Airway hillocks are injury-resistant reservoirs of unique plastic stem cells.

Airway hillocks are stratified epithelial structures of unknown function1. Hillocks persist for months and have a unique population of basal stem cells that express genes associated with barrier function and cell adhesion. Hillock basal stem cells continually replenish overlying squamous barrier cells. They exhibit dramatically higher turnover than the abundant, largely quiescent classic pseudostratified airway epithelium. Hillocks resist a remarkably broad spectrum of injuries, including toxins, infection, acid and physical injury because hillock squamous cells shield underlying hillock basal stem cells from injury. Hillock basal stem cells are capable of massive clonal expansion that is sufficient to resurface denuded airway, and eventually regenerate normal airway epithelium with each of its six component cell types. Hillock basal stem cells preferentially stratify and keratinize in the setting of retinoic acid signalling inhibition, a known cause of squamous metaplasia2,3. Here we show that mouse hillock expansion is the cause of vitamin A deficiency-induced squamous metaplasia. Finally, we identify human hillocks whose basal stem cells generate functional squamous barrier structures in culture. The existence of hillocks reframes our understanding of airway epithelial regeneration. Furthermore, we show that hillocks are one origin of 'squamous metaplasia', which is long thought to be a precursor of lung cancer.

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.

Cystic Fibrosis and the Cells of the Airway Epithelium: What Are Ionocytes and What Do They Do?

Cystic fibrosis (CF) is caused by defects in an anion channel, the cystic fibrosis transmembrane conductance regulator (CFTR). Recently, a new airway epithelial cell type has been discovered and dubbed the pulmonary ionocyte. Unexpectedly, these ionocytes express higher levels of CFTR than any other airway epithelial cell type. However, ionocytes are not the sole CFTR-expressing airway epithelial cells, and CF-associated disease genes are in fact expressed in multiple airway epithelial cell types. The experimental depletion of ionocytes perturbs epithelial physiology in the mouse trachea, but the role of these rare cells in the pathogenesis of human CF remains mysterious. Ionocytes have been described in diverse tissues(kidney and inner ear) and species (frog and fish). We draw on these prior studies to suggest potential roles of airway ionocytes in health and disease. A complete understanding of ionocytes in the mammalian airway will ultimately depend on cell type-specific genetic manipulation.

Longitudinal proteomic analysis of severe COVID-19 reveals survival-associated signatures, tissue-specific cell death, and cell-cell interactions.

Mechanisms underlying severe coronavirus disease 2019 (COVID-19) disease remain poorly understood. We analyze several thousand plasma proteins longitudinally in 306 COVID-19 patients and 78 symptomatic controls, uncovering immune and non-immune proteins linked to COVID-19. Deconvolution of our plasma proteome data using published scRNA-seq datasets reveals contributions from circulating immune and tissue cells. Sixteen percent of patients display reduced inflammation yet comparably poor outcomes. Comparison of patients who died to severely ill survivors identifies dynamic immune-cell-derived and tissue-associated proteins associated with survival, including exocrine pancreatic proteases. Using derived tissue-specific and cell-type-specific intracellular death signatures, cellular angiotensin-converting enzyme 2 (ACE2) expression, and our data, we infer whether organ damage resulted from direct or indirect effects of infection. We propose a model in which interactions among myeloid, epithelial, and T cells drive tissue damage. These datasets provide important insights and a rich resource for analysis of mechanisms of severe COVID-19 disease.

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.

Plasma proteomics reveals tissue-specific cell death and mediators of cell-cell interactions in severe COVID-19 patients.

COVID-19 has caused over 1 million deaths globally, yet the cellular mechanisms underlying severe disease remain poorly understood. By analyzing several thousand plasma proteins in 306 COVID-19 patients and 78 symptomatic controls over serial timepoints using two complementary approaches, we uncover COVID-19 host immune and non-immune proteins not previously linked to this disease. Integration of plasma proteomics with nine published scRNAseq datasets shows that SARS-CoV-2 infection upregulates monocyte/macrophage, plasmablast, and T cell effector proteins. By comparing patients who died to severely ill patients who survived, we identify dynamic immunomodulatory and tissue-associated proteins associated with survival, providing insights into which host responses are beneficial and which are detrimental to survival. We identify intracellular death signatures from specific tissues and cell types, and by associating these with angiotensin converting enzyme 2 (ACE2) expression, we map tissue damage associated with severe disease and propose which damage results from direct viral infection rather than from indirect effects of illness. We find that disease severity in lung tissue is driven by myeloid cell phenotypes and cell-cell interactions with lung epithelial cells and T cells. Based on these results, we propose a model of immune and epithelial cell interactions that drive cell-type specific and tissue-specific damage in severe COVID-19.

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.

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.

Memory of Divisional History Directs the Continuous Process of Primitive Hematopoietic Lineage Commitment.

Hematopoietic stem cells (HSCs) exist in a dormant state and progressively lose regenerative potency as they undergo successive divisions. Why this functional decline occurs and how this information is encoded is unclear. To better understand how this information is stored, we performed RNA sequencing on HSC populations differing only in their divisional history. Comparative analysis revealed that genes upregulated with divisions are enriched for lineage genes and regulated by cell-cycle-associated transcription factors, suggesting that proliferation itself drives lineage priming. Downregulated genes are, however, associated with an HSC signature and targeted by the Polycomb Repressive Complex 2 (PRC2). The PRC2 catalytic subunits Ezh1 and Ezh2 promote and suppress the HSC state, respectively, and successive divisions cause a switch from Ezh1 to Ezh2 dominance. We propose that cell divisions drive lineage priming and Ezh2 accumulation, which represses HSC signature genes to consolidate information on divisional history into memory.

The Human Lung Cell Atlas: A High-Resolution Reference Map of the Human Lung in Health and Disease.

Lung disease accounts for every sixth death globally. Profiling the molecular state of all lung cell types in health and disease is currently revolutionizing the identification of disease mechanisms and will aid the design of novel diagnostic and personalized therapeutic regimens. Recent progress in high-throughput techniques for single-cell genomic and transcriptomic analyses has opened up new possibilities to study individual cells within a tissue, classify these into cell types, and characterize variations in their molecular profiles as a function of genetics, environment, cell-cell interactions, developmental processes, aging, or disease. Integration of these cell state definitions with spatial information allows the in-depth molecular description of cellular neighborhoods and tissue microenvironments, including the tissue resident structural and immune cells, the tissue matrix, and the microbiome. The Human Cell Atlas consortium aims to characterize all cells in the healthy human body and has prioritized lung tissue as one of the flagship projects. Here, we present the rationale, the approach, and the expected impact of a Human Lung Cell Atlas.

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.

Feedback control of pluripotency in embryonic stem cells: Signaling, transcription and epigenetics.

Embryonic stem cells (ESCs) can proliferate and self-renew, maintaining their pluripotency status in vitro for a long period of time. Pluripotent states of ESCs in vitro are supported by a network of signaling, transcriptional and epigenetic regulatory interactions known as the pluripotency gene regulatory network (PGRN). Despite extensive investigation of the network, the exact order of regulatory links and many structural features of the network are still missing. Analysis of published data and literature reveals numerous PGRN components regulating each other in a mutual fashion, thus creating multiple regulatory feedback control circuits. Here we consider possible organizational features of PGRN and describe examples representing known feedback control loops in the context of mouse ESCs. We discuss how the feedback control interactions can contribute to learning behavior and dynamic responses of pluripotency gene network to changing environments.

A Conserved Distal Lung Regenerative Pathway in Acute Lung Injury.

Improved tools have led to a burgeoning understanding of lung regeneration in mice, but it is not yet known how these insights may be relevant to acute lung injury in humans. We report in detail two cases of fulminant idiopathic acute lung injury requiring extracorporeal membrane oxygenation in previously healthy young adults with acute respiratory distress syndrome, one of whom required lung transplantation. Biopsy specimens showed diffuse alveolar injury with a striking paucity of alveolar epithelial regeneration, rare hyaline membranes, and diffuse contiguous airspace lining by macrophages. This novel constellation was termed diffuse alveolar injury with delayed epithelization. In addition, mirroring data from murine models of lung injury/regeneration, peribronchiolar basaloid pods (previously described as squamous metaplasia) and ciliated bronchiolarization were identified in these patients and in 39% of 57 historical cases with diffuse alveolar damage. These findings demonstrate a common and clinically relevant human disease correlate for murine models of severe acute lung injury. Evidence suggests that peribronchiolar basaloid pods and bronchiolarization are related spatially and temporally and likely represent overlapping sequential stages of the response to severe distal airway injury.

Tips from the embryonic lung.

A new source of progenitor cells can now be used to study hidden aspects of human lung development and pediatric lung disease.

Transient HES5 Activity Instructs Mesodermal Cells toward a Cardiac Fate.

