Cardiac cell type-specific gene regulatory programs and disease risk association

A Novel Macrophage Subpopulation Conveys Increased Genetic Risk of Coronary Artery Disease

BACKGROUND

Coronary artery disease (CAD), the leading cause of death worldwide, is influenced by both environmental and genetic factors. Although over 250 genetic risk loci have been identified through genome-wide association studies, the specific causal variants and their regulatory mechanisms are still largely unknown, particularly in disease-relevant cell types such as macrophages.

METHODS

We utilized single-cell RNA-seq and single-cell multiomics approaches in primary human monocyte-derived macrophages to explore the transcriptional regulatory network involved in a critical pathogenic event of coronary atherosclerosis-the formation of lipid-laden foam cells. The relative genetic contribution to CAD was assessed by partitioning disease heritability across different macrophage subpopulations. Meta-analysis of single-cell RNA-seq data sets from 38 human atherosclerotic samples was conducted to provide high-resolution cross-referencing to macrophage subpopulations in vivo.

RESULTS

We identified 18 782 cis-regulatory elements by jointly profiling the gene expression and chromatin accessibility of >5000 macrophages. Integration with CAD genome-wide association study data prioritized 121 CAD-related genetic variants and 56 candidate causal genes. We showed that CAD heritability was not uniformly distributed and was particularly enriched in the gene programs of a novel CD52-hi lipid-handling macrophage subpopulation. These CD52-hi macrophages displayed significantly less lipoprotein accumulation and were also found in human atherosclerotic plaques. We investigated the cis-regulatory effect of a risk variant rs10488763 on FDX1, implicating the recruitment of AP-1 and C/EBP-β in the causal mechanisms at this locus.

CONCLUSIONS

Our results provide genetic evidence of the divergent roles of macrophage subsets in atherogenesis and highlight lipid-handling macrophages as a key subpopulation through which genetic variants operate to influence disease. These findings provide an unbiased framework for functional fine-mapping of genome-wide association study results using single-cell multiomics and offer new insights into the genotype-environment interactions underlying atherosclerotic disease.

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.

Single-Nucleus Multiomic Analyses Identifies Gene Regulatory Dynamics of Phenotypic Modulation in Human Aneurysmal Aortic Root

Aortic root aneurysm is a potentially life-threatening condition that may lead to aortic rupture and is often associated with genetic syndromes, such as Marfan syndrome (MFS). Although studies with MFS animal models have provided valuable insights into the pathogenesis of aortic root aneurysms, this understanding of the transcriptomic and epigenomic landscape in human aortic root tissue remains incomplete. This knowledge gap has impeded the development of effective targeted therapies. Here, this study performs the first integrative analysis of single-nucleus multiomic (gene expression and chromatin accessibility) and spatial transcriptomic sequencing data of human aortic root tissue under healthy and MFS conditions. Cell-type-specific transcriptomic and cis-regulatory profiles in the human aortic root are identified. Regulatory and spatial dynamics during phenotypic modulation of vascular smooth muscle cells (VSMCs), the cardinal cell type, are delineated. Moreover, candidate key regulators driving the phenotypic modulation of VSMC, such as FOXN3, TEAD1, BACH2, and BACH1, are identified. In vitro experiments demonstrate that FOXN3 functions as a novel key regulator for maintaining the contractile phenotype of human aortic VSMCs through targeting ACTA2. These findings provide novel insights into the regulatory and spatial dynamics during phenotypic modulation in the aneurysmal aortic root of humans.

Single-nucleus multi-omic profiling of polyploid heart nuclei identifies fusion-derived cardiomyocytes in the human heart

Understanding the mechanisms of polyploidization in cardiomyocytes is crucial for advancing strategies to stimulate myocardial regeneration. Although endoreplication has long been considered the primary source of polyploid human cardiomyocytes, recent animal work suggests the potential for cardiomyocyte fusion. Moreover, the effects of polyploidization on the genomic-transcriptomic repertoire of human cardiomyocytes have not been studied previously. We applied single-nuclei whole genome sequencing, single nuclei RNA sequencing, and multiome ATAC + gene expression (from the same nuclei) techniques to nuclei isolated from 11 healthy hearts. Utilizing post-zygotic non-inherited somatic mutations occurring during development as "endogenous barcodes," to reconstruct lineage relationships of polyploid cardiomyocytes. Of 482 cardiomyocytes from multiple healthy donor hearts 75.7% can be sorted into several developmental clades marked by one or more somatic single-nucleotide variants (SNVs). At least ~10% of tetraploid cardiomyocytes contain cells from distinct clades, indicating fusion of lineally distinct cells, whereas 60% of higher-ploidy cardiomyocytes contain fused cells from distinct clades. Combined snRNA-seq and snATAC-seq revealed transcriptome and chromatin landscapes of polyploid cardiomyocytes distinct from diploid cardiomyocytes, and show some higher-ploidy cardiomyocytes with transcriptional signatures suggesting fusion between cardiomyocytes and endothelial and fibroblast cells. These observations provide the first evidence for cell and nuclear fusion of human cardiomyocytes, raising the possibility that cell fusion may contribute to developing or maintaining polyploid cardiomyocytes in the human heart.

Transcriptomic and spatial dissection of human ex vivo right atrial tissue reveals proinflammatory microvascular changes in ischemic heart disease

Cardiovascular disease plays a central role in the electrical and structural remodeling of the right atrium, predisposing to arrhythmias, heart failure, and sudden death. Here, we dissect with single-nuclei RNA sequencing (snRNA-seq) and spatial transcriptomics the gene expression changes in the human ex vivo right atrial tissue and pericardial fluid in ischemic heart disease, myocardial infarction, and ischemic and non-ischemic heart failure using asymptomatic patients with valvular disease who undergo preventive surgery as the control group. We reveal substantial differences in disease-associated gene expression in all cell types, collectively suggesting inflammatory microvascular dysfunction and changes in the right atrial tissue composition as the valvular and vascular diseases progress into heart failure. The data collectively suggest that investigation of human cardiovascular disease should expand to all functionally important parts of the heart, which may help us to identify mechanisms promoting more severe types of the disease.

Cardiac GR Mediates the Diurnal Rhythm in Ventricular Arrhythmia Susceptibility

BACKGROUND

Ventricular arrhythmias (VAs) demonstrate a prominent day-night rhythm, commonly presenting in the morning. Transcriptional rhythms in cardiac ion channels accompany this phenomenon, but their role in the morning vulnerability to VAs and the underlying mechanisms are not understood. We investigated the recruitment of transcription factors that underpins transcriptional rhythms in ion channels and assessed whether this mechanism was pertinent to the heart's intrinsic diurnal susceptibility to VA.

