Identification of Hub Genes and Pathways With Immune Cell Infiltration in Cholangiocarcinoma by Weighted Gene Co-expression Network Analysis

doi:10.21203/rs.3.rs-135007/v1

PDF

Research

Identification of Hub Genes and Pathways With Immune Cell Infiltration in Cholangiocarcinoma by Weighted Gene Co-expression Network Analysis

https://doi.org/10.21203/rs.3.rs-135007/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted 29 Dec, 2020

You are reading this latest preprint version

Background: Cholangiocarcinoma (CCA) is the most frequent tumor in biliary tract and the second most common primary tumor of the liver. However, the molecular biomarkers in tumorigenesis of CCA remain unclear. Therefore, we aim to explore the underlying mechanisms of progression and screen for novel prognostic biomarkers and therapeutic targets.

Method: The genes expression sequencing of normal and CCA samples were selected from the Gene Expression Omnibus database (GEO) and the Cancer Genome Atlas (TCGA). The weighted gene co-expression network analysis (WGCNA) was used to build the co-expression network. Gene ontology (GO) as well as Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were applied for the selected genes. The protein-protein interactions of these modules are visualized using cytoscape. Furthermore, the significance of these genes was confirmed by survival analysis. The tumor immune estimation resource (TIMER) was presented to investigate assess the relationship between the hub genes and the immune cells infiltration.

Results: Ten hub genes with CCA development were identified in this study containing CDC20, CCNA2, TOP2A, AURKA, CCNB2, UBE2C, NUSAP1, PRC1, PTTG1 and MCM4. Key genes of CCNB2 and PTTG1 might be potential prognostic biomarkers for CCA. GO analysis indicated the enrichment terms of nuclear division, collagen-containing extracellular matrix and cell adhesion molecule binding. KEGG analysis demonstrated that the cell cycle pathway was the significantly altered pathway. There was a negative correlation between TOP2A, AURKA, CCNB2, PRC1 expression and the infiltration of CD4+T cell, while MCM4 expression was positively associated with the infiltration of neutrophil cells. No significant association between CDC20 levels and CD4+T cell, CD8+T cell, B cell, neutrophil, macrophage, or dendritic cell infiltration in CCA, the same as CCNA2, UBE2C, NUSAP1, PTTG1 respectively.

Conclusion: These candidate genes may involve in the development of CCA. Our results offer novel insights into the etiology, prognosis, and treatment of CCA.

Translational Medicine

cholangiocarcinoma

weighted gene co-expression network analysis

hub gene

pathway

immune cell infiltration

Cholangiocarcinoma (CCA), a group of heterogeneous carcinomas characterized by cholangiocyte differentiation, is the most common biliary tract tumor and the second frequent primary hepatic malignancy^[¹^,²^]. The incidence of CCA has risen to 10-25% of all primary liver malignancies during the past 40 years and accounted for approximately 2% of cancer-related deaths globally per year^[³^,⁴^]. The diagnosis of CCA is challenging because of its atypical symptoms and anatomic location, and most patients are diagnosed with advanced CCA^[⁵^]. At present, the serum biomarkers routinely for CCA diagnosis, such as cancer antigen 125 and carbohydrate antigen 199, are not sensitive or accurate^[⁶^,⁷^]. Overall, the prognosis of CCA is extremely poor, with an average 5-year OS rate of 5 to 10 %^[⁸^]. In non-operable CCA the 5-year OS is less than 5%^[⁹^]. Consequently, the use of novel biomarkers for early detection and monitoring of CCA patients, as well as a better understanding of the molecular mechanisms associated with its development, are critical and may provide new insights for treatment.

In the last decade, gene expression analysis based on the next-generation sequencing (NGS) and high-throughput platforms has been widely used to detect significant expression changes in cancers ^[¹⁰^,¹¹^]. Several studies indicated that different genomic alterations were involved in the development of CCA and might have potential value for its diagnosis and prognosis. Frequent mutations in oncogenes such as KRAS, as well as cancer suppressor genes of TP53 were identified in CCA^[¹²^,¹³^].IDH1 mutation was also founded to be involved in downregulated genes with hypermethylation status in CCA patients^[¹⁴^,¹⁵^]. Although several genes and mechanisms have been demonstrated to be closely implicated in the development of CCA, the comprehensive profile of the whole genes of CCA is still unclear. Weighted gene co-expression network analysis (WGCNA), a powerful biology method, has been successfully used to identify hub genes in many malignancies to describe the co-expression patterns between genes across microarray samples. In this study, to further investigate the molecular mechanisms of CCA development, we used WGCNA to explore the hub genes of CCA and analyze its relationship with immune infiltration.

Data collection and screening

The gene expression profiles GSE26566 (GPL6104 Illumina humanRef-8 v2.0 expression beadchip) were collected from the Gene Expression Omnibus database (GEO, https://www.ncbi.nlm.nih.gov/geo/). This dataset included 6 normal samples and 104 CCA samples. To enhance the robustness, gene expression profile data and survival data of CCA were also downloaded from the Cancer Genome Atlas database (TCGA, https://www.cancer.gov/tcga), including 9 normal and 35 CCA samples. Differentially expressed genes (DEGs) between normal and CCA samples were identified using the R package “limma”.

Weighted Co-expression Network Construction and Module Detection

To further investigate the correlation of gene expression between the normal and CCA samples, we adopted WGCNA package in R software to build a co-expression network for GEO and TCGA databases respectively^[¹⁶^]. A hierarchical clustering tree was constructed and different branches of the tree were represented different gene modules. Then, the adjacency matrix was transformed into topological overlap matrix (TOM). Modules were identified using the dynamic tree cut algorithm based on the TOM-based dissimilarity measure. The pearson correlation examination was performed to explore the relation of module eigengenes (MEs).

Functional Annotation and Pathway Enrichment Analysis

The overlapping of DEGs and the selecting modules of WGCNA co-expression network in GEO and TCGA databases were portrayed using a Venn diagram. The biological function of the overlapping genes was elucidated through functional annotation and pathway enrichment analysis. The gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were annotated and visualized using the “ClusterProfiler package” in R software^[¹⁷^].

Protein-Protein Interaction Network Construction and Analysis

The search tool for the retrieval of interaction genes (STRING, https://string-db.org/) was used to construct the protein–protein interactions(PPI) networks^[¹⁸^]. The overlapping genes were mapped to STRING to evaluate the PPI information with a confidence score > 0.9 as the cut-off standard. Cytoscape is a practical open–source software tool for visually exploring biomolecule interaction networks consisting of proteins, genes and other types of interaction ^[¹⁹^].Topological features of the network were calculated to filter the hub network. The maximal clique centrality (MCC) methods in cytohubba plug-in in cytoscape was used to select the key genes. Finally, the top 10 genes in MCC method were selected as the hub genes^[²⁰^].

Survival Analysis

The prognostic role of hub genes was validated by survival analysis in TCGA database. R package “survival” was tested by log-rank tests and Kaplan-Meier survival curves were plotted ^[²¹^].

Tumor-Infiltrating Immune Cells

The tumor immune estimation resource 2.0 (TIMER 2.0, http://timer.cistrome.org/) database was presented to comprehensively assess the relationship between the hub genes and the immune cells infiltration^[²²^].

Statistical analysis

Statistical analyses were performed using IBM SPSS version 20.0 and the statistical software package R version 4.0.2. Absolute log2 fold-change >1 and P-value < 0.05 were considered statistically significant.

Data processing

The whole work of this study is shown in Figure 1. The gene expression profile dataset GSE26566 was downloaded from the GEO database containing 104 CCA samples and 6 normal samples. Then, 168 significantly up–regulated DEGs and 595 significantly down–regulated DEGs were identified. The TCGA dataset contained of 35 normal samples and 9 CCA samples, and a total of 2808 up-regulated genes and 2467 down-regulated genes were identified. Figure 2A and Figure 2B show the volcano plot of DGEs in GEO and TCGA databases, respectively.

Construction of the weighted gene co-expression network and identification of key modules

Gene co-expression modules were established in the GEO and TCGA cohort by using WGCNA (Figure 3A,3C). In the GEO dataset, the black and green modules were the top two modules closely related to CCA with correlation coefficients at 0.43 and -0.41 respectively (Figure 3B). In the TCGA cohort, it was found that the correlation coefficients of the turquoise modules were greatest (Figure 3D). Consequently, the intersection of the MEgreen module in GEO cohort and MEturquoise in the TCGA cohort combined with DEGs in these two datasets is shown in Figure4A, including 89 genes. These genes were regarded as the essential genes connected with CCA development and were selected for further analysis.

Functional annotation and pathway enrichment analysis

All 89 intersection genes were analyzed by GO and KEGG pathway enrichment analyses. The results demonstrated that the top biological processes associated with these genes was nuclear division. The cellular components most related to these genes was collagen-containing extracellular matrix. Furthermore, the molecular function terms were indicated that these genes were enriched in cell adhesion molecule binding (Figure 4B). In KEGG analysis, cell cycle pathway was the top pathway with significant enrichment (Figure 4C).

PPI analysis and identification of hub genes

The PPI network was constructed by the STRING tool, consisting of 89 nodes and 246 edges (Figure 5A). The hub network extracted from PPI network based on a confidence score > 0.9 (Figure 5B). The result from STRING database was imported in cytoscape for subsequent analysis. The genes that scored in the top 10 by MCC method in cytohubba plus-in were selected as hub genes of CAA. These genes were cell division cycle 20(CDC20), cyclin A2(CCNA2), type IIA topoisomerases(TOP2A), aurora kinase A (AURKA), cyclin B1(CCNB2), ubiquitin-coupled enzyme 2C(UBE2C), nucleolar and spindle-associated protein 1(NUSAP1), polycomb repressor complex 1(PRC1), pituitary tumor transforming gene 1(PTTG1) and minichromosomal maintenance proteins 4 (MCM4), which may play an important role in CCA progression (Figure 5C).

Survival Analysis and Tumor-Infiltrating Immune Cells

The survival analysis results of all hub genes were also validated in TCGA databases. Overexpression of CCNB2 and PTTG1 had significantly predictive effect on the prognosis of CCA, while the other genes had no significant predictive effect (Figure 6). We further explore the correlation between the hub genes expression and immune infiltration using the TIMER database^[²²^]. There was a negative correlation between TOP2A expression and the infiltration of CD4+T cell (Cor=-0.433, p = 9.43e-03), similarly between AURKA (Cor=-0.343, p = 4.35e-02), CCNB2 (Cor=-0.414, p = 1.34e-02), PRC1(Cor=-0.364, p = 3.14e-02) and CD4+T cell. MCM4 expression was positively associated with the infiltration of neutrophil cells (Cor=0.447, p = 7.13e-03). There was no significant association between CDC20 levels and CD4+T cell, CD8+T cell, B cell, neutrophil, macrophage, or dendritic cell infiltration in mesothelioma, the same in CCNA2, UBE2C, NUSAP1, PTTG1 respectively (Figure 7).

CCA is the most common malignancy of bile duct and its incidence is increasing annually^[²³^]. Due to its occult clinical characteristics, difficulties in early diagnosis as well as limited therapeutic approaches, the prognosis of CCA is still poor^[²⁴^]. Although surgical resection was regard as a possible option for certain patients who suffer CCA, the 5-year OS is extremely low^[²⁵^]. Therefore, it is urgent to investigate the pathogenesis of CCA, so as to develop new clinical treatment strategies.

In the present study, a total of 10 genes were identified as hub genes in the development of CCA, including CDC20, CCNA2, TOP2A, AURKA, CCNB2, UBE2C, NUSAP1, PRC1, PTTG1 and MCM4. The genes of CCNA2, CCNB2, CDC20 and UBE2C play a critical role in cell cycle progression. CCNA2 and CCNB2 are both the key members of the cyclin family ^[²⁶^,²⁷^]. Abnormal expression of CCNA2 or CCNB2 has been detected in various types of cancers. It was reported that CCNA2 was dysregulated in lung cancer, gastric cancer and colorectal carcinoma^[^28-30^]. Several studies also pointed the dysregulate of CCNB2 in hepatocellular carcinoma, non-small cell lung cancer and bladder cancer^[^31-33^]. However, the molecular mechanisms of these two genes in CCA remain unclear. CDC20 is an essential co-factor of the anaphase promoting complex (APC) as well as a downstream factor of the spindle assembly checkpoint^[³⁴^,³⁵^]. Previous studies have indicated that CDC20 may represent a promising therapeutic target for human cancers^[³⁶^,³⁷^]. Li et al indicated that CDC20 promoted the resistance of esophageal cancer chemotherapy by regulating E2F1 degradation and thymidylate synthase expression^[³⁸^]. Diane et al demonstrated that CDC20 overexpression increased the ability of human glioblastoma stem-like cells to generate brain tumors in an orthotopic xenograft model and implicated the role of CDC20-APC-targeted strategies in glioblastoma therapy^[³⁹^]. Shi M et al found that CDC20 played an important role in the development of hepatocellular carcinoma by governing prolyl hydroxylase 3 protein^[⁴⁰^]. Thus, exploring the underlying mechanism of CDC20 in the pathogenesis of cholangiocarcinoma and the development of CDC20 inhibitors may contribute to improve the therapeutic efficacy of patients with CCA. UBE2C is a ubiquitin enzyme that is highly expressed in various human cancers, such as NSCLC, ovarian cancers and breast cancers^[⁴¹^,⁴²^]. Overexpression of UBE2C may cause the loss of mitotic spindle checkpoint activity and genomic stability^[⁴³^], suggesting that UBE2C may associated with tumorigenesis of CCA.

DNA topoisomerases, especially TOP2A, demonstrate the ability to regulate the topological states of DNA during transcription. TOP2A has been indicated to be highly expressed in 33% patients with intrahepatic CCA^[⁴⁴^]. Overexpression of TOP2A in CCA was also confirmed in this paper. TOP2A has also been confirmed as a therapeutic target with anticancer and antibacterial effects^[⁴⁵^,⁴⁶^]. Therefore, we can speculate that TOP2A may also be a new therapeutic target for CCA. AURKA belongs to aurora family and plays key roles in regulating mitotic progression and cell proliferation^[⁴⁷^]. Aberrant AURKA activity participates in oncogenic transformation through the development of chromosomal instability and tumor cell heterogeneity^[⁴⁸^,⁴⁹^]. Ding et al has demonstrated that AURKA inhibitor can inhibit cell growth and promote apoptosis of CCA both in vitro and in vivo experiments. These results suggested that AURKA may be an effective target for the treatment of CCA^[⁵⁰^]. NUSAP1 is a microtubule-binding protein involving in chromosomal segregation^[⁵¹^,⁵²^]. NUSAP1 regulates metastasis of cervical carcinoma through Wnt/β-catenin signaling pathway^[⁵³^]. The down-regulation of NUSAP1 can inhibit cell proliferation, migration and invasion in NSCLC, which is related to the regulation of BTG2/PI3K/Akt signaling pathway^[⁵⁴^]. PRC1 is encoded by multiple homologous genes and shows the ability to promote ubiquitination of histone H2A^[⁵⁵^]. Su et al revealed that PRC1 can coordinate stemness with immune suppression and point out the potential clinical application of targeting PRC1 in double-negative prostate cancer^[⁵⁶^]. PTTG1, a ubiquitously expressed transcription factor, is involved in DNA damage repair, angiogenesis, cell differentiation and apoptosis^[⁵⁷^]. It can promote the transcription of genes that are directly or indirectly involved in tumorigenesis^[⁵⁸^]. Therefore, it is important to elucidate the mechanism of PTTG1 involvement in the development of CCA ^[⁵⁹^]. MCM4, a member of MCM family, is required to initiate genome replication in eukaryotes. It may be involved in the formation of replication forks and the recruitment of other DNA replication-related proteins^[⁶⁰^].

In GO analysis, nuclear division, collagen-containing extracellular matrix and cell adhesion molecule binding were the main pathways of gene enrichment. Furthermore, the genes were closely related to the cell cycle process in the KEGG analysis. This result was consistent with the related functions of the hub genes including CCNA2, CCNB2, CDC20 and UBE2C.Cell cycle dysregulation has been confirmed to be associated with the development and progression of various types of cancers, including CCA^[⁶¹^,⁶²^]. Chen et al found that the inhibition of Kelch-like family member 21 (KLHL21) could significantly reduce proliferation, migration and invasion of CCA cells. In addition, KLHL21 knockdown alleviated the activation of the Erk signaling pathway via decreasing the expression of phospho-Erk1/2, suggesting that it might be a potential therapeutic target for CCA^[⁶³^]. Several studies have shown that some drugs, such as telmisartan and cryptotanshinone, can induce cell cycle arrest and apoptosis in CCA and are potential therapeutic agents for CCA^[⁶⁴^,⁶⁵^]. Obviously, more attention should be paied to the cell cycle-related pathways and genes to explore the mechanism of development and treatment of CCA.

The immune system in CCA presents a complex landscape and most CCA is thought to origin from chronic inflammatory of the biliary tract^[⁶⁶^]. In this study, TOP2A, AURKA, PRC1, MCM4 were confirmed to be associated with immune cell infiltration using the TIMER 2.0 algorithm. Understanding the interaction between infiltrating immune cells and tumor cells in CCA may open new opportunities for treatment. Patients with neutrophil distribution in CCA may indicate poor prognosis^[⁶⁷^]. The chronic inflammatory status and the resulting immune state deserve further study and are meaningful for subsequent treatment of cancers. Tumor-associated macrophages (TAMs) play a pivotal role in promoting angiogenesis, tumor cell proliferation and suppression of adaptive immunity^[⁶⁸^]. TAMs drive the occurrence of CCA by producing various soluble mediators, including cytokines, chemokine and reactive nitrogen intermediates^[⁶⁹^]. It has been proven to be a prognostic marker and potential therapeutic intervention in CCA^[⁷⁰^]. Moreover, depletion of the enrichment of CD163+TAMs, CD4+ T cell, as well as CD8+ T cell, are correlated with reducing relapse-free survival and inducing gemcitabine resistance^[⁷¹^]. Mature dendritic cell may enhance CD4+ and CD8+ T cells infiltration into cancers and improve prognosis in patients with CCA^[⁷²^]. B cells have been identified in the tumor infiltrating lymphocytes of CCA, but are rarely observed in patient tissues^[⁷³^,⁷⁴^]. Goeppert et al showed that the total number of tumor-infiltrating B cell is correlated with a longer OS probability in biliary tract cancer. Whether this observation just reflected an occasional phenomenon in the context of immune response activation, or indicated of the true effect of B cell on tumor control should be determined in future studies.

There are still several limitations in this study. First, although a relatively large number of CCA samples were analyzed, more studies with larger samples sizes are needed to verify the findings of this study. Second, although the immune cell infiltration results of the hub genes were evaluated in this study, they were based only on the information from an online database. Consequently, further experiments should be performed to validate the results in this study and to reveal the mechanisms of these hub genes in CCA progression and prognosis.

In summary, by using WGCNA and a series of other bioinformatics approaches, we identified 10 hub genes (CDC20, CCNA2, TOP2A, AURKA, CCNB2, UBE2C, NUSAP1, PRC1, PTTG1 and MCM4), which may contribute to the development and prognosis of CCA. Further studies of the hub genes identified in the present study are urgently required to determine the underlying mechanisms associated with their function in CCA.

CCA, Cholangiocarcinoma; OS, overall survival; NGS, next generation sequencing; WGCNA, Weighted gene co-expression network analysis; GEO, The Gene Expression Omnibus database; TCGA, The Cancer Genome Atlas database; DEGs, differentially expressed genes; TOM, topological overlap matrix; MEs, module eigengenes; GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; STRING, Search the Tool for The Retrieval of Interaction Genes; PPI, Protein–protein interaction; GEPIA, Gene Expression Profiling Interactive Analysis; TIMER, The Tumor Immune Estimation Resource; HPA, The Human Protein Atlas; CCNA2, Cyclin A2; CCNB2, cyclin B1; CDC20, Cell division cycle 20; APC, anaphase promoting complex; TOP2A, type IIA topoisomerases; AURKA, Aurora kinase; UBE2C, Ubiquitin-conjugating enzyme 2C; NUSAP1, Nucleolar and spindle-associated protein 1; PRC1, Polycomb Repressor Complex 1; PTTG1, Pituitary tumor-transforming gene 1. MCM4, Minichromosomal maintenance proteins 4.

Acknowledgements

All the authors indicate no potential conflicts of interests. In addition, Sujing Jiang wants to express her thanks to Miss.Yolo. Sincerely wish her happiness.

Author contributions

YY and SJ conceived the concept and designed the study. SJ and ZL collected and analysised the data. HS and GL prepared the figures. YY and SJ wrote the paper. All authors reviewed and approved the manuscript.

Funding

No funding for this article.

Availability of data and materials

Not applicable.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing Interests

The authors declare that they have no conflicts of interest.

[1] Blechacz B. Cholangiocarcinoma: Current Knowledge and New Developments[J]. Gut Liver, 2017, 11(1): 13-26.

[2] Chan E, Berlin J. Biliary tract cancers: understudied and poorly understood[J]. J Clin Oncol, 2015, 33(16): 1845-8.

[3] Tyson G L, El-Serag H B. Risk factors for cholangiocarcinoma[J]. Hepatology, 2011, 54(1): 173-84.

[4] Razumilava N, Gores G J. Cholangiocarcinoma[J]. Lancet, 2014, 383(9935): 2168-79.

[5] Rizvi S, Gores G J. Pathogenesis, diagnosis, and management of cholangiocarcinoma[J]. Gastroenterology, 2013, 145(6): 1215-29.

[6] Patel A H, Harnois D M, Klee G G, et al. The utility of CA 19-9 in the diagnoses of cholangiocarcinoma in patients without primary sclerosing cholangitis[J]. Am J Gastroenterol, 2000, 95(1): 204-7.

[7] Blechacz B, Gores G J. Cholangiocarcinoma: advances in pathogenesis, diagnosis, and treatment[J]. Hepatology, 2008, 48(1): 308-21.

[8] Doherty B, Nambudiri V E, Palmer W C. Update on the Diagnosis and Treatment of Cholangiocarcinoma[J]. Curr Gastroenterol Rep, 2017, 19(1): 2.

[9] Macias R I R, Banales J M, Sangro B, et al. The search for novel diagnostic and prognostic biomarkers in cholangiocarcinoma[J]. Biochim Biophys Acta Mol Basis Dis, 2018, 1864(4 Pt B): 1468-1477.

[10] Kulasingam V, Diamandis E P. Strategies for discovering novel cancer biomarkers through utilization of emerging technologies[J]. Nat Clin Pract Oncol, 2008, 5(10): 588-99.

[11] Nones K, Patch A M. The Impact of Next Generation Sequencing in Cancer Research[J]. Cancers (Basel), 2020, 12(10).

[12] Chan-On W, Nairismägi M L, Ong C K, et al. Exome sequencing identifies distinct mutational patterns in liver fluke-related and non-infection-related bile duct cancers[J]. Nat Genet, 2013, 45(12): 1474-8.

[13] Zou S, Li J, Zhou H, et al. Mutational landscape of intrahepatic cholangiocarcinoma[J]. Nat Commun, 2014, 5: 5696.

[14] Borger D R, Tanabe K K, Fan K C, et al. Frequent mutation of isocitrate dehydrogenase (IDH)1 and IDH2 in cholangiocarcinoma identified through broad-based tumor genotyping[J]. Oncologist, 2012, 17(1): 72-9.

[15] Zhang X, Miao R, Liu T, et al. IDH1 as a frequently mutated gene has potential effect on exosomes releasement by epigenetically regulating P2RX7 in intrahepatic cholangiocarcinoma[J]. Biomed Pharmacother, 2019, 113: 108774.

[16] Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis[J]. BMC Bioinformatics, 2008, 9: 559.

[17] Yu G, Wang L G, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters[J]. Omics, 2012, 16(5): 284-7.

[18] Szklarczyk D, Gable A L, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets[J]. Nucleic Acids Res, 2019, 47(D1): D607-d613.

[19] Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks[J]. Genome Res, 2003, 13(11): 2498-504.

[20] Bader G D, Hogue C W. An automated method for finding molecular complexes in large protein interaction networks[J]. BMC Bioinformatics, 2003, 4: 2.

[21] Tang Z, Li C, Kang B, et al. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses[J]. Nucleic Acids Res, 2017, 45(W1): W98-w102.

[22] Li T, Fan J, Wang B, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells[J]. Cancer Res, 2017, 77(21): e108-e110.

[23] Torre L A, Bray F, Siegel R L, et al. Global cancer statistics, 2012[J]. CA Cancer J Clin, 2015, 65(2): 87-108.

[24] Blechacz B, Komuta M, Roskams T, et al. Clinical diagnosis and staging of cholangiocarcinoma[J]. Nat Rev Gastroenterol Hepatol, 2011, 8(9): 512-22.

[25] Turgeon M K, Maithel S K. Cholangiocarcinoma: a site-specific update on the current state of surgical management and multi-modality therapy[J]. Chin Clin Oncol, 2020, 9(1): 4.

[26] Zhang S, Tischer T. Cyclin A2 degradation during the spindle assembly checkpoint requires multiple binding modes to the APC/C[J], 2019, 10(1): 3863.

[27] Nam H J, Van Deursen J M. Cyclin B2 and p53 control proper timing of centrosome separation[J]. Nat Cell Biol, 2014, 16(6): 538-49.

[28] Cooper W A, Kohonen-Corish M R, Mccaughan B, et al. Expression and prognostic significance of cyclin B1 and cyclin A in non-small cell lung cancer[J]. Histopathology, 2009, 55(1): 28-36.

[29] Nozoe T, Inutsuka S, Honda M, et al. Clinicopathologic significance of cyclin A expression in colorectal carcinoma[J]. J Exp Clin Cancer Res, 2004, 23(1): 127-33.

[30] Lee Y, Lee C E, Oh S. Pharmacogenomic Analysis Reveals CCNA2 as a Predictive Biomarker of Sensitivity to Polo-Like Kinase I Inhibitor in Gastric Cancer[J], 2020, 12(6).

[31] Li R, Jiang X, Zhang Y, et al. Cyclin B2 Overexpression in Human Hepatocellular Carcinoma is Associated with Poor Prognosis[J]. Arch Med Res, 2019, 50(1): 10-17.

[32] Qian X, Song X, He Y, et al. CCNB2 overexpression is a poor prognostic biomarker in Chinese NSCLC patients[J]. Biomed Pharmacother, 2015, 74: 222-7.

[33] Lei C Y, Wang W, Zhu Y T, et al. The decrease of cyclin B2 expression inhibits invasion and metastasis of bladder cancer[J]. Urol Oncol, 2016, 34(5): 237.e1-10.

[34] Kapanidou M, Curtis N L, Bolanos-Garcia V M. Cdc20: At the Crossroads between Chromosome Segregation and Mitotic Exit[J]. Trends Biochem Sci, 2017, 42(3): 193-205.

[35] Schrock M S, Stromberg B R, Scarberry L, et al. APC/C ubiquitin ligase: Functions and mechanisms in tumorigenesis[J]. Semin Cancer Biol, 2020.

[36] Wang Z, Wan L, Zhong J, et al. Cdc20: a potential novel therapeutic target for cancer treatment[J]. Curr Pharm Des, 2013, 19(18): 3210-4.

[37] Wang L, Zhang J, Wan L, et al. Targeting Cdc20 as a novel cancer therapeutic strategy[J]. Pharmacol Ther, 2015, 151: 141-51.

[38] Li B, Xu W W, Guan X Y, et al. Competitive Binding Between Id1 and E2F1 to Cdc20 Regulates E2F1 Degradation and Thymidylate Synthase Expression to Promote Esophageal Cancer Chemoresistance[J]. Clin Cancer Res, 2016, 22(5): 1243-55.

[39] Mao D D, Gujar A D, Mahlokozera T, et al. A CDC20-APC/SOX2 Signaling Axis Regulates Human Glioblastoma Stem-like Cells[J]. Cell Rep, 2015, 11(11): 1809-21.

[40] Shi M, Dai W Q, Jia R R, et al. APC(CDC20)-mediated degradation of PHD3 stabilizes HIF-1a and promotes tumorigenesis in hepatocellular carcinoma[J]. Cancer Lett, 2020, 496: 144-155.

[41] Guo J, Wu Y, Du J, et al. Deregulation of UBE2C-mediated autophagy repression aggravates NSCLC progression[J]. Oncogenesis, 2018, 7(6): 49.

[42] Kim Y J, Lee G, Han J, et al. UBE2C Overexpression Aggravates Patient Outcome by Promoting Estrogen-Dependent/Independent Cell Proliferation in Early Hormone Receptor-Positive and HER2-Negative Breast Cancer[J]. Front Oncol, 2019, 9: 1574.

[43] Bajaj S, Alam S K, Roy K S, et al. E2 Ubiquitin-conjugating Enzyme, UBE2C Gene, Is Reciprocally Regulated by Wild-type and Gain-of-Function Mutant p53[J]. J Biol Chem, 2016, 291(27): 14231-47.

[44] Potkonjak M, Miura J T, Turaga K K, et al. Intrahepatic cholangiocarcinoma and gallbladder cancer: distinguishing molecular profiles to guide potential therapy[J]. HPB (Oxford), 2015, 17(12): 1119-23.

[45] Delgado J L, Hsieh C M, Chan N L, et al. Topoisomerases as anticancer targets[J]. Biochem J, 2018, 475(2): 373-398.

[46] Hiasa H. DNA Topoisomerases as Targets for Antibacterial Agents[J]. Methods Mol Biol, 2018, 1703: 47-62.

[47] Glover D M, Leibowitz M H, Mclean D A, et al. Mutations in aurora prevent centrosome separation leading to the formation of monopolar spindles[J]. Cell, 1995, 81(1): 95-105.

[48] Maia A R, Van Heesbeen R G, Medema R H. A growing role for Aurora A in chromosome instability[J]. Nat Cell Biol, 2014, 16(8): 739-41.

[49] Goldenson B, Crispino J D. The aurora kinases in cell cycle and leukemia[J]. Oncogene, 2015, 34(5): 537-45.

[50] Ding X, Huang T, Peng C, et al. Therapeutic Rationale to Target Highly Expressed Aurora kinase A Conferring Poor Prognosis in Cholangiocarcinoma[J]. J Cancer, 2020, 11(8): 2241-2251.

[51] Raemaekers T, Ribbeck K, Beaudouin J, et al. NuSAP, a novel microtubule-associated protein involved in mitotic spindle organization[J]. J Cell Biol, 2003, 162(6): 1017-29.

[52] Li C, Xue C, Yang Q, et al. NuSAP governs chromosome oscillation by facilitating the Kid-generated polar ejection force[J]. Nat Commun, 2016, 7: 10597.

[53] Li H, Zhang W, Yan M, et al. Nucleolar and spindle associated protein 1 promotes metastasis of cervical carcinoma cells by activating Wnt/β-catenin signaling[J]. J Exp Clin Cancer Res, 2019, 38(1): 33.

[54] Xu Z, Wang Y, Xiong J, et al. NUSAP1 knockdown inhibits cell growth and metastasis of non-small-cell lung cancer through regulating BTG2/PI3K/Akt signaling[J], 2020, 235(4): 3886-3893.

[55] Wang H, Wang L, Erdjument-Bromage H, et al. Role of histone H2A ubiquitination in Polycomb silencing[J]. Nature, 2004, 431(7010): 873-8.

[56] Su W, Han H H, Wang Y, et al. The Polycomb Repressor Complex 1 Drives Double-Negative Prostate Cancer Metastasis by Coordinating Stemness and Immune Suppression[J]. Cancer Cell, 2019, 36(2): 139-155.e10.

[57] Tong Y, Eigler T. Transcriptional targets for pituitary tumor-transforming gene-1[J]. J Mol Endocrinol, 2009, 43(5): 179-85.

[58] Yoon C H, Kim M J, Lee H, et al. PTTG1 oncogene promotes tumor malignancy via epithelial to mesenchymal transition and expansion of cancer stem cell population[J]. J Biol Chem, 2012, 287(23): 19516-27.

[59] Chen R, Duan J, Li L, et al. mTOR promotes pituitary tumor development through activation of PTTG1[J]. Oncogene, 2017, 36(7): 979-988.

[60] Ishimi Y, Komamura-Kohno Y, Kwon H J, et al. Identification of MCM4 as a target of the DNA replication block checkpoint system[J]. J Biol Chem, 2003, 278(27): 24644-50.

[61] Malumbres M, Barbacid M. Cell cycle, CDKs and cancer: a changing paradigm[J]. Nat Rev Cancer, 2009, 9(3): 153-66.

[62] Briggs C D, Neal C P, Mann C D, et al. Prognostic molecular markers in cholangiocarcinoma: a systematic review[J]. Eur J Cancer, 2009, 45(1): 33-47.

[63] Chen J, Song W, Du Y, et al. Inhibition of KLHL21 prevents cholangiocarcinoma progression through regulating cell proliferation and motility, arresting cell cycle and reducing Erk activation[J]. Biochem Biophys Res Commun, 2018, 499(3): 433-440.

[64] Samukawa E, Fujihara S, Oura K, et al. Angiotensin receptor blocker telmisartan inhibits cell proliferation and tumor growth of cholangiocarcinoma through cell cycle arrest[J]. Int J Oncol, 2017, 51(6): 1674-1684.

[65] Ke F, Wang Z, Song X, et al. Cryptotanshinone induces cell cycle arrest and apoptosis through the JAK2/STAT3 and PI3K/Akt/NFκB pathways in cholangiocarcinoma cells[J]. Drug Des Devel Ther, 2017, 11: 1753-1766.

[66] Augustine M M, Fong Y. Epidemiology and risk factors of biliary tract and primary liver tumors[J]. Surg Oncol Clin N Am, 2014, 23(2): 171-88.

[67] Mao Z Y, Zhu G Q, Xiong M, et al. Prognostic value of neutrophil distribution in cholangiocarcinoma[J]. World J Gastroenterol, 2015, 21(16): 4961-8.

[68] Solinas G, Germano G, Mantovani A, et al. Tumor-associated macrophages (TAM) as major players of the cancer-related inflammation[J]. J Leukoc Biol, 2009, 86(5): 1065-73.

[69] Prakobwong S, Yongvanit P, Hiraku Y, et al. Involvement of MMP-9 in peribiliary fibrosis and cholangiocarcinogenesis via Rac1-dependent DNA damage in a hamster model[J]. Int J Cancer, 2010, 127(11): 2576-87.

[70] Hasita H, Komohara Y, Okabe H, et al. Significance of alternatively activated macrophages in patients with intrahepatic cholangiocarcinoma[J]. Cancer Sci, 2010, 101(8): 1913-9.

[71] Kitano Y, Okabe H, Yamashita Y I, et al. Tumour-infiltrating inflammatory and immune cells in patients with extrahepatic cholangiocarcinoma[J]. Br J Cancer, 2018, 118(2): 171-180.

[72] Takagi S, Miyagawa S, Ichikawa E, et al. Dendritic cells, T-cell infiltration, and Grp94 expression in cholangiocellular carcinoma[J]. Hum Pathol, 2004, 35(7): 881-6.

[73] Goeppert B, Frauenschuh L, Zucknick M, et al. Prognostic impact of tumour-infiltrating immune cells on biliary tract cancer[J]. Br J Cancer, 2013, 109(10): 2665-74.

[74] Kasper H U, Drebber U, Stippel D L, et al. Liver tumor infiltrating lymphocytes: comparison of hepatocellular and cholangiolar carcinoma[J]. World J Gastroenterol, 2009, 15(40): 5053-7.

PDF

Version 1

posted 29 Dec, 2020

You are reading this latest preprint version

Identification of Hub Genes and Pathways With Immune Cell Infiltration in Cholangiocarcinoma by Weighted Gene Co-expression Network Analysis

Status:

Version 1

Abstract

Figures

Introduction

Material And Methods

Results

Discussion

Conclusion

Abbreviations

Declarations

References

Status:

Version 1