Genetic and cultural adaptations underlie the establishment of dairy pastoralism in the Tibetan Plateau.

BACKGROUND: Domestication and introduction of dairy animals facilitated the permanent human occupation of the Tibetan Plateau. Yet the history of dairy pastoralism in the Tibetan Plateau remains poorly understood. Little is known how Tibetans adapted to milk and dairy products. RESULTS: We integrated archeological evidence and genetic analysis to show the picture that the dairy ruminants, together with dogs, were introduced from West Eurasia into the Tibetan Plateau since ~ 3600 years ago. The genetic admixture between the exotic and indigenous dogs enriched the candidate lactase persistence (LP) allele 10974A > G of West Eurasian origin in Tibetan dogs. In vitro experiments demonstrate that - 13838G > A functions as a LP allele in Tibetans. Unlike multiple LP alleles presenting selective signatures in West Eurasians and South Asians, the de novo origin of Tibetan-specific LP allele - 13838G > A with low frequency (~ 6-7%) and absence of selection corresponds - 13910C > T in pastoralists across eastern Eurasia steppe. CONCLUSIONS: Results depict a novel scenario of genetic and cultural adaptations to diet and expand current understanding of the establishment of dairy pastoralism in the Tibetan Plateau.

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.

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.

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.

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.

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.

Identification of the susceptible subpopulations for wide pulse pressure under long-term exposure to ambient particulate matters.

Wide pulse pressure (WPP) is a preclinical indicator for arterial stiffness and cardiovascular diseases. Long-term exposure to ambient particulate matters (PMs) would increase the risk of WPP. Although reducing pollutants emissions and avoiding outdoor activity during a polluted period are effective ways to blunt the adverse effects. Identifying and protecting the susceptible subpopulation is another crucial way to reduce the disease burdens. Therefore, we aimed to identify the susceptible subpopulations of WPP under long-term exposure to PMs. The WPP was defined as pulse pressure over 60 mmHg. Three-year averages of PMs were estimated using random forest approaches. Associations between WPP and PMs exposure were estimated using generalized propensity score weighted logistic regressions. Demographic, socioeconomic characteristics, health-related behaviors, and hematological biomarkers were collected to detect the modification effects on the WPP-PMs associations. Susceptible subpopulations were defined as those with significantly higher risks of WPP under PMs exposures. The PMs-WPP associations were significant with ORs (95%CI) of 1.126 (1.094, 1.159) for PM1, 1.174 (1.140, 1.210) for PM2.5, and 1.111 (1.088, 1.135) for PM10. There were 17 subpopulations more sensitive to WPP under long-term exposure to PMs. The susceptibility was higher in subpopulations with high BMI (Q3-Q4 quartiles), high-intensive physical activity (Q3 or Q4 quartile), insufficient or excessive fruit intake (Q1 or Q5 quartile), insufficient or too long sleep length (<7 or >8 h). Subpopulations with elevated inflammation markers (WBC, LYM, BAS, EOS: Q3-Q4 quartiles) and glucose metabolism indicators (HbA1c, GLU: Q3-Q4 quartiles) were more susceptible. Besides, elder, urban living, low socioeconomic level, and excessive red meat and sodium salt intake were also related to higher susceptibility. Our findings on the susceptibility characteristics would help to develop more targeted disease prevention and therapy strategies. Health resources can be allocated more effectively by putting more consideration to subpopulations with higher susceptibility.

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.

Associations of long-term exposure to ambient air pollution and physical activity with insomnia in Chinese adults.

