Differential methylation in EGLN1 associates with blood oxygen saturation and plasma protein levels in high-altitude pulmonary edema

BACKGROUND: High-altitude (HA, 2500 m) hypoxic exposure evokes a multitude of physiological processes. The hypoxia-sensing genes though influence transcriptional output in disease susceptibility; the exact regulatory mechanisms remain undetermined in high-altitude pulmonary edema (HAPE). Here, we investigated the differential DNA methylation distribution in the two genes encoding the oxygen-sensing HIF-prolyl hydroxylases, prolyl hydroxylase domain protein 2 (PHD2) and factor inhibiting HIF-1α and the consequent contributions to the HAPE pathophysiology. METHODS: Deep sequencing of the sodium bisulfite converted DNA segments of the two genes, Egl nine homolog 1 (EGLN1) and Hypoxia Inducible Factor 1 Subunit Alpha Inhibitor (HIF1AN), was conducted to analyze the differential methylation distribution in three study groups, namely HAPE-patients (HAPE-p), HAPE-free sojourners (HAPE-f) and healthy HA natives (HLs). HAPE-p and HAPE-f were permanent residents of low altitude (< 200 m) of North India who traveled to Leh (3500 m), India, and were recruited through Sonam Norboo Memorial (SNM) hospital, Leh. HLs were permanent residents of altitudes at and above 3500 m. In addition to the high resolution, bisulfite converted DNA sequencing, gene expression of EGLN1 and HIF1AN and their plasma protein levels were estimated. RESULTS: A significantly lower methylation distribution of CpG sites was observed in EGLN1 and higher in HIF1AN (P < 0.01) in HAPE-p compared to the two control groups, HAPE-f and HLs. Of note, differential methylation distribution of a few CpG sites, 231,556,748, 231,556,804, 231,556,881, 231,557,317 and 231,557,329, in EGLN1 were significantly associated with the risk of HAPE (OR = 4.79–10.29; P = 0.048–004). Overall, the methylation percentage in EGLN1 correlated with upregulated plasma PHD2 levels (R = − 0.36, P = 0.002) and decreased peripheral blood oxygen saturation (SpO(2)) levels (R = 0.34, P = 0.004). We also identified a few regulatory SNPs in the DNA methylation region of EGLN1 covering chr1:231,556,683–231,558,443 suggestive of the functional role of differential methylation distribution of these CpG sites in the regulation of the genes and consequently in the HIF-1α signaling. CONCLUSIONS: Significantly lower methylation distribution in EGLN1 and the consequent physiological influences annotated its functional epigenetic relevance in the HAPE pathophysiology. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s13148-022-01338-z.