The Clonal Evolution of Metastatic Osteosarcoma as Shaped by Cisplatin Treatment

SJPedPanel: A pan-cancer gene panel for childhood malignancies

Background

Large scale genomics projects have identified driver alterations for most childhood cancers that provide reliable biomarkers for clinical diagnosis and disease monitoring using targeted sequencing. However, there is lack of a comprehensive panel that matches the list of known driver genes. Here we fill this gap by developing SJPedPanel for childhood cancers.

Results

SJPedPanel covers 5,275 coding exons of 357 driver genes, 297 introns frequently involved in rearrangements that generate fusion oncoproteins, commonly amplified/deleted regions (e.g., MYCN for neuroblastoma, CDKN2A and PAX5 for B-/T-ALL, and SMARCB1 for AT/RT), and 7,590 polymorphism sites for interrogating tumors with aneuploidy, such as hyperdiploid and hypodiploid B-ALL or 17q gain neuroblastoma. We used driver alterations reported from an established real-time clinical genomics cohort (n=253) to validate this gene panel. Among the 485 pathogenic variants reported, our panel covered 417 variants (86%). For 90 rearrangements responsible for oncogenic fusions, our panel covered 74 events (82%). We re-sequenced 113 previously characterized clinical specimens at an average depth of 2,500X using SJPedPanel and recovered 354 (91%) of the 389 reported pathogenic variants. We then investigated the power of this panel in detecting mutations from specimens with low tumor purity (as low as 0.1%) using cell line-based dilution experiments and discovered that this gene panel enabled us to detect ∼80% variants with allele fraction of 0.2%, while the detection rate decreases to ∼50% when the allele fraction is 0.1%. We finally demonstrate its utility in disease monitoring on clinical specimens collected from AML patients in morphologic remission.

Conclusions

SJPedPanel enables the detection of clinically relevant genetic alterations including rearrangements responsible for subtype-defining fusions for childhood cancers by targeted sequencing of ∼0.15% of human genome. It will enhance the analysis of specimens with low tumor burdens for cancer monitoring and early detection.

Does primary tumor resection improve survival for patients with sarcomas of pelvic bones, sacrum, and coccyx who have metastasis at diagnosis ?

BACKGROUND

Recent studies demonstrated that primary tumor resection (PTR) improves survival of patients with metastatic bone sarcomas. However, it remains quite unclear regarding the role of PTR in the treatment of sarcomas of pelvic bones with synchronous metastasis at diagnosis.

METHODS

Using the Surveillance, Epidemiology, and End Results Program, we enrolled a total of 385 patients with sarcomas of pelvic bones, sacrum, and coccyx who have metastasis at initial diagnosis, including 139 patients with osteosarcoma, 176 with Ewing sarcoma, and 70 with chondrosarcoma. Association between PTR and disease-specific survival (DSS) were investigated using the univariable and multivariable Cox regression models. Hazard ratio (HR) and 95% confidence interval (CI) were reported. Representative institutional PTR strategies and clinical outcomes for patients with metastatic pelvic sarcomas from our cancer center were displayed.

RESULTS

The usage rate of PTR was 28.1% (39/139) in osteosarcoma, 13.6% (24/176) in Ewing sarcoma, and 41.4% (29/70) in chondrosarcoma with synchronous metastatic lesions. PTR was not associated with an improved DSS for metastatic pelvic osteosarcoma (HR = 0.686, 95% CI = 0.430 ~ 1.094, P = 0.113) and Ewing sarcoma (HR = 0.580, 95% CI = 0.291 ~ 1.154, P = 0.121). The use of PTR was associated with an improved DSS for metastatic pelvic chondrosarcoma (HR = 0.464, 95% CI = 0.225 ~ 0.954, P = 0.037).

CONCLUSION

Primary lesion resection may provide a survival benefit for metastatic chondrosarcoma, but not for osteosarcoma and Ewing sarcoma of pelvic bones, sacrum, and coccyx. This population-based study recommends an active surgical intervention for metastatic chondrosarcoma while non-surgical treatment for metastatic osteosarcoma and Ewing sarcoma of the pelvis in terms of survival improvement.

Modeling phenotypic heterogeneity towards evolutionarily inspired osteosarcoma therapy

