Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Immunol., 11 February 2026

Sec. Cytokines and Soluble Mediators in Immunity

Volume 17 - 2026 | https://doi.org/10.3389/fimmu.2026.1671540

This article is part of the Research TopicDiscussing the implication of DNA methylation in human diseasesView all 5 articles

Integrative analysis of DNA methylation and inflammatory protein biomarkers in hypertension

Min Zhu&#x;Min ZhuJin Zhang&#x;Jin ZhangJing MaJing MaXiaofeng TangXiaofeng TangYan Wang*Yan Wang*Dingliang ZhuDingliang Zhu
  • Department of Cardiovascular Medicine, Research Center for Hypertension Management and Prevention in Community, Shanghai Key Laboratory of Hypertension, Shanghai Institute of Hypertension, State Key Laboratory of Medical Genomics, Ruijin Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, China

Background and objective: To investigate the relationship between inflammation and hypertension by comparing epigenetic and serum biomarkers of inflammation and their association with target organ damage (TOD).

Methods: The epigenome-wide methylation profiles of peripheral leukocyte DNA from 176 patients with hypertension were analyzed using Illumina Infinium Methylation EPIC BeadChips. The commercially available Olink® Target 96 Inflammation panels were utilized to evaluate markers associated with inflammation. We identified CpG-protein association protein quantitative trait methylation loci (pQTMs) using mixed linear regression, adjusting for potential confounders. The possible roles of the discovered pQTMs were ascertained by Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analyses. The association between pQTMs and TOD was also analyzed.

Results: In our analysis, we found 771 significant associations across 11 biomarkers, with a false discovery rate (FDR) of less than 0.05. Lambda estimates ranged from 0.847 to 1.202. Among these, 39 pQTMs showed an association with urine albumin-creatinine ratio, while 82 pQTMs showed an association with carotid intima-media thickness.

Conclusion: Our findings contribute to the understanding of inflammatory biomarkers associated with alterations in DNA methylation in hypertension. However, these results will need to be validated in future studies.

Introduction

Hypertension ranks as the leading cause of cardiovascular disease and premature death, affecting over one billion adults globally (1). It is thought to be responsible for roughly 54% of strokes and 47% of ischemic heart diseases (2). A number of diseases are largely caused by low-grade inflammation. There is growing evidence that the immune system and an activated inflammatory response play crucial roles in the development of hypertension and are linked to its complications, including myocardial infarction, hemorrhagic stroke, and renal damage (3). It has been demonstrated that the cytokines IL-17, IFN-γ, TNFα, and IL-6 can cause damage to the blood vessels and kidneys, as well as increase sodium retention and systemic vascular resistance (4).

Cytosine-guanine (CpG) dinucleotides are among the most extensively profiled blood-based epigenetic modifications on a large scale. Their methylation (DNAm) can result from exposure to risk factors, such as obesity (5) or smoking (6), impacting regulatory regions of the genome. Additionally, diseases like cancer can lead to global changes in CpG methylation (7). Methylome-wide association studies of protein levels can quantify associations between DNAm at CpGs and protein levels, which are referred to as protein quantitative trait methylation loci (pQTMs) (8). Epigenome-wide association studies (EWAS) have been conducted by several researches on the expression levels of a select group of inflammatory proteins, such as interleukins, IFN-γ, TGF-β, and TNF (9, 10). A meta-analysis of EWAS of CRP revealed multiple inflammation-related CpG sites linked to cardiometabolic phenotypes and coronary heart disease (11). Recently, several studies have focused on the combined effects of multi-omics on clinical phenotypes. Zaghlool et al. performed an EWAS of 1,123 proteins, which identified networks of chronic low-grade inflammatory biomarkers that overlapped with transcriptomics, metabolomics, and clinical endpoints (12). Hillary et al. conducted EWAS on 92 neurological proteins and identified 26 methylation sites associated with nine proteins related to neurological, immunological, and extracellular matrix metabolic pathways (13).

However, research focusing on the influence of inflammatory proteins on changes in DNA methylation patterns in patients with hypertension is scarce. To identify the pathogenic mechanisms underlying inflammation in hypertension, we examined the relationship between serum cytokines detected via Olink Target 96 Inflammation panels and DNA methylation patterns in the peripheral leukocytes of patients with hypertension. Furthermore, we explored their effects on target organ damage (TOD) in hypertension (Figure 1).

Figure 1
Flowchart showing a study involving 176 hypertension patients. Peripheral leukocytes DNA undergoes Methylation EPIC 850K BeadChips analysis on 725,736 CpGs. Plasma is analyzed with the OLINK Inflammatory panel for 74 proteins. Results in pEWAS with 771 pQTMs at FDR<0.05. Of these, 39 pQTMs are associated with UACR, and 82 pQTMs with CIMT.

Figure 1. Study design.

Materials and methods

Study population

The participants in the present study were outpatients with hypertension who were identified in Xinzhuang Community Hospital between May 2012 and December 2015. Of the 176 enrolled patients with hypertension, 110 were males. The average age was 66.2 years, ranging from 50–80 years. Participants were included if they had hypertension, defined as BP ≥ 140/90 mmHg, or were taking antihypertensive medications. Those with previous cardiovascular diseases, including stroke, myocardial infarction, heart failure, or dementia, were excluded from the study. Trained clinicians conducted clinical interviews to gather information on the participants’ history of current smoking, drinking, hypertension, diabetes, and recent drug prescriptions, all within a two-week period. General physical examinations, which included assessments of body height and weight, were conducted by trained nurses. Patients were asked to fast and remove shoes or coats during these examinations. The details of clinical information collection have been described previously (14).

All participants provided written informed consent, and the study protocol was approved by the Ruijin Hospital ethics committee.

Genome-wide DNA methylation profiling analysis

DNA was isolated from peripheral leukocytes using a DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany). The purity and concentration of DNA were estimated using a Nanodrop 2000 spectrophotometer (Thermo). Subsequently, 500 ng of DNA was bisulfite-converted using EZ DNA Methylation Kits (Zymo Research, USA).

Genome-wide DNA methylation profiling analysis was performed using Illumina Infinium Methylation EPIC 850 K BeadChips (Illumina, San Diego, CA, USA). The data were mainly analyzed using the lmer package in R. Probes were filtered out if they failed to hybridize with more than three beads in at least 5% of the samples, or if the detection p-value was > 0.01. Additionally, non-CpG probes, SNP-related probes, multi-hit probes, and probes located on the X and Y chromosomes were excluded. Ultimately, 725,736 probes were retained for subsequent analysis. β value (ranging from 0 to 1) was used to represent the methylation level at CpG sites, according to the formula β = m/(m + u + 100), where m and u represent methylated and unmethylated probe intensities, respectively). Subsequently, the β value was normalized using BMIQ.

Proteomic analysis

Olink® Target 96 Inflammation panels were used to measure plasma inflammatory markers (Olink, Uppsala, Sweden). These panels analyze 92 biomarkers associated with inflammation using a highly sensitive and specific proximity extension assay technology (12). In short, each target protein is matched with a unique complementary DNA barcode and recognized by double antibodies. The Biomark HD (Fluidigm, South San Francisco, CA, USA) is a high-throughput microfluidic real-time PCR device that is used to quantify these barcodes. The resulting cycle threshold (Ct) data were processed for quality control and normalization using internal and external controls. The final assay read-out was given as the normalized protein expression values, which were expressed in arbitrary log2-scale units and corresponded to higher protein levels.

TOD estimation

Serum creatinine levels were measured at the central laboratory of Ruijin Hospital using fasting venous blood. The Modification of Diet in Renal Disease equation was used to determine the estimated glomerular filtration rate (eGFR). Urine chemistry tests were also conducted the same day, and one spot urine specimen was collected in the morning. A Siemens ADVIA 2400 automated chemistry analyzer (Siemens Healthcare Diagnostics, USA) and commercially available kits (Shanghai Kehua Bio-Engineering, China) were used to measure the urine albumin-creatinine ratio (UACR).

Following standard protocols, echocardiographic data were obtained using a commercially available device (E9; GE Healthcare, Milwaukee, WI, USA). LV diameters and septal and posterior wall thicknesses were measured in accordance with the guidelines of the American Society of Echocardiography using 2-dimensionally guided M-mode tracings, with a recording speed of 50–100 cm/s. LV mass was calculated from M-mode echocardiograms using the formula described by Devereux et al. The LV mass index (LVMI) in g/m2 was indexed from the LV mass to the body surface area.

The intima-media thickness (IMT) of the left and right common carotid arteries was measured using high-resolution ultrasound (E9; GE Healthcare, Milwaukee, WI, USA). The distance between the leading edge of the internal lumen and the media adventitia echo was defined as the IMT. Using the larger IMT measurement found in the bilateral common carotid artery, the final IMT value was calculated.

