Epigenetic loci for blood pressure are associated with hypertensive target organ damage in older African Americans from the genetic epidemiology network of Arteriopathy (GENOA) study

Non-coding 886 (nc886/vtRNA2-1), the epigenetic odd duck - implications for future studies

Non-coding 886 (nc886, vtRNA2-1) is the only human polymorphically imprinted gene, in which the methylation status is not determined by genetics. Existing literature regarding the establishment, stability and consequences of the methylation pattern, as well as the nature and function of the nc886 RNAs transcribed from the locus, are contradictory. For example, the methylation status of the locus has been reported to be stable through life and across somatic tissues, but also susceptible to environmental effects. The nature of the produced nc886 RNA(s) has been redefined multiple times, and in carcinogenesis, these RNAs have been reported to have conflicting roles. In addition, due to the bimodal methylation pattern of the nc886 locus, traditional genome-wide methylation analyses can lead to false-positive results, especially in smaller datasets. Herein, we aim to summarize the existing literature regarding nc886, discuss how the characteristics of nc886 give rise to contradictory results, as well as to reinterpret, reanalyse and, where possible, replicate the results presented in the current literature. We also introduce novel findings on how the distribution of the nc886 methylation pattern is associated with the geographical origins of the population and describe the methylation changes in a large variety of human tumours. Through the example of this one peculiar genetic locus and RNA, we aim to highlight issues in the analysis of DNA methylation and non-coding RNAs in general and offer our suggestions for what should be taken into consideration in future analyses.

Development of an epigenetic clock resistant to changes in immune cell composition

Epigenetic clocks are age predictors that use machine-learning models trained on DNA CpG methylation values to predict chronological or biological age. Increases in predicted epigenetic age relative to chronological age (epigenetic age acceleration) are connected to aging-associated pathologies, and changes in epigenetic age are linked to canonical aging hallmarks. However, epigenetic clocks rely on training data from bulk tissues whose cellular composition changes with age. Here, we found that human naive CD8+ T cells, which decrease in frequency during aging, exhibit an epigenetic age 15-20 years younger than effector memory CD8+ T cells from the same individual. Importantly, homogenous naive T cells isolated from individuals of different ages show a progressive increase in epigenetic age, indicating that current epigenetic clocks measure two independent variables, aging and immune cell composition. To isolate the age-associated cell intrinsic changes, we created an epigenetic clock, the IntrinClock, that did not change among 10 immune cell types tested. IntrinClock shows a robust predicted epigenetic age increase in a model of replicative senescence in vitro and age reversal during OSKM-mediated reprogramming.

Map of epigenetic age acceleration: A worldwide analysis

We present a systematic analysis of epigenetic age acceleration based on by far the largest collection of publicly available DNA methylation data for healthy samples (93 datasets, 23 K samples), focusing on the geographic (25 countries) and ethnic (31 ethnicities) aspects around the world. We employed the most popular epigenetic tools for assessing age acceleration and examined their quality metrics and ability to extrapolate to epigenetic data from different tissue types and age ranges different from the training data of these models. In most cases, the models proved to be inconsistent with each other and showed different signs of age acceleration, with the PhenoAge model tending to systematically underestimate and different versions of the GrimAge model tending to systematically overestimate the age prediction of healthy subjects. Referring to data availability and consistency, most countries and populations are still not represented in GEO, moreover, different datasets use different criteria for determining healthy controls. Because of this, it is difficult to fully isolate the contribution of "geography/environment", "ethnicity" and "healthiness" to epigenetic age acceleration. Among the explored metrics, only the DunedinPACE, which measures aging rate, appears to adequately reflect the standard of living and socioeconomic indicators in countries, although it has a limited application to blood methylation data only. Invariably, by epigenetic age acceleration, males age faster than females in most of the studied countries and populations.

Accurate age prediction from blood using a small set of DNA methylation sites and a cohort-based machine learning algorithm

Chronological age prediction from DNA methylation sheds light on human aging, health, and lifespan. Current clocks are mostly based on linear models and rely upon hundreds of sites across the genome. Here, we present GP-age, an epigenetic non-linear cohort-based clock for blood, based upon 11,910 methylomes. Using 30 CpG sites alone, GP-age outperforms state-of-the-art models, with a median accuracy of ∼2 years on held-out blood samples, for both array and sequencing-based data. We show that aging-related changes occur at multiple neighboring CpGs, with implications for using fragment-level analysis of sequencing data in aging research. By training three independent clocks, we show enrichment of donors with consistent deviation between predicted and actual age, suggesting individual rates of biological aging. Overall, we provide a compact yet accurate alternative to array-based clocks for blood, with applications in longitudinal aging research, forensic profiling, and monitoring epigenetic processes in transplantation medicine and cancer.

Integrated epigenome, whole genome sequence and metabolome analyses identify novel multi-omics pathways in type 2 diabetes: a Middle Eastern study

BACKGROUND

T2D is of high prevalence in the middle east and thus studying its mechanisms is of a significant importance. Using 1026 Qatar BioBank samples, epigenetics, whole genome sequencing and metabolomics were combined to further elucidate the biological mechanisms of T2D in a population with a high prevalence of T2D.

METHODS

An epigenome-wide association study (EWAS) with T2D was performed using the Infinium 850K EPIC array, followed by whole genome-wide sequencing SNP-CpG association analysis (> 5.5 million SNPs) and a methylome-metabolome (CpG-metabolite) analysis of the identified T2D sites.

RESULTS

A total of 66 T2D-CpG associations were identified, including 63 novel sites in pathways of fructose and mannose metabolism, insulin signaling, galactose, starch and sucrose metabolism, and carbohydrate absorption and digestion. Whole genome SNP associations with the 66 CpGs resulted in 688 significant CpG-SNP associations comprising 22 unique CpGs (33% of the 66 CPGs) and included 181 novel pairs or pairs in novel loci. Fourteen of the loci overlapped published GWAS loci for diabetes related traits and were used to identify causal associations of HK1 and PFKFB2 with HbA1c. Methylome-metabolome analysis identified 66 significant CpG-metabolite pairs among which 61 pairs were novel. Using the identified methylome-metabolome associations, methylation QTLs, and metabolic networks, a multi-omics network was constructed which suggested a number of metabolic mechanisms underlying T2D methylated genes. 1-palmitoyl-2-oleoyl-GPE (16:0/18:1) - a triglyceride-associated metabolite, shared a common network with 13 methylated CpGs, including TXNIP, PFKFB2, OCIAD1, and BLCAP. Mannonate - a food component/plant shared a common network with 6 methylated genes, including TXNIP, BLCAP, THBS4 and PEF1, pointing to a common possible cause of methylation in those genes. A subnetwork with alanine, glutamine, urea cycle (citrulline, arginine), and 1-carboxyethylvaline linked to PFKFB2 and TXNIP revealed associations with kidney function, hypertension and triglyceride metabolism. The pathway containing STYXL1-POR was associated with a sphingosine-ceramides subnetwork associated with HDL-C and LDL-C and point to steroid perturbations in T2D.

CONCLUSIONS

This study revealed several novel methylated genes in T2D, with their genomic variants and associated metabolic pathways with several implications for future clinical use of multi-omics associations in disease and for studying therapeutic targets.

Association Between DNA Methylation and Blood Pressure: A 5-Year Longitudinal Twin Study

BACKGROUND

Previous EWASs (Epigenome-Wide Association Studies) have reported hundreds of blood pressure (BP) associated 5'-cytosine-phosphate-guanine-3' (CpG) sites. However, their results were inconsistent. Longitudinal observations on the temporal relationship between DNA methylation and BP are lacking.

METHODS

A candidate CpG site association study for BP was conducted on 1072 twins in the Chinese National Twin Registry. PubMed and EMBASE were searched for candidate CpG sites. Cross-lagged models were used to assess the temporal relationship between BP and DNA methylation in 308 twins who completed 2 surveys in 2013 and 2018. Then, the significant cross-lagged associations were validated by adopting the Inference About Causation From Examination of Familial Confounding approach. Finally, to evaluate the cumulative effects of DNA methylation on the progression of hypertension, we established methylation risk scores based on BP-associated CpG sites and performed Markov multistate models.

RESULTS

16 and 20 CpG sites were validated to be associated with systolic BP and diastolic BP, respectively. In the cross-lagged analysis, we detected that methylation of 2 CpG sites could predict subsequent systolic BP, and systolic BP predicted methylation at another 3 CpG sites. For diastolic BP, methylation at 3 CpG sites had significant cross-lagged effects for predicting diastolic BP levels, while the prediction from the opposite direction was observed at one site. Among these, 3 associations were validated in the Inference About Causation From Examination of Familial Confounding analysis. Using the Markov multistate model, we observed that methylation risk scores were associated with the development of hypertension.

CONCLUSIONS

Our findings suggest the significance of DNA methylation in the development of hypertension.

The Interplay of Epigenetic, Genetic, and Traditional Risk Factors on Blood Pressure: Findings from the Health and Retirement Study

The epigenome likely interacts with traditional and genetic risk factors to influence blood pressure. We evaluated whether 13 previously reported DNA methylation sites (CpGs) are associated with systolic (SBP) or diastolic (DBP) blood pressure, both individually and aggregated into methylation risk scores (MRS), in 3070 participants (including 437 African ancestry (AA) and 2021 European ancestry (EA), mean age = 70.5 years) from the Health and Retirement Study. Nine CpGs were at least nominally associated with SBP and/or DBP after adjusting for traditional hypertension risk factors (p < 0.05). MRSSBP was positively associated with SBP in the full sample (β = 1.7 mmHg per 1 standard deviation in MRSSBP; p = 2.7 × 10-5) and in EA (β = 1.6; p = 0.001), and MRSDBP with DBP in the full sample (β = 1.1; p = 1.8 × 10-6), EA (β = 1.1; p = 7.2 × 10-5), and AA (β = 1.4; p = 0.03). The MRS and BP-genetic risk scores were independently associated with blood pressure in EA. The effects of both MRSs were weaker with increased age (pinteraction < 0.01), and the effect of MRSDBP was higher among individuals with at least some college education (pinteraction = 0.02). In AA, increasing MRSSBP was associated with higher SBP in females only (pinteraction = 0.01). Our work shows that MRS is a potential biomarker of blood pressure that may be modified by traditional hypertension risk factors.

Association between alcohol consumption and DNA methylation in blood: a systematic review of observational studies

Aim: We systematically reviewed and evaluated current literature on alcohol consumption and DNA methylation (DNAm) at the genome-wide and probe-wise level in blood of adults. Materials & methods: Five databases (PubMed, Embase, Web of Science, CINAHL and PsycInfo) were searched until 20 December 2020. Studies assessing the effect of alcohol dependence on DNAm were not eligible. Results: 11 cross-sectional studies were included with 88 to 9643 participants. Overall, all studies had a risk of bias criteria unclear or unmet. Epigenome-wide association studies identified between 0 and 5458 differentially methylated positions, and 15 were observed in at least four studies. Conclusion: Potential methylation markers for alcohol consumption have been identified, but further validation in large cohorts is needed.

DNA methylation signatures associated with cardiometabolic risk factors in children from India and The Gambia: results from the EMPHASIS study

BACKGROUND

The prevalence of cardiometabolic disease (CMD) is rising globally, with environmentally induced epigenetic changes suggested to play a role. Few studies have investigated epigenetic associations with CMD risk factors in children from low- and middle-income countries. We sought to identify associations between DNA methylation (DNAm) and CMD risk factors in children from India and The Gambia.

RESULTS

Using the Illumina Infinium HumanMethylation 850 K Beadchip array, we interrogated DNAm in 293 Gambian (7-9 years) and 698 Indian (5-7 years) children. We identified differentially methylated CpGs (dmCpGs) associated with systolic blood pressure, fasting insulin, triglycerides and LDL-Cholesterol in the Gambian children; and with insulin sensitivity, insulinogenic index and HDL-Cholesterol in the Indian children. There was no overlap of the dmCpGs between the cohorts. Meta-analysis identified dmCpGs associated with insulin secretion and pulse pressure that were different from cohort-specific dmCpGs. Several differentially methylated regions were associated with diastolic blood pressure, insulin sensitivity and fasting glucose, but these did not overlap with the dmCpGs. We identified significant cis-methQTLs at three LDL-Cholesterol-associated dmCpGs in Gambians; however, methylation did not mediate genotype effects on the CMD outcomes.

CONCLUSION

This study identified cardiometabolic biomarkers associated with differential DNAm in Indian and Gambian children. Most associations were cohort specific, potentially reflecting environmental and ethnic differences.