optiCall: a robust genotype-calling algorithm for rare, low-frequency and common variants

Polygenic Risk Associations with Clinical Characteristics and Recurrence of Dupuytren Disease

BACKGROUND

Dupuytren disease (DD) is a common complex trait, with varying severity and incompletely understood cause. Genome-wide association studies (GWAS) have identified risk loci. In this article, we examine whether genetic risk profiles of DD in patients are associated with clinical variation and disease severity and with patient genetic risk profiles of genetically correlated traits, including body mass index (BMI), triglycerides, high-density lipoproteins, type 2 diabetes mellitus, and endophenotypes fasting glucose and glycated hemoglobin.

METHODS

The authors used a well-characterized cohort of 1461 DD patients with available phenotypic and genetic data. Phenotype data include age at onset, recurrence, and family history of disease. Polygenic risk scores (PRSs) of DD, BMI, triglycerides, high-density lipoprotein, type 2 diabetes, fasting glucose, and hemoglobin A1c using various significance thresholds were calculated with PRSice using the most recent GWAS summary statistics. Control data from LifeLines were used to determine P value cutoffs for PRS generation explaining most variance.

RESULTS

The PRS for DD was significantly associated with a positive family history for DD, age at onset, disease onset before the age of 50, and recurrence. We also found a significant negative correlation between the PRSs for DD and BMI.

CONCLUSIONS

Although GWAS studies of DD are designed to identify genetic risk factors distinguishing case/control status, we show that the genetic risk profile for DD also explains part of its clinical variation and disease severity. The PRS may therefore aid in accurate prognostication, choosing initial treatment and in personalized medicine in the future.

CLINICAL QUESTION/LEVEL OF EVIDENCE

Risk, III.

Genetic and nongenetic drivers of platelet reactivity in healthy Tanzanian individuals

BACKGROUND

Platelets play a key role in hemostasis, inflammation, and cardiovascular diseases. Platelet reactivity is highly variable between individuals. The drivers of this variability in populations from Sub-Saharan Africa remain largely unknown.

OBJECTIVES

We aimed to investigate the nongenetic and genetic determinants of platelet reactivity in healthy adults living in a rapidly urbanizing area in Northern Tanzania.

METHODS

Platelet activation and reactivity were measured by platelet P-selectin expression and the binding of fibrinogen in unstimulated blood and after ex vivo stimulation with adenosine diphosphate and PAR-1 and PAR-4 ligands. We then analyzed the associations of platelet parameters with host genetic and nongenetic factors, environmental factors, plasma inflammatory markers, and plasma metabolites.

RESULTS

Only a few associations were found between platelet reactivity parameters and plasma inflammatory markers and nongenetic host and environmental factors. In contrast, untargeted plasma metabolomics revealed a large number of associations with food-derived metabolites, including phytochemicals that were previously reported to inhibit platelet reactivity. Genome-wide single-nucleotide polymorphism genotyping identified 2 novel single-nucleotide polymorphisms (rs903650 and rs4789332) that were associated with platelet reactivity at the genome-wide level (P < 5 × 10-8) as well as a number of variants in the PAR4 gene (F2RL3) that were associated with PAR4-induced reactivity.

CONCLUSION

Our study uncovered factors that determine variation in platelet reactivity in a population in East Africa that is rapidly transitioning to an urban lifestyle, including the importance of genetic ancestry and the gradual abandoning of the traditional East African diet.

Sex-biased genetic regulation of inflammatory proteins in the Dutch population

BACKGROUND

Significant differences in immune responses, prevalence or susceptibility of diseases and treatment responses have been described between males and females. Despite this, sex-differentiation analysis of the genetic architecture of inflammatory proteins is largely unexplored. We performed sex-stratified meta-analysis after protein quantitative trait loci (pQTL) mapping using inflammatory biomarkers profiled using targeted proteomics (Olink inflammatory panel) of two population-based cohorts of Europeans.

RESULTS

Even though, around 67% of the pQTLs demonstrated shared effect between sexes, colocalization analysis identified two loci in the males (LINC01135 and ITGAV) and three loci (CNOT10, SRD5A2, and LILRB5) in the females with evidence of sex-dependent modulation by pQTL variants. Furthermore, we identified pathways with relevant functions in the sex-biased pQTL variants. We also showed through cross-validation that the sex-specific pQTLs are linked with sex-specific phenotypic traits.

CONCLUSION

Our study demonstrates the relevance of genetic sex-stratified analysis in the context of genetic dissection of protein abundances among individuals and reveals that, sex-specific pQTLs might mediate sex-linked phenotypes. Identification of sex-specific pQTLs associated with sex-biased diseases can help realize the promise of individualized treatment.

The effect of leptin on trained innate immunity and on systemic inflammation in subjects with obesity

Leptin is associated with cardiometabolic complications of obesity, such as metabolic syndrome and atherosclerosis. In obese men, the presence of metabolic syndrome is associated with higher circulating leptin and interleukin (IL)-6 concentrations and increased monocyte cytokine production capacity. Here, we investigated the effects of leptin on monocyte function and systemic inflammatory markers in obese individuals. We specifically explored whether leptin can induce long-term changes in innate immune function by inducing innate immune memory (also called trained immunity). We exposed human primary monocytes for 24 h to relevant leptin concentrations in vitro and measured cytokine production. In addition, after removing leptin, we incubated monocytes for 5 d in culture medium, and we restimulated them on day 6 to assess cytokine production capacity, phagocytosis, and foam cell formation. Direct stimulation with leptin did not induce cytokine production, but exposure to 50 ng/mL leptin augmented lipopolysaccharide- and R848-induced tumor necrosis factor α (TNF-α) production after 1 wk. In a separate in vivo study in a cohort of 302 obese subjects (body mass index [BMI] >27 kg/m2, 55 to 81 yr), we measured circulating leptin, inflammatory markers, and cytokine production upon ex vivo stimulation of isolated peripheral blood mononuclear cells. Circulating leptin concentrations positively correlated with circulating IL-1β and IL-6, which was more pronounced in men than in women. Four single nucleotide polymorphisms in the leptin gene influenced circulating IL-6 concentrations in men, suggesting a direct effect of leptin on IL-6. In conclusion, in vitro, leptin does not directly stimulate monocytes to produce cytokines, yet induces long-term monocyte hyperresponsiveness, i.e. trained immunity. In obese subjects, leptin is associated with circulating IL-6 in a sex-dependent manner. The underlying mechanisms of the sex-specific effect of leptin on innate immune cells remain to be further investigated.

Multi-omics analysis of innate and adaptive responses to BCG vaccination reveals epigenetic cell states that predict trained immunity

Immune responses are tightly regulated yet highly variable between individuals. To investigate human population variation of trained immunity, we immunized healthy individuals with Bacillus Calmette-Guérin (BCG). This live-attenuated vaccine induces not only an adaptive immune response against tuberculosis but also triggers innate immune activation and memory that are indicative of trained immunity. We established personal immune profiles and chromatin accessibility maps over a 90-day time course of BCG vaccination in 323 individuals. Our analysis uncovered genetic and epigenetic predictors of baseline immunity and immune response. BCG vaccination enhanced the innate immune response specifically in individuals with a dormant immune state at baseline, rather than providing a general boost of innate immunity. This study advances our understanding of BCG's heterologous immune-stimulatory effects and trained immunity in humans. Furthermore, it highlights the value of epigenetic cell states for connecting immune function with genotype and the environment.

Uncovering myocardial infarction genetic signatures using GWAS exploration in Saudi and European cohorts

Genome-wide association studies (GWAS) have yielded significant insights into the genetic architecture of myocardial infarction (MI), although studies in non-European populations are still lacking. Saudi Arabian cohorts offer an opportunity to discover novel genetic variants impacting disease risk due to a high rate of consanguinity. Genome-wide genotyping (GWG), imputation and GWAS followed by meta-analysis were performed based on two independent Saudi Arabian studies comprising 3950 MI patients and 2324 non-MI controls. Meta-analyses were then performed with these two Saudi MI studies and the CardioGRAMplusC4D and UK BioBank GWAS as controls. Meta-analyses of the two Saudi MI studies resulted in 17 SNPs with genome-wide significance. Meta-analyses of all 4 studies revealed 66 loci with genome-wide significance levels of p < 5 × 10-8. All of these variants, except rs2764203, have previously been reported as MI-associated loci or to have high linkage disequilibrium with known loci. One SNP association in Shisa family member 5 (SHISA5) (rs11707229) was evident at a much higher frequency in the Saudi MI populations (> 12% MAF). In conclusion, our results replicated many MI associations, whereas in Saudi-only GWAS (meta-analyses), several new loci were implicated that require future validation and functional analyses.

