Clinical, immunological, and genetic description of a Mexican cohort of patients with DOCK8 deficiency

PURPOSE

We aimed to describe the clinical, immunological, and genetic features of patients with DOCK8 deficiency (DOCK8-Def) in a tertiary care center for children.

METHODS

Retrospective chart review of patients' clinical, immunological, and genetic characteristics with DOCK8-Def. Genetic analysis was performed with targeted- or whole-exome sequencing; we also assessed DOCK8 protein expression and a lymphoproliferation assay and analyzed survival by the Kaplan-Meier method.

RESULTS

We described 11 patients from 8 unrelated kindreds. The median age at symptoms' onset was 10 months (range 1-54 months). The median follow-up time was 53.4 months (4.8-118.8). All patients presented eczema and recurrent sinopulmonary and cutaneous infections. Besides those symptoms, the most frequent manifestations were bronchiectases (8/11), food allergies (6/11), and severe infections (6/11). Infrequent characteristics were detection of CMV in bronchial lavage, C. parvum-driven sclerosing cholangitis, Takayasu vasculitis, neurological syndromes, pulmonary tuberculosis, and lymphomatoid granulomatosis.

CONCLUSION

DOCK8-Def has a broad spectrum of manifestations, including allergy, autoimmunity, inflammation, infection, and cancer. The hallmark of this inborn error of immunity is IEI-associated eczema with eosinophilia and increased IgE. Here, we report six new mutations causing human DOCK8 deficiency and symptoms previously unrecognized to occur in DOCK8-Def. Therefore, an early diagnosis of DOCK8-Def is essential to facilitate an adequate treatment such as HSCT.

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.

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.

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.

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.

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.

Early IFNβ secretion determines variable downstream IL-12p70 responses upon TLR4 activation

The interleukin-12 (IL-12) family comprises the only heterodimeric cytokines mediating diverse functional effects. We previously reported a striking bimodal IL-12p70 response to lipopolysaccharide (LPS) stimulation in healthy donors. Herein, we demonstrate that interferon β (IFNβ) is a major upstream determinant of IL-12p70 production, which is also associated with numbers and activation of circulating monocytes. Integrative modeling of proteomic, genetic, epigenomic, and cellular data confirms IFNβ as key for LPS-induced IL-12p70 and allowed us to compare the relative effects of each of these parameters on variable cytokine responses. Clinical relevance of our findings is supported by reduced IFNβ-IL-12p70 responses in patients hospitalized with acute severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection or chronically infected with hepatitis C (HCV). Importantly, these responses are resolved after viral clearance. Our systems immunology approach defines a better understanding of IL-12p70 and IFNβ in healthy and infected persons, providing insights into how common genetic and epigenetic variation may impact immune responses to bacterial infection.

The genomic signatures of natural selection in admixed human populations

Admixture has been a pervasive phenomenon in human history, extensively shaping the patterns of population genetic diversity. There is increasing evidence to suggest that admixture can also facilitate genetic adaptation to local environments, i.e., admixed populations acquire beneficial mutations from source populations, a process that we refer to as "adaptive admixture." However, the role of adaptive admixture in human evolution and the power to detect it remain poorly characterized. Here, we use extensive computer simulations to evaluate the power of several neutrality statistics to detect natural selection in the admixed population, assuming multiple admixture scenarios. We show that statistics based on admixture proportions, F_adm and LAD, show high power to detect mutations that are beneficial in the admixed population, whereas other statistics, including iHS and F_ST, falsely detect neutral mutations that have been selected in the source populations only. By combining F_adm and LAD into a single, powerful statistic, we scanned the genomes of 15 worldwide, admixed populations for signatures of adaptive admixture. We confirm that lactase persistence and resistance to malaria have been under adaptive admixture in West Africans and in Malagasy, North Africans, and South Asians, respectively. Our approach also uncovers other cases of adaptive admixture, including APOL1 in Fulani nomads and PKN2 in East Indonesians, involved in resistance to infection and metabolism, respectively. Collectively, our study provides evidence that adaptive admixture has occurred in human populations whose genetic history is characterized by periods of isolation and spatial expansions resulting in increased gene flow.

Integrative genetic and immune cell analysis of plasma proteins in healthy donors identifies novel associations involving primary immune deficiency genes

Background

Blood plasma proteins play an important role in immune defense against pathogens, including cytokine signaling, the complement system, and the acute-phase response. Recent large-scale studies have reported genetic (i.e., protein quantitative trait loci, pQTLs) and non-genetic factors, such as age and sex, as major determinants to inter-individual variability in immune response variation. However, the contribution of blood-cell composition to plasma protein heterogeneity has not been fully characterized and may act as a mediating factor in association studies.

Methods

Here, we evaluated plasma protein levels from 400 unrelated healthy individuals of western European ancestry, who were stratified by sex and two decades of life (20–29 and 60–69 years), from the Milieu Intérieur cohort. We quantified 229 proteins by Luminex in a clinically certified laboratory and their levels of variation were analyzed together with 5.2 million single-nucleotide polymorphisms. With respect to non-genetic variables, we included 254 lifestyle and biochemical factors, as well as counts of seven circulating immune cell populations measured by hemogram and standardized flow cytometry.

Results

Collectively, we found 152 significant associations involving 49 proteins and 20 non-genetic variables. Consistent with previous studies, age and sex showed a global, pervasive impact on plasma protein heterogeneity, while body mass index and other health status variables were among the non-genetic factors with the highest number of associations. After controlling for these covariates, we identified 100 and 12 pQTLs acting in cis and trans, respectively, collectively associated with 87 plasma proteins and including 19 novel genetic associations. Genetic factors explained the largest fraction of the variability of plasma protein levels, as compared to non-genetic factors. In addition, blood-cell fractions, including leukocytes, lymphocytes, monocytes, neutrophils, eosinophils, basophils, and platelets, had a larger contribution to inter-individual variability than age and sex and appeared as confounders of specific genetic associations. Finally, we identified new genetic associations with plasma protein levels of five monogenic Mendelian disease genes including two primary immunodeficiency genes (Ficolin-3 and FAS).

Conclusions

Our study identified novel genetic and non-genetic factors associated to plasma protein levels which may inform health status and disease management.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01032-y.

Sporadic occurrence of recent selective sweeps from standing variation in humans as revealed by an approximate Bayesian computation approach

During their dispersals over the last 100,000 years, modern humans have been exposed to a large variety of environments, resulting in genetic adaptation. While genome-wide scans for the footprints of positive Darwinian selection have increased knowledge of genes and functions potentially involved in human local adaptation, they have globally produced evidence of a limited contribution of selective sweeps in humans. Conversely, studies based on machine learning algorithms suggest that recent sweeps from standing variation are widespread in humans, an observation that has been recently questioned. Here, we sought to formally quantify the number of recent selective sweeps in humans, by leveraging approximate Bayesian computation and whole-genome sequence data. Our computer simulations revealed suitable ABC estimations, regardless of the frequency of the selected alleles at the onset of selection and the completion of sweeps. Under a model of recent selection from standing variation, we inferred that an average of 68 (from 56 to 79) and 140 (from 94 to 198) sweeps occurred over the last 100,000 years of human history, in African and Eurasian populations, respectively. The former estimation is compatible with human adaptation rates estimated since divergence with chimps, and reveals numbers of sweeps per generation per site in the range of values estimated in Drosophila. Our results confirm the rarity of selective sweeps in humans and show a low contribution of sweeps from standing variation to recent human adaptation.

Gut microbiome stability and dynamics in healthy donors and patients with non-gastrointestinal cancers

As microbial therapeutics are increasingly being tested in diverse patient populations, it is essential to understand the host and environmental factors influencing the microbiome. Through analysis of 1,359 gut microbiome samples from 946 healthy donors of the Milieu Intérieur cohort, we detail how microbiome composition is associated with host factors, lifestyle parameters, and disease states. Using a genome-based taxonomy, we found biological sex was the strongest driver of community composition. Additionally, bacterial populations shift across decades of life (age 20-69), with Bacteroidota species consistently increased with age while Actinobacteriota species, including Bifidobacterium, decreased. Longitudinal sampling revealed that short-term stability exceeds interindividual differences. By accounting for these factors, we defined global shifts in the microbiomes of patients with non-gastrointestinal tumors compared with healthy donors. Together, these results demonstrated that the microbiome displays predictable variations as a function of sex, age, and disease state. These variations must be considered when designing microbiome-targeted therapies or interpreting differences thought to be linked to pathophysiology or therapeutic response.

Genomic insights into population history and biological adaptation in Oceania

The Pacific region is of major importance for addressing questions regarding human dispersals, interactions with archaic hominins and natural selection processes¹. However, the demographic and adaptive history of Oceanian populations remains largely uncharacterized. Here we report high-coverage genomes of 317 individuals from 20 populations from the Pacific region. We find that the ancestors of Papuan-related ('Near Oceanian') groups underwent a strong bottleneck before the settlement of the region, and separated around 20,000-40,000 years ago. We infer that the East Asian ancestors of Pacific populations may have diverged from Taiwanese Indigenous peoples before the Neolithic expansion, which is thought to have started from Taiwan around 5,000 years ago^2-4. Additionally, this dispersal was not followed by an immediate, single admixture event with Near Oceanian populations, but involved recurrent episodes of genetic interactions. Our analyses reveal marked differences in the proportion and nature of Denisovan heritage among Pacific groups, suggesting that independent interbreeding with highly structured archaic populations occurred. Furthermore, whereas introgression of Neanderthal genetic information facilitated the adaptation of modern humans related to multiple phenotypes (for example, metabolism, pigmentation and neuronal development), Denisovan introgression was primarily beneficial for immune-related functions. Finally, we report evidence of selective sweeps and polygenic adaptation associated with pathogen exposure and lipid metabolism in the Pacific region, increasing our understanding of the mechanisms of biological adaptation to island environments.

Human ancient DNA analyses reveal the high burden of tuberculosis in Europeans over the last 2,000 years

Tuberculosis (TB), usually caused by Mycobacterium tuberculosis bacteria, is the first cause of death from an infectious disease at the worldwide scale, yet the mode and tempo of TB pressure on humans remain unknown. The recent discovery that homozygotes for the P1104A polymorphism of TYK2 are at higher risk to develop clinical forms of TB provided the first evidence of a common, monogenic predisposition to TB, offering a unique opportunity to inform on human co-evolution with a deadly pathogen. Here, we investigate the history of human exposure to TB by determining the evolutionary trajectory of the TYK2 P1104A variant in Europe, where TB is considered to be the deadliest documented infectious disease. Leveraging a large dataset of 1,013 ancient human genomes and using an approximate Bayesian computation approach, we find that the P1104A variant originated in the common ancestors of West Eurasians ∼30,000 years ago. Furthermore, we show that, following large-scale population movements of Anatolian Neolithic farmers and Eurasian steppe herders into Europe, P1104A has markedly fluctuated in frequency over the last 10,000 years of European history, with a dramatic decrease in frequency after the Bronze Age. Our analyses indicate that such a frequency drop is attributable to strong negative selection starting ∼2,000 years ago, with a relative fitness reduction on homozygotes of 20%, among the highest in the human genome. Together, our results provide genetic evidence that TB has imposed a heavy burden on European health over the last two millennia.

Negative selection on human genes underlying inborn errors depends on disease outcome and both the mode and mechanism of inheritance

Genetic variants underlying life-threatening diseases, being unlikely to be transmitted to the next generation, are gradually and selectively eliminated from the population through negative selection. We study the determinants of this evolutionary process in human genes underlying monogenic diseases by comparing various negative selection scores and an integrative approach, CoNeS, at 366 loci underlying inborn errors of immunity (IEI). We find that genes underlying autosomal dominant (AD) or X-linked IEI have stronger negative selection scores than those underlying autosomal recessive (AR) IEI, whose scores are not different from those of genes not known to be disease causing. Nevertheless, genes underlying AR IEI that are lethal before reproductive maturity with complete penetrance have stronger negative selection scores than other genes underlying AR IEI. We also show that genes underlying AD IEI by loss of function have stronger negative selection scores than genes underlying AD IEI by gain of function, while genes underlying AD IEI by haploinsufficiency are under stronger negative selection than other genes underlying AD IEI. These results are replicated in 1,140 genes underlying inborn errors of neurodevelopment. Finally, we propose a supervised classifier, SCoNeS, which predicts better than state-of-the-art approaches whether a gene is more likely to underlie an AD or AR disease. The clinical outcomes of monogenic inborn errors, together with their mode and mechanisms of inheritance, determine the levels of negative selection at their corresponding loci. Integrating scores of negative selection may facilitate the prioritization of candidate genes and variants in patients suspected to carry an inborn error.

Common homozygosity for predicted loss-of-function variants reveals both redundant and advantageous effects of dispensable human genes

Humans homozygous or hemizygous for variants predicted to cause a loss of function (LoF) of the corresponding protein do not necessarily present with overt clinical phenotypes. We report here 190 autosomal genes with 207 predicted LoF variants, for which the frequency of homozygous individuals exceeds 1% in at least one human population from five major ancestry groups. No such genes were identified on the X and Y chromosomes. Manual curation revealed that 28 variants (15%) had been misannotated as LoF. Of the 179 remaining variants in 166 genes, only 11 alleles in 11 genes had previously been confirmed experimentally to be LoF. The set of 166 dispensable genes was enriched in olfactory receptor genes (41 genes). The 41 dispensable olfactory receptor genes displayed a relaxation of selective constraints similar to that observed for other olfactory receptor genes. The 125 dispensable nonolfactory receptor genes also displayed a relaxation of selective constraints consistent with greater redundancy. Sixty-two of these 125 genes were found to be dispensable in at least three human populations, suggesting possible evolution toward pseudogenes. Of the 179 LoF variants, 68 could be tested for two neutrality statistics, and 8 displayed robust signals of positive selection. These latter variants included a known FUT2 variant that confers resistance to intestinal viruses, and an APOL3 variant involved in resistance to parasitic infections. Overall, the identification of 166 genes for which a sizeable proportion of humans are homozygous for predicted LoF alleles reveals both redundancies and advantages of such deficiencies for human survival.

Quantitative genetic analysis deciphers the impact of cis and trans regulation on cell-to-cell variability in protein expression levels

Identifying the factors that shape protein expression variability in complex multi-cellular organisms has primarily focused on promoter architecture and regulation of single-cell expression in cis. However, this targeted approach has to date been unable to identify major regulators of cell-to-cell gene expression variability in humans. To address this, we have combined single-cell protein expression measurements in the human immune system using flow cytometry with a quantitative genetics analysis. For the majority of proteins whose variability in expression has a heritable component, we find that genetic variants act in trans, with notably fewer variants acting in cis. Furthermore, we highlight using Mendelian Randomization that these variability-Quantitative Trait Loci might be driven by the cis regulation of upstream genes. This indicates that natural selection may balance the impact of gene regulation in cis with downstream impacts on expression variability in trans.

Two common disease-associated TYK2 variants impact exon splicing and TYK2 dosage

TYK2 belongs to the JAK protein tyrosine kinase family and mediates signaling of numerous antiviral and immunoregulatory cytokines (type I and type III IFNs, IL-10, IL-12, IL-22, IL-23) in immune and non-immune cells. After many years of genetic association studies, TYK2 is recognized as a susceptibility gene for some inflammatory and autoimmune diseases (AID). Seven TYK2 variants have been associated with AIDs in Europeans, and establishing their causality remains challenging. Previous work showed that a protective variant (P1104A) is hypomorphic and also a risk allele for mycobacterial infection. Here, we have studied two AID-associated common TYK2 variants: rs12720270 located in intron 7 and rs2304256, a non-synonymous variant in exon 8 that causes a valine to phenylalanine substitution (c.1084 G > T, Val362Phe). We found that this amino acid substitution does not alter TYK2 expression, catalytic activity or ability to relay signaling in EBV-B cell lines or in reconstituted TYK2-null cells. Based on in silico predictions that these variants may impact splicing of exon 8, we: i) analyzed TYK2 transcripts in genotyped EBV-B cells and in CRISPR/Cas9-edited cells, ii) measured splicing using minigene assays, and iii) performed eQTL (expression quantitative trait locus) analysis of TYK2 transcripts in primary monocytes and whole blood cells. Our results reveal that the two variants promote the inclusion of exon 8, which, we demonstrate, is essential for TYK2 binding to cognate receptors. In addition and in line with GTEx (Genetic Tissue Expression) data, our eQTL results show that rs2304256 mildly enhances TYK2 expression in whole blood. In all, these findings suggest that these TYK2 variants are not neutral but instead have a potential impact in AID.

The Genomic Impact of European Colonization of the Americas

The human genetic diversity of the Americas has been affected by several events of gene flow that have continued since the colonial era and the Atlantic slave trade. Moreover, multiple waves of migration followed by local admixture occurred in the last two centuries, the impact of which has been largely unexplored. Here, we compiled a genome-wide dataset of ∼12,000 individuals from twelve American countries and ∼6,000 individuals from worldwide populations and applied haplotype-based methods to investigate how historical movements from outside the New World affected (1) the genetic structure, (2) the admixture profile, (3) the demographic history, and (4) sex-biased gene-flow dynamics of the Americas. We revealed a high degree of complexity underlying the genetic contribution of European and African populations in North and South America, from both geographic and temporal perspectives, identifying previously unreported sources related to Italy, the Middle East, and to specific regions of Africa.

Associations of autozygosity with a broad range of human phenotypes

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

A comprehensive assessment of demographic, environmental, and host genetic associations with gut microbiome diversity in healthy individuals

BACKGROUND

The gut microbiome is an important determinant of human health. Its composition has been shown to be influenced by multiple environmental factors and likely by host genetic variation. In the framework of the Milieu Intérieur Consortium, a total of 1000 healthy individuals of western European ancestry, with a 1:1 sex ratio and evenly stratified across five decades of life (age 20-69), were recruited. We generated 16S ribosomal RNA profiles from stool samples for 858 participants. We investigated genetic and non-genetic factors that contribute to individual differences in fecal microbiome composition.

RESULTS

Among 110 demographic, clinical, and environmental factors, 11 were identified as significantly correlated with α-diversity, ß-diversity, or abundance of specific microbial communities in multivariable models. Age and blood alanine aminotransferase levels showed the strongest associations with microbiome diversity. In total, all non-genetic factors explained 16.4% of the variance. We then searched for associations between > 5 million single nucleotide polymorphisms and the same indicators of fecal microbiome diversity, including the significant non-genetic factors as covariates. No genome-wide significant associations were identified after correction for multiple testing. A small fraction of previously reported associations between human genetic variants and specific taxa could be replicated in our cohort, while no replication was observed for any of the diversity metrics.

CONCLUSION

In a well-characterized cohort of healthy individuals, we identified several non-genetic variables associated with fecal microbiome diversity. In contrast, host genetics only had a negligible influence. Demographic and environmental factors are thus the main contributors to fecal microbiome composition in healthy individuals.

TRIAL REGISTRATION

ClinicalTrials.gov identifier NCT01699893.

Human IFN-γ immunity to mycobacteria is governed by both IL-12 and IL-23

Hundreds of patients with autosomal recessive, complete IL-12p40 or IL-12Rβ1 deficiency have been diagnosed over the last 20 years. They typically suffer from invasive mycobacteriosis and, occasionally, from mucocutaneous candidiasis. Susceptibility to these infections is thought to be due to impairments of IL-12-dependent IFN-γ immunity and IL-23-dependent IL-17A/IL-17F immunity, respectively. We report here patients with autosomal recessive, complete IL-12Rβ2 or IL-23R deficiency, lacking responses to IL-12 or IL-23 only, all of whom, unexpectedly, display mycobacteriosis without candidiasis. We show that αβ T, γδ T, B, NK, ILC1, and ILC2 cells from healthy donors preferentially produce IFN-γ in response to IL-12, whereas NKT cells and MAIT cells preferentially produce IFN-γ in response to IL-23. We also show that the development of IFN-γ-producing CD4⁺ T cells, including, in particular, mycobacterium-specific T_H1* cells (CD45RA^-CCR6⁺), is dependent on both IL-12 and IL-23. Last, we show that IL12RB1, IL12RB2, and IL23R have similar frequencies of deleterious variants in the general population. The comparative rarity of symptomatic patients with IL-12Rβ2 or IL-23R deficiency, relative to IL-12Rβ1 deficiency, is, therefore, due to lower clinical penetrance. There are fewer symptomatic IL-23R- and IL-12Rβ2-deficient than IL-12Rβ1-deficient patients, not because these genetic disorders are rarer, but because the isolated absence of IL-12 or IL-23 is, in part, compensated by the other cytokine for the production of IFN-γ, thereby providing some protection against mycobacteria. These experiments of nature show that human IL-12 and IL-23 are both required for optimal IFN-γ-dependent immunity to mycobacteria, both individually and much more so cooperatively.

Tuberculosis and impaired IL-23-dependent IFN-γ immunity in humans homozygous for a common TYK2 missense variant

Inherited IL-12Rβ1 and TYK2 deficiencies impair both IL-12- and IL-23-dependent IFN-γ immunity and are rare monogenic causes of tuberculosis, each found in less than 1/600,000 individuals. We show that homozygosity for the common TYK2 P1104A allele, which is found in about 1/600 Europeans and between 1/1000 and 1/10,000 individuals in regions other than East Asia, is more frequent in a cohort of patients with tuberculosis from endemic areas than in ethnicity-adjusted controls (P = 8.37 × 10^-8; odds ratio, 89.31; 95% CI, 14.7 to 1725). Moreover, the frequency of P1104A in Europeans has decreased, from about 9% to 4.2%, over the past 4000 years, consistent with purging of this variant by endemic tuberculosis. Surprisingly, we also show that TYK2 P1104A impairs cellular responses to IL-23, but not to IFN-α, IL-10, or even IL-12, which, like IL-23, induces IFN-γ via activation of TYK2 and JAK2. Moreover, TYK2 P1104A is properly docked on cytokine receptors and can be phosphorylated by the proximal JAK, but lacks catalytic activity. Last, we show that the catalytic activity of TYK2 is essential for IL-23, but not IL-12, responses in cells expressing wild-type JAK2. In contrast, the catalytic activity of JAK2 is redundant for both IL-12 and IL-23 responses, because the catalytically inactive P1057A JAK2, which is also docked and phosphorylated, rescues signaling in cells expressing wild-type TYK2. In conclusion, homozygosity for the catalytically inactive P1104A missense variant of TYK2 selectively disrupts the induction of IFN-γ by IL-23 and is a common monogenic etiology of tuberculosis.

Homozygosity for TYK2 P1104A underlies tuberculosis in about 1% of patients in a cohort of European ancestry

The human genetic basis of tuberculosis (TB) has long remained elusive. We recently reported a high level of enrichment in homozygosity for the common TYK2 P1104A variant in a heterogeneous cohort of patients with TB from non-European countries in which TB is endemic. This variant is homozygous in ∼1/600 Europeans and ∼1/5,000 people from other countries outside East Asia and sub-Saharan Africa. We report a study of this variant in the UK Biobank cohort. The frequency of P1104A homozygotes was much higher in patients with TB (6/620, 1%) than in controls (228/114,473, 0.2%), with an odds ratio (OR) adjusted for ancestry of 5.0 [95% confidence interval (CI): 1.96-10.31, P = 2 × 10-3]. Conversely, we did not observe enrichment for P1104A heterozygosity, or for TYK2 I684S or V362F homozygosity or heterozygosity. Moreover, it is unlikely that more than 10% of controls were infected with Mycobacterium tuberculosis, as 97% were of European genetic ancestry, born between 1939 and 1970, and resided in the United Kingdom. Had all of them been infected, the OR for developing TB upon infection would be higher. These findings suggest that homozygosity for TYK2 P1104A may account for ∼1% of TB cases in Europeans.

Recent Adaptive Acquisition by African Rainforest Hunter-Gatherers of the Late Pleistocene Sickle-Cell Mutation Suggests Past Differences in Malaria Exposure

The hemoglobin β^S sickle mutation is a textbook case in which natural selection maintains a deleterious mutation at high frequency in the human population. Homozygous individuals for this mutation develop sickle-cell disease, whereas heterozygotes benefit from higher protection against severe malaria. Because the overdominant β^S allele should be purged almost immediately from the population in the absence of malaria, the study of the evolutionary history of this iconic mutation can provide important information about the history of human exposure to malaria. Here, we sought to increase our understanding of the origins and time depth of the β^S mutation in populations with different lifestyles and ecologies, and we analyzed the diversity of HBB in 479 individuals from 13 populations of African farmers and rainforest hunter-gatherers. Using an approximate Bayesian computation method, we estimated the age of the β^S allele while explicitly accounting for population subdivision, past demography, and balancing selection. When the effects of balancing selection are taken into account, our analyses indicate a single emergence of β^S in the ancestors of present-day agriculturalist populations ∼22,000 years ago. Furthermore, we show that rainforest hunter-gatherers have more recently acquired the β^S mutation from the ancestors of agriculturalists through adaptive gene flow during the last ∼6,000 years. Together, our results provide evidence for a more ancient exposure to malarial pressures among the ancestors of agriculturalists than previously appreciated, and they suggest that rainforest hunter-gatherers have been increasingly exposed to malaria during the last millennia.

Genome-wide data from the Bubi of Bioko Island clarifies the Atlantic fringe of the Bantu dispersal

Background

Bioko is one of the few islands that exist around Africa, the most genetically diverse continent on the planet. The native Bantu-speaking inhabitants of Bioko, the Bubi, are believed to have colonized the island about 2000 years ago. Here, we sequenced the genome of thirteen Bubi individuals at high coverage and analysed their sequences in comparison to mainland populations from the Gulf of Guinea.

Results

We found that, genetically, the closest mainland population to the Bubi are Bantu-speaking groups from Angola instead the geographically closer groups from Cameroon. The Bubi possess a lower proportion of rainforest hunter-gatherer (RHG) ancestry than most other Bantu-speaking groups. However, their RHG component most likely came from the same source and could have reached them by gene flow from the mainland after island settlement. By studying identity by descent (IBD) genomic blocks and runs of homozygosity (ROHs), we found evidence for a significant level of genetic isolation among the Bubi, isolation that can be attributed to the island effect. Additionally, as this population is known to have one of the highest malaria incidence rates in the world we analysed their genome for malaria-resistant alleles. However, we were unable to detect any specific selective sweeps related to this disease.

Conclusions

By describing their dispersal to the Atlantic islands, the genomic characterization of the Bubi contributes to the understanding of the margins of the massive Bantu migration that shaped all Sub-Saharan African populations.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5529-0) contains supplementary material, which is available to authorized users.

Exploring the genetic basis of human population differences in DNA methylation and their causal impact on immune gene regulation

BACKGROUND

DNA methylation is influenced by both environmental and genetic factors and is increasingly thought to affect variation in complex traits and diseases. Yet, the extent of ancestry-related differences in DNA methylation, their genetic determinants, and their respective causal impact on immune gene regulation remain elusive.

RESULTS

We report extensive population differences in DNA methylation between 156 individuals of African and European descent, detected in primary monocytes that are used as a model of a major innate immunity cell type. Most of these differences (~ 70%) are driven by DNA sequence variants nearby CpG sites, which account for ~ 60% of the variance in DNA methylation. We also identify several master regulators of DNA methylation variation in trans, including a regulatory hub nearby the transcription factor-encoding CTCF gene, which contributes markedly to ancestry-related differences in DNA methylation. Furthermore, we establish that variation in DNA methylation is associated with varying gene expression levels following mostly, but not exclusively, a canonical model of negative associations, particularly in enhancer regions. Specifically, we find that DNA methylation highly correlates with transcriptional activity of 811 and 230 genes, at the basal state and upon immune stimulation, respectively. Finally, using a Bayesian approach, we estimate causal mediation effects of DNA methylation on gene expression in ~ 20% of the studied cases, indicating that DNA methylation can play an active role in immune gene regulation.

CONCLUSION

Using a system-level approach, our study reveals substantial ancestry-related differences in DNA methylation and provides evidence for their causal impact on immune gene regulation.

Human thymopoiesis is influenced by a common genetic variant within the TCRA-TCRD locus

The thymus is the primary lymphoid organ where naïve T cells are generated; however, with the exception of age, the parameters that govern its function in healthy humans remain unknown. We characterized the variability of thymic function among 1000 age- and sex-stratified healthy adults of the Milieu Intérieur cohort, using quantification of T cell receptor excision circles (TRECs) in peripheral blood T cells as a surrogate marker of thymopoiesis. Age and sex were the only nonheritable factors identified that affect thymic function. TREC amounts decreased with age and were higher in women compared to men. In addition, a genome-wide association study revealed a common variant (rs2204985) within the T cell receptor TCRA-TCRD locus, between the DD2 and DD3 gene segments, which associated with TREC amounts. Strikingly, transplantation of human hematopoietic stem cells with the rs2204985 GG genotype into immunodeficient mice led to thymopoiesis with higher TRECs, increased thymocyte counts, and a higher TCR repertoire diversity. Our population immunology approach revealed a genetic locus that influences thymopoiesis in healthy adults, with potentially broad implications in precision medicine.

The human CIB1-EVER1-EVER2 complex governs keratinocyte-intrinsic immunity to β-papillomaviruses

Patients with epidermodysplasia verruciformis (EV) and biallelic null mutations of TMC6 (encoding EVER1) or TMC8 (EVER2) are selectively prone to disseminated skin lesions due to keratinocyte-tropic human β-papillomaviruses (β-HPVs), which lack E5 and E8. We describe EV patients homozygous for null mutations of the CIB1 gene encoding calcium- and integrin-binding protein-1 (CIB1). CIB1 is strongly expressed in the skin and cultured keratinocytes of controls but not in those of patients. CIB1 forms a complex with EVER1 and EVER2, and CIB1 proteins are not expressed in EVER1- or EVER2-deficient cells. The known functions of EVER1 and EVER2 in human keratinocytes are not dependent on CIB1, and CIB1 deficiency does not impair keratinocyte adhesion or migration. In keratinocytes, the CIB1 protein interacts with the HPV E5 and E8 proteins encoded by α-HPV16 and γ-HPV4, respectively, suggesting that this protein acts as a restriction factor against HPVs. Collectively, these findings suggest that the disruption of CIB1-EVER1-EVER2-dependent keratinocyte-intrinsic immunity underlies the selective susceptibility to β-HPVs of EV patients.

Human genetic variants and age are the strongest predictors of humoral immune responses to common pathogens and vaccines

BACKGROUND

Humoral immune responses to infectious agents or vaccination vary substantially among individuals, and many of the factors responsible for this variability remain to be defined. Current evidence suggests that human genetic variation influences (i) serum immunoglobulin levels, (ii) seroconversion rates, and (iii) intensity of antigen-specific immune responses. Here, we evaluated the impact of intrinsic (age and sex), environmental, and genetic factors on the variability of humoral response to common pathogens and vaccines.

METHODS

We characterized the serological response to 15 antigens from common human pathogens or vaccines, in an age- and sex-stratified cohort of 1000 healthy individuals (Milieu Intérieur cohort). Using clinical-grade serological assays, we measured total IgA, IgE, IgG, and IgM levels, as well as qualitative (serostatus) and quantitative IgG responses to cytomegalovirus, Epstein-Barr virus, herpes simplex virus 1 and 2, varicella zoster virus, Helicobacter pylori, Toxoplasma gondii, influenza A virus, measles, mumps, rubella, and hepatitis B virus. Following genome-wide genotyping of single nucleotide polymorphisms and imputation, we examined associations between ~ 5 million genetic variants and antibody responses using single marker and gene burden tests.

RESULTS

We identified age and sex as important determinants of humoral immunity, with older individuals and women having higher rates of seropositivity for most antigens. Genome-wide association studies revealed significant associations between variants in the human leukocyte antigen (HLA) class II region on chromosome 6 and anti-EBV and anti-rubella IgG levels. We used HLA imputation to fine map these associations to amino acid variants in the peptide-binding groove of HLA-DRβ1 and HLA-DPβ1, respectively. We also observed significant associations for total IgA levels with two loci on chromosome 2 and with specific KIR-HLA combinations.

CONCLUSIONS

Using extensive serological testing and genome-wide association analyses in a well-characterized cohort of healthy individuals, we demonstrated that age, sex, and specific human genetic variants contribute to inter-individual variability in humoral immunity. By highlighting genes and pathways implicated in the normal antibody response to frequently encountered antigens, these findings provide a basis to better understand disease pathogenesis.

TRIALS REGISTRATION

ClinicalTrials.gov , NCT01699893.

Mutations involving the SRY-related gene SOX8 are associated with a spectrum of human reproductive anomalies

SOX8 is an HMG-box transcription factor closely related to SRY and SOX9. Deletion of the gene encoding Sox8 in mice causes reproductive dysfunction but the role of SOX8 in humans is unknown. Here, we show that SOX8 is expressed in the somatic cells of the early developing gonad in the human and influences human sex determination. We identified two individuals with 46, XY disorders/differences in sex development (DSD) and chromosomal rearrangements encompassing the SOX8 locus and a third individual with 46, XY DSD and a missense mutation in the HMG-box of SOX8. In vitro functional assays indicate that this mutation alters the biological activity of the protein. As an emerging body of evidence suggests that DSDs and infertility can have common etiologies, we also analysed SOX8 in a cohort of infertile men (n = 274) and two independent cohorts of women with primary ovarian insufficiency (POI; n = 153 and n = 104). SOX8 mutations were found at increased frequency in oligozoospermic men (3.5%; P < 0.05) and POI (5.06%; P = 4.5 × 10-5) as compared with fertile/normospermic control populations (0.74%). The mutant proteins identified altered SOX8 biological activity as compared with the wild-type protein. These data demonstrate that SOX8 plays an important role in human reproduction and SOX8 mutations contribute to a spectrum of phenotypes including 46, XY DSD, male infertility and 46, XX POI.

IRF4 haploinsufficiency in a family with Whipple's disease

Most humans are exposed to Tropheryma whipplei (Tw). Whipple's disease (WD) strikes only a small minority of individuals infected with Tw (<0.01%), whereas asymptomatic chronic carriage is more common (<25%). We studied a multiplex kindred, containing four WD patients and five healthy Tw chronic carriers. We hypothesized that WD displays autosomal dominant (AD) inheritance, with age-dependent incomplete penetrance. We identified a single very rare non-synonymous mutation in the four patients: the private R98W variant of IRF4, a transcription factor involved in immunity. The five Tw carriers were younger, and also heterozygous for R98W. We found that R98W was loss-of-function, modified the transcriptome of heterozygous leukocytes following Tw stimulation, and was not dominant-negative. We also found that only six of the other 153 known non-synonymous IRF4 variants were loss-of-function. Finally, we found that IRF4 had evolved under purifying selection. AD IRF4 deficiency can underlie WD by haploinsufficiency, with age-dependent incomplete penetrance.

Dispersals and genetic adaptation of Bantu-speaking populations in Africa and North America

Bantu languages are spoken by about 310 million Africans, yet the genetic history of Bantu-speaking populations remains largely unexplored. We generated genomic data for 1318 individuals from 35 populations in western central Africa, where Bantu languages originated. We found that early Bantu speakers first moved southward, through the equatorial rainforest, before spreading toward eastern and southern Africa. We also found that genetic adaptation of Bantu speakers was facilitated by admixture with local populations, particularly for the HLA and LCT loci. Finally, we identified a major contribution of western central African Bantu speakers to the ancestry of African Americans, whose genomes present no strong signals of natural selection. Together, these results highlight the contribution of Bantu-speaking peoples to the complex genetic history of Africans and African Americans.

The Genetic Legacy of the Indian Ocean Slave Trade: Recent Admixture and Post-admixture Selection in the Makranis of Pakistan

From the eighth century onward, the Indian Ocean was the scene of extensive trade of sub-Saharan African slaves via sea routes controlled by Muslim Arab and Swahili traders. Several populations in present-day Pakistan and India are thought to be the descendants of such slaves, yet their history of admixture and natural selection remains largely undefined. Here, we studied the genome-wide diversity of the African-descent Makranis, who reside on the Arabian Sea coast of Pakistan, as well that of four neighboring Pakistani populations, to investigate the genetic legacy, population dynamics, and tempo of the Indian Ocean slave trade. We show that the Makranis are the result of an admixture event between local Baluch tribes and Bantu-speaking populations from eastern or southeastern Africa; we dated this event to ∼300 years ago during the Omani Empire domination. Levels of parental relatedness, measured through runs of homozygosity, were found to be similar across Pakistani populations, suggesting that the Makranis rapidly adopted the traditional practice of endogamous marriages. Finally, we searched for signatures of post-admixture selection at traits evolving under positive selection, including skin color, lactase persistence, and resistance to malaria. We demonstrate that the African-specific Duffy-null blood group-believed to confer resistance against Plasmodium vivax infection-was recently introduced to Pakistan through the slave trade and evolved adaptively in this P. vivax malaria-endemic region. Our study reconstructs the genetic and adaptive history of a neglected episode of the African Diaspora and illustrates the impact of recent admixture on the diffusion of adaptive traits across human populations.

Dual T cell- and B cell-intrinsic deficiency in humans with biallelic RLTPR mutations

In two complementary papers, Casanova, Malissen, and collaborators report the discovery of human RLTPR deficiency, the first primary immunodeficiency of the human CD28 pathway in T cells. Together, the two studies highlight the important and largely (but not completely) overlapping roles of RLTPR in T and B cells of humans and mice.

Combined immunodeficiency (CID) refers to inborn errors of human T cells that also affect B cells because of the T cell deficit or an additional B cell–intrinsic deficit. In this study, we report six patients from three unrelated families with biallelic loss-of-function mutations in RLTPR, the mouse orthologue of which is essential for CD28 signaling. The patients have cutaneous and pulmonary allergy, as well as a variety of bacterial and fungal infectious diseases, including invasive tuberculosis and mucocutaneous candidiasis. Proportions of circulating regulatory T cells and memory CD4⁺ T cells are reduced. Their CD4⁺ T cells do not respond to CD28 stimulation. Their CD4⁺ T cells exhibit a "Th2" cell bias ex vivo and when cultured in vitro, contrasting with the paucity of "Th1," "Th17," and T follicular helper cells. The patients also display few memory B cells and poor antibody responses. This B cell phenotype does not result solely from the T cell deficiency, as the patients’ B cells fail to activate NF-κB upon B cell receptor (BCR) stimulation. Human RLTPR deficiency is a CID affecting at least the CD28-responsive pathway in T cells and the BCR-responsive pathway in B cells.

Genomic Signatures of Selective Pressures and Introgression from Archaic Hominins at Human Innate Immunity Genes

Human genes governing innate immunity provide a valuable tool for the study of the selective pressure imposed by microorganisms on host genomes. A comprehensive, genome-wide study of how selective constraints and adaptations have driven the evolution of innate immunity genes is missing. Using full-genome sequence variation from the 1000 Genomes Project, we first show that innate immunity genes have globally evolved under stronger purifying selection than the remainder of protein-coding genes. We identify a gene set under the strongest selective constraints, mutations in which are likely to predispose individuals to life-threatening disease, as illustrated by STAT1 and TRAF3. We then evaluate the occurrence of local adaptation and detect 57 high-scoring signals of positive selection at innate immunity genes, variation in which has been associated with susceptibility to common infectious or autoimmune diseases. Furthermore, we show that most adaptations targeting coding variation have occurred in the last 6,000-13,000 years, the period at which populations shifted from hunting and gathering to farming. Finally, we show that innate immunity genes present higher Neandertal introgression than the remainder of the coding genome. Notably, among the genes presenting the highest Neandertal ancestry, we find the TLR6-TLR1-TLR10 cluster, which also contains functional adaptive variation in Europeans. This study identifies highly constrained genes that fulfill essential, non-redundant functions in host survival and reveals others that are more permissive to change-containing variation acquired from archaic hominins or adaptive variants in specific populations-improving our understanding of the relative biological importance of innate immunity pathways in natural conditions.

The epigenomic landscape of African rainforest hunter-gatherers and farmers

The genetic history of African populations is increasingly well documented, yet their patterns of epigenomic variation remain uncharacterized. Moreover, the relative impacts of DNA sequence variation and temporal changes in lifestyle and habitat on the human epigenome remain unknown. Here we generate genome-wide genotype and DNA methylation profiles for 362 rainforest hunter-gatherers and sedentary farmers. We find that the current habitat and historical lifestyle of a population have similarly critical impacts on the methylome, but the biological functions affected strongly differ. Specifically, methylation variation associated with recent changes in habitat mostly concerns immune and cellular functions, whereas that associated with historical lifestyle affects developmental processes. Furthermore, methylation variation—particularly that correlated with historical lifestyle—shows strong associations with nearby genetic variants that, moreover, are enriched in signals of natural selection. Our work provides new insight into the genetic and environmental factors affecting the epigenomic landscape of human populations over time.

Genetic and environmental factors affect genome-wide patterns of epigenetic variation. Here, the authors show that while current habitat and historical lifestyle impact the methylome of rainforest hunter-gatherers and sedentary farmers, the biological functions affected and the degree of genetic control differ.

Extensive Admixture and Selective Pressure Across the Sahel Belt

Genome-wide studies of African populations have the potential to reveal powerful insights into the evolution of our species, as these diverse populations have been exposed to intense selective pressures imposed by infectious diseases, diet, and environmental factors. Within Africa, the Sahel Belt extensively overlaps the geographical center of several endemic infections such as malaria, trypanosomiasis, meningitis, and hemorrhagic fevers. We screened 2.5 million single nucleotide polymorphisms in 161 individuals from 13 Sahelian populations, which together with published data cover Western, Central, and Eastern Sahel, and include both nomadic and sedentary groups. We confirmed the role of this Belt as a main corridor for human migrations across the continent. Strong admixture was observed in both Central and Eastern Sahelian populations, with North Africans and Near Eastern/Arabians, respectively, but it was inexistent in Western Sahelian populations. Genome-wide local ancestry inference in admixed Sahelian populations revealed several candidate regions that were significantly enriched for non-autochthonous haplotypes, and many showed to be under positive selection. The DARC gene region in Arabs and Nubians was enriched for African ancestry, whereas the RAB3GAP1/LCT/MCM6 region in Oromo, the TAS2R gene family in Fulani, and the ALMS1/NAT8 in Turkana and Samburu were enriched for non-African ancestry. Signals of positive selection varied in terms of geographic amplitude. Some genomic regions were selected across the Belt, the most striking example being the malaria-related DARC gene. Others were Western-specific (oxytocin, calcium, and heart pathways), Eastern-specific (lipid pathways), or even population-restricted (TAS2R genes in Fulani, which may reflect sexual selection).

Impact of common risk factors of fibrosis progression in chronic hepatitis C

OBJECTIVE

The natural course of chronic hepatitis C varies widely. To improve the profiling of patients at risk of developing advanced liver disease, we assessed the relative contribution of factors for liver fibrosis progression in hepatitis C.

DESIGN

We analysed 1461 patients with chronic hepatitis C with an estimated date of infection and at least one liver biopsy. Risk factors for accelerated fibrosis progression rate (FPR), defined as ≥ 0.13 Metavir fibrosis units per year, were identified by logistic regression. Examined factors included age at infection, sex, route of infection, HCV genotype, body mass index (BMI), significant alcohol drinking (≥ 20 g/day for ≥ 5 years), HIV coinfection and diabetes. In a subgroup of 575 patients, we assessed the impact of single nucleotide polymorphisms previously associated with fibrosis progression in genome-wide association studies. Results were expressed as attributable fraction (AF) of risk for accelerated FPR.

RESULTS

Age at infection (AF 28.7%), sex (AF 8.2%), route of infection (AF 16.5%) and HCV genotype (AF 7.9%) contributed to accelerated FPR in the Swiss Hepatitis C Cohort Study, whereas significant alcohol drinking, anti-HIV, diabetes and BMI did not. In genotyped patients, variants at rs9380516 (TULP1), rs738409 (PNPLA3), rs4374383 (MERTK) (AF 19.2%) and rs910049 (major histocompatibility complex region) significantly added to the risk of accelerated FPR. Results were replicated in three additional independent cohorts, and a meta-analysis confirmed the role of age at infection, sex, route of infection, HCV genotype, rs738409, rs4374383 and rs910049 in accelerating FPR.

CONCLUSIONS

Most factors accelerating liver fibrosis progression in chronic hepatitis C are unmodifiable.

The Milieu Intérieur study - an integrative approach for study of human immunological variance

The Milieu Intérieur Consortium has established a 1000-person healthy population-based study (stratified according to sex and age), creating an unparalleled opportunity for assessing the determinants of human immunologic variance. Herein, we define the criteria utilized for participant enrollment, and highlight the key data that were collected for correlative studies. In this report, we analyzed biological correlates of sex, age, smoking-habits, metabolic score and CMV infection. We characterized and identified unique risk factors among healthy donors, as compared to studies that have focused on the general population or disease cohorts. Finally, we highlight sex-bias in the thresholds used for metabolic score determination and recommend a deeper examination of current guidelines. In sum, our clinical design, standardized sample collection strategies, and epidemiological data analyses have established the foundation for defining variability within human immune responses.

Exploring the occurrence of classic selective sweeps in humans using whole-genome sequencing data sets

Genome-wide scans for selection have identified multiple regions of the human genome as being targeted by positive selection. However, only a small proportion has been replicated across studies, and the prevalence of positive selection as a mechanism of adaptive change in humans remains controversial. Here we explore the power of two haplotype-based statistics--the integrated haplotype score (iHS) and the Derived Intraallelic Nucleotide Diversity (DIND) test--in the context of next-generation sequencing data, and evaluate their robustness to demography and other selection modes. We show that these statistics are both powerful for the detection of recent positive selection, regardless of population history, and robust to variation in coverage, with DIND being insensitive to very low coverage. We apply these statistics to whole-genome sequence data sets from the 1000 Genomes Project and Complete Genomics. We found that putative targets of selection were highly significantly enriched in genic and nonsynonymous single nucleotide polymorphisms, and that DIND was more powerful than iHS in the context of small sample sizes, low-quality genotype calling, or poor coverage. As we excluded genomic confounders and alternative selection models, such as background selection, the observed enrichment attests to the action of recent, strong positive selection. Further support to the adaptive significance of these genomic regions came from their enrichment in functional variants detected by genome-wide association studies, informing the relationship between past selection and current benign and disease-related phenotypic variation. Our results indicate that hard sweeps targeting low-frequency standing variation have played a moderate, albeit significant, role in recent human evolution.

Heterogeneous pattern of selective pressure for PRRT2 in human populations, but no association with autism spectrum disorders

Inherited and de novo genomic imbalances at chromosome 16p11.2 are associated with autism spectrum disorders (ASD), but the causative genes remain unknown. Among the genes located in this region, PRRT2 codes for a member of the synaptic SNARE complex that allows the release of synaptic vesicles. PRRT2 is a candidate gene for ASD since homozygote mutations are associated with intellectual disability and heterozygote mutations cause benign infantile seizures, paroxysmal dyskinesia, or hemiplegic migraine. Here, we explored the contribution of PRRT2 mutations in ASD by screening its coding part in a large sample of 1578 individuals including 431 individuals with ASD, 186 controls and 961 individuals from the human genome Diversity Panel. We detected 24 nonsynonymous variants, 1 frameshift (A217PfsX8) and 1 in-frame deletion of 6 bp (p.A361_P362del). The frameshift mutation was observed in a control with no history of neurological or psychiatric disorders. The p.A361_P362del was observed in two individuals with autism from sub-Saharan African origin. Overall, the frequency of PRRT2 deleterious variants was not different between individuals with ASD and controls. Remarkably, PRRT2 displays a highly significant excess of nonsynonymous (pN) vs synonymous (pS) mutations in Asia (pN/pS = 4.85) and Europe (pN/pS = 1.62) compared with Africa (pN/pS = 0.26; Asia vs Africa: P = 0.000087; Europe vs Africa P = 0.00035; Europe vs Asia P = P = 0.084). We also showed that whole genome amplification performed through rolling cycle amplification could artificially introduce the A217PfsX8 mutation indicating that this technology should not be performed prior to PRRT2 mutation screening. In summary, our results do not support a role for PRRT2 coding sequence variants in ASD, but provide an ascertainment of its genetic variability in worldwide populations that should help researchers and clinicians to better investigate the role of PRRT2 in human diseases.

A genomic portrait of the genetic architecture and regulatory impact of microRNA expression in response to infection

MicroRNAs (miRNAs) are critical regulators of gene expression, and their role in a wide variety of biological processes, including host antimicrobial defense, is increasingly well described. Consistent with their diverse functional effects, miRNA expression is highly context dependent and shows marked changes upon cellular activation. However, the genetic control of miRNA expression in response to external stimuli and the impact of such perturbations on miRNA-mediated regulatory networks at the population level remain to be determined. Here we assessed changes in miRNA expression upon Mycobacterium tuberculosis infection and mapped expression quantitative trait loci (eQTL) in dendritic cells from a panel of healthy individuals. Genome-wide expression profiling revealed that ∼40% of miRNAs are differentially expressed upon infection. We find that the expression of 3% of miRNAs is controlled by proximate genetic factors, which are enriched in a promoter-specific histone modification associated with active transcription. Notably, we identify two infection-specific response eQTLs, for miR-326 and miR-1260, providing an initial assessment of the impact of genotype-environment interactions on miRNA molecular phenotypes. Furthermore, we show that infection coincides with a marked remodeling of the genome-wide relationships between miRNA and mRNA expression levels. This observation, supplemented by experimental data using the model of miR-29a, sheds light on the role of a set of miRNAs in cellular responses to infection. Collectively, this study increases our understanding of the genetic architecture of miRNA expression in response to infection, and highlights the wide-reaching impact of altering miRNA expression on the transcriptional landscape of a cell.

A 1,100-year-old founder effect mutation in IL12B gene is responsible for Mendelian susceptibility to mycobacterial disease in Tunisian patients

Mendelian susceptibility to mycobacterial disease (MSMD) is a rare disorder predisposing apparently healthy individuals to infections caused by weakly virulent mycobacteria such as bacille Calmette-Guerin (BCG), environmental mycobacteria, and poorly virulent Salmonella strains. IL-12p40 deficiency is the first reported human disease due to a cytokine gene defect and is one of the deficiencies that cause MSMD. Nine mutant alleles only have been identified in the IL12B gene, and three of them are recurrent mutations due to a founder effect in specific populations. IL-12p40 deficiency has been identified especially in countries where consanguinity is high and where BCG vaccination at birth is universal. We investigated, in such settings, the clinical, cellular, and molecular features of six IL-12p40-deficient Tunisian patients having the same mutation in IL12B gene (c.298_305del). We found that this mutation is inherited as a common founder mutation arousing ~1,100 years ago. This finding facilitates the development of a preventive approach by genetic counseling and prenatal diagnosis especially in affected families.

Human genetic data reveal contrasting demographic patterns between sedentary and nomadic populations that predate the emergence of farming

Demographic changes are known to leave footprints on genetic polymorphism. Together with the increased availability of large polymorphism data sets, coalescent-based methods allow inferring the past demography of populations from their present-day patterns of genetic diversity. Here, we analyzed both nuclear (20 noncoding regions) and mitochondrial (HVS-I) resequencing data to infer the demographic history of 66 African and Eurasian human populations presenting contrasting lifestyles (nomadic hunter-gatherers, nomadic herders, and sedentary farmers). This allowed us to investigate the relationship between lifestyle and demography and to address the long-standing debate about the chronology of demographic expansions and the Neolithic transition. In Africa, we inferred expansion events for farmers, but constant population sizes or contraction events for hunter-gatherers. In Eurasia, we inferred higher expansion rates for farmers than herders with HVS-I data, except in Central Asia and Korea. Although isolation and admixture processes could have impacted our demographic inferences, these processes alone seem unlikely to explain the contrasted demographic histories inferred in populations with different lifestyles. The small expansion rates or constant population sizes inferred for herders and hunter-gatherers may thus result from constraints linked to nomadism. However, autosomal data revealed contraction events for two sedentary populations in Eurasia, which may be caused by founder effects. Finally, the inferred expansions likely predated the emergence of agriculture and herding. This suggests that human populations could have started to expand in Paleolithic times, and that strong Paleolithic expansions in some populations may have ultimately favored their shift toward agriculture during the Neolithic.

Different selective pressures shape the evolution of Toll-like receptors in human and African great ape populations

The study of the genetic and selective landscape of immunity genes across primates can provide insight into the existing differences in susceptibility to infection observed between human and non-human primates. Here, we explored how selection has driven the evolution of a key family of innate immunity receptors, the Toll-like receptors (TLRs), in African great ape species. We sequenced the 10 TLRs in various populations of chimpanzees and gorillas, and analysed these data jointly with a human data set. We found that purifying selection has been more pervasive in great apes than in humans. Furthermore, in chimpanzees and gorillas, purifying selection has targeted TLRs irrespectively of whether they are endosomal or cell surface, in contrast to humans where strong selective constraints are restricted to endosomal TLRs. These observations suggest important differences in the relative importance of TLR-mediated pathogen sensing, such as that of recognition of flagellated bacteria by TLR5, between humans and great apes. Lastly, we used a population genetics-phylogenetics method that jointly analyses polymorphism and divergence data to detect fine-scale variation in selection pressures at specific codons within TLR genes. We identified different codons at different TLRs as being under positive selection in each species, highlighting that functional variation at these genes has conferred a selective advantage in immunity to infection to specific primate species. Overall, this study showed that the degree of selection driving the evolution of TLRs has largely differed between human and non-human primates, increasing our knowledge on their respective biological contribution to host defence in the natural setting.

Ribosomal protein SA haploinsufficiency in humans with isolated congenital asplenia

Isolated congenital asplenia (ICA) is characterized by the absence of a spleen at birth in individuals with no other developmental defects. The patients are prone to life-threatening bacterial infections. The unbiased analysis of exomes revealed heterozygous mutations in RPSA in 18 patients from eight kindreds, corresponding to more than half the patients and over one-third of the kindreds studied. The clinical penetrance in these kindreds is complete. Expression studies indicated that the mutations carried by the patients-a nonsense mutation, a frameshift duplication, and five different missense mutations-cause autosomal dominant ICA by haploinsufficiency. RPSA encodes ribosomal protein SA, a component of the small subunit of the ribosome. This discovery establishes an essential role for RPSA in human spleen development.

Positive selection in the chromosome 16 VKORC1 genomic region has contributed to the variability of anticoagulant response in humans

VKORC1 (vitamin K epoxide reductase complex subunit 1, 16p11.2) is the main genetic determinant of human response to oral anticoagulants of antivitamin K type (AVK). This gene was recently suggested to be a putative target of positive selection in East Asian populations. In this study, we genotyped the HGDP-CEPH Panel for six VKORC1 SNPs and downloaded chromosome 16 genotypes from the HGDP-CEPH database in order to characterize the geographic distribution of footprints of positive selection within and around this locus. A unique VKORC1 haplotype carrying the promoter mutation associated with AVK sensitivity showed especially high frequencies in all the 17 HGDP-CEPH East Asian population samples. VKORC1 and 24 neighboring genes were found to lie in a 505 kb region of strong linkage disequilibrium in these populations. Patterns of allele frequency differentiation and haplotype structure suggest that this genomic region has been submitted to a near complete selective sweep in all East Asian populations and only in this geographic area. The most extreme scores of the different selection tests are found within a smaller 45 kb region that contains VKORC1 and three other genes (BCKDK, MYST1 (KAT8), and PRSS8) with different functions. Because of the strong linkage disequilibrium, it is not possible to determine if VKORC1 or one of the three other genes is the target of this strong positive selection that could explain present-day differences among human populations in AVK dose requirement. Our results show that the extended region surrounding a presumable single target of positive selection should be analyzed for genetic variation in a wide range of genetically diverse populations in order to account for other neighboring and confounding selective events and the hitchhiking effect.