ORIGINAL RESEARCH article

Front. Cell Dev. Biol., 14 June 2021

Sec. Molecular and Cellular Oncology

Volume 9 - 2021 | https://doi.org/10.3389/fcell.2021.653528

Comprehensive Analysis Reveals the Potential Regulatory Mechanism Between Ub–Proteasome System and Cell Cycle in Colorectal Cancer

  • Department of General Surgery, Zhongshan Hospital, Fudan University, Shanghai, China

Abstract

The ubiquitin (Ub)–proteasome system (UPS) is an important regulatory component in colorectal cancer (CRC), and the cell cycle is also characterized to play a significant role in CRC. In this present study, we firstly identified UPS-associated differentially expressed genes and all the differentially expressed protein-coding genes in CRC through three differential analyses. UPS-associated genes were also further analyzed via survival analysis. A weighted gene co-expression network analysis (WGCNA) was used to identify the cell cycle-associated genes. We used protein–protein interaction (PPI) network to comprehensively mine the potential mechanism of the UPS–cell cycle regulatory axis. Moreover, we constructed a signature based on UPS-associated genes to predict the overall survival of CRC patients. Our research provides a novel insight view of the UPS and cell cycle system in CRC.

Introduction

The ubiquitin (Ub)–proteasome system (UPS) is identified to regulate the cellular protein by ubiquitination modification. Ub is an important component in the UPS. It consists of a 76-amino-acid protein and is usually highly conserved in eukaryotes (). Among the UPS, three significant enzymes function as enzyme cascade to transmit the Ub to the substrate. The three vital enzymes are identified as the Ub-activating enzyme (E1), the Ub-conjugating enzyme (E2), and the Ub-protein ligase (E3) (). The UPS process can be summarized as follows: at first, E1 activates Ub and transmits it to E2 by an adenosine triphosphate-dependent way. Then, E3 mediates the last step by interacting with E2 that is carried with the Ub and recognizing a specific substrate. E3 are usually identified as the most crucial component among the three-enzyme cascade because the interaction between E3 and the substrate is highly specific, and the ubiquitination of the substrate mainly depends on the E3.

The UPS, especially E3, is also proved to contribute a lot in the process of cancers. Dysregulated E3 were widely reported to occur in diverse cancers. For example, HERC3 was reported to mediate the ubiquitination and the degradation of SMAD7 in glioblastoma (). RNF6 was indicated to induce the progression of colorectal cancer (CRC) through mediating ubiquitination of TLE3 (). Moreover, E3 can regulate the downstream substrate including many oncogenes and tumor suppressors. P53 was reported to be degraded by the E3 ligase RING1 (). Cancer is featured as an uncontrolled cell proliferation that is also regulated by many cell cycle-related regulators (; ). These regulators can also be ubiquitination modified and degraded by E3. Thus, the UPS especially E3 may be the key regulators in the proliferation of cancer cells and may also become the therapeutic targets to address the uncontrolled proliferation of cells ().

Colorectal cancer was ranked the top in the aspect of mortality and morbidity among diverse cancer types (). Given the crucial role of E3 in cancers, it is urgent and necessary to research the UPS in CRC. However, the integrated analysis of UPS or the interaction between UPS and cell cycle-related regulators in CRC is still blank. In this present study, we integrated and analyzed the UPS, especially the E3 ligase in CRC; moreover, we constructed a prognostic signature based on UPS-associated regulators and further depicted the potential interactions between UPS and specific cell cycle-related genes in CRC. We provide a novel insight into the UPS and the latent interaction between cell cycle-associated genes and the UPS in the field of CRC.

Materials and Methods

Acquisition and Processing of the Raw Data

The raw microarray data and relevant clinical information of CRC patients that were based on TCGA (The Cancer Genome Atlas) database and the raw microarray data of normal colon samples that were based on the GTEx (Genotype-Tissue Expression) database were downloaded from XENA1. Then, the data from TCGA and data from GTEx were normalized and combined based on the description of the website. The UPS-associated genes were obtained from an article published by . The detailed information of these UPS-associated genes are provided in Supplementary Table 1.

Identification of Differentially Expressed Genes

