ORIGINAL RESEARCH article

Front. Genet., 11 April 2022

Sec. Epigenomics and Epigenetics

Volume 13 - 2022 | https://doi.org/10.3389/fgene.2022.844141

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

    HS

    Hongwei Sun 1

    JX

    Jia Xu 2

    BH

    Bifeng Hu 1

    YL

    Yue Liu 1

    YZ

    Yun Zhai 1

    YS

    Yanyan Sun 1

    HS

    Hongwei Sun 1

    FL

    Fang Li 1

    JW

    Jiamin Wang 2

    AF

    Anqi Feng 2

    YT

    Ying Tang 1*

    JZ

    Jingbo Zhao 2*

  • 1. Department of Neurology, The First Affiliated Hospital of Harbin Medical University, Harbin, China

  • 2. Department of Epidemiology, School of Public Health, Harbin Medical University, Harbin, China

Article metrics

View details

9

Citations

2,8k

Views

1,4k

Downloads

Abstract

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

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

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

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

Introduction

Stroke is a devastating disease due to its high morbidity, disability, and recurrence. The prevalence of stroke is increasing over time, and the affected population is becoming younger (Katan and Luft, 2018). Ischemic stroke (IS) is the dominant subtype of stroke with a proportion more than 80% (Ajoolabady et al., 2021). Many of the environmental and genetic risk factors are associated with IS, and the genetic risk is believed to be in the order of 37.9%. However, genetic variants associated with IS found to date only account for 5–10% of that genetic risk (Bevan et al., 2012; Krupinski et al., 2017), which suggests that more associated heritable risk factors have not yet been discovered. Epigenetic modifications is one of these possible heritable changes (Mahjoubin-Tehran et al., 2021).

Epigenetics refer to chemical modifications of DNA structure without affecting the DNA sequence that may provide a link between environment and gene expression (Majnik and Lane, 2014). Generally, epigenetic modifications of gene expression occur by three main forms: DNA methylation, histone modification and microRNA expression (Schiano et al., 2020; Liu et al., 2020). As a major type of epigenetic process, DNA methylation mainly occurs at the cytosine of a cytosine-phosphate-guanine (CpG) dinucleotide forming 5-methylcytosine, which can inhibit gene expression through transcriptional silencing (Bird, 2007). DNA methylation is a crucial epigenetic mechanism involved in normal and pathological cellular processes (Krupinski et al., 2017).

The current epigenome-wide studies revealed that DNA methylation plays a vital part in the pathogenesis and recurrence of IS (Soriano-Tárraga et al., 2020; Davis Armstrong et al., 2018; Gómez-Úriz et al., 2015; Shen et al., 2019), which have mainly been conducted in American and European populations (Soriano-Tárraga et al., 2020; Davis Armstrong et al., 2018; Gómez-Úriz et al., 2015). An epigenetic study in a Chinese population revealed that hypomethylation of the MTRNR2L8 gene is associated with IS (Shen et al., 2019). It is undeniable that this study is significant; however, it focused on the large-artery atherosclerosis stroke, a subgroup of IS, which cannot fully demonstrate the role of DNA methylation in IS. Thus far, epigenome-wide alterations of IS have not been systematically investigated in a Chinese population.

In our research, we designed a two-stage case-control study to perform integrated analysis of genome-wide DNA methylation profiles to identify novel candidate genes and pathways for IS in a Chinese population.

Materials and Methods

Ethical Approval

We performed this study according to the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments. The study procedure was approved by the Ethics Committee of the First Affiliated Hospital of Harbin Medical University. All participants voluntarily gave written informed consent.

Study Population

In this study, a two-stage case-control DNA methylation study was designed, including discovery and replication analyses. The sample consisted of 192 IS patients and 192 age-(±3 years) and sex-matched controls from the First Affiliated Hospital of Harbin Medical University. Figure 1 demonstrates the study design and selection criteria for cases and controls. The control group was from healthy people who underwent physical examination or patients hospitalized from the First Affiliated Hospital of Harbin Medical University at the same period.

FIGURE 1

FIGURE 1

Two-stage case-control study design and selection criteria for cases and controls. IS, ischemic stroke; DMPs, DNA methylation positions; GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; CpG, cytosine phosphate guanine.

Data and Sample Collection

Data for subjects included questionnaire, laboratory investigations, and clinical characteristics. We collected data using a questionnaire related to demographic information (such as gender, age, and marital status), lifestyle habits (such as tobacco smoking and alcohol drinking), and past medical history (hypertension and type 2 diabetes mellitus) of the participants (Qin et al., 2019).

Laboratory investigations were done at the First Affiliated Hospital of Harbin Medical University. We collected blood samples in tubes with clot activator and gel separator from all participants in the morning after fasting and measured using an auto biochemical analyzer (Biobase BK-600; Shandong, China). The levels of total triglyceride (TG), total cholesterol (TC), high-density lipoprotein (HDL), and low-density lipoprotein (LDL) was collected from all participants.

Clinical characteristics included Trial of Org 10172 in Acute Stroke Treatment (TOAST) type and stroke location. Stratification analysis by TOAST criteria was performed for large-artery atherosclerosis (LAA), cardioembolism (CE), small vessel disease (SVD), other determined cause (ODC), and undetermined (UND) aetiologies. Stroke location was categorized as anterior, posterior, and both anterior and posterior circulation. Anterior circulation was further subdivided as left and right hemisphere based on computerized tomography or magnetic resonance imaging scans.

Two milliliters of blood samples were collected in EDTA tubes from all participants and stored at −80°C until DNA extraction.

Discovery Stage: Genome-Wide DNA Methylation Profiles

Genomic DNA was extracted from the leukocytes in the peripheral blood using an DNA Methylation kit (Zymo, Irvine, CA), following the protocol (Genesky Biotechnologies Inc., Shanghai, China). Genome-wide DNA methylation was assessed using Infinium Human Methylation 850K BeadChip (Illumina Inc., San Diego, CA, United States)). More than 853,000 CpG sites were included in each chip. CpG sites containing documented single-nucleotide polymorphisms and mapping to X and Y chromosomes were removed to avoid potential confounding by single-nucleotide polymorphisms and gender.

Discovery Stage: Bioinformatics Analysis of DNA Methylation Profiles

