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.

A comparative analysis of differentially expressed mRNAs, miRNAs and circRNAs provides insights into the key genes involved in the high-altitude adaptation of yaks.

BACKGROUND: Yaks that inhabit the Tibetan Plateau exhibit striking phenotypic and physiological differences from cattle and have adapted well to the extreme conditions on the plateau. However, the mechanisms used by these animals for the regulation of gene expression at high altitude are not fully understood. RESULTS: Here, we sequenced nine lung transcriptomes of yaks at altitudes of 3400, 4200 and 5000 m, and low-altitude Zaosheng cattle, which is a closely related species, served as controls. The analysis identified 21,764 mRNAs, 1377 circRNAs and 1209 miRNAs. By comparing yaks and cattle, 4975 mRNAs, 252 circRNAs and 75 miRNAs were identified differentially expressed. By comparing yaks at different altitudes, we identified 756 mRNAs, 64 circRNAs and 83 miRNAs that were differentially expressed (fold change ≥2 and P-value < 0.05). The pathways enriched in the mRNAs, circRNAs and miRNAs identified from the comparison of yaks and cattle were mainly associated with metabolism, including 'glycosaminoglycan degradation', 'pentose and glucuronate interconversions' and 'flavone and flavonol biosynthesis', and the mRNAs, circRNAs and miRNAs identified from the comparison of yaks at different altitude gradients were significantly enriched in metabolic pathways and immune and genetic information processing pathways. The core RNAs were identified from the mRNA-miRNA-circRNA networks constructed using the predominant differentially expressed RNAs. The core genes specific to the difference between yaks and cattle were associated with the endoplasmic reticulum and fat deposition, but those identified from the comparison among yaks at different altitude gradients were associated with maintenance of the normal biological functions of cells. CONCLUSIONS: This study enhances our understanding of the molecular mechanisms involved in hypoxic adaptation in yaks and might contribute to improvements in the understanding and prevention of hypoxia-related diseases.

Evaluating genetic diversity and identifying priority conservation for seven Tibetan pig populations in China based on the mtDNA D-loop.

OBJECTIVE: Tibetan pigs, an excellent species unique to China, face serious threats, which in turn affects the development and utilization of the outstanding advantages of plateau hypoxia adaptability and reduces their genetic diversity. Therefore, a discussion of measures to conserve this genetic resource is necessary. The method, based on genetic diversity, genetic divergence and total genetic contribution rate of population, reflects the priority conservation order and varies depending on the three different purposes of conservation. METHODS: We analyzed mitochondrial DNA control region (D-loop) variation in 1,201 individuals from nine Tibetan pig populations across five provinces and downloaded 564 mtDNA D-loop sequences from three indigenous pig breeds in Qinghai, Sichuan, and Yunnan Provinces distributed near the Tibetan pigs. RESULTS: We analyzed three different aspects: Changdu Tibetan pigs have the highest genetic diversity, and from the perspective of genetic diversity, the priority conservation is Changdu Tibetan pigs. Hezuo Tibetan pigs have the highest genetic contribution, so the priority conservation is Hezuo Tibetan pigs in the genetic contribution aspect. Rkaze Tibetan pigs were severely affected by indigenous pig breeds, so if considering from the perspective of introgression, the priority conservation is Rkaze Tibetan pigs. CONCLUSION: This study evaluated genetic diversity and comprehensively assessed conservation priority from three different aspects in nine Tibetan pig populations.

The domestication event of the Tibetan pig revealed to be in the upstream region of the Yellow River based on the mtDNA D-loop.

OBJECTIVE: Evidence from previous reports indicates that pig domestication in East Asia mainly occurred in the Mekong region and the middle and downstream regions of the Yangtze River. Further research identified two new origin centers for domestic pigs in the Tibetan Plateau and the islands of Southeast Asia. However, due to the small sample size of Tibetan pigs, details of the origin and spread of Tibetan pigs has not yet been established. METHODS: We analyzed mitochondrial DNA control region (D-loop) variation in 1,201 individuals from nine Tibetan pig populations across five provinces. Comprehensive Tibetan pig samples were taken to perform the most detailed analysis of Tibetan pigs to date. RESULTS: The result indicate that Rkaze pigs had the lowest level of diversity, while Changdu pigs had the highest diversity. Interestingly, these two populations were both in the Tibetan Plateau area. If we calculate diversity in terms of each province, the Tibetan Plateau area had the lowest diversity, while the Chinese province of Gansu had the highest diversity. Diversity gradient analysis of major haplotypes suggested three domestication centers of Tibetan pigs in the Tibetan Plateau and the Chinese provinces of Gansu and Yunnan. CONCLUSION: We found two new domestication centers for Tibetan pigs. One is in the Chinese province of Gansu, which lies in the upstream region of the Yellow River, and the other is in the Chinese province of Yunnan.

