Blood circulating proteins are confounded readouts of the biological processes that occur in different tissues and organs. Many proteins have been linked to complex disorders and are also under ...substantial genetic control. Here, we investigate the associations between over 1000 blood circulating proteins and body mass index (BMI) in three studies including over 4600 participants. We show that BMI is associated with widespread changes in the plasma proteome. We observe 152 replicated protein associations with BMI. 24 proteins also associate with a genome-wide polygenic score (GPS) for BMI. These proteins are involved in lipid metabolism and inflammatory pathways impacting clinically relevant pathways of adiposity. Mendelian randomization suggests a bi-directional causal relationship of BMI with LEPR/LEP, IGFBP1, and WFIKKN2, a protein-to-BMI relationship for AGER, DPT, and CTSA, and a BMI-to-protein relationship for another 21 proteins. Combined with animal model and tissue-specific gene expression data, our findings suggest potential therapeutic targets further elucidating the role of these proteins in obesity associated pathologies.
Sample collection, processing, storage and isolation methods constitute pre-analytic factors that can influence the quality of samples used in research and clinical practice. With regard to ...biobanking practices, a critical point in the sample's life chain is storage, particularly long-term storage. Since most studies examine the influence of different temperatures (4°C, room temperature) or delays in sample processing on sample quality, there is only little information on the effects of long-term storage at ultra-low (vapor phase of liquid nitrogen) temperatures on biomarker levels. Among these biomarkers, circulating miRNAs hold great potential for diagnosis or prognosis for a variety of diseases, like cancer, infections and chronic diseases, and are thus of high interest in several scientific questions. We therefore investigated the influence of long-term storage on levels of eight circulating miRNAs (miR-103a-3p, miR-191-5p, miR-124-3p, miR-30c-5p, miR-451a, miR-23a-3p, miR-93-5p, miR-24-3p, and miR-33b-5p) from 10 participants from the population-based cohort study KORA. Sample collection took place during the baseline survey S4 and the follow-up surveys F4 and FF4, over a time period spanning from 1999 to 2014. The influence of freeze-thaw (f/t) cycles on miRNA stability was also investigated using samples from volunteers (n = 6). Obtained plasma samples were profiled using Exiqon's miRCURYTM real-time PCR profiling system, and repeated measures ANOVA was used to check for storage or f/t effects. Our results show that detected levels of most of the studied miRNAs showed no statistically significant changes due to storage at ultra-low temperatures for up to 17 years; miR-451a levels were altered due to contamination during sampling. Freeze-thawing of one to four cycles showed an effect only on miR-30c-5p. Our results highlight the robustness of this set of circulating miRNAs for decades of storage at ultra-low temperatures and several freeze-thaw cycles, which makes our findings increasingly relevant for research conducted with biobanked samples.
Despite extensive research on 5-methylcytosine (5mC) in relation to smoking, there has been limited exploration into the interaction between smoking and 5-hydroxymethylcytosine (5hmC). In this study, ...total DNA methylation (5mC+5hmC), true DNA methylation (5mC) and hydroxymethylation (5hmC) levels were profiled utilizing conventional bisulphite (BS) and oxidative bisulphite (oxBS) treatment, measured with the Illumina Infinium Methylation EPIC BeadChip. An epigenome-wide association study (EWAS) of 5mC+5hmC methylation revealed a total of 38,575 differentially methylated positions (DMPs) and 2023 differentially methylated regions (DMRs) associated with current smoking, along with 82 DMPs and 76 DMRs associated with former smoking (FDR-adjusted p < 0.05). Additionally, a focused examination of 5mC identified 33 DMPs linked to current smoking and 1 DMP associated with former smoking (FDR-adjusted p < 0.05). In the 5hmC category, eight DMPs related to current smoking and two DMPs tied to former smoking were identified, each meeting a suggestive threshold (p < 1 × 10−5). The substantial number of recognized DMPs, including 5mC+5hmC (7069/38,575, 2/82), 5mC (0/33, 1/1), and 5hmC (2/8, 0/2), have not been previously reported. Our findings corroborated previously established methylation positions and revealed novel candidates linked to tobacco smoking. Moreover, the identification of hydroxymethylated CpG sites with suggestive links provides avenues for future research.
Studies on the relationship between renal function and the human plasma proteome have identified several potential biomarkers. However, investigations have been conducted largely in European ...populations, and causality of the associations between plasma proteins and kidney function has never been addressed.
A cross-sectional study of 993 plasma proteins among 2882 participants in four studies of European and admixed ancestries (KORA, INTERVAL, HUNT, QMDiab) identified transethnic associations between eGFR/CKD and proteomic biomarkers. For the replicated associations, two-sample bidirectional Mendelian randomization (MR) was used to investigate potential causal relationships. Publicly available datasets and transcriptomic data from independent studies were used to examine the association between gene expression in kidney tissue and eGFR.
In total, 57 plasma proteins were associated with eGFR, including one novel protein. Of these, 23 were additionally associated with CKD. The strongest inferred causal effect was the positive effect of eGFR on testican-2, in line with the known biological role of this protein and the expression of its protein-coding gene (
) in renal tissue. We also observed suggestive evidence of an effect of melanoma inhibitory activity (MIA), carbonic anhydrase III, and cystatin-M on eGFR.
In a discovery-replication setting, we identified 57 proteins transethnically associated with eGFR. The revealed causal relationships are an important stepping stone in establishing testican-2 as a clinically relevant physiological marker of kidney disease progression, and point to additional proteins warranting further investigation.
Dietary intake of n-3 polyunsaturated fatty acids (PUFA) may have a protective effect on the development of cardiovascular diseases, diabetes, depression and cancer, while a high intake of n-6 PUFA ...was often reported to be associated with inflammation-related traits. The effect of PUFAs on health outcomes might be mediated by DNA methylation (DNAm). The aim of our study is to identify the impact of PUFA intake on DNAm in the Cooperative Health Research in the Region of Augsburg (KORA) FF4 cohort and the Leiden Longevity Study (LLS).
DNA methylation levels were measured in whole blood from the population-based KORA FF4 study (N = 1354) and LLS (N = 448), using the Illumina MethylationEPIC BeadChip and Illumina HumanMethylation450 array, respectively. We assessed associations between DNAm and intake of eight and four PUFAs in KORA and LLS, respectively. Where possible, results were meta-analyzed. Below the Bonferroni correction threshold (p < 7.17 × 10
), we identified two differentially methylated positions (DMPs) associated with PUFA intake in the KORA study. The DMP cg19937480, annotated to gene PRDX1, was positively associated with docosahexaenoic acid (DHA) in model 1 (beta: 2.00 × 10
, 95%CI: 1.28 × 10
-2.73 × 10
, P value: 6.98 × 10
), while cg05041783, annotated to gene MARK2, was positively associated with docosapentaenoic acid (DPA) in our fully adjusted model (beta: 9.80 × 10
, 95%CI: 6.25 × 10
-1.33 × 10
, P value: 6.75 × 10
). In the meta-analysis, we identified the CpG site (cg15951061), annotated to gene CDCA7L below Bonferroni correction (1.23 × 10
) associated with eicosapentaenoic acid (EPA) intake in model 1 (beta: 2.00 × 10
, 95% CI: 1.27 × 10
-2.73 × 10
, P value = 5.99 × 10
) and we confirmed the association of cg19937480 with DHA in both models 1 and 2 (beta: 2.07 × 10
, 95% CI: 1.31 × 10
-2.83 × 10
, P value = 1.00 × 10
and beta: 2.19 × 10
, 95% CI: 1.41 × 10
-2.97 × 10
, P value = 5.91 × 10
respectively).
Our study identified three CpG sites associated with PUFA intake. The mechanisms of these sites remain largely unexplored, highlighting the novelty of our findings. Further research is essential to understand the links between CpG site methylation and PUFA outcomes.
We determined the relationships between DNA sequence variation and DNA methylation using blood samples from 3,799 Europeans and 3,195 South Asians. We identify 11,165,559 SNP-CpG associations ...(methylation quantitative trait loci (meQTL), P < 10
), including 467,915 meQTL that operate in trans. The meQTL are enriched for functionally relevant characteristics, including shared chromatin state, High-throuhgput chromosome conformation interaction, and association with gene expression, metabolic variation and clinical traits. We use molecular interaction and colocalization analyses to identify multiple nuclear regulatory pathways linking meQTL loci to phenotypic variation, including UBASH3B (body mass index), NFKBIE (rheumatoid arthritis), MGA (blood pressure) and COMMD7 (white cell counts). For rs6511961 , chromatin immunoprecipitation followed by sequencing (ChIP-seq) validates zinc finger protein (ZNF)333 as the likely trans acting effector protein. Finally, we used interaction analyses to identify population- and lineage-specific meQTL, including rs174548 in FADS1, with the strongest effect in CD8
T cells, thus linking fatty acid metabolism with immune dysregulation and asthma. Our study advances understanding of the potential pathways linking genetic variation to human phenotype.
Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression ...prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.
The difference between an individual's chronological and DNA methylation predicted age (DNAmAge), termed DNAmAge acceleration (DNAmAA), can capture life-long environmental exposures and age-related ...physiological changes reflected in methylation status. Several studies have linked DNAmAA to morbidity and mortality, yet its relationship with kidney function has not been assessed. We evaluated the associations between seven DNAm aging and lifespan predictors (as well as GrimAge components) and five kidney traits (estimated glomerular filtration rate eGFR, urine albumin-to-creatinine ratio uACR, serum urate, microalbuminuria and chronic kidney disease CKD) in up to 9688 European, African American and Hispanic/Latino individuals from seven population-based studies.
We identified 23 significant associations in our large trans-ethnic meta-analysis (p < 1.43E-03 and consistent direction of effect across studies). Age acceleration measured by the Extrinsic and PhenoAge estimators, as well as Zhang's 10-CpG epigenetic mortality risk score (MRS), were associated with all parameters of poor kidney health (lower eGFR, prevalent CKD, higher uACR, microalbuminuria and higher serum urate). Six of these associations were independently observed in European and African American populations. MRS in particular was consistently associated with eGFR (β = - 0.12, 95% CI = - 0.16, - 0.08 change in log-transformed eGFR per unit increase in MRS, p = 4.39E-08), prevalent CKD (odds ratio (OR) = 1.78 1.47, 2.16, p = 2.71E-09) and higher serum urate levels (β = 0.12 0.07, 0.16, p = 2.08E-06). The "first-generation" clocks (Hannum, Horvath) and GrimAge showed different patterns of association with the kidney traits. Three of the DNAm-estimated components of GrimAge, namely adrenomedullin, plasminogen-activation inhibition 1 and pack years, were positively associated with higher uACR, serum urate and microalbuminuria.
DNAmAge acceleration and DNAm mortality predictors estimated in whole blood were associated with multiple kidney traits, including eGFR and CKD, in this multi-ethnic study. Epigenetic biomarkers which reflect the systemic effects of age-related mechanisms such as immunosenescence, inflammaging and oxidative stress may have important mechanistic or prognostic roles in kidney disease. Our study highlights new findings linking kidney disease to biological aging, and opportunities warranting future investigation into DNA methylation biomarkers for prognostic or risk stratification in kidney disease.
Abstract
Educational inequalities in all-cause mortality have been observed for decades. However, the underlying biological mechanisms are not well known. We aimed to assess the role of DNA ...methylation changes in blood captured by epigenetic clocks in explaining these inequalities. Data were from 8 prospective population-based cohort studies, representing 13 021 participants. First, educational inequalities and their portion explained by Horvath DNAmAge, Hannum DNAmAge, DNAmPhenoAge, and DNAmGrimAge epigenetic clocks were assessed in each cohort via counterfactual-based mediation models, on both absolute (hazard difference) and relative (hazard ratio) scales, and by sex. Second, estimates from each cohort were pooled through a random effect meta-analysis model. Men with low education had excess mortality from all causes of 57 deaths per 10 000 person-years (95% confidence interval CI: 38, 76) compared with their more advantaged counterparts. For women, the excess mortality was 4 deaths per 10 000 person-years (95% CI: −11, 19). On the relative scale, educational inequalities corresponded to hazard ratios of 1.33 (95% CI: 1.12, 1.57) for men and 1.15 (95% CI: 0.96, 1.37) for women. DNAmGrimAge accounted for the largest proportion, approximately 50%, of the educational inequalities for men, while the proportion was negligible for women. Most of this mediation was explained by differential effects of unhealthy lifestyles and morbidities of the World Health Organization (WHO) risk factors for premature mortality. These results support DNA methylation-based epigenetic aging as a signature of educational inequalities in life expectancy emphasizing the need for policies to address the unequal social distribution of these WHO risk factors.
The discovery of robust and trans-ethnically replicated DNA methylation markers of metabolic phenotypes, has hinted at a potential role of epigenetic mechanisms in lipid metabolism. However, DNA ...methylation and the lipid compositions and lipid concentrations of lipoprotein sizes have been scarcely studied. Here, we present an epigenome-wide association study (EWAS) (N = 5414 total) of mostly lipid-related metabolic measures, including a fine profiling of lipoproteins. As lipoproteins are the main players in the different stages of lipid metabolism, examination of epigenetic markers of detailed lipoprotein features might improve the diagnosis, prognosis, and treatment of metabolic disturbances.
We conducted an EWAS of leukocyte DNA methylation and 226 metabolic measurements determined by nuclear magnetic resonance spectroscopy in the population-based KORA F4 study (N = 1662) and replicated the results in the LOLIPOP, NFBC1966, and YFS cohorts (N = 3752). Follow-up analyses in the discovery cohort included investigations into gene transcripts, metabolic-measure ratios for pathway analysis, and disease endpoints. We identified 161 associations (p value < 4.7 × 10
), covering 16 CpG sites at 11 loci and 57 metabolic measures. Identified metabolic measures were primarily medium and small lipoproteins, and fatty acids. For apolipoprotein B-containing lipoproteins, the associations mainly involved triglyceride composition and concentrations of cholesterol esters, triglycerides, free cholesterol, and phospholipids. All associations for HDL lipoproteins involved triglyceride measures only. Associated metabolic measure ratios, proxies of enzymatic activity, highlight amino acid, glucose, and lipid pathways as being potentially epigenetically implicated. Five CpG sites in four genes were associated with differential expression of transcripts in blood or adipose tissue. CpG sites in ABCG1 and PHGDH showed associations with metabolic measures, gene transcription, and metabolic measure ratios and were additionally linked to obesity or previous myocardial infarction, extending previously reported observations.
Our study provides evidence of a link between DNA methylation and the lipid compositions and lipid concentrations of different lipoprotein size subclasses, thus offering in-depth insights into well-known associations of DNA methylation with total serum lipids. The results support detailed profiling of lipid metabolism to improve the molecular understanding of dyslipidemia and related disease mechanisms.