Pathogenic variants damage cell composition and single cell transcription in cardiomyopathies

Functional diversity of cardiac macrophages in health and disease

Macrophages make up a substantial portion of the stromal compartment of the heart in health and disease. In the past decade, the origins of these cardiac macrophages have been established as two broad populations derived from either embryonic or definitive haematopoiesis and that can be distinguished by the expression of CC-motif chemokine receptor 2 (CCR2). These cardiac macrophage populations are transcriptionally distinct and have differing cell surface markers and divergent roles in cardiac homeostasis and disease. Embryonic-derived CCR2- macrophages are a tissue-resident population that participates in tissue development, repair and maintenance, whereas CCR2+ macrophages are derived from definitive haematopoiesis and contribute to inflammation and tissue damage. Studies from the past 5 years have leveraged single-cell RNA sequencing technologies to expand our understanding of cardiac macrophage diversity, particularly of the monocyte-derived macrophage populations that reside in the injured and diseased heart. Emerging technologies in spatial transcriptomics have enabled the identification of distinct disease-associated cellular neighbourhoods consisting of macrophages, other immune cells and fibroblasts, highlighting the involvement of macrophages in cell-cell communication. Together, these discoveries lend new insights into the role of specific macrophage populations in the pathogenesis of cardiac disease, which can pave the way for the identification of new therapeutic targets and the development of diagnostic tools. In this Review, we discuss the developmental origin of cardiac macrophages and describe newly identified cell states and associated cellular neighbourhoods in the steady state and injury settings. We also discuss various contributions and effector functions of cardiac macrophages in homeostasis and disease.

Early genetic screening and cardiac intervention in patients with cardiomyopathies in a multidisciplinary clinic

AIMS

Patients with cardiomyopathies are a heterogeneous group of patients who experience high morbidity and mortality. Early cardiac assessment and intervention with access to genetic counselling in a multidisciplinary Cardiomyopathy Clinic may improve outcomes and prevent progression to advanced heart failure.

METHODS AND RESULTS

Our prospective cohort study was conducted at a multidisciplinary Cardiomyopathy Clinic with 421 patients enrolled (42.5% female, median age 58 years), including 224 patients with dilated cardiomyopathy (DCM, 42.9% female, median age 57 years), 72 with hypertrophic cardiomyopathy (HCM, 43.1% female, median age 60 years), 79 with infiltrative cardiomyopathy (65.8% female, median age 70 years) and 46 who were stage A/at risk for genetic cardiomyopathy (54.3% female, median age 36 years). Patients were seen in follow-up at a median of 18 months. A pathogenic/likely pathogenic variant was identified in 28.5% of the total cohort, including 33.3% of the DCM cohort (28% TTN mutations) and 34.1% of the HCM cohort (60% MYBPC3 and 20% MYH7) who underwent genetic testing. The use of angiotensin-converting enzyme inhibitors/angiotensin receptor blockers/angiotensin receptor neprilysin inhibitor (48.3-69.5% of total cohort, P < 0.001), β-blockers (58.4-72.4%, P < 0.001), mineralocorticoid receptor antagonists (33.9-41.4%, P = 0.0014) and sodium/glucose cotransporter-2 inhibitors (5.3-27.9%, P < 0.001) all increased at follow-up. Precision-based therapies were also implemented, including tafamidis for transthyretin amyloidosis (n = 21), enzyme replacement therapy for Fabry disease (n = 14) and mavacamten (n = 4) for HCM. Optimization of medications and devices resulted in improvements in left ventricular ejection fraction (LVEF) from 27% to 43% at follow-up for DCM patients with reduced LVEF at baseline (P < 0.001) and reduction in left ventricular mass index (LVMI) from 156 g/m2 to 128 g/m2 at follow-up for HCM patients with abnormal LVMI at baseline (P = 0.009). Optimization of therapies was associated with stable plasma biomarkers in stage B patients while lowering levels of BNP (619-517.5 pg/mL, P = 0.048), NT-proBNP (777.5-356 ng/L, P < 0.001) and hsTropT (31-22 ng/L, P = 0.005) at follow-up relative to baseline values for stage C patients. Despite stage B patients having overt cardiomyopathy at baseline, stage A and B patients had a similarly high probability of survival (χ2 = 0.204, P = 0.652). The overall cardiovascular mortality rate was low at 1.7% for the cohort (0.5% for stage B and 3.3% for stage C) over a median of 34-month follow-up.

CONCLUSION

Our study demonstrates that a multidisciplinary cardiomyopathy clinic can improve the clinical profiles of patients with diverse genetic cardiomyopathies.

Fibroblast-specific TGF-β signaling mediates cardiac dysfunction, fibrosis, and hypertrophy in obese diabetic mice

AIMS

Transforming growth factor (TGF)-β is up-regulated in the diabetic myocardium and may mediate fibroblast activation. We aimed at examining the role of TGF-β-induced fibroblast activation in the pathogenesis of diabetic cardiomyopathy.

METHODS AND RESULTS

We generated lean and obese db/db mice with fibroblast-specific loss of TbR2, the Type 2 receptor-mediating signaling through all three TGF-β isoforms, and mice with fibroblast-specific Smad3 disruption. Systolic and diastolic function, myocardial fibrosis, and hypertrophy were assessed. Transcriptomic studies and in vitro experiments were used to dissect mechanisms of fibroblast activation. Fibroblast-specific TbR2 loss attenuated systolic and diastolic dysfunction in db/db mice. The protective effects of fibroblast TbR2 loss in db/db mice were associated with attenuated fibrosis and reduced cardiomyocyte hypertrophy, suggesting that in addition to their role in fibrous tissue deposition, TGF-β-stimulated fibroblasts may also exert paracrine actions on cardiomyocytes. Fibroblast-specific Smad3 loss phenocopied the protective effects of fibroblast TbR2 loss in db/db mice. Db/db fibroblasts had increased expression of genes associated with oxidative response (such as Fmo2, encoding flavin-containing monooxygenase 2), matricellular genes (such as Thbs4 and Fbln2), and Lox (encoding lysyl oxidase). Ingenuity pathway analysis (IPA) predicted that neurohumoral mediators, cytokines, and growth factors (such as AGT, TGFB1, and TNF) may serve as important upstream regulators of the transcriptomic profile of diabetic mouse fibroblasts. IPA of scRNA-seq data identified TGFB1, p53, MYC, PDGF-BB, EGFR, and WNT3A/CTNNB1 as important upstream regulators underlying fibroblast activation in db/db hearts. Comparison of the transcriptome of fibroblasts from db/db mice with fibroblast-specific Smad3 loss and db/db Smad3 fl/fl controls identified Thbs4 [encoding thrombospondin-4 (TSP-4), a marker of activated fibroblasts] as a candidate diabetes-induced fibrogenic mediator. However, in vitro experiments showed no significant activating effects of matricellular or intracellular TSP-4 on cardiac fibroblasts.

CONCLUSION

Fibroblast-specific TGF-β/Smad3 signaling mediates ventricular fibrosis, hypertrophy, and dysfunction in Type 2 diabetes.

Myocardial fibrosis from the perspective of the extracellular matrix: Mechanisms to clinical impact

Fibrosis is defined by the excessive accumulation of extracellular matrix (ECM) and constitutes a central pathophysiological process that underlies tissue dysfunction, across organs, in multiple chronic diseases and during aging. Myocardial fibrosis is a key contributor to dysfunction and failure in numerous diseases of the heart and is a strong predictor of poor clinical outcome and mortality. The excess structural and matricellular ECM proteins deposited by cardiac fibroblasts, is found between cardiomyocytes (interstitial fibrosis), in focal areas where cardiomyocytes have died (replacement fibrosis), and around vessels (perivascular fibrosis). Although myocardial fibrosis has important clinical prognostic value, access to cardiac tissue biopsies for histological evaluation is limited. Despite challenges with sensitivity and specificity, cardiac magnetic resonance imaging (CMR) is the most applicable diagnostic tool in the clinic, and the scientific community is currently actively searching for blood biomarkers reflecting myocardial fibrosis, to complement the imaging techniques. The lack of mechanistic insights into specific pro- and anti-fibrotic molecular pathways has hampered the development of effective treatments to prevent or reverse myocardial fibrosis. Development and implementation of anti-fibrotic therapies is expected to improve patient outcomes and is an urgent medical need. Here, we discuss the importance of the ECM in the heart, the central role of fibrosis in heart disease, and mechanistic pathways likely to impact clinical practice with regards to diagnostics of myocardial fibrosis, risk stratification of patients, and anti-fibrotic therapy.

Genome-wide association analysis provides insights into the molecular etiology of dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a leading cause of heart failure and cardiac transplantation. We report a genome-wide association study and multi-trait analysis of DCM (14,256 cases) and three left ventricular traits (36,203 UK Biobank participants). We identified 80 genomic risk loci and prioritized 62 putative effector genes, including several with rare variant DCM associations (MAP3K7, NEDD4L and SSPN). Using single-nucleus transcriptomics, we identify cellular states, biological pathways, and intracellular communications that drive pathogenesis. We demonstrate that polygenic scores predict DCM in the general population and modify penetrance in carriers of rare DCM variants. Our findings may inform the design of genetic testing strategies that incorporate polygenic background. They also provide insights into the molecular etiology of DCM that may facilitate the development of targeted therapeutics.

