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.

Systematic benchmarking of single-cell ATAC-sequencing protocols.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Single-cell multi-omics of mitochondrial DNA disorders reveals dynamics of purifying selection across human immune cells.

Pathogenic mutations in mitochondrial DNA (mtDNA) compromise cellular metabolism, contributing to cellular heterogeneity and disease. Diverse mutations are associated with diverse clinical phenotypes, suggesting distinct organ- and cell-type-specific metabolic vulnerabilities. Here we establish a multi-omics approach to quantify deletions in mtDNA alongside cell state features in single cells derived from six patients across the phenotypic spectrum of single large-scale mtDNA deletions (SLSMDs). By profiling 206,663 cells, we reveal the dynamics of pathogenic mtDNA deletion heteroplasmy consistent with purifying selection and distinct metabolic vulnerabilities across T-cell states in vivo and validate these observations in vitro. By extending analyses to hematopoietic and erythroid progenitors, we reveal mtDNA dynamics and cell-type-specific gene regulatory adaptations, demonstrating the context-dependence of perturbing mitochondrial genomic integrity. Collectively, we report pathogenic mtDNA heteroplasmy dynamics of individual blood and immune cells across lineages, demonstrating the power of single-cell multi-omics for revealing fundamental properties of mitochondrial genetics.

Positional influence on cellular transcriptional identity revealed through spatially segmented single-cell transcriptomics.

Single-cell RNA sequencing (scRNA-seq) is a powerful technique for describing cell states. Identifying the spatial arrangement of these states in tissues remains challenging, with the existing methods requiring niche methodologies and expertise. Here, we describe segmentation by exogenous perfusion (SEEP), a rapid and integrated method to link surface proximity and environment accessibility to transcriptional identity within three-dimensional (3D) disease models. The method utilizes the steady-state diffusion kinetics of a fluorescent dye to establish a gradient along the radial axis of disease models. Classification of sample layers based on dye accessibility enables dissociated and sorted cells to be characterized by transcriptomic and regional identities. Using SEEP, we analyze spheroid, organoid, and in vivo tumor models of high-grade serous ovarian cancer (HGSOC). The results validate long-standing beliefs about the relationship between cell state and position while revealing new concepts regarding how spatially unique microenvironments influence the identity of individual cells within tumors.

Mitochondrial single-cell ATAC-seq for high-throughput multi-omic detection of mitochondrial genotypes and chromatin accessibility.

Natural sequence variation within mitochondrial DNA (mtDNA) contributes to human phenotypes and may serve as natural genetic markers in human cells for clonal and lineage tracing. We recently developed a single-cell multi-omic approach, called 'mitochondrial single-cell assay for transposase-accessible chromatin with sequencing' (mtscATAC-seq), enabling concomitant high-throughput mtDNA genotyping and accessible chromatin profiling. Specifically, our technique allows the mitochondrial genome-wide inference of mtDNA variant heteroplasmy along with information on cell state and accessible chromatin variation in individual cells. Leveraging somatic mtDNA mutations, our method further enables inference of clonal relationships among native ex vivo-derived human cells not amenable to genetic engineering-based clonal tracing approaches. Here, we provide a step-by-step protocol for the use of mtscATAC-seq, including various cell-processing and flow cytometry workflows, by using primary hematopoietic cells, subsequent single-cell genomic library preparation and sequencing that collectively take ~3-4 days to complete. We discuss experimental and computational data quality control metrics and considerations for the extension to other mammalian tissues. Overall, mtscATAC-seq provides a broadly applicable platform to map clonal relationships between cells in human tissues, investigate fundamental aspects of mitochondrial genetics and enable additional modes of multi-omic discovery.

Congenital anemia reveals distinct targeting mechanisms for master transcription factor GATA1.

Master regulators, such as the hematopoietic transcription factor (TF) GATA1, play an essential role in orchestrating lineage commitment and differentiation. However, the precise mechanisms by which such TFs regulate transcription through interactions with specific cis-regulatory elements remain incompletely understood. Here, we describe a form of congenital hemolytic anemia caused by missense mutations in an intrinsically disordered region of GATA1, with a poorly understood role in transcriptional regulation. Through integrative functional approaches, we demonstrate that these mutations perturb GATA1 transcriptional activity by partially impairing nuclear localization and selectively altering precise chromatin occupancy by GATA1. These alterations in chromatin occupancy and concordant chromatin accessibility changes alter faithful gene expression, with failure to both effectively silence and activate select genes necessary for effective terminal red cell production. We demonstrate how disease-causing mutations can reveal regulatory mechanisms that enable the faithful genomic targeting of master TFs during cellular differentiation.

Epigenetic encoding, heritability and plasticity of glioma transcriptional cell states.

