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

Bronchial epithelium epithelial-mesenchymal plasticity forms aberrant basaloid-like cells in vitro

Although epithelial-mesenchymal transition (EMT) is a common feature of fibrotic lung disease, its role in fibrogenesis is controversial. Recently, aberrant basaloid cells were identified in fibrotic lung tissue as a novel epithelial cell type displaying a partial EMT phenotype. The developmental origin of these cells remains unknown. To elucidate the role of EMT in the development of aberrant basaloid cells from the bronchial epithelium, we mapped EMT-induced transcriptional changes at the population and single-cell levels. Human bronchial epithelial cells grown as submerged or air-liquid interface (ALI) cultures with or without EMT induction were analyzed by bulk and single-cell RNA-Sequencing. Comparison of submerged and ALI cultures revealed differential expression of 8,247 protein coding (PC) and 1,621 long noncoding RNA (lncRNA) genes and revealed epithelial cell-type-specific lncRNAs. Similarly, EMT induction in ALI cultures resulted in robust transcriptional reprogramming of 6,020 PC and 907 lncRNA genes. Although there was no evidence for fibroblast/myofibroblast conversion following EMT induction, cells displayed a partial EMT gene signature and an aberrant basaloid-like cell phenotype. The substantial transcriptional differences between submerged and ALI cultures highlight that care must be taken when interpreting data from submerged cultures. This work supports that lung epithelial EMT does not generate fibroblasts/myofibroblasts and confirms ALI cultures provide a physiologically relevant system to study aberrant basaloid-like cells and mechanisms of EMT. We provide a catalog of PC and lncRNA genes and an interactive browser (https://bronc-epi-in-vitro.cells.ucsc.edu/) of single-cell RNA-Seq data for further exploration of potential roles in the lung epithelium in health and lung disease.

Role of Respiratory Epithelial Cells in Allergic Diseases

The airway epithelium provides the first line of defense to the surrounding environment. However, dysfunctions of this physical barrier are frequently observed in allergic diseases, which are tightly connected with pro- or anti-inflammatory processes. When the epithelial cells are confronted with allergens or pathogens, specific response mechanisms are set in motion, which in homeostasis, lead to the elimination of the invaders and leave permanent traces on the respiratory epithelium. However, allergens can also cause damage in the sensitized organism, which can be ascribed to the excessive immune reactions. The tight interaction of epithelial cells of the upper and lower airways with local and systemic immune cells can leave an imprint that may mirror the pathophysiology. The interaction with effector T cells, along with the macrophages, play an important role in this response, as reflected in the gene expression profiles (transcriptomes) of the epithelial cells, as well as in the secretory pattern (secretomes). Further, the storage of information from past exposures as memories within discrete cell types may allow a tissue to inform and fundamentally alter its future responses. Recently, several lines of evidence have highlighted the contributions from myeloid cells, lymphoid cells, stromal cells, mast cells, and epithelial cells to the emerging concepts of inflammatory memory and trained immunity.

Interactive single-cell data analysis using Cellar

Cell type assignment is a major challenge for all types of high throughput single cell data. In many cases such assignment requires the repeated manual use of external and complementary data sources. To improve the ability to uniformly assign cell types across large consortia, platforms and modalities, we developed Cellar, a software tool that provides interactive support to all the different steps involved in the assignment and dataset comparison process. We discuss the different methods implemented by Cellar, how these can be used with different data types, how to combine complementary data types and how to analyze and visualize spatial data. We demonstrate the advantages of Cellar by using it to annotate several HuBMAP datasets from multi-omics single-cell sequencing and spatial proteomics studies. Cellar is open-source and includes several annotated HuBMAP datasets.

Secretory Immunoglobulin A Immunity in Chronic Obstructive Respiratory Diseases

Chronic obstructive pulmonary disease (COPD), asthma and cystic fibrosis (CF) are distinct respiratory diseases that share features such as the obstruction of small airways and disease flare-ups that are called exacerbations and are often caused by infections. Along the airway epithelium, immunoglobulin (Ig) A contributes to first line mucosal protection against inhaled particles and pathogens. Dimeric IgA produced by mucosal plasma cells is transported towards the apical pole of airway epithelial cells by the polymeric Ig receptor (pIgR), where it is released as secretory IgA. Secretory IgA mediates immune exclusion and promotes the clearance of pathogens from the airway surface by inhibiting their adherence to the epithelium. In this review, we summarize the current knowledge regarding alterations of the IgA/pIgR system observed in those major obstructive airway diseases and discuss their implication for disease pathogenesis.

Inflammatory blockade prevents injury to the developing pulmonary gas exchange surface in preterm primates

Perinatal inflammatory stress is associated with early life morbidity and lifelong consequences for pulmonary health. Chorioamnionitis, an inflammatory condition affecting the placenta and fluid surrounding the developing fetus, affects 25 to 40% of preterm births. Severe chorioamnionitis with preterm birth is associated with significantly increased risk of pulmonary disease and secondary infections in childhood, suggesting that fetal inflammation may markedly alter the development of the lung. Here, we used intra-amniotic lipopolysaccharide (LPS) challenge to induce experimental chorioamnionitis in a prenatal rhesus macaque (Macaca mulatta) model that mirrors structural and temporal aspects of human lung development. Inflammatory injury directly disrupted the developing gas exchange surface of the primate lung, with extensive damage to alveolar structure, particularly the close association and coordinated differentiation of alveolar type 1 pneumocytes and specialized alveolar capillary endothelium. Single-cell RNA sequencing analysis defined a multicellular alveolar signaling niche driving alveologenesis that was extensively disrupted by perinatal inflammation, leading to a loss of gas exchange surface and alveolar simplification, with notable resemblance to chronic lung disease in newborns. Blockade of the inflammatory cytokines interleukin-1β and tumor necrosis factor-α ameliorated LPS-induced inflammatory lung injury by blunting stromal responses to inflammation and modulating innate immune activation in myeloid cells, restoring structural integrity and key signaling networks in the developing alveolus. These data provide new insight into the pathophysiology of developmental lung injury and suggest that modulating inflammation is a promising therapeutic approach to prevent fetal consequences of chorioamnionitis.

A single-cell regulatory map of postnatal lung alveologenesis in humans and mice

Ex-utero regulation of the lungs' responses to breathing air and continued alveolar development shape adult respiratory health. Applying single-cell transposome hypersensitive site sequencing (scTHS-seq) to over 80,000 cells, we assembled the first regulatory atlas of postnatal human and mouse lung alveolar development. We defined regulatory modules and elucidated new mechanistic insights directing alveolar septation, including alveolar type 1 and myofibroblast cell signaling and differentiation, and a unique human matrix fibroblast population. Incorporating GWAS, we mapped lung function causal variants to myofibroblasts and identified a pathogenic regulatory unit linked to lineage marker FGF18, demonstrating the utility of chromatin accessibility data to uncover disease mechanism targets. Our regulatory map and analysis model provide valuable new resources to investigate age-dependent and species-specific control of critical developmental processes. Furthermore, these resources complement existing atlas efforts to advance our understanding of lung health and disease across the human lifespan.

Temporal modelling using single-cell transcriptomics

Methods for profiling genes at the single-cell level have revolutionized our ability to study several biological processes and systems including development, differentiation, response programmes and disease progression. In many of these studies, cells are profiled over time in order to infer dynamic changes in cell states and types, sets of expressed genes, active pathways and key regulators. However, time-series single-cell RNA sequencing (scRNA-seq) also raises several new analysis and modelling issues. These issues range from determining when and how deep to profile cells, linking cells within and between time points, learning continuous trajectories, and integrating bulk and single-cell data for reconstructing models of dynamic networks. In this Review, we discuss several approaches for the analysis and modelling of time-series scRNA-seq, highlighting their steps, key assumptions, and the types of data and biological questions they are most appropriate for.

Clustering spatial transcriptomics data

MOTIVATION

Recent advancements in fluorescence in situ hybridization (FISH) techniques enable them to concurrently obtain information on the location and gene expression of single cells. A key question in the initial analysis of such spatial transcriptomics data is the assignment of cell types. To date, most studies used methods that only rely on the expression levels of the genes in each cell for such assignments. To fully utilize the data and to improve the ability to identify novel sub-types, we developed a new method, FICT, which combines both expression and neighborhood information when assigning cell types.

RESULTS

FICT optimizes a probabilistic function that we formalize and for which we provide learning and inference algorithms. We used FICT to analyze both simulated and several real spatial transcriptomics data. As we show, FICT can accurately identify cell types and sub-types, improving on expression only methods and other methods proposed for clustering spatial transcriptomics data. Some of the spatial sub-types identified by FICT provide novel hypotheses about the new functions for excitatory and inhibitory neurons.

AVAILABILITY AND IMPLEMENTATION

FICT is available at: https://github.com/haotianteng/FICT.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The discovAIR project: a roadmap towards the Human Lung Cell Atlas

The Human Cell Atlas (HCA) consortium aims to establish an atlas of all organs in the healthy human body at single-cell resolution to increase our understanding of basic biological processes that govern development, physiology and anatomy, and to accelerate diagnosis and treatment of disease. The Lung Biological Network of the HCA aims to generate the Human Lung Cell Atlas as a reference for the cellular repertoire, molecular cell states and phenotypes, and cell–cell interactions that characterise normal lung homeostasis in healthy lung tissue. Such a reference atlas of the healthy human lung will facilitate mapping the changes in the cellular landscape in disease. The discovAIR project is one of six pilot actions for the HCA funded by the European Commission in the context of the H2020 framework programme. discovAIR aims to establish the first draft of an integrated Human Lung Cell Atlas, combining single-cell transcriptional and epigenetic profiling with spatially resolving techniques on matched tissue samples, as well as including a number of chronic and infectious diseases of the lung. The integrated Human Lung Cell Atlas will be available as a resource for the wider respiratory community, including basic and translational scientists, clinical medicine, and the private sector, as well as for patients with lung disease and the interested lay public. We anticipate that the Human Lung Cell Atlas will be the founding stone for a more detailed understanding of the pathogenesis of lung diseases, guiding the design of novel diagnostics and preventive or curative interventions.

The discovAIR project contributes to the Human Cell Atlas Lung Biological Network by establishing a first draft of the Human Lung Cell Atlas, advancing our insight into the cellular complexity and spatial organisation of the lung in health and diseasehttps://bit.ly/3zX4cad

Refractory neutrophilic asthma and ciliary genes

BACKGROUND

Refractory asthma (RA) remains poorly controlled, resulting in high health care utilization despite guideline-based therapies. Patients with RA manifest higher neutrophilia as a result of increased airway inflammation and subclinical infection, the underlying mechanisms of which remain unclear.

OBJECTIVE

We sought to characterize and clinically correlate gene expression differences between refractory and nonrefractory (NR) asthma to uncover molecular mechanisms driving group distinctions.

METHODS

Microarray gene expression of paired airway epithelial brush and endobronchial biopsy samples was compared between 60 RA and 30 NR subjects. Subjects were hierarchically clustered to identify subgroups of RA, and biochemical and clinical traits (airway inflammatory molecules, respiratory pathogens, chest imaging) were compared between groups. Weighted gene correlation network analysis was used to identify coexpressed gene modules. Module expression scores were compared between groups using linear regression, controlling for age, sex, and body mass index.

RESULTS

Differential gene expression analysis showed upregulation of proneutrophilic and downregulation of ciliary function genes/pathways in RA compared to NR. A subgroup of RA with downregulated ciliary gene expression had increased levels of subclinical infections, airway neutrophilia, and eosinophilia as well as higher chest imaging mucus burden compared to other RA, the dominant differences between RA and NR. Weighted gene correlation network analysis identified gene modules related to ciliary function, which were downregulated in RA and were associated with lower pulmonary function and higher airway wall thickness/inflammation, markers of poorer asthma control.

CONCLUSIONS

Identification of a novel ciliary-deficient subgroup of RA suggests that diminished mucociliary clearance may underlie repeated asthma exacerbations despite adequate treatment, necessitating further exploration of function, mechanism, and therapeutics.

Differential LysoTracker Uptake Defines Two Populations of Distal Epithelial Cells in Idiopathic Pulmonary Fibrosis

Idiopathic pulmonary fibrosis (IPF) is a progressive and fatal degenerative lung disease of unknown etiology. Although in its final stages it implicates, in a reactive manner, all lung cell types, the initial damage involves the alveolar epithelial compartment, in particular the alveolar epithelial type 2 cells (AEC2s). AEC2s serve dual progenitor and surfactant secreting functions, both of which are deeply impacted in IPF. Thus, we hypothesize that the size of the surfactant processing compartment, as measured by LysoTracker incorporation, allows the identification of different epithelial states in the IPF lung. Flow cytometry analysis of epithelial LysoTracker incorporation delineates two populations (Lysohigh and Lysolow) of AEC2s that behave in a compensatory manner during bleomycin injury and in the donor/IPF lung. Employing flow cytometry and transcriptomic analysis of cells isolated from donor and IPF lungs, we demonstrate that the Lysohigh population expresses all classical AEC2 markers and is drastically diminished in IPF. The Lysolow population, which is increased in proportion in IPF, co-expressed AEC2 and basal cell markers, resembling the phenotype of the previously identified intermediate AEC2 population in the IPF lung. In that regard, we provide an in-depth flow-cytometry characterization of LysoTracker uptake, HTII-280, proSP-C, mature SP-B, NGFR, KRT5, and CD24 expression in human lung epithelial cells. Combining functional analysis with extracellular and intracellular marker expression and transcriptomic analysis, we advance the current understanding of epithelial cell behavior and fate in lung fibrosis.

Deciphering the Immune-Tumor Interplay During Early-Stage Lung Cancer Development via Single-Cell Technology

Lung cancer is the leading cause of cancer-related death worldwide. Cancer immunotherapy has shown great success in treating advanced-stage lung cancer but has yet been used to treat early-stage lung cancer, mostly due to lack of understanding of the tumor immune microenvironment in early-stage lung cancer. The immune system could both constrain and promote tumorigenesis in a process termed immune editing that can be divided into three phases, namely, elimination, equilibrium, and escape. Current understanding of the immune response toward tumor is mainly on the "escape" phase when the tumor is clinically detectable. The detailed mechanism by which tumor progenitor lesions was modulated by the immune system during early stage of lung cancer development remains elusive. The advent of single-cell sequencing technology enables tumor immunologists to address those fundamental questions. In this perspective, we will summarize our current understanding and big gaps about the immune response during early lung tumorigenesis. We will then present the state of the art of single-cell technology and then envision how single-cell technology could be used to address those questions. Advances in the understanding of the immune response and its dynamics during malignant transformation of pre-malignant lesion will shed light on how malignant cells interact with the immune system and evolve under immune selection. Such knowledge could then contribute to the development of precision and early intervention strategies toward lung malignancy.

Lack of evidence for increased transcriptional noise in aged tissues

Aging is often associated with a loss of cell type identity that results in an increase in transcriptional noise in aged tissues. If this phenomenon reflects a fundamental property of aging remains an open question. Transcriptional changes at the cellular level are best detected by single-cell RNA sequencing (scRNAseq). However, the diverse computational methods used for the quantification of age-related loss of cellular identity have prevented reaching meaningful conclusions by direct comparison of existing scRNAseq datasets. To address these issues we created Decibel, a Python toolkit that implements side-to-side four commonly used methods for the quantification of age-related transcriptional noise in scRNAseq data. Additionally, we developed Scallop, a novel computational method for the quantification of membership of single cells to their assigned cell type cluster. Cells with a greater Scallop membership score are transcriptionally more stable. Application of these computational tools to seven aging datasets showed large variability between tissues and datasets, suggesting that increased transcriptional noise is not a universal hallmark of aging. To understand the source of apparent loss of cell type identity associated with aging, we analyzed cell type-specific changes in transcriptional noise and the changes in cell type composition of the mammalian lung. No robust pattern of cell type-specific transcriptional noise alteration was found across aging lung datasets. In contrast, age-associated changes in cell type composition of the lung were consistently found, particularly of immune cells. These results suggest that claims of increased transcriptional noise of aged tissues should be reformulated.

Respiratory manifestations in the Ehlers-Danlos syndromes

Persons with the Ehlers-Danlos syndromes (EDS) report a wide range of respiratory symptoms, most commonly shortness of breath, exercise limitation, and cough. Also reported are noisy breathing attributed to asthma, difficulty with deep inhalation, and inspiratory thoracic pain. The literature consists of case reports and small cross-sectional and cohort studies. One case-control study estimated twofold to threefold greater respiratory disease burden among persons with EDS as compared to controls. The differential diagnosis for symptoms is broad. Structural alterations include pectus deformities, scoliosis, recurrent rib subluxations, and tracheobronchomalacia, associated with varying degrees of physiologic impairment. Those with vascular EDS have an increased risk of pneumothorax, intrapulmonary bleeding, cysts, and nonmalignant fibrous nodules. Functional aerodigestive manifestations such as inducible laryngeal obstruction may be misdiagnosed as asthma, with gastro-esophageal dysmotility and reflux as common contributing factors. Inflammatory manifestations include costochondritis, bronchiectasis, and localized respiratory allergic and nonallergic mast cell activation. Cranio-cervical instability can dysregulate respiratory control pathways. There is a need for careful phenotyping using standardized clinical tools and patient-reported outcomes and continuing collaboration with aerodigestive specialists including otolaryngologists and gastroenterologists. Also needed is further evaluation of respiratory symptoms in persons with hypermobility spectrum disorders. Personalized monitoring strategies are invaluable for interpretation and long-term management of respiratory symptoms.

Large-scale integration of single-cell transcriptomic data captures transitional progenitor states in mouse skeletal muscle regeneration

Skeletal muscle repair is driven by the coordinated self-renewal and fusion of myogenic stem and progenitor cells. Single-cell gene expression analyses of myogenesis have been hampered by the poor sampling of rare and transient cell states that are critical for muscle repair, and do not inform the spatial context that is important for myogenic differentiation. Here, we demonstrate how large-scale integration of single-cell and spatial transcriptomic data can overcome these limitations. We created a single-cell transcriptomic dataset of mouse skeletal muscle by integration, consensus annotation, and analysis of 23 newly collected scRNAseq datasets and 88 publicly available single-cell (scRNAseq) and single-nucleus (snRNAseq) RNA-sequencing datasets. The resulting dataset includes more than 365,000 cells and spans a wide range of ages, injury, and repair conditions. Together, these data enabled identification of the predominant cell types in skeletal muscle, and resolved cell subtypes, including endothelial subtypes distinguished by vessel-type of origin, fibro-adipogenic progenitors defined by functional roles, and many distinct immune populations. The representation of different experimental conditions and the depth of transcriptome coverage enabled robust profiling of sparsely expressed genes. We built a densely sampled transcriptomic model of myogenesis, from stem cell quiescence to myofiber maturation, and identified rare, transitional states of progenitor commitment and fusion that are poorly represented in individual datasets. We performed spatial RNA sequencing of mouse muscle at three time points after injury and used the integrated dataset as a reference to achieve a high-resolution, local deconvolution of cell subtypes. We also used the integrated dataset to explore ligand-receptor co-expression patterns and identify dynamic cell-cell interactions in muscle injury response. We provide a public web tool to enable interactive exploration and visualization of the data. Our work supports the utility of large-scale integration of single-cell transcriptomic data as a tool for biological discovery.

Targeting monocytes/macrophages in fibrosis and cancer diseases: Therapeutic approaches

Over almost 140 years since their identification, the knowledge about macrophages has unbelievably evolved. The 'big eaters' from being thought of as simple phagocytic cells have been recognized as master regulators in immunity, homeostasis, healing/repair and organ development. Long considered to originate exclusively from bone marrow-derived circulating monocytes, macrophages have been also demonstrated to be the first immune cells colonizing tissues in the developing embryo and persisting in adult life by self-renewal, as long-lived tissue resident macrophages. Therefore, heterogeneous populations of macrophages with different ontogeny and functions co-exist in tissues. Macrophages act as sentinels of homeostasis and are intrinsically programmed to lead the wound healing and repair processes that occur after injury. However, in certain pathological circumstances macrophages get dysfunctional, and impaired or aberrant macrophage activities become key features of diseases. For instance, in both fibrosis and cancer, that have been defined 'wounds that do not heal', dysfunctional monocyte-derived macrophages overall play a key detrimental role. On the other hand, due to their plasticity these cells can be 're-educated' and exert anti-fibrotic and anti-cancer functions. Therefore macrophages represent an important therapeutic target in both fibrosis and cancer diseases. The current review will illustrate new insights into the role of monocytes/macrophages in these devastating diseases and summarize new therapeutic strategies and applications of macrophage-targeted drug development in their clinical setting.

The cellular composition of the lung lining fluid gradually changes from bronchus to alveolus

Although large advances have recently been made mapping out the cellular composition of lung tissue using single cell sequencing, the composition and distribution of the cellular elements within the lining fluid of the lung has not been extensively studied. Here, we assessed the cellular composition of the lung lining fluid by performing a differential cell analysis on bronchoalveolar lavage fluid (BALF) and epithelial lining fluid (ELF) at four different locations within the lung in post-lung transplantation patients. The percentage of neutrophils and lymphocytes is reduced in more distal regions of the lungs, while the percentage of macrophages increases in these more distal regions. These data provide valuable information to determine which lung lining fluid sampling technique and location is best to use for measuring specific factors and biomarkers, and to increase the understanding of different cell populations in specific lung regions.

A machine learning method for the discovery of minimum marker gene combinations for cell type identification from single-cell RNA sequencing

Single-cell genomics is rapidly advancing our knowledge of the diversity of cell phenotypes, including both cell types and cell states. Driven by single-cell/-nucleus RNA sequencing (scRNA-seq), comprehensive cell atlas projects characterizing a wide range of organisms and tissues are currently underway. As a result, it is critical that the transcriptional phenotypes discovered are defined and disseminated in a consistent and concise manner. Molecular biomarkers have historically played an important role in biological research, from defining immune cell types by surface protein expression to defining diseases by their molecular drivers. Here, we describe a machine learning-based marker gene selection algorithm, NS-Forest version 2.0, which leverages the nonlinear attributes of random forest feature selection and a binary expression scoring approach to discover the minimal marker gene expression combinations that optimally capture the cell type identity represented in complete scRNA-seq transcriptional profiles. The marker genes selected provide an expression barcode that serves as both a useful tool for downstream biological investigation and the necessary and sufficient characteristics for semantic cell type definition. The use of NS-Forest to identify marker genes for human brain middle temporal gyrus cell types reveals the importance of cell signaling and noncoding RNAs in neuronal cell type identity.

Differences and similarities between the upper and lower airway: focusing on innate immunity

The nose is the first respiratory barrier to external pathogens, allergens, pollutants, or cigarette smoke, and vigorous immune responses are triggered when external pathogens come in contact with the nasal epithelium. The mucosal epithelial cells of the nose are essential to the innate immune response against external pathogens and transmit signals that modulate the adaptive immune response. The upper and lower airways share many physiological and immunological features, but there are also numerous differences. It is crucial to understand these differences and their contribution to pathophysiology in order to optimize treatments for inflammatory diseases of the respiratory tract. This review summarizes important differences in the embryological development, histological features, microbiota, immune responses, and cellular subtypes of mucosal epithelial cells of the nose and lungs.

Diagnostic Classification of Interstitial Lung Disease in Clinical Practice

Interstitial lung diseases (ILDs) are challenging to diagnose, requiring integration of multiple complex features that are often difficult to interpret. This article reviews a pragmatic approach to ILD diagnosis and classification, focusing on diagnostic tools and strategies that are used to separate different subtypes and identify the most appropriate management. We discuss the evolution of ILD classification and the contemporary approach that integrates routinely used diagnostic tools in a multidisciplinary discussion. We highlight the increasing importance of taking a multipronged approach to ILD classification that reflects the recent emphasis on disease behavior while also considering etiopathogenesis and morphologic features.

COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure^1-4, but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63⁺ intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

National Heart, Lung, and Blood Institute and Building Respiratory Epithelium and Tissue for Health (BREATH) Consortium Workshop Report: Moving Forward in Lung Regeneration

The National Heart, Lung, and Blood Institute of the National Institutes of Health, together with the Longfonds BREATH consortium, convened a working group to review the field of lung regeneration and suggest avenues for future research. The meeting took place on May 22, 2019, at the American Thoracic Society 2019 conference in Dallas, Texas, United States, and brought together investigators studying lung development, adult stem-cell biology, induced pluripotent stem cells, biomaterials, and respiratory disease. The purpose of the working group was 1) to examine the present status of basic science approaches to tackling lung disease and promoting lung regeneration in patients and 2) to determine priorities for future research in the field.

