Integrative ATAC-seq and RNA-seq analysis of myogenic differentiation of ovine skeletal muscle satellite cell.

Skeletal muscle satellite cells (SMSCs) play an important role in regulating muscle growth and regeneration. Chromatin accessibility allows physical interactions that synergistically regulate gene expression through enhancers, promoters, insulators, and chromatin binding factors. However, the chromatin accessibility altas and its regulatory role in ovine myoblast differentiation is still unclear. Therefore, ATAC-seq and RNA-seq analysis were performed on ovine SMSCs at the proliferation stage (SCG) and differentiation stage (SCD). 17,460 DARs (differential accessibility regions) and 3732 DEGs (differentially expressed genes) were identified. Based on joint analysis of ATAC-seq and RNA-seq, we revealed that PI3K-Akt, TGF-ß and other signaling pathways regulated SMSCs differentiation. We identified two novel candidate genes, FZD5 and MAP2K6, which may affect the proliferation and differentiation of SMSCs. Our data identify potential cis regulatory elements of ovine SMSCs. This study can provide a reference for exploring the mechanisms of the differentiation and regeneration of SMSCs in the future.

Multi-order orbital angular momentum mode generators based on integrated long-period fiber gratings.

We propose integrated long-period fiber gratings (LPFGs) fabricated by a CO2 laser to realize a multi-channel and multi-order orbital angular momentum (OAM) mode generator. The integrated LPFG is inscribed on multiple surfaces of the few-mode fiber (FMF) by rotating the fiber in the opposite direction at an angle θ. By controlling the rotation angle, the number of integrated LPFGs can be set. The selected rotation angle is 43 ∘, which can integrate up to nine LPFGs, i.e., realizing that the number of channels for first-order orbital angular momentum (OAM) mode conversion is nine. The integrated LPFGs fabricated in this method allow a flexible design of channel spacing. In addition, the flexible selection of the integrated grating period achieves the simultaneous generation of multi-channel second-order and third-order OAM mode conversion. The multi-channel and multi-order OAM mode generators have important application in optical communication multiplexing systems and OAM sensing.

Mode conversion in graded-index few-mode fiber via hollow cylindrical long-period fiber gratings.

We demonstrate the realization of mode conversion using hollow cylindrical long-period fiber gratings (LPFGs) inscribed in graded-index few-mode fibers (FMFs) by a femtosecond laser. By precisely shaping the refractive index modulation into a hollow cylindrical structure, we enable efficient coupling from the fundamental mode to high-order modes (LP11, LP02, LP21, LP12, LP31, LP03, and LP22 modes). The method achieves high-efficiency and low-loss mode conversion across various mode groups, marking a first for graded-index FMFs with different grating periods. The influence of the hollow cylindrical LPFG radius on mode conversion efficiency and spectral characteristics is thoroughly investigated. Additionally, incorporating a linear polarizer allowed for distinguishing between LP modes within a mode group, enhancing the tunability of the mode converter. These mode conversion devices have significant potential for applications in mode gain equalization and mode scrambling for mode division multiplexing (MDM) systems.

Porcine ZBED6 regulates growth of skeletal muscle and internal organs via multiple targets.

ZBED6 (zinc finger BED domain containing protein 6) is a transcription factor unique to placental mammals and its interaction with the IGF2 (insulin-like growth factor 2) locus plays a prominent role in the regulation of postnatal skeletal muscle growth. Here, we generated lean Bama miniature pigs by generating ZBED6-knockout (ZBED6-/-) and investigated the mechanism underlying ZBED6 in growth of muscle and internal organs of placental mammals. ZBED6-/- pigs show markedly higher lean mass, lean mass rate, larger muscle fiber area and heavier internal organs (heart and liver) than wild-type (WT) pigs. The striking phenotypic changes of ZBED6-/- pigs coincided with remarkable upregulation of IGF2 mRNA and protein expression across three tissues (gastrocnemius muscle, longissimus dorsi, heart). Despite a significant increase in liver weight, ZBED6-/- pigs show comparable levels of IGF2 expression to those of WT controls. A mechanistic study revealed that elevated methylation in the liver abrogates ZBED6 binding at the IGF2 locus, explaining the unaltered hepatic IGF2 expression in ZBED6-/- pigs. These results indicate that a ZBED6-IGF2-independent regulatory pathway exists in the liver. Transcriptome analysis and ChIP-PCR revealed new ZBED6 target genes other than IGF2, including cyclin dependent kinase inhibitor 1A (CDKN1A) and tsukushi, small leucine rich proteoglycan (TSKU), that regulates growth of muscle and liver, respectively.

