Transcriptional profile of the rat cardiovascular system at single cell resolution.

Background: Despite the critical role of the cardiovascular system, our understanding of its cellular and transcriptional diversity remains limited. We therefore sought to characterize the cellular composition, phenotypes, molecular pathways, and communication networks between cell types at the tissue and sub-tissue level across the cardiovascular system of the healthy Wistar rat, an important model in preclinical cardiovascular research. We obtained high quality tissue samples under controlled conditions that reveal a level of cellular detail so far inaccessible in human studies. Methods and Results: We performed single nucleus RNA-sequencing in 78 samples in 10 distinct regions including the four chambers of the heart, ventricular septum, sinoatrial node, atrioventricular node, aorta, pulmonary artery, and pulmonary veins (PV), which produced an aggregate map of 505,835 nuclei. We identified 26 distinct cell types and additional subtypes, including a number of rare cell types such as PV cardiomyocytes and non-myelinating Schwann cells (NMSCs), and unique groups of vascular smooth muscle cells (VSMCs), endothelial cells (ECs) and fibroblasts (FBs), which gave rise to a detailed cell type distribution across tissues. We demonstrated differences in the cellular composition across different cardiac regions and tissue-specific differences in transcription for each cell type, highlighting the molecular diversity and complex tissue architecture of the cardiovascular system. Specifically, we observed great transcriptional heterogeneities among ECs and FBs. Importantly, several cell subtypes had a unique regional localization such as a subtype of VSMCs enriched in the large vasculature. We found the cellular makeup of PV tissue is closer to heart tissue than to the large arteries. We further explored the ligand-receptor repertoire across cell clusters and tissues, and observed tissue-enriched cellular communication networks, including heightened Nppa - Npr1/2/3 signaling in the sinoatrial node. Conclusions: Through a large single nucleus sequencing effort encompassing over 500,000 nuclei, we broadened our understanding of cellular transcription in the healthy cardiovascular system. The existence of tissue-restricted cellular phenotypes suggests regional regulation of cardiovascular physiology. The overall conservation in gene expression and molecular pathways across rat and human cell types, together with our detailed transcriptional characterization of each cell type, offers the potential to identify novel therapeutic targets and improve preclinical models of cardiovascular disease.

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

BACKGROUND: ZFHX3 (zinc finger homeobox 3), a gene that encodes a large transcription factor, is at the second-most significantly associated locus with atrial fibrillation (AF), but its function in the heart is unknown. This study aims to identify causative genetic variation related to AF at the ZFHX3 locus and examine the impact of Zfhx3 loss on cardiac function in mice. METHODS: CRISPR-Cas9 genome editing, chromatin immunoprecipitation, and luciferase assays in pluripotent stem cell-derived cardiomyocytes were used to identify causative genetic variation related to AF at the ZFHX3 locus. Cardiac function was assessed by echocardiography, magnetic resonance imaging, electrophysiology studies, calcium imaging, and RNA sequencing in mice with heterozygous and homozygous cardiomyocyte-restricted Zfhx3 loss (Zfhx3 Het and knockout, respectively). Human cardiac single-nucleus ATAC (assay for transposase-accessible chromatin)-sequencing data was analyzed to determine which genes in atrial cardiomyocytes are directly regulated by ZFHX3. RESULTS: We found single-nucleotide polymorphism (SNP) rs12931021 modulates an enhancer regulating ZFHX3 expression, and the AF risk allele is associated with decreased ZFHX3 transcription. We observed a gene-dose response in AF susceptibility with Zfhx3 knockout mice having higher incidence, frequency, and burden of AF than Zfhx3 Het and wild-type mice, with alterations in conduction velocity, atrial action potential duration, calcium handling and the development of atrial enlargement and thrombus, and dilated cardiomyopathy. Zfhx3 loss results in atrial-specific differential effects on genes and signaling pathways involved in cardiac pathophysiology and AF. CONCLUSIONS: Our findings implicate ZFHX3 as the causative gene at the 16q22 locus for AF, and cardiac abnormalities caused by loss of cardiac Zfhx3 are due to atrial-specific dysregulation of pathways involved in AF susceptibility. Together, these data reveal a novel and important role for Zfhx3 in the control of cardiac genes and signaling pathways essential for normal atrial function.

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

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

How Does COVID-19 Affect the Heart?

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

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

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

Transcriptome variation in human tissues revealed by long-read sequencing.

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

A Beary Good Genome: Haplotype-Resolved, Chromosome-Level Assembly of the Brown Bear (Ursus arctos).

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

Single-nucleus profiling of human dilated and hypertrophic cardiomyopathy.

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

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

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

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

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

Deep learning enables genetic analysis of the human thoracic aorta.

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

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

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

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

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

COVID-19 and Cardiovascular Disease: From Bench to Bedside.

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

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

BACKGROUND: Atrial fibrillation (AF) often arises from structural abnormalities in the left atria (LA). Annotation of the noncoding genome in human LA is limited, as are effects on gene expression and chromatin architecture. Many AF-associated genetic variants reside in noncoding regions; this knowledge gap impairs efforts to understand the molecular mechanisms of AF and cardiac conduction phenotypes. METHODS: We generated a model of the LA noncoding genome by profiling 7 histone post-translational modifications (active: H3K4me3, H3K4me2, H3K4me1, H3K27ac, H3K36me3; repressive: H3K27me3, H3K9me3), CTCF binding, and gene expression in samples from 5 individuals without structural heart disease or AF. We used MACS2 to identify peak regions (P<0.01), applied a Markov model to classify regulatory elements, and annotated this model with matched gene expression data. We intersected chromatin states with expression quantitative trait locus, DNA methylation, and HiC chromatin interaction data from LA and left ventricle. Finally, we integrated genome-wide association data for AF and electrocardiographic traits to link disease-related variants to genes. RESULTS: Our model identified 21 epigenetic states, encompassing regulatory motifs, such as promoters, enhancers, and repressed regions. Genes were regulated by proximal chromatin states; repressive states were associated with a significant reduction in gene expression (P<2×10-16). Chromatin states were differentially methylated, promoters were less methylated than repressed regions (P<2×10-16). We identified over 15 000 LA-specific enhancers, defined by homeobox family motifs, and annotated several cardiovascular disease susceptibility loci. Intersecting AF and PR genome-wide association studies loci with long-range chromatin conformation data identified a gene interaction network dominated by NKX2-5, TBX3, ZFHX3, and SYNPO2L. CONCLUSIONS: Profiling the noncoding genome provides new insights into the gene expression and chromatin regulation in human LA tissue. These findings enabled identification of a gene network underlying AF; our experimental and analytic approach can be extended to identify molecular mechanisms for other cardiac diseases and traits.

Epigenetic and Transcriptional Networks Underlying Atrial Fibrillation.

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

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

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

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

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