Genetic signatures of high-altitude adaptation in Tibetans

Distinct methylomic signatures of high-altitude acclimatization and adaptation in the Tibetan Plateau

High altitude presents a challenging environment for human settlement. DNA methylation is an essential epigenetic mechanism that responds to environmental stimuli, but its roles in high-altitude short-term acclimatization (STA) and long-term adaptation (LTA) are poorly understood. Here, we conducted a methylome-wide association study involving 687 native highlanders and 299 acclimatized newcomers in the Tibetan Plateau and 462 native lowlanders to identify differentially methylated sites (DMSs) associated with STA or LTA. We identified 93 and 4070 DMSs for STA and LTA, respectively, which had no overlap, showed opposite asymmetric effect size patterns, and resided near genes enriched in distinct biological pathways/processes (e.g., cell cycle for STA and immune diseases and calcium signalling pathway for LTA). Epigenetic clock analysis revealed evidence of accelerated ageing in the acclimatized newcomers compared to the native lowlanders. Our research provides novel insights into epigenetic regulation in relation to high altitude and intervention strategies for altitude-related ageing or illnesses.

Life destiny of erythrocyte in high altitude erythrocytosis: mechanisms underlying the progression from physiological (moderate) to pathological (excessive) high-altitude erythrocytosis

High-altitude polycythemia (HAPC) represents a pathological escalation of the physiological erythrocytosis induced by chronic hypoxia exposure. While moderate erythroid expansion enhances oxygen delivery, HAPC manifests as hematologic disorder characterized by hemoglobin thresholds (≥21 g/dL males; ≥19 g/dL females) and multi-organ complications including microcirculatory thrombosis, right ventricular hypertrophy, and uric acid dysmetabolism. This review critically evaluates the continuum between adaptive and maladaptive polycythemia through multiscale analysis of erythrocyte biology. We integrate genomic predisposition patterns, bone marrow erythroid kinetic studies, and peripheral erythrocyte pathophenotypes revealed by multi-omics profiling (iron-redox proteome, hypoxia-metabolome crosstalk). Current diagnostic limitations are highlighted, particularly the oversimplification of hemoglobin cutoffs that neglect transitional dynamics in erythrocyte turnover. By reconstructing the erythroid life cycle-from hypoxia-sensitive progenitor commitment to senescent cell clearance-we propose a phase transition model where cumulative epigenetic-metabolic derangements overcome homeostatic buffers, triggering pathological erythroid amplification. These insights reframe HAPC as a systems biology failure of erythroid adaptation, informing predictive biomarkers and targeted interventions to preserve hematological homeostasis in hypoxic environments.

Genomic analysis reveals population structure and selection signatures in plateau dairy cattle

BACKGROUND

To solve the problem of an insufficient supply of dairy products in Tibet, work has been carried out to improve native dairy cattle and introduce purebred dairy cattle from low-altitude areas. The harsh environment of the plateau not only severely limits the production performance of high-yielding dairy cattle, such as Holstein and Jersey cattle, but also challenges their survival. The population structure and plateau adaptation mechanism of plateau dairy cattle are rarely reported. In this study, key genes and pathways affecting plateau purebred and crossbred dairy cattle were explored using genetic chip information.

RESULTS

The results showed that the genetic diversity of the Tibet dairy cattle population was higher than that of the native cattle and plains dairy cattle. Purebred Holstein and Jersey cattle in Tibet were genetically closer to dairy cattle in the plains, and crossbred dairy cattle were admixed with more Tibet cattle and Apaijiza cattle. Based on the fixation index (FST), integrated haplotype score (iHS), and cross-population extend haplotype homozygosity (XP-EHH) approaches, 60 and 40 genes were identified in plateau Holstein and Jersey cattle, respectively. A total of 78 and 70 genes were identified in crossbred cattle compared to Holstein and Tibet cattle respectively. These genes are related to cardiac health and development, neuronal development and function, angiogenesis and hematopoietic, pigmentation, growth and development, and immune response.

CONCLUSIONS

Our results provide a glimpse into diverse selection signatures in plateau dairy cattle, which can be used to enhance our understanding of the genomic basis of plateau adaptation in dairy cattle. These results support further research on breeding strategies such as marker-assisted selection and gene editing in plateau dairy cattle populations.

A highland-adaptation variant near MCUR1 reduces its transcription and attenuates erythrogenesis in Tibetans

To identify genomic regions subject to positive selection that might contain genes involved in high-altitude adaptation (HAA), we performed a genome-wide scan by whole-genome sequencing of Tibetan highlanders and Han lowlanders. We revealed a collection of candidate genes located in 30 genomic loci under positive selection. Among them, MCUR1 at 6p23 was a novel pronounced candidate. By single-cell RNA sequencing and comprehensive functional studies, we demonstrated that MCUR1 depletion leads to impairment of erythropoiesis under hypoxia and normoxia. Mechanistically, MCUR1 knockdown reduced mitochondrial Ca2+ uptake and then concomitantly increased cytosolic Ca2+ levels, which thereby reduced erythropoiesis via the CAMKK2-AMPK-mTOR axis. Further, we revealed rs61644582 at 6p23 as an expression quantitative trait locus for MCUR1 and a functional variant that confers an allele-specific transcriptional regulation of MCUR1. Overall, MCUR1-mediated mitochondrial Ca2+ homeostasis is highlighted as a novel regulator of erythropoiesis, deepening our understanding of the genetic mechanism of HAA.

Equitable machine learning counteracts ancestral bias in precision medicine

Gold standard genomic datasets severely under-represent non-European populations, leading to inequities and a limited understanding of human disease. Therapeutics and outcomes remain hidden because we lack insights that could be gained from analyzing ancestrally diverse genomic data. To address this significant gap, we present PhyloFrame, a machine learning method for equitable genomic precision medicine. PhyloFrame corrects for ancestral bias by integrating functional interaction networks and population genomics data with transcriptomic training data. Application of PhyloFrame to breast, thyroid, and uterine cancers shows marked improvements in predictive power across all ancestries, less model overfitting, and a higher likelihood of identifying known cancer-related genes. Validation in fourteen ancestrally diverse datasets demonstrates that PhyloFrame is better able to adjust for ancestry bias across all populations. The ability to provide accurate predictions for underrepresented groups, in particular, is substantially increased. Analysis of performance in the most diverse continental ancestry group, African, illustrates how phylogenetic distance from training data negatively impacts model performance, as well as PhyloFrame's capacity to mitigate these effects. These results demonstrate how equitable artificial intelligence (AI) approaches can mitigate ancestral bias in training data and contribute to equitable representation in medical research.

Gut Metagenome Reveals the Microbiome Signatures in Tibetan and Black Pigs

The harsh conditions of the Qinghai-Tibet Plateau pose significant physiological challenges to local fauna, often resulting in gastrointestinal disorders. However, Tibetan pigs have exhibited remarkable adaptability to the high-altitude stress of the Tibetan Plateau, a phenomenon that remains not fully understood in terms of their gastrointestinal microbiota. This study collected 57 gastrointestinal tract samples from Tibetan pigs (n = 6) and plain black pigs (n = 6) with comparable genetic backgrounds. Samples from the stomach, jejunum, cecum, colon, and rectum, underwent comprehensive metagenomic analysis to elucidate the gut microbiota-related adaptive mechanisms in Tibetan pigs to the extreme high-altitude environment. A predominance of Pseudomonadota was observed within gut microbiome of Tibetan pigs. Significant differences in the microbial composition were also identified across the tested gastrointestinal segments, with 18 genera and 141 species exhibiting differential abundance. Genera such as Bifidobacterium, Megasphaera, Fusobacterium, and Mitsuokella were significantly more abundant in Tibetan pigs than in their lowland counterparts, suggesting specialized adaptations. Network analysis found greater complexity and modularity in the microbiota of Tibetan pigs compared to black pigs, indicating enhanced ecological stability and adaptability. Functional analysis revealed that the Tibetan pig microbiota was particularly enriched with bacterial species involved in metabolic pathways for propionate and butyrate, key short-chain fatty acids that support energy provision under low-oxygen conditions. The enzymatic profiles of Tibetan pigs, characterized by elevated levels of 4-hydroxybutyrate dehydrogenase and glutaconyl-CoA decarboxylase, highlighted a robust fatty acid metabolism and enhanced tricarboxylic acid cycle activity. In contrast, the gut microbiome of plain black pigs showed a reliance on the succinate pathway, with a reduced butyrate metabolism and lower metabolic flexibility. Taken together, these results demonstrate the crucial role of the gastrointestinal microbiota in the adaptation of Tibetan pigs to high-altitude environments by optimizing carbohydrate metabolism and short-chain fatty acid production for efficient energy utilization. This study not only highlights the metabolic benefits conferred by the gut microbiota of Tibetan pigs in extreme environments, but also advances our understanding of the adaptive gastrointestinal mechanisms in plateau-dwelling animals. These insights lay the foundation for exploring metabolic interventions to support health and performance in high-altitude conditions.

