Isthmus progenitor cells contribute to homeostatic cellular turnover and support regeneration following intestinal injury.

The currently accepted intestinal epithelial cell organization model proposes that Lgr5+ crypt-base columnar (CBC) cells represent the sole intestinal stem cell (ISC) compartment. However, previous studies have indicated that Lgr5+ cells are dispensable for intestinal regeneration, leading to two major hypotheses: one favoring the presence of a quiescent reserve ISC and the other calling for differentiated cell plasticity. To investigate these possibilities, we studied crypt epithelial cells in an unbiased fashion via high-resolution single-cell profiling. These studies, combined with in vivo lineage tracing, show that Lgr5 is not a specific ISC marker and that stemness potential exists beyond the crypt base and resides in the isthmus region, where undifferentiated cells participate in intestinal homeostasis and regeneration following irradiation (IR) injury. Our results provide an alternative model of intestinal epithelial cell organization, suggesting that stemness potential is not restricted to CBC cells, and neither de-differentiation nor reserve ISC are drivers of intestinal regeneration.

SMARCAL1 is a dual regulator of innate immune signaling and PD-L1 expression that promotes tumor immune evasion.

Genomic instability can trigger cancer-intrinsic innate immune responses that promote tumor rejection. However, cancer cells often evade these responses by overexpressing immune checkpoint regulators, such as PD-L1. Here, we identify the SNF2-family DNA translocase SMARCAL1 as a factor that favors tumor immune evasion by a dual mechanism involving both the suppression of innate immune signaling and the induction of PD-L1-mediated immune checkpoint responses. Mechanistically, SMARCAL1 limits endogenous DNA damage, thereby suppressing cGAS-STING-dependent signaling during cancer cell growth. Simultaneously, it cooperates with the AP-1 family member JUN to maintain chromatin accessibility at a PD-L1 transcriptional regulatory element, thereby promoting PD-L1 expression in cancer cells. SMARCAL1 loss hinders the ability of tumor cells to induce PD-L1 in response to genomic instability, enhances anti-tumor immune responses and sensitizes tumors to immune checkpoint blockade in a mouse melanoma model. Collectively, these studies uncover SMARCAL1 as a promising target for cancer immunotherapy.

A modular master regulator landscape controls cancer transcriptional identity.

Despite considerable efforts, the mechanisms linking genomic alterations to the transcriptional identity of cancer cells remain elusive. Integrative genomic analysis, using a network-based approach, identified 407 master regulator (MR) proteins responsible for canalizing the genetics of individual samples from 20 cohorts in The Cancer Genome Atlas (TCGA) into 112 transcriptionally distinct tumor subtypes. MR proteins could be further organized into 24 pan-cancer, master regulator block modules (MRBs), each regulating key cancer hallmarks and predictive of patient outcome in multiple cohorts. Of all somatic alterations detected in each individual sample, >50% were predicted to induce aberrant MR activity, yielding insight into mechanisms linking tumor genetics and transcriptional identity and establishing non-oncogene dependencies. Genetic and pharmacological validation assays confirmed the predicted effect of upstream mutations and MR activity on downstream cellular identity and phenotype. Thus, co-analysis of mutational and gene expression profiles identified elusive subtypes and provided testable hypothesis for mechanisms mediating the effect of genetic alterations.

Single-cell protein activity analysis identifies recurrence-associated renal tumor macrophages.

Clear cell renal carcinoma (ccRCC) is a heterogeneous disease with a variable post-surgical course. To assemble a comprehensive ccRCC tumor microenvironment (TME) atlas, we performed single-cell RNA sequencing (scRNA-seq) of hematopoietic and non-hematopoietic subpopulations from tumor and tumor-adjacent tissue of treatment-naive ccRCC resections. We leveraged the VIPER algorithm to quantitate single-cell protein activity and validated this approach by comparison to flow cytometry. The analysis identified key TME subpopulations, as well as their master regulators and candidate cell-cell interactions, revealing clinically relevant populations, undetectable by gene-expression analysis. Specifically, we uncovered a tumor-specific macrophage subpopulation characterized by upregulation of TREM2/APOE/C1Q, validated by spatially resolved, quantitative multispectral immunofluorescence. In a large clinical validation cohort, these markers were significantly enriched in tumors from patients who recurred following surgery. The study thus identifies TREM2/APOE/C1Q-positive macrophage infiltration as a potential prognostic biomarker for ccRCC recurrence, as well as a candidate therapeutic target.

