Identification and Validation of CYBB, CD86, and C3AR1 as the Key Genes Related to Macrophage Infiltration of Gastric Cancer

YTHDF2 favors protumoral macrophage polarization and implies poor survival outcomes in triple negative breast cancer

Patients with triple-negative breast cancer (TNBC) frequently experience resistance to chemotherapy, leading to recurrence. The approach of optimizing anti-tumoral immunological effect is promising in overcoming such resistance, given the heterogeneity and lack of biomarkers in TNBC. In this study, we focused on YTHDF2, an N6-methyladenosine (m6A) RNA-reader protein, in macrophages, one of the most abundant intra-tumoral immune cells. Using single-cell sequencing and ex vivo experiments, we discovered that YTHDF2 significantly promotes pro-tumoral phenotype polarization of macrophages and is closely associated with down-regulated antigen-presentation signaling to other immune cells in TNBC. The in vitro deprivation of YTHDF2 favors anti-tumoral effect. Expressions of multiple transcription factors, especially SPI1, were consistently observed in YTHDF2-high macrophages, providing potential therapeutic targets for new strategies. In conclusion, YTHDF2 in macrophages appears to promote pro-tumoral effects while suppressing immune activity, indicating the treatment targeting YTHDF2 or its transcription factors could be a promising strategy for chemoresistant TNBC.

Necroptosis and immune infiltration in hypertrophic cardiomyopathy: novel insights from bioinformatics analyses

Background

Hypertrophic Cardiomyopathy (HCM), a widespread genetic heart disorder, is largely associated with sudden cardiac fatality. Necroptosis, an emerging type of programmed cell death, plays a fundamental role in several cardiovascular diseases.

Aim

This research utilized bioinformatics analysis to investigate necroptosis's implication in HCM.

Methods

The study retrieved RNA sequencing datasets GSE130036 and GSE141910 from the Gene Expression Omnibus (GEO) database. It detected necroptosis-linked differentially expressed genes (NRDEGs) by reviewing both the gene set for necroptosis and the differently expressed genes (DEGs). The enriched signaling pathway of HCM was assessed using GSEA, while common DEGs were studied through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Concurrently, the Protein-Protein Interaction network (PPI) proved useful for identifying central genes. CIBERSORT facilitated evaluating the correlation between distinct immune cell-type prevalence and NRDEGs by analyzing immune infiltration patterns. Lastly, GSE141910 dataset validated the expression ranks of NRDEGs and immune-cell penetration.

Results

The investigation disclosed significant enrichment and activation of the necroptosis pathway in HCM specimens. Seventeen diverse genes, including CYBB, BCL2, and JAK2 among others, were identified in the process. PPI network scrutiny classified nine of these genes as central genes. Results from GO and KEGG enrichment analyses showed substantial connections of these genes to pathways pertaining to the HIF-1 signaling track, necroptosis, and NOD-like receptor signaling process. Moreover, an imbalance in M2 macrophage cells in HCM samples was observed. Finally, CYBB, BCL2, and JAK2 emerged as vital genes and were validated using the GSE141910 dataset.

Conclusion

These results indicate necroptosis as a probable underlying factor in HCM, with immune cell infiltration playing a part. Additionally, CYBB, BCL2, JAK2 could act as potential biomarkers for recognizing HCM. This information forms crucial insights into the basic mechanisms of HCM and could enhance its diagnosis and management.

Comprehensive analysis of lncRNA-associated ceRNA network reveals novel potential prognostic regulatory axes in glioblastoma multiforme

