Skip to main content
Download PDF

HTRA1 methylation in peripheral blood as a potential marker for the preclinical detection of stroke: a case–control study and a prospective nested case–control study

Clinical Epigenetics volume 14, Article number: 191 (2022) Cite this article

A Correction to this article was published on 29 August 2023

This article has been updated

Abstract

Background

Stroke is the leading cause of mortality in China. DNA methylation has essential roles in multiple diseases, but its association with stroke was barely studied. We hereby explored the association between blood-based HTRA serine protease 1 (HTRA1) methylation and the risk of stroke.

Results

The association was discovered in a hospital-based case–control study (cases/controls = 190:190) and further validated in a prospective nested case–control study including 139 cases who developed stroke within 2 years after recruitment and 144 matched stroke-free controls. We observed stroke-related altered HTRA1 methylation and expression in both case–control study and prospective study. This blood-based HTRA1 methylation was associated with stroke independently from the known risk factors and mostly affected the older population. The prospective results further showed that the altered HTRA1 methylation was detectable 2 years before the clinical determination of stroke and became more robust with increased discriminatory power for stroke along with time when combined with other known stroke-related variables [onset time ≤ 1 year: area under the curve (AUC) = 0.76].

Conclusions

In our study, altered HTRA1 methylation was associated with stroke at clinical and preclinical stages and thus may provide a potential biomarker in the blood for the risk evaluation and preclinical detection of stroke.

Background

Stroke is a leading cause of mortality and a major contributor to neurological disability in low-income and middle-come countries [1, 2]. In 2019, stroke was responsible for 6.55 million deaths and 143.23 million disability-adjusted life years (DALYs) worldwide [2]. In China, there were 3.94 million new stroke cases and 2.19 million deaths with the highest DALYs of 45.9 million lost in 2019 [3].

Stroke is a complex disease attributed to genetic and environmental factors [4]. According to the genome-wide association studies, multiple genetic factors have only explained 5–10% of the risk for stroke [5, 6]. In recent years, attention has been paid to epigenetic factors and their role in the pathophysiological process of stroke [7]. DNA methylation is a crucial epigenetic modification involved in diabetes mellitus [8], atherosclerosis [9] and cardiovascular diseases [10]. DNA methylation can influence the structure of the genome and regulate the process of transcription, which mediates alterations in the expression of specific genes [11, 12]. Candidate approach studies have identified the association between stroke and altered blood-based DNA methylation in LINE-1, ABCB1 and CBS genes, but mainly in case–control studies with small sample sizes [13,14,15]. Our previous study has also revealed that aberrant methylation of ACTB in peripheral blood was associated with the risk of stroke [16].

The HTRA serine protease 1 (encoded by HTRA1 gene) is one of the human serine proteases that participate in diverse physiological processes including the maintenance of mitochondrial homeostasis, cell signaling and apoptosis [17]. Several studies have shown that heterozygous HTRA1 mutations are closely related to cerebral small vessel disease (SVD), small-vessel ischemic stroke and lacunar stroke [18,19,20,21,22,23]. Our previous study also showed that HTRA1 variants might confer the genetic susceptibility to stroke [24]. However, hereditary stroke attributable to HTRA1 gene constitutes only a small proportion of all strokes. So far, no data are available about the association between blood-based HTRA1 methylation and stroke, especially in the prospective studies.

In the present study, the association between peripheral blood-based HTRA1 methylation and the risk of stroke in the Chinese population was evaluated in an independent case–control study consisting of 190 stroke patients and 190 age- and sex-matched controls and was further verified in a prospective nested case–control study including 139 cases who developed stroke within 2 years after recruitment and 144 matched stroke-free controls from a cohort with a total of 11,151 subjects.

Results

Demographic characteristics of the subjects

The demographic and clinical characteristics of subjects in the case–control study and prospective nested case–control study are presented in Table 1. In the case–control study, the stroke cases had a higher proportion of drinking, higher prevalence of hypertension and diabetes, and higher levels of diastolic blood pressure (DBP), but lower levels of high-density lipoprotein cholesterol (HDL-C) than controls (all P < 0.05). No significant difference existed in age, sex, smoking status, systolic blood pressure (SBP), total cholesterol (TC), triglycerides (TG), low-density lipoprotein cholesterol (LDL-C), and fasting plasma glucose between the two groups.

Table 1 Demographic and clinical characteristics of subjects
Full size table

In the prospective nested case–control study, there were no significant differences in age, sex, SBP, DBP, smoking status, drinking status, history of hypertension, diabetes, levels of TC, TG, HDL-C, LDL-C and fasting plasma glucose between stroke cases and controls. The median time between blood draw and initial diagnosis of stroke in cases was 1.32 years.

Association between HTRA1 methylation and stroke in the case–control study and the prospective nested case–control study

To evaluate the association between HTRA1 methylation and the risk of stroke in the case–control study and the prospective nested case–control study, two amplicons, HTRA1-A and HTRA1-B, including 17 measurable CpG sites were semiquantitatively determined by Agena MALDI-TOF mass spectrometry.

In the case–control study (190 cases vs. 190 controls), HTRA1_A_CpG_6 showed significantly lower methylation levels in stroke cases than in the controls [median = 0.10 (interquartile range, IQR = 0.06–0.13) and 0.12 (IQR = 0.10–0.15) for cases and controls, respectively, OR per + 10% methylation = 0.54, 95% CI 0.34–0.86, P = 0.009 by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Table 2]. The other CpG sites in the HTRA1 amplicons showed no significant associations with stroke (Table 2).

Table 2 Methylation difference of HTRA1 between 190 stroke cases and 190 controls in the case–control study
Full size table

