Identification of genes related to white and black plumage formation by RNA-Seq from white and black feather bulbs in ducks

Regulatory gene network for coffee-like color morph of TYRP1 mutant of oujiang color common carp

BACKGROUND

Neither a TYRP1-mediated highly conserved genetic network underlying skin color towards optimum defense nor the pathological tendency of its mutation is well understood. The Oujiang Color Common Carp (Cyprinus carpio var. color) as a model organism, offering valuable insights into genetics, coloration, aquaculture practices, and environmental health. Here, we performed a comparative skin transcriptome analysis on TYRP1 mutant and wild fishes by applying a conservative categorical approach considering different color phenotypes.

RESULTS

Our results reveal that an unusual color phenotype may be sensitized with TYRP1 mutation as a result of upregulating several genes related to an anti-inflammatory autoimmune system in response to the COMT-mediated catecholamine neurotransmitters in the skin. Particularly, catecholamines-derived red/brown, red with blue colored membrane attack complex, and brown/grey colored reduced eumelanin are expected to be aggregated in the regenerated cells.

CONCLUSIONS

It is, thus, concluded that the regenerated cells with catecholamines, membrane attack complex, and eumelanin altogether may contribute to the formation of the unusual (coffee-like) color phenotype in TYRP1 mutant.

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.

Investigation into the association of FNDC1 and ADAMTS12 gene expression with plumage coloration in Muscovy ducks

To elucidate the molecular genetic mechanisms underpinning feather color in Muscovy ducks. A cohort of 100 Muscovy ducks was meticulously selected for this research. Follicular tissues from ducks exhibiting black and white plumage served as the experimental samples. From these tissues, RNA and proteins were extracted for further analysis. The RNA underwent reverse transcription polymerase chain reaction amplification, followed by validation through western blot assays. The data revealed a significant upregulation in the expression of FN domain-containing protein 1 (FNDC1) and ADAMTS12 genes in Muscovy ducks with white plumage traits as opposed to those with black plumage traits. Specifically, individuals with pure white plumage demonstrated a markedly elevated expression of the FNDC1 gene in comparison to their pure black counterparts. Conversely, expression levels of the ADAMTS12 gene were found to be reduced in ducks with pure black plumage relative to those with pure white plumage. Notably, the expression patterns of FNDC1 and ADAMTS12 genes exhibited inconsistencies between mRNA and protein levels. This study offers significant insights into the molecular genetic mechanisms underlying feather color variation in Muscovy ducks. FNDC1 and ADAMTS12 could be considered potential targets for genetic manipulation or selective breeding strategies aimed at achieving specific feather color phenotypes in Muscovy ducks.

The relationship of MITF gene expression and promoter methylation with plumage colour in quail

1. This study focused on the relationship between MITF mRNA expression and plumage colour in quail and the effect of promoter methylation on the expression of MITF mRNA.2. The CDS region of MITF mRNA was cloned by RT-PCR, followed by DNA sequencing. The RT-qPCR method was used to analyse the expression levels of MITF mRNA in dorsal skin tissue in Korean quail and Beijing white quail. The promoter region of the MITF gene was cloned, and the CpG island was predicted by the CpGplot program. The methylation levels of the CpG island were analysed using BS-PCR technology.3. Quail MITF mRNA contains a 1,476 bp complete ORF, which encodes a 492 amino acid residue protein. The MITF protein has no signal peptide or transmembrane region. The expression of MITF mRNA in dorsal tissue of Korean quail was significantly higher than that in Beijing white quail (p < 0.01). Abundant cis-elements and a 346 bp CpG island were found in the promoter region of the MITF gene. The average methylation level of the CpG island was 22 (22%) in Korean quail, and 46 (30%) in Beijing white quail (p < 0.05).4. The hypermethylation of the MITF gene promoter region in Beijing white quail resulted in a decrease in expression level, which was related to white feather colour.

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.

Genome-Wide Association Studies Reveal Candidate Genes Associated with Pigmentation Patterns of Single Feathers of Tianfu Nonghua Ducks

In modern advanced genetics and breeding programs, the study of genes related to pigmentation in ducks is gaining much attention and popularity. Genes and DNA mutation cause variations in the plumage color traits of ducks. Therefore, discovering related genes responsible for different color traits and pigment patterns on each side of the single feathers in Chinese ducks is important for genetic studies. In this study, we collected feather images from 340 ducks and transported them into Image Pro Plus (IPP) 6.0 software to quantify the melanin content in the feathers. Thereafter, a genome-wide association study was conducted to reveal the genes responsible for variations in the feather color trait. The results from this study revealed that the pigmented region was larger in the male ducks as compared to the female ducks. In addition, the pigmented region was larger on the right side of the feather vane than on the left side in both dorsal and ventral feathers, and a positive correlation was observed among the feather color traits. Further, among the annotated genes, WNT3A, DOCK1, RAB1A, and ALDH1A3 were identified to play important roles in the variation in pigmented regions of the various feathers. This study also revealed that five candidate genes, including DPP8, HACD3, INTS14, SLC24A1, and DENND4A, were associated with the color pigment on the dorsal feathers of the ducks. Genes such as PRKG1, SETD6, RALYL, and ZNF704 reportedly play important roles in ventral feather color traits. This study revealed that genes such as WNT3A, DOCK1, RAB1A, and ALDH1A3 were associated with different pigmentation patterns, thereby providing new insights into the genetic mechanisms of single-feather pigmentation patterns in ducks.

The draft genome of the Temminck's tragopan (Tragopan temminckii) with evolutionary implications

BACKGROUND

High-quality genome data of birds play a significant role in the systematic study of their origin and adaptive evolution. The Temminck's tragopan (Tragopan temminckii) (Galliformes, Phasianidae), a larger pheasant, is one of the most abundant and widely distributed species of the genus Tragopan, and was defined as class II of the list of national key protected wild animals in China. The absence of a sequenced genome has restricted previous evolutionary trait studies of this taxa.

RESULTS

The whole genome of the Temminck's tragopan was sequenced using Illumina and PacBio platform, and then de novo assembled and annotated. The genome size was 1.06 Gb, with a contig N50 of 4.17 Mb. A total of 117.22 Mb (11.00%) repeat sequences were identified. 16,414 genes were predicted using three methods, with 16,099 (98.08%) annotated as functional genes based on five databases. In addition, comparative genome analyses were conducted across 12 Galliformes species. The results indicated that T. temminckii was the first species to branch off from the clade containing Lophura nycthemera, Phasianus colchicus, Chrysolophus pictus, Syrmaticus mikado, Perdix hodgsoniae, and Meleagris gallopavo, with a corresponding divergence time of 31.43 million years ago (MYA). Expanded gene families associated with immune response and energy metabolism were identified. Genes and pathways associated with plumage color and feather development, immune response, and energy metabolism were found in the list of positively selected genes (PSGs).