An analysis differentially expressed genes was conducted three times according to the diverse grouping of samples. The three comparative groups were set as follows: GTEx normal colon samples VS. TCGA CRC-adjacent normal colon samples; TCGA CRC-adjacent normal colon samples VS. TCGA CRC samples; and GTEx normal colon samples combined with TCGA CRC-adjacent normal colon samples VS. TCGA CRC samples. Wilcoxon test was utilized to perform the differential analysis. The selection criterion was set as FDR < 0.05.

Survival Analysis

Survival analysis was performed on the UPS-associated genes based on the expression pattern and clinical information [overall survival (OS) information] from TCGA. The Kaplan–Meier plot was used to visualize the results and the median expression of the relevant gene was set as the cutoff. The log-rank test was used to analyze the results, and a P-value < 0.05 was identified as statistically significant.

Gene Set Enrichment Analysis and Single-Sample Gene Set Enrichment Analysis

A gene set enrichment analysis (GSEA) was used to identify the biological pathways that might play significant roles in the process of CRC. The input reference gene sets were all the protein-coding genes from the intersection of the three differential analyses. A single-sample gene set enrichment analysis (ssGSEA) was conducted to calculate the score for individual samples based on a specific reference.

Weighted Gene Co-Expression Network Analysis

A weighted gene co-expression network analysis (WGCNA) was performed through the R package WGCNA in R based on the relevant instructions (). Parameters used in the WGCNA process were set as default. Several cell cycle-related pathways that were statistically significant in GSEA were set as references to calculate the correlation between genes and those pathways. A P-value < 0.05 and a correlation value >0.7 were identified as statistically significant. Before WGCNA, we screened out the differentially expressed and OS-related genes. Firstly, 3 times differential analysis same as previously performed was carried out based on the whole genes, the criterion was set as FDR < 0.05. Then the results were subjected to COX analysis. The criterion was set as P < 0.05 and HR < 1.2 or HR < 0.05.

Construction of Protein–Protein Interaction Network

The protein–protein interaction (PPI) network was drawn based on the STRING2. The input genes were differentially expressed UPS-associated genes obtained from the three differential analyses and the cell cycle-related genes that were obtained from the WGCNA. The PPI network was divided into four parts, each module representing E2 and E3 (E3 adaptor and E3 activity, respectively) and the potential cell cycle-associated substrates. The PPI network was visualized by Cytoscape (3.8.2).

Construction and Internal Validation of a Prognostic Signature Based on Differentially Expressed UPS-Associated Genes

The TCGA patients were randomly divided into two parts at a ratio of 7:3. The 70% of the patients were identified as the internal training set and the 30% of the patients were identified as the internal validation set. A least absolute shrinkage and selection operator (LASSO) analysis was conducted to screen out the variates for further analysis. A multivariate Cox regression analysis was performed to construct the signature. The accuracy of the signature was validated in the internal validation set through the receiver operating characteristic (ROC) curve, risk score analysis, and Kaplan–Meier survival analysis. A P-value < 0.05 was identified as statistically significant.

Results

Differentially Expressed UPS-Associated Genes in CRC

The entire flow of the research is presented in Figure 1. We first identified differentially expressed UPS-associated genes through three differential analyses according to the diverse grouping of the entire patients. We first compared the normal colon samples from GTEx with CRC-adjacent normal samples in TCGA. We then compared the CRC-adjacent normal samples with CRC samples in TCGA. Finally, we combined normal colon samples from GTEx with CRC-adjacent normal samples from TCGA and identified the differentially expressed genes between those samples and CRC samples in TCGA. A heatmap is depicted in Figure 2A. Afterward, we intersected the upregulated genes and the downregulated genes from the three differential analyses, respectively. We got 30 upregulated and 30 downregulated UPS-associated genes, respectively (Figure 2B). The expression landscape of these 60 dysregulated UPS-associated genes in three diverse groups is shown in Figure 2B. The results of the differential analysis when comparing normal colon samples in GTEx plus CRC-adjacent normal samples with CRC samples in TCGA are demonstrated in Table 1. Among these 60 UPS-associated genes, 24 genes were identified as E3. We also used violin plots to depict the expression pattern of 24 E3 in Figure 3. The detailed information of the three differential analyses of the 24 E3 is presented in Supplementary Table 2.

TABLE 1

