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

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.

Gene expression plasticity followed by genetic change during colonization in a high-elevation environment

Phenotypic plasticity facilitates organismal invasion of novel environments, and the resultant phenotypic change may later be modified by genetic change, so called 'plasticity first.' Herein, we quantify gene expression plasticity and regulatory adaptation in a wild bird (Eurasian Tree Sparrow) from its original lowland (ancestral stage), experimentally implemented hypoxia acclimation (plastic stage), and colonized highland (colonized stage). Using a group of co-expressed genes from the cardiac and flight muscles, respectively, we demonstrate that gene expression plasticity to hypoxia tolerance is more often reversed than reinforced at the colonized stage. By correlating gene expression change with muscle phenotypes, we show that colonized tree sparrows reduce maladaptive plasticity that largely associated with decreased hypoxia tolerance. Conversely, adaptive plasticity that is congruent with increased hypoxia tolerance is often reinforced in the colonized tree sparrows. Genes displaying large levels of reinforcement or reversion plasticity (i.e. 200% of original level) show greater genetic divergence between ancestral and colonized populations. Overall, our work demonstrates that gene expression plasticity at the initial stage of high-elevation colonization can be reversed or reinforced through selection-driven adaptive modification.

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.

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.

Microfluidic intestinal organoid-on-a-chip uncovers therapeutic targets by recapitulating oxygen dynamics of intestinal IR injury

Increasing evidence demonstrates that mammals have different reactions to hypoxia with varied oxygen dynamic patterns. It takes ∼24 h for tri-gas incubator to achieve steady cell hypoxia, which fails to recapitulate ultrafast oxygen dynamics of intestinal ischemia/reperfusion (IR) injury. Inspired from the structure of native intestinal villi, we engineered an intestinal organoid chip embedded with engineered artificial microvessels based on co-axial microfluidic technology by using pH-responsive ZIF-8/sodium alginate scaffold. The chip was featured on: (i) eight times the oxygen exchange efficiency compared with the conventional device, tri-gas incubator, (ii) implantation of intestinal organoid reproducing all types of intestinal epithelial cells, and (iii) bio-responsiveness to hypoxia and reoxygenation (HR) by presenting metabolism disorder, inflammatory reaction, and cell apoptosis. Strikingly, it was found for the first time that Olfactomedin 4 (Olfm4) was the most significantly down-regulated gene under a rapid HR condition by sequencing the RNA from the organoids. Mechanistically, OLFM4 played protective functions on HR-induced cell inflammation and tissue damage by inhibiting the NF-kappa B signaling activation, thus it could be used as a therapeutic target. Altogether, this study overcomes the issue of mismatched oxygen dynamics between in vitro and in vivo, and sets an example of next-generation multisystem-interactive organoid chip for finding precise therapeutic targets of IR injury.

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.

Evolutionary origin of genomic structural variations in domestic yaks

Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.

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.

Promising Natural Medicines for the Treatment of High-Altitude Illness

Li Li, Lin Lin, Bo Wen, Peng-cheng Zhao, Da-sheng Liu, Guo-ming Pang, Zi-rong Wang, Yong Tan, and Cheng Lu. Promising natural medicines for the treatment of high-altitude illness. High Alt Med Biol. 24:175-185, 2023.-High-altitude illness (HAI) is a dangerous disease characterized by oxidative stress, inflammatory damage and hemodynamic changes in the body that can lead to severe damage to the lungs, heart, and brain. Natural medicines are widely known for their multiple active ingredients and pharmacological effects, which may be important in the treatment of HAI. In this review, we outline the specific types of HAI and the underlying pathological mechanisms and summarize the currently documented natural medicines applied in the treatment of acute mountain sickness and high-altitude cerebral edema, high-altitude pulmonary edema, chronic mountain sickness, and high-altitude pulmonary hypertension. Their sources, types, and medicinal sites are summarized, and their active ingredients, pharmacological effects, related mechanisms, and potential toxicity are discussed. In conclusion, natural medicines, as an acceptable complementary and alternative strategy with fewer side effects and more long-term application, can provide a reference for developing more natural antialtitude sickness medicines in the future and have good application prospects in HAI treatment.

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.

Divergent contributions of coding and noncoding sequences to initial high-altitude adaptation in passerine birds endemic to the Qinghai-Tibet Plateau

Early events in the evolution of an ancestral lineage can shape the adaptive patterns of descendant species, but the evolutionary mechanisms driving initial adaptation from an ancestor remain largely unexplored. High-altitude adaptations have been extensively explored from the viewpoint of protein-coding genes; however, the contribution of noncoding regions remains relatively neglected. Here, we integrate genomic and transcriptomic data to investigate adaptive evolution in the ancestor of three high-altitude snowfinch species endemic to the Qinghai-Tibet Plateau. Our genome-wide scan for adaptation in the snowfinch ancestor identifies strong adaptation signals in functions of development and metabolism for the coding genes, but in functions of the nervous system development for noncoding regions. This pattern is exclusive to the snowfinch ancestor compared to a control ancestral lineage subject to weak selection. Changes in noncoding regions in the snowfinch ancestor, especially those nearest to coding genes, may be disproportionately associated with the differential expression of genes in the brain tissue compared to other tissues. Extensive gene expression in the brain tissue can be further altered via genetic regulatory networks of transcription factors harbouring potential accelerated regulatory regions (e.g., the development-related transcription factor YEATS4). Altogether, our study provides new evidence concerning how coding and noncoding sequences work through decoupled pathways in initial adaptation to the selective pressure of high-altitude environments. The analysis highlights the idea that noncoding sequences may be promising elements in facilitating the rapid evolution and adaptation to high altitudes.

