Pro-inflammatory alveolar macrophages associated with allograft dysfunction after lung transplantation

Moshkelgosha, Sajad; Wilson, Gavin; Duong, Allen; Andrews, Tallulah; Berra, Gregory; Renaud-Picard, Benjamin; Keshavjee, Shaf; Martinu, Tereza; MacParland, Sonya; Yeung, Jonathan; Juvet, Stephen

Purpose Lung transplant (LT) recipients experience episodes of immune-mediated acute lung allograft dysfunction (ALAD). We have applied single-cell RNA sequencing (scRNAseq) to bronchoalveolar lavage (BAL) cells of stable and ALAD patients to determine key cellular elements in dysfunctional lung allografts. Our particular focus here is on studying alveolar macrophages (AMs) as scRNAseq enables us to elucidate their heterogeneity and possible association with ALAD where our knowledge from cytometry-based assays is very limited. Methods Fresh bronchoalveolar lavage (BAL) cells from 6 LT patients, 3 with stable lung function (3044 ± 1519 cells) and 3 undergoing an episode of ALAD (2593 ± 904 cells) were used for scRNAseq. R Bioconductor and Seurat were used to perform QC, dimensionality reduction, annotation, pathway analysis, and trajectory. Donor and recipient deconvolution was performed using single nucleotide variations. Results Our data revealed that AMs are highly heterogeneous (12 transcriptionally distinct subsets in stable). We identified two AM subsets uniquely represented in ALAD. Based on pathway analysis and the top differentially expressed genes in BAL we annotated them as pro-inflammatory interferon-stimulated genes (ISG) and metallothioneins-mediated inflammatory (MT). Pseudotime analysis suggested that ISG AMs represent an earlier stage of differentiation which may suggest them as monocyte drive macrophages. Our functional analysis on an independent set of BAL samples shows that ALAD samples have significantly higher expression of CXCL10, a marker of ISG AM, as we as higher secretion of pro-inflammatory cytokines. Single nucleotide variation calling algorithm has allowed us to identify macrophages of donor origin and demonstrated that donor AMs are lost with time post-transplant. Conclusion Using scRNAseq, we observed AMs heterogeneity and identified specific subsets that may be associated with allograft dysfunction. Further exploration with scRNAseq will shed light on LT immunobiology and the role of AMs in allograft injury and dysfunction.