A genetic variation map for chicken with 2.8 million single nucleotide polymorphisms



We describe a genetic variation map for the chicken genome containing 2.8 million single nucleotide polymorphisms (SNPs), based on a comparison of the sequences of 3 domestic chickens (broiler, layer, Silkie) to their wild ancestor Red Jungle Fowl (RJF). Subsequent experiments indicate that at least 90% are true SNPs, and at least 70% are common SNPs that segregate in many domestic breeds. Mean nucleotide diversity is about 5 SNP/kb for almost every possible comparison between RJF and domestic lines, between two different domestic lines, and within domestic lines – contrary to the idea that domestic animals are highly inbred relative to their wild ancestors. In fact, most of the SNPs originated prior to domestication, and there is little to no evidence of selective sweeps for adaptive alleles on length scales of greater than 100 kb.

Keywords: chicken, polymorphism, domestication, selection


The generation of a high quality draft sequence for the genome of chicken (Gallus gallus) is an important advance. Chickens are good models for studying the genetic basis of phenotypic traits, because of the extensive diversity among domestic chickens selected for different purposes. Monogenic traits are well-studied, but many interesting traits are complex and determined by an unknown number of genes. Quantitative trait loci (QTLs) have been mapped for a range of traits, including ones for growth, body composition, egg production, antibody response, disease resistance, and behaviour. Determining causative genes is difficult, since each locus controls only a fraction of the phenotypic variance. We will describe a survey of the genetic variation between 3 domestic chickens and their wild ancestor. The 2.8 million single nucleotide polymorphisms (SNPs) that we identified will facilitate mapping of complex traits in many ways. First, improved marker density allows researchers to take advantage of the higher recombination rates in chicken, which are 2.5 to 21 cM/Mb depending on the chromosome, as compared to 1 cM/Mb for human and 0.5 cM/Mb for mouse. The previous linkage map used 2000 markers,, but only 800 of these were microsatellites or SNPs, which are the most useful. More importantly, our new data allow researchers to construct detailed haplotypes that segregate in different QTL crosses. Because any mutation underlying a QTL must once have originated from a single founder animal, haplotype comparisons will facilitate the fine mapping of QTLs. To this end, we conduct a genomewide search for evidence of selection due to domestication, and provide an initial characterization of the expected magnitude of these effects.

Genetic variation and utility

Our experiment is outlined in Figure 1. SNPs are generated by partial sequencing at ¼ coverage for each of 3 domestic breeds (a male broiler, a female layer, and a female Silkie), and comparison of the resultant reads to the 6.6x genome for the wild ancestor of domestic chickens, Red Jungle Fowl (RJF). We expect marked heterozygosity within the 3 domestic lines, but not within RJF because the sequenced bird for the genome project is from a highly inbred line that is essentially homozygous.

An external file that holds a picture, illustration, etc.
Object name is nihms-704-f0001.jpg

SNP discovery experiment. We sampled 3 domestic chickens at 1/4 coverage each and compared the resultant sequence to the 6.6x draft genome of Red Jungle Fowl (RJF). Chicken photographs shown here are provided by Bill Payne (RJF), Paul Hocking (broiler), Leif Andersson (layer), and Ning Yang (Silkie).

Comparing the sequence reads for broiler, layer, and Silkie to the genome of RJF, we identified nearly a million SNPs in each instance, at mean rates of about 5 SNP/kb, as shown in Table 1. Notice that all of the “SNP rates” quoted in this paper are computed as nucleotide diversities (π), and given in units of π×103. After correcting for SNPs detected in more than one line, there are 2,833,578 variant sites, or one potential marker every 374 bp along the 1.06 Gb genome. To assess the reliability of these data, we resequenced 295 SNPs in the same bird in which it was detected (Table S1). As many as 94% of the SNPs were confirmed. However, confirmation rates are sensitive to the functional context (e.g., coding versus non-coding) and SNPs in rare categories are less likely to be confirmed. In fact, only 83% of the non-synonymous SNPs were confirmed. Small indels of a few base pairs in length (mean of 2.3 and median of 1) are detected at rates that are well correlated with the corresponding SNP rates, but smaller by about a factor of 10.

Table 1

Frequency of SNPs in different comparisons of RJF and the 3 domestic chicken lines. In addition, we show comparisons involving 3.8-Mb of finished BAC sequence from another line of the layer (White Leghorn) breed. SNP rates are an estimate of nucleotide diversity (π), as embodied by the effective length, which considers how much of the data is of sufficiently good quality to actually detect SNPs and the probability that overlapping reads might be derived from homologous chromosomes

# of SNPs L(effective) SNP/kb
Wild versus domestic
RJF-Broiler 1,041,948 197,431,517 5.28
RJF-Layer 889,377 170,586,544 5.21
RJF-Silkie 1,217,817 217,841,171 5.59
Between domestic lines
Broiler-Layer 194,605 37,506,800 5.19
Broiler-Silkie 257,849 47,554,311 5.42
Layer-Silkie 246,954 42,682,304 5.79
Within domestic lines
Broiler-Broiler 59,227 13,835,075 4.28
Layer-Layer 40,412 10,863,595 3.72
Silkie-Silkie 83,630 15,253,383 5.48
Compare to layer BACs
RJF-to-BAC 20,925 3,809,567 5.49
BAC-Broiler 4,404 847,456 5.20
BAC-Layer 3,904 740,392 5.27
BAC-Silkie 5,089 925,738 5.50