An organotypic atlas of human vascular cells

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.

Single-cell transcriptomics reveal distinctive patterns of fibroblast activation in heart failure with preserved ejection fraction

Inflammation, fibrosis and metabolic stress critically promote heart failure with preserved ejection fraction (HFpEF). Exposure to high-fat diet and nitric oxide synthase inhibitor N[w]-nitro-l-arginine methyl ester (L-NAME) recapitulate features of HFpEF in mice. To identify disease-specific traits during adverse remodeling, we profiled interstitial cells in early murine HFpEF using single-cell RNAseq (scRNAseq). Diastolic dysfunction and perivascular fibrosis were accompanied by an activation of cardiac fibroblast and macrophage subsets. Integration of fibroblasts from HFpEF with two murine models for heart failure with reduced ejection fraction (HFrEF) identified a catalog of conserved fibroblast phenotypes across mouse models. Moreover, HFpEF-specific characteristics included induced metabolic, hypoxic and inflammatory transcription factors and pathways, including enhanced expression of Angiopoietin-like 4 (Angptl4) next to basement membrane compounds, such as collagen IV (Col4a1). Fibroblast activation was further dissected into transcriptional and compositional shifts and thereby highly responsive cell states for each HF model were identified. In contrast to HFrEF, where myofibroblast and matrifibrocyte activation were crucial features, we found that these cell states played a subsidiary role in early HFpEF. These disease-specific fibroblast signatures were corroborated in human myocardial bulk transcriptomes. Furthermore, we identified a potential cross-talk between macrophages and fibroblasts via SPP1 and TNFɑ with estimated fibroblast target genes including Col4a1 and Angptl4. Treatment with recombinant ANGPTL4 ameliorated the murine HFpEF phenotype and diastolic dysfunction by reducing collagen IV deposition from fibroblasts in vivo and in vitro. In line, ANGPTL4, was elevated in plasma samples of HFpEF patients and particularly high levels associated with a preserved global-longitudinal strain. Taken together, our study provides a comprehensive characterization of molecular fibroblast activation patterns in murine HFpEF, as well as the identification of Angiopoietin-like 4 as central mechanistic regulator with protective effects.

Genome-wide association study reveals mechanisms underlying dilated cardiomyopathy and myocardial resilience

Dilated cardiomyopathy (DCM) is a heart muscle disease that represents an important cause of morbidity and mortality, yet causal mechanisms remain largely elusive. Here, we perform a large-scale genome-wide association study and multitrait analysis for DCM using 9,365 cases and 946,368 controls. We identify 70 genome-wide significant loci, which show broad replication in independent samples and map to 63 prioritized genes. Tissue, cell type and pathway enrichment analyses highlight the central role of the cardiomyocyte and contractile apparatus in DCM pathogenesis. Polygenic risk scores constructed from our genome-wide association study predict DCM across different ancestry groups, show differing contributions to DCM depending on rare pathogenic variant status and associate with systolic heart failure across various clinical settings. Mendelian randomization analyses reveal actionable potential causes of DCM, including higher bodyweight and higher systolic blood pressure. Our findings provide insights into the genetic architecture and mechanisms underlying DCM and myocardial function more broadly.

Screening for Heart Failure: Biomarkers to Detect Heightened Risk in the General Population

PURPOSE OF REVIEW

Heart failure (HF) represents a growing global burden of morbidity and mortality. Identifying individuals at risk for HF development is increasingly important, particularly given the advent of disease-modifying therapies for HF as well as its major risk factors such as obesity actalnd diabetes. We aim to review the key circulating biomarkers associated with future HF which may contribute to HF risk prediction.

RECENT FINDINGS

While current guidelines recommend the use of natriuretic peptides and cardiac troponins in HF risk stratification, there are a diverse array of other emerging protein, metabolic, transcriptomic, and genomic biomarkers of future HF development. These biomarkers not only lend insight into the underlying pathophysiology of HF, which spans inflammation to cardiac fibrosis, but also offer an opportunity to further refine HF risk in addition to established biomarkers. As evolving techniques in molecular biology enable an increased understanding of the complex biologic contributions to HF pathophysiology, there is an important opportunity to construct integrated clinical and multi-omic models to best capture HF risk. Moving forward, future studies should seek to understand the contributions of sex differences, underlying comorbidity burden, and HF subtypes to an individual's HF risk. Further studies are necessary to fully define the clinical utility of biomarker screening approaches to refine HF risk assessment, as well as to link risk assessment directly to preventive strategies for HF.

Preoperative methionine restriction induces perivascular adipose tissue browning and improves vein graft remodeling in male mice

Short-term preoperative methionine restriction (MetR) is a promising translatable strategy to mitigate surgical injury response. However, its application to improve post-interventional vascular remodeling remains underexplored. Here we find that MetR protects from arterial intimal hyperplasia in a focal stenosis model and pathologic vascular remodeling following vein graft surgery in male mice. RNA sequencing reveals that MetR enhances browning in arterial (thoracic aorta) perivascular adipose tissue (PVAT) and induces it in venous (caval vein) PVAT. Specifically, Ppara is highly upregulated in PVAT-adipocytes upon MetR. Furthermore, MetR dampens the postoperative pro-inflammatory response to surgery in PVAT-macrophages in vivo and in vitro. This study shows that the detrimental effects of dysfunctional PVAT on vascular remodeling can be reversed by MetR, and identifies pathways involved in MetR-induced browning of PVAT. Furthermore, we demonstrate the potential of short-term preoperative MetR as a simple intervention to ameliorate vascular remodeling after vascular surgery.

Beyond the Heartbeat: Single-Cell Omics Redefining Cardiovascular Research

PURPOSE OF REVIEW

This review aims to explore recent advances in single-cell omics techniques as applied to various regions of the human heart, illuminating cellular diversity, regulatory networks, and disease mechanisms. We examine the contributions of single-cell transcriptomics, genomics, proteomics, epigenomics, and spatial transcriptomics in unraveling the complexity of cardiac tissues.

RECENT FINDINGS

Recent strides in single-cell omics technologies have revolutionized our understanding of the heart's cellular composition, cell type heterogeneity, and molecular dynamics. These advancements have elucidated pathological conditions as well as the cellular landscape in heart development. We highlight emerging applications of integrated single-cell omics, particularly for cardiac regeneration, disease modeling, and precision medicine, and emphasize the transformative potential of these technologies to advance cardiovascular research and clinical practice.

Targeting immune-fibroblast cell communication in heart failure

Inflammation and tissue fibrosis co-exist and are causally linked to organ dysfunction1,2. However, the molecular mechanisms driving immune-fibroblast cell communication in human cardiac disease remain unexplored and there are at present no approved treatments that directly target cardiac fibrosis3,4. Here we performed multiomic single-cell gene expression, epitope mapping and chromatin accessibility profiling in 45 healthy donor, acutely infarcted and chronically failing human hearts. We identified a disease-associated fibroblast trajectory that diverged into distinct populations reminiscent of myofibroblasts and matrifibrocytes, the latter expressing fibroblast activator protein (FAP) and periostin (POSTN). Genetic lineage tracing of FAP+ fibroblasts in vivo showed that they contribute to the POSTN lineage but not the myofibroblast lineage. We assessed the applicability of experimental systems to model cardiac fibroblasts and demonstrated that three different in vivo mouse models of cardiac injury were superior compared with cultured human heart and dermal fibroblasts in recapitulating the human disease phenotype. Ligand-receptor analysis and spatial transcriptomics predicted that interactions between C-C chemokine receptor type 2 (CCR2) macrophages and fibroblasts mediated by interleukin-1β (IL-1β) signalling drove the emergence of FAP/POSTN fibroblasts within spatially defined niches. In vivo, we deleted the IL-1 receptor on fibroblasts and the IL-1β ligand in CCR2+ monocytes and macrophages, and inhibited IL-1β signalling using a monoclonal antibody, and showed reduced FAP/POSTN fibroblasts, diminished myocardial fibrosis and improved cardiac function. These findings highlight the broader therapeutic potential of targeting inflammation to treat tissue fibrosis and preserve organ function.