In the prospective nested case–control study, comparing 144 controls and all 139 cases who developed stroke within 2 years, the methylation levels of three measurable CpG loci (CpG_3.4, CpG_8.9, and CpG_13.14.15) in the HTRA1-B amplicon were significantly higher in cases than in controls (ORs per + 10% methylation from 1.33 to 1.42, P ≤ 0.034 for all by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Fig. 1A, Additional file 1: Table S1).

Fig. 1
figure 1

Onset-time-dependent association between HTRA1 methylation and stroke in the prospective nested case–control study. The analysis was performed for stroke cases with onset time A < 2 years, B ≤ 1.5 years, C ≤ 1.32 years, and D ≤ 1 year. The P values of all the 17 measureable CpG loci in the HTRA1-A and HTRA1-B amplicons were calculated by logistic regression adjusting for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, and all the P values were transformed by -log10. The dotted lines indicate the thresholds of different P values of 0.05. E and F The OR per + 10% methylation of HTRA1_A_CpG_6 and CpG_3.4, CpG_5, CpG_7, CpG_8.9, CpG_10, CpG_11.12 as well as CpG_13.14.15 in the HTRA1-B amplicon in stroke cases with onset time < 2 years, ≤ 1.5 years, ≤ 1.32 years and ≤ 1 year based on logistic regression analysis adjusting for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C

Full size image

Since the onset time of stroke may have an effect on the methylation levels of HTRA1, we further performed analyses comparing controls and stroke cases stratified by time of stroke onset in the prospective nested case–control study. Interestingly, the differences between the aforementioned three CpG loci became more pronounced in stroke cases with onset time ≤ 1.5 years (ORs per + 10% methylation from 1.43 to 1.58, P ≤ 0.016 for all by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Fig. 1B, Additional file 1: Table S2) and in cases with onset time ≤ 1.32 years (ORs per + 10% methylation from 1.43 to 1.59, P ≤ 0.025 for all by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Fig. 1C, Additional file 1: Table S3), especially in stroke cases with onset time ≤ 1 year (ORs per + 10% methylation from 1.67 to 1.88, P ≤ 0.022 for all by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Fig. 1D, Additional file 1: Table S4). Additionally, CpG_5, CpG_7, CpG_10 and CpG_11.12 in the HTRA1-B amplicon showed borderline associations with stroke in cases with shorter onset time (Fig. 1A–D, Additional file 1: Tables S1–S4). Moreover, stroke cases with onset time ≤ 1.5 years showed significantly lower HTRA1_A_CpG_6 methylation levels than controls (OR per + 10% methylation = 0.39, 95% CI 0.17–0.90, P = 0.027 by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Fig. 1B, Additional file 1: Table S2). Similarly, the differences became more pronounced in stroke cases with onset time ≤ 1.32 years (OR per + 10% methylation = 0.10, 95% CI 0.03–0.42, P = 0.002 by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Fig. 1C, Additional file 1: Table S3), and in cases with onset time ≤ 1 year (OR per + 10% methylation = 0.12, 95% CI 0.02–0.70, P = 0.018 by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Fig. 1D, Additional file 1: Table S4). HTRA1_A_CpG_7 and HTRA1_A_CpG_5 also showed significant association with stroke in cases with onset time ≤ 1.5 years and ≤ 1 year, respectively (P = 0.046 and 0.007 by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Fig. 1B and D Additional file 1: Table S2, Additional file 1: Table S4). Overall, the OR per + 10% methylation showed clear trends of decrease in HTRA1_A_CpG_6 in stroke cases with onset time < 2 years, ≤ 1.5 years, ≤ 1.32 years and ≤ 1 year (Fig. 1E). The OR per + 10% methylation showed clear trends of increase in CpG_3.4, CpG_5, CpG_7, CpG_8.9, CpG_10, CpG_11.12 and CpG_13.14.15 in the HTRA1-B amplicon in stroke cases with onset time < 2 years, ≤ 1.5 years, ≤ 1.32 years and ≤ 1 year (Fig. 1F).

Association between blood-based methylation of HTRA1 and stroke in the older population

Age was reported to play an essential role in the pattern of DNA methylation [25]; we therefore evaluated the association between HTRA1 methylation in the blood and stroke stratified by two major age groups (< 70 and ≥ 70 years). In the case–control study, stroke cases showed significantly lower methylation levels of HTRA1_A_CpG_3 and HTRA1_A_CpG_6 than controls among participants ≥ 70 years old (ORs per + 10% methylation = 0.63 and 0.29, P = 0.022 and 0.008, respectively, by logistic regression adjusted for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Table 3). However, no association between HTRA1 methylation and stroke was observed among subjects < 70 years old (Table 3).

Table 3 Methylation difference of HTRA1 between stroke cases and controls stratified by age in the case–control study
Full size table

In the prospective nested case–control study, the methylation levels of CpG_3.4, CpG_5, CpG_8.9 and CpG_13.14.15 in the HTRA1_B amplicon were significantly higher in stroke cases than in controls among participants ≥ 70 years old (ORs per + 10% methylation from 1.70 to 2.36, P ≤ 0.021 for all by logistic regression adjusting for age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C, Table 4). This association was not observed in subjects < 70 years old in the prospective nested case–control study as well (Table 4).

Table 4 Methylation difference of HTRA1 between stroke cases and controls stratified by age in the prospective nested case–control study
Full size table

Correlation between HTRA1 methylation and stroke-related characteristics