Cross-species integration of single-cell RNA-seq resolved alveolar-epithelial transitional states in idiopathic pulmonary fibrosis

Single-cell transcriptomics analyses of the fibrotic lung uncovered two cell states critical to lung injury recovery in the alveolar epithelium-a reparative transitional cell state in the mouse and a disease-specific cell state (KRT5^-/KRT17⁺) in human idiopathic pulmonary fibrosis (IPF). The murine transitional cell state lies between the differentiation from type 2 (AT2) to type 1 pneumocyte (AT1), and the human KRT5^-/KRT17⁺ cell state may arise from the dysregulation of this differentiation process. We review major findings of single-cell transcriptomics analyses of the fibrotic lung and reanalyzed data from seven single-cell RNA sequencing studies of human and murine models of IPF, focusing on the alveolar epithelium. Our comparative and cross-species single-cell transcriptomics analyses allowed us to further delineate the differentiation trajectories from AT2 to AT1 and AT2 to the KRT5^-/KRT17⁺ cell state. We observed AT1 cells in human IPF retain the transcriptional signature of the murine transitional cell state. Using pseudotime analysis, we recapitulated the differentiation trajectories from AT2 to AT1 and from AT2 to KRT5^-/KRT17⁺ cell state in multiple human IPF studies. We further delineated transcriptional programs underlying cell-state transitions and determined the molecular phenotypes at terminal differentiation. We hypothesize that in addition to the reactivation of developmental programs (SOX4, SOX9), senescence (TP63, SOX4) and the Notch pathway (HES1) are predicted to steer intermediate progenitors to the KRT5^-/KRT17⁺ cell state. Our analyses suggest that activation of SMAD3 later in the differentiation process may explain the fibrotic molecular phenotype typical of KRT5^-/KRT17⁺ cells.

The Known Unknowns of the Immune Response to Coccidioides

Coccidioidomycosis, otherwise known as Valley Fever, is caused by the dimorphic fungi Coccidioides immitis and C. posadasii. While most clinical cases present with self-limiting pulmonary infection, dissemination of Coccidioides spp. results in prolonged treatment and portends higher mortality rates. While the structure, genome, and niches for Coccidioides have provided some insight into the pathogenesis of disease, the underlying immunological mechanisms of clearance or inability to contain the infection in the lung are poorly understood. This review focuses on the known innate and adaptive immune responses to Coccidioides and highlights three important areas of uncertainty and potential approaches to address them. Closing these gaps in knowledge may enable new preventative and therapeutic strategies to be pursued.

Hydroxypropyl Cyclodextrin Improves Amiodarone-induced Aberrant Lipid Homeostasis of Alveolar Cells

Alveolar epithelial type II (AT2) cells secrete pulmonary surfactant via lamellar bodies (LBs). Abnormalities in LBs are associated with pulmonary disorders, including fibrosis. However, high-content screening (HCS) for LB abnormalities is limited by the lack of understanding of AT2 cell functions. In the present study, we have developed LB cells harboring LB-like organelles that secrete surfactant proteins. These cells were more similar to AT2 cells than to parental A549 cells. LB cells recapitulated amiodarone (AMD)-induced LB enlargement, similar to AT2 cells of patients exposed to AMD. To reverse AMD-induced LB abnormalities, we performed HCS of approved drugs and identified 2-hydroxypropyl-β-cyclodextrin (HPβCD), a cyclic oligosaccharide, as a potential therapeutic agent. A transcriptome analysis revealed that HPβCD modulates lipid homeostasis. In addition, HPβCD inhibited AMD-induced LB abnormalities in human induced pluripotent stem cell-derived AT2 cells. Our results demonstrate that LB cells are useful for HCS and suggest that HPβCD is a candidate therapeutic agent for AMD-induced interstitial pneumonia.

Airway organoids as models of human disease

Studies developing and applying organoid technology have greatly increased in volume and visibility over the past decade. Organoids are three-dimensional structures that are established from pluripotent stem cells (PSCs) or adult tissue stem cells (ASCs). They consist of organ-specific cell types that self-organize through cell sorting and spatially restricted lineage commitment to generate architectural and functional characteristics of the tissue of interest. The field of respiratory development and disease has been particularly productive in this regard. Starting from human cells (PSCs or ASCs), models of the two segments of the lung, the airways and the alveoli, can be built. Such organoids allow the study of development, physiology and disease and thus bridge the gap between animal models and clinical studies. This review discusses current developments in the pulmonary organoid field, highlighting the potential and limitations of current models.

A high-resolution cell atlas of the domestic pig lung and an online platform for exploring lung single-cell data

The genetically engineered pig is regarded as an optimal source of organ transplantation for humans and an excellent model for human disease research, given its comparable physiology to human beings. A myriad of single-cell RNA sequencing (scRNA-seq) data on humans has been reported, but such data on pigs are scarce. Here, we apply scRNA-seq technology to study the cellular heterogeneity of 3-month-old pig lungs, generating the single-cell atlas of 13,580 cells covering 16 major cell types. Based on these data, we systematically characterize the similarities and differences in the cellular cross-talk and expression patterns of respiratory virus receptors in each cell type of pig lungs compared with human lungs. Furthermore, we analyze pig lung xenotransplantation barriers and reported the cell-type expression patterns of 10 genes associated with pig-to-human immunobiological incompatibility and coagulation dysregulation. We also investigate the conserved transcription factors (TFs) and their candidate target genes and constructed five conserved TF regulatory networks in the main cell types shared by pig and human lungs. Finally, we present a comprehensive and openly accessible online platform, ScdbLung. Our scRNA-seq atlas of the domestic pig lung and ScdbLung database can guide pig lung research and clinical applicability.

