Search

Interrogation of Oncomine{trade mark, serif} (https://www.oncomine.org/) showed higher expression of eIF4E and eIF5 in breast (top panel) and lung cancer (bottom panel) compared to normal tissue. When these biomarkers were compared for gender, eIF4E and eIF5 expression was proportionately higher in male breast cancer (A, C) but not in lung cancer (B, D), respectively .

Purpose: Breast cancer affects both genders, but is understudied in men. Although still rare, male breast cancer (MBC) is being diagnosed more frequently. Treatments are wholly informed by clinical studies conducted in women, based on assumptions that underlying biology is similar.Experimental Design: A transcriptomic investigation of male and female breast cancer was performed, confirming transcriptomic data in silico. Biomarkers were immunohistochemically assessed in 697 MBCs (n = 477, training; n = 220, validation set) and quantified in pre- and posttreatment samples from an MBC patient receiving everolimus and PI3K/mTOR inhibitor.Results: Gender-specific gene expression patterns were identified. eIF transcripts were upregulated in MBC. eIF4E and eIF5 were negatively prognostic for overall survival alone (log-rank P = 0.013; HR = 1.77, 1.12–2.8 and P = 0.035; HR = 1.68, 1.03–2.74, respectively), or when coexpressed (P = 0.01; HR = 2.66, 1.26–5.63), confirmed in the validation set. This remained upon multivariate Cox regression analysis [eIF4E P = 0.016; HR = 2.38 (1.18–4.8), eIF5 P = 0.022; HR = 2.55 (1.14–5.7); coexpression P = 0.001; HR = 7.04 (2.22–22.26)]. Marked reduction in eIF4E and eIF5 expression was seen post BEZ235/everolimus, with extended survival.Conclusions: Translational initiation pathway inhibition could be of clinical utility in MBC patients overexpressing eIF4E and eIF5. With mTOR inhibitors that target this pathway now in the clinic, these biomarkers may represent new targets for therapeutic intervention, although further independent validation is required. Clin Cancer Res; 23(10); 2575–83. ©2016 AACR.

R.O.C curve analysis was used to determine the cut-off value for each antibody to determine a point of dichotomisation (graphs). Gradation of intensity of immunohistochemical staining (weak, moderate, strong) for eIF1 (A), eIF2 (B), eIF4E (C) and eIF5 (D) is shown in representative MBC TMA cores. Immunoreactivity is shown to the right of each R.O.C curve generated. Each core was individually scored with the proportion of tumour cells stained 0 - 4, followed by the intensity of staining 0 - 3. Combination scores reported were in the range 0 - 7. The Interclass Correlation Coefficient from scorer concordance of each biomarker was; eIF1 0.911 (95%CI 0.769-0.944), eIF2 0.846 (95%CI 0.736-0.910), eIF4E 0.882 (95%CI 0.755-0.913), eIF5 0.865 (95%CI 0.769-0.922). No staining was observed for eIF3 in MBC samples, although this was expressed in positive control tissue for EIF3, colon (E).

Details of the male breast cancer samples used at different stages of the experimental workflow

Relationship between eIF4E, eIF5 and clinico-pathological parameters

Glioblastoma (GBM) brain tumours lackingIDH1mutations (IDHwt) have the worst prognosis of all brain neoplasms. Patients receive surgery and chemoradiotherapy but tumours almost always fatally recur. Using RNAseq data from 107 pairs of pre- and post-standard treatment locally recurrent IDHwt GBM tumours, we identified two responder subtypes based on therapy-driven changes in gene expression. In two thirds of patients a specific subset of genes is up-regulated from primary to recurrence (Up responders) and in one third the same genes are down-regulated (Down responders). Characterisation of the responder subtypes indicates subtype-specific adaptive treatment resistance mechanisms. In Up responders treatment enriches for quiescent proneural GBM stem cells and differentiated neoplastic cells with increased neurotransmitter signalling, whereas Down responders commonly undergo therapy-driven mesenchymal transition. Stratifying GBM tumours by response subtype may lead to more effective treatment. In support of this, modulators of gamma aminobutyric acid (GABA) neurotransmitter signalling differentially sensitise Up and Down responder GBM models to standard treatmentin vitro.

MicroRNAs play an important role in the regulation of mRNA translation and have therapeutic potential in cancer and other diseases. To profile the landscape of microRNAs with significant cytotoxicity in the context of glioblastoma (GBM), we performed a high-throughput screen in adult and pediatric GBM cells using a synthetic oligonucleotide library representing all known human microRNAs. Bioinformatics analysis was used to refine this list and the top seven microRNAs were validated in a larger panel of GBM cells using state-of-the-art in vitro assays. The cytotoxic effect of our most relevant candidate was assessed in a preclinical model. Our screen identified ~100 significantly cytotoxic microRNAs with 70% concordance between cell lines. MicroRNA-1300 (miR-1300) was the most potent and robust candidate. We observed a striking binucleated phenotype in miR-1300 transfected cells due to cytokinesis failure followed by apoptosis. This was also observed in two stem-like patient-derived cultures. We identified the physiological role of miR-1300 as a regulator of endomitosis in megakaryocyte differentiation where blockade of cytokinesis is an essential step. In GBM cells, where miR-1300 is normally not expressed, the oncogene Epithelial Cell Transforming 2 (ECT2) was validated as a direct key target. ECT2 siRNA phenocopied the effects of miR-1300, and ECT2 overexpression led to rescue of miR-1300 induced binucleation. We showed that ectopic expression of miR-1300 led to decreased tumor growth in an orthotopic GBM model. Our screen provides a resource for the neuro-oncology community and identified miR-1300 as a novel regulator of endomitosis with translatable potential for therapeutic application.

Purpose The development of biomarkers and molecularly targeted therapies for patients with Ewing sarcoma (ES) in order to minimise morbidity and improve outcome is urgently needed. Here, we set out to isolate and characterise patient-derived ES primary cell cultures and daughter cancer stem-like cells (CSCs) to identify biomarkers of high-risk disease and candidate therapeutic targets. Methods Thirty-two patient-derived primary cultures were established from treatment-naïve tumours and primary ES-CSCs isolated from these cultures using functional methods. By RNA-sequencing we analysed the transcriptome of ES patient-derived cells (n = 24) and ES-CSCs (n = 11) to identify the most abundant and differentially expressed genes (DEGs). Expression of the top DEG(s) in ES-CSCs compared to ES cells was validated at both RNA and protein levels. The functional and prognostic potential of the most significant gene (neurexin-1) was investigated using knock-down studies and immunohistochemistry of two independent tumour cohorts. Results ES-CSCs were isolated from all primary cell cultures, consistent with the premise that ES is a CSC driven cancer. Transcriptional profiling confirmed that these cells were of mesenchymal origin, revealed novel cell surface targets for therapy that regulate cell-extracellular matrix interactions and identified candidate drivers of progression and relapse. High expression of neurexin-1 and low levels of regulators of its activity, APBA1 and NLGN4X, were associated with poor event-free and overall survival rates. Knock-down of neurexin-1 decreased viable cell numbers and spheroid formation. Conclusions Genes that regulate extracellular interactions, including neurexin-1, are candidate therapeutic targets in ES. High levels of neurexin-1 at diagnosis are associated with poor outcome and identify patients with localised disease that will relapse. These patients could benefit from more intensive or novel treatment modalities. The prognostic significance of neurexin-1 should be validated independently.

AIMS Glioblastoma (GBM) is the most common primary malignant brain tumour in adults. Despite aggressive treatment, a resistant tumour recurs in practically all patients. We therefore aimed to better understand the mechanisms driving this treatment resistance through investigating changes in gene expression across pairs of primary and recurrent GBM tumours. METHOD We generated or acquired bulk tumour RNA sequencing data for primary and first recurrent tumours from 107 patients who received standard treatment. Differential expression analysis between primary and recurrent samples found that the most dysregulated genes were involved in neurodevelopment and neurodifferentiation. We therefore used a publicly available ChIP-seq database to identify DNA binding factors for which binding sites are enriched in the promotors of genes with the largest expression changes from primary to recurrent. RESULTS Jumonji and AT-Rich Interacting Domain 2 (JARID2) was the most strongly enriched for binding to promotors of dysregulated genes. 65 patients showed an up-regulation and 42 showed a down-regulation of genes bound by this protein. The same set of JARID2 bound genes were found to be dysregulated in each direction, and correlated with the largest source of variation between samples in their response to treatment. Further enrichment analyses indicated that ‘Up’ responders may resist treatment through reduced proliferation and increased interaction with the tumour microenvironment, whereas ‘Down’ responders instead rely on a shift to mesenchymal cell states. CONCLUSION These results indicate that GBM tumours can be split into two subtypes that transcriptionally reprogramme in different directions through treatment and may benefit from different treatment approaches.

Summary Glioblastoma (GBM) is an aggressive cancer with a very poor prognosis. Generally viewed as weakly immunogenic, GBM responds poorly to current immunotherapies. To understand this problem more clearly we used a combination of natural killer (NK) cell functional assays together with gene and protein expression profiling to define the NK cell response to GBM and explore immunosuppression in the GBM microenvironment. In addition, we used transcriptome data from patient cohorts to classify GBM according to immunological profiles. We show that glioma stem‐like cells, a source of post‐treatment tumour recurrence, express multiple immunomodulatory cell surface molecules and are targeted in preference to normal neural progenitor cells by natural killer (NK) cells ex vivo. In contrast, GBM‐infiltrating NK cells express reduced levels of activation receptors within the tumour microenvironment, with hallmarks of transforming growth factor (TGF)‐β‐mediated inhibition. This NK cell inhibition is accompanied by expression of multiple immune checkpoint molecules on T cells. Single‐cell transcriptomics demonstrated that both tumour and haematopoietic‐derived cells in GBM express multiple, diverse mediators of immune evasion. Despite this, immunome analysis across a patient cohort identifies a spectrum of immunological activity in GBM, with active immunity marked by co‐expression of immune effector molecules and feedback inhibitory mechanisms. Our data show that GBM is recognized by the immune system but that anti‐tumour immunity is restrained by multiple immunosuppressive pathways, some of which operate in the healthy brain. The presence of immune activity in a subset of patients suggests that these patients will more probably benefit from combination immunotherapies directed against multiple immunosuppressive pathways., Our manuscript challenges the current paradigm that Glioblastoma (GBM) is a poor candidate for immunotherapy, therefore impacting on clinical approaches. We identify NK cells as a target for immunosuppression in GBM, with multiple suppressive pathways within GBM human tumours. We also identify a group of GBM patients with higher immunological activity, who will benefit from immunotherapies that target these immunosuppressive pathways.

Low Temperature Plasma (LTP) generates reactive oxygen and nitrogen species, causing cell death, similarly to radiation. Radiation resistance results in tumour recurrence, however mechanisms of LTP resistance are unknown. LTP was applied to patient-derived prostate epithelial cells and gene expression assessed. A typical global oxidative response (AP-1 and Nrf2 signalling) was induced, whereas Notch signalling was activated exclusively in progenitor cells. Notch inhibition induced expression of prostatic acid phosphatase (PAP), a marker of prostate epithelial cell differentiation, whilst reducing colony forming ability and preventing tumour formation. Therefore, if LTP is to be progressed as a novel treatment for prostate cancer, combination treatments should be considered in the context of cellular heterogeneity and existence of cell type-specific resistance mechanisms.

Background Brain metastases (BrM) develop in 20–40% of cancer patients and represent an unmet clinical need. Limited access of drugs into the brain because of the blood-brain barrier is at least partially responsible for therapeutic failure, necessitating improved drug delivery systems. Methods Green fluorescent protein (GFP)-transduced murine and nontransduced human hematopoietic stem cells (HSCs) were administered into mice (n = 10 and 3). The HSC progeny in mouse BrM and in patient-derived BrM tissue (n = 6) was characterized by flow cytometry and immunofluorescence. Promoters driving gene expression, specifically within the BrM-infiltrating HSC progeny, were identified through differential gene-expression analysis and subsequent validation of a series of promoter-green fluorescent protein-reporter constructs in mice (n = 5). One of the promoters was used to deliver tumor necrosis factor–related apoptosis-inducing ligand (TRAIL) to BrM in mice (n = 17/21 for TRAIL vs control group). Results HSC progeny (consisting mostly of macrophages) efficiently homed to macrometastases (mean [SD] = 37.6% [7.2%] of all infiltrating cells for murine HSC progeny; 27.9% mean [SD] = 27.9% [4.9%] of infiltrating CD45+ hematopoietic cells for human HSC progeny) and micrometastases in mice (19.3–53.3% of all macrophages for murine HSCs). Macrophages were also abundant in patient-derived BrM tissue (mean [SD] = 8.8% [7.8%]). Collectively, this provided a rationale to optimize the delivery of gene therapy to BrM within myeloid cells. MMP14 promoter emerged as the strongest promoter construct capable of limiting gene expression to BrM-infiltrating myeloid cells in mice. TRAIL delivered under MMP14 promoter statistically significantly prolonged survival in mice (mean [SD] = 19.0 [3.4] vs mean [SD] = 15.0 [2.0] days for TRAIL vs control group; two-sided P = .006), demonstrating therapeutic and translational potential of our approach. Conclusions Our study establishes HSC gene therapy using a myeloid cell–specific promoter as a new strategy to target BrM. This approach, with strong translational value, has potential to overcome the blood-brain barrier, target micrometastases, and control multifocal lesions.

The authors wish to make the following corrections to this paper [1] [...]

Dedifferentiation and transdifferentiation are rare and only poorly understood phenomena in cutaneous melanoma. To study this disease more comprehensively we have retrieved 11 primary cutaneous melanomas from our pathology archives showing biphasic features characterized by a conventional melanoma and additional areas of de-/trans-differentiation as defined by a lack of immunohistochemical expression of all conventional melanocytic markers (S-100 protein, SOX10, Melan- A, and HMB-45). The clinical, histologic, and immunohistochemical findings were recorded and follow-up was obtained. The patients were mostly elderly (median: 81 years ; range: 42–86 years) without significant gender predilection, and the sun-exposed skin of the head and neck area was most commonly affected. The tumors were deeply invasive with a mean depth of 7 mm (range: 4–80 mm). The dedifferentiated component showed atypical fibroxanthoma-like features in the majority of cases (7), while additional rhabdomyosarcomatous and epithelial transdifferentiation was noted histologically and/or immunohistochemically in two tumors each. The background conventional melanoma component was of desmoplastic (4), superficial spreading (3), nodular (2), lentigo maligna (1), or spindle cell (1) types. For the seven patients with available follow-up data (median follow-up period of 25 months ; range: 8–36 months), two died from their disease, and three developed metastases. Next- generation sequencing of the cohort revealed somatic mutations of established melanoma drivers including mainly NF1 mutations (5) in the conventional component, which was also detected in the corresponding de-/trans-differentiated component. In summary, the diagnosis of primary cutaneous de-/trans-differentiated melanoma is challenging and depends on the morphologic identification of conventional melanoma. Molecular analysis is diagnostically helpful as the mutated gene profile is shared between the conventional and de-/trans-differentiated components. Importantly, de-/trans-differentiation does not appear to confer a more aggressive behavior.

While oncogenes promote tumorigenesis, they also induce deleterious cellular stresses, such as apoptosis, that cancer cells must combat by coopting adaptive responses. Whether tumor suppressor gene haploinsufficiency leads to such phenomena and their mechanistic basis is unclear. Here, we demonstrate that elevated levels of the anti-apoptotic factor, CASP8 and FADD-like apoptosis regulator (CFLAR), promotes apoptosis evasion in acute myeloid leukemia (AML) cells haploinsufficient for the cut-like homeobox 1 (CUX1) transcription factor, whose loss is associated with dismal clinical prognosis. Genome-wide CRISPR/Cas9 screening identifies CFLAR as a selective, acquired vulnerability in CUX1-deficient AML, which can be mimicked therapeutically using inhibitor of apoptosis (IAP) antagonists in murine and human AML cells. Mechanistically, CUX1 deficiency directly alleviates CUX1 repression of the CFLAR promoter to drive CFLAR expression and leukemia survival. These data establish how haploinsufficiency of a tumor suppressor is sufficient to induce advantageous anti-apoptosis cell survival pathways and concurrently nominate CFLAR as potential therapeutic target in these poor-prognosis leukemias., Cut-like homeobox 1 (CUX1) is a haploinsufficient tumor suppressor commonly inactivated in acute myeloid leukemia and high-risk myelodysplasia. Here, in a genetically modified murine model, the authors show that CUX1 deficiency impairs apoptosis leading to leukemia when combined with mutant Flt3-ITD.

Melanoma represents ~5% of all cutaneous malignancies, yet accounts for the majority of skin cancer deaths due to its propensity to metastasise. To develop new therapies, novel target molecules must to be identified and the accessibility of cell surface proteins makes them attractive targets. Using CRISPR activation technology, we screened a library of guide RNAs targeting membrane protein-encoding genes to identify cell surface molecules whose upregulation enhances the metastatic pulmonary colonisation capabilities of tumour cells in vivo. We show that upregulated expression of the cell surface protein LRRN4CL led to increased pulmonary metastases in mice. Critically, LRRN4CL expression was elevated in melanoma patient samples, with high expression levels correlating with decreased survival. Collectively, our findings uncover an unappreciated role for LRRN4CL in the outcome of melanoma patients and identifies a potential therapeutic target and biomarker., info:eu-repo/semantics/published

Genetic redundancy has evolved as a way for human cells to survive the loss of genes that are single copy and essential in other organisms, but also allows tumours to survive despite having highly rearranged genomes. In this study we CRISPR screen 1191 gene pairs, including paralogues and known and predicted synthetic lethal interactions to identify 105 gene combinations whose co-disruption results in a loss of cellular fitness. 27 pairs influence fitness across multiple cell lines including the paralogues FAM50A/FAM50B, two genes of unknown function. Silencing of FAM50B occurs across a range of tumour types and in this context disruption of FAM50A reduces cellular fitness whilst promoting micronucleus formation and extensive perturbation of transcriptional programmes. Our studies reveal the fitness effects of FAM50A/FAM50B in cancer cells., Systematic screens for synthetic lethal gene pairs represent a powerful approach to define interactions that may be exploited in the clinic. Here, the authors use a dual-guide CRISPR screening approach to screen a curated selection of gene pairs across three cell lines and validate the Fam50a/Fam50b synthetic lethality was in isogenic cell models as well as in an in vivo setting.

Natural killer (NK) cells protect against intracellular infection and cancer. These properties are exploited in oncolytic virus (OV) therapy, where anti-viral responses enhance anti-tumour immunity. We have analysed the mechanism by which reovirus, an oncolytic dsRNA virus, modulates human NK cell activity. Reovirus activates NK cells in a type I interferon (IFN-I) dependent manner, resulting in STAT1 and STAT4 signalling in both CD56dimand CD56brightNK cell subsets. Gene expression profiling revealed the dominance of IFN-I responses and identified induction of genes associated with NK cell cytotoxicity and cell cycle progression, with distinct responses in the CD56dimand CD56brightsubsets. However, reovirus treatment, acting via IFN-I, inhibited NK cell proliferative responses to IL-15 and was associated with reduced AKT signalling.In vivo, human CD56dimand CD56brightNK cells responded with similar kinetics to reovirus treatment, but CD56brightNK cells were transiently lost from the peripheral circulation at the peak of the IFN-I response, suggestive of their redistribution to secondary lymphoid tissue. These results show that reovirus modulates a spectrum of NK cell activityin vivo, encompassing direct action on tumour cells and the regulation of adaptive immunity. Such activity is likely to mirror NK cell responses to natural viral infection.

Targeting the MAPK signaling pathway has transformed the treatment of metastatic melanoma. CRISPR‐Cas9 genetic screens provide a genome‐wide approach to uncover novel genetic dependencies that might serve as therapeutic targets. Here, we analyzed recently reported CRISPR‐Cas9 screens comparing data from 28 melanoma cell lines and 313 cell lines of other tumor types in order to identify fitness genes related to melanoma. We found an average of 1,494 fitness genes in each melanoma cell line. We identified 33 genes, inactivation of which specifically reduced the fitness of melanoma. This set of tumor type‐specific genes includes established melanoma fitness genes as well as many genes that have not previously been associated with melanoma growth. Several genes encode proteins that can be targeted using available inhibitors. We verified that genetic inactivation of DUSP4 and PPP2R2A reduces the proliferation of melanoma cells. DUSP4 encodes an inhibitor of ERK, suggesting that further activation of MAPK signaling activity through its loss is selectively deleterious to melanoma cells. Collectively, these data present a resource of genetic dependencies in melanoma that may be explored as potential therapeutic targets.

Epidemiological studies have identified many environmental agents that appear to significantly increase cancer risk in human populations. By analysis of tumour genomes from mice chronically exposed to one of 20 known or suspected human carcinogens, we reveal that most agents do not generate distinct mutational signatures or increase mutation burden, with most mutations, including driver mutations, resulting from tissue specific endogenous processes. We identify signatures resulting from exposure to cobalt and vinylidene chloride and link distinct human signatures (SBS19 and SBS42) with 1,2,3-Trichloropropane (TCP), a haloalkane and pollutant of drinking water, and find these and other signatures in human tumour genomes. We define the cross-species genomic landscape of tumours induced by an important compendium of agents with relevance to human health.

BACKGROUND: Human prostate cancers display numerous DNA methylation changes compared to normal tissue samples. However, definitive identification of features related to the cells' malignant status has been compromised by the predominance of cells with luminal features in prostate cancers. METHODS: We generated genome-wide DNA methylation profiles of cell subpopulations with basal or luminal features isolated from matched prostate cancer and normal tissue samples. RESULTS: Many frequent DNA methylation changes previously attributed to prostate cancers are here identified as differences between luminal and basal cells in both normal and cancer samples. We also identified changes unique to each of the two cancer subpopulations. Those specific to cancer luminal cells were associated with regulation of metabolic processes, cell proliferation and epithelial development. Within the prostate cancer TCGA dataset, these changes were able to distinguish not only cancers from normal samples, but also organ-confined cancers from those with extraprostatic extensions. Using changes present in both basal and luminal cancer cells, we derived a new 17-CpG prostate cancer signature with high predictive power in the TCGA dataset. CONCLUSIONS: This study demonstrates the importance of comparing phenotypically matched prostate cell populations from normal and cancer tissues to unmask biologically and clinically relevant DNA methylation changes.

Breast cancer can occur in either gender; however, it is rare in men, accounting for

Simple Summary High-grade gliomas remain difficult-to-treat cancers. Novel treatment options include targeting glycogen synthase kinase 3 (GSK-3) to induce cell death but the mode of drug activity remains unknown and combination with conventional treatment including radiotherapy has not been explored. Here, we describe the effect of targeting GSK-3 with the inhibitor AZD2858 in in vitro and in vivo models of glioma. We established that AZD2858 exposure induces mitotic defects leading to cell death in patient-derived glioma cell lines and tumor growth delay in glioma xenografts. Co-administration also enhanced the effect of radiotherapy. We therefore propose AZD2858 as an adjuvant to radiotherapy in high-grade glioma. Abstract Background: Previous data on glycogen synthase kinase 3 (GSK-3) inhibition in cancer models support a cytotoxic effect with selectivity for tumor cells compared to normal tissue but the effect of these inhibitors in glioma has not been widely studied. Here, we investigate their potential as cytotoxics in glioma. Methods: We assessed the effect of pharmacologic GSK-3 inhibition on established (U87, U251) and patient-derived (GBM1, GBM4) glioblastoma (GBM) cell lines using cytotoxicity assays as well as undertaking a detailed investigation of the effect on cell cycle, mitosis, and centrosome biology. We also assessed drug uptake and efficacy of GSK-3 inhibition alone and in combination with radiation in xenograft models. Results: Using the selective GSK-3 inhibitor AZD2858, we demonstrated single agent cytotoxicity in two patient-derived glioma cell lines (GBM1, GBM4) and two established cell lines (U251 and U87) with IC50 in the low micromolar range promoting centrosome disruption, failed mitosis, and S-phase arrest. Glioma xenografts exposed to AZD2858 also showed growth delay compared to untreated controls. Combined treatment with radiation increased the cytotoxic effect of clinical radiation doses in vitro and in orthotopic glioma xenografts. Conclusions: These data suggest that GSK-3 inhibition promotes cell death in glioma through disrupting centrosome function and promoting mitotic failure and that AZD2858 is an effective adjuvant to radiation at clinical doses.

Many traits of cancer progression (e.g., development of metastases or resistance to therapy) are facilitated by tumour evolution: Darwinian selection of subclones with distinct genotypes or phenotypes that enable such progression. Characterising these subclones provide an opportunity to develop drugs to better target their specific properties but requires the accurate identification of somatic mutations shared across multiple spatiotemporal tumours from the same patient. Current best practices for calling somatic mutations are optimised for single samples, and risk being too conservative to identify shared mutations with low prevalence in some samples. We reasoned that datasets from multiple matched tumours can be used for mutual validation and thus propose an adapted two-stage approach: (1) low-stringency mutation calling to identify mutations shared across samples irrespective of the weight of evidence in a single sample; (2) high-stringency mutation calling to further characterise mutations present in a single sample. We applied our approach to three-independent cohorts of paired primary and recurrent glioblastoma tumours, two of which have previously been analysed using existing approaches, and found that it significantly increased the amount of biologically relevant shared somatic mutations identified. We also found that duplicate removal was detrimental when identifying shared somatic mutations. Our approach is also applicable when multiple datasets e.g. DNA and RNA are available for the same tumour.

Recent findings from our group, and the wider community, show that standard treatment does not impose an apparent bottleneck on the clonal evolution of adult glioblastoma (GBM), implying a lack of direct therapeutic opportunity. This does not negate the possibility that multiple treatment-resistance mechanisms co-exist in tumours, repeated across patients, making a combination of targeted therapies a potentially effective approach. We investigated whether treatment resistance may be driven by selection of cellular properties conferred above the level of the genome. Differential expression analysis was performed on 23 pairs of primary and recurrent tumours from patients who received standard treatment and had a local recurrence treated by surgery and second line chemotherapy. This revealed a treatment-induced shift in cell states linked to normal neurodevelopment. The latter is orchestrated by cascades of transcription factors. We, therefore, applied a bespoke gene set enrichment analysis to our paired expression data to investigate whether any factors were implicated in co-regulation of the genes that were altered through therapy. This identified a specific chromatin remodelling machinery, instrumental in normal neurogenesis. We validated our results in an independent cohort of 22 paired GBM samples. Our results suggest that the chromatin remodelling machinery is responsible for determining transcriptional hierarchies in GBM, shown elsewhere to have different treatment sensitivities such that their relative abundances are altered through treatment.

BackgroundMicroRNAs play an important role in the regulation of mRNA translation, and have therapeutic potential in cancer and other diseases.MethodsTo profile the landscape of microRNAs with significant cytotoxicity in the context of glioblastoma (GBM), we performed a high-throughput screen using a synthetic oligonucleotide library representing all known human microRNAs in adult and pediatric GBM cells. Bio-informatics analysis were used to refine this list and the top seven microRNAs were validated in a larger panel of cells by flow-cytometry, and RTqPCR. The downstream mechanism of the strongest and most consistent candidate was investigated by siRNAs, 3’UTR luciferase assays and Western Blotting.ResultsOur screen identified ∼100 significantly cytotoxic microRNAs with 70% concordance between cell lines. MicroRNA-1300 (miR-1300) was the most potent and robust candidate. We observed a striking binucleated phenotype in miR-1300 expressing cells and characterized the mechanism of action as cytokinesis failure followed by apoptosis, which was observed in an extended GBM cell panel including two stem-like patient-derived cultures. We identified the physiological role of miR-1300 as a regulator of endomitosis in megakaryocyte differentiation where blockade of cytokinesis is an essential step. In glioblastoma cells, the oncogene Epithelial Cell Transforming 2 (ECT2) was validated as a direct key target of miR-1300. ECT2 siRNA phenocopied the effects of miR-1300, and its overexpression led to a significant rescue of miR-1300 induced binucleation.ConclusionMiR-1300 was identified as a novel regulator of endomitosis with translatable potential for therapeutic application. The datatasets will be a resource for the neuro-oncology community.Key points (2 or 3 key points 85 characters plus spaces each)70% of cytotoxic microRNAs were shared between adult and pediatric glioblastoma cellsMiR-1300 expression is restricted to endomitosis within megakaryocyte differentiationMiR-1300’s ectopic expression is a potent and promising therapeutic tool in cancerImportance of StudyPrevious functional studies of microRNAs involved in the regulation of glioblastoma cell proliferation and/or survival have focused on adult glioblastoma alone and are restricted to only a few microRNAs at a time. Our study provides the first encompassing landscape of potent cytotoxic microRNAs in pediatric and adult glioblastoma.Not only, does our data provide an invaluable resource for the research community but it also revealed that 70% of microRNAs with significant cytotoxicity were shared by adult and pediatric cells. Finally, we identified and characterized the previously undescribed role of microRNA-1300 in the tight regulation of megakaryocyte differentiation into platelets and how, when expressed outside of this context, miR-1300 consistently causes cytokinesis failure followed by apoptosis, and thus represents a powerful cytotoxic tool with potential for translation towards therapeutic applications.

BACKGROUND Glioblastoma (GBM) tumours recur following standard treatment in almost all cases. We use ‘omics technologies to simultaneously profile pairs of primary and matched recurrent GBM to specifically identify and characterise the cells that resisted treatment, with the aim of determining how to more effectively kill them. MATERIAL AND METHODS We have analysed high coverage RNAseq data from pairs of GBM tumours: primary de novo tumour and matched local recurrence from patients that underwent standard therapy. Our original cohort constituted 23 pairs and our validation cohort was an additional 22 pairs. We also cultured two plates of spheroids directly from a patient’s GBM, treating one with radiation and temozolomide. We monitored growth and captured and sequenced RNA from single cells at two time-points: one week post-treatment when the deviation between untreated and treated spheroid growth curves was most pronounced; and three weeks post-treatment when the growth rate of treated spheroids had recovered. We investigated differential gene expression between primary and recurrent pairs, and single cells pre- and post-treatment, and performed a bespoke per patient gene set enrichment analysis. RESULTS Differential gene expression analysis in 23 tumour pairs indicated a treatment-induced shift in cell states linked to normal neurogenesis and prompted us to develop a novel gene set enrichment analysis approach to identify gene regulatory factors that may orchestrate such a shift. This revealed the significant and universal dysregulation of genes, through therapy, that are targeted by a specific chromatin remodeling machinery. This finding was validated in an independent cohort of 22 further GBM pairs. To understand the therapeutic potential of this finding we must determine whether genes are dysregulated through therapy owing to a) their fixed expression in inherently treatment resistance cells in the primary tumour which get selected during therapy to increase the signal of that profile, or b) changes in expression during the process of cells acquiring treatment resistance. To inspect this, we analysed single cell gene expression data from GBM spheroids pre- and post-treatment. We found that there was significant dysregulation of the genes associated with the chromatin remodeling complex but only at the three-week post-treatment time-point. CONCLUSION Our results indicate that GBM cells are being transcriptionally reprogrammed in response to treatment; the mechanism of which may represent a therapeutic opportunity.