GeneCategoryLogFCP-valueFDRGeneCategoryLogFCP-valueFDR
DYNC1I1E3 adaptor–1.703.18E-991.23E-98DTLE3 adaptor2.101.51E-1441.20E-142
BCL6E3 adaptor–2.331.27E-1025.49E-102WDR43E3 adaptor1.972.88E-1486.27E-146
FBXO17E3 adaptor–1.831.43E-1077.14E-107ENC1E3 adaptor3.312.15E-1511.87E-148
FBXO32E3 adaptor–2.684.24E-1122.55E-111RNF150E3 activity–1.724.11E-1062.00E-105
EML1E3 adaptor–2.143.75E-1122.27E-111PDZRN3E3 activity–2.293.46E-1091.84E-108
DLG4E3 adaptor–1.831.24E-1148.21E-114TRIM2E3 activity1.552.24E-1101.26E-109
WDR86E3 adaptor–1.807.74E-1155.23E-114UBR4E3 activity–2.134.27E-1132.64E-112
FBXO9E3 adaptor–1.717.21E-1195.61E-118RNF128E3 activity3.077.13E-1164.89E-115
DCAF12E3 adaptor1.621.07E-1198.55E-119DTX3LE3 activity1.534.04E-1162.82E-115
SPAG16E3 adaptor–1.615.64E-1224.77E-121CHD3E3 activity–1.582.23E-1161.58E-115
ASB9E3 adaptor1.652.85E-1222.44E-121DTX3E3 activity–2.844.37E-1213.59E-120
DMXL2E3 adaptor–1.796.72E-1235.97E-122PDZRN4E3 activity–2.401.60E-1211.33E-120
SKP1E3 adaptor–2.854.00E-1233.71E-122PHF1E3 activity–1.848.76E-1237.70E-122
CORO6E3 adaptor–3.288.24E-1258.45E-124CHFRE3 activity–1.881.61E-1241.61E-123
ASB2E3 adaptor–2.973.52E-1253.83E-124RNF113AE3 activity1.524.41E-1275.42E-126
WDR34E3 adaptor1.591.04E-1251.17E-124RABGEF1E3 activity–1.927.08E-1288.94E-127
DENND3E3 adaptor–2.532.62E-1273.26E-126TRIM3E3 activity–1.895.51E-1287.05E-127
POC1AE3 adaptor2.271.99E-1292.84E-128HERC1E3 activity–1.995.24E-1286.81E-127
CDC20E3 adaptor3.561.29E-1352.61E-134PDLIM2E3 activity–3.111.10E-1352.29E-134
ZBTB18E3 adaptor1.592.64E-1365.75E-135RNF43E3 activity3.661.07E-1362.44E-135
PAK1IP1E3 adaptor1.763.81E-1379.48E-136PRPF19E3 activity1.734.73E-1381.47E-136
EIF3IE3 adaptor1.682.69E-1377.24E-136DCUN1D5E3 activity1.512.59E-1431.50E-141
SOCS2E3 adaptor–2.062.69E-1377.24E-136UHRF1E3 activity2.097.16E-1444.80E-142
RRP9E3 adaptor2.258.25E-1392.87E-137HERC3E3 activity–1.959.94E-1458.66E-143
WDR4E3 adaptor1.681.03E-1393.90E-138AURKAE3 activity2.922.37E-1462.58E-144
ZBTB33E3 adaptor1.772.90E-1401.20E-138CBX4E3 activity1.981.54E-1462.23E-144
CHAF1BE3 adaptor1.591.06E-1404.60E-139RFWD3E3 activity1.831.45E-1496.32E-147
ZBTB16E3 adaptor–2.556.16E-1413.16E-139UBE2E2E2–1.601.99E-896.14E-89
SKP2E3 adaptor1.893.44E-1431.87E-141UBE2TE22.875.37E-1443.90E-142
CCNFE3 adaptor2.089.83E-1446.12E-142UBE2CE24.271.28E-1483.73E-146

The results of the differential analysis when comparing normal colon samples in GTEx plus CRC-adjacent normal samples with CRC samples in TCGA.

FIGURE 1

FIGURE 2

FIGURE 3

Identification of 12 OS-Related Genes From 60 Dysregulated UPS-Associated Genes

