A Single-Cell Atlas of the Human Healthy Airways

Epithelial stem and progenitor cells of the upper airway

The upper airway acts as a conduit for the passage of air to the respiratory system and is implicated in several chronic diseases. Whilst the cell biology of the distal respiratory system has been described in great detail, less is known about the proximal upper airway. In this review, we describe the relevant anatomy of the upper airway and discuss the literature detailing the identification and roles of the progenitor cells of these regions.

Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches

Introduction

The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing.

Methods

Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors.

Results

Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19.

Discussion

The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

Idiopathic Subglottic Stenosis Is Associated With More Frequent and Abnormal Squamous Metaplasia

OBJECTIVES

Gain insights into the pathophysiology of idiopathic subglottic stenosis (iSGS) by investigating differences in transcriptome of subglottic mucosal tissue between patients with iSGS and controls, and between tracheal and subglottic tissue within patients.

METHODS

RNA sequencing was conducted on biopsied mucosal samples collected from subglottic and tracheal (in-patient control) regions in iSGS patients, and from subglottis in controls. The gene expression differences were validated on a protein level by (1) staining the tissue samples obtained from a second cohort of patients and controls; and (2) in vitro functional assays using primary subglottic epithelial cells from both iSGS patients and healthy donors.

RESULTS

We found 7 upregulated genes in the subglottic region of iSGS patients relative to both the tracheal mucosa and subglottic region of controls. A gene ontology enrichment analysis found that the epithelial cell differentiation and cornification pathways are significant, involving specifically 3 of the genes: involucrin (IVL), small proline rich protein 1B (SPRR1B), and keratin 16 (KRT16). Involvement of these pathways suggests squamous metaplasia of the epithelium. Histological analyses of epithelium in subglottic mucosal biopsies revealed squamous metaplasia in 41% of the samples from iSGS patients and in 25% from controls. Immunohistochemical evaluation of the samples presented with squamous epithelium revealed increased expression of the protein encoded by SPRR1B, hyperproliferative basal cells, shedding of apical layers, and accompanying lesions in iSGS compared to CTRL. Cultured primary subglottic epithelial cells from iSGS patients had higher proliferation rates compared to healthy donors and squamous metaplastic differentiation formed thinner epithelia with increased expression proteins encoded by INV, SPRR1B, and KRT16, suggesting intrinsic dysfunction of basal cells in iSGS.

CONCLUSIONS

Abnormal squamous differentiation of epithelial cells may contribute to the pathogenesis of iSGS. Patients having metaplastic epithelial phenotype may be sensitive to drugs that reverse it to a normal phenotype.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Dictionary learning for integrative, multimodal and scalable single-cell analysis

Mapping single-cell sequencing profiles to comprehensive reference datasets provides a powerful alternative to unsupervised analysis. However, most reference datasets are constructed from single-cell RNA-sequencing data and cannot be used to annotate datasets that do not measure gene expression. Here we introduce 'bridge integration', a method to integrate single-cell datasets across modalities using a multiomic dataset as a molecular bridge. Each cell in the multiomic dataset constitutes an element in a 'dictionary', which is used to reconstruct unimodal datasets and transform them into a shared space. Our procedure accurately integrates transcriptomic data with independent single-cell measurements of chromatin accessibility, histone modifications, DNA methylation and protein levels. Moreover, we demonstrate how dictionary learning can be combined with sketching techniques to improve computational scalability and harmonize 8.6 million human immune cell profiles from sequencing and mass cytometry experiments. Our approach, implemented in version 5 of our Seurat toolkit ( http://www.satijalab.org/seurat ), broadens the utility of single-cell reference datasets and facilitates comparisons across diverse molecular modalities.

Single-cell resolution of human airway epithelial cells exposed to bronchiolitis obliterans-associated chemicals

Bronchiolitis obliterans (BO) is a fibrotic lung disease characterized by progressive luminal narrowing and obliteration of the small airways. In the nontransplant population, inhalation exposure to certain chemicals is associated with BO; however, the mechanisms contributing to disease induction remain poorly understood. This study's objective was to use single-cell RNA sequencing for the identification of transcriptomic signatures common to primary human airway epithelial cells after chemical exposure to BO-associated chemicals-diacetyl or nitrogen mustard-to help explain BO induction. Primary airway epithelial cells were cultured at air-liquid interface and exposed to diacetyl, nitrogen mustard, or control vapors. Cultures were dissociated and sequenced for single-cell RNA. Differential gene expression and functional pathway analyses were compared across exposures. In total, 75,663 single cells were captured and sequenced from all exposure conditions. Unbiased clustering identified 11 discrete phenotypes, including 5 basal, 2 ciliated, and 2 secretory cell clusters. With chemical exposure, the proportion of cells assigned to keratin 5+ basal cells decreased, whereas the proportion of cells aligned to secretory cell clusters increased compared with control exposures. Functional pathway analysis identified interferon signaling and antigen processing/presentation as pathways commonly upregulated after diacetyl or nitrogen mustard exposure in a ciliated cell cluster. Conversely, the response of airway basal cells differed significantly with upregulation of the unfolded protein response in diacetyl-exposed basal cells, not seen in nitrogen mustard-exposed cultures. These new insights provide early identification of airway epithelial signatures common to BO-associated chemical exposures.NEW & NOTEWORTHY Bronchiolitis obliterans (BO) is a devastating fibrotic lung disease of the small airways, or bronchioles. This original manuscript uses single-cell RNA sequencing for identifying common signatures of chemically exposed airway epithelial cells in BO induction. Chemical exposure reduced the proportion of keratin 5+ basal cells while increasing the proportion of keratin 4+ suprabasal cells. Functional pathways contributory to these shifts differed significantly across exposures. These new results highlight similarities and differences in BO induction across exposures.

Combination of CRISPR-Cas9-RNP and Single-Cell RNAseq to Identify Cell State-Specific FOXJ1 Functions in the Human Airway Epithelium

The study of the airway epithelium in vitro is routinely performed using air-liquid culture (ALI) models from nasal or bronchial basal cells. These 3D experimental models allow to follow the regeneration steps of fully differentiated mucociliary epithelium and to study gene function by performing gene invalidation. Recent progress made with CRISPR-based techniques has overcome the experimental difficulty of this approach, by a direct transfection of ribonucleoprotein complexes combining a mix of synthetic small guide RNAs (sgRNAs) and recombinant Cas9. The approach shows more than 95% efficiency and does not require any selection step. A limitation of this approach is that it generates cell populations that contain heterogeneous deletions, which makes the evaluation of invalidation efficiency difficult. We have successfully used Flongle sequencing (Nanopore) to quantify the number of distinct deletions. We describe the use of CRISPR-Cas9 RNP in combination with single-cell RNA sequencing to functionally characterize the impact of gene invalidation in ALI cultures. The complex ecosystem of the airway epithelium, composed of many cell types, makes single-cell approaches particularly relevant to study cell type, or cell state-specific events. This protocol describes the invalidation of FOXJ1 in ALI cultures through the following steps: (1) Establishment of basal cell cultures from nasal turbinates, (2) CRISPR-Cas9 RNP invalidation of FOXJ1, (3) Quantification of FOXJ1 invalidation efficiency by Nanopore sequencing, (4) Dissociation of ALI cultures and single-cell RNAseq, (5) Analysis of single-cell RNAseq data from FOXJ1-invalidated cells.We confirm here that FOXJ1 invalidation impairs the final differentiation step of multiciliated cells and provides a framework to explore other gene functions.

Changes in the Proportion of Each Cell Type After hiPSC-Derived Airway Epithelia Transplantation

No radical treatment is available for the regeneration of dysfunction and defects in airway epithelia. Artificial tracheae made of polypropylene and collagen sponge were used in clinical studies to reconstitute tracheae after resection. For early epithelialization of the luminal surface of the artificial trachea, a model was established, that is, an artificial trachea covered with human-induced pluripotent stem cell-derived airway epithelial cells (hiPSC-AECs) was transplanted into a tracheal defect in an immunodeficient rat. Unlike the cell types of hiPSC-derived cells that are currently used in clinical studies, AECs maintain tissues by proliferation and differentiation of basal cells into various cell types that constitute AECs constantly. Therefore, post-transplantation, the proportion of each cell type, such as ciliated and goblet cells, may change; however, no studies have examined this possibility. In this study, using our hiPSC-AEC-transplanted rat model, we investigated changes in the proportion of each cell type in hiPSC-AECs pre-transplantation and post-transplantation. As a result, the proportion of each cell type changed post-transplantation. The proportion of ciliated, basal, and club cells increased, and the proportion of goblet cells decreased post-transplantation. In addition, the proportion of each cell type in engrafted hiPSC-AECs is more similar to the proportion of each cell type in normal proximal airway tissue than the proportion of each cell type pre-transplantation. The results of this study are useful for the development of therapeutic techniques using hiPSC-AEC transplantation.

Heterogeneity-Preserving Discriminative Feature Selection for Subtype Discovery

The discovery of subtypes is pivotal for disease diagnosis and targeted therapy, considering the diverse responses of different cells or patients to specific treatments. Exploring the heterogeneity within disease or cell states provides insights into disease progression mechanisms and cell differentiation. The advent of high-throughput technologies has enabled the generation and analysis of various molecular data types, such as single-cell RNA-seq, proteomic, and imaging datasets, at large scales. While presenting opportunities for subtype discovery, these datasets pose challenges in finding relevant signatures due to their high dimensionality. Feature selection, a crucial step in the analysis pipeline, involves choosing signatures that reduce the feature size for more efficient downstream computational analysis. Numerous existing methods focus on selecting signatures that differentiate known diseases or cell states, yet they often fall short in identifying features that preserve heterogeneity and reveal subtypes. To identify features that can capture the diversity within each class while also maintaining the discrimination of known disease states, we employed deep metric learning-based feature embedding to conduct a detailed exploration of the statistical properties of features essential in preserving heterogeneity. Our analysis revealed that features with a significant difference in interquartile range (IQR) between classes possess crucial subtype information. Guided by this insight, we developed a robust statistical method, termed PHet (Preserving Heterogeneity) that performs iterative subsampling differential analysis of IQR and Fisher's method between classes, identifying a minimal set of heterogeneity-preserving discriminative features to optimize subtype clustering quality. Validation using public single-cell RNA-seq and microarray datasets showcased PHet's effectiveness in preserving sample heterogeneity while maintaining discrimination of known disease/cell states, surpassing the performance of previous outlier-based methods. Furthermore, analysis of a single-cell RNA-seq dataset from mouse tracheal epithelial cells revealed, through PHet-based features, the presence of two distinct basal cell subtypes undergoing differentiation toward a luminal secretory phenotype. Notably, one of these subtypes exhibited high expression of BPIFA1. Interestingly, previous studies have linked BPIFA1 secretion to the emergence of secretory cells during mucociliary differentiation of airway epithelial cells. PHet successfully pinpointed the basal cell subtype associated with this phenomenon, a distinction that pre-annotated markers and dispersion-based features failed to make due to their admixed feature expression profiles. These findings underscore the potential of our method to deepen our understanding of the mechanisms underlying diseases and cell differentiation and contribute significantly to personalized medicine.

Meta-Analysis of COVID-19 BAL Single-Cell RNA Sequencing Reveals Alveolar Epithelial Transitions and Unique Alveolar Epithelial Cell Fates

Single-cell RNA sequencing (scRNA-seq) of BAL cells has provided insights into coronavirus disease (COVID-19). However, reports have been limited by small patient cohorts. We performed a meta-analysis of BAL scRNA-seq data from healthy control subjects (n = 13) and patients with COVID-19 (n = 20), sourced from six independent studies (167,280 high-quality cells in total). Consistent with the source reports, increases in infiltrating leukocyte subtypes were noted, several with type I IFN signatures and unique gene expression signatures associated with transcellular chemokine signaling. Noting dramatic reductions of inferred NKX2-1 and NR4A1 activity in alveolar epithelial type II (AT-II) cells, we modeled pseudotemporal AT-II-to-AT-I progression. This revealed changes in inferred AT-II cell metabolic activity, increased transitional cells, and a previously undescribed AT-I state. This cell state was conspicuously marked by the induction of genes of the epidermal differentiation complex, including the cornified envelope protein SPRR3 (small proline-rich protein 3), upregulation of multiple KRT (keratin) genes, inferred mitochondrial dysfunction, and cell death signatures including apoptosis and ferroptosis. Immunohistochemistry of lungs from patients with COVID-19 confirmed upregulation and colocalization of KRT13 and SPRR3 in the distal airspaces. Forced overexpression of SPRR3 in human alveolar epithelial cells ex vivo did not activate caspase-3 or upregulate KRT13, suggesting that SPRR3 marks an AT-I cornification program in COVID-19 but is not sufficient for phenotypic changes.

Comparative immunological landscape between pre- and early-stage LUAD manifested as ground-glass nodules revealed by scRNA and scTCR integrated analysis

BACKGROUND

Mechanism underlying the malignant progression of precancer to early-stage lung adenocarcinoma (LUAD) as well as their indolence nature remains elusive.

METHODS

Single-cell RNA sequencing (scRNA) with simultaneous T cell receptor (TCR) sequencing on 5 normal lung tissues, 3 precancerous and 4 early-stage LUAD manifested as pulmonary ground-glass nodules (GGNs) were performed.

RESULTS

Through this integrated analysis, we have delineated five key modules that drive the malignant progression of early-stage LUAD in a disease stage-dependent manner. These modules are related to cell proliferation and metabolism, immune response, mitochondria, cilia, and cell adhesion. We also find that the tumor micro-environment (TME) of early-stage LUAD manifested as GGN are featured with regulatory T (Tregs) cells accumulation with three possible origins, and loss-functional state (decreased clonal expansion and cytotoxicity) of CD8 + T cells. Instead of exhaustion, the CD8 + T cells are featured with a shift to memory phenotype, which is significantly different from the late stage LUAD. Furthermore, we have identified monocyte-derived macrophages that undergo a lipid-phenotype transition and may contribute to the suppressive TME. Intense interaction between stromal cells, myeloid cells including lipid associated macrophages and LAMP3 + DCs, and lymphocytes were also characterized.

CONCLUSIONS

Our work provides new insight into the molecular and cellular mechanism underlying malignant progression of LUAD manifested as GGN, and pave way for novel immunotherapies for GGN. Video Abstract.

Cancer-prone Phenotypes and Gene Expression Heterogeneity at Single-cell Resolution in Cigarette-smoking Lungs

Single-cell RNA sequencing (scRNA-seq) technologies have been broadly utilized to reveal molecular mechanisms of respiratory pathology and physiology at single-cell resolution. Here, we established single-cell meta-analysis (scMeta-analysis) by integrating data from eight public datasets, including 104 lung scRNA-seq samples with clinicopathologic information and designated a cigarette-smoking lung atlas. The atlas revealed early carcinogenesis events and defined the alterations of single-cell transcriptomics, cell population, and fundamental properties of biological pathways induced by smoking. In addition, we developed two novel scMeta-analysis methods: VARIED (Visualized Algorithms of Relationships In Expressional Diversity) and AGED (Aging-related Gene Expressional Differences). VARIED analysis revealed expressional diversity associated with smoking carcinogenesis. AGED analysis revealed differences in gene expression related to both aging and smoking status. The scMeta-analysis paves the way to utilize publicly-available scRNA-seq data and provide new insights into the effects of smoking and into cellular diversity in human lungs, at single-cell resolution.

SIGNIFICANCE

The atlas revealed early carcinogenesis events and defined the alterations of single-cell transcriptomics, cell population, and fundamental properties of biological pathways induced by smoking.

Intercellular Communication in Airway Epithelial Cell Regeneration: Potential Roles of Connexins and Pannexins

Connexins and pannexins are transmembrane proteins that can form direct (gap junctions) or indirect (connexons, pannexons) intercellular communication channels. By propagating ions, metabolites, sugars, nucleotides, miRNAs, and/or second messengers, they participate in a variety of physiological functions, such as tissue homeostasis and host defense. There is solid evidence supporting a role for intercellular signaling in various pulmonary inflammatory diseases where alteration of connexin/pannexin channel functional expression occurs, thus leading to abnormal intercellular communication pathways and contributing to pathophysiological aspects, such as innate immune defense and remodeling. The integrity of the airway epithelium, which is the first line of defense against invading microbes, is established and maintained by a repair mechanism that involves processes such as proliferation, migration, and differentiation. Here, we briefly summarize current knowledge on the contribution of connexins and pannexins to necessary processes of tissue repair and speculate on their possible involvement in the shaping of the airway epithelium integrity.

Integration of risk variants from GWAS with SARS-CoV-2 RNA interactome prioritizes FUBP1 and RAB2A as risk genes for COVID-19

The role of host genetic factors in COVID-19 outcomes remains unclear despite various genome-wide association studies (GWAS). We annotate all significant variants and those variants in high LD (R2 > 0.8) from the COVID-19 host genetics initiative (HGI) and identify risk genes by recognizing genes intolerant nonsynonymous mutations in coding regions and genes associated with cis-expression quantitative trait loci (cis-eQTL) in non-coding regions. These genes are enriched in the immune response pathway and viral life cycle. It has been found that host RNA binding proteins (RBPs) participate in different phases of the SARS-CoV-2 life cycle. We collect 503 RBPs that interact with SARS-CoV-2 RNA concluded from in vitro studies. Combining risk genes from the HGI with RBPs, we identify two COVID-19 risk loci that regulate the expression levels of FUBP1 and RAB2A in the lung. Due to the risk allele, COVID-19 patients show downregulation of FUBP1 and upregulation of RAB2A. Using single-cell RNA sequencing data, we show that FUBP1 and RAB2A are expressed in SARS-CoV-2-infected upper respiratory tract epithelial cells. We further identify NC_000001.11:g.77984833C>A and NC_000008.11:g.60559280T>C as functional variants by surveying allele-specific transcription factor sites and cis-regulatory elements and performing motif analysis. To sum up, our research, which associates human genetics with expression levels of RBPs, identifies FUBP1 and RAB2A as two risk genes for COVID-19 and reveals the anti-viral role of FUBP1 and the pro-viral role of RAB2A in the infection of SARS-CoV-2.

Human airway tuft cells influence the mucociliary clearance through cholinergic signalling

BACKGROUND

Airway tuft cells, formerly called brush cells have long been described only morphologically in human airways. More recent RNAseq studies described a chemosensory cell population, which includes tuft cells, by a distinct gene transcription signature. Yet, until which level in the tracheobronchial tree in native human airway epithelium tuft cells occur and if they function as regulators of innate immunity, e.g., by regulating mucociliary clearance, remained largely elusive.

METHODS

We performed immunohistochemistry, RT-PCR and immunoblotting analyses for various tuft cell markers to confirm the presence of this cell type in human tracheal samples. Immunohistochemistry was conducted to study the distribution of tuft cells along the intrapulmonary airways in humans. We assessed the influence of bitter substances and the taste transduction pathway on mucociliary clearance in mouse and human tracheal samples by measuring particle transport speed.

RESULTS

Tuft cells identified by the expression of their well-established marker POU class 2 homeobox 3 (POU2F3) were present from the trachea to the bronchioles. We identified choline acetyltransferase in POU2F3 expressing cells as well as the transient receptor potential melastatin 5 (TRPM5) channel in a small population of tracheal epithelial cells with morphological appearance of tuft cells. Application of bitter substances, such as denatonium, led to an increase in mucociliary clearance in human tracheal preparations. This was dependent on activation of the TRPM5 channel and involved cholinergic and nitric oxide signalling, indicating a functional role for human tuft cells in the regulation of mucociliary clearance.

CONCLUSIONS

We were able to detect tuft cells in the tracheobronchial tree down to the level of the bronchioles. Moreover, taste transduction and cholinergic signalling occur in the same cells and regulate mucociliary clearance. Thus, tuft cells are potentially involved in the regulation of innate immunity in human airways.

Transplantation of Induced Pluripotent Stem Cell-Derived Airway Epithelia with a Collagen Scaffold into the Nasal Cavity

The nasal cavity is covered with respiratory epithelia, including ciliated cells that eliminate foreign substances trapped in the mucus. In hereditary diseases such as primary ciliary dyskinesia and cystic fibrosis, respiratory epithelial functions are irreversibly impaired; however, no radical treatment has been established yet. Thus, we considered that the transplantation of normal airway epithelia (AE) into the nasal epithelia is one of the strategies that could lead to radical treatment in the future. In our previous study, human induced pluripotent stem cell-derived AE (hiPSC-AE) on the vitrigel membrane were transplanted into the scraped area of the nasal septal mucosa of nude rats. Although human-derived ciliated cells, club cells, and basal cells were observed, they were located in the cysts within the submucosal granulation tissue but not in the nasal mucosal epithelia and the transplanted cells may not contribute to the function of the nasal mucosa with this condition. Therefore, to achieve more functional transplantation, we prepared the graft differently in this study by wrapping the collagen sponge in hiPSC-AE on the vitrigel membrane. As a result, we found the transplanted cells surviving in the nasal mucosal epithelia. These results suggest that hiPSC-AE transplanted into the nasal cavity could be viable in the nasal mucosa. In addition, our method leads to the establishment of nasal mucosa-humanized rats that are used for the development of the drugs and therapeutic methods for hereditary diseases of nasal respiratory epithelia.

The development and characterization of in vivo-like three-dimensional models of bronchial epithelial cell lines

In vitro models of differentiated respiratory epithelium that allow high-throughput screening are an important tool to explore new therapeutics for chronic respiratory diseases. In the present study, we developed in vivo-like three-dimensional (3-D) models of bronchial epithelial cell lines that are commonly used to study chronic lung disease (16HBE14o-, CFBE41o- and CFBE41o- 6.2 WT-CFTR). To this end, cells were cultured on porous microcarrier beads in the rotating wall vessel (RWV) bioreactor, an optimized suspension culture method that allows higher throughput experimentation than other physiologically relevant models. Cell differentiation was compared to conventional two-dimensional (2-D) monolayer cultures and to the current gold standard in the respiratory field, i.e. air-liquid interface (ALI) cultures. Cellular differentiation was assessed in the three model systems by evaluating the expression and localization of markers that reflect the formation of tight junctions (zonula occludens 1), cell polarity (intercellular adhesion molecule 1 at the apical side and collagen IV expression at the basal cell side), multicellular complexity (acetylated α-tubulin for ciliated cells, CC10 for club cells, keratin-5 for basal cells) and mucus production (MUC5AC) through immunostaining and confocal laser scanning microscopy. Results were validated using Western Blot analysis. We found that tight junctions were expressed in 2-D monolayers, ALI cultures and 3-D models for all three cell lines. All tested bronchial epithelial cell lines showed polarization in ALI and 3-D cultures, but not in 2-D monolayers. Mucus secreting goblet-like cells were present in ALI and 3-D cultures of CFBE41o- and CFBE41o- 6.2 WT-CFTR cells, but not in 16HBE14o- cells. For all cell lines, there were no ciliated cells, basal cells, or club cells found in any of the model systems. In conclusion, we developed RWV-derived 3-D models of commonly used bronchial epithelial cell lines and showed that these models are a valuable alternative to ALI cultures, as they recapitulate similar key aspects of the in vivo parental tissue.

A bronchial gene signature specific for severe COPD that is retained in the nose

Introduction

A subset of COPD patients develops advanced disease with severe airflow obstruction, hyperinflation and extensive emphysema. We propose that the pathogenesis in these patients differs from mild-moderate COPD and is reflected by bronchial gene expression. The aim of the present study was to identify a unique bronchial epithelial gene signature for severe COPD patients.

Methods

We obtained RNA sequencing data from bronchial brushes from 123 ex-smokers with severe COPD, 23 with mild-moderate COPD and 23 non-COPD controls. We identified genes specific to severe COPD by comparing severe COPD to non-COPD controls, followed by removing genes that were also differentially expressed between mild-moderate COPD and non-COPD controls. Next, we performed a pathway analysis on these genes and evaluated whether this signature is retained in matched nasal brushings.

Results

We identified 219 genes uniquely differentially expressed in severe COPD. Interaction network analysis identified VEGFA and FN1 as the key genes with the most interactions. Genes were involved in extracellular matrix regulation, collagen binding and the immune response. Of interest were 10 genes (VEGFA, DCN, SPARC, COL6A2, MGP, CYR61, ANXA6, LGALS1, C1QA and C1QB) directly connected to fibronectin 1 (FN1). Most of these genes were lower expressed in severe COPD and showed the same effect in nasal brushings.

Conclusions

We found a unique severe COPD bronchial gene signature with key roles for VEGFA and FN1, which was retained in the upper airways. This supports the hypothesis that severe COPD, at least partly, comprises a different pathology and supports the potential for biomarker development based on nasal brushes in COPD.

Exploring intrinsic variability between cultured nasal and bronchial epithelia in cystic fibrosis

The nasal and bronchial epithelium are unified parts of the respiratory tract that are affected in the monogenic disorder cystic fibrosis (CF). Recent studies have uncovered that nasal and bronchial tissues exhibit intrinsic variability, including differences in mucociliary cell composition and expression of unique transcriptional regulatory proteins which relate to germ layer origin. In the present study, we explored whether intrinsic differences between nasal and bronchial epithelial cells persist in cell cultures and affect epithelial cell functioning in CF. Comparison of air-liquid interface (ALI) differentiated epithelial cells from subjects with CF revealed distinct mucociliary differentiation states of nasal and bronchial cultures. Moreover, using RNA sequencing we identified cell type-specific signature transcription factors in differentiated nasal and bronchial epithelial cells, some of which were already poised for expression in basal progenitor cells as evidenced by ATAC sequencing. Analysis of differentiated nasal and bronchial epithelial 3D organoids revealed distinct capacities for fluid secretion, which was linked to differences in ciliated cell differentiation. In conclusion, we show that unique phenotypical and functional features of nasal and bronchial epithelial cells persist in cell culture models, which can be further used to investigate the effects of tissue-specific features on upper and lower respiratory disease development in CF.

Chromatin accessibility profiling of targeted cell populations with laser capture microdissection coupled to ATAC-seq

Spatially resolved omics technologies reveal context-dependent cellular regulatory networks in tissues of interest. Beyond transcriptome analysis, information on epigenetic traits and chromatin accessibility can provide further insights on gene regulation in health and disease. Nevertheless, compared to the enormous advancements in spatial transcriptomics technologies, the field of spatial epigenomics is much younger and still underexplored. In this study, we report laser capture microdissection coupled to ATAC-seq (LCM-ATAC-seq) applied to fresh frozen samples for the spatial characterization of chromatin accessibility. We first demonstrate the efficient use of LCM coupled to in situ tagmentation and evaluate its performance as a function of cell number, microdissected areas, and tissue type. Further, we demonstrate its use for the targeted chromatin accessibility analysis of discrete contiguous or scattered cell populations in tissues via single-nuclei capture based on immunostaining for specific cellular markers.

Puzzle resolved: CFTR mediates chloride homeostasis by segregating absorption and secretion to different cell types

In the lungs, the cystic fibrosis transmembrane conductance regulator (CFTR) regulates ion transport in surface-airway epithelia and submucosal glands, thus determining airway surface liquid (ASL) volume and mucus hydration. In this issue of the JCI, Lei Lei and colleagues report that the CFTR-rich and barttin/Cl- channel-expressing ionocytes mediate chloride absorption across airway epithelia, whereas the more abundant basal cells and secretory cells mediate chloride secretion. Thus, CFTR-mediated secretion and absorption of chloride ions in the lung are segregated by cell type, which has implications for future molecular therapies for cystic fibrosis lung disease.

CFTR-rich ionocytes mediate chloride absorption across airway epithelia

The volume and composition of a thin layer of liquid covering the airway surface defend the lung from inhaled pathogens and debris. Airway epithelia secrete Cl- into the airway surface liquid through cystic fibrosis transmembrane conductance regulator (CFTR) channels, thereby increasing the volume of airway surface liquid. The discovery that pulmonary ionocytes contain high levels of CFTR led us to predict that ionocytes drive secretion. However, we found the opposite. Elevating ionocyte abundance increased liquid absorption, whereas reducing ionocyte abundance increased secretion. In contrast to other airway epithelial cells, ionocytes contained barttin/Cl- channels in their basolateral membrane. Disrupting barttin/Cl- channel function impaired liquid absorption, and overexpressing barttin/Cl- channels increased absorption. Together, apical CFTR and basolateral barttin/Cl- channels provide an electrically conductive pathway for Cl- flow through ionocytes, and the transepithelial voltage generated by apical Na+ channels drives absorption. These findings indicate that ionocytes mediate liquid absorption, and secretory cells mediate liquid secretion. Segregating these counteracting activities to distinct cell types enables epithelia to precisely control the airway surface. Moreover, the divergent role of CFTR in ionocytes and secretory cells suggests that cystic fibrosis disrupts both liquid secretion and absorption.

FORMATION OF MALIGNANT, METASTATIC SMALL CELL LUNG CANCERS THROUGH OVERPRODUCTION OF cMYC PROTEIN IN TP53 AND RB1 DEPLETED PULMONARY NEUROENDOCRINE CELLS DERIVED FROM HUMAN EMBRYONIC STEM CELLS

We recently described our initial efforts to develop a model for small cell lung cancer (SCLC) derived from human embryonic stem cells (hESCs) that were differentiated to form pulmonary neuroendocrine cells (PNECs), a putative cell of origin for neuroendocrine-positive SCLC. Although reduced expression of the tumor suppressor genes TP53 and RB1 allowed the induced PNECs to form subcutaneous growths in immune-deficient mice, the tumors did not display the aggressive characteristics of SCLC seen in human patients. Here we report that the additional, doxycycline-regulated expression of a transgene encoding wild-type or mutant cMYC protein promotes rapid growth, invasion, and metastasis of these hESC-derived cells after injection into the renal capsule. Similar to others, we find that the addition of cMYC encourages the formation of the SCLC-N subtype, marked by high levels of NEUROD1 RNA. Using paired primary and metastatic samples for RNA sequencing, we observe that the subtype of SCLC does not change upon metastatic spread and that production of NEUROD1 is maintained. We also describe histological features of these malignant, SCLC-like tumors derived from hESCs and discuss potential uses of this model in efforts to control and better understand this recalcitrant neoplasm.

Single-cell RNA sequencing of bronchoscopy specimens: development of a rapid minimal handling protocol

Single-cell RNA sequencing (scRNA-seq) is an important tool for understanding disease pathophysiology, including airway diseases. Currently, the majority of scRNA-seq studies in airway diseases have used invasive methods (airway biopsy, surgical resection), which carry inherent risks and thus present a major limitation to scRNA-seq investigation of airway pathobiology. Bronchial brushing, where the airway mucosa is sampled using a cytological brush, is a viable, less invasive method of obtaining airway cells for scRNA-seq. Here we describe the development of a rapid and minimal handling protocol for preparing single-cell suspensions from bronchial brush specimens for scRNA-seq. Our optimized protocol maximizes cell recovery and cell quality and facilitates large-scale profiling of the airway transcriptome at single-cell resolution.

Single-cell RNA sequencing analysis reveals transcriptional heterogeneity of multiple primary lung cancer

INTRODUCTION

With the advancements in early diagnosis, more and more patients with multiple primary lung cancer (MPLC) have been identified. However, the progression of MPLC involves complex changes in cell composition and metabolic function, which remains largely controversial.

OBJECTIVE

Our study aims to comprehensively reveal the cellular characteristics and inter-cellular connections of MPLC.

METHODS

We performed scRNA-seq from 23 samples of six MPLC patients, combined with bulk whole-exome sequencing. We performed trajectory analysis to investigate the transition of different cell types during the development of MPLC.

RESULTS

A total of 1 67 397 cells were sequenced derived from tumour and adjacent tissues of MPLC patients, and tumour, normal, immune and stromal cells were identified. Two states of epithelial cells were identified, which were associated with immune response and cell death, respectively. Furthermore, both CD8+ naïve and memory T cells participated in the differentiation of CD8+ T cells. The terminal states of CD8+ T cells were exhausted T cells and cytotoxic T cells, which positively regulated cell death and were implicated in the regulation of cytokine production, respectively. Two main subpopulations of B cells with distinct functions were identified, which participate in the regulation of the immune response and antigen presentation, respectively. In addition, we found a specific type of endothelial cells that were abundant in tumour samples, with an increasing trend from normal to tumour samples.

CONCLUSIONS

Our study showed the comprehensive landscape of different cells of MPLC, which might reveal the key cellular mechanisms and, therefore, may provide new insights into the early diagnosis and treatment of MPLC.

Towards Universal Cell Embeddings: Integrating Single-cell RNA-seq Datasets across Species with SATURN

Analysis of single-cell datasets generated from diverse organisms offers unprecedented opportunities to unravel fundamental evolutionary processes of conservation and diversification of cell types. However, inter-species genomic differences limit the joint analysis of cross-species datasets to homologous genes. Here, we present SATURN, a deep learning method for learning universal cell embeddings that encodes genes' biological properties using protein language models. By coupling protein embeddings from language models with RNA expression, SATURN integrates datasets profiled from different species regardless of their genomic similarity. SATURN has a unique ability to detect functionally related genes co-expressed across species, redefining differential expression for cross-species analysis. We apply SATURN to three species whole-organism atlases and frog and zebrafish embryogenesis datasets. We show that cell embeddings learnt in SATURN can be effectively used to transfer annotations across species and identify both homologous and species-specific cell types, even across evolutionarily remote species. Finally, we use SATURN to reannotate the five species Cell Atlas of Human Trabecular Meshwork and Aqueous Outflow Structures and find evidence of potentially divergent functions between glaucoma associated genes in humans and other species.

Spatial transcriptomics delineates molecular features and cellular plasticity in lung adenocarcinoma progression

Indolent (lepidic) and aggressive (micropapillary, solid, and poorly differentiated acinar) histologic subtypes often coexist within a tumor tissue of lung adenocarcinoma (LUAD), but the molecular features associated with different subtypes and their transitions remain elusive. Here, we combine spatial transcriptomics and multiplex immunohistochemistry to elucidate molecular characteristics and cellular plasticity of distinct histologic subtypes of LUAD. We delineate transcriptional reprogramming and dynamic cell signaling that determine subtype progression, especially hypoxia-induced regulatory network. Different histologic subtypes exhibit heterogeneity in dedifferentiation states. Additionally, our results show that macrophages are the most abundant cell type in LUAD, and identify different tumor-associated macrophage subpopulations that are unique to each histologic subtype, which might contribute to an immunosuppressive microenvironment. Our results provide a systematic landscape of molecular profiles that drive LUAD subtype progression, and demonstrate potentially novel therapeutic strategies and targets for invasive lung adenocarcinoma.

Epithelial plasticity and innate immune activation promote lung tissue remodeling following respiratory viral infection

Epithelial plasticity has been suggested in lungs of mice following genetic depletion of stem cells but is of unknown physiological relevance. Viral infection and chronic lung disease share similar pathological features of stem cell loss in alveoli, basal cell (BC) hyperplasia in small airways, and innate immune activation, that contribute to epithelial remodeling and loss of lung function. We show that a subset of distal airway secretory cells, intralobar serous (IS) cells, are activated to assume BC fates following influenza virus infection. Injury-induced hyperplastic BC (hBC) differ from pre-existing BC by high expression of IL-22Ra1 and undergo IL-22-dependent expansion for colonization of injured alveoli. Resolution of virus-elicited inflammation results in BC to IS re-differentiation in repopulated alveoli, and increased local expression of protective antimicrobial factors, but fails to restore normal alveolar epithelium responsible for gas exchange.

A single-cell atlas of West African lungfish respiratory system reveals evolutionary adaptations to terrestrialization

The six species of lungfish possess both lungs and gills and are the closest extant relatives of tetrapods. Here, we report a single-cell transcriptome atlas of the West African lungfish (Protopterus annectens). This species manifests the most extreme form of terrestrialization, a life history strategy to survive dry periods that can last for years, characterized by dormancy and reversible adaptive changes of the gills and lungs. Our atlas highlights the cell type diversity of the West African lungfish, including gene expression consistent with phenotype changes of terrestrialization. Comparison with terrestrial tetrapods and ray-finned fishes reveals broad homology between the swim bladder and lung cell types as well as shared and idiosyncratic changes of the external gills of the West African lungfish and the internal gills of Atlantic salmon. The single-cell atlas presented here provides a valuable resource for further exploration of the respiratory system evolution in vertebrates and the diversity of lungfish terrestrialization.

The contribution of reticular basement membrane proteins to basal airway epithelial attachment, spreading and barrier formation: implications for airway remodeling in asthma

Rationale

In the healthy lung, the pseudostratified conducting airway epithelium is anchored to the reticular basement membrane (RBM) via hemidesmosome junction complexes formed between basal cells and the extracellular matrix (ECM). The RBM within the healthy lung is composed of the ECM proteins laminin and collagen-IV. In patients with asthma, the RBM is remodeled with collagen-I, -III and fibronectin deposition. The goal of this study was to assess the effect of RBM ECM proteins on basal airway epithelial cell attachment, spreading and barrier formation using real-time electrical cell-substrate impedance sensing (ECIS).

Methods

ECIS 8-well arrays were coated with 50 μg/mL of fibronectin, collagen-I, collagen-III, collagen-IV, or laminin and compared to bovine serum albumin (BSA) or uncoated controls. The airway epithelial cell line (1HAEo-) was seeded 40, 50, 60, and 70 k cells/well and continuously monitored over 70 h to assess cell attachment, spreading and barrier formation using high (64 k Hz) and low (500 Hz) frequency resistance and capacitance. Data were analyzed using a one-phase decay model from which half-life (time cells cover half of the electrode area) and rate-constant (cell-spreading rate/h) were determined and a generalized additive mixed effect model (GAMM) was used to assess ECM proteins over the entire experiment.

Results

High-frequency (64 kHz) capacitance measures demonstrated the half-life for 1HAEo-cells to attach was fastest when grown on fibronectin (6.5 h), followed by collagen-I (7.2 h) and collagen-III (8.1 h), compared to collagen-IV (11.3 h), then laminin (13.2 h) compared to BSA (12.4 h) and uncoated (13.9 h) controls. High-frequency (64 kHz) resistance measures demonstrated that the rate of 1HAEo- cell spreading was significantly faster on fibronectin and collagen-I compared to collagen-III, collagen-IV, laminin, BSA and the uncoated control. Low-frequency (500 Hz) resistance measures demonstrated that 1HAEo-cells formed a functional barrier fastest when grown on fibronectin and collagen-I, compared to the other ECM conditions. Lastly, the distance of 1HAEo-cells from the ECM substrates was the smallest when grown on fibronectin reflecting high cell-matrix adhesion.

Conclusion

Airway epithelial cells attach, spread and form a barrier fastest on fibronectin, and collagen-I and these reticular basement membrane ECM proteins may play a protective role in preserving the epithelial barrier during airway remodeling in asthma.

The Dual Role of the Airway Epithelium in Asthma: Active Barrier and Regulator of Inflammation

Chronic airway inflammation is the cornerstone on which bronchial asthma arises, and in turn, chronic inflammation arises from a complex interplay between environmental factors such as allergens and pathogens and immune cells as well as structural cells constituting the airway mucosa. Airway epithelial cells (AECs) are at the center of these processes. On the one hand, they represent the borderline separating the body from its environment in order to keep inner homeostasis. The airway epithelium forms a multi-tiered, self-cleaning barrier that involves an unstirred, discontinuous mucous layer, the dense and rigid mesh of the glycocalyx, and the cellular layer itself, consisting of multiple, densely interconnected cell types. On the other hand, the airway epithelium represents an immunologically highly active tissue once its barrier has been penetrated: AECs play a pivotal role in releasing protective immunoglobulin A. They express a broad spectrum of pattern recognition receptors, enabling them to react to environmental stressors that overcome the mucosal barrier. By releasing alarmins-proinflammatory and regulatory cytokines-AECs play an active role in the formation, strategic orientation, and control of the subsequent defense reaction. Consequently, the airway epithelium is of vital importance to chronic inflammatory diseases, such as asthma.

Intrauterine exposure to nicotine through maternal vaping disrupts embryonic lung and skeletal development via the Kcnj2 potassium channel

Smoking cigarettes during pregnancy is associated with adverse effects on infants including low birth weight, defective lung development, and skeletal abnormalities. Pregnant women are increasingly turning to vaping [use of electronic (e)-cigarettes] as a perceived safer alternative to cigarettes. However, nicotine disrupts fetal development, suggesting that like cigarette smoking, nicotine vaping may be detrimental to the fetus. To test the impact of maternal vaping on fetal lung and skeletal development in mice, pregnant dams were exposed to e-cigarette vapor throughout gestation. At embryonic day (E)18.5, vape exposed litter sizes were reduced, and some embryos exhibited growth restriction compared to air exposed controls. Fetal lungs were collected for histology and whole transcriptome sequencing. Maternally nicotine vaped embryos exhibited histological and transcriptional changes consistent with impaired distal lung development. Embryonic lung gene expression changes mimicked transcriptional changes observed in adult mouse lungs exposed to cigarette smoke, suggesting that the developmental defects may be due to direct nicotine exposure. Fetal skeletons were analyzed for craniofacial and long bone lengths. Nicotine directly binds and inhibits the Kcnj2 potassium channel which is important for bone development. The length of the maxilla, palatal shelves, humerus, and femur were reduced in vaped embryos, which was further exacerbated by loss of one copy of the Kcnj2 gene. Nicotine vapor exposed Kcnj2KO/+ embryos also had significantly lower birth weights than unexposed animals of either genotype. Kcnj2 mutants had severely defective lungs with and without vape exposure, suggesting that potassium channels may be broadly involved in mediating the detrimental developmental effects of nicotine vaping. These data indicate that intrauterine nicotine exposure disrupts fetal lung and skeletal development likely through inhibition of Kcnj2.

The viral fitness and intrinsic pathogenicity of dominant SARS-CoV-2 Omicron sublineages BA.1, BA.2, and BA.5

BACKGROUND

Among the Omicron sublineages that have emerged, BA.1, BA.2, BA.5, and their related sublineages have resulted in the largest number of infections. While recent studies demonstrated that all Omicron sublineages robustly escape neutralizing antibody response, it remains unclear on whether these Omicron sublineages share any pattern of evolutionary trajectory on their replication efficiency and intrinsic pathogenicity along the respiratory tract.

METHODS

We compared the virological features, replication capacity of dominant Omicron sublineages BA.1, BA.2 and BA.5 in the human nasal epithelium, and characterized their pathogenicity in K18-hACE2, A129, young C57BL/6, and aged C57BL/6 mice.

FINDINGS

We found that BA.5 replicated most robustly, followed by BA.2 and BA.1, in the differentiated human nasal epithelium. Consistently, BA.5 infection resulted in higher viral gene copies, infectious viral titres and more abundant viral antigen expression in the nasal turbinates of the infected K18-hACE2 transgenic mice. In contrast, the Omicron sublineages are continuously attenuated in lungs of infected K18-hACE2 and C57BL/6 mice, leading to decreased pathogenicity. Nevertheless, lung manifestations remain severe in Omicron sublineages-infected A129 and aged C57BL/6 mice.

INTERPRETATION

Our results suggested that the Omicron sublineages might be gaining intrinsic replication fitness in the upper respiratory tract, therefore highlighting the importance of global surveillance of the emergence of hyper-transmissive Omicron sublineages. On the contrary, replication and intrinsic pathogenicity of Omicron is suggested to be further attenuated in the lower respiratory tract. Effective vaccination and other precautions should be in place to prevent severe infections in the immunocompromised populations at risk.

FUNDING

A full list of funding bodies that contributed to this study can be found in the Acknowledgements section.

Ionocyte-Specific Regulation of Cystic Fibrosis Transmembrane Conductance Regulator

CFTR (cystic fibrosis transmembrane conductance regulator) is a tightly regulated anion channel that mediates chloride and bicarbonate conductance in many epithelia and in other tissues, but whether its regulation varies depending on the cell type has not been investigated. Epithelial CFTR expression is highest in rare cells called ionocytes. We studied CFTR regulation in control and ionocyte-enriched cultures by transducing bronchial basal cells with adenoviruses that encode only eGFP or FOXI1 (forkhead box I1) + eGFP as separate polypeptides. FOXI1 dramatically increased the number of transcripts for ionocyte markers ASCL3 (Achaete-Scute Family BHLH Transcription Factor 3), BSND, ATP6V1G3, ATP6V0D2, KCNMA1, and CFTR without altering those for secretory (SCGB1A1), basal (KRT5, KRT6, TP63), goblet (MUC5AC), or ciliated (FOXJ1) cells. The number of cells displaying strong FOXI1 expression was increased 7-fold, and there was no evidence for a broad increase in background immunofluorescence. Total CFTR mRNA and protein levels increased 10-fold and 2.5-fold, respectively. Ionocyte-enriched cultures displayed elevated basal current, increased adenylyl cyclase 5 expression, and tonic suppression of CFTR activity by the phosphodiesterase PDE1C, which has not been shown previously to regulate CFTR activity. The results indicate that CFTR regulation depends on cell type and identifies PDE1C as a potential target for therapeutics that aim to increase CFTR function specifically in ionocytes.

Spatial transcriptomics: recent developments and insights in respiratory research

The respiratory system's complex cellular heterogeneity presents unique challenges to researchers in this field. Although bulk RNA sequencing and single-cell RNA sequencing (scRNA-seq) have provided insights into cell types and heterogeneity in the respiratory system, the relevant specific spatial localization and cellular interactions have not been clearly elucidated. Spatial transcriptomics (ST) has filled this gap and has been widely used in respiratory studies. This review focuses on the latest iterative technology of ST in recent years, summarizing how ST can be applied to the physiological and pathological processes of the respiratory system, with emphasis on the lungs. Finally, the current challenges and potential development directions are proposed, including high-throughput full-length transcriptome, integration of multi-omics, temporal and spatial omics, bioinformatics analysis, etc. These viewpoints are expected to advance the study of systematic mechanisms, including respiratory studies.

Human Lung Organoids-A Novel Experimental and Precision Medicine Approach

The global burden of respiratory diseases is very high and still on the rise, prompting the need for accurate models for basic and translational research. Several model systems are currently available ranging from simple airway cell cultures to complex tissue-engineered lungs. In recent years, human lung organoids have been established as highly transferrable three-dimensional in vitro model systems for lung research. For acute infectious and chronic inflammatory diseases as well as lung cancer, human lung organoids have opened possibilities for precise in vitro research and a deeper understanding of mechanisms underlying lung injury and regeneration. Human lung organoids from induced pluripotent stem cells or from adult stem cells of patients' samples introduce tools for understanding developmental processes and personalized medicine approaches. When further state-of-the-art technologies and protocols come into use, the full potential of human lung organoids can be harnessed. High-throughput assays in drug development, gene therapy, and organoid transplantation are current applications of organoids in translational research. In this review, we emphasize novel approaches in translational and personalized medicine in lung research focusing on the use of human lung organoids.

A WNT4- and DKK3-driven canonical to noncanonical Wnt signaling switch controls multiciliogenesis

Multiciliated cells contain hundreds of cilia whose directional movement powers the mucociliary clearance of the airways, a vital host defense mechanism. Multiciliated cell specification requires canonical Wnt signaling, which then must be turned off. Next, ciliogenesis and polarized ciliary orientation are regulated by noncanonical Wnt/planar cell polarity (Wnt/PCP) signaling. The mechanistic relationship between the Wnt pathways is unknown. We show that DKK3, a secreted canonical Wnt regulator and WNT4, a noncanonical Wnt ligand act together to facilitate a canonical to noncanonical Wnt signaling switch during multiciliated cell formation. In primary human airway epithelial cells, DKK3 and WNT4 CRISPR knockout blocks, whereas ectopic expression promotes, multiciliated cell formation by inhibiting canonical Wnt signaling. Wnt4 and Dkk3 single-knockout mice also display defective ciliated cells. DKK3 and WNT4 are co-secreted from basal stem cells and act directly on multiciliated cells via KREMEN1 and FZD6, respectively. We provide a novel mechanism that links specification to cilium biogenesis and polarization for proper multiciliated cell formation.

High-throughput bioprinting of the nasal epithelium using patient-derived nasal epithelial cells

Progenitor human nasal epithelial cells (hNECs) are an essential cell source for the reconstruction of the respiratory pseudostratified columnar epithelium composed of multiple cell types in the context of infection studies and disease modeling. Hitherto, manual seeding has been the dominant method for creating nasal epithelial tissue models through biofabrication. However, this approach has limitations in terms of achieving the intricate three-dimensional (3D) structure of the natural nasal epithelium. 3D bioprinting has been utilized to reconstruct various epithelial tissue models, such as cutaneous, intestinal, alveolar, and bronchial epithelium, but there has been no attempt to use of 3D bioprinting technologies for reconstruction of the nasal epithelium. In this study, for the first time, we demonstrate the reconstruction of the nasal epithelium with the use of primary hNECs deposited on Transwell inserts via droplet-based bioprinting (DBB), which enabled high-throughput fabrication of the nasal epithelium in Transwell inserts of 24-well plates. DBB of progenitor hNECs ranging from one-tenth to one-half of the cell seeding density employed during the conventional cell seeding approach enabled a high degree of differentiation with the presence of cilia and tight-junctions over a 4 weeks air-liquid interface culture. Single cell RNA sequencing of these cultures identified five major epithelial cells populations, including basal, suprabasal, goblet, club, and ciliated cells. These cultures recapitulated the pseudostratified columnar epithelial architecture present in the native nasal epithelium and were permissive to respiratory virus infection. These results denote the potential of 3D bioprinting for high-throughput fabrication of nasal epithelial tissue models not only for infection studies but also for other purposes, such as disease modeling, immunological studies, and drug screening.

Glycosylated extracellular mucin domains protect against SARS-CoV-2 infection at the respiratory surface

Mucins play an essential role in protecting the respiratory tract against microbial infections while also acting as binding sites for bacterial and viral adhesins. The heavily O-glycosylated gel-forming mucins MUC5AC and MUC5B eliminate pathogens by mucociliary clearance. Transmembrane mucins MUC1, MUC4, and MUC16 can restrict microbial invasion at the apical surface of the epithelium. In this study, we determined the impact of host mucins and mucin glycans on epithelial entry of SARS-CoV-2. Human lung epithelial Calu-3 cells express the SARS-CoV-2 entry receptor ACE2 and high levels of glycosylated MUC1, but not MUC4 and MUC16, on their cell surface. The O-glycan-specific mucinase StcE specifically removed the glycosylated part of the MUC1 extracellular domain while leaving the underlying SEA domain and cytoplasmic tail intact. StcE treatment of Calu-3 cells significantly enhanced infection with SARS-CoV-2 pseudovirus and authentic virus, while removal of terminal mucin glycans sialic acid and fucose from the epithelial surface did not impact viral entry. In Calu-3 cells, the transmembrane mucin MUC1 and ACE2 are located to the apical surface in close proximity and StcE treatment results in enhanced binding of purified spike protein. Both MUC1 and MUC16 are expressed on the surface of human organoid-derived air-liquid interface (ALI) differentiated airway cultures and StcE treatment led to mucin removal and increased levels of SARS-CoV-2 replication. In these cultures, MUC1 was highly expressed in non-ciliated cells while MUC16 was enriched in goblet cells. In conclusion, the glycosylated extracellular domains of different transmembrane mucins might have similar protective functions in different respiratory cell types by restricting SARS-CoV-2 binding and entry.

Pseudomonas aeruginosa type IV pili actively induce mucus contraction to form biofilms in tissue-engineered human airways

The opportunistic pathogen Pseudomonas aeruginosa causes antibiotic-recalcitrant pneumonia by forming biofilms in the respiratory tract. Despite extensive in vitro experimentation, how P. aeruginosa forms biofilms at the airway mucosa is unresolved. To investigate the process of biofilm formation in realistic conditions, we developed AirGels: 3D, optically accessible tissue-engineered human lung models that emulate the airway mucosal environment. AirGels recapitulate important factors that mediate host-pathogen interactions including mucus secretion, flow and air-liquid interface (ALI), while accommodating high-resolution live microscopy. With AirGels, we investigated the contributions of mucus to P. aeruginosa biofilm biogenesis in in vivo-like conditions. We found that P. aeruginosa forms mucus-associated biofilms within hours by contracting luminal mucus early during colonization. Mucus contractions facilitate aggregation, thereby nucleating biofilms. We show that P. aeruginosa actively contracts mucus using retractile filaments called type IV pili. Our results therefore suggest that, while protecting epithelia, mucus constitutes a breeding ground for biofilms.

A single-cell atlas of in vitro multiculture systems uncovers the in vivo lineage trajectory and cell state in the human lung

We present an in-depth single-cell atlas of in vitro multiculture systems on human primary airway epithelium derived from normal and diseased lungs of 27 individual donors. Our large-scale single-cell profiling identified new cell states and differentiation trajectories of rare airway epithelial cell types in human distal lungs. By integrating single-cell datasets of human lung tissues, we discovered immune-primed subsets enriched in lungs and organoids derived from patients with chronic respiratory disease. To demonstrate the full potential of our platform, we further illustrate transcriptomic responses to various respiratory virus infections in vitro airway models. Our work constitutes a single-cell roadmap for the cellular and molecular characteristics of human primary lung cells in vitro and their relevance to human tissues in vivo.

Rhinovirus C causes heterogeneous infection and gene expression in airway epithelial cell subsets

Rhinoviruses infect ciliated airway epithelial cells, and rhinoviruses' nonstructural proteins quickly inhibit and divert cellular processes for viral replication. However, the epithelium can mount a robust innate antiviral immune response. Therefore, we hypothesized that uninfected cells contribute significantly to the antiviral immune response in the airway epithelium. Using single-cell RNA sequencing, we demonstrate that both infected and uninfected cells upregulate antiviral genes (e.g. MX1, IFIT2, IFIH1, and OAS3) with nearly identical kinetics, whereas uninfected non-ciliated cells are the primary source of proinflammatory chemokines. Furthermore, we identified a subset of highly infectable ciliated epithelial cells with minimal interferon responses and determined that interferon responses originate from distinct subsets of ciliated cells with moderate viral replication. These findings suggest that the composition of ciliated airway epithelial cells and coordinated responses of infected and uninfected cells could determine the risk of more severe viral respiratory illnesses in children with asthma, chronic obstructive pulmonary disease, and genetically susceptible individuals.

Airway Epithelial-Derived Immune Mediators in COVID-19

The airway epithelium, which lines the conducting airways, is central to the defense of the lungs against inhaled particulate matter and pathogens such as SARS-CoV-2, the virus that causes COVID-19. Recognition of pathogens results in the activation of an innate and intermediate immune response which involves the release of cytokines and chemokines by the airway epithelium. This response can inhibit further viral invasion and influence adaptive immunity. However, severe COVID-19 is characterized by a hyper-inflammatory response which can give rise to clinical presentations including lung injury and lead to acute respiratory distress syndrome, viral pneumonia, coagulopathy, and multi-system organ failure. In response to SARS-CoV-2 infection, the airway epithelium can mount a maladaptive immune response which can delay viral clearance, perpetuate excessive inflammation, and contribute to the pathogenesis of severe COVID-19. In this article, we will review the barrier and immune functions of the airway epithelium, how SARS-CoV-2 can interact with the epithelium, and epithelial-derived cytokines and chemokines and their roles in COVID-19 and as biomarkers. Finally, we will discuss these immune mediators and their potential as therapeutic targets in COVID-19.

Guided construction of single cell reference for human and mouse lung

Accurate cell type identification is a key and rate-limiting step in single-cell data analysis. Single-cell references with comprehensive cell types, reproducible and functionally validated cell identities, and common nomenclatures are much needed by the research community for automated cell type annotation, data integration, and data sharing. Here, we develop a computational pipeline utilizing the LungMAP CellCards as a dictionary to consolidate single-cell transcriptomic datasets of 104 human lungs and 17 mouse lung samples to construct LungMAP single-cell reference (CellRef) for both normal human and mouse lungs. CellRefs define 48 human and 40 mouse lung cell types catalogued from diverse anatomic locations and developmental time points. We demonstrate the accuracy and stability of LungMAP CellRefs and their utility for automated cell type annotation of both normal and diseased lungs using multiple independent methods and testing data. We develop user-friendly web interfaces for easy access and maximal utilization of the LungMAP CellRefs.

Respiratory epithelial cell types, states and fates in the era of single-cell RNA-sequencing

Standalone and consortia-led single-cell atlases of healthy and diseased human airways generated with single-cell RNA-sequencing (scRNA-seq) have ushered in a new era in respiratory research. Numerous discoveries, including the pulmonary ionocyte, potentially novel cell fates, and a diversity of cell states among common and rare epithelial cell types have highlighted the extent of cellular heterogeneity and plasticity in the respiratory tract. scRNA-seq has also played a pivotal role in our understanding of host-virus interactions in coronavirus disease 2019 (COVID-19). However, as our ability to generate large quantities of scRNA-seq data increases, along with a growing number of scRNA-seq protocols and data analysis methods, new challenges related to the contextualisation and downstream applications of insights are arising. Here, we review the fundamental concept of cellular identity from the perspective of single-cell transcriptomics in the respiratory context, drawing attention to the need to generate reference annotations and to standardise the terminology used in literature. Findings about airway epithelial cell types, states and fates obtained from scRNA-seq experiments are compared and contrasted with information accumulated through the use of conventional methods. This review attempts to discuss major opportunities and to outline some of the key limitations of the modern-day scRNA-seq that need to be addressed to enable efficient and meaningful integration of scRNA-seq data from different platforms and studies, with each other as well as with data from other high-throughput sequencing-based genomic, transcriptomic and epigenetic analyses.

Airway transcriptome networks identify susceptibility to frequent asthma exacerbations in children

BACKGROUND

Frequent asthma exacerbators, defined as those experiencing more than 1 hospitalization in a year for an asthma exacerbation, represent an important subgroup of individuals with asthma. However, this group remains poorly defined and understudied in children.

OBJECTIVE

Our aim was to determine the molecular mechanisms underlying asthma pathogenesis and exacerbation frequency.

METHODS

We performed RNA sequencing of upper airway cells from both frequent and nonfrequent exacerbators enrolled in the Ohio Pediatric Asthma Repository.

RESULTS

Through molecular network analysis, we found that nonfrequent exacerbators display an increase in modules enriched for immune system processes, including type 2 inflammation and response to infection. In contrast, frequent exacerbators showed expression of modules enriched for nervous system processes, such as synaptic formation and axonal outgrowth.

CONCLUSION

These data suggest that the upper airway of frequent exacerbators undergoes peripheral nervous system remodeling, representing a novel mechanism underlying pediatric asthma exacerbation.

Lipid-associated macrophages transition to an inflammatory state in human atherosclerosis increasing the risk of cerebrovascular complications

The immune system is integral to cardiovascular health and disease. Targeting inflammation ameliorates adverse cardiovascular outcomes. Atherosclerosis, a major underlying cause of cardiovascular disease (CVD), is conceptualised as a lipid-driven inflammation where macrophages play a non-redundant role. However, evidence emerging so far from single cell atlases suggests a dichotomy between lipid associated and inflammatory macrophage states. Here, we present an inclusive reference atlas of human intraplaque immune cell communities. Combining scRNASeq of human surgical carotid endarterectomies in a discovery cohort with bulk RNASeq and immunohistochemistry in a validation cohort (the Carotid Plaque Imaging Project-CPIP), we reveal the existence of PLIN2hi/TREM1hi macrophages as a toll-like receptor-dependent inflammatory lipid-associated macrophage state linked to cerebrovascular events. Our study shifts the current paradigm of lipid-driven inflammation by providing biological evidence for a pathogenic macrophage transition to an inflammatory lipid-associated phenotype and for its targeting as a new treatment strategy for CVD.

Tetraspanin-8 sequesters syntaxin-2 to control biphasic release propensity of mucin granules

Agonist-mediated stimulated pathway of mucin and insulin release are biphasic in which rapid fusion of pre-docked granules is followed by slow docking and fusion of granules from the reserve pool. Here, based on a cell-culture system, we show that plasma membrane-located tetraspanin-8 sequesters syntaxin-2 to control mucin release. Tetraspanin-8 affects fusion of granules during the second phase of stimulated mucin release. The tetraspanin-8/syntaxin-2 complex does not contain VAMP-8, which functions with syntaxin-2 to mediate granule fusion. We suggest that by sequestering syntaxin-2, tetraspanin-8 prevents docking of granules from the reserve pool. In the absence of tetraspanin-8, more syntaxin-2 is available for docking and fusion of granules and thus doubles the quantities of mucins secreted. This principle also applies to insulin release and we suggest a cell type specific Tetraspanin/Syntaxin combination is a general mechanism regulating the fusion of dense core granules.

Dicoumarol is an effective post-exposure prophylactic for SARS-CoV-2 Omicron infection in human airway epithelium

Repurposing existing drugs to inhibit SARS-CoV-2 infection in airway epithelial cells (AECs) is a quick way to find novel treatments for COVID-19. Computational screening has found dicoumarol (DCM), a natural anticoagulant, to be a potential SARS-CoV-2 inhibitor, but its inhibitory effects and possible working mechanisms remain unknown. Using air-liquid interface culture of primary human AECs, we demonstrated that DCM has potent antiviral activity against the infection of multiple Omicron variants (including BA.1, BQ.1 and XBB.1). Time-of-addition and drug withdrawal assays revealed that early treatment (continuously incubated after viral absorption) of DCM could markedly inhibit Omicron replication in AECs, but DCM did not affect the absorption, exocytosis and spread of viruses or directly eliminate viruses. Mechanistically, we performed single-cell sequencing analysis (a database of 77,969 cells from different airway locations from 10 healthy volunteers) and immunofluorescence staining, and showed that the expression of NAD(P)H quinone oxidoreductase 1 (NQO1), one of the known DCM targets, was predominantly localised in ciliated AECs. We further found that the NQO1 expression level was positively correlated with both the disease severity of COVID-19 patients and virus copy levels in cultured AECs. In addition, DCM treatment downregulated NQO1 expression and disrupted signalling pathways associated with SARS-CoV-2 disease outcomes (e.g., Endocytosis and COVID-19 signalling pathways) in cultured AECs. Collectively, we demonstrated that DCM is an effective post-exposure prophylactic for SARS-CoV-2 infection in the human AECs, and these findings could help physicians formulate novel treatment strategies for COVID-19.

Hallmarks of transcriptional intratumour heterogeneity across a thousand tumours

Each tumour contains diverse cellular states that underlie intratumour heterogeneity (ITH), a central challenge of cancer therapeutics1. Dozens of recent studies have begun to describe ITH by single-cell RNA sequencing, but each study typically profiled only a small number of tumours and provided a narrow view of transcriptional ITH2. Here we curate, annotate and integrate the data from 77 different studies to reveal the patterns of transcriptional ITH across 1,163 tumour samples covering 24 tumour types. Among the malignant cells, we identify 41 consensus meta-programs, each consisting of dozens of genes that are coordinately upregulated in subpopulations of cells within many tumours. The meta-programs cover diverse cellular processes including both generic (for example, cell cycle and stress) and lineage-specific patterns that we map into 11 hallmarks of transcriptional ITH. Most meta-programs of carcinoma cells are similar to those identified in non-malignant epithelial cells, suggesting that a large fraction of malignant ITH programs are variable even before oncogenesis, reflecting the biology of their cell of origin. We further extended the meta-program analysis to six common non-malignant cell types and utilize these to map cell-cell interactions within the tumour microenvironment. In summary, we have assembled a comprehensive pan-cancer single-cell RNA-sequencing dataset, which is available through the Curated Cancer Cell Atlas website, and leveraged this dataset to carry out a systematic characterization of transcriptional ITH.