Skip to main content

Whole blood methylome-derived features to discriminate endocrine hypertension



Arterial hypertension represents a worldwide health burden and a major risk factor for cardiovascular morbidity and mortality. Hypertension can be primary (primary hypertension, PHT), or secondary to endocrine disorders (endocrine hypertension, EHT), such as Cushing's syndrome (CS), primary aldosteronism (PA), and pheochromocytoma/paraganglioma (PPGL). Diagnosis of EHT is currently based on hormone assays. Efficient detection remains challenging, but is crucial to properly orientate patients for diagnostic confirmation and specific treatment. More accurate biomarkers would help in the diagnostic pathway. We hypothesized that each type of endocrine hypertension could be associated with a specific blood DNA methylation signature, which could be used for disease discrimination. To identify such markers, we aimed at exploring the methylome profiles in a cohort of 255 patients with hypertension, either PHT (n = 42) or EHT (n = 213), and at identifying specific discriminating signatures using machine learning approaches.


Unsupervised classification of samples showed discrimination of PHT from EHT. CS patients clustered separately from all other patients, whereas PA and PPGL showed an overall overlap. Global methylation was decreased in the CS group compared to PHT. Supervised comparison with PHT identified differentially methylated CpG sites for each type of endocrine hypertension, showing a diffuse genomic location. Among the most differentially methylated genes, FKBP5 was identified in the CS group. Using four different machine learning methods—Lasso (Least Absolute Shrinkage and Selection Operator), Logistic Regression, Random Forest, and Support Vector Machine—predictive models for each type of endocrine hypertension were built on training cohorts (80% of samples for each hypertension type) and estimated on validation cohorts (20% of samples for each hypertension type). Balanced accuracies ranged from 0.55 to 0.74 for predicting EHT, 0.85 to 0.95 for predicting CS, 0.66 to 0.88 for predicting PA, and 0.70 to 0.83 for predicting PPGL.


The blood DNA methylome can discriminate endocrine hypertension, with methylation signatures for each type of endocrine disorder.


Arterial hypertension affects over a billion people worldwide, with an estimated overall prevalence in adults of around 30–45% [1, 2]. This is one of the major risk factors for multiple cardiovascular and renal disorders and represents one of the most preventable causes of morbidity and mortality. Around 5 to 10% of arterial hypertension cases are estimated to have secondary causes, most commonly due to parenchymal renal disease, renovascular hypertension, obstructive sleep apnea, or endocrine diseases, such as Cushing’s syndrome (CS), primary aldosteronism (PA), and pheochromocytoma/paraganglioma (PPGL) [3, 4]. Diagnostic screening for secondary hypertension is complex and expensive [5] and generally restricted to patients with clinical signs, including hypertension in young adults, sudden worsening of blood pressure in normotensive subjects, and drug-resistant hypertension, among others [3, 4]. However, the prevalence of secondary causes of hypertension is often underestimated and many cases remain unrecognized [6,7,8]. Early detection of endocrine forms of hypertension is crucial to control high blood pressure and prevent hypertension-mediated organ damage and related cardiovascular complications, as well as achieve effective long-term treatment [9]. Moreover, hormone excess may increase individual risk of other consequences beyond hypertension, particularly in the case of CS or PA [10,11,12].

Currently, the diagnosis of endocrine hypertension (EHT) relies on hormonal evaluations, with specific diagnostic algorithms for each type of endocrine disorder. Hormone assays present some limitations, including (i) common borderline values, especially in mild forms of over-secretion, where the diagnosis may be missed or wrongly called; (ii) lack of direct estimation of the individual risk to hormone exposure, despite the important inter-individual variability of hormone excess consequences; (iii) recommendation for multiple hormonal tests, with screening and confirmatory strategies still debated, and varying from one center to the other [13]. New biomarkers could potentially allow to directly measure tissue exposure to hormone excess, with the potential improvement of diagnostic accuracy and prediction of individual susceptibility to different consequences.

Circulating biomarkers in the blood have the advantage of being non-invasive when utilized. Leukocyte DNA methylation is particularly convenient as a biomarker, being a chemically stable yet dynamic biological hallmark, with a key role in epigenetic regulation in both health and disease [14]. Recent epigenome-wide association studies exploring leukocyte DNA methylation in hypertensive patients versus normotensive controls identified different loci associated with blood pressure regulation [15,16,17]. However, it is not yet established whether blood methylation profiles could also discriminate patients with EHT from those with primary hypertension (PHT). Indeed, hormone excess may impact peripheral tissues at the epigenetic level, measurable as hormone-specific methylome signatures. This hypothesis is supported by a recent work, in which we were able to identify a blood methylome signature of CS [18].

In the present study, we explored leukocyte methylation profiles in hypertensive patients. Specifically, we analyzed whole blood methylome signatures in patients with PHT or EHT related to CS, PA, or PPGL.



Blood samples were collected from 255 patients. Patients had been diagnosed either with primary hypertension (PHT, n = 42), or endocrine forms of hypertension (EHT, n = 213), either related to the presence of CS (n = 57), PA (n = 101), or PPGL (n = 55). Each group, except PA, showed a predominance of female patients, and the mean age was lower in patients with EHT than in patients with PHT (p < 0.01; Table 1).

Table 1 Sex and age distribution in the different hypertension types

Whole blood DNA methylome profile in hypertension

