BioMed Research International
Volume 2020 (2020), Article ID 3962953, 7 pages
https://doi.org/10.1155/2020/3962953

Research Article

Integrated Bioinformatics Analysis Reveals Function and Regulatory Network of miR-200b-3p in Endometriosis

Wanxue Hu,¹ Qin Xie,¹ Yicong Xu,¹ Xin Tang,² and Hongbo Zhao1,3

¹Biomedical Engineering Research Center, Kunming Medical University, Kunming 650500, China
²School of Rehabilitation, Kunming Medical University, Kunming 650500, China
³Yunnan Key Laboratory of Stem Cell and Regenerative Medicine, Kunming 650500, China

Correspondence should be addressed to Xin Tang; tangxin8140730@163.com and Hongbo Zhao; zhao.hongbo@hotmail.com

Received 21 April 2020; Revised 1 July 2020; Accepted 16 July 2020; Published 30 July 2020

Guest Editor: Jiangning Song

Copyright © 2020 Wanxue Hu 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.

Abstract.

Objective. MicroRNAs play vital roles in the development of endometriosis. It is reported that miR-200b-3p is downregulated in endometriosis, although its mechanisms in this disease remain still unclear. Therefore, the purpose of this study was to explore the function and potential regulatory network of miR-200b-3p in endometriosis through database analysis. Methods. The endometriosis gene expression profiles were downloaded from the GEO database to screen differentially expressed genes (DEGs). The predicted and validated target genes of miR-200b-3p were obtained from miRWalk and miRTarBase database. Then, a comparison was performed between miR-200b-3p target genes and DEGs. GO enrichment and KEGG pathway analysis of the target genes was performed using clusterProfiler package. STRING was used to predict the protein-protein interaction among the proteins encoded by the target genes. Then, TransmiR, LncBase, StarBase, PROMO, and AnimalTFDB were employed to identify interactive transcription factors and lncRNAs of miR-200b-3p. Results. miR-200b-3p was associated with the transcription factors DNMT1, EZH2, HNF1B, JUN, MYB, ZEB1, and ZEB2 during the pathogenesis of endometriosis. The downstream 110 target genes were involved in the biological processes of positive regulation of MAPK cascade, muscle cell proliferation, organ growth, vasculogenesis, and axon development. KEGG analysis revealed that the main pathways related to miR-200b-3p were microRNAs in cancer, PI3K-Akt signaling pathway, colorectal cancer, and tight junction. In addition, four lncRNAs such as MALAT1, NEAT1, SNHG22, and XIST interacted with miR-200b-3p and were associated with transcription factors FOXP3 and YY1. Conclusion. The predicted target genes and molecular regulatory network of miR-200b-3p in endometriosis not only revealed its biological function but also provided a valuable guideline for further research.

1. Introduction

Endometriosis is a benign gynaecological disease characterized by endometrial glands and stroma located outside the uterine cavity. It affects 10–15% of women in their reproductive age and results in a significantly affected quality of life [1, 2]. Endometriosis symptoms include dyschezia, deep dyspareunia, dysuria, and dysmenorrhea, but the pathogenesis of endometriosis is still largely unknown [3]. Therefore, a better understanding of the molecular mechanism underlying endometriosis progression is of utmost importance to identify novel and more efficient diagnostic predictors for an accurate identification of the disease to use optimal therapeutic strategies.

MicroRNAs (miRNAs) are small, single-strand noncoding RNA molecules containing ~22 nucleotides, which regulate the expression of target genes and control a variety of cellular functions, including apoptosis, proliferation, migration, and invasion [4]. Furthermore, through a combinatorial action of miRNAs, transcription factors (TFs), and long noncoding RNAs (lncRNAs), more complex regulatory networks are often involved in various biological events [5]. Previous microarray studies demonstrated that some miRNAs are aberrantly expressed in endometriosis [6–8]. miR-200b-3p, one of the miRNAs downregulated in endometriosis [9, 10], belongs to the miR-200 family, and it is involved in the regulating of epithelial-to-mesenchymal transition (EMT) in endometriotic cells [11]. However, little is known about the potential function and regulatory network of miR-200b-3p in endometriotic lesions.

Therefore, in this study, an integrated analysis of gene expression profile datasets from GEO database was performed to screen differentially expressed genes (DEGs) in endometriosis. Then, the target genes of miR-200b-3p and the network it regulates were investigated using bioinformatics tools. Altogether, our results showed that miR-200b-3p might play critical roles in the pathology of endometriosis.

2. Materials and Methods

2.1. Expression Profile in Endometriosis

Two endometriosis microarray data sets were downloaded from the GEO website (http://www.ncbi.nlm.nih.gov/geo/). GSE7305 microarray data contained 10 endometriosis tissues and 10 normal endometrial tissues [12]. GSE7307 included 18 endometriosis samples and 23 healthy endometrial samples. In both datasets, the disease samples are ovarian endometriosis, which is the most frequently occurring type. The raw data were first corrected by RMA (Robust Multichip Averaging) algorithm and standardized together with log2 [13]. The DEGs were filtrated using the limma package [14]. The screening criteria were set at a false discovery rate-adjusted value (FDR) <0.05 and fold change >1.5. DEGs that appeared in both sets of data were selected for subsequent analysis.

2.2. Prediction of Target Genes

The miRWalk 2.0 tool (http://zmf.umm.uni-heidelberg.de/mirwalk2/) [15], which is linked to other 11 online prediction databases, was used to predict the putative target genes of miR-200b-3p. Only mRNA targets predicted by four miRNA prediction databases (miRWalk, miRanda, RNA22, and TargetScan), and in common, among them were considered for further analysis. Then, the validated mRNA targets of miR-200b-3p were downloaded from the miRTarBase database [16] and combined with predicted targets. Finally, the target genes of miR-200b-3p in endometriosis were obtained by comparing miR-200b-3p candidate target genes and DEGs from GEO.

2.3. Functional Enrichment Analysis of Target Genes

The clusterProfiler package [17] was used to evaluate GO functional (biological processes) and KEGG pathway analysis on the target genes. FDR <0.05 was used as the cutoff criterion.

2.4. Protein-Protein Interaction (PPI) Analysis of Target Genes

Target genes were submitted into STRING (https://string-db.org/) [18] to analyse the interaction between the encoded proteins. The pairs with combined scores >0.4 were used to construct the PPI network; then, the Cytoscape software was used for visualizing the network [19].

2.5. Prediction of miR-200b-3p TFs

TransmiR v2.0 is a database for TF-miRNA regulations, which covers 623 TFs and 785 miRNAs and contains 3730 literature-curated TF-miRNA regulations [20]. Therefore, the TFs linked with miR-200b-3p were predicted using TransmiR v2.0, and then, the ones included into the DEGs of endometriosis were screened.

2.6. Prediction of the Candidate lncRNAs for miR-200b-3p

LncBase Experimental v2 is an online software for experimental validation and computational derivation of the relationship between miRNA and lncRNA, which was used to predict lncRNAs associated with miR-200b-3p [21]. The StarBase database is based on CLIP-seq data and tumour samples (14 cancer types, >6000 samples), and analyses the interaction between lncRNA, miRNA, ceRNA, RBPs, and mRNA [22]. miRNA-lncRNA module in StarBase was used to screen lncRNAs binding to miR-200b-3p, and the results were compared to the ones obtained with LncBase Experimental v2.

2.7. Prediction of lncRNA Associated TFs

PROMO is a program that predicts TF binding sites in related sequences from a species through the TRANSFAC database [23]. AnimalTFDB is a source providing annotation and prediction of animal TFs from 97 animal genomes and contains 125,135 TF genes and 80,060 transcription cofactor genes [24]. The two databases were used to predict the TFs of lncRNAs and then to take the intersection of the results.

2.8. Construction of the miR-200b-3p Network Regulating Endometriosis

On the basis of the genomic information that was collected, the potential regulatory network of miR-200b-3p in endometriosis was constructed.

3. Results

3.1. Identification of DEGs

As regards DEG analysis, GSE7305 contained 4223 upregulated probes and 3830 downregulated probes in endometriosis, while GSE7307 contained 2360 upregulated probes and 2200 downregulated probes. After comparison, 2105 common upregulated probes corresponding to 1264 genes and 1867 common downregulated probes corresponding to 1269 genes were obtained.

3.2. Identification of miR-200b-3p Target Genes in Endometriosis

A total number of 482 predicted miR-200b-3p target genes that appeared in the four prediction databases were obtained. A total of 186 validated target genes associated with miR-200b-3p were obtained from the mirTarBase. The two sets were taken together, and 634 candidate target genes of miR-200b-3p were obtained and compared with the DEGs obtained from the microarray data analysis. Finally, 110 target genes of miR-200b-3p in endometriosis, including 56 upregulated genes and 54 downregulated genes, were obtained (Figure 1).

Figure 1: Venn diagram of the intersection analysis between miR-200b-3p target genes and DEGs in endometriosis.

3.3. GO and KEGG Pathway Enrichment Analysis of the Target Genes

GO analysis showed that the target genes were mainly enriched in biological processes such as positive regulation of MAPK cascade, muscle cell proliferation, organ growth, vasculogenesis, and axon development (Figure 2). In addition, KEGG pathway analysis indicated that these target genes were mainly involved in pathways of microRNAs in cancer, PI3K-Akt signaling pathway, colorectal cancer, and tight junction (Table 1).

Figure 2: Top 10 significantly enriched biological processes identified by GO analysis of miR-200b-3p target genes.

Table 1: KEGG pathway enrichment analysis on 110 target genes of miR-200b-3p.

Pathway ID
Pathway name
FDR
Included genes

hsa05206
MicroRNAs in cancer
0.0073
BCL2L11, FOXP1, ZEB2, ZFPM2, CCNE2, DNMT1, EGFR, EZH2, VEGFA, ZEB1

hsa04151
PI3K-Akt signaling pathway
0.0111
BCL2L11, FN1, NTRK2, PPP2R1B, CCNE2, EGFR, FGF9, IGF2, MYB, VEGFA

hsa05210
Colorectal cancer
0.0162
BCL2L11, JUN, TCF7L2, EGFR, PMAIP1

hsa05215
Prostate cancer
0.0212
ERG, TCF7L2, CCNE2, EGFR, ZEB1

hsa04530
Tight junction
0.0337
GATA4, JUN, PPP2R1B, MAP3K1, PARD6B, PRKAB1

hsa05165
Human papillomavirus infection
0.0493
FN1, FZD4, PPP2R1B, TCF7L2, CCNE2, EGFR, PARD6B, VEGFA

3.4. PPI Network Construction

STRING was used to extract the proteins expressed by the target genes that can interact with others. A total of 63 target genes were filtered into the PPI network, containing 63 nodes and 155 edges (Figure 3). The NetworkAnalyzer app in Cytoscape was used to calculate the node degree. The top 5 genes were the following: EGFR, EZH2, VEGFA, JUN, and FN1 with a degree >15; thus, they might be regulated by miR-200b-3p and play important roles in the progression of endometriosis.

Figure 3: Protein-protein interaction network analysis of the target genes. The PPI network consisting of 63 nodes and 155 edges was constructed using the target genes. The size of the circle represents the degree of a node.

3.5. TFs Associated to miR-200b-3p

A total of 75 TFs associated with miR-200b-3p were obtained from the TransmiR v2.0 database. Then, DNMT1, EZH2, HNF1B, JUN, MYB, ZEB1, and ZEB2 were found as TFs related to endometriosis by interacting with target genes of miR-200b-3p.

3.6. miR-200b-3p-lncRNAs Association

The lncRNAs associated with miR-200b-3p were predicted by LncBase Experimental v2 and StarBase databases, and the lncRNAs overlapping in both the databases were selected, such as MALAT1, NEAT1, SNHG22, and XIST. Thus, the mutual regulation of these four lncRNAs and miR-200b-3p might play an important role in the development of endometriosis.

3.7. Identification of lncRNA Associated TFs

PROMO and AnimalTFDB were used to predict the TFs that could be associated with four lncRNAs (MALAT1, NEAT1, SNHG22, and XIST) in the upstream 2 kb promoter sequence. The predicted intersections resulted in two TFs, such as FOXP3 and YY1.

3.8. Functional Regulatory Network of miR-200b-3p in Endometriosis

From the above prediction results, a miR-200b-3p network regulating endometriosis was constructed (Figure 4).

Figure 4: Predicted functional regulatory network of miR-200b-3p in endometriosis.

4. Discussion

Although endometriosis is a benign inflammatory disease, some of its characteristics are comparable to the ones of malignant tumours, such as abnormal cell migration, implantation, invasion, and angiogenesis. Several studies indicated that miR-200b-3p is decreased in endometriosis, but its specific pathogenesis in this disease is not yet clear. In this study, the miR-200b-3p network that regulates endometriosis was constructed by integrating microarray data sets and the target pair information from several interactome databases. The four relationships considered in this network were the following: (1) the regulatory pair of miR-200b-3p and target genes, (2) the regulatory pair of miR-200b-3p and TF genes, (3) the regulatory pair of miR-200b-3p and lncRNA genes, and (4) the regulatory pair of lncRNA genes and TF genes. A total of 110 target genes, 7 TF genes, and 4 lncRNAs associated with miR-200b-3p were identified in endometriosis. The functions of the target genes were analysed, and the enriched GO terms were positive regulation of MAPK cascade, muscle cell proliferation, organ growth, vasculogenesis, and axon development. These results suggested that the identified miR-200b-3p target genes in endometriosis and the associated regulators (TFs and lncRNAs) might play an important role in the induction of endometriosis.

By integrating the information from gene expression profiles, the target genes of miRNA can be found more accurately; thus, the genes regulated by miR-200b-3p in endometriosis can be also discovered. A total of 110 target genes of miR-200b-3p were screened, and many of them play important roles in the development of endometriosis. They included FOXP1, GATA6, JUN, and EGFR involved in muscle cell proliferation [25–28]; NTRK2, ZFPM, and VEGFA involved in angiogenesis [29, 30]; and DNMT1, EZH2, FOXP1, VEGFA, ZEB1, and ZEB2 involved in EMT [31, 32]. Recent studies found that EMT plays an important role in the occurrence and development of endometriosis [33]. The miRNA-200 family is a recognized epithelial cell marker and is a major regulator of EMT [34]. Other studies showed that miR-200b-3p downregulates ZEB1 and ZEB2 in endometriosis patients and participates in the feedback regulation of EMT [35], which more easily promotes cell invasion and metastasis. Our results showed that EGFR, EZH2, VEGFA, JUN, and FN1 are important nodes in the target genes PPI network. These results indicated the high accuracy in predicting target genes, and they also proved that miR-200b-3p had important regulatory effects on endometriosis. In addition, some predicted genes related to endometriosis were not discovered during the research process and reported. Although these genes have not been fully studied and verified, they might become important targets in the treatment of endometriosis.

miRNAs not only regulate the development of several diseases by interacting with protein-coding genes but also promote or inhibit cell proliferation, apoptosis, and metastasis by interacting with lncRNAs. The lncRNAs associated with miR-200b-3p were predicted by LncBase and StarBase databases and identified 4 overlapping lncRNAs (MALAT1, NEAT1, SNHG22, and XIST). Therefore, these four lncRNAs might have a mutual regulation with miR-200b-3p. Liang et al. found that the expression of MALAT1 is significantly upregulated in ectopic endometrial tissue and is negatively correlated with the expression of miR-200c [36], suggesting that MALAT1 is a target of miR-200c in primary endometrial stromal cells. Since miR-200c and miR-200b are both belonging to the miRNA-200 family, this could mean that MALATI is a potential target of miR-200b-3p regulating endometrial stromal cells. However, the relationship between these four lncRNAs and miR-200b needs further experimental verification.

The common predicted TFs of the four lncRNAs were FOXP3 and YY1. FOXP3 is essential in the development and function of regulatory T cells. Ectopic Foxp3 expression can phenotypically convert effector T cells to regulatory T cells [37]. A recent research found that endometriosis is related to changes in Treg and Th17 cells in the pelvis, leading to the survival and implantation of ectopic endometrial lesions [38]. The percentage of CD25⁺FOXP3⁺Treg cells is significantly increased in the peritoneal fluid of women with endometriosis [39]. Furthermore, the combination of FCRL3 and FOXP3 genotypes increases the risk of endometriosis, as observed in the Brazilian population samples [40]. YY1 is a TF with diverse and complex biological functions and has a role in controlling X chromosome inactivation (XCI) [41]. YY1 binds directly the Xist 5 region to trigger the activity of the Xist promoter [42], and loss of YY1 prevents Xist upregulation during the initiation and maintenance of X-inactivation [43]. Furthermore, YY1 deletion or disruption in activated female T and B cells disrupts Xist RNA localization to the inactive X with minimal impact on Xist expression, suggesting that YY1 may regulate XCI maintenance through the mediation of Xist RNA localization in lymphocytes [44]. In addition, endometriosis is an estrogen-dependent disease characterized by chronic inflammation and clinically by chronic pelvic pain. Current evidence suggests that mast cells and microglia are considered coordinators of peripheral and central inflammatory processes, and may trigger and sustain the persistent inflammatory process present in the pathogenesis of endometriosis [45]. Han et al. confirmed that YY1 upregulates NEAT1 in the oxygen-glucose deprivation/reoxygenation injury of microglial cells [46].

In conclusion, based on the results of various bioinformatics methods, our hypothesis is that miR-200b-3p is associated with TFs such as DNMT1, EZH2, HNF1B, JUN, MYB, ZEB1, and ZEB2 in the occurrence and development of endometriosis, regulating downstream target genes, and it is also regulated by four lncRNAs. Furthermore, these four lncRNAs interact with the TFs FOXP3 and YY1. Although the prediction of the molecular regulatory network of miR-200b-3p by bioinformatics needs further experiments to be confirmed, these predictions may represent a valuable guideline in the mechanism triggering endometriosis.

Data Availability

The microarray raw data used in this study are freely available from GEO datasets GSE7305 and GSE7307. In addition, the datasets used for network construction were obtained from the tool websites.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

XT and HBZ conceived and designed the study. WXH conducted the bioinformatics analysis and drafted the manuscript. QX and YCX performed the data collection. XT and HBZ revised the important intellectual content of the manuscript. All authors read and approved the final manuscript.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (grant no. 81360336) and the Joint Special Funds for the Department of Science and Technology of Yunnan Province-Kunming Medical University (grant no. 2018FE001(-197)).

References

S. E. Bulun, “Endometriosis,” New England Journal of Medicine, vol. 360, no. 3, pp. 268–279, 2009.

M. Hickey, K. Ballard, and C. Farquhar, “Endometriosis,” BMJ, vol. 348, no. mar19 5, article g1752, 2014.

P. Vercellini, P. Vigano, E. Somigliana, and L. Fedele, “Endometriosis: pathogenesis and treatment,” Nature Reviews: Endocrinology, vol. 10, no. 5, pp. 261–275, 2014.

T. A. Farazi, J. I. Hoell, P. Morozov, and T. Tuschl, “MicroRNAs in human cancer,” Advances in Experimental Medicine and Biology, vol. 774, pp. 1–20, 2013.

S. Arora, R. Rana, A. Chhabra, A. Jaiswal, and V. Rani, “miRNA-transcription factor interactions: a combinatorial regulation of gene expression,” Molecular Genetics and Genomics, vol. 288, no. 3-4, pp. 77–87, 2013.

E. M. C. Ohlsson Teague, K. H. Van der Hoek, M. B. Van der Hoek et al., “MicroRNA-regulated pathways associated with endometriosis,” Molecular Endocrinology, vol. 23, no. 2, pp. 265–275, 2009.

S. M. Hawkins, C. J. Creighton, D. Y. Han et al., “Functional microRNA involved in endometriosis,” Molecular Endocrinology, vol. 25, no. 5, pp. 821–832, 2011.

A. Braza-Boïls, J. Marí-Alexandre, J. Gilabert et al., “MicroRNA expression profile in endometriosis: its relation to angiogenesis and fibrinolytic factors,” Human Reproduction, vol. 29, no. 5, pp. 978–988, 2014.

M. Saare, K. Rekker, T. Laisk-Podar et al., “High-throughput sequencing approach uncovers the miRNome of peritoneal endometriotic lesions and adjacent healthy tissues,” PLoS One, vol. 9, no. 11, article e112630, 2014.

K. Rekker, M. Saare, A. M. Roost et al., “Circulating miR-200-family micro-RNAs have altered plasma levels in patients with endometriosis and vary with blood collection time,” Fertility and Sterility, vol. 104, no. 4, pp. 938–946.e2, 2015.

J. C. Eggers, V. Martino, R. Reinbold et al., “microRNA miR-200b affects proliferation, invasiveness and stemness of endometriotic cells by targeting ZEB1, ZEB2 and KLF4,” Reproductive Biomedicine Online, vol. 32, no. 4, pp. 434–445, 2016.

A. Hever, R. B. Roth, P. Hevezi et al., “Human endometriosis is associated with plasma cells and overexpression of B lymphocyte stimulator,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 30, pp. 12451–12456, 2007.

L. Gautier, L. Cope, B. M. Bolstad, and R. A. Irizarry, “affy--analysis of Affymetrix GeneChip data at the probe level,” Bioinformatics, vol. 20, no. 3, pp. 307–315, 2004.

M. E. Ritchie, B. Phipson, Y. H. Di Wu, C. W. Law, W. Shi, and G. K. Smyth, “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, pp. e47–e47, 2015.

H. Dweep and N. Gretz, “miRWalk2.0: a comprehensive atlas of microRNA-target interactions,” Nature Methods, vol. 12, no. 8, p. 697, 2015.

C. H. Chou, S. Shrestha, C. D. Yang et al., “miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions,” Nucleic Acids Research, vol. 46, no. D1, pp. D296–D302, 2018.

G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.

D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–D613, 2019.

P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.

Z. Tong, Q. Cui, J. Wang, and Y. Zhou, “TransmiR v2.0: an updated transcription factor-microRNA regulation database,” Nucleic Acids Research, vol. 47, no. D1, pp. D253–D258, 2019.

M. D. Paraskevopoulou, I. S. Vlachos, D. Karagkouni et al., “DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts,” Nucleic Acids Research, vol. 44, no. D1, pp. D231–D238, 2016.

J.-H. Li, S. Liu, H. Zhou, L.-H. Qu, and J.-H. Yang, “starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data,” Nucleic Acids Research, vol. 42, no. D1, pp. D92–D97, 2014.

D. Farré, R. Roset, M. Huerta et al., “Identification of patterns in biological sequences at the ALGGEN server: PROMO and MALGEN,” Nucleic Acids Research, vol. 31, no. 13, pp. 3651–3653, 2003.

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.

L. Fu, J. E. Girling, and P. A. Rogers, “Expression of Fox head protein 1 in human eutopic endometrium and endometriosis,” Reproductive Sciences, vol. 15, no. 3, pp. 243–252, 2008.

L. A. Bernardi, M. T. Dyson, H. Tokunaga et al., “The essential role of GATA6 in the activation of estrogen synthesis in endometriosis,” Reproductive Sciences, vol. 26, no. 1, pp. 60–69, 2017.

M. T. Beste, N. Pfaffle-Doyle, E. A. Prentice et al., “Molecular network analysis of endometriosis reveals a role for c-Jun-regulated macrophage activation,” Science Translational Medicine, vol. 6, no. 222, article 222ra16, 2014.

K. Chatterjee, S. Jana, P. DasMahapatra, and S. Swarnakar, “EGFR-mediated matrix metalloproteinase-7 up-regulation promotes epithelial-mesenchymal transition via ERK1-AP1 axis during ovarian endometriosis progression,” The FASEB Journal, vol. 32, no. 8, pp. 4560–4572, 2018.

X. Yu, H. Ren, T. Liu, M. Yong, and H. Zhong, “Expression and significance of ERβ and TrkB in endometriosis,” Clinical and Experimental Obstetrics and Gynecology, vol. 43, no. 1, pp. 75–81, 2016.

A. Vodolazkaia, B. T. Yesilyurt, C. M. Kyama et al., “Vascular endothelial growth factor pathway in endometriosis: genetic variants and plasma biomarkers,” Fertility and Sterility, vol. 105, no. 4, pp. 988–996, 2016.

K. Y. Hsiao, M. H. Wu, N. Chang et al., “Coordination of AUF1 and miR-148a destabilizes DNA methyltransferase 1 mRNA under hypoxia in endometriosis,” Molecular Human Reproduction, vol. 21, no. 12, pp. 894–904, 2015.

Q. Zhang, P. Dong, X. Liu, N. Sakuragi, and S.-W. Guo, “Enhancer of Zeste homolog 2 (EZH2) induces epithelial-mesenchymal transition in endometriosis,” Scientific Reports, vol. 7, no. 1, p. 6804, 2017.

K. Proestling, P. Birner, S. Gamperl et al., “Enhanced epithelial to mesenchymal transition (EMT) and upregulated MYC in ectopic lesions contribute independently to endometriosis,” Reproductive Biology and Endocrinology, vol. 13, no. 1, 2015.

P. S. Mongroo and A. K. Rustgi, “The role of the miR-200 family in epithelial-mesenchymal transition,” Cancer Biology & Therapy, vol. 10, no. 3, pp. 219–222, 2014.

S. M. Park, A. B. Gaur, E. Lengyel, and M. E. Peter, “The miR-200 family determines the epithelial phenotype of cancer cells by targeting the E-cadherin repressors ZEB1 and ZEB2,” Genes & Development, vol. 22, no. 7, pp. 894–907, 2008.

Z. Liang, Y. Chen, Y. Zhao et al., “miR-200c suppresses endometriosis by targeting MALAT1 in vitro and in vivo,” Stem Cell Research & Therapy, vol. 8, no. 1, p. 251, 2017.

S. F. Ziegler, “FOXP3: of mice and men,” Annual Review of Immunology, vol. 24, no. 1, pp. 209–226, 2006.

K. N. Khan, K. Yamamoto, A. Fujishita et al., “Differential levels of regulatory T cells and T-helper-17 cells in women with early and advanced endometriosis,” The Journal of Clinical Endocrinology and Metabolism, vol. 104, no. 10, pp. 4715–4729, 2019.

J. Olkowska-Truchanowicz, K. Bocian, R. B. Maksym et al., “CD4⁺ CD25⁺ FOXP3⁺ regulatory T cells in peripheral blood and peritoneal fluid of patients with endometriosis,” Human Reproduction, vol. 28, no. 1, pp. 119–124, 2013.

C. P. Barbosa, J. S. Teles, T. G. Lerner et al., “Genetic association study of polymorphisms FOXP3 and FCRL3 in women with endometriosis,” Fertility and Sterility, vol. 97, no. 5, pp. 1124–1128, 2012.

N. Xu, M. E. Donohoe, S. S. Silva, and J. T. Lee, “Evidence that homologous X-chromosome pairing requires transcription and Ctcf protein,” Nature Genetics, vol. 39, no. 11, pp. 1390–1396, 2007.

Y. Jeon and J. T. Lee, “YY1 tethers Xist RNA to the inactive X nucleation center,” Cell, vol. 146, no. 1, pp. 119–133, 2011.

M. Makhlouf, J.-F. Ouimette, A. Oldfield, P. Navarro, D. Neuillet, and C. Rougeulle, “A prominent and conserved role for YY1 in Xist transcriptional activation,” Nature communications, vol. 5, no. 1, 2014.

J. Wang, C. M. Syrett, M. C. Kramer, A. Basu, M. L. Atchison, and M. C. Anguera, “Unusual maintenance of X chromosome inactivation predisposes female lymphocytes for increased expression from the inactive X,” Proceedings of the National Academy of Sciences of the United States of America, vol. 113, no. 14, pp. E2029–E2038, 2016.

A. Graziottin, “Inflammation and chronic pelvic pain: a biological trigger for depression in women?” Journal of Depression and Anxiety, vol. 3, no. 1, 2014.

D. Han and Y. Zhou, “YY1-induced upregulation of lncRNA NEAT1 contributes to OGD/R injury-induced inflammatory response in cerebral microglial cells via Wnt/β-catenin signaling pathway,” In vitro cellular & developmental biology - Animal, vol. 55, no. 7, pp. 501–511, 2019.