Abstract:Background: Endometriosis is recognized as an estrogen-dependent inflammation disorder, estimated to affect 8%-15% of women of childbearing age. Currently, the etiology and pathogenesis of endometriosis are not completely clear. Underlying mechanism for endometriosis is still under debate and needs further exploration. The involvement of transcription factors and immune mediations may be involved in the pathophysiological process of endometriosis, but the specific mechanism remains to be explored. This study aims to investigate the underlying molecular mechanisms in endometriosis. Methods: The gene expression profile of endometriosis was obtained from the gene expression omnibus (GEO) database. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were applied to the endometriosis GSE7305 datasets. Cibersort and MCP-counter were used to explore the immune response gene sets, immune response pathway, and immune environment. Differentially expressed genes (DEGs) were identified and screened. Common biological pathways were being investigated using the kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis. Transcription factors were from The Human Transcription Factors. The least absolute shrinkage and selection operator (Lasso) model identified four differential expressions of transcription factors (AEBP1, HOXB6, KLF2, and RORB). Their diagnostic value was calculated by receiver operating characteristic (ROC) curve analysis and validated in the validation cohort (GSE11691, GSE23339). By constructing the interaction network of crucial transcription factors, weighted gene coexpression network analysis (WGCNA) was used to search for key module genes. Metascape was used for enrichment analysis of essential module genes and obtained HOXB6, KLF2. The HOXB6 and KLF2 were further verified as the only two intersection genes according to Support Vector Machine Recursive Feature Elimination (SVM-RFE) and random forest models. We constructed ceRNA (lncRNA-miRNA-mRNA) networks with four potential transcription factors. Finally, we performed molecular docking for goserelin and dienogest with four transcription factors (AEBP1, HOXB6, KLF2, and RORB) to screen potential drug targets. Results: Immune and metabolic pathways were enriched in GSVA and GSEA. In single sample gene set enrichment analysis (ssGSEA), most immune infiltrating cells, immune response gene sets, and immune response pathways are differentially expressed between endometriosis and non-endometriosis. Twenty-seven transcription factors were screened from differentially expressed genes. Most of the twenty-seven transcription factors were correlated with immune infiltrating cells, immune response gene sets and immune response pathways. Furthermore, Adipocyte enhancer binding protein 1 (AEBP1), Homeobox B6 (HOXB6), Kruppel Like Factor 2 (KLF2) and RAR Related Orphan Receptor B (RORB) were selected out from twenty-seven transcription factors. ROC analysis showed that the four genes had a high diagnostic value for endometriosis. In addition, KLF2 and HOXB6 were found to play particularly important roles in multiple modules (String, WGCNA, SVM-RFE, random forest) on the gene interaction network. Using the ceRNA network, we found that NEAT1 may regulate the expressions of AEBP1, HOXB6 and RORB, while X Inactive Specific Transcript (XIST) may control the expressions of HOXB6, RORB and KLF2. Finally, we found that goserelin and dienogest may be potential drugs to regulate AEBP1, HOXB6, KLF2 and RORB through molecular docking. Conclusions: AEBP1, HOXB6, KLF2, and RORB may be potential biomarkers for endometriosis. Two of them, KLF2 and HOXB6, are critical molecules in the gene interaction network of endometriosis. Discovered by molecular docking, AEBP1, HOXB6, KLF2, and RORB are targets for goserelin and dienogest.

Cross-Tissue Regulatory Network Analyses Reveal Novel Susceptibility Genes and Potential Mechanisms for Endometriosis

Identification of Key Genes and Pathways in Endometriosis by Integrated Expression Profiles Analysis.

Pgpm_a_354957 377..392

Identification of MKNK1 and TOP3A As Ovarian Endometriosis Risk-Associated Genes Using Integrative Genomic Analyses and Functional Experiments

Single-cell Transcriptomic Analysis of Endometriosis Provides Insights into Fibroblast Fates and Immune Cell Heterogeneity

Long Non Coding RNA, C8orf49, a Novel Diagnostic and Prognostic Biomarker, Enhances PTEN/FZD4-mediated Cell Growth and Metastasis by Sponging Mir-1323 in Endometriosis.

Therapeutic targets for endometriosis: genome-wide Mendelian randomization and colocalization analyses

A Cross-Study Gene Set Enrichment Analysis Identifies Critical Pathways in Endometriosis

Identification of NFASC and CHL1 As Two Novel Hub Genes in Endometriosis Using Integrated Bioinformatic Analysis and Experimental Verification.

Bioinformatics identification and validation of biomarkers and infiltrating immune cells in endometriosis

Bioinformatic analysis reveals the importance of epithelial-mesenchymal transition in the development of endometriosis

Identification of Key Genes Associated with Endometriosis and Endometrial Cancer by Bioinformatics Analysis

Unraveling pathogenesis, biomarkers and potential therapeutic agents for endometriosis associated with disulfidptosis based on bioinformatics analysis, machine learning and experiment validation

Exploration of the molecular linkage between endometriosis and Crohn disease by bioinformatics methods

Molecular mechanisms underlying endometriosis pathogenesis revealed by bioinformatics analysis of microarray data

Bioinformatics Analysis Identifies Molecular Markers Regulating Development and Progression of Endometriosis and Potential Therapeutic Drugs

Gene expression analysis in endometriosis: Immunopathology insights, transcription factors and therapeutic targets

Identification of candidate genes for endometrial cancer in multi-omics: a Mendelian randomization analysis

Identification of central genes for endometriosis through integration of single-cell RNA sequencing and bulk RNA sequencing analysis

Construction and evaluation of endometriosis diagnostic prediction model and immune infiltration based on efferocytosis-related genes

Telomere-related genes as potential biomarkers to predict endometriosis and immune response: Development of a machine learning-based risk model