Mapping the developing human cardiac endothelium at single-cell resolution identifies MECOM as a regulator of arteriovenous gene expression

Intramyocardial Sprouting Tip Cells Specify Coronary Arterialization

BACKGROUND

The elaborate patterning of coronary arteries critically supports the high metabolic activity of the beating heart. How coronary endothelial cells coordinate hierarchical vascular remodeling and achieve arteriovenous specification remains largely unknown. Understanding the molecular and cellular cues that pattern coronary arteries is crucial to develop innovative therapeutic strategies that restore functional perfusion within the ischemic heart.

METHODS

Single-cell transcriptomics and histological validation were used to delineate heterogeneous transcriptional states of the developing and mature coronary endothelium with a focus on sprouting endothelium and arterial cell specification. Genetic lineage tracing and high-resolution 3-dimensional imaging were used to characterize the origin and mechanisms of coronary angiogenic sprouting, as well as to fate-map selective endothelial lineages. Integration of single-cell transcriptomic data from ischemic adult mouse hearts and human embryonic data served to assess the conservation of transcriptional states across development, disease, and species.

RESULTS

We discover that coronary arteries originate from cells that have previously transitioned through a specific tip cell phenotype. We identify nonoverlapping intramyocardial and subepicardial tip cell populations with differential gene expression profiles and regulatory pathways. Esm1-lineage tracing confirmed that intramyocardial tip cells selectively contribute to coronary arteries and endocardial tunnels, but not veins. Notably, prearterial cells are detected from development stages to adulthood, increasingly in response to ischemic injury, and in human embryos, suggesting that tip cell-to-artery specification is a conserved mechanism.

CONCLUSIONS

A tip cell-to-artery specification mechanism drives arterialization of the intramyocardial plexus and endocardial tunnels throughout life and is reactivated upon ischemic injury. Differential sprouting programs govern the formation and specification of the venous and arterial coronary plexus.

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.

Genetics of Cardiac Aging Implicate Organ-Specific Variation

Heart structure and function change with age, and the notion that the heart may age faster for some individuals than for others has driven interest in estimating cardiac age acceleration. However, current approaches have limited feature richness (heart measurements; radiomics) or capture extraneous data and therefore lack cardiac specificity (deep learning [DL] on unmasked chest MRI). These technical limitations have been a barrier to efforts to understand genetic contributions to age acceleration. We hypothesized that a video-based DL model provided with heart-masked MRI data would capture a rich yet cardiac-specific representation of cardiac aging. In 61,691 UK Biobank participants, we excluded noncardiac pixels from cardiac MRI and trained a video-based DL model to predict age from one cardiac cycle in the 4-chamber view. We then computed cardiac age acceleration as the bias-corrected prediction of heart age minus the calendar age. Predicted heart age explained 71.1% of variance in calendar age, with a mean absolute error of 3.3 years. Cardiac age acceleration was linked to unfavorable cardiac geometry and systolic and diastolic dysfunction. We also observed links between cardiac age acceleration and diet, decreased physical activity, increased alcohol and tobacco use, and altered levels of 239 serum proteins, as well as adverse brain MRI characteristics. We found cardiac age acceleration to be heritable (h2g 26.6%); a genome-wide association study identified 8 loci related to linked to cardiomyopathy (near TTN, TNS1, LSM3, PALLD, DSP, PLEC, ANKRD1 and MYO18B) and an additional 16 loci (near MECOM, NPR3, KLHL3, HDGFL1, CDKN1A, ELN, SLC25A37, PI15, AP3M1, HMGA2, ADPRHL1, PGAP3, WNT9B, UHRF1 and DOK5). Of the discovered loci, 21 were not previously associated with cardiac age acceleration. Mendelian randomization revealed that lower genetically mediated levels of 6 circulating proteins (MSRA most strongly), as well as greater levels of 5 proteins (LXN most strongly) were associated with cardiac age acceleration, as were greater blood pressure and Lp(a). A polygenic score for cardiac age acceleration predicted earlier onset of arrhythmia, heart failure, myocardial infarction, and mortality. These findings provide a thematic understanding of cardiac age acceleration and suggest that heart- and vascular-specific factors are key to cardiac age acceleration, predominating over a more global aging program.

Generation and characterisation of scalable and stable human pluripotent stem cell-derived microvascular-like endothelial cells for cardiac applications

Coronary microvascular disease (CMD) and its progression towards major adverse coronary events pose a significant health challenge. Accurate in vitro investigation of CMD requires a robust cell model that faithfully represents the cells within the cardiac microvasculature. Human pluripotent stem cell-derived endothelial cells (hPSC-ECs) offer great potential; however, they are traditionally derived via differentiation protocols that are not readily scalable and are not specified towards the microvasculature. Here, we report the development and comprehensive characterisation of a scalable 3D protocol enabling the generation of phenotypically stable cardiac hPSC-microvascular-like ECs (hPSC-CMVECs) and cardiac pericyte-like cells. These were derived by growing vascular organoids within 3D stirred tank bioreactors and subjecting the emerging 3D hPSC-ECs to high-concentration VEGF-A treatment (3DV). Not only did this promote phenotypic stability of the 3DV hPSC-ECs; single cell-RNA sequencing (scRNA-seq) revealed the pronounced expression of cardiac endothelial- and microvascular-associated genes. Further, the generated mural cells attained from the vascular organoid exhibited markers characteristic of cardiac pericytes. Thus, we present a suitable cell model for investigating the cardiac microvasculature as well as the endothelial-dependent and -independent mechanisms of CMD. Moreover, owing to their phenotypic stability, cardiac specificity, and high angiogenic potential, the cells described within would also be well suited for cardiac tissue engineering applications.

Early Injury Landscape in Vein Harvest by Single-Cell and Spatial Transcriptomics

BACKGROUND

Vein graft failure following cardiovascular bypass surgery results in significant patient morbidity and cost to the healthcare system. Vein graft injury can occur during autogenous vein harvest and preparation, as well as after implantation into the arterial system, leading to the development of intimal hyperplasia, vein graft stenosis, and, ultimately, bypass graft failure. Although previous studies have identified maladaptive pathways that occur shortly after implantation, the specific signaling pathways that occur during vein graft preparation are not well defined and may result in a cumulative impact on vein graft failure. We, therefore, aimed to elucidate the response of the vein conduit wall during harvest and following implantation, probing the key maladaptive pathways driving graft failure with the overarching goal of identifying therapeutic targets for biologic intervention to minimize these natural responses to surgical vein graft injury.

METHODS

Employing a novel approach to investigating vascular pathologies, we harnessed both single-nuclei RNA-sequencing and spatial transcriptomics analyses to profile the genomic effects of vein grafts after harvest and distension, then compared these findings to vein grafts obtained 24 hours after carotid-carotid vein bypass implantation in a canine model (n=4).

RESULTS

Spatial transcriptomic analysis of canine cephalic vein after initial conduit harvest and distention revealed significant enrichment of pathways (P<0.05) involved in the activation of endothelial cells (ECs), fibroblasts, and vascular smooth muscle cells, namely pathways responsible for cellular proliferation and migration and platelet activation across the intimal and medial layers, cytokine signaling within the adventitial layer, and ECM (extracellular matrix) remodeling throughout the vein wall. Subsequent single-nuclei RNA-sequencing analysis supported these findings and further unveiled distinct EC and fibroblast subpopulations with significant upregulation (P<0.05) of markers related to endothelial injury response and cellular activation of ECs, fibroblasts, and vascular smooth muscle cells. Similarly, in vein grafts obtained 24 hours after arterial bypass, there was an increase in myeloid cell, protomyofibroblast, injury response EC, and mesenchymal-transitioning EC subpopulations with a concomitant decrease in homeostatic ECs and fibroblasts. Among these markers were genes previously implicated in vein graft injury, including VCAN, FBN1, and VEGFC, in addition to novel genes of interest, such as GLIS3 and EPHA3. These genes were further noted to be driving the expression of genes implicated in vascular remodeling and graft failure, such as IL-6, TGFBR1, SMAD4, and ADAMTS9. By integrating the spatial transcriptomics and single-nuclei RNA-sequencing data sets, we highlighted the spatial architecture of the vein graft following distension, wherein activated and mesenchymal-transitioning ECs, myeloid cells, and fibroblasts were notably enriched in the intima and media of distended veins. Finally, intercellular communication network analysis unveiled the critical roles of activated ECs, mesenchymal-transitioning ECs, protomyofibroblasts, and vascular smooth muscle cells in upregulating signaling pathways associated with cellular proliferation (MDK [midkine], PDGF [platelet-derived growth factor], VEGF [vascular endothelial growth factor]), transdifferentiation (Notch), migration (ephrin, semaphorin), ECM remodeling (collagen, laminin, fibronectin), and inflammation (thrombospondin), following distension.

CONCLUSIONS

Vein conduit harvest and distension elicit a prompt genomic response facilitated by distinct cellular subpopulations heterogeneously distributed throughout the vein wall. This response was found to be further exacerbated following vein graft implantation, resulting in a cascade of maladaptive gene regulatory networks. Together, these results suggest that distension initiates the upregulation of pathological pathways that may ultimately contribute to bypass graft failure and presents potential early targets warranting investigation for targeted therapies. This work highlights the first applications of single-nuclei and spatial transcriptomic analyses to investigate venous pathologies, underscoring the utility of these methodologies and providing a foundation for future investigations.

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.

Potential molecular and cellular mechanisms of the effects of cuproptosis-related genes in the cardiomyocytes of patients with diabetic heart failure: a bioinformatics analysis

Background

Diabetes mellitus is an independent risk factor for heart failure, and diabetes-induced heart failure severely affects patients' health and quality of life. Cuproptosis is a newly defined type of programmed cell death that is thought to be involved in the pathogenesis and progression of cardiovascular disease, but the molecular mechanisms involved are not well understood. Therefore, we aimed to identify biomarkers associated with cuproptosis in diabetes mellitus-associated heart failure and the potential pathological mechanisms in cardiomyocytes.

Materials

Cuproptosis-associated genes were identified from the previous publication. The GSE26887 dataset was downloaded from the GEO database.

Methods

The consistency clustering was performed according to the cuproptosis gene expression. Differentially expressed genes were identified using the limma package, key genes were identified using the weighted gene co-expression network analysis(WGCNA) method, and these were subjected to immune infiltration analysis, enrichment analysis, and prediction of the key associated transcription factors. Consistency clustering identified three cuproptosis clusters. The differentially expressed genes for each were identified using limma and the most critical MEantiquewhite4 module was obtained using WGCNA. We then evaluated the intersection of the MEantiquewhite4 output with the three clusters, and obtained the key genes.

Results

There were four key genes: HSDL2, BCO2, CORIN, and SNORA80E. HSDL2, BCO2, and CORIN were negatively associated with multiple immune factors, while SNORA80E was positively associated, and T-cells accounted for a major proportion of this relationship with the immune system. Four enriched pathways were found to be associated: arachidonic acid metabolism, peroxisomes, fatty acid metabolism, and dorsoventral axis formation, which may be regulated by the transcription factor MECOM, through a change in protein structure.

Conclusion

HSDL2, BCO2, CORIN, and SNORA80E may regulate cardiomyocyte cuproptosis in patients with diabetes mellitus-associated heart failure through effects on the immune system. The product of the cuproptosis-associated gene LOXL2 is probably involved in myocardial fibrosis in patients with diabetes, which leads to the development of cardiac insufficiency.

Contribution of VEGF-B-Induced Endocardial Endothelial Cell Lineage in Physiological Versus Pathological Cardiac Hypertrophy

BACKGROUND

Preclinical studies have shown the therapeutic potential of VEGF-B (vascular endothelial growth factor B) in revascularization of the ischemic myocardium, but the associated cardiac hypertrophy and adverse side effects remain a concern. To understand the importance of endothelial proliferation and migration for the beneficial versus adverse effects of VEGF-B in the heart, we explored the cardiac effects of autocrine versus paracrine VEGF-B expression in transgenic and gene-transduced mice.

METHODS

We used single-cell RNA sequencing to compare cardiac endothelial gene expression in VEGF-B transgenic mouse models. Lineage tracing was used to identify the origin of a VEGF-B-induced novel endothelial cell population and adeno-associated virus-mediated gene delivery to compare the effects of VEGF-B isoforms. Cardiac function was investigated using echocardiography, magnetic resonance imaging, and micro-computed tomography.

RESULTS

Unlike in physiological cardiac hypertrophy driven by a cardiomyocyte-specific VEGF-B transgene (myosin heavy chain alpha-VEGF-B), autocrine VEGF-B expression in cardiac endothelium (aP2 [adipocyte protein 2]-VEGF-B) was associated with septal defects and failure to increase perfused subendocardial capillaries postnatally. Paracrine VEGF-B led to robust proliferation and myocardial migration of a novel cardiac endothelial cell lineage (VEGF-B-induced endothelial cells) of endocardial origin, whereas autocrine VEGF-B increased proliferation of VEGF-B-induced endothelial cells but failed to promote their migration and efficient contribution to myocardial capillaries. The surviving aP2-VEGF-B offspring showed an altered ratio of secreted VEGF-B isoforms and developed massive pathological cardiac hypertrophy with a distinct cardiac vessel pattern. In the normal heart, we found a small VEGF-B-induced endothelial cell population that was only minimally expanded during myocardial infarction but not during physiological cardiac hypertrophy associated with mouse pregnancy.

CONCLUSIONS

