Background/Objectives: Osteosarcoma is the most common malignant bone tumor in children, characterized by a high degree of genomic instability, resulting in copy number alterations and genomic rearrangements without disease-defining recurrent mutations. Clinical trials based on molecular characterization have failed to find new effective therapies or improve outcomes over the last 40 years. Methods: To better understand the immune microenvironment of osteosarcoma, we performed single-cell RNA sequencing on six tumor biopsy samples, combined with a previously published cohort of six samples. Additional osteosarcoma samples were profiled using spatial transcriptomics for the validation of discovered subtypes and to add spatial context. Results: Analysis revealed immunosuppressive cells, including myeloid-derived suppressor cells (MDSCs), regulatory and exhausted T cells, and LAMP3+ dendritic cells. Conclusions: Using cell–cell communication modeling, we identified robust interactions between MDSCs and other cells, leading to NF-κB upregulation and an immunosuppressive microenvironment, as well as interactions involving regulatory T cells and osteosarcoma cells that promoted tumor progression and a proangiogenic niche.

• In vitro 3D medulloblastoma models have indistinguishable transcriptomic profiles than the original tumor tissue. • One type of chemically defined media was identified to support four major types of pediatric brain tumor better than serum-containing culture conditions. • Different tumor types show varying spheroid-forming ability depending on stem cell composition in the original tumor cell population. • Different tumor types show different preference to 3D scaffold and/or extracellular matrix for the in vitro growth accompanied by microenvironment-induced gene expression changes.

Supratentorial ependymoma (ST-EPN) is a type of malignant brain tumor mainly seen in children. Since 2014, it has been known that an intrachromosomal fusion C11orf95-RELA is an oncogenic driver in ST-EPN [Parker et al. Nature 506:451–455 (2014); Pietsch et al. Acta Neuropathol 127:609–611 (2014)] but the molecular mechanisms of oncogenesis are unclear. Here we show that the C11orf95 component of the fusion protein dictates DNA binding activity while the RELA component is required for driving the expression of ependymoma-associated genes. Epigenomic characterizations using ChIP-seq and HiChIP approaches reveal that C11orf95-RELA modulates chromatin states and mediates chromatin interactions, leading to transcriptional reprogramming in ependymoma cells. Our findings provide important characterization of the molecular underpinning of C11orf95-RELA fusion and shed light on potential therapeutic targets for C11orf95-RELA subtype ependymoma.

Genome‐wide association studies (GWAS) have identified single‐nucleotide polymorphisms (SNPs) associated with glioma risk on 20q13.33, but the biological mechanisms underlying this association are unknown. We tested the hypothesis that a functional SNP on 20q13.33 impacted the activity of an enhancer, leading to an altered expression of nearby genes. To identify candidate functional SNPs, we identified all SNPs in linkage disequilibrium with the risk‐associated SNP rs2297440 that mapped to putative enhancers. Putative enhancers containing candidate functional SNPs were tested for allele‐specific effects in luciferase enhancer activity assays against glioblastoma multiforme (GBM) cell lines. An enhancer containing SNP rs3761124 exhibited allele‐specific effects on activity. Deletion of this enhancer by CRISPR‐Cas9 editing in GBM cell lines correlated with an altered expression of multiple genes, including STMN3, RTEL1, RTEL1‐TNFRSF6B, GMEB2, and SRMS. Expression quantitative trait loci (eQTL) analyses using nondiseased brain samples, isocitrate dehydrogenase 1 (IDH1) wild‐type glioma, and neurodevelopmental tissues showed STMN3 to be a consistent significant eQTL with rs3761124. RTEL1 and GMEB2 were also significant eQTLs in the context of early CNS development and/or in IDH1 wild‐type glioma. We provide evidence that rs3761124 is a functional variant on 20q13.33 related to glioma/GBM risk that modulates the expression of STMN3 and potentially other genes across diverse cellular contexts.

Osteosarcoma is a primary malignant bone tumor arising from bone-forming mesenchymal cells in children and adolescents. Despite efforts to understand the biology of the disease and identify novel therapeutics, the survival of osteosarcoma patients remains dismal. We have concurrently profiled the copy number and gene expression of 226 osteosarcoma samples as part of the Strategic Partnering to Evaluate Cancer Signatures (SPECS) initiative. Our results demonstrate the heterogeneous landscape of osteosarcoma in younger populations by showing the presence of genome-wide copy number abnormalities occurring both recurrently among samples and in a high frequency. Insulin growth factor receptor 1 (IGF1R) is a receptor tyrosine kinase which binds IGF1 and IGF2 to activate downstream pathways involved in cell apoptosis and proliferation. We identify prevalent amplification of IGF1R corresponding with increased gene expression in patients with poor survival outcomes. Our results substantiate previously tenuously associated copy number abnormalities identified in smaller datasets (13q34+, 20p13+, 4q35-, 20q13.33-), and indicate the significance of high fibroblast growth factor receptor 2 (FGFR2) expression in distinguishing patients with poor prognosis. FGFR2 is involved in cellular proliferation processes such as division, growth and angiogenesis. In summary, our findings demonstrate the prognostic significance of several genes associated with osteosarcoma pathogenesis.

Although familial susceptibility to glioma is known, the genetic basis for this susceptibility remains unidentified in the majority of glioma-specific families. An alternative approach to identifying such genes is to examine cancer pedigrees, which include glioma as one of several cancer phenotypes, to determine whether common chromosomal modifications might account for the familial aggregation of glioma and other cancers. Germline rearrangements in 146 glioma families (from the Gliogene Consortium; http://www.gliogene.org/) were examined using multiplex ligation-dependent probe amplification. These families all had at least 2 verified glioma cases and a third reported or verified glioma case in the same family or 2 glioma cases in the family with at least one family member affected with melanoma, colon, or breast cancer.The genomic areas covering TP53, CDKN2A, MLH1, and MSH2 were selected because these genes have been previously reported to be associated with cancer pedigrees known to include glioma. We detected a single structural rearrangement, a deletion of exons 1-6 in MSH2, in the proband of one family with 3 cases with glioma and one relative with colon cancer. Large deletions and duplications are rare events in familial glioma cases, even in families with a strong family history of cancers that may be involved in known cancer syndromes.