Human history has been shaped by infectious diseases. Human genes, especially host defense genes, have been constantly influenced by the pathogens encountered (Fumagalli and Sironi, 2014; Karlsson et al., 2014; Quintana-Murci and Clark, 2013). Pathogens drive the selection of genetic variants affecting resistance or tolerance to the infection, and heritable variations that increase survival to diseases with high morbidity and mortality will be naturally selected in people before reproductive age (Karlsson et al., 2014). These selection signatures vary with historical period, virulence of the pathogen, and the geographical spread.

Here we investigated the historical evolutionary patterns leading to genetic adaptation in cytokine production and immune-mediated diseases, including infectious, inflammatory, and autoimmune diseases. Cytokine production capacity is a key component of the host defense mechanisms: it induces inflammation, activates phagocytes to eliminate the pathogens and present antigens, and controls induction of T-helper (Th) adaptive immune responses. We have therefore chosen to investigate the evolutionary trajectories of cytokine production capacity in modern human populations during history. To determine the difference in polygenic regulation of diseases and cytokine production capacity, we used data derived from the 500 Functional Genomics (500FG) cohort of the Human Functional Genomics Project (HFGP; http://www.humanfunctionalgenomics.org). The HFGP is an international collaboration aiming to identify the host and environmental factors responsible for the variability of human immune responses in health and disease (Netea et al., 2016). Within the HFGP project, the 500FG study generated a large database of immunological, phenotypic, and multi-omics data from a cohort of 534 individuals of Western-European ancestry, which has been used to integrate the impact of genetic and environmental factors on cytokine production and immune parameters. We subsequently deciphered the factors that influence inter-individual variations in the immune responses against different stimuli (Bakker et al., 2018; Li et al., 2016; Schirmer et al., 2016; Ter Horst et al., 2016).

Peripheral blood mononuclear cells from these individuals were challenged with bacterial, fungal, viral, and non-microbial stimuli, and six cytokines (tumor necrosis factor (TNF)α, interleukin (IL)-1β, IL-6, IL-17, IL-22, and interferon (IFN)γ) were measured at 24 hr or 7 days after stimulation, generating 105 cytokine-stimulation pairs (Figure 1—figure supplement 1 and Supplementary file 1a). The stimulation time intervals were chosen based on extensive studies that showed that the time points used were best suited for assessing monocyte-derived and lymphocyte-derived cytokines per stimulus. Not all the stimuli induce the production of all cytokines; so the selection of the cytokine-stimulus pairs was performed for those pairs for which cytokine production was measurable (Li et al., 2016; Ruschen et al., 1992; Schirmer et al., 2016; Ter Horst et al., 2016; van de Veerdonk et al., 2009). We correlated cytokine production with genetic variant data to obtain cytokine quantitative trait loci (QTLs), which were employed to compute and compare the polygenic risk score (PRS) of the genomes of 827 individuals from different human historical eras (early upper Paleolithic, late upper Paleolithic, Mesolithic, Neolithic, post-Neolithic), which were downloaded from version 37.2 of the compiled dataset containing unimputed published ancient genotypes (https://reich.hms.harvard.edu/downloadable-genotypes-present-day-and-ancient-dna-data-compiled-published-papers) and 250 modern Europeans randomly selected from the European 1000G cohort. The individuals in this cohort present a similar genetic background (Western-European ancestry) and balanced characteristics in terms of age, sex, body mass index (BMI), and habits (Ter Horst et al., 2016), allowing us to represent the natural variability in the immune responses and minimizing the impact of the inter-individual variability in our results. In line with this, cytokines were measured in batches in which a certain cytokine was measured the same day for all the samples, in order to decrease the potential influence of technical variations. Intra-individual variation of cytokine production capacity was found to be limited (Ter Horst et al., 2021,) while inter-individual variability was largely depedent on genetic variants (Bakker et al., 2018; Li et al., 2016).

We investigated the PRS changes over time through constructing linear models and correlation analysis. In order to account for the ancient DNA (aDNA) samples being pseudo-haploid, ambiguous single-nucleotide polymorphisms (SNPs) (A/T and C/G) were excluded when computing PRS to prevent errors due to strand flips. PRS was computed using the most significant QTLs that had a p-value lower than our predetermined threshold for each given trait and removing all variants within a 250-kb window around these variants. Additional PRS models were calculated at varying window sizes (250, 500, 1000 kb) in order to show consistency in the direction of the trends (Figure 1—figure supplement 2). The dosage of these variants was multiplied by their effect size while the dosage of missing variants in a sample was supplemented with the average dosage. Although replacing missing values with the population average does not skew the observed trends over time, it does affect the scale of the data. We opted for this approach, however, due to the scarcity of variants shared by all samples; we focused instead on consistent observations after excluding samples with higher missing genotype rates. Finally, we scaled the PRS to a range of −1 and 1 and correlated the scores of the samples with their respective carbon-dated ages. In order to verify the robustness of our results, we repeated the analysis at multiple threshold combinations for variant missingness and QTL thresholds.

We plotted locally estimated scatterplot smoothing (LOESS) regression models to show the variations in PRS across time for each trait and highlight the time periods around which they drastically change. After plotting the changes before and after the Neolithic revolution at various thresholds for both the missing genotype rate in samples as well as QTL inclusion threshold, only traits that showed consistently significant trends in the same direction were chosen. Additional boxplots and t-tests show the variations in PRS between different pairings besides pre- and post-Neolithic samples. A schematic representation of the steps performed is shown in Figure 1—figure supplement 3. Lastly, we used a separate test to determine whether traits were under selective pressure between any two broad time periods using the trait-associated Wright’s fixation index (Fst). This test shows an increase in the number of traits under selective pressure after the start of the Neolithic and highlights the traits and points in time for which the observed trends in PRS can be partially attributed to selective pressure as opposed to drift (Supplementary file 1b).

Applying the methodology described above, several patterns were apparent (Figure 1). The first overall observation is that the estimation of cytokine production capacity based on PRS shows significant differences between populations in various historical periods, and the strength of evolutionary pressure on cytokine responses was different before and after the Neolithic revolution. We did not observe significant changes in cytokine production capacity between individuals who lived in different historical periods before the Neolithic, whereas strong pressure is apparent after adoption of agriculture and animal domestication in Europe. This different pattern may have resulted from the more limited number of samples available for the older time periods, resulting in lower statistical power, but the presence of some evolutionary pressure also before the Neolithic argues that this is most likely not the full explanation. The development of agriculture and domestication of animals in the Neolithic increased population densities on the one hand and the contact between humans and domesticated animals as a source of pathogens on the other hand. The number of zoonoses increased dramatically (examples being tuberculosis, brucellosis, Q-fever, and influenza), which strongly increased the selective pressure and caused significant adaptations of immunity at the genetic level (Flandroy et al., 2018). Most of the genetic adaptations to pathogens took place in the period since modern humans abandoned their hunting-gathering lifestyle and developed agriculture (Deschamps et al., 2016). In this respect, the strongest changes leading to tolerance (decreased cytokine production) were exerted in the cytokine responses to intracellular zoonotic infections (tuberculosis and Coxiella). In contrast, responses to the extracellular pathogens Staphylococcus aureus and Candida albicans indicate increased resistance, with high production of IL-22 and TNFα, respectively. The increased response to the important fungal pathogen C. albicans after the Neolithic period is validated also at the transcriptional level. Overall, these patterns are reminiscent of the studies showing that human immune responses need to adapt to a new landscape of infectious agents depending on the geographical location and types of microbe encountered (Ferwerda et al., 2007). Such different patterns were most likely encountered also through history.

Figure 1 with 3 supplements see all

Download asset Open asset

Importantly, our results also show significant patterns in changes in the production of specific cytokines during history. The resistance against intracellular pathogens increased after Neolithic with higher IFNγ responses (see Figure 1—figure supplement 1): indeed, it is known that Th1-IFNγ responses are crucial for the host defense against intracellular pathogens such as mycobacteria or Coxiella (Thakur et al., 2019). In addition, the resistance to the extracellular pathogens C. albicans and S. aureus is also increased after this Neolithic era, with TNFα and IFNγ production increasing steadily after. These two cytokines are very well known to be important for anti-Candida and anti-Staphylococcus host defense (MRSA Systems Immunobiology Group et al., 2018; Domínguez-Andrés et al., 2017). On the other hand, a different pattern emerges in relation with the IL1β-IL6-IL17 axis: the production of these cytokines is seen decreasing after Neolithic (see Figure 1a and b). In this context, the decrease through time of poly I:C induction of cytokines, as a model of viral stimulation, is intriguing but potentially very important: many important viruses such as influenza and coronaviruses (severe acute respiratory syndrome (SARS), Middle East respiratory syndrome (MERS), and SARS coronavirus 2 (SARS-CoV-2)) exert life-threatening effects through induction of cytokine-mediated hyperinflammation (also termed ‘cytokine storm’) (Tay et al., 2020); evolutionary processes to curtail these exaggerated responses are thus likely to be protective, and tolerance against viruses becomes a host defense mechanism (Diard and Hardt, 2017).

These evolutionary genetic adaptations to pathogens throughout human history greatly influence the way we respond to multiple diseases in modern times as well. To assess these effects, we calculated the PRS associated with the risk of several highly prevalent immune-mediated diseases. The first focus was on common infectious diseases such as malaria, human immunodeficiency virus-acquired immunodeficiency syndrome (HIV-AIDS), tuberculosis, and chronic viral hepatitis; we calculated the changes in susceptibility to these diseases in the last 50,000 years of human history, based on summary statistics from genome-wide association studies (GWAS) databases available from the literature (Figure 2—figure supplement 1). Our results show that humans are becoming more resistant to these diseases, with the notable exception of tuberculosis, whose risk score remained stable along the period studied (Figure 2). Of note, the QTLs that passed our thresholds were scarce, resulting in a PRS model using a limited number of SNPs, making them sensitive to changes in the p-value cutoff, especially in traits related with infectious diseases. Our results suggest that humans have built up a genetic makeup that made them more resistant to a variety of microbes. The pattern of this adaptation is very interesting as well, with a suggested decrease of susceptibility to malaria especially in the last 10,000 years. The reason for this accelerated resistance after Neolithic might be linked to a higher disease prevalence due to increased population density, as otherwise Plasmodium parasites are known to have circulated in Africa since at least the Paleogene 30 million years ago (Poinar, 2005), and we have likely inherited it from gorillas (Liu et al., 2010). Intriguingly, we also observed a strong decrease in susceptibility to HIV: this is a contemporary pathogen, therefore this signal could be due to common genetic and immune pathways with other infections that were present in human populations. The increased resistance to HIV in Europeans may be derived from selective pressures induced by other pathogens such as Yersinia pestis (Duncan et al., 2005). Our data suggest, on the other hand, that the source of this increased resistance is even older.

Figure 2 with 1 supplement see all

Download asset Open asset

In contrast, the lack of genetic adaptation in susceptibility to tuberculosis is intriguing. This surprising finding may be explained by a concept in which Mycobacterium tuberculosis is at the same time a pathogen and a symbiont, in which latent infection enhances the resistance against other pathogens, and this is why our immune system tolerates mycobacterial presence (Pai et al., 2016). In this regard, individuals with latent tuberculosis exhibit enhanced macrophage functions that may protect against other pathogens through the induction of trained immunity (Joosten et al., 2018). In this context, humanity may not be adapting to tuberculosis because increased resistance against mycobacteria is not evolutionarily advantageous. All in all, these results suggest that the risk of suffering infectious diseases has steadily decreased at least for the last 50,000 years as a result of the selection of genetic variants that confer resistance to infections.

It has been proposed that the increased prevalence of inflammatory and autoimmune diseases is associated with the immune-related alleles that have been positively selected through evolutionary processes to protect against infection; hence, the contrasting differences in the prevalence of autoimmune diseases between populations result from diverse selective pressures (Ramos et al., 2015). In line with this, it has been hypothesized that genetic variants associated with protection against infectious agents are behind the increased prevalence of autoimmune diseases in populations with low pathogen exposure, such as Europeans (Fumagalli et al., 2011; Raj et al., 2013). To study the changing patterns of susceptibility to autoimmune and inflammatory diseases during history, we used publicly available summary statistics from GWAS of digestive tract-related autoimmune and inflammatory diseases and arthritis-related diseases (see Figure 2—figure supplement 1) and calculated the PRS for each of the samples under study. Interestingly, we observed a robust increase of the genetic variants related with the development of inflammatory diseases in the digestive tract after the Neolithic revolution (Figure 3). PRS scores associated with celiac disease, Crohn’s disease, ulcerative colitis, and inflammatory bowel disease were strongly associated with the age of the samples, regardless of the p-value thresholds or the missing genotype rates used for PRS calculation, showing the robustness of these results (see Figure 1—figure supplement 2). The fact that especially intestinal inflammatory pathology is increased after a historical event that fundamentally modified human diet is unlikely to be an accident. Our results are in line with earlier research demonstrating that variants in genes important for immune responses and involved in celiac disease pathophysiology (such as IL-12, IL-18RAP, SH2B3) are under strong positive selection (Zhernakova et al., 2010). The reasons for the selection pressure on these genes are not completely understood, but an advantage for host defense has been suggested (Zhernakova et al., 2010).

Figure 3

Download asset Open asset

In contrast to intestinal inflammation, the PRS of traits linked with juvenile-idiopathic arthritis, rheumatoid arthritis, and multiple sclerosis shows a decrease in genetic susceptibility with the age of the sample after the Neolithic revolution. For pre-Neolithic periods, these patterns had little impact with decreasing PRS for digestive tract diseases and increasing PRS for ankylosing spondylitis and juvenile idiopathic arthritis. A strong decrease in susceptibility to juvenile idiopathic arthritis, rheumatoid arthritis, and multiple sclerosis is seen after the Neolithic period (see Figure 3). This is likely linked to the decreased production of the IL-1/IL-6/IL-17 axis described in Figure 2, which is particularly important in the pathophysiology of these disorders (Akioka, 2019; Mei et al., 2011).

The significant changes in cytokine production and disease susceptibility in European populations after the Neolithic can be due to selective processes on the one hand (as described above), but also due to important demographic changes due to migrations of human communities such as the Anatolians (in Neolithic) or the Yamnaya populations from the Pontic steppe (during the Bronze Age) (Racimo et al., 2020). In this regard, several loci associated with inflammatory diseases displayed a group of alleles linked with Crohn’s disease, celiac disease, and ulcerative colitis in Neolithic Aegeans, the community that spread farming across Europe (Hofmanová et al., 2016), with several of these alleles showing signs of positive selection in modern Europeans (Raj et al., 2013). In addition, the gene expression PRS of several cytokines based on the cis- and trans- [gene] expression quantitative trait locus (eQTLs) from the eQTLGen Consortium (https://www.eqtlgen.org/) displayed a very strong association with time for TNFα after the Neolithic revolution (Figure 4).

Figure 4 with 1 supplement see all

Download asset Open asset

Of note, the availability of samples in the pre-Neolithic dataset is limited compared to the post-Neolithic era. To test the robustness of the trends observed, we recalculated the correlation coefficients and p-values to show the consistency of the results from the down-sampled data of the post-Neolithic samples to match the sample size of the pre-Neolithic samples. This analysis showed that the trajectories observed in our results are consistent regardless of the sample size (Figure 4—figure supplement 1).

In this study we focus on the effects of different pathogens on cytokine production and, through that, on the response to different infections. Human evolution has been influenced by multiple factors such as migration, urbanization, climate, and diet, which also influence the responses to pathogens. The advent of the Neolithic lifestyle was a milestone for human societies: after the Neolithic, the density of human populations and the rate of contact with domesticated animals increased significantly, which increased the emergence of pathogens and the ease with which those could spread. Our data show that while various events and conditions throughout human history have influenced our immune responses to pathogens, the advent of the Neolithic lifestyle was an important turning point for our capacity to respond to pathogens.

Collectively, our results show that the advent of the Neolithic era was a turning point for the evolution of immune-mediated traits in European populations, driving the expansion of alleles that favor the development of tolerance against intracellular pathogens and promote inflammatory responses against extracellular microbes. This is associated with a higher presence of genetic traits related with inflammatory and auto-immune diseases of the digestive tract and a lower number of alleles linked with the development of arthritis. It is important to underline that our results are obtained using European aDNA samples and GWAS summary statistics computed using European cohorts. Therefore, our results would be specific to European populations and should be interpreted with caution for populations from other geographical locations. Further research should compare the trends in different populations that have been exposed to different environments across the planet and clarify the influence of ancestry, time, and rural vs urban lifestyle to shed light on the influence of the infectious environment on genetics and human evolution.

All the data employed in this manuscript have been obtained from publicly available databases.

1. 1. Trynka G
   2. Hunt KA
   3. Bockett NA
   4. Romanos J
   5. Mistry V
   6. Szperl A
   7. Bakker SF
   8. Bardella MT
   9. Bhaw-Rosun L
   10. Castillejo G
   11. Almeida RC
   12. Dias KR
   13. Dubois PC
   14. Duerr RH
   15. Edkins S
   16. Franke L
   17. Fransen K
   18. Gutierrez J
   19. Heap GA
   20. Hrdlickova B
   21. Hunt S
   22. Izzo V
   23. Joosten LA
   24. Langford C
   25. Mazzilli MC
   26. Mein CA
   27. Midah V
   28. Mitrovic M
   29. Mora B

   (2011) NHGRI-EBI GWAS Catalog

   ID 22057235. Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease.

   https://www.ebi.ac.uk/gwas/publications/22057235
2. 1. Liu JZ
   2. Huang H
   3. Ng SC
   4. Alberts R
   5. Takahashi A
   6. Ripke S
   7. Lee JC
   8. Jostins L
   9. Shah T
   10. Abedian S
   11. Cheon JH
   12. Cho J
   13. Dayani NE
   14. Franke L
   15. Fuyuno Y
   16. Hart A
   17. Juyal RC
   18. Juyal G

   (2015) NHGRI-EBI GWAS Catalog

   ID 26192919. Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations.

   https://www.ebi.ac.uk/gwas/publications/26192919
3. 1. Cortes A
   2. Hadler J
   3. Pointon JP
   4. Robinson PC
   5. Karaderi T
   6. Leo P
   7. Cremin K
   8. Pryce K
   9. Harris J
   10. Lee S
   11. Joo KB
   12. Shim SC
   13. Weisman M
   14. Ward M
   15. Zhou X
   16. Garchon HJ
   17. Chiocchia G
   18. Nossent J
   19. Lie BA
   20. Førre Ø
   21. Tuomilehto J
   22. Laiho K
   23. Jiang L
   24. Liu Y
   25. Wu X
   26. Bradbury LA
   27. Elewaut D
   28. Burgos-Vargas R
   29. Stebbings S
   30. Appleton L
   31. Farrah C
   32. Lau J
   33. Kenna TJ
   34. Haroon N
   35. Ferreira MA
   36. Yang J
   37. Mulero J
   38. Fernandez-Sueiro JL
   39. Gonzalez-Gay MA
   40. Lopez-Larrea C
   41. Deloukas P
   42. Donnelly P

   (2013) NHGRI-EBI GWAS Catalog

   ID 23749187. Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci".

   https://www.ebi.ac.uk/gwas/publications/23749187
4. 1. Hinks A
   2. Cobb J
   3. Marion MC
   4. Prahalad S
   5. Sudman M
   6. Bowes J
   7. Martin P
   8. Comeau ME
   9. Sajuthi S
   10. Andrews R
   11. Brown M
   12. Chen WM
   13. Concannon P
   14. Deloukas P
   15. Edkins S
   16. Eyre S
   17. Gaffney PM
   18. Guthery SL
   19. Guthridge JM
   20. Hunt SE
   21. James JA
   22. Keddache M
   23. Moser KL
   24. Nigrovic PA
   25. Onengut-Gumuscu S
   26. Onslow ML
   27. Rosé CD
   28. Rich SS
   29. Steel KJ
   30. Wakeland EK
   31. Wallace CA
   32. Wedderburn LR
   33. Woo P
   34. Bohnsack JF
   35. Haas JP
   36. Glass DN
   37. Langefeld CD
   38. Thomson W
   39. Thompson SD

   (2013) NHGRI-EBI GWAS Catalog

   ID 23603761. Dense genotyping of immune-related disease regions identifies 14 new susceptibility loci for juvenile idiopathic arthritis.

   https://www.ebi.ac.uk/gwas/publications/23603761
5. 1. Okada Y
   2. Wu D
   3. Trynka G
   4. Raj T
   5. Terao C
   6. Ikari K
   7. Kochi Y
   8. Ohmura K
   9. Suzuki A
   10. Yoshida S
   11. Graham RR
   12. Manoharan A
   13. Ortmann W
   14. Bhangale T
   15. Denny JC
   16. Carroll RJ
   17. Eyler AE
   18. Greenberg JD
   19. Kremer JM
   20. Pappas DA
   21. Jiang L
   22. Yin J
   23. Ye L
   24. Su DF
   25. Yang J
   26. Xie G
   27. Keystone E
   28. Westra HJ
   29. Esko T
   30. Metspalu A
   31. Zhou X
   32. Gupta N
   33. Mirel D
   34. Stahl EA
   35. Diogo D
   36. Cui J
   37. Liao K
   38. Guo MH
   39. Myouzen K
   40. Kawaguchi T
   41. Coenen MJ

   (2013) NHGRI-EBI GWAS Catalog

   ID 24390342. Genetics of rheumatoid arthritis contributes to biology and drug discovery.

   https://www.ebi.ac.uk/gwas/publications/24390342
6. 1. Beecham AH
   2. Patsopoulos NA
   3. Xifara DK
   4. Davis MF
   5. Kemppinen A
   6. Cotsapas C
   7. Shah TS
   8. Spencer C
   9. Booth D
   10. Goris A
   11. Oturai A
   12. Saarela J
   13. Fontaine B
   14. Hemmer B
   15. Martin C
   16. Zipp F
   17. D'Alfonso S
   18. Martinelli-Boneschi F
   19. Taylor B
   20. Harbo HF
   21. Kockum I
   22. Hillert J
   23. Olsson T
   24. Ban M
   25. Oksenberg JR
   26. Hintzen R
   27. Barcellos LF
   28. Agliardi C
   29. Alfredsson L
   30. Alizadeh M
   31. Anderson C
   32. Andrews R
   33. Søndergaard HB
   34. Baker A
   35. Band G
   36. Baranzini SE
   37. Barizzone N
   38. Barrett J
   39. Bellenguez C
   40. Bergamaschi L
   41. Bernardinelli L
   42. Berthele A
   43. Biberacher V
   44. Binder TM
   45. Blackburn H
   46. Bomfim IL
   47. Brambilla P
   48. Broadley S
   49. Brochet B
   50. Brundin L
   51. Buck D
   52. Butzkueven H
   53. Caillier SJ
   54. Camu W
   55. Carpentier W
   56. Cavalla P
   57. Celius EG
   58. Coman I
   59. Comi G
   60. Corrado L
   61. Cosemans L
   62. Cournu-Rebeix I
   63. Cree BA
   64. Cusi D
   65. Damotte V
   66. Defer G

   (2013) NHGRI-EBI GWAS Catalog

   ID 24076602. Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis.

   https://www.ebi.ac.uk/gwas/publications/24076602
7. 1. Beecham AH
   2. Patsopoulos NA

   (2019) David Reich Lab published data amh_repo V37.2

   ID V37.2. Downloadable genotypes of present-day and ancient DNA data (compiled from published papers) V37.2.

   https://reich.hms.harvard.edu/allen-ancient-dna-resource-aadr-downloadable-genotypes-present-day-and-ancient-dna-data
8. 1. Bakker OB

   (2018) UK Biobank

   ID uk-biobank. GWAS round 2.

   http://www.nealelab.is/uk-biobank
9. 1. Bakker OB

   (2019) The 1000ibd MOLGENIS database

   ID 500FG/cytokines. 500FG Cytokine QTL data.

   https://molgenis26.gcc.rug.nl/downloads/500FG/cytokines.zip

1. 1. Akioka S

   (2019) Interleukin-6 in juvenile idiopathic arthritis

   Modern Rheumatology 29:275–286.

   https://doi.org/10.1080/14397595.2019.1574697

   • PubMed
   • Google Scholar
2. 1. Bakker OB
   2. Aguirre-Gamboa R
   3. Sanna S
   4. Oosting M
   5. Smeekens SP
   6. Jaeger M
   7. Zorro M
   8. Võsa U
   9. Withoff S
   10. Netea-Maier RT
   11. Koenen H
   12. Joosten I
   13. Xavier RJ
   14. Franke L
   15. Joosten LAB
   16. Kumar V
   17. Wijmenga C
   18. Netea MG
   19. Li Y

   (2018) Integration of multi-omics data and deep phenotyping enables prediction of cytokine responses

   Nature Immunology 19:776–786.

   https://doi.org/10.1038/s41590-018-0121-3

   • PubMed
   • Google Scholar
3. 1. Bycroft C
   2. Freeman C
   3. Petkova D
   4. Band G
   5. Elliott LT
   6. Sharp K
   7. Motyer A
   8. Vukcevic D
   9. Delaneau O
   10. O'Connell J
   11. Cortes A
   12. Welsh S
   13. Young A
   14. Effingham M
   15. McVean G
   16. Leslie S
   17. Allen N
   18. Donnelly P
   19. Marchini J

   (2018) The UK biobank resource with deep phenotyping and genomic data

   Nature 562:203–209.

   https://doi.org/10.1038/s41586-018-0579-z

   • PubMed
   • Google Scholar
4. 1. Deschamps M
   2. Laval G
   3. Fagny M
   4. Itan Y
   5. Abel L
   6. Casanova JL
   7. Patin E
   8. Quintana-Murci L

   (2016) Genomic signatures of selective pressures and introgression from archaic hominins at human innate immunity genes

   The American Journal of Human Genetics 98:5–21.

   https://doi.org/10.1016/j.ajhg.2015.11.014

   • PubMed
   • Google Scholar
5. 1. Diard M
   2. Hardt WD

   (2017) Evolution of bacterial virulence

   FEMS Microbiology Reviews 41:679–697.

   https://doi.org/10.1093/femsre/fux023

   • PubMed
   • Google Scholar
6. 1. Domínguez-Andrés J
   2. Feo-Lucas L
   3. Minguito de la Escalera M
   4. González L
   5. López-Bravo M
   6. Ardavín C

   (2017) Inflammatory Ly6C high monocytes protect against candidiasis through IL-15-Driven NK cell/Neutrophil activation

   Immunity 46:1059–1072.

   https://doi.org/10.1016/j.immuni.2017.05.009

   • PubMed
   • Google Scholar
7. 1. Duncan SR
   2. Scott S
   3. Duncan CJ

   (2005) Reappraisal of the historical selective pressures for the CCR5-Delta32 mutation

   Journal of Medical Genetics 42:205–208.

   https://doi.org/10.1136/jmg.2004.025346

   • PubMed
   • Google Scholar
8. 1. Ferwerda B
   2. McCall MB
   3. Alonso S
   4. Giamarellos-Bourboulis EJ
   5. Mouktaroudi M
   6. Izagirre N
   7. Syafruddin D
   8. Kibiki G
   9. Cristea T
   10. Hijmans A
   11. Hamann L
   12. Israel S
   13. ElGhazali G
   14. Troye-Blomberg M
   15. Kumpf O
   16. Maiga B
   17. Dolo A
   18. Doumbo O
   19. Hermsen CC
   20. Stalenhoef AF
   21. van Crevel R
   22. Brunner HG
   23. Oh DY
   24. Schumann RR
   25. de la Rúa C
   26. Sauerwein R
   27. Kullberg BJ
   28. van der Ven AJ
   29. van der Meer JW
   30. Netea MG

   (2007) TLR4 polymorphisms, infectious diseases, and evolutionary pressure during migration of modern humans

   PNAS 104:16645–16650.

   https://doi.org/10.1073/pnas.0704828104

   • PubMed
   • Google Scholar
9. 1. Flandroy L
   2. Poutahidis T
   3. Berg G
   4. Clarke G
   5. Dao M-C
   6. Decaestecker E
   7. Furman E
   8. Haahtela T
   9. Massart S
   10. Plovier H
   11. Sanz Y
   12. Rook G

   (2018) The impact of human activities and lifestyles on the interlinked Microbiota and health of humans and of ecosystems

   Science of the Total Environment 627:1018–1038.

   https://doi.org/10.1016/j.scitotenv.2018.01.288

   • Google Scholar
10. 1. Fumagalli M
    2. Sironi M
    3. Pozzoli U
    4. Ferrer-Admetlla A
    5. Ferrer-Admettla A
    6. Pattini L
    7. Nielsen R

    (2011) Signatures of environmental genetic adaptation pinpoint pathogens as the main selective pressure through human evolution

    PLOS Genetics 7:e1002355.

    https://doi.org/10.1371/journal.pgen.1002355

    • PubMed
    • Google Scholar
11. 1. Fumagalli M
    2. Sironi M

    (2014) Human genome variability, natural selection and infectious diseases

    Current Opinion in Immunology 30:9–16.

    https://doi.org/10.1016/j.coi.2014.05.001

    • PubMed
    • Google Scholar
12. 1. Hofmanová Z
    2. Kreutzer S
    3. Hellenthal G
    4. Sell C
    5. Diekmann Y
    6. Díez-Del-Molino D
    7. van Dorp L
    8. López S
    9. Kousathanas A
    10. Link V
    11. Kirsanow K
    12. Cassidy LM
    13. Martiniano R
    14. Strobel M
    15. Scheu A
    16. Kotsakis K
    17. Halstead P
    18. Triantaphyllou S
    19. Kyparissi-Apostolika N
    20. Urem-Kotsou D
    21. Ziota C
    22. Adaktylou F
    23. Gopalan S
    24. Bobo DM
    25. Winkelbach L
    26. Blöcher J
    27. Unterländer M
    28. Leuenberger C
    29. Çilingiroğlu Ç
    30. Horejs B
    31. Gerritsen F
    32. Shennan SJ
    33. Bradley DG
    34. Currat M
    35. Veeramah KR
    36. Wegmann D
    37. Thomas MG
    38. Papageorgopoulou C
    39. Burger J

    (2016) Early farmers from across Europe directly descended from neolithic aegeans

    PNAS 113:6886–6891.

    https://doi.org/10.1073/pnas.1523951113

    • PubMed
    • Google Scholar
13. 1. Joosten SA
    2. van Meijgaarden KE
    3. Arend SM
    4. Prins C
    5. Oftung F
    6. Korsvold GE
    7. Kik SV
    8. Arts RJW
    9. van Crevel R
    10. Netea MG
    11. Ottenhoff THM

    (2018) Mycobacterial growth inhibition is associated with trained innate immunity

    Journal of Clinical Investigation 128:1837–1851.

    https://doi.org/10.1172/JCI97508

    • Google Scholar
14. 1. Karlsson EK
    2. Kwiatkowski DP
    3. Sabeti PC

    (2014) Natural selection and infectious disease in human populations

    Nature Reviews Genetics 15:379–393.

    https://doi.org/10.1038/nrg3734

    • PubMed
    • Google Scholar
15. 1. Li Y
    2. Oosting M
    3. Smeekens SP
    4. Jaeger M
    5. Aguirre-Gamboa R
    6. Le KTT
    7. Deelen P
    8. Ricaño-Ponce I
    9. Schoffelen T
    10. Jansen AFM
    11. Swertz MA
    12. Withoff S
    13. van de Vosse E
    14. van Deuren M
    15. van de Veerdonk F
    16. Zhernakova A
    17. van der Meer JWM
    18. Xavier RJ
    19. Franke L
    20. Joosten LAB
    21. Wijmenga C
    22. Kumar V
    23. Netea MG

    (2016) A functional genomics approach to understand variation in cytokine production in humans

    Cell 167:1099–1110.

    https://doi.org/10.1016/j.cell.2016.10.017

    • PubMed
    • Google Scholar
16. 1. Liu W
    2. Li Y
    3. Learn GH
    4. Rudicell RS
    5. Robertson JD
    6. Keele BF
    7. Ndjango JB
    8. Sanz CM
    9. Morgan DB
    10. Locatelli S
    11. Gonder MK
    12. Kranzusch PJ
    13. Walsh PD
    14. Delaporte E
    15. Mpoudi-Ngole E
    16. Georgiev AV
    17. Muller MN
    18. Shaw GM
    19. Peeters M
    20. Sharp PM
    21. Rayner JC
    22. Hahn BH

    (2010) Origin of the human malaria parasite plasmodium falciparum in gorillas

    Nature 467:420–425.

    https://doi.org/10.1038/nature09442

    • PubMed
    • Google Scholar
17. 1. MacArthur J
    2. Bowler E
    3. Cerezo M
    4. Gil L
    5. Hall P
    6. Hastings E
    7. Junkins H
    8. McMahon A
    9. Milano A
    10. Morales J
    11. Pendlington ZM
    12. Welter D
    13. Burdett T
    14. Hindorff L
    15. Flicek P
    16. Cunningham F
    17. Parkinson H

    (2017) The new NHGRI-EBI catalog of published genome-wide association studies (GWAS catalog)

    Nucleic Acids Research 45:D896–D901.

    https://doi.org/10.1093/nar/gkw1133

    • PubMed
    • Google Scholar
18. 1. Mei Y
    2. Pan F
    3. Gao J
    4. Ge R
    5. Duan Z
    6. Zeng Z
    7. Liao F
    8. Xia G
    9. Wang S
    10. Xu S
    11. Xu J
    12. Zhang L
    13. Ye D

    (2011) Increased serum IL-17 and IL-23 in the patient with ankylosing spondylitis

    Clinical Rheumatology 30:269–273.

    https://doi.org/10.1007/s10067-010-1647-4

    • PubMed
    • Google Scholar
19. 1. MRSA Systems Immunobiology Group
    2. Chan LC
    3. Rossetti M
    4. Miller LS
    5. Filler SG
    6. Johnson CW
    7. Lee HK
    8. Wang H
    9. Gjertson D
    10. Fowler VG
    11. Reed EF
    12. Yeaman MR

    (2018) Protective immunity in recurrent Staphylococcus aureus infection reflects localized immune signatures and macrophage-conferred memory

    PNAS 115:E11111–E11119.

    https://doi.org/10.1073/pnas.1808353115

    • PubMed
    • Google Scholar
20. 1. Netea MG
    2. Joosten LAB
    3. Li Y
    4. Kumar V
    5. Oosting M
    6. Smeekens S
    7. Jaeger M
    8. ter Horst R
    9. Schirmer M
    10. Vlamakis H
    11. Notebaart R
    12. Pavelka N
    13. Aguirre-Gamboa RR
    14. Swertz MA
    15. Tunjungputri RN
    16. van de Heijden W
    17. Franzosa EA
    18. Ng A
    19. Graham D
    20. Lassen K
    21. Schraa K
    22. Netea-Maier R
    23. Smit J
    24. de Mast Q
    25. van de Veerdonk F
    26. Kullberg BJ
    27. Tack C
    28. van de Munckhof I
    29. Rutten J
    30. van der Graaf J
    31. Franke L
    32. Hofker M
    33. Jonkers I
    34. Platteel M
    35. Maatman A
    36. Fu J
    37. Zhernakova A
    38. van der Meer JWM
    39. Dinarello CA
    40. van der Ven A
    41. Huttenhouwer C
    42. Koenen H
    43. Joosten I
    44. Xavier RJ
    45. Wijmenga C

    (2016) Understanding human immune function using the resources from the human functional genomics project

    Nature Medicine 22:831–833.

    https://doi.org/10.1038/nm.4140

    • Google Scholar
21. 1. Pai M
    2. Behr MA
    3. Dowdy D
    4. Dheda K
    5. Divangahi M
    6. Boehme CC
    7. Ginsberg A
    8. Swaminathan S
    9. Spigelman M
    10. Getahun H
    11. Menzies D
    12. Raviglione M

    (2016) Tuberculosis

    Nature Reviews Disease Primers 2:76.

    https://doi.org/10.1038/nrdp.2016.76

    • Google Scholar
22. 1. Poinar G

    (2005) Plasmodium dominicana n. sp. (Plasmodiidae: haemospororida) from tertiary dominican amber

    Systematic Parasitology 61:47–52.

    https://doi.org/10.1007/s11230-004-6354-6

    • PubMed
    • Google Scholar
23. 1. Purcell S
    2. Neale B
    3. Todd-Brown K
    4. Thomas L
    5. Ferreira MA
    6. Bender D
    7. Maller J
    8. Sklar P
    9. de Bakker PI
    10. Daly MJ
    11. Sham PC

    (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses

    The American Journal of Human Genetics 81:559–575.

    https://doi.org/10.1086/519795

    • PubMed
    • Google Scholar
24. 1. Quintana-Murci L
    2. Clark AG

    (2013) Population genetic tools for dissecting innate immunity in humans

    Nature Reviews Immunology 13:280–293.

    https://doi.org/10.1038/nri3421

    • PubMed
    • Google Scholar
25. 1. Racimo F
    2. Woodbridge J
    3. Fyfe RM
    4. Sikora M
    5. Sjögren KG
    6. Kristiansen K
    7. Vander Linden M

    (2020) The spatiotemporal spread of human migrations during the european holocene

    PNAS 117:8989–9000.

    https://doi.org/10.1073/pnas.1920051117

    • PubMed
    • Google Scholar
26. 1. Raj T
    2. Kuchroo M
    3. Replogle JM
    4. Raychaudhuri S
    5. Stranger BE
    6. De Jager PL

    (2013) Common risk alleles for inflammatory diseases are targets of recent positive selection

    The American Journal of Human Genetics 92:517–529.

    https://doi.org/10.1016/j.ajhg.2013.03.001

    • PubMed
    • Google Scholar
27. 1. Ramos PS
    2. Shedlock AM
    3. Langefeld CD

    (2015) Genetics of autoimmune diseases: insights from population genetics

    Journal of Human Genetics 60:657–664.

    https://doi.org/10.1038/jhg.2015.94

    • PubMed
    • Google Scholar
28. 1. Ruschen S
    2. Stellberg W
    3. Warnatz H

    (1992) Kinetics of cytokine secretion by mononuclear cells of the blood from rheumatoid arthritis patients are different from those of healthy controls

    Clinical & Experimental Immunology 89:32–37.

    https://doi.org/10.1111/j.1365-2249.1992.tb06873.x

    • PubMed
    • Google Scholar
29. 1. Schirmer M
    2. Smeekens SP
    3. Vlamakis H
    4. Jaeger M
    5. Oosting M
    6. Franzosa EA
    7. Ter Horst R
    8. Jansen T
    9. Jacobs L
    10. Bonder MJ
    11. Kurilshikov A
    12. Fu J
    13. Joosten LAB
    14. Zhernakova A
    15. Huttenhower C
    16. Wijmenga C
    17. Netea MG
    18. Xavier RJ

    (2016) Linking the human gut microbiome to inflammatory cytokine production capacity

    Cell 167:1125–1136.

    https://doi.org/10.1016/j.cell.2016.10.020

    • PubMed
    • Google Scholar
30. 1. Tay MZ
    2. Poh CM
    3. Rénia L
    4. MacAry PA
    5. Ng LFP

    (2020) The trinity of COVID-19: immunity, inflammation and intervention

    Nature Reviews Immunology 20:363–374.

    https://doi.org/10.1038/s41577-020-0311-8

    • PubMed
    • Google Scholar
31. 1. Ter Horst R
    2. Jaeger M
    3. Smeekens SP
    4. Oosting M
    5. Swertz MA
    6. Li Y
    7. Kumar V
    8. Diavatopoulos DA
    9. Jansen AFM
    10. Lemmers H
    11. Toenhake-Dijkstra H
    12. van Herwaarden AE
    13. Janssen M
    14. van der Molen RG
    15. Joosten I
    16. Sweep F
    17. Smit JW
    18. Netea-Maier RT
    19. Koenders M
    20. Xavier RJ
    21. van der Meer JWM
    22. Dinarello CA
    23. Pavelka N
    24. Wijmenga C
    25. Notebaart RA
    26. Joosten LAB
    27. Netea MG

    (2016) Host and environmental factors influencing individual human cytokine responses

    Cell 167:1111–1124.

    https://doi.org/10.1016/j.cell.2016.10.018

    • PubMed
    • Google Scholar
32. 1. Ter Horst R
    2. Jaeger M
    3. van de Wijer L
    4. van der Heijden WA
    5. Janssen AMW
    6. Smeekens SP
    7. Brouwer MAE
    8. van Cranenbroek B
    9. Aguirre-Gamboa R
    10. Netea-Maier RT
    11. van Herwaarden AE
    12. Lemmers H
    13. Dijkstra H
    14. Joosten I
    15. Koenen H
    16. Netea MG
    17. Joosten LAB

    (2021) Seasonal and nonseasonal longitudinal variation of immune function

    The Journal of Immunology 14:2000133.

    https://doi.org/10.4049/jimmunol.2000133

    • Google Scholar
33. 1. Thakur A
    2. Mikkelsen H
    3. Jungersen G

    (2019) Intracellular pathogens: host immunity and microbial persistence strategies

    Journal of Immunology Research 2019:1–24.

    https://doi.org/10.1155/2019/1356540

    • Google Scholar
34. 1. van de Veerdonk FL
    2. Marijnissen RJ
    3. Kullberg BJ
    4. Koenen HJ
    5. Cheng SC
    6. Joosten I
    7. van den Berg WB
    8. Williams DL
    9. van der Meer JW
    10. Joosten LA
    11. Netea MG

    (2009) The macrophage mannose receptor induces IL-17 in response to Candida Albicans

    Cell Host & Microbe 5:329–340.

    https://doi.org/10.1016/j.chom.2009.02.006

    • PubMed
    • Google Scholar
35. 1. Zhernakova A
    2. Elbers CC
    3. Ferwerda B
    4. Romanos J
    5. Trynka G
    6. Dubois PC
    7. de Kovel CG
    8. Franke L
    9. Oosting M
    10. Barisani D
    11. Bardella MT
    12. Joosten LA
    13. Saavalainen P
    14. van Heel DA
    15. Catassi C
    16. Netea MG
    17. Wijmenga C
    18. Finnish Celiac Disease Study Group

    (2010) Evolutionary and functional analysis of celiac risk loci reveals SH2B3 as a protective factor against bacterial infection

    The American Journal of Human Genetics 86:970–977.

    https://doi.org/10.1016/j.ajhg.2010.05.004

    • PubMed
    • Google Scholar

Author details

1. Jorge Domínguez-Andrés

   1. Department of Internal Medicine and Radboud Center for Infectious diseases (RCI), Radboud University Nijmegen Medical Centre, Nijmegen, Netherlands
   2. Radboud Institute for Molecular Life Sciences (RIMLS), RadboudUniversity Medical Center, Nijmegen, Netherlands

   Contribution

   Conceptualization, Investigation, Writing - original draft

   Contributed equally with

   Yunus Kuijpers

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-9091-1961

2. Yunus Kuijpers

   1. Department of Computational Biology for Individualised Infection Medicine, Centre for Individualised Infection Medicine (CiiM), a joint venture between Helmholtz-Centre for Infection Research (HZI) and the Hannover Medical School (MHH), Hannover, Germany
   2. TWINCORE, Centre for Experimental and Clinical Infection Research, a joint venture between Helmholtz-Centre for Infection Research (HZI) and the Hannover Medical School (MHH), Hannover, Germany

   Contribution

   Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing - original draft

   Contributed equally with

   Jorge Domínguez-Andrés

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-5075-3970

3. Olivier B Bakker

   Department of Genetics, University Medical Centre Groningen, Nijmegen, Netherlands

   Contribution

   Resources, Formal analysis, Investigation

   Competing interests

   No competing interests declared

4. Martin Jaeger

   1. Department of Internal Medicine and Radboud Center for Infectious diseases (RCI), Radboud University Nijmegen Medical Centre, Nijmegen, Netherlands
   2. Radboud Institute for Molecular Life Sciences (RIMLS), RadboudUniversity Medical Center, Nijmegen, Netherlands

   Contribution

   Resources, Writing - review and editing

   Competing interests

   No competing interests declared

5. Cheng-Jian Xu

   1. Department of Internal Medicine and Radboud Center for Infectious diseases (RCI), Radboud University Nijmegen Medical Centre, Nijmegen, Netherlands
   2. Department of Computational Biology for Individualised Infection Medicine, Centre for Individualised Infection Medicine (CiiM), a joint venture between Helmholtz-Centre for Infection Research (HZI) and the Hannover Medical School (MHH), Hannover, Germany
   3. TWINCORE, Centre for Experimental and Clinical Infection Research, a joint venture between Helmholtz-Centre for Infection Research (HZI) and the Hannover Medical School (MHH), Hannover, Germany

   Contribution

   Resources, Methodology, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-1586-4672

6. Jos WM Van der Meer

   Department of Internal Medicine and Radboud Center for Infectious diseases (RCI), Radboud University Nijmegen Medical Centre, Nijmegen, Netherlands

   Contribution

   Supervision, Writing - review and editing

   Competing interests

   Senior editor, eLife

7. Mattias Jakobsson

   1. Human Evolution, Department of Organismal Biology, Uppsala University, Uppsala, Sweden
   2. Centre for Anthropological Research, Department of Anthropology and Development Studies, University of Johannesburg, Auckland Park, South Africa

   Contribution

   Resources, Methodology

   Competing interests

   No competing interests declared

8. Jaume Bertranpetit

   Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, Barcelona, Spain

   Contribution

   Resources, Supervision, Writing - review and editing

   Competing interests

   No competing interests declared

9. Leo AB Joosten

   1. Department of Internal Medicine and Radboud Center for Infectious diseases (RCI), Radboud University Nijmegen Medical Centre, Nijmegen, Netherlands
   2. Radboud Institute for Molecular Life Sciences (RIMLS), RadboudUniversity Medical Center, Nijmegen, Netherlands

   Contribution

   Resources, Supervision

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-6166-9830

10. Yang Li

    1. Department of Internal Medicine and Radboud Center for Infectious diseases (RCI), Radboud University Nijmegen Medical Centre, Nijmegen, Netherlands
    2. Radboud Institute for Molecular Life Sciences (RIMLS), RadboudUniversity Medical Center, Nijmegen, Netherlands
    3. Department of Computational Biology for Individualised Infection Medicine, Centre for Individualised Infection Medicine (CiiM), a joint venture between Helmholtz-Centre for Infection Research (HZI) and the Hannover Medical School (MHH), Hannover, Germany
    4. TWINCORE, Centre for Experimental and Clinical Infection Research, a joint venture between Helmholtz-Centre for Infection Research (HZI) and the Hannover Medical School (MHH), Hannover, Germany

    Contribution

    Conceptualization, Resources, Supervision, Funding acquisition, Project administration, Writing - review and editing

    Contributed equally with

    Mihai G Netea

    For correspondence

    yang.li@helmholtz-hzi.de

    Competing interests

    No competing interests declared

    "This ORCID iD identifies the author of this article:" 0000-0003-4022-7341

11. Mihai G Netea

    1. Department of Internal Medicine and Radboud Center for Infectious diseases (RCI), Radboud University Nijmegen Medical Centre, Nijmegen, Netherlands
    2. Radboud Institute for Molecular Life Sciences (RIMLS), RadboudUniversity Medical Center, Nijmegen, Netherlands
    3. Department for Genomics & Immunoregulation, Life and Medical Sciences Institute (LIMES), University of Bonn, Bonn, Germany

    Contribution

    Conceptualization, Resources, Supervision, Funding acquisition, Project administration, Writing - review and editing

    Contributed equally with

    Yang Li

    For correspondence

    Mihai.Netea@radboudumc.nl

    Competing interests

    No competing interests declared

    "This ORCID iD identifies the author of this article:" 0000-0003-2421-6052

• 5,988

  views

• 565

  downloads

• 15

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.