Varied hypoxia adaptation patterns of embryonic brain at different development stages between Tibetan and Dwarf laying chickens

BACKGROUND

Tibetan chickens (Gallus gallus; TBCs), an indigenous breed distributed in the Qinghai-Tibet Plateau, are well adapted to the hypoxic environment. Currently, the molecular genetic basis of hypoxia adaptation in TBCs remains unclear. This study investigated hypoxia adaptation patterns of embryonic brain at different development stages by integrating analysis of the transcriptome with our previously published metabolome data in TBCs and Dwarf Laying Chickens (DLCs), a lowland chicken breed.

RESULTS

During hypoxia, the results revealed that 1334, 578, and 417 differentially expressed genes (DEGs) (|log₂ fold change|>1, p-value < 0.05) on days 8, 12, and 18 of development, respectively between TBCs and DLCs. Gene Ontology (GO) and pathway analyses revealed that DEGs are mainly related to metabolic pathways, vessel development, and immune response under hypoxia. This is consistent with our metabolome data that TBCs have higher energy metabolism than DLCs during hypoxia. Some vital DEGs between TBCs and DLCs, such as EPAS1, VEGFD, FBP1, FBLN5, LDHA, and IL-6 which are involved in the HIF pathway and hypoxia regulation.

CONCLUSION

These results suggest varied adaptation patterns between TBCs and DLCs under hypoxia. Our study provides a basis for uncovering the molecular regulation mechanism of hypoxia adaptation in TBCs and a potential application of hypoxia adaptation research for other animals living on the Qinghai-Tibet Plateau, and may even contribute to the study of brain diseases caused by hypoxia.

Regulation of erythroid differentiation in K562 cells by the EPAS1-IRS2 axis under hypoxic conditions

Red blood cells (RBCs) produced in vitro have the potential to alleviate the worldwide demand for blood transfusion. Hematopoietic cell differentiation and proliferation are triggered by numerous cellular physiological processes, including low oxygen concentration (<5%). In addition, hypoxia inducible factor 2α (HIF-2α) and insulin receptor substrate 2 (IRS2) were found to be involved in the progression of erythroid differentiation. However, the function of the HIF-2α-IRS2 axis in the progression of erythropoiesis is not yet fully understood. Therefore, we used an in vitro model of erythropoiesis generated from K562 cells transduced with shEPAS1 at 5% O₂ in the presence or absence of the IRS2 inhibitor NT157. We observed that erythroid differentiation was accelerated in K562 cells by hypoxia. Conversely, knockdown of EPAS1 expression reduced IRS2 expression and erythroid differentiation. Intriguingly, inhibition of IRS2 could impair the progression of hypoxia-induced erythropoiesis without affecting EPAS1 expression. These findings indicated that the EPAS1-IRS2 axis may be a crucial pathway that regulates erythropoiesis and that drugs targeting this pathway may become promising agents for promoting erythroid differentiation.

Resurrecting the alternative splicing landscape of archaic hominins using machine learning

Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.

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.

Hypoxia induces alterations in tRNA modifications involved in translational control

BACKGROUND

Adaptation to high-altitude hypobaric hypoxia has been shown to require a set of physiological traits enabled by an associated set of genetic modifications, as well as transcriptome regulation. These lead to both lifetime adaptation of individuals to hypoxia at high altitudes and generational evolution of populations as seen for instance in those of Tibet. Additionally, RNA modifications, which are sensitive to environmental exposure, have been shown to play pivotal biological roles in maintaining the physiological functions of organs. However, the dynamic RNA modification landscape and related molecular mechanisms in mouse tissues under hypobaric hypoxia exposure remain to be fully understood. Here, we explore the tissue-specific distribution pattern of multiple RNA modifications across mouse tissues.

RESULTS

By applying an LC-MS/MS-dependent RNA modification detection platform, we identified the distribution of multiple RNA modifications in total RNA, tRNA-enriched fragments, and 17-50-nt sncRNAs across mouse tissues; these patterns were associated with the expression levels of RNA modification modifiers in different tissues. Moreover, the tissue-specific abundance of RNA modifications was sensitively altered across different RNA groups in a simulated high-altitude (over 5500 m) hypobaric hypoxia mouse model with the activation of the hypoxia response in mouse peripheral blood and multiple tissues. RNase digestion experiments revealed that the alteration of RNA modification abundance under hypoxia exposure impacted the molecular stability of tissue total tRNA-enriched fragments and isolated individual tRNAs, such as tRNA^Ala, tRNA^val, tRNA^Glu, and tRNA^Leu. In vitro transfection experiments showed that the transfection of testis total tRNA-enriched fragments from the hypoxia group into GC-2spd cells attenuated the cell proliferation rate and led to a reduction in overall nascent protein synthesis in cells.

CONCLUSIONS