After the participant had rested in the supine position for at least five minutes, Cf-PWV was measured. Using a validated device (SphygmoCor XCEL; AtCor Medical Pty Ltd, Sydney, Australia), the time delay between the right femoral artery and the right common carotid artery’s rapid upstroke was measured. Using a tape measure, the distance between two recording sites on the body’s surface was calculated to estimate the pulse wave’s travel distance. The distance in meters divided by the time delay in seconds was automatically used to compute Cf-PWV.

Utilizing the Mini-Mental State Examination (MMSE), we conducted an oral examination and evaluated cognitive functioning.

Statistical analysis

To identify CpG-protein associations (pQTMs), a mixed linear regression as conducted in R (version 4.0.0), adjusting for age, sex, smoking, slide, array, and estimated cell counts. Cellular heterogeneity due to variations in cell population proportions in whole blood (15) was controlled using cell proportion estimates based on the Houseman algorithm (16). Estimated cell counts were generated for B-lymphocytes, CD4+ and CD8+ T-lymphocytes, natural killer cells, granulocytes, and monocytes using the “estimate cell counts” function in the R package minfi (17, 18). The regression coefficient was estimated using the REML fitting algorithm, and the p-value was estimated using the Satterthwaite method. The association of CpG/protein with the TOD of hypertension (LVMI, IMT, Cf-PWV, UACR, eGFR, and MMSE score) was estimated via linear regression, adjusting for age, sex, and mean arterial BP. To control for multiple testing, we applied false discovery rate (FDR) correction using the Benjamini-Hochberg method.

Results

Key clinical parameters in patients with hypertension

The average BP was 137.2/71.4 mmHg. The average values for LVMI, IMT, Cf-PWV, UACR, eGFR, and MMSE score were 101.9 g/m2, 0.8 mm, 9.4 m/s, 25.3 mg/g, 99.3 ml/min/1.73 m2, and 27.8, respectively (Supplementary Table 1).

EWAS of inflammatory biomarker levels

In total, 74 of the 92 biomarkers passed the quality control and were used in the analysis. Through EWAS, we identified 771 significant associations across the abundance levels of 11 biomarkers (FDR<0.05). The lambda estimates for these associations ranged from 0.847 to 1.202 (Figure 2, Supplementary Figure 1). Among these associations, 557 pQTMs were annotated to known genes. Specifically, 74 were located within the TSS200 region, 81 within TSS1500, and 74 within the 5’UTR (Supplementary Table 2). Notably, the methylation levels at 197 loci showed a positive association with protein abundance, while 360 associations were negative. Twenty regions contained multiple biomarker-associated CpG sites, including those associated with FGF5, FGF23, NTF3, IL10RA, FGF5, and CD274. The degree of methylation of the 13 genes was related to multiple proteins, whereas the level of methylation at loci within PTPRN2 was associated with the abundance of CCL7, FGF23, and FGF5. The top, differentially methylated positions associated with each biomarker in patients with hypertension are shown in Table 1.

Figure 2
Manhattan plot and Q-Q plot for TNFRSF9 gene analysis. The Manhattan plot shows -log10 p-values by chromosome, with significant markers labeled. The Q-Q plot compares observed versus expected p-values, indicating data distribution with slight deviation from the null hypothesis line.

Figure 2. Manhattan plot from the epigenome-wide association studie in 176 patients with hypertension. 725,736 CpG sites were tested for association with the level of TNFRSF9. The black dotted horizontal line represents the FDR< 0.05. QQ plot showing observed vs. expected − log10 (P values) for association at all CpG sites.

Table 1
www.frontiersin.org

Table 1. Top 11 differentially methylated positions associated with each biomarker in patients with hypertension.

pQTMs and TOD in hypertension

