Elucidation and Pharmacologic Targeting of Master Regulator Dependencies in Coexisting Diffuse Midline Glioma Subpopulations.

Diffuse Midline Gliomas (DMGs) are universally fatal, primarily pediatric malignancies affecting the midline structures of the central nervous system. Despite decades of clinical trials, treatment remains limited to palliative radiation therapy. A major challenge is the coexistence of molecularly distinct malignant cell states with potentially orthogonal drug sensitivities. To address this challenge, we leveraged established network-based methodologies to elucidate Master Regulator (MR) proteins representing mechanistic, non-oncogene dependencies of seven coexisting subpopulations identified by single-cell analysis-whose enrichment in essential genes was validated by pooled CRISPR/Cas9 screens. Perturbational profiles of 372 clinically relevant drugs helped identify those able to invert the activity of subpopulation-specific MRs for follow-up in vivo validation. While individual drugs predicted to target individual subpopulations-including avapritinib, larotrectinib, and ruxolitinib-produced only modest tumor growth reduction in orthotopic models, systemic co-administration induced significant survival extension, making this approach a valuable contribution to the rational design of combination therapy.

Tff2 defines transit-amplifying pancreatic acinar progenitors that lack regenerative potential and are protective against Kras-driven carcinogenesis.

While adult pancreatic stem cells are thought not to exist, it is now appreciated that the acinar compartment harbors progenitors, including tissue-repairing facultative progenitors (FPs). Here, we study a pancreatic acinar population marked by trefoil factor 2 (Tff2) expression. Long-term lineage tracing and single-cell RNA sequencing (scRNA-seq) analysis of Tff2-DTR-CreERT2-targeted cells defines a transit-amplifying progenitor (TAP) population that contributes to normal homeostasis. Following acute and chronic injury, Tff2+ cells, distinct from FPs, undergo depopulation but are eventually replenished. At baseline, oncogenic KrasG12D-targeted Tff2+ cells are resistant to PDAC initiation. However, KrasG12D activation in Tff2+ cells leads to survival and clonal expansion following pancreatitis and a cancer stem/progenitor cell-like state. Selective ablation of Tff2+ cells prior to KrasG12D activation in Mist1+ acinar or Dclk1+ FP cells results in enhanced tumorigenesis, which can be partially rescued by adenoviral Tff2 treatment. Together, Tff2 defines a pancreatic TAP population that protects against Kras-driven carcinogenesis.

Ras-dependent activation of BMAL2 regulates hypoxic metabolism in pancreatic cancer.

To identify novel drivers of malignancy in pancreatic ductal adenocarcinoma (PDAC), we employed regulatory network analysis, which calculates the activity of transcription factors and other regulatory proteins based on the integrated expression of their positive and negative target genes. We generated a regulatory network for the malignant epithelial cells of human PDAC using gene expression data from a set of 197 laser capture microdissected human PDAC samples and 45 low-grade precursors, for which we had matched histopathological, clinical, and epidemiological annotation. We then identified the most highly activated and repressed regulatory proteins (e.g. master regulators or MRs) associated with four malignancy phenotypes: precursors vs. PDAC (initiation), low-grade vs. high grade histopathology (progression), survival post resection, and association with KRAS activity. Integrating across these phenotypes, the top MR of PDAC malignancy was found to be BMAL2, a member of the PAS family of bHLH transcription factors. Although the canonical function of BMAL2 is linked to the circadian rhythm protein CLOCK, annotation of BMAL2 target genes highlighted a potential role in hypoxia response. We previously demonstrated that PDAC is hypovascularized and hypoperfused, and here show that PDAC from the genetically engineered KPC model exists in a state of extreme hypoxia, with a partial oxygen pressure of <1mmHg. Given the close homology of BMAL2 to HIF1ß (ARNT) and its potential to heterodimerize with HIF1A and HIF2A, we investigated whether BMAL2 plays a role in the hypoxic response of PDAC. Indeed, BMAL2 controlled numerous hypoxia response genes and could be inhibited following treatment with multiple RAF, MEK, and ERK inhibitors, validating its association with RAS activity. Knockout of BMAL2 in four human PDAC cell lines led to defects in growth and invasion in the setting of hypoxia. Strikingly, BMAL2 null cells failed to induce glycolysis upon exposure to severe hypoxia and this was associated with a loss of expression of the glycolytic enzyme LDHA. Moreover, HIF1A was no longer stabilized under hypoxia in BMAL2 knockout cells. By contrast, HIF2A was hyper-stabilized under hypoxia, indicating a dysregulation of hypoxia metabolism in response to BMAL2 loss. We conclude that BMAL2 is a master regulator of hypoxic metabolism in PDAC, serving as a molecular switch between the disparate metabolic roles of HIF1A- and HIF2A-dependent hypoxia responses.