Role of Csdc2 in Regulating Secondary Hair Follicle Growth in Cashmere Goats.

Cashmere goats possess two types of hair follicles, with the secondary hair follicles producing valuable cashmere fiber used for textiles. The growth of cashmere exhibits a seasonal pattern arising from photoperiod change. Transcription factors play crucial roles during this process. The transcription factor, cold-shock domain, containing C2 (Csdc2) plays a crucial role in modulating cell proliferation and differentiation. Our preceding research indicated that the expression of Csdc2 changes periodically during anagen to telogen. However, the mechanisms of Csdc2 in regulating SHF growth remain unclear. Here, we found that the knockdown of Csdc2 inhibits the proliferation of dermal papilla cells. ChIP-Seq analysis showed that Csdc2 had a unique DNA binding motif in SHFs. Through conjoint analysis of ChIP-Seq and RNA-Seq, we revealed a total of 25 candidate target genes of Csdc2. Notably, we discovered a putative Csdc2 binding site within roundabout guidance receptor 2 (Robo2) on chromosome 1 of the goat genome. Furthermore, qRT-PCR and dual-luciferase reporter assay confirmed Csdc2's positive regulatory influence on Robo2. These findings expand the research field of hair follicle transcriptional regulatory networks, offering insights into molecular breeding strategies to enhance cashmere production in goats.

Broadband linearly polarized mode converter based on over-coupled long-period fiber grating.

We demonstrate the fabrication of over-coupled long-period fiber gratings (LPFGs) in the 1.55-µm and 2-µm wavebands enabling broadband linearly polarized LP11 mode conversion using a CO2 laser. The birefringence of the fiber is caused by on one side laser exposure and increases with the increase of refractive index modulation depth, which realizes the conversion of linearly polarized modes. The mode conversion bandwidth can be significantly increased by using the over-coupled LPFG. The 10-dB bandwidth of the LPFGs with |κ|L values of π/2, 3π/2, and 5π/2 are 33.04, 80.84, and 114.08 nm at 1.55 µm waveband, respectively. The maximum bandwidth of the over-coupled LPFG is 3.79 times higher than that of conventional LPFG. The operating wavelength of the mode converter can be extended to 2.0 µm wavebands and the maximum 10-dB bandwidth reaches 161.32 nm. The proposed broadband linearly polarized mode converters could have potential application in the fields of mode division multiplexing systems, fiber laser systems.

Ultrasensitive magnetic field sensor based on cladding-etched long-period grating.

We demonstrate a high-sensitivity fiber-optic magnetic field sensor, which consists of a cladding-etched long-period fiber grating (LPFG) near the dispersion turning point (DTP) integrated with a magnetic fluid (MF). By reducing the cladding diameter of the LPFG, the fundamental mode is coupled to the lowest order cladding mode (LP0,2) near the DTP, which has a much higher surrounding refractive index sensitivity. Thanks to the excellent magneto-optical characteristics of the MF, the proposed sensor can achieve a magnetic field intensity sensitivity of 44.69 nm/mT in the range of 3-7.4 mT. The minimum magnetic field intensity that can be detected is 0.45 µT due to the 0.02-nm wavelength resolution of the optical spectrum analyzer. The proposed etched DTP-LPFG-based sensor with ultrahigh magnetic field sensitivity could have potential applications in magnetic fields and electrical systems.

Multiple core modes conversion using helical long-period fiber gratings.

We propose and demonstrate the fabrication of an all-fiber mode converter enabling simultaneous generation of multiple high-order core modes, which is realized by inscribing a helical long-period grating (HLPG) in a few-mode fiber (FMF) using a femtosecond laser. Helical refractive index modulation is introduced by continuously irradiating the core region with a highly focused femtosecond laser, while the fiber moves in a spiral path through a three-dimensional translation stage. Mode conversion from the LP01 mode to high-order core modes, including LP11, LP21, LP31, LP02, LP12, and LP41 modes, is achieved by controlling the inscription pitch of the grating. Moreover, first-, second-, third-, and fourth-order orbital angular momentum (OAM) modes can be directly generated using the HLPGs, and multiple OAM modes of different topological charges can be simultaneously excited using a single high diffraction order HLPG. This approach offers a new option for implementing with high-integration high-order mode converters or OAM mode generators.

