Disruption of metazoan gene regulatory networks in cancer alters the balance of co-expression between genes of unicellular and multicellular origins.

BACKGROUND: Metazoans inherited genes from unicellular ancestors that perform essential biological processes such as cell division, metabolism, and protein translation. Multicellularity requires careful control and coordination of these unicellular genes to maintain tissue integrity and homeostasis. Gene regulatory networks (GRNs) that arose during metazoan evolution are frequently altered in cancer, resulting in over-expression of unicellular genes. We propose that an imbalance in co-expression of unicellular (UC) and multicellular (MC) genes is a driving force in cancer. RESULTS: We combine gene co-expression analysis to infer changes to GRNs in cancer with protein sequence conservation data to distinguish genes with UC and MC origins. Co-expression networks created using RNA sequencing data from 31 tumor types and normal tissue samples are divided into modules enriched for UC genes, MC genes, or mixed UC-MC modules. The greatest differences between tumor and normal tissue co-expression networks occur within mixed UC-MC modules. MC and UC genes not commonly co-expressed in normal tissues form distinct co-expression modules seen only in tumors. The degree of rewiring of genes within mixed UC-MC modules increases with tumor grade and stage. Mixed UC-MC modules are enriched for somatic mutations in cancer genes, particularly amplifications, suggesting an important driver of the rewiring observed in tumors is copy number changes. CONCLUSIONS: Our study shows the greatest changes to gene co-expression patterns during tumor progression occur between genes of MC and UC origins, implicating the breakdown and rewiring of metazoan gene regulatory networks in cancer development and progression.

Copy number architectures define treatment-mediated selection of lethal prostate cancer clones.

Despite initial responses to hormone treatment, metastatic prostate cancer invariably evolves to a lethal state. To characterize the intra-patient evolutionary relationships of metastases that evade treatment, we perform genome-wide copy number profiling and bespoke approaches targeting the androgen receptor (AR) on 167 metastatic regions from 11 organs harvested post-mortem from 10 men who died from prostate cancer. We identify diverse and patient-unique alterations clustering around the AR in metastases from every patient with evidence of independent acquisition of related genomic changes within an individual and, in some patients, the co-existence of AR-neutral clones. Using the genomic boundaries of pan-autosome copy number changes, we confirm a common clone of origin across metastases and diagnostic biopsies, and identified in individual patients, clusters of metastases occupied by dominant clones with diverged autosomal copy number alterations. These autosome-defined clusters are characterized by cluster-specific AR gene architectures, and in two index cases are topologically more congruent than by chance (p-values 3.07 × 10-8 and 6.4 × 10-4). Integration with anatomical sites suggests patterns of spread and points of genomic divergence. Here, we show that copy number boundaries identify treatment-selected clones with putatively distinct lethal trajectories.

Comprehensive profiling identifies tumour and immune microenvironmental differences in clinical subsets of cutaneous squamous cell carcinoma.

BACKGROUND: Cutaneous squamous cell carcinomas (cSCCs) are the second most diagnosed skin cancer worldwide; however, little is known about the pathobiological factors that contribute to the diverse clinical outcomes seen. OBJECTIVES: To profile cSCCs comprehensively and identify the pathological processes that contribute to the disparities seen in their clinical behaviour. METHODS: We characterized the genomic, transcriptomic and immunohistochemical profiles of 211 cSCC tumours, including 37 cSCCs from immunocompromised patients. RESULTS: cSCCs from immunocompromised patients were characterized by a lack of B cells in the peritumoral stroma compared with immunocompetent patients. Further, an abundance of a memory B-cell-like population in the peritumoral stroma was associated with a better prognosis in all patients (immunocompetent and immunocompromised), as well as only immunocompetent patients. No differences in genetic -variants, tumour mutational burden or mutational signatures were observed between cSCCs from immunocompetent and immunocompromised patients. Thus, differences in survival between cSCCs from immunocompromised patients and immunocompetent patients are not likely to be driven by tumour genomic factors, but may be associated with differential host immune response. cSCC not from a primary head and neck site had lower tumour mutational burden and exhibited upregulation of the epithelial-mesenchymal transition programme compared with head and neck cSCC. Both factors were implicated with poorer responses to immune checkpoint inhibition, and the latter with poorer survival. CONCLUSIONS: We identified tumour and host immune factors that contribute to the disparate clinical behaviour of cSCC, with broad translational application, including prognostication, treatment prediction to current therapies and the identification of novel anticancer therapy approaches in cSCC.

Spatial analysis with SPIAT and spaSim to characterize and simulate tissue microenvironments.

Spatial proteomics technologies have revealed an underappreciated link between the location of cells in tissue microenvironments and the underlying biology and clinical features, but there is significant lag in the development of downstream analysis methods and benchmarking tools. Here we present SPIAT (spatial image analysis of tissues), a spatial-platform agnostic toolkit with a suite of spatial analysis algorithms, and spaSim (spatial simulator), a simulator of tissue spatial data. SPIAT includes multiple colocalization, neighborhood and spatial heterogeneity metrics to characterize the spatial patterns of cells. Ten spatial metrics of SPIAT are benchmarked using simulated data generated with spaSim. We show how SPIAT can uncover cancer immune subtypes correlated with prognosis in cancer and characterize cell dysfunction in diabetes. Our results suggest SPIAT and spaSim as useful tools for quantifying spatial patterns, identifying and validating correlates of clinical outcomes and supporting method development.

Inhibition of mutant IDH1 promotes cycling of acute myeloid leukemia stem cells.

Approximately 20% of acute myeloid leukemia (AML) patients carry mutations in IDH1 or IDH2 that result in over-production of the oncometabolite D-2-hydroxyglutarate (2-HG). Small molecule inhibitors that block 2-HG synthesis can induce complete morphological remission; however, almost all patients eventually acquire drug resistance and relapse. Using a multi-allelic mouse model of IDH1-mutant AML, we demonstrate that the clinical IDH1 inhibitor AG-120 (ivosidenib) exerts cell-type-dependent effects on leukemic cells, promoting delayed disease regression. Although single-agent AG-120 treatment does not fully eradicate the disease, it increases cycling of rare leukemia stem cells and triggers transcriptional upregulation of the pyrimidine salvage pathway. Accordingly, AG-120 sensitizes IDH1-mutant AML to azacitidine, with the combination of AG-120 and azacitidine showing vastly improved efficacy in vivo. Our data highlight the impact of non-genetic heterogeneity on treatment response and provide a mechanistic rationale for the observed combinatorial effect of AG-120 and azacitidine in patients.

Cystathionine-ß-synthase is essential for AKT-induced senescence and suppresses the development of gastric cancers with PI3K/AKT activation.

Hyperactivation of oncogenic pathways downstream of RAS and PI3K/AKT in normal cells induces a senescence-like phenotype that acts as a tumor-suppressive mechanism that must be overcome during transformation. We previously demonstrated that AKT-induced senescence (AIS) is associated with profound transcriptional and metabolic changes. Here, we demonstrate that human fibroblasts undergoing AIS display upregulated cystathionine-ß-synthase (CBS) expression and enhanced uptake of exogenous cysteine, which lead to increased hydrogen sulfide (H2S) and glutathione (GSH) production, consequently protecting senescent cells from oxidative stress-induced cell death. CBS depletion allows AIS cells to escape senescence and re-enter the cell cycle, indicating the importance of CBS activity in maintaining AIS. Mechanistically, we show this restoration of proliferation is mediated through suppressing mitochondrial respiration and reactive oxygen species (ROS) production by reducing mitochondrial localized CBS while retaining antioxidant capacity of transsulfuration pathway. These findings implicate a potential tumor-suppressive role for CBS in cells with aberrant PI3K/AKT pathway activation. Consistent with this concept, in human gastric cancer cells with activated PI3K/AKT signaling, we demonstrate that CBS expression is suppressed due to promoter hypermethylation. CBS loss cooperates with activated PI3K/AKT signaling in promoting anchorage-independent growth of gastric epithelial cells, while CBS restoration suppresses the growth of gastric tumors in vivo. Taken together, we find that CBS is a novel regulator of AIS and a potential tumor suppressor in PI3K/AKT-driven gastric cancers, providing a new exploitable metabolic vulnerability in these cancers.

Tumor immune microenvironment of primary prostate cancer with and without germline mutations in homologous recombination repair genes.