An expanded universe of cancer targets.

The characterization of cancer genomes has provided insight into somatically altered genes across tumors, transformed our understanding of cancer biology, and enabled tailoring of therapeutic strategies. However, the function of most cancer alleles remains mysterious, and many cancer features transcend their genomes. Consequently, tumor genomic characterization does not influence therapy for most patients. Approaches to understand the function and circuitry of cancer genes provide complementary approaches to elucidate both oncogene and non-oncogene dependencies. Emerging work indicates that the diversity of therapeutic targets engendered by non-oncogene dependencies is much larger than the list of recurrently mutated genes. Here we describe a framework for this expanded list of cancer targets, providing novel opportunities for clinical translation.

Elucidating Compound Mechanism of Action by Network Perturbation Analysis.

Genome-wide identification of the mechanism of action (MoA) of small-molecule compounds characterizing their targets, effectors, and activity modulators represents a highly relevant yet elusive goal, with critical implications for assessment of compound efficacy and toxicity. Current approaches are labor intensive and mostly limited to elucidating high-affinity binding target proteins. We introduce a regulatory network-based approach that elucidates genome-wide MoA proteins based on the assessment of the global dysregulation of their molecular interactions following compound perturbation. Analysis of cellular perturbation profiles identified established MoA proteins for 70% of the tested compounds and elucidated novel proteins that were experimentally validated. Finally, unknown-MoA compound analysis revealed altretamine, an anticancer drug, as an inhibitor of glutathione peroxidase 4 lipid repair activity, which was experimentally confirmed, thus revealing unexpected similarity to the activity of sulfasalazine. This suggests that regulatory network analysis can provide valuable mechanistic insight into the elucidation of small-molecule MoA and compound similarity.

Identification of causal genetic drivers of human disease through systems-level analysis of regulatory networks.

Identification of driver mutations in human diseases is often limited by cohort size and availability of appropriate statistical models. We propose a framework for the systematic discovery of genetic alterations that are causal determinants of disease, by prioritizing genes upstream of functional disease drivers, within regulatory networks inferred de novo from experimental data. We tested this framework by identifying the genetic determinants of the mesenchymal subtype of glioblastoma. Our analysis uncovered KLHL9 deletions as upstream activators of two previously established master regulators of the subtype, C/EBPß and C/EBPδ. Rescue of KLHL9 expression induced proteasomal degradation of C/EBP proteins, abrogated the mesenchymal signature, and reduced tumor viability in vitro and in vivo. Deletions of KLHL9 were confirmed in > 50% of mesenchymal cases in an independent cohort, thus representing the most frequent genetic determinant of the subtype. The method generalized to study other human diseases, including breast cancer and Alzheimer's disease.

The Immune Landscape of Cancer.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

An extensive microRNA-mediated network of RNA-RNA interactions regulates established oncogenic pathways in glioblastoma.

By analyzing gene expression data in glioblastoma in combination with matched microRNA profiles, we have uncovered a posttranscriptional regulation layer of surprising magnitude, comprising more than 248,000 microRNA (miR)-mediated interactions. These include ∼7,000 genes whose transcripts act as miR "sponges" and 148 genes that act through alternative, nonsponge interactions. Biochemical analyses in cell lines confirmed that this network regulates established drivers of tumor initiation and subtype implementation, including PTEN, PDGFRA, RB1, VEGFA, STAT3, and RUNX1, suggesting that these interactions mediate crosstalk between canonical oncogenic pathways. siRNA silencing of 13 miR-mediated PTEN regulators, whose locus deletions are predictive of PTEN expression variability, was sufficient to downregulate PTEN in a 3'UTR-dependent manner and to increase tumor cell growth rates. Thus, miR-mediated interactions provide a mechanistic, experimentally validated rationale for the loss of PTEN expression in a large number of glioma samples with an intact PTEN locus.