Statistical mechanics meets single-cell biology

Single-cell omics is transforming our understanding of cell biology and disease, yet the systems-level analysis and interpretation of single-cell data faces many challenges. In this Perspective, we describe the impact that fundamental concepts from statistical mechanics, notably entropy, stochastic processes and critical phenomena, are having on single-cell data analysis. We further advocate the need for more bottom-up modelling of single-cell data and to embrace a statistical mechanics analysis paradigm to help attain a deeper understanding of single-cell systems biology.

The aging lung: Physiology, disease, and immunity

The population is aging at a rate never seen before in human history. As the number of elderly adults grows, it is imperative we expand our understanding of the underpinnings of aging biology. Human lungs are composed of a unique panoply of cell types that face ongoing chemical, mechanical, biological, immunological, and xenobiotic stress over a lifetime. Yet, we do not fully appreciate the mechanistic drivers of lung aging and why age increases the risk of parenchymal lung disease, fatal respiratory infection, and primary lung cancer. Here, we review the molecular and cellular aspects of lung aging, local stress response pathways, and how the aging process predisposes to the pathogenesis of pulmonary disease. We place these insights into context of the COVID-19 pandemic and discuss how innate and adaptive immunity within the lung is altered with age.

Advanced single-cell technologies to guide the development of bioengineered lungs