Our results reveal that the abundance of RNA modifications for different classes of RNAs under physiological conditions is tissue-specific and responds to hypobaric hypoxia exposure in a tissue-specific manner. Mechanistically, the dysregulation of tRNA modifications under hypobaric hypoxia attenuated the cell proliferation rate, facilitated the sensitivity of tRNA to RNases, and led to a reduction in overall nascent protein synthesis, suggesting an active role of tRNA epitranscriptome alteration in the adaptive response to environmental hypoxia exposure.

A pleiotropic hypoxia-sensitive EPAS1 enhancer is disrupted by adaptive alleles in Tibetans

In Tibetans, noncoding alleles in EPAS1-whose protein product hypoxia-inducible factor 2α (HIF-2α) drives the response to hypoxia-carry strong signatures of positive selection; however, their functional mechanism has not been systematically examined. Here, we report that high-altitude alleles disrupt the activity of four EPAS1 enhancers in one or more cell types. We further characterize one enhancer (ENH5) whose activity is both allele specific and hypoxia dependent. Deletion of ENH5 results in down-regulation of EPAS1 and HIF-2α targets in acute hypoxia and in a blunting of the transcriptional response to sustained hypoxia. Deletion of ENH5 in mice results in dysregulation of gene expression across multiple tissues. We propose that pleiotropic adaptive effects of the Tibetan alleles in EPAS1 underlie the strong selective signal at this gene.

Estimates of genomic heritability and genome-wide association studies for blood parameters in Akkaraman sheep

The aim of this study was to estimate genomic heritability and the impact that genetic backgrounds have on blood parameters in Akkaraman sheep by conducting genome-wide association studies and regional heritability mapping analysis. Genomic heritability estimates for blood parameters ranged from 0.00 to 0.55, indicating that measured phenotypes have a low to moderate heritability. A total of 7 genome- and 13 chromosome-wide significant SNPs were associated with phenotypic changes in 15 blood parameters tested. Accordingly, SCN7A, SCN9A, MYADM-like, CCDC67, ITGA9, MGAT5, SLC19A1, AMPH, NTRK2, MSRA, SLC35F3, SIRT6, CREB3L3, and NAV3 genes as well as three undefined regions (LOC101117887, LOC106991526 and LOC105608461) were suggested as candidates. Most of the identified genes were involved in basic biological processes that are essential to immune system function and cellular growth; specific functions include cellular transport, histone deacetylation, cell differentiation, erythropoiesis, and endocytosis. The top significant SNP for HCT, MCH, and MCHC was found within a genomic region mainly populated by the MYADM-like gene family. This region was previously suggested to be under historical selection pressure in many sheep breeds from various parts of the world. These results have implications on animal breeding program studies due to the effect that the genetic background has on blood parameters, which underlying many productive and wellness related traits.

Microevolutionary mechanism of high-altitude adaptation in Tibetan chicken populations from an elevation gradient

As an indigenous breed, the Tibetan chicken is found in highland regions and shows physiological adaptations to high altitude; however, the genetic changes that determine these adaptations remain elusive. We assumed that the microevolution of the Tibetan chicken occurred from lowland to highland regions with a continuous elevation range. In this study, we analyzed the genome of 188 chickens from lowland areas to the high-altitude regions of the Tibetan plateau with four altitudinal levels. Phylogenetic analysis revealed that Tibetan chickens are significantly different from other altitude chicken populations. Reconstruction of the demographic history showed that the migration and admixture events of the Tibetan chicken occurred at different times. The genome of the Tibetan chicken was also used to analyze positive selection pressure that is associated with high-altitude adaptation, revealing the well-known candidate gene that participates in oxygen binding (HBAD), as well as other novel potential genes (e.g., HRG and ANK2) that are related to blood coagulation and cardiovascular efficiency. Our study provides novel insights regarding the evolutionary history and microevolution mechanisms of the high-altitude adaptation in the Tibetan chicken.

Effective prediction of potential ferroptosis critical genes in clinical colorectal cancer

BACKGROUND

Colon cancer is common worldwide, with high morbidity and poor prognosis. Ferroptosis is a novel form of cell death driven by the accumulation of iron-dependent lipid peroxides, which differs from other programmed cell death mechanisms. Programmed cell death is a cancer hallmark, and ferroptosis is known to participate in various cancers, including colon cancer. Novel ferroptosis markers and targeted colon cancer therapies are urgently needed. To this end, we performed a preliminary exploration of ferroptosis-related genes in colon cancer to enable new treatment strategies.

METHODS

Ferroptosis-related genes in colon cancer were obtained by data mining and screening for differentially expressed genes (DEGs) using bioinformatics analysis tools. We normalized the data across four independent datasets and a ferroptosis-specific database. Identified genes were validated by immunohistochemical analysis of pathological and healthy clinical samples.

RESULTS

We identified DEGs in colon cancer that are involved in ferroptosis. Among these, five core genes were found: ELAVL1, GPX2, EPAS1, SLC7A5, and HMGB1. Bioinformatics analyses revealed that the expression of all five genes, except for EPAS1, was higher in tumor tissues than in healthy tissues.

CONCLUSIONS

The preliminary exploration of the five core genes revealed that they are differentially expressed in colon cancer, playing an essential role in ferroptosis. This study provides a foundation for subsequent research on ferroptosis in colon cancer.

Expression and Variations in EPAS1 Associated with Oxygen Metabolism in Sheep

