Skip to main content
Download PDF

Changes in selective pressures associated with human population expansion may explain metabolic and immune related pathways enriched for signatures of positive selection

BMC Genomics volume 17, Article number: 504 (2016) Cite this article

Abstract

Background

The study of local adaptation processes is a very important research topic in the field of population genomics. There is a particular interest in the study of human populations because they underwent a process of rapid spatial expansion and faced important environmental changes that translated into changes in selective pressures. New mutations may have been selected for in the new environment and previously existing genetic variants may have become detrimental. Immune related genes may have been released from the selective pressure exerted by pathogens in the ancestral environment and new variants may have been positively selected due to pathogens present in the newly colonized habitat. Also, variants that had a selective advantage in past environments may have become deleterious in the modern world due to external stimuli including climatic, dietary and behavioral changes, which could explain the high prevalence of some polygenic diseases such as diabetes and obesity.

Results

We performed an enrichment analysis to identify gene sets enriched for signals of positive selection in humans. We used two genome scan methods, XPCLR and iHS to detect selection using a dense coverage of SNP markers combined with two gene set enrichment approaches. We identified immune related gene sets that could be involved in the protection against pathogens especially in the African population. We also identified the glycolysis & gluconeogenesis gene set, related to metabolism, which supports the thrifty genotype hypothesis invoked to explain the current high prevalence of diseases such as diabetes and obesity. Extending our analysis to the gene level, we found signals for 23 candidate genes linked to metabolic syndrome, 13 of which are new candidates for positive selection.

Conclusions

Our study provides a list of genes and gene sets associated with immunity and metabolic syndrome that are enriched for signals of positive selection in three human populations (Europeans, Africans and Asians). Our results highlight differences in the relative importance of pathogens as drivers of local adaptation in different continents and provide new insights into the evolution and high incidence of metabolic syndrome in modern human populations.

Background

Migration and colonization of new habitats are a major component of human demographic history. Such processes lead to sharp changes in environmental conditions experienced by individuals (soil, water, climate, and pathogens) with concomitant changes in lifestyle and diet, which in turn alter the phenotypic composition of human populations. Changes in environmental conditions can lead to shifts in selective pressures. Thus, new mutations may be advantageous or previous favorable ones may become deleterious [1, 2].

Diseases caused by deleterious mutations are not expected to persist for many generations or to segregate at high frequency. Yet, the prevalence of metabolic syndrome such as diabetes type 2 and obesity has increased dramatically in recent years, especially in industrialized countries [3, 4]. In principle, natural selection should act to eliminate harmful genetic variants but there are many reasons why this is not always the case. One of the most well-known explanations for the high incidence of metabolic syndrome is the “Thrifty Genotype Hypothesis” [5], which posits that alleles that had a selective advantage in past environments became deleterious in the modern world due to external stimuli including climatic, dietary and behavioral changes. A non-adaptive alternative to this hypothesis is the “Drifty Gene Hypothesis”, which contends that early hominins were subject to predation risks that selected against mutations associated with increased weigh. However, the evolution of social behavior, discovery of fire and use of weapons released more recent Homo species from such selective pressure allowing metabolic disease genes to evolve under neutrality [6]. Also, there are many other possible explanations. For example, a deleterious mutation may be driven to high frequency due to linkage disequilibrium with a neighboring locus/region that is under positive selection [7, 8]. Also, an allele might have pleiotropic effects, for example increased disease risk and protection against a bacterial infection [9].

Another good example of changes in selective pressures exerted on humans during their spatial expansion relates to the pathogens they found in the newly colonized habitats. Infectious diseases caused by pathogens are one of the main threats to human populations and there is an increasing interest in the identification of genomic regions associated with parasite resistance or susceptibility [10, 11]. It is well known that the evolutionary arms race between hosts and pathogens can lead to selective pressures on immune related gene sets [11–16]. Recent studies provide several examples of pathogen-driven selection signatures with important implications in terms of human immunology and infectious disease epidemiology [10]. Thus, it is important to investigate if the application of new statistical genetic approaches can uncover additional examples.

The purpose of the present study is to investigate the extent to which colonization of new habitats by humans represent important changes in selective pressures that may explain human phenotypic variation related to common diseases and immunity. In particular, changes in diet are likely to represent important new selective pressures that may explain the high prevalence of diabetes type 2. Additionally, differences in parasite species prevalent in different habitats may also lead to changes in selective pressures on immunity related genes. Thus, we focus on the impact that environmental changes may have on metabolic syndrome and infectious diseases. Until recently, a substantial understanding of the evolution of such complex diseases has been derived from candidate variants that were revealed either by experiments or gene association studies [17, 18]. There is general agreement that complex diseases are inherently polygenic and, therefore, genome scans based solely on SNP data are not well suited for making inferences about their genetic architecture. Individual variants associated with complex traits tend to be weakly selected and, therefore, only the combined effect on gene sets can be strong enough to reach strict cutoff levels of significance. This limitation can be overcome by combining genome scans of SNP data with gene-set enrichment analyses [19]. This approach detects selection signatures at the gene set level and, therefore, is better adapted to detect polygenic selection. For this reason, we first carry out genome scans of selection (XPCLR [20] and iHS [21]) using HapMap SNP data [22] and then we use their output as input of Gene Set Enrichment Analysis (GSEA-Daub et al. [19] and Gowinda [23]). We mainly focus on the gene set level and we refine it to the gene level whenever needed.

We detected gene sets involved in immunity and host-defense interactions in African. We further examined how natural selection has influenced the prevalence of diseases such as diabetes type 2 and obesity, investigating at both the gene and gene set level. Our results indicate that changes in environmental conditions have lead to shifts in selective pressures on genomic regions associated with adaptation to pathogens and with phenotypes associated to metabolic syndrome.

Methods

SNP data

We use SNP data from Hapmap database phase II [22] to test for signatures of positive selection using two genome-scan methods, XPCLR [20] and iHS [21] (the data from HapMap Project are publicly available for unrestricted research use). These methods provide selection scores for each of the SNPs in the dataset, which we later use to perform the enrichment analyses.

The data consists of three populations: Europeans (CEU), Yoruba from Ibadan, Nigeria (YRI), and Han Chinese from Beijing, China and Japanese from Tokyo (CHB + JTP). Recombination rates and genetic coordinates are from assembly NCBI36 [24].

Genome scan methods

Several methods have been developed to detect signals of positive selection. Fst-based [25] and haplotype homozygosity methods [21, 26] are among the most widely used. Although these methods detect genomic regions under positive selection, they do not necessarily detect the same signals. A sensitivity analysis [27] that compared the performance of many recent methods under complex evolutionary scenarios revealed that the Cross Population Composite Likelihood Ratio test (XPCLR) [20] performs the best under both hard and soft sweep scenarios. However, it lacks power to detect very recent (incomplete) selective sweeps, in which case the integrated Haplotype Score (iHS) [21] performs better. More specifically, iHS has the ability to detect selective sweeps at an early stage but lacks power to detect those at intermediate or late stages while XPCLR exhibits the exact opposite behaviour. Another substantial difference between the two genome-scan methods is the underlying theory that they use to identify candidate regions under positive selection. XPCLR is a population differentiation method, which uses the multilocus allele frequency differentiation between two populations; the objective population (under positive selection) and the reference population (under neutrality). On the other hand, iHS is a haplotype-based method that considers the extent of linkage disequilibrium surrounding the core SNP in an isolated population. Therefore, we chose to combine XPCLR and iHS to maximize the power to detect selection signatures at the SNP level.

