ORIGINAL RESEARCH article

Front. Immunol., 13 August 2021

Sec. Cancer Immunity and Immunotherapy

Volume 12 - 2021 | https://doi.org/10.3389/fimmu.2021.702244

Genome-Wide DNA Methylation Pattern in Whole Blood Associated With Primary Intracerebral Hemorrhage

  • 1. Department of Neurology, Xiangya Hospital, Central South University, Changsha, China

  • 2. Department of Geriatrics, Second Xiangya Hospital, Central South University, Changsha, China

  • 3. Department of Neurology, Second Xiangya Hospital, Central South University, Changsha, China

Article metrics

View details

18

Citations

4,2k

Views

2,1k

Downloads

Abstract

Primary intracerebral hemorrhage (ICH) is a significant cause of morbidity and mortality throughout the world. ICH is a multifactorial disease that emerges from interactions among multiple genetic and environmental factors. DNA methylation plays an important role in the etiology of complex traits and diseases. We used the Illumina Infinium Human Methylation 850k BeadChip to detect changes in DNA methylation in peripheral blood samples from patients with ICH and healthy controls to explore DNA methylation patterns in ICH. Here, we compared genomic DNA methylation patterns in whole blood from ICH patients (n = 30) and controls (n = 34). The ICH and control groups showed significantly different DNA methylation patterns at 1530 sites (p-value < 5.92E-08), with 1377 hypermethylated sites and 153 hypomethylated sites in ICH patients compared to the methylation status in healthy controls. A total of 371 hypermethylated sites and 35 hypomethylated sites were in promoters, while 738 hypermethylated sites and 67 hypomethylated sites were in coding regions. Furthermore, the differentially methylated genes between ICH patients and controls were largely related to inflammatory pathways. Abnormalities in the DNA methylation pattern identified in the peripheral blood of ICH patients may play an important role in the development of ICH and warranted further investigation.

Introduction

ICH is a significant healthcare and financial burden on society and the families of patients. In Western populations, approximately 10% of all strokes are due to ICH (1), whereas 18.8%-55.4% of all strokes are due to ICH in China (2). ICH is a complex disease caused by environmental and genetic factors (3). Several studies have shown changes in gene expression in the brain following ICH (4, 5). However, the mechanism by which ICH causes these changes remains unknown. There has been recent increasing interest in the roles of epigenetic mechanisms and the interaction between the genome and environment in human diseases (6). Studies have confirmed that epigenetic modifications affect the regulation of gene expression and development, and epigenetics is predicted to be an important part of future genetic studies of ICH (7). However, there are few reports on epigenetic mechanisms in ICH.

DNA methylation refers to the formation of a covalent bond between the 5’ carbon atom of the cytosine of a DNA CpG (5’-cytosine guanine-3’) dinucleotide and a methyl group by DNA methyltransferase, forming 5-methylcytosine (8). DNA hypermethylation can directly inhibit transcription or indirectly inhibit gene expression through transcriptional silencing (9). Methylated DNA in the human genome is an epigenetic marker of paramount importance for normal development. Changes in DNA methylation are associated with cardiovascular and cerebrovascular diseases, such as coronary heart disease, atherosclerosis, and cerebral infarction (10–12). Hypertension is the most important risk factor for ICH (13). Recent studies have found that the overall DNA methylation level of the peripheral blood of patients with essential hypertension was lower than that of the peripheral blood of healthy people (14). Furthermore, an epigenetic study showed that the level of global DNA hydroxymethylation in mouse brains was decreased after ICH (15). While DNA methylation has been linked to gene silencing, DNA hydroxymethylation is associated with actively transcribed genes (16). Therefore, we hypothesized that changes in DNA methylation play an important role in the development of ICH. DNA methylation exhibits strong tissue specificity, and studies of methylation patterns in peripheral blood samples from patients with many neurological diseases, such as Down Syndrome (17) and Parkinson’s disease (18), have indicated that peripheral blood biomarkers are easily used in clinical practice. However, whether DNA methylation is changed in ICH injury remains unknown. In this study, we identified altered genome-wide DNA methylation patterns in whole blood from patients with ICH for the first time.

Materials and Methods

Subjects

This case-control study recruited consecutive Han Chinese patients with acute ICH diagnosed at the Department of Neurology in Xiangya Hospital in Changsha, between August 2016 and October 2017. The inclusion criteria were as follows: diagnosis of ICH confirmed by computed tomography (CT) and/or magnetic resonance imaging (MRI); typical areas exhibiting ICH, especially hematoma in the basal ganglia; a clear history of hypertension; 18 years of age or older; and signed informed consent. Cases of mixed cerebrovascular disease (postinfarct hemorrhage or infarct after hemorrhage); ICH caused by arteritis, trauma, drug, tumor, cerebrovascular malformation, cerebral amyloid angiopathy, sinus thrombosis or aneurysm; and patients more than 2 weeks after the onset of ICH were excluded (19, 20). Age- and sex-matched healthy control subjects who had a routine physical check-up for the purpose of health maintenance at Xiangya Hospital were enrolled in the study. The control subjects had no symptoms or history of stroke, autoimmune disease, peripheral atherosclerotic disease, coronary artery disease or hematological disease. Finally, 30 patients and 34 controls were enrolled in the study following the procurement of written informed consent from all subjects. All research protocols were approved by the Ethics Committee of Central South University (No.201503225).

