Search

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.

Regulation of transcript structure generates transcript diversity and plays an important role in human disease

Cardiac allograft vasculopathy (CAV) is the leading cause of late allograft failure and mortality after heart transplantation. As current standards of diagnosis and treatment of CAV have significant limitations, understanding cell-specific responses may prove critical for developing improved detection strategies and novel therapeutics. This study is the first to successfully utilize human endomyocardial biopsy (EMB) samples to isolate large numbers of intact nuclei for single-nuclear transcriptomics. These data also lay the groundwork for ongoing experiments to study serial, routinely-collected EMB specimens after heart transplantation to identify novel biomarkers and pathways through which early CAV pathogenesis can be interrupted, thereby prolonging allograft survival.

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.

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.

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.\ud \ud 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.\ud \ud 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.

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.

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.

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.

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 Pitx2 c 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 .

New York State, in particular the New York City metropolitan area, was the early epicenter of the SARS-CoV-2 pandemic in the United States. Similar to initial pandemic dynamics in many metropolitan areas, multiple introductions from various locations appear to have contributed to the swell of positive cases. However, representation and analysis of samples from New York regions outside the greater New York City area were lacking, as were SARS-CoV-2 genomes from the earliest cases associated with the Westchester County outbreak, which represents the first outbreak recorded in New York State. The Wadsworth Center, the public health laboratory of New York State, sought to characterize the transmission dynamics of SARS-CoV-2 across the entire state of New York from March to September with the addition of over 600 genomes from under-sampled and previously unsampled New York counties and to more fully understand the breadth of the initial outbreak in Westchester County. Additional sequencing confirmed the dominance of B.1 and descendant lineages (collectively referred to as B.1.X) in New York State. Community structure, phylogenetic, and phylogeographic analyses suggested that the Westchester outbreak was associated with continued transmission of the virus throughout the state, even after travel restrictions and the on-pause measures of March, contributing to a substantial proportion of the B.1 transmission clusters as of September 30th, 2020.