Considering that inflammatory factors are involved in the damage process of target organs in hypertension, we further evaluated the association of the 771 pQTMs with LVMI, IMT, Cf-PWV, UACR, eGFR, and MMSE scores using Pearson’s correlation and multiple linear regression analysis. Thirty-nine pQTMs were associated with UACR, even after adjusting for risk factors, including age, sex, BMI, current smoking, SBP, and LDL (Table 2). These pQTMs corresponded to four proteins: TNFRSF9, CCL7, FGF23, and IL6, and the expression levels of these five proteins were positively correlated (Figure 3). The methylation levels of most loci were negatively correlated with the abundance of four proteins, except cg04092513, cg10477088, cg10478101, cg15871214, cg19756821, and cg21184859. Of the 39 pQTMs, 33 were related to TNFRSR9. Additionally, 82 pQTMs were found to be associated with CIMT after adjusting for risk factors. Among these, the methylation levels of 59 CpG sites were negatively correlated with protein abundance. These 82 pQTMs encompassed three unique proteins, with NTF3 being associated with 77 CpG sites (Supplementary Table 3, Supplementary Figure 2). No pQTMs were associated with LVMI, Cf-PWV, eGFR, or MMSE scores.

Table 2
www.frontiersin.org

Table 2. Association of the methylation levels of pQTMs with ACR.

Figure 3
Correlation matrix heatmap depicting relationships between multiple variables, with a color gradient from blue to red indicating correlation values ranging from 1 to -1. Variables are labeled along both axes.

Figure 3. Pairwise correlation matrix for 4 proteins abundance with 39 methylation level of pQTMs assocated with ACR.

Functional enrichment analysis

Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the top 771 DMPs correlating with inflammatory factors were performed using the R package clusterProfiler. Under the condition of FDR less than 0.05, the differentially methylated genes were found to be involved in two biological processes: organ development and morphogenesis, as well as the Wnt signaling pathway. Regarding cellular components, enriched terms included the cell body membrane, fibrillar center, glutamatergic synapses, and neuronal cell body membrane (Figure 4). However, no item was found to be enriched in the KEGG pathway analysis with a FDR of less than 0.05.

Figure 4
Scatter plot showing GO enrichment analysis, indicating gene ratios against GO terms. Dots vary in size and color, representing counts and P-values respectively. Two panels display biological processes (BP) and cellular components (CC). Enrichment terms include embryonic development, Wnt signaling pathways, and cell membranes.

Figure 4. GO pathway enrichment analysis was conducted for the 771 DMPs correlated with inflammatory biomarkers using the R package clusterProfiler. The enrichment analysis was performed with a false discovery rate (FDR) < 0.05 as the threshold for statistical significance. The figure displays the top enriched terms in biological process (BP) and cellular component (CC). The size of each bubble represents the number of DMPs enriched in the corresponding term, and the color gradient indicates the FDR value (from blue to red, reflecting increasing enrichment significance).

Discussion

In the present study, we performed the first epigenome-wide analysis of DNA methylation changes associated with the levels of multiple inflammatory proteins in patients with hypertension. Our analysis revealed a total of 771 significant associations with the abundance levels of 11 biomarkers. Notably, among these associations, 39 pQTMs were found to be associated with UACR, while 82 pQTMs showed associations with CIMT. Most associations between DNA methylation and protein levels are related to age (19), blood cell composition (20), and rheumatoid arthritis (21). Additionally, many CpG loci that appear in pQTMs are associated with other plasma proteins, including IGSF3, F11, and PGAM1 (8).

Aberrant inflammatory responses result in tissue damage and are central to the pathogenesis of multiple diseases, including diabetes (22), cardiovascular disorders (23), and cognitive dysfunction (24). The chemokines involved in the pathogenesis of hypertension include monocyte chemoattractant protein-1, MCP-1, CCL2, interferon-inducible protein, interleukin-8, RANTES, fractalkine, along with their receptors. However, the mechanisms involving chemokines and their receptors in hypertension are intricate and not fully elucidated (25). Association studies between DNA methylation and chronic low-grade inflammation have primarily focused on scattered proteins, with CRP being a common target. Wielscher et al. (26) carried out a multi-ethnic EWAS to characterize the DNA methylation signature associated with chronic low-grade inflammation, as measured by CRP. They identified 1,511 independent differentially methylated loci associated with CRP levels. Importantly, they found that an activated CpG signature significantly increased the risk of cardiometabolic diseases and COPD. However, few studies have detected an association between genome-wide methylation alterations and multiple inflammatory proteins, particularly in patients with hypertension. Hillary et al. (27) performed an EWAS on the levels of 70 plasma-derived inflammatory protein biomarkers in healthy, older adults. They identified three CpG sites spread across three proteins (CCL11, IL18R1, and CXCL9) and demonstrated putative causal relationships between IL18R1 and inflammatory bowel disease. No overlap was observed between our results and those reported by Hillary et al. Such inconsistencies may be induced by different disease backgrounds that influence the interaction between immune factors and DNA methylation.

