Common DNA sequence variation influences epigenetic aging in African populations.

Aging is associated with genome-wide changes in DNA methylation in humans, facilitating the development of epigenetic age prediction models. However, most of these models have been trained primarily on European-ancestry individuals, and none account for the impact of methylation quantitative trait loci (meQTL). To address these gaps, we analyzed the relationships between age, genotype, and CpG methylation in 3 understudied populations: central African Baka (n = 35), southern African ‡Khomani San (n = 52), and southern African Himba (n = 51). We find that published prediction methods yield higher mean errors in these cohorts compared to European-ancestry individuals, and find that unaccounted-for DNA sequence variation may be a significant factor underlying this loss of accuracy. We leverage information about the associations between DNA genotype and CpG methylation to develop an age predictor that is minimally influenced by meQTL, and show that this model remains accurate across a broad range of genetic backgrounds. Intriguingly, we also find that the older individuals and those exhibiting relatively lower epigenetic age acceleration in our cohorts tend to carry more epigenetic age-reducing genetic variants, suggesting a novel mechanism by which heritable factors can influence longevity.

Deep estimation of the intensity and timing of natural selection from ancient genomes.

Leveraging past allele frequencies has proven to be key for identifying the impact of natural selection across time. However, this approach suffers from imprecise estimations of the intensity (s) and timing (T) of selection, particularly when ancient samples are scarce in specific epochs. Here, we aimed to bypass the computation of allele frequencies across arbitrarily defined past epochs and refine the estimations of selection parameters by implementing convolutional neural networks (CNNs) algorithms that directly use ancient genotypes sampled across time. Using computer simulations, we first show that genotype-based CNNs consistently outperform an approximate Bayesian computation (ABC) approach based on past allele frequency trajectories, regardless of the selection model assumed and the number of available ancient genotypes. When applying this method to empirical data from modern and ancient Europeans, we replicated the reported increased number of selection events in post-Neolithic Europe, independently of the continental subregion studied. Furthermore, we substantially refined the ABC-based estimations of s and T for a set of positively and negatively selected variants, including iconic cases of positive selection and experimentally validated disease-risk variants. Our CNN predictions support a history of recent positive and negative selection targeting variants associated with host defence against pathogens, aligning with previous work that highlights the significant impact of infectious diseases, such as tuberculosis, in Europe. These findings collectively demonstrate that detecting the footprints of natural selection on ancient genomes is crucial for unravelling the history of severe human diseases.

Conservation of antiviral systems across domains of life reveals immune genes in humans.

Deciphering the immune organization of eukaryotes is important for human health and for understanding ecosystems. The recent discovery of antiphage systems revealed that various eukaryotic immune proteins originate from prokaryotic antiphage systems. However, whether bacterial antiphage proteins can illuminate immune organization in eukaryotes remains unexplored. Here, we use a phylogeny-driven approach to uncover eukaryotic immune proteins by searching for homologs of bacterial antiphage systems. We demonstrate that proteins displaying sequence similarity with recently discovered antiphage systems are widespread in eukaryotes and maintain a role in human immunity. Two eukaryotic proteins of the anti-transposon piRNA pathway are evolutionarily linked to the antiphage system Mokosh. Additionally, human GTPases of immunity-associated proteins (GIMAPs) as well as two genes encoded in microsynteny, FHAD1 and CTRC, are respectively related to the Eleos and Lamassu prokaryotic systems and exhibit antiviral activity. Our work illustrates how comparative genomics of immune mechanisms can uncover defense genes in eukaryotes.

Impact of socioeconomic status on healthy immune responses in humans.

Individuals with low socioeconomic status (SES) are at greater risk of contracting and developing severe disease compared with people with higher SES. Age, sex, host genetics, smoking and cytomegalovirus (CMV) serostatus are known to have a major impact on human immune responses and thus susceptibility to infection. However, the impact of SES on immune variability is not well understood or explored. Here, we used data from the Milieu Intérieur project, a study of 1000 healthy volunteers with extensive demographic and biological data, to examine the effect of SES on immune variability. We developed an Elo-rating system using socioeconomic features such as education, income and home ownership status to objectively rank SES in the 1000 donors. We observed sex-specific SES associations, such as females with a low SES having a significantly higher frequency of CMV seropositivity compared with females with high SES, and males with a low SES having a significantly higher frequency of active smoking compared with males with a high SES. Using random forest models, we identified specific immune genes which were significantly associated with SES in both baseline and immune challenge conditions. Interestingly, many of the SES associations were sex stimuli specific, highlighting the complexity of these interactions. Our study provides a new way of computing SES in human populations that can help identify novel SES associations and reinforces biological evidence for SES-dependent susceptibility to infection. This should serve as a basis for further understanding the molecular mechanisms behind SES effects on immune responses and ultimately disease.

