Whole genome sequencing of Ethiopian highlanders reveals conserved hypoxia tolerance genes

Genomic and physiological mechanisms of high-altitude adaptation in Ethiopian highlanders: a comparative perspective

High-altitude adaptation is a remarkable example of natural selection, yet the genomic and physiological adaptation mechanisms of Ethiopian highlanders remain poorly understood compared to their Andean and Tibetan counterparts. Ethiopian populations, such as the Amhara and Oromo, exhibit unique adaptive strategies characterized by moderate hemoglobin levels and enhanced arterial oxygen saturation, indicating distinct mechanisms of coping with chronic hypoxia. This review synthesizes current genomic insights into Ethiopian high-altitude adaptation, identifying key candidate genes involved in hypoxia tolerance and examining the influence of genetic diversity and historical admixture on adaptive responses. Furthermore, the review highlights significant research gaps, particularly the underrepresentation of Ethiopian populations in global genomic studies, the lack of comprehensive genotype-phenotype analyses, and inconsistencies in research methodologies. Addressing these gaps is crucial for advancing our understanding of the genetic basis of human adaptation to extreme environments and for developing a more complete picture of human physiological resilience. This review offers a comparative perspective with Tibetan and Andean highlanders, emphasizing the need for expanding genomic representation and refining methodologies to uncover the genetic mechanisms underlying high-altitude adaptation in Ethiopian populations.

Physiological and Genetic Basis of High-Altitude Indigenous Animals' Adaptation to Hypoxic Environments

Adaptation is one of the fundamental characteristics of life activities; humans and animals inhabiting high altitudes are well adapted to hypobaric hypoxic environments, and studies on the mechanisms of this adaptation emerged a hundred years ago. Based on these studies, this paper reviews the adaptive changes in hypoxia-sensitive tissues and organs, as well as at the molecular genetic level, such as pulmonary, cardiovascular, O2-consuming tissues, and the hemoglobin and HIF pathway, that occur in animals in response to the challenge of hypobaric hypoxia. High-altitude hypoxia adaptation may be due to the coordinated action of genetic variants in multiple genes and, as a result, adaptive changes in multiple tissues and organs at the physiological and biochemical levels. Unraveling their mechanisms of action can provide a reference for the prevention and treatment of multiple diseases caused by chronic hypoxia.

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.

WGS Data Collections: How Do Genomic Databases Transform Medicine?

As a scientific community we assumed that exome sequencing will elucidate the basis of most heritable diseases. However, it turned out it was not the case; therefore, attention has been increasingly focused on the non-coding sequences that encompass 98% of the genome and may play an important regulatory function. The first WGS-based datasets have already been released including underrepresented populations. Although many databases contain pooled data from several cohorts, recently the importance of local databases has been highlighted. Genomic databases are not only collecting data but may also contribute to better diagnostics and therapies. They may find applications in population studies, rare diseases, oncology, pharmacogenetics, and infectious and inflammatory diseases. Further data may be analysed with Al technologies and in the context of other omics data. To exemplify their utility, we put a highlight on the Polish genome database and its practical application.

Genomic architecture of adaptive radiation and hybridization in Alpine whitefish

Adaptive radiations represent some of the most remarkable explosions of diversification across the tree of life. However, the constraints to rapid diversification and how they are sometimes overcome, particularly the relative roles of genetic architecture and hybridization, remain unclear. Here, we address these questions in the Alpine whitefish radiation, using a whole-genome dataset that includes multiple individuals of each of the 22 species belonging to six ecologically distinct ecomorph classes across several lake-systems. We reveal that repeated ecological and morphological diversification along a common environmental axis is associated with both genome-wide allele frequency shifts and a specific, larger effect, locus, associated with the gene edar. Additionally, we highlight the possible role of introgression between species from different lake-systems in facilitating the evolution and persistence of species with unique trait combinations and ecology. These results highlight the importance of both genome architecture and secondary contact with hybridization in fuelling adaptive radiation.

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.

Genetic differentiation in East African ethnicities and its relationship with endurance running success

Since the 1960s, East African athletes, mainly from Kenya and Ethiopia, have dominated long-distance running events in both the male and female categories. Further demographic studies have shown that two ethnic groups are overrepresented among elite endurance runners in each of these countries: the Kalenjin, from Kenya, and the Oromo, from Ethiopia, raising the possibility that this dominance results from genetic or/and cultural factors. However, looking at the life history of these athletes or at loci previously associated with endurance athletic performance, no compelling explanation has emerged. Here, we used a population approach to identify peaks of genetic differentiation for these two ethnicities and compared the list of genes close to these regions with a list, manually curated by us, of genes that have been associated with traits possibly relevant to endurance running in GWAS studies, and found a significant enrichment in both populations (Kalenjin, P = 0.048, and Oromo, P = 1.6x10^-5). Those traits are mainly related to anthropometry, circulatory and respiratory systems, energy metabolism, and calcium homeostasis. Our results reinforce the notion that endurance running is a systemic activity with a complex genetic architecture, and indicate new candidate genes for future studies. Finally, we argue that a deterministic relationship between genetics and sports must be avoided, as it is both scientifically incorrect and prone to reinforcing population (racial) stereotyping.

Heterozygous Tropomodulin 3 mice have improved lung vascularization after chronic hypoxia

