Gene expression networks in endothelial cellsfrom failing human hearts.

IntroductionChronic heart failure is associated with adverse remodeling of the heart that is typically characterized by cardiomyocyte hypertrophy. This requires the formation of new capillaries to maintain oxygen supply. Insufficient angiogenesis promotes the transition from compensated hypertrophy into heart failure. The aim of the present study was to identify angiogenesis-related gene networks and corresponding regulatory hubs in endothelial cells from failing human hearts.Methods and resultsWe isolated left ventricular endothelial cells from patients with advanced heart failure undergoing left ventricular assist device surgery (n = 15) and healthy organ donors (n = 2) and performed RNA sequencing. Subgroup analysis revealed no impact of co-morbidities on gene expression. In a weighted coexpression network analysis, we found 26 gene clusters, of which nine clusters showed a significant positive or negative correlation with the presence of heart failure. We identified the transcription factors CASZ1 (castor zinc finger 1), ZNF523 (zinc finger protein 523), and NFE2L1 (nuclear factor erythroid 2-related factor 1) as hub genes of a cluster related to angiogenesis. Knock-down of CASZ1, ZNF523, or NFE2L1 in human umbilical vein endothelial cells led to a downregulation of genes from the respective cluster, including CD34 and platelet derived growth factor ß, confirming their regulatory function.ConclusionIn conclusion, we assessed gene networks in endothelial cells and identified transcription factors CASZ1, ZNF532, and NFE2L1 as potential regulators of angiogenesis in failing human hearts. Our study provides insights into the transcriptional regulation of angiogenesis beyond the classical vascular endothelial growth factor signaling pathway.

Atlas of cardiac endothelial cell enhancer elements linking the mineralocorticoid receptor to pathological gene expression.

Endothelial cells play crucial roles in physiology and are increasingly recognized as therapeutic targets in cardiovascular disease. Here, we analyzed the regulatory landscape of cardiac endothelial cells by assessing chromatin accessibility, histone modifications, and 3D chromatin organization and confirmed the functional relevance of enhancer-promoter interactions by CRISPRi-mediated enhancer silencing. We used this dataset to explore mechanisms of transcriptional regulation in cardiovascular disease and compared six different experimental models of heart failure, hypertension, or diabetes. Enhancers that regulate gene expression in diseased endothelial cells were enriched with binding sites for a distinct set of transcription factors, including the mineralocorticoid receptor (MR), a known drug target in heart failure and hypertension. For proof of concept, we applied endothelial cell-specific MR deletion in mice to confirm MR-dependent gene expression and predicted direct MR target genes. Overall, we have compiled here a comprehensive atlas of cardiac endothelial cell enhancer elements that provides insight into the role of transcription factors in cardiovascular disease.

The challenges of research data management in cardiovascular science: a DGK and DZHK position paper-executive summary.

The sharing and documentation of cardiovascular research data are essential for efficient use and reuse of data, thereby aiding scientific transparency, accelerating the progress of cardiovascular research and healthcare, and contributing to the reproducibility of research results. However, challenges remain. This position paper, written on behalf of and approved by the German Cardiac Society and German Centre for Cardiovascular Research, summarizes our current understanding of the challenges in cardiovascular research data management (RDM). These challenges include lack of time, awareness, incentives, and funding for implementing effective RDM; lack of standardization in RDM processes; a need to better identify meaningful and actionable data among the increasing volume and complexity of data being acquired; and a lack of understanding of the legal aspects of data sharing. While several tools exist to increase the degree to which data are findable, accessible, interoperable, and reusable (FAIR), more work is needed to lower the threshold for effective RDM not just in cardiovascular research but in all biomedical research, with data sharing and reuse being factored in at every stage of the scientific process. A culture of open science with FAIR research data should be fostered through education and training of early-career and established research professionals. Ultimately, FAIR RDM requires permanent, long-term effort at all levels. If outcomes can be shown to be superior and to promote better (and better value) science, modern RDM will make a positive difference to cardiovascular science and practice. The full position paper is available in the supplementary materials.

Endothelial deletion of the cytochrome P450 reductase leads to cardiac remodelling.

