Genetic variations in Tibetan populations and high-altitude adaptation at the Himalayas

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.

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.

Comparative mitogenomic analysis of subterranean and surface amphipods (Crustacea, Amphipoda) with special reference to the family Crangonyctidae

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes.

Hypobaric hypoxia drives selection of altitude-associated adaptative alleles in the Himalayan population

Genetic variants play a crucial role in shaping the adaptive phenotypes associated with high-altitude populations. Nevertheless, a comprehensive understanding of the specific impacts of different environments associated with increasing altitudes on the natural selection of these genetic variants remains undetermined. Hence, this study aimed to identify genetic markers responsible for high-altitude adaptation with specific reference to different altitudes, majorly focussing on an altitude elevation range of ∼1500 m and a corresponding decrease of ≥5 % in ambient oxygen availability. We conducted a comprehensive genome-wide investigation (n = 192) followed by a validation study (n = 514) in low-altitude and three high-altitude populations (>2400 m) of Nubra village (NU) (3048 m), Sakti village (SKT) (3812 m), and Tso Moriri village (TK) (4522 m). Extensive genetic analysis identified 86 SNPs that showed significant associations with high-altitude adaptation. Frequency mapping of these SNPs revealed 38 adaptive alleles and specific haplotypes that exhibited a strong linear correlation with increasing altitude. Notably, these SNPs spanned crucial genes, such as ADH6 and NAPG along with the vastly studied genes like EGLN1 and EPAS1, involved in oxygen sensing, metabolism, and vascular homeostasis. Correlation analyses between these adaptive alleles and relevant clinical and biochemical markers provided evidence of their functional relevance in physiological adaptation to hypobaric hypoxia. High-altitude population showed a significant increase in plasma 8-isoPGF2α levels as compared to low-altitude population. Similar observation showcased increased blood pressure in NU as compared to TK (P < 0.0001). In silico analyses further confirmed that these alleles regulate gene expression of EGLN1, EPAS1, COQ7, NAPG, ADH6, DUOXA1 etc. This study provides genetic insights into the effects of hypobaric-hypoxia on the clinico-physiological characteristics of natives living in increasing high-altitude regions. Overall, our findings highlight the synergistic relationship between environment and evolutionary processes, showcasing physiological implications of genetic variants in oxygen sensing and metabolic pathway genes in increasing high-altitude environments.

Comparative Analysis of Chloroplast Genome of Meconopsis (Papaveraceae) Provides Insights into Their Genomic Evolution and Adaptation to High Elevation

The Meconopsis species are widely distributed in the Qinghai-Tibet Plateau, Himalayas, and Hengduan Mountains in China, and have high medicinal and ornamental value. The high diversity of plant morphology in this genus poses significant challenges for species identification, given their propensity for highland dwelling, which makes it a question worth exploring how they cope with the harsh surroundings. In this study, we recently generated chloroplast (cp) genomes of two Meconopsis species, Meconopsis paniculata (M. paniculata) and M. pinnatifolia, and compared them with those of ten Meconopsis cp genomes to comprehend cp genomic features, their phylogenetic relationships, and what part they might play in plateau adaptation. These cp genomes shared a great deal of similarities in terms of genome size, structure, gene content, GC content, and codon usage patterns. The cp genomes were between 151,864 bp and 154,997 bp in length, and contain 133 predictive genes. Through sequence divergence analysis, we identified three highly variable regions (trnD-psbD, ccsA-ndhD, and ycf1 genes), which could be used as potential markers or DNA barcodes for phylogenetic analysis. Between 22 and 38 SSRs and some long repeat sequences were identified from 12 Meconopsis species. Our phylogenetic analysis confirmed that 12 species of Meconopsis clustered into a monophyletic clade in Papaveraceae, which corroborated their intrageneric relationships. The results indicated that M. pinnatifolia and M. paniculata are sister species in the phylogenetic tree. In addition, the atpA and ycf2 genes were positively selected in high-altitude species. The functions of these two genes might be involved in adaptation to the extreme environment in the cold and low CO2 concentration conditions at the plateau.

Functional EPAS1/HIF2A missense variant is associated with hematocrit in Andean highlanders

Hypoxia-inducible factor pathway genes are linked to adaptation in both human and nonhuman highland species. EPAS1, a notable target of hypoxia adaptation, is associated with relatively lower hemoglobin concentration in Tibetans. We provide evidence for an association between an adaptive EPAS1 variant (rs570553380) and the same phenotype of relatively low hematocrit in Andean highlanders. This Andean-specific missense variant is present at a modest frequency in Andeans and absent in other human populations and vertebrate species except the coelacanth. CRISPR-base-edited human cells with this variant exhibit shifts in hypoxia-regulated gene expression, while metabolomic analyses reveal both genotype and phenotype associations and validation in a lowland population. Although this genocopy of relatively lower hematocrit in Andean highlanders parallels well-replicated findings in Tibetans, it likely involves distinct pathway responses based on a protein-coding versus noncoding variants, respectively. These findings illuminate how unique variants at EPAS1 contribute to the same phenotype in Tibetans and a subset of Andean highlanders despite distinct evolutionary trajectories.

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.

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.

