Single-cell RNA-seq reveals altered NK cell subsets and reduced levels of cytotoxic molecules in patients with ankylosing spondylitis.

Ankylosing spondylitis (AS) is an autoimmune disease with unknown aetiology. To unravel the mechanisms mediating AS pathogenesis, we profiled peripheral blood mononuclear cells (PBMCs) from AS patients and healthy subjects using 10X single-cell RNA sequencing. The frequencies of immune cell subsets were evaluated by flow cytometry. NK cells were purified from PBMCs using isolation kit and were examined for gene expression by RT-qPCR. Plasma levels of cytolytic molecules were examined by enzyme-linked immunosorbent assay. Compared to healthy controls, AS patients showed a significant decrease in total NK cells as well as CD56dim NK subset, whereas CD56bright NK cells were increased. Additionally, impaired expression of cytotoxic genes in NK cells of AS patients was observed by bioinformatics algorithm and verified by RT-qPCR and flow cytometry. Consistent with changes in transcriptomics, we found decreased plasma levels of granzymes, but not granulysin, in AS patients. Furthermore, Pearson correlation analysis revealed a negative correlation between plasma GZMB levels and disease activity (r = -0.5275, p = 0.0358). No correlation was observed between plasma cytolytic molecules and biochemical indexes (ESR and CRP). Our findings uncover altered NK cell subsets and cytotoxic profiles in peripheral circulation of AS patients at single-cell resolution.

The efficacy of statins in the treatment of rheumatoid arthritis: A systematic review and meta-analysis.

BACKGROUND: A large body of research has investigated the use of statins in rheumatoid arthritis (RA); however, the therapeutic effects of statins remain uncertain. Thus, we designed a systematic review and meta-analysis to evaluate the role of statins in patients with RA. METHODS: Databases searched to detect clinical randomized controlled trials or clinical controlled trials on the interaction between statins and RA before January 2020 included PubMed, Web of Sciences, Embase, Cochrane Library, CNKI, Wan Fang Database. Efficacy was measured by Disease Activity Score in 28 Joints (DAS28), erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), tenderness of the joint (TJ), swelling of the joint (SJ), and interleukin-6. The level of blood lipid was also evaluated. STATA 12.0 was used for the meta-analysis. The Cochrane method was used for quality assessment. Heterogeneity was considered to determine fixed effects or random effects models. RESULTS: Nineteen clinical trials with a total of 22,906 subjects were included in the meta-analysis. Sixteen studies reported a change in DAS28 after statin treatment. The pooled analysis showed that statins reduced DAS28 in RA patients. Change in ESR after statin treatment was reported in 9 studies. The summary analysis showed that statins lowered ESR in RA patients. Twelve studies reported a change in CRP after statin treatment. The results of the entire analysis showed that statins lowered CRP in RA patients. Seven studies reported a change in TJ after statin treatment. The combined analysis showed that statins reduced TJ at RA patients. Six studies reported changes in IL6 after statin therapy. The results showed that statins failed to reduce IL6 in RA patients. Seven studies reported changes in SJ after statin therapy, which showed that statins failed to reduce SJ in RA patients. We also found that statins can reduce blood lipid levels in RA patients. CONCLUSION: In conclusion, statins were able to reduce DAS28, ESR, CRP, TJ, and blood lipids. It indicated that stains can benefit patients with RA by inhibiting the expression of inflammatory factors and reducing the levels of lipids in the blood. Our study may offer a new perspective on the treatment of RA and provide research ideas for future larger multi-center clinical trials.

WISP2 downregulation inhibits the osteogenic differentiation of BMSCs in congenital scoliosis by regulating Wnt/ß-catenin pathway.

OBJECTIVES: Bone marrow mesenchymal stem cells (BMSCs) are instrumental in bone development, metabolism, and marrow microenvironment homeostasis. Despite this, the relevant effects and mechanisms of BMSCs on congenital scoliosis (CS) remain undefined. Herein, it becomes our focus to reveal the corresponding effects and mechanisms implicated. METHODS: BMSCs from CS patients (hereafter referred as CS-BMSCs) and healthy donors (NC-BMSCs) were observed and identified. Differentially expressed genes in BMSCs were analyzed utilizing scRNA-seq and RNA-seq profiles. The multi-differentiation potential of BMSCs following the transfection or infection was evaluated. The expression levels of factors related to osteogenic differentiation and Wnt/ß-catenin pathway were further determined as appropriate. RESULTS: A decreased osteogenic differentiation ability was shown in CS-BMSCs. Both the proportion of LEPR+ BMSCs and the expression level of WNT1-inducible-signaling pathway protein 2 (WISP2) were decreased in CS-BMSCs. WISP2 knockdown suppressed the osteogenic differentiation of NC-BMSCs, while WISP2 overexpression facilitated the osteogenesis of CS-BMSCs via acting on the Wnt/ß-catenin pathway. CONCLUSIONS: Our study collectively indicates WISP2 knockdown blocks the osteogenic differentiation of BMSCs in CS by regulating Wnt/ß-catenin signaling, thus providing new insights into the aetiology of CS.

Exploring the protective mechanism of baicalin in treatment of atherosclerosis using endothelial cells deregulation model and network pharmacology.

BACKGROUND: Baicalin is a generally available flavonoid with potent biological activity. The present study aimed to assess the underlying mechanism of baicalin in treatment of atherosclerosis (AS) with the help of network pharmacology, molecular docking and experimental validation. METHODS: The target genes of baicalin and AS were identified from public databases, and the overlapping results were considered to be baicalin-AS targets. Core target genes of baicalin were obtained through the PPI network and validated by a clinical microarray dataset (GSE132651). Human aortic endothelial cells (HAECs) were treated with Lipopolysaccharide (LPS) to construct an endothelial injury model. The expression of NOX4 was examined by real-time qPCR and western blot. Flow cytometry was used to detect intracellular levels of reactive oxygen species (ROS). Furthermore, HAECs were transfected with NOX4-specific siRNA and then co-stimulated with baicalin and LPS to investigate whether NOX4 was involved in the anti-oxidative stress effects of baicalin. RESULTS: In this study, baicalin had 45 biological targets against AS. Functional enrichment analysis demonstrated that most targets were involved in oxidative stress. Using the CytoHubba plug-in, we obtained the top 10 genes in the PPI network ranked by the EPC algorithm. Molecular docking and microarray dataset validation indicated that NOX4 may be an essential target of baicalin, and its expression was significantly suppressed in AS samples compared to controls. In endothelial injury model, intervention of HAECs with baicalin increased the expression levels of NOX4 and NOS3 (eNOS), and decreased LPS-induced ROS generation. After inhibition of NOX4, the anti-ROS-generating effect of baicalin was abolished. CONCLUSION: Collectively, we combined network pharmacology and endothelial injury models to investigate the anti-AS mechanism of baicalin. The results demonstrate that baicalin may exert anti-oxidative stress effects by targeting NOX4, providing new mechanisms and insights to baicalin for the treatment of AS.

Identification of the key genes and immune infiltrating cells determined by sex differences in ischaemic stroke through co-expression network module.

Stroke is one of the leading causes of patients' death and long-term disability worldwide, and ischaemic stroke (IS) accounts for nearly 80% of all strokes. Differential genes and weighted gene co-expression network analysis (WGCNA) in male and female patients with IS were compared. The authors used cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT) to analyse the distribution pattern of immune subtypes between male and female patients. In this study, 141 up-regulated and 61 down-regulated genes were gathered and distributed into five modules in response to their correlation degree to clinical traits. The criterion for Gene Ontology (GO) term and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway indicated that detailed analysis had the potential to enhance clinical prediction and to identify the gender-related mechanism. After that, the expression levels of hub genes were measured via the quantitative real-time PCR (qRT-PCR) method. Finally, CCL20, ICAM1 and PTGS2 were identified and these may be some promising targets for sex differences in IS. Besides, the hub genes were further verified by rat experiments. Furthermore, these CIBERSORT results showed that T cells CD8 and Monocytes may be the target for the treatment of male and female patients, respectively.

CD96 Downregulation Promotes the Immune Response of CD4 T Cells and Associates with Ankylosing Spondylitis.

Inhibitory receptors (IRs) play an indispensable role in regulating T cell activation and expansion. This study is aimed at exploring the correlation between IRs and ankylosing spondylitis (AS). Bioinformatics analysis of two datasets (GSE25101 and GSE73754), including 68 AS cases and 36 healthy controls, demonstrated that "T cell receptor signaling pathway" was significantly enriched, and two IRs (CD112R and CD96) were downregulated in AS cases. Real-time Quantitative PCR Detecting System (qPCR) analysis confirmed the decreased expression of CD112R and CD96 in the peripheral blood of AS patients. Flow cytometry demonstrated that the frequency of CD96-positive cells among CD4 T cells in AS patients was significantly reduced and that expressed on the cells was also significantly lower than the healthy controls. In addition, the expression of CD96 was altered on human primary CD4 T cells extracted from 3 healthy volunteers and cocultured with allogeneic dendritic cells (DCs). Also, low expression of CD96 elevated the phosphorylation of ERK in CD4 T cells and increased the level of TNF-α, IL-23, IL-17A, IL-6, and IFN-γ in the cell culture supernatant. These results suggested that CD96 is crucial for the pathogenesis of AS and may be a potential target in the treatment of the disease.

Identification of diagnostic genes and vital microRNAs involved in rheumatoid arthritis: based on data mining and experimental verification.

BACKGROUND: The pathogenesis of rheumatoid arthritis (RA) is complex. This study aimed to identify diagnostic biomarkers and transcriptional regulators that underlie RA based on bioinformatics analysis and experimental verification. MATERIAL AND METHODS: We applied weighted gene co-expression network analysis (WGCNA) to analyze dataset GSE55457 and obtained the key module most relevant to the RA phenotype. We then conducted gene function annotation, gene set enrichment analysis (GSEA) and immunocytes quantitative analysis (CIBERSORT). Moreover, the intersection of differentially expressed genes (DEGs) and genes within the key module were entered into the STRING database to construct an interaction network and to mine hub genes. We predicted microRNA (miRNA) using a web-based tool (miRDB). Finally, hub genes and vital miRNAs were validated with independent GEO datasets, RT-qPCR and Western blot. RESULTS: A total of 367 DEGs were characterized by differential expression analysis. The WGCNA method divided genes into 14 modules, and we focused on the turquoise module containing 845 genes. Gene function annotation and GSEA suggested that immune response and inflammatory signaling pathways are the molecular mechanisms behind RA. Nine hub genes were screened from the network and seven vital regulators were obtained using miRNA prediction. CIBERSORT analysis identified five cell types enriched in RA samples, which were closely related to the expression of hub genes. Through ROC curve and RT-qPCR validation, we confirmed five genes that were specific for RA, including CCL25, CXCL9, CXCL10, CXCL11, and CXCL13. Moreover, we selected a representative gene (CXCL10) for Western blot validation. Vital miRNAs verification showed that only the differences in has-miR-573 and has-miR-34a were statistically significant. CONCLUSION: Our study reveals diagnostic genes and vital microRNAs highly related to RA, which could help improve our understanding of the molecular mechanisms underlying the disorder and provide theoretical support for the future exploration of innovative therapeutic approaches.

Efficacy and Safety of Wenbu Zhibi Granule in Patients with Ankylosing Spondylitis: A Multicenter, Randomized, Double-blind, Placebo-controlled Trial.

Background. Ankylosing spondylitis (AS) is a chronic disease in which the column is the main lesion. It is caused by a combination of genetic and environmental factors, mainly involving the axial skeleton, resulting in column rigidity and difficulty in movement, and there may be different degrees of eye, lung, cardiovascular, kidney, and other organ damage. Long-term treatment lacks in ankylosing spondylitis. Wenbu Zhibi granule (WZG) is a prescription handed down from the history of Chinese medicine for thousands of years, which is used to treat the pain of patients with AS and to prevent the further development of the disease. However, there is no scientific evidence based on clinical trials to evaluate the efficacy and safety of WZG for ankylosing spondylitis. Methods/Design. We will conduct a multicenter, randomized, double-blind, placebo-controlled trial to evaluate the efficacy and safety of the WZG in the treatment of AS. We will randomly assign 100 patients with active AS to two groups, treated for 16 weeks. The primary efficacy endpoint is the proportion of subjects who reached 40% improvement criteria proposed by Assessment of SpondyloArthritis International Society (ASAS40) at 16 weeks from baseline, the secondary efficacy endpoint includes ASAS20 response rate, ASAS partial remission response rate, 5/6 improvement criteria proposed by ASAS (ASAS5/6) response rate, and change in the Spondyloarthritis Research Consortium of Canada (SPARCC) MRI spine score, Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), Bath Ankylosing Spondylitis Functional Index (BASFI), Ankylosing Spondylitis Disease Activity Score (ASDAS), linear Bath Ankylosing Spondylitis Metrology Index (BASMI), ankylosing spondylitis quality of life (ASQoL). In addition, the time points will be set as baseline, 2 weeks, 4 weeks, 8 weeks, 12 weeks, 16 weeks, 24 weeks, and 48 weeks. Discussion. The results of this study will elucidate the efficacy and safety of WZG and provide an appropriate treatment option for patients with AS. Trial registration: ClinicalTrials.gov ID: https://clinicaltrials.gov/ct2/show/ChiCTR2000041010. (Chinese Clinical Trail Registry, Registered 16 December 2020, http://www.chictr.org.cn).

Identification of immune related cells and crucial genes in the peripheral blood of ankylosing spondylitis by integrated bioinformatics analysis.

BACKGROUND: Ankylosing spondylitis (AS) is a progressive rheumatic disease and studies reveal that the immune system is critical for the pathogenesis of AS. In the present study, various bioinformatics analysis methods were comprehensively applied, designed to identify potential key genes and inflammation states of AS. METHODS: The transcriptome profiles of GSE25101 and GSE73754 obtained from the Gene Expression Omnibus (GEO) database were merged for subsequent analyses. The differentially expressed genes (DEGs) were identified using the Bioconductor package Limma and threshold values. Functional enrichment and pathway enrichment analyses were performed using the clusterProfiler package and Gene Set Enrichment Analysis (GSEA). Next, protein-protein interaction (PPI) network of the identified DEGs was constructed by the online database, the Search Tool for the Retrieval of Interacting Genes (STRING), visualization and analysis were performed through Cytoscape software. Subsequently, we applied CIBERSORT algorithm to identify subpopulation proportions of immune cells in peripheral blood samples. Finally, we validated the hub genes with the GSE18781 dataset. Samples were collected from patients to validate gene and protein expression using qRT-PCR and ELISA. RESULTS: A total of 334 DEGs were identified, including 182 upregulated and 152 downregulated DEGs, between AS patients and normal human controls, which were primarily involved in immune response, autophagy, and natural killer cell-mediated cytotoxicity. The most prominent module and candidate biomarkers were identified from the PPI network. Biomarkers were selected for validation and their expressions were significantly decreased in peripheral blood samples which was consistent with transcriptome sequencing results. Nine genes with AUC > 0.70 were considered to be AS hub genes for ROC curve analysis, including GZMA, GZMK, PRF1, GNLY, NKG7, KLRB1, KLRD1, IL2RB and CD247. Furthermore, CIBERSORT results suggest that AS contained a higher proportion of CD8+ T cells, naive CD4+ T cells, neutrophils, and lower levels of gamma delta T cells compared with the normal controls. CONCLUSION: In this study, we identified DEGs combined with their closely related biological functions and propose that granule-associated proteins and immune infiltration maybe involved in the progression of ankylosing spondylitis. These validated hub genes may provide new perspectives for understanding the molecular mechanisms of ankylosing spondylitis.

Comprehensive Bioinformatics Analysis Reveals Hub Genes and Inflammation State of Rheumatoid Arthritis.

Rheumatoid arthritis (RA) is an autoimmune disease characterized by erosive arthritis, which has not been thoroughly cured yet, and standardized treatment is helpful for alleviating clinical symptoms. Here, various bioinformatics analysis tools were comprehensively utilized, aiming to identify critical biomarkers and possible pathogenesis of RA. Three gene expression datasets profiled by microarray were obtained from GEO database. Dataset GSE55235 and GSE55457 were merged for subsequent analyses. We identified differentially expressed genes (DEGs) in RStudio with limma package, performing functional enrichment analysis based on GSEA software and clusterProfiler package. Next, protein-protein interaction (PPI) network was set up through STRING database and Cytoscape. Moreover, CIBERSORT website was used to assess the inflammatory state of RA. Finally, we validated the candidate hub genes with dataset GSE77298. As a result, we identified 106 DEGs (72 upregulated and 34 downregulated genes). Through GO, KEGG, and GSEA analysis, we found that DEGs were mainly involved in immune response and inflammatory signaling pathway. With the help of Cytoscape software and MCODE plug-in, the most prominent subnetwork was screened out, containing 14 genes and 45 edges. For ROC curve analysis, eight genes with AUC >0.80 were considered as hub genes of RA. In conclusion, compared with healthy controls, the DEGs and their closely related biological functions were analyzed, and we held that chemokines and immune cells infiltration promote the progression of rheumatoid arthritis. Targeting the eight biomarkers we identified may be useful for the diagnosis and treatment of rheumatoid arthritis.

Bioinformatics Analysis Reveals Diagnostic Markers and Vital Pathways Involved in Acute Coronary Syndrome.

BACKGROUND: Acute coronary syndrome (ACS) has a high incidence and mortality rate. Early detection and intervention would provide clinical benefits. This study aimed to reveal hub genes, transcription factors (TFs), and microRNAs (miRNAs) that affect plaque stability and provide the possibility for the early diagnosis and treatment of ACS. METHODS: We obtained gene expression matrix GSE19339 for ACS patients and healthy subjects from public database. The differentially expressed genes (DEGs) were screened using Limma package in R software. The biological functions of DEGs were shown by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA). Protein-protein interaction (PPI) network was mapped in Cytoscape, followed by screening of hub genes based on the Molecular Complex Detection (MCODE) plug-in. Functional Enrichment analysis tool (FunRich) and Database for Annotation, Visualization and Integrated Discovery (DAVID) were used to predict miRNAs and TFs, respectively. Finally, GSE60993 expression matrix was chosen to plot receiver operating characteristic (ROC) curves with the aim of further assessing the reliability of our findings. RESULTS: We obtained 176 DEGs and further identified 16 hub genes by MCODE. The results of functional enrichment analysis showed that DEGs mediated inflammatory response and immune-related pathways. Among the predicted miRNAs, hsa-miR-4770, hsa-miR-5195, and hsa-miR-6088 all possessed two target genes, which might be closely related to the development of ACS. Moreover, we identified 11 TFs regulating hub gene transcriptional processes. Finally, ROC curves confirmed three genes with high confidence (area under the curve > 0.9), including VEGFA, SPP1, and VCAM1. CONCLUSION: This study suggests that three genes (VEGFA, SPP1, and VCAM1) were involved in the molecular mechanisms of ACS pathogenesis and could serve as biomarkers of disease progression.

Acupuncture Therapy plus Hyaluronic Acid Injection for Knee Osteoarthritis: A Meta-Analysis of Randomized Controlled Trials.

OBJECTIVE: This meta-analysis aimed to investigate the effectiveness of acupuncture therapy plus hyaluronic acid injection versus hyaluronic acid injection alone for patients with knee osteoarthritis. METHODS: Relevant randomized controlled trials that compared the combined effect of acupuncture therapy and hyaluronic acid injection with hyaluronic acid injection alone for knee osteoarthritis patients were included. 10 studies were included in this meta-analysis, and the relative risk (RR) and weight mean difference (MD) with 95% CI for the Lysholm knee score (LKSS), visual analogue scale (VAS), and effective rate (ER) were evaluated by using RevMan 5.3 software. Besides, the bias assessment of the included studies was evaluated using the Cochrane risk of bias tool, and the GRADE (Grading of Recommendations, Assessment Development, and Evaluation) system was applied to assess the overall quality of the evidence. RESULTS: A total of 10 studies involving 998 participants were included in this study. Compared to hyaluronic acid injection alone, the combined therapy significantly reduced pain on the visual analogue scale (VAS) and improved the ER and knee function on the Lysholm knee score (LKSS). Of these, the pooled LKSS (MD = 8.09, 95% CI = [7.02, 9.16], p < 0.00001, 7 studies) and ER (RR = 1.23, 95% CI 1.15 to 1.31, p < 0.00001, 7 studies) and ER (RR = 1.23, 95% CI 1.15 to 1.31, p < 0.00001, 7 studies) and ER (RR = 1.23, 95% CI 1.15 to 1.31. CONCLUSION: Current evidence suggests that acupuncture therapy combined with hyaluronic acid injection is more effective in alleviating pain, improving the ER and knee function compared with hyaluronic acid injection alone. However, considering the low quality, small size, and high risk of the studies identified in this meta-analysis, more higher methodological quality, rigorously designed randomized controlled trials with large sample sizes are needed to confirm the results.