A Novel Risk Score Model of Lactate Metabolism for Predicting over Survival and Immune Signature in Lung Adenocarcinoma

A machine learning model and identification of immune infiltration for chronic obstructive pulmonary disease based on disulfidptosis-related genes

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a chronic and progressive lung disease. Disulfidptosis-related genes (DRGs) may be involved in the pathogenesis of COPD. From the perspective of predictive, preventive, and personalized medicine (PPPM), clarifying the role of disulfidptosis in the development of COPD could provide a opportunity for primary prediction, targeted prevention, and personalized treatment of the disease.

METHODS

We analyzed the expression profiles of DRGs and immune cell infiltration in COPD patients by using the GSE38974 dataset. According to the DRGs, molecular clusters and related immune cell infiltration levels were explored in individuals with COPD. Next, co-expression modules and cluster-specific differentially expressed genes were identified by the Weighted Gene Co-expression Network Analysis (WGCNA). Comparing the performance of the random forest (RF), support vector machine (SVM), generalized linear model (GLM), and eXtreme Gradient Boosting (XGB), we constructed the ptimal machine learning model.

RESULTS

DE-DRGs, differential immune cells and two clusters were identified. Notable difference in DRGs, immune cell populations, biological processes, and pathway behaviors were noted among the two clusters. Besides, significant differences in DRGs, immune cells, biological functions, and pathway activities were observed between the two clusters.A nomogram was created to aid in the practical application of clinical procedures. The SVM model achieved the best results in differentiating COPD patients across various clusters. Following that, we identified the top five genes as predictor genes via SVM model. These five genes related to the model were strongly linked to traits of the individuals with COPD.

CONCLUSION

Our study demonstrated the relationship between disulfidptosis and COPD and established an optimal machine-learning model to evaluate the subtypes and traits of COPD. DRGs serve as a target for future predictive diagnostics, targeted prevention, and individualized therapy in COPD, facilitating the transition from reactive medical services to PPPM in the management of the disease.

Identification and Immunological Characterization of Cuproptosis Related Genes in Preeclampsia Using Bioinformatics Analysis and Machine Learning

Preeclampsia (PE) is a pregnancy-specific disorder characterized by an unclearly understood pathogenesis and poses a great threat to maternal and fetal safety. Cuproptosis, a novel form of cellular death, has been implicated in the advancement of various diseases. However, the role of cuproptosis and immune-related genes in PE is unclear. The current study aims to elucidate the gene expression matrix and immune infiltration patterns of cuproptosis-related genes (CRGs) in the context of PE. The GSE98224 dataset was obtained from the Gene Expression Omnibus (GEO) database and utilized as the internal training set. Based on the GSE98224 dataset, we explored the differentially expressed cuproptosis related genes (DECRGs) and immunological composition. We identified 10 DECRGs conducted Gene Ontology (GO) function, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, and a protein-protein interaction (PPI) network. Furthermore, patients with PE were categorized into two distinct clusters, and an investigation was conducted to examine the status of immune cell infiltration. Additionally, the application of Weighted Gene Co-expression Network Analysis (WGCNA) was utilized to differentiate modules consisting of co-expressed genes and conduct clustering analysis. The intersecting genes were obtained by intersecting differently expressed genes in PE and PE clusters. The most precise forecasting model was chosen by evaluating the effectiveness of four machine learning models. The ResNet model was established to score the hub genes. The prediction accuracy was assessed by receiver operating characteristic (ROC) curves and an external dataset. We successfully identified five key DECREGs and two pathological clusters in PE, each with distinct immune profiles and biological characteristics. Subsequently, the RF model was deemed the most optimal model for the identification of PE with a large area under the curve (AUC = 0.733). The five genes that ranked highest in the RF machine learning model were considered to be predictor genes. The calibration curve demonstrated a high level of accuracy in aligning the predicted outcomes with the actual outcomes. We validate the ResNet model using the ROC curve with the area under the curve (AUC = 0.82). Cuproptosis and immune infiltration may play an important role in the pathogenesis of PE. The present study elucidated that GSTA4, KCNK5, APLNR, IKZF2, and CAP2 may be potential markers of cuproptosis-associated PE and are considered to play a significant role in the initiation and development of cuproptosis-induced PE.

