Andean and Tibetan Patterns of Adaptation to High Altitude

Objectives: High-altitude hypoxia, or decreased oxygen levels caused by low barometric pressure, challenges the ability of humans to live and reproduce. Despite these challenges, human populations have lived on the Andean Altiplano and the Tibetan Plateau for millennia and exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. We and others have identified natural selection candidate genes and gene regions for these adaptations using dense genome scan data. One gene previously known to be important in cellular oxygen sensing, egl nine homolog 1 (EGLN1), shows evidence of positive selection in both Tibetans and Andeans. Interestingly, the pattern of variation for this gene differs between the two populations. Continued research among Tibetan populations has identified statistical associations between hemoglobin concentration and single nucleotide polymorphism (SNP) genotype at EGLN1 and a second gene, endothelial PAS domain protein 1 (EPAS1). Methods: To measure for the effects of EGLN1 and EPAS1 altitude genotypes on hemoglobin concentration among Andean highlanders, we performed a multiple linear regression analysis of 10 candidate SNPs in or near these two genes. Results: Our analysis did not identify significant associations between EPAS1 or EGLN1 SNP genotypes and hemoglobin concentration in Andeans. Conclusions: These results contribute to our understanding of the unique set of adaptations developed in different highland groups to the hypoxia of high altitude. Overall, the results provide key insights into the patterns of genetic adaptation to high altitude in Andean and Tibetan populations

Possible Positive Selection for an Arsenic-Protective Haplotype in Humans

BACKGROUND: Arsenic in drinking water causes severe health effects. Indigenous people in the South American Andes have likely lived with arsenic-contaminated drinking water for thousands of years. Inhabitants of San Antonio de los Cobres (SAC) in the Argentinean highlands generally carry an AS3MT (the major arsenic-metabolizing gene) haplotype associated with reduced health risks due to rapid arsenic excretion and lower urinary fraction of the monomethylated metabolite. OBJECTIVES: We hypothesized an adaptation to high-arsenic living conditions via a possible positive selection for protective AS3MT variants and compared AS3MT haplotype frequencies among different indigenous groups. METHODS: Indigenous groups we evaluated were a) inhabitants of SAC and villages near Salta in northern Argentina (n = 346), b) three Native American populations from the Human Genome Diversity Project (HGDP; n = 25), and c) five Peruvian populations (n = 97). The last two groups have presumably lower historical exposure to arsenic. RESULTS: We found a significantly higher frequency of the protective AS3MT haplotype in the SAC population (68.7%) compared with the HGDP (14.3%, p < 0.001, Fisher exact test) and Peruvian (50.5%, p < 0.001) populations. Genome-wide micro-satellite (n = 671) analysis showed no detectable level of population structure between SAC and Peruvian populations (measure of population differentiation F-ST = 0.006) and low levels of structure between SAC and HGDP populations (F-ST < 0.055 for all pairs of populations compared). CONCLUSIONS: Because population stratification seems unlikely to explain the differences in AS3MT haplotype frequencies, our data raise the possibility that, during a few thousand years, natural selection for tolerance to the environmental stressor arsenic may have increased the frequency of protective variants of AS3MT. Further studies are needed to investigate this hypothesis

Neurovascular Coupling Remains Intact During Incremental Ascent to High Altitude (4240 m) in Acclimatized Healthy Volunteers

