Identification of polymorphisms in MITF and DCT genes and their associations with plumage colors in Asian duck breeds

Unique feather color characteristics and transcriptome analysis of hair follicles in Liancheng White ducks

Avian feather color is a fascinating trait, and the genetic mechanism of duck plumage formation is still in the preliminary stage. In this study, feather color of Liancheng White ducks was analyzed by determination of melanin content and RNA-seq analysis. In this research, 9 ducks from Mallards (n = 3), Liancheng White (n = 3) and Pekin ducks (n = 3) were used by high performance liquid chromatography (HPLC) and Masson-Fontana staining to reveal the difference of feather melanin content. RNA-seq from 11 hair follicle tissues (1- and 8-wk-old) of Liancheng White ducks (n = 5) and Pekin ducks (n = 7) was used to analyze the candidate genes for the feather melanin synthesis, and Immunofluorescence experiment was used to show the protein expression in 6 black- and white-feathered ducks. Pectorale, skin, liver, fat, brain, heart, kidney, lung, spleen of an 8-wk-old black-feathered Mallard were collected for candidate gene expression. The results showed that the contents of feathers, beak, web melanin in Liancheng White ducks were higher than in Pekin ducks (p < 0.05). Melanin within hair follicles was located in the barb ridge and hair matrix of black feather duck, also we found that TYRP1, TYR, SOX10 genes were differentially expressed between Liancheng White and Pekin ducks (p < 0.05), and these genes were mainly expressed showed in duck skin tissues. This study revealed the unique feather color phenotype of Liancheng White duck shedding light on the transcriptome that underlies it.

Identification of candidate genes associated with primary feathers of tianfu nonghua ducks based on Genome-wide association studies

The primary feathers of ducks have important economic value in the poultry industry. This study quantified the primary feather phenotype of Nonghua ducks, including the primary feathers' length, area, distribution of black spots, and feather symmetry. And genome-wide association analysis was used to screen candidate genes that affect the primary feather traits. The genome-wide association study (GWAS) results identified the genetic region related to feather length (FL) on chromosome 2. Through Linkage disequilibrium (LD) analysis, candidate regions (chr2: 115,246,393-116,501,448 bp) were identified and were further annotated to 5 genes: MRS2, GPLD1, ALDH5A1, KIAA0319, and ATP9B. Secondly, candidate regions related to feather black spots were identified on chromosome 21. Through LD analysis, the candidate regions (chr21: 163,552-2,183,853 bp) were screened and further annotated to 47 genes. Among them, STK4, CCN5, and YWHAB genes were related to melanin-related pathways or pigment deposition, which may be key genes affecting the distribution of black spots on feathers. In addition, we also screened 125 genes on multiple chromosomes that may be related to feather symmetry. Among them, significant SNPs on chromosome 1 were further identified as candidate regions (chr1: 142,118,209-142,223,605 bp) through LD analysis and annotated into 2 genes, TGFBRAP1 and LOC113839965. These results reported the genetic basis of the primary feather from multiple phenotypes, and offered valuable insights into the genetic basis for the growth and development of duck feathers and feather color pattern.

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.

Genome-Wide Association Study Reveals the Genetic Basis of Duck Plumage Colors

Plumage color is an artificially and naturally selected trait in domestic ducks. Black, white, and spotty are the main feather colors in domestic ducks. Previous studies have shown that black plumage color is caused by MC1R, and white plumage color is caused by MITF. We performed a genome-wide association study (GWAS) to identify candidate genes associated with white, black, and spotty plumage in ducks. Two non-synonymous SNPs in MC1R (c.52G>A and c.376G>A) were significantly related to duck black plumage, and three SNPs in MITF (chr13:15411658A>G, chr13:15412570T>C and chr13:15412592C>G) were associated with white plumage. Additionally, we also identified the epistatic interactions between causing loci. Some ducks with white plumage carry the c.52G>A and c.376G>A in MC1R, which also compensated for black and spotty plumage color phenotypes, suggesting that MC1R and MITF have an epistatic effect. The MITF locus was supposed to be an upstream gene to MC1R underlying the white, black, and spotty colors. Although the specific mechanism remains to be further clarified, these findings support the importance of epistasis in plumage color variation in ducks.

Integrated analysis of long non-coding RNAs and mRNA expression profiles identified potential interactions regulating melanogenesis in chicken skin

1. Long non-coding RNAs (lncRNAs) play important roles in various physiological functions. However, the mechanisms underlying the regulation of lncRNAs in melanogenesis remain unclear. To determine the molecular mechanisms involved in skin melanogenesis, the present study depicted the expression profiles of lncRNAs and messenger RNAs (mRNAs) in black- (B group) and white- (W group) skinned chickens using RNA sequencing.2. In total, 373 differentially expressed lncRNAs (DELs; 203 up-regulated and 170 down-regulated) and 253 differentially expressed genes (DEGs; 152 up-regulated and 101 down-regulated) were identified between the B and W groups. A total of eight known melanogenesis-related genes were identified (KIT, TYRP1, DCT (TYRP2), SLC45A2, OCA2, EDNRB2, TRPM1 and RAB38).3. Functional annotation of the co-expressed DEGs and DELs was performed using Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analyses. The co-expressed DEGs were mainly involved in melanogenesis and the co-expressed genes of 117 and 108 DELs were significantly enriched in the melanogenesis and tyrosine metabolism pathways, respectively.4. The DEL-DEG interaction network revealed that three lncRNAs (XR_003072387.1, XR_003075112.1, and XR_003077033.1) and DCT genes may have key roles in regulating melanogenesis in chicken skin. This data provides the groundwork for studying the lncRNA regulatory mechanisms of skin melanogenesis and suggested a new perspective on the modulation of melanogenesis in chicken skin based on a lncRNA-mRNA causal regulatory network.

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.

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.

Biallelic mutations in L-dopachrome tautomerase (DCT) cause infantile nystagmus and oculocutaneous albinism

Infantile nystagmus syndrome (INS) denominates early-onset, involuntary oscillatory eye movements with different etiologies. Nystagmus is also one of the symptoms in oculocutaneus albinism (OCA), a heterogeneous disease mainly caused by defects in melanin synthesis or melanosome biogenesis. Dopachrome tautomerase (DCT, also called TYRP2) together with tyrosinase (TYR) and tyrosin-related protein 1 (TYRP1) is one of the key enzymes in melanin synthesis. Although DCT´s role in pigmentation has been proven in different species, until now only mutations in TYR and TYRP1 have been found in patients with OCA. Detailed ophthalmological and orthoptic investigations identified a consanguineous family with two individuals with isolated infantile nystagmus and one family member with subtle signs of albinism. By whole-exome sequencing and segregation analysis, we identified the missense mutation c.176G > T (p.Gly59Val) in DCT in a homozygous state in all three affected family members. We show that this mutation results in incomplete protein maturation and targeting in vitro compatible with a partial or total loss of function. Subsequent screening of a cohort of patients with OCA (n = 85) and INS (n = 25) revealed two heterozygous truncating mutations, namely c.876C > A (p.Tyr292*) and c.1407G > A (p.Trp469*), in an independent patient with OCA. Taken together, our data suggest that mutations in DCT can cause a phenotypic spectrum ranging from isolated infantile nystagmus to oculocutaneous albinism.

Characterization and Comparative Transcriptomic Analysis of Skeletal Muscle in Pekin Duck at Different Growth Stages Using RNA-Seq

Simple Summary

Skeletal muscle is an important tissue and its development is strictly regulated by genes. In this study, in order to understand the muscle-related gene expression in Pekin duck, RNA-seq was performed to analyze and compare skeletal muscle at different growth stages. Alternative splicing, single nucleotide polymorphisms and insertion–deletions were detected, and 299 novel genes were discovered. MYL4, IGF2BP1, CSRP3, SPP1, KLHL31, LAMB2, LAMA2, ITGB1 and OPN played crucial roles in skeletal muscle development. Oxidative phosphorylation, ECM-receptor interaction, focal adhesion, carbon metabolism, and biosynthesis of amino acids participated in the regulation of skeletal muscle development in Pekin duck. This study provides an important reference for revealing the developmental mechanisms of pectoral and leg muscles in duck.

Abstract

Skeletal muscle, accounting for approximately 50% of body weight, is the largest and most important tissue. In this study, the gene expression profiles and pathways in skeletal muscle of Pekin duck were investigated and compared at embryonic day 17, 21, and 27 and postnatally at 6 months of age. An average of 49,555,936 reads in each sample was obtained from the transcriptome libraries. Over 70.0% of alternative splicing (AS) in each sample was mainly alternative 5′ first exon (transcription start site)—the first exon splicing (TSS) and alternative 3′ last exon (transcription terminal site)—the last exon splicing (TTS), indicating that TSS and TTS were the most common AS event in Pekin ducks, and these AS events were closely related to the regulation of muscle development at different growth stages. The results provided a valuable genomic resource for selective breeding and functional studies of genes. A total of 299 novel genes with ≥2 exons were obtained. There were 294 to 2806 differentially expressed genes (DEGs) in each pairwise comparison of Pekin duck. Notably, 90 DEGs in breast muscle and 9 DEGs in leg muscle were co-expressed at all developmental points. DEGs were validated by qPCR analysis, which confirmed the tendency of the expression. DEGs related to muscle development were involved in biological processes such as “endodermal cell differentiation”, “muscle cell cellular homeostasis”, “skeletal muscle tissue growth” and “skeletal muscle cell differentiation”, and were involved in pathways such as oxidative phosphorylation, ECM-receptor (extracellular matrix receptor) interaction, focal adhesion, carbon metabolism, and biosynthesis of amino acids. Some DEGs, including MYL4, IGF2BP1, CSRP3, SPP1 and KLHL31, as well as LAMB2, LAMA2, ITGB1 and OPN, played crucial roles in muscle growth and development. This study provides valuable information about the expression profile of mRNAs and pathways from duck skeletal muscle at different growth stages, and further functional study of these mRNAs and pathways could provide new ideas for studying the molecular networks of growth and development in duck skeletal muscle.

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.

Integrated Analysis of mRNA Expression, CpG Island Methylation, and Polymorphisms in the MITF Gene in Ducks (Anas platyrhynchos)

Microphthalmia-associated transcription factor (MITF) is a key regulator for the development and function of melanocytes in skin, eye, and plumage pigmentations. Thus, the MITF was selected as a candidate gene associated with plumage coloration in ducks. This study analyzed the mRNA expression, promoter methylation, and polymorphisms in the MITF gene in ducks with different plumage colors (Putian Black, Putian White, Liancheng White, and Longsheng Jade-green). No expression of the MITF melanin-specific isoform (MITF-M) was detected in white feather bulbs. By contrast, the mRNA expression levels of MITF-M were high in black feather bulbs. Bioinformatics analysis showed that two CpG islands were present in the promoter region of the MITF gene. The methylation level of the second CpG island was significantly lower in black feather bulbs than in white feather bulbs. However, the methylation level of the first CpG island was not different among the feather bulbs with various colors except Liancheng White feather bulbs. The methylation status of the whole CpG island significantly and negatively correlated with the mRNA expression of MITF-M (P < 0.05). Furthermore, four novel SNPs (single nucleotide polymorphisms) were identified in the 5′UTR, exon 4, intron 7, and intron 8 of the MITF gene. Allele T in g.39807T>G and allele G in g.40862G>A were the predominant alleles only found in Putian White, whereas the variant A allele in g.32813G>A exhibited a high allele frequency in Liancheng White. Collectively, these results contributed to the understanding of the function of the MITF gene in duck plumage coloration.

The Long Noncoding RNA UCA1 Negatively Regulates Melanogenesis in Melanocytes

The long noncoding RNA UCA1 was first discovered in bladder cancer and is known to regulate the proliferation and migration of melanoma. However, its role in melanogenesis is unclear. In this study, we aimed to explore the role and mechanism of UCA1 in melanogenesis. Our findings showed that the expression of UCA1 was negatively correlated with melanin content in melanocytes and pigmented nevus. Overexpression of UCA1 in melanocytes decreased melanin content and the expression of melanogenesis-related genes, whereas knockdown of UCA1 in melanocytes had the opposite effect. High-throughput sequencing revealed that microphthalmia-associated transcription factor (MITF), an important transcription factor affecting melanogenesis, was also negatively correlated with the expression of UCA1. Furthermore, the transcription factor CRE-binding protein (CREB), which promotes MITF expression, was negatively regulated by UCA1. The cAMP/protein kinase A (PKA), extracellular signal-regulated kinase (ERK), and c-Jun N-terminal kinase (JNK) signaling pathways, which are upstream of the CREB/MITF/melanogenesis axis, were activated or inhibited in response to silencing or enhancing UCA1 expression, respectively. In addition, enhanced UCA1 expression downregulates the expression of melanogenesis-related genes induced by UVB in melanocytes. In conclusion, UCA1 may negatively regulate the CREB/MITF/melanogenesis axis through inhibiting the cAMP/PKA, ERK, and JNK signaling pathways in melanocytes. UCA1 may be a potential therapeutic target for the treatment of pigmented skin diseases.

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.

Whole-genome resequencing reveals signatures of selection and timing of duck domestication

Background

The genetic basis of animal domestication remains poorly understood, and systems with substantial phenotypic differences between wild and domestic populations are useful for elucidating the genetic basis of adaptation to new environments as well as the genetic basis of rapid phenotypic change. Here, we sequenced the whole genome of 78 individual ducks, from two wild and seven domesticated populations, with an average sequencing depth of 6.42X per individual.

Results

Our population and demographic analyses indicate a complex history of domestication, with early selection for separate meat and egg lineages. Genomic comparison of wild to domesticated populations suggests that genes that affect brain and neuronal development have undergone strong positive selection during domestication. Our F_ST analysis also indicates that the duck white plumage is the result of selection at the melanogenesis-associated transcription factor locus.

Conclusions

Our results advance the understanding of animal domestication and selection for complex phenotypic traits.