Detecting past and ongoing natural selection among ethnically Tibetan women at high altitude in Nepal

Genetic characterization of Tibetan pigs adapted to high altitude under natural selection based on a large whole-genome dataset

The Qinghai-Tibet Plateau is a valuable genetic resource pool, and the high-altitude adaptation of Tibetan pigs is a classic example of the adaptive evolution of domestic animals. Here, we report the presence of Darwinian positive selection signatures in Tibetan pigs (TBPs) using 348 genome-wide datasets (127 whole-genome sequence datasets (WGSs) and 221 whole-genome single-nucleotide polymorphism (SNP) chip datasets). We characterized a high-confidence list of genetic signatures related response to high-altitude adaptation in Tibetan pigs, including 4,598 candidate SNPs and 131 candidate genes. Functional annotation and enrichment analysis revealed that 131 candidate genes are related to multiple systems and organs in Tibetan pigs. Notably, eight of the top ten novel genes, RALB, NBEA, LIFR, CLEC17A, PRIM2, CDH7, GK5 and FAM83B, were highlighted and associated with improved adaptive heart functions in Tibetan pigs high-altitude adaptation. Moreover, genome-wide association analysis revealed that 29 SNPs were involved in 13 candidate genes associated with at least one adaptive trait. In particular, among the top ten candidate genes, CLEC17A is related to a reduction in hemoglobin (HGB) in Tibetan pigs. Overall, our study provides a robust SNP/gene list involving genetic adaptation for Tibetan pig high-altitude adaptation, and it will be a valuable resource for future Tibetan pig studies.

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.

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.

Regulation of haemoglobin concentration at high altitude

Lowlanders sojourning for more than 1 day at high altitude (HA) experience a reduction in plasma volume (PV) that increases haemoglobin concentration and thus restores arterial oxygen content. If the sojourn extends over weeks, an expansion of total red cell volume (RCV) occurs and contributes to the haemoconcentration. While the reduction in PV was classically attributed to an increased diuretic fluid loss, recent studies support fluid redistribution, rather than loss, as the underlying mechanism. The fluid redistribution is presumably driven by a disappearance of proteins from the circulation and the resulting reduction in oncotic pressure exerted by the plasma, although the fate of the disappearing proteins remains unclear. The RCV expansion is the result of an accelerated erythropoietic activity secondary to enhanced renal erythropoietin release, but a contribution of other mechanisms cannot be excluded. After return from HA, intravascular volumes return to normal values and the normalisation of RCV might involve selective destruction of newly formed erythrocytes, although this explanation has been strongly challenged by recent studies. In contrast to acclimatised lowlanders, native highlanders originating from the Tibetan and the Ethiopian plateaus present with a normal or only mildly elevated haemoglobin concentration. Genetic adaptations blunting the erythropoietic response to HA exposure have been proposed as an explanation for the absence of more pronounced haemoconcentration in these populations, but new evidence also supports a contribution of a larger than expected PV. The functional significance of the relatively low haemoglobin concentration in Tibetan and Ethiopian highlanders is incompletely understood and warrants further investigation.

Human genetics: An extreme fitness landscape

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

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.

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.

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.

The Oxygen Cascade from Atmosphere to Mitochondria as a Tool to Understand the (Mal)adaptation to Hypoxia

Hypoxia is a life-threatening challenge for about 1% of the world population, as well as a contributor to high morbidity and mortality scores in patients affected by various cardiopulmonary, hematological, and circulatory diseases. However, the adaptation to hypoxia represents a failure for a relevant portion of the cases as the pathways of potential adaptation often conflict with well-being and generate diseases that in certain areas of the world still afflict up to one-third of the populations living at altitude. To help understand the mechanisms of adaptation and maladaptation, this review examines the various steps of the oxygen cascade from the atmosphere to the mitochondria distinguishing the patterns related to physiological (i.e., due to altitude) and pathological (i.e., due to a pre-existing disease) hypoxia. The aim is to assess the ability of humans to adapt to hypoxia in a multidisciplinary approach that correlates the function of genes, molecules, and cells with the physiologic and pathological outcomes. We conclude that, in most cases, it is not hypoxia by itself that generates diseases, but rather the attempts to adapt to the hypoxia condition. This underlies the paradigm shift that when adaptation to hypoxia becomes excessive, it translates into maladaptation.