Structural variants involved in high-altitude adaptation detected using single-molecule long-read sequencing

Structural variants (SVs), accounting for a larger fraction of the genome than SNPs/InDels, are an important pool of genetic variation, enabling environmental adaptations. Here, we perform long-read sequencing data of 320 Tibetan and Han samples and show that SVs are highly involved in high-altitude adaptation. We expand the landscape of global SVs, apply robust models of selection and population differentiation combining SVs, SNPs and InDels, and use epigenomic analyses to predict enhancers, target genes and biological functions. We reveal diverse Tibetan-specific SVs affecting the regulatory circuitry of biological functions, including the hypoxia response, energy metabolism and pulmonary function. We find a Tibetan-specific deletion disrupts a super-enhancer and downregulates EPAS1 using enhancer reporter, cellular knock-out and DNA pull-down assays. Our study expands the global SV landscape, reveals the role of gene-regulatory circuitry rewiring in human adaptation, and illustrates the diverse functional roles of SVs in human biology.

Hypoxia-induced signaling in the cardiovascular system: pathogenesis and therapeutic targets

Hypoxia, characterized by reduced oxygen concentration, is a significant stressor that affects the survival of aerobic species and plays a prominent role in cardiovascular diseases. From the research history and milestone events related to hypoxia in cardiovascular development and diseases, The "hypoxia-inducible factors (HIFs) switch" can be observed from both temporal and spatial perspectives, encompassing the occurrence and progression of hypoxia (gradual decline in oxygen concentration), the acute and chronic manifestations of hypoxia, and the geographical characteristics of hypoxia (natural selection at high altitudes). Furthermore, hypoxia signaling pathways are associated with natural rhythms, such as diurnal and hibernation processes. In addition to innate factors and natural selection, it has been found that epigenetics, as a postnatal factor, profoundly influences the hypoxic response and progression within the cardiovascular system. Within this intricate process, interactions between different tissues and organs within the cardiovascular system and other systems in the context of hypoxia signaling pathways have been established. Thus, it is the time to summarize and to construct a multi-level regulatory framework of hypoxia signaling and mechanisms in cardiovascular diseases for developing more therapeutic targets and make reasonable advancements in clinical research, including FDA-approved drugs and ongoing clinical trials, to guide future clinical practice in the field of hypoxia signaling in cardiovascular diseases.

Structure and Metabolic Characteristics of Intestinal Microbiota in Tibetan and Han Populations of Qinghai-Tibet Plateau and Associated Influencing Factors

Residents of the Qinghai-Tibet Plateau might experience shifts in their gut microbiota composition as a result of the plateau environment. For example, high altitudes can increase the abundance of obligate anaerobic bacteria, decrease the number of aerobic bacteria and facultative anaerobic bacteria, increase probiotics, and decrease pathogenic bacteria. This study aimed to determine the structure and metabolic differences in intestinal microbial communities among the Tibetan and Han populations on the Qinghai-Xizang Plateau and shed light on the factors that influence the abundance of the microbial communities in the gut. The structural characteristics of intestinal microorganisms were detected from blood and fecal samples using 16S rRNA sequencing. Metabolic characteristics were detected using gas chromatography-time-of-flight mass spectrometry (GC-TOFMS). The influencing factors were analyzed using Spearman's correlation analysis. Bacteroides and Bifidobacterium were dominant in the intestinal tract of the Han population, while Bacteroides and Prevotella were dominant in that of the Tibetan population, with marked differences in Pseudomonas, Prevotella, and other genera. Ferulic acid and 4-methylcatechol were the main differential metabolites between the Tibetan and Han ethnic groups. This may be the reason for the different adaptability of Tibetan and Han nationalities to the plateau. Alanine aminotransferase and uric acid also have a high correlation with different bacteria and metabolites, which may play a role. These results reveal notable disparities in the compositions and metabolic characteristics of gut microbial communities in the Tibetan and Han people residing on the Qinghai-Tibet Plateau and may provide insights regarding the mechanism of plateau adaptability.

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

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

Tensor Decomposition-based Feature Extraction and Classification to Detect Natural Selection from Genomic Data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under nonconvex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data although preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx, which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

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.

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.

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

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

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

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

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00117-x.

Searching for molecular hypoxia sensors among oxygen-dependent enzymes

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

Lipid levels in the Jiarong Tibetan's diet at high altitudes: a cross-sectional survey

Despite the ongoing debate on the inconsistent and controversial effects of Tibetan diet on blood lipid levels at high altitude, this cross-sectional study was conducted to analyze the relationship between dietary practices and blood lipid levels among Jiarong Tibetan population. A total of 476 Jiarong Tibetan residents were included, in which basic demographic data, physical activity records, simplified food frequency questionnaire, and biochemical data were collected. Using multivariate logistic regression analysis, the potential associations between the variables were examined, and it was found that fat energy supply ratio increased with the elevation of altitude, while the lipid level showed an inverted U-shaped variation. However, the findings suggested that a diet rich in unsaturated fatty acids might balance the effects of the Tibetan diet on the risk of lipid metabolism disorders. Therefore, it is crucial to concentrate on the fat composition rather than the amount of fat E% intake on the plateau. The results highlighted the importance of investigating the interaction between environment and genes in lipid levels among plateau Tibetan population. However, further large-scale prospective studies are required for better understanding of the complexities involved in dietary practices and their influences on blood lipid levels.

Oxygen stress on age-stage, two-sex life tables and transcriptomic response of diamondback moth (Plutella xylostella)

Elucidating the genetic basis of local adaption is one of the important tasks in evolutionary biology. The Qinghai-Tibet Plateau has the highest biodiversity for an extreme environment worldwide, and provides an ideal natural laboratory to study adaptive evolution. The diamondback moth (DBM), Plutella xylostella, is one of the most devastating pests of the global Brassica industry. A highly heterozygous genome of this pest has facilitated its adaptation to a variety of complex environments, and so provides an ideal model to study fast adaptation. We conducted a pilot study combining RNA-seq with an age-stage, two-sex life table to study the effects of oxygen deprivation on DBM. The developmental periods of all instars were significantly shorter in the hypoxic environment. We compared the transcriptomes of DBM from Fuzhou, Fujian (low-altitude) and Lhasa, Tibet (high-altitude) under hypoxia treatment in a hypoxic chamber. Some DEGs are enriched in pathways associated with DNA replication, such as DNA repair, nucleotide excision repair, base excision repair, mismatch repair and homologous recombination. The pathways with significant changes were associated with metabolism process and cell development. Thus, we assumed that insects could adapt to different environments by regulating their metabolism. Our findings indicated that although adaptive mechanisms to hypoxia in different DBM strains could be similar, DBM individuals from Tibet had superior tolerance to hypoxia compared with those of Fuzhou. Local adaptation of the Tibetan colony was assumed to be responsible for this difference. Our research suggests novel mechanisms of insect responses to hypoxia stress.

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.

Backbone phylogeny and adaptive evolution of Pleurospermum s. l.: New insights from phylogenomic analyses of complete plastome data

Pleurospermum is a taxonomically challenging taxon of Apiaceae, as its circumscription and composition remain controversial for morphological similarities with several related genera, leading to a dispute between Pleurospermum in the broad sense and strict sense. While evidence from previous molecular studies recognized plural branching lineages within the Pleurospermum s. l., it did not support the latest delimitation of Pleurospermum s. str. by only two closely related northern species. So far, no proper delimitation for Pleurospermum has come up, and many of the plural taxa in Pleurospermum s. l. remain unresolved, which may be due to poor phylogenetic resolution yielded barely from ITS sequences. Herein, we newly assembled 40 complete plastomes from 36 species of Pleurospermum s. l. and related genera, 34 of which were first reported and generated a well-resolved backbone phylogeny in a framework of the subfamily Apioideae. From the phylogeny with greatly improved resolution, a total of six well-supported monophyletic lineages within Pleurospermum s. l. were recognized falling in different major clades of Apioideae. Combining morphological characteristics with phylogenetic inference, we suggested to re-delimit the Pleurospermum s. str. by introducing nine species mainly from the Himalayan regions and proposed its boundary features; the remaining species were suggested to be excluded from Pleurospermum to incorporate into their more related taxa being revealed. On this basis, the plastome comparison revealed not only the high conservatism but also the mild differences among lineages in plastome structure and gene evolution. Overall, our study provided a backbone phylogeny essential for further studies of the taxonomically difficult taxa within Pleurospermum s. l.

Tensor decomposition based feature extraction and classification to detect natural selection from genomic data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under non-convex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data while preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx , which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

Genomic evidence reveals intraspecific divergence of the hot-spring snake (Thermophis baileyi), an endangered reptile endemic to the Qinghai-Tibet plateau