Whole-genome blood DNA methylome was determined for the 255 samples, with 731,635 informative CpG sites in all samples. Unsupervised principal component analysis showed a distribution of samples with separation depending on the type of hypertension (Fig. 1a). The main components of variability were associated with the type of hypertension, but also white blood cell count variation and, to a lower extent, with age and sex (Additional file 1: Fig. S1).

Fig. 1
figure 1

Global structure of blood DNA methylation profile in different hypertension types. a Sample projections based on the first two principle components (PC1, PC2) of unsupervised PCA performed on the whole dataset (n = 731,635 CpG sites, n = 255 samples). The center of each group is indicated by the larger white circles. b Boxplot of the most variable CpG sites (n = 48,452 with M-Value SD > 0.4) in the four hypertension types. #CS, PA and PPGL were aggregated into a single EHT type, with random sampling of 55 samples from each group. *Student’s t-test p value < 0.05

The mean methylation level of the most variable CpG sites among all samples (n = 48,452) was not significantly different between the different types of hypertension, except in the group of CS patients, where methylation was decreased (Fig. 1b).

The methylome signature of endocrine hypertension

Comparing the 731,635 informative CpG sites of each type of EHT to PHT, differentially methylated CpG sites were identified in EHT (n = 123,288), CS (n = 197,897), PA (n = 114,988), and PPGL (n = 66,976) versus PHT, respectively (Additional file 2: Table S1). Differentially methylated CpG sites were distributed all over the genome (Fig. 2c–f), and the proportion of hypo- and hyper-methylated CpG sites was not strongly related to CpG genomic density (Fig. 2a), nor to CpG proximity to genes (Fig. 2b). Gene set enrichment analysis of genes associated with the differentially methylated CpG sites in the four comparisons revealed enrichment in pathways related to blood pressure regulation mechanisms and hypertension, including MAPK, Rap1, phospholipase D and calcium signaling pathways [19,20,21,22,23] (FDR < 0.001 in all comparisons; Additional file 3: Table S2).

Fig. 2
figure 2

Characteristics of differentially methylated CpG sites in endocrine hypertension. a Distribution of differentially methylated CpG sites related to CpG genomic density. b Distribution of differentially methylated CpG sites related to genes structure. cf Manhattan plots representing the genomic distribution of differentially methylated CpG sites for each type of endocrine hypertension, in the comparisons CS vs. PHT (c), PA vs. PHT (d), PPGL vs. PHT (e), EHT vs. PHT (f). The most differentially methylated genes (mean –log10(p value) > 11, Limma) are highlighted in black

Beyond the analysis of individual CpG sites, a specific analysis of differentially methylated genes was performed (Table 2 and Additional file 4: Table S3), highlighting the implication of specific genes for each type of EHT (Fig. 2c–f). The most significantly differentially methylated genes were FKBP5 in CS and CD300LG in PA.

Table 2 Top 15 genes associated with the most differentially methylated CpG sites for each endocrine hypertension type

Prediction of endocrine hypertension

Four different machine learning methods—Lasso, Logistic Regression, Random Forest, Support Vector Machine—were used to build a prediction model for each type of endocrine hypertension on subsets of samples (training cohorts), and subsequently tested on remaining samples (validation cohorts). The prediction performance was better for individual types of EHT (CS, PA and PPGL) than for predicting EHT as a whole. Indeed, balanced accuracy against PHT reached 0.95 for CS using SVM, 0.88 for PA using Lasso, 0.83 for PPGL using RF, and 0.74 for EHT using LR (Fig. 3a; selected CpG sites in Additional file 1: Table S4). Misclassified samples (false negatives and false positives) mostly showed a positioning on the global methylome space in regions of overlapping between different hypertension types (Additional file 1: Fig. S2). In the two comparisons with sample numbers imbalance (PA vs. PTH and EHT vs. PHT), down- and up-sampling were used to further explore prediction score, with no major impact on performance (Additional file 1: Fig. S3). The selected sets of CpG sites were specific for each endocrine hypertension type (Fig. 3b). These specific CpG sites, except for just one CpG site in CS, were all present among the differentially methylated CpG sites selected for each type of endocrine hypertension (Additional file 2: Table S1).

Fig. 3
figure 3

Prediction of endocrine hypertension by blood methylome marks. a Heatmap of prediction scores in validation cohorts for each type of endocrine hypertension (CS, PA, PPGL, EHT) versus primary hypertension (PHT), obtained with different machine learning approaches. RF random forests, SVM support vector machine, LR logistic regression, Lasso Least Absolute Shrinkage and Selection Operator. b Venn diagram of top performing features selected for each comparison (prediction score of reference: balanced accuracy)


In this work, we were able to identify a blood methylome signature discriminating different types of EHT -including CS, PA and PPGL- from PHT. These different conditions are currently diagnosed with specific hormone measurements, representing the gold standard [9]. However, hormone assays are not always conclusive, with accuracies far from perfect [24,25,26], and quite expensive. This leads to screening strategies that target hypertensive patients with the highest risk of endocrine etiology [3, 4]. These probabilistic approaches could miss cases compared to a full screening of the entire hypertensive population [27]. Specific non-hormonal biomarkers could potentially help in situations where (i) hormone levels are not applicable—for instance, exogenous glucocorticoid administration; (ii) hormone dosages are not informative, due to situations interfering with hormone assays [28]; (iii) hormone dosages are borderline, not providing a proper classification of patients. Some non-hormonal markers of endocrine hypertension have been proposed, including specifically targeted metabolomics [29], and miRNAs [30, 31]. Blood methylation marks may also be suitable for direct measurement of hormone excess on peripheral tissues. Indeed, DNA methylation is dynamic but stably preserved [14] and could potentially reflect the individual impact of hormone excess over time better than hormonal evaluation, that need to be standardized for time, feeding-status, salt and fluid intake, and co-medication, among others.

Beyond the specific benefits for diagnosis, massive molecular screenings in common diseases like hypertension are expected to improve personalized medicine, through the identification of features predicting response to treatment, or specific complications. Addressing these aims will require specific study designs, focused on follow-up. For these up-coming studies, robustness and statistical power will improve by restricting the number of candidate features to be tested. This prior feature selection will be helped by this work and others, providing lists of CpG sites associated with each type of hypertension or with each type of hypertensive treatment [32].

Of note, the methylome signatures identified here are impacted by white blood cell count. Indeed, neutrophils increase and lymphocytes decrease in case of CS and, to a lesser extent, of PA [33,34,35,36]. However, we showed that white blood cell count was contributing to a limited part of methylome variability. To consider this effect, the methylome signatures provided here were adjusted for white blood cell count. The proportion of neutrophils was used as a unique proxy for this adjustment, given the high correlation with proportions of lymphocytes and lymphocyte subtypes.

Here, we identified blood DNA methylation loci associated with EHT. The blood DNA methylome profile associated with the different types of hormonal excess was also explored, with the identification of genes differentially methylated, such as FKBP5 in CS and CD300LG in PA. FKBP5 encodes a co-chaperone regulating the glucocorticoid receptor activity [37], and its increased expression and decreased methylation in blood have already been described in CS [18, 38]. Moreover, stress-induced epigenetic up-regulation of FKBP5 has been associated with an increased risk of cardiovascular events [39]. CD300LG encodes a vascular endothelial cell adhesion molecule, which is implicated in leukocyte binding and transmigration, and whose expression seems to be differently influenced by the local environment in different tissues [40]. Aldosterone is known to induce endothelial dysfunction [41] and to promote a pro-inflammatory phenotype, with vascular infiltration of monocytes, macrophages, and lymphocytes, as demonstrated in both animal models of hypertension [42, 43] and human vascular endothelial cells [44]. Whether CD300LG possibly contributes to the effect of aldosterone on vascular endothelium and lymphocyte migration, and whether this association represents a specific trait of aldosterone-induced hypertension, remains to be established. One study explored CD300LG variants in a mouse model of hypertension, but with uncertain conclusions [45]. For a majority of other genes, the biological relevance in endocrine hypertension remains to be explored, including the impact of DNA methylation of gene expression.

The hypothesis of a common blood methylome signature of EHT was also addressed. When we compared CS-, PA- and PPGL-related signatures, a limited number of methylation marks were shared, suggesting the prominence of individual endocrine hypertension signatures over a common signature of EHT. However, sample projections based on the variability of methylation levels showed overlap between PA and PPGL, and a tendency of CS to cluster separately. Indeed, secondary aldosteronism is common in PPGL [46] and could in part explain this similarity with PA. To which extent an endocrine-related molecular signature of hypertension shared between CS, PA and PPGL, but distinct from PHT, is induced remains to be determined. By directly comparing EHT as a whole to PHT, we provide here a set of methylome marks, which may facilitate the identification of such molecular signatures.

Whether these methylation marks better reflect the risk of complications related to hormone excess also remains to be elucidated. Of note, the methylation signatures identified here were based on whole methylome analysis. Routine use of such a tool would be limited by cost and complexity. Reliable biomarkers should allow the prediction of a clinical endpoint of interest more easily, at a lower cost and over a shorter time span than the direct measurement of the clinical end point [47]. For methylation marks, this would imply a technology transfer to targeted measurements, such as pyrosequencing, methylation specific-MLPA, methylation specific-high resolution melting analysis, or nanopore sequencing [48, 49]. Setup and validation of suitable targeted assays have to be performed for the identified methylome marks of EHT.