SummaryRegulation of transcript structure generates transcript diversity and plays an important role in human disease. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure. In this paper, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from GTEx tissues and cell lines, complementing the GTEx resource. We identified just under 100,000 new transcripts for annotated genes, and validated the protein expression of a similar proportion of novel and annotated transcripts. We developed a new computational package, LORALS, to analyze genetic effects of rare and common variants on the transcriptome via allele-specific analysis of long reads. We called allele-specific expression and transcript structure events, providing novel 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 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 use this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.

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 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 −16 ). Chromatin states were differentially methylated, promoters were less methylated than repressed regions ( P −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.

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 −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.

The aorta is the largest blood vessel in the body, and enlargement or aneurysm of the aorta can predispose to dissection, an important cause of sudden death. While rare syndromes have been identified that predispose to aortic aneurysm, the common genetic basis for the size of the aorta remains largely unknown. By leveraging a deep learning architecture that was originally developed to recognize natural images, we trained a model to evaluate the dimensions of the ascending and descending thoracic aorta in cardiac magnetic resonance imaging. After manual annotation of just 116 samples, we applied this model to 3,840,140 images from the UK Biobank. We then conducted a genome-wide association study in 33,420 individuals, revealing 68 loci associated with ascending and 35 with descending thoracic aortic diameter, of which 10 loci overlapped. Integration of common variation with transcriptome-wide analyses, rare-variant burden tests, and single nucleus RNA sequencing prioritized SVIL, a gene highly expressed in vascular smooth muscle, that was significantly associated with the diameter of the ascending and descending aorta. A polygenic score for ascending aortic diameter was associated with a diagnosis of thoracic aortic aneurysm in the remaining 391,251 UK Biobank participants who did not undergo imaging (HR = 1.44 per standard deviation; P = 3.7·10−12). Defining the genetic basis of the diameter of the aorta may enable the identification of asymptomatic individuals at risk for aneurysm or dissection and facilitate the prioritization of potential therapeutic targets for the prevention or treatment of aortic aneurysm. Finally, our results illustrate the potential for rapidly defining novel quantitative traits derived from a deep learning model, an approach that can be more broadly applied to biomedical imaging data.

The COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, creates an urgent need for identifying molecular mechanisms that mediate viral entry, propagation, and tissue pathology. Cell membrane bound angiotensin-converting enzyme 2 (ACE2) and associated proteases, transmembrane protease serine 2 (TMPRSS2) and Cathepsin L (CTSL), were previously identified as mediators of SARS-CoV2 cellular entry. Here, we assess the cell type-specific RNA expression of ACE2, TMPRSS2, and CTSL through an integrated analysis of 107 single-cell and single-nucleus RNA-Seq studies, including 22 lung and airways datasets (16 unpublished), and 85 datasets from other diverse organs. Joint expression of ACE2 and the accessory proteases identifies specific subsets of respiratory epithelial cells as putative targets of viral infection in the nasal passages, airways, and alveoli. Cells that co-express ACE2 and proteases are also identified in cells from other organs, some of which have been associated with COVID-19 transmission or pathology, including gut enterocytes, corneal epithelial cells, cardiomyocytes, heart pericytes, olfactory sustentacular cells, and renal epithelial cells. Performing the first meta-analyses of scRNA-seq studies, we analyzed 1,176,683 cells from 282 nasal, airway, and lung parenchyma samples from 164 donors spanning fetal, childhood, adult, and elderly age groups, associate increased levels of ACE2, TMPRSS2, and CTSL in specific cell types with increasing age, male gender, and smoking, all of which are epidemiologically linked to COVID-19 susceptibility and outcomes. Notably, there was a particularly low expression of ACE2 in the few young pediatric samples in the analysis. Further analysis reveals a gene expression program shared by ACE2+TMPRSS2+ cells in nasal, lung and gut tissues, including genes that may mediate viral entry, subtend key immune functions, and mediate epithelial-macrophage cross-talk. Amongst these are IL6, its receptor and co-receptor, IL1R, TNF response pathways, and complement genes. Cell type specificity in the lung and airways and smoking effects were conserved in mice. Our analyses suggest that differences in the cell type-specific expression of mediators of SARS-CoV-2 viral entry may be responsible for aspects of COVID-19 epidemiology and clinical course, and point to putative molecular pathways involved in disease susceptibility and pathogenesis.

Publisher's version (útgefin grein), 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., We provide all investigator and study-specific acknowledgements in Supplementary Note 1, and funding sources in Supplementary Note 2.

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.

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.

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality1,2. We performed multi-ancestry (N=293,051) and European only (N=271,570) genome-wide association (GWAS) meta-analyses for the PR interval, discovering 210 loci of which 149 are novel. Variants at all loci nearly doubled the percentage of heritability explained, from 33.5% to 62.6%. We observed enrichment for genes involved in cardiac muscle development/contraction and the cytoskeleton highlighting key regulation processes for atrioventricular conduction. Additionally, 19 novel loci harbour genes underlying inherited monogenic heart diseases suggesting the role of these genes in cardiovascular pathology in the general population. We showed that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease risk, including distal conduction disease, AF, atrioventricular pre-excitation, non-ischemic cardiomyopathy, and coronary heart disease. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.

SummaryMost disease-associated variants identified by population based genetic studies are non-coding, which compromises finding causative genes and mechanisms. Presumably they interact through looping with nearby genes to modulate transcription. Hi-C provides the most complete and unbiased method for genome-wide identification of potential regulatory interactions, but finding chromatin loops in Hi-C data remains difficult and tissue specific data are limited. We have generated Hi-C data from primary cardiac tissue and developed a method, peakHiC, for sensitive and quantitative loop calling to uncover the human heart regulatory interactome. We identify complex CTCF-dependent and -independent contact networks, with loops between coding and non-coding gene promoters, shared enhancers and repressive sites. Across the genome, enhancer interaction strength correlates with gene transcriptional output and loop dynamics follows CTCF, cohesin and H3K27Ac occupancy levels. Finally, we demonstrate that intersection of the human heart regulatory interactome with cardiovascular disease variants facilitates prioritizing disease-causative genes.

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.

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.

Nonsyndromic Mitral valve prolapse (MVP) is a common degenerative valvular heart disease with severe health consequences, including arrhythmia, heart failure and sudden death. MVP is characterized by excess extracellular matrix secretion and cellular disorganization which leads to bulky valves that are unable to co-apt properly during ventricular systole. However, the triggering mechanisms of this process are mostly unknown. Using pathway enrichment tools applied to GWAS we show that genes at risk loci are involved in biological functions relevant to cell adhesion and migration during cardiac development and in response to shear stress. Through genetic, in silico and in vivo experiments we demonstrates the presence of several genes involved in gene regulation, including GLIS1, a transcription factor that regulates Hedgehog signaling. Our findings define genetic, molecular and cellular mechanisms underlying non-syndromic MVP and implicate disrupted endothelial to mesenchymal transition and cell migration as a potential common cause to this disease.

Recently, we reported a 4-generation family of European ancestry with restrictive cardiomyopathy and a high prevalence of atrial fibrillation.1 By exome sequencing, we identified a missense variant in affected family members in a highly conserved residue (p.V2297M) of FLNC, the gene encoding filamin C. We …

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 2 =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.

Atrial fibrillation (AF) affects more than 33 million individuals worldwide 1 and has a complex heritability 2 . 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.

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.

Atrial fibrillation affects more than 33 million people worldwide and increases the risk of stroke, heart failure, and death(1,2). Fourteen genetic loci have been associated with atrial fibrillation in European and Asian ancestry groups(3-7). To further define the genetic basis of atrial fibrillation, we performed large-scale, trans-ancestry meta-analyses of common and rare variant association studies. The genome-wide association studies (GWAS) included 17,931 individuals with atrial fibrillation and 115,142 referents; the exome-wide association studies (ExWAS) and rare variant association studies (RVAS) involved 22,346 cases and 132,086 referents. We identified 12 new genetic loci that exceeded genome-wide significance, implicating genes involved in cardiac electrical and structural remodeling. Our results nearly double the number of known genetic loci for atrial fibrillation, provide insights into the molecular basis of atrial fibrillation, and may facilitate the identification of new potential targets for drug discovery(8).

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.

Atrial fibrillation (AF) is the most common arrhythmic disorder, and currently affects nearly 3 million Americans, 8.8 million Europeans, and an estimated 30 million individuals worldwide. The clinical risk factors for AF are numerous, with age, sex, hypertension, obesity, and ischemic heart disease among the most prevalent. Over the last ten years, a preponderance of evidence also suggests a large genetic contribution to AF. The earliest report of familial AF dates to the early 1940s1. Since then, it has become apparent that AF in referral populations2,3 and in the community is heritable4,5. Indeed, having a family member with AF is associated with a 40% increased risk for the arrhythmia6. Once the heritability was recognized, traditional genetics techniques for the discovery of rare, monogenic causes of AF were used to identify the initial AF genes. These studies in turn, informed candidate gene screening in AF cohorts. To identify additional sources of heritability for AF, large-scale analyses of common variation through genome wide association studies (GWAS) has recently yielded data identifying risk loci in many regions of the genome. In spite of these advances, the combination of these techniques has, as yet, failed to completely identify the heritability of AF in the population. It is the goal of this review to examine the previous studies on rare variants, address the findings of the recent GWAS studies, and describe future avenues towards defining the heritability of AF.

Objective: Nonsyndromic Mitral valve prolapse (MVP) is a common degenerative valvular heart disease with severe health consequences, including arrhythmia, heart failure and sudden death. MVP is cha...

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 ( SK3T/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 SK3T/T mice (7 of 19 SK3T/T mice). Ambulatory monitoring demonstrated that sudden death was due to heart block and bradyarrhythmias. SK3T/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 SK3T/T mice have slower ventricular conduction compared with WT controls ( SK3T/T vs. WT; 0.45 ± 0.04 vs. 0.60 ± 0.09 mm/ms, P = 0.001). Programmed stimulation in 1-month-old SK3T/T mice demonstrated inducible atrial arrhythmias (50% of SK3T/T vs. 0% of WT mice) and also a shorter atrioventricular nodal refractory period ( SK3T/T vs. WT; 43 ± 6 vs. 52 ± 9 ms, P = 0.02). Three-month-old SK3T/T mice on the other hand displayed a trend towards a more prolonged atrioventricular nodal refractory period ( SK3T/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.

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 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.

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 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 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.

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, Author Summary Fibromuscular Dysplasia (FMD) is a vascular disease characterized by a succession of occlusions and dilatation of medium-sized arteries (e.g renal, carotid or coronary arteries) with important health consequences, mainly resistant hypertension and stroke. FMD is an atypical vascular disease because it is not associated with overweight or dyslipidemia and 80% of patients are early middle aged women. Our genetic study conducted in >1100 patients and >3800 controls demonstrate that a common variant rs9349379 located on chromosome 6 in the phosphatase and actin regulator 1 gene (PHACTR1) increases by ~40% the risk of FMD. This is the first time a genetic risk factor is reported for FMD because it has been longtime considered rare and potentially under a Mendelian mode of inheritance. We also show that rs9349379 correlates with the expression of PHACTR1 in fibroblasts from FMD patients and controls. Interestingly, the same allele that increases the risk of FMD is at risk for cervical artery dissection and migraine, often reported in FMD patients but protective from myocardial infarction and coronary disease, where atherosclerosis is more common. The clear role of PHACTR1 in maintaining vascular well integrity is not fully elucidated. Using a specific antibody we detected PHACTR1 both on endothelial and smooth muscle cells of human FMD and control carotids, which suggests that PHACTR1 may have multiple functions depending on the cell type and the degree of atherosclerosis of the arteries.

Injury to muscle tissue plays a central role in various cardiovascular pathologies. Overexpression of the small heat shock protein Hsp27 protects muscle cells against thermal, oxidative and ischemic stress. However, underlying mechanisms of this protection have not been resolved. A distinctive feature of muscle cells is the stress-induced association of Hsp27 with the sarcomere. The association of Hsp27 with the cytoskeleton, in both muscle and non-muscle cells, is thought to represent interaction with Z-line components or filamentous actin. Here, we examined the association of Hsp27 with myofibrils in adult zebrafish myocardium subjected to hyperthermia and mechanical stretching. Consistent with previously published results, Hsp27 in resting length myofibrils localized to narrowly defined regions, or bands, which colocalized with Z-line markers. However, analysis of stretched myofibrils revealed that the association of Hsp27 with myofibrils was independent of desmin, alpha-actinin, myosin, and filamentous actin. Instead, Hsp27 maintained a consistent relationship with a marker for the titin A/I border over various sarcomeric lengths. Finally, extraction of actin filaments revealed that Hsp27 binds to a component of the remaining sarcomere. Together, these novel data support a mechanism of Hsp27 function where interactions with the titin filament system protect myofibrils from stress-induced degradation.

Background: To date, 17 loci associated with atrial fibrillation (AF) have been identified, and there are 4 regions at 4q25 around PITX2 gene strongly associated with AF. However, the functional variants at all known AF loci remain elusive. Methods and Results: Among the 4 AF loci at 4q25, we focused on one associated region of ~84 kb extending from Chr4: 111,520,926 to 111,604,727 that includes the PITX2 gene itself (defined as all SNPs with r 2 > 0.6 from the sentinel SNP rs1448818). To identify potential functional variants , we analyzed DNase hypersensitivity, histone modification and mammalian conservation from ENCODE and Roadmap Epigenomics Project. Additionally, cap analysis of gene expression (CAGE) data was also investigated. We identified 13 candidate regions, 5 of which showed enhancer activity in zebrafish skeletal muscle and heart. Using luciferase assays in a mouse atrial cell line (HL-1), we then identified one SNP (rs2595104) in which the T risk allele was associated with 31% decreased activity compared to the G non-risk allele (p = 0.007). rs2595104 lies ~10 kb upstream of PITX2c transcriptional start site and is in high linkage disequilibrium (LD) (r 2 = 0.67) with rs1448818. In silico investigation of rs2595104 using UniPROBE and HaploReg database revealed differential binding of each allele to activating enhancer binding protein 2 alpha (TFAP2a). Electrophoretic mobility shift assays showed that TFAP2a bound robustly to the non-risk G allele at rs2595104 but not to the risk T allele. Chromatin immunopricipitation followed by allele-specific qPCR in human embryonic stem cell-derived cardiomyocytes (hESC-CM) further suggested that TFAP2a was preferentially recruited to the non-risk G allele at rs2595104 in an AF-relevant cell type. hESC-CMs with homozygous deletion of rs2595104-containing enhancer using CRSPR-Cas9 expressed less PITX2c (35 ± 6%, p = 0.036), compared to the wild-type hESC-CMs. Finally, mutating the G allele at rs2595104 to T allele by CRSPR-Cas9 led to 21% reduction in PITX2c level in hESC-CMs (p = 0.027). Conclusion: We have found that the AF-associated SNP rs2595104 alters PITX2c expression via interaction with TFAP2a. Such a pathway may ultimately contribute to AF susceptibility at the PITX2 locus.