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

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)(11)(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.

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, lowdensity 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 b value, which was calculated as b= max (signal B, 0)/[max(signal A, 0)+ max (signal B, 0)+ 100]. The b values range from 0-1, indicating completely unmethylated to totally methylated DNA (21).

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 Nshore, 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.

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).

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. b-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.
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 FIGURE 1 | Volcano plot for differential DNA methylation analysis of ICH group and control group. The x-axis shows the DNA methylation difference (delta b), 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.
(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). 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.

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 crossvalidation 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).      Table 3).

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).

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-1phosphate receptor), the expression of which can improve the 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. 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-b1) 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 Ca 2+ , mitogen-activated protein kinase (MAPK) kinase and nuclear factor-kB (NF-kB) 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 signalregulated kinases (ERK1/2) and c-Jun amino-terminal kinases (JNKs). Activation of the MAPK pathway protects against ICHinduced 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-kB 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 Ca 2+ signaling pathway by glycogen synthase kinase-3b (GSK-3b) 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-a (TNF-a) and VEGF play important roles (35)(36)(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, TGFb1, 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 TGFb1 expression dramatic reduces surface KLRK1 expression associated with impaired NK cytotoxicity (39). Previous studies suggest that polymorphisms and haplotypes in the TGFb1 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 TGFb1 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, TGFb1, 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.

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.