IL1β, IL18, NFKB1 and IFNG gene interactions are associated with severity of rheumatoid arthritis: A pilot study

Integration of single-cell and bulk RNA sequencing revealed immune heterogeneity and its association with disease activity in rheumatoid arthritis patients

Rheumatoid arthritis (RA) is a chronic, inflammatory, systemic autoimmune disease characterized by cartilage, bone damage, synovial inflammation, hyperplasia, autoantibody production, and systemic features. To obtain an overall profile of the immune environment in RA patients and its association with clinical features, we performed single-cell transcriptome and T-cell receptor sequencing of mononuclear cells from peripheral blood (PBMC) and synovial fluid (SF) from RA patients, integrated with two large cohorts with bulk RNA sequencing for further validation and investigation. Dendritic cells (DCs) exhibited relatively high functional heterogeneity and tissue specificity in relation to both antigen presentation and proinflammatory functions. Peripheral helper T cells (TPHs) are likely to originate from synovial tissue, undergo activation and exhaustion, and are subsequently released into the peripheral blood. Notably, among all immune cell types, TPHs were found to have the most intense associations with disease activity. In addition, CD8 effector T cells could be clustered into two groups with different cytokine expressions and play distinct roles in RA development. By integrating single-cell data with bulk sequencing from two large cohorts, we identified interactions among TPHs, CD8 cells, CD16 monocytes, and DCs that strongly contribute to the proinflammatory local environment in RA joints. Of note, the swollen 28-joint counts exhibited a more pronounced association with this immune environment compared to other disease activity indexes. The immune environment alternated significantly from PBMCs to SF, which indicated that a series of immune cells was involved in proinflammatory responses in the local joints of RA patients. By integrating single-cell data with two large cohorts, we have uncovered associations between specific immune cell populations and clinical features. This integration provides a rapid and precise methodology for assessing local immune activation, offering valuable insights into the pathophysiological mechanisms at play in RA.

Mechanisms of Tripterygium wilfordii Hook F on treating rheumatoid arthritis explored by network pharmacology analysis and molecular docking

Background

Rheumatoid arthritis (RA) is a chronic inflammatory and disabling disease that imposes significant economic and social costs. Tripterygium wilfordii Hook F (TwHF) has a long history of use in traditional Chinese medicine for treating joint disorders, and it has been shown to be cost-effective in treating RA, but its exact mechanism is unknown.

Objective

The goal of the network pharmacology analysis and molecular docking was to investigate the potential active compounds and associated anti-RA mechanisms of TwHF.

Methods

TCMSP and UniProt databases were searched for active compounds and related targets of TwHF. PharmGKB, DrugBank, OMIM, TTD, and the Human Gene Databases were used to identify RA-related targets. The intersected RA and TwHF targets were entered into the STRING database to create a protein-protein interaction network. R software was used for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Molecular docking technology was used to analyze the optimal effective components from TwHF for docking with the selected target gene.

Results

Following screening and duplicate removal, a total of 51 active compounds and 96 potential targets were chosen. The PPI network revealed that the target proteins are CXCL8, CXCL6, STAT3, STAT1, JUN, PPARG, TP53, IL14, MMP9, VEGFA, RELA, CASP3, PTGS2, IFNG, AKT1, FOS, ICAM1, and MAPK14. The results of the GO enrichment analysis focused primarily on the response to lipopolysaccharide, the response to molecules of bacterial origin, and the response to drugs. The KEGG results indicated that the mechanisms were closely related to lipid and atherosclerosis, chemical carcinogenesis-receptor activation, Kaposi sarcoma-associated, herpesvirus infection, hepatitis B, fluid shear stress and atherosclerosis, IL-17 signaling pathways, Th17-cell differentiation, and so on, all of which are involved in angiogenesis, immune cell chemotaxis, and inflammatory responses. Molecular docking results suggested that triptolide was the appropriate PTGS1, PTGS2, and TNF inhibitors.

Conclusion

