Gene network analysis reveals candidate genes related with the hair follicle development in sheep

Study on the expression patterns of inner root sheath-specific genes in Tan sheep hair follicle during different developmental stages

By analyzing the expression patterns of inner root sheath (IRS) specific genes during different developmental stages of hair follicle (HF) in Tan sheep embryos and at birth, this study aims to reveal the influence of the IRS on crimped wool. Skin tissues from the scapular region of male Tan sheep were collected at 85 days (E85) and 120 days (E120) of fetal development, and at 0 days (D0), 35 days (D35), and 60 days (D60) after birth, with four samples at each stage. Real-time quantitative polymerase chain reaction (RT-qPCR) was employed to determine the relative expression levels of IRS type I keratin genes (KRT25, KRT26, KRT27, KRT28), type II keratin genes (KRT71, KRT72, KRT73, KRT74), and the trichohyalin gene (TCHH) in the skin of Tan sheep at different stages. Results showed that the expression levels of all IRS-specific genes peaked at D0, with the expression of all genes significantly higher than at E85 (P < 0.01), except for KRT73 and TCHH. The expression levels of KRT25, KRT26, and KRT72 were also significantly higher than at E120 (P < 0.01). Furthermore, the expression levels of KRT27, KRT28, KRT71, and KRT74 were significantly higher than both at E120 and D35 (P < 0.01). The expression levels of other genes at different stages showed no significant difference (P > 0.05). Conclusion: The IRS-specific genes exhibit the highest expression levels in Tan sheep at the neonatal stage. The expression levels of KRT71, KRT72, and TCHH, which are consistent with the pattern of wool crimp, may influence the morphology of the IRS and thereby affect the crimp of Tan sheep wool.

Population structure and selection signal analysis of indigenous sheep from the southern edge of the Taklamakan Desert

Analyzing the genetic diversity and selection characteristics of sheep (Ovis aries) holds significant value in understanding their environmental adaptability, enhancing breeding efficiency, and achieving effective conservation and rational utilization of genetic resources. In this study, we utilized Illumina Ovine SNP 50 K BeadChip data from four indigenous sheep breeds from the southern margin of the Taklamakan Desert (Duolang sheep: n = 36, Hetian sheep: n = 74, Kunlun sheep: n = 27, Qira black sheep: n = 178) and three foreign meat sheep breeds (Poll Dorset sheep: n = 105, Suffolk sheep: n = 153, Texel sheep: n = 150) to investigate the population structure, genetic diversity, and genomic signals of positive selection within the indigenous sheep. According to the Principal component analysis (PCA), the Neighbor-Joining tree (NJ tree), and Admixture, we revealed distinct clustering patterns of these seven sheep breeds based on their geographical distribution. Then used Cross Population Extended Haplotype Homozygosity (XP-EHH), Fixation Index (FST), and Integrated Haplotype Score (iHS), we identified a collective set of 32 overlapping genes under positive selection across four indigenous sheep breeds. These genes are associated with wool follicle development and wool traits, desert environmental adaptability, disease resistance, reproduction, and high-altitude adaptability. This study reveals the population structure and genomic selection characteristics in the extreme desert environments of native sheep breeds from the southern edge of the Taklimakan Desert, providing new insights into the conservation and sustainable use of indigenous sheep genetic resources in extreme environments. Additionally, these findings offer valuable genetic resources for sheep and other mammals to adapt to global climate change.

Key genes and metabolites that regulate wool fibre diameter identified by combined transcriptome and metabolome analysis

BACKGROUND

Fibre diameter is an important economic trait of wool fibre. As the fibre diameter decreases, the economic value of wool increases. Therefore, understanding the mechanism of wool fibre diameter regulation is important in improving the value of wool.

RESULTS

In this study, we used non-targeted metabolome and reference transcriptome data to detect differences in metabolites and genes in groups of Alpine Merino sheep with different wool fibre diameter gradients, and integrated metabolome and transcriptome data to identify key genes and metabolites that regulate wool fibre diameter. We found 464 differentially abundant metabolites (DAMs) and 901 differentially expressed genes (DEGs) in four comparisons of groups with different wool fibre diameters. Approximately 25% of the differentially abundant metabolites were lipid and lipid-like molecules. These molecules were predicted to be associated with skin development and keratin filament by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses. Key genes, including COL5A2, COL5A3, CREB3L4, COL1A1, and SFRP4, were identified by gene set enrichment analysis.

CONCLUSIONS

Key genes regulating wool fibre diameter were identified, the effects of lipid molecules on wool performance were investigated, and potential synergies between genes and metabolites were postulated, providing a theoretical framework for fine wool sheep breeding.

Transcriptome meta-analysis reveals the hair genetic rules in six animal breeds and genes associated with wool fineness

Wool plays an irreplaceable role in the lives of livestock and the textile industry. The variety of hair quality and shape leads to the diversity of its functions and applications, and the finer wool has a higher economic value. In this study, 10 coarse and 10 fine ordos fine wool sheep skin samples were collected for RNA-seq, and coarse and fine skin/hair follicle RNA-seq datasets of other five animal breeds were obtained from NCBI. Weighted gene co-expression network analysis showed that the common genes were clustered into eight modules. Similar gene expression patterns in sheep and rabbits with the same wool types, different gene expression patterns in animal species with different hair types, and brown modules were significantly correlated with species and breeds. GO and KEGG enrichment analyses showed that, most genes in the brown module associated with hair follicle development. Hence, gene expression patterns in skin tissues may determine hair morphology in animal. The analysis of differentially expressed genes revealed that 32 highly expressed candidate genes associated with the wool fineness of Ordos fine wool sheep. Among them, KAZALD1 (grey module), MYOC (brown module), C1QTNF6 (brown module), FOS (tan module), ITGAM, MX2, MX1, and IFI6 genes have been reported to be involved in the regulation of the hair follicle cycle or hair loss. Additionally, 12 genes, including KAZALD1, MYOC, C1QTNF6, and FOS, are differentially expressed across various animal breeds and species. The above results suggest that different sheep breeds share a similar molecular regulatory basis of wool fineness. Finally, the study provides a theoretical reference for molecular breeding of sheep breeds as well as for the investigation of the origin and evolution of animal hair.

MiR-23b and miR-133 Cotarget TGFβ2/NOTCH1 in Sheep Dermal Fibroblasts, Affecting Hair Follicle Development

Wool is produced and controlled by hair follicles (HFs). However, little is known about the mechanisms involved in HF development and regulation. Sheep dermal fibroblasts (SDFs) play a key role in the initial stage of HF development. Analyzing the molecular mechanism that regulates early HF development in superfine wool sheep is of great importance for better understanding the HF morphogenesis process and for the breeding of fine wool sheep. Here, we show that two microRNAs (miRNAs) affect the development of HFs by targeting two genes that are expressed by SDFs. Meanwhile, the overexpression and inhibition of oar-miR-23b and oar-miR-133 in SDFs cells and cell proliferation, apoptosis, and migration were further detected using a CCK-8 assay, an Annexin V-FITC assay, a Transwell assay, and flow cytometry. We found that oar-miR-23b, oar-miR-133, and their cotarget genes TGFβ2 and NOTCH1 were differentially expressed during the six stages of HF development in superfine wool sheep. Oar-miR-23b and oar-miR-133 inhibited the proliferation and migration of SDFs and promoted the apoptosis of SDFs through TGFβ2 and NOTCH1. oar-miR-23b and oar-miR-133 inhibited the proliferation and migration of SDFs by jointly targeting TGFβ2 and NOTCH1, thereby inhibiting the development of superfine wool HFs. Our research provides a molecular marker that can be used to guide the breeding of ultrafine wool sheep.

Two mutations at KRT74 and EDAR synergistically drive the fine-wool production in Chinese sheep

INTRODUCTION

Fine-wool sheep are the most common breed used by the wool industry worldwide. Fine-wool sheep have over a three-fold higher follicle density and a 50% smaller fiber diameter than coarse-wool sheep.

OBJECTIVES

This study aims to clarify the underlying genetic basis for the denser and finer wool phenotype in fine-wool breeds.

METHOD

Whole-genome sequences of 140 samples, Ovine HD630K SNP array data of 385 samples, including fine, semi-fine, and coarse wool sheep, as well as skin transcriptomes of nine samples were integrated for genomic selection signature analysis.

RESULTS

Two loci at keratin 74 (KRT74) and ectodysplasin receptor (EDAR) were revealed. Fine-scale analysis in 250 fine/semi-fine and 198 coarse wool sheep narrowed this association to one C/A missense variant of KRT74 (OAR3:133,486,008, P = 1.02E-67) and one T/C SNP in the regulatory region upstream of EDAR (OAR3:61,927,840, P = 2.50E-43). Cellular over-expression and ovine skin section staining assays confirmed that C-KRT74 activated the KRT74 protein and specifically enlarged cell size at the Huxley's layer of the inner root sheath (P < 0.01). This structure enhancement shapes the growing hair shaft into the finer wool than the wild type. Luciferase assays validated that the C-to-T mutation upregulated EDAR mRNA expression via a newly created SOX2 binding site and potentially led to the formation of more hair placodes.

CONCLUSIONS

Two functional mutations driving finer and denser wool production were characterized and offered new targets for genetic breeding during wool sheep selection. This study not only provides a theoretical basis for future selection of fine wool sheep breeds but also contributes to improving the value of wool commodities.

Whole genome sequencing identified genomic diversity and candidated genes associated with economic traits in Northeasern Merino in China

Introduction: Northeast Merino (NMS) is a breed developed in Northeast China during the 1960s for wool and meat production. It exhibits excellent traits such as high wool yield, superior meat quality, rapid growth rate, robust disease resistance, and adaptability to cold climates. However, no studies have used whole-genome sequencing data to investigate the superior traits of NMS. Methods: In this study, we investigated the population structure, genetic diversity, and selection signals of NMS using whole-genome sequencing data from 20 individuals. Two methods (integrated haplotype score and composite likelihood ratio) were used for selection signal analysis, and the Fixation Index was used to explore the selection signals of NMS and the other two breeds, Mongolian sheep and South African meat Merino. Results: The results showed that NMS had low inbreeding levels, high genomic diversity, and a pedigree of both Merino breeds and Chinese local breeds. A total length of 14.09 Mb genomic region containing 287 genes was detected using the two methods. Further exploration of the functions of these genes revealed that they are mainly concentrated in wool production performance (IRF2BP2, MAP3K7, and WNT3), meat production performance (NDUFA9, SETBP1, ZBTB38, and FTO), cold resistance (DNAJC13, LPGAT1, and PRDM16), and immune response (PRDM2, GALNT8, and HCAR2). The selection signals of NMS and the other two breeds annotated 87 and 23 genes, respectively. These genes were also mainly focused on wool and meat production performance. Conclusion: These results provide a basis for further breeding improvement, comprehensive use of this breed, and a reference for research on other breeds.

Modulation of tooth regeneration through opposing responses to Wnt and BMP signals in teleosts

