- 1Reproductive Medicine Center, The Sixth Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
- 2Department of Gynecology and Obstetrics, The Sixth Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
Background: Endometriosis negatively affects fertility, and it is a common disease in assisted reproductive practice. Surgical removal of endometriotic lesions is widely carried out to relieve symptoms and promote fertility. But it is not intensively investigated what changes in the secretory eutopic endometrium of patients with endometriosis after surgery.
Methods: Eighteen patients with stage III/IV endometriosis were included in the study, and they were divided into the untreated group and the treated group (6 vs. 12). Basic clinical data were compared, and transcriptomic data of the secretory eutopic endometrium were analyzed with DESeq2, Cytoscape, ClueGO, CluePedia, and Gene Set Enrichment Analysis (GSEA). CIBERSORT was used to calculate the relative abundance of 22 immune cells in the samples.
Results: We determined 346 differentially expressed genes (DEGs) using DESeq2. These DEGs were used to enrich seven Gene Ontology terms including three associated with immune processes and one correlated to prostaglandin using ClueGO and CluePedia. GSEA enriched 28 Gene Ontology terms in the treated group mainly associated with immune and blood pressure regulation process. Compared to the untreated group, the relative abundance of resting CD4+ memory T cells [0.218 (0.069, 0.334) vs. 0.332 (0.181, 0.429), P = 0.022] and the even less abundant memory B cells [0.001 (0.000, 0.083) vs. 0.033 (0.007, 0.057), P = 0.049] are significantly decreased in the treated group.
Conclusion: Surgical treatment of stage III/IV endometriosis influences some genes and biological processes related to endometrial receptivity, but more evidence is needed.
Introduction
Endometriosis is recognized as a chronic estrogen-dependent complex syndrome and inflammatory disease that affects around 5%~10% of women, characterized by the abnormally existing endometrium-like tissues outside the uterus, commonly located on pelvic organs and tissues (1). Endometriosis significantly disturbs female fertility partly due to abnormal endometrial receptivity since studies in human and animal models show defective implantation (2).
Endometrial receptivity, coined as the ability of the endometrium to allow normal implantation, is fundamental for a successful pregnancy and is affected in patients with endometriosis. Multi-omics were applied to identify potential biomarkers for assessing endometrial receptivity. Clinical treatment of endometriosis improves pregnancy rates, and basic research showed abnormal endometrial milieu. For instance, prokineticin-1 (PROK1), which is important in the vascular function of peri-implantation endometrium and early pregnancy, is significantly less expressed in the eutopic endometrium of patients with endometriosis (3). Other factors that are associated with decidualization and implantation such as L-selectin ligand, ανβ3 integrin, and leukemia inhibitory factor (LIF) are also decreased in the eutopic endometrium (4–6).
Surgery that removes ectopic endometrial tissue is a common measure to relieve symptoms and promote fertility. Some evidence shows that fertility is improved after surgery (7, 8), but the changes in the eutopic endometrium after surgery are not well-illustrated. We hereby conducted this research to preliminarily discover the effect of surgery on the secretory eutopic endometrium.
Materials and methods
Criteria of patients included in the research groups
Six patients in the untreated group were diagnosed with endometriosis by ultrasound or magnetic resonance imaging (MRI), as endometrioma exists and did not receive surgical treatment for endometriosis. Twelve patients in the treated group were confirmed to have endometriosis during laparoscopic or transabdominal surgery in less than 24 months before endometrial sampling, and the endometriotic lesions were removed during surgery. All 18 patients were of stage III/IV according to the revised American Fertility Society (r-AFS) classification. The exclusion criteria are as follows: 1) follicle stimulation hormone (FSH) >12 mIU/ml and anti-Müllerian hormone (AMH) <1.1 ng/ml; 2) abnormal uterine cavity morphology by ultrasound, MRI, or hysteroscopy; 3) signs of malignancies.
Endometrial preparation and sample collection protocol
Patients received 2 mg of estradiol valerate (Progynova, Bayer Pharma AG, Berlin, Germany) twice daily for 28 days starting on the second day of menstrual bleeding. From the 15th day onward, 10 mg of dydrogesterone (Duphaston, Abbott Biologicals B.V., Netherlands) twice daily was applied.
Endometrial sampling was conducted on P+5 (the 5th day of Duphaston administration) using a sterile endometrium sampler kit with a separate package (Type I, Run Ting). The patient was kept in the lithotomy position and received vulvar and vaginal disinfection with the introduction of the disinfected speculum, and a sterile swab was inserted into the cervical canal to avoid contamination. Then, a disposable sterile endometrial sampler (Yikon Inc.) was introduced into the uterine cavity to suck the endometrium. After that, the external surface of the catheter was cleaned with a sterile gauze and half of its content was transferred to cryotubes containing 1.5 ml of tissue preservation solution (XK-039-3, Yikon Inc.) at -20°C for future RNA extraction and PCR test, and the other half was stored in formalin for hematoxylin and eosin staining and immunohistochemistry staining for CD138. Chronic endometritis (CE) is determined as at least one plasma cell is found in the endometrial biopsy using immunohistochemistry staining for CD138 (9, 10).
RNA extraction and RNA sequencing procession
Total RNA was extracted using RNeasy Micro Kit (74004, Qiagen) according to the manufacturer’s instructions. Then, RNA quantitative detection was conducted with Qubit RNA HS Kit (Q32855; Thermo Fisher Scientific). Agilent Bioanalyzer 2100 (Agilent) was used to check the integrity of the extracted total RNA. Only samples whose RNA integrity number (RIN) is greater than 7 were considered qualified samples and used for subsequent testing.
Next, RNA reverse transcription and amplification were conducted using MALBAC® Platinum Single Cell RNA Amplification Kit (KT110700796; Xukang Co., Ltd.). The positive and negative controls were 500 ng of high-quality host total RNA and ultrapure water, respectively. After that, 1 µl cDNA was 10-fold diluted and detected with Agilent Bioanalyzer 2100 (Agilent).
Finally, the library was constructed using gene sequencing and library preparation kit (XK-038, Xukang Co., Ltd.) according to the manufacturer’s instructions. Following purification, the library was quantified with Qubit dsDNA HS kit (Q32584, Invitrogen). According to the quantitative results of Qubit, each sample was taken from 10 ng library and mixed in equal proportions, followed by Qubit quantitative detection again. Pair-end sequencing of the products was performed on the NextSeq CN500 platform (Illumina).
RNA extraction, quantitative real-time PCR, and public database search
Total RNA was extracted from five samples from the untreated group and 12 samples from the treated group using RNeasy Micro Kit (Qiagen, MD, USA; Cat# 74004) according to manufacturer’s instructions. Reverse transcription of first-strand complementary DNAs (cDNAs) was carried out using HiScript III RT SuperMix for qPCR (Vazyme, Nanjing, China, Cat# R323-01). Quantitative PCR (qPCR) was conducted using 2× RealStar Green Power Mixture (Genstar; Cat# A311-101) on a Roche LightCycler 480 II (Roche Diagnostics, Mannheim, Germany). The expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was chosen as an endogenous control. The primers used were as follows:
HOXB2: forward: GATGAAAGAGAAGAAATCCGCC,
reverse: AAGTGGAATTCCTTCTCCAGTT;
GAPDH: forward: GGTCGGAGTCAACGGATTT,
reverse: CCAGCATCGCCCCACTTG.
The Human Protein Atlas (11) is used for the HOXB2 gene expression feature identification.
RNA sequencing analysis procession
Firstly, FastQC v.0.11.9 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) was used to evaluate fastq files, and then the files were processed with trim_galore (http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/) to remove reads containing adapters, more than 10% unknown nucleotides (N), and more than 50% of low-quality (Q-value ≤20) bases. Next, HISAT2 (12) was used to map the processed paired-end reads to the human genome GRCh38 from Gencode v26. FeatureCount (13) was used to generate the count value of gene expression in the form of a matrix. After that, principal component analysis was done with TPM (Transcripts Per Kilobase of exon model per Million mapped reads) value normalized from the matrix. The differentially expressed genes (DEGs) between the two groups were selected by using DESeq2 R package (14) with the P value <0.05 and log fold change (logFC) >1 as the cutoff. Then, functional enrichment analysis of DEGs was performed with ClueGO v.2.5.8 (15) and CluePedia 1.5.8 (16) plugins of Cytoscape v.3.9.1 (17). GSEA 4.2.1 (http://www.gsea-msigdb.org/gsea/index.jsp) was used to conduct functional enrichment analysis of whole gene expression features.
Calculation of endometrial immune cell abundance
CIBERSORT (18) was used to estimate the abundance of 22 types of immune cells in the endometrium.
Statistical analysis
Clinical characteristics of participants, including age, body mass index (BMI), FSH, and AMH, are presented as the mean ± standard deviation (SD) and were compared between groups using the t-test for equality of means. The relative abundance of immune cells is presented as the median (range) and was compared between groups using the nonparametric Mann–Whitney U-test. CE rate is presented as n% (positive/all) and was compared between groups using the Fisher’s exact test. Above statistical analysis is conducted with SPSS 22.0 software (IBM Corporation, Armonk, NY, USA). The relative mRNA expression ratio was compared between groups using the nonparametric Mann–Whitney U-test by GraphPad Prism 9.0.0 software (Bethesda, MD, USA). Statistical significance was set at P value <0.05.
Results
Clinical features of enrolled patients
The age, BMI, FSH, AMH, and CE rates are demonstrated in Table 1.
Transcriptomic features of the secretory eutopic endometrium in between the untreated and treated groups
We first conducted PCA that shows relatively obvious discrimination between the two groups (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/932339/fendo-13-932339-HTML-r2/image_m/fendo-13-932339-g001.jpg)
Figure 1 PCA of the transcriptomic data. The raw expression matrix was normalized to TPM before PCA. PCA, principal component analysis; TPM, Transcripts Per Kilobase of exon model per Million mapped reads; Untreated, the untreated group; Treated, the treated group.
Then, we used DESeq2 to determine the DEGs between the groups. A total of 346 DEGs were filtered with logFC >|1| and P value <0.05. Among these DEGs, 140 and 206 genes are upregulated and downregulated in the treated group, respectively. The top 10 upregulated and downregulated genes with high expression levels are exhibited in Table 2. All 346 DEGs are presented in Supplementary Material 1. Of the 20 top expressed DEGs, we found HOXB2 to be of importance, since it is a downstream gene of type I interferon (IFN) response and we validated its expression by qPCR (Figure 2A). We then searched for the expression feature of HOXB2 in the human protein atlas database, and it turns out that HOXB2 expresses highly in endometrial ciliated cells, natural killer (NK) cells, and T cells (Figure 2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/932339/fendo-13-932339-HTML-r2/image_m/fendo-13-932339-g002.jpg)
Figure 2 The expression of HOXB2. (A) Relative mRNA expression of HOXB2 in the treated (T) and untreated (UT) groups by quantitative real-time PCR. **P value = 0.006. (B) Top 15 cell types with high expression levels of HOXB2. Raw data are from the human protein atlas (www.proteinatlas.org). Untreated, the untreated group; Treated, the treated group.
Functional analysis of the transcriptomic features
Functional analysis derived from ClueGO suggests that “regulation of αβT cell proliferation,” “regulation of natural killer cell mediated immunity,” “regulation of defense response to bacterium,” “regulation of response to wounding,” “epithelial fluid transport,” “nucleoside biphosphate biosynthetic process,” and “prostaglandin transport” are different between the untreated and treated groups (Figure 3A). Of the seven Gene Ontology (GO) terms, three are associated with immune processes. Gene counts, gene ratio, and Bonferroni-corrected P value of each term are shown in Figure 3B.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/932339/fendo-13-932339-HTML-r2/image_m/fendo-13-932339-g003.jpg)
Figure 3 GO enrichment by ClueGO and CluePedia. (A) Network of biological process enriched. (B) Gene counts, gene ratio, and P value after Bonferroni correction of each GO: biological process terms. GO, gene ontology.
To have a broader view of the potential changes based on the whole transcriptomic features, we ran GSEA and enriched 28 GO terms in the treated group, including 17 terms of biological process, 8 terms of molecular function, and 3 terms of cellular component (Table 3). Of these GO terms, “defense response to Gram negative bacterium,” “antibacterial humoral response,” “humoral immune response mediated by circulating immunoglobulin”, “immunological memory process”, “immunoglobulin receptor binding”, “antigen binding”, “peptide antigen binding”, “immunoglobulin complex”, “T cell receptor complex,” and “MHC protein complex” correlate with immune process.
Several terms associated with blood pressure regulation such as “regulation of systemic arterial blood pressure by renin angiotensin”, “serotonin receptor signaling pathway”, “regulation of systemic arterial blood pressure by hormone”, “regulation of systemic arterial blood pressure mediated by a chemical signal”, “regulation of systemic arterial blood pressure by circulatory renin angiotensin”, “vasodilation”, “regulation of systemic arterial blood pressure”, and “serotonin receptor activity” are also enriched.
Calculated immune cell abundance in the secretory eutopic endometrium in between the untreated and treated groups
To further investigate the immunological changes in the endometrium, we used CIBERSORT to estimate the relative percent of 22 types of immune cells. Activated NK cells and resting CD4+ memory T cells comprise the majority of the immune cell population in the secretory eutopic endometrium (Figure 4A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/932339/fendo-13-932339-HTML-r2/image_m/fendo-13-932339-g004.jpg)
Figure 4 Immune cell abundance and comparison of certain immune cells between the groups. (A) CIBERSORT-calculated relative abundance of 22 different immune cells in the 18 samples. (B) Comparison of memory B cells, resting CD4+ memory T cells, and activated NK cells in the treated and untreated groups. *P value <0.05. ns, not significant. Untreated, the untreated group; Treated, the treated group.
Then, the immune cell populations were compared between the groups, three of them are shown in Figure 4B. Activated NK cells are not significantly different [0.424 (0.248, 0.526) in the treated group vs. 0.365 (0.198, 0.568) in the untreated group, P = 0.574]. While the less abundant resting CD4+ memory T cells and the even less abundant memory B cells are significantly decreased in the treated group [0.218 (0.069, 0.334) vs. 0.332 (0.181, 0.429), P = 0.022, and 0.001 (0.000, 0.083) vs. 0.033 (0.007, 0.057), P = 0.049, respectively]. The comparison of all of the immune cells is concluded in Supplementary Material 2.
Discussion
In this study, we applied RNA-seq analysis and CIBERSORT to profile the endometrial environment during the secretory phase after surgical treatment of endometriosis. We found several DEGs and biological processes that could be associated with receptivity. Also, altered function and abundance of certain immune cells are speculated to be related to the possibly differed endometrial microbiota after surgery.
Changes identified at the gene expression level
The expression of certain genes associated with endometriosis and receptivity is altered after surgery.
HOXB2, a relatively less studied member of the homeobox genes in the endometrium is upregulated in type I IFN response (19), and one study pointed out that the response to IFN signaling is dysregulated in the eutopic endometrium of patients with adenomyosis during the secretory phase (20). IFNs bind to specific membrane receptors to exert biological activities. Upregulated transcription of IFN-α in the human endometrium during the window of implantation (21) and the effect of type I IFN-τ on receptivity and implantation in the endometrium of the ruminant animal (22) indicate the importance of IFNs in fertility. IFNs also play an important role in immune modulation, including the adaptive immune response (23). As a downstream target gene of type I IFN, the upregulated HOXB2 expression is a putative implication of the altered IFN signaling in the endometrium of patients with endometriosis after surgery. CREB3L1 that is identified as a target gene of the progesterone receptor is decreased in the eutopic endometrium of patients with endometriosis, and it regulates the phosphorylation of extracellular signal-regulated kinase (ERK)1/2 in the decidualization process (24). The transcription level of CLDN4 is significantly higher in the endometrium of idiopathic infertility and minimal endometriosis, indicating that CLDN4 overexpression might be negatively linked to endometrial receptivity (25). STC1 serves in many processes (26), and it correlates with receptivity markers during the window of implantation (27). In patients with endometriosis, STC1 is upregulated in the mid-secretary eutopic endometrium (28).
In our results, CREB3L1 and HOXB2 are upregulated after surgery, while CLDN4 and STC1 are downregulated. These alterations might be important in remodeling receptivity.
Changes identified at the biological process level
Several biological processes are interconnected and found to be altered after surgery.
Prostaglandins (PGs) are bioactive in tissues and organs. PGE2 is observed to increase during the luteal phase in the animal uterus, and this hints that it might be associated with luteal function and implantation. PGE2 induces LH receptor expression on the corpus luteum through receptor EP2 (29, 30), and it increases the blood supply of the uterus and ovary by increasing a vasodilator called nitric oxide via receptor EP4 (31). In a recent review, the authors hypothesized that PGE2 prevents luteolysis by playing a role in the P4 secretion stimulated by estradiol (E2) and IFN-τ, and the latter is a member of type I IFNs (32). In goat endometrial stromal cells, it has been validated that IFN-τ could increase the ratio of PGE2/PGF2α via JAB1 and the unfolded protein response (UPR) mainly by regulating the amount of PGE2 rather than PGF2α (33). As the increased expression of HOXB2 hints, this process is likely to be the underlying mechanism that alters luteal function.
PGE2 is also found to regulate the CXCR4 expression through the epidermal growth factor receptor (EGFR)-phosphatidyl inositol-3 kinase (PI3K) and ERK1/2 pathway via prostaglandin E receptor 2 (PTGER2) in the endometrium (34–37). One study confirmed that CXCR4 expression is especially increased at the apposition site (38). CXCR4 is vital in embryonic vasculogenesis (39), which might be crucial in placenta attachment. Furthermore, the PI3K/ERK1/2 pathways enhance the growth, proliferation, differentiation, and survival of endometrial cells and embryos (40). PGE2 could also help increase the expression of ανβ3 integrin, thus helping embryo adhesion (41). Fetal–maternal crosstalk between the embryo and the endometrium is essential in successful pregnancy, and PGE2 is widely involved in this process. Our analysis suggests that “prostaglandin transport” is altered after surgery, which might contribute to the altered response of the endometrium to the embryo and luteal function, as PGE2 might function in implantation by means that are mentioned above.
Changes identified at the immune environment level
The endometrium contains different kinds of immune cells and relative molecules, which are important during implantation and pregnancy.
A PGE2–myocyte enhancer factor 2A (MEF2A) axis in type I IFN induction is introduced recently. Research shows that PGE2 interferes with lipopolysaccharide (LPS)-mediated activation of ERK5 that is a known transcriptional partner of MEF2 and thus affects inflammatory gene expression (42). Given the previously mentioned regulation of IFN-τ on PGE2, it can be inferred that the PG and IFN response in the endometrium interact mutually and they maintain balance in an accurate manner, and surgery might influence that regulation process.
Implantation correlates with innate immune cells in the endometrium. Although our analysis did not show significant alteration of several innate immune cells in relative quantity after surgery, “regulation of αβT cell proliferation” and “regulation of natural killer cell mediated immunity” are enriched in the functional enrichment. In patients with endometriosis who obtained successful implantation, uterine NK progenitor cell populations are markedly higher than those in patients who have failed implantation (43). Uterine NK cells secrete angiogenic growth factor, contributing to decidualization and formation of spiral arteries, which makes them important in early pregnancy (44). However, there is evidence that NK cells are not essential for pregnancy (45). The ratio of T helper type 1 (Th1) cytokines over T helper type 2 (Th2) cytokines is related to fetal acceptance. Inflammatory Th1 immunity is dominant at the peri-implantation stage, it could benefit the invading of trophoblasts (46). Interestingly, there is evidence that type I IFN response works in the regulation of NK cells and T cells. In direct and indirect activation of NK cells, type I IFN could activate NK cells through the IFN-activating receptor on NK cells or dendritic cells are first activated by type I IFN, then they activate NK cells by trans-presentation of interleukin (IL)-15 to IL-15 receptors on NK cells, respectively (47). As for T cells, type I IFN could also alter the αβT-cell population by inducing an immune response of dendritic cells via both IL-15 and the IL-15R α-chain (48). Exposure to IFN-α/β results in increased expression of genes associated with cell proliferation and cell survival (49), which could be the mechanism of αβT-cell proliferation. Type I IFN leads to Th1 induction and Th2 restriction, which might contribute to better implantation (50). Due to the lack of adequate subjects, the function of T cells and NK cells after surgery needs to be further elucidated.
Implantation might correlate with adaptive immune cells in the endometrium as well. Memory B cells and resting CD4+ memory T cells are less in the treated group. We also find “antibacterial humoral response,” “humoral immune response mediated by circulating immunoglobulin,” and “immunological memory process” are enriched in the treated group. Higher IL-2 seems to be adverse to pregnancy, since it induces increased NK activity of decidual NK cells (51) that leads to a less tolerant endometrial environment, and application of IL-2 causes fetal development inhibition in mice (52). Surgical treatment could reduce serum IL-2 (53), which might facilitate receptivity. CD4+ memory T cells are typical IL-2-producing cells (54), thus their reduction in the relative population might facilitate a healthier milieu for implantation. Both CD4+ memory T cells and memory B cells are generated in the adaptive immune process, so their abundance could be closely related to the microbes in the uterine cavity. Endometrial microbiota seems to be more diverse in patients with endometriosis (55). One study found a complete absence of Atopobium in the vaginal and cervical microbiota and an increase of potentially pathogenic Gardnerella, Streptococcus, Escherichia, Shigella, and Ureoplasma in the cervical microbiota in III/IV endometriosis (56). Another study categorized the endometrial microbial composition as either Lactobacillus-dominant (LD, >90% Lactobacillus spp.) or non-Lactobacillus-dominant (NLD, <90% Lactobacillus spp. with >10% of other bacteria), and NLD microbiotas are associated with adverse reproductive outcomes (57). Since the female genital tract is consecutively unobstructed and relatively short, the endometrial microbiota possibly has a crosstalk with that in the lower tract. Clinical administration of a broad-spectrum antibiotic significantly decreased several kinds of microbiomes (58). The endometrial microbiota might alter after surgery, since antibiotics were applied during the surgery and thus caused altered adaptive immune cell proportions. Based on these clues, it could be inferred that adaptive immune cells might be a potential marker of endometrial microbiotic milieu due to their interaction with the microbiota, and the microbiota could interfere with endometrial receptivity indirectly by affecting the adaptive immune cells.
Upon these findings and speculation, we hypothesized that surgical treatment might modify immune components and the endometrial microbiota, leading to the alteration of innate and adaptive immune cells. Further investigation of immune factors and endometrial microbiota after surgery is needed to prove this.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: National Genomics Data Center, accession number HRA002724 (https://ngdc.cncb.ac.cn/gsa-human/browse/HRA002724).
Ethics statement
This study was reviewed and approved by The Ethics Committee of the Sixth Affiliated Hospital of Sun Yat-Sen University. All methods were carried out in accordance with the Declaration of Helsinki and the ethical guidelines and regulations of the Ethics Committee of the Sixth Affiliated Hospital of Sun Yat-Sen University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
RX, JP, and HZ proposed the idea and designed the research. ZZ, HL, JZ, JP, and HZ collected clinical data and tissue samples of the patients. RX, PC, ZZ, and CZ analyzed and interpreted the transcriptomic data. RX drafted the manuscript. JP and HZ did the revision. All authors read and approved the final manuscript.
Funding
This work was supported by the National Key Research and Development Program of China (grant no. 2021YFC2700503).
Acknowledgments
We appreciate Xukang Medical·Yikon Genomics for helping with the sequencing process.
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/fendo.2022.932339/full#supplementary-material
References
1. Taylor HS, Kotlyar AM, Flores VA. Endometriosis is a chronic systemic disease: clinical challenges and novel innovations. Lancet (2021) 397:839–52. doi: 10.1016/S0140-6736(21)00389-5
2. Lessey BA, Kim JJ. Endometrial receptivity in the eutopic endometrium of women with endometriosis: it is affected, and let me show you why. Fertil Steril (2017) 108:19–27. doi: 10.1016/j.fertnstert.2017.05.031
3. Tiberi F, Tropea A, Apa R, Romani F, Lanzone A, Marana R. Prokineticin 1 mRNA expression in the endometrium of healthy women and in the eutopic endometrium of women with endometriosis. Fertil Steril (2010) 93:2145–9. doi: 10.1016/j.fertnstert.2009.01.105
4. Lessey BA, Castelbaum AJ, Sawin SW, Buck CA, Schinnar R, Bilker W, et al. Aberrant integrin expression in the endometrium of women with endometriosis. J Clin Endocrinol Metab (1994) 79:643–9. doi: 10.1210/jcem.79.2.7519194
5. Margarit L, Gonzalez D, Lewis PD, Hopkins L, Davies C, Conlan RS, et al. L-selectin ligands in human endometrium: comparison of fertile and infertile subjects. Hum Reprod (2009) 24:2767–77. doi: 10.1093/humrep/dep247
6. Franasiak JM, Holoch KJ, Yuan L, Schammel DP, Young SL, Lessey BA. Prospective assessment of midsecretory endometrial leukemia inhibitor factor expression versus ανβ3 testing in women with unexplained infertility. Fertil Steril (2014) 101:1724–31. doi: 10.1016/j.fertnstert.2014.02.027
7. Zhong C, Gao L, Shu L, Hou Z, Cai L, Huang J, et al. Analysis of IVF/ICSI outcomes in endometriosis patients with recurrent implantation failure: Influence on cumulative live birth rate. Front Endocrinol (Lausanne) (2021) 12:640288. doi: 10.3389/fendo.2021.640288
8. Bafort C, Beebeejaun Y, Tomassetti C, Bosteels J, Duffy JM. Laparoscopic surgery for endometriosis. Cochrane Database Syst Rev (2020) 10:CD011031. doi: 10.1002/14651858.CD011031.pub3
9. Bayer-Garner IB, Nickell JA, Korourian S. Routine syndecan-1 immunohistochemistry aids in the diagnosis of chronic endometritis. Arch Pathol Lab Med (2004) 128:1000–3. doi: 10.5858/2004-128-1000-RSIAIT
10. Kitaya K, Takeuchi T, Mizuta S, Matsubayashi H, Ishikawa T. Endometritis: new time, new concepts. Fertil Steril (2018) 110:344–50. doi: 10.1016/j.fertnstert.2018.04.012
11. Pontén F, Jirström K, Uhlen M. The human protein atlas–a tool for pathology. J Pathol (2008) 216:387–93. doi: 10.1002/path.2440
12. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods (2015) 12:357–60. doi: 10.1038/nmeth.3317
13. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics (2014) 30:923–30. doi: 10.1093/bioinformatics/btt656
14. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol (2014) 15:550. doi: 10.1186/s13059-014-0550-8
15. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics (2009) 25:1091–3. doi: 10.1093/bioinformatics/btp101
16. Bindea G, Galon J, Mlecnik B. CluePedia cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics (2013) 29:661–3. doi: 10.1093/bioinformatics/btt019
17. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res (2003) 13:2498–504. doi: 10.1101/gr.1239303
18. 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:453–7. doi: 10.1038/nmeth.3337
19. Vestergaard AL, Knudsen UB, Munk T, Rosbach H, Martensen PM. Transcriptional expression of type-I interferon response genes and stability of housekeeping genes in the human endometrium and endometriosis. Mol Hum Reprod (2011) 17:243–54. doi: 10.1093/molehr/gaq100
20. Prašnikar E, Kunej T, Gorenjak M, Potočnik U, Kovačič B, Knez J. Transcriptomics of receptive endometrium in women with sonographic features of adenomyosis. Reprod Biol Endocrinol (2022) 20:2. doi: 10.1186/s12958-021-00871-5
21. Li Q, Zhang M, Kumar S, Zhu LJ, Chen D, Bagchi MK, et al. Identification and implantation stage-specific expression of an interferon-alpha-regulated gene in human and rat endometrium. Endocrinology (2001) 142:2390–400. doi: 10.1210/endo.142.6.8101
22. Roberts RM. Interferon-tau, a type 1 interferon involved in maternal recognition of pregnancy. Cytokine Growth Factor Rev (2007) 18:403–8. doi: 10.1016/j.cytogfr.2007.06.010
23. Zivcec M, Spiropoulou CF, Spengler JR. The use of mice lacking type I or both type I and type II interferon responses in research on hemorrhagic fever viruses. part 2: Vaccine efficacy studies. Antiviral Res (2020) 174:104702. doi: 10.1016/j.antiviral.2019.104702
24. Ahn JI, Yoo J-Y, Kim TH, Kim YI, Ferguson SD, Fazleabas AT, et al. cAMP-response element-binding 3-like protein 1 (CREB3L1) is required for decidualization and its expression is decreased in women with endometriosis. Curr Mol Med (2016) 16:276–87. doi: 10.2174/1566524016666160225153659
25. Mikołajczyk M, Wirstlein P, Skrzypczak J. Aberrant claudin-4 transcript levels in eutopic endometrium of women with idiopathic infertility and minimal endometriosis. Ginekol Pol (2013) 84:90–4. doi: 10.17772/gp/1546
26. Bishop A, Cartwright JE, Whitley GS. Stanniocalcin-1 in the female reproductive system and pregnancy. Hum Reprod Update (2021) 27:1098–114. doi: 10.1093/humupd/dmab028
27. Allegra A, Marino A, Coffaro F, Lama A, Rizza G, Scaglione P, et al. Is there a uniform basal endometrial gene expression profile during the implantation window in women who became pregnant in a subsequent ICSI cycle? Hum Reprod (2009) 24:2549–57. doi: 10.1093/humrep/dep222
28. Aghajanova L, Altmäe S, Kasvandik S, Salumets A, Stavreus-Evers A, Giudice LC. Stanniocalcin-1 expression in normal human endometrium and dysregulation in endometriosis. Fertil Steril (2016) 106:681–91.e1. doi: 10.1016/j.fertnstert.2016.05.023
29. Ashley RL, Arreguin-Arevalo JA, Nett TM. Binding characteristics of the ovine membrane progesterone receptor alpha and expression of the receptor during the estrous cycle. Reprod Biol Endocrinol (2009) 7:42. doi: 10.1186/1477-7827-7-42
30. Weems YS, Nett TM, Rispoli LA, Davis TL, Johnson DL, Uchima T, et al. Effects of prostaglandin e and f receptor agonists in vivo on luteal function in ewes. Prostaglandins Other Lipid Mediat (2010) 92:67–72. doi: 10.1016/j.prostaglandins.2010.03.004
31. Hristovska A-M, Rasmussen LE, Hansen PBL, Nielsen SS, Nüsing RM, Narumiya S, et al. Prostaglandin E2 induces vascular relaxation by e-prostanoid 4 receptor-mediated activation of endothelial nitric oxide synthase. Hypertension (2007) 50:525–30. doi: 10.1161/HYPERTENSIONAHA.107.088948
32. Niringiyumukiza JD, Cai H, Xiang W. Prostaglandin E2 involvement in mammalian female fertility: ovulation, fertilization, embryo development and early implantation. Reprod Biol Endocrinol (2018) 16:43. doi: 10.1186/s12958-018-0359-5
33. Yang D, Jiang T, Liu J, Hong J, Lin P, Chen H, et al. Interferon-τ regulates prostaglandin release in goat endometrial stromal cells via JAB1 - unfolded protein response pathway. Theriogenology (2018) 113:237–46. doi: 10.1016/j.theriogenology.2018.03.007
34. Banerjee P, Sapru K, Strakova Z, Fazleabas AT. Chorionic gonadotropin regulates prostaglandin e synthase via a phosphatidylinositol 3-kinase-extracellular regulatory kinase pathway in a human endometrial epithelial cell line: implications for endometrial responses for embryo implantation. Endocrinology (2009) 150:4326–37. doi: 10.1210/en.2009-0394
35. Jabbour HN, Boddy SC. Prostaglandin E2 induces proliferation of glandular epithelial cells of the human endometrium via extracellular regulated kinase 1/2-mediated pathway. J Clin Endocrinol Metab (2003) 88:4481–7. doi: 10.1210/jc.2003-030297
36. Sales KJ, Katz AA, Millar RP, Jabbour HN. Seminal plasma activates cyclooxygenase-2 and prostaglandin E2 receptor expression and signalling in cervical adenocarcinoma cells. Mol Hum Reprod (2002) 8:1065–70. doi: 10.1093/molehr/8.12.1065
37. Sales KJ, Grant V, Catalano RD, Jabbour HN. Chorionic gonadotrophin regulates CXCR4 expression in human endometrium via e-series prostanoid receptor 2 signalling to PI3K-ERK1/2: implications for fetal-maternal crosstalk for embryo implantation. Mol Hum Reprod (2011) 17:22–32. doi: 10.1093/molehr/gaq069
38. Dominguez F, Galan A, Martin JJL, Remohi J, Pellicer A, Simón C. Hormonal and embryonic regulation of chemokine receptors CXCR1, CXCR4, CCR5 and CCR2B in the human endometrium and the human blastocyst. Mol Hum Reprod (2003) 9:189–98. doi: 10.1093/molehr/gag024
39. Nagasawa T, Hirota S, Tachibana K, Takakura N, Nishikawa S, Kitamura Y, et al. Defects of b-cell lymphopoiesis and bone-marrow myelopoiesis in mice lacking the CXC chemokine PBSF/SDF-1. Nature (1996) 382:635–8. doi: 10.1038/382635a0
40. Gentilini D, Busacca M, Di Francesco S, Vignali M, Viganò P, Di Blasio AM. PI3K/Akt and ERK1/2 signalling pathways are involved in endometrial cell migration induced by 17beta-estradiol and growth factors. Mol Hum Reprod (2007) 13:317–22. doi: 10.1093/molehr/gam001
41. Huang X, Liu H, Li R. Prostaglandin E2 promotes BeWo spheroids implantation in RL95-2 cell monolayers. Gynecol Endocrinol (2017) 33:548–52. doi: 10.1080/09513590.2017.1296125
42. Cilenti F, Barbiera G, Caronni N, Iodice D, Montaldo E, Barresi S, et al. A PGE2-MEF2A axis enables context-dependent control of inflammatory gene expression. Immunity (2021) 54:1665–82.e14. doi: 10.1016/j.immuni.2021.05.016
43. Glover LE, Crosby D, Thiruchelvam U, Harmon C, Chorcora CN, Wingfield MB, et al. Uterine natural killer cell progenitor populations predict successful implantation in women with endometriosis-associated infertility. Am J Reprod Immunol (2018) 79:e12817. doi: 10.1111/aji.12817
44. Fu B, Zhou Y, Ni X, Tong X, Xu X, Dong Z, et al. Natural killer cells promote fetal development through the secretion of growth-promoting factors. Immunity (2017) 47:1100–13.e6. doi: 10.1016/j.immuni.2017.11.018
45. Barber EM, Pollard JW. The uterine NK cell population requires IL-15 but these cells are not required for pregnancy nor the resolution of a listeria monocytogenes infection. J Immunol (2003) 171:37–46. doi: 10.4049/jimmunol.171.1.37
46. Wang W, Sung N, Gilman-Sachs A, Kwak-Kim J. T Helper (Th) cell profiles in pregnancy and recurrent pregnancy losses: Th1/Th2/Th9/Th17/Th22/Tfh cells. Front Immunol (2020) 11:2025. doi: 10.3389/fimmu.2020.02025
47. Ml S, Aj L, Aa A. Type I interferon regulation of natural killer cell function in primary and secondary infections. Expert Rev Vaccines (2013) 12:875-84. doi: 10.1586/14760584.2013.814871
48. Mattei F, Schiavoni G, Belardelli F, Tough DF. IL-15 is expressed by dendritic cells in response to type I IFN, double-stranded RNA, or lipopolysaccharide and promotes dendritic cell activation. J Immunol (2001) 167:1179–87. doi: 10.4049/jimmunol.167.3.1179
49. Matikainen S, Sareneva T, Ronni T, Lehtonen A, Koskinen PJ, Julkunen I. Interferon-alpha activates multiple STAT proteins and upregulates proliferation-associated IL-2Ralpha, c-myc, and pim-1 genes in human T cells. Blood (1999) 93:1980–91.
50. Huber JP, Farrar JD. Regulation of effector and memory T-cell functions by type I interferon. Immunology (2011) 132:466–74. doi: 10.1111/j.1365-2567.2011.03412.x
51. Saito S, Umekage H, Nishikawa K, Morii T, Narita N, Enomoto M, et al. Interleukin 4 (IL-4) blocks the IL-2-induced increased in natural killer activity and DNA synthesis of decidual CD16-CD56bright NK cells by inhibiting expression of the IL-2 receptor alpha, beta, and gamma. Cell Immunol (1996) 170:71–7. doi: 10.1006/cimm.1996.0135
52. Tezabwala BU, Johnson PM, Rees RC. Inhibition of pregnancy viability in mice following IL-2 administration. Immunology (1989) 67:115–9.
53. Monsanto SP, Edwards AK, Zhou J, Nagarkatti P, Nagarkatti M, Young SL, et al. Surgical removal of endometriotic lesions alters local and systemic proinflammatory cytokines in endometriosis patients. Fertil Steril (2016) 105:968–77.e5. doi: 10.1016/j.fertnstert.2015.11.047
54. Sathaliyawala T, Kubota M, Yudanin N, Turner D, Camp P, Thome JJC, et al. Distribution and compartmentalization of human circulating and tissue-resident memory T cell subsets. Immunity (2013) 38:187–97. doi: 10.1016/j.immuni.2012.09.020
55. Wessels JM, Domínguez MA, Leyland NA, Agarwal SK, Foster WG. Endometrial microbiota is more diverse in people with endometriosis than symptomatic controls. Sci Rep (2021) 11:18877. doi: 10.1038/s41598-021-98380-3
56. Ata B, Yildiz S, Turkgeldi E, Brocal VP, Dinleyici EC, Moya A, et al. The endobiota study: Comparison of vaginal, cervical and gut microbiota between women with stage 3/4 endometriosis and healthy controls. Sci Rep (2019) 9:2204. doi: 10.1038/s41598-019-39700-6
57. Moreno I, Codoñer FM, Vilella F, Valbuena D, Martinez-Blanch JF, Jimenez-Almazán J, et al. Evidence that the endometrial microbiota has an effect on implantation success or failure. Am J Obstet Gynecol (2016) 215:684–703. doi: 10.1016/j.ajog.2016.09.075
Keywords: endometriosis (EM), surgical treatment, endometrial receptivity, RNA sequencing (RNA-seq), secretory endometrium (SE), immune cell profile
Citation: Xiang R, Chen P, Zeng Z, Liu H, Zhou J, Zhou C, Peng J and Zeng H (2022) Transcriptomic analysis shows that surgical treatment is likely to influence the endometrial receptivity of patients with stage III/IV endometriosis. Front. Endocrinol. 13:932339. doi: 10.3389/fendo.2022.932339
Received: 29 April 2022; Accepted: 04 July 2022;
Published: 29 August 2022.
Edited by:
Sijia Lu, Yikon Genomics, ChinaReviewed by:
Kai-Fai Lee, The University of Hong Kong, Hong Kong SAR, ChinaMin Jin, Zhejiang University, China
Chi Chiu Wang, The Chinese University of Hong Kong, China
Copyright © 2022 Xiang, Chen, Zeng, Liu, Zhou, Zhou, Peng and Zeng. 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: Jintao Peng, pengjt3@mail.sysu.edu.cn; Haitao Zeng, zenghtao@mail.sysu.edu.cn
†These authors have contributed equally to this work and share first authorship