Endothelial PAS domain protein 1 gene (EPAS1) is a member of the HIF gene family. This gene encodes a transcription factor subunit that is involved in the induction of oxygen-regulated genes. Several studies have demonstrated that a mutation in EPAS1 could affect oxygen sensing, polycythemia, and hemoglobin level. However, whether EPAS1 mutation affects sheep oxygen metabolism is still unknown. Therefore, we explored the relationship between the variation of EPAS1 and oxygen metabolism in sheep. In this study, variations in ovine EPAS1 exon 15 were investigated in 332 Tibetan sheep and 339 Hu sheep by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) analysis. In addition, we studied the effect of these variations on blood gas in 176 Tibetan sheep and 231 Hu sheep. Finally, the mRNA expression of EPAS1 in six tissues of Hu sheep and Tibetan sheep living at different altitudes (2500 m, 3500 m, and 4500 m) was analyzed by real-time quantitative PCR (RT-qPCR). Four alleles (A, B, C, and D) were detected, and their distributions highly differed between Tibetan sheep and Hu sheep. In Tibetan sheep, B was the dominant allele, and C and D alleles were rare, whereas all four alleles were common in Hu sheep. Six single nucleotide polymorphisms (SNPs) were identified between the four alleles and one of them was non-synonymous (p.F606L). While studying the blood gas levels in Tibetan sheep and Hu sheep, one variant region was found to be associated with an elevated pO₂ and sO₂, which suggested that variations in EPAS1 are associated with oxygen metabolism in sheep. RT-qPCR results showed that EPAS1 was expressed in the six tissues of Hu sheep and Tibetan sheep at different altitudes. In addition, the expression of EPAS1 in four tissues (heart, liver, spleen, and longissimus dorsi muscle) of Hu sheep was lower than that in Tibetan sheep from three different altitudes, and the expression of EPAS1 was positively correlated with the altitude. These results indicate that the variations and expression of EPAS1 is closely related to oxygen metabolism.

Time Domains of Hypoxia Responses and -Omics Insights

The ability to respond rapidly to changes in oxygen tension is critical for many forms of life. Challenges to oxygen homeostasis, specifically in the contexts of evolutionary biology and biomedicine, provide important insights into mechanisms of hypoxia adaptation and tolerance. Here we synthesize findings across varying time domains of hypoxia in terms of oxygen delivery, ranging from early animal to modern human evolution and examine the potential impacts of environmental and clinical challenges through emerging multi-omics approaches. We discuss how diverse animal species have adapted to hypoxic environments, how humans vary in their responses to hypoxia (i.e., in the context of high-altitude exposure, cardiopulmonary disease, and sleep apnea), and how findings from each of these fields inform the other and lead to promising new directions in basic and clinical hypoxia research.

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.

A highland-adaptation mutation of the Epas1 protein increases its stability and disrupts the circadian clock in the plateau pika

The Qinghai-Tibet Plateau (QTP) harbors hundreds of species well adapted to its extreme conditions, including its low-oxygen (hypoxic) atmosphere. Here, we show that the plateau pika-a keystone mammal of the QTP-lacks robust circadian rhythms. The major form of the plateau pika Epas1 protein includes a 24-residue insert caused by a point mutation at the 5' juncture site of Intron14 and is more stable than other mammalian orthologs. Biochemical studies reveal that an Epas1-Bmal1 complex with lower trans-activation activity occupies the E1/E2 motifs at the promoter of the core-clock gene Per2, thus explaining how an Epas1 mutation-selected in the hypoxic conditions of the QTP-disrupts the molecular clockwork. Importantly, experiments with hypoxic chambers show that mice expressing the plateau pika Epas1 ortholog in their suprachiasmatic nucleus have dysregulated central clocks, and pika Epas1 knockin mice reared in hypoxic conditions exhibit dramatically reduced heart damage compared with wild-type animals.

Physiology and Transcriptomics Analysis Reveal the Contribution of Lungs on High-Altitude Hypoxia Adaptation in Tibetan Sheep

The Tibetan sheep is an indigenous species on the Tibetan plateau with excellent adaptability to high-altitude hypoxia and is distributed at altitudes of 2500–5000 m. The high-altitude hypoxia adaptation of Tibetan sheep requires adaptive reshaping of multiple tissues and organs, especially the lungs. To reveal the mechanisms of adaptation at the tissue and molecular levels in the lungs of Tibetan sheep under hypoxic conditions at different altitudes, we performed light and electron microscopic observations, transcriptomic sequencing, and enzyme-linked immunosorbent assay studies on the lungs of Tibetan sheep from three altitudes (2500, 3500, and 4500 m). The results showed that in addition to continuous increase in pulmonary artery volume, thickness, and elastic fiber content with altitude, Tibetan sheep increase the hemoglobin concentration at an altitude of 3500 m, while they decrease the Hb concentration and increase the surface area of gas exchange and capacity of the blood at an altitude of 4500 m. Other than that, some important differentially expressed genes related to angiogenesis (FNDC1, HPSE, and E2F8), vasomotion and fibrogenesis (GJA4, FAP, COL1A1, COL1A2, COL3A1, and COL14A1), and gas transport (HBB, HBA1, APOLD1, and CHL1) were also identified; these discoveries at the molecular level explain to some extent the physiological findings. In conclusion, the lungs of Tibetan sheep adopt different strategies when adapting to different altitudes, and these findings are valuable for understanding the basis of survival of indigenous species on the Tibetan plateau.

