Revealing Chamber-specific Epigenetics in Cardiovascular Disease

Chamber-specific Epigenetic Analysis Identifies Cardiovascular Disease-associated Regulatory Elements in Cardiomyocytes: A new paper describing chamber-specific epigenetic analysis of human adult cardiomyocytes reports the identification of distal cis-regulatory elements associated with cardiovascular disease.

Cardiovascular Disease and the Epigenetic Status of Gene Regulatory Elements

Combined contacts, contact frequency, and histone tracks in cardiovascular cells
Original traces of chromatin interaction frequency heatmap at 5 kb resolution (Hi-C), virtual viewpoint analysis of PLN promoter interactions (1 kb resolution), gene expression (RNA-seq), mCpG (DNA methylation sequencing), histone modification ChIP-seq and histone code (ChromHMM). Low-methylated regions (LMR) and PIDs are annotated. Gray area highlights the PLN promoter. Hi-C data is merged from 9 biological replicates. CRE, cis-regulatory elements. From figure 1 of Haydar et al.

The onset of cardiovascular diseases occurs alongside widespread epigenetic alterations to the distal cis-regulatory elements (CREs, including gene enhancers) that control the gene expression programs underpinning normal heart development, leading to altered gene expression and cell dysfunction. While we generally lack knowledge regarding CRE activity and function in specific cell types or states, a previously reported comparative analysis of mouse cardiac cells (which included cardiomyocytes) described the cell-type-specific responses of CREs to cardiovascular disease (Lother et al.). At the epigenetic level, a joint analysis of histone modifications and DNA methylation in human cardiomyocytes identified over 100,000 CREs, many of which displayed an enrichment for cardiovascular disease-associated non-coding variants (Gilsbach et al., Tan et al., and Hocker et al.). Genetic variants within CREs may exert cell-type-specific effects, with cardiac arrhythmia and hypertrophy known to be associated with CRE variants in distinct cardiomyocyte populations (Kanemaru et al.). These non-coding variants can alter transcription factor binding to CREs, thereby disrupting enhancer-promoter interactions and altering gene expression ([Ameen et al.](https://doi.org/10.1016/j.cell.2022.11.028

With this research in mind, the laboratory of Ralf Gilsbach recently set out to create a high-resolution genome-wide map of chromatin interactions between CREs and target gene promoters in human adult cardiomyocytes from distinct heart chambers, thereby linking them to gene expression programs and genetic risk factors. Their analyses involved high-resolution high-throughput chromosome conformation capture sequencing (Hi-C; Rao et al.) of healthy atrial and ventricular cardiomyocytes and failing ventricular cardiomyocytes and the integration of additional epigenetic and transcriptomic datasets. Fascinatingly, their new Nature Communications study now provides deeper insights into the chamber-specific regulatory mechanisms controlling human cardiomyocyte function and the interpretation of non-coding variants that affect cardiovascular diseases (Haydar et al.).

Paired-Tag technology from Epigenome Technologies generates joint epigenetic and transcriptomic profiles at single-cell resolution and detects histone modifications and RNA transcripts in individual nuclei with efficiency comparable to single-nucleus RNA-seq/ChIP-seq assays. Could integrating Paired-Tag datasets have provided additional robustness to this study, furthering our understanding of chamber- and cell-specific CREs and their links to cardiovascular disease?

Multi-panel image showing changes between chamber-specific heart cells
a Volcano plot illustrating changes in promoter-interacting domains (PIDs) between NF-LV-CM and NF-LA-CM. Differential PIDs were identified using a Chi-square test (P<0.05), and a replicate-based post-hoc analysis was performed (two-sided unpaired t-test). Chromatin interactions with P<0.05 are considered significant. Interactions significantly enriched in NF-LV-CM (blue points) and, NF-LA-CM (green points) are highlighted. Hi-C data is derived from 9 (NF-LV-CM) or 6 (NF-LA-CM) biological replicates. b Correlation between changes in promoter interactions and gene expression in NF-LV-CM and NF-LA-CM. Genes are grouped according to changes in differential interaction. From figure 3 of Haydar et al.

The authors initially isolated bulk cardiac tissue or cardiomyocyte nuclei from three non-failing left ventricles, generated Hi-C data, constructed genome-wide chromatin state models using ChIP-seq data and ChromHMM (Ernst & Kellis), and identified CREs as regions with low CpG methylation (Gilsbach et al., Stadler et al., and Burger et al.). Analysis of this data revealed the highly cell-type-specific nature of chromatin interactions and 3D genome organization in the human heart, thereby suggesting a requirement for chamber-specific cardiomyocyte analysis to accurately “decode” the distinct regulatory architectures present. They then undertook a similar analysis in human cardiomyocytes from non-failing left atria and left ventricles, and from terminally failing ventricles, to support differential chromatin interaction analysis in healthy vs. diseased chamber-specific cardiomyocytes. An initial analysis indicated the highly conserved nature of higher-order genome structure (topologically associating domains and A/B compartments) across human cardiomyocyte subtypes, despite transcriptional and epigenetic differences.

The next team focused on higher-order chromatin structures to study interactions between gene promoters and CREs (typically enhancer elements), and also integrated chromatin state annotations, DNA methylation profiles, and RNA-seq data. This analysis supported the identification of promoter-interacting CREs in chamber-specific cardiomyocytes, which the team functionally validated using CRISPR interference-mediated perturbation (Qi et al.), a method that enables targeted epigenetic silencing without altering the underlying genomic sequence. The study highlighted links between chamber-specific promoter-CRE interactions and gene expression in cardiomyocytes, and a comparative analysis revealed significant alterations in promoter interactions in diseased cardiomyocytes.

Combined contacts, contact frequency, and histone tracks in cardiovascular cells
Virtual viewpoint analysis of the HEY2 (a) and TBX5 promoters (b) at 1kb resolution. Gene expression (RNA-seq), CpG methylation, chromatin accessibility (ATAC-seq), and histone code are shown. NF-LV-CM and NF-LA-CM tracks are shown in blue and green, respectively. Gray areas highlight the HEY2 and TBX5 promoters and the H3K27ac-marked enhancer. Hi-C data were derived from 9 (NF-LV-CM) or 6 (NF-LA-CM) biological replicates. sgRNA target regions for CRISPRi are annotated. AAV-mediated CRISPRi silencing of HEY2 PID1 in hiPSC-vCM. From figure 3 of Haydar et al.

Cardiac-specific CREs generally display the enrichment of cardiac disease-associated genetic variants (see the original study for a range of references). To address the contribution of non-coding genetic variants to atrial- and ventricular-specific cardiovascular diseases, the authors integrated chamber-specific promoter-CRE interaction data, thereby enhancing variant-to-gene mapping and supporting the discovery of additional disease-associated regulatory loci. This analysis underscored the chamber-specific enrichment of variants associated with heart disease traits in active CREs; for example, the study discovered a more significant prevalence of atrial fibrillation-associated variants in promoter-interacting CREs in atrial as compared to ventricular cardiomyocytes. As such, they highlighted the KCNJ2 gene locus (which encodes a potassium channel associated with ventricular arrhythmia susceptibility) and identified uncharacterized CREs harboring variants that affected QT duration. This study (which also included functional epigenetic silencing) identified enhancer elements harboring QT-duration-associated genetic risk factors that modulated gene expression, altered KCNJ2-dependent channel current, and affected repolarization kinetics in cardiomyocytes, thereby providing mechanistic insights into the genetic basis of QT syndrome (a heart rhythm disorder that causes fast, chaotic heartbeats). Of note, additional CREs highlighted included those interacting with the SCN5A, KCNH219, and KCNH2 gene loci; the authors emphasize that only cell-type-specific epigenome and chromatin-interaction profiling can detect promoter interactions in cardiac tissues and, as such, decode the disease relevance of non-coding regulatory elements.

The Next Steps Towards Decoding Cardiovascular Disease

This paper underscored the general utility of combining CRISPR-based gene silencing with cell- and chamber-type-specific chromatin interaction mapping to identify and explore the function of CREs, exemplified by KCNJ2-associated CREs, which harbor genetic variants associated with the onset of ventricular arrhythmias and whose activity controls gene expression and cellular electrophysiological outcomes. The team now hopes that their maps will help to uncover novel genetic variants driving cardiovascular disease.

The implementation of Paired-Tag technology from Epigenome Technologies, which generates joint epigenetic and transcriptomic profiles at single-cell resolution and detects histone modifications and RNA transcripts in individual nuclei with efficiency comparable to single-nucleus RNA-seq/ChIP-seq assays, has the potential to provide deeper insight into such research aims. What more could the simultaneous single-cell analysis of histone modification and transcriptomic profiles tell us about chamber- and cell-specific CREs and their links to cardiovascular disease?