Single-Cell Transcriptomic Profiling Identifies Molecular Phenotypes of Newborn Human Lung Cells.

While animal model studies have extensively defined the mechanisms controlling cell diversity in the developing mammalian lung, there exists a significant knowledge gap with regards to late-stage human lung development. The NHLBI Molecular Atlas of Lung Development Program (LungMAP) seeks to fill this gap by creating a structural, cellular and molecular atlas of the human and mouse lung. Transcriptomic profiling at the single-cell level created a cellular atlas of newborn human lungs. Frozen single-cell isolates obtained from two newborn human lungs from the LungMAP Human Tissue Core Biorepository, were captured, and library preparation was completed on the Chromium 10X system. Data was analyzed in Seurat, and cellular annotation was performed using the ToppGene functional analysis tool. Transcriptional interrogation of 5500 newborn human lung cells identified distinct clusters representing multiple populations of epithelial, endothelial, fibroblasts, pericytes, smooth muscle, immune cells and their gene signatures. Computational integration of data from newborn human cells and with 32,000 cells from postnatal days 1 through 10 mouse lungs generated by the LungMAP Cincinnati Research Center facilitated the identification of distinct cellular lineages among all the major cell types. Integration of the newborn human and mouse cellular transcriptomes also demonstrated cell type-specific differences in maturation states of newborn human lung cells. Specifically, newborn human lung matrix fibroblasts could be separated into those representative of younger cells (n = 393), or older cells (n = 158). Cells with each molecular profile were spatially resolved within newborn human lung tissue. This is the first comprehensive molecular map of the cellular landscape of neonatal human lung, including biomarkers for cells at distinct states of maturity.

Deciphering Endothelial and Mesenchymal Organ Specification in Vascularized Lung and Intestinal Organoids.

To investigate the co-development of vasculature, mesenchyme, and epithelium crucial for organogenesis and the acquisition of organ-specific characteristics, we constructed a human pluripotent stem cell-derived organoid system comprising lung or intestinal epithelium surrounded by organotypic mesenchyme and vasculature. We demonstrated the pivotal role of co-differentiating mesoderm and endoderm via precise BMP regulation in generating multilineage organoids and gut tube patterning. Single-cell RNA-seq analysis revealed organ specificity in endothelium and mesenchyme, and uncovered key ligands driving endothelial specification in the lung (e.g., WNT2B and Semaphorins) or intestine (e.g., GDF15). Upon transplantation under the kidney capsule in mice, these organoids further matured and developed perfusable human-specific sub-epithelial capillaries. Additionally, our model recapitulated the abnormal endothelial-epithelial crosstalk in patients with FOXF1 deletion or mutations. Multilineage organoids provide a unique platform to study developmental cues guiding endothelial and mesenchymal cell fate determination, and investigate intricate cell-cell communications in human organogenesis and disease. Highlights: BMP signaling fine-tunes the co-differentiation of mesoderm and endoderm.The cellular composition in multilineage organoids resembles that of human fetal organs.Mesenchyme and endothelium co-developed within the organoids adopt organ-specific characteristics.Multilineage organoids recapitulate abnormal endothelial-epithelial crosstalk in FOXF1-associated disorders.

A novel non-recurrent CNV deletion involving TBX4 and leaving TBX2 intact causes congenital alveolar dysplasia.

Congenital alveolar dysplasia (CAD) belongs to rare lethal lung developmental disorders (LLDDs) in neonates, manifesting with acute respiratory failure and pulmonary arterial hypertension refractory to treatment. The majority of CAD cases have been associated with copy-number variant (CNV) deletions at 17q23.1q23.2 or 5p12. Most CNV deletions at 17q23.1q23.2 were recurrent and encompassed two closely located genes, TBX4 and TBX2. In a few CAD cases, intragenic frameshifting deletions or single-nucleotide variants (SNVs) involved TBX4 but not TBX2. Here, we describe a male neonate who died at 27 days of life from acute respiratory failure caused by lung growth arrest along the spectrum of CAD confirmed by histopathological assessment. Trio-based genome sequencing revealed in the proband a novel non-recurrent ~1.07 Mb heterozygous CNV deletion at 17q23.2, encompassing TBX4 that arose de novo on the paternal chromosome. This is the first report of a larger-sized CNV deletion in a CAD patient involving TBX4 and leaving TBX2 intact. Our results, together with previous reports, indicate that perturbations of TBX4, rather than TBX2, cause severe lung phenotypes in humans.

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.

