Soamy Montesino-Goicolea , Lingsong Meng , Asha Rani , Zhiguang Huo , Thomas C. Foster , Roger B. Fillingim , Yenisel Cruz-Almeida
{"title":"基于与膝关节骨关节炎疼痛相关的差异DNA甲基化谱的基因组途径的富集","authors":"Soamy Montesino-Goicolea , Lingsong Meng , Asha Rani , Zhiguang Huo , Thomas C. Foster , Roger B. Fillingim , Yenisel Cruz-Almeida","doi":"10.1016/j.ynpai.2022.100107","DOIUrl":null,"url":null,"abstract":"<div><p>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 <em>minfi</em> (<span>Aryee et al., 2014</span>) 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.</p></div>","PeriodicalId":52177,"journal":{"name":"Neurobiology of Pain","volume":"12 ","pages":"Article 100107"},"PeriodicalIF":0.0000,"publicationDate":"2022-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9755025/pdf/","citationCount":"4","resultStr":"{\"title\":\"Enrichment of genomic pathways based on differential DNA methylation profiles associated with knee osteoarthritis pain\",\"authors\":\"Soamy Montesino-Goicolea , Lingsong Meng , Asha Rani , Zhiguang Huo , Thomas C. Foster , Roger B. Fillingim , Yenisel Cruz-Almeida\",\"doi\":\"10.1016/j.ynpai.2022.100107\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><p>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 <em>minfi</em> (<span>Aryee et al., 2014</span>) 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.</p></div>\",\"PeriodicalId\":52177,\"journal\":{\"name\":\"Neurobiology of Pain\",\"volume\":\"12 \",\"pages\":\"Article 100107\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2022-08-01\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9755025/pdf/\",\"citationCount\":\"4\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Neurobiology of Pain\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S2452073X22000241\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"Medicine\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Neurobiology of Pain","FirstCategoryId":"1085","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S2452073X22000241","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"Medicine","Score":null,"Total":0}
引用次数: 4
摘要
我们的研究旨在确定差异甲基化区域(即多个相邻CpG位点显示差异甲基化的基因组区域)及其丰富的基因组通路与膝关节骨关节炎疼痛(KOA)相关。我们招募了认知健康的中老年(45-85岁)成年人,有(n = 182)和没有(n = 31)自我报告KOA疼痛。我们还从外周血中提取DNA,使用MethylationEPIC阵列进行分析。使用R包minfi (Aryee et al., 2014)进行甲基化数据预处理和质量控制。为了研究受差异甲基化影响的生物学途径,我们使用匠人途径分析(Ingenuity pathway analysis, IPA)进行了途径富集分析,以确定典型途径和上游调节因子。在假设的差异甲基化区域(DMRs, p <0.05)进行IPA分析。无疼痛组与疼痛组在年龄、性别、研究部位上无显著差异(p >0.05)。非西班牙裔黑人在疼痛组的比例过高(p = 0.003)。在raw p <在0.05截断值下,我们共鉴定出19,710个CpG探针,其中包括13,951个高甲基化的CpG探针,在疼痛等级最高的组中,DNA甲基化水平较高。我们还鉴定出5,759个低甲基化的CpG探针,其DNA甲基化水平在疼痛等级较高的疼痛组中较低。IPA显示,与疼痛相关的DMRs在多个通路和上游调节因子中富集。前10个典型途径与免疫应答相关的细胞信号传导过程有关(即抗原呈递、PD-1、PD-L1癌症免疫治疗、B细胞发育、IL-4信号传导、Th1和Th2激活途径以及吞噬体成熟)。在上游调节因子方面,NDUFAF3是最显著的上游调节因子(p = 8.6E-04)。我们的研究结果支持了先前的初步工作,表明免疫系统的表观遗传调控在膝关节疼痛中的重要性,以及未来工作需要了解表观遗传对慢性疼痛的贡献。
Enrichment of genomic pathways based on differential DNA methylation profiles associated with knee osteoarthritis pain
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.