Whole-genome identification of transposable elements reveals the equine repetitive element insertion polymorphism in Chinese horses.

Transposable elements (TEs) are diverse, abundant, and complicated in genomes. They not only can drive the genome evolution process but can also act as special resources for adaptation. However, little is known about the evolutionary processes that shaped horses. In this work, 126 horse assemblages involved in most horse breeds in China were used to investigate the patterns of TE variation for the first time. By using RepeatMasker and melt software, we found that the horse-specific short interspersed repetitive elements family, equine repetitive elements (ERE1), exhibited polymorphisms in horse genomes. Phylogenetic analysis based on these ERE1 loci (minor allele frequency ≥0.05) revealed three major horse groups, namely, those in northern China, southern China, and Qinghai-Tibetan, which mirrors the result determined by SNPs to some extent. The present ERE1 family emerged ~0.26 to 1.77 Mya ago, with an activity peak at ~0.49 Mya, which matches the early stage of the horse lineage and decreases after the divergence of Equus caballus and Equus ferus przewalskii. To detect the functional ERE1(s) associated with adaptation, locus-specific branch length, genome-wide association study, and absolute allele frequency difference analyses were conducted and resulted in two common protein-coding genes annotated by candidate ERE1s. They were clustered into the vascular smooth muscle contraction (p = 0.01, EDNRA) and apelin signalling pathways (p = 0.02, NRF1). Notably, ERE1 insertion into the EDNRA gene showed a higher association with adaptation among southern China horses and other horses in 15 populations and 451 individuals (p = 4.55 e-8). Our results provide a comprehensive understanding of TE variations to analyse the phylogenetic relationships and traits relevant to adaptive evolution in horses.

Comparative analysis of the genetic diversity of the neutral microsatellite loci and second exon of the goat MHC-DQB1 gene.

This study compared and analyzed the genetic diversity and population structure of exon 2 of the DQB1 gene and 13 autosomal neutral microsatellite markers from 14 Chinese goat breeds to explore the potential evolutionary mechanism of the major histocompatibility complex (MHC). A total of 287 haplotypes were constructed from MHC-DQB1 exon 2 from 14 populations, and 82 nucleotide polymorphic sites (SNPs, 31.78%) and 172 heterozygous individuals (79.12%) were identified. The FST values of the microsatellites and MHC-DQB ranged between 0.01831-0.26907 and 0.00892-0.38871, respectively. Furthermore, 14 goat populations showed rich genetic diversity in the microsatellite loci and MHC-DQB1 exon 2. However, the population structure and phylogenetic relationship represented by the two markers were different. Positive selection and Tajima's D test results showed the occurrence of a diversified selection mechanism, which was primarily based on a positive and balancing selection in goat DQB. This study also found that the DQB sequences of bovines exhibited trans-species polymorphism (TSP) among species and families. In brief, this study indicated that positive and balancing selection played a major role in maintaining the genetic diversity of DQB, and TSP of MHC in bovines was common, which enhanced the understanding of the MHC evolution.

Femtosecond laser inscribed parallel long-period fiber gratings for multi-channel core mode conversion.

We propose and demonstrate the inscription of parallel long-period gratings (LPGs) in a few-mode fiber (FMF) using femtosecond lasers. Mode conversion from the fundamental (LP01) mode to high-order core modes, including LP11, LP21, LP31, LP02, and LP12, is achieved by controlling the inscription period of the gratings. Taking advantage of the highly focused femtosecond laser, LPGs with different off-axis offsets were fabricated, and the resonance wavelength and the inscription efficiency of the gratings versus the offset were investigated. Based on the off-axis writing technique and using the femtosecond laser source, we wrote parallel LPGs that contain multi-gratings in a single FMF and achieved a multi-channel core mode converter in a single FMF with flexibility in terms of the resonant wavelength and mode conversion among different modes. This approach offers a new, to the best of our knowledge, option for implementation with high integration, and a multi-channel mode converter, which could find potential applications in FMF multi-wavelength laser systems, and wavelength/mode division multiplex communication systems. Furthermore, these microstructured LPGs integrated into an optical fiber can be used as a multifunctional sensor.

A study using single-locus and multi-locus genome-wide association study to identify genes associated with teat number in Hu sheep.

