Genome-wide profiling of miRNA expression patterns in tubal endometriosis

in Reproduction
Authors:
Hang Qi Department of Obstetrics and Gynecology, International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
Institute of Embryo-Fetal Original Adult Disease, Affiliated to School of Medicine, Shanghai Jiaotong University, Shanghai, China

Search for other papers by Hang Qi in
Current site
Google Scholar
PubMed
Close
,
Guiling Liang Department of Obstetrics and Gynecology, International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
Institute of Embryo-Fetal Original Adult Disease, Affiliated to School of Medicine, Shanghai Jiaotong University, Shanghai, China

Search for other papers by Guiling Liang in
Current site
Google Scholar
PubMed
Close
,
Jin Yu Department of Obstetrics and Gynecology, International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
Institute of Embryo-Fetal Original Adult Disease, Affiliated to School of Medicine, Shanghai Jiaotong University, Shanghai, China

Search for other papers by Jin Yu in
Current site
Google Scholar
PubMed
Close
,
Xiaofeng Wang Department of Obstetrics and Gynecology, International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
Institute of Embryo-Fetal Original Adult Disease, Affiliated to School of Medicine, Shanghai Jiaotong University, Shanghai, China

Search for other papers by Xiaofeng Wang in
Current site
Google Scholar
PubMed
Close
,
Yan Liang Department of Obstetrics and Gynecology, International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
Institute of Embryo-Fetal Original Adult Disease, Affiliated to School of Medicine, Shanghai Jiaotong University, Shanghai, China

Search for other papers by Yan Liang in
Current site
Google Scholar
PubMed
Close
,
Xiaoqing He Department of Obstetrics and Gynecology, International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
Institute of Embryo-Fetal Original Adult Disease, Affiliated to School of Medicine, Shanghai Jiaotong University, Shanghai, China

Search for other papers by Xiaoqing He in
Current site
Google Scholar
PubMed
Close
,
Tienan Feng Clinical Research Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China

Search for other papers by Tienan Feng in
Current site
Google Scholar
PubMed
Close
, and
Jian Zhang Department of Obstetrics and Gynecology, International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiaotong University, Shanghai, China
Institute of Embryo-Fetal Original Adult Disease, Affiliated to School of Medicine, Shanghai Jiaotong University, Shanghai, China

Search for other papers by Jian Zhang in
Current site
Google Scholar
PubMed
Close

Correspondence should be addressed to J Zhang; Email: zhangjian_ipmch@sjtu.edu.cn

*(H Qi and G Liang contributed equally to this work and are co-first authors)

Open access

Sign up for journal news

MicroRNA (miRNA) expression profiles in tubal endometriosis (EM) are still poorly understood. In this study, we analyzed the differential expression of miRNAs and the related gene networks and signaling pathways in tubal EM. Four tubal epithelium samples from tubal EM patients and five normal tubal epithelium samples from uterine leiomyoma patients were collected for miRNA microarray. Bioinformatics analyses, including Ingenuity Pathway Analysis (IPA), Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, were performed. Quantitative real-time polymerase chain reaction (qRT-PCR) validation of five miRNAs was performed in six tubal epithelium samples from tubal EM and six from control. A total of 17 significantly differentially expressed miRNAs and 4343 potential miRNA-target genes involved in tubal EM were identified (fold change >1.5 and FDR-adjusted P value <0.05). IPA indicated connections between miRNAs, target genes and other gynecological diseases like endometrial carcinoma. GO and KEGG analysis revealed that most of the identified genes were involved in the mTOR signaling pathway, SNARE interactions in vesicular transport and endocytosis. We constructed an miRNA-gene-disease network using target gene prediction. Functional analysis showed that the mTOR pathway was connected closely to tubal EM. Our results demonstrate for the first time the differentially expressed miRNAs and the related signal pathways involved in the pathogenesis of tubal EM which contribute to elucidating the pathogenic mechanism of tubal EM-related infertility.

Abstract

MicroRNA (miRNA) expression profiles in tubal endometriosis (EM) are still poorly understood. In this study, we analyzed the differential expression of miRNAs and the related gene networks and signaling pathways in tubal EM. Four tubal epithelium samples from tubal EM patients and five normal tubal epithelium samples from uterine leiomyoma patients were collected for miRNA microarray. Bioinformatics analyses, including Ingenuity Pathway Analysis (IPA), Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, were performed. Quantitative real-time polymerase chain reaction (qRT-PCR) validation of five miRNAs was performed in six tubal epithelium samples from tubal EM and six from control. A total of 17 significantly differentially expressed miRNAs and 4343 potential miRNA-target genes involved in tubal EM were identified (fold change >1.5 and FDR-adjusted P value <0.05). IPA indicated connections between miRNAs, target genes and other gynecological diseases like endometrial carcinoma. GO and KEGG analysis revealed that most of the identified genes were involved in the mTOR signaling pathway, SNARE interactions in vesicular transport and endocytosis. We constructed an miRNA-gene-disease network using target gene prediction. Functional analysis showed that the mTOR pathway was connected closely to tubal EM. Our results demonstrate for the first time the differentially expressed miRNAs and the related signal pathways involved in the pathogenesis of tubal EM which contribute to elucidating the pathogenic mechanism of tubal EM-related infertility.

Introduction

Endometriosis (EM) is the presence of functional endometrial glands and stroma outside the uterine cavity. EM is seen in 5–10% of the general female population, with most affected individuals being of the reproductive age (Woodward et al. 2001). Patients may be asymptomatic or may present with chronic pelvic pain, infertility and dyspareunia (Woodward et al. 2001).

Tubal EM refers to functional endometrial glands and stromal implant involving the fallopian tube and is a part of pelvic EM. Only a few studies have focused on tubal EM. The incidence of tubal EM varies widely and accounts for about 0.29–4% of pelvic EM. (Xue & Zhang 2017). In an earlier study, EM lesions were found in the fallopian tubes of 60% (21/35) EM patients; therefore, tubal involvement in moderate and severe EM is common (Xia et al. 2018). Moreover, EM with tubal EM resulted in lower ciliary beat frequency, lower ciliated cell percentage, weaker muscular contractile activity and lower contraction frequency (Xia et al. 2018). Thus, tubal EM may cause infertility in women by affecting the function of the fallopian tube (Xia et al. 2018).

miRNAs are short non-coding RNA molecules (22–25 nucleotides long), which are transcribed in the nucleus as primary miRNA and undergo a series of maturation steps utilizing the sequential action of the endonucleases, Drosha and Dicer, to attain functional competence (Taganov et al. 2007). miRNAs have diverse functions, including morphogenesis, tissue maintenance, cell growth, differentiation, apoptosis and metabolism (Ouellet et al. 2006, Song & Tuan 2006). Aberrant miRNA expression is associated with a number of human diseases (Singh et al. 2008), including benign gynecological conditions, gynecological malignancies and fertility disorders of the human female reproductive tract. Furthermore, miRNAs can impose epigenetic effects on the pathogenesis of EM (Zondervan et al. 2018). Previous studies have identified the miRNA expression profiles of eutopic and ectopic endometrium from women with EM (Ohlsson Teague et al. 2009) and differentially expressed miRNAs between women with EM and those without EM (Jia et al. 2013, Okamoto et al. 2015). By performing bioinformatics analyses on miRNA and mRNA microarray data, novel miRNA-associated molecular pathways that are likely to play an important role in the pathophysiology of this disorder were identified (Hull et al. 2008, Ohlsson Teague et al. 2009). Additionally, animal models of EM and microarray methods have identified numerous individual EM-associated transcripts that encode the proteins involved in broad signaling pathways. These pathways mediate inflammation, tissue remodeling, apoptosis, cellular proliferation, angiogenesis and wound healing in endometriotic tissues (Eyster et al. 2007, Flores et al. 2007, Hever et al. 2007, Borghese et al. 2008, Hull et al. 2008). Nevertheless, little is known about the expression profiles and functions of miRNAs in tubal EM.

In this study, the miRNA expression profiles in tubal EM were examined by microarray analysis. Differentially expressed miRNAs between tubal EM and controls were identified, followed by the identification of target genes using an integrated method. These differentially expressed miRNAs may play a potential role in the pathogenesis and development of tubal EM.

Materials and methods

Clinical tissue specimens

Seven premenopausal patients with tubal EM (case group) and seven premenopausal patients suffering from uterine leiomyoma (control group) were recruited. Tubal EM and normal fallopian tubes in the control group were confirmed by histopathological examinations and were reviewed by at least two experienced pathologists. This study was approved by the Institutional Ethics Committee of the International Peace Maternity and Child Health Hospital in Shanghai, China (No. GKLW 2015-34) and written informed consents was obtained from all participants. From June 2016 to March 2017, tubal epithelium specimens were obtained from surgically resected fallopian tubes. Specimens were immediately preserved in liquid nitrogen and were transported from operation room to the lab for experiments within 10 min. Patient demographics and clinical pathological data were also collected.

RNA extraction and quality control

Total RNA was extracted and purified using mirVana™ miRNA Isolation Kit without phenol (Cat # AM1561, Ambion), following the manufacturer’s instructions and checked for RNA integrity by RNA integrity number (RIN) using an Agilent Bioanalyzer 2100 (Agilent Technologies). The initial sample for chip experiment using total RNA was examined using a NanoDrop nd-2000 spectrophotometer and an Agilent Bioanalyzer 2100 (Agilent Technologies).

MiRNA microarray

miRNA in total RNA was labeled by miRNA Complete Labeling and Hyb Kit (Cat # 5190-0456, Agilent Technologies) following the manufacturer’s instructions. Each slide was hybridized with 100 ng Cy3-labeled RNA in a hybridization oven (Cat # G2545A, Agilent Technologies) at 55°C and 0.06  g for 20 h. The Agilent miRNA Complete Labeling and Hyb Kit provides all the needed components for sample labeling and hybridization preparation. These include calf intestinal alkaline phosphatase (CIP), 10× calf intestinal phosphatase buffer, T4 RNA ligase, 10× T4 RNA ligase buffer, DMSO, nuclease-free water, cyanine 3-pCp, 10× GE-blocking agent and 2X Hi-RPM hybridization buffer. After hybridization, slides were washed in staining dishes (Cat # 121, Thermos Shandon, Waltham, MA, USA) with Gene-Expression Wash Buffer Kit (Cat #5188-5327, Agilent Technologies).

The slides were then scanned using an Agilent Microarray Scanner (Cat # G2565CA, Agilent Technologies) and the data were extracted using Agilent Feature Extraction Software (Agilent Technologies) with default settings.

Bio-informatic analyses to predict miRNA-target genes and network

To further investigate the functional roles of miRNAs, putative targets of miRNAs were predicted by TargetScan (http://www.targetscan.org/), and for more accuracy, the Ingenuity Pathway Analysis software (IPA software, Qiagen) was also used for the prediction. Besides, IPA identified the miRNA-target gene-disease networks, which connected tubal EM with other similar diseases. Genes were then subjected to enrichment analysis of Gene Ontology (GO) functions, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and DAVID (https://david.ncifcrf.gov). GO analysis was performed to explore the functional roles of target genes in terms of biological processes (BP) and molecular functions (MFs). Pathway analysis is a functional analysis for mapping genes to KEGG. The FDR-adjusted P value denotes the significance of the pathway correlated to the conditions. IPA showed how miRNAs participate in the mTOR signaling pathway.

Quantitative real-time PCR

The expression level of five miRNAs was measured by qRT-PCR in six paired tubal epitheliums (six tubal EM and six control women). Total RNA was isolated using TRIzol (Invitrogen), and cDNA was prepared using the bulge-loop™ miRNA qRT-PCR Starter kit. miRNA qPCR analysis was conducted to determine the miRNA expression using the SYBR Green Mix system (RiboBio, Guangzhou, China) in 384-well plates. The amplification thermal cycling conditions were an initial 95°C for 10 min followed by 40 cycles of 95°C for 2 s, 60°C for 20 s and 70°C for 10 s. After PCR, a dissociation curve was constructed at 70°C for 5 s for detection of PCR product specificity. The bulge-loop™ miRNA qRT-PCR Primer Sets (one RT primer and a pair of qPCR primers for each set) specific for miR-21-3p, miR-223-3p, miR-3663-3p, miR-450a-5p and miR-503-5p were designed by RiboBio. The relative expression level of miRNAs was normalized to U6 RNA and calculated using the comparative Ct method (2−∆∆CT).

Statistical analysis

Raw data were normalized using the Quantile algorithm, included in the R package AgiMicroRna (López-Romero, P. BMC Genomics. 2011). The obtained array data were normalized to external and internal controls to remove system-related variations. The statistical significance of the difference was estimated using the Student ’s t-test, fold change >1.5 and FDR-adjusted P values <0.05 were considered statistically significant. RT-qPCR data were analyzed using GraphPad Prism 7 (GraphPad Software Inc.).

Results

Participant demographics and clinical pathological data

After initial sample examinations, nine qualified samples (four patients and five controls) were identified for subsequent miRNA chip experiments. The participants’ demographics and clinical pathological data are summarized in Table 1. N302837-BL, N129922-BL, N305896-BL and N308632-BL were included in the case group and N305643-DZ, N305044-DZ, N306224-DZ, N303534-DZ and N306515-DZ were included in the control group. The participants varied in age from 42 to 47 years and were in secretory phase of their menstrual cycles, as estimated by the number of days since their last menstrual cycle and the endometrial histology. The P value of unpaired Student’s t-test for age was 0.53 indicating that there was no significant difference in age between the case group and the control group.

Table 1

Demographics and clinical pathological data of participants.

Sample number Age Reproductive history (TPAL) Tubal EM Pathologic diagnosis Menstrual cycle
N302837-BL 43 0-0-0-0 Diffuse adenomyosis of the uterus and tubal endometriosis Secretory phase
N129922-BL 47 0-0-0-0 Endometriotic cysts of the ovary and tubal endometriosis Secretory phase
N305896-BL 45 1-0-1-1 Endometriotic cysts of the ovary and tubal endometriosis Secretory phase
N308632-BL 42 1-0-0-1 Endometriotic cysts of the ovary and tubal endometriosis Secretory phase
N305643-DZ 43 1-0-2-1 × Uterine leiomyoma Secretory phase
N305044-DZ 45 1-0-1-1 × Uterine leiomyoma Secretory phase
N306224-DZ 46 2-0-1-2 × Uterine leiomyoma Secretory phase
N303534-DZ 45 1-0-1-4 × Uterine leiomyoma Secretory phase
N306515-DZ 46 1-0-1-1 × Uterine leiomyoma Secretory phase

TPAL, term infants, preterm infants, abortuses, living children.

Normalizing the raw data

Normalization is an essential step in microarray data preprocessing. It aims to remove non-biologic effects resulting from the experimental process so that biologic effects can be accurately identified (Qin & Satagopan 2009). After homogenization and logarithmic transformation of the original data from single fluorescence chip, boxplots were used to evaluate the data to visually compare overall signal distribution between multiple samples (Fig. 1A). The overall distribution and dispersion of the hybridization signal value data of each chip were estimated using boxplots. The boxplots showed that the distributions of miRNAs for the nine samples were similar after normalization (Fig. 1A).

Figure 1
Figure 1

(A) Box plots showing the distribution of miRNAs in the nine samples. The majority of the miRNAs on the chip were not expressed in the sample or the expression level was very low. The signal is focused on a relatively flat casing, and a small number of miRNAs with a relatively high signal are located in the upper third quartile. (B) Hierarchical clustering of miRNAs. The samples are in columns and miRNAs in rows. BL represents case group and DZ represents the control group. The red columns represent upregulation in the case group, and the blue columns represent downregulation in the case group. The brightness of the color indicates the degree of upregulation or downregulation. (C, D, E, F and G) Verification of five differentially expressed miRNAs related to tubal EM identified by miRNA-microarray. The expression levels of the five differentially expressed miRNAs in tubal epithelium of tubal EM (n = 6) and control women (n = 6). (C, D and E) are upregulated, and F-G are downregulated in transcription analysis and qRT-PCR. Data are presented as mean ± s.e. P value was determined by unpaired Student’s t-test. *P < 0.05, **0.01 < P < 0.05 and ***P < 0.01, compared to the corresponding controls.

Citation: Reproduction 157, 6; 10.1530/REP-18-0631

Identification of differential miRNA profiles

A total of 2532 miRNAs were scanned during microarray analysis. After normalization of the raw data, fold change and Student’s t-test were used to screen the differentially expressed genes; the selection criteria were as follows: fold change >1.5 and FDR-adjusted P value <0.05; finally, 17 miRNAs were found significantly differentially expressed including 4 upregulated and 13 downregulated miRNAs in tubal EM (shown in Supplementary Table 1, see section on supplementary data given at the end of this article). To be more intuitive, differentially expressed miRNAs with statistical significance between the two groups were identified using hierarchical clustering (Fig. 1B).

Validation of differentially expressed miRNAs with qPCR

To confirm the miRNA–microarray results, three upregulated miRNAs and two downregulated miRNAs with fold change >2 were identified with qRT-PCR. The results of qRT-PCR (Fig. 1C, D, E, F and G) showed that the expression of miRNAs was consistent with the miRNA–microarray results.

IPA of miRNAs

Based on an analysis of predicted gene targets using the IPA software, one upregulated miRNA hsa-miR-223-3p and two downregulated miRNAs hsa-miR-450a -5p and hsa-miR-503-5p were predicted to top 70 target genes involved in the pathogenesis of EM. These miRNAs were selected because their expression profiling showed a fold change >2 which was confirmed by qRT-PCR. Besides, IPA showed, that among the miR-223-3p target genes, 15 were associated with endometrioid carcinoma, 17 with female genital tract adenocarcinoma and 22 with endometrial carcinoma (Fig. 2A). Among the miR-450a-5p target genes, 22 were associated with endometrioid carcinoma, 22 with endometrial carcinoma, 24 with female genital tract adenocarcinoma, 27 with female genital tract cancer and 27 with tumorigenesis of the reproductive tract (Fig. 2B). Among the miR-503-5p target genes, 41 were associated with endometrioid carcinoma, 43 with endometrial carcinoma, 46 with female genital tract adenocarcinoma, 52 with tumorigenesis of the reproductive tract and 52 were associated with female genital tract cancer (Fig. 2B). More details of the miRNA–gene disease network are shown in Supplementary Table 2.

Figure 2
Figure 2

Inverse relationship of differentially expressed miRNAs with target genes and gynecological diseases. The central blue comb-like sign represents miRNAs, the red rectangular, triangular, trapezoidal and oval boxes, represent their target genes; the blue cross-shaped boxes represent gynecological diseases and the discontinuous lines connecting them represent the interaction. (A) Network of miR-223-3P-downstream gene-gynecological diseases. (B) Network of miR-450a-5p, miR-503-5p-downstream gene-gynecological diseases.

Citation: Reproduction 157, 6; 10.1530/REP-18-0631

Integrating miRNAs to functional pathways

TargetScan predicted 4343 target genes on which we performed functional annotation including GO and KEGG analyses. The results of the GO analysis of these genes according to the values in the enrichment score under BP and MF are shown in Fig. 3A. In BP, the top two processes were regulation of nuclear-transcribed mRNA poly(A) tail shortening and positive regulation of nuclear-transcribed mRNA poly(A) tail shortening. In MF, there were type-2 fibroblast growth factor receptor binding and copper ion transmembrane transporter activity. Other GO categories were also enriched such as regulation of nuclear-transcribed mRNA catabolic process and positive regulation of nuclear-transcribed mRNA catabolic process (Fig. 3A). KEGG enrichment analysis of differentially expressed genes can enrich pathways which are significantly different and is suitable for finding biological regulatory pathways with significant differences under experimental conditions. Count numbers higher than two and FDR-adjusted P value less than 0.05 were chosen as the cut-off criterion. Five hundred and ten KEGG pathways were identified in all the genes (data not shown). The top 30 KEGG pathways are shown in Fig. 3B. KEGG results showed that pathways in cancer were associated with 85 genes, which was the largest association. The interesting candidate pathways included the MAPK signaling pathway and Wnt signaling pathway, which has been reported to be involved in the pathogenesis of EM. Moreover, SNARE interactions in vesicular transport and mTOR signaling pathway were the top three pathways in terms of enriching factor.

Figure 3
Figure 3

(A) GO analysis of predicted genes according to the values in the enrichment score under the themes of BP and MF. Color represents transitional values about-log10 of q-value, size represents the number of genes enrichment, shape represents GO categories. BP, biological process; MF, molecular function. The horizontal axis represents the proportion of genetic enrichment. GO enrichment analysis classifies genes according to their different functions. (B) KEGG enrichment analysis of the differentially expressed genes shown as a scatter plot. In this figure, the enrichment degree of KEGG was measured by the number of genes enriched to this pathway, Rich factor and q-value. The horizontal coordinate is a Rich factor. Greater the value of the rich factor, greater is the degree of enrichment. Rich factor represents a number of all annotated genes under this pathway term. The vertical coordinate is the pathway term with high enrichment. Q-value is the value of p after multiple verifications, and the value ranges from 0 to 1 is expressed in colors. The size of the dot indicates the number of different genes under this term – larger the dot, higher the number of genes.

Citation: Reproduction 157, 6; 10.1530/REP-18-0631

3500 miRNA downstream target genes were predicted by IPA. DAVID was used to explore the KEGG pathway enrichment terms. The top 19 of pathways are shown in Fig. 4A. We focused our attention on to mTOR signaling mechanisms because both TargetScan and IPA-predicted genes enriched in this pathway. Remarkably, four genes at different levels of the mTOR signaling pathway were targeted by miR-450a-5p and miR-503-5p (Fig. 4B and C). Based on these results, we propose that the downregulated miRNAs may elicit the constitutive activation of the mTOR signaling pathway.

Figure 4
Figure 4

(A) Top 19 of the annotated pathways for the IPA-predicted genes. (B) miRNAs downregulated in tubal EM target multiple genes in the mTOR pathway. Shown here are the two miRNAs significantly downregulated in tubal EM and their predicted target genes that participate in the mTOR signaling pathway. Circle represents the miRNAs. Figure developed from the output of IPA. (C) Details of interaction between miRNAs and molecules involved in the mTOR signaling pathway.

Citation: Reproduction 157, 6; 10.1530/REP-18-0631

Discussion

We have seen numerous advances in the field of EM in the past few decades; however, only a few studies have focused on the pathogenesis, diagnosis and treatment of tubal EM. Till now, the gold standard of diagnosis of EM was laparoscopic visualization, ideally with histological verification (Zondervan et al. 2018). In EM, miRNAs play an emerging important role and many researches have explored the clinical applications of miRNA-based therapeutics (Panir et al. 2018). MicroRNA has also been studied as a potential diagnostic marker of EM (Nisenblat et al. 2016) since the presence of altered miRNA profiles in the plasma was reported in several types of solid tumors, such as ovarian cancer (Taylor & Gercel-Taylor 2008), and other diseases. Additionally, since a single miRNA targets multiple downstream genes, miRNA dysregulation can affect several pathways involved in the development of diseases. Unfortunately, the expression profile and the role of miRNAs in tubal EM is still unknown. We performed miRNA–microarray screening followed by qRT-PCR validation, and functional analysis to explore their roles in tubal EM. Bioinformatic analysis showed that specific miRNAs might promote the occurrence and development of tubal EM by taking part in SNARE interactions in vesicular transport, the mTOR signaling pathway and endocytosis. Tubal EM may lead to infertility, but the pathogenic mechanisms of it is still unclear. By regulating target genes and participating in mTOR signaling pathway, these differentially expressed miRNAs play potential role in the pathogenesis of it. Altered miRNAs and gene expression profiling affected micro-environment in fallopian tube and then may cause infertility, though the tubal epithelium showed no significant damage. Our study emphasizes miRNAs and related signaling pathways in the tubal epithelium of tubal EM which provides new clue to the pathogenesis of tubal EM.

In this study, we identified 17 differentially expressed miRNAs in the tubal epithelium from tubal EM. What we found was expected because miRNAs play a vital role in the regulation of diverse molecular pathways as a generic mechanism and as a core disease phenotype (Ben-Hamo & Efroni 2015), which are strongly linked with the pathogenesis of EM (Teague et al. 2010, Hawkins et al. 2011, Braza-Boils et al. 2014). Previous studies found that miRNA-21 was involved in tissue repair of EM and miRNA-223 was associated with angiogenesis of EM (Teague et al. 2010). miRNA-21 was overexpressed and was found to directly participate in the regulation and activation of a number of EM-related signaling pathways such as Wnt and TGF-β (Wang et al. 2012). miRNA-21 overexpression has also been found to contribute to cell proliferation by targeting PTEN in endometrioid endometrial cancer (Qin et al. 2012). It was found that the miR-503 expression in endometriotic cyst stromal cells was downregulated and that miR-503 inhibited cell proliferation, VEGF-A production and induced apoptosis and G0/G1 cell-cycle arrest of these cells (Hirakawa et al. 2016). A reduced expression of miR-503 has also been identified in endometrial cancer (Xu et al. 2013). These findings were consistent with our results that miRNA-21-3p and miRNA-223-3p were upregulated, whereas miRNA-503-3p was downregulated in tubal EM compared to the controls. These suggest that these differentially expressed miRNAs may play a potential role in the pathogenesis and development of tubal EM.

To explore how these miRNAs participate in tubal EM pathogenesis, an miRNA-target gene-disease network was generated, which identified the potential connections between miRNAs and their downstream genes. IPA identified 78 target genes (e.g. ALDOB and ZNF695) which were associated with endometrioid carcinoma, 79 (e.g. ARSK and ELL2) with female genital tract cancer, 79 (e.g. ARSK, CDK5RAP1 and HS6ST2) with tumorigenesis of reproductive tract, 86 (e.g. ARSK, ELL2 and TOX) with female genital tract adenocarcinoma, and 87 (e.g. ARSK, CERKL and ZNF484) with endometrial carcinoma. Most target genes were cancer-related. Previous studies showed that HS6ST2 was involved in ovarian cancer (Cole et al. 2014) and pancreatic cancer (Song et al. 2011). CDK5RAP1 was associated with tumor growth in breast cancer (Wang et al. 2015). Though tubal EM is not a malignant disease, some biological behaviors of it are similar to cancer, including invasion, adhesion, angiogenesis and proliferation. Therefore, we conclude that these target genes are associated with the clinical phenotype of tubal EM and the 17 differentially expressed miRNAs may regulate the pathogenesis of tubal EM by these target genes.

In this study, bioinformatic analysis and prediction of the downstream genes of miRNAs were conducted on miRNAs isolated from the fallopian tube epithelium from tubal EM patients and control patients. The five identified miRNAs were predicted to participate in a mTOR signaling pathway which is related to autophagy and cancer. Some studies have reported that in women with EM, this pathway is overexpressed and may significantly modulate survival, the proliferation of endometriotic cells and angiogenesis in endometriotic implants (Lee & Kim 2014). The mTOR pathway is associated with cell viability and proliferation and many kinases within this pathway may be overactive in endometriotic cells (McKinnon et al. 2016). Besides, it has been suggested that mTOR inhibitors can be used in the treatment of EM. Our results suggest that the downregulated miRNAs may participate in the pathogenesis of tubal EM by activating mTOR signaling pathway.

However, certain limitations in our study should be considered when interpreting the results – first, the number of samples was small, so, recruiting more participants is important for future studies; second, target genes related to biological functions, MFs and signaling pathways were obtained from existing databases. These findings must be further confirmed by mechanistic studies performed in tubal EM cell lines and animal experiments; third, the relationship between differentially expressed miRNAs and tubal reproductive functions needs further confirmation; fourth, patients with uterine leiomyoma were recruited as control group, which affects the accuracy of our results; therefore, we will select appropriate controls in further studies.

Overall, based on the microarray data, we found for the first time, an miRNA expression profiling and an associated miRNA-gene-disease network in tubal EM. We revealed that 17 miRNAs and 4343 genes were differentially expressed in the tubal EM patients compared to the control group. Most of the genes targeted by the five identified miRNAs appear to be cancer-related and anticipate primarily in the mTOR signaling pathway. Further studies on the pathogenesis of tubal EM are necessary to confirm the five identified miRNAs as pathogenic factors of tubal EM. Our results suggest that it is feasible to recruit more patients with tubal EM to study miRNAs and their target genes within tissues in order to provide novel idea to the pathogenesis of tubal EM-related infertility and to identify molecular targets for the diagnosis and treatment of tubal EM.

Supplementary data

This is linked to the online version of the paper at https://doi.org/10.1530/REP-18-0631.

Declaration of interest

The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.

Funding

This work was supported by Shanghai Jiaotong University (grant number YG2016ZD07).

Acknowledgements

The authors thank Elsevier Webshop (http://webshop.elsevier.com) for its linguistic assistance during the preparation of this manuscript.

References

  • Ben-Hamo R & Efroni S 2015 MicroRNA regulation of molecular pathways as a generic mechanism and as a core disease phenotype. Oncotarget 6 15941604. (https://doi.org/10.18632/oncotarget.2734)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • Borghese B, Mondon F, Noel JC, Fayt I, Mignot TM, Vaiman D & Chapron C 2008 Gene expression profile for ectopic versus eutopic endometrium provides new insights into endometriosis oncogenic potential. Molecular Endocrinology 22 25572562. (https://doi.org/10.1210/me.2008-0322)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Braza-Boils A, Mari-Alexandre J, Gilabert J, Sanchez-Izquierdo D, Espana F, Estelles A & Gilabert-Estelles J 2014 MicroRNA expression profile in endometriosis: its relation to angiogenesis and fibrinolytic factors. Human Reproduction 29 978988. (https://doi.org/10.1093/humrep/deu019)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Cole CL, Rushton G, Jayson GC & Avizienyte E 2014 Ovarian cancer cell heparan sulfate 6-O-sulfotransferases regulate an angiogenic program induced by heparin-binding epidermal growth factor (EGF)-like growth factor/EGF receptor signaling. Journal of Biological Chemistry 289 1048810501. (https://doi.org/10.1074/jbc.M113.534263)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Eyster KM, Klinkova O, Kennedy V & Hansen KA 2007 Whole genome deoxyribonucleic acid microarray analysis of gene expression in ectopic versus eutopic endometrium. Fertility and Sterility 88 15051533. (https://doi.org/10.1016/j.fertnstert.2007.01.056)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Flores I, Rivera E, Ruiz LA, Santiago OI, Vernon MW & Appleyard CB 2007 Molecular profiling of experimental endometriosis identified gene expression patterns in common with human disease. Fertility and Sterility 87 11801199. (https://doi.org/10.1016/j.fertnstert.2006.07.1550)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Hawkins SM, Creighton CJ, Han DY, Zariff A, Anderson ML, Gunaratne PH & Matzuk MM 2011 Functional microRNA involved in endometriosis. Molecular Endocrinology 25 821832. (https://doi.org/10.1210/me.2010-0371)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Hever A, Roth RB, Hevezi P, Marin ME, Acosta JA, Acosta H, Rojas J, Herrera R, Grigoriadis D & White E et al. 2007 Human endometriosis is associated with plasma cells and overexpression of B lymphocyte stimulator. PNAS 104 1245112456. (https://doi.org/10.1073/pnas.0703451104)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Hirakawa T, Nasu K, Abe W, Aoyagi Y, Okamoto M, Kai K, Takebayashi K & Narahara H 2016 miR-503, a microRNA epigenetically repressed in endometriosis, induces apoptosis and cell-cycle arrest and inhibits cell proliferation, angiogenesis, and contractility of human ovarian endometriotic stromal cells. Human Reproduction 31 25872597. (https://doi.org/10.1093/humrep/dew217)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Hull ML, Escareno CR, Godsland JM, Doig JR, Johnson CM, Phillips SC, Smith SK, Tavare S, Print CG & Charnock-Jones DS 2008 Endometrial-peritoneal interactions during endometriotic lesion establishment. American Journal of Pathology 173 700715. (https://doi.org/10.2353/ajpath.2008.071128)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Jia SZ, Yang Y, Lang J, Sun P & Leng J 2013 Plasma miR-17-5p, miR-20a and miR-22 are down-regulated in women with endometriosis. Human Reproduction 28 322330. (https://doi.org/10.1093/humrep/des413)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Lee II & Kim JJ 2014 Influence of AKT on progesterone action in endometrial diseases. Biology of Reproduction 91 63. (https://doi.org/10.1095/biolreprod.114.119255)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • McKinnon BD, Kocbek V, Nirgianakis K, Bersinger NA & Mueller MD 2016 Kinase signalling pathways in endometriosis: potential targets for non-hormonal therapeutics. Human Reproduction Update 22 382403. (https://doi.org/10.1093/humupd/dmv060)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Nisenblat V, Bossuyt PM, Shaikh R, Farquhar C, Jordan V, Scheffers CS, Mol BW, Johnson N & Hull ML 2016 Blood biomarkers for the non-invasive diagnosis of endometriosis. Cochrane Database of Systematic Reviews CD012179. (https://doi.org/10.1002/14651858.CD012179)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • Ohlsson Teague EM, Van Der Hoek KH, Van Der Hoek MB, Perry N, Wagaarachchi P, Robertson SA, Print CG & Hull LM 2009 MicroRNA-regulated pathways associated with endometriosis. Molecular Endocrinology 23 265275. (https://doi.org/10.1210/me.2008-0387)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Okamoto M, Nasu K, Abe W, Aoyagi Y, Kawano Y, Kai K, Moriyama M & Narahara H 2015 Enhanced miR-210 expression promotes the pathogenesis of endometriosis through activation of signal transducer and activator of transcription 3. Human Reproduction 30 632641. (https://doi.org/10.1093/humrep/deu332)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Ouellet DL, Perron MP, Gobeil LA, Plante P & Provost P 2006 MicroRNAs in gene regulation: when the smallest governs it all. Journal of Biomedicine and Biotechnology 2006 69616. (https://doi.org/10.1155/JBB/2006/69616)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • Panir K, Schjenken JE, Robertson SA & Hull ML 2018 Non-coding RNAs in endometriosis: a narrative review. Human Reproduction Update 24 497515. (https://doi.org/10.1093/humupd/dmy014)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Qin LX & Satagopan JM 2009 Normalization method for transcriptional studies of heterogeneous samples--simultaneous array normalization and identification of equivalent expression. Statistical Applications in Genetics and Molecular Biology 8 Article 10. (https://doi.org/10.2202/1544-6115.1339)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • Qin X, Yan L, Zhao X, Li C & Fu Y 2012 microRNA-21 overexpression contributes to cell proliferation by targeting PTEN in endometrioid endometrial cancer. Oncology Letters 4 12901296. (https://doi.org/10.3892/ol.2012.896)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Singh SK, Pal Bhadra M, Girschick HJ & Bhadra U 2008 MicroRNAs--micro in size but macro in function. FEBS Journal 275 49294944. (https://doi.org/10.1111/j.1742-4658.2008.06624.x)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Song L & Tuan RS 2006 MicroRNAs and cell differentiation in mammalian development. Birth Defects Research Part C: Embryo Today 78 140149. (https://doi.org/10.1002/bdrc.20070)

  • Song K, Li Q, Peng YB, Li J, Ding K, Chen LJ, Shao CH, Zhang LJ & Li P 2011 Silencing of hHS6ST2 inhibits progression of pancreatic cancer through inhibition of Notch signalling. Biochemical Journal 436 271282. (https://doi.org/10.1042/BJ20110297)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Taganov KD, Boldin MP & Baltimore D 2007 MicroRNAs and immunity: tiny players in a big field. Immunity 26 133137. (https://doi.org/10.1016/j.immuni.2007.02.005)

  • Taylor DD & Gercel-Taylor C 2008 MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecologic Oncology 110 1321. (https://doi.org/10.1016/j.ygyno.2008.04.033)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Teague EM, Print CG & Hull ML 2010 The role of microRNAs in endometriosis and associated reproductive conditions. Human Reproduction Update 16 142165. (https://doi.org/10.1093/humupd/dmp034)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Wang T, Zhang L, Shi C, Sun H, Wang J, Li R, Zou Z, Ran X & Su Y 2012 TGF-β-induced miR-21 negatively regulates the antiproliferative activity but has no effect on EMT of TGF-β in HaCaT cells. International Journal of Biochemistry and Cell Biology 44 366376. (https://doi.org/10.1016/j.biocel.2011.11.012)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Wang H, Wei L, Li C, Zhou J & Li Z 2015 CDK5RAP1 deficiency induces cell cycle arrest and apoptosis in human breast cancer cell line by the ROS/JNK signaling pathway. Oncology Reports 33 10891096. (https://doi.org/10.3892/or.2015.3736)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Woodward PJ, Sohaey R & Mezzetti TP 2001 Endometriosis: radiologic-pathologic correlation. RadioGraphics 21 193216. (https://doi.org/10.1148/radiographics.21.1.g01ja14193)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Xia W, Zhang D, Ouyang J, Liang Y, Zhang H, Huang Z, Liang G, Zhu Q, Guan X & Zhang J 2018 Effects of pelvic endometriosis and adenomyosis on ciliary beat frequency and muscular contractions in the human fallopian tube. Reproductive Biology and Endocrinology 16 48. (https://doi.org/10.1186/s12958-018-0361-y)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Xu YY, Wu HJ, Ma HD, Xu LP, Huo Y & Yin LR 2013 MicroRNA-503 suppresses proliferation and cell-cycle progression of endometrioid endometrial cancer by negatively regulating cyclin D1. FEBS Journal 280 37683779. (https://doi.org/10.1111/febs.12365)

    • Crossref
    • PubMed
    • Search Google Scholar
    • Export Citation
  • Xue H & Zhang J 2017 Research status of tubal endometriosis. Progress in Obstetrics and Gynecology 26 98100.

  • Zondervan KT, Becker CM, Koga K, Missmer SA, Taylor RN & Viganò P 2018 Endometriosis. Nature Reviews Disease Primers 4 9. (https://doi.org/10.1038/s41572-018-0008-5)

 

  • Collapse
  • Expand
  • (A) Box plots showing the distribution of miRNAs in the nine samples. The majority of the miRNAs on the chip were not expressed in the sample or the expression level was very low. The signal is focused on a relatively flat casing, and a small number of miRNAs with a relatively high signal are located in the upper third quartile. (B) Hierarchical clustering of miRNAs. The samples are in columns and miRNAs in rows. BL represents case group and DZ represents the control group. The red columns represent upregulation in the case group, and the blue columns represent downregulation in the case group. The brightness of the color indicates the degree of upregulation or downregulation. (C, D, E, F and G) Verification of five differentially expressed miRNAs related to tubal EM identified by miRNA-microarray. The expression levels of the five differentially expressed miRNAs in tubal epithelium of tubal EM (n = 6) and control women (n = 6). (C, D and E) are upregulated, and F-G are downregulated in transcription analysis and qRT-PCR. Data are presented as mean ± s.e. P value was determined by unpaired Student’s t-test. *P < 0.05, **0.01 < P < 0.05 and ***P < 0.01, compared to the corresponding controls.

  • Inverse relationship of differentially expressed miRNAs with target genes and gynecological diseases. The central blue comb-like sign represents miRNAs, the red rectangular, triangular, trapezoidal and oval boxes, represent their target genes; the blue cross-shaped boxes represent gynecological diseases and the discontinuous lines connecting them represent the interaction. (A) Network of miR-223-3P-downstream gene-gynecological diseases. (B) Network of miR-450a-5p, miR-503-5p-downstream gene-gynecological diseases.

  • (A) GO analysis of predicted genes according to the values in the enrichment score under the themes of BP and MF. Color represents transitional values about-log10 of q-value, size represents the number of genes enrichment, shape represents GO categories. BP, biological process; MF, molecular function. The horizontal axis represents the proportion of genetic enrichment. GO enrichment analysis classifies genes according to their different functions. (B) KEGG enrichment analysis of the differentially expressed genes shown as a scatter plot. In this figure, the enrichment degree of KEGG was measured by the number of genes enriched to this pathway, Rich factor and q-value. The horizontal coordinate is a Rich factor. Greater the value of the rich factor, greater is the degree of enrichment. Rich factor represents a number of all annotated genes under this pathway term. The vertical coordinate is the pathway term with high enrichment. Q-value is the value of p after multiple verifications, and the value ranges from 0 to 1 is expressed in colors. The size of the dot indicates the number of different genes under this term – larger the dot, higher the number of genes.

  • (A) Top 19 of the annotated pathways for the IPA-predicted genes. (B) miRNAs downregulated in tubal EM target multiple genes in the mTOR pathway. Shown here are the two miRNAs significantly downregulated in tubal EM and their predicted target genes that participate in the mTOR signaling pathway. Circle represents the miRNAs. Figure developed from the output of IPA. (C) Details of interaction between miRNAs and molecules involved in the mTOR signaling pathway.