Osteosarcoma is the most common bone sarcoma in children and young adults. While universally delivered, chemotherapy only benefits roughly half of patients with localized disease. Increasingly, intratumoral heterogeneity is recognized as a source of therapeutic resistance. In this study, we develop and evaluate an in vitro model of osteosarcoma heterogeneity based on phenotype and genotype. Cancer cell populations vary in their environment-specific growth rates and in their sensitivity to chemotherapy. We present the genotypic and phenotypic characterization of an osteosarcoma cell line panel with a focus on co-cultures of the most phenotypically divergent cell lines, 143B and SAOS2. Modest environmental (pH, glutamine) or chemical perturbations dramatically shift the success and composition of cell lines. We demonstrate that in nutrient rich culture conditions 143B outcompetes SAOS2. But, under nutrient deprivation or conventional chemotherapy, SAOS2 growth can be favored in spheroids. Importantly, when the simplest heterogeneity state is evaluated, a two-cell line coculture, perturbations that affect the faster growing cell line have only a modest effect on final spheroid size. Thus the only evaluated therapies to eliminate the spheroids were by switching therapies from a first strike to a second strike. This extensively characterized, widely available system, can be modeled and scaled to allow for improved strategies to anticipate resistance in osteosarcoma due to heterogeneity.

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease.

SIGNIFICANCE

The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

Transcriptional profiling of canine osteosarcoma identifies prognostic gene expression signatures with translational value for humans

Canine osteosarcoma is increasingly recognized as an informative model for human osteosarcoma. Here we show in one of the largest clinically annotated canine osteosarcoma transcriptional datasets that two previously reported, as well as de novo gene signatures devised through single sample Gene Set Enrichment Analysis (ssGSEA), have prognostic utility in both human and canine patients. Shared molecular pathway alterations are seen in immune cell signaling and activation including TH1 and TH2 signaling, interferon signaling, and inflammatory responses. Virtual cell sorting to estimate immune cell populations within canine and human tumors showed similar trends, predominantly for macrophages and CD8+ T cells. Immunohistochemical staining verified the increased presence of immune cells in tumors exhibiting immune gene enrichment. Collectively these findings further validate naturally occurring osteosarcoma of the pet dog as a translationally relevant patient model for humans and improve our understanding of the immunologic and genomic landscape of the disease in both species.

Subclonal somatic copy number alterations emerge and dominate in recurrent osteosarcoma

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma.

Exploiting active nuclear import for efficient delivery of Auger electron emitters into the cell nucleus

BACKGROUND

The most attractive features of Auger electrons (AEs) in cancer therapy are their extremely short range and sufficiently high linear energy transfer (LET) for a majority of them. The cytotoxic effects of AE emitters can be realized only in close vicinity to sensitive cellular targets and they are negligible if the emitters are located outside the cell. The nucleus is considered the compartment most sensitive to high LET particles. Therefore, the use of AE emitters could be most useful in specific recognition of a cancer cell and delivery of AE emitters into its nucleus.

PURPOSE

This review describes the studies aimed at developing effective anticancer agents for the delivery of AE emitters to the nuclei of target cancer cells. The use of peptide-based conjugates, nanoparticles, recombinant proteins, and other constructs for AE emitter targeted intranuclear delivery as well as their advantages and limitations are discussed.

CONCLUSION

Transport from the cytoplasm to the nucleus along with binding to the cancer cell is one of the key stages in the delivery of AE emitters; therefore, several constructs for exploitation of this transport have been developed. The transport is carried out through a nuclear pore complex (NPC) with the use of specific amino acid nuclear localization sequences (NLS) and carrier proteins named importins, which are located in the cytosol. Therefore, the effectiveness of NLS-containing delivery constructs designed to provide energy-dependent transport of AE emitter into the nuclei of cancer cells also depends on their efficient entry into the cytosol of the target cell.

How Genetics and Genomics Advances Are Rewriting Pediatric Cancer Research and Clinical Care

In the last two decades, thanks to the data that have been obtained from the Human Genome Project and the development of next-generation sequencing (NGS) technologies, research in oncology has produced extremely important results in understanding the genomic landscape of pediatric cancers, which are the main cause of death during childhood. NGS has provided significant advances in medicine by detecting germline and somatic driver variants that determine the development and progression of many types of cancers, allowing a distinction between hereditary and non-hereditary cancers, characterizing resistance mechanisms that are also related to alterations of the epigenetic apparatus, and quantifying the mutational burden of tumor cells. A combined approach of next-generation technologies allows us to investigate the numerous molecular features of the cancer cell and the effects of the environment on it, discovering and following the path of personalized therapy to defeat an "ancient" disease that has had victories and defeats. In this paper, we provide an overview of the results that have been obtained in the last decade from genomic studies that were carried out on pediatric cancer and their contribution to the more accurate and faster diagnosis in the stratification of patients and the development of new precision therapies.