BackgroundGlioblastoma (GBM) is a fatal and incurable brain cancer with a dismal prognosis. In order to impact on this disease, we need to understand how infiltrating, non resectable tumour cells resist chemoradiation and facilitate disease recurrence. To this end, we generated or acquired bulk tumour RNA sequencing data from 45 paired primary and locally recurrent GBM tumours (split into original and validation cohorts) from patients that received standard treatment. We also generated DNA methylation profiles for 9 pairs and sequenced RNA from single cells isolated from a patient derived GBM spheroid model at different timepoints following in vitro chemoradiation.ResultsWe have identified a set of genes with Jumonji and AT-Rich Interacting Domain 2 (JARID2) binding sites in their promoters that are universally dysregulated in post-standard treatment recurrent GBMs compared to the primary tumour. The direction of dysregulation is patient-dependent and not associated with differential promoter DNA methylation. Our in vitro experiments suggest that this dysregulation occurs dynamically following treatment as opposed to resulting from selection of cells with specific expression profiles.ConclusionJARID2 is an accessory protein to a chromatin remodeling complex, responsible for histone modifications observed during cell state transitions in both normal brain and GBM. We propose that JARID2 facilitates GBM recurrence following treatment by indirect transcriptional reprogramming of surviving cells in whichever manner is needed to reproduce the phenotypic heterogeneity required for tumour regrowth in vivo. The mechanism of this reprogramming may present a therapeutic vulnerability for more effective treatment of GBM.

The immune response to melanoma improves the survival in untreated patients and predicts the response to immune checkpoint blockade. Here, we report genetic and environmental predictors of the immune response in a large primary cutaneous melanoma cohort. Bioinformatic analysis of 703 tumor transcriptomes was used to infer immune cell infiltration and to categorize tumors into immune subgroups, which were then investigated for association with biological pathways, clinicopathologic factors, and copy number alterations. Three subgroups, with “low”, “intermediate”, and “high” immune signals, were identified in primary tumors and replicated in metastatic tumors. Genes in the low subgroup were enriched for cell-cycle and metabolic pathways, whereas genes in the high subgroup were enriched for IFN and NF-κB signaling. We identified high MYC expression partially driven by amplification, HLA-B downregulation, and deletion of IFNγ and NF-κB pathway genes as the regulators of immune suppression. Furthermore, we showed that cigarette smoking, a globally detrimental environmental factor, modulates immunity, reducing the survival primarily in patients with a strong immune response. Together, these analyses identify a set of factors that can be easily assessed that may serve as predictors of response to immunotherapy in patients with melanoma. Significance: These findings identify novel genetic and environmental modulators of the immune response against primary cutaneous melanoma and predict their impact on patient survival. See related commentary by Anichini, p. 2457

Tumor stem cells and malignant multicellular networks have been separately implicated in the therapeutic resistance of Glioblastoma Multiforme (GBM), the most aggressive type of brain cancer in adults. We show that small molecule inhibition of RHO-associated serine/threonine kinase (ROCKi) significantly promoted the outgrowth of neurite-like cell projections in cultures of heterogeneous patient-derived GBM stem-like cells. These projections formed de novo-induced cellular network (iNet) ‘webs’, which regressed after withdrawal of ROCKi. Connected cells within the iNet web exhibited long range calcium signal transmission, and significant lysosomal and mitochondrial trafficking. In contrast to their less-connected vehicle control counterparts, iNet cells remained viable and proliferative after high-dose radiation. These findings demonstrate a link between ROCKi-regulated cell projection dynamics and the formation of radiation-resistant multicellular networks. Our study identifies means to reversibly induce iNet webs ex vivo, and may thereby accelerate future studies into the biology of GBM cellular networks.