Understanding how and why species evolve requires knowledge on intraspecific divergence. In this study, we examined intraspecific divergence in the endangered hot-spring snake (Thermophis baileyi), an endemic species on the Qinghai-Tibet Plateau (QTP). Whole-genome resequencing of 58 sampled individuals from 15 populations was performed to identify the drivers of intraspecific divergence and explore the potential roles of genes under selection. Our analyses resolved three groups, with major intergroup admixture occurring in regions of group contact. Divergence probably occurred during the Pleistocene as a result of glacial climatic oscillations, Yadong-Gulu rift, and geothermal fields differentiation, while complex gene flow between group pairs reflected a unique intraspecific divergence pattern on the QTP. Intergroup fixed loci involved selected genes functionally related to divergence and local adaptation, especially adaptation to hot spring microenvironments in different geothermal fields. Analysis of structural variants, genetic diversity, inbreeding, and genetic load indicated that the hot-spring snake population has declined to a low level with decreased diversity, which is important for the conservation management of this endangered species. Our study demonstrated that the integration of demographic history, gene flow, genomic divergence genes, and other information is necessary to distinguish the evolutionary processes involved in speciation.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Regional variation in growth status: A cross-sectional survey among Tibetan adolescents living at three different high altitudes in Tibet, China

OBJECTIVE

To explore the regional and age variation in growth status and age-by-site interaction effect on the growth status among Chinese Tibetan adolescents at different altitudes in Tibet, China.

METHODS

The research was conducted in three regions of Tibet, China: Nyingchi (average altitude 3100 m), Lhasa (average altitude 3650 m), and Nagqu (average altitude 4500 m). A total of 3817 Chinese Tibetan adolescents aged 12-18 years were tested for height, weight, chest circumference, and waist circumference. One-way ANOVA was used to compare the growth status of Chinese Tibetan adolescents. Two-way ANOVA was used to explore the age-by-site interaction effect on the growth status of Chinese Tibetan adolescents.

RESULTS

The height, weight, and chest circumference of Chinese Tibetan adolescents in Nagqu are the lowest among the three cities. Age and site have an interaction effect influence on the growth status of Chinese Tibetan adolescents (p < .01).

CONCLUSIONS

The growth indicators (height, weight, chest circumference, WC) of Chinese Tibetan adolescents differed with altitudes. Policies to improve the growth status of Chinese Tibetan adolescents in Nagqu are urgently needed.

Review of the relationship and underlying mechanisms between the Qinghai-Tibet plateau and host intestinal flora

The intestinal microbial community is the largest ecosystem in the human body, in which the intestinal flora plays a dominant role and has a wide range of biological functions. However, it is vulnerable to a variety of factors, and exposure to extreme environments at high altitudes, as seen on the Qinghai-Tibet plateau, may cause changes in the structure and function of the host intestinal flora. Conversely, the intestinal flora can help the host adapt to the plateau environment through a variety of ways. Herein, we review the relationship and underlying mechanism between the host intestinal flora and the plateau environment by discussing the characteristics of the plateau environment, its influence on the intestinal flora, and the important role of the intestinal flora in host adaptation to the plateau environment. This review aimed to provide a reference for maintaining the health of the plateau population.

Risk factors and SCN5A-H558R polymorphism for atrial fibrillation in Tibetans living at different altitudes

Several studies have found associations of genes with atrial fibrillation (AF), including SCN5A-H558R. However, there are limited data of these associations among populations living at different altitudes. We investigated the relationship between the SCN5A-H558R polymorphism and AF in Tibetans living at different altitudes in Qinghai, China. General clinical and genotype data were obtained from 72 patients with AF and 109 non-AF (NAF) individuals at middle altitudes, and from 102 patients with AF and 143 NAF individuals at high altitudes. Multifactor logistic regression was performed to determine associations and AF risk factors. SCN5A-H558R genotypes differed significantly between the AF and NAF groups (P < .0125) and the G allele was an independent AF risk factor (P < .05) at both altitudes, with no significant differences according to altitude (P > .0125). At middle altitudes, age, red blood cell distribution width (RDW-SD), left atrial internal diameter (LAD), and G allele were independent AF risk factors. At high altitudes, age, smoking, hypertension, RDW-SD, free triiodothyronine, LAD, and G allele were independent AF risk factors (P < .05). The G allele of SCN5A-H558R might be an independent risk factor of AF both high and middle altitude, but there are some differences in other clinical risk factors of AF.

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.

Comparative Analysis of Chloroplast Genomes within Saxifraga (Saxifragaceae) Takes Insights into Their Genomic Evolution and Adaption to the High-Elevation Environment

