Clinically relevant treatment of PDX models reveals patterns of neuroblastoma chemoresistance

Advances in the approaches used to repurpose drugs for neuroblastoma

INTRODUCTION

Neuroblastoma (NB) remains a challenging pediatric malignancy with limited treatment options, particularly for high-risk cases. Drug repurposing offers a convenient and cost-effective strategy for treating rare diseases like NB. Using existing drugs with known safety profiles accelerates the availability of new treatments, reduces development costs, and mitigates risks, offering hope for improved patient outcomes in challenging conditions.

AREAS COVERED

This review provides an overview of the advances in approaches used to repurpose drugs for NB therapy. The authors discuss strategies employed in drug repurposing, including computational and experimental methods, and rational drug design, highlighting key examples of repurposed drugs with promising clinical results. Additionally, the authors examine the challenges and opportunities associated with drug repurposing in NB and discuss future directions and potential areas for further research.

EXPERT OPINION

The fact that only one new drug has been approved in the last 30 years for the treatment of neuroblastoma plus a significant proportion of high-risk NB patients that remain uncurable, evidences the need for new fast and cost-effective alternatives. Drug repurposing may accelerate the treatment development process while reducing expenses and risks. This approach can swiftly bring effective NB therapies to market, enhancing survival rates and patient quality of life.

NBAtlas: A harmonized single-cell transcriptomic reference atlas of human neuroblastoma tumors

Neuroblastoma, a rare embryonic tumor arising from neural crest development, is responsible for 15% of pediatric cancer-related deaths. Recently, several single-cell transcriptome studies were performed on neuroblastoma patient samples to investigate the cell of origin and tumor heterogeneity. However, these individual studies involved a small number of tumors and cells, limiting the conclusions that could be drawn. To overcome this limitation, we integrated seven single-cell or single-nucleus datasets into a harmonized cell atlas covering 362,991 cells across 61 patients. We use this atlas to decipher the transcriptional landscape of neuroblastoma at single-cell resolution, revealing associations between transcriptomic profiles and clinical outcomes within the tumor compartment. In addition, we characterize the complex immune-cell landscape and uncover considerable heterogeneity among tumor-associated macrophages. Finally, we showcase the utility of our atlas as a resource by expanding it with additional data and using it as a reference for data-driven cell-type annotation.

Early evolutionary branching across spatial domains predisposes to clonal replacement under chemotherapy in neuroblastoma

Neuroblastoma (NB) is one of the most lethal childhood cancers due to its propensity to become treatment resistant. By spatial mapping of subclone geographies before and after chemotherapy across 89 tumor regions from 12 NBs, we find that densely packed territories of closely related subclones present at diagnosis are replaced under effective treatment by islands of distantly related survivor subclones, originating from a different most recent ancestor compared to lineages dominating before treatment. Conversely, in tumors that progressed under treatment, ancestors of subclones dominating later in disease are present already at diagnosis. Chemotherapy treated xenografts and cell culture models replicate these two contrasting scenarios and show branching evolution to be a constant feature of proliferating NB cells. Phylogenies based on whole genome sequencing of 505 individual NB cells indicate that a rich repertoire of parallel subclones emerges already with the first oncogenic mutations and lays the foundation for clonal replacement under treatment.

MELATONIN ENHANCES TEMOZOLOMIDE-INDUCED APOPTOSIS IN GLIOBLASTOMA AND NEUROBLASTOMA CELLS

BACKGROUND

The combination of temozolomide (TMZ) and paclitaxel (PTX) is the most commonly used chemotherapy regimen for glioblastoma, but there is no specific treatment for neuroblastoma due to the acquired multidrug resistance. Approximately half of treated glioblastoma patients develop resistance to TMZ and experience serious side effects. Melatonin (MEL), a multifunctional hormone long known for its antitumor effects, has a great advantage in combination cancer therapy thanks to its ability to affect tumors differently than normal cells.

AIM

This study aims to evaluate the in vitro inhibitory effects of MEL in combination with TMZ on cancer cell viability and to elucidate the underlying mechanisms in the glioblastoma and neuroblastoma cell lines.

MATERIALS AND METHODS

C6 (Rattus norvegicus) and N1E-115 (Mus musculus) cancer cell lines and C8-D1A (mice) healthy cell lines were used. Cell proliferation was evaluated using the MTT test. IC50 values were determined by probit analysis. Two concentrations of TMZ (IC50 and 1/2 IC50) were used to induce cytotoxicity in the C6 and N1E-115 cell lines, both alone and in combination with PXT and MEL (all at IC50). The viable, dead, and apoptotic cells were determined by image-based cytometry using Annexin V/PI staining. The gene expression related to signaling pathways was assessed by the quantitative reverse transcription polymerase chain reaction (qRT-PCR), and key proteins were identified by the Western blot analysis.

