Human lung single-cell atlas reveals the complexity and diversity of aberrant cellular populations in pulmonary fibrosis.

Background: Cellular diversity of the lung endothelium has not been systematically characterized in humans. We provide a reference atlas of human lung endothelial cells (ECs) to facilitate a better understanding of the phenotypic diversity and composition of cells comprising the lung endothelium. Methods: We reprocessed human control single-cell RNA sequencing (scRNAseq) data from 6 datasets. EC populations were characterized through iterative clustering with subsequent differential expression analysis. Marker genes were validated by fluorescent microscopy and in situ hybridization. scRNAseq of primary lung ECs cultured in vitro was performed. The signaling network between different lung cell types was studied. For cross-species analysis or disease relevance, we applied the same methods to scRNAseq data obtained from mouse lungs or from human lungs with pulmonary hypertension. Results: Six lung scRNAseq datasets were reanalyzed and annotated to identify >15 000 vascular EC cells from 73 individuals. Differential expression analysis of EC revealed signatures corresponding to endothelial lineage, including panendothelial, panvascular, and subpopulation-specific marker gene sets. Beyond the broad cellular categories of lymphatic, capillary, arterial, and venous ECs, we found previously indistinguishable subpopulations; among venous EC, we identified 2 previously indistinguishable populations: pulmonary–venous ECs (COL15A1 ^neg ) localized to the lung parenchyma and systemic–venous ECs (COL15A1 ^pos ) localized to the airways and the visceral pleura; among capillary ECs, we confirmed their subclassification into recently discovered aerocytes characterized by EDNRB , SOSTDC1 , and TBX2 and general capillary EC. We confirmed that all 6 endothelial cell types, including the systemic–venous ECs and aerocytes, are present in mice and identified endothelial marker genes conserved in humans and mice. Ligand-receptor connectome analysis revealed important homeostatic crosstalk of EC with other lung resident cell types. scRNAseq of commercially available primary lung ECs demonstrated a loss of their native lung phenotype in culture. scRNAseq revealed that endothelial diversity is maintained in pulmonary hypertension. Our article is accompanied by an online data mining tool ( www.LungEndothelialCellAtlas.com ). Conclusions: Our integrated analysis provides a comprehensive and well-crafted reference atlas of ECs in the normal lung and confirms and describes in detail previously unrecognized endothelial populations across a large number of humans and mice.

To develop a systems biology model of fibrosis progression within the human lung we performed RNA sequencing and microRNA analysis on 95 samples obtained from 10 idiopathic pulmonary fibrosis (IPF) and 6 control lungs. Extent of fibrosis in each sample was assessed by microCT-measured alveolar surface density (ASD) and confirmed by histology. Regulatory gene expression networks were identified using linear mixed-effect models and dynamic regulatory events miner (DREM). Differential gene expression analysis identified a core set of genes increased or decreased before fibrosis was histologically evident that continued to change with advanced fibrosis. DREM generated a systems biology model (www.sb.cs.cmu.edu/IPFReg) that identified progressively divergent gene expression tracks with microRNAs and transcription factors that specifically regulate mild or advanced fibrosis. We confirmed model predictions by demonstrating that expression of POU2AF1, previously unassociated with lung fibrosis but proposed by the model as regulator, is increased in B lymphocytes in IPF lungs and that POU2AF1-knockout mice were protected from bleomycin-induced lung fibrosis. Our results reveal distinct regulation of gene expression changes in IPF tissue that remained structurally normal compared with moderate or advanced fibrosis and suggest distinct regulatory mechanisms for each stage.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive, and often fatal disorder. Two U.S. Food and Drug Administration–approved antifibrotic drugs, nintedanib and pirfenidone, slow the rate of decline in lung function, but responses are variable and side effects are common. Objectives: Using an in silico data-driven approach, we identified a robust connection between the transcriptomic perturbations in IPF disease and those induced by saracatinib, a selective Src kinase inhibitor originally developed for oncological indications. Based on these observations, we hypothesized that saracatinib would be effective at attenuating pulmonary fibrosis. Methods: We investigated the antifibrotic efficacy of saracatinib relative to nintedanib and pirfenidone in three preclinical models: 1) in vitro in normal human lung fibroblasts; 2) in vivo in bleomycin and recombinant Ad-TGF-β (adenovirus transforming growth factor-β) murine models of pulmonary fibrosis; and 3) ex vivo in mice and human precision-cut lung slices from these two murine models as well as patients with IPF and healthy donors. Measurements and Main Results: In each model, the effectiveness of saracatinib in blocking fibrogenic responses was equal or superior to nintedanib and pirfenidone. Transcriptomic analyses of TGF-β–stimulated normal human lung fibroblasts identified specific gene sets associated with fibrosis, including epithelial–mesenchymal transition, TGF-β, and WNT signaling that was uniquely altered by saracatinib. Transcriptomic analysis of whole-lung extracts from the two animal models of pulmonary fibrosis revealed that saracatinib reverted many fibrogenic pathways, including epithelial–mesenchymal transition, immune responses, and extracellular matrix organization. Amelioration of fibrosis and inflammatory cascades in human precision-cut lung slices confirmed the potential therapeutic efficacy of saracatinib in human lung fibrosis. Conclusions: These studies identify novel Src-dependent fibrogenic pathways and support the study of the therapeutic effectiveness of saracatinib in IPF treatment.

Rationale: Changes in peripheral blood cell populations have been observed, but not detailed, at single-cell resolution in idiopathic pulmonary fibrosis (IPF). Objectives: We sought to provide an atlas of the changes in the peripheral immune system in stable and progressive IPF. Methods: Peripheral blood mononuclear cells (PBMCs) from patients with IPF and control subjects were profiled using 10× chromium 5′ single-cell RNA sequencing. Flow cytometry was used for validation. Protein concentrations of regulatory T cells (Tregs) and monocyte chemoattractants were measured in plasma and lung homogenates from patients with IPF and control subjects. Measurements and Main Results: Thirty-eight PBMC samples from 25 patients with IPF and 13 matched control subjects yielded 149,564 cells that segregated into 23 subpopulations. Classical monocytes were increased in patients with progressive and stable IPF compared with control subjects (32.1%, 25.2%, and 17.9%, respectively; P < 0.05). Total lymphocytes were decreased in patients with IPF versus control subjects and in progressive versus stable IPF (52.6% vs. 62.6%, P = 0.035). Tregs were increased in progressive versus stable IPF (1.8% vs. 1.1% of all PBMCs, P = 0.007), although not different than controls, and may be associated with decreased survival (P = 0.009 in Kaplan-Meier analysis; and P = 0.069 after adjusting for age, sex, and baseline FVC). Flow cytometry analysis confirmed this finding in an independent cohort of patients with IPF. The fraction of Tregs out of all T cells was also increased in two cohorts of lung single-cell RNA sequencing. CCL22 and CCL18, ligands for CCR4 and CCR8 Treg chemotaxis receptors, were increased in IPF. Conclusions: The single-cell atlas of the peripheral immune system in IPF reveals an outcome-predictive increase in classical monocytes and Tregs, as well as evidence for a lung–blood immune recruitment axis involving CCL7 (for classical monocytes) and CCL18/CCL22 (for Tregs).

Idiopathic pulmonary fibrosis (IPF) is a progressive and ultimately fatal disease. Recent findings have shown a marked metabolic reprogramming associated with changes in mitochondrial homeostasis and autophagy during pulmonary fibrosis. The microRNA-33 (miR-33) family of microRNAs (miRNAs) encoded within the introns of sterol regulatory element binding protein ( SREBP ) genes are master regulators of sterol and fatty acid (FA) metabolism. miR-33 controls macrophage immunometabolic response and enhances mitochondrial biogenesis, FA oxidation, and cholesterol efflux. Here, we show that miR-33 levels are increased in bronchoalveolar lavage (BAL) cells isolated from patients with IPF compared with healthy controls. We demonstrate that specific genetic ablation of miR-33 in macrophages protects against bleomycin-induced pulmonary fibrosis. The absence of miR-33 in macrophages improves mitochondrial homeostasis and increases autophagy while decreasing inflammatory response after bleomycin injury. Notably, pharmacological inhibition of miR-33 in macrophages via administration of anti–miR-33 peptide nucleic acids (PNA-33) attenuates fibrosis in different in vivo and ex vivo mice and human models of pulmonary fibrosis. These studies elucidate a major role of miR-33 in macrophages in the regulation of pulmonary fibrosis and uncover a potentially novel therapeutic approach to treat this disease.