The cytochrome P450 reductase (POR) transfers electrons to all microsomal cytochrome P450 enzymes (CYP450) thereby driving their activity. In the vascular system, the POR/CYP450 system has been linked to the production of epoxyeicosatrienoic acids (EETs) but also to the generation of reactive oxygen species. In cardiac myocytes (CMs), EETs have been shown to modulate the cardiac function and have cardioprotective effects. The functional importance of the endothelial POR/CYP450 system in the heart is unclear and was studied here using endothelial cell-specific, inducible knockout mice of POR (ecPOR-/-). RNA sequencing of murine cardiac cells revealed a cell type-specific expression of different CYP450 homologues. Cardiac endothelial cells mainly expressed members of the CYP2 family which produces EETs, and of the CYP4 family that generates omega fatty acids. Tamoxifen-induced endothelial deletion of POR in mice led to cardiac remodelling under basal conditions, as shown by an increase in heart weight to body weight ratio and an increased CM area as compared to control animals. Endothelial deletion of POR was associated with a significant increase in endothelial genes linked to protein synthesis with no changes in genes of the oxidative stress response. CM of ecPOR-/- mice exhibited attenuated expression of genes linked to mitochondrial function and an increase in genes related to cardiac myocyte contractility. In a model of pressure overload (transverse aortic constriction, TAC with O-rings), ecPOR-/- mice exhibited an accelerated reduction in cardiac output (CO) and stroke volume (SV) as compared to control mice. These results suggest that loss of endothelial POR along with a reduction in EETs leads to an increase in vascular stiffness and loss in cardioprotection, resulting in cardiac remodelling.

HIF1α-AS1 is a DNA:DNA:RNA triplex-forming lncRNA interacting with the HUSH complex.

DNA:DNA:RNA triplexes that are formed through Hoogsteen base-pairing of the RNA in the major groove of the DNA duplex have been observed in vitro, but the extent to which these interactions occur in cells and how they impact cellular functions remains elusive. Using a combination of bioinformatic techniques, RNA/DNA pulldown and biophysical studies, we set out to identify functionally important DNA:DNA:RNA triplex-forming long non-coding RNAs (lncRNA) in human endothelial cells. The lncRNA HIF1α-AS1 was retrieved as a top hit. Endogenous HIF1α-AS1 reduces the expression of numerous genes, including EPH Receptor A2 and Adrenomedullin through DNA:DNA:RNA triplex formation by acting as an adapter for the repressive human silencing hub complex (HUSH). Moreover, the oxygen-sensitive HIF1α-AS1 is down-regulated in pulmonary hypertension and loss-of-function approaches not only result in gene de-repression but also enhance angiogenic capacity. As exemplified here with HIF1α-AS1, DNA:DNA:RNA triplex formation is a functionally important mechanism of trans-acting gene expression control.

DNA topoisomerase inhibition with the HIF inhibitor acriflavine promotes transcription of lncRNAs in endothelial cells.

The transcription factor hypoxia-inducible factor 1 (HIF1) is an important driver of cancer and is therefore an attractive drug target. Acriflavine (ACF) has been suggested to inhibit HIF1, but its mechanism of action is unknown. Here we investigated the interaction of ACF with DNA and long non-coding RNAs (lncRNAs) and its function in human endothelial cells. ACF promoted apoptosis and reduced proliferation, network formation, and angiogenic capacity. It also induced changes in gene expression, as determined by RNA sequencing (RNA-seq), which could not be attributed to specific inhibition of HIF1. A similar response was observed in murine lung endothelial cells. Although ACF increased and decreased a similar number of protein-coding genes, lncRNAs were preferentially upregulated under normoxic and hypoxic conditions. An assay for transposase accessibility with subsequent DNA sequencing (ATAC-seq) demonstrated that ACF induced strong changes in chromatin accessibility at lncRNA promoters. Immunofluorescence showed displacement of DNA:RNA hybrids. Such effects might be due to ACF-mediated topoisomerase inhibition, which was indeed the case, as reflected by DNA unwinding assays. Comparison with other acridine derivatives and topoisomerase inhibitors suggested that the specific function of ACF is an effect of acridinium-class compounds. This study demonstrates that ACF inhibits topoisomerases rather than HIF specifically and that it elicits a unique expression response of lncRNAs.

Nuclear receptor activation shapes spatial genome organization essential for gene expression control: lessons learned from the vitamin D receptor.

Spatial genome organization is tightly controlled by several regulatory mechanisms and is essential for gene expression control. Nuclear receptors are ligand-activated transcription factors that modulate physiological and pathophysiological processes and are primary pharmacological targets. DNA binding of the important loop-forming insulator protein CCCTC-binding factor (CTCF) was modulated by 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3). We performed CTCF HiChIP assays to produce the first genome-wide dataset of CTCF long-range interactions in 1,25(OH)2D3-treated cells, and to determine whether dynamic changes of spatial chromatin interactions are essential for fine-tuning of nuclear receptor signaling. We detected changes in 3D chromatin organization upon vitamin D receptor (VDR) activation at 3.1% of all observed CTCF interactions. VDR binding was enriched at both differential loop anchors and within differential loops. Differential loops were observed in several putative functional roles including TAD border formation, promoter-enhancer looping, and establishment of VDR-responsive insulated neighborhoods. Vitamin D target genes were enriched in differential loops and at their anchors. Secondary vitamin D effects related to dynamic chromatin domain changes were linked to location of downstream transcription factors in differential loops. CRISPR interference and loop anchor deletion experiments confirmed the functional relevance of nuclear receptor ligand-induced adjustments of the chromatin 3D structure for gene expression regulation.

Specificities of Gßγ subunits for the SNARE complex before and after stimulation of α2a-adrenergic receptors.

Ligand binding to G protein­coupled receptors (GPCRs), such as the α2a-adrenergic receptor (α2aAR), results in the activation of heterotrimeric G proteins, which consist of functionally distinct Gα subunits and Gßγ dimers. α2aAR-dependent inhibition of synaptic transmission regulates functions such as spontaneous locomotor activity, anesthetic sparing, and working memory enhancement and requires the soluble NSF attachment protein receptor (SNARE) complex, a Gßγ effector. To understand how the Gßγ-SNARE complex underlies the α2aAR-dependent inhibition of synaptic transmission, we examined the specificity of Gßγ subunits for the SNARE complex in adrenergic neurons, in which auto-α2aARs respond to epinephrine released from these neurons, and nonadrenergic neurons, in which hetero-α2aARs respond to epinephrine released from other neurons. We performed a quantitative, targeted multiple reaction monitoring proteomic analysis of Gß and Gγ subunits bound to the SNARE complex in synaptosomes from mouse brains. In the absence of stimulation of auto-α2aARs, Gß1 and Gγ3 interacted with the SNARE complex. However, Gß1, Gß2, and Gγ3 were found in the complex when auto-α2aARs were activated by epinephrine. Further understanding of the specific usage of distinct Gßγ subunits in vivo may provide insights into the homeostatic regulation of synaptic transmission and the mechanisms of dysfunction that occur in neurological diseases.

Sequential Defects in Cardiac Lineage Commitment and Maturation Cause Hypoplastic Left Heart Syndrome.

BACKGROUND: Complex molecular programs in specific cell lineages govern human heart development. Hypoplastic left heart syndrome (HLHS) is the most common and severe manifestation within the spectrum of left ventricular outflow tract obstruction defects occurring in association with ventricular hypoplasia. The pathogenesis of HLHS is unknown, but hemodynamic disturbances are assumed to play a prominent role. METHODS: To identify perturbations in gene programs controlling ventricular muscle lineage development in HLHS, we performed whole-exome sequencing of 87 HLHS parent-offspring trios, nuclear transcriptomics of cardiomyocytes from ventricles of 4 patients with HLHS and 15 controls at different stages of heart development, single cell RNA sequencing, and 3D modeling in induced pluripotent stem cells from 3 patients with HLHS and 3 controls. RESULTS: Gene set enrichment and protein network analyses of damaging de novo mutations and dysregulated genes from ventricles of patients with HLHS suggested alterations in specific gene programs and cellular processes critical during fetal ventricular cardiogenesis, including cell cycle and cardiomyocyte maturation. Single-cell and 3D modeling with induced pluripotent stem cells demonstrated intrinsic defects in the cell cycle/unfolded protein response/autophagy hub resulting in disrupted differentiation of early cardiac progenitor lineages leading to defective cardiomyocyte subtype differentiation/maturation in HLHS. Premature cell cycle exit of ventricular cardiomyocytes from patients with HLHS prevented normal tissue responses to developmental signals for growth, leading to multinucleation/polyploidy, accumulation of DNA damage, and exacerbated apoptosis, all potential drivers of left ventricular hypoplasia in absence of hemodynamic cues. CONCLUSIONS: Our results highlight that despite genetic heterogeneity in HLHS, many mutations converge on sequential cellular processes primarily driving cardiac myogenesis, suggesting novel therapeutic approaches.

A hierarchical regulatory network analysis of the vitamin D induced transcriptome reveals novel regulators and complete VDR dependency in monocytes.

The transcription factor vitamin D receptor (VDR) is the high affinity nuclear target of the biologically active form of vitamin D3 (1,25(OH)2D3). In order to identify pure genomic transcriptional effects of 1,25(OH)2D3, we used VDR cistrome, transcriptome and open chromatin data, obtained from the human monocytic cell line THP-1, for a novel hierarchical analysis applying three bioinformatics approaches. We predicted 75.6% of all early 1,25(OH)2D3-responding (2.5 or 4 h) and 57.4% of the late differentially expressed genes (24 h) to be primary VDR target genes. VDR knockout led to a complete loss of 1,25(OH)2D3-induced genome-wide gene regulation. Thus, there was no indication of any VDR-independent non-genomic actions of 1,25(OH)2D3 modulating its transcriptional response. Among the predicted primary VDR target genes, 47 were coding for transcription factors and thus may mediate secondary 1,25(OH)2D3 responses. CEBPA and ETS1 ChIP-seq data and RNA-seq following CEBPA knockdown were used to validate the predicted regulation of secondary vitamin D target genes by both transcription factors. In conclusion, a directional network containing 47 partly novel primary VDR target transcription factors describes secondary responses in a highly complex vitamin D signaling cascade. The central transcription factor VDR is indispensable for all transcriptome-wide effects of the nuclear hormone.

Proximity to injury, but neither number of nuclei nor ploidy define pathological adaptation and plasticity in cardiomyocytes.

The adult mammalian heart consists of mononuclear and binuclear cardiomyocytes (CMs) with various ploidies. However, it remains unclear whether a variation in ploidy or number of nuclei is associated with distinct functions and injury responses in CMs, including regeneration. Therefore, we investigated transcriptomes and cellular as well as nuclear features of mononucleated and binucleated CMs in adult mouse hearts with and without injury. To be able to identify the role of ploidy we analyzed control and failing human ventricular CMs because human CMs show a larger and disease-sensitive degree of polyploidization. Using transgenic Myh6-H2BmCh to identify mononucleated and binucleated mouse CMs, we found that cellular volume and RNA content were similar in both. On average nuclei of mononuclear CMs showed a 2-fold higher ploidy, as compared to binuclear CMs indicating that most mononuclear CMs are tetraploid. After myocardial infarction mononucleated and binucleated CMs in the border zone of the lesion responded with hypertrophy and corresponding changes in gene expression, as well as a low level of induction of cell cycle gene expression. Human CMs allowed us to study a wide range of polyploidy spanning from 2n to 16n. Notably, basal as well as pathological gene expression signatures and programs in failing CMs proved to be independent of ploidy. In summary, gene expression profiles were induced in proximity to injury, but independent of number of nuclei or ploidy levels in CMs.

Diabetes changes gene expression but not DNA methylation in cardiac cells.

BACKGROUND: Diabetes mellitus is a worldwide epidemic that causes high mortality due to cardiovascular complications, in particular heart failure. Diabetes is associated with profound pathophysiological changes in the heart. The aim of this study was to investigate the impact of diabetes on gene expression and DNA methylation in cardiac cells. METHODS AND RESULTS: Transcriptome analysis of heart tissue from mice with streptozotocin-induced diabetes revealed only 39 genes regulated, whereas cell type-specific analysis of the diabetic heart was more sensitive and more specific than heart tissue analysis and revealed a total of 3205 differentially regulated genes in five cell types. Whole genome DNA methylation analysis with basepair resolution of distinct cardiac cell types identified highly specific DNA methylation signatures of genic and regulatory regions. Interestingly, despite marked changes in gene expression, DNA methylation remained stable in streptozotocin-induced diabetes. Integrated analysis of cell type-specific gene expression enabled us to assign the particular contribution of single cell types to the pathophysiology of the diabetic heart. Finally, analysis of gene regulation revealed ligand-receptor pairs as potential mediators of heterocellular interaction in the diabetic heart, with fibroblasts and monocytes showing the highest degree of interaction. CONCLUSION: In summary, cell type-specific analysis reveals differentially regulated gene programs that are associated with distinct biological processes in diabetes. Interestingly, despite these changes in gene expression, cell type-specific DNA methylation signatures of genic and regulatory regions remain stable in diabetes. Analysis of heterocellular interactions in the diabetic heart suggest that the interplay between fibroblasts and monocytes is of pivotal importance.

Congenital heart disease risk loci identified by genome-wide association study in European patients.

Genetic factors undoubtedly affect the development of congenital heart disease (CHD) but still remain ill defined. We sought to identify genetic risk factors associated with CHD and to accomplish a functional analysis of SNP-carrying genes. We performed a genome-wide association study (GWAS) of 4034 White patients with CHD and 8486 healthy controls. One SNP on chromosome 5q22.2 reached genome-wide significance across all CHD phenotypes and was also indicative for septal defects. One region on chromosome 20p12.1 pointing to the MACROD2 locus identified 4 highly significant SNPs in patients with transposition of the great arteries (TGA). Three highly significant risk variants on chromosome 17q21.32 within the GOSR2 locus were detected in patients with anomalies of thoracic arteries and veins (ATAV). Genetic variants associated with ATAV are suggested to influence the expression of WNT3, and the variant rs870142 related to septal defects is proposed to influence the expression of MSX1. We analyzed the expression of all 4 genes during cardiac differentiation of human and murine induced pluripotent stem cells in vitro and by single-cell RNA-Seq analyses of developing murine and human hearts. Our data show that MACROD2, GOSR2, WNT3, and MSX1 play an essential functional role in heart development at the embryonic and newborn stages.

ZNF354C is a transcriptional repressor that inhibits endothelial angiogenic sprouting.

Zinc finger proteins (ZNF) are a large group of transcription factors with diverse functions. We recently discovered that endothelial cells harbour a specific mechanism to limit the action of ZNF354C, whose function in endothelial cells is unknown. Given that ZNF354C has so far only been studied in bone and tumour, its function was determined in endothelial cells. ZNF354C is expressed in vascular cells and localises to the nucleus and cytoplasm. Overexpression of ZNF354C in human endothelial cells results in a marked inhibition of endothelial sprouting. RNA-sequencing of human microvascular endothelial cells with and without overexpression of ZNF354C revealed that the protein is a potent transcriptional repressor. ZNF354C contains an active KRAB domain which mediates this suppression as shown by mutagenesis analysis. ZNF354C interacts with dsDNA, TRIM28 and histones, as observed by proximity ligation and immunoprecipitation. Moreover, chromatin immunoprecipitation revealed that the ZNF binds to specific endothelial-relevant target-gene promoters. ZNF354C suppresses these genes as shown by CRISPR/Cas knockout and RNAi. Inhibition of endothelial sprouting by ZNF354C is dependent on the amino acids DV and MLE of the KRAB domain. These results demonstrate that ZNF354C is a repressive transcription factor which acts through a KRAB domain to inhibit endothelial angiogenic sprouting.

Galaxy HiCExplorer 3: a web server for reproducible Hi-C, capture Hi-C and single-cell Hi-C data analysis, quality control and visualization.

The Galaxy HiCExplorer provides a web service at https://hicexplorer.usegalaxy.eu. It enables the integrative analysis of chromosome conformation by providing tools and computational resources to pre-process, analyse and visualize Hi-C, Capture Hi-C (cHi-C) and single-cell Hi-C (scHi-C) data. Since the last publication, Galaxy HiCExplorer has been expanded considerably with new tools to facilitate the analysis of cHi-C and to provide an in-depth analysis of Hi-C data. Moreover, it supports the analysis of scHi-C data by offering a broad range of tools. With the help of the standard graphical user interface of Galaxy, presented workflows, extensive documentation and tutorials, novices as well as Hi-C experts are supported in their Hi-C data analysis with Galaxy HiCExplorer.

Author Correction: The in vivo specificity of synaptic Gß and Gγ subunits to the α2a adrenergic receptor at CNS synapses.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

miR-128a Acts as a Regulator in Cardiac Development by Modulating Differentiation of Cardiac Progenitor Cell Populations.

MicroRNAs (miRs) appear to be major, yet poorly understood players in regulatory networks guiding cardiogenesis. We sought to identify miRs with unknown functions during cardiogenesis analyzing the miR-profile of multipotent Nkx2.5 enhancer cardiac progenitor cells (NkxCE-CPCs). Besides well-known candidates such as miR-1, we found about 40 miRs that were highly enriched in NkxCE-CPCs, four of which were chosen for further analysis. Knockdown in zebrafish revealed that only miR-128a affected cardiac development and function robustly. For a detailed analysis, loss-of-function and gain-of-function experiments were performed during in vitro differentiations of transgenic murine pluripotent stem cells. MiR-128a knockdown (1) increased Isl1, Sfrp5, and Hcn4 (cardiac transcription factors) but reduced Irx4 at the onset of cardiogenesis, (2) upregulated Isl1-positive CPCs, whereas NkxCE-positive CPCs were downregulated, and (3) increased the expression of the ventricular cardiomyocyte marker Myl2 accompanied by a reduced beating frequency of early cardiomyocytes. Overexpression of miR-128a (4) diminished the expression of Isl1, Sfrp5, Nkx2.5, and Mef2c, but increased Irx4, (5) enhanced NkxCE-positive CPCs, and (6) favored nodal-like cardiomyocytes (Tnnt2+, Myh6+, Shox2+) accompanied by increased beating frequencies. In summary, we demonstrated that miR-128a plays a so-far unknown role in early heart development by affecting the timing of CPC differentiation into various cardiomyocyte subtypes.

α2A-adrenergic heteroreceptors are required for stress-induced reinstatement of cocaine conditioned place preference.

The α2a-adrenergic receptor (α2a-AR) agonist guanfacine has been investigated as a potential treatment for substance use disorders. While decreasing stress-induced reinstatement of cocaine seeking in animal models and stress-induced craving in human studies, guanfacine has not been reported to decrease relapse rates. Although guanfacine engages α2a-AR autoreceptors, it also activates excitatory Gi-coupled heteroreceptors in the bed nucleus of the stria terminalis (BNST), a key brain region in driving stress-induced relapse. Thus, BNST α2a-AR heteroreceptor signaling might decrease the beneficial efficacy of guanfacine. We aimed to determine the role of α2a-AR heteroreceptors and BNST Gi-GPCR signaling in stress-induced reinstatement of cocaine conditioned place preference (CPP) and the effects of low dose guanfacine on BNST activity and stress-induced reinstatement. We used a genetic deletion strategy and the cocaine CPP procedure to first define the contributions of α2a-AR heteroreceptors to stress-induced reinstatement. Next, we mimicked BNST Gi-coupled α2a-AR heteroreceptor signaling using a Gi-coupled designer receptor exclusively activated by designer drug (Gi-DREADD) approach. Finally, we evaluated the effects of low-dose guanfacine on BNST cFOS immunoreactivity and stress-induced reinstatement. We show that α2a-AR heteroreceptor deletion disrupts stress-induced reinstatement and that BNST Gi-DREADD activation is sufficient to induce reinstatement. Importantly, we found that low-dose guanfacine does not increase BNST activity, but prevents stress-induced reinstatement. Our findings demonstrate a role for α2a-AR heteroreceptors and BNST Gi-GPCR signaling in stress-induced reinstatement of cocaine CPP and provide insight into the impact of dose on the efficacy of guanfacine as a treatment for stress-induced relapse of cocaine use.