- Research article
- Open Access
Comparison of linkage disequilibrium and haplotype diversity on macro- and microchromosomes in chicken
BMC Genetics volume 10, Article number: 86 (2009)
The chicken (Gallus gallus), like most avian species, has a very distinct karyotype consisting of many micro- and a few macrochromosomes. While it is known that recombination frequencies are much higher for micro- as compared to macrochromosomes, there is limited information on differences in linkage disequilibrium (LD) and haplotype diversity between these two classes of chromosomes. In this study, LD and haplotype diversity were systematically characterized in 371 birds from eight chicken populations (commercial lines, fancy breeds, and red jungle fowl) across macro- and microchromosomes. To this end we sampled four regions of ~1 cM each on macrochromosomes (GGA1 and GGA2), and four 1.5 -2 cM regions on microchromosomes (GGA26 and GGA27) at a high density of 1 SNP every 2 kb (total of 889 SNPs).
At a similar physical distance, LD, haplotype homozygosity, haploblock structure, and haplotype sharing were all lower for the micro- as compared to the macrochromosomes. These differences were consistent across populations. Heterozygosity, genetic differentiation, and derived allele frequencies were also higher for the microchromosomes. Differences in LD, haplotype variation, and haplotype sharing between populations were largely in line with known demographic history of the commercial chicken. Despite very low levels of LD, as measured by r2 for most populations, some haploblock structure was observed, particularly in the macrochromosomes, but the haploblock sizes were typically less than 10 kb.
Differences in LD between micro- and macrochromosomes were almost completely explained by differences in recombination rate. Differences in haplotype diversity and haplotype sharing between micro- and macrochromosomes were explained by differences in recombination rate and genotype variation. Haploblock structure was consistent with demography of the chicken populations, and differences in recombination rates between micro- and macrochromosomes. The limited haploblock structure and LD suggests that future whole-genome marker assays will need 100+K SNPs to exploit haplotype information. Interpretation and transferability of genetic parameters will need to take into account the size of chromosomes in chicken, and, since most birds have microchromosomes, in other avian species as well.
Accurate characterization of linkage disequilibrium (LD) within and between populations is important in domestic animal genetic studies because LD underlies all forms of mapping studies  and is an important parameter in the design of marker panels. A recent study on LD in the chicken (Gallus gallus)  showed substantial differences in the extent of LD between populations, which is consistent with pronounced differences observed between breeds of other domesticated species [3, 4]. The extent of LD in broilers estimated by SNP genotyping was shown to be very limited with r2 values below 0.3, on average, for distances between 0.25 and 1 cM . Studies in layers suggest higher extent of LD compared to broilers, although studies on LD in chicken other than broilers are sparse [5, 6]).
Despite low levels of LD in broilers, Andreescu et al  showed that there is significant overlap in the LD of marker pairs between different populations. This result suggests that haplotypes are shared between commercial broiler lines. Estimating haplotype sharing between populations is of interest because it can aid in predicting transferability of genetic parameters, such as genomic estimated breeding values (GEBV)  or QTL, from one population to another. Characterization of LD and haplotype sharing has been achieved in humans by making extensive haplotype maps of different populations [8, 9]. These maps confirmed the organization of haplotypes in so called haploblocks. The elucidation of genome-wide haploblock structure has been beneficial for designing SNP genotyping assays that capture the maximum amount of haplotype diversity in humans . The design of a genome-wide haplotype map for chicken would be equally beneficial for the design of high-density genotyping assays.
In most birds, the genome is organized into a few very large chromosomes and many very small chromosomes[11, 12]. The best studied bird genome to date, that of the chicken, has five so called macrochromosomes (GGA1-5) that range in size from 50 to 200 Mb, five intermediate chromosomes (GGA6-10) ranging from 20 to 40 Mb, and 28 microchromosomes (GGA11-38) that average ~12 Mb. The smallest of the microchromosomes were estimated to be less than 5 Mb.
The microchromosomes have structural differences compared to the larger chromosomes, such as higher GC content. Intergenic distances on the microchromosomes are also lower as well as the average size of the introns, resulting in a much higher gene density compared to the macrochromosomes . Recombination rates in microchromosomes (50 - 100 kb/cM) are much higher compared to the macrochromosomes (~300 kb/cM), possibly resulting from the requirement of at least one chiasma per chromosome per meiosis, and possibly facilitated by a higher density of cohesin binding sites[14, 15].
The higher recombination rate on microchromosomes is expected to reduce LD compared to macrochromosomes. This conclusion was ascertained by Aerts et al for one microchromosome indicating that LD may not be the same throughout the chicken genome. However, a systematic survey of differences in haplotype structure between micro- and macrochromosomes in birds has not been done. Recombination rate was also shown to correlate positively with nucleotide diversity in chicken . The degree of haplotype sharing is expected to decline with increasing recombination rate, which should result in lower transferability of genetic parameters for microchromosomes. Other measures such as haploblock structure, are also expected to be affected by differences in recombination rate. All SNP based studies so far have used SNP densities that were insufficient to ascertain LD and haplotype structure in chicken.
We applied a focused high-density SNP typing strategy, sampling 1 SNP every 2 kb, to quantify the differences in LD, haplotype diversity, haplotype sharing and haploblock structure between the micro- and macrochromosomes in chicken. For a good representation of both types of chromosomes, we sampled four regions of ~1 cM each on macrochromosomes (GGA1 and GGA2), and four 1.5 -2 cM regions on microchromosomes (GGA26 and GGA27). To test the generality of our conclusions, a broad and diverse set of populations were genotyped including all the important commercial types (white egg and brown egg layers, sire and dam broilers), as well as two traditional Dutch fancy breeds and a wild chicken population.
Commercial and non-commercial populations were surveyed to provide a broad view of LD patterns and haplotypic diversity in chicken. Unrelated individuals were taken from each population. For the commercial populations, representatives of all major types were selected: one white egg layer (E2), one brown egg layer (B2), two female or dam broiler lines (one of a closed line - E5, and one of an open line - A3), and one male or sire broiler line (E3). The commercial lines were provided by Hendrix Genetics; E2, E5 and E3 lines were surveyed previously . Furthermore, two traditional Dutch breeds (Owl Bearded, AvDiv_09, and Frisian Fowl, AvDiv_10), and one wild chicken population (Gallus gallus spadiceus, AvDiv_101) were sampled; these populations were previously surveyed in the AvianDiv project . We also included seven Gallus lafayetii (Ceylon Jungle Fowl, SRI) specimens, from Sri Lanka, for outgroup comparison and determining the ancestral state of the SNPs. Additional details on the populations and sampling can be found in Table 1.
SNP selection and typing
SNPs were selected from dbSNP 125 and mapped to Gallus gallus build 2.1. Four chromosomes, two macro- (GGA1 and GGA2) and two microchromosomes (GGA26 and GGA27) were each surveyed at two regions, for a total of eight regions. The size of each region was ~300 kb (~1 cM) for the macrochromosomes, and ~150 kb (~1.5-2 cM) for the microchromosomes [14, 15]. Regions were selected for having sufficient SNP information to allow selection of one SNP per 2 kb and good reliability of the assembly. Furthermore, regions on the same chromosome were chosen to be far apart, to minimize effects of hitchhiking due to linkage. One SNP was selected per 2 kb, on average, with a total of 889 SNPs (Additional File 1).
Genotyping was performed using the GoldenGate/Sentrix Array technology from Illumina , according to manufacturer's protocols. The 889 SNPs were part of a larger 1536-plex assay.
LD was calculated as pairwise r2 and D' values using Haploview 4.0  for each of the populations and for each of the genomic regions. For each population, only markers with a 75% or higher genotyping success, without significant deviation from Hardy-Weinberg (p < 0.001), and with a minor allele frequency (MAF) greater than 5% were included in the analysis.
where LD ijk is the observed LD for marker pair i of population j in region k, d ijk is the distance in bp for marker pair i of population j in genomic region k, β jk is the coefficient that describes the decline of LD with distance for population j in genomic region k and e ijk is a random residual. For each genomic region within population was estimated using the nls function in the R environment http://www.r-project.org/.
Population and genomic region effects on LD extent were tested using LDc ijk , which is the distance corrected and variance stabilized LD for marker pair i in genomic region k and breed j and it was estimated using and obtained with equation 1 :
Differences in LD between genomic regions and populations were analyzed by testing their significance when included as fixed effect in a linear model .
Effective population sizes were estimated by transforming physical distances to genetic distances (300 kb/cM for the macrochromosomes, and 65 kb/cM for the microchromosomes ). Past effective population sizes were calculated using a sliding window and transforming physical distances to genetic distances to estimate the number of past generations as [22–24]
General population statistics
Observed heterozygosity and MAF were calculated with custom Perl scripts. The ancestral state of the SNPs was determined by assessing the state of the SNP in the seven Ceylon Jungle Fowl samples, and the allele present in the same monomorphic state in all seven animals was inferred to be the ancestral state for G. gallus. Tests for the effect of genomic region on heterozygosity and on ancestral frequency were done using a linear model.
Genetic distances were calculated based on allele frequencies using the Gendist program, and a neighbor joining tree was constructed using the Neighbor program from the Phylip package .
Haplotype and haploblock analysis
Haploview was used for inferring haplotypes and haplotype frequencies using the '-blocks' option , using custom Perl scripts to generate block definitions and to collect haplotypes and haplotype frequencies for a sliding window (scripts available upon request from HJM). Haplotype homozygosity (HH) was calculated as the sum of products of haplotype frequencies . Haploblock structure was determined using two haploblock rules, the Gabriel rule  and the 4 gamete rule, as implemented in Haploview 4.0 , based on haplotypes with 5% or higher occurrence.
Haplotype sharing between two populations was calculated as the number of haplotypes that were shared by both populations divided by the average number of haplotypes in these populations. The average number of haplotypes was calculated by taking the sum of the haplotype count of both populations divided by two. Only haplotypes with greater than 5% occurrence in each of the populations were considered.
Of 889 SNPs assayed, 806 were successfully genotyped. A 90% success rate is comparable to previous data sets generated for chicken on the same platform . The numbers of SNPs successfully genotyped per region can be found in Additional File 1. Of these 806 SNPs, 91% had a MAF greater than 5% in at least one population. The wild chicken population (AvDiv101) had the highest percentage of polymorphic markers (70%). The dam broiler line A3 was a close second with slightly less than 70% of SNPs polymorphic. The white egg layer line was clearly the population with the smallest number of polymorphic markers with less than 35% of SNPs informative (Additional File 2).
Heterozygosity was lower (P < 0.005) for the macrochromosomes (0.193) compared to the microchromosomes (0.207) when analyzed across all populations. Within most populations, heterozygosity was found to vary between chromosomal regions; for example, a three-fold difference was observed between GGA26 and GGA27 in the white egg layer population (Table 1).
Of the SNPs genotyped in Ceylon Jungle Fowl, 729 (90%) were successfully genotyped and 669 (92%) were not polymorphic in the sample of seven individuals. The alleles present in Ceylon Jungle Fowl were putatively inferred as the ancestral allele in G. gallus. Taken over all populations, ancestral frequencies were slightly lower in the microchromosomes (0.669 vs. 0.680 for macrochromosomes), but the difference was not significant. Population allele frequencies were skewed towards the ancestral state. The highest ancestral frequency was observed in the wild chicken population, and the lowest in the white egg layer and the Dutch breeds (Table 1, Additional File 3). Taken over all populations, ancestral frequencies were slightly lower in the microchromosomes (0.669, compared to 0.680 for macrochromosomes), but the difference was not significant.
Across all populations, LD for the microchromosomes was significantly lower compared to the macrochromosomes (P < 0.0001), and fitted values for LD were consistently lower for the microchromosomes for all populations. Differences in LD between micro- and macrochromosomes from the global fit to the Sved equation resulted in a 2.8× lower estimated recombination rate for the latter (Table 2). Observed values were almost consistently lower for the microchromosomes, although a few local exceptions were observed (e.g. dam broiler E5 showed somewhat higher LD at the microchromosomes at short distances).
Observed values for r2 were never < 0.1 even at distances of ~1 cM, contrary to the fit to the Sved equation predicted, and the r2 values never reached high values (> 0.9) even at very small marker distances. This resulted in a much flatter observed versus predicted LD curves. Observed values for D' showed similar trends as r2, but were always much higher. Average D' never was < 0.5 at the macrochromosomes for any population even at 250 kb (Figure 1).
The white egg layer showed the highest LD of all populations, followed by the brown layer and the dam broiler line E5. The wild chicken had, in general, the lowest extent of LD, closely followed by dam broiler line A3 and the sire broiler line. The differences between the breeds were expected to depend on their effective population sizes, which were estimated to be between 40 and 1200 for the domesticated chicken populations, and > 2000 for the wild chicken population, based on all the macrochromosomes and estimated across all marker distances. Based on microchromosomes, Ne was estimated systematically lower for all populations, on average 1.6 times lower (Table 2). When marker distance was taken into account to allow estimation of Ne for a given point in the population history, a continuous reduction in Ne was observed for all populations (Additional File 4).
The proportion of the regions captured by haploblocks was consistently lower for the microchromosomes compared to the macrochromosomes (P < 0.001, sign test). Haploblock sizes also tended to be smaller for the microchromosomes (Figure 2). Only in layers were more than 10% of the SNPs on the microchromosomes captured in blocks > 40 kb. Congruence in haploblock boundaries - both for micro- and macrochromosomes - between populations was very limited (Additional File 5).
Haploblock structure varied widely between populations (Figure 2). The white egg layer had large (> 100 kb) blocks covering 71% of the regions. In the brown egg layer and dam broiler line E5, block structure was still considerable with well over 40% of regions in blocks larger than 10 kb. There were pronounced differences between the results from different methods of block inference. In general, the Gabriel method inferred far fewer blocks that tended to be larger compared to the 4 Gamete Rule (Additional File 6). Nevertheless, the overall observations on differences between populations, and between micro- and macrochromosomes were consistent between the two methods.
Haplotype homozygosity and evidence for selective sweep
Haplotype homozygosity (HH) measured over a sliding window with bin sizes of 10 SNPs (~20 kb) ranged from 0.11 to 1 (Figure 3). Over all populations, the microchromosomes showed consistently lower average HH, with the exception of dam broiler line E5. For the macrochromosomes a relatively small number (1-7) of haplotypes accounted for the vast majority of haplotype diversity (> 90%) in most domesticated populations while in the microchromosomes a larger number of haplotypes tended to explain a smaller part of the variation (Additional File 7). The white egg layer displayed extended regions with only a single haplotype, while HH was between 0.1 and 0.2 in Red Jungle Fowl for large parts of the regions covered in this study. All other populations showed a wide range of HH between and within regions.
Haplotype sharing on the microchromosomes was substantially lower (P < 0,0001, sign test) compared to the macrochromosomes (Figure 4C, Table 3), between 25 and 50% of sharing in macrochromosomes for haplotypes of the same physical size. Genetic distances between populations were also consistently larger (Figure 4B, P < 0.0001, sign test) based on genotypes derived from microchromosomal SNPs, but the difference was less pronounced compared to haplotype sharing.
Haplotype sharing on the macrochromosomes was highest between the two dam broiler lines (Table 3). It was also high between the sire broiler line and the dam broiler lines, and between the two Dutch traditional breeds. The two Dutch traditional breeds exhibited haplotype sharing with the white egg layer. Patterns of sharing were consistent with genetic distances based on genotype data (Figure 4A).
Our aim was to study differences in LD and haplotype variation between micro- and macrochromosomes, using a very high marker density across populations reflecting commercial diversity as well as fancy breeds and wild chicken. Chromosomal regions were chosen to represent the maximum range in size of macro- (GGA1 and GGA2) and microchromosomes (GGA26 and GGA27); the two microchromosomes were among the smallest well-assembled chromosomes available within the current genome build . Selection of SNPs was based solely on position (with the requirement of having 1 SNP every 2 kb), and thus systematic bias due to SNP selection was unlikely. Populations were chosen to reflect variation in the degree of polymorphism, and hence expected LD, to the widest possible extent, with the white egg layer at the lower end and Red Jungle Fowl at the upper end . Extent of LD in chicken has been studied before [2, 5, 6] but these studies were limited in numbers of markers, marker density, population sampling or sampling across chromosomes to accurately and comprehensively asses LD to the same degree as the present study.
Based on higher recombination rates in microchromosomes compared to macrochromosome differences in LD and haplotype variation were expected, but measures of these differences have not been previously reported. We found LD, HH, haploblock structure, and haplotype sharing all consistently lower for microchromosomes compared to macrochromosomes when measured using physical distance. A direct effect of recombination on these measures comes from changing the relationship from physical distance to genetic distance. From the fit of LD to the Sved equation , and assuming that Ne is the same for all chromosomes, the recombination rate was estimated to be on average 2.8 times higher at the microchromosomes (Table 2). This difference is less than the expected 4.5 times higher recombination rate for microchromosomes compared to macrochromosomes . The recombination frequency for the microchromosomes based on LD, therefore, appears to be systematically underestimated for all populations. While regional differences in recombination frequency are expected, currently no recombination map is available providing information at the scale of the present study (< 1 cM scale), not even for the macrochromosomes. For the smallest microchromosomes current recombination maps are even less detailed.
The inferred rate of 2.8× smaller recombination rate for macrochromosomes compared to microchromosomes, which is inconsistent with previous estimates (~4.5×, ), is due to a bias in the analysis from fitting the Sved equation across the same physical distance in micro- and macro chromosomes. LD at different distances has been shown to relate to effective population sizes at different numbers of generations by 1/(2 c), where c is the median distance between markers in Morgan . By performing local fits to the data, using SNP distance bins that are similar in genetic rather than physical distances, the systematic difference in Ne between micro- and macrochromosomes disappears. For most populations, past population sizes derived from both classes of chromosomes become quite similar when measured against genetic distance (Additional File 4).
Since Ne does not seem to deviate systematically once distances are properly corrected for differences in recombination rate, the main explanation for observed differences in heterozygosity, genotype differentiation, and derived allele frequencies in the microchromosomes is higher mutation rate. Higher heterozygosity is known to be positively correlated to recombination rate [16, 28], although the mechanism is not fully understood. We found derived allele frequency to be slightly higher on the microchromosomes, which suggests a higher evolutionary rate. A higher evolutionary rate for microchromosomes has been found before in a comparison between chicken and turkey macro- and microchromosomes . Higher levels of differentiation could result from increased background directional selection for higher GC content in the microchromosomes. The effect of directional selection would have the same effect as a smaller effective population size. Since there is no evidence for differences in Ne a higher mutation rate seems to be the best explanation for higher genetic differentiation on the microchromosomes.
The Sved equation assumes a static population size . However, the fact that the observed values or r2 (Figure 1) show more of a flat line compared to the expected values of r2 based on the fit to the global Sved equation is an indication of declining population size[22, 24]. Fits based on local inter-marker distance-bins reveal declining effective population sizes as shown in Additional File 4. Differences in LD and derived effective population sizes are largely consistent with known population histories, with white egg layers known to be more inbred than other breeds, while most of the commercial broiler lines are considered outbred [18, 30, 17]. Nevertheless, the dam broiler E5 has been a closed line for many generations (AV, unpublished results), which explains higher LD and HH in this population. The decline of effective sizes for the eight populations is consistent with earlier findings of substantial loss of allelic variation in domesticated chicken , reemphasizing the concern to maintain genetic diversity in this species.
In humans, markers diagnostic for haplotypes, so called tag SNPs, are often transferable between populations because of haplotype sharing and populations having common haploblock boundaries . In the chicken, haploblock boundaries show little overlap between populations, and haplotype sharing between populations is low. This difference between the two species could be the result of differences in demography, with the block-like structure of haplotype variation in humans being the result of population expansion in the past 10+ thousand of years originating from a population with an effective size of thousands to tens of thousands at most . Conversely, the present study finds evidence for population contraction in chicken, which is consistent with the relatively small number of long haplotypes and levels of haplotype homozygosity. These long current haplotypes are expected to be a mosaic of a much higher diversity of small past haplotypes, similar to what is observed in dogs . The ancient small haplotypes that make up today' s longer haplotypes, therefore, do not result in a very high r2 (unless inbreeding becomes very high and only a very small number of haplotypes remain). They do, however, result in high D' as the inbreeding erodes away many of the possible - and previously existing - haplotypes in a population. D' will more often result in high LD values when only part of all possible haplotypes are present compared to r2[32, 33]. Since the block construction methods applied here were based on D' it was not surprising to find considerable block structure - albeit often found in small blocks.
The block structure in the genomes of layers may be exploited to make genome-wide marker assays with 10,000 to 20,000 well chosen tag SNPs that would cover around 70% of the genome, supplemented by a similar number of SNPs to survey the remaining ~30%. Since block structure is mostly at the scale of < 10 kb for the more outbred broiler populations, and LD (measured as r2) near 0.2 at a similar scale, the number of informative SNPs would need to be > 100,000. However, as tag SNPs are probably not highly transferable between commercial populations, a general purpose assay might even need many more markers than 100 K.
Understanding sharing of haplotypes between populations is of further importance as it determines the success of transferring genetic parameters from one population to another . The present study confirms the findings of Andreescu et al  in that high overlap in haplotypes between broilers exists. However, it appears to only exist between closely related populations. Transferability of marker information between more distantly related populations may be problematic. For the microchromosomes, haplotype sharing is very small even among the broilers, showing that population-to-population transferability of marker information should be treated differently for micro- and macrochromosomes at the same physical scale. Since genotype differentiation is also systematically higher for the microchromosomes, differences in haplotype sharing are likely the result of both increased mutation rate and recombination frequency for the microchromosomes.
Patterns of LD, haplotype variation, and haplotype sharing, as well as genotype variation and genotype differentiation, are all different in the microchromosomes compared to macrochromosomes in chicken. While differences in LD are congruent with differences in recombination rate, differences in haplotype differentiation may be partly explained by an increased genotype differentiation. Differences in genotype differentiation seem best explained by a higher mutation rate for the microchromosomes. It is vital that whole-genome studies in chicken take these differences into account, both in the genotype assay design phase, as well as in interpretation and application of results. Because most birds have microchromosomes it is likely that the findings presented in this study are relevant to a wider group of avian species.
Andersson L, Georges M: Domestic-animal genomics: deciphering the genetics of complex traits. Nat Rev Genet. 2004, 5: 202-212. 10.1038/nrg1294.
Aerts J, Megens HJ, Veenendaal T, Ovcharenko I, Crooijmans R, Gordon L, Stubbs L, Groenen M: Extent of linkage disequilibrium in chicken. Cytogenet Genome Res. 2007, 117: 338-45. 10.1159/000103196.
Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, Clamp M, Chang JL, Kulbokas EJ, Zody MC, Mauceli E, Xie X, Breen M, Wayne RK, et al: Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005, 438: 803-19. 10.1038/nature04338.
Amaral AJ, Megens H, Crooijmans RPMA, Heuven HCM, Groenen MAM: Linkage Disequilibrium Decay and Haplotype Block Structure in the Pig. Genetics. 2008, 179: 569-579. 10.1534/genetics.107.084277.
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. 10.1534/genetics.107.082206.
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-1181. 10.1534/genetics.105.040782.
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-1512. 10.1534/genetics.107.084301.
International HapMap Consortium: A haplotype map of the human genome. Nature. 2005, 437: 1299-320. 10.1038/nature04226.
International HapMap Consortium: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449: 851-61. 10.1038/nature06258.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296: 2225-9. 10.1126/science.1069424.
Griffin DK, Robertson LBW, Tempest HG, Skinner BM: The evolution of the avian genome as revealed by comparative molecular cytogenetics. Cytogenet Genome Res. 2007, 117: 64-77. 10.1159/000103166.
Nishida C, Ishijima J, Kosaka A, Tanabe H, Habermann FA, Griffin DK, Matsuda Y: Characterization of chromosome structures of Falconinae (Falconidae, Falconiformes, Aves) by chromosome painting and delineation of chromosome rearrangements during their differentiation. Chromosome Res. 2008, 16: 171-81. 10.1007/s10577-007-1210-6.
Chicken genome consortium: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.
Groenen MA, Cheng HH, Bumstead N, Benkel BF, Briles WE, Burke T, Burt DW, Crittenden LB, Dodgson J, Hillel J, Lamont S, de Leon AP, Soller M, Takahashi H, Vignal A: A consensus linkage map of the chicken genome. Genome Res. 2000, 10: 137-47.
Groenen MAM, Wahlberg P, Foglio M, Cheng HH, Megens H, Crooijmans RPMA, Besnier F, Lathrop M, Muir WM, Wong GK, Gut I, Andersson L: A high-density SNP-based linkage map of the chicken genome reveals sequence features correlated with recombination rate. Genome Res. 2009, 19: 510-9. 10.1101/gr.086538.108.
Webster MT, Axelsson E, Ellegren H: Strong Regional Biases in Nucleotide Substitution in the Chicken Genome. Mol Biol Evol. 2006, 23: 1203-1216. 10.1093/molbev/msk008.
Muir WM, Wong GK, Zhang Y, Wang J, Groenen MAM, Crooijmans RPMA, Megens H, Zhang H, Okimoto R, Vereijken A, Jungerius A, Albers GAA, Lawley CT, Delany ME, MacEachern S, Cheng HH: Genome-wide assessment of worldwide chicken SNP genetic diversity indicates significant absence of rare alleles in commercial breeds. Proc Natl Acad Sci USA. 2008, 105: 17312-7. 10.1073/pnas.0806569105.
Hillel J, Groenen MA, Tixier-Boichard M, Korol AB, David L, Kirzhner VM, Burke T, Barre-Dirie A, Crooijmans RP, Elo K, Feldman MW, Freidlin PJ, Mäki-Tanila A, Oortwijn M, Thomson P, Vignal A, Wimmers K, Weigend S: Biodiversity of 52 chicken populations assessed by microsatellite typing of DNA pools. Genet Sel Evol. 2003, 35: 533-557. 10.1186/1297-9686-35-6-533.
Fan JB, Oliphant A, Shen R, Kermani BG, Garcia F, Gunderson KL, Hansen M, Steemers F, Butler SL, Deloukas P, Galver L, Hunt S, McBride C, Bibikova M, Rubano T, Chen J, Wickham E, Doucet D, Chang W, Campbell D, Zhang B, Kruglyak S, Bentley D, Haas J, Rigault P, Zhou L, Stuelpnagel J, Chee MS: Highly parallel SNP genotyping. Cold Spring Harb Symp Quant Biol. 2003, 6869-78.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-5. 10.1093/bioinformatics/bth457.
Sved JA: Linkage disequilibrium and homozygosity of chromosome segments in finite populations. Theor Popul Biol. 1971, 2: 125-41. 10.1016/0040-5809(71)90011-6.
Hayes BJ, Visscher PM, McPartlan HC, Goddard ME: Novel multilocus measure of linkage disequilibrium to estimate past effective population size. Genome Res. 2003, 13: 635-43. 10.1101/gr.387103.
Tenesa A, Navarro P, Hayes BJ, Duffy DL, Clarke GM, Goddard ME, Visscher PM: Recent human effective population size estimated from linkage disequilibrium. Genome Res. 2007, 17: 520-6. 10.1101/gr.6023607.
Villa-Angulo R, Matukumalli L, Gill C, Choi J, Van Tassell C, Grefenstette J: High-resolution haplotype block structure in the cattle genome. BMC Genetics. 2000, 10: 19-10.1186/1471-2156-10-19.
Felsenstein : PHYLIP (Phylogeny Inference Package), version 3.66. 2008, Distributed by the author. Department of Genome Sciences, University of Washington, Seattle
Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, Ramachandran S, Cann HM, Barsh GS, Feldman M, Cavalli-Sforza LL, Myers RM: Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008, 319: 1100-4. 10.1126/science.1153717.
Muir W, Wong G, Zhang Y, Wang J, Groenen M, Crooijmans R, Megens H, Zhang H, McKay J, McLeod S, Okimoto R, Fulton J, Settar P, O'Sullivan N, Vereijken A, Jungerius A, Albers G, Taylor Lawley C, Delany M, Cheng H: Review of the Initial Validation and Characterization of a 3 K Chicken SNP Array. World's Poultry Science Journal. 2008, 64: 219-226. 10.1017/S0043933908000019.
Fang L, Ye J, Li N, Zhang Y, Li S, Wong G, Wang J: Positive correlation between recombination rate and nucleotide diversity is shown under domestication selection in the chicken genome. Chinese Science Bulletin. 2008, 53: 746-750. 10.1007/s11434-008-0159-y.
Axelsson E, Webster MT, Smith NG, Burt DW, Ellegren H: Comparison of the chicken and turkey genomes reveals a higher rate of nucleotide divergence on microchromosomes than macrochromosomes. Genome Res. 2005, 15: 120-125. 10.1101/gr.3021305.
Twito T, Weigend S, Blum S, Granevitze Z, Feldman MW, Perl-Treves R, Lavi U, Hillel J: Biodiversity of 20 chicken breeds assessed by SNPs located in gene regions. Cytogenet Genome Res. 2007, 117: 319-26. 10.1159/000103194.
de Bakker PIW, Burtt NP, Graham RR, Guiducci C, Yelensky R, Drake JA, Bersaglieri T, Penney KL, Butler J, Young S, Onofrio RC, Lyon HN, Stram DO, Haiman CA, Freedman ML, Zhu X, Cooper R, Groop L, Kolonel LN, Henderson BE, Daly MJ, Hirschhorn JN, Altshuler D: Transferability of tag SNPs in genetic association studies in multiple populations. Nat Genet. 2006, 38: 1298-303. 10.1038/ng1899.
Devlin B, Risch N: A Comparison of Linkage Disequilibrium Measures for Fine-Scale Mapping. Genomics. 1995, 29: 311-322. 10.1006/geno.1995.9003.
VanLiere JM, Rosenberg NA: Mathematical properties of the r2 measure of linkage disequilibrium. Theor Popul Biol. 2008, 74: 130-137. 10.1016/j.tpb.2008.05.006.
We would like to thank Andreia Amaral for aid in statistical analysis and Coline van Moorsel for comments on the manuscript. This work was supported in part by United States Department of Agriculture (USDA) National Research Initiative Competitive Grants Program (NRICGP) grant no. 2004-05434.
The authors declare that they have no competing interests.
HJM carried out genotyping and statistical analyses and drafted the manuscript. JB helped in interpreting results and drafting the manuscript. MG and RC coordinated the research. RC prepared DNA samples. RJ designed the marker assay. WM and HC helped in drafting and revising the manuscript. HK and AC aided in statistical analyses. RAV, OH, and PS provided DNA samples. All authors read, edited and approved the final manuscript.
Electronic supplementary material
Additional file 1: An overview of the SNP selection. Two regions were sampled per chromosome, each region was ~1 cM in size. SNPs were selected with a spacing of ~2 kb. Coordinates of the regions are based on chicken genome build WASHU2. (PDF 25 KB)
Additional file 4: Estimates of past effective popution size. Estimates from microchromosomes were comparable to estimates from macrochromosomes. (PDF 417 KB)
Additional file 5: Haploblock structure in the eight selected regions for eight populations derived from the Gabriel method . Top and middle sections show two regions for GGA 1 and GGA 2 (macrochromosomes). Bottom section shows the four regions derived from the two microchromosomes (GGA26 and GGA27). 1 = white egg layer; 2 = brown egg layer; 3 = dam broiler E5 (closed line); 4 = sire broiler line; 5 = dam broiler A3 (open line); 6 = Frisian Fowl; 7 = Owl Bearded; 8 = Red Jungle Fowl. (PDF 198 KB)
Additional file 6: Proportion of macro- and microchromosomes capture in haploblocks of different size. Block definitions were according to the 4 Gamete Rule . (PDF 192 KB)
Additional file 7: estimated number of major haplotypes (haplotypes occurring with a frequency > 5% in the populations studied). The second column for each of the populations shows the percentage of haplotype diversity explained by the major haplotypes. (PDF 36 KB)
About this article
Cite this article
Megens, HJ., Crooijmans, R.P., Bastiaansen, J.W. et al. Comparison of linkage disequilibrium and haplotype diversity on macro- and microchromosomes in chicken. BMC Genet 10, 86 (2009). https://doi.org/10.1186/1471-2156-10-86
- Linkage Disequilibrium
- Recombination Rate
- Effective Population Size
- High Linkage Disequilibrium
- Haplotype Sharing