CONCLUSIONS

A genome draft of the Temminck's tragopan was reported, genome feature and comparative genome analysis were described, and genes and pathways related to plumage color and feather development, immune response, and energy metabolism were identified. The genomic data of the Temminck's tragopan considerably contribute to the genome evolution and phylogeny of the genus Tragopan and the whole Galliformes species underlying ecological adaptation strategies.

Embryonic Development of Parthenogenetic and Sexual Eggs in Lower Termites

Worldwide, termites are one of few social insects. In this research, the stages of embryonic development in the parthenogenetic and sexual eggs of Reticulitermes aculabialis and R. flaviceps were observed and described. In R. flaviceps, the egg development of the FF and FM groups happened during the early phases of development, whereas in R. aculabialis, this appeared mainly during the late phase of development. The variance in the number of micropyles between the R. flaviceps FF colony type and the R. aculabialis FF colony type was statistically significant. Five stages of egg development were found in both types of R. aculabialis but only the sexual eggs of R. flaviceps. In R. flaviceps, 86% of the parthenogenetic eggs stopped growing during the blastoderm development, with the yolk cell assembling frequently in the center of the egg. According to the results of the single-cell transcriptome sequencing, we investigated the egg-to-larval expression level of genes (pka, map2k1, mapk1/3, hgk, mkp, and pax6) and indicated that the levels of essential gene expression in RaFF were considerably higher than in RfFF (p < 0.05). We also discovered that the oocyte cleavage rate in the FF colony type was considerably lower in R. flaviceps compared to R. aculabialis, which gave rise to a smaller number of mature oocytes in R. flaviceps. During ovulation in both species, oocytes underwent activation and one or two cleavage events, but the development of unfertilized eggs ceased in R. flaviceps. It was shown that termite oocyte and embryonic development were heavily influenced by genes with significant expressions. Results from the databases KEGG, COG, and GO unigenes revealed the control of numerous biological processes. This study is the first to complete a database of parthenogenetic and sexual eggs of R. flaviceps and R. aculabialis.

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.

Coordinated microRNA/mRNA Expression Profiles Reveal Unique Skin Color Regulatory Mechanisms in Chinese Giant Salamander (Andrias davidianus)

The Chinese giant salamander (Andrias davidianus) has been increasingly popular in the aquaculture market in China in recent years. In the breeding process of Andrias davidianus, we found that some albino individuals were extremely rare and could not be inherited stably, which severely limits their commercialization in the aquaculture market. In this study, we performed transcriptome and small RNA (sRNA) sequencing analyses in the skin samples of wild-type (WT) and albino (AL) Andrias davidianus. In total, among 5517 differentially expressed genes (DEGs), 2911 DEGs were down-regulated in AL, including almost all the key genes involved in melanin formation. A total of 25 miRNAs were differentially expressed in AL compared to WT, of which 17 were up-regulated. Through the integrated analysis, no intersection was found between the target genes of the differentially expressed miRNAs and the key genes for melanin formation. Gene Ontology (GO) and KEGG pathway analyses on DEGs showed that these genes involved multiple processes relevant to melanin synthesis and the key signal pathway MAPK. Interestingly, the transcription factors SOX10 and PAX3 and the Wnt signaling pathway that play a key role in other species were not included, while the other two transcription factors in the SOX family, SOX21 and SOX7, were included. After analyzing the key genes for melanin formation, it was interesting to note an alternative splicing form of the MITF in WT and a critical mutation of the SLC24A5 gene in AL, which might be the main reason for the skin color change of Andrias davidianus. The results contributed to understanding the molecular mechanism of skin pigmentation in Andrias davidianus and accelerating the acquisition process of individuals with specific body colors by genetic means.

Association analysis of melanophilin (MLPH) gene expression and polymorphism with plumage color in quail

We explore the relationship between the melanophilin (MLPH) gene and quail plumage color and provide a reference for subsequent quail plumage color breeding. In this experiment, real-time quantitative PCR (RT-qPCR) technology was used to analyze the relative mRNA expression levels of Korean quail (maroon) and Beijing white quail embryos at different developmental stages. Two single-nucleotide polymorphisms (SNPs) in the MLPH gene were screened based on the RNA-sequencing (RNA-Seq) data of skin tissues of Korean quail and Beijing white quail during the embryonic stage. Kompetitive allele-specific PCR (KASP) technology was used for genotyping in the resource population, and correlation analysis was carried out with the plumage color traits of quail. Finally, bioinformatics was used to predict the effects of these two SNPs on the structure and function of the encoded protein. The results showed that the expression level of the MLPH gene during embryonic development of Beijing white quail was significantly higher than that of Korean quail ( P < 0.01 ). The frequency distribution of the three genotypes (CC, CA and AA) of the Beijing white quail at the c.1807C  >  A mutation site was significantly different from that of the Korean quail ( P < 0.01 ). The frequency distribution of the three genotypes (GG, GA and AA) of the Beijing white quail at the c.2129G  >  A mutation site was significantly different from that of the Korean quail ( P < 0.01 ). And there was a significant correlation between the c.1807C  >  A mutation site and the white plumage phenotype. Bioinformatics showed that SNP1 (c.1807C  >  A) was a neutral mutation and that SNP2 (c.2129G  >  A) was a deleterious mutation. The prediction of protein conservation showed that the mutation sites of coding proteins R603S and G710D caused by SNP1 (c.1807C  >  A) and SNP2 (c.2129G  >  A) were highly conserved.

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.

Involvement of B-aat1 and Cbs in regulating mantle pigmentation in the Pacific oyster (Crassostrea gigas)

BACKGROUND

Shell color formation is an important physiological process in bivalves, the molecular genetic basis has potential application in bivalve aquaculture, but there is still remaining unclear about this issue. The cystine/glutamate transporter (Slc7a11) and cystathionine beta-synthase (Cbs) are integral genes in pheomelanin synthesis pathway, which is vital to skin pigmentation.

METHODS AND RESULTS

Here, the sequences of b (0, +) -type amino acid transporter 1 (B-aat1) and Cbs in Pacific oyster (Crassostrea gigas) (CgB-aat1, CgCbs) were characterized. Phylogenetically, the deduced amino acid sequences of CgB-aat1 and CgCbs both possessed conserved features. Genes were both ubiquitously expressed in six tested tissues with more abundant expression level in central mantle. Besides, the polyclonal antibodies of CgB-aat1, CgCbs, CgTyr, and CgTyrp2 were successfully prepared. Immunofluorescence analysis revealed that CgB-aat1 and CgCbs proteins were both expressed in gill rudiments of eyed-larvae and concentrated mainly in cytoplasm of epithelial cell and nerve axons in mantle. Additionally, after CgB-aat1 or CgCbs silencing, expressions at mRNA and protein levels of CgB-aat1 and CgCbs involved in pheomelanin synthesis were significantly suppressed, and CgTyr, CgTyrp1 and CgTyrp2 related to eumelanin synthesis were also down-regulated but no apparent differences, respectively. Moreover, micrographic examination found less brown-granules at mantle edge in CgB-aat1 interference group.

CONCLUSION

These results implied that pheomelanin synthesis was possible induced by CgB-aat1-CgTyr-CgCbs axis, and it played an essential role on mantle pigmentation in the oysters. These findings provide the useful genetic knowledge and enrich the physiological information for the shell color formation in bivalve aquaculture.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Draft genome assemblies of four manakins

Manakins are a family of small suboscine passerine birds characterized by their elaborate courtship displays, non-monogamous mating system, and sexual dimorphism. This family has served as a good model for the study of sexual selection. Here we present genome assemblies of four manakin species, including Cryptopipo holochlora, Dixiphia pipra (also known as Pseudopipra pipra), Machaeropterus deliciosus and Masius chrysopterus, generated by Single-tube Long Fragment Read (stLFR) technology. The assembled genome sizes ranged from 1.10 Gb to 1.19 Gb, with average scaffold N50 of 29 Mb and contig N50 of 169 Kb. On average, 12,055 protein-coding genes were annotated in the genomes, and 9.79% of the genomes were annotated as repetitive elements. We further identified 75 Mb of Z-linked sequences in manakins, containing 585 to 751 genes and an ~600 Kb pseudoautosomal region (PAR). One notable finding from these Z-linked sequences is that a possible Z-to-autosome/PAR reversal could have occurred in M. chrysopterus. These de novo genomes will contribute to a deeper understanding of evolutionary history and sexual selection in manakins.

Measurement(s)	whole genome sequencing
Technology Type(s)	BGISEQ-500 Sequencing
Sample Characteristic - Organism	Pipridae

The evolutionary history and mechanistic basis of female ornamentation in a tropical songbird

Ornamentation, such as the showy plumage of birds, is widespread among female vertebrates, yet the evolutionary pressures shaping female ornamentation remain uncertain. In part this is due to a poor understanding of the mechanistic route to ornamentation in females. To address this issue, we evaluated the evolutionary history of ornament expression in a tropical passerine bird, the White-shouldered Fairywren, whose females, but not males, strongly vary between populations in occurrence of ornamented black-and-white plumage. We first use phylogenomic analysis to demonstrate that female ornamentation is derived and that female ornamentation evolves independently of changes in male plumage. We then use exogenous testosterone in a field experiment to induce partial ornamentation in naturally unornamented females. By sequencing the transcriptome of experimentally induced ornamented and natural feathers, we identify genes expressed during ornament production and evaluate the degree to which female ornamentation in this system is associated with elevated testosterone, as is common in males. We reveal that some ornamentation in females is linked to testosterone and that sexes differ in ornament-linked gene expression. Lastly, using genomic outlier analysis we identify a candidate melanogenesis gene that lies in a region of high genomic divergence among populations that is also differentially expressed in feather follicles of different female plumages. Taken together, these findings are consistent with sex-specific selection favoring the evolution of female ornaments and demonstrate a key role for testosterone in generating population divergence in female ornamentation through gene regulation. More broadly, our work highlights similarities and differences in how ornamentation evolves in the sexes.

Genome-Wide Analysis Identifies Candidate Genes Encoding Feather Color in Ducks

Comparative population genomics and genome-wide association studies (GWAS) offer opportunities to discover human-driven detectable signatures within the genome. From the point of view of evolutionary biology, the identification of genes associated with the domestication of traits is of interest for the elucidation of the selection of these traits. To this end, an F2 population of ducks, consisting of 275 ducks, was genotyped using a whole genome re-sequence containing 12.6 Mb single nucleotide polymorphisms (SNPs) and four plumage colors. GWAS was used to identify the candidate and potential SNPs of four plumage colors in ducks (white, spot, grey, and black plumage). In addition, FST and genetic diversity (π ratio) were used to screen signals of the selective sweep, which relate to the four plumage colors. Major genomic regions associated with white, spotted, and black feathers overlapped with their candidate selection regions, whereas no such overlap was observed with grey plumage. In addition, MITF and EDNRB2 are functional candidate genes that contribute to white and black plumage due to their indirect involvement in the melanogenesis pathway. This study provides new insights into the genetic factors that may influence the diversity of plumage color.

Sex Identification of Feather Color in Geese and the Expression of Melanin in Embryonic Dorsal Skin Feather Follicles

In production practice, we have found that the gray and black down on the backs of the Holdobaggy goslings is usually darker in females than in males. Melanin is the key pigment affecting the color of poultry plumage. Therefore, to determine whether the darkness of the dorsal plumage of the Holdobaggy goslings is related to sex, we study the melanin in the feather follicles of the dorsal skin during the embryonic period. The feather follicle structure and melanin distribution on the dorsal surface of the goose embryo is observed by HE staining and melanin-specific staining. The melanin content in the feather follicles of the dorsal skin of goslings is determined by ELISA. The results showed that the melanin content is higher in female geese than in males (p < 0.05). In addition, we also analyze the mRNA and protein expression levels of melanin-related genes (TYRP1 and ASIP) by quantitative real-time PCR and Western blotting analysis. The results show that the mRNA expression level of TYRP1 is significantly higher in the females’ dorsal skin feather follicles (p < 0.05), while the mRNA expression level of ASIP is significantly higher in the dorsal skin feather follicles of male geese (p < 0.05). In conclusion, the difference between males and females in the color of the black feathers on the dorsal track of the Holdobaggy goslings is verified, and it is feasible to identify the sex by the initial plumage color.

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.

Transcriptome Analysis Reveals Genes Associated With Sexual Dichromatism of Head Feather Color in Mallard

Sexual dimorphism of feather color is typical in mallards, in which drakes exhibit green head feathers, while females show dull head feather color. We showed that more melanosomes deposited in the males' head's feather barbules than females and further form a two-dimensional hexagonal lattice, which conferred the green feather coloration of drakes. Additionally, transcriptome analysis revealed that some essential melanin biosynthesis genes were highly expressed in feather follicles during the development of green feathers, contributing to melanin deposition. We further identified 18 candidate differentially expressed genes, which may affect the sharp color differences between the males' head feathers, back feathers, and the females' head feathers. TYR and TYRP1 genes are associated with melanin biosynthesis directly. Their expressions in the males' head feather follicles were significantly higher than those in the back feather follicles and females' head feather follicles. Most clearly, the expression of TYRP1 was 256 and 32 times higher in the head follicles of males than in those of the female head and the male back, respectively. Hence, TYR and TYRP1 are probably the most critical candidate genes in DEGs. They may affect the sexual dimorphism of head feather color by cis-regulation of some transcription factors and the Z-chromosome dosage effect.

Comparative Transcriptome Analysis Identifies Candidate Genes Related to Black-Spotted Pattern Formation in Spotted Scat (Scatophagus argus)

Simple Summary

Spotted scat (Scatophagus argus) is a commercially important marine aquaculture and ornamental fish species in China and East Asian countries. There are dozens of black spots on each side of the body, and the caudal fin, which is yellow and black, is appreciated in ornamental fish markets. To explore the genetic mechanisms of its pattern formation, we found 2357 differentially expressed genes (DEGs) by comparing the transcriptome in the black-spotted skin, non-spotted skin and caudal fin in S. argus. The results will expand our knowledge about the molecular mechanism of important genes and pathways associated with pigment pattern formation and provide a certain theoretical basis for the molecular breeding in S. argus.

Abstract

Spotted scat (Scatophagus argus) is an economically important marine aquaculture and ornamental fish species in Asia, especially in southeast China. In this study, skin transcriptomes of S. argus were obtained for three types of skin, including black-spotted skin (A), non-spotted skin (B) and caudal fin (C). A total of nine complementary DNA (cDNA) libraries were obtained by Illumina sequencing. Bioinformatics analysis revealed that 1358, 2086 and 487 genes were differentially expressed between A and B, A and C, and B and C, respectively. The results revealed that there were 134 common significantly differentially expressed genes (DEGs) and several key genes related to pigment synthesis and pigmentation, including tyrp1, mitf, pmel, slc7a2, tjp1, hsp70 and mart-1. Of these, some DEGs were associated with pigmentation-related Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, such as tyrosine metabolism, melanogenesis, the Wnt signaling pathway and the mitogen-activated protein kinase (MAPK) signaling pathway. The results will facilitate understanding the molecular mechanisms of skin pigmentation differentiation in S. argus and provide valuable information for skin coloration, especially the formation of spotted patterns on other marine fish species.

Analysis of MC1R, MITF, TYR, TYRP1, and MLPH Genes Polymorphism in Four Rabbit Breeds with Different Coat Colors

Simple Summary

Coat color is an important breed characteristic and economic trait for rabbits, and it is regulated by a few genes. In this study, the gene frequencies of some pigmentation genes were investigated in four Chinese native rabbit breeds with different coat colors. A total of 14 genetic variants were detected in the gene fragments of MC1R, MITF, TYR, TYRP1, and MLPH genes, and there was low-to-moderate polymorphism in the populations. The gene frequency showed significant differences among the four rabbit populations. The above results suggest that these genetic variations play an important role in regulating the coat color of rabbits. This study will provide potential molecular markers for the breeding of coat color traits in rabbits.

Abstract

Pigmentation genes such as MC1R, MITF, TYR, TYRP1, and MLPH play a major role in rabbit coat color. To understand the genotypic profile underlying coat color in indigenous Chinese rabbit breeds, portions of the above-mentioned genes were amplified and variations in them were analyzed by DNA sequencing. Based on the analysis of 24 Tianfu black rabbits, 24 Sichuan white rabbits, 24 Sichuan gray rabbits, and 24 Fujian yellow rabbits, two indels in MC1R, three SNPs in MITF, five SNPs (single nucleotide polymorphisms) in TYR, one SNP in TYRP1, and three SNPs in MLPH were discovered. These variations have low-to-moderate polymorphism, and there are significant differences in their distribution among the different breeds (p < 0.05). These results provide more information regarding the genetic background of these native rabbit breeds and reveal their high-quality genetic resources.

Comparative transcriptomics reveals the molecular genetic basis of pigmentation loss in Sinocyclocheilus cavefishes

Cave-dwelling animals evolve distinct troglomorphic traits, such as loss of eyes, skin pigmentation, and augmentation of senses following long-term adaptation to perpetual darkness. However, the molecular genetic mechanisms underlying these phenotypic variations remain unclear. In this study, we conducted comparative histology and comparative transcriptomics study of the skin of eight Sinocyclocheilus species (Cypriniformes: Cyprinidae) that included surface- and cave-dwelling species. We analyzed four surface and four cavefish species by using next-generation sequencing, and a total of 802,798,907 clean reads were generated and assembled into 505,495,009 transcripts, which contributed to 1,037,334 unigenes. Bioinformatic comparisons revealed 10,629 and 6,442 significantly differentially expressed unigenes between four different surface-cave fish groups. Further, tens of differentially expressed genes (DEGs) potentially related to skin pigmentation were identified. Most of these DEGs (including GNAQ, PKA, NRAS, and p38) are downregulated in cavefish species. They are involved in key signaling pathways of pigment synthesis, such as the melanogenesis, Wnt, and MAPK pathways. This trend of downregulation was confirmed through qPCR experiments. This study will deepen our understanding of the formation of troglomorphic traits in cavefishes.

Transcriptome analysis of feather follicles reveals candidate genes and pathways associated with pheomelanin pigmentation in chickens

Yellow plumage is common in chickens, especially in breeds such as the Huiyang Bearded chicken, which is indigenous to China. We evaluated plumage colour distribution in F1, F2, and F3 populations of an Huiyang Bearded chicken × White Leghorn chicken cross, the heredity of the yellow plumage trait was distinguished from that of the gold plumage and other known plumage colours. Microscopic analysis of the feather follicles indicated that pheomelanin particles were formed in yellow but not in white feathers. To screen genes related to formation of the pheomelanin particles, we generated transcriptome data from yellow and white feather follicles from 7- and 11-week-old F3 chickens using RNA-seq. We identified 27 differentially expressed genes (DEGs) when comparing the yellow and white feather follicles. These DEGs were enriched in the Gene Ontology classes ‘melanosome’ and ‘melanosome organization’ related to the pigmentation process. Down-regulation of TYRP1, DCT, PMEL, MLANA, and HPGDS, verified using quantitative reverse transcription PCR, may lead to reduced eumelanin and increased pheomelanin synthesis in yellow plumage. Owing to the presence of the Dominant white locus, both white and yellow plumage lack eumelanin, and white feathers showed no pigments. Our results provide an understanding of yellow plumage formation in chickens.

Identification of key microRNAs affecting melanogenesis of breast muscle in Muchuan black-boned chickens by RNA sequencing

1. Melanin content is considered an important indicator of meat quality in black-boned chickens, which have a high market value. To understand the complex physiological processes underlying muscle melanogenesis in this chicken, differentially expressed miRNAs (DEMs) were detected between black muscle (BM) and white muscle (WM) of chickens using high-throughput sequencing technology. Six small RNA libraries were constructed, and more than 16.75 million clean reads were obtained for each library. 2. A total of 582 known miRNAs and 65 novel miRNAs were identified from the six chicken sequence libraries. A total of 19 DEMs were identified between the two groups, of which nine were upregulated and 10 were downregulated. Furthermore, the DEMs were predicted to target 572 genes. 3. Certain DEMs (such as miR-204, miR-133b, and miR-12 229-3p) and their target genes may play an important role in muscle melanogenesis of chickens. These findings provide a foundation for clarifying the miRNA regulatory mechanisms involved in muscle pigmentation in avian species.

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.

Morphological Characteristics and Comparative Transcriptome Analysis of Three Different Phenotypes of Pristella maxillaris

Pristella maxillaris is known as the X-ray fish based on its translucent body. However, the morphological characteristics and the molecular regulatory mechanisms of these translucent bodies are still unknown. In this study, the following three phenotypes, a black-and-gray body color or wild-type (WT), a silvery-white body color defined as mutant I (MU1), and a fully transparent body with a visible visceral mass named as mutant II (MU2), were investigated to analyze their chromatophores and molecular mechanisms. The variety and distribution of pigment cells in the three phenotypes of P. maxillaris significantly differed by histological assessment. Three types of chromatophores (melanophores, iridophores, and xanthophores) were observed in the WT, whereas MU1 fish were deficient in melanophores, and MU2 fish lacked melanophores and iridophores. Transcriptome sequencing of the skin and peritoneal tissues of P. maxillaris identified a total of 166,089 unigenes. After comparing intergroup gene expression levels, more than 3,000 unigenes with significantly differential expression levels were identified among three strains. Functional annotation and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the differentially expressed genes (DEGs) identified a number of candidates melanophores and iridophores genes that influence body color. Some DEGs that were identified using transcriptome analysis were confirmed by quantitative real-time PCR. This study serves as a global survey of the morphological characteristics and molecular mechanism of different body colors observed in P. maxillaris and thus provides a valuable theoretical foundation for the molecular regulation of the transparent phenotype.

Integrated analysis of RNA-seq and microRNA-seq depicts miRNA-mRNA networks involved in stripe patterns of Botia superciliaris skin

Botia superciliaris, an endemic cobitid fish in China, is widely accepted by Chinese consumers because its edibility. Recently, the black and yellow stripes of B. superciliaris skin have made this species increasingly popular as a novel ornamental fish. However, the genetic basis of the stripe patterns in B. superciliaris skin has not been extensively studied. In this study, Illumina sequencing was employed to identify the mRNAs and miRNAs involved in stripe pattern formation in B. superciliaris skin. A total of 147.25 and 155.15 million (M) high-quality transcriptome reads were generated from three black and yellow skin libraries respectively, which resulted in 159,327 unigenes that were used as reference sequences. A total of 3169 genes exhibited significantly differential expression patterns (fold-change ≥ 2 or ≤ 0.5 and q ≤ 0.05), including 1891 upregulated genes (59.67%) and 1278 downregulated genes (40.33%) in black vs yellow skin. These genes were enriched in 50 GO terms and 10 KEGG pathways (q ≤ 0.05), including melanogenesis, with 21 upregulated genes and 5 downregulated genes in black vs yellow skin. Based on the zebrafish genome, miRNA-seq identified a total of 355 miRNAs, which included 38 novel miRNAs. Furthermore, 87 differentially expressed miRNAs including 50 upregulated and 37 downregulated miRNAs were identified in different color skin (fold-change ≥ 2 or ≤ 0.5 and q ≤ 0.05). Then, target prediction revealed a variety of putative target genes; differentially expressed mRNAs and miRNAs patterns of high-throughput sequencing were validated in 5 mRNAs and miR-217-5p by qRT-PCR. In vivo tests and dual-luciferase reporter assay revealed that overexpression of miR-217-5p can inhibit pheomelanin formation by targeting Zgc. In this study, a comparative analysis was conducted to profile the transcriptome of black and yellow skin for B. superciliaris, and we detected key genes and important miRNAs involved in the B. superciliaris skin pigmentation process. These results will enhance understanding of molecular mechanisms underlying skin pigmentation and facilitate molecular-assisted selection of highly valued skin colors.

Identifying Hub Genes for Heat Tolerance in Water Buffalo (Bubalus bubalis) Using Transcriptome Data

Heat stress has a detrimental effect on the physiological and production performance of buffaloes. Elucidating the underlying mechanisms of heat stress is challenging, therefore identifying candidate genes is urgent and necessary. We evaluated the response of buffaloes (n = 30) to heat stress using the physiological parameters, ELISA indexes, and hematological parameters. We then performed mRNA and microRNA (miRNA) expression profiles analysis between heat tolerant (HT, n = 4) and non-heat tolerant (NHT, n = 4) buffaloes, as well as the specific modules, significant genes, and miRNAs related to the heat tolerance identified using the weighted gene co-expression network analysis (WGCNA). The results indicated that the buffaloes in HT had a significantly lower rectal temperature (RT) and respiratory rate (RR) and displayed a higher plasma heat shock protein (HSP70 and HSP90) and cortisol (COR) levels than those of NHT buffaloes. Differentially expressed analysis revealed a total of 753 differentially expressed genes (DEGs) and 16 differentially expressed miRNAs (DEmiRNAs) were identified between HT and NHT. Using the WGCNA analysis, these DEGs assigned into 5 modules, 4 of which were significantly correlation with the heat stress indexes. Interestingly, 158 DEGs associated with heat tolerance in the turquoise module were identified, 35 of which were found within the protein-protein interaction network. Several hub genes (IL18RAP, IL6R, CCR1, PPBP, IL1B, and IL1R1) were identified that significantly enriched in the Cytokine-cytokine receptor interaction. The findings may help further elucidate the underlying mechanisms of heat tolerance in buffaloes.

In situ quantification of individual mRNA transcripts in melanocytes discloses gene regulation of relevance to speciation

Functional validation of candidate genes involved in adaptation and speciation remains challenging. Here, we exemplify the utility of a method quantifying individual mRNA transcripts in revealing the molecular basis of divergence in feather pigment synthesis during early-stage speciation in crows. Using a padlock probe assay combined with rolling circle amplification, we quantified cell-type-specific gene expression in the histological context of growing feather follicles. Expression of Tyrosinase Related Protein 1 (TYRP1), Solute Carrier Family 45 member 2 (SLC45A2) and Hematopoietic Prostaglandin D Synthase (HPGDS) was melanocyte-limited and significantly reduced in follicles from hooded crow, explaining the substantially lower eumelanin content in grey versus black feathers. The central upstream Melanocyte Inducing Transcription Factor (MITF) only showed differential expression specific to melanocytes - a feature not captured by bulk RNA-seq. Overall, this study provides insight into the molecular basis of an evolutionary young transition in pigment synthesis, and demonstrates the power of histologically explicit, statistically substantiated single-cell gene expression quantification for functional genetic inference in natural populations.

Transcriptomic Analysis of Coding Genes and Non-Coding RNAs Reveals Complex Regulatory Networks Underlying the Black Back and White Belly Coat Phenotype in Chinese Wuzhishan Pigs

Coat color is one of the most important characteristics for distinguishing Chinese indigenous pig breeds. In Wuzhishan pigs, the animals have black on the back and white on the abdomen. However, the molecular genetic basis of this phenotype is unclear. In this study, we used high-throughput RNA sequencing to compare expression profiles of coding and non-coding RNAs from white and black skin samples obtained from individual Wuzhishan pigs. The expression profiling revealed that 194 lncRNAs (long non-coding RNAs), 189 mRNAs (messenger RNAs), and 162 miRNAs (microRNAs) had significantly different levels of expression (|log₂ fold change| > 1, p-value < 0.05) in white and black skin. Compared to RNA levels in black skin, white skin had higher levels of expression of 185 lncRNAs, 181 mRNAs, and 23 miRNAs and lower levels of expression of 9 lncRNAs, 8 mRNAs, and 139 miRNAs. Functional analysis suggested that the differentially expressed transcripts are involved in biological processes such as melanin biosynthesis, pigmentation and tyrosine metabolism. Several key genes involved in melanogenesis, including MLANA, PMEL, TYR, TYRP1, DTC, TRPM1 and CAMK2A, had significantly different levels of expression in the two skin tissues. Potential lncRNA–miRNA–gene interactions were also examined. A total of 15 lncRNAs, 11 miRNAs and 7 genes formed 23 lncRNA–miRNA–gene pairs, suggesting that complex regulatory networks of coding and non-coding genes underlie the coat color trait in Wuzhishan pigs. Our study provides a foundation for understanding how lncRNA, miRNA and genes interact to regulate coat color in black-back/white-belly pigs. We also constructed lncRNA–miRNA–gene interaction networks to elucidate the complex molecular mechanisms underlying skin physiology and melanogenesis. The results extend our knowledge about the diversity of coat color among different domestic animals and provide a foundation for studying novel mechanisms that control coat color in Chinese indigenous pigs.

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.

Targeted disruption of tyrosinase causes melanin reduction in Carassius auratus cuvieri and its hybrid progeny

The white crucian carp (Carassius auratus cuvieri, WCC) not only is one of the most economically important fish in Asia, characterized by strong reproductive ability and rapid growth rates, but also represents a good germplasm to produce hybrid progenies with heterosis. Gene knockout technique provides a safe and acceptant way for fish breeding. Achieving gene knockout in WCC and its hybrid progeny will be of great importance for both genetic studies and hybridization breeding. Tyrosinase (TYR) is a key enzyme in melanin synthesis. Depletion of tyr in zebrafish and mice results in mosaic pigmentation or total albinism. Here, we successfully used CRISPR-Cas9 to target tyr in WCC and its hybrid progeny (WR) derived from the cross of WCC (♀) and red crucian carp (Carassius auratus red var., RCC, ♂). The level of TYR protein was significantly reduced in mutant WCC. Both the mutant WCC and the mutant WR showed different degrees of melanin reduction compared with the wild-type sibling control fish, resulting from different mutation efficiency ranging from 60% to 90%. In addition, the transcriptional expression profiles of a series of pivotal pigment synthesis genes, i.e. tyrp1, mitfa, mitfb, dct and sox10, were down-regulated in tyr-CRISPR WCC, which ultimately caused a reduction in melanin synthesis. These results demonstrated that tyr plays a key role in melanin synthesis in WCC and WR, and CRISPR-Cas9 is an effective tool for modifying the genome of economical fish. Furthermore, the tyr-CRISPR models could be valuable in understanding fundamental mechanisms of pigment formation in non-model fish.

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.

Transcriptome profile analysis identifies candidate genes for the melanin pigmentation of breast muscle in Muchuan black-boned chicken

Melanin-based coloration in the meat of black-boned chicken is a major economic issue in China. Variation in the pigmentation (hypopigmentation) of chicken muscle causes direct economic losses every year. To determine the molecular mechanisms involved in the melanogenesis of muscle tissue, this study used high-throughput sequencing to compare differences in the transcriptome between black (BM) and white (WM) chicken breast muscles. We constructed 6 cDNA libraries from BM and WM groups in Muchuan black-boned chickens. A comparison between the BM and WM groups revealed 264 differentially expressed genes, of which 152 were upregulated, whereas 112 were downregulated in black muscle. Gene ontology and a Kyoto Encyclopedia of Genes and Genomes pathway analysis identified several differentially enriched biological functions and processes of the 2 muscles. Seven promising candidate genes [PMEL, Ras-related protein RAB29, and 5 solute carrier superfamily genes: SLC6A9, SLC38A4, SLC22A5, SLC35F3, and SLC16A3] may play an important role in the melanogenesis of chicken muscle. Our data provide a valuable resource for identifying genes whose functions are critical for muscle melanogenesis, and will assist studies of the molecular mechanisms of melanogenesis regulation in chicken muscle.

Feather follicles transcriptome profiles in Bashang long-tailed chickens with different plumage colors