SnRNA-seq reveals differential functional transcriptional pathway alterations in three mutant types of dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a leading cause of heart failure, characterized by ventricular dilation, thinning of the ventricular walls, and systolic dysfunction in either the left or both ventricles, often accompanied by fibrosis. Human cardiac tissue is composed of various cell types, including cardiomyocytes (CMs), fibroblasts (FBs), endothelial cells (ECs), macrophages, lymphocytes and so on. In DCM patients, these cells frequently undergo functional and phenotypic changes, contributing to contractile dysfunction, inflammation, fibrosis, and cell death, thereby increasing the risk of heart failure. This study focuses on DCM patients with mutations (LMNA, RBM20, and TTN) and analyzes functional changes in subpopulations of four cardiac cell types. The study involves functional annotation of subpopulations within each cell type and explores the association between gene mutations and specific functions and pathways. Additionally, the SCENIC method is employed of a particular cell subpopulation with significant functional importance, aiming to identify key transcriptional regulators in specific cell states. By analyzing the expression levels of ligand-receptor pairs in vCM4, vFB2, EC5.0, T cells, and NK cells across the DCM mutant genotypes, we predicted their signaling pathways and communications. This research provides insights into the molecular mechanisms of DCM and potential therapeutic targets.

Evolution of translational control and the emergence of genes and open reading frames in human and non-human primate hearts

Evolutionary innovations can be driven by changes in the rates of RNA translation and the emergence of new genes and small open reading frames (sORFs). In this study, we characterized the transcriptional and translational landscape of the hearts of four primate and two rodent species through integrative ribosome and transcriptomic profiling, including adult left ventricle tissues and induced pluripotent stem cell-derived cardiomyocyte cell cultures. We show here that the translational efficiencies of subunits of the mitochondrial oxidative phosphorylation chain complexes IV and V evolved rapidly across mammalian evolution. Moreover, we discovered hundreds of species-specific and lineage-specific genomic innovations that emerged during primate evolution in the heart, including 551 genes, 504 sORFs and 76 evolutionarily conserved genes displaying human-specific cardiac-enriched expression. Overall, our work describes the evolutionary processes and mechanisms that have shaped cardiac transcription and translation in recent primate evolution and sheds light on how these can contribute to cardiac development and disease.

Transcriptomic Profiling of Human Myocardium at Sudden Death to Define Vulnerable Substrate for Lethal Arrhythmias

BACKGROUND

While some chronic pathological substrates for sudden cardiac death (SCD) are well known (eg, coronary artery disease and left ventricular [LV] dysfunction), the acute vulnerable myocardial state predisposing to fatal arrhythmia remains a critical barrier to near-term SCD prevention.

OBJECTIVES

This study sought to define the distinct myocardial transcriptomic profile of autopsy-defined arrhythmic sudden deaths, compared to nonarrhythmic sudden deaths and trauma deaths, to determine the acute vulnerable state in the hours to days before SCD.

METHODS

We used autopsy to adjudicate arrhythmic from nonarrhythmic causes in 1,265 sudden deaths in San Francisco County from 2011 to 2018. We performed a degradation-tolerant transcriptomic evaluation of LVs sampled at the time of SCD from 245 consented cases using a curated panel of 448 gene probes with known or hypothesized association with SCD.

RESULTS

The targeted transcriptome of arrhythmic (n = 129) vs nonarrhythmic (n = 90 nonarrhythmic sudden deaths + 26 trauma deaths) LV samples revealed 31 differentially up-regulated and 36 down-regulated genes (adjusted P < 0.05) related to the collagen-containing extracellular matrix (up-regulation of FAP, FMOD, and LTBP2), regulation of ion transport (up-regulation of KCNA5 and KCNN3 and down-regulation of KCNJ8, KCNK1, and KCNJ5), and contraction (down-regulation of MYH6). Fibrosis-related genes showed the highest magnitude increased expression in arrhythmic vs nonarrhythmic deaths and vs published transcriptomes from end-stage heart failure. After molecular stratification by known markers for mature (COL1A1, COL1A2, COL3A1) and active (POSTN, MEOX1) fibrosis, cases with the highest expression of both had the largest proportion of arrhythmic cause of death (n = 27 of 36 [75%]) vs cases with low expression of both markers (n = 87 of 181 [38%]) (P = 0.006) or vs mature only (n = 10 of 14 [71%]) or active only (n = 5 of 14 [36%]). Activated fibroblast gene expression signature was enriched in arrhythmic female vs arrhythmic male cases, among other sex-specific differences in ion-channel and myosin (up-regulation of SCN4B, SCN8A, and KCNAB1 in females and KCNJ4 and MYH7B in males) expression.

CONCLUSIONS

RNA profiling of the myocardium at SCD identifies active fibrosis, undetectable by conventional clinical methods, in the presence of fixed scar and selected ion-channel dysregulation (more pronounced among female cases) as an acute vulnerable substrate for fatal arrhythmias. These findings may represent novel directions to identify patients at elevated near-term risk for SCD and critical pathways for intervention to reduce acute lethal arrhythmias.

From Cell to Gene: Deciphering the Mechanism of Heart Failure With Single-Cell Sequencing

Heart failure (HF) is a prevalent cardiovascular disease with significant morbidity and mortality rates worldwide. Due to the intricate structure of the heart, diverse cell types, and the complex pathogenesis of HF, further in-depth investigation into the underlying mechanisms  is required. The elucidation of the heterogeneity of cardiomyocytes and the intercellular communication network is particularly important. Traditional high-throughput sequencing methods provide an average measure of gene expression, failing to capture the "heterogeneity" between cells and impacting the accuracy of gene function knowledge. In contrast, single-cell sequencing techniques allow for the amplification of the entire genome or transcriptome at the individual cell level, facilitating the examination of gene structure and expression with unparalleled precision. This approach offers valuable insights into disease mechanisms, enabling the identification of changes in cellular components and gene expressions during hypertrophy associated with HF. Moreover, it reveals distinct cell populations and their unique roles in the HF microenvironment, providing a comprehensive understanding of the cellular landscape that underpins HF pathogenesis. This review focuses on the insights provided by single-cell sequencing techniques into the mechanisms underlying HF and discusses the challenges encountered in current cardiovascular research.

Identification of macrophage driver genes in fibrosis caused by different heart diseases based on omics integration

BACKGROUND

Myocardial fibrosis, a hallmark of heart disease, is closely associated with macrophages, yet the genetic pathophysiology remains incompletely understood. In this study, we utilized integrated single-cell transcriptomics and bulk RNA-seq analysis to investigate the relationship between macrophages and myocardial fibrosis across omics integration.

METHODS

We examined and curated existing single-cell data from dilated cardiomyopathy (DCM), ischemic cardiomyopathy (ICM), myocardial infarction (MI), and heart failure (HF), and analyzed the integrated data using cell communication, transcription factor identification, high dimensional weighted gene co-expression network analysis (hdWGCNA), and functional enrichment to elucidate the drivers of macrophage polarization and the macrophage-to-myofibroblast transition (MMT). Additionally, we assessed the accuracy of single-cell data from the perspective of driving factors, cell typing, anti-fibrosis performance of left ventricular assist device (LVAD). Candidate drugs were screened using L1000FWD.

RESULTS

All four heart diseases exhibit myocardial fibrosis, with only MI showing an increase in macrophage proportions. Macrophages participate in myocardial fibrosis through various fibrogenic molecules, especially evident in DCM and MI. Abnormal RNA metabolism and dysregulated transcription are significant drivers of macrophage-mediated fibrosis. Furthermore, profibrotic macrophages exhibit M1 polarization and increased MMT. In HF patients, those responding to LVAD therapy showed a significant decrease in driver gene expression, M1 polarization, and MMT. Drug repurposing identified cinobufagin as a potential therapeutic agent.

CONCLUSION

Using integrated single-cell transcriptomics, we identified the drivers of macrophage-mediated myocardial fibrosis in four heart diseases and confirmed the therapeutic effect of LVAD on improving HF with single-cell accuracy, providing novel insights into the diagnosis and treatment of myocardial fibrosis.

Identification of potential therapeutic targets for nonischemic cardiomyopathy in European ancestry: an integrated multiomics analysis

BACKGROUND

Nonischemic cardiomyopathy (NISCM) is a clinical challenge with limited therapeutic targets. This study aims to identify promising drug targets for NISCM.

METHODS

We utilized cis-pQTLs from the deCODE study, which includes data from 35,559 Icelanders, and SNPs from the FinnGen study, which includes data from 1,754 NISCM cases and 340,815 controls of Finnish ancestry. Mendelian randomization (MR) analysis was performed to estimate the causal relationship between circulating plasma protein levels and NISCM risk. Proteins with significant associations underwent false discovery rate (FDR) correction, followed by Bayesian colocalization analysis. The expression of top two proteins, LILRA5 and NELL1, was further analyzed using various NISCM datasets. Descriptions from the Human Protein Atlas (HPA) validated protein expression. The impact of environmental exposures on LILRA5 was assessed using the Comparative Toxicogenomics Database (CTD), and molecular docking identified the potential small molecule interactions.

RESULTS

MR analysis identified 255 circulating plasma proteins associated with NISCM, with 16 remaining significant after FDR correction. Bayesian colocalization analysis identified LILRA5 and NELL1 as significant, with PP.H4 > 0.8. LILRA5 has a protective effect (OR = 0.758, 95% CI, 0.670-0.857) while NELL1 displays the risk effect (OR = 1.290, 95% CI, 1.199-1.387) in NISCM. Decreased LILRA5 expression was found in NISCM such as diabetic, hypertrophic, dilated, and inflammatory cardiomyopathy, while NELL1 expression increased in hypertrophic cardiomyopathy. HPA data indicated high LILRA5 expression in neutrophils, macrophages and endothelial cells within normal heart and limited NELL1 expression. Immune infiltration analysis revealed decreased neutrophil in diabetic cardiomyopathy. CTD analysis identified several small molecules that affect LILRA5 mRNA expression. Among these, Estradiol, Estradiol-3-benzoate, Gadodiamide, Topotecan, and Testosterone were found to stably bind to the LILRA5 protein at the conserved VAL-15 or THR-133 residues in the Ig-like C2 domain.

CONCLUSION

Based on European Ancestry Cohort, this study reveals that LILRA5 and NELL1 are potential therapeutic targets for NISCM, with LILRA5 showing particularly promising prospects in diabetic cardiomyopathy. Several small molecules interact with LILRA5, implying potential clinical implication.

Hyperactivation of ATF4/TGF-β1 signaling contributes to the progressive cardiac fibrosis in Arrhythmogenic cardiomyopathy caused by DSG2 Variant

BACKGROUND

Arrhythmogenic cardiomyopathy (ACM) is an inherited cardiomyopathy characterized with progressive cardiac fibrosis and heart failure. However, the exact mechanism driving the progression of cardiac fibrosis and heart failure in ACM remains elusive. This study aims to investigate the underlying mechanisms of progressive cardiac fibrosis in ACM caused by newly identified Desmoglein-2 (DSG2) variation.

METHODS

We identified homozygous DSG2F531C variant in a family with 8 ACM patients using whole-exome sequencing and generated Dsg2F536C knock-in mice. Neonatal and adult mouse ventricular myocytes isolated from Dsg2F536C knock-in mice were used. We performed functional, transcriptomic and mass spectrometry analyses to evaluate the mechanisms of ACM caused by DSG2F531C variant.

RESULTS

All eight patients with ACM were homozygous for DSG2F531C variant. Dsg2F536C/F536C mice displayed cardiac enlargement, dysfunction, and progressive cardiac fibrosis in both ventricles. Mechanistic investigations revealed that the variant DSG2-F536C protein underwent misfolding, leading to its recognition by BiP within the endoplasmic reticulum, which triggered endoplasmic reticulum stress, activated the PERK-ATF4 signaling pathway and increased ATF4 levels in cardiomyocytes. Increased ATF4 facilitated the expression of TGF-β1 in cardiomyocytes, thereby activating cardiac fibroblasts through paracrine signaling and ultimately promoting cardiac fibrosis in Dsg2F536C/F536C mice. Notably, inhibition of the PERK-ATF4 signaling attenuated progressive cardiac fibrosis and cardiac systolic dysfunction in Dsg2F536C/F536C mice.

CONCLUSIONS

Hyperactivation of the ATF4/TGF-β1 signaling in cardiomyocytes emerges as a novel mechanism underlying progressive cardiac fibrosis in ACM. Targeting the ATF4/TGF-β1 signaling may be a novel therapeutic target for managing ACM.

LIANA+ provides an all-in-one framework for cell-cell communication inference

The growing availability of single-cell and spatially resolved transcriptomics has led to the development of many approaches to infer cell-cell communication, each capturing only a partial view of the complex landscape of intercellular signalling. Here we present LIANA+, a scalable framework built around a rich knowledge base to decode coordinated inter- and intracellular signalling events from single- and multi-condition datasets in both single-cell and spatially resolved data. By extending and unifying established methodologies, LIANA+ provides a comprehensive set of synergistic components to study cell-cell communication via diverse molecular mediators, including those measured in multi-omics data. LIANA+ is accessible at https://github.com/saezlab/liana-py with extensive vignettes ( https://liana-py.readthedocs.io/ ) and provides an all-in-one solution to intercellular communication inference.

Pathophysiological insights into HFpEF from studies of human cardiac tissue

Heart failure with preserved ejection fraction (HFpEF) is a major, worldwide health-care problem. Few therapies for HFpEF exist because the pathophysiology of this condition is poorly defined and, increasingly, postulated to be diverse. Although perturbations in other organs contribute to the clinical profile in HFpEF, altered cardiac structure, function or both are the primary causes of this heart failure syndrome. Therefore, studying myocardial tissue is fundamental to improve pathophysiological insights and therapeutic discovery in HFpEF. Most studies of myocardial changes in HFpEF have relied on cardiac tissue from animal models without (or with limited) confirmatory studies in human cardiac tissue. Animal models of HFpEF have evolved based on theoretical HFpEF aetiologies, but these models might not reflect the complex pathophysiology of human HFpEF. The focus of this Review is the pathophysiological insights gained from studies of human HFpEF myocardium. We outline the rationale for these studies, the challenges and opportunities in obtaining myocardial tissue from patients with HFpEF and relevant comparator groups, the analytical approaches, the pathophysiological insights gained to date and the remaining knowledge gaps. Our objective is to provide a roadmap for future studies of cardiac tissue from diverse cohorts of patients with HFpEF, coupling discovery biology with measures to account for pathophysiological diversity.

Single-cell analysis of chromatin and expression reveals age- and sex-associated alterations in the human heart

Sex differences and age-related changes in the human heart at the tissue, cell, and molecular level have been well-documented and many may be relevant for cardiovascular disease. However, how molecular programs within individual cell types vary across individuals by age and sex remains poorly characterized. To better understand this variation, we performed single-nucleus combinatorial indexing (sci) ATAC- and RNA-Seq in human heart samples from nine donors. We identify hundreds of differentially expressed genes by age and sex and find epigenetic signatures of variation in ATAC-Seq data in this discovery cohort. We then scale up our single-cell RNA-Seq analysis by combining our data with five recently published single nucleus RNA-Seq datasets of healthy adult hearts. We find variation such as metabolic alterations by sex and immune changes by age in differential expression tests, as well as alterations in abundance of cardiomyocytes by sex and neurons with age. In addition, we compare our adult-derived ATAC-Seq profiles to analogous fetal cell types to identify putative developmental-stage-specific regulatory factors. Finally, we train predictive models of cell-type-specific RNA expression levels utilizing ATAC-Seq profiles to link distal regulatory sequences to promoters, quantifying the predictive value of a simple TF-to-expression regulatory grammar and identifying cell-type-specific TFs. Our analysis represents the largest single-cell analysis of cardiac variation by age and sex to date and provides a resource for further study of healthy cardiac variation and transcriptional regulation at single-cell resolution.

Assessing the impact of transcriptomics data analysis pipelines on downstream functional enrichment results

Transcriptomics is widely used to assess the state of biological systems. There are many tools for the different steps, such as normalization, differential expression, and enrichment. While numerous studies have examined the impact of method choices on differential expression results, little attention has been paid to their effects on further downstream functional analysis, which typically provides the basis for interpretation and follow-up experiments. To address this, we introduce FLOP, a comprehensive nextflow-based workflow combining methods to perform end-to-end analyses of transcriptomics data. We illustrate FLOP on datasets ranging from end-stage heart failure patients to cancer cell lines. We discovered effects not noticeable at the gene-level, and observed that not filtering the data had the highest impact on the correlation between pipelines in the gene set space. Moreover, we performed three benchmarks to evaluate the 12 pipelines included in FLOP, and confirmed that filtering is essential in scenarios of expected moderate-to-low biological signal. Overall, our results underscore the impact of carefully evaluating the consequences of the choice of preprocessing methods on downstream enrichment analyses. We envision FLOP as a valuable tool to measure the robustness of functional analyses, ultimately leading to more reliable and conclusive biological findings.

Unraveling cardiomyocyte responses and intercellular communication alterations in primary carnitine deficiency cardiomyopathy via single-nucleus RNA sequencing

Background

Primary Carnitine Deficiency (PCD) is a potentially life-threatening autosomal recessive monogenic disorder arising from mutations in the organic cation transporter 2 (OCTN2) gene. Dilated cardiomyopathy (DCM) is a prevalent symptom associated with this condition, and episodes of metabolic disturbance may lead to sudden death. However, the pathogenic mechanism remains unclear. Here, we sought to investigate the response of cardiomyocytes and alterations in the intercellular communication in individuals with PCD DCM.

Methods

The GSE211650 dataset was downloaded. Subsequently, modular analysis was performed using hdWGCNA. SCENIC was employed for transcription factor analysis. Monocle2 and SCP were applied to conduct trajectory inference and characterize dynamic features. CellChat was used to investigate intercellular interactions.

Results

OCTN2-deficient cardiomyocytes displayed transcriptomic alterations indicative of reduced contractility, developmental abnormalities, and fibrosis. The reduced expression of genes encoding troponin, myosin, and calcium ion transporters may underlie the observed decrease in contractility. Suppressed Wnt signaling and downregulated transcription factors associated with myocardial development suggest potential developmental disturbances in cardiomyocytes. Growth arrest-specific 6 (GAS6) secreted by TNNC1 high cardiomyocytes is implicated in myocardial inflammation and fibrosis. Macrophages-derived secreted phosphoprotein 1 (SPP1) promotes the activation of fibroblasts. Furthermore, there was a reduction in neuronal genes in the OCTN2-deficient group.

Conclusions

Our research has unveiled, for the first time, the responses of cardiomyocytes and alterations in the intercellular communication in PCD DCM, offering valuable insights for the precision treatment of this condition.

Altered Inflammatory State and Mitochondrial Function Identified by Transcriptomics in Paediatric Congenital Heart Patients Prior to Surgical Repair

Congenital heart disease (CHD) remains the most common birth defect, with surgical intervention required in complex cases. Right ventricle (RV) function is known to be a major predictor of sustained cardiac health in these patients; thus, by elucidating the divergent profiles between CHD and the control through tissue analysis, this study aims to identify new avenues of investigation into the mechanisms surrounding reduced RV function. Transcriptomic profiling, in-silico deconvolution and functional network analysis were conducted on RV biopsies, identifying an increase in the mitochondrial dysfunction genes RPPH1 and RMPR (padj = 4.67 × 10-132, 2.23 × 10-107), the cytotoxic T-cell markers CD8a, LAGE3 and CD49a (p = 0.0006, p < 0.0001, and p = 0.0118) and proinflammatory caspase-1 (p = 0.0055) in CHD. Gene-set enrichment identified mitochondrial dysfunctional pathways, predominately changes within oxidative phosphorylation processes. The negative regulation of mitochondrial functions and metabolism was identified in the network analysis, with dysregulation of the mitochondrial complex formation. A histological analysis confirmed an increase in cellular bodies in the CHD RV tissue and positive staining for both CD45 and CD8, which was absent in the control. The deconvolution of bulk RNAseq data suggests a reduction in CD4+ T cells (p = 0.0067) and an increase in CD8+ T cells (p = 0.0223). The network analysis identified positive regulation of the immune system and cytokine signalling clusters in the inflammation functional network, as there were lymphocyte activation and leukocyte differentiation. Utilising RV tissue from paediatric patients undergoing CHD cardiac surgery, this study identifies dysfunctional mitochondrial pathways and an increase in inflammatory T-cell presence prior to reparative surgery.

Deciphering the role of SMARCA4 in cardiac disorders: Insights from single-cell studies on dilated cardiomyopathy and coronary heart disease

BACKGROUND

Dilated cardiomyopathy (DCM) and coronary heart disease (CHD) stand as two of the foremost causes of mortality. However, the comprehensive comprehension of the regulatory mechanisms governing DCM and CHD remains limited, particularly from the vantage point of single-cell transcriptional analysis.

METHOD

We used the GSE121893 dataset from the GEO database, analyzing single-cell expressions with tools like DropletUtils, Seurat, and Monocle. We also utilized the GSVA package for comparing gene roles in DCM and CHD, Finally, we conducted qRT-PCR and Western blot analyses to measure the expression levels of SMARCA4, Col1A1, Col3A1 and α-SMA, and the role of SMARCA4 on fibroblasts were explored by EdU and Transwell assay.

RESULTS

Our analysis identified six cell types in heart tissue, with fibroblasts showing the most interaction with other cells. DEGs in fibroblasts were linked to muscle development and morphogenesis. Pseudotime analysis revealed the dynamics of fibroblast changes in both the normal and disease groups and many transcription factors (TFs) potentially involved in this process. Among these TFs, SMARCA4 which was translated into protein BRG1, showed the most significantly difference. In vivo experiments have demonstrated that SMARCA4 indeed promoted fibroblasts proliferation and migration.

CONCLUSION

This study provides a clearer understanding of cell-type dynamics in heart diseases, emphasizing the role of fibroblasts and the significance of SMARCA4 in their function. Our results offer insights into the cellular mechanisms underlying DCM and CHD, potentially guiding future therapeutic strategies.

Building and analyzing metacells in single-cell genomics data

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

LMNA -Related Dilated Cardiomyopathy: Single-Cell Transcriptomics during Patient-derived iPSC Differentiation Support Cell type and Lineage-specific Dysregulation of Gene Expression and Development for Cardiomyocytes and Epicardium-Derived Cells with Lamin A/C Haploinsufficiency

LMNA -Related Dilated Cardiomyopathy (DCM) is an autosomal-dominant genetic condition with cardiomyocyte and conduction system dysfunction often resulting in heart failure or sudden death. The condition is caused by mutation in the Lamin A/C ( LMNA ) gene encoding Type-A nuclear lamin proteins involved in nuclear integrity, epigenetic regulation of gene expression, and differentiation. Molecular mechanisms of disease are not completely understood, and there are no definitive treatments to reverse progression or prevent mortality. We investigated possible mechanisms of LMNA -Related DCM using induced pluripotent stem cells derived from a family with a heterozygous LMNA c.357-2A>G splice-site mutation. We differentiated one LMNA mutant iPSC line derived from an affected female (Patient) and two non-mutant iPSC lines derived from her unaffected sister (Control) and conducted single-cell RNA sequencing for 12 samples (4 Patient and 8 Control) across seven time points: Day 0, 2, 4, 9, 16, 19, and 30. Our bioinformatics workflow identified 125,554 cells in raw data and 110,521 (88%) high-quality cells in sequentially processed data. Unsupervised clustering, cell annotation, and trajectory inference found complex heterogeneity: ten main cell types; many possible subtypes; and lineage bifurcation for Cardiac Progenitors to Cardiomyocytes (CM) and Epicardium-Derived Cells (EPDC). Data integration and comparative analyses of Patient and Control cells found cell type and lineage differentially expressed genes (DEG) with enrichment to support pathway dysregulation. Top DEG and enriched pathways included: 10 ZNF genes and RNA polymerase II transcription in Pluripotent cells (PP); BMP4 and TGF Beta/BMP signaling, sarcomere gene subsets and cardiogenesis, CDH2 and EMT in CM; LMNA and epigenetic regulation and DDIT4 and mTORC1 signaling in EPDC. Top DEG also included: XIST and other X-linked genes, six imprinted genes: SNRPN , PWAR6 , NDN , PEG10 , MEG3 , MEG8 , and enriched gene sets in metabolism, proliferation, and homeostasis. We confirmed Lamin A/C haploinsufficiency by allelic expression and Western blot. Our complex Patient-derived iPSC model for Lamin A/C haploinsufficiency in PP, CM, and EPDC provided support for dysregulation of genes and pathways, many previously associated with Lamin A/C defects, such as epigenetic gene expression, signaling, and differentiation. Our findings support disruption of epigenomic developmental programs as proposed in other LMNA disease models. We recognized other factors influencing epigenetics and differentiation; thus, our approach needs improvement to further investigate this mechanism in an iPSC-derived model.

Multi-level transcriptomic analysis of LMNA -related dilated cardiomyopathy identifies disease-driving processes

LMNA- related dilated cardiomyopathy ( LMNA -DCM) is one of the most severe forms of DCM. The incomplete understanding of the molecular disease mechanisms results in lacking treatment options, leading to high mortality amongst patients. Here, using an inducible, cardiomyocyte-specific lamin A/C depletion mouse model, we conducted a comprehensive transcriptomic study, combining both bulk and single nucleus RNA sequencing, and spanning LMNA -DCM disease progression, to identify potential disease drivers. Our refined analysis pipeline identified 496 genes already misregulated early in disease. The expression of these genes was largely driven by disease specific cardiomyocyte sub-populations and involved biological processes mediating cellular response to DNA damage, cytosolic pattern recognition, and innate immunity. Indeed, DNA damage in LMNA -DCM hearts was significantly increased early in disease and correlated with reduced cardiomyocyte lamin A levels. Activation of cytosolic pattern recognition in cardiomyocytes was independent of cGAS, which is rarely expressed in cardiomyocytes, but likely occurred downstream of other pattern recognition sensors such as IFI16. Altered gene expression in cardiac fibroblasts and immune cell infiltration further contributed to tissue-wide changes in gene expression. Our transcriptomic analysis further predicted significant alterations in cell-cell communication between cardiomyocytes, fibroblasts, and immune cells, mediated through early changes in the extracellular matrix (ECM) in the LMNA -DCM hearts. Taken together, our work suggests a model in which nuclear damage in cardiomyocytes leads to activation of DNA damage responses, cytosolic pattern recognition pathway, and other signaling pathways that activate inflammation, immune cell recruitment, and transcriptional changes in cardiac fibroblasts, which collectively drive LMNA -DCM pathogenesis.

Revisiting Cardiac Biology in the Era of Single Cell and Spatial Omics

Throughout our lifetime, each beat of the heart requires the coordinated action of multiple cardiac cell types. Understanding cardiac cell biology, its intricate microenvironments, and the mechanisms that govern their function in health and disease are crucial to designing novel therapeutical and behavioral interventions. Recent advances in single-cell and spatial omics technologies have significantly propelled this understanding, offering novel insights into the cellular diversity and function and the complex interactions of cardiac tissue. This review provides a comprehensive overview of the cellular landscape of the heart, bridging the gap between suspension-based and emerging in situ approaches, focusing on the experimental and computational challenges, comparative analyses of mouse and human cardiac systems, and the rising contextualization of cardiac cells within their niches. As we explore the heart at this unprecedented resolution, integrating insights from both mouse and human studies will pave the way for novel diagnostic tools and therapeutic interventions, ultimately improving outcomes for patients with cardiovascular diseases.

Patterns of Unwanted Biological and Technical Expression Variation Among 49 Human Tissues

Tissue gene expression studies are impacted by biological and technical sources of variation, which can be broadly classified into wanted and unwanted variation. The latter, if not addressed, results in misleading biological conclusions. Methods have been proposed to reduce unwanted variation, such as normalization and batch correction. A more accurate understanding of all causes of variation could significantly improve the ability of these methods to remove unwanted variation while retaining variation corresponding to the biological question of interest. We used 17,282 samples from 49 human tissues in the Genotype-Tissue Expression data set (v8) to investigate patterns and causes of expression variation. Transcript expression was transformed to z-scores, and only the most variable 2% of transcripts were evaluated and clustered based on coexpression patterns. Clustered gene sets were assigned to different biological or technical causes based on histologic appearances and metadata elements. We identified 522 variable transcript clusters (median: 11 per tissue) among the samples. Of these, 63% were confidently explained, 16% were likely explained, 7% were low confidence explanations, and 14% had no clear cause. Histologic analysis annotated 46 clusters. Other common causes of variability included sex, sequencing contamination, immunoglobulin diversity, and compositional tissue differences. Less common biological causes included death interval (Hardy score), disease status, and age. Technical causes included blood draw timing and harvesting differences. Many of the causes of variation in bulk tissue expression were identifiable in the Tabula Sapiens data set of single-cell expression. This is among the largest explorations of the underlying sources of tissue expression variation. It uncovered expected and unexpected causes of variable gene expression and demonstrated the utility of matched histologic specimens. It further demonstrated the value of acquiring meaningful tissue harvesting metadata elements to use for improved normalization, batch correction, and analysis of both bulk and single-cell RNA-seq data.

Identification of Potential lncRNA-miRNA-mRNA Regulatory Network Contributing to Arrhythmogenic Right Ventricular Cardiomyopathy

Arrhythmogenic right ventricular cardiomyopathy (ARVC) can lead to sudden cardiac death and life-threatening heart failure. Due to its high fatality rate and limited therapies, the pathogenesis and diagnosis biomarker of ARVC needs to be explored urgently. This study aimed to explore the lncRNA-miRNA-mRNA competitive endogenous RNA (ceRNA) network in ARVC. The mRNA and lncRNA expression datasets obtained from the Gene Expression Omnibus (GEO) database were used to analyze differentially expressed mRNA (DEM) and lncRNA (DElnc) between ARVC and non-failing controls. Differentially expressed miRNAs (DEmiRs) were obtained from the previous profiling work. Using starBase to predict targets of DEmiRs and intersecting with DEM and DElnc, a ceRNA network of lncRNA-miRNA-mRNA was constructed. The DEM and DElnc were validated by real-time quantitative PCR in human heart tissue. Protein-protein interaction network and weighted gene co-expression network analyses were used to identify hub genes. A logistic regression model for ARVC diagnostic prediction was established with the hub genes and their ceRNA pairs in the network. A total of 448 DEMs (282 upregulated and 166 downregulated) were identified, mainly enriched in extracellular matrix and fibrosis-related GO terms and KEGG pathways, such as extracellular matrix organization and collagen fibril organization. Four mRNAs and two lncRNAs, including COL1A1, COL5A1, FBN1, BGN, XIST, and LINC00173 identified through the ceRNA network, were validated by real-time quantitative PCR in human heart tissue and used to construct a logistic regression model. Good ARVC diagnostic prediction performance for the model was shown in both the training set and the validation set. The potential lncRNA-miRNA-mRNA regulatory network and logistic regression model established in our study may provide promising diagnostic methods for ARVC.

Mast cells: a novel therapeutic avenue for cardiovascular diseases?

Mast cells are tissue-resident immune cells strategically located in different compartments of the normal human heart (the myocardium, pericardium, aortic valve, and close to nerves) as well as in atherosclerotic plaques. Cardiac mast cells produce a broad spectrum of vasoactive and proinflammatory mediators, which have potential roles in inflammation, angiogenesis, lymphangiogenesis, tissue remodelling, and fibrosis. Mast cells release preformed mediators (e.g. histamine, tryptase, and chymase) and de novo synthesized mediators (e.g. cysteinyl leukotriene C4 and prostaglandin D2), as well as cytokines and chemokines, which can activate different resident immune cells (e.g. macrophages) and structural cells (e.g. fibroblasts and endothelial cells) in the human heart and aorta. The transcriptional profiles of various mast cell populations highlight their potential heterogeneity and distinct gene and proteome expression. Mast cell plasticity and heterogeneity enable these cells the potential for performing different, even opposite, functions in response to changing tissue contexts. Human cardiac mast cells display significant differences compared with mast cells isolated from other organs. These characteristics make cardiac mast cells intriguing, given their dichotomous potential roles of inducing or protecting against cardiovascular diseases. Identification of cardiac mast cell subpopulations represents a prerequisite for understanding their potential multifaceted roles in health and disease. Several new drugs specifically targeting human mast cell activation are under development or in clinical trials. Mast cells and/or their subpopulations can potentially represent novel therapeutic targets for cardiovascular disorders.

The Macrophage Landscape Across the Lifespan of a Human Cardiac Allograft

BACKGROUND

Much of our knowledge of organ rejection after transplantation is derived from rodent models.

METHODS

We used single-nucleus RNA sequencing to investigate the inflammatory myocardial microenvironment in human pediatric cardiac allografts at different stages after transplantation. We distinguished donor- from recipient-derived cells using naturally occurring genetic variants embedded in single-nucleus RNA sequencing data.

RESULTS

Donor-derived tissue resident macrophages, which accompany the allograft into the recipient, are lost over time after transplantation. In contrast, monocyte-derived macrophages from the recipient populate the heart within days after transplantation and form 2 macrophage populations: recipient MP1 and recipient MP2. Recipient MP2s have cell signatures similar to donor-derived resident macrophages; however, they lack signatures of pro-reparative phagocytic activity typical of donor-derived resident macrophages and instead express profibrotic genes. In contrast, recipient MP1s express genes consistent with hallmarks of cellular rejection. Our data suggest that recipient MP1s activate a subset of natural killer cells, turning them into a cytotoxic cell population through feed-forward signaling between recipient MP1s and natural killer cells.

CONCLUSIONS

Our findings reveal an imbalance of donor-derived and recipient-derived macrophages in the pediatric cardiac allograft that contributes to allograft failure.

Fast skeletal myosin binding protein-C expression exacerbates dysfunction in heart failure

During heart failure, gene and protein expression profiles undergo extensive compensatory and pathological remodeling. We previously observed that fast skeletal myosin binding protein-C (fMyBP-C) is upregulated in diseased mouse hearts. While fMyBP-C shares significant homology with its cardiac paralog, cardiac myosin binding protein-C (cMyBP-C), there are key differences that may affect cardiac function. However, it is unknown if the expression of fMyBP-C expression in the heart is a pathological or compensatory response. We aim to elucidate the cardiac consequence of either increased or knockout of fMyBP-C expression. To determine the sufficiency of fMyBP-C to cause cardiac dysfunction, we generated cardiac-specific fMyBP-C over-expression mice. These mice were further crossed into a cMyBP-C null model to assess the effect of fMyBP-C in the heart in the complete absence of cMyBP-C. Finally, fMyBP-C null mice underwent transverse aortic constriction (TAC) to define the requirement of fMyBP-C during heart failure development. We confirmed the upregulation of fMyBP-C in several models of cardiac disease, including the use of lineage tracing. Low levels of fMyBP-C caused mild cardiac remodeling and sarcomere dysfunction. Exclusive expression of fMyBP-C in a heart failure model further exacerbated cardiac pathology. Following 8 weeks of TAC, fMyBP-C null mice demonstrated greater protection against heart failure development. Mechanistically, this may be due to the differential regulation of the myosin super-relaxed state. These findings suggest that the elevated expression of fMyBP-C in diseased hearts is a pathological response. Targeted therapies to prevent upregulation of fMyBP-C may prove beneficial in the treatment of heart failure.

Significance Statement

Recently, the sarcomere - the machinery that controls heart and muscle contraction - has emerged as a central target for development of cardiac therapeutics. However, there remains much to understand about how the sarcomere is modified in response to disease. We recently discovered that a protein normally expressed in skeletal muscle, is present in the heart in certain settings of heart disease. How this skeletal muscle protein affects the function of the heart remained unknown. Using genetically engineered mouse models to modulate expression of this skeletal muscle protein, we determined that expression of this skeletal muscle protein in the heart negatively affects cardiac performance. Importantly, deletion of this protein from the heart could improve heart function suggesting a possible therapeutic avenue.

PBX/Knotted 1 homeobox-2 (PKNOX2) is a novel regulator of myocardial fibrosis

Much effort has been made to uncover the cellular heterogeneities of human hearts by single-nucleus RNA sequencing. However, the cardiac transcriptional regulation networks have not been systematically described because of the limitations in detecting transcription factors. In this study, we optimized a pipeline for isolating nuclei and conducting single-nucleus RNA sequencing targeted to detect a higher number of cell signal genes and an optimal number of transcription factors. With this unbiased protocol, we characterized the cellular composition of healthy human hearts and investigated the transcriptional regulation networks involved in determining the cellular identities and functions of the main cardiac cell subtypes. Particularly in fibroblasts, a novel regulator, PKNOX2, was identified as being associated with physiological fibroblast activation in healthy hearts. To validate the roles of these transcription factors in maintaining homeostasis, we used single-nucleus RNA-sequencing analysis of transplanted failing hearts focusing on fibroblast remodelling. The trajectory analysis suggested that PKNOX2 was abnormally decreased from fibroblast activation to pathological myofibroblast formation. Both gain- and loss-of-function in vitro experiments demonstrated the inhibitory role of PKNOX2 in pathological fibrosis remodelling. Moreover, fibroblast-specific overexpression and knockout of PKNOX2 in a heart failure mouse model induced by transverse aortic constriction surgery significantly improved and aggravated myocardial fibrosis, respectively. In summary, this study established a high-quality pipeline for single-nucleus RNA-sequencing analysis of heart muscle. With this optimized protocol, we described the transcriptional regulation networks of the main cardiac cell subtypes and identified PKNOX2 as a novel regulator in suppressing fibrosis and a potential therapeutic target for future translational studies.

Histone H1.0 couples cellular mechanical behaviors to chromatin structure

Tuning of genome structure and function is accomplished by chromatin-binding proteins, which determine the transcriptome and phenotype of the cell. Here we investigate how communication between extracellular stress and chromatin structure may regulate cellular mechanical behaviors. We demonstrate that histone H1.0, which compacts nucleosomes into higher-order chromatin fibers, controls genome organization and cellular stress response. We show that histone H1.0 has privileged expression in fibroblasts across tissue types and that its expression is necessary and sufficient to induce myofibroblast activation. Depletion of histone H1.0 prevents cytokine-induced fibroblast contraction, proliferation and migration via inhibition of a transcriptome comprising extracellular matrix, cytoskeletal and contractile genes, through a process that involves locus-specific H3K27 acetylation. Transient depletion of histone H1.0 in vivo prevents fibrosis in cardiac muscle. These findings identify an unexpected role of linker histones to orchestrate cellular mechanical behaviors, directly coupling force generation, nuclear organization and gene transcription.

A Multimodal Omics Framework to Empower Target Discovery for Cardiovascular Regeneration

Ischaemic heart disease is a global healthcare challenge with high morbidity and mortality. Early revascularisation in acute myocardial infarction has improved survival; however, limited regenerative capacity and microvascular dysfunction often lead to impaired function and the development of heart failure. New mechanistic insights are required to identify robust targets for the development of novel strategies to promote regeneration. Single-cell RNA sequencing (scRNA-seq) has enabled profiling and analysis of the transcriptomes of individual cells at high resolution. Applications of scRNA-seq have generated single-cell atlases for multiple species, revealed distinct cellular compositions for different regions of the heart, and defined multiple mechanisms involved in myocardial injury-induced regeneration. In this review, we summarise findings from studies of healthy and injured hearts in multiple species and spanning different developmental stages. Based on this transformative technology, we propose a multi-species, multi-omics, meta-analysis framework to drive the discovery of new targets to promote cardiovascular regeneration.

AAV9:PKP2 improves heart function and survival in a Pkp2-deficient mouse model of arrhythmogenic right ventricular cardiomyopathy

BACKGROUND

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a familial cardiac disease associated with ventricular arrhythmias and an increased risk of sudden cardiac death. Currently, there are no approved treatments that address the underlying genetic cause of this disease, representing a significant unmet need. Mutations in Plakophilin-2 (PKP2), encoding a desmosomal protein, account for approximately 40% of ARVC cases and result in reduced gene expression.

METHODS

Our goal is to examine the feasibility and the efficacy of adeno-associated virus 9 (AAV9)-mediated restoration of PKP2 expression in a cardiac specific knock-out mouse model of Pkp2.

RESULTS

We show that a single dose of AAV9:PKP2 gene delivery prevents disease development before the onset of cardiomyopathy and attenuates disease progression after overt cardiomyopathy. Restoration of PKP2 expression leads to a significant extension of lifespan by restoring cellular structures of desmosomes and gap junctions, preventing or halting decline in left ventricular ejection fraction, preventing or reversing dilation of the right ventricle, ameliorating ventricular arrhythmia event frequency and severity, and preventing adverse fibrotic remodeling. RNA sequencing analyses show that restoration of PKP2 expression leads to highly coordinated and durable correction of PKP2-associated transcriptional networks beyond desmosomes, revealing a broad spectrum of biological perturbances behind ARVC disease etiology.

CONCLUSIONS

We identify fundamental mechanisms of PKP2-associated ARVC beyond disruption of desmosome function. The observed PKP2 dose-function relationship indicates that cardiac-selective AAV9:PKP2 gene therapy may be a promising therapeutic approach to treat ARVC patients with PKP2 mutations.

Microfluidic platform for omics analysis on single cells with diverse morphology and size: A review

BACKGROUND

Microfluidic techniques have emerged as powerful tools in single-cell research, facilitating the exploration of omics information from individual cells. Cell morphology is crucial for gene expression and physiological processes. However, there is currently a lack of integrated analysis of morphology and single-cell omics information. A critical challenge remains: what platform technologies are the best option to decode omics data of cells that are complex in morphology and size?

RESULTS

This review highlights achievements in microfluidic-based single-cell omics and isolation of cells based on morphology, along with other cell sorting methods based on physical characteristics. Various microfluidic platforms for single-cell isolation are systematically presented, showcasing their diversity and adaptability. The discussion focuses on microfluidic devices tailored to the distinct single-cell isolation requirements in plants and animals, emphasizing the significance of considering cell morphology and cell size in optimizing single-cell omics strategies. Simultaneously, it explores the application of microfluidic single-cell sorting technologies to single-cell sequencing, aiming to effectively integrate information about cell shape and size.

SIGNIFICANCE AND NOVELTY

The novelty lies in presenting a comprehensive overview of recent accomplishments in microfluidic-based single-cell omics, emphasizing the integration of different microfluidic platforms and their implications for cell morphology-based isolation. By underscoring the pivotal role of the specialized morphology of different cells in single-cell research, this review provides robust support for delving deeper into the exploration of single-cell omics data.

Brooklyn plots to identify co-expression dysregulation in single cell sequencing

Altered open chromatin regions, impacting gene expression, is a feature of some human disorders. We discovered it is possible to detect global changes in genomically-related adjacent gene co-expression within single cell RNA sequencing (scRNA-seq) data. We built a software package to generate and test non-randomness using 'Brooklyn plots' to identify the percent of genes significantly co-expressed from the same chromosome in ∼10 MB intervals across the genome. These plots establish an expected low baseline of co-expression in scRNA-seq from most cell types, but, as seen in dilated cardiomyopathy cardiomyocytes, altered patterns of open chromatin appear. These may relate to larger regions of transcriptional bursting, observable in single cell, but not bulk datasets.

Comparison of the negative effect of remimazolam and propofol on cardiac contractility: Analysis of a randomised parallel-group trial and a preclinical ex vivo study

Remimazolam is a newly developed ultra-short-acting benzodiazepine that exerts sedative effects. This study aimed to clarify the effects of remimazolam on cardiac contractility. In a randomised-parallel group trial, haemodynamic parameters were compared between propofol (n = 11) and remimazolam (n = 12) groups during the induction of general anaesthesia in patients undergoing non-cardiac surgery. In a preclinical study, the direct effects of remimazolam on cardiac contractility were also evaluated using isolated rat hearts. RNA sequence data obtained from rat and human hearts were analysed to assess the expression patterns of the cardiac γ-aminobutyric acid type A (GABAA ) receptor subunits. In a clinical study, the proportional change of the maximum rate of arterial pressure rise was milder during the study period in the remimazolam group (propofol: -52.6 [10.2] (mean [standard deviation])% vs. remimazolam: -39.7% [10.5%], p = 0.007). In a preclinical study, remimazolam did not exert a negative effect on left ventricle developed pressure, whereas propofol did exert a negative effect after bolus administration of a high dose (propofol: -26.9% [3.5%] vs. remimazolam: -1.1 [6.9%], p < 0.001). Analysis of the RNA sequence revealed a lack of γ subunits, which are part of the major benzodiazepine binding site of the GABAA receptor, in rat and human hearts. These results indicate that remimazolam does not have a direct negative effect on cardiac contractility, which might contribute to its milder effect on cardiac contractility during the induction of general anaesthesia. The expression patterns of cardiac GABAA receptor subunits might be associated with the unique pharmacokinetics of benzodiazepines in the heart.