A single-cell view on host immune transcriptional response to in vivo BCG-induced trained immunity

Bacillus Calmette-Guérin (BCG) vaccination is a prototype model for the study of trained immunity (TI) in humans, and results in a more effective response of innate immune cells upon stimulation with heterologous stimuli. Here, we investigate the heterogeneity of TI induction by single-cell RNA sequencing of immune cells collected from 156 samples. We observe that both monocytes and CD8⁺ T cells show heterologous transcriptional responses to lipopolysaccharide, with an active crosstalk between these two cell types. Furthermore, the interferon-γ pathway is crucial in BCG-induced TI, and it is upregulated in functional high responders. Data-driven analyses and functional experiments reveal STAT1 to be one of the important transcription factors for TI shared in all identified monocyte subpopulations. Finally, we report the role of type I interferon-related and neutrophil-related TI transcriptional programs in patients with sepsis. These findings provide comprehensive insights into the importance of monocyte heterogeneity during TI in humans.

Genome-Wide Association Study Reveals CLEC7A and PROM1 as Potential Regulators of Paracoccidioides brasiliensis-Induction of Cytokine Production in Peripheral Blood Mononuclear Cells

Paracoccidioidomycosis (PCM) is a systemic mycosis caused by fungi of the genus Paracoccidioides and the different clinical forms of the disease are associated with the host immune responses. Quantitative trait loci mapping analysis was performed to assess genetic variants associated with mononuclear-cells-derived cytokines induced by P. brasiliensis on 158 individuals. We identified the rs11053595 SNP, which is present in the CLEC7A gene (encodes the Dectin-1 receptor) and the rs62290169 SNP located in the PROM1 gene (encodes CD133) associated with the production of IL-1β and IL-22, respectively. Functionally, the blockade of the dectin-1 receptor abolished the IL-1β production in P. brasiliensis-stimulated PBMCs. Moreover, the rs62290169-GG genotype was associated with higher frequency of CD38+ Th1 cells in PBMCs cultured with P. brasiliensis yeasts. Therefore, our research indicates that the CLEC7A and PROM1 genes are important for the cytokine response induced by P. brasiliensis and may influence the Paracoccidioidomycosis disease outcome.

Genetic regulators of cytokine responses upon BCG vaccination in children from West Africa

Genetic variation is a key factor influencing cytokine production capacity, but which genetic loci regulate cytokine production before and after vaccination, particularly in African population is unknown. Here, we aimed to identify single-nucleotide polymorphisms (SNPs) controlling cytokine responses (cQTLs) after microbial stimulation in infants of West-African ancestry, comprising of low-birth-weight neonates randomized to bacillus Calmette-Guérin (BCG) vaccine-at-birth (intervention) or to the usual delayed BCG (control). Genome-wide cytokine QTL mapping revealed 12 independent cQTLs loci, of which the LINC01082-LINC00917 locus influenced more than half of the cytokine-stimulation pairs assessed. Furthermore, nine distinct cQTLs were found among infants randomized to BCG. Functional validation confirmed that several complement genes affect cytokine response after BCG vaccination. We observed a limited overlap of common cQTLs between the West-African infants and cohorts of Western European individuals. These data reveal strong population-specific genetic effects on cytokine production and may indicate new opportunities for therapeutic intervention and vaccine development in African populations.

Genetic associations with resistance to Meloidogyne enterolobii in guava (Psidium sp.) using cross-genera SNPs and comparative genomics to Eucalyptus highlight evolutionary conservation across the Myrtaceae

Tropical fruit tree species constitute a yet untapped supply of outstanding diversity of taste and nutritional value, barely developed from the genetics standpoint, with scarce or no genomic resources to tackle the challenges arising in modern breeding practice. We generated a de novo genome assembly of the Psidium guajava, the super fruit “apple of the tropics”, and successfully transferred 14,268 SNP probesets from Eucalyptus to Psidium at the nucleotide level, to detect genomic loci linked to resistance to the root knot nematode (RKN) Meloidogyne enterolobii derived from the wild relative P. guineense. Significantly associated loci with resistance across alternative analytical frameworks, were detected at two SNPs on chromosome 3 in a pseudo-assembly of Psidium guajava genome built using a syntenic path approach with the Eucalyptus grandis genome to determine the order and orientation of the contigs. The P. guineense-derived resistance response to RKN and disease onset is conceivably triggered by mineral nutrients and phytohormone homeostasis or signaling with the involvement of the miRNA pathway. Hotspots of mapped resistance quantitative trait loci and functional annotation in the same genomic region of Eucalyptus provide further indirect support to our results, highlighting the evolutionary conservation of genomes across genera of Myrtaceae in the adaptation to pathogens. Marker assisted introgression of the resistance loci mapped should accelerate the development of improved guava cultivars and hybrid rootstocks.

Genetic determinants of fungi-induced ROS production are associated with the risk of invasive pulmonary aspergillosis

Reactive oxygen species (ROS) are an essential component of the host defense against fungal infections. However, little is known about how common genetic variation affects ROS-mediated antifungal host defense. In the present study, we investigated the genetic factors that regulate ROS production capacity in response to the two human fungal pathogens: Candida albicans and Aspergillus fumigatus. We investigated fungal-stimulated ROS production by immune cells isolated from a population-based cohort of approximately 200 healthy individuals (200FG cohort), and mapped ROS-quantitative trait loci (QTLs). We identified several genetic loci that regulate ROS levels (P < 9.99 × 10-6), with some of these loci being pathogen-specific, and others shared between the two fungi. These ROS-QTLs were investigated for their influence on the risk of invasive pulmonary aspergillosis (IPA) in a disease relevant context. We stratified hematopoietic stem-cell transplant (HSCT) recipients based on the donor's SNP genotype and tested their impact on the risk of IPA. We identified rs4685368 as a ROS-QTL locus that was significantly associated with an increased risk of IPA after controlling for patient age and sex, hematological malignancy, type of transplantation, conditioning regimen, acute graft-versus-host-disease grades III-IV, and antifungal prophylaxis. Collectively, this data provides evidence that common genetic variation can influence ROS production capacity, and, importantly, the risk of developing IPA among HSCT recipients. This evidence warrants further research for patient stratification based on the genetic profiling that would allow the identifications of patients at high-risk for an invasive fungal infection, and who would benefit the most from a preventive strategy.

Impact of Human Genetic Variation on C-Reactive Protein Concentrations and Acute Appendicitis

Background

Acute appendicitis is one of the most common abdominal emergencies worldwide. Both environmental and genetic factors contribute to the disease. C-reactive protein (CRP) is an important biomarker in the diagnosis of acute appendicitis. CRP concentrations are significantly affected by genetic variation. However, whether such genetic variation is causally related to appendicitis risk remains unclear. In this study, the causal relationship between single-nucleotide polymorphisms (SNPs) associated with circulating CRP concentrations and the risk and severity of acute appendicitis was investigated.

Methods

CRP concentrations in serum of appendicitis patients (n = 325) were measured. Appendicitis was categorized as complicated/uncomplicated and gangrenous/non-gangrenous. Imputed SNP data (n = 287) were generated. A genome-wide association study (GWAS) on CRP concentrations and appendicitis severity was performed. Intersection and colocalization of the GWAS results were performed with appendicitis and CRP-associated loci from the Pan-UKBB cohort. A functional-genomics approach to prioritize genes was employed.

Results

Thirteen percent of significant CRP quantitative trait loci (QTLs) that were previously identified in a large cohort of healthy individuals were replicated in our small patient cohort. Significant enrichment of CRP-QTLs in association with appendicitis was observed. Among these shared loci, the two top loci at chromosomes 1q41 and 8p23.1 were characterized. The top SNP at chromosome 1q41 is located within the promoter of H2.0 Like Homeobox (HLX) gene, which is involved in blood cell differentiation, and liver and gut organogeneses. The expression of HLX is increased in the appendix of appendicitis patients compared to controls. The locus at 8p23.1 contains multiple genes, including cathepsin B (CTSB), which is overexpressed in appendix tissue from appendicitis patients. The risk allele of the top SNP in this locus also increases CTSB expression in the sigmoid colon of healthy individuals. CTSB is involved in collagen degradation, MHC class II antigen presentation, and neutrophil degranulation.

Conclusions

The results of this study prioritize HLX and CTSB as potential causal genes for appendicitis and suggest a shared genetic mechanism between appendicitis and CRP concentrations.

A genome-wide functional genomics approach uncovers genetic determinants of immune phenotypes in type 1 diabetes

Background:

The large inter-individual variability in immune-cell composition and function determines immune responses in general and susceptibility o immune-mediated diseases in particular. While much has been learned about the genetic variants relevant for type 1 diabetes (T1D), the pathophysiological mechanisms through which these variations exert their effects remain unknown.

Methods:

Blood samples were collected from 243 patients with T1D of Dutch descent. We applied genetic association analysis on >200 immune-cell traits and >100 cytokine production profiles in response to stimuli measured to identify genetic determinants of immune function, and compared the results obtained in T1D to healthy controls.

Results:

Genetic variants that determine susceptibility to T1D significantly affect T cell composition. Specifically, the CCR5+ regulatory T cells associate with T1D through the CCR region, suggesting a shared genetic regulation. Genome-wide quantitative trait loci (QTLs) mapping analysis of immune traits revealed 15 genetic loci that influence immune responses in T1D, including 12 that have never been reported in healthy population studies, implying a disease-specific genetic regulation.

Conclusions:

This study provides new insights into the genetic factors that affect immunological responses in T1D.

Funding:

This work was supported by an ERC starting grant (no. 948207) and a Radboud University Medical Centre Hypatia grant (2018) to YL and an ERC advanced grant (no. 833247) and a Spinoza grant of the Netherlands Association for Scientific Research to MGN CT received funding from the Perspectief Biomarker Development Center Research Programme, which is (partly) financed by the Netherlands Organisation for Scientific Research (NWO). AJ was funded by a grant from the European Foundation for the Study of Diabetes (EFSD/AZ Macrovascular Programme 2015). XC was supported by the China Scholarship Council (201706040081).

Every year around the world, over 100,000 people are diagnosed with type 1 diabetes. This disease develops when the immune system mistakenly destroys the cells that produce a hormone called insulin, leaving affected individuals unable to regulate their blood sugar levels. Type 1 diabetes patients must rely on regular injections of manufactured insulin to survive.

The composition and activity of the human immune system is under genetic control, and people with certain changes in their genes are more susceptible than others to develop type 1 diabetes. Previous studies have identified around 60 locations in the human DNA (known as loci) associated with the condition, but it remains unclear how these loci influence the immune system and whether diabetes will emerge.

Chu, Janssen, Koenen et al. explored how variations in genetic information can influence the composition of the immune system, and the type of molecules it releases to perform its role. To do so, blood samples from 243 individuals of Dutch descent with type 1 diabetes were collected, and genetic associations were investigated.

The results revealed that a major type of immune actors known as T cells are under the control of genetic factors associated with type 1 diabetes susceptibility. For instance, a specific type of T cells showed shared genetic control with type 1 diabetes. In addition, 15 loci were identified that influenced immune responses in the patients. Among those, 12 have never been reported to be involved in immune responses in healthy people, implying that these regions might only regulate the immune system of individuals with type 1 diabetes and other similar disorders. Finally, Chu, Janssen, Koenen et al. propose 11 genes within the identified loci as potential targets for new diabetes medication. These results represent an important resource for researchers exploring the genetic and immune basis of type 1 diabetes, and they could open new avenues for drug development.

Differences in thrombin and plasmin generation potential between East African and Western European adults: The role of genetic and non-genetic factors

BACKGROUND

Geographic variability in coagulation across populations and their determinants are poorly understood.

OBJECTIVE

To compare thrombin (TG) and plasmin (PG) generation parameters between healthy Tanzanian and Dutch individuals, and to study associations with inflammation and different genetic, host and environmental factors.

METHODS

TG and PG parameters were measured in 313 Tanzanians of African descent living in Tanzania and 392 Dutch of European descent living in the Netherlands and related to results of a dietary questionnaire, circulating inflammatory markers, genotyping, and plasma metabolomics.

RESULTS

Tanzanians exhibited an enhanced TG and PG capacity, compared to Dutch participants. A higher proportion of Tanzanians had a TG value in the upper quartile with a PG value in the lower/middle quartile, suggesting a relative pro-coagulant state. Tanzanians also displayed an increased normalized thrombomodulin sensitivity ratio, suggesting reduced sensitivity to protein C. In Tanzanians, PG parameters (lag time and TTP) were associated with seasonality and food-derived plasma metabolites. The Tanzanians had higher concentrations of pro-inflammatory cytokines, which correlated strongly with TG and PG parameters. There was limited overlap in genetic variation associated with TG and PG parameters between the two cohorts. Pathway analysis of genetic variants in the Tanzanian cohort revealed multiple immune pathways that were enriched with TG and PG traits, confirming the importance of co-regulation between coagulation and inflammation.

CONCLUSIONS

Tanzanians have an enhanced TG and PG potential compared to Dutch individuals, which may relate to differences in inflammation, genetics and diet. These observations highlight the importance of better understanding of the geographic variability in coagulation across populations.

Single-Cell RNA Sequencing of Peripheral Blood Mononuclear Cells From Pediatric Coeliac Disease Patients Suggests Potential Pre-Seroconversion Markers

Celiac Disease (CeD) is a complex immune disorder involving villous atrophy in the small intestine that is triggered by gluten intake. Current CeD diagnosis is based on late-stage pathophysiological parameters such as detection of specific antibodies in blood and histochemical detection of villus atrophy and lymphocyte infiltration in intestinal biopsies. To date, no early onset biomarkers are available that would help prevent widespread villous atrophy and severe symptoms and co-morbidities. To search for novel CeD biomarkers, we used single-cell RNA sequencing (scRNAseq) to investigate PBMC samples from 11 children before and after seroconversion for CeD and 10 control individuals matched for age, sex and HLA-genotype. We generated scRNAseq profiles of 9559 cells and identified the expected major cellular lineages. Cell proportions remained stable across the different timepoints and health conditions, but we observed differences in gene expression profiles in specific cell types when comparing patient samples before and after disease development and comparing patients with controls. Based on the time when transcripts were differentially expressed, we could classify the deregulated genes as biomarkers for active CeD or as potential pre-diagnostic markers. Pathway analysis showed that active CeD biomarkers display a transcriptional profile associated with antigen activation in CD4+ T cells, whereas NK cells express a subset of biomarker genes even before CeD diagnosis. Intersection of biomarker genes with CeD-associated genetic risk loci pinpointed genetic factors that might play a role in CeD onset. Investigation of potential cellular interaction pathways of PBMC cell subpopulations highlighted the importance of TNF pathways in CeD. Altogether, our results pinpoint genes and pathways that are altered prior to and during CeD onset, thereby identifying novel potential biomarkers for CeD diagnosis in blood.

A functional genomics approach in Tanzanian population identifies distinct genetic regulators of cytokine production compared to European population

Humans exhibit remarkable interindividual and interpopulation immune response variability upon microbial challenges. Cytokines play a vital role in regulating inflammation and immune responses, but dysregulation of cytokine responses has been implicated in different disease states. Host genetic factors were previously shown to significantly impact cytokine response heterogeneity mainly in European-based studies, but it is unclear whether these findings are transferable to non-European individuals. Here, we aimed to identify genetic variants modulating cytokine responses in healthy adults of East African ancestry from Tanzania. We leveraged both cytokine and genetic data and performed genome-wide cytokine quantitative trait loci (cQTLs) mapping. The results were compared with another cohort of healthy adults of Western European ancestry via direct overlap and functional enrichment analyses. We also performed meta-analyses to identify cQTLs with congruent effect direction in both populations. In the Tanzanians, cQTL mapping identified 80 independent suggestive loci and one genome-wide significant locus (TBC1D22A) at chromosome 22; SNP rs12169244 was associated with IL-1b release after Salmonella enteritidis stimulation. Remarkably, the identified cQTLs varied significantly when compared to the European cohort, and there was a very limited percentage of overlap (1.6% to 1.9%). We further observed ancestry-specific pathways regulating induced cytokine responses, and there was significant enrichment of the interferon pathway specifically in the Tanzanians. Furthermore, contrary to the Europeans, genetic variants in the TLR10-TLR1-TLR6 locus showed no effect on cytokine response. Our data reveal both ancestry-specific effects of genetic variants and pathways on cytokine response heterogeneity, hence arguing for the importance of initiatives to include diverse populations into genomics research.

Mitochondrial Genome Study Identifies Association Between Primary Open-Angle Glaucoma and Variants in MT-CYB, MT-ND4 Genes and Haplogroups

Background and purpose: Primary open-angle glaucoma (POAG) is an optic neuropathy characterized by death of retinal ganglion cells and atrophy of the optic nerve head. The susceptibility of the optic nerve to damage has been shown to be mediated by mitochondrial dysfunction. In this study, we aimed to determine a possible association between mitochondrial SNPs or haplogroups and POAG.

Methods: Mitochondrial DNA single nucleotide polymorphisms (mtSNPs) were genotyped using the Illumina Infinium Global Screening Array-24 (GSA) 700K array set. Genetic analyses were performed in a POAG case-control study involving the cohorts, Groningen Longitudinal Glaucoma Study-Lifelines Cohort Study and Amsterdam Glaucoma Study, including 721 patients and 1951 controls in total. We excluded samples not passing quality control for nuclear genotypes and samples with low call rate for mitochondrial variation. The mitochondrial variants were analyzed both as SNPs and haplogroups. These were determined with the bioinformatics software HaploGrep, and logistic regression analysis was used for the association, as well as for SNPs.

Results: Meta-analysis of the results from both cohorts revealed a significant association between POAG and the allele A of rs2853496 [odds ratio (OR) = 0.64; p = 0.006] within the MT-ND4 gene, and for the T allele of rs35788393 (OR = 0.75; p = 0.041) located in the MT-CYB gene. In the mitochondrial haplogroup analysis, the most significant p-value was reached by haplogroup K (p = 1.2 × 10⁻⁰⁵), which increases the risk of POAG with an OR of 5.8 (95% CI 2.7–13.1).

Conclusion: We identified an association between POAG and polymorphisms in the mitochondrial genes MT-ND4 (rs2853496) and MT-CYB (rs35788393), and with haplogroup K. The present study provides further evidence that mitochondrial genome variations are implicated in POAG. Further genetic and functional studies are required to substantiate the association between mitochondrial gene polymorphisms and POAG and to define the pathophysiological mechanisms of mitochondrial dysfunction in glaucoma.

An integrative genomics approach identifies KDM4 as a modulator of trained immunity

Innate immune cells are able to build memory characteristics via a process termed “trained immunity.” Host factors that influence the magnitude of the individual trained immunity response remain largely unknown. Using an integrative genomics approach, our study aimed to prioritize and understand the role of specific genes in trained immunity responses. In vitro‐induced trained immunity responses were assessed in two independent population‐based cohorts of healthy individuals, the 300 Bacillus Calmette‐Guérin (300BCG; n = 267) and 200 Functional Genomics (200FG; n = 110) cohorts from the Human Functional Genomics Project. Genetic loci that influence cytokine responses upon trained immunity were identified by conducting a meta‐analysis of QTLs identified in the 300BCG and 200FG cohorts. From the identified QTL loci, we functionally validated the role of PI3K‐Akt signaling pathway and two genes that belong to the family of Siglec receptors (Siglec‐5 and Siglec‐14). Furthermore, we identified the H3K9 histone demethylases of the KDM4 family as major regulators of trained immunity responses. These data pinpoint an important role of metabolic and epigenetic processes in the regulation of trained immunity responses, and these findings may open new avenues for vaccine design and therapeutic interventions.

Single-cell transcriptomic profiles reveal changes associated with BCG-induced trained immunity and protective effects in circulating monocytes

Bacillus Calmette-Guérin (BCG) vaccine is one of the most widely used vaccines worldwide. In addition to protection against tuberculosis, BCG confers a degree of non-specific protection against other infections by enhancing secondary immune responses to heterologous pathogens, termed "trained immunity." To better understand BCG-induced immune reprogramming, we perform single-cell transcriptomic measurements before and after BCG vaccination using secondary immune stimulation with bacterial lipopolysaccharide (LPS). We find that BCG reduces systemic inflammation and identify 75 genes with altered LPS responses, including inflammatory mediators such as CCL3 and CCL4 that have a heightened response. Co-expression analysis reveals that gene modules containing these cytokines lose coordination after BCG. Other modules exhibit increased coordination, including several humanin nuclear isoforms that we confirm induce trained immunity in vitro. Our results link in vivo BCG administration to single-cell transcriptomic changes, validated in human genetics experiments, and highlight genes that are putatively responsible for non-specific protective effects of BCG.

Integration of metabolomics, genomics, and immune phenotypes reveals the causal roles of metabolites in disease

BACKGROUND

Recent studies highlight the role of metabolites in immune diseases, but it remains unknown how much of this effect is driven by genetic and non-genetic host factors.

RESULT

We systematically investigate circulating metabolites in a cohort of 500 healthy subjects (500FG) in whom immune function and activity are deeply measured and whose genetics are profiled. Our data reveal that several major metabolic pathways, including the alanine/glutamate pathway and the arachidonic acid pathway, have a strong impact on cytokine production in response to ex vivo stimulation. We also examine the genetic regulation of metabolites associated with immune phenotypes through genome-wide association analysis and identify 29 significant loci, including eight novel independent loci. Of these, one locus (rs174584-FADS2) associated with arachidonic acid metabolism is causally associated with Crohn's disease, suggesting it is a potential therapeutic target.

CONCLUSION

This study provides a comprehensive map of the integration between the blood metabolome and immune phenotypes, reveals novel genetic factors that regulate blood metabolite concentrations, and proposes an integrative approach for identifying new disease treatment targets.

oxLDL-Induced Trained Immunity Is Dependent on Mitochondrial Metabolic Reprogramming

Following brief exposure to endogenous atherogenic particles, such as oxidized low-density lipoprotein (oxLDL), monocytes/macrophages can adopt a long-term pro-inflammatory phenotype, which is called trained immunity. This mechanism might contribute to the chronic low-grade inflammation that characterizes atherosclerosis. In this study, we aim to elucidate immunometabolic pathways that drive oxLDL-induced trained immunity. Primary isolated human monocytes were exposed to oxLDL for 24 h, and after five days stimulated with LPS to measure the cytokine production capacity. RNA-sequencing revealed broad increases in genes enriched in mitochondrial pathways after 24 h of oxLDL exposure. Further omics profiling of oxLDL-trained macrophages via intracellular metabolomics showed an enrichment for tricarboxylic acid (TCA) cycle metabolites. Single cell analysis revealed that oxLDL-trained macrophages contain larger mitochondria, potentially likely linked to increased oxidative phosphorylation (OXPHOS) activity. Co-incubation with pharmacological blockers of OXPHOS inhibited oxLDL-induced trained immunity. The relevance of OXPHOS was confirmed in a cohort of 243 healthy subjects showing that genetic variation in genes coding for enzymes relevant to OXPHOS correlated with the capacity of monocytes to be trained with oxLDL. Interestingly, OXPHOS appears to play an important role in the increased cytokine hyperresponsiveness by oxLDL-trained macrophages. The TCA-cycle can also be fuelled by glutamine and free fatty acids, and pharmacological blockade of these pathways could prevent oxLDL-induced trained immunity. This study demonstrates that the mitochondria of oxLDL-trained macrophages undergo changes to their function and form with OXPHOS being an important mechanism for trained immunity, which could unveil novel pharmacological targets to prevent atherogenesis.

The role of sirtuin 1 on the induction of trained immunity

Sirtuin 1 (SIRT1) has been described to modify immune responses by modulation of gene transcription. As transcriptional reprogramming is the molecular substrate of trained immunity, a de facto innate immune memory, we investigated the role of SIRT1 in the induction of trained immunity. We identified various SIRT1 genetic single nucleotide polymorphisms affecting innate and adaptive cytokine production of human peripheral blood mononuclear cells (PBMCs) in response to various stimuli on the one hand, and in vitro induction of trained immunity on the other hand. Furthermore, inhibition of SIRT1 upregulated pro-inflammatory innate cytokine production upon stimulation of PBMCs. However, inhibition of SIRT1 in vitro had no effect on cytokine responses upon induction of trained immunity, while activation of SIRT1 mildly modified trained immunity responses. In conclusion, SIRT1 modifies innate cytokine production by PBMCs in response to various microbes, but has only a secondary role for BCG and β-glucan-induced trained immunity responses.

Glutathione Metabolism Contributes to the Induction of Trained Immunity

The innate immune system displays heterologous memory characteristics, which are characterized by stronger responses to a secondary challenge. This phenomenon termed trained immunity relies on epigenetic and metabolic rewiring of innate immune cells. As reactive oxygen species (ROS) production has been associated with the trained immunity phenotype, we hypothesized that the increased ROS levels and the main intracellular redox molecule glutathione play a role in the induction of trained immunity. Here we show that pharmacological inhibition of ROS in an in vitro model of trained immunity did not influence cell responsiveness; the modulation of glutathione levels reduced pro-inflammatory cytokine production in human monocytes. Single nucleotide polymorphisms (SNPs) in genes involved in glutathione metabolism were found to be associated with changes in pro-inflammatory cytokine production capacity upon trained immunity. Also, plasma glutathione concentrations were positively associated with ex vivo IL-1β production, a biomarker of trained immunity, produced by monocytes of BCG-vaccinated individuals. In conclusion, glutathione metabolism is involved in the induction of trained immunity, and future studies are warranted to explore its functional consequences in human diseases.

Urban living in healthy Tanzanians is associated with an inflammatory status driven by dietary and metabolic changes

Sub-Saharan Africa currently experiences an unprecedented wave of urbanization, which has important consequences for health and disease patterns. This study aimed to investigate and integrate the immune and metabolic consequences of rural or urban lifestyles and the role of nutritional changes associated with urban living. In a cohort of 323 healthy Tanzanians, urban as compared to rural living was associated with a pro-inflammatory immune phenotype, both at the transcript and protein levels. We identified different food-derived and endogenous circulating metabolites accounting for these differences. Serum from urban dwellers induced reprogramming of innate immune cells with higher tumor necrosis factor production upon microbial re-stimulation in an in vitro model of trained immunity. These data demonstrate important shifts toward an inflammatory phenotype associated with an urban lifestyle and provide new insights into the underlying dietary and metabolic factors, which may affect disease epidemiology in sub-Sahara African countries.

Transethnic analysis of the human leukocyte antigen region for ulcerative colitis reveals not only shared but also ethnicity-specific disease associations

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.

de Bree LC, de Mast Q, van de Veerdonk FL, Joosten LA, Li Y, Netea MG, Xu C. The genetic risk for COVID-19 severity is associated with defective innate immune responses.. [DOI: 10.1101/2020.11.10.20229203]

Recent genome-wide association studies (GWASs) of COVID-19 patients of European ancestry have identified genetic loci significantly associated with disease severity (1). Here, we employed the detailed clinical, immunological and multi-omics dataset of the Human Functional Genomics Projects (HFGP) to explore the physiological significance of the host genetic variants that influence susceptibility to severe COVID-19. A genomics investigation intersected with functional characterization of individuals with high genetic risk for severe COVID-19 susceptibility identified several major patterns: i. a large impact of genetically determined innate immune responses in COVID-19, with increased susceptibility for severe disease in individuals with defective monocyte-derived cytokine production; ii. genetic susceptibility related to ABO blood groups is probably mediated through the von Willebrand factor (VWF) and endothelial dysfunction. We further validated these identified associations at transcript and protein levels by using independent disease cohorts. These insights allow a physiological understanding of genetic susceptibility to severe COVID-19, and indicate pathways that could be targeted for prevention and therapy.

One Sentence summary

In this study, we explore the physiological significance of the genetic variants associated with COVID-19 severity using detailed clinical, immunological and multi-omics data from large cohorts. Our findings allow a physiological understanding of genetic susceptibility to severe COVID-19, and indicate pathways that could be targeted for prevention and therapy.

The influence of human genetic variation on early transcriptional responses and protective immunity following immunization with Rotarix vaccine in infants in Ho Chi Minh City in Vietnam: A study protocol for an open single-arm interventional trial

Background: Rotavirus (RoV) remains the leading cause of acute gastroenteritis in infants and children aged under five years in both high- and low-middle-income countries (LMICs). In LMICs, RoV infections are associated with substantial mortality. Two RoV vaccines (Rotarix and Rotateq) are widely available for use in infants, both of which have been shown to be highly efficacious in Europe and North America. However, for unknown reasons, these RoV vaccines have markedly lower efficacy in LMICs. We hypothesize that poor RoV vaccine efficacy across in certain regions may be associated with genetic heritability or gene expression in the human host. Methods/design: We designed an open-label single-arm interventional trial with the Rotarix RoV vaccine to identify genetic and transcriptomic markers associated with generating a protective immune response against RoV. Overall, 1,000 infants will be recruited prior to Expanded Program on Immunization (EPI) vaccinations at two months of age and vaccinated with oral Rotarix vaccine at two and three months, after which the infants will be followed-up for diarrheal disease until 18 months of age. Blood sampling for genetics, transcriptomics, and immunological analysis will be conducted before each Rotarix vaccination, 2-3 days post-vaccination, and at each follow-up visit (i.e. 6, 12 and 18 months of age). Stool samples will be collected during each diarrheal episode to identify RoV infection. The primary outcome will be Rotarix vaccine failure events (i.e. symptomatic RoV infection despite vaccination), secondary outcomes will be antibody responses and genotypic characterization of the infection virus in Rotarix failure events. Discussion: This study will be the largest and best powered study of its kind to be conducted to date in infants, and will be critical for our understanding of RoV immunity, human genetics in the Vietnam population, and mechanisms determining RoV vaccine-mediated protection. Registration: ClinicalTrials.gov, ID: NCT03587389. Registered on 16 July 2018.

Benefits and limitations of genome-wide association studies

Genome-wide association studies (GWAS) involve testing genetic variants across the genomes of many individuals to identify genotype-phenotype associations. GWAS have revolutionized the field of complex disease genetics over the past decade, providing numerous compelling associations for human complex traits and diseases. Despite clear successes in identifying novel disease susceptibility genes and biological pathways and in translating these findings into clinical care, GWAS have not been without controversy. Prominent criticisms include concerns that GWAS will eventually implicate the entire genome in disease predisposition and that most association signals reflect variants and genes with no direct biological relevance to disease. In this Review, we comprehensively assess the benefits and limitations of GWAS in human populations and discuss the relevance of performing more GWAS.

Joint sequencing of human and pathogen genomes reveals the genetics of pneumococcal meningitis

Streptococcus pneumoniae is a common nasopharyngeal colonizer, but can also cause life-threatening invasive diseases such as empyema, bacteremia and meningitis. Genetic variation of host and pathogen is known to play a role in invasive pneumococcal disease, though to what extent is unknown. In a genome-wide association study of human and pathogen we show that human variation explains almost half of variation in susceptibility to pneumococcal meningitis and one-third of variation in severity, identifying variants in CCDC33 associated with susceptibility. Pneumococcal genetic variation explains a large amount of invasive potential (70%), but has no effect on severity. Serotype alone is insufficient to explain invasiveness, suggesting other pneumococcal factors are involved in progression to invasive disease. We identify pneumococcal genes involved in invasiveness including pspC and zmpD, and perform a human-bacteria interaction analysis. These genes are potential candidates for the development of more broadly-acting pneumococcal vaccines.

Latent autoimmunity across disease-specific boundaries in at-risk first-degree relatives of SLE and RA patients

BACKGROUND

Autoimmune disease prevention requires tools to assess an individual's risk of developing a specific disease. One tool is disease-associated autoantibodies, which accumulate in an asymptomatic preclinical period. However, patients sometimes exhibit autoantibodies associated with a different disease classification. When and how these alternative autoantibodies first appear remain unknown. This cross-sectional study characterizes alternative autoimmunity, and associated genetic and environmental factors, in unaffected first-degree relatives (FDRs) of patients, who exhibit increased future risk for the same disease.

METHODS

Samples (n = 1321) from disease-specific autoantibody-positive (aAb+) systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and type 1 diabetes (T1D) patients; and unaffected aAb+ and autoantibody-negative (aAb-) SLE and RA FDRs were tested for SLE, RA, and T1D aAbs, as well as anti-tissue transglutaminase, anti-cardiolipin and anti-thyroperoxidase. FDR SLE and RA genetic risk scores (GRS) were calculated.

FINDINGS

Alternative autoimmunity occurred in SLE patients (56%) and FDRs (57·4%), RA patients (32·6%) and FDRs (34·8%), and T1D patients (43%). Expanded autoimmunity, defined as autoantibodies spanning at least two other diseases, occurred in 18·5% of SLE patients, 16·4% of SLE FDRs, 7·8% of RA patients, 5·3% of RA FDRs, and 10·8% of T1D patients. SLE FDRs were more likely to have alternative (odds ratio [OR] 2·44) and expanded (OR 3·27) autoimmunity than RA FDRs. Alternative and expanded autoimmunity were associated with several environmental exposures. Alternative autoimmunity was associated with a higher RA GRS in RA FDRs (OR 1·41), and a higher SLE GRS in aAb+ RA FDRs (OR 1·87), but not in SLE FDRs.

INTERPRETATION

Autoimmunity commonly crosses disease-specific boundaries in systemic (RA, SLE) and organ-specific (T1D) autoimmune diseases. Alternative autoimmunity is more common in SLE FDRs than RA FDRs, and is influenced by genetic and environmental factors. These findings have substantial implications for preclinical disease pathogenesis and autoimmune disease prevention studies. FUND: NIH U01AI101981, R01AR051394, U19AI082714, P30AR053483, P30GM103510, U54GM104938, U01AI101934, R01AI024717, U01AI130830, I01BX001834, & U01HG008666.

The 1000IBD project: multi-omics data of 1000 inflammatory bowel disease patients; data release 1

BACKGROUND

Inflammatory bowel disease (IBD) is a chronic complex disease of the gastrointestinal tract. Patients with IBD can experience a wide range of symptoms, but the pathophysiological mechanisms that cause these individual differences in clinical presentation remain largely unknown. In consequence, IBD is currently classified into subtypes using clinical characteristics. If we are to develop a more targeted treatment approach, molecular subtypes of IBD need to be discovered that can be used as new drug targets. To achieve this, we need multiple layers of molecular data generated from the same IBD patients.

CONSTRUCTION AND CONTENT

We initiated the 1000IBD project ( https://1000ibd.org ) to prospectively follow more than 1000 IBD patients from the Northern provinces of the Netherlands. For these patients, we have collected a uniquely large number of phenotypes and generated multi-omics profiles. To date, 1215 participants have been enrolled in the project and enrolment is on-going. Phenotype data collected for these participants includes information on dietary and environmental factors, drug responses and adverse drug events. Genome information has been generated using genotyping (ImmunoChip, Global Screening Array and HumanExomeChip) and sequencing (whole exome sequencing and targeted resequencing of IBD susceptibility loci), transcriptome information generated using RNA-sequencing of intestinal biopsies and microbiome information generated using both sequencing of the 16S rRNA gene and whole genome shotgun metagenomic sequencing.

UTILITY AND DISCUSSION

All molecular data generated within the 1000IBD project will be shared on the European Genome-Phenome Archive ( https://ega-archive.org , accession no: EGAS00001002702). The first data release, detailed in this announcement and released simultaneously with this publication, will contain basic phenotypes for 1215 participants, genotypes of 314 participants and gut microbiome data from stool samples (315 participants) and biopsies (107 participants) generated by tag sequencing the 16S gene. Future releases will comprise many more additional phenotypes and -omics data layers. 1000IBD data can be used by other researchers as a replication cohort, a dataset to test new software tools, or a dataset for applying new statistical models.

CONCLUSIONS

We report on the establishment and future development of the 1000IBD project: the first comprehensive multi-omics dataset aimed at discovering IBD biomarker profiles and treatment targets.

Investigating multiple sclerosis genetic susceptibility on the founder population of east-central Sardinia via association and linkage analysis of immune-related loci

BACKGROUND

A wealth of single-nucleotide polymorphisms (SNPs) responsible for multiple sclerosis (MS) susceptibility have been identified; however, they explain only a fraction of MS heritability.

OBJECTIVES

We contributed to discovery of new MS susceptibility SNPs by studying a founder population with high MS prevalence.

METHODS

We analyzed ImmunoChip data from 15 multiplex families and 94 unrelated controls from the Nuoro Province, Sardinia, Italy. We tested each SNP for both association and linkage with MS, the linkage being explored in terms of identity-by-descent (IBD) sharing excess and using gene dropping to compute a corresponding empirical p-value. By targeting regions that are both associated and in linkage with MS, we increase chances of identifying interesting genomic regions.

RESULTS

We identified 486 MS-associated (p < 1 × 10-4) and 18,426 MS-linked (p < 0.05) SNPs. A total of 111 loci were both linked and associated with MS, 18 of them pointing to 14 non-major histocompatibility complex (MHC) genes, and 93 of them located in the MHC region.

CONCLUSION

We discovered new suggestive signals and confirmed some previously identified ones. We believe this to represent a significant step toward an understanding of the genetic basis of MS.

Detecting past and ongoing natural selection among ethnically Tibetan women at high altitude in Nepal

Adaptive evolution in humans has rarely been characterized for its whole set of components, i.e. selective pressure, adaptive phenotype, beneficial alleles and realized fitness differential. We combined approaches for detecting polygenic adaptations and for mapping the genetic bases of physiological and fertility phenotypes in approximately 1000 indigenous ethnically Tibetan women from Nepal, adapted to high altitude. The results of genome-wide association analyses and tests for polygenic adaptations showed evidence of positive selection for alleles associated with more pregnancies and live births and evidence of negative selection for those associated with higher offspring mortality. Lower hemoglobin level did not show clear evidence for polygenic adaptation, despite its strong association with an EPAS1 haplotype carrying selective sweep signals.

Investigation of common, low-frequency and rare genome-wide variation in anorexia nervosa

Anorexia nervosa (AN) is a complex neuropsychiatric disorder presenting with dangerously low body weight, and a deep and persistent fear of gaining weight. To date, only one genome-wide significant locus associated with AN has been identified. We performed an exome-chip based genome-wide association studies (GWAS) in 2158 cases from nine populations of European origin and 15 485 ancestrally matched controls. Unlike previous studies, this GWAS also probed association in low-frequency and rare variants. Sixteen independent variants were taken forward for in silico and de novo replication (11 common and 5 rare). No findings reached genome-wide significance. Two notable common variants were identified: rs10791286, an intronic variant in OPCML (P=9.89 × 10-6), and rs7700147, an intergenic variant (P=2.93 × 10-5). No low-frequency variant associations were identified at genome-wide significance, although the study was well-powered to detect low-frequency variants with large effect sizes, suggesting that there may be no AN loci in this genomic search space with large effect sizes.

Interplay of host genetics and gut microbiota underlying the onset and clinical presentation of inflammatory bowel disease

OBJECTIVE

Patients with IBD display substantial heterogeneity in clinical characteristics. We hypothesise that individual differences in the complex interaction of the host genome and the gut microbiota can explain the onset and the heterogeneous presentation of IBD. Therefore, we performed a case-control analysis of the gut microbiota, the host genome and the clinical phenotypes of IBD.

DESIGN

Stool samples, peripheral blood and extensive phenotype data were collected from 313 patients with IBD and 582 truly healthy controls, selected from a population cohort. The gut microbiota composition was assessed by tag-sequencing the 16S rRNA gene. All participants were genotyped. We composed genetic risk scores from 11 functional genetic variants proven to be associated with IBD in genes that are directly involved in the bacterial handling in the gut: NOD2, CARD9, ATG16L1, IRGM and FUT2.

RESULTS

Strikingly, we observed significant alterations of the gut microbiota of healthy individuals with a high genetic risk for IBD: the IBD genetic risk score was significantly associated with a decrease in the genus Roseburia in healthy controls (false discovery rate 0.017). Moreover, disease location was a major determinant of the gut microbiota: the gut microbiota of patients with colonic Crohn's disease (CD) is different from that of patients with ileal CD, with a decrease in alpha diversity associated to ileal disease (p=3.28×10^-13).

CONCLUSIONS

We show for the first time that genetic risk variants associated with IBD influence the gut microbiota in healthy individuals. Roseburia spp are acetate-to-butyrate converters, and a decrease has already been observed in patients with IBD.

The NOD2 receptor is crucial for immune responses towards New World Leishmania species

American Tegumentary Leishmaniasis is a chronic infection caused by Leishmania protozoan. It is not known whether genetic variances in NOD-like receptor (NLR) family members influence the immune response towards Leishmania parasites and modulate intracellular killing. Using functional genomics, we investigated whether genetic variants in NOD1 or NOD2 influence the production of cytokines by human PBMCs exposed to Leishmania. In addition, we examined whether recognition of Leishmania by NOD2 contributes to intracellular killing. Polymorphisms in the NOD2 gene decreased monocyte- and lymphocyte-derived cytokine production after stimulation with L. amazonensis or L. braziliensis compared to individuals with a functional NOD2 receptor. The phagolysosome formation is important for Leishmania-induced cytokine production and upregulation of NOD2 mRNA expression. NOD2 is crucial to control intracellular infection caused by Leishmania spp. NOD2 receptor is important for Leishmania recognition, the control of intracellular killing, and the induction of innate and adaptive immune responses.

The effects of signal transducer and activator of transcription three mutations on human platelets

Involvement of signal transducer and activator of transcription 3 (STAT3) in inflammation is well known. Recently, a role for STAT3 in platelet activation and platelet production has been suggested. Platelets exhibit important immune functions and engagement of STAT3 in platelet physiology may link inflammation and hemostasis. This study investigated the effects of STAT3 loss-of-function mutations and single nucleotide polymorphisms (SNPs) in STAT3 on glycoprotein VI (GPVI)-mediated platelet activation and platelet numbers in humans. Two cohorts were studied. The first cohort concerned patients with STAT3 loss-of-function mutations. Platelet numbers were investigated in eight patients and GPVI-mediated platelet activation was functionally tested in four patients. Additional experiments were performed to investigate underlying mechanisms. The second cohort concerned 334 healthy volunteers and investigated the consequences of SNPs in STAT3 on GPVI-mediated platelet activation and platelet numbers. Platelet activation was lower in STAT3 loss-of-function patients at baseline and after stimulation of the GPVI receptor, reflected by decreased P-selectin expression. This was independent of gene transcription. Blockade of the adenosine di-phosphate (ADP) pathway resulted in a further decrease of P-selectin expression, particularly in STAT3 loss-of-function patients. In contrast, the SNPs in STAT3 did not influence GPVI-mediated platelet activation. Also, platelet numbers were not affected by STAT3 loss-of-function mutations, nor was there an association with the SNPs. In conclusion, STAT3 signaling does not seem to play a major role in thrombopoiesis. We confirm that STAT3 is involved in GPVI-mediated platelet activation in humans, independent of gene transcription. GPVI-mediated platelet activation is highly dependent on secondary ADP release. Our findings suggest that STAT3 modulation may affect inflammation, hemostasis, and their interaction.

Transancestral mapping and genetic load in systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is an autoimmune disease with marked gender and ethnic disparities. We report a large transancestral association study of SLE using Immunochip genotype data from 27,574 individuals of European (EA), African (AA) and Hispanic Amerindian (HA) ancestry. We identify 58 distinct non-HLA regions in EA, 9 in AA and 16 in HA (∼50% of these regions have multiple independent associations); these include 24 novel SLE regions (P<5 × 10⁻⁸), refined association signals in established regions, extended associations to additional ancestries, and a disentangled complex HLA multigenic effect. The risk allele count (genetic load) exhibits an accelerating pattern of SLE risk, leading us to posit a cumulative hit hypothesis for autoimmune disease. Comparing results across the three ancestries identifies both ancestry-dependent and ancestry-independent contributions to SLE risk. Our results are consistent with the unique and complex histories of the populations sampled, and collectively help clarify the genetic architecture and ethnic disparities in SLE.

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong ethnic and gender bias. In a transancestral genetic association study, Langefeld et al. identify 24 novel regions associated with risk to lupus and propose a cumulative hits hypothesis for loci conferring risk to SLE.

Fine-mapping inflammatory bowel disease loci to single-variant resolution

Inflammatory bowel diseases are chronic gastrointestinal inflammatory disorders that affect millions of people worldwide. Genome-wide association studies have identified 200 inflammatory bowel disease-associated loci, but few have been conclusively resolved to specific functional variants. Here we report fine-mapping of 94 inflammatory bowel disease loci using high-density genotyping in 67,852 individuals. We pinpoint 18 associations to a single causal variant with greater than 95% certainty, and an additional 27 associations to a single variant with greater than 50% certainty. These 45 variants are significantly enriched for protein-coding changes (n = 13), direct disruption of transcription-factor binding sites (n = 3), and tissue-specific epigenetic marks (n = 10), with the last category showing enrichment in specific immune cells among associations stronger in Crohn's disease and in gut mucosa among associations stronger in ulcerative colitis. The results of this study suggest that high-resolution fine-mapping in large samples can convert many discoveries from genome-wide association studies into statistically convincing causal variants, providing a powerful substrate for experimental elucidation of disease mechanisms.

Dense genotyping of immune-related loci identifies HLA variants associated with increased risk of collagenous colitis

OBJECTIVE

Collagenous colitis (CC) is a major cause of chronic non-bloody diarrhoea, particularly in the elderly female population. The aetiology of CC is unknown, and still poor is the understanding of its pathogenesis. This possibly involves dysregulated inflammation and immune-mediated reactions in genetically predisposed individuals, but the contribution of genetic factors to CC is underinvestigated. We systematically tested immune-related genes known to impact the risk of several autoimmune diseases for their potential CC-predisposing role.

DESIGN

Three independent cohorts of histologically confirmed CC cases (N=314) and controls (N=4299) from Sweden and Germany were included in a 2-step association analysis. Immunochip and targeted single nucleotide polymorphism (SNP) genotype data were produced, respectively, for discovery and replication purposes. Classical human leucocyte antigen (HLA) variants at 2-digit and 4-digit resolution were obtained via imputation from single marker genotypes. SNPs and HLA variants passing quality control filters were tested for association with CC with logistic regression adjusting for age, sex and country of origin.

RESULTS

Forty-two markers gave rise to genome-wide significant association signals, all contained within the HLA region on chromosome 6 (best p=4.2×10^-10 for SNP rs4143332). Among the HLA variants, most pronounced risk effects were observed for 8.1 haplotype alleles including DQ2.5, which was targeted and confirmed in the replication data set (p=2.3×10^-11; OR=2.06; 95% CI (1.67 to 2.55) in the combined analysis).

CONCLUSIONS

HLA genotype associates with CC, thus implicating HLA-related immune mechanisms in its pathogenesis.

Association of IFIH1 and pro-inflammatory mediators: Potential new clues in SLE-associated pathogenesis

Antiviral defenses are inappropriately activated in systemic lupus erythematosus (SLE) and association between SLE and the antiviral helicase gene, IFIH1, is well established. We sought to extend the previously reported association of pathogenic soluble mediators and autoantibodies with mouse Mda5 to its human ortholog, IFIH1. To better understand the role this gene plays in human lupus, we assessed association of IFIH1 variants with soluble mediators and autoantibodies in 357 European-American SLE patients, first-degree relatives, and unrelated, unaffected healthy controls. Association between each of 135 genotyped SNPs in IFIH1 and four lupus-associated plasma mediators, IL-6, TNF-α, IFN-β, and IP-10, were investigated via linear regression. No significant associations were found to SNPs orthologous to those identified in exon 13 of the mouse. However, outside of this region there were significant associations between IL-6 and rs76162067 (p = 0.008), as well as IP-10 and rs79711023 (p = 0.003), located in a region of IFIH1 previously shown to directly influence MDA-5 mediated IP-10 and IL-6 secretion. SLE patients and FDRs carrying the minor allele for rs79711023 demonstrated lower levels of IP-10, while only FDRs carrying the minor allele for rs76162067 demonstrated an increased level of IL-6. This would suggest that the change in IP-10 is genotypically driven, while the change in IL-6 may be reflective of SLE transition status. These data suggest that IFIH1 may contribute to SLE pathogenesis via altered inflammatory mechanisms.

Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease

Genetic association studies have identified 215 risk loci for inflammatory bowel disease, thereby uncovering fundamental aspects of its molecular biology. We performed a genome-wide association study of 25,305 individuals and conducted a meta-analysis with published summary statistics, yielding a total sample size of 59,957 subjects. We identified 25 new susceptibility loci, 3 of which contain integrin genes that encode proteins in pathways that have been identified as important therapeutic targets in inflammatory bowel disease. The associated variants are correlated with expression changes in response to immune stimulus at two of these genes (ITGA4 and ITGB8) and at previously implicated loci (ITGAL and ICAM1). In all four cases, the expression-increasing allele also increases disease risk. We also identified likely causal missense variants in a gene implicated in primary immune deficiency, PLCG2, and a negative regulator of inflammation, SLAMF8. Our results demonstrate that new associations at common variants continue to identify genes relevant to therapeutic target identification and prioritization.

Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease

Genetic association studies have identified 215 risk loci for inflammatory bowel disease, thereby uncovering fundamental aspects of its molecular biology. We performed a genome-wide association study of 25,305 individuals and conducted a meta-analysis with published summary statistics, yielding a total sample size of 59,957 subjects. We identified 25 new susceptibility loci, 3 of which contain integrin genes that encode proteins in pathways that have been identified as important therapeutic targets in inflammatory bowel disease. The associated variants are correlated with expression changes in response to immune stimulus at two of these genes (ITGA4 and ITGB8) and at previously implicated loci (ITGAL and ICAM1). In all four cases, the expression-increasing allele also increases disease risk. We also identified likely causal missense variants in a gene implicated in primary immune deficiency, PLCG2, and a negative regulator of inflammation, SLAMF8. Our results demonstrate that new associations at common variants continue to identify genes relevant to therapeutic target identification and prioritization.

Functional and Genomic Architecture of Borrelia burgdorferi-Induced Cytokine Responses in Humans

Despite the importance of immune variation for the symptoms and outcome of Lyme disease, the factors influencing cytokine production during infection with the causal pathogen Borrelia burgdorferi remain poorly understood. Borrelia infection-induced monocyte- and T cell-derived cytokines were profiled in peripheral blood from two healthy human cohorts of Western Europeans from the Human Functional Genomics Project. Both non-genetic and genetic host factors were found to influence Borrelia-induced cytokine responses. Age strongly impaired IL-22 responses, and genetic studies identified several independent QTLs that impact Borrelia-induced cytokine production. Genetic, transcriptomic, and functional validation studies revealed an important role for HIF-1α-mediated glycolysis in the cytokine response to Borrelia. HIF-1α pathway activation and increase in glycolysis-derived lactate was confirmed in Lyme disease patients. In conclusion, functional genomics approaches reveal the architecture of cytokine production induced by Borrelia infection of human primary leukocytes and suggest a connection between cellular glucose metabolism and Borrelia-induced cytokine production.

The effect of host genetics on the gut microbiome

The gut microbiome is affected by multiple factors, including genetics. In this study, we assessed the influence of host genetics on microbial species, pathways and gene ontology categories, on the basis of metagenomic sequencing in 1,514 subjects. In a genome-wide analysis, we identified associations of 9 loci with microbial taxonomies and 33 loci with microbial pathways and gene ontology terms at P < 5 × 10^-8. Additionally, in a targeted analysis of regions involved in complex diseases, innate and adaptive immunity, or food preferences, 32 loci were identified at the suggestive level of P < 5 × 10^-6. Most of our reported associations are new, including genome-wide significance for the C-type lectin molecules CLEC4F-CD207 at 2p13.3 and CLEC4A-FAM90A1 at 12p13. We also identified association of a functional LCT SNP with the Bifidobacterium genus (P = 3.45 × 10^-8) and provide evidence of a gene-diet interaction in the regulation of Bifidobacterium abundance. Our results demonstrate the importance of understanding host-microbe interactions to gain better insight into human health.

Inter-individual variability and genetic influences on cytokine responses to bacteria and fungi

Little is known about the inter-individual variation of cytokine responses to different pathogens in healthy individuals. To systematically describe cytokine responses elicited by distinct pathogens and to determine the effect of genetic variation on cytokine production, we profiled cytokines produced by peripheral blood mononuclear cells from 197 individuals of European origin from the 200 Functional Genomics (200FG) cohort in the Human Functional Genomics Project (http://www.humanfunctionalgenomics.org), obtained over three different years. We compared bacteria- and fungi-induced cytokine profiles and found that most cytokine responses were organized around a physiological response to specific pathogens, rather than around a particular immune pathway or cytokine. We then correlated genome-wide single-nucleotide polymorphism (SNP) genotypes with cytokine abundance and identified six cytokine quantitative trait loci (QTLs). Among them, a cytokine QTL at the NAA35-GOLM1 locus markedly modulated interleukin (IL)-6 production in response to multiple pathogens and was associated with susceptibility to candidemia. Furthermore, the cytokine QTLs that we identified were enriched among SNPs previously associated with infectious diseases and heart diseases. These data reveal and begin to explain the variability in cytokine production by human immune cells in response to pathogens.

Combined role of vitamin D status and CYP24A1 in the transition to systemic lupus erythematosus

OBJECTIVE

We examined whether measures of vitamin D were associated with transitioning to systemic lupus erythematosus (SLE) in individuals at risk for SLE.

METHODS

436 individuals who reported having a relative with SLE but who did not have SLE themselves were evaluated at baseline and again an average of 6.3 (±3.9) years later. Fifty-six individuals transitioned to SLE (≥4 cumulative American College of Rheumatology criteria). 25-Hydroxyvitamin D (25[OH]D) levels were measured by ELISA. Six single-nucleotide polymorphisms in four vitamin D genes were genotyped. Generalised estimating equations, adjusting for correlation within families, were used to test associations between the vitamin D variables and the outcome of transitioning to SLE.

RESULTS

Mean baseline 25[OH]D levels (p=0.42) and vitamin D supplementation (p=0.65) were not different between those who did and did not transition to SLE. Vitamin D deficiency (25[OH]D <20 ng/mL) was greater in those who transitioned compared with those who did not transition to SLE (46% vs 33%, p=0.05). The association between 25[OH]D and SLE was modified by CYP24A1 rs4809959, where for each additional minor allele increased 25[OH]D was associated with decreased SLE risk: zero minor alleles (adjusted OR: 1.03, CI 0.98 to 1.09), one minor allele (adjusted OR: 1.01, CI 0.97 to 1.05) and two minor alleles (adjusted OR: 0.91, CI 0.84 to 0.98). Similarly, vitamin D deficiency significantly increased the risk of transitioning to SLE in those with two minor alleles at rs4809959 (adjusted OR: 4.90, CI 1.33 to 18.04).

CONCLUSIONS

Vitamin D status and CYP24A1 may have a combined role in the transition to SLE in individuals at increased genetic risk for SLE.

Genetic Evaluation of Schizophrenia Using the Illumina HumanExome Chip

INTRODUCTION

Schizophrenia is a genetically heterogeneous disorder that is associated with several common and rare genetic variants. As technology involved, cost advantages of chip based genotyping was combined with information about rare variants, resulting in the Infinium HumanExome Beadchip. Using this chip, a sample of 493 patients with schizophrenia or schizoaffective disorder and 484 healthy controls was genotyped.

RESULTS

From the initial 242901 SNVs, 88306 had at least one minor allele and passed quality control. No variant reached genomewide-significant results (p<10(-8)). The SNP with the lowest p-value was rs1230345 in WISP3 (p = 3.05*10(-6)), followed by rs9311525 in CACNA2D3 (p = 1.03*10(-5)) and rs1558557 (p = 3.85*10(-05)) on chromosome 7. At the gene level, 3 genes were of interest: WISP3, on chromosome 6q21, a signally protein from the extracellular matrix. A second candidate gene is CACNA2D3, a regulator of the intracerebral calcium pathway. A third gene is TNFSF10, associated with p53 mediated apoptosis.