Sex-biased regulatory changes in the placenta of native highlanders contribute to adaptive fetal development.

Compared with lowlander migrants, native Tibetans have a higher reproductive success at high altitude though the underlying mechanism remains unclear. Here, we compared the transcriptome and histology of full-term placentas between native Tibetans and Han migrants. We found that the placental trophoblast shows the largest expression divergence between Tibetans and Han, and Tibetans show decreased immune response and endoplasmic reticulum stress. Remarkably, we detected a sex-biased expression divergence, where the male-infant placentas show a greater between-population difference than the female-infant placentas. The umbilical cord plays a key role in the sex-biased expression divergence, which is associated with the higher birth weight of the male newborns of Tibetans. We also identified adaptive histological changes in the male-infant placentas of Tibetans, including larger umbilical artery wall and umbilical artery intima and media, and fewer syncytial knots. These findings provide valuable insights into the sex-biased adaptation of human populations, with significant implications for medical and genetic studies of human reproduction.

High Arterial Oxygen Saturation in the Acclimatized Lowlanders Living at High Altitude.

Blood oxygen saturation (SpO2) is a key indicator of oxygen availability in the body. It is known that a low SpO2 at high altitude is associated with morbidity and mortality risks due to physiological hypoxemia. Previously, it was proposed that the lowlander immigrants living at high altitude should have a lower SpO2 level compared to the highlander natives, but this proposal has not been rigorously tested due to the lack of data from the lowlander immigrants living at high altitude. In this study, we compared arterial oxygen saturation of 5929 Tibetan natives and 1034 Han Chinese immigrants living at altitudes ranging from 1120 m to 5020 m. Unexpectedly, the Han immigrants had a higher SpO2 than the Tibetan natives at the same high altitudes. At the same time, there is a higher prevalence of chronic mountain sickness in Han than in Tibetans at the same altitude. This result suggests that the relatively higher SpO2 level of the acclimatized Han is associated with a physiological cost, and the SpO2 level of Tibetans tends to be sub-optimal. Consequently, SpO2 alone is not a robust indicator of physiological performance at high altitude. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-023-00117-x.

Polygenic adaptation leads to a higher reproductive fitness of native Tibetans at high altitude.

The adaptation of Tibetans to high-altitude environments has been studied extensively. However, the direct assessment of evolutionary adaptation, i.e., the reproductive fitness of Tibetans and its genetic basis, remains elusive. Here, we conduct systematic phenotyping and genome-wide association analysis of 2,252 mother-newborn pairs of indigenous Tibetans, covering 12 reproductive traits and 76 maternal physiological traits. Compared with the lowland immigrants living at high altitudes, indigenous Tibetans show better reproductive outcomes, reflected by their lower abortion rate, higher birth weight, and better fetal development. The results of genome-wide association analyses indicate a polygenic adaptation of reproduction in Tibetans, attributed to the genomic backgrounds of both the mothers and the newborns. Furthermore, the EPAS1-edited mice display higher reproductive fitness under chronic hypoxia, mirroring the situation in Tibetans. Collectively, these results shed new light on the phenotypic pattern and the genetic mechanism of human reproductive fitness in extreme environments.

Deep phenotyping of 11,880 highlanders reveals novel adaptive traits in native Tibetans.

Tibetans are the ideal population to study genetic adaptation in extreme environments. Here, we performed systematic phenotyping of 11,880 highlanders, covering 133 quantitative traits of 13 organ systems. We provided a comprehensive phenotypic atlas by comparing altitude adaptation and altitude acclimatization. We found the differences between adaptation and acclimatization are quantitative rather than qualitative, with a whole-system "blunted effect" seen in the adapted Tibetans. We characterized twelve different functional changes between adaptation and acclimatization. More importantly, we established a landscape of adaptive phenotypes of indigenous Tibetans, including 45 newly identified Tibetan adaptation-nominated traits, involving specific changes of Tibetans in internal organ state, metabolism, eye morphology, and skin pigmentation. In addition, we observed a sex-biased pattern between altitude acclimatization and adaptation. The generated atlas of phenotypic landscape provides new insights into understanding of human adaptation to high-altitude environments, and it serves as a valuable blueprint for future medical and physiological studies.

Large-scale genome sequencing redefines the genetic footprints of high-altitude adaptation in Tibetans.