Single-cell RNA sequencing has revealed extensive transcriptional cell state diversity in cancer, often observed independently of genetic heterogeneity, raising the central question of how malignant cell states are encoded epigenetically. To address this, here we performed multiomics single-cell profiling-integrating DNA methylation, transcriptome and genotype within the same cells-of diffuse gliomas, tumors characterized by defined transcriptional cell state diversity. Direct comparison of the epigenetic profiles of distinct cell states revealed key switches for state transitions recapitulating neurodevelopmental trajectories and highlighted dysregulated epigenetic mechanisms underlying gliomagenesis. We further developed a quantitative framework to directly measure cell state heritability and transition dynamics based on high-resolution lineage trees in human samples. We demonstrated heritability of malignant cell states, with key differences in hierarchal and plastic cell state architectures in IDH-mutant glioma versus IDH-wild-type glioblastoma, respectively. This work provides a framework anchoring transcriptional cancer cell states in their epigenetic encoding, inheritance and transition dynamics.

Single-cell meta-analysis of SARS-CoV-2 entry genes across tissues and demographics.

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.

Skin-resident innate lymphoid cells converge on a pathogenic effector state.

Tissue-resident innate lymphoid cells (ILCs) help sustain barrier function and respond to local signals. ILCs are traditionally classified as ILC1, ILC2 or ILC3 on the basis of their expression of specific transcription factors and cytokines1. In the skin, disease-specific production of ILC3-associated cytokines interleukin (IL)-17 and IL-22 in response to IL-23 signalling contributes to dermal inflammation in psoriasis. However, it is not known whether this response is initiated by pre-committed ILCs or by cell-state transitions. Here we show that the induction of psoriasis in mice by IL-23 or imiquimod reconfigures a spectrum of skin ILCs, which converge on a pathogenic ILC3-like state. Tissue-resident ILCs were necessary and sufficient, in the absence of circulatory ILCs, to drive pathology. Single-cell RNA-sequencing (scRNA-seq) profiles of skin ILCs along a time course of psoriatic inflammation formed a dense transcriptional continuum-even at steady state-reflecting fluid ILC states, including a naive or quiescent-like state and an ILC2 effector state. Upon disease induction, the continuum shifted rapidly to span a mixed, ILC3-like subset also expressing cytokines characteristic of ILC2s, which we inferred as arising through multiple trajectories. We confirmed the transition potential of quiescent-like and ILC2 states using in vitro experiments, single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) and in vivo fate mapping. Our results highlight the range and flexibility of skin ILC responses, suggesting that immune activities primed in healthy tissues dynamically adapt to provocations and, left unchecked, drive pathological remodelling.

Massively parallel single-cell mitochondrial DNA genotyping and chromatin profiling.

Natural mitochondrial DNA (mtDNA) mutations enable the inference of clonal relationships among cells. mtDNA can be profiled along with measures of cell state, but has not yet been combined with the massively parallel approaches needed to tackle the complexity of human tissue. Here, we introduce a high-throughput, droplet-based mitochondrial single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), a method that combines high-confidence mtDNA mutation calling in thousands of single cells with their concomitant high-quality accessible chromatin profile. This enables the inference of mtDNA heteroplasmy, clonal relationships, cell state and accessible chromatin variation in individual cells. We reveal single-cell variation in heteroplasmy of a pathologic mtDNA variant, which we associate with intra-individual chromatin variability and clonal evolution. We clonally trace thousands of cells from cancers, linking epigenomic variability to subclonal evolution, and infer cellular dynamics of differentiating hematopoietic cells in vitro and in vivo. Taken together, our approach enables the study of cellular population dynamics and clonal properties in vivo.

Myocyte Specific Upregulation of ACE2 in Cardiovascular Disease: Implications for SARS-CoV-2 mediated myocarditis.

Coronavirus disease 2019 (COVID-19) is a global pandemic caused by a novel severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). SARS-CoV-2 infection of host cells occurs predominantly via binding of the viral surface spike protein to the human angiotensin-converting enzyme 2 (ACE2) receptor. Hypertension and pre-existing cardiovascular disease are risk factors for morbidity from COVID-19, and it remains uncertain whether the use of angiotensin converting enzyme inhibitors (ACEi) or angiotensin receptor blockers (ARB) impacts infection and disease. Here, we aim to shed light on this question by assessing ACE2 expression in normal and diseased human myocardial samples profiled by bulk and single nucleus RNA-seq.

SARS-CoV-2 Receptor ACE2 Is an Interferon-Stimulated Gene in Human Airway Epithelial Cells and Is Detected in Specific Cell Subsets across Tissues.

