Gene Program Research Articles

Single-cell profiling aligns CD56bright and cytomegalovirus-induced adaptive natural killer cells to a naïve-memory relationship

Development of antigen-specific memory upon pathogen exposure is a hallmark of the adaptive immune system. While natural killer (NK) cells are considered part of the innate immune system, humans exposed to the chronic viral pathogen cytomegalovirus (CMV) often possess a distinct NK cell population lacking in individuals who have not been exposed, termed “adaptive” NK cells. To identify the “naïve” population from which this “memory” population derives, we performed phenotypic, transcriptional, and functional profiling of NK cell subsets. We identified immature precursors to the Adaptive NK cells that are equally present in both CMV+ and CMV- individuals, resolved an Adaptive transcriptional state distinct from most mature NK cells and sharing a common gene program with the immature CD56bright population, and demonstrated retention of proliferative capacity and acquisition of superior IFNγ production in the Adaptive population. Furthermore, we distinguish the CD56bright and Adaptive NK populations by expression of the transcription factor CXXC5, positioning these memory NK cells at the inflection point between innate and adaptive lymphocytes.

HMX3 is a critical vulnerability in MECOM-negative KMT2A::MLLT3 acute myelomonocytic leukemia.

KMT2A::MLLT3 acute myelomonocytic leukemia (AML) comes in two clinically and biologically different subtypes. One is characterized by inferior outcome, older age, and MECOM oncogene expression. The other is mainly observed in children and young adults, associates with better clinical outcome, but lacks MECOM. To identify cell fate determining transcription factors downstream of KMT2A::MLLT3, we applied a bioinformatic algorithm that integrates gene and enhancer expression from primary MECOM-positive and -negative KMT2A::MLLT3 AML samples. This identified MECOM to be most influential in the MECOM-positive group, while neuronal transcription factor HMX3 was most influential in the MECOM-negative group. In large AML cohorts, HMX3 expression associated with a unique gene expression profile, younger age (p < 0.002) and KMT2A-rearranged and KAT6A-CREBBP leukemia (p < 0.00001). HMX3 was not expressed in other major genetic risk groups and healthy blood cells. RNA-sequencing analyses following forced HMX3 expression in healthy CD34+ cells and its silencing in KMT2A::MLT3 cells showed that HMX3 drives cancer-associated E2F and MYC gene programs (p < 0.001). HMX3 expression in healthy CD34+ cells blocked monocytic but not granulocytic colony formation. Strikingly, HMX3 silencing in KMT2A::MLLT3 patient cells resulted in cell cycle arrest, monocytic differentiation and apoptosis. Thus, the neuronal transcription factor HMX3 is a leukemia-specific vulnerability in KMT2A::MLLT3 AML.

Differential activity of MAPK signalling defines fibroblast subtypes in pancreatic cancer

Fibroblast heterogeneity is increasingly recognised across cancer conditions. Given their important contribution to disease progression, mapping fibroblasts’ heterogeneity is critical to devise effective anti-cancer therapies. Cancer-associated fibroblasts (CAFs) represent the most abundant cell population in pancreatic ductal adenocarcinoma (PDAC). Whether CAF phenotypes are differently specified by PDAC cell lineages remains to be elucidated. Here, we reveal an important role for the MAPK signalling pathway in defining PDAC CAF phenotypes. We show that epithelial MAPK activity promotes the myofibroblastic differentiation of CAFs by sustaining the expression and secretion of TGF-β1. We integrate single-cell profiling of post-perturbation transcriptional responses from mouse models with cellular and spatial profiles of human tissues to define a MAPKhigh CAF (mapCAF) phenotype. We show that this phenotype associates with basal-like tumour cells and reduced frequency of CD8+ T cells. In addition to elevated MAPK activity, this mapCAF phenotype is characterized by TGF-β signalling, hypoxia responsive signatures, and immunoregulatory gene programs. Furthermore, the mapCAF signature is enriched in myofibroblastic CAFs from various cancer conditions and correlates with reduced response to immune checkpoint inhibition in melanoma. Altogether, our data expand our knowledge on CAF phenotype heterogeneity and reveal a potential strategy for targeting myofibroblastic CAFs in vivo.

First person – Ayyappa Raja Desingu Rajan

ABSTRACT First Person is a series of interviews with the first authors of a selection of papers published in Disease Models &amp; Mechanisms, helping researchers promote themselves alongside their papers. Ayyappa Raja Desingu Rajan is first author on ‘ Generation of a zebrafish neurofibromatosis model via inducible knockout of nf2a/b’, published in DMM. Ayyappa Raja is a postdoctoral researcher in the lab of Marianne Bronner at California Institute of Technology, Pasadena, CA, USA. Using the neural crest as a model, Ayyappa Raja investigates lineage decisions that contribute to specific tissue types, such as the meninges, and how developmental gene programs can go array in cancers such as meningiomas and neurofibromatosis.

The universal role of adaptive transcription in health and disease.

In animals, adaptive transcription is a crucial mechanism to connect environmental stimulation to changes in gene expression and subsequent organism remodeling. Adaptive transcriptional programs involving molecules such as CREB, SRF, MEF2, FOS, and EGR1 are central to a wide variety of organism functions, including learning and memory, immune system plasticity, and muscle hypertrophy, and their activation increases cellular resilience and prevents various diseases. Yet, they also form the basis for many maladaptive processes and are involved in the progression of addiction, depression, cancer, cardiovascular disorders, autoimmune conditions, and metabolic dysfunction among others and are thus prime examples for mediating the adaptation-maladaptation dilemma. They are implicated in the therapeutic effects of major treatment modalities such as antidepressants and can have negative effects on treatment, for example, contributing to therapy resistance in cancer. This review examines the universal role of adaptive transcription as a mechanism for the induction of adaptive cell state transitions in health and disease and explores how many medical disorders can be conceptualized as caused by errors in cellular adaptation goals. It also considers the underlying principles in the basic structure of adaptive gene programs such as their division into a core and a directional program. Finally, it analyses how one might best reprogram cells via targeting of adaptive transcription in combination with complex stimulation patterns to leverage endogenous cellular reprogramming dynamics and achieve optimal health of the whole organism.

Functional profiling of murine glioma models highlights targetable immune evasion phenotypes

Cancer-intrinsic immune evasion mechanisms and pleiotropy are a barrier to cancer immunotherapy. This is apparent in certain highly fatal cancers, including high-grade gliomas and glioblastomas (GBM). In this study, we evaluated two murine syngeneic glioma models (GL261 and CT2A) as preclinical models for human GBM using functional genetic screens, single-cell transcriptomics and machine learning approaches. Through CRISPR genome-wide co-culture killing screens with various immune cells (cytotoxic T cells, natural killer cells, and macrophages), we identified three key cancer-intrinsic evasion mechanisms: NFκB signaling, autophagy/endosome machinery, and chromatin remodeling. Additional fitness screens identified dependencies in murine gliomas that partially recapitulated those seen in human GBM (e.g., UFMylation). Our single-cell analyses showed that different glioma models exhibited distinct immune infiltration patterns and recapitulated key immune gene programs observed in human GBM, including hypoxia, interferon, and TNF signaling. Moreover, in vivo orthotopic tumor engraftment was associated with phenotypic shifts and changes in proliferative capacity, with murine tumors recapitulating the intratumoral heterogeneity observed in human GBM, exhibiting propensities for developmental- and mesenchymal-like phenotypes. Notably, we observed common transcription factors and cofactors shared with human GBM, including developmental (Nfia and Tcf4), mesenchymal (Prrx1 and Wwtr1), as well as cycling-associated genes (Bub3, Cenpa, Bard1, Brca1, and Mis18bp1). Perturbation of these genes led to reciprocal phenotypic shifts suggesting intrinsic feedback mechanisms that balance in vivo cellular states. Finally, we used a machine-learning approach to identify two distinct immune evasion gene programs, one of which represents a clinically-relevant phenotype and delineates a subpopulation of stem-like glioma cells that predict response to immune checkpoint inhibition in human patients. This comprehensive characterization helps bridge the gap between murine glioma models and human GBM, providing valuable insights for future therapeutic development.