BACKGROUND: Tibetans are genetically adapted to high-altitude environments. Though many studies have been conducted, the genetic basis of the adaptation remains elusive due to the poor reproducibility for detecting selective signatures in the Tibetan genomes. RESULTS: Here, we present whole-genome sequencing (WGS) data of 1001 indigenous Tibetans, covering the major populated areas of the Qinghai-Tibetan Plateau in China. We identify 35 million variants, and more than one-third of them are novel variants. Utilizing the large-scale WGS data, we construct a comprehensive map of allele frequency and linkage disequilibrium and provide a population-specific genome reference panel, referred to as 1KTGP. Moreover, with the use of a combined approach, we redefine the signatures of Darwinian-positive selection in the Tibetan genomes, and we characterize a high-confidence list of 4320 variants and 192 genes that have undergone selection in Tibetans. In particular, we discover four new genes, TMEM132C, ATP13A3, SANBR, and KHDRBS2, with strong signals of selection, and they may account for the adaptation of cardio-pulmonary functions in Tibetans. Functional annotation and enrichment analysis indicate that the 192 genes with selective signatures are likely involved in multiple organs and physiological systems, suggesting polygenic and pleiotropic effects. CONCLUSIONS: Overall, the large-scale Tibetan WGS data and the identified adaptive variants/genes can serve as a valuable resource for future genetic and medical studies of high-altitude populations.

Seasonality and Sex-Biased Fluctuation of Birth Weight in Tibetan Populations.

Birth weight (BW) is a key determinant of infant mortality. Previous studies have reported seasonal fluctuation of BW. However, the responsible environmental factors remain disputable. High-altitude environment provides a great opportunity to test the current hypotheses due to its distinctive climate conditions. We collected BW data of ~ 9000 Tibetan singletons born at Lhasa (elevation: 3660 m) from 2014 to 2018. Using regression models, we analyzed BW seasonality of highland Tibetans. Multivariate models with meteorological factors as independent variables were employed to examine responsible environmental factors accounting for seasonal variation. We compared BW, low-BW prevalence and sex ratio between highland and lowland populations, and we observed a significant seasonal pattern of BW in Tibetans, with a peak in winter and a trough in summer. Notably, there is a marked sex-biased pattern of BW seasonality (more striking in males than in females). Sunlight exposure in the 3rd trimester and barometric pressure exposure in the 2nd trimester are significantly correlated with BW, and the latter can be explained by seasonal change of oxygen partial pressure. In particular, due to the male-biased BW seasonality, we found a more serious BW reduction and higher prevalence of low-BW in males, and a skewed sex ratio in highlanders. The infant BW of highland Tibetans has a clear pattern of seasonality. The winter BW is larger than the summer BW, due to the longer sunlight exposure during the late-trimester. Male infants are more sensitive to hypoxia than female infants during the 2nd trimester, leading to more BW reduction and higher mortality. Supplementary Information: The online version contains supplementary material available at 10.1007/s43657-021-00038-7.

Molecular mechanisms detected in yak lung tissue via transcriptome-wide analysis provide insights into adaptation to high altitudes.

Due to their long-term colonization of and widespread distribution in plateau environments, yaks can serve as an ideal natural animal model for the adaptive evolution of other plateau species, including humans. Some studies reported that the lung and heart are two key organs that show adaptive transcriptional changes in response to high altitudes, and most of the genes that show differential expression in lung tissue across different altitudes display nonlinear regulation. To explore the molecular mechanisms that are activated in yak lung tissue in response to hypoxia, the mRNAs, lncRNAs and miRNAs of lung tissue from 9 yaks living at three different altitudes (3400 m, 4200 m and 5000 m), with three repetitions per altitude, were sequenced. Two Zaosheng cattle from 1500 m were selected as low-altitude control. A total of 21,764 mRNAs, 14,168 lncRNAs and 1209 miRNAs (305 known and 904 novel miRNAs) were identified. In a comparison of yaks and cattle, 4975 mRNAs, 3326 lncRNAs and 75 miRNAs were differentially expressed. A total of 756 mRNAs, 346 lncRNAs and 83 miRNAs were found to be differentially expressed among yaks living at three different altitudes (fold change ≥ 2 and P-value < 0.05). The differentially expressed genes between yaks and cattle were functionally enriched in long-chain fatty acid metabolic process and protein processing, while the differentially expressed genes among yaks living at three different altitudes were enriched in immune response and the cell cycle. Furthermore, competing endogenous RNA (ceRNA) networks were investigated to illustrate the roles of ceRNAs in this process, the result was also support the GO and KEGG analysis. The present research provides important genomic insights for discovering the mechanisms that are activated in response to hypoxia in yak lung tissue.

De novo assembly of a Tibetan genome and identification of novel structural variants associated with high-altitude adaptation.

