Cancer is one of the leading causes of death, and there is great interest in understanding the underlying molecular mechanisms involved in the pathogenesis and progression of individual tumors. We used systems-level approaches to analyze the genome-wide transcriptome of the protein-coding genes of 17 major cancer types with respect to clinical outcome. A general pattern emerged: Shorter patient survival was associated with up-regulation of genes involved in cell growth and with down-regulation of genes involved in cellular differentiation. Using genome-scale metabolic models, we show that cancer patients have widespread metabolic heterogeneity, highlighting the need for precise and personalized medicine for cancer treatment. All data are presented in an interactive open-access database (www.proteinatlas.org/pathology) to allow genome-wide exploration of the impact of individual proteins on clinical outcomes.

Hepatocellular carcinoma (HCC) is one of the most frequent forms of liver cancer, and effective treatment methods are limited due to tumor heterogeneity. There is a great need for comprehensive approaches to stratify HCC patients, gain biological insights into subtypes, and ultimately identify effective therapeutic targets. We stratified HCC patients and characterized each subtype using transcriptomics data, genome-scale metabolic networks and network topology/controllability analysis. This comprehensive systems-level analysis identified three distinct subtypes with substantial differences in metabolic and signaling pathways reflecting at genomic, transcriptomic, and proteomic levels. These subtypes showed large differences in clinical survival associated with altered kynurenine metabolism, WNT/β-catenin–associated lipid metabolism, and PI3K/AKT/mTOR signaling. Integrative analyses indicated that the three subtypes rely on alternative enzymes (e.g., ACSS1/ACSS2/ACSS3, PKM/PKLR, ALDOB/ALDOA, MTHFD1L/MTHFD2/MTHFD1) to catalyze the same reactions. Based on systems-level analysis, we identified 8 to 28 subtype-specific genes with pivotal roles in controlling the metabolic network and predicted that these genes may be targeted for development of treatment strategies for HCC subtypes by performing in silico analysis. To validate our predictions, we performed experiments using HepG2 cells under normoxic and hypoxic conditions and observed opposite expression patterns between genes expressed in high/moderate/low-survival tumor groups in response to hypoxia, reflecting activated hypoxic behavior in patients with poor survival. In conclusion, our analyses showed that the heterogeneous HCC tumors can be stratified using a metabolic network-driven approach, which may also be applied to other cancer types, and this stratification may have clinical implications to drive the development of precision medicine.