Most vertebrate species undergo tooth replacement throughout adult life. This process is marked by the shedding of existing teeth and the regeneration of tooth organs. However, little is known about the genetic circuitry regulating tooth replacement. Here, we tested whether fish orthologs of genes known to regulate mammalian hair regeneration have effects on tooth replacement. Using two fish species that demonstrate distinct modes of tooth regeneration, threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), we found that transgenic overexpression of four different genes changed tooth replacement rates in the direction predicted by a hair regeneration model: Wnt10a and Grem2a increased tooth replacement rate, whereas Bmp6 and Dkk2 strongly inhibited tooth formation. Thus, similar to known roles in hair regeneration, Wnt and BMP signals promote and inhibit regeneration, respectively. Regulation of total tooth number was separable from regulation of replacement rates. RNA sequencing of stickleback dental tissue showed that Bmp6 overexpression resulted in an upregulation of Wnt inhibitors. Together, these data support a model in which different epithelial organs, such as teeth and hair, share genetic circuitry driving organ regeneration.

Association analysis for SNPs of LIPE and ITGB4 genes with cashmere production performance, body measurement traits and milk production traits in Liaoning cashmere goats

Liaoning cashmere goat (LCG) is one of the excellent cashmere goat breeds in China. Because of its larger size, better cashmere, and better cashmere production performance, people pay special attention to it. This article mainly studied the relationship between SNP loci of LIPE gene and ITGB4 gene and milk production, cashmere production and body measurement traits of LCGs. We further identified potential SNP loci by PCR-Seq polymorphism detection and gene sequence comparison of LIPE and ITGB4 genes. Further, we use SPSS and SHEsis software to analyze their relationship to production performance. The consequence indicated that CC genotype of LIPE gene T16409C locus was dominant genotype in milk production and cashmere production, while CT genotype of LIPE gene T16409C locus was dominant in body size. The CT genotype of C168T locus of ITGB4 gene is the dominant genotype of body type and cashmere production, while the dominant genotype of milk production is TT genotype. Through joint analysis, in haploid combinations, H1H2:CCCT is the dominant haplotype combination in cashmere fineness. H3H4:TTCT is a dominant haplotype combination of milk production traits and body measurement traits. These dominant genotypes can provide a reliable basis for the study of production performance of LCG.

Proteome Analysis of Alpine Merino Sheep Skin Reveals New Insights into the Mechanisms Involved in Regulating Wool Fiber Diameter

Wool fiber is a textile material that is highly valued based on its diameter, which is crucial in determining its economic value. To analyze the molecular mechanisms regulating wool fiber diameter, we used a Data-independent acquisition-based quantitative proteomics approach to analyze the skin proteome of Alpine Merino sheep with four fiber diameter ranges. From three contrasts of defined groups, we identified 275, 229, and 190 differentially expressed proteins (DEPs). Further analysis using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways revealed that pathways associated with cyclic adenosine monophosphate and peroxisome proliferator-activated receptor signaling are relevant to wool fiber diameter. Using the K-means method, we investigated the DEP expression patterns across wool diameter ranges. Using weighted gene co-expression network analysis, we identified seven key proteins (CIDEA, CRYM, MLX, TPST2, GPD1, GOPC, and CAMK2G) that may be involved in regulating wool fiber diameter. Our findings provide a theoretical foundation for identifying DEPs and pathways associated with wool fiber diameter in Alpine Merino sheep to enable a better understanding of the molecular mechanisms underlying the genetic regulation of wool fiber quality.

Construction of a neural network diagnostic model for renal fibrosis and investigation of immune infiltration characteristics

Background

Recently, the incidence rate of renal fibrosis has been increasing worldwide, greatly increasing the burden on society. However, the diagnostic and therapeutic tools available for the disease are insufficient, necessitating the screening of potential biomarkers to predict renal fibrosis.

Methods

Using the Gene Expression Omnibus (GEO) database, we obtained two gene array datasets (GSE76882 and GSE22459) from patients with renal fibrosis and healthy individuals. We identified differentially expressed genes (DEGs) between renal fibrosis and normal tissues and analyzed possible diagnostic biomarkers using machine learning. The diagnostic effect of the candidate markers was evaluated using receiver operating characteristic (ROC) curves and verified their expression using Reverse transcription quantitative polymerase chain reaction (RT-qPCR). The CIBERSORT algorithm was used to determine the proportions of 22 types of immune cells in patients with renal fibrosis, and the correlation between biomarker expression and the proportion of immune cells was studied. Finally, we developed an artificial neural network model of renal fibrosis.

Results

Four candidate genes namely DOCK2, SLC1A3, SOX9 and TARP were identified as biomarkers of renal fibrosis, with the area under the ROC curve (AUC) values higher than 0.75. Next, we verified the expression of these genes by RT-qPCR. Subsequently, we revealed the potential disorder of immune cells in the renal fibrosis group through CIBERSORT analysis and found that immune cells were highly correlated with the expression of candidate markers.

Conclusion

DOCK2, SLC1A3, SOX9, and TARP were identified as potential diagnostic genes for renal fibrosis, and the most relevant immune cells were identified. Our findings provide potential biomarkers for the diagnosis of renal fibrosis.

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.

Genomics of adaptive evolution in the woolly mammoth

Ancient genomes provide a tool to investigate the genetic basis of adaptations in extinct organisms. However, the identification of species-specific fixed genetic variants requires the analysis of genomes from multiple individuals. Moreover, the long-term scale of adaptive evolution coupled with the short-term nature of traditional time series data has made it difficult to assess when different adaptations evolved. Here, we analyze 23 woolly mammoth genomes, including one of the oldest known specimens at 700,000 years old, to identify fixed derived non-synonymous mutations unique to the species and to obtain estimates of when these mutations evolved. We find that at the time of its origin, the woolly mammoth had already acquired a broad spectrum of positively selected genes, including ones associated with hair and skin development, fat storage and metabolism, and immune system function. Our results also suggest that these phenotypes continued to evolve during the last 700,000 years, but through positive selection on different sets of genes. Finally, we also identify additional genes that underwent comparatively recent positive selection, including multiple genes related to skeletal morphology and body size, as well as one gene that may have contributed to the small ear size in Late Quaternary woolly mammoths.

Integrated analysis of miRNAs and mRNA profiling reveals the potential roles of miRNAs in sheep hair follicle development

BACKGROUND

Merino sheep exhibit high wool production and excellent wool quality. The fleece of Merino sheep is predominantly composed of wool fibers grown from hair follicles (HFs). The HF is a complex biological system involved in a dynamic process governed by gene regulation, and gene expression is regulated by microRNAs (miRNAs). miRNA inhibits posttranscriptional gene expression by specifically binding to target messenger RNA (mRNA) and plays an important role in regulating gene expression, the cell cycle and biological development sequences. The purpose of this study was to examine mRNA and miRNA binding to identify key miRNAs and target genes related to HF development. This will provide new and important insights into fundamental mechanisms that regulate cellular activity and cell fate decisions within and outside of the skin.

RESULTS

We analyzed miRNA data in skin tissues collected from 18 Merino sheep on four embryonic days (E65, E85, E105 and E135) and two postnatal days (D7 and D30) and identified 87 differentially expressed miRNAs (DE-miRNAs). These six stages were further divided into two longer developmental stages based on heatmap cluster analysis, and the results showed that DE-mRNAs in Stage A were closely related to HF morphogenesis. A coanalysis of Stage A DE-mRNAs and DE-miRNAs revealed that 9 DE-miRNAs and 17 DE-mRNAs presented targeting relationships in Stage A. We found that miR-23b and miR-133 could target and regulate ACVR1B and WNT10A. In dermal fibroblasts, the overexpression of miR-133 significantly reduced the mRNA and protein expression levels of ACVR1B. The overexpression of miR-23b significantly reduced the mRNA and protein expression levels of WNT10A.

CONCLUSION

This study provides a new reference for understanding the molecular basis of HF development and lays a foundation for further improving sheep HF breeding. miRNAs and target genes related to hair follicular development were found, which provided a theoretical basis for molecular breeding for the culture of fine-wool sheep.