Individual variations and sex differences in hemodynamics and percutaneous arterial oxygen saturation (SpO2) in Tibetan highlanders of Tsarang in the Mustang district of Nepal

Background

Many studies have indicated specific low-hemoglobin (Hb) adaptation to high altitude in the Tibetan population, but studies focusing on physiological variations within this population are limited. This study aimed to investigate the relationships between SpO₂ and related factors, including individual variations and sex differences, to assess the generality of high-altitude adaptation in the Tibetan population of Tsarang.

Methods

The participants were 31 male and 41 female community-dwelling people aged ≥18 years living in Tsarang, in the Mustang district of Nepal. Height, weight, SpO₂, Hb concentration, finger temperature, heart rate, and blood pressure were measured. Lifestyle information was obtained by interview.

Results

Men had significantly higher systolic blood pressure (p = 0.002) and Hb (p < 0.001) than women. There was no significant correlation between SpO₂ and other parameters in men. In women, SpO₂ was negatively correlated with heart rate (p = 0.036), Hb (p = 0.004), and finger temperature (p = 0.037). In multiple regression analysis, a higher SpO₂ was marginally correlated with lower age (β = −0.109, p = 0.086) and higher Hb (β = 0.547, p = 0.053) in men. In women, higher SpO₂ was significantly correlated with lower heart rate (β = −0.045, p = 0.036) and Hb (β = −0.341, p = 0.018). Mean hemoglobin (95% confidence interval) was 13.6 g/dl (13.1–14.0 g/dl), which is lower than that found previously in Andeans and almost equal to that in Japanese lowlanders measured using the same device. In some participants of both sexes, hemoglobin was >17.0 g/dl.

Conclusion

Higher SpO₂ was marginally correlated with younger age and higher Hb in men and with lower heart rate and lower Hb in women. Hemoglobin concentration was similar to that found previously in lowlanders, but higher in some individuals. These results indicate individual variation and sex differences in the hemodynamics of high-altitude adaptation in Tibetan highlanders of Tsarang, as well as low-Hb adaptation to high altitude equal to that of other Tibetans.

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.

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.

How placenta promotes the successful reproduction in high-altitude populations: a transcriptome comparison between adaptation and acclimatization

As the best adapted high altitude population, Tibetans feature a relatively high offspring survival rate. Genome-wide studies have identified hundreds of candidate SNPs related to high altitude adaptation of Tibetans, although most of them have unknown functional relevance. To explore the mechanisms behind successful reproduction at high altitudes, we compared the placental transcriptomes of Tibetans, sea level Hans (SLHan), and Han immigrants (ImHan). Among the three populations, placentas from ImHan showed a hyperactive gene expression pattern. Their increased activation demonstrates a hypoxic stress response similar to sea level individuals experiencing hypoxic conditions. Unlike ImHan, Tibetan placentas were characterized by the significant up-regulation of placenta-specific genes, and the activation of autophagy and the tricarboxylic acid (TCA) cycle. Certain conserved hypoxia response functions, including the antioxidant system and angiogenesis, were activated in both ImHan and Tibetans, but mediated by different genes. The coherence of specific transcriptome features linked to possible genetic contribution was observed in Tibetans. Furthermore, we identified a novel Tibetan-specific EPAS1 isoform with a partial deletion at exon six, which may be involved in the adaption to hypoxia through the EPAS1-centred gene network in the placenta. Overall, our results show that the placenta grants successful pregnancies in Tibetans by strengthening the natural functions of the placenta itself. On the other hand, the placenta of ImHan was in an inhabiting time-dependent acclimatization process representing a common hypoxic stress response pattern.

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

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

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.

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.