To screen functional DNA methylation positions (DMPs), DMPs located in promoter regions (upstream or 5′untranslated region, 5′UTR), with |methylation difference|>0.1 and p < 0.05 were selected. We performed Gene ontology (GO, http://www.geneontology.org) and the Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.kegg.jp/) pathway enrichment analysis to clarify the function and biological pathways of differential methylation loci. Through GO enrichment analysis, the differentially methylated genes were classified according to cellular component, molecular function and biological process. In GO and KEGG enrichment analysis, p < 0.05 was significantly enriched by differential methylation loci-related genes.

Replication Stage: MethylTargetTM Assay

Considering the p-values and biological functions of the genes, we selected 8 candidate genes from GO analytics top-ranked pathway including 10 CpG sites to validate 850 K BeadChips results using MethylTargetTM assay, which was a multi-targeted CpG methylation analysis method based on next-generation sequencing (Wan et al., 2021). As the 10 CpG sites in 850 K BeadChips results are not top-ranked loci and FDR p values of the 10 CpG sites are not statistically significant, we replicated 10 CpGs in a small replication sample (20 cases:20 matched controls, 8 samples from DNA methylation chip analysis and 32 samples from validation analysis). Moreover, we further qualified the methylation level of promoter regions in 8 candidate genes, containing 17 CpG islands and 308 CpG sites, using a large replication sample (188 cases:188 individually matched controls). The locations of the 17 CpG islands in 8 genes are shown in Figure 2. The primers were carefully designed (Additional File S1).

FIGURE 2

FIGURE 2

Diagram of the structure of the 8 genes and 17 CpG islands. The line represents the gene, and double arrows represent CpG islands. CpG, cytosine phosphate guanine; TSS, transcription start sites; Chr, chromosome.

Validation of the mRNA Expression of Differentially Methylated Genes

The dataset GSE22255 was acquired from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) and was utilized as the validation sample. Gene expression profiling was performed in peripheral blood of 20 IS patients and 20 sex- and age-matched controls. From this dataset, the downloaded data format was MINIML. Box plots of expression were drawn by the R software package ggplot2, and principal component analysis graphs were drawn by the R software package ggord. The Wilcoxon rank-sum test was used to detect whether the mRNA expression of differentially methylated genes was different between stroke patients and controls.

Statistical Analysis

Quantitative variables were reported as the mean and standard deviation, while qualitative variables were expressed as frequencies and percentages (%). The normality of quantitative variables was assessed using the Shapiro–Wilk test.

First, we compared the differences of baseline demographic characteristics, lifestyle habits, medical history, and lipid levels between IS cases and controls. The differences of the categorical variables were measured with a chi-square test. We performed a paired t test for the normal distributed continuous variables. Non-parametric test was used for the skewed distributed continuous variables.

Second, the methylation status was compared between IS cases and marched controls using conditional logistic regression analysis. Smoking, drinking, medical history of diabetes mellitus and hypertension, and lipid levels were added as adjusted factors. The association of the methylation level of genes was assessed using Pearson’s correlation. All subjects were classified into four groups (Q1, Q2, Q3, and Q4) based on the quartile value of methylation level. Q1 group having the lowest methylation level is regarded as the reference group, and all other groups (Q2 having 25–50% of the values, Q3 having 50–75%, and Q4 having 75–100%) were compared with Q1 group (Qin et al., 2019).

The methylation level of each CpG site was measured as the percentage of the methylated cytosines over total tested cytosines. The average methylation level of all CpG sites was worked out as the methylation level of the DNA segments and genes. Results were considered statistically significant when the p values were less than 0.05. All statistical analyses were performed with the SAS version 9.4 for Windows (SAS Institute, Inc., Cary, NC).

Results

Characteristics of Participants

The distribution of characteristics for the participants are presented in Table 1. Compared with controls, IS patients in validation stage had a higher prevalence of smoking, drinking and hypertension (p < 0.0001), and higher level of triglyceride (p = 0.0288).

TABLE 1

DNA methylation chip analysisPDNA methylation validation analysisP
IS (n = 4)Control (n = 4)IS (n = 188)Control (n = 188)
Age65.2397 ± 2.656165.9904 ± 2.92580.717162.2155 ± 10.666062.1466 ± 10.57390.9499
Female2 (50%)2 (50%)1.00044 (23.40%)44 (23.40%)1.000
Smoking (%)1 (25%)1 (25%)1.000106 (56.38%)41 (21.81%)<0.0001*
Drinking (%)1 (25%)1 (25%)1.00091 (48.40%)39 (20.74%)<0.0001*
Hypertension (%)0 (0%)0 (0%)1.000112 (59.57%)69 (36.70%)<0.0001*
Diabetes (%)0 (0%)0 (0%)1.00035 (18.62)26 (13.83%)0.2080
Total cholesterol (mmol/L)5.4900 ± 0.86765.1325 ± 0.27490.46205.4598 ± 4.56054.8362 ± 1.05030.0691
Triglyceride (mmol/L)1.7350 ± 0.89161.4225 ± 1.01740.66041.9978 ± 1.37281.7139 ± 1.12170.0288*
High-density lipoprotein (mmol/L)1.1650 ± 0.28831.5425 ± 0.59160.29501.2963 ± 0.95801.2990 ± 0.40340.9715
Low-density lipoprotein (mmol/L)3.5500 ± 0.44972.9450 ± 0.36920.08283.2304 ± 0.79742.9960 ± 1.83160.1089

Characteristics of participants.

*Statistically significant difference (p < 0.05). IS, ischemic stroke.

Discovery Stage: Epigenome-Wide Association Study of IS

In the discovery analysis, 2656 DMPs exhibited differences (p < 0.05, |methylation difference|>0.1) in DNA methylation between the two groups. The results are shown in another unpublished article (Xu J, et al.). After screening, a total of 462 functional DMPs (located in the promoter region, |methylation difference| > 0.1 and p < 0.05) corresponding to 373 genes were selected (Additional File S2). Among them, the majority of the DMPs were in CpG island (26.4%) and N-shore (28.6%). In addition, 87.7% (405 of 462) CpG cites corresponding to 327 genes were found to be hypomethylated in IS cases compared to controls.

Discovery Stage: Bioinformatics Analysis of DNA Methylation Profiles

We screened and finally included 12 KEGG pathways (p < 0.05, Figure 3). A total of 315 significant GO terms (p < 0.05) were selected. The top 10 GO items were shown in Table 2. Among them, the “calcium-dependent cell-cell adhesion via plasma membrane cell adhesion molecules” item showed the most significant fold enrichment. There are 8 genes in this item including CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB6/PCDHB9, and the 8 genes were selected for replication.

FIGURE 3

FIGURE 3

Top 12 enriched KEGG pathways of differential methylation loci-related genes. The dot size represents the number of genes in the pathway, and the dot color represents the enrichment p value. KEGG, Kyoto Encyclopedia of Genes and Genomes; GABA, gamma-aminobutyric acid.

TABLE 2

IDDescriptionGenesP
BP
 GO:0016339calcium-dependent cell-cell adhesion via plasma membrane cell adhesion moleculesPCDHB14/PCDHB11/PCDHB3/PCDHB16/PCDHB9/CDH2/PCDHB10/PCDHB6p < 0.0001
 GO:0007156homophilic cell adhesion via plasma membrane adhesion moleculesDSC3/PCDHGA4/PCDHA6/PCDHB14/PCDHGA11/PCDHB11/PCDHB3/PCDHGA1/PCDHA1/PCDHB16/PCDHB9/CDH2/PCDHB7/PCDHB10/PCDHB6p < 0.0001
 GO:0098742cell-cell adhesion via plasma-membrane adhesion moleculesDSC3/PCDHGA4/PCDHA6/CBLN1/PCDHB14/PCDHGA11/PCDHB11/PCDHB3/PCDHGA1/PCDHA1/PCDHB16/PCDHB9/CDH2/PCDHB7/PCDHB10/CD164/PCDHB6p < 0.0001
 GO:0007416synapse assemblyCBLN1/NRG1/PCDHB14/THBS2/PCDHB11/PCDHB3/PCDHB16/PCDHB9/CDH2/NRXN2/PCDHB7/PCDHB10/PCDHB6p < 0.0001
 GO:0050808synapse organizationANK3/CBLN1/NRG1/PCDHB14/THBS2/CTNND2/CACNB2/PCDHB11/PCDHB3/PCDHB16/PCDHB9/CDH2/NRXN2/PCDHB7/PCDHB10/PCDHB6p < 0.0001
 GO:0048863stem cell differentiationPSMA8/NRG1/KIT/EDN3/TP73/MYOCD/TAL1/PSMD1/PSMB2/SFRP1/SEMA3C/MSX1/YTHDF2p = 0.0008
CC
 GO:0045211postsynaptic membranePDLIM4/ANK3/GRIK4/GRIK1/CBLN1/TANC1/CHRM2/ABI1/GABRA5/SSPN/GABRG1/GRM7/NETO1/GABRB3p < 0.0001
 GO:0098794postsynapseADD2/PDLIM4/ANK3/GRIK4/GRIK1/CBLN1/CALB1/TANC1/CHRM2/ABI1/GABRA5/SSPN/GABRG1/GRM7/NETO1/CDH2/KCNN2/GABRB3p = 0.0004
 GO:0097060synaptic membranePDLIM4/ANK3/GRIK4/GRIK1/CBLN1/TANC1/CHRM2/ABI1/GABRA5/SSPN/GABRG1/GRM7/NETO1/GABRB3p = 0.0006
MF
 GO:0005338nucleotide-sugar transmembrane transporter activitySLC35D2/SLC35C1/SLC35A1p = 0.0005

Top 10 enriched GO pathways of differential methylation loci-related genes.

GO, gene ontology; BP, biological process; CC, cellular component; MF, molecular function.

Replication Stage

In the replication stage, 5 of 10 CpG sites in 8 candidate genes from 850K BeadChip results were verified to be statistically different (p < 0.05). The trend of the methylation level of the rest 5 CpG sites with no significant differences is consistent with the 850k BeadChip results (Additional File S3), suggesting that the 850k results are reliable. The genes screened by GO analysis are worthy of further verification.

In 188 paired samples, 203 of 308 CpG sites from 17 CpG islands in the promoter regions of 8 candidate genes were significantly different (p < 0.05) between the IS and control group (Additional File S4). Consistent with the 850k BeadChip results, the methylation level of all 203 CpG sites was lower in the IS group than in the control group. 13 of 17 CpG islands were differentially methylated between IS patients and controls (p < 0.05, Additional File S5). 7 genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9) showed differences (p < 0.05) between the two groups even if adjusting for smoking, drinking, history of hypertension and diabetes, and blood lipid levels, and adjust FDR p values of 6 genes (CDH2/PCDHB10/PCDHB11/PCDHB16/PCDHB3/PCDHB9) were statistically significant (p < 0.05, Table 3). The methylation levels of the 7 genes were highly correlated in both the case group and the control group (p < 0.05, Additional File S6). Besides, we further split the samples by sex (male or female). There were 7 differentially methylated genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9) in the male subgroup and 4 (CDH2/PCDHB10/PCDHB11/PCDHB9) in the female subgroup, which is shown in Table 4.

TABLE 3

TargetGroup differencep valueAdjust p valueAdjust FDR p value
CDH2−0.0127p < 0.0001*p < 0.0001*p = 0.0001*
PCDHB10−0.0228p < 0.0001*p < 0.0001*p = 0.0002*
PCDHB11−0.0232p < 0.0001*p = 0.0001*p = 0.0006*
PCDHB14−0.0193p = 0.0002*p = 0.0498*p = 0.0940
PCDHB16−0.0162p = 0.0117*p = 0.0148*p = 0.0315*
PCDHB3−0.0158p = 0.0016*p = 0.0089*p = 0.0217*
PCDHB6−0.0087p = 0.3376p = 0.3399p = 0.5268
PCDHB9−0.0251p < 0.0001*p = 0.0002*p = 0.0008*

Validation of 8 candidate genes in the ischemic stroke group and control group.

Adjusted factors: smoking, drinking, previous medical history of hypertension and diabetes mellitus, and plasma lipid levels (total triglyceride, total cholesterol, high-density lipoprotein, and low-density lipoprotein). *Statistically significant difference (p < 0.05).

TABLE 4

TargetGroup differencep valueAdjust p value
Male (case:control = 144:144)
CDH2−0.0118p < 0.0001*p = 0.0006*
PCDHB10−0.0211p = 0.0001*p < 0.0001*
PCDHB11−0.0216p = 0.0005*p = 0.0003*
PCDHB14−0.0201p = 0.0002*p = 0.0139*
PCDHB16−0.0153p = 0.0327*p = 0.0070*
PCDHB3−0.0139p = 0.0081*p = 0.0038*
PCDHB6−0.0038p = 0.7111p = 0.3113
PCDHB9−0.0230p = 0.0003*p = 0.0003*
Female (case:control = 44:44)
CDH2−0.0158p = 0.0014*p = 0.0299*
PCDHB10−0.0283p = 0.0103*p = 0.0270*
PCDHB11−0.0283p = 0.0129*p = 0.0160*
PCDHB14−0.0164p = 0.2108p = 0.0562
PCDHB16−0.0192p = 0.1810p = 0.2095
PCDHB3−0.0220p = 0.0875p = 0.0907
PCDHB6−0.0247p = 0.2093p = 0.3205
PCDHB9−0.0317p = 0.0124*p = 0.0195*

Stratified analysis by gender.

Adjusted factors: smoking, drinking, previous medical history of hypertension and diabetes mellitus, and plasma lipid levels (total triglyceride, total cholesterol, high-density lipoprotein, and low-density lipoprotein). *Statistically significant difference (p < 0.05).

Stratified Analysis by IS Subtypes and Location

We split the case samples by IS subtypes (TOAST classification), stratified by LAA (n = 67), CE (n = 13), SVD (n = 64), ODC (n = 5), and UND (n = 39). However, no analysis was conducted in the CE and ODC groups because of the small sample size. We also did not perform analysis in the UND group as it is hard to identify the cause of IS. There were 8 differentially methylated genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB6/PCDHB9) in the LAA subgroup and 1 (CDH2) in the SVD subgroup after adjusting for smoking, drinking, history of hypertension and diabetes, and blood lipid levels (p < 0.05, Table 5). Besides, the case samples were seperated according to the location of the lesion (anterior circulation-left hemisphere, n = 72; anterior circulation-right hemisphere, n = 63; posterior circulation, n = 48; and both anterior and posterior circulation, n = 5). We also did not perform analysis in the “both anterior and posterior circulation” group because of the small sample size. Four genes (CDH2/PCDHB10/PCDHB11/PCDHB9) in the anterior circulation-left hemisphere subgroup, 8 genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB6/PCDHB9) in the anterior circulation-right hemisphere subgroup and 4 genes (CDH2/PCDHB10/PCDHB11/PCDHB9) in the posterior circulation subgroup showed differences after adjusting for smoking, drinking, history of hypertension and diabetes, and blood lipid levels (p < 0.05, Table 6).

TABLE 5

TargetGroup differencep valueAdjust p value
LAA (case:control = 67:67)
CDH2−0.0138p = 0.0016*p = 0.0059*
PCDHB10−0.0249p = 0.0015*p = 0.0028*
PCDHB11−0.0264p = 0.0035*p = 0.0034*
PCDHB14−0.0273p = 0.0016*p = 0.0070*
PCDHB16−0.0204p = 0.0597p = 0.0043*
PCDHB3−0.0182p = 0.0247*p = 0.0035*
PCDHB6−0.0146p = 0.3849p = 0.0434*
PCDHB9−0.0292p = 0.0014*p = 0.0015*
SVD (case:control = 64:64)
CDH2−0.0101p = 0.0022*p = 0.0343*
PCDHB10−0.0190p = 0.0322*p = 0.1069
PCDHB11−0.0195p = 0.0502p = 0.1195
PCDHB14−0.0131p = 0.1816p = 0.3445
PCDHB16−0.0151p = 0.1758p = 0.3626
PCDHB3−0.0093p = 0.3182p = 0.5970
PCDHB6−0.0031p = 0.8273p = 0.9188
PCDHB9−0.0117p = 0.2385p = 0.4120

Stratified analysis by ischemic stroke subtypes.

Adjusted factors: smoking, drinking, previous medical history of hypertension and diabetes mellitus, and plasma lipid levels (total triglyceride, total cholesterol, high-density lipoprotein, and low-density lipoprotein). *Statistically significant difference (p < 0.05). LAA, large-artery atherosclerosis; SVD, small vessel disease.

TABLE 6

TargetGroup differencep valueAdjust p value
Anterior circulation-left hemisphere (case:control = 72:72)
CDH2−0.0133p = 0.0005*p = 0.0015*
PCDHB10−0.0200p = 0.0201*p = 0.0278*
PCDHB11−0.0191p = 0.0419*p = 0.0400*
PCDHB14−0.0155p = 0.0997p = 0.2873
PCDHB16−0.0111p = 0.2982p = 0.1031
PCDHB3−0.0093p = 0.3007p = 0.1201
PCDHB60.0034p = 0.8128p = 0.3847
PCDHB9−0.0237p = 0.0214*p = 0.0144*
Anterior circulation-right hemisphere (case:control = 63:63)
CDH2−0.0103p = 0.0005*p = 0.0058*
PCDHB10−0.0249p = 0.0026*p = 0.0037*
PCDHB11−0.0257p = 0.0075*p = 0.0060*
PCDHB14−0.0219p = 0.0067*p = 0.0470*
PCDHB16−0.0204p = 0.0644p = 0.0243*
PCDHB3−0.0224p = 0.0051*p = 0.0076*
PCDHB6−0.0248p = 0.1181p = 0.0276*
PCDHB9−0.0257p = 0.0052*p = 0.0090*
Posterior circulation (case:control = 48:48)
CDH2−0.0174p = 0.0022*p = 0.0363*
PCDHB10−0.0243p = 0.0066*p = 0.0113*
PCDHB11−0.0265p = 0.0064*p = 0.0196*
PCDHB14−0.0187p = 0.0557p = 0.0717
PCDHB16−0.0190p = 0.1315p = 0.1532
PCDHB3−0.0172p = 0.0661p = 0.1380
PCDHB6−0.0112p = 0.5626p = 0.9945
PCDHB9−0.0268p = 0.0075*p = 0.0205*

Stratified analysis by ischemic stroke location.

Adjusted factors: smoking, drinking, previous medical history of hypertension and diabetes mellitus, and plasma lipid levels (total triglyceride, total cholesterol, high-density lipoprotein, and low-density lipoprotein). *Statistically significant difference (p < 0.05).

The Association Between DNA Methylation in 7 Differentially Methylated Genes and IS

The methylation level of each gene was grouped into quartiles. The DNA methylation level of the 7 differentially methylated genes in Q4 groups was associated with higher risk of IS than Q1 groups (p < 0.05), which are listed in Table 7.

TABLE 7

TargetOR95%CIp valueAdjust p value
CDH2
 Q2 vs. Q18.8443.776–20.710p < 0.0001*p < 0.0001*
 Q3 vs. Q16.0922.593–14.309p < 0.0001*p < 0.0001*
 Q4 vs. Q124.6159.518–63.660p < 0.0001*p < 0.0001*
PCDHB10
 Q2 vs. Q12.6341.363–5.089p = 0.0040*p = 0.0004*
 Q3 vs. Q13.6271.775–7.413p = 0.0004*p = 0.0001*
 Q4 vs. Q17.7973.707–16.399p < 0.0001*p < 0.0001*
PCDHB11
 Q2 vs. Q12.6941.371–5.296p = 0.0040*p = 0.0016*
 Q3 vs. Q14.1031.982–8.497p = 0.0001*p = 0.0003*
 Q4 vs. Q16.3462.963–13.592p < 0.0001*p < 0.0001*
PCDHB14
 Q2 vs. Q11.4020.774–2.540p = 0.2653p = 0.3518
 Q3 vs. Q14.0352.051–7.937p < 0.0001*p = 0.0055*
 Q4 vs. Q14.0482.053–7.983p < 0.0001*p = 0.0028*
PCDHB16
 Q2 vs. Q12.8621.537–5.327p = 0.0009*p = 0.0007*
 Q3 vs. Q11.3190.671–2.595p = 0.4221p = 0.2949
 Q4 vs. Q13.3551.743–6.456p = 0.0003*p = 0.0002*
PCDHB3
 Q2 vs. Q11.3900.764–2.529p = 0.2813p = 0.0897
 Q3 vs. Q12.0321.125–3.671p = 0.0187*p = 0.0055*
 Q4 vs. Q12.9651.561–5.630p = 0.0009*p = 0.0016*
PCDHB9
 Q2 vs. Q12.4021.221–4.725p = 0.0111*p = 0.0008*
 Q3 vs. Q13.5001.799–6.807p = 0.0002*p = 0.0010*
 Q4 vs. Q15.3362.656–10.720p < 0.0001*p < 0.0001*

Associations between DNA methylation in 7 differentially methylated genes and risk of ischemic stroke.

Adjusted factors: smoking, drinking, previous medical history of hypertension and diabetes mellitus, and plasma lipid levels (total triglyceride, total cholesterol, high-density lipoprotein, and low-density lipoprotein). *Statistically significant difference (p < 0.05). OR odds ratio; CI, confidence interval.

Biomarker Potential of the 7 Differentially Methylated Genes for IS

We tested the biomarker potential of the 7 differentially methylated genes using receiver operating characteristic (ROC) analysis. The sensitivity, specificity, accuracy and area under the curve (AUC) in CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9 genes was determined based on methylation values of replication stage (Table 8, Figure 4A). To enhance the diagnostic ability of these differentially methylated genes in IS, we tested the performance of the polygenic methylation model. The sensitivity, specificity, accuracy and AUC of polygenic methylation model was shown in Table 8 and Figure 4B using ROC analysis. The polygenic methylation model performed better than each of the individual genes, especially with higher specificity (0.8723), sensitivity (0.883), accuracy (0.8777), and AUC (0.9384).

TABLE 8

TargetAdjust AUCAdjust SensitivityAdjust SpecificityAdjust Accuracy
CDH20.90670.8830.78720.8351
PCDHB100.90150.86170.79790.8298
PCDHB110.89240.82980.81910.8245
PCDHB140.86790.81910.81910.8191
PCDHB160.8740.82450.81910.8218
PCDHB30.87930.82450.82450.8245
PCDHB90.8910.86170.78720.8245
Polygenic model0.93840.8830.87230.8777

ROC analysis of 7 differentially methylated genes and polygenic methylation model.

Adjusted factors: smoking, drinking, previous medical history of hypertension and diabetes mellitus, and plasma lipid levels (total triglyceride, total cholesterol, high-density lipoprotein, and low-density lipoprotein). AUC area under the curve; ROC receiver operating characteristic.

FIGURE 4

FIGURE 4

ROC for 7 differentially methylated genes and polygenic methylation model. (A) ROC curve measuring sensitivity (Y axis) and 1-specificity (X axis) of methylation of CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9 genes as biomarkers for discriminating IS from controls. (B) ROC curve of the polygenic methylation model. AUC, area under the curve; ROC, receiver operating characteristic; CI, confidence interval.

Validation of the mRNA Expression of Differentially Methylated Genes

In the replication stage, there were 7 differentially hypomethylated genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9) in the IS group. Generally, a decrease in DNA methylation leads to abnormal gene transcription, resulting in the upregulation of gene expression. We validated the expression levels of these 7 genes in the GSE22255 dataset. The expression of two genes (PCDHB9 and PCDHB11) increased, while the methylation level of the same 2 genes decreased in IS patients compared to controls (p < 0.05, Figure 5). The expression levels of the other three genes (CDH2/PCDHB14/PCDHB16) also increased in the IS group, but there were no significant differences between the two groups (p > 0.05, Additional Files S7–11).

FIGURE 5

FIGURE 5

Validation of the mRNA expression of PCDHB9 and PCDHB11 in GSE22255. (A): Methylation level of the PCDHB9 gene in the control and case (IS) groups; (B): Methylation level of the PCDHB11 gene in the control and case (IS) groups; (C): mRNA expression of the PCDHB9 gene in the control and case (IS) groups; (D): mRNA expression of the PCDHB11 gene in the control and case (IS) groups. *Statistically significant difference (p < 0.05). ****Statistically significant difference (p < 0.0001).

Discussion

Certainly, IS is a heterogeneous disease with a genetic predisposition (Acosta et al., 2021). Emerging evidence has revealed the importance of epigenetic regulations in IS, particularly DNA methylation (Gómez-Úriz et al., 2015; Davis Armstrong et al., 2018; Shen et al., 2019; Soriano-Tárraga et al., 2020). Our study presented genome-wide alterations in DNA methylation in IS patients and controls by identifying a total of 462 functional DMPs corresponding to 373 annotated genes and revealed that hypomethylated sites were eightfold more numerous than DNA hypermethylation sites in IS cases, demonstrating that the hypomethylated modification was predominant in IS.

For DNA methylation profiling, GO analysis showed that these genes form interconnected networks involved in the calcium-dependent cell-cell adhesion via plasma membrane cell adhesion molecules, homophilic cell adhesion via plasma membrane adhesion molecules, cell-cell adhesion via plasma-membrane adhesion molecules and other networks. Among them, the most enriched term in GO analysis was “calcium-dependent cell-cell adhesion via plasma membrane cell adhesion molecules” item including 8 candidate genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB6/PCDHB9), which were related to the cadherins.

Cadherins are a superfamily of calcium-dependent adhesion molecules mainly involved in tissue and embryonic cell development (Schaarschuch and Hertel, 2018). They can be classified into several subfamilies such as classic cadherins (CDHs), protocadherins (PCDHs), desmosomal cadherins and so on (Pancho et al., 2020). The PCDH gene family comprises three gene clusters (PCDHA, PCDHB, and PCDHG, respectively) (Peek et al., 2017; Canzio and Maniatis, 2019). PCDHs are highly expressed in the central nervous system in neurons, astrocytes, pericytes and brain microvascular endothelial cells, and mediate various developmental processes, including synaptic maintenance, neuronal survival, and spatial patterning of axons and dendrites (Gabbert et al., 2020; Mancini et al., 2020; Ing-Esteves et al., 2018). PCDHs have been reported to be relevant to neurodegenerative diseases, epilepsy, schizophrenia and mood disorders (Schaarschuch and Hertel, 2018; Langfelder et al., 2016; Li et al., 2017; Pederick et al., 2018; Lencz et al., 2021; Flaherty and Maniatis, 2020). However, the association between PCDHs and IS is rarely reported. Ana et al. found altered PCDHGA3 gene expression was strongly associated with reduced stroke volume (Ortega et al., 2016). Some studies revealed that cadherins are involved in the regulation of angiogenesis and inflammation in endothelial cells and lack of PCDHs may lead to the destruction of the blood-brain barrier (Gabbert et al., 2020; Nanes et al., 2012), which are the mechanisms of IS. However, no direct evidence of the association between PCDHB genes and IS has been shown in previous studies. The semi-stochastic expression of PCDH genes is regulated by DNA methylation (Hirayama and Yagi, 2017). In our results, the methylation level of PCDHB (PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9) genes was lower in IS cases compared with controls, which suggested the methylation status of PCDHB genes may be involved in the pathogenesis of IS and be potential diagnostic biomarkers for IS.

Like the PCDHs, classic CDHs have been implicated in neurulation, brain development, and regulation of synaptic function (Polanco et al., 2021; de Agustín-Durán et al., 2021; Sanes and Zipursky, 2020). They are also involved in the formation of atherosclerotic plaques (He et al., 2017), which is an important clinical feature of IS. One of the most important members of CDHs is the neuronal cadherin (N-cadherin, CDH2) (Schaarschuch and Hertel, 2018; Martinez-Garay, 2020). So far, there has been no direct evidence to show the relationship between CDH2 gene and stroke. A Study about myocardial infarction revealed overexpression of CDH2, CDH12, PCDH17, and PCDH18 in myocardial infarction vascular smooth muscle cells compared with controls (Derda et al., 2018). In this study we found decreased methylation level of the promoter region in CDH2 gene. Generally, hypomethylation in promotor is believed to upregulate gene transcription, which is consistent with the previous study of myocardial infarction in a way. This is the first time we have found the correlation between the methylation of CDH2 gene and IS.

There are also some interesting findings in our study. First, when we split the samples by sex, we found fewer differentially methylated genes (4 genes) in females, which we thought might be because the sample size of women is much smaller than that of men. What’s more, compared with SVD-control groups, there are more differentially methylated genes in LAA-control groups. We thought it might hint that the LAA subtype may be more susceptible to epigenetic regulation, and further studies are needed. Lastly, there existed more differentially methylated genes in the right hemisphere-control groups than the other two paired groups, which need further studies in a large sample size.

The peripheral blood is a good choice for epigenetic research of IS, as it is easy to obtain with minor invasion (Qin et al., 2019). Moreover, IS is a disease related to the vasculature and interrupting blood supply to the brain. Liu et al. found that stroke patients had lower methylation level of the ACTB gene in blood (Liu et al., 2021), which revealed that peripheral blood could identify the methylation aberrations associated with stroke.

There were certain limitations in our study. First, gene expression was verified in a public dataset, but not explored in our blood samples. Our study aimed to describe the global DNA methylation patterns in IS and to explore potential diagnostic biomarkers for IS in a Chinese population. Although gene and protein expression was not the main purpose of our research, it will be carried out in follow-up studies. In addition, we couldn’t detect the methylation level of more CpG sites due to the limited funds. However, they will be the directions for our future study.

Conclusion

The present study demonstrated the changes in genome-wide DNA methylation between IS cases and controls and identified 7 novel DNA methylation genes (CDH2/PCDHB10/PCDHB11/PCDHB14/PCDHB16/PCDHB3/PCDHB9) related to IS in the replication stage. These data may provide preliminary evidence for further exploring the role of DNA methylation in IS.

Statements

Data availability statement

The datasets presented in this article are not readily available because of ethical restrictions. Requests to access the datasets should be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Harbin Medical University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

Study design: JZ and YT. Blood collection: HS, JX, BH and YL. Data check: YZ, YS, HS, FL, JW and AF. Manuscript drafting: HS. Manuscript edition: HS. Administrative and material support: JZ and YT. Study supervision: JZ and YT. All authors read and approved the final manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant number 81771508); Scientific research innovation fund of the First Affiliated Hospital of Harbin Medical University (Grant number 2020B05).

Acknowledgments

We thank all participants for their cooperation throughout the whole study. We acknowledged the technical support from the Shanghai Genesky Biotechnology Company (Shanghai, China).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.844141/full#supplementary-material

Abbreviations

CDHs, cadherins; CE, cardioembolism; CpG, cytosine phosphate guanine; DMPs, DNA methylation positions; GO, gene ontology; HDL, high-density lipoprotein; IS, ischemic stroke; KEGG, Kyoto Encyclopedia of Genes and Genomes; LAA, large-artery atherosclerosis; LDL, low-density lipoprotein; ODC, other determined cause; PCDHs, protocadherins; ROC, receiver operating characteristic; SVD, small vessel disease; TC, total cholesterol; TG, total triglyceride; TOAST, Trial of ORG 10172 in Acute Stroke Treatment; UND, undetermined; 5′UTR, 5′untranslated region; AUC, area under the curve.

References

  • 1

    AcostaJ. N.SzejkoN.FalconeG. J. (2021). Mendelian Randomization in Stroke: A Powerful Approach to Causal Inference and Drug Target Validation. Front. Genet.12, 683082. 10.3389/fgene.2021.683082

  • 2

    AjoolabadyA.WangS.KroemerG.PenningerJ. M.UverskyV. N.PraticoD. (2021). Targeting Autophagy in Ischemic Stroke: From Molecular Mechanisms to Clinical Therapeutics. Pharmacol. Ther.225, 107848. 10.1016/j.pharmthera.2021.107848

  • 3

    BevanS.TraylorM.Adib-SamiiP.MalikR.PaulN. L. M.JacksonC.et al (2012). Genetic Heritability of Ischemic Stroke and the Contribution of Previously Reported Candidate Gene and Genomewide Associations. Stroke43, 31613167. 10.1161/STROKEAHA.112.665760

  • 4

    BirdA. (2007). Perceptions of Epigenetics. Nature447, 396398. 10.1038/nature05913

  • 5

    CanzioD.ManiatisT. (2019). The Generation of a Protocadherin Cell-Surface Recognition Code for Neural Circuit Assembly. Curr. Opin. Neurobiol.59, 213220. 10.1016/j.conb.2019.10.001

  • 6

    Davis ArmstrongN. M.ChenW. M.BrewerM. S.WilliamsS. R.SaleM. M.WorrallB. B.et al (2018). Epigenome-Wide Analyses Identify Two Novel Associations with Recurrent Stroke in the Vitamin Intervention for Stroke Prevention Clinical Trial. Front. Genet.9, 358. 10.3389/fgene.2018.00358

  • 7

    de Agustín-DuránD.Mateos-WhiteI.Fabra-BeserJ.Gil-SanzC. (2021). Stick Around: Cell–Cell Adhesion Molecules during Neocortical Development. Cells10, 118. 10.3390/cells10010118

  • 8

    DerdaA. A.WooC. C.WongsurawatT.RichardsM.LeeC. N.KofidisT.et al (2018). Gene Expression Profile Analysis of Aortic Vascular Smooth Muscle Cells Reveals Upregulation of Cadherin Genes in Myocardial Infarction Patients. Physiol. Genomics50, 648657. 10.1152/physiolgenomics.00042.2017

  • 9

    FlahertyE.ManiatisT. (2020). The Role of Clustered Protocadherins in Neurodevelopment and Neuropsychiatric Diseases. Curr. Opin. Genet. Develop.65, 144150. 10.1016/j.gde.2020.05.041

  • 10

    GabbertL.DillingC.MeybohmP.BurekM. (2020). Deletion of Protocadherin Gamma C3 Induces Phenotypic and Functional Changes in Brain Microvascular Endothelial Cells In Vitro. Front. Pharmacol.11, 590144. 10.3389/fphar.2020.590144

  • 11

    Gómez-ÚrizA. M.MilagroF. I.MansegoM. L.CorderoP.AbeteI.De ArceA.et al (2015). Obesity and Ischemic Stroke Modulate the Methylation Levels of KCNQ1 in white Blood Cells. Hum. Mol. Genet.24, 14321440. 10.1093/hmg/ddu559

  • 12

    HeX.LiD. R.CuiC.WenL. J. (2017). Clinical Significance of Serum MCP-1 and VE-Cadherin Levels in Patients with Acute Cerebral Infarction. Eur. Rev. Med. Pharmacol. Sci.21, 804808.

  • 13

    HirayamaT.YagiT. (2017). Regulation of Clustered Protocadherin Genes in Individual Neurons. Semin. Cel Develop. Biol.69, 122130. 10.1016/j.semcdb.2017.05.026

  • 14

    Ing-EstevesS.KostadinovD.MarochaJ.SingA. D.JosephK. S.LaboulayeM. A.et al (2018). Combinatorial Effects of Alpha- and Gamma-Protocadherins on Neuronal Survival and Dendritic Self-Avoidance. J. Neurosci.38, 2713. 10.1523/JNEUROSCI.3035-17.2018

  • 15

    KatanM.LuftA. (2018). Global Burden of Stroke. Semin. Neurol.38, 208211. 10.1055/s-0038-1649503

  • 16

    KrupinskiJ.CarreraC.MuiñoE.TorresN.Al-BaradieR.CullellN.et al (2017). DNA Methylation in Stroke. Update of Latest Advances. Comput. Struct. Biotechnol. J.16, 15. 10.1016/j.csbj.2017.12.001

  • 17

    LangfelderP.CantleJ. P.ChatzopoulouD.WangN.GaoF.Al-RamahiI.et al (2016). Integrated Genomics and Proteomics Define Huntingtin CAG Length–dependent Networks in Mice. Nat. Neurosci.19, 623633. 10.1038/nn.4256

  • 18

    LenczT.YuJ.KhanR. R.FlahertyE.CarmiS.LamM.et al (2021). Novel Ultra-rare Exonic Variants Identified in a Founder Population Implicate Cadherins in Schizophrenia. Neuron109, 14651478. 10.1016/j.neuron.2021.03.004

  • 19

    LiY.ChenZ.GaoY.PanG.ZhengH.ZhangY.et al (2017). Synaptic Adhesion Molecule Pcdh-Γc5 Mediates Synaptic Dysfunction in Alzheimer's Disease. J. Neurosci.37, 9259. 10.1523/jneurosci.1051-17.2017

  • 20

    LiuC.YinQ.LiM.FanY.ShenC.YangR. (2021). ACTB Methylation in Blood as a Potential Marker for the Pre-clinical Detection of Stroke: A Prospective Nested Case-Control Study. Front. Neurosci.15, 644943. 10.3389/fnins.2021.644943

  • 21

    LiuX.FanB.ChoppM.ZhangZ. (2020). Epigenetic Mechanisms Underlying Adult Post Stroke Neurogenesis. Int. J. Mol. Sci.21, 6179. 10.3390/ijms21176179

  • 22

    Mahjoubin-TehranM.RezaeiS.JesmaniA.BirangN.MorshediK.KhanbabaeiH.et al (2021). New Epigenetic Players in Stroke Pathogenesis: From Non-coding RNAs to Exosomal Non-coding RNAs. Biomed. Pharmacother.140, 111753. 10.1016/j.biopha.2021.111753

  • 23

    MajnikA. V.LaneR. H. (2014). Epigenetics: where Environment, Society and Genetics Meet. Epigenomics6, 14. 10.2217/epi.13.83

  • 24

    ManciniM.BassaniS.PassafaroM. (2020). Right Place at the Right Time: How Changes in Protocadherins Affect Synaptic Connections Contributing to the Etiology of Neurodevelopmental Disorders. Cells9, 2711. 10.3390/cells9122711

  • 25

    Martinez-GarayI. (2020). Molecular Mechanisms of Cadherin Function during Cortical Migration. Front. Cel. Dev. Biol.8, 588152. 10.3389/fcell.2020.588152

  • 26

    NanesB. A.Chiasson-MacKenzieC.LoweryA. M.IshiyamaN.FaundezV.IkuraM.et al (2012). p120-catenin Binding Masks an Endocytic Signal Conserved in Classical Cadherins. J. Cel Biol199, 365380. 10.1083/jcb.201205029

  • 27

    OrtegaA.Gil-CayuelaC.TarazónE.García-ManzanaresM.MonteroJ. A.CincaJ.et al (2016). New Cell Adhesion Molecules in Human Ischemic Cardiomyopathy. PCDHGA3 Implications in Decreased Stroke Volume and Ventricular Dysfunction. PLoS One11, e0160168. 10.1371/journal.pone.0160168

  • 28

    PanchoA.AertsT.MitsogiannisM. D.SeuntjensE. (2020). Protocadherins at the Crossroad of Signaling Pathways. Front. Mol. Neurosci.13, 117. 10.3389/fnmol.2020.00117

  • 29

    PederickD. T.RichardsK. L.PiltzS. G.KumarR.Mincheva-TashevaS.MandelstamS. A.et al (2018). Abnormal Cell Sorting Underlies the Unique X-Linked Inheritance of PCDH19 Epilepsy. Neuron97, 5966. e55. 10.1016/j.neuron.2017.12.005

  • 30

    PeekS. L.MahK. M.WeinerJ. A. (2017). Regulation of Neural Circuit Formation by Protocadherins. Cell Mol. Life Sci.74, 41334157. 10.1007/s00018-017-2572-3

  • 31

    PolancoJ.Reyes-VigilF.WeisbergS. D.DhimitrukaI.BrusésJ. L. (2021). Differential Spatiotemporal Expression of Type I and Type II Cadherins Associated with the Segmentation of the Central Nervous System and Formation of Brain Nuclei in the Developing Mouse. Front. Mol. Neurosci.14, 633719. 10.3389/fnmol.2021.633719

  • 32

    QinX.LiJ.WuT.WuY.TangX.GaoP.et al (2019). Overall and Sex-specific Associations between Methylation of the ABCG1 and APOE Genes and Ischemic Stroke or Other Atherosclerosis-Related Traits in a Sibling Study of Chinese Population. Clin. Epigenetics11, 189. 10.1186/s13148-019-0784-0

  • 33

    SanesJ. R.ZipurskyS. L. (2020). Synaptic Specificity, Recognition Molecules, and Assembly of Neural Circuits. Cell181, 536556. 10.1016/j.cell.2020.04.008

  • 34

    SchaarschuchA.HertelN. (2018). Expression Profile of N-Cadherin and Protocadherin-19 in Postnatal Mouse Limbic Structures. J. Comp. Neurol.526, 663680. 10.1002/cne.24359

  • 35

    SchianoC.BenincasaG.FranzeseM.Della MuraN.PaneK.SalvatoreM.et al (2020). Epigenetic-sensitive Pathways in Personalized Therapy of Major Cardiovascular Diseases. Pharmacol. Ther.210, 107514. 10.1016/j.pharmthera.2020.107514

  • 36

    ShenY.PengC.BaiQ.DingY.YiX.DuH.et al (2019). Epigenome-Wide Association Study Indicates Hypomethylation of MTRNR2L8 in Large-Artery Atherosclerosis Stroke. Stroke50, 13301338. 10.1161/STROKEAHA.118.023436

  • 37

    Soriano-TárragaC.LazcanoU.Giralt-SteinhauerE.Avellaneda-GómezC.OisÁ.Rodríguez-CampelloA.et al (2020). Identification of 20 Novel Loci Associated with Ischaemic Stroke. Epigenome-wide Association Study. Epigenetics15, 988997. 10.1080/15592294.2020.1746507

  • 38

    WanS.LiuL.RenB.QuM.WuH.JiangW.et al (2021). DNA Methylation Patterns in the HLA-DPB1 and PDCD1LG2 Gene Regions in Patients with Autoimmune Thyroiditis from Different Water Iodine Areas. Thyroid31, 17411748. 10.1089/thy.2021.0221

Summary

Keywords

ischemic stroke, epigenome-wide association study, DNA methylation, gene ontology, cadherins

Citation

Sun H, Xu J, Hu B, Liu Y, Zhai Y, Sun Y, Sun H, Li F, Wang J, Feng A, Tang Y and Zhao J (2022) Association of DNA Methylation Patterns in 7 Novel Genes With Ischemic Stroke in the Northern Chinese Population. Front. Genet. 13:844141. doi: 10.3389/fgene.2022.844141

Received

27 December 2021

Accepted

18 March 2022

Published

11 April 2022

Volume

13 - 2022

Edited by

Cunyou Zhao, Southern Medical University, China

Reviewed by

Yi Huang, Zhejiang University School of Medicine, China

Yingying Xu, The Second Hospital of Shandong University, China

Updates

Copyright

*Correspondence: Jingbo Zhao, ; Ying Tang,

This article was submitted to Epigenomics and Epigenetics, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics