Discovering the Immune Cell Landscape in IPF by Whole Lung Single Cell Transcriptomics

Single cell transcriptomes of 27 IPF, 13 COPD and 14 control lungs were established. Cells were clustered into major cell types, then the immune cells consisting of 116,057 myeloid, 4,742 B-, and 19,165 T-& NK-cells were analyzed and subclustered. Subclustering revealed IPF-associated unique gene expression profiles subpopulations that were specific to the IPF lung. For example, we could identify an IPF-specific macrophage subpopulation which showed a pro-fibrotic profile characterized by SPP1, CHIT1, CCL18, SPARC, LGMN, CTSK, and MMP9 gene expression. Antagonistically, NK- and cytotoxic T-cells revealed a significantly increased expression of the protease inhibitors A2M and PZP in the IPF lun

Genetic and geographic influence on phenotypic variation in European sarcoidosis patients

IntroductionSarcoidosis is a highly variable disease in terms of organ involvement, type of onset and course. Associations of genetic polymorphisms with sarcoidosis phenotypes have been observed and suggest genetic signatures.MethodsAfter obtaining a positive vote of the competent ethics committee we genotyped 1909 patients of the deeply phenotyped Genetic-Phenotype Relationship in Sarcoidosis (GenPhenReSa) cohort of 31 European centers in 12 countries with 116 potentially disease-relevant single-nucleotide polymorphisms (SNPs). Using a meta-analysis, we investigated the association of relevant phenotypes (acute vs. sub-acute onset, phenotypes of organ involvement, specific organ involvements, and specific symptoms) with genetic markers. Subgroups were built on the basis of geographical, clinical and hospital provision considerations.ResultsIn the meta-analysis of the full cohort, there was no significant genetic association with any considered phenotype after correcting for multiple testing. In the largest sub-cohort (Serbia), we confirmed the known association of acute onset with TNF and reported a new association of acute onset an HLA polymorphism. Multi-locus models with sets of three SNPs in different genes showed strong associations with the acute onset phenotype in Serbia and Lublin (Poland) demonstrating potential region-specific genetic links with clinical features, including recently described phenotypes of organ involvement.DiscussionThe observed associations between genetic variants and sarcoidosis phenotypes in subgroups suggest that gene–environment-interactions may influence the clinical phenotype. In addition, we show that two different sets of genetic variants are permissive for the same phenotype of acute disease only in two geographic subcohorts pointing to interactions of genetic signatures with different local environmental factors. Our results represent an important step towards understanding the genetic architecture of sarcoidosis

Integrated Single-Cell Atlas of Endothelial Cells of the Human Lung

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

Single-cell multi-omics reveals dyssynchrony of the innate and adaptive immune system in progressive COVID-19.

Dysregulated immune responses against the SARS-CoV-2 virus are instrumental in severe COVID-19. However, the immune signatures associated with immunopathology are poorly understood. Here we use multi-omics single-cell analysis to probe the dynamic immune responses in hospitalized patients with stable or progressive course of COVID-19, explore V(D)J repertoires, and assess the cellular effects of tocilizumab. Coordinated profiling of gene expression and cell lineage protein markers shows that S100

Chronic lung diseases are associated with gene expression programs favoring SARS-CoV-2 entry and severity.

Patients with chronic lung disease (CLD) have an increased risk for severe coronavirus disease-19 (COVID-19) and poor outcomes. Here, we analyze the transcriptomes of 611,398 single cells isolated from healthy and CLD lungs to identify molecular characteristics of lung cells that may account for worse COVID-19 outcomes in patients with chronic lung diseases. We observe a similar cellular distribution and relative expression of SARS-CoV-2 entry factors in control and CLD lungs. CLD AT2 cells express higher levels of genes linked directly to the efficiency of viral replication and the innate immune response. Additionally, we identify basal differences in inflammatory gene expression programs that highlight how CLD alters the inflammatory microenvironment encountered upon viral exposure to the peripheral lung. Our study indicates that CLD is accompanied by changes in cell-type-specific gene expression programs that prime the lung epithelium for and influence the innate and adaptive immune responses to SARS-CoV-2 infection

Gene coexpression networks reveal novel molecular endotypes in alpha-1 antitrypsin deficiency

Background Alpha-1 antitrypsin deficiency (AATD) is a genetic condition that causes early onset pulmonary emphysema and airways obstruction. The complete mechanisms via which AATD causes lung disease are not fully understood. To improve our understanding of the pathogenesis of AATD, we investigated gene expression profiles of bronchoalveolar lavage (BAL) and peripheral blood mononuclear cells (PBMCs) in AATD individuals. Methods We performed RNA-Seq on RNA extracted from matched BAL and PBMC samples isolated from 89 subjects enrolled in the Genomic Research in Alpha-1 Antitrypsin Deficiency and Sarcoidosis (GRADS) study. Subjects were stratified by genotype and augmentation therapy. Supervised and unsupervised differential gene expression analyses were performed using Weighted Gene Co-expression Network Analysis (WGCNA) to identify gene profiles associated with subjects’ clinical variables. The genes in the most significant WGCNA module were used to cluster AATD individuals. Gene validation was performed by NanoString nCounter Gene Expression Assay. Result We observed modest effects of AATD genotype and augmentation therapy on gene expression. When WGCNA was applied to BAL transcriptome, one gene module, ME31 (2312 genes), correlated with the highest number of clinical variables and was functionally enriched with numerous immune T-lymphocyte related pathways. This gene module identified two distinct clusters of AATD individuals with different disease severity and distinct PBMC gene expression patterns. Conclusions We successfully identified novel clusters of AATD individuals where severity correlated with increased immune response independent of individuals’ genotype and augmentation therapy. These findings may suggest the presence of previously unrecognised disease endotypes in AATD that associate with T-lymphocyte immunity and disease severity