Lung organogenesis requires precise timing and coordination to effect spatial organization and function of the parenchymal cells. To provide a systematic broad-based view of the mechanisms governing the dynamic alterations in parenchymal cells over crucial periods of development, we performed a single-cell RNA-sequencing time-series yielding 102,571 epithelial, endothelial and mesenchymal cells across nine time points from embryonic day 12 to postnatal day 14 in mice. Combining computational fate-likelihood prediction with RNA in situ hybridization and immunofluorescence, we explore lineage relationships during the saccular to alveolar stage transition. The utility of this publicly searchable atlas resource (www.sucrelab.org/lungcells) is exemplified by discoveries of the complexity of type 1 pneumocyte function and characterization of mesenchymal Wnt expression patterns during the saccular and alveolar stages – wherein major expansion of the gas-exchange surface occurs. We provide an integrated view of cellular dynamics in epithelial, endothelial and mesenchymal cell populations during lung organogenesis.

Large-scale changes in the structure and cellular makeup of the distal lung are a hallmark of pulmonary fibrosis (PF), but the spatial contexts that contribute to disease pathogenesis have remained uncertain. Using image-based spatial transcriptomics, we analyzed the gene expression of 1.6 million cells from 35 unique lungs. Through complementary cell-based and innovative cell-agnostic analyses, we characterized the localization of PF-emergent cell types, established the cellular and molecular basis of classical PF histopathologic features and identified a diversity of distinct molecularly defined spatial niches in control and PF lungs. Using machine learning and trajectory analysis to segment and rank airspaces on a gradient of remodeling severity, we identified compositional and molecular changes associated with progressive distal lung pathology, beginning with alveolar epithelial dysregulation and culminating with changes in macrophage polarization. Together, these results provide a unique, spatially resolved view of PF and establish methods that could be applied to other spatial transcriptomic studies.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) novel coronavirus 2019 (COVID-19) global pandemic has led to millions of cases and hundreds of thousands of deaths. While older adults appear at high risk for severe disease, hospitalizations and deaths due to SARS-CoV-2 among children have been relatively rare. Integrating single-cell RNA sequencing (scRNA-seq) of developing mouse lung with temporally resolved immunofluorescence in mouse and human lung tissue, we found that expression of SARS-CoV-2 Spike protein primer TMPRSS2 was highest in ciliated cells and type I alveolar epithelial cells (AT1), and TMPRSS2 expression increased with aging in mice and humans. Analysis of autopsy tissue from fatal COVID-19 cases detected SARS-CoV-2 RNA most frequently in ciliated and secretory cells in airway epithelium and AT1 cells in peripheral lung. SARS-CoV-2 RNA was highly colocalized in cells expressing TMPRSS2. Together, these data demonstrate the cellular spectrum infected by SARS-CoV-2 in lung epithelium and suggest that developmental regulation of TMPRSS2 may underlie the relative protection of infants and children from severe respiratory illness.

The continued emergence and spread of infectious agents is of great concern, and systems biology approaches to infectious disease research can advance our understanding of host–pathogen relationships and facilitate the development of new therapies and vaccines. Molecular characterization of infectious samples outside of appropriate biosafety containment can take place only subsequent to pathogen inactivation. Herein, we describe a modified Folch extraction using chloroform/methanol that facilitates the molecular characterization of infectious samples by enabling simultaneous pathogen inactivation and extraction of proteins, metabolites, and lipids for subsequent mass spectrometry-based multi-omics measurements. This single-sample metabolite, protein and lipid extraction (MPLEx) method resulted in complete inactivation of clinically important bacterial and viral pathogens with exposed lipid membranes, including Yersinia pestis, Salmonella Typhimurium, and Campylobacter jejuni in pure culture, and Yersinia pestis, Campylobacter jejuni, and West Nile, MERS-CoV, Ebola, and influenza H7N9 viruses in infection studies. In addition, >99% inactivation, which increased with solvent exposure time, was also observed for pathogens without exposed lipid membranes including community-associated methicillin-resistant Staphylococcus aureus, Clostridium difficile spores and vegetative cells, and adenovirus type 5. The overall pipeline of inactivation and subsequent proteomic, metabolomic, and lipidomic analyses was evaluated using a human epithelial lung cell line infected with wild-type and mutant influenza H7N9 viruses, thereby demonstrating that MPLEx yields biomaterial of sufficient quality for subsequent multi-omics analyses. Based on these experimental results, we believe that MPLEx will facilitate systems biology studies of infectious samples by enabling simultaneous pathogen inactivation and multi-omics measurements from a single specimen with high success for pathogens with exposed lipid membranes.

Campylobacter jejuni, a zoonotic pathogen that frequently colonizes poultry, possesses two Microbial Surface Components Recognizing Adhesive Matrix Molecule(s) (MSCRAMMs) termed CadF and FlpA that bind to the glycoprotein fibronectin (FN). Previous to this study, it was not known whether the CadF and FlpA proteins were functionally redundant or if both were required to potentiate host cell binding and signaling processes. We addressed these questions by generating a complete repertoire of cadF and flpA mutants and complemented isolates, and performing multiple phenotypic assays. Both CadF and FlpA were found to be necessary for the maximal binding of C. jejuni to FN and to host cells. In addition, both CadF and FlpA are required for the delivery of the C. jejuni Cia effector proteins into the cytosol of host target cells, which in turn activates the MAPK signaling pathway (Erk 1/2) that is required for the C. jejuni invasion of host cells. These data demonstrate the non-redundant and bi-functional nature of these two C. jejuni FN-binding proteins. Taken together, the C. jejuni CadF and FlpA adhesins facilitate the binding of C. jejuni to the host cells, permit delivery of effector proteins into the cytosol of a host target cell, and aid in the rewiring of host cell signaling pathways to alter host cell behavior.

Lung injury in preterm infants leads to structural and functional respiratory deficits, with a risk for bronchopulmonary dysplasia (BPD) that in its most severe form is accompanied by pulmonary hypertension (PH). To identify potential cellular and molecular drivers of BPD in humans, we performed single-cell RNA sequencing of preterm infant lungs with evolving BPD and BPD + PH compared to term infants. Examination of endothelial cells reveals a unique, aberrant capillary cell-state in BPD + PH defined by ANKRD1 expression. Within the alveolar parenchyma in infants with BPD/BPD + PH, predictive signaling analysis identifies surprising deficits in the semaphorin guidance-cue pathway, with decreased expression of pro-angiogenic transcription factor FOXF1. Loss of semaphorin signaling is replicated in a murine BPD model and in humans with causal FOXF1 mutations for alveolar capillary dysplasia (ACDMPV), suggesting a mechanistic link between developmental programs underlying BPD and ACDMPV and uncovering a critical role for semaphorin signaling in normal lung development.