Based on the clinical information (OS) in TCGA, we conducted a survival analysis through Kaplan–Meier plots, and the log-rank test was utilized as the analytical method. We identified 12 genes from the 60 dysregulated UPS-associated genes as OS-related genes (statistically significant) (Figure 4). They were BCL6, CCNF, CHAF1B, CORO6, DLG4, HERC3, PAK1IP1, PHF1, POC1A, RNF113A, TRIM2, and UBE2E2.

FIGURE 4

Identification of Potential Significant Biological Pathways in the Process of CRC

We identified dysregulated protein-coding genes through three differential analyses. The comparative groups were also set as described in the “Materials and Methods” section. We also intersected the upregulated genes and downregulated genes, respectively. Finally, we obtained 3,332 upregulated genes and 3,635 downregulated genes (Figure 5A). GSEA was performed based on these genes. The result of GSEA indicated that cell cycle might be the significant process during the progression of CRC because MYC_TARGETS, E2F_TARGETS, G2M_CHECKPOINT, and P53_PATHWAY were important regulator during the process of cell cycle (Figure 5B). ssGSEA was also conducted, and the cell cycle-related biological pathways obtained from the GSEA were set as reference.

FIGURE 5

Identification of Potential Cell Cycle-Associated Substrates

Before WGCNA, we screened out the differentially expressed and OS-related genes as described before. We performed WGCNA as mentioned in MATERIALS AND METHODS. As depicted in Figure 5C, the suitable soft threshold (power) was set as eight. Through the WGCNA process, we identified six modules, and the cell cycle-related biological pathways obtained from the GSEA based on the all dysregulated protein-coding genes were set as reference. Then, MEbrown was identified as the most correlated module with cell cycle-related pathways. There are a total of 97 genes in MEbrown and they were identified as potential cell cycle-associated substrates (Figures 5D,E).

Identification of the PPI Network

Based on the dysregulated 60 UPS-associated and the potential cell cycle-associated substrates obtained from WGCNA, we conducted the construction of the PPI network. The entire network is shown in Figure 6, and the details of the network are demonstrated in Supplementary Table 3. We divided the entire network into four modules as dysregulated E2, dysregulated E3 adaptor, dysregulated E3 activity, and potential cell cycle-associated substrates. Detailed information of genes involved in the PPI is shown in Supplementary Table 4.

FIGURE 6

Construction of the Prognostic Signature Based on 60 Dysregulated UPS-Associated Genes

We randomly divided the entire TCGA CRC patients into two parts at a ratio of 7:3; moreover, LASSO and multivariate cox regression analyses were performed to construct the signature based on the 70% of patients. The results of the univariate Cox analysis based on 60 differentially expressed UPS-associated genes are presented in Table 2. The efficiency was evaluated through a risk score analysis, Kaplan–Meier survival analysis, and ROC. The efficiency was also validated in the 30% of the patients through equivalent methods (Figures 7A–C). The formula was risk score = exp of CCNF × (−0.545) + exp of ZBTB18 × (−0.502) + exp of RNF113A × 0.687 + exp of CHD3 × 0.347 + exp of PHF1 × 0.401 + exp of BCL6 × 0.368 + exp of ZBTB16 × (−1.216) + exp of SKP1 × (−0.588). Then, the cutoff between high risk and low risk was 1.282. The areas under the curve (AUC) of ROC at 1, 3, and 5 years in the internal training set were 0.715, 0.71, and 0.805, respectively. The AUC of ROC at 1, 3, and 5 years in the internal test set were 0.715, 0.621, and 0.655, respectively. In the entire TCGA dataset, the AUC of ROC at 1, 3, and 5 years were 0.713, 0.677, and 0.757, respectively.

TABLE 2

GeneHROR (95% CI high)OR (95% CI low)P-valueGeneHROR (95% CI high)OR (95% CI low)P-value
PHF11.811.322.470.00SKP20.840.621.140.27
DLG41.641.172.300.00ASB90.910.751.090.30
CCNF0.650.480.880.01DTL0.850.621.170.32
BCL61.411.101.820.01DENND31.240.791.930.34
WDR861.931.153.250.01FBXO171.120.881.420.37
TRIM20.740.580.940.01UBE2T0.890.691.160.40
CHD31.421.071.880.01PDZRN31.090.871.360.47
CORO62.541.175.490.02DTX3L0.890.631.260.52
ZBTB180.720.540.950.02FBXO90.840.501.410.52
UBE2E21.291.021.640.03ASB21.080.831.410.55
RNF430.850.730.990.03PDLIM21.120.781.600.55
DCAF120.580.350.960.03EML11.110.791.550.55
DTX31.291.011.650.04UBE2C0.940.761.170.59
CHAF1B0.720.521.000.05RNF1280.950.781.160.61
POC1A0.780.611.000.05ENC10.940.711.240.65
RNF113A1.360.991.860.06SPAG160.940.711.240.67
AURKA0.790.601.040.09RRP90.950.721.250.71
WDR430.730.501.080.11CHFR1.060.771.450.73
SKP10.730.491.090.13CBX41.050.791.400.73
DCUN1D50.760.541.080.13UBR40.930.631.390.73
PAK1IP10.780.551.090.14SOCS20.950.681.340.78
RFWD30.720.461.130.16HERC10.940.611.450.79
PDZRN41.280.891.830.18WDR341.030.821.300.79
DYNC1I11.260.891.780.19PRPF191.040.711.540.84
CDC200.870.711.070.19FBXO321.020.831.250.85
HERC30.760.501.150.19TRIM31.030.721.490.85
WDR40.820.611.110.20DMXL21.030.721.470.88
RNF1501.280.831.980.26RABGEF11.040.571.880.90
EIF3I0.840.621.140.26ZBTB330.990.721.360.95
UHRF10.860.661.120.27ZBTB161.000.541.840.99

The results of univariate Cox analysis based on 60 differentially expressed UPS-associated genes.

FIGURE 7

Discussion

Colorectal cancer is identified as one digestive system cancer type that has high morbidity and mortality (). Although the treatment of CRC progresses a lot during recent years, the potential regulatory mechanism of CRC is not clearly elucidated yet. The cell cycle is a crucial component in biological development. The cell cycle can be summarized into four phases: in the S phase, DNA replication is frequent; in the M phase, a single cell can divide into two daughter cells. There are also two gap phases between S and M, which can be characterized as G1 and G2. It is widely recognized that G1 is a phase in which a cell can be sensitive to stimulation of growth. G2 is the phase after S and specific for cell entering mitosis. Moreover, cancer is characterized as uncontrolled cell proliferation (; ). Uncontrolled cancer cell proliferation is also commonly observed in CRC. found that miR-1258 could regulate the CRC cell proliferation via regulating the cell cycle. Moreover, there was also a prognostic signature based on cell cycle specific in colon cancer (). These studies implied that the cell cycle also has a significant function in the process of CRC. However, the detailed mechanism of the regulation of the cell cycle in CRC is not clear until now, and according to the vital role of the cell cycle in CRC, it is urgent to demonstrate the potential regulatory mechanism of the cell cycle in CRC.

The Ub–UPS is proved to be involved in the regulation of the process of many cancers. It can be summarized into many regulatory mechanisms. Firstly, the proteins that are involved in the UPS can be oncogenes or tumor suppressors. For instance, HERC3 was reported to be a tumor suppressor via regulating SMAD7 in glioblastoma (). RNF6 was reported to induce the progression of CRC through mediating ubiquitination of TLE3 (). Secondly, many oncogenes or tumor suppressors can be regulated by E3, and P53 was indicated to be degraded by the E3 ligase RING1 (). Given the important role of UPS in cancer, it is also urgent to comprehensively analyze the regulatory mechanism of UPS in CRC.

Given the important role of cell cycle and UPS in CRC, we comprehensively analyzed differentially expressed UPS-associated genes in CRC through three differential analyses. Moreover, we discovered that the cell cycle is one of the most important biological processes in the progression of CRC. We also used WGCNA to identify some cell cycle-associated genes that are specific in CRC. Furthermore, we used the differentially expressed UPS-associated genes and cell cycle-associated genes to construct the PPI network. Afterward, we used UPS-associated genes to construct the OS prognostic signature in CRC with relative considerable AUC.