Dissecting genomes of multiple yak populations: unveiling ancestry and high-altitude adaptation through whole-genome resequencing analysis

The present study was undertaken to elucidate the population structure and differentiation of Indian yak from Chinese and wild cohorts on genome-wide scale by identifying the selection sweeps and genomic basis of their adaptation across different comparisons while analyzing whole genome sequencing (WGS) data using latest bioinformatics tools. The study included 105 individuals from three distinct yak populations i.e., Indian yak (n = 29); Chinese yak (n = 61) and wild yak (n = 15), hypothesized to be related along the evolutionary timescale. Efficient variant calling and quality control in GATK and PLINK programs resulted in around 1 million (1,002,970) high-quality (LD-independent) SNPs with an average genotyping rate of 96.55%. The PCA, ADMIXTURE and TREEMIX analysis revealed stratification of the yak groups into three distinct clusters. The empirical distribution pattern of minor allele frequency (MAF) of SNPs on genome-wide scale was also elucidated for three yak cohorts revealing unique distribution across five different bins. The selection signature analysis revealed candidate genes that are important for the adaptation of Indian yak against harsh environmental conditions in their habitats. Under iHS analysis, several genes were identified to be under selection pressure in Indian yak including ABCA12, EXOC1, JUNB, KLF1, PRDX2, NANOS3, RFX1, RFX2, and CACNG7. On the other hand, across population analysis revealed the genes like NR2F2, OSBPL10, CIDEC, WFIKKN2, ADCY, THSD7A, ADGRB3, TRPC1, VASH2, and ABHD5 to be part of selective sweeps under these comparisons. A total of 53 genes were found common between intra- and inter-population selection signature analysis of Indian yak. Notably, the genes harbouring the SNPs under selection pressure were significant for adaptation traits including lipidogenesis, energy metabolism, thermogenesis, hair follicle formation, oxidation-reduction reactions, hypoxia and reproduction. These genes may be evaluated as candidate genes for livestock adaptation to harsh environmental conditions and to further the research and application in the present era of climate change.

Human dispersal into East Eurasia: ancient genome insights and the need for research on physiological adaptations

Humans have long pondered their genesis. The answer to the great question of where Homo sapiens come from has evolved in conjunction with biotechnologies that have allowed us to more brightly illuminate our distant past. The "Multiregional Evolution" model was once the hegemonic theory of Homo sapiens origins, but in the last 30 years, it has been supplanted by the "Out of Africa" model. Here, we review the major findings that have resulted in this paradigmatic shift. These include hominin brain expansion, classical insight from the mitochondrial genome (mtDNA) regarding the timing of the divergence point between Africans and non-Africans, and next-generation sequencing (NGS) of the Neanderthal and Denisovan genomes. These findings largely bolstered the "Out of Africa" model, although they also revealed a small degree of introgression of the Neanderthal and Denisovan genomes into those of non-African Homo sapiens. We also review paleogenomic studies for which migration route, north or south, early migrants to East Eurasia most likely traversed. Whichever route was taken, the migrants moved to higher latitudes, which necessitated adaptation for lower light conditions, colder clines, and pro-adipogenic mechanisms to counteract food scarcity. Further genetic and epigenetic investigations of these physiological adaptations constitute an integral aspect of the story of human origins and human migration to East Asia.

Transcriptomic analysis of iPSC-derived endothelium reveals adaptations to high altitude hypoxia in energy metabolism and inflammation

Tibetan adaptation to high-altitude hypoxia remains a classic example of Darwinian selection in humans. Amongst Tibetan populations, alleles in the EPAS1 gene - whose protein product, HIF-2α, is a central regulator of the hypoxia response - have repeatedly been shown to carry some of the strongest signals of positive selection in humans. However, selective sweep signals alone may only account for some of the phenotypes that differentiate high-altitude adapted populations from closely related lowlanders. Therefore, there is a pressing need to functionally probe adaptive alleles and their impact at both the locus-specific and genome-wide levels and across cell types to uncover the full range of beneficial traits. To this end, we established a library of induced pluripotent stem cells (iPSCs) derived from Tibetan and Han Chinese individuals, a robust model system allowing precise exploration of allelic effects on transcriptional responses, and we differentiated them into vascular endothelium. Using this system, we focus first on a hypoxia-dependent enhancer (ENH5) that contributes to the regulation of EPAS1 to investigate its locus-specific effects in endothelium. Then, to cast a wider net, we harness the same experimental system to compare the transcriptome of Tibetan and Han Chinese cells in hypoxia and find evidence that angiogenesis, energy metabolism and immune pathways differ between these two populations with different histories of long-term residence at high altitude. Coupled with evidence of polygenic adaptations targeting the same pathways, these results suggests that the observed transcriptional differences between the two populations were shaped by natural selection.

Selection on the vascular-remodeling BMPER gene is associated with altitudinal adaptation in an insular lizard

High altitude imposes several extreme constraints on life, such as low oxygen pressure and high levels of ultraviolet radiation, which require specialized adaptations. Many studies have focused on how endothermic vertebrates respond to these challenging environments, but there is still uncertainty on how ectotherms adapt to these conditions. Here, we used whole-genome sequencing of low-altitude (100-600 m) and high-altitude (3,550 m) populations of the wide-ranging Tenerife lizard Gallotia galloti to uncover signatures of selection for altitudinal adaptation. The studied populations show reduced differentiation, sharing similar patterns of genetic variation. Selective sweep mapping suggests that signatures of adaptation to high altitude are not widespread across the genome, clustering in a relatively small number of genomic regions. One of these regions contains BMPER, a gene involved with vascular remodeling, and that has been associated with hypoxia-induced angiogenic response. By genotyping samples across 2 altitudinal transects, we show that allele frequency changes at this locus are not gradual, but rather show a well-defined shift above ca. 1,900 m. Transcript and protein structure analyses on this gene suggest that putative selection likely acts on noncoding variation. These results underline how low oxygen pressure generates the most consistent selective constraint in high-altitude environments, to which vertebrates with vastly contrasting physiological profiles need to adapt in the context of ongoing climate change.

Variations in HBA gene contribute to high-altitude hypoxia adaptation via affected O2 transfer in Tibetan sheep

