Precision Oncology Beyond Genomics: The Future Is Here-It Is Just Not Evenly Distributed.

Cancer is a multifactorial disease with increasing incidence. There are more than 100 different cancer types, defined by location, cell of origin, and genomic alterations that influence oncogenesis and therapeutic response. This heterogeneity between tumors of different patients and also the heterogeneity within the same patient's tumor pose an enormous challenge to cancer treatment. In this review, we explore tumor heterogeneity on the longitudinal and the latitudinal axis, reviewing current and future approaches to study this heterogeneity and their potential to support oncologists in tailoring a patient's treatment regimen. We highlight how the ideal of precision oncology is reaching far beyond the knowledge of genetic variants to inform clinical practice and discuss the technologies and strategies already available to improve our understanding and management of heterogeneity in cancer treatment. We will focus on integrating multi-omics technologies with suitable in vitro models and their proficiency in mimicking endogenous tumor heterogeneity.

Inhibiting phosphoglycerate dehydrogenase counteracts chemotherapeutic efficacy against MYCN-amplified neuroblastoma.

Here we sought metabolic alterations specifically associated with MYCN amplification as nodes to indirectly target the MYCN oncogene. Liquid chromatography-mass spectrometry-based proteomics identified seven proteins consistently correlated with MYCN in proteomes from 49 neuroblastoma biopsies and 13 cell lines. Among these was phosphoglycerate dehydrogenase (PHGDH), the rate-limiting enzyme in de novo serine synthesis. MYCN associated with two regions in the PHGDH promoter, supporting transcriptional PHGDH regulation by MYCN. Pulsed stable isotope-resolved metabolomics utilizing 13 C-glucose labeling demonstrated higher de novo serine synthesis in MYCN-amplified cells compared to cells with diploid MYCN. An independence of MYCN-amplified cells from exogenous serine and glycine was demonstrated by serine and glycine starvation, which attenuated nucleotide pools and proliferation only in cells with diploid MYCN but did not diminish these endpoints in MYCN-amplified cells. Proliferation was attenuated in MYCN-amplified cells by CRISPR/Cas9-mediated PHGDH knockout or treatment with PHGDH small molecule inhibitors without affecting cell viability. PHGDH inhibitors administered as single-agent therapy to NOG mice harboring patient-derived MYCN-amplified neuroblastoma xenografts slowed tumor growth. However, combining a PHGDH inhibitor with the standard-of-care chemotherapy drug, cisplatin, revealed antagonism of chemotherapy efficacy in vivo. Emergence of chemotherapy resistance was confirmed in the genetic PHGDH knockout model in vitro. Altogether, PHGDH knockout or inhibition by small molecules consistently slows proliferation, but stops short of killing the cells, which then establish resistance to classical chemotherapy. Although PHGDH inhibition with small molecules has produced encouraging results in other preclinical cancer models, this approach has limited attractiveness for patients with neuroblastoma.

Reflection of neuroblastoma intratumor heterogeneity in the new OHC-NB1 disease model.

Accurate modeling of intratumor heterogeneity presents a bottleneck against drug testing. Flexibility in a preclinical platform is also desirable to support assessment of different endpoints. We established the model system, OHC-NB1, from a bone marrow metastasis from a patient diagnosed with MYCN-amplified neuroblastoma and performed whole-exome sequencing on the source metastasis and the different models and passages during model development (monolayer cell line, 3D spheroid culture and subcutaneous xenograft tumors propagated in mice). OHC-NB1 harbors a MYCN amplification in double minutes, 1p deletion, 17q gain and diploid karyotype, which persisted in all models. A total of 80-540 single-nucleotide variants (SNVs) was detected in each sample, and comparisons between the source metastasis and models identified 34 of 80 somatic SNVs to be propagated in the models. Clonal reconstruction using the combined copy number and SNV data revealed marked clonal heterogeneity in the originating metastasis, with four clones being reflected in the model systems. The set of OHC-NB1 models represents 43% of somatic SNVs and 23% of the cellularity in the originating metastasis with varying clonal compositions, indicating that heterogeneity is partially preserved in our model system.