Notch signaling plays a role in specifying a cardiac fate but the downstream effectors remain unknown. In this study we implicate the Notch downstream effector HES5 in cardiogenesis. We show transient Hes5 expression in early mesoderm of gastrulating embryos and demonstrate, by loss and gain-of-function experiments in mouse embryonic stem cells, that HES5 favors cardiac over primitive erythroid fate. Hes5 overexpression promotes upregulation of the cardiac gene Isl1, while the hematopoietic regulator Scl is downregulated. Moreover, whereas a pulse of Hes5 instructs cardiac commitment, sustained expression after lineage specification impairs progression of differentiation to contracting cardiomyocytes. These findings establish a role for HES5 in cardiogenesis and provide insights into the early cardiac molecular network.

A dual molecular analogue tuner for dissecting protein function in mammalian cells.

Loss-of-function studies are fundamental for dissecting gene function. Yet, methods to rapidly and effectively perturb genes in mammalian cells, and particularly in stem cells, are scarce. Here we present a system for simultaneous conditional regulation of two different proteins in the same mammalian cell. This system harnesses the plant auxin and jasmonate hormone-induced degradation pathways, and is deliverable with only two lentiviral vectors. It combines RNAi-mediated silencing of two endogenous proteins with the expression of two exogenous proteins whose degradation is induced by external ligands in a rapid, reversible, titratable and independent manner. By engineering molecular tuners for NANOG, CHK1, p53 and NOTCH1 in mammalian stem cells, we have validated the applicability of the system and demonstrated its potential to unravel complex biological processes.

Single-Cell Analyses of ESCs Reveal Alternative Pluripotent Cell States and Molecular Mechanisms that Control Self-Renewal.

Analyses of gene expression in single mouse embryonic stem cells (mESCs) cultured in serum and LIF revealed the presence of two distinct cell subpopulations with individual gene expression signatures. Comparisons with published data revealed that cells in the first subpopulation are phenotypically similar to cells isolated from the inner cell mass (ICM). In contrast, cells in the second subpopulation appear to be more mature. Pluripotency Gene Regulatory Network (PGRN) reconstruction based on single-cell data and published data suggested antagonistic roles for Oct4 and Nanog in the maintenance of pluripotency states. Integrated analyses of published genomic binding (ChIP) data strongly supported this observation. Certain target genes alternatively regulated by OCT4 and NANOG, such as Sall4 and Zscan10, feed back into the top hierarchical regulator Oct4. Analyses of such incoherent feedforward loops with feedback (iFFL-FB) suggest a dynamic model for the maintenance of mESC pluripotency and self-renewal.

Tbx3 Controls Dppa3 Levels and Exit from Pluripotency toward Mesoderm.

Tbx3, a member of the T-box family, plays important roles in development, stem cells, nuclear reprogramming, and cancer. Loss of Tbx3 induces differentiation in mouse embryonic stem cells (mESCs). However, we show that mESCs exist in an alternate stable pluripotent state in the absence of Tbx3. In-depth transcriptome analysis of this mESC state reveals Dppa3 as a direct downstream target of Tbx3. Also, Tbx3 facilitates the cell fate transition from pluripotent cells to mesoderm progenitors by directly repressing Wnt pathway members required for differentiation. Wnt signaling regulates differentiation of mESCs into mesoderm progenitors and helps to maintain a naive pluripotent state. We show that Tbx3, a downstream target of Wnt signaling, fine tunes these divergent roles of Wnt signaling in mESCs. In conclusion, we identify a signaling-TF axis that controls the exit of mESCs from a self-renewing pluripotent state toward mesoderm differentiation.

Tex10 Coordinates Epigenetic Control of Super-Enhancer Activity in Pluripotency and Reprogramming.

Super-enhancers (SEs) are large clusters of transcriptional enhancers that are co-occupied by multiple lineage-specific transcription factors driving expression of genes that define cell identity. In embryonic stem cells (ESCs), SEs are highly enriched for the core pluripotency factors Oct4, Sox2, and Nanog. In this study, we sought to dissect the molecular control mechanism of SE activity in pluripotency and reprogramming. Starting from a protein interaction network surrounding Sox2, we identified Tex10 as a key pluripotency factor that plays a functionally significant role in ESC self-renewal, early embryo development, and reprogramming. Tex10 is enriched at SEs in a Sox2-dependent manner and coordinates histone acetylation and DNA demethylation at SEs. Tex10 activity is also important for pluripotency and reprogramming in human cells. Our study therefore highlights Tex10 as a core component of the pluripotency network and sheds light on its role in epigenetic control of SE activity for cell fate determination.