Microglia and monocyte-derived macrophages drive progression of pediatric high-grade gliomas and are transcriptionally shaped by histone mutations.

Pediatric high-grade gliomas (pHGGs), including hemispheric pHGGs and diffuse midline gliomas (DMGs), harbor mutually exclusive tumor location-specific histone mutations. Using immunocompetent de novo mouse models of pHGGs, we demonstrated that myeloid cells were the predominant infiltrating non-neoplastic cell population. Single-cell RNA sequencing (scRNA-seq), flow cytometry, and immunohistochemistry illustrated the presence of heterogeneous myeloid cell populations shaped by histone mutations and tumor location. Disease-associated myeloid (DAM) cell phenotypes demonstrating immune permissive characteristics were identified in murine and human pHGG samples. H3.3K27M DMGs, the most aggressive DMG, demonstrated enrichment of DAMs. Genetic ablation of chemokines Ccl8 and Ccl12 resulted in a reduction of DAMs and an increase in lymphocyte infiltration, leading to increased survival of tumor-bearing mice. Pharmacologic inhibition of chemokine receptors CCR1 and CCR5 resulted in extended survival and decreased myeloid cell infiltration. This work establishes the tumor-promoting role of myeloid cells in DMG and the potential therapeutic opportunities for targeting them.

Hypoxic niches attract and sequester tumor-associated macrophages and cytotoxic T cells and reprogram them for immunosuppression.

Glioblastoma (GBM), a highly lethal brain cancer, is notorious for immunosuppression, but the mechanisms remain unclear. Here, we documented a temporospatial patterning of tumor-associated myeloid cells (TAMs) corresponding to vascular changes during GBM progression. As tumor vessels transitioned from the initial dense regular network to later scant and engorged vasculature, TAMs shifted away from perivascular regions and trafficked to vascular-poor areas. This process was heavily influenced by the immunocompetence state of the host. Utilizing a sensitive fluorescent UnaG reporter to track tumor hypoxia, coupled with single-cell transcriptomics, we revealed that hypoxic niches attracted and sequestered TAMs and cytotoxic T lymphocytes (CTLs), where they were reprogrammed toward an immunosuppressive state. Mechanistically, we identified chemokine CCL8 and cytokine IL-1ß as two hypoxic-niche factors critical for TAM trafficking and co-evolution of hypoxic zones into pseudopalisading patterns. Therefore, perturbation of TAM patterning in hypoxic zones may improve tumor control.

A revised airway epithelial hierarchy includes CFTR-expressing ionocytes.

The airways of the lung are the primary sites of disease in asthma and cystic fibrosis. Here we study the cellular composition and hierarchy of the mouse tracheal epithelium by single-cell RNA-sequencing (scRNA-seq) and in vivo lineage tracing. We identify a rare cell type, the Foxi1+ pulmonary ionocyte; functional variations in club cells based on their location; a distinct cell type in high turnover squamous epithelial structures that we term 'hillocks'; and disease-relevant subsets of tuft and goblet cells. We developed 'pulse-seq', combining scRNA-seq and lineage tracing, to show that tuft, neuroendocrine and ionocyte cells are continually and directly replenished by basal progenitor cells. Ionocytes are the major source of transcripts of the cystic fibrosis transmembrane conductance regulator in both mouse (Cftr) and human (CFTR). Knockout of Foxi1 in mouse ionocytes causes loss of Cftr expression and disrupts airway fluid and mucus physiology, phenotypes that are characteristic of cystic fibrosis. By associating cell-type-specific expression programs with key disease genes, we establish a new cellular narrative for airways disease.

Novel microenvironment-based classification of intrahepatic cholangiocarcinoma with therapeutic implications.

OBJECTIVE: The diversity of the tumour microenvironment (TME) of intrahepatic cholangiocarcinoma (iCCA) has not been comprehensively assessed. We aimed to generate a novel molecular iCCA classifier that incorporates elements of the stroma, tumour and immune microenvironment ('STIM' classification). DESIGN: We applied virtual deconvolution to transcriptomic data from ~900 iCCAs, enabling us to devise a novel classification by selecting for the most relevant TME components. Murine models were generated through hydrodynamic tail vein injection and compared with the human disease. RESULTS: iCCA is composed of five robust STIM classes encompassing both inflamed (35%) and non-inflamed profiles (65%). The inflamed classes, named immune classical (~10%) and inflammatory stroma (~25%), differ in oncogenic pathways and extent of desmoplasia, with the inflammatory stroma showing T cell exhaustion, abundant stroma and KRAS mutations (p<0.001). Analysis of cell-cell interactions highlights cancer-associated fibroblast subtypes as potential mediators of immune evasion. Among the non-inflamed classes, the desert-like class (~20%) harbours the lowest immune infiltration with abundant regulatory T cells (p<0.001), whereas the hepatic stem-like class (~35%) is enriched in 'M2-like' macrophages, mutations in IDH1/2 and BAP1, and FGFR2 fusions. The remaining class (tumour classical: ~10%) is defined by cell cycle pathways and poor prognosis. Comparative analysis unveils high similarity between a KRAS/p19 murine model and the inflammatory stroma class (p=0.02). The KRAS-SOS inhibitor, BI3406, sensitises a KRAS-mutant iCCA murine model to anti-PD1 therapy. CONCLUSIONS: We describe a comprehensive TME-based stratification of iCCA. Cross-species analysis establishes murine models that align closely to human iCCA for the preclinical testing of combination strategies.

Loss of DNA methyltransferase activity in primed human ES cells triggers increased cell-cell variability and transcriptional repression.

Maintenance of pluripotency and specification towards a new cell fate are both dependent on precise interactions between extrinsic signals and transcriptional and epigenetic regulators. Directed methylation of cytosines by the de novo methyltransferases DNMT3A and DNMT3B plays an important role in facilitating proper differentiation, whereas DNMT1 is essential for maintaining global methylation levels in all cell types. Here, we generated single-cell mRNA expression data from wild-type, DNMT3A, DNMT3A/3B and DNMT1 knockout human embryonic stem cells and observed a widespread increase in cellular and transcriptional variability, even with limited changes in global methylation levels in the de novo knockouts. Furthermore, we found unexpected transcriptional repression upon either loss of the de novo methyltransferase DNMT3A or the double knockout of DNMT3A/3B that is further propagated upon differentiation to mesoderm and ectoderm. Taken together, our single-cell RNA-sequencing data provide a high-resolution view into the consequences of depleting the three catalytically active DNMTs in human pluripotent stem cells.

Transcription factor binding dynamics during human ES cell differentiation.

Pluripotent stem cells provide a powerful system to dissect the underlying molecular dynamics that regulate cell fate changes during mammalian development. Here we report the integrative analysis of genome-wide binding data for 38 transcription factors with extensive epigenome and transcriptional data across the differentiation of human embryonic stem cells to the three germ layers. We describe core regulatory dynamics and show the lineage-specific behaviour of selected factors. In addition to the orchestrated remodelling of the chromatin landscape, we find that the binding of several transcription factors is strongly associated with specific loss of DNA methylation in one germ layer, and in many cases a reciprocal gain in the other layers. Taken together, our work shows context-dependent rewiring of transcription factor binding, downstream signalling effectors, and the epigenome during human embryonic stem cell differentiation.

Toward a unified physical model of nucleosome patterns flanking transcription start sites.

Recent genome-wide maps of nucleosome positions in different eukaryotes revealed patterns around transcription start sites featuring a nucleosome-free region flanked by a periodic modulation of the nucleosome density. For Saccharomyces cerevisiae, the average in vivo pattern was previously shown to be quantitatively described by a "nucleosome gas" model based on the statistical positioning mechanism. However, this simple physical description is challenged by the fact that the pattern differs quantitatively between species and by recent experiments that appear incompatible with statistical positioning, indicating important roles for chromatin remodelers. We undertake a data-driven search for a unified physical model to describe the nucleosome patterns of 12 yeast species and also consider an extension of the model to capture remodeling effects. We are led to a nucleosome gas that takes into account nucleosome breathing, i.e., transient unwrapping of nucleosomal DNA segments. This known biophysical property of nucleosomes rationalizes a "pressure"-induced dependence of the effective nucleosome size that is suggested by the data. By fitting this model to the data, we find an average energy cost for DNA unwrapping consistent with previous biophysical experiments. Although the available data are not sufficient to reconstruct chromatin remodeling mechanisms, a minimal model extension by one mechanism yields an "active nucleosome gas" that can rationalize the behavior of systems with reduced histone-DNA ratio and remodeler knockouts. We therefore establish a basis for a physical description of nucleosome patterns that can serve as a null model for sequence-specific effects at individual genes and in models of transcription regulation.

Hypoxia drives shared and distinct transcriptomic changes in two invasive glioma stem cell lines.

Glioblastoma (GBM) is the most common primary malignant cancer of the central nervous system. Insufficient oxygenation (hypoxia) has been linked to GBM invasion and aggression, leading to poor patient outcomes. Hypoxia induces gene expression for cellular adaptations. However, GBM is characterized by high intertumoral (molecular subtypes) and intratumoral heterogeneity (cell states), and it is not well understood to what extent hypoxia triggers patient-specific gene responses and cellular diversity in GBM. Here, we surveyed eight patient-derived GBM stem cell lines for invasion phenotypes in 3D culture, which identified two GBM lines showing increased invasiveness in response to hypoxia. RNA-seq analysis of the two patient GBM lines revealed a set of shared hypoxia response genes concerning glucose metabolism, angiogenesis, and autophagy, but also a large set of patient-specific hypoxia-induced genes featuring cell migration and anti-inflammation, highlighting intertumoral diversity of hypoxia responses in GBM. We further applied the Shared GBM Hypoxia gene signature to single cell RNA-seq datasets of glioma patients, which showed that hypoxic cells displayed a shift towards mesenchymal-like (MES) and astrocyte-like (AC) states. Interestingly, in response to hypoxia, tumor cells in IDH-mutant gliomas displayed a strong shift to the AC state, whereas tumor cells in IDH-wildtype gliomas mainly shifted to the MES state. This distinct hypoxia response of IDH-mutant gliomas may contribute to its more favorable prognosis. Our transcriptomic studies provide a basis for future approaches to better understand the diversity of hypoxic niches in gliomas.

A cellular and spatial atlas of TP53 -associated tissue remodeling in lung adenocarcinoma.

TP53 is the most frequently mutated gene across many cancers and is associated with shorter survival in lung adenocarcinoma (LUAD). To define how TP53 mutations affect the LUAD tumor microenvironment (TME), we constructed a multi-omic cellular and spatial tumor atlas of 23 treatment-naïve human lung tumors. We found that TP53 -mutant ( TP53 mut ) malignant cells lose alveolar identity and upregulate highly proliferative and entropic gene expression programs consistently across resectable LUAD patient tumors, genetically engineered mouse models, and cell lines harboring a wide spectrum of TP53 mutations. We further identified a multicellular tumor niche composed of SPP1 + macrophages and collagen-expressing fibroblasts that coincides with hypoxic, pro-metastatic expression programs in TP53 mut tumors. Spatially correlated angiostatic and immune checkpoint interactions, including CD274 - PDCD1 and PVR - TIGIT , are also enriched in TP53 mut LUAD tumors, which may influence response to checkpoint blockade therapy. Our methodology can be further applied to investigate mutation-specific TME changes in other cancers.

Myeloid progenitor dysregulation fuels immunosuppressive macrophages in tumors.

Monocyte-derived macrophages (mo-macs) drive immunosuppression in the tumor microenvironment (TME) and tumor-enhanced myelopoiesis in the bone marrow (BM) fuels these populations. Here, we performed paired transcriptome and chromatin analysis over the continuum of BM myeloid progenitors, circulating monocytes, and tumor-infiltrating mo-macs in mice and in patients with lung cancer to identify myeloid progenitor programs that fuel pro-tumorigenic mo-macs. Analyzing chromatin accessibility and histone mark changes, we show that lung tumors prime accessibility for Nfe2l2 (NRF2) in BM myeloid progenitors as a cytoprotective response to oxidative stress. NRF2 activity is sustained and increased during monocyte differentiation into mo-macs in the lung TME to regulate oxidative stress, in turn promoting metabolic adaptation, resistance to cell death, and contributing to immunosuppressive phenotype. NRF2 genetic deletion and pharmacological inhibition significantly reduced mo-macs' survival and immunosuppression in the TME, enabling NK and T cell therapeutic antitumor immunity and synergizing with checkpoint blockade strategies. Altogether, our study identifies a targetable epigenetic node of myeloid progenitor dysregulation that sustains immunoregulatory mo-macs in the TME.

Single cell view of tumor microenvironment gradients in pleural mesothelioma.

Immunotherapies have shown great promise in pleural mesothelioma (PM), yet most patients still do not achieve significant clinical response, highlighting the importance of improving understanding of the tumor microenvironment (TME). Here, we utilized high-throughput, single-cell RNA-sequencing to de novo identify 54 expression programs and construct a comprehensive cellular catalogue of the PM TME. We found four cancer-intrinsic programs associated with poor disease outcome and a novel fetal-like, endothelial cell population that likely responds to VEGF signaling and promotes angiogenesis. Throughout cellular compartments, we observe substantial difference in the TME associated with a cancer-intrinsic sarcomatoid signature, including enrichment in fetal-like endothelial cells, CXCL9+ macrophages, cytotoxic, exhausted, and regulatory T cells, which we validated using imaging and bulk deconvolution analyses on independent cohorts. Finally, we show, both computationally and experimentally, that NKG2A-HLA-E interaction between NK and tumor cells represents an important new therapeutic axis in PM, especially for epithelioid cases.

The role of nucleosome positioning in the evolution of gene regulation.

Chromatin organization plays a major role in gene regulation and can affect the function and evolution of new transcriptional programs. However, it can be difficult to decipher the basis of changes in chromatin organization and their functional effect on gene expression. Here, we present a large-scale comparative genomic analysis of the relationship between chromatin organization and gene expression, by measuring mRNA abundance and nucleosome positions genome-wide in 12 Hemiascomycota yeast species. We found substantial conservation of global and functional chromatin organization in all species, including prominent nucleosome-free regions (NFRs) at gene promoters, and distinct chromatin architecture in growth and stress genes. Chromatin organization has also substantially diverged in both global quantitative features, such as spacing between adjacent nucleosomes, and in functional groups of genes. Expression levels, intrinsic anti-nucleosomal sequences, and trans-acting chromatin modifiers all play important, complementary, and evolvable roles in determining NFRs. We identify five mechanisms that couple chromatin organization to evolution of gene regulation and have contributed to the evolution of respiro-fermentation and other key systems, including (1) compensatory evolution of alternative modifiers associated with conserved chromatin organization, (2) a gradual transition from constitutive to trans-regulated NFRs, (3) a loss of intrinsic anti-nucleosomal sequences accompanying changes in chromatin organization and gene expression, (4) re-positioning of motifs from NFRs to nucleosome-occluded regions, and (5) the expanded use of NFRs by paralogous activator-repressor pairs. Our study sheds light on the molecular basis of chromatin organization, and on the role of chromatin organization in the evolution of gene regulation.

Monocyte depletion enhances neutrophil influx and proneural to mesenchymal transition in glioblastoma.

Myeloid cells comprise the majority of immune cells in tumors, contributing to tumor growth and therapeutic resistance. Incomplete understanding of myeloid cells response to tumor driver mutation and therapeutic intervention impedes effective therapeutic design. Here, by leveraging CRISPR/Cas9-based genome editing, we generate a mouse model that is deficient of all monocyte chemoattractant proteins. Using this strain, we effectively abolish monocyte infiltration in genetically engineered murine models of de novo glioblastoma (GBM) and hepatocellular carcinoma (HCC), which show differential enrichment patterns for monocytes and neutrophils. Eliminating monocyte chemoattraction in monocyte enriched PDGFB-driven GBM invokes a compensatory neutrophil influx, while having no effect on Nf1-silenced GBM model. Single-cell RNA sequencing reveals that intratumoral neutrophils promote proneural-to-mesenchymal transition and increase hypoxia in PDGFB-driven GBM. We further demonstrate neutrophil-derived TNF-a directly drives mesenchymal transition in PDGFB-driven primary GBM cells. Genetic or pharmacological inhibiting neutrophils in HCC or monocyte-deficient PDGFB-driven and Nf1-silenced GBM models extend the survival of tumor-bearing mice. Our findings demonstrate tumor-type and genotype dependent infiltration and function of monocytes and neutrophils and highlight the importance of targeting them simultaneously for cancer treatments.

A Mutation-driven oncofetal regression fuels phenotypic plasticity in colorectal cancer.

Targeting cancer stem cells (CSCs) is crucial for effective cancer treatment 1 . However, the molecular mechanisms underlying resistance to LGR5 + CSCs depletion in colorectal cancer (CRC) 2,3 remain largely elusive. Here, we unveil the existence of a primitive cell state dubbed the oncofetal (OnF) state, which works in tandem with the LGR5 + stem cells (SCs) to fuel tumor evolution in CRC. OnF cells emerge early during intestinal tumorigenesis and exhibit features of lineage plasticity. Normally suppressed by the Retinoid X Receptor (RXR) in mature SCs, the OnF program is triggered by genetic deletion of the gatekeeper APC. We demonstrate that diminished RXR activity unlocks an epigenetic circuity governed by the cooperative action of YAP and AP1, leading to OnF reprogramming. This high-plasticity state is inherently resistant to conventional chemotherapies and its adoption by LGR5 + CSCs enables them to enter a drug-tolerant state. Furthermore, through phenotypic tracing and ablation experiments, we uncover a functional redundancy between the OnF and stem cell (SC) states and show that targeting both cellular states is essential for sustained tumor regression in vivo . Collectively, these findings establish a mechanistic foundation for developing effective combination therapies with enduring impact on CRC treatment.

A deep lung cell atlas reveals cytokine-mediated lineage switching of a rare cell progenitor of the human airway epithelium.

The human airway contains specialized rare epithelial cells whose roles in respiratory disease are not well understood. Ionocytes express the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR), while chemosensory tuft cells express asthma-associated alarmins. However, surprisingly, exceedingly few mature tuft cells have been identified in human lung cell atlases despite the ready identification of rare ionocytes and neuroendocrine cells. To identify human rare cell progenitors and define their lineage relationship to mature tuft cells, we generated a deep lung cell atlas containing 311,748 single cell RNA-Seq (scRNA-seq) profiles from discrete anatomic sites along the large and small airways and lung lobes of explanted donor lungs that could not be used for organ transplantation. Of 154,222 airway epithelial cells, we identified 687 ionocytes (0.45%) that are present in similar proportions in both large and small airways, suggesting that they may contribute to both large and small airways pathologies in CF. In stark contrast, we recovered only 3 mature tuft cells (0.002%). Instead, we identified rare bipotent progenitor cells that can give rise to both ionocytes and tuft cells, which we termed tuft-ionocyte progenitor cells (TIP cells). Remarkably, the cycling fraction of these TIP cells was comparable to that of basal stem cells. We used scRNA-seq and scATAC-seq to predict transcription factors that mark this novel rare cell progenitor population and define intermediate states during TIP cell lineage transitions en route to the differentiation of mature ionocytes and tuft cells. The default lineage of TIP cell descendants is skewed towards ionocytes, explaining the paucity of mature tuft cells in the human airway. However, Type 2 and Type 17 cytokines, associated with asthma and CF, diverted the lineage of TIP cell descendants in vitro , resulting in the differentiation of mature tuft cells at the expense of ionocytes. Consistent with this model of mature tuft cell differentiation, we identify mature tuft cells in a patient who died from an asthma flare. Overall, our findings suggest that the immune signaling pathways active in asthma and CF may skew the composition of disease-relevant rare cells and illustrate how deep atlases are required for identifying physiologically-relevant scarce cell populations.

A paracrine circuit of IL-1ß/IL-1R1 between myeloid and tumor cells drives genotype-dependent glioblastoma progression.

Monocytes and monocyte-derived macrophages (MDMs) from blood circulation infiltrate glioblastoma (GBM) and promote growth. Here, we show that PDGFB-driven GBM cells induce the expression of the potent proinflammatory cytokine IL-1ß in MDM, which engages IL-1R1 in tumor cells, activates the NF-κB pathway, and subsequently leads to induction of monocyte chemoattractant proteins (MCPs). Thus, a feedforward paracrine circuit of IL-1ß/IL-1R1 between tumors and MDM creates an interdependence driving PDGFB-driven GBM progression. Genetic loss or locally antagonizing IL-1ß/IL-1R1 leads to reduced MDM infiltration, diminished tumor growth, and reduced exhausted CD8+ T cells and thereby extends the survival of tumor-bearing mice. In contrast to IL-1ß, IL-1α exhibits antitumor effects. Genetic deletion of Il1a/b is associated with decreased recruitment of lymphoid cells and loss-of-interferon signaling in various immune populations and subsets of malignant cells and is associated with decreased survival time of PDGFB-driven tumor-bearing mice. In contrast to PDGFB-driven GBM, Nf1-silenced tumors have a constitutively active NF-κB pathway, which drives the expression of MCPs to recruit monocytes into tumors. These results indicate local antagonism of IL-1ß could be considered as an effective therapy specifically for proneural GBM.

Intratumoral dendritic cell-CD4+ T helper cell niches enable CD8+ T cell differentiation following PD-1 blockade in hepatocellular carcinoma.

Despite no apparent defects in T cell priming and recruitment to tumors, a large subset of T cell rich tumors fail to respond to immune checkpoint blockade (ICB). We leveraged a neoadjuvant anti-PD-1 trial in patients with hepatocellular carcinoma (HCC), as well as additional samples collected from patients treated off-label, to explore correlates of response to ICB within T cell-rich tumors. We show that ICB response correlated with the clonal expansion of intratumoral CXCL13+CH25H+IL-21+PD-1+CD4+ T helper cells ("CXCL13+ TH") and Granzyme K+ PD-1+ effector-like CD8+ T cells, whereas terminally exhausted CD39hiTOXhiPD-1hiCD8+ T cells dominated in nonresponders. CD4+ and CD8+ T cell clones that expanded post-treatment were found in pretreatment biopsies. Notably, PD-1+TCF-1+ (Progenitor-exhausted) CD8+ T cells shared clones mainly with effector-like cells in responders or terminally exhausted cells in nonresponders, suggesting that local CD8+ T cell differentiation occurs upon ICB. We found that these Progenitor CD8+ T cells interact with CXCL13+ TH within cellular triads around dendritic cells enriched in maturation and regulatory molecules, or "mregDC". These results suggest that discrete intratumoral niches that include mregDC and CXCL13+ TH control the differentiation of tumor-specific Progenitor exhasuted CD8+ T cells following ICB.

An atlas of late prenatal human neurodevelopment resolved by single-nucleus transcriptomics.

Late prenatal development of the human neocortex encompasses a critical period of gliogenesis and cortical expansion. However, systematic single-cell analyses to resolve cellular diversity and gliogenic lineages of the third trimester are lacking. Here, we present a comprehensive single-nucleus RNA sequencing atlas of over 200,000 nuclei derived from the proliferative germinal matrix and laminating cortical plate of 15 prenatal, non-pathological postmortem samples from 17 to 41 gestational weeks, and 3 adult controls. This dataset captures prenatal gliogenesis with high temporal resolution and is provided as a resource for further interrogation. Our computational analysis resolves greater complexity of glial progenitors, including transient glial intermediate progenitor cell (gIPC) and nascent astrocyte populations in the third trimester of human gestation. We use lineage trajectory and RNA velocity inference to further characterize specific gIPC subpopulations preceding both oligodendrocyte (gIPC-O) and astrocyte (gIPC-A) lineage differentiation. We infer unique transcriptional drivers and biological pathways associated with each developmental state, validate gIPC-A and gIPC-O presence within the human germinal matrix and cortical plate in situ, and demonstrate gIPC states being recapitulated across adult and pediatric glioblastoma tumors.

Single-nucleus cross-tissue molecular reference maps toward understanding disease gene function.

Understanding gene function and regulation in homeostasis and disease requires knowledge of the cellular and tissue contexts in which genes are expressed. Here, we applied four single-nucleus RNA sequencing methods to eight diverse, archived, frozen tissue types from 16 donors and 25 samples, generating a cross-tissue atlas of 209,126 nuclei profiles, which we integrated across tissues, donors, and laboratory methods with a conditional variational autoencoder. Using the resulting cross-tissue atlas, we highlight shared and tissue-specific features of tissue-resident cell populations; identify cell types that might contribute to neuromuscular, metabolic, and immune components of monogenic diseases and the biological processes involved in their pathology; and determine cell types and gene modules that might underlie disease mechanisms for complex traits analyzed by genome-wide association studies.

Author Correction: A single-cell and single-nucleus RNA-Seq toolbox for fresh and frozen human tumors.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.