Single Cell Multiomics Identifies Cells and Genetic Networks Underlying Alveolar Capillary Dysplasia.

Rationale: Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a lethal developmental disorder of lung morphogenesis caused by insufficiency of FOXF1 (forkhead box F1) transcription factor function. The cellular and transcriptional mechanisms by which FOXF1 deficiency disrupts human lung formation are unknown. Objectives: To identify cell types, gene networks, and cell-cell interactions underlying the pathogenesis of ACDMPV. Methods: We used single-nucleus RNA and assay for transposase-accessible chromatin sequencing, immunofluorescence confocal microscopy, and RNA in situ hybridization to identify cell types and molecular networks influenced by FOXF1 in ACDMPV lungs. Measurements and Main Results: Pathogenic single-nucleotide variants and copy-number variant deletions involving the FOXF1 gene locus in all subjects with ACDMPV (n = 6) were accompanied by marked changes in lung structure, including deficient alveolar development and a paucity of pulmonary microvasculature. Single-nucleus RNA and assay for transposase-accessible chromatin sequencing identified alterations in cell number and gene expression in endothelial cells (ECs), pericytes, fibroblasts, and epithelial cells in ACDMPV lungs. Distinct cell-autonomous roles for FOXF1 in capillary ECs and pericytes were identified. Pathogenic variants involving the FOXF1 gene locus disrupt gene expression in EC progenitors, inhibiting the differentiation or survival of capillary 2 ECs and cell-cell interactions necessary for both pulmonary vasculogenesis and alveolar type 1 cell differentiation. Loss of the pulmonary microvasculature was associated with increased VEGFA (vascular endothelial growth factor A) signaling and marked expansion of systemic bronchial ECs expressing COL15A1 (collagen type XV α 1 chain). Conclusions: Distinct FOXF1 gene regulatory networks were identified in subsets of pulmonary endothelial and fibroblast progenitors, providing both cellular and molecular targets for the development of therapies for ACDMPV and other diffuse lung diseases of infancy.

Single-cell multiomic analysis identifies a HOX-PBX gene network regulating the survival of lymphangioleiomyomatosis cells.

Lymphangioleiomyomatosis (LAM) is a rare, progressive lung disease that predominantly affects women. LAM cells carry TSC1/TSC2 mutations, causing mTORC1 hyperactivation and uncontrolled cell growth. mTORC1 inhibitors stabilize lung function; however, sustained efficacy requires long-term administration, and some patients fail to tolerate or respond to therapy. Although the genetic basis of LAM is known, mechanisms underlying LAM pathogenesis remain elusive. We integrated single-cell RNA sequencing and single-nuclei ATAC-seq of LAM lungs to construct a gene regulatory network controlling the transcriptional program of LAM cells. We identified activation of uterine-specific HOX-PBX transcriptional programs in pulmonary LAMCORE cells as regulators of cell survival depending upon HOXD11-PBX1 dimerization. Accordingly, blockage of HOXD11-PBX1 dimerization by HXR9 suppressed LAM cell survival in vitro and in vivo. PBX1 regulated STAT1/3, increased the expression of antiapoptotic genes, and promoted LAM cell survival in vitro. The HOX-PBX gene network provides promising targets for treatment of LAM/TSC mTORC1-hyperactive cancers.

EMC3 regulates mesenchymal cell survival via control of the mitotic spindle assembly.

