Hematological and pulmonary responses to high altitude in Quechuas: a multivariate approach

The genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture

Western South America was one of the worldwide cradles of civilization. The well-known Inca Empire was the tip of the iceberg of an evolutionary process that started 11,000 to 14,000 years ago. Genetic data from 18 Peruvian populations reveal the following: 1) The between-population homogenization of the central southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes, and Amazonia accompanied cultural and socioeconomic interactions revealed by archeology. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower, and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia. 2) The genetic homogenization between the populations of the arid Andes is not only due to migrations during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the Wari Empire (600 to 1,000 years before present). 3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes, a putative enhancer in HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, a noncoding gene related to cardiovascular function) and rs269868-C/Ser1067 in DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes and, in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T and B lymphocytes in viral-host interaction.

The history behind the mosaic of the Americas

Focusing on literature published in 2018-2020, we review inferences about: (i) how ancient DNA is contributing to clarify the peopling of the Americas and the dispersal of its first inhabitants, (ii) how the interplay between environmental diversity and culture has influenced the genetic structure and adaptation of Andean and Amazon populations, (iii) how genetics has contributed to our understanding of the Pre-Columbian Tupi expansion in Eastern South America, (iv) the subcontinental origins and dynamics of Post-Columbian admixture in the Americas, and finally, (v) episodes of adaptive natural selection in the American continent, particularly in the high altitudes of the Andes.

The increase in hemoglobin concentration with altitude varies among human populations

Decreased oxygen availability at high altitude requires physiological adjustments allowing for adequate tissue oxygenation. One such mechanism is a slow increase in the hemoglobin concentration ([Hb]) resulting in elevated [Hb] in high-altitude residents. Diagnosis of anemia at different altitudes requires reference values for [Hb]. Our aim was to establish such values based on published data of residents living at different altitudes by applying meta-analysis and multiple regressions. Results show that [Hb] is increased in all high-altitude residents. However, the magnitude of increase varies among the regions analyzed and among ethnic groups within a region. The highest increase was found in residents of the Andes (1 g/dL/1000 m), but this increment was smaller in all other regions of the world (0.6 g/dL/1000 m). While sufficient data exist for adult males and females showing that sex differences in [Hb] persist with altitude, data for infants, children, and pregnant women are incomplete preventing such analyses. Because WHO reference values were originally based on [Hb] of South American people, we conclude that individual reference values have to be defined for ethnic groups to reliably diagnose anemia and erythrocytosis in high-altitude residents. Future studies need to test their applicability for children of different ages and pregnant women.

CFTR Modulator Use Is Associated with Higher Hemoglobin Levels in Individuals with Cystic Fibrosis

RATIONALE

Understanding how cystic fibrosis transmembrane conductance regulator (CFTR) modulators influence comorbid conditions like anemia is of great interest to the cystic fibrosis community.

OBJECTIVES

To test the hypothesis that CFTR modulators are associated with higher hemoglobin (Hgb) levels.

METHODS

Annualized Hgb and other laboratory, demographic, and anthropometric data were abstracted from the U.S. CF Foundation Patient Registry for adult and pediatric registrants before and after therapy with ivacaftor (IVA) or lumacaftor/ivacaftor (LUM/IVA) between January 2010 and December 2016. Univariate and multivariate linear mixed models were used to examine the effect of IVA on Hgb in patients with G551D-CFTR, and the effect of LUM/IVA on Hgb in F508del-CFTR homozygotes. Linear regression was used to characterize change in mean Hgb over time.

RESULTS

A total of 1,347 registrants (707 males and 640 females) with G551D-CFTR and 12,582 F508del-CFTR homozygotes (6,640 males and 5,942 females) who had never undergone lung transplant and had contemporaneous data regarding Hgb and CFTR modulator use were identified. IVA was associated with average Hgb increases of 0.54 gm/dl (95% confidence interval [CI], 0.39-0.69; P < 0.0001) and 0.18 gm/dl (95% CI, 0.01-0.35; P = 0.037) for males and females, respectively, with G551D-CFTR. LUM/IVA was associated with average Hgb increases of 0.58 gm/dl (95% CI, 0.48-0.68; P < 0.0001) and 0.26 gm/dl (95% CI, 0.20-0.33; P < 0.0001) for male and female F508del-CFTR homozygotes, respectively. In multivariate models, IVA positively affected Hgb in males but not females, and LUM/IVA positively affected Hgb in both sexes.

