Identification of a Tibetan-specific mutation in the hypoxic gene EGLN1 and its contribution to high-altitude adaptation

Skin colour: A window into human phenotypic evolution and environmental adaptation

As modern humans ventured out of Africa and dispersed around the world, they faced novel environmental challenges that led to geographic adaptations including skin colour. Over the long history of human evolution, skin colour has changed dramatically, showing tremendous diversity across different geographical regions, for example, the majority of individuals from the expansive lands of Africa have darker skin, whereas the majority of people from Eurasia exhibit lighter skin. What adaptations did lighter skin confer upon modern humans as they migrated from Africa to Eurasia? What genetic mechanisms underlie the diversity of skin colour observed in different populations? In recent years, scientists have gradually gained a deeper understanding of the interactions between pigmentation gene and skin colour through population-based genomic studies of different groups around the world, particularly in East Asia and Africa. In this review, we summarize our current understanding of 26 skin colour-related pigmentation genes and 48 SNPs that influence skin colour. Important pigmentation genes across three major populations are described in detail: MFSD12, SLC24A5, PDPK1 and DDB1/CYB561A3/TMEM138 influence skin colour in African populations; OCA2, KITLG, SLC24A2, GNPAT and PAH are key to the evolution of skin pigmentation in East Asian populations; and SLC24A5, SLC45A2, TYR, TYRP1, ASIP, MC1R and IRF4 significantly contribute to the lightening of skin colour in European populations. We summarized recent findings in genomic studies of skin colour in populations that implicate diverse geographic environments, local adaptation among populations, gene flow and multi-gene interactions as factors influencing skin colour diversity.

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.

The cardiovascular system at high altitude: A bibliometric and visualization analysis

BACKGROUND

When exposed to high-altitude environments, the cardiovascular system undergoes various changes, the performance and mechanisms of which remain controversial.

AIM

To summarize the latest research advancements and hot research points in the cardiovascular system at high altitude by conducting a bibliometric and visualization analysis.

METHODS

The literature was systematically retrieved and filtered using the Web of Science Core Collection of Science Citation Index Expanded. A visualization analysis of the identified publications was conducted employing CiteSpace and VOSviewer.

RESULTS

A total of 1674 publications were included in the study, with an observed annual increase in the number of publications spanning from 1990 to 2022. The United States of America emerged as the predominant contributor, while Universidad Peruana Cayetano Heredia stood out as the institution with the highest publication output. Notably, Jean-Paul Richalet demonstrated the highest productivity among researchers focusing on the cardiovascular system at high altitude. Furthermore, Peter Bärtsch emerged as the author with the highest number of cited articles. Keyword analysis identified hypoxia, exercise, acclimatization, acute and chronic mountain sickness, pulmonary hypertension, metabolism, and echocardiography as the primary research hot research points and emerging directions in the study of the cardiovascular system at high altitude.

CONCLUSION

Over the past 32 years, research on the cardiovascular system in high-altitude regions has been steadily increasing. Future research in this field may focus on areas such as hypoxia adaptation, metabolism, and cardiopulmonary exercise. Strengthening interdisciplinary and multi-team collaborations will facilitate further exploration of the pathophysiological mechanisms underlying cardiovascular changes in high-altitude environments and provide a theoretical basis for standardized disease diagnosis and treatment.

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.

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.

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.

Hypoxia modeling techniques: A review

Hypoxia is the main cause and effect of a large number of diseases, including the most recent one facing the world, the coronavirus disease (COVID-19). Hypoxia is divided into short-term, long-term, and periodic, it can be the result of diseases, climate change, or living and traveling in the high mountain regions of the world. Since each type of hypoxia can be a cause and a consequence of various physiological changes, the methods for modeling these hypoxias are also different. There are many techniques for modeling hypoxia under experimental conditions. The most common animal for modeling hypoxia is a rat. Hypoxia models (hypoxia simulations) in rats are a tool to study the effect of various conditions on the oxygen supply of the body. These models can provide a necessary information to understand hypoxia and also provide effective treatment, highlighting the importance of various reactions of the body to hypoxia. The main parameters when choosing a model should be reproducibility and the goal that the scientist wants to achieve. Hypoxia in rats can be reproduced both ways exogenously and endogenously. The reason for writing this review was the aim to systematize the models of rats available in the literature in order to facilitate their selection by scientists. The relative strengths and limitations of each model need to be identified and understood in order to evaluate the information obtained from these models and extrapolate these results to humans to develop the necessary generalizations. Despite these problems, animal models have been and remain vital to understanding the mechanisms involved in the development and progression of hypoxia. The eligibility criteria for the selected studies was a comprehensive review of the methods and results obtained from the studies. This made it possible to make generalizations and give recommendations on the application of these methods. The review will assist scientists in choosing an appropriate hypoxia simulation method, as well as assist in interpreting the results obtained with these methods.

Cultural and demic co-diffusion of Tubo Empire on Tibetan Plateau

A high point of Tibetan Plateau (TP) civilization, the expansive Tubo Empire (618-842 AD) wielded great influence across ancient western China. However, whether the Tubo expansion was cultural or demic remains unclear due to sparse ancient DNA sampling. Here, we reported ten ancient genomes at 0.017- to 0.867-fold coverages from the Dulan site with typical Tubo archaeological culture dating to 1308-1130 BP. Nine individuals from three different grave types have close relationship with previously reported ancient highlanders from the southwestern Himalayas and modern core-Tibetan populations. A Dulan-related Tubo ancestry contributed overwhelmingly (95%-100%) to the formation of modern Tibetans. A genetic outlier with dominant Eurasian steppe-related ancestry suggesting a potential population movement into the Tubo-controlled regions from Central Asia. Together with archeological evidence from burial styles and customs, our study suggested the impact of the Tubo empire on the northeast edge of the TP involved both cultural and demic diffusion.

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

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

The HMOX2 polymorphism contributes to the carotid body chemoreflex in European sea-level residents by regulating hypoxic ventilatory responses

This study investigates whether a functional single nucleotide polymorphism of HMOX2 (heme oxygenase-2) (rs4786504 T>C) is involved in individual chemosensitivity to acute hypoxia, as assessed by ventilatory responses, in European individuals. These responses were obtained at rest and during submaximal exercise, using a standardized and validated protocol for exposure to acute normobaric hypoxia. Carriers of the ancestral T allele (n = 44) have significantly lower resting and exercise hypoxic ventilatory responses than C/C homozygous carriers (n = 40). In the literature, a hypoxic ventilatory response threshold to exercise has been identified as an independent predictor of severe high altitude-illness (SHAI). Our study shows that carriers of the T allele have a higher risk of SHAI than carriers of the mutated C/C genotype. Secondarily, we were also interested in COMT (rs4680 G > A) polymorphism, which may be indirectly involved in the chemoreflex response through modulation of autonomic nervous system activity. Significant differences are present between COMT genotypes for oxygen saturation and ventilatory responses to hypoxia at rest. In conclusion, this study adds information on genetic factors involved in individual vulnerability to acute hypoxia and supports the critical role of the ≪ O2 sensor ≫ - heme oxygenase-2 - in the chemosensitivity of carotid bodies in Humans.

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.

Differential methylation in EGLN1 associates with blood oxygen saturation and plasma protein levels in high-altitude pulmonary edema

BACKGROUND

High-altitude (HA, 2500 m) hypoxic exposure evokes a multitude of physiological processes. The hypoxia-sensing genes though influence transcriptional output in disease susceptibility; the exact regulatory mechanisms remain undetermined in high-altitude pulmonary edema (HAPE). Here, we investigated the differential DNA methylation distribution in the two genes encoding the oxygen-sensing HIF-prolyl hydroxylases, prolyl hydroxylase domain protein 2 (PHD2) and factor inhibiting HIF-1α and the consequent contributions to the HAPE pathophysiology.

