Endothelin Receptor B2 (EDNRB2) Gene Is Associated with Spot Plumage Pattern in Domestic Ducks (Anas platyrhynchos)

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.


Introduction
In birds, plumage color is crucial to attract opposite sex individuals for mating and to avoid predators [1,2]. Plumage coloration variants include differences in shades of basic colors (e.g. dilution), hue (e.g. black, white, yellow) and patterns (e.g. spotting, barring) [3,4]. Biosynthesis of eumelanin and pheomelanin in melanocyte is responsible for melanin-based coloration of feathers [5]. Melanocytes differentiate from undifferentiated precursors called melanoblasts, which are derived from the neural crest cells(NCC). Melanoblasts migrate from the neural crest to the epidermis and into developing feather follicles in birds [6]. To date, various genes that affect melanocyte differentiation, proliferation, migration, survival, morphology, structure and function have been shown to affect pigmentation [2,7,8].
Aberrations in development of melanocytes resulting in white spotting phenotype has been described in several species, such as dog [4], mouse [9], alpaca [10], horse [11], rabbit [12], chicken [13] and Japanese quail [14]. Up to now, endothelin 3(EDN3) gene and endothelin receptor B(EDNRB) gene have been well studied to cause pronounced white spotting phenotype [13]. Paracrine expression of EDN3 was shown to coordinate localized coat color differences between tabby marking and wild type cats [15]. Mutations in EDNRB gene were found in horse with lethal white foal syndrome [8] and piebald mice [16]. Endothelin receptor B subtype 2 (EDNRB2) gene is a paralog of the EDNRB gene [17]. EDNRB2 gene encodes a seven-transmembrane domain G-protein-coupled receptor EDNRB2 and participates in melanoblast differentiation and migration [18][19][20][21]. To date, EDNRB2 has been investigated in chicken [13,17], quail [14], fish, monotremes (platypus) [22] and frog [23], while it has been lost in therians lineages(marsupials and placentals) with the rise of the therian sex chromosomes [22]. A study by Miwa et al. [14] showed that an Arg332His amino acid change in EDNRB2 was associated with the panda plumage color mutation in Japanese quail(Coturnix coturnix). Recently, Kinoshita et al. [13] reported that a Cys244Phe mutation in EDNRB2 was associated with the mo w mutation, while an Arg332His mutation was associated with the mo mutation in chicken.
As described in livestock, many plumage color mutants of ducks are produced by the combined effect of controlled breeding and selection pressures from domestication [12,24]. These mutants provide opportunities for us to unravel mechanisms that controlling the inheritance of white spotting phenotype. In our previous research [24], we reported the interaction of alleles at two different loci determined the black(grey) plumage color in duck. We also found that plumage pattern of the ducks displayed great variation. Some white ducks had black patches on the head, back, tail and wing. In the current study, this kind of plumage pattern, which is similar to the panda Japanese quail described by Tsudzukis et al. [14,25] was termed as spot phenotype (Fig 1).
The aim of this study was to investigate causative mutation of the spot phenotype in duck. EDNRB2 gene was selected as the most likely candidate gene due to its functional importance in melanocyte development and phenotypic similarity between spot duck and panda Japanese quail [14]. We characterized the transcript and genomic structure of EDNRB2 gene in duck, performed mutation screening, investigated its association with plumage pattern, and carried out mating tests to validate the results of the association study. In addition, expression patterns of marker genes were applied to investigate the probably mechanism of the spot plumage pattern. We present strong evidence that spot plumage pattern in duck is caused by a single nucleotide substitution of the EDNRB2 gene.