Information related to gender, age, body mass index (BMI), systolic blood pressure (SBP), diastolic blood pressure (DBP), fasting blood glucose (FBG), total cholesterol, triglyceride characteristics, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol was recorded on admission. The hematoma volumes on admission to hospital were determined from base-line computed tomography scans (SOMATOM® Definition Flash; Siemens Healthcare, Erlangen, Germany) by the ABC/2 method for all patients. ICH patients and/or their caregivers were contacted and interviewed at 6 months after index ICH by telephone. We defined favorable outcome as achieving and retaining during entire follow-up an mRS ≤ 3.

Whole Blood DNA Extraction

Whole blood DNA was extracted from the blood samples of patients and control subjects using a commercial kit (Beijing Adly Biological Company, Beijing, China) according to the manufacturer’s instructions. One microliter of DNA was diluted 10 times with distilled water and shaken evenly, and the DNA concentration (µg/µl) was measured using a NanoDrop spectrophotometer; DNA samples with an OD260/280 ratio from 1.7 to 2.0 were used for further experimentations. The integrity of all DNA samples was assessed by agarose gel electrophoresis and ethidium bromide staining. The main band of the sample electrophoresis is clear, not less than 10kb, no obvious degradation.

BeadChip DNA Methylation Assay

An Illumina Human Methylation 850k BeadChip Kit was used to determine the pattern and extent of DNA methylation in the samples. DNA samples were processed by a series of steps according to the manufacturer’s protocols: DNA denaturation, whole genome amplification, fragmentation, precipitation, resuspension and hybridization. Signals (grayscale) from the final hybridization products were captured with the iScan system (Illumina, Inc.). The extent of methylation was determined by the β value, which was calculated as β= max (signal B, 0)/[max(signal A, 0)+ max (signal B, 0)+ 100]. The β values range from 0–1, indicating completely unmethylated to totally methylated DNA (21).

Validation of Methylation by Pyrosequencing

After analysis by Methylation 850k BeadChip, 6 CpG loci methylated to different extents (cg08909806, cg26620021, cg20460771, cg24507266, cg05799811, cg26934362) were selected. Then, we selected blood samples from 15 ICH cases and 15 sex- and age-matched control cases for bisulfite conversion of the target region using a PyroMark PCR kit (Qiagen, 59104). Pyrosequencing was performed using primers designed by PyroMark Assay Design 2.0 software (Qiagen, CA, USA). Pyrosequencing was performed using PyroMark Q96/48 ID (Qiagen, CA, USA).

Data and Bioinformatics Analyses