Deciphering the lncRNA-associated competitive endogenous RNA (ceRNA) network is essential in decoding glioblastoma multiforme (GBM) pathogenesis by regulating miRNA availability and controlling mRNA stability. This study aimed to explore novel biomarkers for GBM by constructing a lncRNA-miRNA-mRNA network. A ceRNA network in GBM was constructed using lncRNA, mRNA and miRNA expression profiles from the TCGA and GEO datasets. Seed nodes were identified by protein-protein interaction (PPI) network analysis of deregulated-mRNAs (DEmRNAs) in the ceRNA network. A lncRNA-miRNA-seed network was constructed by mapping the seed nodes into the preliminary ceRNA network. The impact of the seed nodes on the overall survival (OS) of patients was assessed by the GSCA database. Functional enrichment analysis of the deregulated-lncRNAs (DElncRNA) in the ceRNA network and genes interacting with OS-related genes in the PPI network were performed. Finally, the positive correlation between seed nodes and their associated lncRNAs and the expression level of these molecules in GBM tissue compared with normal samples was validated using the GEPIA database. Our analyzes revealed that three novel regulatory axes AL161785.1/miR-139-5p/MS4A6A, LINC02611/miR-139-5p/MS4A6A and PCED1B-AS1/miR-433-3p/MS4A6A may play essential roles in GBM pathogenesis. MS4A6A is upregulated in GBM and closely associated with shorter survival time of patients. We also identified that MS4A6A expression positively correlates with genes related to tumour-associated macrophages, which induce macrophage infiltration and immune suppression. The functional enrichment analysis demonstrated that DElncRNAs are mainly involved in neuroactive ligand-receptor interaction, calcium/MAPK signalling pathway, ribosome, GABAergic/Serotonergic/Glutamatergic synapse and immune system process. In addition, genes related to MS4A6A contribute to immune and inflammatory-related biological processes. Our findings provide novel insights to understand the ceRNA regulation in GBM and identify novel prognostic biomarkers or therapeutic targets.

Exosome as a crucial communicator between tumor microenvironment and gastric cancer (Review)

Gastric cancer is one of the most common malignancies and has relatively high morbidity and mortality rates. Exosomes are nanoscale extracellular vesicles that originate from a diverse array of cells and may be found throughout various bodily fluids. These vesicles are endogenous nanocarriers in their natural state with the unique ability to transport lipids, proteins, DNA and RNA. Exosomes contain DNA, RNA, proteins, lipids and other bioactive components that have crucial roles in the transmission of information and regulation of cell activities in gastric cancer. This paper begins with an exploration of the composition, formation and release mechanisms of exosomes. Subsequently, the role of exosomes in the tumor microenvironment is reviewed in terms of the immune cell population, nonimmune cell population and other factors. Finally, the current status and challenges of exosome‑based research on the progression, diagnosis and therapeutic methods of gastric cancer are summarized. This holistic review offers insight that may guide future research directions for exosomes and potentially pave the way for novel therapeutic interventions in the management of gastric cancer.

Construct dysregulated miRNA-mRNA interaction networks to conjecture possible pathogenesis for Stomach adenocarcinomas

BACKGROUND

Post-transcriptional regulation of mRNA induced by microRNA is known crucial in tumor occurrence, progression, and metastasis. This study aims at identifying significant miRNA-mRNA axes for stomach adenocarcinomas (STAD).

METHOD

RNA expression profiles were collected from The Cancer Genome Atlas (TCGA) and GEO database for screening differently expressed RNAs and miRNAs (DE-miRNAs/DE-mRNAs). Functional enrichment analysis was conducted with Hiplot and DAVID-mirPath. Connectivity MAP was applied in compounds prediction. MiRNA-mRNA axes were forecasted by TarBase and MiRTarBase. Real-time reverse transcription polymerase chain reaction (RT-qPCR) of stomach specimen verified these miRNA-mRNA pairs. Diagnosis efficacy of miRNA-mRNA interactions was measured by Receiver operation characteristic curve and Decision Curve Analysis. Clinical and survival analysis were also carried out. CIBERSORT and ESTIMATE was employed for immune microenvironment measurement.

RESULT

Totally 228 DE-mRNAs (105 upregulated and 123 downregulated) and 38 DE-miRNAs (22 upregulated and 16 downregulated) were considered significant. TarBase and MiRTarBase identified 18 miRNA-mRNA pairs, 12 of which were verified in RT-qPCR. The network of miR-301a-3p/ELL2 and miR-1-3p/ANXA2 were established and verified in external validation. The model containing all 4 signatures showed better diagnosis ability. Via interacting with M0 macrophage and resting mast cell, these miRNA-mRNA axes may influence tumor microenvironment.

CONCLUSION

This study established a miRNA-mRNA network via bioinformatic analysis and experiment validation for STAD.

Establishment and Characterization of Cell Lines from Canine Metastatic Osteosarcoma

