The multifaceted functional role of DNA methylation in immune-mediated rheumatic diseases

DEEP-EP: Identification of epigenetic protein by ensemble residual convolutional neural network for drug discovery

Epigenetic proteins (EP) play a role in the progression of a wide range of diseases, including autoimmune disorders, neurological disorders, and cancer. Recognizing their different functions has prompted researchers to investigate them as potential therapeutic targets and pharmacological targets. This paper proposes a novel deep learning-based model that accurately predicts EP. This study introduces a novel deep learning-based model that accurately predicts EP. Our approach entails generating two distinct datasets for training and evaluating the model. We then use three distinct strategies to transform protein sequences to numerical representations: Dipeptide Deviation from Expected Mean (DDE), Dipeptide Composition (DPC), and Group Amino Acid (GAAC). Following that, we train and compare the performance of four advanced deep learning models algorithms: Ensemble Residual Convolutional Neural Network (ERCNN), Generative Adversarial Network (GAN), Convolutional Neural Network (CNN), and Gated Recurrent Unit (GRU). The DDE encoding combined with the ERCNN model demonstrates the best performance on both datasets. This study demonstrates deep learning's potential for precisely predicting EP, which can considerably accelerate research and streamline drug discovery efforts. This analytical method has the potential to find new therapeutic targets and advance our understanding of EP activities in disease.

Epigenetic regulation of 15-lipoxygenase-1 expression in human chondrocytes by promoter methylation

OBJECTIVE AND DESIGN

15-Lipoxygenase-1 (15-LOX-1) catalyzes the biosynthesis of many anti-inflammatory and immunomodulatory lipid mediators and was reported to have protective properties in several inflammatory conditions, including osteoarthritis (OA). This study was designed to evaluate the expression of 15-LOX-1 in cartilage from normal donors and patients with OA, and to determine whether it is regulated by DNA methylation.

METHODS

Cartilage samples were obtained at autopsy from normal knee joints and from OA-affected joints at the time of total knee joint replacement surgery. The expression of 15-LOX-1 was evaluated using real-time polymerase chain reaction (PCR). The role of DNA methylation in 15-LOX-1 expression was assessed using the DNA methyltransferase inhibitor 5-Aza-2'-desoxycytidine (5-Aza-dC). The effect of CpG methylation on 15-LOX-1 promoter activity was evaluated using a CpG-free luciferase vector. The DNA methylation status of the 15-LOX-1 promoter was determined by pyrosequencing.

RESULTS

Expression of 15-LOX-1 was upregulated in OA compared to normal cartilage. Treatment with 5-Aza-dC increased 15-LOX-1 mRNA levels in chondrocytes, and in vitro methylation decreased 15-LOX-1 promoter activity. There was no difference in the methylation status of the 15-LOX-1 gene promoter between normal and OA cartilage.

CONCLUSION

The expression level of 15-LOX-1 was elevated in OA cartilage, which may be part of a repair process. The upregulation of 15-LOX-1 in OA cartilage was not associated with the methylation status of its promoter, suggesting that other mechanisms are involved in its upregulation.

Targeting post-translational modifications of Foxp3: a new paradigm for regulatory T cell-specific therapy

A healthy immune system is pivotal for the hosts to resist external pathogens and maintain homeostasis; however, the immunosuppressive tumor microenvironment (TME) damages the anti-tumor immunity and promotes tumor progression, invasion, and metastasis. Recently, many studies have found that Foxp3+ regulatory T (Treg) cells are the major immunosuppressive cells that facilitate the formation of TME by promoting the development of various tumor-associated cells and suppressing the activity of effector immune cells. Considering the role of Tregs in tumor progression, it is pivotal to identify new therapeutic drugs to target and deplete Tregs in tumors. Although several studies have developed strategies for targeted deletion of Treg to reduce the TME and support the accumulation of effector T cells in tumors, Treg-targeted therapy systematically affects the Treg population and may lead to the progression of autoimmune diseases. It has been understood that, nevertheless, in disease conditions, Foxp3 undergoes several definite post-translational modifications (PTMs), including acetylation, glycosylation, phosphorylation, ubiquitylation, and methylation. These PTMs not only elevate or mitigate the transcriptional activity of Foxp3 but also affect the stability and immunosuppressive function of Tregs. Various studies have shown that pharmacological targeting of enzymes involved in PTMs can significantly influence the PTMs of Foxp3; thus, it may influence the progression of cancers and/or autoimmune diseases. Overall, this review will help researchers to understand the advances in the immune-suppressive mechanisms of Tregs, the post-translational regulations of Foxp3, and the potential therapeutic targets and strategies to target the Tregs in TME to improve anti-tumor immunity.

