MAPK/ERK signaling in gliomas modulates interferon responses, T cell recruitment, microglia phenotype, and immune checkpoint blockade efficacy.

Background: Glioblastoma (GB) remains a formidable challenge in neuro-oncology, with immune checkpoint blockade (ICB) showing limited efficacy in unselected patients. We previously recently established that MAPK/ERK signaling is associated with overall survival following anti-PD-1 and anti-CTLA-4 treatment in recurrent GB. However, the causal relationship between MAPK/ERK signaling and susceptibility to ICB, as well as the mechanisms underlying this association, remain poorly understood. Method: We conducted in vivo kinome-wide CRISPR/Cas9 screenings in murine gliomas to identify key regulators of susceptibility to anti-PD-1 and CD8+ T cell responses and performed survival studies to validate the most relevant genes. Additionally, paired single cell RNA-sequencing (scRNA-seq) with p-ERK staining, spatial transcriptomics on GB samples, and ex-vivo slice culture of a BRAFV600E mutant GB tumor treated with BRAFi/MEKi were used to determine the causal relationship between MAPK signaling, tumor cell immunogenicity, and modulation of microglia phenotype. Results: CRISPR/Cas9 screens identified the MAPK pathway, particularly the RAF-MEK-ERK pathway, as the most critical modulator of glioma susceptibility to CD8+ T cells, and anti-PD-1 across all kinases. Experimentally-induced ERK phosphorylation in gliomas enhanced survival with ICB treatment, led to durable anti-tumoral immunity upon re-challenge and memory T cell infiltration in long-term survivors. Elevated p-ERK in glioma cells correlated with increased interferon responses, antigen presentation and T cell infiltration in GB. Moreover, spatial transcriptomics and scRNA-seq analysis revealed the modulation of interferon responses by the MAPK/ERK pathway in BRAFV600E human GB cells with ERK1/2 knockout and in slice cultures of human BRAFV600E GB tissue. Notably, BRAFi/MEKi treatment disrupted the interaction between tumor cells and tumor-associated macrophages/microglia in slice cultures from BRAFV600E mutant GB. Conclusion: Our data indicate that the MAPK/ERK pathway is a critical regulator of GB cell susceptibility to anti-tumoral immunity, modulating interferon responses, and antigen-presentation in glioma cells, as well as tumor cell interaction with microglia. These findings not only elucidate the mechanistic underpinnings of immunotherapy resistance in GB but also highlight the MAPK/ERK pathway as a promising target for enhancing immunotherapeutic efficacy in this challenging malignancy.

Tumor detection by analysis of both symmetric- and hemi-methylation of plasma cell-free DNA.

Aberrant DNA methylation patterns have been used for cancer detection. However, DNA hemi-methylation, present at about 10% CpG dinucleotides, has been less well studied. Here we show that a majority of differentially hemi-methylated regions (DHMRs) in liver tumor DNA or plasma cells free (cf) DNA do not overlap with differentially methylated regions (DMRs) of the same samples, indicating that DHMRs could serve as independent biomarkers. Furthermore, we analyzed the cfDNA methylomes of 215 samples from individuals with liver or brain cancer and individuals without cancer (controls), and trained machine learning models using DMRs, DHMRs or both. The models incorporated with both DMRs and DHMRs show a superior performance compared to models trained with DMRs or DHMRs, with AUROC being 0.978, 0.990, and 0.983 in distinguishing control, liver and brain cancer, respectively, in a validation cohort. This study supports the potential of utilizing both DMRs and DHMRs for multi-cancer detection.

Fc-enhanced anti-CTLA-4, anti-PD-1, doxorubicin, and ultrasound-mediated BBB opening: A novel combinatorial immunotherapy regimen for gliomas.

BACKGROUND: Glioblastoma is a highly aggressive brain cancer that is resistant to conventional immunotherapy strategies. Botensilimab, an Fc-enhanced anti-CTLA-4 antibody (FcE-aCTLA-4), has shown durable activity in "cold" and immunotherapy-refractory cancers. METHOD: We evaluated the efficacy and immune microenvironment phenotype of a mouse analogue of FcE-aCTLA-4 in treatment-refractory preclinical models of glioblastoma, both as a monotherapy and in combination with doxorubicin delivered via low-intensity pulsed ultrasound and microbubbles (LIPU/MB). Additionally, we studied 4 glioblastoma patients treated with doxorubicin, anti-PD-1 with concomitant LIPU/MB to investigate the novel effect of doxorubicin modulating FcγR expressions in tumor associated macrophages/microglia (TAMs). RESULTS: FcE-aCTLA-4 demonstrated high-affinity binding to FcγRIV, the mouse ortholog of human FcγRIIIA, which was highly expressed in TAMs in human glioblastoma, most robustly at diagnosis. Notably, FcE-aCTLA-4 mediated selective depletion of intra-tumoral regulatory T cells (Tregs) via TAM-mediated phagocytosis, while sparing peripheral Tregs. Doxorubicin, a chemotherapeutic drug with immunomodulatory functions, was found to upregulate FcγRIIIA on TAMs in glioblastoma patients who received doxorubicin and anti-PD-1 with concomitant LIPU/MB. In murine models of immunotherapy-resistant gliomas, a combinatorial regimen of FcE-aCTLA-4, anti-PD-1, and doxorubicin with LIPU/MB, achieved a 90% cure rate, that was associated robust infiltration of activated CD8+ T cells and establishment of immunological memory as evidenced by rejection upon tumor rechallenge. CONCLUSION: Our findings demonstrate that FcE-aCTLA-4 promotes robust immunomodulatory and anti-tumor effects in murine gliomas and is significantly enhanced when combined with anti-PD-1, doxorubicin, and LIPU/MB. We are currently investigating this combinatory strategy in a clinical trial (clinicaltrials.gov NCT05864534).

GET: a foundation model of transcription across human cell types.

Transcriptional regulation, involving the complex interplay between regulatory sequences and proteins, directs all biological processes. Computational models of transcription lack generalizability to accurately extrapolate in unseen cell types and conditions. Here, we introduce GET, an interpretable foundation model designed to uncover regulatory grammars across 213 human fetal and adult cell types. Relying exclusively on chromatin accessibility data and sequence information, GET achieves experimental-level accuracy in predicting gene expression even in previously unseen cell types. GET showcases remarkable adaptability across new sequencing platforms and assays, enabling regulatory inference across a broad range of cell types and conditions, and uncovering universal and cell type specific transcription factor interaction networks. We evaluated its performance on prediction of regulatory activity, inference of regulatory elements and regulators, and identification of physical interactions between transcription factors. Specifically, we show GET outperforms current models in predicting lentivirus-based massive parallel reporter assay readout with reduced input data. In fetal erythroblasts, we identify distal (>1Mbp) regulatory regions that were missed by previous models. In B cells, we identified a lymphocyte-specific transcription factor-transcription factor interaction that explains the functional significance of a leukemia-risk predisposing germline mutation. In sum, we provide a generalizable and accurate model for transcription together with catalogs of gene regulation and transcription factor interactions, all with cell type specificity.

Understanding variants of unknown significance: the computational frontier.

The rapid advancement of sequencing technologies has led to the identification of numerous mutations in cancer genomes, many of which are variants of unknown significance (VUS). Computational models are increasingly being used to predict the functional impact of these mutations, in both coding and noncoding regions. Integration of these models with emerging genomic datasets will refine our understanding of mutation effects and guide clinical decision making. Future advancements in modeling protein interactions and transcriptional regulation will further enhance our ability to interpret VUS. Periodic incorporation of these developments into VUS reclassification practice has the potential to significantly improve personalized cancer care.

Systematic Characterization of p53-Regulated Long Noncoding RNAs across Human Cancers Reveals Remarkable Heterogeneity among Different Tumor Types.

The p53 tumor suppressor protein, a sequence-specific DNA binding transcription factor, regulates the expression of a large number of genes, in response to various forms of cellular stress. Although the protein coding target genes of p53 have been well studied, less is known about its role in regulating long noncoding genes and their functional relevance to cancer. Here we report the genome-wide identification of a large set (>1,000) of long noncoding RNAs (lncRNA), which are putative p53 targets in a colon cancer cell line and in human patient datasets from five different common types of cancer. These lncRNAs have not been annotated by other studies of normal unstressed systems. In the colon cancer cell line, a high proportion of these lncRNAs are uniquely induced by different chemotherapeutic agents that activate p53, whereas others are induced by more than one agent tested. Further, subsets of these lncRNAs independently predict overall and disease-free survival of patients across the five different common cancer types. Interestingly, both genetic alterations and patient survival associated with different lncRNAs are unique to each cancer tested, indicating extraordinary tissue-specific variability in the p53 noncoding response. The newly identified noncoding p53 target genes have allowed us to construct a classifier for tumor diagnosis and prognosis. IMPLICATIONS: Our results not only identify myriad p53-regulated long noncoding (lncRNA), they also reveal marked drug-induced, as well as tissue- and tumor-specific heterogeneity in these putative p53 targets and our findings have enabled the construction of robust classifiers for diagnosis and prognosis.

Immune Response following FLASH and Conventional Radiation in Diffuse Midline Glioma.

PURPOSE: Diffuse midline glioma (DMG) is a fatal tumor traditionally treated with radiation therapy (RT) and previously characterized as having a noninflammatory tumor immune microenvironment (TIME). FLASH is a novel RT technique using ultra-high dose rate that is associated with decreased toxicity and effective tumor control. However, the effect of FLASH and conventional (CONV) RT on the DMG TIME has not yet been explored. METHODS AND MATERIALS: Here, we performed single-cell RNA sequencing (scRNA-seq) and flow cytometry on immune cells isolated from an orthotopic syngeneic murine model of brainstem DMG after the use of FLASH (90 Gy/sec) or CONV (2 Gy/min) dose-rate RT and compared to unirradiated tumor (SHAM). RESULTS: At day 4 post-RT, FLASH exerted similar effects as CONV in the predominant microglial (MG) population, including the presence of two activated subtypes. However, at day 10 post-RT, we observed a significant increase in the type 1 interferon α/ß receptor (IFNAR+) in MG in CONV and SHAM compared to FLASH. In the non-resident myeloid clusters of macrophages (MACs) and dendritic cells (DCs), we found increased type 1 interferon (IFN1) pathway enrichment for CONV compared to FLASH and SHAM by scRNA-seq. We observed this trend by flow cytometry at day 4 post-RT in IFNAR+ MACs and DCs, which equalized by day 10 post-RT. DMG control and murine survival were equivalent between RT dose rates. CONCLUSIONS: Our work is the first to map CONV and FLASH immune alterations of the DMG TIME with single-cell resolution. Although DMG tumor control and survival were similar between CONV and FLASH, we found that changes in immune compartments differed over time. Importantly, although both RT modalities increased IFN1, we found that the timing of this response was cell-type and dose-rate dependent. These temporal differences, particularly in the context of tumor control, warrant further study.

Meta-analyses of mouse and human prostate single-cell transcriptomes reveal widespread epithelial plasticity in tissue regression, regeneration, and cancer.

Recent advances in single-cell RNA-sequencing (scRNA-seq) technology have facilitated studies of cell states and plasticity in tissue maintenance and cancer, including in the prostate. Here we present meta-analyses of multiple new and published scRNA-seq datasets to establish reference cell type classifications for the normal mouse and human prostate. Our analyses demonstrate transcriptomic similarities between epithelial cell states in the normal prostate, in the regressed prostate after androgen-deprivation, and in primary prostate tumors. During regression in the mouse prostate, all epithelial cells shift their expression profiles towards a proximal periurethral (PrU) state, demonstrating an androgen-dependent plasticity that is restored to normal during androgen restoration and regeneration. In the human prostate, we find progressive rewiring of transcriptional programs across epithelial cell types in benign prostate hyperplasia and treatment-naïve prostate cancer. Notably, we detect copy number variants predominantly within Luminal Acinar cells in prostate tumors, suggesting a bias in their cell type of origin, as well as a larger field of transcriptomic alterations in non-tumor cells. Finally, we observe that Luminal Acinar tumor cells in treatment-naïve prostate cancer display heterogeneous androgen receptor (AR) signaling activity, including a split between high-AR and low-AR profiles with similarity to PrU-like states. Taken together, our analyses of cellular heterogeneity and plasticity provide important translational insights into the origin and treatment response of prostate cancer.

Virome Data Explorer: A web resource to longitudinally explore respiratory viral infections, their interactions with other pathogens and host transcriptomic changes in over 100 people.

Viral respiratory infections are an important public health concern due to their prevalence, transmissibility, and potential to cause serious disease. Disease severity is the product of several factors beyond the presence of the infectious agent, including specific host immune responses, host genetic makeup, and bacterial coinfections. To understand these interactions within natural infections, we designed a longitudinal cohort study actively surveilling respiratory viruses over the course of 19 months (2016 to 2018) in a diverse cohort in New York City. We integrated the molecular characterization of 800+ nasopharyngeal samples with clinical data from 104 participants. Transcriptomic data enabled the identification of respiratory pathogens in nasopharyngeal samples, the characterization of markers of immune response, the identification of signatures associated with symptom severity, individual viruses, and bacterial coinfections. Specific results include a rapid restoration of baseline conditions after infection, significant transcriptomic differences between symptomatic and asymptomatic infections, and qualitatively similar responses across different viruses. We created an interactive computational resource (Virome Data Explorer) to facilitate access to the data and visualization of analytical results.

Smoother: a unified and modular framework for incorporating structural dependency in spatial omics data.

Spatial omics technologies can help identify spatially organized biological processes, but existing computational approaches often overlook structural dependencies in the data. Here, we introduce Smoother, a unified framework that integrates positional information into non-spatial models via modular priors and losses. In simulated and real datasets, Smoother enables accurate data imputation, cell-type deconvolution, and dimensionality reduction with remarkable efficiency. In colorectal cancer, Smoother-guided deconvolution reveals plasma cell and fibroblast subtype localizations linked to tumor microenvironment restructuring. Additionally, joint modeling of spatial and single-cell human prostate data with Smoother allows for spatial mapping of reference populations with significantly reduced ambiguity.

A pan-cancer analysis implicates human NKIRAS1 as a tumor-suppressor gene.

The NF-κB family of transcription factors and the Ras family of small GTPases are important mediators of proproliferative signaling that drives tumorigenesis and carcinogenesis. The κB-Ras proteins were previously shown to inhibit both NF-κB and Ras activation through independent mechanisms, implicating them as tumor suppressors with potentially broad relevance to human cancers. In this study, we have used two mouse models to establish the relevance of the κB-Ras proteins for tumorigenesis. Additionally, we have utilized a pan-cancer bioinformatics analysis to explore the role of the κB-Ras proteins in human cancers. Surprisingly, we find that the genes encoding κB-Ras 1 (NKIRAS1) and κB-Ras 2 (NKIRAS2) are rarely down-regulated in tumor samples with oncogenic Ras mutations. Reduced expression of human NKIRAS1 alone is associated with worse prognosis in at least four cancer types and linked to a network of genes implicated in tumorigenesis. Our findings provide direct evidence that loss of NKIRAS1 in human tumors that do not carry oncogenic RAS mutations is associated with worse clinical outcomes.

Aberrant pace of cortical neuron development in brain organoids from patients with 22q11.2 deletion syndrome and schizophrenia.

Adults and children afflicted with the 22q11.2 deletion syndrome (22q11.2DS) exhibit cognitive, social, and emotional impairments, and are at significantly heightened risk for schizophrenia (SCZ). The impact of this deletion on early human brain development, however, has remained unclear. Here we harness organoid models of the developing human cerebral cortex, cultivated from subjects with 22q11.2DS and SCZ, as well as unaffected control samples, to identify cell-type-specific developmental abnormalities arising from this genomic lesion. Leveraging single-cell RNA-sequencing in conjunction with experimental validation, we find that the loss of genes within the 22q11.2 locus leads to a delayed development of cortical neurons. This compromised development was reflected in an elevated proportion of actively proliferating neural progenitor cells, coupled with a decreased fraction of more mature neurons. Furthermore, we identify perturbed molecular imprints linked to neuronal maturation, observe the presence of sparser neurites, and note a blunted amplitude in glutamate-induced Ca2+ transients. The aberrant transcription program underlying impaired development contains molecular signatures significantly enriched in neuropsychiatric genetic liability. MicroRNA profiling and target gene investigation suggest that microRNA dysregulation may drive perturbations of genes governing the pace at which maturation unfolds. Using protein-protein interaction network analysis we define complementary effects stemming from additional genes residing within the deleted locus. Our study uncovers reproducible neurodevelopmental and molecular alterations due to 22q11.2 deletions. These findings have the potential to facilitate disease modeling and promote the pursuit of therapeutic interventions.

Radiation therapy promotes unsaturated fatty acids to maintain survival of glioblastoma.

Radiation therapy (RT) is essential for the management of glioblastoma (GBM). However, GBM frequently relapses within the irradiated margins, thus suggesting that RT might stimulate mechanisms of resistance that limits its efficacy. GBM is recognized for its metabolic plasticity, but whether RT-induced resistance relies on metabolic adaptation remains unclear. Here, we show in vitro and in vivo that irradiated GBM tumors switch their metabolic program to accumulate lipids, especially unsaturated fatty acids. This resulted in an increased formation of lipid droplets to prevent endoplasmic reticulum (ER) stress. The reduction of lipid accumulation with genetic suppression and pharmacological inhibition of the fatty acid synthase (FASN), one of the main lipogenic enzymes, leads to mitochondrial dysfunction and increased apoptosis of irradiated GBM cells. Combination of FASN inhibition with focal RT improved the median survival of GBM-bearing mice. Supporting the translational value of these findings, retrospective analysis of the GLASS consortium dataset of matched GBM patients revealed an enrichment in lipid metabolism signature in recurrent GBM compared to primary. Overall, these results demonstrate that RT drives GBM resistance by generating a lipogenic environment permissive to GBM survival. Targeting lipid metabolism might be required to develop more effective anti-GBM strategies.

Mechanisms of insertions at a DNA double-strand break.

Insertions and deletions (indels) are common sources of structural variation, and insertions originating from spontaneous DNA lesions are frequent in cancer. We developed a highly sensitive assay called insertion and deletion sequencing (Indel-seq) to monitor rearrangements in human cells at the TRIM37 acceptor locus that reports indels stemming from experimentally induced and spontaneous genome instability. Templated insertions, which derive from sequences genome wide, require contact between donor and acceptor loci, require homologous recombination, and are stimulated by DNA end-processing. Insertions are facilitated by transcription and involve a DNA/RNA hybrid intermediate. Indel-seq reveals that insertions are generated via multiple pathways. The broken acceptor site anneals with a resected DNA break or invades the displaced strand of a transcription bubble or R-loop, followed by DNA synthesis, displacement, and then ligation by non-homologous end joining. Our studies identify transcription-coupled insertions as a critical source of spontaneous genome instability that is distinct from cut-and-paste events.

ATR kinase supports normal proliferation in the early S phase by preventing replication resource exhaustion.

The ATR kinase, which coordinates cellular responses to DNA replication stress, is also essential for the proliferation of normal unstressed cells. Although its role in the replication stress response is well defined, the mechanisms by which ATR supports normal cell proliferation remain elusive. Here, we show that ATR is dispensable for the viability of G0-arrested naïve B cells. However, upon cytokine-induced proliferation, Atr-deficient B cells initiate DNA replication efficiently, but by mid-S phase they display dNTP depletion, fork stalling, and replication failure. Nonetheless, productive DNA replication and dNTP levels can be restored in Atr-deficient cells by suppressing origin firing, such as partial inhibition of CDC7 and CDK1 kinase activities. Together, these findings indicate that ATR supports the proliferation of normal unstressed cells by tempering the pace of origin firing during the early S phase to avoid exhaustion of dNTPs and importantly also other replication factors.

ATR kinase supports normal proliferation in the early S phase by preventing replication resource exhaustion.

The ATR kinase, which coordinates cellular responses to DNA replication stress, is also essential for the proliferation of normal unstressed cells. Although its role in the replication stress response is well defined, the mechanisms by which ATR supports normal cell proliferation remain elusive. Here, we show that ATR is dispensable for the viability of G0-arrested naïve B cells. However, upon cytokine-induced proliferation, Atr-deficient B cells initiate DNA replication efficiently in early S phase, but by mid-S phase they display dNTP depletion, fork stalling, and replication failure. Nonetheless, productive DNA replication can be restored in Atr-deficient cells by pathways that suppress origin firing, such as downregulation of CDC7 and CDK1 kinase activities. Together, these findings indicate that ATR supports the proliferation of normal unstressed cells by tempering the pace of origin firing during the early S phase to avoid exhaustion of dNTPs and other replication factors.

A transposase-derived gene required for human brain development.

DNA transposable elements and transposase-derived genes are present in most living organisms, including vertebrates, but their function is largely unknown. PiggyBac Transposable Element Derived 5 (PGBD5) is an evolutionarily conserved vertebrate DNA transposase-derived gene with retained nuclease activity in cells. Vertebrate brain development is known to be associated with prominent neuronal cell death and DNA breaks, but their causes and functions are not well understood. Here, we show that PGBD5 contributes to normal brain development in mice and humans, where its deficiency causes disorder of intellectual disability, movement and seizures. In mice, Pgbd5 is required for the developmental induction of post-mitotic DNA breaks and recurrent somatic genome rearrangements in neurons. Together, these studies nominate PGBD5 as the long-hypothesized neuronal DNA nuclease required for brain function in mammals.

Inverse agonists of retinoic acid receptor/retinoid X receptor signaling as lineage-specific antitumor agents against human adenoid cystic carcinoma.

BACKGROUND: Adenoid cystic carcinoma (ACC) is a lethal malignancy of exocrine glands, characterized by the coexistence within tumor tissues of 2 distinct populations of cancer cells, phenotypically similar to the myoepithelial and ductal lineages of normal salivary epithelia. The developmental relationship linking these 2 cell types, and their differential vulnerability to antitumor treatments, remains unknown. METHODS: Using single-cell RNA sequencing, we identified cell-surface markers (CD49f, KIT) that enabled the differential purification of myoepithelial-like (CD49fhigh/KITneg) and ductal-like (CD49flow/KIT+) cells from patient-derived xenografts (PDXs) of human ACCs. Using prospective xenotransplantation experiments, we compared the tumor-initiating capacity of the 2 cell types and tested whether one could differentiate into the other. Finally, we searched for signaling pathways with differential activation between the 2 cell types and tested their role as lineage-specific therapeutic targets. RESULTS: Myoepithelial-like cells displayed higher tumorigenicity than ductal-like cells and acted as their progenitors. Myoepithelial-like and ductal-like cells displayed differential expression of genes encoding for suppressors and activators of retinoic acid signaling, respectively. Agonists of retinoic acid receptor (RAR) or retinoid X receptor (RXR) signaling (all-trans retinoic acid, bexarotene) promoted myoepithelial-to-ductal differentiation, whereas suppression of RAR/RXR signaling with a dominant-negative RAR construct abrogated it. Inverse agonists of RAR/RXR signaling (BMS493, AGN193109) displayed selective toxicity against ductal-like cells and in vivo antitumor activity against PDX models of human ACC. CONCLUSIONS: In human ACCs, myoepithelial-like cells act as progenitors of ductal-like cells, and myoepithelial-to-ductal differentiation is promoted by RAR/RXR signaling. Suppression of RAR/RXR signaling is lethal to ductal-like cells and represents a new therapeutic approach against human ACCs.

Checkpoint kinase 1/2 inhibition potentiates anti-tumoral immune response and sensitizes gliomas to immune checkpoint blockade.

Whereas the contribution of tumor microenvironment to the profound immune suppression of glioblastoma (GBM) is clear, tumor-cell intrinsic mechanisms that regulate resistance to CD8 T cell mediated killing are less understood. Kinases are potentially druggable targets that drive tumor progression and might influence immune response. Here, we perform an in vivo CRISPR screen to identify glioma intrinsic kinases that contribute to evasion of tumor cells from CD8 T cell recognition. The screen reveals checkpoint kinase 2 (Chek2) to be the most important kinase contributing to escape from CD8 T-cell recognition. Genetic depletion or pharmacological inhibition of Chek2 with blood-brain-barrier permeable drugs that are currently being evaluated in clinical trials, in combination with PD-1 or PD-L1 blockade, lead to survival benefit in multiple preclinical glioma models. Mechanistically, loss of Chek2 enhances antigen presentation, STING pathway activation and PD-L1 expression in mouse gliomas. Analysis of human GBMs demonstrates that Chek2 expression is inversely associated with antigen presentation and T-cell activation. Collectively, these results support Chek2 as a promising target for enhancement of response to immune checkpoint blockade therapy in GBM.