MANOCCA: a robust and computationally efficient test of covariance in high-dimension multivariate omics data.

Multivariate analysis is becoming central in studies investigating high-throughput molecular data, yet, some important features of these data are seldom explored. Here, we present MANOCCA (Multivariate Analysis of Conditional CovAriance), a powerful method to test for the effect of a predictor on the covariance matrix of a multivariate outcome. The proposed test is by construction orthogonal to tests based on the mean and variance and is able to capture effects that are missed by both approaches. We first compare the performances of MANOCCA with existing correlation-based methods and show that MANOCCA is the only test correctly calibrated in simulation mimicking omics data. We then investigate the impact of reducing the dimensionality of the data using principal component analysis when the sample size is smaller than the number of pairwise covariance terms analysed. We show that, in many realistic scenarios, the maximum power can be achieved with a limited number of components. Finally, we apply MANOCCA to 1000 healthy individuals from the Milieu Interieur cohort, to assess the effect of health, lifestyle and genetic factors on the covariance of two sets of phenotypes, blood biomarkers and flow cytometry-based immune phenotypes. Our analyses identify significant associations between multiple factors and the covariance of both omics data.

Smoking changes adaptive immunity with persistent effects.

Individuals differ widely in their immune responses, with age, sex and genetic factors having major roles in this inherent variability1-6. However, the variables that drive such differences in cytokine secretion-a crucial component of the host response to immune challenges-remain poorly defined. Here we investigated 136 variables and identified smoking, cytomegalovirus latent infection and body mass index as major contributors to variability in cytokine response, with effects of comparable magnitudes with age, sex and genetics. We find that smoking influences both innate and adaptive immune responses. Notably, its effect on innate responses is quickly lost after smoking cessation and is specifically associated with plasma levels of CEACAM6, whereas its effect on adaptive responses persists long after individuals quit smoking and is associated with epigenetic memory. This is supported by the association of the past smoking effect on cytokine responses with DNA methylation at specific signal trans-activators and regulators of metabolism. Our findings identify three novel variables associated with cytokine secretion variability and reveal roles for smoking in the short- and long-term regulation of immune responses. These results have potential clinical implications for the risk of developing infections, cancers or autoimmune diseases.

Standardized high-dimensional spectral cytometry protocol and panels for whole blood immune phenotyping in clinical and translational studies.

Flow cytometry is the method of choice for immunophenotyping in the context of clinical, translational, and systems immunology studies. Among the latter, the Milieu Intérieur (MI) project aims at defining the boundaries of a healthy immune response to identify determinants of immune response variation. MI used immunophenotyping of a 1000 healthy donor cohort by flow cytometry as a principal outcome for immune variance at steady state. New generation spectral cytometers now enable high-dimensional immune cell characterization from small sample volumes. Therefore, for the MI 10-year follow up study, we have developed two high-dimensional spectral flow cytometry panels for deep characterization of innate and adaptive whole blood immune cells (35 and 34 fluorescent markers, respectively). We have standardized the protocol for sample handling, staining, acquisition, and data analysis. This approach enables the reproducible quantification of over 182 immune cell phenotypes at a single site. We have applied the protocol to discern minor differences between healthy and patient samples and validated its value for application in immunomonitoring studies. Our protocol is currently used for characterization of the impact of age and environmental factors on peripheral blood immune phenotypes of >400 donors from the initial MI cohort.

Immune-response 3'UTR alternative polyadenylation quantitative trait loci contribute to variation in human complex traits and diseases.

Genome-wide association studies (GWASs) have identified thousands of non-coding variants that are associated with human complex traits and diseases. The analysis of such GWAS variants in different contexts and physiological states is essential for deciphering the regulatory mechanisms underlying human disease. Alternative polyadenylation (APA) is a key post-transcriptional modification for most human genes that substantially impacts upon cell behavior. Here, we mapped 9,493 3'-untranslated region APA quantitative trait loci in 18 human immune baseline cell types and 8 stimulation conditions (immune 3'aQTLs). Through the comparison between baseline and stimulation data, we observed the high responsiveness of 3'aQTLs to immune stimulation (response 3'aQTLs). Co-localization and mendelian randomization analyses of immune 3'aQTLs identified 678 genes where 3'aQTL are associated with variation in complex traits, 27.3% of which were derived from response 3'aQTLs. Overall, these analyses reveal the role of immune 3'aQTLs in the determination of complex traits, providing new insights into the regulatory mechanisms underlying disease etiologies.

Unveiling Candida albicans intestinal carriage in healthy volunteers: the role of micro- and mycobiota, diet, host genetics and immune response.

Candida albicans is a commensal yeast present in the gut of most healthy individuals but with highly variable concentrations. However, little is known about the host factors that influence colonization densities. We investigated how microbiota, host lifestyle factors, and genetics could shape C. albicans intestinal carriage in 695 healthy individuals from the Milieu Intérieur cohort. C. albicans intestinal carriage was detected in 82.9% of the subjects using quantitative PCR. Using linear mixed models and multiway-ANOVA, we explored C. albicans intestinal levels with regard to gut microbiota composition and lifestyle factors including diet. By analyzing shotgun metagenomics data and C. albicans qPCR data, we showed that Intestinimonas butyriciproducens was the only gut microbiota species whose relative abundance was negatively correlated with C. albicans concentration. Diet is also linked to C. albicans growth, with eating between meals and a low-sodium diet being associated with higher C. albicans levels. Furthermore, by Genome-Wide Association Study, we identified 26 single nucleotide polymorphisms suggestively associated with C. albicans colonization. In addition, we found that the intestinal levels of C. albicans might influence the host immune response, specifically in response to fungal challenge. We analyzed the transcriptional levels of 546 immune genes and the concentration of 13 cytokines after whole blood stimulation with C. albicans cells and showed positive associations between the extent of C. albicans intestinal levels and NLRP3 expression, as well as secreted IL-2 and CXCL5 concentrations. Taken together, these findings open the way for potential new interventional strategies to curb C. albicans intestinal overgrowth.

Dissecting human population variation in single-cell responses to SARS-CoV-2.

Humans display substantial interindividual clinical variability after SARS-CoV-2 infection1-3, the genetic and immunological basis of which has begun to be deciphered4. However, the extent and drivers of population differences in immune responses to SARS-CoV-2 remain unclear. Here we report single-cell RNA-sequencing data for peripheral blood mononuclear cells-from 222 healthy donors of diverse ancestries-that were stimulated with SARS-CoV-2 or influenza A virus. We show that SARS-CoV-2 induces weaker, but more heterogeneous, interferon-stimulated gene activity compared with influenza A virus, and a unique pro-inflammatory signature in myeloid cells. Transcriptional responses to viruses display marked population differences, primarily driven by changes in cell abundance including increased lymphoid differentiation associated with latent cytomegalovirus infection. Expression quantitative trait loci and mediation analyses reveal a broad effect of cell composition on population disparities in immune responses, with genetic variants exerting a strong effect on specific loci. Furthermore, we show that natural selection has increased population differences in immune responses, particularly for variants associated with SARS-CoV-2 response in East Asians, and document the cellular and molecular mechanisms by which Neanderthal introgression has altered immune functions, such as the response of myeloid cells to viruses. Finally, colocalization and transcriptome-wide association analyses reveal an overlap between the genetic basis of immune responses to SARS-CoV-2 and COVID-19 severity, providing insights into the factors contributing to current disparities in COVID-19 risk.

Integrated analysis of whole blood oxylipin and cytokine responses after bacterial, viral, and T cell stimulation reveals new immune networks.

Oxylipins are major immunomodulating mediators, yet studies of inflammation focus mainly on cytokines. Here, using a standardized whole-blood stimulation system, we characterized the oxylipin-driven inflammatory responses to various stimuli and their relationships with cytokine responses. We performed a pilot study in 25 healthy individuals using 6 different stimuli: 2 bacterial stimuli (LPS and live BCG), 2 viral stimuli (vaccine-grade poly I:C and live H1N1 attenuated influenza), an enterotoxin superantigen and a Null control. All stimuli induced a strong production of oxylipins but most importantly, bacterial, viral, and T cell immune responses show distinct oxylipin signatures. Integration of the oxylipin and cytokine responses for each condition revealed new immune networks improving our understanding of inflammation regulation. Finally, the oxylipin responses and oxylipin-cytokine networks were compared in patients with active tuberculosis or with latent infection. This revealed different responses to BCG but not LPS stimulation highlighting new regulatory pathways for further investigations.

High-throughput droplet-based analysis of influenza A virus genetic reassortment by single-virus RNA sequencing.

The segmented RNA genome of influenza A viruses (IAVs) enables viral evolution through genetic reassortment after multiple IAVs coinfect the same cell, leading to viruses harboring combinations of eight genomic segments from distinct parental viruses. Existing data indicate that reassortant genotypes are not equiprobable; however, the low throughput of available virology techniques does not allow quantitative analysis. Here, we have developed a high-throughput single-cell droplet microfluidic system allowing encapsulation of IAV-infected cells, each cell being infected by a single progeny virion resulting from a coinfection process. Customized barcoded primers for targeted viral RNA sequencing enabled the analysis of 18,422 viral genotypes resulting from coinfection with two circulating human H1N1pdm09 and H3N2 IAVs. Results were highly reproducible, confirmed that genetic reassortment is far from random, and allowed accurate quantification of reassortants including rare events. In total, 159 out of the 254 possible reassortant genotypes were observed but with widely varied prevalence (from 0.038 to 8.45%). In cells where eight segments were detected, all 112 possible pairwise combinations of segments were observed. The inclusion of data from single cells where less than eight segments were detected allowed analysis of pairwise cosegregation between segments with very high confidence. Direct coupling analysis accurately predicted the fraction of pairwise segments and full genotypes. Overall, our results indicate that a large proportion of reassortant genotypes can emerge upon coinfection and be detected over a wide range of frequencies, highlighting the power of our tool for systematic and exhaustive monitoring of the reassortment potential of IAVs.

Genetic adaptation to pathogens and increased risk of inflammatory disorders in post-Neolithic Europe.

Ancient genomics can directly detect human genetic adaptation to environmental cues. However, it remains unclear how pathogens have exerted selective pressures on human genome diversity across different epochs and affected present-day inflammatory disease risk. Here, we use an ancestry-aware approximate Bayesian computation framework to estimate the nature, strength, and time of onset of selection acting on 2,879 ancient and modern European genomes from the last 10,000 years. We found that the bulk of genetic adaptation occurred after the start of the Bronze Age, <4,500 years ago, and was enriched in genes relating to host-pathogen interactions. Furthermore, we detected directional selection acting on specific leukocytic lineages and experimentally demonstrated that the strongest negatively selected candidate variant in immunity genes, lipopolysaccharide-binding protein (LBP) D283G, is hypomorphic. Finally, our analyses suggest that the risk of inflammatory disorders has increased in post-Neolithic Europeans, possibly because of antagonistic pleiotropy following genetic adaptation to pathogens.

Unravelling the determinants of human health in French Polynesia: the MATAEA project.

Background: French Polynesia is a French overseas collectivity in the Southeast Pacific, comprising 75 inhabited islands across five archipelagoes. The human settlement of the region corresponds to the last massive migration of humans to empty territories, but its timeline is still debated. Despite their recent population history and geographical isolation, inhabitants of French Polynesia experience health issues similar to those of continental countries. Modern lifestyles and increased longevity have led to a rise in non-communicable diseases (NCDs) such as obesity, diabetes, hypertension, and cardiovascular diseases. Likewise, international trade and people mobility have caused the emergence of communicable diseases (CDs) including mosquito-borne and respiratory diseases. Additionally, chronic pathologies including acute rheumatic fever, liver diseases, and ciguatera, are highly prevalent in French Polynesia. However, data on such diseases are scarce and not representative of the geographic fragmentation of the population. Objectives: The present project aims to estimate the prevalence of several NCDs and CDs in the population of the five archipelagoes, and identify associated risk factors. Moreover, genetic analyses will contribute to determine the sequence and timings of the peopling history of French Polynesia, and identify causal links between past genetic adaptation to island environments, and present-day susceptibility to certain diseases. Methods: This cross-sectional survey is based on the random selection of 2,100 adults aged 18-69 years and residing on 18 islands from the five archipelagoes. Each participant answered a questionnaire on a wide range of topics (including demographic characteristics, lifestyle habits and medical history), underwent physical measurements (height, weight, waist circumference, arterial pressure, and skin pigmentation), and provided biological samples (blood, saliva, and stool) for biological, genetic and microbiological analyses. Conclusion: For the first time in French Polynesia, the present project allows to collect a wide range of data to explore the existence of indicators and/or risk factors for multiple pathologies of public health concern. The results will help health authorities to adapt actions and preventive measures aimed at reducing the incidence of NCDs and CDs. Moreover, the new genomic data generated in this study, combined with anthropological data, will increase our understanding of the peopling history of French Polynesia. Clinical trial registration: https://clinicaltrials.gov/, identifier: NCT06133400.