A Patient-to-Model-to-Patient (PMP) Cancer Drug Discovery Protocol for Identifying and Validating Therapeutic Agents Targeting Tumor Regulatory Architecture.

The current Achilles heel of cancer drug discovery is the inability to forge precise and predictive connections among mechanistic drivers of the cancer cell state, therapeutically significant molecular targets, effective drugs, and responsive patient subgroups. Although advances in molecular biology have helped identify molecular markers and stratify patients into molecular subtypes, these associational strategies typically fail to provide a mechanistic rationale to identify cancer vulnerabilities. Recently, integrative systems biology methodologies have been used to reverse engineer cellular networks and identify master regulators (MRs), proteins whose activity is both necessary and sufficient to implement phenotypic states under physiological and pathological conditions, which are organized into highly interconnected regulatory modules called tumor checkpoints. Because of their functional relevance, MRs represent ideal pharmacological targets and biomarkers. Here, we present a six-step patient-to-model-to-patient protocol that employs computational and experimental methodologies to reconstruct and interrogate the regulatory logic of human cancer cells for identifying and therapeutically targeting the tumor checkpoint with novel as well as existing pharmacological agents. This protocol systematically identifies, from specific patient tumor samples, the MRs that comprise the tumor checkpoint. Then, it identifies in vitro and in vivo models that, by recapitulating the patient's tumor checkpoint, constitute the appropriate cell lines and xenografts to further elucidate the tissue context-specific drug mechanism of action (MOA) and permit precise, biomarker-based preclinical validations of drug efficacy. The combination of determination of a drug's context-specific MOA and precise identification of patients' tumor checkpoints provides a personalized, mechanism-based biomarker to enrich prospective clinical trials with patients likely to respond. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC.

EZH2-Mediated H3K27me3 Targets Transcriptional Circuits of Neuronal Differentiation.

The Polycomb Repressive Complex 2 (PRC2) plays important roles in the epigenetic regulation of cellular development and differentiation through H3K27me3-dependent transcriptional repression. Aberrant PRC2 activity has been associated with cancer and neurodevelopmental disorders, particularly with respect to the malfunction of sits catalytic subunit EZH2. Here, we investigated the role of the EZH2-mediated H3K27me3 apposition in neuronal differentiation. We made use of a transgenic mouse model harboring Ezh2 conditional KO alleles to derive embryonic stem cells and differentiate them into glutamatergic neurons. Time course transcriptomics and epigenomic analyses of H3K27me3 in absence of EZH2 revealed a significant dysregulation of molecular networks affecting the glutamatergic differentiation trajectory that resulted in: (i) the deregulation of transcriptional circuitries related to neuronal differentiation and synaptic plasticity, in particular LTD, as a direct effect of EZH2 loss and (ii) the appearance of a GABAergic gene expression signature during glutamatergic neuron differentiation. These results expand the knowledge about the molecular pathways targeted by Polycomb during glutamatergic neuron differentiation.

A cell-of-origin epigenetic tracer reveals clinically distinct subtypes of high-grade serous ovarian cancer.

BACKGROUND: High-grade serous ovarian cancer (HGSOC) is a major unmet need in oncology. The remaining uncertainty on its originating tissue has hampered the discovery of molecular oncogenic pathways and the development of effective therapies. METHODS: We used an approach based on the retention in tumors of a DNA methylation trace (OriPrint) that distinguishes the two putative tissues of origin of HGSOC, the fimbrial (FI) and ovarian surface epithelia (OSE), to stratify HGSOC by several clustering methods, both linear and non-linear. The identified tumor subtypes (FI-like and OSE-like HGSOC) were investigated at the RNAseq level to stratify an in-house cohort of macrodissected HGSOC FFPE samples to derive overall and disease-free survival and identify specific transcriptional alterations of the two tumor subtypes, both by classical differential expression and weighted correlation network analysis. We translated our strategy to published datasets and verified the co-occurrence of previously described molecular classification of HGSOC. We performed cytokine analysis coupled to immune phenotyping to verify alterations in the immune compartment associated with HGSOC. We identified genes that are both differentially expressed and methylated in the two tumor subtypes, concentrating on PAX8 as a bona fide marker of FI-like HGSOC. RESULTS: We show that: - OriPrint is a robust DNA methylation tracer that exposes the tissue of origin of HGSOC. - The tissue of origin of HGSOC is the main determinant of DNA methylation variance in HGSOC. - The tissue of origin is a prognostic factor for HGSOC patients. - FI-like and OSE-like HGSOC are endowed with specific transcriptional alterations that impact patients' prognosis. - OSE-like tumors present a more invasive and immunomodulatory phenotype, compatible with its worse prognostic impact. - Among genes that are differentially expressed and regulated in FI-like and OSE-like HGSOC, PAX8 is a bona fide marker of FI-like tumors. CONCLUSIONS: Through an integrated approach, our work demonstrates that both FI and OSE are possible origins for human HGSOC, whose derived subtypes are both molecularly and clinically distinct. These results will help define a new roadmap towards rational, subtype-specific therapeutic inroads and improved patients' care.

Single-Cell Genomics Reveals a Novel Cell State During Smooth Muscle Cell Phenotypic Switching and Potential Therapeutic Targets for Atherosclerosis in Mouse and Human.

BACKGROUND: Smooth muscle cells (SMCs) play significant roles in atherosclerosis via phenotypic switching, a pathological process in which SMC dedifferentiation, migration, and transdifferentiation into other cell types. Yet how SMCs contribute to the pathophysiology of atherosclerosis remains elusive. METHODS: To reveal the trajectories of SMC transdifferentiation during atherosclerosis and to identify molecular targets for disease therapy, we combined SMC fate mapping and single-cell RNA sequencing of both mouse and human atherosclerotic plaques. We also performed cell biology experiments on isolated SMC-derived cells, conducted integrative human genomics, and used pharmacological studies targeting SMC-derived cells both in vivo and in vitro. RESULTS: We found that SMCs transitioned to an intermediate cell state during atherosclerosis, which was also found in human atherosclerotic plaques of carotid and coronary arteries. SMC-derived intermediate cells, termed "SEM" cells (stem cell, endothelial cell, monocyte), were multipotent and could differentiate into macrophage-like and fibrochondrocyte-like cells, as well as return toward the SMC phenotype. Retinoic acid (RA) signaling was identified as a regulator of SMC to SEM cell transition, and RA signaling was dysregulated in symptomatic human atherosclerosis. Human genomics revealed enrichment of genome-wide association study signals for coronary artery disease in RA signaling target gene loci and correlation between coronary artery disease risk alleles and repressed expression of these genes. Activation of RA signaling by all-trans RA, an anticancer drug for acute promyelocytic leukemia, blocked SMC transition to SEM cells, reduced atherosclerotic burden, and promoted fibrous cap stability. CONCLUSIONS: Integration of cell-specific fate mapping, single-cell genomics, and human genetics adds novel insights into the complexity of SMC biology and reveals regulatory pathways for therapeutic targeting of SMC transitions in atherosclerotic cardiovascular disease.

The Host Cell ViroCheckpoint: Identification and Pharmacologic Targeting of Novel Mechanistic Determinants of Coronavirus-Mediated Hijacked Cell States.