BeadStudio Methylation Module v3.2 software (Illumina, Inc.) (http://www.illumina.com) was used to analyze differential DNA methylation sites. Differential methylation sites with p-value<5.92 E-08 according to the Bonferroni correction method were screened to reduce the occurrence of class II statistical errors and ensure the accuracy of the research results. The R (https://www.r-project.org/, version 4.1.0) was used to perform all data analysis. Methylation of gene promoters was defined as methylation in the TSS1500 (upstream), TSS200 (upstream), 5’ UTR and first exon regions. Methylation involving CpG island structures was defined as methylation of a CpG island or the N-shore, S-shore, N-shelf and S-shelf regions, with shores occurring 0–2 kb from a CpG island and shelves occurring 2–4 kb from a CpG island. Differentially methylated genes were selected using the DAVID bioinformatics database (http://david.abcc.ncifcrf.gov/home.jsp/). Pathways enriched in differentially methylated genes were determined using the Gene Ontology (GO) (http://geneontology.org/) and KEGG (http://www.genome.jp/kegg/) databases, and GO terms and KEGG pathways with a Bonferroni corrected p-value of less than 0.05 were reported. Used Pearson’s chi-squared test to compare the distribution of hypomethylation and hypermethylation sites with all probes. The relationship between DNA methylation and clinical outcomes/hematoma volumes was analyzed by Pearson correlation analysis. Measured data are expressed as the mean ± SD. Differences in counted data were statistically analyzed with the chi-square test, and differences in measured data were analyzed by Student’s t-test. p-value < 0.05 indicated statistical significance.

Results

Characteristics of Cases and Controls

The total sample population of 64 individuals comprised 30 (aged 41–75 years, the mean age is 54.7) unrelated patients diagnosed with ICH and 34 (aged 33–74 years, the mean age is 54.9) age- and sex-matched unrelated controls (Table 1).

Table 1

Patients n=30Controls n=34p
Age54.90 ± 7.9154.71 ± 9.650.931
Sex (M/F)19/1123/110.717
SBP/mmHg159.7 ± 23.82138.53 ± 23.590.000*
DBP/mmHg93.7 ± 11.9886.97 ± 12.760.036*
hypertension30/3016/340.000*
Diabetes1/300/340.283
Smoking11/3010/340.537
Drinking9/309/340.754
TC/(mmol/L)4.97 ± 1.164.72 ± 0.810.396
TG/(mmol/L)1.96 ± 0.602.07 ± 1.310.665
LDL-C/(mmo l/L)3.20 ± 0.902.65 ± 0.730.021
HDL-C/(mmol/L)1.11 ± 0.301.19 ± 0.260.245
Hematoma volumes/ml16.95 ± 11.88/

Demographic and clinical characteristics of ICH and control group.

Data presented as mean± SD or n (%).

SBP, systolic blood pressure; DBP, diastolic blood pressure; TC, total cholesterol; TG, triglycerides; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol.

*P<0.05.

Overall Changes in Blood Genomic DNA Methylation in Patients With ICH

To assess differences in methylation between the ICH cases and the control group, we conducted genome-wide analysis of DNA methylation in samples from these groups using the Illumina Infinium Human Methylation 850k BeadChip, which allows methylation-specific hybridization to an array of 844668 DNA methylation sites spanning the entire human genome. β-values for sites exhibiting differential DNA methylation between controls (n = 34) and patients with ICH (n= 30) are reported. Statistical analysis showed a total of 1530 loci showing differential DNA methylation between the two groups (defined using cutoff p-values < 5.92 E-08) (Figure 1). These substantially methylated sites appear to be randomly distributed on 22 chromosomes (Figure 2). Of these differential methylation sites, 1377 sites were hypermethylated, and 153 sites were hypomethylated (yielding a ratio of 9) compared to the methylation status of those sites in the control samples (Figure 3). The 10 methylation sites with the most significant differences in methylation level are listed in Table 2.

Figure 1

Figure 1

Volcano plot for differential DNA methylation analysis of ICH group and control group. The x-axis shows the DNA methylation difference (delta β), the y-axis shows the – log10 p-value of each CpG site. The dashed line indicates statistically significant (P <5.92 E-08). Red represents hypermethylation sites, blue represents hypomethylation sites, and gray represents no significant.

Figure 2

Figure 2

Differential methylation site (P <5.92 E-08) circos diagram. The distribution of 1530 differential DNA methylation sites on 22 chromosomes, the outer layer is the physical location of the chromosome; the middle layer is the chromosome segment; the inner layer is the methylation of the site, red indicates hypermethylation, and green indicates hypomethylation.

Figure 3

Figure 3

Heat map generated from clustering-analysis of microarray data illustrating differentially methylated DNA sites in blood whole genome in the ICH patient group (n = 30) relative to the control group (n = 34). Blue and red represent levels of hypomethylation and hypermethylation, respectively, whereas white indicates no change in methylation, relative to control.

Table 2

TargetIDdelta βP valueCHRGROUPGENE_NAME
Top ten hypermethylated sitescg029024120.1954921575.03851E-0918BodyATP9B
cg269343620.1482529418.89275E-105BodyLPCAT1
cg144071790.1410039221.26452E-107BodyADAP1
cg057998110.1363352947.13005E-091BodyCD247
cg156932990.1340058827.95408E-0912TSS1500LOC360030
cg276420020.1324666672.04617E-0811BodyCLMP
cg179722130.1301666671.36499E-081BodyS1PR1
cg135288730.1295411761.0542E-096BodyGMDS
cg206738670.1292470597.13057E-099BodyMOB3B
cg234748900.1291450983.7247E-0812BodyBCL7A
Top ten hypomethylated sitescg23351010-0.13464.09969E-083TSS1500;BodyFLJ10213;PPP4R2
cg03330490-0.134525498.20857E-0910BodyHHEX
cg08301869-0.1294372555.27767E-084BodyTBC1D14
cg11546475-0.1244764718.93678E-0911BodyCAT
cg07801620-0.1224568637.25833E-0981stExon;Body;5’UTRPLEC1
cg23749353-0.1120901965.59553E-0912BodyCORO1C
cg05756780-0.1119274514.17303E-081BodyIL6R
cg09447621-0.1103450985.38782E-0811BodyB3GAT1
cg16284789-0.1075333334.42562E-081BodyIL10
cg26246740-0.10349.43931E-0910TSS1500WDFY4

Top ten hyper- and hypo- methylated sites in whole blood genome in ICH.

Biological information was available for 1157 sites of the 1530 positive sites, and there were 5 sites that overlap with genetic variant regions (22). Further analysis of DNA functional domains among the hypermethylated sites revealed 414 sites located in noncoding regions (27%), 738 sites in coding regions (47%), 377 sites in promoters (24%) and 34 in the termini of 3’ untranslated regions (UTRs) (2%). Among the hypomethylated sites, 45 sites were in noncoding intergenic regions (28%), 76 sites were in coding regions (46%), 36 sites were in promoters (22%), and 7 (4%) were in the end of a 3’ UTR. There was no significant difference in the distribution of hypomethylation and hypermethylation sites compared with all probes on the 850k Illumina array (Figure 4).

Figure 4

Figure 4

Analysis of differentially methylated DNA sites in reference to genomic structure domains in ICH. (A) The distribution of hypermethylation sites, (B) the distribution of hypomethylation sites. (C) Compared the distribution of hypomethylation and hypermethylation sites with all probes on the Illumina array.

Among the 1530 positive sites, 406 were in promoter regions (Figure 5), and among these sites, 371, accounting for 91.38% of the total, were hypermethylated, and 35, accounting for 8.62% of the total, were hypomethylated. Among the 406 methylation sites in promoter regions, 42 were located in a CpG island, 46 were in a north (N)-shore, 59 were in a south (S)-shore, 227 were in an open area, 14 were in a north (N)-shelf, and 18 were in a south (S)-shelf.

Figure 5

Figure 5

Heat map generated from clustering-analysis of microarray data illustrating differentially methylated DNA sites in the promoter region.

Validation of Methylation by Pyrosequencing

Because previous studies found a good correlation between percent methylation determined by Illumina array and pyrosequencing (23), we randomly selected 6 differential DNA methylation sites for cross-validation with pyrosequencing. Data on both the differential DNA methylation sites and the mean signals for these sites collected by Illumina Infinium Human Methylation 850k BeadChip and pyrosequencing were analyzed. Data on individual sites from Illumina 850k BeadChip array and pyrosequencing were highly correlated, and the means of the differential DNA methylation sites ranged from 0.899 to 0.971 (p-value<0.001) (Figures 6A–F).

Figure 6

Figure 6

(A–F) Methylation chip and pyrosequencing correlation scatter plot Shown are degrees of methylation in 6 CpG loci (cg08909806, cg26620021, cg20460771, cg24507266, cg05799811, cg26934362) reported by Illumina Infinium Human Methylation 850k BeadChip (Y axis, ratio) and pyrosequencing (X axis, percentage) assays. Blue and red dots in each panel plot values of methylation at a given site among individual patients with ICH (n=15) and controls (n=15).

Ontological Profiling of Differentially Methylated Genes in the Blood of Patients With ICH

The 1157 sites with available biological data corresponded to 949 genes. The corresponding biological functions of these genes were obtained from the Database for Annotation, Visualization, and Integrated Discovery (DAVID) database. All genes were enriched in 69 cellular component (CC) terms, and gene enrichment in the following 7 terms was significant (Bonferroni corrected p-value<0.05): (1) GO:0005829, cytosol; (2) GO:0005925, focal adhesion; (3) GO:0015629, actin cytoskeleton; (4) GO:0016020, membrane; (5) GO:0005737, cytoplasm; (6) GO:0005886, plasma membrane; and (7) GO:0016023, cytoplasmic, membrane-bounded vesicle. Among 179 biological process (BP) terms, the following 4 terms were significantly enriched in differentially methylated genes (Bonferroni corrected p-value<0.05): (1) GO:0007165, signal transduction; (2) GO:0043547, positive regulation of GTPase activity; (3) GO:0035556, intracellular signal transduction; and (4) GO:0016477, cell migration. Among the molecular function (MF) terms, 63 terms were enriched in the differentially methylated genes, and 1 of these terms was significantly enriched (Bonferroni corrected p-value<0.05): GO:0005515, protein binding (Figure 7).

Figure 7

Figure 7

Partial GO pathway enrichment analysis bubble chart. The vertical axis indicates the name of the GO pathway, the horizontal axis is the enrichment factor, and the size of the circle indicates the number of differential methylation-related genes involved in the pathway, The color response of each channel takes p value of – log10.

Genomic Network Profiling of Differential Methylated Genes in the Blood of Patients With ICH

Through bioinformatics assessment against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/), a total of 90 genomic networks/pathways were predicted to exhibit biological interplay between/among differentially methylated genes in ICH, and the 5 following networks/pathways were significant (Bonferroni corrected p-value<0.05): (1) hsa04660, T cell receptor signaling pathway; (2) hsa04650, natural killer (NK) cell-mediated cytotoxicity; (3) hsa04662, B cell receptor signaling pathway; (4) hsa04071, sphingolipid signaling pathway; and (5) hsa05205, proteoglycans in cancer (Table 3).

