- 1Department of Gynecology, Southern Medical University Affiliated Maternal & Child Health Hospital of Foshan, Foshan, Guangdong, China
- 2Southern Medical University, Graduate School, Guangzhou, Guangdong, China
- 3Department of Gynecology, Shenzhen People’s Hospital, Shenzhen, Guangdong, China
Background: Endometriosis (EMs), a common gynecological disorder, adversely affects the quality of life of females. The pathogenesis of EMs has not been elucidated and the diagnostic methods for EMs have limitations. This study aimed to identify potential molecular biomarkers for the diagnosis and treatment of EMs.
Methods: Differential gene expression (DEG) and functional enrichment analyses were performed using the R language. WGCNA, Random Forest, SVM-REF and LASSO methods were used to identify core immune genes. The CIBERSORT algorithm was then used to analyse the differences in immune cell infiltration and to explore the correlation between immune cells and core genes. In addition, the extent of immune cell infiltration and the expression of immune core genes were investigated using single-cell RNA (scRNA) sequencing data. Finally, we performed molecular docking of three core genes with dienogest and goserelin to screen for potential drug targets.
Results: DEGs enriched in immune response, angiogenesis and estrogen processes. CXCL12, ROBO3 and SCG2 were identified as core immune genes. RT-PCR confirmed that the expression of CXCL12 and SCG2 was significantly upregulated in 12Z cells compared to hESCs cells. ROC curves showed high diagnostic value for these genes. Abnormal immune cell distribution, particularly increased macrophages, was observed in endometriosis. CXCL12, ROBO3 and SCG2 correlated with immune cell levels. Molecular docking suggested their potential as drug targets.
Conclusion: This study investigated the correlation between EMs and the immune system and identified potential immune-related biomarkers. These findings provided valuable insights for developing clinically relevant diagnostic and therapeutic strategies for EMs.
Introduction
Endometriosis (EMs), which is one of the prevalent gynecological disorderrs among women of reproductive age, is characterized by the occurrence of endometrial tissues (glands and stroma) outside the uterine cavity. Globally, EMs affects approximately 10% of women in their reproductive years with the highest incidence observed among women aged 25–35 years (1). Approximately 190 million women suffer from EMs worldwide (2). The prevalence rates of EMs in the asymptomatic, infertile, and pelvic pain groups are 2%–11%, 5%–50%, and 5%–21%, respectively. In symptomatic adolescents, the prevalence rate of EMs is in the range of 49%–75% (3). The symptoms of EMs include pain, development of masses and nodules, and infertility, which adversely affect the quality of life of females. The incidence of co-morbidities in females with EMs is higher than that in females without EMs. The infertility incidence rate is as high as 40%–50% in patients with EMs (2). The risk of early recurrent spontaneous abortion in pregnant women with EMs is significantly higher than that in healthy pregnant women (4).Furthermore, patients with EMs are at an increased risk of developing clear cell and endometrioid ovarian cancer, non-Hodgkin’s lymphoma, thyroid cancer, and endometrial cancer (3, 5, 6). Treatments for EMs include surgery, drug therapy, and assisted reproductive technologies. However, these treatments are not effective in some patients. In cases where the treatments are effective, approximately 50% of patients experience recurrence within five years of treatment (7) The average time lag between the onset of symptoms and the definitive diagnosis of EMs is 6–7 years (8). Currently used diagnostic methods for EMs are associated with advantages and disadvantages. Thus, there is a need to develop a combined diagnostic approach to reduce the delay and increase the precision of EMs diagnosis. Additionally, preventive strategies, as well as therapeutic strategies, for EMs must be developed to alleviate the suffering of patients. Furthermore, the molecular biomarkers for the diagnosis and treatment of EMs must be identified.
Several studies have examined the etiology of EMs although it has not been completely elucidated. Various hypotheses have been proposed for the underlying mechanisms and pathological determinants of EMs. In particular, the pivotal role of immune factors in EMs development has been the recent research focus (9–11). The occurrence of EMs is directly influenced by dysfunctional immune mechanisms and related factors (12, 13). Patients with EMs exhibit alterations in the peritoneal immune microenvironment and are characterized by augmented peritoneal fluid volume, compromised phagocytic function of peritoneal macrophages, reduced cytotoxicity of natural killer (NK) cells, variations in the quantity and functional properties of T lymphocytes and B lymphocytes, and upregulated levels of pro-inflammatory and chemotactic cytokines originating from macrophages and other cellular sources (9). This aberrant peritoneal immune environment can facilitate the immune escape of ectopic endometrial tissue, providing immune conditions for the implantation, invasion, and infiltration of endometriotic lesions and consequently inducing the development of EMs (14). Limited studies have examined the mechanisms underlying this process, emphasizing the need for further research in this domain. The elucidation of the correlation between immunological dysregulation and EMs will enable the identification of novel therapeutic targets and diagnostic modalities and contribute to alleviating the suffering of patients with EMs.
In this study, the microarray datasets were retrieved from the Gene Expression Omnibus (GEO) database (15) (GSE141549) and analyzed using R language. EMs was closely associated with pathways related to the immune system (including Antigen processing and presentation, T helper cell 17 (Th17 cell) differentiation, Th1 and Th2 cell differentiation, Leukocyte transendothelial migration), immune disorders (Rheumatoid arthritis and Systemic lupus erythematosus), and drug metabolism-cytochrome P450. Additionally, three potential immune-related biomarkers (C-X-C Motif Chemokine Ligand 12 (CXCL12), Secretogranin II (SCG2), Roundabout Guidance Receptor 3 (ROBO3)) for EMs were identified and a diagnostic model based on these biomarkers was developed. The expression levels CXCL12, SCG2, and ROBO3 were validated at the cellular level. Furthermore, the cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm and single-cell RNA sequencing (scRNA-seq) analyses were performed to examine the correlation between these biomarkers and immune cell infiltration. The findings of this study provided valuable insights for the development of diagnostic and therapeutic strategies for EMs.
Materials and methods
Data collection and processing
The GSE141549 (16), and GSE7305 (17) datasets were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) and used as the test cohort and validation datasets, respectively (Table 1). Next, scRNA-seq data (GSE179640 (18)), comprising the data of 9 eutopic endometrium samples (GSM6102537, GSM6102540, GSM6102543, GSM6102546, GSM6102549, GSM6102551, GSM6102554, GSM6102555, and GSM6102560) and 12 EMs samples (GSM6102536, GSM6102539, GSM6102542, GSM6102545, GSM6102548, GSM6102550, GSM6102552, GSM6102553, GSM6102556, GSM6102559, GSM6102561, and GSM6102562) were downloaded. The raw data from the GSE141549 dataset were processed using the “limma” package (19) in the R environment to standardize the data and remove the data of genes that were missing, duplicated, and highly downregulated. differential expression analysis was performed to identify the differentially expressed genes (DEGs) between the EMs and control groups. Volcano plots and heatmaps were used to visually depict the distribution and expression of DEGs. The DEGs were identified based on the following criteria: an absolute fold change greater than 1 (|log2FC|>1) and an adjusted p-value less than 0.05 (adj.P.Val<0.05). Next, the DEGs were subjected to functional enrichment analysis using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (20) modules of the OmicShare tool (https://www.omicshare.com/tools/).
Screening of the critical immune genes
The immune gene set was sourced from the IMMPORT website (http://www.immport.org/, Supplementary File 1). A scale-free weighted gene co-expression network of the GSE141549 dataset was constructed using the “Weighted Gene Coexpression Network Analysis (WGCNA)” toolkit (21) in R to identify co-expressed genes and modules associated with EMs. The intersection of the gene set within the module exhibiting the closest association with EMs, the DEGs, and the immune gene set facilitated the identification of key immune genes.
Key immune gene screening and diagnostic model construction
Random forest (RF) model, support vector machine-recursive feature elimination (SVM-RFE) algorithm, and least absolute shrinkage and selection operator regression (LASSO) regression were used to analyze the critical immune genes and screen the corresponding diagnostic genes. Key immune genes associated with EMs were identified using the intersection set. To examine the diagnostic efficacy of these genes, we calculated the area under the receiver operating characteristic (ROC) curve (AUC) for each key gene was calculated. The expression and diagnostic potential of these genes were validated using the GSE7305 datasets. Additionally, we conducted gene set enrichment analysis (GSEA) was performed to evaluate the potential functional roles of these genes in the context of EMs. Finally, a diagnostic column line graph was generated using the “Rms” package (22) in R to predict the occurrence of EMs. The consistency index (C-index) and calibration curve set decision curve analysis (DCA) were used to assess the predictive power and clinical application of the model.
Immune infiltration analysis
The CIBERSORT algorithm (23) was used to investigate differential immune cell infiltration between the EMs and control groups. Additionally, the correlation between immune cells and the key genes was examined using Spearman’s rank correlation coefficient.
scRNA-seq analysis
The “Seurat” package in R was used to analyze the scRNA-seq data. The data of low-quality cells, which were determined based on the following criteria, were censored from the analysis: cells expressing less than 200 genes and those expressing less than 3 genes; cells with < 1000 genes; cells with < 200 or > 10000 Unique Molecular Identifiers (UMIs) <200 or >10000; mitochondrial genes expressed in more than 25% of cells and ribosomal genes expressed in more than 20% of cells. Next, the “NormalizeData” function in Seurat was used to perform data normalization, followed by batch correction and dimensionality reduction using Harmony. The top 2000 highly variable genes were identified. Cell clusters were identified using 20 principal components (PCs) at resolution of 0.7, resulting in the subgrouping of all cells into 28 clusters. The DEGs for each cluster were identified using the “FindAllMarkers” function in Seurat with the following parameters: “min.pct = 0.25” and “logfc.threshold = 0.25”.
Single-cell annotation was performed using singleR. For accurate annotation, manual correction was performed, referencing established markers from published literature. The following marker genes were used for annotation: Platelet Derived Growth Factor Receptor Beta (PDGFRB) for (epithelial cells marker), Platelet And Endothelial Cell Adhesion Molecule 1 (PECAM1) and CD34 Molecule (CD34) for (endothelial cells marker), Collagen Type I Alpha 1 Chain (COL1A1), Collagen Type III Alpha 1 Chain (COL3A1), Collagen Type I Alpha 2 Chain (COL1A2), Fibroblast Growth Factor 7 (FGF7),and Membrane Metalloendopeptidasee (MME) for (fibroblasts markers), CD79a Molecule (CD68), CD19 Molecule (CD19), and Membrane Spanning 4-Domains A1 (MS4A1) for (B cells markers), CD3d Molecule (CD3D), CD3e Molecule (CD3E), and CD3g Molecule (CD3G) for (T cells markers), T Cell Receptor Delta Constant (TRDC), Killer Cell Lectin Like Receptor C1 (KLRC1), Killer Cell Lectin Like Receptor D1 (KLRD1), Natural Killer Cell Granule Protein 7 (NKG7), Killer Cell Lectin Like Receptor B1 (KLRB1), Neural Cell Adhesion Molecule 1 (NCAM1/CD56), and Killer Cell Lectin Like Receptor G1 (KLRG1) for (NK cells markers), CD14 Molecule (CD14) for Monocytes, CD68 Molecule (CD68) for (M1 macrophages marker), Colony Stimulating Factor 1 Receptor (CSF1R/CD115), CD163 Molecule (CD163), Macrophage Scavenger Receptor 1 (MSR1/CD204), and Mannose Receptor C-Type 1 (MRC1/CD206) for (M2 macrophages marker), Integrin Subunit Alpha X (ITGAX/CD11C) and CD1c Molecule (CD1C) for (dendritic cells marker), Membrane Spanning 4-Domains A2 (MS4A2), Glycoprotein IX Platelet (GP9/CD42a), and Glycoprotein Ib Platelet Subunit Alpha (GP1BA/CD42B) for (mast cells marker). Finally, a histogram of the content ratio of each cell, the distribution of EMs-related key immune genes expressed in single cells and gene differences was plotted.
Protein-ligand interaction analysis
The three-dimensional structure images of the three core target proteins were obtained from Uniprot (24) (https://www.uniprot.org/). The structures of dienogest and goserelin were downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov/) in the structure data file (SDF) format. Then OpenBabel software was used to transform the obtained SDF structures into mol2 structures. The protein receptors and small molecule ligands were imported into AutoDock software for molecular docking. The docking results were plotted using PyMol software.Binding pairs with binding energy values of < −1.2 kcal/mol were considered to exhibit good binding. The binding energy values were inversely proportional to the docking ability.
Quantitative real-time polymerase chain reaction
The eutopic stromal cell line (human embryonic stem cells (hESCs)) and the immortalized endometriotic epithelial cell line (12Z cells) were procured from ABM (Richmond, BC, Canada). Cells were cultured in Dulbecco’s modified Eagle’s medium/F12 (DMEM/F12; Gibco) supplemented with 10% fetal bovine serum (FBS; Gibco) and 1% penicillin-streptomycin (PS; Gibco) at 37°C and 5% CO2 in an incubator. Total RNA was extracted from hESCs and 12Z cells using the cell total RNA isolation kit (Foregene) and RNA extraction kits, following the manufacturer’s instructions. The RNA samples were subjected to reverse transcription using the 5× PrimeScript RT Master Mix (Takara). The resulting complementary DNA was quantified using the TB Green® Premix Ex Taq™ II (Takara) with the CFX96 touch q-PCR system (Bio-Rad Laboratories,Inc.). The relative expression levels of the target gene were determined using the 2−ΔΔCt. ACTB served as the internal reference. The experiments were repeated thrice for robustness. The primer sequences used in qRT-PCR analysis are shown in Supplementary Table 1.
Statistical analysis
All statistical analyses were performed using R software (version 4.1.2). The differential core genes expression and immune cell infiltration levels were compared using the Wilcoxon test. Spearman analysis was performed to evaluate the correlation analysis. Differences were considered significant at p < 0.05.
Results
DEG screening and functional enrichment analysis
After the standardization of the microarray data (GSE141549), 647 DEGs(247 downregulated genes and 399 upregulated genes) were identified using the “limma” package (Supplementary File 2) (Figures 1A, B). The biological functions of these DEGs were examined using GO annotation and KEGG pathway enrichment analyses. The DEGs were significantly enriched in immune-related processes, encompassing antigen processing and presentation, Th17 cell differentiation, Th1 and Th2 cell differentiation, and leukocyte transendothelial migration. Furthermore, the DEGs were associated with the pathways related to immune diseases (such as rheumatoid arthritis and systemic lupus erythematosus) and drug metabolism-cytochrome P450 (Figure 1C; Supplementary Table 2). The DEGs were enriched in different GO terms as follows: biological process terms, noteworthy enrichments were observed in critical processes, including blood vessel development, angiogenesis, regulation of the immune system process, regulation of leukocyte chemotaxis, and positive regulation of the immune system process; molecular function terms, chemokine receptor binding, cytokine activity, T cell receptor binding, estrogen receptor binding, and macrophage migration inhibitory factor binding; cellular component terms, extracellular structures and pathways related to macrophage migration (Figure 1D; Supplementary Table 3). The results of KEGG and GO revealed that EMs is associated with immune responses, vascular development, and estrogen-related processes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g001.jpg)
Figure 1 Screening of differentially expressed genes (DEGs) and functional enrichment analysis. (A) Heatmap of DEGs. (B) Volcano plot of DEGs. (C) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs. (D) Gene Ontology (GO) enrichment analysis of DEGs.
Weighted gene co-expression network construction
The data of samples in the GSE141549 database were clustered, and the soft threshold was set to 9 (Supplementary Figure 1B). A co-expression network was established when the R2 value was > 0.9, indicating high average connectivity (Figure 2A). Among the 14 gene modules identified, the MEmagenta module showed the highest correlation (0.82) with EMs (Figure 2B). The MEmagenta module comprised 407 genes (Supplementary File 3). The intersection of genes within the MEmagenta module DEGs, and immune-related genes revealed the following 20 critical immune genes associated with EMs: Bone Marrow Stromal Cell Antigen 2 (BST2), C-C Motif Chemokine Ligand 19 (CCL19), C-C Motif Chemokine Ligand 21 (CCL21), CXCL12, Fatty Acid Binding Protein 4 (FABP4), Growth Hormone Receptor (GHR), LIM Zinc Finger Domain Containing 1 (LIMS1), LDL Receptor Related Protein 1 (LRP1), Latent Transforming Growth Factor Beta Binding Protein 2 (LTBP2), Nerve Growth Factor (NGF), Platelet Derived Growth Factor Receptor Like (PDGFRL), Phospholipase A2 Group IIA (PLA2G2A), Prostaglandin F Receptor (PTGFR), ROBO3, Sphingosine-1-Phosphate Receptor 1 (S1PR1), SCG2, Semaphorin 3C (SEMA3C), Transforming Growth Factor Beta 3 (TGFB3), Transforming Growth Factor Beta Receptor 2 (TGFBR2), and Wnt Family Member 5A (WNT5A) (Figure 2C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g002.jpg)
Figure 2 Construction of weighted gene co-expression (WGCNA) network. (A) Screening the co -expression module. (B) Heatmap of module-trait correlations (red and blue indicate positive and negative correlations, respectively). (C) Venn plot of key module genes versus DEGs and immune genes.
Selection and validation of three core immune genes (CXCL12, ROBO3, and SCG2)
Immunological biomarkers with diagnostic significance were identified using three machine learning algorithms. The RF model identified seven genes (SCG2, ROBO3, FABP4, S1PR1, LRP1, CCL19, and CXCL12) based on the following criterion: importance scores > 4 (Figures 3A, B). Next, the SVM-RFE algorithm identified 16 genes (SCG2, ROBO3, S1PR1, FABP4, NGF, LIMS1, PLA2G2A, LRP1, TGFBR2, CCL19, PDGFRL, CXCL12, GHR, PTGFR, CCL21, and TGFB3) (Figure 3C). Finally, the LASSO regression analysis yielded six genes (TGFBR2, SCG2, ROBO3, CXCL12, NGF, and LIMS1) from the optimal model (Figures 3D, E). The intersection of these genes using a Venn diagram revealed three robust core biomarkers (CXCL12, ROBO3, and SCG2) (Figure 3F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g003.jpg)
Figure 3 Multiple machine learning models were used to identify core immune genes (CXCL12, ROBO3, and SCG2) in the GSE141549 dataset. (A, B) The genes were selected based on the random forest (RF) model. (C) The support vector machine (SVM) algorithm was used to select the genes. (D, E) Least absolute shrinkage and selection operator (LASSO) regression analysis was used to select the genes. (F) Three key immune genes identified from the Venn diagram.
Compared with those in the endometrium samples, the CXCL12, ROBO3, and SCG2 levels were significantly upregulated in the EMs samples (Figure 4A). The expression patterns of these three genes inEMs were validated using the GSE7305 dataset. Consistently, the CXCL12, ROBO3, and SCG2 were significantly upregulated in the EMs samples in the GSE7305 dataset (Figure 4C). qRT-PCR analysis was performed to further validate the dysregulation of CXCL12, ROBO3, and SCG2 EMs. The CXCL12 and SCG2 levels in 12Z cells were significantly upregulated when compared with those in hESCs (Figures 4E, F). Meanwhile, ROBO3 was downregulated in both 12Z cells and hESCs to an undetectable level. Furthermore, the diagnostic value of these genes was evaluated using ROC curve analysis. The AUC values for CXCL12, ROBO3, and SCG2 were almost all > 0.9 (Figures 4B, D). These findings indicate that CXCL12, ROBO3, and SCG2 are potential novel diagnostic biomarkers for EMs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g004.jpg)
Figure 4 Validating the expression and diagnostic value of core immune genes. (A, C) The expression levels of CXCL12, ROBO3, and SCG2 in the test (GSE141549) and validation cohorts (GSE7305). (B, D) Receiver operating characteristic (ROC) curves for evaluating the diagnostic values of CXCL12, ROBO3, and SCG2 in the test (GSE141549) and validation cohorts (GSE7305). (E, F) The expression levels of CXCL12 and SCG2 in 12Z cells and hESCs cells.
The pathways in which the three core immune genes were enriched in GSEA were similar to those in which DEGs were enriched (Figures 1C, D). These pathways included systemic lupus erythematosus, steroid hormone biosynthesis, leukocyte transendothelial migration, drug metabolism other enzymes, cytokine cytokine receptor interaction, and the Notch signaling pthway (Supplementary Figures 2A–C). Thus, CXCL12, ROBO3, and SCG2 are involved in essential immune responses and processes related to estrogen metabolism. In the validation set GSE7305, CXCL12, ROBO3, and SCG2 were significantly enriched in immune-related pathways, including cytokine-cytokine receptor interaction, chemokine signaling pathway, leukocyte transendothelial migration, systemic lupus erythematosus, and MAPK signaling pathway (Supplementary Figures 2D–F).
Establishment and testing of a diagnostic column line graph
The Rms package was used to construct an EMs diagnostic column line graph (Figure 5A). The calibration curve revealed minimal deviation between the actual event risk and the predicted risk, indicating the high accuracy of the model (Figure 5B). Furthermore, DCA revealed that the clinical net benefit of the diagnostic column line graph was higher than that of all other strategies (Figure 5C). As the high-risk threshold increased from 0.4 to 1, the “Number high risk” curve gradually overlapped with the “Number high risk with event” curve (Figure 5D). Additionally, the diagnostic column line graph increased high AUC values in the test cohort (GSE141549), and validation cohort (GSE7305) (0.981, 1.000, respectively; Figures 5E, F). These results validate the enhanced predictive performance of the diagnostic column line graph.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g005.jpg)
Figure 5 Construction and validation of the endometriosis diagnostic column line graph. (A) Diagnostic column line graph was used to predict the occurrence of endometriosis. (B) Calibration curve to examine the predictive power of the diagnostic column line graph. (C) Decision curve analysis (DCA) to evaluate the predictive power of the diagnostic column line graph. (D) Clinical impact curve to assess the predictive power of the diagnostic column line graph. (E, F) Receiver operating characteristic (ROC) curves to assess the clinical value of the diagnostic column line graph in the test (GSE141549) and validation cohorts (GSE7305).
Immune cell infiltration analysis
The CIBERSORT algorithm was used to estimate the infiltration abundance of 22 immune cell types in the samples of the GSE141549 dataset (Figure 6A). Comparative analysis revealed that the infiltration abundance of B cells memory, CD8+ T cells, resting CD4+ T cells memory, activated CD4+ T cells memory, T cells gamma delta, monocytes, M1 macrophages, M2 macrophages, and resting mast cells in the EMs samples were upregulated when compared with those in the endometrium tissue samples. In contrast, the infiltration abundances of plasma cells, regulatory T cells (Tregs), resting NK cells, activated NK cells, and activated dendritic cells were downregulated in the EMs samples (P<0.05, Figure 6B). The expression levels of CXCL12, ROBO3, and SCG2 were significantly and positively correlated with the abundances of B cells memory, M1 macrophages, M2 macrophages, resting mast cells, and T cells gamma delta. Conversely, the expression levels of CXCL12, ROBO3, and SCG2 were negatively correlated with the abundances of activated dendritic cells, activated NK cells, resting NK cells, plasma cells, and Tregs (Figure 7). These findings provided valuable insights into the correlation between immune cell infiltration and the expression patterns of CXCL12, ROBO3, and SCG2 in the context of EMs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g006.jpg)
Figure 6 Analysis of immune landscape related to endometriosis. Heatmap (A) and violin plot (B) showing the distribution of 22 immune cells in endometrium samples and endometriosis samples of the GSE141549 dataset. *P < 0.05; **P < 0.01; ***P < 0.001;****P < 0.0001; -, non-significant (P > 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g007.jpg)
Figure 7 Correlation between infiltrating immune cells and core immune genes (CXCL12 (A), ROBO3 (B), and SCG2 (C)).
scRNA profiling in EMs
The scRNA-seq of tissues obtained from subjects with EMs and eutopic endometrium, encompassing a total of 102,095 cells before filtering (60,663 cells from patients with EMs and 41,432 cells from controls (eutopic endometrium)), were analyzed. After stringent data quality control measures, a total of 36,881 cells were retained (18,206 cells from patients with EMs and 18,675 cells from controls) (Supplementary Figure 3A). The data were subjected to t-distributed stochastic neighbor embedding (t-SNE) analysis. The clusters were annotated using single-cell singleR annotation. To ensure accuracy, manual annotation correction based on published markers was performed (Supplementary Figure 3F). Various cell types, including fibroblasts, epithelial cells, endothelial cells, T cells, B cells, NK cells, dendritic cells, monocytes, M1 macrophages, M2 macrophages, and dematopoietic stem cell_-granulocyte colony-stimulating factor (HSC_-G-CSF) (Figures 8A–C). The proportions of different cell subtypes in each sample and tissue type are shown in Figure 8D and Table 2. Compared with those in the endometrium group, the proportions of M1 macrophages (10.65% vs. 5.64%), M2 macrophages (11.31% vs. 4.56%), B cells (1.55% vs. 0.81%), T cells (8.29% vs. 6.38%), and dendritic cells (2.41% vs. 1.15%)) were upregulated and the proportions of NK cells (14.64% vs. 16.05%) were downregulated in the EMs group. These results are consistent with those of immune cell infiltration analysis. Next, this study focused on the distribution of core genes (CXCL12, ROBO3, and SCG2) in t-SNE plots. Based on the t-SNE plot, CXCL12 and ROBO3 were predominantly expressed in the fibroblasts, epithelial cells, and M2 macrophage subclusters (Figure 8E). Furthermore, the expression levels of CXCL12, ROBO3, and SCG2 were upregulated in the EMs group, with CXCL12 and ROBO3 exhibiting significant differential expression (Figure 8F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g008.jpg)
Figure 8 Visualization plots of single-cell RNA sequencing (scRNA-seq) data of endometriosis samples. t-Distributed stochastic neighbor embedding (tSNE) plots showing integrated analysis of eutopic endometrium and endometriosis. Cells were colored according to cell types (A) or tissue types (B). (C) Violin plot of marker genes in different cell type. (D) Proportions of different T cell subtypes in each sample (left) or different tissue types (right). (E) The expression pattern of CXCL12, ROBO3, and SCG2 were representedd in the t-SNE plot. (F) The gene expression levels of CXCL12, ROBO3, and SCG2 in the scRNA-seq dataset (GSE179640). ****P < 0.0001; ns, non-significant (P > 0.05).
Molecular docking of three core immune genes
Dienogest and goserelin have been used for the treatment of EMs in several countries and regions worldwide. Molecular docking assesses the binding ability of the molecules to the target in terms of binding energy. The binding energy value was inversely proportional to the stability of the binding conformation. The three core genes could bind well to dienogest and goserelin, indicating that they have a favorable affinity (Figure 9; Supplementary Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1288263/fimmu-14-1288263-HTML/image_m/fimmu-14-1288263-g009.jpg)
Figure 9 Molecular docking patterns of dienogest or goserelin with target proteins. ALA, alanine; ASP, asparticacid; ARG, arginine; GLN, glutarnine; GLU, glutamicacid; GLY, glycine; PHE, phenylalanine; RPO, proline; SER, serine; THR, threonine.
Discussion
The societal costs of EMs, which adversely affects the quality of life of women and their families, are similar to those of chronic diseases, such as diabetes and rheumatoid arthritis (25). Although the exact mechanisms underlying EMs development have not been elucidated, the development of EMs is hypothesized to involve the interaction between multiple factors, especially the immune components (26). The understanding of these mechanisms is limited. Hence, early detection and intervention are critical for enhancing the quality of life of patients with EMs. Additionally, the diagnosis and treatment of EMs can be further improved by identifying early diagnostic biomarkers.
This study identified 647 DEGs in EMs. Functional enrichment analysis revealed that the DEGs were enriched in pathways related to immunity, angiogenesis, and estrogen. Leveraging WGCNA and various machine learning algorithms were used to identify three core immune-related diagnostic biomarkers (CXCL12, ROBO3, and SCG2). The roles of these core immune genes in EMs were examined using GSEA. The results of GSEA were consistent with those of enrichment analysis of DEGs, which revealed the enrichment of pathways related to the immune system (such as leukocyte transendothelial migration), immune disorders (such as systemic lupus erythematosus), and estrogen metabolism (such as steroid hormone biosynthesis, drug metabolism other enzymes).
These findings indicated the critical role of the immune system. Thus, the CIBERSORT algorithm was used to comprehensively analyze the infiltration levels of 22 immune cell types. Compared with those in the healthy endometrium tissue samples to Endometriosis samples, the infiltration B cells memory, CD8+ T cells, resting CD4+ T cells memory, activated CD4+ T cells memory, T cells gamma delta, monocytes, M1 macrophages, M2 macrophages, and resting mast cells were significantly upregulated and the infiltration levels of plasma cells, Tregs, resting NK cells, activated NK cells, and activated dendritic cells downregulated in the EMs samples. scRNA sequencing data analysis revealed 11 distinct cell type clusters. The levels of these 11 clusters varied between the endometrium group and the EMs group. Compared with those in the endometrium group, the proportions of M1 macrophages, M2 macrophages, B cells, T cells, and dendritic cells were upregulated and the proportions of NK cells were downregulated in the EMs group. Furthermore, the expression levels of CXCL12, ROBO3, and SCG2 were correlated with the infiltration levels of multiple immune cells. These findings indicate that CXCL12, ROBO3, and SCG2 are potential therapeutic targets for EMs.
Genetic and environmental factors are key factors involved in the development of EMs. In particular, immune dysfunction is involved during the entire pathological process of EMs (2, 27). Suryawanshi et al. (28) reported he presence of an immune environment within the peritoneal cavity of patients with EMs that was similar to tumor-like inflammatory environment. This milieu fosters the implantation, neovascularization, and proliferation of ectopic endometrial tissues, thereby facilitating the progression of EMs. Consistently, this study demonstrated the significant enrichment of pathways associated with the immune system (e.g. antigen processing and presentation, Th17 cell differentiation, Th1 and Th2 cell differentiation, and leukocyte transendothelial migration) and immune diseases (e.g., rheumatoid arthritis and systemic lupus erythematosus) in the EMs group (Figures 1C, D). These findings confirm the presence of an immune environment and an inflammation-prone environment in patients with EMs. Consistent with these findings, Chang KK et al. (29) revealed the effect of interleukin-10 (IL-10) and Th17 differentiation on the growth and invasion of endometriotic cells, consequently contributing to the progressive development and aggravation of the EMs. Additionally, several clinical studies have reported that systemic lupus erythematosus and rheumatoid arthritis increase the risk of EMs, suggesting that they have mutual influencing factors with similar immunological and inflammatory responses (5, 30, 31).
In this study, the distribution of most immune cells was dysregulated in EMs, which is consistent with the findings of previous studies (14, 32). In particular, EMs was associated with decreased cytotoxicity of NK cells, an increased abundance and cytotoxicity of macrophages, aberrant activation of T cells and B cells, and enhanced density and number of mast cells and dendritic cells. The activity of NK cells may be suppressed in endometriotic lesions facilitating the evasion of ectopic endometrial cells from immune surveillance and promoting their survival and implantation (33–35). Macrophages are the primary immune cells that produce pro-inflammatory chemokines and are the main source of neurovascularization (36). The disruption of the dynamic equilibrium between M1 and M2 macrophage phenotypes may be involved in the pathogenesis of EMs. Among women with EMs, M2 macrophages are enriched in the peritoneal environment (37, 38). M2 macrophages can exert pathogenic effects by mediating immune suppression, promoting neovascularization and tissue remodeling, and supporting the formation of endometriotic lesions and the fibrotic repair process (39–41). The peritoneal fluid microenvironment in patients with EMs specifically induces the differentiation of monocytes into macrophages rather than dendritic cells (42). This is consistent with the findings of this study. These findings suggest the predominance of macrophages relative to dendritic cells.
This study comprehensively demonstrated that CXCL12, ROBO3, and SCG2 are potential key biomarkers for EMs. The diagnostic values of CXCL12, ROBO3, and SCG2 were established using ROC curve analysis. The levels of CXCL12 in endometriotic lesions were upregulated when compared with those in healthy endometrial tissues, which was consistent with the findings of previous studies. CXCL12, which belongs the CXC chemokine family, is extensively expressed in diverse cells and tissues. CXCL12 effectively binds to CXC motif chemokine receptors 4/7 (CXCR4/7) and subsequently activates multiple signaling pathways (43–45). Additionally, CXCL12 is involved in cell homing, angiogenesis, immune and inflammatory responses, tumor infiltration, and metastasis. Recent studies have reported that CXCL12 and its receptors directly affect the proliferation, migration, invasion, and angiogenesis of endometrial stem cells in normal and ectopic endometrial tissue, promoting the occurrence and progression of endometriotic lesions (46, 47). Additionally, the CXCL12-CXCR4 interaction significantly inhibits the proliferation, migration, and invasion of ectopic endometrial cells (48).CXCL12 is significantly upregulated in endometriotic lesions and enhances the chemotactic activity of bone marrow stem cells (49). These findings are consistent with the results of this study as macrophages and mast cells differentiate from bone marrow stem cells. In this study, the proportions of M1 macrophages, M2 macrophages, and resting mast cells were positively correlated with CXCL12 expression. Consequently, CXCL12 may induce inflammatory reactions by recruiting immune cells to the ectopic lesions, contributing to the formation and development of endometriotic lesions. However, the expression patterns and functional roles of ROBO3 and SCG2 in EMs have not been examined. Limited studies have examined the function of ROBO3, a member of the roundabout guidance receptor (ROBO) family and a transmembrane protein mainly expressed in neuronal cells. The axon guidance ligand (SLIT)/ROBO signaling pathway is reported to play an important role in embryonic development, nervous system development, angiogenesis, and chemotaxis of inflammatory cells (50–52). Denk et al. performed cellular studies and reported that ROBO3 polymorphisms increase the susceptibility to rheumatoid arthritis. Compared with those in healthy individuals, the synovial fibroblast levels of ROBO3, pro-inflammatory cytokines, and cartilage degradation were significantly upregulated in patients with rheumatoid arthritis (53). Recent immunological studies have reported that autoimmune dysfunction may be involved in the development of EMs. Several clinical studies have confirmed that patients with EMs are associated with significantly increased incidence rates of autoimmune diseases, such as systemic lupus erythematosus and rheumatoid arthritis and that disease pathogenesis, clinical symptoms, and disease progression and regression are correlated (30, 31). This study demonstrated that the pathways related to systemic lupus erythematosus, chemokine signaling, and leukocyte transendothelial migration were significantly enriched in the high-ROBO3 expression group. Additionally, the infiltration levels of M2 macrophages and M1 macrophages were upregulated and the infiltration levels of activated dendritic cells, activated NK cells, resting NK cells, and Tregs were downregulated in the high-ROBO3 expression group. ROBO3 was hypothesized to contribute to aberrant immune responses, inflammation, and tissue damage that drive the progression of EMs. Secretoneurin (SN), which is the bioactive peptide product of SCG2 cleavage, is involved in neurotransmitter release, leukocyte migration, and angiogenesis (54, 55). SCG2/SN stimulates ovarian angiogenesis in response to human chorionic gonadotropin via a granulocyte signaling mechanism downstream of the luteinizing hormone receptor (56). In particular, SCG2/SN can prevent endothelial cell apoptosis by activating the vascular endothelial growth factor and MAPK signaling pathways and promote endothelial cell proliferation, migration, and angiogenesis (55, 57, 58), and inducing macrophage accumulation (59). In tumor tissues, SCG2 upregulation may promote the polarization of M0 macrophages to M2 macrophages and their subsequent differentiation into tumor-associated macrophages, enhancing tumor cell invasion and metastasis and angiogenesis and suppressing immune activity (60). In this study, SCG2 was upregulated in EMs tissues. Additionally, cytokine-cytokine receptor interaction, leukocyte transendothelial migration, and MAPK signaling pathway were significantly enriched in the high-SCG2 expression group. SCG2 expression was significantly and positively correlated with the abundances of M2 macrophages and M1 macrophages and negatively correlated with the abundances of activated NK cells and resting NK cells. The upregulation of SCG2 was hypothesized to affect the immune system and promote angiogenesis. These changes impair the ability of the immune system to effectively clear ectopic cells and enable adequate blood supply, promoting the survival and development of ectopic tissues. However, further studies are needed to determine the exact roles of CXCL12, ROBO3, and SCG2 in the pathogenesis of EMs.
Molecular docking analysis examines intermolecular interactions by evaluating the binding energies of small molecules in a network with protein ligands, enabling the mining of potential therapeutic targets. GSEA revealed that steroid hormone biosynthesis was significantly enriched in EMs tissues from the high-CXCL12 and high-ROBO3 expression groups. Dienogest (a novel and highly potent progestin) and goserelin (a gonadotropin-releasing hormone agonist (GnRH-a)) inhibit ovarian steroidogenesis through different mechanisms, downregulate the serum levels of endogenous steroids, and promote ectopic endometrial tissue decidualization, and then atrophy (61).The results of molecular docking analysis revealed that CXCL12, ROBO3, and SCG2 exhibited good binding affinity with dienogest and goserelin, suggesting that they are potential targets for dienogest and goserelin. However, further studies are needed to elucidate the specific regulatory mechanisms.
Conclusion
This study identified potential diagnostic biomarkers (CXCL12, ROBO3, and SCG2) for EMs and demonstrated the role of immune dysregulation in the EMs pathogenesis. These findings offer valuable insights for developing early diagnostic and targeted therapeutic strategies for EMs.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
Ethical approval was not required for the studies on humans in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used.
Author contributions
CZ: Conceptualization, Methodology, Software, Visualization, Writing – original draft. MF: Writing – original draft. YC: Writing – original draft. SL: Writing – original draft. YZ: Software, Visualization, Writing – original draft. JC: Funding acquisition, Software, Visualization, Writing – review & editing. RZ: Conceptualization, Writing – review & editing. XH: Conceptualization, Funding acquisition, Project administration, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by grants from the Foshan Transnatural Cavity Surgery Development and Innovation Engineering Technology Research Center (2020001003507) and the 2021 Foshan Self-funded Science and Technology Innovation Projects (2220001004062).
Acknowledgments
We thank the authors who provided the GSE141549, GSE7305, and GSE179640 public datasets. We thank Bullet Edits Limited for the linguistic editing and proofreading of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1288263/full#supplementary-material
References
1. Wei Y, Liang Y, Lin H, Dai Y, Yao S. Autonomic nervous system and inflammation interaction in endometriosis-associated pain. J Neuroinflamm (2020) 17(1):80. doi: 10.1186/s12974-020-01752-1
2. Zondervan KT, Becker CM, Missmer SA. Endometriosis. N Engl J Med (2020) 382(13):1244–56. doi: 10.1056/NEJMra1810764
3. Shafrir AL, Farland LV, Shah DK, Harris HR, Kvaskoff M, Zondervan K, et al. Risk for and consequences of endometriosis: A critical epidemiologic review. Best Pract Res Clin Obstet Gynaecol (2018) 51:1–15. doi: 10.1016/j.bpobgyn.2018.06.001
4. Saraswat L, Ayansina DT, Cooper KG, Bhattacharya S, Miligkos D, Horne AW, et al. Pregnancy outcomes in women with endometriosis: A national record linkage study. Bjog (2017) 124(3):444–52. doi: 10.1111/1471-0528.13920
5. Kvaskoff M, Mu F, Terry KL, Harris HR, Poole EM, Farland L, et al. Endometriosis: A high-risk population for major chronic diseases? Hum Reprod Update (2015) 21(4):500–16. doi: 10.1093/humupd/dmv013
6. Wang C, Liang Z, Liu X, Zhang Q, Li S. The association between endometriosis, tubal ligation, hysterectomy and epithelial ovarian cancer: meta-analyses. Int J Environ Res Public Health (2016) 13(11):1138. doi: 10.3390/ijerph13111138
7. Becker CM, Gattrell WT, Gude K, Singh SS. Reevaluating response and failure of medical treatment of endometriosis: A systematic review. Fertil Steril (2017) 108(1):125–36. doi: 10.1016/j.fertnstert.2017.05.004
8. Bedaiwy MA, Alfaraj S, Yong P, Casper R. New developments in the medical treatment of endometriosis. Fertil Steril (2017) 107(3):555–65. doi: 10.1016/j.fertnstert.2016.12.025
9. Chen S, Liu Y, Zhong Z, Wei C, Liu Y, Zhu X. Peritoneal immune microenvironment of endometriosis: role and therapeutic perspectives. Front Immunol (2023) 14:1134663. doi: 10.3389/fimmu.2023.1134663
10. Giudice LC, Kao LC. Endometriosis. Lancet (2004) 364(9447):1789–99. doi: 10.1016/s0140-6736(04)17403-5
11. Symons LK, Miller JE, Kay VR, Marks RM, Liblik K, Koti M, et al. The immunopathophysiology of endometriosis. Trends Mol Med (2018) 24(9):748–62. doi: 10.1016/j.molmed.2018.07.004
12. Wu MY, Ho HN. The role of cytokines in endometriosis. Am J Reprod Immunol (2003) 49(5):285–96. doi: 10.1034/j.1600-0897.2003.01207.x
13. Giacomini E, Minetto S, Li Piani L, Pagliardini L, Somigliana E, Viganò P. Genetics and inflammation in endometriosis: improving knowledge for development of new pharmacological strategies. Int J Mol Sci (2021) 22(16):9033. doi: 10.3390/ijms22169033
14. Vallvé-Juanico J, Houshdaran S, Giudice LC. The endometrial immune environment of women with endometriosis. Hum Reprod Update (2019) 25(5):564–91. doi: 10.1093/humupd/dmz018
15. NCBI Resource Coordinators. Database resources of the national center for biotechnology information. Nucleic Acids Res (2016) 44(D1):D7–19. doi: 10.1093/nar/gkv1290
16. Gabriel M, Fey V, Heinosalo T, Adhikari P, Rytkönen K, Komulainen T, et al. A relational database to identify differentially expressed genes in the endometrium and endometriosis lesions. Sci Data (2020) 7(1):284. doi: 10.1038/s41597-020-00623-x
17. Hever A, Roth RB, Hevezi P, Marin ME, Acosta JA, Acosta H, et al. Human endometriosis is associated with plasma cells and overexpression of B lymphocyte stimulator. Proc Natl Acad Sci U S A (2007) 104(30):12451–6. doi: 10.1073/pnas.0703451104
18. Tan Y, Flynn WF, Sivajothi S, Luo D, Bozal SB, Davé M, et al. Single-cell analysis of endometriosis reveals a coordinated transcriptional programme driving immunotolerance and angiogenesis across eutopic and ectopic tissues. Nat Cell Biol (2022) 24(8):1306–18. doi: 10.1038/s41556-022-00961-5
19. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for Rna-sequencing and microarray studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007
20. Yang H. Comprehensive analysis of the expression and clinical significance of a ferroptosis-related genome in ovarian serous cystadenocarcinoma: A study based on Tcga data. Oncologie (2022) 24(4):835–63. doi: 10.32604/oncologie.2022.026447
21. Langfelder P, Horvath S. Wgcna: an R package for weighted correlation network analysis. BMC Bioinf (2008) 9:559. doi: 10.1186/1471-2105-9-559
22. Deng M, Yin Y, Zhang Q, Zhou X, Hou G. Identification of inflammation-related biomarker lp-pla2 for patients with Copd by comprehensive analysis. Front Immunol (2021) 12:670971. doi: 10.3389/fimmu.2021.670971
23. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods (2015) 12(5):453–7. doi: 10.1038/nmeth.3337
24. UniProt Consortium. Uniprot: the universal protein knowledgebase in 2023. Nucleic Acids Res (2023) 51(D1):D523–d31. doi: 10.1093/nar/gkac1052
25. Zondervan KT, Becker CM, Koga K, Missmer SA, Taylor RN, Viganò P. Endometriosis. Nat Rev Dis Primers (2018) 4(1):9. doi: 10.1038/s41572-018-0008-5
26. Abramiuk M, Grywalska E, Małkowska P, Sierawska O, Hrynkiewicz R, Niedźwiedzka-Rystwej P. The role of the immune system in the development of endometriosis. Cells (2022) 11(13):2028. doi: 10.3390/cells11132028
27. Wang Y, Nicholes K, Shih IM. The origin and pathogenesis of endometriosis. Annu Rev Pathol (2020) 15:71–95. doi: 10.1146/annurev-pathmechdis-012419-032654
28. Suryawanshi S, Huang X, Elishaev E, Budiu RA, Zhang L, Kim S, et al. Complement pathway is frequently altered in endometriosis and endometriosis-associated ovarian cancer. Clin Cancer Res (2014) 20(23):6163–74. doi: 10.1158/1078-0432.Ccr-14-1338
29. Chang KK, Liu LB, Jin LP, Zhang B, Mei J, Li H, et al. Il-27 triggers il-10 production in th17 cells via a C-Maf/Rorγt/Blimp-1 signal to promote the progression of endometriosis. Cell Death Dis (2017) 8(3):e2666. doi: 10.1038/cddis.2017.95
30. Porpora MG, Scaramuzzino S, Sangiuliano C, Piacenti I, Bonanni V, Piccioni MG, et al. High prevalence of autoimmune diseases in women with endometriosis: A case-control study. Gynecol Endocrinol (2020) 36(4):356–9. doi: 10.1080/09513590.2019.1655727
31. Shigesi N, Kvaskoff M, Kirtley S, Feng Q, Fang H, Knight JC, et al. The association between endometriosis and autoimmune diseases: A systematic review and meta-analysis. Hum Reprod Update (2019) 25(4):486–503. doi: 10.1093/humupd/dmz014
32. Li W, Lin A, Qi L, Lv X, Yan S, Xue J, et al. Immunotherapy: A promising novel endometriosis therapy. Front Immunol (2023) 14:1128301. doi: 10.3389/fimmu.2023.1128301
33. Wu MY, Yang JH, Chao KH, Hwang JL, Yang YS, Ho HN. Increase in the expression of killer cell inhibitory receptors on peritoneal natural killer cells in women with endometriosis. Fertil Steril (2000) 74(6):1187–91. doi: 10.1016/s0015-0282(00)01592-2
34. Reis JL, Rosa NN, Ângelo-Dias M, Martins C, Borrego LM, Lima J. Natural killer cell receptors and endometriosis: A systematic review. Int J Mol Sci (2022) 24(1):331. doi: 10.3390/ijms24010331
35. Kang YJ, Jeung IC, Park A, Park YJ, Jung H, Kim TD, et al. An increased level of il-6 suppresses Nk cell activity in peritoneal fluid of patients with endometriosis via regulation of Shp-2 expression. Hum Reprod (2014) 29(10):2176–89. doi: 10.1093/humrep/deu172
36. Jeljeli M, Riccio LGC, Chouzenoux S, Moresi F, Toullec L, Doridot L, et al. Macrophage immune memory controls endometriosis in mice and humans. Cell Rep (2020) 33(5):108325. doi: 10.1016/j.celrep.2020.108325
37. Miller JE, Ahn SH, Marks RM, Monsanto SP, Fazleabas AT, Koti M, et al. Il-17a modulates peritoneal macrophage recruitment and M2 polarization in endometriosis. Front Immunol (2020) 11:108. doi: 10.3389/fimmu.2020.00108
38. Bacci M, Capobianco A, Monno A, Cottone L, Di Puppo F, Camisa B, et al. Macrophages are alternatively activated in patients with endometriosis and required for growth and vascularization of lesions in a mouse model of disease. Am J Pathol (2009) 175(2):547–56. doi: 10.2353/ajpath.2009.081011
39. Zhou WJ, Yang HL, Shao J, Mei J, Chang KK, Zhu R, et al. Anti-inflammatory cytokines in endometriosis. Cell Mol Life Sci (2019) 76(11):2111–32. doi: 10.1007/s00018-019-03056-x
40. Abomaray F, Gidlöf S, Götherström C. Mesenchymal stromal cells are more immunosuppressive in vitro if they are derived from endometriotic lesions than from eutopic endometrium. Stem Cells Int (2017) 2017:3215962. doi: 10.1155/2017/3215962
41. Artemova D, Vishnyakova P, Khashchenko E, Elchaninov A, Sukhikh G, Fatkhudinov T. Endometriosis and cancer: exploring the role of macrophages. Int J Mol Sci (2021) 22(10):5196. doi: 10.3390/ijms22105196
42. Na YJ, Jin JO, Lee MS, Song MG, Lee KS, Kwak JY. Peritoneal fluid from endometriosis patients switches differentiation of monocytes from dendritic cells to macrophages. J Reprod Immunol (2008) 77(1):63–74. doi: 10.1016/j.jri.2007.03.013
43. Walentowicz-Sadlecka M, Sadlecki P, Bodnar M, Marszalek A, Walentowicz P, Sokup A, et al. Stromal derived factor-1 (Sdf-1) and its receptors Cxcr4 and Cxcr7 in endometrial cancer patients. PLoS One (2014) 9(1):e84629. doi: 10.1371/journal.pone.0084629
44. Leconte M, Chouzenoux S, Nicco C, Chéreau C, Arkwright S, Santulli P, et al. Role of the Cxcl12-Cxcr4 axis in the development of deep rectal endometriosis. J Reprod Immunol (2014) 103:45–52. doi: 10.1016/j.jri.2013.12.121
45. Singh AK, Arya RK, Trivedi AK, Sanyal S, Baral R, Dormond O, et al. Chemokine receptor trio: Cxcr3, Cxcr4 and Cxcr7 crosstalk via Cxcl11 and Cxcl12. Cytokine Growth Factor Rev (2013) 24(1):41–9. doi: 10.1016/j.cytogfr.2012.08.007
46. Sahin Ersoy G, Zolbin MM, Cosar E, Moridi I, Mamillapalli R, Taylor HS. Cxcl12 promotes stem cell recruitment and uterine repair after injury in Asherman's syndrome. Mol Ther Methods Clin Dev (2017) 4:169–77. doi: 10.1016/j.omtm.2017.01.001
47. Hufnagel D, Li F, Cosar E, Krikun G, Taylor HS. The role of stem cells in the etiology and pathophysiology of endometriosis. Semin Reprod Med (2015) 33(5):333–40. doi: 10.1055/s-0035-1564609
48. Ruiz A, Ruiz L, Colón-Caraballo M, Torres-Collazo BJ, Monteiro JB, Bayona M, et al. Pharmacological blockage of the Cxcr4-Cxcl12 axis in endometriosis leads to contrasting effects in proliferation, migration, and invasion. Biol Reprod (2018) 98(1):4–14. doi: 10.1093/biolre/iox152
49. Moridi I, Mamillapalli R, Cosar E, Ersoy GS, Taylor HS. Bone marrow stem cell chemotactic activity is induced by elevated Cxcl12 in endometriosis. Reprod Sci (2017) 24(4):526–33. doi: 10.1177/1933719116672587
50. McConnell RE, Edward van Veen J, Vidaki M, Kwiatkowski AV, Meyer AS, Gertler FB. A requirement for filopodia extension toward slit during robo-mediated axon repulsion. J Cell Biol (2016) 213(2):261–74. doi: 10.1083/jcb.201509062
51. Klagsbrun M. Eichmann A. A role for axon guidance receptors and ligands in blood vessel development and tumor angiogenesis. Cytokine Growth Factor Rev (2005) 16(4-5):535–48. doi: 10.1016/j.cytogfr.2005.05.002
52. Wu JY, Feng L, Park HT, Havlioglu N, Wen L, Tang H, et al. The neuronal repellent slit inhibits leukocyte chemotaxis induced by chemotactic factors. Nature (2001) 410(6831):948–52. doi: 10.1038/35073616
53. Denk AE, Kaufmann S, Stark K, Schedel J, Lowin T, Schubert T, et al. Slit3 inhibits robo3-induced invasion of synovial fibroblasts in rheumatoid arthritis. Arthritis Res Ther (2010) 12(2):R45. doi: 10.1186/ar2955
54. Montero-Hadjadje M, Vaingankar S, Elias S, Tostivint H, Mahata SK, Anouar Y. Chromogranins a and B and secretogranin ii: evolutionary and functional aspects. Acta Physiol (Oxf) (2008) 192(2):309–24. doi: 10.1111/j.1748-1716.2007.01806.x
55. Kirchmair R, Gander R, Egger M, Hanley A, Silver M, Ritsch A, et al. The neuropeptide secretoneurin acts as a direct angiogenic cytokine in vitro and in vivo. Circulation (2004) 109(6):777–83. doi: 10.1161/01.Cir.0000112574.07422.C1
56. Hannon PR, Duffy DM, Rosewell KL, Brännström M, Akin JW, Curry TE Jr. Ovulatory induction of Scg2 in human, nonhuman primate, and rodent granulosa cells stimulates ovarian angiogenesis. Endocrinology (2018) 159(6):2447–58. doi: 10.1210/en.2018-00020
57. Albrecht-Schgoer K, Schgoer W, Holfeld J, Theurl M, Wiedemann D, Steger C, et al. The angiogenic factor secretoneurin induces coronary angiogenesis in a model of myocardial infarction by stimulation of vascular endothelial growth factor signaling in endothelial cells. Circulation (2012) 126(21):2491–501. doi: 10.1161/circulationaha.111.076950
58. Luo MJ, Rao SS, Tan YJ, Yin H, Hu XK, Zhang Y, et al. Fasting before or after Wound Injury Accelerates Wound Healing through the Activation of Pro-Angiogenic Smoc1 and Scg2. Theranostics (2020) 10(8):3779–92. doi: 10.7150/thno.44115
59. Liu W, Wang F, Zhao M, Fan Y, Cai W, Luo M. The neuropeptide secretoneurin exerts a direct effect on arteriogenesis in vivo and in vitro. Anat Rec (Hoboken) (2018) 301(11):1917–27. doi: 10.1002/ar.23929
60. Kim J, Bae JS. Tumor-associated macrophages and neutrophils in tumor microenvironment. Mediators Inflammation (2016) 2016:6058147. doi: 10.1155/2016/6058147
Keywords: endometriosis, single-cell RNA sequencing, CIBERSORT, immune, molecular docking
Citation: Zhou C, Feng M, Chen Y, Lv S, Zhang Y, Chen J, Zhang R and Huang X (2023) Unraveling immunotherapeutic targets for endometriosis: a transcriptomic and single-cell analysis. Front. Immunol. 14:1288263. doi: 10.3389/fimmu.2023.1288263
Received: 05 September 2023; Accepted: 31 October 2023;
Published: 16 November 2023.
Edited by:
Issam El Naqa, University of Michigan, United StatesReviewed by:
Arsalan Jalili, Royan Institute for Stem Cell Biology and Technology (RI-SCBT), IranKeyue Ma, Zai Lab (China), China
Copyright © 2023 Zhou, Feng, Chen, Lv, Zhang, Chen, Zhang and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiaobin Huang, xiaobinhuangcn@gmail.com; Rujian Zhang, erichy@126.com; Jiebo Chen, 351006256@qq.com
†These authors have contributed equally to this work and share first authorship