key: cord-0019809-v8733ubv authors: Kao, Tzu-Jen; Wu, Chung-Che; Phan, Nam Nhut; Liu, Yen-Hsi; Ta, Hoang Dang Khoa; Anuraga, Gangga; Wu, Yung-Fu; Lee, Kuen-Haur; Chuang, Jian-Ying; Wang, Chih-Yang title: Prognoses and genomic analyses of proteasome 26S subunit, ATPase (PSMC) family genes in clinical breast cancer date: 2021-07-30 journal: Aging (Albany NY) DOI: 10.18632/aging.203345 sha: 0e9a5de143051f76ce46e751a35cf47d45cb1ad9 doc_id: 19809 cord_uid: v8733ubv Breast cancer is a complex disease, and several processes are involved in its development. Therefore, potential therapeutic targets need to be discovered for these patients. Proteasome 26S subunit, ATPase gene (PSMC) family members are well reported to be involved in protein degradation. However, their roles in breast cancer are still unknown and need to be comprehensively researched. Leveraging publicly available databases, such as cBioPortal and Oncomine, for high-throughput transcriptomic profiling to provide evidence-based targets for breast cancer is a rapid and robust approach. By integrating the aforementioned databases with the Kaplan–Meier plotter database, we investigated potential roles of six PSMC family members in breast cancer at the messenger RNA level and their correlations with patient survival. The present findings showed significantly higher expression profiles of PSMC2, PSMC3, PSMC4, PSMC5, and PSMC6 in breast cancer compared to normal breast tissues. Besides, positive correlations were also revealed between PSMC family genes and ubiquinone metabolism, cell cycle, and cytoskeletal remodeling. Meanwhile, we discovered that high levels of PSMC1, PSMC3, PSMC4, PSMC5, and PSMC6 transcripts were positively correlated with poor survival, which likely shows their importance in breast cancer development. Collectively, PSMC family members have the potential to be novel and essential prognostic biomarkers for breast cancer development. In 2020, breast cancer accounted for 30% of all types of cancer in women in the United States. Expressions of the estrogen receptor (ER), progesterone receptor, and human epidermal growth factor receptor (HER)-2 are used to subgroup breast cancer cases. Currently, salvage therapy for breast cancer patients includes fulvestrant (selective ER downregulators) [1, 2] , cyclin-dependent kinase 4/6 inhibitors [3] , aromatase inhibitors combined with everolimus (a mammalian analog of rapamycin which acts as a mammalian target of rapamycin (mTOR) inhibitor) [4] , and histone deacetylase (HDAC) inhibitors [5] . High expression of B-cell lymphoma 2 was detected in nearly 70% of metastatic breast cancer patients, and treatment with a selective inhibitor improved apoptosis in a preclinical model of breast cancer [6, 7] . Meanwhile, proteasome 26S subunit ATPase (PSMC), proteasome 20S subunit beta (PSMB), GATA-binding protein, serine/threonine kinase, and matrix metallopeptidase family genes, signal transducer and activator of transcription (STAT), Notch, and phosphatidylinositol 3kinase (PI3K) were reported to be causes of those alterations [8] [9] [10] [11] . The PSMC family is comprised of six members, namely PSMC1, PSMC2, PSMC3, PSMC4, PSMC5, and PSMC6, which partially constitute formation of the 19S regulatory complex. This complex plays an important role in regulating the 26S proteasome, which in turn, catalyzes the unfolding and translocation of substrates into the 20S proteasome. In addition, members of the PSMC gene family, except for PSMC3, are known to cause N-CoR degradation [12] . Previous studies showed that PSMC6 promotes osteoblast apoptosis and cancer cell proliferation, while PSMC2 inhibits apoptosis. Furthermore, PSMC6 also inhibited activation of the PI3K/AKT signaling pathway in an animal model of ovariectomy-induced osteoporosis [13] [14] [15] . PSMC5 participates in degradation of Tln1 and angiogenesis [16] . In hepatocellular carcinoma cells, knockdown of PSMC3IP resulted in suppression of xenograft proliferation and tumorigenesis [17] . High-throughput technologies are widely used as systemic approaches to explore differences in expressions of thousands of genes for both biological and genomics systems [18] [19] [20] . It is well recognized that many upregulated and downregulated genes are associated with oncogenic or tumor-suppressive functions in cancer development [21] [22] [23] [24] [25] [26] . Nevertheless, a holistic approach to exploring messenger (m)RNA levels of the entire PSMC family in breast cancer has not been conducted. Therefore, in the present study, we analyzed all available mRNA data from public breast cancer databases, comparing datasets from breast cancer patients with those from normal tissues. We also predicted interactive networks and gene regulatory networks related to the PSMC family to determine potential biomarkers. A meta-analysis approach was adopted to screen downstream molecules associated with PSMC genes. Based on our analysis, PSMC family members and their downstream-regulated genes are potential candidates for new therapeutic targets in breast cancer progression. Previous studies identified six PSMC family members in Homo species, and some of these genes play crucial roles in cancer progression. Oncomine platform contained a total of 392 unique analyses for PSMC1 expression, and PSMC1 had significant in 13 of 392 unique analyses. PSMC2 had significant in 55 of 433 unique analyses, PSMC3 had significant in 28 of 421 unique analyses, PSMC4 had significant in 71 of 432 unique analyses, PSMC5 had significant in 24 of 420 unique analyses, PSMC6 had significant in 28 of 445 unique analyses ( Figure 1 ). However, a meta-analysis is needed to clarify gene signatures of PSMC family members in breast cancer. According to our results from an Oncomine analysis of mRNA expressions of PSMC2, PSMC3, PSMC4, PSMC5, and PSMC6, these members are highly upregulated in breast cancer tissues; therefore, we chose breast cancer to perform further bioinformatics analyses ( Figure 1 ). Furthermore, in the METABRIC database, expressions of PSMC members in breast cancer tissues were significantly higher than those in normal tissues; p values ranged from 2.16E-45 to 0.023 for PSMC1, 1.37E-29 to 0.016 for PSMC2, 3.18E-21 to 0.001 for PSMC3, 1.28E-53 to 0.018 for PSMC4, 7.02E-36 to 0.041 for PSMC5, and 9.03E-12 to 0.039 for PSMC6 (Supplementary Table 1) . Meanwhile, to further explore gene expressions of the entire PSMC family in breast cancer, we compared transcript levels of different breast cancer subtypes, such as the triple-negative, HER-transcripts were positively correlated with poor survival, which likely shows their importance in breast cancer development. Collectively, PSMC family members have the potential to be novel and essential prognostic biomarkers for breast cancer development. AGING 2, and luminal subtypes, relative to normal breast tissues, in TCGA database (Supplementary Figure 1) . Interestingly, we discovered that PSMC genes were overexpressed in a subtype-specific manner: specifically, PSMC1, PSMC2, PSMC3, PSMC4 were highly expressed in the triple-negative subtype, PSMC5 in HER-2, and PSMC6 in luminal cancer. These results suggest oncogenic effects of PSMC family genes on tumor progression. GEPIA2 datasets were used to analyze mRNA levels of PSMC members in breast cancer tissues compared to normal tissues. We found that levels of all six PSMCs were upregulated in breast cancer tissues relative to normal breast tissues (Figure 2A-2H) . Additionally, analysis of the CCLE dataset (https://www.broadinstitute.org/ccle) also showed differential expressions of PSMC family members in breast cell lines ( Figure 3 ). We analyzed genes coexpressed with PSMC1 in the Perou Breast 2 dataset from the Oncomine platform. We found that PSMC1 was positively correlated with MYL9, PCOLCE, ANXA6, DVL1, MAPKAPK2, PROC, FTL, ZNF358, CHRNA1, COL5A1, MPP1, PDE6A, COL1A1, MIA, SH2B3, COL6A1, and BGN. We used the Landemaine dataset to analyze genes coexpressed with PSMC2 and found that its expression was positively correlated with SMAD5, found that its expression was positively correlated with CLPX, CCDC9OB, FAM18B2, C60rf62, ZBTB33, PYROXD1, CDC42SE2, COMMD6, LOC401397, CAPZAI, TPRKB, GABPA, MATR3 , ZDHHC20, SCOC, and COPS2 ( Figure 4A ). Additionally, associations among PSMC1, PSMC2, PSMC3, PSMC4, PSMC5, and PSMC6 were also analyzed using the GEPIA dataset. Specifically, PSMC1 was positively correlated with PSMC2 (R=0. 41 AGING PSMC6 (R=0.15, p<0.05). PSMC4 was positively correlated with PSMC5 (R=0.021, p<0.05) and PSMC6 (R=0.068, p<0.05). Finally, PSMC5 was positively correlated with PSMC6 (R=0.13, p<0.05) ( Figure 4B ). Meanwhile, we obtained similar results from the cBioPortal and the Cytoscaped and METABRIC databases, which revealed that the six PSMC members were correlated with cell cycle-related genes (Supplementary Figure 2 ). In addition, expressions of PSMC family members were also correlated with immune infiltration profiles in breast cancer, as analyzed with the Tumor Immune Estimation Resource (TIMER; cistrome.shinyapps.io/timer) tool. Expression of each PSMC gene was associated with tumor purity and markers of different types of immune cells (Supplementary Figure 3) . Since PSMC family genes were differentially expressed in samples from breast cancer patients, we next explored the potential roles of these genes in human breast cancer tissues, correlating their expressions with other potential biomarkers related to molecular subtypes of breast cancer. To determine expressions of PSMC AGING family members and their clinical relevance, we analyzed protein expressions of individual PSMC members in clinical specimens from the Human Protein Atlas. The data demonstrated that PSMC1-6 presented moderate protein expressions, and PSMC2, PSMC3, and PSMC5 were highly expressed in certain clinical tissues from breast cancer specimens ( Figure 5 ). The Kaplan-Meier plotter database also AGING showed that PSMC1, PSMC3, PSMC4, PSMC5, and PSMC6 had high expression levels in breast cancer tissues may have oncogenic roles in breast cancer progression. High transcription levels of PSMC1, PSMC3, PSMC4, PSMC5, and PSMC6 predicted poor survival, whereas PSMC2 did not show the same pattern ( Figure 6 ). First, to explore the universally regulated pathways of the entire PSMC family, GeneGo Metacore was leveraged to investigate downstream networks according to coexpression patterns of PSMC genes. An AGING analysis on the GeneGo Metacore platform demonstrated that genes coexpressed with the six PSMC genes participated in biological processes related to cancer progression. MetaCore can be used to construct downstream networks associated with biological processes from uploaded genes. By uploading PSMC coexpressed genes from the METABRIC database into the Metacore platform, we found that several cancer progression-related pathways were correlated with genes of the PSMC family (Supplementary Figures 4, 5 and Supplementary Table 2), including "Cytoskeleton remodeling_Regulation of actin cytoskeleton organization by the kinase effectors of Rho GTPases", "Cell cycle_Role of APC in cell cycle regulation", "Cell cycle_Chromosome condensation in prometaphase", "Cell cycle_ Nucleocytoplasmic transport of CDK/Cyclins", and "Transcription_Role of heterochromatin protein 1 family in transcriptional silencing" (Figure 7) . Next, the STRING platform was used to externally validate and search for potential protein-protein interactions (PPIs). The resulting network with a core cluster contained all of the genes associated with cancer progression and metastasis ( Figure 8 ). Next, we explored whether individual genes of the PSMC family regulate specific pathways and networks in breast cancer development. We obtained coexpression profiles for PSMC1 from TCGA and METABRIC breast cancer datasets. Afterward, GeneGo Metacore annotations of each biological process suggested that genes coexpressed with PSMC1 were involved inG-protein-coupled receptor (GPCR)-and apoptosis-related pathways and networks such as AGING "Chemotaxis_Lysophosphatidic acid signaling via GPCRs", "Development_Positive regulation of WNT/Beta-catenin signaling in the cytoplasm", and "Apoptosis and survival_Regulation of apoptosis by mitochondrial proteins", and play essential roles in breast cancer (Supplementary Figure 6 and Supplementary Table 3) . PSMC2-related genes were involved in Wnt-and hypoxia-related pathways and networks such as "Development_Negative regulation of WNT/Beta-catenin signaling in the cytoplasm" and "Transcription_HIF-1 targets", which may be involved in breast cancer (Supplementary Figure 7 and Supplementary Table 4 ). Genes coexpressed with PSMC3 participated in processes of cytoskeleton-and organization-related pathways and networks such as "Cytoskeleton remodeling_Regulation of actin cytoskeleton organization by the kinase effectors of Rho GTPases" (Supplementary Figure 8 and Supplementary Table 5 ). Genes associated with the proteasome 26S subunit, ATPase (PSMC) family were uploaded to the STRING platform to establish the network. Using k-means clustering, the network was further separated into different clusters. AGING PSMC4-related genes were involved in mitogen-activated protein kinase (MAPK)-and inflammation-related pathways and networks such as "Signal trans-duction_CXCR4 signaling via MAPKs cascades" and "Signal transduction_Angiotensin II/AGTR1 signaling via Notch, Beta-catenin and NF-κB pathways", which may participate in breast cancer (Supplementary Figure 9 and Supplementary Table 6 ). Genes found to be coexpressed with PSMC5 were involved in oxidative stress-and cell adhesion-related pathways and networks such as "Oxidative stress_ROS-induced cellular signaling" and "Cell adhesion_Tight junctions" (Supplementary Figure 10 and Supplementary Table 7) . PSMC6-coexpressed genes were involved in calcium-and hormone-related pathways and networks such as "Signal transduction_Calcium-mediated signaling" and "Reproduction_Gonadotropin-releasing hormone (GnRH) signaling ", which could participate in breast cancer (Supplementary Figure 11 and Supplementary Table 8 ). Breast cancer has the highest prevalence rate compared to other types of cancer, particularly in females. Despite several years of extraordinary efforts to increase our knowledge of tumor biology and improve surgical treatments and chemotherapies, prognoses of advanced breast cancer patients have not improved [27] [28] [29] [30] [31] [32] . Therefore, it is very important to investigate new diagnostic tools and novel biomarkers that can allow us to refine patient prognoses and investigate effective interventions. Most genes of the PSMC family are upregulated in many types and subtypes of cancer. PSMC members were proven to be involved in tumor progression. For example, overexpression of PFN1 is associated with PSMC1 in the MDA-MB-231 triple-negative breast cancer cell line and may involve multiple mechanisms for cancer progression [33] . PSMC2 is highly expressed in pancreatic cancer, and PSMC2-knockdown significantly decreased cell proliferation. PSMC3 was identified as a crucial node in a PPI network in glioma cells [34] . PSMC4 and PSMC5 contribute to prostate tumorigenesis [35] . Additionally, PSMC4 was identified as one of the best biomarkers for endometrial cancer [36] , and risk model construction revealed that it is also a prognostic marker for the same type of cancer [37] . PSMC5 acts as a novel regulator and is involved in the extracellular signal-regulated kinase 1/2 signaling pathway [38] , and it was observed to have a relatively higher cytoplasmic expression pattern in most cancer types [39] . Both the PSMC6 and MAPK8 genes were upregulated in melanosis coli patients [40] , and CRISPR Genome-Wide Screening demonstrated that the PSMC6 subunit is an important and sensitive target for bortezomib in multiple myeloma cells [41] . Since the roles of PSMC family members in breast cancer are poorly described, the present findings show their importance, by providing preliminary clues for prospective studies in breast cancer research. Findings from the current study are in line with previous reports on the roles of PSMC genes in cancer. Both mRNA and protein levels of PSMC2, PSMC3, PSMC4, PSMC5, and PSMC6 were significant in cancer tissues, and PSMC1, PSMC3, PSMC4, PSMC5, and PSMC6 overexpression was associated with poor prognoses of breast cancer patients. Meanwhile, to further clarify its role of PSMC family genes in pan-cancer, we also used GEPIA2 database to confirm that PSMC family genes had prognostic value in these integrated analyses (Supplementary Figure 12) , as well as the PSMC family gene expressions in CCLE database (Supplementary Figure 13) , The coexpression analysis revealed the positive correlative roles of PSMC family genes in cytoskeletal remodeling and CDK/cyclins, as well as cell cyclerelated pathways and networks, which are consistent with previous studies. To some extent, this is the first report on both mRNA and protein expressions of the PSMC family in cancer cell lines and tissues, together with their associations with breast cancer patient survival. Collectively, by integrating multiple high-throughput databases, our study uncovered that PSMC genes have prognostic and predictive value in breast cancer. To comprehensively provide a complete picture of the PSMC members not only in breast but also in other types and subtypes of cancer, our results can be used as hints for further examination of this family, and possibly they can serve as novel biomarkers and potential prognostic indicators in breast cancer. Oncomine (https://www.oncomine.org/), a wellknown high-throughput database for mRNA, was used to query expressions of PSMC family members [42] . In brief, each gene symbol of the PSMC family was used to search for expression levels in 20 types of cancer relative to matched normal-type samples. The search thresholds included a multiple of change of 2fold, a p value of <0.01, and a gene ranking in the top 10%. The search output displayed the number of datasets that satisfied the above thresholds among all unique analyses. Upregulated genes in the generated dataset were displayed in a red gradient, decreasing with the top-ranked percentage, while downregulated genes in the generated dataset were presented in a blue gradient that decreased with the gene ranking. A gene could be upregulated in one dataset and downregulated in another depending on the search thresholds and study parameters. It is useful to examine discrepancies in gene expressions among studies, and plots of breast cancer subtypes were conducted with the ggpubr package in R environment as we previously described [43] [44] [45] [46] . In addition to investigating mRNA expressions of PSMC family members in cancer tissues from the Oncomine database, we further searched for their expression levels in cell lines via the CCLE database (https://portals.broadinstitute.org/ccle) [47] . The CCLE is a high-throughput web-based tool with large numbers of human cancer cell lines (n = 1457) and unique datasets (n = 136,488). An RNA sequencing method was selected to search for expressions of PSMC family members in 60 breast cancer cell lines, and results were plotted with default settings as we previously described [48] [49] [50] . To acquire coexpression patterns of PSMC family members in the METABRIC and cBioPortal databases [51] , a Venn diagram was created using InteractiVenn (http://www.interactivenn.net/). Then, 1588 coexpressed genes were further uploaded to Gene Ontology for pathway and network analyses using the MetaCore platform (https://portal.genego.com/), a functional annotation platform for exploring the biological significance behind a large list of genes. Statistical significance as the boundary criterion was set to p<0.05, as we previously described [52] [53] [54] . Together with investigating mRNA expression levels, we concomitantly performed searches for PPI networks of PSMC family members based on coexpressed genes using the STRING database. STRING has protein data comprising 24.6 million proteins in more than 5000 organisms, resulting in more than 2 billion interactions [55] . We selected the k-means clustering algorithm to classify target proteins into different clusters. To understand how mRNA expression levels of PSMC gene family members affected relapse-free survival (RFS) and distant metastasis-free survival (DMFS) of breast cancer patients, we performed a survival analysis using the Kaplan-Meier plotter database (https://kmplot.com/) [56] . Meanwhile, we also assessed the prognostic value of PSMC gene family members for pan-cancer analysis in GEPIA2, which contains RNA sequencing expression data from different types of tumors as well as the normal samples from the TCGA [57, 58] . All default settings inKaplan-Meier were selected for our analysis, namely survival curves, log-rank p values, and hazard ratios with 95% confidence intervals (CIs). To further analyze the infiltration level of immune cells, we applied TIMER 2.0 (http://timer.compgenomics.org) across 31 cancer types comprising of more 10,000 samples [59, 60] . The differences between normal and tumor in mRNA expression of PSMC genes were obtained using DiffExp module. We then selected B cells, T cells clusters including CD4+ and CD8+ together with neutrophils, macrophages, and dendritic cells for our analysis. The oral selective oestrogen receptor degrader (SERD) AZD9496 is comparable to fulvestrant in antagonising ER and circumventing endocrine resistance Challenges in Treating Estrogen Receptor-Positive Breast Cancer Targeting the cyclindependent kinases (CDK) 4/6 in estrogen receptorpositive breast cancers Overcoming acquired resistance to anticancer therapy: focus on the PI3K/AKT/mTOR pathway A phase II study of the histone deacetylase inhibitor vorinostat combined with tamoxifen for the treatment of patients with hormone therapy-resistant breast cancer Targeting BCL-2 with the BH3 mimetic ABT-199 in estrogen receptor-positive breast cancer A Phase Ib Dose-Escalation and Expansion Study of the BCL2 Inhibitor Venetoclax Combined with Tamoxifen in ER and BCL2-Positive Metastatic Breast Cancer Cell cycle-dependent localization of the proteasome to chromatin The 26S proteasome is a multifaceted target for anti-cancer therapies MMP9 modulates the metastatic cascade and immune landscape for breast cancer anti-metastatic therapy PSMB5 plays a dual role in cancer development and immunosuppression A corepressor/coactivator exchange complex required for transcriptional activation by nuclear receptors and other regulated transcription factors PSMC6 promotes osteoblast apoptosis through inhibiting PI3K/AKT signaling pathway activation in ovariectomyinduced osteoporosis mouse model PSMC2 is Upregulated in Pancreatic Cancer and Promotes Cancer Cell Proliferation and Inhibits Apoptosis PSMC2 is upregulated in osteosarcoma and regulates osteosarcoma cell proliferation, apoptosis and migration Inhibition of AZIN2-sv induces neovascularization and improves prognosis after myocardial infarction by blocking ubiquitin-dependent talin1 degradation and activating the Akt pathway Knockdown of PSMC3IP suppresses the proliferation and xenografted tumorigenesis of hepatocellular carcinoma cell Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles Distinct roles and differential expression levels of Wnt5a mRNA isoforms in colorectal cancer cells Increased expression of secreted frizzled related protein 1 (SFRP1) predicts ampullary adenocarcinoma recurrence Microarray technology: beyond transcript profiling and genotype analysis Comprehensive Analysis of Prognostic and Genetic Signatures for General Transcription Factor III (GTF3) in Clinical Colorectal Cancer Patients Using Bioinformatics Approaches Analysis of LAGEs Family Gene Signature and Prognostic Relevance in Breast Cancer Characterizing the Relapse Potential in Different Luminal Subtypes of Breast Cancers with Functional Proteomics Bioinformatics Data Mining Repurposes the JAK2 (Janus Kinase 2) Inhibitor Fedratinib for Treating Pancreatic Ductal Adenocarcinoma by Reversing the KRAS (Kirsten Rat Sarcoma 2 Viral Oncogene Homolog)-Driven Gene Signature Sequential Interferon β-Cisplatin Treatment Enhances the Surface Exposure of Calreticulin in Cancer Cells via an Interferon Regulatory Factor 1-Dependent Manner Sorafenib Inhibits Ribonucleotide Reductase Regulatory Subunit M2 (RRM2) in Hepatocellular Carcinoma Cells Profilin-1 overexpression in MDA-MB-231 breast cancer cells is associated with alterations in proteomics biomarkers of cell proliferation, survival, and motility as revealed by global proteomics analyses Caspaserelated apoptosis genes in gliomas by RNA-seq and bioinformatics analysis Transcription alterations of members of the ubiquitin-proteasome network in prostate carcinoma Selection of Endogenous Control Reference Genes for Studies on Type 1 or Type 2 Endometrial Cancer A novel risk score system of immune genes associated with prognosis in endometrial cancer Spatial control of Shoc2-scaffold-mediated ERK1/2 signaling requires remodeling activity of the ATPase PSMC5 Using Quantitative Seroproteomics to Identify Antibody Biomarkers in Pancreatic Cancer Identification of candidate biomarkers associated with apoptosis in melanosis coli: GNG5, LPAR3, MAPK8, and PSMC6 CRISPR Genome-Wide Screening Identifies Dependence on the Proteasome Subunit PSMC6 for Bortezomib Sensitivity in Multiple Myeloma ONCOMINE: a cancer microarray database and integrated data-mining platform Single-cell RNA sequencing reveals gene expression signatures of breast cancer-associated endothelial cells Argininosuccinate lyase interacts with cyclin A2 in cytoplasm and modulates growth of liver tumor cells Gene signatures and prognostic analyses of the Tob/BTG pituitary tumor-transforming gene (PTTG) family in clinical breast cancer patients The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity Knockdown of serine/threonine-protein kinase 24 promotes tumorigenesis and myeloid-derived suppressor cell expansion in an orthotopic immunocompetent gastric cancer animal model Mutation of the PTCH1 gene predicts recurrence of breast cancer The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data Gene signatures and potential therapeutic targets of amino acid metabolism in estrogen receptorpositive breast cancer CCDC167 as a potential therapeutic target and regulator of cell cycle-related networks in breast cancer Gene signatures of SARS-CoV/SARS-CoV-2-infected ferret lungs in short-and long-term models STRING v11: proteinprotein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells Approaches to treat immune hot, altered and cold tumours with combination immunotherapies Rac1-related 2 Cell cycle_Chromosome condensation in prometaphase 1.24E-06 CAP-G, Cyclin A, CAP-G/G2 Transcription_Role of heterochromatin protein 1 (HP1) family in transcriptional silencing 1 Cell cycle_Regulation of G1/S transition (part 1) 9 TGF-beta receptor type II Rb protein, Cyclin D3, Cyclin D, c-Fos 9 Abnormalities in cell cycle in SCLC 1 Rb protein Prolactin/JAK2 signaling in breast cancer 4 MAD2a, CENP-H 13 Transport_The role of AVP in regulation of Aquaporin 2 and renal water reabsorption 7 Fos 14 Immune response_Function of MEF2 in T lymphocytes 7 MEF2, HDAC5, MEF2C, Calcineurin A (beta) 15 Cell cycle progression in Prostate Cancer 1 Immune response_ETV3 effect on CSF1-promoted macrophage differentiation 1 The authors especially thank Mr. Daniel P. Chamberlin from the Office of Research and Development at Taipei Medical University for professional English editing. Bioinformatics analyses and data mining were conducted by the Bioinformatics Core Facility at Taipei Medical University. The authors declare that they have no conflicts of interest. The study was supported by grants from Taipei Medical University Hospital (107TMU-TMUH-02 to C.C.W. and J.Y.C.) and the Ministry of Science and Technology of Taiwan (110-2636-B-038-004 to J.Y.C). Please browse Full Text version to see the data of Supplementary Tables 3-8 . Table 7 . Pathway analysis of genes coexpressed with proteasome 26S subunit, ATPase 5 (PSMC5) from public breast cancer databases using the MetaCore database (with p<0.01 set as the cutoff value). Supplementary Table 8 . Pathway analysis of genes coexpressed with proteasome 26S subunit, ATPase 6 (PSMC6) from public breast cancer databases using the MetaCore database (with p<0.01 set as the cutoff value).