SummaryTumor stem cells and malignant multicellular networks have been separately implicated in the therapeutic resistance of Glioblastoma Multiforme (GBM), the most aggressive type of brain cancer in adults. We show that small molecule inhibition of RHO-associated serine/threonine kinase (ROCKi) significantly promoted the outgrowth of neurite-like cell projections in cultures of heterogeneous patient-derived GBM stem-like cells. These projections formedde novo-induced cellular network (iNet) ‘webs’, which regressed after withdrawal of ROCKi. Connected cells within the iNet web exhibited long range calcium signal transmission, and significant lysosomal and mitochondrial trafficking. In contrast to their less-connected vehicle control counterparts, iNet cells remained viable and proliferative after high-dose radiation. These findings demonstrate a link between ROCKi-regulated cell projection dynamics and the formation of radiation-resistant multicellular networks. Our study identifies means to reversibly induce iNet websex vivo, and may thereby accelerate future studies into the biology of GBM cellular networks.

Pharmacological inhibition of uncontrolled cell growth with small-molecule inhibitors is a potential strategy for treating glioblastoma multiforme (GBM), the most malignant primary brain cancer. We showed that the synthetic small-molecule KHS101 promoted tumor cell death in diverse GBM cell models, independent of their tumor subtype, and without affecting the viability of noncancerous brain cell lines. KHS101 exerted cytotoxic effects by disrupting the mitochondrial chaperone heat shock protein family D member 1 (HSPD1). In GBM cells, KHS101 promoted aggregation of proteins regulating mitochondrial integrity and energy metabolism. Mitochondrial bioenergetic capacity and glycolytic activity were selectively impaired in KHS101-treated GBM cells. In two intracranial patient-derived xenograft tumor models in mice, systemic administration of KHS101 reduced tumor growth and increased survival without discernible side effects. These findings suggest that targeting of HSPD1-dependent metabolic pathways might be an effective strategy for treating GBM.