Table 3

CategoryTermCountBonferroni P value GenesFold Enrichment
KEGG_PATHWAYhsa04660:T cell receptor signaling pathway236.05E-07PTPRC, PTPN6, ITK, TNF, VAV3, CD8A, MAP2K2, NFKBIB, PIK3CD, CD247, RAF1, IL10, NCK2, PRKCQ, MAPK1, FYN, GSK3B, PAK4, ZAP70, PIK3R5, PPP3CA, NFATC2, AKT24.5922908
KEGG_PATHWAYhsa04650:Natural killer cell mediated cytotoxicity243.33E-06PRKCA, PRF1, ITGAL, PTPN6, TNF, VAV3, MAP2K2, CD247, PIK3CD, KLRK1, RAF1, FASLG, GZMB, MAPK1, TNFRSF10D, FYN, ZAP70, KLRC4-KLRK1, PIK3R5, PPP3CA, NFATC2, KLRD1, TYROBP, SH3BP24.045667447
KEGG_PATHWAYhsa04662:B cell receptor signaling pathway159.76E-04PTPN6, VAV3, MAP2K2, NFKBIB, PIK3CD, RAF1, MAPK1, GSK3B, CD81, PIK3AP1, PIK3R5, CD79B, PPP3CA, NFATC2, AKT24.470755694
KEGG_PATHWAYhsa04071:Sphingolipid signaling pathway209.93E-04PRKCA, TNF, SGPP1, MAP2K2, PPP2R5C, PIK3CD, RAF1, MAPK10, SGMS1, ADORA1, TRADD, TNFRSF1A, MAPK1, S1PR1, FYN, S1PR5, PIK3R5, PPP2R5E, PPP2R2C, AKT23.427579365
KEGG_PATHWAYhsa05205:Proteoglycans in cancer250.007635914TNF, CAMK2G, RPS6KB2, FASLG, TGFB1, CTNNB1, IGF1R, CD44, PRKACA, PIK3R5, GPC1, CAMK2A, AKT2, PRKCA, PTPN6, ARHGEF1, MAP2K2, PIK3CD, HSPG2, RAF1, DDX5, ITPR1, MAPK1, SMO, RRAS22.570684524

Differential methylation-associated gene KEGG pathway TOP5 in ICH.

The Relationship Between DNA Methylation Patterns and Clinical Outcomes, Hematoma Volumes, Immunological Homeostasis

The top 10 methylation sites correlated with clinical outcomes and hematoma volumes were showed in the Table 4. In order to explore the role of immunological homeostasis in the ICH methylation pattern, we screened out the DNA methylation sites in the promoter region related to clinical prognosis, and found the DNA methylation sites involved in immunological homeostasis, including cg24270157(IL15), cg23343408 (TLR5), which are negatively correlated with the mRS scores (Figure 8).

Table 4