BACKGROUND: Aberrations in homologous recombination repair (HRR) genes are emerging as important biomarkers for personalized treatment in prostate cancer (PCa). HRR deficiency (HRD) could affect the tumor immune microenvironment (TIME), potentially contributing to differential responses to poly ADP-ribose polymerase (PARP) inhibitors and immune checkpoint inhibitors. Spatial distribution of immune cells in a range of cancers identifies novel disease subtypes and is related to prognosis. In this study we aimed to determine the differences in the TIME of PCa with and without germline (g) HRR mutations. METHODS: We performed gene expression analysis, multiplex immunohistochemistry of T and B cells and quantitative spatial analysis of PCa samples from 36 patients with gHRD and 26 patients with sporadic PCa. Samples were archival tumor tissue from radical prostatectomies with the exception of one biopsy. Results were validated in several independent cohorts. RESULTS: Although the composition of the T cell and B cells was similar in the tumor areas of gHRD-mutated and sporadic tumors, the spatial profiles differed between these cohorts. We describe two T-cell spatial profiles across primary PCa, a clustered immune spatial (CIS) profile characterized by dense clusters of CD4+ T cells closely interacting with PD-L1+ cells, and a free immune spatial (FIS) profile of CD8+ cells in close proximity to tumor cells. gHRD tumors had a more T-cell inflamed microenvironment than sporadic tumors. The CIS profile was mainly observed in sporadic tumors, whereas a FIS profile was enriched in gHRD tumors. A FIS profile was associated with lower Gleason scores, smaller tumors and longer time to biochemical recurrence and metastasis. CONCLUSIONS: gHRD-mutated tumors have a distinct immune microenvironment compared with sporadic tumors. Spatial profiling of T-cells provides additional information beyond T-cell density and is associated with time to biochemical recurrence, time to metastasis, tumor size and Gleason scores.

Adaptive translational reprogramming of metabolism limits the response to targeted therapy in BRAFV600 melanoma.

Despite the success of therapies targeting oncogenes in cancer, clinical outcomes are limited by residual disease that ultimately results in relapse. This residual disease is often characterized by non-genetic adaptive resistance, that in melanoma is characterised by altered metabolism. Here, we examine how targeted therapy reprograms metabolism in BRAF-mutant melanoma cells using a genome-wide RNA interference (RNAi) screen and global gene expression profiling. Using this systematic approach we demonstrate post-transcriptional regulation of metabolism following BRAF inhibition, involving selective mRNA transport and translation. As proof of concept we demonstrate the RNA processing kinase U2AF homology motif kinase 1 (UHMK1) associates with mRNAs encoding metabolism proteins and selectively controls their transport and translation during adaptation to BRAF-targeted therapy. UHMK1 inactivation induces cell death by disrupting therapy induced metabolic reprogramming, and importantly, delays resistance to BRAF and MEK combination therapy in multiple in vivo models. We propose selective mRNA processing and translation by UHMK1 constitutes a mechanism of non-genetic resistance to targeted therapy in melanoma by controlling metabolic plasticity induced by therapy.

Reprogrammed mRNA translation drives resistance to therapeutic targeting of ribosome biogenesis.

Elevated ribosome biogenesis in oncogene-driven cancers is commonly targeted by DNA-damaging cytotoxic drugs. Our previous first-in-human trial of CX-5461, a novel, less genotoxic agent that specifically inhibits ribosome biogenesis via suppression of RNA polymerase I (Pol I) transcription, revealed single-agent efficacy in refractory blood cancers. Despite this clinical response, patients were not cured. In parallel, we demonstrated a marked improvement in the in vivo efficacy of CX-5461 in combination with PI3K/AKT/mTORC1 pathway inhibitors. Here, we reveal the molecular basis for this improved efficacy observed in vivo, which is associated with specific suppression of translation of mRNAs encoding regulators of cellular metabolism. Importantly, acquired resistance to this cotreatment is driven by translational rewiring that results in dysregulated cellular metabolism and induction of a cAMP-dependent pathway critical for the survival of blood cancers including lymphoma and acute myeloid leukemia. Our studies thus identify key molecular mechanisms underpinning the response of blood cancers to selective inhibition of ribosome biogenesis and define metabolic vulnerabilities that will facilitate the rational design of more effective regimens for Pol I-directed therapies.

CX-5461 activates the DNA damage response and demonstrates therapeutic efficacy in high-grade serous ovarian cancer.

Acquired resistance to PARP inhibitors (PARPi) is a major challenge for the clinical management of high grade serous ovarian cancer (HGSOC). Here, we demonstrate CX-5461, the first-in-class inhibitor of RNA polymerase I transcription of ribosomal RNA genes (rDNA), induces replication stress and activates the DNA damage response. CX-5461 co-operates with PARPi in exacerbating replication stress and enhances therapeutic efficacy against homologous recombination (HR) DNA repair-deficient HGSOC-patient-derived xenograft (PDX) in vivo. We demonstrate CX-5461 has a different sensitivity spectrum to PARPi involving MRE11-dependent degradation of replication forks. Importantly, CX-5461 exhibits in vivo single agent efficacy in a HGSOC-PDX with reduced sensitivity to PARPi by overcoming replication fork protection. Further, we identify CX-5461-sensitivity gene expression signatures in primary and relapsed HGSOC. We propose CX-5461 is a promising therapy in combination with PARPi in HR-deficient HGSOC and also as a single agent for the treatment of relapsed disease.

A functional genetic screen defines the AKT-induced senescence signaling network.

Exquisite regulation of PI3K/AKT/mTORC1 signaling is essential for homeostatic control of cell growth, proliferation, and survival. Aberrant activation of this signaling network is an early driver of many sporadic human cancers. Paradoxically, sustained hyperactivation of the PI3K/AKT/mTORC1 pathway in nontransformed cells results in cellular senescence, which is a tumor-suppressive mechanism that must be overcome to promote malignant transformation. While oncogene-induced senescence (OIS) driven by excessive RAS/ERK signaling has been well studied, little is known about the mechanisms underpinning the AKT-induced senescence (AIS) response. Here, we utilize a combination of transcriptome and metabolic profiling to identify key signatures required to maintain AIS. We also employ a whole protein-coding genome RNAi screen for AIS escape, validating a subset of novel mediators and demonstrating their preferential specificity for AIS as compared with OIS. As proof of concept of the potential to exploit the AIS network, we show that neurofibromin 1 (NF1) is upregulated during AIS and its ability to suppress RAS/ERK signaling facilitates AIS maintenance. Furthermore, depletion of NF1 enhances transformation of p53-mutant epithelial cells expressing activated AKT, while its overexpression blocks transformation by inducing a senescent-like phenotype. Together, our findings reveal novel mechanistic insights into the control of AIS and identify putative senescence regulators that can potentially be targeted, with implications for new therapeutic options to treat PI3K/AKT/mTORC1-driven cancers.

Somatic mutations in early metazoan genes disrupt regulatory links between unicellular and multicellular genes in cancer.

Extensive transcriptional alterations are observed in cancer, many of which activate core biological processes established in unicellular organisms or suppress differentiation pathways formed in metazoans. Through rigorous, integrative analysis of genomics data from a range of solid tumors, we show many transcriptional changes in tumors are tied to mutations disrupting regulatory interactions between unicellular and multicellular genes within human gene regulatory networks (GRNs). Recurrent point mutations were enriched in regulator genes linking unicellular and multicellular subnetworks, while copy-number alterations affected downstream target genes in distinctly unicellular and multicellular regions of the GRN. Our results depict drivers of tumourigenesis as genes that created key regulatory links during the evolution of early multicellular life, whose dysfunction creates widespread dysregulation of primitive elements of the GRN. Several genes we identified as important in this process were associated with drug response, demonstrating the potential clinical value of our approach.

How the evolution of multicellularity set the stage for cancer.