Despite the rich variety in plumage color found in nature, genetic studies on how feather follicles affect pigmentation are often limited to animals that have black and white pigment. To test how gene expression influences plumage color, transcriptomes of chicken feather follicles with white, black, hemp, reed catkins, silvery grey, and landscape plumage colors were generated using Illumina sequencing. We generated six RNA-Seq libraries with over 25 million paired-end clean reads per library with percentage of paired-end clean reads ranging from 96.73 to 96.98%. 78% of the reads mapped to the chicken genome, and approximately 70% of the reads were mapped to exons and 6% mapped to introns. Transcriptomes of feather follicles producing hemp and land plumage were similar, but these two showed moderate differences compared with gray and reed colored plumage. The black and white follicle transcriptomes were most divergent from the other colors. We identified several candidate genes, including GPNMB, PMEL, TYRP1, GPR143, OCA2, SOX10, SLC45A2, KRT75, and TYR. All of these genes are known to induce pigment formation in mice. White feathers result from the lack of pigment formation, and our results suggest that the white chickens due to the recessive insertion mutation of TYR. The formation of black area size and color depth may be due to the expression levels of GPNMB, PMEL, TYRP1, GPR143, OCA2, SOX10, SLC45A2, KRT75, and TYR. The GO analysis of the differentially expressed genes (DEGs) revealed that DEGs in our transcriptome analysis were enriched in cytoskeleton and cell structure related pathways. The black plumage transcriptome showed significant differences in melanogenesis, tyrosine metabolism, and riboflavin metabolism compared with transcriptomes of other plumage colors. The transcriptome profiles of the different chicken plumage colors provide a valuable resource to understand how gene expression influences plumage color, and will be an important resource for identifying candidate genes in breeding programs.

Comparative liver transcriptome analysis in ducklings infected with duck hepatitis A virus 3 (DHAV-3) at 12 and 48 hours post-infection through RNA-seq

Duck hepatitis A virus 3 (DHAV-3), the only member of the novel genus Avihepatovirus, in the family Picornaviridae, can cause significant economic losses for duck farms in China. Reports on the pathogenicity and the antiviral molecular mechanisms of the lethal DHAV-3 strain in ducklings are inadequate and remain poorly understood. We conducted global gene expression profiling and screened differentially expressed genes (DEG) of duckling liver tissues infected with lethal DHAV-3. There were 1643 DEG and 8979 DEG when compared with mock ducklings at 12 hours post-infection (hpi) and at 48 hpi, respectively. Gene pathway analysis of DEG highlighted mainly biological processes involved in metabolic pathways, host immune responses, and viral invasion. The results may provide valuable information for us to explore the pathogenicity of the virulent DHAV-3 strain and to improve our understanding of host–virus interactions.

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

OBJECTIVE

The aim of this study was to investigate the effect of single nucleotide polymorphisms (SNPs) of the melanogenesis associated transcription factor (MITF) and dopachrome tautomerase (DCT) genes on plumage coloration in Asian native duck breeds. MITF encodes a protein for microphthalmia-associated transcription factor, which regulates the development and function of melanocytes for pigmentation of skin, hair, and eyes. Among the tyrosinase-related family genes, DCT is a pigment cell-specific gene that plays important roles in the melanin synthesis pathway and the expression of skin, feather, and retina color.

METHODS

Five Asian duck varieties (black Korean native, white Korean native, commercial Peking, Nageswari, and Bangladeshi Deshi white ducks) were investigated to examine the polymorphisms associated with plumage colors. Among previously identified SNPs, three synonymous SNPs and one indel of MITF and nine SNPs in exon regions of DCT were genotyped. The allele frequencies for SNPs of the black and white plumage color populations were estimated and Fisher's exact test was conducted to assess the association between the allele frequencies of these two populations.

RESULTS

Two synonymous SNPs (c.114T>G and c.147T>C) and a 14-bp indel (GCTGCAAAC AGATG) in intron 7 of MITF were significantly associated with the black- and white-colored breeds (p<0.001). One non-synonymous SNP [c.938A>G (p.His313Arg)] in DCT, was highly significantly associated (p<0.001) and a synonymous SNP (c.753A>G) was significantly associated (p<0.05) with black and white color plumage in the studied duck populations.

CONCLUSION

The results of this study provide a basis for further investigations of the associations between polymorphisms and plumage color phenotypes in Asian duck breeds.

Comparative Transcriptome Analysis of Mink (Neovison vison) Skin Reveals the Key Genes Involved in the Melanogenesis of Black and White Coat Colour

Farmed mink (Neovison vison) is one of the most important fur-bearing species worldwide, and coat colour is a crucial qualitative characteristic that contributes to the economic value of the fur. To identify additional genes that may play important roles in coat colour regulation, Illumina/Solexa high-throughput sequencing technology was used to catalogue the global gene expression profiles in mink skin with two different coat colours (black and white). RNA-seq analysis indicated that a total of 12,557 genes were differentially expressed in black versus white minks, with 3,530 genes up-regulated and 9,027 genes down-regulated in black minks. Significant differences were not observed in the expression of MC1R and TYR between the two different coat colours, and the expression of ASIP was not detected in the mink skin of either coat colour. The expression levels of KITLG, LEF1, DCT, TYRP1, PMEL, Myo5a, Rab27a and SLC7A11 were validated by qRT-PCR, and the results were consistent with RNA-seq analysis. This study provides several candidate genes that may be associated with the development of two coat colours in mink skin. These results will expand our understanding of the complex molecular mechanisms underlying skin physiology and melanogenesis in mink and will provide a foundation for future studies.

Comparative Transcriptome Analysis Identifies Candidate Genes Related to Skin Color Differentiation in Red Tilapia

Red tilapia is becoming more popular for aquaculture production in China in recent years. However, the pigmentation differentiation in genetic breeding is the main problem limiting its development of commercial red tilapia culture and the genetic basis of skin color variation is still unknown. In this study, we conducted Illumina sequencing of transcriptome on three color variety red tilapia. A total of 224,895,758 reads were generated, resulting in 160,762 assembled contigs that were used as reference contigs. The contigs of red tilapia transcriptome had hits in the range of 53.4% to 86.7% of the unique proteins of zebrafish, fugu, medaka, three-spined stickleback and tilapia. And 44,723 contigs containing 77,423 simple sequence repeats (SSRs) were identified, with 16,646 contigs containing more than one SSR. Three skin transcriptomes were compared pairwise and the results revealed that there were 148 common significantly differentially expressed unigenes and several key genes related to pigment synthesis, i.e. tyr, tyrp1, silv, sox10, slc24a5, cbs and slc7a11, were included. The results will facilitate understanding the molecular mechanisms of skin pigmentation differentiation in red tilapia and accelerate the molecular selection of the specific strain with consistent skin colors.

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.

Transcriptomic analyses of regenerating adult feathers in chicken

Background