There is pressing urgency to understand the pathogenesis of the severe acute respiratory syndrome coronavirus clade 2 (SARS-CoV-2), which causes the disease COVID-19. SARS-CoV-2 spike (S) protein binds angiotensin-converting enzyme 2 (ACE2), and in concert with host proteases, principally transmembrane serine protease 2 (TMPRSS2), promotes cellular entry. The cell subsets targeted by SARS-CoV-2 in host tissues and the factors that regulate ACE2 expression remain unknown. Here, we leverage human, non-human primate, and mouse single-cell RNA-sequencing (scRNA-seq) datasets across health and disease to uncover putative targets of SARS-CoV-2 among tissue-resident cell subsets. We identify ACE2 and TMPRSS2 co-expressing cells within lung type II pneumocytes, ileal absorptive enterocytes, and nasal goblet secretory cells. Strikingly, we discovered that ACE2 is a human interferon-stimulated gene (ISG) in vitro using airway epithelial cells and extend our findings to in vivo viral infections. Our data suggest that SARS-CoV-2 could exploit species-specific interferon-driven upregulation of ACE2, a tissue-protective mediator during lung injury, to enhance infection.

Transcriptional States and Chromatin Accessibility Underlying Human Erythropoiesis.

Human erythropoiesis serves as a paradigm of physiologic cellular differentiation. This process is also of considerable interest for better understanding anemias and identifying new therapies. Here, we apply deep transcriptomic and accessible chromatin profiling to characterize a faithful ex vivo human erythroid differentiation system from hematopoietic stem and progenitor cells. We reveal stage-specific transcriptional states and chromatin accessibility during various stages of erythropoiesis, including 14,260 differentially expressed genes and 63,659 variably accessible chromatin peaks. Our analysis suggests differentiation stage-predominant roles for specific master regulators, including GATA1 and KLF1. We integrate chromatin profiles with common and rare genetic variants associated with erythroid cell traits and diseases, finding that variants regulating different erythroid phenotypes likely act at variable points during differentiation. In addition, we identify a regulator of terminal erythropoiesis, TMCC2, more broadly illustrating the value of this comprehensive analysis to improve our understanding of erythropoiesis in health and disease.

Lineage Tracing in Humans Enabled by Mitochondrial Mutations and Single-Cell Genomics.

Lineage tracing provides key insights into the fate of individual cells in complex organisms. Although effective genetic labeling approaches are available in model systems, in humans, most approaches require detection of nuclear somatic mutations, which have high error rates, limited scale, and do not capture cell state information. Here, we show that somatic mutations in mtDNA can be tracked by single-cell RNA or assay for transposase accessible chromatin (ATAC) sequencing. We leverage somatic mtDNA mutations as natural genetic barcodes and demonstrate their utility as highly accurate clonal markers to infer cellular relationships. We track native human cells both in vitro and in vivo and relate clonal dynamics to gene expression and chromatin accessibility. Our approach should allow clonal tracking at a 1,000-fold greater scale than with nuclear genome sequencing, with simultaneous information on cell state, opening the way to chart cellular dynamics in human health and disease.

Lipid-coated nanocapillaries for DNA sensing.

We report a simple and efficient way to accomplish the chemical modification of glass nanopores by means of lipid self-assembly. Lipid coating improves the success rate of these glass nanopores as biosensors to detect λ-DNA.

α(V)ß(3) integrin-targeted PLGA-PEG nanoparticles for enhanced anti-tumor efficacy of a Pt(IV) prodrug.

Targeted delivery of therapeutics to tumor neovasculature is potentially a powerful approach for selective cancer treatment. Integrins are heterodimeric transmembrane proteins involved in cell adhesion and cell signaling, and their expression is commonly upregulated in cancers and inflammatory diseases. The α(v)ß(3) integrin is differentially upregulated on angiogenic endothelial cells as well as on many cancer cells. Here we demonstrate the differential targeting of cisplatin prodrug-encapsulated poly(d,l-lactic-co-glycolic acid)-block-polyethylene glycol (PLGA-PEG) nanoparticles (NPs) to the α(v)ß(3) integrin on cancer cells using the cyclic pentapeptide c(RGDfK). Cisplatin is one of the most widely used anticancer drugs, and approaches that can improve its therapeutic index are of broad importance. The RGD-targeted Pt(IV)-encapsulated NPs displayed enhanced cytotoxicity as compared to cisplatin administered in its conventional dosage form in model prostate and breast cancer epithelial cells in vitro. Cytotoxicities were also elevated in comparison to those of previously reported systems, a small molecule Pt(IV)-RGD conjugate and a Pt(IV) nanoscale coordination polymer carrying RGD moieties. This result encouraged us also to evaluate the anticancer effect of the new construct in an animal model. The RGD-targeted PLGA-PEG NPs were more efficacious and better tolerated by comparison to cisplatin in an orthotopic human breast cancer xenograft model in vivo.