METHODS AND RESULTS

Assay for transposase-accessible chromatin with sequencing performed in mouse ventricular myocyte nuclei at the beginning of the animals' inactive (ZT0) and active (ZT12) periods revealed differentially accessible chromatin sites annotating to rhythmically transcribed ion channels and distinct transcription factor binding motifs in these regions. Notably, motif enrichment for the glucocorticoid receptor (GR; transcriptional effector of corticosteroid signaling) in open chromatin profiles at ZT12 was observed, in line with the well-recognized ZT12 peak in circulating corticosteroids. Molecular, electrophysiological, and in silico biophysically-detailed modeling approaches demonstrated GR-mediated transcriptional control of ion channels (including Scn5a underlying the cardiac Na+ current, Kcnh2 underlying the rapid delayed rectifier K+ current, and Gja1 responsible for electrical coupling) and their contribution to the day-night rhythm in the vulnerability to VA. Strikingly, both pharmacological block of GR and cardiomyocyte-specific genetic knockout of GR blunted or abolished ion channel expression rhythms and abolished the ZT12 susceptibility to pacing-induced VA in isolated hearts.

CONCLUSIONS

Our study registers a day-night rhythm in chromatin accessibility that accompanies diurnal cycles in ventricular myocytes. Our approaches directly implicate the cardiac GR in the myocyte excitability rhythm and mechanistically link the ZT12 surge in glucocorticoids to intrinsic VA propensity at this time.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

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.

Gender-specific genetic and epigenetic signatures in cardiovascular disease

Cardiac sex differences represent a pertinent focus in pursuit of the long-awaited goal of personalized medicine. Despite evident disparities in the onset and progression of cardiac pathology between sexes, historical oversight has led to the neglect of gender-specific considerations in the treatment of patients. This oversight is attributed to a predominant focus on male samples and a lack of sex-based segregation in patient studies. Recognizing these sex differences is not only relevant to the treatment of cisgender individuals; it also holds paramount importance in addressing the healthcare needs of transgender patients, a demographic that is increasingly prominent in contemporary society. In response to these challenges, various agencies, including the National Institutes of Health, have actively directed their efforts toward advancing our comprehension of this phenomenon. Epigenetics has proven to play a crucial role in understanding sex differences in both healthy and disease states within the heart. This review presents a comprehensive overview of the physiological distinctions between males and females during the development of various cardiac pathologies, specifically focusing on unraveling the genetic and epigenetic mechanisms at play. Current findings related to distinct sex-chromosome compositions, the emergence of gender-biased genetic variations, and variations in hormonal profiles between sexes are highlighted. Additionally, the roles of DNA methylation, histone marks, and chromatin structure in mediating pathological sex differences are explored. To inspire further investigation into this crucial subject, we have conducted global analyses of various epigenetic features, leveraging data previously generated by the ENCODE project.

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.

Genome-wide association studies highlight novel risk loci for septal defects and left-sided congenital heart defects

BACKGROUND

Congenital heart defects (CHD) are structural defects of the heart affecting approximately 1% of newborns. They exhibit low penetrance and non-Mendelian patterns of inheritance as varied and complex traits. While genetic factors are known to play an important role in the development of CHD, the specific genetics remain unknown for the majority of patients. To elucidate the underlying genetic risk, we performed a genome wide association study (GWAS) of CHDs in general and specific CHD subgroups using the FinnGen Release 10 (R10) (N > 393,000), followed by functional fine-mapping through eQTL and co-localization analyses using the GTEx database.

RESULTS

We discovered three genome-wide significant loci associated with general CHD. Two of them were located in chromosome 17: 17q21.32 (rs2316327, intronic: LRRC37A2, Odds ratio (OR) [95% Confidence Interval (CI)] = 1.17[1.12-1.23], p = 1.5 × 10-9) and 17q25.3 (rs1293973611, nearest: BAHCC1, OR[95%CI] = 4.48[2.80-7.17], p = 7.0 × 10-10), respectively, and in addition to general CHD, the rs1293973611 locus was associated with the septal defect subtype. The third locus was in band 1p21.2 (rs35046143, nearest: PALMD, OR[95%CI] = 1.15[1.09-1.21], p = 7.1 × 10-9), and it was associated with general CHD and left-sided lesions. In the subgroup analysis, two additional loci were associated with septal defects (rs75230966 and rs6824295), and one with left-sided lesions (rs1305393195). In the eQTL analysis the variants rs2316327 (general CHD), and rs75230966 (septal defects) both located in 17q21.32 (with a LD r2 of 0.41) were both predicted to significantly associate with the expression of WNT9B in the atrial appendage tissue category. This effect was further confirmed by co-localization analysis, which also implicated WNT3 expression in the atrial appendage. A meta-analysis of general CHD together with the UK Biobank (combined N = 881,678) provided a different genome-wide significant locus in LRRC37A2; rs16941382 (OR[95%CI] = 1.15[1.11-1.20], p = 1.5 × 10-9) which is in significant LD with rs2316327.

CONCLUSIONS

Our results of general CHD and different CHD subcategories identified a complex risk locus on chromosome 17 near BAHCC1 and LRRC37A2, interacting with the genes WNT9B, WNT3 and MYL4, may constitute potential novel CHD risk associated loci, warranting future experimental tests to determine their role.

Spatially organized cellular communities form the developing human heart

The heart, which is the first organ to develop, is highly dependent on its form to function1,2. However, how diverse cardiac cell types spatially coordinate to create the complex morphological structures that are crucial for heart function remains unclear. Here we integrated single-cell RNA-sequencing with high-resolution multiplexed error-robust fluorescence in situ hybridization to resolve the identity of the cardiac cell types that develop the human heart. This approach also provided a spatial mapping of individual cells that enables illumination of their organization into cellular communities that form distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, which support their specialization according to the cellular ecosystem and anatomical region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities using in vivo conditional genetic mouse models and in vitro human pluripotent stem cell systems revealed multicellular signalling pathways that orchestrate the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. These detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodelling the human heart offer new insights into structural heart diseases and the engineering of complex multicellular tissues for human heart repair.

Atrial fibrillation: pathophysiology, genetic and epigenetic mechanisms

