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

Transcriptome analysis of sexual dimorphism in dorsal down coloration in goslings

BACKGROUND

In day-old Hungarian white goose goslings, there is a noticeable difference in dorsal down coloration between males and females, with females having darker dorsal plumage and males having lighter plumage. The ability to autosex day-old goslings based on their dorsal down coloration is important for managing them efficiently and planning their nutrition in the poultry industry. The aim of this study was to determine the biological and genetic factors underlying this difference in dorsal down colorationthrough histological analysis, biochemical assays, transcriptomic profiling, and q‒PCR analysis.

RESULTS

Tissue analysis and biochemical assays revealed that compared with males, 17-day-old embryos and day-old goslings of female geese exhibited a greater density of melanin-containing feather follicles and a greater melanin concentration in these follicles during development. Both female and male goslings had lower melanin concentrations in their dorsal skin compared to 17-day-old embryos. Transcriptome analysis identified a set of differentially expressed genes (DEGs) (MC1R, TYR, TYRP1, DCT and MITF) associated with melanogenesis pathways that were downregulated or silenced specifically in the dorsal skin of day-old goslings compared to 17-day-old embryos, affecting melanin synthesis in feather follicles. Additionally, two key genes (MC1R and MITF) associated with feather coloration showed differences between males and females, with females having higher expression levels correlated with increased melanin synthesis and darker plumage.

CONCLUSION

The expression of multiple melanogenesis genes determines melanin synthesis in goose feather follicles. The dorsal down coloration of day-old Hungarian white goose goslings shows sexual dimorphism, likely due to differences in the expression of the MC1R and MITF genes between males and females. These results could help us better understand why male and female goslings exhibit different plumage patterns.

Identification of Differentially Expressed Genes and microRNAs in the Gray and White Feather Follicles of Shitou Geese

The Shitou goose, a highly recognized indigenous breed with gray plumage originating from Chaozhou Raoping in Guangdong Province, China, is renowned for being the largest goose species in the country. Notably, during the pure breeding process of Shitou geese, approximately 2% of the offspring in each generation unexpectedly exhibited white plumage. To better understand the mechanisms underlying white plumage color formation in Shitou geese, we conducted a comparative transcriptome analysis between white and gray feather follicles, aiming to identify key genes and microRNAs that potentially regulate white plumage coloration in this unique goose breed. Our results revealed a number of pigmentation genes, encompassing TYR, TYRP1, EDNRB2, MLANA, SOX10, SLC45A2, GPR143, TRPM1, OCA2, ASIP, KIT, and SLC24A5, which were significantly down-regulated in the white feather follicles of Shitou geese. Among these genes, EDNRB2 and KIT emerged as the most promising candidate genes for white plumage coloration in Shitou geese. Additionally, our analysis also uncovered 46 differentially expressed miRNAs. Of these, miR-144-y may play crucial roles in the regulation of feather pigmentation. Furthermore, the expression of novel-m0086-5p, miR-489-y, miR-223-x, miR-7565-z, and miR-3535-z exhibits a significant negative correlation with the expression of pigmentation genes including TYRP1, EDNRB2, MLANA, SOX10, TRPM1, and KIT, suggesting these miRNAs may indirectly regulate the expression of these genes, thereby influencing feather color. Our findings provide valuable insights into the genetic mechanisms underlying white plumage coloration in Shitou geese and contribute to the broader understanding of avian genetics and coloration research.

Genetic insights into fiber quality, coat color and adaptation in Changthangi and Muzzafarnagri sheep: A comparative skin transcriptome analysis

Changthangi sheep, which inhabit the high-altitude regions of Ladakh, are known for their fine fiber production and are characterized by grey skin and either black or white coats. In contrast, Muzzafarnagri sheep from the plains of Uttar Pradesh produce coarse wool and have white skin and coats. We conducted comparative global gene expression profiling on four biological replicates of skin from each breed. Notably, our analysis identified 149 up-regulated genes and 2,139 down-regulated genes in Changthangi sheep compared to Muzzafarnagri sheep, with a p-adjusted value (padj) of ≤0.05 and a Log2 fold change of ≥1.5. Gene Ontology analysis of the up-regulated genes revealed an enrichment of terms related to melanin biosynthesis and developmental pigmentation. Additionally, enriched KEGG pathways included tyrosine metabolism and metabolic pathways. Among the melanogenesis-related genes that exhibited higher expression in Changthangi sheep were TYR, TYRP1, DCT, SLC45A2, PMEL, MLANA, and OCA2. These findings confirm melanin's role in both the animals' black coat color and UV protection at high-altitude. Furthermore, we observed more pronounced expression of genes related to fiber quality, namely KRTAP6, KRTAP7, KRTAP13, and KRTAP2, in the fine wool-producing sheep from Ladakh. The results of the RNA sequencing were validated using real-time PCR on 10 genes governing fiber quality and coat color, with ACTB and PPIB serving as reference genes. In conclusion, our comparative skin transcriptome analysis of Changthangi and Muzzafarnagri sheep sheds light on the genetic differences associated with distinct phenotypic traits and environmental adaptability, offering valuable insights into the underlying mechanisms.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

Stripes and loss of color in ball pythons (Python regius) are associated with variants affecting endothelin signaling

Color patterns in nonavian reptiles are beautifully diverse, but little is known about the genetics and development of these patterns. Here, we investigated color patterning in pet ball pythons (Python regius), which have been bred to show color phenotypes that differ dramatically from the wildtype form. We report that several color phenotypes in pet animals are associated with putative loss-of-function variants in the gene encoding endothelin receptor EDNRB1: (1) frameshift variants in EDNRB1 are associated with conversion of the normal mottled color pattern to skin that is almost fully white, (2) missense variants affecting conserved sites of the EDNRB1 protein are associated with dorsal, longitudinal stripes, and (3) substitutions at EDNRB1 splice donors are associated with subtle changes in patterning compared to wildtype. We propose that these phenotypes are caused by loss of specialized color cells (chromatophores), with loss ranging from severe (fully white) to moderate (dorsal striping) to mild (subtle changes in patterning). Our study is the first to describe variants affecting endothelin signaling in a nonavian reptile and suggests that reductions in endothelin signaling in ball pythons can produce a variety of color phenotypes, depending on the degree of color cell loss.

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

Background: 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. Methods: 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. Results: Seven genome-wide significant single nucleotide polymorphisms (SNPs, Bonferroni-adjusted p-value <8.03 × 10^-7) within the introns of the genes STARD9, ZNF106, SLC7A5, and BANP genes were associated with the primary feather color trait. Twenty-two genome-wide suggestive SNPs (Bonferroni-adjusted p-value <1.61 × 10^-5) of 17 genes (besides ZNF106 and SLC7A5) 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 STARD9, SLC7A5, BANP, LOC101798015, and IPMK were involved pigmentation and follicle development, especially, STARD9 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. Conclusion: The preliminarily findings suggested candidate genes and regions, and the genetic basis of primary feather color trait in a female duck.

Genome-Wide Population Structure Analysis and Genetic Diversity Detection of Four Chinese Indigenous Duck Breeds from Fujian Province

Simple Summary

The aim of this study was to conduct a genome-wide comparative analysis of four indigenous Chinese duck breeds (Jinding, Liancheng white, Putian black, and Shanma ducks) from Fujian Province, to understand their genetic diversity and population structure. Population parameters showed that the four indigenous breeds were separated groups. Five genomic regions are presented as hotspots of autozygosity among these indigenous duck breeds, with candidate genes involved in muscle growth, pigmentation, and neuroregulation. Genomic information may play a vital role in the improvement of conservation strategies.

Abstract

The assessment of population genetic structure is the basis for understanding the genetic information of indigenous breeds and is important for the protection and management of indigenous breeds. However, the population genetic differentiation of many local breeds still remains unclear. Here, we performed a genome-wide comparative analysis of Jinding, Liancheng white, Putian black, and Shanma ducks based on the genomic sequences using RAD sequencing to understand their population structure and genetic diversity. The population parameters showed that there were obvious genetic differences among the four indigenous breeds, which were separated groups. Among them, Liancheng white and Shanma ducks may come from the same ancestor because the phylogenetic tree forms three tree trunks. In addition, during the runs of homozygosity (ROH), we found that the average inbreeding coefficient of Liancheng white and Putian black ducks was the lowest and the highest, respectively. Five genomic regions were considered to be the hotspots of autozygosity among these indigenous duck breeds, and the candidate genes involved a variety of potential variations, such as muscle growth, pigmentation, and neuroregulation. These findings provide insights into the further improvement and conservation of Fujian duck breeds.

Using comparative genomics to detect mutations regulating plumage variations in graylag (A. anser) and swan geese (A. cygnoides)

Although graylag geese (A. anser) showed similar plumages of white, grey, and white with grey patches compared to those in swan geese (A. cygnoides), it was believed the substantial molecular mechanism for plumage variations were different. To date, studies on genes responsible for diverse plumages among graylag geese were limited and causal mutations remain unknown. In this study, genomes from 57 individuals belonging to six breeds showing different plumages were sequenced at ∼10X depth. Firstly, the allele frequency differences (AFD) of variants on the scaffold394 (NW_013185915.1) between grey and white goose breeds (A. anser) was calculated and a genomic region between 768,290-779,889 bp was detected to carry candidate variants associated with plumages, including one SNP (g. 775,151G > T, ∼18.6 kb upstream of EDNRB2) found to be fixed in white geese. This region was overlapped with the one detected by the haplotype-based sweep analysis, in which significant signals defined a candidate region of 736,610-820,622 bp on the same scaffold. Results from the transcriptomic data showed that expression levels of EDNRB2 and many other melanogenesis-related genes were significantly decreased among white geese compared to that in grey geese, especially at late embryonic stages (>E15). Modifications at transcriptional levels might result in abnormal melanocyte developments and thus the white plumages when they grow up. In addition, a frameshift mutation (C > -) in exon4 of MLANA gene on scaffold176 (NW_013185876.1) was suggested as the causal mutation for sex-linked dilution phenotype in graylag geese although this requires more demonstration experiments. Together with observed white plumages caused by EDNRB2 mutations in coding regions among swan geese and chicken, our study provided new examples to study the parallel evolution.

Genome-Wide Association Study for Screening and Identifying Potential Shin Color Loci in Ducks

Shin color diversity is a widespread phenomenon in birds. In this study, ducks were assessed to identify candidate genes for yellow, black, and spotted tibiae. For this purpose, we performed whole-genome resequencing of an F2 population consisting of 275 ducks crossed between Runzhou crested-white ducks and Cherry Valley ducks. We obtained 12.6 Mb of single nucleotide polymorphism (SNP) data, and the three shin colors were subsequently genotyped. Genome-wide association studies (GWASs) were performed to identify candidate and potential SNPs for the three shin colors. According to the results, 2947 and 3451 significant SNPs were associated with black and yellow shins, respectively, and six potential SNPs were associated with spotted shins. Based on the SNP annotations, the MITF, EDNRB2, POU family members, and the SLC superfamily were the candidate genes regulating pigmentation. In addition, the isoforms of EDNRB2, TYR, TYRP1, and MITF-M were significantly different between the black and yellow tibiae. MITF and EDNRB2 may have synergistic roles in the regulation of melanin synthesis, and their mutations may lead to phenotypic differences in the melanin deposition between individuals. This study provides new insights into the genetic factors that may influence tibia color diversity in birds.

Transcriptome analysis of genes potentially associated with white and black plumage formation in Chinese indigenous ducks (Anas platyrhynchos)

1. Plumage colour is an important recognisable characteristic of duck (Anas platyrhynchos), but the coloration mechanisms remain largely unknown. To elucidate the molecular mechanisms underlying the formation of black and white plumage, the following study applied RNA sequencing (RNA-Seq) to catalogue the global gene expression profiles in the duck feather bulbs of black and white colours.2. Black feather bulbs were collected from Putian Black ducks (B-PTB) and black Longsheng Jade-green ducks (B-LS), while white feather bulbs were collected from Putian White ducks (W-PTW), Putian Black ducks (W-PTB) and Longsheng Jade-green ducks (W-LS). Sixteen cDNA libraries were constructed and sequenced for transcriptional analysis. Three comparison groups were employed to analyse differentially expressed genes (DEGs), including W-PTB versus B-PTB, W-PTW versus B-PTB and W-LS versus B-LS.3. The results showed 180 DEGs between W-PTB and B-PTB, 303 DEGs between W-PTW and B-PTB, and 108 DEGs between W-LS and B-LS. Further analysis showed that 18 DEGs were directly involved in the pigmentation process and melanogenesis signalling pathway. Additionally, the distribution of DEGs varied amongst groups whereby ASIP appeared only in the W-LS versus B-LS group, GNAI1 and ZEB2 appeared only in the W-PTW versus B-PTB group, and KITLG, EDN3 and FZD4 appeared only in W-PTB versus B-PTB.4. The findings suggested that the mechanism of feather albinism may differ between duck breeds. This study provided new information for discovering genes that are important for feather pigmentation and helps elucidate molecular mechanisms involved in black and white plumage in ducks.

Two Genomic Loci Control Three Eye Colors in the Domestic Pigeon (Columba livia)

The iris of the eye shows striking color variation across vertebrate species, and may play important roles in crypsis and communication. The domestic pigeon (Columba livia) has three common iris colors, orange, pearl (white), and bull (dark brown), segregating in a single species, thereby providing a unique opportunity to identify the genetic basis of iris coloration. We used comparative genomics and genetic mapping in laboratory crosses to identify two candidate genes that control variation in iris color in domestic pigeons. We identified a nonsense mutation in the solute carrier SLC2A11B that is shared among all pigeons with pearl eye color, and a locus associated with bull eye color that includes EDNRB2, a gene involved in neural crest migration and pigment development. However, bull eye is likely controlled by a heterogeneous collection of alleles across pigeon breeds. We also found that the EDNRB2 region is associated with regionalized plumage depigmentation (piebalding). Our study identifies two candidate genes for eye colors variation, and establishes a genetic link between iris and plumage color, two traits that vary widely in the evolution of birds and other vertebrates.

New insights into genetics underlying of plumage color

Plumage color can be considered as a social signal in chickens and a breeding identification tool among breeders. The relationship between plumage color and trait groups of immunity, growth and fertility is still a controversial issue. This research aimed to determine the genome-wide additive and epistatic variants affecting plumage color variation in chickens using the chicken Illumina 60k high-density SNP array. Two scenarios of genome-wide additive association studies using all SNPs and independent SNPs were carried out. To perform epistatic association analysis, the LD pruning approach was used to reduce the complexity of the analysis. We detected seven novel significant loci using all of the SNPs in the model and 14 SNPs using the LD pruning approach associated with plumage color. Moreover, 89 significantly associated SNP-SNP interactions (P-value <10^-6 ) distributed in 25 chromosomes were identified, indicating that all of the signals together putatively influence the quantitative variation of plumage color. By annotating genes relevant to top SNPs, we have distinguished 18 potential candidate genes comprising HNF4beta, CKMT1B, TBC1D22A, RPL8, CACNA2D1, FZD4, SGMS1, IRF8, OPTN, LOC420362, TRABD, OvoDA1, DAD1, USP6, RBM12B, MIR1772, MIR1709 and MIR6696 and also 89 putative gene-gene combinations responsible for plumage color variation in chickens. Furthermore, several KEGG pathways including metabolic pathway, cytokine-cytokine receptor interaction, focal adhesion, melanogenesis, glycosaminoglycan biosynthesis-keratan sulfate and sphingolipid metabolism were enriched in the gene-set analysis. The results indicated that plumage color is a highly polygenic trait which, in turn, can be affected by multiple coding genes, regulatory genes and gene-gene epistasis interactions. In addition to genes with additive effects, epistatic genes with tiny individual effect sizes but significant effects in a pair have the potential to control plumage coloration in chickens.

Genome-wide association analysis reveals that EDNRB2 causes a dose-dependent loss of pigmentation in ducks

Background

Birds have various plumage color patterns, and spot is a common phenotype. Herein, we conducted genome-wide association studies (GWAS) in a population of 225 ducks with different sized black spots to reveal the genetic basis of this phenomenon.

Results

First, we quantified the black spot phenotype within the duck population. The results showed that the uncolored area of the body surface first appeared on the ventral side. With increasing duck age, the area of the black spots was highly conserved across the whole body surface. The GWAS results identified a 198 kb (Chr4: 10,149,651 bp to 10,348,068 bp) genetic region that was significantly associated with the black spot phenotype. The conditional GWAS and linkage disequilibrium (LD) analysis further narrowed the ultimate candidate region to 167 kb (Chr4: 10,180,939 bp to 10,348,068 bp). A key gene regulating melanoblast migration and differentiation, EDNRB2 (Endothelin B receptor-like), was found in the candidate region and having significant mRNA expression level changes in embryonic duck skin tissue with different spot sizes. The significant SNPs (single nucleotide polymorphisms) associated with the EDNRB2 gene were annotated, and two mutations (Chr4: 10,180,939 T > C and Chr4: 10,190,671 A > T) were found to result in the loss of binding sites for two trans-factors, XBP1 and cMYB. The phenotypic effect of these two mutations suggested that they can regulate the size of black spots in a dose-dependent manner, and Chr4: 10,180,939 T > C was the major allele locus.

Conclusions

Our results revealed that EDNRB2 was the gene responsible for the variation in duck body surface spot size. Chr4: 10,180,939 T > C was the major allele that explained 49.5 % (dorsal side) and 32.9 % (ventral side) of the variation in duck body surface spot size, while 32.1 % (dorsal side) and 19.1 % (ventral side) of the variation could be explained by Chr4: 10,190,671 A > T. The trans-factor prediction also suggested that XBP1 and cMYB have the potential to interact with EDNRB2, providing new insights into the mechanism of action of these genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07719-7.

A complex genetic architecture in zebrafish relatives Danio quagga and D. kyathit underlies development of stripes and spots

Vertebrate pigmentation is a fundamentally important, multifaceted phenotype. Zebrafish, Danio rerio, has been a valuable model for understanding genetics and development of pigment pattern formation due to its genetic and experimental tractability, advantages that are shared across several Danio species having a striking array of pigment patterns. Here, we use the sister species D. quagga and D. kyathit, with stripes and spots, respectively, to understand how natural genetic variation impacts phenotypes at cellular and organismal levels. We first show that D. quagga and D. kyathit phenotypes resemble those of wild-type D. rerio and several single locus mutants of D. rerio, respectively, in a morphospace defined by pattern variation along dorsoventral and anteroposterior axes. We then identify differences in patterning at the cellular level between D. quagga and D. kyathit by repeated daily imaging during pattern development and quantitative comparisons of adult phenotypes, revealing that patterns are similar initially but diverge ontogenetically. To assess the genetic architecture of these differences, we employ reduced-representation sequencing of second-generation hybrids. Despite the similarity of D. quagga to D. rerio, and D. kyathit to some D. rerio mutants, our analyses reveal a complex genetic basis for differences between D. quagga and D. kyathit, with several quantitative trait loci contributing to variation in overall pattern and cellular phenotypes, epistatic interactions between loci, and abundant segregating variation within species. Our findings provide a window into the evolutionary genetics of pattern-forming mechanisms in Danio and highlight the complexity of differences that can arise even between sister species. Further studies of natural genetic diversity underlying pattern variation in D. quagga and D. kyathit should provide insights complementary to those from zebrafish mutant phenotypes and more distant species comparisons.

Pigment patterns of fishes are diverse and function in a wide range of behaviors. Common pattern themes include stripes and spots, exemplified by the closely related minnows Danio quagga and D. kyathit, respectively. We show that these patterns arise late in development owing to alterations in the development and arrangements of pigment cells. In the closely related model organism zebrafish (D. rerio) single genes can switch the pattern from stripes to spots. Yet, we show that pattern differences between D. quagga and D. kyathit have a more complex genetic basis, depending on multiple genes and interactions between these genes. Our findings illustrate the importance of characterizing naturally occurring genetic variants, in addition to laboratory induced mutations, for a more complete understanding of pigment pattern development and evolution.

Transcriptome Reveals Multi Pigmentation Genes Affecting Dorsoventral Pattern in Avian Body

Certain animals exhibit a special dorsoventral pattern with a lighter ventral side compared to the dorsal one and this phenomenon was preserved in the long-term evolution process. Birds also retain this trait. Recently, Inaba et al. (2019) found that ASIP (agouti signal protein) regulated interconversion between different melanocyte types leads to dorsal stripe pattern, which may partly explain the birds' dorsoventral plumage color difference. In this study, we used the embryo samples of LBM (light brown mottling) ducks (Anas platyrhynchos) with white ventral and dark dorsal body parts to investigate the mechanism of dorsoventral color variation. Firstly, melanin deposition process of duck embryos was investigated. The result indicated that E13 and E16 were the active stages of melanin synthesis. Moreover, the melanin deposition on the dorsum of LBM ducks was higher than that on the ventral side throughout. Then, RNA-seq was conducted for the dorsal and ventral skin tissues from E7 (early), E13 (middle) and E19 (late) of LBM ducks. Expression pattern analysis showed that the mRNA expression of most melanin synthesis related genes were at the highest level at E13, which was consistent with the section analysis. A correlation was found between melanogenesis pathway and dorsoventral color difference by co-expression analysis. In the DEG (differentially expressed gene) analysis, we added the dorsal skin transcriptome of embryonic white and black duck of same subspecies (Anas platyrhynchos domestica) for horizontal comparison. The results showed that 8 melanogenesis related genes (TYR, TYRP1, MLANA, RAB38, OCA2, TSPAN10, MC1R, and MSLN) were the common DEGs (Differential expressed genes) in the comparisons of body parts and breeds suggesting that the underlying molecular regulatory mechanism of dorsoventral plumage color difference may be similar to that of albino and melanic duck, which were caused by the different expression of multiple genes in melanin synthesis pathway. In addition, the molecular regulation of melanin synthesis pathway in the dorsal and ventral side of LBM ducks was analyzed. In this pathway, ASIP, MC1R, TYR, and TYRP1 have differential mRNA expression. ASIP, as an upstream gene in this pathway, was likely to play a decisive role in determining the dorsoventral plumage pattern.

Combined transcriptomics and proteomics forecast analysis for potential genes regulating the Columbian plumage color in chickens

BACKGROUND

Coloration is one of the most recognizable characteristics in chickens, and clarifying the coloration mechanisms will help us understand feather color formation. "Yufen I" is a commercial egg-laying chicken breed in China that was developed by a three-line cross using lines H, N and D. Columbian plumage is a typical feather character of the "Yufen I" H line. To elucidate the molecular mechanism underlying the pigmentation of Columbian plumage, this study utilizes high-throughput sequencing technology to compare the transcriptome and proteome differences in the follicular tissue of different feathers, including the dorsal neck with black and white striped feather follicles (Group A) and the ventral neck with white feather follicles (Group B) in the "Yufen I" H line.

RESULTS

In this study, we identified a total of 21,306 genes and 5,203 proteins in chicken feather follicles. Among these, 209 genes and 382 proteins were differentially expressed in two locations, Group A and Group B, respectively. A total of 8 differentially expressed genes (DEGs) and 9 differentially expressed proteins (DEPs) were found to be involved in the melanogenesis pathway. Additionally, a specifically expressed MED23 gene and a differentially expressed GNAQ protein were involved in melanin synthesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis mapped 190 DEGs and 322 DEPs to 175 and 242 pathways, respectively, and there were 166 pathways correlated with both DEGs and DEPs. 49 DEPs/DEGs overlapped and were enriched for 12 pathways. Transcriptomic and proteomic analyses revealed that the following pathways were activated: melanogenesis, cardiomyocyte adrenergic, calcium and cGMP-PKG. The expression of DEGs was validated by real-time quantitative polymerase chain reaction (qRT-PCR) that produced results similar to those from RNA-seq. In addition, we found that the expression of the MED23, FZD10, WNT7B and WNT11 genes peaked at approximately 8 weeks in the "Yufen I" H line, which is consistent with the molting cycle. As both groups showed significant differences in terms of the expression of the studied genes, this work opens up avenues for research in the future to assess their exact function in determining plumage color.

CONCLUSION

Common DEGs and DEPs were enriched in the melanogenesis pathway. MED23 and GNAQ were also reported to play a crucial role in melanin synthesis. In addition, this study is the first to reveal gene and protein variations in in the "Yufen I" H line during Columbian feather color development and to discover principal genes and proteins that will aid in functional genomics studies in the future. The results of the present study provide a significant conceptual basis for the future breeding schemes with the "Yufen I" H line and provide a basis for research on the mechanisms of feather pigmentation.

Transcriptome Analysis of the Breast Muscle of Xichuan Black-Bone Chickens Under Tyrosine Supplementation Revealed the Mechanism of Tyrosine-Induced Melanin Deposition

The Xichuan black-bone chicken, which is a rare local chicken species in China, is an important genetic resource of black-bone chickens. Tyrosine can affect melanin production, but the molecular mechanism underlying tyrosine-induced melanin deposition in Xichuan black-bone chickens is poorly understood. Here, the blackness degree and melanin content of the breast muscle of Xichuan black-bone chickens fed a basic diet with five levels of added tyrosine (i.e., 0.2, 0.4, 0.6, 0.8, and 1.0%; these groups were denoted test groups I-V, respectively) were assessed, and the results showed that 0.8% tyrosine was the optimal level of added tyrosine. Moreover, the effects of tyrosine supplementation on the proliferation and tyrosinase content of melanocytes in Xichuan black-bone chickens were evaluated. The results revealed a dose-dependent relationship between tyrosine supplementation and melanocyte proliferation. In addition, 417 differentially expressed genes (DEGs), including 160 upregulated genes and 257 downregulated genes, were identified in a comparative analysis of the transcriptome profiles constructed using the pooled total RNA from breast muscle tissues of the control group and test group IV, respectively (fold change ≥2.0, P < 0.05). These DEGs were mainly involved in melanogenesis, the calcium signaling pathway, the Wnt signaling pathway, the mTOR signaling pathway, and vascular smooth muscle contraction. The pathway analysis of the DEGs identified some key genes associated with pigmentation, such as DCT and EDNRB2. In summary, the melanin content of breast muscle could be markedly enhanced by adding an appropriate amount of tyrosine to the diet of Xichuan black-bone chickens, and the EDNRB2-mediated molecular regulatory network could play a key role in the biological process of tyrosine-induced melanin deposition. These results have deepened the understanding of the molecular regulatory mechanism of melanin deposition in black-bone chickens and provide a basis for the regulation of nutrition and genetic breeding associated with melanin deposition in Xichuan black-bone chickens.

Comparative transcriptome analysis reveals the genetic basis of coat color variation in Pashmina goat

The genetics of coat color variation remains a classic area. Earlier studies have focused on a limited number of genes involved in color determination; however, the complete set of trait determinants are still not well known. In this study, we used high-throughput sequencing technology to identify and characterize intricate interactions between genes that cause complex coat color variation in Changthangi Pashmina goats, producer of finest and costly commercial animal fiber. We systematically identified differentially expressed mRNAs and lncRNAs from black, brown and white Pashmina goat skin samples by using RNA-sequencing technique. A pairwise comparison of black, white and brown skin samples yielded 2479 significantly dysregulated genes (2422 mRNA and 57 lncRNAs). Differentially expressed genes were enriched in melanin biosynthesis, melanocyte differentiation, developmental pigmentation, melanosome transport activities GO terms. Our analysis suggested the potential role of lncRNAs on color coding mRNAs in cis and trans configuration. We have also developed online data repository as a component of the study to provide a central location for data access, visualization and interpretation accessible through http://pcd.skuastk.org/.

The recessive C locus in the MITF gene plays a key regulatory role in the plumage colour pattern of duck (Anas platyrhynchos)

1. The c/c alleles are responsible for the white plumage colour of ducks; however, the gene corresponding to this locus is still unclear. In order to identify the locus-related candidate gene associated with duck's plumage colour pattern, it was necessary to analyse the whole genome resequencing data. 2. A total of 929,465 SNPs in chromosome 13 and 1,688 SNPs in the region of the Microphthalmia-Associated Transcription Factor (MITF) gene were identified from whole genome resequencing data analysis. After construction of an F_ST plot from chromosome 13, MITF was highlighted as a candidate gene, possessing the highest F_ST value (0.811) on chromosome 13. 3. Six novel SNPs were discovered, located in the intronic region of the MITF gene. F2 progeny of Kaiya × Liancheng ducks (N = 1,061) were selected for genotyping by the Restriction Fragment Length Polymorphism (RFLP) technique. Association analysis using Haploview software was used for validation of the results. 4. Association results between SNPs and phenotypes showed significant association with corresponding phenotypes. All the significantly associated SNPs were located in the identified candidate gene. 5. The identified candidate gene provided novel information which is important in marker-assisted selection and breeding of duck and for the investigation of the C locus recessive white genetic mechanisms underlying plumage colour pattern.

Evolution of Endothelin signaling and diversification of adult pigment pattern in Danio fishes

Fishes of the genus Danio exhibit diverse pigment patterns that serve as useful models for understanding the genes and cell behaviors underlying the evolution of adult form. Among these species, zebrafish D. rerio exhibit several dark stripes of melanophores with sparse iridophores that alternate with light interstripes of dense iridophores and xanthophores. By contrast, the closely related species D. nigrofasciatus has an attenuated pattern with fewer melanophores, stripes and interstripes. Here we demonstrate species differences in iridophore development that presage the fully formed patterns. Using genetic and transgenic approaches we identify the secreted peptide Endothelin-3 (Edn3)—a known melanogenic factor of tetrapods—as contributing to reduced iridophore proliferation and fewer stripes and interstripes in D. nigrofasciatus. We further show the locus encoding this factor is expressed at lower levels in D. nigrofasciatus owing to cis-regulatory differences between species. Finally, we show that functions of two paralogous loci encoding Edn3 have been partitioned between skin and non-skin iridophores. Our findings reveal genetic and cellular mechanisms contributing to pattern differences between these species and suggest a model for evolutionary changes in Edn3 requirements for pigment patterning and its diversification across vertebrates.

Neural crest derived pigment cells generate the spectacular variation in skin pigment patterns among vertebrates. Mammals and birds have just a single skin pigment cell, the melanocyte, whereas ectothermic vertebrates have several pigment cells including melanophores, iridophores and xanthophores, that together organize into a diverse array of patterns. In the teleost zebrafish, Danio rerio, an adult pattern of stripes depends on interactions between pigment cell classes and between pigment cells and their tissue environment. The close relative D. nigrofasciatus has fewer stripes and prior analyses suggested a difference between these species that lies extrinsic to the pigment cells themselves. A candidate for mediating this difference is Endothelin-3 (Edn3), essential for melanocyte development in warm-blooded animals, and required by all three classes of pigment cells in an amphibian. We show that Edn3 specifically promotes iridophore development in Danio, and that differences in Edn3 expression contribute to differences in iridophore complements, and striping, between D. rerio and D. nigrofasciatus. Our study reveals a novel function for Edn3 and provides new insights into how changes in gene expression yield morphogenetic outcomes to effect diversification of adult form.

Genome-wide association analysis identifies potential regulatory genes for eumelanin pigmentation in chicken plumage

Plumage color in chicken is determined by the proportion of eumelanin and pheomelanin pigmentation. As the main ingredient in plumage melanin, eumelanin plays a key role in the dark black, brown and grey coloration. However, very few studies have been performed to identify the related genes and mutations on a genome-wide scale. Herein, a resource family consisting of one backcross population and two F2 cross populations between a black roster and Yukou Brown I parent stockbreed was constructed for identification of genes related to eumelanin pigmentation. Chickens with eumelanin in their plumage were classified as the case group, and the rest were considered the control group. A genome-wide association study of this phenotype and genotypes using Affymetrix 600K HD SNP arrays in this F2 family revealed 13 significantly associated SNPs and in 10 separate genes on chromosomes 1, 2, 3 and 5. Based on previous studies in model species, we inferred that genes, including NUAK family kinase 1 (NUAK1) and sonic hedgehog (SHH), may play roles in the development of neural crest cells or melanoblasts during the embryonic period, which may also affect the eumelanin pigmentation. Our results facilitate the understanding of the genetic basis of eumelanin pigmentation in chicken plumage.

A gene expression study of dorso-ventrally restricted pigment pattern in adult fins of Neolamprologus meeli, an African cichlid species

Fish color patterns are among the most diverse phenotypic traits found in the animal kingdom. Understanding the molecular and cellular mechanisms that control in chromatophore distribution and pigmentation underlying this diversity is a major goal in developmental and evolutionary biology, which has predominantly been pursued in the zebrafish model system. Here, we apply results from zebrafish work to study a naturally occurring color pattern phenotype in the fins of an African cichlid species from Lake Tanganyika. The cichlid fish Neolamprologus meeli displays a distinct dorsal color pattern, with black and white stripes along the edges of the dorsal fin and of the dorsal half of the caudal fin, corresponding with differences in melanophore density. To elucidate the molecular mechanisms controlling the differences in dorsal and ventral color patterning in the fins, we quantitatively assessed the expression of 15 candidate target genes involved in adult zebrafish pigmentation and stripe formation. For reference gene validation, we screened the expression stability of seven widely expressed genes across the investigated tissue samples and identified tbp as appropriate reference. Relative expression levels of the candidate target genes were compared between the dorsal, striped fin regions and the corresponding uniform, grey-colored regions in the anal and ventral caudal fin. Dorso-ventral expression differences, with elevated levels in both white and black stripes, were observed in two genes, the melanosome protein coding gene pmel and in igsf11, which affects melanophore adhesion, migration and survival. Next, we predicted potential shared upstream regulators of pmel and igsf11. Testing the expression patterns of six predicted transcriptions factors revealed dorso-ventral expression difference of irf1 and significant, negative expression correlation of irf1 with both pmel and igsf11. Based on these results, we propose pmel, igsf11 and irf1 as likely components of the genetic mechanism controlling distinct dorso-ventral color patterns in N. meeli fins.

Agouti signalling protein (ASIP) gene: molecular cloning, sequence characterisation and tissue distribution in domestic goose

Agouti signalling protein (ASIP) is an endogenous antagonist of melanocortin-1 receptor (MC1R) and is involved in the regulation of pigmentation in mammals. The objective of this study was to identify and characterise the ASIP gene in domestic goose. The goose ASIP cDNA consisted of a 44-nucleotide 5'-terminal untranslated region (UTR), a 390-nucleotide open-reading frame (ORF) and a 45-nucleotide 3'-UTR. The length of goose ASIP genomic DNA was 6176 bp, including three coding exons and two introns. Bioinformatic analysis indicated that the ORF encodes a protein of 130 amino-acid residues with a molecular weight of 14.88 kDa and an isoelectric point of 9.73. Multiple sequence alignments and phylogenetic analysis showed that the amino-acid sequence of ASIP was conserved in vertebrates, especially in the avian species. RT-qPCR showed that the goose ASIP mRNA was differentially expressed in the pigment deposition tissues, including eye, foot, feather follicle, skin of the back, as well as in skin of the abdomen. The expression level of the ASIP gene in skin of the abdomen was higher than that in skin of the back. Those findings will contribute to further understanding the functions of the ASIP gene in geese plumage colouring.