Glioblastoma cell ability to adapt their functioning to microenvironment changes is a source of the extensive intra-tumor heterogeneity characteristic of this devastating malignant brain tumor. A systemic view of the metabolic pathways underlying glioblastoma cell functioning states is lacking. We analyzed public single cell RNA-sequencing data from glioblastoma surgical resections, which offer the closest available view of tumor cell heterogeneity as encountered at the time of patients’ diagnosis. Unsupervised analyses revealed that information dispersed throughout the cell transcript repertoires encoded the identity of each tumor and masked information related to cell functioning states. Data reduction based on an experimentally-defined signature of transcription factors overcame this hurdle. It allowed cell grouping according to their tumorigenic potential, regardless of their tumor of origin. The approach relevance was validated using independent datasets of glioblastoma cell and tissue transcriptomes, patient-derived cell lines and orthotopic xenografts. Overexpression of genes coding for amino acid and lipid metabolism enzymes involved in anti-oxidative, energetic and cell membrane processes characterized cells with high tumorigenic potential. Modeling of their expression network highlighted the very long chain polyunsaturated fatty acid synthesis pathway at the core of the network. Expression of its most downstream enzymatic component, ELOVL2, was associated with worsened patient survival, and required for cell tumorigenic properties in vivo. Our results demonstrate the power of signature-driven analyses of single cell transcriptomes to obtain an integrated view of metabolic pathways at play within the heterogeneous cell landscape of patient tumors.

Cell motility is critical for tumor malignancy. Metabolism being an obligatory step in shaping cell behavior, we looked for metabolic weaknesses shared by motile cells across the diverse genetic contexts of patients’ glioblastoma. Computational analyses of single-cell transcriptomes from thirty patients’ tumors isolated cells with high motile potential and highlighted their metabolic specificities. These cells were characterized by enhanced mitochondrial load and oxidative stress coupled with mobilization of the cysteine metabolism enzyme 3-Mercaptopyruvate sulfurtransferase (MPST). Functional assays with patients’ tumor-derived cells and -tissue organoids, and genetic and pharmacological manipulations confirmed that the cells depend on enhanced ROS production and MPST activity for their motility. MPST action involved protection of protein cysteine residues from damaging hyperoxidation. Its knockdown translated in reduced tumor burden, and a robust increase in mice survival. Starting from cell-by-cell analyses of the patients’ tumors, our work unravels metabolic dependencies of cell malignancy maintained across heterogeneous genomic landscapes.