BACKGROUND: Air pollution is a potential environmental risk for sleep disturbance. However, the evidence is very limited in China. On the other hand, physical activity (PA) is a preventive behavior that can improve insomnia, but whether PA mitigates the negative impact of air pollution on insomnia is unknown. METHODS: We obtained data from the baseline of China Multi-Ethnic Cohort (CMEC) survey, and examined the association between air pollution and insomnia, as well as PA's modification effect of on this association. We included 70,668 respondents and assessed insomnia by self-reported symptoms collected using electronic questionnaires. Using satellite data, we estimated the residence-specified, three-year average PM1, PM2.5, PM10 (particulate matter with aerodynamic diameters of ≤1 µm, ≤2.5 µm and 10 µm, respectively), O3 (ozone), and NO2 (nitrogen dioxide) concentrations. We established the associations between air pollutants and insomnia through logistic regression. We evaluated the modification impact of total and domain-specific PA (leisure, occupation, housework, transportation) by introducing an interaction term. RESULTS: Positive associations were observed between long-term exposure to PM1, PM2.5, PM10, and O3 and insomnia symptoms, with ORs (95% CI) of 1.09 (1.03-1.16), 1.11 (1.07-1.15), 1.07 (1.05-1.10) and 1.15 (1.11-1.20), respectively. As total PA increased, the ORs of air pollution for insomnia tended to decrease and then rise. We observed varying modification effects of domain-specific PA. With an increase in leisure PA, the ORs for PM2.5 and PM10 significantly declined. However, increased ORs of air pollutants were related to insomnia among participants with higher levels of occupational and housework PA. CONCLUSION: Long-term exposure to higher concentrations of PM1, PM2.5, PM10, and O3 increases the risk of insomnia symptoms. Moderate to high levels of leisure PA alleviate the harmful effects of air pollution on insomnia, while high levels of occupation and housework PA intensify such effects.

Chromatin accessibility landscape and regulatory network of high-altitude hypoxia adaptation.

High-altitude adaptation of Tibetans represents a remarkable case of natural selection during recent human evolution. Previous genome-wide scans found many non-coding variants under selection, suggesting a pressing need to understand the functional role of non-coding regulatory elements (REs). Here, we generate time courses of paired ATAC-seq and RNA-seq data on cultured HUVECs under hypoxic and normoxic conditions. We further develop a variant interpretation methodology (vPECA) to identify active selected REs (ASREs) and associated regulatory network. We discover three causal SNPs of EPAS1, the key adaptive gene for Tibetans. These SNPs decrease the accessibility of ASREs with weakened binding strength of relevant TFs, and cooperatively down-regulate EPAS1 expression. We further construct the downstream network of EPAS1, elucidating its roles in hypoxic response and angiogenesis. Collectively, we provide a systematic approach to interpret phenotype-associated noncoding variants in proper cell types and relevant dynamic conditions, to model their impact on gene regulation.

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.

The Transcriptomic Landscape of Yaks Reveals Molecular Pathways for High Altitude Adaptation.

Yak is one of the largest native mammalian species at the Himalayas, the highest plateau area in the world with an average elevation of >4,000 m above the sea level. Yak is well adapted to high altitude environment with a set of physiological features for a more efficient blood flow for oxygen delivery under hypobaric hypoxia. Yet, the genetic mechanism underlying its adaptation remains elusive. We conducted a cross-tissue, cross-altitude, and cross-species study to characterize the transcriptomic landscape of domestic yaks. The generated multi-tissue transcriptomic data greatly improved the current yak genome annotation by identifying tens of thousands novel transcripts. We found that among the eight tested tissues (lung, heart, kidney, liver, spleen, muscle, testis, and brain), lung and heart are two key organs showing adaptive transcriptional changes and >90% of the cross-altitude differentially expressed genes in lung display a nonlinear regulation. Pathways related to cell survival and proliferation are enriched, including PI3K-Akt, HIF-1, focal adhesion, and ECM-receptor interaction. These findings, in combination with the comprehensive transcriptome data set, are valuable to understanding the genetic mechanism of hypoxic adaptation in yak.

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.

Sherpas share genetic variations with Tibetans for high-altitude adaptation.

BACKGROUND: Sherpas, a highlander population living in Khumbu region of Nepal, are well known for their superior climbing ability in Himalayas. However, the genetic basis of their adaptation to high-altitude environments remains elusive. METHODS: We collected DNA samples of 582 Sherpas from Nepal and Tibetan Autonomous Region of China, and we measured their hemoglobin levels and degrees of blood oxygen saturation. We genotyped 29 EPAS1 SNPs, two EGLN1 SNPs and the TED polymorphism (3.4 kb deletion) in Sherpas. We also performed genetic association analysis among these sequence variants with phenotypic data. RESULTS: We found similar allele frequencies on the tested 32 variants of these genes in Sherpas and Tibetans. Sherpa individuals carrying the derived alleles of EPAS1 (rs113305133, rs116611511 and rs12467821), EGLN1 (rs186996510 and rs12097901) and TED have lower hemoglobin levels when compared with those wild-type allele carriers. Most of the EPAS1 variants showing significant association with hemoglobin levels in Tibetans were replicated in Sherpas. CONCLUSION: The shared sequence variants and hemoglobin trait between Sherpas and Tibetans indicate a shared genetic basis for high-altitude adaptation, consistent with the proposal that Sherpas are in fact a recently derived population from Tibetans and they inherited adaptive variants for high-altitude adaptation from their Tibetan ancestors.

