To elucidate the genes involved in the formation of white and black plumage in ducks, RNA from white and black feather bulbs of an F₂ population were analyzed using RNA-Seq. A total of 2,642 expressed sequence tags showed significant differential expression between white and black feather bulbs. Among these tags, 186 matched 133 annotated genes that grouped into 94 pathways. A number of genes controlling melanogenesis showed differential expression between the two types of feather bulbs. This differential expression was confirmed by qPCR analysis and demonstrated that Tyr (Tyrosinase) and Tyrp1 (Tyrosinase-related protein-1) were expressed not in W-W (white feather bulb from white dorsal plumage) and W-WB (white feather bulb from white-black dorsal plumage) but in B-B (black feather bulb from black dorsal plumage) and B-WB (black feather bulb from white-black dorsal plumage) feather bulbs. Tyrp2 (Tyrosinase-related protein-2) gene did not show expression in the four types of feather bulbs but expressed in retina. C-kit (The tyrosine kinase receptor) expressed in all of the samples but the relative mRNA expression in B-B or B-WB was approximately 10 fold higher than that in W-W or W-WB. Additionally, only one of the two Mitf isoforms was associated with plumage color determination. Downregulation of c-Kit and Mitf in feather bulbs may be the cause of white plumage in the duck.

Endothelin receptor B subtype 2 (EDNRB2) is a seven-transmembrane G-protein coupled receptor. In this study, we investigated EDNRB2 gene as a candidate gene for duck spot plumage pattern according to studies of chicken and Japanese quail. The entire coding region was cloned by the reverse transcription polymerase chain reaction (RT-PCR). Sequence analysis showed that duck EDNRB2 cDNA contained a 1311bp open reading frame and encoded a putative protein of 436 amino acids residues. The transcript shared 89%-90% identity with the counterparts in other avian species. A phylogenetic tree based on amino acid sequences showed that duck EDNRB2 was evolutionary conserved in avian clade. The entire coding region of EDNRB2 were sequenced in 20 spot and 20 non-spot ducks, and 13 SNPs were identified. Two of them (c.940G>A and c.995G>A) were non-synonymous substitutions, and were genotyped in 647 ducks representing non-spot and spot phenotypes. The c.995G>A mutation, which results in the amino acid substitution of Arg332His, was completely associated with the spot phenotype: all 152 spot ducks were carriers of the AA genotype and the other 495 individuals with non-spot phenotype were carriers of GA or GG genotype, respectively. Segregation in 17 GA×GG and 22 GA×GA testing combinations confirmed this association since the segregation ratios and genotypes of the offspring were in agreement with the hypothesis. In order to investigate the underlying mechanism of the spot phenotype, MITF gene was used as cell type marker of melanocyte progenitor cells while TYR and TYRP1 gene were used as cell type markers of mature melanocytes. Transcripts of MITF, TYR and TYRP1 gene with expected size were identified in all pigmented skin tissues while PCR products were not obtained from non-pigmented skin tissues. It was inferred that melanocytes are absent in non-pigmented skin tissues of spot ducks.