Glioblastoma (GBM) constitutes the most common and aggressive primary brain tumor. To better understand how GBM evolves we analyzed longitudinal genomic and transcriptomic data of 114 patients. The analysis reveals a highly branched evolutionary pattern in which 63% of patients experience expression-based subtype changes. The branching pattern together with estimates of evolutionary rates suggest that the relapse associated clone typically preexisted years before diagnosis. 15% of tumors present hypermutations at relapse in highly expressed genes with a clear mutational signature. We find that 11% of recurrent tumors harbor mutations in LTBP4, a protein binding to TGF-β. Silencing LTBP4 in GBM cells leads to TGF-β activity suppression and decreased proliferation. In IDH1-wild-type recurrent GBM, high LTBP4 expression is associated with worse prognosis, highlighting the TGF-β pathway as a potential therapeutic target in GBM.

Precision medicine in cancer proposes that genomic characterization of tumors can inform personalized targeted therapies1–5. This proposition, however, is complicated by spatial and temporal heterogeneity6–14. Here we study genomic and expression profiles across 127 multi-sector or longitudinal specimens from 52 glioblastoma (GBM) patients. Using bulk and single-cell data, we find that samples from the same tumor mass share genomic and expression signatures, while geographically separated multifocal tumors and/or long-term recurrent tumors are seeded from different clones. Chemical screening of patient-derived glioma cells (PDCs) shows that therapeutic response is associated to genetic similarity, and multifocal tumors enriched with PIK3CA mutations have a heterogeneous drug response pattern. Importantly, we show that targeting truncal events is more efficacious in reducing tumor burden. In summary, this work demonstrates that evolutionary inference from integrated genomic analysis in multi-sector biopsies can inform targeted therapeutic interventions for GBM patients.

No abstract

Pancreatic ductal adenocarcinoma (PDAC) is characterized by the presence of abundant desmoplastic stroma primarily composed of cancer-associated fibroblasts (CAFs). It is generally accepted that CAFs stimulate tumor progression and might be implicated in drug resistance and immunosuppression. Here, we have compared the transcriptional profile of PDGFRα CAFs isolated from genetically engineered mouse PDAC tumors with that of normal pancreatic fibroblasts to identify genes potentially implicated in their protumorigenic properties. We report that the most differentially expressed gene, , a member of the serum amyloid A (SAA) apolipoprotein family, is a key mediator of the protumorigenic activity of PDGFRα CAFs. Whereas -competent CAFs stimulate the growth of tumor cells in an orthotopic model,-null CAFs inhibit tumor growth. Saa3 also plays a role in the cross talk between CAFs and tumor cells. Ablation of in pancreatic tumor cells makes them insensitive to the inhibitory effect of-null CAFs. As a consequence, germline ablation of does not prevent PDAC development in mice. The protumorigenic activity of Saa3 in CAFs is mediated by Mpp6, a member of the palmitoylated membrane protein subfamily of the peripheral membrane-associated guanylate kinases (MAGUK). Finally, we interrogated whether these observations could be translated to a human scenario. Indeed,, the ortholog of murine Saa3, is overexpressed in human CAFs. Moreover, high levels of in the stromal component correlate with worse survival. These findings support the concept that selective inhibition of SAA1 in CAFs may provide potential therapeutic benefit to PDAC patients.

SF3B1, which encodes an essential spliceosomal protein, is frequently mutated in myelodysplastic syndromes (MDS) and many cancers. However, the defect of mutant SF3B1 is unknown. Here, we analyzed RNA sequencing data from MDS patients and confirmed that SF3B1 mutants use aberrant 3 0 splice sites. To elucidate the underlying mechanism, we purified complexes containing either wild-type or the hotspot K700E mutant SF3B1 and found that levels of a poorly studied spliceosomal protein, SUGP1, were reduced in mutant spliceosomes. Strikingly, SUGP1 knockdown completely recapitulated the splicing errors, whereas SUGP1 overexpression drove the protein, which our data suggest plays an important role in branchsite recognition, into the mutant spliceosome and partially rescued splicing. Other hotspot SF3B1 mutants showed similar altered splicing and diminished interaction with SUGP1. Our study demonstrates that SUGP1 loss is a common defect of spliceosomes with disease-causing SF3B1 mutations and, because this defect can be rescued, suggests possibilities for therapeutic intervention.(A) Distribution of cryptic 3 0 ss around the associated canonical 3 0 ss. Red and black histograms indicate 1,145 cryptic 3 0 ss more frequently used (q value < 0.05) in six mutant SF3B1 MDS patient samples and 186 cryptic 3 0 ss in nine WT SF3B1 samples, respectively. (B) Hierarchical clustering and heatmap analysis of the 627 cryptic 3 0 ss differentially used in six mutant versus nine WT SF3B1 samples (q value < 0.05 and cryptic 3 0 ss closer than 100 nt upstream of the associated canonical 3 0 ss). Each row represents one cryptic 3 0 ss and each column one MDS patient sample. Z scores in the matrix represent normalized percent-spliced-in values. The color bars above the heatmap indicate the SF3B1 mutations. (C) Volcano plot representation of genes associated with the 169 cryptic 3 0 ss differentially used in mutant versus WT SF3B1 samples (q value < 0.05, closer than 50 nt upstream of the associated canonical 3 0 ss and more than 15 supporting reads averaged over mutant SF3B1 samples). The horizontal axis shows the PSI difference between mutant and WT SF3B1 samples, and the vertical axis shows the significance. Genes selected for further experimental validation are highlighted in red.

