Next Article in Journal
Using the Colloidal Method to Prepare Au Catalysts for the Alkylation of Aniline by Benzyl Alcohol
Next Article in Special Issue
The Brain Metabolome Is Modified by Obesity in a Sex-Dependent Manner
Previous Article in Journal
The Cytoprotective Role of Autophagy in Response to BRAF-Targeted Therapies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Epigenetic Marks as Predictors of Metabolic Response to Bariatric Surgery: Validation from an Epigenome Wide Association Study

by
Carolina Gutiérrez-Repiso
1,2,3,*,
Antonio Cantarero-Cuenca
4,
Andrés González-Jiménez
4,
Teresa Linares-Pineda
1,2,
Nerea Peña-Montero
1,2,
Luis Ocaña-Wilhelmi
5,6,
Francisco J. Tinahones
1,2,3,7,† and
Sonsoles Morcillo
1,2,3,*,†
1
Unidad de Gestión Clínica de Endocrinología y Nutrición, Hospital Clínico Universitario Virgen de la Victoria, Campus de Teatinos s/n, 29010 Málaga, Spain
2
Instituto de Investigación Biomédica de Málaga y Plataforma en Nanomedicina-IBIMA Plataforma BIONAND, Campus de Teatinos s/n, 29010 Málaga, Spain
3
Centro de Investigación Biomédica en Red de Fisiopatología de la Obesidad y la Nutrición (CIBERobn), Instituto de Salud Carlos III, 28029 Madrid, Spain
4
ECAI Bioinformática, Instituto de Investigación Biomédica de Málaga y Plataforma en Nanomedicina-IBIMA Plataforma BIONAND, 29590 Málaga, Spain
5
Unidad de Gestión Clínica de Cirugía General, Digestiva y Trasplantes, Hospital Universitario Virgen de la Victoria, 29590 Málaga, Spain
6
Departamento de Especialidades Quirúrgicas, Bioquímica e Inmunología, Universidad de Málaga, 29590 Málaga, Spain
7
Departamento de Medicina y Dermatología, Universidad de Málaga, 29590 Málaga, Spain
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(19), 14778; https://doi.org/10.3390/ijms241914778
Submission received: 30 August 2023 / Revised: 26 September 2023 / Accepted: 28 September 2023 / Published: 30 September 2023
(This article belongs to the Special Issue Molecular Insights in Obesity and Metabolism)

Abstract

:
Little is known about the potential role of epigenetic marks as predictors of the resolution of obesity-related comorbidities after bariatric surgery. In this study, 20 patients were classified according to the metabolic improvement observed 6 months after sleeve gastrectomy, based on the diagnosis of metabolic syndrome, into responders if metabolic syndrome reversed after bariatric surgery (n = 10) and non-responders if they had metabolic syndrome bariatric surgery (n = 10). Blood DNA methylation was analyzed at both study points using the Infinium Methylation EPIC Bead Chip array-based platform. Twenty-six CpG sites and their annotated genes, which were previously described to be associated with metabolic status, were evaluated. Cg11445109 and cg19469447 (annotated to Cytochrome P450 2E1 (CYP2E1) gene) were significantly more hypomethylated in the responder group than in the non-responder group at both study points, whilst cg25828445 (annotated to Nucleolar Protein Interacting With The FHA Domain Of MKI67 Pseudogene 3 (NIFKP3) gene) showed to be significantly more hypermethylated in the non-responder group compared to the responder group at both study points. The analysis of the methylation sites annotated to the associated genes showed that CYP2E1 had 40% of the differentially methylated CpG sites, followed by Major Histocompatibility Complex, Class II, DR Beta 1 (HLA-DRB1) (33.33%) and Zinc Finger Protein, FOG Family Member 2 (ZFPM2) (26.83%). Cg11445109, cg19469447 and cg25828445 could have a role in the prediction of metabolic status and potential value as biomarkers of response to bariatric surgery.

1. Introduction