Despite the advancements in treatments for other cancers, the outcomes for osteosarcoma (OSA) patients have not improved in the past forty years, especially in metastatic patients. Moreover, the major cause of death in OSA patients is due to metastatic lesions. In the current study, we report on the establishment of three cell lines derived from metastatic canine OSA patients and their transcriptome as compared to normal canine osteoblasts. All the OSA cell lines displayed significant upregulation of genes in the epithelial mesenchymal transition (EMT) pathway, and upregulation of key cytokines such as CXCL8, CXCL10 and IL6. The two most upregulated genes are MX1 and ISG15. Interestingly, ISG15 has recently been identified as a potential therapeutic target for OSA. In addition, there is notable downregulation of cell cycle control genes, including CDKN2A, CDKN2B and THBS1. At the protein level, p16INK4A, coded by CDKN2A, was undetectable in all the canine OSA cell lines, while expression of the tumor suppressor PTEN was variable, with one cell line showing complete absence and others showing low levels of expression. In addition, the cells express a variety of actionable genes, including KIT, ERBB2, VEGF and immune checkpoint genes. These findings, similar to those reported in human OSA, point to some genes that can be used for prognosis, targeted therapies and novel drug development for both canine and human OSA patients.

Bioinformatics analysis and prediction of Alzheimer's disease and alcohol dependence based on Ferroptosis-related genes

Background

Alzheimer's disease (AD) is a neurodegenerative disease whose origins have not been universally accepted. Numerous studies have demonstrated the relationship between AD and alcohol dependence; however, few studies have combined the origins of AD, alcohol dependence, and programmed cell death (PCD) to analyze the mechanistic relationship between the development of this pair of diseases. We demonstrated in previous studies the relationship between psychiatric disorders and PCD, and in the same concerning neurodegeneration-related AD, we found an interesting link with the Ferroptosis pathway. In the present study, we explored the bioinformatic interactions between AD, alcohol dependence, and Ferroptosis and tried to elucidate and predict the development of AD from this aspect.

Methods

We selected the Alzheimer's disease dataset GSE118553 and alcohol dependence dataset GSE44456 from the Gene Expression Omnibus (GEO) database. Ferroptosis-related genes were gathered through Gene Set Enrichment Analysis (GSEA), Kyoto Encyclopedia of Genes and Genomes (KEGG), and relevant literature, resulting in a total of 88 related genes. For the AD and alcohol dependence datasets, we conducted Limma analysis to identify differentially expressed genes (DEGs) and performed functional enrichment analysis on the intersection set. Furthermore, we used ferroptosis-related genes and the DEGs to perform machine learning crossover analysis, employing Least Absolute Shrinkage and Selection Operator (LASSO) regression to identify candidate immune-related central genes. This analysis was also used to construct protein-protein interaction networks (PPI) and artificial neural networks (ANN), as well as to plot receiver operating characteristic (ROC) curves for diagnosing AD and alcohol dependence. We analyzed immune cell infiltration to explore the role of immune cell dysregulation in AD. Subsequently, we conducted consensus clustering analysis of AD using three relevant candidate gene models and examined the immune microenvironment and functional pathways between different subgroups. Finally, we generated a network of gene-gene interactions and miRNA-gene interactions using Networkanalyst.

Results

The crossover of AD and alcohol dependence DEG contains 278 genes, and functional enrichment analysis showed that both AD and alcohol dependence were strongly correlated with Ferroptosis, and then crossed them with Ferroptosis-related genes to obtain seven genes. Three candidate genes were finally identified by machine learning to build a diagnostic prediction model. After validation by ANN and PPI analysis, ROC curves were plotted to assess the diagnostic value of AD and alcohol dependence. The results showed a high diagnostic value of the predictive model. In the immune infiltration analysis, functional metabolism and immune microenvironment of AD patients were significantly associated with Ferroptosis. Finally, analysis of target genes and miRNA-gene interaction networks showed that hsa-mir-34a-5p and has-mir-106b-5p could simultaneously regulate the expression of both CYBB and ACSL4.

Conclusion

We obtained a diagnostic prediction model with good effect by comprehensive analysis, and validation of ROC in AD and alcohol dependence data sets showed good diagnostic, predictive value for both AD (AUC 0. 75, CI 0.91-0.60), and alcohol dependence (AUC 0.81, CI 0.95-0.68). In the consensus clustering grouping, we identified variability in the metabolic and immune microenvironment between subgroups as a likely cause of the different prognosis, which was all related to Ferroptosis function. Finally, we discovered that hsa-mir-34a-5p and has-mir-106b-5p could simultaneously regulate the expression of both CYBB and ACSL4.