Neurovascular coupling (NVC) is the temporal link between neuronal metabolic activity and regional cerebral blood flow (CBF), supporting adequate delivery of nutrients. Exposure to high altitude (HA) imposes several stressors, including hypoxia and hypocapnia, which modulate cerebrovascular tone in an antagonistic fashion. Whether these contrasting stressors and subsequent adaptations affect NVC during incremental ascent to HA is unclear. The aim of this study was to assess whether incremental ascent to HA influences the NVC response. Given that CBF is sensitive to changes in arterial blood gasses, in particular PaCO2, we hypothesized that the vasoconstrictive effect of hypocapnia during ascent would decrease the NVC response. 10 healthy study participants (21.7 ± 1.3 years, 23.57 ± 2.00 kg/m2, mean ± SD) were recruited as part of a research expedition to HA in the Nepal Himalaya. Resting posterior cerebral artery velocity (PCAv), arterial blood gasses (PaO2, SaO2, PaCO2, [HCO3-], base excess and arterial blood pH) and NVC response of the PCA were measured at four pre-determined locations: Calgary/Kathmandu (1045/1400 m, control), Namche (3440 m), Deboche (3820 m) and Pheriche (4240 m). PCAv was measured using transcranial Doppler ultrasound. Arterial blood draws were taken from the radial artery and analyzed using a portable blood gas/electrolyte analyzer. NVC was determined in response to visual stimulation (VS; Strobe light; 6 Hz; 30 s on/off × 3 trials). The NVC response was averaged across three VS trials at each location. PaO2, SaO2, and PaCO2 were each significantly decreased at 3440, 3820, and 4240 m. No significant differences were found for pH at HA (P &gt; 0.05) due to significant reductions in [HCO3-] (P &lt; 0.043). As expected, incremental ascent to HA induced a state of hypoxic hypocapnia, whereas normal arterial pH was maintained due to renal compensation. NVC was quantified as the delta (Δ) PCAv from baseline for mean PCAv, peak PCAv and total area under the curve (ΔPCAv tAUC) during VS. No significant differences were found for Δmean, Δpeak or ΔPCAv tAUC between locations (P &gt; 0.05). NVC remains remarkably intact during incremental ascent to HA in healthy acclimatized individuals. Despite the array of superimposed stressors associated with ascent to HA, CBF and NVC regulation may be preserved coincident with arterial pH maintenance during acclimatization

Cardiorespiratory hysteresis during incremental high altitude ascent-descent quantifies the magnitude of ventilatory acclimatization

Maintenance of arterial blood gases is achieved through sophisticated regulation of ventilation, mediated by central and peripheral chemoreflexes. Respiratory chemoreflexes are important during exposure to high altitude due to the competing influence of hypoxia and hypoxic hyperventilation‐mediated hypocapnia on steady‐state ventilatory drive. Inter‐individual variability exists in ventilatory acclimatization to high altitude, potentially affecting the development of acute mountain sickness (AMS). We aimed to quantify ventilatory acclimatization to high altitude by comparing differential ascent and descent values (i.e. hysteresis) in steady‐state cardiorespiratory variables. We hypothesized that (a) the hysteresis area formed by cardiorespiratory variables during ascent and descent would quantify the magnitude of ventilatory acclimatization, and (b) larger hysteresis areas would be associated with lower AMS symptom scores during ascent. In 25 healthy, Diamox‐free trekkers ascending to and descending from 5160 m, cardiorespiratory hysteresis was measured in the pressure of end‐tidal (PET)CO2, peripheral oxygen saturation (SpO2), minute ventilation (V̇E), chemoreceptor stimulus index (SI; PETCO2/SpO2) and the calculated steady‐state chemoreflex drive (SS‐CD; V̇E/SI) using portable devices (capnograph, peripheral pulse oximeter and respirometer, respectively). AMS symptoms were assessed daily using the Lake Louise Questionnaire. We found that (a) ascent‐descent hysteresis was present in all cardiorespiratory variables, (b) SS‐CD is a valid metric for tracking ventilatory acclimatization to high altitude and (c) highest AMS scores during ascent were significantly, moderately and inversely‐correlated to SS‐CD hysteresis magnitude (rs = ‐0.408, P = 0.043). We propose that ascent‐descent hysteresis is a novel and feasible way to quantify ventilatory acclimatization in trekkers during high altitude exposure

Identifying Signatures of Natural Selection in Tibetan and Andean Populations Using Dense Genome Scan Data

High-altitude hypoxia (reduced inspired oxygen tension due to decreased barometric pressure) exerts severe physiological stress on the human body. Two high-altitude regions where humans have lived for millennia are the Andean Altiplano and the Tibetan Plateau. Populations living in these regions exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. Although these responses have been well characterized physiologically, their underlying genetic basis remains unknown. We performed a genome scan to identify genes showing evidence of adaptation to hypoxia. We looked across each chromosome to identify genomic regions with previously unknown function with respect to altitude phenotypes. In addition, groups of genes functioning in oxygen metabolism and sensing were examined to test the hypothesis that particular pathways have been involved in genetic adaptation to altitude. Applying four population genetic statistics commonly used for detecting signatures of natural selection, we identified selection-nominated candidate genes and gene regions in these two populations (Andeans and Tibetans) separately. The Tibetan and Andean patterns of genetic adaptation are largely distinct from one another, with both populations showing evidence of positive natural selection in different genes or gene regions. Interestingly, one gene previously known to be important in cellular oxygen sensing, EGLN1 (also known as PHD2), shows evidence of positive selection in both Tibetans and Andeans. However, the pattern of variation for this gene differs between the two populations. Our results indicate that several key HIF-regulatory and targeted genes are responsible for adaptation to high altitude in Andeans and Tibetans, and several different chromosomal regions are implicated in the putative response to selection. These data suggest a genetic role in high-altitude adaption and provide a basis for future genotype/phenotype association studies necessary to confirm the role of selection-nominated candidate genes and gene regions in adaptation to altitude

Development of a Panel of Genome-Wide Ancestry Informative Markers to Study Admixture Throughout the Americas

Most individuals throughout the Americas are admixed descendants of Native American, European, and African ancestors. Complex historical factors have resulted in varying proportions of ancestral contributions between individuals within and among ethnic groups. We developed a panel of 446 ancestry informative markers (AIMs) optimized to estimate ancestral proportions in individuals and populations throughout Latin America. We used genome-wide data from 953 individuals from diverse African, European, and Native American populations to select AIMs optimized for each of the three main continental populations that form the basis of modern Latin American populations. We selected markers on the basis of locus-specific branch length to be informative, well distributed throughout the genome, capable of being genotyped on widely available commercial platforms, and applicable throughout the Americas by minimizing within-continent heterogeneity. We then validated the panel in samples from four admixed populations by comparing ancestry estimates based on the AIMs panel to estimates based on genome-wide association study (GWAS) data. The panel provided balanced discriminatory power among the three ancestral populations and accurate estimates of individual ancestry proportions (R2>0.9 for ancestral components with significant between-subject variance). Finally, we genotyped samples from 18 populations from Latin America using the AIMs panel and estimated variability in ancestry within and between these populations. This panel and its reference genotype information will be useful resources to explore population history of admixture in Latin America and to correct for the potential effects of population stratification in admixed samples in the region

Hydrologic controls on basin-scale distribution of benthic invertebrates

Streamflow variability is a major determinant of basin-scale distributions of benthic invertebrates. Here we present a novel procedure based on a probabilistic approach aiming at a spatially explicit quantitative assessment of benthic invertebrate abundance as derived from near-bed flow variability. Although the proposed approach neglects ecological determinants other than hydraulic ones, it is nevertheless relevant in view of its implications on the predictability of basin-scale patterns of organisms. In the present context, aquatic invertebrates are considered, given that they are widely employed as sensitive indicators of fluvial ecosystem health and human-induced perturbations. Moving from the analytical characterization of site-specific probability distribution functions of streamflow and bottom shear stress, we achieve a spatial extension to an entire stream network. Bottom shear stress distributions, coupled with habitat suitability curves derived from field studies, are used to produce maps of invertebrate suitability to shear stress conditions. Therefore, the proposed framework allows one to inspect the possible impacts on river ecology of human-induced perturbations of streamflow variability. We apply this framework to an Austrian river network for which rainfall and streamflow time series, river network hydraulic properties, and local information on invertebrate abundance for a limited number of sites are available. A comparison between observed species density versus modeled suitability to shear stress is also presented. Although the proposed strategy focuses on a single controlling factor and thus represents an ecological minimal model, it allows derivation of important implications for water resource management and fluvial ecosystem protection. Key Points Hydrologic variability is a major control of invertebrate habitat suitability New analytical basin-scale approach for pdfs of ecohydrological key features Austrian river basin used for ecohydrological data-model compariso