Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation

Lysophospholipid acyltransferase-mediated formation of saturated glycerophospholipids maintained cell membrane integrity for hypoxic adaptation

Adaptation to hypoxia has attracted much public interest because of its clinical significance. However, hypoxic adaptation in the body is complicated and difficult to fully explore. To explore previously unknown conserved mechanisms and key proteins involved in hypoxic adaptation in different species, we first used a yeast model for mechanistic screening. Further multi-omics analyses in multiple species including yeast, zebrafish and mice revealed that glycerophospholipid metabolism was significantly involved in hypoxic adaptation with up-regulation of lysophospholipid acyltransferase (ALE1) in yeast, a key protein for the formation of dipalmitoyl phosphatidylcholine [DPPC (16:0/16:0)], which is a saturated phosphatidylcholine. Importantly, a mammalian homolog of ALE1, lysophosphatidylcholine acyltransferase 1 (LPCAT1), enhanced DPPC levels at the cell membrane and exhibited the same protective effect in mammalian cells under hypoxic conditions. DPPC supplementation effectively attenuated growth restriction, maintained cell membrane integrity and increased the expression of epidermal growth factor receptor under hypoxic conditions, but unsaturated phosphatidylcholine did not. In agreement with these findings, DPPC treatment could also repair hypoxic injury of intestinal mucosa in mice. Taken together, ALE1/LPCAT1-mediated DPPC formation, a key pathway of glycerophospholipid metabolism, is crucial for cell viability under hypoxic conditions. Moreover, we found that ALE1 was also involved in glycolysis to maintain sufficient survival conditions for yeast. The present study offers a novel approach to understanding lipid metabolism under hypoxia and provides new insights into treating hypoxia-related diseases.

The transition from genomics to phenomics in personalized population health

Modern health care faces several serious challenges, including an ageing population and its inherent burden of chronic diseases, rising costs and marginal quality metrics. By assessing and optimizing the health trajectory of each individual using a data-driven personalized approach that reflects their genetics, behaviour and environment, we can start to address these challenges. This assessment includes longitudinal phenome measures, such as the blood proteome and metabolome, gut microbiome composition and function, and lifestyle and behaviour through wearables and questionnaires. Here, we review ongoing large-scale genomics and longitudinal phenomics efforts and the powerful insights they provide into wellness. We describe our vision for the transformation of the current health care from disease-oriented to data-driven, wellness-oriented and personalized population health.

Oncogenic ETS fusions promote DNA damage and proinflammatory responses via pericentromeric RNAs in extracellular vesicles

Aberrant expression of the E26 transformation-specific (ETS) transcription factors characterizes numerous human malignancies. Many of these proteins, including EWS:FLI1 and EWS:ERG fusions in Ewing sarcoma (EwS) and TMPRSS2:ERG in prostate cancer (PCa), drive oncogenic programs via binding to GGAA repeats. We report here that both EWS:FLI1 and ERG bind and transcriptionally activate GGAA-rich pericentromeric heterochromatin. The respective pathogen-like HSAT2 and HSAT3 RNAs, together with LINE, SINE, ERV, and other repeat transcripts, are expressed in EwS and PCa tumors, secreted in extracellular vesicles (EVs), and are highly elevated in plasma of patients with EwS with metastatic disease. High human satellite 2 and 3 (HSAT2,3) levels in EWS:FLI1- or ERG-expressing cells and tumors were associated with induction of G2/M checkpoint, mitotic spindle, and DNA damage programs. These programs were also activated in EwS EV-treated fibroblasts, coincident with accumulation of HSAT2,3 RNAs, proinflammatory responses, mitotic defects, and senescence. Mechanistically, HSAT2,3-enriched cancer EVs induced cGAS-TBK1 innate immune signaling and formation of cytosolic granules positive for double-strand RNAs, RNA-DNA, and cGAS. Hence, aberrantly expressed ETS proteins derepress pericentromeric heterochromatin, yielding pathogenic RNAs that transmit genotoxic stress and inflammation to local and distant sites. Monitoring HSAT2,3 plasma levels and preventing their dissemination may thus improve therapeutic strategies and blood-based diagnostics.

Surface and Global Proteome Analyses Identify ENPP1 and Other Surface Proteins as Actionable Immunotherapeutic Targets in Ewing Sarcoma

PURPOSE

Ewing sarcoma is the second most common bone sarcoma in children, with 1 case per 1.5 million in the United States. Although the survival rate of patients diagnosed with localized disease is approximately 70%, this decreases to approximately 30% for patients with metastatic disease and only approximately 10% for treatment-refractory disease, which have not changed for decades. Therefore, new therapeutic strategies are urgently needed for metastatic and refractory Ewing sarcoma.

EXPERIMENTAL DESIGN

This study analyzed 19 unique Ewing sarcoma patient- or cell line-derived xenografts (from 14 primary and 5 metastatic specimens) using proteomics to identify surface proteins for potential immunotherapeutic targeting. Plasma membranes were enriched using density gradient ultracentrifugation and compared with a reference standard of 12 immortalized non-Ewing sarcoma cell lines prepared in a similar manner. In parallel, global proteome analysis was carried out on each model to complement the surfaceome data. All models were analyzed by Tandem Mass Tags-based mass spectrometry to quantify identified proteins.

RESULTS

The surfaceome and global proteome analyses identified 1,131 and 1,030 annotated surface proteins, respectively. Among surface proteins identified, both approaches identified known Ewing sarcoma-associated proteins, including IL1RAP, CD99, STEAP1, and ADGRG2, and many new cell surface targets, including ENPP1 and CDH11. Robust staining of ENPP1 was demonstrated in Ewing sarcoma tumors compared with other childhood sarcomas and normal tissues.

CONCLUSIONS

Our comprehensive proteomic characterization of the Ewing sarcoma surfaceome provides a rich resource of surface-expressed proteins in Ewing sarcoma. This dataset provides the preclinical justification for exploration of targets such as ENPP1 for potential immunotherapeutic application in Ewing sarcoma. See related commentary by Bailey, p. 934.

Current updates in sarcoma biomarker discovery: emphasis on next-generation sequencing-based methods

Soft tissue sarcomas comprise a heterogeneous group of neoplasms. Although soft tissue malignancies make up only 2% of adult cancers, classification based on histomorphology presents a diagnostic challenge. Characterisation of soft tissue sarcomas by molecular analysis is rapidly evolving to improve diagnostic accuracy and develop targeted therapies. This review highlights the advances in molecular techniques, including current next-generation sequencing-based assays (fusion detection by RNA sequencing, targeted/whole exome sequencing, microRNA profiling), as well as emerging methods (liquid biopsies, DNA methylation profiling, single-cell molecular profiling and next-generation immunohistochemistry) for future clinical applications.

Single cell RNA-sequencing of Ewing sarcoma tumors demonstrates transcriptional heterogeneity and clonal evolution

Ewing sarcoma is the second most common bone cancer in children, accounting for 2% of pediatric cancer diagnoses. Patients who present with metastatic disease at the time of diagnosis have a dismal prognosis, compared to the >70% 5-year survival of those with localized disease. Here, we utilized single cell RNA-sequencing to characterize the transcriptional landscape of primary Ewing sarcoma tumors and surrounding tumor microenvironment (TME). Copy-number analysis identified subclonal evolution within patients even prior to treatment. Primary tumor samples demonstrate a heterogenous transcriptional landscape with several conserved gene expression programs, including those composed of genes related to proliferation and EWS targets. We also were able to identify the composition of the TME and molecularly dissect the transcriptional profile of circulating tumor cells in peripheral blood at the time of diagnosis.

Genetic analyses identify evidence for a causal relationship between Ewing sarcoma and hernias

Knowledge of Ewing sarcoma (EWS) risk factors is exceedingly limited; however, multiple small, independent studies have suggested a possible connection between hernia and EWS. By leveraging hernia summary statistics from the UK Biobank and a recently published genome-wide association study of EWS (733 EWS cases and 1,346 controls), we conducted a genetic investigation of the relationship of 5 hernia types (diaphragmatic, inguinal, umbilical, femoral, and ventral) and EWS. We discovered a positive causal relationship between inguinal hernia and EWS (OR 1.27, 95% confidence interval [CI] 1.01-1.59, and p = 0.041) through Mendelian randomization analysis. Further analyses suggested shared pathways through three genes: HMGA2, LOX, and FBXW7. Diaphragmatic hernia showed a stronger causal relationship with EWS among all of the hernia types (OR 2.26, 95% CI 1.30-3.95, p = 0.004), but no statistically significant local correlation pattern was observed. No evidence of a causal or genetic relationship was observed between EWS and the other three hernia types, including umbilical hernia, despite a previous report indicating an OR as high as 3.3. The finding of our genetic analysis provided additional support to the hypothesis that EWS and hernias may share a common origin.

Illuminating phenotypic drug responses of sarcoma cells to kinase inhibitors by phosphoproteomics

Kinase inhibitors (KIs) are important cancer drugs but often feature polypharmacology that is molecularly not understood. This disconnect is particularly apparent in cancer entities such as sarcomas for which the oncogenic drivers are often not clear. To investigate more systematically how the cellular proteotypes of sarcoma cells shape their response to molecularly targeted drugs, we profiled the proteomes and phosphoproteomes of 17 sarcoma cell lines and screened the same against 150 cancer drugs. The resulting 2550 phenotypic profiles revealed distinct drug responses and the cellular activity landscapes derived from deep (phospho)proteomes (9-10,000 proteins and 10-27,000 phosphorylation sites per cell line) enabled several lines of analysis. For instance, connecting the (phospho)proteomic data with drug responses revealed known and novel mechanisms of action (MoAs) of KIs and identified markers of drug sensitivity or resistance. All data is publicly accessible via an interactive web application that enables exploration of this rich molecular resource for a better understanding of active signalling pathways in sarcoma cells, identifying treatment response predictors and revealing novel MoA of clinical KIs.