Validity and reliability of the Tibetan version of s-EMBU for measuring parenting styles.

Parenting style experienced during childhood has profound effects on children's futures. Scales developed in other countries have never been validated in the Tibetan context. The present study aimed to examine the construct validity and reliability of a Tibetan translation of the 23-item short form of the Egna Minnen Beträffande Uppfostran [One's Memories of Upbringing] (s-EMBU) and to test the correlation between the parenting styles of fathers and mothers. A cross-sectional study was conducted in a sample of 847 students aged 12-21 years from Lhasa, Tibet, during September and October 2015 with a participation rate of 97.7%. The Tibetan translation of self-completed s-EMBU was administered. Confirmatory factor analysis was employed to test the scale's validity on the first half of the sample and was then cross-validated with the second half of the sample. The final model consisted of six factors: three (rejection, emotional warmth, and overprotection) for each parent, equality constrained on factor loadings, factor correlations, and error variance between father and mother. Father-mother correlation coefficients ranged from 0.81 to 0.86, and the level of consistency ranged from 0.62 to 0.82. Thus, the slightly modified s-EMBU is suitable for use in the Tibetan culture where both the father and the mother have consistent parenting styles.

Down-Regulation of EPAS1 Transcription and Genetic Adaptation of Tibetans to High-Altitude Hypoxia.

Tibetans are well adapted to the hypoxic environments at high altitude, yet the molecular mechanism of this adaptation remains elusive. We reported comprehensive genetic and functional analyses of EPAS1, a gene encoding hypoxia inducible factor 2α (HIF-2α) with the strongest signal of selection in previous genome-wide scans of Tibetans. We showed that the Tibetan-enriched EPAS1 variants down-regulate expression in human umbilical endothelial cells and placentas. Heterozygous EPAS1 knockout mice display blunted physiological responses to chronic hypoxia, mirroring the situation in Tibetans. Furthermore, we found that the Tibetan version of EPAS1 is not only associated with the relatively low hemoglobin level as a polycythemia protectant, but also is associated with a low pulmonary vasoconstriction response in Tibetans. We propose that the down-regulation of EPAS1 contributes to the molecular basis of Tibetans' adaption to high-altitude hypoxia.

HMOX2 Functions as a Modifier Gene for High-Altitude Adaptation in Tibetans.

Tibetans are well adapted to high-altitude environments. Among the adaptive traits in Tibetans, the relatively low hemoglobin level is considered a blunted erythropoietic response to hypoxic challenge. Previously, EPAS1 and EGLN1, the major upstream regulators in the hypoxic pathway, were reportedly involved in the hemoglobin regulation in Tibetans. In this study, we report a downstream gene (HMOX2) involved in heme catabolism, which harbors potentially adaptive variants in Tibetans. We first resequenced the entire genomic region (45.6 kb) of HMOX2 in Tibetans, which confirmed the previously suspected signal of positive selection on HMOX2 in Tibetans. Subsequent association analyses of hemoglobin levels in two independent Tibetan populations (a total of 1,250 individuals) showed a male-specific association between the HMOX2 variants and hemoglobin levels. Tibetan males with the derived C allele at rs4786504:T>C displayed lower hemoglobin level as compared with the T allele carriers. Furthermore, our in vitro experiments indicated that the C allele of rs4786504 could increase the expression of HMOX2, presumably leading to a more efficient breakdown of heme that may help maintain a relatively low hemoglobin level at high altitude. Collectively, we propose that HMOX2 contributes to high-altitude adaptation in Tibetans by functioning as a modifier in the regulation of hemoglobin metabolism.