The molecular mechanisms leading to high altitude pulmonary hypertension (HAPH) remains poorly understood. We previously analyzed the whole genome sequence of Kyrgyz highland population and identified eight genomic intervals having a potential role in HAPH. Tropomodulin 3 gene (TMOD3) which encodes a protein that binds and caps the pointed ends of actin filaments and inhibits cell migration, was one of the top candidates. Here we systematically sought additional evidence to validate the functional role of TMOD3. In-silico analysis reveals that some of the SNPs in HAPH associated genomic intervals were positioned in a regulatory region that could result in alternative splicing of TMOD3. In order to functionally validate the role of TMOD3 in HAPH, we exposed Tmod3-/+ mice to 4 weeks of constant hypoxia, i.e. 10% O2 and analyzed both functional (hemodynamic measurements) and structural (angiography) parameters related to HAPH. The hemodynamic measurements, such as right ventricular systolic pressure, a surrogate measure for pulmonary arterial systolic pressure, and right ventricular contractility (RV- ± dP/dt), increases with hypoxia did not separate between Tmod3-/+ and control mice. Remarkably, there was a significant increase in the number of lung vascular branches and total length of pulmonary vascular branches (p < 0.001) in Tmod3-/+ after 4 weeks of constant hypoxia as compared to controls. Notably, the Tmod3-/+ endothelial cells migration was also significantly higher than that from the wild-type littermates. Our results indicate that, under chronic hypoxia, lower levels of Tmod3 play an important role in the maintenance or neo-vascularization of pulmonary arteries.

Multiple mechanisms drive genomic adaptation to extreme O2 levels in Drosophila melanogaster

To detect the genomic mechanisms underlying evolutionary dynamics of adaptation in sexually reproducing organisms, we analyze multigenerational whole genome sequences of Drosophila melanogaster adapting to extreme O₂ conditions over an experiment conducted for nearly two decades. We develop methods to analyze time-series genomics data and predict adaptive mechanisms. Here, we report a remarkable level of synchronicity in both hard and soft selective sweeps in replicate populations as well as the arrival of favorable de novo mutations that constitute a few asynchronized sweeps. We additionally make direct experimental observations of rare recombination events that combine multiple alleles on to a single, better-adapted haplotype. Based on the analyses of the genes in genomic intervals, we provide a deeper insight into the mechanisms of genome adaptation that allow complex organisms to survive harsh environments.

The genomic details of adaptation to extreme environments remain challenging to characterize. Using new methods to analyze flies experimentally evolved to survive extreme O₂ conditions, the authors find a surprising level of synchronicity in selective sweeps, de novo mutations and adaptive recombination events.

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

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

Positive selection in admixed populations from Ethiopia

BACKGROUND

In the process of adaptation of humans to their environment, positive or adaptive selection has played a main role. Positive selection has, however, been under-studied in African populations, despite their diversity and importance for understanding human history.

RESULTS

Here, we have used 119 available whole-genome sequences from five Ethiopian populations (Amhara, Oromo, Somali, Wolayta and Gumuz) to investigate the modes and targets of positive selection in this part of the world. The site frequency spectrum-based test SFselect was applied to idfentify a wide range of events of selection (old and recent), and the haplotype-based statistic integrated haplotype score to detect more recent events, in each case with evaluation of the significance of candidate signals by extensive simulations. Additional insights were provided by considering admixture proportions and functional categories of genes. We identified both individual loci that are likely targets of classic sweeps and groups of genes that may have experienced polygenic adaptation. We found population-specific as well as shared signals of selection, with folate metabolism and the related ultraviolet response and skin pigmentation standing out as a shared pathway, perhaps as a response to the high levels of ultraviolet irradiation, and in addition strong signals in genes such as IFNA, MRC1, immunoglobulins and T-cell receptors which contribute to defend against pathogens.

CONCLUSIONS

Signals of positive selection were detected in Ethiopian populations revealing novel adaptations in East Africa, and abundant targets for functional follow-up.

Cross-Species Insights Into Genomic Adaptations to Hypoxia

Over millions of years, vertebrate species populated vast environments spanning the globe. Among the most challenging habitats encountered were those with limited availability of oxygen, yet many animal and human populations inhabit and perform life cycle functions and/or daily activities in varying degrees of hypoxia today. Of particular interest are species that inhabit high-altitude niches, which experience chronic hypobaric hypoxia throughout their lives. Physiological and molecular aspects of adaptation to hypoxia have long been the focus of high-altitude populations and, within the past decade, genomic information has become increasingly accessible. These data provide an opportunity to search for common genetic signatures of selection across uniquely informative populations and thereby augment our understanding of the mechanisms underlying adaptations to hypoxia. In this review, we synthesize the available genomic findings across hypoxia-tolerant species to provide a comprehensive view of putatively hypoxia-adaptive genes and pathways. In many cases, adaptive signatures across species converge on the same genetic pathways or on genes themselves [i.e., the hypoxia inducible factor (HIF) pathway). However, specific variants thought to underlie function are distinct between species and populations, and, in most cases, the precise functional role of these genomic differences remains unknown. Efforts to standardize these findings and explore relationships between genotype and phenotype will provide important clues into the evolutionary and mechanistic bases of physiological adaptations to environmental hypoxia.

Mitochondrial DNA genomes revealed different patterns of high-altitude adaptation in high-altitude Tajiks compared with Tibetans and Sherpas

High-altitude Tajiks (HA-Tajiks), Tibetans and Sherpas are three groups of high-altitude native people in China. The differences in the mtDNA genome between the three populations and the role of the mtDNA genome in the high-altitude adaptation of HA-Tajiks were seldom investigated. In this study, 80 HA-Tajiks were enrolled, and their whole mtDNA genomes were sequenced. The haplogroup of each subject was determined by comparison to the revised Cambridge Reference Sequence (rCRS). Ten additional populations from East Asia and Central Asia, including Tibetans and Sherpas, were selected as references. The top haplogroup was U, followed by H, T and J. Principle component analysis and genetic distance analysis indicated that HA-Tajiks showed a close relationship with Wakhi Tajiks, Pamiri Tajiks and Sarikoli Tajiks, indicating that they should be considered one nation scattered around the Pamirs. The difference in the mtDNA genome between HA-Tajiks and Sherpas was significantly greater than that between HA-Tajiks and Tibetans. Among the 13 genes related to the OXPHOS pathway encoded by the mtDNA genome, HA-Tajiks showed more significant differences in ND3 and CYTB compared to Tibetans. Compared to Sherpas, HA-Tajiks showed more significant differences in ND1, ND2, COX1, ATP8, ATP6, ND3, ND4L, ND4, ND5 and CYTB. The associated functional changes and underlying molecular mechanisms should be explored by molecular and biochemical investigations in further studies.

