Epigenome-Wide Analyses Identify Two Novel Associations With Recurrent Stroke in the Vitamin Intervention for Stroke Prevention Clinical Trial

Genome-wide analysis of DNA methylation and its relationship with serum homocysteine levels in patients with hypertension

BACKGROUND

Homocysteine (Hcy) is an independent risk factor for cardiovascular diseases, and elevated plasma Hcy levels could aggravate vascular injury in hypertension. Hyperhomocysteinemia can change the methylation status of global DNA and specific genes. In the present study, we aim to examine the comprehensive influence of Hcy levels on DNA methylation status in patients with hypertension.

METHODS

Epigenome-wide methylation profiles of the peripheral leukocyte DNA of 218 patients with hypertension were analyzed using the Illumina Infinium Methylation EPIC BeadChip. Differentially methylated positions (DMPs) associated with serum Hcy levels were identified by mixed linear regression with the adjustment of potential confounders. Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis were conducted to determine the potential functions of the identified DMPs. The association between the methylation level of DMPs and carotid-femoral pulse wave velocity (Cf-PWV) was also analyzed.

RESULTS

Five DMPs at cg13169662, cg03179312, cg21976560, cg25262698, and cg09433843 showed significant association with serum Hcy levels (false discovery rate-corrected P  < 0.05). An additional six CpG sites met the threshold for suggestive significance ( P  < 1 × 10 -6 ), among which three DMPs (cg25781123, cg26463106, and cg06679221) were annotated to THUMPD3 . Furthermore, the methylation levels of cg13169662 and cg25262698 (RPRD1A) were significantly associated with Cf-PWV.

CONCLUSION

Our results suggest that Hcy could induce DNA methylation alteration in patients with hypertension. Further functional research is warranted to elucidate the concrete role of DMPs in hypertension.

Plasma Homocysteine (Hcy) Concentration Functions as a Predictive Biomarker of SPECT-Evaluated Post-Ischemic Hyperperfusion in Acute Ischemic Stroke

Introduction

Homocysteine (Hcy) concentration has been reported to be associated with ischemic stroke. In this study, we aimed to investigate the potential of plasma Hcy in the prediction of post-ischemic hyperperfusion in AIS patients, which was diagnosed with the single-photon emission computed tomography (SPECT) method.

Methods

A total of 112 ischemic stroke patients were recruited in this study. According to whether the patients were subjected to post-ischemic hyperperfusion, all recruited subjects were divided into a post-ischemic hyperperfusion (+) group (N=48) and post-ischemic hyperperfusion (-) group (N=64). The basic demographical data, clinicopathological data and laboratory biochemical data were collected and compared. Level of homocysteine (Hcy) and cystatin-C (Cys-C) and their potential as predictive biomarker are also investigated.

Results

No significant differences were spotted between the post-ischemic hyperperfusion group (+) and post-ischemic hyperperfusion (-) group in respect to the basic demographical and clinicopathological data. And the serum Hcy levels were lower in the post-ischemic hyperperfusion (+) group. Moreover, ROC analysis indicated significant relationships between Hcy levels and the onset of post-ischemic hyperperfusion.

Conclusion

In conclusion, we validated that the plasma Hcy concentration can be used as a predictive biomarker of SPECT-evaluated post-ischemic hyperperfusion in patients suffering from acute ischemic stroke.

DNA methylation and cardiovascular disease in humans: a systematic review and database of known CpG methylation sites

BACKGROUND

Cardiovascular disease (CVD) is the leading cause of death worldwide and considered one of the most environmentally driven diseases. The role of DNA methylation in response to the individual exposure for the development and progression of CVD is still poorly understood and a synthesis of the evidence is lacking.

RESULTS