Tibetan sheep are indigenous to the Qinghai-Xizang Plateau. Owing to the harsh hypoxic environment in this plateau, the hemoglobin (Hb) protein in Tibetan sheep has undergone adaptive changes over time. Hb is primarily responsible for transporting O2 and CO2 between the lungs and other tissues of the body. The α subunit of Hb, encoded by the HBA gene, is a crucial component of the protein. However, whether variations in the HBA gene sequence affect the adaptation of Tibetan sheep to high-altitude hypoxia remains unclear. In this study, we sequenced the HBA gene and identified three single nucleotide polymorphisms (SNPs). These SNPs were genotyped in Tibetan and Hu sheep using Kompetitive Allele-Specific PCR (KASP). The results showed that the frequencies of the AT genotype and H1H2 haplotype were higher in Tibetan sheep than in Hu sheep. Individuals with the AT genotype exhibited higher P50 levels, whereas those with the H1H2 haplotype exhibited lower PO2 and SaO2 levels. The higher P50 levels indicated that O2 was more readily released from oxygenated Hb into the tissues, with the lower PO2 and SaO2 levels facilitating this process. These findings indicate that variations in the HBA gene sequence contribute to enhancing O2 transfer efficiency in Tibetan sheep.

Archaic introgression contributed to shape the adaptive modulation of angiogenesis and cardiovascular traits in human high-altitude populations from the Himalayas

It is well established that several Homo sapiens populations experienced admixture with extinct human species during their evolutionary history. Sometimes, such a gene flow could have played a role in modulating their capability to cope with a variety of selective pressures, thus resulting in archaic adaptive introgression events. A paradigmatic example of this evolutionary mechanism is offered by the EPAS1 gene, whose most frequent haplotype in Himalayan highlanders was proved to reduce their susceptibility to chronic mountain sickness and to be introduced in the gene pool of their ancestors by admixture with Denisovans. In this study, we aimed at further expanding the investigation of the impact of archaic introgression on more complex adaptive responses to hypobaric hypoxia evolved by populations of Tibetan/Sherpa ancestry, which have been plausibly mediated by soft selective sweeps and/or polygenic adaptations rather than by hard selective sweeps. For this purpose, we used a combination of composite-likelihood and gene network-based methods to detect adaptive loci in introgressed chromosomal segments from Tibetan WGS data and to shortlist those presenting Denisovan-like derived alleles that participate to the same functional pathways and are absent in populations of African ancestry, which are supposed to do not have experienced Denisovan admixture. According to this approach, we identified multiple genes putatively involved in archaic introgression events and that, especially as regards TBC1D1, RASGRF2, PRKAG2, and KRAS, have plausibly contributed to shape the adaptive modulation of angiogenesis and of certain cardiovascular traits in high-altitude Himalayan peoples. These findings provided unprecedented evidence about the complexity of the adaptive phenotype evolved by these human groups to cope with challenges imposed by hypobaric hypoxia, offering new insights into the tangled interplay of genetic determinants that mediates the physiological adjustments crucial for human adaptation to the high-altitude environment.

Higher oxygen content and transport characterize high-altitude ethnic Tibetan women with the highest lifetime reproductive success

We chose the "natural laboratory" provided by high-altitude native ethnic Tibetan women who had completed childbearing to examine the hypothesis that multiple oxygen delivery traits were associated with lifetime reproductive success and had genomic associations. Four hundred seventeen (417) women aged 46 to 86 y residing at ≥3,500 m in Upper Mustang, Nepal, provided information on reproductive histories, sociocultural factors, physiological measurements, and DNA samples for this observational cohort study. Simultaneously assessing multiple traits identified combinations associated with lifetime reproductive success measured as the number of livebirths. Women with the most livebirths had distinctive hematological and cardiovascular traits. A hemoglobin concentration near the sample mode and a high percent of oxygen saturation of hemoglobin raised arterial oxygen concentration without risking elevated blood viscosity. We propose ongoing stabilizing selection on hemoglobin concentration because extreme values predicted fewer livebirths and directional selection favoring higher oxygen saturation because higher values had more predicted livebirths. EPAS1, an oxygen homeostasis locus with strong signals of positive natural selection and a high frequency of variants occurring only among populations indigenous to the Tibetan Plateau, associated with hemoglobin concentration. High blood flow into the lungs, wide left ventricles, and low hypoxic heart rate responses aided effective convective oxygen transport to tissues. Women with physiologies closer to unstressed, low altitude values had the highest lifetime reproductive success. This example of ethnic Tibetan women residing at high altitudes in Nepal links reproductive fitness with trait combinations increasing oxygen delivery under severe hypoxic stress and demonstrates ongoing natural selection.

Whole Genome Analysis Reveals Evolutionary History and Introgression Events in Bale Monkeys

Background/Objective: The Bale monkey (Chlorocebus djamdjamensis) is a threatened primate species endemic to Ethiopia and, in contrast to other members of the genus Chlorocebus, lives at high altitudes and feeds mainly on bamboo. Two populations of the species are present, one in continuous bamboo forest (CF) in the eastern part of the species' range, and the other in fragmented forest (FF) in the western part. Based on mitochondrial DNA and phenotypic characteristics, previous studies have suggested introgression by parapatric congeners into the FF population but not into the CF population. The objective of this study was to gain insights into the evolutionary history of Bale monkeys and their potential genetic adaptations to high altitudes and for bamboo consumption. Methods: We sequenced the whole genomes of individuals from both populations and compared their genomes with those of the other five Chlorocebus species. We applied phylogenetic methods and conducted population demographic simulations to elucidate their evolutionary history. A genome-wide analysis was conducted to assess gene flow and identify mutations potentially associated with adaptations to high altitudes and for bamboo metabolism. Results: Our analyses revealed Bale monkeys as the sister clade to Chlorocebus aethiops and showed that gene flow occurred between C. aethiops and FF but not between C. aethiops and CF. In addition, we detected non-synonymous mutations in genes potentially associated with the adaptation to high altitudes (EPAS1) in both populations and with the adaptation for bamboo metabolism (TAS2R16, MPST, and TST) mainly in the CF population. Conclusions: Our study provides insights into the evolutionary history of a threatened primate species and reveals the genetic basis for its adaptions to unique environments and for diet specialization.

Physiological and Genetic Basis of High-Altitude Indigenous Animals' Adaptation to Hypoxic Environments

Adaptation is one of the fundamental characteristics of life activities; humans and animals inhabiting high altitudes are well adapted to hypobaric hypoxic environments, and studies on the mechanisms of this adaptation emerged a hundred years ago. Based on these studies, this paper reviews the adaptive changes in hypoxia-sensitive tissues and organs, as well as at the molecular genetic level, such as pulmonary, cardiovascular, O2-consuming tissues, and the hemoglobin and HIF pathway, that occur in animals in response to the challenge of hypobaric hypoxia. High-altitude hypoxia adaptation may be due to the coordinated action of genetic variants in multiple genes and, as a result, adaptive changes in multiple tissues and organs at the physiological and biochemical levels. Unraveling their mechanisms of action can provide a reference for the prevention and treatment of multiple diseases caused by chronic hypoxia.

Ancestral origins and post-admixture adaptive evolution of highland Tajiks

It remains debatable how many genes and how various the mechanisms are behind human adaptation to extreme environments, such as high altitudes. Despite extensive studies on Tibetans, Andeans and Ethiopians, new insights are expected to be provided with careful analysis of underrepresented highlanders living in a different geographical region, such as the Tajiks, who reside on the Pamir Plateau at an average altitude exceeding 4000 meters. Moreover, genetic admixture, as we observed in the current whole-genome deep-sequencing study of Xinjiang Tajiks (XJT), offers a unique opportunity to explore how admixture may facilitate adaptation to high-altitude environments. Compared with other extensively studied highlanders, XJT showed pronounced admixture patterns: most of their ancestry are derived from West Eurasians (34.5%-48.3%) and South Asians (21.4%-40.0%), and some minor ancestry from East Asians and Siberians (3.62%-17.5%). The greater genetic diversity in XJT than in their ancestral source populations provides a genetic basis for their adaptation to high-altitude environments. The admixture gain of functional adaptive components from ancestral populations could facilitate adaptation to high-altitude environments. Specifically, admixture-facilitated adaptation was strongly associated with skin-related candidate genes that respond to UV radiation (e.g. HERC2 and BNC2) and cardiovascular-system-related genes (e.g. MPI and BEST1). Notably, no adaptive variants of genes showing outstanding natural selection signatures in the Tibetan or Andean highlanders were identified in XJT, including EPAS1 and EGLN1, indicating that a different set of genes contributed to XJT's survival on the Pamir Plateau, although some genes underlying natural selection in XJT have been previously reported in other highlanders. Our results highlight the unique genetic adaptations in XJT and propose that admixture may play a vital role in facilitating high-altitude adaptation. By introducing and elevating diversity, admixture likely induces novel genetic factors that contribute to the survival of populations in extreme environments like the highlands.

Protective alleles and precision healthcare in crewed spaceflight

Common and rare alleles are now being annotated across millions of human genomes, and omics technologies are increasingly being used to develop health and treatment recommendations. However, these alleles have not yet been systematically characterized relative to aerospace medicine. Here, we review published alleles naturally found in human cohorts that have a likely protective effect, which is linked to decreased cancer risk and improved bone, muscular, and cardiovascular health. Although some technical and ethical challenges remain, research into these protective mechanisms could translate into improved nutrition, exercise, and health recommendations for crew members during deep space missions.

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.

Novel insight into the genetic signatures of altitude adaptation related body composition in Tibetans

Background

The Tibetan population residing in high-altitude (HA) regions has adapted to extreme hypoxic environments. However, there is limited understanding of the genetic basis of body compositions in Tibetan population adapted to HA.

Methods

We performed a genome-wide association study (GWAS) to identify genetic variants associated with HA and HA-related body composition traits. A total of 755,731 single nucleotide polymorphisms (SNPs) were genotyped using the precision medicine diversity array from 996 Tibetan college students. T-tests and Pearson correlation analysis were used to estimate the association between body compositions and altitude. The mixed linear regression identified the SNPs significantly associated with HA and HA-related body compositions. LASSO regression was used to screen for important SNPs in HA and body compositions.

Results

Significant differences were observed in lean body mass (LBW), muscle mass (MM), total body water (TBW), standard weight (SBW), basal metabolic rate (BMR), total protein (TP), and total inorganic salt (Is) in different altitudes stratification. We identified three SNPs in EPAS1 (rs1562453, rs7589621 and rs7583392) that were significantly associated with HA (p < 5 × 10-7). GWAS analysis of 7 HA-related body composition traits, we identified 14 SNPs for LBM, 11 SNPs for TBW, 15 SNPs for MM, 16 SNPs for SBW, 9 SNPs for BMR, 12 SNPs for TP, and 26 SNPs for Is (p < 5.0 × 10-5).

Conclusion

These findings provide insight into the genetic basis of body composition in Tibetan college students adapted to HA, and lay the foundation for further investigation into the molecular mechanisms underlying HA adaptation.

Polymorphic transposable elements provide new insights on high-altitude adaptation in the Tibetan Plateau

Several studies demonstrated that populations living in the Tibetan plateau are genetically and physiologically adapted to high-altitude conditions, showing genomic signatures ascribable to the action of natural selection. However, so far most of them relied solely on inferences drawn from the analysis of coding variants and point mutations. To fill this gap, we focused on the possible role of polymorphic transposable elements in influencing the adaptation of Tibetan and Sherpa highlanders. To do so, we compared high-altitude and middle/low-lander individuals of East Asian ancestry by performing in silico analyses and differentiation tests on 118 modern and ancient samples. We detected several transposable elements associated with high altitude, which map genes involved in cardiovascular, hematological, chem-dependent and respiratory conditions, suggesting that metabolic and signaling pathways taking part in these functions are disproportionately impacted by the effect of environmental stressors in high-altitude individuals. To our knowledge, our study is the first hinting to a possible role of transposable elements in the adaptation of Tibetan and Sherpa highlanders.

Immune consequences of exercise in hypoxia: A narrative review

Immune outcomes are key mediators of many health benefits of exercise and are determined by exercise type, dose (frequency/duration, intensity), and individual characteristics. Similarly, reduced availability of ambient oxygen (hypoxia) modulates immune functions depending on the hypoxic dose and the individual capacity to respond to hypoxia. How combined exercise and hypoxia (e.g., high-altitude training) sculpts immune responses is not well understood, although such combinations are becoming increasingly popular. Therefore, in this paper, we summarize the impact on immune responses of exercise and of hypoxia, both independently and together, with a focus on specialized cells in the innate and adaptive immune system. We review the regulation of the immune system by tissue oxygen levels and the overlapping and distinct immune responses related to exercise and hypoxia, then we discuss how they may be modulated by nutritional strategies. Mitochondrial, antioxidant, and anti-inflammatory mechanisms underlie many of the adaptations that can lead to improved cellular metabolism, resilience, and overall immune functions by regulating the survival, differentiation, activation, and migration of immune cells. This review shows that exercise and hypoxia can impair or complement/synergize with each other while regulating immune system functions. Appropriate acclimatization, training, and nutritional strategies can be used to avoid risks and tap into the synergistic potentials of the poorly studied immune consequences of exercising in a hypoxic state.

Genetic control of DNA methylation is largely shared across European and East Asian populations

DNA methylation is an ideal trait to study the extent of the shared genetic control across ancestries, effectively providing hundreds of thousands of model molecular traits with large QTL effect sizes. We investigate cis DNAm QTLs in three European (n = 3701) and two East Asian (n = 2099) cohorts to quantify the similarities and differences in the genetic architecture across populations. We observe 80,394 associated mQTLs (62.2% of DNAm probes with significant mQTL) to be significant in both ancestries, while 28,925 mQTLs (22.4%) are identified in only a single ancestry. mQTL effect sizes are highly conserved across populations, with differences in mQTL discovery likely due to differences in allele frequency of associated variants and differing linkage disequilibrium between causal variants and assayed SNPs. This study highlights the overall similarity of genetic control across ancestries and the value of ancestral diversity in increasing the power to detect associations and enhancing fine mapping resolution.

Patient Sex and Origin Influence Distribution of Driver Genes and Clinical Presentation of Paraganglioma

Context

Sexual and ancestral differences in driver gene prevalence have been described in many cancers but have not yet been investigated in pheochromocytoma and paraganglioma (PPGL).

Objective

This study aims to assess whether sex and ancestry influence prevalence of PPGL driver genes and clinical presentation.

Methods

We conducted a retrospective analysis of patients with PPGL considering studies from 2010 onwards that included minimal data of type of disease, sex, mutated gene, and country of origin. Additional features were recorded when available (age, tumor location, bilateral or multifocal, somatic or germline, and metastatic disease).

Results

We included 2162 patients: 877 in Europe and 757 in Asia. Males presented more often with germline pathogenic variants (PVs) in genes activating hypoxia pathways (P = .0006) and had more often sympathetic paragangliomas (P = .0005) and metastasis (P = .0039). On the other hand, females with PPGLs due to MAX PVs were diagnosed later than males (P = .0378) and more often developed metastasis (P = .0497). European but not Asian females presented more often with PPGLs due to PVs in genes related to kinase signaling (P = .0052), particularly RET and TMEM127. Contrary to experiences from Europe, Asian patients with PPGL due to PVs in kinase signaling genes NF1, HRAS, and FGFR1 showed a high proportion of sympathetic tumors, while European patients almost exclusively had adrenal tumors (P < .005).

Conclusion

Personalized management of patients with PPGL might benefit from considering sexual and ancestral differences. Further studies with better clinically aligned cohorts from various origins are required to better dissect ancestral influences on PPGL development.

Evolution and emergence of Mycobacterium tuberculosis

Tuberculosis (TB) remains one of the deadliest infectious diseases in human history, prevailing even in the 21st century. The causative agents of TB are represented by a group of closely related bacteria belonging to the Mycobacterium tuberculosis complex (MTBC), which can be subdivided into several lineages of human- and animal-adapted strains, thought to have shared a last common ancestor emerged by clonal expansion from a pool of recombinogenic Mycobacterium canettii-like tubercle bacilli. A better understanding of how MTBC populations evolved from less virulent mycobacteria may allow for discovering improved TB control strategies and future epidemiologic trends. In this review, we highlight new insights into the evolution of mycobacteria at the genus level, describing different milestones in the evolution of mycobacteria, with a focus on the genomic events that have likely enabled the emergence and the dominance of the MTBC. We also review the recent literature describing the various MTBC lineages and highlight their particularities and differences with a focus on host preferences and geographic distribution. Finally, we discuss on putative mechanisms driving the evolution of tubercle bacilli and mycobacteria in general, by taking the mycobacteria-specific distributive conjugal transfer as an example.

An Anthropometric Study of the Morphologic Facial Index of Tibetan Youth in Tibet

The purpose of this study was to understand sex differences and variations in facial indices among Tibetans and to create and evaluate anthropometric data on facial morphology. The study population consisted of 476 native Tibetans (241 males and 235 females) aged 18 to 24 years. The means and SD facial width was 133.53±7.31 mm for males and 133.95±8.10 mm for females; the difference between the sexes was not statistically significant. The means and SD facial height was 107.68±5.76 mm for males and 111.95±14.28 mm for females; the difference between the sexes was statistically significant (u=-8.394, P=0.000). The morphologic facial index was 80.86±5.82 (means±SD) for males and 83.91±11.90 (means±SD) for females; the difference between the sexes was statistically significant (u=-6.581, P=0.000). The proportion of the Tibetan male facial shape was hypereuryprosopic (45.6%) > euryprosopic (31.1%) > mesoprosopic (18.7%) > leptoprosopic (3.3%) > hyperleptoprosopic (1.2%). The proportion of the Tibetan female facial shape was hypereuryprosopic (25.5%) > mesoprosopic (22.6%) > euryprosopic (21.7%) > leptoprosopic (17.4%) > hyperleptoprosopic (12.8%). Facial width was positively correlated with height (male r=0.306, P=0.000; female r=0.144, P=0.027), weight (r=0.470, P=0.000 for males; r=0.337, P=0.000 for females), and BMI (r=0.378, P=0.000 for males; r=0.291, P=0.000 for females). Facial height was positively correlated with height (r=0.329, P=0.000 for males; r=0.137, P=0.035 for females) and weight (r=0.391, P=0.000 for males; r=0.170, P=0.009 for females). Facial height was positively correlated with BMI in Tibetan males (r=0.293, P=0.000), but no significant correlation was found in Tibetan females. The morphologic facial index of Tibetans was positively correlated with age (r=0.183, P=0.004 for males; r=0.171, P=0.009 for females). The results indicated that Tibetan youth in Tibet have a predominantly hypereuryprosopic facial shape and that facial features are related to age, height, and weight. Some common facial morphology features exist among the Tibet Tibetans, northeastern Indians, and Nepalese in the 3 different regions of the Sino-Tibetan language family. The data from this study provide basic information for the study of Tibetans in the fields of physical anthropology, forensic medicine, maxillofacial surgery, and plastic surgery.

Genomic structural variation is associated with hypoxia adaptation in high-altitude zokors

Zokors, an Asiatic group of subterranean rodents, originated in lowlands and colonized high-elevational zones following the uplift of the Qinghai-Tibet plateau about 3.6 million years ago. Zokors live at high elevation in subterranean burrows and experience hypobaric hypoxia, including both hypoxia (low oxygen concentration) and hypercapnia (elevated partial pressure of CO2). Here we report a genomic analysis of six zokor species (genus Eospalax) with different elevational ranges to identify structural variants (deletions and inversions) that may have contributed to high-elevation adaptation. Based on an assembly of a chromosome-level genome of the high-elevation species, Eospalax baileyi, we identified 18 large inversions that distinguished this species from congeners native to lower elevations. Small-scale structural variants in the introns of EGLN1, HIF1A, HSF1 and SFTPD of E. baileyi were associated with the upregulated expression of those genes. A rearrangement on chromosome 1 was associated with altered chromatin accessibility, leading to modified gene expression profiles of key genes involved in the physiological response to hypoxia. Multigene families that underwent copy-number expansions in E. baileyi were enriched for autophagy, HIF1 signalling and immune response. E. baileyi show a significantly larger lung mass than those of other Eospalax species. These findings highlight the key role of structural variants underlying hypoxia adaptation of high-elevation species in Eospalax.

Deciphering climate resilience in Indian cattle breeds by selection signature analyses

The signature of selection is a crucial concept in evolutionary biology that refers to the pattern of genetic variation which arises in a population due to natural selection. In the context of climate adaptation, the signature of selection can reveal the genetic basis of adaptive traits that enable organisms to survive and thrive in changing environmental conditions. Breeds living in diverse agroecological zones exhibit genetic "footprints" within their genomes that mirror the influence of climate-induced selective pressures, subsequently impacting phenotypic variance. It is assumed that the genomes of animals residing in these regions have been altered through selection for various climatic adaptations. These regions are known as signatures of selection and can be identified using various summary statistics. We examined genotypic data from eight different cattle breeds (Gir, Hariana, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, and Tharparkar) that are adapted to diverse regional climates. To identify selection signature regions in this investigation, we used four intra-population statistics: Tajima's D, CLR, iHS, and ROH. In this study, we utilized Bovine 50 K chip data and four genome scan techniques to assess the genetic regions of positive selection for high-temperature adaptation. We have also performed a genome-wide investigation of genetic diversity, inbreeding, and effective population size in our target dataset. We identified potential regions for selection that are likely to be caused by adverse climatic conditions. We observed many adaptation genes in several potential selection signature areas. These include genes like HSPB2, HSPB3, HSP20, HSP90AB1, HSF4, HSPA1B, CLPB, GAP43, MITF, and MCHR1 which have been reported in the cattle populations that live in varied climatic regions. The findings demonstrated that genes involved in disease resistance and thermotolerance were subjected to intense selection. The findings have implications for marker-assisted breeding, understanding the genetic landscape of climate-induced adaptation, putting breeding and conservation programs into action.

Hypoxia Inducible Factor pathway proteins in high-altitude mammals

Humans and other mammals inhabit hypoxic high-altitude locales. In many of these species, genes under positive selection include ones in the Hypoxia Inducible Factor (HIF) pathway. One is PHD2 (EGLN1), which encodes for a key oxygen sensor. Another is HIF2A (EPAS1), which encodes for a PHD2-regulated transcription factor. Recent studies have provided insights into mechanisms for these high-altitude alleles. These studies have (i) shown that selection can occur on nonconserved, unstructured regions of proteins, (ii) revealed that high altitude-associated amino acid substitutions can have differential effects on protein-protein interactions, (iii) provided evidence for convergent evolution by different molecular mechanisms, and (iv) suggested that mutations in different genes can complement one another to produce a set of adaptive phenotypes.

Meta-analysis identifying gut microbial biomarkers of Qinghai-Tibet Plateau populations and the functionality of microbiota-derived butyrate in high-altitude adaptation

The extreme environmental conditions of a plateau seriously threaten human health. The relationship between gut microbiota and human health at high altitudes has been extensively investigated. However, no universal gut microbiota biomarkers have been identified in the plateau population, limiting research into gut microbiota and high-altitude adaptation. 668 16s rRNA samples were analyzed using meta-analysis to reduce batch effects and uncover microbiota biomarkers in the plateau population. Furthermore, the robustness of these biomarkers was validated. Mendelian randomization (MR) results indicated that Tibetan gut microbiota may mediate a reduced erythropoietic response. Functional analysis and qPCR revealed that butyrate may be a functional metabolite in high-altitude adaptation. A high-altitude rat model showed that butyrate reduced intestinal damage caused by high altitudes. According to cell experiments, butyrate may downregulate hypoxia-inducible factor-1α (HIF-1α) expression and blunt cellular responses to hypoxic stress. Our research found universally applicable biomarkers and investigated their potential roles in promoting human health at high altitudes.

Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside

Rhodiola L. is a genus that has undergone rapid radiation in the mid-Miocene and may represent a typic case of adaptive radiation. Many species of Rhodiola have also been widely used as an important adaptogen in traditional medicines for centuries. However, a lack of high-quality chromosome-level genomes hinders in-depth study of its evolution and biosynthetic pathway of secondary metabolites. Here, we assembled two chromosome-level genomes for two Rhodiola species with different chromosome number and sexual system. The assembled genome size of R. chrysanthemifolia (2n = 14; hermaphrodite) and R. kirilowii (2n = 22; dioecious) were of 402.67 and 653.62 Mb, respectively, with approximately 57.60% and 69.22% of transposable elements (TEs). The size difference between the two genomes was mostly due to proliferation of long terminal repeat-retrotransposons (LTR-RTs) in the R. kirilowii genome. Comparative genomic analysis revealed possible gene families responsible for high-altitude adaptation of Rhodiola, including a homolog of plant cysteine oxidase 2 gene of Arabidopsis thaliana (AtPCO2), which is part of the core molecular reaction to hypoxia and contributes to the stability of Group VII ethylene response factors (ERF-VII). We found extensive chromosome fusion/fission events and structural variations between the two genomes, which might have facilitated the initial rapid radiation of Rhodiola. We also identified candidate genes in the biosynthetic pathway of salidroside. Overall, our results provide important insights into genome evolution in plant rapid radiations, and possible roles of chromosome fusion/fission and structure variation played in rapid speciation.

Longitudinal study on blood and biochemical indexes of Tibetan and Han in high altitude area

Objective

This study aims to review the blood routine and biochemical indicators of the plateau population for three consecutive years, and analyze the impact of the plateau on these blood indicators of the Tibetan population and the Han immigrant population.

Method

These parameters were extracted from the Laboratory Department of Ali District People's Hospital in Tibet from January 2019 to December 2021, including blood routine, liver and kidney function, blood lipids, myocardial enzyme spectrum, and rheumatic factor indicators. Changes in these parameters were analyzed over 3 consecutive years according to inclusion and exclusion criteria.

Result

A total of 114 Tibetans and 93 Hans participated in the study. These parameters were significantly different between Tibetan and Han populations. Red blood cells (RBC), hemoglobin (HGB), hematocrit (HCT), mean hemoglobin content (MCH), mean corpuscular hemoglobin concentration (MCHC), white blood cells (WBC), lymphocytes (LYMPH) and monocytes (MONO) were significantly higher in Hans than Tibetans (p < 0.05). Biochemically, total bilirubin (TBIL), direct bilirubin (DBIL), albumin (ALB), urea nitrogen (Urea), creatinine (Cr), uric acid (UA), glucose (GLU), triglycerides (TG) and creatine kinase isoenzyme (CKMB) were significantly higher in Hans than Tibetans; aspartate aminotransferase (AST), glutamyl transpeptidase (GGT), alkaline phosphatase (ALP), antistreptolysin (ASO), and C-reactive protein (CRP) were significantly higher in Tibetans than Hans (p < 0.05). There were no obvious continuous upward or downward trend of the parameters for 3 consecutive years.

Conclusion

In high-altitude areas, Han immigrants have long-term stress changes compared with Tibetans. The main differences are reflected in the blood system, liver and kidney functions, etc., which provide basic data for further research on the health status of plateau populations.

Human genetics: An extreme fitness landscape

A new study aims to identify how genetic and physiological adaptations to altitude affect pregnancy, childbirth and neonatal health in one of the most extreme environments on Earth, the Tibetan Plateau.

Genetic Origins and Adaptive Evolution of the Deng People on the Tibetan Plateau

The Tibetan Plateau is populated by diverse ethnic groups, but most of them are underrepresented in genomics studies compared with the Tibetans (TIB). Here, to gain further insight into the genetic diversity and evolutionary history of the people living in the Tibetan Plateau, we sequenced 54 whole genomes of the Deng people with high coverage (30-60×) and analyzed the data together with that of TIB and Sherpas, as well as 968 ancient Asian genomes and available archaic and modern human data. We identified 17.74 million novel single-nucleotide variants from the newly sequenced genomes, although the Deng people showed reduced genomic diversity and a relatively small effective population size. Compared with the other Tibetan highlander groups which are highly admixed, the Deng people are dominated by a sole ancestry that could be traced to some ancient northern East Asian populations. The divergence between Deng and Tibetan people (∼4,700-7,200 years) was more recent than that between highlanders and the Han Chinese (Deng-HAN, ∼9,000-14,000 years; TIB-HAN, 7,200-10,000 years). Adaptive genetic variants (AGVs) identified in the Deng are only partially shared with those previously reported in the TIB like HLA-DQB1, whereas others like KLHL12 were not reported in TIB. In contrast, the top candidate genes harboring AGVs as previously identified in TIB, like EPAS1 and EGLN1, do not show strong positive selection signals in Deng. Interestingly, Deng also showed a different archaic introgression scenario from that observed in the TIB. Our results suggest that convergent adaptation might be prevalent on the Tibetan Plateau.

Natural Selection Signatures in the Hondo and Ryukyu Japanese Subpopulations

Natural selection signatures across Japanese subpopulations are under-explored. Here we conducted genome-wide selection scans with 622,926 single nucleotide polymorphisms for 20,366 Japanese individuals, who were recruited from the main-islands of Japanese Archipelago (Hondo) and the Ryukyu Archipelago (Ryukyu), representing two major Japanese subpopulations. The integrated haplotype score (iHS) analysis identified several signals in one or both subpopulations. We found a novel candidate locus at IKZF2, especially in Ryukyu. Significant signals were observed in the major histocompatibility complex region in both subpopulations. The lead variants differed and demonstrated substantial allele frequency differences between Hondo and Ryukyu. The lead variant in Hondo tags HLA-A*33:03-C*14:03-B*44:03-DRB1*13:02-DQB1*06:04-DPB1*04:01, a haplotype specific to Japanese and Korean. While in Ryukyu, the lead variant tags DRB1*15:01-DQB1*06:02, which had been recognized as a genetic risk factor for narcolepsy. In contrast, it is reported to confer protective effects against type 1 diabetes and human T lymphotropic virus type 1-associated myelopathy/tropical spastic paraparesis. The FastSMC analysis identified 8 loci potentially affected by selection within the past 20-150 generations, including 2 novel candidate loci. The analysis also showed differences in selection patterns of ALDH2 between Hondo and Ryukyu, a gene recognized to be specifically targeted by selection in East Asian. In summary, our study provided insights into the selection signatures within the Japanese and nominated potential sources of selection pressure.

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.

The putatively high-altitude adaptation of macaque monkeys: Evidence from the fecal metabolome and gut microbiome

Animals living in high-altitude environments, such as the Tibetan Plateau, must face harsh environmental conditions (e.g., hypoxia, cold, and strong UV radiation). These animals' physiological adaptations (e.g., increased red cell production and turnover rate) might also be associated with the gut microbial response. Bilirubin is a component of red blood cell turnover or destruction and is excreted into the intestine and reduced to urobilinoids and/or urobilinogen by gut bacteria. Here, we found that the feces of macaques living in high-altitude regions look significantly browner (with a high concentration of stercobilin, a component from urobilinoids) than those living in low-altitude regions. We also found that gut microbes involved in urobilinogen reduction (e.g., beta-glucuronidase) were enriched in the high-altitude mammal population compared to the low-altitude population. Moreover, the spatial-temporal change in gut microbial function was more profound in the low-altitude macaques than in the high-altitude population, which might be attributed to profound changes in food resources in the low-altitude regions. Therefore, we conclude that a high-altitude environment's stress influences living animals and their symbiotic microbiota.

Genome-wide DNA methylation landscape of four Chinese populations and epigenetic variation linked to Tibetan high-altitude adaptation

DNA methylation (DNAm) is one of the major epigenetic mechanisms in humans and is important in diverse cellular processes. The variation of DNAm in the human population is related to both genetic and environmental factors. However, the DNAm profiles have not been investigated in the Chinese population of diverse ethnicities. Here, we performed double-strand bisulfite sequencing (DSBS) for 32 Chinese individuals representing four major ethnic groups including Han Chinese, Tibetan, Zhuang, and Mongolian. We identified a total of 604,649 SNPs and quantified DNAm at more than 14 million CpGs in the population. We found global DNAm-based epigenetic structure is different from the genetic structure of the population, and ethnic difference only partially explains the variation of DNAm. Surprisingly, non-ethnic-specific DNAm variations showed stronger correlation with the global genetic divergence than these ethnic-specific DNAm. Differentially methylated regions (DMRs) among these ethnic groups were found around genes in diverse biological processes. Especially, these DMR-genes between Tibetan and non-Tibetans were enriched around high-altitude genes including EPAS1 and EGLN1, suggesting DNAm alteration plays an important role in high-altitude adaptation. Our results provide the first batch of epigenetic maps for Chinese populations and the first evidence of the association of epigenetic changes with Tibetans' high-altitude adaptation.

Maternal genetic history of ancient Tibetans over the past 4000 years

The settlement of the Tibetan Plateau epitomizes human adaptation to a high-altitude environment that poses great challenges to human activity. Here, we reconstruct a 4000-year maternal genetic history of Tibetans using 128 ancient mitochondrial genome data from 37 sites in Tibet. The phylogeny of haplotypes M9a1a, M9a1b, D4g2, G2a'c, and D4i show that ancient Tibetans share the most recent common ancestor with ancient Middle and Upper Yellow River populations around the Early and Middle Holocene. In addition, the connections between Tibetans and Northeastern Asians vary over the past 4000 years, with a stronger matrilineal connection between the two during 4000 BP-3000 BP, and a weakened connection after 3000 BP, that are coincident with climate change, followed by a reinforced connection after the Tubo period (1400 BP-1100 BP). Besides, an over 4000-year matrilineal continuity is observed in some of the maternal lineages. We also find the maternal genetic structure of ancient Tibetans is correlated to the geography and interactions between ancient Tibetans and ancient Nepal and Pakistan populations. Overall, the maternal genetic history of Tibetans can be characterized as a long-term matrilineal continuity with frequent internal and external population interactions that are dynamically shaped by geography, climate changes, as well as historical events.

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.

Associations between high-altitude adaptation and risk of cardiovascular diseases: a bidirectional Mendelian randomization study

High-altitude adaptation (HAA) was reported to be significantly associated with reduced risks for multiple cardiovascular diseases (CVDs). However, the causality and direction of the associations are largely uncharacterized. We aimed to examine the potential causal relationships between HAA and six types of CVD, including coronary artery disease (CAD), cerebral aneurysm, ischemic stroke, peripheral artery disease, arrhythmia and atrial fibrillation. We obtained the summary data from largest available genome-wide association study of HAA and six types of CVD. Two-sample bidirectional Mendelian randomization (MR) analyses were performed to infer the causality between them. In the sensitivity analyses, MR-Egger regression analyses and MR-Pleiotropy RESidual Sum and Outlier (MR-PRESSO) global analyses were used to assess the pleiotropic effects; Cochran's Q tests were used to test the heterogeneity by inverse variance-weighted (IVW) and MR-Egger methods; and the leave-one-out analyses were used to examine whether some single nucleotide polymorphisms (SNPs) could influence the results independently. The MR main analyses showed that the genetically instrumented HAA was significantly causally associated with the reduced risks of CAD (odds ratio [OR] = 0.029; 95% confidence interval [CI] = 0.004-0.234; P = 8.64 × 10-4). In contrast, there was no statistically significant relationship between CVDs and HAA. Our findings provide evidence for the causal effects of HAA on the reduced risks of CAD. However, there is no causality of CVDs on HAA. These findings might be helpful in developing the prevention and intervention strategies for CAD.

Massively parallel sequencing of 74 microhaplotypes and forensic characteristics in three Chinese Sino-Tibetan populations

Microhaplotype (MH), as an emerging type of forensic genetic marker in recent years, has the potential to support multiple forensic applications, especially for mixture deconvolution and biogeographic ancestry inference. Herein, we investigated the genotype data of 74 MHs included in a novel MH panel, the Ion AmpliSeq MH-74 Plex Microhaplotype Research Panel, in three Chinese Sino-Tibetan populations (Han, Tibetan, and Yi) using the Ion Torrent semiconductor sequencing. The sequencing performance, allele frequencies, effective number of alleles (Ae), informativeness (In), and forensic parameters were subsequently estimated and calculated. In addition, principal component analysis (PCA) and structure analysis were performed to explore the population relationships among the three populations and the ancestry component distribution. Overall, this novel MH panel is robust and reliable, and has an excellent sequencing performance. The Ae values ranged from 1.0126 to 7.0855 across all samples, and 75.68 % of MHs had Ae values >2.0000. Allele frequencies at some loci varied considerably among the three studied populations, and the mean In value was 0.0195. Moreover, the genetic affinity between Tibetans and Yis was closer than that between Tibetans and Hans. The aforementioned results suggest that the Ion AmpliSeq MH-74 Plex Microhaplotype Research Panel is highly polymorphic in three investigated populations and could be used as an effective tool for human forensics. Although these 74 MHs have demonstrated the competency in continental population stratification, a higher resolution for distinguishing intracontinental subpopulations and a more comprehensive database with sufficient reference population data still remain to be accomplished.

Potential role of the gut microbiota of bumblebee Bombus pyrosoma in adaptation to high-altitude habitats

The gut microbiota affects the health and overall fitness of bumblebees. It can enhance the host's ecological range by leveraging their metabolic capacities. However, the diversity of the gut microbiota and adaptive functional evolution in high-altitude regions remain unclear. To explore how the gut microbiota helps the host adapt to high-altitude environments, we analyzed the differences in diversity and function of the gut microbiota between high- and low-altitude regions through full-length 16S rRNA sequencing. Our results show that high-altitude regions have a lower abundance of Fructobacillus and Saccharibacter compared to low-altitude regions. Additionally, some individuals in low-altitude regions were invaded by opportunistic pathogens. The gut microbiota in high-altitude regions has a greater number of pathways involved in "Protein digestion and absorption" and "Biosynthesis of amino acids," while fewer carbohydrate pathways are involved in "digestion and absorption" and "Salmonella infection." Our finding suggests that plateau hosts typically reduce energy metabolism and enhance immunity in response to adverse environments. Correspondingly, the gut microbiota also makes changes, such as reducing carbohydrate degradation and increasing protein utilization in response to the host. Additionally, the gut microbiota regulates their abundance and function to help the host adapt to adverse high-altitude environments.

Searching for molecular hypoxia sensors among oxygen-dependent enzymes

The ability to sense and respond to changes in cellular oxygen levels is critical for aerobic organisms and requires a molecular oxygen sensor. The prototypical sensor is the oxygen-dependent enzyme PHD: hypoxia inhibits its ability to hydroxylate the transcription factor HIF, causing HIF to accumulate and trigger the classic HIF-dependent hypoxia response. A small handful of other oxygen sensors are known, all of which are oxygen-dependent enzymes. However, hundreds of oxygen-dependent enzymes exist among aerobic organisms, raising the possibility that additional sensors remain to be discovered. This review summarizes known and potential hypoxia sensors among human O2-dependent enzymes and highlights their possible roles in hypoxia-related adaptation and diseases.

High-Altitude Andean H194R HIF2A Allele Is a Hypomorphic Allele

