Stochastic epigenetic mutation profiles as biomarkers of clinical activity in juvenile idiopathic arthritis: a multi-omic machine learning approach for gene prioritization.

Abstract

BACKGROUND: Juvenile idiopathic arthritis (JIA) is a rare autoimmune disease arising from a complex interplay between genetic and environmental factors. Epigenetic modifications such as DNA methylation (DNAm) have been described as potential mediators in gene-environment interactions, contributing to immune system dysregulation. Emerging evidence suggests that DNAm profiles also predict therapeutic responses in autoimmune diseases. This study aims to identify epigenetic biomarkers and epigenetic-driven gene expression changes associated with JIA clinical activity. METHODS: We reanalyzed a publicly available dataset of 44 JIA patients, with whole-genome DNAm and gene expression from CD4 + T cells measured at two points: at anti-TNF therapy withdrawal (T0) and eight months later (Tend). At Tend, 30 patients maintained inactive disease (ID) while 14 did not (NO ID). We investigated differences between ID and NO ID patients in the epigenetic mutation load and various epigenetic clocks through linear regression models, and prioritized genomic regions with significantly higher number of epimutations in NO ID patients through machine learning. RESULTS: We found a higher mutation load in NO ID than ID patients, both at T0 and at Tend, with the differences at Tend reaching statistical significance (p = 0.02). In contrast, we found no evidence of association between epigenetic clocks and JIA clinical activity. Using a multi-omic approach, we identified a List of candidate epigenetically-driven differentially expressed genes, 80 up-regulated and 77 down-regulated, in NO ID patients. Finally, comparing our candidate gene list with the Connectivity Map database, we identified new candidate potential therapeutic targets. Key findings were validated in independent datasets: DNAm profiles from CD4 + T cells (56 JIA patients, 57 controls) and transcriptomic data from PBMCs of JIA patients with active or inactive disease, confirming dysregulation of pathways such as TNF-α signaling via NF-kB and TGF-β signaling among others. CONCLUSIONS: We described a significant association of epigenetic mutations with JIA clinical activity, indicating that epigenetic changes might precede clinical symptoms and may serve as biomarkers for early disease monitoring. Further, our results shed light on biomolecular mechanisms of JIA, supporting the development of more effective treatments.