Genome-wide analysis reflects novel 5-hydroxymethylcytosines implicated in diabetic nephropathy and the biomarker potential

Abstract

Aim: Diabetic nephropathy (DN) has become the most common cause of end-stage renal disease in most countries for patients with type 2 diabetes (T2D). Elucidating novel epigenetic contributors to DN can not only enhance our understanding of this complex disorder but also lay the foundation for developing more effective monitoring tools and preventive interventions in the future, thus contributing to our ultimate goal of improving patient care.

Methods: 5-hydroxymethylcytosines (5hmC)-Seal, a highly selective chemical labeling technique, was used to profile genome-wide 5hmC, a stable cytosine modification type marking gene activation, in circulating cell-free DNA (cfDNA) samples from a cohort of patients recruited at Zhongnan Hospital, including T2D patients with nephropathy (DN, n = 12), T2D patients with non-DN vascular complications (non-DN, n = 29), and T2D patients without any complication (controls, n = 14). Differential analysis was performed to find DN-associated 5hmC features, followed by the exploration of biomarker potential of 5hmC in cfDNA for DN using a machine learning approach.

Results: Genome-wide analyses of 5hmC in cfDNA detected 427 and 336 differential 5hmC modifications associated with DN, compared with non-DN individuals and controls, and suggested relevant pathways such as NOD-like receptor signaling pathway and tyrosine metabolism. Our exploration using a machine learning approach revealed an exploratory model comprised of ten 5hmC genes showing the possibility to distinguish DN from non-DN individuals or controls.

Conclusion: Genome-wide analysis suggests the possibility of exploiting novel 5hmC in patient-derived cfDNA as a non-invasive tool for monitoring DN in high-risk T2D patients in the future.

Keywords

Type 2 diabetes / nephropathy / epigenetics / 5-hydroxymethylcytosine / cfDNA

References

[1] Chen J,Liu D.Exosomes: advances, development and potential therapeutic strategies in diabetic nephropathy.Metabolism2021;122:154834

[2] American Diabetes Association. 11. Microvascular complications and foot care: standards of medical care in diabetes-2021.Diabetes Care2021;44:S151-67

[3] Lecamwasam A,Saffery R.Potential for novel biomarkers in diabetes-associated chronic kidney disease: epigenome, metabolome, and gut microbiome.Biomedicines2020;8:341 PMCID:PMC7555227

[4] Bakdash K,Annam A,Kondo K.Complications of percutaneous renal biopsy.Semin Intervent Radiol2019;36:97-103 PMCID:PMC6531025

[5] Zhang W.Towards clinical implementation of circulating cell-free DNA in precision medicine.J Transl Genet Genom2019;3:11 PMCID:PMC6774639

[6] Zampieri M,Barchetta I.Increased PARylation impacts the DNA methylation process in type 2 diabetes mellitus.Clin Epigenetics2021;13:114 PMCID:PMC8130175

[7] Zeng C,Wang J,Hou L.Application of the high-throughput TAB-Array for the discovery of novel 5-hydroxymethylcytosine biomarkers in pancreatic ductal adenocarcinoma.Epigenomes2019;3:16 PMCID:PMC6693877

[8] Yuan EF,Cheng L.Hyperglycemia affects global 5-methylcytosine and 5-hydroxymethylcytosine in blood genomic DNA through upregulation of SIRT6 and TETs.Clin Epigenetics2019;11:63 PMCID:PMC6466651

[9] Zeng C,Zhang Z,Zhang W.Towards precision medicine: advances in 5-hydroxymethylcytosine cancer biomarker discovery in liquid biopsy.Cancer Commun (Lond)2019;39:12 PMCID:PMC6440138

[10] Song CX,Fu Y.Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine.Nat Biotechnol2011;29:68-72 PMCID:PMC3107705

[11] Li W,Lu X.5-Hydroxymethylcytosine signatures in circulating cell-free DNA as diagnostic biomarkers for human cancers.Cell Res2017;27:1243-57 PMCID:PMC5630683

[12] Cai J,Zhang Z.Genome-wide mapping of 5-hydroxymethylcytosines in circulating cell-free DNA as a non-invasive approach for early detection of hepatocellular carcinoma.Gut2019;68:2195-205 PMCID:PMC6872444

[13] Cai J,Hua W.An integrative analysis of genome-wide 5-hydroxymethylcytosines in circulating cell-free DNA detects noninvasive diagnostic markers for gliomas.Neurooncol Adv2021;3:vdab049 PMCID:PMC8209591

[14] Chiu BC,You Q.Alterations of 5-hydroxymethylation in circulating cell-free DNA reflect molecular distinctions of subtypes of non-Hodgkin lymphoma.NPJ Genom Med2021;6:11 PMCID:PMC7878492

