The Tibetan sheep is an indigenous species of the Tibetan plateau and has been well adapted to high-altitude hypoxia. In comparison to lowland sheep breeds, the blood gas indicators have changed and the HIFs signaling pathway is activated in Tibetan sheep. These phenotypic and genetic alterations in Tibetan sheep are thought to be an important basis for adaptation to high-altitude hypoxia and variation in genes encoding the subunits that make up HIFs, such as HIF-1α can affect blood gas indicators. In this study, exons 9, 10, 12 of the HIF-1α gene were sequenced to find variations and 3 SNPs were detected, and these 3 SNPs were genotyped by KASP in 341 Hu sheep and 391 Tibetan sheep. In addition, 197 Hu sheep, 160 Tibetan sheep and 12 Gansu alpine merino sheep were used for blood gas indicators analysis. The results showed significant differences between the blood gas indicators of high-altitude breeds (Tibetan sheep and Gansu alpine merino sheep) and low-altitude breeds (Hu sheep), implying that the differences in blood gas indicators are mainly caused by differences in altitude. The haplotype combinations H2H3 and H1H3 were more frequent in the Tibetan sheep population, H2H3 increases O₂ carrying capacity by increasing hematocrit and hemoglobin concentrations; H1H3 makes O₂ dissociate more readily from oxyhemoglobin by decreasing partial pressure of oxygen and oxygen saturation. These results suggest that variants at the HIF-1α promote the ability of oxygen utilization in Tibetan sheep, which may underpin the survival and reproduction of Tibetan sheep on the Tibetan plateau.

Lymphatic enhancer factor 1 (Lef1) and distal-less homeobox 3 (Dlx3) are the transcription factors involved in regulating hair follicle development in mice, goats, and other animals. Their deletion can lead to hair follicle deficiency. In this study, hematoxylin–eosin staining (HE), real-time quantitative PCR (RT-qPCR), immunohistochemistry, and immunofluorescence were used to analyze the expression, location, and biological functions of Lef1 and Dlx3 in the lateral skin of Gansu Alpine Merino aged 1, 30, 60, and 90 days. The results revealed that the number of hair follicles decreased with age and was significantly higher at 1 day than in the other three age groups (p < 0.05). The mRNA levels of Lef1 and Dlx3 in the skin of 30-day old Gansu Alpine Merino were significantly higher than those in the other three age groups (p < 0.05). Protein expression of Lef1 and Dlx3 was lowest at 1 day (p < 0.05) and peaked at 60 days. Lef1 and Dlx3 exhibited a high density and strong positive expression in the dermal papillae; additionally, Dlx3 exhibited a high density and strong positive expression in the inner and outer root sheaths. Collectively, Lef1 and Dlx3 may facilitate the maturation of secondary hair follicles, which is mainly achieved through the dermal papillae and inner and outer root sheaths.

The Tibetan sheep has an intricate mechanism of adaptation to low oxygen levels, which is influenced by both genetic and environmental factors. The heart plays a crucial role in the adaptation of Tibetan sheep to hypoxia. In the present study, we utilized transcriptomic and proteomic technologies to comprehensively analyze and identify the long non-coding RNAs (lncRNAs), genes, proteins, pathways, and gene ontology (GO) terms associated with hypoxic adaptation in Tibetan sheep at three different altitudes (2500 m, 3500 m, and 4500 m). By integrating the differentially expressed (DE) lncRNA target genes, differentially expressed proteins (DEPs), and differentially expressed genes (DEGs), we were able to identify and characterize the mechanisms underlying hypoxic adaptation in Tibetan sheep. Through this integration, we identified 41 shared genes/proteins, and functional enrichment analyses revealed their close association with lipid metabolism, glycolysis/gluconeogenesis, and angiogenesis. Additionally, significant enrichment was observed in important pathways such as the PPAR signaling pathway, glycolysis/gluconeogenesis, the oxoacid metabolic process, and angiogenesis. Furthermore, the co-expression network of lncRNAs and mRNAs demonstrated that lncRNAs (MSTRG.4748.1, ENSOART00020025894, and ENSOART00020036371) may play a pivotal role in the adaptation of Tibetan sheep to the hypoxic conditions of the plateau. In conclusion, this study expands the existing database of lncRNAs and proteins in Tibetan sheep, and these findings may serve as a reference for the prevention of altitude sickness in humans.

Keratin-related proteins (KAPs) are structural components of wool fibers and are thought to play a key role in regulating the physical and mechanical properties of fibers. Among all KAP genes (KRTAPs), KRTAP6 gene family (KRTAP6-1, KRTAP6-2, KRTAP6-3, KRTAP6-4, and KRTAP6-5) is a very important member with high polymorphism and notable association with some wool traits. In this study, we used real-time fluorescence quantitative PCR (RT-qPCR) and in situ hybridization to investigate spatiotemporal expression of KRTAP6s. The results revealed that KRTAP6 family genes were significantly expressed during anagen compared to other stages (p < 0.05). And it was found the five genes were expressed predominantly in the dermal papillae, inner and outer root sheaths, and showed a distinct spatiotemporal expression pattern. Also, it was found that KRTAP6-1 and KRTAP6-5 mRNA expression was negatively correlated with wool mean fiber diameter (MFD) and mean staple strength (MSS) (p < 0.05). In summary, the KRTAP6 family genes share a similar spatiotemporal expression pattern. And KRTAP6-1 and KRTAP6-5 may regulate the MFD and MSS of Gansu Alpine fine-wool sheep wool by changing the expression.

Changes in keratin gene expression and spatiotemporal regulation determine the compositional content and cellular localization of wool keratin, thereby affecting wool traits. Therefore, keratin gene family member 32 (KRT32) was selected for a study using RT-qPCR, immunofluorescence, and penta-primer amplification refractory mutation system (PARMS) techniques. The results showed that KRT32 mRNA was highly expressed in the skin and localized to the inner root sheath (IRS), outer root sheath (ORS) and dermal papilla (DP). Sequencing results identified eight SNPs in KRT32, and association analyses revealed that the variations were significantly associated with multiple traits in wool (p < 0.05), including MFD, CF and MFC. The constructed haplotype combination H₂H₃ has higher CF and smaller MFD than other haplotype combination (p < 0.05). In conclusion, KRT32 can be used as a candidate gene for molecular genetic improvement of wool in Gansu Alpine Fine-wool sheep.

Homobox C13 (Hoxc13) is an important transcription factor in hair follicle cycle development, and its deletion had been found in a variety of animals leading to abnormal hair growth and disruption of the hair follicle system. In this study, we used immunofluorescence, immunohistochemistry, real-time fluorescence quantitative PCR (RT-qPCR), and Kompetitive Allele-Specific PCR (KASP) genotyping to investigate molecular genetic characteristics of the Hoxc13 gene in Gansu alpine fine-wool sheep. The results revealed that Hoxc13 was significantly expressed during both the anagen and catagen phases (p < 0.05). It was found to be highly expressed predominantly in the dermal papillae and the inner and outer root sheaths, showing a distinct spatiotemporal expression pattern. Two single nucleotide polymorphisms (SNPs) in the exon 1 of Hoxc13, both the individual locus genotypes and the combined haplotypes were found to be correlated with wool length (p < 0.05). It was determined the mutations led to changes in mRNA expression, in which higher expression of this gene was related with longer wool length. In summary, this unique spatiotemporal expression pattern of the Hoxc13 gene may regulate the wool length of Gansu alpine fine-wool sheep, which can be used as a molecular genetic marker for wool traits and thus improve the breed.