Skip to main content

Comparing gene expression in deep infiltrating endometriosis with adenomyosis uteri: evidence for dysregulation of oncogene pathways



The pathogenesis of deep infiltrating endometriosis (DIE) is poorly understood. It is considered a benign disease but has histologic features of malignancy, such as local invasion or gene mutations. Moreover, it is not clear whether its invasive potential is comparable to that of adenomyosis uteri (FA), or whether it has a different biological background. Therefore, the aim of this study was to molecularly characterize the gene expression signatures of both diseases in order to gain insight into the common or different underlying pathomechanisms and to provide clues to pathomechanisms of tumor development based on these diseases.


In this study, we analyzed formalin-fixed and paraffin-embedded tissue samples from two independent cohorts. One cohort involved 7 female patients with histologically confirmed FA, the other cohort 19 female patients with histologically confirmed DIE. The epithelium of both entities was microdissected in a laser-guided fashion and RNA was extracted. We analyzed the expression of 770 genes using the nCounter expression assay human PanCancer (Nanostring Technology).


In total, 162 genes were identified to be significantly down-regulated (n = 46) or up-regulated (n = 116) in DIE (for log2-fold changes of < 0.66 or > 1.5 and an adjusted p-value of < 0.05) compared to FA. Gene ontology and KEGG pathway analysis of increased gene expression in DIE compared to FA revealed significant overlap with genes upregulated in the PI3K pathway and focal adhesion signaling pathway as well as other solid cancer pathways. In FA, on the other hand, genes of the RAS pathway showed significant expression compared to DIE.


DIE and FA differ significantly at the RNA expression level: in DIE the most expressed genes were those belonging to the PI3K pathway, and in FA those belonging to the RAS pathway.


Deep infiltrating endometriosis (DIE) and adenomyosis (FA) are very common benign gynecological conditions in women of childbearing age: 10–15% of women undergoing laparoscopy for benign reasons are generally found to have endometriosis [1], and in patients with infertility the incidence is even higher at 30–50% [2]. DIE is thought to occur in at least 20% of women with pelvic endometriosis [3]. FA affects 19.5% of women of childbearing age [4] [5]. In hysterectomy specimens, the incidence of diagnosed FA is 10–35% [6].

Both conditions are defined by the presence of endometrial glands and stroma either in the myometrium or outside the uterus. Originally DIE was even called “adenomyosis externa” [7].

There is a major discussion, if FA and especially DIE are related diseases or not, as in women with both entities the phenotype of FA appears to be related to the severity of endometriosis, particularly as women with DIE had a significantly higher frequency of focal adenomyosis in the external myometrium than patients with ovarian endometriosis [8]. Clinically however, the two conditions differ markedly. FA is commonly associated with dysmenorrhea, infertility, repeated implantation failure and pregnancy loss [9,10,11]. Multiparity and previous uterine surgery are discussed as risk factors for this condition [12,13,14]. DIE, on the other hand, is known to cause extensive adhesions up to complete obliteration of the Douglas space, as well as constriction of affected organs such as the bowel and bladder, and thus can cause not only dysmenorrhea but also chronic abdominal pain and, in severe cases, bowel obstruction up to the point of an ileus.

Although histologically benign, both DIE and FA are characterized by their propensity for local tissue invasion and resistance to apoptosis [15]. Notably, DIE has been described as a “benign tumor” [16]. Recent work using next-generation sequencing (NGS) has demonstrated driver mutations in cancer associated genes such as PIK3CA, ARID1A, PPP2R1A and KRAS in both, ovarian endometrioma [17] and DIE [15].

The presence of PIK3CA- or KRAS-mutated clones in histologically normal uterine endometrium in endometriosis [17] but also in patients without endometriosis has also been demonstrated [18], so the theory of the cellular origin of endometriosis requires further investigation. In contrast, the discovery of identical mutations in the KRAS gene in coexisting adenomyotic and endometriotic lesions in several patients [11], supports the theory of a common pathogenesis of adenomyosis uteri and endometriosis and a common molecular mechanism in these diseases [11, 19].

In this study, we aimed to further characterize the molecular mechanisms involved in FA and especially DIE to find molecular similarities and differences in both diseases. To this end, we analyzed cancer-related signaling pathways at the gene expression level using a nanostring gene panel encompassing the major signaling pathways of carcinogenesis using epithelial cells of FA and DIE. Any relevant genomic changes at the DNA level should be reflected in their gene expression and provide insight into the common or different underlying pathomechanisms of the two diseases and provide clues to the pathomechanisms of tumor development in these diseases and treatment options.

Patients and methods

Study population

For this study, we collected formalin-fixed and paraffin-embedded (FFPE) tissue samples for the analysis of two independent cohorts of patients with DIE or FA. Patients underwent surgery at the University Hospital, Heidelberg or cooperating clinics and the samples were histologically examined and assessed at the Dept. of Pathology of the University Hospital, Heidelberg between 2003 and 2018. Clinical records and histology were reviewed. Exclusion criteria were histological indications of cancer or dysplasia or a lesion size too small to gain sufficient material for further analysis. The samples were provided by the Tissue Bank at the National Center for Tumor Diseases (Heidelberg, Germany) in concordance with the Ethics Committee of the University of Heidelberg (approval No. S-362/2017).

Staining and laser microdissection (LMD)

For the RNA extraction, FFPE tissue blocks from FA and DIE were selected after reviewing all original tissue slides and were recut for hematoxylin & eosin sections, to be used for reference and to determine the lesion size. RNA was extracted using 10–20 FFPE slides for each entity.

For mounting on Zeiss 1.0 PEN slides (Carl Zeiss, Oberkochen, Germany) and for better adhesion of the tissue to the membrane, slides were irradiated with UV light/254nm for 30 min before, FFPE tissue blocks were cut at 8 μm thickness and incubated overnight at 37 °C. They were dewaxed in xylene (100%), rehydrated through decreasing concentrations of ethanol (100, 95, 75%), stained in 1% cresyl violet acetate (Sigma-Aldrich, Taufkirchen, Germany) and again dehydrated in increasing ethanol concentrations (75, 95, 100%). After that, tissue sections were dried and stored at 4 °C. Using a ZEISS PALM LMD laser capture microdissection unit, regions of interest (epithelium of the adenomyosis or epithelium of deep infiltrating endometriosis) were microdissected. The isolated tissue fragments corresponded to an area of approximately 20.000.000 μm² for each sample. They were collected in AdhesiveCap 500 opaque tubes (Carl Zeiss) and stored at − 20 °C until further processing (Fig. 1).

Fig. 1
figure 1

Light microscopy of deep infiltrating endometriosis of the rectum (A) and of adenomyosis uteri (B) stained with cresylviolet (10x). In each case the second image in the row shows the marking of the tissue for laser dissection [2]. The third image shows the tissue after laser dissection [3]

RNA isolation

Extraction of the total RNA from microdissected tissue samples was performed using the AllPrep DNA/RNA FFPE Kit (Qiagen, Venlo, the Netherlands) according to the manufacturer’s protocol. They were quantified with the Nanodrop ND-1000 spectrophotometer (NanoDrop Technologies, Rockland, DE, USA).

Gene expression analysis

We analyzed the expression of 770 genes (Codeset: Human PanCancer Pathways) and hybridisation counts were measured using the nCounter technology (both by Nanostring™ Technology, Seattle, WA). A minimum of approximately 50 ng of total RNA was used. Hybridization time per cartridge was 16 h before measurement. According to the manufacturer’s protocol, the examined genes were attached to specific tag sequences and hybridized for 16 h at 65 °C to a capture/reporter probe pair equipped with a fluorescent barcode. These gene-specific barcodes were then detected by the nCounter Digital Analyzer providing count of genes. No cases were excluded.

The 770 genes codeset included 730 genes from 13 canonical pathways (e.g., cell cycle, chromatin modeling, apoptosis, MAPK, and PI3K) and 40 housekeeping genes [20]. The raw data were pre-analyzed for consistency using the manufacturer’s software (nSolver version 4.0). The geNorm pairwise variation statistics was used for stepwise selection of normalization genes from the housekeeping genes [21]. Six genes with minimal pairwise variation statistics were finally selected for normalization (TLK2, VPS33B, TMUB2, C10orf76, SLC4A1AP, ERCC3).


Differential expression analysis was carried out using a linear data model in limma [22, 23], and nominal p-values were corrected for multiple comparisons using Benjamini and Hochberg’s method [24]. All genes with an adjusted false discovery rate (FDR) of p < 0.05 and fold change of < 0.66 or > 1.5 were considered differentially expressed. Differentially expressed genes (DEGs) were subjected to functional annotation and clusterization using DAVID Bioinformatics Resources (version 6.8, [25, 26]) after conversion of gene symbols to Entrez IDs and uploading to DAVID using the “RDAVIDWebService” BioConductor library [27]. Basal cytokeratin co-regulated genes were identified using DAVID analysis in the “Biological Process” category and the KEGG pathway enrichment function with a significance threshold of 0.05. The p-values of selected GO terms were corrected using Benjamini-Hochberg correction and described as adjusted p-values [24]. Otherwise, differences between samples were tested using Wilcoxon signed-rank test, and correlation was tested using Spearman’s rank correlation test, and p-value of 0.05 was considered significant. All statistical calculations were done using R version 4.0 [28]. For visualization, the R packages ComplexHeatmap and Ggplot2 were utilized.


Study population

The study cohort included 19 female patients with DIE. DIE lesions were defined as histologically confirmed endometriosis with infiltrative growth with a depth of more than 5 mm into the wall of pelvic organs, e.g. in the bowel or bladder (Tables 1 and 2). The control group included seven patients with FA and infiltration of less than half of myometrium (superficial) or more than one half of myometrium (deep adenomyosis) (Tables 1 and 2).

Table 1 Cohort description including age (years), symptoms/presentation in clinic, site of manifestation, and severity according to rASRM/Enzian classification for each of the 26 patients enrolled in the study
Table 2 Overview of clinical and pathological patient data

The mean age of the patients suffering from FA was 48.4 years, as compared to 33.4 years for patients suffering from DIE. Patients with adenomyosis had undergone surgery due to bleeding disorders, while DIE was resected for different reasons mainly including acute and chronic pain or incipient intestinal obstruction. The GI-tract (rectum, recto-sigmoid, colon and ileum, 15 cases) was mostly affected by DIE, followed by the vagina (7 cases) and the bladder (1 case), including overlapping sites. FA included cases with superficial and deep infiltration.

Dysregulated genes in DIE vs. adenomyosis uteri

A mathematical model was constructed for the analysis of differential gene expression in FA and DIE. When using a threshold for fold changes (FC) < 0.66 or > 1.5, a total of 162 genes were identified that were up- or downregulated (adj. p < 0.05). This analysis included significantly more genes with upregulation in DIE (116 genes), as compared to 46 genes with relative downregulation in DIE, compared to adenomyosis (FC < 0.66). When using a stricter threshold of significance (p < 0.001 was used), 15 genes were upregulated in DIE, and only one gene (FDZ2) was upregulated in adenomyosis uteri (Table 3; Fig. 2). With regards to the functional properties, no specific pathway could be assigned, and therefore a separate gene ontology analysis was performed (see below).

Table 3 Genes with significant relative upregulation (n = 15) or downregulation (n = 1) in DIE vs. adenomyosis (> 1.5 fold, Benjamini–Hochberg adjusted P ≤ 0.001)
Fig. 2
figure 2

Volcano plot showing genes with most significant dysregulation in DIE vs. adenomyosis (adj. p-adjusted < 0.001 and log2 FC > 1.5). 15 genes were highly significantly upregulated in DIE, and one gene (FZD2) was downregulated (log2 FC < 0.66).

In order to relate dysregulated genes to clinical characteristics of adenomyosis and DIE, an unsupervised heatmap was constructed. Here, clustering revealed a clear separation of DIE and adenomyosis cases with only one DIE case clustering within adenomyosis (Fig. 3). But generally, samples from DIE had generally higher pathway activity scores than samples from adenomyosis (FA). In this analysis, one larger gene group with upregulation in DIE could be separated from a smaller gene group with upregulation in adenomyosis (Fig. 3). However, in this clustering no correlation of gene expression with clinical characteristics (organ, BMI, depth of adenomyosis) was evident.

Fig. 3
figure 3

Heatmap of supervised hierarchical clustering of differential genes (adj. p < 0.05, log2 fold change > 1.5 or < 0.66 and adjusted p-value < 0.05, n = 162) for the FA versus DIE group. This includes 116 upregulated genes in DIE and 46 genes with upregulation in adenomyosis. Gene expression with cases of adenomyosis (FA) is clearly distinct from cases with deep infiltration endometriosis (DIE).

Pathways with activation in adenomyosis uteri and DIE

In order to analyze the functional properties of dysregulated genes in both diseases we performed gene ontology analysis using the KEGG pathway analysis. This analysis revealed upregulating of several signaling pathways in DIE, and interestingly, the PIK3CA pathway was most significantly upregulated. Other gene ontology groups included pathways involved in virus infection, focal adhesion, endocrine resistance and malignancy (Fig. 4a).

Fig. 4
figure 4

(a) KEGG enrichment analyses for differentially upregulated genes (p < 0.05, n = 116), 12 most significantly upregulated pathways are shown. This analysis revealed upregulating of signaling pathways in DIE, most significantly the PI3K pathway, but also pathways involved in virus infection, focal adhesion, endocrine resistance and malignancy. (b) Same KEGG analysis, for differentially downregulated genes (p < 0.05, n = 46), 4 most significant pathways (p < 0.0001) are shown. Here, RAS, PI3K-AKT, RAP1 and Calcium signaling pathways are significant

The identification of virus infection pathways was an unexpected finding, but further analysis revealed that these virus-pathway associated genes had 12 genes in common, and 11 of these genes were also common to the PIK3CA pathway. Therefore, it is believed that the PIK3CA pathway upregulation is the root cause for showing virus-pathway associated pathways in this analysis, and that this analysis does not point to a virus related cause of DIE. Upregulated pathways in adenomyosis included RAS, PI3K-AKT, RAP1 and calcium signaling pathways (Fig. 4b).


To our knowledge, this is the first study comparing isolated epithelium cells of deep infiltrating endometriosis with epithelium of adenomyosis uteri using the nanostring technology. This is of particular interest, as adenomyosis and endometriosis lesions are often surrounded by many stromal and inflammatory cells, which cause some blurring of studies at the molecular level [29]. Anglesio et al., Inoue et al. and Moore et al. have recently shown that the somatic mutation occurs in the epithelial component of DIE [15], as well as of FA [11] and in the histologically normal endometrium of healthy patients [18]. In the latter studies, analysis of the laser microdissected epithelium has been shown to yield promising results [11, 15, 17, 18, 30]. Therefore, this method was also used in the current work.

In our studies, we show for the first time that the epithelium of DIE and FA differ significantly at the RNA expression level. Interestingly, these differences in RNA expression between both entities are independent of the site of DIE manifestation and the body mass index.

Looking first at the KEGG analysis of genes whose expression is elevated in DIE compared to FA, it appears that the PI3K pathway is significantly activated. This is consistent with the results of whole-exome studies by Anglesio 2017 and Suda 2018, which detected somatic driver mutations at the DNA level in the PIK3CA gene in the epithelium of DIE lesions [15] and in the epithelium of ovarian endometriosis, as well as in the eutopic endometrium of healthy patients [17], among others. Since most PIK3CA mutations in cancers show gain of function and growth advantages [31], it was concluded that the presence of the same mutations in endometriotic epithelial cells has functional significance in the pathogenesis of the disease [17, 32]. We show that this may also be reflected in the activation of the pathway at the RNA level. Previous work has also indicated dysregulation of the PI3K pathway in endometriosis: Yin et al. 2012 demonstrated an increase in pAKT, albeit in stromal cells from endometriomas, compared to cells from the eutopic endometrium of healthy women, and Guo et al. 2015 also showed that phosphorylated mTor is increased in ectopic endometrial lesions compared to eutopic endometrium from endometriosis patients [33, 34]. Also, in a recent work by Madanes, the authors are able to demonstrate increased expression of PI3K, reduced expression of PTEN, and increased levels of pAkt in the ectopic and eutopic endometrium of patients with peritoneal endometriosis [35].

The PI3K-AKT-mTOR pathway is one of the most frequently dysregulated signaling pathways in carcinoma diseases [36]. Its significant activation in DIE compared to FA may explain the different behavior of the two entities under study.

We believe, that the fact, that the KEGG analysis revealed activation of virus associated gene groups (HPV, KSAHV and CMV) is due to the large intersection of genes activated in these gene groups with genes that are significantly activated in the PI3K pathway, rather than a virus-related cause of DIE, as discussed in detail above.

The focal adhesion pathway plays an essential role in cell motility, cell proliferation and cell differentiation. Its increased activation in DIE compared to FA could also explain the more progressive behavior of DIE. Interestingly, in a recent analysis of the proteome of the eutopic endometrium of endometriosis patients, Méar et al. also demonstrated an increased activation of the PI3K pathway and the focal adhesion pathway compared to healthy controls [37].

The KEGG analysis of genes that are downregulated in DIE and upregulated in adenomyosis compared to DIE shows that the RAS pathway in particular is upregulated in adenomyosis. This fits well with the findings of Inoue et al., who demonstrated a mutation in the KRAS gene with consecutive activation of KRAS in 37.1% of adenomyosis cases [11]. They detected a mutation in PI3KCA in only two of 70 patients, which may explain the increased expression of PI3K pathway-associated genes in DIE in our study compared with FA. In contrast, in patients with both adenomyosis and endometriosis lesions, Inoue et al. were able to detect the same KRAS mutation in both lesions [11]. By comparing both tissues in our study and looking at the relative gene expressions of both entities, the absolute activation of the RAS pathway is shown to be lower in endometriosis patients, which would explain why it is not detected in the KEGG analysis. In direct comparison of the two entities at the level of gene expression however, the PI3K pathway appears to be the dominant pathway in DIE and the RAS pathway in FA.

This could explain, for example, the different sensitivity of the two entities to certain therapeutic approaches. While DIE responds well to therapy with the progestin dienogest, progesterone resistance is often described in adenomyosis patients, which according to Inoue et al. is due to the KRAS mutations [11, 32].

If we look at the changes of single genes, it is striking that Frizzled class receptor 2 (FDZ2) is the only gene that is significantly upregulated in adenomyosis compared to DIE. FZD2 is discussed as an important trigger of TGF-ß induced epithelial–mesenchymal transition (EMT) [38] and cell migration [39]. Accordingly, the induction of EMT and ultimate fibrosis by TGF-β1 appears to play a critical role in the pathogenesis of adenomyosis [40]. In addition, EMT promoted by FZD2 also plays an important role in the metastasis of endometrial cancer [41], therefore suggesting that the invasive behavior of epithelial cells in FA has a cancerogenic aspect.

The two genes upregulated most in DIE compared to FA, Heat Shock Protein Family B (Small) Member 1 (HSPB1) and Calpain 2 (CAPN2), have not yet been described in endometriosis, despite the fact, that they play a role in proliferation and invasion in various solid tumors [42,43,44,45] and may serve this same function in DIE.

Catenin Beta 1 (CTNNB1), on the other hand, which is likewise upregulated in DIE epithelial cells compared to cells from adenomyosis, is discussed as a key factor in the regulation of proliferation and invasion of endometriosis [46, 47]. Since both papers demonstrate an upregulation of CTNNB1 in endometrial stromal cells in endometriosis lesions, we can raise the question of an additional important role for CTNNB1 action in the epithelium of endometriosis lesions.

Limitations of this study may include the following: [1] The group of patients suffering from deep infiltrating endometriosis is younger than the group of patients with adenomyosis who underwent hysterectomy. However, some of the younger patients were also on therapy with a GnRH analogue, which hormonally corresponds to a menopausal status and therefore attenuates any age-related differences. [2] By comparing adenomyosis and DIE without comparison to normal endometrium, we can only show relative expression differences, but not absolute differences compared to healthy tissue. Nevertheless, especially in view of the frequent co-occurrence of the two entities and the presumed common molecular origin, we consider our results at the gene expression level to be of further value, particularly with regard to possible different therapeutic approaches. [3] The sample size is relatively small with 7 patients with adenomyosis and 19 patients with DIE. Further studies with a larger number of patients to evaluate the possible influence of age, therapy concept or site of manifestation of DIE would be useful.


Deep infiltrating endometriosis and adenomyosis uteri differ significantly at the RNA expression level: for deep infiltrating endometriosis, the genes most expressed were those belonging to the PI3K pathway, and for adenomyosis, those belonging to the RAS pathway.

Data Availability

Data or material are available on reasonable request.



AT-rich interaction domain1A


Body mass index


chromosome 10, open-reading frame 76


Calpain 2




Catenin Beta 1


differentially expressed genes


deep infiltrating endometriosis


deoxyribonucleic acid


epithelial mesenchymal transition


ERCC excision repair 3


Adenomyosis uteri


Fold change


False discorvery rate


formalin-fixed and paraffin-embedded


Frizzled class receptor 2


Gastrointestinal tract


Gene ontology


Human papillomavirus


Heat Shock Protein Family B (Small) Member 1


KRAS Proto-oncogene, GTPase


kaposi sarcoma associated herpesvirus infection


laser microdissection


next-generation sequencing


phosphorylated AKT Serine/Threonine Kinase


Phosphatidylinositol 3-kinase


Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Alpha


Protein phosphatase 2, regulatory subunit A


Phosphatase and Tensin Homolog


Ras-related protein 1


Rat sarcoma


Revised American Society for Reproductive Medicine classification


Ribonucleic acid


Solute Carrier Family 4 Member 1 Adaptor Protein


Transforming growth factor ß


Tousled Like Kinase 2


Transmembrane and Ubiquitin Like Domain Containing 2


Vacuolar Protein Sorting-Associated Protein 33B


  1. Giudice LC. Clinical practice. Endometriosis. N Engl J Med. 2010;362(25):2389–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. 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 

  3. Chapron C, Fauconnier A, Vieira M, Barakat H, Dousset B, Pansini V et al. Anatomical distribution of deeply infiltrating endometriosis: surgical implications and proposition for a classification.Human Reproduction. 2003;18:157–61.

  4. Devlieger R, D’Hooghe T, Timmerman D. Uterine adenomyosis in the infertility clinic. Hum Reprod Update. 2003;9(2):139–47.

    Article  PubMed  Google Scholar 

  5. Garcia-Solares J, Donnez J, Donnez O, Dolmans MM. Pathogenesis of uterine adenomyosis: invagination or metaplasia? Fertil Steril. 2018;109(3):371–9.

    Article  PubMed  Google Scholar 

  6. Yeniel O, Cirpan T, Ulukus M, Ozbal A, Gundem G, Ozsener S, et al. Adenomyosis: prevalence, risk factors, symptoms and clinical findings. Clin Exp Obstet Gynecol. 2007;34(3):163–7.

    CAS  PubMed  Google Scholar 

  7. Cornillie FJ, Oosterlynck D, Lauweryns JM, Koninckx PR. Deeply infiltrating pelvic endometriosis: histology and clinical significance. Fertil Steril. 1990;53(6):978–83.

    Article  CAS  PubMed  Google Scholar 

  8. Chapron C, Tosti C, Marcellin L, Bourdon M, Lafay-Pillet MC, Millischer AE, et al. Relationship between the magnetic resonance imaging appearance of adenomyosis and endometriosis phenotypes. Hum Reprod. 2017;32(7):1393–401.

    Article  PubMed  Google Scholar 

  9. Vannuccini S, Clifton VL, Fraser IS, Taylor HS, Critchley H, Giudice LC, et al. Infertility and reproductive disorders: impact of hormonal and inflammatory mechanisms on pregnancy outcome. Hum Reprod Update. 2016;22(1):104–15.

    Article  CAS  PubMed  Google Scholar 

  10. Hashimoto A, Iriyama T, Sayama S, Nakayama T, Komatsu A, Miyauchi A, et al. Adenomyosis and adverse perinatal outcomes: increased risk of second trimester miscarriage, preeclampsia, and placental malposition. J Matern Fetal Neonatal Med. 2018;31(3):364–9.

    Article  PubMed  Google Scholar 

  11. Inoue S, Hirota Y, Ueno T, Fukui Y, Yoshida E, Hayashi T, et al. Uterine adenomyosis is an oligoclonal disorder associated with KRAS mutations. Nat Commun. 2019;10(1):5785.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Parazzini F, Vercellini P, Panazza S, Chatenoud L, Oldani S, Crosignani PG. Risk factors for adenomyosis. Hum Reprod. 1997;12(6):1275–9.

    Article  CAS  PubMed  Google Scholar 

  13. Vercellini P, Vigano P, Somigliana E, Daguati R, Abbiati A, Fedele L. Adenomyosis: epidemiological factors. Best Pract Res Clin Obstet Gynaecol. 2006;20(4):465–77.

    Article  PubMed  Google Scholar 

  14. Templeman C, Marshall SF, Ursin G, Horn-Ross PL, Clarke CA, Allen M, et al. Adenomyosis and endometriosis in the California Teachers Study. Fertil Steril. 2008;90(2):415–24.

    Article  PubMed  Google Scholar 

  15. Anglesio MS, Papadopoulos N, Ayhan A, Nazeran TM, Noe M, Horlings HM, et al. Cancer-Associated mutations in endometriosis without Cancer. N Engl J Med. 2017;376(19):1835–48.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Gordts S, Koninckx P, Brosens I. Pathogenesis of deep endometriosis. Fertil Steril. 2017;108(6):872–85. e1.

    Article  PubMed  Google Scholar 

  17. Suda K, Nakaoka H, Yoshihara K, Ishiguro T, Tamura R, Mori Y, et al. Clonal expansion and diversification of Cancer-Associated mutations in endometriosis and normal endometrium. Cell Rep. 2018;24(7):1777–89.

    Article  CAS  PubMed  Google Scholar 

  18. Moore L, Leongamornlert D, Coorens THH, Sanders MA, Ellis P, Dentro SC, et al. mutational Landsc normal Hum endometrial epithelium. 2020;580(7805):640–6.

    CAS  Google Scholar 

  19. Koninckx PR, Ussia A, Adamyan L, Wattiez A, Gomel V, Martin DC. Pathogenesis of endometriosis: the genetic/epigenetic theory. Fertil Steril. 2019;111(2):327–40.

    Article  CAS  PubMed  Google Scholar 

  20. Cesano A. nCounter((R)) PanCancer Immune Profiling Panel (NanoString Technologies, Inc., Seattle, WA). J Immunother Cancer. 2015;3:42.

  21. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3(7):RESEARCH0034.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Smyth GK. Limma: Linear Models for microarray data. In: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W, editors. Bioinformatics and Computational Biology solutions using R and Bioconductor. Statistics for Biology and Health. Berlin, Heidelberg, New York: Springer; 2005. pp. 397–420.

    Google Scholar 

  23. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Benjamini Y, Hochberg Y. Controlling the false Discovery Rate - a practical and powerful Approach to multiple testing. J R Stat Soc B. 1995;57(1):289–300.

    Google Scholar 

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

    Article  PubMed  Google Scholar 

  26. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.

    Article  PubMed  Google Scholar 

  27. Fresno C, Fernandez EA. RDAVIDWebService: a versatile R interface to DAVID. Bioinformatics. 2013;29(21):2810–1.

    Article  CAS  PubMed  Google Scholar 

  28. R Core Team. R: A language and environment for statistical computing. Date of access: 20.05.2020. Available from:

  29. 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):382–403.

    Article  CAS  PubMed  Google Scholar 

  30. Matsuzaki S, Canis M, Mage G. Use of laser capture microdissection in studying hormone-dependent diseases: endometriosis. Methods Mol Biol. 2009;590:295–306.

    Article  CAS  PubMed  Google Scholar 

  31. Gymnopoulos M, Elsliger M-A, Vogt PK. Rare cancer-specific mutations in PIK3CA show gain of function. Proceedings of the National Academy of Sciences. 2007;104.

  32. Bulun SE, Yildiz S, Adli M, Wei J-J. Adenomyosis pathogenesis: insights from next-generation sequencing.Human Reproduction Update. 2021;27:1086–97.

  33. Yin X, Pavone ME, Lu Z, Wei J, Kim JJ. Increased Activation of the PI3K/AKT Pathway Compromises Decidualization of Stromal Cells from Endometriosis.The Journal of Clinical Endocrinology \& Metabolism. 2012;97:E35–E43.

  34. Guo J, Gao J, Yu X, Luo H, Xiong X, Huang O. Expression of DJ-1 and mTOR in Eutopic and Ectopic Endometria of Patients with Endometriosis and Adenomyosis.Gynecologic and Obstetric Investigation. 2015;79:195–200.

  35. Madanes D, Bilotas MA, Bastón JI, Singla JJ, Meresman GF, Barañao RI et al. PI3K/AKT pathway is altered in the endometriosis patient’s endometrium and presents differences according to severity stage. Gynecol Endocrinol. 2020;36:436–40.

  36. Janku F, Yap TA, Meric-Bernstam F. Targeting the PI3K pathway in cancer: are we making headway?Nature Reviews Clinical Oncology. 2018;15:273–91.

  37. Méar L, Com E, Fathallah K, Guillot L, Lavigne R, Guével B et al. The Eutopic Endometrium Proteome in Endometriosis Reveals Candidate Markers and Molecular Mechanisms of Physiopathology.Diagnostics. 2022;12:419.

  38. Tuluhong D, Chen T, Wang J, Zeng H, Li H, Dunzhu W, et al. FZD2 promotes TGF-beta-induced epithelial-to-mesenchymal transition in breast cancer via activating notch signaling pathway. Cancer Cell Int. 2021;21(1):199.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Gujral TS, Chan M, Peshkin L, Sorger PK, Kirschner MW, MacBeath G. A noncanonical Frizzled2 pathway regulates epithelial-mesenchymal transition and metastasis. Cell. 2014;159(4):844–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Kay N, Huang CY, Shiu LY, Yu YC, Chang Y, Suen JL, et al. The Effects of Anti-TGF-beta1 on epithelial-mesenchymal transition in the pathogenesis of adenomyosis. Reprod Sci. 2020;27(9):1698–706.

    Article  CAS  PubMed  Google Scholar 

  41. Bian Y, Chang X, Liao Y, Wang J, Li Y, Wang K, et al. Promotion of epithelial-mesenchymal transition by Frizzled2 is involved in the metastasis of endometrial cancer. Oncol Rep. 2016;36(2):803–10.

    Article  CAS  PubMed  Google Scholar 

  42. Calderwood SK, Gong J. Heat shock proteins promote Cancer: it’s a Protection Racket. Trends Biochem Sci. 2016;41(4):311–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Alberti G, Vergilio G, Paladino L, Barone R, Cappello F, Conway de Macario E et al. The Chaperone System in Breast Cancer: Roles and Therapeutic Prospects of the Molecular Chaperones Hsp27, Hsp60, Hsp70, and Hsp90.Int J Mol Sci. 2022;23(14).

  44. Li P, Miao C, Liang C, Shao P, Wang Z, Li J. Silencing CAPN2 expression inhibited castration-resistant prostate Cancer cells Proliferation and Invasion via AKT/mTOR Signal Pathway. Biomed Res Int. 2017;2017:2593674.

    PubMed  PubMed Central  Google Scholar 

  45. Chen B, Tang J, Guo YS, Li Y, Chen ZN, Jiang JL. Calpains are required for invasive and metastatic potentials of human HCC cells. Cell Biol Int. 2013;37(7):643–52.

    Article  CAS  PubMed  Google Scholar 

  46. Quattrone F, Sanchez AM, Pannese M, Hemmerle T, Vigano P, Candiani M, et al. The targeted delivery of interleukin 4 inhibits development of endometriotic lesions in a mouse model. Reprod Sci. 2015;22(9):1143–52.

    Article  CAS  PubMed  Google Scholar 

  47. Wang X, Wu P, Zeng C, Zhu J, Zhou Y, Lu Y et al. Long Intergenic Non-Protein Coding RNA 02381 Promotes the Proliferation and Invasion of Ovarian Endometrial Stromal Cells through the miR-27b-3p/CTNNB1 Axis.Genes (Basel). 2022;13(3).

Download references


We would like to thank the working group of Prof. Dr. rer. nat. Stefanie Rössler for the many fruitful discussions and support during the evaluation as well as MTA Jutta Scheurer for the great support during the LMD. We would also like to thank Professor Stenzinger’s research group for continuous support of the project.


Open Access funding enabled and organized by Projekt DEAL. A.M. has received financial support for the project from the Kußmaul and Walter-Erb Foundation. The other authors received no financial support for the research and authorship of this article. For the publication fee we acknowledge financial support by Deutsche Forschungsgemeinschaft within the funding programme „Open Access Publikationskosten“ as well as by Heidelberg University

Author information

Authors and Affiliations



The study conception and design was conducted by A.M., P.S. and A.G. Material preparation, data collection and analysis were performed by A.M, F.K., H.O., M.K., J.J. and P.S. Supervision of the project was conducted by A.G. and T.S. Financial support and equipment were made available by T.S. The first draft of the manuscript was written by A.M., P.S. and A.G. and all authors commented on the following versions of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to A. Marshall.

Ethics declarations

Ethics approval and consent to participate

Patient consent was obtained at the time of hospitalization for the use of diagnostic materials for research. The tissue samples were provided by the Tissue Bank at the National Center for Tumor Diseases (Heidelberg, Germany) in concordance with the Ethics Committee of the University of Heidelberg (approval No. S-362/2017).

Consent for publication

All authors read and approved the manuscript in the present form.

Competing interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Marshall, A., Kommoss, K.F., Ortmann, H. et al. Comparing gene expression in deep infiltrating endometriosis with adenomyosis uteri: evidence for dysregulation of oncogene pathways. Reprod Biol Endocrinol 21, 33 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Gene expression analysis
  • Deep infiltrating endometriosis
  • Adenomyosis
  • PI3K pathway
  • RAS pathway