Multiple datasets to explore the tumor microenvironment of cutaneous squamous cell carcinoma

Machine learning and weighted gene co-expression network analysis identify a three-gene signature to diagnose rheumatoid arthritis

Background

Rheumatoid arthritis (RA) is a systemic immune-related disease characterized by synovial inflammation and destruction of joint cartilage. The pathogenesis of RA remains unclear, and diagnostic markers with high sensitivity and specificity are needed urgently. This study aims to identify potential biomarkers in the synovium for diagnosing RA and to investigate their association with immune infiltration.

Methods

We downloaded four datasets containing 51 RA and 36 healthy synovium samples from the Gene Expression Omnibus database. Differentially expressed genes were identified using R. Then, various enrichment analyses were conducted. Subsequently, weighted gene co-expression network analysis (WGCNA), random forest (RF), support vector machine-recursive feature elimination (SVM-RFE), and least absolute shrinkage and selection operator (LASSO) were used to identify the hub genes for RA diagnosis. Receiver operating characteristic curves and nomogram models were used to validate the specificity and sensitivity of hub genes. Additionally, we analyzed the infiltration levels of 28 immune cells in the expression profile and their relationship with the hub genes using single-sample gene set enrichment analysis.

Results

Three hub genes, namely, ribonucleotide reductase regulatory subunit M2 (RRM2), DLG-associated protein 5 (DLGAP5), and kinesin family member 11 (KIF11), were identified through WGCNA, LASSO, SVM-RFE, and RF algorithms. These hub genes correlated strongly with T cells, natural killer cells, and macrophage cells as indicated by immune cell infiltration analysis.

Conclusion

RRM2, DLGAP5, and KIF11 could serve as potential diagnostic indicators and treatment targets for RA. The infiltration of immune cells offers additional insights into the underlying mechanisms involved in the progression of RA.

Gene signature based on glycolysis is closely related to immune infiltration of patients with osteoarthritis

BACKGROUND

Osteoarthritis (OA) is a degenerative arthritis with high levels of clinical heterogeneity. Aberrant metabolism such as shifting from oxidative phosphorylation to glycolysis is a response to changes in the inflammatory microenvironment of OA. Therefore, there is a pressing need to identify novel glycolysis regulators during OA progression.

METHODS

We systematically studied glycolysis patterns mediated by 141 glycolysis regulators in 74 human synovial samples and discussed the characteristics of the immune microenvironment modified by glycolysis. The random forest (RF) method was applied to screen candidate hub glycolysis regulators in OA. RT-qPCR was performed to validate these key regulators. Then distinct glycolysis patterns were identified, and systematic correlation between these glycolysis patterns and immune cell infiltration was analyzed. The glycolysis score was constructed to quantify glycolysis patterns together with immune infiltration of individual OA patient.

RESULTS

56 glycolysis-related differentially expressed genes (DEGs) were identified between OA and non-OA samples. STC1, VEGFA, KDELR3, DDIT4 and PGAM1 were selected as candidate genes to predict the probability of OA. Two glycolysis patterns in OA were identified. Glycolysis cluster A with higher glycolysis score was related to an inflamed phenotype.

CONCLUSIONS

Taken together, our results established a glycolysis-based genetic signature for OA, guided in-depth studies on the metabolic mechanism of OA, and facilitated to explore new clinical treatment strategies.

Identification of four-gene signature to diagnose osteoarthritis through bioinformatics and machine learning methods

BACKGROUND

Although osteoarthritis (OA) is one of the most prevalent joint disorders, effective biomarkers to diagnose OA are still unavailable. This study aimed to acquire some key synovial biomarkers (hub genes) and analyze their correlation with immune infiltration in OA.

METHODS

Gene expression profiles and clinical characteristics of OA and healthy synovial samples were retrieved from the Gene Expression Omnibus (GEO) database. Hub genes for OA were mined based on a combination of weighted gene co-expression network analysis (WGCNA), the least absolute shrinkage and selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE), and random forest (RF) algorithms. A diagnostic nomogram model for OA prediction was developed based on the hub genes. Receiver operating characteristic curves (ROC) were performed to confirm the abnormal expression of hub genes in the experimemtal and validation datasets. qRT-PCR using patients' samples were conducted as well. In addition, the infiltration level of 28 immune cells in the expression profile and their relationship with hub genes were analyzed using single-sample GSEA (ssGSEA).

RESULTS

4 hub genes (ZBTB16, TNFSF11, SCRG1 and KDELR3) were obtained by WGCNA, lasso, SVM-RFE, RF algorithms as potential biomarkers for OA. The immune infiltration analyses revealed that hub genes were most correlated with regulatory T cell and natural killer cell.

CONCLUSION

A machine learning model to diagnose OA based on ZBTB16, TNFSF11, SCRG1 and KDELR3 using synovial tissue was constructed, providing theoretical foundation and guideline for diagnostic and treatment targets in OA.

Identification and validation of BCL6 and VEGFA as biomarkers and ageing patterns correlating with immune infiltrates in OA progression

Osteoarthritis (OA), the most common type of arthritis, is a complex biological response caused by cartilage wear and synovial inflammation that links biomechanics and inflammation. The progression of OA correlates with a rise in the number of senescent cells in multiple joint tissues. However, the mechanisms by which senescent cells and their involvement with immune infiltration promote OA progression are not fully understood. The gene expression profiles and clinical information of OA and healthy control synovial tissue samples were retrieved from the Gene Expression Omnibus database, and then differential analysis of senescence regulators between OA and normal samples was performed. The random forest (RF) was used to screen candidate senescence regulators to predict the occurrence of OA. The reverse transcription quantitative real-time PCR experiments at tissue's level was performed to confirm these biomarkers. Moreover, two distinct senescence patterns were identified and systematic correlation between these senescence patterns and immune cell infiltration was analyzed. The senescence score and senescence gene clusters were constructed to quantify senescence patterns together with immune infiltration of individual OA patient. 73 senescence differentially expressed genes were identified between OA patients and normal controls. The RF method was utilized to build an OA risk model based on two senescence related genes: BCL6 and VEGFA. Next, two distinct aging patterns were determined in OA synovial samples. Most patients from senescence cluster A were further classified into gene cluster B and high senescence score group correlated with a non-inflamed phenotype, whereas senescence cluster B were classified into gene cluster A and low senescence score group correlated with an inflamed phenotype. Our study revealed that senescence played an important role in in OA synovial inflammation. Evaluating the senescence patterns of individuals with OA will contribute to enhancing our cognition of immune infiltration characterization, providing novel diagnostic and prognostic biomarkers, and guiding more effective immunotherapy strategies.