Selection signatures for high-altitude adaptation in ruminants

High-altitude areas are important socio-economical habitats with ruminants serving as a major source of food and commodities for humans. Living at high altitude, however, is extremely challenging, predominantly due to the exposure to hypoxic conditions, but also because of cold temperatures and limited feed for livestock. To survive in high-altitude environments over the long term, ruminants have evolved adaptation strategies, e.g. physiological and morphological modifications, which allow them to cope with these harsh conditions. Identification of such selection signatures in ruminants may contribute to more informed breeding decisions, and thus improved productivity. Moreover, studying the genetic background of altitude adaptation in ruminants provides insights into a common molecular basis across species and thus a better understanding of the physiological basis of this adaptation. In this paper, we review the major effects of high altitude on the mammalian body and highlight some of the most important short-term (coping) and genetically evolved (adaptation) physiological modifications. We then discuss the genetic architecture of altitude adaptation and target genes that show evidence of being under selection based on recent studies in various species, with a focus on ruminants. The yak is presented as an interesting native species that has adapted to the high-altitude regions of Tibet. Finally, we conclude with implications and challenges of selection signature studies on altitude adaptation in general. We found that the number of studies on genetic mechanisms that enable altitude adaptation in ruminants is growing, with a strong focus on identifying selection signatures, and hypothesise that the investigation of genetic data from multiple species and regions will contribute greatly to the understanding of the genetic basis of altitude adaptation.

The Genetic Architecture of Chronic Mountain Sickness in Peru

Chronic mountain sickness (CMS) is a pathological condition resulting from chronic exposure to high-altitude hypoxia. While its prevalence is high in native Andeans (>10%), little is known about the genetic architecture of this disease. Here, we performed the largest genome-wide association study (GWAS) of CMS (166 CMS patients and 146 controls living at 4,380 m in Peru) to detect genetic variants associated with CMS. We highlighted four new candidate loci, including the first CMS-associated variant reaching GWAS statistical significance (rs7304081; P = 4.58 × 10⁻⁹). By looking at differentially expressed genes between CMS patients and controls around these four loci, we suggested AEBP2, CAST, and MCTP2 as candidate CMS causal genes. None of the candidate loci were under strong natural selection, consistent with the observation that CMS affects fitness mainly after the reproductive years. Overall, our results reveal new insights on the genetic architecture of CMS and do not provide evidence that CMS-associated variants are linked to a strong ongoing adaptation to 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'.

Novel insight into the genetic basis of high-altitude pulmonary hypertension in Kyrgyz highlanders

The Central Asian Kyrgyz highland population provides a unique opportunity to address genetic diversity and understand the genetic mechanisms underlying high-altitude pulmonary hypertension (HAPH). Although a significant fraction of the population is unaffected, there are susceptible individuals who display HAPH in the absence of any lung, cardiac or hematologic disease. We report herein the analysis of the whole-genome sequencing of healthy individuals compared with HAPH patients and other controls (total n = 33). Genome scans reveal selection signals in various regions, encompassing multiple genes from the first whole-genome sequences focusing on HAPH. We show here evidence of three candidate genes MTMR4, TMOD3 and VCAM1 that are functionally associated with well-known molecular and pathophysiological processes and which likely lead to HAPH in this population. These processes are (a) dysfunctional BMP signaling, (b) disrupted tissue repair processes and (c) abnormal endothelial cell function. Whole-genome sequence of well-characterized patients and controls and using multiple statistical tools uncovered novel candidate genes that belong to pathways central to the pathogenesis of HAPH. These studies on high-altitude human populations are pertinent to the understanding of sea level diseases involving hypoxia as a main element of their pathophysiology.

Mitochondrial dysfunction in iPSC-derived neurons of subjects with chronic mountain sickness

Patients with chronic mountain sickness (CMS) suffer from hypoxemia, erythrocytosis, and numerous neurologic deficits. Here we used induced pluripotent stem cell (iPSC)-derived neurons from both CMS and non-CMS subjects to study CMS neuropathology. Using transmission electron microscopy, we report that CMS neurons have a decreased mitochondrial volume density, length, and less cristae membrane surface area. Real-time PCR confirmed a decreased mitochondrial fusion gene optic atrophy 1 (OPA1) expression. Immunoblot analysis showed an accumulation of the short isoform of OPA1 (S-OPA1) in CMS neurons, which have reduced ATP levels under normoxia and increased lactate dehydrogenase (LDH) release and caspase 3 activation after hypoxia. Improving the balance between the long isoform of OPA1 and S-OPA1 in CMS neurons increased the ATP levels and attenuated LDH release under hypoxia. Our data provide initial evidence for altered mitochondrial morphology and function in CMS neurons, and reveal increased cell death under hypoxia due in part to altered mitochondrial dynamics. NEW & NOTEWORTHY Induced pluripotent stem cell-derived neurons from chronic mountain sickness (CMS) subjects have altered mitochondrial morphology and dynamics, and increased sensitivity to hypoxic stress. Modification of OPA1 can attenuate cell death after hypoxic treatment, providing evidence that altered mitochondrial dynamics play an important role in increased vulnerability under stress in CMS neurons.

Cardiac-specific knockout and pharmacological inhibition of Endothelin receptor type B lead to cardiac resistance to extreme hypoxia

