MYCN controls an alternative RNA splicing program in high-risk metastatic neuroblastoma.

The molecular mechanisms underlying the aggressive behavior of MYCN driven neuroblastoma (NBL) is under intense investigation; however, little is known about the impact of this family of transcription factors on the splicing program. Here we used high-throughput RNA sequencing to systematically study the expression of RNA isoforms in stage 4 MYCN-amplified NBL, an aggressive subtype of metastatic NBL. We show that MYCN-amplified NBL tumors display a distinct gene splicing pattern affecting multiple cancer hallmark functions. Six splicing factors displayed unique differential expression patterns in MYCN-amplified tumors and cell lines, and the binding motifs for some of these splicing factors are significantly enriched in differentially-spliced genes. Direct binding of MYCN to promoter regions of the splicing factors PTBP1 and HNRNPA1 detected by ChIP-seq demonstrates that MYCN controls the splicing pattern by direct regulation of the expression of these key splicing factors. Furthermore, high expression of PTBP1 and HNRNPA1 was significantly associated with poor overall survival of stage4 NBL patients (p ≤ 0.05). Knocking down PTBP1, HNRNPA1 and their downstream target PKM2, an isoform of pro-tumor-growth, result in repressed growth of NBL cells. Therefore, our study reveals a novel role of MYCN in controlling global splicing program through regulation of splicing factors in addition to its well-known role in the transcription program. These findings suggest a therapeutically potential to target the key splicing factors or gene isoforms in high-risk NBL with MYCN-amplification.

Comprehensive genomic analysis of rhabdomyosarcoma reveals a landscape of alterations affecting a common genetic axis in fusion-positive and fusion-negative tumors.

UNLABELLED: Despite gains in survival, outcomes for patients with metastatic or recurrent rhabdomyosarcoma remain dismal. In a collaboration between the National Cancer Institute, Children's Oncology Group, and Broad Institute, we performed whole-genome, whole-exome, and transcriptome sequencing to characterize the landscape of somatic alterations in 147 tumor/normal pairs. Two genotypes are evident in rhabdomyosarcoma tumors: those characterized by the PAX3 or PAX7 fusion and those that lack these fusions but harbor mutations in key signaling pathways. The overall burden of somatic mutations in rhabdomyosarcoma is relatively low, especially in tumors that harbor a PAX3/7 gene fusion. In addition to previously reported mutations in NRAS, KRAS, HRAS, FGFR4, PIK3CA, and CTNNB1, we found novel recurrent mutations in FBXW7 and BCOR, providing potential new avenues for therapeutic intervention. Furthermore, alteration of the receptor tyrosine kinase/RAS/PIK3CA axis affects 93% of cases, providing a framework for genomics-directed therapies that might improve outcomes for patients with rhabdomyosarcoma. SIGNIFICANCE: This is the most comprehensive genomic analysis of rhabdomyosarcoma to date. Despite a relatively low mutation rate, multiple genes were recurrently altered, including NRAS, KRAS, HRAS, FGFR4, PIK3CA, CTNNB1, FBXW7, and BCOR. In addition, a majority of rhabdomyosarcoma tumors alter the receptor tyrosine kinase/RAS/PIK3CA axis, providing an opportunity for genomics-guided intervention.

Massively parallel sequencing reveals an accumulation of de novo mutations and an activating mutation of LPAR1 in a patient with metastatic neuroblastoma.

Neuroblastoma is one of the most genomically heterogeneous childhood malignances studied to date, and the molecular events that occur during the course of the disease are not fully understood. Genomic studies in neuroblastoma have showed only a few recurrent mutations and a low somatic mutation burden. However, none of these studies has examined the mutations arising during the course of disease, nor have they systemically examined the expression of mutant genes. Here we performed genomic analyses on tumors taken during a 3.5 years disease course from a neuroblastoma patient (bone marrow biopsy at diagnosis, adrenal primary tumor taken at surgical resection, and a liver metastasis at autopsy). Whole genome sequencing of the index liver metastasis identified 44 non-synonymous somatic mutations in 42 genes (0.85 mutation/MB) and a large hemizygous deletion in the ATRX gene which has been recently reported in neuroblastoma. Of these 45 somatic alterations, 15 were also detected in the primary tumor and bone marrow biopsy, while the other 30 were unique to the index tumor, indicating accumulation of de novo mutations during therapy. Furthermore, transcriptome sequencing on the 3 tumors demonstrated only 3 out of the 15 commonly mutated genes (LPAR1, GATA2, and NUFIP1) had high level of expression of the mutant alleles, suggesting potential oncogenic driver roles of these mutated genes. Among them, the druggable G-protein coupled receptor LPAR1 was highly expressed in all tumors. Cells expressing the LPAR1 R163W mutant demonstrated a significantly increased motility through elevated Rho signaling, but had no effect on growth. Therefore, this study highlights the need for multiple biopsies and sequencing during progression of a cancer and combinatorial DNA and RNA sequencing approach for systematic identification of expressed driver mutations.