XPCLR analyses

We carried out XPCLR analyses of all possible population pairs from the Hapmap Phase II dataset. In all of the comparisons, the first population was considered as the objective population (population under positive selection) and the second one was the reference population (non-selected population). Thus, the population pairs considered were: CEU - YRI, YRI - CEU, CEU - CHB + JTP, CHB + JTP - CEU, CHB + JTP - YRI and YRI - CHB + JTP. A previous study that evaluated the performance of several genome scan methods [27] shows that the FDR of XPCLR increases substantially once the SNP reaches fixation. Therefore, we only used SNPs that were polymorphic in both populations in order to minimize false positives. Thus, there were 2,179,305 SNPs for population pairs including CEU and CHB + JTP, 2,203,610 SNPs for those that considered CEU and YRI and 2,118,211 SNPs for pairs including YRI and CHB + JTP. In principle, since humans originate in Africa, the YRI population should be considered as the reference population when looking for signatures of positive selection due to changed environmental conditions. However, we also carried out analyses with YRI as the objective population, as variants that were positively selected in Africa may be released from selective pressure in newly colonized habitats (Europe and Asia). We evaluated grid points every 200 bp along the genome and evaluated a window size of 50 SNPs around each point. We normalized the SNP scores across the whole genome such that they have zero mean and unit variance. This was performed separately for each of the pairwise comparisons conducted in the dataset.

iHS analyses

iHS scores for the three populations were derived for 1,131,343 SNPs for CEU, 1,264,811 for YRI and 976,006 for CHB + JTP. iHS is based on the decay in Extended Haplotype Homozygosity (EHH) to the right and left of a focal SNP, which is measured by the integrated Haplotype Homozygosity (iHH). More precisely, the EHH is calculated for all SNPs to the right and left of the focal SNP until EHH is equal to or larger than 0.05. However, if the focal SNP is close to the end of the chromosome or the start of a gap >200 kb, there may be no SNP for which EHH reaches or is below that value, in which case the iHS statistic cannot be calculated.

iHS performs better when the frequency of the selected allele is low (~0.1 and ~0.3) [27]. As in the case of the XPCLR analyses, iHS scores were also normalized. More details about the genome scan methods are given in Additional file 1: Text 1.

The results of the two genome-scan methods (XPCLR and iHS) for all the chromosomes in the HapMap Phase II database are available under request.

Gene set enrichment analysis (GSEA)

After we calculated the selection scores (from XPCLR and iHS) for each of the SNPs in the dataset for all the different populations, we performed the GSEA. We used Daub et al. [19] approach and Gowinda to detect significant gene sets (more details are presented in Additional file 1: Text 2). This procedure was performed twice, using a) the XPCLR scores and b) the iHS scores.

Gene data

We downloaded the human protein coding gene data from the Entrez NCBI database from assembly GRCh38 on May 5/2014 and we extracted the start and the end position of 27,121 genes. For all genes that possess more than one reference sequence, we took the utmost start and end position of each of them (longest isoform).

Assignment of SNPs to genes

To match the SNP positions from assembly NCBI36 to gene positions (start/end) from assembly GRCh38, we converted the SNP positions from assembly NCBI36 to assembly GRCh38 using two tools: a) Remap tool from NCBI [28] and b) the liftOver tool from UCSC [29]. Conversion with the Remap tool was made on 5/3/2014 and conversion with liftOver on 10/3/2014. The SNPs that were mapped to the same location with both tools were kept for the final enrichment analysis.

Gene sets

We downloaded all the human gene sets from the Biosystems database [30] on 6/5/2014. The NCBI Biosystem database includes gene sets from KEGG, REACTOME, WikiPathways, Pathway Interaction Database (PID) and BIOCYC. We discarded gene sets that contain less than 10 genes [31]. To avoid redundancy due to duplicated gene sets, we merged gene sets with a similarity larger than 95 %. After these corrections, we were left with 928 gene sets out of the 2682 initial ones to continue the analysis.

Gene set enrichment analyses (GSEA) methods

Both Daub et al. [19] approach and Gowinda calculate the significance of gene sets based on permutation tests. They avoid biases due to differences in the length of genes and the overlap of the genes between gene sets. Nevertheless, the two methodologies are expected to give different results as they differ in several aspects, the most important being in the calculation of the gene set selection scores. Daub et al. [19] uses the selection scores from all SNPs in the analysis, estimating the sum of selection scores of all the genes in the gene set. On the other hand, Gowinda focuses only on the selection scores of the candidate SNPs, which are those that are significant for the chosen cut-off significance threshold (usually the 1- or 5-percentile of the total number of SNPs used in the genome scan). Other important differences include the permutation procedure, the correction of overlapping genes, and the correction for multiple testing. More details about the two GSEA approaches are given in Additional file 1: Text 2. To determine the significance in the real dataset, we used a q-value smaller or equal to 0.09 for both GSEA approaches.

Gene enrichment analysis for metabolic syndrome

To evaluate more in detail if the “Thrifty genotype hypothesis” can help explain the high incidence of diabetes 2, obesity, and metabolic syndrome in modern human populations we identified all genes that have been associated with these diseases using the Bio4j tool and the STRING database (see below). We then determined whether or not these genes were enriched for signatures of positive selection using the results of the XPCLR and iHS genome scans. A threshold of 1 % was used to determine significance of the gene scores.

Bio4j

Using the Bio4j tool [32], a graph-based platform that integrates semantically biological data from 6 different databases (Uniprot KB –SwissProt + Trembl, Gene Ontology, UniRef (50,90,100), NCBI Taxonomy, and Expasy Enzyme DB), we identified a total of 683 genes associated with obesity, diabetes 2 and metabolic syndrome. More details for Bio4j analysis are available in the Additional file 1: Text 3.1.

STRING database

We also used the Search Tool for the Retrieval of Interacting Genes (STRING) database [33] to find genes that interact with the candidate genes identified using Bio4j. The Protein-Protein Interaction (PPI) networks can reveal additional genes associated with metabolic syndrome, which could then be evaluated for signatures of positive selection. More details for STRING analysis are given in the Additional file 1: Text 3.2.

Results and Discussion

In this study, we use the Cross Population Composite Likelihood Ratio (XPCLR) [20] and the integrated Haplotype Score (iHS) [21] to uncover the broadest range of positive selection signatures in the human genome. Numerous gene sets were identified to be involved in local adaptation but we focus on those related to immune responses and metabolic syndrome.

We found a total 161 gene sets enriched for signals of positive selection; 15 identified by Daub et al. [19] approach and 152 by Gowinda, with six that were identified by both. Figure 1 summarizes the results of the four GSEAs and the common gene sets that were found among them. An interesting observation is the large number of significant gene sets based on the iHS scores and Gowinda’s enrichment analysis (140). Using XPCLR, we performed six pairwise analyses of the three populations (Europeans, Chinese/Japanese, Yoruba) to investigate all the possible selection signals.

Fig. 1
figure 1

Illustration of the 167 gene sets that were detected with the four combinations of the genome scan and GSEA methods. Each circle shows the methods that were used with the number of candidate gene sets. The overlap among them is also shown

Full size image

When iHS scores were used as input of the enrichment analysis, we found no significant gene sets for the CHB + JTP population with Gowinda and only one with Daub et al. [19] approach. In contrast, we identified a total of 140 different significant gene sets for the YRI population. Finally, the number of enriched gene sets obtained for the CEU population was four only (Table 1 and Additional file 2: Tables S1 and S2).