The concept of metabolic syndrome encompasses a wide range of metabolic abnormalities, including abdominal obesity, insulin resistance, hypertension and hyperlipidaemia, that result in a high risk of developing cardiovascular disease [1]. Although the number of different definitions of metabolic syndrome leads to discrepancies in its prevalence rates, the incidence and prevalence of metabolic syndrome go hand in hand with that of obesity and type 2 diabetes [2].
Lifestyle and environmental factors, in addition to genetic and epigenetic factors, are involved in the pathogenesis of metabolic syndrome [3]. Epigenetic mechanisms have been proposed to be a relevant factor involved in the development of metabolic diseases including metabolic syndrome, with epigenetics emerging as a promising tool for the management of metabolic diseases.
One of the most common epigenetic alterations associated with metabolic diseases is DNA methylation. Promoter methylation has usually been associated with a decrease in gene expression, whilst intragenic methylation has been related to increased gene expression [4].
Previous studies have investigated the relationship between epigenetic marks, DNA methylation-based, and metabolic syndrome as a whole or with separate components of the metabolic syndrome [5,6].
Furthermore, DNA methylation patterns respond to environmental factors, including weight loss intervention. Bariatric surgery is the most effective and durable treatment for weight loss and the amelioration of obesity-related comorbidities, including metabolic syndrome. Epigenetic factors have been suggested to be involved in the mechanisms related to metabolic improvement after bariatric surgery and the impact of bariatric surgery on DNA methylation levels has been shown in different tissues such as blood, adipose tissue, skeletal muscle and liver, among others [7]. Although DNA methylation changes have been associated with metabolic improvement observed after weight loss interventions [8,9], most of the studies have focused on evaluating the reversion of obesity-related epigenetic marks [10], and only a few studies have analyzed the potential role of epigenetic marks as predictors of the weight loss achievement and resolution of obesity-related comorbidities after bariatric surgery [11,12], although these studies have evaluated the methylation profiles of selected genes, warranting more studies to investigate the role of epigenetic marks in the prediction of bariatric surgery outcomes.
In a previous longitudinal study of the methylation status in a metabolically healthy obese population with a 11-year follow-up, we showed 26 differentially methylated CpG sites that could have a predictive role of the stability/worsening of metabolically healthy obese phenotype at long term [13]. The aim of the present study is to validate these epigenetics marks in a group of patients undergoing bariatric surgery and their association with the metabolic improvement observed after bariatric surgery, in order to identify those epigenetic marks that could predict the metabolic response to the bariatric surgery.

2. Results

The anthropometric and biochemical variables of the patients included in the study are shown in Table 1. Before surgery, triglyceride levels were significantly higher in the MetS non-responder group. There were no statistically significant differences between groups in the rest of variables studied at baseline. As was expected, patients improved their metabolic status 6 months after surgery; statistically significant differences between groups are shown in Table 1.

2.1. Validation of Differentially Methylated CpG Sites

Of the 26 significantly differentially methylated CpG sites previously described [13], three of them (cg11445109, cg19469447, cg25828445) were shown to be differentially methylated at baseline and 6 months after surgery. Both cg11445109 and cg19469447 (annotated to CYP2E1 gene) were hypomethylated, with their methylation levels significantly lower in the MetS responder group compared to the MetS non-responder group both before surgery and six months after bariatric surgery, whilst cg25828445 (annotated to NIFKP3 gene) was shown to be hypermethylated, with significantly higher methylation levels in the MeS non-responder group compared to the MetS responder group, both before surgery and 6 months after bariatric surgery (Figure 1). The methylation levels of the 26 CpG sites are shown in Supplementary Figure S1.
In a previous study based in a metabolically healthy obese population with a 11-year follow-up [13], both cg11445109 and cg19469447 were also hypomethylated, with their methylation levels being significantly lower in patients with metabolically healthy phenotype during 11-year follow-up compared to those who transitioned to unhealthy phenotype during the follow-up, whilst cg25828445 was significantly hypermethylated in those obese patients who transitioned to an unhealthy phenotype (Figure 1).

2.2. Differentially Methylated Genes

These 26 CpG sites were related to a total of 17 genes and two pseudogenes. A more in-depth analysis was performed on these genes. We collected all CpG sites in each of these genes from the UCSC genome browser and checked whether they were detected in the Methylation EPIC Bead. The detected CpG sites described in each gene were analyzed to investigate to what extent these genes presented multiple differentially methylated CpG sites. Fourteen of the nineteen genes identified showed multiple differentially methylated CpG sites. The gene with the highest percentage of methylated sites was CYP2E1 (40%), followed by HLA-DRB1 (33.33%) and ZFPM2 (26.83%) (Table 2).

2.3. Associations between DNA Methylation and Other Anthropometric and Biochemical Variables