Prognostic value of a lactate metabolism gene signature in lung adenocarcinoma and its associations with immune checkpoint blockade therapy response

Lung adenocarcinoma (LUAD) is a study that examines the prognostic value of lactate metabolism genes in tumor cells, which are associated with clinical prognosis. We analyzed the expression and clinical data for LUAD from The Cancer Genome Atlas database, using the GSE68465 dataset from the Gene Expression Omnibus and the MSigDB database. LASSO Cox regression and stepwise Cox regression were used to identify the optimal lactate metabolism gene signature. Differences in immune infiltration, tumor mutation burden (TMB), and response to immune checkpoint blockade (ICB) therapy were evaluated between groups. LASSO and Cox regression analyses showed an eight-lactate metabolism gene signature for model construction in both TCGA cohort and GSE68465 data, with higher survival outcomes in high-risk groups. The lactate metabolism risk score had an independent prognostic value (hazard ratio: 2.279 [1.652-3.146], P < .001). Immune cell infiltration differed between the risk groups, such as CD8+ T cells, macrophages, dendritic cells, mast cells, and neutrophils. The high-risk group had higher tumor purity, lower immune and stromal scores, and higher TMB. High-risk samples had high tumor immune dysfunction and exclusion (TIDE) scores and low cytolytic activity (CYT) scores, indicating a poor response to ICB therapy. Similarly, most immune checkpoint molecules, immune inhibitors/stimulators, and major histocompatibility complex (MHC) molecules were highly expressed in the high-risk group. The 8-lactate metabolism gene-based prognostic model predicts patient survival, immune infiltration, and ICB response in patients with LUAD, driving the development of therapeutic strategies to target lactate metabolism.

Identification of immune infiltration and PANoptosis-related molecular clusters and predictive model in Alzheimer's disease based on transcriptome analysis

This study aims to explore the expression profile of PANoptosis-related genes (PRGs) and immune infiltration in Alzheimer's disease (AD). Based on the Gene Expression Omnibus database, this study investigated the differentially expressed PRGs and immune cell infiltration in AD and explored related molecular clusters. Gene set variation analysis (GSVA) was used to analyze the expression of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes in different clusters. Weighted gene co-expression network analysis was utilized to find co-expressed gene modules and core genes in the network. By analyzing the intersection genes in random forest, support vector machine, generalized linear model, and extreme gradient boosting (XGB), the XGB model was determined. Eventually, the first five genes (Signal Transducer and Activator of Transcription 3, Tumor Necrosis Factor (TNF) Receptor Superfamily Member 1B, Interleukin 4 Receptor, Chloride Intracellular Channel 1, TNF Receptor Superfamily Member 10B) in XGB model were selected as predictive genes. This research explored the relationship between PANoptosis and AD and established an XGB learning model to evaluate and screen key genes. At the same time, immune infiltration analysis showed that there were different immune infiltration expression profiles in AD.

Identification of copper death-associated molecular clusters and immunological profiles for lumbar disc herniation based on the machine learning

Lumbar disc herniation (LDH) is a common clinical spinal disorder, yet its etiology remains unclear. We aimed to explore the role of cuproptosis-related genes (CRGs) and identify potential diagnostic biomarkers. Our analysis involved interrogating the GSE124272 and GSE150408 datasets for differential gene expression profiles associated with CRGs and immune characteristics. Molecular clustering was performed on LDH samples, followed by expression and immune infiltration analyses. Using the WGCNA algorithm, specific genes within CRG clusters were identified. After selecting the most predictive genes from the optimal model, four machine learning models were constructed and validated. This study identified nine CRGs associated with copper-regulated cell death. Two copper-containing molecular clusters linked to death were detected in LDH samples. Elevated expression and immune infiltration levels were found in LDH patients, particularly in CRG cluster C2. Utilizing XGB, five genes were identified for constructing a diagnostic model, achieving an area under the curve values of 0.715. In conclusion, this research provides valuable insights into the association between LDH and copper-regulated cell death, alongside proposing a promising predictive model.