Dimethyl fumarate ameliorates chronic graft-versus-host disease by inhibiting Tfh differentiation via Nrf2.

Chronic graft-versus-host disease (cGVHD), characterized by chronic tissue inflammation and fibrosis involving multiple organs, remains a major complication after allogeneic hematopoietic stem cell transplantation (allo-HSCT). Dimethyl fumarate (DMF) is an anti-inflammatory drug approved for the treatment of multiple sclerosis and psoriasis. We previously reported that DMF effectively inhibits acute GVHD (aGVHD) while preserving the graft-versus-leukemia effect. However, the role of DMF in cGVHD progression remains unknown. Here, we found that DMF administration significantly suppresses follicular helper T cell (Tfh) differentiation, and germinal center formation and alleviates disease severity in different murine cGVHD models. Mechanistically, DMF treatment downregulates IL-21 transcription by activation of Nrf2, thus orchestrating Tfh-related gene programs both in mice and humans. The inhibitory role of DMF on Tfh cell differentiation was diminished in Nrf2 deficient T cells. Importantly, the therapeutic potential of DMF in clinical cGVHD has been validated in human data whereby DMF effectively reduces IL-21 production and Tfh cell generation in peripheral blood mononuclear cells from active cGVHD patients and further attenuates xenograft GVHD. Collectively, our findings reveal that DMF potently inhibits cGVHD development by repressing Tfh cell differentiation via Nrf2, paving the way for the treatment of cGVHD in the clinic.

Cardiac contractility modulation: from molecular patterns to tailored treatment in heart failure subgroups.

Cardiac contractility modulation (CCM) signals are non-excitatory signals that are applied during the myocyte's absolute refractory period. These signals have been demonstrated to have an inotropic effect without increasing myocardial oxygen consumption. This has been observed in both preclinical animal studies and randomized clinical trials. CCM influences the expression of various genes that are abnormally expressed in heart failure: it reverses fetal myocyte gene programming associated with heart failure and regulates the expression of genes associated with calcium cycling and myocardial contractile machinery. Clinical investigations have primarily focused on patients with heart failure and normal QRS duration where CCM has demonstrated its safety and effectiveness in reducing heart failure-related hospitalizations, as well as improving symptoms, functional capacity, and overall quality of life. Currently, for individuals experiencing symptomatic heart failure with an ejection fraction ranging from 25% to 45% and a QRS duration of less than 130 ms, who are not suitable candidates for cardiac resynchronization therapy, CCM offers a viable treatment option. Even though promising results in specific HF subgroups have been published, further studies are needed to understand the role of CCM in tailored treatment for heart failure. Moreover, the role of multimodality imaging in lead placement and prognostic stratification in CCM patients should be further investigated. This review aims to summarize the main pathophysiological evidence related to the use of CCM and to highlight its role as a possible additional weapon in tailored treatment for specific subgroups of patients with heart failure.

Melanin-like nanoparticles slow cyst growth in ADPKD by dual inhibition of oxidative stress and CREB.

Melanin-like nanoparticles (MNPs) have recently emerged as valuable agents in antioxidant therapy due to their excellent biocompatibility and potent capacity to scavenge various reactive oxygen species (ROS). However, previous studies have mainly focused on acute ROS-related diseases, leaving a knowledge gap regarding their potential in chronic conditions. Furthermore, apart from their well-established antioxidant effects, it remains unclear whether MNPs target other intracellular molecular pathways. In this study, we synthesized ultra-small polyethylene glycol-incorporated Mn2+-chelated MNP (MMPP). We found that MMPP traversed the glomerular filtration barrier and specifically accumulated in renal tubules. Autosomal dominant polycystic kidney disease (ADPKD) is a chronic genetic disorder closely associated with increased oxidative stress and featured by the progressive enlargement of cysts originating from various segments of the renal tubules. Treatment with MMPP markedly attenuated oxidative stress levels, inhibited cyst growth, thereby improving renal function. Interestingly, we found that MMPP effectively inhibits a cyst-promoting gene program downstream of the cAMP-CREB pathway, a crucial signaling pathway implicated in ADPKD progression. Mechanistically, we observed that MMPP directly binds to the bZIP DNA-binding domain of CREB, leading to competitive inhibition of CREB's DNA binding ability and subsequent reduction in CREB target gene expression. In summary, our findings identify an intracellular target of MMPP and demonstrate its potential for treating ADPKD by simultaneously targeting oxidative stress and CREB transcriptional activity.

Modeling gene interactions in polygenic prediction via geometric deep learning.

Polygenic risk score (PRS) is a widely-used approach for predicting individuals' genetic risk of complex diseases, playing a pivotal role in advancing precision medicine. Traditional PRS methods, predominantly following a linear structure, often fall short in capturing the intricate relationships between genotype and phenotype. In this study, we present PRS-Net, an interpretable geometric deep learning-based framework that effectively models the nonlinearity of biological systems for enhanced disease prediction and biological discovery. PRS-Net begins by deconvoluting the genome-wide PRS at the single-gene resolution, and then explicitly encapsulates gene-gene interactions leveraging a graph neural network (GNN) for genetic risk prediction, enabling a systematic characterization of molecular interplay underpinning diseases. An attentive readout module is introduced to facilitate model interpretation. Extensive tests across multiple complex traits and diseases demonstrate the superior prediction performance of PRS-Net compared to conventional PRS methods. The interpretability of PRS-Net further enhances the identification of disease-relevant genes and gene programs. PRS-Net provides a potent tool for concurrent genetic risk prediction and biological discovery for complex diseases.

Cutaneous T cell lymphoma atlas reveals malignant TH2 cells supported by a B cell-rich tumor microenvironment

Cutaneous T cell lymphoma (CTCL) is a potentially fatal clonal malignancy of T cells primarily affecting the skin. The most common form of CTCL, mycosis fungoides, can be difficult to diagnose, resulting in treatment delay. We performed single-cell and spatial transcriptomics analysis of skin from patients with mycosis fungoides-type CTCL and an integrated comparative analysis with human skin cell atlas datasets from healthy and inflamed skin. We revealed the co-optation of T helper 2 (TH2) cell-immune gene programs by malignant CTCL cells and modeling of the tumor microenvironment to support their survival. We identified MHC-II+ fibroblasts and dendritic cells that can maintain TH2 cell-like tumor cells. CTCL tumor cells are spatially associated with B cells, forming tertiary lymphoid structure-like aggregates. Finally, we validated the enrichment of B cells in CTCL and its association with disease progression across three independent patient cohorts. Our findings provide diagnostic aids, potential biomarkers for disease staging and therapeutic strategies for CTCL.

Redistribution of PU.1 partner transcription factor RUNX1 binding secures cell survival during leukemogenesis

Transcription factors (TFs) orchestrating lineage-development often control genes required for cellular survival. However, it is not well understood how cells survive when such TFs are lost, for example in cancer. PU.1 is an essential TF for myeloid fate, and mice with downregulated PU.1 levels develop acute myeloid leukemia (AML). Combining a multi-omics approach with a functional genetic screen, we reveal that PU.1-downregulated cells fundamentally change their survival control from cytokine-driven pathways to overexpression of an autophagy-predominated stem cell gene program, for which we also find evidence in human AML. Control of this program involves redirected chromatin occupancy of the PU.1 partner TF Runx1 to a lineage-inappropriate binding site repertoire. Hence, genomic reallocation of TF binding upon loss of a partner TF can act as a pro-oncogenic failsafe mechanism by sustaining cell survival during leukemogenesis.

Abstract 4115669: Single-cell and Spatial Transcriptomics Identified Fatty Acid-binding Proteins Controlling Endothelial Glycolytic and Arterial Programming in Pulmonary Hypertension

We are seeking to elucidate the endothelial fatty acid metabolism role for obstructive vascular remodeling in the pathogenesis of PAH. In this study, we bred a severe mouse model of PH Egln1 Tie2Cre mice with Fabp45 -/- mice to generate Egln1 Tie2Cre / Fabp45 -/- mice. We applied single-cell RNA sequencing (scRNA-seq) to profile the pulmonary cells in Egln1 Tie2Cre mice and Egln1 Tie2Cre / Fabp45 -/- mice. Human hPAEC from idiopathic PAH patients and healthy donors were used to measure fatty acid-binding protein 4 and 5 (FABP4 and FABP5) expression via sc-RNA seq and spatial transcriptomics. siRNA-mediated knockdown of FABP4 and FABP5 and lentivirus-mediated FABP4 and 5 overexpression were performed to study cell proliferation, apoptosis, and glycolysis. Our scRNA-seq analysis demonstrated that both FABP4 and 5 were highly induced in the ECs of Egln1 Tie2Cre mice. PAECs from IPAH patients also showed higher expression of FABP4 and 5. Knockdown of FABP4-5 reduced EC proliferation, starvation-induced Caspase 3/7 activity, and glycolysis. Overexpression of FABP4-5 promoted EC glycolysis and proliferation. Genetic deletion of Fabp4 and 5 in Egln1 Tie2Cre mice exhibited a reduction of RVSP, RV hypertrophy, and reduction of EC glycolysis gene programming compared to Egln1 Tie2Cre mice. In conclusion, we found that FABP4 and 5 control EC glycolysis and contribute to the development of PAH.

STEM-14. INTEGRATED ANALYSIS OF GLIOBLASTOMA PROGENITOR SUBTYPES REVEALS A TUMOR-DERIVED STEM CELL POPULATION WITH NEURAL AND VASCULAR POTENTIAL

Abstract Glioblastoma (GBM) exhibits a vast heterogeneity in cell types and states. Recent research highlights the reactivation of neurodevelopmental programs in GBM, contributing to this heterogeneity and, ultimately, therapy resistance. To characterize the landscape of progenitor cell types in human primary GBM, we compiled a meta-atlas of single-cell transcriptomes from seven published studies. Using a novel bioinformatic method for gene program analysis, we derived GBM gene expression meta-modules represented across multiple tumors. Our analyses validated previously described GBM gene programs and revealed novel tumor progenitor subtypes. Among these subtypes, we identified a tumor-derived stem cell population exhibiting a mixed vascular identity and transcriptional programs reminiscent of neural stem cells, which we identified as a putative neurovascular progenitor (NVP). To our knowledge, this population has not been previously described in GBM. We validated the existence of NVPs in situ using immunofluorescence from primary patient GBM samples. We then performed an enrichment for NVPs using fluorescence-activated cell sorting direct from patient tumors, followed by single-cell RNA sequencing. Examination of data from traditional model systems, including gliomaspheres and mouse GBM models, confirmed the preservation of NVP identity across systems. In vivo analysis of NVPs in multiple mouse models of GBM elucidate the influence of NVPs on GBM cell type composition, highlighting their role as multipotent progenitors. Finally, we utilized single-cell lineage barcoding to trace the progeny of patient-derived NVP cells in a novel organoid transplantation model, revealing NVP differentiation into vascular-like and neural-like cells. These data indicate that NVP cells play crucial roles in tumor plasticity and may represent a valuable therapeutic target for future exploration.

EPCO-43. DECONVOLUTION OF GLIOBLASTOMA RNA-SEQ REVEALS IMMUNE SUBPOPULATIONS PREDICTING THERAPEUTIC RESPONSE

Abstract Bulk RNA-seq is now among the standard of care for patients with newly diagnosed glioblastoma (GBM), but whether this data can predict therapeutic response is unclear. In silico deconvolution of bulk RNA-seq data allows for quantification of cell types comprising the tumor and its microenvironment. Our study tests whether cell type proportions predict therapeutic response in patients with GBM and high-grade glioma (HGG) enrolled in trials of investigational agents. 127 patients with GBM/HGG underwent post-resection/biopsy whole exome RNA-seq and reads were post-processed using standard computational pipelines. In silico deconvolution was performed using GBMDeconvoluteR, which generates cellular proportions comprising the RNA-seq reads. Multivariate Cox regression models were used to examine each cell proportion as a predictor of progression-free (PFS) and overall (OS) survival. Statistical significance was ascertained using Bonferroni correction for multiple comparisons. Among 54 patients treated with a PD-L1 inhibitor (nivolumab) and VEGF inhibitor (bevacizumab), the proportions of NK cells (p=0.0024), radial glial cells (p=0.0010), and mural cells (p=0.0025), was associated with increased PFS, while proportion of endothelial cells (p=0.0008) was associated with decreased PFS. Astrocyte (p=0.0022), macrophage (p=0.0015), NK cell (p&amp;lt;0.0001), and mural cell (p=0.0004) proportions were associated with longer OS, while endothelial cell proportion (p=0.0001) was associated with shorter OS. Among 21 patients treated with a tyrosine kinase inhibitor (ibrutinib) combined with radiation and temozolomide (RT/TMZ), cell proportions assessed at the time of diagnosis did not predict PFS or OS. In 52 patients with HGG treated with a JAK inhibitor (ruxolitinib) with RT/TMZ, only mural cell proportion was associated with a longer OS (p=0.0017). Our results highlight how bulk RNA-seq data can be used to predict therapeutic response. In GBM patients treated with immunotherapies, the proportion of NK cells may predict response. Further work will examine which gene programs are predictive of therapeutic response.

EPCO-08. EPIGENETIC THERAPY POTENTIATES TRANSPOSABLE ELEMENT TRANSCRIPTION TO CREATE TUMOR-ENRICHED ANTIGENS IN GLIOBLASTOMA

Abstract Epigenetic treatment has been studied as a means to augment anti-tumor immune cell activity by increasing the expression of antiviral response genes. Previous studies, including our own, have revealed an additional consequence of epigenetic treatment: reactivation of transposable elements (TEs). Normally, most TEs are epigenetically repressed in healthy somatic tissues. However, epigenetic treatment can reactivate TEs as previously unannotated promoters, particularly in rapidly proliferating cells. These reactivated TEs could generate unique proteins that may act as novel antigens for immunotherapy. To study epigenetic treatment as a novel method to augment immunotherapy in cancers with a low number of somatic mutation-derived antigens, we chose glioblastoma as a model. Specifically, we investigated if drug-induced TE expression increases the antigen repertoire upon combinatorial treatment with DNMT inhibitor (Decitabine) and HDAC inhibitor (Panobinostat), agents currently being tested for clinical use in cancers, including CNS tumors. Using patient-derived primary glioblastoma stem cell (GSC) lines and control cell lines (astrocytes and fibroblasts), we first profiled treatment-induced changes in the epigenetic landscape using WGBS-seq and ATAC-seq. Then, using nanoCAGE-seq, both short-read and long-read RNA-seq, we identified treatment-induced, TE-derived transcripts that are predominantly expressed in glioblastoma cells along with upregulation of antiviral gene programs. Finally, using LC-MS/MS, we validated that tumor-enriched reactivated TEs, which mostly originate from the LTR class, encode HLA-I presented antigens in primary GSCs. Notably, epigenetic therapy also reactivated many TEs in proliferating control cell lines, and thus we provide evidence that targeted approaches like CRISPRa might reduce the side effects of activating nonspecific TE loci. Overall, our findings suggest that the therapeutic utilization of epigenetic treatment-induced, TE-derived antigens holds great promise but also provides a cautionary tale regarding the careful selection of tumor-enriched antigen targets.

PATH-29. LONGITUDINAL EPIGENOMIC EVOLUTION OF IDH-WILDTYPE GLIOBLASTOMA

Abstract Glioblastoma, IDH-wildtype remains a molecularly heterogeneous disease that is refractory to current therapies. To investigate how longitudinal epigenomic evolution of glioblastoma drives tumor progression and treatment resistance, we performed comprehensive histopathologic assessment, targeted genomic sequencing, and genome-wide DNA methylation profiling of paired initial and recurrent IDH-wildtype glioblastoma samples from 98 patients, as well as a separate longitudinal cohort of IDH-mutant astrocytomas for comparison. Primary treatment-naïve IDH-wildtype glioblastomas had more globally hypomethylated genomes compared to primary IDH-mutant astrocytomas. At time of recurrence, IDH-mutant astrocytomas uniformly became more hypomethylated similar to the global methylation levels in IDH-wildtype glioblastomas. In contrast, recurrent IDH-wildtype glioblastomas were heterogeneous with subsets becoming more globally hypermethylated (27%), hypomethylated (22%), or remaining stable (51%). Accordingly, a substantial proportion (52%) of recurrent glioblastomas underwent methylation class switching compared to their matched primary tumor, with a predominance of the RTK2 epigenetic subtype at initial surgery and the Mesenchymal epigenetic subtype at recurrence. Among approximately 850,000 interrogated CpG sites, we identified only 110 and 153 sites that consistently became more hypermethylated or hypomethylated between initial and recurrent IDH-wildtype glioblastomas, which regulate genes involved in neuronal morphogenesis, oligodendrocyte differentiation, and myelination. In contrast, we identified 260 and 18,527 CpG sites that consistently became more hypermethylated or hypomethylated between initial and recurrent IDH-mutant astrocytomas, which regulate genes involved in neuronal stem-like differentiation, cell cycle control, and other biologic processes. Finally, we developed a DNA methylation evolution signature incorporating change in mean beta-methylation values from initial to recurrent tumor specimens across 347 critical CpG sites that significantly correlated with clinical outcomes for patients with IDH-wildtype glioblastoma. In summary, IDH-wildtype glioblastomas demonstrate heterogeneous epigenomic evolution patterns distinct from IDH-mutant astrocytomas. These changes in epigenetic patterns may identify unique subgroups of glioblastoma with differential treatment response and inform future therapeutic studies targeting differentially methylated gene programs.

Zymogen granule protein 16B (ZG16B) is a druggable epigenetic target to modulate the mammary extracellular matrix.

High tissue density of the mammary gland is considered a pro-tumorigenic factor, hence suppressing the stimuli that induce matrix buildup carries the potential for cancer interception. We found that in non-malignant mammary epithelial cells the combination of the chemopreventive agents bexarotene (Bex) and carvedilol (Carv) suppresses the zymogen granule protein 16B (ZG16B, PAUF) through an interaction of ARID1A with a proximal enhancer. Bex + Carv also reduced ZG16B levels invivo in normal breast tissue and MDA-MB231 tumor xenografts. The relevance of ZG16B is underscored by ongoing clinical trials targeting ZG16B in pancreatic cancers, but its role in breast cancer development is unclear. In immortalized mammary epithelial cells, secreted recombinant ZG16B stimulated mitogenic kinase phosphorylation, detachment and mesenchymal characteristics, and promoted proliferation, motility and clonogenic growth. Highly concerted induction of specific laminin, collagen and integrin isoforms indicated a shift in matrix properties toward increased density and cell-matrix interactions. Exogenous ZG16B alone blocked Bex + Carv-mediated control of cell growth and migration, and antagonized Bex + Carv-induced gene programs regulating cell adhesion and migration. In breast cancer cells ZG16B induced colony formation and anchorage-independent growth, and stimulated migration in a PI3K/Akt-dependent manner. In contrast, Bex + Carv inhibited colony formation, reduced Ki67 levels, ZG16B expression and glucose uptake in MDA-MB231 xenografts. These data establish ZG16B as a druggable pro-tumorigenic target in breast cell transformation and suggest a key role of the matrisome network in rexinoid-dependent antitumor activity.

Induced hepatocyte-like cells derived from adipose-derived stem cells alleviates liver injury in mice infected with Echinococcus Multilocularis

Accumulating evidence has shown that adipose-derived stem cells (ADSCs) have the potential to differentiate into hepatic lineages, which are ideal engraftments for tissue-engineered repair. In this study, we investigated the potential of transplanted induced hepatocyte-like cells (iHEPs) in treating hepatic alveolar echinococcosis and describe an efficient three-step protocol for the generation of iHEPs in vitro from ADSCs. The expression of hepatocyte lineage markers was assessed and iHEPs function was evaluated by Periodic acid-Schiff staining. iHEPs were intravenously transplanted into mice infected with Echinococcus multilocularis. Histopathological analysis and liver function tests were used to assess therapeutic effects. The iHEPs exhibit morphological features and a glycogen storage function similar to those of mature hepatocytes and demonstrate an upregulation in hepatic gene programs with increasing induction time. Following transplantation, iHEPs were observed surrounding the metacestode lesions in the liver parenchyma of E. multilocularis-infected mice. iHEPs transplantation effectively restored liver function and improved liver injury in the infected mice. Additionally, we observed significant activation of the Wnt/β-catenin signaling pathway in the livers of infected mice transplanted with iHEPs. Our results provide evidence that iHEPs transplantation can alleviate E. multilocularis-induced liver injury, potentially creating new avenues for treating liver injury in end-stage hepatic alveolar echinococcosis.

Contractile regulation of the glucocorticoid-sensitive transcriptome in young and aged skeletal muscle.

Elevated glucocorticoids alter the skeletal muscle transcriptome to induce a myopathy characterized by muscle atrophy, muscle weakness, and decreased metabolic function. These effects are more likely to occur and be more severe in aged muscles. Resistance exercise can blunt the development of glucocorticoid myopathy in young muscle, but the potential to oppose the signals initiating myopathy in aged muscle is unknown. To answer this, young (4-mo-old) and aged (24-to 25-mo-old) male C57BL/6 mice were randomized to receive either an intraperitoneal (IP) injection of dexamethasone (DEX; 2 mg/kg) or saline as a control. Two hours postinjections, the tibialis anterior (TA) muscles of mice were subjected to unilateral high-force contractions. Muscles were harvested 4 h later. The glucocorticoid- and contraction-sensitive genes were determined by RNA sequencing. The number of glucocorticoid-sensitive genes was similar between young and aged muscle. Contractions opposed changes to more glucocorticoid-sensitive genes in aged muscle, with this outcome primarily occurring when hormone levels were elevated. Glucocorticoid-sensitive gene programs opposed by contractions were primarily related to metabolism in young mice and muscle size regulation and inflammation in aged mice. In silico analysis implied peroxisome proliferator-activated receptor gamma-1 (PPARG) contributed to the contraction-induced opposition of glucocorticoid-sensitive genes in aged muscle. Increasing PPARG expression in the TA of aged mice using adeno-associated virus serotype 9 partially counteracted the glucocorticoid-induced reduction in runt-related transcription factor 1 (Runx1) mRNA content, recapitulating the effects observed by contractions. Overall, these data contribute to our understanding of the contractile regulation of the glucocorticoid transcriptome in aged skeletal muscle.NEW & NOTEWORTHY We establish the extent to which muscle contractions oppose changes to the glucocorticoid-sensitive transcriptome in both young and aged muscle. We also identify peroxisome proliferator-activated receptor gamma (PPARG) as a transcription factor likely contributing to contraction-induced opposition to the glucocorticoid transcriptome in aged muscle. Overall, these data contribute to our understanding of the contractile regulation of the glucocorticoid transcriptome.