Elevated Mutational Age in Blood of Children Treated for Cancer Contributes to Therapy-Related Myeloid Neoplasms

Childhood cancer survivors are confronted with various chronic health conditions like therapy-related malignancies. However, it is unclear how exposure to chemotherapy contributes to the mutation burden and clonal composition of healthy tissues early in life. Here, we studied mutation accumulation in hematopoietic stem and progenitor cells (HSPC) before and after cancer treatment of 24 children. Of these children, 19 developed therapy-related myeloid neoplasms (t-MN). Posttreatment HSPCs had an average mutation burden increase comparable to what treatment-naïve cells accumulate during 16 years of life, with excesses up to 80 years. In most children, these additional mutations were induced by clock-like processes, which are also active during healthy aging. Other patients harbored mutations that could be directly attributed to treatments like platinum-based drugs and thiopurines. Using phylogenetic inference, we demonstrate that most t-MN in children originate after the start of treatment and that leukemic clones become dominant during or directly after chemotherapy exposure.

SIGNIFICANCE

Our study shows that chemotherapy increases the mutation burden of normal blood cells in cancer survivors. Only few drugs damage the DNA directly, whereas in most patients, chemotherapy-induced mutations are caused by processes similar to those present during normal aging. This article is highlighted in the In This Issue feature, p. 1825.

Therapeutic and prognostic insights from the analysis of cancer mutational signatures

The somatic mutations in each cancer genome are caused by multiple mutational processes, each of which leaves a characteristic imprint (or 'signature'), potentially caused by specific etiologies or exposures. Deconvolution of these signatures offers a glimpse into the evolutionary history of individual tumors. Recent work has shown that mutational signatures may also yield therapeutic and prognostic insights, including the identification of cell-intrinsic signatures as biomarkers of drug response and prognosis. For example, mutational signatures indicating homologous recombination deficiency are associated with poly(ADP)-ribose polymerase (PARP) inhibitor sensitivity, whereas APOBEC-associated signatures are associated with ataxia telangiectasia and Rad3-related kinase (ATR) inhibitor sensitivity. Furthermore, therapy-induced mutational signatures implicated in cancer progression have also been uncovered, including the identification of thiopurine-induced TP53 mutations in leukemia. In this review, we explore the various ways mutational signatures can reveal new therapeutic and prognostic insights, thus extending their traditional role in identifying disease etiology.

Therapeutic vulnerability to PARP1,2 inhibition in RB1-mutant osteosarcoma

Loss-of-function mutations in the RB1 tumour suppressor are key drivers in cancer, including osteosarcoma. RB1 loss-of-function compromises genome-maintenance and hence could yield vulnerability to therapeutics targeting such processes. Here we demonstrate selective hypersensitivity to clinically-approved inhibitors of Poly-ADP-Polymerase1,2 inhibitors (PARPi) in RB1-defective cancer cells, including an extended panel of osteosarcoma-derived lines. PARPi treatment results in extensive cell death in RB1-defective backgrounds and prolongs survival of mice carrying human RB1-defective osteosarcoma grafts. PARPi sensitivity is not associated with canonical homologous recombination defect (HRd) signatures that predict PARPi sensitivity in cancers with BRCA1,2 loss, but is accompanied by rapid activation of DNA replication checkpoint signalling, and active DNA replication is a prerequisite for sensitivity. Importantly, sensitivity in backgrounds with natural or engineered RB1 loss surpasses that seen in BRCA-mutated backgrounds where PARPi have established clinical benefit. Our work provides evidence that PARPi sensitivity extends beyond cancers identifiable by HRd and advocates PARP1,2 inhibition as a personalised strategy for RB1-mutated osteosarcoma and other cancers.

The phosphoinositide 3-kinase inhibitor ZSTK474 increases the susceptibility of osteosarcoma cells to oncolytic vesicular stomatitis virus VSVΔ51 via aggravating endoplasmic reticulum stress

Blockage of phosphoinositide 3-kinase (PI3K)/protein kinase B (Akt) signal pathway is effective to increase the cytotoxic effects of oncolytic virus on cancer cells, but the detailed mechanisms are still largely unknown. Based on this, the present study managed to investigate the anti-tumor effects of PI3K inhibitor ZSTK474 and oncolytic vesicular stomatitis virus VSVΔ51 combination treatments on osteosarcoma (OS) in vitro and in vivo. Specifically, ZSTK474 aggravated the inhibiting effects of VSVΔ51 on osteosarcoma development by triggering endoplasmic reticulum (ER)-stress mediated apoptotic cell death. Mechanistically, either ZSTK474 or VSVΔ51 alone had limited effects on cell viability in osteosarcoma cells, while ZSTK474 and VSVΔ51 combination treatments significantly induced osteosarcoma cell apoptosis. Interestingly, VSVΔ51 increased the expression levels of IRE1α and p-PERK to initiate ER stress in osteosarcoma cells, which were aggravated by co-treating cells with ZSTK474. Next, the promoting effects of ZSTK474-VSVΔ51 combined treatment on osteosarcoma cell death were abrogated by the ER-stress inhibitor 4-phenyl butyric acid (4-PBA), indicating that ZSTK474 enhanced the cytotoxic effects of VSVΔ51 on osteosarcoma cells in an ER-stress dependent manner. Finally, the xenograft tumor-bearing mice models were established, and the results showed that ZSTK474-VSVΔ51 combined treatment synergistically hindered tumorigenesis of osteosarcoma cells in vivo. Taken together, our data suggested that ZSTK474 was a novel agent to enhance the cytotoxic effects of VSVΔ51 on osteosarcoma by aggravating ER-stress, and the present study might provide alternative therapy treatments for osteosarcoma in clinic.

The developmental origin of cancers defines basic principles of cisplatin resistance

Cisplatin-based chemotherapy has been used for more than four decades as a standard therapeutic option in several tumor entities. However, being a multifaceted and heterogeneous phenomenon, inherent or acquired resistance to cisplatin remains a major obstacle during the treatment of several solid malignancies and inevitably results in disease progression. Hence, we felt there was an urgent need to evaluate common mechanisms between multifarious cancer entities to identify patient-specific therapeutic strategies. We found joint molecular and (epi)genetic resistance mechanisms and specific cisplatin-induced mutational signatures that depended on the developmental origin (endo-, meso-, ectoderm) of the tumor tissue. Based on the findings of thirteen tumor entities, we identified three resistance groups, where Group 1 (endodermal origin) prominently indicates NRF2-pathway activation, Group 2 (mesodermal origin, primordial germ cells) shares elevated DNA repair mechanisms and decreased apoptosis induction, and Group 3 (ectodermal and paraxial mesodermal origin) commonly presents deregulated apoptosis induction and alternating pathways as the main cisplatin-induced resistance mechanisms. This review further proposes potential and novel therapeutic strategies to improve the outcome of cisplatin-based chemotherapy.

Combination Therapy of Cisplatin and other Agents for Osteosarcoma: A Review

Background:

Osteosarcoma is the most common type of primary bone tumor in children and adolescents, which is associated with rapid progression and poor prognosis. Multimodal therapy is the most common approach utilized for osteosarcoma management, such as the application of chemotherapy in combination with surgery or radiation therapy. Cisplatin is one of the predominantly used chemotherapeutic agents for osteosarcoma. Optimally, it is employed in combination with other chemotherapeutic drugs along with surgery or radiation therapy. Despite the availability of numerous treatment approaches, the patient survival rate has not definitively improved over the past three decades.

Methods:

We have summarized all findings regarding the combination of cisplatin with other chemotherapeutic agents as well as with phytochemical compounds.

Results:

A combination of cisplatin with a phytochemical compound synergistically enhances the killing effect of cisplatin on osteosarcoma cells with fewer side effects compared to combination with other chemotherapeutic agents.

Conclusion:

Conclusively, a combination of cisplatin with selected chemotherapeutic drugs has been shown to be effective. However, the unchanged survival rate has posed an urge to search for a new combination regimen. As a collaborative effort to substantiate the therapeutic efficacy, the combination with phytochemical compounds shows a promising response both in vitro as well as in the preclinical study.

Data-driven approaches to advance research and clinical care for pediatric cancer

Pediatric cancer is a rare disease with a distinct etiology and mutational landscape compared with adult cancer. Multi-omics profiling of retrospective and prospective cohorts coupled with innovative computational analysis have been instrumental in uncovering mechanisms of tumorigenesis and drug resistance that are now informing pediatric cancer clinical therapy. In this review we present the major data resources of pediatric cancer and actionable insights into pediatric cancer etiology stemming from the identification of oncogenic gene fusions, mutational signature analysis, systems biology, cancer predisposition and survivorship studies - that have led to improved clinical diagnosis, discovery of new drug-targets, pharmacological therapy, and screening for genetic predisposition. Ultimately, integration of large-scale omics datasets generated through international collaboration is required to maximize the power of data-driven approaches to advance pediatric cancer research informing clinical therapy.

A comparative analysis of the mutagenicity of platinum-containing chemotherapeutic agents reveals direct and indirect mutagenic mechanisms

Platinum-based drugs are a mainstay of cancer chemotherapy. However, their mutagenic effect can increase tumour heterogeneity, contribute to the evolution of treatment resistance and also induce secondary malignancies. We coupled whole genome sequencing with phenotypic investigations on two cell line models to compare the magnitude and examine the mechanism of mutagenicity of cisplatin, carboplatin and oxaliplatin. Cisplatin induced significantly more base substitution mutations than carboplatin or oxaliplatin when used at equitoxic concentrations on human TK6 or chicken DT40 cells, and also induced the highest number of short insertions and deletions. The analysis of base substitution spectra revealed that all three tested platinum drugs elicit both a direct mutagenic effect at purine dinucleotides, and an indirect effect of accelerating endogenous mutagenic processes, whereas the direct mutagenic effect appeared to correlate with the level of DNA damage caused as assessed through histone H2AX phosphorylation and single-cell agarose gel electrophoresis, the indirect mutagenic effects were equal. The different mutagenicity and DNA-damaging effect of equitoxic platinum drug treatments suggest that DNA damage independent mechanisms significantly contribute to their cytotoxicity. Thus, the comparatively high mutagenicity of cisplatin should be taken into account in the design of chemotherapeutic regimens.

SequencErr: measuring and suppressing sequencer errors in next-generation sequencing data

Background

There is currently no method to precisely measure the errors that occur in the sequencing instrument/sequencer, which is critical for next-generation sequencing applications aimed at discovering the genetic makeup of heterogeneous cellular populations.

Results

We propose a novel computational method, SequencErr, to address this challenge by measuring the base correspondence between overlapping regions in forward and reverse reads. An analysis of 3777 public datasets from 75 research institutions in 18 countries revealed the sequencer error rate to be ~ 10 per million (pm) and 1.4% of sequencers and 2.7% of flow cells have error rates > 100 pm. At the flow cell level, error rates are elevated in the bottom surfaces and > 90% of HiSeq and NovaSeq flow cells have at least one outlier error-prone tile. By sequencing a common DNA library on different sequencers, we demonstrate that sequencers with high error rates have reduced overall sequencing accuracy, and removal of outlier error-prone tiles improves sequencing accuracy. We demonstrate that SequencErr can reveal novel insights relative to the popular quality control method FastQC and achieve a 10-fold lower error rate than popular error correction methods including Lighter and Musket.

Conclusions

Our study reveals novel insights into the nature of DNA sequencing errors incurred on DNA sequencers. Our method can be used to assess, calibrate, and monitor sequencer accuracy, and to computationally suppress sequencer errors in existing datasets.

St. Jude Cloud: A Pediatric Cancer Genomic Data-Sharing Ecosystem

Effective data sharing is key to accelerating research to improve diagnostic precision, treatment efficacy, and long-term survival in pediatric cancer and other childhood catastrophic diseases. We present St. Jude Cloud (https://www.stjude.cloud), a cloud-based data-sharing ecosystem for accessing, analyzing, and visualizing genomic data from >10,000 pediatric patients with cancer and long-term survivors, and >800 pediatric sickle cell patients. Harmonized genomic data totaling 1.25 petabytes are freely available, including 12,104 whole genomes, 7,697 whole exomes, and 2,202 transcriptomes. The resource is expanding rapidly, with regular data uploads from St. Jude's prospective clinical genomics programs. Three interconnected apps within the ecosystem-Genomics Platform, Pediatric Cancer Knowledgebase, and Visualization Community-enable simultaneously performing advanced data analysis in the cloud and enhancing the Pediatric Cancer knowledgebase. We demonstrate the value of the ecosystem through use cases that classify 135 pediatric cancer subtypes by gene expression profiling and map mutational signatures across 35 pediatric cancer subtypes. SIGNIFICANCE: To advance research and treatment of pediatric cancer, we developed St. Jude Cloud, a data-sharing ecosystem for accessing >1.2 petabytes of raw genomic data from >10,000 pediatric patients and survivors, innovative analysis workflows, integrative multiomics visualizations, and a knowledgebase of published data contributed by the global pediatric cancer community.This article is highlighted in the In This Issue feature, p. 995.

Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations

Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes. To better understand neuroblastoma pathogenesis, here we analyze whole-genome, whole-exome and/or transcriptome data from 702 neuroblastoma samples. Forty percent of samples harbor at least one recurrent driver gene alteration and most aberrations, including MYCN, ATRX, and TERT alterations, differ in frequency by age. MYCN alterations occur at median 2.3 years of age, TERT at 3.8 years, and ATRX at 5.6 years. COSMIC mutational signature 18, previously associated with reactive oxygen species, is the most common cause of driver point mutations in neuroblastoma, including most ALK and Ras-activating variants. Signature 18 appears early and is continuous throughout disease evolution. Signature 18 is enriched in neuroblastomas with MYCN amplification, 17q gain, and increased expression of mitochondrial ribosome and electron transport-associated genes. Recurrent FGFR1 variants in six patients, and ALK N-terminal structural alterations in five samples, identify additional patients potentially amenable to precision therapy.

Pharmacogenomics for the efficacy of platinum-based chemotherapy: Old drugs, new integrated perspective

Platinum-based chemotherapy remains the cornerstone of treatment for many malignancies. However, although therapeutic efficiency varies greatly among individuals, there is a lack of pharmacogenomic biomarkers that can be used in clinical settings to identify chemosensitive patients and allow stratification. With the development of high-throughput screening techniques and systems biology approaches, a growing body of evidence has shown that platinum resistance is a multifactorial, multi-dimensional, dynamic process incorporating genetic background, tumor evolution and gut microbes. This review critically summarizes potential pharmacogenomic biomarkers for predicting the efficacy of platinum drugs and provides a comprehensive, time-varying perspective that integrates multiple markers.

Extensive Clonal Branching Shapes the Evolutionary History of High-Risk Pediatric Cancers

Darwinian evolution of tumor cells remains underexplored in childhood cancer. We here reconstruct the evolutionary histories of 56 pediatric primary tumors, including 24 neuroblastomas, 24 Wilms tumors, and 8 rhabdomyosarcomas. Whole-genome copy-number and whole-exome mutational profiling of multiple regions per tumor were performed, followed by clonal deconvolution to reconstruct a phylogenetic tree for each tumor. Overall, 88% of the tumors exhibited genetic variation among primary tumor regions. This variability typically emerged through collateral phylogenetic branching, leading to spatial variability in the distribution of more than 50% (96/173) of detected diagnostically informative genetic aberrations. Single-cell sequencing of 547 individual cancer cells from eight solid pediatric tumors confirmed branching evolution to be a fundamental underlying principle of genetic variation in all cases. Strikingly, cell-to-cell genetic diversity was almost twice as high in aggressive compared with clinically favorable tumors (median Simpson index of diversity 0.45 vs. 0.88; P = 0.029). Similarly, a comparison of multiregional sampling data from a total of 274 tumor regions showed that new phylogenetic branches emerge at a higher frequency per sample and carry a higher mutational load in high-risk than in low-risk tumors. Timelines based on spatial genetic variation showed that the mutations most influencing relapse risk occur at initiation of clonal expansion in neuroblastoma and rhabdomyosarcoma, whereas in Wilms tumor, they are late events. Thus, from an evolutionary standpoint, some high-risk childhood cancers are born bad, whereas others grow worse over time. SIGNIFICANCE: Different pediatric cancers with a high risk of relapse share a common generic pattern of extensively branching evolution of somatic mutations. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/7/1512/F1.large.jpg.

Therapy-induced mutations drive the genomic landscape of relapsed acute lymphoblastic leukemia

To study the mechanisms of relapse in acute lymphoblastic leukemia (ALL), we performed whole-genome sequencing of 103 diagnosis-relapse-germline trios and ultra-deep sequencing of 208 serial samples in 16 patients. Relapse-specific somatic alterations were enriched in 12 genes (NR3C1, NR3C2, TP53, NT5C2, FPGS, CREBBP, MSH2, MSH6, PMS2, WHSC1, PRPS1, and PRPS2) involved in drug response. Their prevalence was 17% in very early relapse (<9 months from diagnosis), 65% in early relapse (9-36 months), and 32% in late relapse (>36 months) groups. Convergent evolution, in which multiple subclones harbor mutations in the same drug resistance gene, was observed in 6 relapses and confirmed by single-cell sequencing in 1 case. Mathematical modeling and mutational signature analysis indicated that early relapse resistance acquisition was frequently a 2-step process in which a persistent clone survived initial therapy and later acquired bona fide resistance mutations during therapy. In contrast, very early relapses arose from preexisting resistant clone(s). Two novel relapse-specific mutational signatures, one of which was caused by thiopurine treatment based on in vitro drug exposure experiments, were identified in early and late relapses but were absent from 2540 pan-cancer diagnosis samples and 129 non-ALL relapses. The novel signatures were detected in 27% of relapsed ALLs and were responsible for 46% of acquired resistance mutations in NT5C2, PRPS1, NR3C1, and TP53. These results suggest that chemotherapy-induced drug resistance mutations facilitate a subset of pediatric ALL relapses.

Genomic Complexity of Osteosarcoma and Its Implication for Preclinical and Clinical Targeted Therapies

Osteosarcoma is a genomically complex disease characterized by few recurrent single-nucleotide mutations or in-frame fusions. In contrast, structural alterations, including copy number changes, chromothripsis, kataegis, loss of heterozygosity (LOH), and other large-scale genomic alterations, are frequent and widespread across the osteosarcoma genome. These observed structural alterations lead to activation of oncogenes and loss of tumor suppressors which together contribute to oncogenesis. To date, few targeted therapies for osteosarcoma have been identified. It is likely that effectiveness of targeted therapies will vary greatly in subsets of tumors with distinct key driver events. Model systems which can recapitulate the genetic heterogeneity of this disease are needed to test this hypothesis. One possible approach is to use patient-derived xenograft (PDX) models characterized with regards to their similarity to the human tumor samples from which they were derived. Here we review evidence pointing to the genomic complexity of osteosarcoma and how this is reflected in available model systems. We also review the current state of preclinical testing for targeted therapies using these models.

Relapsed Osteosarcoma Trial Concepts to Match the Complexity of the Disease

Osteosarcoma relapses not only herald a very poor prognosis but also opportunities to treat this genetically diverse complex cancer in new ways. This review will attempt to show that the field is a rapidly evolving one in which not only cytotoxic agents but also local control strategies and the immune system can be harnessed to improve the prognosis of relapsed patients. The molecular heterogeneity and the difficulty of effectively treating most common patterns of relapse with surgery and/or radiation (lung and/or bone metastases) have been responsible for a wide variety of approaches to learning whether agents are active against osteosarcoma. This chapter will highlight past, current, and potential future approaches to provide more effective systemic therapy for the problem of recurrent metastases of osteosarcoma. These include single-agent trials with a wide variety of agents, radiopharmaceuticals, and immune therapies. Finally, how such efforts are integrated into more effective local control strategies is also discussed.

Ubiquitin-like protein FAT10 promotes osteosarcoma growth by modifying the ubiquitination and degradation of YAP1

Osteosarcoma is a common malignancy of the bone tissue. The rapid growth exhibited by this cancer is a primary challenge in its treatment. In many types of cancers, FAT10, a ubiquitin-like protein, is involved in several biological activities, especially cell proliferation. Herein, we demonstrate that FAT10 plays a vital role in tumorigenesis and is overexpressed in tumor tissues compared to its expression in adjacent normal tissues. Functional assays revealed that knockdown of FAT10 expression significantly repressed the proliferation of osteosarcoma in vitro and in vivo. Furthermore, our results indicate that FAT10 exhibits oncogenic functions by regulating the level of YAP1, a key protein of the Hippo/YAP signaling pathway, and a significant positive correlation exists between the levels of FAT10 and YAP1. Further analysis showed that FAT10-induced growth of osteosarcoma cells is dependent on YAP1. Mechanistically, FAT10 stabilizes YAP1 expression by regulating its ubiquitination and degradation. Taken together, our results link the two drivers of cell growth in osteosarcoma and reveal a novel pathway for FAT10 regulation. We provide new evidence for the biological and clinical significance of FAT10 as a potential biomarker for osteosarcoma.

The mutational footprints of cancer therapies

Some cancer therapies damage DNA and cause mutations both in cancer and healthy cells of the patient. Therapy-induced mutations may underlie some of the long-term and late side effects of treatments, such as mental disabilities, organ toxicities and secondary neoplasms. Currently we ignore the mutation burden caused by different cancer treatments. Here we identify mutational signatures, or footprints of six widely-used anti-cancer therapies across more than 3,500 metastatic tumors originating from different organs. These include previously known and new mutational signatures generated by platinum-based drugs, and a novel signature of nucleoside metabolic inhibitors. Exploiting these mutational footprints, we estimate the contribution of different treatments to the mutation burden of tumors and their risk of contributing coding and potential driver mutations in the genome. The mutational footprints identified here allow for precisely assessing the mutational risk of different cancer therapies to understand their long-term side effects.

Osteosarcoma in the Post Genome Era: Preclinical Models and Approaches to Identify Tractable Therapeutic Targets

PURPOSE OF REVIEW

Osteosarcoma (OS) is the most common cancer of bone, yet is classified as a rare cancer. Treatment and outcomes for OS have not substantively changed in several decades. While the decoding of the OS genome greatly advanced the understanding of the mutational landscape of OS, immediately actionable therapeutic targets were not apparent. Here we describe recent preclinical models that can be leveraged to identify, test, and prioritize therapeutic candidates.

RECENT FINDINGS

The generation of multiple high fidelity murine models of OS, the spontaneous disease that arises in pet dogs, and the establishment of a diverse collection of patient-derived OS xenografts provide a robust preclinical platform for OS. These models enable evidence to be accumulated across multiple stages of preclinical evaluation. Chemical and genetic screening has identified therapeutic targets, often demonstrating cross species activity. Clinical trials in both PDX models and in canine OS have effectively tested new therapies for prioritization. Improving clinical outcomes in OS has proven elusive. The integrated target discovery and testing possible through a cross species platform provides validation of a putative target and may enable the rigorous evaluation of new therapies in models where endpoints can be rapidly assessed.

Tumor heterogeneity and acquired drug resistance in FGFR2-fusion-positive cholangiocarcinoma through rapid research autopsy

Cholangiocarcinoma is a highly aggressive and lethal malignancy, with limited treatment options available. Recently, FGFR inhibitors have been developed and utilized in FGFR-mutant cholangiocarcinoma; however, resistance often develops and the genomic determinants of resistance are not fully characterized. We completed whole-exome sequencing (WES) of 11 unique tumor samples obtained from a rapid research autopsy on a patient with FGFR-fusion-positive cholangiocarcinoma who initially responded to the pan-FGFR inhibitor, INCB054828. In vitro studies were carried out to characterize the novel FGFR alteration and secondary FGFR2 mutation identified. Multisite WES and analysis of tumor heterogeneity through subclonal inference identified four genetically distinct cancer cell populations, two of which were only observed after treatment. Additionally, WES revealed an FGFR2 N549H mutation hypothesized to confer resistance to the FGFR inhibitor INCB054828 in a single tumor sample. This hypothesis was corroborated with in vitro cell-based studies in which cells expressing FGFR2-CLIP1 fusion were sensitive to INCB054828 (IC50 value of 10.16 nM), whereas cells with the addition of the N549H mutation were resistant to INCB054828 (IC50 value of 1527.57 nM). Furthermore, the FGFR2 N549H secondary mutation displayed cross-resistance to other selective FGFR inhibitors, but remained sensitive to the nonselective inhibitor, ponatinib. Rapid research autopsy has the potential to provide unprecedented insights into the clonal evolution of cancer throughout the course of the disease. In this study, we demonstrate the emergence of a drug resistance mutation and characterize the evolution of tumor subclones within a cholangiocarcinoma disease course.