Among the UPS, E3 ligase is undoubtedly the most important component because the interaction between E3 and the substrate is highly specific, and the ubiquitination of the substrate mainly depends on the E3. In this research, we found a total of 60 differentially expressed UPS-associated genes, among them, 24 genes were identified as E3 ligases. Among these 24 differentially expressed E3 ligases, many were proved to have relevant regulatory roles in the process of CRC. For example, AURKA was upregulated by ARID3 in CRC (). CBX4 was reported to involve in the process of long non-coding RNA RAMS11 regulating the metastasis of CRC (). ASB8 was reported to be controlled by miR-452 in CRC cells (). FBXO45 could be regulated by RP11 through the Siah1-Fbxo45/Zeb1 axis (). TRIM27 was also reported to be an oncogene in CRC (). Previous researches also confirmed that our analysis was reliable. Combing the results of survival analysis and the results of differential analysis, we found that HERC3 is an E3 ligase that owns the same trend of survival analysis and differential analysis, indicating the potential research value of HERC3 in CRC.

Through three differential analyses, we identified 3,323 upregulated genes and 3,635 downregulated genes in CRC. GSEA indicated that the cell cycle is an important component in CRC. Via WGCNA, we also identified cell cycle-associated genes specific in CRC. Moreover, the PPI network based on UPS-associated genes and cell cycle-associated genes provided many latent research orientations for the mechanism of uncontrolled cell proliferation in CRC.

Finally, we constructed an OS-associated signature based on the 60 differentially expressed UPS-associated genes. However, the signature still lacks enough validation, and this limitation is also a novel research direction for us. Compared with other signatures in CRC, our signature can predict the prognosis in colon cancer and rectal cancer with a considerate AUC; other signatures mainly focus on one cancer type only (colon cancer or rectal cancer) (,,,, ). Among these genes involved in the signature, some of them were previously reported to perform latent jobs in the progression of cancer. ZBTB18 was reported to be upregulated by circTP63 and further promote hepatocellular carcinoma progression (). RNF113A was revealed to promote the proliferation, migration, and invasion in esophageal squamous cell carcinoma (). SKP1 was demonstrated to be involved in an axis to promote bladder cancer proliferation and is controlled by circGLIS3 ().

In conclusion, we used bioinformatic analysis to reveal the potential regulatory mechanism between UPS-associated genes and potential cell cycle-related substrates specific in CRC. Besides, we constructed a prognostic signature based on the UPS-associated genes. Our research provides a novel insight of the UPS and cell cycle system in CRC.

Statements

Data availability statement

The raw microarray data and relevant clinical information of CRC patients that were based on TCGA (The Cancer Genome Atlas) database, and raw microarray data of normal colon samples that were based on GTEx (The Genotype-Tissue Expression) database were downloaded from XENA (http://xena.ucsc.edu/).

Ethics statement

The studies involving human participants were reviewed and approved by the data was obtained from TCGA and GTEx database. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

ZZ, JX, and LR designed and conducted the study. ZZ and QF wrote the manuscript. JC and WT helped to improve and design the study. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.653528/full#supplementary-material

Supplementary Table 1

Detailed information of UPS-associated genes.

Supplementary Table 2

Detailed information of three differential analyses of the 24 E3.

Supplementary Table 3

Detailed interactions in the PPI network.

Supplementary Table 4

Detailed information of genes involved in the PPI.

References

Summary

Keywords

Ub–proteasome system, cell cycle, signature, colorectal cancer, bioinformatic

Citation

Zhang Z, Chen J, Tang W, Feng Q, Xu J and Ren L (2021) Comprehensive Analysis Reveals the Potential Regulatory Mechanism Between Ub–Proteasome System and Cell Cycle in Colorectal Cancer. Front. Cell Dev. Biol. 9:653528. doi: 10.3389/fcell.2021.653528

Received

14 January 2021

Accepted

18 March 2021

Published

14 June 2021

Volume

9 - 2021

Edited by

Mingyan Zhu, Affiliated Hospital of Nantong University, China

Reviewed by

Hao Qian, Shanghai Jiao Tong University, China; Yuhua Shi, Yancheng Third People’s Hospital, China

Updates

Copyright

*Correspondence: Qingyang Feng, Jianmin Xu, Li Ren,

These authors have contributed equally to this work

This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Cell and Developmental Biology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics