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.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.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 (COL15A1neg) localized to the lung parenchyma and systemic-venous ECs (COL15A1pos) 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).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.

Abstract Post-acute infection syndromes may develop after acute viral disease 1 . Infection with SARS-CoV-2 can result in the development of a post-acute infection syndrome known as long COVID. Individuals with long COVID frequently report unremitting fatigue, post-exertional malaise, and a variety of cognitive and autonomic dysfunctions 2–4 . However, the biological processes that are associated with the development and persistence of these symptoms are unclear. Here 275 individuals with or without long COVID were enrolled in a cross-sectional study that included multidimensional immune phenotyping and unbiased machine learning methods to identify biological features associated with long COVID. Marked differences were noted in circulating myeloid and lymphocyte populations relative to the matched controls, as well as evidence of exaggerated humoral responses directed against SARS-CoV-2 among participants with long COVID. Furthermore, higher antibody responses directed against non-SARS-CoV-2 viral pathogens were observed among individuals with long COVID, particularly Epstein–Barr virus. Levels of soluble immune mediators and hormones varied among groups, with cortisol levels being lower among participants with long COVID. Integration of immune phenotyping data into unbiased machine learning models identified the key features that are most strongly associated with long COVID status. Collectively, these findings may help to guide future studies into the pathobiology of long COVID and help with developing relevant biomarkers.

Abstract Single-cell technologies have transformed our understanding of human tissues. Yet, studies typically capture only a limited number of donors and disagree on cell type definitions. Integrating many single-cell datasets can address these limitations of individual studies and capture the variability present in the population. Here we present the integrated Human Lung Cell Atlas (HLCA), combining 49 datasets of the human respiratory system into a single atlas spanning over 2.4 million cells from 486 individuals. The HLCA presents a consensus cell type re-annotation with matching marker genes, including annotations of rare and previously undescribed cell types. Leveraging the number and diversity of individuals in the HLCA, we identify gene modules that are associated with demographic covariates such as age, sex and body mass index, as well as gene modules changing expression along the proximal-to-distal axis of the bronchial tree. Mapping new data to the HLCA enables rapid data annotation and interpretation. Using the HLCA as a reference for the study of disease, we identify shared cell states across multiple lung diseases, including SPP1 + profibrotic monocyte-derived macrophages in COVID-19, pulmonary fibrosis and lung carcinoma. Overall, the HLCA serves as an example for the development and use of large-scale, cross-dataset organ atlases within the Human Cell Atlas.

Abstract In this study, we leveraged a combination of single cell RNAseq, cytometry by time of flight (CyTOF), and flow cytometry to study the biology of a unique macrophage population in pulmonary fibrosis. Using the profiling data from 312,928 cells derived from 32 idiopathic pulmonary fibrosis (IPF), 29 healthy control and 18 chronic obstructive pulmonary disease (COPD) lungs, we identified an expanded population of macrophages in IPF that have a unique transcriptional profile associated with pro-fibrotic signature. These macrophages attain a hybrid transitional state between alveolar and interstitial macrophages, are enriched with biological processes of pro-fibrotic immune cells, and express novel surface markers and genes that have not been previously reported. We then applied single cell CyTOF to simultaneously measure 37 markers to precisely phenotype the uniquely expanded macrophage subset in IPF lungs. The SPADE algorithm independently identified an expanded macrophage cluster, and validated CD84 and CD36 as novel surface markers that highly label this cluster. Using a separate validation cohort, we confirmed an increase in CD84 ++ CD36 ++ macrophage population in IPF compared to control and COPD lungs by flow cytometry. Further, using the signature from the IPF-specific macrophages and the LINCS drug database, we predicted small molecules that could reverse the signature of IPF-specific macrophages, and validated two molecules, CRT and Cucur, using THP-1 derived human macrophages and precision-cut lung slices (PCLS) from IPF patients. Utilizing a multi-dimensional translational approach, our work identified a novel and targetable population of macrophages found in end-stage pulmonary fibrosis. One Sentence Summary Single cell RNAseq, CyTOF, and flow cytometry reveal the presence of an aberrant macrophage population in pulmonary fibrosis

Abstract Single-cell RNA-sequencing data can revolutionize our understanding of the patterns of cell-cell and ligand-receptor connectivity that influence the function of tissues and organs. However, the quantification and visualization of these patterns are major computational and epistemological challenges. Here, we present Connectome , a software package for R which facilitates rapid calculation, and interactive exploration, of cell-cell signaling network topologies contained in single-cell RNA-sequencing data. Connectome can be used with any reference set of known ligand-receptor mechanisms. It has built-in functionality to facilitate differential and comparative connectomics, in which complete mechanistic networks are quantitatively compared between systems. Connectome includes computational and graphical tools designed to analyze and explore cell-cell connectivity patterns across disparate single-cell datasets. We present approaches to quantify these topologies and discuss some of the biologic theory leading to their design.

Abstract Summary Recent years have seen the release of several toolsets that reveal cell-cell interactions from single-cell data. However, all existing approaches leverage mean celltype gene expression values, and do not preserve the single-cell fidelity of the original data. Here, we present NICHES ( N iche I nteractions and C ommunication H eterogeneity in E xtracellular S ignaling), a tool to explore extracellular signaling at the truly single-cell level. NICHES allows embedding of ligand-receptor signal proxies to visualize heterogeneous signaling archetypes within cell clusters, between cell clusters, and across experimental conditions. When applied to spatial transcriptomic data, NICHES can be used to reflect local cellular microenvironment. NICHES can operate with any list of ligand-receptor signaling mechanisms and is compatible with existing single-cell packages and pseudotime techniques. NICHES is also a user friendly and extensible program, allowing rapid analysis of cell-cell signaling at single-cell resolution. Availability and implementation NICHES is an open-source software implemented in R under academic free license v3.0 and it is available at github.com/msraredon/NICHES. Use-case vignettes are available at https://msraredon.github.io/NICHES/ . Contact michasam.raredon@yale.edu ; yuval.kluger@yale.edu

ABSTRACT Idiopathic pulmonary fibrosis (IPF) is a fatal disease with limited treatment options. In this study we focus on the profibrotic properties of airway basal cells (ABC) obtained from patients with IPF (IPF-ABC). Single cell RNA sequencing of bronchial brushes revealed extensive reprogramming of IPF-ABC towards a KRT17 high PTEN low dedifferentiated cell type. In the 3D organoid model, compared to ABC obtained from healthy volunteers, IPF-ABC give rise to more bronchospheres, de novo bronchial structures resembling lung developmental processes, induce fibroblast proliferation and extracellular matrix deposition in co-culture. Intratracheal application of IPF-ABC into minimally injured lungs of Rag2 -/- or NRG mice causes severe fibrosis, remodeling of the alveolar compartment, and formation of honeycomb cyst-like structures. Connectivity MAP analysis of scRNA seq of bronchial brushings suggested that gene expression changes in IPF-ABC can be reversed by SRC inhibition. After demonstrating enhanced SRC expression and activity in these cells, and in IPF lungs, we tested the effects of saracatinib, a potent SRC inhibitor previously studied in humans. We demonstrated that saracatinib modified in-vitro and in-vivo the profibrotic changes observed in our 3D culture system and novel mouse xenograft model.

Abstract Animal tissues are comprised of diverse cell types. However, the mechanisms controlling the number of each cell type within tissue compartments remain poorly understood. Here, we report that different cell types utilize distinct strategies to control population numbers. Proliferation of fibroblasts, stromal cells important for tissue integrity, is limited by space availability. In contrast, proliferation of macrophages, innate immune cells involved in defense, repair, and homeostasis, is constrained by growth factor availability. Examination of density-dependent gene expression in fibroblasts revealed that Hippo and TGF- β target genes are both regulated by cell density. We found YAP1, the transcriptional co-activator of the Hippo signaling pathway, directly regulates expression of Csf1 , the lineage-specific growth factor for macrophages, through an enhancer of Csf1 that is specifically active in fibroblasts. Activation of YAP1 in fibroblasts elevates Csf1 expression and is sufficient to increase the number of macrophages at steady state. Our data also suggest that expression programs in fibroblasts that change with density may result from sensing of mechanical force through actin-dependent mechanisms. Altogether, we demonstrate that two different modes of population control are connected and coordinated to regulate cell numbers of distinct cell types. Sensing of the tissue environment may serve as a general strategy to control tissue composition. Significance Statement Collections of distinct cell types constitute animal tissues. To perform their unique functions, each cell type must exist in the correct number and proportion in a given tissue compartment. However, many of the mechanisms regulating and coordinating cell population sizes remain enigmatic. Our study characterizes two different modes of population size control, utilized by two ubiquitous cell types, macrophages and fibroblasts. Macrophage populations are more sensitive to the presence of growth factors in the environment and fibroblasts are more sensitive to space limitations. Intriguingly, space-sensing mechanisms in fibroblasts directly control the production of growth factor for macrophages and thus macrophage numbers. This link suggests a mechanism by which macrophage compartment size is controlled by stromal cells according to the microenvironment.

Abstract Background Despite its importance in health and disease, the cellular diversity of the lung endothelium has not been systematically characterized in humans. Here 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, both in health and disease. Methods We reprocessed control single cell RNA sequencing (scRNAseq) data from five datasets of whole lungs that were used for the analysis of pan-endothelial markers, we later included a sixth dataset of sorted control EC for the vascular subpopulation analysis. EC populations were characterized through iterative clustering with subsequent differential expression analysis. Marker genes were validated by immunohistochemistry and in situ hybridization. Signaling network between different lung cell types was studied using connectomic analysis. For cross species analysis we applied the same methods to scRNAseq data obtained from mouse lungs. Results The six lung scRNAseq datasets were reanalyzed and annotated to identify over 15,000 vascular EC cells from 73 individuals. Differential expression analysis of EC revealed signatures corresponding to endothelial lineage, including pan-endothelial, pan-vascular 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 two 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 EC we confirmed their subclassification into recently discovered aerocytes characterized by EDNRB, SOSTDC1 and TBX2 and general capillary EC. We confirmed that all six endothelial cell types, including the systemic-venous EC 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. Our manuscript is accompanied by an online data mining tool ( www.LungEndothelialCellAtlas.com ). Conclusion Our integrated analysis provides the comprehensive and well-crafted reference atlas of lung endothelial cells in the normal lung and confirms and describes in detail previously unrecognized endothelial populations across a large number of humans and mice.

Abstract While inhibition of T cell co-inhibitory receptors has revolutionized cancer therapy, the mechanisms governing their expression on human T cells have not been elucidated. Type 1 interferon (IFN-I) modulates T cell immunity in viral infection, autoimmunity, and cancer, and may facilitate induction of T cell exhaustion in chronic viral infection 1,2 . Here we show that IFN-I regulates co-inhibitory receptors expression on human T cells, inducing PD-1/TIM-3/LAG-3 while surprisingly inhibiting TIGIT expression. High-temporal-resolution mRNA profiling of IFN-I responses enabled the construction of dynamic transcriptional regulatory networks uncovering three temporal transcriptional waves. Perturbation of key transcription factors on human primary T cells revealed both canonical and non-canonical IFN-I transcriptional regulators, and identified unique regulators that control expression of co-inhibitory receptors. To provide direct in vivo evidence for the role of IFN-I on co-inhibitory receptors, we then performed single cell RNA-sequencing in subjects infected with SARS-CoV-2, where viral load was strongly associated with T cell IFN-I signatures. We found that the dynamic IFN-I response in vitro closely mirrored T cell features with acute IFN-I linked viral infection, with high LAG3 and decreased TIGIT expression. Finally, our gene regulatory network identified SP140 as a key regulator for differential LAG3 and TIGIT expression. The construction of co-inhibitory regulatory networks induced by IFN-I with identification of unique transcription factors controlling their expression may provide targets for enhancement of immunotherapy in cancer, infectious diseases, and autoimmunity.