Genome-Wide Selective Signatures Reveal Candidate Genes Associated with Hair Follicle Development and Wool Shedding in Sheep

Transcriptomic analysis of regulatory mechanisms in the telogen-anagen transition of ovine hair follicles

BACKGROUND

Dorper sheep are celebrated for their fast maturation and superior meat quality, with some shedding their wool each spring. Wool shedding occurs naturally due to the hair follicle (HF) cycle, but its regulatory mechanisms remain unclear and need further investigation.

RESULTS

In this study, shedding and non-shedding sheep were selected from the same Dorper flock. Skin samples were collected in September of the first year and January and March of the following years. RNA sequencing was performed on these samples. Principal component analysis (PCA) was used to assess the results. A total of 2536 differentially expressed genes (DEGs) were identified. Using a clustering heatmap and fuzzy clustering analysis three distinct gene expression patterns were identified: A pattern (high expression in anagen), T1 pattern, and T2 pattern (high expression in telogen). For each pattern, differentially expressed genes (DEGs) were analyzed through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Combining this with pathway expression analysis, six A-pattern and fourteen T-pattern pathways linked to telogen-anagen transition in the HF cycle were identified. Networks of key pathways were then constructed. Additionally, key genes were identified in the telogen-anagen transition, including one A-pattern gene and seven T-pattern (T1, 1; T2, 6) genes, using the Maximal Clique Centrality (MCC) tool in Cytoscape. Predicted transcription factors (TFs) involved in key pathways, such as LEF and STAT5B, were identified. Finally, RNA-seq results were confirmed by RT-qPCR.

CONCLUSION

This study highlights critical genes and pathways in the telogen-anagen transition, and transcriptome sequencing along with bioinformatics analysis provides new insights into the regulatory mechanisms of the HF cycle and development.

Genetic Diversity and Selection Signal Analysis of Hu Sheep Based on SNP50K BeadChip

This research is designed to examine the genetic diversity and kinship among Hu sheep, as well as to discover genes associated with crucial economic traits. A selection of 50 unrelated adult male Hu sheep underwent genotyping with the SNP50K BeadChip. Seven indicators of genetic diversity were assessed based on high-quality SNP data: effective population size (Ne), polymorphic information content (PIC), polymorphic marker ratio (PN), expected heterozygosity (He), observed heterozygosity (Ho), effective number of alleles, and minor allele frequency (MAF). Plink software was employed to compute the IBS genetic distance matrix and detect runs of homozygosity (ROHs), while the G matrix and principal component analysis were performed using GCTA software. Selective sweep analysis was carried out using ROH, Pi, and Tajima's D methodologies. This study identified a total of 64,734 SNPs, of which 56,522 SNPs remained for downstream analysis after quality control. The population displayed relatively high genetic diversity. The 50 Hu sheep were ultimately grouped into 12 distinct families, with families 6, 8, and 10 having the highest numbers of individuals, each consisting of 6 sheep. Furthermore, a total of 294 ROHs were detected, with the majority having lengths between 1 and 5 Mb, and the inbreeding coefficient FROH was 0.01. In addition, 41, 440, and 994 candidate genes were identified by ROH, Pi, and Tajima's D methods, respectively, with 3 genes overlapping (BMPR1B, KCNIP4, and FAM13A). These results offer valuable insights for future Hu sheep breeding, genetic assessment, and population management.

Population structure and selective signature of Kirghiz sheep by Illumina Ovine SNP50 BeadChip

Objective

By assessing the genetic diversity and associated selective traits of Kirghiz sheep (KIR), we aim to uncover the mechanisms that contribute to sheep's adaptability to the Pamir Plateau environment.

Methods

This study utilized Illumina Ovine SNP50 BeadChip data from KIR residing in the Pamir Plateau, Qira Black sheep (QBS) inhabiting the Taklamakan Desert, and commonly introduced breeds including Dorper sheep (DOR), Suffolk sheep (SUF), and Hu sheep (HU). The data was analyzed using principal component analysis, phylogenetic analysis, population admixture analysis, kinship matrix analysis, linkage disequilibrium analysis, and selective signature analysis. We employed four methods for selective signature analysis: fixation index (Fst), cross-population extended homozygosity (XP-EHH), integrated haplotype score (iHS), and nucleotide diversity (Pi). These methods aim to uncover the genetic mechanisms underlying the germplasm resources of Kirghiz sheep, enhance their production traits, and explore their adaptation to challenging environmental conditions.