Most antiviral agents are designed to target virus-specific proteins and mechanisms rather than the host cell proteins that are critically dysregulated following virus-mediated reprogramming of the host cell transcriptional state. To overcome these limitations, we propose that elucidation and pharmacologic targeting of host cell Master Regulator proteins-whose aberrant activities govern the reprogramed state of coronavirus-infected cells-presents unique opportunities to develop novel mechanism-based therapeutic approaches to antiviral therapy, either as monotherapy or as a complement to established treatments. Specifically, we propose that a small module of host cell Master Regulator proteins (ViroCheckpoint) is hijacked by the virus to support its efficient replication and release. Conventional methodologies are not well suited to elucidate these potentially targetable proteins. By using the VIPER network-based algorithm, we successfully interrogated 12h, 24h, and 48h signatures from Calu-3 lung adenocarcinoma cells infected with SARS-CoV, to elucidate the time-dependent reprogramming of host cells and associated Master Regulator proteins. We used the NYS CLIA-certified Darwin OncoTreat algorithm, with an existing database of RNASeq profiles following cell perturbation with 133 FDA-approved and 195 late-stage experimental compounds, to identify drugs capable of virtually abrogating the virus-induced Master Regulator signature. This approach to drug prioritization and repurposing can be trivially extended to other viral pathogens, including SARS-CoV-2, as soon as the relevant infection signature becomes available.

Cross-Species Single-Cell Analysis of Pancreatic Ductal Adenocarcinoma Reveals Antigen-Presenting Cancer-Associated Fibroblasts.

Cancer-associated fibroblasts (CAF) are major players in the progression and drug resistance of pancreatic ductal adenocarcinoma (PDAC). CAFs constitute a diverse cell population consisting of several recently described subtypes, although the extent of CAF heterogeneity has remained undefined. Here we use single-cell RNA sequencing to thoroughly characterize the neoplastic and tumor microenvironment content of human and mouse PDAC tumors. We corroborate the presence of myofibroblastic CAFs and inflammatory CAFs and define their unique gene signatures in vivo. Moreover, we describe a new population of CAFs that express MHC class II and CD74, but do not express classic costimulatory molecules. We term this cell population "antigen-presenting CAFs" and find that they activate CD4+ T cells in an antigen-specific fashion in a model system, confirming their putative immune-modulatory capacity. Our cross-species analysis paves the way for investigating distinct functions of CAF subtypes in PDAC immunity and progression. SIGNIFICANCE: Appreciating the full spectrum of fibroblast heterogeneity in pancreatic ductal adenocarcinoma is crucial to developing therapies that specifically target tumor-promoting CAFs. This work identifies MHC class II-expressing CAFs with a capacity to present antigens to CD4+ T cells, and potentially to modulate the immune response in pancreatic tumors.See related commentary by Belle and DeNardo, p. 1001.This article is highlighted in the In This Issue feature, p. 983.

The ion channels and transporters gene expression profile indicates a shift in excitability and metabolisms during malignant progression of Follicular Lymphoma.

The definition of the gene expression profile of genes encoding Ion Channels and Transporters (ICT-GEP) represents a novel and attracting aspect in cancer. We determined the ICT-GEP of Follicular Lymphoma (FL), and compared it with that of the more aggressive Diffuse Large B Cell Lymphoma (DLBCL). cDNA microarray data were collected both from patients enrolled for this study, and from public datasets. In FL the ICT-GEP indicated the overexpression of both the K+ channel encoding gene KCNN4, and SLC2A1, which encodes the Glut1 glucose transporter. SLC2A1 turned out to represent the hub of a functional network, connecting channels and transporters in FL. Relapsed FL patients were characterised by 38 differentially expressed ICT genes, among which ATP9A, SLC2A1 and KCNN4 were under-expressed, indicating a down-regulation of both excitability and glycolysis. A completely different profile of K+ channel encoding genes emerged in DLBCL accompanied by the over-expression of the fatty acid transporter-encoding gene SLC27A1 as well as of the metabolism regulator NCoR1. This indicates a change in excitability and a shift towards an oxidative metabolism in DLBCL. Overall, the ICT-GEP may contribute to identifying novel lymphoma biomarkers related to excitability and metabolic pathways, with particular relevance for drug resistant, relapsed FL.

Experimental microdissection enables functional harmonisation of pancreatic cancer subtypes.