The multiple teats trait is common in many species of mammals and is considered related to lactation ability in swine. However, in Hu sheep, related gene research is still relatively limited. In this study, a genome-wide association study was used to identify genetic markers and genes related to the number of teats in the Hu sheep population, a native Chinese sheep breed. A single marker method and several multi-locus methods were utilized. A total of 61 SNPs were found to be related to the number of teats. Among these, 11 SNPs and one SNP were consistently detected by two and three multi-locus models respectively. Four SNPs were concordantly identified between the single marker and multi-locus methods. We also performed quantitative real-time PCR testing of these identified candidate genes, identifying three genes with significantly different expression. Our study suggested that the LHFP, DPYSL2, and TDP-43 genes may be related to the number of teats in sheep. The combination of single and multi-locus GWAS detected additional SNPs not found with only one model. Our results provide new and important insights into the genetic mechanisms of the mammalian multiparous teat phenotype. These findings may be useful for future breeding and understanding the genetics of sheep and other livestock.

A deletion variant within the FGF5 gene in goats is associated with gene expression levels and cashmere growth.

The FGF5 gene has been associated with the regulation of fibre length in mammals, including cashmere goats. A deletion variant at ~14 kb downstream of the FGF5 gene showed significant divergence between cashmere and non-cashmere goats in previous studies. In this study, we designed specific primers to genotype the deletion variant. The results of gel electrophoresis and Sanger sequencing revealed that a 507-bp deletion mutation is located at 95 454 685-95 455 191 of chromosome 6 in goats. Genotyping data from a large panel of 288 goats showed that the deletion at the FGF5 gene locus appeared to be associated with cashmere length. The deletion variant was close to fixation (frequency 0.97) in cashmere goats. Furthermore, electrophoretic mobility shift assays for evaluating DNA-protein interaction and mRNA expression levels of FGF5 suggested that the deletion variant may serve as a cis-acting element by specifically binding transcription factors to mediate quantitative changes in FGF5 mRNA expression. Our study illustrates how a structural mutation of the FGF5 gene has contributed to the cashmere growth phenotype in domestic goats. The deletion mutation within the FGF5 gene could potentially serve as a molecular marker of cashmere growth in cashmere goat breeding.

Genome-Wide Analysis of Circular RNAs Reveals circCHRNG Regulates Sheep Myoblast Proliferation via miR-133/SRF and MEF2A Axis.

As relatively new members of the non-coding RNA family, circRNAs play important roles in a variety of biological processes. However, the temporal expression pattern and the function of circRNAs during sheep skeletal muscle development remains unclear. This study aimed to identify circRNAs related to sheep skeletal muscle development and explore their roles in myoblast proliferation. The circRNA expression profiles of longissimus dorsi of sheep from F90, L30, and A3Y were obtained by the RNA-seq method. The function and mechanisms of the novel circCHRNG in muscle satellite cell proliferation were explored using CCK-8 assay, Western blot, qPCR, and dual-luciferase reporter assay. We identified 12,375 circRNAs, including 476, 133, and 233 DEcircRNAs found among three comparative groups. KEGG results showed that DEcircRNAs were enriched in muscle contraction, the regulation of cell proliferation, and the AMPK, insulin, and PI3K-Akt signaling pathways. Notably, a novel circRNA, termed circRNA CHRNG, acts as a miR-133 sponge to promote skeletal muscle satellite cell proliferation. Our study provides a systematic description of circRNAs of ovine skeletal muscle across fetal, lamb, and adult stages. GO and KEGG analyses showed that DEcircRNAs were enriched in multiple pathways associated with muscle development, such as the PI3K-Akt and AMPK signaling pathways. In addition, we propose that circCHRNG acts as a miR-133 sponge to upregulate the expression levels of SRF and MEF2A, thereby promoting myoblast proliferation.

Comprehensive analysis of lncRNAs involved in skeletal muscle development in ZBED6-knockout Bama pigs.

