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.

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 Mb, respectively. 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, the brown bear is highly syntenic with the polar bear, and our phylogenetic analysis of 7,246 single-copy BUSCOs 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.SignificanceBrown bears (Ursus arctos) are the most widespread, large terrestrial carnivore on the planet and represent an interesting example of speciation through hybridization, as well as a medical model for sedentary lifestyle-related disease. Although a previous genome for a brown bear has been published, the reported contig N50 was low (only ∼530 kb), despite being scaffolded into putative chromosomes. Genomes of this quality limit the accuracy of analyses which rely on long contiguous stretches of the genome to be assembled (such as with many demographic analyses) as well as attempts at connecting genotype to phenotype (such as in association analyses). In order to support studies on both the complex hybridization history of the brown bear and investigations into medically-relevant phenotypes, we generated a fully-phased, chromosome-level assembly from a male grizzly bear. The genome has a total size of 2.47 Gb and 90% of the genome is contained in 36 scaffolds, roughly corresponding to one autosome per scaffold. This high-quality genome will enable studies across a variety of disciplines, including conservation, evolution, and medicine.

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.

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.

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.

Introduction: Approximately 50% of all heart failure (HF) diagnoses can be classified as HF with preserved ejection fraction (HFpEF), which has no FDA approved therapies. HFpEF is more prevalent in females compared to males, but the underlying mechanisms for the development as a sex-based disorder are unknown. We previously described how slow progressive pressure overload (PO) in male felines recapitulates the HFpEF phenotype but have not investigated the female phenotype. Hypothesis: Females will develop a less severe HFpEF phenotype compared to males under the same pathological stress. Methods & Results: Male (m) and female (f) domestic short felines (age 2mo) underwent aortic constriction (m: n=11; f: n=10) using a customized pre-shaped band or a sham procedure (m: n=7; f: n=7). Before surgery (baseline), there was no difference in body weight (BW) between groups and lung compliance was not different. Echocardiography revealed no significant difference in the ratio of left atrium to aortic root (LA/Ao), LA ejection fraction (LA EF), left ventricle (LV) ejection fraction, LV wall-thickness, and E/A ratio. By 4mo post-surgery, both males and females had developed cardiac dysfunction and decreased lung compliance. At this time, females had significantly smaller BWs than males. Despite the difference in BW, LV wall thickness and changes in E/A ratio were similar in both sexes in banded vs. shams. Importantly, LV EF did not change in any group. There was a decrease in LA EF and increased LA/Ao in all banded animals. Invasive hemodynamics at 4mo post-surgery showed no differences between sexes for the systolic pressure gradient generated by aortic banding. Banded males had higher LV end-diastolic pressure vs. banded females (m: 15.0±2.7mmHg; f: 8.1±1.9mmHg). However, there was a trend towards prolongation of tau and lower dp/dt min in banded females, suggestive of abnormal relaxation. There were no differences between banded males and females in heart weight/BW or cardiomyocyte cross-sectional area and both developed fibrosis. Conclusions: Exposure of male and female felines to PO resulted in similar cardiac hypertrophy, fibrosis, and decreased lung compliance. Females had lower LVEDP than males suggesting they may be protected from diastolic dysfunction.

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 .

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.

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.

Heart failure encompasses a heterogeneous set of clinical features that converge on impaired cardiac contractile function

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.

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.

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.

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.

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.

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.

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

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

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.

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.

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

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.

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.

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.

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.