Comprehensive analysis of disulfidptosis-related genes reveals the effect of disulfidptosis in ulcerative colitis

Ulcerative colitis (UC) is a chronic inflammatory condition of the intestinal tract. Various programmed cell death pathways in the intestinal mucosa are crucial to the pathogenesis of UC. Disulfidptosis, a recently identified form of programmed cell death, has not been extensively reported in the context of UC. This study evaluated the expression of disulfidptosis-related genes (DRGs) in UC through public databases and assessed disulfide accumulation in the intestinal mucosal tissues of UC patients and dextran sulfate sodium (DSS)-induced colitis mice via targeted metabolomics. We utilized various bioinformatics techniques to identify UC-specific disulfidptosis signature genes, analyze their potential functions, and investigate their association with immune cell infiltration in UC. The mRNA and protein expression levels of these signature genes were confirmed in the intestinal mucosa of DSS-induced colitis mice and UC patients. A total of 24 DRGs showed differential expression in UC. Our findings underscore the role of disulfide stress in UC. Four UC-related disulfidptosis signature genes-SLC7A11, LRPPRC, NDUFS1, and CD2AP-were identified. Their relationships with immune infiltration in UC were analyzed using CIBERSORT, and their expression levels were validated by quantitative real-time PCR and western blotting. This study provides further insights into their potential functions and explores their links to immune infiltration in UC. In summary, disulfidptosis, as a type of programmed cell death, may significantly influence the pathogenesis of UC by modulating the homeostasis of the intestinal mucosal barrier.

Identification of a novel ADCC-related gene signature for predicting the prognosis and therapy response in lung adenocarcinoma

BACKGROUND

Previous studies have largely neglected the role of ADCC in LUAD, and no study has systematically compiled ADCC-associated genes to create prognostic signatures.

METHODS

In this study, 1564 LUAD patients, 2057 NSCLC patients, and more than 5000 patients with various cancer types from diverse cohorts were included. R package ConsensusClusterPlus was utilized to classify patients into different subtypes. A number of machine-learning algorithms were used to construct the ADCCRS. GSVA and ClusterProfiler were used for enrichment analyses, and IOBR was used to quantify immune cell infiltration level. GISTIC2.0 and maftools were used to analyze the CNV and SNV data. The Oncopredict package was used to predict drug information based on the GDSC1. Three immunotherapy cohorts were used to evaluate patient response to immunotherapy. The Seurat package was used to process single-cell data, the AUCell package was used to calculate cells' geneset activity scores, and the Scissor algorithm was used to identify ADCCRS-associated cells.

RESULTS

Through unsupervised clustering, two distinct subtypes of LUAD were identified, each exhibiting distinct clinical characteristics. The ADCCRS, consisted of 16 genes, was constructed by integrated machine-learning methods. The prognostic power of ADCCRS was validated in 28 independent datasets. Further, ADCCRS shows better predictive abilities than 102 previously published signatures in predicting LUAD patients' survival. A nomogram incorporating ADCCRS and clinical features was constructed, demonstrating high predictive performance. ADCCRS positively correlates with patients' gene mutation, and integrated analysis of bulk and single-cell transcriptome data revealed the association of ADCCRS with TME modulators. Cells representing high-ADCCRS phenotype exhibited more malignant features. LUAD patients with high ADCCRS levels exhibited sensitivity to chemotherapy and targeted therapy, while displaying resistance to immunotherapy. In pan-cancer analysis, ADCCRS still exhibited significant prognostic value and was found to be a risk factor for most cancer patients.

