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

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.

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.

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

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

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.

Identifying Candidate Genes for Litter Size and Three Morphological Traits in Youzhou Dark Goats Based on Genome-Wide SNP Markers

This study aimed to reveal the potential genetic basis for litter size, coat colour, black middorsal stripe and skin colour by combining genome-wide association analysis (GWAS) and selection signature analysis and ROH detection within the Youzhou dark (YZD) goat population (n = 206) using the Illumina GoatSNP54 BeadChip. In the GWAS, we identified one SNP (snp54094-scaffold824-899720) on chromosome 11 for litter size, two SNPs on chromosome 26 (snp11508-scaffold142-1990450, SORCS3) and chromosome 12 (snp55048-scaffold842-324525, LOC102187779) for coat colour and one SNP on chromosome 18 (snp56013-scaffold873-22716, TCF25) for the black middorsal stripe. In contrast, no SNPs were identified for skin colour. In selection signature analysis, 295 significant iHS genomic regions with a mean |iHS| score > 2.66, containing selection signatures encompassing 232 candidate genes were detected. In particular, 43 GO terms and one KEGG pathway were significantly enriched in the selected genes, which may contribute to the excellent environmental adaptability and characteristic trait formation during the domestication of YZD goats. In ROH detection, we identified 4446 ROH segments and 282 consensus ROH regions, among which nine common genes overlapped with those detected using the iHS method. Some known candidate genes for economic traits such as reproduction (TSHR, ANGPT4, CENPF, PIBF1, DACH1, DIS3, CHST1, COL4A1, PRKD1 and DNMT3B) and development and growth (TNPO2, IFT80, UCP2, UCP3, GHRHR, SIM1, CCM2L, CTNNA3 and CTNNA1) were revealed by iHS and ROH detection. Overall, this study is limited by the small population size, which affects the results of GWAS to a certain extent. Nevertheless, our findings could provide the first overview of the genetic mechanism underlying these important traits and provide novel insights into the future conservation and utilisation of Chinese goat germplasm resources.

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.