Eukaryotic cells transit through the cell cycle to produce two daughter cells. Dysregulation of the cell cycle leads to cell death or tumorigenesis. Herein, we found a subunit of the ER membrane complex, EMC3, as a key regulator of cell cycle. Conditional deletion of Emc3 in mouse embryonic mesoderm led to reduced size and patterning defects of multiple organs. Emc3 deficiency impaired cell proliferation, causing spindle assembly defects, chromosome mis-segregation, cell cycle arrest at G2/M, and apoptosis. Upon entry into mitosis, mesenchymal cells upregulate EMC3 protein levels and localize EMC3 to the mitotic centrosomes. Further analysis indicated that EMC3 works together with VCP to tightly regulate the levels and activity of Aurora A, an essential factor for centrosome function and mitotic spindle assembly: while overexpression of EMC3 or VCP degraded Aurora A, their loss led to increased Aurora A stability but reduced Aurora A phosphorylation in mitosis.

Molecular Function and Contribution of TBX4 in Development and Disease.

Over the past decade, recognition of the profound impact of the TBX4 (T-box 4) gene, which encodes a member of the evolutionarily conserved family of T-box-containing transcription factors, on respiratory diseases has emerged. The developmental importance of TBX4 is emphasized by the association of TBX4 variants with congenital disorders involving respiratory and skeletal structures; however, the exact role of TBX4 in human development remains incompletely understood. Here, we discuss the developmental, tissue-specific, and pathological TBX4 functions identified through human and animal studies and review the published TBX4 variants resulting in variable disease phenotypes. We also outline future research directions to fill the gaps in our understanding of TBX4 function and of how TBX4 disruption affects development.

PRDM3/16 Regulate Chromatin Accessibility Required for NKX2-1 Mediated Alveolar Epithelial Differentiation and Function.

Differential chromatin accessibility accompanies and mediates transcriptional control of diverse cell fates and their differentiation during embryogenesis. While the critical role of NKX2-1 and its transcriptional targets in lung morphogenesis and pulmonary epithelial cell differentiation is increasingly known, mechanisms by which chromatin accessibility alters the epigenetic landscape and how NKX2-1 interacts with other co-activators required for alveolar epithelial cell differentiation and function are not well understood. Here, we demonstrate that the paired domain zinc finger transcriptional regulators PRDM3 and PRDM16 regulate chromatin accessibility to mediate cell differentiation decisions during lung morphogenesis. Combined deletion of Prdm3 and Prdm16 in early lung endoderm caused perinatal lethality due to respiratory failure from loss of AT2 cell function. Prdm3/16 deletion led to the accumulation of partially differentiated AT1 cells and loss of AT2 cells. Combination of single cell RNA-seq, bulk ATAC-seq, and CUT&RUN demonstrated that PRDM3 and PRDM16 enhanced chromatin accessibility at NKX2-1 transcriptional targets in peripheral epithelial cells, all three factors binding together at a multitude of cell-type specific cis-active DNA elements. Network analysis demonstrated that PRDM3/16 regulated genes critical for perinatal AT2 cell differentiation, surfactant homeostasis, and innate host defense. Lineage specific deletion of PRDM3/16 in AT2 cells led to lineage infidelity, with PRDM3/16 null cells acquiring partial AT1 fate. Together, these data demonstrate that NKX2-1-dependent regulation of alveolar epithelial cell differentiation is mediated by epigenomic modulation via PRDM3/16.

LungMAP Portal Ecosystem: Systems-Level Exploration of the Lung.

An improved understanding of the human lung necessitates advanced systems models informed by an ever-increasing repertoire of molecular omics, cellular, imaging, and pathological datasets. To centralize and standardize information across broad lung research efforts we expanded the LungMAP.net website into a new gateway portal. This portal connects a broad spectrum of research networks, bulk and single-cell multi-omics data and a diverse collection of image data that span mammalian lung development, and disease. The data are standardized across species and technologies using harmonized data and metadata models that leverage recent advances including those from the Human Cell Atlas, diverse ontologies, and the LungMAP CellCards initiative. To cultivate future discoveries, we have aggregated a diverse collection of single-cell atlases for multiple species (human, rhesus, mouse), to enable consistent queries across technologies, cohorts, age, disease, and drug treatment. These atlases are provided as independent and integrated queryable datasets, with an emphasis on dynamic visualization, figure generation, re-analysis, cell-type curation, and automated reference-based classification of user-provided single-cell genomics datasets (Azimuth). As this resource grows, we intend to increase the breadth of available interactive interfaces, supported data types, data portals and datasets from LungMAP and external research efforts.

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.

Proteomic Analysis of Human Lung Development.

Rationale: The current understanding of human lung development derives mostly from animal studies. Although transcript-level studies have analyzed human donor tissue to identify genes expressed during normal human lung development, protein-level analysis that would enable the generation of new hypotheses on the processes involved in pulmonary development are lacking. Objectives: To define the temporal dynamic of protein expression during human lung development. Methods: We performed proteomics analysis of human lungs at 10 distinct times from birth to 8 years to identify the molecular networks mediating postnatal lung maturation. Measurements and Main Results: We identified 8,938 proteins providing a comprehensive view of the developing human lung proteome. The analysis of the data supports the existence of distinct molecular substages of alveolar development and predicted the age of independent human lung samples, and extensive remodeling of the lung proteome occurred during postnatal development. Evidence of post-transcriptional control was identified in early postnatal development. An extensive extracellular matrix remodeling was supported by changes in the proteome during alveologenesis. The concept of maturation of the immune system as an inherent part of normal lung development was substantiated by flow cytometry and transcriptomics. Conclusions: This study provides the first in-depth characterization of the human lung proteome during development, providing a unique proteomic resource freely accessible at Lungmap.net. The data support the extensive remodeling of the lung proteome during development, the existence of molecular substages of alveologenesis, and evidence of post-transcriptional control in early postnatal development.

YAP regulates alveolar epithelial cell differentiation and AGER via NFIB/KLF5/NKX2-1.

Ventilation is dependent upon pulmonary alveoli lined by two major epithelial cell types, alveolar type-1 (AT1) and 2 (AT2) cells. AT1 cells mediate gas exchange while AT2 cells synthesize and secrete pulmonary surfactants and serve as progenitor cells which repair the alveoli. We developed transgenic mice in which YAP was activated or deleted to determine its roles in alveolar epithelial cell differentiation. Postnatal YAP activation increased epithelial cell proliferation, increased AT1 cell numbers, and caused indeterminate differentiation of subsets of alveolar cells expressing atypical genes normally restricted to airway epithelial cells. YAP deletion increased expression of genes associated with mature AT2 cells. YAP activation enhanced DNA accessibility in promoters of transcription factors and motif enrichment analysis predicted target genes associated with alveolar cell differentiation. YAP participated with KLF5, NFIB, and NKX2-1 to regulate AGER. YAP plays a central role in a transcriptional network that regulates alveolar epithelial differentiation.

Therapeutic Potential of Endothelial Progenitor Cells in Pulmonary Diseases.

Compromised alveolar development and pulmonary vascular remodeling are hallmarks of pediatric lung diseases such as bronchopulmonary dysplasia (BPD) and alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV). Although advances in surfactant therapy, corticosteroids, and antiinflammatory drugs have improved clinical management of preterm infants, those who suffer with severe vascular complications still lack viable treatment options. Paucity of the alveolar capillary network in ACDMPV causes respiratory distress and leads to mortality in a vast majority of infants with ACDMPV. The discovery of endothelial progenitor cells (EPCs) in 1997 brought forth the paradigm of postnatal vasculogenesis and hope for promoting vascularization in fragile patient populations, such as those with BPD and ACDMPV. The identification of diverse EPC populations, both hematopoietic and nonhematopoietic in origin, provided a need to identify progenitor cell-selective markers that are linked to progenitor properties needed to develop cell-based therapies. Focusing on the future potential of EPCs for regenerative medicine, this review will discuss various aspects of EPC biology, beginning with the identification of hematopoietic, nonhematopoietic, and tissue-resident EPC populations. We will review knowledge related to cell surface markers, signature gene expression, and key transcriptional regulators and will explore the translational potential of EPCs for cell-based therapy for BPD and ACDMPV. The ability to produce pulmonary EPCs from patient-derived induced pluripotent stem cells in vitro holds promise for restoring vascular growth and function in the lungs of patients with pediatric pulmonary disorders.

Heterogeneity in Human Induced Pluripotent Stem Cell-derived Alveolar Epithelial Type II Cells Revealed with ABCA3/SFTPC Reporters.

Alveolar epithelial type 2 cells (AEC2s), the facultative progenitors of lung alveoli, are typically identified through the use of the canonical markers, SFTPC and ABCA3. Self-renewing AEC2-like cells have been generated from human induced pluripotent stem cells (iPSCs) through the use of knock-in SFTPC fluorochrome reporters. However, developmentally, SFTPC expression onset begins in the fetal distal lung bud tip and thus is not specific to mature AEC2s. Furthermore, SFTPC reporters appear to identify only those iPSC-derived AEC2s (iAEC2s) expressing the highest SFTPC levels. Here, we generate an ABCA3 knock-in GFP fusion reporter (ABCA3:GFP) that enables the purification of iAEC2s while allowing visualization of lamellar bodies, organelles associated with AEC2 maturation. Using an SFTPCtdTomato and ABCA3:GFP bifluorescent line for in vitro distal lung-directed differentiation, we observe later onset of ABCA3:GFP expression and broader identification of the subsequently emerging iAEC2 population based on ABCA3:GFP expression compared with SFTPCtdTomato expression. Comparing ABCA3:GFP/SFTPCtdTomato double-positive with ABCA3:GFP single-positive (SP) cells by RNA sequencing and functional studies reveals iAEC2 cellular heterogeneity with both populations functionally processing surfactant proteins but the SP cells exhibiting faster growth kinetics, increased clonogenicity, increased expression of progenitor markers, lower levels of SFTPC expression, and lower levels of AEC2 maturation markers. Over time, we observe that each population (double-positive and SP) gives rise to the other and each can serve as the parents of indefinitely self-renewing iAEC2 progeny. Our results indicate that iAEC2s are a heterogeneous population of cells with differing proliferation versus maturation properties, the majority of which can be tracked and purified using the ABCA3:GFP reporter or surrogate cell surface proteins, such as SLC34A2 and CPM.

VEGF receptor 2 (KDR) protects airways from mucus metaplasia through a Sox9-dependent pathway.

Mucus-secreting goblet cells are the dominant cell type in pulmonary diseases, e.g., asthma and cystic fibrosis (CF), leading to pathologic mucus metaplasia and airway obstruction. Cytokines including IL-13 are the major players in the transdifferentiation of club cells into goblet cells. Unexpectedly, we have uncovered a previously undescribed pathway promoting mucous metaplasia that involves VEGFa and its receptor KDR. Single-cell RNA sequencing analysis coupled with genetic mouse modeling demonstrates that loss of epithelial VEGFa, KDR, or MEK/ERK kinase promotes excessive club-to-goblet transdifferentiation during development and regeneration. Sox9 is required for goblet cell differentiation following Kdr inhibition in both mouse and human club cells. Significantly, airway mucous metaplasia in asthmatic and CF patients is also associated with reduced KDR signaling and increased SOX9 expression. Together, these findings reveal an unexpected role for VEGFa/KDR signaling in the defense against mucous metaplasia, offering a potential therapeutic target for this common airway pathology.