Evidence for shared genetic risk factors between lymphangioleiomyomatosis and pulmonary function

Introduction

Lymphangioleiomyomatosis (LAM) is a rare low-grade metastasising disease characterised by cystic lung destruction. The genetic basis of LAM remains incompletely determined, and the disease cell-of-origin is uncertain. We analysed the possibility of a shared genetic basis between LAM and cancer, and LAM and pulmonary function.

Methods

The results of genome-wide association studies of LAM, 17 cancer types and spirometry measures (forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC), FEV1/FVC ratio and peak expiratory flow (PEF)) were analysed for genetic correlations, shared genetic variants and causality. Genomic and transcriptomic data were examined, and immunodetection assays were performed to evaluate pleiotropic genes.

Results

There were no significant overall genetic correlations between LAM and cancer, but LAM correlated negatively with FVC and PEF, and a trend in the same direction was observed for FEV1. 22 shared genetic variants were uncovered between LAM and pulmonary function, while seven shared variants were identified between LAM and cancer. The LAM-pulmonary function shared genetics identified four pleiotropic genes previously recognised in LAM single-cell transcriptomes: ADAM12, BNC2, NR2F2 and SP5. We had previously associated NR2F2 variants with LAM, and we identified its functional partner NR3C1 as another pleotropic factor. NR3C1 expression was confirmed in LAM lung lesions. Another candidate pleiotropic factor, CNTN2, was found more abundant in plasma of LAM patients than that of healthy women.

Conclusions

This study suggests the existence of a common genetic aetiology between LAM and pulmonary function.

Comparison of sequencing data processing pipelines and application to underrepresented African human populations

Background

Population genetic studies of humans make increasing use of high-throughput sequencing in order to capture diversity in an unbiased way. There is an abundance of sequencing technologies, bioinformatic tools and the available genomes are increasing in number. Studies have evaluated and compared some of these technologies and tools, such as the Genome Analysis Toolkit (GATK) and its “Best Practices” bioinformatic pipelines. However, studies often focus on a few genomes of Eurasian origin in order to detect technical issues. We instead surveyed the use of the GATK tools and established a pipeline for processing high coverage full genomes from a diverse set of populations, including Sub-Saharan African groups, in order to reveal challenges from human diversity and stratification.

Results

We surveyed 29 studies using high-throughput sequencing data, and compared their strategies for data pre-processing and variant calling. We found that processing of data is very variable across studies and that the GATK “Best Practices” are seldom followed strictly. We then compared three versions of a GATK pipeline, differing in the inclusion of an indel realignment step and with a modification of the base quality score recalibration step. We applied the pipelines on a diverse set of 28 individuals. We compared the pipelines in terms of count of called variants and overlap of the callsets. We found that the pipelines resulted in similar callsets, in particular after callset filtering. We also ran one of the pipelines on a larger dataset of 179 individuals. We noted that including more individuals at the joint genotyping step resulted in different counts of variants. At the individual level, we observed that the average genome coverage was correlated to the number of variants called.

Conclusions

We conclude that applying the GATK “Best Practices” pipeline, including their recommended reference datasets, to underrepresented populations does not lead to a decrease in the number of called variants compared to alternative pipelines. We recommend to aim for coverage of > 30X if identifying most variants is important, and to work with large sample sizes at the variant calling stage, also for underrepresented individuals and populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04407-x.

Repeatability of adaptive traits among ethnic Tibetan highlanders

OBJECTIVES

Connecting traits to biological pathways and genes relies on stable observations. Researchers typically determine traits once, expecting careful study protocols to yield measurements free of noise. This report examines that expectation with test-retest repeatability analyses for traits used regularly in research on adaptation to high-altitude hypoxia, often in settings without climate control.

METHODS

Two hundred ninety-one ethnic Tibetan women residing from 3500 to 4200 m in Upper Mustang District, Nepal, provided three observations of hemoglobin concentration, percent of oxygen saturation of hemoglobin, and pulse by noninvasive pulse oximetry under conditions designed to minimize environmental noise.

RESULTS

High-intraclass correlation coefficients and low within-subject coefficients of variation reflected consistent measurements. Percent of oxygen saturation had the highest intraclass correlation coefficient and the smallest within-subject coefficient of variability; measurement noise occurred mainly in the lower values. Hemoglobin concentration and pulse presented slightly higher within-subject coefficients of variation; measurement noise occurred across the range of values. The women had performed the same measurements 7 years earlier using the same devices and protocol. The sample means and SD observed across 7 years differed little. Hemoglobin concentration increased substantially after menopause.

CONCLUSIONS

Analyzing repeatability features of traits may improve our interpretation of statistical analyses and detection of variation from measurement or biology. The high levels of measurement repeatability and biological stability support the continued use of these robust traits for investigating human adaptation in this altitude range.

Fetal growth, high altitude, and evolutionary adaptation: a new perspective

Residence at high altitude is consistently associated with low birthweight among placental mammals. This reduction in birthweight influences long-term health trajectories for both the offspring and mother. However, the physiological processes that contribute to fetal growth restriction at altitude are still poorly understood, and thus our ability to safely intervene remains limited. One approach to identify the factors that mitigate altitude-dependent fetal growth restriction is to study populations that are protected from fetal growth restriction through evolutionary adaptations (e.g., high altitude-adapted populations). Here, we examine human gestational physiology at high altitude from a novel evolutionary perspective that focuses on patterns of physiological plasticity, allowing us to identify 1) the contribution of specific physiological systems to fetal growth restriction and 2) the mechanisms that confer protection in highland-adapted populations. Using this perspective, our review highlights two general findings: first, that the beneficial value of plasticity in maternal physiology is often dependent on factors more proximate to the fetus; and second, that our ability to understand the contributions of these proximate factors is currently limited by thin data from altitude-adapted populations. Expanding the comparative scope of studies on gestational physiology at high altitude and integrating studies of both maternal and fetal physiology are needed to clarify the mechanisms by which physiological responses to altitude contribute to fetal growth outcomes. The relevance of these questions to clinical, agricultural, and basic research combined with the breadth of the unknown highlight gestational physiology at high altitude as an exciting niche for continued work.

Enhanced Placental Mitochondrial Respiration in Tibetan Women at High Altitude

Living at high altitudes is extremely challenging as it entails exposure to hypoxia, low temperatures, and high levels of UV radiation. However, the Tibetan population has adapted to such conditions on both a physiological and genetic level over 30,000–40,000 years. It has long been speculated that fetal growth restriction is caused by abnormal placental development. We previously demonstrated that placentas from high-altitude Tibetans were protected from oxidative stress induced by labor compared to those of European descent. However, little is known about how placental mitochondria change during high-altitude adaptation. In this study, we aimed to uncover the mechanism of such adaptation by studying the respiratory function of the placental mitochondria of high-altitude Tibetans, lower-altitude Tibetans, and lower-altitude Chinese Han. We discovered that mitochondrial respiration was greater in high-altitude than in lower-altitude Tibetans in terms of OXPHOS via complexes I and I+II, ETS_max capacity, and non-phosphorylating respiration, whereas non-ETS respiration, LEAK/ETS, and OXPHOS via complex IV did not differ. Respiration in lower-altitude Tibetans and Han was similar for all tested respiratory states. Placentas from high-altitude Tibetan women were protected from acute ischemic/hypoxic insult induced by labor, and increased mitochondrial respiration may represent an acute response that induces mitochondrial adaptations.

Gene expression plasticity and desert adaptation in house mice

Understanding how organisms adapt to new environments is a key problem in evolution, yet it remains unclear whether phenotypic plasticity generally facilitates or hinders this process. Here we studied evolved and plastic responses to water-stress in lab-born descendants of wild house mice (Mus musculus domesticus) collected from desert and non-desert environments and measured gene expression and organismal phenotypes under control and water-stressed conditions. After many generations in the lab, desert mice consumed significantly less water than mice from other localities, indicating that this difference has a genetic basis. Under water-stress, desert mice maintained more weight than non-desert mice, and exhibited differences in blood chemistry related to osmoregulatory function. Gene expression in the kidney revealed evolved differences between mice from different environments as well as plastic responses between hydrated and dehydrated mice. Desert mice showed reduced expression plasticity under water-stress compared to non-desert mice. Importantly, non-desert mice under water-stress generally showed shifts toward desert-like expression, consistent with adaptive plasticity. Finally, we identify several co-expression modules linked to phenotypes of interest. These findings provide evidence for local adaptation after a recent invasion and suggest that adaptive plasticity may have facilitated colonization of the desert environment.

HIF-1 directly induces TET3 expression to enhance 5-hmC density and induce erythroid gene expression in hypoxia

In mammalian cells, cytosines found within cytosine guanine dinucleotides can be methylated to 5-methylcytosine (5-mC) by DNA methyltransferases and further oxidized by the Ten-eleven translocation dioxygenase (TET) enzymes to 5-hydroxymethylcytosine (5-hmC). We have previously shown that hematopoietic stem and progenitor cells (HSPCs) with TET2 mutations have aberrant 5-hmC distribution and less erythroid differentiation potential. However, these experiments were performed under standard tissue culture conditions with 21% oxygen (O2), whereas HSPCs in human bone marrow reside in ∼1% O2. Therefore, to model human erythropoiesis more accurately, we compared 5-hmC distribution and gene expression in hypoxic vs normoxic conditions. Despite TET enzymes having limited O2 as a substrate in hypoxia, 5-hmC peaks were more numerous and pronounced than in normoxia. Among the TET genes, TET3 was upregulated specifically in hypoxia. We identified 2 HIF-1 binding sites in TET3 by chromatin immunoprecipitation of HIF-1α followed by sequencing, and TET3 upregulation was abrogated with deletion of both sites, indicating that TET3 is a direct HIF-1 target. Finally, we showed that loss of one or both of these HIF-1 binding sites in K562 cells disrupted erythroid differentiation in hypoxia and lowered cell viability. This work provides a molecular link between O2 availability, epigenetic modification of chromatin, and erythroid differentiation.

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.

Human genetic admixture

Throughout human history, large-scale migrations have facilitated the formation of populations with ancestry from multiple previously separated populations. This process leads to subsequent shuffling of genetic ancestry through recombination, producing variation in ancestry between populations, among individuals in a population, and along the genome within an individual. Recent methodological and empirical developments have elucidated the genomic signatures of this admixture process, bringing previously understudied admixed populations to the forefront of population and medical genetics. Under this theme, we present a collection of recent PLOS Genetics publications that exemplify recent progress in human genetic admixture studies, and we discuss potential areas for future work.

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.

Unique hemoglobin dynamics in female Tibetan highlanders

BACKGROUND

Tibetan highlanders have adapted to hypoxic environments through the development of unique mechanisms that suppress an increase in hemoglobin (Hb) concentration even in high-altitude areas. Hb concentrations generally decrease with increasing age. However, in the highlands, chronic altitude sickness is known to occur in the elderly population. To investigate how aging in a hypoxic environment affects Hb levels in Tibetan highlanders, we focused on the Mustang people, who live above 3500 m. We tried to clarify the pure relationship between aging and Hb levels in a hypoxic environment.

RESULTS

We found that the Hb concentration increased with increasing age in females but not in males. Multivariate analysis showed that age, pulse pressure, the poverty index, and vascular diameter were strongly correlated with the Hb concentration.

CONCLUSIONS

We found unique Hb dynamics among female Tibetan highlanders. As seen in these Hb dynamics, there may be sex-based differences in the adaptive mechanism in Tibetan highlanders.

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

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

Population History and Altitude-Related Adaptation in the Sherpa

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

The overlooked significance of plasma volume for successful adaptation to high altitude in Sherpa and Andean natives

In contrast to Andean natives, high-altitude Tibetans present with a lower hemoglobin concentration that correlates with reproductive success and exercise capacity. Decades of physiological and genomic research have assumed that the lower hemoglobin concentration in Himalayan natives results from a blunted erythropoietic response to hypoxia (i.e., no increase in total hemoglobin mass). In contrast, herein we test the hypothesis that the lower hemoglobin concentration is the result of greater plasma volume, rather than an absence of increased hemoglobin production. We assessed hemoglobin mass, plasma volume and blood volume in lowlanders at sea level, lowlanders acclimatized to high altitude, Himalayan Sherpa, and Andean Quechua, and explored the functional relevance of volumetric hematological measures to exercise capacity. Hemoglobin mass was highest in Andeans, but also was elevated in Sherpa compared with lowlanders. Sherpa demonstrated a larger plasma volume than Andeans, resulting in a comparable total blood volume at a lower hemoglobin concentration. Hemoglobin mass was positively related to exercise capacity in lowlanders at sea level and in Sherpa at high altitude, but not in Andean natives. Collectively, our findings demonstrate a unique adaptation in Sherpa that reorientates attention away from hemoglobin concentration and toward a paradigm where hemoglobin mass and plasma volume may represent phenotypes with adaptive significance at high altitude.

Convergent evolution in human and domesticate adaptation to high-altitude environments

Humans and their domestic animals have lived and thrived in high-altitude environments worldwide for thousands of years. These populations have developed a number of adaptations to survive in a hypoxic environment, and several genomic studies have been conducted to identify the genes that drive these adaptations. Here, we discuss the various adaptations and genetic variants that have been identified as adaptive in human and domestic animal populations and the ways in which convergent evolution has occurred as these populations have adapted to high-altitude environments. We found that human and domesticate populations have adapted to hypoxic environments in similar ways. Specific genes and biological pathways have been involved in high-altitude adaptation for multiple populations, although the specific variants differ between populations. Additionally, we found that the gene EPAS1 is often a target of selection in hypoxic environments and has been involved in multiple adaptive introgression events. High-altitude environments exert strong selective pressures, and human and animal populations have evolved in convergent ways to cope with a chronic lack of oxygen.

This article is part of the theme issue ‘Convergent evolution in the genomics era: new insights and directions'.

Evidence of Polygenic Adaptation to High Altitude from Tibetan and Sherpa Genomes

Although Tibetans and Sherpa present several physiological adjustments evolved to cope with selective pressures imposed by the high-altitude environment, especially hypobaric hypoxia, few selective sweeps at a limited number of hypoxia related genes were confirmed by multiple genomic studies. Nevertheless, variants at these loci were found to be associated only with downregulation of the erythropoietic cascade, which represents an indirect aspect of the considered adaptive phenotype. Accordingly, the genetic basis of Tibetan/Sherpa adaptive traits remains to be fully elucidated, in part due to limitations of selection scans implemented so far and mostly relying on the hard sweep model.

In order to overcome this issue, we used whole-genome sequence data and several selection statistics as input for gene network analyses aimed at testing for the occurrence of polygenic adaptation in these high-altitude Himalayan populations. Being able to detect also subtle genomic signatures ascribable to weak positive selection at multiple genes of the same functional subnetwork, this approach allowed us to infer adaptive evolution at loci individually showing small effect sizes, but belonging to highly interconnected biological pathways overall involved in angiogenetic processes.

Therefore, these findings pinpointed a series of selective events neglected so far, which likely contributed to the augmented tissue blood perfusion observed in Tibetans and Sherpa, thus uncovering the genetic determinants of a key biological mechanism that underlies their adaptation to high altitude.