For over 10,000 years, Andeans have resided at high altitude where the partial pressure of oxygen challenges human survival. Recent studies have provided evidence for positive selection acting in Andeans on the HIF2A (also known as EPAS1) locus, which encodes for a central transcription factor of the hypoxia-inducible factor pathway. However, the precise mechanism by which this allele might lead to altitude-adaptive phenotypes, if any, is unknown. By analyzing whole genome sequencing data from 46 high-coverage Peruvian Andean genomes, we confirm evidence for positive selection acting on HIF2A and a unique pattern of variation surrounding the Andean-specific single nucleotide variant (SNV), rs570553380, which encodes for an H194R amino acid substitution in HIF-2α. Genotyping the Andean-associated SNV rs570553380 in a group of 299 Peruvian Andeans from Cerro de Pasco, Peru (4,338 m), reveals a positive association with increased fraction of exhaled nitric oxide, a marker of nitric oxide biosynthesis. In vitro assays show that the H194R mutation impairs binding of HIF-2α to its heterodimeric partner, aryl hydrocarbon receptor nuclear translocator. A knockin mouse model bearing the H194R mutation in the Hif2a gene displays decreased levels of hypoxia-induced pulmonary Endothelin-1 transcripts and protection against hypoxia-induced pulmonary hypertension. We conclude the Andean H194R HIF2A allele is a hypomorphic (partial loss of function) allele.

Gut microbiome signatures of extreme environment adaption in Tibetan pig

Tibetan pigs (TPs) can adapt to the extreme environments in the Tibetan plateau implicated by their self-genome signals, but little is known about roles of the gut microbiota in the host adaption. Here, we reconstructed 8210 metagenome-assembled genomes from TPs (n = 65) living in high-altitude and low-altitude captive pigs (87 from China-CPs and 200 from Europe-EPs) that were clustered into 1050 species-level genome bins (SGBs) at the threshold of 95% average nucleotide identity. 73.47% of SGBs represented new species. The gut microbial community structure analysis based on 1,048 SGBs showed that TPs was significantly different from low-altitude captive pigs. TP-associated SGBs enabled to digest multiple complex polysaccharides, including cellulose, hemicellulose, chitin and pectin. Especially, we found TPs showed the most common enrichment of phyla Fibrobacterota and Elusimicrobia, which were involved in the productions of short- and medium-chain fatty acids (acetic acid, butanoate and propanoate; octanomic, decanoic and dodecanoic acids), as well as in the biosynthesis of lactate, 20 essential amino acids, multiple B vitamins (B1, B2, B3, B5, B7 and B9) and cofactors. Unexpectedly, Fibrobacterota solely showed powerful metabolic capacity, including the synthesis of acetic acid, alanine, histidine, arginine, tryptophan, serine, threonine, valine, B2, B5, B9, heme and tetrahydrofolate. These metabolites might contribute to host adaptation to high-altitude, such as energy harvesting and resistance against hypoxia and ultraviolet radiation. This study provides insights into understanding the role of gut microbiome played in mammalian high-altitude adaptation and discovers some potential microbes as probiotics for improving animal health.

Dynamic Profile of the Yak Mammary Transcriptome during the Lactation Cycle

The objective of this study was to assess the transcriptome of the mammary tissue of four yaks during the whole lactation cycle. For this purpose, biopsies of the mammary gland were performed at -30, -15, 1, 15, 30, 60, 120, 180, and 240 days relative to parturition (d). The transcriptome analysis was performed using a commercial bovine microarray platform and the results were analyzed using several bioinformatic tools. The statistical analysis using an overall false discovery rate ≤ 0.05 for the effect of whole lactation and p < 0.05 for each comparison identified >6000 differentially expressed genes (DEGs) throughout lactation, with a large number of DEGs observed at the onset (1 d vs. -15 d) and at the end of lactation (240 d vs. 180 d). Bioinformatics analysis revealed a major role of genes associated with BTA3, BTA4, BTA6, BTA9, BTA14, and BTA28 in lactation. Functional analysis of DEG underlined an overall induction of lipid metabolism, suggesting an increase in triglycerides synthesis, likely regulated by PPAR signaling. The same analysis revealed an induction of amino acid metabolism and secretion of protein, with a concomitant decrease in proteasome, indicating a major role of amino acid handling and reduced protein degradation in the synthesis and secretion of milk proteins. Glycan biosynthesis was induced for both N-glycan and O-glycan, suggesting increased glycan content in the milk. The cell cycle and immune response, especially antigen processing and presentation, were strongly inhibited during lactation, suggesting that morphological changes are minimized during lactation, while the mammary gland prevents immune hyper-response. Transcripts associated with response to radiation and low oxygen were enriched in the down-regulated DEG affected by the stage of lactation. Except for this last finding, the functions affected by the transcriptomic adaptation to lactation in mammary tissue of yak are very similar to those observed in dairy cows.

Paleoproteomic evidence reveals dairying supported prehistoric occupation of the highland Tibetan Plateau

The extreme environments of the Tibetan Plateau offer considerable challenges to human survival, demanding novel adaptations. While the role of biological and agricultural adaptations in enabling early human colonization of the plateau has been widely discussed, the contribution of pastoralism is less well understood, especially the dairy pastoralism that has historically been central to Tibetan diets. Here, we analyze ancient proteins from the dental calculus (n = 40) of all human individuals with sufficient calculus preservation from the interior plateau. Our paleoproteomic results demonstrate that dairy pastoralism began on the highland plateau by ~3500 years ago. Patterns of milk protein recovery point to the importance of dairy for individuals who lived in agriculturally poor regions above 3700 m above sea level. Our study suggests that dairy was a critical cultural adaptation that supported expansion of early pastoralists into the region's vast, non-arable highlands, opening the Tibetan Plateau up to widespread, permanent human occupation.

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.

Chromosome-level genome assembly and population genomics of Mongolian racerunner (Eremias argus) provide insights into high-altitude adaptation in lizards

BACKGROUND

Although the extreme environmental adaptation of organisms is a hot topic in evolutionary biology, genetic adaptation to high-altitude environment remains poorly characterized in ectothermic animals. Squamates are among the most diverse terrestrial vertebrates, with tremendous ecological plasticity and karyotype diversity, and are a unique model system to investigate the genetic footprints of adaptation.

RESULTS

We report the first chromosome-level assembly of the Mongolian racerunner (Eremias argus) and our comparative genomics analyses found that multiple chromosome fissions/fusions events are unique to lizards. We further sequenced the genomes of 61 Mongolian racerunner individuals that were collected from altitudes ranging from ~ 80 to ~ 2600 m above sea level (m.a.s.l.). Population genomic analyses revealed many novel genomic regions under strong selective sweeps in populations endemic to high altitudes. Genes embedded in those genomic regions are mainly associated with energy metabolism and DNA damage repair pathways. Moreover, we identified and validated two substitutions of PHF14 that may enhance the lizards' tolerance to hypoxia at high altitudes.

CONCLUSIONS

Our study reveals the molecular mechanism of high-altitude adaptation in ectothermic animal using lizard as a research subject and provides a high-quality lizard genomic resource for future research.

Genomic evidence for adaptation to tuberculosis in the Andes before European contact

Most studies focusing on human high-altitude adaptation in the Andean highlands have thus far been focused on Peruvian populations. We present high-coverage whole genomes from Indigenous people living in the Ecuadorian highlands and perform multi-method scans to detect positive natural selection. We identified regions of the genome that show signals of strong selection to both cardiovascular and hypoxia pathways, which are distinct from those uncovered in Peruvian populations. However, the strongest signals of selection were related to regions of the genome that are involved in immune function related to tuberculosis. Given our estimated timing of this selection event, the Indigenous people of Ecuador may have adapted to Mycobacterium tuberculosis thousands of years before the arrival of Europeans. Furthermore, we detect a population collapse that coincides with the arrival of Europeans, which is more severe than other regions of the Andes, suggesting differing effects of contact across high-altitude populations.