The most common cause of death due to COVID-19 remains respiratory failure. Yet, our understanding of the precise cellular and molecular changes underlying lung alveolar damage is limited. Here, we integrate single cell transcriptomic data of COVID-19 and donor lung tissue with spatial transcriptomic data stratifying histopathological stages of diffuse alveolar damage. We identify changes in cellular composition across progressive damage, including waves of molecularly distinct macrophages and depletion of epithelial and endothelial populations. Predicted markers of pathological states identify immunoregulatory signatures, including IFN-alpha and metallothionein signatures in early damage, and fibrosis-related collagens in late damage. Furthermore, we predict a fibrinolytic shutdown via endothelial upregulation of SERPINE1/PAI-1. Cell-cell interaction analysis revealed macrophage-derived SPP1/osteopontin signalling as a key regulator during early steps of alveolar damage. These results provide a comprehensive, spatially resolved atlas of alveolar damage progression in COVID-19, highlighting the cellular mechanisms underlying pro-inflammatory and pro-fibrotic pathways in severe disease. Here the authors characterise the cellular and molecular progression of lung alveolar damage in severe COVID-19 patients using integrated histopathology and cell atlassing, pinpointing a role for macrophage SPP1 signalling to vasculature in this process.

The human vascular system, comprising endothelial cells (ECs) and mural cells, covers a vast surface area in the body, providing a critical interface between blood and tissue environments. Functional differences exist across specific vascular beds, but their molecular determinants across tissues remain largely unknown. In this study, we integrated single-cell transcriptomics data from 19 human organs and tissues and defined 42 vascular cell states from approximately 67,000 cells (62 donors), including angiotypic transitional signatures along the arterial endothelial axis from large to small caliber vessels. We also characterized organotypic populations, including splenic littoral and blood-brain barrier ECs, thus clarifying the molecular profiles of these important cell states. Interrogating endothelial-mural cell molecular crosstalk revealed angiotypic and organotypic communication pathways related to Notch, Wnt, retinoic acid, prostaglandin and cell adhesion signaling. Transcription factor network analysis revealed differential regulation of downstream target genes in tissue-specific modules, such as those of FOXF1 across multiple lung vascular subpopulations. Additionally, we make mechanistic inferences of vascular drug targets within different vascular beds. This open-access resource enhances our understanding of angiodiversity and organotypic molecular signatures in human vascular cells, and has therapeutic implications for vascular diseases across tissues.

Recurrent myocardial ischemia can lead to left ventricular (LV) dysfunction in patients with coronary artery disease (CAD). In this observational cohort study, we assessed for chronic metabolomic and transcriptomic adaptations within LV myocardium of patients undergoing coronary artery bypass grafting. During surgery, paired transmural LV biopsies were acquired on the beating heart from regions with and without evidence of inducible ischemia on preoperative stress perfusion cardiovascular magnetic resonance. From 33 patients, 63 biopsies were acquired, compared to analysis of LV samples from 11 donor hearts. The global myocardial adenosine triphosphate (ATP):adenosine diphosphate (ADP) ratio was reduced in patients with CAD as compared to donor LV tissue, with increased expression of oxidative phosphorylation (OXPHOS) genes encoding the electron transport chain complexes across multiple cell types. Paired analyses of biopsies obtained from LV segments with or without inducible ischemia revealed no significant difference in the ATP:ADP ratio, broader metabolic profile or expression of ventricular cardiomyocyte genes implicated in OXPHOS. Differential metabolite analysis suggested dysregulation of several intermediates in patients with reduced LV ejection fraction, including succinate. Overall, our results suggest that viable myocardium in patients with stable CAD has global alterations in bioenergetic and transcriptional profile without large regional differences between areas with or without inducible ischemia.

Common genetic variants across individuals modulate the cellular response to pathogens and are implicated in diverse immune pathologies, yet how they dynamically alter the response upon infection is not well understood. Here, we triggered antiviral responses in human fibroblasts from 68 healthy donors, and profiled tens of thousands of cells using single-cell RNA-sequencing. We developed GASPACHO (GAuSsian Processes for Association mapping leveraging Cell HeterOgeneity), a statistical approach designed to identify nonlinear dynamic genetic effects across transcriptional trajectories of cells. This approach identified 1,275 expression quantitative trait loci (local false discovery rate 10%) that manifested during the responses, many of which were colocalized with susceptibility loci identified by genome-wide association studies of infectious and autoimmune diseases, including the OAS1 splicing quantitative trait locus in a COVID-19 susceptibility locus. In summary, our analytical approach provides a unique framework for delineation of the genetic variants that shape a wide spectrum of transcriptional responses at single-cell resolution.

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system1. The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG+ and IgA+ plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs.