Skip to main content
Download PDF

Integration analysis of microRNA and mRNA paired expression profiling identifies deregulated microRNA-transcription factor-gene regulatory networks in ovarian endometriosis

Reproductive Biology and Endocrinology volume 16, Article number: 4 (2018) Cite this article

Abstract

Background

The etiology and pathophysiology of endometriosis remain unclear. Accumulating evidence suggests that aberrant microRNA (miRNA) and transcription factor (TF) expression may be involved in the pathogenesis and development of endometriosis. This study therefore aims to survey the key miRNAs, TFs and genes and further understand the mechanism of endometriosis.

Methods

Paired expression profiling of miRNA and mRNA in ectopic endometria compared with eutopic endometria were determined by high-throughput sequencing techniques in eight patients with ovarian endometriosis. Binary interactions and circuits among the miRNAs, TFs, and corresponding genes were identified by the Pearson correlation coefficients. miRNA-TF-gene regulatory networks were constructed using bioinformatic methods. Eleven selected miRNAs and TFs were validated by quantitative reverse transcription-polymerase chain reaction in 22 patients.

Results

Overall, 107 differentially expressed miRNAs and 6112 differentially expressed mRNAs were identified by comparing the sequencing of the ectopic endometrium group and the eutopic endometrium group. The miRNA-TF-gene regulatory network consists of 22 miRNAs, 12 TFs and 430 corresponding genes. Specifically, some key regulators from the miR-449 and miR-34b/c cluster, miR-200 family, miR-106a-363 cluster, miR-182/183, FOX family, GATA family, and E2F family as well as CEBPA, SOX9 and HNF4A were suggested to play vital regulatory roles in the pathogenesis of endometriosis.

Conclusion

Integration analysis of the miRNA and mRNA expression profiles presents a unique insight into the regulatory network of this enigmatic disorder and possibly provides clues regarding replacement therapy for endometriosis.

Background

Endometriosis is characterized by the presence of endometrial tissue (endometrial glandular and stromal) abnormally outside the uterine cavity [1]. In women of reproductive age, there is a 6–10% incidence rate, which may even reach 35–50% in women who are enduring pelvic pain and/or infertility [2, 3]. Despite its high prevalence, the etiology and pathophysiology of endometriosis are still not fully understood. It is widely accepted that endometriosis is a multifactorial, polygenic disorder [4]. A large amount of molecular aberrations exist between the ectopic endometrium (EC) (endometriotic lesions) and the eutopic endometrium (EU), which could potentially explain the mechanism of the abnormal growth of EU outside the uterus.

microRNAs (miRNAs) are endogenous, single-stranded, 18–22 nucleotide RNAs that mainly inhibit gene expression at the post-transcriptional level [5]. They negatively regulate the translation of target messenger RNAs (mRNAs) via the sequence specific recognition of the “seed sequence” and repress translation and/or degrade the target mRNAs according to the degree of complementary nucleotides [6]. Transcription factors (TFs) are a group of proteins that regulate gene expression by binding to specific DNA sequences at the transcriptional level [7]. Depending on its function, the transcription of an adjacent gene is either activated or repressed [8]. As two vital gene regulatory molecules, miRNAs and TFs can share common target genes and exert effects on each other [9]. As such, the miRNA expression levels can be regulated by TFs, and the mRNAs encoding TFs can be inhibited by miRNAs [10]. In addition, miRNAs and TFs can form a feed-back loop (FBL) or feed-forward loop (FFL) in cooperation to fine tune gene expression [11]. The combinatorial regulation of TFs and miRNAs is suggested to play vital roles in various biological processes and disease pathogenesis [12,13,14,15].

Accumulating evidence suggests that aberrant miRNA and TF expression may be involved in the pathogenesis and development of endometriosis [16,17,18,19,20,21,22,23,24,25]. Although differences in miRNA expression profiling between EC and EU from women with endometriosis have previously been reported, almost no previous studies have focused on the integrative analysis of miRNA-mRNA interactions and solved the complex problem of the identification of miRNA-TF-gene regulatory networks. With the aim of further understanding the abnormal molecule changes in the occurrence and development of endometriosis, we used high-throughput sequencing techniques to detect miRNA and mRNA expression profiling in paired EC and EU. Joint analysis of the miRNA-TF-gene regulatory network can help us survey the key genomic factors and understand the mechanism of endometriosis at the molecular level.

Methods

Patient samples

The study protocol was approved by the Local Ethical Committee of the Chinese PLA General Hospital, and each patient was required to provide written consent.

Sixty tissue samples (30 paired EC and EU) from 30 patients with ovarian endometriosis were obtained from the Chinese PLA general hospital. Among them, eight paired EC and EU samples were selected and prepared for small RNA sequencing and mRNA sequencing. The remaining samples were used for validation. All of the patients were confirmed to have endometriosis by histological examination and diagnosed as being at a moderate to severe (III-IV) stage by the revised American Fertility Society (r-AFS) classification during laparoscopic surgery [26]. EC samples were obtained during laparoscopy, and EU samples were obtained through curettage before the laparoscopic procedure. Only patients in the secretory phase of the menstrual cycle, which was confirmed by the method of Noyes et al. [27], and without any hormonal treatment history, were included in the study. Clinical characteristics of the patients are listed in Additional file 1.

Tissue processing, RNA extraction and quality control

All tissue samples were divided into two parts: one half was fixed and prepared for pathological examination and the other half was placed in RNA-later solution (Sigma Aldrich, Poland) at 4 °C for 24 h and subsequently transferred to −80 °C until further use. Total RNA was extracted using a single-step acid guanidinium thiocyanate-phenol-chloroform method [28]. The quality and purity of RNA were examined by a Nanodrop 8000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). RNA integrity was analyzed using the RNA 6000 Nano Kit and Small RNA Kit with the Bioanalyzer 2100 (Agilent, Santa Clara, CA). Samples with an absorbance wavelength ratio (A260/A280) ≥ 1.9 and an RNA integrity number ≥ 8 were included.

Small RNA sequencing (small RNA-seq) and data analysis