The Influence of Environmental Hypoxia on Hemostasis—A Systematic Review

Humans have been ascending to high altitudes for centuries, with a growing number of professional- and leisure-related sojourns occurring in this millennium. A multitude of scientific reports on hemostatic disorders at high altitude suggest that hypoxia is an independent risk factor. However, no systematic analysis of the influence of environmental hypoxia on coagulation, fibrinolysis and platelet function has been performed. To fill this gap, we performed a systematic literature review, including only the data of healthy persons obtained during altitude exposure (&lt;60 days). The results were stratified by the degree of hypoxia and sub-categorized into active and passive ascents and sojourns. Twenty-one studies including 501 participants were included in the final analysis. Since only one study provided relevant data, no conclusions regarding moderate altitudes (1,500–2,500 m) could be drawn. At high altitude (2,500–5,400 m), only small pathophysiological changes were seen, with a possible impact of increasing exercise loads. Elevated thrombin generation seems to be balanced by decreased platelet activation. Viscoelastic methods do not support increased thrombogenicity, with fibrinolysis being unaffected by high altitude. At extreme altitude (5,400–8,850 m), the limited data showed activation of coagulation in parallel with stimulation of fibrinolysis. Furthermore, multiple confounding variables at altitude, like training status, exercise load, fluid status and mental stress, prevent definitive conclusions being drawn on the impact of hypoxia on hemostasis. Thus, we cannot support the hypothesis that hypoxia triggers hypercoagulability and increases the risk of thromboembolic disorders, at least in healthy sojourners.

A Ferroptosis-Related Gene Prognostic Index Associated With Biochemical Recurrence and Radiation Resistance for Patients With Prostate Cancer Undergoing Radical Radiotherapy

Background: Ferroptosis is a new type of programmed cell death which has been reported to be involved in the development of various cancers. In this study, we attempted to explore the possible links between ferroptosis and prostate cancer (PCa), and a novel ferroptosis-related gene prognostic index (FGPI) was constructed to predict biochemical recurrence (BCR) and radiation resistance for PCa patients undergoing radical radiotherapy (RRT). Moreover, the tumor immune microenvironment (TME) of PCa was analyzed.

Methods: We merged four GEO datasets by removing batch effects. All analyses were conducted with R version 3.6.3 and its suitable packages. Cytoscape 3.8.2 was used to establish a network of transcriptional factor and competing endogenous RNA.

Results: We established the FGPI based on ACSL3 and EPAS1. We observed that FGPI was an independent risk factor of BCR for PCa patients (HR: 3.03; 95% CI: 1.68–5.48), consistent with the result of internal validation (HR: 3.44; 95% CI: 1.68–7.05). Furthermore, FGPI showed high ability to identify radiation resistance (AUC: 0.963; 95% CI: 0.882–1.00). LncRNA PART1 was significantly associated with BCR and might modulate the mRNA expression of EPAS1 and ACSL3 through interactions with 60 miRNAs. Gene set enrichment analysis indicated that FGPI was enriched in epithelial–mesenchymal transition, allograft rejection, TGF beta signaling pathway, and ECM receptor interaction. Immune checkpoint and m6A analyses showed that PD-L2, CD96, and METTL14 were differentially expressed between BCR and no BCR groups, among which CD96 was significantly associated with BCR-free survival (HR: 1.79; 95% CI: 1.06–3.03). We observed that cancer-related fibroblasts (CAFs), macrophages, stromal score, immune score, estimate score, and tumor purity were differentially expressed between BCR and no BCR groups and closely related to BCR-free survival (HRs were 2.17, 1.79, 2.20, 1.93, 1.92, and 0.52 for cancer-related fibroblasts, macrophages, stromal score, immune score, estimate score, and tumor purity, respectively). Moreover, cancer-related fibroblasts (coefficient: 0.20), stromal score (coefficient: 0.14), immune score (coefficient: 0.14), estimate score (coefficient: 0.15), and tumor purity (coefficient: −0.15) were significantly related to FGPI, among which higher positive correlation between cancer-related fibroblasts and FGPI was observed.

Conclusion: We found that FGPI based on ACSL3 and EPAS1 might be used to predict BCR and radiation resistance for PCa patients. CD96 and PD-L2 might be a possible target for drug action. Besides, we highlighted the importance of immune evasion in the process of BCR.

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.

High Altitude Pregnancies and Vascular Dysfunction: Observations From Latin American Studies

An estimated human population of 170 million inhabit at high-altitude (HA, above 2,500 m). The potential pathological effects of HA hypobaric hypoxia during gestation have been the focus of several researchers around the world. The studies based on the Himalayan and Central/South American mountains are particularly interesting as these areas account for nearly 70% of the HA world population. At present, studies in human and animal models revealed important alterations in fetal development and growth at HA. Moreover, vascular responses to chronic hypobaria in the pregnant mother and her fetus may induce marked cardiovascular impairments during pregnancy or in the neonatal period. In addition, recent studies have shown potential long-lasting postnatal effects that may increase cardiovascular risk in individuals gestated under chronic hypobaria. Hence, the maternal and fetal adaptive responses to hypoxia, influenced by HA ancestry, are vital for a better developmental and cardiovascular outcome of the offspring. This mini-review exposes and discusses the main determinants of vascular dysfunction due to developmental hypoxia at HA, such as the Andean Mountains, at the maternal and fetal/neonatal levels. Although significant advances have been made from Latin American studies, this area still needs further investigations to reveal the mechanisms involved in vascular dysfunction, to estimate complications of pregnancy and postnatal life adequately, and most importantly, to determine potential treatments to prevent or treat the pathological effects of being developed under chronic hypobaric hypoxia.

Analysis of whole exome sequencing in severe mental illness hints at selection of brain development and immune related genes

Evolutionary trends may underlie some aspects of the risk for common, non-communicable disorders, including psychiatric disease. We analyzed whole exome sequencing data from 80 unique individuals from India coming from families with two or more individuals with severe mental illness. We used Population Branch Statistics (PBS) to identify variants and genes under positive selection and identified 74 genes as candidates for positive selection. Of these, 20 were previously associated with Schizophrenia, Alzheimer’s disease and cognitive abilities in genome wide association studies. We then checked whether any of these 74 genes were involved in common biological pathways or related to specific cellular or molecular functions. We found that immune related pathways and functions related to innate immunity such as antigen binding were over-represented. We also evaluated for the presence of Neanderthal introgressed segments in these genes and found Neanderthal introgression in a single gene out of the 74 candidate genes. However, the introgression pattern indicates the region is unlikely to be the source for selection. Our findings hint at how selection pressures in individuals from families with a history of severe mental illness may diverge from the general population. Further, it also provides insights into the genetic architecture of severe mental illness, such as schizophrenia and its link to immune factors.

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.

Identification of a miRNA-mRNA Regulatory Networks in Placental Tissue Associated With Tibetan High Altitude Adaptation

The Tibetan population has lived and successfully reproduced at high altitude for many generations. Studies have shown that Tibetans have various mechanisms for protection against high-altitude hypoxia, which are probably due, at least in part, to placental adaptation. However, comprehensive in silico analyses of placentas in Tibetans are lacking. We performed a microarray-based comparative transcriptome analysis of 10 Tibetan women from Yushu, Qinghai, CHN (∼3,780 m) and 10 European women living in Leadville, CO, United States (∼3,100 m) for less than three generations. Expression of HIF-1α, STAT3, EGFR, HSP5A, XBP1, and ATF6A mRNA was less in the Tibetan placentas as compared with European placentas. A total of 38 miRNAs were involved in regulating these genes. Differentially expressed genes were enriched for HIF1α signaling pathways, protein processing in the endoplasmic reticulum, PI3K-AKT signaling pathways, and MAPK signaling pathways. Based on the transcriptome profiles, the Tibetan population was distinct from the European population; placental tissues from the Tibetan population are lacking hypoxic responses, and “passivation” occurs in response to hypoxic stress. These results provide insights into the molecular signature of adaptation to high altitudes in these two populations.

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.

Perinatal Hypoxemia and Oxygen Sensing

The development of the control of breathing begins in utero and continues postnatally. Fetal breathing movements are needed for establishing connectivity between the lungs and central mechanisms controlling breathing. Maturation of the control of breathing, including the increase of hypoxia chemosensitivity, continues postnatally. Insufficient oxygenation, or hypoxia, is a major stressor that can manifest for different reasons in the fetus and neonate. Though the fetus and neonate have different hypoxia sensing mechanisms and respond differently to acute hypoxia, both responses prevent deviations to respiratory and other developmental processes. Intermittent and chronic hypoxia pose much greater threats to the normal developmental respiratory processes. Gestational intermittent hypoxia, due to maternal sleep-disordered breathing and sleep apnea, increases eupneic breathing and decreases the hypoxic ventilatory response associated with impaired gasping and autoresuscitation postnatally. Chronic fetal hypoxia, due to biologic or environmental (i.e. high-altitude) factors, is implicated in fetal growth restriction and preterm birth causing a decrease in the postnatal hypoxic ventilatory responses with increases in irregular eupneic breathing. Mechanisms driving these changes include delayed chemoreceptor development, catecholaminergic activity, abnormal myelination, increased astrocyte proliferation in the dorsal respiratory group, among others. Long-term high-altitude residents demonstrate favorable adaptations to chronic hypoxia as do their offspring. Neonatal intermittent hypoxia is common among preterm infants due to immature respiratory systems and thus, display a reduced drive to breathe and apneas due to insufficient hypoxic sensitivity. However, ongoing intermittent hypoxia can enhance hypoxic sensitivity causing ventilatory overshoots followed by apnea; the number of apneas is positively correlated with degree of hypoxic sensitivity in preterm infants. Chronic neonatal hypoxia may arise from fetal complications like maternal smoking or from postnatal cardiovascular problems, causing blunting of the hypoxic ventilatory responses throughout at least adolescence due to attenuation of carotid body fibers responses to hypoxia with potential roles of brainstem serotonin, microglia, and inflammation, though these effects depend on the age in which chronic hypoxia initiates. Fetal and neonatal intermittent and chronic hypoxia are implicated in preterm birth and complicate the respiratory system through their direct effects on hypoxia sensing mechanisms and interruptions to the normal developmental processes. Thus, precise regulation of oxygen homeostasis is crucial for normal development of the respiratory control network. © 2021 American Physiological Society. Compr Physiol 11:1653-1677, 2021.

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.

Hypoxia-inducible factor signaling in pulmonary hypertension