Bone morphogenic protein-4 availability in the cardiac microenvironment controls inflammation and fibrosis in autoimmune myocarditis

Myocarditis is an inflammatory heart disease that leads to loss of cardiomyocytes and frequently precipitates fibrotic remodeling of the myocardium, culminating in heart failure. However, the molecular mechanisms underlying immune cell control and maintenance of tissue integrity in the inflamed cardiac microenvironment remain elusive. In this study, we found that bone morphogenic protein-4 (BMP4) gradients maintain cardiac tissue homeostasis by single-cell transcriptomics analyses of inflamed murine and human myocardial tissues. Cardiac BMP pathway dysregulation was reflected by reduced BMP4 serum concentration in patients with myocarditis. Restoration of BMP signaling by antibody-mediated neutralization of the BMP inhibitors gremlin-1 and gremlin-2 ameliorated T cell-induced myocardial inflammation in mice. Moreover, progression to inflammatory cardiomyopathy was blocked through the reduction of fibrotic remodeling and preservation of cardiomyocyte integrity. These results unveil the BMP4-gremlin axis as a druggable pathway for the treatment of myocardial inflammation, limiting the severe sequelae of cardiac fibrosis and heart failure.

Enrichment on steps, not genes, improves inference of differentially expressed pathways

Enrichment analysis is frequently used in combination with differential expression data to investigate potential commonalities amongst lists of genes and generate hypotheses for further experiments. However, current enrichment analysis approaches on pathways ignore the functional relationships between genes in a pathway, particularly OR logic that occurs when a set of proteins can each individually perform the same step in a pathway. As a result, these approaches miss pathways with large or multiple sets because of an inflation of pathway size (when measured as the total gene count) relative to the number of steps. We address this problem by enriching on step-enabling entities in pathways. We treat sets of protein-coding genes as single entities, and we also weight sets to account for the number of genes in them using the multivariate Fisher's noncentral hypergeometric distribution. We then show three examples of pathways that are recovered with this method and find that the results have significant proportions of pathways not found in gene list enrichment analysis.

Detection of gene mutation in the prognosis of a patient with arrhythmogenic right ventricular cardiomyopathy: a case report

BACKGROUND

Arrhythmogenic right ventricular cardiomyopathy (ARVC), or more recently known as arrhythmogenic cardiomyopathy (ACM), is an heritable disorder of the myocardium characterized by progressive fibrofatty replacement the heart muscle and risk of ventricular arrhythmias and sudden cardiac death (SCD). We report a case study to demonstrate the role of gene mutation detection in risk stratification for primary prevention of SCD in a young patient diagnosed with ARVC.

CASE PRESENTATION

A 15-year-old Asian (Vietnamese) male patient with no history of documented tachyarrhythmia or syncope and a family history of potential SCD was admitted due to palpitations. Clinical findings and work-up including cardiac magnetic resonance imaging (MRI) were highly suggestive of ARVC. Gene sequencing was performed for SCD risk stratification, during which PKP2 gene mutation was found. Based on the individualized risk stratification, an ICD was implanted for primary prevention of SCD. At 6 months post ICD implantation, the device detected and successfully delivered an appropriate shock to terminate an episode of potentially fatal ventricular arrhythmia. ICD implantation was therefore proven to be appropriate in this patient.

CONCLUSIONS

While gene mutations are known to be an important factor in the diagnosis of ARVC according to the 2010 Task Force Criteria and recent clinical guidelines, their role in risk stratification of SCD remains controversial. Our case demonstrated that when used with other clinical factors and family history, this information could be helpful in identifying appropriate indication for ICD implantation.

Activated fibroblasts drive cellular interactions in end-stage pediatric hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is a relatively rare but debilitating diagnosis in the pediatric population and patients with end-stage HCM require heart transplantation. In this study, we performed single-nucleus RNA sequencing on pediatric HCM and control myocardium. We identified distinct underling cellular processes in pediatric, end-stage HCM in cardiomyocytes, fibroblasts, endothelial cells, and myeloid cells, compared to controls. Pediatric HCM was enriched in cardiomyocytes exhibiting "stressed" myocardium gene signatures and underlying pathways associated with cardiac hypertrophy. Cardiac fibroblasts exhibited clear activation signatures and heightened downstream processes associated with fibrosis, more so than adult counterparts. There was notable depletion of tissue-resident macrophages, and increased vascular remodeling in endothelial cells. Our analysis provides the first single nuclei analysis focused on end-stage pediatric HCM.

Spatial multi-omics: novel tools to study the complexity of cardiovascular diseases

Spatial multi-omic studies have emerged as a promising approach to comprehensively analyze cells in tissues, enabling the joint analysis of multiple data modalities like transcriptome, epigenome, proteome, and metabolome in parallel or even the same tissue section. This review focuses on the recent advancements in spatial multi-omics technologies, including novel data modalities and computational approaches. We discuss the advancements in low-resolution and high-resolution spatial multi-omics methods which can resolve up to 10,000 of individual molecules at subcellular level. By applying and integrating these techniques, researchers have recently gained valuable insights into the molecular circuits and mechanisms which govern cell biology along the cardiovascular disease spectrum. We provide an overview of current data analysis approaches, with a focus on data integration of multi-omic datasets, highlighting strengths and weaknesses of various computational pipelines. These tools play a crucial role in analyzing and interpreting spatial multi-omics datasets, facilitating the discovery of new findings, and enhancing translational cardiovascular research. Despite nontrivial challenges, such as the need for standardization of experimental setups, data analysis, and improved computational tools, the application of spatial multi-omics holds tremendous potential in revolutionizing our understanding of human disease processes and the identification of novel biomarkers and therapeutic targets. Exciting opportunities lie ahead for the spatial multi-omics field and will likely contribute to the advancement of personalized medicine for cardiovascular diseases.

Myocardial B cells have specific gene expression and predicted interactions in Dilated Cardiomyopathy and Arrhythmogenic Right Ventricular Cardiomyopathy

Introduction

Growing evidence from animal models indicates that the myocardium hosts a population of B cells that play a role in the development of cardiomyopathy. However, there is minimal data on human myocardial B cells in the context of cardiomyopathy.

Methods

We integrated single-cell and single-nuclei datasets from 45 healthy human hearts, 70 hearts with dilated cardiomyopathy (DCM), and 8 hearts with Arrhythmogenic Right Ventricular Cardiomyopathy (ARVC). Interactions between B cells and other cell types were investigated using the CellChat Package. Differential gene expression analysis comparing B cells across conditions was performed using DESeq2. Pathway analysis was performed using Ingenuity, KEGG, and GO pathways analysis.

Results

We identified 1,100 B cells, including naive B cells and plasma cells. B cells showed an extensive network of interactions within the healthy myocardium that included outgoing signaling to macrophages, T cells, endothelial cells, and pericytes, and incoming signaling from endothelial cells, pericytes, and fibroblasts. This niche relied on ECM-receptor, contact, and paracrine interaction; and changed significantly in the context of cardiomyopathy, displaying disease-specific features. Differential gene expression analysis showed that in the context of DCM both naive and plasma B cells upregulated several pathways related to immune activation, including upregulation of oxidative phosphorylation, upregulation of leukocyte extravasation, and, in naive B cells, antigen presentation.

Discussion

The human myocardium contains naive B cells and plasma cells, integrated into a diverse and dynamic niche that has distinctive features in healthy myocardium, DCM, and ARVC. Naive myocardial-associated B cells likely contribute to the pathogenesis of human DCM.

CardiOmics signatures reveal therapeutically actionable targets and drugs for cardiovascular diseases

Cardiovascular diseases are the leading cause of death worldwide, with heart failure being a complex condition that affects millions of individuals. Single-nucleus RNA sequencing has recently emerged as a powerful tool for unraveling the molecular mechanisms behind cardiovascular diseases. This cutting-edge technology enables the identification of molecular signatures, intracellular networks, and spatial relationships among cardiac cells, including cardiomyocytes, mast cells, lymphocytes, macrophages, lymphatic endothelial cells, endocardial cells, endothelial cells, epicardial cells, adipocytes, fibroblasts, neuronal cells, pericytes, and vascular smooth muscle cells. Despite these advancements, the discovery of essential therapeutic targets and drugs for precision cardiology remains a challenge. To bridge this gap, we conducted comprehensive in silico analyses of single-nucleus RNA sequencing data, functional enrichment, protein interactome network, and identification of the shortest pathways to physiological phenotypes. This integrated multi-omics analysis generated CardiOmics signatures, which allowed us to pinpoint three therapeutically actionable targets (ADRA1A1, PPARG, and ROCK2) and 15 effective drugs, including adrenergic receptor agonists, adrenergic receptor antagonists, norepinephrine precursors, PPAR receptor agonists, and Rho-associated kinase inhibitors, involved in late-stage cardiovascular disease clinical trials.