Most of the 39 pQTMs related to UACR were linked to TNFRSF9, with 34 out of 39 being specifically related to this protein. TNFRSF9, also known as CD137 or 4-1BB, belongs to a group of costimulatory immune receptors and is a member of the tumor necrosis factor receptor (TNF-R) superfamily. It is preferentially found on activated T cells, regulatory T cells (Tregs), and innate immune cells (28). TNFRSF9 is implicated in various immunomodulatory responses related to kidney disease, as well as biological processes like cell survival, proliferation, and death. In a study by Nano et al. (29), TNFRSF9 was associated with baseline UACR and an increased risk of incident chronic kidney disease, independent of other cardiorenal risk factors. A similar study by Niewczas et al. found a kidney risk inflammatory signature consisting of 17 inflammatory proteins enriched for TNF-R superfamily members in participants from three different cohorts with type 1 and type 2 diabetes (30), which significantly increased the 10-year risk of end-stage renal disease. The level of TNFRSF9 was positively related to UACR (β = 0.4, p = 3.73e-08), consistent with the previous findings. Among the 34 CpGs associated with TNFRSF9, 25 loci were annotated as known genes, wherein fat atypical cadherin 1 (FAT1) was associated with UACR levels in genome-wide association studies (31) and family studies (32). FAT1 is a calcium ion-dependent adhesion protein that belongs to the fat family. Notably, FAT1 regulates cell migration and growth through specific protein–protein interactions in its cytoplasmic tail (33). Moreover, FAT1 is abundantly expressed in the glomerular podocytes; specific deletion of FAT1 in mouse podocytes can induce abnormal development of the glomerular filtration barrier, reduce cell adhesion and migration of fibroblasts and podocytes, and lead to the disappearance of foot processes and slits (31). Renal biopsies of patients with membranous nephropathy after hematopoietic stem cell transplantation also revealed granular staining for FAT1 along the glomerular basement membrane (34). In a study on the pan-cancer landscape of T-cell heterogeneity, Zheng et al. reported that FAT1 mutations were positively correlated with TNFRSF9+ Treg cell frequencies. However, the influence of TNFRSF9 on the methylation of FAT1 and its combined effect on UACR requires further study.

Our study has several strengths and limitations. To the best of our knowledge, this is the first study to investigate the association between peripheral blood leukocyte DNA methylation and the levels of multiple inflammatory proteins in patients with hypertension using an 850k chip and an Olink inflammation panel. These results should also be interpreted with consideration of their limitations. First, owing to the relatively small sample size, the risk of false-positive results should be considered, necessitating further replication studies with larger sample sizes. However, to mitigate the impact of this limitation, we applied the Benjamini–Hochberg method to control for FDR. Second, given the descriptive nature of the study, it cannot establish a definitive cause-and-effect relationship. Additional prospective studies are therefore required to confirm our findings. Third, our study included patients with hypertension, which may hinder the applicability of our conclusions to the general population. Nevertheless, previous research has indicated that the disease background can influence methylation changes induced by inflammatory proteins (35, 36). Therefore, targeting DNA methylation changes related to inflammatory protein levels, especially in hypertension, may reveal specific underlying mechanisms. Additionally, we found that pQTMs in hypertension were also related to UACR and CIMT, the surrogate markers of TOD, which offered insights into the underlying function.

In conclusion, we integrated DNA methylation data and inflammatory protein abundance in patients with hypertension and identified 771 novel pQTMs, some of which were related to UACR and CIMT. Our findings contribute to the understanding of the inflammatory factors associated with DNA methylation alterations in hypertension. However, our results will need to be validated in future studies.

Data availability statement

The data presented in the study are deposited in NCBI's Gene Expression Omnibus (GEO) repository, accession number GSE314180 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE314180).

Ethics statement

The studies involving humans were approved by Ruijin Hospital ethics committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

MZ: Methodology, Writing – original draft, Formal analysis. YW: Writing – review & editing, Data curation, Conceptualization, Writing – original draft, Funding acquisition. JZ: Writing – original draft, Data curation. JM: Writing – original draft, Methodology. XT: Writing – original draft, Methodology. DZ: Writing – original draft, Conceptualization.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This project was supported by grants from National Research Center for Translational Medicine (Shanghai) grant NRCTM(SH)-2021-05, Shanghai Municipal Health Bureau Foundation (202540024), the National Natural Science Foundation of China (81670213). The funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript. The datasets analyzed during the current study are not publicly available due to the consent obtained, which does not grant permission for deposition but are available from the corresponding author on reasonable request.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

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

Supplementary material

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

References

1. Mills KT, Stefanescu A, and He J. The global epidemiology of hypertension. Nat Rev Nephrol. (2020) 16:223–37. doi: 10.1038/s41581-019-0244-2

PubMed Abstract | Crossref Full Text | Google Scholar

2. Lawes CM, Vander Hoorn S, Rodgers A, and International Society of, H. Global burden of blood-pressure-related disease, 2001. Lancet. (2008) 371:1513–8. doi: 10.1016/S0140-6736(08)60655-8

PubMed Abstract | Crossref Full Text | Google Scholar

3. Zhang Z, Zhao L, Zhou X, Meng X, and Zhou X. Role of inflammation, immunity, and oxidative stress in hypertension: New insights and potential therapeutic targets. Front Immunol. (2022) 13:1098725. doi: 10.3389/fimmu.2022.1098725

PubMed Abstract | Crossref Full Text | Google Scholar

4. McMaster WG, Kirabo A, Madhur MS, and Harrison DG. Inflammation, immunity, and hypertensive end-organ damage. Circ Res. (2015) 116:1022–33. doi: 10.1161/CIRCRESAHA.116.303697

PubMed Abstract | Crossref Full Text | Google Scholar

5. Sharma NK, Comeau ME, Montoya D, Pellegrini M, Howard TD, Langefeld CD, et al. Integrative analysis of glucometabolic traits, adipose tissue DNA methylation, and gene expression identifies epigenetic regulatory mechanisms of insulin resistance and obesity in African Americans. Diabetes. (2020) 69:2779–93. doi: 10.2337/db20-0117

PubMed Abstract | Crossref Full Text | Google Scholar

6. Tsai PC, Glastonbury CA, Eliot MN, Bollepalli S, Yet I, Castillo-Fernandez JE, et al. Smoking induces coordinated DNA methylation and gene expression changes in adipose tissue with consequences for metabolic health. Clin Epigenet. (2018) 10:126. doi: 10.1186/s13148-018-0558-0

PubMed Abstract | Crossref Full Text | Google Scholar

7. Skvortsova K, Stirzaker C, and Taberlay P. The DNA methylation landscape in cancer. Essays Biochem. (2019) 63:797–811. doi: 10.1042/EBC20190037

PubMed Abstract | Crossref Full Text | Google Scholar

8. Gadd DA, Hillary RF, McCartney DL, Shi L, Stolicyn A, Robertson NA, et al. Integrated methylome and phenome study of the circulating proteome reveals markers pertinent to brain health. Nat Commun. (2022) 13:4670. doi: 10.1038/s41467-022-32319-8

PubMed Abstract | Crossref Full Text | Google Scholar

9. Verschoor CP, McEwen LM, Kohli V, Wolfson C, Bowdish DM, Raina P, et al. The relation between DNA methylation patterns and serum cytokine levels in community-dwelling adults: a preliminary study. BMC Genet. (2017) 18:57. doi: 10.1186/s12863-017-0525-3

PubMed Abstract | Crossref Full Text | Google Scholar

10. Liang L, Willis-Owen SAG, Laprise C, Wong KCC, Davies GA, Hudson TJ, et al. An epigenome-wide association study of total serum immunoglobulin E concentration. Nature. (2015) 520:670–4. doi: 10.1038/nature14125

PubMed Abstract | Crossref Full Text | Google Scholar

11. Ligthart S, Marzi C, Aslibekyan S, Mendelson MM, Conneely KN, Tanaka T, et al. DNA methylation signatures of chronic low-grade inflammation are associated with complex diseases. Genome Biol. (2016) 17:255. doi: 10.1186/s13059-016-1119-5

PubMed Abstract | Crossref Full Text | Google Scholar

12. Zaghlool SB, Kuhnel B, Elhadad MA, Kader S, Halama A, Thareja G, et al. Epigenetics meets proteomics in an epigenome-wide association study with circulating blood plasma protein traits. Nat Commun. (2020) 11:15. doi: 10.1038/s41467-019-13831-w

PubMed Abstract | Crossref Full Text | Google Scholar

13. Hillary RF, McCartney DL, Harris SE, Stevenson AJ, Seeboth A, Zhang Q, et al. Genome and epigenome wide studies of neurological protein biomarkers in the Lothian Birth Cohort 1936. Nat Commun. (2019) 10:3160.

PubMed Abstract | Google Scholar

14. Zhu M, An D, Zhang J, Tang X, Wang Y, and Zhu D. Genome-wide analysis of DNA methylation and its relationship with serum homocysteine levels in patients with hypertension. J Hypertens. (2023) 41:1626–33. doi: 10.1097/HJH.0000000000003515

PubMed Abstract | Crossref Full Text | Google Scholar

15. Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlen SE, Greco D, et al. Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PloS One. (2012) 7:e41361. doi: 10.1371/journal.pone.0041361

PubMed Abstract | Crossref Full Text | Google Scholar

16. Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinf. (2012) 13:86. doi: 10.1186/1471-2105-13-86

PubMed Abstract | Crossref Full Text | Google Scholar

17. Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. (2014) 30:1363–9. doi: 10.1093/bioinformatics/btu049

PubMed Abstract | Crossref Full Text | Google Scholar

18. Jaffe AE and Irizarry RA. Accounting for cellular heterogeneity is critical in epigenome-wide association studies. Genome Biol. (2014) 15:R31. doi: 10.1186/gb-2014-15-2-r31

PubMed Abstract | Crossref Full Text | Google Scholar

19. Mulder RH, Neumann A, Cecil CAM, Walton E, Houtepen LC, Simpkin AJ, et al. Epigenome-wide change and variation in DNA methylation in childhood: trajectories from birth to late adolescence. Hum Mol Genet. (2021) 30:119–34. doi: 10.1093/hmg/ddaa280

PubMed Abstract | Crossref Full Text | Google Scholar

20. Islam SA, Goodman SJ, MacIsaac JL, Obradovic J, Barr RG, Boyce WT, et al. Integration of DNA methylation patterns and genetic variation in human pediatric tissues help inform EWAS design and interpretation. Epigenet Chromatin. (2019) 12:1. doi: 10.1186/s13072-018-0245-6

PubMed Abstract | Crossref Full Text | Google Scholar

21. Liu Y, Aryee MJ, Padyukov L, Fallin MD, Hesselberg E, Runarsson A, et al. Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis. Nat Biotechnol. (2013) 31:142–7. doi: 10.1038/nbt.2487

PubMed Abstract | Crossref Full Text | Google Scholar

22. Schmidt MI, Duncan BB, Sharrett AR, Lindberg G, Savage PJ, Offenbacher S, et al. Markers of inflammation and prediction of diabetes mellitus in adults (Atherosclerosis Risk in Communities study): a cohort study. Lancet. (1999) 353:1649–52. doi: 10.1016/S0140-6736(99)01046-6

PubMed Abstract | Crossref Full Text | Google Scholar

23. Danesh J, Wheeler JG, Hirschfield GM, Eda S, Eiriksdottir G, Rumley A, et al. C-reactive protein and other circulating markers of inflammation in the prediction of coronary heart disease. N Engl J Med. (2004) 350:1387–97. doi: 10.1056/NEJMoa032804

PubMed Abstract | Crossref Full Text | Google Scholar

24. Kim J, Yoon S, Lee S, Hong H, Ha E, Joo Y, et al. A double-hit of stress and low-grade inflammation on functional brain network mediates posttraumatic stress symptoms. Nat Commun. (2020) 11:1898. doi: 10.1038/s41467-020-15655-5

PubMed Abstract | Crossref Full Text | Google Scholar

25. Rudemiller NP and Crowley SD. The role of chemokines in hypertension and consequent target organ damage. Pharmacol Res. (2017) 119:404–11. doi: 10.1016/j.phrs.2017.02.026

PubMed Abstract | Crossref Full Text | Google Scholar

26. Wielscher M, Mandaviya PR, Kuehnel B, Joehanes R, Mustafa R, Robinson O, et al. DNA methylation signature of chronic low-grade inflammation and its role in cardio-respiratory diseases. Nat Commun. (2022) 13:2408. doi: 10.1038/s41467-022-29792-6

PubMed Abstract | Crossref Full Text | Google Scholar

27. Hillary RF, Trejo-Banos D, Kousathanas A, McCartney DL, Harris SE, Stevenson AJ, et al. Multi-method genome- and epigenome-wide studies of inflammatory protein levels in healthy older adults. Genome Med. (2020) 12:60. doi: 10.1186/s13073-020-00754-1

PubMed Abstract | Crossref Full Text | Google Scholar

28. Stoll A, Bruns H, Fuchs M, Volkl S, Nimmerjahn F, Kunz M, et al. CD137 (4-1BB) stimulation leads to metabolic and functional reprogramming of human monocytes/macrophages enhancing their tumoricidal activity. Leukemia. (2021) 35:3482–96. doi: 10.1038/s41375-021-01287-1

PubMed Abstract | Crossref Full Text | Google Scholar

29. Nano J, Schottker B, Lin JS, Huth C, Ghanbari M, Garcia PM, et al. Novel biomarkers of inflammation, kidney function and chronic kidney disease in the general population. Nephrol Dial Transplant. (2022) 37:1916–26. doi: 10.1093/ndt/gfab294

PubMed Abstract | Crossref Full Text | Google Scholar

30. Niewczas MA, Pavkov ME, Skupien J, Smiles A, Md Dom ZI, Wilson JM, et al. A signature of circulating inflammatory proteins and development of end-stage renal disease in diabetes. Nat Med. (2019) 25:805–13. doi: 10.1038/s41591-019-0415-5

PubMed Abstract | Crossref Full Text | Google Scholar

31. Lahrouchi N, George A, Ratbi I, Schneider R, Elalaoui SC, Moosa S, et al. Homozygous frameshift mutations in FAT1 cause a syndrome characterized by colobomatous-microphthalmia, ptosis, nephropathy and syndactyly. Nat Commun. (2019) 10:1180. doi: 10.1038/s41467-019-08547-w

PubMed Abstract | Crossref Full Text | Google Scholar

32. Gee HY, Sadowski CE, Aggarwal PK, Porath JD, Yakulov TA, Schueler M, et al. FAT1 mutations cause a glomerulotubular nephropathy. Nat Commun. (2016) 7:10822. doi: 10.1038/ncomms10822

PubMed Abstract | Crossref Full Text | Google Scholar

33. Sadeqzadeh E, de Bock CE, and Thorne RF. Sleeping giants: emerging roles for the fat cadherins in health and disease. Med Res Rev. (2014) 34:190–221. doi: 10.1002/med.21286

PubMed Abstract | Crossref Full Text | Google Scholar

34. Sethi S, Madden B, Casal Moura M, Nasr SH, Klomjit N, Gross L, et al. Hematopoietic stem cell transplant-Membranous nephropathy is associated with protocadherin FAT1. J Am Soc Nephrol. (2022) 33:1033–44. doi: 10.1681/ASN.2021111488

PubMed Abstract | Crossref Full Text | Google Scholar

35. Mori KM, McElroy JP, Weng DY, Chung S, Fadda P, Reisinger SA, et al. Lung mitochondrial DNA copy number, inflammatory biomarkers, gene transcription and gene methylation in vapers and smokers. EBioMedicine. (2022) 85:104301. doi: 10.1016/j.ebiom.2022.104301

PubMed Abstract | Crossref Full Text | Google Scholar

36. Zhao N, Ruan M, Koestler DC, Lu J, Salas LA, Kelsey KT, et al. Methylation-derived inflammatory measures and lung cancer risk and survival. Clin Epigenet. (2021) 13:222. doi: 10.1186/s13148-021-01214-2

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: DNA methylation, epigenome-wide analysis, hypertension, inflammatory biomarkers, target organ damage

Citation: Zhu M, Zhang J, Ma J, Tang X, Wang Y and Zhu D (2026) Integrative analysis of DNA methylation and inflammatory protein biomarkers in hypertension. Front. Immunol. 17:1671540. doi: 10.3389/fimmu.2026.1671540

Received: 23 July 2025; Accepted: 19 January 2026; Revised: 15 January 2026;
Published: 11 February 2026.

Edited by:

Viviana Barra, University of Palermo, Italy

Reviewed by:

Ying Chen, Kunming Medical University, China
Sujith Dassanayaka, University of Louisville, United States
Tiar Masykuroh Pratamawati, Swadaya Gunung Djati University, Indonesia

Copyright © 2026 Zhu, Zhang, Ma, Tang, Wang and Zhu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yan Wang, eWFud2FuZ3NoamRAMTYzLmNvbQ==

These authors have contributed equally to this work

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.