Identification and Analysis of Potential Immune-Related Biomarkers in Endometriosis

He, Yanan; Li, Jixin; Qu, Yanjun; Sun, Liyuan; Zhao, Xibo; Wu, Han; Zhang, Guangmei

doi:https://doi.org/10.1155/2023/2975581

Journal of Immunology Research

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Ethical Approval Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2023 | Article ID 2975581 | https://doi.org/10.1155/2023/2975581

Identification and Analysis of Potential Immune-Related Biomarkers in Endometriosis

Yanan He,¹Jixin Li,¹Yanjun Qu,¹Liyuan Sun,¹Xibo Zhao,¹Han Wu,¹and Guangmei Zhang¹

Academic Editor: Amar Singh

Received26 Oct 2022

Revised03 Dec 2022

Accepted06 Dec 2022

Published10 Jan 2023

Abstract

Background. Endometriosis is an inflammatory gynecological disease leading to deep pelvic pain, dyspareunia, and infertility. The pathophysiology of endometriosis is complex and depends on a variety of biological processes and pathways. Therefore, there is an urgent need to identify reliable biomarkers for early detection and accurate diagnosis to predict clinical outcomes and aid in the early intervention of endometriosis. We screened transcription factor- (TF-) immune-related gene (IRG) regulatory networks as potential biomarkers to reveal new molecular subgroups for the early diagnosis of endometriosis. Methods. To explore potential therapeutic targets for endometriosis, the Gene Expression Omnibus (GEO), Immunology Database and Analysis Portal (ImmPort), and TF databases were used to obtain data related to the recognition of differentially expressed genes (DEGs), differentially expressed IRGs (DEIRGs), and differentially expressed TFs (DETFs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the DETFs and DEIRGs. Then, DETFs and DEIRGs were further validated in the external datasets of GSE51981 and GSE1230103. Then, we used quantitative real-time polymerase chain reaction (qRT-PCR) to verify the hub genes. Simultaneously, the Pearson correlation analysis and protein-protein interaction (PPI) analyses were used to indicate the potential mechanisms of TF-IRGs at the molecular level and obtain hub IRGs. Finally, the receiver operating characteristic (ROC) curve analysis was used to assess the diagnostic value of the hub IRGs. Results. We screened a total of 94 DETFs and 121 DEIRGs in endometriosis. Most downregulated DETFs showed decreased expression in the endometria of moderate/severe endometriosis patients. The top-ranked upregulated DEIRGs were upregulated in the endometra of infertile women. Functional analysis showed that DETFs and DEIRGs may be involved in the biological behaviors and pathways of endometriosis. The TF-IRG PPI network was successfully constructed. Compared with the control group, high C3, VCAM1, ITGB2, and C3AR1 expression had statistical significance in endometriosis among the hub DEIRGs. They also showed higher sensitivity and specificity by ROC analysis for the diagnosis of endometriosis. Finally, compared with controls, C3 and VCAM1 were highly expressed in endometriosis tissue samples. In addition, they also showed high specificity and sensitivity for diagnosing endometriosis. Conclusion. Overall, we discovered the TF-IRG regulatory network and analyzed 4 hub IRGs that were closely related to endometriosis, which contributes to the diagnosis of endometriosis. Additionally, we verified that DETFs or DEIRGs were associated with the clinicopathological features of endometriosis, and external datasets also confirmed the hub IRGs. Finally, C3 and VCAM1 were highly expressed in endometriosis tissue samples compared with controls and may be potential biomarkers of endometriosis, which are helpful for the early diagnosis of endometriosis.

1. Introduction

Endometriosis is an inflammatory gynecological disease characterized by the presence of endometrial tissues outside the uterus [1]. It affects approximately 10% of females in their reproductive years leading to a heavy financial burden on patients [2]. The typical clinical symptoms are chronic pelvic pain, dysmenorrhea, dyspareunia, and infertility, causing a decrease in patients’ quality of life [3]. In addition, endometriosis surgery is the second most common surgery in premenopausal women. The occurrence and development of endometriosis are usually slow and are related to the local inflammatory response, proliferation, invasion, and angiogenesis of endometrial cells [4]. There are many theories about the etiology of endometriosis, but no exact theory can explain the pathogenesis of endometriosis [5]. Among the theories, the most prevalent is Sampson’s theory of retrograde menstruation proposed in 1921. However, there are many arguments against this theory [6]. Because 90% of women have menstrual reflux, only 10% of women suffer from endometriosis. Although our understanding of endometriosis is growing, the exact molecular mechanisms underlying this tumor-like disease are still far from being understood. The pathophysiology of the occurrence and development of endometriosis is complex and depends on a variety of biological processes and pathways [3]. Therefore, there is an urgent need to determine reliable diagnostic biomarkers to predict early diagnosis and clinical severity.

The immune system plays a major role in survival in the pelvic microenvironment, including causing immune tolerance, depressing immunosurveillance, and escaping phagocytosis by immune cells [7]. Previous studies have indicated that immune-related genes (IRGs) play an important role in the complex regulatory network of tumors [8], and they have been explored to indicate the development of tumor immunity and the pathophysiological mechanisms of tumors [9], such as ovarian cancer. Emerging evidence has shown that women with endometriosis not only have a changed immune status of the endometrium but also have an altered peripheral immune system [10]. Consistent with the changes in the peritoneal environment of endometriosis, a large number of immune cells, inflammatory factors, and relevant cytokines have also been recruited to contribute to the abnormal immune environment in endometriosis [11, 12]. Nevertheless, the relationship between IRGs and the diagnosis of endometriosis patients is not clear, and further research is needed. This study is aimed at exploring the potential correlation between the onset of endometriosis and IRGs, which are potential molecular biomarkers to reveal new molecular subgroups for the early diagnosis of endometriosis.

Some transcription factors (TFs) are closely related to IRGs and can also regulate the function of IRGs in some diseases. Aberrant TF-IRGs could influence the various processes of tumor development. Additionally, the differential expression of TFs and their downstream target genes has been found to be related to the progression of endometriosis. Previous studies have shown that IRGs act as important regulators in diverse pathological processes. Therefore, studying the role of IRGs and their related molecular mechanisms in endometriosis is crucial, which is beneficial for exploring the pathogenesis of endometriosis and detecting more effective potential diagnostic markers.

2. Materials and Methods

2.1. Preparation and Processing of TF and IRG Data in Endometriosis

We searched two online TF datasets and downloaded 1665 TFs from the Human Transcription Factor Database (Human TFDB) [13] and 1639 TFs from the Human Transcription Factors Database [14]. The 1508 TFs obtained by the intersection of these two TF databases were used in our research on endometriosis. We constructed a diagnostic signature from the IRGs downloaded from the Immunology Database and Analysis Portal (ImmPort) database (http://www.immport.org) [15]. We used the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) to analyze gene expression datasets. A total of 1871 series of endometriosis were retrieved from the database. We selected five GEO datasets (GSE7305, GSE7307, GSE51981, GSE1230103, and GSE23339) after filtering.

We matched the gene symbols of the data with the corresponding GEO platforms (GPL). In total, 10 endometriosis cases and 10 control samples were obtained from GSE7305, while 23 endometriosis patients’ specimens and 18 control specimens were acquired from GSE7307. Both two expression microarrays were based on the GPL570 ((HG-U133_Plus_2) Affymetrix Human Genome U133 Plus 2.0 Array) platform. Moreover, GSE51981 and GSE120103, based on the GPL570 ((HG-U133_Plus_2) Affymetrix Human Genome U133 Plus 2.0 Array) and GPL6480 (Agilent-014850 Whole Huma Genome Microarray 4x44K G4112F) platform, respectively, were chosen for further validation. All of the data are freely available online.

2.2. Identification of DETFs and DEIRGs in Endometriosis

We selected two GEO gene expression datasets (GSE7305 and GSE7307) and divided the above data into the endometriosis group and the control group. First, the differentially expressed genes (DEGs) between the endometriosis and control samples were identified using the GEO2R online analysis tool (https://www.ncbi.nlm.nih.gov/geo/geo2r/), and the genes met the cutoff criteria based on the criteria of and [16]. Statistical analysis was carried out on each dataset, and the Venn diagram webtool (http://bioinformatics.psb.ugent.be/webtools/Venn/) was used to identify the intersection. The DEGs were then intersected with TFs from both databases and IRGs from the ImmPort database to obtain coupregulated differentially expressed TFs (DETFs), codownregulated DETFs, coupregulated differentially expressed IRGs (DEIRGs), and codownregulated DEIRGs.

2.3. GO Enrichment and KEGG Pathway Analyses of DETFs and DEIRGs

In this study, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of DETFs and DEIRGs were performed using the ClusterProfiler R package (version 3.18.0) [17], and was considered statistically significant. The R package can automate the process of biological term classification and the enrichment analysis of gene clusters to unravel the biological meaning behind a large list of genes. GO analysis can be mainly classified into three domains: biological process (BP), molecular function (MF), and cellular component (CC). The version of R used in our research was 4.0.3.

2.4. PPI Network Construction and Hub Gene Identification

We used the Search Tool for the Retrieval of Interacting Genes (STRING) database (http://string-db.org/) to analyze protein-protein interaction (PPI) networks. It was essential to display the molecular mechanisms of key activities in endometriosis. To investigate the potential PPI relationships, the previously identified DETFs and DEIRGs were mapped to the STRING database. The PPIs were extracted with a combined . Subsequently, the visualized PPI network was constructed by Cytoscape software (version 3.7.1). The Hmisc R package (version 4.4.2) (https://hbiostat.org/R/doc/sintro.pdf) was utilized to test the correlations between DETFs and DEIRGs with the cutoff criteria set as correlation and . The Molecular Complex Detection (MCODE) plugin of Cytoscape software was utilized to recognize the most prominent clustering modules. Functional enrichment analysis of the genes in individual modules was achieved by DAVID, an online tool (https://david.ncifcrf.gov/), with a significance threshold of . Nodes with higher connectivity tend to be more important for maintaining the stability of the entire network. Therefore, cytoHubba, a plugin in Cytoscape, was used to screen out hub genes.

2.5. Collection of Human Tissues

Ectopic endometrium tissues were collected from chocolate cyst in endometriosis to identify hub genes expression (). Endometrium tissue from patients undergoing surgery for uterine fibroids served as a control group (). All patients in our study with or without endometriosis had no menstrual disorders. Patients who had received hormone therapy or other serious diseases were not included in this study. All tissue samples obtained were approved by the Ethics Commission of Harbin Medical University (202106).

2.6. Reverse Transcription and Quantitative Real-Time Polymerase Chain Reaction

12 endometriosis samples and 12 controls were frozen in liquid nitrogen, and total RNA was extracted using the TRIzol^® reagent (15596026, America). The relative expression of VCAM1, ITGB2, C3AR1, and C3 mRNA was normalized to GAPDH, and calculated using the method (). The total RNA was used only if the A260/280 ratio of the absorbances was between 1.8 and 2.2 when measured by spectrophotometry. Reverse transcription was performed at 42°C (15 min) followed by 95°C (3 min), then in a 10 ml SYBR reaction system using the Talent qPCR PreMix (FP209-02, China) with 1 cycle of 95°C for 3 minutes, and 40 cycles of 95°C for 5 seconds and 60°C for 15 seconds. We identified the target mRNA sequences with ideal melting curves and sizes. Sequences of the primers are shown in Table 1.

2.7. Statistical Analysis

The method used to compare DETF expression in different severity groups of endometriosis was unpaired Student’s test. At the same time, the comparison method between the expression of DEIRGs in infertile and fertile endometriosis was unpaired Student’s test. The pROC R package (version 1.18.0) [18] was used to evaluate the sensitivity and specificity of DETFs and DEIRGs in the diagnosis of endometriosis. was considered statistically significant in our study.

3. Results

3.1. Identification of DETFs and DEIRGs in Endometriosis

We obtained 1508 intersecting TFs from Human TFDB and the Human Transcription Factors database (Figure 1(a)). Subsequently, we chose gene expression datasets from the GEO datasets. GSE7305 and GSE7307 were selected to identify DEGs because both included the endometrial samples with or without endometriosis. Based on the criteria of and |log2FC| ≥1 [13], a total of 1141 DEGs from GSE7305 and GSE7307 were acquired by the GEO2R analysis tool, including 525 upregulated genes and 616 downregulated genes (Figures 1(b) and 1(c)). The DEGs were visualized by volcano plots in GSE7305 and GSE7307 (Figures 2(a) and 2(b)). Then, the results were intersected with 1508 TFs, identifying the 94 DETFs (35 upregulated DETFs and 59 downregulated DETFs) (Figures 2(c) and 2(d)). Similarly, 1793 unique IRGs were downloaded from the ImmPort database, and 111 DEIRGs (80 upregulated genes and 31 downregulated genes) were obtained from the intersection of the IRGs and the DEGs (Figures 2(e) and 2(f)). The DETFs and DEIRGs were visualized by a heatmap in GSE7305 and GSE7307, and there was a clear division between the endometriosis and control groups (Figures 3–6).

(a)

(b)

(c)

(a)

(b)

(c)

(d)

(e)

(f)

3.2. Validation of the DETFs in Different Severity Groups of Endometriosis

To verify the credibility and applicability of the DETFs, we selected the external dataset GSE51981, which contained endometriosis samples with different severities. We obtained the 40 most highly expressed DETFs (14 upregulated DETFs and 26 downregulated DETFs) in 10 randomly selected samples from the minimal/mild group and moderate/severe group in GSE7305 and GSE7307. As shown in Figure 7(a), we found that a series of upregulated DETFs were still highly expressed in the moderate/severe group, and these upregulated DETFs could predict the severity of endometriosis by receiver operating characteristic (ROC) curve analysis () (Figure 8(a)). Downregulated DETFs had a more pronounced advantage in predicting endometriosis severity than upregulated DETFs, and most of the downregulated DETFs remained expressed at low levels in the moderate/severe group () (Figure 7(b)). In addition, the ROC curve provided powerful evidence to support this view, with area under the curve (AUC) values all over 0.7(Figure 8(b)).

(a)

(b)

(a)

(b)

3.3. The Expression of DEIRGs in Women with Endometriosis with or without Infertility

The GSE120103 dataset was chosen for subsequent validation because it included infertile and fertile females with endometriosis, and we obtained the top 40 DEIRGs expressed in it. Interestingly, most of the upregulated DEIRGs were increased in infertile females with endometriosis, while some downregulated DEIRGs were expressed at low levels in the endometria of infertile women () (Figures 9(a) and 9(b)). For this result, we conducted ROC analysis to evaluate the values of DEIRGs in the diagnosis of endometriosis, and the AUC also verified the sensitivity and specificity of DEIRGs with (Figures 10(a) and 10(b)).

(a)

(b)

(a)

(b)

3.4. Functional Enrichment Analysis of DETFs and DEIRGs

To indicate the biological properties of DETFs and DEIRGs, functional analysis was performed, including GO functional and KEGG encrichment analyses. The enriched GO terms were divided into BP, CC, and MF ontologies. The GO functional enrichment results of the DETFs were mainly enriched in the BP ontology. BP analysis showed that the DETFs were significantly enriched in reproductive structure or system development, cell fate commitment, and anterior/posterior pattern specification. For the cell component, the DETFs were enriched in the transcription regulator complex and nuclear speck RNA polymerase II transcription regulator complex. The MF ontology of DETFs was mainly related to ligand-activated transcription factor activity, RNA polymerase II-specific DNA-binding transcription factor binding, and nuclear receptor binding or activity (Figure 11(a)). For the KEGG analysis, the DETFs were mainly enriched in the signaling pathways associated with transcriptional misregulation in cancer and the Notch signaling pathway (Figure 11(b)).

(a)

(b)

(c)

(d)

Likewise, the DEIRGs were also enriched in the regulation chemotaxis, lipase or phospholipase activity, cytoplasmic vesicle lumen, external side of plasma membrane, cytokine activity, and nuclear receptor activity (Figure 11(c)). The KEGG pathways of the DEIRGs were mainly enriched in viral protein interactions with cytokines and cytokine receptors, the PI3K-Akt signaling pathway, and the MAPK signaling pathway (Figure 11(d)).

3.5. PPI Network Construction and Pearson’s Correlation Analysis

Protein interactions between the DETFs and DEIRGs were constructed using the STRING online database, and the PPI network was constructed using Cytoscape. Five subnetworks were recognized. Therefore, we obtained TF-IRG regulatory networks containing 44 nodes and 73 edges to accurately illustrate the regulatory relationships between the DETFs and DEIRGs (Figure 12(a)). Pearson’s correlation analysis was used to analyze the TF-IRG regulatory network, and most nodes were correlated with other nodes at the expression level with (Figure 12(b)). The MCODE plugin of Cytoscape was used to complement the module analysis, with the corresponding modules shown in Figure 13. Furthermore, the most significantly enriched functional modules were those related to complement and coagulation cascades, Staphylococcus aureus infection, proteoglycans in cancer, focal adhesion, and the Rap1 signaling pathway (Table 2).

(a)

(b)

(c)

Subsequently, we used the cytoHubba plugin of Cytoscape to identify hub genes according to the three most important topological features in network analysis, including degree, betweenness, and closeness. We then ranked the top ten nodes for each set of the three different topological measurements (Table 3). As a result, we obtained five nodes (CXCL2, C3, VCAM1, ITGB2, and C3AR1) in all three of the lists (Figure 12(c)). These five DEIRGs can therefore be considered hub genes in the regulatory network.

3.6. Identification and Validation of Hub IRGs

As shown in GSE7305 and GSE7307, the expression of each hub IRG was significantly higher in the endometriosis group than in the control group in the box plot () (Figures 14(a) and 14(b)). The GSE23339 dataset (GPL6102 Illumina human-6 v2.0 expression beadchip) and the publicly accessible ENDOMET Turku Endometriosis Database were also used to verify the DEIRGs (Figures 14(c) and 14(d)). However, in the additional database validation of GSE23339, CXCL2 was not statistically significant. Thus, four hub genes were obtained (C3, VCAM1, ITGB2, and C3AR1). In view of the above results, we verified the above four target genes by qRT-PCR, and the highly expressed C3 and VCAM1 were statistically significant in endometriosis, and the AUCs were 0.96 and 0.76 (Figure 14(e)). However, ITGB2 and C3AR1 were not statistically significant (Figure S1). The AUCs calculated from ROC analysis in GSE7305 (Figure 15(a)), GSE7307 (Figure 15(b)), and GSE23339 (Figure 15(c)) were used to evaluate the diagnostic value of endometriosis. The AUC values of all four hub DEIRGs were over 80%, which meant that the hub DEIRGs played a critical role as novel biomarkers for endometriosis.

(a)

(b)

(c)

(d)

(e)

(a)

(b)

(c)

4. Discussion

Endometriosis is a common benign gynecological disorder characterized by immunity, inflammation, and hormone dependence. Previous studies on endometriosis have mainly focused on TFs [19] or IRGs [20], and there have been few studies on TF-related IRGs in endometriosis. TF-related IRGs not only function in immunity regulation but can also be used as prognostic biomarkers and play a key role in the development of cancer [21]. The exact molecular mechanisms of endometriosis are still unclear, and the current treatments are limited. Therefore, the discovery of new therapeutic targets and potential diagnostic biomarkers remains a research focus.

With the rapid development of the high-throughput methods and data analysis of various databases, Bohler et al. have focused on bioinformatics analysis, which can also serve as the basis for molecular biology experiments of endometriosis [22]. This study mainly analyzed DETFs and DEIRGs in endometriosis by bioinformatics methods and analyzed the expression of DETFs and DEIRGs in the endometria of women with different disease severities and infertility statuses. In addition, the enrichment analysis and networks were performed and constructed on DETFs and DEIRGs to discover valuable TFs, IRGs, and related endometriosis pathways.

Among patients with different severities of endometriosis, upregulated DETFs were more highly expressed in patients with moderate/severe endometriosis. However, downregulated DETFs presented the opposite trend. We believe that this phenomenon may be related to the degree of macrophage infiltration of different severities of endometriosis. Compared with the r-AFS stage I-II of endometriosis, the proportion of M2 macrophages in stage III-IV endometriosis was higher, suggesting that the degree of M2 macrophage infiltration was related to the severity of endometriosis [23]. M2 macrophages have anti-inflammatory properties, promote wound healing, promote fibrosis, and enable the immune escape of ectopic endometrium [24]. In addition, M2 macrophages produce matrix metalloproteinases (MMPs), such as MMP 9, which promote the ectopic growth and progression of endometrial cells by degrading the extracellular matrix and enhancing intercellular adhesion [25, 26]. At the same time, the results of DETF enrichment analysis also suggested that they may affect the activation of macrophages in endometriosis, such as the Notch pathway, which was consistent with the findings of previous research [27].

In fertile or infertile women with endometriosis, most of the upregulated DEIRGs were also highly expressed in infertile women, whereas only 4 downregulated DEIRGs were significantly expressed at a low level in the infertile group. Immunity, inflammation, and DEGs have important implications in infertile patients. These factors may affect the expression of IRGs, and that they and IRGs coregulate to influence susceptibility in patients with endometriosis-related infertility. A previous study found that BDNF (Met) single-nucleotide polymorphism, an IRG, was associated with endometriosis-related infertility women, suggesting that low levels of BDNF may be responsible for poor in vitro fertilization (IVF) outcomes in infertile patients with the BDNF (Met/Met) genotype [28]. Yin et al. studied another IRG, PTX3, which is also associated with endometriosis-related infertility [29]. In our study, IGF1 had lower expression levels in infertile women with endometriosis, and it had the ability to attenuate oocyte and embryo development resulting in endometriotic infertility, as reported in the study of Ding et al. [30]. At the same time, ESR1 showed low expression, and related studies suggested that ESR1 can affect the possibility of pregnancy in infertile patients with endometriosis [31].

In this study, the GO and KEGG enrichment analyses showed that the IRGs were mainly related to immune-related functions and pathways, such as the external side of the plasma membrane, cytokine activity and nuclear receptor activity, the PI3K-Akt signaling pathway, and the MAPK signaling pathway. In general, the external side of the plasma membrane plays an important role in endometriosis immunity. Antigens derived from the plasma membrane might directly assay reactive autoantibodies to indicate the immunoreactivity of endometriosis severity [32]. Cytokine activity, such as proinflammatory cytokines (IL-1β, IL-6) and tumor growth factor-beta (TNF-β), plays an important role in evading immune surveillance and predicting the disease severity of endometriosis [33, 34]. In addition, upregulation of MAPK subfamilies promoted the occurrence of endometriosis by influencing the function of various cytokines, including IL-6 and IL-8 [35]. The PI3K-Akt and MAPK pathways are interconnected with each other [36, 37]. The activation of the PI3K-Akt signaling pathway and the ERK-related intracellular MAPK signaling pathway was correlated with endometriosis [36], and both were shown to be involved in the immunity [38]. In addition, our research further identified the involvement of IRGs in the regulation of the PI3K-Akt signaling pathway and MAPK signaling pathway in endometriosis. A previous study showed that upregulation of the adaptor protein SHC1 had the ability to activate the PI3K-Akt and/or MAPK pathways in endometriosis samples [39]. The activation of the PI3K-Akt and MAPK pathways was associated with the immune-related pathway, nuclear factor-κB (NF-κB) signaling pathway in endometriosis cells [40]. Therefore, the study of IRGs in endometriosis is essential. These findings shed light on the screening of new potential biomarkers and the early diagnosis of endometriosis.

A previous study reported a TF-targeted gene network indicating the onset of endometriosis [41]. To further investigate the possible underlying molecular regulatory mechanisms, a TF-IRG network was constructed to study the mechanism of endometriosis. A total of 39 IRGs (RND3, PLK2, AURKA, RCAN1, EZH2, etc.) were selected to analyze the TF-IRG PPI network. Five DETFs (IRF6, EGR1, FOSB, JUNB, and MECOM) were connected with several IRGs. The genes in the PPI network were closely linked and cross regulated with each other. For instance, in the multifunctional network, IRF6/BST2 was involved in the regulation of immunity. Currently, there are many related studies on IRF6 and BST2 in immunity. Aberrant DNA methylation of IRF6 and BST2 in CD4+ T cells induced autoimmune responses [42]. Meanwhile, Figure 13 shows the predicted binding sites of IRF6 and BST2, which suggests that IRF6 and BST2 may function through mutual regulation. In a human papillomavirus type 16 (HPV16) study of host immunity, inhibition of IRF6 was responsible for immune escape from HPV16 blocking IL-1β secretion [43]. In this study, we found that BST2/CD317 in combination with TLR agonists specifically presented Ag by plasmacytoid dendritic cells in vivo, which contributed to the strong cellular and humoral immune responses [44]. However, the detailed mechanism of the PPI network should be elucidated in the future. Our findings provide an informatics basis for future research in this direction.

In this research, we mainly aimed to construct an IRG-related diagnostic model, which was established based on DEGs. ROC analysis revealed that four IRGs can be used as potential biomarkers of endometriosis, which also demonstrated the feasibility in terms of the AUC, a signal for endometriosis occurrence. Recently, C3 was considered a candidate diagnostic biomarker of endometrosis, and its expression was correlated with the engraftment of the endometriotic cysts [45]. The overexpression of VCAM-1 on the peritoneum of endometriosis had been revealed by Schutt et al. [46]. The increased expression of ITGB2 had been previously reported in endometriosis tissues compared with normal tissues [47], and high C3AR1 expression might be used as a diagnostic factor for the endometriosis-associated malignant phenotype [5]. In this research, two hub DEIRGs (C3 and VCAM1) with diagnostic value were obtained. However, this research had some limitations. First, the applicability of the diagnostic model needs to be validated in a larger sample population in future studies. Second, we will continue to complete the molecular mechanism study on the role of IRGs in endometriosis.

5. Conclusion

The TF-IRG network could be used to present novel prospective molecular mechanisms underlying the development of tumors [45]. However, studies of the regulatory mechanisms underlying TFs and IRGs in endometriosis are still in progress. In our study, IRGs were used to construct a diagnostic model to predict the onset of endometriosis patients by bioinformatics analysis. ROC analysis confirmed that the diagnostic value of hub genes (C3 and VCAM1) was clinically feasible. Additionally, the TF-IRG regulatory network broadened the horizon for research concerning the pathogenesis of endometriosis.

Data Availability

All of the data we used in this study were publicly available as described in the methods section and can be found in online Github page: https://github.com/zgm19661026/zgm19661026.git.

Ethical Approval

This research has been conducted using publicly available datasets, and no ethical approval was required. All tissue samples obtained were approved by the Ethics Commission of Harbin Medical University (202106).

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

GMZ designed all the study. YNH and JXL performed the data processing and experimental analysis. YJQ, LYS, XBZ, and HW drafted the manuscript. All authors reviewed and approved the final version of the manuscript. Yanan He and Jixin Li contributed equally to this work.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (81971359), Natural Science Foundation of Heilongjiang Province (LH2019H027), the Heilongjiang Postdoctoral Program Foundation (LBH-Z19085), and the Outstanding Young Medical Talents Training Fund project of the First Affiliated Hospital of Harbin Medical University (HYD2020YQ0021).

Supplementary Materials

Figure S1: the differences in the mRNA expression levels of ITGB2 and C3AR1 by qRT-PCR between the endometriosis tissues and the controls. (A) ITGB2. (B) C3AR1. (Supplementary Materials)

References

T. T. Tapmeier, N. Rahmioglu, J. Lin et al., “Neuropeptide S receptor 1 is a nonhormonal treatment target in endometriosis,” Science Translational Medicine, vol. 13, no. 608, 2021.
View at: Publisher Site | Google Scholar
H. S. Taylor, A. M. Kotlyar, and V. A. Flores, “Endometriosis is a chronic systemic disease: clinical challenges and novel innovations,” The Lancet, vol. 397, no. 10276, pp. 839–852, 2021.
View at: Publisher Site | Google Scholar
K. T. Zondervan, C. M. Becker, and S. A. Missmer, “Endometriosis,” The New England Journal of Medicine, vol. 382, no. 13, pp. 1244–1256, 2020.
View at: Publisher Site | Google Scholar
S. C. Lin, H. C. Lee, C. T. Hsu et al., “Targeting anthrax toxin receptor 2 ameliorates endometriosis progression,” Theranostics, vol. 9, no. 3, pp. 620–632, 2019.
View at: Publisher Site | Google Scholar
C. DeAngelo, M. B. Tarasiewicz, A. Strother et al., “Endometriosis: a malignant fingerprint,” Journal of Cancer Research and Therapeutic Oncology, vol. 8, no. 2, 2020.
View at: Google Scholar
J. E. A. Sampson, “Peritoneal endometriosis due to the menstrual dissemination of endometrial tissue into the peritoneal cavity,” American Journal of Obstetrics & Gynecology, vol. 14, no. 4, pp. 422–469, 1927.
View at: Publisher Site | Google Scholar
M. O. Perez, R. Santos, and G. S. Novaes, “Endometriosis and autoimmunity,” Revista Brasileira de Medicina, vol. 70, no. 3, pp. 66–69, 2013.
View at: Google Scholar
W. H. Fridman, L. Zitvogel, C. Saut’es-Fridman, and G. Kroemer, “The immune contexture in cancer prognosis and treatment,” Nature Reviews Clinical Oncology, vol. 14, no. 12, pp. 717–734, 2017.
View at: Publisher Site | Google Scholar
S. Shen, G. Wang, R. Zhang et al., “Development and validation of an immune gene-set based prognostic signature in ovarian cancer,” eBioMedicine, vol. 40, pp. 318–326, 2019.
View at: Publisher Site | Google Scholar
J. Vallvé-Juanico, S. Houshdaran, and L. C. Giudice, “The endometrial immune environment of women with endometriosis,” Human Reproduction Update, vol. 25, no. 5, pp. 564–591, 2019.
View at: Publisher Site | Google Scholar
P. T. K. Saunders and A. W. Horne, “Endometriosis: etiology, pathobiology, and therapeutic prospects,” Cell, vol. 184, no. 11, pp. 2807–2824, 2021.
View at: Publisher Site | Google Scholar
T. Zhang, C. De Carolis, G. C. W. Man, and C. C. Wang, “The link between immunity, autoimmunity and endometriosis: a literature update,” Autoimmunity Reviews, vol. 17, no. 10, pp. 945–955, 2018.
View at: Publisher Site | Google Scholar
H. Hu, Y. R. Miao, L. H. Jia, Q. Y. Yu, Q. Zhang, and A. Y. Guo, “AnimalTFDB 3.0: a comprehensive resource for annotation and prediction of animal transcription factors,” Nucleic Acids Research, vol. 47, no. D1, pp. D33–D38, 2019.
View at: Publisher Site | Google Scholar
S. A. Lambert, A. Jolma, L. F. Campitelli et al., “The human transcription factors,” Cell, vol. 172, no. 4, pp. 650–665, 2018.
View at: Publisher Site | Google Scholar
S. Bhattacharya, S. Andorf, L. Gomes et al., “ImmPort: disseminating data to the public for the future of immunology,” Immunologic Research, vol. 58, no. 2-3, pp. 234–239, 2014.
View at: Publisher Site | Google Scholar
Y. Xue, P. Xu, S. Xu et al., “Peptidomic analysis of endometrial tissue from patients with ovarian endometriosis,” Cellular Physiology and Biochemistry, vol. 47, no. 1, pp. 107–118, 2018.
View at: Publisher Site | Google Scholar
G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
X. Robin, N. Turck, A. Hainard et al., “pROC: an open-source package for R and S+ to analyze and compare ROC curves,” BMC Bioinformatics, vol. 12, no. 1, p. 77, 2011.
View at: Publisher Site | Google Scholar
D. Mathew, J. A. Drury, A. J. Valentijn, O. Vasieva, and D. K. Hapangama, “In silico, in vitro and in vivo analysis identifies a potential role for steroid hormone regulation of FOXD3 in endometriosis-associated genes,” Human Reproduction, vol. 31, no. 2, p. dev307, 2015.
View at: Google Scholar
H. C. Bohler, C. Gercel-Taylor, B. A. Lessey, and D. D. Taylor, “Endometriosis markers: immunologic alterations as diagnostic indicators for endometriosis,” Reproductive Sciences, vol. 14, no. 6, pp. 595–604, 2007.
View at: Publisher Site | Google Scholar
Y. L. Sun, Y. Zhang, Y. C. Guo, Z. H. Yang, and Y. C. Xu, “A prognostic model based on the immune-related genes in colon adenocarcinoma,” International Journal of Medical Sciences, vol. 17, no. 13, pp. 1879–1896, 2020.
View at: Publisher Site | Google Scholar
V. Singh, S. H. Ahn, and C. Tayade, “Biomarkers in endometriosis: challenges and opportunities,” Fertility and Sterility, vol. 107, no. 3, 2017.
View at: Publisher Site | Google Scholar
Q. Zhong, F. Yang, X. Chen, J. Li, C. Zhong, and S. Chen, “Patterns of immune infiltration in endometriosis and their relationship to r-AFS stages,” Frontiers in Genetics, vol. 12, article 631715, 2021.
View at: Publisher Site | Google Scholar
M. Hesketh, K. B. Sahin, Z. E. West, and R. Z. Murray, “Macrophage phenotypes regulate scar formation and chronic wound healing,” International Journal of Molecular Sciences, vol. 18, no. 7, p. 1545, 2017.
View at: Publisher Site | Google Scholar
C. Atri, F. Z. Guerfali, and D. Laouini, “Role of human macrophage polarization in inflammation during infectious diseases,” International Journal of Molecular Sciences, vol. 19, no. 6, p. 1801, 2018.
View at: Publisher Site | Google Scholar
H. Liu, J. Wang, H. Wang et al., “Correlation between matrix metalloproteinase-9 and endometriosis,” International Journal of Clinical and Experimental Pathology, vol. 8, no. 10, pp. 13399–13404, 2015.
View at: Google Scholar
W. Chen, Y. Liu, J. Chen et al., “The notch signaling pathway regulates macrophage polarization in liver diseases,” International Immunopharmacology, vol. 99, article 107938, 2021.
View at: Publisher Site | Google Scholar
Q. Y. Zhang, Q. Guan, Y. Wang et al., “BDNF Val66Met polymorphism is associated with stage III-IV endometriosis and poor in vitro fertilization outcome,” Human Reproduction, vol. 27, no. 6, pp. 1668–1675, 2012.
View at: Publisher Site | Google Scholar
Y. Yin, Y. Mao, A. Liu et al., “Insufficient cumulus expansion and poor oocyte retrieval in endometriosis-related infertile women,” Reproductive Sciences, vol. 28, no. 5, pp. 1412–1420, 2021.
View at: Publisher Site | Google Scholar
G.-L. Ding, X.-J. Chen, Q. Luo, M.-Y. Dong, N. Wang, and H.-F. Huang, “Attenuated oocyte fertilization and embryo development associated with altered growth factor/signal transduction induced by endometriotic peritoneal fluid,” Fertility and Sterility, vol. 93, no. 8, pp. 2538–2544, 2010.
View at: Publisher Site | Google Scholar
M. Lamp, M. Peters, E. Reinmaa et al., “Polymorphisms in ESR1, ESR2 and HSD17B1 genes are associated with fertility status in endometriosis,” Gynecological Endocrinology, vol. 27, no. 6, pp. 425–433, 2011.
View at: Publisher Site | Google Scholar
H. Ho, M. Wu, and Y. Yang, “Peritoneal cellular immunity and endometriosis,” American Journal of Reproductive Immunology, vol. 38, no. 6, pp. 400–412, 2011.
View at: Google Scholar
D. I. Sokolov, N. G. Solodovnikova, O. V. Pavlov, D. A. Niauri, N. N. Volkov, and S. A. Sel'kov, “Study of cytokine profile and angiogenic potential of peritoneal fluid in patients with external genital endometriosis,” Bulletin of Experimental Biology and Medicine, vol. 140, no. 5, pp. 541–544, 2005.
View at: Publisher Site | Google Scholar
A. Pizzo, F. M. Salmeri, F. V. Ardita, V. Sofo, M. Tripepi, and S. Marsico, “Behaviour of cytokine levels in serum and peritoneal fluid of women with endometriosis,” Gynecologic and Obstetric Investigation, vol. 54, no. 2, pp. 82–87, 2003.
View at: Publisher Site | Google Scholar
O. Yoshino, Y. Osuga, Y. Hirota et al., “Possible pathophysiological roles of mitogen-activated protein kinases (MAPKs) in endometriosis,” American Journal of Reproductive Immunology, vol. 52, no. 5, pp. 306–311, 2004.
View at: Publisher Site | Google Scholar
G. Bora and A. Yaba, “The role of mitogen-activated protein kinase signaling pathway in endometriosis,” The Journal of Obstetrics and Gynaecology Research, vol. 47, no. 5, pp. 1610–1623, 2021.
View at: Publisher Site | Google Scholar
A. Makker, M. M. Goel, V. Das, and A. Agarwal, “PI3K-Akt-mTOR and MAPK signaling pathways in polycystic ovarian syndrome, uterine leiomyomas and endometriosis: an update,” Gynecological Endocrinology, vol. 28, no. 3, pp. 175–181, 2012.
View at: Publisher Site | Google Scholar
L. Wang, X. Zhang, G. Wu et al., “Streptococcus pneumoniae aminopeptidase N contributes to bacterial virulence and elicits a strong innate immune response through MAPK and PI3K/AKT signaling,” Journal of Microbiology, vol. 58, no. 4, pp. 330–339, 2020.
View at: Publisher Site | Google Scholar
H. Honda, F. F. Barrueto, J. Gogusev, D. D. Im, and P. J. Morin, “Serial analysis of gene expression reveals differential expression between endometriosis and normal endometrium. Possible roles for AXL and SHC1 in the pathogenesis of endometriosis,” Reproductive Biology and Endocrinology, vol. 6, no. 1, p. 59, 2008.
View at: Publisher Site | Google Scholar
H. Zhang, X. Zhao, S. Liu, J. Li, Z. Wen, and M. Li, “17βE₂ promotes cell proliferation in endometriosis by decreasing PTEN via NFκB-dependent pathway,” Molecular and Cellular Endocrinology, vol. 317, no. 1-2, pp. 31–43, 2010.
View at: Publisher Site | Google Scholar
S. Cong, Q. Guo, Y. Cheng et al., “Identification and analyzation of differentially expressed transcription factors in endometriosis,” Frontiers in Molecular Biosciences, vol. 7, article 614427, 2021.
View at: Publisher Site | Google Scholar
K. A. Koelsch, R. Webb, M. Jeffries et al., “Functional characterization of the MECP2/IRAK1 lupus risk haplotype in human T cells and a human MECP2 transgenic mouse,” Journal of Autoimmunity, vol. 41, pp. 168–174, 2013.
View at: Publisher Site | Google Scholar
M. Ainouze, P. Rochefort, P. Parroche et al., “Human papillomavirus type 16 antagonizes IRF6 regulation of IL-1β,” PLoS Pathogens, vol. 14, no. 8, article e1007158, 2018.
View at: Publisher Site | Google Scholar
J. Loschko, A. Schlitzer, D. Dudziak et al., “Antigen delivery to plasmacytoid dendritic cells via BST2 induces protective T cell-mediated immunity,” Journal of Immunology, vol. 186, no. 12, pp. 6718–6725, 2011.
View at: Publisher Site | Google Scholar
C. Agostinis, S. Zorzet, A. Balduit et al., “The inflammatory feed-forward loop triggered by the complement component C3 as a potential target in endometriosis,” Frontiers in Immunology, vol. 12, article 693118, 2021.
View at: Publisher Site | Google Scholar
A. K. Schutt, K. A. Atkins, J. K. Slack-Davis, and D. W. Stovall, “VCAM-1 on peritoneum and α4β1 integrin in endometrium and their implications in endometriosis,” International Journal of Gynecological Pathology, vol. 34, no. 1, pp. 85–89, 2015.
View at: Publisher Site | Google Scholar
J. Sundqvist, K. L. Andersson, G. Scarselli, K. Gemzell-Danielsson, and P. G. L. Lalitkumar, “Expression of adhesion, attachment and invasion markers in eutopic and ectopic endometrium: a link to the aetiology of endometriosis,” Human Reproduction, vol. 27, no. 9, pp. 2737–2746, 2012.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2023 Yanan He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

840

Downloads

718

Citations