Comprehensive bioinformatics analysis of acquired progesterone resistance in endometrial cancer cell line

Background Progesterone resistance is a problem in endometrial carcinoma, and its underlying molecular mechanisms remain poorly understood. The aim of this study was to elucidate the molecular mechanisms of progesterone resistance and to identify the key genes and pathways mediating progesterone resistance in endometrial cancer using bioinformatics analysis. Methods We developed a stable MPA (medroxyprogesterone acetate)-resistant endometrial cancer cell subline named IshikawaPR. Microarray analysis was used to identify differentially expressed genes (DEGs) from triplicate samples of Ishikawa and IshikawaPR cells. PANTHER, DAVID and Metascape were used to perform gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and cBioPortal for progesterone receptor (PGR) coexpression analysis. GEO microarray (GSE17025) was utilized for validation. The protein–protein interaction network (PPI) and modular analyses were performed using Metascape and Cytoscape. Further validation were performed by real-time polymerase chain reaction (RT-PCR). Results In total, 821 DEGs were found and further analyzed by GO, KEGG pathway enrichment and PPI analyses. We found that lipid metabolism, immune system and inflammation, extracellular environment-related processes and pathways accounted for a significant portion of the enriched terms. PGR coexpression analysis revealed 7 PGR coexpressed genes (ANO1, SOX17, CGNL1, DACH1, RUNDC3B, SH3YL1 and CRISPLD1) that were also dramatically changed in IshikawaPR cells. Kaplan–Meier survival statistics revealed clinical significance for 4 out of 7 target genes. Furthermore, 8 hub genes and 4 molecular complex detections (MCODEs) were identified. Conclusions Using microarray and bioinformatics analyses, we identified DEGs and determined a comprehensive gene network of progesterone resistance. We offered several possible mechanisms of progesterone resistance and identified therapeutic and prognostic targets of progesterone resistance in endometrial cancer. Electronic supplementary material The online version of this article (10.1186/s12967-019-1814-6) contains supplementary material, which is available to authorized users.


Background
Endometrial cancer (EC) is the most common gynecologic malignancy in developed countries. According to cancer statistics in China, there were 634 per 100,000 newly diagnosed and 21.8 per 100,000 mortalities in 2017. EC had a significant upward trend in age-standardized incidence rates [1].
Approximately 80% of EC cases are type I endometrial cancer, which are highly associated with prolonged unopposed estrogen action and insufficient progesterone. Common progestin therapy drugs including medroxyprogesterone acetate (MPA) and megestrol acetate (MA), have been used as a conservative treatment for stage 1A low-grade EC patients who desire to preserve fertility or patients with comorbidities who are unsuitable candidates for operation. Recurrent endometrial cancer patients can also receive hormonal therapy. Some studies have indicated that synthetic progesterone is effective for stage 1A EC [2][3][4][5][6][7]. However, approximately 30% of early stage 1A EC patients never responded to treatment or only exhibit a temporary response when treated as primary therapy [2,7]. For advanced cases, the response rate is only 20-40% [8]. Although 70% of early patients respond to treatment, 57% of these patients experience recurrence. When treated by MPA again after diagnosis of recurrence, 63% were nonresponsive [7]. All of the above findings indicate that de novo or acquired progestin resistance is a major clinical problem.
Several studies have investigated the mechanism of progesterone resistance in EC, but the precise mechanism remains unknown. When progesterone resistance occurs, progestin can promote proliferation and invasiveness rather than inhibit growth and promote apoptosis of cancer cells [9]. Previous mechanistic explanations include imbalance of ER and PR subtypes [10], overexpression of EGFR and activation of TGF-EGFR signaling [11][12][13], activation of the PI3K/AKT/mTOR pathway [14,15], dysregulation of Nrf2-survivin and overexpression of survivin [16]. We previously developed a stable MPAresistant Ishikawa cell and confirmed that SIRT1/FoxO1/ SREBP-1 acts as a pathway targeting PR that is involved in the development of progesterone resistance in endometrial cancer cell [17]. However, current research does not sufficient explain the key genes that cause the down regulation of PR in progesterone resistance and molecular network involved in progesterone resistance.
Microarray analysis is a high-throughput approach that has been used more than 10 years. Bioinformatics analysis of microarray data, including gene clustering, gene ontology and pathway analysis, has shown great significance to identify potential key genes and pathways within a complex disease or biological process [18]. No microarray profiling study of progesterone resistance within endometrial cancer has been performed. In this work, we performed microarray analyses to comprehensively analyze the expression profile of progesterone resistance. In total, 821 DEGs, 7 PGR coexpressed genes, 8 hub genes and 4 MCODEs were found and further analyzed by GO, pathway enrichment and PPI analyses. We found that lipid metabolism, the immune system and inflammatory response and the extracellular environment may be crucial in progesterone resistance. These hub genes may help us identify novel biomarkers and treatment targets for progesterone resistance in the future.

Cell culture
Ishikawa EC cells were purchased from American Type Culture Collection (ATCC; Rockvile, MD, US). MPAresistant Ishikawa cell, which is referred to as IshikawaPR was established as previously described [17]. Ishikawa and IshikawaPR were routinely grown in RPMI 1640 (HyClone, USA) containing 10% fetal bovine serum (FBS) at 37 °C in a 5% CO 2 humidified atmosphere. IshikawaPR cell was routinely cultured in 10 μM MPA to maintain resistance.

Microarrays and bioinformatics analysis
Microarray analysis was performed using triplicate samples of parental Ishikawa and IshikawaPR cells. Total RNAs were prepared and subjected to Agilent Human 4 × 44 K Gene Expression Microarrays according to the manufacturer's instructions (Agilent) after quality control assessment. Sample labeling and microarray hybridization were performed using the Agilent One-Color Microarray-Based Gene Expression Analysis protocol (Agilent Technology). Array images were acquired and analyzed using the Agilent Feature Extraction software (version 11.0.1.1). Normalization and subsequent analytical data processing were performed using the Gene-Spring GX v12.1 software package. Genes for which at least 3 out of 6 samples exhibited detection were chosen for further data analysis. The raw data.tar has been submitted to GEO (Series GSE121367 https ://www.ncbi.nlm. nih.gov/geo/query /acc.cgi?&acc=GSE12 1367) and will be released until Dec 31, 2019.
Differential expressed genes (DEG) with statistical significance were identified through volcano plot filtering. The thresholds for DEG were |logFC| ≥ 4.0 and p value < 0.05 (Additional file 1:

PGR coexpression analysis and validation by bioinformatics
Assessment of the coexpression genes of PGR was performed using the cBioPortal database (http://www. cbiop ortal .org). The data obtained were RNA-Seq data from TCGA database that included 549 endometrial cancer tissues. Pearson's correlation score and Spearman score (≥ 0.3 was considered positively correlated and ≤ − 0.3 was considered negatively correlated with PGR) were used to select PGR coexpressed genes. To predict the target genes that were changed in Ishi-kawaPR, we use FunRich to identify the overlapping genes between DEGs and PGR coexpressed genes. For validation of target genes, the gene expression profile result, GSE17025, deposited by Day et al. [19] was used. The gene expression profile has 91 Stage I endometrial cancer patients and 12 postmenopausal healthy tissues. Then we calculated Pearson and Spearman score between target genes and PGR and compared the expression of target genes among healthy and different type of endometrial cancer tissues (p < 0.05 as cut-off criterion) using GraphPad Prism 7.0 and SPSS 22.0 software. Kaplan-Meier curves for target genes were generated with the online tool Kaplan-Meier Plotter (http://www.kmplo t.com/). A total of 542 RNAseq data samples of uterine corpus endometrial carcinoma were interrogated. The patients were split into 2 groups (high vs. low) based on the expression level.

Validation of target genes by real-time PCR
Total RNA was extracted from Ishikawa and Ishi-kawaPR cells using TRIzol reagent (Invitrogen, Carlsbad, CA) according to the manufacture' instructions. Total RNA (3 μg) was reverse transcribed using the M-MLV reverse transcriptase (Cat no. C28025-011, Invitrogen, China). Then RNA expression level of detected genes were quantified using an ABI Prism 7500 Sequence Detection System (Applied Biosystems, USA) with SYBR Green Master Mix (Takara, Japan) in a 20 µl reaction mixture. The primers were synthesized by Sangon Biotech Corporation and presented in Additional file 2: Table S2. Experiments were repeated in triplicate.

Protein-protein interaction (PPI) network construction analysis
STRING and Cytoscape were used to establish a PPI network, and proteins with degree > 1 were selected. The network analyzers CentiScape and MCODE of Cytoscape software were used to analyze the topology property of the network. Genes with a degree of connectivity > 15 were defined as hub genes. MCODEs were extracted when the Node score cut-off was 0.2 and K-core was 2.

Identification of DEGs from parental and IshikawaPR cells
We performed microarray analysis to identify differentially expressed genes between parental and IshikawaPR cells (Fig. 1). Volcano plots displayed the distribution of the 27,831 expressed genes (Fig. 1b). Using |LogFC| > 4 and p-value < 0.05 as cut-off criteria, 821 DEGs were extracted, including 453 upregulated and 368 downregulated genes (Fig. 1c, Additional file 1: Tables S3, Additional file 3: Table S2). In total, 760 genes were coding RNA, and 61 genes were noncoding RNA, including 30 long noncoding RNA (LncRNA), 19 noncoding RNA (NR) and 12 uncharacterized RNA. Using the Morpheus website, we developed a clustering heatmap of the DEGs (Fig. 1c). Interestingly, based on the heatmap of genes involved in lipid metabolism and biosynthetic processes, we found that lipid metabolism process-related genes, including fatty-acid biosynthetic processes, lipid translocation, regulation of lipid metabolic processes and lipid transport, were mostly upregulated, lipid biosynthetic process-related genes were downregulated in IshikawaPR cell (Fig. 2, Additional file 4: Table S4). RNA expression level changes in lipid metabolism-related genes in Ishi-kawaPR cells indicated that reduced lipid biosynthesis and increased metabolism may participate in the mechanism of progesterone resistance.

Gene ontology analysis of DEGs
To gain more biological insight, we performed gene oncology (GO) enrichment analysis using multiple online databases, including PANTHER, DAVID, and Metascape. The DEGs were classified into three functional groups: biological process (BP), molecular function (MF) and cellular component (CC). The most enriched BP functions were cellular processes (311 genes) and metabolic processes (200 genes). For MF, binding (210 genes) and catalytic activity (160 genes) were the most enriched. In the clusters of CC, cell part (170 genes) and organelle (106 genes) genes were the most enriched (Fig. 3a). As shown in Table 1 and Fig. 3b, in the biological process group, upregulated genes were mainly enriched in anterior/ posterior pattern specification, response to retinoic acid, extracellular matrix organization, antigen processing and presentation of peptide antigens via MHC class I, downregulated genes were mainly enriched in palate development, bicellular tight junction assembly, extracellular matrix organization, and angiogenesis. In the molecular function group, upregulated genes were mainly enriched in sequence-specific DNA binding. In the cellular component group, upregulated genes mainly enriched in extracellular space, extracellular region and the cell surface, downregulated genes were mainly enriched in lateral plasma membrane, ruffle membrane and bicellular tight junction. These results showed that DEGs were mainly enriched in extracellular environment, cell junction, and immune system.

Encoded protein class and KEGG pathway enrichment analysis of DEGs
To further comprehend the proteins function of those DEGs, we analyzed the protein classes of DEGs (Fig. 4a). DEGs encoded proteins were mainly distributed among transcription factors, hydrolases, enzyme modulators. In addition, the number of proteins associated with energy and lipid metabolism was also high, including fatty acid desaturase, lipase, lipoxygenase, regulator of metabolic enzymes, etc. This result was consistent with the GO analysis and DEGs identification analysis, suggesting that lipid metabolism and the immune system may play an important role in the mechanism of progesterone resistance. KEGG pathway enrichment of the 821 DEGs was conducted using the Panther database. Twenty-nine pathways were enriched with a criterion of p < 0.05 (Table 2, Fig. 4b). Complement and coagulation cascades, signaling pathways regulating pluripotency of stem cells and Staphylococcus aureus infection pathway were highly enriched in upregulated genes. Transcriptional misregulation in cancer, Rap1 signaling pathway and cancer related pathways were the most obvious pathways in downregulated genes. Among the 29 pathways, 11 pathways were involved in the immune system. Five pathways were related to extracellular environment and epithelia mesenchymal transition (EMT). This analysis suggested that extracellular environment, EMT and immune system may be critical in the development of progesterone resistance.
To determine the clinical significance of 7 target genes, Kaplan-Meier survival statistics were generated for a large cohort of endometrial cancer. In total, data from 542 endometrial cancer patients were interrogated and hazard ratios (HR) and p-values for statistical significance were determined. The data are summarized in Table 3. Interestingly, low expression of the 7 target genes correlated with poor prognosis especially ANO1, SOX17, CGNL1 and SH3YL1 (Fig. 5b, Additional file 7: Figure S1). These 7 target genes should be further studied to explore their association with PGR expression, which might expose the mechanism of progesterone resistance as well as downregulation of PR during the development of endometrial cancer.

Protein-protein interaction network (PPI) and modular analysis
The protein-protein interaction (PPI) network was generated using Metascape and visualized with Cytoscape 3.5.1. Proteins with degree > 1 were selected. In total, 352 nodes (42.9% of all 821 DEGs) and 592 PPI relationships were obtained (Fig. 7). Eight genes with a degree of connectivity > 15 were defined as hub genes for progesterone resistance (Table 4, Fig. 7a). According to the degree rank, the eight hub genes included HSPA1A, EEF1A2, AR, POU5F1, C3, SYK, LPAR1 and NMU. The 8 hub genes interact directly with 112 DEGs. As the most intensive hub gene, HSPA1A interacts with 15 upregulated and 15 downregulated genes. In addition, these hub genes could interact with each other. AR could interact with three hub genes (HSPA1A, EEF1A2, SYK), HSPA1A, EEF1A2, C3, LPAR1 and NMU could interact with other 2 hub genes. These results suggest that these hub genes might play an important role in progesterone resistance and should be further studied.
Ten densely connected regions in the networks were identified using the MCODE clustering algorithm by Metascape [20]. Four most significant MCODEs were extracted when the node score cut-off was 0.2 and K-core was 2 (Fig. 7b). GO enrichment analysis of the 4 MCODEs-related genes showed that these genes were mainly associated with the following biological processes terms: antigen processing and presentation, immune response, G-protein coupled receptor signaling pathway and chemotaxis. Regarding molecular function terms, these genes are mainly enriched in steroid hydroxylase activity, heme binding, oxidoreductase activity, aromatase activity and peptide antigen binding. Regarding cellular component terms, the top terms were MHC class I protein complex, intracellular and integral component of plasma membrane. The top 15 GO enrichment terms are presented in Table 5. KEGG pathway enrichment results showed that these genes were significantly enriched in antigen processing  and presentation, phagosome, graft-versus-host disease, allograft rejection, and some microbials infection (Table 6). These results were consistent with the GO and KEGG pathway enrichment analysis of DEGs, suggesting that the immune system and lipid metabolism play an extremely important role in progesterone resistance and should be further studied.

Discussion
Endometrial carcinoma (EC) is the most common gynecologic malignancy. De novo or acquired progesterone resistance is a major clinical problem. Our study is the first comprehensive investigation to explore the mechanism of progesterone resistance by microarray analysis. We developed an endometrial cancer cell subline that exhibits stabilized resistance to MPA and performed microarray analysis to identify DEGs. In total, 821 DEGs, 7 PGR coexpressed genes, 8 hub genes and 4 MCODEs were obtained. Using GO and KEGG pathway enrichment analysis, protein class analysis, PGR coexpression identification and validation, PPI and modular analysis of DEGs, we further explored the following possible mechanisms involved in progesterone resistance.

Identification of predictive biomarkers for progesterone treatment response
Previous research suggested several prognostic markers for response to progesterone treatment, including progesterone receptor (PR) [21], survivin [22], Fas/FasL [23], Nrf2 and AKR1C1 [24], DUSP6 [25] and DHCR24 [26]. These findings were all based on exiting research and lack a systematic view of the molecular network of progesterone resistance.
Our study is based on microarray analysis of 27,831 expressed genes, representing a more comprehensive and integrated exploration of the internal mechanism. Predictive biomarkers could identify patients who would benefit from the treatment. In total, 821 DEGs were identified from microarray analysis, including 760 coding RNA, indicating that these gene coding proteins might serve as protein biomarkers. By detecting these predictive biomarkers before treatment, patients who will be insensitive to progesterone can be screened, and can be recommended for other treatments, avoiding disease progression. In addition, the development of targeted drugs against drug resistance will increase the sensitivity to progesterone and help to provide benefit from conservative treatment. Additional validation studies are needed to verify the biomarkers value and optimal filter conditions.

Potential mechanistic hypotheses of PR dysregulation
The presence of the progesterone receptor (PR) is the precondition for progesterone response and PR is a predictive marker for response of progesterone [21]. Progesterone binds to its receptor PR-A and PR-B, subsequently inhibiting tumor growth and promoting tumor apoptosis by regulating downstream genes. Constant stimulation of progesterone reduced the expression of PGR [10] and promoted the development of drug resistance. Thus, downregulation of PR especially PRB must be involved in progesterone resistance. However, the molecular mechanism of PGR dysfunction remains unclear.
A study defined the genome-wide PR cistrome in the murine uterus using Chip-seq to identify novel PR target gene in circadian rhythm reveled that Sox17 is a novel mediator of progesterone signaling in normal endometrium [27]. After acute progesterone treatment, Sox17 is upregulated at both the mRNA and protein levels. In our study, PGR coexpression analysis showed that Sox17 was negatively correlated with PGR and downregulated in IshikawaPR cell. These results suggest that when resistance occurred, Sox17 was downregulated rather than upregulated in normal endometrium treated with progesterone. Kaplan-Meier survival statistics revealed that Table 3 Distinguishing and validation of 7 PGR co-expressed genes among DEGs PGR co-expression showed correlation analysis among target genes and PGR using Cbioportal website. The third column shows the result of microarray analysis between Ishikawa and IshikawaPR cell. Negative means they were down-regulated in IshikawaPR cell. The fourth column is the PGR correlation analysis of target genes using GSE17025 dataset. Survival curve was performed using Kaplan-Meier plotter. FC, fold change its higher expression in endometrial cancer patients is associated with better survival. Targeting Sox17 may represent a new approach to reverse resistance. In a recent study, Huizhe Wu [28] reveals that the calcium-activated chloride channel Ano1 promotes cell proliferation in ER-positive, PR-positive, and HER2-negative breast cancer cell, but inhibits cell growth in ER-negative, PR-negative, and HER2-negative breast cancer cells. These results suggest that Ano1 may differentially regulate cell proliferation depending on the status of steroid receptors. In breast cancer patients treated with tamoxifen, overexpression of the selective estrogen-receptor modulator, Ano1 is associated with good prognosis in PR-positive patients [29]. Similar to our study, Ano1 was . Font size of gene symbol determined by degree of connectivity. The width of the edge was determined by the combined score of the PPI relationship. Hub genes were HSPA1A, EEF1A2, AR, POU5F1, C3, SYK, LPAR1, NMU. b 4 MCODEs identification in the PPI network with MCODE score ≥ 4. The label size was determined according to the degree of connectivity. Yellow represented seed nodes. Red nodes were up-regulated genes and blue nodes were down-regulated genes in IshikawaPR cell downregulated in IshikawaPR cell, which is a PR-negative cell. Positive correlation was noted between Ano1 and PGR. High Ano1 expression is essential for survival and its expression was altered in papillary serous tumor and endometrioid tumor compared with healthy endometrium. Strategies to activate Ano1 are needed to promote the response to progesterone.
PR must first recruit coregulators with intrinsic histone and DNA-modifying activities and then regulate target gene. Many coregulators of hormonal receptor have been identified and Src homology 3 domain containing, Ysc84like 1 (SH3YL1) is a novel coregulator that bind to the polyproline domain of androgen receptor (AR) [30] in prostate cancer cell. In our study, SH3YL1 was downregulated in resistant cell and positively correlated with PGR. SH3YL1 expression was considerably in endometrial cancer tissues compared with inactive endometrium. Whether SH3YL1 is also the coregulator of PR remains to be investigated.
Other PGR coexpressed genes and their encoded proteins were all potential PR targets and should be further studied. Nrf2 (Nuclear factor erythroid related factor 2), overexpression of which results in progesterone resistance in endometrial cancer cells [24], binds to the SLC40A1 promoter and transcriptionally suppressed SLC40A1 expression [31]. SLC40A1 subsequently promotes Nrf2 expression. In addition, SLC40A1 is associated with cisplatin resistance in ovarian cancer. Consistent with our results, Nrf2 was upregulated and SLC40A1 was downregulated in IshikawaPR cells. Interestingly, SLC40A1 was positively correlated with PGR and exhibited reduced expression in endometrial cancer tissues than normal endometrium, indicating a potential relationship between SLC40A1 and PGR dysfunction as well as a mechanism of endometrial cancer proliferation.
Furthermore, NPAS3, CCDC146, RUNDC3B, SEMA3D, CRISPLD1, CGNL1 and CDS1 were all negatively correlated with PGR and have not been thoroughly studied. Thus, further understanding their function within progesterone resistance is required.

Immune system and inflammatory response may be associated with progesterone resistance
This study first found that the immune system and inflammatory response are potentially associated with the development of progesterone resistance. Immunomodulation in endometrial cancer and chemoresistance have been studied for years but have never been studied in the context of progesterone resistance in endometrial cancer. The immune system has pro-and antitumorigenic functions in the endometrium. By forming an integral mucosal immune system, normal endometrial epithelial cells could secrete defensin [32], which exhibits immediate antimicrobial function and the ability to active the adaptive immune system by attracting T cells and dendritic cells (DC) [33]. Endometrial epithelia cells are antigen-presenting cells that express major histocompatibility complex (MHC) and non-classical MHC class I molecules and human leukocyte antigen G (HLA-G) is down regulated in endometrial cancer [34]. In our study, HLA-G was downregulated in IshikawaPR cell and a MCODEs gene. HLA-G maybe offer protective function to avoid NK cell lysis. Hormones influence the inflammatory environment. Estrogen increases the expression of inflammatory cytokines released from murine and uterine cells but progesterone decreases their production [35]. A previous study demonstrated that the inflammatory responses is increased in mouse uterus lacking PGR, and infiltration of leukocytes and extensive tissue remodeling also increase [36]. Alterations in steroid receptors expression leads to an increase in estrogen to progesterone signaling, subsequently increasing inflammatory mediators and promoting tumor growth [37]. Many inflammatory response and immune system related GO terms were significantly enriched in our study. Within the hub genes, six of the top 20 GO enrichment terms of MCODEs genes were involved in immune response and immune systems, especially antigen process and presentation of peptide antigens. Genes involved in these terms included TAP1 (ATP-binding cassette) and TAP2, which are transporters associated with antigen processing and are major contributors to the multidrug resistance phenotype [38]. In addition, CXCL1 and CXCL8 are chemokines that are elevated in endometrial adenocarcinoma [39,40]. However in progesterone resistant cells, these genes are downregulated, requiring further study.
Another potential mechanism of progesterone resistance involves the G protein-coupled estrogen receptor 1 (GPER1), which mediates immunomodulation. GPER1 was upregulated in progesterone resistant cell in our study and is also a MCODEs genes. Given that estrogen anti-inflammatory activity is mediated by estrogen receptor ERα36 and GPR30/GPER1 [41], the increase in GPER1 expression in progesterone-resistant cells may hasten the inflammatory response and promote the proliferation effect by estrogen. Our study provides supportive evidence for the novel theory that the immune system and inflammatory response play vital role in progesterone resistance of endometrial cancer. These findings may provide new strategies to block the immune response to reverse the progesterone resistance in the future.

The activation of lipid metabolism is vital in progesterone resistance
In cancer cells, previous studies declare that metabolism alterations including aerobic glycolysis and de novo lipid biosynthesis promote the transformation of cells from a noncancerous to cancerous state [42]. Given that fatty acids provide twice as much ATP as glucose and are preferred nutrients for storage, they serve as an energy source when the requirement for ATP production increases in cancer cells [43]. Several studies have demonstrated the vital role of fatty acid metabolism in endometrial cancer progression. We previously reported that SREBP1 promoted fatty acid metabolism by regulating the transcription of FASN and played an important role in the tumorigenesis of endometrial cancer [44]. In addition, in endometrial cancer cell, SREBP1/FASN was suppressed during the proliferation suppression and apoptosis induced by progesterone [45]. In this study, the expression of genes involved in fatty acid biosynthetic processes, lipid metabolic processes, lipid translocation and regulation of lipid metabolic process related genes were all increased. However, genes involved in the lipid biosynthetic process were decreased. Consistent with KEGG pathway analysis, elevation of fatty acid synthesis and lipid metabolism were extremely important in progesterone resistance.

Conclusions
We offer a comprehensive analysis of gene expression profile in progesterone resistant endometrial cancer cell and raise some novel, potential mechanistic hypotheses of progesterone resistance. We identified 821 DEGs, 7 PGR coexpressed genes, 8 hub genes and 4 MCODEs that are significantly enriched in immune system, inflammatory response, lipid metabolism and extracellular environment. These findings provide comprehensive information for further study how to promote the response of progesterone treatment in endometrial carcinoma.