Dissociation protocols used for sarcoma tissues bias the transcriptome observed in single-cell and single-nucleus RNA sequencing.

BACKGROUND: Single-cell RNA-seq has emerged as an innovative technology used to study complex tissues and characterize cell types, states, and lineages at a single-cell level. Classification of bulk tumors by their individual cellular constituents has also created new opportunities to generate single-cell atlases for many organs, cancers, and developmental models. Despite the tremendous promise of this technology, recent evidence studying epithelial tissues and diverse carcinomas suggests the methods used for tissue processing, cell disaggregation, and preservation can significantly bias gene expression and alter the observed cell types. To determine whether sarcomas - tumors of mesenchymal origin - are subject to the same technical artifacts, we profiled patient-derived tumor explants (PDXs) propagated from three aggressive subtypes: osteosarcoma (OS), Ewing sarcoma (ES), desmoplastic small round cell tumor (DSRCT). Given the rarity of these sarcoma subtypes, we explored whether single-nuclei RNA-seq from more widely available archival frozen specimens could accurately be identified by gene expression signatures linked to tissue phenotype or pathognomonic fusion proteins. RESULTS: We systematically assessed dissociation methods across different sarcoma subtypes. We compared gene expression from single-cell and single-nucleus RNA-sequencing of 125,831 whole-cells and nuclei from ES, DSRCT, and OS PDXs. We detected warm dissociation artifacts in single-cell samples and gene length bias in single-nucleus samples. Classic sarcoma gene signatures were observed regardless of the dissociation method. In addition, we showed that dissociation method biases could be computationally corrected. CONCLUSIONS: We highlighted transcriptional biases, including warm dissociation and gene-length biases, introduced by the dissociation method for various sarcoma subtypes. This work is the first to characterize how the dissociation methods used for sc/snRNA-seq may affect the interpretation of the molecular features in sarcoma PDXs.

Enzalutamide induces cytotoxicity in desmoplastic small round cell tumor independent of the androgen receptor.

Desmoplastic Small Round Cell Tumor (DSRCT) is a rare, pediatric cancer caused by the EWSR1::WT1 fusion protein. DSRCT predominantly occurs in males, which comprise 80-90% of the patient population. While the reason for this male predominance remains unknown, one hypothesis is that the androgen receptor (AR) plays a critical role in DSRCT and elevated testosterone levels in males help drive tumor growth. Here, we demonstrate that AR is highly expressed in DSRCT relative to other fusion-driven sarcomas and that the AR antagonists enzalutamide and flutamide reduce DSRCT growth. However, despite these findings, which suggest an important role for AR in DSRCT, we show that DSRCT cell lines form xenografts in female mice at the same rate as male mice and AR depletion does not significantly alter DSRCT growth in vitro. Further, we find that AR antagonists reduce DSRCT growth in cells depleted of AR, establishing an AR-independent mechanism of action. These findings suggest that AR dependence is not the reason for male predominance in DSRCT and that AR-targeted therapies may provide therapeutic benefit primarily through an AR-independent mechanism that requires further elucidation.

Mapping the Single-cell Differentiation Landscape of Osteosarcoma.

PURPOSE: The genetic intratumoral heterogeneity observed in human osteosarcomas (OS) poses challenges for drug development and the study of cell fate, plasticity, and differentiation, processes linked to tumor grade, cell metastasis, and survival. EXPERIMENTAL DESIGN: To pinpoint errors in OS differentiation, we transcriptionally profiled 31,527 cells from a tissue-engineered model that directsMSCs toward adipogenic and osteoblastic fates. Incorporating pre-existing chondrocyte data, we applied trajectory analysis and non-negative matrix factorization (NMF) to generate the first human mesenchymal differentiation atlas. RESULTS: This 'roadmap' served as a reference to delineate the cellular composition of morphologically complex OS tumors and quantify each cell's lineage commitment. Projecting a bulk RNA-seq OS dataset onto this roadmap unveiled a correlation between a stem-like transcriptomic phenotype and poorer survival outcomes. CONCLUSIONS: Our study quantifies OS differentiation and lineage, a prerequisite to better understanding lineage-specific differentiation bottlenecks that might someday be targeted therapeutically.

Mapping the Single-cell Differentiation Landscape of Osteosarcoma.

The genetic and intratumoral heterogeneity observed in human osteosarcomas (OS) poses challenges for drug development and the study of cell fate, plasticity, and differentiation, processes linked to tumor grade, cell metastasis, and survival. To pinpoint errors in OS differentiation, we transcriptionally profiled 31,527 cells from a tissue-engineered model that directs MSCs toward adipogenic and osteoblastic fates. Incorporating pre-existing chondrocyte data, we applied trajectory analysis and non-negative matrix factorization (NMF) to generate the first human mesenchymal differentiation atlas. This 'roadmap' served as a reference to delineate the cellular composition of morphologically complex OS tumors and quantify each cell's lineage commitment. Projecting these signatures onto a bulk RNA-seq OS dataset unveiled a correlation between a stem-like transcriptomic phenotype and poorer survival outcomes. Our study takes the critical first step in accurately quantifying OS differentiation and lineage, a prerequisite to better understanding global differentiation bottlenecks that might someday be targeted therapeutically.

Enzalutamide Induces Cytotoxicity in Desmoplastic Small Round Cell Tumor Independent of the Androgen Receptor.

Desmoplastic Small Round Cell Tumor (DSRCT) is a rare, pediatric cancer caused by the EWSR1::WT1 fusion protein. DSRCT predominantly occurs in males, which comprise 80-90% of the patient population. While the reason for this male predominance remains unknown, one hypothesis is that the androgen receptor (AR) plays a critical role in DSRCT and elevated testosterone levels in males help drive tumor growth. Here, we demonstrate that AR is highly expressed in DSRCT relative to other fusion-driven sarcomas and that the AR antagonists enzalutamide and flutamide reduce DSRCT growth. However, despite these findings, which suggest an important role for AR in DSRCT, we show that DSRCT cell lines form xenografts in female mice at the same rate as male mice and AR depletion does not significantly alter DSRCT growth in vitro. Further, we find that AR antagonists reduce DSRCT growth in cells depleted of AR, establishing an AR-independent mechanism of action. These findings suggest that AR dependence is not the reason for male predominance in DSRCT and that AR-targeted therapies may provide therapeutic benefit primarily through an AR-independent mechanism that requires further elucidation.

Targeting the IGF/PI3K/mTOR pathway and AXL/YAP1/TAZ pathways in primary bone cancer.

Primary bone cancers (PBC) belong to the family of mesenchymal tumors classified based on their cellular origin, extracellular matrix, genetic regulation, and epigenetic modification. The three major PBC types, Ewing sarcoma, osteosarcoma, and chondrosarcoma, are frequently aggressive tumors, highly metastatic, and typically occur in children and young adults. Despite their distinct origins and pathogenesis, these sarcoma subtypes rely upon common signaling pathways to promote tumor progression, metastasis, and survival. The IGF/PI3K/mTOR and AXL/YAP/TAZ pathways, in particular, have gained significant attention recently given their ties to oncogenesis, cell fate and differentiation, metastasis, and drug resistance. Naturally, these pathways - and their protein constituents - have caught the eye of the pharmaceutical industry, and a wide array of small molecule inhibitors and antibody drug-conjugates have emerged. Here, we review how the IGF/PI3K/mTOR and AXL/YAP/TAZ pathways promote PBC and highlight the drug candidates under clinical trial investigation.

Correlation of nuclear pIGF-1R/IGF-1R and YAP/TAZ in a tissue microarray with outcomes in osteosarcoma patients.

Osteosarcoma (OS) is a genetically diverse bone cancer that lacks a consistent targetable mutation. Recent studies suggest the IGF/PI3K/mTOR pathway and YAP/TAZ paralogs regulate cell fate and proliferation in response to biomechanical cues within the tumor microenvironment. How this occurs and their implication upon osteosarcoma survival, remains poorly understood. Here, we show that IGF-1R can translocate into the nucleus, where it may act as part of a transcription factor complex. To explore the relationship between YAP/TAZ and total and nuclear phosphorylated IGF-1R (pIGF-1R), we evaluated sequential tumor sections from a 37-patient tissue microarray by confocal microscopy. Next, we examined the relationship between stained markers, clinical disease characteristics, and patient outcomes. The nuclear to cytoplasmic ratios (N:C ratio) of YAP and TAZ strongly correlated with nuclear pIGF-1R (r = 0.522, p = 0.001 for each pair). Kaplan-Meier analyses indicated that nuclear pIGF-1R predicted poor overall survival, a finding confirmed in the Cox proportional hazards model. Though additional investigation in a larger prospective study will be required to validate the prognostic accuracy of these markers, our results may have broad implications for the new class of YAP, TAZ, AXL, or TEAD inhibitors that have reached early phase clinical trials this year.

Multi-site desmoplastic small round cell tumors are genetically related and immune-cold.

Desmoplastic small round cell tumor (DSRCT) is a highly aggressive soft tissue sarcoma that is characterized by the EWSR1-WT1 fusion protein. Patients present with hundreds of tumor implants in their abdominal cavity at various sites. To determine the genetic relatedness among these sites, exome and RNA sequencing were performed on 22 DSRCT specimens from 14 patients, four of whom had specimens from various tissue sites. Multi-site tumors from individual DSRCT patients had a shared origin and were highly related. Other than the EWSR1-WT1 fusion, very few secondary cancer gene mutations were shared among the sites. Among these, ARID1A, was recurrently mutated, which corroborates findings by others in DSRCT patients. Knocking out ARID1A in JN-DSRCT cells using CRISPR/CAS9 resulted in significantly lower cell proliferation and increased drug sensitivity. The transcriptome data were integrated using network analysis and drug target database information to identify potential therapeutic opportunities in EWSR1-WT1-associated pathways, such as PI3K and mTOR pathways. Treatment of JN-DSRCT cells with the PI3K inhibitor alpelisib and mTOR inhibitor temsirolimus reduced cell proliferation. In addition, the low mutation burden was associated with an immune-cold state in DSRCT. Together, these data reveal multiple genomic and immune features of DSRCT and suggest therapeutic opportunities in patients.

