Lung Cell Atlases in Health and Disease.

The human lung cellular portfolio, traditionally characterized by cellular morphology and individual markers, is highly diverse, with over 40 cell types and a complex branching structure highly adapted for agile airflow and gas exchange. While constant during adulthood, lung cellular content changes in response to exposure, injury, and infection. Some changes are temporary, but others are persistent, leading to structural changes and progressive lung disease. The recent advance of single-cell profiling technologies allows an unprecedented level of detail and scale to cellular measurements, leading to the rise of comprehensive cell atlas styles of reporting. In this review, we chronical the rise of cell atlases and explore their contributions to human lung biology in health and disease.

Single Cell Analysis Reveals Novel Immune Perturbations in Fibrotic Hypersensitivity Pneumonitis.

Rationale: Fibrotic hypersensitivity pneumonitis is a debilitating interstitial lung disease driven by incompletely understood immune mechanisms. Objectives: To elucidate immune aberrations in fibrotic hypersensitivity pneumonitis in single-cell resolution. Methods: Single-cell 5' RNA sequencing was conducted on peripheral blood mononuclear cells and bronchoalveolar lavage cells obtained from 45 patients with fibrotic hypersensitivity pneumonitis, 63 idiopathic pulmonary fibrosis, 4 non-fibrotic hypersensitivity pneumonitis, and 36 healthy controls in the United States and Mexico. Analyses included differential gene expression (Seurat), transcription factor activity imputation (DoRothEA-VIPER), and trajectory analyses (Monocle3/Velocyto-scVelo-CellRank). Measurements and Main Results: Overall, 501,534 peripheral blood mononuclear cells from 110 patients and controls and 88,336 bronchoalveolar lavage cells from 19 patients were profiled. Compared to controls, fibrotic hypersensitivity pneumonitis has elevated classical monocytes (adjusted-p=2.5e-3) and are enriched in CCL3hi/CCL4hi and S100Ahi classical monocytes (adjusted-p<2.2e-16). Trajectory analyses demonstrate that S100Ahi classical monocytes differentiate into SPP1hi lung macrophages associated with fibrosis. Compared to both controls and idiopathic pulmonary fibrosis, fibrotic hypersensitivity pneumonitis patient cells are significantly enriched in GZMhi cytotoxic T cells. These cells exhibit transcription factor activities indicative of TGFß and TNFα/NFκB pathways. These results are publicly available at https://ildimmunecellatlas.org. Conclusions: Single-cell transcriptomics of fibrotic hypersensitivity pneumonitis patients uncovered novel immune perturbations, including previously undescribed increases in GZMhi cytotoxic CD4+ and CD8+ T cells - reflecting this disease's unique inflammatory T-cell driven nature - as well as increased S100Ahi and CCL3hi/CCL4hi classical monocytes also observed in idiopathic pulmonary fibrosis. Both cell populations may guide the development of new biomarkers and therapeutic interventions.

Microenvironmental sensing by fibroblasts controls macrophage population size.

Animal tissues comprise 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 coactivator 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.

Single-cell transcriptomic analysis of human pleura reveals stromal heterogeneity and informs in vitro models of mesothelioma.

The pleural lining of the thorax regulates local immunity, inflammation and repair. A variety of conditions, both benign and malignant, including pleural mesothelioma, can affect this tissue. A lack of knowledge concerning the mesothelial and stromal cells comprising the pleura has hampered the development of targeted therapies. Here, we present the first comprehensive single-cell transcriptomic atlas of the human parietal pleura and demonstrate its utility in elucidating pleural biology. We confirm the presence of known universal fibroblasts and describe novel, potentially pleural-specific, fibroblast subtypes. We also present transcriptomic characterisation of multiple in vitro models of benign and malignant mesothelial cells, and characterise these through comparison with in vivo transcriptomic data. While bulk pleural transcriptomes have been reported previously, this is the first study to provide resolution at the single-cell level. We expect our pleural cell atlas will prove invaluable to those studying pleural biology and disease. It has already enabled us to shed light on the transdifferentiation of mesothelial cells, allowing us to develop a simple method for prolonging mesothelial cell differentiation in vitro.

iDESC: identifying differential expression in single-cell RNA sequencing data with multiple subjects.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) technology has enabled assessment of transcriptome-wide changes at single-cell resolution. Due to the heterogeneity in environmental exposure and genetic background across subjects, subject effect contributes to the major source of variation in scRNA-seq data with multiple subjects, which severely confounds cell type specific differential expression (DE) analysis. Moreover, dropout events are prevalent in scRNA-seq data, leading to excessive number of zeroes in the data, which further aggravates the challenge in DE analysis. RESULTS: We developed iDESC to detect cell type specific DE genes between two groups of subjects in scRNA-seq data. iDESC uses a zero-inflated negative binomial mixed model to consider both subject effect and dropouts. The prevalence of dropout events (dropout rate) was demonstrated to be dependent on gene expression level, which is modeled by pooling information across genes. Subject effect is modeled as a random effect in the log-mean of the negative binomial component. We evaluated and compared the performance of iDESC with eleven existing DE analysis methods. Using simulated data, we demonstrated that iDESC had well-controlled type I error and higher power compared to the existing methods. Applications of those methods with well-controlled type I error to three real scRNA-seq datasets from the same tissue and disease showed that the results of iDESC achieved the best consistency between datasets and the best disease relevance. CONCLUSIONS: iDESC was able to achieve more accurate and robust DE analysis results by separating subject effect from disease effect with consideration of dropouts to identify DE genes, suggesting the importance of considering subject effect and dropouts in the DE analysis of scRNA-seq data with multiple subjects.

VISTA (PD-1H) Is a Crucial Immune Regulator to Limit Pulmonary Fibrosis.

VISTA (V domain immunoglobulin suppressor of T cell activation, also called PD-1H [programmed death-1 homolog]), a novel immune regulator expressed on myeloid and T lymphocyte lineages, is upregulated in mouse and human idiopathic pulmonary fibrosis (IPF). However, the significance of VISTA and its therapeutic potential in regulating IPF has yet to be defined. To determine the role of VISTA and its therapeutic potential in IPF, the expression profile of VISTA was evaluated from human single-cell RNA sequencing data (IPF Cell Atlas). Inflammatory response and lung fibrosis were assessed in bleomycin-induced experimental pulmonary fibrosis models in VISTA-deficient mice compared with wild-type littermates. In addition, these outcomes were evaluated after VISTA agonistic antibody treatment in the wild-type pulmonary fibrosis mice. VISTA expression was increased in lung tissue-infiltrating monocytes of patients with IPF. VISTA was induced in the myeloid population, mainly circulating monocyte-derived macrophages, during bleomycin-induced pulmonary fibrosis. Genetic ablation of VISTA drastically promoted pulmonary fibrosis, and bleomycin-induced fibroblast activation was dependent on the interaction between VISTA-expressing myeloid cells and fibroblasts. Treatment with VISTA agonistic antibody reduced fibrotic phenotypes accompanied by the suppression of lung innate immune and fibrotic mediators. In conclusion, these results suggest that VISTA upregulation in pulmonary fibrosis may be a compensatory mechanism to limit inflammation and fibrosis, and stimulation of VISTA signaling using VISTA agonists effectively limits the fibrotic innate immune landscape and consequent tissue fibrosis. Further studies are warranted to test VISTA as a novel therapeutic target for the IPF treatment.

Single-cell RNA-seq uncovers cellular heterogeneity and provides a signature for paediatric sleep apnoea.

BACKGROUND: Obstructive sleep apnoea (OSA) is a highly prevalent disease and a major cause of systemic inflammation leading to neurocognitive, behavioural, metabolic and cardiovascular dysfunction in children and adults. However, the impact of OSA on the heterogeneity of circulating immune cells remains to be determined. METHODS: We applied single-cell transcriptomics analysis (scRNA-seq) to identify OSA-induced changes in transcriptional landscape in peripheral blood mononuclear cell (PBMC) composition, which uncovered severity-dependent differences in several cell lineages. Furthermore, a machine-learning approach was used to combine scRNAs-seq cell-specific markers with those differentially expressed in OSA. RESULTS: scRNA-seq demonstrated OSA-induced heterogeneity in cellular composition and enabled the identification of previously undescribed cell types in PBMCs. We identified a molecular signature consisting of 32 genes, which distinguished OSA patients from various controls with high precision (area under the curve 0.96) and accuracy (93% positive predictive value and 95% negative predictive value) in an independent PBMC bulk RNA expression dataset. CONCLUSION: OSA deregulates systemic immune function and displays a molecular signature that can be assessed in standard cellular RNA without the need for pre-analytical cell separation, thereby making the assay amenable to application in a molecular diagnostic setting.

CD38 Mediates Lung Fibrosis by Promoting Alveolar Epithelial Cell Aging.

Rationale: A prevailing paradigm recognizes idiopathic pulmonary fibrosis (IPF) originating from various alveolar epithelial cell (AEC) injuries, and there is a growing appreciation of AEC aging as a key driver of the pathogenesis. Despite this progress, it is incompletely understood what main factor(s) contribute to the worsened alveolar epithelial aging in lung fibrosis. It remains a challenge how to dampen AEC aging and thereby mitigate the disease progression. Objectives: To determine the role of AEC CD38 (cluster of differentiation 38) in promoting cellular aging and lung fibrosis. Methods: We used single-cell RNA sequencing, real-time PCR, flow cytometry, and Western blotting. Measurements and Main Results: We discovered a pivotal role of CD38, a cardinal nicotinamide adenine dinucleotide (NAD) hydrolase, in AEC aging and its promotion of lung fibrosis. We found increased CD38 expression in IPF lungs that inversely correlated with the lung functions of patients. CD38 was primarily located in the AECs of human lung parenchyma and was markedly induced in IPF AECs. Similarly, CD38 expression was elevated in the AECs of fibrotic lungs of young mice and further augmented in those of old mice, which was in accordance with a worsened AEC aging phenotype and an aggravated lung fibrosis in the old animals. Mechanistically, we found that CD38 elevation downregulated intracellular NAD, which likely led to the aging promoting impairment of the NAD-dependent cellular and molecular activities. Furthermore, we demonstrated that genetic and pharmacological inactivation of CD38 improved these NAD dependent events and ameliorated bleomycin-induced lung fibrosis. Conclusions: Our study suggests targeting alveolar CD38 as a novel and effective therapeutic strategy to treat this pathology.

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

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 (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). 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.

Fibroblasts positive for meflin have anti-fibrotic properties in pulmonary fibrosis.

The prognosis of elderly individuals with idiopathic pulmonary fibrosis (IPF) remains poor. Fibroblastic foci, in which aggregates of proliferating fibroblasts and myofibroblasts are involved, are the pathological hallmark lesions in IPF to represent focal areas of active fibrogenesis. Fibroblast heterogeneity in fibrotic lesions hampers the discovery of the pathogenesis of pulmonary fibrosis. Therefore, to determine the pathogenesis of IPF, identification of functional fibroblasts is warranted. The aim of this study was to determine the role of fibroblasts positive for meflin, identified as a potential marker for mesenchymal stromal cells, during the development of pulmonary fibrosis.We characterised meflin-positive cells in a single-cell atlas established by single-cell RNA sequencing (scRNA-seq)-based profiling of 243 472 cells from 32 IPF lungs and 29 normal lung samples. We determined the role of fibroblasts positive for meflin using bleomycin (BLM)-induced pulmonary fibrosis.scRNA-seq combined with in situ RNA hybridisation identified proliferating fibroblasts positive for meflin in fibroblastic foci, not dense fibrosis, of fibrotic lungs in IPF patients. A BLM-induced lung fibrosis model for meflin-deficient mice showed that fibroblasts positive for meflin had anti-fibrotic properties to prevent pulmonary fibrosis. Although transforming growth factor-ß-induced fibrogenesis and cell senescence with the senescence-associated secretory phenotype were exacerbated in fibroblasts via the repression or lack of meflin, these were inhibited in meflin-deficient fibroblasts with meflin reconstitution.These findings provide evidence to show the biological importance of meflin expression on fibroblasts and myofibroblasts in the active fibrotic region of pulmonary fibrosis.

Transcriptomics of bronchoalveolar lavage cells identifies new molecular endotypes of sarcoidosis.

BACKGROUND: Sarcoidosis is a multisystem granulomatous disease of unknown origin with a variable and often unpredictable course and pattern of organ involvement. In this study we sought to identify specific bronchoalveolar lavage (BAL) cell gene expression patterns indicative of distinct disease phenotypic traits. METHODS: RNA sequencing by Ion Torrent Proton was performed on BAL cells obtained from 215 well-characterised patients with pulmonary sarcoidosis enrolled in the multicentre Genomic Research in Alpha-1 Antitrypsin Deficiency and Sarcoidosis (GRADS) study. Weighted gene co-expression network analysis and nonparametric statistics were used to analyse genome-wide BAL transcriptome. Validation of results was performed using a microarray expression dataset of an independent sarcoidosis cohort (Freiburg, Germany; n=50). RESULTS: Our supervised analysis found associations between distinct transcriptional programmes and major pulmonary phenotypic manifestations of sarcoidosis including T-helper type 1 (Th1) and Th17 pathways associated with hilar lymphadenopathy, transforming growth factor-ß1 (TGFB1) and mechanistic target of rapamycin (MTOR) signalling with parenchymal involvement, and interleukin (IL)-7 and IL-2 with airway involvement. Our unsupervised analysis revealed gene modules that uncovered four potential sarcoidosis endotypes including hilar lymphadenopathy with increased acute T-cell immune response; extraocular organ involvement with PI3K activation pathways; chronic and multiorgan disease with increased immune response pathways; and multiorgan involvement, with increased IL-1 and IL-18 immune and inflammatory responses. We validated the occurrence of these endotypes using gene expression, pulmonary function tests and cell differentials from Freiburg. CONCLUSION: Taken together, our results identify BAL gene expression programmes that characterise major pulmonary sarcoidosis phenotypes and suggest the presence of distinct disease molecular endotypes.

Single-Cell Transcriptional Archetypes of Airway Inflammation in Cystic Fibrosis.

Rationale: Cystic fibrosis (CF) is a life-shortening, multisystem hereditary disease caused by abnormal chloride transport. CF lung disease is driven by innate immune dysfunction and exaggerated inflammatory responses that contribute to tissue injury. To define the transcriptional profile of this airway immune dysfunction, we performed the first single-cell transcriptome characterization of CF sputum.Objectives: To define the transcriptional profile of sputum cells and its implication in the pathogenesis of immune function and the development of CF lung disease.Methods: We performed single-cell RNA sequencing of sputum cells from nine subjects with CF and five healthy control subjects. We applied novel computational approaches to define expression-based cell function and maturity profiles, herein called transcriptional archetypes.Measurements and Main Results: The airway immune cell repertoire shifted from alveolar macrophages in healthy control subjects to a predominance of recruited monocytes and neutrophils in CF. Recruited lung mononuclear phagocytes were abundant in CF and were separated into the following three archetypes: activated monocytes, monocyte-derived macrophages, and heat shock-activated monocytes. Neutrophils were the most prevalent in CF, with a dominant immature proinflammatory archetype. Although CF monocytes exhibited proinflammatory features, both monocytes and neutrophils showed transcriptional evidence of abnormal phagocytic and cell-survival programs.Conclusions: Our findings offer an opportunity to understand subject-specific immune dysfunction and its contribution to divergent clinical courses in CF. As we progress toward personalized applications of therapeutic and genomic developments, we hope this inflammation-profiling approach will enable further discoveries that change the natural history of CF lung disease.

Dedifferentiated early postnatal lung myofibroblasts redifferentiate in adult disease.

Alveolarization ensures sufficient lung surface area for gas exchange, and during bulk alveolarization in mice (postnatal day [P] 4.5-14.5), alpha-smooth muscle actin (SMA)+ myofibroblasts accumulate, secrete elastin, and lay down alveolar septum. Herein, we delineate the dynamics of the lineage of early postnatal SMA+ myofibroblasts during and after bulk alveolarization and in response to lung injury. SMA+ lung myofibroblasts first appear at ∼ P2.5 and proliferate robustly. Lineage tracing shows that, at P14.5 and over the next few days, the vast majority of SMA+ myofibroblasts downregulate smooth muscle cell markers and undergo apoptosis. Of note, ∼8% of these dedifferentiated cells and another ∼1% of SMA+ myofibroblasts persist to adulthood. Single cell RNA sequencing analysis of the persistent SMA- cells and SMA+ myofibroblasts in the adult lung reveals distinct gene expression profiles. For instance, dedifferentiated SMA- cells exhibit higher levels of tissue remodeling genes. Most interestingly, these dedifferentiated early postnatal myofibroblasts re-express SMA upon exposure of the adult lung to hypoxia or the pro-fibrotic drug bleomycin. However, unlike during alveolarization, these cells that re-express SMA do not proliferate with hypoxia. In sum, dedifferentiated early postnatal myofibroblasts are a previously undescribed cell type in the adult lung and redifferentiate in response to injury.

Noninvasive assessment of the lung inflammation-fibrosis axis by targeted imaging of CMKLR1.

Precision management of fibrotic lung diseases is challenging due to their diverse clinical trajectories and lack of reliable biomarkers for risk stratification and therapeutic monitoring. Here, we validated the accuracy of CMKLR1 as an imaging biomarker of the lung inflammation-fibrosis axis. By analyzing single-cell RNA sequencing datasets, we demonstrated CMKLR1 expression as a transient signature of monocyte-derived macrophages (MDMφ) enriched in patients with idiopathic pulmonary fibrosis (IPF). Consistently, we identified MDMφ as the major driver of the uptake of CMKLR1-targeting peptides in a murine model of bleomycin-induced lung fibrosis. Furthermore, CMKLR1-targeted positron emission tomography in the murine model enabled quantification and spatial mapping of inflamed lung regions infiltrated by CMKLR1-expressing macrophages and emerged as a robust predictor of subsequent lung fibrosis. Last, high CMKLR1 expression by bronchoalveolar lavage cells identified an inflammatory endotype of IPF with poor survival. Our investigation supports the potential of CMKLR1 as an imaging biomarker for endotyping and risk stratification of fibrotic lung diseases.

Dendritic Cell - Fibroblast Crosstalk via TLR9 and AHR Signaling Drives Lung Fibrogenesis.

Idiopathic pulmonary fibrosis (IPF) is characterized by progressive scarring and loss of lung function. With limited treatment options, patients succumb to the disease within 2-5 years. The molecular pathogenesis of IPF regarding the immunologic changes that occur is poorly understood. We characterize a role for non-canonical aryl-hydrocarbon receptor signaling (ncAHR) in dendritic cells (DCs) that leads to production of IL-6 and IL-17, promoting fibrosis. TLR9 signaling in myofibroblasts is shown to regulate production of TDO2 which converts tryptophan into the endogenous AHR ligand kynurenine. Mice with augmented ncAHR signaling were created by crossing floxed AHR exon-2 deletion mice (AHR Δex2 ) with mice harboring a CD11c-Cre. Bleomycin was used to study fibrotic pathogenesis. Isolated CD11c+ cells and primary fibroblasts were treated ex-vivo with relevant TLR agonists and AHR modulating compounds to study how AHR signaling influenced inflammatory cytokine production. Human datasets were also interrogated. Inhibition of all AHR signaling rescued fibrosis, however, AHR Δex2 mice treated with bleomycin developed more fibrosis and DCs from these mice were hyperinflammatory and profibrotic upon adoptive transfer. Treatment of fibrotic fibroblasts with TLR9 agonist increased expression of TDO2. Study of human samples corroborate the relevance of these findings in IPF patients. We also, for the first time, identify that AHR exon-2 floxed mice retain capacity for ncAHR signaling.

Rational engineering of lung alveolar epithelium.

Engineered whole lungs may one day expand therapeutic options for patients with end-stage lung disease. However, the feasibility of ex vivo lung regeneration remains limited by the inability to recapitulate mature, functional alveolar epithelium. Here, we modulate multimodal components of the alveolar epithelial type 2 cell (AEC2) niche in decellularized lung scaffolds in order to guide AEC2 behavior for epithelial regeneration. First, endothelial cells coordinate with fibroblasts, in the presence of soluble growth and maturation factors, to promote alveolar scaffold population with surfactant-secreting AEC2s. Subsequent withdrawal of Wnt and FGF agonism synergizes with tidal-magnitude mechanical strain to induce the differentiation of AEC2s to squamous type 1 AECs (AEC1s) in cultured alveoli, in situ. These results outline a rational strategy to engineer an epithelium of AEC2s and AEC1s contained within epithelial-mesenchymal-endothelial alveolar-like units, and highlight the critical interplay amongst cellular, biochemical, and mechanical niche cues within the reconstituting alveolus.

Emergence of division of labor in tissues through cell interactions and spatial cues.

Most cell types in multicellular organisms can perform multiple functions. However, not all functions can be optimally performed simultaneously by the same cells. Functions incompatible at the level of individual cells can be performed at the cell population level, where cells divide labor and specialize in different functions. Division of labor can arise due to instruction by tissue environment or through self-organization. Here, we develop a computational framework to investigate the contribution of these mechanisms to division of labor within a cell-type population. By optimizing collective cellular task performance under trade-offs, we find that distinguishable expression patterns can emerge from cell-cell interactions versus instructive signals. We propose a method to construct ligand-receptor networks between specialist cells and use it to infer division-of-labor mechanisms from single-cell RNA sequencing (RNA-seq) and spatial transcriptomics data of stromal, epithelial, and immune cells. Our framework can be used to characterize the complexity of cell interactions within tissues.

A hybrid machine learning and regression method for cell type deconvolution of spatial barcoding-based transcriptomic data.

Spatial barcoding-based transcriptomic (ST) data require cell type deconvolution for cellular-level downstream analysis. Here we present SDePER, a hybrid machine learning and regression method, to deconvolve ST data using reference single-cell RNA sequencing (scRNA-seq) data. SDePER uses a machine learning approach to remove the systematic difference between ST and scRNA-seq data (platform effects) explicitly and efficiently to ensure the linear relationship between ST data and cell type-specific expression profile. It also considers sparsity of cell types per capture spot and across-spots spatial correlation in cell type compositions. Based on the estimated cell type proportions, SDePER imputes cell type compositions and gene expression at unmeasured locations in a tissue map with enhanced resolution. Applications to coarse-grained simulated data and four real datasets showed that SDePER achieved more accurate and robust results than existing methods, suggesting the importance of considering platform effects, sparsity and spatial correlation in cell type deconvolution.

A Multi-omic Analysis of the Human Lung Reveals Distinct Cell Specific Aging and Senescence Molecular Programs.

Age is a major risk factor for lung disease. To understand the mechanisms underlying this association, we characterized the changing cellular, genomic, transcriptional, and epigenetic landscape of lung aging using bulk and single-cell RNAseq (scRNAseq) data. Our analysis revealed age-associated gene networks that reflected hallmarks of aging, including mitochondrial dysfunction, inflammation, and cellular senescence. Cell type deconvolution revealed age-associated changes in the cellular composition of the lung: decreased alveolar epithelial cells and increased fibroblasts and endothelial cells. In the alveolar microenvironment, aging is characterized by decreased AT2B cells and reduced surfactant production, a finding that was validated by scRNAseq and IHC. We showed that a previously reported senescence signature, SenMayo, captures cells expressing canonical senescence markers. SenMayo signature also identified cell-type specific senescence-associated co-expression modules that have distinct molecular functions, including ECM regulation, cell signaling, and damage response pathways. Analysis of somatic mutations showed that burden was highest in lymphocytes and endothelial cells and was associated with high expression of senescence signature. Finally, aging and senescence gene expression modules were associated with differentially methylated regions, with inflammatory markers such as IL1B, IL6R, and TNF being significantly regulated with age. Our findings provide new insights into the mechanisms underlying lung aging and may have implications for the development of interventions to prevent or treat age-related lung diseases.