Spatially resolved multiomics of human cardiac niches.

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.

Mapping interindividual dynamics of innate immune response at single-cell resolution.

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.

Assessment of left ventricular tissue mitochondrial bioenergetics in patients with stable coronary artery disease.

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.

Pathogenic variants damage cell composition and single cell transcription in cardiomyopathies.

Pathogenic variants in genes that cause dilated cardiomyopathy (DCM) and arrhythmogenic cardiomyopathy (ACM) convey high risks for the development of heart failure through unknown mechanisms. Using single-nucleus RNA sequencing, we characterized the transcriptome of 880,000 nuclei from 18 control and 61 failing, nonischemic human hearts with pathogenic variants in DCM and ACM genes or idiopathic disease. We performed genotype-stratified analyses of the ventricular cell lineages and transcriptional states. The resultant DCM and ACM ventricular cell atlas demonstrated distinct right and left ventricular responses, highlighting genotype-associated pathways, intercellular interactions, and differential gene expression at single-cell resolution. Together, these data illuminate both shared and distinct cellular and molecular architectures of human heart failure and suggest candidate therapeutic targets.