TargetIDRP valueCHRGROUPGENE_NAME
Top 10 methylation sites correlated with clinical outcomescg041285630.7599391.1E-066TSS1500AMD1
cg21112259-0.732674.2E-06125’UTR;BodyCSAD
cg07624428-0.723446.3E-0616
cg137499390.7219656.7E-0611TSS1500SLC25A22
cg108434260.7142269.3E-061
cg03770703-0.710111.1E-051TSS1500NMNAT2
cg27084073-0.709931.1E-0518
cg24764979-0.704841.4E-0516TSS1500;5’UTR;1stExonGRIN2A
cg216696530.7039481.4E-0512
cg13110300-0.701061.6E-056
Top 10 methylation sites correlated with hematoma volumescg07877097-0.802079.9E-086
cg01174235-0.735883.6E-067BodyWBSCR17
cg169091240.7272685.3E-0612TSS200NR2C1
cg20342930-0.723616.2E-063BodySPATA16
cg080490600.7207117.1E-063TSS200TMEM115
cg236847110.7179268.0E-0611stExon;5’UTRISG15
cg22716633-0.713829.5E-064BodyPPP2R2C
cg00617889-0.712571.0E-055
cg04662567-0.710241.1E-056
cg26260819-0.709111.2E-0511Body;1stExonCRY2

Top 10 methylation sites correlated with clinical outcomes and hematoma volumes.

Figure 8

Figure 8

The DNA methylation sites involved in immunological homeostasis and mRS correlation scatter plot. degrees of methylation in cg24270157, cg23343408 are negatively correlated with the 3 months mRS score (R=-0.58, R=-0.59), red dots indicate the favorable outcome patients and the blue means poor outcome.

Discussion

This is the first report to identify an altered DNA methylation pattern in the whole blood of ICH patients. Using the Illumina Infinium Human Methylation 850k BeadChip, 844668 DNA methylation sites were tested. Differentially methylated sites with a Bonferroni corrected p-value < 0.05 were considered significant. We identified 1530 differential DNA methylation sites, among which 1377 sites were significantly hypermethylated, and 153 sites were hypomethylated (yielding a ratio of 9) in ICH patients relative to their methylation in the healthy control group. The methylation of sites around promoters has the strongest inhibitory effect on transcription (24). In this study, among the 1530 positive sites, 406 were located in promoter regions; among these 406 sites, 371 were hypermethylated, and 35 were hypomethylated in patients with ICH, suggesting predominantly increased methylation at these sites. Such changes could enhance the expression of the corresponding genes. This study identified the genome-wide DNA methylation pattern of peripheral blood from patients with ICH and found that the main differences in methylation between ICH patients and normal controls are hypermethylation.

The bioinformatics profile of our data indicates that the differentially methylated regions are correlated to 949 known functional genes. The top ten hypermethylated and top ten hypomethylated genes (Table 2), such as S1PR1 (sphingosine-1-phosphate receptor), the expression of which can improve the clinical outcome of ICH, are related to the pathophysiology of ICH (25). The levels of 5-hydroxymethylcytosine in the CpG-rich regions of the AKT2(serine/threonine kinase 2), PDPK1 (3-phosphoinositide dependent protein kinase 1)and VEGF (ascular endothelial growth factor) genes were significantly decreased and at their minimum levels, suggesting the coordinated upregulation of proinflammatory/anti-inflammatory signaling pathways and neuronal signaling systems in ICH mouse brains; furthermore, changes in the expression of these genes contribute to cell death after ICH (15). This confirms the involvement of epigenetic mechanisms such as DNA hydroxymethylation in pathological changes following ICH. Our current data revealed the hypermethylation of the promoters of the AKT2 and PDPK1 genes in ICH patients relative to their methylation in the control group. Studies have shown that increased expression of the TSPO (translocator protein) gene after ICH may be an intrinsic mechanism that prevents an increase in the inflammatory response (26). The polymorphism rs10940495 in the IL6ST (interleukin 6 signal transducer) gene is associated with functional outcome 6 months after ICH (27), transforming growth factor beta 1(TGF-β1) modulates microglia-mediated neuroinflammation after ICH and promotes functional recovery (28). The relationship between these genes and ICH needs further research.

KEGG pathway analysis showed that these genes form interconnected networks involved in the T cell receptor signaling pathway, NK cell (natural killer cell)-mediated cytotoxicity, the B cell receptor signaling pathway, the sphingolipid signaling pathway and other networks. The T cell receptor signaling network, which is a proinflammatory signaling pathway, contains three major signaling pathways, the Ca2+, mitogen-activated protein kinase (MAPK) kinase and nuclear factor-κB (NF-κB) signaling pathways, that mobilize transcription factors that are critical for gene expression and for T cell growth and differentiation (29). The MAPK family is divided into three subfamilies: p38 MAPKs, extracellular signal-regulated kinases (ERK1/2) and c-Jun amino-terminal kinases (JNKs). Activation of the MAPK pathway protects against ICH-induced secondary brain injury (30). In the present study, genes in the MAPK pathway (RAF1, MAPK1 and MAP2K2) were differentially methylated in ICH patients relative to their methylation status in the control group.

NF-κB signaling pathways plays a key role in immunological homeostasis (31). Previous research indicated a relationship between NF-kB and the pathobiology of perilesional cell death after ICH. In our study, it was found that the methylation pattern of partially DNA methylation sites located in the promoter region was negatively correlated with the follow-up mRS scores of ICH. And these methylation site-associated genes (IL15, TLR5) are related to NF-kB signaling pathways (32, 33). We speculate that hypomethylation of the promoter region of these genes may affect their expression and participate in immunological homeostasis of ICH, which is worth-study in the future.

The PI3K/Akt signaling pathway can suppress the Ca2+ signaling pathway by glycogen synthase kinase-3β (GSK-3β) in the T cell receptor signaling networks, and statins can increase brain-derived neurotrophic factor (BDNF) and VEGF expression, which activated the PI3K/Akt-mediated signaling signaling pathway and improved neurological function in an autologous ICH model (34). Among the various cytokines that mediate inflammatory reactions in ICH, the proinflammatory cytokines IL, tumor necrosis factor-α (TNF-α) and VEGF play important roles (35–37). In addition, altered methylation in the T cell receptor signaling networks may change the corresponding inflammatory pathways in ICH.

In the present study, genes connected with the NK pathway (KLRK1, TGFβ1, ITGAL) were differentially methylated in ICH patients relative to their methylation status in the control group. KLRK1, the receptor for UL16-binding proteins, is expressed on primary NK cells (38), and TGFβ1 expression dramatic reduces surface KLRK1 expression associated with impaired NK cytotoxicity (39). Previous studies suggest that polymorphisms and haplotypes in the TGFβ1 gene are associated with cerebral infarction and ICH (40, 41). Therefore, far more research on the relationship between changes in the methylation of KLRK1 and TGFβ1 with ICH are needed. ITGAL gene, which encodes CD11a, is important in adhesive interactions between T cells and other cells of the immune system in the NK pathway. Furthermore, ITGAL gene promoter methylation and chromatin structure may contribute to the tissue-specific expression of CD11a (42). However, the relationship between changes in ITGAL methylation and ICH is unclear.

There are several limitations to this study. First, we used DNA from peripheral blood leukocytes, representing groups of different cells with different epigenetic profiles. The purpose of this study is to find out some potential epigenetic bio-marks of ICH, future research on more specific types of samples (such as NK-cells, CD8+ T-cells, Neutrophils), which by sorting of subsets by flow cytometry, may provide us with more detailed results. Secondly, the sample size needs to be increased to improve the reliability of the research in the future. Studies have found that primary damage following ICH is mainly the result of the mass effect of hematoma (43). Secondary damage in the brain following ICH occurs through the induction of hematoma toxicity, oxidative stress, cerebral edema and inflammation. Genes associated with the differentially methylated sites identified in this study are primarily involved in inflammatory pathways. DNA methylation is speculated to be involved in the inflammatory response-related mechanism of ICH, which provides new ideas for epigenetic research on ICH. Therefore, further evidence to support proteomics-related studies is still needed.

Conclusion

We identified 1530 differentially methylated sites between patients with ICH and normal controls for the first time with the Human Methylation 850k BeadChip. Some of the differentially methylated sites identified in the present study are associated with promoters and genes (such as RAF1, MAPK1, MAP2K2, AKT2, KLRK1, TGFβ1, ITGAL) important for the biological and metabolic regulation of many different inflammatory pathways. Therefore, DNA methylation may play an important role in the pathogenesis of ICH. We hope that future research will combine genomics and proteomics studies to further elucidate the mechanism of DNA methylation in ICH.

Funding

This study is funded by the National Science & Technology Foundational Resource Investigation Program of China (Grant No. 2018FY100900), the National Natural Science Foundation of China (Grant No.81571151), the Natural Science Foundation of Hunan Province, China (Grant No. 2016JJ2164) and the Natural Science Foundation of Hunan Province, China (Grant No. 2017JJ2356).

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.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI GEO GSE179759.

Ethics statement

The studies involving human participants were reviewed and approved by Medical Ethics Committee of Xiangya Hospital Central South University (201503225). The patients/participants provided their written informed consent to participate in this study.

Author contributions

The co-authors YuZ and HL contributed equally to this article in clinical recruitment,drafted the manuscript and statistical analysis. LZ and YiZ designed the study. SW, MX, JL and LC collected the samples. RC, XW and LW contributed to the reagents, materials, and analysis tools. SW, YS and WX provided technical assistance. All authors contributed to the article and approved the submitted version.

Acknowledgments

The authors thank the patients and their caregivers who took part in the study, as well as the staff at Neurology Department of Xiangya Hospital.

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.

The handling editor declared a shared affiliation with the authors at time of review.