We further analyzed the associations between anthropometric and biochemical variables and the methylation levels of differentially methylated CpG sites at both study points. Before surgery, cg11445109 and cg19469447 (annotated to CYP2E1 gene) were correlated positively with hip circumference (r = 0.48; p = 0.04 and r = 0.48; p = 0.004, respectively) and diastolic blood pressure (r = 0.49; p = 0.003 and r = 0.47; p = 0.004, respectively). Six months after surgery, cg11445109 and cg19469447 (CYP2E1 gene) were correlated positively with triglycerides (r = 0.61; p = 0.003 and r = 0.62; p = 0.003, respectively), cholesterol (r = 0.58; p = 0.007 and r = 0.65; p = 0.001, respectively) and LDL-cholesterol levels (r = 0.48; p = 0.03 and r = 0.54; p = 0.01, respectively). Cg25828445 (annotated to NIFKP3 gene) showed no association with the anthropometric and biochemical variables included in the study (Figure 2).
The associations found with the rest of the differentially methylated CpG sites from each gene are shown in Figure 2. It should be highlighted that CpG sites annotated to CYP2E1 (cg00321709, cg13315147, cg21024264 and cg24530264) correlated positively with triglycerides, cholesterol and LDL-cholesterol 6 months after bariatric surgery (Figure 2).

3. Discussion

In a previous study [13], we found 26 CpG sites associated with the stability/worsening of the metabolically healthy phenotype in a group of patients with obesity followed up for 11 years. In the present study, we studied these CpG sites in a group of patients with metabolic syndrome who underwent bariatric surgery and analyzed their association with metabolic outcomes observed after surgery. In this study, we confirm the potential role of cg11445109, cg19469447 and cg25828445 in the metabolic homeostasis and their association with metabolic health.
Our results showed that cg11445109 and cg19469447 (annotated to CYP2E1 gene) were significantly more hypomethylated in patients who improved their metabolic status compared to those who remained with metabolic syndrome after bariatric surgery. Indeed, CYP2E1 showed multiple differentially methylated sites in patients who improved their metabolic status after bariatric surgery compared to those who did not. In a previous study [13], we showed that the hypomethylation of both sites was associated with the stability of metabolically healthy phenotype in a group of patients followed up for 11 years.
CYP2E1 is a member of the CYP superfamily, mainly expressed in hepatocytes, which is involved in the metabolism of several low-molecular weight xenobiotics such as ethanol or benzene, and endogenous compounds such as fatty acids and bile acids [14]. CYP2E1 has been shown to be a relevant contributor to ROS production. In animal models, its expression has been associated with an increase in oxidative stress and apoptosis [15].
CYP2E1 has been related to metabolic diseases; its activity has been shown to be increased in patients with T2DM and obesity [16,17] and several polymorphisms have been significantly associated with the risk of T2DM [18].
In animal models, CYP2E1 deletion has been shown to induce adipose tissue browning, to increase whole-body energy expenditure and to improve glucose homeostasis [19]. CYP2E1 deficiency has been suggested to activate PPARα, which is involved in lipid homeostasis and energy regulation [20]. Moreover, our results show that several CpG sites annotated to CYP2E1 were correlated with markers of lipid levels at both study points. Previous studies have shown that overexpression of CYP2E1 is associated with plasma lipid levels and liver steatosis [21]. These findings point out CYP2E1 as a promising target to develop new strategies to fight against obesity and related diseases.
Our results show that cg25828445, annotated to NIFKP3 gene, was significantly hypermethylated in patients who remained with metabolic syndrome after bariatric surgery. In a previous study [13], we showed that cg25828445 was significantly hypermethylated in those patients who transitioned to metabolically unhealthy phenotype in a group of patients followed up for 11 years.
To the best of our knowledge, there is little information about NIFKP3. NIFKP3 is a processed pseudogene that is not encoded into a functional protein. Processed pseudogenes are a result of the reintegration of a gene due to retrotransposition. Pseudogenes have been suggested to have a regulatory role in the cells, and RNA transcripts of pseudogenes have been suggested to play a role as siRNA, miRNA-like and antisense RNA, which can regulate target gene expression [22].
The search for non-invasive biomarkers based on blood DNA methylation has increased in recent decades, becoming a huge challenge in metabolic diseases such as metabolic syndrome, a complex condition multifactorial in origin that involves multiple tissues and organs that are inaccessible without surgery. The search of biomarkers based on DNA methylation is an intricate task. DNA methylation is widely known to be influenced by genetic, physiological, and environmental factors. The gut microbiota, which has been shown to play an important role in host metabolism, and the use of chronic medication (i.e., metformin) are some of the factors whose actions have been suggested to be mediated by epigenetic modifications, which deserves further investigation.
This study highlights the potential role of DNA methylation as a predictive tool for metabolic outcomes after bariatric surgery. Further research is warranted to clarify the potential use of DNA methylation patterns as non-invasive biomarkers to predict metabolic improvement associated with bariatric surgery [11,23]. The identification of epigenetic biomarkers that might be useful in the prediction of metabolic improvement observed after bariatric surgery could help obesity management and development of appropriate treatment strategies.
Our study also presents some limitations. We could not perform RNA analysis to associate DNA methylation status to gene expression. Another limitation of the study is the small simple size. The use of peripherical blood could be considered a limitation of this study, although the main aim of this study is the identification of potential biomarkers. Whole blood is characterized by a high cellular heterogeneity, and methylation is a tissue-specific molecular mechanism. However, when it comes to the prediction role of a biomarker, this does not have to be involved in the causation of a disease, and may not have a mechanistic role in the disease.
However, as a strength of the study, we have confirmed the association between three epigenetic marks (cg11445109, cg19469447, cg25828445) and the metabolic status in a different population. Previously, we have shown the potential use of these CpG sites as markers of stability of metabolically healthy phenotype sites in patients with obesity followed up for 11 years with no weight loss intervention. On the other hand, in the present study, we have validated these CpG sites as potential markers of bariatric surgery response in patients with morbid obesity.

