- Research article
A genome-wide scan of selective sweeps in two broiler chicken lines divergently selected for abdominal fat content
BMC Genomicsvolume 13, Article number: 704 (2012)
Genomic regions controlling abdominal fatness (AF) were studied in the Northeast Agricultural University broiler line divergently selected for AF. In this study, the chicken 60KSNP chip and extended haplotype homozygosity (EHH) test were used to detect genome-wide signatures of AF.
A total of 5357 and 5593 core regions were detected in the lean and fat lines, and 51 and 57 reached a significant level (P<0.01), respectively. A number of genes in the significant core regions, including RB1, BBS7, MAOA, MAOB, EHBP1, LRP2BP, LRP1B, MYO7A, MYO9A and PRPSAP1, were detected. These genes may be important for AF deposition in chickens.
We provide a genome-wide map of selection signatures in the chicken genome, and make a contribution to the better understanding the mechanisms of selection for AF content in chickens. The selection for low AF in commercial breeding using this information will accelerate the breeding progress.
The linkage disequilibrium (LD) is important in livestock genetics for its key role in genomic selection  and detecting the causal mutations of economically important traits [2–6]. Based on the LD information, there are two main methods to detect genes underlying phenotypic variation, including one from phenotype to genome and another one from genome to phenotype. The first method is performed by targeting particular candidate genes or by quantitative trait loci (QTL) mapping and positional cloning of QTL. In the second method, patterns of LD in populations that are incompatible with the hypothesis of genetic neutrality are identified, and these patterns are selection signatures . The aim of the second method is to identify artificial selections by statistically evaluating the genomic data .
Allele frequencies underlying selection are expected to change. A neutral mutation will take many generations until the mutated allele reaches a high or low population frequency. In this case, the LD between the mutation and its neighboring loci will be degraded because of the recombination in every generation . The frequency of a novel mutation will increase or decrease more rapidly than the neutral mutation because it is underlying artificial selection, so that the surrounding conserved haplotype was long [9, 10]. This is the background of the extended haplotype homozygosity (EHH) statistic method used to detect selection signatures . There are also many other methods to detect selective sweeps from DNA sequence data, including the Tajima’s D and Fay and Wu’s H-test  for selected mutations, measuring large allele-frequency differences among populations by FST, and the integrated Haplotype Score (iHS) , which is an extension of the EHH statistic . Among these methods, the EHH test is particularly useful [7, 11]. The EHH test is used to detect artificial selections according to the characteristics of haplotypes within a single population, and do not require the genotype of the ancestor . Furthermore, the EHH test is less sensitive to ascertainment bias than other approaches, so it was designed to work with SNP rather than sequencing data [7, 16].
The broilers used in this study were selected for eleven generations and genomic regions controlling AF deposition are expected to exhibit signatures of selective sweep. The aim of this study was to identify the selection signatures underlying the artificial selection for AF in chicken and to investigate the genes important for AF deposition.
All animal work was conducted according to the guidelines for the care and use of experimental animals established by the Ministry of Science and Technology of the People’s Republic of China (Approval number: 2006–398) and approved by the Laboratory Animal Management Committee of Northeast Agricultural University.
DNA samples and data preparation
Broilers used in this study were from two Northeast Agricultural University broiler lines divergently selected for AF content (NEAUHLF). The two lines have been selected since 1996 using AF percentage (%AFW or AFP) and plasma very low-density lipoprotein (VLDL) concentration as selection criteria . The two lines were selected for 11 generations and the AFP changes over the 11 generations are shown in Figure 1. A total of 475 individuals from generation 11 of NEAUHLF were used in this study.
Genotyping was carried out using the Illumina chicken 60K SNP chip containing a total of 57636 SNPs. Markers were filtered to exclude loci with unknown positions, monomorphic loci and loci with a minor allele frequency <0.05.
The haplotype and LD analysis
The fastPHASE  (http://depts.washington.edu/fphase/download/) was used to reconstruct the haplotypes for every chromosome using the default parameters. The reconstructed haplotypes were inserted into HAPLOVIEW v4.1  to estimate LD statistics based on pairwise r2 and to construct the blocking pattern in the candidate regions of interest to enable selection signature analysis.
The EHH test
The “core region” was defined as the region in the genome characterized by the strong LD among SNPs involving a set of “core haplotypes” . The Sweep v.1.1 (http://www.soft82.com/get/download/windows/sweep/) was used to identify the core regions . The algorithm defined a pair of SNPs to be in strong LD if the upper 95% confidence bound of D’ is between 0.70 and 0.98 . The program was set to select core regions with at least two SNPs. EHH was defined as the probability that two randomly chosen haplotypes carrying the candidate core haplotype were homozygous for the entire interval spanning the core region to a given locus . The EHH test  was based on one of the core haplotype vs. other haplotypes in the same position. The “Relative Extended Haplotype Homozygosity” (REHH) statistic corrects EHH for the variability in recombination rates . It was computed by EHHt / ; with defined as the decay of EHH on all other core haplotypes combined. The REHH value was used in the current study to determine the selection signatures. To determine the significance of REHH values, the haplotypes were ordered into 20 bins according to their frequencies . The REHH values of each haplotype in a candidate region were compared with all equally frequent haplotypes and the P-values were obtained . The significant selection signatures were defined as P<0.01.
Markers and core haplotypes
A total of 43034 SNPs on 28 autosomes in chickens were included in the selection signature analysis (Table 1). These markers covered 950.68 Mb of the genome, with an average of 22.09 kb between adjacent markers.
For the SNPs analyzed in this study, the average minor allele frequency was 0.29 ± 0.13. A summary of genome-wide markers and core haplotype distribution in the data set is shown in Table 1. A total of 5357 and 5593 core regions spanning 549523.91 kb and 480784.79 kb of the genome, respectively, in the lean and fat lines were detected (Table 1). Mean core region length was estimated as 102.58±37.24 kb and 85.96±26.65 kb, with a maximum of 2288.64 kb and 2191.34 kb in the lean and fat lines, respectively (Table 1). Chromosome 1 was the largest chromosome in chickens, and it had the largest haplotypic structures in the genome, which covered 110644.43 kb and 105728.03 kb in the lean and fat lines, respectively. For each chromosome, the proportion of length covered by core regions vs. total length, as well as the number of SNPs forming core regions vs. the total number of SNPs, are shown in Table 1. The distribution of the size of core regions is shown in Figure 2. Overall, 25069 and 22180 SNPs in the lean and fat lines, respectively, participated in forming core regions, with a range of 2 to 19 SNPs per core.
Whole genome selection signatures
For all 5357 and 5593 core regions in the lean and fat lines, respectively, a total of 44822 and 46775 EHH tests, with an average of 8.37 and 8.36 tests per core region, were calculated. To find outlying core haplotypes, we calculated REHH at 1 Mb distances both on the upstream and downstream sides. Figure 3 shows the distribution of REHH values vs. haplotype frequencies in the lean and fat lines, respectively. Corresponding P-values are indicated by different colored symbols. The –log10 of the P-values associated with REHH against the chromosomal position was plotted to visualize the chromosomal distribution of outlying core haplotypes with frequency <25% (Figure 4). The results indicated that these selection signals were not uniformly distributed across all chromosomes, with a substantial overrepresentation on chromosomes 1, 2, 3 and 4.
The genome-wide statistics of the selection signature test, including the number of tests and outlying core haplotypes for each chromosome, are presented in Table 2. Of 16677 and 18346 tests on core haplotypes with frequency ≥0.25, there were 51 and 57 tests with P<0.01 in the lean and fat lines, respectively. There were 153 and 251 tests with P<0.05 in the lean and fat lines, respectively.
The conformity of the distribution of Tukey’s outliers was examined, with outlying core haplotypes defined at the threshold level of 0.01. Figure 5 displays box plots of the distribution of –log10 (P-values) within each bin of core haplotype frequency. The results indicated that the extreme outliers appear in the small haplotype frequencies bins.
Mapping selection signatures to genes
A summary of statistics for 51 and 57 positively selected core regions with P<0.01 of the REHH tests in the lean and fat lines, respectively, is presented in Table 3. Corresponding genes were identified by aligning the core positions with the chicken genome sequence (Table 3). The full genes names were from Ensembl online (http://www.ensembl.org/index.html). A total of 66 and 46 genes in the core regions were detected in the lean and fat lines, respectively, including RB1 (retinoblastoma 1), BBS7 (Bardet-Biedl syndrome 7), MAOA (monoamine oxidase A), MAOB (monoamine oxidase B), EHBP1 (EH domain binding protein 1), LRP2BP (LRP2 binding protein), LRP1B (low-density lipoprotein receptor-related protein 1B), MYO7A (myosin VIIA), MYO9A (myosin IXA) and PRPSAP1 (phosphoribosyl pyrophosphate synthetase-associated protein 1). The haplotype analysis of these genes revealed that the haplotype frequencies were significantly different (P<0.01) between the two lines (Table 4).
Mapping selection signatures to QTLs
The chicken QTL database available online (http://www.animalgenome.org/cgi-bin/QTLdb/GG/index) was explored to identify any overlapping of the core regions with significant REHH P-values (P<0.01) and published QTLs in chickens. The approximate positions of the overlapping QTLs for each core region are listed in Table 5. There were many overlaps between the core regions with significant REHH P-values (P<0.01) and published QTLs for AF content in chickens.
Selective sweep is used to detect genomic regions with reduced variation in allele frequency in any population experiencing divergent selection for specific traits. Here, we determined the feasibility of the selective sweep approach for finding genes important for AF deposition in chickens. The long-range haplotype test was employed, which detects selection signature by measuring the characteristics of haplotypes within the lean and fat lines divergently selected for AF content. There were 5357 and 5593 core regions in the lean and fat lines, respectively. When comparing the average marker spacing with mean core length and number of SNPs forming cores, we revealed that core regions are more likely to appear in regions with higher marker density.
The selection signatures on the whole genome were calculated, and a subset of putative core regions with significant REHH P-values (P<0.01) was identified. The genes in these core regions were detected and 10 genes, including RB1, BBS7, MAOA, MAOB, EHBP1, LRP2BP, LRP1B, MYO7A, MYO9A and PRPSAP1, were important for fatness. Among these 10 important genes, seven genes, including RB1, BBS7, MAOA, MAOB, EHBP1, LRP2BP and LRP1B, were all in the QTL regions reported previously for AF in chickens (Table 5). Although the other three genes, including MYO7A, MYO9A and PRPSAP1, were not in the QTL regions, these genes were also important for the AF deposition.
The known functions of these 10 genes were analyzed and the results indicated that they were likely to be linked with fatness. The RB1 gene regulates the C/EBP-DNA-binding activity during 3T3-L1 adipogenesis and plays a key role in adipocyte differentiation [40, 41].
The BBS7 gene is a member of the Bardet-Biedl syndrome (BBS) family. BBS is a pleiotropic genetic disorder characterized by obesity, photoreceptor degeneration, polydactyly, hypogenitalism, renal abnormalities, and developmental delay . BBS is recognized to be a genetically heterogeneous autosomal recessive disorder mapped to eight loci . Positional cloning and candidate genes identified six BBS genes, including BBS1, BBS2, BBS4, BBS6, BBS7, and BBS8. These BBS genes may be important for obesity.
The MAOA and MAOB are two enzymes important for dopamine production. The dopamine levels influence the risk of obesity and MAOA and MOAB may be implicated in human obesity .
The EHBP1 gene is required for insulin-stimulated GLUT4 movements . Insulin stimulates glucose transport in adipose tissues by recruiting intracellular membrane vesicles containing the glucose transporter GLUT4 to the plasma membrane . The mechanisms involved in the biogenesis of these vesicles and their translocation to the cell surface were studied and the results indicated that EHD1 and EHBP1 are required for perinuclear localization of GLUT4, and the loss of EHBP1 disrupts insulin-regulated GLUT4 recycling in cultured adipocytes . This indicates that the EHBP1 gene may be important in adipocyte differentiation.
The LRP2BP and LRP1B genes are two members of the low-density lipoprotein receptor family that participates in a wide range of physiological processes, including the regulation of lipid metabolism, protection against atherosclerosis, neurodevelopment, and transport of nutrients and vitamins .
The MYO7A and MYO9A are two myosin genes. A spontaneous mutant mouse line, Myo7ash1-6J, was used to study the function of the MYO7A gene, and the result indicated that the mutant male homozygous mice displayed decreased body weight and body fat . The MYO9A gene was in the BBS4 region of chromosome 15q22-q23 , which might be important for obesity.
The PRPSAP1 gene is named as phosphoribosyl pyrophosphate synthetase-associated protein 1. The results of differentially expressed genes associated with insulin resistance indicate that PRPSAP1 gene is associated with percentage of body fat .
The associations of these 10 genes with obesity or lipid metabolism were mainly in humans and mice. Because of the high conservation of these genes between humans, mice and chickens, the 10 genes might also be important for AF deposition in chickens.
Our results provide a genome-wide map of selection signatures in two chicken lines divergently selected for AF content. There were 51 and 57 core regions showing significant P-values (P<0.01) of selection signatures in the lean and fat lines, respectively. In these core regions there were a number of important genes, including RB1, BBS7, MAOA, MAOB, EHBP1, LRP2BP, LRP1B, MYO7A, MYO9A and PRPSAP1. These genes are important for AF deposition in chickens.
Abdominal fat percentage
Abdominal fat weight
Extended haplotype homozygosity
Integrated Haplotype Score
Northeast Agricultural University
Northeast Agricultural University broiler lines divergently selected for abdominal fat content
National Research Council
Quantitative trait loci
Relative Extended Haplotype Homozygosity
Single nucleotide polymorphism
Very low-density lipoprotein.
Meuwissen TH, Hayes BJ, Goddard ME: Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001, 157 (4): 1819-1829.
Pollinger JP, Bustamante CD, Fledel-Alon A, Schmutz S, Gray MM, Wayne RK: Selective sweep mapping of genes with large phenotypic effects. Genome Res. 2005, 15 (12): 1809-1819.
Daetwyler HD, Schenkel FS, Sargolzaei M, Robinson JA: A genome scan to detect quantitative trait loci for economically important traits in Holstein cattle using two methods and a dense single nucleotide polymorphism map. J Dairy Sci. 2008, 91 (8): 3225-3236.
Hayes BJ, Lien S, Nilsen H, Olsen HG, Berg P, Maceachern S, Potter S, Meuwissen TH: The origin of selection signatures on bovine chromosome 6. Anim Genet. 2008, 39 (2): 105-111.
Hayes BJ, Chamberlain AJ, Maceachern S, Savin K, McPartlan H, MacLeod I, Sethuraman L, Goddard ME: A genome map of divergent artificial selection between Bos taurus dairy cattle and Bos taurus beef cattle. Anim Genet. 2009, 40 (2): 176-184.
Prasad A, Schnabel RD, McKay SD, Murdoch B, Stothard P, Kolbehdari D, Wang Z, Taylor JF, Moore SS: Linkage disequilibrium and signatures of selection on chromosomes 19 and 29 in beef and dairy cattle. Anim Genet. 2008, 39 (6): 597-605.
Qanbari S, Pimentel EC, Tetens J, Thaller G, Lichtner P, Sharifi AR, Simianer H: A genome-wide scan for signatures of recent selection in Holstein cattle. Anim Genet. 2010, 41 (4): 377-389.
Kimura M: The Neutral Theory of Molecular Evolution. 1983, Cambridge, MA: Cambridge University Press
Smith JM, Haigh J: The hitch-hiking effect of a favourable gene. Genet Res. 1974, 23 (1): 23-35.
Nielsen R: Molecular signatures of natural selection. Annu Rev Genet. 2005, 39: 197-218.
Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, Gabriel SB, Platko JV, Patterson NJ, McDonald GJ, Ackerman HC, Campbell SJ, Altshuler D, Cooper R, Kwiatkowski D, Ward R, Lander ES: Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002, 419 (6909): 832-837.
Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.
Fay JC, Wu CI: Hitchhiking under positive Darwinian selection. Genetics. 2000, 155 (3): 1405-1413.
Akey JM, Zhang G, Zhang K, Jin L, Shriver MD: Interrogating a high-density SNP map for signatures of natural selection. Genome Res. 2002, 12 (12): 1805-1814.
Voight BF, Kudaravalli S, Wen X, Pritchard JK: A map of recent positive selection in the human genome. PLoS Biol. 2006, 4 (3): e72-
Tang K, Thornton KR, Stoneking M: A new approach for using genome scans to detect recent positive selection in the human genome. PLoS Biol. 2007, 5 (7): e171-
Guo L, Sun B, Shang Z, Leng L, Wang Y, Wang N, Li H: Comparison of adipose tissue cellularity in chicken lines divergently selected for fatness. Poult Sci. 2011, 90 (9): 2024-2034.
Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006, 78 (4): 629-644.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265.
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 (5576): 2225-2229.
Jennen DG, Vereijken AL, Bovenhuis H, Crooijmans RM, van der Poel JJ, Groenen MA: Confirmation of quantitative trait loci affecting fatness in chickens. Genet Sel Evol. 2005, 37 (2): 215-228.
Ikeobi CO, Woolliams JA, Morrice DR, Law A, Windsor D, Burt DW, Hocking PM: Quantitative trait loci affecting fatness in the chicken. Anim Genet. 2002, 33 (6): 428-435.
Liu X, Zhang H, Li H, Li N, Zhang Y, Zhang Q, Wang S, Wang Q, Wang H: Fine-mapping quantitative trait loci for body weight and abdominal fat traits: effects of marker density and sample size. Poult Sci. 2008, 87 (7): 1314-1319.
McElroy JP, Kim JJ, Harry DE, Brown SR, Dekkers JC, Lamont SJ: Identification of trait loci affecting white meat percentage and other growth and carcass traits in commercial broiler chickens. Poult Sci. 2006, 85 (4): 593-605.
Nadaf J, Pitel F, Gilbert H, Duclos MJ, Vignoles F, Beaumont C, Vignal A, Porter TE, Cogburn LA, Aggrey SE, Simon J, Le Bihan-Duval E: QTL for several metabolic traits map to loci controlling growth and body composition in an F2 intercross between high- and low-growth chicken lines. Physiol Genomics. 2009, 38 (3): 241-249.
Lagarrigue S, Pitel F, Carré W, Abasht B, Le Roy P, Neau A, Amigues Y, Sourdioux M, Simon J, Cogburn L, Aggrey S, Leclercq B, Vignal A, Douaire M: Mapping quantitative trait loci affecting fatness and breast muscle weight in meat-type chicken lines divergently selected on abdominal fatness. Genet Sel Evol. 2006, 38 (1): 85-97.
Park HB, Jacobsson L, Wahlberg P, Siegel PB, Andersson L: QTL analysis of body composition and metabolic traits in an intercross between chicken lines divergently selected for growth. Physiol Genomics. 2006, 25 (2): 216-223.
Campos RL, Nones K, Ledur MC, Moura AS, Pinto LF, Ambo M, Boschiero C, Ruy DC, Baron EE, Ninov K, Altenhofen CA, Silva RA, Rosário MF, Burt DW, Coutinho LL: Quantitative trait loci associated with fatness in a broiler-layer cross. Anim Genet. 2009, 40 (5): 729-736.
Atzmon G, Blum S, Feldman M, Lavi U, Hillel J: Detection of agriculturally important QTLs in chickens and analysis of the factors affecting genotyping strategy. Cytogenet Genome Res. 2007, 117 (1–4): 327-337.
Jennen DG, Vereijken AL, Bovenhuis H, Crooijmans RP, Veenendaal A, van der Poel JJ, Groenen MA: Detection and localization of quantitative trait loci affecting fatness in broilers. Poult Sci. 2004, 83 (3): 295-301.
Abasht B, Pitel F, Lagarrigue S, Le Bihan-Duval E, Le Roy P, Demeure O, Vignoles F, Simon J, Cogburn L, Aggrey S, Vignal A, Douaire M: Fatness QTL on chicken chromosome 5 and interaction with sex. Genet Sel Evol. 2006, 38 (3): 297-311.
Ankra-Badu GA, Shriner D, Le Bihan-Duval E, Mignon-Grasteau S, Pitel F, Beaumont C, Duclos MJ, Simon J, Porter TE, Vignal A, Cogburn LA, Allison DB, Yi N, Aggrey SE: Mapping main, epistatic and sex-specific QTL for body composition in a chicken population divergently selected for low or high growth rate. BMC Genomics. 2010, 11: 107-
Atzmon G, Blum S, Feldman M, Cahaner A, Lavi U, Hillel J: QTLs detected in a multigenerational resource chicken population. J Hered. 2008, 99 (5): 528-538.
Zhou H, Deeb N, Evock-Clover CM, Ashwell CM, Lamont SJ: Genome-wide linkage analysis to identify chromosomal regions affecting phenotypic traits in the chicken. II. Body composition. Poult Sci. 2006, 85 (10): 1712-1721.
Nones K, Ledur MC, Ruy DC, Baron EE, Melo CM, Moura AS, Zanella EL, Burt DW, Coutinho LL: Mapping QTLs on chicken chromosome 1 for performance and carcass traits in a broiler x layer cross. Anim Genet. 2006, 37 (2): 95-100.
Uemoto Y, Sato S, Odawara S, Nokata H, Oyamada Y, Taguchi Y, Yanai S, Sasaki O, Takahashi H, Nirasawa K, Kobayashi E: Genetic mapping of quantitative trait loci affecting growth and carcass traits in F2 intercross chickens. Poult Sci. 2009, 88 (3): 477-482.
Le Mignon G, Pitel F, Gilbert H, Le Bihan-Duval E, Vignoles F, Demeure O, Lagarrigue S, Simon J, Cogburn LA, Aggrey SE, Douaire M, Le Roy P: A comprehensive analysis of QTL for abdominal fat and breast muscle weights on chicken chromosome 5 using a multivariate approach. Anim Genet. 2009, 40 (2): 157-164.
Li H, Deeb N, Zhou H, Mitchell AD, Ashwell CM, Lamont SJ: Chicken quantitative trait loci for growth and body composition associated with transforming growth factor-beta genes. Poult Sci. 2003, 82 (3): 347-356.
Hu G, Wang SZ, Wang ZP, Li YM, Li H: Genetic epistasis analysis of 10 peroxisome proliferator-activated receptor γ-correlated genes in broiler lines divergently selected for abdominal fat content. Poult Sci. 2010, 89 (11): 2341-2350.
Fajas L, Egler V, Reiter R, Hansen J, Kristiansen K, Debril MB, Miard S, Auwerx J: The retinoblastoma-histone deacetylase 3 complex inhibits PPARgamma and adipocyte differentiation. Dev Cell. 2002, 3 (6): 903-910.
Cole KA, Harmon AW, Harp JB, Patel YM: Rb regulates C/EBPbeta-DNA-binding activity during 3T3-L1 adipogenesis. Am J Physiol Cell Physiol. 2004, 286 (2): C349-C354.
Sheffield VC: Use of isolated populations in the study of a human obesity syndrome, the Bardet-Biedl syndrome. Pediatr Res. 2004, 55 (6): 908-911.
Need AC, Ahmadi KR, Spector TD, Goldstein DB: Obesity is associated with genetic variants that alter dopamine availability. Ann Hum Genet. 2006, 70 (Pt 3): 293-303.
Guilherme A, Soriano NA, Furcinitti PS, Czech MP: Role of EHD1 and EHBP1 in perinuclear sorting and insulin-regulated GLUT4 recycling in 3T3-L1 adipocytes. J Biol Chem. 2004, 279 (38): 40062-40075.
May P, Woldt E, Matz RL, Boucher P: The LDL receptor-related protein (LRP) family: an old family of proteins with new physiological functions. Ann Med. 2007, 39 (3): 219-228.
Gibson F, Walsh J, Mburu P, Varela A, Brown KA, Antonio M, Beisel KW, Steel KP, Brown SD: A type VII myosin encoded by the mouse deafness gene shaker-1. Nature. 1995, 374 (6517): 62-64.
Gorman SW, Haider NB, Grieshammer U, Swiderski RE, Kim E, Welch JW, Searby C, Leng S, Carmi R, Sheffield VC, Duhl DM: The cloning and developmental expression of unconventional myosin IXA (MYO9A) a gene in the Bardet-Biedl syndrome (BBS4) region at chromosome 15q22-q23. Genomics. 1999, 59 (2): 150-160.
Lee YH, Tokraks S, Pratley RE, Bogardus C, Permana PA: Identification of differentially expressed genes in skeletal muscle of non-diabetic insulin-resistant and insulin-sensitive Pima Indians by differential display PCR. Diabetologia. 2003, 46 (11): 1567-1575.
The authors would like to acknowledge the members of the Poultry Breeding Group of the College of Animal Science and Technology in Northeast Agricultural University for managing the birds and collecting data. This research was supported by the China Agriculture Research System (No. CARS-42), National 863 Project of China (No. 2011AA100301), National 973 Project of China (No. 2009CB941604) and Program for Innovation Research Team in University of Heilongjiang Province (No. 2010td02).
There are no potential competing interests related to this manuscript.
HZ contributed to, conceived and designed the experiments, participated in the interpretation of the data, and drafted and wrote the manuscript. SZW participated in the design of the study and interpretation of the data, and contributed to writing the manuscript. ZPW participated in the design of the study and contributed to writing the manuscript. YD carried out analysis and interpretation of the data. NW participated in the design of the study and contributed to the analysis of the data. XXH participated in the analysis and interpretation of the data. YDZ participated in the design of the study and contributed to the analysis of the data. YXW participated in the design of the study. LL contributed reagents/materials/analysis tools. ZQT contributed reagents/materials/analysis tools. HL co-led the conception and design of the study, participated in the interpretation of the data, and contributed to writing the manuscript. All authors submitted comments on drafts, and read and approved the final manuscript.