Summary Tumour evolution results in progressive cancer phenotypes such as metastatic spread and treatment resistance. To better treat cancers, we must characterize tumour evolution and the genetic events that confer progressive phenotypes. This is facilitated by high coverage genome or exome sequencing. However, the best approach by which, or indeed whether, these data can be used to accurately model and interpret underlying evolutionary dynamics is yet to be confirmed. Establishing this requires sequencing data from appropriately heterogeneous tumours in which the exact trajectory and combination of events occurring throughout its evolution are known. We therefore developed HeteroGenesis: a tool to generate realistically evolved tumour genomes, which can be sequenced using weighted-Wessim (w-Wessim), an in silico exome sequencing tool that we have adapted from previous methods. HeteroGenesis simulates more complex and realistic heterogeneous tumour genomes than existing methods, can model different evolutionary dynamics, and enables the creation of multi-region and longitudinal data. Availability and implementation HeteroGenesis and w-Wessim are freely available under the GNU General Public Licence from https://github.com/GeorgetteTanner, implemented in Python and supported on linux and MS Windows. Supplementary information Supplementary data are available at Bioinformatics online.

Mutations in pre-mRNA processing factors (PRPFs) cause 40% of autosomal dominant retinitis pigmentosa (RP), but it is unclear why mutations in ubiquitously expressed PRPFs cause retinal disease. To understand the molecular basis of this phenotype, we have generated RP type 11 (PRPF31-mutated) patient-specific retinal organoids and retinal pigment epithelium (RPE) from induced pluripotent stem cells (iPSC). Impaired alternative splicing of genes encoding pre-mRNA splicing proteins occurred in patient-specific retinal cells and Prpf31 /- mouse retinae, but not fibroblasts and iPSCs, providing mechanistic insights into retinal-specific phenotypes of PRPFs. RPE was the most affected, characterised by loss of apical-basal polarity, reduced trans-epithelial resistance, phagocytic capacity, microvilli, and cilia length and incidence. Disrupted cilia morphology was observed in patient-derived-photoreceptors that displayed progressive features associated with degeneration and cell stress. In situ gene-editing of a pathogenic mutation rescued key structural and functional phenotypes in RPE and photoreceptors, providing proof-of-concept for future therapeutic strategies.

SummaryMutations in pre-mRNA processing factors (PRPFs) cause 40% of autosomal dominant retinitis pigmentosa (RP), but it is unclear why mutations in ubiquitously expressed PRPFs cause retinal disease. To understand the molecular basis of this phenotype, we have generated RP type 11 (PRPF31-mutated) patient-specific retinal organoids and retinal pigment epithelium (RPE) from induced pluripotent stem cells (iPSC). Impaired alternative splicing of genes encoding pre-mRNA splicing proteins occurred in patient-specific retinal cells and Prpf31+/− mouse retinae, but not fibroblasts and iPSCs, providing mechanistic insights into retinal-specific phenotypes of PRPFs. RPE was the most affected, characterised by loss of apical-basal polarity, reduced trans-epithelial resistance, phagocytic capacity, microvilli, and cilia length and incidence. Disrupted cilia morphology was observed in patient-derived-photoreceptors that displayed progressive features associated with degeneration and cell stress. In situ gene-editing of a pathogenic mutation rescued key structural and functional phenotypes in RPE and photoreceptors, providing proof-of-concept for future therapeutic strategies.eTOCPRPF31 is a ubiquitously expressed pre-mRNA processing factor that when mutated causes autosomal dominant RP. Using a patient-specific iPSC approach, Buskin and Zhu et al. show that retinal-specific defects result from altered splicing of genes involved in the splicing process itself, leading to impaired splicing, loss of RPE polarity and diminished phagocytic ability as well as reduced cilia incidence and length in both photoreceptors and RPE.HighlightsSuccessful generation of iPSC-derived RPE and photoreceptors from four RP type 11 patientsRPE cells express the mutant PRPF31 protein and show the lowest expression of wildtype proteinPRPF31 mutations result in altered splicing of genes involved in pre-mRNA splicing in RPE and retinal organoidsPrpf31 haploinsufficiency results in altered splicing of genes involved in pre-mRNA splicing in mouse retinaRPE cells display loss of polarity, reduced barrier function and phagocytosisPhotoreceptors display shorter and fewer cilia and degenerative featuresRPE cells display most abnormalities suggesting they might be the primary site of pathogenesisIn situ gene editing corrects the mutation and rescues key phenotypes

Bladder cancer incurs a higher lifetime treatment cost than other cancers due to frequent recurrence of non-invasive disease. Improved prognostic biomarkers and localised therapy are needed for this large patient group. We defined two major genomic subtypes of primary stage Ta tumors that showed differential risk of recurrence. The higher risk subtype was characterised by loss of 9q including TSC1, increased KI67 labelling index, upregulated glycolysis, DNA repair, mTORC1 signaling, features of the unfolded protein response and altered cholesterol homeostasis. Comparison with muscle-invasive bladder cancer mutation profiles revealed lower overall mutation rates and more frequent mutations in RHOB and chromatin modifier genes. More mutations in the histone-lysine demethylase KDM6A were present in non-invasive tumors from females than males.

Pharmacological inhibition of uncontrolled cell growth with small molecule inhibitors is a potential strategy against glioblastoma multiforme (GBM), the most malignant primary brain cancer. Phenotypic profiling of the neurogenic small molecule KHS101 revealed the chemical induction of lethal cellular degradation in molecularly-diverse GBM cells, independent of their tumor subtype, whereas non-cancerous brain cells remained viable. Mechanism-of-action (MOA) studies showed that KHS101 specifically bound and inhibited the mitochondrial chaperone HSPD1. In GBM but not non-cancerous brain cells, KHS101 elicited the aggregation of an enzymatic network that regulates energy metabolism. Compromised glycolysis and oxidative phosphorylation (OXPHOS) resulted in the metabolic energy depletion in KHS101-treated GBM cells. Consistently, KHS101 induced key mitochondrial unfolded protein response factor DDIT3 in vitro and in vivo, and significantly reduced intracranial GBM xenograft tumor growth upon systemic administration, without discernible side effects. These findings suggest targeting of HSPD1-dependent oncometabolic pathways as an anti-GBM therapy.

Many traits of cancer progression (e.g., development of metastases or resistance to therapy) are facilitated by tumour evolution: Darwinian selection of subclones with distinct genotypes or phenotypes that enable such progression. Characterising these subclones provide an opportunity to develop drugs to better target their specific properties but requires the accurate identification of somatic mutations shared across multiple spatiotemporal tumours from the same patient. Current best practices for calling somatic mutations are optimised for single samples, and risk being too conservative to identify shared mutations with low prevalence in some samples. We reasoned that datasets from multiple matched tumours can be used for mutual validation and thus propose an adapted two-stage approach: (1) low-stringency mutation calling to identify mutations shared across samples irrespective of the weight of evidence in a single sample; (2) high-stringency mutation calling to further characterise mutations present in a single sample. We applied our approach to three-independent cohorts of paired primary and recurrent glioblastoma tumours, two of which have previously been analysed using existing approaches, and found that it significantly increased the amount of biologically relevant shared somatic mutations identified. We also found that duplicate removal was detrimental when identifying shared somatic mutations. Our approach is also applicable when multiple datasets e.g. DNA and RNA are available for the same tumour.

A collaboration between cancer biologists and academic software engineers has been exploring the development of an agent-based simulator to inform and support work on the dynamics of cell proliferation in the study of prostate disorders. The research has influenced and been informed by the CoSMoS project. This paper presents the simulation project (which is not yet complete). We reflect on the reality of following CoSMoS principles; we describe the domain exploration and show how software modelling approaches (here, Petri nets, state diagrams) can be used to express both biological and software models. We explore fitness for purpose and consider ways to present a fitness argument. We consider issues in choosing simulation media and mapping from domain models through to code. The implementation emphasis is on traceability to support reuse and extension of the simulator, as well as demonstrable fitness for purpose. Initial work on calibration is presented. We discuss the calibration results, that both support and challenge the design and assumptions captured in the domain modelling and development activities.