Chicken autosomes are sorted by size into 5 large macrochromosomes (GGA1-5), 5 intermediate chromosomes (GGA6-10), and 28 microchromosomes (GGA11-38). SNP and indel rates are independent of chromosome size, as shown in Figure 2. GGA16 is the sole exception, because it contains the highly variable MHC. This result is surprising, as recombination rates on microchromosomes are much higher than on macrochromosomes and studies in other organisms exhibit a positive correlation between recombination rates and polymorphism rates. We expect that higher gene densities on microchromosomes likely counteract the effect of higher recombination rates.

An external file that holds a picture, illustration, etc.
Object name is nihms-704-f0002.jpg

SNP and indel rates versus chromosome number. We excluded all sequences with “random” chromosome positions. Because of the assembly problems on W, it is not shown. The rates are computed as an average of all 3 domestic lines.

SNP rates between and within chicken lines can be determined from the overlaps between reads. Table 1 demonstrates that almost every pairwise combination gives a SNP rate of just over 5 SNP/kb, except for broiler-broiler and layer-layer, which show about 4 SNP/kb, as expected since the sequenced broiler and layer are from closed breeding lines. To ensure that there are no confounding factors from the single read nature of our data, or the complexities of the overlap analysis, we used comparisons to 3.8 Mb of finished BAC sequence of a different White Leghorn from the same breed but not the same line as the layer sequenced herein. 15 chromosomes were sampled, and the results confirm our rates of 5 SNP/kb. In another study of 15 kb of introns in 25 birds from 10 divergent breeds of domestic chickens, an autosomal rate of 6.5 SNP/kb was reported.

To quantify SNP and indel rate variation versus functional context, we considered three gene sets representing 3868 confirmed mRNA transcripts, 995 chicken orthologs of human disease genes, and 17,709 Ensembl annotations from the RJF analysis. Complete details for all 3 lines are tabulated in the supplements (Table S2). An excerpt for broiler is shown in Table 2. Within genes defined by mRNA transcripts, the SNP rates are 3.5, 2.1, 5.7, and 3.4 SNP/kb in 5′-UTR, coding exon, intron, and 3′-UTR regions respectively. In coding regions, indel rates are 43 times smaller than SNP rates. Ka/Ks is 0.098, similar to what is typically seen in vertebrate comparisons. We also studied “conserved non-coding regions” from the RJF analysis. SNP rates are similar to those of coding exons, but indel rates are intermediate to those of coding exons and UTRs, which supports the notion that these regions are functional, but may not encode proteins.

Table 2

RJF-Broiler polymorphisms

Frequency of sequence polymorphisms between RJF and broiler, decomposed by functional context based on three non-redundant gene sets of 3868 confirmed mRNA transcripts, 995 chicken orthologs of known human disease genes, and 17,709 Ensembl annotations. Human-chicken motifs are conserved sequences that exhibit no evidence of being genic in origin. Gene regions are subdivided into 5′-UTR, coding exon, intron, and 3′-UTR. Ka and Ks indicate non-synonymous and synonymous rates

SNP/kb Indel/kb # of SNP # of Indel
Confirmed mRNA transcripts
 5′-UTR 3.45 0.46 203 27
 coding region 2.11 0.05 1,772 41
  non-synonymous (Ka) 0.73
  synonymous (Ks) 7.44
 introns 5.70 0.52 86,586 7,915
 3′-UTR 3.40 0.42 1,946 243
Human disease genes
 coding region 2.74 0.04 1,005 15
  non-synonymous (Ka) 1.10
  synonymous (Ks) 9.40
 introns 5.36 0.49 27,768 2,553
Ensembl (final version 040427)
 5′-UTR 4.22 0.37 616 54
 coding region 2.71 0.06 12,229 276
  non-synonymous (Ka) 1.17
  synonymous (Ks) 8.28
 introns 5.64 0.52 367,361 33,869
 3′-UTR 3.92 0.43 2,130 236
Human-chicken motifs 2.41 0.25 3,636 379
Genomewide average 5.28 0.48 1,041,948 94,578

Utility of these SNPs depends on their frequency of occurrence in commonly used chicken populations. Hence, we typed 125 SNPs (including coding and non-coding SNPs, randomly distributed across the chicken genome) in 10 unrelated individuals from each of 9 divergent lines representing an assortment of European breeds. This collection includes commercial broiler and layer breeds, standardized breeds selected for their morphological traits, and an unselected breed from Iceland (Table S3). Both alleles segregated in 73% of 1113 successful marker-line combinations (out of 1125 possible combinations). Averaged minor allele frequency is 27%, but it decreases to 20% if marker-line combinations where one of the two alleles is fixed are included. This indicates that a majority of the SNPs are common variants that predate the divergence of modern breeds. Only 12% of the markers had a minor allele frequency of less than 10% in the 90 animals tested.

We demonstrate by example how these data can be used to target specific genome regions. Details of our experiments are in Supplement E (Examples). First, we consider a body weight related QTL on GGA4 that was previously mapped to a 150 cM interval,. After a year of effort, where every known microsatellite (>50) was tested, 26 informative markers were developed. Further progress would have required the laborious sequencing of multiple chickens to find additional polymorphisms in this target region. With the SNP map, we selected 47 random broiler-layer SNPs, and ABI SNPlex assays were developed to type an experimental F2 cross (n = 466). 28 (60%) of these SNPs were informative, but none had breed specific alleles, confirming that most variations predate domestication. In just one month, we doubled the number of markers, and resolved the initial QTL into two QTLs that affect body weight at 3 and 9 weeks of age.

In addition to providing markers for fine mapping, these SNPs are a rich source of candidate polymorphisms for the causative differences underlying important traits. As an example, candidate genes for disease resistance often include TGF-β,, cytokines, and the major histocompatibility complex (MHC). We thus identified 40 SNPs from the SNP map in the coding or promoter regions of 12 cytokine genes. When typed against 8 inbred layer lines, 32 of these SNPs were informative. Cytokine genes on GGA13, including IL4 and IL13, two genes that are expressed in T helper-2 (Th2) cells, drive antibody response. Four of the six SNPs that were polymorphic among lines were in IL4 and IL13, and these SNPs were fixed for different alleles in lines N and 15I, which show differential antibody response to vaccination. These SNPs therefore allow us to test whether the IL4 and IL13 loci directly determine the observed differential antibody response.

Domestication and selection

Domestic animals are useful models of phenotypic evolution under selection. The challenge is to find not only those loci that determine phenotypic differences, but also the causative alleles. We adopt two approaches, searching for evidence of selective sweeps, and for non-synonymous amino acid substitutions at highly conserved sites. One example of a selective sweep is the IGF2 locus in pigs. Given the available data, determining the exact haplotype structure is difficult, because blocks of shared alleles can be erroneously disrupted by heterozygosity of the domestic lines and by sequencing errors. However, we can still search for the local reductions in heterozygosity that accompany selective sweep, as long as we are mindful of the sequencing error rate.

We did 3-way comparisons of RJF and all possible combinations of two domestic lines. Given the limited coverage of the latter, we only examined 100 kb segments with at least 10 SNP sites where each qualifying site must have read coverage from every line. In practice, these segments contained an average of 25 to 28 SNPs. Then, we computed how often 80% or more of the SNP sites are identical in the two domestic lines but different in RJF. In Table S4, we show that 0.4 to 1.5% of the segments qualified. However, when we searched for shared alleles between RJF and one of the domestic lines, 1.2 to 2.6% of the segments qualified. Heterozygosity of the domestic lines is more of a confounding factor in searching for blocks of shared alleles between two domestic lines, versus between RJF and one domestic line. This could explain the difference, but if so, then heterozygosity of the domestic lines is the dominant factor in determining such blocks of shared alleles, not selective sweeps. Hence, selective sweeps that occurred before the divergence of modern domestic breeds must have left behind footprints that are much smaller than 100 kb. This would however be entirely consistent with the historically large effective population size of domestic chickens, and the reported high recombination rates.

For a glimpse of the true haplotype patterns, one can compare the aforementioned 3.8 Mb of finished BAC sequence, from the second layer line (L2), to the genome of RJF. These results are overlaid alongside the primary SNP data set in Figure 3. Short RJF-type fragments can be seen in all 4 lines. Shared domestic-type fragments can also be seen, but at sizes of 5 to 15 kb. This is consistent with our inability to detect footprints of selective sweeps at length scales of 100 kb and suggests that a better choice is 10 kb. However, our data are insufficient for such a genomewide analysis.

An external file that holds a picture, illustration, etc.
Object name is nihms-704-f0003.jpg

Detailed haplotype patterns in 3 regions, each covered by 2 overlapping BACs from the second layer line (L2). The primary SNP data are labeled B (broiler), L1 (layer), and S (Silkie). All comparisons are to RJF, and we show only those sites where a SNP is identified in at least one of the 4 lines. Hence, the horizontal scale is linear in the number of SNP sites, but non-linear for size. BLUE colors indicate where a particular line agrees with RJF, while RED colors indicate where it does not. Overlapping BACs on GGA1 and GGA7, but not GGA14, are clearly from different haplotypes.

It has been hypothesized that loss-of-function (LOF) mutations have accumulated in domestic animals, as the result of relaxed purifying selection and selection for adaptive benefits. An example of the latter is the deletion in the myostatin gene in cattle selected for muscularity. Such deletions are rare, and so we looked for non-synonymous SNPs at highly conserved sites using SIFT. Every substitution is thus classified as being likely to affect function (intolerant) or not (tolerant). For genes defined by mRNA transcripts, 26% of testable SNPs are intolerant, although only 11% are intolerant if we restrict this to high confidence assessments (Table S5). Usually, it is the domestic allele that is intolerant, but we would emphasize that intolerant SNPs are rare, and only 59% were confirmed by PCR resequencing. Given that the domestic allele is represented by a single read, as opposed to 6.6 for the wild allele, much of this effect is likely due to sequencing errors. However, we noticed the same effect in 424 non-synonymous SNPs that we identified from an analysis of 330,000 ESTs, where every allele was seen in two or more ESTs. We conclude that the LOF hypothesis remains intriguing, but any effect is likely to be small.

