Comprehensive Analysis of Hub Telomere-Related Genes and Synovial Immune Characteristics in Osteoarthritis.

Objectives: Idiopathic pulmonary arterial hypertension (IPAH) is a rare but severe lung disorder, which may lead to heart failure and early mortality. However, little is known about the etiology of IPAH. Thus, the present study aimed to establish the differentially expressed genes (DEGs) between IPAH and normal tissues, which may serve as potential prognostic markers in IPAH. Furthermore, we utilized a versatile computational method, CIBERSORT to identify immune cell infiltration characteristics in IPAH.Materials and Methods: The GSE117261 and GSE48149 datasets were obtained from the Gene Expression Omnibus database. The GSE117261 dataset was adopted to screen DEGs between IPAH and the control groups with the criterion of |log2 fold change| ≥ 1, adjusted P < 0.05, and to further explore their potential biological functions via Gene Ontology analysis, Kyoto Encyclopedia of Genes and Genomes Pathway analysis, and Gene Set Enrichment Analysis. Moreover, the support vector machine (SVM)-recursive feature elimination and the least absolute shrinkage and selection operator regression model were performed jointly to identify the best potential biomarkers. Then we built a regression model based on these selected variables. The GSE48149 dataset was used as a validation cohort to appraise the diagnostic efficacy of the SVM classifier by receiver operating characteristic (ROC) analysis. Finally, immune infiltration was explored by CIBERSORT in IPAH. We further analyzed the correlation between potential biomarkers and immune cells.Results: In total, 75 DEGs were identified; 40 were downregulated, and 35 genes were upregulated. Functional enrichment analysis found a significantly enrichment in heme binding, inflammation, chemokines, cytokine activity, and abnormal glycometabolism. HBB, RNASE2, S100A9, and IL1R2 were identified as the best potential biomarkers with an area under the ROC curve (AUC) of 1 (95%CI = 0.937–1.000, specificity = 100%, sensitivity = 100%) in the discovery cohort and 1(95%CI = 0.805–1.000, specificity = 100%, sensitivity = 100%) in the validation cohort. Moreover, immune infiltration analysis by CIBERSORT showed a higher level of CD8+ T cells, resting memory CD4+ T cells, gamma delta T cells, M1 macrophages, resting mast cells, as well as a lower level of naïve CD4+ T cells, monocytes, M0 macrophages, activated mast cells, and neutrophils in IPAH compared with the control group. In addition, HBB, RNASE2, S100A9, and IL1R2 were correlated with immune cells.Conclusion: HBB, RNASE2, S100A9, and IL1R2 were identified as potential biomarkers to discriminate IPAH from the control. There was an obvious difference in immune infiltration between patient with IPAH and normal groups.

This study aimed to find potential diagnostic markers for osteoarthritis (OA) and analyze the role of immune cells infiltration in this pathology. We used OA datasets from the Gene Expression Omnibus database. First, R software was used to identify differentially expressed genes (DEGs) and perform functional correlation analysis. Then least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine-recursive feature elimination algorithms were used to screen and verify the diagnostic markers of OA. Finally, CIBERSORT was used to evaluate the infiltration of immune cells in OA tissues, and the correlation between diagnostic markers and infiltrating immune cells was analyzed. A total of 458 DEGs were screened in this study. GRB10 and E2F3 (AUC = 0.962) were identified as diagnostic markers of OA. Immune cell infiltration analysis found that resting mast cells, T regulatory cells, CD4 memory resting T cells, activated NK cells, and eosinophils may be involved in the OA process. In addition, GRB10 was correlated with NK resting cells, naive CD4 + T cells, and M1 macrophages, while E2F3 was correlated with resting mast cells. In conclusion, GRB10 and E2F3 can be used as diagnostic markers of osteoarthritis, and immune cell infiltration plays an important role in the occurrence and progression of OA.