Table 1 Summary of the results of the enrichment analysis using iHS scores as the baseline
Full size table

When XPCLR scores were used as input for the enrichment analysis, we identified a total of 17 significant gene sets across all population pairs and enrichment analyses approaches (Table 2; Additional file 2: Tables S3 and S4). Using Daub et al. [19] approach, we found no significant gene sets for the population pairs CEU-CHB + JTP and YRI-CEU, while with Gowinda, we found no significant gene sets for the CEU-YRI population pair. Interestingly, when YRI was the objective population and CHB + JTP was the reference one, we detected seven gene sets, more than in the comparisons where YRI was the reference. We obtained similar results when YRI was the objective population and CEU was the reference, but only when using Gowinda.

Table 2 Summary of the results of the enrichment analysis using XPCLR scores as the baseline
Full size table

Overall, these results suggest that selective pressures are much more important in Africa than in European and Asian populations. We posit that many of the enriched gene sets were under selection in Africa but were released from selection in the newly colonized habitats (Europe and Asia).

In what follows, we first present the results of the GSEAs based on XPCLR and then the results of analyses based on iHS mainly for the Yoruba population. We also explore candidates of positive selection that could be related to metabolism and therefore have been involved in the evolution of diseases such as diabetes type 2, obesity and metabolic syndrome.

Infectious diseases

As individuals migrated to new environments (e.g. climate, diet), they were forced to adapt, and natural selection acted to favor variants that were advantageous both in pathogens and hosts. Both environmental factors and pathogen load must have had a great impact on the human genetic variation and infectious disease susceptibility. Based on all the analyses carried out, we identify a total 35 gene sets related to immune protection and host defense interactions; seven of them were detected with XPCLR and 28 with iHS (details about the populations that the gene sets were detected are presented in Additional file 2: Table S5).

XPCLR-based GSEA

When Africa was compared with the Asian population (YRI: objective population), we observed increased selection pressures on host-defense mechanisms that we describe below in more detail (Additional file 2: Table S5).

Of primary importance is the Cell surface interactions at the vascular wall gene set, which was detected with both GSEA approaches. Platelets and leukocytes interact with the endothelium as a protective response of the host to infection by bacteria, pathogens or injury. Platelets can also enhance their activation through the generation of signals by the alphaIIbbeta3 integrin gene set [34], also significant in our analysis. Key regulator of host defense is the T-Cell-Receptor (TCR) gene set, whose main function is the induction of responses to the cell nucleus to isolate and destroy the malignant cells [35, 36]. Also significant was the chemokine gene set, which is directly involved in both innate and adaptive immune responses, as secreted chemokines and their receptors mediate leukocyte recruitment in the site of infection [37]. Xu et al. (2014) suggest that CXCRs, main receptors of the chemokine subfamily CXC, are also enriched for positive selection [38].

Although four immune mediated gene sets were detected in the pairwise comparison between YRI-CHB + JTP, in the opposite comparison (CHB + JTP-YRI) we detected only one, the apoptotic gene set, involved in programmed cell death. The dysfunction of genes involved in apoptosis can lead to inflammatory diseases and cancer [39–42]. da Fonseca et al. (2010) shows that positive selection signals on the apoptotic genes in mammals are mainly the result of evolutionary forces between host and pathogens [43].

These observations overall suggest that gene sets related to host defense interactions in the immune system were under selective pressures in Africa. However, they were released from such pressures once humans migrated to Asia, presumably because many of the pathogens (bacteria, viruses and other infectious pathogens) were absent in the newly colonized continent.

iHS-based GSEA

Interestingly, numerous gene sets (147) were detected as enriched for signatures of positive selection in the YRI population using iHS. Here we mainly focus on the immune- and host defense- related ones.

Ten gene sets associated with different phases of Human Immunodeficiency Virus (HIV) and HIV-1 infection were detected in the African population (Additional file 2: Table S5). HIV is among the most threatening infectious diseases in Africa and its prevalence accounts for 71 % of HIV worldwide [44]. The genetic diversity in HIV patients is believed to be the result of the virus evolution in relation to the human immune system [45, 46]. Positive selection was also detected in previous studies on the human genome of HIV patients [47] as well as on HIV genome sequences [45, 48, 49], indicating selective pressures on both the virus and the host.

Three other significant gene sets, the Salmonella, Shigellosis, and Helicobacter pylori (H. pylori) infection gene sets are clearly involved in the defense against pathogens. These infections are usually caused by environmental contamination (e.g. food or water) [50–52]. Bacteria continuously evolve to find new strategies to deceive the host barriers [53]. Positive selection on the specific bacteria strains [54–56] is consistent with an evolutionary arm race between bacteria and the human host immune system. Salmonella and Shigella use the same mechanism, known as the trigger model, to invade organisms [57], whereas H. pylori adheres to and reproduces in epithelial cells [58]. Interestingly, we also detected selection pressures on the invasion of bacteria through epithelial cells gene set. Epithelial cells are involved in antimicrobial host immune defense and induce inflammatory responses to prevent the pathogens from invading [59].

The Toll-Like-Receptors (10 different cascades were significant) and RIG-1 receptors (2 different significant pathways) play a critical role to introduce innate immune responses. They are associated with the recognition and defense of pathogenic microorganisms, viruses and bacteria [60, 61]. Two Interleukin gene sets, a group of cytokines that secret signals to regulate the immune response, were also identified as candidates of positive selection.

Overall, our results indicate that selective pressures on host-defense mechanisms exerted by pathogens in Africa became much less important once humans colonized other continents because of dramatic changes in the environment experienced by colonizers. In agreement with previous observations, selection is associated with protective immunity (e.g. TLRs- RIG-1 receptors, interleukin) against infectious diseases (e.g. bacterial infections, HIV) [19, 61–66].

Obesity – metabolic syndrome – diabetes

Obesity, diabetes and metabolic syndrome are chronic diseases that have been thoroughly studied. A large effort has been made to find variants, genes and gene sets that contribute to their pathogenesis. Although many studies have identified multiple genetic factors, few have examined whether or not they are enriched for signatures of positive selection. Past selective pressures could help explain the high prevalence of metabolic syndrome in modern human populations [67–70]. More precisely, changes in environmental conditions, nutrition, or lifestyle may induce radical changes in selective patterns; genes initially subject to positive selection in ancient humans may become deleterious under the environmental conditions experienced by modern human societies. Potential examples are the insulin receptor substrate-1 (IRS1) and TCF7L2 genes that have already been found to be under positive selection [71, 72]. Nevertheless, given the likely polygenic architecture of metabolic syndrome, it is important to carry out gene set enrichment analyses in order to identify new candidate gene sets. In our study, we found two gene sets that could be associated with metabolic syndrome; the Glycolysis and gluconeogenesis and the Signal Attenuation gene sets.

The Glycolysis and gluconeogenesis gene set was detected with XPCLR and Gowinda when the European was the objective population and the Asian was the reference population. In principle, glycolysis and gluconeogenesis are two metabolic processes that have opposite effects. Glycolysis degrades glucose in periods of fasting, diet or intense exercise and gluconeogenesis produces glucose in periods of feeding. The role of glycolysis and gluconeogenesis gene set is well defined in the pathogenesis of diabetes. Glycolysis controls the hepatic glucose production [73], which causes hyperglycemia, a symptom of diabetes. It is also involved in cases of insulin resistance and deficiency [74], which is usually observed in type 2 diabetes [75]. Felig et al. (1974) also indicates that glucose precursors are increased in obese population [76], although the clear role of this gene set in obesity is still controversial [77]. A few studies have already investigated the possibility of glycolysis and gluconeogenesis gene set as targets for diabetes therapies [78–80].