Lung transplantation remains the only viable option for individuals suffering from end-stage lung failure. However, a number of current limitations exist including a continuing shortage of suitable donor lungs and immune rejection following transplantation. To address these concerns, engineering a decellularized biocompatible lung scaffold from cadavers reseeded with autologous lung cells to promote tissue regeneration is being explored. Proof-of-concept transplantation of these bioengineered lungs into animal models has been accomplished. However, these lungs were incompletely recellularized with resulting epithelial and endothelial leakage and insufficient basement membrane integrity. Failure to repopulate lung scaffolds with all of the distinct cell populations necessary for proper function remains a significant hurdle for the progression of current engineering approaches and precludes clinical translation. Advancements in 3D bioprinting, lung organoid models, and microfluidic device and bioreactor development have enhanced our knowledge of pulmonary lung development, as well as important cell-cell and cell-matrix interactions, all of which will help in the path to a bioengineered transplantable lung. However, a significant gap in knowledge of the spatiotemporal interactions between cell populations as well as relative quantities and localization within each compartment of the lung necessary for its proper growth and function remains. This review will provide an update on cells currently used for reseeding decellularized scaffolds with outcomes of recent lung engineering attempts. Focus will then be on how data obtained from advanced single-cell analyses, coupled with multiomics approaches and high-resolution 3D imaging, can guide current lung bioengineering efforts for the development of fully functional, transplantable lungs.

Integrative analysis of cell state changes in lung fibrosis with peripheral protein biomarkers

The correspondence of cell state changes in diseased organs to peripheral protein signatures is currently unknown. Here, we generated and integrated single-cell transcriptomic and proteomic data from multiple large pulmonary fibrosis patient cohorts. Integration of 233,638 single-cell transcriptomes (n = 61) across three independent cohorts enabled us to derive shifts in cell type proportions and a robust core set of genes altered in lung fibrosis for 45 cell types. Mass spectrometry analysis of lung lavage fluid (n = 124) and plasma (n = 141) proteomes identified distinct protein signatures correlated with diagnosis, lung function, and injury status. A novel SSTR2+ pericyte state correlated with disease severity and was reflected in lavage fluid by increased levels of the complement regulatory factor CFHR1. We further discovered CRTAC1 as a biomarker of alveolar type-2 epithelial cell health status in lavage fluid and plasma. Using cross-modal analysis and machine learning, we identified the cellular source of biomarkers and demonstrated that information transfer between modalities correctly predicts disease status, suggesting feasibility of clinical cell state monitoring through longitudinal sampling of body fluid proteomes.

Overlapping and distinct features of viral and allergen immunity in the human lung

Immunity in the human respiratory tract is provided by a diverse range of tissue-resident cells, including specialized epithelial and macrophage populations and a network of innate and innate-like lymphocytes, such as natural killer cells, innate lymphoid cells, and invariant T cells. Lung-resident memory T and B cells contribute to this network following initial exposure to antigenic stimuli. This review explores how advances in the study of human immunology have shaped our understanding of this resident immune network and its response to two of the most commonly encountered inflammatory stimuli in the airways: viruses and allergens. It discusses the many ways in which pathogenic infection and allergic inflammation mirror each other, highlighting the key checkpoints at which they diverge and how this can result in a lifetime of allergic exacerbation versus protective anti-viral immunity.

Versatile workflow for cell type-resolved transcriptional and epigenetic profiles from cryopreserved human lung