A systematic review of articles examining measurements of DNA cytosine methylation in CVD was conducted in accordance with PRISMA (preferred reporting items for systematic reviews and meta-analyses) guidelines. The search yielded 5,563 articles from PubMed and CENTRAL databases. From 99 studies with a total of 87,827 individuals eligible for analysis, a database was created combining all CpG-, gene- and study-related information. It contains 74,580 unique CpG sites, of which 1452 CpG sites were mentioned in ≥ 2, and 441 CpG sites in ≥ 3 publications. Two sites were referenced in ≥ 6 publications: cg01656216 (near ZNF438) related to vascular disease and epigenetic age, and cg03636183 (near F2RL3) related to coronary heart disease, myocardial infarction, smoking and air pollution. Of 19,127 mapped genes, 5,807 were reported in ≥ 2 studies. Most frequently reported were TEAD1 (TEA Domain Transcription Factor 1) and PTPRN2 (Protein Tyrosine Phosphatase Receptor Type N2) in association with outcomes ranging from vascular to cardiac disease. Gene set enrichment analysis of 4,532 overlapping genes revealed enrichment for Gene Ontology molecular function "DNA-binding transcription activator activity" (q = 1.65 × 10^-11) and biological processes "skeletal system development" (q = 1.89 × 10^-23). Gene enrichment demonstrated that general CVD-related terms are shared, while "heart" and "vasculature" specific genes have more disease-specific terms as PR interval for "heart" or platelet distribution width for "vasculature." STRING analysis revealed significant protein-protein interactions between the products of the differentially methylated genes (p = 0.003) suggesting that dysregulation of the protein interaction network could contribute to CVD. Overlaps with curated gene sets from the Molecular Signatures Database showed enrichment of genes in hemostasis (p = 2.9 × 10^-6) and atherosclerosis (p = 4.9 × 10^-4).

CONCLUSION

This review highlights the current state of knowledge on significant relationship between DNA methylation and CVD in humans. An open-access database has been compiled of reported CpG methylation sites, genes and pathways that may play an important role in this relationship.

Altered methylation pattern in EXOC4 is associated with stroke outcome: an epigenome-wide association study

BACKGROUND AND PURPOSE

The neurological course after stroke is highly variable and is determined by demographic, clinical and genetic factors. However, other heritable factors such as epigenetic DNA methylation could play a role in neurological changes after stroke.

METHODS

We performed a three-stage epigenome-wide association study to evaluate DNA methylation associated with the difference between the National Institutes of Health Stroke Scale (NIHSS) at baseline and at discharge (ΔNIHSS) in ischaemic stroke patients. DNA methylation data in the Discovery (n = 643) and Replication (n = 62) Cohorts were interrogated with the 450 K and EPIC BeadChip. Nominal CpG sites from the Discovery (p value < 10-06) were also evaluated in a meta-analysis of the Discovery and Replication cohorts, using a random-fixed effect model. Metabolic pathway enrichment was calculated with methylGSA. We integrated the methylation data with 1305 plasma protein expression levels measured by SOMAscan in 46 subjects and measured RNA expression with RT-PCR in a subgroup of 13 subjects. Specific cell-type methylation was assessed using EpiDISH.

RESULTS

The meta-analysis revealed an epigenome-wide significant association in EXOC4 (p value = 8.4 × 10-08) and in MERTK (p value = 1.56 × 10-07). Only the methylation in EXOC4 was also associated in the Discovery and in the Replication Cohorts (p value = 1.14 × 10-06 and p value = 1.3 × 10-02, respectively). EXOC4 methylation negatively correlated with the long-term outcome (coefficient = - 4.91) and showed a tendency towards a decrease in EXOC4 expression (rho = - 0.469, p value = 0.091). Pathway enrichment from the meta-analysis revealed significant associations related to the endocytosis and deubiquitination processes. Seventy-nine plasma proteins were differentially expressed in association with EXOC4 methylation. Pathway analysis of these proteins showed an enrichment in natural killer (NK) cell activation. The cell-type methylation analysis in blood also revealed a differential methylation in NK cells.

CONCLUSIONS

DNA methylation of EXOC4 is associated with a worse neurological course after stroke. The results indicate a potential modulation of pathways involving endocytosis and NK cells regulation.

DNA Methylation and Ischemic Stroke Risk: An Epigenome-Wide Association Study

BACKGROUND

Ischemic stroke (IS) risk heritability is partly explained by genetics. Other heritable factors, such as epigenetics, could explain an unknown proportion of the IS risk. The objective of this study is to evaluate DNA methylation association with IS using epigenome-wide association studies (EWAS).

