RESUMO
Our study aimed to identify differentially methylated regions (i.e., genomic region where multiple adjacent CpG sites show differential methylation) and their enriched genomic pathways associated with knee osteoarthritis pain (KOA). We recruited cognitively healthy middle to older aged (age 45-85) adults with (n = 182) and without (n = 31) self-reported KOA pain. We also extracted DNA from peripheral blood that was analyzed using MethylationEPIC arrays. The R package minfi (Aryee et al., 2014) was used to perform methylation data preprocessing and quality control. To investigate biological pathways impacted by differential methylation, we performed pathway enrichment analysis using Ingenuity Pathway Analysis (IPA) to identify canonical pathways and upstream regulators. Annotated genes within ± 5 kb of the putative differentially methylated regions (DMRs, p < 0.05) were subjected to the IPA analysis. There was no significant difference in age, sex, study site between no pain and pain group (p > 0.05). Non-Hispanic black individuals were overrepresented in the pain group (p = 0.003). At raw p < 0.05 cutoff, we identified a total of 19,710 CpG probes, including 13,951 hypermethylated CpG probes, for which DNA methylation level was higher in the groups with highest pain grades. We also identified 5,759 hypomethylated CpG probes for which DNA methylation level was lower in the pain groups with higher pain grades. IPA revealed that pain-related DMRs were enriched across multiple pathways and upstream regulators. The top 10 canonical pathways were linked to cellular signaling processes related to immune responses (i.e., antigen presentation, PD-1, PD-L1 cancer immunotherapy, B cell development, IL-4 signaling, Th1 and Th2 activation pathway, and phagosome maturation). Moreover, in terms of upstream regulators, NDUFAF3 was the most significant (p = 8.6E-04) upstream regulator. Our findings support previous preliminary work suggesting the importance of epigenetic regulation of the immune system in knee pain and the need for future work to understand the epigenetic contributions to chronic pain.
RESUMO
Our study aimed to identify differentially methylated CpGs/regions and their enriched genomic pathways associated with underlying chronic musculoskeletal pain in older individuals. We recruited cognitively healthy older adults with (n = 20) and without (n = 9) self-reported musculoskeletal pain and collected DNA from peripheral blood that was analyzed using MethylationEPIC arrays. We identified 31,739 hypermethylated CpG and 10,811 hypomethylated CpG probes (ps ≤ 0.05). All CpG probes were clustered into 5966 regions, among which 600 regions were differentially methylated at p ≤ 0.05 level, including 294 hypermethylated regions and 306 hypomethylated regions (differentially methylated regions). Ingenuity pathway enrichment analysis revealed that the pain-related differentially methylated regions were enriched across multiple pathways. The top 10 canonical pathways were linked to cellular signaling processes related to immune responses (i.e. antigen presentation, programed cell death 1 receptor/PD-1 ligand 1, interleukin-4, OX40 signaling, T cell exhaustion, and apoptosis) and gamma-aminobutyric acid receptor signaling. Further, Weighted Gene Correlation Network Analysis revealed a comethylation network module in the pain group that was not preserved in the control group, where the hub gene was the cyclic adenosine monophosphate-dependent transcription factor ATF-2. Our preliminary findings provide new epigenetic insights into the role of aberrant immune signaling in musculoskeletal pain in older adults while further supporting involvement of dysfunctional GABAergic signaling mechanisms in chronic pain. Our findings need to be urgently replicated in larger cohorts as they may serve as a basis for developing and targeting future interventions.