Our findings provide an essential role and basis for further immune inflammatory studies into the molecular mechanisms of TwHF and PTGS1, PTGS2, and TNF inhibitor development in RA.

Interferon-gamma (IFN-γ) intronic variant (rs2430561) is a risk factor for systemic lupus erythematosus: Observation from a meta-analysis

BACKGROUND

The role of interferon-gamma (IFN-γ) in autoimmune disorders has been well documented. Elevated levels of IFN-γ are observed in rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) and are linked with disease severity. Single nucleotide polymorphism in the intronic region of the IFN-γ gene (+874 T>A rs2430561) has been associated with susceptibility to the development of RA and SLE; however, the reports remained contradictories. We conducted a meta-analysis using earlier published articles to reach a valid conclusion on the role of IFN-γ polymorphism (+874 T>A) in autoimmune diseases.

MATERIALS AND METHODS

Various online databases such as PubMed, Google Scholar, Science Direct, and Scopus were searched to find eligible reports for inclusion in the present analysis. Two independent authors extracted eligible studies and data. The meta-analysis was performed by comprehensive meta-analysis software (CMA) v.3.1. Trial sequential analysis was performed to test whether enough case-control studies have already been conducted worldwide to reach a valid observation.

RESULTS

Six published reports on the role of IFN-γ +874 T>A in SLE and four in RA were found after searching various databases. However, out of those six studies in SLE, in one study, the distribution of genotypes was not following the hardy-Weinberg equilibrium. In RA, three studies were deviated out of four reports. Thus, a total of five studies comprising 1440 SLE patients and 1748 controls were considered for the present meta-analysis. Meta-analysis showed a significant association between IFN-γ +874 T>A variants with susceptibility to SLE (homozygous comparison: p = 0.036, OR = 1.592, heterozygous model: p = 0.042, OR = 1.507, dominant model: p = 0.002, OR = 1.309).

CONCLUSIONS

IFN-γ +874 T>A variant is associated with predisposition to SLE development.

Analysis of the Whole-Genome Sequences from an Equus Parent-Offspring Trio Provides Insight into the Genomic Incompatibilities in the Hybrid Mule

Interspecific hybridization often shows negative effects on hybrids. However, only a few multicellular species, limited to a handful of plants and animals, have shown partial genetic mechanisms by which hybridization leads to low fitness in hybrids. Here, to explore the outcome of combining the two genomes of a horse and donkey, we analyzed the whole-genome sequences from an Equus parent-offspring trio using Illumina platforms. We generated 41.39× and 46.21× coverage sequences for the horse and mule, respectively. For the donkey, a 40.38× coverage sequence was generated and stored in our laboratory. Approximately 24.86 million alleles were discovered that varied from the reference genome. Single nucleotide polymorphisms were used as polymorphic markers for assigning alleles to their parental genomic inheritance. We identified 25,703 Mendelian inheritance error single nucleotide polymorphisms in the mule genome that were not inherited from the parents through Mendelian inheritance. A total of 555 de novo single nucleotide polymorphisms were also identified. The rate of de novo single nucleotide polymorphisms was 2.21 × 10^-7 in the mule from the Equus parent-offspring trio. This rate is obviously higher than the natural mutation rate for Equus, which is also consistent with the previous hypothesis that interracial crosses may have a high mutation rate. The genes associated with these single nucleotide polymorphisms are mainly involved in immune processes, DNA repair, and cancer processes. The results of the analysis of three genomes from an Equus parent-offspring trio improved our knowledge of the consequences of the integration of parental genomes in mules.

A Preliminary Inquiry Into the Potential Mechanism of Huang-Lian-Jie-Du Decoction in Treating Rheumatoid Arthritis via Network Pharmacology and Molecular Docking