METHODS

 We performed a two-stage EWAS comprising 1,156 subjects. Differentially methylated positions (DMPs) and differentially methylated regions (DMRs) were assessed using the Infinium 450K and EPIC BeadChip in the discovery cohort (252 IS and 43 controls). Significant DMPs were replicated in an independent cohort (618 IS and 243 controls). Stroke subtype associations were also evaluated. Differentially methylated cell-type (DMCT) was analyzed in the replicated CpG sites using EpiDISH. We additionally performed pathway enrichment analysis and causality analysis with Mendelian randomization for the replicated CpG sites.

RESULTS

 A total of 957 CpG sites were epigenome-wide-significant (p ≤ 10^-7) in the discovery cohort, being CpG sites in the top signals (logFC = 0.058, p = 2.35 × 10^-22; logFC = 0.035, p = 3.22 × 10^-22, respectively). ZFHX3 and MAP3K1 were among the most significant DMRs. In addition, 697 CpG sites were replicated considering Bonferroni-corrected p-values (p < 5.22 × 10^-5). All the replicated DMPs were associated with risk of cardioembolic, atherothrombotic, and undetermined stroke. The DMCT analysis demonstrated that the significant associations were driven by natural killer cells. The pathway enrichment analysis showed overrepresentation of genes belonging to certain pathways including oxidative stress. ZFHX3 and MAP3K1 methylation was causally associated with specific stroke-subtype risk.

CONCLUSION

 Specific DNA methylation pattern is causally associated with IS risk. These results could be useful for specifically predicting stroke occurrence and could potentially be evaluated as therapeutic targets.

Association of DNA Methylation Patterns in 7 Novel Genes With Ischemic Stroke in the Northern Chinese Population

Background: Ischemic stroke is a highly complex disorder. This study aims to identify novel methylation changes in ischemic stroke.

Methods: We carried out an epigenome-wide study of ischemic stroke using an Infinium HumanMethylation 850K array (cases:controls = 4:4). 10 CpG sites in 8 candidate genes from gene ontology analytics top-ranked pathway were selected to validate 850K BeadChip results (cases:controls = 20:20). We further qualified the methylation level of promoter regions in 8 candidate genes (cases:controls = 188:188). Besides, we performed subgroup analysis, dose-response relationship and diagnostic prediction polygenic model of candidate genes.

Results: In the discovery stage, we found 462 functional DNA methylation positions to be associated with ischemic stroke. Gene ontology analysis highlighted the “calcium-dependent cell-cell adhesion via plasma membrane cell adhesion molecules” item, including 8 candidate genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB6/PCDHB9). In the replication stage, we identified 5 differentially methylated loci in 20 paired samples and 7 differentially methylated genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9) in 188 paired samples. Subgroup analysis showed that the methylation level of above 7 genes remained significantly different in the male subgroup, large-artery atherosclerosis subgroup and right hemisphere subgroup. The methylation level of each gene was grouped into quartiles, and Q4 groups of the 7 genes were associated with higher risk of ischemic stroke than Q1 groups (p &lt; 0.05). Besides, the polygenic model showed high diagnostic specificity (0.8723), sensitivity (0.883), and accuracy (0.8777).

Conclusion: Our results demonstrate that DNA methylation plays a crucial part in ischemic stroke. The methylation of these 7 genes may be potential diagnostic biomarker for ischemic stroke.

DNA methylation analyses identify an intronic ZDHHC6 locus associated with time to recurrent stroke in the Vitamin Intervention for Stroke Prevention (VISP) clinical trial