Structural variants (SVs) may play important roles in human adaptation to extreme environments such as high altitude but have been under-investigated. Here, combining long-read sequencing with multiple scaffolding techniques, we assembled a high-quality Tibetan genome (ZF1), with a contig N50 length of 24.57 mega-base pairs (Mb) and a scaffold N50 length of 58.80 Mb. The ZF1 assembly filled 80 remaining N-gaps (0.25 Mb in total length) in the reference human genome (GRCh38). Markedly, we detected 17 900 SVs, among which the ZF1-specific SVs are enriched in GTPase activity that is required for activation of the hypoxic pathway. Further population analysis uncovered a 163-bp intronic deletion in the MKL1 gene showing large divergence between highland Tibetans and lowland Han Chinese. This deletion is significantly associated with lower systolic pulmonary arterial pressure, one of the key adaptive physiological traits in Tibetans. Moreover, with the use of the high-quality de novo assembly, we observed a much higher rate of genome-wide archaic hominid (Altai Neanderthal and Denisovan) shared non-reference sequences in ZF1 (1.32%-1.53%) compared to other East Asian genomes (0.70%-0.98%), reflecting a unique genomic composition of Tibetans. One such archaic hominid shared sequence-a 662-bp intronic insertion in the SCUBE2 gene-is enriched and associated with better lung function (the FEV1/FVC ratio) in Tibetans. Collectively, we generated the first high-resolution Tibetan reference genome, and the identified SVs may serve as valuable resources for future evolutionary and medical studies.

Multiple Domestication Centers Revealed by the Geographical Distribution of Chinese Native Pigs.

Previous studies have shown that Southeast Asian pigs were independently domesticated from local wild boars. However, the domestication of Chinese native pigs remains a subject of debate. In the present study, phylogenetic analysis of Chinese native pigs was performed by screening for haplotypes inferred from a phylogenetic tree of pig mitochondrial DNA (mtDNA) sequences based on sequence-specific mutations. A total of 2466 domestic pigs formed 124 haplotypes and were assigned to four clades. Clade A comprised pigs distributed mainly in the Qinghai-Tibet Plateau and its surrounding areas; these pigs clustered into three groups. The pigs of clade B were mainly from the Mekong River Basin in Yunnan Province and had been exposed to genetic infiltration from European populations. Clade C comprised pigs mainly from the middle and lower reaches of the Yangtze River. The pigs of clade D were distributed mainly at the intersection of Yunnan, Sichuan, and Gansu provinces east of the Hengduan Mountains (YSGH). Compared with wild boar, at least three domestication centers and one expansion center of pigs in China were detected. Among the four centers detected, two were for Tibetan pigs and were in the Qinghai-Tibet Plateau and at the YSGH intersection, and the other two were in the Mekong River Basin in Yunnan Province and the middle and downstream regions of the Yangtze River.

GCH1 plays a role in the high-altitude adaptation of Tibetans.

Tibetans are well adapted to high-altitude hypoxia. Previous genome-wide scans have reported many candidate genes for this adaptation, but only a few have been studied. Here we report on a hypoxia gene ( GCH1, GTP-cyclohydrolase I), involved in maintaining nitric oxide synthetase (NOS) function and normal blood pressure, that harbors many potentially adaptive variants in Tibetans. We resequenced an 80.8 kb fragment covering the entire gene region of GCH1 in 50 unrelated Tibetans. Combined with previously published data, we demonstrated many GCH1 variants showing deep divergence between highlander Tibetans and lowlander Han Chinese. Neutrality tests confirmed a signal of positive Darwinian selection on GCH1 in Tibetans. Moreover, association analysis indicated that the Tibetan version of GCH1 was significantly associated with multiple physiological traits in Tibetans, including blood nitric oxide concentration, blood oxygen saturation, and hemoglobin concentration. Taken together, we propose that GCH1 plays a role in the genetic adaptation of Tibetans to high altitude hypoxia.

EP300 contributes to high-altitude adaptation in Tibetans by regulating nitric oxide production.

The genetic adaptation of Tibetans to high altitude hypoxia likely involves a group of genes in the hypoxic pathway, as suggested by earlier studies. To test the adaptive role of the previously reported candidate gene EP300 (histone acetyltransferase p300), we conducted resequencing of a 108.9 kb gene region of EP300 in 80 unrelated Tibetans. The allele-frequency and haplotype-based neutrality tests detected signals of positive Darwinian selection on EP300 in Tibetans, with a group of variants showing allelic divergence between Tibetans and lowland reference populations, including Han Chinese, Europeans, and Africans. Functional prediction suggested the involvement of multiple EP300 variants in gene expression regulation. More importantly, genetic association tests in 226 Tibetans indicated significant correlation of the adaptive EP300 variants with blood nitric oxide (NO) concentration. Collectively, we propose that EP300 harbors adaptive variants in Tibetans, which might contribute to high-altitude adaptation through regulating NO production.