4. Material and Methods

4.1. Subjects

This study was carried out at the Virgen de la Victoria University Hospital (Málaga, Spain). The inclusion criteria were patients with morbid obesity and metabolic syndrome who underwent laparoscopic sleeve gastrectomy procedure between January 2020 and December 2020 and attended both visits (before surgery and 6 months after surgery). Metabolic syndrome was established using the National Cholesterol Education Program (NCEP) Adult Treatment Panel III (NCEP ATPIII criteria) [24]. Patients were excluded if they had cardiovascular disease, acute inflammatory disease, infectious disease, or had complications after the surgery. From the patients who met the criteria, we selected the first 10 patients who were no longer patients with metabolic syndrome 6 months after the surgery (MetS responder group) and the first 10 patients who continued with metabolic syndrome 6 months after the surgery (MetS non-responder group).
Weight and height measurements were made before surgery and 6 months after surgery. Body mass index (BMI) was calculated as weight (kg)/height2 (m2). Blood pressure was measured twice using a sphygmomanometer with an interval of 5 min between measurements, and the average of the two measurements was calculated.
Blood samples were collected after a 10–12 h fast before surgery and 6 months after surgery. Serum was separated and immediately frozen at −80 °C. Blood glucose was measured using the glucose oxidase method (Bayer, Leverkusen, Germany). Total cholesterol, triglycerides and high-density lipoprotein cholesterol levels were measured using enzymatic methods. Insulin was determined using radioimmunoassay (DIASource ImmunoAssays SA, Nivelles, Belgium). The homeostasis model assessment of insulin resistance (HOMA-IR) was calculated as follows: fasting insulin (µIU/mL) × fasting glucose (mmol/L)/22.5.

4.2. DNA Methylation Assay

DNA was extracted from peripheral blood using the QIAmp DNA Blood Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. DNA concentration was quantified with a Qubit 3.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA) using Qubit dsDNA HS Assay Kit Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA) After quantification, a total of 500 ng of genomic DNA was bisulfite-treated using a Zymo EZ-96 DNA Methylation™ Kit (Zymo Research Corp, Irvine, CA, USA) and was purified using a DNA-Clean-Up Kit (Zymo Research Corp, Irvine, CA, USA).
Over 850,000 methylation sites were interrogated with the Infinium Methylation EPIC Bead Chip Kit (Illumina, San Diego, CA, USA) following the Infinium HD Assay Methylation protocol, and raw data were obtained from iS (Illumina) software (v2.0).

4.3. Methylation Data Analysis

We used statistical programming language R 3.5.1 accessed on 1 March 2021 (https://www.r-project.org) to perform the methylation data analysis. Raw data files (idat files) were read with the minfi package [25] to calculate raw β-values. Normal-exponential out-of-band (NOOB) normalization [26] was used to correct the background. Probes located at sexual chromosomes or near SNPs were removed from the analysis. Low quality probes (those with a detection p-value > 0.01 in at least 10% of samples) were also removed.
Finally, beta-mixture quantile (BMIQ) normalization [27] was applied to correct for the two different bead designs in the microarrays. For the differential methylation analysis, we transformed β-values to M-values.
In the current study, we evaluate the methylation levels of those previously identified CpG sites and neighboring CpG sites related to the same genes (Supplementary Table S1) [13].

4.4. Statistical Analysis

Differences in anthropometric and biochemical variables were analyzed with the Mann–Whitney U test for continuous data and the chi-square test for categorial data using R software (3.5.1). Data were expressed as the mean ± standard deviation. Values were statistically significant when p < 0.05.
Differences between groups for each CpG site were analyzed using the Mann–Whitney U test. Correlations between DNA methylation levels, and anthropometric and biochemicals variables were analyzed using Spearman’s rank correlation coefficients.

5. Conclusions

In conclusion, we confirmed that cg11445109, cg19469447 and cg25828445 sites could have a potential role in the prediction of metabolic status. Identifying non-invasive blood-based biomarkers for metabolic syndrome may aid in diagnosing and screening those patients at risk of developing metabolic diseases. Importantly, these biomarkers may help in the management and treatment of patients with metabolic disease.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241914778/s1.

Author Contributions

S.M. and F.J.T. designed and coordinated the study. F.J.T. and L.O.-W. selected the patients. L.O.-W. performed bariatric surgery procedure. T.L.-P. and N.P.-M. collected and processed the samples and conducted research. A.C.-C. and A.G.-J. implemented the bioinformatic workflow and analyzed the data. C.G.-R. analyzed the data and wrote the manuscript. S.M. and F.J.T. were involved in writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

CGR is supported by Miguel Servet program from Instituto de Salud Carlos III (CP20/00066), TMLP was supported by a grant from Instituto de Salud Carlos III (FI19/00178). NPM is supported by a grant from “Programa Estatal de Promoción del Talento y su empleabilidad 2021” (POEJ-00039-18-IBIMA). SM is supported by Nicolas Monardes program from the Consejería de Salud de la Junta de Andalucía (RC-0008-2021). This work was supported in part by a grant from the Instituto de Salud Carlos III (PI15-01350). This study has been co-funded by FEDER funds (“A way to make Europe”). Funding for open access charge: Universidad de Málaga/CBUA. CIBER Fisiopatología de la Obesidad y Nutrición (CIBEROBN) is ISCIII project.

Institutional Review Board Statement

The study was performed in accordance with the ethical standards of the Declaration of Helsinki and was reviewed and approved by the Biomedical Research Ethic Coordinator Committee of Andalucía (CCEIBA) (Project code: PI15-01350, date of approval: 22 March 2018).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data sets used and analyzed in this study are available from the corresponding authors on request.

Conflicts of Interest

Authors declare that they have no conflict of interest.

References

  1. Virani, S.S.; Alonso, A.; Benjamin, E.J.; Bittencourt, M.S.; Callaway, C.W.; Carson, A.P.; Chamberlain, A.M.; Chang, A.R.; Cheng, S.; Delling, F.N.; et al. Heart Disease and Stroke Statistics-2020 Update: A Report from the American Heart Association. Circulation 2020, 141, E139–E596. [Google Scholar] [CrossRef] [PubMed]
  2. Fahed, G.; Aoun, L.; Bou Zerdan, M.; Allam, S.; Bou Zerdan, M.; Bouferraa, Y.; Assi, H.I. Metabolic Syndrome: Updates on Pathophysiology and Management in 2021. Int. J. Mol. Sci. 2022, 23, 786. [Google Scholar] [CrossRef] [PubMed]
  3. Ussar, S.; Griffin, N.W.; Bezy, O.; Fujisaka, S.; Vienberg, S.; Softic, S.; Deng, L.; Bry, L.; Gordon, J.I.; Kahn, C.R. Interactions between Gut Microbiota, Host Genetics and Diet Modulate the Predisposition to Obesity and Metabolic Syndrome. Cell Metab. 2015, 22, 516–530. [Google Scholar] [CrossRef] [PubMed]
  4. Jones, P.A. Functions of DNA methylation: Islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 2012, 13, 484–492. [Google Scholar] [CrossRef]
  5. Nuotio, M.-L.; Pervjakova, N.; Joensuu, A.; Karhunen, V.; Hiekkalinna, T.; Milani, L.; Kettunen, J.; Järvelin, M.-R.; Jousilahti, P.; Metspalu, A.; et al. An epigenome-wide association study of metabolic syndrome and its components. Sci. Rep. 2020, 10, 20567. [Google Scholar] [CrossRef]
  6. Akinyemiju, T.; Do, A.N.; Patki, A.; Aslibekyan, S.; Zhi, D.; Hidalgo, B.; Tiwari, H.K.; Absher, D.; Geng, X.; Arnett, D.K.; et al. Epigenome-wide association study of metabolic syndrome in African-American adults. Clin. Epigenetics 2018, 10, 49. [Google Scholar] [CrossRef]
  7. ElGendy, K.; Malcomson, F.C.; Bradburn, D.M.; Mathers, J.C. Effects of bariatric surgery on DNA methylation in adults: A systematic review and meta-analysis. Surg. Obes. Relat. Dis. 2020, 16, 128–136. [Google Scholar] [CrossRef]
  8. Faenza, M.; Benincasa, G.; Docimo, L.; Nicoletti, G.F.; Napoli, C. Clinical epigenetics and restoring of metabolic health in severely obese patients undergoing batriatric and metabolic surgery. Updates Surg. 2022, 74, 431–438. [Google Scholar] [CrossRef]
  9. Fraszczyk, E.; Luijten, M.; Spijkerman, A.M.W.; Snieder, H.; Wackers, P.F.K.; Bloks, V.W.; Nicoletti, C.F.; Nonino, C.B.; Crujeiras, A.B.; Buurman, W.A.; et al. The effects of bariatric surgery on clinical profile, DNA methylation, and ageing in severely obese patients. Clin. Epigenetics 2020, 12, 14. [Google Scholar] [CrossRef]
  10. Morcillo, S.; Macías-González, M.; Tinahones, F.J. The Effect of Metabolic and Bariatric Surgery on DNA Methylation Patterns. Curr. Atheroscler. Rep. 2017, 19, 40. [Google Scholar] [CrossRef]
  11. Nicoletti, C.F.; Nonino, C.B.; de Oliveira, B.A.P.; Pinhel, M.A.d.S.; Mansego, M.L.; Milagro, F.I.; Zulet, M.A.; Martinez, J.A. DNA Methylation and Hydroxymethylation Levels in Relation to Two Weight Loss Strategies: Energy-Restricted Diet or Bariatric Surgery. Obes. Surg. 2016, 26, 603–611. [Google Scholar] [CrossRef] [PubMed]
  12. Coppedè, F.; Seghieri, M.; Stoccoro, A.; Santini, E.; Giannini, L.; Rossi, C.; Migliore, L.; Solini, A. DNA methylation of genes regulating appetite and prediction of weight loss after bariatric surgery in obese individuals. J. Endocrinol. Investig. 2019, 42, 37–44. [Google Scholar] [CrossRef] [PubMed]
  13. Gutiérrez-Repiso, C.; Linares-Pineda, T.M.; Gonzalez-Jimenez, A.; Aguilar-Lineros, F.; Valdés, S.; Soriguer, F.; Rojo-Martínez, G.; Tinahones, F.J.; Morcillo, S. Epigenetic Biomarkers of Transition from Metabolically Healthy Obesity to Metabolically Unhealthy Obesity Phenotype: A Prospective Study. Int. J. Mol. Sci. 2021, 22, 10417. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, J.; Jiang, S.; Wang, J.; Renukuntla, J.; Sirimulla, S.; Chen, J. A comprehensive review of cytochrome P450 2E1 for xenobiotic metabolism. Drug Metab. Rev. 2019, 51, 178–195. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, W.; Lu, D.; Dong, W.; Zhang, L.; Zhang, X.; Quan, X.; Ma, C.; Lian, H.; Zhang, L. Expression of CYP2E1 increases oxidative stress and induces apoptosis of cardiomyocytes in transgenic mice. FEBS J. 2011, 278, 1484–1492. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Z.; Hall, S.D.; Maya, J.F.; Li, L.; Asghar, A.; Gorski, J.C. Diabetes mellitus increases the in vivo activity of cytochrome P450 2E1 in humans. Br. J. Clin. Pharmacol. 2003, 55, 77–85. [Google Scholar] [CrossRef]
  17. van Rongen, A.; Välitalo, P.A.J.; Peeters, M.Y.M.; Boerma, D.; Huisman, F.W.; van Ramshorst, B.; van Dongen, E.P.A.; Anker, J.N.v.D.; Knibbe, C.A.J. Morbidly Obese Patients Exhibit Increased CYP2E1-Mediated Oxidation of Acetaminophen. Clin. Pharmacokinet. 2016, 55, 833–847. [Google Scholar] [CrossRef]
  18. Zhang, L.; Tang, J.; Wang, Y.; Wang, X.; Wang, F. Association of CYP7A1 and CYP2E1 Polymorphisms with Type 2 Diabetes in the Chinese Han Populations. Pharmgenomics Pers. Med. 2022, 15, 843–855. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Zhang, L.; Xu, P.; Qin, X.; Wang, P.; Cheng, Y.; Yao, B.; Wang, X. Cytochrome P450 2E1 gene knockout or inhibition prevents obesity induced by high-fat diet via regulating energy expenditure. Biochem. Pharmacol. 2022, 202, 115160. [Google Scholar] [CrossRef]
  20. Zhang, Y.; Yan, T.; Wang, T.; Liu, X.; Hamada, K.; Sun, D.; Sun, Y.; Yang, Y.; Wang, J.; Takahashi, S.; et al. Crosstalk between CYP2E1 and PPAR α substrates and agonists modulate adipose browning and obesity. Acta Pharm. Sin. B 2022, 12, 2224–2238. [Google Scholar] [CrossRef]
  21. Kathirvel, E.; Morgan, K.; French, S.W.; Morgan, T.R. Overexpression of liver-specific cytochrome P4502E1 impairs hepatic insulin signaling in a transgenic mouse model of nonalcoholic fatty liver disease. Eur. J. Gastroenterol. Hepatol. 2009, 21, 973–983. [Google Scholar] [CrossRef] [PubMed]
  22. Witek, J.; Mohiuddin, S.S. Biochemistry, Pseudogenes; StatPearls: St. Petersburg, FL, USA, 2022. [Google Scholar]
  23. Keller, M.; Meir, A.Y.; Bernhart, S.H.; Gepner, Y.; Shelef, I.; Schwarzfuchs, D.; Tsaban, G.; Zelicha, H.; Hopp, L.; Müller, L.; et al. DNA methylation signature in blood mirrors successful weight-loss during lifestyle interventions: The CENTRAL trial. Genome Med. 2020, 12, 18. [Google Scholar] [CrossRef]
  24. Cleeman, J.I. Executive summary of the third report of the National Cholesterol Education Program (NCEP) expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (adult treatment panel III). J. Am. Med. Assoc. 2001, 285, 2486–2497. [Google Scholar] [CrossRef]
  25. Aryee, M.J.; Jaffe, A.E.; Corrada-Bravo, H.; Ladd-Acosta, C.; Feinberg, A.P.; Hansen, K.D.; Irizarry, R.A. Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 2014, 30, 1363–1369. [Google Scholar] [CrossRef] [PubMed]
  26. Triche, T.J.; Weisenberger, D.J.; Van Den Berg, D.; Laird, P.W.; Siegmund, K.D. Low-level processing of Illumina Infinium DNA Methylation BeadArrays. Nucleic Acids Res. 2013, 41, e90. [Google Scholar] [CrossRef] [PubMed]
  27. Teschendorff, A.E.; Marabita, F.; Lechner, M.; Bartlett, T.; Tegner, J.; Gomez-Cabrero, D.; Beck, S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics 2013, 29, 189–196. [Google Scholar] [CrossRef]
Figure 1. Methylation levels of the differentially methylated CpG sites. MetS responder: patients with metabolic syndrome who improved their metabolic status 6 months after bariatric surgery. MetS non-responder: patients with metabolic syndrome 6 months after surgery. Stable MHO: metabolically healthy obese patients during 11-year follow up. Unstable MHO: metabolically healthy obese patients who developed metabolic complications at 11-year follow-up.
Figure 1. Methylation levels of the differentially methylated CpG sites. MetS responder: patients with metabolic syndrome who improved their metabolic status 6 months after bariatric surgery. MetS non-responder: patients with metabolic syndrome 6 months after surgery. Stable MHO: metabolically healthy obese patients during 11-year follow up. Unstable MHO: metabolically healthy obese patients who developed metabolic complications at 11-year follow-up.
Ijms 24 14778 g001
Figure 2. Heatmap correlation plot between methylation levels and anthropometric and biochemical variables at both study points.
Figure 2. Heatmap correlation plot between methylation levels and anthropometric and biochemical variables at both study points.
Ijms 24 14778 g002
Table 1. Anthropometric and biochemical variables of the patients included in the study.
Table 1. Anthropometric and biochemical variables of the patients included in the study.
MetS Non-Responder at BaselineMetS Responder at BaselineMetS Non-Responder 6 m after SurgeryMetS Responder 6 m after Surgery
Sex (M/F)1/93/7
Age (years)49.7 ± 10.2947.7 ± 5.49
Weight (kg)126.45 ± 22.91125.32 ± 20.56101.29 ± 15.9194.59 ± 12.35
BMI (kg/m2)51.54 ± 10.3445.9 ± 5.0341.29 ± 7.2534.73 ± 2.82 #
Glucose (mg/dL)113.4 ± 18.77117.1 ± 28.8596.60 ± 14.8589.1 ± 7.92
Insulin (µUI/mL)19.84 ± 9.0127.91 ± 18.6310.80 ± 2.7811.39 ± 11.18
HOMA-IR5.6 ± 2.958.21 ± 5.552.53 ± 0.632.54 ± 2.61
C-peptide (ng/mL)4.59 ± 1.474.31 ± 1.373.14 ± 0.592.52 ± 0.81 #
Cholesterol (mg/dL)188.3 ± 25.18175.6 ± 37.76207.50 ± 37.24174.3 ± 29.77 #
Triglycerides (mg/dL)197.3 ± 45.59137.4 ± 70.45 *153.60 ± 21.0681.30 ± 19.39 #
HDL-cholesterol (mg/dL)38.6 ± 4.1945.8 ± 8.9543.20 ± 3.7358.10 ± 12.83 #
LDL-cholesterol (mg/dL)111.18 ± 25.37104.92 ± 31.39131.26 ± 35.2699.94 ± 24.88
AST (U/L)28.5 ± 19.6925.4 ± 13.3817.56 ± 7.1215.70 ± 3.77
ALT (U/L)44.7 ± 26.0639.3 ± 22.8127.70 ± 8.9021.40 ± 6.62
GGT (U/L)46.8 ± 31.8345.1 ± 17.4636.00 ± 27.5826.00 ± 14.53
SBP (mm Hg)145.11 ± 24.4137.8 ± 21.17141.0 ± 17.45124.4 ± 14.51 #
DBP (mm Hg)87.11 ± 13.1483.4 ± 12.6387.78 ± 11.5578.6 ± 11.85
Diabetes treatment (y/n)7/37/34/61/9
HTA treatment (y/n)7/38/27/34/6
GGT: Gamma-glutamyl transferase. AST: aspartate aminotransferase. ALT: alanine aminotransferase. SBP: systolic blood pressure. DBP: diastolic blood pressure. * p < 0.05 between MetS non-responder and MetS responder group at baseline. # p < 0.05 between MetS non-responder and MetS responder group 6 months after surgery.
Table 2. CpG sites differentially methylated in the identified genes.
Table 2. CpG sites differentially methylated in the identified genes.
GeneKEGG OrthologyDifferentially Methylated CpG SitesTotal Number of CpG Sites%
RAP1GAP2Signaling and cellular processes218823.86
MUC2Membrane trafficking116018.33
ZFPM2Transcription factor114126.83
COG5Membrane trafficking86312.70
CYP2E1Metabolism82040.00
PLCZ1Membrane trafficking62920.69
SLC1A1Signaling and cellular processes62524.00
FRMD1Signaling and cellular processes54211.90
HLA-DRB1Signaling and cellular processes51533.33
ZBTB45Transcription factor52718.52
HLA-DQ2Signaling and cellular processes4498.16
HOOK3Signaling and cellular processes. Membrane trafficking43511.43
DENND1AMembrane trafficking31816.67
TNFRSF4Signaling and cellular processes21315.38
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gutiérrez-Repiso, C.; Cantarero-Cuenca, A.; González-Jiménez, A.; Linares-Pineda, T.; Peña-Montero, N.; Ocaña-Wilhelmi, L.; Tinahones, F.J.; Morcillo, S. Epigenetic Marks as Predictors of Metabolic Response to Bariatric Surgery: Validation from an Epigenome Wide Association Study. Int. J. Mol. Sci. 2023, 24, 14778. https://doi.org/10.3390/ijms241914778

AMA Style

Gutiérrez-Repiso C, Cantarero-Cuenca A, González-Jiménez A, Linares-Pineda T, Peña-Montero N, Ocaña-Wilhelmi L, Tinahones FJ, Morcillo S. Epigenetic Marks as Predictors of Metabolic Response to Bariatric Surgery: Validation from an Epigenome Wide Association Study. International Journal of Molecular Sciences. 2023; 24(19):14778. https://doi.org/10.3390/ijms241914778

Chicago/Turabian Style

Gutiérrez-Repiso, Carolina, Antonio Cantarero-Cuenca, Andrés González-Jiménez, Teresa Linares-Pineda, Nerea Peña-Montero, Luis Ocaña-Wilhelmi, Francisco J. Tinahones, and Sonsoles Morcillo. 2023. "Epigenetic Marks as Predictors of Metabolic Response to Bariatric Surgery: Validation from an Epigenome Wide Association Study" International Journal of Molecular Sciences 24, no. 19: 14778. https://doi.org/10.3390/ijms241914778

APA Style

Gutiérrez-Repiso, C., Cantarero-Cuenca, A., González-Jiménez, A., Linares-Pineda, T., Peña-Montero, N., Ocaña-Wilhelmi, L., Tinahones, F. J., & Morcillo, S. (2023). Epigenetic Marks as Predictors of Metabolic Response to Bariatric Surgery: Validation from an Epigenome Wide Association Study. International Journal of Molecular Sciences, 24(19), 14778. https://doi.org/10.3390/ijms241914778

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop