Search

Summary Microbiome influences multiple human systems, but its effects on gene methylation is unknown. We investigated the relations between gene methylation in blood and the abundance of common gut bacteria profiled by 16s rRNA gene sequencing in two population-based Dutch cohorts: LifeLines-Deep (LLD, n = 616, discovery) and the Netherlands Twin Register (NTR, n = 296, replication). In LLD, we also explored microbial pathways using data generated by shotgun metagenomic sequencing (n = 683). Methylation in both cohorts was profiled in blood samples using the Illumina 450K array. Discovery and replication analysis identified two independent CpGs associated with the genus Eggerthella: cg16586104 (Pmeta−analysis = 3.21 × 10−11) and cg12234533 (Pmeta−analysis = 4.29 × 10−10). We also show that microbiome can mediate the effect of environmental factors on host gene methylation. In this first association study linking epigenome to microbiome, we found and replicated the associations of two CpGs to the abundance of genus Eggerthella and identified microbiome as a mediator of the exposome. These associations are observational and suggest further investigation in larger and longitudinal set-ups.

Abstract Genetic variants identified through genome-wide association studies (GWAS) are typically non-coding, exerting small regulatory effects on downstream genes. However, which downstream genes are ultimately impacted and how they confer risk remains mostly unclear. By contrast, variants that cause rare Mendelian diseases are often coding and have a more direct impact on disease development. Here we demonstrate that common and rare genetic diseases can be linked by studying how common disease-associated variants influence gene co-expression in 57 different tissues and cell types. We implemented this method in a framework called Downstreamer and applied it to 88 GWAS traits. We find that predicted downstream “genes” are enriched with Mendelian disease genes, e.g. key genes for height are enriched for genes that cause skeletal abnormalities and Ehlers–Danlos syndromes. 78% of these key genes are located outside of GWAS loci, suggesting that they result from complex trans regulation rather than being impacted by disease-associated variants in cis. Based on our findings, we discuss the challenges in reconstructing gene regulatory networks and provide a roadmap to improve the identification of these highly connected genes linked to common traits and diseases.

Abstract Background The plasma metabolome reflects the physiological state of various biological processes and can serve as a proxy for disease risk. Plasma metabolite variation, influenced by genetic and epigenetic mechanisms, can also affect the cellular microenvironment and blood cell epigenetics. The interplay between the plasma metabolome and the blood cell epigenome remains elusive. In this study, we performed an epigenome-wide association study (EWAS) of 1183 plasma metabolites in 693 participants from the LifeLines-DEEP cohort and investigated the causal relationships in DNA methylation–metabolite associations using bidirectional Mendelian randomization and mediation analysis. Results After rigorously adjusting for potential confounders, including genetics, we identified five robust associations between two plasma metabolites (l-serine and glycine) and three CpG sites located in two independent genomic regions (cg14476101 and cg16246545 in PHGDH and cg02711608 in SLC1A5) at a false discovery rate of less than 0.05. Further analysis revealed a complex bidirectional relationship between plasma glycine/serine levels and DNA methylation. Moreover, we observed a strong mediating role of DNA methylation in the effect of glycine/serine on the expression of their metabolism/transport genes, with the proportion of the mediated effect ranging from 11.8 to 54.3%. This result was also replicated in an independent population-based cohort, the Rotterdam Study. To validate our findings, we conducted in vitro cell studies which confirmed the mediating role of DNA methylation in the regulation of PHGDH gene expression. Conclusions Our findings reveal a potential feedback mechanism in which glycine and serine regulate gene expression through DNA methylation.

Abstract Background The genetic background of cancer remains complex and challenging to integrate. Many somatic mutations within genes are known to cause and drive cancer, while genome-wide association studies (GWAS) of cancer have revealed many germline risk factors associated with cancer. However, the overlap between known somatic driver genes and positional candidate genes from GWAS loci is surprisingly small. We hypothesised that genes from multiple independent cancer GWAS loci should show tissue-specific co-regulation patterns that converge on cancer-specific driver genes. Results We studied recent well-powered GWAS of breast, prostate, colorectal and skin cancer by estimating co-expression between genes and subsequently prioritising genes that show significant co-expression with genes mapping within susceptibility loci from cancer GWAS. We observed that the prioritised genes were strongly enriched for cancer drivers defined by COSMIC, IntOGen and Dietlein et al. The enrichment of known cancer driver genes was most significant when using co-expression networks derived from non-cancer samples of the relevant tissue of origin. Conclusion We show how genes within risk loci identified by cancer GWAS can be linked to known cancer driver genes through tissue-specific co-expression networks. This provides an important explanation for why seemingly unrelated sets of genes that harbour either germline risk factors or somatic mutations can eventually cause the same type of disease.

Abstract Recent innovations in single-cell RNA-sequencing (scRNA-seq) provide the technology to investigate biological questions at cellular resolution. Pooling cells from multiple individuals has become a common strategy, and droplets can subsequently be assigned to a specific individual by leveraging their inherent genetic differences. An implicit challenge with scRNA-seq is the occurrence of doublets—droplets containing two or more cells. We develop Demuxafy, a framework to enhance donor assignment and doublet removal through the consensus intersection of multiple demultiplexing and doublet detecting methods. Demuxafy significantly improves droplet assignment by separating singlets from doublets and classifying the correct individual.

Summary: Background: 1H-NMR metabolomics and DNA methylation in blood are widely known biomarkers predicting age-related physiological decline and mortality yet exert mutually independent mortality and frailty signals. Methods: Leveraging multi-omics data in four Dutch population studies (N = 5238, ∼40% of which male) we investigated whether the mortality signal captured by 1H-NMR metabolomics could guide the construction of DNA methylation-based mortality predictors. Findings: We trained DNA methylation-based surrogates for 64 metabolomic analytes and found that analytes marking inflammation, fluid balance, or HDL/VLDL metabolism could be accurately reconstructed using DNA-methylation assays. Interestingly, a previously reported multi-analyte score indicating mortality risk (MetaboHealth) could also be accurately reconstructed. Sixteen of our derived surrogates, including the MetaboHealth surrogate, showed significant associations with mortality, independent of relevant covariates. Interpretation: The addition of our metabolic analyte-derived surrogates to the well-established epigenetic clock GrimAge demonstrates that our surrogates potentially represent valuable mortality signal. Funding: BBMRI-NL, X-omics, VOILA, Medical Delta, NWO, ERC.

Summary: Human induced pluripotent stem cell (hiPSC)-derived intestinal organoids are valuable tools for researching developmental biology and personalized therapies, but their closed topology and relative immature state limit applications. Here, we use organ-on-chip technology to develop a hiPSC-derived intestinal barrier with apical and basolateral access in a more physiological in vitro microenvironment. To replicate growth factor gradients along the crypt-villus axis, we locally expose the cells to expansion and differentiation media. In these conditions, intestinal epithelial cells self-organize into villus-like folds with physiological barrier integrity, and myofibroblasts and neurons emerge and form a subepithelial tissue in the bottom channel. The growth factor gradients efficiently balance dividing and mature cell types and induce an intestinal epithelial composition, including absorptive and secretory lineages, resembling the composition of the human small intestine. This well-characterized hiPSC-derived intestine-on-chip system can facilitate personalized studies on physiological processes and therapy development in the human small intestine.

Abstract Expression quantitative trait loci (eQTL) offer insights into the regulatory mechanisms of trait-associated variants, but their effects often rely on contexts that are unknown or unmeasured. We introduce PICALO, a method for hidden variable inference of eQTL contexts. PICALO identifies and disentangles technical from biological context in heterogeneous blood and brain bulk eQTL datasets. These contexts are biologically informative and reproducible, outperforming cell counts or expression-based principal components. Furthermore, we show that RNA quality and cell type proportions interact with thousands of eQTLs. Knowledge of hidden eQTL contexts may aid in the inference of functional mechanisms underlying disease variants.

Summary: Encounters with pathogens and other molecules can imprint long-lasting effects on our immune system, influencing future physiological outcomes. Given the wide range of microbes to which humans are exposed, their collective impact on health is not fully understood. To explore relations between exposures and biological aging and inflammation, we profiled an antibody-binding repertoire against 2,815 microbial, viral, and environmental peptides in a population cohort of 1,443 participants. Utilizing antibody-binding as a proxy for past exposures, we investigated their impact on biological aging, cell composition, and inflammation. Immune response against cytomegalovirus (CMV), rhinovirus, and gut bacteria relates with telomere length. Single-cell expression measurements identified an effect of CMV infection on the transcriptional landscape of subpopulations of CD8 and CD4 T-cells. This examination of the relationship between microbial exposures and biological aging and inflammation highlights a role for chronic infections (CMV and Epstein-Barr virus) and common pathogens (rhinoviruses and adenovirus C).

IntroductionIn the evolving healthcare landscape, precision medicine's rise necessitates adaptable doctoral training. The European Union has recognized this and promotes the development of international, training-focused programmes called Innovative Training Networks (ITNs). In this article, we introduce TranSYS, an ITN focused on educating the next generation of precision medicine researchers. In an ambition to go beyond describing the consortium goals, our article explores two key aspects of ITNs: the training and collaboration.MethodsUsing self-report questionnaires, we evaluate the scientific, professional, and personal growth of ESRs over the duration of the ITN and investigate whether this can be linked to network activities.ResultsOur quantitative analysis approach reveals substantial improvements in scientific, professional, and social skills among young researchers facilitated by the engagement in this interdisciplinary network. We provide case studies underlining the advantages of collaborative environments, featuring innovative scientific exchange within TranSYS.DiscussionWhile challenging, ITNs foster positive growth in young researchers, yet exhibit weaknesses such as balancing stakeholder interests and partner commitment. We believe this study may benefit a variety of stakeholders, from prospective ITN creators to industry partners, to design better sustainable training networks going forward.

Summary: Background: Smoking impacts DNA methylation, but data are lacking on smoking-related differential methylation by sex or dietary intake, recent smoking cessation (

Abstract Background Expression quantitative trait loci (eQTL) studies show how genetic variants affect downstream gene expression. Single-cell data allows reconstruction of personalized co-expression networks and therefore the identification of SNPs altering co-expression patterns (co-expression QTLs, co-eQTLs) and the affected upstream regulatory processes using a limited number of individuals. Results We conduct a co-eQTL meta-analysis across four scRNA-seq peripheral blood mononuclear cell datasets using a novel filtering strategy followed by a permutation-based multiple testing approach. Before the analysis, we evaluate the co-expression patterns required for co-eQTL identification using different external resources. We identify a robust set of cell-type-specific co-eQTLs for 72 independent SNPs affecting 946 gene pairs. These co-eQTLs are replicated in a large bulk cohort and provide novel insights into how disease-associated variants alter regulatory networks. One co-eQTL SNP, rs1131017, that is associated with several autoimmune diseases, affects the co-expression of RPS26 with other ribosomal genes. Interestingly, specifically in T cells, the SNP additionally affects co-expression of RPS26 and a group of genes associated with T cell activation and autoimmune disease. Among these genes, we identify enrichment for targets of five T-cell-activation-related transcription factors whose binding sites harbor rs1131017. This reveals a previously overlooked process and pinpoints potential regulators that could explain the association of rs1131017 with autoimmune diseases. Conclusion Our co-eQTL results highlight the importance of studying context-specific gene regulation to understand the biological implications of genetic variation. With the expected growth of sc-eQTL datasets, our strategy and technical guidelines will facilitate future co-eQTL identification, further elucidating unknown disease mechanisms.

Here, the authors present a TWAS framework OTTERS that adapts multiple polygenic risk score methods to estimate eQTL weights from summary-level eQTL data. Both simulation and real studies show OTTERS is powerful across a wide range of genetic architectures.

Circulating lipids can influence immune cell function, which could have implications for inflammatory diseases such as rheumatoid arthritis. Here, the authors use Mendelian randomization to identify genes whose expression is influenced by triglyceride levels in blood, implicating genes involved in lipid metabolism and allergic response.

BackgroundDuring the COVID-19 pandemic, many healthcare workers faced extreme working conditions and were at higher risk of infection with the coronavirus. These circumstances may have led to mental health problems, such as anxiety, among healthcare workers. Most studies that examined anxiety among healthcare workers during the COVID-19 pandemic were cross-sectional and focused on the first months of the pandemic only. Therefore, this study aimed to investigate the longitudinal association between working in healthcare and anxiety during a long-term period (i.e., 18 months) of the COVID-19 pandemic.MethodsData were used from online questionnaires of the Lifelines COVID-19 prospective cohort with 22 included time-points (March 2020–November 2021). In total, 2,750 healthcare workers and 9,335 non-healthcare workers were included. Anxiety was assessed with questions from the Mini-International Neuropsychiatric Interview, and an anxiety sum score (0–7) was calculated. Negative binomial generalized estimating equations (GEE), adjusted for demographic, work and health covariates, were used to examine the association between working in healthcare and anxiety.ResultsAnxiety sum scores over time during the COVID-19 pandemic were similar for healthcare workers and non-healthcare workers. No differences between the anxiety sum scores of healthcare workers and non-healthcare workers were found [incidence rate ratio (IRR) = 0.97, 95% CI = 0.91–1.04].ConclusionThis study did not find differences between healthcare workers and non-healthcare in perceived anxiety during the COVID-19 pandemic.

Abstract Background The genetic underpinning of sexual dimorphism is very poorly understood. The prevalence of many diseases differs between men and women, which could be in part caused by sex-specific genetic effects. Nevertheless, only a few published genome-wide association studies (GWAS) were performed separately in each sex. The reported enrichment of expression quantitative trait loci (eQTLs) among GWAS-associated SNPs suggests a potential role of sex-specific eQTLs in the sex-specific genetic mechanism underlying complex traits. Methods To explore this scenario, we combined sex-specific whole blood RNA-seq eQTL data from 3447 European individuals included in BIOS Consortium and GWAS data from UK Biobank. Next, to test the presence of sex-biased causal effect of gene expression on complex traits, we performed sex-specific transcriptome-wide Mendelian randomization (TWMR) analyses on the two most sexually dimorphic traits, waist-to-hip ratio (WHR) and testosterone levels. Finally, we performed power analysis to calculate the GWAS sample size needed to observe sex-specific trait associations driven by sex-biased eQTLs. Results Among 9 million SNP-gene pairs showing sex-combined associations, we found 18 genes with significant sex-biased cis-eQTLs (FDR 5%). Our phenome-wide association study of the 18 top sex-biased eQTLs on >700 traits unraveled that these eQTLs do not systematically translate into detectable sex-biased trait-associations. In addition, we observed that sex-specific causal effects of gene expression on complex traits are not driven by sex-specific eQTLs. Power analyses using real eQTL- and causal-effect sizes showed that millions of samples would be necessary to observe sex-biased trait associations that are fully driven by sex-biased cis-eQTLs. Compensatory effects may further hamper their detection. Conclusions Our results suggest that sex-specific eQTLs in whole blood do not translate to detectable sex-specific trait associations of complex diseases, and vice versa that the observed sex-specific trait associations cannot be explained by sex-specific eQTLs.

Abstract Single-cell RNA sequencing (scRNA-seq) has revolutionized the study of the cellular landscape of organs. Most single-cell protocols require fresh material, which limits sample size per experiment, and consequently, introduces batch effects. This is especially true for samples acquired through complex medical procedures, such as intestinal mucosal biopsies. Moreover, the tissue dissociation procedure required for obtaining single cells is a major source of noise; different dissociation procedures applied to different compartments of the tissue induce artificial gene expression differences between cell subsets. To overcome these challenges, we have developed a one-step dissociation protocol and demonstrated its use on cryopreserved gut mucosal biopsies. Using flow cytometry and scRNA-seq analysis, we compared this one-step dissociation protocol with the current gold standard, two-step collagenase digestion, and an adaptation of a recently published alternative, three-step cold-active Bacillus licheniformus protease digestion. Both cell viability and cell type composition were comparable between the one-step and two-step collagenase dissociation, with the former being more time-efficient. The cold protease digestion resulted in equal cell viability, but better preserves the epithelial cell types. Consequently, to analyze the rarer cell types, such as glial cells, larger total biopsy cell numbers are required as input material. The multi-step protocols affected cell types spanning multiple compartments differently. In summary, we show that cryopreserved gut mucosal biopsies can be used to overcome the logistical challenges and batch effects in large scRNA-seq studies. Furthermore, we demonstrate that using cryopreserved biopsies digested using a one-step collagenase protocol enables large-scale scRNA-seq, FACS, organoid generation and intraepithelial lymphocyte expansion.

Deep molecular and phenotypic data highlights the links of human telomere lengths from six different blood cells with genetics, parental phenotypes mediated by epigenetic signals and expression changes at the single cell level.

Not just differential gene expression but also differential gene regulation in immune cells account for individual differences in the immune response. Authors show here by single-cell RNA-sequencing of peripheral blood mononuclear cells from a large cohort of genetically diverse individuals that gene expression and regulatory changes in these cells depend on the context of and interactions between cell types, genetics, type of pathogen and time after exposure.

Summary: Genetics have nominated many schizophrenia risk genes and identified convergent signals between schizophrenia and neurodevelopmental disorders. However, functional interpretation of the nominated genes in the relevant brain cell types is often lacking. We executed interaction proteomics for six schizophrenia risk genes that have also been implicated in neurodevelopment in human induced cortical neurons. The resulting protein network is enriched for common variant risk of schizophrenia in Europeans and East Asians, is down-regulated in layer 5/6 cortical neurons of individuals affected by schizophrenia, and can complement fine-mapping and eQTL data to prioritize additional genes in GWAS loci. A sub-network centered on HCN1 is enriched for common variant risk and contains proteins (HCN4 and AKAP11) enriched for rare protein-truncating mutations in individuals with schizophrenia and bipolar disorder. Our findings showcase brain cell-type-specific interactomes as an organizing framework to facilitate interpretation of genetic and transcriptomic data in schizophrenia and its related disorders.

Abstract Background Epigenetic clocks use DNA methylation (DNAm) levels of specific sets of CpG dinucleotides to accurately predict individual chronological age. A popular application of these clocks is to explore whether the deviation of predicted age from chronological age is associated with disease phenotypes, where this deviation is interpreted as a potential biomarker of biological age. This wide application, however, contrasts with the limited insight in the processes that may drive the running of epigenetic clocks. Results We perform a functional genomics analysis on four epigenetic clocks, including Hannum’s blood predictor and Horvath’s multi-tissue predictor, using blood DNA methylome and transcriptome data from 3132 individuals. The four clocks result in similar predictions of individual chronological age, and their constituting CpGs are correlated in DNAm level and are enriched for similar histone modifications and chromatin states. Interestingly, DNAm levels of CpGs from the clocks are commonly associated with gene expression in trans. The gene sets involved are highly overlapping and enriched for T cell processes. Further analysis of the transcriptome and methylome of sorted blood cell types identifies differences in DNAm between naive and activated T and NK cells as a probable contributor to the clocks. Indeed, within the same donor, the four epigenetic clocks predict naive cells to be up to 40 years younger than activated cells. Conclusions The ability of epigenetic clocks to predict chronological age involves their ability to detect changes in proportions of naive and activated immune blood cells, an established feature of immuno-senescence. This finding may contribute to the interpretation of associations between clock-derived measures and age-related health outcomes.

Both gene expression and protein concentrations are regulated by genetic variants. Exploring the regulation of both eQTLs and pQTLs simultaneously in a context- and cell-type dependent manner may help to unravel mechanistic basis for genetic regulation of pQTLs. Here, we performed meta-analysis of Candida albicans-induced pQTLs from two population-based cohorts and intersected the results with Candida-induced cell-type specific expression association data (eQTL). This revealed systematic differences between the pQTLs and eQTL, where only 35% of the pQTLs significantly correlated with mRNA expressions at single cell level, indicating the limitation of eQTLs use as a proxy for pQTLs. By taking advantage of the tightly co-regulated pattern of the proteins, we also identified SNPs affecting protein network upon Candida stimulations. Colocalization of pQTLs and eQTLs signals implicated several genomic loci including MMP-1 and AMZ1. Analysis of Candida-induced single cell gene expression data implicated specific cell types that exhibit significant expression QTLs upon stimulation. By highlighting the role of trans-regulatory networks in determining the abundance of secretory proteins, our study serve as a framework to gain insights into the mechanisms of genetic regulation of protein levels in a context-dependent manner.

Summary: Polygenic risk scores (PRSs) have been widely explored in precision medicine. However, few studies have thoroughly investigated their best practices in global populations across different diseases. We here utilized data from Global Biobank Meta-analysis Initiative (GBMI) to explore methodological considerations and PRS performance in 9 different biobanks for 14 disease endpoints. Specifically, we constructed PRSs using pruning and thresholding (P + T) and PRS-continuous shrinkage (CS). For both methods, using a European-based linkage disequilibrium (LD) reference panel resulted in comparable or higher prediction accuracy compared with several other non-European-based panels. PRS-CS overall outperformed the classic P + T method, especially for endpoints with higher SNP-based heritability. Notably, prediction accuracy is heterogeneous across endpoints, biobanks, and ancestries, especially for asthma, which has known variation in disease prevalence across populations. Overall, we provide lessons for PRS construction, evaluation, and interpretation using GBMI resources and highlight the importance of best practices for PRS in the biobank-scale genomics era.

Many genetic loci have been identified to be associated with kidney disease, but the molecular mechanisms are not well understood. Here, the authors perform epigenome-wide association studies on kidney function measures to identify epigenetic marks and pathways involved in kidney function.

IntroductionWorking from home during the COVID-19 pandemic has been associated both with physical inactivity and musculoskeletal pain. However, it has not been examined whether physical activity and sedentary behavior are underlying mechanisms in the association between working from home and musculoskeletal pain. Therefore, we examined their mediating role in this association.MethodsData were used from 24 questionnaire rounds of the Lifelines COVID-19 cohort (March 2020–January 2022). Longitudinal information on work situation (location, home, hybrid), physical activity, sedentary behavior, and musculoskeletal pain was collected among 28,586 workers. Analysis of physical activity/sedentary behavior as mediators of the association between working from home and musculoskeletal pain was performed using multilevel structural equation modeling.ResultsHome workers more often had pain in the upper back [odds ratio (OR) = 1.17, 95%-confidence interval (CI) = 1.02–1.34] and arm, neck, and/or shoulder (ANS) (OR = 1.32, 95%-CI = 1.19–1.47) than location workers. Furthermore, home workers were more often sedentary for >9 h per work day than location workers (OR = 2.82, 95%-CI = 2.56–3.09), and being more sedentary was associated with musculoskeletal pain (upper back: OR = 1.17, 95%-CI = 1.06–1.30; ANS: OR = 1.25, 95%-CI = 1.16–1.34). Corresponding indirect effects were OR = 1.18 (95%-CI = 1.04–1.33) and OR = 1.26 (95%-CI = 1.12–1.35). No indirect effect was found for physical activity. Similar indirect effects were observed for hybrid workers.ConclusionHome and hybrid workers were more likely to have pain in the upper musculoskeletal system during the COVID-19 pandemic than location workers, which was partly mediated by increased sedentary behavior, but not by reduced physical activity. Measures to reduce sedentary time in home workers may contribute to preventing musculoskeletal pain.

Summary: Copy-number variations (CNV) are believed to play an important role in a wide range of complex traits, but discovering such associations remains challenging. While whole-genome sequencing (WGS) is the gold-standard approach for CNV detection, there are several orders of magnitude more samples with available genotyping microarray data. Such array data can be exploited for CNV detection using dedicated software (e.g., PennCNV); however, these calls suffer from elevated false-positive and -negative rates. In this study, we developed a CNV quality score that weights PennCNV calls (pCNVs) based on their likelihood of being true positive. First, we established a measure of pCNV reliability by leveraging evidence from multiple omics data (WGS, transcriptomics, and methylomics) obtained from the same samples. Next, we built a predictor of omics-confirmed pCNVs, termed omics-informed quality score (OQS), using only PennCNV software output parameters. Promisingly, OQS assigned to pCNVs detected in close family members was up to 35% higher than the OQS of pCNVs not carried by other relatives (p

Abstract 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.

While coffee and tea consumption has been associated with risk of diseases, their mechanisms of action remain elusive. Here the authors present a large EWAS on coffee and tea consumption in cohorts of European and African-American ancestries, finding that coffee consumption is associated with differential DNA methylation levels at multiple CpGs.

Abstract Background Aging is a multifactorial process that affects multiple tissues and is characterized by changes in homeostasis over time, leading to increased morbidity. Whole blood gene expression signatures have been associated with aging and have been used to gain information on its biological mechanisms, which are still not fully understood. However, blood is composed of many cell types whose proportions in blood vary with age. As a result, previously observed associations between gene expression levels and aging might be driven by cell type composition rather than intracellular aging mechanisms. To overcome this, previous aging studies already accounted for major cell types, but the possibility that the reported associations are false positives driven by less prevalent cell subtypes remains. Results Here, we compared the regression model from our previous work to an extended model that corrects for 33 additional white blood cell subtypes. Both models were applied to whole blood gene expression data from 3165 individuals belonging to the general population (age range of 18–81 years). We evaluated that the new model is a better fit for the data and it identified fewer genes associated with aging (625, compared to the 2808 of the initial model; P ≤ 2.5⨯10−6). Moreover, 511 genes (~ 18% of the 2808 genes identified by the initial model) were found using both models, indicating that the other previously reported genes could be proxies for less abundant cell types. In particular, functional enrichment of the genes identified by the new model highlighted pathways and GO terms specifically associated with platelet activity. Conclusions We conclude that gene expression analyses in blood strongly benefit from correction for both common and rare blood cell types, and recommend using blood-cell count estimates as standard covariates when studying whole blood gene expression.

Abstract Background Expression quantitative trait loci (eQTL) studies are used to interpret the function of disease-associated genetic risk factors. To date, most eQTL analyses have been conducted in bulk tissues, such as whole blood and tissue biopsies, which are likely to mask the cell type-context of the eQTL regulatory effects. Although this context can be investigated by generating transcriptional profiles from purified cell subpopulations, current methods to do this are labor-intensive and expensive. We introduce a new method, Decon2, as a framework for estimating cell proportions using expression profiles from bulk blood samples (Decon-cell) followed by deconvolution of cell type eQTLs (Decon-eQTL). Results The estimated cell proportions from Decon-cell agree with experimental measurements across cohorts (R ≥ 0.77). Using Decon-cell, we could predict the proportions of 34 circulating cell types for 3194 samples from a population-based cohort. Next, we identified 16,362 whole-blood eQTLs and deconvoluted cell type interaction (CTi) eQTLs using the predicted cell proportions from Decon-cell. CTi eQTLs show excellent allelic directional concordance with eQTL (≥ 96–100%) and chromatin mark QTL (≥87–92%) studies that used either purified cell subpopulations or single-cell RNA-seq, outperforming the conventional interaction effect. Conclusions Decon2 provides a method to detect cell type interaction effects from bulk blood eQTLs that is useful for pinpointing the most relevant cell type for a given complex disease. Decon2 is available as an R package and Java application ( https://github.com/molgenis/systemsgenetics/tree/master/Decon2 ) and as a web tool ( www.molgenis.org/deconvolution ).

Here, the authors provide an in-depth study of the metabolome in two large population-based prospective cohorts and identify 32 microbial traits associated with various metabolic biomarkers and specific lipoprotein subfractions, providing insights into the role of microbiota in influencing host lipid levels.

Abstract Background A large number of analysis strategies are available for DNA methylation (DNAm) array and RNA-seq datasets, but it is unclear which strategies are best to use. We compare commonly used strategies and report how they influence results in large cohort studies. Results We tested the associations of DNAm and RNA expression with age, BMI, and smoking in four different cohorts (n = ~ 2900). By comparing strategies against the base model on the number and percentage of replicated CpGs for DNAm analyses or genes for RNA-seq analyses in a leave-one-out cohort replication approach, we find the choice of the normalization method and statistical test does not strongly influence the results for DNAm array data. However, adjusting for cell counts or hidden confounders substantially decreases the number of replicated CpGs for age and increases the number of replicated CpGs for BMI and smoking. For RNA-seq data, the choice of the normalization method, gene expression inclusion threshold, and statistical test does not strongly influence the results. Including five principal components or excluding correction of technical covariates or cell counts decreases the number of replicated genes. Conclusions Results were not influenced by the normalization method or statistical test. However, the correction method for cell counts, technical covariates, principal components, and/or hidden confounders does influence the results.

Inbreeding depression has been observed in many different species, but in humans a systematic analysis has been difficult so far. Here, analysing more than 1.3 million individuals, the authors show that a genomic inbreeding coefficient (FROH) is associated with disadvantageous outcomes in 32 out of 100 traits tested.

Abstract Background Identification of imprinted genes, demonstrating a consistent preference towards the paternal or maternal allelic expression, is important for the understanding of gene expression regulation during embryonic development and of the molecular basis of developmental disorders with a parent-of-origin effect. Combining allelic analysis of RNA-Seq data with phased genotypes in family trios provides a powerful method to detect parent-of-origin biases in gene expression. Results We report findings in 296 family trios from two large studies: 165 lymphoblastoid cell lines from the 1000 Genomes Project and 131 blood samples from the Genome of the Netherlands (GoNL) participants. Based on parental haplotypes, we identified > 2.8 million transcribed heterozygous SNVs phased for parental origin and developed a robust statistical framework for measuring allelic expression. We identified a total of 45 imprinted genes and one imprinted unannotated transcript, including multiple imprinted transcripts showing incomplete parental expression bias that was located adjacent to strongly imprinted genes. For example, PXDC1, a gene which lies adjacent to the paternally expressed gene FAM50B, shows a 2:1 paternal expression bias. Other imprinted genes had promoter regions that coincide with sites of parentally biased DNA methylation identified in the blood from uniparental disomy (UPD) samples, thus providing independent validation of our results. Using the stranded nature of the RNA-Seq data in lymphoblastoid cell lines, we identified multiple loci with overlapping sense/antisense transcripts, of which one is expressed paternally and the other maternally. Using a sliding window approach, we searched for imprinted expression across the entire genome, identifying a novel imprinted putative lncRNA in 13q21.2. Overall, we identified 7 transcripts showing parental bias in gene expression which were not reported in 4 other recent RNA-Seq studies of imprinting. Conclusions Our methods and data provide a robust and high-resolution map of imprinted gene expression in the human genome.

A genetic diagnosis remains unattainable for many individuals with a rare disease because of incomplete knowledge about the genetic basis of many diseases. Here, the authors present the web-based tool GADO (GeneNetwork Assisted Diagnostic Optimization) that uses public RNA-seq data for prioritization of candidate genes.

Recent contributions highlighted gender differences in the mental health consequences of COVID-19 lockdowns. However, their cross-sectional designs cannot differentiate between pre-existing gender differences and differences induced by lockdowns. Estimating fixed-effects models using longitudinal data from the Lifelines biobank and cohort study with repeated mental health measurements throughout the lockdown, we overcome this caveat. Significant gender differences in mental health during the lockdown were found, where women experienced more depression symptoms and disorders and men experienced more anxiety symptoms and disorders stemming from the lockdown. Policymakers need to keep in mind that the COVID-19 lockdowns have different effects on mental health for men and women.

Circulatory inflammatory proteins play a significant role in anti-Candida host immune defence. However, little is known about the genetic variation that contributes to the variability of inflammatory responses in response to C. albicans. To systematically characterize inflammatory responses in Candida infection, we profiled 91 circulatory inflammatory proteins in peripheral blood mononuclear cells (PBMCs) stimulated with C. albicans yeast isolated from 378 individuals of European origin from the 500 Functional Genomics (500FG) cohort of the Human Functional Genomics Project (HFGP) and Lifelines Deep cohort. To identify the genetic factors that determine variation in inflammatory protein responses, we correlated genome-wide single nucleotide polymorphism (SNP) genotypes with protein abundance (protein quantitative trait loci, pQTLs) produced by the Candida-stimulated PBMCs. Furthermore, we investigated whether differences in survival of candidaemia patients can be explained by modulating levels of inflammatory proteins. We identified five genome-wide significant pQTLs that modulate IL-8, MCP-2, MMP-1, and CCL3 in response to C. albicans. In addition, our genetic analysis suggested that GADD45G from rs10114707 locus that reached genome-wide significance could be a potential core gene that regulates a cytokine network upon Candida infection. Last but not least, we observed that a trans-pQTL marked from SNP rs7651677 at chromosome 3 that influences urokinase plasminogen activator (uPA) is strongly associated with patient survival (Psurvival = 3.52 x 10-5, OR 3). Overall, our genetic analysis showed that genetic variation determines the abundance of circulatory proteins in response to Candida infection.

Educational attainment is a key behavioural measure in studies of cognitive and physical health, and socioeconomic status. We measured DNA methylation at 410,746 CpGs (N = 4152) and identified 58 CpGs associated with educational attainment at loci characterized by pleiotropic functions shared with neuronal, immune and developmental processes. Associations overlapped with those for smoking behaviour, but remained after accounting for smoking at many CpGs: Effect sizes were on average 28% smaller and genome-wide significant at 11 CpGs after adjusting for smoking and were 62% smaller in never smokers. We examined sources and biological implications of education-related methylation differences, demonstrating correlations with maternal prenatal folate, smoking and air pollution signatures, and associations with gene expression in cis, dynamic methylation in foetal brain, and correlations between blood and brain. Our findings show that the methylome of lower-educated people resembles that of smokers beyond effects of their own smoking behaviour and shows traces of various other exposures.

Purpose The Lifelines COVID-19 cohort was set up to assess the psychological and societal impacts of the COVID-19 pandemic and investigate potential risk factors for COVID-19 within the Lifelines prospective population cohort.Participants Participants were recruited from the 140 000 eligible participants of Lifelines and the Lifelines NEXT birth cohort, who are all residents of the three northern provinces of the Netherlands. Participants filled out detailed questionnaires about their physical and mental health and experiences on a weekly basis starting in late March 2020, and the cohort consists of everyone who filled in at least one questionnaire in the first 8 weeks of the project.Findings to date >71 000 unique participants responded to the questionnaires at least once during the first 8 weeks, with >22 000 participants responding to seven questionnaires. Compiled questionnaire results are continuously updated and shared with the public through the Corona Barometer website. Early results included a clear signal that younger people living alone were experiencing greater levels of loneliness due to lockdown, and subsequent results showed the easing of anxiety as lockdown was eased in June 2020.Future plans Questionnaires were sent on a (bi)weekly basis starting in March 2020 and on a monthly basis starting July 2020, with plans for new questionnaire rounds to continue through 2020 and early 2021. Questionnaire frequency can be increased again for subsequent waves of infections. Cohort data will be used to address how the COVID-19 pandemic developed in the northern provinces of the Netherlands, which environmental and genetic risk factors predict disease susceptibility and severity and the psychological and societal impacts of the crisis. Cohort data are linked to the extensive health, lifestyle and sociodemographic data held for these participants by Lifelines, a 30-year project that started in 2006, and to data about participants held in national databases.

Objectives Studies in clinical settings showed a potential relationship between socioeconomic status (SES) and lifestyle factors with COVID-19, but it is still unknown whether this holds in the general population. In this study, we investigated the associations of SES with self-reported, tested and diagnosed COVID-19 status in the general population.Design, setting, participants and outcome measures Participants were 49 474 men and women (46±12 years) residing in the Northern Netherlands from the Lifelines cohort study. SES indicators and lifestyle factors (i.e., smoking status, physical activity, alcohol intake, diet quality, sleep time and TV watching time) were assessed by questionnaire from the Lifelines Biobank. Self-reported, tested and diagnosed COVID-19 status was obtained from the Lifelines COVID-19 questionnaire.Results There were 4711 participants who self-reported having had a COVID-19 infection, 2883 participants tested for COVID-19, and 123 positive cases were diagnosed in this study population. After adjustment for age, sex, lifestyle factors, body mass index and ethnicity, we found that participants with low education or low income were less likely to self-report a COVID-19 infection (OR [95% CI]: low education 0.78 [0.71 to 0.86]; low income 0.86 [0.79 to 0.93]) and be tested for COVID-19 (OR [95% CI]: low education 0.58 [0.52 to 0.66]; low income 0.86 [0.78 to 0.95]) compared with high education or high income groups, respectively.Conclusion Our findings suggest that the low SES group was the most vulnerable population to self-reported and tested COVID-19 status in the general population.

Celiac disease (CeD) is a complex T cell-mediated enteropathy induced by gluten. Although genome-wide association studies have identified numerous genomic regions associated with CeD, it is difficult to accurately pinpoint which genes in these loci are most likely to cause CeD. We used four different in silico approaches—Mendelian randomization inverse variance weighting, COLOC, LD overlap, and DEPICT—to integrate information gathered from a large transcriptomics dataset. This identified 118 prioritized genes across 50 CeD-associated regions. Co-expression and pathway analysis of these genes indicated an association with adaptive and innate cytokine signaling and T cell activation pathways. Fifty-one of these genes are targets of known drug compounds or likely druggable genes, suggesting that our methods can be used to pinpoint potential therapeutic targets. In addition, we detected 172 gene combinations that were affected by our CeD-prioritized genes in trans. Notably, 41 of these trans-mediated genes appear to be under control of one master regulator, TRAF-type zinc finger domain containing 1 (TRAFD1), and were found to be involved in interferon (IFN)γ signaling and MHC I antigen processing/presentation. Finally, we performed in vitro experiments in a human monocytic cell line that validated the role of TRAFD1 as an immune regulator acting in trans. Our strategy confirmed the role of adaptive immunity in CeD and revealed a genetic link between CeD and IFNγ signaling as well as with MHC I antigen processing, both major players of immune activation and CeD pathogenesis.

Background: Exposure to environmental endocrine disrupting chemicals (EDCs) may play an important role in the epidemic of metabolic diseases. Epigenetic alterations may functionally link EDCs with gene expression and metabolic traits. Objectives: We aimed to evaluate metabolic-related effects of the exposure to endocrine disruptors including five parabens, three bisphenols, and 13 metabolites of nine phthalates as measured in 24-hour urine on epigenome-wide DNA methylation. Methods: A blood-based epigenome-wide association study was performed in 622 participants from the Lifelines DEEP cohort using Illumina Infinium HumanMethylation450 methylation data and EDC excretions in 24-hour urine. Out of the 21 EDCs, 13 compounds were detected in >75% of the samples and, together with bisphenol F, were included in these analyses. Furthermore, we explored the putative function of identified methylation markers and their correlations with metabolic traits. Results: We found 20 differentially methylated cytosine-phosphate-guanines (CpGs) associated with 10 EDCs at suggestive p-value

Coronavirus disease 2019 (COVID-19) shows a wide variation in expression and severity of symptoms, from very mild or no symptoms, to flu-like symptoms, and in more severe cases, to pneumonia, acute respiratory distress syndrome, and even death. Large differences in outcome have also been observed between males and females. The causes for this variability are likely to be multifactorial, and to include genetics. The SARS-CoV-2 virus responsible for the infection depends on two human genes: the human receptor angiotensin converting enzyme 2 (ACE2) for cell invasion, and the serine protease TMPRSS2 for S protein priming. Genetic variation in these two genes may thus modulate an individual's genetic predisposition to infection and virus clearance. While genetic data on COVID-19 patients is being gathered, we carried out a phenome-wide association scan (PheWAS) to investigate the role of these genes in other human phenotypes in the general population. We examined 178 quantitative phenotypes including cytokines and cardio-metabolic biomarkers, as well as usage of 58 medications in 36,339 volunteers from the Lifelines population cohort, in relation to 1,273 genetic variants located in or near ACE2 and TMPRSS2. While none reached our threshold for significance, we observed several interesting suggestive associations. For example, single nucleotide polymorphisms (SNPs) near the TMPRSS2 genes were associated with thrombocytes count (p = 1.8 × 10−5). SNPs within the ACE2 gene were associated with (1) the use of angiotensin II receptor blockers (ARBs) combination therapies (p = 5.7 × 10−4), an association that is significantly stronger in females (pdiff = 0.01), and (2) with the use of non-steroid anti-inflammatory and antirheumatic products (p = 5.5 × 10−4). While these associations need to be confirmed in larger sample sizes, they suggest that these variants could play a role in diseases such as thrombocytopenia, hypertension, and chronic inflammation that are often observed in the more severe COVID-19 cases. Further investigation of these genetic variants in the context of COVID-19 is thus promising for better understanding of disease variability. Full results are available at https://covid19research.nl.

Candida bloodstream infection, i.e. candidemia, is the most frequently encountered life-threatening fungal infection worldwide, with mortality rates up to almost 50%. In the majority of candidemia cases, Candida albicans is responsible. Worryingly, a global increase in the number of patients who are susceptible to infection (e.g. immunocompromised patients), has led to a rise in the incidence of candidemia in the last few decades. Therefore, a better understanding of the anti-Candida host response is essential to overcome this poor prognosis and to lower disease incidence. Here, we integrated genome-wide association studies with bulk and single-cell transcriptomic analyses of immune cells stimulated with Candida albicans to further our understanding of the anti-Candida host response. We show that differential expression analysis upon Candida stimulation in single-cell expression data can reveal the important cell types involved in the host response against Candida. This confirmed the known major role of monocytes, but more interestingly, also uncovered an important role for NK cells. Moreover, combining the power of bulk RNA-seq with the high resolution of single-cell RNA-seq data led to the identification of 27 Candida-response QTLs and revealed the cell types potentially involved herein. Integration of these response QTLs with a GWAS on candidemia susceptibility uncovered a potential new role for LY86 in candidemia susceptibility. Finally, experimental follow-up confirmed that LY86 knockdown results in reduced monocyte migration towards the chemokine MCP-1, thereby implying that this reduced migration may underlie the increased susceptibility to candidemia. Altogether, our integrative systems genetics approach identifies previously unknown mechanisms underlying the immune response to Candida infection.

Abstract Only a small fraction of patients respond to the drug prescribed to treat their disease, which means that most are at risk of unnecessary exposure to side effects through ineffective drugs. This inter-individual variation in drug response is driven by differences in gene interactions caused by each patient’s genetic background, environmental exposures, and the proportions of specific cell types involved in disease. These gene interactions can now be captured by building gene regulatory networks, by taking advantage of RNA velocity (the time derivative of the gene expression state), the ability to study hundreds of thousands of cells simultaneously, and the falling price of single-cell sequencing. Here, we propose an integrative approach that leverages these recent advances in single-cell data with the sensitivity of bulk data to enable the reconstruction of personalized, cell-type- and context-specific gene regulatory networks. We expect this approach will allow the prioritization of key driver genes for specific diseases and will provide knowledge that opens new avenues towards improved personalized healthcare.