Pulmonary hypertension (PH) is characterized by pulmonary artery remodeling that can subsequently culminate in right heart failure and premature death. Emerging evidence suggests that hypoxia-inducible factor (HIF) signaling plays a fundamental and pivotal role in the pathogenesis of PH. This Review summarizes the regulation of HIF isoforms and their impact in various PH subtypes, as well as the elaborate conditional and cell-specific knockout mouse studies that brought the role of this pathway to light. We also discuss the current preclinical status of pan- and isoform-selective HIF inhibitors, and propose new research areas that may facilitate HIF isoform-specific inhibition as a novel therapeutic strategy for PH and right heart failure.

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.

Pulmonary Hypertension in Acute and Chronic High Altitude Maladaptation Disorders

Alveolar hypoxia is the most prominent feature of high altitude environment with well-known consequences for the cardio-pulmonary system, including development of pulmonary hypertension. Pulmonary hypertension due to an exaggerated hypoxic pulmonary vasoconstriction contributes to high altitude pulmonary edema (HAPE), a life-threatening disorder, occurring at high altitudes in non-acclimatized healthy individuals. Despite a strong physiologic rationale for using vasodilators for prevention and treatment of HAPE, no systematic studies of their efficacy have been conducted to date. Calcium-channel blockers are currently recommended for drug prophylaxis in high-risk individuals with a clear history of recurrent HAPE based on the extensive clinical experience with nifedipine in HAPE prevention in susceptible individuals. Chronic exposure to hypoxia induces pulmonary vascular remodeling and development of pulmonary hypertension, which places an increased pressure load on the right ventricle leading to right heart failure. Further, pulmonary hypertension along with excessive erythrocytosis may complicate chronic mountain sickness, another high altitude maladaptation disorder. Importantly, other causes than hypoxia may potentially underlie and/or contribute to pulmonary hypertension at high altitude, such as chronic heart and lung diseases, thrombotic or embolic diseases. Extensive clinical experience with drugs in patients with pulmonary arterial hypertension suggests their potential for treatment of high altitude pulmonary hypertension. Small studies have demonstrated their efficacy in reducing pulmonary artery pressure in high altitude residents. However, no drugs have been approved to date for the therapy of chronic high altitude pulmonary hypertension. This work provides a literature review on the role of pulmonary hypertension in the pathogenesis of acute and chronic high altitude maladaptation disorders and summarizes current knowledge regarding potential treatment options.

[A Review of High-altitude Hypoxia Adaptation and Hypoxic Solid Tumor]

Historically, the Cambrian explosion was a major life evolution event caused by changes of natural environmental oxygen concentration. The use of oxygen was part of the basic survival instinct of higher life, which evolved a complex regulation system in response to variant levels of oxygen concentration. Hypoxia is one of the typical environmental characteristics in plateau areas. After long-term natural selection in hypoxic conditions, numerous species living in plateau areas have evolved unique mechanisms adapted to hypoxia. Recent studies have found that there are some similarities in adaptation to hypoxia between the animals in highland and different types of human solid tumor cells. Herein, we will summarize recent findings about the hypoxia adaptation evolution in high-altitude animals and the characteristics of hypoxic solid tumors, especially the reactive oxygen species responses in hypoxic solid tumors. We believe that deciphering the underlying molecular mechanisms involved in hypoxia adaptation in highland will facilitate the identification of new genes or biomarkers critical for research on hypoxic solid tumors in the future.

HAHmiR.DB: a server platform for high-altitude human miRNA-gene coregulatory networks and associated regulatory circuits

Around 140 million people live in high-altitude (HA) conditions! and even a larger number visit such places for tourism, adventure-seeking or sports training. Rapid ascent to HA can cause severe damage to the body organs and may lead to many fatal disorders. During induction to HA, human body undergoes various physiological, biochemical, hematological and molecular changes to adapt to the extreme environmental conditions. Several literature references hint that gene-expression-regulation and regulatory molecules like miRNAs and transcription factors (TFs) control adaptive responses during HA stress. These biomolecules are known to interact in a complex combinatorial manner to fine-tune the gene expression and help in controlling the molecular responses during this stress and ultimately help in acclimatization. High-Altitude Human miRNA Database (HAHmiR.DB) is a unique, comprehensive and curated collection of miRNAs that have been experimentally validated to be associated with HA stress, their level of expression in different altitudes, fold change, experiment duration, biomarker association, disease and drug association, tissue-specific expression level, Gene Ontology (GO) and Kyoto Encyclopaedia of Gene and Genomes (KEGG) pathway associations. As a server platform, it also uniquely constructs and analyses interactive miRNA-TF-gene coregulatory networks and extracts regulatory circuits/feed-forward loops (FFLs). These regulatory circuits help to offer mechanistic insights into complex regulatory mechanisms during HA stress. The server can also build these regulatory networks between two and more miRNAs of the database and also identify the regulatory circuits from this network. Hence, HAHmiR.DB is the first-of-its-kind database in HA research, which is a reliable platform to explore, compare, analyse and retrieve miRNAs associated with HA stress, their coregulatory networks and FFL regulatory-circuits. HAHmiR.DB is freely accessible at http://www.hahmirdb.in.

Vascular homeostasis at high-altitude: role of genetic variants and transcription factors