Osteoarthritis (OA) is one of the most common refractory degenerative joint diseases worldwide. Synovitis is believed to drive joint cartilage destruction during OA pathogenesis. Cuproptosis is a novel form of copper-induced cell death. However, few studies have examined the correlations between cuproptosis-related genes (CRGs), immune infiltration, and synovitis. Therefore, we analyzed CRGs in synovitis during OA. Microarray datasets (GSE55235, GSE55457, GSE12021, GSE82107 and GSE176308) were downloaded from the Gene Expression Omnibus database. Next, we conducted differential and subtype analyses of CRGs across synovitis. Immune infiltration and correlation analyses were performed to explore the association between CRGs and immune cell abundance in synovitis. Finally, single-cell RNA-seq profiling was performed using the GSE176308 dataset to investigate the expression of CRGs in the various cell clusters. We found that the expression of five CRGs (FDX1, LIPT1, PDHA1, PDHB, and CDKN2A) was significantly increased in the OA synovium. Moreover, abundant and various types of immune cells infiltrated the synovium during OA, which was correlated with the expression of CRGs. Additionally, single-cell RNA-seq profiling revealed that the cellular composition of the synovium was complex and that their proportions varied greatly as OA progressed. The expression of CRGs differed across various cell types in the OA synovium. The current study predicted that cuproptosis may be involved in the pathogenesis of synovitis. The five screened CRGs (FDX1, LIPT1, PDHA1, PDHB, and CDKN2A) could be explored as candidate biomarkers or therapeutic targets for OA synovitis.

Gliomas are the most prevalent primary brain tumors, and patients typically exhibit poor prognoses. Increasing evidence suggests that telomere maintenance mechanisms play a crucial role in glioma development. However, the prognostic value of telomere-related genes in glioma remains uncertain. This study aimed to construct a prognostic model of telomere-related genes and further elucidate the potential association between the two. We acquired RNA-seq data for low-grade glioma (LGG) and glioblastoma (GBM), along with corresponding clinical information from The Cancer Genome Atlas (TCGA) database, and normal brain tissue data from the Genotype-Tissue Expression (GTEX) database for differential analysis. Telomere-related genes were obtained from TelNet. Initially, we conducted a differential analysis on TCGA and GTEX data to identify differentially expressed telomere-related genes, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses on these genes. Subsequently, univariate Cox analysis and log-rank tests were employed to obtain prognosis-related genes. Least absolute shrinkage and selection operator (LASSO) regression analysis and multivariate Cox regression analysis were sequentially utilized to construct prognostic models. The model's robustness was demonstrated using receiver operating characteristic (ROC) curve analysis, and multivariate Cox regression of risk scores for clinical characteristics and prognostic models were calculated to assess independent prognostic factors. The aforementioned results were validated using the Chinese Glioma Genome Atlas (CGGA) dataset. Finally, the CIBERSORT algorithm analyzed differences in immune cell infiltration levels between high- and low-risk groups, and candidate genes were validated in the Human Protein Atlas (HPA) database. Differential analysis yielded 496 differentially expressed telomere-related genes. GO and KEGG pathway analyses indicated that these genes were primarily involved in telomere-related biological processes and pathways. Subsequently, a prognostic model comprising ten telomere-related genes was constructed through univariate Cox regression analysis, log-rank test, LASSO regression analysis, and multivariate Cox regression analysis. Patients were stratified into high-risk and low-risk groups based on risk scores. Kaplan-Meier (K-M) survival analysis revealed worse outcomes in the high-risk group compared to the low-risk group, and establishing that this prognostic model was a significant independent prognostic factor for glioma patients. Lastly, immune infiltration analysis was conducted, uncovering notable differences in the proportion of multiple immune cell infiltrations between high- and low-risk groups, and eight candidate genes were verified in the HPA database. This study successfully constructed a prognostic model of telomere-related genes, which can more accurately predict glioma patient prognosis, offer potential targets and a theoretical basis for glioma treatment, and serve as a reference for immunotherapy through immune infiltration analysis.

