SUMMARY
OBJECTIVE: Juvenile idiopathic arthritis is the most common chronic rheumatic disease of autoimmune origin in childhood. While altered DNA methylation has been associated with adult autoimmune rheumatic diseases, studies on juvenile idiopathic arthritis remain limited. This study evaluated levels of 5-methylcytosine and 5-hydroxymethylcytosine, important epigenetic biomarkers, in juvenile idiopathic arthritis and its subtypes.
METHODS: This single-center cross-sectional study included patients diagnosed with juvenile idiopathic arthritis aged 1 month to 18 years. Children with systemic or autoimmune diseases and comorbidities were excluded. The study comprised 58 juvenile idiopathic arthritis patients and 35 age-gender-matched healthy controls presenting with non-rheumatological complaints. Demographic and clinical characteristics, including juvenile idiopathic arthritis subtypes, were recorded. The 5-methylcytosine and 5-hydroxymethylcytosine levels were assessed quantitatively using an enzyme-linked immunosorbent assay kit. The relationships between 5-methylcytosine and 5-hydroxymethylcytosine and juvenile idiopathic arthritis and its subtypes were evaluated statistically.
RESULTS: The 5-hydroxymethylcytosine level was significantly higher, whereas the 5-methylcytosine and 5-methylcytosine/5-hydroxymethylcytosine ratios were significantly lower in the patient group than in the control group. No significant differences between different juvenile idiopathic arthritis subtypes and disease activity status were found in 5-hydroxymethylcytosine and 5-methylcytosine levels or 5-methylcytosine/5-hydroxymethylcytosine ratios (p>0.05). The 5-methylcytosine/5-hydroxymethylcytosine ratio, with an optimal cut-off value of ≤28.05, demonstrated a sensitivity of 70.69% and a specificity of 77.14% in distinguishing juvenile idiopathic arthritis patients from controls (p<0.001).
CONCLUSION: 5-Hydroxymethylcytosine and 5-methylcytosine can be potential biomarkers for juvenile idiopathic arthritis and its subtypes, indicating their roles in juvenile idiopathic arthritis development regarding molecular mechanisms and novel therapeutic approaches.
KEYWORDS:
Juvenile idiopathic arthritis; Epigenetic; 5-mC; 5-methylcytosine; DNA methylation; 5-hydroxymethyl cytosine; ELISA; DNA
INTRODUCTION
Juvenile idiopathic arthritis (JIA) is a complex childhood autoimmune rheumatic disease that affects children aged 6 months to 16 years1. Literature data suggest that several immune system components, such as inflammatory cells, including neutrophils, plasma cells, dendritic cells, and activated T-cells, are involved in the pathogenesis of JIA. Although the exact cause of JIA is unknown, family studies indicate that it is a complex disease in which environmental and genetic variables play significant roles2.
Global DNA methylation via the insertion of a methyl at position 5 of cytosine in cytosine-phosphate-guanine (CpG) dinucleotides, resulting in the formation of 5-methylcytosine (5-mC), which affects gene expression, chromatin remodeling, and genome integrity, might have a role in understanding physiological mechanisms and disease conditions3,4. Transformation of 5-mC into 5-hydroxymethylcytosine (5-hmC) may influence various intracellular functions5. Imbalances in 5-mC and 5-hmC markers cause transcriptional dysregulation of enhancers and promoters, altering the cell's fate6. Aberrant DNA methylation is also implicated in several autoimmune diseases, including rheumatoid arthritis7.
This study aimed to investigate the 5-mC and 5-hmC levels in JIA patients and healthy controls and to analyze the impact of JIA subtypes and disease activity on methylation biomarkers.
METHODS
Study design
The study was designed as a single-center cross-sectional study. The protocol was approved by the Non-Interventional Clinical Research Ethics Committee of Samsun Training and Research Hospital (GOKAEK 2024/23/16). The study was conducted according to the ethical considerations outlined in the Declaration of Helsinki. Informed consent was obtained from the parents or legal guardians of the children included in the study. Children's consent was also obtained to ensure that they understood the purpose and procedures of the study.
Population and sample
The study population comprised patients between 1 month and 18 years diagnosed with JIA according to the International League of Associations for Rheumatology criteria8. Patient data were retrospectively collected from two centers: the Department of Pediatric Rheumatology at Başakşehir Çam and Sakura City Hospital in Istanbul and Samsun University Training and Research Hospital in Samsun, Turkey. The study was conducted under ethical approval obtained from Samsun Training and Research Hospital. Children with any systemic or autoimmune diseases and comorbidities were excluded from the study. In the end, the patient group consisted of 58 patients with JIA, and the control group consisted of 35 age- and gender-matched healthy children who presented to the outpatient clinics of the same hospital with complaints unrelated to rheumatological conditions.
Data collection
Children's demographic and clinical characteristics, including disease duration, disease activity status, JIA subtype classification, and biological treatments and inflammatory marker levels (C-reactive protein [CRP] and erythrocyte sedimentation rate [ESR]), were recorded.
DNA isolation
Three milliliters of blood were collected from all participants, regardless of disease status, into tubes containing ethylenediaminetetraacetic acid upon their admission to the pediatric rheumatology outpatient clinics. Blood samples from healthy controls were obtained during the same study period.
Genomic DNA was extracted from blood samples taken from the participants using ELK Biotechnology's Blood/Cell/Tissue Genomic DNA Extraction Kit (Cat no: EP007), following the manufacturer's instructions. The extracted DNA samples were stored at −20°C.
5-Methylcytosine measurements
For 5-mC global DNA methylation, the Abcam Global DNA Methylation Assay Kit (5-Methyl Cytosine, Colorimetric, cat no: ab233486) was used. The rate of methylated DNA (5-mC) in the total DNA was calculated using the following formula: 5-mC%=sample's optical density (OD)-negative control's OD×100%: slope ×S.
5-Hydroxymethylcytosine measurements
For 5-hmC global DNA hydroxymethylation, the Abcam Global DNA Hydroxymethylation Assay Kit (5-hydroxymethyl Cytosine, Colorimetric, cat no: ab233487) was used. The percentage of hydroxymethylated DNA (5-hmC) in the total DNA was calculated using the following formula: 5-hmC%=sample's OD-negative control's OD×100%: slope ×S.
Statistical analysis
Statistical analyses were conducted using the JAMOVI software (Version 2.3.28, The Jamovi Project, 2023, https://www.jamovi.org) and JASP (Version 0.19.0, Jeffreys’ Amazing Statistics Program, 2024, https://jasp-stats.org). Given the sample sizes in our JIA (n=58) and control (n=35) groups, we employed the Kolmogorov-Smirnov test for the JIA group and the Shapiro-Wilk test for the control group to assess normality distribution, complemented by visual examination of histograms and Q-Q plots.
The results of the statistical analyses were expressed using descriptive statistics, that is, median values for quantitative variables and numbers for categorical variables.
The differences between the 5-mC global methylation and 5-hmC global hydroxymethylation results of the patient and control groups were assessed using the independent samples median test and Fisher's exact test, where necessary. The relationships between global methylation, hydroxymethylation levels, and clinical parameters were evaluated using the Spearman rank correlation (ρ) test. Probability (p) statistics of <0.05 were deemed to indicate statistical significance.
To evaluate the 5-mC/5-hmC ratio's diagnostic potential in distinguishing JIA patients from controls, receiver operating characteristic (ROC) curve analysis was performed using DeLong's method. The calculation of area under the curve (AUC), sensitivity, specificity, and optimal cut-off value was determined by the Youden index (J=sensitivity+specificity—1), where a larger value indicates better diagnostic performance.
Height, weight, and body mass index (BMI) Z-scores were calculated using the 2000 Centers for Disease Control and Prevention (CDC) growth charts as reference values.
RESULTS
The distribution of the demographic and clinical characteristics by the patient and control groups is shown in Table 1. Eight patients (13.8%) had active disease at the time of blood sampling. Accordingly, the most common (34.5%) JIA subtype was oligoarticular JIA, and uveitis was present in four (6.9%) patients. In terms of biological treatments JIA patients received, adalimumab was used most frequently (58.6%), followed by etanercept (20.7%) and tocilizumab (10.3%). No statistically significant differences were observed between patient and control groups regarding height Z-scores (p=0.210), weight Z-scores (p=0.183), and BMI Z-scores (p=0.157) based on CDC reference values.
The distribution of biomarker and inflammatory parameters by the patient and control groups is shown in Table 1. The 5-hmC level was significantly higher, whereas the 5-mC level and 5-mC/5-hmC ratio were significantly lower in the patient group compared to controls (p<0.001, p=0.038, and p<0.001, respectively) (Figure 1). Regarding inflammatory markers, CRP level was significantly higher in patients (p<0.001), while ESR showed no significant difference between groups (p>0.05).
Comparison of DNA methylation and hydroxymethylation levels between juvenile idiopathic arthritis patients and healthy controls. Box plots showing (A) 5-hydroxymethylcytosine levels, (B) 5-methylcytosine levels, and (C) 5-methylcytosine/5-hydroxymethylcytosine ratios in juvenile idiopathic arthritis patients (n=58) versus healthy controls (n=35). Juvenile idiopathic arthritis patients demonstrated significantly higher 5-hydroxymethylcytosine levels (p<0.001), significantly lower 5-methylcytosine levels (p=0.038), and significantly lower 5-methylcytosine/5-hydroxymethylcytosine ratios (p<0.001) compared to controls. Box plots display median values with interquartile ranges; individual data points represent outliers. Statistical significance was determined using independent samples median test.
There was no significant difference in 5-hmC and 5-mC levels and, thus, in the 5-mC/5-hmC ratio between patients with different JIA subtypes (p>0.05) (Table 2).
Comparison of 5-methylcytosine, 5-hydroxymethylcytosine levels, and 5-methylcytosine/5-hydroxymethylcytosine ratio according to disease activity in juvenile idiopathic arthritis patients.
As shown in Table 2, no significant differences were observed in 5-hmC levels, 5-mC levels, or 5-mC/5-hmC ratios between patients with active and inactive JIA (p>0.05 for all comparisons).
The ROC curve analysis revealed that the 5-mC/5-hmC ratio demonstrated good diagnostic performance in distinguishing JIA patients from controls, with an AUC of 0.760 (95%CI 0.660–0.842, p<0.001), sensitivity of 70.69%, and specificity of 77.14% at the optimal cut-off value of ≤28.05.
DISCUSSION
JIA involves synovial tissue containing various inflammatory cells, including neutrophils, plasma cells, dendritic cells, and most activated T cells2. It has been suggested that the highly activated Th1-type phenotype of T cells in chronically inflamed joints of children with JIA may reflect specific recruitment events that contribute to the polarization of these cells.9. These T-cells produce tumor necrosis factor beta, interleukin 2, and interferon-γ. Th1-dominant responses have been linked to several other autoimmune conditions, including psoriasis, inflammatory bowel disease, rheumatoid arthritis (RA), and type 1 diabetes mellitus2. Although the etiology of autoimmune disorders involves a complex genetic predisposition10, environmental variables also play significant roles in the development of autoimmune diseases, as demonstrated in genetically identical monozygotic twins11.
Examples of epigenetic mechanisms and elements affected by hereditary and environmental variables include DNA methylation, post-transcriptional modifications of histones, and microribonucleic acids12. These mechanisms and elements control gene transcription rates, modify chromatin architecture, and govern accessibility to transcriptional regulatory factors. The most extensively researched and thoroughly described epigenetic alteration mechanism is DNA methylation. The DNA methyltransferase (DNMT)-mediated enzymatic transfer of a methyl group from S-adenosyl-L-methionine to the CpG dinucleotide results in in the formation of 5-mC. The CpG dinucleotides, which are not typical and severely methylated in the bulk genome, are typically unmethylated in CpG island clusters, which are collections of CpG dinucleotides. Methylation of CpG islands in gene promoters is linked to the long-term suppression of gene expression. On the other hand, tissue-specific and variable methylation preferentially occurs outside of CpG islands. In a study investigating the transcript levels of types DNMT1, DNMT3A, and DNMT3B in peripheral blood mononuclear cells to assess the importance of epigenetic changes in the pathogenesis of JIA, Ghavidel et al. found that DNMT1 and DNMT3A expressions were significantly lower in JIA patients than in healthy control subjects (7- and 5.5-fold, respectively), and this difference was even more significant in the case of young male JIA patients13.
Several inflammatory and autoimmune disorders have been linked to global DNA hypomethylation, which involves abnormal gene and ribosomal RNA expression that likely plays a role in the disease's pathophysiology14. Autoimmune diseases have not been linked to other DNA hypomethylation-related alterations, such as genomic instability, mutations, or the use of cryptic promoters15. It has been demonstrated that systemic lupus erythematosus (SLE) cluster of differentiation 4 (CD4+) T cells have global DNA hypomethylation, contributing to the over-reactivity of T cells in SLE patients16,17. SLE CD4+ T cells have lower levels of DNMT1, an enzyme that maintains DNA methylation, indicating that passive DNA demethylation contributes to the pathophysiology of SLE18. Fibroblast-like synoviocytes (FLS) in RA have been shown to exhibit DNA hypomethylation. Overexpression of genes essential to the disease process has been linked to hypomethylation at particular CpG sites in FLS. The active phenotype typical of FLS develops following 5-azacytidine-induced demethylation, highlighting the effects of DNA hypomethylation in FLS. However, there are still many unresolved aspects of DNA methylation in RA. DNA methylation in the gene encoding the pro-inflammatory cytokine interleukin-32 is decreased in CD4+ T cells of JIA patients7. All these findings suggest that unmethylated cytokine production increases and upregulates immune system hyperactivity, potentially leading to autoimmune diseases. Hence, decreased DNA methylation is likely linked to the pathogenesis of JIA.
Based on these literature data, we investigated the peripheral blood 5-mC and 5-hmC levels in JIA patients and the relationships between DNA methylation and hydroxymethylation and JIA subtypes. Consequently, we found that the 5-hmC level was significantly higher. In contrast, the 5-mC level and, thus, the 5-mC/5-hmC ratio were significantly lower in the patients with JIA than in healthy control subjects, suggesting that epigenetic changes play a role in JIA pathogenesis.
Our ROC curve analysis highlighted the potential diagnostic utility of the 5-mC/5-hmC ratio as a novel biomarker for JIA, demonstrating good discriminatory power between patients and healthy controls. We observed balanced sensitivity and specificity, suggesting its potential value as a supplementary diagnostic tool. These findings are particularly noteworthy given the current challenges in early JIA diagnosis and the need for objective biomarkers. The relatively high specificity suggests that the 5-mC/5-hmC ratio could be especially valuable in confirming disease when used with traditional clinical and laboratory parameters. However, the moderate sensitivity indicates that this marker alone may not be sufficient for screening. Future longitudinal studies should investigate whether this epigenetic ratio could also serve as a predictive marker for disease progression or treatment response, potentially contributing to more personalized therapeutic approaches in JIA management.
We also found that these epigenetic modifications were associated with elevated CRP levels. However, we did not find any significant relationship between these epigenetic changes and JIA subtypes, which may be attributed to the insufficient patients in the JIA subgroups.
The limitations of the study
The primary limitation of our study was that the cross-sectional design may not have allowed for the establishment of causal relationships between epigenetic changes and disease progression. Secondly, insufficient patients in the JIA subgroups may have limited our ability to detect JIA subtype-specific methylation differences. Our study does not exclude the need to investigate gene-specific methylation alterations. Lastly, the lack of pre-treatment methylation data and the differences between the treatment regimens received by patients prevented us from drawing any conclusions about the impact of current medical treatments on methylation patterns. Given these limitations of our study, caution is advised when interpreting its results.
CONCLUSION
Although the etiology of JIA is not fully known, the interaction between genes and environment is likely to be important, as in other autoimmune diseases. In this context, a critical finding of our study is that the 5-hmC levels were elevated in patients with JIA, highlighting a key aspect of epigenetic mechanisms influencing inflammatory processes in JIA. In conclusion, 5-hmC and 5-mC can be potential biomarkers of JIA and its subtypes. The study's findings regarding the roles of 5-hmC and 5-mC in the development of JIA pave the way for further studies to elucidate the pathogenesis of JIA and develop new treatment options.
DATA AVAILABILITY STATEMENT
The datasets generated and/or analyzed during the current study are available from the corresponding author upon reasonable request.
REFERENCES
-
1 Chavez-Valencia RA, Chiaroni-Clarke RC, Martino DJ, Munro JE, Allen RC, Akikusa JD, et al. The DNA methylation landscape of CD4+ T cells in oligoarticular juvenile idiopathic arthritis. J Autoimmun. 2018;86:29-38. https://doi.org/10.1016/j.jaut.2017.09.010
» https://doi.org/10.1016/j.jaut.2017.09.010 -
2 Prahalad S, Glass DN. A comprehensive review of the genetics of juvenile idiopathic arthritis. Pediatr Rheumatol Online J. 2008;6:11. https://doi.org/10.1186/1546-0096-6-11
» https://doi.org/10.1186/1546-0096-6-11 -
3 Xu W, Wang F, Yu Z, Xin F. Epigenetics and cellular metabolism. Genet Epigenet. 2016;8:43-51. https://doi.org/10.4137/GEG.S32160
» https://doi.org/10.4137/GEG.S32160 -
4 Steenbergen RD, Snijders PJ, Heideman DA, Meijer CJ. Clinical implications of (epi)genetic changes in HPV-induced cervical precancerous lesions. Nat Rev Cancer. 2014;14(6):395-405. https://doi.org/10.1038/nrc3728
» https://doi.org/10.1038/nrc3728 -
5 Wang J, Tang J, Lai M, Zhang H. 5-Hydroxymethylcytosine and disease. Mutat Res Rev Mutat Res. 2014;762:167-75. https://doi.org/10.1016/j.mrrev.2014.09.003
» https://doi.org/10.1016/j.mrrev.2014.09.003 -
6 Laird A, Thomson JP, Harrison DJ, Meehan RR. 5-hydroxymethylcytosine profiling as an indicator of cellular state. Epigenomics. 2013;5(6):655-69. https://doi.org/10.2217/epi.13.69
» https://doi.org/10.2217/epi.13.69 -
7 Meyer B, Chavez RA, Munro JE, Chiaroni-Clarke RC, Akikusa JD, Allen RC, et al. DNA methylation at IL32 in juvenile idiopathic arthritis. Sci Rep. 2015;5:11063. https://doi.org/10.1038/srep11063
» https://doi.org/10.1038/srep11063 - 8 Petty RE, Southwood TR, Manners P, Baum J, Glass DN, Goldenberg J, et al. International League of Associations for Rheumatology classification of juvenile idiopathic arthritis: second revision, Edmonton, 2001. J Rheumatol. 2004;31(2):390-2. PMID: 14760812
- 9 Wedderburn LR, Robinson N, Patel A, Varsani H, Woo P. Selective recruitment of polarized T cells expressing CCR5 and CXCR3 to the inflamed joints of children with juvenile idiopathic arthritis. Arthritis Rheum. 2000;43(4):765-74. https://doi.org/10.1002/1529-0131(200004)43:4<765::AID-ANR7>3.0.CO;2-B
-
10 Floreani A, Leung PS, Gershwin ME. Environmental basis of autoimmunity. Clin Rev Allergy Immunol. 2016;50(3):287-300. https://doi.org/10.1007/s12016-015-8493-8
» https://doi.org/10.1007/s12016-015-8493-8 -
11 Ballestar E. Epigenetics lessons from twins: prospects for autoimmune disease. Clin Rev Allergy Immunol. 2010;39(1):30-41. https://doi.org/10.1007/s12016-009-8168-4
» https://doi.org/10.1007/s12016-009-8168-4 -
12 Liberman N, Wang SY, Greer EL. Transgenerational epigenetic inheritance: from phenomena to molecular mechanisms. Curr Opin Neurobiol. 2019;59:189-206. https://doi.org/10.1016/j.conb.2019.09.012
» https://doi.org/10.1016/j.conb.2019.09.012 -
13 Ghavidel AA, Shiari R, Hassan-Zadeh V, Farivar S. The expression of DNMTs is dramatically decreased in peripheral blood mononuclear cells of male patients with juvenile idiopathic arthritis. Allergol Immunopathol (Madr). 2020;48(2):182-6. https://doi.org/10.1016/j.aller.2019.08.003
» https://doi.org/10.1016/j.aller.2019.08.003 - 14 Javierre BM, Hernando H, Ballestar E. Environmental triggers and epigenetic deregulation in autoimmune disease. Discov Med. 2011;12(67):535-45. PMID: 22204770
-
15 Nakano K, Whitaker JW, Boyle DL, Wang W, Firestein GS. DNA methylome signature in rheumatoid arthritis. Ann Rheum Dis. 2013;72(1):110-7. https://doi.org/10.1136/annrheumdis-2012-201526
» https://doi.org/10.1136/annrheumdis-2012-201526 -
16 Wu H, Zhao M, Chang C, Lu Q. The real culprit in systemic lupus erythematosus: abnormal epigenetic regulation. Int J Mol Sci. 2015;16(5):11013-33. https://doi.org/10.3390/ijms160511013
» https://doi.org/10.3390/ijms160511013 -
17 Zhang Y, Zhao M, Sawalha AH, Richardson B, Lu Q. Impaired DNA methylation and its mechanisms in CD4(+)T cells of systemic lupus erythematosus. J Autoimmun. 2013;41:92-9. https://doi.org/10.1016/j.jaut.2013.01.005
» https://doi.org/10.1016/j.jaut.2013.01.005 -
18 Lei W, Luo Y, Lei W, Luo Y, Yan K, Zhao S, et al. Abnormal DNA methylation in CD4+ T cells from patients with systemic lupus erythematosus, systemic sclerosis, and dermatomyositis. Scand J Rheumatol. 2009;38(5):369-74. https://doi.org/10.1080/03009740902758875
» https://doi.org/10.1080/03009740902758875
Edited by
-
Scientifıc Editor:
José Maria Soares Júnior https://orcid.org/0000-0003-0774-9404