Neoplastic growth and many of the hallmark properties of cancer are driven by the disruption of molecular networks established during the emergence of multicellularity. Regulatory pathways and molecules that evolved to impose regulatory constraints upon networks established in earlier unicellular organisms enabled greater communication and coordination between the diverse cell types required for multicellularity, but also created liabilities in the form of points of vulnerability in the network that when mutated or dysregulated facilitate the development of cancer. These factors are usually overlooked in genomic analyses of cancer, but understanding where vulnerabilities to cancer lie in the networks of multicellular species would provide important new insights into how core molecular processes and gene regulation change during tumourigenesis. We describe how the evolutionary origins of genes influence their roles in cancer, and how connections formed between unicellular and multicellular genes that act as key regulatory hubs for normal tissue homeostasis can also contribute to malignant transformation when disrupted. Tumours in general are characterised by increased dependence on unicellular processes for survival, and major dysregulation of the control structures imposed on these processes during the evolution of multicellularity. Mounting molecular evidence suggests altered interactions at the interface between unicellular and multicellular genes play key roles in the initiation and progression of cancer. Furthermore, unicellular network regions activated in cancer show high degrees of robustness and plasticity, conferring increased adaptability to tumour cells by supporting effective responses to environmental pressures such as drug exposure. Examining how the links between multicellular and unicellular regions get disrupted in tumours has great potential to identify novel drivers of cancer, and to guide improvements to cancer treatment by identifying more effective therapeutic strategies. Recent successes in targeting unicellular processes by novel compounds underscore the logic of such approaches. Further gains could come from identifying genes at the interface between unicellular and multicellular processes and manipulating the communication between network regions of different evolutionary ages.

Altered interactions between unicellular and multicellular genes drive hallmarks of transformation in a diverse range of solid tumors.

Tumors of distinct tissues of origin and genetic makeup display common hallmark cellular phenotypes, including sustained proliferation, suppression of cell death, and altered metabolism. These phenotypic commonalities have been proposed to stem from disruption of conserved regulatory mechanisms evolved during the transition to multicellularity to control fundamental cellular processes such as growth and replication. Dating the evolutionary emergence of human genes through phylostratigraphy uncovered close association between gene age and expression level in RNA sequencing data from The Cancer Genome Atlas for seven solid cancers. Genes conserved with unicellular organisms were strongly up-regulated, whereas genes of metazoan origin were primarily inactivated. These patterns were most consistent for processes known to be important in cancer, implicating both selection and active regulation during malignant transformation. The coordinated expression of strongly interacting multicellularity and unicellularity processes was lost in tumors. This separation of unicellular and multicellular functions appeared to be mediated by 12 highly connected genes, marking them as important general drivers of tumorigenesis. Our findings suggest common principles closely tied to the evolutionary history of genes underlie convergent changes at the cellular process level across a range of solid cancers. We propose altered activity of genes at the interfaces between multicellular and unicellular regions of human gene regulatory networks activate primitive transcriptional programs, driving common hallmark features of cancer. Manipulation of cross-talk between biological processes of different evolutionary origins may thus present powerful and broadly applicable treatment strategies for cancer.

ProNGF derived from rat sciatic nerves downregulates neurite elongation and axon specification in PC12 cells.

Several reports have shown that a sciatic nerve conditioned media (CM) causes neuronal-like differentiation in PC12 cells. This differentiation is featured by neurite outgrowth, which are exclusively dendrites, without axon or sodium current induction. In previous studies, our group reported that the CM supplemented with a generic inhibitor for tyrosine kinase receptors (k252a) enhanced the CM-induced morphological differentiation upregulating neurite outgrowth, axonal formation and sodium current elicitation. Sodium currents were also induced by depletion of endogenous precursor of nerve growth factorr (proNGF) from the CM (pNGFd-CM). Given that sodium currents, neurite outgrowth and axon specification are important features of neuronal differentiation, in the current manuscript, first we investigated if proNGF was hindering the full PC12 cell neuronal-like differentiation. Second, we studied the effects of exogenous wild type (pNGFwt) and mutated (pNGFmut) proNGF isoforms over sodium currents and whether or not their addition to the pNGFd-CM would prevent sodium current elicitation. Third, we investigated if proNGF was exerting its negative regulation through the sortilin receptor, and for this, the proNGF action was blocked with neurotensin (NT), a factor known to compete with proNGF for sortilin. Thereby, here we show that pNGFd-CM enhanced cell differentiation, cell proportion with long neurites, total neurite length, induced axonal formation and sodium current elicitation. Interestingly, treatment of PC12 cells with wild type or mutated proNGF isoforms elicited sodium currents. Supplementing pNGFd-CM with pNGFmut reduced 35% the sodium currents. On the other hand, pNGFd-CM+pNGFwt induced larger sodium currents than pNGFd-CM. Finally, treatments with CM supplemented with NT showed that sortilin was mediating proNGF negative regulation, since its blocking induced similar effects than the pNGFd-CM treatment. Altogether, our results suggest that proNGF within the CM, is one of the main inhibitors of full neuronal differentiation, acting through sortilin receptor.