Identification of an early survival prognostic gene signature for localized osteosarcoma patients.

Osteosarcoma is the most prevalent bone tumor in pediatric patients. Neoadjuvant chemotherapy has improved osteosarcoma patient survival, however the 5-year survival rate for localized osteosarcoma is 75% with a 30-50% recurrence rate. We, therefore, sought to identify a prognostic gene signature which could predict poor prognosis in localized osteosarcoma patients. Using the TARGET osteosarcoma transcriptomic dataset, we identified a 13-hub gene signature associated with overall survival and time to death of localized osteosarcoma patients, with the high-risk group showing a 22% and the low-risk group showing 100% overall survival. Furthermore, network analysis identified five modules of co-expressed genes that significantly correlated with survival, and identified 65 pathways enriched across 3 modules, including Hedgehog signaling, which includes 2 of the 13 genes, IHH and GLI1. Subsequently, we demonstrated that GLI antagonists inhibited growth of a recurrent localized PDX-derived cell line with elevated IHH and GLI1 expression, but not a non-relapsed cell line with low pathway activation. Finally, we show that our signature outperforms previously reported signatures in predicting poor prognosis and death within 3 years in patients with localized osteosarcoma.

SOD1 is a synthetic lethal target in PPM1D-mutant leukemia cells.

The DNA damage response is critical for maintaining genome integrity and is commonly disrupted in the development of cancer. PPM1D (protein phosphatase, Mg2+/Mn2+ dependent 1D) is a master negative regulator of the response; gain-of-function mutations and amplifications of PPM1D are found across several human cancers making it a relevant pharmacologic target. Here, we used CRISPR/Cas9 screening to identify synthetic-lethal dependencies of PPM1D, uncovering superoxide dismutase-1 (SOD1) as a potential target for PPM1D-mutant cells. We revealed a dysregulated redox landscape characterized by elevated levels of reactive oxygen species and a compromised response to oxidative stress in PPM1D-mutant cells. Altogether, our results demonstrate the protective role of SOD1 against oxidative stress in PPM1D-mutant leukemia cells and highlight a new potential therapeutic strategy against PPM1D-mutant cancers.

SOD1 is a synthetic-lethal target in PPM1D-mutant leukemia cells.

The DNA damage response is critical for maintaining genome integrity and is commonly disrupted in the development of cancer. PPM1D (protein phosphatase Mg2+/Mn2+-dependent 1D) is a master negative regulator of the response; gain-of-function mutations and amplifications of PPM1D are found across several human cancers making it a relevant pharmacological target. Here, we used CRISPR/Cas9 screening to identify synthetic-lethal dependencies of PPM1D, uncovering superoxide dismutase-1 (SOD1) as a potential target for PPM1D-mutant cells. We revealed a dysregulated redox landscape characterized by elevated levels of reactive oxygen species and a compromised response to oxidative stress in PPM1D-mutant cells. Altogether, our results demonstrate a role for SOD1 in the survival of PPM1D-mutant leukemia cells and highlight a new potential therapeutic strategy against PPM1D-mutant cancers.

Functional genomic analysis of non-canonical DNA regulatory elements of the aryl hydrocarbon receptor.

The aryl hydrocarbon receptor (AHR) is a ligand-dependent transcription factor that binds DNA and regulates genes in response to halogenated and polycyclic aromatic hydrocarbons. AHR also regulates the development and function of the liver and the immune system. In the canonical pathway, AHR binds a consensus DNA sequence, termed the xenobiotic response element (XRE), recruits protein coregulators, and regulates target gene expression. Emerging evidence suggests that AHR may regulate gene expression via an additional pathway, by binding to a non-consensus DNA sequence termed the non-consensus XRE (NC-XRE). The prevalence of NC-XRE motifs in the genome is not known. Studies using chromatin immunoprecipitation and reporter genes provide indirect evidence of AHR-NC-XRE interactions, but direct evidence for an AHR-NCXRE interaction that regulates transcription in a natural genomic context is lacking. Here, we analyzed AHR binding to NC-XRE DNA on a genome-wide scale in mouse liver. We integrated ChIP-seq and RNA-seq data and identified putative AHR target genes with NC-XRE motifs in regulatory regions. We also performed functional genomics at a single locus, the mouse Serpine1 gene. Deleting NC-XRE motifs from the Serpine1 promoter reduced the upregulation of Serpine1 by TCDD, an AHR ligand. We conclude that AHR upregulates Serpine1 via NC-XRE DNA. NC-XRE motifs are prevalent throughout regions of the genome where AHR binds. Taken together, our results suggest that AHR regulates genes via NC-XRE motifs. Our results will also improve our ability to identify AHR target genes and their physiologic relevance.

Modeling Ewing Sarcoma Lung Metastasis.

Ewing Sarcoma (EwS) is the second most common malignant bone tumor in adolescents and young adults. The single-most powerful predictor of outcome in EwS is presence of metastatic burden at the time of diagnosis. Patients with metastatic Ewing Sarcoma have an abysmal 5-year survival rate of 10%-25%, which has not changed over the past 30-40 years. Thus, unraveling underlying mechanisms of EwS metastasis are imperative for developing effective therapeutic measures. Investigations towards this goal are limited by the lack of reliable genetically engineered mouse models and specialized metastatic models. Using two established cell lines, A673 and TC71, we generated lung specific metastatic cell lines by serial orthotopic intra-tibial injection followed by isolation of cells from lung metastases. The lung metastatic lines generated exhibit distinct differential molecular signatures from the parental cells when analyzed using a multi-omics approach. These signatures overlapped with EwS patient primary bone and metastatic lung specimens supporting the clinical relevance of these preclinical models of EwS. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Intra-Tibial injection in NSG mice Basic Protocol 2: Development and characterization of lung metastatic cell line.

MYC regulates CSF1 expression via microRNA 17/20a to modulate tumor-associated macrophages in osteosarcoma.

Osteosarcoma (OS) is the most common primary bone tumor of childhood. Approximately 20%-30% of OSs carry amplification of chromosome 8q24, which harbors the oncogene c-MYC and correlates with a poor prognosis. To understand the mechanisms that underlie the ability of MYC to alter both the tumor and its surrounding tumor microenvironment (TME), we generated and molecularly characterized an osteoblast-specific Cre-Lox-Stop-Lox-c-MycT58A p53fl/+ knockin genetically engineered mouse model (GEMM). Phenotypically, the Myc-knockin GEMM had rapid tumor development with a high incidence of metastasis. MYC-dependent gene signatures in our murine model demonstrated significant homology to the human hyperactivated MYC OS. We established that hyperactivation of MYC led to an immune-depleted TME in OS demonstrated by the reduced number of leukocytes, particularly macrophages. MYC hyperactivation led to the downregulation of macrophage colony-stimulating factor 1, through increased microRNA 17/20a expression, causing a reduction of macrophage population in the TME of OS. Furthermore, we developed cell lines from the GEMM tumors, including a degradation tag-MYC model system, which validated our MYC-dependent findings both in vitro and in vivo. Our studies utilized innovative and clinically relevant models to identify a potentially novel molecular mechanism through which MYC regulates the profile and function of the OS immune landscape.

Integrative metabolomics and transcriptomics analysis reveals novel therapeutic vulnerabilities in lung cancer.

BACKGROUND: Non-small cell lung cancer (NSCLC) comprises the majority (~85%) of all lung tumors, with lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC) being the most frequently diagnosed histological subtypes. Multi-modal omics profiling has been carried out in NSCLC, but no studies have yet reported a unique metabolite-related gene signature and altered metabolic pathways associated with LUAD and LUSC. METHODS: We integrated transcriptomics and metabolomics to analyze 30 human lung tumors and adjacent noncancerous tissues. Differential co-expression was used to identify modules of metabolites that were altered between normal and tumor. RESULTS: We identified unique metabolite-related gene signatures specific for LUAD and LUSC and key pathways aberrantly regulated at both transcriptional and metabolic levels. Differential co-expression analysis revealed that loss of coherence between metabolites in tumors is a major characteristic in both LUAD and LUSC. We identified one metabolic onco-module gained in LUAD, characterized by nine metabolites and 57 metabolic genes. Multi-omics integrative analysis revealed a 28 metabolic gene signature associated with poor survival in LUAD, with six metabolite-related genes as individual prognostic markers. CONCLUSIONS: We demonstrated the clinical utility of this integrated metabolic gene signature in LUAD by using it to guide repurposing of AZD-6482, a PI3Kß inhibitor which significantly inhibited three genes from the 28-gene signature. Overall, we have integrated metabolomics and transcriptomics analyses to show that LUAD and LUSC have distinct profiles, inferred gene signatures with prognostic value for patient survival, and identified therapeutic targets and repurposed drugs for potential use in NSCLC treatment.

LncRNA PVT-1 promotes osteosarcoma cancer stem-like properties through direct interaction with TRIM28 and TSC2 ubiquitination.