Feathers have diverse forms with hierarchical branching patterns and are an excellent model for studying the development and evolution of morphological traits. The complex structure of feathers allows for various types of morphological changes to occur. The genetic basis of the structural differences between different parts of a feather and between different types of feather is a fundamental question in the study of feather diversity, yet there is only limited relevant information for gene expression during feather development.

Results

We conducted transcriptomic analysis of five zones of feather morphologies from two feather types at different times during their regeneration after plucking. The expression profiles of genes associated with the development of feather structure were examined. We compared the gene expression patterns in different types of feathers and different portions of a feather and identified morphotype-specific gene expression patterns. Many candidate genes were identified for growth control, morphogenesis, or the differentiation of specific structures of different feather types.

Conclusion

This study laid the ground work for studying the evolutionary origin and diversification of feathers as abundant data were produced for the study of feather morphogenesis. It significantly increased our understanding of the complex molecular and cellular events in feather development processes and provided a foundation for future studies on the development of other skin appendages.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1966-6) contains supplementary material, which is available to authorized users.

Identification and profiling of microRNAs associated with white and black plumage pigmentation in the white and black feather bulbs of ducks by RNA sequencing

MicroRNAs (miRNAs) play important roles in many biological processes by regulating gene expression at the post-transcriptional level. However, the mechanism by which specific miRNAs may regulate plumage pigmentation has remained largely elusive. In this study, we sequenced miRNAs using Solexa sequencing and then performed a detailed analysis of their expression profiles between the black and white feather bulbs of ducks from Cui Hei, Kaiya, Liancheng pure breeds and a Kaiya-Liancheng F2 population. mirdeep2 software identified 121 conserved and eight novel miRNAs. Five differentially expressed miRNAs between the two tissues types were also identified by degseq software. Notably, miR-204 was predominantly expressed in black feather bulbs. To further validate the sequencing data, we applied stem-loop quantitative PCR of ten known miRNAs based on the identified sequences. Furthermore, in exploring the temporal expression pattern of miR-204, we performed profiling in nine duck tissues. The targets of these miRNAs were predicted using a PITA algorithm and were later grouped based on Gene Ontology and KEGG pathway analysis using the DAVID website. The melanogenesis pathway was among the identified signalling pathways, implying key roles of these miRNAs in plumage pigmentation. Expression analysis of the target genes in the melanogenesis pathways was also performed. This study provides the foundation for subsequent studies on the prospective practical role for such miRNAs in post-transcriptional gene regulation linked to plumage pigmentation.

Transcriptome analysis of adiposity in domestic ducks by transcriptomic comparison with their wild counterparts

Excessive adiposity is a major problem in the duck industry, but its molecular mechanisms remain unknown. Genetic comparisons between domestic and wild animals have contributed to the exploration of genetic mechanisms responsible for many phenotypic traits. Significant differences in body fat mass have been detected between domestic and wild ducks. In this study, we used the Peking duck and Anas platyrhynchos as the domestic breed and wild counterpart respectively and performed a transcriptomic comparison of abdominal fat between the two breeds to comprehensively analyze the transcriptome basis of adiposity in ducks. We obtained approximately 350 million clean reads; assembled 61 250 transcripts, including 23 699 novel ones; and identified alternative 5' splice sites, alternative 3' splice sites, skipped exons and retained intron as the main alternative splicing events. A differential expression analysis between the two breeds showed that 753 genes exhibited differential expression. In Peking ducks, some lipid metabolism-related genes (IGF2, FABP5, BMP7, etc.) and oncogenes (RRM2, AURKA, CYR61, etc.) were upregulated, whereas genes related to tumor suppression and immunity (TNFRSF19, TNFAIP6, IGSF21, NCF1, etc.) were downregulated, suggesting adiposity might closely associate with tumorigenesis in ducks. Furthermore, 280 576 single-nucleotide variations were found differentiated between the two breeds, including 8641 non-synonymous ones, and some of the non-synonymous ones were found enriched in genes involved in lipid-associated and immune-associated pathways, suggesting abdominal fat of the duck undertakes both a metabolic function and immune-related function. These datasets enlarge our genetic information of ducks and provide valuable resources for analyzing mechanisms underlying adiposity in ducks.

A Comparative Transcriptome Analysis between Wild and Albino Yellow Catfish (Pelteobagrus fulvidraco)

Body colours are important and striking features for individual survival and reproductive success, in particular in vertebrates where mating behaviour and mate preference may be strongly influenced by non-normal phenotypes. Pigmentation disorders may be generated by disruption of one or many independent genes as well as by environmental factors. The first discovery of albino yellow catfish (Pelteobagrus fulvidraco Richardson) with golden skin colour from fish farms in China provides us valuable material to study the molecular mechanism underlying the abnormalities of pigmentation. In this study, transcriptome sequencing of fin tissues corresponding to the distinct body colours, wild type and mutant albino yellow catfish, were performed using Illumina sequencing technology. Based on next-generation sequencing technology and de novo assembly, we generated a transcriptome of P. fulvidraco. A number of genes differentially expressed between the wild types and albinos were identified, suggesting their contribution to the different phenotypes and fitness. However, non-synonymous mutations result from single nucleotide substitutions residing in coding regions may not contribute to such differences. Based on the high-throughput expression data generated for the two different types of P. fulvidraco, we found that alterations of expression pattern may be more common than non-synonymous mutations. The transcriptome of P. fulvidraco will be an invaluable resource for subsequent comparative genomics and evolutionary analyses of this economically important fish.

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.

Feathers as a source of RNA for genomic studies in avian species

Dioxins and dioxin-like chemicals (DLCs) cause a suite of adverse effects in terrestrial species. Most of the adverse effects occur subsequent to binding to the aryl hydrocarbon receptor. Avian species vary in their sensitivity to the effects of DLCs and current research indicates that this is mediated by variations in the amino acid sequence within the ligand binding domain (LBD) of the aryl hydrocarbon receptor 1 (AHR1). Eighty-eight avian species have been classified into three broad categories of sensitivity, based on the amino acid variations within the AHR1 LBD: sensitive type 1 (Ile324_Ser380), moderately sensitive type 2 (Ile324_Ala380), and relatively insensitive type 3 (Val324_Ala380). Risk assessment of avian species can be complicated due to the variability in sensitivity among species. A predictive tool for selecting the priority species at a given site would have broad implications for the risk assessment community. We present a method for AHR1 genotyping using plucked feathers as a source of RNA. The method is extremely robust, requires minimal sample processing and handling, and eliminates the need for blood sampling or tissue collection from the species of interest. Using this method we were able to determine the amino acid sequence of the AHR LBD of three avian species: the chicken, the herring gull, and the zebra finch, and to categorize them based on the identity of amino acids at key sites within the LBD.