Aberrant DNA methylation profiles have been implicated in numerous cardiovascular diseases; however, few studies have investigated how these epigenetic modifications contribute to stroke recurrence. The aim of this study was to identify methylation loci associated with the time to recurrent cerebro- and cardiovascular events in individuals of European and African descent. DNA methylation profiles were generated for 180 individuals from the Vitamin Intervention for Stroke Prevention clinical trial using Illumina HumanMethylation 450K BeadChip microarrays, resulting in beta values for 470,871 autosomal CpG sites. Ethnicity-stratified survival analyses were performed using Cox Proportional Hazards regression models for associations between each methylation locus and the time to recurrent stroke or composite vascular event. Results were validated in the Vall d’Hebron University Hospital cohort from Barcelona, Spain. Network analyses of the methylation loci were generated using weighted gene coexpression network analysis. Primary analysis identified four significant loci, cg04059318, ch.2.81927627R, cg03584380, and cg24875416, associated with time to recurrent stroke. Secondary analysis identified three loci, cg00076998, cg16758041, and cg02365967, associated with time to composite vascular endpoint. Locus cg03584380, which is located in an intron of ZDHHC6, was replicated in the Vall d’Hebron University Hospital cohort. The results from this study implicate the degree of methylation at cg03584380 is associated with the time of recurrence for stroke or composite vascular events across two ethnically diverse groups. Furthermore, modules of loci were associated with clinical traits and blood biomarkers including previous number of strokes, prothrombin fragments 1 + 2, thrombomodulin, thrombin-antithrombin complex, triglyceride levels, and tissue plasminogen activator. Ultimately, these loci could serve as potential epigenetic biomarkers that could identify at-risk individuals in recurrence-prone populations.

Multi-omic analysis of stroke recurrence in African Americans from the Vitamin Intervention for Stroke Prevention (VISP) clinical trial

African Americans endure a nearly two-fold greater risk of suffering a stroke and are 2–3 times more likely to die from stroke compared to those of European ancestry. African Americans also have a greater risk of recurrent stroke and vascular events, which are deadlier and more disabling than incident stroke. Stroke is a multifactorial disease with both heritable and environmental risk factors. We conducted an integrative, multi-omic study on 922 plasma metabolites, 473,864 DNA methylation loci, and 556 variants from 50 African American participants of the Vitamin Intervention for Stroke Prevention clinical trial to help elucidate biomarkers contributing to recurrent stroke rates in this high risk population. Sixteen metabolites, including cotinine, N-delta-acetylornithine, and sphingomyelin (d17:1/24:1) were identified in t-tests of recurrent stroke outcome or baseline smoking status. Serum tricosanoyl sphingomyelin (d18:1/23:0) levels were significantly associated with recurrent stroke after adjusting for covariates in Cox Proportional Hazards models. Weighted Gene Co-expression Network Analysis identified moderate correlations between sphingolipid markers and clinical traits including days to recurrent stroke. Integrative analyses between genetic variants in sphingolipid pathway genes identified 29 nominal associations with metabolite levels in a one-way analysis of variance, while epigenomic analyses identified xenobiotics, predominately smoking-associated metabolites and pharmaceutical drugs, associated with methylation profiles. Taken together, our results suggest that metabolites, specifically those associated with sphingolipid metabolism, are potential plasma biomarkers for stroke recurrence in African Americans. Furthermore, genetic variation and DNA methylation may play a role in the regulation of these metabolites.

Detecting methylation signatures in neurodegenerative disease by density-based clustering of applications with reducing noise

There have been numerous genetic and epigenetic datasets generated for the study of complex disease including neurodegenerative disease. However, analysis of such data often suffers from detecting the outliers of the samples, which subsequently affects the extraction of the true biological signals involved in the disease. To address this critical issue, we developed a novel framework for identifying methylation signatures using consecutive adaptation of a well-known outlier detection algorithm, density based clustering of applications with reducing noise (DBSCAN) followed by hierarchical clustering. We applied the framework to two representative neurodegenerative diseases, Alzheimer's disease (AD) and Down syndrome (DS), using DNA methylation datasets from public sources (Gene Expression Omnibus, GEO accession ID: GSE74486). We first applied DBSCAN algorithm to eliminate outliers, and then used Limma statistical method to determine differentially methylated genes. Next, hierarchical clustering technique was applied to detect gene modules. Our analysis identified a methylation signature comprising 21 genes for AD and a methylation signature comprising 89 genes for DS, respectively. Our evaluation indicated that these two signatures could lead to high classification accuracy values (92% and 70%) for these two diseases. In summary, this framework will be useful to better detect outlier-free genetic and epigenetic signatures in various complex diseases and their developmental stages.

Comparison of DNA Methylation Profiles of Hemostatic Genes between Liver Tissue and Peripheral Blood within Individuals