To explore the relationship between the blood-based HTRA1 methylation and stroke-related characteristics, the subjects including both stroke cases and controls with available data were interpreted. In the case–control study, weak correlations were observed between the methylation levels of a few CpG sites in the HTRA1 amplicons and hypertension, HDL-C, current smoking status, TC as well as TG levels (Additional file 1: Table S5). None of the 17 CpG loci showed any correlation with drinking, diabetes and levels of LDL-C (Additional file 1: Table S5).

In the prospective nested case–control study, the methylation level of HTRA1_A_CpG_2 was positively correlated with current drinking status (Spearman rho = 0.199, P = 0.001, Additional file 1: Table S6), while the methylation levels of CpG_3, CpG_4 and CpG_5 in HTRA1_A amplicon were inversely correlated with hypertension (Spearman rho = − 0.121, − 0.263 and − 0.198, respectively, Additional file 1: Table S6). No correlations were observed between individual CpG sites and smoking, diabetes, levels of TC, TG, HDL-C and LDL-C (Additional file 1: Table S6).

Expression of HTRA1 in peripheral blood leukocytes

To understand the correlation between the methylation and expression of HTRA1 gene, RNA was isolated from peripheral blood leukocytes of 96 stroke cases and 96 controls in the case–control study. The median of relative expression levels of HTRA1 in the leukocytes of stroke cases was 1.68-fold lower than in the controls (HTRA1 relative expression, median and IQR of cases: 0.56 (0.25–1.42), median and IQR of controls: 0.94 (0.36–2.52), P = 0.023, Fig. 2A).

Fig. 2
figure 2

The expression of HTRA1 in the peripheral blood leukocytes in controls and stroke cases. A The relative expression of HTRA1 between stroke cases and controls in the case–control study. B The relative expression of HTRA1 between stroke cases and controls in the prospective nested case–control study. The dots represent the relative expression levels of HTRA1 in each individual. The upper line and lower line represent upper quartile and lower quartile, and the medium line represents the median. The P values were calculated by Mann–Whitney U tests

Full size image

Next, we plan to measure the expression of HTRA1 in peripheral blood leukocytes in a prospective study. However, no RNA materials are available from the nested case–control study mentioned above as well as in the materials and methods. Thus, we analyzed blood RNA from 48 stroke cases who developed stroke within 1 year after the enrollment and 48 matched stroke-free controls from another prospective study in Jurong (Additional file 1: Table S7). The median of relative expression levels of HTRA1 in the leukocytes of stroke cases was 1.24-fold lower than in the controls (HTRA1 relative expression, median and IQR of cases: 0.76 (0.49–1.34), median and IQR of controls: 0.94 (0.64–1.46), Fig. 2B), although the P value did not reach significance (P = 0.054) probably due to the limited sample size.

HTRA1 methylation as a marker for the early detection of stroke

To estimate the potential clinical utility of HTRA1 methylation as a marker for the early detection of stroke, ROC curve analyses were performed based on the methylation data generated from the prospective nested case–control study using all the analyzed CpG sites and stroke-related variables including age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C by logistic regression. The combination of all the analyzed CpG sites and stroke-related variables showed better discriminatory power for differentiating stroke cases at any time of onset than stroke-related variables alone (Fig. 3A–D). For stroke cases with shorter onset time, the HTRA1 methylation and stroke-related variables showed increasing discriminatory power for differentiating stroke cases from healthy controls (area under the curve (AUC) = 0.66, 0.72 and 0.72 for stroke cases with onset time < 2 years, ≤ 1.5 years and ≤ 1.32 years, respectively, Fig. 3A–C). Notably, the combination of HTRA1 methylation and stroke-related variables had robust power to distinguish stroke cases with onset time ≤ 1 year from controls (AUC = 0.76, 95% CI 0.67–0.84, Fig. 3D).

Fig. 3
figure 3

Methylation level of HTRA1 in the peripheral blood DNA as a marker for the early detection of stroke. ROC curve analyses for the discriminatory power of HTRA1 methylation to distinguish stroke cases with onset time A < 2 years, B ≤ 1.5 years, C ≤ 1.32 years, and D ≤ 1 year from controls. The ROC analyses were performed by logistic regression using all the analyzed CpG sites and stroke-related variables including age, sex, smoking, drinking, hypertension, diabetes, TC, TG, HDL-C and LDL-C

Full size image

Discussion

In our study, the stroke-related altered HTRA1 methylation in blood cells was observed in both hospital-based case–control study and prospective nested case–control study. The prospective data further disclosed that the changed HTRA1 methylation was detectable 2 years preclinical of stroke and became more pronounced 1.5 years and 1 year before the clinical determination of stroke. These findings provided pronounced evidence that differential DNA methylation signatures in the peripheral blood could be a potential biomarker for the risk evaluation and preclinical detection of stroke.

HTRA1, a primarily secreted serine protease, degrades a number of substrates that involve in a variety of disease processes such as neurodegeneration [26], age-related macular degeneration [27], and arthritis [28]. The loss of HTRA1 function has been reported to cause cerebral autosomal recessive arteriopathy with subcortical infarcts and leukoencephalopathy (CARASIL), an inherited form of cerebral SVD characterized by an early-onset stroke in young adults [29]. To date, little is yet known about the epigenetic impact of HTRA1 on the initiation and progression of vascular disorders including stroke. This study was novel in reporting a significant association between aberrant blood-based HTRA1 and increased risk of stroke, and this association was more pronounced when the baseline was closer to the onset time of stroke. The prospective results suggested that the methylation level of HTRA1 in blood DNA might be a preclinical marker for stroke. Moreover, significantly decreased expression of HTRA1 was observed in the stroke cases, as well as in the preclinical stroke cases, indicating that the increased HTRA1 methylation in blood cells might be biologically meaningful. Future prospective multicenter studies are warranted to verify these findings.

Many studies disclosed that the DNA methylation in human peripheral leukocytes could vary with age [30, 31]. In our case–control study and prospective study, a significant association between methylation of HTRA1 and stroke was mainly observed in the older population, which was in accordance with the fact that the elderly men have a higher risk of stroke [32], and further indicated that HTRA1 methylation might be an age-dependent factor for the risk of stroke. Alternatively, the methylation pattern might be affected by lifestyle such as alcohol drinking [33] and other stroke-related factors like hypertension [34]. In our prospective nested case–control study, we observed that the methylation level of HTRA1_A_CpG_2 was positively correlated with current drinking status, while the methylation levels of CpG_3, CpG_4 and CpG_5 in the HTRA1_A amplicon were inversely correlated with hypertension. The specific mechanisms of DNA methylation and alcohol drinking as well as hypertension remained unclear. Further studies with larger sample sizes are needed to investigate the complex interactions between lifestyle, environmental factors, and DNA methylation throughout the progression of stroke.

Despite the prevalence of stroke, risk prediction and early detection remain challenging. Early detection of stroke is vital for providing timely treatment through thrombolysis, improving clinical outcomes and reducing resulting complications, and leading to better survival results and less disease burden [35, 36]. At present, magnetic resonance imaging (MRI) and computed tomography (CT) provide efficient tools for stroke diagnosis, but are only efficient at the onset of stroke. In this study, the blood-based HTRA1 methylation can objectively predict stroke a year in advance with an area under the curve of 0.76. We should admit that the HTRA1 methylation alone is not sufficient as a screening tool, the combination of multiple biomarkers including blood-based DNA methylation profiles and others as well as stroke-related environmental factors might be a solution for the screening and prediction of stroke in preclinic.

The long-term follow-up of stroke incidence and prospectively collected samples with careful uniform processing and storage in this study could exclude the possibility of observed differences between cases and controls owning to potential treatment or processing effects. Although with a 2-year follow-up in a cohort of 11,151 subjects, only 139 stroke cases were identified after the enrollment. This observed HTRA1 methylation-related risk of stroke was still based on a relatively small sample size. The association between HTRA1 methylation and stroke in the case–control study was not completely consistent with the prospective nested case–control study when each CpG was considered. This may also be attributed to the limited sample size. Further prospective studies with larger sample sizes are needed for validation and may also provide additional insights.

Conclusions

In conclusion, this study revealed a significant correlation between altered HTRA1 methylation and stroke in preclinic, thus suggesting the potential of blood-based DNA methylation for the risk evaluation and even prevention of stroke. Further validation studies are warranted to ensure the generalizability of our findings and evaluate the potential usage of HTRA1 methylation as a clinical biomarker for stroke.

Methods

Study population

Hospital-based case–control study

The case–control study included 190 stroke patients and 190 controls matched for age and sex. The patients (mean age of 64.99 ± 9.85 years) were recruited from the Affiliated Hospital of Xuzhou Medical University, Affiliated Jiangning Hospital of Nanjing Medical University and Jurong People's Hospital from January 2019 to August 2020. The clinical diagnosis of stroke was based on the clinical examination and confirmed by MRI or CT. The diagnostic information and clinical examination data were obtained from hospital medical records. The demographic information including age, sex, nationality, status of smoking and status of drinking and medical history of patients were collected using standardized questionnaires. The controls (mean age of 65.89 ± 10.12 years) were recruited from a cohort study conducted from November 2018 to January 2019 in Jurong City, Jiangsu Province, and those who remained stroke-free during the 2-year follow-up were selected. Blood samples were collected after 8 h from the last meal or after an overnight fast to measure the levels of fasting plasma glucose, TC, TG, HDL-C, and LDL-C. The individuals were considered as diabetes cases in the presence of fasting plasma glucose ≥ 7.0 mmol/L or a self-reported diabetes history according to the diagnostic criteria for diabetes suggested by the American Diabetes Association. Additionally, blood samples of all the cases and controls were processed following the same protocol and stored at − 80 °C after drawing.

Prospective nested case–control study

Subjects for this nested case–control study were selected from an ongoing population-based cohort study conducted in Jurong City, Jiangsu Province, which has been described in our previous study [16]. Briefly, a total of 11,151 subjects aged 18 years or older were recruited in this cohort from October to November of 2015, and all the subjects were reported to be stroke-free at baseline. Detailed information on demography and status of smoking and status of drinking as well as family history were collected via questionnaires, along with the collection of blood samples and anthropometric measurements at enrollment. The incidence of stroke was identified by records from the local hospitals and health service centers in the community. All the 139 subjects who developed a stroke within 2 years after the enrollment were included in this study as cases (mean age of 67.64 ± 9.51 years), and 150 participants who remained stroke-free during the 2.71-year follow-up time were selected as controls. After excluding 6 low-quality clotted blood samples, 144 controls (mean age of 67.77 ± 9.11 years) who well matched the stroke patients by age and sex were included.

DNA processing and bisulfite conversion

The Genomic DNA Extraction Kit (Rebece, Nanjing, China) was used to extract human peripheral blood cells’ genomic DNA as described previously [37]. Then the extracted DNA was further bisulfate-converted by EZ-96 DNA Methylation Gold Kit (Zymo Research, Orange County, USA) according to the manufacturer’s protocol. After bisulfite treatment, all non-methylated cytosine (C) bases were converted to uracil (U), whereas all methylated C bases were kept unchanged.

MALDI-TOF mass spectrometry

DNA methylation of HTRA1 in the peripheral blood was quantitatively measured using MALDI-TOF mass spectrometry (Agena Bioscience, San Diego, California, USA) as described by Yang et al. [38, 39]. Briefly, the bisulfite-treated genomic DNA was amplified by bisulfite-specific primers. The primers and sequences of polymerase chain reaction (PCR) amplicons (HTRA1-A and HTRA1-B) are presented in Additional file 1: Fig. S1. The HTRA1-A amplicon in the gene promoter region is located at 664 bp upstream of the nearby CpG island. The HTRA1-B amplicon in the first intron is located at 120 bp downstream of the nearby CpG island. Alternatively, these two amplicons are located at the so-called CpG island shores where differential methylation occurs frequently [40, 41]. There are no SNPs located at the primer regions or overlapped with any of the CpG sites in the amplicon. The PCR products were processed following shrimp alkaline phosphatase cleanup, T cleavage and Resin Clean steps. The final products were transferred to a SpectroCHIP G384 by a Nanodispenser RS1000 apparatus (Agena, USA). The chips were detected by MassARRAY spectrometry. The quantitative methylation levels of each CpG site or aggregate of multiple CpG sites were collected by EpiTYPER v1.2 software. For each batch of MassARRAY analyses, stroke case-and-control pairs were randomly arranged, and the samples were treated and analyzed in parallel in all the processes.

Quantitative real-time PCR

Total RNA was extracted from peripheral blood leucocytes using the QIAamp RNA Blood Mini Kit (Qiagen, Catalog no. 52304, Germany) and then reverse-transcribed into cDNA using the PrimeScript™ RT Reagent Kit (Takara, RR047A, Japan). Quantitative real-time PCR was performed to determine the expressions of HTRA1 gene and the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as an endogenous control using 2 × SYBR Green qPCR Master Mix (Bimake, B21202, USA). The relative expression of HTRA1 for each sample was normalized by GAPDH according to the 2−ΔΔct method.

Statistical analysis

All statistical analyses were conducted in SPSS version 25.0 (SPSS Inc., Chicago, USA). The methylation levels of each CpG site and gene expression data with non-Gaussian distribution were expressed as median (interquartile range). Differences between cases and controls were analyzed with Mann–Whitney U tests. The association between the level of HTRA1 methylation and stroke was assessed using logistic regression models with covariates-adjusted odds ratios (ORs) and 95% confidence intervals (CIs). The correlation between HTRA1 methylation and stroke-related characteristics was evaluated by Spearman’s rank correlation coefficients. Receiver operating characteristic (ROC) curve analysis was performed to assess the discriminatory power of methylation levels. A two‐tailed P < 0.05 was considered statistically significant.

Availability of data and materials

The datasets reported in the current study are available from the corresponding author upon reasonable request.

Change history

Abbreviations

AUC:

Area under the curve

CARASIL:

Cerebral autosomal recessive arteriopathy with subcortical infarcts and leukoencephalopathy

CIs:

Confidence intervals

CpG:

Cytosine–phosphate–guanine

CT:

Computed tomography

DBP:

Diastolic blood pressure

HDL-C:

High-density lipoprotein cholesterol

HTRA1:

HTRA serine protease 1

LDL-C:

Low-density lipoprotein cholesterol

MRI:

Magnetic resonance imaging

ORs:

Odds ratios

PCR:

Polymerase chain reaction

ROC:

Receiver operating characteristic

SBP:

Systolic blood pressure

SNP:

Single nucleotide polymorphism

SVD:

Small vessel disease

TC:

Total cholesterol

TG:

Triglycerides

TGF-β:

Transforming growth factor-β

References

  1. Kim J, Thayabaranathan T, Donnan GA, Howard G, Howard VJ, Rothwell PM, Feigin V, Norrving B, Owolabi M, Pandian J, Liu L, Cadilhac DA, Thrift AG. Global stroke statistics 2019. Int J Stroke. 2020;15(8):819–38.

    Article  PubMed  Google Scholar 

  2. Thayabaranathan T, Kim J, Cadilhac D, Thrift AG, Donnan G, Howard G, Howard VJ, Rothwell PM, Feigin VL, Norrving B, Owolabi M, Pandian JD, Liu L, Olaiya MT. Global stroke statistics 2022. Int J Stroke. 2022. https://doi.org/10.1177/17474930221123175.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Ma Q, Li R, Wang L, Yin P, Wang Y, Yan C, Ren Y, Qian Z, Vaughn MG, McMillin SE, Hay SI, Naghavi M, Cai M, Wang C, Zhang Z, Zhou M, Lin H, Yang Y. Temporal trend and attributable risk factors of stroke burden in China, 1990–2019: an analysis for the Global Burden of Disease Study 2019. Lancet Public Health. 2021;6(12):e897–906.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Benjamin EJ, Blaha MJ, Chiuve SE, Cushman M, Das SR, Deo R, de Ferranti SD, Floyd J, Fornage M, Gillespie C, Isasi CR, Jimenez MC, Jordan LC, Judd SE, Lackland D, Lichtman JH, Lisabeth L, Liu S, Longenecker CT, Mackey RH, Matsushita K, Mozaffarian D, Mussolino ME, Nasir K, Neumar RW, Palaniappan L, Pandey DK, Thiagarajan RR, Reeves MJ, Ritchey M, Rodriguez CJ, Roth GA, Rosamond WD, Sasson C, Towfighi A, Tsao CW, Turner MB, Virani SS, Voeks JH, Willey JZ, Wilkins JT, Wu JH, Alger HM, Wong SS, Muntner P, American Heart Association Statistics Committee, Stroke Statistics Subcommittee. Heart disease and stroke statistics-2017 update: a report from the American Heart Association. Circulation. 2017;135(10):e146–603.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Bevan S, Traylor M, Adib-Samii P, Malik R, Paul NL, Jackson C, Farrall M, Rothwell PM, Sudlow C, Dichgans M, Markus HS. Genetic heritability of ischemic stroke and the contribution of previously reported candidate gene and genomewide associations. Stroke. 2012;43(12):3161–7.

    Article  CAS  PubMed  Google Scholar 

  6. Malik R, Dichgans M. Challenges and opportunities in stroke genetics. Cardiovasc Res. 2018;114(9):1226–40.

    CAS  PubMed  Google Scholar 

  7. Felling RJ, Song H. Epigenetic mechanisms of neuroplasticity and the implications for stroke recovery. Exp Neurol. 2015;268:37–45.

    Article  CAS  PubMed  Google Scholar 

  8. Zhou Z, Sun B, Li X, Zhu C. DNA methylation landscapes in the pathogenesis of type 2 diabetes mellitus. Nutr Metab (Lond). 2018;15:47.

    Article  PubMed  Google Scholar 

  9. Zaina S, Heyn H, Carmona FJ, Varol N, Sayols S, Condom E, Ramirez-Ruz J, Gomez A, Goncalves I, Moran S, Esteller M. DNA methylation map of human atherosclerosis. Circ Cardiovasc Genet. 2014;7(5):692–700.

    Article  CAS  PubMed  Google Scholar 

  10. Liu G, Bin P, Wang T, Ren W, Zhong J, Liang J, Hu CA, Zeng Z, Yin Y. DNA methylation and the potential role of methyl-containing nutrients in cardiovascular diseases. Oxid Med Cell Longev. 2017;2017:1670815.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Nicoglou A, Merlin F. Epigenetics: a way to bridge the gap between biological fields. Stud Hist Philos Biol Biomed Sci. 2017;66:73–82.

    Article  PubMed  Google Scholar 

  12. Bird A. Perceptions of epigenetics. Nature. 2007;447(7143):396–8.

    Article  CAS  PubMed  Google Scholar 

  13. Yang J, Zhou JS, Zhao YX, Yang ZH, Zhao HD, Zhang YD, Zou JJ. ABCB1 hypomethylation is associated with decreased antiplatelet effects of clopidogrel in Chinese ischemic stroke patients. Pharmazie. 2015;70(2):97–102.

    CAS  PubMed  Google Scholar 

  14. Wang C, Xu G, Wen Q, Peng X, Chen H, Zhang J, Xu S, Zhang C, Zhang M, Ma J, Hui Z, Wu G, Ma M. CBS promoter hypermethylation increases the risk of hypertension and stroke. Clinics (Sao Paulo). 2019;74:e630.

    Article  PubMed  Google Scholar 

  15. Zhou S, Zhang Y, Wang L, Zhang Z, Cai B, Liu K, Zhang H, Dai M, Sun L, Xu X, Cai H, Liu X, Lu G, Xu G. CDKN2B methylation is associated with carotid artery calcification in ischemic stroke patients. J Transl Med. 2016;14(1):333.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Liu C, Yin Q, Li M, Fan Y, Shen C, Yang R. ACTB methylation in blood as a potential marker for the pre-clinical detection of stroke: a prospective nested case-control study. Front Neurosci. 2021;15:644943.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Zurawa-Janicka D, Wenta T, Jarzab M, Skorko-Glonek J, Glaza P, Gieldon A, Ciarkowski J, Lipinska B. Structural insights into the activation mechanisms of human HtrA serine proteases. Arch Biochem Biophys. 2017;621:6–23.

    Article  CAS  PubMed  Google Scholar 

  18. Wardlaw JM, Smith C, Dichgans M. Mechanisms of sporadic cerebral small vessel disease: insights from neuroimaging. Lancet Neurol. 2013;12(5):483–97.

    Article  PubMed  Google Scholar 

  19. Verdura E, Herve D, Scharrer E, Amador Mdel M, Guyant-Marechal L, Philippi A, Corlobe A, Bergametti F, Gazal S, Prieto-Morin C, Beaufort N, Le Bail B, Viakhireva I, Dichgans M, Chabriat H, Haffner C, Tournier-Lasserve E. Heterozygous HTRA1 mutations are associated with autosomal dominant cerebral small vessel disease. Brain. 2015;138(Pt 8):2347–58.

    Article  PubMed  Google Scholar 

  20. Nozaki H, Kato T, Nihonmatsu M, Saito Y, Mizuta I, Noda T, Koike R, Miyazaki K, Kaito M, Ito S, Makino M, Koyama A, Shiga A, Uemura M, Sekine Y, Murakami A, Moritani S, Hara K, Yokoseki A, Kuwano R, Endo N, Momotsu T, Yoshida M, Nishizawa M, Mizuno T, Onodera O. Distinct molecular mechanisms of HTRA1 mutants in manifesting heterozygotes with CARASIL. Neurology. 2016;86(21):1964–74.

    Article  CAS  PubMed  Google Scholar 

  21. Lee YC, Chung CP, Chao NC, Fuh JL, Chang FC, Soong BW, Liao YC. Characterization of heterozygous HTRA1 mutations in Taiwanese patients with cerebral small vessel disease. Stroke. 2018;49(7):1593–601.

    Article  CAS  PubMed  Google Scholar 

  22. Rannikmae K, Sivakumaran V, Millar H, Malik R, Anderson CD, Chong M, Dave T, Falcone GJ, Fernandez-Cadenas I, Jimenez-Conde J, Lindgren A, Montaner J, O’Donnell M, Pare G, Radmanesh F, Rost NS, Slowik A, Soderholm M, Traylor M, Pulit SL, Seshadri S, Worrall BB, Woo D, Markus HS, Mitchell BD, Dichgans M, Rosand J, Sudlow CLM, Stroke Genetics Network (SiGN), METASTROKE Collaboration, International Stroke Genetics Consortium (ISGC). COL4A2 is associated with lacunar ischemic stroke and deep ICH: meta-analyses among 21,500 cases and 40,600 controls. Neurology. 2017;89(17):1829–39.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Mishra A, Chauhan G, Violleau MH, Vojinovic D, Jian X, Bis JC, Li S, Saba Y, Grenier-Boley B, Yang Q, Bartz TM, Hofer E, Soumare A, Peng F, Duperron MG, Foglio M, Mosley TH, Schmidt R, Psaty BM, Launer LJ, Boerwinkle E, Zhu Y, Mazoyer B, Lathrop M, Bellenguez C, Van Duijn CM, Ikram MA, Schmidt H, Longstreth WT, Fornage M, Seshadri S, Joutel A, Tzourio C, Debette S. Association of variants in HTRA1 and NOTCH3 with MRI-defined extremes of cerebral small vessel disease in older subjects. Brain. 2019;142(4):1009–23.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Tian Y, Tang W, Yang S, Zhao Y, Chen Y, Zhao X, Liu C, Chen X, Shen C. HTRA1 variants and the interaction with smoking confer the genetic susceptibility to ischemic stroke. Int J Med Sci. 2021;18(8):1840–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Horvath S, Zhang Y, Langfelder P, Kahn RS, Boks MP, van Eijk K, van den Berg LH, Ophoff RA. Aging effects on DNA methylation modules in human brain and blood tissue. Genome Biol. 2012;13(10):R97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Tennstaedt A, Popsel S, Truebestein L, Hauske P, Brockmann A, Schmidt N, Irle I, Sacca B, Niemeyer CM, Brandt R, Ksiezak-Reding H, Tirniceriu AL, Egensperger R, Baldi A, Dehmelt L, Kaiser M, Huber R, Clausen T, Ehrmann M. Human high temperature requirement serine protease A1 (HTRA1) degrades tau protein aggregates. J Biol Chem. 2012;287(25):20931–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Yang Z, Camp NJ, Sun H, Tong Z, Gibbs D, Cameron DJ, Chen H, Zhao Y, Pearson E, Li X, Chien J, Dewan A, Harmon J, Bernstein PS, Shridhar V, Zabriskie NA, Hoh J, Howes K, Zhang K. A variant of the HTRA1 gene increases susceptibility to age-related macular degeneration. Science. 2006;314(5801):992–3.

    Article  CAS  PubMed  Google Scholar 

  28. Tiaden AN, Richards PJ. The emerging roles of HTRA1 in musculoskeletal disease. Am J Pathol. 2013;182(5):1482–8.

    Article  CAS  PubMed  Google Scholar 

  29. Hara K, Shiga A, Fukutake T, Nozaki H, Miyashita A, Yokoseki A, Kawata H, Koyama A, Arima K, Takahashi T, Ikeda M, Shiota H, Tamura M, Shimoe Y, Hirayama M, Arisato T, Yanagawa S, Tanaka A, Nakano I, Ikeda S, Yoshida Y, Yamamoto T, Ikeuchi T, Kuwano R, Nishizawa M, Tsuji S, Onodera O. Association of HTRA1 mutations and familial ischemic cerebral small-vessel disease. N Engl J Med. 2009;360(17):1729–39.

    Article  CAS  PubMed  Google Scholar 

  30. Jung SE, Shin KJ, Lee HY. DNA methylation-based age prediction from various tissues and body fluids. BMB Rep. 2017;50(11):546–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Chen X, Zhao C, Zhao Z, Wang H, Fang Z. Specific glioma prognostic subtype distinctions based on DNA methylation patterns. Front Genet. 2019;10:786.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Chao TF, Lip GY, Liu CJ, Tuan TC, Chen SJ, Wang KL, Lin YJ, Chang SL, Lo LW, Hu YF, Chen TJ, Chiang CE, Chen SA. Validation of a modified CHA2DS2-VASc score for stroke risk stratification in Asian patients with atrial fibrillation: a nationwide cohort study. Stroke. 2016;47(10):2462–9.

    Article  CAS  PubMed  Google Scholar 

  33. Liu C, Marioni RE, Hedman AK, Pfeiffer L, Tsai PC, Reynolds LM, Just AC, Duan Q, Boer CG, Tanaka T, Elks CE, Aslibekyan S, Brody JA, Kuhnel B, Herder C, Almli LM, Zhi D, Wang Y, Huan T, Yao C, Mendelson MM, Joehanes R, Liang L, Love SA, Guan W, Shah S, McRae AF, Kretschmer A, Prokisch H, Strauch K, Peters A, Visscher PM, Wray NR, Guo X, Wiggins KL, Smith AK, Binder EB, Ressler KJ, Irvin MR, Absher DM, Hernandez D, Ferrucci L, Bandinelli S, Lohman K, Ding J, Trevisi L, Gustafsson S, Sandling JH, Stolk L, Uitterlinden AG, Yet I, Castillo-Fernandez JE, Spector TD, Schwartz JD, Vokonas P, Lind L, Li Y, Fornage M, Arnett DK, Wareham NJ, Sotoodehnia N, Ong KK, van Meurs JBJ, Conneely KN, Baccarelli AA, Deary IJ, Bell JT, North KE, Liu Y, Waldenberger M, London SJ, Ingelsson E, Levy D. A DNA methylation biomarker of alcohol consumption. Mol Psychiatry. 2018;23(2):422–33.

    Article  CAS  PubMed  Google Scholar 

  34. Han L, Liu Y, Duan S, Perry B, Li W, He Y. DNA methylation and hypertension: emerging evidence and challenges. Brief Funct Genom. 2016;15(6):460–9.

    CAS  Google Scholar 

  35. Ducas RA, Labos C, Allen D, Golian M, Jeyaraman M, Lys J, Mann A, Copstein L, Vokey S, Rabbani R, Zarychanski R, Abou-Setta AM, Menkis AH. Association of pre-hospital ECG administration with clinical outcomes in ST-segment myocardial infarction: a systematic review and meta-analysis. Can J Cardiol. 2016;32(12):1531–41.

    Article  PubMed  Google Scholar 

  36. Saver JL. Time is brain-quantified. Stroke. 2006;37(1):263–6.

    Article  PubMed  Google Scholar 

  37. Yin Q, Yang X, Li L, Xu T, Zhou W, Gu W, Ma F, Yang R. The association between breast cancer and blood-based methylation of S100P and HYAL2 in the Chinese population. Front Genet. 2020;11:977.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Yang R, Pfutze K, Zucknick M, Sutter C, Wappenschmidt B, Marme F, Qu B, Cuk K, Engel C, Schott S, Schneeweiss A, Brenner H, Claus R, Plass C, Bugert P, Hoth M, Sohn C, Schmutzler R, Bartram CR, Burwinkel B. DNA methylation array analyses identified breast cancer-associated HYAL2 methylation in peripheral blood. Int J Cancer. 2015;136(8):1845–55.

    Article  CAS  PubMed  Google Scholar 

  39. Yang R, Stocker S, Schott S, Heil J, Marme F, Cuk K, Chen B, Golatta M, Zhou Y, Sutter C, Wappenschmidt B, Schmutzler R, Bugert P, Qu B, Bartram CR, Sohn C, Schneeweiss A, Burwinkel B. The association between breast cancer and S100P methylation in peripheral blood by multicenter case-control studies. Carcinogenesis. 2017;38(3):312–20.

    Article  CAS  PubMed  Google Scholar 

  40. Doi A, Park IH, Wen B, Murakami P, Aryee MJ, Irizarry R, Herb B, Ladd-Acosta C, Rho J, Loewer S, Miller J, Schlaeger T, Daley GQ, Feinberg AP. Differential methylation of tissue- and cancer-specific CpG island shores distinguishes human induced pluripotent stem cells, embryonic stem cells and fibroblasts. Nat Genet. 2009;41(12):1350–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, Ji H, Potash J, Sabunciyan S, Feinberg AP. The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009;41(2):178–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank Yu Liu, Min Liu and Haifeng Xu for their contribution to the study. We would also like to express our gratitude to all of the participants of the contributing studies to each of these datasets.