CONCLUSIONS

ADCCRS offers a critical prognostic insight for patients with LUAD, shedding light on the tumor microenvironment and forecasting treatment responsiveness.

Identification of cuproptosis-related gene clusters and immune cell infiltration in major burns based on machine learning models and experimental validation

Introduction

Burns are a global public health problem. Major burns can stimulate the body to enter a stress state, thereby increasing the risk of infection and adversely affecting the patient's prognosis. Recently, it has been discovered that cuproptosis, a form of cell death, is associated with various diseases. Our research aims to explore the molecular clusters associated with cuproptosis in major burns and construct predictive models.

Methods

We analyzed the expression and immune infiltration characteristics of cuproptosis-related factors in major burn based on the GSE37069 dataset. Using 553 samples from major burn patients, we explored the molecular clusters based on cuproptosis-related genes and their associated immune cell infiltrates. The WGCNA was utilized to identify cluster-specific genes. Subsequently, the performance of different machine learning models was compared to select the optimal model. The effectiveness of the predictive model was validated using Nomogram, calibration curves, decision curves, and an external dataset. Finally, five core genes related to cuproptosis and major burn have been was validated using RT-qPCR.

Results

In both major burn and normal samples, we determined the cuproptosis-related genes associated with major burns through WGCNA analysis. Through immune infiltrate profiling analysis, we found significant immune differences between different clusters. When K=2, the clustering number is the most stable. GSVA analysis shows that specific genes in cluster 2 are closely associated with various functions. After identifying the cross-core genes, machine learning models indicate that generalized linear models have better accuracy. Ultimately, a generalized linear model for five highly correlated genes was constructed, and validation with an external dataset showed an AUC of 0.982. The accuracy of the model was further verified through calibration curves, decision curves, and modal graphs. Further analysis of clinical relevance revealed that these correlated genes were closely related to time of injury.

Conclusion

This study has revealed the intricate relationship between cuproptosis and major burns. Research has identified 15 cuproptosis-related genes that are associated with major burn. Through a machine learning model, five core genes related to cuproptosis and major burn have been selected and validated.

Identification of cuproptosis-related molecular classification and characteristic genes in ulcerative colitis

Ulcerative colitis (UC) is a refractory inflammatory disease with imbalances in intestinal mucosal homeostasis. Cuproptosis serves as newly identified programmed cell death (PCD) form involved in UC. In the study, UC-related datasets were extracted from the Gene Expression Omnibus (GEO) database. A comparison of UC patients and healthy controls identified 11 differentially expressed cuproptosis-related genes (DE-CRGs), where FDX1, LIAS, and DLAT were differentially expressed in UC groups from the mouse models and clinical samples, with their expression correlating with disease severity. By comprehending weighted gene co-expression network analysis (WGCNA) and differential expression analysis, the key genes common to the module genes relevant to different cuproptosis-related clusters and differentially expressed genes (DEGs) both in different clusters and patients with and without UC were identified using several bioinformatic analysis. Furthermore, the mRNA levels of four characteristic genes with diagnostic potential demonstrated significant decrease in both mouse models and clinical UC samples. Our discoveries offer a theoretical foundation for cuproptosis effect in UC.

Identification and validation of eight lysosomes-related genes signatures and correlation with immune cell infiltration in lung adenocarcinoma

Lung cancer is the leading cause of cancer-related death. Lysosomes are key degradative compartments that maintain protein homeostasis. In current study, we aimed to construct a lysosomes-related genes signature to predict the overall survival (OS) of patients with Lung Adenocarcinoma (LUAD). Differentially expressed lysosomes-related genes (DELYs) were analyzed using The Cancer Genome Atlas (TCGA-LUAD cohort) database. The prognostic risk signature was identified by Least Absolute Shrinkage and Selection Operator (LASSO)-penalized Cox proportional hazards regression and multivariate Cox analysis. The predictive performance of the signature was assessed by Kaplan-Meier curves and Time-dependent receiver operating characteristic (ROC) curves. Gene set variant analysis (GSVA) was performed to explore the potential molecular biological function and signaling pathways. ESTIMATE and single sample gene set enrichment analysis (ssGSEA) were applied to estimate the difference of tumor microenvironment (TME) between the different risk subtypes. An eight prognostic genes (ACAP3, ATP8B3, BTK, CAV2, CDK5R1, GRIA1, PCSK9, and PLA2G3) signature was identified and divided patients into high-risk and low-risk groups. The prognostic signature was an independent prognostic factor for OS (HR > 1, p < 0.001). The molecular function analysis suggested that the signature was significantly correlated with cancer-associated pathways, including angiogenesis, epithelial mesenchymal transition, mTOR signaling, myc-targets. The low-risk patients had higher immune cell infiltration levels than high-risk group. We also evaluated the response to chemotherapeutic, targeted therapy and immunotherapy in high- and low-risk patients with LUAD. Furthermore, we validated the expression of the eight gene expression in LUAD tissues and cell lines by qRT-PCR. LYSscore signature provide a new modality for the accurate diagnosis and targeted treatment of LUAD and will help expand researchers' understanding of new prognostic models.

Identification of immune infiltration and cuproptosis-related molecular clusters in tuberculosis

Background

Tuberculosis (TB) is an infectious disease caused by Mycobacterium tuberculosis (Mtb) infection. Cuproptosis is a novel cell death mechanism correlated with various diseases. This study sought to elucidate the role of cuproptosis-related genes (CRGs) in TB.

Methods

Based on the GSE83456 dataset, we analyzed the expression profiles of CRGs and immune cell infiltration in TB. Based on CRGs, the molecular clusters and related immune cell infiltration were explored using 92 TB samples. The Weighted Gene Co-expression Network Analysis (WGCNA) algorithm was utilized to identify the co-expression modules and cluster-specific differentially expressed genes. Subsequently, the optimal machine learning model was determined by comparing the performance of the random forest (RF), support vector machine (SVM), generalized linear model (GLM), and eXtreme Gradient Boosting (XGB). The predictive performance of the machine learning model was assessed by generating calibration curves and decision curve analysis and validated in an external dataset.

Results

11 CRGs were identified as differentially expressed cuproptosis genes. Significant differences in immune cells were observed in TB patients. Two cuproptosis-related molecular clusters expressed genes were identified. Distinct clusters were identified based on the differential expression of CRGs and immune cells. Besides, significant differences in biological functions and pathway activities were observed between the two clusters. A nomogram was generated to facilitate clinical implementation. Next, calibration curves were generated, and decision curve analysis was conducted to validate the accuracy of our model in predicting TB subtypes. XGB machine learning model yielded the best performance in distinguishing TB patients with different clusters. The top five genes from the XGB model were selected as predictor genes. The XGB model exhibited satisfactory performance during validation in an external dataset. Further analysis revealed that these five model-related genes were significantly associated with latent and active TB.

Conclusion

Our study provided hitherto undocumented evidence of the relationship between cuproptosis and TB and established an optimal machine learning model to evaluate the TB subtypes and latent and active TB patients.

Identification of copper death-associated molecular clusters and immunological profiles in rheumatoid arthritis

Objective

An analysis of the relationship between rheumatoid arthritis (RA) and copper death-related genes (CRG) was explored based on the GEO dataset.

Methods

Based on the differential gene expression profiles in the GSE93272 dataset, their relationship to CRG and immune signature were analysed. Using 232 RA samples, molecular clusters with CRG were delineated and analysed for expression and immune infiltration. Genes specific to the CRGcluster were identified by the WGCNA algorithm. Four machine learning models were then built and validated after selecting the optimal model to obtain the significant predicted genes, and validated by constructing RA rat models.

Results

The location of the 13 CRGs on the chromosome was determined and, except for GCSH. LIPT1, FDX1, DLD, DBT, LIAS and ATP7A were expressed at significantly higher levels in RA samples than in non-RA, and DLST was significantly lower. RA samples were significantly expressed in immune cells such as B cells memory and differentially expressed genes such as LIPT1 were also strongly associated with the presence of immune infiltration. Two copper death-related molecular clusters were identified in RA samples. A higher level of immune infiltration and expression of CRGcluster C2 was found in the RA population. There were 314 crossover genes between the 2 molecular clusters, which were further divided into two molecular clusters. A significant difference in immune infiltration and expression levels was found between the two. Based on the five genes obtained from the RF model (AUC = 0.843), the Nomogram model, calibration curve and DCA also demonstrated their accuracy in predicting RA subtypes. The expression levels of the five genes were significantly higher in RA samples than in non-RA, and the ROC curves demonstrated their better predictive effect. Identification of predictive genes by RA animal model experiments was also confirmed.

Conclusion

This study provides some insight into the correlation between rheumatoid arthritis and copper mortality, as well as a predictive model that is expected to support the development of targeted treatment options in the future.

Characterization of lung adenocarcinoma based on immunophenotyping and constructing an immune scoring model to predict prognosis

Background: Lung cancer poses great threat to human health, and lung adenocarcinoma (LUAD) is the main subtype. Immunotherapy has become first line therapy for LUAD. However, the pathogenic mechanism of LUAD is still unclear. Methods: We scored immune-related pathways in LUAD patients using single sample gene set enrichment analysis (ssGSEA) algorithm, and further identified distinct immune-related subtypes through consistent clustering analysis. Next, immune signatures, Kaplan-Meier survival analysis, copy number variation (CNV) analysis, gene methylation analysis, mutational analysis were used to reveal differences between subtypes. pRRophetic method was used to predict the response to chemotherapeutic drugs (half maximal inhibitory concentration). Then, weighted gene co-expression network analysis (WGCNA) was performed to screen hub genes. Significantly, we built an immune score (IMscore) model to predict prognosis of LUAD. Results: Consensus clustering analysis identified three LUAD subtypes, namely immune-Enrich subtype (Immune-E), stromal-Enrich subtype (Stromal-E) and immune-Deprived subtype (Immune-D). Stromal-E subtype had a better prognosis, as shown by Kaplan-Meier survival analysis. Higher tumor purity and lower immune cell scores were found in the Immune-D subtype. CNV analysis showed that homologous recombination deficiency was lower in Stromal-E and higher in Immune-D. Likewise, mutational analysis found that the Stromal-E subtype had a lower mutation frequency in TP53 mutations. Difference in gene methylation (ZEB2, TWIST1, CDH2, CDH1 and CLDN1) among three subtypes was also observed. Moreover, Immune-E was more sensitive to traditional chemotherapy drugs Cisplatin, Sunitinib, Crizotinib, Dasatinib, Bortezomib, and Midostaurin in both the TCGA and GSE cohorts. Furthermore, a 6-gene signature was constructed to predicting prognosis, which performed better than TIDE score. The performance of IMscore model was successfully validated in three independent datasets and pan-cancer.

Construction of the model for predicting prognosis by key genes regulating EGFR-TKI resistance

Background: Previous studies have suggested that patients with lung adenocarcinoma (LUAD) will significantly benefit from epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKI). However, many LUAD patients will develop resistance to EGFR-TKI. Thus, our study aims to develop models to predict EGFR-TKI resistance and the LUAD prognosis. Methods: Two Gene Expression Omnibus (GEO) datasets (GSE31625 and GSE34228) were used as the discovery datasets to find the common differentially expressed genes (DEGs) in EGFR-TKI resistant LUAD profiles. The association of these common DEGs with LUAD prognosis was investigated in The Cancer Genome Atlas (TCGA) database. Moreover, we constructed the risk score for prognosis prediction of LUAD by LASSO analysis. The performance of the risk score for predicting LUAD prognosis was calculated using an independent dataset (GSE37745). A random forest model by risk score genes was trained in the training dataset, and the diagnostic ability for distinguishing sensitive and EGFR-TKI resistant samples was validated in the internal testing dataset and external testing datasets (GSE122005, GSE80344, and GSE123066). Results: From the discovery datasets, 267 common upregulated genes and 374 common downregulated genes were identified. Among these common DEGs, there were 59 genes negatively associated with prognosis, while 21 genes exhibited positive correlations with prognosis. Eight genes (ABCC2, ARL2BP, DKK1, FUT1, LRFN4, PYGL, SMNDC1, and SNAI2) were selected to construct the risk score signature. In both the discovery and independent validation datasets, LUAD patients with the higher risk score had a poorer prognosis. The nomogram based on risk score showed good performance in prognosis prediction with a C-index of 0.77. The expression levels of ABCC2, ARL2BP, DKK1, LRFN4, PYGL, SMNDC1, and SNAI2 were positively related to the resistance of EGFR-TKI. However, the expression level of FUT1 was favorably correlated with EGFR-TKI responsiveness. The RF model worked wonderfully for distinguishing sensitive and resistant EGFR-TKI samples in the internal and external testing datasets, with predictive area under the curves (AUC) of 0.973 and 0.817, respectively. Conclusion: Our investigation revealed eight genes associated with EGFR-TKI resistance and provided models for EGFR-TKI resistance and prognosis prediction in LUAD patients.

Immune micro-environment and drug analysis of peritoneal endometriosis based on epithelial-mesenchymal transition classification

BACKGROUND

Epithelial-mesenchymal transition (EMT) is a complex event that drives polar epithelial cells transform from adherent cells to motile mesenchymal cells, in which are involved immune cells and stroma cells. EMT plays crucial roles in migration and invasion of endometriosis. The interaction of endometrial implants with the surrounding peritoneal micro-environment probably affects the development of peritoneal endometriosis. To date, very few studies have been carried out on peritoneal endometriosis sub-type classification and micro-environment analysis based on EMT. The purpose of this study is to investigate the potential application of EMT-based classification in precise diagnosis and treatment of peritoneal endometriosis.

METHOD

Based on EMT hallmark genes, 76 peritoneal endometriosis samples were classified into two clusters by consistent cluster classification. EMT scores, which calculated by Z score of 8 epithelial cell marker genes and 8 mesenchymal cell marker genes, were compared in two clusters. Then, immune scores and the abundances of corresponding immune cells, stroma scores and the abundances of corresponding stroma cells were analyzed by the "xCell" package. Futhermore, a diagnostic model was constructed based on 9 diagnostic markers which related to immune score and stroma score by Lasso-Logistic regression analysis. Finally, based on EMT classification, a total of 8 targeted drugs against two clusters were screened out by drug susceptibility analysis via "pRRophetic" package.

RESULTS

Hallmark epithelial-mesenchymal transition was the mainly enriched pathway of differentially expressed genes between peritoneal endometriosis tissues and endometrium tissues. Compared with cluster 2, EMT score and the abundances of most infiltrating stroma cell were significantly higher, while the abundances of most infiltrating immune cells were dramatically less. The diagnostic model could accurately distinguish cluster 1 from cluster 2. Pathway analysis showed drug candidates targeting cluster 1 mainly act on the IGF-1 signaling pathway, and drug candidates targeting cluster 2 mainly block the EGFR signaling pathway.

CONCLUSION

In peritoneal endometriosis, EMT was probably promoted by stroma cell infiltration and inhibited by immune cell infiltration. Besides, our study highlighted the potential uses of the EMT classification in the precise diagnosis and treatment of peritoneal endometriosis.