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.

Common and Specific Gene Signatures among Three Different Endometriosis Subtypes

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

Elevated Immunoreactivity of RANTES and CCR1 Correlate with the Severity of Stages and Dysmenorrhea in Women with Deep Infiltrating Endometriosis.

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

Single-cell profiling identifies distinct hormonal, immunologic, and inflammatory signatures of endometriosis-constituting cells

Identification of Differentially Expressed Genes and Signaling Pathways Involved in Endometriosis by Integrated Bioinformatics Analysis.

Identification of Endometriosis-Associated Genes and Pathways Based on Bioinformatic Analysis

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

Integrated bioinformatic analysis reveals the gene signatures, epigenetic roles, and regulatory networks in endometriosis

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

Analysis of Differences in the Transcriptomic Profiles of Eutopic and Ectopic Endometriums in Women with Ovarian Endometriosis

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

In Search of Key Genes Associated with Endometriosis Using Bioinformatics Approach.

Bioinformatical analysis of the key differentially expressed genes and associations with immune cell infiltration in development of endometriosis

Identification and Analyzation of Differentially Expressed Transcription Factors in Endometriosis

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

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

Identification of genes and pathways associated with endometriosis using bioinformatic analysis