The androgen receptor is a therapeutic target in desmoplastic small round cell sarcoma.

Desmoplastic small round cell tumor (DSRCT) is an aggressive, usually incurable sarcoma subtype that predominantly occurs in post-pubertal young males. Recent evidence suggests that the androgen receptor (AR) can promote tumor progression in DSRCTs. However, the mechanism of AR-induced oncogenic stimulation remains undetermined. Herein, we demonstrate that enzalutamide and AR-directed antisense oligonucleotides (AR-ASO) block 5α-dihydrotestosterone (DHT)-induced DSRCT cell proliferation and reduce xenograft tumor burden. Gene expression analysis and chromatin immunoprecipitation sequencing (ChIP-seq) were performed to elucidate how AR signaling regulates cellular epigenetic programs. Remarkably, ChIP-seq revealed novel DSRCT-specific AR DNA binding sites adjacent to key oncogenic regulators, including WT1 (the C-terminal partner of the pathognomonic fusion protein) and FOXF1. Additionally, AR occupied enhancer sites that regulate the Wnt pathway, neural differentiation, and embryonic organ development, implicating AR in dysfunctional cell lineage commitment. Our findings have direct clinical implications given the widespread availability of FDA-approved androgen-targeted agents used for prostate cancer.

Transcriptional activators YAP/TAZ and AXL orchestrate dedifferentiation, cell fate, and metastasis in human osteosarcoma.

Osteosarcoma (OS) is a molecularly heterogeneous, aggressive, poorly differentiated pediatric bone cancer that frequently spreads to the lung. Relatively little is known about phenotypic and epigenetic changes that promote lung metastases. To identify key drivers of metastasis, we studied human CCH-OS-D OS cells within a previously described rat acellular lung (ACL) model that preserves the native lung architecture, extracellular matrix, and capillary network. This system identified a subset of cells-termed derived circulating tumor cells (dCTCs)-that can migrate, intravasate, and spread within a bioreactor-perfused capillary network. Remarkably, dCTCs highly expressed epithelial-to-mesenchymal transition (EMT)-associated transcription factors (EMT-TFs), such as ZEB1, TWIST, and SOX9, which suggests that they undergo cellular reprogramming toward a less differentiated state by coopting the same epigenetic machinery used by carcinomas. Since YAP/TAZ and AXL tightly regulate the fate and plasticity of normal mesenchymal cells in response to microenvironmental cues, we explored whether these proteins contributed to OS metastatic potential using an isogenic pair of human OS cell lines that differ in AXL expression. We show that AXL inhibition significantly reduced the number of MG63.2 pulmonary metastases in murine models. Collectively, we present a laboratory-based method to detect and characterize a pure population of dCTCs, which provides a unique opportunity to study how OS cell fate and differentiation contributes to metastatic potential. Though the important step of clinical validation remains, our identification of AXL, ZEB1, and TWIST upregulation raises the tantalizing prospect that EMT-TF-directed therapies might expand the arsenal of therapies used to combat advanced-stage OS.

A Human Organotypic Microfluidic Tumor Model Permits Investigation of the Interplay between Patient-Derived Fibroblasts and Breast Cancer Cells.

Tumor-stroma interactions significantly influence cancer cell metastasis and disease progression. These interactions are partly comprised of the cross-talk between tumor and stromal fibroblasts, but the key molecular mechanisms within the cross-talk that govern cancer invasion are still unclear. Here, we adapted our previously developed microfluidic device as a 3D in vitro organotypic model to mechanistically study tumor-stroma interactions by mimicking the spatial organization of the tumor microenvironment on a chip. We cocultured breast cancer and patient-derived fibroblast cells in 3D tumor and stroma regions, respectively, and combined functional assessments, including cancer cell migration, with transcriptome profiling to unveil the molecular influence of tumor-stroma cross-talk on invasion. This led to the observation that cancer-associated fibroblasts (CAF) enhanced invasion in 3D by inducing expression of a novel gene of interest, glycoprotein nonmetastatic B (GPNMB), in breast cancer cells, resulting in increased migration speed. Importantly, knockdown of GPNMB blunted the influence of CAF on enhanced cancer invasion. Overall, these results demonstrate the ability of our model to recapitulate patient-specific tumor microenvironments to investigate the cellular and molecular consequences of tumor-stroma interactions. SIGNIFICANCE: An organotypic model of tumor-stroma interactions on a microfluidic chip reveals that CAFs promote invasion by enhancing expression of GPNMB in breast cancer cells.