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.

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 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 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.

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.

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.

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.

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.

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.

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

Single-cell genomics is essential to chart tumor ecosystems. Although single-cell RNA-Seq (scRNA-Seq) profiles RNA from cells dissociated from fresh tumors, single-nucleus RNA-Seq (snRNA-Seq) is needed to profile frozen or hard-to-dissociate tumors. Each requires customization to different tissue and tumor types, posing a barrier to adoption. Here, we have developed a systematic toolbox for profiling fresh and frozen clinical tumor samples using scRNA-Seq and snRNA-Seq, respectively. We analyzed 216,490 cells and nuclei from 40 samples across 23 specimens spanning eight tumor types of varying tissue and sample characteristics. We evaluated protocols by cell and nucleus quality, recovery rate and cellular composition. scRNA-Seq and snRNA-Seq from matched samples recovered the same cell types, but at different proportions. Our work provides guidance for studies in a broad range of tumors, including criteria for testing and selecting methods from the toolbox for other tumors, thus paving the way for charting tumor atlases.

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.

Analysis of chromatin accessibility uncovers TEAD1 as a regulator of migration in human glioblastoma.

The intrinsic drivers of migration in glioblastoma (GBM) are poorly understood. To better capture the native molecular imprint of GBM and its developmental context, here we isolate human stem cell populations from GBM (GSC) and germinal matrix tissues and map their chromatin accessibility via ATAC-seq. We uncover two distinct regulatory GSC signatures, a developmentally shared/proliferative and a tumor-specific/migratory one in which TEAD1/4 motifs are uniquely overrepresented. Using ChIP-PCR, we validate TEAD1 trans occupancy at accessibility sites within AQP4, EGFR, and CDH4. To further characterize TEAD's functional role in GBM, we knockout TEAD1 or TEAD4 in patient-derived GBM lines using CRISPR-Cas9. TEAD1 ablation robustly diminishes migration, both in vitro and in vivo, and alters migratory and EMT transcriptome signatures with consistent downregulation of its target AQP4. TEAD1 overexpression restores AQP4 expression, and both TEAD1 and AQP4 overexpression rescue migratory deficits in TEAD1-knockout cells, implicating a direct regulatory role for TEAD1-AQP4 in GBM migration.

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.

Global delay in nascent strand DNA methylation.

Cytosine methylation is widespread among organisms and essential for mammalian development. In line with early postulations of an epigenetic role in gene regulation, symmetric CpG methylation can be mitotically propagated over many generations with extraordinarily high fidelity. Here, we combine BrdU labeling and immunoprecipitation with genome-wide bisulfite sequencing to explore the inheritance of cytosine methylation onto newly replicated DNA in human cells. Globally, we observe a pronounced lag between the copying of genetic and epigenetic information in embryonic stem cells that is reconsolidated within hours to accomplish faithful mitotic transmission. Populations of arrested cells show a global reduction of lag-induced intermediate CpG methylation when compared to proliferating cells, whereas sites of transcription factor engagement appear cell-cycle invariant. Alternatively, the cancer cell line HCT116 preserves global epigenetic heterogeneity independently of cell-cycle arrest. Taken together, our data suggest that heterogeneous methylation largely reflects asynchronous proliferation, but is intrinsic to actively engaged cis-regulatory elements and cancer.

Author Correction: Global delay in nascent strand DNA methylation.

Following online publication of this article, the Gene Expression Omnibus records corresponding to accession codes GSM2406773, MN-d6, and GSM2406772, MN-d14, listed in the data availability statement were deleted. The data are now available under accession codes GSM3039355, WGBS_hESC_WT_D6_R4 (MN day 6), and GSM3039351, WGBS_hESC_WT_D14_R4 (MN day 14), and the data availability statement has been updated with the new accession codes in the HTML and PDF versions of the article.

Induced Pluripotent Stem Cell Differentiation Enables Functional Validation of GWAS Variants in Metabolic Disease.

Genome-wide association studies (GWAS) have highlighted a large number of genetic variants with potential disease association, but functional analysis remains a challenge. Here we describe an approach to functionally validate identified variants through differentiation of induced pluripotent stem cells (iPSCs) to study cellular pathophysiology. We collected peripheral blood cells from Framingham Heart Study participants and reprogrammed them to iPSCs. We then differentiated 68 iPSC lines into hepatocytes and adipocytes to investigate the effect of the 1p13 rs12740374 variant on cardiometabolic disease phenotypes via transcriptomics and metabolomic signatures. We observed a clear association between rs12740374 and lipid accumulation and gene expression in differentiated hepatocytes, in particular, expression of SORT1, CELSR2, and PSRC1, consistent with previous analyses of this variant using other approaches. Initial investigation of additional SNPs also highlighted correlations with gene expression. These findings suggest that iPSC-based population studies hold promise as tools for the functional validation of GWAS variants.

A comparison of genetically matched cell lines reveals the equivalence of human iPSCs and ESCs.

The equivalence of human induced pluripotent stem cells (hiPSCs) and human embryonic stem cells (hESCs) remains controversial. Here we use genetically matched hESC and hiPSC lines to assess the contribution of cellular origin (hESC vs. hiPSC), the Sendai virus (SeV) reprogramming method and genetic background to transcriptional and DNA methylation patterns while controlling for cell line clonality and sex. We find that transcriptional and epigenetic variation originating from genetic background dominates over variation due to cellular origin or SeV infection. Moreover, the 49 differentially expressed genes we detect between genetically matched hESCs and hiPSCs neither predict functional outcome nor distinguish an independently derived, larger set of unmatched hESC and hiPSC lines. We conclude that hESCs and hiPSCs are molecularly and functionally equivalent and cannot be distinguished by a consistent gene expression signature. Our data further imply that genetic background variation is a major confounding factor for transcriptional and epigenetic comparisons of pluripotent cell lines, explaining some of the previously observed differences between genetically unmatched hESCs and hiPSCs.