BackgroundAutophagy plays a critical role in the progression of osteoarthritis (OA), mainly by regulating inflammatory and immune responses. However, the underlying mechanisms remain unclear. This study aimed to investigate the potential relevance of autophagy-related genes (ARGs) associated with infiltrating immune cells in OA.MethodsGSE114007, GSE169077, and ARGs were obtained from the Gene Expression Omnibus (GEO) database and the Human Autophagy database. R software was used to identify the differentially expressed autophagy-related genes (DEARGs) in OA. Functional enrichment and protein–protein interaction (PPI) analyses were performed to explore the role of DEARGs in OA cartilage, and then Cytoscape was utilized to screen hub ARGs. Single-sample gene set enrichment analysis (ssGSEA) was used to conduct immune infiltration analysis and evaluate the potential correlation of key ARGs and immune cell infiltration. Then, the expression levels of hub ARGs in OA were further verified by the GSE169077 and qRT-PCR. Finally, Western blotting and immunohistochemistry were used to validate the final hub ARGs.ResultsA total of 24 downregulated genes and five upregulated genes were identified, and these genes were enriched in autophagy, mitophagy, and inflammation-related pathways. The intersection results identified nine hub genes, namely, CDKN1A, DDIT3, FOS, VEGFA, RELA, MAP1LC3B, MYC, HSPA5, and HSPA8. GSE169077 and qRT-PCR validation results showed that only four genes, CDKN1A, DDT3, MAP1LC3B, and MYC, were consistent with the bioinformatics analysis results. Western blotting and immunohistochemical (IHC) showed that the expression of these four genes was significantly downregulated in the OA group, which is consistent with the qPCR results. Immune infiltration correlation analysis indicated that DDIT3 was negatively correlated with immature dendritic cells in OA, and FOS was positively correlated with eosinophils.ConclusionCDKN1A, DDIT3, MAP1LC3B, and MYC were identified as ARGs that were closely associated with immune infiltration in OA cartilage. Among them, DDIT3 showed a strong negative correlation with immature dendritic cells. This study found that the interaction between ARGs and immune cell infiltration may play a crucial role in the pathogenesis of OA; however, the specific interaction mechanism needs further research to be clarified. This study provides new insights to further understand the molecular mechanisms of immunity involved in the process of OA by autophagy.

This study aimed to identify hub genes associated with mitophagy involved in the pathogenesis and progression of diabetic foot ulcer (DFU), and to characterize their immune cell infiltration features and single-cell expression profiles. DFU-related datasets (GSE80178, GSE68183) were retrieved from the GEO database. Subsequently, differentially expressed genes (DEGs) were identified via limma analysis, followed by gene set enrichment analysis (GSEA) to assess gene function enrichment. Identified DEGs were intersected with mitophagy-related genes. Machine learning (ML) algorithms were further employed to identify hub genes. Additionally, immune cell infiltration was examined via the CIBERSORT algorithm, and the correlation between the identified genes and immune infiltration was investigated. Finally, hub genes identified were validated via the single-cell RNA sequencing dataset GSE165816, and further validated using RT-PCR and Western blot (WB) assays. Two hub genes, ANO6 and ALDH2, were identified and found to be significantly downregulated in the skin tissues of patients with DFU. Receiver operating characteristic (ROC) analysis demonstrated robust diagnostic potential (ANO6, AUC = 0.833, ALDH2, AUC = 0.806). Immune cell infiltration analysis demonstrated notable differences between the DFU and normal groups in naïve B cells, monocytes, resting mast cells, γδT cells, and regulatory T cells (Tregs). The findings were further validated through single-cell RNA sequencing (scRNA-seq) analysis and experimental studies, which confirmed the downregulation of ANO6 and ALDH2 in DFU tissues. Two mitophagy-related hub genes, ANO6 and ALDH2, were identified and validated as being significantly downregulated in DFU. Both genes demonstrated diagnostic potential and showed an association with immune cell infiltration. These findings suggest that mitophagy dysfunction may contribute to the pathophysiology of DFU, potentially through the dysregulation of inflammatory pathways and immune responses. While the results provide valuable insights into DFU and its management, further studies with larger cohorts and deeper exploration of mechanistic links to inflammation are necessary to translate these findings into therapeutic strategies.

The poliovirus receptor (PVR) is a member of the immunoglobulin superfamily (Ig SF) and is essential for the promotion of cancer cell proliferation and invasion. However, the correlation between PVR expression and prognosis as well as immune infiltration in hepatocellular carcinoma (HCC) remains unclear. The expression level of PVR was quantified using the Tumor and Tumor Immunity Evaluation Resource (TIMER) and Sangerbox. The Gene Expression Omnibus (GEO) database was used to validate the PVR expression. The receiver operating characteristic (ROC) curve was used to evaluate the feasibility of using PVR as a differentiating factor according to the area under curve (AUC) score. A PVR binding protein network was built using the STRING tool. An enrichment analysis using the R package clusterProfiler was used to explore the potential function of PVR. Immune infiltration analysis was calculated with ESTIMATE algorithms. We also assessed the correlation between PVR expression and immune infiltration by the single-sample Gene Set Enrichment Analysis (ssGSEA) method from the R package GSVA and TIMER database. The results showed that PVR was commonly overexpressed in multiple types of tumors including HCC. The data of GSE64041 confirmed the same result. The ROC curve suggested that PVR could be a potential diagnostic biomarker. Additionally, high mRNA expression of PVR in HCC was significantly correlated with poor overall survival (OS) and relapse free survival (RFS). Results also indicated correlations between PVR mRNA expression with the level of infiltration immune cells including B cells, CD8+ T cells, cytotoxic cells, DCs, CD56dim NK cells, pDCs, and Th2 cells. Furthermore, the PVR level was significantly correlated with immune markers for immunosuppressive cells in HCC. In conclusion, PVR might be an important regulator of tumor immune cell infiltration and a valuable prognostic biomarker in HCC. However, additional work is needed to fully elucidate the underlying mechanisms.

Osteoarthritis (OA) and rheumatoid arthritis (RA) are two common diseases that result in limb disability and a decrease in quality of life. The major symptoms of OA and RA are pain, swelling, stiffness, and malformation of joints, and each disease also has unique characteristics. To compare the pathological mechanisms of OA and RA via weighted correlation network analysis (WGCNA) and immune infiltration analysis and find potential diagnostic and pharmaceutical targets for the treatment of OA and RA. The gene expression profiles of ten OA and ten RA synovial tissue samples were downloaded from the Gene Expression Omnibus (GEO) database (GSE55235). After obtaining differentially expressed genes (DEGs) via GEO2R, WGCNA was conducted using an R package, and modules and genes that were highly correlated with OA and RA were identified. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and protein-protein interaction (PPI) network analyses were also conducted. Hub genes were identified using the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape software. Immune infiltration analysis was conducted using the Perl program and CIBERSORT software. Two hundred ninety-nine DEGs, 24 modules, 16 GO enrichment terms, 6 KEGG pathway enrichment terms, 10 hub genes (CXCL9, CXCL10, CXCR4, CD27, CD69, CD3D, IL7R, STAT1, RGS1, and ISG20), and 8 kinds of different infiltrating immune cells (plasma cells, CD8 T cells, activated memory CD4 T cells, T helper follicular cells, M1 macrophages, Tregs, resting mast cells, and neutrophils) were found to be involved in the different pathological mechanisms of OA and RA. Inflammation-associated genes were the top differentially expressed hub genes between OA and RA, and their expression was downregulated in OA. Genes associated with lipid metabolism may have upregulated expression in OA. In addition, immune cells that participate in the adaptive immune response play an important role in RA. OA mainly involves immune cells that are associated with the innate immune response.

Objectives: This study aimed to identify novel biomarkers for osteoarthritis (OA) and explore potential pathological immune cell infiltration. Methods: We identified differentially expressed genes (DEGs) between OA and normal synovial tissues using the limma package in R, and performed enrichment analyses to understand the functions and enriched pathways of DEGs. Weighted gene co-expression network analysis (WGCNA) and distinct machine-learning algorithms were then used to identify hub modules and candidate biomarkers. We assessed the diagnostic value of the candidate biomarkers using receiver operating characteristic (ROC) analysis. We then used the CIBERSORT algorithm to analyze immune cell infiltration patterns, and the Wilcoxon test to screen out hub immune cells that might affect OA occurrence. Finally, the expression levels of hub biomarkers were confirmed by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Results: We identified 102 up-regulated genes and 110 down-regulated genes. The functional enrichment analysis results showed that DEGs are enriched mainly in immune response pathways. Combining the results of the algorithms and ROC analysis, we identified GUCA1A and NELL1 as potential diagnostic biomarkers for OA, and validated their diagnosibility using an external dataset. Construction of a TF-mRNA-miRNA network enabled prediction of potential candidate compounds targeting hub biomarkers. Immune cell infiltration analyses revealed the expression of hub biomarkers to be correlated with CD8 T cells, memory B cells, M0/M2 macrophages, resting mast cells and resting dendritic cells. qRT-PCR results showed both GUCA1A and NELL1 were significantly increased in OA samples (p < 0.01). All validations are consistent with the microarray hybridization, indicating that GUCA1A and NELL1 may be involved in the pathogenesis of OA. Conclusion: The findings suggest that GUCA1A and NELL1, closely related to OA occurrence and progression, represent new OA candidate markers, and that immune cell infiltration plays a significant role in the progression of OA.

Chemotaxis of cultured human mesenchymal stem cells to osteoarthritic cartilage and synovium

Immune infiltration plays a pivotal role in the pathogenesis of ischemic stroke. A novel form of cell death known as disulfidptosis has emerged in recent studies. However, there is currently a lack of research investigating the regulatory mechanism of disulfidptosis-related genes in immune infiltration during ischemic stroke. Using machine learning methods, we identified candidate key disulfidptosis-related genes (DRGs). Subsequently, we performed an analysis of immune cell infiltration to investigate the dysregulation of immune cells in the context of ischemic stroke. We assessed their diagnostic value by employing receiver operating characteristic (ROC) curves. To gain further insights, we conducted functional enrichment analyses to elucidate the signaling pathways associated with these seven DRGs. We identified two distinct subclusters based on the expression patterns of these seven DRGs. The unique roles of these subclusters were further evaluated through KEGG analysis and immune infiltration studies. Furthermore, we validated the expression profiles of these seven DRGs using both single-cell datasets and external datasets. Lastly, molecular docking was performed to explore potential drugs for the treatment of ischemic stroke. We identified seven DRGs. The seven DRGs are related to immune cells. Additionally, these seven DRGs also demonstrate potential diagnostic value in ischemic stroke. Functional enrichment analysis highlighted pathways such as platelet aggregation and platelet activation. Two subclusters related to disulfidptosis were defined, and functional enrichment analysis of their differentially expressed genes (DEGs) primarily involved pathways like cytokine-cytokine receptor interaction. Single-cell analysis indicated that these seven DRGs were primarily distributed among immune cell types. Molecular docking results suggested that genistein might be a potential therapeutic drug. This study has opened up new avenues for exploring the causes of ischemic stroke and developing potential therapeutic targets.

To build an innovative telomere-associated scoring model to predict prognosis and treatment responsiveness in acute myeloid leukemia (AML). AML is a highly heterogeneous malignant hematologic disorder with a poor prognosis. While telomere maintenance is frequently observed in tumors, investigations into telomere-related genes (TRGs) in AML remain limited. This study aimed to identify prognostic TRGs using the least absolute shrinkage and selection operator (LASSO) Cox regression and multivariate Cox regression, evaluate their predictive value, explore the association between TRG scores and immune cell infiltration, and assess the sensitivity of high-scoring AML patients to chemotherapeutic agents. Univariate Cox regression analysis was conducted on the TCGA cohort to identify prognostic TRGs and to develop the TRG scoring model using LASSO-Cox and multivariate Cox regression. Validation was performed on the GSE37642 cohort. Immune cell infiltration patterns were assessed through computational analysis, and the sensitivity to chemotherapeutic agents was evaluated. Thirteen prognostic TRGs were identified, and a seven-TRG scoring model (including NOP10, OBFC1, PINX1, RPA2, SMG5, MAPKAPK5, and SMN1) was developed. Higher TRG scores were associated with a poorer prognosis, as confirmed in the GSE37642 cohort, and remained an independent prognostic factor even after adjusting for other clinical characteristics. The high-score group was characterized by elevated infiltration of B cells, T helper cells, natural killer cells, tumor-infiltrating lymphocytes, regulatory T (Treg) cells, M2 macrophages, neutrophils, and monocytes, along with reduced infiltration of gamma delta T cells, CD4- T cells, and resting mast cells. Moreover, high infiltration of M2 macrophages and Tregs was associated with poor overall survival compared to low infiltration. Notably, high-risk AML patients were resistant to Erlotinib, Parthenolide, and Nutlin-3a, but sensitive to AC220, Midostaurin, and Tipifarnib. Additionally, using RT-qPCR, we observed significantly higher expression of two model genes, OBFC1 and SMN1, in AML tissues compared to control tissues. This innovative TRG scoring model demonstrates considerable predictive value for AML patient prognosis, offering valuable insights for optimizing treatment strategies and personalized medicine approaches. The identified TRGs and associated scoring models could aid in risk stratification and guide tailored therapeutic interventions in AML patients.

BackgroundThe mechanisms whereby CYFIP2 acts in tumor development and drives immune infiltration have been poorly explored. Thus, this study aimed to identifying the role of CYFIP2 in tumors and immune response.MethodsIn this study, we first explored expression patterns, diagnostic role and prognostic value of CYFIP2 in cancers, particularly in lung adenocarcinoma (LUAD). Then, we performed functional enrichment, genetic alterations, DNA methylation analysis, and immune cell infiltration analysis of CYFIP2 to uncover its potential mechanisms involved in immune microenvironment.ResultsWe found that CYFIP2 significantly differentially expressed in different tumors including LUAD compared with normal tissues. Furthermore, CYFIP2 was found to be significantly correlated with clinical parameters in LUAD. According to the diagnostic and survival analysis, CYFIP2 may be employed as a potential diagnostic and prognostic biomarker. Moreover, genetic alterations revealed that mutation of CYFIP2 was the main types of alterations in different cancers. DNA methylation analysis indicated that CYFIP2 mRNA expression correlated with hypomethylation. Afterwards, functional enrichment analysis uncovered that CYFIP2 was involved in tumor-associated and immune-related pathways. Immune infiltration analysis indicated that CYFIP2 was significantly correlated with immune cells infiltration. In particular, CYFIP2 was strongly linked with immune microenvironment scores. Additionally, CYFIP2 exhibited a significant relationship with immune regulators and immune-related genes including chemokines, chemokines receptors, and MHC genes.ConclusionOur results suggested that CYFIP2 may serve as a prognostic cancer biomarker for determining prognosis and might be a promising therapeutic strategy for tumor immunotherapy.

SummaryObjectiveThe aim of the study was to characterise the protein complement of synovial fluid (SF) in health and osteoarthritis (OA) using liquid chromatography mass spectrometry (LC-MS/MS) following peptide-based depletion of high abundance proteins.DesignSF was used from nine normal and nine OA Thoroughbred horses. Samples were analysed with LC-MS/MS using a NanoAcquity™ LC coupled to an LTQ Orbitrap Velos. In order to enrich the lower-abundance protein fractions protein equalisation was first undertaken using ProteoMiner™. Progenesis-QI™ LC-MS software was used for label-free quantification. In addition immunohistochemistry, western blotting and mRNA expression analysis was undertaken on selected joint tissues.ResultsThe number of protein identifications was increased by 33% in the ProteoMiner™ treated SF compared to undepleted SF. A total of 764 proteins (462 with≥2 significant peptides) were identified in SF. A subset of 10 proteins were identified which were differentially expressed in OA SF. S100-A10, a calcium binding protein was upregulated in OA and validated with western blotting and immunohistochemistry. Several new OA specific peptide fragments (neopeptides) were identified.ConclusionThe protein equalisation method compressed the dynamic range of the synovial proteins identifying the most comprehensive SF proteome to date. A number of proteins were identified for the first time in SF which may be involved in the pathogenesis of OA. We identified a distinct set of proteins and neopeptides that may act as potential biomarkers to distinguish between normal and OA joints.

Objective The present study identified novel renal tubular biomarkers that may influence the diagnosis and treatment of focal segmental glomerulosclerosis (FSGS) based on immune infiltration. Methods Three FSGS microarray datasets, GSE108112, GSE133288 and GSE121211, were downloaded from the Gene Expression Omnibus (GEO) database. The R statistical software limma package and the combat function of the sva package were applied for preprocessing and to remove the batch effects. Differentially expressed genes (DEGs) between 120 FSGS and 15 control samples were identified with the limma package. Disease Ontology (DO) pathway enrichment analysis was conducted with statistical R software to search for related diseases. Gene set enrichment analysis (GSEA) was used to interpret the gene expression data and it revealed many common biological pathways. A protein-protein interaction (PPI) network was built using the Search Tool for the Retrieval of Interacting Genes (STRING) database, and hub genes were identified by the Cytoscape (version 3.7.2) plug-in CytoHubba. The plug-in Molecular Complex Detection (MCODE) was used to screen hub modules of the PPI network in Cytoscape, while functional analysis of the hub genes and hub nodes involved in the submodule was performed by ClusterProfiler. The least absolute shrinkage and selection operator (LASSO) regression and support vector machine recursive feature elimination (SVM-RFE) analysis were used to screen characteristic genes and build a logistic regression model. Receiver operating characteristic (ROC) curve analyses were used to investigate the logistic regression model and it was then validated by an external dataset GSE125779, which contained 8 FSGS samples and 8 healthy subjects. Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) was used to calculate the immune infiltration of FSGS samples. Results We acquired 179 DEGs, 79 genes with downregulated expression (44.1%) and 100 genes with upregulated expression (55.9%), in the FSGS samples. The DEGs were significantly associated with arteriosclerosis, kidney disease and arteriosclerotic cardiovascular disease. GSEA revealed that these gene sets were significantly enriched in allograft rejection signaling pathways and activation of immune response in biological processes. Fifteen genes were demonstrated to be hub genes by PPI, and three submodules were screened by MCODE linked with FSGS. Analysis by machine learning methodologies identified nuclear receptor subfamily 4 group A member 1 (NR4A1) and dual specificity phosphatase 1 (DUSP1) as sensitive tubular renal biomarkers in the diagnosis of FSGS, and they were selected as hub genes, as well as hub nodes which were enriched in the MAPK signaling pathway. Immune cell infiltration analysis revealed that the genetic biomarkers were both correlated with activated mast cells, which may amplify FSGS biological processes. Conclusion DUSP1 and NR4A1 were identified as sensitive potential biomarkers in the diagnosis of FSGS. Activated mast cells have a decisive effect on the occurrence and development of FSGS through tubular lesions and tubulointerstitial inflammation, and they are expected to become therapeutic targets in FSGS.