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.

Transcriptional and epigenetic dynamics during specification of human embryonic stem cells.

Differentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

Class IIa HDAC4 and HDAC7 cooperatively regulate gene transcription in Th17 cell differentiation.

Class II histone deacetylases (HDACs) are important in regulation of gene transcription during T cell development. However, our understanding of their cell-specific functions is limited. In this study, we reveal that class IIa Hdac4 and Hdac7 (Hdac4/7) are selectively induced in transcription, guiding the lineage-specific differentiation of mouse T-helper 17 (Th17) cells from naive CD4+ T cells. Importantly, Hdac4/7 are functionally dispensable in other Th subtypes. Mechanistically, Hdac4 interacts with the transcription factor (TF) JunB, facilitating the transcriptional activation of Th17 signature genes such as Il17a/f. Conversely, Hdac7 collaborates with the TF Aiolos and Smrt/Ncor1-Hdac3 corepressors to repress transcription of Th17 negative regulators, including Il2, in Th17 cell differentiation. Inhibiting Hdac4/7 through pharmacological or genetic methods effectively mitigates Th17 cell-mediated intestinal inflammation in a colitis mouse model. Our study uncovers molecular mechanisms where HDAC4 and HDAC7 function distinctively yet cooperatively in regulating ordered gene transcription during Th17 cell differentiation. These findings suggest a potential therapeutic strategy of targeting HDAC4/7 for treating Th17-related inflammatory diseases, such as ulcerative colitis.

An essential role for Cmtr2 in mammalian embryonic development.

CMTR2 is an mRNA cap methyltransferase with poorly understood physiological functions. It catalyzes 2'-O-ribose methylation of the second transcribed nucleotide of mRNAs, potentially serving to mark RNAs as "self" to evade the cellular innate immune response. Here we analyze the consequences of Cmtr2 deficiency in mice. We discover that constitutive deletion of Cmtr2 results in mouse embryos that die during mid-gestation, exhibiting defects in embryo size, placental malformation and yolk sac vascularization. Endothelial cell deletion of Cmtr2 in mice results in vascular and hematopoietic defects, and perinatal lethality. Detailed characterization of the constitutive Cmtr2 KO phenotype shows an activation of the p53 pathway and decreased proliferation, but no evidence of interferon pathway activation. In summary, our study reveals the essential roles of Cmtr2 in mammalian cells beyond its immunoregulatory function.

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.

Dissecting neural differentiation regulatory networks through epigenetic footprinting.

Models derived from human pluripotent stem cells that accurately recapitulate neural development in vitro and allow for the generation of specific neuronal subtypes are of major interest to the stem cell and biomedical community. Notch signalling, particularly through the Notch effector HES5, is a major pathway critical for the onset and maintenance of neural progenitor cells in the embryonic and adult nervous system. Here we report the transcriptional and epigenomic analysis of six consecutive neural progenitor cell stages derived from a HES5::eGFP reporter human embryonic stem cell line. Using this system, we aimed to model cell-fate decisions including specification, expansion and patterning during the ontogeny of cortical neural stem and progenitor cells. In order to dissect regulatory mechanisms that orchestrate the stage-specific differentiation process, we developed a computational framework to infer key regulators of each cell-state transition based on the progressive remodelling of the epigenetic landscape and then validated these through a pooled short hairpin RNA screen. We were also able to refine our previous observations on epigenetic priming at transcription factor binding sites and suggest here that they are mediated by combinations of core and stage-specific factors. Taken together, we demonstrate the utility of our system and outline a general framework, not limited to the context of the neural lineage, to dissect regulatory circuits of differentiation.

Automated, high-throughput derivation, characterization and differentiation of induced pluripotent stem cells.

Induced pluripotent stem cells (iPSCs) are an essential tool for modeling how causal genetic variants impact cellular function in disease, as well as an emerging source of tissue for regenerative medicine. The preparation of somatic cells, their reprogramming and the subsequent verification of iPSC pluripotency are laborious, manual processes limiting the scale and reproducibility of this technology. Here we describe a modular, robotic platform for iPSC reprogramming enabling automated, high-throughput conversion of skin biopsies into iPSCs and differentiated cells with minimal manual intervention. We demonstrate that automated reprogramming and the pooled selection of polyclonal pluripotent cells results in high-quality, stable iPSCs. These lines display less line-to-line variation than either manually produced lines or lines produced through automation followed by single-colony subcloning. The robotic platform we describe will enable the application of iPSCs to population-scale biomedical problems including the study of complex genetic diseases and the development of personalized medicines.

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

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.