[15] Chiu BC,You Q.Prognostic implications of 5-hydroxymethylcytosines from circulating cell-free DNA in diffuse large B-cell lymphoma.Blood Adv2019;3:2790-9 PMCID:PMC6784517

[16] Yang Y,Lu X.5-Hydroxymethylcytosines in Circulating Cell-Free DNA Reveal Vascular Complications of Type 2 Diabetes.Clin Chem2019;65:1414-25 PMCID:PMC7055673

[17] Han L,Lu X.Alterations of 5-hydroxymethylcytosines in circulating cell-free DNA reflect retinopathy in type 2 diabetes.Genomics2021;113:79-87

[18] American Diabetes Association. Standards of medical care in diabetes-2017: summary of revisions.Diabetes Care2017;40:S4-5

[19] Spanopoulos D,Zaccardi F.Temporal variation of renal function in people with type 2 diabetes mellitus: a retrospective UK clinical practice research datalink cohort study.Diabetes Obes Metab2019;21:1817-23 PMCID:PMC6767485

[20] Han D,Shih AH.A highly sensitive and robust method for genome-wide 5hmc profiling of rare cell populations.Mol Cell2016;63:711-9 PMCID:PMC4992443

[21] Harrow J,Gonzalez JM.GENCODE: the reference human genome annotation for the ENCODE project.Genome Res2012;22:1760-74 PMCID:PMC3431492

[22] Liao Y,Shi W.FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features.Bioinformatics2014;30:923-30

[23] Kundaje A,Ernst J.Roadmap epigenomics consortiumintegrative analysis of 111 reference human epigenomes.Nature2015;518:317-30 PMCID:PMC4530010

[24] Shannon P,Ozier O.Cytoscape: a software environment for integrated models of biomolecular interaction networks.Genome Res2003;13:2498-504 PMCID:PMC403769

[25] Wu T,Xu S.ClusterProfiler 4.0: a universal enrichment tool for interpreting Omics data. I.nnovation (N Y)2021;2:100141 PMCID:PMC8454663

[26] Szklarczyk D,Nastou KC.The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets.Nucleic Acids Res2021;49:D605-12 PMCID:PMC7779004

[27] Subramanian A,Mootha VK.Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proc Natl Acad Sci102:15545-50 PMCID:PMC1239896

[28] Yu G,Han Y.ClusterProfiler: an R package for comparing biological themes among gene clusters.OMICS2012;16:284-7 PMCID:PMC3339379

[29] Kanehisa M,Sato Y,Tanabe M.KEGG for integration and interpretation of large-scale molecular data sets.Nucleic Acids Res2012;40:D109-14 PMCID:PMC3245020

[30] Zou H.Regularization and variable selection via the elastic net.J Royal Statistical Soc B2005;67:301-20

[31] Tang YL,Zeng ZG.Gene expression-based analysis identified NTNG1 and HGF as biomarkers for diabetic kidney disease.Medicine (Baltimore)2020;99:e18596 PMCID:PMC6946191

[32] Dafinger C,Borgal L.Targeted deletion of the AAA-ATPase Ruvbl1 in mice disrupts ciliary integrity and causes renal disease and hydrocephalus. E.xp Mol Med2018;50:1-17 PMCID:PMC6026120

[33] Huang F,Guo F.FoxO1-mediated inhibition of STAT1 alleviates tubulointerstitial fibrosis and tubule apoptosis in diabetic kidney disease.EBioMedicine2019;48:491-504 PMCID:PMC6838438

[34] Anders HJ.NOD-like and Toll-like receptors or inflammasomes contribute to kidney disease in a canonical and a non-canonical manner.Kidney Int2013;84:225-8

[35] Qu W,Li M,Li L.Revealing the underlying mechanism of diabetic nephropathy viewed by microarray analysis.Exp Clin Endocrinol Diabetes2015;123:353-9

[36] Nunemaker CS,Verrilli GM,Upadhye A.Increased serum CXCL1 and CXCL5 are linked to obesity, hyperglycemia, and impaired islet function.J Endocrinol2014;222:267-76 PMCID:PMC4135511

[37] Geng XD,Feng Z.Identification of key genes and pathways in diabetic nephropathy by bioinformatics analysis.J Diabetes Investig2019;10:972-84 PMCID:PMC6626994

[38] Zhang Z,Zhang W.Considerations before the application of 5-hydroxymethylation levels of long non-coding RNAs for non-invasive cancer diagnosis.Extracell Vesicles Circ Nucl Acids2022;3:10-3 PMCID:PMC8797161

[39] Moen EL,Mu W.Genome-wide variation of cytosine modifications between European and African populations and the implications for complex traits.Genetics2013;194:987-96 PMCID:PMC3730924