Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep.

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

The ribosomal chaperone NACA recruits PHD2 to cotranslationally modify HIF-α.

Prolyl hydroxylase domain protein 2 (PHD2)-catalyzed modification of hypoxia-inducible factor (HIF)-α is a key event in oxygen sensing. We previously showed that the zinc finger of PHD2 binds to a Pro-Xaa-Leu-Glu (PXLE) motif. Here, we show that the zinc finger binds to this motif in the ribosomal chaperone nascent polypeptide complex-α (NACA). This recruits PHD2 to the translation machinery to cotranslationally modify HIF-α. Importantly, this cotranslational modification is enhanced by a translational pause sequence in HIF-α. Mice with a knock-in Naca gene mutation that abolishes the PXLE motif display erythrocytosis, a reflection of HIF pathway dysregulation. In addition, human erythrocytosis-associated mutations in the zinc finger of PHD2 ablate interaction with NACA. Tibetans, who have adapted to the hypoxia of high altitude, harbor a PHD2 variant that we previously showed displays a defect in zinc finger binding to p23, a PXLE-containing HSP90 cochaperone. We show here that Tibetan PHD2 maintains interaction with NACA, thereby showing differential interactions with PXLE-containing proteins and providing an explanation for why Tibetans are not predisposed to erythrocytosis.

Strengthening the argument for a large Greater India.

The Sangdanlin section in southern Tibet represents a geologic Rosetta stone to constrain the initiation of the India-Asia collision from its sedimentary and paleomagnetic records. However, geoscientists have arrived at fundamentally divergent interpretations surrounding the age of the strata and its paleomagnetic record. Here, we report paleontologic, petrographic, and paleomagnetic data from the Sangdanlin section that recognize the sequence as a thrust complex containing interlaced Barremian-Albian (Early Cretaceous) and Paleocene strata, each separated by thrust faults. Recognizing two complexly interwoven formations of distinctly different ages contradicts a continuous stratigraphic superposition. Assigning an Early Cretaceous, instead of Paleocene, age to the units collected for paleomagnetic data revises paleogeographic models thereby supporting a large (2,000 to 3,000 km) extent of Greater India, with collision initiating at 55 ± 5 Ma in the western Himalayas. A contiguous plate in the Neotethys Ocean precludes that Asia's southern margin was built through a succession of accreted terrains.

Tectonic and orbital forcing of the South Asian monsoon in central Tibet during the late Oligocene.

The modern pattern of the Asian monsoon is thought to have formed around the Oligocene/Miocene transition and is generally attributed to Himalaya-Tibetan Plateau (H-TP) uplift. However, the timing of the ancient Asian monsoon over the TP and its response to astronomical forcing and TP uplift remains poorly known because of the paucity of well-dated high-resolution geological records from the TP interior. Here, we present a precession-scale cyclostratigraphic sedimentary section of 27.32 to 23.24 million years ago (Ma) during the late Oligocene epoch from the Nima Basin to show that the South Asian monsoon (SAM) had already advanced to the central TP (32°N) at least by 27.3 Ma, which is indicated by cyclic arid-humid fluctuations based on environmental magnetism proxies. A shift of lithology and astronomically orbital periods and amplified amplitude of proxy measurements as well as a hydroclimate transition around 25.8 Ma suggest that the SAM intensified at ~25.8 Ma and that the TP reached a paleoelevation threshold for enhancing the coupling between the uplifted plateau and the SAM. Orbital short eccentricity-paced precipitation variability is argued to be mainly driven by orbital eccentricity-modulated low-latitude summer insolation rather than glacial-interglacial Antarctic ice sheet fluctuations. The monsoon data from the TP interior provide key evidence to link the greatly enhanced tropical SAM at 25.8 Ma with TP uplift rather than global climate change and suggest that SAM's northward expansion to the boreal subtropics was dominated by a combination of tectonic and astronomical forcing at multiple timescales in the late Oligocene epoch.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep.

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

Tibetan mesenchymal stem cell-derived exosomes alleviate pulmonary vascular remodeling in hypoxic pulmonary hypertension rats.

Hypoxic pulmonary hypertension (HPH) is characterized by progressive pulmonary vasoconstriction, vascular remodeling, and right ventricular hypertrophy, causing right heart failure. This study aimed to investigate the therapeutic effects of exosomes from Tibetan umbilical cord mesenchymal stem cells on HPH via the TGF-ß1/Smad2/3 pathway, comparing them with exosomes from Han Chinese individuals. An HPH rat model was established in vivo, and a hypoxia-induced injury in the rat pulmonary artery smooth muscle cells (rPASMCs) was simulated in vitro. Exosomes from human umbilical cord mesenchymal stem cells were administered to HPH model rats or added to cultured rPASMCs. The therapeutic effects of Tibetan-mesenchymal stem cell-derived exosomes (Tibetan-MSC-exo) and Han-mesenchymal stem cell-derived exosomes (Han-MSC-exo) on HPH were investigated through immunohistochemistry, western blotting, EdU, and Transwell assays. The results showed that Tibetan-MSC-exo significantly attenuated pulmonary vascular remodeling and right ventricular hypertrophy in HPH rats compared with Han-MSC-exo. Tibetan-MSC-exo demonstrated better inhibition of hypoxia-induced rPASMCs proliferation and migration. Transcriptome sequencing revealed upregulated genes (Nbl1, Id2, Smad6, and Ltbp1) related to the TGFß pathway. Nbl1 knockdown enhanced hypoxia-induced rPASMCs proliferation and migration, reversing Tibetan-MSC-exo-induced downregulation of TGFß1 and p-Smad2/3. Furthermore, TGFß1 overexpression hindered the therapeutic effects of Tibetan-MSC-exo and Han-MSC-exo on hypoxic injury. These findings suggest that Tibetan-MSC-exo favors HPH treatment better than Han-MSC-exo, possibly through the modulation of the TGFß1/Smad2/3 pathway via Nbl1.

Pulsed rise and growth of the Tibetan Plateau to its northern margin since ca. 30 Ma.

The onset of mountain building along margins of the Tibetan Plateau provides a key constraint on the processes by which the high topography in Eurasia formed. Although progressive expansion of thickened crust underpins most models, several studies suggest that the northern extent of the plateau was established early, soon after the collision between India and Eurasia at ca. 50 Ma. This inference relies heavily on the age and provenance of Cenozoic sediments preserved in the Qaidam basin. Here, we present evidence in the northern plateau for a considerably younger inception and evolution of the Qaidam basin, based on magnetostratigraphies combined with detrital apatite fission-track ages that date the basin fills to be from ca. 30 to 4.8 Ma. Detrital zircon-provenance analyses coupled with paleocurrents reveal that two-stage growth of the Qilian Shan in the northeastern margin of the Tibetan Plateau began at ca. 30 and at 10 Ma, respectively. Evidence for ca. 30 and 10 to 15 Ma widespread synchronous deformation throughout the Tibetan Plateau and its margins suggests that these two stages of outward growth may have resulted from the removal of mantle lithosphere beneath different portions of the Tibetan Plateau.

Genetic adaptation of skin pigmentation in highland Tibetans.

Strong ultraviolet (UV) radiation at high altitude imposes a serious selective pressure, which may induce skin pigmentation adaptation of indigenous populations. We conducted skin pigmentation phenotyping and genome-wide analysis of Tibetans in order to understand the underlying mechanism of adaptation to UV radiation. We observe that Tibetans have darker baseline skin color compared with lowland Han Chinese, as well as an improved tanning ability, suggesting a two-level adaptation to boost their melanin production. A genome-wide search for the responsible genes identifies GNPAT showing strong signals of positive selection in Tibetans. An enhancer mutation (rs75356281) located in GNPAT intron 2 is enriched in Tibetans (58%) but rare in other world populations (0 to 18%). The adaptive allele of rs75356281 is associated with darker skin in Tibetans and, under UVB treatment, it displays higher enhancer activities compared with the wild-type allele in in vitro luciferase assays. Transcriptome analyses of gene-edited cells clearly show that with UVB treatment, the adaptive variant of GNPAT promotes melanin synthesis, likely through the interactions of CAT and ACAA1 in peroxisomes with other pigmentation genes, and they act synergistically, leading to an improved tanning ability in Tibetans for UV protection.

Seismic constraints and geodynamic implications of differential lithosphere-asthenosphere flow revealed in East Asia.

During the last 50 Ma, the East Asian continent has been a zone of massive continental collision and lithospheric deformation. While the consequences of this for Asian surface and lithospheric deformation have been intensively studied over the past 4 decades, the relationships between lithospheric deformation and underlying asthenospheric flow have been more difficult to constrain. Here we present a high resolution 3-D azimuthal anisotropy model for the northeastern Tibetan Plateau and its eastward continuation based on surface-wave tomography and shear-wave splitting measurements. This model shows that eastward lateral flow of asthenosphere beneath the northeastern Tibetan Plateau is being blocked by thick Ordos and Sichuan cratonic keels. The damming effect of these keels induces flow to first rotate around the Ordos keel and then transition into strong east-west flow beneath the thinner lithosphere that forms the lithospheric suture between the two cratonic keels. We further find that asthenosphere flow directions can differ from those of overlying lithosphere, with the asthenosphere neither being passively dragged by overlying lithosphere, nor being able to drag the overlying plate to mimic its subsurface flow. Finally, the region of eastward-channeled asthenospheric flow from Tibet underlies a belt of stronger intracontinental deformation in eastern China.

Limited underthrusting of India below Tibet: 3He/4He analysis of thermal springs locates the mantle suture in continental collision.

During continent­continent collision, does the downgoing continental plate underplate far inboard of the collisional boundary or does it subduct steeply into the mantle, and how is this geometry manifested in the mantle flow field? We test conflicting models for these questions for Earth's archetypal continental collision forming the Himalaya and Tibetan Plateau. Air-corrected helium isotope data (3He/4He) from 225 geothermal springs (196 from our group, 29 from the literature) delineate a boundary separating a Himalayan domain of only crustal helium from a Tibetan domain with significant mantle helium. This 1,000-km-long boundary is located close to the Yarlung-Zangbo Suture (YZS) in southern Tibet from 80 to 92°E and is interpreted to overlie the "mantle suture" where cold underplated Indian lithosphere is juxtaposed at >80 km depth against a sub-Tibetan incipiently molten asthenospheric mantle wedge. In southeastern Tibet, the mantle suture lies 100 km south of the YZS, implying delamination of the mantle lithosphere from the Indian crust. This helium-isotopic boundary helps resolve multiple, mutually conflicting seismological interpretations. Our synthesis of the combined data locates the northern limit of Indian underplating beneath Tibet, where the Indian plate bends to steeper dips or breaks off beneath a (likely thin) asthenospheric wedge below Tibetan crust, thereby defining limited underthrusting for the Tibetan continental collision.

Rupture process of the 2021 M7.4 Maduo earthquake and implication for deformation mode of the Songpan-Ganzi terrane in Tibetan Plateau.

The deformation mode of the Tibetan Plateau is of crucial importance for understanding its construction and extrusion processes, as well as for the assessment of regional earthquake potential. Block motion and viscous flow models have been proposed to describe the deformation field but are not fully supported by modern geophysical observations. The 2021 Mw 7.4 Maduo earthquake, which occurred inside the Songpan-Ganzi terrane (SGT) in central-east Tibet, provides a chance to evaluate the associated deformation mode of the region. We conduct a joint inversion for this earthquake and resolve a bilateral rupture process, which is characterized by super- and subshear rupture velocities, respectively. We interpret this distinct rupture behavior to be the result of the respective slip concentration depths of the two ruptured segments. We analyze geological, seismic, and geodetic evidence and find that the SGT upper crust shows distributed shear deformation and distinct transverse anisotropy, which are associated with folded structures originating from compression of the paleo-Tethys ocean accretional prism realigned by following shear deformation. The SGT receives lateral shear loading from its NS boundary and accommodates a right-step sinistral motion across the terrane boundary faults. The unique tectonic setting of the SGT defines locations and behaviors of internal faulting and strong earthquakes such as the 2021 Maduo earthquake, with the latter occurring on slow-moving faults at intervals of several thousands of years.