Complexity of lung microenvironment and changes in cellular composition during disease make it exceptionally hard to understand molecular mechanisms driving development of chronic lung diseases. Although recent advances in cell type-resolved approaches hold great promise for studying complex diseases, their implementation relies on local access to fresh tissue, as traditional tissue storage methods do not allow viable cell isolation. To overcome these hurdles, we developed a versatile workflow that allows storage of lung tissue with high viability, permits thorough sample quality check before cell isolation, and befits sequencing-based profiling. We demonstrate that cryopreservation enables isolation of multiple cell types from both healthy and diseased lungs. Basal cells from cryopreserved airways retain their differentiation ability, indicating that cellular identity is not altered by cryopreservation. Importantly, using RNA sequencing and EPIC Array, we show that gene expression and DNA methylation signatures are preserved upon cryopreservation, emphasizing the suitability of our workflow for omics profiling of lung cells. Moreover, we obtained high-quality single-cell RNA-sequencing data of cells from cryopreserved human lungs, demonstrating that cryopreservation empowers single-cell approaches. Overall, thanks to its simplicity, our workflow is well suited for prospective tissue collection by academic collaborators and biobanks, opening worldwide access to viable human tissue.

Single-cell RNA sequencing of human liver reveals hepatic stellate cell heterogeneity

Background & Aims

The multiple vital functions of the human liver are performed by highly specialised parenchymal and non-parenchymal cells organised in complex collaborative sinusoidal units. Although crucial for homeostasis, the cellular make-up of the human liver remains to be fully elucidated. Here, single-cell RNA-sequencing was used to unravel the heterogeneity of human liver cells, in particular of hepatocytes (HEPs) and hepatic stellate cells (HSCs).

Method

The transcriptome of ~25,000 freshly isolated human liver cells was profiled using droplet-based RNA-sequencing. Recently published data sets and RNA in situ hybridisation were integrated to validate and locate newly identified cell populations.

Results

In total, 22 cell populations were annotated that reflected the heterogeneity of human parenchymal and non-parenchymal liver cells. More than 20,000 HEPs were ordered along the portocentral axis to confirm known, and reveal previously undescribed, zonated liver functions. The existence of 2 subpopulations of human HSCs with unique gene expression signatures and distinct intralobular localisation was revealed (i.e. portal and central vein-concentrated GPC3⁺ HSCs and perisinusoidally located DBH⁺ HSCs). In particular, these data suggest that, although both subpopulations collaborate in the production and organisation of extracellular matrix, GPC3⁺ HSCs specifically express genes involved in the metabolism of glycosaminoglycans, whereas DBH⁺ HSCs display a gene signature that is reminiscent of antigen-presenting cells.

Conclusions

This study highlights metabolic zonation as a key determinant of HEP transcriptomic heterogeneity and, for the first time, outlines the existence of heterogeneous HSC subpopulations in the human liver. These findings call for further research on the functional implications of liver cell heterogeneity in health and disease.

Lay summary

This study resolves the cellular landscape of the human liver in an unbiased manner and at high resolution to provide new insights into human liver cell biology. The results highlight the physiological heterogeneity of human hepatic stellate cells.

•

A cell atlas from the near-native transcriptome of >25,000 human liver cells is presented.

•

Hepatocytes were ordered along the portocentral axis to reveal previously undescribed gene expression patterns and zonated liver functions.

•

Two subpopulations of human hepatic stellate cells (HSCs) are reported, characterised by different spatial distribution in the native tissue.

•

Characteristic gene signatures of HSC subpopulations are suggestive of far-reaching functional differences.

Network Analysis and Transcriptome Profiling Identify Autophagic and Mitochondrial Dysfunctions in SARS-CoV-2 Infection

Analyzing host cells' transcriptional response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection will help delineate biological processes underlying viral pathogenesis. First, analysis of expression profiles of lung cell lines A549 and Calu3 revealed upregulation of antiviral interferon signaling genes in response to all three SARS-CoV-2, MERS-CoV, or influenza A virus (IAV) infections. However, perturbations in expression of genes involved in inflammatory, mitochondrial, and autophagy processes were specifically observed in SARS-CoV-2-infected cells. Next, a validation study in infected human nasopharyngeal samples also revealed perturbations in autophagy and mitochondrial processes. Specifically, mTOR expression, mitochondrial ribosomal, mitochondrial complex I, lysosome acidification, and mitochondrial fission promoting genes were concurrently downregulated in both infected cell lines and human samples. SARS-CoV-2 infection impeded autophagic flux either by upregulating GSK3B in lung cell lines or by downregulating autophagy genes, SNAP29, and lysosome acidification genes in human samples, contributing to increased viral replication. Therefore, drugs targeting lysosome acidification or autophagic flux could be tested as intervention strategies. Finally, age-stratified SARS-CoV-2-positive human data revealed impaired upregulation of chemokines, interferon-stimulated genes, and tripartite motif genes that are critical for antiviral signaling. Together, this analysis has revealed specific aspects of autophagic and mitochondrial function that are uniquely perturbed in SARS-CoV-2-infected host cells.

CiliOPD: a ciliopathy-associated COPD endotype

The pathophysiology of chronic obstructive pulmonary disease (COPD) relies on airway remodelling and inflammation. Alterations of mucociliary clearance are a major hallmark of COPD caused by structural and functional cilia abnormalities. Using transcriptomic databases of whole lung tissues and isolated small airway epithelial cells (SAEC), we comparatively analysed cilia-associated and ciliopathy-associated gene signatures from a set of 495 genes in 7 datasets including 538 non-COPD and 508 COPD patients. This bio-informatics approach unveils yet undescribed cilia and ciliopathy genes associated with COPD including NEK6 and PROM2 that may contribute to the pathology, and suggests a COPD endotype exhibiting ciliopathy features (CiliOPD).

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

Epithelial Cells and Inflammation in Pulmonary Wound Repair