References

  • 1

    ZiaEEngstrMGSvenssonPJNorrvingBPessah-RasmussenH. Three-Year Survival and Stroke Recurrence Rates in Patients With Primary Intracerebral Hemorrhage. Stroke (2009) 40:3567–73. doi: 10.1161/STROKEAHA.109.556324

  • 2

    YangQDNiuQZhouYHLiuYHXuHWGuWPet al. Incidence of Cerebral Hemorrhage in the Changsha Community. A Prospective Study From 1986 to 2000. Cerebrovasc Dis (2004) 17:303–13. doi: 10.1159/000077341

  • 3

    LiuBZhangLYangQ. Genetics of Intracerebral Hemorrhage: Insights From Candidate Gene Approaches. Neurol India (2012) 60:3–8. doi: 10.4103/0028-3886.93581

  • 4

    CarmichaelSTVespaPMSaverJLCoppolaGGeschwindDHStarkmanSet al. Genomic Profiles of Damage and Protection in Human Intracerebral Hemorrhage. J Cereb Blood Flow Metab (2008) 28:1860–75. doi: 10.1038/jcbfm.2008.77

  • 5

    RosellAVilaltaAGarcía-BerrocosoTFernández-CadenasIMontanerJ. Brain Perihematoma Genomic Profile Following Spontaneous Human Intracerebral Hemorrhage. PloS One (2012) 6:e16750. doi: 10.1371/journal.pone.0016750

  • 6

    PortelaAEstellerM. Epigenetic Modifications and Human Disease. Nat Biotechnol (2010) 28:1057–68. doi: 10.1038/nbt.1685

  • 7

    CarpenterAMSinghIPGandhiCDPrestigiacomoCJ. Genetic Risk Factors for Spontaneous Intracerebral Haemorrhage. Nat Rev Neurol (2015) 12:40–9. doi: 10.1038/nrneurol.2015.226

  • 8

    AttwoodJTYungRLRichardsonBC. DNA Methylation and the Regulation of Gene Transcription. Cell Mol Life Sci Cmls (2002) . 59:241–57. doi: 10.1007/s00018-002-8420-z

  • 9

    QureshiIAMehlerMF. Epigenetic Mechanisms Underlying Human Epileptic Disorders and the Process of Epileptogenesis. Neurobiol Dis (2010) 39:53–60. doi: 10.1016/j.nbd.2010.02.005

  • 10

    SharmaPKumarJGargGKumarAPatowaryAKarthikeyanGet al. Detection of Altered Global DNA Methylation in Coronary Artery Disease Patients. DNA Cell Biol (2008) 27:357–65. doi: 10.1089/dna.2007.0694

  • 11

    LeeJParkJHYanBCKimIHChoGJeoungDet al. Effects of Transient Cerebral Ischemia on the Expression of DNA Methyltransferase 1 in the Gerbil Hippocampal CA1 Region. Neurochemical Res (2013) 38:74–81. doi: 10.1007/s11064-012-0890-2

  • 12

    ZainaSHeynHCarmonaFJVarolNSayolsSCondomEet al. DNA Methylation Map of Human Atherosclerosis. Circ: Cardiovasc Genet (2014) 7:692–700. doi: 10.1161/CIRCGENETICS.113.000441

  • 13

    O'DonnellMJChinSLRangarajanSXavierDLiuLZhangHet al. Global and Regional Effects of Potentially Modifiable Risk Factors Associated With Acute Stroke in 32 Countries (INTERSTROKE): A Case-Control Study. Lancet (2016) 388:761–75. doi: 10.1016/S0140-6736(16)30506-2

  • 14

    SmolarekIWyszkoEBarciszewskaAM. Global DNA Methylation Changes in Blood of Patients With Essential Hypertension. Med Sci Monitor (2010) 16:R149–55. doi: 10.1016/j.mehy.2009.07.026

  • 15

    TangYHanSAsakawaTLuoYHanXXiaoBet al. Effects of Intracerebral Hemorrhage on 5-Hydroxymethylcytosine Modification in Mouse Brains. Neuropsychiatr Dis Treat (2016) 12:617–24. doi: 10.2147/NDT.S97456

  • 16

    FiczGBrancoMRSeisenbergerSSantosFKruegerFHoreTAet al. Dynamic Regulation of 5-Hydroxymethylcytosine in Mouse ES Cells and During Differentiation. Nature (2011) 473:398–402. doi: 10.1038/nature10008

  • 17

    BacaliniMGGentiliniDBoattiniAGiampieriEPirazziniCGiulianiCet al. Identification of a DNA Methylation Signature in Blood Cells From Persons With Down Syndrome. Aging (2015) 7:82–93. doi: 10.18632/aging.100715

  • 18

    MasliahEDumaopWGalaskoDDesplatsP. Distinctive Patterns of DNA Methylation Associated With Parkinson Disease. Epigenetics: Off J DNA Methylation Soc (2013) 8:1030–8. doi: 10.4161/epi.25865

  • 19

    WooDSauerbeckLRKisselaBMKhouryJCSzaflarskiJPGebelJet al. Genetic and Environmental Risk Factors for Intracerebral Hemorrhage: Preliminary Results of a Population-Based Study * Editorial Comment: Preliminary Results of a Population-Based Study. Stroke (2002) 33:1190–6. doi: 10.1161/01.STR.0000014774.88027.22

  • 20

    ZengYMaMLiuBXiaJXuHLiuYet al. Association Between ECE1 Gene Polymorphisms and Risk of Intracerebral Haemorrhage. J Int Med Res (2016) 44:444–52. doi: 10.1177/0300060516635385

  • 21

    LongHYFengLKangJLuoZHXiaoWBLongLLet al. Blood DNA Methylation Pattern Is Altered in Mesial Temporal Lobe Epilepsy. Sci Rep (2017) 7:43810. doi: 10.1038/srep43810

  • 22

    PidsleyRZotenkoEPetersTJLawrenceMGClarkSJ. Critical Evaluation of the Illumina Methylationepic Beadchip Microarray for Whole-Genome Dna Methylation Profiling. Genome Biol (2016) 17:208. doi: 10.1186/s13059-016-1066-1

  • 23

    RajendramRFerreiraJCGrafodatskayaDChoufaniSChiangTPuSet al. Assessment of Methylation Level Prediction Accuracy in Methyl-DNA Immunoprecipitation and Sodium Bisulfite Based Microarray Platforms. Epigenetics (2011) 6:410–5. doi: 10.4161/epi.6.4.14763

  • 24

    SmithZDChanMMMikkelsenTSGuHGnirkeARegevAet al. A Unique Regulatory Phase of DNA Methylation in the Early Mammalian Embryo. Nature (2012) 484:339–44. doi: 10.1038/nature10960

  • 25

    SunNShenYHanWShiKWoodKFuYet al. Selective Sphingosine-1-Phosphate Receptor 1 Modulation Attenuates Experimental Intracerebral Hemorrhage. Stroke (2016) 47:1899–906. doi: 10.1161/STROKEAHA.115.012236

  • 26

    BonsackFAlleyneCHSukumari-RameshS. Augmented Expression of TSPO After Intracerebral Hemorrhage: A Role in Inflammation? J Neuroinflamm (2016) 13:151. doi: 10.1186/s12974-016-0619-2

  • 27

    El HusseiniNHoffmanBMBennettERLiYWilliamson TaylorRAHaileyCEet al. Association of IL6ST (Gp130) Polymorphism With Functional Outcome Following Spontaneous Intracerebral Hemorrhage. J Stroke Cerebrovascular Dis (2018) 27:125–31. doi: 10.1016/j.jstrokecerebrovasdis.2017.08.017

  • 28

    TaylorRAChangCFGoodsBAHammondMDGroryBMAiYet al. TGF-β1 Modulates Microglial Phenotype and Promotes Recovery After Intracerebral Hemorrhage. J Clin Invest (2016) 1:280–92. doi: 10.1172/JCI88647

  • 29

    BrownlieRJZamoyskaR. T Cell Receptor Signalling Networks: Branched, Diversified and Bounded. Nat Rev Immunol (2013) 13:257–69. doi: 10.1038/nri3403

  • 30

    WenZMeiBLiHDouYTianXShenMet al. P2X7 Participates in Intracerebral Hemorrhage-Induced Secondary Brain Injury in Rats via MAPKs Signaling Pathways. Neurochemical Res (2017) 42:2372–83. doi: 10.1007/s11064-017-2257-1

  • 31

    AnderssonUTraceyKJ. Reflex Principles of Immunological Homeostasis. Annu Rev Immunol (2012) 30:313–35. doi: 10.1146/annurev-immunol-020711-075015

  • 32

    SmithMFMitchellALiGDingSFitzmauriceAMRyanKet al. TLR2 and TLR5, But Not TLR4, Are Required for Helicobacter Pylori-Induced NF-kB Activation and Chemokine Expression by Epithelial Cells. J Biol Chem (2003) 35:33552–60. doi: 10.1074/jbc.M305536200

  • 33

    YanWFanWChenCWuYFanZChenJet al. IL-15 Up-Regulates the MMP-9 Expression Levels and Induces Inflammatory Infiltration of Macrophages in Polymyositis Through Regulating the NF-kB Pathway. Gene (2016) 1:137–47. doi: 10.1016/j.gene.2016.06.055

  • 34

    YangDHanYZhangJChoppMSeyfriedDM. Statins Enhance Expression of Growth Factors and Activate the PI3K/Akt-Mediated Signaling Pathway After Experimental Intracerebral Hemorrhage. World J Neurosci (2012) 02:74–80. doi: 10.4236/wjns.2012.22011

  • 35

    MracskoEVeltkampR. Neuroinflammation After Intracerebral Hemorrhage. Front Cell Neurosci (2014) 8:388. doi: 10.3389/fncel.2014.00388

  • 36

    YangGShaoG. Elevated Serum IL-11, TNF α, and VEGF Expressions Contribute to the Pathophysiology of Hypertensive Intracerebral Hemorrhage (HICH). Neurological Sci (2016) 37:1253–9. doi: 10.1007/s10072-016-2576-z

  • 37

    ZhengJSunJYangLZhaoBFanZ. The Potential Role of Vascular Endothelial Growth Factor as a New Biomarker in Severe Intracerebral Hemorrhage. J Clin Lab Anal (2017) 31:e22076. doi: 10.1002/jcla.22076

  • 38

    SutherlandCLChalupnyNJSchooleyKVandenBosTKubinMCosmanD. UL16-Binding Proteins, Novel MHC Class I-Related Proteins, Bind to NKG2D and Activate Multiple Signaling Pathways in Primary NK Cells. J Immunol (2002) 2:671–9. doi: 10.4049/jimmunol.168.2.671

  • 39

    LeeJCLeeKMKimDWHeoDS. Elevated TGF-Beta1 Secretion and Down-Modulation of NKG2D Underlies Impaired NK Cytotoxicity in Cancer Patients. J Immunol (2004) 172:7335–40. doi: 10.4049/jimmunol.172.12.7335

  • 40

    TaoHChenGChengGShanX. The Haplotype of the Tgfβ1 Gene Associated With Cerebral Infarction in Chinese. Can J Neurological Sci (2012) 39:626–31. doi: 10.1017/S0317167100015365

  • 41

    KumarPKumarAMisraSSagarRFarooqMKumariRet al. Association of Transforming Growth Factor-β1 Gene C509T, G800A and T869C Polymorphisms With Intracerebral Hemorrhage in North Indian Population: A Case–Control Study. Neurological Sci (2016) 37:353–9. doi: 10.1007/s10072-015-2426-4

  • 42

    LuQRayDGutschDRichardsonB. Effect of DNA Methylation and Chromatin Structure on ITGAL Expression. Blood (2002) 99:4503–8. doi: 10.1182/blood.v99.12.4503

  • 43

    QureshiAIMendelowADHanleyDF. Intracerebral Haemorrhage. Lancet (2009) 373:1632–44. doi: 10.1016/S0140-6736(09)60371-8

Summary

Keywords

intracerebral hemorrhage, epigenome-wide association study, epigenetics, DNA methylation, stroke

Citation

Zhang Y, Long H, Wang S, Xiao W, Xiong M, Liu J, Chen L, Chen R, Wei X, Shu Y, Zeng Y and Zhang L (2021) Genome-Wide DNA Methylation Pattern in Whole Blood Associated With Primary Intracerebral Hemorrhage. Front. Immunol. 12:702244. doi: 10.3389/fimmu.2021.702244

Received

29 April 2021

Accepted

28 July 2021

Published

13 August 2021

Volume

12 - 2021

Edited by

Mingzhu Yin, Central South University, China

Reviewed by

Wen Jiang, Fourth Military Medical University, China; Suyue Pan, Southern Medical University, China

Updates

Copyright

*Correspondence: Le Zhang,

†These authors have contributed equally to this work

This article was submitted to Cancer Immunity and Immunotherapy, a section of the journal Frontiers in Immunology

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