Paracrine and autocrine secretions of VEGF-B induce expansion of a specific endocardium-derived endothelial cell population with distinct angiogenic markers. However, autocrine VEGF-B signaling fails to promote VEGF-B-induced endothelial cell migration and contribution to myocardial capillaries, predisposing to septal defects and inducing a mismatch between angiogenesis and myocardial growth, which results in pathological cardiac hypertrophy.

Disruption of maternal vascular remodeling by a fetal endoretrovirus-derived gene in preeclampsia

BACKGROUND

Preeclampsia, one of the most lethal pregnancy-related diseases, is associated with the disruption of uterine spiral artery remodeling during placentation. However, the early molecular events leading to preeclampsia remain unknown.

RESULTS

By analyzing placentas from preeclampsia, non-preeclampsia, and twin pregnancies with selective intrauterine growth restriction, we show that the pathogenesis of preeclampsia is attributed to immature trophoblast and maldeveloped endothelial cells. Delayed epigenetic reprogramming during early extraembryonic tissue development leads to generation of excessive immature trophoblast cells. We find reduction of de novo DNA methylation in these trophoblast cells results in selective overexpression of maternally imprinted genes, including the endoretrovirus-derived gene PEG10 (paternally expressed gene 10). PEG10 forms virus-like particles, which are transferred from the trophoblast to the closely proximate endothelial cells. In normal pregnancy, only a low amount of PEG10 is transferred to maternal cells; however, in preeclampsia, excessive PEG10 disrupts maternal vascular development by inhibiting TGF-beta signaling.

CONCLUSIONS

Our study reveals the intricate epigenetic mechanisms that regulate trans-generational genetic conflict and ultimately ensure proper maternal-fetal interface formation.

Stem cell-derived vessels-on-chip for cardiovascular disease modeling

The metabolic syndrome is accompanied by vascular complications. Human in vitro disease models are hence required to better understand vascular dysfunctions and guide clinical therapies. Here, we engineered an open microfluidic vessel-on-chip platform that integrates human pluripotent stem cell-derived endothelial cells (SC-ECs). The open microfluidic design enables seamless integration with state-of-the-art analytical technologies, including single-cell RNA sequencing, proteomics by mass spectrometry, and high-resolution imaging. Beyond previous systems, we report SC-EC maturation by means of barrier formation, arterial toning, and high nitric oxide synthesis levels under gravity-driven flow. Functionally, we corroborate the hallmarks of early-onset atherosclerosis with low sample volumes and cell numbers under flow conditions by determining proteome and secretome changes in SC-ECs stimulated with oxidized low-density lipoprotein and free fatty acids. More broadly, our organ-on-chip platform enables the modeling of patient-specific human endothelial tissue and has the potential to become a general tool for animal-free vascular research.

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.

Cationic cholesterol-dependent LNP delivery to lung stem cells, the liver, and heart

Adding a cationic helper lipid to a lipid nanoparticle (LNP) can increase lung delivery and decrease liver delivery. However, it remains unclear whether charge-dependent tropism is universal or, alternatively, whether it depends on the component that is charged. Here, we report evidence that cationic cholesterol-dependent tropism can differ from cationic helper lipid-dependent tropism. By testing how 196 LNPs delivered mRNA to 22 cell types, we found that charged cholesterols led to a different lung:liver delivery ratio than charged helper lipids. We also found that combining cationic cholesterol with a cationic helper lipid led to mRNA delivery in the heart as well as several lung cell types, including stem cell-like populations. These data highlight the utility of exploring charge-dependent LNP tropism.

Integrated single-nuclei and spatial transcriptomic analysis reveals propagation of early acute vein harvest and distension injury signaling pathways following arterial implantation

Background

Vein graft failure (VGF) following cardiovascular bypass surgery results in significant patient morbidity and cost to the healthcare system. Vein graft injury can occur during autogenous vein harvest and preparation, as well as after implantation into the arterial system, leading to the development of intimal hyperplasia, vein graft stenosis, and, ultimately, bypass graft failure. While previous studies have identified maladaptive pathways that occur shortly after implantation, the specific signaling pathways that occur during vein graft preparation are not well defined and may result in a cumulative impact on VGF. We, therefore, aimed to elucidate the response of the vein conduit wall during harvest and following implantation, probing the key maladaptive pathways driving graft failure with the overarching goal of identifying therapeutic targets for biologic intervention to minimize these natural responses to surgical vein graft injury.

Methods

Employing a novel approach to investigating vascular pathologies, we harnessed both single-nuclei RNA-sequencing (snRNA-seq) and spatial transcriptomics (ST) analyses to profile the genomic effects of vein grafts after harvest and distension, then compared these findings to vein grafts obtained 24 hours after carotid-cartoid vein bypass implantation in a canine model (n=4).

Results