Results

The test results unveiled potential selective signals associated with adaptive traits and growth characteristics in sheep under harsh environmental conditions, and annotated the corresponding genes accordingly. These genes encompass various functionalities such as adaptations associated with plateau, cold, and arid environment (ETAA1, UBE3D, TLE4, NXPH1, MAT2B, PPARGC1A, VEGFA, TBX15 and PLXNA4), wool traits (LMO3, TRPS1, EPHA5), body size traits (PLXNA2, EFNA5), reproductive traits (PPP3CA, PDHA2, NTRK2), and immunity (GATA3).

Conclusion

Our study identified candidate genes associated with the production traits and adaptation to the harsh environment of the Pamir Plateau in Kirghiz sheep. These findings provide valuable resources for local sheep breeding programs. The objective of this study is to offer valuable insights for the sustainable development of the Kirghiz sheep industry.

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.

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.

Comprehensive analysis of the expression profiles of mRNA, lncRNA, circRNA, and miRNA in primary hair follicles of coarse sheep fetal skin

BACKGROUND

The Qinghai Tibetan sheep, a local breed renowned for its long hair, has experienced significant deterioration in wool characteristics due to the absence of systematic breeding practices. Therefore, it is imperative to investigate the molecular mechanisms underlying follicle development in order to genetically enhance wool-related traits and safeguard the sustainable utilization of valuable germplasm resources. However, our understanding of the regulatory roles played by coding and non-coding RNAs in hair follicle development remains largely elusive.

RESULTS

A total of 20,874 mRNAs, 25,831 circRNAs, 4087 lncRNAs, and 794 miRNAs were annotated. Among them, we identified 58 DE lncRNAs, 325 DE circRNAs, 924 DE mRNAs, and 228 DE miRNAs during the development of medullary primary hair follicle development. GO and KEGG functional enrichment analyses revealed that the JAK-STAT, TGF-β, Hedgehog, PPAR, cGMP-PKG signaling pathway play crucial roles in regulating fibroblast and epithelial development during skin and hair follicle induction. Furthermore, the interactive network analysis additionally identified several crucial mRNA, circRNA, and lncRNA molecules associated with the process of primary hair follicle development. Ultimately, by investigating DEmir's role in the ceRNA regulatory network mechanism, we identified 113 circRNA-miRNA pairs and 14 miRNA-mRNA pairs, including IGF2BP1-miR-23-x-novel-circ-01998-MSTRG.7111.3, DPT-miR-370-y-novel-circ-005802-MSTRG.14857.1 and TSPEAR-oar-miR-370-3p-novel-circ-005802- MSTRG.10527.1.

CONCLUSIONS

Our study offers novel insights into the distinct expression patterns of various transcription types during hair follicle morphogenesis, establishing a solid foundation for unraveling the molecular mechanisms that drive hair development and providing a scientific basis for selectively breeding desirable wool-related traits in this specific breed.

Screening of hair follicle telogen-associated circRNAs in sheep and construction of their ceRNA network

Sheep breeds with hair-shedding traits have many advantages over non-shedding sheep breeds, not only because of reduced shearing labor and feeding management costs but also because it reduces in vitro parasites and improves adaptability to summer heat stress. The wool of Dorper sheep naturally sheds in spring due to the periodic growth of hair follicles. CircRNAs primarily regulate the morphogenesis of hair follicles through the ceRNA mechanism. In this study, five 2-year-old Dorper ewes with extreme hair-shedding phenotype (S) and three Dorper ewes with non-shedding (N) phenotype were selected for subsequent analyses. For RNA extraction, skin tissues were collected on 27th September 2019 (S1, N1), 3rd January 2020 (S2, N2), and 17th March 2020 (S3, N3), which were then subjected to RNA-seq. RNA-seq technology revealed 20,185 novel circRNAs in the hair follicles of Dorper sheep. Among them, 1450 circRNAs were differentially expressed (DE). Clustering heatmap and expression pattern analyses were performed on DE circRNAs, which indicated 78 circRNAs with T pattern (Telogen, highly expressed in telogen), and the source genes for candidate circRNAs were further screened by functional enrichment analysis, which identified 13 crucial genes enriched in pathways associated with hair follicle development. Additionally, a ceRNA regulatory network comprising 4 circRNAs, 11 miRNAs, and 13 target genes was constructed. Overall, this study screened circRNAs that may be associated with the telogen phase of hair follicles in sheep, providing a relevant theoretical basis for wool shedding in sheep and for breeding Dorper sheep with automatic wool shedding.

Phenotype transition from wild mouflon to domestic sheep

The domestication of animals started around 12,000 years ago in the Near East region. This "endless process" is characterized by the gradual accumulation of changes that progressively marked the genetic, phenotypic and physiological differences between wild and domesticated species. The main distinctive phenotypic characteristics are not all directly attributable to the human-mediated selection of more productive traits. In the last decades, two main hypotheses have been proposed to clarify the emergence of such a set of phenotypic traits across a variety of domestic species. The first hypothesis relates the phenotype of the domesticated species to an altered thyroid hormone-based signaling, whereas the second one relates it to changes in the neural crest cells induced by selection of animals for tameness. These two hypotheses are not necessarily mutually exclusive since they may have contributed differently to the process over time and space. The adaptation model induced by domestication can be adopted to clarify some aspects (that are still controversial and debated) of the long-term evolutionary process leading from the wild Neolithic mouflon to the current domestic sheep. Indeed, sheep are among the earliest animals to have been domesticated by humans, around 12,000 years ago, and since then, they have represented a crucial resource in human history. The aim of this review is to shed light on the molecular mechanisms and the specific genomic variants that underlie the phenotypic variability between sheep and mouflon. In this regard, we carried out a critical review of the most recent studies on the molecular mechanisms that are most accredited to be responsible for coat color and phenotype, tail size and presence of horns. We also highlight that, in such a complicate context, sheep/mouflon hybrids represent a powerful and innovative model for studying the mechanism by which the phenotypic traits related to the phenotypic responses to domestication are inherited. Knowledge of these mechanisms could have a significant impact on the selection of more productive breeds. In fact, as in a journey back in time of animal domestication, the genetic traits of today's domestic species are being progressively and deliberately shaped according to human needs, in a direction opposite to that followed during domestication.

Genome-Wide Association Study of Fiber Diameter in Alpacas

The aim of this study was the identification of candidate genomic regions associated with fiber diameter in alpacas. DNA samples were collected from 1011 female Huacaya alpacas from two geographical Andean regions in Peru (Pasco and Puno), and three alpaca farms within each region. The samples were genotyped using an Affymetrix Custom Alpaca genotyping array containing 76,508 SNPs. After the quality controls, 960 samples and 51,742 SNPs were retained. Three association study methodologies were performed. The GWAS based on a linear model allowed us to identify 11 and 35 SNPs (-log10(p-values) > 4) using information on all alpacas and alpacas with extreme values of fiber diameter, respectively. The haplotype and marker analysis method allowed us to identify nine haplotypes with standardized haplotype heritability higher than six standard deviations. The selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) allowed us to identify 180 SNPs with XP-EHH values greater than |3|. Four candidate regions with adjacent SNPs identified via two association methods of analysis are located on VPA6, VPA9, VPA29 and one chromosomally unassigned scaffold. This study represents the first analysis of alpaca whole genome association with fiber diameter, using a recently assembled alpaca SNP microarray.

Detection of Candidate Genes Associated with Fecundity through Genome-Wide Selection Signatures of Katahdin Ewes

One of the strategies to genetically improve reproductive traits, despite their low inheritability, has been the identification of candidate genes. Therefore, the objective of this study was to detect candidate genes associated with fecundity through the fixation index (F_ST) and runs of homozygosity (ROH) of selection signatures in Katahdin ewes. Productive and reproductive records from three years were used and the genotypes (OvineSNP50K) of 48 Katahdin ewes. Two groups of ewes were identified to carry out the genetic comparison: with high fecundity (1.3 ± 0.03) and with low fecundity (1.1 ± 0.06). This study shows for the first time evidence of the influence of the CNOT11, GLUD1, GRID1, MAPK8, and CCL28 genes in the fecundity of Katahdin ewes; in addition, new candidate genes were detected for fecundity that were not reported previously in ewes but that were detected for other species: ANK2 (sow), ARHGAP22 (cow and buffalo cow), GHITM (cow), HERC6 (cow), DPF2 (cow), and TRNAC-GCA (buffalo cow, bull). These new candidate genes in ewes seem to have a high expression in reproduction. Therefore, future studies are needed focused on describing the physiological basis of changes in the reproductive behavior influenced by these genes.

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.