Funding

This work was funded by the National Natural Science Foundation of China (Grant No. 81573232) and by the research funding from Nanjing Medical University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Author notes

  1. Chunlan Liu and Mengxia Li contributed equally to this work

Authors and Affiliations

  1. Department of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, 101 Longmian Avenue, Jiangning, Nanjing, 211166, China

    Chunlan Liu, Mengxia Li, Qiming Yin, Chong Shen & Rongxi Yang

  2. Division of Clinical Epidemiology, Affiliated Geriatric Hospital of Nanjing Medical University, Nanjing, 211166, China

    Yao Fan

Authors
  1. Chunlan Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Mengxia Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Qiming Yin
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yao Fan
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Chong Shen
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Rongxi Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

CL and ML analyzed the data and drafted the article. QY and YF collected the data. CS and RY performed the conception or design of the work, data interpretation and critical revision of the article. All authors reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Chong Shen or Rongxi Yang.

Ethics declarations

Ethics approval and consent to participate

The studies were approved by the Ethics Committee of Nanjing Medical University. All the recruited participants were Chinese Han and provided written informed consent.

Consent for publication

All authors give consent for the publication of the manuscript.

Competing interests

Rongxi Yang, Chong Shen, Chunlan Liu and Qiming Yin are inventors of a provisional patent application relating to the subject matter of this manuscript and therefore declare potential competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Table S1. Methylation difference of HTRA1 between 144 controls and 139 stroke cases with onset time < 2 years in the prospective nested casecontrol study. Table S2. Methylation difference of HTRA1 between 144 controls and 91 stroke cases with onset time ≤ 1.5 years in the prospective nested case-control study. Table S3. Methylation difference of HTRA1 between 144 controls and 67 stroke cases with onset time ≤ 1.32 years in the prospective nested case-control study. Table S4. Methylation difference of HTRA1 between 144 controls and 35 stroke cases with onset time ≤ 1 year in the prospective nested case-control study. Table S5. Correlation between HTRA1 methylation and stroke-related characteristics in the case-control study. Table S6. Correlation between HTRA1 methylation and stroke-related characteristics in the prospective nested case-control study. Table S7. Description of the subjects for RNA analysis. Fig. S1 Amplicon and primer design for MassARRAY methylation analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, C., Li, M., Yin, Q. et al. HTRA1 methylation in peripheral blood as a potential marker for the preclinical detection of stroke: a case–control study and a prospective nested case–control study. Clin Epigenet 14, 191 (2022). https://doi.org/10.1186/s13148-022-01418-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-022-01418-0

Keywords

Clinical Epigenetics

ISSN: 1868-7083

Contact us