Spatial transcriptomic analysis of canine cephalic vein after initial conduit harvest and distention revealed significant enrichment of pathways (P < 0.05) involved in the activation of endothelial cells (ECs), fibroblasts (FBs), and vascular smooth muscle cells (VSMCs), namely pathways responsible for cellular proliferation and migration and platelet activation across the intimal and medial layers, cytokine signaling within the adventitial layer, and extracellular matrix (ECM) remodeling throughout the vein wall. Subsequent snRNA-seq analysis supported these findings and further unveiled distinct EC and FB subpopulations with significant upregulation (P < 0.00001) of markers related to endothelial injury response and cellular activation of ECs, FBs, and VSMCs. Similarly, in vein grafts obtained 24 hours after arterial bypass, there was an increase in myeloid cell, protomyofibroblast, injury-response EC, and mesenchymal-transitioning EC subpopulations with a concomitant decrease in homeostatic ECs and fibroblasts. Among these markers were genes previously implicated in vein graft injury, including VCAN (versican), FBN1 (fibrillin-1), and VEGFC (vascular endothelial growth factor C), in addition to novel genes of interest such as GLIS3 (GLIS family zinc finger 3) and EPHA3 (ephrin-A3). These genes were further noted to be driving the expression of genes implicated in vascular remodeling and graft failure, such as IL-6, TGFBR1, SMAD4, and ADAMTS9. By integrating the ST and snRNA-seq datasets, we highlighted the spatial architecture of the vein graft following distension, wherein activated and mesenchymal-transitioning ECs, myeloid cells, and FBs were notably enriched in the intima and media of distended veins. Lastly, intercellular communication network analysis unveiled the critical roles of activated ECs, mesenchymal transitioning ECs, protomyofibroblasts, and VSMCs in upregulating signaling pathways associated with cellular proliferation (MDK, PDGF, VEGF), transdifferentiation (Notch), migration (ephrin, semaphorin), ECM remodeling (collagen, laminin, fibronectin), and inflammation (thrombospondin), following distension.

Conclusions

Vein conduit harvest and distension elicit a prompt genomic response facilitated by distinct cellular subpopulations heterogeneously distributed throughout the vein wall. This response was found to be further exacerbated following vein graft implantation, resulting in a cascade of maladaptive gene regulatory networks. Together, these results suggest that distension initiates the upregulation of pathological pathways that may ultimately contribute to bypass graft failure and presents potential early targets warranting investigation for targeted therapies. This work highlights the first applications of single-nuclei and spatial transcriptomic analyses to investigate venous pathologies, underscoring the utility of these methodologies and providing a foundation for future investigations.

Control of coronary vascular cell fate in development and regeneration

The coronary vasculature consists of a complex hierarchal network of arteries, veins, and capillaries which collectively function to perfuse the myocardium. However, the pathways controlling the temporally and spatially restricted mechanisms underlying the formation of this vascular network remain poorly understood. In recent years, the increasing use and refinement of transgenic mouse models has played an instrumental role in offering new insights into the cellular origins of the coronary vasculature, as well as identifying a continuum of transitioning cell states preceding the full maturation of the coronary vasculature. Coupled with the emergence of single cell RNA sequencing platforms, these technologies have begun to uncover the key regulatory factors mediating the convergence of distinct cellular origins to ensure the formation of a collectively functional, yet phenotypically diverse, vascular network. Furthermore, improved understanding of the key regulatory factors governing coronary vessel formation in the embryo may provide crucial clues into future therapeutic strategies to reactivate these developmentally functional mechanisms to drive the revascularisation of the ischaemic adult heart.

Molecular and Spatial Signatures of Mouse Embryonic Endothelial Cells at Single-Cell Resolution

BACKGROUND

Mature endothelial cells (ECs) are heterogeneous, with subtypes defined by tissue origin and position within the vascular bed (ie, artery, capillary, vein, and lymphatic). How this heterogeneity is established during the development of the vascular system, especially arteriovenous specification of ECs, remains incompletely characterized.

METHODS

We used droplet-based single-cell RNA sequencing and multiplexed error-robust fluorescence in situ hybridization to define EC and EC progenitor subtypes from E9.5, E12.5, and E15.5 mouse embryos. We used trajectory inference to analyze the specification of arterial ECs (aECs) and venous ECs (vECs) from EC progenitors. Network analysis identified candidate transcriptional regulators of arteriovenous differentiation, which we tested by CRISPR (clustered regularly interspaced short palindromic repeats) loss of function in human-induced pluripotent stem cells undergoing directed differentiation to aECs or vECs (human-induced pluripotent stem cell-aECs or human-induced pluripotent stem cell-vECs).

