Rheumatoid arthritis–associated DNA methylation sites in peripheral blood mononuclear cells

Epigenome-850K-wide profiling reveals peripheral blood differential methylation in term low birth weight

Aim: We investigate the genome-wide DNA methylation (DNAm) patterns of term low birth weight (TLBW) neonates. Methods: In the discovery phase, we assayed 32 samples (TLBW/control:16/16) using the EPIC 850k BeadChip Array. Targeted pyrosequencing of in 60 samples (TLBW/control:28/32) using targeted pyrosequencing during the replication phase. Results: The 850K array identified TLBW-associated 144 differentially methylated positions (DMPs) and 149 DMRs. Nearly 77% DMPs exhibited hypomethylation, located in the opensea and gene body regions. The most significantly enriched pathway in KEGG is sphingolipid metabolism (hsa00600), and the genes GALC and SGMS1 related to this pathway both show hypomethylation. Conclusion: Our analysis provides evidence of genome-wide DNAm alterations in TLBW. Further investigations are needed to elucidate the functional significance of these DNAm changes.

Identification of clinical characteristics biomarkers for rheumatoid arthritis through targeted DNA methylation sequencing

OBJECTIVES

Rheumatoid arthritis (RA) is a complex disease with a challenging diagnosis, especially in seronegative patients. The aim of this study is to investigate whether the methylation sites associated with the overall immune response in RA can assist in clinical diagnosis, using targeted methylation sequencing technology on peripheral venous blood samples.

METHODS

The study enrolled 241 RA patients, 30 osteoarthritis patients (OA), and 30 healthy volunteers control (HC). Fifty significant cytosine guanine (CG) sites between undifferentiated arthritis and RA were selected and analyzed using targeted DNA methylation sequencing. Logistic regression models were used to establish diagnostic models for different clinical features of RA, and six machine learning methods (logit model, random forest, support vector machine, adaboost, naive bayes, and learning vector quantization) were used to construct clinical diagnostic models for different subtypes of RA. Least absolute shrinkage and selection operator regression and detrended correspondence analysis were utilized to screen for important CGs. Spearman correlation was used to calculate the correlation coefficient.

RESULTS

The study identified 16 important CG sites, including tumor necrosis factort receptor associated factor 5 (TRAF5) (chr1:211500151), mothers against decapentaplegic homolog 3 (SMAD3) (chr15:67357339), tumor endothelial marker 1 (CD248) (chr11:66083766), lysosomal trafficking regulator (LYST) (chr1:235998714), PR domain zinc finger protein 16 (PRDM16) (chr1:3307069), A-kinase anchoring protein 10 (AKAP10) (chr17:19850460), G protein subunit gamma 7 (GNG7) (chr19:2546620), yes1 associated transcriptional regulator (YAP1) (chr11:101980632), PRDM16 (chr1:3163969), histone deacetylase complex subunit sin3a (SIN3A) (chr15:75747445), prenylated rab acceptor protein 2 (ARL6IP5) (chr3:69134502), mitogen-activated protein kinase kinase kinase 4 (MAP3K4) (chr6:161412392), wnt family member 7A (WNT7A) (chr3:13895991), inhibin subunit beta B (INHBB) (chr2:121107018), deoxyribonucleic acid replication helicase/nuclease 2 (DNA2) (chr10:70231628) and chromosome 14 open reading frame 180 (C14orf180) (chr14:105055171). Seven CG sites showed abnormal changes between the three groups (P < 0.05), and 16 CG sites were significantly correlated with common clinical indicators (P < 0.05). Diagnostic models constructed using different CG sites had an area under the receiver operating characteristic curve (AUC) range of 0.64-0.78 for high-level clinical indicators of high clinical value, with specificity ranging from 0.42 to 0.77 and sensitivity ranging from 0.57 to 0.88. The AUC range for low-level clinical indicators of high clinical value was 0.63-0.72, with specificity ranging from 0.48 to 0.74 and sensitivity ranging from 0.72 to 0.88. Diagnostic models constructed using different CG sites showed good overall diagnostic accuracy for the four subtypes of RA, with an accuracy range of 0.61-0.96, a balanced accuracy range of 0.46-0.94, and an AUC range of 0.46-0.94.

CONCLUSIONS

This study identified potential clinical diagnostic biomarkers for RA and provided novel insights into the diagnosis and subtyping of RA. The use of targeted deoxyribonucleic acid (DNA) methylation sequencing and machine learning methods for establishing diagnostic models for different clinical features and subtypes of RA is innovative and can improve the accuracy and efficiency of RA diagnosis.

Identification of Poly(ADP-ribose) Polymerase 9 (PARP9) as a Potent Suppressor for Mycobacterium tuberculosis Infection

ADP-ribosylation is a reversible and dynamic post-translational modification mediated by ADP-ribosyltransferases (ARTs). Poly(ADP-ribose) polymerases (PARPs) are an important family of human ARTs. ADP-ribosylation and PARPs have crucial functions in host-pathogen interaction, especially in viral infections. However, the functions and potential molecular mechanisms of ADP-ribosylation and PARPs in Mycobacterium infection remain unknown. In this study, bioinformatics analysis revealed significantly changed expression levels of several PARPs in tuberculosis patients compared to healthy individuals. Moreover, the expression levels of these PARPs returned to normal following tuberculosis treatment. Then, the changes in the expression levels of PARPs during Mycobacterium infection were validated in Tohoku Hospital Pediatrics-1 (THP1)-induced differentiated macrophages infected with Mycobacterium model strains bacillus Calmette-Guérin (BCG) and in human lung adenocarcinoma A549 cells infected with Mycobacterium smegmatis (Ms), respectively. The mRNA levels of PARP9, PARP10, PARP12, and PARP14 were most significantly increased during infection, with corresponding increases in protein levels, indicating the possible biological functions of these PARPs during Mycobacterium infection. In addition, the biological function of host PARP9 in Mycobacterium infection was further studied. PARP9 deficiency significantly increased the infection efficiency and intracellular proliferation ability of Ms, which was reversed by the reconstruction of PARP9. Collectively, this study updates the understanding of changes in PARP expression during Mycobacterium infection and provides evidence supporting PARP9 as a potent suppressor for Mycobacterium infection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00112-2.

DiseaseNet: a transfer learning approach to noncommunicable disease classification

As noncommunicable diseases (NCDs) pose a significant global health burden, identifying effective diagnostic and predictive markers for these diseases is of paramount importance. Epigenetic modifications, such as DNA methylation, have emerged as potential indicators for NCDs. These have previously been exploited in other contexts within the framework of neural network models that capture complex relationships within the data. Applications of neural networks have led to significant breakthroughs in various biological or biomedical fields but these have not yet been effectively applied to NCD modeling. This is, in part, due to limited datasets that are not amenable to building of robust neural network models. In this work, we leveraged a neural network trained on one class of NCDs, cancer, as the basis for a transfer learning approach to non-cancer NCD modeling. Our results demonstrate promising performance of the model in predicting three NCDs, namely, arthritis, asthma, and schizophrenia, for the respective blood samples, with an overall accuracy (f-measure) of 94.5%. Furthermore, a concept based explanation method called Testing with Concept Activation Vectors (TCAV) was used to investigate the importance of the sample sources and understand how future training datasets for multiple NCD models may be improved. Our findings highlight the effectiveness of transfer learning in developing accurate diagnostic and predictive models for NCDs.

Methylation changes of liver DNA during the formation of gallstones

Aim: To explore the overall methylation changes in liver tissues during the formation of gallstones, as well as the key pathways and genes involved in the process. Methods: Reduced-representation bisulfite sequencing and RNA sequencing were conducted on the liver tissues of mice with gallstones and control normal mice. Results: A total of 8705 differentially methylated regions in CpG and 1410 differentially expressed genes were identified. The joint analysis indicated that aberrant DNA methylation may be associated with dysregulated gene expression in key pathways such as cholesterol metabolism and bile secretion. Conclusion: We propose for the first time that methylation changes in some key pathway genes in liver tissue may be involved in the formation of gallstones.

Quantitative relationships of FAM50B and PTCHD3 methylation with reduced intelligence quotients in school aged children exposed to lead: Evidence from epidemiological and in vitro studies

At present, the application of DNA methylation (DNAm) biomarkers in environmental health risk assessment (EHRA) is more challenging due to the unclearly quantitative relationship between them. We aimed to explore the role of FAM50B and PTCHD3 at the level of signaling pathways, and establish the quantitative relationship between them and children's intelligence quotients (IQs). DNAm of target regions was measured in multiple cell models and was compared with the human population data. Then the dose-response relationships of lead exposure with neurotoxicity and DNAm were established by benchmark dose (BMD) model, followed by potential signaling pathway screening. Results showed that there was a quantitative linear relationship between children's IQs and FAM50B/PTCHD3 DNAm (DNAm between 51.40 % - 78.78 % and 31.41 % - 74.19 % for FAM50B and PTCHD3, respectively), and this relationship was more significant when children's IQs > 90. The receiver operating characteristic (ROC) and calibration curves showed that FAM50B/PTCHD3 DNAm had a satisfying accuracy and consistency in predicting children's IQs, which was confirmed by sensitivity analysis of gender and CpG site grouping data. In cell experiments, there was also a quantitative linear relationship between FAM50B DNAm and reactive oxygen species (ROS) production, which was mediated by PI3K-AKT signaling pathway. In addition, the lead BMD of ROS was close to that of FAM50B DNAm, suggesting that FAM50B DNAm was a suitable biomarker for the risk assessments of adverse outcomes induced by lead. Taken collectively, these results suggest that FAM50B/PTCHD3 can be applied to EHRA and the prevention/intervention of adverse effects of lead on children's IQs.