Respiratory diseases are frequently characterised by epithelial injury, airway inflammation, defective tissue repair, and airway remodelling. This may occur in a subacute or chronic context, such as asthma and chronic obstructive pulmonary disease, or occur acutely as in pathogen challenge and acute respiratory distress syndrome (ARDS). Despite the frequent challenge of lung homeostasis, not all pulmonary insults lead to disease. Traditionally thought of as a quiescent organ, emerging evidence highlights that the lung has significant capacity to respond to injury by repairing and replacing damaged cells. This occurs with the appropriate and timely resolution of inflammation and concurrent initiation of tissue repair programmes. Airway epithelial cells are key effectors in lung homeostasis and host defence; continual exposure to pathogens, toxins, and particulate matter challenge homeostasis, requiring robust defence and repair mechanisms. As such, the epithelium is critically involved in the return to homeostasis, orchestrating the resolution of inflammation and initiating tissue repair. This review examines the pivotal role of pulmonary airway epithelial cells in initiating and moderating tissue repair and restitution. We discuss emerging evidence of the interactions between airway epithelial cells and candidate stem or progenitor cells to initiate tissue repair as well as with cells of the innate and adaptive immune systems in driving successful tissue regeneration. Understanding the mechanisms of intercellular communication is rapidly increasing, and a major focus of this review includes the various mediators involved, including growth factors, extracellular vesicles, soluble lipid mediators, cytokines, and chemokines. Understanding these areas will ultimately identify potential cells, mediators, and interactions for therapeutic targeting.

Rho-Kinase 1/2 Inhibition Prevents Transforming Growth Factor-β-Induced Effects on Pulmonary Remodeling and Repair

Transforming growth factor (TGF)-β-induced myofibroblast transformation and alterations in mesenchymal-epithelial interactions contribute to chronic lung diseases such as chronic obstructive pulmonary disease (COPD), asthma and pulmonary fibrosis. Rho-associated coiled-coil-forming protein kinase (ROCK) consists as two isoforms, ROCK1 and ROCK2, and both are playing critical roles in many cellular responses to injury. In this study, we aimed to elucidate the differential role of ROCK isoforms on TGF-β signaling in lung fibrosis and repair. For this purpose, we tested the effect of a non-selective ROCK 1 and 2 inhibitor (compound 31) and a selective ROCK2 inhibitor (compound A11) in inhibiting TGF-β-induced remodeling in lung fibroblasts and slices; and dysfunctional epithelial-progenitor interactions in lung organoids. Here, we demonstrated that the inhibition of ROCK1/2 with compound 31 represses TGF-β-driven actin remodeling as well as extracellular matrix deposition in lung fibroblasts and PCLS, whereas selective ROCK2 inhibition with compound A11 did not. Furthermore, the TGF-β induced inhibition of organoid formation was functionally restored in a concentration-dependent manner by both dual ROCK 1 and 2 inhibition and selective ROCK2 inhibition. We conclude that dual pharmacological inhibition of ROCK 1 and 2 counteracts TGF-β induced effects on remodeling and alveolar epithelial progenitor function, suggesting this to be a promising therapeutic approach for respiratory diseases associated with fibrosis and defective lung repair.

Regulation of Immune Responses by Nonhematopoietic Cells in Asthma

Nonhematopoietic cells are emerging as important contributors to many inflammatory diseases, including allergic asthma. Recent advances have led to a deeper understanding of how these cells interact with traditional immune cells, thereby modulating their activities in both homeostasis and disease. In addition to their well-established roles in gas exchange and barrier function, lung epithelial cells express an armament of innate sensors that can be triggered by various inhaled environmental agents, leading to the production of proinflammatory molecules. Advances in cell lineage tracing and single-cell RNA sequencing have expanded our knowledge of rare, but immunologically important nonhematopoietic cell populations. In parallel with these advances, novel reverse genetic approaches are revealing how individual genes in different lung-resident nonhematopoietic cell populations contribute to the initiation and maintenance of asthma. This knowledge is already revealing new pathways that can be selectively targeted to treat distinct forms of asthma.

Distinct disease severity between children and older adults with COVID-19: Impacts of ACE2 expression, distribution, and lung progenitor cells

BACKGROUND

Children and older adults with coronavirus disease 2019 (COVID-19) display a distinct spectrum of disease severity yet the risk factors aren't well understood. We sought to examine the expression pattern of angiotensin-converting enzyme 2 (ACE2), the cell-entry receptor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and the role of lung progenitor cells in children and older patients.

METHODS

We retrospectively analysed clinical features in a cohort of 299 patients with COVID-19. The expression and distribution of ACE2 and lung progenitor cells were systematically examined using a combination of public single-cell RNA-seq datasets, lung biopsies, and ex vivo infection of lung tissues with SARS-CoV-2 pseudovirus in children and older adults. We also followed up patients who had recovered from COVID-19.

RESULTS

Compared with children, older patients (> 50 yrs.) were more likely to develop into serious pneumonia with reduced lymphocytes and aberrant inflammatory response (p = 0.001). The expression level of ACE2 and lung progenitor cell markers were generally decreased in older patients. Notably, ACE2 positive cells were mainly distributed in the alveolar region, including SFTPC positive cells, but rarely in airway regions in the older adults (p < 0.01). The follow-up of discharged patients revealed a prolonged recovery from pneumonia in the older (p < 0.025).

CONCLUSION

Compared to children, ACE2 positive cells are generally decreased in older adults and mainly presented in the lower pulmonary tract. The lung progenitor cells are also decreased. These risk factors may impact disease severity and recovery from pneumonia caused by SARS-Cov-2 infection in older patients.

Regulation of immune responses by the airway epithelial cell landscape

The community of cells lining our airways plays a collaborative role in the preservation of immune homeostasis in the lung and provides protection from the pathogens and pollutants in the air we breathe. In addition to its structural attributes that provide effective mucociliary clearance of the lower airspace, the airway epithelium is an immunologically active barrier surface that senses changes in the airway environment and interacts with resident and recruited immune cells. Single-cell RNA-sequencing is illuminating the cellular heterogeneity that exists in the airway wall and has identified novel cell populations with unique molecular signatures, trajectories of differentiation and diverse functions in health and disease. In this Review, we discuss how our view of the airway epithelial landscape has evolved with the advent of transcriptomic approaches to cellular phenotyping, with a focus on epithelial interactions with the local neuronal and immune systems.