Linking the gut microbiome to host DNA methylation by a discovery and replication epigenome-wide association study.

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

NMR metabolomics-guided DNA methylation mortality predictors.

Background: H-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 H-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.

Epigenome-wide association study on the plasma metabolome suggests self-regulation of the glycine and serine pathway through DNA methylation.

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.

An iPSC-derived small intestine-on-chip with self-organizing epithelial, mesenchymal, and neural cells.

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.

Antibody signatures against viruses and microbiome reflect past and chronic exposures and associate with aging and inflammation.

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.

Multiomic analyses uncover immunological signatures in acute and chronic coronary syndromes.

Acute and chronic coronary syndromes (ACS and CCS) are leading causes of mortality. Inflammation is considered a key pathogenic driver of these diseases, but the underlying immune states and their clinical implications remain poorly understood. Multiomic factor analysis (MOFA) allows unsupervised data exploration across multiple data types, identifying major axes of variation and associating these with underlying molecular processes.

Genetic Correlates of Biological Aging and the Influence on Prediction of Mortality.

Longevity and disease-free survival are influenced by a combination of genetics and lifestyle. Biological age (BioAge), a measure of aging based on composite biomarkers, may outperform chronological age in predicting health and longevity. This study investigated the relationship between genetic risks, lifestyle factors, and delta age (Δage), estimated as the difference between biological and chronological age.

PICALO: principal interaction component analysis for the identification of discrete technical, cell-type, and environmental factors that mediate eQTLs.

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.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively.

Semi-automatic translation of medicine usage data (in Dutch, free-text) from Lifelines COVID-19 questionnaires to ATC codes.

The mapping of human-entered data to codified data formats that can be analysed is a common problem across medical research and health care. To identify risk and protective factors for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) susceptibility and coronavirus disease 2019 (COVID-19) severity, frequent questionnaires were sent out to participants of the Lifelines Cohort Study starting 30 March 2020. Because specific drugs were suspected COVID-19 risk factors, the questionnaires contained multiple-choice questions about commonly used drugs and open-ended questions to capture all other drugs used.

Identification of genetic variants that impact gene co-expression relationships using large-scale single-cell data.

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.

OTTERS: a powerful TWAS framework leveraging summary-level reference data.

Most existing TWAS tools require individual-level eQTL reference data and thus are not applicable to summary-level reference eQTL datasets. The development of TWAS methods that can harness summary-level reference data is valuable to enable TWAS in broader settings and enhance power due to increased reference sample size. Thus, we develop a TWAS framework called OTTERS (Omnibus Transcriptome Test using Expression Reference Summary data) that adapts multiple polygenic risk score (PRS) methods to estimate eQTL weights from summary-level eQTL reference data and conducts an omnibus TWAS.

Integration of -induced single-cell gene expression data and secretory protein concentrations reveal genetic regulators of inflammation.

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 -induced pQTLs from two population-based cohorts and intersected the results with -induced cell-type specific expression association data (eQTL).

KidneyNetwork: using kidney-derived gene expression data to predict and prioritize novel genes involved in kidney disease.

Genetic testing in patients with suspected hereditary kidney disease may not reveal the genetic cause for the disorder as potentially pathogenic variants can reside in genes that are not yet known to be involved in kidney disease. We have developed KidneyNetwork, that utilizes tissue-specific expression to inform candidate gene prioritization specifically for kidney diseases. KidneyNetwork is a novel method constructed by integrating a kidney RNA-sequencing co-expression network of 878 samples with a multi-tissue network of 31,499 samples.

Lipid-induced transcriptomic changes in blood link to lipid metabolism and allergic response.

Immune cell function can be altered by lipids in circulation, a process potentially relevant to lipid-associated inflammatory diseases including atherosclerosis and rheumatoid arthritis. To gain further insight in the molecular changes involved, we here perform a transcriptome-wide association analysis of blood triglycerides, HDL cholesterol, and LDL cholesterol in 3229 individuals, followed by a systematic bidirectional Mendelian randomization analysis to assess the direction of effects and control for pleiotropy. Triglycerides are found to induce transcriptional changes in 55 genes and HDL cholesterol in 5 genes.