Our result are consistent with the hypothesis that in the past, the glycolysis and gluconeogenesis gene set was subject to positive selective pressures associated with adaptation to potentially long fasting periods. But it then may have become maladaptive due to extreme changes in the environment experienced by humans in modern societies. We suggest that these results could be explained by the fact that recent negative selection should be pronounced in Europeans but not in Asians presumably because of differences in diet and lifestyle (e.g. food, increased sedentarily and limited exercise). For example, Europeans tend to consume more meat than Asians [81], and it has been shown that increased consumption of meat increases the incidence of diabetes 2 [82].

We also detected the signal attenuation gene set in the European population, which includes 15 insulin related genes. This gene set was identified in the comparison with the Africans (CEU-YRI pair) using XPCLR and Daub et al. [19] approach. We note that the significance of this gene set may be due mostly to the top-scoring outlier gene, DOK1, which plays an important role in insulin receptor binding [83]. Nevertheless, insulin related genes closely interact with DOK1 [84], according to the STRING database (Additional file 1: Figure S2). Among them, the IRS1 and IRS2 genes play an important role in the pathogenesis of diabetes and obesity [85–88]. Our analysis did not detect the IRS- genes presumably because of the stringent cut- off level that we used, however, a previous study has shown that IRS1 is under positive selection [72].

Lastly, we identified the autoimmune thyroid disease in the pairwise comparison between Europeans (objective) and Yoruba, which is an inflammatory disease in the thyroid gland where blood cells and antibodies attack healthy tissue. Interestingly, it has also been associated with the pathogenesis of most types of diabetes [89] and its prevalence shows a ~4 % increase in diabetic populations compared to non-diabetic. Our findings agree with Jin et al. (2012) [90], that also detected the auto-immune thyroid disease and allograft rejection gene sets enriched for positive selection in populations from Sub-Saharan Africa and North America using an Fst-based genome scan method.

Having identified gene sets associated with metabolic syndrome and enriched for signatures of positive selection, we next investigated if individual genes potentially associated with this disease, based on current knowledge, were also enriched for signatures of positive selection. To this end we used Bio4j to extract the genes that are associated to metabolic syndrome and then used STRING to identify genes that interact closely with them.

Using this approach, we identified 23 genes associated with metabolic syndrome and enriched for signatures of positive selection; 17 of them were detected with the XPCLR-based analysis and six other were detected with iHS (Additional file 2: Tables S6 and S7 respectively). Ten of them (SLC19A2, SLC2A10, UCP2, SLC27A4, BLK, ESP15, EGFR, SOCS1, TSC, and PTH) have already been suggested to be under positive selection [15, 91–101]. That leaves us with 13 new genes associated directly (six) and indirectly (seven) with obesity, diabetes or metabolic syndrome and also subject to natural selection. Interestingly, the KDM3A (or Jhdm2a) gene also enriched for positive selection is a key regulator for obesity and metabolic syndrome in mice [102, 103], but it has not been studied thoroughly in human populations. Nine of the 23 genes play a potential risk role in the development of diabetes, obesity or metabolic syndrome. Five of them (GATA6, CP, LARS2, MRAP2 and SGIP1) belong to the 13 novel genes that were identified in the current analyses. Four in total have a potential protective effect, with the PIK3CB gene to be a novel candidate for positive selection, detected in the comparison of Asians (objective) with Europeans (reference). PIK3CB was previously shown to be associated with protection against insulin resistance [104]. The rest of the genes have an indirect association without clear evidence for their direct implication in metabolic syndrome (Additional file 2: Tables S6 and S7).

Comparing our results with two previous analyses on diabetes [68, 69], we found that TSPAN8 and NOTCH2 genes were also detected in our analyses, while some others (ADAMTS9, JAZF1, GRB14 and PRC1) genes were just below the threshold we used. Figure 2 presents the genes that show evidence of positive selection, and are risk or protective factors for metabolic syndrome. Surprisingly, we do not detect the TCF7L2 gene, unlike many studies [71, 105]. Nevertheless, the TBL1XR1 and CUL1 genes that closely interact with TCF7L2 gene (Additional file 1: Figure S3; derived from STRING database) are under positive selection in our analysis.

Fig. 2
figure 2

Strength of positive selection (XPCLR scores in the y-axis) on genes that are associated with metabolic syndrome. The blue line illustrates the 1 % cut off for each pair. The population pairs where these genes were detected are shown on the top. The first population is always the objective and the second one is the reference. Genes with no star symbol were suggested to be susceptible to metabolic syndrome, with one star, they need further investigation for their association with metabolic syndrome and with two stars, they might offer protection against metabolic syndrome

Full size image

Conclusions

In this study, we were particularly interested in uncovering selection signatures due to changed environmental conditions found by expanding human populations as they left the African continent. We hypothesized that such changes can influence gene sets involved in immune response to pathogens or associated with diseases such as diabetes, obesity and metabolic syndrome, which could be connected with changes in life-style or diet.

It is widely accepted that after the first exodus from Africa, humans confronted changed environmental conditions and had to face new pathogens [106, 107]. Therefore, migrants were likely to encounter new selective pressures. The large number of positively selected gene sets that we found in the African population, mainly using the iHS genome scan method, suggests that the most common scenario involved a release from selective pressures present in the African continent but absent in Europe and Asia. Indeed, our results show that selective pressures on immune related- and host-defense gene sets are detected in Africa but not in Eurasia. This result is consistent with the idea that pathogens may be more abundant in Africa than in other continents [108].

It is important to note that the signature of selection identified by iHS (long stretches of homozygosity) is similar to that left by purifying selection acting on conserved regions [109]. In this regard, we observed that gene sets uncovered by both Daub et al. [19] approach and Gowinda when using the iHS selection scores (spliceosome, cell cycle, mitotic and DNA repair) are conserved across several species. The distribution of the iHS scores of some of the conserved pathways (Additional file 1: Figure S3), does not give a clear pattern whether selection in the pathways are mainly driven by selection on the ancestral or the derived allele. Thus, although iHS is widely used to detect positive selection, we posit that several of the gene sets identified by iHS might in fact be subject to purifying selection instead of local adaptation. However, according to previous studies [109, 110], this happens mainly in non-coding regions, unlike the gene sets detected here that are in coding regions.

We investigated thoroughly the gene sets and the genes that could be associated with metabolic syndrome, mainly diabetes and obesity, whose prevalence increased in recent years. We examine diabetes and obesity as a single entity because they are intimately associated. More specifically, obesity is considered one of the most important factors in the pathogenesis of diabetes [111]. Our results suggest that metabolic syndrome could in part be explained by the thrifty genotype hypothesis. Even though, there is a controversy on whether or not insulin resistance increased metabolic efficiency in the past [112], our result is in agreement with the reduced selection pressure on the Glycolysis and Gluconeogenesis gene set that was observed in the Asian population (see above). The process of gluconeogenesis is controlled by insulin hormone but in the case of insulin resistance, it fails to stop glucose production [74]. Therefore, there might be an association between the deficiency of the insulin receptor and increased levels of nutrients that promoted fat storage in early humans. Lastly, two potentially protective genes, EGFR and SOCS1, were identified with iHS, which suggests that these genes may be undergoing recent selective sweeps associated with the high prevalence of metabolic syndrome in modern humans.