Saxifraga species are widely distributed in alpine and arctic regions in the Northern hemisphere. Highly morphological diversity within this genus brings great difficulties for species identification, and their typical highland living properties make it interesting how they adapt to the extreme environment. Here, we newly generated the chloroplast (cp) genomes of two Saxifraga species and compared them with another five Saxifraga cp genomes to understand the characteristics of cp genomes and their potential roles in highland adaptation. The genome size, structure, gene content, GC content, and codon usage pattern were found to be highly similar. Cp genomes ranged from 146,549 bp to 151,066 bp in length, most of which comprised 130 predicted genes. Yet, due to the expansion of IR regions, the second copy of rps19 in Saxifraga stolonifera was uniquely kept. Through sequence divergence analysis, we identified seven hypervariable regions and detected some signatures of regularity associated with genetic distance. We also identified 52 to 89 SSRs and some long repeats among seven Saxifraga species. Both ML and BI phylogenetic analyses confirmed that seven Saxifraga species formed a monophyletic clade in the Saxifragaceae family, and their intragenus relationship was also well supported. Additionally, the ndhI and ycf1 genes were considered under positive selection in species inhabiting relatively high altitudes. Given the conditions of intense light and low CO₂ concentration in the highland, the products of these two genes might participate in the adaptation to the extreme environment.

Chromosome-level Genome Assembly of the High-altitude Leopard (Panthera pardus) Sheds Light on Its Environmental Adaptation

The leopard (Panthera pardus) has the largest natural distribution from low- to high-altitude areas of any wild felid species, but recent studies have revealed that leopards have disappeared from large areas, probably owing to poaching, a decline of prey species, and habitat degradation. Here, we reported the chromosome-scale genome assembly of the high-altitude leopard (HL) based on nanopore sequencing and high-throughput chromatin conformation capture (Hi-C) technology. Panthera genomes revealed similar repeat composition, and there was an appreciably conserved synteny between HL and the other two Panthera genomes. Divergence time analysis based on the whole genomes revealed that the HL and the low-altitude leopard differentiate from a common ancestor ∼2.2 Ma. Through comparative genomics analyses, we found molecular genetic signatures that may reflect high-altitude adaptation of the HL. Three HL-specific missense mutations were detected in two positively selected genes, that is, ITGA7 (Ala112Gly, Asp113Val, and Gln115Pro) and NOTCH2 (Ala2398Ser), which are likely to be associated with hypoxia adaptation. The chromosome-level genome of the HL provides valuable resources for the investigation of high-altitude adaptation and protection management of the vulnerable leopard.

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.

Human adaptation to high altitude: a review of convergence between genomic and proteomic signatures

Both genomics- and proteomics-based investigations have identified several essential genes, proteins, and pathways that may facilitate human adaptive genotype/phenotype in a population-specific manner. This comprehensive review provides an up-to-date list of genes and proteins identified for human adaptive responses to high altitudes. Genomics studies for indigenous high-altitude populations like Tibetans, Andeans, Ethiopians, and Sherpas have identified 169 genes under positive natural selection. Similarly, global proteomics studies have identified 258 proteins (± 1.2-fold or more) for Tibetan, Sherpa, and Ladakhi highlanders. The primary biological processes identified for genetic signatures include hypoxia-inducible factor (HIF)-mediated oxygen sensing, angiogenesis, and erythropoiesis. In contrast, major biological processes identified for proteomics signatures include 14–3-3 mediated sirtuin signaling, integrin-linked kinase (ILK), phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT), and integrin signaling. Comparing genetic and protein signatures, we identified 7 common genes/proteins (HBB/hemoglobin subunit beta, TF/serotransferrin, ANGPTL4/angiopoietin-related protein 4, CDC42/cell division control protein 42 homolog, GC/vitamin D-binding protein, IGFBP1/insulin-like growth factor-binding protein 1, and IGFBP2/insulin-like growth factor-binding protein 2) involved in crucial molecular functions like IGF-1 signaling, LXR/RXR activation, ferroptosis signaling, iron homeostasis signaling and regulation of cell cycle. Our combined multi-omics analysis identifies common molecular targets and pathways for human adaptation to high altitude. These observations further corroborate convergent positive selection of hypoxia-responsive molecular pathways in humans and advocate using multi-omics techniques for deciphering human adaptive responses to high altitude.

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.

What to do with low O2: Redox adaptations in vertebrates native to hypoxic environments

