Skin transcriptome profiles associated with coat color in sheep

Wnt signaling pathway and retinoic acid signaling pathway involved in delamination and migration of chicken trunk NCCs and contributing to HVP phenotype

Hyperpigmentation of the visceral peritoneum (HVP) is a hereditary trait that significantly affects the carcass quality in bearded chickens, yet its molecular mechanisms remain unclear. This study utilized data-independent acquisition proteomics to analyze the protein expression profiles of black peritoneum (B), faded peritoneum (F), and normal peritoneum (N) in bearded chickens at 40 and 120 d of age. Combined with histopathological and functional enrichment analyses, we revealed the regulatory network underlying HVP formation. Results indicated that the melanin content was significantly elevated in HVP samples, without accompanying inflammatory responses or tumor characteristics, suggesting that its formation is driven by developmental abnormalities. A total of 9,375 high-confidence proteins were identified through proteomics, with differentially abundant proteins at 40 d of age (219 proteins) primarily enriched in ribosomal function, tyrosine metabolism, and melanin synthesis pathways. In comparison, at 120 d of age (246 proteins), they were enriched in transcription regulation and chromatin remodeling pathways. The abnormal expression of key co-expressed proteins DHRS3 and DACT1 suggests that the dysregulation of retinoic acid (RA) and the Wnt signaling pathway may promote the directed differentiation of melanocytes by regulating neural crest cells (NCCs). The reduced abundance of the chondroitin sulfate proteoglycan, VCAN, weakened the peritoneal barrier function, whereas estradiol accelerated melanin synthesis via hormonal microenvironmental regulation. Furthermore, the formation of HVP led to a reprogramming of energy metabolism, reduced fat deposition, and a downregulation of immune-related molecules, implying that pigment deposition may weaken the chicken immune response. This study systematically elucidates the molecular mechanisms of HVP and provides potential targets for molecular breeding of HVP.

Genetic advancements and future directions in ruminant livestock breeding: from reference genomes to multiomics innovations

Ruminant livestock provide a rich source of products, such as meat, milk, and wool, and play a critical role in global food security and nutrition. Over the past few decades, genomic studies of ruminant livestock have provided valuable insights into their domestication and the genetic basis of economically important traits, facilitating the breeding of elite varieties. In this review, we summarize the main advancements for domestic ruminants in reference genome assemblies, population genomics, and the identification of functional genes or variants for phenotypic traits. These traits include meat and carcass quality, reproduction, milk production, feed efficiency, wool and cashmere yield, horn development, tail type, coat color, environmental adaptation, and disease resistance. Functional genomic research is entering a new era with the advancements of graphical pangenomics and telomere-to-telomere (T2T) gap-free genome assembly. These advancements promise to improve our understanding of domestication and the molecular mechanisms underlying economically important traits in ruminant livestock. Finally, we provide new perspectives and future directions for genomic research on ruminant genomes. We suggest how ever-increasing multiomics datasets will facilitate future studies and molecular breeding in livestock, including the potential to uncover novel genetic mechanisms underlying phenotypic traits, to enable more accurate genomic prediction models, and to accelerate genetic improvement programs.

Transcriptome-based screening and validation of key genes for wool color in cashmere goats

BACKGROUND

Colored wool from cashmere goats is increasingly popular among consumers, but the transcriptomic differences between coat colors are poorly understood.

OBJECTIVES

This study aimed to screen for coat color regulation-associated genes in cashmere goats to ascertain their underlying molecular mechanisms.

METHODS

Transcriptomic sequencing of skin tissues from black (BC), brown (YC), and white cashmere (WC) goats was performed. Immunohistochemistry and western blotting were used to validate SLC24A4 and DCT expression, two essential genes identified for coat color determination.

RESULTS

We identified 6,518 differentially expressed genes (DEGs) in the BC vs. WC group (3,919 upregulated, 2,599 downregulated). Next, 5,593 DEGs were identified in the YC vs. WC group (3,629 upregulated, 1,964 downregulated). Finally, 4,538 DEGs were expressed in both groups, with 1,980 and 1,055 DEGs exclusively expressed in either group. Functions and pathways associated with hair color were enriched, including melanosomes, melanocyte migration, melanin biosynthesis processes and functions, and melanogenesis pathways. TYRP1, SLC24A4, PMEL, OCA2, and DCT were significantly upregulated in BC goat skin, while ASIP was significantly upregulated in YC skin. Additionally, KIT, POMC, SLC24A5, Wnt3a, and EDN3 were DEGs for different coat colors. Immunohistochemistry revealed SLC24A4 and DCT expression in dermal papillae, inner and outer root sheaths, and the hair follicle matrix. Western blotting showed that SLC24A4 protein levels were highest in BC goat skin. DCT protein levels were also highest in BC goat skin, albeit not significantly.

CONCLUSION

These results further our understanding of coat color regulation in cashmere goats, establishing a foundation for their molecular breeding.

Skin proteomic screening and functional analysis of differential proteins associated with coat color in sheep (Ovis aries)

OBJECTIVE

Coat color is an important characteristic and economic trait in domestic sheep. In this study, we explored the potential mechanisms and the signaling pathways involved in coat color regulation for sheep.

METHODS

Isobaric tags for relative and absolute quantification (iTRAQ) technology was used to catalog global protein expression profiles in skin of sheep with black versus white coat color. Immunofluorescence was used to observe the expression localization of differential protein. Western blot and quantitative real time polymerase chain reaction (qRT-PCR) were used to evaluate their role in the coat color formation of sheep.

RESULTS

A total of 136 differential proteins were obtained in different coat colors, including 101 up-regulated and 35 down-regulated. Pigmentation function entries were enriched through gene ontology annotation. Tyrosine metabolism and platelet activation signaling pathway were extracted by Kyoto encyclopedia of genes and genomes analysis. Apolipoprotein A-1 (APOA1) and fibrinogen alpha chain (FGA) were found to be critical differential proteins by the interaction of differential proteins in the direct-interaction network diagram. Strikingly, twenty candidate differential proteins were screened, from which beta-actin (ACTB) protein showed higher expression in white sheep skin, while albumin (ALB), APOA1, MAOA (amine oxidase) and FGA proteins showed higher expression in black sheep skin, which was validated by immunofluorescence, western blot, and qRT-PCR.

CONCLUSION

This study identified several novel proteins that may be involved in the coat color formation of sheep. The white and black sheep skin proteome profiles obtained provide a valuable resource for future research to understand the network of protein expression controlling skin physiology and melanogenesis in sheep.

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep

BACKGROUND

Managing genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright's fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs).

RESULTS

Katahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms.

CONCLUSIONS

In this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Transcriptomic and Metabolomic Analyses Reveal Molecular Regulatory Networks for Pigmentation Deposition in Sheep

Domestic animals have multiple phenotypes of skin and coat color, which arise from different genes and their products, such as proteins and metabolites responsible with melanin deposition. However, the complex regulatory network of melanin synthesis remains to be fully unraveled. Here, the skin and tongue tissues of Liangshan black sheep (black group) and Liangshan semi-fine-wool sheep (pink group) were collected, stained with hematoxylin-eosin (HE) and Masson-Fontana, and the transcriptomic and metabolomic data were further analyzed. We found a large deposit of melanin granules in the epidermis of the black skin and tongue. Transcriptome and metabolome analysis identified 744 differentially expressed genes (DEGs) and 443 differentially expressed metabolites (DEMs) between the pink and black groups. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses revealed the DEGs and DEMs were mainly enriched in the pathways of secondary metabolic processes, melanin biosynthesis processes, melanin metabolism processes, melanosome membranes, pigment granule membranes, melanosome, tyrosine metabolism, and melanogenesis. Notably, we revealed the gene ENSARG00020006042 may be a family member of YWHAs and involved in regulating melanin deposition. Furthermore, several essential genes (TYR, TYRP1, DCT, PMEL, MLANA, SLC45A2) were significantly associated with metabolite prostaglandins and compounds involved in sheep pigmentation. These findings provide new evidence of the strong correlation between prostaglandins and related compounds and key genes that regulate sheep melanin synthesis, furthering our understanding of the regulatory mechanisms and molecular breeding of pigmentation in sheep.

Selection signatures of wool color in Gangba sheep revealed by genome-wide SNP discovery

BACKGROUND

Gangba sheep as a famous breed of Tibetan sheep, its wool color is mainly white and black. Gangba wool is economically important as a high-quality raw material for Tibetan blankets and Tibetan serge. However, relatively few studies have been conducted on the wool color of Tibetan sheep.

RESULTS

To fill this research gap, this study conducted an in-depth analysis of two populations of Gangba sheep (black and white wool color) using whole genome resequencing to identify genetic variation associated with wool color. Utilizing PCA, Genetic Admixture, and N-J Tree analyses, the present study revealed a consistent genetic relationship and structure between black and white wool colored Gangba sheep populations, which is consistent with their breed history. Analysis of selection signatures using multiple methods (FST, π ratio, Tajima's D), 370 candidate genes were screened in the black wool group (GBB vs GBW); among them, MC1R, MLPH, SPIRE2, RAB17, SMARCA4, IRF4, CAV1, USP7, TP53, MYO6, MITF, MC2R, TET2, NF1, JAK1, GABRR1 genes are mainly associated with melanin synthesis, melanin delivery, and distribution. The enrichment results of the candidate genes identified 35 GO entries and 19 KEGG pathways associated with the formation of the black phenotype. 311 candidate genes were screened in the white wool group (GBW vs GBB); among them, REST, POU2F1, ADCY10, CCNB1, EP300, BRD4, GLI3, and SDHA genes were mainly associated with interfering with the differentiation of neural crest cells into melanocytes, affecting the proliferation of melanocytes, and inhibiting melanin synthesis. 31 GO entries and 22 KEGG pathways were associated with the formation of the white phenotype.

CONCLUSIONS

This study provides important information for understanding the genetic mechanism of wool color in Gangba, and provides genetic knowledge for improving and optimizing the wool color of Tibetan sheep. Genetic improvement and selective breeding to produce wool of specific colors can meet the demand for a diversity of wool products in the Tibetan wool textile market.

Uncovering the genetic diversity and adaptability of Butuo Black Sheep through whole-genome re-sequencing

The Butuo Black Sheep (BBS) is well-known for its ability to thrive at high altitudes, resist diseases, and produce premium-quality meat. Nonetheless, there is insufficient data regarding its genetic diversity and population-specific Single nucleotide polymorphisms (SNPs). This paper centers on the genetic diversity of (BBS). The investigation conducted a whole-genome resequencing of 33 BBS individuals to recognize distinct SNPs exclusive to BBS. The inquiry utilized bioinformatic analysis to identify and explain SNPs and pinpoint crucial mutation sites. The findings reveal that reproductive-related genes (GHR, FSHR, PGR, BMPR1B, FST, ESR1), lipid-related genes (PPARGC1A, STAT6, DGAT1, ACACA, LPL), and protein-related genes (CSN2, LALBA, CSN1S1, CSN1S2) were identified as hub genes. Functional enrichment analysis showed that genes associated with reproduction, immunity, inflammation, hypoxia, PI3K-Akt, and AMPK signaling pathways were present. This research suggests that the unique ability of BBS to adapt to low oxygen levels in the plateau environment may be owing to mutations in a variety of genes. This study provides valuable insights into the genetic makeup of BBS and its potential implications for breeding and conservation efforts. The genes and SPNs identified in this study could serve as molecular markers for BBS.

m6A mRNA Methylation Analysis Provides Novel Insights into Pigmentation in Sheep Skin

N6-methyladenosine (m⁶A) is the most universal post-transcriptional modification of mRNA which may play important roles in verious species. However, the potential roles of m⁶A in the pigmentation of skin are not completely understood. To explore the role of m6A modification in pigmentation of sheep skin, we used MeRIP-seq and RNA-seq to profile the skin transcriptome in black and white coat color (n=3). Our results showed that an average of 7701 m⁶A peaks were obtained for all samples and the average length was 305.89 bp. The GGACUU sequence was the most enrichment motif and shared in black skin and white skin. The m⁶A peaks were mainly enriched in the CDS, 3'UTR and 5'UTR, especially in CDS region near the stop codon of the transcript. 235 significantly differential peaks were found in black skin vs. white skin. The KEGG signaling pathways of downregulated and upregulated m⁶A peaks were mainly enriched in AGE-RAGE signaling pathway in diabetic complications, Viral carcinogenesis, Transcriptional misregulation in cancer, ABC transporters, Basal transcription factors and Thyroid hormone synthesis (P value <0.05). For RNA-seq, 71 differently expressed genes (DEGs) were scanned in black skin vs. white skin. DEGs were significantly enriched in tyrosine metabolism, melanogenesis, neuroactive ligand-receptor interaction pathway (P value <0.05). Combined m⁶A-seq and RNA-seq analysis showed that the hyper-up genes and hypo-up genes were both enriched in ErbB signaling pathway (P value <0.05). In conclusion, it provide a basis for further research into the functions of m⁶A methylation modifications in pigmentation.

Whole Genome Resequencing Reveals Selection Signals Related to Wool Color in Sheep

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

Mitochondrial calcium signaling mediated transcriptional regulation of keratin filaments is a critical determinant of melanogenesis

Mitochondria are versatile organelles that regulate several physiological functions. Many mitochondria-controlled processes are driven by mitochondrial Ca2+ signaling. However, role of mitochondrial Ca2+ signaling in melanosome biology remains unknown. Here, we show that pigmentation requires mitochondrial Ca2+ uptake. In vitro gain and loss of function studies demonstrated that Mitochondrial Ca2+ Uniporter (MCU) is crucial for melanogenesis while the MCU rheostats, MCUb and MICU1 negatively control melanogenesis. Zebrafish and mouse models showed that MCU plays a vital role in pigmentation in vivo. Mechanistically, MCU controls activation of transcription factor NFAT2 to induce expression of three keratins (keratin 5, 7 and 8), which we report as positive regulators of melanogenesis. Interestingly, keratin 5 in turn modulates mitochondrial Ca2+ uptake thereby this signaling module acts as a negative feedback loop that fine-tunes both mitochondrial Ca2+ signaling and melanogenesis. Mitoxantrone, an FDA approved drug that inhibits MCU, decreases physiological melanogenesis. Collectively, our data demonstrates a critical role for mitochondrial Ca2+ signaling in vertebrate pigmentation and reveal the therapeutic potential of targeting MCU for clinical management of pigmentary disorders. Given the centrality of mitochondrial Ca2+ signaling and keratin filaments in cellular physiology, this feedback loop may be functional in a variety of other pathophysiological conditions.

Transcriptome Profile Analysis Identifies Candidate Genes for the Melanin Pigmentation of Skin in Tengchong Snow Chickens

Tengchong Snow chickens are one of the most precious, black-boned chickens in Yunnan province and usually produce black meat. However, we found a small number of white meat traits in the chicken population during feeding. In order to determine the pattern of melanin deposition and the molecular mechanism of formation in the Tengchong Snow chicken, we measured the luminance value (L value) and melanin content in the skin of black meat chickens (Bc) and white meat chickens (Wc) using a color colorimeter, ELISA kit, and enzyme marker. The results showed that the L value of skin tissues in black meat chickens was significantly lower than that of white meat chickens, and the L value of skin tissues gradually increased with an increase in age. The melanin content of skin tissues in black meat chickens was higher than that of white meat chickens, and melanin content in the skin tissues gradually decreased with an increase in age, but this difference was not significant (p > 0.05); the L value of skin tissues in black meat chickens was negatively correlated with melanin content, and the correlation coefficient was mostly above -0.6. In addition, based on the phenotypic results, we chose to perform the comparative transcriptome profiling of skin tissues at 90 days of age. We screened a total of 44 differential genes, of which 32 were upregulated and 12 were downregulated. These DEGs were mainly involved in melanogenesis, tyrosine metabolism and RNA transport. We identified TYR, DCT, and EDNRB2 as possible master effector genes for skin pigmentation in Tengchong Snow black meat chickens through DEGs analysis. Finally, we measured the mRNA of TYR, DCT, MC1R, EDNRB2, GPR143, MITF, and TYRP1 genes through a quantitative real-time polymerase chain reaction (qPCR) and found that the mRNA of all the above seven genes decreased with increasing age. In conclusion, our study initially constructed an evaluation system for the black-boned traits of Tengchong Snow chickens and found key candidate genes regulating melanin deposition, which could provide an important theoretical basis for the selection and breeding of black-boned chickens.

Analysis of reproduction-related transcriptomes on pineal-hypothalamic-pituitary-ovarian tissues during estrus and anestrus in Tan sheep

Seasonal estrus is an important factor limiting the fertility of some animals such as sheep. Promoting estrus in the anestrus season is one of the major ways in improving the fecundity of seasonally breeding animals. The pineal-hypothalamus-pituitary-ovary (PHPO) axis plays a decisive role in regulating animal reproduction. However, the molecular mechanisms by which the PHPO axis regulates seasonal reproduction in animals are not well understood, especially in Tan sheep. To this end, we collected pineal, hypothalamus, pituitary and ovary tissues from Tan sheep during estrus and anestrus for RNA-Sequencing, and performed bioinformatics analysis on the entire regulatory axis of the pineal-hypothalamic-pituitary-ovary (PHPO). The results showed that 940, 1,638, 750, and 971 DEGs (differentially expressed genes, DEGs) were identified in pineal, hypothalamus, pituitary and ovary, respectively. GO analysis showed that DEGs from PHPO axis-related tissues were mainly enriched in "biological processes" such as transmembrane transport, peptide and amide biosynthesis and DNA synthesis. Meanwhile, KEGG enrichment analysis showed that the bile acid secretion pathway and the neuroactive ligand-receptor interaction pathway were significantly enriched. Additionally, four potential candidate genes related to seasonal reproduction (VEGFA, CDC20, ASPM, and PLCG2) were identified by gene expression profiling and protein-protein interaction (PPI) analysis. These findings will contribute to be better understanding of seasonal reproduction regulation in Tan sheep and will serve as a useful reference for molecular breeding of high fertility Tan sheep.

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.

Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes

BACKGROUND

After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.

REVIEW

Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.

CONCLUSIONS

Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.

De novo assembly transcriptome analysis reveals the genes associated with body color formation in the freshwater ornamental shrimps Neocaridina denticulate sinensis

The body color of Neocaridina denticulate sinensis is a compelling phenotypic trait, in which a cascade of carotenoid metabolic processes plays an important role. The study was conducted to compare the transcriptome of cephalothoraxes among three pigmentation phenotypes (red, blue, and chocolate) of N. denticulate sinensis. The purpose of this study was to explore the candidate genes associated with different colors of N. denticulate sinensis. Nine cDNA libraries in three groups were constructed from the cephalothoraxes of shrimps. After assembly, 75022 unigenes were obtained in total with an average length of 1026 bp and N50 length of 1876 bp. There were 45977, 25284, 23605, 21913 unigenes annotated in the Nr, Swissprot, KOG, and KEGG databases, respectively. Differential expression analysis revealed that there were 829, 554, and 3194 differentially expressed genes (DEGs) in RD vs BL, RD vs CH, and BL vs CH, respectively. These DEGs may play roles in the absorption, transport, and metabolism of carotenoids. We also emphasized that electron transfer across the inner mitochondrial membrane (IMM) was a key process in pigment metabolism. In addition, a total of 6328 simple sequence repeats (SSRs) were also detected in N. denticulate sinensis. The results laid a solid foundation for further research on the molecular mechanism of integument pigmentation in the crustacean and contributed to developing more attractive aquatic animals.

Evolutionary selection of alleles in the melanophilin gene that impacts on prostate organ function and cancer risk

Background and objectives

Several hundred inherited genetic variants or SNPs that alter the risk of cancer have been identified through genome-wide association studies. In populations of European ancestry, these variants are mostly present at relatively high frequencies. To gain insight into evolutionary origins, we screened a series of genes and SNPs linked to breast or prostate cancer for signatures of historical positive selection.

Methodology

We took advantage of the availability of the 1000 genome data and we performed genomic scans for positive selection in five different Caucasian populations as well as one African reference population. We then used prostate organoid cultures to provide a possible functional explanation for the interplay between the action of evolutionary forces and the disease risk association.

Results

Variants in only one gene showed genomic signatures of positive, evolutionary selection within Caucasian populations melanophilin (MLPH). Functional depletion of MLPH in prostate organoids, by CRISPR/Cas9 mutation, impacted lineage commitment of progenitor cells promoting luminal versus basal cell differentiation and on resistance to androgen deprivation.

Conclusions and implications

The MLPH variants influencing prostate cancer risk may have been historically selected for their adaptive benefit on skin pigmentation but MLPH is highly expressed in the prostate and the derivative, positively selected, alleles decrease the risk of prostate cancer. Our study suggests a potential functional mechanism via which MLPH and its genetic variants could influence risk of prostate cancer, as a serendipitous consequence of prior evolutionary benefits to another tissue.

Lay Summary

We screened a limited series of genomic variants associated with breast and prostate cancer risk for signatures of historical positive selection. Variants within the melanophilin (MLPH) gene fell into this category. Depletion of MLPH in prostate organoid cultures, suggested a potential functional mechanism for impacting on cancer risk, as a serendipitous consequence of prior evolutionary benefits to another tissue.

In Silico Prediction of Transcription Factor Collaborations Underlying Phenotypic Sexual Dimorphism in Zebrafish (Danio rerio)

The transcriptional regulation of gene expression in higher organisms is essential for different cellular and biological processes. These processes are controlled by transcription factors and their combinatorial interplay, which are crucial for complex genetic programs and transcriptional machinery. The regulation of sex-biased gene expression plays a major role in phenotypic sexual dimorphism in many species, causing dimorphic gene expression patterns between two different sexes. The role of transcription factor (TF) in gene regulatory mechanisms so far has not been studied for sex determination and sex-associated colour patterning in zebrafish with respect to phenotypic sexual dimorphism. To address this open biological issue, we applied bioinformatics approaches for identifying the predicted TF pairs based on their binding sites for sex and colour genes in zebrafish. In this study, we identified 25 (e.g., STAT6-GATA4; JUN-GATA4; SOX9-JUN) and 14 (e.g., IRF-STAT6; SOX9-JUN; STAT6-GATA4) potentially cooperating TFs based on their binding patterns in promoter regions for sex determination and colour pattern genes in zebrafish, respectively. The comparison between identified TFs for sex and colour genes revealed several predicted TF pairs (e.g., STAT6-GATA4; JUN-SOX9) are common for both phenotypes, which may play a pivotal role in phenotypic sexual dimorphism in zebrafish.

Comparative transcriptome and histological analyses provide insights into the skin pigmentation in Minxian black fur sheep (Ovis aries)

Background

Minxian black fur (MBF) sheep are found in the northwestern parts of China. These sheep have developed several special traits. Skin color is a phenotype subject to strong natural selection and diverse skin colors are likely a consequence of differences in gene regulation.

Methods

Skin structure, color differences, and gene expression (determined by RNA sequencing) were evaluated the Minxian black fur and Small-tail Han sheep (n = 3 each group), which are both native Chinese sheep breeds.

Results

Small-tail Han sheep have a thicker skin and dermis than the Minxian black fur sheep (P < 0.01); however, the quantity of melanin granules is greater (P < 0.01) in Minxian black fur sheep with a more extensive distribution in skin tissue and hair follicles. One hundred thirty-three differentially expressed genes were significantly associated with 37 ontological terms and two critical KEGG pathways for pigmentation (“tyrosine metabolism” and “melanogenesis” pathways). Important genes from those pathways with known involvement in pigmentation included OCA2 melanosomal transmembrane protein (OCA2), dopachrome tautomerase (DCT), tyrosinase (TYR) and tyrosinase related protein (TYRP1), melanocortin 1 receptor (MC1R), and premelanosome protein (PMEL). The results from our histological and transcriptome analyses will form a foundation for additional investigation into the genetic basis and regulation of pigmentation in these sheep breeds.

Preliminary comparison of skin transcriptome from sheep with different wool fibre diameters

Context Wool is one of the most important animal fibres for the textile industry, and its diameter directly affects its economic value. However, the molecular mechanisms underlying wool fibre diameter (FD) in sheep have not been fully elucidated. Aims The aims of the work were to make an initial comparison of skin transcriptomes from sheep with wool of high and low FD, and to identify key genes affecting FD. Methods High-throughput RNA-Seq technology was employed to explore the skin transcriptome, using three sheep with fine wool (FD &lt;21.0 μm) and three sheep with coarse wool (FD &gt;27.0 μm). Key results We obtained 28607228 bp of clean sequence data, 78.9% (±3.8%) of which uniquely aligned to the reference genome across the six samples. In total, 19914 mRNA transcripts were expressed (FPKM &gt;0) in the six skin samples, among which were certain well-known genes involved in the skin–hair cycle, such as KRTAP7-1, KRT14, Wnt10b, Wnt2b, β-catenin and FGF5. Furthermore, 467 genes were significantly differentially expressed between the fine-wool and coarse-wool groups, including 21 genes with upregulated and 446 genes with downregulated expression in the sheep with lower FD. These differentially expressed genes were particularly enriched in the gene ontology processes related to lipid metabolism, skin development, differentiation and immune function (P &lt; 0.05). The biological processes were involved in collagen catabolism, negative regulation of macromolecule metabolism, steroid hormone stimulation and lipid metabolism. A significant Kyoto Encyclopedia of Genes and Genomes pathway involving the metabolism of lipids and lipoproteins was also enriched, revealing that lipid metabolism might be one of the key factors affecting FD. The expression of these differentially expressed genes that were involved in the metabolism of lipids and lipoproteins pathway was verified by quantitative real-time PCR (qPCR). The correlation between the mRNA expression level from qPCR and RNA-Seq data was 0.999 (P &lt; 0.001). Conclusions The 467 differentially expressed genes, especially those involved in lipid metabolism and immune function, may play key roles in wool follicle metabolism and the expression of wool FD. Implications This study provided valuable data for future studies aimed at elucidating the mechanisms that underlie wool follicle metabolism and wool FD. The work may also have implications for studies of the human hair follicle.

Genetics and Evolution of Mammalian Coat Pigmentation

The diversity of mammalian coat colors, and their potential adaptive significance, have long fascinated scientists as well as the general public. The recent decades have seen substantial improvement in our understanding of their genetic bases and evolutionary relevance, revealing novel insights into the complex interplay of forces that influence these phenotypes. At the same time, many aspects remain poorly known, hampering a comprehensive understanding of these phenomena. Here we review the current state of this field and indicate topics that should be the focus of additional research. We devote particular attention to two aspects of mammalian pigmentation, melanism and pattern formation, highlighting recent advances and outstanding challenges, and proposing novel syntheses of available information. For both specific areas, and for pigmentation in general, we attempt to lay out recommendations for establishing novel model systems and integrated research programs that target the genetics and evolution of these phenotypes throughout the Mammalia.

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.

Comparative proteomics reveals genetic mechanisms underlying secondary hair follicle development in fine wool sheep during the fetal stage

The aim of this study was to investigate the genetic mechanisms underlying wool production by characterizing the skin protein profile and determining the proteomic changes that occur as a consequence of development in wool-producing sheep using a label-free proteomics approach. Samples were collected at four stages during gestation (87, 96, 102, and 138 days), and every two consecutive stages were statistically compared (87 versus 96, 96 versus 102, and 102 versus 138 days). We identified 227 specific proteins in the sheep proteome that were present in all four stages, and 123 differentially abundant proteins (DAPs). We also observed that the microstructure of the secondary follicles changed significantly during the development of the fetal skin hair follicle. The screened DAPs were strictly related to metabolic and skin development pathways, and were associated with pathways such as the glycolysis/gluconeogenesis. These analyses indicated that the wool production of fine wool sheep is regulated via a variety of pathways. These findings provide an important resource that can be used in future studies of the genetic mechanisms underlying wool traits in fine wool sheep, and the identified DAPs should be further investigated as candidate markers for predicting wool traits in sheep. SIGNIFICANCE: Wool quality (fiber diameter, length, etc.) is an important economic trait of fine wool sheep that is determined by secondary follicle differentiation and re-differentiation. Secondary follicles of fine wool sheep developed from a bud (87 days), and underwent differentiation (96 days) and rapid growth (102 days) until maturity (138 days) during gestation. Comparative analysis based on differential proteomics of these four periods could provide a better understanding of the wool growth mechanism of fine wool sheep and offer novel strategies for improving fine wool quality by breeding.

Expression and distribution of bone morphogenetic protein 4 and its antagonist Noggin in the skin of Kazakh sheep (Ovis aries) with a white and brown coat color

The natural coat color is an important trait of vertebrate animals. For example, the coat color can help avoid harm to human beings caused by chemical dyeing, and it has economic significance for domestic animals. The bone morphogenetic protein 4 (BMP4) and its antagonist Noggin can regulate pigmentation and the generation of coat color in mice; thus, they may also regulate the coat color of Kazakh sheep. To gain mechanistic insight into this possibility, we determined the relative expression levels of BMP4 and Noggin in the skin of white and brown Kazakh sheep by quantitative real-time polymerase chain reaction (qPCR) and western blotting analysis. The localization of BMP4 and Noggin were detected by immunohistochemistry (IHC). The results of qPCR and western blot analysis demonstrated that the relative expression levels of BMP4 and Noggin in the skin of brown Kazakh sheep were significantly higher than those in white Kazakh sheep. Our IHC results showed that the BMP4 protein was expressed in the epidermis and root sheath of the Kazakh sheep skin. The Noggin protein was expressed in the epidermis, root sheath, hair shaft, and dermal papilla of the Kazakh sheep skin. These results provide a theoretical basis for additional studies regarding the association and mechanism of BMP4 and Noggin in coat-color formation in Kazakh sheep. These results may provide new methods for developing treatment strategies for pigmentation disorders and diseases.

Transcriptome analysis identifies candidate genes associated with skin color variation in Triplophysa siluroides

In vertebrates, skin pigmentation is the most diverse phenotypic trait, and it is produced by a complex biological process that is often genetically controlled. Recently, two different colors (the typical brown and orange varieties) of Triplophysa siluroides, a species restricted to Yellow River drainage in China, were discovered. In the present study, the skin, brain and liver transcriptomes of T. siluroides of both colors were sequenced to search for genes related to skin pigmentation. Transcriptome sequencing generated 1,484,197,774 clean reads, resulting in a total of 222.6 Gb of sequence. The reads were assembled into 470,788 unigenes with a mean length of 1550 bp and an N50 size of 2944 bp. Functional annotation of the unigene dataset showed that 214,507, 304,161, 112,886, 179,074, 180,064, 184,837 and 82,081 unigenes were significantly matched to entries in the Nr protein, Nt, KO, Swiss-Prot, Pfam, GO and KOG databases, respectively. A differential expression analysis revealed that 2774, 3552 and 1529 unigenes were upregulated and 2720, 2663 and 1103 unigenes were downregulated in the skin, brain and liver of orange-skinned T. siluroides, respectively. Several genes that play key roles in pigmentation, i.e., Agouti, Slc45a2, Cbs, Mift and Slc7a11, showed significantly differential expression between brown and orange fish. In addition, we detected 158,863 simple sequence repeats (SSRs) in the T. siluroides transcriptome, and a total of 201,338 single-nucleotide polymorphisms (SNPs) were discovered in the different transcriptomes. The present results will facilitate further study of the molecular mechanisms of skin pigmentation and marker-assisted breeding of fish with valuable skin colors.

Transcriptional Differences of Coding and Non-Coding Genes Related to the Absence of Melanocyte in Skins of Bama Pig

Skin is the body’s largest organ, and the main function of skin is to protect underlying organs from possible external damage. Melanocytes play an important role in skin pigmentation. The Bama pig has a “two-end-black” phenotype with different coat colors across skin regions, e.g., white skin (without melanocytes) and black skin (with melanocytes), which could be a model to investigate skin-related disorders, specifically loss of melanocytes. Here, we generated expression profiles of mRNAs and long noncoding RNAs in Bama pig skins with different coat colors. In total, 14,900 mRNAs and 7549 lncRNAs were expressed. Overall, 2338 mRNAs/113 lncRNAs with FDR-adjusted p-value ≤ 0.05 were considered to be differentially expressed (DE) mRNAs/lncRNAs, with 1305 down-regulated mRNAs and 1033 up-regulated mRNAs in white skin with|log₂(fold change)| > 1. The genes down-regulated in white skin were associated with pigmentation, melanocyte–keratinocyte interaction, and keratin, while up-regulated ones were mainly associated with cellular energy metabolisms. Furthermore, those DE lncRNAs were predicted to be implicated in pigmentation, keratin synthesis and cellular energy metabolism. In general, this study provides insight into the transcriptional difference involved in melanocyte-loss-induced keratinocyte changes and promotes the Bama pig as a biomedical model in skin research.

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.

Expression and tissue distribution analysis of Angiotensin II in sheep (Ovis aries) skins associated with white and black coat colors

Angiotensin II (AngII) regulates pigment synthesis by tyrosinase in melanocytes. To evaluate the association between AngII and coat color formation, we detected the expression distribution of AngII in white and black sheep skins by LC-ESI-MS/MS, western blot, quantitative real-time-PCR (qPCR) and distribution of AngII by immunohistochemistry.Liquid chromatography-electrospray ionization tandem MS (LC-ESI-MS/MS) results showed that AngII was found in white and black skin tissues of sheep. Western blot results verified the LC-ESI-MS/MS results and suggested that AngII was expressed at significantly higher levels in black sheep skins compared with the white sheep skins. Quantitative real time PCR (qRT-PCR) results also revealed that the expression level of AngII mRNA was higher in black sheep skins than that in white sheep skins. Immunohistochemical analysis further demonstrated that AngII protein was localized in the hair bulb and outer root sheath of hair follicle in sheep. In summary, protein and transcripts exhibited the same expression pattern in white and black sheep skins. Furthermore, the expressions of AngII in the hair bulb and outer root sheath of black sheep were stronger than those in white sheep. These results suggested that AngII functions in sheep coat color regulation and offer a novel insight for further investigation on the role of AngII in the coat color formation in sheep.

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

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

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.

Contrasting patterns of coding and flanking region evolution in mammalian keratin associated protein-1 genes

Mammalian genomes contain a number of duplicated genes, and sequence identity between these duplicates can be maintained by purifying selection. However, between-duplicate recombination can also maintain sequence identity between copies, resulting in a pattern known as concerted evolution where within-genome repeats are more similar to each other than to orthologous repeats in related species. Here we investigated the tandemly-repeated keratin-associated protein 1 (KAP1) gene family, KRTAP1, which encodes proteins that are important components of hair and wool in mammals. Comparison of eutherian mammal KRTAP1 gene repeats within and between species shows a strong pattern of concerted evolution. However, in striking contrast to the coding regions of these genes, we find that the flanking regions have a divergent pattern of evolution. This contrast in evolutionary pattern transitions abruptly near the start and stop codons of the KRTAP1 genes. We reveal that this difference in evolutionary patterns is not explained by conventional purifying selection, nor is it likely a consequence of codon adaptation or reverse transcription of KRTAP1-n mRNA. Instead, the evidence suggests that these contrasting patterns result from short-tract gene conversion events that are biased to the KRTAP1 coding region by selection and/or differential sequence divergence. This work demonstrates the power that gene conversion has to finely shape the evolution of repetitive genes, and provides another distinctive pattern of contrasting evolutionary outcomes that results from gene conversion. A greater emphasis on exploring the evolution of multi-gene eukaryotic families will reveal how common different contrasting evolutionary patterns are in gene duplicates.

De Novo assembly and comparative transcriptome analyses of purple and green morphs of Apostichopus japonicus during body wall pigmentation process

Pigmentation processes provide a traceable and relevant trait for understanding key issues in evolutionary biology such as adaptation, speciation and the maintenance of balanced polymorphisms. The sea cucumber Apostichopus japonicus, which has nutritive and medical properties, is considered the most valuable commercial species in many parts of Asia. Compared with the green morph, the purple morph is rare and has great appeal to consumers. However, little is currently known about the molecular mechanism of body color formation in A. japonicus, even in echinoderm. Here, we employ illumina sequencing to examine expression patterns of the gene network underlying body wall development in purple and green morphs of A. japonicus. Overall, the number of down-regulated genes in the green morph was significantly more than in the purple morph during the pigmentation stage. We observed dynamic expression patterns of a large number of pigment, regulation and growth genes from the "Melanogenesis", "Melanoma", "Wnt signaling pathway", "Notch signaling pathway", "epithelium development", "epidermal growth factor receptor binding","growth factor activity" and "growth", including contrasting expression patterns of these genes in green and purple morph. This study provides comprehensive lists of differentially expressed genes during body wall development in the green and purple morphs, revealing potential candidate genes that may be involved in regulating body color formation and polymorphism. These data will provide valuable information for future genetic studies on sea cucumbers elucidating the molecular mechanisms underlying pigmentation, and may support the culture of desirable color morphs.

Genetic background of coat colour in sheep

Abstract. The coat colour of animals is an extremely important trait that affects their behaviour and is decisive for survival in the natural environment. In farm animal breeding, as a result of the selection of a certain coat colour type, animals are characterized by a much greater variety of coat types. This makes them an appropriate model in research in this field. A very important aspect of the coat colour types of farm animals is distinguishing between breeds and varieties based on this trait. Furthermore, for the sheep breeds which are kept for skins and wool, coat/skin colour is an important economic trait. Until now the study of coat colour inheritance in sheep proved the dominance of white colour over pigmented/black coat or skin and of black over brown. Due to the current knowledge of the molecular basis of ovine coat colour inheritance, there is no molecular test to distinguish coat colour types in sheep although some are available for other species, such as cattle, dogs, and horses. Understanding the genetic background of variation in one of the most important phenotypic traits in livestock would help to identify new genes which have a great effect on the coat colour type. Considering that coat colour variation is a crucial trait for discriminating between breeds (including sheep), it is important to broaden our knowledge of the genetic background of pigmentation. The results may be used in the future to determine the genetic pattern of a breed. Until now, identified candidate genes that have a significant impact on colour type in mammals mainly code for factors located in melanocytes. The proposed candidate genes code for the melanocortin 1 receptor (MC1R), agouti signaling protein (ASIP), tyrosinase-related protein 1 (TYRP1), microphthalmia-associated transcription factor MITF, and v-kit Hardy–Zuckerman 4 feline sarcoma viral oncogene homologue (KIT). However, there is still no conclusive evidence of established polymorphisms for specific coat colour types in sheep.

An examination of melanogenic traits and <i>TYRP1</i> polymorphism in Nanping and Romney Marsh sheep breeds

Abstract. The effects of mutations of the gene for tyrosinase-related protein 1 (TYRP1) on the black muscles and coat color in Nanping black-boned sheep were investigated. Tyrosinase activity and melanin content in plasma were measured and compared in three random groups of sheep: Nanping black-boned (101 heads), Nanping normal (106 heads) and Romney Marsh sheep (82 heads, Ovis aries). Eight exons and their partial flanking regions of the TYRP1 gene were amplified. Six intronic mutations and six exonic polymorphisms including two non-synonymous mutations [c.203C > T (p.A68V) and c.1202T > C (p.V401A)] were identified. Using a bi-directional polymerase chain reaction allele-specific amplification (bi-PASA) of the mutation c.203C > T it was shown that the frequencies of allele C in the Nanping black-boned, Nanping normal and Romney Marsh sheep were respectively 0.955, 0.967 and 0.744. For the mutation c.1202T > C, the frequencies of allele T in the three populations of sheep were respectively 0.777, 0.745 and 0.793 as measured using the single-strand conformation polymorphism. When the data from sheep of all three populations with the CC genotype of SNP c.203C > T were pooled, it was found that there was significantly higher (P < 0.05) tyrosinase activity, content of alkali-soluble melanin and ratio of eumelanin : total melanin than in the plasma of sheep with the CT and TT genotypes. This was not so within each of the three groups of sheep. No significant effect of the TRYP1 genotype on coat color was found. Further studies will be necessary to determine the cause of the black traits in Nanping black-boned sheep.

Evolution of Trichocyte Keratin Associated Proteins

The major components of hair are keratins and keratin associated proteins (KRTAPs). KRTAPs form the interfilamentous matrix between intermediate filament bundles through extensive disulfide bond cross-linking with the numerous cysteine residues in hair keratins. A variable number of approximately100-180 genes compose the KRTAP gene family in mammals. KRTAP gene family members present a typical pattern of concerted evolution, and its evolutionary features are consistent with the evolution of mammalian hair. KRATP genes might be more important in determining the structure of cashmere fibers in domestic mammals like sheep and goats. KRTAP gene variants thus should provide information for improved wool by sheep and goat breeding.

A comparison of transcriptomic patterns measured in the skin of Chinese fine and coarse wool sheep breeds

We characterised wool traits, and skin gene expression profiles of fine wool Super Merino (SM) and coarse wool Small Tail Han (STH) sheep. SM sheep had a significantly higher total density of wool follicles, heavier fleeces, finer fibre diameter, and increased crimp frequency, staple length and wool grease (lanolin) production. We found 435 genes were expressed at significantly different levels in the skin of the two breeds (127 genes more highly in SM and 308 genes more highly in STH sheep). Classification of the genes more highly expressed in SM sheep revealed numerous lipid metabolic genes as well as genes encoding keratins, keratin-associated proteins, and wool follicle stem cell markers. In contrast, mammalian epidermal development complex genes and other genes associated with skin cornification and muscle function were more highly expressed in STH sheep. Genes identified in this study may be further evaluated for inclusion in breeding programs, or as targets for therapeutic or genetic interventions, aimed at altering wool quality or yield. Expression of the lipid metabolic genes in the skin of sheep may be used as a novel trait with the potential to alter the content or properties of lanolin or the fleece.

Functional Role of Cyclin-Dependent Kinase 5 in the Regulation of Melanogenesis and Epidermal Structure

The mammalian integumentary system plays important roles in body homeostasis, and dysfunction of melanogenesis or epidermal development may lead to a variety of skin diseases, including melanoma. Skin pigmentation in humans and coat color in fleece-producing animals are regulated by many genes. Among them, microphthalmia-associated transcription factor (MITF) and paired-box 3 (PAX3) are at the top of the cascade and regulate activities of many important melanogenic enzymes. Here, we report for the first time that cyclin-dependent kinase 5 (Cdk5) is an essential regulator of MITF and PAX3. Cdk5 knockdown in mice causes a lightened coat color, a polarized distribution of melanin and hyperproliferation of basal keratinocytes. Reduced expression of Keratin 10 (K10) resulting from Cdk5 knockdown may be responsible for an abnormal epidermal structure. In contrast, overexpression of Cdk5 in sheep (Ovis aries) only produces brown patches on a white background, with no other observable abnormalities. Collectively, our findings show that Cdk5 has an important functional role in the regulation of melanin production and transportation and in normal development of the integumentary system.

Regulatory pathway analysis of coat color genes in Mongolian horses

BACKGROUND

Studies on the molecular genetics of horse skin pigmentation have typically focused on very few genes and proteins. In this study, we used Illumina sequencing to determine the global gene expression profiles in horses with white-colored coats and those with black-colored coats, with the goal of identifying novel genes that could regulate horse coat color.

RESULTS

Genes encoding ribosomal-associated proteins were highly expressed in horse skin. We found a total of 231 unigenes that were differentially expressed between horses with white coats and horses with black coats; 119 were down-regulated, and 112 were up-regulated. Many of the up-regulated genes in black horses, such as genes related to tyrosine metabolism, may directly regulate dark coat color. Keratin genes, MIA family genes, fatty acid-related genes, and melanoma-associated genes were also differentially regulated, which suggests that they may play important roles in coat color formation.

CONCLUSIONS

These findings show that the transcription profiles from white and black horse skin provide useful information to understand the genetics underlying the control of skin melanin synthesis in horses, which may enhance our knowledge of human skin diseases, such as melanoma and albinism.

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.

Gene expression profiling in Pekin duck embryonic breast muscle

Lean-type Pekin duck is a breed gained through long-term selection and great effort has been exerted to understand the mechanisms underlying increased muscle yields. However, the genes involved in Pekin duck embryonic breast muscle development have not been explored to date. In this study, we investigated gene expression profiles in Pekin Duck embryonic breast muscle at hatched day 13 (E13), E19, and E27 using RNA-seq. In total, we produced 519,312,178 raw reads resulting in 497,348,158 high-quality reads after filtering. The mapping, distribution of reads along annotated genes, and consistency across replicates demonstrates the high quality of the RNA-seq data used in this study, allowing us to continue with the downstream analysis. Significantly fewer differentially expressed genes (DEGs) were identified between E13 and E19 (203 DEGs) compared to E27 and E19 (2,797 DEGs). Many DEGs highly expressed in E19 are involved in metabolic processes and cell division. KEGG analysis showed many pathways associated with fat development were significantly enriched for DEGs highly expressed in E27. These results provide a basis for the further investigation of the mechanisms involved in Pekin duck embryonic breast muscle development.

The expression of KRT2 and its effect on melanogenesis in alpaca skins

In order to investigate the effects of the keratin 2 (KRT2) on alpaca melanocyte in vivo and vitro, the immunohistochemistry (IHC), quantitative real-time PCR (qPCR), Western blot, and alpaca melanocytes transfection methods were used. The results showed that mRNA and protein expression of KRT2 was highly expressed in brown skin in comparison with that in white skin. Moreover, we found that KRT2 was expressed in alpaca melanocytes in vitro by immunocytochemistry. After transfection with KRT2 in alpaca melanocytes, the relative mRNA and protein expression of KRT2, microphthalmia-associtated transcription factor (MITF), tyrosinase (TYR) and tyrosinase-related protein 1 (TYRP1) in alpaca skin melanocytes was increased with significant differences; a further result was the increase of melanin production. The results suggested that KRT2 functions in alpaca hair color formation, which offered an essential theoretical basis for further exploration of the role of melanogenesis.

Comparative transcriptome analyses of seven anurans reveal functions and adaptations of amphibian skin

Animal skin, which is the tissue that directly contacts the external surroundings, has evolved diverse functions to adapt to various environments. Amphibians represent the transitional taxon from aquatic to terrestrial life. Exploring the molecular basis of their skin function and adaptation is important to understand the survival and evolutionary mechanisms of vertebrates. However, comprehensive studies on the molecular mechanisms of skin functions in amphibians are scarce. In this study, we sequenced the skin transcriptomes of seven anurans belonging to three families and compared the similarities and differences in expressed genes and proteins. Unigenes and pathways related to basic biological processes and special functions, such as defense, immunity, and respiration, were enriched in functional annotations. A total of 108 antimicrobial peptides were identified. The highly expressed genes were similar in species of the same family but were different among families. Additionally, the positively selected orthologous groups were involved in biosynthesis, metabolism, immunity, and defense processes. This study is the first to generate extensive transcriptome data for the skin of seven anurans and provides unigenes and pathway candidates for further studies on amphibian skin function and adaptation.

Characterization and expression of soluble guanylate cyclase in skins and melanocytes of sheep

The study reported the characterization of soluble guanylate cyclase (sGC) with the size of CDS of 1860bp, encoding a protein of 620 amino acids and containing several conserved functional domains including HNOB, HNOBA, and CHD. Quantitative real time PCR analysis of sGC showed that the expression of sGC mRNA is higher (∼5 fold) in white sheep skin relative to black sheep skin with significant difference (P<0.01). Using a rabbit polyclonal anti-sGC antibody, an immune reactive band corresponding to sheep sGC protein was detected in the skin samples by Western blotting analysis, and the expression of sGC protein was significantly higher in white sheep skin compared to black sheep skin (P<0.01). Immunohistochemical analysis revealed that sGC protein was localized in cytoplasm and intercellular substance of upper hair papilla in hair follicles of white sheep skin, but the protein was localized in cytoplasm and intercellular substance of lower hair bulb and outer root sheath cells in hair follicles of black sheep skin. The immunocytochemical analysis revealed that sGC was expressed in melanocytes in vitro of sheep skin. Over expression of sGC in melanocytes resulted in decreased expression of key melanogenic genes including microphthalmia transcription factor (MITF), tyrosinase (TYR), tyrosinase related protein 1(TYRTP1), and tyrosinase related protein 2(TYRP2) both at mRNA and protein level. Moreover, the melanocytes was capable of producing cGMP and cAMP. The observed differential expression and localization of sGC in sheep skins and melanocytes and the capability of producing cGMP and cAMP, which suggested a potential role for this gene in hair color regulation.

Comparative Transcriptome Analysis of Fetal Skin Reveals Key Genes Related to Hair Follicle Morphogenesis in Cashmere Goats

Cashmere goat skin contains two types of hair follicles (HF): primary hair follicles (PHF) and secondary hair follicles (SHF). Although multiple genetic determinants associated with HF formation have been identified, the molecules that determine the independent morphogenesis of HF in cashmere goats remain elusive. The growth and development of SHF directly influence the quantity and quality of cashmere production. Here, we report the transcriptome profiling analysis of nine skin samples from cashmere goats using 60- and 120-day-old embryos (E60 and E120, respectively), as well as newborns (NB), through RNA-sequencing (RNA-seq). HF morphological changes indicated that PHF were initiated at E60, with maturation from E120, while differentiation of SHF was identified at E120 until formation of cashmere occurred after birth (NB). The RNA-sequencing analysis generated over 20.6 million clean reads from each mRNA library. The number of differentially expressed genes (DEGs) in E60 vs. E120, E120 vs. NB, and E60 vs. NB were 1,024, 0 and 1,801, respectively, indicating that no significant differences were found at transcriptomic levels between E120 and NB. Key genes including B4GALT4, TNC, a-integrin, and FGFR1, were up-regulated and expressed in HF initiation from E60 to E120, while regulatory genes such as GPRC5D, PAD3, HOXC13, PRR9, VSIG8, LRRC15, LHX2, MSX-2, and FOXN1 were up-regulated and expressed in HF keratinisation and hair shaft differentiation from E120 and NB to E60. Several genes belonging to the KRT and KRTAP gene families were detected throughout the three HF developmental stages. The transcriptional trajectory analyses of all DEGs indicated that immune privilege, glycosaminoglycan biosynthesis, extracellular matrix receptor interaction, and growth factor receptors all played dominant roles in the epithelial-mesenchymal interface and HF formation. We found that the Wnt, transforming growth factor-beta/bone morphogenetic protein, and Notch family members played vital roles in HF differentiation and maturation. The DEGs we found could be attributed to the generation and development of HF, and thus will be critically important for improving the quantity and quality of fleece production in animals for fibres.

Characterization and comparative analyses of transcriptomes for in vivo and in vitro produced peri-implantation conceptuses and endometria from sheep

An increasing number of reports indicate that in vitro fertilization (IVF) is highly associated with long‑term side effects on embryonic and postnatal development, and can sometimes result in embryonic implant failure. While high‑throughput gene expression analysis has been used to explore the mechanisms underlying IVF-induced side effects on embryonic development, little is known about the effects of IVF on conceptus–endometrial interactions during the peri-implantation period. Using sheep as a model, we performed a comparative transcriptome analysis between in vivo (IVO; in vivo fertilized followed by further development in the uterus) and in vitro produced (IVP; IVF with further culture in the incubator) conceptuses, and the caruncular and intercaruncular areas of the ovine endometrium. We identified several genes that were differentially expressed between the IVO and IVP groups on day 17, when adhesion between the trophoblast and the uterine luminal epithelium begins in sheep. By performing Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, we found that, in the conceptus, differentially expressed genes (DEGs) were associated mainly with functions relating to cell binding and the cell cycle. In the endometrial caruncular area, DEGs were involved in cell adhesion/migration and apoptosis, and in the intercaruncular area, they were significantly enriched in pathways of signal transduction and transport. Thus, these DEGs are potential candidates for further exploring the mechanism underlying IVF/IVP-induced embryonic implant failure that occurs due to a loss of interaction between the conceptus and endometrium during the peri-implantation period.

Identification of differentially expressed Gnαs and Gnα11 in sheep (Ovis aries) skins associated with white and black coat colors

Guanine nucleotide-binding protein subunit alpha-s (Gnαs) and guanine nucleotide-binding protein subunit alpha-11 (Gnα11) play an important role in coat color formation. To improve our understanding of Gnαs and Gnα11 expression levels and pattern in the skin of black sheep (Ovis aries) compared with white sheep, we analyzed the expression levels through quantitative real time PCR (qPCR) and Western blot, immunohistochemistry and immunofluorescence. qPCR and Western blot results suggested that Gnαs and Gnα11 were significantly expressed at high levels in black sheep skin compared with the white sheep skin. Gnα11 expression was higher than Gnαs expression in both skin colors, transcripts and protein exhibited the same expression pattern in white and black sheep skins. Immunohistochemical results revealed that Gnαs and Gnα11 were localized in the outer root sheath of hair follicle in sheep. Furthermore, the expression of Gnα11 in outer root sheath of hair follicle was stronger than that in Gnαs. Immunofluorescencence further demonstrated that signals of Gnαs and Gnα11 were detected in outer root sheath and hair papilla. These results provide a novel insight into the role of Gnαs and Gnα11 in the regulation of sheep coat color.