For small RNA sequencing, 3 μg of total RNA was used for library preparation as previously reported with little modification [29]. The purified libraries were sequenced on an Illumina HiSeq 2000 platform by Annoroad Genome (Beijing, China). The detailed protocol is described in Additional file 2. In brief, the raw data were processed with Python scripts to ensure its quality and were then filtered by Q30 statistics. The clean data were mapped to the Ensemble database (GRCh37) (http://grch37.ensembl.org/index.html) [30] using Bowtie (v1.01) [31] and compared with miRBase (Release 21) [32] to identify mature miRNAs. For each sample, the count and RPM (reads per million total reads) value of the miRNAs were collected. DEGseq (v1.18.0) [33] was used for differential expression analysis with the parameters of differentially expressed miRNAs (DEM) set with a false discovery rate (FDR) < 0.05 and |log2 fold change (FC)| ≥ 1.

mRNA sequencing (mRNA-seq) and data analysis

For mRNA sequencing, 3 μg of RNA per sample was prepared for library preparation. The Ribo-Zero Gold Kit (Epicentre, USA) and NEB Next Ultra RNA Library Prep Kit (NEB, Ipswich, USA) were used for rRNA removal and library construction following the manufacturer’s protocols. For high-throughput sequencing, paired-end 150-bp sequencing of the cDNAs was performed on an Illumina HiSeq4000 system (Illumina, USA) conducted by Annoroad Genome (Beijing, China). The raw data were processed with Perl scripts to ensure the quality of the data used in subsequent analyses (details in Additional file 2). Bowtie2 (v2.2.3) [34] was used to build the genome index, and clean data were mapped to the human genome (GRCh37) (http://grch37.ensembl.org/index.html) using TopHat (v2.0.12) [35]. The read counts of each gene were counted by HTSeq (v0.6.0) [36], and the reads per kilobase of one gene per million reads (RPKM) were calculated to estimate the expression level of genes in each sample. DEGseq (v1.18.0) [33] was used to analyze differentially expressed genes (DEGs) with parameters of FDR < 0.05 and |log2FC| ≥ 1.

Identification of miRNA/TF/gene interactions and a regulatory network

The miRNA and mRNA expression profiles were analyzed by MAGIA2 (http://gencomp.bio.unipd.it/magia2/) to identify miRNA-TF-gene binary interactions and circuits [37]. The miRNA-target interactions were predicted by the TargetScan database [38] with a z-score ≥ 0.7. TF-miRNA interactions were identified from the mirGen2.0 [39] and TransmiR [40] databases, and TF-gene interactions were obtained from the ‘TFBS conserved’ track of the UCSC genome annotation for human (version hg19) (http://genome.ucsc.edu/) database, restricting the z-score to ≥ 3. Two types of mixed regulatory circuits were also identified: (i) FFL, which is a TF that regulates both a given miRNA and their common target gene and (ii) FBL, which is a miRNA that regulates both a given TF and their common regulated gene. The Pearson correlation coefficient (PCC) was used to measure the relationships with a PCC > 0.6 and a P-value < 0.05. The integrative miRNA-TF-gene regulatory network was visualized by Cytoscape (v3.2) [41].

Functional analysis

To explore the functional roles of DEGs in the miRNA-TF-gene networks, we used DAVID, which integrates the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases to analyze biological functions [42]. The enrichment values of GO terms and KEGG pathways were implemented by the hypergeometric test, and q (adjusted as p-value) < 0.05 was considered to be significantly enriched.

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) validation

Six miRNAs (miR-34c-5p, miR-106a-5p, miR-182-5p, miR-200a-3p, miR-449b-5p, and miR-615-3p) and five TFs (CEBPA, FOXC1, E2F1, GATA1, and HNF4A) were selected for validation analysis. U6 was selected as the miRNA endogenous control, and GAPDH was used as the mRNA endogenous control. Primer sequences are listed in Additional file 3. miRNA cDNA synthesis was conducted with the miRNA cDNA Kit (ComWin Biotech, Beijing), and total cDNA synthesis was conducted by a RevertAid™ First Strand cDNA Synthesis Kit (Thermo, USA). Relative miRNA expression was determined according to the miRNA Real-Time PCR Assay Kit (ComWin Biotech, Beijing), and relative mRNA expression was determined according to the THUNDERBIRD™ SYBR qPCR Mix (TOYOBO, Japan). qRT-PCR was performed on an ABI PRISM 7500 (Applied Biosystems, CA, USA). Relative expression was calculated using the ABI PRISM 7500 version 2.0.6 software (Applied Biosystems, CA, USA) with the 2−ΔΔCt method [43].

Results

Small RNA-seq analysis and DEM profiling

The small RNA-seq generated 216,333,583 reads with an average of 13,520,849 reads per sample. The FastQC quality test showed that 182,420,380 of the reads had a Q-score ≥ 30; these reads were considered in further analyses. A proportion of clean reads was mapped to the miRBase database (v21.0) for each sample and is listed in Additional file 4. Of the 1574 miRNAs analyzed by small RNA sequencing, 1163 with a normalized count value greater than 1 in at least one sample in both the EC and EU groups were further analyzed. Unsupervised hierarchical clustering showed that the overall miRNA profile of EC was wholly different from that of EU (Fig.1a). Using the selection criteria mentioned above, 107 DEMs were found, including 41 up-regulated miRNAs and 66 down-regulated miRNAs, between EC and paired EU. Among them, miR-514b-3p was the most up-regulated, with an FC of 330.31, whereas miR-375 was the most down-regulated, with an FC of 154.52 (Additional file 5).

Fig. 1
figure 1

Heat map showing differentially expressed miRNAs and mRNAs from ectopic endometria compared with paired eutopic endometria. Rows represent miRNAs (a) or mRNAs (b); columns represent profiled samples. The relative expression is depicted according to the color scale. Red indicates up-regulation; blue indicates down-regulation. EU means eutopic endometrium and EC means ectopic endometrium

Full size image

mRNA-seq analysis and DEG profiling

mRNA-seq generated 1,748,534,484 reads with an average of 109,283,405 reads per sample. The FastQC quality test showed that 1,721,709,042 (98.46%) reads reached Q30; these reads were considered for further analyses. Of these reads, 93.17% were mapped to hg19, and 98.56% were uniquely aligned. The detailed filtering and mapping data are listed in Additional file 6. Unsupervised hierarchical clustering analysis showed that the transcriptome profiling of the EC group was significantly different from that of the EU group (Fig. 1b). Using the same criteria that were used for the miRNAs, a total of 6112 DEGs with 3228 up-regulated mRNA transcripts and 2884 down-regulated mRNA transcripts were identified. The most up- and down-regulated mRNAs were IGHV3–49 (ENSG00000211965) and FXYD4 (ENSG00000150201), with FCs of 3253.41 and 1024.16, respectively. Additional file 7 lists the top 25 up- and down-regulated mRNAs identified by our sequencing.

Identification of miRNA/TF/gene interactions and regulatory networks

One hundred seven DEMs and 6112 DEGs identified in the previous steps were considered for the miRNA/TF/gene binary interaction and regulatory network analysis.

According to the standard of absolute PCC > 0.6, correlation analysis revealed 379 paired miRNA-target regulatory interactions between 15 miRNAs (miR-20b-5p, miR-30d-5p, miR-34c-5p, miR-106a-5p, miR-133b, miR-141-3p, miR-182-5p, miR-183-5p, miR-196b-5p, miR-200a-3p, miR-200b-3p, miR-200c-3p, miR-363-3p, miR-449a, and miR-449b-5p), five TFs (GATA6, GATA2, FOXO1, FOXC1, and CEBPA) and 223 genes. Although miRNAs usually suppress target mRNAs, positively paired miRNA-mRNA correlations are occasionally observed under certain conditions [44]. For this reason, we identified all of the differentially expressed miRNA-gene pairs, including 174 negative correlating pairs and 205 positive correlating pairs. In the case of TF-gene interactions, 310 regulatory pairs involving seven TFs (E2F1, CEBPA, FOXC1, FOXL1, FOXO1, GATA1, and SOX9) and 261 corresponding genes were identified. For TF-miRNA interactions, a total of 23 TF-miRNA pairs between seven TFs (E2F1, E2F2, FOXJ1, FOXL1, GATA1, HNF4A, and SOX9) and 19 miRNAs (miR-20b-5p, miR-34b-3p, miR-34b-5p, miR-34c-5p, miR-106a-5p, miR-133b, miR-182-5p, miR-183-5p, miR-196b-5p, miR-200a-3p, miR-200b-3p, miR-202-3p, miR-216a-5p, miR-224-5p, miR-363-3p, miR-449a, miR-449b-5p, miR-449c-5p, and miR-708-5p) were involved in the comparison. Moreover, we revealed five mixed circuits, including one FFL (FOXL1/miR-133b/PRDM1) and four FBLs (miR-182-5p/FOXO1/NUAK1, miR-182-5p/FOXO1/SCN9A, miR-133b/FOXC1/EYA4, and miR-133b/FOXC1/LHFP).

Based on the relationships of binary interactions and circuits, we constructed the miRNA-TF-gene regulatory network. The final integrated network contained 22 miRNAs, 12 TFs and 430 genes (Table 1, Fig. 2). To identify key regulators in the integrated networks, we also calculated the node degree of the involved miRNAs and TFs. The node degree is a centrality measure method that directly counts the links of one node. The few nodes that have a high degree in the network are called hubs, which reflect their importance in a large number of interactions [45]. As a result, the node degrees of the selected miRNAs and TFs ranged from 1 to 117 in our data. The top five key regulators in the miRNA-TF-gene network of endometriosis were FOXC1, FOXO1, miR-182-5p, miR-106a-5p, and CEBPA, which were directly connected to 117, 64, 57, 49, and 49 corresponding targets, respectively (Table 1).

Table 1 The key regulators and their targets identified in the integrated network in endometriosis
Full size table
Fig. 2
figure 2

The integrated endometriosis-related miRNA-transcription factor-gene regulation network. Red represents miRNA; yellow represents transcription factor; green represents target gene

Full size image

Functional analysis

Functional enrichment analysis of all of the target genes involved in the networks provided us with an overall clue regarding their functional roles in endometriosis. GO analysis revealed that the target genes in the GO cellular component (GO-CC) category were enriched in 16 terms, such as the plasma membrane, proteinaceous extracellular matrix and Z disc. The GO molecular function (GO-MF) category was enriched in six terms, including sequence-specific DNA binding, transcriptional activator activity and cell adhesion molecule binding. In the biological process category of GO terms (GO-BP), 20 terms included cell adhesion, whereas embryonic skeletal system morphogenesis, extracellular matrix organization and angiogenesis were significantly enriched (Fig. 3a). KEGG pathway analysis demonstrated that the DEG targets were associated with 27 pathways, such as the calcium-signaling pathway, cAMP signaling pathway, and transcriptional misregulation in cancer pathway (Fig. 3b).

Fig.3
figure 3

a GO terms enriched by differentially expressed genes that involved in regulatory network. BP represents biological processes; CC represents cellular components; MF represents molecular functions. b KEGG pathways enriched by differentially expressed genes that involved in regulatory network

Full size image

Validation by qRT-PCR

The qRT-PCR results were consistent with the sequencing results in that all six selected miRNAs (miR-34c-5p, miR-106a-5p, miR-182-5p, miR-200a-3p, miR-449b-5p, and miR-615-3p) and five TFs (CEBPA, FOXC1, E2F1, GATA1, and HNF4A) were differentially expressed with the same trend (up- or down-regulated) of sequencing and reached statistical significance (Fig. 4).

Fig. 4
figure 4

qRT-PCR analysis of the six miRNAs expressions (a) and five transcription factors expressions (b) in ectopic endometria (n = 22) and paired eutopic endometria (n = 22). Grey column represents ectopic endometrium; white column represents eutopic endometrium. Significant fold changes in ectopic vs. eutopic endometrium are marked by *, P ≤ 0.01; ** P ≤ 0.001. The error bar denotes SEM (standard error of the mean)

Full size image

Discussion

Several microarray studies were conducted to obtain DEMs between EC and EU from women with endometriosis [18, 46, 47], but overlapping data are limited, which suggests that the molecular changes caused by endometriosis might be more complex than we assumed. Additionally, miRNA target prediction methods based on sequence complementarity have been used to annotate the functions of DEMs in these studies, which may contain a high false positive rate and false negative rate [48]. Therefore, the functions of most of the identified DE miRNAs in endometriosis remain unclear due to a lack of definite target genes and confounding downstream validations. Network analysis allows structured grouping of miRNAs, TFs and target genes; thus, network construction is essential for identifying critical molecular changes and understanding their functional roles. To our knowledge, this is the first study to establish concomitant miRNA-TF-gene binary interactions and regulatory networks between the EC and the EU from patients with endometriosis. At first, our study revealed 107 DEMs and 6112 DEGs by sequencing. The large differential expression profiling was refined by PCC to 22 functionally DE miRNAs, 12 TFs and 430 target genes. The network included some well-studied candidate miRNAs and genes (e.g., miR-183-5p, FOXO1, and VEGFA) but also identified many novel objects (e.g., the miR-449 family and miR-106a-363 cluster). Here, we discuss these key regulators in the network, with the hope of providing valuable information for understanding the molecular mechanism of endometriosis.

The highly conserved miR-449 family encoded by the CDC20B gene consists of three members: miR-449a, miR-449b, and miR-449c. Accumulating evidence suggests that the miR-449 family plays an important role in cell proliferation, migration and invasion and might serve as a tumor suppressor in miRNAs [49]. Consistently, miR-449 members were absent or down-regulated in various human malignancies, such as liver cancer [50], breast cancers [51], lung cancer [52], and ovarian clear cell carcinoma [53]. In the current study, we reported the first significant down-regulation of miR-449a, miR-449b-5p, and miR-449c-5p in ectopic endometrial tissue. Considering that endometriosis shares many characteristics with malignant diseases, such as invasion and immigration, it is very likely that the lower expression of miR-449 plays an important role in the pathogenesis of endometriosis. miR-34b and miR-34c are classified into the same cluster and share related functions with miR-449 due to their structural similarities [54, 55]. Much attention has been focused on miR-34 therapeutics for cancer [56, 57]. In addition, miR-34 mimics have been used in several preclinical studies and a phase I clinical trial, and they have shown potential as anticancer therapeutics [58]. Our study identified the down-regulation of miR-34b-3p, which is consistent with the results from Teague et al. [16]. Moreover, we also identified two down-regulated miRNAs, namely, miR-34b-5p and miR-34c-5p, which have been reported by Burney et al. in the eutopic endometrium compared with normal endometrium [17]. Their results suggested that down-regulation of miR-34b and miR-34c was involved in the initiation and maintenance of endometriosis. Given the therapeutic role of miR-34 in cancer and its structural similarity to miR-449 [54], mimics from the miR-449 and miR-34b/c cluster could also be potential miRNA replacement therapies in endometriosis. We therefore suggest that further studies focus on both the miR-34 and miR-449 families by analyzing their function and target genes in endometriosis.

The miR-200 family is another important set of miRNAs that was down-regulated in the present study and has been reported to be reduced in many solid tumors [59]. The miR-200 family contains two groups, miR-200a/b/429 and miR-200c/141, that regulate a cohort of target genes involved in epithelial-mesenchymal transition, metastasis and angiogenesis [60], which are important events in the pathogenesis of endometriosis [61]. Previous studies revealed deregulated miR-200 family members in endometriotic lesions compared with eutopic endometrium [16, 18] and in eutopic endometrium compared with normal endometrium [19]. Circulating miR-200a and miR-141 have also been proposed as candidate biomarkers for endometriosis [62]. The present study is not the first to identify lower levels of miR-141-3p, miR-200a-3p, miR-200b-3p and miR-200c-3p in endometriosis. However, our report confirmed the previous results and identified a series of corresponding genes of these miRNAs, which may be involved in many physical or pathological processes in endometriosis. For example, JAZF1, which is regulated by miR-200b/c, is a proto-oncogene known to be affected by chromosomal translocation in endometrial stromal tumors [63]. In the study by Yotova et al., JAZF1 was hyper-methylated and down-regulated in stromal cells from ovarian ectopic endometria compared with normal endometria [25]. As another example, TP73, which is regulated by miR-200b, is part of the TP53 family and plays a role in the cell cycle, apoptosis and infertility [64], and no previous research has revealed its functions in endometriosis.

miR-363, miR-20b, and miR-106a are transcribed from the miR-106a-363 cluster [65]. All of these miRNAs were down-regulated in ectopic lesions according to our results. Intriguingly, these miRNAs have dual functions either as tumor suppressors or as oncogenic miRNAs. A high level of miR-363-3p suppresses the proliferation of human hepatocellular carcinoma cells [66] and decreases the metastasis of neuroblastoma cells [67], whereas the knock-down of miR-363-3p was reported to suppress carcinogenesis in gastric cancer cells [68]. Li et al. demonstrated that miR-20b acts as an oncogene by contributing to breast tumorigenesis [69], but Hong et al. found that miR-20b was markedly down-regulated and served as a tumor suppressor miRNA in papillary thyroid carcinoma [70]. Li et al. identified increased miR-106a as playing an oncogenic role in pancreatic cancer [71]. In contrast, Zhi et al. reported that miR-106a could inhibit proliferation and promote apoptosis in astrocytoma cells [72]. As no study has previously reported the functions of miR-363-3p, miR-20b and miR-106a in endometriosis, we infer that the decreased expression of these three miRNAs might have tumor suppressive-like functions in the initiation, progression, and metastasis of endometriotic lesions. Additionally, many corresponding target genes, such as ADAMTS family members, which are also known tumor suppressor genes, are down-regulated, consistent with the present study, which confirmed our inference indirectly.

miR-182-5p, which belongs to the miR-183/96/182 cluster, contains the highest node degree in the identified DEMs and shows a 5.3-fold down-regulation in EC versus EU. This finding is consistent with the results obtained from the microarray screening and qRT-PCR validation by Filigheddu et al. [18] Similar to miR-182-5p, miR-183-5p expression in ectopic endometrium was 5.7-fold down-regulated in our study and was confirmed by previously published reports [47, 73]. Additionally, functional analysis indicated that the down-regulation of miR-183 can inhibit apoptosis and enhance the invasive ability of endometrial stromal cells [47]. These findings indicate that aberrant miR-182/183 expression is part of the epigenetic mechanism of pathogenesis and the development of endometriosis.

Based on the integrative regulatory network, a series of TFs also attracted our attention, including FOX family members, GATA family members, and E2F family members, as well as CEBPA, SOX9 and HNF4A. FOX (fork-head box) proteins are a family of transcription factors involved in many pathologic and physiological processes, such as cell growth, proliferation, and longevity [74]. The FOX subfamilies, such as FOXO, FOXC and FOXP, play different roles in different types of diseases [75]. Here, we identified four important FOX subfamily members, FOXO1, FOXC1, FOXL1 and FOXJ1, which cooperatively regulate the integrative networks in endometriosis. FOXO1 serves as an anti-oncogene in various malignancies through diverse mechanisms, such as promoting apoptosis and facilitating DNA repair [76] and was decreased by 3.5-fold in endometriotic tissue in the present study. Previous data showed the same trend with lower FOXO1 levels expressed in endometriosis, which might contribute to an overactive PI3K/AKT pathway, but the underlying mechanisms are not well understood [77]. Additionally, two FBLs, miR-182-5p/FOXO1/NUAK1 and miR-182-5p/FOXO1/SCN9A, were observed and their functions are expected to be explored in the future. In recent years, an elevated expression of FOXC1 has been detected in different tumors, such as esophageal cancer [78], hepatocellular carcinoma [79] and osteosarcoma [80]. The current views on its function and molecular mechanisms reveal that FOXC1 as an oncogene plays multiple roles in cancer progression and metastasis by affecting different targets [78,79,80,81]. The present study evaluated FOXC1 expression by high-throughput sequencing and qRT-PCR, and both results indicated that FOXC1 increased approximately sevenfold in EC compared with EU. Furthermore, FOXC1 was identified as a direct target of miR-133b and shares the same targets including EYA4 and LHFP.

GATA1, GATA2 and GATA6 belong to the GATA family, which mainly plays roles in morphogenesis and organogenesis [82]. Our study observed that GATA1 could increase the miR-202-3p expression level and act on 40 targets in the endometriotic lesion. Although the cause of the high expression of GATA1 is not yet clear, Hawkins et al. reported that miR-202-3p was up-regulated in EC compared to normal endometria [19]. Moreover, we suggest that GATA2 and GATA6, as target TFs of the miR-200 family, might be involved in the pathogenesis of endometriosis together with the miR-200 family.

E2F1 is the most important TF in the E2F family, which is known to integrate cell cycle progression with DNA repair, replication, and multiple checkpoints [83]. Interestingly, the function of E2F1 in tumors is a paradox; both positive and negative effects of E2F1 on tumorigenesis have been observed [84]. The present study showed down-regulation of the E2F1 target miR-449 family members, which suggests that their combination possibly initiates endometriosis. E2F2, another important family member of E2F, is also involved in regulating the cell cycle and is associated with cancer growth [85]. A previous study reported that E2F2 can promote the development of liver cancer [86] and can also inhibit the proliferation of prostate cancer cell lines [87]. In our joint analysis, we observed that the down-regulated E2F2 is consistent with a lower expression of miR-196b-5p. Before our research, Abe et al. found that miR-196b-5p was down-regulated in ectopic endometria compared with normal endometria by microarray. Their subsequent functional experiments revealed that miR-196b can inhibit proliferation and induce apoptosis in endometriotic stromal cells [21], which led us to hypothesize that E2F2 might also play important roles in endometriosis by targeting miR-196b.

CEBPA is an important TF controlling hematopoietic differentiation and homeostasis, gene mutation and aberrant expression, which may contribute to hematological system diseases [88]. We observed that highly expressed CEBPA in EC targets 49 corresponding genes and that like FOXO1, was targeted by miR-182-5p. In a recent study, Wang et al. reported that miR-182 can repress CEBPA in hepatocellular carcinoma and is involved in the AKT signaling pathway [89]. This negative correlation between miR-182 and CEBPA is consistent with our result; however, their specific mechanism and involvement in the downstream signaling pathway in endometriosis need to be further confirmed.

Notably, the integrative data revealed that HNF4A targets 12 miRNAs and thus might be a miRNA “hub” in endometriosis. HNF4A is a critical factor in liver function and digestive diseases [90]. The mechanisms of HNF4A and its target miRNAs have not been fully elucidated yet; it is known that HNF4A can directly up-regulate some specific miRNAs to stabilize the hepatocyte phenotype [91]. This unexpected finding requires further validation to identify the underlying mechanism mediated by HNF4A in endometriosis.

GO and KEGG enrichment analysis delineated 42 functional terms and 27 important signaling pathways, most of which are consistent with the current knowledge on endometriosis. For example, the most significant term and pathway are plasma membrane (GO: 0005886) and the PI3K-Akt signaling pathway, which enriched 145 and 17 DEGs, respectively. Ectopic lesions are exposed to a unique peritoneal microenvironment that is characterized by elevated levels of hormones, inflammation, oxidative stress and iron [92]. These abnormal molecules can influence cell survival and regulate cellular functions through binding to cell membrane receptors, and then lead to the subsequent cascading activation of kinase signaling pathways [93]. Therefore, the plasma membrane term and PI3K-Akt signaling pathway are clearly relevant for endometriosis. The PI3K-Akt pathway, which integrates a variety of extracellular signals and regulates various cellular functions, including cell growth, differentiation, transformation and survival, is being investigated as a therapeutic target in other diseases and thus may also represent a target for endometriosis treatment [94, 95]. In addition, our enrichment analysis revealed some novel dysregulated signaling pathways, such as the Rap1 signaling pathway and the retrograde endocannabinoid signaling pathway, and their functions should be further investigated.

Despite our novel results, the limitations of this study should also be mentioned. First, the number of sequenced samples was relatively inadequate, and only non-hormonal treated patients with moderate to severe stage disease within the secretory phase were enrolled. Endometriosis is a complex disease that is influenced by multiple factors, such as the menstrual cycle, CA125 levels and various subtypes. These differences might impact the expression profile. Although normalized samples are helpful to determine the key pathogenic molecules, our study was most likely underpowered in discovering DEMs and DEGs with subtle expression differences due to the limited sample size. Second, the present results are preliminary and descriptive. Integrative analysis of paired miRNA and mRNA profiling data cannot entirely exclude false positive results. In addition, only some of the identified miRNAs and TFs were validated by qRT-PCR. Therefore, other experiments, such as cross-linking immunoprecipitation or functional experimental validation, are needed to holistically validate the interactions identified by our method. Furthermore, the positive findings that were not included in the network might contain valuable information and are worthy of further discussion. Third, various computational methods have been devised to measure the miRNA-mRNA regulatory relationships in recent years, such as correlation analysis, regression models, Bayesian network learning and causal inference techniques [96]. However, which one is the best method remains an open question in practical use. Here, we used the typical Pearson correlation coefficient to analyze the relationship between DEMs and DEGs, which might lead to some differences from other methods and an inevitable loss of information.

Conclusion

In summary, we provide the first characterization of miRNA-TF-gene co-expression network in paired miRNA and mRNA expression profiling of endometriosis. By establishing the regulatory network in miRNA and mRNA expression profiling, our present study confirms and significantly extends the results of prior studies, thereby defining some crucial miRNAs, TFs and genes involved in the pathogenesis of endometriosis. The identified interactions could also have implications toward targeted therapeutic strategies for endometriosis. Further in-depth functional studies are encouraged to confirm our results.

Abbreviations

BP:

Biological process

CC:

GO cellular component

DEG:

Differentially expressed gene

DEM:

Differentially expressed miRNA

EC:

Ectopic endometrium

EU:

Eutopic endometrium

FBL:

Feed-back loop

FC:

Fold change

FDR:

False discovery rate

FFL:

Feed-forward loop

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

MF:

Molecular function

miRNA:

microRNA

mRNA:

messenger RNA

PCC:

Pearson correlation coefficient

qRT-PCR:

Quantitative reverse transcription-polymerase chain reaction

r-AFS:

Revised American fertility society

RPKM:

Reads per kilobase of one gene per million reads

RPM:

Reads per million total reads

small RNA-seq:

Small RNA sequencing

TF:

Transcription factors

References

  1. Kennedy S, Bergqvist A, Chapron C, D'Hooghe T, Dunselman G, Greb R, Hummelshoj L, Prentice A, Saridogan E. Endometriosis ESIGf, Endometrium guideline development G. ESHRE guideline for the diagnosis and treatment of endometriosis. Hum Reprod. 2005;20(10):2698–704.

    Article  PubMed  Google Scholar 

  2. Stern RC, Dash R, Bentley RC, Snyder MJ, Haney AF, Robboy SJ. Malignancy in endometriosis: frequency and comparison of ovarian and extraovarian types. Int J Gynecol Pathol. 2001;20(2):133–9.

    Article  CAS  PubMed  Google Scholar 

  3. Meuleman C, Vandenabeele B, Fieuws S, Spiessens C, Timmerman D, D'Hooghe T. High prevalence of endometriosis in infertile women with normal ovulation and normospermic partners. Fertil Steril. 2009;92(1):68–74.

    Article  PubMed  Google Scholar 

  4. Montgomery GW, Nyholt DR, Zhao ZZ, Treloar SA, Painter JN, Missmer SA, Kennedy SH, Zondervan KT. The search for genes contributing to endometriosis risk. Hum Reprod Update. 2008;14(5):447–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.

    Article  CAS  PubMed  Google Scholar 

  6. Eulalio A, Huntzinger E, Izaurralde E. Getting to the root of miRNA-mediated gene silencing. Cell. 2008;132(1):9–14.

    Article  CAS  PubMed  Google Scholar 

  7. Latchman DS. Transcription factors: an overview. Int J Biochem Cell Biol. 1997;29(12):1305–12.

    Article  CAS  PubMed  Google Scholar 

  8. Vaquerizas JM, Kummerfeld SK, Teichmann SA, Luscombe NM. A census of human transcription factors: function, expression and evolution. Nat Rev Genet. 2009;10(4):252–63.

    Article  CAS  PubMed  Google Scholar 

  9. Hobert O. Gene regulation by transcription factors and microRNAs. Science. 2008;319(5871):1785–6.

    Article  CAS  PubMed  Google Scholar 

  10. Arora S, Rana R, Chhabra A, Jaiswal A, Rani V. miRNA-transcription factor interactions: a combinatorial regulation of gene expression. Mol Gen Genomics. 2013;288(3–4):77–87.

    Article  CAS  Google Scholar 

  11. Cui Q, Yu Z, Pan Y, Purisima EO, Wang E. MicroRNAs preferentially target the genes with high transcriptional regulation complexity. Biochem Biophys Res Commun. 2007;352(3):733–8.

    Article  CAS  PubMed  Google Scholar 

  12. Shalgi R, Lieber D, Oren M, Pilpel Y. Global and local architecture of the mammalian microRNA-transcription factor regulatory network. PLoS Comput Biol. 2007;3(7):e131.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Marson A, Levine SS, Cole MF, Frampton GM, Brambrink T, Johnstone S, Guenther MG, Johnston WK, Wernig M, Newman J, Calabrese JM, Dennis LM, Volkert TL, Gupta S, Love J, Hannett N, Sharp PA, Bartel DP, Jaenisch R, Young RA. Connecting microRNA genes to the core transcriptional regulatory circuitry of embryonic stem cells. Cell. 2008;134(3):521–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Sadeghi M, Ranjbar B, Ganjalikhany MR, MK F, Schmitz U, Wolkenhauer O, Gupta SK. MicroRNA and transcription factor gene regulatory network analysis reveals key regulatory elements associated with prostate cancer progression. PLoS One. 2016;11(12):e0168760.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Feng F, Wu J, Gao Z, Yu S, Cui Y. Screening the key microRNAs and transcription factors in prostate cancer based on microRNA functional synergistic relationships. Medicine (Baltimore). 2017;96(1):e5679.

    Article  CAS  Google Scholar 

  16. Ohlsson Teague EM, Van der Hoek KH, Van der Hoek MB, Perry N, Wagaarachchi P, Robertson SA, Print CG, Hull LM. MicroRNA-regulated pathways associated with endometriosis. Mol Endocrinol. 2009;23(2):265–75.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Burney RO, Hamilton AE, Aghajanova L, Vo KC, Nezhat CN, Lessey BA, Giudice LC. MicroRNA expression profiling of eutopic secretory endometrium in women with versus without endometriosis. Mol Hum Reprod. 2009;15(10):625–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Filigheddu N, Gregnanin I, Porporato PE, Surico D, Perego B, Galli L, Patrignani C, Graziani A, Surico N. Differential expression of microRNAs between eutopic and ectopic endometrium in ovarian endometriosis. J Biomed Biotechnol. 2010;2010:369549.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Hawkins SM, Creighton CJ, Han DY, Zariff A, Anderson ML, Gunaratne PH, Matzuk MM. Functional microRNA involved in endometriosis. Mol Endocrinol. 2011;25(5):821–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Laudanski P, Charkiewicz R, Kuzmicki M, Szamatowicz J, Charkiewicz A, Niklinski J. MicroRNAs expression profiling of eutopic proliferative endometrium in women with ovarian endometriosis. Reprod Biol Endocrinol. 2013;11:78.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Abe W, Nasu K, Nakada C, Kawano Y, Moriyama M, Narahara H. miR-196b targets c-myc and Bcl-2 expression, inhibits proliferation and induces apoptosis in endometriotic stromal cells. Hum Reprod. 2013;28(3):750–61.

    Article  CAS  PubMed  Google Scholar 

  22. Braza-Boils A, Mari-Alexandre J, Gilabert J, Sanchez-Izquierdo D, Espana F, Estelles A, Gilabert-Estelles J. MicroRNA expression profile in endometriosis: its relation to angiogenesis and fibrinolytic factors. Hum Reprod. 2014;29(5):978–88.

    Article  CAS  PubMed  Google Scholar 

  23. Saare M, Rekker K, Laisk-Podar T, Soritsa D, Roost AM, Simm J, Velthut-Meikas A, Samuel K, Metsalu T, Karro H, Soritsa A, Salumets A, Peters M. High-throughput sequencing approach uncovers the miRNome of peritoneal endometriotic lesions and adjacent healthy tissues. PLoS One. 2014;9(11):e112630.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Laudanski P, Charkiewicz R, Tolwinska A, Szamatowicz J, Charkiewicz A, Niklinski J. Profiling of selected MicroRNAs in proliferative Eutopic Endometrium of women with ovarian endometriosis. Biomed Res Int. 2015;2015:760698.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Yotova I, Hsu E, Do C, Gaba A, Sczabolcs M, Dekan S, Kenner L, Wenzl R, Tycko B. Epigenetic alterations affecting transcription factors and signaling pathways in Stromal cells of endometriosis. PLoS One. 2017;12(1):e0170859.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Canis M, Donnez JG, Guzick, DS, Halme JK, Rock JA, Schenken RS, Vernon MW. Revised American Society for Reproductive Medicine classification of endometriosis: 1996. Fertil Steril. 1997;67(5):817–21.

  27. Noyes RW, Hertig AT, Rock J. Dating the endometrial biopsy. Am J Obstet Gynecol. 1975;122(2):262–3.

    Article  CAS  PubMed  Google Scholar 

  28. Chomczynski P, Sacchi N. The single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction: twenty-something years on. Nat Protoc. 2006;1(2):581–5.

    Article  CAS  PubMed  Google Scholar 

  29. Li H, Wu B, Geng J, Zhou J, Zheng R, Chai J, Li F, Peng J, Jiang S. Integrated analysis of miRNA/mRNA network in placenta identifies key factors associated with labor onset of large white and Qingping sows. Sci Rep. 2015;5:13074.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Aken BL, Ayling S, Barrell D, Clarke L, Curwen V, Fairley S, Fernandez Banet J, Billis K, Garcia Giron C, Hourlier T, Howe K, Kahari A, Kokocinski F, Martin FJ, Murphy DN, Nag R, Ruffier M, Schuster M, Tang YA, Vogel JH, White S, Zadissa A, Flicek P, Searle SM. The Ensembl gene annotation system. Database (Oxford). 2016;2016

  31. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(Database issue):D68–73.

    Article  CAS  PubMed  Google Scholar 

  33. Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.

    Article  PubMed  Google Scholar 

  34. Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Anders S, Pyl PT, Huber W. HTSeq--a python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.

    Article  CAS  PubMed  Google Scholar 

  37. Bisognin A, Sales G, Coppe A, Bortoluzzi S, Romualdi C. MAGIA(2): from miRNA and genes expression data integrative analysis to microRNA-transcription factor mixed regulatory circuits (2012 update). Nucleic Acids Res. 2012;40(Web Server issue):W13–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007;27(1):91–105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Friard O, Re A, Taverna D, De Bortoli M, Cora D. CircuitsDB: a database of mixed microRNA/transcription factor feed-forward regulatory circuits in human and mouse. BMC Bioinformatics. 2010;11:435.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Wang J, Lu M, Qiu C, Cui Q. TransmiR: a transcription factor-microRNA regulation database. Nucleic Acids Res. 2010;38(Database issue):D119–22.

    Article  CAS  PubMed  Google Scholar 

  41. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  Google Scholar 

  43. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  PubMed  Google Scholar 

  44. Nunez YO, Truitt JM, Gorini G, Ponomareva ON, Blednov YA, Harris RA, Mayfield RD. Positively correlated miRNA-mRNA regulatory networks in mouse frontal cortex during early stages of alcohol dependence. BMC Genomics. 2013;14:725.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Barabasi AL, Oltvai ZN. Network biology: understanding the cell's functional organization. Nat Rev Genet. 2004;5(2):101–13.

    Article  CAS  PubMed  Google Scholar 

  46. Teague EM, Print CG, Hull ML. The role of microRNAs in endometriosis and associated reproductive conditions. Hum Reprod Update. 2010;16(2):142–65.

    Article  PubMed  Google Scholar 

  47. Shi XY, Gu L, Chen J, Guo XR, Shi YL. Downregulation of miR-183 inhibits apoptosis and enhances the invasive potential of endometrial stromal cells in endometriosis. Int J Mol Med. 2014;33(1):59–67.

    Article  CAS  PubMed  Google Scholar 

  48. Rajewsky N. microRNA target predictions in animals. Nat Genet. 2006;38(Suppl):S8–13.

    Article  CAS  PubMed  Google Scholar 

  49. Bou Kheir T, Futoma-Kazmierczak E, Jacobsen A, Krogh A, Bardram L, Hother C, Gronbaek K, Federspiel B, Lund AH, Friis-Hansen L. miR-449 inhibits cell proliferation and is down-regulated in gastric cancer. Mol Cancer. 2011;10:29.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Sandbothe M, Buurman R, Reich N, Greiwe L, Vajen B, Gurlevik E, Schaffer V, Eilers M, Kuhnel F, Vaquero A, Longerich T, Roessler S, Schirmacher P, Manns MP, Illig T, Schlegelberger B, Skawran B. The microRNA-449 family inhibits TGF-beta-mediated liver cancer cell migration by targeting SOX4. J Hepatol. 2017;

  51. Zhang Z, Wang J, Gao R, Yang X, Zhang Y, Li J, Zhang J, Zhao X, Xi C, Lu X. Downregulation of microRNA-449 promotes migration and invasion of breast cancer cells by targeting tumor protein D52 (TPD52). Oncol Res. 2016;

  52. Luo W, Huang B, Li Z, Li H, Sun L, Zhang Q, Qiu X, Wang E. MicroRNA-449a is downregulated in non-small cell lung cancer and inhibits migration and invasion by targeting c-met. PLoS One. 2013;8(5):e64759.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Jang SG, Yoo CW, Park SY, Kang S, Kim HK. Low expression of miR-449 in gynecologic clear cell carcinoma. Int J Gynecol Cancer. 2014;24(9):1558–63.

    Article  PubMed  Google Scholar 

  54. Wu J, Bao J, Kim M, Yuan S, Tang C, Zheng H, Mastick GS, Xu C, Yan W. Two miRNA clusters, miR-34b/c and miR-449, are essential for normal brain development, motile ciliogenesis, and spermatogenesis. Proc Natl Acad Sci U S A. 2014;111(28):E2851–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Hermeking H. The miR-34 family in cancer and apoptosis. Cell Death Differ. 2010;17(2):193–9.

    Article  CAS  PubMed  Google Scholar 

  56. Wang R, Ma J, Wu Q, Xia J, Miele L, Sarkar FH, Wang Z. Functional role of miR-34 family in human cancer. Curr Drug Targets. 2013;14(10):1185–91.

    Article  CAS  PubMed  Google Scholar 

  57. Misso G, Di Martino MT, De Rosa G, Farooqi AA, Lombardi A, Campani V, Zarone MR, Gulla A, Tagliaferri P, Tassone P, Caraglia M. Mir-34: a new weapon against cancer? Mol Ther Nucleic Acids. 2014;3:e194.

    Article  CAS  PubMed  Google Scholar 

  58. Rupaimoole R, Slack FJ. MicroRNA therapeutics: towards a new era for the management of cancer and other diseases. Nat Rev Drug Discov. 2017;

  59. Feng X, Wang Z, Fillmore R, Xi Y. MiR-200, a new star miRNA in human cancer. Cancer Lett. 2014;344(2):166–73.

    Article  CAS  PubMed  Google Scholar 

  60. Matsuzaki S, Darcha C. Epithelial to mesenchymal transition-like and mesenchymal to epithelial transition-like processes might be involved in the pathogenesis of pelvic endometriosis. Hum Reprod. 2012;27(3):712–21.

    Article  CAS  PubMed  Google Scholar 

  61. Proestling K, Birner P, Gamperl S, Nirtl N, Marton E, Yerlikaya G, Wenzl R, Streubel B, Husslein H. Enhanced epithelial to mesenchymal transition (EMT) and upregulated MYC in ectopic lesions contribute independently to endometriosis. Reprod Biol Endocrinol. 2015;13:75.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Rekker K, Saare M, Roost AM, Kaart T, Soritsa D, Karro H, Soritsa A, Simon C, Salumets A, Peters M. Circulating miR-200-family micro-RNAs have altered plasma levels in patients with endometriosis and vary with blood collection time. Fertil Steril. 2015;104(4):938–46. e2

    Article  CAS  PubMed  Google Scholar 

  63. Li H, Wang J, Mor G, Sklar J. A neoplastic gene fusion mimics trans-splicing of RNAs in normal human cells. Science. 2008;321(5894):1357–61.

    Article  CAS  PubMed  Google Scholar 

  64. Vagnini LD, Renzi A, Oliveira-Pelegrin GR, Canas Mdo C, Petersen CG, Mauri AL, Oliveira JB, Baruffi RL, Cavagna M, Franco Junior JG. The TP73 gene polymorphism (rs4648551, a>G) is associated with diminished ovarian reserve. PLoS One. 2015;10(3):e0120048.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Mendell JT. miRiad roles for the miR-17-92 cluster in development and disease. Cell. 2008;133(2):217–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Han H, Sun D, Li W, Shen H, Zhu Y, Li C, Chen Y, Lu L, Li W, Zhang J, Tian Y, Li Y. A c-Myc-MicroRNA functional feedback loop affects hepatocarcinogenesis. Hepatology. 2013;57(6):2378–89.

    Article  CAS  PubMed  Google Scholar 

  67. Qiao J, Lee S, Paul P, Theiss L, Tiao J, Qiao L, Kong A, Chung DH. miR-335 and miR-363 regulation of neuroblastoma tumorigenesis and metastasis. Surgery. 2013;154(2):226–33.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Hsu KW, Wang AM, Ping YH, Huang KH, Huang TT, Lee HC, Lo SS, Chi CW, Yeh TS. Downregulation of tumor suppressor MBP-1 by microRNA-363 in gastric carcinogenesis. Carcinogenesis. 2014;35(1):208–17.

    Article  CAS  PubMed  Google Scholar 

  69. Li D, Ilnytskyy Y, Kovalchuk A, Khachigian LM, Bronson RT, Wang B, Kovalchuk O. Crucial role for early growth response-1 in the transcriptional regulation of miR-20b in breast cancer. Oncotarget. 2013;4(9):1373–87.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Hong S, Yu S, Li J, Yin Y, Liu Y, Zhang Q, Guan H, Li Y, Xiao H. MiR-20b displays tumor-suppressor functions in papillary thyroid carcinoma by regulating the MAPK/ERK signaling pathway. Thyroid. 2016;26(12):1733–43.

    Article  CAS  PubMed  Google Scholar 

  71. Li P, Xu Q, Zhang D, Li X, Han L, Lei J, Duan W, Ma Q, Wu Z, Wang Z. Upregulated miR-106a plays an oncogenic role in pancreatic cancer. FEBS Lett. 2014;588(5):705–12.

    Article  CAS  PubMed  Google Scholar 

  72. Zhi F, Zhou G, Shao N, Xia X, Shi Y, Wang Q, Zhang Y, Wang R, Xue L, Wang S, Wu S, Peng Y, Yang Y. miR-106a-5p inhibits the proliferation and migration of astrocytoma cells and promotes apoptosis by targeting FASTK. PLoS One. 2013;8(8):e72390.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Chen J, Gu L, Ni J, Hu P, Hu K, Shi YL. MiR-183 regulates ITGB1P expression and promotes invasion of endometrial Stromal cells. Biomed Res Int. 2015;2015:340218.

    PubMed  PubMed Central  Google Scholar 

  74. Zaret KS, Carroll JS. Pioneer transcription factors: establishing competence for gene expression. Genes Dev. 2011;25(21):2227–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Hannenhalli S, Kaestner KH. The evolution of fox genes and their role in development and disease. Nat Rev Genet. 2009;10(4):233–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Fu Z, Tindall DJ. FOXOs, cancer and regulation of apoptosis. Oncogene. 2008;27(16):2312–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Yin X, Pavone ME, Lu Z, Wei J, Kim JJ. Increased activation of the PI3K/AKT pathway compromises decidualization of stromal cells from endometriosis. J Clin Endocrinol Metab. 2012;97(1):E35–43.

    Article  CAS  PubMed  Google Scholar 

  78. Myatt SS, Wang J, Monteiro LJ, Christian M, Ho KK, Fusi L, Dina RE, Brosens JJ, Ghaem-Maghami S, Lam EW. Definition of microRNAs that repress expression of the tumor suppressor gene FOXO1 in endometrial cancer. Cancer Res. 2010;70(1):367–77.

    Article  CAS  PubMed  Google Scholar 

  79. Zhu X, Wei L, Bai Y, Wu S, Han S. FoxC1 promotes epithelial-mesenchymal transition through PBX1 dependent transactivation of ZEB2 in esophageal cancer. Am J Cancer Res. 2017;7(8):1642–53.

    PubMed  PubMed Central  Google Scholar 

  80. Xia L, Huang W, Tian D, Zhu H, Qi X, Chen Z, Zhang Y, Hu H, Fan D, Nie Y, Wu K. Overexpression of forkhead box C1 promotes tumor metastasis and indicates poor prognosis in hepatocellular carcinoma. Hepatology. 2013;57(2):610–24.

    Article  CAS  PubMed  Google Scholar 

  81. Deng L, Liu T, Zhang B, Wu H, Zhao J, Chen J. Forkhead box C1 is targeted by microRNA-133b and promotes cell proliferation and migration in osteosarcoma. Exp Ther Med. 2017;14(4):2823–30.

    Article  PubMed  PubMed Central  Google Scholar 

  82. Myatt SS, Lam EW. The emerging roles of forkhead box (fox) proteins in cancer. Nat Rev Cancer. 2007;7(11):847–59.

    Article  CAS  PubMed  Google Scholar 

  83. Tremblay JJ, Viger RS. GATA factors differentially activate multiple gonadal promoters through conserved GATA regulatory elements. Endocrinology. 2001;142(3):977–86.

    Article  CAS  PubMed  Google Scholar 

  84. Ren B, Cam H, Takahashi Y, Volkert T, Terragni J, Young RA, Dynlacht BD. E2F integrates cell cycle progression with DNA repair, replication, and G(2)/M checkpoints. Genes Dev. 2002;16(2):245–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Johnson DG. The paradox of E2F1: oncogene and tumor suppressor gene. Mol Carcinog. 2000;27(3):151–7.

    Article  CAS  PubMed  Google Scholar 

  86. Fujiwara K, Yuwanita I, Hollern DP, Andrechek ER. Prediction and genetic demonstration of a role for activator E2Fs in Myc-induced tumors. Cancer Res. 2011;71(5):1924–32.

    Article  CAS  PubMed  Google Scholar 

  87. Zhan L, Huang C, Meng XM, Song Y, XQ W, Miu CG, Zhan XS, Li J. Promising roles of mammalian E2Fs in hepatocellular carcinoma. Cell Signal. 2014;26(5):1075–81.

    Article  CAS  PubMed  Google Scholar 

  88. Dong Q, Meng P, Wang T, Qin W, Qin W, Wang F, Yuan J, Chen Z, Yang A, Wang H. MicroRNA let-7a inhibits proliferation of human prostate cancer cells in vitro and in vivo by targeting E2F2 and CCND2. PLoS One. 2010;5(4):e10147.

    Article  PubMed  PubMed Central  Google Scholar 

  89. Zhang P, Iwasaki-Arai J, Iwasaki H, Fenyus ML, Dayaram T, Owens BM, Shigematsu H, Levantini E, Huettner CS, Lekstrom-Himes JA, Akashi K, Tenen DG. Enhancement of hematopoietic stem cell repopulating capacity and self-renewal in the absence of the transcription factor C/EBP alpha. Immunity. 2004;21(6):853–63.

    Article  CAS  PubMed  Google Scholar 

  90. Wang C, Ren R, Hu H, Tan C, Han M, Wang X, Zheng Y. MiR-182 is up-regulated and targeting Cebpa in hepatocellular carcinoma. Chin J Cancer Res. 2014;26(1):17–29.

    CAS  PubMed  PubMed Central  Google Scholar 

  91. Weltmeier F, Borlak J. A high resolution genome-wide scan of HNF4alpha recognition sites infers a regulatory gene network in colon cancer. PLoS One. 2011;6(7):e21667.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Cicchini C, Amicone L, Alonzi T, Marchetti A, Mancone C, Tripodi M. Molecular mechanisms controlling the phenotype and the EMT/MET dynamics of hepatocyte. Liver Int. 2015;35(2):302–10.

    Article  PubMed  Google Scholar 

  93. Goud PT, Goud AP, Joshi N, Puscheck E, Diamond MP, Abu-Soud HM. Dynamics of nitric oxide, altered follicular microenvironment, and oocyte quality in women with endometriosis. Fertil Steril. 2014;102(1):151–9. e5

    Article  CAS  PubMed  Google Scholar 

  94. Little AS, Smith PD, Cook SJ. Mechanisms of acquired resistance to ERK1/2 pathway inhibitors. Oncogene. 2013;32(10):1207–15.

    Article  CAS  PubMed  Google Scholar 

  95. McKinnon BD, Kocbek V, Nirgianakis K, Bersinger NA, Mueller MD. Kinase signalling pathways in endometriosis: potential targets for non-hormonal therapeutics. Hum Reprod Update. 2016;22(3)

  96. Le TD, Zhang J, Liu L, Li J. Ensemble methods for MiRNA target prediction from expression data. PLoS One. 2015;10(6):e0131627.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank all patients for their participation.

Funding

This study was funded by the National Natural Science Foundation of China (grant number: 81571411).

Availability of data and materials

The primary sequencing datasets are available from the GEO database (GSE105764 for mRNA sequencing and GSE105765 for small RNA sequencing). The detailed data after analyzed are included in this published article and its additional files.

Author information

Author notes

    Authors and Affiliations

    1. Department of Gynecology and Obstetrics, People’s Liberation Army (PLA) Medical School, Chinese PLA General Hospital, Beijing, 100853, China

      Luyang Zhao, Chenglei Gu, Mingxia Ye, Zhe Zhang, Li’an Li, Wensheng Fan & Yuanguang Meng

    2. Department of Gynecology and Obstetrics, Peking University People’s Hospital, Beijing, China

      Luyang Zhao

    3. Department of Gynecology and Obstetrics, the 309th Hospital of Chinese PLA, Beijing, China

      Chenglei Gu

    Authors
    1. Luyang Zhao
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Chenglei Gu
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Mingxia Ye
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Zhe Zhang
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Li’an Li
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Wensheng Fan
      View author publications

      You can also search for this author in PubMed Google Scholar

    7. Yuanguang Meng
      View author publications

      You can also search for this author in PubMed Google Scholar

    Contributions

    WSF and YGM conceived and designed the experiments. LYZ, CLG, MXY and ZZ performed the experiments. LYZ, CLG and LLA analyzed the data. LYZ and CLG wrote the paper. WSF revised the manuscript for important intellectual content. YGM supervised the manuscript writing. All authors read and approved the final manuscript.

    Corresponding authors

    Correspondence to Wensheng Fan or Yuanguang Meng.

    Ethics declarations

    Ethics approval and consent to participate

    The study protocol was approved by the Local Ethical Committee of Chinese PLA General Hospital (No. 20140098).

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Publisher’s Note

    Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

    Additional files

    Additional file 1:

    Clinical characteristics of 30 enrolled patients with ovarian endometriosis. (DOCX 19 kb)

    Additional file 2:

    Supplementary methods introduction for small RNA-seq and mRNA seq. (DOCX 21 kb)

    Additional file 3:

    Primers had been used in qRT-PCR (DOCX 14 kb)

    Additional file 4:

    Summary of the Small RNA sequencing data after filtering and mapping (DOCX 15 kb)

    Additional file 5:

    List of 107 differentially expressed miRNAs in ectopic endometria compared with eutopic endometria (DOCX 17 kb)

    Additional file 6:

    Summary of the mRNA sequencing data after filtering and mapping (DOCX 16 kb)

    Additional file 7:

    Top 25 up-regulated and down-regulated mRNAs in ectopic endometria compared with paired eutopic endometria in ovarian endometriosis. (DOCX 15 kb)

    Rights and permissions

    Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

    Reprints and permissions

    About this article

    Check for updates. Verify currency and authenticity via CrossMark

    Cite this article

    Zhao, L., Gu, C., Ye, M. et al. Integration analysis of microRNA and mRNA paired expression profiling identifies deregulated microRNA-transcription factor-gene regulatory networks in ovarian endometriosis. Reprod Biol Endocrinol 16, 4 (2018). https://doi.org/10.1186/s12958-017-0319-5

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/s12958-017-0319-5

    Keywords

    Reproductive Biology and Endocrinology

    ISSN: 1477-7827

    Contact us