RESULTS

MTT assay showed that the combination of TMZ and MEL significantly reduces the viability of both glioblastoma and neuroblastoma cells compared to the vehicle-treated controls. Notably, MEL combined with 1/2 IC50 TMZ showed a significant death rate of cancer cells compared to controls and PTX. According to qRT-PCR data, the TMZ + MEL combination resulted in the upregulation of the genes of antioxidative enzymes (Sod1 and Sod2) and DNA repair genes (Mlh1, Exo1, and Rad18) in both cell lines. Moreover, the levels of Nfkb1 and Pik3cg were significantly reduced following the TMZ + MEL treatment. The combination of MEL with TMZ also enhanced the cell cycle arrest and increased the expression of p53 and pro-apoptotic proteins (Bax and caspase-3), while significantly decreasing the expression of anti-apoptotic protein Bcl-2.

CONCLUSIONS

Our findings indicate that the combination of MEL with a low dose of TMZ may serve as an upstream inducer of apoptosis. This suggests the potential development of a novel selective therapeutic strategy as an alternative to TMZ for the treatment of both glioblastoma and neuroblastoma.

Targeting SWI/SNF ATPases reduces neuroblastoma cell plasticity

Tumor cell heterogeneity defines therapy responsiveness in neuroblastoma (NB), a cancer derived from neural crest cells. NB consists of two primary subtypes: adrenergic and mesenchymal. Adrenergic traits predominate in NB tumors, while mesenchymal features becomes enriched post-chemotherapy or after relapse. The interconversion between these subtypes contributes to NB lineage plasticity, but the underlying mechanisms driving this phenotypic switching remain unclear. Here, we demonstrate that SWI/SNF chromatin remodeling complex ATPases are essential in establishing an mesenchymal gene-permissive chromatin state in adrenergic-type NB, facilitating lineage plasticity. Targeting SWI/SNF ATPases with SMARCA2/4 dual degraders effectively inhibits NB cell proliferation, invasion, and notably, cellular plasticity, thereby preventing chemotherapy resistance. Mechanistically, depletion of SWI/SNF ATPases compacts cis-regulatory elements, diminishes enhancer activity, and displaces core transcription factors (MYCN, HAND2, PHOX2B, and GATA3) from DNA, thereby suppressing transcriptional programs associated with plasticity. These findings underscore the pivotal role of SWI/SNF ATPases in driving intrinsic plasticity and therapy resistance in neuroblastoma, highlighting an epigenetic target for combinational treatments in this cancer.

Characterizing Relationships between T-cell Inflammation and Outcomes in Patients with High-Risk Neuroblastoma According to Mesenchymal and Adrenergic Signatures

Recent insights have identified adrenergic (ADRN) and mesenchymal (MES) cell lineages as distinct biologic cell types and T-cell inflammation as a prognostic marker in neuroblastoma. We hypothesized that elucidating unique and overlapping aspects of these biologic features could serve as novel biomarkers for informing ongoing efforts to improve therapeutic approaches for children with high-risk neuroblastoma. We identified lineage-specific, single-stranded super-enhancers to define ADRN and MES specific genes. Publicly available RNA-seq of diagnostic tumor biopsies was used in Discovery and Validation cohorts. Each tumor was assigned a relative MES score and T-cell inflammation (TCI) score. Survival was assessed using the Kaplan-Meier method, and differences were assessed by the log-rank test. Inflammation scores were correlated with MES scores and anticorrelated with MYCN-amplification in both cohorts. Among patients with high-risk, ADRN tumors, those with TCI tumors had superior overall survival to those with non-inflamed tumors. A similar, but nonsignificant, trend was observed in the Validation cohort. Conversely, there was no difference according to TCI status in the MES cohort in either the Discover or Validation cohorts. High-inflammation scores were correlated with improved survival in some patients with high-risk, ADRN but not MES neuroblastoma. Our findings bolster support for further developing T-cell-based and immunotherapy-based approaches for children with high-risk neuroblastoma of varying MES and ADRN expression.

SIGNIFICANCE

Adrenergic (ADRN) and mesenchymal (MES) lineages are distinct biologic cell types in neuroblastoma. We defined ADRN and MES specific genes and found that high-risk, ADRN tumors harboring elevated T-cell inflammation signatures had superior overall survival. Our findings bolster support for further developing immunotherapy-based approaches for children with high-risk neuroblastoma.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs

BACKGROUND

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models.

RESULTS

Here, we generate single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We develop an unsupervised machine learning approach ("automatic consensus nonnegative matrix factorization" (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirm a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly, however, this weak-mesenchymal-like program is maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 h, suggesting an uncharacterized therapy-escape mechanism.

CONCLUSIONS

Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed an unsupervised machine learning approach ('automatic consensus nonnegative matrix factorization' (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirmed a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly however, this weak-mesenchymal-like program was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

Preclinical Anticipation of On- and Off-Target Resistance Mechanisms to Anti-Cancer Drugs: A Systematic Review

The advent of targeted therapies has led to tremendous improvements in treatment options and their outcomes in the field of oncology. Yet, many cancers outsmart precision drugs by developing on-target or off-target resistance mechanisms. Gaining the ability to resist treatment is the rule rather than the exception in tumors, and it remains a major healthcare challenge to achieve long-lasting remission in most cancer patients. Here, we discuss emerging strategies that take advantage of innovative high-throughput screening technologies to anticipate on- and off-target resistance mechanisms before they occur in treated cancer patients. We divide the methods into non-systematic approaches, such as random mutagenesis or long-term drug treatment, and systematic approaches, relying on the clustered regularly interspaced short palindromic repeats (CRISPR) system, saturated mutagenesis, or computational methods. All these new developments, especially genome-wide CRISPR-based screening platforms, have significantly accelerated the processes for identification of the mechanisms responsible for cancer drug resistance and opened up new avenues for future treatments.

Two bullets in the gun: combining immunotherapy with chemotherapy to defeat neuroblastoma by targeting adrenergic-mesenchymal plasticity

Neuroblastoma (NB) is a childhood tumor that originates in the peripheral sympathetic nervous system and is responsible for 15% of cancer-related deaths in the pediatric population. Despite intensive multimodal treatment, many patients with high-risk NB relapse and develop a therapy-resistant tumor. One of the phenomena related to therapeutic resistance is intratumor heterogeneity resulting from the adaptation of tumor cells in response to different selective environmental pressures. The transcriptional and epigenetic profiling of NB tissue has recently revealed the existence of two distinct cellular identities in the NB, termed adrenergic (ADRN) and mesenchymal (MES), which can spontaneously interconvert through epigenetic regulation. This phenomenon, known as tumor plasticity, has a major impact on cancer pathogenesis. The aim of this review is to describe the peculiarities of these two cell states, and how their plasticity affects the response to current therapeutic treatments, with special focus on the immunogenic potential of MES cells. Furthermore, we will discuss the opportunity to combine immunotherapy with chemotherapy to counteract NB phenotypic interconversion.

T-cell inflammation is prognostic of survival in patients with high-risk neuroblastoma enriched for an adrenergic signature

Purpose

T-cell inflammation (TCI) has been shown to be a prognostic marker in neuroblastoma, a tumor comprised of cells that can exist in two epigenetic states, adrenergic (ADRN) and mesenchymal (MES). We hypothesized that elucidating unique and overlapping aspects of these biologic features could serve as novel biomarkers.

Patients and Methods

We detected lineage-specific, single-stranded super-enhancers defining ADRN and MES specific genes. Publicly available neuroblastoma RNA-seq data from GSE49711 (Cohort 1) and TARGET (Cohort 2) were assigned MES, ADRN, and TCI scores. Tumors were characterized as MES (top 33%) or ADRN (bottom 33%), and TCI (top 67% TCI score) or non-inflamed (bottom 33% TCI score). Overall survival (OS) was assessed using the Kaplan-Meier method, and differences were assessed by the log-rank test.

Results

We identified 159 MES genes and 373 ADRN genes. TCI scores were correlated with MES scores (R=0.56, p<0.001 and R=0.38, p<0.001) and anticorrelated with MYCN -amplification (R=-0.29, p<0.001 and -0.18, p=0.03) in both cohorts. Among Cohort 1 patients with high-risk, ADRN tumors (n=59), those with TCI tumors (n=22) had superior OS to those with non-inflammed tumors (n=37) (p=0.01), though this comparison did not reach significance in Cohort 2. TCI status was not associated with survival in patients with high-risk MES tumors in either cohort.

Conclusions

High inflammation scores were correlated with improved survival in some high-risk patients with, ADRN but not MES neuroblastoma. These findings have implications for approaches to treating high-risk neuroblastoma.

Using 3D-bioprinted models to study pediatric neural crest-derived tumors

The use of three-dimensional (3D) bioprinting has remained at the forefront of tissue engineering and has recently been employed for generating bioprinted solid tumors to be used as cancer models to test therapeutics. In pediatrics, neural crest-derived tumors are the most common type of extracranial solid tumors. There are only a few tumor-specific therapies that directly target these tumors, and the lack of new therapies remains detrimental to improving the outcomes for these patients. The absence of more efficacious therapies for pediatric solid tumors, in general, may be due to the inability of the currently employed preclinical models to recapitulate the solid tumor phenotype. In this study, we utilized 3D bioprinting to generate neural crest-derived solid tumors. The bioprinted tumors consisted of cells from established cell lines and patient-derived xenograft tumors mixed with a 6% gelatin/1% sodium alginate bioink. The viability and morphology of the bioprints were analyzed via bioluminescence and immunohisto chemistry, respectively. We compared the bioprints to traditional twodimensional (2D) cell culture under conditions such as hypoxia and therapeutics. We successfully produced viable neural crest-derived tumors that retained the histology and immunostaining characteristics of the original parent tumors. The bioprinted tumors propagated in culture and grew in orthotopic murine models. Furthermore, compared to cells grown in traditional 2D culture, the bioprinted tumors were resistant to hypoxia and chemotherapeutics, suggesting that the bioprints exhibited a phenotype that is consistent with that seen clinically in solid tumors, thus potentially making this model superior to traditional 2D culture for preclinical investigations. Future applications of this technology entail the potential to rapidly print pediatric solid tumors for use in high-throughput drug studies, expediting the identification of novel, individualized therapies.

Analysis of Serial Neuroblastoma PDX Passages in Mice Allows the Identification of New Mediators of Neuroblastoma Aggressiveness

Neuroblastoma is a neural crest cell-derived pediatric tumor characterized by high inter- and intra-tumor heterogeneity, and by a poor outcome in advanced stages. Patient-derived xenografts (PDXs) have been shown to be useful models for preserving and expanding original patient biopsies in vivo, and for studying neuroblastoma biology in a more physiological setting. The maintenance of genetic, histologic, and phenotypic characteristics of the original biopsy along serial PDX passages in mice is a major concern regarding this model. Here we analyze consecutive PDX passages in mice, at both transcriptomic and histological levels, in order to identify potential changes or highlight similarities to the primary sample. We studied temporal changes using mRNA and miRNA expression and correlate those with neuroblastoma aggressiveness using patient-derived databases. We observed a shortening of tumor onset and an increase in proliferative potential in the PDXs along serial passages. This behavior correlates with changes in the expression of genes related to cell proliferation and neuronal differentiation, including signaling pathways described as relevant for neuroblastoma malignancy. We also identified new genes and miRNAs that can be used to stratify patients according to survival, and which could be potential new players in neuroblastoma aggressiveness. Our results highlight the usefulness of the PDX neuroblastoma model and reflect phenotypic changes that might be occurring in the mouse environment. These findings could be useful for understanding the progression of tumor aggressiveness in this pathology.

Macrophage infiltration promotes regrowth in MYCN-amplified neuroblastoma after chemotherapy

Despite aggressive treatment, the 5-year event-free survival rate for children with high-risk neuroblastoma is <50%. While most high-risk neuroblastoma patients initially respond to treatment, often with complete clinical remission, many eventually relapse with therapy-resistant tumors. Novel therapeutic alternatives that prevent the recurrence of therapy-resistant tumors are urgently needed. To understand the adaptation of neuroblastoma under therapy, we analyzed the transcriptomic landscape in 46 clinical tumor samples collected before (PRE) or after (POST) treatment from 22 neuroblastoma patients. RNA sequencing revealed that many of the top-upregulated biological processes in POST MYCN amplified (MNA⁺) tumors compared to PRE MNA⁺ tumors were immune-related, and there was a significant increase in numerous genes associated with macrophages. The infiltration of macrophages was corroborated by immunohistochemistry and spatial digital protein profiling. Moreover, POST MNA⁺ tumor cells were more immunogenic compared to PRE MNA⁺ tumor cells. To find support for the macrophage-induced outgrowth of certain subpopulations of immunogenic tumor cells following treatment, we examined the genetic landscape in multiple clinical PRE and POST tumor samples from nine neuroblastoma patients revealing a significant correlation between an increased amount of copy number aberrations (CNA) and macrophage infiltration in POST MNA⁺ tumor samples. Using an in vivo neuroblastoma patient-derived xenograft (PDX) chemotherapy model, we further show that inhibition of macrophage recruitment with anti-CSF1R treatment prevents the regrowth of MNA⁺ tumors following chemotherapy. Taken together, our work supports a therapeutic strategy for fighting the relapse of MNA⁺ neuroblastoma by targeting the immune microenvironment.