- Primary research
- Open Access
ESR1 as a recurrence-related gene in intrahepatic cholangiocarcinoma: a weighted gene coexpression network analysis
Cancer Cell International volume 21, Article number: 225 (2021)
Intrahepatic cholangiocarcinoma (iCCA) is the second most common malignant hepatic tumor and has a high postoperative recurrence rate and a poor prognosis. The key roles of most tumor recurrence-associated molecules in iCCA remain unclear. This study aimed to explore hub genes related to the postsurgical recurrence of iCCA.
Differentially expressed genes (DEGs) between iCCA samples and normal liver samples were screened from The Cancer Genome Atlas (TCGA) database and used to construct a weighted gene coexpression network. Module-trait correlations were calculated to identify the key module related to recurrence in iCCA patients. Genes in the key module were subjected to functional enrichment analysis, and candidate hub genes were filtered through coexpression and protein–protein interaction (PPI) network analysis. Validation studies were conducted to detect the “real” hub gene. Furthermore, the biological functions and the underlying mechanism of the real hub gene in iCCA tumorigenesis and progression were determined via in vitro experiments.
A total of 1019 DEGs were filtered and used to construct four coexpression modules. The red module, which showed the highest correlations with the recurrence status, family history, and day to death of patients, was identified as the key module. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses demonstrated that genes in the red module were enriched in genes and pathways related to tumorigenesis and tumor progression. We performed validation studies and identified estrogen receptor 1 (ESR1), which significantly impacted the prognosis of iCCA patients, as the real hub gene related to the recurrence of iCCA. The in vitro experiments demonstrated that ESR1 overexpression significantly suppressed cell proliferation, migration, and invasion, whereas ESR1 knockdown elicited opposite effects. Further investigation into the mechanism demonstrated that ESR1 acts as a tumor suppressor by inhibiting the JAK/STAT3 signaling pathway.
ESR1 was identified as the real hub gene related to the recurrence of iCCA that plays a critical tumor suppressor role in iCCA progression. ESR1 significantly impacts the prognosis of iCCA patients and markedly suppresses cholangiocarcinoma cell proliferation, migration and invasion by inhibiting JAK/STAT3 signaling pathway.
Intrahepatic cholangiocarcinoma (iCCA) is the second most common malignant hepatic tumor after hepatocellular carcinoma (HCC) [1,2,3]. iCCA is a devastating disease with a poor prognosis and a high postoperative recurrence rate of 40–80% [3, 4]. Fewer than one-third of iCCA patients who undergo curative-intent surgery survive beyond 5 years after surgery [4, 5]. The age-specific mortality rate of this disease has increased in almost all countries across all continents, albeit at different rates . The rising incidence and high mortality of iCCA have aroused growing worldwide concerns, especially over the last two decades . The significant increase in the mortality rate of this primary hepatobiliary cancer coincides with the rapidly growing interest of clinicians and investigators [3, 7]. It is therefore an urgent task to clarify the mechanisms underlying the rapid progression of iCCA and identify novel molecular biomarkers for the recurrence and prognosis of this disease. Risk factors contributing to the prognosis of iCCA include cirrhosis, chronic viral hepatitis, excess alcohol use, diabetes, and obesity. The similarities between the risk factors of iCCA and HCC suggest that common pathobiological pathways may exist in all primary liver parenchymal tumors . Numerous genetic alterations occur during iCCA development and progression. These alterations include mutational activation involving IDH1/2, KRAS, BRAS, and EGFR; mutational inactivation involving ARID1A, BAP1, and PBRM1; and changes in tyrosine kinase (TK) fusion proteins, including FGFR2 and ROS [2, 9, 10]. Tumor development and progression are complex and involve progressive and multifactorial processes . However, the molecular pathogenesis of iCCA remains unclear.
With the discovery and development of high-throughput data analysis, screening for key information has become the basis for subsequent research. Weighted gene coexpression network analysis (WGCNA) is an effective bioinformatics strategy commonly used to explore the complex relationships between genes and related clinical traits. This method has been increasingly used in bioinformatics to analyze gene expression microarray profiles [12,13,14]. Genes that are highly coexpressed with other genes are connected in networks and then grouped into modules. Each module exhibits highly connected network regions, and genes in the same module may share similar biological regulatory functions. The relationship between modules and clinical traits can also be analyzed. Hub genes in hub modules are potential biomarkers for prognosis or therapy and can thus be selected for further validation.
In this study, we used WGCNA for the first time to analyze the clinical traits and mRNA expression profiles of iCCA patients obtained from The Cancer Genome Atlas (TCGA) database in an attempt to identify hub genes associated with iCCA recurrence after surgery. Our research determined that ESR1 plays a critical tumor suppressor role in iCCA. iCCA patients with low ESR1 expression had a poorer prognosis than those with high ESR1 expression. These findings could be beneficial for the assessment of the malignant potential of iCCA and provide therapeutic methods for this cancer.
Materials and methods
Data and tissue collections
The mRNA expression profiles and related clinical traits of iCCA patients were obtained from the TCGA database (https://cancergenome.nih.gov/). The TCGA–CHOL project consists of seven hilar/perihilar cholangiocarcinoma samples, two distal cholangiocarcinoma samples, 30 iCCA samples, and eight adjacent normal samples. The RNA-Seq data of 30 iCCA samples and eight corresponding adjacent normal samples were selected for analysis. Another 30 paired iCCA tumor tissues and normal liver tissues were collected from Eastern Hepatobiliary Surgery Hospital as the validation cohort. Informed consent was obtained from all the patients, and all procedures performed were in accordance with the Helsinki Declaration of 1975.
Screening for differentially expressed genes (DEGs)
DEGs between iCCA samples and adjacent normal samples were screened using the “DESeq2” R package, in which the thresholds were set to false discovery rate (FDR) < 0.01 and |log2FoldChange|> 1.0.
The quality of the RNA-Seq data was evaluated before the construction of the weighted coexpression network. Samples were deemed defective if they were distant from the majority in clustering by average linkage (defective sample = 0.75). The remaining samples were selected for follow-up analysis. The “WGCNA” package in R was used to construct a scale-free coexpression network for the DEGs that met the screening criteria. We considered that selection of the appropriate soft-thresholding power for a scale-free coexpression network is crucial and separately evaluating thresholding powers from 1 to 10 to identify the most appropriate value. Adjacencies among all DEGs were calculated using the selected power function. The adjacency matrix was subsequently transformed into a topological overlap matrix (TOM) to calculate dissimilarity (1-TOM). Modules were identified through a dynamic tree cut method, and the minimum size of each module was set to 30 genes. MEDissTres was set as 0.24 to merge similar modules.
Identification of significant clinical modules
Gene significance (GS), module significance (MS) and the module eigengene (ME) were assessed. GS was defined as the log10 transformation of the P-value in the linear regression between gene expression and clinical traits (GS = lgP). MS was defined as the mean GS for all the genes in a module. The ME was defined as the first principal component in the principal component analysis for a given gene module. The module with the maximum absolute MS among all modules is generally regarded as related to a clinical trait. In addition, the correlation between the MEs and clinical traits was calculated. The module with the highest absolute MS was considered the key module related to a given clinical trait and was selected for further analysis.
Functional enrichment analysis of the key modules
The Database for Annotation, Visualization, and Integrated Discovery (DAVID, https://david.ncifcrf.gov/summary.jsp) was used for GO enrichment analysis. The gene list for the key module was uploaded to DAVID, and Gene Ontology (GO) enrichment results and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were obtained. A P-value < 0.05 was deemed statistically significant.
Identification of candidate hub genes
The hub genes in the hub module were initially selected using the following criteria: absolute value of GS for an interesting clinical trait > 0.2 and module membership in the hub module > 0.8. The genes screened with these criteria comprised the Color_Hub gene list. All genes from the hub module (P-value of GS < 0.05, P-value of MM < 0.05) were simultaneously uploaded to the Search Tool for the Retrieval of Interacting Genes database (STRING). A protein–protein interaction (PPI) network was constructed based on a minimum required interaction confidence score > 0.4. The result of the PPI network analysis was downloaded and visualized again using Cytoscape software. The hub genes in the PPI network were screened using the cytoHubba  application in Cytoscape software, which was designed to identify key nodes in networks. The key genes were screened using cytoHubba to generate the PPI_Hub gene list. Genes common to the Color_Hub and PPI_Hub gene lists were regarded as candidate hub genes for further validation and analysis.
Validation of real hub genes
Thirty pairs of iCCA and adjacent nontumor tissues were obtained from the Eastern Hepatobiliary Surgery Hospital (Shanghai, China). The mRNA expression levels of the screened candidate hub genes were validated by qRT-PCR. In addition, differences in the protein levels of the candidate genes between iCCA specimens and normal liver specimens were investigated by immunohistochemistry (IHC). Finally, the candidate hub genes that significantly influenced the disease-free survival (DFS) of patients were selected as the “real” hub genes.
Protein was extracted from cells with RIPA lysis buffer supplemented with a protease inhibitor. Protein extracts (30 μg) were subjected to 4–20% gradient SDS/PAGE and blotted onto a nitrocellulose (NC) membrane. The membrane was then blocked with 5% nonfat milk in TBS (Sigma, T8793), incubated with the primary antibody overnight at 4 °C and incubated with the secondary antibody; then, the membrane was washed with TBST three times, and finally visualized using an enhanced chemiluminescence (ECL) detection kit (Thermo Fisher Scientific, USA). Quantification was performed with ImageJ (Version 2.0.0).
For IHC analysis, tissue microarray (TMA) sections were incubated with anti-ESR1 antibody (1:400 dilution). IHC staining was scored by two independent pathologists who were blinded to the clinical characteristics of the patients. The scoring system was based on the intensity and extent of staining. Staining intensity was classified as 0 (negative), 1 (weak), 2 (moderate), or 3 (strong); staining extent was dependent on the percentage of positive cells (among 200 assessed cells) and was classified as 0 (< 5%), 1 (5–25%), 2 (26–50%), 3 (51–75%), or 4 (> 75%). According to the staining intensity and staining extent scores, the IHC result was classified as 0–1, negative (−); 2–4, weakly positive ( +); 5–8, moderately positive (+ +), and 9–12, strongly positive (+ + +).
Cells were cultured at 37 °C with a 5% CO2 atmosphere. HiBEpiC, RBE, HCCC-9810 and FRH0201 cells were cultured in DMEM (Gibco, USA) supplemented with 10% FBS (Gibco), and QBC939 cells were cultured in RPMI 1640 (Gibco) supplemented with 10% FBS (Gibco). At 80–90% confluence, the cells were trypsinized and replated at 5 × 105 cells/75 cm2.
We generated an ESR1-overexpressing HCCC-9810 cell line using a lentiviral system and an ESR1-knockdown RBE cell line using siRNAs targeting ESR1. The corresponding negative control cell lines were also generated using an empty lentivirus or a nontargeting control siRNA. Our cell transfection experiments were performed using Lipofectamine 3000 (Invitrogen, USA) according to the manufacturer’s instructions. The efficiency of cell transfection was determined by Western blotting after 48 h.
Cell proliferation assay
Cell proliferation ability was investigated using the Cell Counting Kit-8 (CCK-8; Beyotime Institute of Biotechnology, Shanghai, China) assay. Approximately 5000 cells were plated in 96-well plates for 12 h. Then, the supernatant was removed, and CCK-8 solution (Thermo Fisher Scientific) was added to each well and further incubated for 24, 48, and 72 h at 37 °C. Finally, the absorbance of each well at 450 nm was detected by a spectrophotometer (Thermo Fisher Scientific, Waltham, USA).
Colony formation assay
Cells at a density of 200 cells/well were plated in six-well plates. After culturing for 14 days, the generated colonies were fixed with methanol for 15 min, stained with a hematoxylin solution and photographed using a microscope (Olympus, Japan). The assays were conducted in triplicate.
Cell migration and invasion assay
Cell migration and invasion assays were performed using Transwell chambers (Corning, USA). Approximately 1 × 105 transfected cells were inoculated in the upper chamber with serum-free medium (with a permeable membrane) and cultured for 24 h. Medium with 10% FBS was added into the lower chamber as a chemical inducer. After 12 h of incubation, cells that migrated to the bottom of the membrane were fixed with 4% formaldehyde for 10 min. After staining and fixation with 0.1% crystal violet, the cells were observed and counted under a microscope.
All statistical analyses were performed using R 3.6.1. For quantitative values according to the normal distribution, Student’s t-test was used for comparison. For nonnormally distributed quantitative data, the Mann–Whitney U test was utilized for comparison. DFS was assessed with the Kaplan–Meier method, and the differences between the groups were compared by the log-rank test. P < 0.05 was considered statistically significant.
Screening of DEGs
RNA-Seq data from the TCGA–CHOL dataset were downloaded from the TCGA database. Seven hilar/perihilar cholangiocarcinoma tissues and two distal cholangiocarcinoma tissues were excluded in accordance with the postoperative pathological results. The RNA-Seq data for 30 iCCA tissues and eight adjacent normal tissues were retained for subsequent analysis. The clinicopathological characteristics of the aforementioned 30 iCCA patients are listed in Table 1. After data preprocessing, 1019 DEGs (184 upregulated and 835 downregulated) were screened using the threshold criteria of FDR < 0.01 and |log2FoldChange|> 1.0 (Additional file 1: Table S1).
Construction of the weighted gene coexpression network
All DEGs were used to construct a weighted gene coexpression network. One sample was excluded because it was identified as defective according to the given criteria. The 29 remaining iCCA samples were clustered using the “flashClust” R package. The results of the clustering and relative trait analyses are illustrated in Fig. 1a. Then, thresholding powers from 1 to 10 were tested separately to identify the most appropriate soft-thresholding power. As shown in Fig. 1b, a power value of 8 was ultimately selected for constructing a scale-free network with a scale-free R2 of 0.87 and a mean connectivity of 10.12 (Additional file 2: Figure S1A). Seven coexpression modules were constructed after excluding the gray module. The MEDissThres value was set as 0.24 to merge similar modules, leaving four modules (Fig. 1c, Additional file 2: Figure S1B). The network heatmap of the four modules is shown in Fig. 1d, and the results showed high correlations between the intramodular genes. The black, blue, brown and red modules included 406, 171, 132 and 231 genes, respectively (Additional file 3: Table S2).
Identification of the key module
The interaction relationship between the four modules was calculated. All these modules independently validated each other, indicating a high level of independence between the four modules (Additional file 2: Figure S1C). In addition, the correlation between the module eigengenes and clinical traits was calculated (Fig. 1e). The red module was selected as the key module because it had the three highest correlations with the day to death, family history and recurrence status. We were particularly interested in the relationships between the recurrence status of the patients and the expression of hub genes in the red module, given its potential significance in the prognosis of iCCA patients. Furthermore, univariable COX regression for DFS was performed, and the results revealed that ESR1 expression was an independent prognostic factor for better DFS of iCCA patients in both cohorts (Table 2).
Functional enrichment analysis of genes in the red module
All genes in the red module were uploaded to DAVID for functional enrichment analysis. All significantly enriched GO terms (P < 0.05) are plotted in Fig. 2a. The GO enrichment analysis showed that the genes in the red module were significantly enriched in terms of cofactor binding, GTPase binding, negative regulation of protein phosphorylation, second messenger-mediated signaling, and epithelial cell proliferation. Most of these terms were related to tumorigenesis and tumor microenvironment regulation. Then, KEGG pathway enrichment of genes in the red module was performed (Fig. 2b). The genes in the red module were significantly enriched in pathways that were closely related to malignancy. These pathways included the Rap1 signaling pathway, cGMP-PKG signaling pathway, and N-glycan biosynthesis. In summary, these results indicate that genes in the red module are closely associated with tumor development and disease progression.
Identification of candidate hub genes in the red module
The candidate hub genes in the red module were initially selected based on the following criteria: absolute value of GS for recurrence > 0.2 and MM > 0.8 (Fig. 2c). Forty genes were selected and comprised the Red_Hub gene list (Fig. 2e, Additional file 4: Table S3). A PPI network was constructed using the appropriate cutoff for the defective samples of confidence > 0.4 (Additional file 2: Figure S1D). The key genes in the PPI network were screened by employing the cytoHubba algorithm. The 40 genes were filtered and formed the PPI_Hub gene list (Fig. 2d). The genes common to the two gene lists were selected as candidate hub genes for predicting recurrence. Finally, 11 genes (Fig. 2e), namely, ESR1, A2M, TEK, KDR, S1PR1, HGF, RAPGEF4, STAB2, EPAS1, TAB2, and VDR, were selected. Further validation experiments were conducted for these candidate hub genes.
Validation of real hub genes in the red module
Three approaches were used to screen the real hub genes. First, qRT-PCR was adopted to measure the relative transcript abundance of the 11 candidate hub genes between 30 iCCA tumor tissues and adjacent nontumor tissues from patients in the validation cohort. Four genes were excluded from this step because the difference in the relative mRNA expression levels between the tumor tissues and adjacent nontumor tissues was not statistically significant (P-value > 0.05). Seven candidate hub genes (ESR1, A2M, RAPGEF4, S1PR1, HGF, STAB2, and VDR) were retained for the following validation (Additional file 2: Figure S1E). The IHC results from the Human Protein Atlas database demonstrated that only three genes (ESR1, A2M, and RAPGEF4) showed expression differences at the protein level between the iCCA specimens and normal liver specimens (Fig. 3a, Additional file 2: Figure S1E). Then, the Table-S4 values of ESR1, A2M, and RAPGEF4 were determined using receiver operating characteristic (ROC) analysis and Youden index calculation. ESR1 was filtered and deemed the real hub gene that has a significant impact on the DFS of the patients (Fig. 3b, Additional file 5: Table S4).
The protein expression level of ESR1 in cholangiocarcinoma cell lines
To detect the protein expression level of ESR1 in cholangiocarcinoma cell lines (QBC939, RBE, HCCC-9810, and FRH0201) and a normal biliary epithelial cell line (HiBEpiC), a Western blot assay was conducted. The results showed that the ESR1 expression level in the HCCC-9810 cell line was the lowest and that the RBE cell line was the highest among the four cholangiocarcinoma cell lines (Fig. 3c). Therefore, we conducted gain-of-function experiments in the HCCC-9810 cell line and loss-of-function experiments in the RBE cell line. As shown in Fig. 3d, the protein expression level of ESR1 in the Si-ESR1 group was significantly lower than that in the NC-ESR1 group in the RBE cell line (P < 0.0001), and the protein expression level of ESR1 in the ESR1 overexpression (Ov-ESR1) group was significantly higher than that in the NC-ESR1 group in the HCCC-9810 cell line (P < 0.0001).
The effect of ESR1 on iCCA cell proliferation, migration and invasion
To determine the effect of ESR1 on iCCA cell proliferation, CCK-8 and colony forming assays were performed. ESR1 overexpression significantly inhibited cell proliferation, while ESR1 knockdown promoted cell proliferation (Fig. 3e, g). To further determine the effect of ESR1 on cell migration and invasion, a Transwell assay was performed. ESR1 overexpression markedly inhibited the migration and invasion abilities of HCCC-9810 cells, and knockdown of ESR1 promoted the migration and invasion abilities of RBE cells (P < 0.001, Fig. 3f). All these findings demonstrate that ESR1 expression inhibits the proliferation, migration and invasion of iCCA cells.
ESR1 functions as a tumor suppressor by inhibiting the JAK/STAT3 signaling pathway
To further investigate the potential antitumor mechanism of ESR1, we explored whether ESR1 inhibits the proliferation, migration and invasion of iCCA cells through the JAK/STAT3 signaling pathway. As the experiments illustrated, ESR1 overexpression resulted in a significant reduction in the expression of p-JAK1 and p-STAT3 in the HCCC-9810 cell line, while the expression levels of JAK1 and STAT3 did not change significantly. Our results illustrated that ESR1 may act as a tumor suppressor by inhibiting the JAK/STAT3 signaling pathway.
iCCA is the second most common primary liver cancer and has a high likelihood of postoperative recurrence. The incidence of iCCA is on the rise worldwide . High-throughput genomic analysis is a promising tool that has considerable clinical applications in medical oncology [17, 18]. The initiation and progression of iCCA are complex, given that they involve various genetic changes . Although many genes have been selected for the diagnosis and treatment of iCCA, novel potential biomarkers are still required to improve the understanding of tumor progression, especially tumor recurrence after surgery. In this study, WGCNA was used to screen candidate biomarkers associated with the recurrence of iCCA.
First, based on the TCGA-CHOL dataset, 1019 DEGs were filtered out after dataset preprocessing and were used to construct a weighted gene coexpression network. Although our threshold for filtering DEGs was higher than that in previous similar studies [13, 17], it proved to be effective. After network construction, module identification, and assessment of the relations between the modules and clinical traits, we found that the red module was associated with tumor recurrence. Finally, after a series of hub gene identification and validation studies, ESR1 was selected as the real hub gene in our analysis and was significantly associated with recurrence.
ESR1 encodes an estrogen receptor. Estrogen and its receptors are essential for sexual development and reproductive function. Previous studies have shown that ESR1 is closely related to the occurrence and development of various urogenital cancers [20,21,22,23], especially breast cancer [24, 25] and endometrial cancer [26,27,28]. Recent studies have reported the important role of ESR1 in non-small-cell lung cancer (NSCLC) [29, 30] and bladder cancer . In addition, another study demonstrated that the methylation ratio of ESR1 was closely associated with old age and smoking history in patients with NSCLC .
Recently, with the rapid development of genome-wide technologies, we have gained deep insight into the molecular mechanisms of iCCA, and more potential biomarkers have been identified for clinical use. For example, Chen et al. found that TGF-β1 was an independent risk factor for early tumor recurrence in iCCA . In addition, a study found that iCCA patients with high levels of CYFRA21-1 had worse 3-year recurrence-free survival rates than those with low levels (DFS: 25.0% vs. 76.2%, P < 0.01), which suggests that CYFRA21-1 is associated with tumor recurrence . Other studies showed that low expression of SMAD4, a tumor suppressor protein, increased the rates of lymph node recurrence in iCCA . High epidermal growth factor receptor (EGFR) expression was also identified as a risk factor for poor overall survival (P = 0.0006) and recurrence (P = 0.0335) in iCCA patients .
In 2020, using data from the Gene Expression Omnibus (GEO) database, Qin et al. found that ESR1 was downregulated in iCCA patients. ESR1 and its corresponding miRNA, hsa-miR-26a-5p, might be novel biomarkers for prognosis . However, Qin et al. failed to validate their conclusions through in vitro experiments and clinical tissue samples. Moreover, ESR1 was identified through an indirect method in Qin’s study. In this study, we found that ESR1 is a hub gene in iCCA that is significantly associated with recurrence through WGCNA and PPI network analysis. Then, a series of in vitro experiments demonstrated that ESR1 negatively regulates iCCA cell proliferation, migration and invasion. Furthermore, we performed in-depth studies to understand the mechanism and illustrated that ESR1 may act as a tumor suppressor by inhibiting the JAK/STAT3 signaling pathway. Our results suggest that the roles of ESR1 in human cancers are more extensive than previously believed and deepen our understanding of ESR1 function in human cancers.
Our study has several limitations. First, ESR1 was identified as a hub gene through analysis of 29 iCCA patients from the TCGA database, and the validation cohort was also small, with only 30 patients; thus, the statistical power of the calculations in our study is low. More iCCA patients should be used to validate our conclusions in the future. Second, our study was a retrospective study, so the significance and robustness of the results and recurrence-related hub genes should be validated in prospective cohorts. Last, the clinical applicability of ESR1 as a biomarker is limited due to a lack of adequate reproducibility in daily clinical practice.
In conclusion, our research revealed that ESR1 plays a critical role as a tumor suppressor in iCCA. ESR1 significantly impacts the prognosis of iCCA patients and markedly inhibits cell proliferation, migration and invasion. ESR1 may exert a tumor suppressor function by inhibiting the JAK/STAT3 pathway. The findings of the present study may help researchers gain a better understanding of the mechanism of iCCA development and provide a novel target for the treatment of iCCA in the future.
ESR1 was identified as a recurrence-related gene of iCCA via WGCNA and PPI network analysis.
Low expression of ESR1 predicts a poor prognosis in iCCA patients.
The expression of ESR1 inhibits the proliferation, migration and invasion of iCCA cells.
ESR1 may exert a tumor suppressor function by inhibiting the JAK/STAT3 signaling pathway.
Availability of data and materials
All data generated or analyzed during this study are included in the additional files.
Differentially expressed genes
The Cancer Genome Atlas
The Kyoto Encyclopedia of Genes and Genomes
The Database for Annotation, Visualization, and Integrated Discovery
The search tool for the retrieval of interacting genes database
Gene expression profiling interactive analysis
Sodium dodecyl sulfate
Polyacrylamide gel electrophoresis
Non-small-cell lung cancer
Kelley RK, Bridgewater J, Gores GJ, Zhu AX. Systemic therapies for intrahepatic cholangiocarcinoma. J Hepatol. 2020;72(2):353–63.
Moeini A, Sia D, Bardeesy N, Mazzaferro V, Llovet J. Molecular pathogenesis and targeted therapies for intrahepatic cholangiocarcinoma. Clin Cancer Res. 2016;22(2):291–300.
Sirica A, Gores G, Groopman J, Selaru F, Strazzabosco M, Wei Wang X, Zhu A. Intrahepatic cholangiocarcinoma: continuing challenges and translational advances. Hepatology. 2019;69(4):1803–15.
Wirth TC, Vogel A. Surveillance in cholangiocellular carcinoma. Best Pract Res Clin Gastroenterol. 2016;30(6):987–99.
Mavros MN, Economopoulos KP, Alexiou VG, Pawlik TM. Treatment and prognosis for patients with intrahepatic cholangiocarcinoma: systematic review and meta-analysis. JAMA Surg. 2014;149(6):565–74.
Saleh M, Virarkar M, Bura V, Valenzuela R, Javadi S, Szklaruk J, Bhosale P. Intrahepatic cholangiocarcinoma: pathogenesis, current staging, and radiological findings. Abdominal Radiol. 2020;45:3662–80.
Mejia JC, Pasko J. Primary liver cancers: intrahepatic cholangiocarcinoma and hepatocellular carcinoma. Surg Clin North Am. 2020;100(3):535–49.
Palmer W, Patel T. Are common factors involved in the pathogenesis of primary liver cancers? A meta-analysis of risk factors for intrahepatic cholangiocarcinoma. J Hepatol. 2012;57(1):69–76.
Martínez-Jiménez F, Muiños F, Sentís I, Deu-Pons J, Reyes-Salazar I, Arnedo-Pac C, Mularoni L, Pich O, Bonet J, Kranas H, et al. A compendium of mutational cancer driver genes. Nat Rev Cancer. 2020;20:555–72.
Patel T. New insights into the molecular pathogenesis of intrahepatic cholangiocarcinoma. J Gastroenterol. 2014;49(2):165–72.
Miller KD, Fidler-Benaoudia M, Keegan TH, Hipp HS, Jemal A, Siegel RL. Cancer statistics for adolescents and young adults, 2020. CA Cancer J Clin. 2020;70:443–59.
Gu H, Yang M, Guo J, Zhang C, Lin L, Liu Y, Wei R. Identification of the biomarkers and pathological process of osteoarthritis: weighted gene co-expression network analysis. Front Physiol. 2019;10:275.
Qian G, Chen L, Wu C, Dan H, Xiao Y, Wang X. Co-expression network analysis of biomarkers for adrenocortical carcinoma. Front Genet. 2018;9:1–13.
Wang J, Zhang J. Identification of biomarkers of chromophobe renal cell carcinoma by weighted gene co-expression network analysis. Cancer Cell Int. 2018;18:206.
Chin C, Chen S, Wu H, Ho C, Ko M, Lin C. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8(Suppl 4):S11.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
Tang J, Lu M, Cui Q, Zhang D, Kong D, Liao X, Ren J, Gong Y, Wu G. Overexpression of ASPM, CDC20, and TTK confer a poorer prognosis in breast cancer identified by gene co-expression network analysis. Front Oncol. 2019;9:310.
Fabris L, Andersen JB, Fouassier L. Intrahepatic cholangiocarcinoma: A single-cell resolution unraveling the complexity of the tumor microenvironment. J Hepatol. 2020;73:1007–9.
Yu Y, Liu Q, Li W, Qu Y, Zhang Y, Liu T. Identification of a novel EHBP1-MET fusion in an intrahepatic cholangiocarcinoma responding to crizotinib. Oncologist. 2020;25:1005–8.
Brokken L, Lundberg-Giwercman Y, Rajpert De-Meyts E, Eberhard J, Ståhl O, Cohn-Cedermark G, Daugaard G, Arver S, Giwercman A. Association of polymorphisms in genes encoding hormone receptors ESR1, ESR2 and LHCGR with the risk and clinical features of testicular germ cell cancer. Mol Cell Endocrinol. 2012;351(2):279–85.
Fu C, Dong W, Wang A, Qiu G. The influence of ESR1 rs9340799 and ESR2 rs1256049 polymorphisms on prostate cancer risk. Tumour Biol. 2014;35(8):8319–28.
Kostova M, Brennen W, Lopez D, Anthony L, Wang H, Platz E, Denmeade S. PSA-alpha-2-macroglobulin complex is enzymatically active in the serum of patients with advanced prostate cancer and can degrade circulating peptide hormones. Prostate. 2018;78(11):819–29.
McIntyre M, Kantoff P, Stampfer M, Mucci L, Parslow D, Li H, Gaziano J, Abe M, Ma J. Prostate cancer risk and ESR1 TA, ESR2 CA repeat polymorphisms. Cancer Epidemiol Biomark Prev. 2007;16(11):2233–6.
Nagel A, Szade J, Iliszko M, Elzanowska J, Welnicka-Jaskiewicz M, Skokowski J, Stasilojc G, Bigda J, Sadej R, Zaczek A, et al. Clinical and biological significance of gene alteration and estrogen receptors isoforms expression in breast cancer patients. Int J Mol Sci. 2019;20(8):1881.
Zinger L, Merenbakh-Lamin K, Klein A, Elazar A, Journo S, Boldes T, Pasmanik-Chor M, Spitzer A, Rubinek T, Wolf I. Ligand-binding domain-activating mutations of ESR1 rewire cellular metabolism of breast cancer cells. Clin Cancer Res. 2019;25(9):2900–14.
Holst F, Hoivik E, Gibson W, Taylor-Weiner A, Schumacher S, Asmann Y, Grossmann P, Trovik J, Necela B, Thompson E, et al. Recurrent hormone-binding domain truncated ESR1 amplifications in primary endometrial cancers suggest their implication in hormone independent growth. Sci Rep. 2016;6:25521.
Kurz S, Thieme R, Amberg R, Groth M, Jahnke H, Pieroh P, Horn L, Kolb M, Huse K, Platzer M, et al. The anti-tumorigenic activity of A2M-A lesson from the naked mole-rat. PLoS ONE. 2017;12(12):e0189514.
O’Mara T, Glubb D, Painter J, Cheng T, Dennis J, Attia J, Holliday E, McEvoy M, Scott R, Ashton K, et al. Comprehensive genetic assessment of the ESR1 locus identifies a risk region for endometrial cancer. Endocr Relat Cancer. 2015;22(5):851–61.
Lin Q, Geng J, Ma K, Yu J, Sun J, Shen Z, Bao G, Chen Y, Zhang H, He Y, et al. RASSF1A, APC, ESR1, ABCB1 and HOXC9, but not p16INK4A, DAPK1, PTEN and MT1G genes were frequently methylated in the stage I non-small cell lung cancer in China. J Cancer Res Clin Oncol. 2009;135(12):1675–84.
Suga Y, Miyajima K, Oikawa T, Maeda J, Usuda J, Kajiwara N, Ohira T, Uchida O, Tsuboi M, Hirano T, et al. Quantitative p16 and ESR1 methylation in the peripheral blood of patients with non-small cell lung cancer. Oncol Rep. 2008;20(5):1137–42.
Ge Q, Lu M, Ju L, Qian K, Wang G, Wu C, Liu X, Xiao Y, Wang X. miR-4324-RACGAP1-STAT3-ESR1 feedback loop inhibits proliferation and metastasis of bladder cancer. Int J Cancer. 2019;144(12):3043–55.
Chen Y, Ma L, He Q, Zhang S, Zhang C, Jia W. TGF-beta1 expression is associated with invasion and metastasis of intrahepatic cholangiocarcinoma. Biol Res. 2015;48:26.
Uenishi T, Yamazaki O, Tanaka H, Takemura S, Yamamoto T, Tanaka S, Nishiguchi S, Kubo S. Serum cytokeratin 19 fragment (CYFRA21-1) as a prognostic factor in intrahepatic cholangiocarcinoma. Ann Surg Oncol. 2008;15(2):583–9.
Chuang SC, Lee KT, Tsai KB, Sheen PC, Nagai E, Mizumoto K, Tanaka M. Immunohistochemical study of DPC4 and p53 proteins in gallbladder and bile duct cancers. World J Surg. 2004;28(10):995–1000.
Yoshikawa D, Ojima H, Iwasaki M, Hiraoka N, Kosuge T, Kasai S, Hirohashi S, Shibata T. Clinicopathological and prognostic significance of EGFR, VEGF, and HER2 expression in cholangiocarcinoma. Br J Cancer. 2008;98(2):418–25.
Qin X, Song Y. Bioinformatics analysis identifies the estrogen receptor 1 (ESR1) Gene and hsa-miR-26a-5p as potential prognostic biomarkers in patients with intrahepatic cholangiocarcinoma. Med Sci Monit. 2020;26:e921815.
We would like to thank the biotrainee team for technical advice in data analysis.
This study was supported by grants from the State Key Project on Infectious Diseases (2018ZX10723204-001); the Joint Tackling Project of Emerging Frontier Technologies in Shanghai Hospitals in 2017 (SHDC12017122); the Key Breakthrough Direction of Special Disease of the Navy Medical University Investigative-oriented Department and Investigative Type Doctor Construction Project in 2017 (Kui Wang and Yong Xia); and the Establishment of Guidelines for Surgical Treatment of Hepatolithiasis (201502014).
Ethics approval and consent to participate
This study was carried out based on the ethics approval of Eastern Hepatobiliary Surgery Hospital. The clinical analysis was performed according to the principles of the Helsinki Declaration.
Consent for publication
All authors agree to publish.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The list of the 1019 DEGs in 29 iCCA samples from TCGA.
Analysis of the mean connectivity for various soft-thresholding powers (b). Mean connectivity (y-axis) is a strictly decreasing function of the power b (x-axis). (B) Clustering of MEs. MEDissTres was set as 0.24 to merge similar modules. The red and yellow modules were merged into the new red module, which was selected as the hub module in our research. (C) Heatmap of the adjacencies in the hub gene network. (D) PPI network of all genes in the red module. (E) qRT-PCR results for ESR1, A2M, S1PR1, HGF, VDR, RAPGEF4 and STAB2. (F) IHC results for A2M and RAPGEF4 in iCCA tissues and normal liver tissues based on the Human Protein Atlas database. A2M and RAPGEF4 were significantly overexpressed in normal liver tissue vs. iCCA tissue.
GS and MM in four integrated modules.
The genes in red_hub list and PPI_hub list.
Validation of real hub genes.
About this article
Cite this article
Li, F., Chen, Q., Yang, Y. et al. ESR1 as a recurrence-related gene in intrahepatic cholangiocarcinoma: a weighted gene coexpression network analysis. Cancer Cell Int 21, 225 (2021). https://doi.org/10.1186/s12935-021-01929-5
- Intrahepatic cholangiocarcinoma
- Weighted gene coexpression network analysis