Reactive oxygen species (ROS) are important cellular signalling molecules but sudden changes in redox balance can be deleterious to cells and lethal to the whole organism. ROS production is inherently linked to environmental oxygen availability and many species live in variable oxygen environments that can range in both severity and duration of hypoxic exposure. Given the importance of redox homeostasis to cell and animal viability, it is not surprising that early studies in species adapted to various hypoxic niches have revealed diverse strategies to limit or mitigate deleterious ROS changes. Although research in this area is in its infancy, patterns are beginning to emerge in the suites of adaptations to different hypoxic environments. This review focuses on redox adaptations (i.e., modifications of ROS production and scavenging, and mitigation of oxidative damage) in hypoxia-tolerant vertebrates across a range of hypoxic environments. In general, evidence suggests that animals adapted to chronic lifelong hypoxia are in homeostasis, and do not encounter major oxidative challenges in their homeostatic environment, whereas animals exposed to seasonal chronic anoxia or hypoxia rapidly downregulate redox balance to match a hypometabolic state and employ robust scavenging pathways during seasonal reoxygenation. Conversely, animals adapted to intermittent hypoxia exposure face the greatest degree of ROS imbalance and likely exhibit enhanced ROS-mitigation strategies. Although some progress has been made, research in this field is patchy and further elucidation of mechanisms that are protective against environmental redox challenges is imperative for a more holistic understanding of how animals survive hypoxic environments.

Denisovans and Homo sapiens on the Tibetan Plateau: dispersals and adaptations

Recent archaeological discoveries suggest that both archaic Denisovans and Homo sapiens occupied the Tibetan Plateau earlier than expected. Genetic studies show that a pulse of Denisovan introgression was involved in the adaptation of Tibetan populations to high-altitude hypoxia. These findings challenge the traditional view that the plateau was one of the last places on earth colonized by H. sapiens and warrant a reappraisal of the population history of this highland. Here, we integrate archaeological and genomic evidence relevant to human dispersal, settlement, and adaptation in the region. We propose two testable models to address the peopling of the plateau in the broader context of H. sapiens dispersal and their encounters with Denisovans in Asia.

Lung transcriptome analysis for the identification of genes involved in the hypoxic adaptation of plateau pika (Ochotona curzoniae)

The plateau pika, a typical hypoxia-tolerant mammal lives 3000-5000 m above sea level on the Qinghai-Tibet Plateau, has acquired many physiological and morphological characteristics and strategies in its adaptation to sustained, high-altitude hypoxia. Blunted hypoxic pulmonary vasoconstriction is one such strategy, but the genes involved in this strategy have not been elucidated. Here, we investigated the genes involved and their expression profiles in the lung transcriptome of plateau pikas subjected to different hypoxic conditions (using low-pressure oxygen cabins). A slight, right ventricular hypertrophy was observed in pikas of the control group (altitude: 3200 m) vs. those exposed to 5000 m altitude conditions for one week. Our assembly identified 67,774 genes; compared with their expression in the control animals, 866 and 8364 genes were co-upregulated and co-downregulated, respectively, in pikas subjected to 5000 m altitude conditions for 1 and 4 w. We elucidated pathways that were associated with pulmonary vascular arterial pressure, including vascular smooth muscle contraction, HIF-1 signalling, calcium signalling, cGMP-PKG signalling, and PI3K-Akt signalling based on the differentially expressed genes; the top-100 pathway enrichments were found between the control group and the group exposed to 5000 m altitude conditions for 4 w. The mRNA levels of 18 candidate gene showed that more than 83% of genes were expressed and the number of transcriptome The up-regulated genes were EPAS1, Hbα, iNOS, CX40, CD31, PPM1B, HIF-1α, MYLK, Pcdh12, Surfactant protein B, the down-regulated genes were RYR2, vWF, RASA1, CLASRP, HIF-3α. Our transcriptome data are a valuable resource for future genomic studies on plateau pika.

Targeting HIF-2α in the Tumor Microenvironment: Redefining the Role of HIF-2α for Solid Cancer Therapy

Inadequate oxygen supply, or hypoxia, is characteristic of the tumor microenvironment and correlates with poor prognosis and therapeutic resistance. Hypoxia leads to the activation of the hypoxia-inducible factor (HIF) signaling pathway and stabilization of the HIF-α subunit, driving tumor progression. The homologous alpha subunits, HIF-1α and HIF-2α, are responsible for mediating the transcription of a multitude of critical proteins that control proliferation, angiogenic signaling, metastasis, and other oncogenic factors, both differentially and sequentially regulating the hypoxic response. Post-translational modifications of HIF play a central role in its behavior as a mediator of transcription, as well as the temporal transition from HIF-1α to HIF-2α that occurs in response to chronic hypoxia. While it is evident that HIF-α is highly dynamic, HIF-2α remains vastly under-considered. HIF-2α can intensify the behaviors of the most aggressive tumors by adapting the cell to oxidative stress, thereby promoting metastasis, tissue remodeling, angiogenesis, and upregulating cancer stem cell factors. The structure, function, hypoxic response, spatiotemporal dynamics, and roles in the progression and persistence of cancer of this HIF-2α molecule and its EPAS1 gene are highlighted in this review, alongside a discussion of current therapeutics and future directions.

The genomic basis of high-elevation adaptation in wild house mice (Mus musculus domesticus) from South America

Understanding the genetic basis of environmental adaptation in natural populations is a central goal in evolutionary biology. The conditions at high elevation, particularly the low oxygen available in the ambient air, impose a significant and chronic environmental challenge to metabolically active animals with lowland ancestry. To understand the process of adaptation to these novel conditions and to assess the repeatability of evolution over short timescales, we examined the signature of selection from complete exome sequences of house mice (Mus musculus domesticus) sampled across two elevational transects in the Andes of South America. Using phylogenetic analysis, we show that house mice colonized high elevations independently in Ecuador and Bolivia. Overall, we found distinct responses to selection in each transect and largely nonoverlapping sets of candidate genes, consistent with the complex nature of traits that underlie adaptation to low oxygen availability (hypoxia) in other species. Nonetheless, we also identified a small subset of the genome that appears to be under parallel selection at the gene and SNP levels. In particular, three genes (Col22a1, Fgf14, and srGAP1) bore strong signatures of selection in both transects. Finally, we observed several patterns that were common to both transects, including an excess of derived alleles at high elevation, and a number of hypoxia-associated genes exhibiting a threshold effect, with a large allele frequency change only at the highest elevations. This threshold effect suggests that selection pressures may increase disproportionately at high elevations in mammals, consistent with observations of some high-elevation diseases in humans.

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

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

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s43657-021-00038-7.

Genetic and immune changes in Tibetan high-altitude populations contribute to biological adaptation to hypoxia

BACKGROUND

Tibetans have lived at very high altitudes for thousands of years, and have a distinctive suite of physiological traits that enable them to tolerate environmental hypoxia. Expanding awareness and knowledge of the differences in hematology, hypoxia-associated genes, immune system of people living at different altitudes and from different ethnic groups may provide evidence for the prevention of mountain sickness.

METHOD

Ninety-five Han people at mid-altitude, ninety-five Tibetan people at high-altitude and ninety-eight Han people at high-altitude were recruited. Red blood cell parameters, immune cells, the contents of cytokines, hypoxia-associated gene single nucleotide polymorphisms (SNPs) were measured.

RESULTS

The values of Hematocrit (HCT), Mean cell volume (MCV) and Mean cell hemoglobin (MCH) in red blood cell, immune cell CD19⁺ B cell number, the levels of cytokines Erb-B2 receptor tyrosine kinase 3 (ErbB3) and Tumor necrosis factor receptor II (TNF-RII) and the levels of hypoxia-associated factors Hypoxia inducible factor-1α (HIF-1α), Hypoxia inducible factor-2α (HIF-2α) and HIF prolyl 4-hydroxylase 2 (PHD2) were decreased, while the frequencies of SNPs in twenty-six Endothelial PAS domain protein 1 (EPAS1) and Egl-9 family hypoxia inducible factor 1 (EGLN1) were increased in Tibetan people at high-altitude compared with that of Han peoples at high-altitude. Furthermore, compared with mid-altitude individuals, high-altitude individuals showed lower blood cell parameters including Hemoglobin concentration (HGB), HCT, MCV and MCH, higher Mean cell hemoglobin concentration (MCHC), lower immune cells including CD19⁺ B cells, CD4⁺ T cells and CD4/CD8 ratio, higher immune cells containing CD8⁺ T cells and CD16/56NK cells, decreased Growth regulated oncogene alpha (GROa), Macrophage inflammatory protein 1 beta (MIP-1b), Interleukin-8 (IL-8), and increased Thrombomodulin, downregulated hypoxia-associated factors including HIF1α, HIF2α and PHD2, and higher frequency of EGLN1 rs2275279.

CONCLUSIONS

These results indicated that biological adaption to hypoxia at high altitude might have been mediated by changes in immune cells, cytokines, and hypoxia-associated genes during the evolutionary history of Tibetan populations. Furthermore, different responses to high altitude were observed in different ethnic groups, which may provide a useful knowledge to improve the protection of high-altitude populations from mountain sickness.

Genetic characterization of the highlander Tibetan population from Qinghai-Tibet Plateau revealed by X chromosomal STRs

Tibetans are considered an East Asian ethnic group and primarily live in the high Tibetan plateau, the western Sichuan and Yunnan mountains of central and southern China, and areas throughout the Himalayas and around the Tibetan plateau. These people exhibit rare molecular machinery that allows them to adapt to hypoxic environments in the Qinghai-Tibet Plateau and make them a potential candidate for providing insights related to medical genetic, molecular medicine and human population studies. In the current study, we have genotyped 549 individuals with Investigator Argus X-12 Kit. For 12 X-STRs, a total of 174 unique alleles were found, among them DXS10134 and DXS10135 were the most polymorphic loci. All of the loci were in Hardy-Weinberg Equilibrium (HWE). The numbers of observed haplotypes in Highlander Tibetans males were 161,112, 96 and 108, respectively, whereas haplotype diversities (HD) were 0.9959, 0.9880, 0.9809 and 0.9873, respectively. The combined discrimination power for males (PDm) was 0.999 999 99701 and for females (PDf) was 0.999 999 999 999 9958. This study represents an extensive report on X chromosomal STR markers variation in the Highlander Tibetans population for forensic applications and population genetic studies.

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.

High-quality evergreen azalea genome reveals tandem duplication-facilitated low-altitude adaptability and floral scent evolution

Azalea belongs to Rhododendron, which is one of the largest genera of flowering plants and is well known for the diversity and beauty in its more than 1000 woody species. Rhododendron contains two distinct groups: the most high-altitude and a few low-altitude species; however, the former group is difficult to be domesticated for urban landscaping, and their evolution and adaptation are little known. Rhododendron ovatum has broad adaptation in low-altitude regions but possesses evergreen characteristics like high-altitude species, and it has floral fragrance that is deficient in most cultivars. Here we report the chromosome-level genome assembly of R. ovatum, which has a total length of 549 Mb with scaffold N50 of 41 Mb and contains 41 264 predicted genes. Genomic micro-evolutionary analysis of R. ovatum in comparison with two high-altitude Rhododendron species indicated that the expansion genes in R. ovatum were significantly enriched in defence responses, which may account for its adaptability in low altitudes. The R. ovatum genome contains much more terpene synthase genes (TPSs) compared with the species that lost floral fragrance. The subfamily b members of TPS are involved in the synthesis of sesquiterpenes as well as monoterpenes and play a major role in flora scent biosynthesis and defence responses. Tandem duplication is the primary force driving expansion of defence-responsive genes for extensive adaptability to the low-altitude environments. The R. ovatum genome provides insights into low-altitude adaptation and gain or loss of floral fragrance for Rhododendron species, which are valuable for alpine plant domestication and floral scent breeding.

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.

FGF5 and EPAS1 gene polymorphisms are associated with high-altitude adaptation in Nepalese goat breeds

Several studies have reported the gene polymorphisms associated with high-altitude adaptation in goats. The FGF5 gene is a regulator in the hair-growth and a SNP c.-253G>A located within 5'UTR has been reported to cause long-haired phenotype. The EPAS1 gene is a transcription factor for various genes that have hypoxia-adaptive functions and a nonsynonymous SNP (Q579L) located in exon 5 has been reported to be associated with the mean corpuscular hemoglobin concentration. Nepal has large difference in altitudes in the north-south direction and four indigenous goat breeds are bred depending on the altitude. We used a total of 130 animals in Nepal, Chyangra (n = 37), Sinhal (n = 24), Khari (n = 33), and Terai (n = 36), and genotyped these two gene polymorphisms to compare the gene frequencies among the breeds and investigate the associations between breeding altitudes and allele frequencies. The genotyping results revealed that the mutant allele frequency in both polymorphisms tended to increase, as the breeding altitude of each population increased. In addition, correlation coefficients showed a relatively strong positive correlation between the breeding altitude and the mutant allele frequencies (r = 0.87 in FGF5 and r = 0.68 in EPAS1). These results suggested that both polymorphisms would significantly contribute to the high-altitude adaptation in Nepalese goat breeds.

Genomic and functional evidence reveals convergent evolution in fishes on the Tibetan Plateau

High-altitude environments are strong drivers of adaptive evolution in endemic organisms. However, little is known about the genetic mechanisms of convergent adaptation among different lineages, especially in fishes. There are three independent fish groups on the Tibetan Plateau: Tibetan Loaches, Schizothoracine fishes and Glyptosternoid fishes; all are well adapted to the harsh environmental conditions. They represent an excellent example of convergent evolution but with an unclear genetic basis. We used comparative genomic analyses between Tibetan fishes and fishes from low altitudes and detected genomic signatures of convergent evolution in fishes on the Tibetan Plateau. The Tibetan fishes exhibited genome-wide accelerated evolution in comparison with a control set of fishes from low altitudes. A total of 368 positively selected genes were identified in Tibetan fishes, which were enriched in functional categories related to energy metabolism and hypoxia response. Widespread parallel amino acid substitutions were detected among the Tibetan fishes and a subset of these substitutions occurred in positively selected genes associated with high-altitude adaptation. Functional assays suggested that von Hippel-Lindau (VHL) tumour suppressor genes from Tibetan fishes enhance hypoxia-inducible factor (HIF) activity convergently under hypoxia compared to low-altitude fishes. The results provide genomic and functional evidence supporting convergent genetic mechanisms for high-altitude adaptation in fishes on the Tibetan Plateau.

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.