Genome-wide association study for the primary feather color trait in a native Chinese duck.

To reveal candidate genes and the molecular genetic mechanism underlying primary feather color trait in ducks, a genome-wide association study (GWAS) for the primary feather color trait was performed based on the genotyping-by-sequencing (GBS) technology for a native Chinese female duck, Longyan Shan-ma ducks. Blood genomic DNA from 314 female Longyan Shan-ma duck were genotyped using GBS technology. A GWAS for the primary feather color trait with genome variations was performed using an univariate linear mixed model based on all SNPs in autosomes. Seven genome-wide significant single nucleotide polymorphisms (SNPs, Bonferroni-adjusted -value <8.03 × 10) within the introns of the genes , , , and genes were associated with the primary feather color trait. Twenty-two genome-wide suggestive SNPs (Bonferroni-adjusted -value <1.61 × 10) of 17 genes (besides and ) were also identified. Seven SNPs were located at one 0.22 Mb region (38.65-38.87 Mb) on chromosome 5, and six SNPs were located at one 0.31 Mb region (19.53-19.84 Mb) on chromosome 11. The functions of , 5, , , and were involved pigmentation and follicle development, especially, upregulated expression in black feather (haplotype-CCCC) bulb tissue compared with in pockmarked feather (haplotype-TGTT) bulb tissue, implicating these genes as candidate genes for primary feather color trait. The preliminarily findings suggested candidate genes and regions, and the genetic basis of primary feather color trait in a female duck.