BACKGROUND: The mutation of insulin-like growth factor 2 (IGF2 mutation) that a single-nucleotide substitution (G→A) in the third intron of IGF2 abrogates the interaction with zinc finger BED-type containing 6 (ZBED6) and leads to increased muscle mass in pigs. IGF2 mutation knock-in (IGF2 KI) and ZBED6 knockout (ZBED6 KO) lead to changes in IGF2 expression and increase muscle mass in mice and pigs. Long noncoding RNAs (lncRNAs) may participate in numerous biological processes, including skeletal muscle development. However, the role of the ZBED6-lncRNA axis in skeletal muscle development is poorly characterized. RESULTS: In this study, we assembled transcriptomes using RNA-seq data published in previous studies by our group and identified 11,408 known lncRNAs and 2269 potential lncRNAs in seven tissues, heart, longissimus dorsi, gastrocnemius muscle, liver, spleen, lung and kidney, of ZBED6 KO (lean mass model) and WT Bama pigs. ZBED6 affected the expression of 1570 lncRNAs (differentially expressed lncRNAs [DE-lncRNAs]; log2-fold change ≥ 1, nominal p-value ≤ 0.05) in the seven examined tissues. The expressed lncRNAs (FPKM > 0.1) exhibited tissue-specific patterns in WT pigs. Specifically, 3410 lncRNAs were expressed exclusively in only one tissue. Potential functions of lncRNAs were indirectly predicted by searching their target cis- and trans-regulated protein-coding genes. LncRNAs with tissue-specific expression influence numerous genes related to tissue functions. Weighted gene coexpression network analysis (WGCNA) of 1570 DE-lncRNAs between WT and ZBED6 KO pigs was used to define the following six lncRNA modules specific to different tissues: skeletal muscle, heart, lung, spleen, kidney and liver modules. Furthermore, by conjoint analysis of longissimus dorsi data (tissue-specific expression, muscle module and DE-lncRNAs) and ChIP-PCR revealed NONSUSG002145.1 (adjusted p-values = 0.044), which is coexpressed with the IGF2 gene and binding with ZBED6, may play important roles in ZBED6 KO pig skeletal muscle development. CONCLUSIONS: These findings indicate that the identified lncRNAs may play essential roles in tissue function and regulate the mechanism of ZBED6 action in skeletal muscle development in pigs. To our knowledge, this is the first study describing lncRNAs in ZBED6 KO pigs. These results may open new research directions leading to a better understanding of the global functions of ZBED6 and of lncRNA functions in skeletal muscle development in pigs.

Wavelength-tunable L-band mode-locked fiber laser using a long-period fiber grating.

We demonstrate an L-band wavelength-tunable passively mode-locked fiber laser using a single long-period fiber grating (LPFG) as a narrow-band optical attenuator (NBOA). Through bending the LPFG, the central wavelength can be continuously tuned from 1582.02 to 1597.29 nm, while the output power only varies from 1.465 to 1.057 mW, approximately a rate of 22 µW/nm variation. This is the first time that LPFG is functioned as a NBOA in mode-locked fiber lasers, showing the great advantage of less impact on output power variation reduction. Besides, the total cavity length is 5.08 m, which is the shortest length yet reported in wavelength-tunable mode-locked fiber lasers. The wavelength tuning could also be realized at harmonic mode locking with tuning range of 14.69 nm under 5th harmonic.

[Design and Clinical Application of Split Memory Alloy Sternum Plate].

OBJECTIVE: The design and development of split memory alloy sternum bone plate are discussed, and the effect of split memory alloy sternum bone plate internal fixation in the treatment of sternal fractures are analysed. METHODS: The structure of the product is designed according to the anatomy and physiological characteristics of human bones, and the cross section shape of the product is designed according to the cross section shape of human bones. Internal fixation is effective in the treatment of sternal fracture. RESULTS: The split memory alloy sternal plate was successfully designed and developed, and all the patients with sternal fractures treated by internal fixation were clinically healed, the hospitalization and fracture healing time were significantly shortened, and no obvious complications occurred. CONCLUSIONS: The application of split memory alloy sternal plate internal fixation in the treatment of sternal fracture has the advantages of small trauma, simple operation, safety, reliable fixation, good histocompatibility and less complications, and is conducive to promoting fracture healing and respiratory function improvement.

Trend analysis of the role of circular RNA in goat skeletal muscle development.