ETS1, a target gene of the EWSR1::FLI1 fusion oncoprotein, regulates the expression of the focal adhesion protein TENSIN3

The mechanistic basis for the metastasis of Ewing sarcomas remains poorly understood, as these tumors harbor few mutations beyond the chromosomal translocation that initiates the disease. Instead, the epigenome of Ewing sarcoma (EWS) cells reflects the regulatory state of genes associated with the DNA binding activity of the fusion oncoproteins EWSR1::FLI1 or EWSR1::ERG. In this study, we examined the EWSR1::FLI1/ERG's repression of transcription factor genes, concentrating on those that exhibit a broader range of expression in tumors than in EWS cell lines. Focusing on one of these target genes, ETS1, we detected EWSR1::FLI1 binding and an H3K27me3 repressive mark at this locus. Depletion of EWSR1::FLI1 results in ETS1's binding of promoter regions, substantially altering the transcriptome of EWS cells, including the upregulation of the gene encoding TENSIN3 (TNS3), a focal adhesion protein. EWS cell lines expressing ETS1 (CRISPRa) exhibited increased TNS3 expression and enhanced movement compared to control cells. The cytoskeleton of control cells and ETS1-activated EWS cell lines also differed. Specifically, control cells exhibited a distributed vinculin signal and a network-like organization of F-actin. In contrast, ETS1-activated EWS cells showed an accumulation of vinculin and F-actin towards the plasma membrane. Interestingly, the phenotype of ETS1-activated EWS cell lines depleted of TNS3 resembled the phenotype of the control cells. Critically, these findings have clinical relevance as TNS3 expression in EWS tumors positively correlates with that of ETS1.

Discovering chromatin dysregulation induced by protein-coding perturbations at scale

Although population-scale databases have expanded to millions of protein-coding variants, insight into variant mechanisms has not kept pace. We present PROD-ATAC, a high-throughput method for discovering the effects of protein-coding variants on chromatin. A pooled library of variants is expressed in a disease-agnostic cell line, and single-cell ATAC resolves each variant's effect on chromatin. Using PROD-ATAC, we characterized the effects of >100 oncofusions (a class of cancer-causing chimeric proteins) and controls and revealed that pioneer activity is a common feature of fusions spanning an enormous range of fusion frequencies. Further, fusion-induced dysregulation can be context-agnostic as observed mechanisms often overlapped with cancer and cell-type specific prior knowledge. We also showed that gain-of-function pioneering is common among oncofusions. This work provides a global view of fusion-induced chromatin. We uncovered convergent mechanisms among disparate oncofusions and shared modes of dysregulation across different cancers. PROD-ATAC is generalizable to any set of protein-coding variants.

YBX1-interacting small RNAs and RUNX2 can be blocked in primary bone cancer using CADD522

Primary bone cancer (PBC) comprises several subtypes each underpinned by distinctive genetic drivers. This driver diversity produces novel morphological features and clinical behaviour that serendipitously makes PBC an excellent metastasis model. Here, we report that some transfer RNA-derived small RNAs termed tRNA fragments (tRFs) perform as a constitutive tumour suppressor mechanism by blunting a potential pro-metastatic protein-RNA interaction. This mechanism is reduced in PBC progression with a gradual loss of tRNAGly^TCC cleavage into 5' end tRF-Gly^TCC when comparing low-grade, intermediate-grade and high-grade patient tumours. We detected recurrent activation of miR-140 leading to upregulated RUNX2 expression in high-grade patient tumours. Both tRF-Gly^TCC and RUNX2 share a sequence motif in their 3' ends that matches the YBX1 recognition site known to stabilise pro-metastatic mRNAs. Investigating some aspects of this interaction network, gain- and loss-of-function experiments using small RNA mimics and antisense LNAs, respectively, showed that ectopic tRF-Gly^TCC reduced RUNX2 expression and dispersed 3D micromass architecture in vitro. iCLIP sequencing revealed YBX1 physical binding to the 3' UTR of RUNX2. The interaction between YBX1, tRF-Gly^TCC and RUNX2 led to the development of the RUNX2 inhibitor CADD522 as a PBC treatment. CADD522 assessment in vitro revealed significant effects on PBC cell behaviour. In xenograft mouse models, CADD522 as a single agent without surgery significantly reduced tumour volume, increased overall and metastasis-free survival and reduced cancer-induced bone disease. Our results provide insight into PBC molecular abnormalities that have led to the identification of new targets and a new therapeutic.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.