Although mutations in the gene encoding the RNA splicing factor SF3B1 are frequent in multiple cancers, their functional effects and therapeutic dependencies are poorly understood. Here, we characterize 98 tumors and 12 isogenic cell lines harboring SF3B1 hotspot mutations, identifying hundreds of cryptic 3′ splice sites common and specific to different cancer types. Regulatory network analysis revealed that the most common SF3B1 mutation activates MYC via effects conserved across human and mouse cells. SF3B1 mutations promote decay of transcripts encoding the protein phosphatase 2A (PP2A) subunit PPP2R5A, increasing MYC S62 and BCL2 S70 phosphorylation which, in turn, promotes MYC protein stability and impair apoptosis, respectively. Genetic PPP2R5A restoration or pharmacologic PP2A activation impaired SF3B1-mutant tumorigenesis, elucidating a therapeutic approach to aberrant splicing by mutant SF3B1.SigNiFiCANCe: Here, we identify that mutations in SF3B1, the most commonly mutated splicing factor gene across cancers, alter splicing of a specific subunit of the PP2A serine/threonine phosphatase complex to confer post-translational MYC and BCL2 activation, which is therapeutically intervenable using an FDA-approved drug.

Here, we identify iPLA2β as a critical regulator for p53-driven ferroptosis upon reactive oxygen species (ROS)-induced stress. The calcium-independent phospholipase iPLA2β is known to cleave acyl tails from the glycerol backbone of lipids and release oxidized fatty acids from phospholipids. We found that iPLA2β-mediated detoxification of peroxidized lipids is sufficient to suppress p53-driven ferroptosis upon ROS-induced stress, even in GPX4-null cells. Moreover, iPLA2β is overexpressed in human cancers; inhibition of endogenous iPLA2β sensitizes tumor cells to p53-driven ferroptosis and promotes p53-dependent tumor suppression in xenograft mouse models. These results demonstrate that iPLA2β acts as a major ferroptosis repressor in a GPX4-independent manner. Notably, unlike GPX4, loss of iPLA2β has no obvious effect on normal development or cell viability in normal tissues but iPLA2β plays an essential role in regulating ferroptosis upon ROS-induced stress. Thus, our study suggests that iPLA2β is a promising therapeutic target for activating ferroptosis-mediated tumor suppression without serious toxicity concerns.

Plasmablastic lymphoma (PBL) is an aggressive B-cell non-Hodgkin lymphoma associated with immunodeficiency in the context of human immunodeficiency virus (HIV) infection or iatrogenic immunosuppression. While a rare disease in general, the incidence is dramatically increased in regions of the world with high HIV prevalence. The molecular pathogenesis of this disease is poorly characterized. Here, we defined the genomic features of PBL in a cohort of 110 patients from South Africa (15 by whole-exome sequencing and 95 by deep targeted sequencing). We identified recurrent mutations in genes of the JAK-STAT signaling pathway, including STAT3 (42%), JAK1 (14%), and SOCS1 (10%), leading to its constitutive activation. Moreover, 24% of cases harbored gain-of-function mutations in RAS family members (NRAS and KRAS). Comparative analysis with other B-cell malignancies uncovered PBL-specific somatic mutations and transcriptional programs. We also found recurrent copy number gains encompassing the CD44 gene (37%), which encodes for a cell surface receptor involved in lymphocyte activation and homing, and was found expressed at high levels in all tested cases, independent of genetic alterations. These findings have implications for the understanding of the pathogenesis of this disease and the development of personalized medicine approaches. SIGNIfICANCe: Plasmablastic lymphoma is a poorly studied and extremely aggressive tumor. Here we define the genomic landscape of this lymphoma in HIV-positive individuals from South Africa and identify pervasive mutations in JAK-STAT3 and RAS-MAPK signaling pathways. These data offer a genomic framework for the design of improved treatment strategies targeting these circuits.