Received Date: September 07, 2018; Accepted Date: September 22, 2018; Published Date: September 29, 2018
Citation: Kazmi N, Elliott HR, Burrows K, Tillin T, Hughes AD, et al. (2018) Associations between High Blood Pressure and DNA Methylation. J Clin Epigenet 4:16. doi: 10.21767/2472-1158.100098
Background: High blood pressure is one of the major risk factors for cardiovascular disease and is influenced by both environmental and genetic factors. Epigenetic processes such as DNA methylation potentially mediate the relationship between genetic factors, the environment and cardiovascular disease. Despite an increased risk of cardiovascular disease in South Asians, it is not clear whether associations between blood pressure and DNA methylation differ between Europeans and South Asians.
Methods and Findings: In this study, we performed an epigenome-wide association study and differentially methylated region analysis to identify DNA methylation sites and regions that were associated with systolic blood pressure, diastolic blood pressure and hypertension. We analyzed samples from 364 European and 348 South Asian men from the Southall and Brent REvisited cohort, measuring DNA methylation from blood using the Illumina Infinium® HumanMethylation450 BeadChip. One CpG site was found to be associated with DBP in trans-ancestry analyses (i.e. both ethnic groups combined), while in Europeans seven CpG sites were associated with DBP. No associations were identified between DNA methylation and either SBP or hypertension. Comparison of effect sizes between South Asian and European EWAS for DBP, SBP and hypertension revealed little agreement between analyses. DMR analysis identified several DMRs including regions with known relationships with CVD and its risk factors.
Conclusion: This study identified differentially methylated sites and regions associated with blood pressure and revealed ethnic differences in these associations. These findings may uncover new molecular pathways which may explain differences in disease risk experienced by South Asians.
High blood pressure (BP) results from abnormalities in the control systems that normally regulate blood pressure . It is one of the strongest risk factors for cardiovascular disease (CVD) which is the leading cause of death worldwide [2,3]. People of South Asian descent have increased risk of CVD compared to Europeans [4,5]. For example, in South Asians living in the United Kingdom, death rate from stroke is between 20% and 25% greater than the rest of the population . It is also established that associations between SBP or DBP and stroke are stronger in South Asians than Europeans . DNA methylation is an example of an environmentally responsive, mitotically stable, epigenetic mark that is associated with biological processes leading to high blood pressure and stroke [7-9]. Candidate gene analysis in cell-line and animal studies has demonstrated a role of DNA methylation in the pathogenesis of hypertension [10-14]. For one of these genes, HSD11B2, methylation has also been associated with hypertension in humans . One recent study found associations between systolic and diastolic BP and DNA methylation among participants of European, Hispanic and African American decent . We aimed to identify DNA methylation associated with SBP, DBP and hypertension in peripheral blood of European and South Asian men in data collected using the Infinium HumanMethylation450 (HM450) BeadChip. We planned to analyse all samples together in a trans-ancestry analysis, then to conduct analysis in each ethnic group separately. We hypothesised that BP associated epigenetic marks would differ between South Asian and European groups, highlighting potential mechanisms explaining the disparity in CVD and stroke risk between the two ethnicities.
SABRE (Southall and Brent REvisited) is a population-based cohort including 4857 people of European, South Asian and African Caribbean origin aged 40 to 69 living in West London, UK . Peripheral blood samples were collected from the Southall participants at baseline (1988-91) for DNA extraction.
In the current analysis, 800 (400 European and 400 South Asian) samples from the SABRE cohort were randomly selected from available baseline samples of good DNA quality. We lost some samples during quality control procedures and some samples were discarded due to missing information. After exclusion, we remained with 712 (364 European and 348 South Asian) individuals. These individuals did not have known diabetes or coronary heart disease at baseline and were stratified by fouryear age group and ethnicity. All individuals were male. The SABRE study predominantly focused on the recruitment of men , for that reason epigenetic analyses were restricted to male participants.
Ethnicity in the SABRE cohort was assigned based on grandparental origins from participant questionnaire. Blood pressure was measured on one occasion (the average of 2 consecutive readings) in the baseline research clinic as described previously .
All participants gave written informed consent. Approval for the baseline study was obtained from Ealing, Hounslow and Spelthorne, Parkside and University College London research ethnics committees. Characteristics of participants are shown in Table 1.
|SBP||122.7 (10.6)||122.0 (10.7)||123.5 (10.3)||0.24|
|DBP||78.4 (17.4)||76.8 (18.1)||80.0 (16.6)||8.3 × 10-5|
|51 (7.1)||52 (7.2)||51 (7.0)||0.06|
|25.8 (3.5)||26.0 (3.3)||25.7 (3.6)||0.3|
|Smokers (ever smokers)||361||112||266||<0.001*|
|Social class (Manual)||473||264||231||0.15*|
P-value: t-test for differences between European and South Asian groups and *P-value=Fisher’s exact test.
Table 1: Distributions of study characteristics included in the EWAS analysis of SBP, DBP and hypertension.
Traits of interest
We investigated SBP, DBP and hypertension as our traits of interest. Hypertension was defined as occurrence of SBP ≥ 140 mm Hg and/or DBP ≥ 90 mm Hg, or receiving medication for hypertension as described previously . The BP protocol was based on the INTERSALT study protocol  and was followed very carefully.
Models were adjusted for age (years), body mass index (BMI) (kg/m2), ethnicity (European or South Asian for trans-ancestry analyses), smoking status (never or ever smoking) and social class (manual or non-manual occupation). Age, ethnicity, smoking status and social class of participants were collected from questionnaires. BMI was calculated from clinic measures of height and weight.
DNA methylation and pre-processing
DNA was bisulfite converted using the Zymo EZ DNA MethylationTM kit (Zymo, Irvine, CA). Following conversion, DNA methylation was measured using the HM450 BeadChip in line with standard protocols at the University of Bristol, UK .
Samples failing QC (average probe detection p value ≥ 0.01) were repeated and if unsuccessful excluded from further analysis. BeadChip intensity data were converted to β-values using the minfi package  in the R statistical programming language. Methylation beta values range from 0 (no cytosine methylation) to 1 (completely cytosine methylated). Raw beta values were normalised using the Functional Normalisation method of the minfi package . We excluded control probes (n=65), any probes with a detection p-value > 0.05 in more than 5% samples, non-CpG probes, polymorphic probes (defined as SNP-overlapping probes, probes with a SNP at the target CpG site, or probes with a SNP at the base next to the target CpG) and probes with a minor-allele frequency (MAF) ≥ 5%; based on UCSC common SNPs track for dbSNP build 137. We further excluded probes that are considered as cross-hybridizing . We applied this stringent CpG filtering because polymorphic and cross-hybridizing probes can interfere with accurate detection of methylation levels. After excluding these features, 402,331 probes remained for the analysis.
Estimation of cell counts
Cell count estimates were derived using the reference-based Houseman method  in the R minfi package  using the Reinius et al. dataset as reference . This method estimates the relative proportions of six white blood cell subtypes (CD4+ T-lymphocytes, CD8+ T-lymphocytes and NK (natural killer) cells, B-lymphocytes, monocytes and granulocytes).
Epigenome-wide association study (EWAS) analysis
We considered a trans-ancestry analysis as the primary model because this provided maximal statistical power by enabling us to include all participants. Our primary hypothesis was that the potential epigenetic mechanisms of elevated blood pressure would be different across these two ethnicities. However, we postulated that these differences may be in the magnitude of effect at key loci, rather than being completely different loci.
Three sets of EWAS analyses were run to identify CpG sites associated with SBP, DBP or hypertension. Normalised, untransformed beta-values, which are on a scale of 0 (completely unmethylated) to 1 (completely methylated) were utilised. Multiple linear regression models were run to evaluate the association between traits of interest and DNA methylation. Analyses were run for each ethnic subgroup and with all samples combined in a trans-ancestry analysis. To remove potential batch effects, principal components (PCs) were generated from methylation beta values using Principal Component Analysis (PCA) [25-27]. The first four PCs were included as covariates in each of the EWAS. These PCs were not associated with exposure and captured 27% of the variance.
Two models were run for each trait of interest: i) an unadjusted model and ii) a model adjusted for confounders (age, BMI, ethnicity (where appropriate), smoking status, social class), estimated cell counts and PCs.
The genomic inflation factor (lambda) was computed and Manhattan plots were generated to compare the genome wide distribution of p-values in EWAS. Multiple testing was accounted for using the false discovery rate (FDR) procedure by Benjamini and Hochberg (BH) . CpG sites with FDR-corrected p-value<0.05 were considered to be associated with the trait of interest. We considered the fully adjusted trans-ancestry models as the primary analysis models.
Differentially methylated region analysis
In addition to EWAS analyses, differentially methylated region (DMR) analyses in relation to SBP, DBP and hypertension were conducted separately using the R package DMRcate . In the DMR analysis, normalised, untransformed beta-values were used and the models were adjusted for confounders, estimated cell counts and PCs. DMRcate groups associated probes into separate DMRs if the gap between nucleotides is ≥1000 base pairs. P-values of associations were adjusted for multiple testing using the BH method.
Data used in these analyses are available from SABRE on request.
Trans-ancestry EWAS analyses
SBP: In the unadjusted EWAS investigating the association between DNA methylation and SBP we identified eight CpG sites with FDR corrected p-value<0.05 (Supplementary Table S1). These CpG sites were cg06639320 near FHL2, cg06192883 in the gene body of MYO5C, cg16867657 near ELOVL2, cg14361627 and cg08097417 near KLF14, cg00494544 in the gene body of PPP1R2, cg24590708 in the gene body of MYO5C and cg22454769 near FHL2. Higher SBP was associated with higher methylation levels at these eight CpG sites. Amongst associated CpG sites, the difference in methylation per unit SBP was 0.03% to 0.08%. After adjustment for potential confounders these associations were markedly attenuated and no longer associated when observing an FDR corrected p-value<0.05 (Supplementary Table S2). DBP: In the unadjusted EWAS investigating the association between DNA methylation and DBP we identified six CpG sites, with FDR corrected p-value<0.05 (Supplementary Table S3). These sites were cg05575921 in the gene body of AHRR, cg12803068 and cg22132788 in the gene body of MYO1G, cg03636183 in the gene body of F2RL3, cg09935388 in the gene body of GFI1 and cg06355652 near CD59. DBP was associated with both higher and lower methylation levels at these CpG sites and the absolute difference in methylation per unit DBP was 0.01% to 0.26%. After adjustment for potential confounders, the associations identified in the unadjusted EWAS were markedly attenuated and no longer associated when observing an FDR corrected p-value<0.05. However, one additional CpG site (cg07598370 near OR5AP2) was identified in the adjusted EWAS, with FDR corrected p-value<0.05 (Figure 1 and Supplementary Table S4). DBP was associated with lower methylation at this CpG site; 0.1% decrease in DNA methylation per 1 mmHg increase in DBP. Hypertension: In the unadjusted EWAS investigating the association between DNA methylation and hypertension, we identified three CpG sites, with FDR corrected p-value<0.05 (Supplementary Table S5). These CpG sites are cg00494544 in the gene body of PPP1R2, cg18120259 in the gene body of LOC100132354 and cg10461004 near LOC285548. Hypertension was associated with both hitgher and lower methylation levels at these CpG sites and the absolute difference in methylation was 0.3% to 2.9% between cases and controls. After adjustment for potential confounders, the associations identified in the unadjusted EWAS were markedly attenuated and no longer associated when observing an FDR corrected p-value<0.05 (Supplementary Table S6).
SBP: In both unadjusted and fully adjusted EWAS investigating the association between SBP and DNA methylation, no CpG sites were associated with SBP after correction for multiple testing (Supplementary Tables S7 and S8).
DBP: In the unadjusted EWAS investigating the association between DNA methylation and DBP, we identified eleven CpG sites with FDR corrected p-value<0.05 (Supplementary Table S9). These CpG sites are cg16241714 near CEBPD, cg05575921 in the gene body of AHRR, cg00006122 near C12orf44, cg04751533 in the gene body of AFAP1, cg13799986 in the gene body of TMEM141, cg06726740 in the gene body of CACNA1C, cg05284742 in the gene body of ITPK1, cg06355652 near CD59, cg04202002 in an intragenic region, cg13691892 near PNPLA7/ MRPL41 and cg25648203 in the gene body of AHRR. DBP was associated with both higher and lower methylation levels at these CpG sites and the absolute difference in methylation per unit DBP was 0.02% to 0.4%. In the fully adjusted EWAS, six of the eleven CpG sites: cg16241714, cg00006122, cg04751533, cg13691892, cg13799986 and cg06355652 were identified with FDR corrected p-value<0.05 (Supplementary Table S10). Their direction of association was consistent with the unadjusted EWAS. In addition to these six CpG sites, cg18901940 near MYPN was found with FDR corrected p-value<0.05. DBP was associated with lower methylation at this CpG site; 0.08% decrease in DNA methylation per 1 mmHg increase in DBP (Figures 2,3 and Supplementary Table S10).
Hypertension: In both unadjusted and fully adjusted EWAS investigating the association between hypertension and DNA methylation, no CpG sites were associated with hypertension after correction for multiple testing (Supplementary Tables S11 and S12).
South Asian EWAS
SBP: In the unadjusted EWAS investigating the association between DNA methylation and SBP, we identified one CpG site (cg07963349 in the gene body of GALR2), with FDR corrected p-value<0.05 (Supplementary Table S13). In the fully adjusted EWAS, this CpG site was no longer associated after multiple corrections but the direction of association was consistent between models (Supplementary Table S14).
DBP: In both unadjusted and fully adjusted EWAS investigating the association between DBP and DNA methylation, no CpG site was associated following correction for multiple testing (Supplementary Tables S15 and S16).
Hypertension: In both unadjusted and fully adjusted EWAS investigating the association between hypertension and DNA methylation, no CpG site was associated following correction for multiple testing (Supplementary Tables S17 and S18).
The genomic inflation factor, lambda, for all fully adjusted models is provided in Table 2.
Table 2: Confounder and cell counts adjusted levels of genomic inflation (Lambda).
We directly compared the direction and magnitude of associations from the fully adjusted EWAS of SBP, DBP and hypertension between trans-ancestry, Europeans and South Asians in a pairwise comparison. We found moderate consistency between the trans-ancestry and Europeans EWAS (goodness-of-fit R2=0.46, R2=0.48 and R2=0.45 for SBP, DBP and hypertension respectively (Supplementary Figure S1)). The goodness of fit was improved when we compared EWAS of trans-ancestry with South Asians (goodness-of-fit R2=0.48, R2=0.50 and R2=0.56 for SBP, DBP and hypertension respectively (Supplementary Figure S1)). There was little consistency between the EWAS of Europeans and South Asians (goodness-of-fit R2=0 for SBP, DBP and hypertension (Supplementary Figure S1)).
Differentially methylated region analysis
DMR analyses were carried out to identify regions of DNA methylation that were associated with SBP, DBP and hypertension. As in the EWAS analyses, DMR analysis was conducted for the full trans-ancestry group and additionally for European and South Asian sub-groups. Trans-ancestry DMR analysis identified 395 regions (mapped to 325 annotated genes) for SBP, 237 regions (mapped to 160 annotated genes) for DBP and 0 DMRs for hypertension using a FDR-corrected p-value<0.05 (Supplementary Tables S19-S24). Twelve DMRs annotated to genes were common between trans-ancestry SBP and transancestry DBP at FDR p-value<0.05 (Table 3). In Europeans, 313 DMRS for SBP, 58 for DBP and 0 for hypertension were identified using a FDR-corrected p-value<0.05. Nine DMRs annotated to genes were common between European SBP and European DBP at FDR p-value<0.05 (Table 3). In South Asians, 60 DMRs for SBP, 88 DBP and 0 hypertension DMRs were identified using a FDR-corrected p-value<0.05. Thirty-five DMRs annotated to genes were common between South Asian SBP and South Asian DBP at FDR p-value<0.05 (Table 3).
|Overlap for trans-ancestry||Overlap for Europeans||Overlap for South Asians|
Table 3: Overlap of DMRs between fully adjusted DMR analysis of SBP and DBP for trans-ancestry, Europeans and South Asians respectively.
Known genetic variants in DMRs
Two of the twelve DMRs common in trans-ancestry SBP and DBP analysis, TFAP2D and HLX, contain SNPs previously associated with blood pressure [29,30]. The genetic variant in the PDZD2 DMR was associated with myocardial infarction . HLX was also identified in DMR analysis of SBP and DBP in Europeans (Table 6). PRDM16 was found to be common to SBP and DBP in Europeans its genetic variants were previously associated with dilated cardiomyopathy . PRKAG2 was found to be common to SBP and DBP in South Asians, and has been previously associated with hypertrophic cardiomyopathy  and chronic kidney disease . Other identified genes previously reported in association to CVD-related GWAS were ELL2 (GWAS of BP , insulin  and glucose ), TRPS1 (GWAS of BP ), PID1 (GWAS of stroke , lung function  and chronic obstructive pulmonary disease ) and WIPF1 (GWAS of resting heart rate ).
We investigated the association between SBP, DBP and hypertension and DNA methylation measured in peripheral blood of European and South Asian men combined and then across individual ethnicities using the HM450 BeadChip array. In the trans-ancestry fully adjusted EWAS, we found DBP was associated with methylation at one CpG site (cg07598370 near OR5AP2) at FDR-corrected p-value<0.05. The genetic variant near the olfactory receptor, family 5, subfamily AP, member 2 (OR5AP2) has been reported associated with hematological phenotypes  and olfactory receptors regulate blood pressure by the kidney via renal expression . SBP and hypertension were not associated with DNA methylation after adjustment for confounders, estimated cell counts and PCs. EWAS were also conducted in European and South Asian groups separately. In Europeans fully adjusted EWAS, seven CpG sites were associated with DBP with FDR-corrected p-value<0.05. In South Asians fully adjusted EWAS, SBP, DBP and hypertension were not associated with DNA methylation after multiple testing corrections. Several of the initial associations observed in the unadjusted model of DBP were noted to be documented loci responsive to tobacco smoking CpG sites, for example, (cg05575921 , cg12803068 , cg03636183 , cg22132788  and cg09935388 ), hence their attenuation on adjustment for smoking was predictable. This highlights the capacity of DNA methylation to index exposure to relevant risk factors.
We found seven CpG sites associated with DBP in Europeans. One CpG site (cg04751533) was in the gene body of AFAP (actin filament-associated protein), AFAP an action-binding and crosslinking protein is enriched in SRC and phorbol ester induced podosomes . Podosomes are specialized plasma-membrane actin-based microdomains and have been suggested to play a role in arterial vessel remodeling . C12orf44 (cg00006122) encodes autophagy-related protein 101 that is also known as ATG101. ATG101 is a newly found important autophagy related protein likewise many other proteins. Autophagy plays a key role in pulmonary vascular remodeling via regulation of apoptosis and hyperproliferation of pulmonary arterial endothelial cells . ATG101 is an essential gene for the initiation of autophagy and may be involved in endothelial cell growth through regulation of autophagy in pulmonary hypertension . CEBPD (cg16241714) encodes CCAAT/enhancer binding protein delta. CEBPD is involved in regulation of apoptosis and cell proliferation and probably acts as tumour suppressor . CD59 (CD59, molecule (CD59 blood group); cg06355652) has been associated with paroxysmal nocturnal hemogloinuria . Blood coagulation and cell surface receptor signalling are among its biological pathways [51,52]. In apolipoprotein E knockout mice, CD59, slowed down the development of atherosclerotic vulnerable plaque and retarded the progress of atherosclerosis . Component of the sarcomere protein, MYPN (myopalladin; cg18901940), chains together nebulin (skeletal muscle) and nebulette (cardiac muscle) to alpha-actinin, at the Z lines . It has been associated with dilated cardiomyopathy [55,56].
A large EWAS study of SBP and DBP was recently conducted in amongst individuals of European, African American and Hispanic/Latino ancestry (n=17,101) . The study reported 31 associations of DNA methylation with SBP or DBP. Of these, cg19693031 (near TXNIP) and cg18120259 (in gene body LOC100132354), found within the top hundred CpG sites of the trans-ancestry EWAS of SBP in our analyses, were among the 31 previously reported associations . The direction of effect was consistent with the previous analysis and the magnitude of association was slightly larger in our study. The above reported cg18120259 was also the top CpG site in our trans-ancestry EWAS of hypertension, the direction of effect was consistent but the magnitude of association was stronger in our study. One additional CpG site, cg06690548 (in gene body SLC7A11, among top 100 CpG sites in our study), was found in common between the trans-ancestry EWAS of hypertension and the previous study with same direction of effect and the magnitude of association was stronger in the current study. These three CpG sites were also replicated previously . Although the current study is smaller in size, we found evidence of some overlap between our results and a recent EWAS of blood pressure . Our study is the first blood pressure EWAS to our knowledge that has included South Asians, offering the chance to compare results from European and South Asian individuals.
There was no consistency in the magnitude and direction of associations comparing Europeans to South Asians, suggesting that peripheral blood DNA methylation patterns may differ between Europeans and South Asians in relation to blood pressure. This may reflect the fact that DNA methylation could index exposure to a different suite of risk factors in the two ethnic groups, or that different mechanisms contribute to the pathogenesis of hypertension and its related phenotypic traits. Of note, the South Asian participants in the SABRE study are first generation migrants, arriving in the UK as young adults. The potential early life and developmental antecedents of hypertension and blood pressure will therefore be considerably different between the two ethnic groups. This may explain to some extent the lack of consistency in methylation variable loci observed between the two ethnic groups. However, the study has limited statistical power due to the relatively modest sample size for these analyses.
Where associations were observed, the effect sizes were modest in size between cases and controls at the identified CpG sites. Such differences are unlikely to have profound biological consequences but may in turn exert a polygenic-like effect, altering disease risk or trait characteristics by small amounts. Further work is required to understand the functional consequences of such subtle shifts in DNA methylation.
We carried out DMR analysis and identified a large number of DMRs for SBP and DBP in the trans-ancestry, European and South Asian sub groups. The analyses found support for some of the DMRs for CVD related traits in the literature including BP , myocardial infarction , dilated cardiomyopathy  and stroke .
The strengths of this study include the study design where participants were from two different ethnicities; European and South Asian that were analyzed collectively and individually using robust statistical methods. Additionally, the utilization of HM450 arrays provided good coverage of the genome in terms of known annotated genes (although in total only covers <2% of all CpGs). The relatively modest sample size and the utilization of only male participants are among the limitations of this work.
In conclusion, we identified associations between methylation and DBP across trans-ancestry and European-specific analyses. Lack of associations identified in South Asian specific analyses indicates that the associations between methylation and blood pressure may be different between European and South Asian populations.
The SABRE study was funded at baseline by the Medical Research Council, Diabetes UK, and the British Heart Foundation. At follow-up the study was funded by the Wellcome Trust and the British Heart Foundation. Methylation analysis in the SABRE cohort was supported by a Wellcome Trust Enhancement grant (082464/Z/07/C). AH and NC received support from the Medical Research Council (Programme Code MC_UU_12019/1) and the National Institute for Health Research University College London Hospitals Biomedical Research Centre. The SABRE study team also acknowledges the support of the National Institute of Health Research Clinical Research Network (NIHR CRN). NK, HRE, TRG and CLR are supported by the Medical Research Council Integrative Epidemiology Unit at the University of Bristol (MC_UU_12013/2 and MC_UU_12013/8).
The authors declare that they have no competing interests. TRG receives funding from GlaxoSmithKline, Biogen and Sanofi for research unrelated to the work presented here. The SABRE study team also acknowledges the support of the National Institute of Health Research Clinical Research Network (NIHR CRN).
All participants gave written informed consent. Approval for the baseline study was obtained from Ealing, Hounslow and Spelthorne, Parkside and University College London research ethnics committees.
Data used in these analyses are available from SABRE on request.
NK, HRE, NC and CLR designed the study. NK conducted the statistical analysis and wrote the manuscript with input from all authors. Correspondence and material requests should be addressed to NK ([email protected]). All authors read and approved the final manuscript.
All Published work is licensed under a Creative Commons Attribution 4.0 International License
Copyright © 2019 All rights reserved. iMedPub LTD Last revised : March 23, 2019