Positive selection in admixed populations from Ethiopia
BMC Genetics volume 21, Article number: 108 (2020)
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.
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.
Signals of positive selection were detected in Ethiopian populations revealing novel adaptations in East Africa, and abundant targets for functional follow-up.
Genetic and archaeological data demonstrate that Africa is the origin of anatomically modern humans [1,2,3,4], and that populations outside Africa derive from an Out-of-Africa (OoA) migration some 60,000 years ago [5,6,7,8,9]. African populations are genetically more diverse than any other human population, holding the highest amount of genetic variation, low linkage disequilibrium (LD), and deep population structure [10,11,12,13]. They also carry high cultural and phenotypic diversity, speak almost one-third of the world’s languages , live in a wide variety of environments including deserts, tropical rainforests and mountain highlands, and follow many subsistence strategies, including pastoralism, agriculture and hunter-gathering . However, African populations are underrepresented in big genetic projects such as the HGDP , 1000 Genomes Project  and HapMap . Consequently, not only is our understanding of the evolutionary processes that shape human diversity and adaptation limited, but also medical studies are prone to falter when African populations are included, due either to the fact that the single nucleotide polymorphisms (SNP) used are ascertained mainly in Eurasian populations, or to the lower LD found in all African populations . Additional African-specific studies are needed to counterbalance this historical bias [13, 19].
Ethiopian populations lie geographically near a possible embarkation point of the OoA migration [20, 21], exhibit high linguistic diversity encompassing three branches of the Afroasiatic language family (Omotic, Semitic, Cushitic) and also the Nilotic language family, and inhabit environments from lowland to highland. Previous genotyping studies have found a strong match between linguistic and genetic structures, and revealed admixture between Ethiopian (principally Afroasiatic) and OoA populations (most likely from West Asia) around 2600 years ago, contributing about half of the ancestry of some present-day populations in what has been called a “back to Africa” migration [11, 22, 23]. A 4500-year old ancient Ethiopian fossil, Mota, does not show this West Asian backflow , and provides direct insights into the earlier genetic make-up.
Because African populations have adapted to a variety of environments and subsistence strategies, it is crucial to conduct natural selection studies in order to observe how selective pressures shaped their genomes and understand both our evolutionary history as a species and the population-specific local adaptations to these circumstances. Given the diverse features of African populations, we could expect to find a considerable number of signals of local adaptation. A number of approaches have been established to detect positive selection [25,26,27], and a few signals of adaptive selection in Africans have been reported. Some of the most well-known cases involve malaria resistance, driven by genes such as Glucose-6-phosphate dehydrogenase (G6PD) and the Duffy antigen protein . There is also evidence of high-altitude adaptation in Ethiopians living in the highlands, as well of recent positive selection for lactase persistence in eastern African pastoralists [9, 29,30,31]. However, genome-wide analyses of adaptive selection footprints have often reported fewer signals in Africans than in OoA populations [32,33,34], or failed to find adaptive selection in Africans, some arguing that neutral simulations demonstrate that the tails of the empirical distributions contain mainly false positive signals , meaning that demographic events (bottlenecks, population structure and expansions), rather than selection, dominate the results . Thus African populations offer a challenge in recognizing events of adaptive selection in the genome.
In addition, the power and false-positive rates of positive selection tests in recently-admixed populations have only been addressed in a few studies. In a study of African-Americans using real and simulated genetic data, recent admixture did not result in an increase of false positive rates for site frequency spectrum-based tests, but in general the power decreased . In contrast, in some cases when the selective pressure was very strong, studying the admixed population could provide more power to detect selection than the ancestral population because the signature of derived alleles around the fixed selected site was lost in the ancestral population, but admixture made them polymorphic again producing a signature that is easier to detect. Studies with Latin American , Tibetan , Malagasy  and South Asian  populations have found potential admixture-mediated adaptive regions using this methodology, although some controversy exists since another Latin American study did not find evidence of directional selection after admixture .
These examples reveal that detecting positive selection is far from trivial. Most positive selection tests assume a simple model of a hard sweep, where a mutation arises and spreads rapidly in a population until fixation, carrying the adjacent neutral variation with it . However, the relative importance of hard versus soft sweeps in explaining the adaptation of different human populations is debated, and all forms of selection need to be considered [17, 44,45,46,47].
Here, we analyse previously-generated whole-genome sequences of 119 Ethiopians from five different populations  covering a wide geographical range and belonging to four different linguistic groups (Nilotic, Omotic, Cushitic and Semitic) (Table 1 and Supplementary Fig. 1 in Additional File 1). We provide new information about the adaptive processes that these populations have undergone by first detecting the regions of the genome that have been selected, and then interpreting the biological meaning and context of these adaptations.
SFselect selection analysis
In order to detect selective sweeps, we first analysed the data using the site frequency spectrum-based test SFselect  to identify old events of selection in the five Ethiopian populations (Supplementary Fig. 2, Additional File 1). This approach generates a score for each 30 kilobase (kb) window in the genome. We assessed the statistical significance of the scores by defining a critical value of the test, after performing extensive neutral simulations (see methods), as corresponding to the 99.99th percentile of the neutral distribution (see Methods); the threshold is different for each population. Our simulations were based on a three-population demographic model representing Africans, Europeans and Asians , adding an admixture event between Africans and Europeans. We calculated two different thresholds, one for an unadmixed African population (here, the Gumuz) and the second for an admixed African population (here, the four Afroasiatic populations) (Supplementary Fig. 3 in Additional File 1 and Supplementary Table 1 in Additional file 2). After applying the relevant threshold to each of the five populations, we obtained windows considered as putative candidates for adaptive selection. The number of significant windows is shown in Table 2a.
To interpret these windows, we annotated the protein-coding genes that intersected them (Supplementary Table 2, Additional file 2). Many of the signals were shared between Afroasiatic populations (Table 2a and Fig. 1a), as expected from their genetic similarity and shared environment (see Additional File 3 and Supplementary Figs. 4 and 5 in Additional File 1 for a short demographic analysis of the studied populations). The Amhara and Oromo populations shared the highest number of signals (79), whereas the Gumuz shared the least (from 37 to 41). We found many examples of shared signals of selection by all five East African populations (Supplementary Table 3 in Additional file 2 and Fig. 1). We discuss illustrative examples of shared and population-specific signals here, and further examples in the Additional File 3.
One of the top-scoring windows in all populations (Amhara 4.8, Oromo 4.7, Somali 5.2, Wolayta 5.1, Gumuz 3.7) contains genes including FOLR1 and FOLR2 (Fig. 2a), members of the folic acid receptor family. Members of this gene family bind folic acid and its reduced derivatives, and transport 5-methyltetrahydrofolate into cells. The gene products are secreted proteins that either anchor to membranes via a glycosyl-phosphatidylinositol linkage or exist in a soluble form. Mutations in these genes have been associated with neurodegeneration due to cerebral folate transport deficiency; supplementation of folic acid is usually recommended for pregnant women to avoid neural tube defects during foetal development . Folate is also essential in DNA synthesis, survival and growth of the malaria parasite, so antifolate antimalarial drugs are widely used in the treatment of malaria . To our knowledge, this is the first study that finds this gene cluster to be under selection. The fact that we found this window under selection in all populations, together with the important functions of these genes especially during development, indicates that these genes have probably played a pivotal role during the evolutionary history of East Africans and possibly in general within the human species. We discovered, and discuss below, other selection signals related directly or indirectly to folic acid metabolism.
Another example of a top-scoring window among all five populations does not directly overlap with any gene, but the very strong signal lies downstream of the gene ZNF473 (Fig. 2b and Supplementary Table 2 in Additional file 2). This is an interesting region since it has been described as under long-term balancing selection in African populations, and that has been recently targeted by positive selection in Eurasian populations .
A signal shared among the Amhara, Somali and Wolayta populations contains MRC1 (Supplementary Table 2a, c, d in Additional file 2 and Fig. 2c). MRC1 (also known as CD206) encodes a mannose receptor that is part of the C-type lectin superfamily and plays important roles in both adaptive and innate immune systems such as clearance of endogenous molecules and antigen presentation . MRC1 is an endocytic receptor that can bind to numerous endogenous and exogenous molecules and is mainly expressed in macrophages, dendritic cells and nonvascular epithelium . Numerous studies have shown that the C-type lectin-like domain (CTLD) of MRC1 can bind to viruses (HIV, Dengue, HBV), fungi (Candida albicans) and bacteria (Mycobacterium tuberculosis) . It has also been shown that MRC1 can internalize antigens that can then be processed for cross-presentation in antigen-presenting cells  and that it directly interacts with and inhibits CD45 on the T-cell surface resulting in impaired cytotoxic activity of T-cells and antigen-specific T-cell tolerance . This inhibitory effect of T cells by MRC1 has been proposed as a possible therapeutic strategy to downregulate the excessive immune response of autoimmune diseases. In fact, variants in MRC1 have been associated with asthma  and sarcoidosis . In addition, variants of MRC1 have been associated with susceptibility to leprosy in Vietnamese and Brazilian patients  and to pulmonary tuberculosis in Chinese patients . This example introduces a second recurring theme, of selection on defence-related genes, which will be encountered further below.
We also detected population-specific signals of positive selection, and a particularly strong signal was found in the Amhara population, where the 30 kb window containing IFNA14, IFNA16 and IFNA17 showed a very high and statistically significant SFselect score of 3.9. These genes are members of the Interferon Alpha gene family (Fig. 2d and Supplementary Table 2a in Additional file 2); Interferon Alpha is produced in virus-infected leukocytes and has antiviral activity. It has been shown in vitro that IFNA17 is three times more efficient against Hepatitis C than IFNA2A, which is the most effective current treatment . Moreover, polymorphisms in IFNA17 have been associated with a 3.6-fold increased risk for Crimean-Congo Haemorrhagic Fever development . These interferon genes provide further examples of selection on likely defence against pathogens. Other categories of population-specific signal are discussed in the Additional File 3.
iHS selection analysis captures recent events of selection
We next used the linkage-disequilibrium-based test iHS  to detect recent events of selection in the five Ethiopian populations (Supplementary Fig. 6 in Additional File 1). We analysed mean iHS scores in windows of 30 kb that passed the critical value defined after performing extensive neutral simulations (see Methods). We set a restrictive threshold at the 99.99th percentile of the neutral distribution and only signals that passed this threshold were considered as candidates for adaptive selection (Supplementary Table 1 in Additional file 2 and Supplementary Fig. 3 in Additional File 1).
The number of significant windows per population was low, and the five populations shared some windows under recent positive selection (Table 2b, Supplementary Table 4 in Additional file 2 and Fig. 3a). Amhara and Oromo shared the highest number of windows (12 each) with Wolayta, while the lowest number of window shared (four each) was between Somali and Wolayta, and Somali and Gumuz. The Gumuz in general shared the lowest number of windows with the rest of the populations. The lower numbers of significant windows and shared windows from the iHS analysis compared with the SFselect analysis could be because the populations split quite recently, so there has been little time for selection signals to build up.
Although the intersection of signals between all five populations is modest, we do find some strong shared signals (Supplementary Table 4 in Additional file 2). OTOA (otoancorin) shows high and significant mean iHS scores in all populations (except Wolayta, which is close) and variants with significant p-values. Specifically, one of the top variants in all populations (rs370153558) show p-values of 1.8 × 10− 6, 5 × 10− 8, 3.4 × 10− 5, 1 × 10− 6, 9 × 10− 8 for Amhara, Oromo, Somali, Wolayta and Gumuz respectively (Fig. 4a). All the variants in the OTOA gene found under strong selection lie in the intron 21. The protein encoded by OTOA is expressed on the apical surface of epithelial cells in the sensory organs of the inner ear. Mutations in OTOA have been found in Palestinians and Pakistanis to be causative for autosomal recessive deafness 22 [65, 66]. Hearing is a rapidly-evolving phenotype in humans  and thus the likely target of selection here.
Several examples of selection signals were related to ultraviolet (UV) protection or skin pigmentation, and thus indirectly to the folic acid metabolism discussed above. One of these was UVRAG, where we find a signal in all Afroasiatic populations. The strongest signal is found in the Oromo and Somali with mean iHS scores of 2.6 and 3.4 in the region and specific intronic variants such as rs10899132 (Somali p = 10− 6, Oromo p = 10− 4) are found although no clear functional predictions are yet described (Supplementary Table 5 in Additional file 2 and Fig. 4b). The Amhara also show significant mean iHS scores in other windows containing UVRAG (mean iHS score of 2.65 and rs7117696 variant with p = 5 × 10− 5). The lack of a signal from SFselect in this region of the genome also supports the idea of recent selection. UVRAG plays an essential role in protecting cells from UV-induced DNA damage by activating the nucleotide excision repair pathway . In addition, it acts as an autophagic tumour suppressor that is mutated in common human cancers .
A second candidate, shared between Wolayta and Gumuz, is BNC2 with mean iHS scores per window 2.79 and 2.86 (Fig. 4c) and specific variants such as rs113571602 with significant p-values (4 × 10− 7 and 2 × 10− 7 respectively). Again, all the highest scoring variants fall in an intron of BNC2, pointing towards a putative regulatory change of gene expression. This gene codes for a DNA-binding zinc-finger protein that acts as an mRNA-processing enzyme and a transcription factor . It is expressed in melanocytes and keratinocytes and variants have been associated with skin colour, where higher expression levels correspond to darker skin . Interestingly, BNC2 has been found to lie in an adaptive introgressed region from Neanderthals to Europeans , but the signal of selection in Ethiopia lies outside the reported introgressed region.
The third example in this category is found in a region containing ZRANB3 with statistically significant mean iHS scores of 3.32 and variants with significant p-values such as rs11892059 (p = 1.8 × 10− 6) (Supplementary Table 5c in Additional file 2, Fig. 4d). ZRANB3 is an annealing helicase, fork remodeller and structure-specific nuclease; its deficiency can cause genome instability and hypersensitivity to diverse DNA damaging agents such as UV radiation . This region has previously been reported as a putative selection candidate in the Maasai population  but the authors did not link the signal to adaptation related to UV radiation because this variant is in linkage disequilibrium with the well-known lactase (LCT) gene which has many times been reported to be under selection in several populations [29, 75]. In Ethiopians, the signal is clearly in ZRANB3 and not in LCT. ZRANB3 has also been found under selection among black Tibetan wild boars, providing more evidence for its important function to maintain genomic stability against the high UV radiation found in the Tibetan Plateau .
Many examples of selection signals related to defence were also found. Among these was a window in the Wolayta showing a statistically significant mean iHS of 3.14 and many variants with p < 10− 5 in the upstream region of IFNL1 and the intergenic location between IFNL1 and IFNL2. Signals in the 99.9th percentile are also found in the Amhara and Gumuz (Supplementary Table 5d in Additional file 2 and Fig. 4e). The Interferon-λ family or type III IFNs has three members (IFNL1, IFNL2, IFNL3). These genes play a critical role in antiviral, antiproliferative, antitumor and immune responses . These responses often overlap with IFN-α functions such as MHC class I antigen expression and induction of antiviral cascades. Some of the antiviral activities of INFLs target hepatitis B and C virus, cytomegalovirus, influenza A virus, coronaviruses, encephalomyocarditis virus, intestinal infection viruses (noroviruses and rotaviruses) and human immunodeficiency virus by inducing the expression of antiviral proteins by the infected cells . Clinical trials against hepatitis C virus have tested PEGylated IFNL1 and showed a better or equal effectiveness than PEGylated IFN-α with less extrahepatic adverse effects . Since humans are very frequently exposed to viruses of low pathogenicity, and IFN-λ mostly targets mucosal epithelial cells, the function of type III IFNs could be to protect from infections without triggering the severe inflammation and tissue damage that type I IFNs often produce in the long term . Hence IFN-λ could have been a good target of local adaptation to newly encountered pathogens. Additional signals related to defence in TPCN1, CHUK, THEMIS and TRAV are discussed in the Additional File 3.
Finally, a signal specific to Gumuz was found in PKD2L1 (Fig. 4f), with high iHS scoring variants such as rs74154621 and rs74154622 (iHS score p < 10− 8) that are both whole blood eQTLs with a normalized effect size of − 0.669 and p < 10− 4 according to the GTEx portal. We also find several non-synonymous changes with a high derived allele frequency in the Gumuz (rs17112895 and rs7909153 both at a frequency of 0.70). PKD2L1 belongs to the TRPP subfamily of ion channels that are characterised by large extracellular domains . It is expressed in a subset of taste receptor cells in specific taste areas in mice  and has been identified as a candidate for sour taste in mammals . In humans, two patients with sour taste ageusia have been reported and neither had detectable PKD2L1 transcripts, indicating a potential role of PKD2L1 in human sour taste . Sour taste is one of the five basic tastes and although other tastes have a clear evolutionary purpose (sweet indicates carbohydrate rich food, salty taste sodium, bitter potentially poisonous and umami protein rich), sour tasting remains unexplored in humans. One of the main hypotheses of the evolutionary sour tasting function is that it could warn against the acidic ingestion of rotten or immature fruit . Further signals of selection from other functional categories such as skin pigmentation (BLOC1S2) and one in an RNA gene (NSUN3), were also detected (Additional File 3).
Effect of admixture on detecting ancient and recent selection
The power and false-positive rates of positive selection tests in admixed populations have only been addressed in a few studies [9, 37,38,39,40]. To provide further support for our selection analyses, we have investigated whether similar results could be obtained without the West Asian ancestry genetic component among the Afroasiatic populations. For that purpose, we masked the West Asian component from our data, keeping only the East African component (see Methods). Given that, on average, almost half of the genome was masked by this procedure, we merged all four Afroasiatic populations in a single meta-population and re-ran the positive selection tests used previously. Principal Component Analysis (PCA) of the retained East African component confirmed the high similarity between the East African component of the Afroasiatic populations, supporting the combined meta-analysis of all the individuals (Supplementary Fig. 7 in Additional File 1).
The comparisons between the top 20 signals of the SFselect analysis between each single population and the merged East African component show a high similarity between the two analyses (Fig. 1b). In contrast, the overlaps of the iHS analyses were not as strong (Fig. 3b). This last result could be due because of the breaking down of the Ethiopian haplotypes by the ancestry switches that occurred after the West Asian admixture in the area or because of the nature of iHS that detects recent selection more likely to be specific to each population.
Enrichment of west Asian ancestry in windows under selection
The masked West Asian component measures the proportion of West Asian ancestry in each population (Table 3). The Amhara and Oromo populations have the highest amount (54 and 51%, respectively), Wolayta and Somali show 43 and 44%, respectively, while in contrast the Gumuz show the low amount of 0.7%. These values agree with previous estimates . To detect for any enrichment of West Asian ancestry in windows under selection, the same calculation of the proportion of the West Asian component was performed among regions under positive selection (99.99th percentile after neutral simulations and 1% extreme scores) for all populations, for both iHS and SFselect. This analysis revealed a general increase of West Asian ancestry among the regions putatively under selection found with both the SFselect and iHS tests, with similar percentages for the two tests (Table 3). A resampling analysis shows that the difference is highly significant (p < 10− 5, see Methods): there is thus an overall enrichment of West Asian ancestry in regions under selection. It is worth mentioning that this enrichment of West Asian ancestry is not a source of false positive signals in our analysis given the results obtained when we analysed the effect of admixture on detecting adaptive selection (see above).
Unbalanced ancestry regions
Previous studies have used ancestral component proportions to detect regions with a strong ancestry imbalance that could potentially have positive or negative effects on the fitness of admixed populations . The admixture event between Ethiopian and West Asian populations is dated to 2500–3000 years ago [7, 11], meaning that under a neutral model, we would expect the percentage of the West Asian ancestry component to be evenly distributed across the genome. Therefore we report regions with significant deviations from the expected distribution of West Asian component in several populations that could be candidates of adaptation.
A good example is a long stretch of chromosome 17 with an extreme 95% of African ancestry spanning more than 0.5 Megabases (Mb) in all Afroasiatic populations. Moreover, in this region we find the CRHR1 (Corticotrophin Releasing Hormone Receptor 1) gene with a high SFselect score of 3 (significant after simulations) in the Amhara, and 2.22 and 2.1 in the Somali and Gumuz, respectively (close to significance). There are other genes in the region such as KANSL1 and MAPT (Fig. 5a). We also find in all Afroasiatic populations, except Wolayta, an excess of African ancestry in windows under selection containing among other genes FADS1 and FADS2, two enzymes that participate in the omega-3 and omega-6 biosynthesis and found to be under positive selection in other human populations (Fig. 5b) [86, 87]. In Oromo and Wolayta, high African ancestry (77 and 81% respectively) and high iHS scores (top scoring SNPs with p < 10− 4) were found around the immunoglobulin heavy variable 1–8 genes (IGHV1–8), central to defence (Fig. 5c).
Signatures of polygenic adaptation through functional enrichment analysis
Functional enrichment analysis can be used to understand the biological functions of groups of genes, in this case those that have putatively been under positive selection. For this analysis, we listed the genes contained in windows with scores higher than the empirical 99.5 percentile, either for SFselect or mean iHS. We relaxed the thresholds of significance since we are trying to detect loci contributing to polygenic selection and a biological term was considered significant if the p-value after a Benjamini-Hochbert correction was below an alpha value of 0.05 (Table 4). Details of the significant terms and associated genes in each population, and selection tests, can be found in the Additional file 3. Many of the biological categories significant in all populations are related to immune responses and defence (Table 4). Folate metabolism is also a recurrent function found in many populations, as well as for calcium homeostasis related functions. Finally, muscle development function also appears in several populations.
All in all, the enrichment analysis reinforced our previous analyses of selection, again highlighting several of the main adaptations that Ethiopian populations have undergone.
In this study, we have found new gene candidates under adaptive selection in populations from East Africa. We have been able, by performing extensive simulations, to assess the significance of our candidate adaptive selection signals. We have provided evidence for both old and recent selective sweeps, and both shared and population-specific signals of selection, while accounting for any effect of admixture. Our work has also highlighted the genetic similarity among Afroasiatic Ethiopian populations since many of the old signals of selection are shared between them.
Selection analysis in recently admixed populations is of special interest as the adaptation process may maintain pre-admixture adaptations or use one of the components as the genetic background for new adaptations. It is thus of interest to compare the more ancient (likely pre-admixture) and more recent (post-admixture and population specific) adaptations. The site-frequency-based test SFselect captures ancient and shared selection events before the gene flow from West Asia into Africa and thus before the admixture of the West Asian and the East African components. Conversely, iHS captures recent selection that probably happened after admixture, and that is often population-specific either for the Nilotic (Gumuz) or Afroasiatic populations (Amhara, Oromo, Wolayta, Somali).
We have found that folate metabolism appears to have been crucial for Ethiopian populations, a trait that is new as an adaptation . Specifically, we have identified the genes FOLR1, FOLR2 and DHFRL1 (see Additional File 3) as candidates of adaptive selection, while the functional enrichment analysis also highlighted folate metabolism as a main function potentially under selection, and many genes related to skin pigmentation or UV protection were picked out. Folate is crucial for DNA biosynthesis, methylation and repair and its deficiency can cause fatal birth defects and hence can directly affect reproductive success. Sufficient folate is associated with a 72% reduced risk of neural tube defects  and it is known that folate deficiency severely challenges the nucleotide excision repair mechanism needed to remove UV induced DNA photoproducts . Ethiopia experiences very high ultraviolet radiation, which has consequences that include severe DNA damage and impaired genome integrity. It has been hypothesized that under high UVB and UVA radiation, dark skin pigmentation has been selected in order to avoid folate photolysis (the “vitamin D-folate hypothesis”) [90, 91]. Among the recent selective sweeps were many on genes involved in UV radiation response and pigmentation. In the Afroasiatic populations, we have found a region containing the UVRAG gene that activates the nucleotide excision repair pathway when there is UV-induced damage in cells. We have also found as selection candidates BNC2 (among Amhara, Oromo, Wolayta and Gumuz) whose high expression is associated with dark skin colour, and BLOC1S2, encoding a subunit of the complex BLOC-1 that produces strong pigmentation phenotypes in mice and Hermansky-Pudlak syndrome in humans and also many functional enrichment categories related to UV responses. In addition, we have found ZRANB3 in the Somali population, where deficiency causes genome instability and hypersensitivity to DNA damaging including UV radiation and has been found to be under selection in black boars from the Tibetan plateau (an area also challenged with a high UV radiation) . Thus, there is strong evidence in our study pointing towards folate and pigmentation related adaptations.
The environmental changes and migrations that humans have often experienced have made immunological adaptations a key process during human evolution. Our study gives further insights into these immune-related adaptations in East Africa where the major causes of death are due to infections (HIV, tuberculosis, malaria and other acute lower respiratory infections). For example, we have found in Amhara a region containing IFNA genes that encode for interferon alpha (pivotal for antiviral responses) and a region that the Amhara, Somali and Wolayta share in common containing MRC1 (an endocytic receptor involved in adaptive and innate immune responses). Most importantly, we have found in the Gumuz population regions under potentially recent adaptive selection containing genes belonging to the immunoglobulin heavy constant and variable chains and to the T-cell receptor alpha variable locus.
We used two positive selection methods to detect signals of adaptation in the genomes of five Ethiopian populations. Although we have been able to highlight potentially adaptive regions through computational methods and elaborated on the possible biological implications that could have been pivotal for adaptation in East Africa, this is just a first step towards a better understanding of human adaptive evolution. Our work provides the foundation for further functional studies that are necessary to fully understand the adaptive phenotype behind the footprints of positive selection in Ethiopian genomes.
The dataset comprised five East African populations (Amhara, Oromo, Somali, Wolayta and Gumuz) with 24 individuals each from Pagani et al. 2015. One Wolayta individual was excluded from all subsequent analysis due to a high degree of relatedness (data not shown). Additional samples from the 1000 Genomes Project  and a set of 100 Egyptian samples also from Pagani et al. 2015 were included in PCA and ADMIXTURE analyses. The genome assembly of the data is GRCh37 (hg19). A summary of the dataset is shown in Table 1.
PCA and ADMIXTURE
The Principal Component Analysis was performed with smartpca from the Eigensoft 6.0.1 software . All individuals from the dataset were used to perform the worldwide PCA. For the local PCA, all Ethiopians were included, plus a random subset of 24 West African Yoruba (YRI) and 24 European (CEU) individuals from 1000 Genomes Project . We applied a general filter requiring minor allele frequency higher than 0.05. The PCA of the West Asian masked samples was done with the lsqproject mode that is suitable when the samples have large amounts of missing data.
Population structure analysis was performed with the ADMIXTURE software  on a reduced set of 13 populations, 24 individuals per population (with the exception of 23 Wolayta). Variants were pruned using the PLINK software  with parameters --indep 50 5 2 to remove the effect of linkage disequilibrium.
SFselect and iHS
SFselect is a machine-learning site frequency spectrum-based method to detect adaptive selection in polymorphism data , which is available at . The program was developed using supervised learning (support vector machines) trained with extensive forward population simulations. The simulations were performed under a neutral scenario and under a positive selection scenario where a selected allele experienced 200 different combinations of the parameters s (selection coefficient) and τ (time under selection). SFselect shows high power to detect positive selection compared to other tests based on the site frequency spectrum. The minimum sample size of 46 chromosomes per population provides with enough accuracy to make adaptive selection inferences based on the site frequency spectrum . In this study, we used the general support vector machine trained model of SFselect and applied the test by dividing the whole genome into 30 kb windows with 5 kb overlap between windows.
We used the linkage disequilibrium-based test iHS  to detect recent events of selection in the five Ethiopian populations. The sample size per population is of minimum 46 chromosomes, which according to  provides with enough power to detect signals of positive selection with iHS (minimum of 40 chromosomes is recommended). We used the physical positions to calculate iHS since there is no specific genetic map for these populations. We used the software rehh 2.0  to calculate iHS for all the variants with a minor allele frequency higher than 0.05 and excluded a variant from the calculation if a 20 kb gap was found when calculating EHHs, as they may produce biases. In addition to the iHS score per SNP, we also calculated the mean iHS score (average of iHS scores across SNPs), and the maximum iHS value and –log10(p-value) of a SNP in each 30 kb window; these windows were the same as in the SFselect analysis.
To annotate the protein-coding gene content of windows, we used bedtools 2.24.0  to intersect windows with the hg19 gene annotations from RefSeq. To annotate individual variants, ANNOVAR  was used.
Masking was performed as described previously . African and West Asian ancestries of the Ethiopian individuals were deconvoluted using PCAdmix on 20-SNP windows. The CEU and Gumuz populations were used as surrogate sources for the West Asian component and East African component respectively. The West Asian ancestry was subsequently masked.
After the masking procedure, the proportion of West Asian ancestry in a population was estimated by averaging the proportion of masked data across each SNP. For a specific 30 kb window, the same calculation was done but only including SNPs falling in the window. Consecutive 30 kb windows under selection were merged when calculating the West Asian component proportions. A resampling analysis was used to test if the general increase of West Asian component ancestry among the significant 99.99th percentile SFselect and iHS windows was significant. We sampled the number of selected windows 105 times from the genome-wide windows and calculated the mean West Asian component ancestry in each to obtain a distribution of means. The values obtained for windows were compared with this distribution.
To test whether demographic events could mimic the genomic patterns expected from adaptation, we performed extensive simulations using a simple demographic model that captures the key elements to define the critical values for each of the tests. We used the sequence simulator SLiM  to generate samples of the human neutral demography. A demographic model adapted from  was used, adding a simple model of admixture between a sub-Saharan population and an OoA like population  (as a proxy for the West Asians) 2600 years ago (Supplementary Fig. 8 in Additional File 1). For simplicity, the Amhara population was used as example to model the admixture event common to all Afroasiatic populations, using a West Asian admixture proportion of 0.54.
We next checked the validity of the model by comparing the derived site frequency spectra from the real and simulated data (Supplementary Fig. 9 in Additional File 1). The main differences between real and simulated data were seen among the singletons: a deficit of singletons was observed in the real data due to the low coverage, but otherwise the differences are very small, meaning that our model fits our data well.
There is an increase of extreme SFselect and iHS scores in our real data (Supplementary Fig. 3 in Additional File 1). The 99.99th percentile SFselect score thresholds after the neutral simulations for the Gumuz and the Afroasiatic populations are 2.24 and 2.31, respectively (Supplementary Table 1 in Additional file 2). For iHS, we calculated after the neutral simulations the 99.99th percentile of both the per SNP p-value distribution and the 30-kb window of the mean absolute iHS scores (for an easy comparison with SFselect). We found that for the SNP-based analysis, the 99.99th percentiles per SNP were 3.88 and 3.62 for Gumuz and the Afroasiatic populations, respectively. The window analysis set the 99.99th percentile thresholds at 2.54 and 2.53.
Functional enrichment analysis
To understand the biological functions that may have been under positive selection, we used ClueGo , a Cytoscape  plug-in that integrates Gene Ontology, KEGG pathways Reactome and Wikipathway databases to map groups of genes to specific functions. ClueGo enables visualisation in a functionally grouped annotation network, a pie graph showing the group leading terms (most significant term among a group) and a histogram with all significant terms after p-adjustment (< 0.05 after Benjamini-Hochbert correction) and their number of genes from the analysed cluster found in our list of genes. In this case, we used the genes falling among the top 99.5 percentile of SFselect and mean iHS scores. All information about significant terms and associated genes for each population and selection tests can be found in the Additional file 3.
Availability of data and materials
The Ethiopian datasets analysed during the current study are available in the European Genome-phenome Archive (EGA) repository with accession number EGAS00001000238.
Utah residents with Northern and Western European Ancestry
Out of Africa
Principal Component Analysis
Single nucleotide polymorphism
Yoruba from Ibadan, Nigeria
Hublin J-J, Ben-Ncer A, Bailey SE, Freidline SE, Neubauer S, Skinner MM, et al. New fossils from Jebel Irhoud, Morocco and the pan-African origin of Homo sapiens. Nature. 2017;546:289–92.
Ramachandran S, Deshpande O, Roseman CC, Rosenberg NA, Feldman MW, Cavalli-Sforza LL. Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa. Proc Natl Acad Sci. 2005;102:15942–7.
Nielsen R, Akey JM, Jakobsson M, Pritchard JK, Tishkoff S, Willerslev E. Tracing the peopling of the world through genomics. Nature. 2017;541:302–10.
Schlebusch CM, Jakobsson M. Tales of Human Migration, Admixture, and Selection in Africa. Annu Rev Genomics Hum Genet. 2018;19:annurev-genom-083117-021759.
Campbell MC, Tishkoff SA. The evolution of human genetic and phenotypic variation in Africa. Curr Biol. 2010;20:R166–73.
Gravel S, Henn BM, Gutenkunst RN, Indap AR, Marth GT, Clark AG, et al. Demographic history and rare allele sharing among human populations. Proc Natl Acad Sci. 2011;108:11983–8.
Pagani L, Schiffels S, Gurdasani D, Danecek P, Scally A, Chen Y, et al. Tracing the route of modern humans out of Africa by using 225 human genome sequences from Ethiopians and Egyptians. Am J Hum Genet. 2015;96:986–91.
Schiffels S, Durbin R. Inferring human population size and separation history from multiple genome sequences. Nat Genet. 2014;46:919–25.
Huerta-Sánchez E, DeGiorgio M, Pagani L, Tarekegn A, Ekong R, Antao T, et al. Genetic signatures reveal high-altitude adaptation in a set of Ethiopian populations. Mol Biol Evol. 2013;30:1877–88.
Kreager P, Winney B, Ulijaszek S, Capelli C. Population in the human sciences: concepts, models, evidence: Oxford University Press; 2015.
Pagani L, Kivisild T, Tarekegn A, Ekong R, Plaster C, Gallego Romero I, et al. Ethiopian genetic diversity reveals linguistic stratification and complex influences on the Ethiopian gene pool. Am J Hum Genet. 2012;91:83–96.
Dobon B, Hassan HY, Laayouni H, Luisi P, Ricaño-Ponce I, Zhernakova A, et al. The genetics of East African populations: a Nilo-Saharan component in the African genetic landscape. Sci Rep. 2015;5:9996.
Gurdasani D, Carstensen T, Tekola-Ayele F, Pagani L, Tachmazidou I, Hatzikotoulas K, et al. The African genome variation project shapes medical genetics in Africa. Nature. 2015;517:327–32.
Eberhard D, Simons G. Fenning C. Languages of the World: Ethnologue; 2020. http://www.ethnologue.com.
Campbell MC, Tishkoff SA. African genetic diversity: implications for human demographic history, modern human origins, and complex disease mapping. Annu Rev Genomics Hum Genet. 2008;9:403–33.
Cann HM, de Toma C, Cazes L, Legrand M-F, Morel V, Piouffre L, et al. A human genome diversity cell line panel. Science. 2002;296:261–2.
Gibbs RA, Boerwinkle E, Doddapaneni H, Han Y, Korchina V, Kovar C, et al. A global reference for human genetic variation. Nature. 2015;526:68–74.
Gibbs RA, Belmont JW, Hardenbol P, Willis TD, Yu F, Zhang H, et al. The international HapMap project. Nature. 2003;426:789–96.
Sherman RM, Forman J, Antonescu V, Puiu D, Daya M, Rafaels N, et al. Assembly of a pan-genome from deep sequencing of 910 humans of African descent. Nat Genet. 2019;51:30–5.
Melé M, Javed A, Pybus M, Zalloua P, Haber M, Comas D, et al. Recombination gives a new insight in the effective population size and the history of the Old World human populations. Mol Biol Evol. 2012;29:25–30.
Quintana-Murci L, Semino O, Bandelt H-J, Passarino G, McElreavey K, Santachiara-Benerecetti AS. Genetic evidence of an early exit of Homo sapiens sapiens from Africa through eastern Africa. Nat Genet. 1999;23:437–41.
Olivieri A, Achilli A, Pala M, Battaglia V, Fornarino S, Al-Zahery N, et al. The mtDNA legacy of the Levantine early upper Palaeolithic in Africa. Science. 2006;314:1767–70.
Pickrell JK, Patterson N, Loh P-R, Lipson M, Berger B, Stoneking M, et al. Ancient west Eurasian ancestry in southern and eastern Africa. Proc Natl Acad Sci U S A. 2014;111:2632–7.
Gallego Llorente M, Jones ER, Eriksson A, Siska V, Arthur KW, Arthur JW, et al. Ancient Ethiopian genome reveals extensive Eurasian admixture throughout the African continent. Science (80- ). 2015;350:820–2.
Fan S, Hansen MEB, Lo Y, Tishkoff SA. Going global by adapting local: a review of recent human adaptation. Science. 2016;354:54–9.
Wilde S, Timpson A, Kirsanow K, Kaiser E, Kayser M, Unterländer M, et al. Direct evidence for positive selection of skin, hair, and eye pigmentation in Europeans during the last 5,000 y. Proc Natl Acad Sci. 2014;111:4832–7.
Mathieson I, Lazaridis I, Rohland N, Mallick S, Patterson N, Roodenberg SA, et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature. 2015;528:499–503.
Kwiatkowski DP. How malaria has affected the human genome and what human genetics can teach us about malaria. Am J Hum Genet. 2005;77:171–92.
Tishkoff SA, Reed FA, Ranciaro A, Voight BF, Babbitt CC, Silverman JS, et al. Convergent adaptation of human lactase persistence in Africa and Europe. Nat Genet. 2007;39:31–40.
Scheinfeldt LB, Soi S, Thompson S, Ranciaro A, Woldemeskel D, Beggs W, et al. Genetic adaptation to high altitude in the Ethiopian highlands. Genome Biol. 2012;13:R1.
Udpa N, Ronen R, Zhou D, Liang J, Stobdan T, Appenzeller O, et al. Whole genome sequencing of Ethiopian highlanders reveals conserved hypoxia tolerance genes. Genome Biol. 2014;15:R36.
Storz JF, Payseur BA, Nachman MW. Genome scans of DNA variability in humans reveal evidence for selective sweeps outside of Africa. Mol Biol Evol. 2004;21:1800–11.
Coop G, Pickrell JK, Novembre J, Kudaravalli S, Li J, Absher D, et al. The role of geography in human adaptation. PLoS Genet. 2009;5:e1000500.
Pybus M, Luisi P, Dall’Olio GM, Uzkudun M, Laayouni H, Bertranpetit J, et al. Hierarchical boosting: a machine-learning framework to detect and classify hard selective sweeps in human populations. Bioinformatics. 2015;31:btv493.
Granka JM, Henn BM, Gignoux CR, Kidd JM, Bustamante CD, Feldman MW. Limited evidence for classic selective sweeps in African populations. Genetics. 2012;192:1049–64.
Hofer T, Ray N, Wegmann D, Excoffier L. Large allele frequency differences between human continental groups are more likely to have occurred by drift during range expansions than by selection. Ann Hum Genet. 2009;73:95–108.
Lohmueller KE, Bustamante CD, Clark AG. Detecting directional selection in the presence of recent admixture in African-Americans. Genetics. 2011;187:823–35.
Jin W, Xu S, Wang H, Yu Y, Shen Y, Wu B, et al. Genome-wide detection of natural selection in African Americans pre- and post-admixture. Genome Res. 2012;22:519–27.
Jeong C, Alkorta-Aranburu G, Basnyat B, Neupane M, Witonsky DB, Pritchard JK, et al. Admixture facilitates genetic adaptations to high altitude in Tibet. Nat Commun. 2014;5:3281.
Pierron D, Heiske M, Razafindrazaka H, Pereda-Loth V, Sanchez J, Alva O, et al. Strong selection during the last millennium for African ancestry in the admixed population of Madagascar. Nat Commun. 2018;9:1–9.
Yelmen B, Mondal M, Marnetto D, Pathak AK, Montinaro F, Gallego Romero I, et al. Ancestry-specific analyses reveal differential demographic histories and opposite selective pressures in modern south Asian populations. Mol Biol Evol. 2019;36:1628–42.
Bhatia G, Tandon A, Patterson N, Aldrich MC, Ambrosone CB, Amos C, et al. Genome-wide scan of 29,141 African Americans finds no evidence of directional selection since admixture. Am J Hum Genet. 2014;95:437–44.
Vitti JJ, Grossman SR, Sabeti PC. Detecting natural selection in genomic data. Annu Rev Genet. 2013;47:97–120.
Hermisson J, Pennings PS. Soft sweeps: molecular population genetics of adaptation from standing genetic variation. Genetics. 2005;169:2335–52.
Hermisson J, Pennings PS. Soft sweeps and beyond: understanding the patterns and probabilities of selection footprints under rapid adaptation. Methods Ecol Evol. 2017;8:700–16.
Hamblin MT, Di Rienzo A. Detection of the signature of natural selection in humans: evidence from the Duffy blood group locus. Am J Hum Genet. 2000;66:1669–79.
Pritchard JK, Pickrell JK, Coop G. The genetics of human adaptation: hard sweeps, soft sweeps, and polygenic adaptation. Curr Biol. 2010;20:R208–15.
Ronen R, Udpa N, Halperin E, Bafna V. Learning natural selection from the site frequency spectrum. Genetics. 2013;195:181–93.
Jouganous J, Long W, Ragsdale AP, Gravel S. Inferring the joint demographic history of multiple populations: beyond the diffusion approximation. Genetics. 2017;206:1549–67.
Prevention of neural tube defects: Results of the Medical Research Council Vitamin Study. Lancet. 1991;338:131–7.
Metz J. Folic Acid Metabolism and Malaria. Food Nutr Bull. 2007;28 4_suppl4:S540–9.
De Filippo C, Key FM, Ghirotto S, Benazzo A, Meneu JR, Weihmann A, et al. Recent selection changes in human genes under Long-term balancing selection. Mol Biol Evol. 2016;33:1435–47.
Gazi U, Martinez-Pomares L. Influence of the mannose receptor in host immune responses. Immunobiology. 2009;214:554–61.
East L, Isacke CM. The mannose receptor family. Biochim Biophys Acta Gen Subj. 2002;1572:364–86.
Martinez-Pomares L. The mannose receptor. J Leukoc Biol. 2012;92:1177–86.
Burgdorf S, Kautz A, Bohnert V, Knolle PA, Kurts C. Distinct Pathways of Antigen Uptake and Intracellular Routing in CD4 and CD8 T Cell Activation. Science (80- ). 2007;316:612–6.
Schuette V, Embgenbroich M, Ulas T, Welz M, Schulte-Schrepping J, Draffehn AM, et al. Mannose receptor induces T-cell tolerance via inhibition of CD45 and up-regulation of CTLA-4. Proc Natl Acad Sci. 2016;113:10649–54.
Hattori T, Konno S, Hizawa N, Isada A, Takahashi A, Shimizu K, et al. Genetic variants in the mannose receptor gene (MRC1) are associated with asthma in two independent populations. Immunogenetics. 2009;61:731–8.
Hattori T, Konno S, Takahashi A, Isada A, Shimizu K, Shimizu K, et al. Genetic variants in mannose receptor gene (MRC1) confer susceptibility to increased risk of sarcoidosis. BMC Med Genet. 2010;11:151.
Alter A, De Léséleuc L, Van Thuc N, Thai VH, Huong NT, Ba NN, et al. Genetic and functional analysis of common MRC1 exon 7 polymorphisms in leprosy susceptibility. Hum Genet. 2010;127:337–48.
Zhang X, Li X, Zhang W, Wei L, Jiang T, Chen Z, et al. The novel human MRC1 gene polymorphisms are associated with susceptibility to pulmonary tuberculosis in Chinese Uygur and Kazak populations. Mol Biol Rep. 2013;40:5073–83.
Dubois A, François C, Descamps V, Fournier C, Wychowski C, Dubuisson J, et al. Enhanced anti-HCV activity of interferon alpha 17 subtype. Virol J. 2009;6:70.
Elaldi N, Yilmaz M, Bagci B, Yelkovan I, Bagci G, Gozel MG, et al. Relationship between IFNA1 , IFNA5 , IFNA10 , and IFNA17 gene polymorphisms and Crimean-Congo hemorrhagic fever prognosis in a Turkish population range. J Med Virol 2016;88:1159–1167.
Voight BF, Kudaravalli S, Wen X, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006;4:0446–58.
Zwaenepoel I, Mustapha M, Leibovici M, Verpy E, Goodyear R, Liu XZ, et al. Otoancorin, an inner ear protein restricted to the interface between the apical surface of sensory epithelia and their overlying acellular gels, is defective in autosomal recessive deafness DFNB22. Proc Natl Acad Sci. 2002;99:6240–5.
Lee K, Chiu I, Santos-Cortez R, Basit S, Khan S, Azeem Z, et al. Novel OTOA mutations cause autosomal recessive non-syndromic hearing impairment in Pakistani families. Clin Genet. 2013;84:294–6.
Clark AG, Glanowski S, Nielsen R, Thomas PD, Kejariwal A, Todd MA, et al. Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios. Science. 2003;302:1960–3.
Yang Y, Quach C, Liang C. Autophagy modulator plays a part in UV protection. Autophagy. 2016;12:1677–8.
He S, Zhao Z, Yang Y, O’Connell D, Zhang X, Oh S, et al. Truncating mutation in the autophagy gene UVRAG confers oncogenic properties and chemosensitivity in colorectal cancers. Nat Commun. 2015;6:7839.
Visser M, Palstra R-J, Kayser M. Human skin color is influenced by an intergenic DNA polymorphism regulating transcription of the nearby BNC2 pigmentation gene. Hum Mol Genet. 2014;23:5750–62.
Chahal HS, Lin Y, Ransohoff KJ, Hinds DA, Wu W, Dai H-J, et al. Genome-wide association study identifies novel susceptibility loci for cutaneous squamous cell carcinoma. Nat Commun. 2016;7:12048.
Vernot B, Akey JM. Resurrecting Surviving Neandeltal Linages from Modern Human Genomes. Science (80- ). 2014;343 February:1017–21.
Poole LA, Cortez D. Functions of SMARCAL1, ZRANB3, and HLTF in maintaining genome stability. Crit Rev Biochem Mol Biol. 2017;52:696–714.
Ferrer-Admetlla A, Liang M, Korneliussen T, Nielsen R. On detecting incomplete soft or hard selective sweeps using haplotype structure. Mol Biol Evol. 2014;31:1275–91.
Enattah NS, Sahi T, Savilahti E, Terwilliger JD, Peltonen L, Järvelä I. Identification of a variant associated with adult-type hypolactasia. Nat Genet. 2002;30:233–7.
Li M, Tian S, Jin L, Zhou G, Li Y, Zhang Y, et al. Genomic analyses identify distinct patterns of selection in domesticated pigs and Tibetan wild boars. Nat Genet. 2013;45:1431–8.
Pestka S, Krause CD, Sarkar D, Walter MR, Shi Y, Fisher PB. Interleukin-10 and related cytokines and receptors. Annu Rev Immunol. 2004;22:929–79.
Li M, Liu X, Zhou Y, Su SB. Interferon- s: the modulators of antivirus, antitumor, and immune responses. J Leukoc Biol. 2009;86:23–32.
Muir AJ, Arora S, Everson G, Flisiak R, George J, Ghalib R, et al. A randomized phase 2b study of peginterferon lambda-1a for the treatment of chronic HCV infection. J Hepatol. 2014;61:1238–46.
Wack A, Terczyńska-Dyla E, Hartmann R. Guarding the frontiers: the biology of type III interferons. Nat Immunol. 2015;16:802–9.
Su Q, Hu F, Liu Y, Ge X, Mei C, Yu S, et al. Cryo-EM structure of the polycystic kidney disease-like channel PKD2L1. Nat Commun. 2018;9:1–12.
Ishimaru Y, Inada H, Kubota M, Zhuang H, Tominaga M, Matsunami H. Transient receptor potential family members PKD1L3 and PKD2L1 form a candidate sour taste receptor. Proc Natl Acad Sci. 2006;103:12569–74.
Huang AL, Chen X, Hoon MA, Chandrashekar J, Guo W, Tränkner D, et al. The cells and logic for mammalian sour taste detection. Nature. 2006;442:934–8.
Huque T, Cowart BJ, Dankulich-Nagrudny L, Pribitkin EA, Bayley DL, Spielman AI, et al. Sour ageusia in two individuals implicates ion channels of the ASIC and PKD families in human sour taste perception at the anterior tongue. PLoS One. 2009;4.
DeSimone JA, Heck GL, DeSimone SK. Active ion transport in dog tongue: a possible role in taste. Science. 1981;214:1039–41.
Fumagalli M, Moltke I, Grarup N, Racimo F, Bjerregaard P, Jørgensen ME, et al. Greenlandic Inuit show genetic signatures of diet and climate adaptation. Science. 2015;349:1343–7.
Ameur A, Enroth S, Johansson Å, Zaboli G, Igl W, Johansson ACV, et al. Genetic adaptation of fatty-acid metabolism: a human-specific haplotype increasing the biosynthesis of Long-chain Omega-3 and Omega-6 fatty acids. Am J Hum Genet. 2012;90:809–20.
Arciero E, Biagini SA, Chen Y, Xue Y, Luiselli D, Tyler-Smith C, et al. Genes regulated by vitamin D in bone cells are positively selected in East Asians. PLoS One. 2015;10:e0146072.
Han J, Colditz GA, Hunter DJ. Polymorphisms in the MTHFR and VDR genes and skin cancer risk. Carcinogenesis. 2007;28:390–7.
Jablonski NG, Chaplin G. Human skin pigmentation as an adaptation to UV radiation. Proc Natl Acad Sci. 2010;107 Supplement_2:8962–8.
Jones P, Lucock M, Veysey M, Beckett E. The vitamin D−Folate hypothesis as an evolutionary model for skin pigmentation: an update and integration of current ideas. Nutrients. 2018;10.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904–9.
Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19:1655–64.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Ronen R. Learning natural selection from the site frequency Spectrum. 2015. https://github.com/rronen/SFselect.
Robinson JD, Coffman AJ, Hickerson MJ, Gutenkunst RN. Sampling strategies for frequency spectrum-based population genomic inference. BMC Evol Biol. 2014;14.
Sabeti PC, Reich DE, Higgins JM, Levine HZP, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419:832–7.
Pickrell JK, Coop G, Novembre J, Kudaravalli S, Li JZ, Absher D, et al. Signals of recent positive selection in a worldwide sample of human populations. Genome Res. 2009;19:826–37.
Gautier M, Klassmann A, Vitalis R. Rehh 2.0: a reimplementation of the R package rehh to detect positive selection from haplotype structure. Mol Ecol Resour. 2017;17:78–90.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38:e164.
Haller BC, Messer PW. SLiM 2: flexible, interactive forward genetic simulations. Mol Biol Evol. 2017;34:230–40.
Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. 2009;25:1091–3.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
The authors thank the comments of the anonymous referees that have made possible the improvement of the manuscript.
About this supplement
This article has been published as part of BMC Genetics Volume 21 Supplement 1, 2020: Selected Topics in “Systems Biology and Bioinformatics” - 2019: genetics. The full contents of the supplement are available online at https://bmcgenet.biomedcentral.com/articles/supplements/volume-21-supplement-1.
This study has been possible thanks to the F.P.I. grant BES-2014-068994 to SW, and grant BFU2016–77961-P (AEI/FEDER, UE) both awarded by the Agencia Estatal de Investigación (MINECO, Spain) and with the support of Secretaria d’Universitats i Recerca del Departament d’Economia i Coneixement de la Generalitat de Catalunya (GRC 2017 SGR 702). Part of the “Unidad de Excelencia María de Maeztu”, funded by the AEI (CEX2018–000792-M). YX and CTS are supported by Wellcome Trust (098051), LP is supported by the European Union through the European Regional Development Fund Project No. 2014–2020.4.01.16–0024, MOBTT53. Publication costs were funded by Wellcome (grant 098051). The funding body has not played any role in the design of the study and collection, analysis, and interpretation of data and in writing of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Location of the five sampled populations. Supplementary Fig. 2. Genome-wide Manhattan plots of SFselect scores of the five populations of study. Supplementary Fig. 3. Density plots of SFselect and iHS scores of neutral and real data. Supplementary Fig. 4. Principal component analysis of the five East African populations. Supplementary Fig. 5. ADMIXTURE analysis of the Ethiopian samples and a set of worldwide populations. Supplementary Fig. 6. Genome-wide Manhattan plots of the –log10(p-value) of iHS of the five populations of study. Supplementary Fig. 7. PCA of the masked East African samples with a set of Europeans and Africans. Supplementary Fig. 8. Schematic representation of the demographic model used to simulate neutral sequences. Supplementary Fig. 9. Relative site frequency spectrum of Afroasiatic and Gumuz populations.
The 99.99th and 99.90th percentile thresholds of SFselect and iHS calculated after the neutral simulations. Supplementary Table 2. SFselect positive selection signals found in the five populations of study. Supplementary Table 3. SFselect positive selection signals found in the five populations of study. Supplementary Table 4. iHS positive selection signals found in the five populations of study. Supplementary Table 5. iHS positive selection signals shared among the five populations of study.
. Information about additional examples of shared and population-specific signals of positive selection.
About this article
Cite this article
Walsh, S., Pagani, L., Xue, Y. et al. Positive selection in admixed populations from Ethiopia. BMC Genet 21 (Suppl 1), 108 (2020). https://doi.org/10.1186/s12863-020-00908-5