DNA methylation has become increasingly recognized in the etiology of complex diseases, including thrombotic disorders. Blood is often collected in epidemiological studies for genotyping and has recently also been used to examine DNA methylation in epigenome-wide association studies. DNA methylation patterns are often tissue-specific, thus, peripheral blood may not accurately reflect the methylation pattern in the tissue of relevance. Here, we collected paired liver and blood samples concurrently from 27 individuals undergoing liver surgery. We performed targeted bisulfite sequencing for a set of 35 hemostatic genes primarily expressed in liver to analyze DNA methylation levels of >10,000 cytosine-phosphate-guanine (CpG) dinucleotides. We evaluated whether DNA methylation in blood could serve as a proxy for DNA methylation in liver at individual CpGs. Approximately 30% of CpGs were nonvariable and were predominantly hypo- (<25%) or hypermethylated (>70%) in both tissues. While blood can serve as a proxy for liver at these CpGs, the low variability renders these unlikely to explain phenotypic differences. We therefore focused on CpG sites with variable methylation levels in liver. The level of blood-liver tissue correlation varied widely across these variable CpGs; moderate correlations (0.5 ≤ r < 0.75) were detected for 6% and strong correlations (r ≥ 0.75) for a further 4%. Our findings indicate that it is essential to study the concordance of DNA methylation between blood and liver at individual CpGs. This paired blood-liver dataset is intended as a resource to aid interpretation of blood-based DNA methylation results.

Impact of homocysteine levels on clinical outcome in patients with acute ischemic stroke receiving intravenous thrombolysis therapy

Background

The purpose of this study was to retrospectively assess the potential correlation between clinical outcomes and homocysteine (Hcy) levels in acute ischemic stroke (AIS) patients after recombinant tissue plasminogen activator (rtPA) treatment.

Methods

AIS patients treated by rtPA were enrolled between September 2018 and March 2019 in the Stroke Center (Department of Neurology and Neurosurgery), Shanghai Tenth People’s Hospital, Tongji University School of Medicine. Demographics, baseline and clinical characteristics, and modified Rankin Scale (mRS) score after three months from the onset were retrospectively analyzed. Then we compared data about demographics, baseline and clinical characteristics between patients with favorable (mRS score 0–2) and unfavorable (mRS score 3–6) outcomes.

Results

Among 141 patients, 36 patients had poor outcome, for an incidence of 25.53%. Univariate analysis showed that higher Hcy levels (OR = 1.07, 95% CI [1.02–1.12]), older age (OR = 1.06, 95% CI [1.02–1.10]), longer door to needle time (DNT) (OR = 1.03, 95% CI [1.01–1.05]), higher D-Dimer levels (OR = 1.33, 95% CI [1.03–1.71]), and higher National Institutes of Health Stroke Scale (NIHSS) score before treatment (OR = 1.21, 95% CI [1.08–1.35]) were each associated with poor outcome. Also, without internal carotid artery plaque (OR = 0.30, 95% CI [0.10–0.92]) showed a protective effect on patients’ clinical outcome. Patients with higher levels of Hcy decline also showed an increased risk of poor outcome for AIS patients obtaining rtPA treatment (Non-adjusted: OR = 1.07, 95% CI [1.02–1.12]; Adjust model I adjusts for demographics (age, male): OR = 1.06, 95% CI [1.02–1.11]; Adjust model II adjusts for hospital care factors (onset to treatment, DNT): OR = 1.08, 95% CI [1.03–1.13]; Adjust model III adjusts for health and stroke factors (INR, D-Dimer, HGB, NIHSS score before treatment, smoking, drinking, hypertension, diabetes, coronary disease, hyperlipidemia, previous stroke, atrial fibrillation, hemorrhagic transformation, internal carotid artery plaque): OR = 1.06, 95% CI [1.02–1.11]). The results are very stable in all three models constructed.

Conclusion

The results of this study indicate that increased Hcy level independently predicts unfavorable outcome in AIS patients accepting thrombolytic therapy. However, the contribution of Hcy to the outcome, although significant, is relatively small and perhaps not clinically significant when all the other confounders are considered.