RESULTS

From the single-cell transcriptomes of 7682 E9.5 to E15.5 ECs, we identified 19 EC subtypes, including Etv2+Bnip3+ EC progenitors. Spatial transcriptomic analysis of 15 448 ECs provided orthogonal validation of these EC subtypes and established their spatial distribution. Most embryonic ECs were grouped by their vascular-bed types, while ECs from the brain, heart, liver, and lung were grouped by their tissue origins. Arterial (Eln, Dkk2, Vegfc, and Egfl8), venous (Fam174b and Clec14a), and capillary (Kcne3) marker genes were identified. Compared with aECs, embryonic vECs and capillary ECs shared fewer markers than their adult counterparts. Early capillary ECs with venous characteristics functioned as a branch point for differentiation of aEC and vEC lineages.

CONCLUSIONS

Our results provide a spatiotemporal map of embryonic EC heterogeneity at single-cell resolution and demonstrate that the diversity of ECs in the embryo arises from both tissue origin and vascular-bed position. Developing aECs and vECs share common venous-featured capillary precursors and are regulated by distinct transcriptional regulatory networks.

Cell atlas of the Atlantic salmon spleen reveals immune cell heterogeneity and cell-specific responses to bacterial infection

The spleen is a conserved secondary lymphoid organ that emerged in parallel to adaptive immunity in early jawed vertebrates. Recent studies have applied single cell transcriptomics to reveal the cellular composition of spleen in several species, cataloguing diverse immune cell types and subpopulations. In this study, 51,119 spleen nuclei transcriptomes were comprehensively investigated in the commercially important teleost Atlantic salmon (Salmo salar L.), contrasting control animals with those challenged with the bacterial pathogen Aeromonas salmonicida. We identified clusters of nuclei representing the expected major cell types, namely T cells, B cells, natural killer-like cells, granulocytes, mononuclear phagocytes, endothelial cells, mesenchymal cells, erythrocytes and thrombocytes. We discovered heterogeneity within several immune lineages, providing evidence for resident macrophages and melanomacrophages, infiltrating monocytes, several candidate dendritic cell subpopulations, and B cells at distinct stages of differentiation, including plasma cells and an igt + subset. We provide evidence for twelve candidate T cell subsets, including cd4+ T helper and regulatory T cells, one cd8+ subset, three γδT subsets, and populations double negative for cd4 and cd8. The number of genes showing differential expression during the early stages of Aeromonas infection was highly variable across immune cell types, with the largest changes observed in macrophages and infiltrating monocytes, followed by resting mature B cells. Our analysis provides evidence for a local inflammatory response to infection alongside B cell maturation in the spleen, and upregulation of ccr9 genes in igt + B cells, T helper and cd8+ cells, and monocytes, consistent with the recruitment of immune cell populations to the gut to deal with Aeromonas infection. Overall, this study provides a new cell-resolved perspective of the immune actions of Atlantic salmon spleen, highlighting extensive heterogeneity hidden to bulk transcriptomics. We further provide a large catalogue of cell-specific marker genes that can be leveraged to further explore the function and structural organization of the salmonid immune system.

Transcriptional regulators of arterial and venous identity in the developing mammalian embryo

The complex and hierarchical vascular network of arteries, veins, and capillaries features considerable endothelial heterogeneity, yet the regulatory pathways directing arteriovenous specification, differentiation, and identity are still not fully understood. Recent advances in analysis of endothelial-specific gene-regulatory elements, single-cell RNA sequencing, and cell lineage tracing have both emphasized the importance of transcriptional regulation in this process and shed considerable light on the mechanism and regulation of specification within the endothelium. In this review, we discuss recent advances in our understanding of how endothelial cells acquire arterial and venous identity and the role different transcription factors play in this process.

Single-nuclei multiomic analyses identify human cardiac lymphatic endothelial cells associated with coronary arteries in the epicardium

Cardiac lymphatic vessels play important roles in fluid homeostasis, inflammation, disease, and regeneration of the heart. The developing cardiac lymphatics in human fetal hearts are closely associated with coronary arteries, similar to those in zebrafish hearts. We identify a population of cardiac lymphatic endothelial cells (LECs) that reside in the epicardium. Single-nuclei multiomic analysis of the human fetal heart reveals the plasticity and heterogeneity of the cardiac endothelium. Furthermore, we find that VEGFC is highly expressed in arterial endothelial cells and epicardium-derived cells, providing a molecular basis for the arterial association of cardiac lymphatic development. Using a cell-type-specific integrative analysis, we identify a population of cardiac lymphatic endothelial cells marked by the PROX1 and the lymphangiocrine RELN and enriched in binding motifs of erythroblast transformation specific (ETS) variant (ETV) transcription factors. We report the in vivo molecular characterization of human cardiac lymphatics and provide a valuable resource to understand fetal heart development.