Huang-Lian-Jie-Du decoction (HLJDD) has been widely applied to treat inflammation-associated diseases for thousands of years in China. However, the concrete molecular mechanism of HLJDD in the treatment of rheumatoid arthritis (RA) remains unclear. In this work, network pharmacology and molecular docking were applied to preliminarily analyze the potential active ingredients, drug targets, and related pathways of HLJDD on treating RA. A total of 102 active compounds with corresponding 189 targets were identified from HLJDD, and 41 common targets were further identified by intersecting with RA-related targets. Functional enrichment analysis was performed to screen the biological pathways associated with RA. Ten hub targets were further identified through constructing the protein–protein interaction (PPI) network of common targets, which were mainly enriched in the interleukin-17 (IL-17) signaling pathway, tumor necrosis factor (TNF) signaling pathway, and Toll-like receptor signaling pathway. Furthermore, a complex botanical drugs-ingredients-hub-targets-disease network was successfully constructed. The molecular docking results exhibited that these vital ingredients of HLJDD had a stable binding to the hub targets. Among these ingredients, quercetin (MOL000098) was the most common molecule with stable binding to all the targets, and PTGS2 was considered the most important target with multiple regulations by the most active ingredients. In vitro, we successfully validated the inhibitory role of quercetin in the cellular proliferation of human RA fibroblast-like synoviocyte cell line (MH7A cells). These findings indicated that the potential mechanisms of HLJDD for RA treatment might be attributed to inhibiting the immune-inflammatory response, reducing the release of chemokines, and alleviating the destruction of extracellular matrix (ECM) in the synovial compartment.

The influence of the -94 Ins/Del ATTG polymorphism of NFkB on the anti-CCP antibody levels in patients with rheumatoid arthritis

Rheumatoid arthritis (RA) is an autoimmune disease characterized by an inflammatory process that affects mainly synovial tissue in joints, and by the production of cyclic citrullinated peptides (anti-CCP) antibodies. In the inflammatory process the regulation of the nuclear factor kappa B (NFkB) transcription factor activation is a key point in the production of inflammatory cytokines. On the other hand, polymorphisms in several genes could contribute to the promotion of the inflammatory process observed in RA, and the association of the rs28362491 polymorphism in the NFkB gene with RA has been studied in different population. Therefore, it could be one of the interest targets to analyze their association with RA in a Mexican population.This is a case-control study to determine the influence of rs28362491 in the NFkB gene on RA and on clinical features of this disease, such as anti-CCP antibody levels, Disease Activity Score, and Health Assessment Questionnaire-Disability Index.The genotype of rs28362491 in the NFkB gene was determined in 140 RA patients and 135 healthy controls using the polymerase chain reaction-restriction fragment length polymorphism method with the enzyme PflMI. The following clinical variables were also determined: anti-CCP levels, Disease Activity Score, and Spanish version of the Health Assessment Questionnaire Disability-Index.Although no association of the polymorphism as a risk/protection factor with RA was found, the RA patients who carried the Ins/Ins genotype showed higher anti-CCP levels, while those with the Del/Del genotype showed higher Spanish version of the Health Assessment Questionnaire-Disability Index levels, compared to the other genotypes.The NFkB -94 Ins/Del ATTG (rs28362491) polymorphism is, therefore, associated with higher levels of anti-CCP antibodies, though no significant association as a risk or protection factor in RA cases was identified.

LncRNA ZNF667-AS1 alleviates rheumatoid arthritis by sponging miR-523-3p and inactivating the JAK/STAT signalling pathway

BACKGROUND

Rheumatoid arthritis (RA) is an autoimmune disease, which compromises the synovial membrane resulting in chronic inflammation. Increasing evidence has demonstrated that long non-coding RNAs (lncRNAs) are implicated in the pathogenesis of RA. This study investigated the role of lncRNA ZNF667-AS1 in RA progression.

METHODS

Synovial tissues and fibroblast-like synoviocytes (FLSs) were obtained from patients with RA. Gene expression was measured using RT-qPCR. Chondrocytes were treated with lipopolysaccharide (LPS) to establish in vitro models of OA. Cell counting kit-8 (CCK-8), western blot, and enzyme-linked immunosorbent assay (ELISA) were used to examine the proliferation and inflammatory cytokine production in chondrocytes. Animal models of OA were established in SD rats. Peripheral blood mononuclear cells (PBMCs) were isolated from the OA rats. Flow cytometry was used to measure the changes of the inflammatory T-helper cell 17 (Th17) cells. The relationship between ZNF667-AS1 and miR-523-3p was verified by luciferase reporter assay.

RESULTS

ZNF667-AS1 was downregulated in RA-FLSs and LPS-stimulated chondrocytes. ZNF667-AS1 overexpression significantly promoted cell proliferation and inhibited the production of IL-6, IL-17 and TNF-α in LPS-stimulated chondrocytes. Additionally, ZNF667-AS1 overexpression reduced the generation of CD4 + IL-17+ cells. In mechanism, ZNF667-AS1 acted a sponge for miR-523-3p. MiR-523-3p overexpression reversed the ZNF667-AS1-mediated regulation of cell proliferation and inflammation. Furthermore, miR-523-3p overexpression abolished the inhibitory effects of ZNF667-AS1 on the JAK/STAT signalling activation.

CONCLUSION

ZNF667-AS1 exerts protective effects during RA development by sponging miR-523-3p and inactivating the JAK/STAT signalling.

NFKB1 promoter -94 insertion/deletion ATTG polymorphism (rs28362491) is associated with severity and disease progression of rheumatoid arthritis through interleukin-6 levels modulation in Egyptian patients

OBJECTIVE

Rheumatoid arthritis (RA) is an inflammatory autoimmune disorder, which can cause progressive and functional disability. Previous data suggests that some inflammatory cytokines are dysregulated in patients with RA. Polymorphisms in the NFKB1 gene were studied in different populations with RA. Specific studies showed that the NFKB1 promoter -94ins/delATTG (rs28362491) polymorphism appears to be correlated with alterations in the IL-6 expression and may lead to disease development. We aimed to evaluate the association between the NFKB1 -94ins/delATTG polymorphism and biochemical, and clinical markers for severity of RA in Egyptian patients.

METHODS

Study subjects included 196 RA patients from the Egyptian population. NFKB1 -94ins/delATTG polymorphism was genotyped by real-time PCR using the TaqMan assay. Concentrations of plasma IL-6 were assessed using the ELISA method.

RESULTS

The frequencies of (del/del + ins/del) genotype in cases with erosive arthritis were significantly increased as compared to cases with non-erosive arthritis (63.0% vs. 47.7%, OR = 1.86, 95% CI: 1.05-3.30, p: 0.043). Carriers of del allele had high activity and severity markers compared with those of ins/ins genotype. The del allele was significantly associated with higher IL-6 levels in a dose-dependent manner. Plasma levels of IL-6 were significantly higher in the del/del (41.4 ± 16.2 pg/ml) and ins/del (19.1 ± 12.4 pg/ml) genotype when compared with the ins/ins genotype (11.4 ± 4.21 pg/ml). In a multivariate analysis of variance, including confounding factors associated with higher IL-6 levels (RF, disease duration, and DAS28), the NFKB1 -94ins/delATTG polymorphism retained its role. Logistic regression analyses revealed that high IL-6 plasma levels independently associated with an increased risk of presenting erosive RA, while -94ins/delATTG polymorphism has no direct association with the progression of erosive arthritis.

CONCLUSION

Our data indicate that the NFKB1 -94ins/delATTG polymorphism contributes to the severity and progression of RA through IL-6 levels modulation in Egyptian patients. Key Points • Carriers of del allele had high activity and severity markers compared with those of ins/ins genotype. • In RA patients, the del allele was significantly associated with higher IL-6 levels in a dose-dependent manner. • IL-6 plasma levels are independently associated with an increased risk of presenting erosive arthritis. • The NFKB1 -94ins/delATTG polymorphism contributes to the severity and progression of RA through IL-6 levels modulation in Egyptian patients.