Transcriptional profile of the rat cardiovascular system at single cell resolution

Background

Despite the critical role of the cardiovascular system, our understanding of its cellular and transcriptional diversity remains limited. We therefore sought to characterize the cellular composition, phenotypes, molecular pathways, and communication networks between cell types at the tissue and sub-tissue level across the cardiovascular system of the healthy Wistar rat, an important model in preclinical cardiovascular research. We obtained high quality tissue samples under controlled conditions that reveal a level of cellular detail so far inaccessible in human studies.

Methods and Results

We performed single nucleus RNA-sequencing in 78 samples in 10 distinct regions including the four chambers of the heart, ventricular septum, sinoatrial node, atrioventricular node, aorta, pulmonary artery, and pulmonary veins (PV), which produced an aggregate map of 505,835 nuclei. We identified 26 distinct cell types and additional subtypes, including a number of rare cell types such as PV cardiomyocytes and non-myelinating Schwann cells (NMSCs), and unique groups of vascular smooth muscle cells (VSMCs), endothelial cells (ECs) and fibroblasts (FBs), which gave rise to a detailed cell type distribution across tissues. We demonstrated differences in the cellular composition across different cardiac regions and tissue-specific differences in transcription for each cell type, highlighting the molecular diversity and complex tissue architecture of the cardiovascular system. Specifically, we observed great transcriptional heterogeneities among ECs and FBs. Importantly, several cell subtypes had a unique regional localization such as a subtype of VSMCs enriched in the large vasculature. We found the cellular makeup of PV tissue is closer to heart tissue than to the large arteries. We further explored the ligand-receptor repertoire across cell clusters and tissues, and observed tissue-enriched cellular communication networks, including heightened Nppa - Npr1 / 2 / 3 signaling in the sinoatrial node.

Conclusions

Through a large single nucleus sequencing effort encompassing over 500,000 nuclei, we broadened our understanding of cellular transcription in the healthy cardiovascular system. The existence of tissue-restricted cellular phenotypes suggests regional regulation of cardiovascular physiology. The overall conservation in gene expression and molecular pathways across rat and human cell types, together with our detailed transcriptional characterization of each cell type, offers the potential to identify novel therapeutic targets and improve preclinical models of cardiovascular disease.

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.

Cell-Specific Mechanisms in the Heart of COVID-19 Patients

From the onset of the pandemic, evidence of cardiac involvement in acute COVID-19 abounded. Cardiac presentations ranged from arrhythmias to ischemia, myopericarditis/myocarditis, ventricular dysfunction to acute heart failure, and even cardiogenic shock. Elevated serum cardiac troponin levels were prevalent among hospitalized patients with COVID-19; the higher the magnitude of troponin elevation, the greater the COVID-19 illness severity and in-hospital death risk. Whether these consequences were due to direct SARS-CoV-2 infection of cardiac cells or secondary to inflammatory responses steered early cardiac autopsy studies. SARS-CoV-2 was reportedly detected in endothelial cells, cardiac myocytes, and within the extracellular space. However, findings were inconsistent and different methodologies had their limitations. Initial autopsy reports suggested that SARS-CoV-2 myocarditis was common, setting off studies to find and phenotype inflammatory infiltrates in the heart. Nonetheless, subsequent studies rarely detected myocarditis. Microthrombi, cardiomyocyte necrosis, and inflammatory infiltrates without cardiomyocyte damage were much more common. In vitro and ex vivo experimental platforms have assessed the cellular tropism of SARS-CoV-2 and elucidated mechanisms of viral entry into and replication within cardiac cells. Data point to pericytes as the primary target of SARS-CoV-2 in the heart. Infection of pericytes can account for the observed pericyte and endothelial cell death, innate immune response, and immunothrombosis commonly observed in COVID-19 hearts. These processes are bidirectional and synergistic, rendering a definitive order of events elusive. Single-cell/nucleus analyses of COVID-19 myocardial tissue and isolated cardiac cells have provided granular data about the cellular composition and cell type-specific transcriptomic signatures of COVID-19 and microthrombi-positive COVID-19 hearts. Still, much remains unknown and more in vivo studies are needed. This review seeks to provide an overview of the current understanding of COVID-19 cardiac pathophysiology. Cell type-specific mechanisms and the studies that provided such insights will be highlighted. Given the unprecedented pace of COVID-19 research, more mechanistic details are sure to emerge since the writing of this review. Importantly, our current knowledge offers significant clues about the cardiac pathophysiology of long COVID-19, the increased postrecovery risk of cardiac events, and thus, the future landscape of cardiovascular disease.

How Does COVID-19 Affect the Heart?

PURPOSE OF REVIEW

Cardiac consequences occur in both acute COVID-19 and post-acute sequelae of COVID-19 (PASC). Here, we highlight the current understanding about COVID-19 cardiac effects, based upon clinical, imaging, autopsy, and molecular studies.

RECENT FINDINGS

COVID-19 cardiac effects are heterogeneous. Multiple, concurrent cardiac histopathologic findings have been detected on autopsies of COVID-19 non-survivors. Microthrombi and cardiomyocyte necrosis are commonly detected. Macrophages often infiltrate the heart at high density but without fulfilling histologic criteria for myocarditis. The high prevalences of microthrombi and inflammatory infiltrates in fatal COVID-19 raise the concern that recovered COVID-19 patients may have similar but subclinical cardiac pathology. Molecular studies suggest that SARS-CoV-2 infection of cardiac pericytes, dysregulated immunothrombosis, and pro-inflammatory and anti-fibrinolytic responses underlie COVID-19 cardiac pathology. The extent and nature by which mild COVID-19 affects the heart is unknown. Imaging and epidemiologic studies of recovered COVID-19 patients suggest that even mild illness confers increased risks of cardiac inflammation, cardiovascular disorders, and cardiovascular death. The mechanistic details of COVID-19 cardiac pathophysiology remain under active investigation. The ongoing evolution of SARS-CoV-2 variants and vast numbers of recovered COVID-19 patients portend a burgeoning global cardiovascular disease burden. Our ability to prevent and treat cardiovascular disease in the future will likely depend on comprehensive understanding of COVID-19 cardiac pathophysiologic phenotypes.

Sex-specific responses to slow progressive pressure overload in a large animal model of HFpEF

Approximately 50% of all heart failure (HF) diagnoses can be classified as HF with preserved ejection fraction (HFpEF). HFpEF is more prevalent in females compared with males, but the underlying mechanisms are unknown. We previously showed that pressure overload (PO) in male felines induces a cardiopulmonary phenotype with essential features of human HFpEF. The goal of this study was to determine if slow progressive PO induces distinct cardiopulmonary phenotypes in females and males in the absence of other pathological stressors. Female and male felines underwent aortic constriction (banding) or sham surgery after baseline echocardiography, pulmonary function testing, and blood sampling. These assessments were repeated at 2 and 4 mo postsurgery to document the effects of slow progressive pressure overload. At 4 mo, invasive hemodynamic studies were also performed. Left ventricle (LV) tissue was collected for histology, myofibril mechanics, extracellular matrix (ECM) mass spectrometry, and single-nucleus RNA sequencing (snRNAseq). The induced pressure overload (PO) was not different between sexes. PO also induced comparable changes in LV wall thickness and myocyte cross-sectional area in both sexes. Both sexes had preserved ejection fraction, but males had a slightly more robust phenotype in hemodynamic and pulmonary parameters. There was no difference in LV fibrosis and ECM composition between banded male and female animals. LV snRNAseq revealed changes in gene programs of individual cell types unique to males and females after PO. Based on these results, both sexes develop cardiopulmonary dysfunction but the phenotype is somewhat less advanced in females.NEW & NOTEWORTHY We performed a comprehensive assessment to evaluate the effects of slow progressive pressure overload on cardiopulmonary function in a large animal model of heart failure with preserved ejection fraction (HFpEF) in males and females. Functional and structural assessments were performed at the organ, tissue, cellular, protein, and transcriptional levels. This is the first study to compare snRNAseq and ECM mass spectrometry of HFpEF myocardium from males and females. The results broaden our understanding of the pathophysiological response of both sexes to pressure overload. Both sexes developed a robust cardiopulmonary phenotype, but the phenotype was equal or a bit less robust in females.

A beary good genome: Haplotype-resolved, chromosome-level assembly of the brown bear (Ursus arctos)

The brown bear (Ursus arctos) is the second largest and most widespread extant terrestrial carnivore on Earth and has recently emerged as a medical model for human metabolic diseases. Here, we report a fully phased chromosome-level assembly of a male North American brown bear built by combining Pacific Biosciences (PacBio) HiFi data and publicly available Hi-C data. The final genome size is 2.47 Gigabases (Gb) with a scaffold and contig N50 length of 70.08 and 43.94 Megabases (Mb), respectively. Benchmarking Universal Single-Copy Ortholog (BUSCO) analysis revealed that 94.5% of single copy orthologs from Mammalia were present in the genome (the highest of any ursid genome to date). Repetitive elements accounted for 44.48% of the genome and a total of 20,480 protein coding genes were identified. Based on whole genome alignment to the polar bear, the brown bear is highly syntenic with the polar bear, and our phylogenetic analysis of 7,246 single-copy orthologs supports the currently proposed species tree for Ursidae. This highly contiguous genome assembly will support future research on both the evolutionary history of the bear family and the physiological mechanisms behind hibernation, the latter of which has broad medical implications.

Transcriptome variation in human tissues revealed by long-read sequencing

Regulation of transcript structure generates transcript diversity and plays an important role in human disease1-7. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure8-16. In this Article, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from Genotype-Tissue Expression (GTEx) tissues and cell lines, complementing the GTEx resource. We identified just over 70,000 novel transcripts for annotated genes, and validated the protein expression of 10% of novel transcripts. We developed a new computational package, LORALS, to analyse the genetic effects of rare and common variants on the transcriptome by allele-specific analysis of long reads. We characterized allele-specific expression and transcript structure events, providing new insights into the specific transcript alterations caused by common and rare genetic variants and highlighting the resolution gained from long-read data. We were able to perturb the transcript structure upon knockdown of PTBP1, an RNA binding protein that mediates splicing, thereby finding genetic regulatory effects that are modified by the cellular environment. Finally, we used this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.

Single-nucleus profiling of human dilated and hypertrophic cardiomyopathy

Heart failure encompasses a heterogeneous set of clinical features that converge on impaired cardiac contractile function1,2 and presents a growing public health concern. Previous work has highlighted changes in both transcription and protein expression in failing hearts3,4, but may overlook molecular changes in less prevalent cell types. Here we identify extensive molecular alterations in failing hearts at single-cell resolution by performing single-nucleus RNA sequencing of nearly 600,000 nuclei in left ventricle samples from 11 hearts with dilated cardiomyopathy and 15 hearts with hypertrophic cardiomyopathy as well as 16 non-failing hearts. The transcriptional profiles of dilated or hypertrophic cardiomyopathy hearts broadly converged at the tissue and cell-type level. Further, a subset of hearts from patients with cardiomyopathy harbour a unique population of activated fibroblasts that is almost entirely absent from non-failing samples. We performed a CRISPR-knockout screen in primary human cardiac fibroblasts to evaluate this fibrotic cell state transition; knockout of genes associated with fibroblast transition resulted in a reduction of myofibroblast cell-state transition upon TGFβ1 stimulation for a subset of genes. Our results provide insights into the transcriptional diversity of the human heart in health and disease as well as new potential therapeutic targets and biomarkers for heart failure.

Increased atrial effectiveness of flecainide conferred by altered biophysical properties of sodium channels

Atrial fibrillation (AF) affects over 1% of the population and is a leading cause of stroke and heart failure in the elderly. A feared side effect of sodium channel blocker therapy, ventricular pro-arrhythmia, appears to be relatively rare in patients with AF. The biophysical reasons for this relative safety of sodium blockers are not known. Our data demonstrates intrinsic differences between atrial and ventricular cardiac voltage-gated sodium currents (INa), leading to reduced maximum upstroke velocity of action potential and slower conduction, in left atria compared to ventricle. Reduced atrial INa is only detected at physiological membrane potentials and is driven by alterations in sodium channel biophysical properties and not by NaV1.5 protein expression. Flecainide displayed greater inhibition of atrial INa, greater reduction of maximum upstroke velocity of action potential, and slowed conduction in atrial cells and tissue. Our work highlights differences in biophysical properties of sodium channels in left atria and ventricles and their response to flecainide. These differences can explain the relative safety of sodium channel blocker therapy in patients with atrial fibrillation.

Deep learning enables genetic analysis of the human thoracic aorta

Enlargement or aneurysm of the aorta predisposes to dissection, an important cause of sudden death. We trained a deep learning model to evaluate the dimensions of the ascending and descending thoracic aorta in 4.6 million cardiac magnetic resonance images from the UK Biobank. We then conducted genome-wide association studies in 39,688 individuals, identifying 82 loci associated with ascending and 47 with descending thoracic aortic diameter, of which 14 loci overlapped. Transcriptome-wide analyses, rare-variant burden tests and human aortic single nucleus RNA sequencing prioritized genes including SVIL, which was strongly associated with descending aortic diameter. A polygenic score for ascending aortic diameter was associated with thoracic aortic aneurysm in 385,621 UK Biobank participants (hazard ratio = 1.43 per s.d., confidence interval 1.32-1.54, P = 3.3 × 10-20). Our results illustrate the potential for rapidly defining quantitative traits with deep learning, an approach that can be broadly applied to biomedical images.

Clinico-histopathologic and single nuclei RNA sequencing insights into cardiac injury and microthrombi in critical COVID-19

Acute cardiac injury is prevalent in critical COVID-19 and associated with increased mortality. Its etiology remains debated, as initially presumed causes--- myocarditis and cardiac necrosis--- have proven uncommon. To elucidate the pathophysiology of COVID-19-associated cardiac injury, we conducted a prospective study of the first 69 consecutive COVID-19 decedents at Columbia University Irving Medical Center in New York City. Of six acute cardiac histopathologic features, microthrombi was the most commonly detected amongst our cohort (n=48, 70%). We tested associations of cardiac microthrombi with biomarkers of inflammation, cardiac injury, and fibrinolysis and with in-hospital antiplatelet therapy, therapeutic anticoagulation, and corticosteroid treatment, while adjusting for multiple clinical factors, including COVID-19 therapies. Higher peak erythrocyte sedimentation rate and c-reactive protein were independently associated with increased odds of microthrombi, supporting an immunothrombotic etiology. Using single nuclei RNA-sequencing analysis on 3 patients with and 4 patients without cardiac microthrombi, we discovered an enrichment of pro-thrombotic/anti-fibrinolytic, extracellular matrix remodeling, and immune-potentiating signaling amongst cardiac fibroblasts in microthrombi-positive, relative to microthrombi-negative, COVID-19 hearts. Non-COVID-19 non-failing hearts were used as reference controls. Our study identifies a specific transcriptomic signature in cardiac fibroblasts as a salient feature of microthrombi-positive COVID-19 hearts. Our findings warrant further mechanistic study as cardiac fibroblasts may represent a potential therapeutic target for COVID-19-associated cardiac microthrombi.

SnRNA sequencing defines signaling by RBC-derived extracellular vesicles in the murine heart

In a unique model of fluorescent based mapping of EV recipient cells, RBC-EVs were found to signal to cardiac cells and regulate gene expression in a model of ischemic heart failure.

Extracellular vesicles (EVs) mediate intercellular signaling by transferring their cargo to recipient cells, but the functional consequences of signaling are not fully appreciated. RBC-derived EVs are abundant in circulation and have been implicated in regulating immune responses. Here, we use a transgenic mouse model for fluorescence-based mapping of RBC-EV recipient cells to assess the role of this intercellular signaling mechanism in heart disease. Using fluorescent-based mapping, we detected an increase in RBC-EV–targeted cardiomyocytes in a murine model of ischemic heart failure. Single cell nuclear RNA sequencing of the heart revealed a complex landscape of cardiac cells targeted by RBC-EVs, with enrichment of genes implicated in cell proliferation and stress signaling pathways compared with non-targeted cells. Correspondingly, cardiomyocytes targeted by RBC-EVs more frequently express cellular markers of DNA synthesis, suggesting the functional significance of EV-mediated signaling. In conclusion, our mouse model for mapping of EV-recipient cells reveals a complex cellular network of RBC-EV–mediated intercellular communication in ischemic heart failure and suggests a functional role for this mode of intercellular signaling.

Molecular Pathophysiology of Cardiac Injury and Cardiac Microthrombi in Fatal COVID-19: Insights from Clinico-histopathologic and Single Nuclei RNA Sequencing Analyses

Cardiac injury is associated with critical COVID-19, yet its etiology remains debated. To elucidate the pathogenic mechanisms of COVID-19-associated cardiac injury, we conducted a single-center prospective cohort study of 69 COVID-19 decedents. Of six cardiac histopathologic features, microthrombi was the most commonly detected (n=48, 70%). We tested associations of cardiac microthrombi with biomarkers of inflammation, cardiac injury, and fibrinolysis and with in-hospital antiplatelet therapy, therapeutic anticoagulation, and corticosteroid treatment, while adjusting for multiple clinical factors, including COVID-19 therapies. Higher peak ESR and CRP during hospitalization were independently associated with higher odds of microthrombi. Using single nuclei RNA-sequence analysis, we discovered an enrichment of pro-thrombotic/anti-fibrinolytic, extracellular matrix remodeling, and immune-potentiating signaling amongst cardiac fibroblasts in microthrombi-positive COVID-19 hearts relative to microthrombi-negative COVID-19. Non-COVID-19 non-failing hearts were used as reference controls. Our cumulative findings identify the specific transcriptomic changes in cardiac fibroblasts as salient features of COVID-19-associated cardiac microthrombi.

COVID-19 and Cardiovascular Disease: From Bench to Bedside

A pandemic of historic impact, coronavirus disease 2019 (COVID-19) has potential consequences on the cardiovascular health of millions of people who survive infection worldwide. Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2), the etiologic agent of COVID-19, can infect the heart, vascular tissues, and circulating cells through ACE2 (angiotensin-converting enzyme 2), the host cell receptor for the viral spike protein. Acute cardiac injury is a common extrapulmonary manifestation of COVID-19 with potential chronic consequences. This update provides a review of the clinical manifestations of cardiovascular involvement, potential direct SARS-CoV-2 and indirect immune response mechanisms impacting the cardiovascular system, and implications for the management of patients after recovery from acute COVID-19 infection.

Single-cell meta-analysis of SARS-CoV-2 entry genes across tissues and demographics

Angiotensin-converting enzyme 2 (ACE2) and accessory proteases (TMPRSS2 and CTSL) are needed for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cellular entry, and their expression may shed light on viral tropism and impact across the body. We assessed the cell-type-specific expression of ACE2, TMPRSS2 and CTSL across 107 single-cell RNA-sequencing studies from different tissues. ACE2, TMPRSS2 and CTSL are coexpressed in specific subsets of respiratory epithelial cells in the nasal passages, airways and alveoli, and in cells from other organs associated with coronavirus disease 2019 (COVID-19) transmission or pathology. We performed a meta-analysis of 31 lung single-cell RNA-sequencing studies with 1,320,896 cells from 377 nasal, airway and lung parenchyma samples from 228 individuals. This revealed cell-type-specific associations of age, sex and smoking with expression levels of ACE2, TMPRSS2 and CTSL. Expression of entry factors increased with age and in males, including in airway secretory cells and alveolar type 2 cells. Expression programs shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues included genes that may mediate viral entry, key immune functions and epithelial-macrophage cross-talk, such as genes involved in the interleukin-6, interleukin-1, tumor necrosis factor and complement pathways. Cell-type-specific expression patterns may contribute to the pathogenesis of COVID-19, and our work highlights putative molecular pathways for therapeutic intervention.

Epigenetic Analyses of Human Left Atrial Tissue Identifies Gene Networks Underlying Atrial Fibrillation

BACKGROUND

Atrial fibrillation (AF) often arises from structural abnormalities in the left atria (LA). Annotation of the noncoding genome in human LA is limited, as are effects on gene expression and chromatin architecture. Many AF-associated genetic variants reside in noncoding regions; this knowledge gap impairs efforts to understand the molecular mechanisms of AF and cardiac conduction phenotypes.

METHODS

We generated a model of the LA noncoding genome by profiling 7 histone post-translational modifications (active: H3K4me3, H3K4me2, H3K4me1, H3K27ac, H3K36me3; repressive: H3K27me3, H3K9me3), CTCF binding, and gene expression in samples from 5 individuals without structural heart disease or AF. We used MACS2 to identify peak regions (P<0.01), applied a Markov model to classify regulatory elements, and annotated this model with matched gene expression data. We intersected chromatin states with expression quantitative trait locus, DNA methylation, and HiC chromatin interaction data from LA and left ventricle. Finally, we integrated genome-wide association data for AF and electrocardiographic traits to link disease-related variants to genes.

RESULTS

Our model identified 21 epigenetic states, encompassing regulatory motifs, such as promoters, enhancers, and repressed regions. Genes were regulated by proximal chromatin states; repressive states were associated with a significant reduction in gene expression (P<2×10-16). Chromatin states were differentially methylated, promoters were less methylated than repressed regions (P<2×10-16). We identified over 15 000 LA-specific enhancers, defined by homeobox family motifs, and annotated several cardiovascular disease susceptibility loci. Intersecting AF and PR genome-wide association studies loci with long-range chromatin conformation data identified a gene interaction network dominated by NKX2-5, TBX3, ZFHX3, and SYNPO2L.

CONCLUSIONS

Profiling the noncoding genome provides new insights into the gene expression and chromatin regulation in human LA tissue. These findings enabled identification of a gene network underlying AF; our experimental and analytic approach can be extended to identify molecular mechanisms for other cardiac diseases and traits.

Myocyte-Specific Upregulation of ACE2 in Cardiovascular Disease: Implications for SARS-CoV-2-Mediated Myocarditis

Supplemental Digital Content is available in the text.

Transcriptional and Cellular Diversity of the Human Heart

BACKGROUND

The human heart requires a complex ensemble of specialized cell types to perform its essential function. A greater knowledge of the intricate cellular milieu of the heart is critical to increase our understanding of cardiac homeostasis and pathology. As recent advances in low-input RNA sequencing have allowed definitions of cellular transcriptomes at single-cell resolution at scale, we have applied these approaches to assess the cellular and transcriptional diversity of the nonfailing human heart.

METHODS

Microfluidic encapsulation and barcoding was used to perform single nuclear RNA sequencing with samples from 7 human donors, selected for their absence of overt cardiac disease. Individual nuclear transcriptomes were then clustered based on transcriptional profiles of highly variable genes. These clusters were used as the basis for between-chamber and between-sex differential gene expression analyses and intersection with genetic and pharmacologic data.

RESULTS

We sequenced the transcriptomes of 287 269 single cardiac nuclei, revealing 9 major cell types and 20 subclusters of cell types within the human heart. Cellular subclasses include 2 distinct groups of resident macrophages, 4 endothelial subtypes, and 2 fibroblast subsets. Comparisons of cellular transcriptomes by cardiac chamber or sex reveal diversity not only in cardiomyocyte transcriptional programs but also in subtypes involved in extracellular matrix remodeling and vascularization. Using genetic association data, we identified strong enrichment for the role of cell subtypes in cardiac traits and diseases. Intersection of our data set with genes on cardiac clinical testing panels and the druggable genome reveals striking patterns of cellular specificity.

CONCLUSIONS

Using large-scale single nuclei RNA sequencing, we defined the transcriptional and cellular diversity in the normal human heart. Our identification of discrete cell subtypes and differentially expressed genes within the heart will ultimately facilitate the development of new therapeutics for cardiovascular diseases.

Epigenetic and Transcriptional Networks Underlying Atrial Fibrillation

Genome-wide association studies have uncovered over a 100 genetic loci associated with atrial fibrillation (AF), the most common arrhythmia. Many of the top AF-associated loci harbor key cardiac transcription factors, including PITX2, TBX5, PRRX1, and ZFHX3. Moreover, the vast majority of the AF-associated variants lie within noncoding regions of the genome where causal variants affect gene expression by altering the activity of transcription factors and the epigenetic state of chromatin. In this review, we discuss a transcriptional regulatory network model for AF defined by effector genes in Genome-wide association studies loci. We describe the current state of the field regarding the identification and function of AF-relevant gene regulatory networks, including variant regulatory elements, dose-sensitive transcription factor functionality, target genes, and epigenetic states. We illustrate how altered transcriptional networks may impact cardiomyocyte function and ionic currents that impact AF risk. Last, we identify the need for improved tools to identify and functionally test transcriptional components to define the links between genetic variation, epigenetic gene regulation, and atrial function.

Multi-ancestry GWAS of the electrocardiographic PR interval identifies 202 loci underlying cardiac conduction

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N = 293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.

Myocyte Specific Upregulation of ACE2 in Cardiovascular Disease: Implications for SARS-CoV-2 mediated myocarditis

Coronavirus disease 2019 (COVID-19) is a global pandemic caused by a novel severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). SARS-CoV-2 infection of host cells occurs predominantly via binding of the viral surface spike protein to the human angiotensin-converting enzyme 2 (ACE2) receptor. Hypertension and pre-existing cardiovascular disease are risk factors for morbidity from COVID-19, and it remains uncertain whether the use of angiotensin converting enzyme inhibitors (ACEi) or angiotensin receptor blockers (ARB) impacts infection and disease. Here, we aim to shed light on this question by assessing ACE2 expression in normal and diseased human myocardial samples profiled by bulk and single nucleus RNA-seq.

Identification of Functional Variant Enhancers Associated With Atrial Fibrillation

RATIONALE

Genome-wide association studies have identified a large number of common variants (single-nucleotide polymorphisms) associated with atrial fibrillation (AF). These variants are located mainly in noncoding regions of the genome and likely include variants that modulate the function of transcriptional regulatory elements (REs) such as enhancers. However, the actual REs modulated by variants and the target genes of such REs remain to be identified. Thus, the biological mechanisms by which genetic variation promotes AF has thus far remained largely unexplored.

OBJECTIVE

To identify REs in genome-wide association study loci that are influenced by AF-associated variants.

METHODS AND RESULTS

We screened 2.45 Mbp of human genomic DNA containing 12 strongly AF-associated loci for RE activity using self-transcribing active regulatory region sequencing and a recently generated monoclonal line of conditionally immortalized rat atrial myocytes. We identified 444 potential REs, 55 of which contain AF-associated variants (P<10^-8). Subsequently, using an adaptation of the self-transcribing active regulatory region sequencing approach, we identified 24 variant REs with allele-specific regulatory activity. By mining available chromatin conformation data, the possible target genes of these REs were mapped. To define the physiological function and target genes of such REs, we deleted the orthologue of an RE containing noncoding variants in the Hcn4 (potassium/sodium hyperpolarization-activated cyclic nucleotide-gated channel 4) locus of the mouse genome. Mice heterozygous for the RE deletion showed bradycardia, sinus node dysfunction, and selective loss of Hcn4 expression.

CONCLUSIONS

We have identified REs at multiple genetic loci for AF and found that loss of an RE at the HCN4 locus results in sinus node dysfunction and reduced gene expression. Our approach can be broadly applied to facilitate the identification of human disease-relevant REs and target genes at cardiovascular genome-wide association studies loci.

Long-range Pitx2c enhancer-promoter interactions prevent predisposition to atrial fibrillation

Genome-wide association studies found that increased risk for atrial fibrillation (AF), the most common human heart arrhythmia, is associated with noncoding sequence variants located in proximity to PITX2 Cardiomyocyte-specific epigenomic and comparative genomics uncovered 2 AF-associated enhancers neighboring PITX2 with varying conservation in mice. Chromosome conformation capture experiments in mice revealed that the Pitx2c promoter directly contacted the AF-associated enhancer regions. CRISPR/Cas9-mediated deletion of a 20-kb topologically engaged enhancer led to reduced Pitx2c transcription and AF predisposition. Allele-specific chromatin immunoprecipitation sequencing on hybrid heterozygous enhancer knockout mice revealed that long-range interaction of an AF-associated region with the Pitx2c promoter was required for maintenance of the Pitx2c promoter chromatin state. Long-range looping was mediated by CCCTC-binding factor (CTCF), since genetic disruption of the intronic CTCF-binding site caused reduced Pitx2c expression, AF predisposition, and diminished active chromatin marks on Pitx2 AF risk variants located at 4q25 reside in genomic regions possessing long-range transcriptional regulatory functions directed at PITX2.

Genome-Wide Association Study-Driven Gene-Set Analyses, Genetic, and Functional Follow-Up Suggest GLIS1 as a Susceptibility Gene for Mitral Valve Prolapse

Background Mitral valve prolapse (MVP) is a common heart valve disease, the most frequent indication for valve repair or replacement. MVP is characterized by excess extracellular matrix secretion and cellular disorganization, which leads to bulky valves that are unable to coapt correctly during ventricular systole resulting in mitral regurgitation, and it is associated with sudden cardiac death. Here we aim to characterize globally the biological mechanisms underlying genetic susceptibility to MVP to better characterize its triggering mechanisms. Methods We applied i-GSEA4GWAS and DEPICT, two pathway enrichment tools to MVP genome-wide association studies. We followed-up the association with MVP in an independent dataset of cases and controls. This research was conducted using the UK Biobank Resource. Immunohistochemistry staining for Glis1 (GLIS family zinc finger 1) was conducted in developing heart of mice. Knockdown of Glis1 using morpholinos was performed in zebrafish animals 72 hours postfertilization. Results We show that genes at risk loci are involved in biological functions relevant to actin filament organization, cytoskeleton biology, and cardiac development. The enrichment for positive regulation of transcription, cell proliferation, and migration motivated the follow-up of GLIS1, a transcription factor from the Krüppel-like zinc finger family. In combination with previously available data, we now report a genome-wide significant association with MVP (odds ratio, 1.20; P=4.36×10-10), indicating that Glis1 is expressed during embryonic development predominantly in nuclei of endothelial and interstitial cells of mitral valves in mouse. We also show that Glis1 knockdown causes atrioventricular regurgitation in developing hearts in zebrafish. Conclusions Our findings define globally molecular and cellular mechanisms underlying common genetic susceptibility to MVP and implicate established and unprecedented mechanisms. Through the GLIS1 association and function, we point at regulatory functions during cardiac development as common mechanisms to mitral valve degeneration.

Genetic Reduction in Left Ventricular Protein Kinase C-α and Adverse Ventricular Remodeling in Human Subjects

BACKGROUND

Inhibition of PKC-α (protein kinase C-α) enhances contractility and cardioprotection in animal models, but effects in humans are unknown. Genotypes at rs9912468 strongly associate with PRKCA expression in the left ventricle, enabling genetic approaches to measure effects of reduced PKC-α in human populations.

METHODS AND RESULTS

We analyzed the cis expression quantitative trait locus for PRKCA marked by rs9912468 using 313 left ventricular specimens from European Ancestry patients. The forward strand minor allele (G) at rs9912468 is associated with reduced PKC-α transcript abundance (1.7-fold reduction in minor allele homozygotes, P=1×10^-41). This association was cardiac specific in expression quantitative trait locus data sets that span 16 human tissues. Cardiac epigenomic data revealed a predicted enhancer in complete (R²=1.0) linkage disequilibrium with rs9912468 within intron 2 of PRKCA. We cloned this region and used reporter constructs to verify cardiac-specific enhancer activity in vitro in cardiac and noncardiac cells and in vivo in zebrafish. The PRKCA enhancer contains 2 common genetic variants and 4 haplotypes; the haplotype correlated with the rs9912468 PKC-α-lowering allele (G) showed lowest activity. In contrast to previous reports in animal models, the PKC-α-lowering allele is associated with adverse left ventricular remodeling (higher mass, larger diastolic dimension), reduced fractional shortening, and higher risk of dilated cardiomyopathy in human populations.

CONCLUSIONS

These findings support PKC-α as a regulator of the human heart but suggest that PKC-α inhibition may adversely affect the left ventricle depending on timing and duration. Pharmacological studies in human subjects are required to discern potential benefits and harms of PKC-α inhibitors as an approach to treat heart disease.

Association Between Titin Loss-of-Function Variants and Early-Onset Atrial Fibrillation

IMPORTANCE

Atrial fibrillation (AF) is the most common arrhythmia affecting 1% of the population. Young individuals with AF have a strong genetic association with the disease, but the mechanisms remain incompletely understood.

OBJECTIVE

To perform large-scale whole-genome sequencing to identify genetic variants related to AF.

DESIGN, SETTING, AND PARTICIPANTS

The National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine Program includes longitudinal and cohort studies that underwent high-depth whole-genome sequencing between 2014 and 2017 in 18 526 individuals from the United States, Mexico, Puerto Rico, Costa Rica, Barbados, and Samoa. This case-control study included 2781 patients with early-onset AF from 9 studies and identified 4959 controls of European ancestry from the remaining participants. Results were replicated in the UK Biobank (346 546 participants) and the MyCode Study (42 782 participants).

EXPOSURES

Loss-of-function (LOF) variants in genes at AF loci and common genetic variation across the whole genome.

MAIN OUTCOMES AND MEASURES

Early-onset AF (defined as AF onset in persons <66 years of age). Due to multiple testing, the significance threshold for the rare variant analysis was P = 4.55 × 10-3.

RESULTS

Among 2781 participants with early-onset AF (the case group), 72.1% were men, and the mean (SD) age of AF onset was 48.7 (10.2) years. Participants underwent whole-genome sequencing at a mean depth of 37.8 fold and mean genome coverage of 99.1%. At least 1 LOF variant in TTN, the gene encoding the sarcomeric protein titin, was present in 2.1% of case participants compared with 1.1% in control participants (odds ratio [OR], 1.76 [95% CI, 1.04-2.97]). The proportion of individuals with early-onset AF who carried a LOF variant in TTN increased with an earlier age of AF onset (P value for trend, 4.92 × 10-4), and 6.5% of individuals with AF onset prior to age 30 carried a TTN LOF variant (OR, 5.94 [95% CI, 2.64-13.35]; P = 1.65 × 10-5). The association between TTN LOF variants and AF was replicated in an independent study of 1582 patients with early-onset AF (cases) and 41 200 control participants (OR, 2.16 [95% CI, 1.19-3.92]; P = .01).

CONCLUSIONS AND RELEVANCE

In a case-control study, there was a statistically significant association between an LOF variant in the TTN gene and early-onset AF, with the variant present in a small percentage of participants with early-onset AF (the case group). Further research is necessary to understand whether this is a causal relationship.

Novel Mutation in FLNC (Filamin C) Causes Familial Restrictive Cardiomyopathy

BACKGROUND

Restrictive cardiomyopathy (RCM) is a rare cardiomyopathy characterized by impaired diastolic ventricular function resulting in a poor clinical prognosis. Rarely, heritable forms of RCM have been reported, and mutations underlying RCM have been identified in genes that govern the contractile function of the cardiomyocytes.

METHODS AND RESULTS

We evaluated 8 family members across 4 generations by history, physical examination, electrocardiography, and echocardiography. Affected individuals presented with a pleitropic syndrome of progressive RCM, atrioventricular septal defects, and a high prevalence of atrial fibrillation. Exome sequencing of 5 affected members identified a single novel missense variant in a highly conserved residue of FLNC (filamin C; p.V2297M). FLNC encodes filamin C-a protein that acts as both a scaffold for the assembly and organization of the central contractile unit of striated muscle and also as a mechanosensitive signaling molecule during cell migration and shear stress. Immunohistochemical analysis of FLNC localization in cardiac tissue from an affected family member revealed a diminished localization at the z disk, whereas traditional localization at the intercalated disk was preserved. Stem cell-derived cardiomyocytes mutated to carry the effect allele had diminished contractile activity when compared with controls.

CONCLUSION

We have identified a novel variant in FLNC as pathogenic variant for familial RCM-a finding that further expands on the genetic basis of this rare and morbid cardiomyopathy.

Multi-ethnic genome-wide association study for atrial fibrillation

Atrial fibrillation (AF) affects more than 33 million individuals worldwide1 and has a complex heritability2. We conducted the largest meta-analysis of genome-wide association studies (GWAS) for AF to date, consisting of more than half a million individuals, including 65,446 with AF. In total, we identified 97 loci significantly associated with AF, including 67 that were novel in a combined-ancestry analysis, and 3 that were novel in a European-specific analysis. We sought to identify AF-associated genes at the GWAS loci by performing RNA-sequencing and expression quantitative trait locus analyses in 101 left atrial samples, the most relevant tissue for AF. We also performed transcriptome-wide analyses that identified 57 AF-associated genes, 42 of which overlap with GWAS loci. The identified loci implicate genes enriched within cardiac developmental, electrophysiological, contractile and structural pathways. These results extend our understanding of the biological pathways underlying AF and may facilitate the development of therapeutics for AF.

Diminished PRRX1 Expression Is Associated With Increased Risk of Atrial Fibrillation and Shortening of the Cardiac Action Potential

BACKGROUND

Atrial fibrillation (AF) affects over 33 million individuals worldwide. Genome-wide association studies have identified at least 30 AF loci, but the mechanisms through which individual variants lead to altered disease risk have remained unclear for the majority of these loci. At the 1q24 locus, we hypothesized that the transcription factor PRRX1 could be a strong candidate gene as it is expressed in the pulmonary veins, a source of AF in many individuals. We sought to identify the molecular mechanism, whereby variation at 1q24 may lead to AF susceptibility.

METHODS AND RESULTS

We sequenced a ≈158 kb region encompassing PRRX1 in 962 individuals with and without AF. We identified a broad region of association with AF at the 1q24 locus. Using in silico prediction and functional validation, we identified an enhancer that interacts with the promoter of PRRX1 in cells of cardiac lineage. Within this enhancer, we identified a single-nucleotide polymorphism, rs577676, which alters enhancer activity in a mouse atrial cell line and in embryonic zebrafish and differentially regulates PRRX1 expression in human left atria. We found that suppression of PRRX1 in human embryonic stem cell-derived cardiomyocytes and embryonic zebrafish resulted in shortening of the atrial action potential duration, a hallmark of AF.

CONCLUSIONS

We have identified a functional genetic variant that alters PRRX1 expression, ultimately resulting in electrophysiological alterations in atrial myocytes that may promote AF.

Common and Rare Coding Genetic Variation Underlying the Electrocardiographic PR Interval

BACKGROUND

Electrical conduction from the cardiac sinoatrial node to the ventricles is critical for normal heart function. Genome-wide association studies have identified more than a dozen common genetic loci that are associated with PR interval. However, it is unclear whether rare and low-frequency variants also contribute to PR interval heritability.

METHODS

We performed large-scale meta-analyses of the PR interval that included 83 367 participants of European ancestry and 9436 of African ancestry. We examined both common and rare variants associated with the PR interval.

RESULTS

We identified 31 genetic loci that were significantly associated with PR interval after Bonferroni correction (P<1.2×10-6), including 11 novel loci that have not been reported previously. Many of these loci are involved in heart morphogenesis. In gene-based analysis, we found that multiple rare variants at MYH6 (P=5.9×10-11) and SCN5A (P=1.1×10-7) were associated with PR interval. SCN5A locus also was implicated in the common variant analysis, whereas MYH6 was a novel locus.

CONCLUSIONS

We identified common variants at 11 novel loci and rare variants within 2 gene regions that were significantly associated with PR interval. Our findings provide novel insights to the current understanding of atrioventricular conduction, which is critical for cardiac activity and an important determinant of health.

Gain-of-function mutations in GATA6 lead to atrial fibrillation

BACKGROUND

The genetic basis of atrial fibrillation (AF) and congenital heart disease remains incompletely understood.

OBJECTIVE

We sought to determine the causative mutation in a family with AF, atrial septal defects, and ventricular septal defects.

METHODS

We evaluated a pedigree with 16 family members, 1 with an atrial septal defect, 1 with a ventricular septal defect, and 3 with AF; we performed whole exome sequencing in 3 affected family members. Given that early-onset AF was prominent in the family, we then screened individuals with early-onset AF, defined as an age of onset <66 years, for mutations in GATA6. Variants were functionally characterized using reporter assays in a mammalian cell line.

RESULTS

Exome sequencing in 3 affected individuals identified a conserved mutation, R585L, in the transcription factor gene GATA6. In the Massachusetts General Hospital Atrial Fibrillation (MGH AF) Study, the mean age of AF onset was 47.1 ± 10.9 years; 79% of the participants were men; and there was no evidence of structural heart disease. We identified 3 GATA6 variants (P91S, A177T, and A543G). Using wild-type and mutant GATA6 constructs driving atrial natriuretic peptide promoter reporter, we found that 3 of the 4 variants had a marked upregulation of luciferase activity (R585L: 4.1-fold, P < .0001; P91S: 2.5-fold, P = .0002; A177T; 1.7-fold, P = .03). In addition, when co-overexpressed with GATA4 and MEF2C, GATA6 variants exhibited upregulation of the alpha myosin heavy chain and atrial natriuretic peptide reporter activity.

CONCLUSION

Overall, we found gain-of-function mutations in GATA6 in both a family with early-onset AF and atrioventricular septal defects as well as in a family with sporadic, early-onset AF.

A Functional Variant Associated with Atrial Fibrillation Regulates PITX2c Expression through TFAP2a

The most significantly associated genetic locus for atrial fibrillation (AF) is in chromosomal region 4q25, where four independent association signals have been identified. Although model-system studies suggest that altered PITX2c expression might underlie the association, the link between specific variants and the direction of effect on gene expression remains unknown for all four signals. In the present study, we analyzed the AF-associated region most proximal to PITX2 at 4q25. First, we identified candidate regulatory variants that might confer AF risk through a combination of mammalian conservation, DNase hypersensitivity, and histone modification from ENCODE and the Roadmap Epigenomics Project, as well as through in vivo analysis of enhancer activity in embryonic zebrafish. Within candidate regions, we then identified a single associated SNP, rs2595104, which displayed dramatically reduced enhancer activity with the AF risk allele. CRISPR-Cas9-mediated deletion of the rs2595104 region and editing of the rs2595104 risk allele in human stem-cell-derived cardiomyocytes resulted in diminished PITX2c expression in comparison to that of the non-risk allele. This differential activity was mediated by activating enhancer binding protein 2 alpha (TFAP2a), which bound robustly to the non-risk allele at rs2595104, but not to the risk allele, in cardiomyocytes. In sum, we found that the AF-associated SNP rs2595104 altered PITX2c expression via interaction with TFAP2a. Such a pathway could ultimately contribute to AF susceptibility at the PITX2 locus associated with AF.

Mutation of a common amino acid in NKX2.5 results in dilated cardiomyopathy in two large families

BACKGROUND

The genetic basis for dilated cardiomyopathy (DCM) can be difficult to determine, particularly in familial cases with complex phenotypes. Next generation sequencing may be useful in the management of such cases.

METHODS

We report two large families with pleiotropic inherited cardiomyopathy. In addition to DCM, the phenotypes included atrial and ventricular septal defects, cardiac arrhythmia and sudden death. Probands underwent whole exome sequencing to identify potentially causative variants.

RESULTS

Each whole exome sequence yielded over 18,000 variants. We identified distinct mutations affecting a common amino acid in NKX2.5. Segregation analysis of the families support the pathogenic role of these variants.

CONCLUSION

Our study emphasizes the utility of next generation sequencing in identifying causative mutations in complex inherited cardiac disease. We also report a novel pathogenic NKX2.5 mutation.

PHACTR1 Is a Genetic Susceptibility Locus for Fibromuscular Dysplasia Supporting Its Complex Genetic Pattern of Inheritance

Fibromuscular dysplasia (FMD) is a nonatherosclerotic vascular disease leading to stenosis, dissection and aneurysm affecting mainly the renal and cerebrovascular arteries. FMD is often an underdiagnosed cause of hypertension and stroke, has higher prevalence in females (~80%) but its pathophysiology is unclear. We analyzed ~26K common variants (MAF>0.05) generated by exome-chip arrays in 249 FMD patients and 689 controls. We replicated 13 loci (P<10-4) in 402 cases and 2,537 controls and confirmed an association between FMD and a variant in the phosphatase and actin regulator 1 gene (PHACTR1). Three additional case control cohorts including 512 cases and 669 replicated this result and overall reached the genomic level of significance (OR = 1.39, P = 7.4×10-10, 1,154 cases and 3,895 controls). The top variant, rs9349379, is intronic to PHACTR1, a risk locus for coronary artery disease, migraine, and cervical artery dissection. The analyses of geometrical parameters of carotids from ~2,500 healthy volunteers indicate higher intima media thickness (P = 1.97×10-4) and wall to lumen ratio (P = 0.002) in rs9349379-A carriers, suggesting indices of carotid hypertrophy previously described in carotids of FMD patients. Immunohistochemistry detected PHACTR1 in endothelium and smooth muscle cells of FMD and normal human carotids. The expression of PHACTR1 by genotypes in primary human fibroblasts showed higher expression in rs9349379-A carriers (N = 86, P = 0.003). Phactr1 knockdown in zebrafish resulted in dilated vessels indicating subtle impaired vascular development. We report the first susceptibility locus for FMD and provide evidence for a complex genetic pattern of inheritance and indices of shared pathophysiology between FMD and other cardiovascular and neurovascular diseases.

Discovery and validation of sub-threshold genome-wide association study loci using epigenomic signatures

Genetic variants identified by genome-wide association studies explain only a modest proportion of heritability, suggesting that meaningful associations lie 'hidden' below current thresholds. Here, we integrate information from association studies with epigenomic maps to demonstrate that enhancers significantly overlap known loci associated with the cardiac QT interval and QRS duration. We apply functional criteria to identify loci associated with QT interval that do not meet genome-wide significance and are missed by existing studies. We demonstrate that these 'sub-threshold' signals represent novel loci, and that epigenomic maps are effective at discriminating true biological signals from noise. We experimentally validate the molecular, gene-regulatory, cellular and organismal phenotypes of these sub-threshold loci, demonstrating that most sub-threshold loci have regulatory consequences and that genetic perturbation of nearby genes causes cardiac phenotypes in mouse. Our work provides a general approach for improving the detection of novel loci associated with complex human traits.

Discovery and validation of sub-threshold genome-wide association study loci using epigenomic signatures

Genetic variants identified by genome-wide association studies explain only a modest proportion of heritability, suggesting that meaningful associations lie 'hidden' below current thresholds. Here, we integrate information from association studies with epigenomic maps to demonstrate that enhancers significantly overlap known loci associated with the cardiac QT interval and QRS duration. We apply functional criteria to identify loci associated with QT interval that do not meet genome-wide significance and are missed by existing studies. We demonstrate that these 'sub-threshold' signals represent novel loci, and that epigenomic maps are effective at discriminating true biological signals from noise. We experimentally validate the molecular, gene-regulatory, cellular and organismal phenotypes of these sub-threshold loci, demonstrating that most sub-threshold loci have regulatory consequences and that genetic perturbation of nearby genes causes cardiac phenotypes in mouse. Our work provides a general approach for improving the detection of novel loci associated with complex human traits.

DOI:http://dx.doi.org/10.7554/eLife.10557.001

Most complex traits are governed by a large number of genetic contributors, each playing only a modest effect. This makes it difficult to identify the genetic variants that increase disease risk, hindering the discovery of new drug targets and the development of new therapeutics.

To overcome this limitation in discovery power, the field of human genetics has traditionally sought increasingly large groups, or cohorts, of afflicted and non-afflicted individuals. Studies of large cohorts are a powerful approach for discovering new disease genes, but such groups are often impractical and sometimes impossible to obtain. However, it has become possible to complement the genetic evidence found in disease association studies with biological evidence of the effects of disease-associated genetic variants.

Wang et al. focus specifically on genetic sites, or loci, that do not affect protein sequence but instead affect the non-coding control regions. These are known as enhancer elements, as they can enhance the expression of nearby genes. These loci constitute the majority of disease regions, and thus are extremely important, but their discovery has been hindered by our relatively poor understanding of the human genome.

Chemical modifications known as epigenomic marks are indicative of enhancer regions. By studying the factors that affect heart rhythm, Wang et al. show that specific combinations of epigenomic marks are enriched in known trait-associated regions. This knowledge was then used to prioritize the further investigation of genetic regions that genome-wide association studies had only weakly linked to heart rhythm alterations. Wang et al. directly confirmed that genetic differences in “sub-threshold” regions indeed alter the activity of these regulatory regions in human heart cells. Furthermore, mutating or perturbing the predicted target genes of the sub-threshold enhancers caused heart defects in mouse and zebrafish.

Wang et al. have demonstrated that epigenome maps can help to distinguish which sub-threshold regions from genome-wide association studies are more likely to contribute to a disease. This allows for the discovery of new disease genes with much smaller cohorts than would be needed otherwise, thus speeding up the development of new therapeutics by many years.

DOI:http://dx.doi.org/10.7554/eLife.10557.002

Common variation in atrial fibrillation: navigating the path from genetic association to mechanism

Atrial fibrillation (AF) is the most common cardiac arrhythmia with well-established clinical and genetic risk components. Genome-wide association studies (GWAS) have identified 17 independent susceptibility signals for AF at 14 genomic regions, but the mechanisms through which these loci confer risk to AF remain largely undefined. This problem is not unique to AF, as the field of functional genomics, which attempts to bridge this gap from genotype to phenotype, has only uncovered the mechanisms for a handful of GWAS loci. Recent functional genomic studies have made great strides towards translating genetic discoveries to an underlying mechanism, but the large-scale application of these techniques to AF has remain limited. These advances, as well as the continued unresolved challenges for both common variation in AF and the functional genomics field in general, will be the subject of the following review.

Integrating genetic, transcriptional, and functional analyses to identify 5 novel genes for atrial fibrillation

BACKGROUND

Atrial fibrillation (AF) affects >30 million individuals worldwide and is associated with an increased risk of stroke, heart failure, and death. AF is highly heritable, yet the genetic basis for the arrhythmia remains incompletely understood.

METHODS AND RESULTS

To identify new AF-related genes, we used a multifaceted approach, combining large-scale genotyping in 2 ethnically distinct populations, cis-eQTL (expression quantitative trait loci) mapping, and functional validation. Four novel loci were identified in individuals of European descent near the genes NEURL (rs12415501; relative risk [RR]=1.18; 95% confidence interval [CI], 1.13-1.23; P=6.5×10(-16)), GJA1 (rs13216675; RR=1.10; 95% CI, 1.06-1.14; P=2.2×10(-8)), TBX5 (rs10507248; RR=1.12; 95% CI, 1.08-1.16; P=5.7×10(-11)), and CAND2 (rs4642101; RR=1.10; 95% CI, 1.06-1.14; P=9.8×10(-9)). In Japanese, novel loci were identified near NEURL (rs6584555; RR=1.32; 95% CI, 1.26-1.39; P=2.0×10(-25)) and CUX2 (rs6490029; RR=1.12; 95% CI, 1.08-1.16; P=3.9×10(-9)). The top single-nucleotide polymorphisms or their proxies were identified as cis-eQTLs for the genes CAND2 (P=2.6×10(-19)), GJA1 (P=2.66×10(-6)), and TBX5 (P=1.36×10(-5)). Knockdown of the zebrafish orthologs of NEURL and CAND2 resulted in prolongation of the atrial action potential duration (17% and 45%, respectively).

CONCLUSIONS

We have identified 5 novel loci for AF. Our results expand the diversity of genetic pathways implicated in AF and provide novel molecular targets for future biological and pharmacological investigation.

Emerging directions in the genetics of atrial fibrillation

Atrial fibrillation (AF) is the most common arrhythmia and is associated with increased morbidity. As the population ages and the prevalence of AF continues to rise, the socioeconomic consequences of AF will become increasingly burdensome. Although there are well-defined clinical risk factors for AF, a significant heritable component is also recognized. To identify the molecular basis for the heritability of AF, investigators have used a combination of classical Mendelian genetics, candidate gene screening, and genome-wide association studies. However, these avenues have, as yet, failed to define the majority of the heritability of AF. The goal of this review is to describe the results from both candidate gene and genome-wide studies, as well as to outline potential future avenues for creating a more complete understanding of AF genetics. Ultimately, a more comprehensive view of the genetic underpinnings for AF will lead to the identification of novel molecular pathways and improved risk prediction of this complex arrhythmia.

A novel trafficking-defective HCN4 mutation is associated with early-onset atrial fibrillation

BACKGROUND

Atrial fibrillation (AF) is the most common arrhythmia, and a recent genome-wide association study identified the hyperpolarization-activated cyclic nucleotide-gated channel 4 (HCN4) as a novel AF susceptibility locus. HCN4 encodes for the cardiac pacemaker channel, and HCN4 mutations are associated with familial sinus bradycardia and AF.

OBJECTIVE

The purpose of this study was to determine whether novel variants in the coding region of HCN4 contribute to the susceptibility for AF.

METHODS

We sequenced the coding region of HCN4 for novel variants from 527 cases with early-onset AF from the Massachusetts General Hospital AF Study and 443 referents from the Framingham Heart Study. We used site-directed mutagenesis, cellular electrophysiology, immunocytochemistry, and confocal microscopy to functionally characterize novel variants.

RESULTS

We found the frequency of novel coding HCN4 variants was 2-fold greater for individuals with AF (7 variants) compared to the referents (3 variants). We determined that of the 7 novel HCN4 variants in our AF cases, 1 (p.Pro257Ser, located in the amino-terminus adjacent to the first transmembrane spanning domain) did not traffic to cell membrane, whereas the remaining 6 were not functionally different from wild type. In addition, the 3 novel variants in our referents did not alter function compared to wild-type. Coexpression studies showed that the p.Pro257Ser mutant channel failed to colocalize with the wild-type HCN4 channel on the cell membrane.

CONCLUSION

Our findings are consistent with HCN4 haploinsufficiency as the likely mechanism for early-onset AF in the p.Pro257Ser carrier.

Overexpression of KCNN3 results in sudden cardiac death

BACKGROUND

A recent genome-wide association study identified a susceptibility locus for atrial fibrillation at the KCNN3 gene. Since the KCNN3 gene encodes for a small conductance calcium-activated potassium channel, we hypothesized that overexpression of the SK3 channel increases susceptibility to cardiac arrhythmias.

METHODS AND RESULTS

We characterized the cardiac electrophysiological phenotype of a mouse line with overexpression of the SK3 channel. We generated homozygote (SK3(T/T)) and heterozygote (SK3(+/T)) mice with overexpression of the channel and compared them with wild-type (WT) controls. We observed a high incidence of sudden death among SK3(T/T) mice (7 of 19 SK3(T/T) mice). Ambulatory monitoring demonstrated that sudden death was due to heart block and bradyarrhythmias. SK3(T/T) mice displayed normal body weight, temperature, and cardiac function on echocardiography; however, histological analysis demonstrated that these mice have abnormal atrioventricular node morphology. Optical mapping demonstrated that SK3(T/T) mice have slower ventricular conduction compared with WT controls (SK3(T/T) vs. WT; 0.45 ± 0.04 vs. 0.60 ± 0.09 mm/ms, P = 0.001). Programmed stimulation in 1-month-old SK3(T/T) mice demonstrated inducible atrial arrhythmias (50% of SK3(T/T) vs. 0% of WT mice) and also a shorter atrioventricular nodal refractory period (SK3(T/T) vs. WT; 43 ± 6 vs. 52 ± 9 ms, P = 0.02). Three-month-old SK3(T/T) mice on the other hand displayed a trend towards a more prolonged atrioventricular nodal refractory period (SK3(T/T) vs. WT; 61 ± 1 vs. 52 ± 6 ms, P = 0.06).

CONCLUSION

Overexpression of the SK3 channel causes an increased risk of sudden death associated with bradyarrhythmias and heart block, possibly due to atrioventricular nodal dysfunction.

Abstract 285: Overexpression of SK3 in a Murine Model Results in Sudden Cardiac Death

Background: Recently, the gene responsible for encoding SK3, KCNN3, was implicated in a genome-wide association study as a susceptibility locus for atrial fibrillation. SK3 is one of a family of three small conductance, voltage-independent, calcium-activated potassium channels which is abundantly expressed in cardiac tissues. Therefore, we hypothesized that transcriptional misregulation of the KCNN3 gene could result in arrhythmogenesis.

Methods: We characterized a mouse line which had constitutive ubiquitous overexpression of Kcnn3 using ambulatory cardiac rhythm monitoring, optical mapping and in vivo electrophysiology studies. Both homozygote (SK3 T/T ) and heterozygote (SK3 T/+ ) animals were compared to their wild-type (SK3 +/+ ) littermates for assessment of cardiac phenotypes.

Results: SK3 T/T mice displayed increased expression of SK3, as assessed by qRT-PCR, in all four chambers of the heart, although highest increase of expression was observed in the ventricles. Examinations of SK3 T/T mice revealed no gross morphological changes in the myocardium, no observable cardiac fibrosis, and normal echocardiograms when compared to SK3 +/+ mice. However, pups from SK3 T/+ crosses did not produce expected Mendelian ratios (SK3 +/+ = 17%), and 7 out of 19 SK3 +/+ mice died suddenly by 3 months of age, whereas all of the SK3 T/+ or SK3 +/+ survived. To address potential mechanisms, we performed ambulatory monitoring on mice beginning at approximately one month of age. Of the 6 SK3 +/+ mice, 4 died suddenly, whereas neither of the wild-type controls expired. The cardiac rhythm recorded at the time of death in three of the four mice was heart block followed by severe bradycardia. Additionally, during the period of ambulatory monitoring, SK3 T/T mice had a lower mean heart rate compared to the SK3 +/+ control mice, and SK3 T/T mice also displayed frequent episodes of atrioventricular dissociation, both at rest and during periods of activity. SK3 T/T mice also displayed more pronounced variability of the heart rate and the PR interval. Optical mapping revealed slower ventricular conduction velocity.

Conclusion: These data suggest a mechanism whereby overexpression of SK3 leads to sudden death due to bradyarrhythmias and heart block.

HSF1 is essential for the resistance of zebrafish eye and brain tissues to hypoxia/reperfusion injury

Ischemia and subsequent reperfusion (IR) produces injury to brain, eye and other tissues, contributing to the progression of important clinical pathologies. The response of cells to IR involves activation of several signaling pathways including those activating hypoxia and heat shock responsive transcription factors. However, specific roles of these responses in limiting cell damage and preventing cell death after IR have not been fully elucidated. Here, we have examined the role of heat shock factor 1 (HSF1) in the response of zebrafish embryos to hypoxia and subsequent return to normoxic conditions (HR) as a model for IR. Heat shock preconditioning elevated heat shock protein expression and protected zebrafish embryo eye and brain tissues against HR-induced apoptosis. These effects were inhibited by translational suppression of HSF1 expression. Reduced expression of HSF1 also increased cell death in brain and eye tissues of embryos subjected to hypoxia and reperfusion without prior heat shock. Surprisingly, reduced expression of HSF1 had only a modest effect on hypoxia-induced expression of Hsp70 and no effect on hypoxia-induced expression of Hsp27. These results establish the zebrafish embryo as a model for the study of ischemic injury in the brain and eye and reveal a critical role for HSF1 in the response of these tissues to HR. Our results also uncouple the role of HSF1 expression from that of Hsp27, a well characterized heat shock protein considered essential for cell survival after hypoxia. Alternative roles for HSF1 are considered.

Hsp27 is persistently expressed in zebrafish skeletal and cardiac muscle tissues but dispensable for their morphogenesis

Constitutive expression of Hsp27 has been demonstrated in vertebrate embryos, especially in developing skeletal and cardiac muscle. Results of several previous studies have indicated that Hsp27 could play a role in the development of these tissues. For example, inhibition of Hsp27 expression has been reported to cause defective development of mammalian myoblasts in vitro and frog embryos in vivo. In contrast, transgenic mice lacking Hsp27 develop normally. Here, we examined the distribution of Hsp27 protein in developing and adult zebrafish and effects of suppressing Hsp27 expression using phosphorodiamidate morpholino oligonucleotides (PMO) on zebrafish development. Consistent with our previous analysis of hsp27 messenger RNA expression, we detected the protein Hsp27 in cardiac, smooth, and skeletal muscle of both embryonic and adult zebrafish. However, embryos lacking detectable Hsp27 after injection of antisense hsp27 PMO exhibited comparable heart beat rates to that of control embryos and cardiac morphology was indistinguishable in the presence or absence of Hsp27. Loss of Hsp27 also had no effect on the structure of the skeletal muscle myotomes in the developing embryo. Finally, embryos injected with antisense hsp27 and scrambled control PMO displayed equal motility. We conclude that Hsp27 is dispensable for zebrafish morphogenesis but could play a role in long-term maintenance of heart and muscle tissues.