Clinical features and shared mechanisms of chronic gastritis and osteoporosis

Chronic gastritis (CG) and osteoporosis (OP) are common and occult diseases in the elderly and the relationship of these two diseases have been increasingly exposed. We aimed to explore the clinical characteristics and shared mechanisms of CG patients combined with OP. In the cross-sectional study, all participants were selected from BEYOND study. The CG patients were included and classified into two groups, namely OP group and non-OP group. Univariable and multivariable logistic regression methods were used to evaluate the influencing factors. Furthermore, CG and OP-related genes were obtained from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified using the GEO2R tool and the Venny platform. Protein-protein interaction information was obtained by inputting the intersection targets into the STRING database. The PPI network was constructed by Cytoscape v3.6.0 software again, and the key genes were screened out according to the degree value. Gene function enrichment of DEGs was performed by Webgestalt online tool. One hundred and thirty CG patients were finally included in this study. Univariate correlation analysis showed that age, gender, BMI and coffee were the potential influencing factors for the comorbidity (P < 0.05). Multivariate Logistic regression model found that smoking history, serum PTH and serum β-CTX were positively correlated with OP in CG patients, while serum P1NP and eating fruit had an negative relationship with OP in CG patients. In studies of the shared mechanisms, a total of 76 intersection genes were identified between CG and OP, including CD163, CD14, CCR1, CYBB, CXCL10, SIGLEC1, LILRB2, IGSF6, MS4A6A and CCL8 as the core genes. The biological processes closely related to the occurrence and development of CG and OP mainly involved Ferroptosis, Toll-like receptor signaling pathway, Legionellosis and Chemokine signaling pathway. Our study firstly identified the possible associated factors with OP in the patients with CG, and mined the core genes and related pathways that could be used as biomarkers or potential therapeutic targets to reveal the shared mechanisms.

Bioinformatics analysis and verification of key candidate genes influencing the pathogenesis of chronic rhinosinusitis with nasal polyps

OBJECTIVES

Chronic rhinosinusitis (CRS) with nasal polyps (CRSwNP) is a prominent public health issue. Furthermore, the prognosis of eosinophilic CRSwNP is poor, with a high recurrence rate. The underlying molecular mechanisms of eosinophilic CRSwNP remain unclear. Therefore, in this study, we sought to determine the crucial genes underlying eosinophil infiltration in eosinophilic CRSwNP pathogenesis.

METHODS

We used the Gene Expression Omnibus database (GEO) (GSE36830 and GSE23552 datasets) to mine gene expression profiles of CRSwNP patients and normal subjects. Differentially expressed genes (DEGs) between normal and CRSwNP tissues were identified and subjected to the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses. Co-expression networks were established using a weighted gene co-expression network analysis (WGCNA) and single-sample gene set enrichment analysis (GSEA). Protein-protein interaction networks were developed to detect functional protein modules. Based on the common DEGs, candidate miRNAs and related lncRNAs were predicted using the mirTarBase and StarBase databases. Finally, we generated immune cell subtypes of CRSwNP.

RESULTS

A total of 146 DEGs were identified. Of these, 131 genes were upregulated, whereas 15 were downregulated. GO analysis indicated that DEGs primarily participated in leukocyte chemotaxis and migration as well as cell chemotaxis. KEGG pathway analysis suggested that DEGs participated in the interactions between cytokines and viral proteins, osteoclast differentiation, and cytokine-cytokine receptor interactions. Real-time quantitative polymerase chain reaction analysis showed that Complement C5a Receptor 1 (C5AR1), C-C Motif Chemokine Receptor 3 (CCR3), Complement C3a Receptor 1 (C3AR1), and C-C Motif Chemokine Ligand 13 (CCL13) expression levels were significantly upregulated in nasal polyps, whereas C-C Motif Chemokine Ligand 4 (CCL4) expression levels were significantly downregulated.

CONCLUSIONS

The candidate genes identified in this study may influence the activation and accumulation of eosinophils, cell chemotaxis, and inflammatory responses, thereby potentially representing molecular targets for future studies of CRSwNP.

Identification and Validation of a Novel Ferroptotic Prognostic Genes-Based Signature of Clear Cell Renal Cell Carcinoma

Simple Summary

Clear cell renal cell carcinoma (ccRCC) is one of the leading types of kidney malignancy and is closely related to ferroptosis that is an iron-dependent regulated cell death with lipid peroxide accumulation. A signature of nine ferroptotic genes was identified as an independent prognostic factor via construction in The Cancer Genome Atlas (TCGA) database and validation in the ArrayExpress database. This signature could successfully divide patients into low- and high-risk groups to predict survival rate. Compared with the other eight genes, glutaminase 2 (GLS2) played a crucial role during erastin-induced ferroptosis in ACHN and Caki-1 cells. It was discovered for the first time that GLS2 might be a ferroptotic suppressor in ccRCC.

Abstract

Renal cell carcinoma (RCC), as one of the primary urological malignant neoplasms, shows poor survival, and the leading pathological type of RCC is clear cell RCC (ccRCC). Differing from other cell deaths (such as apoptosis, necroptosis, pyroptosis, and autophagy), ferroptosis is characterized by iron-dependence, polyunsaturated fatty acid oxidization, and lipid peroxide accumulation. We analyzed the ferroptosis database (FerrDb V2), Gene Expression Omnibus database, The Cancer Genome Atlas database, and the ArrayExpress database. Nine genes that were differentially expressed and related to prognosis were involved in the ferroptotic prognostic model via the least absolute shrinkage and selection operator Cox regression analysis, which was established in ccRCC patients from the kidney renal clear cell carcinoma (KIRC) cohort in TCGA database, and validated in ccRCC patients from the E-MTAB-1980 cohort in the ArrayExpress database. The signature could be an independent prognostic factor for ccRCC, and high-risk patients showed worse overall survival. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes were utilized to investigate the potential mechanisms. The nine genes in ccRCC cells with erastin or RSL3 treatment were validated to find the crucial gene. The glutaminase 2 (GLS2) gene was upregulated during ferroptosis in ccRCC cells, and cells with GLS2 shRNA displayed lower survival, a lower glutathione level, and a high lipid peroxide level, which illustrated that GLS2 might be a ferroptotic suppressor in ccRCC.

Expression profiles and functions of ferroptosis-related genes in intimal hyperplasia induced by carotid artery ligation in mice

Intimal hyperplasia (IH) is a prominent pathological event that occurs during in-stent restenosis and atherosclerosis. Ferroptosis, characterized by iron-dependent and lipid peroxidation, has become the recent focus of studies on the occurrence and progress of cardiovascular diseases. However, there are few studies on ferroptosis and IH. Therefore, we aimed to identify and validate ferroptosis-related markers in IH to explore new possibilities for IH diagnosis and treatment. The IH microarray dataset (GSE182291) was downloaded from the Gene Expression Omnibus (GEO) database and ferroptosis-related genes (FRGs) were obtained from the FerrDb databases. The differentially expressed genes (DEGs) were analyzed using the GEO2R. Overlapping was performed to identify the ferroptosis-related DEGs among the DEGs and FRGs. Then, clustering, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, and protein–protein interaction (PPI) analyses were performed. Subsequently, the hub genes were identified using Cytoscape and hub gene–transcription factors and hub gene–microRNA networks were constructed. Finally, real-time qPCR (RT-qPCR) and immunohistochemistry (IHC) were used to verify the mRNA and protein levels of the hub FRGs in IH. Thirty-four FRGs showing significantly different expression were identified from a total of 1,197 DEGs 2 days after ligation; 31 FRGs were selected from a total of 1,556 DEGs 14 days after ligation. The GO and KEGG analyses revealed that these 34 ferroptosis-related DEGs identified 2 days after ligation were mainly enriched in the basolateral plasma membrane, ferroptosis, lipid and atherosclerosis, and IL-17 signaling pathways. The 31 ferroptosis-related DEGs in endometrial hyperplasia identified 14 days after ligation were mainly enriched in response to oxidative stress, ferroptosis, tumor necrosis factor signaling pathway, and lipid and atherosclerosis. Five hub FRGs (Il1b, Ptgs2, Cybb, Cd44, and Tfrc) were identified using PPI networks; four hub FRGs (Il1b, Ptgs2, Cybb, and Cd44) were validated to be upregulated 2 and 14 days after ligation using RT-qPCR and show significantly different expression 14 days after ligation via IHC. Our findings verify the expression of hub DEGs related to ferroptosis in IH and elucidate the potential relationship between ferroptosis and IH, providing more evidence about the vital role of ferroptosis in IH.

Construction of a ferroptosis-associated circRNA-miRNA-mRNA network in age-related macular degeneration

Age-related macular degeneration (AMD) is a leading cause of severe vision impairment in the aging population. However, the underlying molecular mechanism remains unclear. Ferroptosis is a novel non-apoptotic programmed cell death pathway, that contributes to AMD. In addition, non-coding RNA-led epigenetic profile was identified in the regulation of AMD progression. Considering that non-coding RNAs are vital regulators of ferroptosis-related genes in various pathological events, we explored and constructed a ferroptosis-associated circRNA-miRNA-mRNA network in AMD. Differential expression of fourteen ferroptosis-associated genes were identified based on our microarray analysis and the FerrDb tool at the threshold of P < 0.05 and log₂|fold change| ≥ 1, which were subsequently validated by the public datasets. We further screened eight miRNAs via public datasets and the miRNet database. Based on these eight miRNAs, 23 circRNAs were mined using the Starbase tool. Taking all these together, we obtained a ferroptosis-related network with 414 pairs of circRNA-miRNA-mRNA, which are potential targets in future AMD treatments.

Identification of potential predictive biomarkers and biological pathways and the correction with immune infiltration in the activation of Crohn's disease

Crohn's disease (CD), a subtype of inflammatory bowel disease (IBD), has increasing prevalence in the world. Due to the lack of cure strategy, most patients with CD develop progressive disease companying with a series of serious complications. Therefore, exploring molecular mechanism differences between active and inactive CD will help in the screening of predict markers and therapeutic targets. In this study, we analyzed differentially expressed genes (DEGs) and molecular pathways through between active and inactive CD patients. In addition, the abundance of 22 immune cell types were assessed by using the CIBERSORT. The hub DEGs were screened out by the CytoHubba in Cytoscape, followed by the least absolute shrinkage and selection operator (LASSO) regression. Finally, the clinical predictive model was constructed by binary logistic regression model. The diagnostic efficacy was tested by receiver operating characteristic (ROC) curve and verified in independent datasets. The results showed that there were 137 DEGs between the active and inactive CD. Most of them were involved in regulating the immunity process. In addition, the decreased abundance of CD8 T cells and the increased abundance of M0, M1 macrophages, and neutrophils were closely related to CD activation. CXCL9, C3AR1, IL1B, and TLR4 were the hub gene and can be applied to the prediction of CD activation. Our results provided important targets for the prediction of CD activation and the selection of therapeutic targets.

The Roles of Blood Lipid-Metabolism Genes in Immune Infiltration Could Promote the Development of IDD

Objectives: Intervertebral disc degeneration is a progressive and chronic disease, usually manifesting as low back pain. This study aimed to screen effective biomarkers for medical practice as well as figuring out immune infiltration situations between circulation and intervertebral discs.

Methods: Gene expression profiles of GSE124272 was included for differentially analysis, WGCNA and immune infiltration analysis from GEO database, and other GSE series were used as validation datasets. A series of validation methods were conducted to verify the robustness of hub genes, such as principal component analysis, machine learning models, and expression verification. Lastly, nomogram was established for medical practice.

Results: 10 genes were commonly screened via combination of DEGs, WGCNA analysis and lipid metabolism related genes. Furthermore, 3 hub gens CYP27A1, FAR2, CYP1B1 were chosen for subsequent analysis based on validation of different methods. GSEA analysis discovered that neutrophil extracellular traps formation and NOD-like receptor signaling pathway was activated during IDD. Immune infiltration analysis demonstrated that the imbalance of neutrophils and γδT cells were significantly correlated with IDD progression. Nomogram was established based on CYP27A1, FAR2, CYP1B1 and age, the calibration plot confirmed the stability of our model.

Conclusion: CYP27A1, FAR2, CYP1B1 were considered as hub lipid metabolism related genes (LMRGs) in the development of IDD, which were regarded as candidate diagnostic biomarkers especially in circulation. The effects are worth expected in the early diagnosis of IDD through detecting these genes in blood.