Epigenome-wide association study identifies Vogt-Koyanagi-Harada disease-associated DNA methylation loci in Chinese

DNA methylation is one of the important epigenetic mechanisms for modulating gene expression. By performing a genome-wide methylation association analysis of whole peripheral blood from 60 Vogt-Koyanagi-Harada disease (VKH) patients and 60 healthy controls, we depicted the global DNA methylation status of VKH disease. Further pyrosequencing validation in 160 patients and 159 controls identified 3 aberrant CpG sites in HLA gene regions including cg04026937 and cg18052547 (located in HLA-DRB1 region), and cg13778567 (HLA-DQA1). We also identified 9 aberrant CpG sites in non-HLA gene regions including cg13979407, cg21075643, cg24290586, cg10135747 and cg22707857 (BTNL2), cg22155039 (NOTCH4), cg02605387 (TNXB), cg06255004 (AGPAT2) and cg18855195 (RIBC2). Increased mRNA levels of BTNL2, NOTCH4 and TNXB were identified in VKH patients when compared with healthy controls, consistent with the hypomethylated CpG status in these gene regions. Moreover, seven aberrantly methylated CpG sites may serve as a diagnostic marker for VKH disease (AUC = 84.95%, 95%CI: 79.49%-90.41%).

DNA Methylation in Autoimmune Thyroid Disease

Graves disease and Hashimoto disease form part of the spectrum of autoimmune thyroid disease (AITD), to which genetic and environmental factors are recognized contributors. Epigenetics provides a potential link between environmental influences, gene expression, and thyroid autoimmunity. DNA methylation (DNAm) is the best studied epigenetic process, and global hypomethylation of leukocyte DNA is reported in several autoimmune disorders. This review summarizes the current understanding of DNAm in AITD. Targeted DNAm studies of blood samples from AITD patients have reported differential DNAm in the promoter regions of several genes implicated in AITD, including TNF, IFNG, IL2RA, IL6, ICAM1, and PTPN22. In many cases, however, the findings await replication and are unsupported by functional studies to support causal roles in AITD pathogenesis. Furthermore, thyroid hormones affect DNAm, and in many studies confounding by reverse causation has not been considered. Recent studies have shown that DNAm patterns in candidate genes including ITGA6, PRKAA2, and DAPK1 differ between AITD patients from regions with different iodine status, providing a potential mechanism for associations between iodine and AITD. Research focus in the field is moving from candidate gene studies to an epigenome-wide approach. Genome-wide methylation studies of AITD patients have demonstrated multiple differentially methylated positions, including some in immunoregulatory genes such as NOTCH1, HLA-DRB1, TNF, and ICAM1. Large, epigenome-wide studies are required to elucidate the pathophysiological role of DNAm in AITD, with the potential to provide novel diagnostic and prognostic biomarkers as well as therapeutic targets.

DNA Methylation Expression Profile of Blood Heat Syndrome and Blood Stasis Syndrome in TCM Psoriasis

Objective

Traditional Chinese medicine (TCM) emphasizes treatment based on syndrome differentiation. This study aimed to clarify the characteristics of DNA methylation expression profiles in peripheral blood mononuclear cells (PBMCs) in patients with psoriasis and analyze the differences in these profiles among different TCM syndromes of psoriasis in order to provide a material basis for the diversity of these syndromes.

Methods

Blood samples were collected from 32 participants, including 14 patients with psoriatic blood heat syndrome (BHS), 12 patients with psoriatic blood stasis syndrome (BSS), and 6 healthy controls. PBMCs were extracted and subjected to DNA quality inspection. An Illumina Human Methylation 850k chip was used to sequence each group of samples. According to gene annotation classification together with CpG island annotation classification, the differentially methylated regions between sample groups were screened, while Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were applied to perform functional analyses of DMGs. Finally, the DMGs closely correlating with psoriatic severity were screened using Spearman's correlation analysis.

Results

Compared with normal controls, patients with psoriasis showed an overall trend of hypermethylation. In psoriasis, the differential methylation probes were mainly distributed on gene body region on the genome, while those in CpG regions were mainly distributed in CpG islands. Compared with healthy controls, the overall trends in methylation were similar in psoriatic BHS and BSS patients compared to healthy controls. However, bioinformatic analysis revealed different functions of DMGs. We also found that the methylation levels of TRIM14 and PRDM16 were closely correlated with PASI scores and could serve as potential biomarkers to assess the severity of psoriasis.

Conclusions

Our study, for the first time, indicated the possible involvement of DNA methylation in regulating the characteristics of TCM syndromes of psoriasis, providing a new direction for research into TCM psoriatic syndromes.

The hypermethylation of FOXP3 gene as an epigenetic marker for the identification of arsenic poisoning risk

Background and Purpose: Arsenic exposure can lead to skin lesions and multiple organ damage, which are not easily reversible and for which there is no effective therapeutics. Identification of reliable epigenetic markers is essential for early recognition of arsenic poisoning risk. Anomalous DNA methylation of immune homeostasis regulator FOXP3 is a critical mechanism for triggering arsenic poisoning. This study aims to explore the value of FOXP3 methylation in the identification of arsenic poisoning risk.Methods: 88 arsenic poisoning subjects and 41 references were recruited. Urinary arsenic contents and FOXP3 methylation in PBLCs was measured by ICP-MS and pyrosequencing, respectively.Results: The results showed that the elevated FOXP3 methylation in PBLCs were associated with the increased levels of urinary arsenic and were positively associated with the increased risk of arsenic poisoning and its progression. The result of mediation analysis revealed that 24.3% of the effect of arsenic exposure on the risk of arsenic poisoning was mediated by increased FOXP3 methylation. Additionally, we constructed a nomogram model with FOXP3 methylation as an epigenetic predictor to assess the probability of individual arsenic poisoning. The model showed a robust ability in the discrimination of arsenic poisoning risk, with an area under receiver operating characteristics curve of 0.897(0.845-0.949) and more than 70% accuracy. The calibration curves and the Harrell concordance index showed that the consistency rate between the probability predicted by the nomogram model and the actual probability is 89.7%.Conclusions: Taken together, we found the great potential of FOXP3 methylation for the identification of arsenic poisoning risk and provided a new approach to the application of epigenetic markers in accurately quantifying the risk of adverse outcomes.

Comprehensive analysis of epigenetic modifications and immune-cell infiltration in tissues from patients with systemic lupus erythematosus

Aim: To explore potential abnormal epigenetic modifications and immune-cell infiltration in tissues from systemic lupus erythematosus (SLE) patients. Materials & methods: To utilize bioinformatics analysis and 'wet lab' methods to identify and verify differentially expressed genes in multiple targeted organs in SLE. Results: Seven key genes, IFI44, IFI44L, IFIT1, IFIT3, PLSCR1, RSAD2 and OAS2, which are regulated by epigenetics and may be involved in the pathogenesis of SLE, are identified by combined long noncoding RNA-miRNA-mRNA network analysis and DNA methylation analysis. The results of quantitative reverse transcription PCR, immunohistochemistry and DNA methylation analysis confirmed the potential of these genes as biomarkers. Conclusion: This study reveals the potential mechanisms in SLE from epigenetic modifications and immune-cell infiltration, providing diagnostic biomarkers and therapeutic targets for SLE.

Epigenome association study for DNA methylation biomarkers in buccal and monocyte cells for female rheumatoid arthritis

Genetics (i.e., mutations) has been assumed to be the major factor in rheumatoid arthritis (RA) etiology, but accounts for a minority of the variance in disease risk for RA. In contrast to genetics, the environment can have dramatic impacts on epigenetics that associate with disease etiology. The current study used buccal cells and purified blood monocytes from two different clinical cohorts involving Caucasian or African American female populations with or without arthritis. The differential DNA methylation regions (DMRs) between the control and RA populations were identified with an epigenome-wide association study. The DMRs (i.e., epimutations) identified in the buccal cells and monocytes were found to be distinct. The DMR associated genes were identified and many have previously been shown to be associated with arthritis. Observations demonstrate DNA methylation epimutation RA biomarkers are cell type specific and similar findings were observed with the two racial background populations. Rheumatoid arthritis susceptibility epigenetic diagnosis appears feasible and may improve the clinical management of RA and allowpreventative medicine considerations.

Advanced genomics and clinical phenotypes in psoriatic arthritis

Psoriatic Arthritis (PsA) is a complex polygenic inflammatory disease showing a variable musculoskeletal involvement in patients with skin psoriasis. PsA coexist in 25-40 % of patients with the dermatological manifestations, but PsA may also predate the appearance of psoriasis. Nonetheless, the immunopathogenesis of psoriasis and PsA manifest significant similarities, with a major role of the individual susceptibility in both cases. Genome wide association studies (GWAS) identified several genes/loci associated with the risk to develop PsA, both dependent and independent of psoriasis. The major challenge is thus represented by the need to translate the identification of functional polymorphisms and other genetics findings into biological mechanisms along with the identification of novel putative drug targets. A functional genomics approach aims to increase GWAS power and recent evidence supports the use of a multilayer process, including eQTL, methylome, chromatin conformation analysis and genome editing to discover novel genes that can be affected by disease-associated variants, such as PsA. The available data have considered PsA as a unique homogeneous clinical entity while the clinical experience supports a wide variability of skin and joint manifestations coexisting in diverse patients with different mechanisms underlying the musculoskeletal and dermatological domains. A better discrimination of the patient features is encouraged by the limited data on functional genomics. We provide herein a review of the latest findings on PsA functional genomics highlighting the exciting developments in the field and how these might lead to a better understanding of gene regulation underpinning disease mechanisms and ultimately refine clinical phenotyping.

A Variant of sNASP Exacerbates Lymphocyte Subset Disorder and Nephritis in a Spontaneous Lupus Model Sle1.Yaa Mouse

A variant of somatic nuclear autoantigenic sperm protein (sNASP) was identified from the murine lupus susceptibility locus Sle2c1 by whole exome sequencing (WES). Previous studies have shown that mutant sNASP could synergize with the Fas^lpr mutation in exacerbating autoimmunity and aggravating end-organ inflammation. In the current study, the sNASP mutation was introduced into Sle1.Yaa mice to detect whether it has a synergistic effect with Sle1 or Yaa loci. As expected, compared with Sle1.Yaa mice, Sle1.Yaa.ΔsNASP mice showed enlarged lymph nodes, aggravated renal inflammation, and shortened survival time. The proportions of CD3⁺ T cells, activated CD19⁺CD86⁺ B cells, Th1 cells in the spleen and lymph nodes, and Th17 cells in lymph nodes in Sle1.Yaa.ΔsNASP mice were increased compared to those in Sle1.Yaa mice. The levels of IFN-γ and TNF-α in the serum of Sle1.Yaa.ΔsNASP mice were higher than those of Sle1.Yaa mice. The above results show that mutant sNASP can interact with different lupus susceptibility genes and promote the disease process of systemic lupus erythematosus.

Identification of circular RNAs hsa_circ_0140271 in peripheral blood mononuclear cells as a novel diagnostic biomarker for female rheumatoid arthritis

Background

Rheumatoid arthritis (RA) is a chronic autoimmune disease, which commonly affects women. Accumulating evidence shows that differentially expressed circular RNAs (circRNAs) play crucial roles in the progress of RA. However, the roles of circRNAs in female RA remains unclear. This study explores potential role and diagnostic value of hsa_circ_0140271 from peripheral blood mononuclear cells (PBMC) in female RA.

Methods

Differential expression of circRNAs was determined by RNA-sequencing in PBMC from 4 healthy controls (HC) and 4 RA patients, and we further measured the level of hsa_circ_0140271 in a validation cohort consisting of 47 RA and 47 HC via RT-qPCR. Besides, correlation studies with clinical variables were also examined. What’s more, we performed bioinformatics analysis to predict the potential role of hsa_circ_0140271.

Results

PBMC expression of hsa_circ_0140271 of female RA was significantly higher than that of female HC, and it was positively correlated with antistreptolysin (ASO). Furthermore, the receiver operating characteristic (ROC) curve indicated that hsa_circ_0140271 could distinguish female RA from female HC and female patients with ankylosing spondylitis (AS) or osteoarthritis (OA). Besides, the combined diagnosis anti-cyclic citrullinated peptide (Anti-CCP) + hsa_circ_0140271 could improve diagnostic accuracy with an area under the curve (AUC) of 0.818 to compared with Anti-CCP. Furthermore, KEGG pathway enrichment analysis indicated hsa_circ_0140271 may act as microRNA sponge and participate in fatty acid metabolism pathways.

Conclusion

Hsa_circ_0140271 was likely to be used as a promising diagnostic biomarker for female RA; it may act as microRNA sponge to regulate fatty acid metabolism pathways in RA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13018-021-02794-8.

Ten Years of EWAS

Epigenome-wide association study (EWAS) has been applied to analyze DNA methylation variation in complex diseases for a decade, and epigenome as a research target has gradually become a hot topic of current studies. The DNA methylation microarrays, next-generation, and third-generation sequencing technologies have prepared a high-quality platform for EWAS. Here, the progress of EWAS research is reviewed, its contributions to clinical applications, and mainly describe the achievements of four typical diseases. Finally, the challenges encountered by EWAS and make bold predictions for its future development are presented.

Long non‑coding RNA expression profiles identify lncRNA‑XLOC_I2_006631 as a potential novel blood biomarker for Hashimoto's thyroiditis

Long non‑coding RNAs (lncRNAs) have been increasingly recognized as important immune checkpoints involved in the pathogenesis of autoimmune diseases. However, the exact role of lncRNAs in Hashimoto's thyroiditis (HT) has been rarely studied. The aim of the present study was to investigate the role of lncRNAs and the potential biomarkers in HT, a total of 33 patients with HT and 32 healthy volunteers were enrolled in the present study, and five patients and five healthy controls were investigated using next generation sequencing. A total of 218 dysregulated lncRNAs, including 94 upregulated and 124 downregulated lncRNAs, were identified and examined in the peripheral blood mononuclear cells (PBMCs) from patients with HT. The majority of the lncRNAs were intergenic and exonic (66.06%). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis demonstrated that abnormally expressed lncRNAs were enriched in the 'NF‑kB expression', in the 'TGF‑β signaling pathway' and in the 'JAK‑STAT signaling pathway', which are associated with the immunopathogenic mechanisms of HT. In total, three lncRNAs (LOC729737, XLOC_I2_006631 and BC041964) were validated and had a trend identical to that detected by the sequencing results. The expression of lncRNA‑XLOC_I2_006631 was upregulated and was positively correlated with the serum concentrations of anti‑thyroperoxidase antibody in patients with HT. Methyl‑CpG‑binding protein 2 (MECP2) was identified as the potential regulatory gene of lncRNA‑XLOC_I2_006631 using a prediction program. The expression of MECP2 was increased and was positively correlated with the elevated expression levels of lncRNA‑XLOC_I2_006631 and anti‑thyroperoxidase antibody in patients with HT. Furthermore, lncRNA‑XLOC_I2_006631 was able to regulate MECP2 expression in vitro. Receiver operating characteristic curve analysis suggested that lncRNA‑XLOC_I2_006631 has a potential diagnostic value. Collectively, the present results indicated the important role of dysregulated lncRNAs in HT and demonstrated that lncRNA‑XLOC_I2_006631 functioned as a positive regulator of MECP2 expression, suggesting a potential mechanism. Thus, lncRNA‑XLOC_I2_006631 may be used as a biomarker of HT.