Atrial fibrillation (AF) is the most common supraventricular arrhythmia affecting up to 1% of the general population. Its prevalence dramatically increases with age and could reach up to ∼10% in the elderly. The management of AF is a complex issue that is object of extensive ongoing basic and clinical research, it depends on its genetic and epigenetic causes, and it varies considerably geographically and also according to the ethnicity. Mechanistically, over the last decade, Genome Wide Association Studies have uncovered over 100 genetic loci associated with AF, and have shown that European ancestry is associated with elevated risk of AF. These AF-associated loci revolve around different types of disturbances, including inflammation, electrical abnormalities, and structural remodeling. Moreover, the discovery of epigenetic regulatory mechanisms, involving non-coding RNAs, DNA methylation and histone modification, has allowed unravelling what modifications reshape the processes leading to arrhythmias. Our review provides a current state of the field regarding the identification and functional characterization of AF-related genetic and epigenetic regulatory networks, including ethnic differences. We discuss clear and emerging connections between genetic regulation and pathophysiological mechanisms of AF.

Age-related structural and functional changes of the intracardiac nervous system

BACKGROUND

Although aging is known to be associated with an increased incidence of both atrial and ventricular arrhythmias, there is limited knowledge about how Schwann cells (SC) and the intracardiac nervous system (iCNS) remodel with age. Here we investigate the differences in cardiac SC, parasympathetic nerve fibers, and muscarinic acetylcholine receptor M2 (M2R) expression in young and old mice. Additionally, we examine age-related changes in cardiac responses to sympathomimetic and parasympathomimetic drugs.

METHODS AND RESULTS

Lower SC density, lower SC proliferation and fewer parasympathetic nerve fibers were observed in cardiac and, as a control sciatic nerves from old (20-24 months) compared to young mice (2-3 months). In old mice, chondroitin sulfate proteoglycan 4 (CSPG4) was increased in sciatic but not cardiac nerves. Expression of M2R was lower in ventricular myocardium and ventricular conduction system from old mice compared to young mice, while no significant difference was seen in M2R expression in sino-atrial or atrio-ventricular node pacemaker tissue. Heart rate was slower and PQ intervals were longer in Langendorff-perfused hearts from old mice. Ventricular tachycardia and fibrillation were more frequently observed in response to carbachol administration in hearts from old mice versus those from young mice.

CONCLUSIONS

On the background of reduced presence of SC and parasympathetic nerve fibers, and of lower M2R expression in ventricular cardiomyocytes and conduction system of aged hearts, the propensity of ventricular arrhythmogenesis upon parasympathomimetic drug application is increased. Whether this is caused by an increase in heterogeneity of iCNS structure and function remains to be elucidated.

Epigenetics

Epigenetics is the study of heritable changes to the genome and gene expression patterns that are not caused by direct changes to the DNA sequence. Examples of these changes include posttranslational modifications to DNA-bound histone proteins, DNA methylation, and remodeling of nuclear architecture. Collectively, epigenetic changes provide a layer of regulation that affects transcriptional activity of genes while leaving DNA sequences unaltered. Sequence variants or mutations affecting enzymes responsible for modifying or sensing epigenetic marks have been identified in patients with congenital heart disease (CHD), and small-molecule inhibitors of epigenetic complexes have shown promise as therapies for adult heart diseases. Additionally, transgenic mice harboring mutations or deletions of genes encoding epigenetic enzymes recapitulate aspects of human cardiac disease. Taken together, these findings suggest that the evolving field of epigenetics will inform our understanding of congenital and adult cardiac disease and offer new therapeutic opportunities.

The benign nature and rare occurrence of cardiac myxoma as a possible consequence of the limited cardiac proliferative/ regenerative potential: a systematic review

BACKGROUND

Cardiac Myxoma is a primary tumor of heart. Its origins, rarity of the occurrence of primary cardiac tumors and how it may be related to limited cardiac regenerative potential, are not yet entirely known. This study investigates the key cardiac genes/ transcription factors (TFs) and signaling pathways to understand these important questions.

METHODS

Databases including PubMed, MEDLINE, and Google Scholar were searched for published articles without any date restrictions, involving cardiac myxoma, cardiac genes/TFs/signaling pathways and their roles in cardiogenesis, proliferation, differentiation, key interactions and tumorigenesis, with focus on cardiomyocytes.

RESULTS

The cardiac genetic landscape is governed by a very tight control between proliferation and differentiation-related genes/TFs/pathways. Cardiac myxoma originates possibly as a consequence of dysregulations in the gene expression of differentiation regulators including Tbx5, GATA4, HAND1/2, MYOCD, HOPX, BMPs. Such dysregulations switch the expression of cardiomyocytes into progenitor-like state in cardiac myxoma development by dysregulating Isl1, Baf60 complex, Wnt, FGF, Notch, Mef2c and others. The Nkx2-5 and MSX2 contribute predominantly to both proliferation and differentiation of Cardiac Progenitor Cells (CPCs), may possibly serve roles based on the microenvironment and the direction of cell circuitry in cardiac tumorigenesis. The Nkx2-5 in cardiac myxoma may serve to limit progression of tumorigenesis as it has massive control over the proliferation of CPCs. The cardiac cell type-specific genetic programming plays governing role in controlling the tumorigenesis and regenerative potential.

CONCLUSION

The cardiomyocytes have very limited proliferative and regenerative potential. They survive for long periods of time and tightly maintain the gene expression of differentiation genes such as Tbx5, GATA4 that interact with tumor suppressors (TS) and exert TS like effect. The total effect such gene expression exerts is responsible for the rare occurrence and benign nature of primary cardiac tumors. This prevents the progression of tumorigenesis. But this also limits the regenerative and proliferative potential of cardiomyocytes. Cardiac Myxoma develops as a consequence of dysregulations in these key genes which revert the cells towards progenitor-like state, hallmark of CM. The CM development in carney complex also signifies the role of TS in cardiac cells.

Precise identification of cell states altered in disease using healthy single-cell references

Joint analysis of single-cell genomics data from diseased tissues and a healthy reference can reveal altered cell states. We investigate whether integrated collections of data from healthy individuals (cell atlases) are suitable references for disease-state identification and whether matched control samples are needed to minimize false discoveries. We demonstrate that using a reference atlas for latent space learning followed by differential analysis against matched controls leads to improved identification of disease-associated cells, especially with multiple perturbed cell types. Additionally, when an atlas is available, reducing control sample numbers does not increase false discovery rates. Jointly analyzing data from a COVID-19 cohort and a blood cell atlas, we improve detection of infection-related cell states linked to distinct clinical severities. Similarly, we studied disease states in pulmonary fibrosis using a healthy lung atlas, characterizing two distinct aberrant basal states. Our analysis provides guidelines for designing disease cohort studies and optimizing cell atlas use.