Oxygen plays a central role in cardiac energy metabolism. At high altitude where the ambient oxygen level is low, we found EDNRB is associated with human hypoxia adaptation. Our subsequent study in global heterozygous knockout mice (Ednrb^-/+) revealed that cardiac function was conserved in these mice when exposed to extreme hypoxia. The major goal of this study was (i) to determine the functional role of cardiomyocyte EdnrB in maintaining cardiac function under hypoxic stress and (ii) to validate the phenotypes we detected in Ednrb^-/+ mice using EDNRB blockers. Unlike the global knockouts, cardiac-specific heterozygote (EdnrB^flox/+) and homozygote (EdnrB^flox/flox) EdnrB knockout mice were phenotypically normal. When treated with graded low levels of oxygen (10% and 5% O₂), both EdnrB^flox/+ and EdnrB^flox/flox were hypoxia tolerant. The cardiac indexes at 10% and 5% O₂ for EdnrB^flox/+ were significantly higher and lactate levels were significantly lower when compared to the cre-negative controls (P < 0.05). Simultaneously, mice treated with BQ-788 (EDNRB-specific blocker) had a significantly higher cardiac index (P < 0.005) and significantly lower lactate levels (P < 0.0001) than in control mice. A similar result was obtained with mice treated with Bosentan (non-specific). These data indicate that a lower level or complete lack of EdnrB in the cardiomyocytes significantly improves cardiac performance under extreme hypoxia, a novel role of cardiomyocyte EdnrB in the regulation of cardiac function. Furthermore, this rescue under extreme hypoxia can also be achieved using EDNRB-specific pharmacological agents, e.g., BQ-788. This systematically confirms, both genetically and pharmacologically, the protective role of a lower EDNRB under extreme hypoxia stress.

KEY MESSAGES

Under normal condition, cardiomyocytes-specific EdnrB knockout mice, both heterozygote and homozygote, are phenotypically normal. Under hypoxic condition, a lower level or complete deletion of cardiomyocyte EdnrB conserves cardiac function by maintaining high cardiac index. Similarly, mice treated with both specific (BQ-788) and non-specific (Bosentan) EDNRB blockers are tolerant to hypoxia by maintaining better cardiac function. The oxygen perfusion under extreme hypoxia is better in the mice with lower EDNRB, as depicted by lower lactate level at 5% oxygen. Our current study systematically confirms, both genetically and pharmacologically, the protective role of a lower EDNRB under extreme hypoxia stress. Overall, it supports our hypothesis that studies on human hypoxia adaptation provide new insight to common disease pathogenesis and treatments.

New Insights into the Genetic Basis of Monge's Disease and Adaptation to High-Altitude

Human high-altitude (HA) adaptation or mal-adaptation is explored to understand the physiology, pathophysiology, and molecular mechanisms that underlie long-term exposure to hypoxia. Here, we report the results of an analysis of the largest whole-genome-sequencing of Chronic Mountain Sickness (CMS) and nonCMS individuals, identified candidate genes and functionally validated these candidates in a genetic model system (Drosophila). We used PreCIOSS algorithm that uses Haplotype Allele Frequency score to separate haplotypes carrying the favored allele from the noncarriers and accordingly, prioritize genes associated with the CMS or nonCMS phenotype. Haplotypes in eleven candidate regions, with SNPs mostly in nonexonic regions, were significantly different between CMS and nonCMS subjects. Closer examination of individual genes in these regions revealed the involvement of previously identified candidates (e.g., SENP1) and also unreported ones SGK3, COPS5, PRDM1, and IFT122 in CMS. Remarkably, in addition to genes like SENP1, SGK3, and COPS5 which are HIF-dependent, our study reveals for the first time HIF-independent gene PRDM1, indicating an involvement of wider, nonHIF pathways in HA adaptation. Finally, we observed that down-regulating orthologs of these genes in Drosophila significantly enhanced their hypoxia tolerance. Taken together, the PreCIOSS algorithm, applied on a large number of genomes, identifies the involvement of both new and previously reported genes in selection sweeps, highlighting the involvement of multiple hypoxia response systems. Since the overwhelming majority of SNPs are in nonexonic (and possibly regulatory) regions, we speculate that adaptation to HA necessitates greater genetic flexibility allowing for transcript variability in response to graded levels of hypoxia.

Transcriptome Analysis of Circulating PBMCs to Understand Mechanism of High Altitude Adaptation in Native Cattle of Ladakh Region

Ladakhi cattle is native population of Leh and Ladakh region and constantly exposed to hypobaric hypoxia over many generations. In present study, transcriptome signatures of cattle from Ladakh region (~5500 m) and Sahiwal cattle from tropical regions were evaluated using Agilent 44 K microarray chip. The top up-regulated genes in Ladakhi cows were INHBC, ITPRI, HECA, ABI3, GPR171, and HIF-1α involved in hypoxia and stress response. In Sahiwal cows, the top up-regulated genes eEF1A1, GRO1, CXCL2, DEFB3 and BOLA-DQA3 were associated with immune function and inflammatory response indicating their strong immune potential to combat the pathogens prevalent in the tropical conditions. The molecular pathways highly impacted were MAPK signaling, ETC, apoptosis, TLR signaling and NF- kB signaling pathway indicating signatures of adaptive evolution of these two cattle types in response to diverse environments. Further, qPCR analysis revealed increased expression of DEGs such as HIF-1, EPAS-1, VEGFA, NOS2, and GLUT-1/SLC2A1 in cattle types from high altitude suggesting their pivotal role in association with high altitude adaptation. Based on data generated, native cattle of Ladakh region was found to be genetically distinct from native cattle adapted to the tropical region of India.

Exome Sequencing Provides Evidence of Polygenic Adaptation to a Fat-Rich Animal Diet in Indigenous Siberian Populations

Siberia is one of the coldest environments on Earth and has great seasonal temperature variation. Long-term settlement in northern Siberia undoubtedly required biological adaptation to severe cold stress, dramatic variation in photoperiod, and limited food resources. In addition, recent archeological studies show that humans first occupied Siberia at least 45,000 years ago; yet our understanding of the demographic history of modern indigenous Siberians remains incomplete. In this study, we use whole-exome sequencing data from the Nganasans and Yakuts to infer the evolutionary history of these two indigenous Siberian populations. Recognizing the complexity of the adaptive process, we designed a model-based test to systematically search for signatures of polygenic selection. Our approach accounts for stochasticity in the demographic process and the hitchhiking effect of classic selective sweeps, as well as potential biases resulting from recombination rate and mutation rate heterogeneity. Our demographic inference shows that the Nganasans and Yakuts diverged ∼12,000-13,000 years ago from East-Asian ancestors in a process involving continuous gene flow. Our polygenic selection scan identifies seven candidate gene sets with Siberian-specific signals. Three of these gene sets are related to diet, especially to fat metabolism, consistent with the hypothesis of adaptation to a fat-rich animal diet. Additional testing rejects the effect of hitchhiking and favors a model in which selection yields small allele frequency changes at multiple unlinked genes.

Epigenomics and human adaptation to high altitude

Over the past decade, major technological and analytical advancements have propelled efforts toward identifying the molecular mechanisms that govern human adaptation to high altitude. Despite remarkable progress with respect to the identification of adaptive genomic signals that are strongly associated with the "hypoxia-tolerant" physiological characteristics of high-altitude populations, many questions regarding the fundamental biological processes underlying human adaptation remain unanswered. Vital to address these enduring questions will be determining the role of epigenetic processes, or non-sequence-based features of the genome, that are not only critical for the regulation of transcriptional responses to hypoxia but heritable across generations. This review proposes that epigenomic processes are involved in shaping patterns of adaptation to high altitude by influencing adaptive potential and phenotypic variability under conditions of limited oxygen supply. Improved understanding of the interaction between genetic, epigenetic, and environmental factors holds great promise to provide deeper insight into the mechanisms underlying human adaptive potential, and clarify its implications for biomedical research.

Reduced Cardiac Calcineurin Expression Mimics Long-Term Hypoxia-Induced Heart Defects in Drosophila

BACKGROUND

Hypoxia is often associated with cardiopulmonary diseases, which represent some of the leading causes of mortality worldwide. Long-term hypoxia exposures, whether from disease or environmental condition, can cause cardiomyopathy and lead to heart failure. Indeed, hypoxia-induced heart failure is a hallmark feature of chronic mountain sickness in maladapted populations living at high altitude. In a previously established Drosophila heart model for long-term hypoxia exposure, we found that hypoxia caused heart dysfunction. Calcineurin is known to be critical in cardiac hypertrophy under normoxia, but its role in the heart under hypoxia is poorly understood.

METHODS AND RESULTS

In the present study, we explore the function of calcineurin, a gene candidate we found downregulated in the Drosophila heart after lifetime and multigenerational hypoxia exposure. We examined the roles of 2 homologs of Calcineurin A, CanA14F, and Pp2B in the Drosophila cardiac response to long-term hypoxia. We found that knockdown of these calcineurin catalytic subunits caused cardiac restriction under normoxia that are further aggravated under hypoxia. Conversely, cardiac overexpression of Pp2B under hypoxia was lethal, suggesting that a hypertrophic signal in the presence of insufficient oxygen supply is deleterious.

CONCLUSIONS

Our results suggest a key role for calcineurin in cardiac remodeling during long-term hypoxia with implications for diseases of chronic hypoxia, and it likely contributes to mechanisms underlying these disease states.

High-altitude adaptation in humans: from genomics to integrative physiology

About 1.2 to 33% of high-altitude populations suffer from Monge's disease or chronic mountain sickness (CMS). Number of factors such as age, sex, and population of origin (older, male, Andean) contribute to the percentage reported from a variety of samples. It is estimated that there are around 83 million people who live at altitudes > 2500 m worldwide and are at risk for CMS. In this review, we focus on a human "experiment in nature" in various high-altitude locations in the world-namely, Andean, Tibetan, and Ethiopian populations that have lived under chronic hypoxia conditions for thousands of years. We discuss the adaptive as well as mal-adaptive changes at the genomic and physiological levels. Although different genes seem to be involved in adaptation in the three populations, we can observe convergence at genetic and signaling, as well as physiological levels. What is important is that we and others have shown that lessons learned from the genes mined at high altitude can be helpful in better understanding and treating diseases that occur at sea level. We discuss two such examples: EDNRB and SENP1 and their role in cardiac tolerance and in the polycythemic response, respectively.

Cic Loss Promotes Gliomagenesis via Aberrant Neural Stem Cell Proliferation and Differentiation

Inactivating mutations in the transcriptional repression factor Capicua (CIC) occur in approximately 50% of human oligodendrogliomas, but mechanistic links to pathogenesis are unclear. To address this question, we generated Cic-deficient mice and human oligodendroglioma cell models. Genetic deficiency in mice resulted in a partially penetrant embryonic or perinatal lethal phenotype, with the production of an aberrant proliferative neural population in surviving animals. In vitro cultured neural stem cells derived from Cic conditional knockout mice bypassed an EGF requirement for proliferation and displayed a defect in their potential for oligodendrocyte differentiation. Cic is known to participate in gene suppression that can be relieved by EGFR signal, but we found that cic also activated expression of a broad range of EGFR-independent genes. In an orthotopic mouse model of glioma, we found that Cic loss potentiated the formation and reduced the latency in tumor development. Collectively, our results define an important role for Cic in regulating neural cell proliferation and lineage specification, and suggest mechanistic explanations for how CIC mutations may impact the pathogenesis and therapeutic targeting of oligodendroglioma. Cancer Res; 77(22); 6097-108. ©2017 AACR.

Measuring high-altitude adaptation

High altitudes (>8,000 ft or 2,500 m) provide an experiment of nature for measuring adaptation and the physiological processes involved. Studies conducted over the past ~25 years in Andeans, Tibetans, and, less often, Ethiopians show varied but distinct O₂ transport traits from those of acclimatized newcomers, providing indirect evidence for genetic adaptation to high altitude. Short-term (acclimatization, developmental) and long-term (genetic) responses to high altitude exhibit a temporal gradient such that, although all influence O₂ content, the latter also improve O₂ delivery and metabolism. Much has been learned concerning the underlying physiological processes, but additional studies are needed on the regulation of blood flow and O₂ utilization. Direct evidence of genetic adaptation comes from single-nucleotide polymorphism (SNP)-based genome scans and whole genome sequencing studies that have identified gene regions acted upon by natural selection. Efforts have begun to understand the connections between the two with Andean studies on the genetic factors raising uterine blood flow, fetal growth, and susceptibility to Chronic Mountain Sickness and Tibetan studies on genes serving to lower hemoglobin and pulmonary arterial pressure. Critical for future studies will be the selection of phenotypes with demonstrable effects on reproductive success, the calculation of actual fitness costs, and greater inclusion of women among the subjects being studied. The well-characterized nature of the O₂ transport system, the presence of multiple long-resident populations, and relevance for understanding hypoxic disorders in all persons underscore the importance of understanding how evolutionary adaptation to high altitude has occurred.NEW & NOTEWORTHY Variation in O₂ transport characteristics among Andean, Tibetan, and, when available, Ethiopian high-altitude residents supports the existence of genetic adaptations that improve the distribution of blood flow to vital organs and the efficiency of O₂ utilization. Genome scans and whole genome sequencing studies implicate a broad range of gene regions. Future studies are needed using phenotypes of clear relevance for reproductive success for determining the mechanisms by which naturally selected genes are acting.

Genomic and phenotypic signatures of climate adaptation in an Anolis lizard

Integrated knowledge on phenotype, physiology, and genomic adaptations is required to understand the effects of climate on evolution. The functional genomic basis of organismal adaptation to changes in the abiotic environment, its phenotypic consequences, and its possible convergence across vertebrates are still understudied. In this study, we use a comparative approach to verify predicted gene functions for vertebrate thermal adaptation with observed functions underlying repeated genomic adaptations in response to elevation in the lizard Anolis cybotes. We establish a direct link between recurrently evolved phenotypes and functional genomics of altitude-related climate adaptation in three highland and lowland populations in the Dominican Republic. We show that across vertebrates, genes contained in this interactome are expressed within the brain, the endocrine system, and during development. These results are relevant to elucidate the effect of global climate change across vertebrates and might aid in furthering insight into gene-environment relationships under disturbances to homeostasis.

Physiological Changes to the Cardiovascular System at High Altitude and Its Effects on Cardiovascular Disease

Riley, Callum James, and Matthew Gavin. Physiological changes to the cardiovascular system at high altitude and its effects on cardiovascular disease. High Alt Med Biol. 18:102-113, 2017.-The physiological changes to the cardiovascular system in response to the high altitude environment are well understood. More recently, we have begun to understand how these changes may affect and cause detriment to cardiovascular disease. In addition to this, the increasing availability of altitude simulation has dramatically improved our understanding of the physiology of high altitude. This has allowed further study on the effect of altitude in those with cardiovascular disease in a safe and controlled environment as well as in healthy individuals. Using a thorough PubMed search, this review aims to integrate recent advances in cardiovascular physiology at altitude with previous understanding, as well as its potential implications on cardiovascular disease. Altogether, it was found that the changes at altitude to cardiovascular physiology are profound enough to have a noteworthy effect on many forms of cardiovascular disease. While often asymptomatic, there is some risk in high altitude exposure for individuals with certain cardiovascular diseases. Although controlled research in patients with cardiovascular disease was largely lacking, meaning firm conclusions cannot be drawn, these risks should be a consideration to both the individual and their physician.

Senp1 drives hypoxia-induced polycythemia via GATA1 and Bcl-xL in subjects with Monge's disease

Azad and collaborators propose that Senp1 drives excessive erythropoiesis in high-altitude Andean dwellers suffering from chronic mountain sickness.

In this study, because excessive polycythemia is a predominant trait in some high-altitude dwellers (chronic mountain sickness [CMS] or Monge’s disease) but not others living at the same altitude in the Andes, we took advantage of this human experiment of nature and used a combination of induced pluripotent stem cell technology, genomics, and molecular biology in this unique population to understand the molecular basis for hypoxia-induced excessive polycythemia. As compared with sea-level controls and non-CMS subjects who responded to hypoxia by increasing their RBCs modestly or not at all, respectively, CMS cells increased theirs remarkably (up to 60-fold). Although there was a switch from fetal to adult HgbA0 in all populations and a concomitant shift in oxygen binding, we found that CMS cells matured faster and had a higher efficiency and proliferative potential than non-CMS cells. We also established that SENP1 plays a critical role in the differential erythropoietic response of CMS and non-CMS subjects: we can convert the CMS phenotype into that of non-CMS and vice versa by altering SENP1 levels. We also demonstrated that GATA1 is an essential downstream target of SENP1 and that the differential expression and response of GATA1 and Bcl-xL are a key mechanism underlying CMS pathology.

The African diaspora: history, adaptation and health

The trans-Atlantic slave trade brought millions of Africans to the New World. Advances in genomics are providing novel insights into the history and health of Africans and the diasporan populations. Recent examples reviewed here include the unraveling of substantial hunter-gatherer and 'Eurasian' admixtures across sub-Saharan Africa, expanding our understanding of ancestral African genetics; the global ubiquity of mixed ancestry; the revealing of African ancestry in Latin Americans that likely derived from the slave trade; and understanding of the ancestral backgrounds of APOL1 and LPL found to influence kidney disease and lipid levels, respectively, providing specific insights into disease etiology and health disparities.

Adaptive genetic changes related to haemoglobin concentration in native high-altitude Tibetans

NEW FINDINGS

What is the topic of this review? Tibetans have genetic adaptations that are hypothesized to underlie the distinct set of traits they exhibit at altitude. What advances does it highlight? Several adaptive signatures in the same genomic regions have been identified among Tibetan populations resident throughout the Qinghai-Tibetan Plateau. Many highland Tibetans exhibit a haemoglobin concentration within the range expected at sea level, and this trait is associated with putatively adaptive regions harbouring the hypoxia-inducible factor pathway genes EGLN1, EPAS1 and PPARA. Precise functional variants at adaptive loci and relationships to physiological traits, beyond haemoglobin concentration, are currently being examined in this population. Some native Tibetan, Andean and Ethiopian populations have lived at altitudes ranging from 3000 to >4000 m above sea level for hundreds of generations and exhibit distinct combinations of traits at altitude. It was long hypothesized that genetic factors contribute to adaptive differences in these populations, and recent advances in genomics provide evidence that some of the strongest signatures of positive selection in humans are those identified in Tibetans. Many of the top adaptive genomic regions highlighted thus far harbour genes related to hypoxia sensing and response. Putatively adaptive copies of three hypoxia-inducible factor pathway genes, EPAS1, EGLN1 and PPARA, are associated with sea-level range, rather than elevated, haemoglobin concentration observed in many Tibetans at high altitude, and recent studies provide insight into some of the precise adaptive variants, timing of adaptive events and functional roles. While several studies in highland Tibetans have converged on a few hypoxia-inducible factor pathway genes, additional candidates have been reported in independent studies of Tibetans located throughout the Qinghai-Tibetan Plateau. Various aspects of adaptive significance have yet to be identified, integrated, and fully explored. Given the rapid technological advances and interdisciplinary efforts in genomics, physiology and molecular biology, careful examination of Tibetans and comparisons with other distinctively adapted highland populations will provide valuable insight into evolutionary processes and models for both basic and clinical research.

Genetics of human origin and evolution: high-altitude adaptations

High altitude, defined as elevations lying above 2500m sea level, challenges human survival and reproduction. This environment provides a natural experimental design wherein specific populations, Andeans, Ethiopians, and Tibetans, have lived in a chronic hypoxia state for millennia. These human groups have overcome the low ambient oxygen tension of high elevation via unique physiologic and genetic adaptations. Genomic studies have identified several genes that underlie high-altitude adaptive phenotypes, many of which are central components of the Hypoxia Inducible Factor (HIF) pathway. Further study of mechanisms governing the adaptive changes responsible for high-altitude adaptation will contribute to our understanding of the molecular basis of evolutionary change and assist in the functional annotation of the human genome.

Anchored pseudo-de novo assembly of human genomes identifies extensive sequence variation from unmapped sequence reads

The human genome reference (HGR) completion marked the genomics era beginning, yet despite its utility universal application is limited by the small number of individuals used in its development. This is highlighted by the presence of high-quality sequence reads failing to map within the HGR. Sequences failing to map generally represent 2-5 % of total reads, which may harbor regions that would enhance our understanding of population variation, evolution, and disease. Alternatively, complete de novo assemblies can be created, but these effectively ignore the groundwork of the HGR. In an effort to find a middle ground, we developed a bioinformatic pipeline that maps paired-end reads to the HGR as separate single reads, exports unmappable reads, de novo assembles these reads per individual and then combines assemblies into a secondary reference assembly used for comparative analysis. Using 45 diverse 1000 Genomes Project individuals, we identified 351,361 contigs covering 195.5 Mb of sequence unincorporated in GRCh38. 30,879 contigs are represented in multiple individuals with ~40 % showing high sequence complexity. Genomic coordinates were generated for 99.9 %, with 52.5 % exhibiting high-quality mapping scores. Comparative genomic analyses with archaic humans and primates revealed significant sequence alignments and comparisons with model organism RefSeq gene datasets identified novel human genes. If incorporated, these sequences will expand the HGR, but more importantly our data highlight that with this method low coverage (~10-20×) next-generation sequencing can still be used to identify novel unmapped sequences to explore biological functions contributing to human phenotypic variation, disease and functionality for personal genomic medicine.

Linking between genetic structure and geographical distance: Study of the maternal gene pool in the Ethiopian population

Background The correlation between genetics and geographical distance has already been examined through the study of the dispersion of human populations, especially in terms of uniparental genetic markers. Aim The present work characterises, at the level of the mitochondrial DNA (mtDNA), two new samples of Amhara and Oromo populations from Ethiopia to evaluate the possible pattern of distribution for mtDNA variation and to test the hypothesis of the Isolation-by-Distance (IBD) model among African, European and Middle-Eastern populations. Subjects and methods This study analysed 173 individuals belonging to two ethnic groups of Ethiopia, Amhara and Oromo, by assaying HVS-I and HVS-II of mtDNA D-loop and informative coding region SNPs of mtDNA. Results The analysis suggests a relationship between genetic and geographic distances, affirming that the mtDNA pool of Africa, Europe and the Middle East might be coherent with the IBD model. Moreover, the mtDNA gene pools of the Sub-Saharan African and Mediterranean populations were very different. Conclusion In this study the pattern of mtDNA distribution, beginning with the Ethiopian plateau, was tested in the IBD model. It could be affirmed that, on a continent scale, the mtDNA pool of Africa, Europe and the Middle East might fall under the IBD model.

Altitude Adaptation: A Glimpse Through Various Lenses

Simonson, Tatum S. Altitude adaptation: A glimpse through various lenses. High Alt Med Biol 16:125-137, 2015.--Recent availability of genome-wide data from highland populations has enabled the identification of adaptive genomic signals. Some of the genomic signals reported thus far among Tibetan, Andean, and Ethiopian are the same, while others appear unique to each population. These genomic findings parallel observations conveyed by decades of physiological research: different continental populations, resident at high altitude for hundreds of generations, exhibit a distinct composite of traits at altitude. The most commonly reported signatures of selection emanate from genomic segments containing hypoxia-inducible factor (HIF) pathway genes. Corroborative evidence for adaptive significance stems from associations between putatively adaptive gene copies and sea-level ranges of hemoglobin concentration in Tibetan and Amhara Ethiopians, birth weights and metabolic factors in Andeans and Tibetans, maternal uterine artery diameter in Andeans, and protection from chronic mountain sickness in Andean males at altitude. While limited reports provide mechanistic insights thus far, efforts to identify and link precise genetic variants to molecular, physiological, and developmental functions are underway, and progress on the genomics front continues to provide unprecedented movement towards these goals. This combination of multiple perspectives is necessary to maximize our understanding of orchestrated biological and evolutionary processes in native highland populations, which will advance our understanding of both adaptive and non-adaptive responses to hypoxia.

Sports genetics moving forward: lessons learned from medical research

Sports genetics can take advantage of lessons learned from human disease genetics. By righting past mistakes and increasing scientific rigor, we can magnify the breadth and depth of knowledge in the field. We present an outline of challenges facing sports genetics in the light of experiences from medical research. Sports performance is complex, resulting from a combination of a wide variety of different traits and attributes. Improving sports genetics will foremost require analyses based on detailed phenotyping. To find widely valid, reproducible common variants associated with athletic phenotypes, study sample sizes must be dramatically increased. One paradox is that in order to confirm relevance, replications in specific populations must be undertaken. Family studies of athletes may facilitate the discovery of rare variants with large effects on athletic phenotypes. The complexity of the human genome, combined with the complexity of athletic phenotypes, will require additional metadata and biological validation to identify a comprehensive set of genes involved. Analysis of personal genetic and multiomic profiles contribute to our conceptualization of precision medicine; the same will be the case in precision sports science. In the refinement of sports genetics it is essential to evaluate similarities and differences between sexes and among ethnicities. Sports genetics to date have been hampered by small sample sizes and biased methodology, which can lead to erroneous associations and overestimation of effect sizes. Consequently, currently available genetic tests based on these inherently limited data cannot predict athletic performance with any accuracy.

Shared Genetic Signals of Hypoxia Adaptation in Drosophila and in High-Altitude Human Populations

The ability to withstand low oxygen (hypoxia tolerance) is a polygenic and mechanistically conserved trait that has important implications for both human health and evolution. However, little is known about the diversity of genetic mechanisms involved in hypoxia adaptation in evolving populations. We used experimental evolution and whole-genome sequencing in Drosophila melanogaster to investigate the role of natural variation in adaptation to hypoxia. Using a generalized linear mixed model we identified significant allele frequency differences between three independently evolved hypoxia-tolerant populations and normoxic control populations for approximately 3,800 single nucleotide polymorphisms. Around 50% of these variants are clustered in 66 distinct genomic regions. These regions contain genes that are differentially expressed between hypoxia-tolerant and normoxic populations and several of the differentially expressed genes are associated with metabolic processes. Additional genes associated with respiratory and open tracheal system development also show evidence of directional selection. RNAi-mediated knockdown of several candidate genes’ expression significantly enhanced survival in severe hypoxia. Using genomewide single nucleotide polymorphism data from four high-altitude human populations—Sherpas, Tibetans, Ethiopians, and Andeans, we found that several human orthologs of the genes under selection in flies are also likely under positive selection in all four high-altitude human populations. Thus, our results indicate that selection for hypoxia tolerance can act on standing genetic variation in similar genes and pathways present in organisms diverged by hundreds of millions of years.

Intracellular PAF-Acetylhydrolase Type I

Platelet-activating factor (PAF) is a phospholipid mediator whose synthesis and degradation depend on specific sets of enzymes. PAF-acetylhydrolase (PAF-AH) hydrolyzes the acetyl moiety of PAF at its sn-2 position and thereby inactivates it. PAF-AH Ib, originally identified in brain, exists in the cytoplasm of many (probably all) types of mammalian cells and tissues. PAF-AH Ib consists of three subunits (α1, α2, and β), in which the α subunits provide the catalytic activity. The finding that the β subunit is the product of the causative gene for Miller-Dieker lissencephaly led to extensive analyses of PAF-AH Ib subunits in the field of cell biology and neurobiology. More than 20 molecules are known to bind to PAF-AH Ib subunits, and PAF-AH Ib has been implicated in neuronal development, neuronal functions, Alzheimer's disease, bipolar disorder, cancer, spermatogenesis, and tolerance to hypoxia. However, in almost all of these cases, how the catalytic activity is involved and the identity of the most important substrate of this enzyme are unclear. In this chapter, the structure and functions of PAF-AH Ib and its subunit proteins are summarized and their contributions to human diseases are discussed.

Predicting Carriers of Ongoing Selective Sweeps without Knowledge of the Favored Allele

Methods for detecting the genomic signatures of natural selection have been heavily studied, and they have been successful in identifying many selective sweeps. For most of these sweeps, the favored allele remains unknown, making it difficult to distinguish carriers of the sweep from non-carriers. In an ongoing selective sweep, carriers of the favored allele are likely to contain a future most recent common ancestor. Therefore, identifying them may prove useful in predicting the evolutionary trajectory—for example, in contexts involving drug-resistant pathogen strains or cancer subclones. The main contribution of this paper is the development and analysis of a new statistic, the Haplotype Allele Frequency (HAF) score. The HAF score, assigned to individual haplotypes in a sample, naturally captures many of the properties shared by haplotypes carrying a favored allele. We provide a theoretical framework for computing expected HAF scores under different evolutionary scenarios, and we validate the theoretical predictions with simulations. As an application of HAF score computations, we develop an algorithm (PreCIOSS: Predicting Carriers of Ongoing Selective Sweeps) to identify carriers of the favored allele in selective sweeps, and we demonstrate its power on simulations of both hard and soft sweeps, as well as on data from well-known sweeps in human populations.

Methods for detecting the genomic signatures of natural selection have been heavily studied, and they have been successful in identifying genomic regions under positive selection. However, methods that detect positive selective sweeps do not typically identify the favored allele, or even the haplotypes carrying the favored allele. The main contribution of this paper is the development and analysis of a new statistic (the HAF score), assigned to individual haplotypes. Using both theoretical analyses and simulations, we describe how the HAF scores differ for carriers and non-carriers of the favored allele, and how they change dynamically during a selective sweep. We also develop an algorithm, PreCIOSS, for separating carriers and non-carriers. Our tool has broad applicability as carriers of the favored allele are likely to contain a future most recent common ancestor. Therefore, identifying them may prove useful in predicting the evolutionary trajectory—for example, in contexts involving drug-resistant pathogen strains or cancer subclones.