- Research
- Open access
- Published:
Associations of NPPA promoter true methylation and hydroxymethylation with ischemic stroke and its functional outcome
Epigenetics Communications volume 4, Article number: 6 (2024)
Abstract
Background
Atrial natriuretic peptide (ANP) has been associated with ischemic stroke (IS), but the underlying molecular mechanisms are not clear. The objective of this study was to examine whether DNA methylation and hydroxymethylation in the coding gene of ANP (NPPA) were associated with IS, as well as its functional outcome.
Methods and results
DNA methylation and hydroxymethylation in NPPA promoter region were quantified by targeted bisulfite sequencing and APOBEC-coupled epigenetic sequencing, respectively, in 615 IS patients and 610 age- and sex-matched healthy controls. True methylation was calculated as the quantified methylation level minus the quantified hydroxymethylation level. The functional outcome of IS was evaluated by the modified Rankin Scale (mRS) at a 3-month follow-up visit after onset. Multiple testing was controlled by the false discovery rate approach. Of the nine CpG sites assayed, hypomethylation at eight CpGs was not only associated with an increased risk of having IS but also predicted a higher mRS score after onset (all q < 0.05) and hydroxymethylation at one CpG located at Chr1:11908348 was positively associated with IS (OR = 1.39, 95%CI:1.15–1.68, q < 0.05), after adjusting for covariates and multiple testing. While, hydroxymethylation levels at none of the CpGs assayed were significantly associated with the mRS score.
Conclusions
DNA methylation and hydroxymethylation in NPPA promoter were associated with the risk of IS and its functional outcome in Chinese adults.
Background
Atrial natriuretic peptide (ANP), a hormone synthesized and released mainly from myocytes with pleiotropic cardiovascular and metabolic properties [1], plays a critical role in hemodynamic homeostasis through natriuresis, diuresis, and vasodilation [2] and cardiovascular remodeling due to its antiproliferative, antifibrotic, and antiangiogenic functions [3]. All of these processes constitute the physiopathology of ischemic stroke (IS) [4]. Admittedly, circulating ANP levels have been associated with IS [5] as well as related cardiovascular disorders, e.g., hypertension [6], heart failure [7], and myocardial infarction [8]. In addition to IS development, plasma ANP levels could also predict 90-day mortality and functional outcomes after IS onset [9, 10]. However, the clinical implementation of ANP in the management of heart failure with ANP-targeted agents, such as sacubitril/valsartan [11] and recombinant carperitide [12, 13] has encountered significant challenges. These include the occurrence of hypotension, angioedema, and fatalities in hospital settings. A better understanding of the molecular mechanisms underlying the association between ANP and the occurrence and development of IS may help its drug development and application in the prevention and control of this debilitating disorder.
As an interface between the fixed genome and dynamic environment, epigenetic modifications such as DNA methylation regulate gene expression without changing the DNA sequence [14], thereby participating in the mechanisms of complex phenotypes. For example, epigenome-wide association studies and candidate gene studies have identified many methylation markers of IS [15,16,17], as well as its prognosis [18, 19], with a lack of genome-widely significant epigenetic markers related to ANP. Considering the contribution of genetic variants of NPPA, the coding gene of ANP, to the susceptibility of IS [20, 21], we hypothesized that DNA methylation in the NPPA promoter region may explain, at least partially, the molecular mechanisms underlying the role of ANP in IS development and prognosis. Further, our previous studies have already found that NPPA promoter hypomethylation predicted an increased risk of cardiovascular disease (CVD) [22], central obesity [23], and hypertension [24].
DNA methylation is a dynamic process in response to the environment. DNA methyltransferase binds a methyl group to cytosine 5’ carbon atom at a CpG site forming 5-methylcytosine (5mC) [25]. When DNA methylation occurs at CpG islands of the gene promoter region, it is normally associated with the repression of gene expression [26]. Under the action of the TET enzyme, 5mC can be further oxidized to 5-hydroxymethyl cytosine (5hmC) by incorporating a single oxygen atom from molecular oxygen into its substrate, 5-methylcytosine [27]. This process is hydroxymethylation, a known demethylation intermediate that may upregulate gene expression [28]. The demethylation process also includes the oxidation of 5hmC to 5-formylcytosine (5fC), 5-carboxylcytosine (5caC), and further conversion to an unmodified cytosine by replication-dependent dilution or thymine DNA glycosylase (TDG)-dependent base excision repair [29]. Of these epigenetic modifications, 5mC and 5hmC are the most common but have conflicting functions on the regulation of gene expression [30]. Nevertheless, DNA methylation levels are currently quantified based on bisulfite sequencing which cannot discriminate between 5mC and 5hmC [31]. Our previous study has found that NPPA promoter methylation assayed using such a method was associated with IS [32], but the methylation level assayed actually reported a composite level including both 5mC and 5hmC and may therefore be prone to information bias. Moreover, mice experiments have demonstrated that 5hmC could protect the brain from ischemic injury [33, 34]. DNA hydroxymethylation of the NPPA gene may also be associated with the development and prognosis of IS, but this association was scarcely studied in humans. Therefore, we quantified the levels of DNA methylation and hydroxymethylation simultaneously at identical CpG sites in the NPPA promoter region for patients with IS and age- and sex-matched healthy controls. We then calculated true methylation levels by subtracting the hydroxymethylation level from the methylation level at corresponding site to examine whether NPPA promoter true methylation and hydroxymethylation are associated with IS, as well as its functional outcome after IS.
Results
Baseline characteristics of study participants
A total of 615 patients with IS (mean aged 62 years, 47% males) and 610 healthy controls (mean aged 62 years, 48% males) were included in the current study. Their baseline characteristics are presented in Table 1. As expected, patients with IS had more metabolic risk factors including hypertension, diabetes, obesity, hyperglycemia, and dyslipidemia than healthy controls (all P < 0.05). Additionally, the difference in education level was also statistically significant between two groups (P < 0.05).
Association between NPPA promoter true methylation and IS
As shown in eTable 1, eight of the nine CpG sites assayed were hypomethylated, as suggested by true methylation levels, in patients with IS, in comparison to healthy controls (all q < 0.05). The single CpG association analysis using binary logistic regression with log2-transformed methylation level as the independent variable showed that true methylation levels at eight CpG sites (except for CpG9) were negatively associated with the risk of having IS (all P < 0.05), after adjusting for age, sex, education level, current smoking, current drinking, BMI, LDL-C, HDL-C, hypertension, and diabetes (Fig. 1). After further correction of multiple testing, these associations were still statistically significant (all q < 0.05).
A forest plot showing the single-CpG associations of true methylation and hydroxymethylation of NPPA promoter with ischemic stroke. The black squares and line segments illustrate the ORs and 95% CIs of having ischemic stroke associated with DNA methylation and hydroxymethylation levels (after log2-transformation), adjusting for age, sex, education level, current smoking, current drinking, BMI, LDL-C, HDL-C, hypertension, and diabetes. q value indicates the significance level after correction of multiple testing by false discovery rate approach. CI: confidence interval; OR: odds ratio
By combing P-values of the single CpG associations, wTPM revealed significant gene-based associations between true methylation of NPPA promoter and IS (P < 0.05), suggesting that DNA methylation of the nine CpG sites may be jointly associated with IS. The impact of the mean of true methylation levels on the risk of IS was visualized in Fig. 2a, the mean methylation level was also significantly associated with IS in the same direction as the single CpG associations (OR = 0.94, 95%CI:0.91–0.97, P < 0.001). In subgroup analyses, the association between the mean of true methylation levels of the nine CpG sites and IS was more pronounced among elders, females, non-smokers, non-drinkers, and those with a low level of education (eTable 2), which may have practical guiding significance for the clinical application of epigenetic markers. Further, we did not find a statistically significant group difference regarding metabolic risk factors of IS, such as obesity and hypertension, which have been associated with NPPA methylation in our previous studies. The results of the subgroup analysis suggested a robust association between NPPA true methylation and IS and this association may not be driven by metabolic factors.
Cubic spline curves visualizing the impacts of the average levels of true methylation and hydroxymethylation in the NPPA promoter region on ischemic stroke. Red lines indicate odds ratios (solid) and their 95% confidence intervals (dotted) of ischemic stroke associated with the mean levels of true methylation (a) and hydroxymethylation (b), adjusting for age, sex, education level, current smoking, current drinking, BMI, LDL-C, HDL-C, hypertension, and diabetes. The histogram shows the distribution of participants according to increasing levels of true methylation and hydroxymethylation. CI: confidence interval; OR: odds ratio
Association between NPPA promoter hydroxymethylation and IS
In contrast to methylation, the median levels of hydroxymethylation at the nine CpG sites assayed were generally higher in patients with IS than that in healthy controls, although no statistically significant group-difference was found (eTable 1). After adjusting for covariates, hydroxymethylation level at one CpG site (CpG2, located in Chr1:11908348) was positively associated with IS (OR = 1.39, 95%CI: 1.15–1.68, P < 0.05 after log2-transformation, Fig. 1). This association persisted after further correction of multiple testing (q < 0.05). The wTPM revealed a significant gene-based association (P < 0.05), suggesting that NPPA promoter hydroxymethylation as a whole was significantly associated with IS. Nevertheless, we failed to observe a significant association between the mean of hydroxymethylation levels of the nine CpG sites and IS (Fig. 2b and eTable 3).
Association between NPPA promoter true methylation and the functional outcome of IS
Of the 615 IS patients with available true methylation data of NPPA promoter at baseline, 17 patients were lost of follow-up and the 598 patients remained were included in the prospective analysis to examine whether true methylation of NPPA promoter was associated with the functional outcome of IS. Ordered logistic regression found that the true methylation levels at eight CpG sites assayed at baseline were significantly associated with a lower risk of having a higher mRS score at the 3-month visit, after multivariate adjustment and multiple testing correction (all q < 0.05, Table 2). To ease data interpretation, binary logistic regression with poor functional outcome defined as mRS ≥ 3 as the dependent variable found similar associations (eFigure 1). Gene-based association between true methylation of NPPA promoter and mRS (all P < 0.05) was also significant, suggesting that DNA methylation of the nine CpG sites may be jointly associated with functional outcome of IS.
Association between NPPA promoter hydroxymethylation and the functional outcome of IS
DNA hydroxymethylation levels at none of the CpG sites assayed were significantly associated with the mRS (Table 2) and the poor functional outcome (eFigure 1) at the 3-month follow-up visit. In addition, no gene-based association was found that NPPA promoter hydroxymethylation as a whole was jointly associated with functional outcome of IS (P = 0.453).
Discussion
In a large sample of patients with IS and healthy controls, we simultaneously examined for the first time whether DNA true methylation and hydroxymethylation in NPPA promoter were associated with IS as well as its functional outcome. We found that NPPA promoter hypomethylation was not only associated with an increased risk of having IS but also predicted poorer neurologic recovery at 3-month follow up visit after onset and hydroxymethylation was elevated in patients with IS but failed to predict the functional outcome of IS, independent of conventional risk factors. These results suggest that DNA methylation and hydroxymethylation in NPPA promoter may participate in the development and prognosis of IS. The identified CpG sites could be used as epigenetic markers and even therapeutic targets for IS.
In line with our study, the role of ANP and its coding gene in IS has been extensively studied. For example, previous studies have identified ANP as a potential biomarker for early detection of IS [35] and a therapeutic target for disease due to its biological functions of natriuresis, diuresis, vasorelaxation, and inhibition of renin and aldosterone secretion [2, 36]. Researches on the role of ANP in the pathogenesis of cerebrovascular diseases have been conducted in great depth, yet the clinical application of ANP-related drugs is still very limited [11,12,13]. Genetic polymorphisms in the NPPA gene could not only determine the circulating levels of ANP [37] but also indicate the susceptibility of IS [20, 21]. A genetic association study including 2660 participants revealed that the NPPA rs5063 was associated with a reduced risk for having IS [20]. Another study conducted in a Korean population found that rs5065 and rs5067 of the NPPA gene were associated with the susceptibility for cardioembolic stroke [21]. As an interface between the fixed genome and dynamic environment, DNA methylation in the NPPA gene may affect its function and subsequent ANP synthesis and excretion. DNA methylation in NPPA promoter may explain the underlying molecular mechanisms of the relationship between ANP and IS as well as its functional outcome. Importantly, since such epigenetic modifications are reversible processes, they could assist the development of new therapeutic approaches for the treatment of human diseases [38].
In recent years, with the continuous development of technology, the correlation between DNA methylation markers and IS as well as its functional outcome have been gradually identified. In line with our study, animal studies have found that dysregulated DNA methylation could cause the aggravation of central ischemic injury in mice [39, 40]. In humans, methylation markers at some candidate genes such as SLC6A4 [41] and CRY1 [42] were related to the long-term adverse outcomes in patients with IS. A comprehensive review suggested that and DNA methylation at hundreds of methylation markers, in addition to global methylation, covering a wide range of genes but none was related to the NPPA gene, were associated with IS [43]. However, DNA methylation levels reported by these studies were measured utilizing bisulfite treatment-based methods, which cannot distinguish between the methylated and hydroxymethylated forms of cytosine [44]. In our study, we calculated true methylation levels by quantifying hydroxymethylation levels at the same CpG sites. We are the first, to the best of our knowledge, to examine and found that DNA true methylation in NPPA promoter was associated with IS and its functional outcome. In our study, participants with IS were tended to have lower true methylation levels, and hypomethylation was associated with increased risk of having IS and poor neurologic recovery. It is well known that methylation status at the gene promoter region is inversely correlated to the gene expression level generally [45], and many studies have also found that participants with IS had higher levels of serum ANP than those without and the severity of the disease could be distinguished in the light of ANP levels [46,47,48]. Our results were consistent with the findings of existing studies, though the causal association between NPPA promoter methylation and ANP expression level has not been verified yet.
As an important member of epigenetic modifications, DNA hydroxymethylation has conflicting functions of regulating gene expression compared with DNA methylation so that we cannot ignore its role in the occurrence and development of disease. Therefore, we further explored the associations of NPPA promoter hydroxymethylation with IS and its functional outcome. The prevailing evidence has shown that DNA hydroxymethylation is involved in the pathogenesis and prognosis of IS while the elucidation of the link between them is still in its infancy, especially in humans. In animal models of focal ischemia, TET2/3, the enzyme that catalyzed the formation of 5hmC, knockdown further increased brain degeneration in mice, demonstrating a role of 5hmC in endogenous protection against stroke [33, 49]. In humans, there was only a small sample of 24 participants showed a significant increase in the abundance of 5hmC in blood samples from patients with IS [33]. In our study, we found that NPPA promoter hydroxymethylation levels at only one CpG site (CpG2 located at Chr1:11908348) seemed to be positively associated with the risk of having IS, revealing the possibility of DNA hydroxymethylation in NPPA promoter as a potential biomarker and a therapeutic target for IS.
In this study, it can be seen that methylation levels and hydroxymethylation levels differed substantially between CpG sites assayed, and their effects on the IS were also different. Of the nine CpG sites assayed, one CpG site (CpG2 located at Chr1:11908348) exhibited a strong correlation with IS, as both methylation level and hydroxymethylation level of this site were significantly associated with IS. This CpG site is highly expected to be a precise intervention target for IS if the relevant technology is feasible in the future. Nevertheless, the methylation and hydroxymethylation levels at these CpG sites were highly correlated. In addition, the contribution of methylation or hydroxymethylation at a single CpG site to the complex phenotypes was generally small [50], and some CpG sites may be overlooked due to a lack of statistical significance before and after multiple testing corrections. Thus, we further tested the combined effects of multiple CpG sites and found that NPPA promoter methylation and hydroxymethylation as a whole were jointly associated with IS. Our research suggested that testing the combined effects of multiple CpG sites was another important approach in epigenetic analysis for complex diseases.
Several strengths in our study deserved to be mentioned. First, careful and systematic analyses were made to examine the associations of true methylation and hydroxymethylation in NPPA promoter with IS and its functional outcome. In addition, associations of NPPA promoter methylation and hydroxymethylation with pathogenesis and prognosis of IS were taken into consideration simultaneously in this article, providing a comprehensive basis for further development of the clinical application of ANP. Some limitations also characterized the present study. First, an observational study is incapable of deducing whether aberrant methylation and hydroxymethylation in NPPA promoter is a risk factor, outcome or just a concomitant phenomenon of IS. Second, although our results remained significant after adjusting for a variety of known risk factors, residual confounding may still exist. Third, the generalizability of our results to other ethnic populations or younger populations is unknown, since our findings were derived from Chinese adults whose cardiovascular profiles could be different from other populations with different ethnic backgrounds or other age groups. Fourth, although our results supported that DNA methylation and hydroxymethylation in NPPA promoter might explain the molecular mechanisms underlying the relationship between ANP and IS as well as its functional outcome, whether or to what extent the contribution of ANP to IS could be explained is unknown due to lacking data on circulating ANP levels. In our previous study [22], we examined the associations between methylation levels at the same CpG sites of the NPPA gene and serum proANP, but it is with regret that no CpG sites reached statistical significance. Although the association between NPPA promoter methylation and ANP has not been confirmed yet, many studies have already found that the variants of the NPPA gene were associated with circulating levels of ANP.
Conclusions
Levels of DNA methylation and hydroxymethylation in NPPA promoter was associated with the risk of IS and its functional outcome in Chinese adults. Dysregulated methylation and hydroxymethylation in the NPPA gene may participate in the development and prognosis of IS. Due to epigenetic modifications are modifiable, the CpG sites identified may be epigenetic markers and even potential therapeutic targets for IS.
Methods
Study participants
As described in Fig. 3, cases of IS and their healthy controls included in the current case-control study were selected from the China Antihypertensive Trial in Ischemic Stroke (CATIS) and the Prevention of Metabolic syndrome and multi-metabolic disorders Study (PMMS), respectively. The detailed methods of the selection of study participants have been detailed elsewhere [32]. In brief, the CATIS was a multicenter randomized controlled clinical trial designed to test the effect of blood pressure reduction on the prognosis of IS. A total of 4071 patients with IS confirmed by brain tomography or MAGNETIC resonance imaging were recruited. Among them, 3013 patients had DNA samples and 1000 patients were randomly selected as candidate cases. A total of 1000 healthy controls were randomly selected by frequency matching based on age (± 4 years) and sex with cases from 3999 participants with available DNA samples in the PMMS, which was a community-based prospective cohort study of CVD and its risk factors in Chinese adults. Of the 2000 candidate participants whose DNA samples were assayed, 615 IS cases and 610 healthy controls with successful quantification of hydroxymethylation and further calculation of true methylation in NPPA promoter were finally included in the current analysis, after excluding samples that failed in quantification of methylation and hydroxymethylation.
A flowchart illustrating the study design and selection of study participants. A case-control analysis was performed in 615 patients with ischemic stroke and 610 healthy controls to examine the associations of true methylation and hydroxymethylation of NPPA promoter with ischemic stroke. Leveraging the follow-up data for the patients with ischemic stroke, a prospective analysis was performed to examine whether true methylation and hydroxymethylation of NPPA promoter predicted the functional outcome after ischemic stroke onset
Standard protocol approvals
The protocols of the CATIS were approved by the institutional review boards at Tulane University in the United States and Soochow University in China, as well as ethical committees at the participating hospitals. The protocols of the PMMS were approved by the Soochow University Ethics Committee. Written informed consent was obtained from all study participants or their immediate family members in both the CATIS and the PMMS.
Quantification of NPPA promoter methylation
As illustrated in Fig. 4, targeted bisulfite sequencing [51] was applied to quantify DNA methylation levels in the promoter region of the NPPA gene. The targeted region was selected from Chr1:11908380 to Chr1:11907840, relative to the transcriptional start site − 540 bp to -270 bp. Primers were carefully designed based on the targeted sequence coordinated in the Genome Reference Consortium Human Build 37 (GRCh37). The genomic DNA, extracted from peripheral blood mononuclear cells, was first bisulfite treated using the EZ DNA Methylation-Gold Kit (Zymo Research, Inc., CA, USA) and then amplified by polymerase chain reaction (PCR). Using the paired-end sequencing protocol according to the manufacturer’s guidelines, the PCR products were ultimately sequenced by Illumina Hiseq 2000 (Illumina, Inc., CA, USA). After bisulfite treatment, unmethylated cytosines were deaminated and converted to uracil while methylated cytosines, containing traces of hydroxymethylated cytosines, were unchanged. The methylation level at each CpG site was calculated as the percentage of methylated alleles over the total alleles. In the aspect of quality control, we took the percentage of converted non-CpG cytosines as the conversion rate of each sample and excluded samples whose conversion rate was lower than 98%. In addition, samples with an average coverage of less than 20× were further excluded.
A schematic illustration of the targeted sequence, primers, and quantification of DNA methylation and hydroxymethylation. A total of 9 CpG sites (red colored) were finally assayed in the NPPA gene promoter (-540 ~ -276Â bp from TSS). ACE-seq: APOBEC-coupled epigenetic sequencing; BS-seq: Bisulfite sequencing; TSS: transcriptional start site; 5mC: 5-methylcytosine; 5hmC: 5-hydroxymethyl cytosine; 5gmC: 5hmc labeled with glucose; C: cytosine; U: uracil; T: thymine
Quantification of NPPA promoter hydroxymethylation
Using the same DNA samples, we additionally quantified the hydroxymethylation levels at the same CpG sites with APOBEC-coupled sequencing (Fig. 4). With 0.2% synthetic full methylated and full hydroxymethylated control DNA spiked in, the genomic DNA was treated using T4 β-glucosyltransferase (β-GT, NEB) and APOBEC3A. β-GT was used for selective labeling of hydroxymethylated cytosines with glucose which can prevent them from being deaminated by APOBEC3A. As a result, this enzymatic conversion converts cytosines and methylated cytosines into uracils and leaves hydroxymethylated cytosines unchanged. Using the paired-end sequencing protocol according to the manufacturer’s guidelines, the enzyme-treated samples were then sequenced by Illumina Hiseq 2000 (Illumina, Inc., CA, United States). The hydroxymethylation level at each CpG site was calculated as the percentage of the hydroxymethylated alleles. The same procedures of quality control as used in methylation quantification were applied.
Calculation of NPPA promoter true methylation
In the methods of assessment for methylation, bisulfite conversion does not allow for distinguishing between methylation and hydroxymethylation. As a result, the amounts of methylated alleles detected in fact comprise hydroxymethylated alleles. Therefore, the true methylation level at each CpG site could be calculated as the corresponding quantified methylation level minus the quantified hydroxymethylation level.
Assessment of functional outcome of IS
The modified Rankin Scale (mRS) was applied to assess the functional outcome of patients with IS at 3 months after onset. It is rated on a scale of 0–6, with 0–2 indicating no symptoms to a slight disability, 3–5 indicating moderate to severe disability, and 6 indicating death. The mRS score was evaluated by trained neurologic physicians blinded to the baseline characteristics of the patients with IS.
Collection of risk factors for IS
As detailed in eMethods (supplementary data), demographic data (age, sex, education level), lifestyle factors (current smoking, current drinking), disease history, body mass index (BMI), blood pressure, fasting glucose, and blood lipids, including total cholesterol, triglycerides, low-density lipoprotein cholesterol (LDL-C), and high-density lipoprotein cholesterol (HDL-C), were collected for all study participants. For patients with IS, subtypes including thrombotic, embolism, lacunar infarction, and other subtypes were assessed based on the symptoms and imaging data [52]. The severity of IS was assessed by trained neurologists using the National Institutes of Health Stroke Scale (NIHSS) [53].
Statistical analysis
All statistical analyses were performed using R version 4.1.3. The baseline characteristics of study participants were presented in patients with IS and healthy controls, respectively. Log2-transformation was applied to maximize the normality of epigenetic data and the generated values were used in logistic regression analyses.
Analysis of the associations of NPPA true methylation and hydroxymethylation with IS
To examine whether true methylation of NPPA promoter was associated with IS, we first conducted a single CpG association analysis followed by a gene-based association analysis. In specific, the single CpG association analysis was performed by constructing a binary logistic regression model in which IS (y/n) was the dependent variable and true methylation level at each CpG site was the independent variable, adjusting for age, sex, education level, current smoking, current drinking, BMI, LDL-C, HDL-C, hypertension, and diabetes. Multiple testing was controlled by adjusting for the total number of CpG loci tested using the false discovery rate (FDR) approach, and an FDR-adjusted P value (i.e., q value) of less than 0.05 was considered statistically significant. Based on the results of the single CpG association analysis, the weighted truncated product method (wTPM) [54] was used to test the gene-based association between true methylation at multiple CpG sites and IS. This method combined P values of all CpGs that reached a preselected threshold (e.g., raw P < 0.1 in this study). The regression coefficient of each individual CpG methylation was included as weights in the wTPM statistic. To visualization of the effect of DNA methylation on IS, the mean of true methylation levels at multiple CpG sites was considered a surrogate of the true methylation level of the target region and was also used to examine the gene-based association by constructing a restricted cubic spline regression model with 4 knots (5th, 35th, 65th, and 95th percentiles) [55]. The association between DNA hydroxymethylation and IS was similarly examined.
Analysis of the associations of NPPA true methylation and hydroxymethylation with the functional outcome of IS
The associations of true methylation and hydroxymethylation of NPPA promoter with the functional outcome of IS were similarly examined. In brief, an ordered logistic regression model with mRS score at 3-month follow-up as an ordinal response variable was constructed. In addition to the variables listed above, ischemic stroke subtypes, hours from onset to hospitalization, NIHSS score, and treatment group at admission were also adjusted for in this model. To ease data interpretation, the associations of true methylation and hydroxymethylation of NPPA promoter with poor functional outcome defined as mRS score ≥ 3 were also examined by constructing a binary logistic regression model.
Data availability
The data that supports the finding of this study was publicly available from Dryad: https://doiorg.publicaciones.saludcastillayleon.es/10.5061/dryad.h18931zvn. This link is permanent and will be active once the paper is published. Before publication, we can access the data using the temporary URL below: https://datadryad.org/stash/share/if0mcxAb7wcCTZNgkI690VuxPDY-k5LweCABJJnTbf0.
Abbreviations
- ANP:
-
Atrial natriuretic peptide
- ACE-seq:
-
APOBEC-coupled epigenetic sequencing
- BMI:
-
Body mass index
- BS-seq:
-
Bisulfite sequencing
- CATIS:
-
China Antihypertensive Trial in Acute Ischemic Stroke
- CI:
-
Confidence interval
- CVD:
-
Cardiovascular disease
- DBP:
-
Diastolic blood pressure
- FDR:
-
False discovery rate
- GRCh37:
-
Genome Reference Consortium Human Build 37
- HDL:
-
High-density lipoprotein
- IS:
-
Ischemic stroke
- LDL:
-
Low-density lipoprotein
- mRS:
-
Modified Rankin Scale
- NIHSS:
-
National Institutes of Health Stroke Scale
- OR:
-
Odds ratio
- PMMS:
-
Prevention of Metabolic syndrome and multi-metabolic disorders Study
- SBP:
-
Systolic blood pressure
- TDG:
-
Thymine DNA glycosylase
- TSS:
-
Transcriptional start site
- wTPM:
-
Weighted truncated product method
- 5caC:
-
5-carboxylcytosine
- 5fC:
-
5-formylcytosine
- 5hmC:
-
5-hydroxymethyl cytosine
- 5mC:
-
5-methylcytosine
- β-GT:
-
β-glucosyltransferase
References
Levin ER, Gardner DG, Samson WK. Natriuretic peptides. N Engl J Med. 1998;339(5):321–8. PubMed PMID: 9682046.
Potter LR, Yoder AR, Flora DR, Antos LK, Dickey DM. Natriuretic peptides: their structures, receptors, physiologic functions and therapeutic applications. Handb Exp Pharmacol2009. pp. 341 – 66.
Rubattu S, Sciarretta S, Valenti V, Stanzione R, Volpe M. Natriuretic peptides: an update on bioactivity, potential therapeutic use, and implication in cardiovascular diseases. Am J Hypertens. 2008;21(7):733–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/ajh.2008.174. Epub 2008/05/10.
Song W, Wang H, Wu Q. Atrial natriuretic peptide in cardiovascular biology and disease (NPPA). Gene. 2015;569(1):1–6. 029. PubMed PMID: 26074089.
Berntsson J, Smith JG, Nilsson PM, Hedblad B, Melander O, Engström G. Pro-atrial natriuretic peptide and prediction of atrial fibrillation and stroke: the Malmö Preventive Project. Eur J Prev Cardiol. 2017;24(8):788–95. PubMed PMID: 28195503; PubMed Central PMCID: PMCQ1.
Macheret F, Heublein D, Costello-Boerrigter LC, Boerrigter G, McKie P, Bellavia D, et al. Human hypertension is characterized by a lack of activation of the antihypertensive cardiac hormones ANP and BNP. J Am Coll Cardiol. 2012;60(16):1558–65. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jacc.2012.05.049. Epub 2012/10/13.
Volpe M, Francia P, Tocci G, Rubattu S, Cangianiello S, Elena Rao MA, et al. Prediction of long-term survival in chronic heart failure by multiple biomarker assessment: a 15-year prospective follow-up study. Clin Cardiol. 2010;33(11):700–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/clc.20813. PubMed PMID: 21089115.
Sabatine MS, Morrow DA, de Lemos JA, Omland T, Sloan S, Jarolim P, et al. Evaluation of multiple biomarkers of cardiovascular stress for risk prediction and guiding medical therapy in patients with stable coronary disease. Circulation. 2012;125(2):233–40. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/CIRCULATIONAHA.111.063842. PubMed PMID: 22179538.
Barbato E, Bartunek J, Marchitti S, Mangiacapra F, Stanzione R, Delrue L, et al. NT-proANP circulating level is a prognostic marker in stable ischemic heart disease. Int J Cardiol. 2012;155(2):311–2. PubMed PMID: 22177588.
Katan M, Fluri F, Schuetz P, Morgenthaler NG, Zweifel C, Bingisser R, et al. Midregional pro-atrial natriuretic peptide and outcome in patients with acute ischemic stroke. J Am Coll Cardiol. 2010;56(13):1045–53. PubMed PMID: 20846604.
Murphy SP, Prescott MF, Camacho A, Iyer SR, Maisel AS, Felker GM, et al. Atrial Natriuretic Peptide and treatment with Sacubitril/Valsartan in heart failure with reduced ejection fraction. JACC Heart Fail. 2021;9(2):127–36. PubMed PMID: 33189632.
Suwa M, Seino Y, Nomachi Y, Matsuki S, Funahashi K. Multicenter prospective investigation on efficacy and safety of carperitide for acute heart failure in the ‘real world’ of therapy. Circ J. 2005;69(3):283–90. PubMed PMID: 15731532.
Nomura F, Kurobe N, Mori Y, Hikita A, Kawai M, Suwa M, et al. Multicenter prospective investigation on efficacy and safety of carperitide as a first-line drug for acute heart failure syndrome with preserved blood pressure: COMPASS: Carperitide effects observed through monitoring dyspnea in Acute Decompensated Heart failure study. Circ J. 2008;72(11):1777–86. PubMed PMID: 18832779.
Zhu H, Wu LF, Mo XB, Lu X, Tang H, Zhu XW, et al. Rheumatoid arthritis-associated DNA methylation sites in peripheral blood mononuclear cells. Ann Rheum Dis. 2019;78(1):36–42. https://doiorg.publicaciones.saludcastillayleon.es/10.1136/annrheumdis-2018-213970. Epub 2018/10/10.
Soriano-Tarraga C, Lazcano U, Giralt-Steinhauer E, Avellaneda-Gomez C, Ois A, Rodriguez-Campello A, et al. Identification of 20 novel loci associated with ischaemic stroke. Epigenome-wide association study. Epigenetics. 2020;15(9):988–97. PubMed PMID: 32202197; PubMed Central PMCID: PMCPMC7518691.
Cullell N, Soriano-Tarraga C, Gallego-Fabrega C, Carcel-Marquez J, Torres-Aguila NP, Muino E, et al. DNA methylation and ischemic stroke risk: an Epigenome-wide Association study. Thromb Haemost. 2022;122(10):1767–78. https://doiorg.publicaciones.saludcastillayleon.es/10.1055/s-0042-1749328. Epub 2022/06/20.
Zhao L, Chen X, Zhou S, Lin Z, Yu X, Huang Y. DNA methylation of AHCY may increase the risk of ischemic stroke. Bosn J Basic Med Sci. 2020;20(4):471–6. https://doiorg.publicaciones.saludcastillayleon.es/10.17305/bjbms.2020.4535. Epub 2020/02/06.
Santoro M, Siotto M, Germanotta M, Mastrorosa A, Papadopoulou D, Aprile I. Association Study of (5-HTTLPR) polymorphism and its promoter methylation with Rehabilitation Outcome in patients with Subacute Stroke. Genes. 2021;12(4):579. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/genes12040579. PubMed PMID: 33923526.
Santoro M, Siotto M, Germanotta M, Bray E, Mastrorosa A, Galli C, et al. BDNF rs6265 polymorphism and its methylation in patients with Stroke Undergoing Rehabilitation. Int J Mol Sci. 2020;21(22):8438. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms21228438. PubMed PMID: 33182716; PubMed Central PMCID: PMCQ1.
Ou W, Liu X, Shen Y, Li J, He L, Yuan Y, et al. Association of CVD candidate gene polymorphisms with ischemic stroke and cerebral hemorrhage in Chinese individuals. PLoS ONE. 2014;9(8):e105516. https://doiorg.publicaciones.saludcastillayleon.es/10.1371/journal.pone.0105516. PubMed PMID: 25144711.
Kim JH, Jang BH, Go HY, Park S, Shin YC, Kim SH, et al. Potential association between frequent nonsynonymous variant of NPPA and cardioembolic stroke. DNA Cell Biol. 2012;31(6):993–1000. https://doiorg.publicaciones.saludcastillayleon.es/10.1089/dna.2011.1463. Epub 2012/03/10.
Chen L, Li J, Zhang M, Zhang Q, Wu L, Lu Y, et al. Serum atrial natriuretic peptide, promoter methylation, and Cardiovascular Disease: a 10-year Follow-Up study in Chinese adults. Glob Heart. 2022;17(1):27. https://doiorg.publicaciones.saludcastillayleon.es/10.5334/gh.1116. PubMed PMID: 35586748.
Li J, Zhu J, Zhang Q, Chen L, Ma S, Lu Y, et al. NPPA promoter hypomethylation predicts central obesity development: a prospective longitudinal study in Chinese adults. Obes Facts. 2022;15(2):257–70. 10.1159/000521295. PubMed PMID: 34875662.
Li J, Zhu J, Ren L, Ma S, Shen B, Yu J, et al. Association between NPPA promoter methylation and hypertension: results from Gusu cohort and replication in an independent sample. Clin Epigenetics. 2020;12(1):133. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s13148-020-00927-0. PubMed PMID: 32883357.
Okano M, Bell DW, Haber DA, Li E. DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell. 1999;99(3):247–57. PubMed PMID: 10555141.
Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16(1):6–21. PubMed PMID: 11782440.
Ito S, D’Alessio AC, Taranova OV, Hong K, Sowers LC, Zhang Y. Role of Tet proteins in 5mC to 5hmC conversion, ES-cell self-renewal and inner cell mass specification. Nature. 2010;466(7310):1129–33. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nature09303. PubMed PMID: 20639862.
Branco MR, Ficz G, Reik W. Uncovering the role of 5-hydroxymethylcytosine in the epigenome. Nat Rev Genet. 2011;13(1):7–13. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nrg3080. PubMed PMID: 22083101.
Wu X, Zhang Y, TET-mediated active. DNA demethylation: mechanism, function and beyond. Nat Rev Genet. 2017;18(9):517 – 34. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nrg.2017.33. PubMed PMID: 28555658.
Ponnaluri VKC, Ehrlich KC, Zhang G, Lacey M, Johnston D, Pradhan S, et al. Association of 5-hydroxymethylation and 5-methylation of DNA cytosine with tissue-specific gene expression. Epigenetics. 2017;12(2):123–38. PubMed PMID: 27911668; PubMed Central PMCID: PMCQ2.
De Borre M, Branco MR. Oxidative Bisulfite Sequencing: An Experimental and Computational Protocol. Methods Mol Biol2021. pp. 333 – 48.
Peng H, Fan Y, Li J, Zheng X, Zhong C, Zhu Z, et al. Neurol Genet. 2022;8(3):e679. 0000000000000679. PubMed PMID: 35620136; PubMed Central PMCID: PMCPMC9128040. DNA Methylation of the Natriuretic Peptide System Genes and Ischemic Stroke: Gene-Based and Gene Set Analyses.
Morris-Blanco KC, Kim T, Lopez MS, Bertogliat MJ, Chelluboina B, Vemuganti R. Induction of DNA hydroxymethylation protects the Brain after Stroke. Stroke. 2019;50(9):2513–21. https://doiorg.publicaciones.saludcastillayleon.es/10.1161/STROKEAHA.119.025665. Epub 2019/07/23.
Miao Z, He Y, Xin N, Sun M, Chen L, Lin L, et al. Altering 5-hydroxymethylcytosine modification impacts ischemic brain injury. Hum Mol Genet. 2015;24(20):5855–66. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/hmg/ddv307. PubMed PMID: 26231219.
Xing DG, Zhang DY, Wang ZF, Ding DL, Wang J, Wang YJ. Correlations of ANP genetic polymorphisms and serum levels with ischemic stroke risk: a meta-analysis. Genet Test Mol Biomarkers. 2014;18(5):349–56. https://doiorg.publicaciones.saludcastillayleon.es/10.1089/gtmb.2013.0498. Epub 2014/03/25.
Pandey KN. Guanylyl cyclase / atrial natriuretic peptide receptor-A: role in the pathophysiology of cardiovascular regulation. Can J Physiol Pharmacol. 2011;89(8):557–73. https://doiorg.publicaciones.saludcastillayleon.es/10.1139/y11-054. PubMed PMID: 21815745.
Salo PP, Havulinna AS, Tukiainen T, Raitakari O, Lehtimaki T, Kahonen M, et al. Genome-Wide Association Study Implicates Atrial Natriuretic Peptide rather than B-Type Natriuretic Peptide in the regulation of blood pressure in the General Population. Circ Cardiovasc Genet. 2017;10(6):e001713. PubMed PMID: 29237677; PubMed Central PMCID: PMCPMC6072381.
Stanzione R, Cotugno M, Bianchi F, Marchitti S, Forte M, Volpe M, et al. Pathogenesis of ischemic stroke: role of epigenetic mechanisms. Genes (Basel). 2020;11(1):89. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/genes11010089. Epub 2020/01/17.
Hu CJ, Chen SD, Yang DI, Lin TN, Chen CM, Huang TH, et al. Promoter region methylation and reduced expression of thrombospondin-1 after oxygen-glucose deprivation in murine cerebral endothelial cells. J Cereb Blood Flow Metab. 2006;26(12):1519–26. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/sj.jcbfm.9600304. Epub 2006/03/30.
Hu Z, Zhong B, Tan J, Chen C, Lei Q, Zeng L. The emerging role of epigenetics in cerebral ischemia. Mol Neurobiol. 2017;54(3):1887–905. PubMed PMID: 26894397.
Kang HJ, Lee EH, Kim JW, Kim S-W, Shin I-S, Kim J-T, et al. Association of SLC6A4 methylation with long-term outcomes after stroke: focus on the interaction with suicidal ideation. Sci Rep. 2021;11(1):2710. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-021-81854-9. PubMed PMID: 33526821.
Cantone L, Tobaldini E, Favero C, Albetti B, Sacco RM, Torgano G, et al. Particulate Air Pollution, clock gene methylation, and stroke: effects on Stroke Severity and Disability. Int J Mol Sci. 2020;21(9):3090. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/ijms21093090. PubMed PMID: 32349365.
Zeng M, Zhen J, Zheng X, Qiu H, Xu X, Wu J et al. The Role of DNA Methylation in Ischemic Stroke: A Systematic Review. Front Neurol. 2020/11/17 ed2020. p. 566124.
Rustad SR, Papale LA, Alisch RS, editors. DNA Methylation and Hydroxymethylation and Behavior. Curr Top Behav Neurosci. 2019/08/09 ed2019. pp. 51–82.
Moore LD, Le T, Fan G. DNA methylation and its basic function. Neuropsychopharmacology. 2013;38(1):23–38. 10.1038. Epub 2012/07/12.
Estrada V, Téllez MJ, Moya J, Fernández-Durango R, Egido J, Fernández Cruz AF. High plasma levels of endothelin-1 and atrial natriuretic peptide in patients with acute ischemic stroke. Am J Hypertens. 1994;7(12):1085–9. PubMed PMID: 7702803.
Sato Y, Maruoka H, Honda Y, Hachiya N, Oizumi K. Plasma concentrations of atrial natriuretic peptide in cardioembolic stroke with atrial fibrillation. Kurume Med J. 1995;42(2):71–7. PubMed PMID: 7564167.
Wang TJ, Larson MG, Levy D, Benjamin EJ, Leip EP, Omland T, et al. Plasma natriuretic peptide levels and the risk of cardiovascular events and death. N Engl J Med. 2004;350(7):655–63. PubMed PMID: 14960742.
Ma X, Zhang M, Yan R, Wu H, Yang B, Miao Z. β2SP/TET2 complex regulates gene 5hmC modification after cerebral ischemia. J Cell Mol Med. 2021;25(24):11300–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/jcmm.17060. PubMed PMID: 34799994.
Breton CV, Marsit CJ, Faustman E, Nadeau K, Goodrich JM, Dolinoy DC, et al. Small-magnitude effect sizes in epigenetic end points are important in children’s Environmental Health studies: the Children’s Environmental Health and Disease Prevention Research Center’s epigenetics Working Group. Environ Health Perspect. 2017;125(4):511–26. PubMed PMID: 28362264; PubMed Central PMCID: PMCQ1.
Pu W, Wang C, Chen S, Zhao D, Zhou Y, Ma Y et al. Targeted bisulfite sequencing identified a panel of DNA methylation-based biomarkers for esophageal squamous cell carcinoma (ESCC). Clin Epigenetics. 2017/12/23 ed2017. p. 129.
Adams HP, Bendixen BH, Kappelle LJ, Biller J, Love BB, Gordon DL, et al. Classification of subtype of acute ischemic stroke. Definitions for use in a multicenter clinical trial. TOAST. Trial of Org 10172 in Acute Stroke Treatment. Stroke. 1993;24(1):35–41. PubMed PMID: 7678184; PubMed Central PMCID: PMCQ1.
Kwah LK, Diong J. National Institutes of Health Stroke Scale (NIHSS). J Physiother. 2014;60(1):61. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jphys.2013.12.012. PubMed PMID: 24856948.
Zaykin DV, Zhivotovsky LA, Westfall PH, Weir BS. Truncated product method for combining P-values. Genet Epidemiol. 2002;22(2):170–85. PubMed PMID: 11788962.
Desquilbet L, Mariotti F. Dose-response analyses using restricted cubic spline functions in public health research. Stat Med. 2010;29(9):1037-57. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/sim.3841. PubMed PMID: 20087875.
Acknowledgements
We gratefully acknowledge the cooperation and participation of our study. We especially thank the clinical staff at all participating hospitals for their support and contribution to this project. Without their contribution, this research would not have been possible. DNA methylation and hydroxymethylation were quantified by Genesky Biotechnologies.
Funding
This study was supported by the National Natural Science Foundation of China (NO. 82173596, 81903384, and 81872690), the Suzhou Municipal Science and Technology Bureau (NO. SYS2020091), Suzhou Key Technologies of Prevention and Control of Major Diseases and Infectious Diseases (NO. GWZX202001), and a Project of the Priority Academic Program Development of Jiangsu Higher Education Institutions.
Author information
Authors and Affiliations
Contributions
JY and LS performed the statistical data analysis and drafted the manuscript. HP and YH developed the concept of the study design and contributed to draft the manuscript. JZ, LC, WF, YL, XZ, JJ (Jiexiang), YJ, and MZ obtained the clinical data and critically reviewed the manuscript. JJ (Jun), YH, and LS contributed to the interpretation of the results. All authors contributed to draft the final versions of the manuscript. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval
The protocols of the CATIS were approved by the institutional review boards at Tulane University in the United States and Soochow University in China, as well as ethical committees at the participating hospitals. The protocols of the PMMS were approved by the Soochow University Ethics Committee. Written informed consent was obtained from all study participants or their immediate family members in both the CATIS and the PMMS.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Yao, J., Song, L., Jiang, J. et al. Associations of NPPA promoter true methylation and hydroxymethylation with ischemic stroke and its functional outcome. Epigenetics Commun. 4, 6 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s43682-024-00029-5
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s43682-024-00029-5