MEF2B mutations lead to deregulated expression of the oncogene BCL6 in diffuse large B cell lymphoma.

MEF2B encodes a transcriptional activator and is mutated in ∼11% of diffuse large B cell lymphomas (DLBCLs) and ∼12% of follicular lymphomas (FLs). Here we found that MEF2B directly activated the transcription of the proto-oncogene BCL6 in normal germinal-center (GC) B cells and was required for DLBCL proliferation. Mutation of MEF2B resulted in enhanced transcriptional activity of MEF2B either through disruption of its interaction with the corepressor CABIN1 or by rendering it insensitive to inhibitory signaling events mediated by phosphorylation and sumoylation. Consequently, the transcriptional activity of Bcl-6 was deregulated in DLBCLs with MEF2B mutations. Thus, somatic mutations of MEF2B may contribute to lymphomagenesis by deregulating BCL6 expression, and MEF2B may represent an alternative target for blocking Bcl-6 activity in DLBCLs.

Entinostat in patients with relapsed or refractory abdominal neuroendocrine tumors.

BACKGROUND: Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) are rare neoplasms with an increasing annual incidence and prevalence. Many are metastatic at presentation or recur following surgical resection and require systemic therapy, for which somatostatin analogs such as octreotide or lanreotide comprise typical first-line therapies. Nonetheless, treatment options remain limited. Epigenetic processes such as histone modifications have been implicated in malignant transformation and progression. In this study, we evaluated the anti-proliferative effects of a histone deacetylase (HDAC) inhibitor, entinostat, which was computationally predicted to show anti-cancer activity, as confirmed in in vitro and in vivo models of GEP-NETs. METHODS: This was a phase II study to evaluate the efficacy and safety of entinostat in patients with relapsed or refractory abdominal NETs. The primary objective was to estimate the objective response rate to entinostat. Additionally, with each patient as his/her own control we estimated the rates of tumor growth prior to enrollment on study and while receiving entinostat. Patients received 5 mg entinostat weekly until disease progression or intolerable toxicity. The dose could be changed to 10 mg biweekly for patients who did not experience grade ≥ 2 treatment-related adverse events (AEs) in cycle 1, but was primarily administered at the starting 5 mg weekly dose. RESULTS: The study enrolled only 5 patients due to early termination by the drug sponsor. The first patient that enrolled had advanced disease and died within days of enrollment before follow-up imaging due to a grade 5 AE unrelated to study treatment and was considered non-evaluable. Best RECIST response for the remaining 4 patients was stable disease (SD) with time on study of 154+, 243, 574, and 741 days. With each patient as his/her own control, rates of tumor growth on entinostat were markedly reduced with rates 20%, 33%, 54%, and 68% of the rates prior to enrollment on study. Toxicities possibly or definitely related to entinostat included grade 2/3 neutrophil count decrease [2/4 (50%)/ 2/4 (50%)], grade 3 hypophosphatemia [1/4, (25%)], grade 1/2 fatigue [1/4 (25%)/ 2/4 (50%)], and other self-limiting grade 1/2 AEs. CONCLUSION: In the treatment of relapsed or refractory abdominal NETs, entinostat 5 mg weekly led to prolonged SD and reduced the rate of tumor growth by 32% to 80% with an acceptable safety profile (ClinicalTrials.gov Identifier: NCT03211988).

Inhibition of the autocrine IL-6-JAK2-STAT3-calprotectin axis as targeted therapy for HR-/HER2+ breast cancers.

HER2-positive (HER2(+)) breast adenocarcinomas are a heterogeneous group in which hormone receptor (HR) status influences therapeutic decisions and patient outcome. By combining genome-wide RNAi screens with regulatory network analysis, we identified STAT3 as a critically activated master regulator of HR(-)/HER2(+) tumors, eliciting tumor dependency in these cells. Mechanistically, HR(-)/HER2(+) cells secrete high levels of the interleukin-6 (IL-6) cytokine, inducing the activation of STAT3, which in turn promotes a second autocrine stimulus to increase S100A8/9 complex (calprotectin) production and secretion. Increased calprotectin levels activate signaling pathways involved in proliferation and resistance. Importantly, we demonstrated that inhibition of the IL-6-Janus kinase 2 (JAK2)-STAT3-calprotectin axis with FDA-approved drugs, alone and in combination with HER2 inhibitors, reduced the tumorigenicity of HR(-)/HER2(+) breast cancers, opening novel targeted therapeutic opportunities.

NaRnEA: An Information Theoretic Framework for Gene Set Analysis.

Gene sets are being increasingly leveraged to make high-level biological inferences from transcriptomic data; however, existing gene set analysis methods rely on overly conservative, heuristic approaches for quantifying the statistical significance of gene set enrichment. We created Nonparametric analytical-Rank-based Enrichment Analysis (NaRnEA) to facilitate accurate and robust gene set analysis with an optimal null model derived using the information theoretic Principle of Maximum Entropy. By measuring the differential activity of ~2500 transcriptional regulatory proteins based on the differential expression of each protein's transcriptional targets between primary tumors and normal tissue samples in three cohorts from The Cancer Genome Atlas (TCGA), we demonstrate that NaRnEA critically improves in two widely used gene set analysis methods: Gene Set Enrichment Analysis (GSEA) and analytical-Rank-based Enrichment Analysis (aREA). We show that the NaRnEA-inferred differential protein activity is significantly correlated with differential protein abundance inferred from independent, phenotype-matched mass spectrometry data in the Clinical Proteomic Tumor Analysis Consortium (CPTAC), confirming the statistical and biological accuracy of our approach. Additionally, our analysis crucially demonstrates that the sample-shuffling empirical null models leveraged by GSEA and aREA for gene set analysis are overly conservative, a shortcoming that is avoided by the newly developed Maximum Entropy analytical null model employed by NaRnEA.

Adult enteric Dclk1-positive glial and neuronal cells reveal distinct responses to acute intestinal injury.

Intestinal ganglionic cells in the adult enteric nervous system (ENS) are continually exposed to stimuli from the surrounding microenvironment and need at times to respond to disturbed homeostasis following acute intestinal injury. The kinase DCLK1 and intestinal Dclk1-positive cells have been reported to contribute to intestinal regeneration. Although Dclk1-positive cells are present in adult enteric ganglia, their cellular identity and response to acute injury have not been investigated in detail. Here, we reveal the presence of distinct Dclk1-tdTom+/CD49b+ glial-like and Dclk1-tdTom+/CD49b- neuronal cell types in adult myenteric ganglia. These ganglionic cells demonstrate distinct patterns of tracing over time yet show a similar expansion in response to elevated serotonergic signaling. Interestingly, Dclk1-tdTom+ glial-like and neuronal cell types appear resistant to acute irradiation injury-mediated cell death. Moreover, Dclk1-tdTom+/CD49b+ glial-like cells show prominent changes in gene expression profiles induced by injury, in contrast to Dclk1-tdTom+/CD49b- neuronal cell types. Finally, subsets of Dclk1-tdTom+/CD49b+ glial-like cells demonstrate prominent overlap with Nestin and p75NTR and strong responses to elevated serotonergic signaling or acute injury. These findings, together with their role in early development and their neural crest-like gene expression signature, suggest the presence of reserve progenitor cells in the adult Dclk1 glial cell lineage.NEW & NOTEWORTHY The kinase DCLK1 identifies glial-like and neuronal cell types in adult murine enteric ganglia, which resist acute injury-mediated cell death yet differ in their cellular response to injury. Interestingly, Dclk1-labeled glial-like cells show prominent transcriptional changes in response to injury and harbor features reminiscent of previously described enteric neural precursor cells. Our data thus add to recently emerging evidence of reserve cellular plasticity in the adult enteric nervous system.

An ID2-dependent mechanism for VHL inactivation in cancer.

Mechanisms that maintain cancer stem cells are crucial to tumour progression. The ID2 protein supports cancer hallmarks including the cancer stem cell state. HIFα transcription factors, most notably HIF2α (also known as EPAS1), are expressed in and required for maintenance of cancer stem cells (CSCs). However, the pathways that are engaged by ID2 or drive HIF2α accumulation in CSCs have remained unclear. Here we report that DYRK1A and DYRK1B kinases phosphorylate ID2 on threonine 27 (Thr27). Hypoxia downregulates this phosphorylation via inactivation of DYRK1A and DYRK1B. The activity of these kinases is stimulated in normoxia by the oxygen-sensing prolyl hydroxylase PHD1 (also known as EGLN2). ID2 binds to the VHL ubiquitin ligase complex, displaces VHL-associated Cullin 2, and impairs HIF2α ubiquitylation and degradation. Phosphorylation of Thr27 of ID2 by DYRK1 blocks ID2-VHL interaction and preserves HIF2α ubiquitylation. In glioblastoma, ID2 positively modulates HIF2α activity. Conversely, elevated expression of DYRK1 phosphorylates Thr27 of ID2, leading to HIF2α destabilization, loss of glioma stemness, inhibition of tumour growth, and a more favourable outcome for patients with glioblastoma.

The miR-424(322)/503 cluster orchestrates remodeling of the epithelium in the involuting mammary gland.

The mammary gland is a very dynamic organ that undergoes continuous remodeling. The critical regulators of this process are not fully understood. Here we identify the microRNA cluster miR-424(322)/503 as an important regulator of epithelial involution after pregnancy. Through the generation of a knockout mouse model, we found that regression of the secretory acini of the mammary gland was compromised in the absence of miR-424(322)/503. Mechanistically, we show that miR-424(322)/503 orchestrates cell life and death decisions by targeting BCL-2 and IGF1R (insulin growth factor-1 receptor). Furthermore, we demonstrate that the expression of this microRNA cluster is regulated by TGF-ß, a well-characterized regulator of mammary involution. Overall, our data suggest a model in which activation of the TGF-ß pathway after weaning induces the transcription of miR-424(322)/503, which in turn down-regulates the expression of key genes. Here, we unveil a previously unknown, multilayered regulation of epithelial tissue remodeling coordinated by the microRNA cluster miR-424(322)/503.

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.

Fimepinostat (CUDC-907) in patients with relapsed/refractory diffuse large B cell and high-grade B-cell lymphoma: report of a phase 2 trial and exploratory biomarker analyses.

Fimepinostat (CUDC-907), a first-in-class oral small-molecule inhibitor of histone deacetylase and phosphatidylinositol 3-kinase, demonstrated efficacy in a phase 1 study of patients with relapsed/refractory (R/R) diffuse large and high-grade B-cell lymphomas (DLBCL/HGBL), particularly those with increased MYC protein expression and/or MYC gene rearrangement/copy number gain (MYC-altered disease). Therefore, a phase 2 study of fimepinostat was conducted in this patient population with 66 eligible patients treated. The primary end-point of overall response (OR) rate for patients with MYC-IHC ≥40% (n = 46) was 15%. Subsequently, exploratory pooled analyses were performed including patients treated on both the phase 1 and 2 studies based upon the presence of MYC-altered disease as well as a biomarker identified by Virtual Inference of Protein activity by Enriched Regulon analysis (VIPER). For these patients with MYC-altered disease (n = 63), the overall response (OR) rate was 22% with seven responding patients remaining on treatment for approximately two years or longer, and VIPER yielded a three-protein biomarker classification with positive and negative predictive values of ≥85%. Prolonged durations of response were achieved by patients with MYC-altered R/R DLBCL/HGBL treated with single-agent fimepinostat. Combination therapies and/or biomarker-based patient selection strategies may lead to higher response rates in future clinical trials.