Novel DNA methylome biomarkers associated with adalimumab response in rheumatoid arthritis patients

Background and aims

Rheumatoid arthritis (RA) patients are currently treated with biological agents mostly aimed at cytokine blockade, such as tumor necrosis factor-alpha (TNFα). Currently, there are no biomarkers to predict therapy response to these agents. Here, we aimed to predict response to adalimumab (ADA) treatment in RA patients using DNA methylation in peripheral blood (PBL).

Methods

DNA methylation profiling on whole peripheral blood from 92 RA patients before the start of ADA treatment was determined using Illumina HumanMethylationEPIC BeadChip array. After 6 months, treatment response was assessed according to the European Alliance of Associations for Rheumatology (EULAR) criteria for disease activity. Patients were classified as responders (Disease Activity Score in 28 Joints (DAS28) < 3.2 or decrease of 1.2 points) or as non-responders (DAS28 > 5.1 or decrease of less than 0.6 points). Machine learning models were built through stability-selected gradient boosting to predict response prior to ADA treatment with predictor DNA methylation markers.

Results

Of the 94 RA patients, we classified 49 and 43 patients as responders and non-responders, respectively. We were capable of differentiating responders from non-responders with a high performance (area under the curve (AUC) 0.76) using a panel of 27 CpGs. These classifier CpGs are annotated to genes involved in immunological and pathophysiological pathways related to RA such as T-cell signaling, B-cell pathology, and angiogenesis.

Conclusion

Our findings indicate that the DNA methylome of PBL provides discriminative capabilities in discerning responders and non-responders to ADA treatment and may therefore serve as a tool for therapy prediction.

Methylation and transcriptomic expression profiles of HUVEC in the oxygen and glucose deprivation model and its clinical implications in AMI patients

The obstructed coronary artery undergoes a series of pathological changes due to ischemic-hypoxic shocks during acute myocardial infarction (AMI). However, the altered DNA methylation levels in endothelial cells under these conditions and their implication for the etiopathology of AMI have not been investigated in detail. This study aimed to explore the relationship between DNA methylation and pathologically altered gene expression profile in human umbilical vein endothelial cells (HUVECs) subjected to oxygen-glucose deprivation (OGD), and its clinical implications in AMI patients. The Illumina Infinium MethylationEPIC BeadChip assay was used to explore the genome-wide DNA methylation profile using the Novaseq6000 platform for mRNA sequencing in 3 pairs of HUVEC-OGD and control samples. GO and KEGG pathway enrichment analyses, as well as correlation, causal inference test (CIT), and protein-protein interaction (PPI) analyses identified 22 hub genes that were validated by MethylTarget sequencing as well as qRT-PCR. ELISA was used to detect four target molecules associated with the progression of AMI. A total of 2,524 differentially expressed genes (DEGs) and 22,148 differentially methylated positions (DMPs) corresponding to 6,642 differentially methylated genes (DMGs) were screened (|Δβ|>0.1 and detection p < 0.05). After GO, KEGG, correlation, CIT, and PPI analyses, 441 genes were filtered. qRT-PCR confirmed the overexpression of VEGFA, CCL2, TSP-1, SQSTM1, BCL2L11, and TIMP3 genes, and downregulation of MYC, CD44, BDNF, GNAQ, RUNX1, ETS1, NGFR, MME, SEMA6A, GNAI1, IFIT1, and MEIS1. DNA fragments BDNF_1_ (r = 0.931, p < 0.0001) and SQSTM1_2_NEW (r = 0.758, p = 0.0043) were positively correlated with the expressions of corresponding genes, and MYC_1_ (r = -0.8245, p = 0.001) was negatively correlated. Furthermore, ELISA confirmed TNFSF10 and BDNF were elevated in the peripheral blood of AMI patients (p = 0.0284 and p = 0.0142, respectively). Combined sequencing from in vitro cellular assays with clinical samples, aiming to establish the potential causal chain of the causal factor (DNA methylation) - mediator (mRNA)-cell outcome (endothelial cell ischemic-hypoxic injury)-clinical outcome (AMI), our study identified promising OGD-specific genes, which provided a solid basis for screening fundamental diagnostic and prognostic biomarkers of coronary endothelial cell injury of AMI. Moreover, it furnished the first evidence that during ischemia and hypoxia, the expression of BNDF was regulated by DNA methylation in endothelial cells and elevated in peripheral blood.

Demethylation of CDKN2A in systemic lupus erythematosus and rheumatoid arthritis: a blood biomarker for diagnosis and assessment of disease activity

INTRODUCTION/OBJECTIVES

Considering the phenotypic and serological heterogeneity of systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA), significant challenges may intervene with the precise diagnosis. In this regard, numerous studies have shown that changes in DNA methylation levels can be used to distinguish between healthy individuals and those with SLE and RA, as well as to predict disease activity and prognosis.

METHODS

In the current study, we evaluated quantitative methylation level of CDKN2A promoter in peripheral blood mononuclear cells (PBMCs) of SLE and RA patients, and healthy controls by methylation-quantification of endonuclease-resistant DNA (MethyQESD), a bisulfite conversion-independent method.

RESULTS

Our findings revealed an excessive hypomethylation of CDKN2A in SLE and RA patients compared to healthy individuals (P < 0.001). Besides, by determining efficient cutoff value, the specificity of CDKN2A for correct diagnosis of healthy subjects was measured to be 77.30% and the sensitivity for SLE and RA diagnosis were 81.33%, and 72%, respectively. Furthermore, CDKN2A methylation level was shown to be positively associated with C3 and C4 levels and negatively associated with anti‑dsDNA concentration (P < 0.001). Moreover, a statistically significant difference in the DNA methylation levels of CDKN2A promoter was identified between SLE cases with age of ≤ 18 and patients with > 18 years of age (P = 0.025).

CONCLUSION

Our findings demonstrated that CDKN2A methylation levels in PBMCs of SLE and RA patients could be used as a promising diagnostic biomarker. The significant association between hypomethylation of CDKN2A promoter and disease activity factors in SLE patients, is suggesting that CDKN2A hypomethylation could be used as an alternative biomarker for assessment of disease activity. Key Points • Several studies have reported increased expression of CDKN2A in SLE and RA suggesting that it may be involved in the pathogenesis of these disorders. • CDKN2A hypomethylation has been implicated in different autoimmune diseases. • Our findings demonstrated that CDKN2A methylation levels in PBMCs of SLE and RA patients could be used as a promising diagnostic and prognostic biomarker.

Differences of RUNX2 gene promoter methylation and transcription level in ankylosing spondylitis

BACKGROUND AND OBJECTIVE

Ankylosing spondylitis is a refractory immune disease that seriously affects the life and work of patients. Epigenetic modifications, especially DNA methylation, have become a research hotspot in complex diseases. We aim to explore the changes in runt-related transcription factor 2 (RUNX2) gene promoter methylation and transcription level in AS.

METHOD

We detected the RUNX2 gene promoter methylation in 83 AS patients and 83 healthy controls (HCs), then inspected the mRNA difference of RUNX2 between 30 AS patients and 30 HCs by the quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR).

RESULTS

The RUNX2 gene promoter was hypomethylated in AS patients compared to HCs (p < .001). The research involved 4 CpG regions and 74 CpG sites of RUNX2, of which CpG-2, CpG-4 regions, and 18 CpG sites have been differentially methylated. The CpG-4 island methylation was negatively correlated with C-reactive protein (p < .05) in AS patients. In the qRT-PCR validation phase, the mRNA level of RUNX2 in AS patients was significantly higher than HCs (p < .05), and in AS patients who were treated with biologics, the methylation level of CpG-2 island showed a negative correlation to mRNA (p < .05). ROC results indicated that RUNX2 methylation and its transcription level have good potential to distinguish AS patients from HCs.

CONCLUSION

The RUNX2 gene promoter was hypomethylated in AS patients. Meanwhile, the qRT-PCR verified the up-regulated expression on the transcription level, suggesting the abnormal methylation of RUNX2 contributes to the pathogenesis of AS.

A multiplex blood-based assay targeting DNA methylation in PBMCs enables early detection of breast cancer

The immune system can monitor tumor development, and DNA methylation is involved in the body's immune response to tumors. In this work, we investigate whether DNA methylation alterations in peripheral blood mononuclear cells (PBMCs) could be used as markers for early detection of breast cancer (BC) from the perspective of tumor immune alterations. We identify four BC-specific methylation markers by combining Infinium 850 K BeadChips, pyrosequencing and targeted bisulfite sequencing. Based on the four methylation markers in PBMCs of BC, we develop an efficient and convenient multiplex methylation-specific quantitative PCR assay for the detection of BC and validate its diagnostic performance in a multicenter cohort. This assay was able to distinguish early-stage BC patients from normal controls, with an AUC of 0.940, sensitivity of 93.2%, and specificity of 90.4%. More importantly, this assay outperformed existing clinical diagnostic methods, especially in the detection of early-stage and minimal tumors.

DNA Methylation Signatures of Response to Conventional Synthetic and Biologic Disease-Modifying Antirheumatic Drugs (DMARDs) in Rheumatoid Arthritis

Rheumatoid arthritis (RA) is a complex condition that displays heterogeneity in disease severity and response to standard treatments between patients. Failure rates for conventional, target synthetic, and biologic disease-modifying rheumatic drugs (DMARDs) are significant. Although there are models for predicting patient response, they have limited accuracy, require replication/validation, or for samples to be obtained through a synovial biopsy. Thus, currently, there are no prediction methods approved for routine clinical use. Previous research has shown that genetics and environmental factors alone cannot explain the differences in response between patients. Recent studies have demonstrated that deoxyribonucleic acid (DNA) methylation plays an important role in the pathogenesis and disease progression of RA. Importantly, specific DNA methylation profiles associated with response to conventional, target synthetic, and biologic DMARDs have been found in the blood of RA patients and could potentially function as predictive biomarkers. This review will summarize and evaluate the evidence for DNA methylation signatures in treatment response mainly in blood but also learn from the progress made in the diseased tissue in cancer in comparison to RA and autoimmune diseases. We will discuss the benefits and challenges of using DNA methylation signatures as predictive markers and the potential for future progress in this area.

The Association of Methylation Level in the CYP39A1 Gene with High Altitude Pulmonary Edema in the Chinese Population

Background

High altitude pulmonary edema (HAPE) is still the most common fatal disease at high altitudes. DNA methylation proceeds with an important role in HAPE progression. This study was designed to investigate the association between CYP39A1 methylation and HAPE.

Methods

Peripheral blood samples were enrolled from 106 participants (53 HAPE patients and 53 healthy subjects) to study the association of CYP39A1 methylation with HAPE. DNA methylation site in the promoter region of CYP39A1 was detected by Sequenom MassARRAY EpiTYPER platform.

Results

Probability analysis showed that the methylation probabilities of CYP39A1_1_CpG_5 and CYP39A1_3_CpG_21 are significant differences between the cases and controls (p< 0.05). The methylation level analysis indicated that CYP39A1_1_CpG_2.3.4, CYP39A1_5_CpG_6.7, and CYP39A1_5_CpG_9.10 were higher methylation in HAPE compared to the controls (p< 0.05). CYP39A1_3_CpG_21 and CYP39A1_4_CpG_3 exhibited a lower methylation level in HAPE than that in the controls (p< 0.05). The association analysis given that CYP39A1_1_CpG_2.3.4 (OR 2.56, p= 0.035), CYP39A1_5_CpG_6.7 (OR 3.99, p= 0.003), CYP39A1_5_CpG_9.10 (OR 3.99, p= 0.003), CYP39A1_5_CpG_16.17.18 (OR 2.53, p= 0.033), and CYP39A1_5_CpG_20 (OR 3.05, p= 0.031) are associated with an increased risk of HAPE. Whereas CYP39A1_1_CpG_5 (OR 0.33, p= 0.016) and CYP39A1_3_CpG_21 (OR 0.18, p= 0.005) have a protective role in HAPE. Besides, age-stratification analysis showed that CYP39A1_1_CpG_5 (OR 0.16, p= 0.014) and CYP39A1_3_CpG_21 (OR 0.08, p= 0.023) had a protective impact on HAPE in people aged ≤32 years. CYP39A1_5_CpG_6.7 (OR 6.70, p= 0.008) and CYP39A1_5_CpG_9.10 (OR 6.70, p= 0.008) were related to an increased susceptibility to HAPE aged >32 years. Moreover, the diagnostic value of CYP39A1_3_CpG_21 (AUC = 0.712, p< 0.001) was significantly better than other CpG sites.

Conclusion

The methylation level of CYP39A1 was associated with a risk of HAPE in the Chinese population, which provided new perspective for preventing and diagnosing of HAPE.

Current trends in epigenetic, cellular and molecular pathways in management of rheumatoid arthritis

Rheumatoid arthritis is a systemic chronic polyarticular autoimmune disorder of joints and joint membrane mainly affecting feet and hands. The pathological manifestation of the disease includes infiltration of immune cells, hyperplasia of the lining of synovium, formation of pannus and bone and cartilage destruction. If left untreated, the appearance of small focal necrosis, adhesion of granulation, and formation of fibrous tissue on the surface of articular cartilage is noted. The disease primarily affects nearly 1% of the population globally, women being more affected than men with a ratio 2:1 and can initiate regardless of any age. The synovial fibroblast in rheumatoid arthritis individuals exhibits an aggressive phenotype which upregulates the manifestation of protooncogenes, adhesive compounds, inflammatory cytokines and matrix-deteriorating enzymes. Apart from the inflammatory effects of cytokines, chemokines are also noted to induce swelling and pain in arthritic individuals by residing in synovial membrane and forming pannus. The current treatment of rheumatoid arthritis includes treatment with non-steroidal anti-inflammatory drugs, disease-modifying antirheumatic drugs, treatment with biologics such as inhibitors of TNF-α, interleukins, platelet activating factor, etc. which provides significant relief from symptoms and aids in management of the disease. The current review highlights the pathogenesis involved in the onset of rheumatoid arthritis and also covers epigenetic, cellular and molecular parameters associated with it to aid better and advanced therapeutic approaches for management of the debilitating disease.

Differential CpG DNA methylation of peripheral B cells, CD4+ T cells, and salivary gland tissues in IgG4-related disease

OBJECTIVES

Immunoglobulin-G4-related disease (IgG4-RD) is a distinct systemic autoimmune-mediated disease manifesting as chronic inflammation and tissue fibrosis. Since the role of DNA methylation in the pathogenesis of IgG4-RD is still unclear, we conduct this study to investigate epigenetic modifications in IgG4-RD.

METHODS

A genome-wide DNA methylation study was conducted with B cells, CD4⁺ T cells, and salivary gland tissues from IgG4-RD patients and matched controls by using the Illumina HumanMethylation 850K BeadChip. We further performed pyrosequencing and immunohistochemistry assays to validate the methylation status of some targets of interest.

RESULTS

We identified differentially methylated CpG sites including 44 hypomethylated and 166 hypermethylated differentially methylated probes (DMPs) in B cells and 260 hypomethylated and 112 hypermethylated DMPs in CD4⁺ T cells from 10 IgG4-RD patients compared with 10 healthy controls. We also identified 36945 hypomethylated and 78380 hypermethylated DMPs in salivary gland tissues of 4 IgG4-RD patients compared with 4 controls. DPM2 (cg21181453), IQCK (cg10266221), and ABCC13 (cg05699681, cg04985582) were hypermethylated and MBP (cg18455083) was hypomethylated in B cells, CD4⁺ T cells, and salivary gland tissues of IgG4-RD patients. We also observed the hypomethylated HLA-DQB2 in CD4⁺ T cells from IgG4-RD patients. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of DMPs in salivary gland tissues of IgG4-RD patients revealed enrichment of pathways involved in the regulation of immune cell responses and fibrosis.

CONCLUSION

This is the first DNA methylation study in peripheral B cells, CD4⁺ T cells, and salivary gland tissues from IgG4-RD patients. Our findings highlighted the role of epigenetic modification of DNA methylation and identified several genes and pathways possibly involved in IgG4-RD pathogenesis.

Seropositivity-Dependent Association between LINE-1 Methylation and Response to Methotrexate Therapy in Early Rheumatoid Arthritis Patients

Background: Methotrexate (MTX) is considered the first choice among disease-modifying anti-rheumatic drugs (DMARDs) for rheumatoid arthritis (RA) treatment. However, response to it varies as approximately 40% of the patients do not respond and would lose the most effective period of treatment time. Therefore, having a predictive biomarker before starting MTX treatment is of utmost importance. Methylation of long interspersed nucleotide element-1 (LINE-1) is generally considered a surrogate marker for global genomic methylation, which has been reported to associate with disease activity after MTX therapy. Methods: We performed a prospective study on 273 naïve early RA (ERA) patients who were treated with MTX, followed up to 12 months, and classified according to their therapy response. The baseline LINE-1 methylation levels in peripheral blood mononuclear cells (PBMC) of cases were assessed by bisulfite pyrosequencing. Results: Baseline LINE-1 methylation level per se turned out not to predict the response to the therapy, nor did age, sex, body mass index, or smoking status. However, if cases were stratified according to positivity to rheumatoid factor (RF) and anti-citrullinated protein antibody (ACPA) or seronegativity, we observed an opposite association between baseline LINE-1 methylation levels and optimal response to MTX therapy among responders. The best response to MTX therapy was associated with hypermethylated LINE-1 among double-positive ERA cases (p-value: 0.002) and with hypomethylated LINE-1 in seronegative ERA patients (p-value: 0.01). Conclusion: The LINE-1 methylation level in PBMCs of naïve ERA cases associates with the degree of response to MTX therapy in an opposite way depending on the presence of RF and ACPA antibodies. Our results suggest LINE-1 methylation level as a new epigenetic biomarker for predicting the degree of response to MTX in both double-positive and seronegative ERA patients.

Single-cell RNA sequencing analysis of T helper cell differentiation and heterogeneity

Single-cell transcriptomics has emerged as a powerful tool to investigate cells' biological landscape and focus on the expression profile of individual cells. Major advantage of this approach is an analysis of highly complex and heterogeneous cell populations, such as a specific subpopulation of T helper cells that are known to differentiate into distinct subpopulations. The need for distinguishing the specific expression profile is even more important considering the T cell plasticity. However, importantly, the universal pipelines for single-cell analysis are usually not sufficient for every cell type. Here, the aims are to analyze the diversity of T cell phenotypes employing classical in vitro cytokine-mediated differentiation of human T cells isolated from human peripheral blood by single-cell transcriptomic approach with support of labelled antibodies and a comprehensive bioinformatics analysis using combination of Seurat, Nebulosa, GGplot and others. The results showed high expression similarities between Th1 and Th17 phenotype and very distinct Th2 expression profile. In a case of Th2 highly specific marker genes SPINT2, TRIB3 and CST7 were expressed. Overall, our results demonstrate how donor difference, Th plasticity and cell cycle influence the expression profiles of distinct T cell populations. The results could help to better understand the importance of each step of the analysis when working with T cell single-cell data and observe the results in a more practical way by using our analyzed datasets.

Identification of methylation signatures and rules for predicting the severity of SARS-CoV-2 infection with machine learning methods

Patients infected with SARS-CoV-2 at various severities have different clinical manifestations and treatments. Mild or moderate patients usually recover with conventional medical treatment, but severe patients require prompt professional treatment. Thus, stratifying infected patients for targeted treatment is meaningful. A computational workflow was designed in this study to identify key blood methylation features and rules that can distinguish the severity of SARS-CoV-2 infection. First, the methylation features in the expression profile were deeply analyzed by a Monte Carlo feature selection method. A feature list was generated. Next, this ranked feature list was fed into the incremental feature selection method to determine the optimal features for different classification algorithms, thereby further building optimal classifiers. These selected key features were analyzed by functional enrichment to detect their biofunctional information. Furthermore, a set of rules were set up by a white-box algorithm, decision tree, to uncover different methylation patterns on various severity of SARS-CoV-2 infection. Some genes (PARP9, MX1, IRF7), corresponding to essential methylation sites, and rules were validated by published academic literature. Overall, this study contributes to revealing potential expression features and provides a reference for patient stratification. The physicians can prioritize and allocate health and medical resources for COVID-19 patients based on their predicted severe clinical outcomes.

Disease mechanisms in preclinical rheumatoid arthritis: A narrative review

In the last decades, the concept of preclinical rheumatoid arthritis (RA) has become established. In fact, the discovery that disease mechanisms start years before the onset of clinical RA has been one of the major recent insights in the understanding of RA pathogenesis. In accordance with the complex nature of the disease, preclinical events extend over several sequential phases. In a genetically predisposed host, environmental factors will further increase susceptibility for incident RA. In the initial steps of preclinical disease, immune disturbance mechanisms take place outside the joint compartment, namely in mucosal surfaces, such as the lung, gums or gut. Herein, the persistent immunologic response to altered antigens will lead to breach of tolerance and trigger autoimmunity. In a second phase, the immune response matures and is amplified at a systemic level, with epitope spreading and widening of the autoantibody repertoire. Finally, the synovial and bone compartment are targeted by specific autoantibodies against modified antigens, initiating a local inflammatory response that will eventually culminate in clinically evident synovitis. In this review, we discuss the elaborate disease mechanisms in place during preclinical RA, providing a broad perspective in the light of current evidence.

Interferon-α-mediated therapeutic resistance in early rheumatoid arthritis implicates epigenetic reprogramming

OBJECTIVES

An interferon (IFN) gene signature (IGS) is present in approximately 50% of early, treatment naive rheumatoid arthritis (eRA) patients where it has been shown to negatively impact initial response to treatment. We wished to validate this effect and explore potential mechanisms of action.

METHODS

In a multicentre inception cohort of eRA patients (n=191), we examined the whole blood IGS (MxA, IFI44L, OAS1, IFI6, ISG15) with reference to circulating IFN proteins, clinical outcomes and epigenetic influences on circulating CD19+ B and CD4+ T lymphocytes.

RESULTS

We reproduced our previous findings demonstrating a raised baseline IGS. We additionally showed, for the first time, that the IGS in eRA reflects circulating IFN-α protein. Paired longitudinal analysis demonstrated a significant reduction between baseline and 6-month IGS and IFN-α levels (p<0.0001 for both). Despite this fall, a raised baseline IGS predicted worse 6-month clinical outcomes such as increased disease activity score (DAS-28, p=0.025) and lower likelihood of a good EULAR clinical response (p=0.034), which was independent of other conventional predictors of disease activity and clinical response. Molecular analysis of CD4+ T cells and CD19+ B cells demonstrated differentially methylated CPG sites and dysregulated expression of disease relevant genes, including PARP9, STAT1, and EPSTI1, associated with baseline IGS/IFNα levels. Differentially methylated CPG sites implicated altered transcription factor binding in B cells (GATA3, ETSI, NFATC2, EZH2) and T cells (p300, HIF1α).

CONCLUSIONS

Our data suggest that, in eRA, IFN-α can cause a sustained, epigenetically mediated, pathogenic increase in lymphocyte activation and proliferation, and that the IGS is, therefore, a robust prognostic biomarker. Its persistent harmful effects provide a rationale for the initial therapeutic targeting of IFN-α in selected patients with eRA.

Geniposide restricts angiogenesis in experimentary arthritis via inhibiting Dnmt1-mediated PTEN hypermethylation

Neovascularization in rheumatoid arthritis (RA) is a key bridge between malignant proliferative synovial tissue and pannus. In view of previous studies on the efficacy of Geniposide (GE) in experimentary arthritis, the purpose of this study was to investigate the possible mechanism of GE inhibiting angiogenesis by regulating the gene of phosphate and tension homology deleted on chromosome ten (PTEN). In this study, human umbilical vein endothelial cells (HUVEC) and adjuvant arthritis (AA) rat models were performed to research in vitro and in vivo. The results showed that GE treatment significantly reduced synovitis and angiogenesis in AA rats, which may be associated with the increased expression of PTEN with GE treatment. Meanwhile, the hypermethylation of PTEN accompanied by the over-expression of DNA methyltransferases (Dnmts) was demonstrated in TNF-α-induced HUVEC and AA rats. Knockdown of Dnmt1 by Dnmt1- siRNA significantly inhibited the tube formation of HUVEC in vitro. GE significantly restricted the angiogenesis of HUVEC by inhibiting DNA methylation, which was attributed to the down-regulation of Dnmt1 rather than Dnmt3a and Dnmt3b. The anti-angiogenesis effect of GE was further verified in AA model by the inhibition of Dnmt1. These results indicate that GE exhibited anti-angiogenesis effects in experimentary arthritis by inhibiting Dnmt1-mediated PTEN gene hypermethylation, which may brings new insights for the prevention and research of RA.

DNA Methylation and mRNA Expression of B7-H3 Gene in Ankylosing Spondylitis: A Case-Control Study

BACKGROUND

Ankylosing spondylitis (AS) is a common inflammatory arthritis, with a high prevalence in patients in their mid-20s. Its pathogenesis is not well understood; however, genetic factors likely play a critical role. Epigenetic DNA changes may be involved in the pathogenesis of AS. In this study, we explored the methylation and transcription levels of the B7-H3 gene and its association with AS in an eastern Chinese Han population.

METHODS

Peripheral blood of AS patients and healthy controls was used to extract genomic DNA and B7-H3 methylation levels were analyzed using sodium bisulfite followed by multiplex polymerase chain reaction. SPSS software was used to determine the statistical significance of the results.

RESULTS

Hypomethylation of the promoter of the B7-H3 gene was observed in AS patients, whereas the B7-H3 gene expression was significantly enhanced in AS patients.

CONCLUSION

Epigenetic modifications of B7-H3 were associated with susceptibility to AS. Hypomethylation of the B7-H3 promoter, which leads to B7-H3 overexpression, may be involved in the pathogenesis of AS.

Mitochondrial Dysfunction and Oxidative Stress in Rheumatoid Arthritis

Control of excessive mitochondrial oxidative stress could provide new targets for both preventive and therapeutic interventions in the treatment of chronic inflammation or any pathology that develops under an inflammatory scenario, such as rheumatoid arthritis (RA). Increasing evidence has demonstrated the role of mitochondrial alterations in autoimmune diseases mainly due to the interplay between metabolism and innate immunity, but also in the modulation of inflammatory response of resident cells, such as synoviocytes. Thus, mitochondrial dysfunction derived from several danger signals could activate tricarboxylic acid (TCA) disruption, thereby favoring a vicious cycle of oxidative/mitochondrial stress. Mitochondrial dysfunction can act through modulating innate immunity via redox-sensitive inflammatory pathways or direct activation of the inflammasome. Besides, mitochondria also have a central role in regulating cell death, which is deeply altered in RA. Additionally, multiple evidence suggests that pathological processes in RA can be shaped by epigenetic mechanisms and that in turn, mitochondria are involved in epigenetic regulation. Finally, we will discuss about the involvement of some dietary components in the onset and progression of RA.

DNA Methylation of the Natriuretic Peptide System Genes and Ischemic Stroke: Gene-Based and Gene Set Analyses

Background and Objectives

The natriuretic peptide (NP) system has been considered an important regulator for ischemic stroke (IS) with a limited clinical implication. A better understanding of the underlying molecular mechanisms is urgent. Here, we aimed to examine the role of DNA methylation of NP system genes in IS.

Methods

DNA methylation at promoter regions of 4 core NP system genes, e.g., CORIN, FURIN, NPPA, and NPPB, was measured by targeted bisulfite sequencing in 853 patients with IS and 918 controls. We first examined the association between DNA methylation at each single CpG and IS, followed by gene-based and gene set analyses to examine the joint associations of DNA methylation at multiple CpGs in a gene or all 4 genes as a pathway with IS.

Results

After control of covariates and multiple testing, DNA methylation at 19 of the 36 assayed CpGs was individually associated with IS at q < 0.05. Higher average methylation levels at the targeted regions of CORIN (odds ratio [OR] = 0.64, 95% confidence interval [CI]: 0.56–0.73), FURIN (OR = 0.78, 95% CI: 0.69–0.88), and NPPA (OR = 0.78, 95% CI: 0.69–0.88) were associated with a lower odds of IS (all q < 0.05). The truncated product method revealed the same gene-based associations (all q < 0.05) and found that DNA methylation at all 4 NP system genes together was jointly associated with IS (p = 0.0001).

Discussion

DNA methylation at NP system genes was downregulated in patients with IS. Our results may unravel a molecular mechanism underlying the regulating effect of the NP system on IS and highlight the relevance of testing the joint effect of multiple CpGs in the epigenetic analysis.

Gene Interaction Network Analysis Reveals IFI44L as a Drug Target in Rheumatoid Arthritis and Periodontitis

Simple Summary

In spite of substantial investigation, the biological link between periodontitis and rheumatoid arthritis remains unexplained. This study intended to correlate periodontitis and rheumatoid arthritis gene expression patterns to find shared targets for both the disease. We identified the differentially expressed genes (DEGs) in periodontitis and rheumatoid arthritis. The network was built by integrating DEGs and ranking the genes using GeneMANIA. FINDSITEcomb2.0 was used to find a possible inhibitor for the top-ranked gene. Further, the binding effectiveness and protein-ligand complex stability were then determined by molecular docking and molecular dynamics. The network analysis showed IFI44L as a highly ranking gene implicated in most immunological pathways. A virtual screening of 6507 compounds revealed vemurafenib as the best candidate for the IFI44L target. Molecular docking and molecular dynamics modelling revealed the stability of the IFI44L-vemurafenib complex, which suggest IFI44L is potential target and vemurafenib could be the better candidate to treat both diseases.

Abstract

Objective: Despite extensive research on periodontitis and rheumatoid arthritis, the underlying molecular connectivity between these condition remains largely unknown. This research aimed to integrate periodontitis and rheumatoid arthritis gene expression profiles to identify interconnecting genes and focus to develop a common lead molecule against these inflammatory conditions. Materials and Methods: Differentially expressed genes (DEGs) of periodontitis and rheumatoid arthritis were identified from the datasets retrieved from the Gene Expression Omnibus database. The network was constructed by merging DEGs, and the interconnecting genes were identified and ranked using GeneMANIA. For the selected top ranked gene, the potential inhibitor was searched using FINDSITE^comb2.0. Subsequently, the molecular docking and molecular dynamics were performed to determine the binding efficiency and protein-ligand complex stability, respectively. Results: From the network analysis, IFN-induced protein 44-like (IFI44L) was identified as a top ranked gene involved in most of the immunological pathway. With further virtual screening of 6507 molecules, vemurafenib was identified to be the best fit against the IFI44L target. The binding energy and stability of IFI44L with vemurafenib were investigated using molecular docking and molecular dynamics simulation. Docking results show binding energy of −7.7 Kcal/mol, and the simulation results show stability till 100 ns. Conclusions: The identified IFI44L may represent a common drug target for periodontitis and rheumatoid arthritis. Vemurafenib could be a potent anti-inflammatory drug for both diseases.

Epigenetic Regulation of Immune and Inflammatory Responses in Rheumatoid Arthritis

Purpose

Rheumatoid arthritis (RA) is a disease associated with multiple factors. Epigenetics can affect gene expression without altering the DNA sequence. In this study, we aimed to comprehensively analyze epigenetic regulation in RA.

Methods

Using the Gene Expression Omnibus database, we identified a methylation chip, RNA-sequencing, and miRNA microarray for RA. First, we searched for DNA methylation, genes, and miRNAs associated with RA using differential analysis. Second, we determined the regulatory networks for RA-specific methylation, miRNA, and m6A using cross-analysis. Based on these three regulatory networks, we built a comprehensive epigenetic regulatory network and identified hub genes.

Results

Using a differential analysis, we identified 16,852 differentially methylated sites, 4877 differentially expressed genes, and 32 differentially expressed miRNAs. The methylation-expression regulatory network was mainly associated with the PI3K-Akt and T-cell receptor signaling pathways. The miRNA expression regulatory network was mainly related to the MAPK and chemokine signaling pathways. M6A regulatory network was mainly associated with the MAPK signaling pathway. Additionally, five hub genes were identified in the epigenetic regulatory network: CHD3, SETD1B, FBXL19, SMARCA4, and SETD1A. Functional analysis revealed that these five genes were associated with immune cells and inflammatory responses.

Conclusion

We constructed a comprehensive epigenetic network associated with RA and identified core regulatory genes. This study provides a new direction for future research on the epigenetic mechanisms of RA.

piRNA-6426 increases DNMT3B-mediated SOAT1 methylation and improves heart failure

PURPOSE

Previous studies found that piRNAs could participate in disease progression by regulating DNA methylation, but there are few reports on their roles in heart failure (HF).

METHODS

The level of piRNA-6426 in the venous blood of HF patients and volunteers was detected by RT-qPCR. Hypoxia-induced cardiomyocytes were transfected with lentiviral-mediated piRNA-6426 overexpression vector (LV-piRNA-6426) or together with LV-DNMT3B, and then cell viability and apoptosis, glucose uptake, ROS production, LDH activity and secretion of inflammatory factors were detected. Also, cardiomyocytes were transfected with LV-piRNA-6426, sh-piRNA-6426 or sh-SOAT1, as well as LV-piRNA-6426 or together with LV-DNMT3B or sh-DNMT3B. The interaction between piRNA-6426 and methyltransferase 3B (DNMT3B) was detected with RNA immunoprecipitation (RIP). And the methylation level of sterol o-acyltransferase 1 (SOAT1) and the enrichment of DNMT3B in the SOAT1 promoter were detected with Methylation-specific PCR (MSP) and ChIP assays. Then a HF rat model constructed with coronary artery occlusion method was injected with LV-piRNA-6426, and heart function index and infarcted area of rat heart were detected.

RESULTS

piRNA-6426 expression was decreased in the blood of HF patients. LV-piRNA-6426 transfection increased the enrichment of DNMT3B in SOAT1 promoter, thereby inhibiting the expression level of SOAT1, and decreased hypoxia-induced oxidative stress and inflammation in cardiomyocytes, while sh-piRNA-6426 transfection had the opposite effect. And LV-DNMT3B transfection enhanced the effect of LV-piRNA-6426 transfection on SOAT1 expression and cardiomyocyte dysfunction. Injection of LV-piRNA-6426 significantly inhibited the heart dysfunction of rats.

CONCLUSIONS

piRNA-6426 overexpression inhibits hypoxia-induced cardiomyocyte dysfunction and HF by promoting DNMT3B-mediated methylation of SOAT1 promoter.

Epigenetic Regulation in the Pathogenesis of Rheumatoid Arthritis

Rheumatoid arthritis (RA) is an autoimmune disease. The etiology of RA remains undetermined and the pathogenesis is complex. There remains a paucity of ideal therapeutic drugs and treatment strategies. The epigenetic modifications affect and regulate the function and characteristics of genes through mechanisms, including DNA methylation, histone modification, chromosome remodeling, and RNAi, thereby exerting a significant impact on the living state of the body. Recently, the phenomenon of epigenetic modification in RA has garnered growing research interest. The application of epigenetically modified methods is the frontier field in the research of RA pathogenesis. This review highlights the research on the pathogenesis of RA based on epigenetic modification in the recent five years, thereby suggesting new methods and strategies for the diagnosis and treatment of RA.

Association of DNA Methylation in Blood Pressure-Related Genes With Ischemic Stroke Risk and Prognosis

Background

A genome-wide association study identified 12 genetic loci influencing blood pressure and implicated a role of DNA methylation. However, the relationship between methylation and ischemic stroke has not yet been clarified. We conducted a large-sample sequencing study to identify blood leukocyte DNA methylations as novel biomarkers for ischemic stroke risk and prognosis based on previously identified genetic loci.

Methods

Methylation levels of 17 genes were measured by sequencing in 271 ischemic stroke cases and 323 controls, and the significant associations were validated in another independent sample of 852 cases and 925 controls. The associations between methylation levels and ischemic stroke risk and prognosis were evaluated.

Results

Methylation of AMH, C17orf82, HDAC9, IGFBP3, LRRC10B, PDE3A, PRDM6, SYT7 and TBX2 was significantly associated with ischemic stroke. Compared to participants without any hypomethylated targets, the odds ratio (OR) (95% confidence interval, CI) for those with 9 hypomethylated genes was 1.41 (1.33–1.51) for ischemic stroke. Adding methylation levels of the 9 genes to the basic model of traditional risk factors significantly improved the risk stratification for ischemic stroke. Associations between AMH, HDAC9, IGFBP3, PDE3A and PRDM6 gene methylation and modified Rankin Scale scores were significant after adjustment for covariates. Lower methylation levels of AMH, C17orf82, PRDM6 and TBX2 were significantly associated with increased 3-month mortality. Compared to patients without any hypomethylated targets, the OR (95% CI) for those with 4 hypomethylated targets was 1.12 (1.08–1.15) for 3-month mortality (P = 2.28 × 10−10).

Conclusion

The present study identified blood leukocyte DNA methylations as potential factors affecting ischemic stroke risk and prognosis among Han Chinese individuals.

Identification of novel rheumatoid arthritis-associated MiRNA-204-5p from plasma exosomes

Rheumatoid arthritis (RA) is an autoimmune disease characterized by infiltration of immune cells in the synovium. However, the crosstalk of immune cells and synovial fibroblasts is still largely unknown. Here, global miRNA screening in plasma exosomes was carried out with a custom microarray (RA patients vs. healthy controls = 9:9). A total of 14 exosomal miRNAs were abnormally expressed in the RA patients. Then, downregulated expression of exosomal miR-204-5p was confirmed in both the replication (RA patients vs. healthy controls = 30:30) and validation groups (RA patients vs. healthy controls = 56:60). Similar to the findings obtained in humans, a decreased abundance of exosomal miR-204-5p was observed in mice with collagen-induced arthritis (CIA). Furthermore, Spearman correlation analysis indicated that plasma exosomal miR-204-5p expression was inversely correlated with disease parameters of RA patients, such as rheumatoid factor, erythrocyte sedimentation rate, and C-reactive protein. In vitro, our data showed that human T lymphocytes released exosomes containing large amounts of miR-204-5p, which can be transferred into synovial fibroblasts, inhibiting cell proliferation. Overexpression of miR-204-5p in synovial fibroblasts suppressed synovial fibroblast activation by targeting genes related to cell proliferation and invasion. In vivo assays found that administration of lentiviruses expressing miR-204-5p markedly alleviated the disease progression of the mice with CIA. Collectively, this study identified a novel RA-associated plasma exosomal miRNA-204-5p that mediates the communication between immune cells and synovial fibroblasts and can be used as a potential biomarker for RA diagnosis and treatment.

A microRNA that is significantly reduced in joint tissues in rheumatoid arthritis could provide a therapeutic target and act as a biomarker for disease progression. In rheumatoid arthritis, immune cells release exosomes, tiny vesicles containing microRNA and proteins that are transferred to cells in the synovium, the connective tissue lining the inside of the joint capsule. This transfer of molecules influences synovial cell activity. Shu-Feng Lei and Fei-Yan Deng at the Medical School of Soochow University, Suzhou, China, and co-workers identifed exosomal microRNAs present in rheumatoid arthritis, and examined their effect on synovial cells. Levels of one exosomal microRNA, miR-204-5p, were significantly lower in patient samples and mice models, inversely correlating with disease severity. The team believe that chronic inflammation may suppress levels of miR-204-5p. Treatment boosting microRNA levels in mice models slowed disease progression.

The Role of Oxidative Stress in Epigenetic Changes Underlying Autoimmunity

Significance: Epigenetic dysregulation plays an important role in the pathogenesis and development of autoimmune diseases. Oxidative stress is associated with autoimmunity and is also known to alter epigenetic mechanisms. Understanding the interplay between oxidative stress and epigenetics will provide insights into the role of environmental triggers in the development of autoimmunity in genetically susceptible individuals. Recent Advances: Abnormal DNA and histone methylation patterns in genes and pathways involved in interferon and tumor necrosis factor signaling, cellular survival, proliferation, metabolism, organ development, and autoantibody production have been described in autoimmunity. Inhibitors of DNA and histone methyltransferases showed potential therapeutic effects in animal models of autoimmune diseases. Oxidative stress can regulate epigenetic mechanisms via effects on DNA damage repair mechanisms, cellular metabolism and the local redox environment, and redox-sensitive transcription factors and pathways. Critical Issues: Studies looking into oxidative stress and epigenetics in autoimmunity are relatively limited. The number of available longitudinal studies to explore the role of DNA methylation in the development of autoimmune diseases is small. Future Directions: Exploring the relationship between oxidative stress and epigenetics in autoimmunity will provide clues for potential preventative measures and treatment strategies. Inception cohorts with longitudinal follow-up would help to evaluate epigenetic marks as potential biomarkers for disease development, progression, and treatment response in autoimmunity. Antioxid. Redox Signal. 36, 423-440.

Promoter DNA Methylation in GWAS-Identified Genes as Potential Functional Elements for Blood Pressure: An Observational and Mendelian Randomization Study

Genome-wide association studies have identified numerous genetic loci for blood pressure (BP). However, the relationships of functional elements inside these loci with BP are not fully understood. This study represented an effort to determine if promoter DNA methylations inside BP-associated loci were associated with BP.We conducted a cross-sectional study investigating the association between promoter DNA methylations of 10 candidate genes and BP in 1,241 Chinese individuals. Twenty-one genomic fragments in the CpG Islands were sequenced. The associations of methylation levels with BP and hypertension were assessed in regression models. Mendelian randomization (MR) analysis was then applied to find supporting evidence for the identified associations.A total of 413 DNA methylation sites were examined in an observational study. Methylation levels of 24 sites in PRDM6, IGFBP3, SYT7, PDE3A, TBX2 and C17orf82 were significantly associated with BP. Methylation levels of PRDM6 and SYT7 were significantly associated with hypertension. Methylation levels of five sites (including cg06713098) in IGFBP3 were significantly associated with DBP. MR analysis found associations between the methylation levels of six CpG sites (cg06713098, cg14228300, cg23193639, cg21268650, cg10677697 and cg04812164) around the IGFBP3 promoter and DBP. Methylation levels of cg14228300 and cg04812164 were associated with SBP. By further applying several MR methods we showed that the associations may not be due to pleiotropy. Association between IGFBP3 mRNA levels in blood cells and BP was also found in MR analysis. This study identified promoter methylation as potential functional element for BP. The identified methylations may be involved in the regulatory pathway linking genetic variants to BP.

NPPA Promoter Hypomethylation Predicts Central Obesity Development: A Prospective Longitudinal Study in Chinese Adults

INTRODUCTION

Atrial natriuretic peptide plays a potential role in obesity with unclear molecular mechanisms. The objective of this study was to examine the association between its coding gene (natriuretic peptide A [NPPA]) methylation and obesity.

METHODS

Peripheral blood DNA methylation of NPPA promoter was quantified at baseline by targeted bisulfite sequencing for 2,497 community members (mean aged 53 years, 38% men) in the Gusu cohort. Obesity was repeatedly assessed by body mass index (BMI) and waist circumference (WC) at baseline and follow-up examinations. The cross-sectional, longitudinal, and prospective associations between NPPA promoter methylation and obesity were examined.

RESULTS

Of the 9 CpG loci assayed, DNA methylation levels at 6 CpGs were significantly lower in participants with central obesity than those without (all p < 0.05 for permutation test). These CpG methylation levels at baseline were also inversely associated with dynamic changes in BMI or WC during follow-up (all p < 0.05 for permutation test). After an average 4 years of follow-up, hypermethylation at the 6 CpGs (CpG2 located at Chr1:11908348, CpG3 located at Chr1:11908299, CpG4 located at Chr1:11908200, CpG5 located at Chr1:11908182, CpG6 located at Chr1:11908178, and CpG8 located at Chr1:11908165) was significantly associated with a lower risk of incident central obesity (all p < 0.05 for permutation test).

CONCLUSIONS

Hypomethylation at NPPA promoter was associated with increased future risk of central obesity in Chinese adults. Aberrant DNA methylation of the NPPA gene may participate in the mechanisms of central obesity.

[Methylation status and expression of TWEAK gene promoter region in peripheral blood of patients with rheumatoid arthritis]

OBJECTIVE

To explore the relationship between tumor necrosis factor like weak inducer of apoptosis (TWEAK) gene and the pathogenesis of rheumatoid arthritis (RA) by detecting the DNA methylation level, mRNA expression level and serum protein concentration of TWEAK gene in peripheral blood.

METHODS

The MassARRAY method was used to detect the DNA methylation level of the TWEAK gene in the peripheral blood of 112 RA patients and 86 matched healthy volunteers. The real-time quantitative polymerase chain reaction method was used to detect the mRNA expression level of the TWEAK gene in the peripheral blood of the subjects. The enzyme-linked immunosorbent assay method was used to detect the serum TWEAK protein concentration of the subjects. The TWEAK gene DNA methylation level, mRNA expression level and serum protein concentration between the RA group and the healthy control group were compared, and the relationship between it and the degree of disease activity analyzed.

RESULTS

The overall DNA methylation level of TWEAK gene and the DNA methylation levels of CpG_11, CpG_17.18.19.20, CpG_40.41.42 site in the RA group were higher than those in the healthy control group (P=0.002, P=0.01, P=0.006, P=0.002, respectively). The DNA methylation level of CpG_55.56 site in the high disease activity group was higher than that in the medium and low disease activity group (P=0.041). The expression level of TWEAK gene mRNA in the peripheral blood of the RA group was lower than that of the healthy control group (P=0.023). The expression level of TWEAK gene mRNA in the high disease activity group was lower than that in the medium and low disease activity group (P=0.035). The serum TWEAK protein concentration of the RA group was not significantly different from that of the healthy control group (P=0.508), but it was positively correlated with the mRNA expression level (r=0.482, P < 0.001).

CONCLUSION

The TWEAK gene is closely related to the onset and progression of RA, and its hypermethylation state may be one of the epigenetic mechanisms regulating its low mRNA expression, and it can be used as one of the important indicators for clinical monitoring and evaluation of RA.

Etiology and Risk Factors for Rheumatoid Arthritis: A State-of-the-Art Review

Rheumatoid arthritis (RA) is the most common systemic inflammatory rheumatic disease. It is associated with significant burden at the patient and societal level. Extensive efforts have been devoted to identifying a potential cause for the development of RA. Epidemiological studies have thoroughly investigated the association of several factors with the risk and course of RA. Although a precise etiology remains elusive, the current understanding is that RA is a multifactorial disease, wherein complex interactions between host and environmental factors determine the overall risk of disease susceptibility, persistence and severity. Risk factors related to the host that have been associated with RA development may be divided into genetic; epigenetic; hormonal, reproductive and neuroendocrine; and comorbid host factors. In turn, environmental risk factors include smoking and other airborne exposures; microbiota and infectious agents; diet; and socioeconomic factors. In the present narrative review, aimed at clinicians and researchers in the field of RA, we provide a state-of-the-art overview of the current knowledge on this topic, focusing on recent progresses that have improved our comprehension of disease risk and development.

Fetal malnutrition is associated with impairment of endogenous melatonin synthesis in pineal via hypermethylation of promoters of protein kinase C alpha and cAMP response element-binding

This study investigated whether and how fetal malnutrition would influence endogenous melatonin synthesis, and whether such effect of fetal malnutrition would transmit to the next generation. We enrolled 2466 participants and 1313 of their offspring. The urine 6-hydroxymelatonin sulfate and serum melatonin rhythm were measured. Methylation microarray detection and bioinformatics analysis were performed to identify hub methylated sites. Additionally, rat experiment was performed to elucidate mechanisms. The participants with fetal malnutrition had lower 6-hydroxymelatonin sulfate (16.59 ± 10.12 μg/24 hours vs 24.29 ± 11.99 μg/24 hours, P < .001) and arear under curve of melatonin rhythm (67.11 ± 8.16 pg/mL vs 77.11 ± 8.04 pg/mL, P < .001). We identified 961 differentially methylated sites, in which the hub methylated sites were locating on protein kinase C alpha (PRKCA) and cAMP response element-binding protein (CREB1) promoters, mediating the association of fetal malnutrition with impaired melatonin secretion. However, such effects were not observed in the offspring (all P > .05). Impaired histomorphology of pineal, decreased melatonin in serum, pineal, and pinealocyte were also found in the in vivo and in vitro experiments (P < .05 for the differences of the indicators). Hypermethylation of 10 CpG sites on the PRKCA promoter and 8 CpG sites on the CREB1 promoter were identified (all P < .05), which down-regulated PRKCA and CREB1 expressions, leading to decreased expression of AANAT, and then resulting in the impaired melatonin synthesis. Collectively, fetal malnutrition can impair melatonin synthesis through hypermethylation of PRKCA and CREB1 promoters, and such effects cannot be transmitted to the next generation.

Ankylosing Spondylitis Patients Display Aberrant ERAP1 Gene DNA Methylation and Expression

BACKGROUND

Endoplasmic reticulum aminopeptidase 1 (ERAP1) is known to participate in the pathogenesis of ankylosing spondylitis (AS). This study aimed to evaluate the relationship between promoter methylation and mRNA levels of ERAP1 and AS susceptibility.

METHODS

DNA methylation levels of 100 AS patients and 100 healthy controls (HCs) were tested using a targeted bisulfite sequencing assay. To verify the results of DNA methylation, mRNA levels of ERAP1 were measured in 20 AS patients and HCs used quantitative real-time reverse transcription-polymerase chain reaction.

RESULTS

The DNA methylation levels of two CpG islands containing 31 loci in ERAP1 promoter were measured. ERAP1_1 (P< .001) and ERAP1_2 (P< .001) islands were significantly hypermethylated in AS patients compared with HCs. In the verification study, the mRNA levels of ERAP1 were significantly decreased in AS patients. The ROC curve analysis showed that the sensitivity, specificity and area under curve were 0.717, 0.737, and 0.779 of differential methylated CpG loci of ERAP1 for AS diagnosis. In AS patients, the methylation levels of EARP1 were associated with family history, non-steroidal anti-inflammatory drugs use, X-ray classification, and clinical manifestations.

CONCLUSIONS

Our study demonstrated that the ERAP1 gene is significantly hypermethylated, and mRNA levels of EARP1 decreased, in AS patients. Our findings suggested that the aberrant methylation of ERAP1 promoter may be involved in the pathogenesis of AS and could be considered as a diagnostic tool and therapeutic target of AS.Abbreviations AS: Ankylosing Spondylitis; AUC: Area Under Curve; BASDAI: Bath Ankylosing Spondylitis Disease Activity Index; BASFI: Bath Ankylosing Spondylitis Functional Index; CI: Confidence Interval; CpG: Cytosine-guanine Dinucleotide; CRP: C-reactive Protein; ERAP1: Endoplasmic Reticulum Aminopeptidase 1; ESR: Erythrocyte Sedimentation Rate; EWAS: Epigenome-Wide Association Study; HLA: Human Leukocyte Antigen; OR: Odds Ratio; PCR: Polymerase Chain Reaction; ROC: Receiver Operating Characteristic; NSAIDs: Non-Steroidal Anti-Inflammatory Drugs.

Epigenetic and transcriptional control of interferon-β

The three classes of interferons (IFNs) share the ability to inhibit viral replication, activating cell transcriptional programs that regulate both innate and adaptive responses to viral and intracellular bacterial challenge. Due to their unique potency in regulating viral replication, and their association with numerous autoimmune diseases, the tightly orchestrated transcriptional regulation of IFNs has long been a subject of intense investigation. The protective role of early robust IFN responses in the context of infection with SARS-CoV-2 has further underscored the relevance of these pathways. In this viewpoint, rather than focusing on the downstream effects of IFN signaling (which have been extensively reviewed elsewhere), we will summarize the historical and current understanding of the stepwise assembly and function of factors that regulate IFNβ enhancer activity (the "enhanceosome") and highlight opportunities for deeper understanding of the transcriptional control of the ifnb gene.

Emerging epigenetic targets in rheumatoid arthritis

Rheumatoid arthritis is a complex disorder that is characterized by irreversible and progressive destructions of joints, but its exact etiology remains mainly unknown. The occurrence and the progression of the disease entirely depend on environmental and genetic factors. In recent years, various epigenetic changes involving DNA methylation, histone modification, miRNA, X-chromosome inactivation, bromodomain, sirtuin, and many others were identified that were found to be linked to the activation and the aggressive phenotype in rheumatoid arthritis. Epigenetics is found to be one of the root causes, which brings changes in the heritable phenotype and is not determined by changes in the DNA sequences and understanding these epigenetic mechanisms and the pathogenesis of the disease can help in understanding the disease and various other possible ways for its control and/or prevention. The various epigenetic modification occurring are reversible and can be modulated by drugs, diet, and environmental factors. This article focuses on various epigenetic factors involved in the pathogenesis of rheumatoid arthritis. Further, various epigenetic therapies that might be successful in inhibiting these epigenetic modifications are summarized. Several therapeutic agents alter the epigenetic modifications occurring in various diseases and many of the epigenetic therapies are under pre-clinical and clinical trial. However, exploring these epigenetic prognostic biomarkers would give a broader perspective and provide more ideas and knowledge regarding the process and pathways through which the diseases occur, and also combining various therapeutic agents would show more beneficial and synergistic effects.

PTEN Methylation Promotes Inflammation and Activation of Fibroblast-Like Synoviocytes in Rheumatoid Arthritis

Rheumatoid arthritis (RA) is characterized by a tumor-like expansion of the synovium and subsequent destruction of adjacent articular cartilage and bone. In our previous work we showed that phosphatase and tension homolog deleted on chromosome 10 (PTEN) contributes to the activation of fibroblast-like synoviocytes (FLS) in adjuvant-induced arthritis (AIA), but the underlying mechanism is not unknown. In this study, we show that PTEN is downregulated while DNA methyltransferase (DNMT)1 is upregulated in FLS from RA patients and a rat model of AIA. DNA methylation of PTEN was increased by administration of tumor necrosis factor (TNF)-α in FLS of RA patients, as determined by chromatin immunoprecipitation and methylation-specific PCR. Treatment with the methylation inhibitor 5-azacytidine suppressed cytokine and chemokine release and FLS activation in vitro and alleviated paw swelling in vivo. PTEN overexpression reduced inflammation and activation of FLS via protein kinase B (AKT) signaling in RA, and intra-articular injection of PTEN-expressing adenovirus into the knee of AIA rats markedly reduced inflammation and paw swelling. Thus, PTEN methylation promotes the inflammation and activation of FLS in the pathogenesis of RA. These findings provide insight into the molecular basis of articular cartilage destruction in RA, and indicate that therapeutic strategies that prevent PTEN methylation may an effective treatment.

Type I Interferons in the Pathogenesis and Treatment of Autoimmune Diseases

Type I interferons (IFN-Is) are a very important group of cytokines that are produced by innate immune cells but also act on adaptive immune cells. IFN-Is possess antiviral, antitumor, and anti-proliferative effects, as well are associated with the initiation and maintenance of autoimmune disorders. Studies have shown that aberrantly expressed IFN-Is and/or type I IFN-inducible gene signatures in the serum or tissues of patients with autoimmune disorders are linked to their pathogenesis, clinical manifestations, and disease activity. Type I interferonopathies with mutations in genes impacting the type I IFN signaling pathway have shown symptoms and characteristics similar to those of systemic lupus erythematosus (SLE). Furthermore, both interventions in animal models and clinical trials of therapies targeting the type I IFN signaling pathway have shown efficacy in the treatment of autoimmune diseases. Our review aims to summarize the functions and targeted therapies (as well as clinical trials) of IFN-Is in both adult and pediatric autoimmune diseases, such as SLE, pediatric SLE (pSLE), rheumatoid arthritis (RA), juvenile idiopathic arthritis (JIA), juvenile dermatomyositis (JDM), Sjögren syndrome (SjS), and systemic sclerosis (SSc), discussing the potential abnormal regulation of transcription factors and epigenetic modifications and providing a potential mechanism for pathogenesis and therapeutic strategies for future clinical use.

The DNA methylation Profile of Undifferentiated Arthritis Patients Anticipates their Subsequent Differentiation to Rheumatoid Arthritis

OBJECTIVE

Undifferentiated arthritis (UA) is the term used to cover all the cases of arthritis that do not fit a specific diagnosis. A significant percentage of UA patients progress to rheumatoid arthritis (RA), others to a different definite rheumatic disease, and the rest undergo spontaneous remission. Therapeutic intervention in patients with UA can delay or halt disease progression and its long-term consequences. It is therefore of inherent interest to identify those UA patients with a high probability of progressing to RA who would benefit from early appropriate therapy. We hypothesized that alterations in the DNA methylation profiles of immune cells may inform on the genetically- or environmentally-determined status of patients and potentially discriminate between disease subtypes.

METHODS

In this study, we performed DNA methylation profiling of a UA patient cohort, in which progression into RA occurs for a significant proportion of the patients.

RESULTS

We find differential DNA methylation in UA patients compared to healthy controls. Most importantly, our analysis identifies a DNA methylation signature characteristic of those UA cases that differentiate to RA. We demonstrate that the methylome of peripheral mononuclear cells can be used to anticipate the evolution of UA to RA, and that this methylome is associated with a number of inflammatory pathways and transcription factors. Finally, we design a machine-learning strategy for DNA methylation-based classification that predicts the differentiation of UA patients towards RA.

CONCLUSION

DNA methylation profiling provides a good predictor of UA-to-RA progression to anticipate targeted treatments and improve clinical management.

Association analysis of juvenile idiopathic arthritis genetic susceptibility factors in Estonian patients

Background

Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic condition of childhood. Genetic association studies have revealed several JIA susceptibility loci with the strongest effect size observed in the human leukocyte antigen (HLA) region. Genome-wide association studies have augmented the number of JIA-associated loci, particularly for non-HLA genes. The aim of this study was to identify new associations at non-HLA loci predisposing to the risk of JIA development in Estonian patients.

Methods

We performed genome-wide association analyses in an entire JIA case–control sample (All-JIA) and in a case–control sample for oligoarticular JIA, the most prevalent JIA subtype. The entire cohort was genotyped using the Illumina HumanOmniExpress BeadChip arrays. After imputation, 16,583,468 variants were analyzed in 263 cases and 6956 controls.

Results

We demonstrated nominal evidence of association for 12 novel non-HLA loci not previously implicated in JIA predisposition. We replicated known JIA associations in CLEC16A and VCTN1 regions in the oligoarticular JIA sample. The strongest associations in the All-JIA analysis were identified at PRKG1 (P = 2,54 × 10⁻⁶), LTBP1 (P = 9,45 × 10⁻⁶), and ELMO1 (P = 1,05 × 10⁻⁵). In the oligoarticular JIA analysis, the strongest associations were identified at NFIA (P = 5,05 × 10⁻⁶), LTBP1 (P = 9,95 × 10⁻⁶), MX1 (P = 1,65 × 10⁻⁵), and CD200R1 (P = 2,59 × 10⁻⁵).

Conclusion

This study increases the number of known JIA risk loci and provides additional evidence for the existence of overlapping genetic risk loci between JIA and other autoimmune diseases, particularly rheumatoid arthritis. The reported loci are involved in molecular pathways of immunological relevance and likely represent genomic regions that confer susceptibility to JIA in Estonian patients.

Key Points • Juvenile idiopathic arthritis (JIA) is the most common childhood rheumatic disease with heterogeneous presentation and genetic predisposition. • Present genome-wide association study for Estonian JIA patients is first of its kind in Northern and Northeastern Europe. • The results of the present study increase the knowledge about JIA risk loci replicating some previously described associations, so adding weight to their relevance and describing novel loci. • The study provides additional evidence for the existence of overlapping genetic risk loci between JIA and other autoimmune diseases, particularly rheumatoid arthritis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10067-021-05756-x.

Association of methylation level and transcript level in TRAF5 gene with ankylosing spondylitis: a case-control study

To explore the association between methylation level and transcript level of TNF receptor-associated factor 5 (TRAF5) gene with ankylosing spondylitis (AS) in Chinese Han population. Methylation and mRNA expression level of the TRAF5 gene were tested in 98 patients and 98 healthy controls. Among the 21 CpG sites, methylation levels at eight sites were significantly different between AS patients and healthy controls. However, only three sites remained significantly different after the correction by the Benjamini-Hochberg method. Compared with controls, the CpG island of TRAF5 gene promoter was highly methylated in AS patients, and the relative mRNA expression level of TRAF5 was significantly reduced in AS patients. And the mRNA level was negatively correlated with the methylation level of TRAF5 gene in AS patients (r_s = -0.453, P < 0.001). Subgroup analyses indicated that there was no significant difference in the level of methylation between groups of different status of HLA-B27 and medications in AS patients. Multiple linear regression showed that disease-modifying antirheumatic drugs could reduce methylation levels of AS patients after adjusting for the effects of other drugs. In conclusion, the hypermethylation of the TRAF5 might contribute to the pathogenesis of AS, but many open questions remain.

Evaluation of the Genetic Association and Methylation of Immune Response Pathway Genes with the Risk of Chronic Periodontitis in the Uighur Population

Aim: The aim of this study was to explore the possible associations between single nucleotide polymorphisms (SNPs) and DNA methylation levels of seven genes in the inflammatory response pathway with susceptibility to chronic periodontitis (CP) among the Uighur population of the Xinjiang Autonomous Region of China. Methods: A total of 444 eligible subjects (279 CP patients and 165 healthy controls) were enrolled in the study. Genomic DNA was obtained from gingival tissue for genotyping eight SNPs and performing methylation measurements of seven genes. Results: SNP rs2070745 in the formyl peptide receptor 1 (FPR1) gene achieved statistical significance in a standard allelic association analysis for CP (p = 0.02). The frequency of the rs2070745 minor allele G was higher in the cases than in controls (0.367 vs. 0.291). Additionally, rs2070745 was significantly associated with CP under the dominant genetic model (p = 0.03). Using logistic regression analysis, rs2070745 was found to be consistently associated with CP under the additive dominant model, and this association remained significant after covariates were taken into account [odds ratio (OR) = 1.49 (1.09-2.05), p = 0.014; OR = 1.58 (1.04-2.40), p = 0.031, respectively]. No significant gene-gene interactions were identified. Although we did not find a polymorphism in interleukin 6 (IL6) associated with CP in our study, the methylation level of a CpG island region located within the promoter region of IL6 was significantly less in CP patients compared with controls (p < 0.05). Conclusions: The genetic polymorphism rs2070745 in FPR1 and the methylation level of the promoter region of IL6 might be associated with CP in the Uighur population of China.

Interaction between early-life pet exposure and methylation pattern of ADAM33 on allergic rhinitis among children aged 3-6 years in China

Background

Recent research has pointed out the important roles of epigenetic modifications in the development and persistence of allergic rhinitis (AR), especially in relation to DNA methylation of disease-associated genes. We investigated whether AR susceptibility genes were epigenetically regulated, and whether methylation modulation of these genes in response to early-life environment could be a molecular mechanism underlying the risk for AR onset in a cohort of children aged 3–6 years in China.

Methods

Peripheral blood mononuclear cell (PBMC) samples were collected from 130 children patients, aged 3–6 years and diagnosed with AR; and 154 matched controls to detect promoter methylation in 25 AR susceptibility genes with the MethylTarget approach. Methylation levels were compared for each CpG site, each amplified region, and each gene. In addition, the relationship among DNA methylation, early-life environmental risk factors and AR onset were assessed.

Results

Maternal allergic history (P = 0.0390) and pet exposure (P = 0.0339) were significantly associated with increased AR risk. Differential methylation analyses were successfully performed for 507 CpG sites, 34 amplified regions and 17 genes and significant hypomethylation was observed in the promoter region of ADAM33 in AR patients [multiple test-corrected (FDR) P-value < 0.05]. Spearman correlation analysis revealed that the hypomethylation of ADAM33 was significantly associated with higher eosinophil counts (Spearman’s ρ: − 0.187, P-value = 0.037). According to the results of the multiple regression analysis, after adjusting for cofounders, the interaction of early-life pet exposure with methylation level of ADAM33 increased the risk for AR onset 1.423 times more in children (95% CI = 0.0290–4.109, P-value = 0.005).

Conclusion

This study provides evidence that early-life pet exposure and low methylation level of ADAM33 increase AR risk in children, and the interaction between pet exposure and methylation level of ADAM33 may play an important role in the development of AR.

Epigenome wide association study of response to methotrexate in early rheumatoid arthritis patients

Aim

To identify differentially methylated positions (DMPs) and regions (DMRs) that predict response to Methotrexate (MTX) in early rheumatoid arthritis (RA) patients.

Materials and methods

DNA from baseline peripheral blood mononuclear cells was extracted from 72 RA patients. DNA methylation, quantified using the Infinium MethylationEPIC, was assessed in relation to response to MTX (combination) therapy over the first 3 months.

Results

Baseline DMPs associated with response were identified; including hits previously described in RA. Additionally, 1309 DMR regions were observed. However, none of these findings were genome-wide significant. Likewise, no specific pathways were related to response, nor could we replicate associations with previously identified DMPs.

Conclusion

No baseline genome-wide significant differences were identified as biomarker for MTX (combination) therapy response; hence meta-analyses are required.