Gene expression analysis identifies new candidate genes associated with the development of black skin spots in Corriedale sheep

Novel Insights to Assess Climate Resilience in Goats Using a Holistic Approach of Skin-Based Advanced NGS Technologies

A novel study was conducted to elucidate heat-stress responses on a number of hair- and skin-based traits in two indigenous goat breeds using a holistic approach that considered a number of phenotypic and genomic variables. The two goat breeds, Kanni Aadu and Kodi Aadu, were subjected to a simulated heat-stress study using the climate chambers. Four groups consisting of six goats each (KAC, Kaani Aadu control; KAH, Kanni Aadu heat stress; KOC, Kodi Aadu control; and KOH, Kodi Aadu heat stress) were considered for the study. The impact of heat stress on caprine skin tissue along with a comparative assessment of the thermal resilience of the two goat breeds was assessed. The variables considered were hair characteristics, hair cortisol, hair follicle quantitative PCR (qPCR), sweating (sweating rate and active sweat gland measurement), skin histometry, skin-surface infrared thermography (IRT), skin 16S rRNA V3-V4 metagenomics, skin transcriptomics, and skin bisulfite sequencing. Heat stress significantly influenced the hair fiber characteristics (fiber length) and hair follicle qPCR profile (Heat-shock protein 70 (HSP70), HSP90, and HSP110). Significantly higher sweating rate, activated sweat gland number, skin epithelium, and sweat gland number (histometry) were observed in heat stressed goats. The skin microbiota was also observed to be significantly altered due to heat stress, with a relatively higher alteration being noticed in Kanni Aadu goats than in Kodi Aadi goats. Furthermore, the transcriptomics and epigenetics analysis also pointed towards the significant impact of heat stress at the cellular and molecular levels in caprine skin tissue. The higher proportion of differentially expressed genes (DEGs) along with higher differentially methylated regions (DMRs) in Kanni Aadu goats due to heat stress when compared to Kodi Aadu goats pointed towards the better resilience of the latter breed. A number of established skin, adaptation, and immune-response genes were also observed to be significantly expressed/methylated. Additionally, the influence of heat stress at the genomic level was also predicted to result in significant functional alterations. This novel study thereby highlights the impact of heat stress on the caprine skin tissue and also the difference in thermal resilience exhibited by the two indigenous goat breeds, with Kodi Aadu goats being more resilient.

Genome-wide association analysis of nine reproduction and morphological traits in three goat breeds from Southern China

OBJECTIVE

This study aimed to investigate the significant single nucleotide polymorphisms (SNPs) and genes associated with nine reproduction and morphological traits in three breed populations of Chinese goats.

METHODS

The genome-wide association of nine reproduction and morphological traits (litter size, nipple number, wattle, skin color, coat color, black dorsal line, beard, beard length, and hind leg hair) were analyzed in three Chinese native goat breeds (n = 336) using an Illumina Goat SNP50 Beadchip.

RESULTS

A total of 17 genome-wide or chromosome-wide significant SNPs associated with one reproduction trait (litter size) and six morphological traits (wattle, coat color, black dorsal line, beard, beard length, and hind leg hair) were identified in three Chinese native goat breeds, and the candidate genes were annotated. The significant SNPs and corresponding putative candidate genes for each trait are as follows: two SNPs located on chromosomes 6 (CSN3) and 24 (TCF4) for litter size trait; two SNPs located on chromosome 9 (KATNA1) and 1 (UBASH3A) for wattle trait; three SNPs located on chromosome 26 (SORCS3), 24 (DYM), and 20 (PDE4D) for coat color trait; two SNPs located on chromosome 18 (TCF25) and 15 (CLMP) for black dorsal line trait; four SNPs located on chromosome 8, 2 (PAX3), 5 (PIK3C2G), and 28 (PLA2G12B and OIT3) for beard trait; one SNP located on chromosome 18 (KCNG4) for beard length trait; three SNPs located on chromosome 17 (GLRB and GRIA2), 28 (PGBD5), and 4 for hind leg hair trait. In contrast, there were no SNPs identified for nipple number and skin color.

CONCLUSION

The significant SNPs or genes identified in this study provided novel insights into the genetic mechanism underlying important reproduction and morphological traits of three local goat breeds in Southern China as well as further potential applications for breeding goats.

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.

Identification of Potential Blind-Side Hypermelanosis-Related lncRNA–miRNA–mRNA Regulatory Network in a Flatfish Species, Chinese Tongue Sole (Cynoglossus semilaevis)

Blind-side hypermelanosis has emerged as a major concern in commercial rearing environments of the flatfish aquaculture industry. To date, the underlying molecular mechanisms are not well understood. To fill this gap, in this study, whole transcriptomic sequencing and analyses were performed using normal skins and hypermelanic skins of the blind side of Chinese tongue sole (Cynoglossus semilaevis). Differentially expressed long non-coding RNAs (DElncRNAs), miRNAs (DEmiRNAs), and differentially expressed genes as well as their competing endogenous RNA (ceRNA) networks were identified. A total of 34 DElncRNAs, 226 DEmiRNAs, and 610 DEGs were identified. Finally, lncRNA–miRNA–mRNA regulatory networks (involving 29 DElncRNAs, 106 DEmiRNAs, and 162 DEGs) associated with blind-side hypermelanosis were constructed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of 162 DEGs in ceRNA networks identified DEGs (e.g., oca2, mc1r, and ihhb) in pigmentation-related biological processes and DEGs (e.g., ca4, glul, and fut9) in nitrogen metabolism, glycosphingolipid biosynthesis, and folate biosynthesis pathways, as well as their corresponding DElncRNAs and DEmiRNAs to potentially play key regulatory roles in blind-side hypermelanosis. In conclusion, this is the first study on the ceRNA regulatory network associated with blind-side hypermelanosis in flatfish. These new findings expand the spectrum of non-coding regulatory mechanisms underpinning blind-side hypermelanosis, which facilitates the further exploration of molecular regulatory mechanisms of malpigmentation in flatfish.

The role of KLF4 in melanogenesis and homeostasis in sheep melanocytes

KLF4 expression has been associated with hair color in mammals and has also been found to regulate melanoma cell growth. Here, we assessed the influence of KLF4 on coat color formation and melanocytes. We found that KLF4 was highly expressed in the black skin of sheep both at the mRNA and protein levels compared with white skin. KLF4 immunostaining further showed that KLF4 protein was mainly expressed in epidermal, outer root, and hair bulb regions. In sheep melanocytes, the proliferation of melanocytes was inhibited by KLF4 overexpression and this decrease in cell proliferation was coupled with induction of the S phase, cell cycle arrest, and apoptosis. In vitro cell migration assays showed that KLF4 suppressed cell migration. In addition, KLF4 overexpression significantly increased melanin production and pigment-related gene expression. Collectively, our findings show that KLF4 is important for coat color formation and melanocyte homeostasis.

Discovery of genes and proteins possibly regulating mean wool fibre diameter using cDNA microarray and proteomic approaches

Wool fibre diameter (WFD) is one of the wool traits with higher economic impact. However, the main genes specifically regulating WFD remain unidentified. In this current work we have used Agilent Sheep Gene Expression Microarray and proteomic technology to investigate the gene expression patterns of body side skin, bearing more wool, in Aohan fine wool sheep, a Chinese indigenous breed, and compared them with that of small tail Han sheep, a sheep bread with coarse wool. Microarray analyses showed that most of the genes likely determining wool diameter could be classified into a few categories, including immune response, regulation of receptor binding and growth factor activity. Certain gene families might play a role in hair growth regulation. These include growth factors, immune cytokines, solute carrier families, cellular respiration and glucose transport amongst others. Proteomic analyses also identified scores of differentially expressed proteins.

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.

Skin transcriptome profiles associated with black- and white-coated regions in Boer and Macheng black crossbred goats

To increase the current understanding of the gene-expression profiles in different skin regions associated with different coat colors and identify key genes for the regulation of color patterns in goats, we used the Illumina RNA-Seq method to compare the skin transcriptomes of the black- and white-coated regions containing hair follicles from the Boer and Macheng Black crossbred goat, which has a black head and a white body. Six cDNA libraries derived from skin samples of the white-coated region (n = 3) and black-coated region (n = 3) were constructed from three full-sib goats. On average, we obtained approximately 76.5 and 73.5 million reads for skin samples from black- and white-coated regions, respectively, of which 75.39% and 76.05% were covered in the genome database. A total of 165 differentially expressed genes (DEGs) were detected between these two color regions, among which 110 were upregulated and 55 were downregulated in the skin samples of white- vs. black-coated regions. The results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that some of these DEGs may play an important role in controlling the pigmentation of skin or hair follicles. We identified three key DEGs, i.e., Agouti, DCT, and TYRP1, in the pathway related to melanogenesis in the different skin regions of the crossbred goat. DCT and TYRP1 were downregulated and Agouti was upregulated in the skin of the white-coated region, suggesting a lack of mature melanocytes in this region and that Agouti might play a key developmental role in color-pattern formation. All data sets (Gene Expression Omnibus) are available via public repositories. In addition, MC1R was genotyped in 200 crossbred goats with a black head and neck. Loss-of-function mutations in MC1R as well as homozygosity for the mutant alleles were widely found in this population. The MC1R gene did not seem to play a major role in determining the black head and neck in our crossbred goats. Our study provides insights into the transcriptional regulation of two distinct coat colors, which might serve as a key resource for understanding coat color pigmentation in goats. The region-specific expression of Agouti may be associated with the distribution of pigments across the body in Boer and Macheng Black crossbred goats.

Exploring differentially expressed genes associated with coat color in goat skin using RNA-seq

Fur color in domestic goats is an important, genetically determined characteristic that is associated with economic value. This study was designed to perform a comprehensive expression profiling of genes expressed in the skin tissues from Laiwu Black goat and Lubei White goat. Comparisons of black and white goat skin transcriptomes revealed 102 differentially expressed genes (DEGs), of which 38 were upregulated and 64 downregulated in black skin compared with white skin. Among the DEGs, we identified six genes involved in pigmentation, including agouti signaling protein (ASIP), CAMP responsive element binding protein 3-like 1 (CREB3L1), dopachrome tautomerase (DCT), premelanosome protein (PMEL), transient receptor potential cation channel subfamily M member 1 (TRPM1), and tyrosinase-related protein 1 (TYRP1). Notably, there were no significant differences in the expression of melanocortin 1 receptor, microphthalmia-associated transcription factor, tyrosinase, and KIT proto-oncogene receptor tyrosine kinase between the black and white skin samples, whereas ASIP expression was detected only in white skin. PMEL, TRPM1, TYRP1, and DCT showed higher expression in black goat skin, but ASIP and CREB3L1 had higher expression in white goat skin. Quantitative polymerase chain reaction results for PMEL, TRPM1, DCT, TYRP1, and CREB3L1 expression were consistent with those for RNA-seq. These results will expand our understanding of the complex molecular mechanisms of skin physiology and melanogenesis in goats, and provide a foundation for future studies.

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.

Identification of genes and proteins associated with anagen wool growth

Identifying genes of major effect for wool growth would offer strategies for improving the quality and increasing the yield of fine wool. In this study, we employed the Agilent Sheep Gene Expression Microarray and proteomic technology to investigate the gene expression patterns of body side skin (more wool growing) in Aohan fine wool sheep (a Chinese indigenous breed) in comparison with groin skin (no wool growing) at the anagen stage of the wool follicle. A microarray study revealed that 4772 probes were differentially expressed, including 2071 upregulated and 2701 downregulated probes, in the comparisons of body side skin vs. groin skin (S/G). The microarray results were verified by means of quantitative PCR. A total of 1099 probes were assigned to unique genes/transcripts. The number of distinct genes/transcripts (annotated) was 926, of which 352 were upregulated and 574 were downregulated. In S/G, 13 genes were upregulated by more than 10 fold, whereas 60 genes were downregulated by more than 10 fold. Further analysis revealed that the majority of the genes possibly related to the wool growth could be assigned to categories including regulation of cell division, intermediate filament, cytoskeletal part and growth factor activity. Several potential gene families may participate in hair growth regulation, including fibroblast growth factors, transforming growth factor-β, WNTs, insulin-like growth factor, vascular endothelial growth factors and so on. Proteomic analysis also revealed 196 differentially expressed protein points, of which 121 were identified as single protein points.

Skin transcriptome profiles associated with skin color in chickens

Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.

Identification of skin-expressed genes possibly associated with wool growth regulation of Aohan fine wool sheep

Background

Sheep are valuable resources for the animal fibre industry. Therefore, identifying genes which regulate wool growth would offer strategies for improving the quality of fine wool. In this study, we employed Agilent sheep gene expression microarray and proteomic technology to compare the gene expression patterns of the body side (hair-rich) and groin (hairless) skins of Aohan fine wool sheep (a Chinese indigenous breed).

Results

Comparing the body side to the groin skins (S/G) of Aohan fine wool sheep, the microarray study revealed that 1494 probes were differentially expressed, including 602 more highly expressed and 892 less highly expressed probes. The microarray results were verified by means of quantitative PCR. Cluster analysis could distinguish the body side skin and the groin skin. Based on the Database for Annotation, Visualization and Integrated Discovery (DAVID), 38 of the differentially expressed genes were classified into four categories, namely regulation of receptor binding, multicellular organismal process, protein binding and macromolecular complex. Proteomic study revealed that 187 protein spots showed significant (p < 0.05) differences in their respective expression levels. Among them, 46 protein entries were further identified by MALDI-TOF/MS analyses.

Conclusions

Microarray analysis revealed thousands of differentially expressed genes, many of which were possibly associated with wool growth. Several potential gene families might participate in hair growth regulation. Proteomic analysis also indentified hundreds of differentially expressed proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0144-1) contains supplementary material, which is available to authorized users.

Differential expression of genes and proteins associated with wool follicle cycling

Sheep are valuable resources for the wool industry. Wool growth of Aohan fine wool sheep has cycled during different seasons in 1 year. Therefore, identifying genes that control wool growth cycling might lead to ways for improving the quality and yield of fine wool. In this study, we employed Agilent sheep gene expression microarray and proteomic technology to compare the gene expression patterns of the body side skins at August and December time points in Aohan fine wool sheep (a Chinese indigenous breed). Microarray study revealed that 2,223 transcripts were differentially expressed, including 1,162 up-regulated and 1,061 down-regulated transcripts, comparing body side skin at the August time point to the December one (A/D) in Aohan fine wool sheep. Then seven differentially expressed genes were selected to validated the reliability of the gene chip data. The majority of the genes possibly related to follicle development and wool growth could be assigned into the categories including regulation of receptor binding, extracellular region, protein binding and extracellular space. Proteomic study revealed that 84 protein spots showed significant differences in expression levels. Of the 84, 63 protein spots were upregulated and 21 were downregulated in A/D. Finally, 55 protein points were determined through MALDI-TOF/MS analyses. Furthermore, the regulation mechanism of hair follicle might resemble that of fetation.

Skin transcriptome profiles associated with coat color in sheep

Background

Previous molecular genetic studies of physiology and pigmentation of sheep skin have focused primarily on a limited number of genes and proteins. To identify additional genes that may play important roles in coat color regulation, Illumina sequencing technology was used to catalog global gene expression profiles in skin of sheep with white versus black coat color.

Results

There were 90,006 and 74,533 unigenes assembled from the reads obtained from white and black sheep skin, respectively. Genes encoding for the ribosomal proteins and keratin associated proteins were most highly expressed. A total of 2,235 known genes were differentially expressed in black versus white sheep skin, with 479 genes up-regulated and 1,756 genes down-regulated. A total of 845 novel genes were differentially expressed in black versus white sheep skin, consisting of 107 genes which were up-regulated (including 2 highly expressed genes exclusively expressed in black sheep skin) and 738 genes that were down-regulated. There was also a total of 49 known coat color genes expressed in sheep skin, from which 13 genes showed higher expression in black sheep skin. Many of these up-regulated genes, such as DCT, MATP, TYR and TYRP1, are members of the components of melanosomes and their precursor ontology category.

Conclusion

The white and black sheep skin transcriptome profiles obtained provide a valuable resource for future research to understand the network of gene expression controlling skin physiology and melanogenesis in sheep.

Genome-wide transcriptome analysis between small-tail Han sheep and the Surabaya fur sheep using high-throughput RNA sequencing

The small-tail Han sheep and the Surabaya fur sheep are two local breeds in north China, which are characterized by high-fecundity and low-prolificacy breed respectively. Significant genetic differences between these two breeds have provided increasing interests in the identification and utilization of major prolificacy genes in these sheep. High prolificacy is a complex trait, and it is difficult to comprehensively identify the candidate genes related to this trait using the single molecular biology technique. To understand the molecular mechanisms of fecundity and provide more information about high prolificacy candidate genes in high- and low-fecundity sheep, we explored the utility of next-generation sequencing technology in this work. A total of 1.8 Gb sequencing reads were obtained and resulted in more than 20 000 contigs that averaged ∼300 bp in length. Ten differentially expressed genes were further verified by quantitative real-time RT-PCR to confirm the reliability of RNA-seq results. Our work will provide a basis for the future research of the sheep reproduction.