Our analyses uncovered several new gene sets with immune- and glucose- related functions that may be subject to positive selection. Moreover, it presents a substantial number of novel genes associated with diabetes and obesity that are enriched for signatures of positive selection. Overall, our study brings new insight into the emergence and evolutionary history of infectious diseases and metabolic syndrome. In the era of unlimited public data resources and cutting-edge approaches to analyze them, understanding the way that polygenic selection has shaped the human genetic diversity could lead towards better prevention and treatment of these diseases.

Abbreviations

CEU, Europeans; CHB + JTP, Han Chinese from Beijing, China and Japanese from Tokyo; GSEA, gene set enrichment analysis; iHS, integrated haplotype score [20]; XPCLR, Cross Population Composite Likelihood Ratio (XPCLR) [23]; YRI, Yoruba from Ibadan, Nigeria.

References

  1. Innan H, Kim Y. Pattern of polymorphism after strong artificial selection in a domestication event. Proc Natl Acad Sci U S A. 2004;101:10667–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Przeworski M, Coop G, Wall JD. The signature of positive selection on standing genetic variation. Evolution. 2005;59:2312–23.

    Article  PubMed  Google Scholar 

  3. Armstrong GL, Conn LA, Pinner RW. Trends in infectious disease mortality in the United States during the 20th century. Jama. 1999;281:61–6.

    Article  CAS  PubMed  Google Scholar 

  4. Chan RS, Woo J. Prevention of overweight and obesity: how effective is the current public health approach. Int J Environ Res Public Health. 2010;7:765–83.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Neel JV. Diabetes Mellitus: A “Thrifty” Genotype Rendered Detrimental by “Progress”? Am J Hum Genet. 1962;14(4):353–62.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Speakman JR. Evolutionary perspectives on the obesity epidemic: adaptive, maladaptive, and neutral viewpoints. Annu Rev Nutr. 2013;33:289–317.

    Article  CAS  PubMed  Google Scholar 

  7. Barreiro LB, Quintana-Murci L. From evolutionary genetics to human immunology: how selection shapes host defence genes. Nat Rev Genet. 2010;11:17–30.

    Article  CAS  PubMed  Google Scholar 

  8. Holopainen PM, Partanen JA. Technical note: linkage disequilibrium and disease-associated CTLA4 gene polymorphisms. J Immunol. 2001;167:2457–8.

    Article  CAS  PubMed  Google Scholar 

  9. Zhernakova A, Elbers CC, Ferwerda B, Romanos J, Trynka G, Dubois PC, de Kovel CG, Franke L, Oosting M, Barisani D, et al. Evolutionary and functional analysis of celiac risk loci reveals SH2B3 as a protective factor against bacterial infection. Am J Hum Genet. 2010;86:970–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Fumagalli M, Sironi M. Human genome variability, natural selection and infectious diseases. Curr Opin Immunol. 2014;30:9–16.

    Article  CAS  PubMed  Google Scholar 

  11. Karlsson EK, Kwiatkowski DP, Sabeti PC. Natural selection and infectious disease in human populations. Nat Rev Genet. 2014;15:379–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Baker C, Antonovics J. Evolutionary determinants of genetic variation in susceptibility to infectious diseases in humans. PLoS One. 2012;7, e29089.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Corona E, Dudley JT, Butte AJ. Extreme evolutionary disparities seen in positive selection across seven complex diseases. PLoS One. 2010;5, e12236.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Karlsson EK, Harris JB, Tabrizi S, Rahman A, Shlyakhter I, Patterson N, O’Dushlaine C, Schaffner SF, Gupta S, Chowdhury F, et al. Natural selection in a bangladeshi population from the cholera-endemic ganges river delta. Sci Transl Med. 2013;5:192ra186.

    Article  Google Scholar 

  15. Ramos PS, Shaftman SR, Ward RC, Langefeld CD. Genes associated with SLE are targets of recent positive selection. Autoimmune Dis. 2014;2014:203435.

    PubMed  PubMed Central  Google Scholar 

  16. Sironi M, Cagliani R, Forni D, Clerici M. Evolutionary insights into host-pathogen interactions from mammalian sequence data. Nat Rev Genet. 2015;16:224–36.

    Article  CAS  PubMed  Google Scholar 

  17. Estrada K, Aukrust I, Bjorkhaug L, Burtt NP, Mercader JM, Garcia-Ortiz H, Huerta-Chagoya A, Moreno-Macias H, Walford G, et al. Association of a low-frequency variant in HNF1A with type 2 diabetes in a Latino population. Jama. 2014;311:2305–14.

    Article  PubMed  CAS  Google Scholar 

  18. Scuteri A, Sanna S, Chen WM, Uda M, Albai G, Strait J, Najjar S, Nagaraja R, Orru M, Usala G, et al. Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits. PLoS Genet. 2007;3, e115.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Daub JT, Hofer T, Cutivet E, Dupanloup I, Quintana-Murci L, Robinson-Rechavi M, Excoffier L. Evidence for polygenic adaptation to pathogens in the human genome. Mol Biol Evol. 2013;30:1544–58.

    Article  CAS  PubMed  Google Scholar 

  20. Chen H, Patterson N, Reich D. Population differentiation as a test for selective sweeps. Genome Res. 2010;20:393–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Voight BF, Kudaravalli S, Wen X, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006;4, e72.

    Article  PubMed  PubMed Central  Google Scholar 

  22. International HapMap C, Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, et al. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007;449:851–61.

    Article  CAS  Google Scholar 

  23. Kofler R, Schlotterer C. Gowinda: unbiased analysis of gene set enrichment for genome-wide association studies. Bioinformatics. 2012;28:2084–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. International HapMap. 2007. http://hapmap.ncbi.nlm.nih.gov/downloads/recombination/2008-03_rel22_B36/.

  25. Foll M, Gaggiotti O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics. 2008;180:977–93.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, Xie X, Byrne EH, McCarroll SA, Gaudet R, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007;449:913–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Vatsiou AI, Bazin E, and Gaggiotti OE. Detection of selective sweeps in structured populations: a comparison of recent methods. Molecular ecology. 2016;25(1):89-103.

  28. NCBI Genome Remapping Service. http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/genome/tools/remap.

  29. Lift Genome Annotation. https://genome.ucsc.edu/cgi-bin/hgLiftOver.

  30. Geer LY, Marchler-Bauer A, Geer RC, Han L, He J, He S, Liu C, Shi W, Bryant SH. The NCBI Biosystems database. Nucleic Acids Res. 2010;38:D492–6.

    Article  CAS  PubMed  Google Scholar 

  31. Tintle NL, Borchers B, Brown M, Bekmetjev A. Comparing gene set analysis methods on single-nucleotide polymorphism data from Genetic Analysis Workshop 16. BMC Proc. 2009;3 Suppl 7:S96.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Pareja-Tobes P, Tobes R, Manrique M, Pareja E, and Pareja-Tobes E. Bio4j: a high-performance cloud-enabled graph-based data platform. 2015. [Database website: http://bio4j.com/].

  33. Snel B, Lehmann G, Bork P, Huynen MA. STRING: a web-server to retrieve and display the repeatedly occurring neighbourhood of a gene. Nucleic Acids Res. 2000;28:3442–4. Database website: http://string-db.org/.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Parise LV. Integrin alpha(IIb)beta(3) signaling in platelet adhesion and aggregation. Curr Opin Cell Biol. 1999;11:597–601.

    Article  CAS  PubMed  Google Scholar 

  35. Almena M, Andrada E, Liebana R, Merida I. Diacylglycerol metabolism attenuates T-cell receptor signaling and alters thymocyte differentiation. Cell Death Dis. 2013;4, e912.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Datta S, Milner JD. Altered T-cell receptor signaling in the pathogenesis of allergic disease. J Allergy Clin Immunol. 2011;127:351–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Wong MM, Fish EN. Chemokines: attractive mediators of the immune response. Semin Immunol. 2003;15:5–14.

    Article  CAS  PubMed  Google Scholar 

  38. Xu T, Zhu Z, Sun Y, Ren L, Wang R. Characterization and expression of the CXCR1 and CXCR4 in miiuy croaker and evolutionary analysis shows the strong positive selection pressures imposed in mammal CXCR1. Dev Comp Immunol. 2014;44:133–44.

    Article  CAS  PubMed  Google Scholar 

  39. Bai L, Wang S. Targeting apoptosis pathways for new cancer therapeutics. Annu Rev Med. 2014;65:139–55.

    Article  CAS  PubMed  Google Scholar 

  40. Behl C. Apoptosis and Alzheimer’s disease. J Neural Transm. 2000;107:1325–44.

    Article  CAS  PubMed  Google Scholar 

  41. Munoz LE, van Bavel C, Franz S, Berden J, Herrmann M, van der Vlag J. Apoptosis in the pathogenesis of systemic lupus erythematosus. Lupus. 2008;17:371–5.

    Article  CAS  PubMed  Google Scholar 

  42. Pitrak DL. Apoptosis and Its Role in Neutrophil Dysfunction in AIDS. Oncologist. 1997;2:121–4.

    PubMed  Google Scholar 

  43. da Fonseca RR, Kosiol C, Vinar T, Siepel A, Nielsen R. Positive selection on apoptosis related genes. FEBS Lett. 2010;584:469–76.

    Article  PubMed  CAS  Google Scholar 

  44. World Health Organization. http://www.who.int/gho/hiv/en/.

  45. Rambaut A, Posada D, Crandall KA, Holmes EC. The causes and consequences of HIV evolution. Nat Rev Genet. 2004;5:52–61.

    Article  CAS  PubMed  Google Scholar 

  46. Sharp PM, Hahn BH. Origins of HIV and the AIDS pandemic. Cold Spring Harb Perspect Med. 2011;1:a006841.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Chen L, Perlina A, Lee CJ. Positive selection detection in 40,000 human immunodeficiency virus (HIV) type 1 sequences automatically identifies drug resistance and positive fitness mutations in HIV protease and reverse transcriptase. J Virol. 2004;78:3722–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. de Oliveira T, Salemi M, Gordon M, Vandamme AM, van Rensburg EJ, Engelbrecht S, Coovadia HM, Cassol S. Mapping sites of positive selection and amino acid diversification in the HIV genome: an alternative approach to vaccine design? Genetics. 2004;167:1047–58.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Snoeck J, Fellay J, Bartha I, Douek DC, Telenti A. Mapping of positive selection sites in the HIV-1 genome in the context of RNA and protein structural constraints. Retrovirology. 2011;8:87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Al-Sulami AA, Al-Edani TA, Al-Abdula AA. Culture Method and PCR for the Detection of Helicobacter pylori in Drinking Water in Basrah Governorate Iraq. Gastroenterol Res Pract. 2012;2012:245167.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Lyle CH, Annandale CH, Gouws J, Morley PS. Comparison of two culture techniques used to detect environmental contamination with Salmonella enterica in a large-animal hospital. J S Afr Vet Assoc. 2015;86:E1–5.

    Article  PubMed  Google Scholar 

  52. Qureshi T, Saeed A, Usmanghani K. Report: prevalence of shigellosis in three different areas of Karachi. Pak J Pharm Sci. 2014;27:381–8.

    PubMed  Google Scholar 

  53. Ribet D, Cossart P. How bacterial pathogens colonize their hosts and invade deeper tissues. Microbes Infect. 2015;17:173–83.

    Article  CAS  PubMed  Google Scholar 

  54. Lan R, Stevenson G, Reeves PR. Comparison of two major forms of the Shigella virulence plasmid pINV: positive selection is a major force driving the divergence. Infect Immun. 2003;71:6298–306.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Soto-Giron MJ, Ospina OE, Massey SE. Elevated levels of adaption in Helicobacter pylori genomes from Japan; a link to higher incidences of gastric cancer? Evol Med Public Health. 2015;2015:88–105.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Soyer Y, Orsi RH, Rodriguez-Rivera LD, Sun Q, Wiedmann M. Genome wide evolutionary analyses reveal serotype specific patterns of positive selection in selected Salmonella serotypes. BMC Evol Biol. 2009;9:264.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  57. Cossart P, Sansonetti PJ. Bacterial invasion: the paradigms of enteroinvasive pathogens. Science. 2004;304:242–8.

    Article  CAS  PubMed  Google Scholar 

  58. Chu YT, Wang YH, Wu JJ, Lei HY. Invasion and multiplication of Helicobacter pylori in gastric epithelial cells and implications for antibiotic resistance. Infect Immun. 2010;78:4157–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Kato A, Schleimer RP. Beyond inflammation: airway epithelial cells are at the interface of innate and adaptive immunity. Curr Opin Immunol. 2007;19:711–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Akira S, Uematsu S, Takeuchi O. Pathogen recognition and innate immunity. Cell. 2006;124:783–801.

    Article  CAS  PubMed  Google Scholar 

  61. Loo YM, Gale MJ. Immune signaling by RIG-I-like receptors. Immunity. 2011;34:680–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Areal H, Abrantes J, Esteves PJ. Signatures of positive selection in Toll-like receptor (TLR) genes in mammals. BMC Evol Biol. 2011;11:368.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Babik W, Dudek K, Fijarczyk A, Pabijan M, Stuglik M, Szkotak R, Zielinski P. Constraint and adaptation in newt toll-like receptor genes. Genome Biol Evol. 2015;7:81–95.

    Article  CAS  Google Scholar 

  64. Barreiro LB, Ben-Ali M, Quach H, Laval G, Patin E, Pickrell JK, Bouchier C, Tichit M, Neyrolles O, Gicquel B, et al. Evolutionary dynamics of human Toll-like receptors and their different contributions to host defense. PLoS Genet. 2009;5, e1000562.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  65. Darfour-Oduro KA, Megens HJ, Roca AL, Groenen MA, Schook LB. Adaptive Evolution of Toll-Like Receptors (TLRs) in the Family Suidae. PLoS One. 2015;10, e0124069.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Lemos de Matos A, McFadden G, Esteves PJ. Positive evolutionary selection on the RIG-I-like receptor genes in mammals. PLoS One. 2013;8:e81864.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  67. Al-Daghri NM, Cagliani R, Forni D, Alokail MS, Pozzoli U, Alkharfy KM, Sabico S, Clerici M, Sironi M. Mammalian NPC1 genes may undergo positive selection and human polymorphisms associate with type 2 diabetes. BMC Med. 2012;10:140.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Ayub Q, Moutsianas L, Chen Y, Panoutsopoulou K, Colonna V, Pagani L, Prokopenko I, Ritchie GR, Tyler-Smith C, McCarthy MI, et al. Revisiting the thrifty gene hypothesis via 65 loci associated with susceptibility to type 2 diabetes. Am J Hum Genet. 2014;94:176–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Pickrell JK, Coop G, Novembre J, Kudaravalli S, Li JZ, Absher D, Srinivasan BS, Barsh GS, Myers RM, Feldman MW, et al. Signals of recent positive selection in a worldwide sample of human populations. Genome Res. 2009;19:826–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Segurel L, Austerlitz F, Toupance B, Gautier M, Kelley JL, Pasquet P, Lonjou C, Georges M, Voisin S, Cruaud C, et al. Positive selection of protective variants for type 2 diabetes from the Neolithic onward: a case study in Central Asia. Eur J Hum Genet. 2013;21:1146–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Helgason A, Palsson S, Thorleifsson G, Grant SF, Emilsson V, Gunnarsdottir S, Adeyemo A, Chen Y, Chen G, Reynisdottir I, et al. Refining the impact of TCF7L2 gene variants on type 2 diabetes and adaptive evolution. Nat Genet. 2007;39:218–25.

    Article  CAS  PubMed  Google Scholar 

  72. Yoshiuchi I. Evidence of selection at insulin receptor substrate-1 gene loci. Acta Diabetol. 2013;50:775–9.

    Article  CAS  PubMed  Google Scholar 

  73. Guo X, Li H, Xu H, Woo S, Dong H, Lu F, Lange A, and Wu C. Glycolysis in the control of blood glucose homeostasis. 2012;2(4):358–67.

  74. Edgerton DS, Ramnanan CJ, Grueter CA, Johnson KM, Lautz M, Neal DW, Williams PE, Cherrington AD. Effects of insulin on the metabolic control of hepatic gluconeogenesis in vivo. Diabetes. 2009;58(12):2766–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Taylor R. Insulin resistance and type 2 diabetes. Diabetes. 2012;61(4):778–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Felig P, Wahren J, Hendler R, Brundin T. Splanchnic glucose and amino acid metabolism in obesity. J Clin Invest. 1974;53(2):582–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Gastaldelli A, Baldi S, Pettiti M, Toschi E, Camastra S, Natali A, Landau BR, Ferrannini E. Influence of obesity and type 2 diabetes on gluconeogenesis and glucose output in humans: a quantitative study. Diabetes. 2000;49(8):1367–73.

    Article  CAS  PubMed  Google Scholar 

  78. Foretz M, Hebrard S, Leclerc J, Zarrinpashneh E, Soty M, Mithieux G, Sakamoto K, Andreelli F, Viollet B. Metformin inhibits hepatic gluconeogenesis in mice independently of the LKB1/AMPK pathway via a decrease in hepatic energy state. J Clin Invest. 2010;120:2355–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Guoa X, Lia H, Xua H, Wooa S, Dongb H, Lub F, Langec A, Wu C. Glycolysis in the control of blood glucose homeostasis. Acta Pharm Sin B. 2012;2:358–67.

    Article  CAS  Google Scholar 

  80. Tousoulis D, Koniari K, Antoniades C, Miliou A, Noutsou M, Nikolopoulou A, Papageorgiou N, Marinou K, Stefanadi E, Stefanadis C. Impact of 6 weeks of treatment with low-dose metformin and atorvastatin on glucose-induced changes of endothelial function in adults with newly diagnosed type 2 diabetes mellitus: A single-blind study. Clin Ther. 2010;32:1720–8.

    Article  CAS  PubMed  Google Scholar 

  81. Current Worldwide Annual Meat consumption per capita. http://chartsbin.com/view/12730.

  82. Song Y, Manson JE, Buring JE, Liu S. A prospective study of red meat consumption and type 2 diabetes in middle-aged and elderly women: the women’s health study. Diabetes Care. 2004;27:2108–15.

    Article  CAS  PubMed  Google Scholar 

  83. Hosooka T, Noguchi T, Kotani K, Nakamura T, Sakaue H, Inoue H, Ogawa W, Tobimatsu K, Takazawa K, Sakai M, et al. Dok1 mediates high-fat diet-induced adipocyte hypertrophy and obesity through modulation of PPAR-gamma phosphorylation. Nat Med. 2008;14:188–93.

    Article  CAS  PubMed  Google Scholar 

  84. Kahn SE, Hull RL, Utzschneider KM. Mechanisms linking obesity to insulin resistance and type 2 diabetes. Nature. 2006;444:840–6.

    Article  CAS  PubMed  Google Scholar 

  85. Brady MJ. IRS2 takes center stage in the development of type 2 diabetes. J Clin Invest. 2004;114:886–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Kovacs P, Hanson RL, Lee YH, Yang X, Kobes S, Permana PA, Bogardus C, Baier LJ. The role of insulin receptor substrate-1 gene (IRS1) in type 2 diabetes in Pima Indians. Diabetes. 2003;52:3005–9.

    Article  CAS  PubMed  Google Scholar 

  87. Rung J, Cauchi S, Albrechtsen A, Shen L, Rocheleau G, Cavalcanti-Proenca C, Bacot F, Balkau B, Belisle A, Borch-Johnsen K, et al. Genetic variant near IRS1 is associated with type 2 diabetes, insulin resistance and hyperinsulinemia. Nat Genet. 2009;41:1110–5.

    Article  CAS  PubMed  Google Scholar 

  88. Wang T, Kusudo T, Takeuchi T, Yamashita Y, Kontani Y, Okamatsu Y, Saito M, Mori N, Yamashita H. Evodiamine inhibits insulin-stimulated mTOR-S6K activation and IRS1 serine phosphorylation in adipocytes and improves glucose tolerance in obese/diabetic mice. PLoS One. 2013;8, e83264.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  89. Hage M, Zantout MS, Azar ST. Thyroid disorders and diabetes mellitus. J Thyroid Res. 2011;2011:439463.

    Article  PubMed  PubMed Central  Google Scholar 

  90. Jin W, Xu S, Wang H, Yu Y, Shen Y, Wu B, Jin L. Genome-wide detection of natural selection in African Americans pre- and post-admixture. Genome Res. 2012;22:519–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Arbiza L, Dopazo J, Dopazo H. Positive selection, relaxation, and acceleration in the evolution of the human and chimp genome. PLoS Comput Biol. 2006;2, e38.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  92. Catlett IM, Hedrick SM. Suppressor of cytokine signaling 1 is required for the differentiation of CD4+ T cells. Nat Immunol. 2005;6:715–21.

    Article  CAS  PubMed  Google Scholar 

  93. Hancock AM, Clark VJ, Qian Y, Di Rienzo A. Population genetic analysis of the uncoupling proteins supports a role for UCP3 in human cold resistance. Mol Biol Evol. 2011;28:601–14.

    Article  CAS  PubMed  Google Scholar 

  94. Hesselson SE, Matsson P, Shima JE, Fukushima H, Yee SW, Kobayashi Y, Gow JM, Ha C, Ma B, Poon A, et al. Genetic variation in the proximal promoter of ABC and SLC superfamilies: liver and kidney specific expression and promoter activity predict variation. PLoS One. 2009;4, e6942.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  95. Laland KN, Odling-Smee J, Myles S. How culture shaped the human genome: bringing genetics and the human sciences together. Nat Rev Genet. 2010;11:137–48.

    Article  CAS  PubMed  Google Scholar 

  96. Lao O, de Gruijter JM, van Duijn K, Navarro A, Kayser M. Signatures of positive selection in genes associated with human skin pigmentation as revealed from analyses of single nucleotide polymorphisms. Ann Hum Genet. 2007;71:354–69.

    Article  CAS  PubMed  Google Scholar 

  97. Li DF, Liu WB, Liu JF, Yi GQ, Lian L, Qu LJ, Li JY, Xu GY, Yang N. Whole-genome scan for signatures of recent selection reveals loci associated with important traits in White Leghorn chickens. Poult Sci. 2012;91:1804–12.

    Article  CAS  PubMed  Google Scholar 

  98. Luisi P, Alvarez-Ponce D, Dall’Olio GM, Sikora M, Bertranpetit J, Laayouni H. Network-level and population genetics analysis of the insulin/TOR signal transduction pathway across human populations. Mol Biol Evol. 2012;29:1379–92.

    Article  CAS  PubMed  Google Scholar 

  99. Morgan CC, Shakya K, Webb A, Walsh TA, Lynch M, Loscher CE, Ruskin HJ, O’Connell MJ. Colon cancer associated genes exhibit signatures of positive selection at functionally significant positions. BMC Evol Biol. 2012;12:114.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Patrick MS, Oda H, Hayakawa K, Sato Y, Eshima K, Kirikae T, Iemura S, Shirai M, Abe T, Natsume T, et al. Gasp, a Grb2-associating protein, is critical for positive selection of thymocytes. Proc Natl Acad Sci U S A. 2009;106:16345–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  101. Yu XJ, Zheng HK, Wang J, Wang W, Su B. Detecting lineage-specific adaptive evolution of brain-expressed genes in human using rhesus macaque as outgroup. Genomics. 2006;88:745–51.

    Article  CAS  PubMed  Google Scholar 

  102. Inagaki T, Tachibana M, Magoori K, Kudo H, Tanaka T, Okamura M, Naito M, Kodama T, Shinkai Y, Sakai J. Obesity and metabolic syndrome in histone demethylase JHDM2a-deficient mice. Genes Cells. 2009;14:991–1001.

    Article  CAS  PubMed  Google Scholar 

  103. Tateishi K, Okada Y, Kallin EM, Zhang Y. Role of Jhdm2a in regulating metabolic gene expression and obesity resistance. Nature. 2009;458:757–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  104. Clement K, Le Stunff C, Meirhaeghe A, Dechartres A, Ferrieres J, Basdevant A, Boitard C, Amouyel P, Bougneres P. In obese and non-obese adults, the cis-regulatory rs361072 promoter variant of PIK3CB is associated with insulin resistance not with type 2 diabetes. Mol Genet Metab. 2009;96:129–32.

    Article  CAS  PubMed  Google Scholar 

  105. Gloyn AL, Braun M, Rorsman P. Type 2 diabetes susceptibility gene TCF7L2 and its role in beta-cell function. Diabetes. 2009;58:800–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  106. Balaresque PL, Ballereau SJ, Jobling MA. Challenges in human genetic diversity: demographic history and adaptation. Hum Mol Genet. 2007;16(Spec No 2):R134–9.

    Article  CAS  PubMed  Google Scholar 

  107. Comas I, Coscolla M, Luo T, Borrell S, Holt KE, Kato-Maeda M, Parkhill J, Malla B, Berg S, Thwaites G, et al. Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans. Nat Genet. 2013;45:1176–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  108. Cagliani R, Sironi M. Pathogen-driven selection in the human genome. Int J Evol Biol. 2013;2013:204240.

    Article  PubMed  PubMed Central  Google Scholar 

  109. Jackson BC, Campos JL, and Zeng K. The effects of purifying selection on patterns of genetic differentiation between Drosophila melanogaster populations. Heredity. 2015;114(2):163-74.

  110. Casillas S, Barbadilla A, Bergman CM. Purifying selection maintains highly conserved noncoding sequences in Drosophila. Mol Biol Evol. 2007;24:2222–34.

    Article  CAS  PubMed  Google Scholar 

  111. England PH. Adult obesity and type 2 diabetes. 2014

  112. Watve MG, Yajnik CS. Evolutionary origins of insulin resistance: a behavioral switch hypothesis. BMC Evol Biol. 2007;7:61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by the Marie-Curie Initial Training Network INTERCROSSING (European Commission FP7). OEG was further supported by the MASTS pooling initiative (The Marine Alliance for Science and Technology for Scotland). The authors thank Raquel Tobes for extracting the genes from the Bio4j database.

Availability of data and materials

All the genomic data used for the current analysis are freely available to be downloaded from the HapMap Project. The tool used for the gene set enrichment analysis is freely available in https://github.com/INTERCROSSING/SEL-GSEA. Additional methodological details and results are available as Additional files 1 and 2 an in the end of the main manuscript. Further information is available under request.

Authors’ contributions

OEG designed the study. AIV wrote the scripts and carried out the analysis. OEG has supervised the whole implementation and paper writing. All authors read and approved the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

  1. Laboratoire d’Écologie Alpine (LECA), Univesrity Joseph Fourier, 2233 Rue de la Piscine, 38041, Grenoble, Cedex 9, France

    Alexandra I. Vatsiou, Eric Bazin & Oscar E. Gaggiotti

  2. Scottish Oceans Institute, East Sands, University of St Andrews, St Andrews, KY16 8LB, Scotland, UK

    Alexandra I. Vatsiou & Oscar E. Gaggiotti

  3. Oh no sequences! Research group, Era7Bioinformatics, Plaza de Campo Verde, 3, 18001, Granada, Spain

    Alexandra I. Vatsiou

Authors
  1. Alexandra I. Vatsiou
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Eric Bazin
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Oscar E. Gaggiotti
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Alexandra I. Vatsiou.

Additional files

Additional file 1:

Supplementary Material. (DOCX 1057 kb)

Additional file 2: Table S1.

Significant gene sets found using iHS scores and Daub et al. [19] approach. Threshold is set to <0.09. Table S2. Significant gene sets found using iHS scores and Gowinda approach. Threshold is set to <0.09. Table S3. Significant gene sets found using XPCLR scores and Daub et al. [19] approach. Threshold is set to <0.09. For all the population pairs, the first population is the objective one and the second, the reference. Table S4. Significant gene sets found using XPCLR scores and Gowinda approach. Threshold is set to 0.09. For all the population pairs, the first population is the objective one and the second, the reference. Table S5. Immunity related gene sets detected with the GSEA approaches. Q-value threshold is set to <0.09. Table S6. 17 significant genes related to obesity, diabetes and metabolic syndrome that were found to be under positive selection with XPCLR and iHS analysis using the list of genes derived from Bio4j. Some of them (indicated with *) have been detected in previous studies to be under positive selection, too. The threshold was calculated based on the 1 % cut off level. Genes are categorized in groups of potential risk factors, potential protective and indirect associations. Table S7. Significant genes related to obesity, diabetes and metabolic syndrome that we found to be under positive selection from XPCLR and iHS analysis using the Protein-Protein Interaction (PPI) networks from the STRING database. Five of them (indicated with *) have been detected in previous studies to be under positive selection. (XLSX 29 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Vatsiou, A.I., Bazin, E. & Gaggiotti, O.E. Changes in selective pressures associated with human population expansion may explain metabolic and immune related pathways enriched for signatures of positive selection. BMC Genomics 17, 504 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-016-2783-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-016-2783-2

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us