Some of the experimentally confirmed SIFT intolerant SNPs could be functionally important. We show one example in Figure 4, from the ornithine transcarbamylase (OTC) gene. It substitutes glycine in RJF to arginine in layer and broiler. This SNP is identical to the G188R substitution associated with hyperammonaemia in humans. Resequencing of additional domestic birds revealed a high frequency for the intolerant allele in both White Leghorns (p=0.65, n=20) and in broilers (p=0.75, n=6). In mammals, OTC is expressed in the liver and catalyzes the second step of the urea cycle. Chicken OTC is expressed in the kidney and exhibits a low enzymatic activity, with substantial variability among breeds. Preservation and sequence conservation of OTC, along with all other enzymes in the urea cycle, was unexpected because avian species excrete uric acid (not urea) as their primary component of nitrogenous waste, and were believed to be lacking a functional urea cycle. The deleterious nature of human G188R makes this an attractive candidate for phenotypic studies of avian-specific adaptations in the urea cycle.

An external file that holds a picture, illustration, etc.
Object name is nihms-704-f0004.jpg

Multi-species alignments for ornithine transcarbamylase (OTC), indicating non-synonymous substitutions relative to human protein. SIFT intolerant position is indicated by site number and bold-faced lettering. WT=wild type. MUT=mutant.


This analysis has provided the first global assessment of nucleotide diversity for a domestic animal in comparison to a representative of its wild ancestor. The small number of birds sequenced is compensated for by the vast number of sites examined. We detected surprisingly little difference in diversity in comparisons between RJF and domestic lines, between different domestic lines, and within domestic lines. The total rates are typically 5 SNP/kb, with the only exception being a slight reduction to 4 SNP/kb in broiler and layer lines that are maintained as closed breeding populations. In comparison, 5 SNP/kb is 6 to 7-fold larger than humans and domestic dogs, 3-fold larger than gorillas, but similar to the diversity between different mouse subspecies.

Most of the nucleotide diversity observed between and within domestic lines must have originated prior to the domestication of chickens 5,000 to 10,000 years ago. Given a neutral substitution rate of 1.8×10-9 sites per year for galliform birds, we estimate that a coalescence time of 1.4 million years would be required to account for the observed rates of 5 SNP/kb. Considering that the rates observed between RJF and domestic lines are not much higher than those between domestic lines, it would seem that domestication has not resulted in a substantial genomewide loss of diversity, as would be expected had a severe population bottleneck occurred. This is important, because it contradicts the assertion that animal domestication began from a small number of individuals in a restricted geographic region. That is still a possible scenario for the very earliest phases of domestication, but if so, our data imply that subsequent crossing with the wild ancestor (in the first thousand years until more developed breeds were established) restored this diversity. Nevertheless, extensive diversity is consistent with the ongoing improvements in agricultural traits that have been achieved over the last 80 years, in layer and broiler lines.

The most important application for this SNP map will be in analysis of QTLs and other genetic traits. Although the density of markers far exceeds what is needed for initial mapping, the principal challenge is not in the detection of linkage but in the identification of genes underlying QTLs. By itself, our SNP map is not adequate. It must be combined with novel strategies and novel resources (like mapping populations specifically designed for fine mapping). The essential problem is the lack of a one-to-one relationship between genotype and phenotype, as the latter is influenced by multiple genetic and environmental factors. This can be overcome, in experimental and domestic animals, by progeny testing and segregation analysis, which permit detailed characterization of haplotypes associated with different QTL alleles, and may eventually lead to the identification of the underlying causative mutations. This SNP map will facilitate fine mapping.

As an example, the major Growth1 QTL on GGA1 explains about one third of the difference between RJF and White Leghorn in adult body weight and egg weight. Initial mapping assigned this locus a ∼20 cM confidence interval. Selective back-crossing using sires that have recombinant chromosomes, and QTL analysis using subsequent intercross generations, are currently employed to refine the localization to a few cM, expected to be less than ∼1 Mb. This also establishes a collection of chromosomes of known QTL status. Our SNP map can then be used for haplotype analysis, assuming that the White Leghorns share a chromosomal segment – identical by descent (IBD) – with the causative mutation. The small haplotype blocks detected in this study underscore the need for a larger number of SNPs to identify such IBD segments. Although these small blocks may require greater marker density and more recombinants to identify the causative haplotype, less effort will be required to resolve the actual QTL alleles once the haplotype is found.

Materials and methods

Our broiler and layer lines are from European breeds with dramatic differences in meat and egg production traits. This specialization started only during the first half of the 20th century. The sequenced male White Cornish-type broiler is from a closed breeding population commonly used in the production of commercial meat-type hybrids (Aviagen, Newbridge, Scotland); effective population size is about 800. The female White Leghorn layer is from a closed line developed at Swedish University of Agricultural Sciences; its effective population size has been 60 to 80 birds for the past 30 years. The Chinese Silkie is used in meat/egg production and traditional Chinese medicine. Selection intensity has been low, and the sequenced female is from a large outbred population.

DNA was extracted from erythrocytes of a single bird, sheared by sonication, and size fractionated via agarose gels. Fragments of 3-kb size were ligated to SmaI-cut blunt-ended pUC18 plasmid vectors. Single colonies were grown overnight, and plasmids were extracted by an alkaline lysis protocol. Sequences were read from both ends of the insert, with vector primers and Amersham MegaBACE 1000 capillary sequencers. Roughly one million reads were generated for each bird. For broiler, layer, and Silkie, we got a total of 841,790, 841,555, and 870,556 successful reads, whose Q20 lengths add to 380,729,199-bp, 372,263,344-bp, and 397,831,117-bp, respectively.

To minimize sequencing errors, we use the Phred quality, Q,. This is related to the single base error rate by the equation: -10×log10(Q). We use more stringent thresholds than normal, with Q>25 for the variant site and Q>20 in both flanking 5-bp regions. For an insertion-deletion (indel), the variant site in the shorter allele is given the quality of its two flanking bases. We originally found many artifactual deletions relative to RJF, which upon a closer examination of the sequence reads were due to doublet peaks that got called as singlet peaks. This is an unavoidable flaw of the base caller software. Hence, we raised the indel thresholds to Q30 and Q25. We must still advise caution, and to that end, indels in simple repeats are flagged and none are counted in our summary tables.

Paralog confusion is detected in the course of the genome level BlastN search that determines where the read is supposed to go. Once this is known, the detailed alignments are done within CrossMatch. Analysis of the RJF genome shows that recent segmental duplications typically agree to 2%. When the best and second best BlastN hits were more than 2% apart, and the best hit was not to a known segmental duplication, the best hit was taken. When either rule was violated, clone-end pairs information was used to resolve the ambiguity. Every alignment had to incorporate 80% of the read. Mapped back to the RJF genome, the amount of usable data for broiler, layer, and Silkie covered 190,513,980-bp, 165,154,746-bp, and 210,214,479-bp respectively.

Polymorphism rates are normalized to the length of the sequence on which we can detect SNPs. To correct for heterozygosity within a line, we compute nucleotide diversity using the approximation:π=K/n1i=1Li, where K is the number of variant sites found by sequencing n chromosomes in a region of length L. When comparing RJF to one of the 3 domestic lines, n can only be 2 or 3, and it is a stochastic variable, because there is a 50% chance that any two overlapping reads are from the same chromosome. When there are m overlapping reads, the denominator is L2m1(1+(2m11)(1+12)). We then sum over all possible regions, with different L and m for each region, to get what we call the “effective length”. Similar considerations are used to compute SNP rates within a line, except that n is 1 or 2, and as a result, the denominator becomes L2m1(2m11).

We compute gene context relative to 5 different data sets. The first 3 are based on experimentally derived genes and the last 2 are based on computer annotations. Riken1 is a set of 1758 full-length cDNAs taken from bursal B-cells of a two week old CB inbred. GenBank refers to 1178 chicken genes with “complete CDS” designation, downloaded as version 2003-12-15. BBSRC is a set of 1184 cDNAs, taken from a larger group of 18,034 cDNAs, which are full-length using a TBlastX mapping to vertebrate Refseq and BlastX mapping to SWALL. Merging all 3 data sets, we have 3868 non-redundant genes. For the detailed gene models, we do a genome level search in BLAT and use SIM4 to compute the exon-intron boundaries. The last two data sets are for 995 chicken orthologs of human disease genes and 17,709 non-redundant Ensembl genes.

Additional details are in Supplement M (Methods).

Supplementary Material




Beijing Institute of Genomics of Chinese Academy of Sciences Gallus gallus SNP discovery and analysis was supported by Chinese Academy of Sciences (KSCX2-SW-223), State Development Planning Commission, Ministry of Science and Technology (2002AA104250; 2004AA231050; 2001AA231061; 2001AA231101), National Natural Science Foundation of China (30200163; 90208019), Beijing Municipal Government, Zhejiang Provincial Government, Hangzhou Municipal Government, Zhejiang University, and China National Grid. Some equipment and reagents were provided by Wellcome Trust and Sanger Institute of the UK. Recent segmental duplications were analyzed by G. Cheng and E.E. Eichler. Riken1 cDNAs were provided by R. Caldwell and J.M. Buerstedde. Noncoding conserved motifs were analyzed by J. Taylor and W. Miller. Washington University School of Medicine Gallus gallus sequence generation was supported by National Human Genome Research Institute. Uppsala University HE was supported by Swedish Research Council, Knut and Alice Wallenberg Foundation, and Royal Academy of Sciences. LA was supported by Wallenberg Consortium North, Foundation for Strategic Research, and Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning. Institute for Animal Health PK, NB, JRY, and JK were supported by BBSRC. Iowa State University SJL was supported by Hatch Act and State of Iowa. Skeletal data for ISU resource population was collected by C. Ashwell and A. Mitchell. Roslin Institute PMH, AL, DJK, and DWB were supported by BBSRC. SNP genotyping was partially funded by Cobb-Vantress. USDA-ARS Avian Disease and Oncology Laboratory J. Kenyon and N. Evenson provided technical assistance. University of Oxford CPP was supported by UK Medical Research Council. University of Manchester Institute of Science and Technology SJH was supported by BBSRC. University of Sheffield SAW was supported by BBSRC.

We dedicate this paper to Nat Bumstead, who died during preparation of the manuscript. Nat was recognised as a major figure in researching the genetics of disease resistance in poultry. He worked tirelessly to realise the sequence of the chicken genome, which led in part to this consortium.

International Chicken Polymorphism Map Consortium

†International Chicken Polymorphism Map Consortium

(Group contributions are listed by their order of appearance in the manuscript)

Polymorphism discovery and analysis:Beijing Institute of Genomics of Chinese Academy of Sciences: Gane Ka-Shu Wong13,*, Bin Liu1,*, Jun Wang1,2,*, Yong Zhang1,4,*, Xu Yang1,*, Zengjin Zhang1, Qingshun Meng1, Jun Zhou1, Dawei Li1, Jingjing Zhang1, Peixiang Ni1, Songgang Li1,4, Longhua Ran5, Heng Li1,6, Jianguo Zhang1, Ruiqiang Li1, Shengting Li1, Hongkun Zheng1, Wei Lin1, Guangyuan Li1, Xiaoling Wang1, Wenming Zhao1, Jun Li1, Chen Ye1, Mingtao Dai1, Jue Ruan1, Yan Zhou2, Yuanzhe Li1, Ximiao He1, Yunze Zhang1, Jing Wang1,4, Xiangang Huang1, Wei Tong1, Jie Chen1, Jia Ye1,2, Chen Chen1, Ning Wei1, Guoqing Li1, Le Dong1, Fengdi Lan1, Yongqiao Sun1, Zhenpeng Zhang1, Zheng Yang1, Yingpu Yu2, Yanqing Huang1, Dandan He1, Yan Xi1, Dong Wei1, Qiuhui Qi1, Wenjie Li1, Jianping Shi1, Miaoheng Wang1, Fei Xie1, Jianjun Wang1, Xiaowei Zhang1, Pei Wang1, Yiqiang Zhao7, Ning Li7, Ning Yang7, Wei Dong1, Songnian Hu1, Changqing Zeng1, Weimou Zheng1,6, Bailin Hao1,6

Genome sequence of Red Jungle Fowl:Washington University School of Medicine: LaDeana W. Hillier8, Shiaw-Pyng Yang8, Wesley C. Warren8, Richard K. Wilson8

Molecular evolution:Uppsala University: Mikael Brandström9, Hans Ellegren9

Population genotyping, BAC sequences and haplotypes:Wageningen University: Richard P.M.A. Crooijmans10, Jan J. van der Poel10, Henk Bovenhuis10, Martien A.M. Groenen10; Lawrence Livermore National Laboratory: Ivan Ovcharenko11,12, Laurie Gordon11,13, Lisa Stubbs11; DOE Joint Genome Institute: Susan Lucas13, Tijana Glavina13, Andrea Aerts13

Examples of application to complex traits:Institute for Animal Health: Pete Kaiser14, Lisa Rothwell14, John R. Young14, Sally Rogers14, Brian A. Walker14, Andy van Hateren14, Jim Kaufman14, Nat Bumstead14; Iowa State University: Susan J. Lamont15, Huaijun Zhou15; Roslin Institute: Paul M. Hocking16, David Morrice16, Dirk-Jan de Koning16, Andy Law16, Neil Bartley16, David W. Burt16; USDA-ARS Avian Disease and Oncology Laboratory: Henry Hunt17, Hans H. Cheng17

Domestication and selection:Uppsala University: Ulrika Gunnarsson18, Per Wahlberg18, Leif Andersson18,19,; Karolinska Institutet: Ellen Kindlund20, Martti T. Tammi20,21, Björn Andersson20

Human disease genes:University of Oxford: Caleb Webber22, Chris P. Ponting22

EST-based SNP data:University of Manchester Institute of Science and Technology: Ian M. Overton23, Paul E Boardman23, Haizhou Tang23, Simon J. Hubbard23; University of Sheffield: Stuart A Wilson24

Scientific management:Beijing Institute of Genomics of Chinese Academy of Sciences: Jun Yu1,2, Jian Wang1,2, HuanMing Yang1,2,

1Beijing Institute of Genomics of Chinese Academy of Sciences, Beijing Genomics Institute, Beijing Proteomics Institute, Beijing 101300, China

2James D. Watson Institute of Genome Sciences of Zhejiang University, Hangzhou Genomics Institute, Key Laboratory of Bioinformatics of Zhejiang Province, Hangzhou 310007, China

3UW Genome Center, Department of Medicine, University of Washington, Seattle, WA 98195, USA

4College of Life Sciences, Peking University, Beijing 100871, China

5Beijing North Computation Center, Beijing 100091, China

6The Institute of Theoretical Physics Chinese Academy of Sciences, Beijing 100080, China

7China Agricultural University, Beijing 100094, China

8Genome Sequencing Center, Washington University School of Medicine, Campus Box 8501, 4444 Forest Park Avenue, St. Louis, MO 63108, USA

9Department of Evolutionary Biology, Evolutionary Biology Centre, Uppsala University, Norbyvägen 18D, SE-752 34 Uppsala, Sweden

10Animal Breeding and Genetics Group, Wageningen University, Marijkewg 40, 6709 PG Wageningen, The Netherlands

11Genome Biology Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA

12Energy, Environment, Biology and Institutional Computing, Lawrence Livermore National Laboratory

13DOE Joint Genome Institute, Walnut Creek, CA 94598, USA

14Institute for Animal Health, Compton, Berkshire RG20 7NN, UK

15Department of Animal Science, Iowa State Univeristy, Ames, IA 50011, USA

16Roslin Institute (Edinburgh), Roslin, Midlothian EH25 9PS, UK

17USDA-ARS Avian Disease and Oncology Laboratory, 3606 E. Mount Hope Rd., East Lansing, MI 48823, USA

18Department of Medical Biochemistry and Microbiology, Uppsala University, Box 597, SE-751 24 Uppsala, Sweden

19Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, SE-751 24 Uppsala, Sweden

20Center for Genomics and Bioinformatics, Karolinska Institutet, SE-171 77 Stockholm, Sweden

21Departments of Biological Sciences and Biochemistry, National University of Singapore, Singapore

22MRC Functional Genetics Unit, University of Oxford, Department of Human Anatomy and Genetics, South Parks Road, Oxford OX1 3QX, UK

23Department of Biomolecular Sciences, University of Manchester Institute of Science and Technology, PO Box 88, Manchester M60 1QD, UK

24Department of Molecular Biology and Biotechnology, University of Sheffield, Firth Court, Western Bank, Sheffield S10 2TN, UK

Corresponding authors: Gane Ka-Shu Wong nc.gro.scimoneg@wskg, Leif Andersson es.uu.mibmi@nossredna.fiel, HuanMing Yang nc.gro.scimoneg@gnayh.

*These authors contributed equally to this work.


The individual SNPs were deposited at GenBank/dbSNP with submitted SNP (ss) number ranges: 24821291 to 24922086, 24922088 to 26161960, 26161962 to 28446123, and 28452569 to 28452598. They may also be found at http://chicken.genomics.org.cn, the UCSC genome browser, and the Ensembl genome browser. Access to raw sequencing traces is being provided through the NCBI Trace Archive.


1. The Chicken Genome Sequencing Consortium The sequence of the chicken genome, Gallus gallus. Nature. (companion paper) []
2. Pisenti JM, et al. Avian genetic resources at risk: an assessment and proposal for conservation of genetic stocks in the USA and Canada. Avian Poult. Biol. Rev. 2001;12:1–102. []
3. Dodgson JB, Romanov MN. Use of chicken models for the analysis of human disease. In: Dracopoli NC, et al., editors. Current Protocols in Human Genetics. Hoboken: John Wiley & Sons; 2004. pp. 15.5.1–11. []
4. Nicholas FW. Online Mendelian Inheritance in Animals (OMIA): a comparative knowledgebase of genetic disorders and other familial traits in non-laboratory animals. Nucleic Acids Res. 2003;31:275–277. http://www.angis.org.au/Databases/BIRX/omia. [PMC free article] [PubMed] []
5. ChickAce database from the Animal Science Group of the Wageningen University and Research Center; https://acedb.asg.wur.nl. []
6. Groenen MA, et al. A consensus linkage map of the chicken genome. Genome Res. 2000;10:137–147. [PMC free article] [PubMed] []
7. Groenen MA, Crooijmans RP. Structural genomics: integrating linkage, physical and sequence maps. In: Muir WM, Aggrey SE, editors. Poultry Genetics, Breeding and Biotechnology. Wallingford: CABI Publishing; 2003. pp. 497–536. []
8. Vignal A, Milan D, SanCristobal M, Eggen A. A review on SNP and other types of molecular markers and their use in animal genetics. Genet. Sel. Evol. 2002;34:275–305. [PMC free article] [PubMed] []
9. Andersson L, Georges M. Domestic-animal genomics: deciphering the genetics of complex traits. Nat. Rev. Genet. 2004;5:202–212. [PubMed] []
10. There are only 20-kb of aligned sequence on GGA16, and if we were to remove it, the total SNP rate would only change by 0.02%.
11. Begun DJ, Aquadro CF. Levels of naturally occurring DNA polymorphism correlate with recombination rates in D. melanogaster. Nature. 1992;356:519–520. [PubMed] []
12. Nachman MW. Single nucleotide polymorphisms and recombination rate in humans. Trends Genet. 2001;(9):481–485. [PubMed] []
13. Crooijmans RP, Vrebalov J, Dijkhof RJ, van der Poel JJ, Groenen MA. Two-dimensional screening of the Wageningen chicken BAC library. Mamm. Genome. 2000;11:360–363. [PubMed] []
14. Sundstrom H, Webster MT, Ellegren H. Reduced variation on the chicken Z chromosome. Genetics. 2004;167:377–385. [PMC free article] [PubMed] []
15. Ikeobi CO, et al. Quantitative trait loci for muscling in a broiler layer cross. Livest. Prod. Sci. 2004;87:143–151. []
16. Sewalem A, et al. Mapping of quantitative trait loci for body weight at three, six, and nine weeks of age in a broiler layer cross. Poult. Sci. 2002;81:1775–1781. [PubMed] []
17. Li H, et al. Chicken quantitative trait loci for growth and body composition associated with transforming growth factor-β genes. Poult. Sci. 2003;82:347–356. [PubMed] []
18. Zhou H, Li H, Lamont SJ. Genetic markers associated with antibody response kinetics in adult chickens. Poult. Sci. 2003;82:699–708. [PubMed] []
19. Gallagher G, Eskdale J, Bidwell JL. Cytokine genetics – polymorphisms, functional variations and disease associations. In: Thomson AW, Lotze MT, editors. The Cytokine Handbook. 4th Edition London: Academic Press; 2003. pp. 19–55. []
20. Bumstead N, et al. EU Project FAIR3 PL96-1502 New Molecular Approaches for Improved Poultry Vaccines. Compton: Institute for Animal Health; 2000. []
21. Maynard-Smith J, Haigh J. The hitch-hiking effect of a favourable gene. Genet. Res. 1974;23:23–35. [PubMed] []
22. Van Laere AS, et al. regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig. Nature. 2003;425:832–836. [PubMed] []
23. Olson MV. When less is more: gene loss as an engine of evolutionary change. Am. J. Hum. Genet. 1999;64:18–23. [PMC free article] [PubMed] []
24. Grobet L, et al. A deletion in the bovine myostatin gene causes the double- muscled phenotype in cattle. Nat. Genet. 1997;17:71–74. [PubMed] []
25. Ng PC, Henikoff S. Predicting deleterious amino acid substitutions. Genome Res. 2001;11:863–874. http://blocks.fhcrc.org/sift/SIFT.html. [PMC free article] [PubMed] []
26. Gilbert-Dussardier B, et al. Partial duplication [dup. TCAC (178)] and novel point mutations (T125M, G188R, A209V, and H302L) of the ornithine transcarbamylase gene in congenital hyperammonemia. Hum. Mutat. 1996;8:74–76. [PubMed] []
27. Tamir H, Ratner S. Enzymes of arginine metabolism in chicks. Arch. Biochem. Biophys. 1963;102:249–258. [PubMed] []
28. Sachidanandam R, et al. A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms. Nature. 2001;409:928–933. [PubMed] []
29. Parker HG, et al. Genetic structure of the purebred domestic dog. Science. 2004;304:1160–1164. [PubMed] []
30. Yu N, Jensen-Seaman MI, Chemnick L, Ryder O, Li WH. Nucleotide diversity in gorillas. Genetics. 2004;166:1375–1383. [PMC free article] [PubMed] []
31. Lindblad-Toh K, et al. Large-scale discovery and genotyping of single- nucleotide polymorphisms in the mouse. Nat. Genet. 2000;24:381–386. [PubMed] []
32. Axelsson E, Smith NG, Sundstrom H, Berlin S, Ellegren H. Male-biased mutation rate and divergence in autosomal, Z-linked and W-linked introns of chicken and turkey. Mol. Biol. Evol. 2004;21:1538–1547. [PubMed] []
33. Mason IL, editor. Evolution of Domesticated Animals. New York: Longman, Inc.; 1984. []
34. Arthur JA, Albers GA. Industrial perspective on problems and issues associated with poultry breeding. In: Muir WM, Aggrey SE, editors. Poultry Genetics, Breeding and Biotechnology. Wallingford: CABI Publishing; 2003. pp. 1–12. []
35. Kerje S, et al. The twofold difference in adult size between the red junglefowl and White Leghorn chickens is largely explained by a limited number of QTLs. Anim. Genet. 2003;34:264–274. [PubMed] []
36. Crawford RD, editor. Poultry Breeding and Genetics. New York: Elsevier Science; 1990. []
37. Liljedahl LE, Kolstad N, Sorensen P, Maijala K. Scandinavian selection and cross-breeding experiment with laying hens. 1. Background and general outline. Acta Agricult. Scand. 1979;29:273–285. []
38. Niu D, et al. The origin and genetic diversity of Chinese native chicken breeds. Biochem. Genet. 2002;40:163–174. [PubMed] []
39. Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998;8:175–185. [PubMed] []
40. Ewing B, Green P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998;8:186–194. [PubMed] []
41. Altshuler D, et al. A SNP map of the human genome generated by reduced representation shotgun sequencing. Nature. 2000;407:513–516. [PubMed] []
42. Green P. CrossMatch is the underlying alignment tool for the Phrap assembly software. at http://www.phrap.org.
43. Cargill M, et al. Characterization of single-nucleotide polymorphisms in coding regions of human genes. Nat. Genet. 1999;22:231–238. [PubMed] []
44. Caldwell R, et al. A large collection of bursal full-length cDNA sequences to facilitate gene function analysis. Genome Biol. (companion issue) []
45. Hubbard SJ, et al. Transcriptome analysis for the chicken based on 19,626 finished cDNA sequences and 485,337 expressed sequence tags. Genome Res. (companion issue) [PMC free article] [PubMed] []
46. Kent WJ. BLAT – the BLAST-like alignment tool. Genome Res. 2002;12:656–664. http://www.genome.ucsc.edu/cgi-bin/hgBlat. [PMC free article] [PubMed] []
47. Florea L, Hartzell G, Zhang Z, Rubin GM, Miller W. A computer program for aligning a cDNA sequence with a genomic DNA sequence. Genome Res. 1998;8:967–974. http://globin.cse.psu.edu/html/docs/sim4.html. [PMC free article] [PubMed] []
48. Wang J, et al. ChickVD: a sequence variation database for the chicken genome. Nucleic Acids Res. 2005 Jan; [PMC free article] [PubMed] []


Please enter your comment!
Please enter your name here