OBJECTIVE: Pancreatic ductal adenocarcinoma (PDA) has among the highest stromal fractions of any cancer and this has complicated attempts at expression-based molecular classification. The goal of this work is to profile purified samples of human PDA epithelium and stroma and examine their respective contributions to gene expression in bulk PDA samples. DESIGN: We used laser capture microdissection (LCM) and RNA sequencing to profile the expression of 60 matched pairs of human PDA malignant epithelium and stroma samples. We then used these data to train a computational model that allowed us to infer tissue composition and generate virtual compartment-specific expression profiles from bulk gene expression cohorts. RESULTS: Our analysis found significant variation in the tissue composition of pancreatic tumours from different public cohorts. Computational removal of stromal gene expression resulted in the reclassification of some tumours, reconciling functional differences between different cohorts. Furthermore, we established a novel classification signature from a total of 110 purified human PDA stroma samples, finding two groups that differ in the extracellular matrix-associated and immune-associated processes. Lastly, a systematic evaluation of cross-compartment subtypes spanning four patient cohorts indicated partial dependence between epithelial and stromal molecular subtypes. CONCLUSION: Our findings add clarity to the nature and number of molecular subtypes in PDA, expand our understanding of global transcriptional programmes in the stroma and harmonise the results of molecular subtyping efforts across independent cohorts.

Polycomb dysregulation in gliomagenesis targets a Zfp423-dependent differentiation network.

Malignant gliomas constitute one of the most significant areas of unmet medical need, owing to the invariable failure of surgical eradication and their marked molecular heterogeneity. Accumulating evidence has revealed a critical contribution by the Polycomb axis of epigenetic repression. However, a coherent understanding of the regulatory networks affected by Polycomb during gliomagenesis is still lacking. Here we integrate transcriptomic and epigenomic analyses to define Polycomb-dependent networks that promote gliomagenesis, validating them both in two independent mouse models and in a large cohort of human samples. We find that Polycomb dysregulation in gliomagenesis affects transcriptional networks associated with invasiveness and de-differentiation. The dissection of these networks uncovers Zfp423 as a critical Polycomb-dependent transcription factor whose silencing negatively impacts survival. The anti-gliomagenic activity of Zfp423 requires interaction with the SMAD proteins within the BMP signalling pathway, pointing to a novel synergic circuit through which Polycomb inhibits BMP signalling.

Deregulation of ion channel and transporter encoding genes in pediatric gliomas.

Brain tumors, including the majority gliomas, are the leading cause of cancer-related death in children. World Health Organization has divided pediatric brain tumors into different grades and, based upon cDNA microarray data identifying gene expression profiles (GEPs), it has become evident in the last decade that the various grades involve different types of genetic alterations. However, it is not known whether ion channel and transporter genes, intimately involved in brain functioning, are associated with such GEPs. We determined the GEPs in an available cohort of 10 pediatric brain tumors initially by comparing the data obtained from four primary tumor samples and corresponding short-term cultures. The correspondence between the two types of samples was statistically significant. We then performed bioinformatic analyses on those samples (a total of nine) which corresponded to tumors of glial origin, either tissues or cell cultures, depending on the best "RNA integrity number." We used R software to evaluate the genes which were differentially expressed (DE) in gliomas compared with normal brain. Applying a p-value below 0.01 and fold change ≥4, led to identification of 2284 DE genes. Through a Functional Annotation Analysis (FAA) using the NIH-DAVID software, the DE genes turned out to be associated mainly with: immune/inflammatory response, cell proliferation and survival, cell adhesion and motility, neuronal phenotype, and ion transport. We have shown that GEPs of pediatric brain tumors can be studied using either primary tumor samples or short-term cultures with similar results. From FAA, we concluded that, among DE genes, pediatric gliomas show a strong deregulation of genes related to ion channels and transporters.

Deterministic and stochastic aspects of VEGF-A production and the cooperative behavior of tumoral cell colony.

A model is proposed to study the process of hypoxia-induced angiogenesis in cancer cells. The model accounts for the role played by the vascular endothelial growth factor (VEGF)-A in regulating the oxygen intake. VEGF-A is dynamically controlled by the HIF-1α concentration. If not degraded, HIF-1α can bind to the subunit termed HIF-1ß and so experience translocation to the nucleus, to exert its proper transcriptional activity. The delicate balance between these opposing tendencies translates into the emergence of distinct macroscopic behaviors in terms of the associated molecular concentrations that we here trace back to normoxia, hypoxia and death regimes. These aspects are firstly analyzed with reference to the ideal mean-field scenario. Stochastic fluctuations are also briefly discussed and shown to seed a cooperative interaction among cellular units, competing for the same oxygen reservoir.