Ethics statement
All animal experiments were carried out according to protocols (No. 5 proclaim of the Standing Committee of Hubei People's Congress) approved by the Standing Committee of Hubei People's Congress, and the ethics committee of Huazhong Agricultural University, China. The approved permit number for this study is ''HBAC2010113". All efforts were made to minimize the number of animals used in this study and their suffering.

Animals
A three-generation intercross between white Kaiya and white Liancheng has been generated in our previous study [24]. In order to verify the association between EDNRB2 gene and spot phenotype in duck found in current research and remove the subjective bias, we deliberately designed the mating tests according to single-blind trials [26]. White ducks were used as parental ducks since their plumage pattern were not expressed due to the absence of melanin. After genotyping of the Arg332His mutation site, four mating combinations were crossed to produce 443 grey ducks: (i) two black males with genotype GA were crossed with ten white females with genotype GG, and fourteen white females with genotype GA; (ii) two white males with genotype GA were crossed with ten black females with genotype GG, and fourteen black females with genotype GA. Collection of eggs and incubation were the same as described by Gong et al. [24]. After eliminating full-sib families with less than five numbers, segregation data of plumage pattern phenotypes were analyzed with the chi-square test by the IBM SPSS Statistics Software(Version 19.0, New York, USA). All the ducks were reared in cages in a semi-open house and subjected to conventional management conditions. High quality digital pictures were taken of all one-day-old and sixty-day-old ducks (DSC-W50, Sony, Japan and ST90, Samsung, Korean) and used for phenotypic classification.

Samples
All blood samples were collected from the F2 generation [24] comprising 647 colored ducks (non-spot and spot plumage pattern) and stored at -20°C. Tissue samples for RNA extraction were collected from skin tissues of three spot and three non-spot ducks. For each duck, six skin specimens were collected from different regions of the body (i.e., mantle, rump, breast, abdomen, the regions located at the proximal part of the wings and the distal part of the wings) after feather being plucked in each duck(S2 Fig). All skin specimens were immediately frozen in liquid nitrogen and stored at -80°C. DNA extraction, total RNA extraction and cDNA synthesis DNA samples were extracted from blood samples by the phenol-chloroform method, and store at 4°C for use [27]. Total RNA was extracted from duck skin specimens using TRIzol Reagent (Invitrogen, San Diego, CA, USA) according to the manufacturer's protocol. After RNA extraction, about 1μg DNase-treated total RNA was transcribed into cDNA using the ReverTraAce-α first strand cDNA synthesis M-MLV reverse transcriptase kit (TOYOBO, Osaka, Japan).
Sequencing of the EDNRB2 gene, phylogenetic analyses and identification of polymorphisms According to the EDNRB2 mRNA sequence of chicken and Japanese quail (GenBank accession numbers NM_204120.1, AB275309.1), we designed three pairs of primers to amplify the complete coding sequence(CDS) of duck EDNRB2 gene by the Primer Premier 5.0 software (http:// www.premierbiosoft.com). Amplification reactions were followed PCR protocols in a volume of 15 μl, using 100-200ng cDNA, 0.3μM of each primer, 10×EasyTaq Buffer, 2.0 mM MgSO 4 , 0.2 mM dNTPs, 0.5 unit EasyTaq DNA Polymerase (TransGen Biotech, Beijing, China) and ddH 2 O. Reactions were performed in a thermal cycler (Applied Biosystems, Foster City, CA, USA). PCR amplification condition was as follows: 5 min at 94°C; 35 amplification cycles of 30s at 94°C, 30s at annealing temperature, 40s at 72°C; 7 min at 72°C. PCR products were examined on 1.5% agarose gel electrophoresis and purified using the Gel Extraction Kit, then cloned into the PEASY-T1 cloning vector (TransGen Biotech, Beijing, China) and sequenced commercially (Sangon, Shanghai, China).
According to the obtained genomic sequence of duck EDNRB2 gene, seven pairs of primers were designed to identify the polymorphisms of this gene (S1 Table). 20 spot and 20 non-spot DNA samples were applied to construct two DNA poolings, respectively. PCR products were amplified from DNA pools and sequenced commercially (Sangon, Shanghai, China). All sequences were visually edited, assembled and aligned with the SeqMan procedure of DNASTAR software (http://www.dnastar.com) and Clustal W software (http://www.ebi.ac.uk/Tools/msa/ clustalw2/). Primer Premier 5.0 software was used to translate the nucleotide into protein. The protein secondary structure predictions were performed with the online website ExPASy (http://www.expasy.org/vg/index/Protein) and SMART (http://smart.embl-heidelberg.de/). Phylogeny of endothelin receptor protein sequences was obtained with the Neighbour Joining tree option of the MEGA6.0 software [28]. All sequences used for Phylogeny and multiple sequence alignment were obtained from GenBank database.

Genotyping of EDNRB2 gene and association analysis
Two pairs of primers were designed for genotyping of the 647 colored ducks of F2 generation [24] and 430 progeny of the mating tests using PCR-RFLP methods (S1 Table). PCR conditions were the same as those previously described. The restriction endonucleases NlaIII (MBI Fermentas, Hanover, MD, USA) and SfaNI (NEB, Ipswich, Massachusetts, USA) were used for the identification of the c.940G> A and c.995G> A substitution, respectively. 5 μl of PCR product was digested overnight in a standard restriction digestion protocol using 3 units of restriction enzyme at 37°C (NlaIII) and 65°C (SfaNI). Digested products were visualized on 2% (NlaIII) and 3.5% (SfaNI) agarose gel, respectively. For the c.940G>A substitution, using a PCR product amplified from genomic DNA with primer pairs NlaIII-F/R, the AA genotype is digested into two fragments of sizes 213bp and 110bp, the GA genotype is digested into three fragments of sizes 323bp, 213bp and 110bp while the GG genotype remains intact at 323bp. For the c.995G>A substitution, using a PCR product amplified from genomic DNA with primer pairs SfaNI-F/R, the GG genotype is digested into three fragments of sizes 240bp, 111bp and 23bp, the GA genotype is digested into four fragments of sizes 240bp, 134bp, 111bp and 23bp while the AA genotype is digested into two fragments of sizes 240bp and 134bp (S3 Fig). Test of significance for association between genotypes and duck plumage pattern phenotypes were conducted using Fisher 0 s exact test for 2×3 contigency tables by IBM SPSS Statistics Software (Version 19.0, New York, USA).

Expression of MITF, TYR and TYRP1 genes in skin tissues
To investigate the mechanism of the spot phenotype produced by the mutation of EDNRB2 gene, we assessed the expression levels of three marker genes that are involved in melanocyte development. The isform M of MITF gene is expressed exclusively in melanocyte and required for development of NC-derived melanocyte [29][30][31][32]. Therefore, it was selected as the marker gene for melanocyte progenitor cells. The tyrosinase gene family members (e.g. TYR, TYRP1) are pigment cell-specific genes express in differentiated melanocyte [31][32][33]. Consequently, TYR and TYRP1 gene were selected as the marker genes for melanocyte.
Four pairs of primers were designed for detecting these three marker genes and the reference gene (β-actin) by the Primer Premier 5.0 software. PCR reactions were carried out using primer pairs β-actin-F/R, MITF-F/R, TYR-F/R and TYRP1-F/R (S1 Table). cDNA samples of the 36 skin specimens mentioned previously were used as the template. Amplification conditions were as follows: 5 min at 94°C; 35 amplification cycles of 20s at 94°C, 20s at annealing temperature, 20s at 72°C; 5min at 72°C. All reactions contained 3 replicates. All PCR products were analyzed on a 2% agarose gel electrophoresis.

Duck EDNRB2 gene sequences
In our study, we amplified and sequenced 1362bp of the duck EDNRB2 gene (GenBank: KP192480). This transcript is composed of a 1311bp open reading frame (ORF) and a 51bp 5'untranslated region (UTR), which is 90% and 89% identical to the chicken and Japanese quail sequence (NM_204120.1, AB275309.1), respectively. The ORF of duck EDNRB2 gene theoretically would translate into 436 amino acids (AJL35291.1), which is 95% and 95% identical to EDNRB2 of Japanese quail and chicken (BAF37676.1, NP_989451.1), respectively. The molecular weight of duck EDNRB2 is 49.02kDa and the isoelectric point(pI) is 9.01. Duck EDNRB2 is a typical membrane protein and contains seven transmembrane domains which locate at the region of 93-115, 127-149, 164-186, 206-228, 264-286, 314-333 and 353-373, respectively. Amplification products of genomic DNA were assembled and displayed a fragment of 6079bp (GenBank: KP203838) that contained the coding sequence(CDS) and intronic sequences of the EDNRB2 gene. Exon boundaries were elucidated and numbered by comparative alignment of chicken exon sequence versus duck cDNA sequence according to the GT-AG rule.

The phylogeny of EDNRB2
In order to investigate the evolutionary relationships among various species of endothelin receptors(EDNRs), a phylogeny tree was constructed based on the deduced 436aa duck EDNRB2 protein and other EDNRs. According to the results, these endothelin receptors were divided into three subgroups, EDNRB, EDNRB2(EDNRC) and EDNRA (Fig 2). This phylogeny tree indicated that the deduced duck EDNRB2 protein is closer to avian species than to the platypus, frog and fish. Duck EDNRB2 protein showed a closer phylogenic relationship with EDNRB2 than EDNRB and EDNA receptors of the mentioned species.

Association between EDNRB2 genotypes and spot plumage pattern mutation
According to the PCR-RFLP (Fig 5) and association analyses, all 152 spot individuals demonstrated a homozygous genotype AA (240/134bp) for the c.995G>A mutation, while all 495 non-spot individuals had either homozygous genotype GG (240/111/23bp) or heterozygous genotype GA (240/134/111/23bp). This mutation showed a perfect co-segregation with the spot phenotype (P<0.0001, Fisher 0 s exact test). Moreover, c.940G>A and c. 995G>A mutations were in complete linkage phase (Table 1).

Mating tests
Later on, mating tests were carried out to confirm the co-segregation of the spot locus with the c.995G>A mutation (S1 Fig). A total of 430 individuals in 39 full-sib families were analyzed ( Table 2). In mating combinations GA×GG Ã and GA×GG # , it was assumed that all of the offspring would be non-spot plumage pattern. Data in Table 2 show that 191 progeny from 17 families were all non-spot plumage pattern. Thus, the segregation ratios of the non-spot and spot phenotype were in accordance with the expected ratio 1: 0. This ratio is based on the hypothesis that spot and non-spot phenotype are controlled by an autosomal allele and non-spot is dominant to spot phenotype. In mating combination GA×GA Ã , 80 exhibited non-spot phenotype and 25 exhibited spot phenotype. In mating combination GA×GA # , 98 exhibited nonspot phenotype and 36 exhibited spot phenotype. Both segregation ratios were in agreement with the expected 3:1 ration based on our hypothesis (P<0.80 and P<0.70, respectively). These results showed that the spot phenotype is recessive to the non-spot phenotype and AA genotype will show the spot phenotype.
All the ducks in Table 2 were genotyped to verify our hypothesis. As expected, the genotyping data of the c.995G>A and c.940G>A mutation were in accordance with our prediction. For the c.995G>A mutation, all 61 spot individuals demonstrated a homozygous genotype AA, while 171 non-spot individuals had homozygous genotype GG and 198 non-spot individuals had heterozygous genotype GA. For the c.940G>A mutation, all 61 spot individuals demonstrated a homozygous genotype GG, while 171 non-spot individuals had homozygous genotype AA and 198 non-spot individuals had heterozygous genotype GA (S3 Table).

Expression of MITF, TYR and TYRP1 genes in skin tissues
The expression patterns of the marker genes. MITF, TYR and TYRP1 genes were found to express differentially in non-spot and spot ducks (Fig 6). Non-spot ducks (n = 3) displayed a robust TYR, TYRP1 and MITF expression and these genes were expressed in all samples regardless of their location. These results demonstrated that both melanocyte progenitor cells and melanocytes were distributed evenly in non-spot ducks. Interestingly, MITF, TYR and TYRP1 gene were exclusively expressed in pigmented skin tissues which located at the rump and wings in spot ducks, while they were not detected in non-pigmented skin tissues (the region located at mantle, breast, abdomen, distal part of the wings). These results indicated that both melanocyte progenitor cells and melanocytes were only found in pigmented skin tissues, while none of them were found in non-pigmented skin tissues in spot ducks. In brief, transcripts of MITF, TYR and TYRP1 gene with expected size were identified in all pigmented skin tissues while no PCR product was obtained from non-pigmented skin tissues.

Discussion
In this study, we reported the genomic structure of the EDNRB2 gene, its full length transcripts in skin tissues and polymorphism of sequences between spot and non-spot ducks. The obtained nucleotide sequence shares about 89%-90% identical with EDNRB2 gene in the chicken and Japanese quail. The phylogeny tree showed the deduced protein of duck EDNRB2 share high amino acid identities with chicken and Japanese quail. We can infer that the biological function of EDNRB2 is conservative during evolution. The closer relationship in evolution and highly conservative structure suggest that this gene may has similar functions in avian species [34,35]. The finding of EDNRB2 gene in duck also give insight into the fact that this gene is present in aves while it was lost in the course of sex chromosome evolution in therian mammals [22].
In chicken, 4 non-synonymous mutations (Thr8Ile, Thr15Ala, Cys244Phe and Arg332His) and 5 synonymous mutations were detected [13]. In Japanese quail, 1 non-synonymous mutation (Arg332His) and 8 synonymous mutations were detected [14]. In the present study, 11 synonymous mutations and 2 non-synonymous mutations (Val314Met and Arg332His) were detected in the coding region of duck EDNRB2. We have not found the above synonymous mutations reported in chicken and Japanese quail, but identified 11 novel synonymous mutations [13,14]. Val314Met and Arg332His mutations were genotyped in 647 ducks representing nonspot and spot phenotypes. The valine to methionine change at the highly conserved position 314 is a novel mutation. However, this mutation is present only in non-spot ducks, so it is not responsible for the spot phenotype. We hypothesis that the Val314Met mutation has no or little influence on function of the protein. Another possible explanation is that phenotypic  expression of this mutation is influenced by epistasis mutations or other related genes including EphB2 [18,36]. Arg332His substitution of EDNRB2 has previously been reported to be associated with panda mutant in Japanese quail and mottled mutant in chicken, and it is located at the sixth transmembrane region [13,14]. In the present study, it is the most probable cause for the spot  Spot Plumage Caused by EDNRB2 in Duck phenotype in ducks since it was only detected in all 152 spot mutant ducks of the F2 generation [24] and 61 spot ducks of the mating tests while not detected in all non-spot ducks. The results of our mating tests provided convincing evidence that this mutation was causal mutation. To date, it is the third time that the same mutation was identified in avian species and this mutational convergence suggest that evolution is highly repeatable [37]. Therefore, this mutation is highly interesting both from a genetic and a structural perspective. It is noticed that no destructive or even lethal effects of this gene were observed in duck and the identical situation is found in chicken and quail mutants [13,14]. We can speculate that the Arg332His substitution found in EDNRB2 is a neutral mutation which is not harmful for surviving or reproduction in spite of possible strong selective pressure in the wild [1,38]. Variation in EDNRB2 sequences is able to help us to infer phylogenetic and phylogeographic relationships among organisms, assess the adaptive evolution of plumage color and understand gene evolution and evolutionary history of aves [39,40]. Although ducks diverged from the chicken approximately 90-100 million years ago [41], Arg332His substitution was found both in duck and chicken. It is postulated that this substitution of aves came from their ancestral population or introgression. This may occurred after aves' divergence from Reptilia [37,41,42]. Further investigation of the EDNRB2 in other 10,500 living birds and other divergent taxa (e.g. medaka, platpus) will afford opportunities to answer whether the Arg332His substitution is conserved throughout evolution [22,43]. Indeed, it seems reasonable to assume that the Arg residue at 332 position is important for protein function. The sixth transmembrane of G protein-coupled receptors (GPCRs) has been suggested to play a critical role in signaling and ligand selectivity [44,45]. We speculate that this mutation may change the 3-dimensional structure of the EDNRB2 protein(e.g. the flexibility of the third cytoplasmic loop) thus the ability of EDNRB2 to act as a EDN3 receptor is altered. One possible mechanism is that this mutation enhance the steric hindrances of the sixth transmembrane movement and reduce receptor activation [13,44,46,47]. Further genetic and pharmacological studies may help us to evaluate possible functional differences between the Arg332His change and to understand the mechanism of the convergent phenotypes caused by the Arg332His change [1,48].
The phenotypic outcomes of EDNRB2 mutations can be partly explained based on the current understanding of its expression pattern and function. The expression of EDNRB2 was significantly higher in Fibromelanosis skin tissue as compared to wild-type skin [49]. In panda mutant Japanese quail and mottled mutant chicken, the relative expression level of EDNRB2 mRNA from skin was lower than in the wild type bird [13,14]. Harris et al. [18] reported that in EDNRB2-siSTRIKE-treated chicken embryos, the dorsal migration of melanoblast was reduced as a result of a severe reduction in EDNRB2 expression. Therefore, it is postulated that normal expression level of EDNRB2 is indispensable to normal pigmentation. The mechanism of action of the Arg332His mutation in duck remains to be established, but one possible explanation is that this mutation also down-regulates the expression of EDNRB2 in a similar manner. Reduced expression of EDNRB2 is supposed to cause the defective in melanocyte development in spot ducks.
In order to test our hypothesis, we used RT-PCR to investigate the expression patterns of three marker genes (MITF, TYR and TYRP1) in spot and non-spot ducks. Our results suggested that unpigmented areas of the spot ducks completely lack both melanocyte progenitor cells and melanocytes. Since melanocyte progenitor cells and melanocytes are derived from the melanoblasts, we hypothesize that this situation may be caused by the early elimination of melanoblasts based on other studies [32,50]. This is supported by the recent studies on the EDN3/ EDNRB2 signaling which is required for melanoblast migration in EDNRB2 gene-conserved animals. EDNRB2 is considered to be a powerful tool for studying development of melanophore development in Xenopus and fate restrictions of premigratory avian neural crest [20,23].
In Xenopus, EDNRB2 has an important role in melanoblast migration [23]. In chickens and quails, EDNRB2 is upregulated in neural crest cells before the initiation of dorsolateral migration. Expression of EDNRB2 in chick help the melanoblasts to overcome the high levels of repulsive cues in the dorsolateral environment, and it is required for melanoblasts to enter the migratory pathway and to differentiate [17][18][19]. Taking the evidence together, it is plausible that EDNRB2 participates in melanocyte development in duck. In the future study, it is necessary to investigate the underlying mechanism of the lack of melanocyte in non-pigmented areas, to identify the determinants including EDNRB2 in dorsolateral pathfinding of melanoblasts in duck and compare the results with that obtained in chicken and quail. This can be investigated by analyzing expression patterns of EDN3, EDNRB2, MITF, TYR, TYRP1 as well as some other related genes that may be affected during the development stage of duck embryos.
Researches in the candidate genes of pigmentation are able to provide insights into human diseases, the functions and network of pigmentation related genes, and other coat color mutants with similar molecular alterations [6,13]. Given the recent advances in surgical manipulation, cell culture, cell marking techniques, transgenesis by electroporation, retrovirally mediate gene transfer and chimeras technique which have been successfully applied in study of chicken and Japanese quail embryo, rapid progress in research of duck embryo is likely [51][52][53][54][55]. All of the mentioned above suggest that it is indeed possible to use ducks with pigmentation mutations as a model to learn the development of melanocyte. Any mechanistic insights gained in duck could also be used to understand the role of genes involved in pigmentation.

Conclusion
In this study, we confirmed the significant role of EDNRB2 mutations in spot ducks and identified the causative mutation. The c.995G>A mutation results in the amino acid change Arg332His, was completely associated with the spot phenotype in duck. Our results of mating tests confirmed this association. It was inferred that melanocytes are absent in non-pigmented skin tissues of spot ducks by detected expression patterns of marker genes involved in melanocyte development. As the first study to detect polymorphism in EDNRB2 gene of duck, the present demonstration of the spot mutation may help us to clarify the genetics mechanism of coloration in ducks. Future study may explore the functional difference between the Arg332His change and evolution of EDNRB2 gene.
Supporting Information S1 Checklist. Completed "The ARRIVE Guidelines Checklist" for reporting animal data in this manuscript.