Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.

Avian egg coloration is shaped by natural selection, but its genetic basis remains unclear. Here, we used genome‐wide association analysis and identity by descent to finely map green egg colour to a 179‐kb region of Chr4 based on the resequencing of 352 ducks (Anas platyrhynchos) from a segregating population resulting from the mating of Pekin ducks (white‐shelled eggs) and mallards (green‐shelled eggs). We further narrowed the candidate region to a 30‐kb interval by comparing genome divergence in seven indigenous duck populations. Among the genes located in the finely mapped region, only one transcript of the ABCG2 gene (XM_013093252.2) exhibited higher uterine expression in green‐shelled individuals than in white‐shelled individuals, as supported by transcriptome data from four populations. ABCG2 has been reported to encode a protein that functions as a membrane transporter for biliverdin. Sanger sequencing of the whole 30‐kb candidate region (Chr4: 47.41–47.44 Mb) and a plasmid reporter assay helped to identify a single nucleotide polymorphism (Chr4: 47,418,074 G>A) located in a conserved predicted promoter region whose variation may alter ABCG2 transcription activity. We provide a useful molecular marker for duck breeding and contribute data to the research on ecological evolution based on egg colour patterns among birds.

The egg-laying rate is an important indicator for evaluating fertility of poultry. In order to better understand the laying mechanism of Muscovy ducks, gene expression profiles and pathways of ovarian tissues in high- and low-laying black (BH and BL) and white Muscovy ducks (WH and WL) during the peak production period were performed by using RNA-seq. The total number of reads produced for each ovarian sample ranged from 44,344,070 to 47,963,328. A total of 113, 619 and 87 differentially expressed genes (DEGs) were identified in BH-vs-WH, BL-vs-BH and BL-vs-WL, respectively. Among them, 54, 356 and 49 genes were up regulated and 59, 263 and 38 genes were down regulated. In addition, there were only 10 up-regulated genes in WL-vs-WH. In the comparison of DEGs in black and white Muscovy ducks, two co-expressed DEG genes were detected between BH-vs-WH and BL-vs-WL and seven DEGs were co-expressed between BL-vs-BH and WL-vs-WH. The RNA-Seq data were confirmed to be reliable by qPCR. Numerous DEGs known to be involved in ovarian development were identified, including TGFβ2, NGFR, CEBPD, CPEB2, POSTN, SMOC1, FGF18, EFNA5 and SDC4. Gene Ontology (GO) annotations indicated that DEGs related to ovarian development were mainly enriched in biological processes of “circadian sleep/wake cycle process,” “negative regulation of transforming growth factor-β secretion,” “positive regulation of calcium ion transport” in BH-vs-WH and “cell surface receptor signaling pathway,” “Notch signaling pathway” and “calcium ion transport” in BL-vs-BH. Besides, “steroid biosynthetic process,” “granulosa cell development” and “egg coat formation” were mainly enriched in BL-vs-WL and “reproduction,” “MAPK cascade” and “mitotic cell cycle” were mainly enriched in WL-vs-WH. KEGG pathway analysis showed that the PI3K-Akt signaling pathway and ovarian steroidogenesis were the most enriched in Muscovy duck ovary transcriptome data. This work highlights potential genes and pathways that may affect ovarian development in Muscovy duck.