Skip to main content
Download PDF

A cross-study gene set enrichment analysis identifies critical pathways in endometriosis

Reproductive Biology and Endocrinology volume 7, Article number: 94 (2009) Cite this article

Abstract

Background

Endometriosis is an enigmatic disease. Gene expression profiling of endometriosis has been used in several studies, but few studies went further to classify subtypes of endometriosis based on expression patterns and to identify possible pathways involved in endometriosis. Some of the observed pathways are more inconsistent between the studies, and these candidate pathways presumably only represent a fraction of the pathways involved in endometriosis.

Methods

We applied a standardised microarray preprocessing and gene set enrichment analysis to six independent studies, and demonstrated increased concordance between these gene datasets.

Results

We find 16 up-regulated and 19 down-regulated pathways common in ovarian endometriosis data sets, 22 up-regulated and one down-regulated pathway common in peritoneal endometriosis data sets. Among them, 12 up-regulated and 1 down-regulated were found consistent between ovarian and peritoneal endometriosis. The main canonical pathways identified are related to immunological and inflammatory disease. Early secretory phase has the most over-represented pathways in the three uterine cycle phases. There are no overlapping significant pathways between the dataset from human endometrial endothelial cells and the datasets from ovarian endometriosis which used whole tissues.

Conclusion

The study of complex diseases through pathway analysis is able to highlight genes weakly connected to the phenotype which may be difficult to detect by using classical univariate statistics. By standardised microarray preprocessing and GSEA, we have increased the concordance in identifying many biological mechanisms involved in endometriosis. The identified gene pathways will shed light on the understanding of endometriosis and promote the development of novel therapies.

Background

Endometriosis is defined as the presence of endometrium-like tissue in sites outside the uterine cavity and occurs in 6-10% of women in the general population [1]. The main clinical features are chronic pelvic pain, pain during intercourse, and infertility [2]. As cellular and molecular mechanisms involved in endometriosis are still uncovered, the classification of this disease evolved from a local disorder to a complex, chronic systemic disease [3]. Despite extensive researches, the etiology of endometriosis remains obscure. Gene expression profiling has been used in several studies of endometriosis, in which from a few to hundreds differentially expressed genes were identified [4–17]. For previously identified genes, their roles in the pathogenesis of endometriosis are further discussed. But it is hard to interpret individual genes on a list with many significant genes.

A common challenge in the analysis of genome wide expression no longer lies in obtaining gene expression profiles, but rather in interpreting the results to gain insights into biological mechanisms [18]. Pathway analysis of microarray data evaluates gene expression profiles of a priori defined biological pathways in association with a phenotype of interest. Recently gene expression patterns were further used in the classification of subtypes of endometriosis as well as in the identification of the pathways involved in endometriosis [4, 13–16]. So far the observed pathways were discordant between the studies that suggest that these previously identified pathways only represent a fraction of the pathways involved in endometriosis.

Currently the most well-known and widely used approach to gene set analysis, the Gene Set Enrichment Analysis (GSEA) method was introduced by Mootha et al. [19], which was used to identify pre-defined gene sets which exhibited significant differences in expression between samples from normal and patients. The methodology was subsequently refined by Subramanian et al. [18]. The algorithms calculate the statistical significance of the expression changes across groups or pathways rather than individual gene, thus allowing identification of groups or pathways most strongly affected by the observed expression changes. The analysis based on a group of relevant genes instead of on an individual gene increases the likelihood for investigators to identify the critical functional processes under the biological phenomena under study. GSEA is likely to be more powerful than conventional single-gene methods in the study of complex diseases in which many genes make subtle contributions [20].

In a single data set the GSEA will generally not result in significant findings beyond major pathways. Here we will use standardised microarray preprocessing and GSEA with comprehensive expression profiles in an attempt to find greater data convergence and provide a systematic insight into the pathways altered during endometriosis pathogenesis.

Methods

Datasets

We searched GEO http://www.ncbi.nlm.nih.gov/geo/, and ArrayExpress http://www.ebi.ac.uk/arrayexpress/ for the gene expression profiling studies related to endometriosis disease. Data were included in our re-analysis if they met the following conditions: 1) the data is in genome-wide; 2) comparison was conducted between endometriosis patients and controls; 3) complete microarray raw or normalized data are available.

Finally six public gene expression data sets were involved in our study, which assessed endometriosis transcripts on a genome-wide basis. In data set GSE7307, total 677 samples from more than 90 distinct tissue types were processed, but only the profiles related to endometriosis and eutopic endometrium were considered here. The data generated from human endometrial endothelial cells by Sha et al. [4] were also included in our combined re-analysis to compare with the whole endometriosis tissues data sets. The related information about these datasets, such as the microarray platform, sample type, sample size, is listed in Table 1.

Table 1 Characteristics of datasets included in the studies.
Full size table

Data Preprocessing

Data preprocessing was performed using software packages developed in version 2.4.0 of Bioconductor [21] and R version 2.9.0 [22]. Each Affymetrix dataset was background adjusted, normalized and log2 probe-set intensities calculated using the Robust Multichip Averaging (RMA) algorithm in affy package [23, 24], and the Codelink arrays normalizations performed in GSE5108 were retained. Genes which cannot be mapped to any KEGG pathway identified were excluded from the further analysis. The interquartile range (IQR) was used as a measure of variability. From the resulting distribution of IQR values for all genes we set a cut-off so as to exclude values under 0.5. When multiple probe sets target one gene, the probe set with largest variability was kept. Pathway analysis was performed separately in each data set.

Gene set enrichment analysis of pathways

GSEA implemented in the Category package (version 2.10.1 [25]). The goal of GSEA is to determine whether the members of a gene set S randomly distributed throughout the entire reference gene list L or are primarily found at the top or bottom. One of the advantages of GSEA is the relative robustness to noise and outliers in the data. The gene sets represented by less than 10 genes were excluded. The t-statistic mean of the genes was computed in each pathway. Using a permutation test with 1000 times, the significantly changed pathways were identified with p-value ≤ 0.05.

Results

For the studies which used multiple locations (ovarian and peritoneal) or uterine cycle phases (proliferative to secretory) of endometriosis, each type or phase was treated as a separate data set. These six studies provided 9 case-control data sets including 74 endometriosis cases and 74 controls. Common GSEA method was applied to the 9 datasets. For individual analysis, we obtained the significant pathways and the genes included [See Additional file 1, 2, 3 and 4]. The analysis results were summarised in table 2. We postulated that the pathways and genes that appear consistently as different expressed in multiple studies are more likely to be important in endometriosis. To look for such convergence we compared the GSEA results.

Table 2 Summary of each data set used in the re-analysis and the number of differentially expressed pathways.
Full size table

Common significant pathways in ovarian endometriosis

Endometriosis is most commonly localized in the ovaries. We first compared the ovarian endometriosis data sets for their convergence and reproducibility. The results for the ovarian endometriosis data sets were presented in Table 3. Based on the permutation p-values, the GSEA reported 24 & 31 (Eyster data set), 56 & 59 (Hever data set), and 36 & 44 (Roth data set) up- and down- regulated pathways. The overlaps among them are presented in Figure 1, with 16 up-regulated and 19 down-regulated pathway detected in common. Among the up-regulated list, the main canonical pathways affected are related to immunological and inflammatory diseases including Asthma, Autoimmune thyroid disease, Systemic lupus erythematosus, Allograft rejection, Graft-versus-host disease and Type I diabetes mellitus. Cytokine-cytokine receptor interaction and Cell adhesion molecules (CAMs) belongs to signaling molecules and interaction. PPAR signaling pathway is part of endocrine system. The most down-regulated pathways are metabolism and repair related pathways. And we extracted the collective genes in each pathway from the 3 data sets [See Additional file 5].

Figure 1
figure 1

Significant pathways identified and overlap between the 3 ovarian endometriosis data sets. For each data set, GSEA implanted to generate p-values for each pathway and using a permutation test with 1000 times, obtained the significant pathways with p-value cut-off of ≤ 0.05. (a) GSEA detected 24 (Eyster data set), 56 (Hever data set) and 36 (Roth data set) up-regulated pathways and 16 are common. (b) GSEA detected 31 (Eyster data set), 59 (Hever data set) and 44(Roth data set) down-regulated pathways and 19 are common.

Full size image
Table 3 Common pathway categories identified by GSEA of the 3 ovarian endometriosis datasets
Full size table

Common significant pathways in peritoneal endometriosis

We also analyzed the data sets from peritoneal endometriosis which is another important endometriosis type. 23 pathways were commonly significantly regulated in the two peritoneal endometriosis data sets, summarised in table 4 and Figure 2. Of these 23 pathways, 12 up- and one down- regulation pathways were also identified from the three ovarian endometriosis data sets. The significances of four pathways (Drug metabolism - cytochrome P450, Metabolism of xenobiotics by cytochrome P450, Olfactory transduction, Toll-like receptor signaling pathway) were found from the data of Eyster et al. [15], but not in the data of Hull et al. [17]. Collective genes in each pathway were extracted [See Additional file 6].

Figure 2
figure 2

Significant pathways identified and overlap between the 2 peritoneal endometriosis data sets. For each data set, GSEA implanted to generate p-values for each pathway and using a permutation test with 1000 times, obtained the significant pathways with p-value cut-off of ≤ 0.05. (a) GSEA detected 57 (Eyster data set) and 29 (Hull data set) up-regulated pathways and 22 are common. (b) GSEA detected 11 (Eyster data set) and 7 (Hull data set) down-regulated pathways and 1 is common.

Full size image
Table 4 Common pathway categories identified by GSEA of the 2 peritoneal endometriosis datasets
Full size table

Differential pathways between the timing of the uterine cycle

Burney et al. conducted global gene expression analysis of endometrium from women with and without moderate/severe stage endometriosis and compared the gene expression signatures across various phases of the menstrual cycle [14]. Specimens were classified as proliferative (PE, d 8-14), early secretory (ESE, d 15-18), midsecretory (MSE, d 19-23). They did not examine later than Day 23 in the menstrual cycle. Each phase had endometriosis more than one type. We re-analyzed their microarray CEL files and each phase was treated as a separate data set. In proliferative phase, 4 pathways were up-regulated and 3 pathways were down-regulated. Early secretory phase had the most over-represented pathways, 12 were up-regulated and 21 were down-regulated. There is no significant pathway up-regulated at p ≤ 0.05 in midsecretory phase while 22 down-regulated pathways were identified [See Additional file 3]. The overlap of pathways among these phases is very low (Figure 3).

Figure 3
figure 3

Significant pathways identified and overlap between the phases of uterine cycle. For each phase, GSEA implanted to generate p-values for each pathway and using a permutation test with 1000 times, obtained the significant pathways with p-value cut-off of ≤ 0.05. (a) GSEA detected 4 (Proliferative phase), 12 (Early secretory phase) and 0 (Mid secretory phase) up-regulated pathways and no overlap. (a) GSEA detected 3 (Proliferative phase), 21 (Early secretory phase) and 22 (Mid secretory phase) down-regulated pathways and the overlap is very low.

Full size image

Differential expression in human endometrial endothelial cells

Sha et al. selected the eutopic endothelial compartment as the subject for exploring the differential expression profile between endometriosis patients and normal controls [4]. Our re-analysis revealed 44 up-regulated and 1 down-regulated pathways in the dataset [See Additional file 4]. There is no overlap compared with the common list from ovarian endometriosis data sets using the whole tissue as sample.

Discussion

Endometriosis is an enigmatic disease. No existing single theory can explain all cases of endometriosis. The genome-wide microarrays are very powerful because they allow the identification of gene families or pathways that change in concert in a disease state comprehensively. Biologically relevant inference should therefore be reproducible across laboratories. For single gene analysis, different statistical methods and different datasets examining the same biological condition may lead to significant discrepancies [26]. Pathways analysis applied to different datasets yields interesting common results, diminishing the large discrepancies observed in direct comparisons of lists of differentially expressed genes obtained from different datasets. Therefore, the study of complex diseases through pathway analysis is able to highlight genes weakly connected to the phenotype which may be difficult to detect by using classical univariate statistics.

We have performed gene set enrichment analysis of six independent publicly available gene expression data sets to understand in depth the common biological mechanisms involved in endometriosis. Our study compared the gene expression between lesion locations (ovarian vs. peritoneal), phases of the uterine cycle (proliferative to midsecretory) and cell types (endometrial endothelial cells vs. whole tissue), as well as overall eutopic versus ectopic endometrium. The transcriptomes of eutopic endometrium and ectopic endometrial lesions suggest that ovarian endometriosis and peritoneal disease are different disorders [13]. Our findings suggest that most of the pathways impacted the ovarian and peritoneal endometriosis are consistent. Many of differentially expressed pathways found in this study have already been reported to be involved in endometriosis pathogenesis. Here, this discussion presents several of the differentially expressed pathways and hypotheses regarding the role of these pathways in endometriosis.

Most significant of the common up-regulated pathways are involved in immune system and immune disorders. It has been widely documented that endometriosis, as an inflammatory disease, induces an immune response, leading to both cellular and humoral immune changes [27, 28]. The association between endometriosis and immune disorders were literature supported [29]. Also some studies concluded that women with endometriosis do not have a higher risk of having asthma, systemic lupus erythematosus and Sjögren's syndrome than other subjects [30, 31]. Our GSEA results showed that expression of Asthma, Graft-versus-host disease, Autoimmune thyroid disease, Allograft rejection, Systemic lupus erythematosus and Type I diabetes mellitus pathways are the significantly imbalanced between endometriosis and eutopic endometrium. We found that human leukocyte antigen (HLA) genes are critical genes in these pathways. HLA are key components of the major histocompatibility complex (MHC), which is involved in immune cell signalling processes such as T-cell activation. People with certain HLA antigens are more likely to develop certain autoimmune diseases, such as Type I Diabetes, Ankylosing spondylitis, Celiac Disease, Systemic Lupus Erythematosus, Myasthenia Gravis and Sjögren's syndrome et al [32].

Cytokine-cytokine receptor interaction and Cell adhesion molecules (CAMs) included in GSEA were up-regulated in endometriosis. Cell adhesion molecules are (glyco) proteins expressed on the cell surface and play a critical role in a wide array of biologic processes that include hemostasis, the immune response, inflammation, embryogenesis, and development of neuronal tissue. Clinical observations and in vitro experiments imply that endometriotic cells are invasive and able to metastasize. Analogous to tumour metastasis, it is likely that cell adhesion molecules are central for the invasion and metastasis of endometriotic cells. The expression of some integrins is aberrant in endometriotic lesions compared to eutopic endometrium [33]. Cytokines are key mediators of intercellular communication within the immune system. Several cytokines including interleukin (IL)-1, 6, 8, 10, tumor necrosis factor (TNF)-α, and vascular endothelial growth factor (VEGF) were reported to be increased in the peritoneal fluid (PF) of women with endometriosis [34–41]. Peroxisome proliferator-activated receptors (PPARs) signaling pathway is up-regulated according to GSEA. PPAR are nuclear hormone receptors that are activated by fatty acids and their derivatives. PPAR-γ is present in human ovarian cells. Activation of PPAR-γ enhances steroidogenesis via activation of StAR protein and leads to the activation of insulin-signaling pathways [42].

The expression patterns of ER (estrogen receptors) and PR (progesterone receptors) in endometriotic lesions are different from those in the eutopic endometrium. Endometriosis is an estrogen-dependent disease [43]. Studies of hormone-ligand binding assays and enzyme immunoassays showed a consistent reduction in the content of ER and PR in endometriotic implants [44–47]. Androgen and estrogen metabolism pathway appeared in most of our re-analysis results. Oxidative phosphorylation pathway possibly affect oocyte quality, fertilization rate, and further embryo development [48], is down-regulated in our analysis.

Burney et al. studied proliferative, early secretory and mid-secretory eutopic endometrium (up to Day 23), from women with endometriosis and controls. They found that endometrial gene expression differed most, between these groups, in the early secretory phase (Days 15-18). They found far fewer differences in the mid-secretory phase where no transcripts were found to be up- or down-regulated 4-fold. The result is consistent with the findings of other studies [12, 49, 50]. Corresponding to their result, there is no significant up-regulated pathway in mid-secretory phase by our GSEA. The molecular phenotype of mid-secretory, eutopic endometrium from women with endometriosis and from controls appears to be very similar [50]. Current efforts to develop minimally invasive diagnostic tests for the presence of endometriosis and also tests to distinguish minimal/mild and moderate/severe disease, by sampling the endometrium, should be focused on the early secretory phase of the menstrual cycle [50].

We hypothesize that all cell types in the endometriotic lesion contribute to the pathology of the disease. Matsuzaki et al. compared global gene expression in eutopic endometrium, from controls and patients with deep endometriosis, at various time points throughout the menstrual cycle. They found no genes were up- or down-regulated in all phases of the cycle, in either tissue compartments [12]. None of the genes from their study that had been identified as differentially expressed in either the stromal or epithelial compartments was shown to be differentially expressed. This may be due to the relative contribution that the epithelial and stromal transcriptomes make to whole tissue gene expression. Our GSEA result showed that the significant pathways from the human endometrial endothelial cells had low overlap with the list from ovarian endometriosis data sets used whole tissue.

Conclusion

The pathogenesis of endometriosis is likely multifactorial. A deeper understanding of the mechanisms of these diseases can be reached by focusing on deregulation of gene sets or pathways rather than on individual genes. By standardised microarray preprocessing and GSEA, we have increased the concordance to identify many biological mechanisms are involved in endometriosis which are novel in terms of their connection to endometriosis (as mined from the existing literature). More studies about the specific role and interactions of the genes included in related pathways are needed to improve the understanding of endometriosis.

References

  1. Eskenazi B, Warner ML: Epidemiology of endometriosis. Obstet Gynecol Clin North Am. 1997, 24 (2): 235-258. 10.1016/S0889-8545(05)70302-8.

    Article  CAS  PubMed  Google Scholar 

  2. Giudice LC, Kao LC: Endometriosis. Lancet. 2004, 364 (9447): 1789-1799. 10.1016/S0140-6736(04)17403-5.

    Article  PubMed  Google Scholar 

  3. Bulun SE: Endometriosis. N Engl J Med. 2009, 360 (3): 268-279. 10.1056/NEJMra0804690.

    Article  CAS  PubMed  Google Scholar 

  4. Sha G, Wu D, Zhang L, Chen X, Lei M, Sun H, Lin S, Lang J: Differentially expressed genes in human endometrial endothelial cells derived from eutopic endometrium of patients with endometriosis compared with those from patients without endometriosis. Hum Reprod. 2007, 22 (12): 3159-3169. 10.1093/humrep/dem266.

    Article  CAS  PubMed  Google Scholar 

  5. Kato N, Sasou S, Motoyama T: Expression of hepatocyte nuclear factor-1beta (HNF-1beta) in clear cell tumors and endometriosis of the ovary. Mod Pathol. 2006, 19 (1): 83-89. 10.1038/modpathol.3800492.

    Article  CAS  PubMed  Google Scholar 

  6. Eyster KM, Boles AL, Brannian JD, Hansen KA: DNA microarray analysis of gene expression markers of endometriosis. Fertil Steril. 2002, 77 (1): 38-42. 10.1016/S0015-0282(01)02955-7.

    Article  PubMed  Google Scholar 

  7. Lebovic DI, Baldocchi RA, Mueller MD, Taylor RN: Altered expression of a cell-cycle suppressor gene, Tob-1, in endometriotic cells by cDNA array analyses. Fertil Steril. 2002, 78 (4): 849-854. 10.1016/S0015-0282(02)03319-8.

    Article  PubMed  Google Scholar 

  8. Arimoto T, Katagiri T, Oda K, Tsunoda T, Yasugi T, Osuga Y, Yoshikawa H, Nishii O, Yano T, Taketani Y, Nakamura Y: Genome-wide cDNA microarray analysis of gene-expression profiles involved in ovarian endometriosis. Int J Oncol. 2003, 22 (3): 551-560.

    CAS  PubMed  Google Scholar 

  9. Kao LC, Germeyer A, Tulac S, Lobo S, Yang JP, Taylor RN, Osteen K, Lessey BA, Giudice LC: Expression profiling of endometrium from women with endometriosis reveals candidate genes for disease-based implantation failure and infertility. Endocrinology. 2003, 144 (7): 2870-2881. 10.1210/en.2003-0043.

    Article  CAS  PubMed  Google Scholar 

  10. Fan Y, Chen BL, Ma XD, Su MQ: Detection of expression of endometriosis-related cytokine and their receptor genes by cDNA microarray technique. Xi Bao Yu Fen Zi Mian Yi Xue Za Zhi. 2005, 21 (4): 489-492.

    CAS  PubMed  Google Scholar 

  11. Matsuzaki S, Canis M, Pouly JL, Dechelotte P, Okamura K, Mage G: The macrophage stimulating protein/RON system: a potential novel target for prevention and treatment of endometriosis. Mol Hum Reprod. 2005, 11 (5): 345-349. 10.1093/molehr/gah162.

    Article  CAS  PubMed  Google Scholar 

  12. Matsuzaki S, Canis M, Vaurs-Barriere C, Boespflug-Tanguy O, Dastugue B, Mage G: DNA microarray analysis of gene expression in eutopic endometrium from patients with deep endometriosis using laser capture microdissection. Fertil Steril. 2005, 84 (Suppl 2): 1180-1190. 10.1016/j.fertnstert.2005.04.041.

    Article  CAS  PubMed  Google Scholar 

  13. Wu Y, Kajdacsy-Balla A, Strawn E, Basir Z, Halverson G, Jailwala P, Wang Y, Wang X, Ghosh S, Guo SW: Transcriptional characterizations of differences between eutopic and ectopic endometrium. Endocrinology. 2006, 147 (1): 232-246. 10.1210/en.2005-0426.

    Article  CAS  PubMed  Google Scholar 

  14. Burney R, Talbi S, Hamilton A, Vo K, Nyegaard M, Nezhat C, Lessey B, Giudice L: Gene Expression Analysis of Endometrium Reveals Progesterone Resistance and Candidate Susceptibility Genes in Women with Endometriosis. Endocrinology. 2007, 148 (8): 3814-10.1210/en.2006-1692.

    Article  CAS  PubMed  Google Scholar 

  15. Eyster KM, Klinkova O, Kennedy V, Hansen KA: Whole genome deoxyribonucleic acid microarray analysis of gene expression in ectopic versus eutopic endometrium. Fertil Steril. 2007, 88 (6): 1505-1533. 10.1016/j.fertnstert.2007.01.056.

    Article  CAS  PubMed  Google Scholar 

  16. Hever A, Roth R, Hevezi P, Marin M, Acosta J, Acosta H, Rojas J, Herrera R, Grigoriadis D, White E: Human endometriosis is associated with plasma cells and overexpression of B lymphocyte stimulator. Proceedings of the National Academy of Sciences. 2007, 104 (30): 12451-10.1073/pnas.0703451104.

    Article  CAS  Google Scholar 

  17. Hull ML, Escareno CR, Godsland JM, Doig JR, Johnson CM, Phillips SC, Smith SK, Tavare S, Print CG, Charnock-Jones DS: Endometrial-peritoneal interactions during endometriotic lesion establishment. Am J Pathol. 2008, 173 (3): 700-715. 10.2353/ajpath.2008.071128.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, Houstis N, Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B, Lander ES, Hirschhorn JN, Altshuler D, Groop LC: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.

    Article  CAS  PubMed  Google Scholar 

  20. Jing Shi, Walker MG: Gene set enrichment analysis (GSEA) for interpreting gene expression profiles. Current Bioinformatics. 2007, 2 (2): 133-137. 10.2174/157489307780618231.

    Article  Google Scholar 

  21. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.

    Article  PubMed Central  PubMed  Google Scholar 

  22. R Development Core Team: R: A language and environment for statistical computing. 2009, R Foundation for Statistical Computing, Vienna, Austria, ISBN 3-900051-07-0, [http://www.R-project.org/]

    Google Scholar 

  23. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.

    Article  PubMed  Google Scholar 

  24. Gautier L, Cope L, Bolstad BM, Irizarry RA: affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20 (3): 307-315. 10.1093/bioinformatics/btg405.

    Article  CAS  PubMed  Google Scholar 

  25. Using Categories to Model Genomic Data. 2009, [http://www.bioconductor.org/packages/2.4/bioc/vignettes/Category/inst/doc/Category.pdf]

  26. Manoli T, Gretz N, Grone H-J, Kenzelmann M, Eils R, Brors B: Group testing for pathway analysis improves comparability of different microarray datasets. Bioinformatics. 2006, 22 (20): 2500-2506. 10.1093/bioinformatics/btl424.

    Article  CAS  PubMed  Google Scholar 

  27. Sinaii N, Cleary SD, Ballweg ML, Nieman LK, Stratton P: High rates of autoimmune and endocrine disorders, fibromyalgia, chronic fatigue syndrome and atopic diseases among women with endometriosis: a survey analysis. Hum Reprod. 2002, 17 (10): 2715-2724. 10.1093/humrep/17.10.2715.

    Article  CAS  PubMed  Google Scholar 

  28. Lebovic DI, Mueller MD, Taylor RN: Immunobiology of endometriosis. Fertility and Sterility. 2001, 75 (1): 1-10. 10.1016/S0015-0282(00)01630-7.

    Article  CAS  PubMed  Google Scholar 

  29. Berbic M, Schulke L, Markham R, Tokushige N, Russell P, Fraser IS: Macrophage expression in endometrium of women with and without endometriosis. Hum Reprod. 2009, 24 (2): 325-332. 10.1093/humrep/den393.

    Article  CAS  PubMed  Google Scholar 

  30. Ferrero S, Petrera P, Colombo BM, Navaratnarajah R, Parisi M, Anserini P, Remorgida V, Ragni N: Asthma in women with endometriosis. Hum Reprod. 2005, 20 (12): 3514-3517. 10.1093/humrep/dei263.

    Article  CAS  PubMed  Google Scholar 

  31. Matorras R, Ocerin I, Unamuno M, Nieto A, Peiro E, Burgos J, Exposito A: Prevalence of endometriosis in women with systemic lupus erythematosus and Sjogren's syndrome. Lupus. 2007, 16 (9): 736-740. 10.1177/0961203307081339.

    Article  CAS  PubMed  Google Scholar 

  32. de Bakker PIW, McVean G, Sabeti PC, Miretti MM, Green T, Marchini J, Ke XY, Monsuur AJ, Whittaker P, Delgado M, Morrison J, Richardson A, Walsh EC, Gao XJ, Galver L, Hart J, Hafler DA, Pericak-Vance M, Todd JA, Daly MJ, Trowsdale J, Wijmenga C, Vyse TJ, Beck S, Murray SS, Carrington M, Gregory S, Deloukas P, Rioux JD: A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC. Nature Genet. 2006, 38 (10): 1166-1172. 10.1038/ng1885.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Starzinski-Powitz A, Handrow-Metzmacher H, Kotzian S: The putative role of cell adhesion molecules in endometriosis: can we learn from tumour metastasis?. Molecular Medicine Today. 1999, 5 (7): 304-309. 10.1016/S1357-4310(99)01497-5.

    Article  CAS  PubMed  Google Scholar 

  34. Wu MY, Ho HN, Chen SU, Chao KH, Chen CD, Yang YS: Increase in the production of interleukin-6, interleukin-10, and interleukin-12 by lipopolysaccharide-stimulated peritoneal macrophages from women with endometriosis. Am J Reprod Immunol. 1999, 41 (1): 106-111.

    Article  CAS  PubMed  Google Scholar 

  35. Ho HN, Wu MY, Chao KH, Chen CD, Chen SU, Yang YS: Peritoneal interleukin-10 increases with decrease in activated CD4+ T lymphocytes in women with endometriosis. Hum Reprod. 1997, 12 (11): 2528-2533. 10.1093/humrep/12.11.2528.

    Article  CAS  PubMed  Google Scholar 

  36. Ho HN, Wu MY, Chao KH, Chen CD, Chen SU, Chen HF, Yang YS: Decrease in interferon gamma production and impairment of T-lymphocyte proliferation in peritoneal fluid of women with endometriosis. Am J Obstet Gynecol. 1996, 175 (5): 1236-1241. 10.1016/S0002-9378(96)70034-0.

    Article  CAS  PubMed  Google Scholar 

  37. Keenan JA, Chen TT, Chadwell NL, Torry DS, Caudle MR: IL-1 beta, TNF-alpha, and IL-2 in peritoneal fluid and macrophage-conditioned media of women with endometriosis. Am J Reprod Immunol. 1995, 34 (6): 381-385.

    Article  CAS  PubMed  Google Scholar 

  38. Harada T, Yoshioka H, Yoshida S, Iwabe T, Onohara Y, Tanikawa M, Terakawa N: Increased interleukin-6 levels in peritoneal fluid of infertile patients with active endometriosis. Am J Obstet Gynecol. 1997, 176 (3): 593-597. 10.1016/S0002-9378(97)70553-2.

    Article  CAS  PubMed  Google Scholar 

  39. Punnonen J, Teisala K, Ranta H, Bennett B, Punnonen R: Increased levels of interleukin-6 and interleukin-10 in the peritoneal fluid of patients with endometriosis. Am J Obstet Gynecol. 1996, 174 (5): 1522-1526. 10.1016/S0002-9378(96)70600-2.

    Article  CAS  PubMed  Google Scholar 

  40. Hornung D, Bentzien F, Wallwiener D, Kiesel L, Taylor RN: Chemokine bioactivity of RANTES in endometriotic and normal endometrial stromal cells and peritoneal fluid. Mol Hum Reprod. 2001, 7 (2): 163-168. 10.1093/molehr/7.2.163.

    Article  CAS  PubMed  Google Scholar 

  41. Ponce C, Torres M, Galleguillos C, Sovino H, Boric MA, Fuentes A, Johnson MC: Nuclear Factor kappa B (NF{kappa}B) pathway and Interleukin-6 (IL-6) are affected in eutopic endometrium of women with endometriosis. Reproduction. 2009

    Google Scholar 

  42. Seto-Young D, Avtanski D, Strizhevsky M, Parikh G, Patel P, Kaplun J, Holcomb K, Rosenwaks Z, Poretsky L: Interactions among Peroxisome Proliferator Activated Receptor-{gamma}, Insulin Signaling Pathways, and Steroidogenic Acute Regulatory Protein in Human Ovarian Cells. J Clin Endocrinol Metab. 2007, 92 (6): 2232-2239. 10.1210/jc.2006-1935.

    Article  CAS  PubMed  Google Scholar 

  43. Bulun SE, Yang S, Fang Z, Gurates B, Tamura M, Sebastian S: Estrogen production and metabolism in endometriosis. Ann N Y Acad Sci. 2002, 955: 75-85. 10.1111/j.1749-6632.2002.tb02767.x.

    Article  CAS  PubMed  Google Scholar 

  44. Tamaya T, Motoyama T, Ohono Y, Ide N, Tsurusaki T, Okada H: Steroid receptor levels and histology of endometriosis and adenomyosis. Fertil Steril. 1979, 31 (4): 396-400.

    CAS  PubMed  Google Scholar 

  45. Janne O, Kauppila A, Kokko E, Lantto T, Ronnberg L, Vihko R: Estrogen and progestin receptors in endometriosis lesions: comparison with endometrial tissue. Am J Obstet Gynecol. 1981, 141 (5): 562-566.

    Article  CAS  PubMed  Google Scholar 

  46. Lyndrup J, Thorpe S, Glenthoj A, Obel E, Sele V: Altered progesterone/estrogen receptor ratios in endometriosis. A comparative study of steroid receptors and morphology in endometriosis and endometrium. Acta Obstet Gynecol Scand. 1987, 66 (7): 625-629. 10.3109/00016348709022068.

    Article  CAS  PubMed  Google Scholar 

  47. Bergqvist A, Ferno M: Estrogen and progesterone receptors in endometriotic tissue and endometrium: comparison according to localization and recurrence. Fertil Steril. 1993, 60 (1): 63-68.

    CAS  PubMed  Google Scholar 

  48. Yamashita Y, Asano M, Morishima S, Fujino K, Terai Y, Ohmichi M: Mitochondrial gene expression in granulosa cells of severe endometriosis with in vitro fertilization and embryo transfer. Fertil Steril. 2007, 88 (6): 1703-1705. 10.1016/j.fertnstert.2007.01.111.

    Article  CAS  PubMed  Google Scholar 

  49. Kao LC, Germeyer A, Tulac S, Lobo S, Yang JP, Taylor RN, Osteen K, Lessey BA, Giudice LC: Expression Profiling of Endometrium from Women with Endometriosis Reveals Candidate Genes for Disease-Based Implantation Failure and Infertility. Endocrinology. 2003, 144 (7): 2870-2881. 10.1210/en.2003-0043.

    Article  CAS  PubMed  Google Scholar 

  50. Sherwin JRA, Sharkey AM, Mihalyi A, Simsa P, Catalano RD, D'Hooghe TM: Global gene analysis of late secretory phase, eutopic endometrium does not provide the basis for a minimally invasive test of endometriosis. Hum Reprod. 2008, 23 (5): 1063-1068. 10.1093/humrep/den078.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Natural Science Foundation of China (grant numbers 30871782 and 30671492), National High Technology Research and Development Program of China (863 Program) (grant number 2006AA10Z1E3).

Author information

Authors and Affiliations

  1. School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, PR China

    Hongbo Zhao, Qishan Wang, Chunyan Bai, Kan He & Yuchun Pan

  2. Shanghai Key Lab for Veterinary Biotechnology, 200240, Shanghai, PR China

    Hongbo Zhao, Qishan Wang, Chunyan Bai, Kan He & Yuchun Pan

Authors
  1. Hongbo Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Qishan Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Chunyan Bai
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Kan He
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yuchun Pan
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Yuchun Pan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HBZ and YCP conceived and designed the study. HBZ contributed the data analysis and drafted the manuscript. KH performed data collection. QSW and CYB mined the literature and manuscript drafting. All authors read and approved the final manuscript.

Electronic supplementary material

12958_2009_566_MOESM1_ESM.xls

Additional file 1: Table 5: Pathway analysis of ovarian endometriosis data sets. The data provided represent the list of significant pathways identified by GSEA in the 3 ovarian endometriosis data sets. (XLS 110 KB)

12958_2009_566_MOESM2_ESM.xls

Additional file 2: Table 6: Pathway analysis of peritoneal endometriosis data sets. The data provided represent the list of significant pathways identified by GSEA in the 2 peritoneal endometriosis data sets. (XLS 61 KB)

12958_2009_566_MOESM3_ESM.xls

Additional file 3: Table 7: Pathway analysis between the timing of the uterine cycle. The data provided represent the list of significant pathways identified by GSEA between the timing of the uterine cycle. (XLS 33 KB)

12958_2009_566_MOESM4_ESM.xls

Additional file 4: Table 8: Pathway analysis in the specific phenotype of HEECs from endometriosis. The data provided represent the list of significant pathways identified by GSEA involved in the specific phenotype of HEECs from endometriosis. (XLS 35 KB)

12958_2009_566_MOESM5_ESM.xls

Additional file 5: Table 9: Common significant pathways in ovarian endometriosis data sets. The data provided represent the list of common significant pathways from GSEA of the 3 ovarian endometriosis datasets. (XLS 28 KB)

12958_2009_566_MOESM6_ESM.xls

Additional file 6: Table 10: Common significant pathways in peritoneal endometriosis data sets. The data provided represent the list of common significant pathways from GSEA of the 2 peritoneal endometriosis datasets. (XLS 22 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Zhao, H., Wang, Q., Bai, C. et al. A cross-study gene set enrichment analysis identifies critical pathways in endometriosis. Reprod Biol Endocrinol 7, 94 (2009). https://doi.org/10.1186/1477-7827-7-94

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1477-7827-7-94

Keywords

Reproductive Biology and Endocrinology

ISSN: 1477-7827

Contact us