ABSTRACT
Background: High heterogeneity is an essential feature of malignant tumors. This study aims to reveal the drivers of hepatocellular carcinoma heterogeneity for prognostic stratification and to guide individualized treatment. Methods: Omics data and clinical data for two HCC cohorts were derived from the Cancer Genome Atlas (TCGA) and the International Cancer Genome Atlas (ICGC), respectively. CNV data and methylation data were downloaded from the GSCA database. GSVA was used to estimate the transcriptional activity of KEGG pathways, and consensus clustering was used to categorize the HCC samples. The pRRophetic package was used to predict the sensitivity of samples to anticancer drugs. TIMER, MCPcounter, quanTIseq, and TIDE algorithms were used to assess the components of TME. LASSO and COX analyses were used to establish a prognostic gene signature. The biological role played by genes in HCC cells was confirmed by in vitro experiments. Results: We classified HCC tissues into two categories based on the activity of prognostic pathways. Among them, the transcriptional profile of cluster A HCC is similar to that of normal tissue, dominated by cancer-suppressive metabolic pathways, and has a better prognosis. In contrast, cluster B HCC is dominated by high proliferative activity and has significant genetic heterogeneity. Meanwhile, cluster B HCC is often poorly differentiated, has a high rate of serum AFP positivity, is prone to microvascular invasion, and has shorter overall survival. In addition, we found that mutations, copy number variations, and aberrant methylation were also crucial drivers of the differences in heterogeneity between the two HCC subtypes. Meanwhile, the TME of the two HCC subtypes is also significantly different, which offers the possibility of precision immunotherapy for HCC patients. Finally, based on the prognostic value of molecular subtypes, we developed a gene signature that could accurately predict patients' OS. The riskscore quantified by the signature could evaluate the heterogeneity of HCC and guide clinical treatment. Finally, we confirmed through in vitro experiments that RFPL4B could promote the progression of Huh7 cells. Conclusion: The molecular subtypes we identified effectively exposed the heterogeneity of HCC, which is important for discovering new effective therapeutic targets.
ABSTRACT
BACKGROUND: Peroxisomes are pivotal metabolic organelles that exist in almost all eukaryote cells. A reduction in numbers and enzymatic activities of peroxisomes was found in colon adenocarcinomas. However, the role of peroxisomes or the peroxisome pathway in colorectal cancer (CRC) is not defined. METHODS: In the current study, a peroxisome score was calculated to indicate the activity of the peroxisome pathway using gene set variant analysis based on transcriptomic datasets. CIBERSORTx was chosen to infer enriched immune cells for tumors among subgroups. The SubMap algorithm was applied to predict its sensitivity to immunotherapy. RESULTS: The patients with a relatively low peroxisome score and high level of T-cell immunoglobulin and mucin domain 3 (TIM-3) presented the worse overall survival than others. Moreover, low peroxisome scores were associated with high infiltration of lymphocytes and poor prognosis in those CRC patients. Thus, a PERLowTIM3High CRC risk subpopulation was identified and characterized by high immune infiltration. The results also showed that CD8 T cells and macrophages highly infiltrated tumors of the PERLowTIM3High group, regardless of consortium molecular subtype and microsatellite instability status. This subgroup had the highest tumor mutational burden and overexpression of immune checkpoint genes. Further, the PERLowTIM3High group showed a higher probability of responding to programmed cell death protein-1-based immunotherapy. In addition, genes involved in peroxisomal metabolic processes in CRC were also investigated since peroxisome is a rather pleiotropic and highly metabolic organelle in cell. The results indicated that only those genes involved in fatty acid alpha oxidation could be used to stratify CRC patients as similar as peroxisome pathway genes. CONCLUSIONS: We revealed the favorable prognostic value of the peroxisome pathway in CRC and provided a new CRC stratification based on peroxisomes and TIM3, which might be helpful for CRC diagnostics and personalized treatment.