Identification of Cell-Specific Differential DNA Methylation Associated With Methotrexate Treatment Response in Rheumatoid Arthritis

DNA and RNA Methylation in Rheumatoid Arthritis-A Narrative Review

Rheumatoid arthritis (RA) is a progressive autoimmune disease leading to structural and functional joint damage and, eventually, to physical disability. The pathogenesis of the disease is highly complex and involves interactions between fibroblast-like synoviocytes (FLSs) and immune cells, which stimulate the secretion of pro-inflammatory factors, leading to chronic inflammation. In recent years, studies have demonstrated the importance of epigenetics in RA. Specifically, epigenetic alterations have been suggested to serve as diagnostic and treatment biomarkers, while epigenetic mechanisms are thought to be involved in the pathogenesis of RA. Epigenetic regulators coordinate gene expression, and in the case of inflammatory diseases, they regulate the expression of a broad range of inflammatory molecules. In this review, we discuss current evidence on the involvement of DNA and RNA methylation in RA.

DMRIntTk: Integrating different DMR sets based on density peak clustering

BACKGROUND

Identifying differentially methylated regions (DMRs) is a basic task in DNA methylation analysis. However, due to the different strategies adopted, different DMR sets will be predicted on the same dataset, which poses a challenge in selecting a reliable and comprehensive DMR set for downstream analysis.

RESULTS

Here, we develop DMRIntTk, a toolkit for integrating DMR sets predicted by different methods on a same dataset. In DMRIntTk, the genome is segmented into bins, and the reliability of each DMR set at different methylation thresholds is evaluated. Then, the bins are weighted based on the covered DMR sets and integrated into final DMRs using a density peak clustering algorithm. To demonstrate the practicality of DMRIntTk, it was applied to different scenarios, including tissues with relatively large methylation differences, cancer tissues versus normal tissues with medium methylation differences, and disease tissues versus normal tissues with subtle methylation differences. Our results show that DMRIntTk can effectively trim regions with small methylation differences from the original DMR sets and thereby enriching the proportion of DMRs with larger methylation differences. In addition, the overlap analysis suggests that the integrated DMR sets are quite comprehensive, and functional analyses indicate the integrated disease-related DMRs are significantly enriched in biological pathways associated with the pathological mechanisms of the diseases. A comparative analysis of the integrated DMR set versus each original DMR set further highlights the superiority of DMRIntTk, demonstrating the unique biological insights it can provide.

CONCLUSIONS

Conclusively, DMRIntTk can help researchers obtain a reliable and comprehensive DMR set from many prediction methods.

Pre-diagnosis blood DNA methylation profiling of twin pairs discordant for breast cancer points to the importance of environmental risk factors

BACKGROUND

Assessment of breast cancer (BC) risk generally relies on mammography, family history, reproductive history, and genotyping of major mutations. However, assessing the impact of environmental factors, such as lifestyle, health-related behavior, or external exposures, is still challenging. DNA methylation (DNAm), capturing both genetic and environmental effects, presents a promising opportunity. Previous studies have identified associations and predicted the risk of BC using DNAm in blood; however, these studies did not distinguish between genetic and environmental contributions to these DNAm sites. In this study, associations between DNAm and BC are assessed using paired twin models, which control for shared genetic and environmental effects, allowing testing for associations between DNAm and non-shared environmental exposures and behavior.

RESULTS

Pre-diagnosis blood samples of 32 monozygotic (MZ) and 76 dizygotic (DZ) female twin pairs discordant for BC were collected at the mean age of 56.0 years, with the mean age at diagnosis 66.8 years and censoring 75.2 years. We identified 212 CpGs (p < 6.4*10-8) and 15 DMRs associated with BC risk across all pairs using paired Cox proportional hazard models. All but one of the BC risks associated with CpGs were hypomethylated, and 198/212 CpGs had their DNAm associated with BC risk independent of genetic effects. According to previous literature, at least five of the top CpGs were related to estrogen signaling. Following a comprehensive two-sample Mendelian randomization analysis, we found evidence supporting a dual causal impact of DNAm at cg20145695 (gene body of NXN, rs480351) with increased risk for estrogen receptor positive BC and decreased risk for estrogen receptor negative BC.

CONCLUSION

While causal effects of DNAm on BC risk are rare, most of the identified CpGs associated with the risk of BC appear to be independent of genetic effects. This suggests that DNAm could serve as a valuable biomarker for environmental risk factors for BC, and may offer potential benefits as a complementary tool to current risk assessment procedures.

Characteristics of T-Cell Receptor Repertoire for Differential Response to Methotrexate Treatment for Rheumatoid Arthritis

BACKGROUND

Methotrexate (MTX) serves as the initial treatment for rheumatoid arthritis (RA). However, a substantial proportion of RA patients, estimated between 30% and 50%, do not respond positively to MTX. While the T-cell receptor (TCR) is crucial for the immune response during RA, its role in differentiating MTX responsiveness has not been thoroughly investigated.

METHODS

This study used next-generation sequencing to analyze the TCR β-chain complementary determining region sequences in peripheral blood mononuclear cells obtained from RA patients before MTX treatment. This study aimed to compare the characteristics of the TCR repertoire between the MTX responder and non-responder groups.

RESULTS

The study identified a significant difference in the TRBV6-6 gene (p = .003) concerning MTX treatment response. Additionally, a significant difference was found in the number of "3" nucleotide deletions at the 5'Jdels site (p = .023) in the VDJ rearrangement.

CONCLUSION

These findings suggest distinct TCR repertoire characteristics between MTX responder and non-responder groups among RA patients. This discovery offers new insights into understanding the variable responses of RA patients to MTX therapy.

Single-cell RNA-Seq analysis reveals cell subsets and gene signatures associated with rheumatoid arthritis disease activity

Rheumatoid arthritis (RA) management leans toward achieving remission or low disease activity. In this study, we conducted single-cell RNA sequencing (scRNA-Seq) of peripheral blood mononuclear cells (PBMCs) from 36 individuals (18 patients with RA and 18 matched controls, accounting for age, sex, race, and ethnicity), to identify disease-relevant cell subsets and cell type-specific signatures associated with disease activity. Our analysis revealed 18 distinct PBMC subsets, including an IFN-induced transmembrane 3-overexpressing (IFITM3-overexpressing) IFN-activated monocyte subset. We observed an increase in CD4+ T effector memory cells in patients with moderate-high disease activity (DAS28-CRP ≥ 3.2) and a decrease in nonclassical monocytes in patients with low disease activity or remission (DAS28-CRP < 3.2). Pseudobulk analysis by cell type identified 168 differentially expressed genes between RA and matched controls, with a downregulation of proinflammatory genes in the γδ T cell subset, alteration of genes associated with RA predisposition in the IFN-activated subset, and nonclassical monocytes. Additionally, we identified a gene signature associated with moderate-high disease activity, characterized by upregulation of proinflammatory genes such as TNF, JUN, EGR1, IFIT2, MAFB, and G0S2 and downregulation of genes including HLA-DQB1, HLA-DRB5, and TNFSF13B. Notably, cell-cell communication analysis revealed an upregulation of signaling pathways, including VISTA, in both moderate-high and remission-low disease activity contexts. Our findings provide valuable insights into the systemic cellular and molecular mechanisms underlying RA disease activity.

An Elevated IL10 mRNA Combined with Lower TNFA mRNA Level in Active Rheumatoid Arthritis Peripheral Blood

We aimed to investigate the expression of pro-inflammatory cytokine genes TNFA, IL6, IL12B, IL23, IL18 and immunoregulatory genes FOXP3, TGFB1, and IL10 in the peripheral blood of patients with rheumatoid arthritis (RA) at messenger ribonucleic acid (mRNA) level. The total RNA was isolated from peripheral blood samples. Real-time quantitative PCR was used to perform TaqMan-based assays to quantify mRNAs from 8 target genes. IL23A was upregulated (1.7-fold), whereas IL6 (5-fold), FOXP3 (4-fold), and IL12B (2.56-fold) were downregulated in patients compared to controls. In addition, we found a strong positive correlation between the expression of FOXP3 and TNFA and a moderate correlation between FOXP3 and TGFB1. These data showed the imbalance of the T helper (Th) 1/Th17/ T regulatory (Treg) axis at a systemic level in RA. In cases with active disease, the IL10 gene expression was approximately 2-fold higher; in contrast, the expression of FOXP3 was significantly decreased (3.38-fold). The main part of patients with higher disease activity expressed upregulation of IL10 and downregulation of TNFA. Different disease activity cohorts could be separated based on IL10, TNFA and IL12B expression combinations. In conclusion, our results showed that active disease is associated with an elevated IL10 and lower TNFA mRNA level in peripheral blood cells of RA patients.

Role of IFN-α in Rheumatoid Arthritis

PURPOSE OF REVIEW

Type 1 interferons (IFN-I) are of increasing interest across a wide range of autoimmune rheumatic diseases. Historically, research into their role in rheumatoid arthritis (RA) has been relatively neglected, but recent work continues to highlight a potential contribution to RA pathophysiology.

RECENT FINDINGS

We emphasise the importance of disease stage when examining IFN-I in RA and provide an overview on how IFN-I may have a direct role on a variety of relevant cellular functions. We explore how clinical trajectory may be influenced by increased IFN-I signalling, and also, the limitations of scores composed of interferon response genes. Relevant environmental triggers and inheritable RA genetic risk relating to IFN-I signalling are explored with emphasis on intriguing data potentially linking IFN-I exposure, epigenetic changes, and disease relevant processes. Whilst these data cumulatively illustrate a likely role for IFN-I in RA, they also highlight the knowledge gaps, particularly in populations at risk for RA, and suggest directions for future research to both better understand IFN-I biology and inform targeted therapeutic strategies.

Cell-type specific EWAS identifies genes involved in HIV pathogenesis and oncogenesis among people with HIV infection

Epigenome-wide association studies (EWAS) of heterogenous blood cells have identified CpG sites associated with chronic HIV infection, which offer limited knowledge of cell-type specific methylation patterns associated with HIV infection. Applying a computational deconvolution method validated by capture bisulfite DNA methylation sequencing, we conducted a cell type-based EWAS and identified differentially methylated CpG sites specific for chronic HIV infection among five immune cell types in blood: CD4+ T-cells, CD8+ T-cells, B cells, Natural Killer (NK) cells, and monocytes in two independent cohorts (N total =1,134). Differentially methylated CpG sites for HIV-infection were highly concordant between the two cohorts. Cell-type level meta-EWAS revealed distinct patterns of HIV-associated differential CpG methylation, where 67% of CpG sites were unique to individual cell types (false discovery rate, FDR <0.05). CD4+ T-cells had the largest number of HIV-associated CpG sites (N=1,472) compared to any other cell type. Genes harboring statistically significant CpG sites are involved in immunity and HIV pathogenesis (e.g. CX3CR1 in CD4+ T-cells, CCR7 in B cells, IL12R in NK cells, LCK in monocytes). More importantly, HIV-associated CpG sites were overrepresented for hallmark genes involved in cancer pathology ( FDR <0.05) (e.g. BCL family, PRDM16, PDCD1LGD, ESR1, DNMT3A, NOTCH2 ). HIV-associated CpG sites were enriched among genes involved in HIV pathogenesis and oncogenesis such as Kras-signaling, interferon-α and -γ, TNF-α, inflammatory, and apoptotic pathways. Our findings are novel, uncovering cell-type specific modifications in the host epigenome for people with HIV that contribute to the growing body of evidence regarding pathogen-induced epigenetic oncogenicity, specifically on HIV and its comorbidity with cancers.