A potential interest of selected methylome marks could be the integration in a targeted multi-omic chip for the diagnosis of endocrine hypertension (ENSAT-HT project; Whether these methylation marks could positively impact the prediction of a multi-modal classifier remains to be determined, especially for sensitive detection of the rarest forms of endocrine hypertension, namely CS and PPGL.

This work aimed to develop machine learning models with methylome data to predict endocrine hypertension. This approach might be seen as too simple regarding the pathophysiological complexity of these conditions, and the extensive and complex diagnosis workup currently needed to call these conditions. However, machine learning could help in leading a step forward toward a robust and accurate screening tool and could provide a unique opportunity to understand the complex relationships within an omic dataset and a complex phenotype. On the other side, this data-driven strategy may be impacted by biases, including recruitment biases—in expert centers and optimal conditions-, and representation biases—considering the rare prevalence of CS or PPGL compared to the other types of hypertension. This latter risk is mitigated here by the study design splitting samples into independent training and validation cohorts. Finally, the heterogeneity of sample sizes for the different types of endocrine hypertension may impact the machine learning process. However, the group of EHT was balanced in terms of endocrine hypertension types when analyzed as a whole, and optimized strategies were used for unbalanced comparisons, including balanced accuracy measurement, or down- and up-sampling of the training cohorts.


Endocrine hypertension can be identified by blood DNA methylation markers, with specific methylation signatures for each type of endocrine disorder.


Patients and samples

A total of 255 blood samples were collected in eight different centers, as part of the European ENSAT-HT study ( The cohort included patients with PHT (n = 42) and patients with different types of EHT, namely CS (n = 57), PA (n = 101), or PPGL (n = 55) (Additional file 5: Table S5). The diagnosis was made according to the current guidelines for screening and management of each specific disease [4, 50,51,52]. Diagnosis of PHT also required the exclusion of endocrine hypertension and other secondary causes (renal disease, pharmacological cause and obstructive sleep apnea syndrome) as well as the exclusion of patients with low-renin hypertension. Patients with uncertain diagnosis, those with pregnancy, severe comorbidities (including heart failure, chronic kidney disease, active malignancy) were also excluded.

Whole-genome DNA methylation measurement

Leukocyte DNA was extracted from EDTA blood samples, using the DNA Isolation kit for Mammalian Blood (Roche, Basil, Switzerland). DNA quality was assessed on a Genomic DNA ScreenTape system (Agilent, Santa Clara, CA, US) and quantified using a Qubit 3.0 Fluorometer (Thermofisher, Waltham, MA, US). DNA was treated by bisulfite and then hybridized to the Infinium MethylationEPIC BeadChip (Illumina, San Diego, CA, US; ~ 865,000 sites), starting from 500 ng of DNA. All experiments were performed following the manufacturer’s instructions at the P3S Post-Genomic Platform of Sorbonne University (Paris, France).

Bioinformatics and statistics

All samples passed the quality controls provided by the Genome Studio software (v. 2011.1; Illumina). Data were exported in Intensity Data (IDAT) format and then processed using the minfi package (v. 1.32.0) [53] in the R software environment (v. 3.6.3) ( An additional quality control of samples was performed based on probe detection quality, confirming the good quality of each sample (mean detection p value < 0.05).

Data were normalized using the stratified quantile normalization procedure implemented in the “preprocessQuantile” minfi function [54] and the methylation score for each CpG probe was extracted as a β-value. Then, the ChAMP package (v. 2.16.1) was used to filter the probes [55]. A total of 731,635 CpG sites passed the following criteria: detection p value < 0.01, presence of the targeted CpG site, absence of frequent SNPs in the probe, single hybridization hit, autosomal target. Of note, filtering the probes before normalization versus after normalization did not significantly impact the values (data not shown).

The significant components of variation in the dataset were assessed using the singular value decomposition method (SVD) for methylation data [56], and a detected batch effect (Slide) was corrected using the “ComBat” method [57], as implemented in the ChAMP package.

White blood cell count of subpopulations (neutrophils, lymphocytes B, lymphocytes T4, lymphocytes T8, lymphocytes NK, monocytes) were estimated by the reference-based “RefbaseEWAS” method [58] implemented in the ChAMP package (Additional file 6: Table S6). Since neutrophils were the most represented cell type in all samples, and since the proportions of neutrophils and lymphocytes were anti-correlated (Pearson’s r = − 0.96; Additional file 1: Fig. S1c, d), the estimated proportion of neutrophils was chosen as the unique proxy reflecting variations in white blood cell count.

M-values, used for statistical analysis, were calculated from β-values (log2 ratio of the intensities of methylated vs. unmethylated probes) using the lumi package (v. 2.36.0) [59].

The global data structure was assessed on β-values by principal component analysis (PCA), using all the CpG probes. The most variable CpG probes were selected based on their Standard Deviation (SD cutoff: 0.4, n = 48,452 probes) among all samples. Differentially methylated CpG sites were identified using the Limma package (v. 3.40.6) [60] for each of the following comparison: EHT vs. PHT, CS vs. PHT, PA vs. PHT, PPGL vs. PHT. For the EHT vs. PHT comparison, the EHT group was obtained by a random sampling of CS, PA and PPGL (n = 55 samples for each type of endocrine hypertension), to avoid imbalance between the endocrine hypertension types. The estimated neutrophil count and age were included as covariates and CpG sites were selected based on an adjusted p value < 0.05 (Benjamin-Hochberg correction procedure; Additional file 2: Table S1). Gene set enrichment analysis of genes associated with differentially methylated CpG sites was performed using the “gometh” method (KEGG chosen as collection of pathways to test) implemented in the missmethyl package (v. 1.18.0) [61], adjusting for the number of CpG sites associated to each gene [62] (Additional file 3: Table S2). For each endocrine hypertension type, differentially methylated genes were selected when they presented at least three differentially methylated CpG sites in each comparison and a − log10(p value) > 5 (Additional file 4: Table S3).

For predicting endocrine hypertension, the same four comparisons described above were considered for supervised learning (CS vs. PHT, PA vs. PHT, PPGL vs. PHT, and EHT vs. PHT). For each comparison, training (80%) and validation (20%) samples were randomly selected by maintaining the initial proportion of samples in each hypertension type in each comparison (Additional file 1: Table S7). Starting from the most variable CpG sites previously selected (n = 48,452 probes with M-value SD > 0.4), the most discriminating CpG sites for each comparison (n = 200 for CS vs. PHT, n = 141 for PA vs. PHT, n = 70 for PPGL vs. PHT, and n = 135 for EHT vs. PHT) were pre-selected from the training sets using the boruta package (v. 0.3) [63] (Additional file 1: Table S8). For each comparison, four different classifier methods were then used: (i) penalized Lasso (Least Absolute Shrinkage and Selection Operator) regression with tenfold cross-validation, using the glmnet package (v. 4.0–2) [64]; (ii) Support Vector Machine (SVM) [65]; iii) Random Forest (RF) [66]; iv) Logistic Regression (LR) [67] with Fast Correlation-Based Filter (FCBF) [68] using the orange toolbox (v. 3.30.1) [69] (Additional file 1: Table S9). The trained models were tested on the validation sets for each comparison, and the prediction performance was evaluated using balanced accuracy, sensitivity, specificity, F1, Kappa, and AUC scores. The model training was further evaluated using up-sampling and down-sampling (70).

Group comparisons were performed using Student’s t test for variables normally distributed, or chi-square test for binary categorical variables. All tests were computed in the R software environment.

Availability of data and materials

The methylome dataset generated and analyzed during the current study is available in the EMBL-EBI BioStudies repository (Reference Number: S-BSST831).


  1. Chow CK, Teo KK, Rangarajan S, Islam S, Gupta R, Avezum A, et al. Prevalence, awareness, treatment, and control of hypertension in rural and urban communities in high-, middle-, and low-income countries. JAMA. 2013;310(9):959–68.

    Article  CAS  PubMed  Google Scholar 

  2. NCD Risk Factor Collaboration (NCD-RisC). Worldwide trends in blood pressure from 1975 to 2015: a pooled analysis of 1479 population-based measurement studies with 19·1 million participants. Lancet. 2017;389(10064):37–55.

  3. Whelton PK, Carey RM, Aronow WS, Casey DE, Collins KJ, Dennison Himmelfarb C, et al. 2017 ACC/AHA/AAPA/ABC/ACPM/AGS/APhA/ASH/ASPC/NMA/PCNA guideline for the prevention, detection, evaluation, and management of high blood pressure in adults: a report of the American college of cardiology/American heart association task force on clinical practice guidelines. Hypertension. 2018;71(6):e13-115.

    CAS  PubMed  Google Scholar 

  4. Williams B, Mancia G, Spiering W, Agabiti Rosei E, Azizi M, Burnier M, et al. 2018 ESC/ESH Guidelines for the management of arterial hypertension: The Task Force for the management of arterial hypertension of the European Society of Cardiology (ESC) and the European Society of Hypertension (ESH). Eur Heart J. 2018;39(33):3021–104.

    Article  PubMed  Google Scholar 

  5. Hirsch JS, Hong S. The demystification of secondary hypertension: diagnostic strategies and treatment algorithms. Curr Treat Options Cardiovasc Med. 2019;21(12):90.

    Article  PubMed  Google Scholar 

  6. Börgel J, Springer S, Ghafoor J, Arndt D, Duchna HW, Barthel A, et al. Unrecognized secondary causes of hypertension in patients with hypertensive urgency/emergency: prevalence and co-prevalence. Clin Res Cardiol. 2010;99(8):499–506.

    Article  PubMed  Google Scholar 

  7. Brown JM, Siddiqui M, Calhoun DA, Carey RM, Hopkins PN, Williams GH, et al. The unrecognized prevalence of primary aldosteronism: a cross-sectional study. Ann Intern Med. 2020;173(1):10–20.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Kotliar C, Obregón S, Koretzky M, Botto F, Di Leva A, Boscaro M, et al. Improved identification of secondary hypertension: use of a systematic protocol. Ann Transl Med. 2018;6(15):293.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Rossi GP, Bisogni V, Rossitto G, Maiolino G, Cesari M, Zhu R, et al. Practice recommendations for diagnosis and treatment of the most common forms of secondary hypertension. High Blood Press Cardiovasc Prev. 2020;27(6):547–60.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Redheuil A, Blanchard A, Pereira H, Raissouni Z, Lorthioir A, Soulat G, et al. Aldosterone-related myocardial extracellular matrix expansion in hypertension in humans: a proof-of-concept study by cardiac magnetic resonance. JACC Cardiovasc Imaging. 2020;13(10):2149–59.

    Article  PubMed  Google Scholar 

  11. Stowasser M, Gordon RD. Primary aldosteronism: changing definitions and new concepts of physiology and pathophysiology both inside and outside the kidney. Physiol Rev. 2016;96(4):1327–84.

    Article  CAS  PubMed  Google Scholar 

  12. Pivonello R, Isidori AM, De Martino MC, Newell-Price J, Biller BMK, Colao A. Complications of Cushing’s syndrome: state of the art. Lancet Diabetes Endocrinol. 2016;4(7):611–29.

    Article  CAS  PubMed  Google Scholar 

  13. de Freminville JB, Amar L. How to explore an endocrine cause of hypertension. J Clin Med. 2022;11(2):420.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Schübeler D. Function and information content of DNA methylation. Nature. 2015;517(7534):321–6.

    Article  PubMed  Google Scholar 

  15. Huang Y, Ollikainen M, Muniandy M, Zhang T, van Dongen J, Hao G, et al. Identification, heritability, and relation with gene expression of novel DNA methylation loci for blood pressure. Hypertension. 2020;76(1):195–205.

    Article  CAS  PubMed  Google Scholar 

  16. Kazmi N, Elliott HR, Burrows K, Tillin T, Hughes AD, Chaturvedi N, et al. Associations between high blood pressure and DNA methylation. PLoS ONE. 2020;15(1): e0227728.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Richard MA, Huan T, Ligthart S, Gondalia R, Jhun MA, Brody JA, et al. DNA methylation analysis identifies loci for blood pressure regulation. Am J Hum Genet. 2017;101(6):888–902.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Armignacco R, Jouinot A, Bouys L, Septier A, Lartigue T, Neou M, et al. Identification of glucocorticoid-related molecular signature by whole blood methylome analysis. Eur J Endocrinol. 2022;186(2):297–308.

    Article  CAS  PubMed  Google Scholar 

  19. Nelson RK, Ya-Ping J, Gadbery J, Abedeen D, Sampson N, Lin RZ, et al. Phospholipase D2 loss results in increased blood pressure via inhibition of the endothelial nitric oxide synthase pathway. Sci Rep. 2017;7(1):9112.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Lakshmikanthan S, Zieba BJ, Ge ZD, Momotani K, Zheng X, Lund H, et al. Rap1b in smooth muscle and endothelium is required for maintenance of vascular tone and normal blood pressure. Arterioscler Thromb Vasc Biol. 2014;34(7):1486–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Wilson C, Zhang X, Buckley C, Heathcote HR, Lee MD, McCarron JG. Increased vascular contractility in hypertension results from impaired endothelial calcium signaling. Hypertension. 2019;74(5):1200–14.

    Article  CAS  PubMed  Google Scholar 

  22. Touyz RM, Alves-Lopes R, Rios FJ, Camargo LL, Anagnostopoulou A, Arner A, et al. Vascular smooth muscle contraction in hypertension. Cardiovasc Res. 2018;114(4):529–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Lu Y, Sun X, Peng L, Jiang W, Li W, Yuan H, et al. Angiotensin II-Induced vascular remodeling and hypertension involves cathepsin L/V- MEK/ERK mediated mechanism. Int J Cardiol. 2020;1(298):98–106.

    Article  Google Scholar 

  24. Bansal V, Asmar NE, Selman WR, Arafah BM. Pitfalls in the diagnosis and management of Cushing’s syndrome. Neurosurg Focus. 2015;38(2):E4.

    Article  PubMed  Google Scholar 

  25. Rossi GP. Primary aldosteronism: JACC state-of-the-art review. J Am Coll Cardiol. 2019;74(22):2799–811.

    Article  CAS  PubMed  Google Scholar 

  26. Ilias I, Thomopoulos C. Addressing delays in the diagnosis of pheochromocytoma/paraganglioma. Expert Rev Endocrinol Metab. 2019;14(5):359–63.

    Article  CAS  PubMed  Google Scholar 

  27. Ruhle BC, White MG, Alsafran S, Kaplan EL, Angelos P, Grogan RH. Keeping primary aldosteronism in mind: deficiencies in screening at-risk hypertensives. Surgery. 2019;165(1):221–7.

    Article  PubMed  Google Scholar 

  28. Lopez AG, Fraissinet F, Lefebvre H, Brunel V, Ziegler F. Pharmacological and analytical interference in hormone assays for diagnosis of adrenal incidentaloma. Ann Endocrinol (Paris). 2019;80(4):250–8.

    Article  Google Scholar 

  29. Erlic Z, Reel P, Reel S, Amar L, Pecori A, Larsen CK, et al. Targeted metabolomics as a tool in discriminating endocrine from primary hypertension. J Clin Endocrinol Metab. 2021;106(4):1111–28.

    Article  PubMed  Google Scholar 

  30. MacKenzie SM, Saunders H, van Kralingen JC, Robertson S, Riddell A, Zennaro MC, et al. Circulating microRNAs as diagnostic markers in primary aldosteronism. Cancers (Basel). 2021;13(21):5312.

    Article  CAS  Google Scholar 

  31. Vetrivel S, Zhang R, Engel M, Altieri B, Braun L, Osswald A, et al. Circulating microRNA expression in cushing’s syndrome. Front Endocrinol (Lausanne). 2021;12: 620012.

    Article  Google Scholar 

  32. Nuotio ML, Sánez Tähtisalo H, Lahtinen A, Donner K, Fyhrquist F, Perola M, et al. Pharmacoepigenetics of hypertension: genome-wide methylation analysis of responsiveness to four classes of antihypertensive drugs using a double-blind crossover study design. Epigenetics. 2022;25:1–14.

    Google Scholar 

  33. Masri-Iraqi H, Robenshtok E, Tzvetov G, Manistersky Y, Shimon I. Elevated white blood cell counts in Cushing’s disease: association with hypercortisolism. Pituitary. 2014;17(5):436–40.

    Article  PubMed  Google Scholar 

  34. Nakagawa M, Terashima T, D’yachkova Y, Bondy GP, Hogg JC, van Eeden SF. Glucocorticoid-induced granulocytosis: contribution of marrow release and demargination of intravascular granulocytes. Circulation. 1998;98(21):2307–13.

    Article  CAS  PubMed  Google Scholar 

  35. Libianto R, Hu J, Chee MR, Hoo J, Lim YY, Shen J, et al. A Multicenter study of neutrophil-to-lymphocyte ratio in primary aldosteronism. J Endocr Soc. 2020;4(12):bvaa153.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Gilet A, Zou F, Boumenir M, Frippiat JP, Thornton SN, Lacolley P, et al. Aldosterone up-regulates MMP-9 and MMP-9/NGAL expression in human neutrophils through p38, ERK1/2 and PI3K pathways. Exp Cell Res. 2015;331(1):152–63.

    Article  CAS  PubMed  Google Scholar 

  37. Fries GR, Gassen NC, Rein T. The FKBP51 glucocorticoid receptor co-chaperone: regulation, function, and implications in health and disease. Int J Mol Sci. 2017;18(12):E2614.

    Article  PubMed  Google Scholar 

  38. Bancos I, Hatipoglu BA, Yuen KCJ, Chandramohan L, Chaudhari S, Moraitis AG. Evaluation of FKBP5 as a cortisol activity biomarker in patients with ACTH-dependent Cushing syndrome. J Clin Transl Endocrinol. 2021;1(24): 100256.

    Google Scholar 

  39. Zannas AS, Jia M, Hafner K, Baumert J, Wiechmann T, Pape JC, et al. Epigenetic upregulation of FKBP5 by aging and stress contributes to NF-κB–driven inflammation and cardiovascular risk. PNAS. 2019;116(23):11370–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Umemoto E, Takeda A, Jin S, Luo Z, Nakahogi N, Hayasaka H, et al. Dynamic changes in endothelial cell adhesion molecule nepmucin/CD300LG expression under physiological and pathological conditions. PLoS ONE. 2013;8(12): e83681.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Liu G, Yin GS, Tang JY, Ma DJ, Ru J, Huang XH. Endothelial dysfunction in patients with primary aldosteronism: a biomarker of target organ damage. J Hum Hypertens. 2014;28(12):711–5.

    Article  CAS  PubMed  Google Scholar 

  42. Rocha R, Rudolph AE, Frierdich GE, Nachowiak DA, Kekec BK, Blomme EAG, et al. Aldosterone induces a vascular inflammatory phenotype in the rat heart. Am J Physiol Heart Circ Physiol. 2002;283(5):H1802-1810.

    Article  CAS  PubMed  Google Scholar 

  43. Sun Y, Zhang J, Lu L, Chen SS, Quinn MT, Weber KT. Aldosterone-induced inflammation in the rat heart : role of oxidative stress. Am J Pathol. 2002;161(5):1773–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Caprio M, Newfell BG, la Sala A, Baur W, Fabbri A, Rosano G, et al. Functional mineralocorticoid receptors in human vascular endothelial cells regulate intercellular adhesion molecule-1 expression and promote leukocyte adhesion. Circ Res. 2008;102(11):1359–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Støy J, Grarup N, Hørlyck A, Ibsen L, Rungby J, Poulsen PL, et al. Blood pressure levels in male carriers of Arg82Cys in CD300LG. PLoS ONE. 2014;9(10): e109646.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Mao JJ, Baker JE, Rainey WE, Young WF Jr, Bancos I. Concomitant pheochromocytoma and primary aldosteronism: a case series and literature review. J Endoc Soc. 2021;5(8):e b107.

  47. Aronson JK, Ferner RE. Biomarkers-A General Review. Curr Protoc Pharmacol. 2017;76:9231–92317.

    Article  Google Scholar 

  48. Simpson JT, Workman RE, Zuzarte PC, David M, Dursi LJ, Timp W. Detecting DNA cytosine methylation using nanopore sequencing. Nat Methods. 2017;14(4):407–10.

    Article  CAS  PubMed  Google Scholar 

  49. García-Giménez JL, Seco-Cervera M, Tollefsbol TO, Romá-Mateo C, Peiró-Chova L, Lapunzina P, et al. Epigenetic biomarkers: Current strategies and future challenges for their use in the clinical laboratory. Crit Rev Clin Lab Sci. 2017;54(7–8):529–50.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Funder JW, Carey RM, Mantero F, Murad MH, Reincke M, Shibata H, et al. The management of primary aldosteronism: case detection, diagnosis, and treatment: an endocrine society clinical practice guideline. J Clin Endocrinol Metab. 2016;101(5):1889–916.

    Article  CAS  PubMed  Google Scholar 

  51. Lenders JWM, Duh QY, Eisenhofer G, Gimenez-Roqueplo AP, Grebe SKG, Murad MH, et al. Pheochromocytoma and paraganglioma: an endocrine society clinical practice guideline. J Clin Endocrinol Metab. 2014;99(6):1915–42.

    Article  CAS  PubMed  Google Scholar 

  52. Nieman LK, Biller BMK, Findling JW, Newell-Price J, Savage MO, Stewart PM, et al. The diagnosis of Cushing’s syndrome: an Endocrine Society Clinical Practice Guideline. J Clin Endocrinol Metab. 2008;93(5):1526–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. 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(10):1363–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Touleimat N, Tost J. Complete pipeline for Infinium(®) Human Methylation 450K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics. 2012;4(3):325–41.

    Article  CAS  PubMed  Google Scholar 

  55. Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, et al. ChAMP: 450k chip analysis methylation pipeline. Bioinformatics. 2014;30(3):428–30.

    Article  CAS  PubMed  Google Scholar 

  56. Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Gayther SA, Apostolidou S, et al. An epigenetic signature in peripheral blood predicts active ovarian cancer. PLoS ONE. 2009;4(12): e8274.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.

    Article  PubMed  Google Scholar 

  58. 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 Bioinform. 2012;13(1):86.

    Article  Google Scholar 

  59. Du P, Kibbe WA, Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008;24(13):1547–8.

    Article  CAS  PubMed  Google Scholar 

  60. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7): e47.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Phipson B, Maksimovic J, Oshlack A. missMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics. 2016;32(2):286–8.

    Article  CAS  PubMed  Google Scholar 

  62. Geeleher P, Hartnett L, Egan LJ, Golden A, Raja Ali RA, Seoighe C. Gene-set analysis is severely biased when applied to genome-wide methylation data. Bioinformatics. 2013;29(15):1851–7.

    Article  CAS  PubMed  Google Scholar 

  63. Kursa MB, Rudnicki WR. Feature selection with the boruta package. J Stat Softw. 2010;16(36):1–13.

    Google Scholar 

  64. Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33(1):1–22.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Chang CC, Lin CJ. LIBSVM: a library for support vector machines. ACM Trans Intell Syst Technol. 2011;2(3):271–2727.

    Article  Google Scholar 

  66. Breiman L. Random forests. Mach Learn. 2001;45(1):5–32.

    Article  Google Scholar 

  67. McCullagh P, Nelder JA. Generalized linear models. 2nd ed. Boca Raton: Routledge; 2019. p. 532.

    Book  Google Scholar 

  68. Yu L, Liu H. Feature Selection for High-Dimensional Data: Proceedings, Twentieth International Conference on Machine Learning. Fawcett T, Mishra N, editors. Proceedings, Twentieth International Conference on Machine Learning. 2003 Dec 1;856–63.

  69. Demšar J, Curk T, Erjavec A, Gorup Č, Hočevar T, Milutinovič M, et al. Orange: data mining toolbox in python. J Mach Learn Res. 2013;14(1):2349–53.

    Google Scholar 

  70. Kuhn M, Johnson K. Applied Predictive Modeling. ISBN: 978-1-4614-6849-3. Springer Nature, 2013.

Download references


We thank Gisèle Bonne, Badreddine Mohand Oumoussa and Abiba Doukani for their technical supply at the P3S Post-Genomic Platform of Sorbonne University, as well as Gabriel Valeix and Meriama Saidi for their clinical research assistance at Cochin hospital (APHP, Paris, France). We thank Kerstin Schaefer and Stephanie Zopp for their research assistance at LMU hospital (Ludwig-Maximilians-University, Munich, Germany), and Thomas Baumgartner at the Universitätsspital Zürich (Zurich, Switzerland).


This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 633983, the Programme Hospitalier de Recherche Clinique “CompliCushing” (PHRC AOM 12-002-0064), the Agence Nationale pour la Recherche (ANR-18-CE14-0008-01), the Else Kröner-Fresenius Stiftung (2012_A103 and 2015_A228 to MR) and the Deutsche Forschungsgemeinschaft (DFG) within the CRC/Transregio 205/1 “The Adrenal: Central Relay in Health and Disease” (to MR, FB, AR). The work was further supported by the Clinical Research Priority Program of the University of Zurich for the CRPP HYRENE (to FB).

Author information

Authors and Affiliations



RA, PSR, SR, GA contributed to conceptualization; RA, PSR, SR, GA contributed to methodology; RA, PSR, SR contributed to Software; CKL, LB, LB, AR, MK, FBS, LA, AB, APGR, AP, AJ, PD, ED, SMM, GPR, LL, FC, CS, PM, TAW, AP, SM, FB, MR, MCZ, JB, GA contributed to clinical data and sample collection; RA, PSR, SR, GA contributed to visualization; LA, AB, APGR, AP, ED, GPR, LL, CS, PM, FB, MR, MCZ, JB, EJ, GA contributed to project administration and ethical aspects management; RA, CG, KP contributed to samples handling and genomic data generation; MCZ, EJ, GA contributed to supervision; RA, PSR, SR, AJ, AS, EJ, GA contributed to bioinformatics and statistical analyses; RA, GA contributed to original manuscript preparation. All authors reviewed and approved the final manuscript.

Corresponding authors

Correspondence to Roberta Armignacco or Guillaume Assié.

Ethics declarations

Ethics approval and consent to participate

Signed informed consent for molecular analysis of blood samples and for access to clinical data was obtained from all patients, and the study was approved by the Ethic Committee of each participating center.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Components of variation in the whole methylome dataset. a) Scree plot representing the percentage of explained variability by the first five principal components of PCA performed on the whole dataset (n=731,635 CpG sites, n=255 samples). b) Singular value decomposition (SVD) plot assessing the correlation between the first five significant components of variation in the dataset and other biological factors of interest (Hypertension type –PHT, CS, PA, PPGL-, age, sex, cell composition–neutrophils used as the unique proxy). c) Correlation between the proportion of neutrophils and of lymphocytes. d) Correlation plot between the proportion of neutrophils and of the other different cell subtypes. Fig. S2. Misclassified sample positioning in the global structure of blood DNA methylation. Samples with discrepant methylome prediction and hormonal status are indicated by the white squares. Fig. S3. Heatmap of prediction scores in validation cohorts for primary aldosteronism (PA) and endocrine hypertension (EHT) versus primary hypertension (PHT), obtained with different machine learning approaches, after up-sampling and down-sampling the training cohort. RF random forests, SVM support vector machine; LR logistic regression; Lasso Least Absolute Shrinkage and Selection Operator. Table S4. Top performing methylome features. List of CpG sites selected for predicting each type of endocrine hypertension, based on the best performing method for each comparison. Table S7. Training/Validation split of samples. Table S8. Parameters chosen for Boruta package. For details, see Table S9. Parameters chosen for different machine learning models for classification task.

Additional file 2: Table S1.

Differentially methylated CpG sites in endocrine hypertension.

Additional file 3: Table S2.

Enriched signaling pathways in endocrine hypertension.

Additional file 4: Table S3.

Significant differentially methylated genes in endocrine hypertension.

Additional file 5: Table S5.

Sample characteristics.

Additional file 6: Table S6.

Methylome-based estimation of blood cell composition.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Armignacco, R., Reel, P.S., Reel, S. et al. Whole blood methylome-derived features to discriminate endocrine hypertension. Clin Epigenet 14, 142 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Endocrine hypertension
  • Whole blood methylome
  • Circulating biomarker