Use of δ18Oatm in dating a Tibetan ice core record of Holocene/Late Glacial climate.

Ice cores from the northwestern Tibetan Plateau (NWTP) contain long records of regional climate variability, but refrozen meltwater and dust in these cores has hampered development of robust timescales. Here, we introduce an approach to dating the ice via the isotopic composition of atmospheric O2 in air bubbles (δ18Oatm), along with annual layer counting and radiocarbon dating. We provide a robust chronology for water isotope records (δ18Oice and d-excess) from three ice cores from the Guliya ice cap in the NWTP. The measurement of δ18Oatm, although common in polar ice core timescales, has rarely been used on ice cores from low-latitude, high-altitude glaciers due to (1) low air pressure, (2) the common presence of refrozen melt that adds dissolved gases and reduces the amount of air available for analysis, and (3) the respiratory consumption of molecular oxygen (O2) by micro-organisms in the ice, which fractionates the δ18O of O2 from the atmospheric value. Here, we make corrections for melt and respiration to address these complications. The resulting records of water isotopes from the Guliya ice cores reveal climatic variations over the last 15,000 y, the timings of which correspond to those observed in independently dated lake and speleothem records and confirm that the Guliya ice cap existed before the Holocene. The millennial-scale drivers of δ18Oice are complex and temporally variable; however, Guliya δ18Oice values since the mid-20th century are the highest since the beginning of the Holocene and have increased with regional air temperature.

Stress adaptation in Tibetan cashmere goats is governed by inherent metabolic differences and manifested through variable cashmere phenotypes.

Tibetan cashmere goats are not only served as a valuable model for studying adaptation to hypoxia and high-altitude conditions but also playing a pivotal role in bolstering local economies through the provision of premium quality cashmere yarn. In this study, we performed an integration and network analysis of metabolomic, transcriptomic and proteomic to elucidate the role of differentially expressed genes, important metabolites, and relevant cellular and metabolic pathways between the fine (average 12.04 ± 0.03 µm of mean fiber diameter) and coarse cashmere (average 14.88 ± 0.05 µm of mean fber diameter) producing by Tibetan cashmere goats. We identified a distinction of 56 and 71 differential metabolites (DMs) between the F and C cashmere groups under positive and negative ion modes, respectively. The KEGG pathway enrichment analysis of these DMs highlighted numerous pathways predominantly involved in amino acid and protein metabolism, as indicated by the finding that the most impactful pathway was the mammalian target of rapamycin (mTOR) signalling pathway. In the F group, we identified a distinctive metabolic profile where amino acid metabolites including serine, histidine, asparagine, glutamic acid, arginine, valine, aspartic acid, tyrosine, and methionine were upregulated, while lysine, isoleucine, glutamine, tryptophan, and threonine were downregulated. The regulatory network and gene co-expression network revealed crucial genes, metabolites, and metabolic pathways. The integrative omics analysis revealed a high enrichment of several pathways, notably encompassing protein digestion and absorption, sphingolipid signalling, and the synaptic vesicle cycle. Within the sphere of our integrative analysis, DNMT3B was identified as a paramount gene, intricately associated with significant proteins such as HMCN1, CPB2, GNG12, and LRP1. Our present study delineated the molecular underpinnings governing the variations in cashmere characteristics by conducting comprehensive analyses across metabolomic, transcriptomic, and proteomic dimensions. This research provided newly insights into the mechanisms regulating cashmere traits and facilitated the advancement of selective breeding programs aimed at cultivating high-quality superfine Tibetan cashmere goats.

Metabolomics and transcriptomics of embryonic livers reveal hypoxia adaptation of Tibetan chickens.

BACKGROUND: Exploring the hypoxia adaptation mechanism of Tibetan chicken is of great significance for revealing the survival law of Tibetan chicken and plateau animal husbandry production. To investigate the hypoxia adaptation of Tibetan chickens (TBCs), an integrative metabolomic-transcriptomic analysis of the liver on day 18 of embryonic development was performed. Dwarf laying chickens (DLCs), a lowland breed, were used as a control. RESULTS: A total of 1,908 metabolites were identified in both TBCs and DLCs. Energy metabolism and amino acid metabolism related differentially regulated metabolites (DRMs) were significantly enriched under hypoxia. Important metabolic pathways including the TCA cycle and arginine and proline metabolism were screened; PCK1, SUCLA2, and CPS1 were found to be altered under hypoxic conditions. In addition, integrated analysis suggested potential differences in mitochondrial function, which may play a crucial role in the study of chicken oxygen adaptation. CONCLUSIONS: These results suggest that hypoxia changed the gene expression and metabolic patterns of embryonic liver of TBCs compared to DLCs. Our study provides a basis for uncovering the molecular regulation mechanisms of hypoxia adaptation in TBCs with the potential application of hypoxia adaptation research for other animals living on the Qinghai-Tibet plateau, and may even contribute to the study of diseases caused by hypoxia.

Genome-wide association studies of body size traits in Tibetan sheep.

BACKGROUND: Elucidating the genetic variation underlying phenotypic diversity will facilitate improving production performance in livestock species. The Tibetan sheep breed in China holds significant historical importance, serving as a fundamental pillar of Qinghai's animal husbandry sector. The Plateau-type Tibetan sheep, comprising 90% of the province's population, are characterized by their tall stature and serve as the primary breed among Tibetan sheep. In contrast, Zhashijia sheep exhibit larger size and superior meat quality. These two species provide an excellent model for elucidating the genetic basis of body size variation. Therefore, this study aims to conduct a comprehensive genome-wide association study on these two Tibetan sheep breeds to identify single nucleotide polymorphism loci and regulatory genes that influence body size traits in Tibetan sheep. RESULT: In this study, the phenotypic traits of body weight, body length, body height, chest circumference, chest depth, chest width, waist angle width, and pipe circumference were evaluated in two Tibetan sheep breeds: Plateau-type sheep and Zhashijia Tibetan sheep. Whole genome sequencing generated 48,215,130 high-quality SNPs for genome-wide association study. Four methods were applied and identified 623 SNPs significantly associated with body size traits. The significantly associated single nucleotide polymorphisms identified in this study are located near or within 111 candidate genes. These genes exhibit enrichment in the cAMP and Rap1 signaling pathways, significantly affecting animal growth, and body size. Specifically, the following genes were associated: ASAP1, CDK6, FRYL, NAV2, PTPRM, GPC6, PTPRG, KANK1, NTRK2 and ADCY8. CONCLUSION: By genome-wide association study, we identified 16 SNPs and 10 candidate genes associated with body size traits in Tibetan sheep, which hold potential for application in genomic selection breeding programs in sheep. Identifying these candidate genes will establish a solid foundation for applying molecular marker-assisted selection in sheep breeding and improve our understanding of body size control in farmed animals.

Effect of supplementing lysins and methionine to low-protein diets on growth performance, hepatic antioxidant capacity, immune status, and glycolytic activity of tibetan sheep.

Reducing the levels of dietary protein is an effective nutritional approach in lowering feed cost and nitrogen emissions in ruminants. The purpose of this study was to evaluate the effects of dietary Lys/Met ratio in a low protein diet (10%, dry matter basis) on the growth performance and hepatic function (antioxidant capacity, immune status, and glycolytic activity) in Tibetan lambs. Ninety two-month-old rams with an average weight of 15.37 ± 0.92 kg were randomly assigned to LP-L (dietary Lys/Met = 1:1), LP-M (dietary Lys/Met = 2:1) and LP-H (dietary Lys/Met = 3:1) treatments. The trial was conducted over 100 d, including 10 d of adaption to the diets. Hepatic phenotypes, antioxidant capacity, immune status, glycolytic activity and gene expression profiling was detected after the conclusion of the feeding trials. The results showed that the body weight was higher in the LP-L group when compared to those on the LP-M group (P < 0.05). In addition, the activities of the catalase (CAT) and glutathione peroxidase (GSH-Px) in the LP-L group were significantly increased compared with the LP-M group (P < 0.05), while the malondialdehyde (MDA) levels in LP-H group were significantly decreased (P < 0.05). Compared with LP-H group, both hepatic glycogen (P < 0.01) and lactate dehydrogenase (LDH) (P < 0.05) were significantly elevated in LP-L group. For the LP-L group, the hepatocytes were arranged radially with the central vein in the center, and hepatic plates exhibited tight arrangement. Transcriptome analysis identified 29, 179, and 129 differentially expressed genes (DEGs) between the LP-M vs. LP-L, LP-H vs. LP-M, and LP-H vs. LP-L groups, respectively (Q-values < 0.05 and |log2Fold Change| > 1). Gene Ontology (GO) and correlation analyses showed that in the LP-L group, core genes (C1QA and JUNB) enriched in oxidoreductase activity were positively correlated with antioxidant indicators, while the MYO9A core gene enriched in the immune response was positively associated with immune indicators, and core genes enriched in molecular function (PDK3 and PDP2) were positively correlated with glycolysis indicators. In summary, low-protein diet with a low Lys/Met ratio (1:1) could reduce the hepatic oxidative stress and improve the glycolytic activity by regulating the expression of related genes of Tibetan sheep.

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.

Selection signatures of wool color in Gangba sheep revealed by genome-wide SNP discovery.

BACKGROUND: Gangba sheep as a famous breed of Tibetan sheep, its wool color is mainly white and black. Gangba wool is economically important as a high-quality raw material for Tibetan blankets and Tibetan serge. However, relatively few studies have been conducted on the wool color of Tibetan sheep. RESULTS: To fill this research gap, this study conducted an in-depth analysis of two populations of Gangba sheep (black and white wool color) using whole genome resequencing to identify genetic variation associated with wool color. Utilizing PCA, Genetic Admixture, and N-J Tree analyses, the present study revealed a consistent genetic relationship and structure between black and white wool colored Gangba sheep populations, which is consistent with their breed history. Analysis of selection signatures using multiple methods (FST, π ratio, Tajima's D), 370 candidate genes were screened in the black wool group (GBB vs GBW); among them, MC1R, MLPH, SPIRE2, RAB17, SMARCA4, IRF4, CAV1, USP7, TP53, MYO6, MITF, MC2R, TET2, NF1, JAK1, GABRR1 genes are mainly associated with melanin synthesis, melanin delivery, and distribution. The enrichment results of the candidate genes identified 35 GO entries and 19 KEGG pathways associated with the formation of the black phenotype. 311 candidate genes were screened in the white wool group (GBW vs GBB); among them, REST, POU2F1, ADCY10, CCNB1, EP300, BRD4, GLI3, and SDHA genes were mainly associated with interfering with the differentiation of neural crest cells into melanocytes, affecting the proliferation of melanocytes, and inhibiting melanin synthesis. 31 GO entries and 22 KEGG pathways were associated with the formation of the white phenotype. CONCLUSIONS: This study provides important information for understanding the genetic mechanism of wool color in Gangba, and provides genetic knowledge for improving and optimizing the wool color of Tibetan sheep. Genetic improvement and selective breeding to produce wool of specific colors can meet the demand for a diversity of wool products in the Tibetan wool textile market.

Whole-genome SNP allele frequency differences between Tibetan and Large white pigs reveal genes associated with skeletal muscle growth.

BACKGROUND: The skeletal muscle growth rate and body size of Tibetan pigs (TIB) are lower than Large white pigs (LW). However, the underlying genetic basis attributing to these differences remains uncertain. To address this knowledge gap, the present study employed whole-genome sequencing of TIB (slow growth) and LW (fast growth) individuals, and integrated with existing NCBI sequencing datasets of TIB and LW individuals, enabling the identification of a comprehensive set of genetic variations for each breed. The specific and predominant SNPs in the TIB and LW populations were detected by using a cutoff value of 0.50 for SNP allele frequency and absolute allele frequency differences (△AF) between the TIB and LW populations. RESULTS: A total of 21,767,938 SNPs were retrieved from 44 TIB and 29 LW genomes. The analysis detected 2,893,106 (13.29%) and 813,310 (3.74%) specific and predominant SNPs in the TIB and LW populations, and annotated to 24,560 genes. Further GO analysis revealed 291 genes involved in biological processes related to striated and/or skeletal muscle differentiation, proliferation, hypertrophy, regulation of striated muscle cell differentiation and proliferation, and myoblast differentiation and fusion. These 291 genes included crucial regulators of muscle cell determination, proliferation, differentiation, and hypertrophy, such as members of the Myogenic regulatory factors (MRF) (MYOD, MYF5, MYOG, MYF6) and Myocyte enhancer factor 2 (MEF2) (MEF2A, MEF2C, MEF2D) families, as well as muscle growth inhibitors (MSTN, ACVR1, and SMAD1); KEGG pathway analysis revealed 106 and 20 genes were found in muscle growth related positive and negative regulatory signaling pathways. Notably, genes critical for protein synthesis, such as MTOR, IGF1, IGF1R, IRS1, INSR, and RPS6KA6, were implicated in these pathways. CONCLUSION: This study employed an effective methodology to rigorously identify the potential genes associated with skeletal muscle development. A substantial number of SNPs and genes that potentially play roles in the divergence observed in skeletal muscle growth between the TIB and LW breeds were identified. These findings offer valuable insights into the genetic underpinnings of skeletal muscle development and present opportunities for enhancing meat production through pig breeding.

Genomic Variation, Population History, and Long-Term Genetic Adaptation to High Altitudes in Tibetan Partridge (Perdix hodgsoniae).

Species residing across elevational gradients display adaptations in response to environmental changes such as oxygen availability, ultraviolet radiation, and temperature. Here, we study genomic variation, gene expression, and long-term adaptation in Tibetan Partridge (Perdix hodgsoniae) populations residing across the elevational gradient of the Tibetan Plateau. We generated a high-quality draft genome and used it to carry out downstream population genomic and transcriptomic analysis. The P. hodgsoniae populations residing across various elevations were genetically distinct, and their phylogenetic clustering was consistent with their geographic distribution. We identified possible evidence of gene flow between populations residing in <3,000 and >4,200 m elevation that is consistent with known habitat expansion of high-altitude populations of P. hodgsoniae to a lower elevation. We identified a 60 kb haplotype encompassing the Estrogen Receptor 1 (ESR1) gene, showing strong genetic divergence between populations of P. hodgsoniae. We identified six single nucleotide polymorphisms within the ESR1 gene fixed for derived alleles in high-altitude populations that are strongly conserved across vertebrates. We also compared blood transcriptome profiles and identified differentially expressed genes (such as GAPDH, LDHA, and ALDOC) that correlated with differences in altitude among populations of P. hodgsoniae. These candidate genes from population genomics and transcriptomics analysis were enriched for neutrophil degranulation and glycolysis pathways, which are known to respond to hypoxia and hence may contribute to long-term adaptation to high altitudes in P. hodgsoniae. Our results highlight Tibetan Partridges as a useful model to study molecular mechanisms underlying long-term adaptation to high altitudes.