Evolution of cytokine production capacity in ancient and modern European populations
Abstract
As our ancestors migrated throughout different continents, natural selection increased the presence of alleles advantageous in the new environments. Heritable variations that alter the susceptibility to diseases vary with the historical period, the virulence of the infections, and their geographical spread. In this study we built polygenic scores for heritable traits that influence the genetic adaptation in the production of cytokines and immune-mediated disorders, including infectious, inflammatory, and autoimmune diseases, and applied them to the genomes of several ancient European populations. We observed that the advent of the Neolithic was a turning point for immune-mediated traits in Europeans, favoring those alleles linked with the development of tolerance against intracellular pathogens and promoting inflammatory responses against extracellular microbes. These evolutionary patterns are also associated with an increased presence of traits related to inflammatory and auto-immune diseases.
Introduction
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).
Results and discussion
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.
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.
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).
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).
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.
Materials and methods
Cohort selection
Request a detailed protocolAncient DNA genotype data were downloaded from version 37.2 of the published aDNA genotype database, compiled by and available on the David Reich Lab website (https://reich.hms.harvard.edu/downloadable-genotypes-present-day-and-ancient-dna-data-compiled-published-papers). The aDNA samples consisted of pseudo-haploid genotype data. This was due to the low genotyping coverage. Samples with variant missingness above 96% were filtered out using Plink (Purcell et al., 2007). This was done in order to remove outliers with extremely low coverage. Only samples within Europe were used for this study, and these samples were selected based on their geographic location, that is latitude (within 35 and 70 degrees north) and longitude (within 10 degrees west and 40 degrees east). Samples without a carbon-dated age were also filtered out. We also selected 250 European samples from the 1000 genomes project phase 3. Only variants present in both the ancient samples and the modern samples were retained. This resulted in a dataset of 827 ancient samples and 250 modern samples containing 1,233,013 variants.
Carbon-dated sample origin and geographical location
Request a detailed protocolBoth carbon-dated age of origin and latitudinal and longitudinal data were available for these 827 ancient European samples. Broad time periods were assigned to these samples with the early upper Paleolithic era for all samples originating from before 25,000 years before the common era standardized to 1950 (BCE). The late upper Paleolithic era follows until 11,000 BCE. The Mesolithic era ranges from 11,000 to 5500 BCE. The Neolithic era ranges from 8500 to 3900 BCE, and the post-Neolithic era ranges from 5000 BCE and more recent ages. Using the geographical data in combination with archeological clues and the genetic data, the broad time period of origin was also available for samples that were dated to a point in time with overlapping broad time periods. This allowed the samples to be classified as either early upper Paleolithic, late upper Paleolithic, Mesolithic, Neolithic, or post-Neolithic. The sample age of the 250 modern European samples was set to 0.
Summary statistics of GWAS and cytokine QTLs
Request a detailed protocolSummary statistics for complex traits were obtained from the UK Biobank (Bycroft et al., 2018) and the GWAS catalog (MacArthur et al., 2017) last accessed on March 29, 2020. The stimulated cytokine response summary statistics from the 500FG cohort of the HFGP were used (Li et al., 2016). Some complex traits had multiple different sets of summary statistics available. In these cases, the data, which were more recent and used bigger cohorts that were either of European or mixed (European and Asian) ancestry, were selected. The variants of these summary statistics were then filtered by only keeping bi-allelic variants. Most aDNA genotypes available are pseudo-haploid as a consequence of their lower sample quality. We excluded ambiguous SNPs (A/T and C/G) in order to prevent errors due to strand flips present in these pseudo-haploid samples.
Polygenic risk score calculation
Request a detailed protocolPolygenic risk scores were then calculated by first intersecting the filtered variants from the summary statistics with the variants present in the DNA samples. Starting at the most significant variant, all variants within a 250-kb window around that variant were excluded until no variants remained. We then multiplied the dosage of these variants with the effect size and these values were summed. If a variant is missing in a sample, the dosage is substituted with the average genotyped dosage for that variant within the entire dataset. This way the PRS is not skewed in any specific direction. The formula for this is described below with the score S being the weighted sum of a variant’s dosage Xn multiplied by its associated weight or βn calculated using m variants.
Relation between PRS and carbon-dated sample age
Request a detailed protocolWe constructed piecewise linear models for each trait by separating the samples into two groups. These two groups consisted of all samples preceding the Neolithic era and those of the Neolithic era and later, respectively. We correlated PRS with the carbon-dated age of our samples. We then multiplied the -log10 of the Pearson correlation p-values with the sign of the correlation coefficients.
In addition, we plotted LOESS regression models to highlight the change in PRS at each point in time independent of any predefined breakpoint between historical periods. We also performed a group-based comparison using Student’s t-test. We compared pre-Neolithic, Neolithic, post-Neolithic, and Mesolithic samples with their respective adjacent historical periods to show the difference between other historical transitions besides the pre- and post-Neolithic.
Selection test
Request a detailed protocolWe tested whether traits were under selection or if observed changes were due to the genetic drift between adjacent time periods. We performed a two-tailed test using the mean Fst calculated with trait-specific SNPs between two adjacent periods and a reference distribution of 10,000 random linkage disequilibrium (LD)- and minor allele frequency (MAF)-matched mean Fst scores calculated using an equal amount of SNPs (Schirmer et al., 2016). Bonferroni correction was performed to account for multiple testing.
Robustness of results
Request a detailed protocolIn order to test the robustness of our results we calculated PRS using multiple different p-value thresholds for QTL inclusion. We used p-value thresholds from 10−3 to 10−8 for the complex traits obtained through GWAS catalog and the UK Biobank. The thresholds used for the stimulated cytokine responses ranged from 10−3 to 10−6. We also calculated PRS using different variant missingness thresholds. This means we removed samples with a variant missingness rate higher than 96, 90, 80, or 70%. All of the results from the piecewise linear models were then used to create a heatmap depicting the consistency and robustness of our observed correlations.
Additionally, various window sizes were used for clumping the QTLs, and LD-based clumping was also performed excluding variants with an LD greater than 0.2 compared to our lead SNP within a window. In order to see whether our observations were due to sample imbalances between the pre-Neolithic period and the later periods, samples originating from the Neolithic period and later were randomly down-sampled to the same number of samples as the pre-Neolithic samples. Correlation coefficients between PRS and sample age were then recalculated for the Neolithic and younger samples and compared to the coefficients obtained using all Neolithic and younger samples.
Data availability
All the data employed in this manuscript have been obtained from publicly available databases.
-
NHGRI-EBI GWAS CatalogID 22057235. Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease.
-
NHGRI-EBI GWAS CatalogID 26192919. Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations.
-
NHGRI-EBI GWAS CatalogID 23749187. Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci".
-
NHGRI-EBI GWAS CatalogID 23603761. Dense genotyping of immune-related disease regions identifies 14 new susceptibility loci for juvenile idiopathic arthritis.
-
NHGRI-EBI GWAS CatalogID 24390342. Genetics of rheumatoid arthritis contributes to biology and drug discovery.
-
NHGRI-EBI GWAS CatalogID 24076602. Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis.
-
David Reich Lab published data amh_repo V37.2ID V37.2. Downloadable genotypes of present-day and ancient DNA data (compiled from published papers) V37.2.
-
The 1000ibd MOLGENIS databaseID 500FG/cytokines. 500FG Cytokine QTL data.
References
-
Interleukin-6 in juvenile idiopathic arthritisModern Rheumatology 29:275–286.https://doi.org/10.1080/14397595.2019.1574697
-
Genomic signatures of selective pressures and introgression from archaic hominins at human innate immunity genesThe American Journal of Human Genetics 98:5–21.https://doi.org/10.1016/j.ajhg.2015.11.014
-
Evolution of bacterial virulenceFEMS Microbiology Reviews 41:679–697.https://doi.org/10.1093/femsre/fux023
-
Reappraisal of the historical selective pressures for the CCR5-Delta32 mutationJournal of Medical Genetics 42:205–208.https://doi.org/10.1136/jmg.2004.025346
-
The impact of human activities and lifestyles on the interlinked Microbiota and health of humans and of ecosystemsScience of the Total Environment 627:1018–1038.https://doi.org/10.1016/j.scitotenv.2018.01.288
-
Human genome variability, natural selection and infectious diseasesCurrent Opinion in Immunology 30:9–16.https://doi.org/10.1016/j.coi.2014.05.001
-
Mycobacterial growth inhibition is associated with trained innate immunityJournal of Clinical Investigation 128:1837–1851.https://doi.org/10.1172/JCI97508
-
Natural selection and infectious disease in human populationsNature Reviews Genetics 15:379–393.https://doi.org/10.1038/nrg3734
-
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
-
Increased serum IL-17 and IL-23 in the patient with ankylosing spondylitisClinical Rheumatology 30:269–273.https://doi.org/10.1007/s10067-010-1647-4
-
PLINK: a tool set for whole-genome association and population-based linkage analysesThe American Journal of Human Genetics 81:559–575.https://doi.org/10.1086/519795
-
Population genetic tools for dissecting innate immunity in humansNature Reviews Immunology 13:280–293.https://doi.org/10.1038/nri3421
-
Common risk alleles for inflammatory diseases are targets of recent positive selectionThe American Journal of Human Genetics 92:517–529.https://doi.org/10.1016/j.ajhg.2013.03.001
-
Genetics of autoimmune diseases: insights from population geneticsJournal of Human Genetics 60:657–664.https://doi.org/10.1038/jhg.2015.94
-
The trinity of COVID-19: immunity, inflammation and interventionNature Reviews Immunology 20:363–374.https://doi.org/10.1038/s41577-020-0311-8
-
Seasonal and nonseasonal longitudinal variation of immune functionThe Journal of Immunology 14:2000133.https://doi.org/10.4049/jimmunol.2000133
-
Intracellular pathogens: host immunity and microbial persistence strategiesJournal of Immunology Research 2019:1–24.https://doi.org/10.1155/2019/1356540
-
Evolutionary and functional analysis of celiac risk loci reveals SH2B3 as a protective factor against bacterial infectionThe American Journal of Human Genetics 86:970–977.https://doi.org/10.1016/j.ajhg.2010.05.004
Article and author information
Author details
Funding
European Commission (833247)
- Mihai G Netea
European Commission (948207)
- Yang Li
ZonMw (09150161910024)
- Jorge Domínguez-Andrés
Netherlands Organisation for Scientific Research (Spinoza)
- Mihai G Netea
Radboud University Medical Center (Hypatia Grant (2018))
- Yang Li
Netherlands Organisation for Scientific Research (09150161910024)
- Jorge Domínguez-Andrés
MINECO (PID2019-110933GB-I00 (AEI/FEDER UE))
- Jaume Bertranpetit
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
MGN was supported by an ERC Advanced Grant (833247) and a Spinoza Grant of the Netherlands Organization for Scientific Research. YL was supported by an ERC Starting Grant (948207) and the Radboud University Medical Centre Hypatia Grant (2018) for Scientific Research. JD-A is supported by The Netherlands Organization for Scientific Research (VENI grant 09150161910024). JB was supported by PID2019-110933GB-I00 (AEI/FEDER, UE) MINECO, Spain.
Copyright
© 2021, Domínguez-Andrés et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 6,196
- views
-
- 597
- downloads
-
- 18
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Genetics and Genomics
Repurposing of pleiotropic factors during execution of diverse cellular processes has emerged as a regulatory paradigm. Embryonic development in metazoans is controlled by maternal factors deposited in the egg during oogenesis. Here, we explore maternal role(s) of Caspar (Casp), the Drosophila orthologue of human Fas-associated factor-1 (FAF1) originally implicated in host-defense as a negative regulator of NF-κB signaling. Maternal loss of either Casp or it’s protein partner, transitional endoplasmic reticulum 94 (TER94) leads to partial embryonic lethality correlated with aberrant centrosome behavior, cytoskeletal abnormalities, and defective gastrulation. Although ubiquitously distributed, both proteins are enriched in the primordial germ cells (PGCs), and in keeping with the centrosome problems, mutant embryos display a significant reduction in the PGC count. Moreover, the total number of pole buds is directly proportional to the level of Casp. Consistently, it’s ‘loss’ and ‘gain’ results in respective reduction and increase in the Oskar protein levels, the master determinant of PGC fate. To elucidate this regulatory loop, we analyzed several known components of mid-blastula transition and identify the translational repressor Smaug, a zygotic regulator of germ cell specification, as a potential critical target. We present a detailed structure-function analysis of Casp aimed at understanding its novel involvement during PGC development.
-
- Genetics and Genomics
Recent studies have revealed a role for zinc in insulin secretion and glucose homeostasis. Randomized placebo-controlled zinc supplementation trials have demonstrated improved glycemic traits in patients with type II diabetes (T2D). Moreover, rare loss-of-function variants in the zinc efflux transporter SLC30A8 reduce T2D risk. Despite this accumulated evidence, a mechanistic understanding of how zinc influences systemic glucose homeostasis and consequently T2D risk remains unclear. To further explore the relationship between zinc and metabolic traits, we searched the exome database of the Regeneron Genetics Center-Geisinger Health System DiscovEHR cohort for genes that regulate zinc levels and associate with changes in metabolic traits. We then explored our main finding using in vitro and in vivo models. We identified rare loss-of-function (LOF) variants (MAF <1%) in Solute Carrier Family 39, Member 5 (SLC39A5) associated with increased circulating zinc (p=4.9 × 10-4). Trans-ancestry meta-analysis across four studies exhibited a nominal association of SLC39A5 LOF variants with decreased T2D risk. To explore the mechanisms underlying these associations, we generated mice lacking Slc39a5. Slc39a5-/- mice display improved liver function and reduced hyperglycemia when challenged with congenital or diet-induced obesity. These improvements result from elevated hepatic zinc levels and concomitant activation of hepatic AMPK and AKT signaling, in part due to zinc-mediated inhibition of hepatic protein phosphatase activity. Furthermore, under conditions of diet-induced non-alcoholic steatohepatitis (NASH), Slc39a5-/- mice display significantly attenuated fibrosis and inflammation. Taken together, these results suggest SLC39A5 as a potential therapeutic target for non-alcoholic fatty liver disease (NAFLD) due to metabolic derangements including T2D.