Cis-Cardio: A comprehensive analysis platform for cardiovascular-relavant cis-regulation in human and mouse

Cis-regulatory elements are important molecular switches in controlling gene expression and are regarded as determinant hubs in the transcriptional regulatory network. Collection and processing of large-scale cis-regulatory data are urgent to decipher the potential mechanisms of cardiovascular diseases from a cis-regulatory element aspect. Here, we developed a novel web server, Cis-Cardio, which aims to document a large number of available cardiovascular-related cis-regulatory data and to provide analysis for unveiling the comprehensive mechanisms at a cis-regulation level. The current version of Cis-Cardio catalogs a total of 45,382,361 genomic regions from 1,013 human and mouse epigenetic datasets, including ATAC-seq, DNase-seq, Histone ChIP-seq, TF/TcoF ChIP-seq, RNA polymerase ChIP-seq, and Cohesin ChIP-seq. Importantly, Cis-Cardio provides six analysis tools, including region overlap analysis, element upstream/downstream analysis, transcription regulator enrichment analysis, variant interpretation, and protein-protein interaction-based co-regulatory analysis. Additionally, Cis-Cardio provides detailed and abundant (epi-) genetic annotations in cis-regulatory regions, such as super-enhancers, enhancers, transcription factor binding sites (TFBSs), methylation sites, common SNPs, risk SNPs, expression quantitative trait loci (eQTLs), motifs, DNase I hypersensitive sites (DHSs), and 3D chromatin interactions. In summary, Cis-Cardio is a valuable resource for elucidating and analyzing regulatory cues of cardiovascular-specific cis-regulatory elements. The platform is freely available at http://www.licpathway.net/Cis-Cardio/index.html.

Single-cell genomics improves the discovery of risk variants and genes of atrial fibrillation

Genome-wide association studies (GWAS) have linked hundreds of loci to cardiac diseases. However, in most loci the causal variants and their target genes remain unknown. We developed a combined experimental and analytical approach that integrates single cell epigenomics with GWAS to prioritize risk variants and genes. We profiled accessible chromatin in single cells obtained from human hearts and leveraged the data to study genetics of Atrial Fibrillation (AF), the most common cardiac arrhythmia. Enrichment analysis of AF risk variants using cell-type-resolved open chromatin regions (OCRs) implicated cardiomyocytes as the main mediator of AF risk. We then performed statistical fine-mapping, leveraging the information in OCRs, and identified putative causal variants in 122 AF-associated loci. Taking advantage of the fine-mapping results, our novel statistical procedure for gene discovery prioritized 46 high-confidence risk genes, highlighting transcription factors and signal transduction pathways important for heart development. In summary, our analysis provides a comprehensive map of AF risk variants and genes, and a general framework to integrate single-cell genomics with genetic studies of complex traits.

Loss of the Atrial Fibrillation-Related Gene, Zfhx3, Results in Atrial Dilation and Arrhythmias

BACKGROUND

ZFHX3 (zinc finger homeobox 3), a gene that encodes a large transcription factor, is at the second-most significantly associated locus with atrial fibrillation (AF), but its function in the heart is unknown. This study aims to identify causative genetic variation related to AF at the ZFHX3 locus and examine the impact of Zfhx3 loss on cardiac function in mice.

METHODS

CRISPR-Cas9 genome editing, chromatin immunoprecipitation, and luciferase assays in pluripotent stem cell-derived cardiomyocytes were used to identify causative genetic variation related to AF at the ZFHX3 locus. Cardiac function was assessed by echocardiography, magnetic resonance imaging, electrophysiology studies, calcium imaging, and RNA sequencing in mice with heterozygous and homozygous cardiomyocyte-restricted Zfhx3 loss (Zfhx3 Het and knockout, respectively). Human cardiac single-nucleus ATAC (assay for transposase-accessible chromatin)-sequencing data was analyzed to determine which genes in atrial cardiomyocytes are directly regulated by ZFHX3.

RESULTS

We found single-nucleotide polymorphism (SNP) rs12931021 modulates an enhancer regulating ZFHX3 expression, and the AF risk allele is associated with decreased ZFHX3 transcription. We observed a gene-dose response in AF susceptibility with Zfhx3 knockout mice having higher incidence, frequency, and burden of AF than Zfhx3 Het and wild-type mice, with alterations in conduction velocity, atrial action potential duration, calcium handling and the development of atrial enlargement and thrombus, and dilated cardiomyopathy. Zfhx3 loss results in atrial-specific differential effects on genes and signaling pathways involved in cardiac pathophysiology and AF.

CONCLUSIONS

Our findings implicate ZFHX3 as the causative gene at the 16q22 locus for AF, and cardiac abnormalities caused by loss of cardiac Zfhx3 are due to atrial-specific dysregulation of pathways involved in AF susceptibility. Together, these data reveal a novel and important role for Zfhx3 in the control of cardiac genes and signaling pathways essential for normal atrial function.

Uncovering Protein Networks in Cardiovascular Proteomics

Biological networks have been widely used in many different diseases to identify potential biomarkers and design drug targets. In the present review, we describe the main computational techniques for reconstructing and analyzing different types of protein networks and summarize the previous applications of such techniques in cardiovascular diseases. Existing tools are critically compared, discussing when each method is preferred such as the use of co-expression networks for functional annotation of protein clusters and the use of directed networks for inferring regulatory associations. Finally, we are presenting examples of reconstructing protein networks of different types (regulatory, co-expression, and protein-protein interaction networks). We demonstrate the necessity to reconstruct networks separately for each cardiovascular tissue type and disease entity and provide illustrative examples of the importance of taking into consideration relevant post-translational modifications. Finally, we demonstrate and discuss how the findings of protein networks could be interpreted using single-cell RNA-sequencing data.

Interpreting non-coding disease-associated human variants using single-cell epigenomics

Genome-wide association studies (GWAS) have linked hundreds of thousands of sequence variants in the human genome to common traits and diseases. However, translating this knowledge into a mechanistic understanding of disease-relevant biology remains challenging, largely because such variants are predominantly in non-protein-coding sequences that still lack functional annotation at cell-type resolution. Recent advances in single-cell epigenomics assays have enabled the generation of cell type-, subtype- and state-resolved maps of the epigenome in heterogeneous human tissues. These maps have facilitated cell type-specific annotation of candidate cis-regulatory elements and their gene targets in the human genome, enhancing our ability to interpret the genetic basis of common traits and diseases.

The heterocellular heart: identities, interactions, and implications for cardiology

The heterocellular nature of the heart has been receiving increasing attention in recent years. In addition to cardiomyocytes as the prototypical cell type of the heart, non-myocytes such as endothelial cells, fibroblasts, or immune cells are coming more into focus. The rise of single-cell sequencing technologies enables  identification of ever more subtle differences and has reignited the question of what defines a cell's identity. Here we provide an overview of the major cardiac cell types, describe their roles in homeostasis, and outline recent findings on non-canonical functions that may be of relevance for cardiology. We highlight modes of biochemical and biophysical interactions between different cardiac cell types and discuss the potential implications of the heterocellular nature of the heart for basic research and therapeutic interventions.

Exploring the genetic basis of coronary artery disease using functional genomics

Genome-wide Association Studies (GWAS) have identified more than 300 loci associated with coronary artery disease (CAD), defining the genetic risk map of the disease. However, the translation of the association signals into biological-pathophysiological mechanisms constitute a major challenge. Through a group of examples of studies focused on CAD, we discuss the rationale, basic principles and outcomes of the main methodologies implemented to prioritize and characterize causal variants and their target genes. Additionally, we highlight the strategies as well as the current methods that integrate association and functional genomics data to dissect the cellular specificity underlying the complexity of disease mechanisms. Despite the limitations of existing approaches, the increasing knowledge generated through functional studies helps interpret GWAS maps and opens novel avenues for the clinical usability of association data.

Enhlink infers distal and context-specific enhancer-promoter linkages

Enhancers play a crucial role in regulating gene expression and their functional status can be queried with cell type precision using using single-cell (sc)ATAC-seq. To facilitate analysis of such data, we developed Enhlink, a novel computational approach that leverages single-cell signals to infer linkages between regulatory DNA sequences, such as enhancers and promoters. Enhlink uses an ensemble strategy that integrates cell-level technical covariates to control for batch effects and biological covariates to infer robust condition-specific links and their associated p-values. It can integrate simultaneous gene expression and chromatin accessibility measurements of individual cells profiled by multi-omic experiments for increased specificity. We evaluated Enhlink using simulated and real scATAC-seq data, including those paired with physical enhancer-promoter links enumerated by promoter capture Hi-C and with multi-omic scATAC-/RNA-seq data we generated from the mouse striatum. These examples demonstrated that our method outperforms popular alternative strategies. In conjunction with eQTL analysis, Enhlink revealed a putative super-enhancer regulating key cell type-specific markers of striatal neurons. Taken together, our analyses demonstrate that Enhlink is accurate, powerful, and provides features that can lead to novel biological insights.

Three-dimensional genome rewiring in loci with human accelerated regions

Human accelerated regions (HARs) are conserved genomic loci that evolved at an accelerated rate in the human lineage and may underlie human-specific traits. We generated HARs and chimpanzee accelerated regions with an automated pipeline and an alignment of 241 mammalian genomes. Combining deep learning with chromatin capture experiments in human and chimpanzee neural progenitor cells, we discovered a significant enrichment of HARs in topologically associating domains containing human-specific genomic variants that change three-dimensional (3D) genome organization. Differential gene expression between humans and chimpanzees at these loci suggests rewiring of regulatory interactions between HARs and neurodevelopmental genes. Thus, comparative genomics together with models of 3D genome folding revealed enhancer hijacking as an explanation for the rapid evolution of HARs.

Multiomics Network Medicine Approaches to Precision Medicine and Therapeutics in Cardiovascular Diseases

Cardiovascular diseases (CVD) are the leading cause of death worldwide and display complex phenotypic heterogeneity caused by many convergent processes, including interactions between genetic variation and environmental factors. Despite the identification of a large number of associated genes and genetic loci, the precise mechanisms by which these genes systematically influence the phenotypic heterogeneity of CVD are not well understood. In addition to DNA sequence, understanding the molecular mechanisms of CVD requires data from other omics levels, including the epigenome, the transcriptome, the proteome, as well as the metabolome. Recent advances in multiomics technologies have opened new precision medicine opportunities beyond genomics that can guide precise diagnosis and personalized treatment. At the same time, network medicine has emerged as an interdisciplinary field that integrates systems biology and network science to focus on the interactions among biological components in health and disease, providing an unbiased framework through which to integrate systematically these multiomics data. In this review, we briefly present such multiomics technologies, including bulk omics and single-cell omics technologies, and discuss how they can contribute to precision medicine. We then highlight network medicine-based integration of multiomics data for precision medicine and therapeutics in CVD. We also include a discussion of current challenges, potential limitations, and future directions in the study of CVD using multiomics network medicine approaches.

Single-cell transcriptomes in the heart: when every epigenome counts

The response of an organ to stimuli emerges from the actions of individual cells. Recent cardiac single-cell RNA-sequencing studies of development, injury, and reprogramming have uncovered heterogeneous populations even among previously well-defined cell types, raising questions about what level of experimental resolution corresponds to disease-relevant, tissue-level phenotypes. In this review, we explore the biological meaning behind this cellular heterogeneity by undertaking an exhaustive analysis of single-cell transcriptomics in the heart (including a comprehensive, annotated compendium of studies published to date) and evaluating new models for the cardiac function that have emerged from these studies (including discussion and schematics that depict new hypotheses in the field). We evaluate the evidence to support the biological actions of newly identified cell populations and debate questions related to the role of cell-to-cell variability in development and disease. Finally, we present emerging epigenomic approaches that, when combined with single-cell RNA-sequencing, can resolve basic mechanisms of gene regulation and variability in cell phenotype.

In Vivo Dissection of Chamber-Selective Enhancers Reveals Estrogen-Related Receptor as a Regulator of Ventricular Cardiomyocyte Identity

BACKGROUND

Cardiac chamber-selective transcriptional programs underpin the structural and functional differences between atrial and ventricular cardiomyocytes (aCMs and vCMs). The mechanisms responsible for these chamber-selective transcriptional programs remain largely undefined.

METHODS

We nominated candidate chamber-selective enhancers (CSEs) by determining the genome-wide occupancy of 7 key cardiac transcription factors (GATA4, MEF2A, MEF2C, NKX2-5, SRF, TBX5, TEAD1) and transcriptional coactivator P300 in atria and ventricles. Candidate enhancers were tested using an adeno-associated virus-mediated massively parallel reporter assay. Chromatin features of CSEs were evaluated by performing assay of transposase accessible chromatin sequencing and acetylation of histone H3 at lysine 27-HiChIP on aCMs and vCMs. CSE sequence requirements were determined by systematic tiling mutagenesis of 29 CSEs at 5 bp resolution. Estrogen-related receptor (ERR) function in cardiomyocytes was evaluated by Cre-loxP-mediated inactivation of ERRα and ERRγ in cardiomyocytes.

RESULTS

We identified 134 066 and 97 506 regions reproducibly occupied by at least 1 transcription factor or P300, in atria or ventricles, respectively. Enhancer activities of 2639 regions bound by transcription factors or P300 were tested in aCMs and vCMs by adeno-associated virus-mediated massively parallel reporter assay. This identified 1092 active enhancers in aCMs or vCMs. Several overlapped loci associated with cardiovascular disease through genome-wide association studies, and 229 exhibited chamber-selective activity in aCMs or vCMs. Many CSEs exhibited differential chromatin accessibility between aCMs and vCMs, and CSEs were enriched for aCM- or vCM-selective acetylation of histone H3 at lysine 27-anchored loops. Tiling mutagenesis of 29 CSEs identified the binding motif of ERRα/γ as important for ventricular enhancer activity. The requirement of ERRα/γ to activate ventricular CSEs and promote vCM identity was confirmed by loss of the vCM gene profile in ERRα/γ knockout vCMs.

CONCLUSIONS

We identified 229 CSEs that could be useful research tools or direct therapeutic gene expression. We showed that chamber-selective multi-transcription factor, P300 occupancy, open chromatin, and chromatin looping are predictive features of CSEs. We found that ERRα/γ are essential for maintenance of ventricular identity. Finally, our gene expression, epigenetic, 3-dimensional genome, and enhancer activity atlas provide key resources for future studies of chamber-selective gene regulation.

Genetic architecture of spatial electrical biomarkers for cardiac arrhythmia and relationship with cardiovascular disease

The 3-dimensional spatial and 2-dimensional frontal QRS-T angles are measures derived from the vectorcardiogram. They are independent risk predictors for arrhythmia, but the underlying biology is unknown. Using multi-ancestry genome-wide association studies we identify 61 (58 previously unreported) loci for the spatial QRS-T angle (N = 118,780) and 11 for the frontal QRS-T angle (N = 159,715). Seven out of the 61 spatial QRS-T angle loci have not been reported for other electrocardiographic measures. Enrichments are observed in pathways related to cardiac and vascular development, muscle contraction, and hypertrophy. Pairwise genome-wide association studies with classical ECG traits identify shared genetic influences with PR interval and QRS duration. Phenome-wide scanning indicate associations with atrial fibrillation, atrioventricular block and arterial embolism and genetically determined QRS-T angle measures are associated with fascicular and bundle branch block (and also atrioventricular block for the frontal QRS-T angle). We identify potential biology involved in the QRS-T angle and their genetic relationships with cardiovascular traits and diseases, may inform future research and risk prediction.

Effect of tissue-grouped regulatory variants associated to type 2 diabetes in related secondary outcomes

Genome-wide association studies have identified over five hundred loci that contribute to variation in type 2 diabetes (T2D), an established risk factor for many diseases. However, the mechanisms and extent through which these loci contribute to subsequent outcomes remain elusive. We hypothesized that combinations of T2D-associated variants acting on tissue-specific regulatory elements might account for greater risk for tissue-specific outcomes, leading to diversity in T2D disease progression. We searched for T2D-associated variants acting on regulatory elements and expression quantitative trait loci (eQTLs) in nine tissues. We used T2D tissue-grouped variant sets as genetic instruments to conduct 2-Sample Mendelian Randomization (MR) in ten related outcomes whose risk is increased by T2D using the FinnGen cohort. We performed PheWAS analysis to investigate whether the T2D tissue-grouped variant sets had specific predicted disease signatures. We identified an average of 176 variants acting in nine tissues implicated in T2D, and an average of 30 variants acting on regulatory elements that are unique to the nine tissues of interest. In 2-Sample MR analyses, all subsets of regulatory variants acting in different tissues were associated with increased risk of the ten secondary outcomes studied on similar levels. No tissue-grouped variant set was associated with an outcome significantly more than other tissue-grouped variant sets. We did not identify different disease progression profiles based on tissue-specific regulatory and transcriptome information. Bigger sample sizes and other layers of regulatory information in critical tissues may help identify subsets of T2D variants that are implicated in certain secondary outcomes, uncovering system-specific disease progression.

Fine mapping spatiotemporal mechanisms of genetic variants underlying cardiac traits and disease

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

The adapted Activity-By-Contact model for enhancer-gene assignment and its application to single-cell data

MOTIVATION

Identifying regulatory regions in the genome is of great interest for understanding the epigenomic landscape in cells. One fundamental challenge in this context is to find the target genes whose expression is affected by the regulatory regions. A recent successful method is the Activity-By-Contact (ABC) model which scores enhancer-gene interactions based on enhancer activity and the contact frequency of an enhancer to its target gene. However, it describes regulatory interactions entirely from a gene's perspective, and does not account for all the candidate target genes of an enhancer. In addition, the ABC model requires two types of assays to measure enhancer activity, which limits the applicability. Moreover, there is neither implementation available that could allow for an integration with transcription factor (TF) binding information nor an efficient analysis of single-cell data.

RESULTS

We demonstrate that the ABC score can yield a higher accuracy by adapting the enhancer activity according to the number of contacts the enhancer has to its candidate target genes and also by considering all annotated transcription start sites of a gene. Further, we show that the model is comparably accurate with only one assay to measure enhancer activity. We combined our generalized ABC model with TF binding information and illustrated an analysis of a single-cell ATAC-seq dataset of the human heart, where we were able to characterize cell type-specific regulatory interactions and predict gene expression based on TF affinities. All executed processing steps are incorporated into our new computational pipeline STARE.

AVAILABILITY AND IMPLEMENTATION

The software is available at https://github.com/schulzlab/STARE.

CONTACT

marcel.schulz@em.uni-frankfurt.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Lineage-specific regulatory changes in hypertrophic cardiomyopathy unraveled by single-nucleus RNA-seq and spatial transcriptomics

Hypertrophic cardiomyopathy (HCM) is the most common cardiac genetic disorder characterized by cardiomyocyte hypertrophy and cardiac fibrosis. Pathological cardiac remodeling in the myocardium of HCM patients may progress to heart failure. An in-depth elucidation of the lineage-specific changes in pathological cardiac remodeling of HCM is pivotal for the development of therapies to mitigate the progression. Here, we performed single-nucleus RNA-seq of the cardiac tissues from HCM patients or healthy donors and conducted spatial transcriptomic assays on tissue sections from patients. Unbiased clustering of 55,122 nuclei from HCM and healthy conditions revealed 9 cell lineages and 28 clusters. Lineage-specific changes in gene expression, subpopulation composition, and intercellular communication in HCM were discovered through comparative analyses. According to the results of pseudotime ordering, differential expression analysis, and differential regulatory network analysis, potential key genes during the transition towards a failing state of cardiomyocytes such as FGF12, IL31RA, and CREB5 were identified. Transcriptomic dynamics underlying cardiac fibroblast activation were also uncovered, and potential key genes involved in cardiac fibrosis were obtained such as AEBP1, RUNX1, MEOX1, LEF1, and NRXN3. Using the spatial transcriptomic data, spatial activity patterns of the candidate genes, pathways, and subpopulations were confirmed on patient tissue sections. Moreover, we showed experimental evidence that in vitro knockdown of AEBP1 could promote the activation of human cardiac fibroblasts, and overexpression of AEBP1 could attenuate the TGFβ-induced activation. Our study provided a comprehensive analysis of the lineage-specific regulatory changes in HCM, which laid the foundation for targeted drug development in HCM.

Characterizing cis-regulatory elements using single-cell epigenomics

Cell type-specific gene expression patterns and dynamics during development or in disease are controlled by cis-regulatory elements (CREs), such as promoters and enhancers. Distinct classes of CREs can be characterized by their epigenomic features, including DNA methylation, chromatin accessibility, combinations of histone modifications and conformation of local chromatin. Tremendous progress has been made in cataloguing CREs in the human genome using bulk transcriptomic and epigenomic methods. However, single-cell epigenomic and multi-omic technologies have the potential to provide deeper insight into cell type-specific gene regulatory programmes as well as into how they change during development, in response to environmental cues and through disease pathogenesis. Here, we highlight recent advances in single-cell epigenomic methods and analytical tools and discuss their readiness for human tissue profiling.

Single-cell epigenome analysis reveals age-associated decay of heterochromatin domains in excitatory neurons in the mouse brain

Loss of heterochromatin has been implicated as a cause of pre-mature aging and age-associated decline in organ functions in mammals; however, the specific cell types and gene loci affected by this type of epigenetic change have remained unclear. To address this knowledge gap, we probed chromatin accessibility at single-cell resolution in the brains, hearts, skeletal muscles, and bone marrows from young, middle-aged, and old mice, and assessed age-associated changes at 353,126 candidate cis-regulatory elements (cCREs) across 32 major cell types. Unexpectedly, we detected increased chromatin accessibility within specific heterochromatin domains in old mouse excitatory neurons. The gain of chromatin accessibility at these genomic loci was accompanied by the cell-type-specific loss of heterochromatin and activation of LINE1 elements. Immunostaining further confirmed the loss of the heterochromatin mark H3K9me3 in the excitatory neurons but not in inhibitory neurons or glial cells. Our results reveal the cell-type-specific changes in chromatin landscapes in old mice and shed light on the scope of heterochromatin loss in mammalian aging.

Human-gained heart enhancers are associated with species-specific cardiac attributes

The heart, a vital organ which is first to develop, has adapted its size, structure and function in order to accommodate the circulatory demands for a broad range of animals. Although heart development is controlled by a relatively conserved network of transcriptional/chromatin regulators, how the human heart has evolved species-specific features to maintain adequate cardiac output and function remains to be defined. Here, we show through comparative epigenomic analysis the identification of enhancers and promoters that have gained activity in humans during cardiogenesis. These cis-regulatory elements (CREs) are associated with genes involved in heart development and function, and may account for species-specific differences between human and mouse hearts. Supporting these findings, genetic variants that are associated with human cardiac phenotypic/disease traits, particularly those differing between human and mouse, are enriched in human-gained CREs. During early stages of human cardiogenesis, these CREs are also gained within genomic loci of transcriptional regulators, potentially expanding their role in human heart development. In particular, we discovered that gained enhancers in the locus of the early human developmental regulator ZIC3 are selectively accessible within a subpopulation of mesoderm cells which exhibits cardiogenic potential, thus possibly extending the function of ZIC3 beyond its conserved left-right asymmetry role. Genetic deletion of these enhancers identified a human gained enhancer that was required for not only ZIC3 and early cardiac gene expression at the mesoderm stage but also cardiomyocyte differentiation. Overall, our results illuminate how human gained CREs may contribute to human-specific cardiac attributes, and provide insight into how transcriptional regulators may gain cardiac developmental roles through the evolutionary acquisition of enhancers.

Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Advances in Single-Cell Multi-Omics and Application in Cardiovascular Research

With the development of ever more powerful and versatile high-throughput sequencing techniques and innovative ways to capture single cells, mapping the multicellular tissues at the single-cell level is becoming routine practice. However, it is still challenging to depict the epigenetic landscape of a single cell, especially the genome-wide chromatin accessibility, histone modifications, and DNA methylation. We summarize the most recent methodologies to profile these epigenetic marks at the single-cell level. We also discuss the development and advancement of several multi-omics sequencing technologies from individual cells. Advantages and limitations of various methods to compare and integrate datasets obtained from different sources are also included with specific practical notes. Understanding the heart tissue at single-cell resolution and multi-modal levels will help to elucidate the cell types and states involved in physiological and pathological events during heart development and disease. The rich information produced from single-cell multi-omics studies will also promote the research of heart regeneration and precision medicine on heart diseases.

Single-nucleus chromatin accessibility profiling highlights regulatory mechanisms of coronary artery disease risk

Coronary artery disease (CAD) is a complex inflammatory disease involving genetic influences across cell types. Genome-wide association studies (GWAS) have identified over 200 loci associated with CAD, where the majority of risk variants reside in noncoding DNA sequences impacting cis-regulatory elements (CREs). Here, we applied single-nucleus ATAC-seq to profile 28,316 nuclei across coronary artery segments from 41 patients with varying stages of CAD, which revealed 14 distinct cellular clusters. We mapped ~320,000 accessible sites across all cells, identified cell type-specific elements, transcription factors, and prioritized functional CAD risk variants. . We identified elements in smooth muscle cell (SMC) transition states (e.g. fibromyocytes) and functional variants predicted to alter SMC and macrophage-specific regulation of MRAS (3q22) and LIPA (10q23), respectively. We further nominated key driver transcription factors such as PRDM16 and TBX2. Together, this single nucleus atlas provides a critical step towards interpreting regulatory mechanisms across the continuum of CAD risk.

Genomic enhancers in cardiac development and disease

The Human Genome Project marked a major milestone in the scientific community as it unravelled the ~3 billion bases that are central to crucial aspects of human life. Despite this achievement, it only scratched the surface of understanding how each nucleotide matters, both individually and as part of a larger unit. Beyond the coding genome, which comprises only ~2% of the whole genome, scientists have realized that large portions of the genome, not known to code for any protein, were crucial for regulating the coding genes. These large portions of the genome comprise the 'non-coding genome'. The history of gene regulation mediated by proteins that bind to the regulatory non-coding genome dates back many decades to the 1960s. However, the original definition of 'enhancers' was first used in the early 1980s. In this Review, we summarize benchmark studies that have mapped the role of cardiac enhancers in disease and development. We highlight instances in which enhancer-localized genetic variants explain the missing link to cardiac pathogenesis. Finally, we inspire readers to consider the next phase of exploring enhancer-based gene therapy for cardiovascular disease.

A single-cell atlas of chromatin accessibility in the human genome

Current catalogs of regulatory sequences in the human genome are still incomplete and lack cell type resolution. To profile the activity of gene regulatory elements in diverse cell types and tissues in the human body, we applied single-cell chromatin accessibility assays to 30 adult human tissue types from multiple donors. We integrated these datasets with previous single-cell chromatin accessibility data from 15 fetal tissue types to reveal the status of open chromatin for ∼1.2 million candidate cis-regulatory elements (cCREs) in 222 distinct cell types comprised of >1.3 million nuclei. We used these chromatin accessibility maps to delineate cell-type-specificity of fetal and adult human cCREs and to systematically interpret the noncoding variants associated with complex human traits and diseases. This rich resource provides a foundation for the analysis of gene regulatory programs in human cell types across tissues, life stages, and organ systems.

Improved epicardial cardiac fibroblast generation from iPSCs

Since the initial isolation of human embryonic stem cells and subsequent discovery of reprogramming methods for somatic cells, thousands of protocols have been developed to create each of the hundreds of cell types found in-vivo with significant focus on disease-prone systems, e.g., cardiovascular. Robust protocols exist for many of these cell types, except for cardiac fibroblasts (CF). Very recently, several competing methods have been developed to generate these cells through a developmentally conserved epicardial pathway. Such methods generate epicardial cells, but here we report that prolonged exposure to growth factors such as bFGF induces fibroblast spindle-like morphology and similar chromatin architecture to primary CFs. Media conditions for growth and assays are provided, as well as suggestions for seeding densities and timepoints for protein harvest of extracellular matrix. We demonstrate marker expression and matrix competency of resultant cells as shown next to primary human cardiac fibroblasts. These methods provide additional guidance to the original protocol and result in an increasingly stable phenotype.

ExpressHeart: Web Portal to Visualize Transcriptome Profiles of Non-Cardiomyocyte Cells

Unveiling the molecular features in the heart is essential for the study of heart diseases. Non-cardiomyocytes (nonCMs) play critical roles in providing structural and mechanical support to the working myocardium. There is an increasing amount of single-cell RNA-sequencing (scRNA-seq) data characterizing the transcriptomic profiles of nonCM cells. However, no tool allows researchers to easily access the information. Thus, in this study, we develop an open-access web portal, ExpressHeart, to visualize scRNA-seq data of nonCMs from five laboratories encompassing three species. ExpressHeart enables comprehensive visualization of major cell types and subtypes in each study; visualizes gene expression in each cell type/subtype in various ways; and facilitates identifying cell-type-specific and species-specific marker genes. ExpressHeart also provides an interface to directly combine information across datasets, for example, generating lists of high confidence DEGs by taking the intersection across different datasets. Moreover, ExpressHeart performs comparisons across datasets. We show that some homolog genes (e.g., Mmp14 in mice and mmp14b in zebrafish) are expressed in different cell types between mice and zebrafish, suggesting different functions across species. We expect ExpressHeart to serve as a valuable portal for investigators, shedding light on the roles of genes on heart development in nonCM cells.

Variants in a cis-regulatory element of TBX1 in conotruncal heart defect patients impair GATA6-mediated transactivation

Background

TBX1 (T-box transcription factor 1) is a major candidate gene that likely contributes to the etiology of velo-cardio-facial syndrome/DiGeorge syndrome (VCFS/DGS). Although the haploinsufficiency of TBX1 in both mice and humans results in congenital cardiac malformations, little has been elucidated about its upstream regulation. We aimed to explore the transcriptional regulation and dysregulation of TBX1.

Methods

Different TBX1 promoter reporters were constructed. Luciferase assays and electrophoretic mobility shift assays (EMSAs) were used to identify a cis-regulatory element within the TBX1 promoter region and its trans-acting factor. The expression of proteins was identified by immunohistochemistry and immunofluorescence. Variants in the cis-regulatory element were screened in conotruncal defect (CTD) patients. In vitro functional assays were performed to show the effects of the variants found in CTD patients on the transactivation of TBX1.

Results

We identified a cis-regulatory element within intron 1 of TBX1 that was found to be responsive to GATA6 (GATA binding protein 6), a transcription factor crucial for cardiogenesis. The expression patterns of GATA6 and TBX1 overlapped in the pharyngeal arches of human embryos. Transfection experiments and EMSA indicated that GATA6 could activate the transcription of TBX1 by directly binding with its GATA cis-regulatory element in vitro. Furthermore, sequencing analyses of 195 sporadic CTD patients without the 22q11.2 deletion or duplication identified 3 variants (NC_000022.11:g.19756832C > G, NC_000022.11:g.19756845C > T, and NC_000022.11:g. 19756902G > T) in the non-coding cis-regulatory element of TBX1. Luciferase assays showed that all 3 variants led to reduced transcription of TBX1 when incubated with GATA6.

Conclusions

Our findings showed that TBX1 might be a direct transcriptional target of GATA6, and variants in the non-coding cis-regulatory element of TBX1 disrupted GATA6-mediated transactivation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-021-01981-4.