MECOM Deficiency: from Bone Marrow Failure to Impaired B-Cell Development

MECOM deficiency is a recently identified inborn error of immunity and inherited bone marrow failure syndrome caused by haploinsufficiency of the hematopoietic transcription factor MECOM. It is unique among inherited bone marrow failure syndromes, many of which present during later childhood or adolescence, because of the early age of onset and severity of the pancytopenia, emphasizing the importance and gene dose dependency of MECOM during hematopoiesis. B-cell lymphopenia and hypogammaglobulinemia have been described in a subset of patients with MECOM deficiency. While the mechanisms underlying the B-cell deficiency are currently unknown, recent work has provided mechanistic insights into the function of MECOM in hematopoietic stem cell (HSC) maintenance. MECOM binds to regulatory enhancers that control the expression of a network of genes essential for HSC maintenance and self-renewal. Heterozygous mutations, as seen in MECOM-deficient bone marrow failure, lead to dysregulated MECOM network expression. Extra-hematopoietic manifestations of MECOM deficiency, including renal and cardiac anomalies, radioulnar synostosis, clinodactyly, and hearing loss, have been reported. Individuals with specific genotypes have some of the systemic manifestations with isolated mild thrombocytopenia or without hematologic abnormalities, highlighting the tissue specificity of mutations in some MECOM domains. Those infants with MECOM-associated bone marrow failure require HSC transplantation for survival. Here, we review the expanding cohort of patient phenotypes and accompanying genotypes resulting in MECOM deficiency, and the proposed mechanisms underlying MECOM regulation of human HSC maintenance and B-cell development.

ScRNA-seq and spatial transcriptomics: exploring the occurrence and treatment of coronary-related diseases starting from development

Single-cell RNA sequencing (scRNA-seq) is a new technology that can be used to explore molecular changes in complex cell clusters at the single-cell level. Single-cell spatial transcriptomic technology complements the cell-space location information lost during single-cell sequencing. Coronary artery disease is an important cardiovascular disease with high mortality rates. Many studies have explored the physiological development and pathological changes in coronary arteries from the perspective of single cells using single-cell spatial transcriptomic technology. This article reviews the molecular mechanisms underlying coronary artery development and diseases as revealed by scRNA-seq combined with spatial transcriptomic technology. Based on these mechanisms, we discuss the possible new treatments for coronary diseases.

The tissue-specific transcriptional landscape underlines the involvement of endothelial cells in health and disease

Endothelial cells (ECs) that line vascular and lymphatic vessels are being increasingly recognized as important to organ function in health and disease. ECs participate not only in the trafficking of gases, metabolites, and cells between the bloodstream and tissues but also in the angiocrine-based induction of heterogeneous parenchymal cells, which are unique to their specific tissue functions. The molecular mechanisms regulating EC heterogeneity between and within different tissues are modeled during embryogenesis and become fully established in adults. Any changes in adult tissue homeostasis induced by aging, stress conditions, and various noxae may reshape EC heterogeneity and induce specific transcriptional features that condition a functional phenotype. Heterogeneity is sustained via specific genetic programs organized through the combinatory effects of a discrete number of transcription factors (TFs) that, at the single tissue-level, constitute dynamic networks that are post-transcriptionally and epigenetically regulated. This review is focused on outlining the TF-based networks involved in EC specialization and physiological and pathological stressors thought to modify their architecture.

Multiomics uncover the proinflammatory role of Trpm4 deletion after myocardial infarction in mice

Upon myocardial infarction (MI), ischemia-induced cell death triggers an inflammatory response responsible for removing necrotic material and inducing tissue repair. TRPM4 is a Ca2+-activated ion channel permeable to monovalent cations. Although its role in cardiomyocyte-driven hypertrophy and arrhythmia post-MI has been established, no study has yet investigated its role in the inflammatory process orchestrated by endothelial cells, immune cells, and fibroblasts. This study aims to assess the role of TRPM4 in 1) survival and cardiac function, 2) inflammation, and 3) healing post-MI. We performed ligation of the left coronary artery or sham intervention on 154 Trpm4 WT or KO mice under isoflurane anesthesia. Survival and echocardiographic functions were monitored up to 5 wk. We collected serum during the acute post-MI phase to analyze proteomes and performed single-cell RNA sequencing on nonmyocytic cells of hearts after 24 and 72 h. Lastly, we assessed chronic fibrosis and angiogenesis. We observed no significant differences in survival or cardiac function, even though our proteomics data showed significantly decreased tissue injury markers (i.e., creatine kinase M and VE-cadherin) in KO serum after 12 h. On the other hand, inflammation, characterized by serum amyloid P component in the serum, higher number of recruited granulocytes, inflammatory monocytes, and macrophages, as well as expression of proinflammatory genes, was significantly higher in KO. This correlated with increased chronic cardiac fibrosis and angiogenesis. Since inflammation and fibrosis are closely linked to adverse remodeling, future therapeutic attempts at inhibiting TRPM4 will need to assess these parameters carefully before proceeding with translational studies.NEW & NOTEWORTHY Deletion of Trpm4 increases markers of cardiac and systemic inflammation within the first 24 h after MI, while inducing an earlier fibrotic transition at 72 h and more overall chronic fibrosis and angiogenesis at 5 wk. The descriptive, robust, and methodologically broad approach of this study sheds light on an important caveat that will need to be taken into account in all future therapeutic attempts to inhibit TRPM4 post-MI.

Multiancestry Genome-Wide Association Study of Aortic Stenosis Identifies Multiple Novel Loci in the Million Veteran Program

BACKGROUND

Calcific aortic stenosis (CAS) is the most common valvular heart disease in older adults and has no effective preventive therapies. Genome-wide association studies (GWAS) can identify genes influencing disease and may help prioritize therapeutic targets for CAS.

METHODS

We performed a GWAS and gene association study of 14 451 patients with CAS and 398 544 controls in the Million Veteran Program. Replication was performed in the Million Veteran Program, Penn Medicine Biobank, Mass General Brigham Biobank, BioVU, and BioMe, totaling 12 889 cases and 348 094 controls. Causal genes were prioritized from genome-wide significant variants using polygenic priority score gene localization, expression quantitative trait locus colocalization, and nearest gene methods. CAS genetic architecture was compared with that of atherosclerotic cardiovascular disease. Causal inference for cardiometabolic biomarkers in CAS was performed using Mendelian randomization and genome-wide significant loci were characterized further through phenome-wide association study.

RESULTS

We identified 23 genome-wide significant lead variants in our GWAS representing 17 unique genomic regions. Of the 23 lead variants, 14 were significant in replication, representing 11 unique genomic regions. Five replicated genomic regions were previously known risk loci for CAS (PALMD, TEX41, IL6, LPA, FADS) and 6 were novel (CEP85L, FTO, SLMAP, CELSR2, MECOM, CDAN1). Two novel lead variants were associated in non-White individuals (P<0.05): rs12740374 (CELSR2) in Black and Hispanic individuals and rs1522387 (SLMAP) in Black individuals. Of the 14 replicated lead variants, only 2 (rs10455872 [LPA], rs12740374 [CELSR2]) were also significant in atherosclerotic cardiovascular disease GWAS. In Mendelian randomization, lipoprotein(a) and low-density lipoprotein cholesterol were both associated with CAS, but the association between low-density lipoprotein cholesterol and CAS was attenuated when adjusting for lipoprotein(a). Phenome-wide association study highlighted varying degrees of pleiotropy, including between CAS and obesity at the FTO locus. However, the FTO locus remained associated with CAS after adjusting for body mass index and maintained a significant independent effect on CAS in mediation analysis.

CONCLUSIONS

We performed a multiancestry GWAS in CAS and identified 6 novel genomic regions in the disease. Secondary analyses highlighted the roles of lipid metabolism, inflammation, cellular senescence, and adiposity in the pathobiology of CAS and clarified the shared and differential genetic architectures of CAS with atherosclerotic cardiovascular diseases.

Vascular endothelial cell development and diversity

Vascular endothelial cells form the inner layer of blood vessels where they have a key role in the development and maintenance of the functional circulatory system and provide paracrine support to surrounding non-vascular cells. Technical advances in the past 5 years in single-cell genomics and in in vivo genetic labelling have facilitated greater insights into endothelial cell development, plasticity and heterogeneity. These advances have also contributed to a new understanding of the timing of endothelial cell subtype differentiation and its relationship to the cell cycle. Identification of novel tissue-specific gene expression patterns in endothelial cells has led to the discovery of crucial signalling pathways and new interactions with other cell types that have key roles in both tissue maintenance and disease pathology. In this Review, we describe the latest findings in vascular endothelial cell development and diversity, which are often supported by large-scale, single-cell studies, and discuss the implications of these findings for vascular medicine. In addition, we highlight how techniques such as single-cell multimodal omics, which have become increasingly sophisticated over the past 2 years, are being utilized to study normal vascular physiology as well as functional perturbations in disease.