Osteosarcoma, the most common pediatric bone tumor, is an aggressive heterogeneous malignancy defined by complex chromosomal aberrations. Overall survival rates remain at ~70%, but patients with chemoresistant or metastatic disease have extremely poor outcomes of <30%. A subgroup of tumors harbor amplification of chromosome 8q24.2 and increased expression of the oncogenic long noncoding RNA (lncRNA) Plasmacytoma Variant Translocation-1 (PVT-1), which is associated with an extremely poor clinical prognosis. This study demonstrates that PVT-1 is critical for osteosarcoma tumor-initiation potential. Chromatin Hybridization by RNA Purification analysis identified Tripartite-Motif Containing Family 28 (TRIM28) as a novel PVT-1 binding partner. Mechanistically, co-immunoprecipitation studies showed the PVT-1/TRIM28 complex binds and increases SUMOylation of phosphatidylinositol 3-kinase catalytic subunit type 3 (Vps34), which leads to enhanced ubiquitination and degradation of tumor suppressor complex 2 (TSC2), thus contributing to increased self-renewal and stem cell phenotypes. Furthermore, we identified that osteosarcoma cells with increased PVT-1 have enhanced sensitivity to the SUMOylation inhibitor, TAK-981. Altogether, this study elucidated a role for PVT-1 in the enhancement of cancer stem-like behaviors, including migration and invasion, in osteosarcoma, and identified the novel PVT-1/TRIM28 axis signaling cascade as a potential therapeutic target for osteosarcoma treatment.

K-Ras and p53 mouse model with molecular characteristics of human rhabdomyosarcoma and translational applications.

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in children, with overall long-term survival rates of ∼65-70%. Thus, additional molecular insights and representative models are critical for identifying and evaluating new treatment modalities. Using MyoD-Cre-mediated introduction of mutant K-RasG12D and perturbations in p53, we developed a novel genetically engineered mouse model (GEMM) for RMS. The anatomic sites of primary RMS development recapitulated human disease, including tumors in the head, neck, extremities and abdomen. We confirmed RMS histology and diagnosis through Hematoxylin and Eosin staining, and positive immunohistochemical staining for desmin, myogenin, and phosphotungstic acid-Hematoxylin. Cell lines from GEMM tumors were established with the ability to engraft in immunocompetent mice with comparable histological and staining features as the primary tumors. Tail vein injection of cell lines had high metastatic potential to the lungs. Transcriptomic analyses of p53R172H/K-RasG12D GEMM-derived tumors showed evidence of high molecular homology with human RMS. Finally, pre-clinical use of these murine RMS lines showed similar therapeutic responsiveness to chemotherapy and targeted therapies as human RMS cell lines.

MYCN-driven fatty acid uptake is a metabolic vulnerability in neuroblastoma.

Neuroblastoma (NB) is a childhood cancer arising from sympatho-adrenal neural crest cells. MYCN amplification is found in half of high-risk NB patients; however, no available therapies directly target MYCN. Using multi-dimensional metabolic profiling in MYCN expression systems and primary patient tumors, we comprehensively characterized the metabolic landscape driven by MYCN in NB. MYCN amplification leads to glycerolipid accumulation by promoting fatty acid (FA) uptake and biosynthesis. We found that cells expressing amplified MYCN depend highly on FA uptake for survival. Mechanistically, MYCN directly upregulates FA transport protein 2 (FATP2), encoded by SLC27A2. Genetic depletion of SLC27A2 impairs NB survival, and pharmacological SLC27A2 inhibition selectively suppresses tumor growth, prolongs animal survival, and exerts synergistic anti-tumor effects when combined with conventional chemotherapies in multiple preclinical NB models. This study identifies FA uptake as a critical metabolic dependency for MYCN-amplified tumors. Inhibiting FA uptake is an effective approach for improving current treatment regimens.

p21-activated kinases as viable therapeutic targets for the treatment of high-risk Ewing sarcoma.

Ewing sarcoma (ES) is the second most common bone tumor in children and young adults. Unfortunately, there have been minimal recent advancements in improving patient outcomes, especially in metastatic and recurrent diseases. In this study, we investigated the biological role of p21-activated kinases (PAKs) in ES, and the ability to therapeutically target them in high-risk disease. Via informatics analysis, we established the inverse association of PAK1 and PAK4 expression with clinical stage and outcome in ES patients. Through expression knockdown and small-molecule inhibition of PAKs, utilizing FRAX-597, KPT-9274, and PF-3758309 in multiple ES cell lines and patient-derived xenograft models, we further explored the role of PAKs in ES tumor growth and metastatic capabilities. In vitro studies in several ES cell lines indicated that diminishing PAK1 and PAK4 expression reduces tumor cell viability, migratory, and invasive properties. In vivo studies using PAK4 inhibitors, KPT-9274 and PF-3758309 demonstrated significant inhibition of primary and metastatic tumor formation, while transcriptomic analysis of PAK4-inhibitor-treated tumors identified concomitant suppression of Notch, ß-catenin, and hypoxia-mediated signatures. In addition, the analysis showed enrichment of anti-tumor immune regulatory mechanisms, including interferon (IFN)-É£ and IFN-α responses. Altogether, our molecular and pre-clinical studies are the first to establish a critical role for PAKs in ES development and progression, and consequently as viable therapeutic targets for the treatment of high-risk ES in the near future.