METHODS

Deep sequencing of the sodium bisulfite converted DNA segments of the two genes, Egl nine homolog 1 (EGLN1) and Hypoxia Inducible Factor 1 Subunit Alpha Inhibitor (HIF1AN), was conducted to analyze the differential methylation distribution in three study groups, namely HAPE-patients (HAPE-p), HAPE-free sojourners (HAPE-f) and healthy HA natives (HLs). HAPE-p and HAPE-f were permanent residents of low altitude (< 200 m) of North India who traveled to Leh (3500 m), India, and were recruited through Sonam Norboo Memorial (SNM) hospital, Leh. HLs were permanent residents of altitudes at and above 3500 m. In addition to the high resolution, bisulfite converted DNA sequencing, gene expression of EGLN1 and HIF1AN and their plasma protein levels were estimated.

RESULTS

A significantly lower methylation distribution of CpG sites was observed in EGLN1 and higher in HIF1AN (P < 0.01) in HAPE-p compared to the two control groups, HAPE-f and HLs. Of note, differential methylation distribution of a few CpG sites, 231,556,748, 231,556,804, 231,556,881, 231,557,317 and 231,557,329, in EGLN1 were significantly associated with the risk of HAPE (OR = 4.79-10.29; P = 0.048-004). Overall, the methylation percentage in EGLN1 correlated with upregulated plasma PHD2 levels (R = - 0.36, P = 0.002) and decreased peripheral blood oxygen saturation (SpO₂) levels (R = 0.34, P = 0.004). We also identified a few regulatory SNPs in the DNA methylation region of EGLN1 covering chr1:231,556,683-231,558,443 suggestive of the functional role of differential methylation distribution of these CpG sites in the regulation of the genes and consequently in the HIF-1α signaling.

CONCLUSIONS

Significantly lower methylation distribution in EGLN1 and the consequent physiological influences annotated its functional epigenetic relevance in the HAPE pathophysiology.

Genomic insight into the population history and biological adaptations of high-altitude Tibetan highlanders in Nagqu

Tibetan, one of the largest indigenous populations living in the high-altitude region of the Tibetan Plateau (TP), has developed a suite of physiological adaptation strategies to cope with the extreme highland environment in TP. Here, we reported genome-wide SNP data from 48 Kham-speaking Nagqu Tibetans and analyzed it with published data from 1,067 individuals in 167 modern and ancient populations to characterize the detailed Tibetan subgroup history and population substructure. Overall, the patterns of allele sharing and haplotype sharing suggested (1) the relatively genetic homogeny between the studied Nagqu Tibetans and ancient Nepalese as well as present-day core Tibetans from Lhasa, Nagqu, and Shigatse; and (2) the close relationship between our studied Kham-speaking Nagqu Tibetans and Kham-speaking Chamdo Tibetans. The fitted qpAdm models showed that the studied Nagqu Tibetans could be fitted as having the main ancestry from late Neolithic upper Yellow River millet farmers and deeply diverged lineages from Southern East Asians (represented by Upper Paleolithic Guangxi_Longlin and Laos_Hoabinhian), and a non-neglectable western Steppe herder-related ancestry (∼3%). We further scanned the candidate genomic regions of natural selection for our newly generated Nagqu Tibetans and the published core Tibetans via FST, iHS, and XP-EHH tests. The genes overlapping with these regions were associated with essential human biological functions such as immune response, enzyme activity, signal transduction, skin development, and energy metabolism. Together, our results shed light on the admixture and evolutionary history of Nagqu Tibetan populations.

A Study on the Molecular Mechanism of High Altitude Heart Disease in Children

Objective

High altitude heart disease (HAHD) is a common pediatric disease in high altitude areas. It usually occurs in people who have lived for a long time or have lived for more than 2500m above sea level. Its common inducement is respiratory tract infection. The clinical differential diagnosis is difficult because the symptoms of HAHD are similar to those of congenital heart disease; Due to the limitation of medical conditions, many patients are in the state of losing follow-up or not seeking medical treatment, resulting in poor prognosis of HAHD and becoming a high-altitude disease with high mortality. Clarifying the molecular mechanism of HAHD, developing early molecular screening technology and accurate treatment methods of HAHD are the key to improve the ability of prevention and treatment of HAHD.

Methods

First, the literature in the PubMed and CNKI databases were screened based on keywords and abstracts. Then, the literature for the study was identified based on the fitness between the content of the literature, the research objectives, and the timeliness of the literature. Finally, a systematic molecular mechanism of HAHD was established by investigating the literature and sorting out the genetic adaptations of Tibetan populations compared with low-altitude populations that migrated to the plateau.

Results

With the investigation of the 48 papers screened, it was found that genes capable of enhancing the hypoxic ventilatory response and resistance to pulmonary hypertension were all correlated with the hypoxia-inducible factor (HIF) pathway, consisting mainly of three pathways, HIF-1α, HIF-2α, and NO.

Conclusion

The low prevalence of HAHD in Tibetan aboriginal children was mainly due to the genetic adaptation of the Tibetan population to the high altitude environment, which coordinated the cellular response to hypoxia by regulating the downstream hypoxia control genes in the HIF pathway.

Ancient genomes from the Himalayas illuminate the genetic history of Tibetans and their Tibeto-Burman speaking neighbors

Present-day Tibetans have adapted both genetically and culturally to the high altitude environment of the Tibetan Plateau, but fundamental questions about their origins remain unanswered. Recent archaeological and genetic research suggests the presence of an early population on the Plateau within the past 40 thousand years, followed by the arrival of subsequent groups within the past 10 thousand years. Here, we obtain new genome-wide data for 33 ancient individuals from high elevation sites on the southern fringe of the Tibetan Plateau in Nepal, who we show are most closely related to present-day Tibetans. They derive most of their ancestry from groups related to Late Neolithic populations at the northeastern edge of the Tibetan Plateau but also harbor a minor genetic component from a distinct and deep Paleolithic Eurasian ancestry. In contrast to their Tibetan neighbors, present-day non-Tibetan Tibeto-Burman speakers living at mid-elevations along the southern and eastern margins of the Plateau form a genetic cline that reflects a distinct genetic history. Finally, a comparison between ancient and present-day highlanders confirms ongoing positive selection of high altitude adaptive alleles.

Comparative Transcriptome Analyses of Gayal (Bos frontalis), Yak (Bos grunniens), and Cattle (Bos taurus) Reveal the High-Altitude Adaptation

Gayal and yak are well adapted to their local high-altitude environments, yet the transcriptional regulation difference of the plateau environment among them remains obscure. Herein, cross-tissue and cross-species comparative transcriptome analyses were performed for the six hypoxia-sensitive tissues from gayal, yak, and cattle. Gene expression profiles for all single-copy orthologous genes showed tissue-specific expression patterns. By differential expression analysis, we identified 3,020 and 1,995 differentially expressed genes (DEGs) in at least one tissue of gayal vs. cattle and yak vs. cattle, respectively. Notably, we found that the adaptability of the gayal to the alpine canyon environment is highly similar to the yak living in the Qinghai-Tibet Plateau, such as promoting red blood cell development, angiogenesis, reducing blood coagulation, immune system activation, and energy metabolism shifts from fatty acid β-oxidation to glycolysis. By further analyzing the common and unique DEGs in the six tissues, we also found that numerous expressed regulatory genes related to these functions are unique in the gayal and yak, which may play important roles in adapting to the corresponding high-altitude environment. Combined with WGCNA analysis, we found that UQCRC1 and COX5A are the shared differentially expressed hub genes related to the energy supply of myocardial contraction in the heart-related modules of gayal and yak, and CAPS is a shared differential hub gene among the hub genes of the lung-related module, which is related to pulmonary artery smooth muscle contraction. Additionally, EDN3 is the unique differentially expressed hub gene related to the tracheal epithelium and pulmonary vasoconstriction in the lung of gayal. CHRM2 is a unique differentially expressed hub gene that was identified in the heart of yak, which has an important role in the autonomous regulation of the heart. These results provide a basis for further understanding the complex transcriptome expression pattern and the regulatory mechanism of high-altitude domestication of gayal and yak.

Adaptation of mammals to hypoxia

Oxygen plays a pivotal role in the metabolism and activities of mammals. However, oxygen is restricted in some environments-subterranean burrow systems or habitats at high altitude or deep in the ocean-and this could exert hypoxic stresses such as oxidative damage on organisms living in these environments. In order to cope with these stresses, organisms have evolved specific strategies to adapt to hypoxia, including changes in physiology, gene expression regulation, and genetic mutations. Here, we review how mammals have adapted to the three high-altitude plateaus of the world, the limited oxygen dissolved in deep water habitats, and underground tunnels, with the aim of better understanding the adaptation of mammals to hypoxia.

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

BACKGROUND

Yaks that inhabit the Tibetan Plateau exhibit striking phenotypic and physiological differences from cattle and have adapted well to the extreme conditions on the plateau. However, the mechanisms used by these animals for the regulation of gene expression at high altitude are not fully understood.

RESULTS

Here, we sequenced nine lung transcriptomes of yaks at altitudes of 3400, 4200 and 5000 m, and low-altitude Zaosheng cattle, which is a closely related species, served as controls. The analysis identified 21,764 mRNAs, 1377 circRNAs and 1209 miRNAs. By comparing yaks and cattle, 4975 mRNAs, 252 circRNAs and 75 miRNAs were identified differentially expressed. By comparing yaks at different altitudes, we identified 756 mRNAs, 64 circRNAs and 83 miRNAs that were differentially expressed (fold change ≥2 and P-value < 0.05). The pathways enriched in the mRNAs, circRNAs and miRNAs identified from the comparison of yaks and cattle were mainly associated with metabolism, including 'glycosaminoglycan degradation', 'pentose and glucuronate interconversions' and 'flavone and flavonol biosynthesis', and the mRNAs, circRNAs and miRNAs identified from the comparison of yaks at different altitude gradients were significantly enriched in metabolic pathways and immune and genetic information processing pathways. The core RNAs were identified from the mRNA-miRNA-circRNA networks constructed using the predominant differentially expressed RNAs. The core genes specific to the difference between yaks and cattle were associated with the endoplasmic reticulum and fat deposition, but those identified from the comparison among yaks at different altitude gradients were associated with maintenance of the normal biological functions of cells.

CONCLUSIONS

This study enhances our understanding of the molecular mechanisms involved in hypoxic adaptation in yaks and might contribute to improvements in the understanding and prevention of hypoxia-related diseases.

Phenotypic and genomic adaptations to the extremely high elevation in plateau zokor (Myospalax baileyi)

The evolutionary outcomes of high elevation adaptation have been extensively described. However, whether widely distributed high elevation endemic animals adopt uniform mechanisms during adaptation to different elevational environments remains unknown, especially with respect to extreme high elevation environments. To explore this, we analysed the phenotypic and genomic data of seven populations of plateau zokor (Myospalax baileyi) along elevations ranging from 2,700 to 4,300 m. Based on whole-genome sequencing data and demographic reconstruction of the evolutionary history, we show that two populations of plateau zokor living at elevations exceeding 3,700 m diverged from other populations nearly 10,000 years ago. Further, phenotypic comparisons reveal stress-dependent adaptation, as two populations living at elevations exceeding 3,700 m have elevated ratios of heart mass to body mass relative to other populations, and the highest population (4,300 m) displays alterations in erythrocytes. Correspondingly, genomic analysis of selective sweeps indicates that positive selection might contribute to the observed phenotypic alterations in these two extremely high elevation populations, with the adaptive cardiovascular phenotypes of both populations possibly evolving under the functional constrains of their common ancestral population. Taken together, phenotypic and genomic evidence demonstrates that heterogeneous stressors impact adaptations to extreme elevations and reveals stress-dependent and genetically constrained adaptation to hypoxia, collectively providing new insights into the high elevation adaptation.

Genomes reveal selective sweeps in kiang and donkey for high-altitude adaptation

Over the last several hundred years, donkeys have adapted to high-altitude conditions on the Tibetan Plateau. Interestingly, the kiang, a closely related equid species, also inhabits this region. Previous reports have demonstrated the importance of specific genes and adaptive introgression in divergent lineages for adaptation to hypoxic conditions on the Tibetan Plateau. Here, we assessed whether donkeys and kiangs adapted to the Tibetan Plateau via the same or different biological pathways and whether adaptive introgression has occurred. We assembled a de novo genome from a kiang individual and analyzed the genomes of five kiangs and 93 donkeys (including 24 from the Tibetan Plateau). Our analyses suggested the existence of a strong hard selective sweep at the EPAS1 locus in kiangs. In Tibetan donkeys, however, another gene, i.e., EGLN1, was likely involved in their adaptation to high altitude. In addition, admixture analysis found no evidence for interspecific gene flow between kiangs and Tibetan donkeys. Our findings indicate that despite the short evolutionary time scale since the arrival of donkeys on the Tibetan Plateau, as well as the existence of a closely related species already adapted to hypoxia, Tibetan donkeys did not acquire adaptation via admixture but instead evolved adaptations via a different biological pathway.

Adaptive introgression of the beta-globin cluster in two Andean waterfowl

Introgression of beneficial alleles has emerged as an important avenue for genetic adaptation in both plant and animal populations. In vertebrates, adaptation to hypoxic high-altitude environments involves the coordination of multiple molecular and cellular mechanisms, including selection on the hypoxia-inducible factor (HIF) pathway and the blood-O2 transport protein hemoglobin (Hb). In two Andean duck species, a striking DNA sequence similarity reflecting identity by descent is present across the ~20 kb β-globin cluster including both embryonic (HBE) and adult (HBB) paralogs, though it was yet untested whether this is due to independent parallel evolution or adaptive introgression. In this study, we find that identical amino acid substitutions in the β-globin cluster that increase Hb-O2 affinity have likely resulted from historical interbreeding between high-altitude populations of two different distantly-related species. We examined the direction of introgression and discovered that the species with a deeper mtDNA divergence that colonized high altitude earlier in history (Anas flavirostris) transferred adaptive genetic variation to the species with a shallower divergence (A. georgica) that likely colonized high altitude more recently possibly following a range shift into a novel environment. As a consequence, the species that received these β-globin variants through hybridization might have adapted to hypoxic conditions in the high-altitude environment more quickly through acquiring beneficial alleles from the standing, hybrid-origin variation, leading to faster evolution.

Gain-of-function Tibetan PHD2D4E;C127S variant suppresses monocyte function: A lesson in inflammatory response to inspired hypoxia

Background

We have previously described an evolutionarily selected Tibetan prolyl hydroxylase-2 (PHD2^D4E;C127S) variant that degrades the hypoxia-inducible factor (HIFα) more efficiently and protects these highlanders from hypoxia-triggered elevation in haemoglobin concentration. High altitude is known to cause acute mountain sickness (AMS) and high-altitude pulmonary edema (HAPE) in a section of rapidly ascending non-acclimatised lowlanders. These morbidities are often accompanied by inflammatory response and exposure to hypobaric hypoxia is presumed to be the principal causative agent. We have investigated whether PHD2^D4E;C127S variant is associated with prevention of hypoxia-mediated inflammatory milieu in Tibetan highlanders and therefore identify a potential target to regulate inflammation.

Methods

We genotyped the Tibetans using DNA isolated from whole blood. Thereafter immunophenotying was performed on PBMCs from homozygous PHD2^D4E;C127S and PHD2^WT individuals using flow cytometry. RNA isolated from these individuals was used to evaluate the peripheral level of important transcripts associated with immune as well as hypoxia response employing the nCounter technology. The ex-vivo findings were validated by generating monocytic cell lines (U937 cell line) expressing PHD2^D4E;C127S and PHD2^WT variants post depletion of endogenous PHD2. We had also collected whole blood samples from healthy travellers and travellers afflicted with AMS and HAPE to evaluate the significance of our ex-vivo and in vitro findings. Hereafter, we also attempted to resolve hypoxia-triggered inflammation in vitro as well as in vivo by augmenting the function of PHD2 using alpha-ketoglutarate (αKG), a co-factor of PHD2.

Findings

We report that homozygous PHD2^D4E;C127S highlanders harbour less inflammatory and patrolling monocytes in circulation as compared to Tibetan PHD2^WT highlanders. In response to in vitro hypoxia, secretion of IL6 and IL1β from PHD2^D4E;C127S monocytes, and their chemotactic response compared to the PHD2^WT are compromised, corresponding to the down-modulated expression of related signalling molecules RELA, JUN, STAT1, ATF2 and CXCR4. We verified these functional outcomes in monocytic U937 cell line engineered to express PHD2^D4E;C127S and confirmed the down-modulation of the signalling molecules at protein level under hypoxia. In contrast, non-Tibetan sojourners with AMS and HAPE at high altitude (3,600 m above sea level) displayed significant increase in these inflammatory parameters. Our data henceforth underline the role of gain-of-function of PHD2 as the rate limiting factor to harness hyper-activation of monocytes in hypoxic environment. Therefore upon pre-treatment with αKG, we observed diminished inflammatory response of monocytes in vitro and reduction in leukocyte infiltration to the lungs in mice exposed to normobaric hypoxia.

Interpretation

Our report suggests that gain-of-function PHD2 ^D4E;C127S variant can therefore protect against inflammation elicited by hypobaric hypoxia. Augmentation of PHD2 activity therefore may be an important method to alleviate inflammatory response to inspired hypoxia.

Funding

This study is supported by the Department of Biotechnology, Government of India.

Characterization of structural variation in Tibetans reveals new evidence of high-altitude adaptation and introgression

BACKGROUND

Structural variation (SV) acts as an essential mutational force shaping the evolution and function of the human genome. However, few studies have examined the role of SVs in high-altitude adaptation and little is known of adaptive introgressed SVs in Tibetans so far.

RESULTS

Here, we generate a comprehensive catalog of SVs in a Chinese Tibetan (n = 15) and Han (n = 10) population using nanopore sequencing technology. Among a total of 38,216 unique SVs in the catalog, 27% are sequence-resolved for the first time. We systematically assess the distribution of these SVs across repeat sequences and functional genomic regions. Through genotyping in additional 276 genomes, we identify 69 Tibetan-Han stratified SVs and 80 candidate adaptive genes. We also discover a few adaptive introgressed SV candidates and provide evidence for a deletion of 335 base pairs at 1p36.32.

CONCLUSIONS

Overall, our results highlight the important role of SVs in the evolutionary processes of Tibetans' adaptation to the Qinghai-Tibet Plateau and provide a valuable resource for future high-altitude adaptation studies.

Effects of circadian variation, lifestyle and environment on hematological parameters: A narrative review

The complete blood count (CBC) is the most widely prescribed laboratory test. It plays a key role in screening, diagnosing, and monitoring a variety of medical disorders. Preanalytical and analytical variables are responsible for more than 50% of laboratory errors that may lead to spurious CBC results. The effects of blood sampling, transport, storage, and analytical errors on hematological parameters have been well described. Circadian variation and changes in lifestyle and environment can also affect blood cells. It has been extensively studied in the past, but highly variable methodology and the presence of confounding factors have provided scattered and inconsistent results. We have investigated the literature to define the impact of circadian variation, modification of the sleep-wake cycle, acute and chronic exercise, eating habits, alcohol, tobacco, drugs of abuse, high-altitude, heat/cold exposure, and air pollution on CBC results. The affected cell type along with the intensity and duration of changes are detailed for each condition. We aim at providing a comprehensive overview of which situations may induce clinically significant changes and have to be taken into account by healthcare professionals before considering a hematological parameter as pathological and requesting complementary tests.

KLF4, a Key Regulator of a Transitive Triplet, Acts on the TGF-β Signaling Pathway and Contributes to High-Altitude Adaptation of Tibetan Pigs

Tibetan pigs are native mammalian species on the Tibetan Plateau that have evolved distinct physiological traits that allow them to tolerate high-altitude hypoxic environments. However, the genetic mechanism underlying this adaptation remains elusive. Here, based on multitissue transcriptional data from high-altitude Tibetan pigs and low-altitude Rongchang pigs, we performed a weighted correlation network analysis (WGCNA) and identified key modules related to these tissues. Complex network analysis and bioinformatics analysis were integrated to identify key genes and three-node network motifs. We found that among the six tissues (muscle, liver, heart, spleen, kidneys, and lungs), lung tissue may be the key organs for Tibetan pigs to adapt to hypoxic environment. In the lung tissue of Tibetan pigs, we identified KLF4, BCL6B, EGR1, EPAS1, SMAD6, SMAD7, KDR, ATOH8, and CCN1 genes as potential regulators of hypoxia adaption. We found that KLF4 and EGR1 genes might simultaneously regulate the BCL6B gene, forming a KLF4-EGR1-BCL6B complex. This complex, dominated by KLF4, may enhance the hypoxia tolerance of Tibetan pigs by mediating the TGF-β signaling pathway. The complex may also affect the PI3K-Akt signaling pathway, which plays an important role in angiogenesis caused by hypoxia. Therefore, we postulate that the KLF4-EGR1-BCL6B complex may be beneficial for Tibetan pigs to survive better in the hypoxia environments. Although further molecular experiments and independent large-scale studies are needed to verify our findings, these findings may provide new details of the regulatory architecture of hypoxia-adaptive genes and are valuable for understanding the genetic mechanism of hypoxic adaptation in mammals.

High-Altitude Adaptation: Mechanistic Insights from Integrated Genomics and Physiology

Population genomic analyses of high-altitude humans and other vertebrates have identified numerous candidate genes for hypoxia adaptation, and the physiological pathways implicated by such analyses suggest testable hypotheses about underlying mechanisms. Studies of highland natives that integrate genomic data with experimental measures of physiological performance capacities and subordinate traits are revealing associations between genotypes (e.g., hypoxia-inducible factor gene variants) and hypoxia-responsive phenotypes. The subsequent search for causal mechanisms is complicated by the fact that observed genotypic associations with hypoxia-induced phenotypes may reflect second-order consequences of selection-mediated changes in other (unmeasured) traits that are coupled with the focal trait via feedback regulation. Manipulative experiments to decipher circuits of feedback control and patterns of phenotypic integration can help identify causal relationships that underlie observed genotype–phenotype associations. Such experiments are critical for correct inferences about phenotypic targets of selection and mechanisms of adaptation.

Physiological Genomics of Adaptation to High-Altitude Hypoxia

Population genomic studies of humans and other animals at high altitude have generated many hypotheses about the genes and pathways that may have contributed to hypoxia adaptation. Future advances require experimental tests of such hypotheses to identify causal mechanisms. Studies to date illustrate the challenge of moving from lists of candidate genes to the identification of phenotypic targets of selection, as it can be difficult to determine whether observed genotype-phenotype associations reflect causal effects or secondary consequences of changes in other traits that are linked via homeostatic regulation. Recent work on high-altitude models such as deer mice has revealed both plastic and evolved changes in respiratory, cardiovascular, and metabolic traits that contribute to aerobic performance capacity in hypoxia, and analyses of tissue-specific transcriptomes have identified changes in regulatory networks that mediate adaptive changes in physiological phenotype. Here we synthesize recent results and discuss lessons learned from studies of high-altitude adaptation that lie at the intersection of genomics and physiology.

Roles of HIF and 2-Oxoglutarate-Dependent Dioxygenases in Controlling Gene Expression in Hypoxia

Simple Summary

Hypoxia—reduction in oxygen availability—plays key roles in both physiological and pathological processes. Given the importance of oxygen for cell and organism viability, mechanisms to sense and respond to hypoxia are in place. A variety of enzymes utilise molecular oxygen, but of particular importance to oxygen sensing are the 2-oxoglutarate (2-OG) dependent dioxygenases (2-OGDs). Of these, Prolyl-hydroxylases have long been recognised to control the levels and function of Hypoxia Inducible Factor (HIF), a master transcriptional regulator in hypoxia, via their hydroxylase activity. However, recent studies are revealing that such dioxygenases are involved in almost all aspects of gene regulation, including chromatin organisation, transcription and translation.

Abstract

Hypoxia—reduction in oxygen availability—plays key roles in both physiological and pathological processes. Given the importance of oxygen for cell and organism viability, mechanisms to sense and respond to hypoxia are in place. A variety of enzymes utilise molecular oxygen, but of particular importance to oxygen sensing are the 2-oxoglutarate (2-OG) dependent dioxygenases (2-OGDs). Of these, Prolyl-hydroxylases have long been recognised to control the levels and function of Hypoxia Inducible Factor (HIF), a master transcriptional regulator in hypoxia, via their hydroxylase activity. However, recent studies are revealing that dioxygenases are involved in almost all aspects of gene regulation, including chromatin organisation, transcription and translation. We highlight the relevance of HIF and 2-OGDs in the control of gene expression in response to hypoxia and their relevance to human biology and health.

Effect of EGLN1 Genetic Polymorphisms on Hemoglobin Concentration in Andean Highlanders

The physiological characteristics of Andean natives living at high altitudes have been investigated extensively, with many studies reporting that Andean highlanders have a higher hemoglobin (Hb) concentration than other highlander populations. It has previously been reported that positive natural selection has acted independently on the egl-9 family hypoxia inducible factor 1 (EGLN1) gene in Tibetan and Andean highlanders and is related to Hb concentration in Tibetans. However, no study has yet revealed the genetic determinants of Hb concentration in Andeans even though several single-nucleotide polymorphisms (SNPs) in EGLN1 have previously been examined. Therefore, we explored the relationship between hematological measurements and tag SNPs designed to cover the whole EGLN1 genomic region in Andean highlanders living in Bolivia. Our findings indicated that haplotype frequencies estimated from the EGLN1 SNPs were significantly correlated with Hb concentration in the Bolivian highlanders. Moreover, we found that an Andean-dominant haplotype related to high Hb level may have expanded rapidly in ancestral Andean highlander populations. Analysis of genotype data in an ~436.3 kb genomic region containing EGLN1 using public databases indicated that the population structure based on EGLN1 genetic markers in Andean highlanders was largely different from that in other human populations. This finding may be related to an intrinsic or adaptive physiological characteristic of Andean highlanders. In conclusion, the high Hb concentrations in Andean highlanders can be partly characterized by EGLN1 genetic variants.

Hypoxia-inducible factors not only regulate but also are myeloid-cell treatment targets

Hypoxia describes limited oxygen availability at the cellular level. Myeloid cells are exposed to hypoxia at various bodily sites and even contribute to hypoxia by consuming large amounts of oxygen during respiratory burst. Hypoxia-inducible factors (HIFs) are ubiquitously expressed heterodimeric transcription factors, composed of an oxygen-dependent α and a constitutive β subunit. The stability of HIF-1α and HIF-2α is regulated by oxygen-sensing prolyl-hydroxylases (PHD). HIF-1α and HIF-2α modify the innate immune response and are context dependent. We provide a historic perspective of HIF discovery, discuss the molecular components of the HIF pathway, and how HIF-dependent mechanisms modify myeloid cell functions. HIFs enable myeloid-cell adaptation to hypoxia by up-regulating anaerobic glycolysis. In addition to effects on metabolism, HIFs control chemotaxis, phagocytosis, degranulation, oxidative burst, and apoptosis. HIF-1α enables efficient infection defense by myeloid cells. HIF-2α delays inflammation resolution and decreases antitumor effects by promoting tumor-associated myeloid-cell hibernation. PHDs not only control HIF degradation, but also regulate the crosstalk between innate and adaptive immune cells thereby suppressing autoimmunity. HIF-modifying pharmacologic compounds are entering clinical practice. Current indications include renal anemia and certain cancers. Beneficial and adverse effects on myeloid cells should be considered and could possibly lead to drug repurposing for inflammatory disorders.

Cytoprotection by a naturally occurring variant of ATP5G1 in Arctic ground squirrel neural progenitor cells

Many organisms in nature have evolved mechanisms to tolerate severe hypoxia or ischemia, including the hibernation-capable Arctic ground squirrel (AGS). Although hypoxic or ischemia tolerance in AGS involves physiological adaptations, little is known about the critical cellular mechanisms underlying intrinsic AGS cell resilience to metabolic stress. Through cell survival-based cDNA expression screens in neural progenitor cells, we identify a genetic variant of AGS Atp5g1 that confers cell resilience to metabolic stress. Atp5g1 encodes a subunit of the mitochondrial ATP synthase. Ectopic expression in mouse cells and CRISPR/Cas9 base editing of endogenous AGS loci revealed causal roles of one AGS-specific amino acid substitution in mediating cytoprotection by AGS ATP5G1. AGS ATP5G1 promotes metabolic stress resilience by modulating mitochondrial morphological change and metabolic functions. Our results identify a naturally occurring variant of ATP5G1 from a mammalian hibernator that critically contributes to intrinsic cytoprotection against metabolic stress.

When animals hibernate, they lower their body temperature and metabolism to conserve the energy they need to withstand cold harsh winters. One such animal is the Arctic ground squirrel, an extreme hibernator that can drop its body temperatures to below 0°C. This hibernation ability means the cells of Arctic ground squirrels can survive severe shortages of blood and oxygen. But, it is unclear how their cells are able to endure this metabolic stress.

To answer this question, Singhal, Bai et al. studied the cells of Arctic ground squirrels for unique features that might make them more durable to stress. Examining the genetic code of these resilient cells revealed that Arctic ground squirrels may have a variant form of a protein called ATP5G1. This protein is found in a cellular compartment called the mitochondria, which is responsible for supplying energy to the rest of the cell and therefore plays an important role in metabolic processes.

Singhal, Bai et al. found that when this variant form of ATP5G1 was introduced into the cells of mice, their mitochondria was better at coping with stress conditions, such as low oxygen, low temperature and poisoning. Using a gene editing tool to selectively substitute some of the building blocks, also known as amino acids, that make up the ATP5G1 protein revealed that improvements to the mitochondria were caused by switching specific amino acids. However, swapping these amino acids, which presumably affects the role of ATP5G1, did not completely remove the cells’ resilience to stress. This suggests that variants of other genes and proteins may also be involved in providing protection.

These findings provide the first evidence of a protein variant that is responsible for protecting cells during the metabolic stress conditions caused by hibernation. The approach taken by Singhal, Bai et al. could be used to identify and study other proteins that increase resilience to metabolic stress. These findings could help develop new treatments for diseases caused by a limited blood supply to human organs, such as a stroke or heart attack.

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Hypoxia-induced pulmonary hypertension-Utilizing experiments of nature

An increase in pulmonary artery pressure is a common observation in adult mammals exposed to global alveolar hypoxia. It is considered a maladaptive response that places an increased workload on the right ventricle. The mechanisms initiating and maintaining the elevated pressure are of considerable interest in understanding pulmonary vascular homeostasis. There is an expectation that identifying the key molecules in the integrated vascular response to hypoxia will inform potential drug targets. One strategy is to take advantage of experiments of nature, specifically, to understand the genetic basis for the inter-individual variation in the pulmonary vascular response to acute and chronic hypoxia. To date, detailed phenotyping of highlanders has focused on haematocrit and oxygen saturation rather than cardiovascular phenotypes. This review explores what we can learn from those studies with respect to the pulmonary circulation. LINKED ARTICLES: This article is part of a themed issue on Risk factors, comorbidities, and comedications in cardioprotection. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v178.1/issuetoc.

Sexual reproduction and genetic polymorphism within the cosmopolitan marine diatom Pseudo-nitzschia pungens

Different clades belonging to the cosmopolitan marine diatom Pseudo-nitzschia pungens appear to be present in different oceanic environments, however, a ‘hybrid zone’, where populations of different clades interbreed, has also been reported. Many studies have investigated the sexual reproduction of P. pungens, focused on morphology and life cycle, rather than the role of sexual reproduction in mixing the genomes of their parents. We carried out crossing experiments to determine the sexual compatibility/incompatibility between different clades of P. pungens, and examined the genetic polymorphism in the ITS2 region. Sexual reproduction did not occur only between clades II and III under any of experimental temperature conditions. Four offspring strains were established between clade I and III successfully. Strains established from offspring were found interbreed with other offspring strains as well as viable with their parental strains. We confirmed the hybrid sequence patterns between clades I and III and found novel sequence types including polymorphic single nucleotide polymorphisms (SNPs) in the offspring strains. Our results implicate that gene exchange and mixing between different clades are still possible, and that sexual reproduction is a significant ecological strategy to maintain the genetic diversity within this diatom species.

EPAS1 regulates proliferation of erythroblasts in chronic mountain sickness

Excessive erythrocytosis (EE) is a characteristic of chronic mountain sickness (CMS). Currently, the pathogenesis of CMS remains unclear. This study was intended to investigate the role of EPAS1 in the proliferation of erythroblasts in CMS. Changes of HIF-1α and EPAS1/HIF-2α in the bone marrow erythroblasts of 21 patients with CMS and 14 control subjects residing at the same altitudes were determined by RT-qPCR and western blotting. We also developed a lentiviral vector, Lv-EPAS1/sh-EPAS1, to over-express/silence EPAS1 in K562 cells. Cells cycle and proliferation were detected by flow cytometry. Transcriptome analyses were carried out on Illumina. CMS patients showed a higher expression of EPAS1/HIF-2α in the bone marrow erythroblasts than those of controls. Variations in EPAS1 expression in CMS patients were positively correlated with RBC levels, and negatively correlated with SaO₂. Over-expressing of EPAS1 in K562 cells accelerated the erythroid cells cycle progression and promoted the erythroid cells proliferation-and vice versa. Transcriptome data indicated that proliferation-related DEGs were significantly enriched in EPAS1 overexpression/silencing K562 cells. Our results suggest that EPAS1 might participate in the pathogenesis of EE by regulating the proliferation of erythroblasts.

Tibetan PHD2, an allele with loss-of-function properties

Tibetans have adapted to the chronic hypoxia of high altitude and display a distinctive suite of physiologic adaptations, including augmented hypoxic ventilatory response and resistance to pulmonary hypertension. Genome-wide studies have consistently identified compelling genetic signatures of natural selection in two genes of the Hypoxia Inducible Factor pathway, PHD2 and HIF2A The product of the former induces the degradation of the product of the latter. Key issues regarding Tibetan PHD2 are whether it is a gain-of-function or loss-of-function allele, and how it might contribute to high-altitude adaptation. Tibetan PHD2 possesses two amino acid changes, D4E and C127S. We previously showed that in vitro, Tibetan PHD2 is defective in its interaction with p23, a cochaperone of the HSP90 pathway, and we proposed that Tibetan PHD2 is a loss-of-function allele. Here, we report that additional PHD2 mutations at or near Asp-4 or Cys-127 impair interaction with p23 in vitro. We find that mice with the Tibetan Phd2 allele display augmented hypoxic ventilatory response, supporting this loss-of-function proposal. This is phenocopied by mice with a mutation in p23 that abrogates the PHD2:p23 interaction. Hif2a haploinsufficiency, but not the Tibetan Phd2 allele, ameliorates hypoxia-induced increases in right ventricular systolic pressure. The Tibetan Phd2 allele is not associated with hemoglobin levels in mice. We propose that Tibetans possess genetic alterations that both activate and inhibit selective outputs of the HIF pathway to facilitate successful adaptation to the chronic hypoxia of high altitude.

Neural network correlates of high-altitude adaptive genetic variants in Tibetans: A pilot, exploratory study

Although substantial progress has been made in the identification of genetic substrates underlying physiology, neuropsychology, and brain organization, the genotype–phenotype associations remain largely unknown in the context of high‐altitude (HA) adaptation. Here, we related HA adaptive genetic variants in three gene loci (EGLN1, EPAS1, and PPARA) to interindividual variance in a set of physiological characteristics, neuropsychological tests, and topological attributes of large‐scale structural and functional brain networks in 135 indigenous Tibetan highlanders. Analyses of individual HA adaptive single‐nucleotide polymorphisms (SNPs) revealed that specific SNPs selectively modulated physiological characteristics (erythrocyte level, ratio between forced expiratory volume in the first second to forced vital capacity, arterial oxygen saturation, and heart rate) and structural network centrality (the left anterior orbital gyrus) with no effects on neuropsychology or functional brain networks. Further analyses of genetic adaptive scores, which summarized the overall degree of genetic adaptation to HA, revealed significant correlations only with structural brain networks with respect to local interconnectivity of the whole networks, intermodule communication between the right frontal and parietal module and the left occipital module, nodal centrality in several frontal regions, and connectivity strength of a subnetwork predominantly involving in intramodule edges in the right temporal and occipital module. Moreover, the associations were dependent on gene loci, weight types, or topological scales. Together, these findings shed new light on genotype–phenotype interactions under HA hypoxia and have important implications for developing new strategies to optimize organism and tissue responses to chronic hypoxia induced by extreme environments or diseases.

Population Structure, and Selection Signatures Underlying High-Altitude Adaptation Inferred From Genome-Wide Copy Number Variations in Chinese Indigenous Cattle

Copy number variations (CNVs) have been demonstrated as crucial substrates for evolution, adaptation and breed formation. Chinese indigenous cattle breeds exhibit a broad geographical distribution and diverse environmental adaptability. Here, we analyzed the population structure and adaptation to high altitude of Chinese indigenous cattle based on genome-wide CNVs derived from the high-density BovineHD SNP array. We successfully detected the genome-wide CNVs of 318 individuals from 24 Chinese indigenous cattle breeds and 37 yaks as outgroups. A total of 5,818 autosomal CNV regions (683 bp-4,477,860 bp in size), covering ~14.34% of the bovine genome (UMD3.1), were identified, showing abundant CNV resources. Neighbor-joining clustering, principal component analysis (PCA), and population admixture analysis based on these CNVs support that most Chinese cattle breeds are hybrids of Bos taurus taurus (hereinafter to be referred as Bos taurus) and Bos taurus indicus (Bos indicus). The distribution patterns of the CNVs could to some extent be related to the geographical backgrounds of the habitat of the breeds, and admixture among cattle breeds from different districts. We analyzed the selective signatures of CNVs positively involved in high-altitude adaptation using pairwise Fst analysis within breeds with a strong Bos taurus background (taurine-type breeds) and within Bos taurus×Bos indicus hybrids, respectively. CNV-overlapping genes with strong selection signatures (at top 0.5% of Fst value), including LETM1 (Fst = 0.490), TXNRD2 (Fst = 0.440), and STUB1 (Fst = 0.420) within taurine-type breeds, and NOXA1 (Fst = 0.233), RUVBL1 (Fst = 0.222), and SLC4A3 (Fst=0.154) within hybrids, were potentially involved in the adaptation to hypoxia. Thus, we provide a new profile of population structure from the CNV aspects of Chinese indigenous cattle and new insights into high-altitude adaptation in cattle.

Association of EGLN1 gene with high aerobic capacity of Peruvian Quechua at high altitude

Highland native Andeans have resided at altitude for millennia. They display high aerobic capacity (VO₂max) at altitude, which may be a reflection of genetic adaptation to hypoxia. Previous genomewide (GW) scans for natural selection have nominated Egl-9 homolog 1 gene (EGLN1) as a candidate gene. The encoded protein, EGLN1/PHD2, is an O₂ sensor that controls levels of the Hypoxia Inducible Factor-α (HIF-α), which regulates the cellular response to hypoxia. From GW association and analysis of covariance performed on a total sample of 429 Peruvian Quechua and 94 US lowland referents, we identified 5 EGLN1 SNPs associated with higher VO₂max (L⋅min^-1 and mL⋅min^-1⋅kg^-1) in hypoxia (rs1769793, rs2064766, rs2437150, rs2491403, rs479200). For 4 of these SNPs, Quechua had the highest frequency of the advantageous (high VO₂max) allele compared with 25 diverse lowland comparison populations from the 1000 Genomes Project. Genotype effects were substantial, with high versus low VO₂max genotype categories differing by ∼11% (e.g., for rs1769793 SNP genotype TT = 34.2 mL⋅min^-1⋅kg^-1 vs. CC = 30.5 mL⋅min^-1⋅kg^-1). To guard against spurious association, we controlled for population stratification. Findings were replicated for EGLN1 SNP rs1769793 in an independent Andean sample collected in 2002. These findings contextualize previous reports of natural selection at EGLN1 in Andeans, and support the hypothesis that natural selection has increased the frequency of an EGLN1 causal variant that enhances O₂ delivery or use during exercise at altitude in Peruvian Quechua.

Whole mitochondrial genome analysis of highland Tibetan ethnicity using massively parallel sequencing

Mitochondrial DNA (mtDNA) is a key player in numerous multifaceted and intricate biological processes and plays a pivotal role in dissecting the peopling of different populations, due to its maternally inherited property and comparatively high mutation rate. In this study, 119 Tibetan individuals from the Muli Tibetan Autonomous County of China (average altitude above 3,000 m) were employed in mitochondrial genome (mitogenome) sequencing by massively parallel sequencing (MPS) techniques using the Precision ID mtDNA Whole Genome Panel on an Ion S5XL system. The dataset presented 88 distinct haplotypes, resulting in the haplotype diversity of 0.9909. The majority of haplotypes were assigned to East Asian lineages and the distribution of haplogroups of Muli Tibetan significantly differed from reference Tibetan populations. The maximum parsimony phylogeny reconstructed by 119 newly generated mitogenomes revealed 12 major Muli Tibetan lineages. Intriguingly, a Sherpa-specific sub-haplogroup A15c1 with the lack of mutations at 4216 and 15,924 was discerned in our dataset, which suggested that the maternal gene pool of Sherpas may derive from Tibetan populations. The shared haplogroups between Muli Tibetan and lowland Han Chinese hinted that these lineages may derive from non-Tibetans and have already differentiated before their arrival on the Tibetan Plateau. Furthermore, extensive pairwise population comparisons displayed that Muli Tibetan had a closer genetic relationship with ethnically or linguistically close Nyingtri Tibetan, Nyingtri Lhoba and Chamdo Tibetan populations. Genetic affinity was also observed between the Muli Tibetan and North Han Chinese. Collectively, the results generated in this study enriched the existing forensic mtDNA database and raised additional interest in the application of whole mitogenome sequencing in forensic investigations.

Prioritizing natural-selection signals from the deep-sequencing genomic data suggests multi-variant adaptation in Tibetan highlanders

Human genetic adaptation to high altitudes (>2500 m) has been extensively studied over the last few years, but few functional adaptive genetic variants have been identified, largely owing to the lack of deep-genome sequencing data available to previous studies. Here, we build a list of putative adaptive variants, including 63 missense, 7 loss-of-function, 1,298 evolutionarily conserved variants and 509 expression quantitative traits loci. Notably, the top signal of selection is located in TMEM247, a transmembrane protein-coding gene. The Tibetan version of TMEM247 harbors one high-frequency (76.3%) missense variant, rs116983452 (c.248C > T; p.Ala83Val), with the T allele derived from archaic ancestry and carried by >94% of Tibetans but absent or in low frequencies (<3%) in non-Tibetan populations. The rs116983452-T is strongly and positively correlated with altitude and significantly associated with reduced hemoglobin concentration (p = 5.78 × 10-5), red blood cell count (p = 5.72 × 10-7) and hematocrit (p = 2.57 × 10-6). In particular, TMEM247-rs116983452 shows greater effect size and better predicts the phenotypic outcome than any EPAS1 variants in association with adaptive traits in Tibetans. Modeling the interaction between TMEM247-rs116983452 and EPAS1 variants indicates weak but statistically significant epistatic effects. Our results support that multiple variants may jointly deliver the fitness of the Tibetans on the plateau, where a complex model is needed to elucidate the adaptive evolution mechanism.

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.

PGG.SNV: understanding the evolutionary and medical implications of human single nucleotide variations in diverse populations

Despite the tremendous growth of the DNA sequencing data in the last decade, our understanding of the human genome is still in its infancy. To understand the implications of genetic variants in the light of population genetics and molecular evolution, we developed a database, PGG.SNV ( https://www.pggsnv.org ), which gives much higher weight to previously under-investigated indigenous populations in Asia. PGG.SNV archives 265 million SNVs across 220,147 present-day genomes and 1018 ancient genomes, including 1009 newly sequenced genomes, representing 977 global populations. Moreover, estimation of population genetic diversity and evolutionary parameters is available in PGG.SNV, a unique feature compared with other databases.

Population History and Altitude-Related Adaptation in the Sherpa

The first ascent of Mount Everest by Tenzing Norgay and Sir Edmund Hillary in 1953 brought global attention to the Sherpa people and human performance at altitude. The Sherpa inhabit the Khumbu Valley of Nepal, and are descendants of a population that has resided continuously on the Tibetan plateau for the past ∼25,000 to 40,000 years. The long exposure of the Sherpa to an inhospitable environment has driven genetic selection and produced distinct adaptive phenotypes. This review summarizes the population history of the Sherpa and their physiological and genetic adaptation to hypoxia. Genomic studies have identified robust signals of positive selection across EPAS1, EGLN1, and PPARA, that are associated with hemoglobin levels, which likely protect the Sherpa from altitude sickness. However, the biological underpinnings of other adaptive phenotypes such as birth weight and the increased reproductive success of Sherpa women are unknown. Further studies are required to identify additional signatures of selection and refine existing Sherpa-specific adaptive phenotypes to understand how genetic factors have underpinned adaptation in this population. By correlating known and emerging signals of genetic selection with adaptive phenotypes, we can further reveal hypoxia-related biological mechanisms of adaptation. Ultimately this work could provide valuable information regarding treatments of hypoxia-related illnesses including stroke, heart failure, lung disease and cancer.

Insights into hypoxic adaptation in Tibetan chicken embryos from comparative proteomics

Tibetan chicken, an indigenous breed, is highly adapted to the extreme environment of the Qinghai-Tibet Plateau. It serves as a model organism to identify genetic differences between hypoxia-adapted and lowland breeds. However, the mechanisms underlying hypoxia adaptation are yet unclear. This study aimed to identify differently abundant proteins (DAPs) and elucidate the mechanisms involved in hypoxic adaptation in the Tibetan chicken. In this study, we obtained proteome data for the embryonic heart tissues of Tibetan and Chahua chickens incubated under hypoxia (TCH and CHH) and normoxia (TCN and CHN) using isobaric tags for relative and absolute quantitation (iTRAQ) technology. We identified 4210 proteins from 53,352 unique peptides in the heart tissue of chicken embryos. Pairwise TCH vs. CHH, TCH vs. TCN, CHH vs. CHN, and TCN vs. CHN comparisons revealed 118, 176, 103, and 162 differently abundant proteins, respectively. Several key proteins (EGLN1, MAP2K2, MYLK, QARS, NOTCH2, and MYH7) and pathways (glutathione metabolism, PPAR signaling pathway, and vascular smooth muscle contraction) were identified and considered important candidates for high-altitude adaptation in Tibetan chicken. This study provides novel insights into the chicken embryonic heart tissue and furthers the current understanding of the mechanisms of survival among animals in high-altitude environments.

Functional Assays to Screen and Dissect Genomic Hits: Doubling Down on the National Investment in Genomic Research

The National Institutes of Health have made substantial investments in genomic studies and technologies to identify DNA sequence variants associated with human disease phenotypes. The National Heart, Lung, and Blood Institute has been at the forefront of these commitments to ascertain genetic variation associated with heart, lung, blood, and sleep diseases and related clinical traits. Genome-wide association studies, exome- and genome-sequencing studies, and exome-genotyping studies of the National Heart, Lung, and Blood Institute-funded epidemiological and clinical case-control studies are identifying large numbers of genetic variants associated with heart, lung, blood, and sleep phenotypes. However, investigators face challenges in identification of genomic variants that are functionally disruptive among the myriad of computationally implicated variants. Studies to define mechanisms of genetic disruption encoded by computationally identified genomic variants require reproducible, adaptable, and inexpensive methods to screen candidate variant and gene function. High-throughput strategies will permit a tiered variant discovery and genetic mechanism approach that begins with rapid functional screening of a large number of computationally implicated variants and genes for discovery of those that merit mechanistic investigation. As such, improved variant-to-gene and gene-to-function screens-and adequate support for such studies-are critical to accelerating the translation of genomic findings. In this White Paper, we outline the variety of novel technologies, assays, and model systems that are making such screens faster, cheaper, and more accurate, referencing published work and ongoing work supported by the National Heart, Lung, and Blood Institute's R21/R33 Functional Assays to Screen Genomic Hits program. We discuss priorities that can accelerate the impressive but incomplete progress represented by big data genomic research.

Genetic structure and forensic characteristics of Tibeto-Burman-speaking Ü-Tsang and Kham Tibetan Highlanders revealed by 27 Y-chromosomal STRs

Culturally diverse Tibetans (Ü-Tsang, Kham and Ando) harboring a unique molecular mechanism that allows them to successfully adapt to hypoxic environments in the Qinghai-Tibet Plateau have been a subject of great interest in medical genetics, linguistics, archeology and forensic science. However, forensic characteristics and genetic variations of the Y-chromosomal 27-marker haplotype included in the Yfiler Plus system in the Ü-Tsang and Kham Tibeto-Burman-speaking Tibetans remain unexplored. Thus, we genotyped 27 Y-STRs in 230 Shigatse Ü-Tsang Tibetans (SUT) and 172 Chamdo Kham Tibetans (CKT) to investigate the forensic characterization and genetic affinity of Chinese Tibetan Highlanders. The haplotype diversities were 0.999962028 in SUT and 0.999796002 in CKT. Forensic diversity measures indicated that this 27-Y-STR amplification system is appropriate for routine forensic applications, such as identifying and separating unrelated males in deficiency paternity cases, male disaster victims and missing person identification and determining male components in sexual assault cases. Moreover, the genetic relationships among 63 worldwide populations (16,282 individuals), 16 Asian populations, and 21 Chinese populations were analyzed and reconstructed using principal component analysis, multidimensional scaling plots and a phylogenetic tree. Considerable genetic differences were observed between Tibetan populations and other geographically/ethnically diverse populations (Han Chinese). Our studied SUT and CKT have a genetically closer relationship with Gansu Ando Tibetans than with other Asians. In total, our analyses indicated that subpopulation structures exist among Asian and Chinese populations, and population-specific reference databases should be established for forensic applications.

Thin Air, Thick Vessels: Historical and Current Perspectives on Hypoxic Pulmonary Hypertension

The association between pulmonary hypertension (PH) and hypoxia is well-established, with two key mechanistic processes, hypoxic pulmonary vasoconstriction and hypoxia-induced vascular remodeling, driving changes in pulmonary arterial pressure. In contrast to other forms of pulmonary hypertension, the vascular changes induced by hypoxia are reversible, both in humans returning to sea-level from high altitude and in animal models. This raises the intriguing possibility that the molecular drivers of these hypoxic processes could be targeted to modify pulmonary vascular remodeling in other contexts. In this review, we outline the history of research into PH and hypoxia, before discussing recent advances in our understanding of this relationship at the molecular level, focussing on the role of the oxygen-sensing transcription factors, hypoxia inducible factors (HIFs). Emerging links between HIF and vascular remodeling highlight the potential utility in inhibiting this pathway in pulmonary hypertension and raise possible risks of activating this pathway using HIF-stabilizing medications.

UNVEILing connections between genotype, phenotype, and fitness in natural populations

Understanding the links between genetic variation and fitness in natural populations is a central goal of evolutionary genetics. This monumental task spans the fields of classical and molecular genetics, population genetics, biochemistry, physiology, developmental biology, and ecology. Advances to our molecular and developmental toolkits are facilitating integrative approaches across these traditionally separate fields, providing a more complete picture of the genotype-phenotype map in natural and non-model systems. Here, we summarize research presented at the first annual symposium of the UNVEIL Network, an NSF-funded collaboration between the University of Montana and the University of Nebraska, Lincoln, which took place from the 1st to the 3rd of June, 2018. We discuss how this body of work advances basic evolutionary science, what it implies for our ability to predict evolutionary change, and how it might inform novel conservation strategies.

Natural selection on TMPRSS6 associated with the blunted erythropoiesis and improved blood viscosity in Tibetan pigs

Tibetan pigs, indigenous to Tibetan plateau, are well adapted to hypoxia. So far, there have been not any definitively described genes and functional sites responsible for hypoxia adaptation for the Tibetan pig. The whole genome-wide association studies in human suggested that genetic variations in TMPRSS6 was associated with hemoglobin concentration (HGB) and red cell counts (RBC). Here we conducted resequencing of the nearly entire genomic region (40.1 kb) of the candidate gene TMPRSS6 in 40 domestic pigs and 40 wild boars along continuous altitudes and identified 708 SNPs, in addition to an indel (CGTG/----) in the intron 10. We conduct the CGTG indel in 838 domestic pigs, both the CGTG deletion frequency and the pairwise r² linkage disequilibrium showed an increase with elevated altitudes, suggesting that TMPRSS6 has been under Darwinian positive selection. As the conserved core sequence of hypoxia-response elements (HREs), the deletion of CGTG in Tibetan pigs decreased the expression levels of TMPRSS6 mRNA and protein in the liver revealed by real-time quantitative PCR and western blot, respectively. We compared domestic pigs and Tibetan pigs living continuous altitudes, found that the blood-related traits with the increase of altitude, however, the HGB did not increase with the elevation in Tibetan pigs. Genotype association analysis results dissected a genetic effect on reducing HGB by 13.25 g/L in Gongbo'gyamda Tibetan pigs, decreasing mean corpuscular volume (MCV) by 4.79 fl in Diqing Tibetan pigs. In conclusion, the CGTG deletion of TMPRSS6 resulted in lower HGB and smaller MCV, which could reflect a blunting erythropoiesis and improving blood viscosity as well as erythrocyte deformability. It remains to be determined whether a blunting of erythropoiesis for TMPRSS6 or others genetic effects are the physiological adaptations among Tibetan pigs.