CIITA promoter polymorphism impairs monocytes HLA-DR expression in patients with septic shock.

Low monocyte (m)HLA-DR expression is associated with mortality in sepsis. G-286A∗rs3087456 polymorphism in promoter III of HLA class II transactivator (CIITA), the master regulator of HLA, has been associated with autoimmune diseases but its role in sepsis has never been demonstrated. In 203 patients in septic shock, GG genotype was associated with 28-day mortality and mHLA-DR remained low whereas it increased in patients with AA or AG genotype. In ex vivo cells, mHLA-DR failed to augment in GG in comparison with AG or AA genotype on exposure to IFN-γ. Promoter III transcript levels were similar in control monocytes regardless of genotype and exposure to IFN-γ. Promoter III activity was decreased in GG genotype in monocyte cell line but restored after stimulation with IFN-γ. Hereby, we demonstrated that G-286A∗rs3087456 significantly impact mHLA-DR expression in patients with septic shock in part through CIITA promoter III activity, that can be rescued using IFN-γ.

The immune factors driving DNA methylation variation in human blood.

Epigenetic changes are required for normal development, yet the nature and respective contribution of factors that drive epigenetic variation in humans remain to be fully characterized. Here, we assessed how the blood DNA methylome of 884 adults is affected by DNA sequence variation, age, sex and 139 factors relating to life habits and immunity. Furthermore, we investigated whether these effects are mediated or not by changes in cellular composition, measured by deep immunophenotyping. We show that DNA methylation differs substantially between naïve and memory T cells, supporting the need for adjustment on these cell-types. By doing so, we find that latent cytomegalovirus infection drives DNA methylation variation and provide further support that the increased dispersion of DNA methylation with aging is due to epigenetic drift. Finally, our results indicate that cellular composition and DNA sequence variation are the strongest predictors of DNA methylation, highlighting critical factors for medical epigenomics studies.

The genomic landscape of contemporary western Remote Oceanians.

The Vanuatu archipelago served as a gateway to Remote Oceania during one of the most extensive human migrations to uninhabited lands ∼3,000 years ago. Ancient DNA studies suggest an initial settlement by East Asian-related peoples that was quickly followed by the arrival of Papuan-related populations, leading to a major population turnover. Yet there is uncertainty over the population processes and the sociocultural factors that have shaped the genomic diversity of ni-Vanuatu, who present nowadays among the world's highest linguistic and cultural diversity. Here, we report new genome-wide data for 1,433 contemporary ni-Vanuatu from 29 different islands, including 287 couples. We find that ni-Vanuatu derive their East Asian- and Papuan-related ancestry from the same source populations and descend from relatively synchronous, sex-biased admixture events that occurred ∼1,700-2,300 years ago, indicating a peopling history common to the whole archipelago. However, East Asian-related ancestry proportions differ markedly across islands, suggesting that the Papuan-related population turnover was geographically uneven. Furthermore, we detect Polynesian ancestry arriving ∼600-1,000 years ago to Central and South Vanuatu in both Polynesian-speaking and non-Polynesian-speaking populations. Last, we provide evidence for a tendency of spouses to carry similar genetic ancestry, when accounting for relatedness avoidance. The signal is not driven by strong genetic effects of specific loci or trait-associated variants, suggesting that it results instead from social assortative mating. Altogether, our findings provide an insight into both the genetic history of ni-Vanuatu populations and how sociocultural processes have shaped the diversity of their genomes.