CONCLUSIONS

IVA and LUM/IVA use are both associated with higher Hgb levels in patients with CF.

Anthropometric and hematological tests to diagnose nutritional deficiencies in schoolchildren of indigenous communities living in the Andean region of Ecuador

ABSTRACT Objective To carry out the anthropometric and biometric-hematological assessments in schoolchildren of the Andean region of Ecuador, in order to improve the diagnosis of nutritional deficiencies. Methods The study has been carried out in the San Juan School (Chimborazo, Ecuador), located at 3,240m of altitude, to 36 children of 5 and 6 years old. Anthropometric analyses (weight, height and body mass index), and hematocrit and hemoglobin concentrations were measured. The hemoglobin measurement was evaluated considering the normal value and the one adapted to the altitude of the area. Results The schoolchildren showed high prevalence of stunting (44%). The values of hematocrit (.=0.001) and hemoglobin (.=0.003) were higher in girls. It should be highlighted that using the normal value of hemoglobin, anemia was not detected. However almost a fifth of the schoolchildren studied were diagnosed with anemia when we applied the correction factors adapted to the altitude. Conclusion The use of correction factors adapted to the altitude is considered essential to do the hematology test in populations that live in high altitude in order to avoid a false diagnosis. Moreover, it is necessary to establish the environmental factors related to the stunted growth of this population of the Andean region.

A decade of rapid change: Biocultural influences on child growth in highland Peru

OBJECTIVES

In the past decade many areas of Peru have been undergoing extreme environmental, economic, and cultural change. In the highland hamlet of Chugurpampa, La Libertad, climate change has ruined harvests and led to frequent periods of migration to the coast in search of livelihood. This biocultural research examines how the changes could be affecting the growth of children who maintain residence in the highlands.

METHODS

Clinical records from the early 2000s were compared to those from the early 2010s. Charts were randomly selected to record anthropometric data, netting a sample of 75 children ages 0-60 months of age. Analysis of covariance was run to compare mean stature, weight, and BMI between cohorts. Percentage of children who fall below the -2 threshold for z-scores for height and weight were compared by age and cohort.

RESULTS

A significant secular trend in growth was found, with children born more recently larger than those born a decade before. The effect is most notable in the first year of life, with the growth advantage attenuated by the age of 3 for height and age 4 for weight. While children were unlikely to be stunted from 0 to 3 years of age, 44% of the later cohort were stunted and 11% were underweight from 4 to 5 years of age.

CONCLUSIONS

Three possible explanations for the rapid shift are entertained: more time spent on the coast during gestation and early childhood, which may attenuate the effect of hypoxia on child growth; dietary change; and increased use of biomedicine.

Surname-inferred Andean ancestry is associated with child stature and limb lengths at high altitude in Peru, but not at sea level

OBJECTIVES

Native Andean ancestry gives partial protection from reduced birthweight at high altitude in the Andes compared with European ancestry. Whether Andean ancestry is also associated with body proportions and greater postnatal body size at altitude is unknown. Therefore, we tested whether a greater proportion of Andean ancestry is associated with stature and body proportions among Peruvian children at high and low altitude.

METHODS

Height, head circumference, head-trunk height, upper and lower limb lengths, and tibia, ulna, hand and foot lengths, were measured in 133 highland and 169 lowland children aged 6 months to 8.5 years. For highland and lowland groups separately, age-sex-adjusted anthropometry z scores were regressed on the number of indigenous parental surnames as a proxy for Andean ancestry, adjusting for potential confounders (maternal age and education, parity, altitude [highlands only]).

RESULTS

Among highland children, greater Andean ancestry was negatively associated with stature and tibia, ulna, and lower limb lengths, independent of negative associations with greater altitude for these measurements. Relationships were strongest for tibia length: each additional Andean surname or 1,000 m increase at altitude among highland children was associated with 0.18 and 0.65 z score decreases in tibia length, respectively. Anthropometry was not significantly associated with ancestry among lowland children.

CONCLUSIONS

Greater Andean ancestry is associated with shorter stature and limb measurements at high but not low altitude. Gene-environment interactions between high altitude and Andean ancestry may exacerbate the trade-off between chest dimensions and stature that was proposed previously, though we could not test this directly.

Bayesian inferences suggest that Amazon Yunga Natives diverged from Andeans less than 5000 ybp: implications for South American prehistory

BACKGROUND

Archaeology reports millenary cultural contacts between Peruvian Coast-Andes and the Amazon Yunga, a rainforest transitional region between Andes and Lower Amazonia. To clarify the relationships between cultural and biological evolution of these populations, in particular between Amazon Yungas and Andeans, we used DNA-sequence data, a model-based Bayesian approach and several statistical validations to infer a set of demographic parameters.

RESULTS

We found that the genetic diversity of the Shimaa (an Amazon Yunga population) is a subset of that of Quechuas from Central-Andes. Using the Isolation-with-Migration population genetics model, we inferred that the Shimaa ancestors were a small subgroup that split less than 5300 years ago (after the development of complex societies) from an ancestral Andean population. After the split, the most plausible scenario compatible with our results is that the ancestors of Shimaas moved toward the Peruvian Amazon Yunga and incorporated the culture and language of some of their neighbors, but not a substantial amount of their genes. We validated our results using Approximate Bayesian Computations, posterior predictive tests and the analysis of pseudo-observed datasets.

CONCLUSIONS

We presented a case study in which model-based Bayesian approaches, combined with necessary statistical validations, shed light into the prehistoric demographic relationship between Andeans and a population from the Amazon Yunga. Our results offer a testable model for the peopling of this large transitional environmental region between the Andes and the Lower Amazonia. However, studies on larger samples and involving more populations of these regions are necessary to confirm if the predominant Andean biological origin of the Shimaas is the rule, and not the exception.

The Andean adaptive toolkit to counteract high altitude maladaptation: genome-wide and phenotypic analysis of the Collas

During their migrations out of Africa, humans successfully colonised and adapted to a wide range of habitats, including extreme high altitude environments, where reduced atmospheric oxygen (hypoxia) imposes a number of physiological challenges. This study evaluates genetic and phenotypic variation in the Colla population living in the Argentinean Andes above 3500 m and compares it to the nearby lowland Wichí group in an attempt to pinpoint evolutionary mechanisms underlying adaptation to high altitude hypoxia. We genotyped 730,525 SNPs in 25 individuals from each population. In genome-wide scans of extended haplotype homozygosity Collas showed the strongest signal around VEGFB, which plays an essential role in the ischemic heart, and ELTD1, another gene crucial for heart development and prevention of cardiac hypertrophy. Moreover, pathway enrichment analysis showed an overrepresentation of pathways associated with cardiac morphology. Taken together, these findings suggest that Colla highlanders may have evolved a toolkit of adaptative mechanisms resulting in cardiac reinforcement, most likely to counteract the adverse effects of the permanently increased haematocrit and associated shear forces that characterise the Andean response to hypoxia. Regulation of cerebral vascular flow also appears to be part of the adaptive response in Collas. These findings are not only relevant to understand the evolution of hypoxia protection in high altitude populations but may also suggest new avenues for medical research into conditions where hypoxia constitutes a detrimental factor.

A divergent Artiodactyl MYADM-like repeat is associated with erythrocyte traits and weight of lamb weaned in domestic sheep

A genome-wide association study (GWAS) was performed to investigate seven red blood cell (RBC) phenotypes in over 500 domestic sheep (Ovis aries) from three breeds (Columbia, Polypay, and Rambouillet). A single nucleotide polymorphism (SNP) showed genome-wide significant association with increased mean corpuscular hemoglobin concentration (MCHC, P = 6.2×10⁻¹⁴) and genome-wide suggestive association with decreased mean corpuscular volume (MCV, P = 2.5×10⁻⁶). The ovine HapMap project found the same genomic region and the same peak SNP has been under extreme historical selective pressure, demonstrating the importance of this region for survival, reproduction, and/or artificially selected traits. We observed a large (>50 kb) variant haplotype sequence containing a full-length divergent artiodactyl MYADM-like repeat in strong linkage disequilibrium with the associated SNP. MYADM gene family members play roles in membrane organization and formation in myeloid cells. However, to our knowledge, no member of the MYADM gene family has been identified in development of morphologically variant RBCs. The specific RBC differences may be indicative of alterations in morphology. Additionally, erythrocytes with altered morphological structure often exhibit increased structural fragility, leading to increased RBC turnover and energy expenditure. The divergent artiodactyl MYADM-like repeat was also associated with increased ewe lifetime kilograms of lamb weaned (P = 2×10⁻⁴). This suggests selection for normal RBCs might increase lamb weights, although further validation is required before implementation in marker-assisted selection. These results provide clues to explain the strong selection on the artiodactyl MYADM-like repeat locus in sheep, and suggest MYADM family members may be important for RBC morphology in other mammals.

Body size, composition, and blood pressure of high-altitude Quechua from the Peruvian Central Andes (Huancavelica, 3,680 m)

Although much information is available about the effects of high altitude on physiological characteristics, less is know about its effect on body composition. In the present study, anthropometric and body composition variables were investigated in a sample of 77 adult Quechua males from the Peruvian Central Andes (Huancavelica, 3,680 m). The subjects are shorter in relation to body weight than other ethnic groups, whereas body proportions are macrocormic (indicating a long trunk relative to the lower extremities), with intermediate values of the acromial-iliac index. All skinfold thicknesses are low (approximately 15th percentiles of NHANES reference values for the triceps and subscapular skinfolds), but tend to be higher than in the other Quechua populations. Similar results are obtained when percentage fat is estimated. Somatotypes are dominant in mesomorphy with very low ectomorphy. Comparison with a sample of high-altitude Kirghiz (3,200 m), previously studied with the same methods, shows higher values in the Peruvian sample for all variables related to adiposity. The presence of low adiposity in the Quechua population could be associated with stresses of the high-altitude environment. Mean values of blood pressure are very low and there is no correlation with age.

Lung volume, chest size, and hematological variation in low-, medium-, and high-altitude central Asian populations

To evaluate adaptive responses to high-altitude environment, we examined three groups of healthy adult males from Central Asia: 94 high-altitude (HA) Kirghiz subjects (3,200 m above sea level); 114 middle-altitude (MA) Kazakh subjects (2,100 m), and 90 low-altitude (LA) Kirghiz subjects (900 m). Data on chest size (chest perimeter and chest diameter), lung volume (forced expiratory volume (FEV) and forced expiratory volume in 1 sec (FEV1)), and hematological parameters (hemoglobin, erythrocytes, hematocrit, and SaO(2)) are discussed. The results show that 1) chest shape is less flat in the samples living at higher altitude. In the HA sample, chest perimeter is lower but chest excursion is high. 2) In the highlanders, forced vital capacity (FVC) and FEV1 are no higher than in the other samples, even when corrected for stature and body weight. The negative correlation between FVC-FEV1 and age decreases with increasing altitude. 3) The HA and MA samples have higher values of hemoglobin, erythrocytes, and hematocrit. The HA sample has lower SaO(2) and higher arterial oxygen content than the LA sample. No association between hematocrit and age was detected in the four samples. The results indicate that the high-altitude Kirghiz present features of developmental acclimatization to hypobaric hypoxia which are also strongly influenced by other major high-altitude environmental stresses.

Genetic structure of Quechua-speakers of the Central Andes and geographic patterns of gene frequencies in South Amerindian populations

A sample of 141 Quechua-speaking individuals of the population of Tayacaja, in the Peruvian Central Andes, was typed for the following 16 genetic systems: ABO, Rh, MNSs, P, Duffy, AcP1, EsD, GLOI, PGM1, AK, 6-PGD, Hp, Gc, Pi, C3, and Bf. The genetic structure of the population was analyzed in relation to the allele frequencies available for other South Amerindian populations, using a combination of multivariate and multivariable techniques. Spatial autocorrelation analysis was performed independently for 13 alleles to identify patterns of gene flow in South America as a whole and in more specific geographic regions. We found a longitudinal cline for the AcP1*a and EsD*1 alleles which we interpreted as the result of an ancient longitudinal expansion of a putative ancestral population of modern Amerindians. Monmonnier's algorithm, used to identify areas of sharp genetic discontinuity, suggested a clear east-west differentiation of native South American populations, which was confirmed by analysis of the distribution of genetic distances. We suggest that this pattern of genetic structures is the consequence of the independent peopling of western and eastern South America or to low levels of gene flow between these regions, related to different environmental and demographic histories.