BACKGROUND: Circular RNA (circRNA) is produced during the splicing of mRNA (in addition to linear splicing) and is part of the gene regulatory network. The temporal expression patterns the different developmental stages were inseparable from these molecules' function. RESULTS: Skeletal muscles of Anhui white goat (AWG) across seven fetal to postnatal development stages were sequenced and 21 RNA sequencing libraries were constructed. We thereby identified 9090 circRNAs and analyzed their molecular properties, temporal expression patterns, and potential functions at the different stages. CircRNAs showed complexities and diversity of formation as the same host gene produces multiple isoforms of these nucleic acids with different expression profiles. The differential expression of 2881 circRNAs (DECs, P < 0.05) was identified and four were randomly selected and validated by qPCR. Moreover, 1118 DECs under strict selected (SDECs, |log2FC| > 2 and P-adj value < 0.01) showed 4 expression trends (Clusters 0, 19, 16 and 18). Cluster 0 molecules had increasing expression at all stages with effects on muscle through metabolism, regulation of enzyme activity, and biosynthesis. Cluster 16 circRNAs had high expression in the early and late stages and are involved in "Wnt signaling pathway", "AMPK signaling pathway" and others. Cluster 18 molecules were mainly expressed at F120 and participate in "cytoskeletal protein binding", "Notch signaling pathway" and so on. Cluster 19 circRNAs were down-regulated at all stages and related to muscle structure and development. Lastly, the SDECs divided the period of skeletal muscle development into three transitional stages: stage 1 (F45 to F90), which related to muscle satellite cell proliferation and muscle fiber structure; stage 2 (F90 to B1), in which the attachment of the cytoplasmic surface to the actin cytoskeleton initiates; and stage 3, which involved the "cGMP-PKG signaling pathway". Moreover, the paraffin sections messages also validated that there are three transitional stages of skeletal muscle development. CONCLUSION: Our current study provides a catalog of goat muscle-related circRNAs that can stratify skeletal muscle development fetus 45 days to newborn 90 days into three developmental stages. These findings better our understanding of functional transitions during mammalian muscle development.

Genomic analysis of worldwide sheep breeds reveals PDGFD as a major target of fat-tail selection in sheep.

BACKGROUND: Fat tail is a unique trait in sheep acquired during domestication. Several genomic analyses have been conducted in sheep breeds from limited geographic origins to identify the genetic factors underlying this trait. Nevertheless, these studies obtained different candidates. The results of these regional studies were easily biased by the breed structures. RESULTS: To minimize the bias and distinguish the true candidates, we used an extended data set of 968 sheep representing 18 fat-tailed breeds and 14 thin-tailed breeds from around the world, and integrated two statistical tests to detect selection signatures, including Genetic Fixation Index (FST) and difference of derived allele frequency (ΔDAF). The results showed that platelet derived growth factor D (PDGFD) exhibited the highest genetic differentiation between fat- and thin-tailed sheep breeds. Analysis of sequence variation identified that a 6.8-kb region within the first intron of PDGFD is likely the target of positive selection and contains regulatory mutation(s) in fat-tailed sheep. Histological and gene expression analyses demonstrated that PDGFD expression is associated with maturation and hemostasis of adipocytes. Further retrospective analysis of public transcriptomic datasets revealed that PDGFD expression is down-regulated during adipogenesis in both human and mouse, and is higher in fat tissues of obese individuals than that in lean individuals. CONCLUSIONS: These results reveal that PDGFD is the predominant factor for the fat tail phenotype in sheep by contributing to adiopogenesis and maintaining the hemostasis of mature adipocytes. This study provides insights into the selection of fat-tailed sheep and has important application to animal breeding, as well as obesity-related human diseases.

EPAS1 Gain-of-Function Mutation Contributes to High-Altitude Adaptation in Tibetan Horses.

High altitude represents some of the most extreme environments worldwide. The genetic changes underlying adaptation to such environments have been recently identified in multiple animals but remain unknown in horses. Here, we sequence the complete genome of 138 domestic horses encompassing a whole altitudinal range across China to uncover the genetic basis for adaptation to high-altitude hypoxia. Our genome data set includes 65 lowland animals across ten Chinese native breeds, 61 horses living at least 3,300 m above sea level across seven locations along Qinghai-Tibetan Plateau, as well as 7 Thoroughbred and 5 Przewalski's horses added for comparison. We find that Tibetan horses do not descend from Przewalski's horses but were most likely introduced from a distinct horse lineage, following the emergence of pastoral nomadism in Northwestern China ∼3,700 years ago. We identify that the endothelial PAS domain protein 1 gene (EPAS1, also HIF2A) shows the strongest signature for positive selection in the Tibetan horse genome. Two missense mutations at this locus appear strongly associated with blood physiological parameters facilitating blood circulation as well as oxygen transportation and consumption in hypoxic conditions. Functional validation through protein mutagenesis shows that these mutations increase EPAS1 stability and its hetero dimerization affinity to ARNT (HIF1B). Our study demonstrates that missense mutations in the EPAS1 gene provided key evolutionary molecular adaptation to Tibetan horses living in high-altitude hypoxic environments. It reveals possible targets for genomic selection programs aimed at increasing hypoxia tolerance in livestock and provides a textbook example of evolutionary convergence across independent mammal lineages.