High-altitude pulmonary edema occurs most frequently in non-acclimatized low landers on exposure to altitude ≥2500 m. High-altitude pulmonary edema is a complex condition that involves perturbation of signaling pathways in vasoconstrictors, vasodilators, anti-diuretics, and vascular growth factors. Genetic variations are instrumental in regulating these pathways and evidence is accumulating for a role of epigenetic modification in hypoxic responses. This review focuses on the crosstalk between high-altitude pulmonary edema-associated genetic variants and transcription factors, comparing high-altitude adapted and high-altitude pulmonary edema-afflicted subjects. This approach might ultimately yield biomarker information both to understand and to design therapies for high-altitude adaptation.

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

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

De Novo Transcriptomic and Metabolomic Analyses Reveal the Ecological Adaptation of High-Altitude Bombus pyrosoma

Bombus pyrosoma is one of the most abundant bumblebee species in China, with a distribution range of very varied geomorphology and vegetation, which makes it an ideal pollinator species for research into high-altitude adaptation. Here, we sequenced and assembled transcriptomes of B. pyrosoma from the low-altitude North China Plain and the high-altitude Tibet Plateau. Subsequent comparative analysis of de novo transcriptomes from the high- and low-altitude groups identified 675 common upregulated genes (DEGs) in the high-altitude B. pyrosoma. These genes were enriched in metabolic pathways and corresponded to enzyme activities involved in energy metabolism. Furthermore, according to joint analysis with comparative metabolomics, we suggest that the metabolism of coenzyme A (CoA) and the metabolism and transport of energy resources contribute to the adaptation of high-altitude B. pyrosoma. Meanwhile, we found many common upregulated genes enriched in the Toll and immune deficiency (Imd)signaling pathways that act as important immune defenses in insects, and hypoxia and cold temperatures could induce the upregulation of immune genes in insects. Therefore, we suppose that the Toll and Imd signaling pathways also participated in the high-altitude adaptation of B. pyrosoma. Like other organisms, we suggest that the high-altitude adaptation of B. pyrosoma is controlled by diverse mechanisms.

Genome-Wide Association Study Using Individual Single-Nucleotide Polymorphisms and Haplotypes for Erythrocyte Traits in Alpine Merino Sheep

Adaptation to high-altitude hypoxia is essential for domestic animals, such as yak, Tibetan chicken, and Tibetan sheep, living on high plateaus, as it ensures efficient oxygen absorption and utilization. Red blood cells are the primary medium for transporting oxygen in the blood. However, little is known about the genetic mechanism of erythrocyte traits. Genome-wide association studies (GWASs) based on single markers or haplotypes have identified potential mechanisms for genetic variation and quantitative traits. To identify loci associated with erythrocyte traits, we performed a GWAS based on the method of the single marker and haplotype in 498 Alpine Merino sheep for six erythrocyte traits: red blood cell count (RBC), hemoglobin (HGB), hematocrit (HCT), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), and RBC volume distribution width coefficient of variation (RWD_CV). Forty-two significant single-nucleotide polymorphisms (SNPs) associated with the six erythrocyte traits were detected by means of a single-marker GWAS, and 34 significant haplotypes associated with five erythrocyte traits were detected by means of haplotype analysis. We identified six genes (DHCR24, SPATA9, FLI1, PLCB1, EFNB2, and SH2B3) as potential genes of interest via gene function annotations, location, and expression variation. In particular, FLI1 and PLCB1 were associated with hematopoiesis and erythropoiesis, respectively. These results provide a theoretical basis for analyzing erythrocyte traits and high-altitude hypoxia adaptation in Alpine Merino sheep and will be a useful reference for future studies of plateau-dwelling livestock.

Transcriptomic modifications in developmental cardiopulmonary adaptations to chronic hypoxia using a murine model of simulated high-altitude exposure

Mechanisms driving adaptive developmental responses to chronic high-altitude (HA) exposure are incompletely known. We developed a novel rat model mimicking the human condition of cardiopulmonary adaptation to HA starting at conception and spanning the in utero and postnatal timeframe. We assessed lung growth and cardiopulmonary structure and function and performed transcriptome analyses to identify mechanisms facilitating developmental adaptations to chronic hypoxia. To generate the model, breeding pairs of Sprague-Dawley rats were exposed to hypobaric hypoxia (equivalent to 9,000 ft elevation). Mating, pregnancy, and delivery occurred in hypoxic conditions. Six weeks postpartum, structural and functional data were collected in the offspring. RNA-Seq was performed on right ventricle (RV) and lung tissue. Age-matched breeding pairs and offspring under room air (RA) conditions served as controls. Hypoxic rats exhibited significantly lower body weights and higher hematocrit levels, alveolar volumes, pulmonary diffusion capacities, RV mass, and RV systolic pressure, as well as increased pulmonary artery remodeling. RNA-Seq analyses revealed multiple differentially expressed genes in lungs and RVs from hypoxic rats. Although there was considerable similarity between hypoxic lungs and RVs compared with RA controls, several upstream regulators unique to lung or RV were identified. We noted a pattern of immune downregulation and regulation patterns of immune and hormonal mediators similar to the genome from patients with pulmonary arterial hypertension. In summary, we developed a novel murine model of chronic hypoxia exposure that demonstrates functional and structural phenotypes similar to human adaptation. We identified transcriptomic alterations that suggest potential mechanisms for adaptation to chronic HA.

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.