Search

Background: The pathogenic missense variant p.G125R in TBX5 (T-box transcription factor 5) causes Holt–Oram syndrome (also known as hand–heart syndrome) and early onset of atrial fibrillation. Revealing how an altered key developmental transcription factor modulates cardiac physiology in vivo will provide unique insights into the mechanisms underlying atrial fibrillation in these patients. Methods: We analyzed ECGs of an extended family pedigree of Holt–Oram syndrome patients. Next, we introduced the TBX5-p.G125R variant in the mouse genome ( Tbx5 G125R ) and performed electrophysiologic analyses (ECG, optical mapping, patch clamp, intracellular calcium measurements), transcriptomics (single-nuclei and tissue RNA sequencing), and epigenetic profiling (assay for transposase-accessible chromatin using sequencing, H3K27ac [histone H3 lysine 27 acetylation] CUT&RUN [cleavage under targets and release under nuclease sequencing]). Results: We discovered high incidence of atrial extra systoles and atrioventricular conduction disturbances in Holt–Oram syndrome patients. Tbx5 G125R/+ mice were morphologically unaffected and displayed variable RR intervals, atrial extra systoles, and susceptibility to atrial fibrillation, reminiscent of TBX5-p.G125R patients. Atrial conduction velocity was not affected but systolic and diastolic intracellular calcium concentrations were decreased and action potentials were prolonged in isolated cardiomyocytes of Tbx5 G125R/+ mice compared with controls. Transcriptional profiling of atria revealed the most profound transcriptional changes in cardiomyocytes versus other cell types, and identified over a thousand coding and noncoding transcripts that were differentially expressed. Epigenetic profiling uncovered thousands of TBX5-p.G125R-sensitive, putative regulatory elements (including enhancers) that gained accessibility in atrial cardiomyocytes. The majority of sites with increased accessibility were occupied by Tbx5. The small group of sites with reduced accessibility was enriched for DNA-binding motifs of members of the SP (specificity protein) and KLF (Krüppel-like factor) families of transcription factors. These data show that Tbx5-p.G125R induces changes in regulatory element activity, alters transcriptional regulation, and changes cardiomyocyte behavior, possibly caused by altered DNA binding and cooperativity properties. Conclusions: Our data reveal that a disease-causing missense variant in TBX5 induces profound changes in the atrial transcriptional regulatory network and epigenetic state in vivo, leading to arrhythmia reminiscent of those seen in human TBX5-p.G125R variant carriers.

Rationale: Atrial fibrillation (AF) is the most common cardiac arrhythmia diagnosed in clinical practice. Genome-wide association studies have identified AF-associated common variants across 100+ genomic loci, but the mechanism underlying the impact of these variant loci on AF susceptibility in vivo has remained largely undefined. One such variant region, highly associated with AF, is found at 1q24, close to PRRX1 , encoding the paired-related homeobox 1 transcription factor. Objective: To identify the mechanistic link between the variant region at 1q24 and AF predisposition. Methods and Results: The mouse orthologue of the noncoding variant genomic region ( R1A ) at 1q24 was deleted using CRISPR genome editing. Among the genes sharing the topologically associated domain with the deleted R1A region ( Kifap3 , Prrx1 , Fmo2 , and Prrc2c ), only the broadly expressed gene Prrx1 was downregulated in mutants, and only in cardiomyocytes. Expression and epigenetic profiling revealed that a cardiomyocyte lineage-specific gene program ( Mhrt , Myh6 , Rbm20 , Tnnt2 , Ttn , and Ckm ) was upregulated in R1A −/− atrial cardiomyocytes and that Mef2 (myocyte enhancer factor)-binding motifs were significantly enriched at differentially accessible chromatin sites. Consistently, Prrx1 suppressed Mef2-activated enhancer activity in HL-1 cells. Mice heterozygous or homozygous for the R1A deletion were susceptible to atrial arrhythmia induction, had atrial conduction slowing and more irregular RR intervals. Isolated R1A −/− mouse left atrial cardiomyocytes showed lower action potential upstroke velocities and sodium current, as well as increased systolic and diastolic calcium concentrations compared with controls. Conclusions: The noncoding AF variant region at 1q24 modulates Prrx1 expression in cardiomyocytes. Cardiomyocyte-specific reduction of Prrx1 expression upon deletion of the noncoding region leads to a profound induction of a cardiac lineage-specific gene program and to propensity for AF. These data indicate that AF-associated variants in humans may exert AF predisposition through reduced PRRX1 expression in cardiomyocytes.

Genome-wide association studies have identified noncoding variants near TBX3 that are associated with PR interval and QRS duration, suggesting that subtle changes in TBX3 expression affect atrioventricular conduction system function. To explore whether and to what extent the atrioventricular conduction system is affected by Tbx3 dose reduction, we first characterized electrophysiological properties and morphology of heterozygous Tbx3 mutant (Tbx3(+/−)) mouse hearts. We found PR interval shortening and prolonged QRS duration, as well as atrioventricular bundle hypoplasia after birth in heterozygous mice. The atrioventricular node size was unaffected. Transcriptomic analysis of atrioventricular nodes isolated by laser capture microdissection revealed hundreds of deregulated genes in Tbx3(+/−) mutants. Notably, Tbx3(+/−) atrioventricular nodes showed increased expression of working myocardial gene programs (mitochondrial and metabolic processes, muscle contractility) and reduced expression of pacemaker gene programs (neuronal, Wnt signaling, calcium/ion channel activity). By integrating chromatin accessibility profiles (ATAC sequencing) of atrioventricular tissue and other epigenetic data, we identified Tbx3-dependent atrioventricular regulatory DNA elements (REs) on a genome-wide scale. We used transgenic reporter assays to determine the functionality of candidate REs near Ryr2, an up-regulated chamber-enriched gene, and in Cacna1g, a down-regulated conduction system-specific gene. Using genome editing to delete candidate REs, we showed that a strong intronic bipartite RE selectively governs Cacna1g expression in the conduction system in vivo. Our data provide insights into the multifactorial Tbx3-dependent transcriptional network that regulates the structure and function of the cardiac conduction system, which may underlie the differences in PR duration and QRS interval between individuals carrying variants in the TBX3 locus.

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.

Heart development and rhythm control are highly Tbx5 dosage-sensitive. TBX5 haploinsufficiency causes congenital conduction disorders, whereas increased expression levels of TBX5 in human heart samples has been associated with atrial fibrillation. We deleted the conserved mouse orthologues of two independent AF-associated genomic regions in the Tbx5 locus, one intronic (RE(int)-/-) and one downstream of Tbx5 (RE(down)-/-). In both lines we observed a modest (30%) increase of Tbx5 in the postnatal atria. To gain insight into the effects of slight dosage increase in vivo, we investigated the atrial transcriptional, epigenetic and electrophysiological properties of both lines. We observed induction of genes involved in development, ion transport and conduction, increased action potential duration and increased susceptibility to atrial arrhythmias. We identified an AF-associated variant in the human intronic regulatory region that increases transcriptional activity. Expression of the AF-associated transcription factor Prrx1 was induced in RE(int)-/- cardiomyocytes. We found that some of the transcriptional and functional changes in the atria caused by increased Tbx5 expression were normalized when reducing cardiac Prrx1 expression in RE(int)-/- mice, indicating an interaction between these two AF genes. We conclude that modest increases in expression of dose-dependent transcription factors, caused by common regulatory variants, significantly impact on the cardiac gene regulatory network and disease susceptibility.

Disease-associated genetic variants that lie in non-coding regions found by genome-wide association studies are thought to alter the functionality of transcription regulatory elements and target gene expression. To uncover causal genetic variants, variant regulatory elements and their target genes, here we cross-reference human transcriptomic, epigenomic and chromatin conformation datasets. Of 104 genetic variant regions associated with atrial fibrillation candidate target genes are prioritized. We optimize EMERGE enhancer prediction and use accessible chromatin profiles of human atrial cardiomyocytes to more accurately predict cardiac regulatory elements and identify hundreds of sub-threshold variants that co-localize with regulatory elements. Removal of mouse homologues of atrial fibrillation-associated regions in vivo uncovers a distal regulatory region involved in Gja1 (Cx43) expression. Our analyses provide a shortlist of genes likely affected by atrial fibrillation-associated variants and provide variant regulatory elements in each region that link genetic variation and target gene regulation, helping to focus future investigations.

PURPOSE: Rare genetic variants in KDR, encoding the vascular endothelial growth factor receptor 2 (VEGFR2), have been reported in patients with tetralogy of Fallot (TOF). However, their role in disease causality and pathogenesis remains unclear. METHODS: We conducted exome sequencing in a familial case of TOF and large-scale genetic studies, including burden testing, in >1,500 patients with TOF. We studied gene-targeted mice and conducted cell-based assays to explore the role of KDR genetic variation in the etiology of TOF. RESULTS: Exome sequencing in a family with two siblings affected by TOF revealed biallelic missense variants in KDR. Studies in knock-in mice and in HEK 293T cells identified embryonic lethality for one variant when occurring in the homozygous state, and a significantly reduced VEGFR2 phosphorylation for both variants. Rare variant burden analysis conducted in a set of 1,569 patients of European descent with TOF identified a 46-fold enrichment of protein-truncating variants (PTVs) in TOF cases compared to controls (P = 7 × 10-11). CONCLUSION: Rare KDR variants, in particular PTVs, strongly associate with TOF, likely in the setting of different inheritance patterns. Supported by genetic and in vivo and in vitro functional analysis, we propose loss-of-function of VEGFR2 as one of the mechanisms involved in the pathogenesis of TOF. ispartof: GENETICS IN MEDICINE vol:23 issue:10 pages:1952-1960 ispartof: location:United States status: published

[Figure: see text].

Genetic variants inThe expression ofWe found that cardiomyocytes of the atria, sinoatrial node, and ventricular conduction system express a short transcript comprising the last 7 exons of the gene (Genetic variants in and around

Background: Genetic variants in SCN10A , encoding the neuronal voltage-gated sodium channel Na V 1.8, are strongly associated with atrial fibrillation, Brugada syndrome, cardiac conduction velocities, and heart rate. The cardiac function of SCN10A has not been resolved, however, and diverging mechanisms have been proposed. Here, we investigated the cardiac expression of SCN10A and the function of a variant-sensitive intronic enhancer previously linked to the regulation of SCN5A , encoding the major essential cardiac sodium channel Na V 1.5. Methods: The expression of SCN10A was investigated in mouse and human hearts. With the use of CRISPR/Cas9 genome editing, the mouse intronic enhancer was disrupted, and mutant mice were characterized by transcriptomic and electrophysiological analyses. The association of genetic variants at SCN5A-SCN10A enhancer regions and gene expression were evaluated by genome-wide association studies single-nucleotide polymorphism mapping and expression quantitative trait loci analysis. Results: We found that cardiomyocytes of the atria, sinoatrial node, and ventricular conduction system express a short transcript comprising the last 7 exons of the gene ( Scn10a-short ). Transcription occurs from an intronic enhancer-promoter complex, whereas full-length Scn10a transcript was undetectable in the human and mouse heart. Expression quantitative trait loci analysis revealed that the genetic variants in linkage disequilibrium with genetic variant rs6801957 in the intronic enhancer associate with SCN10A transcript levels in the heart. Genetic modification of the enhancer in the mouse genome led to reduced cardiac Scn10a-short expression in atria and ventricles, reduced cardiac sodium current in atrial cardiomyocytes, atrial conduction slowing and arrhythmia, whereas the expression of Scn5a , the presumed enhancer target gene, remained unaffected. In patch-clamp transfection experiments, expression of Scn10a-short –encoded Na V 1.8-short increased Na V 1.5-mediated sodium current. We propose that noncoding genetic variation modulates transcriptional regulation of Scn10a-short in cardiomyocytes that impacts Na V 1.5-mediated sodium current and heart rhythm. Conclusions: Genetic variants in and around SCN10A modulate enhancer function and expression of a cardiac-specific SCN10A-short transcript. We propose that noncoding genetic variation modulates transcriptional regulation of a functional C-terminal portion of Na V 1.8 in cardiomyocytes that impacts on Na V 1.5 function, cardiac conduction velocities, and arrhythmia susceptibility.

Rationale: The development and function of the pacemaker cardiomyocytes of the sinoatrial node (SAN), the leading pacemaker of the heart, are tightly controlled by a conserved network of transcription factors, including TBX3 (T-box transcription factor 3), ISL1 (ISL LIM homeobox 1), and SHOX2 (short stature homeobox 2). Yet, the regulatory DNA elements (REs) controlling target gene expression in the SAN pacemaker cells have remained undefined. Objective: Identification of the regulatory landscape of human SAN-like pacemaker cells and functional assessment of SAN-specific REs potentially involved in pacemaker cell gene regulation. Methods and Results: We performed Assay for Transposase-Accessible Chromatin using sequencing on human pluripotent stem cell–derived SAN-like pacemaker cells and ventricle-like cells and identified thousands of putative REs specific for either human cell type. We validated pacemaker cell–specific elements in the SHOX2 and TBX3 loci. CRISPR-mediated homozygous deletion of the mouse ortholog of a noncoding region with candidate pacemaker-specific REs in the SHOX2 locus resulted in selective loss of Shox2 expression from the developing SAN and embryonic lethality. Putative pacemaker-specific REs were identified up to 1 Mbp upstream of TBX3 in a region close to MED13L harboring variants associated with heart rate recovery after exercise. The orthologous region was deleted in mice, which resulted in selective loss of expression of Tbx3 from the SAN and (cardiac) ganglia and in neonatal lethality. Expression of Tbx3 was maintained in other tissues including the atrioventricular conduction system, lungs, and liver. Heterozygous adult mice showed increased SAN recovery times after pacing. The human REs harboring the associated variants robustly drove expression in the SAN of transgenic mouse embryos. Conclusions: We provided a genome-wide collection of candidate human pacemaker-specific REs, including the loci of SHOX2 , TBX3 , and ISL1 , and identified a link between human genetic variants influencing heart rate recovery after exercise and a variant RE with highly conserved function, driving SAN expression of TBX3 .

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.

Purpose: Rare genetic variants in KDR, encoding the vascular endothelial growth factor receptor 2 (VEGFR2), have been reported in patients with tetralogy of Fallot (TOF). However, their role in disease causality and pathogenesis remains unclear. Methods: We conducted exome sequencing in a familial case of TOF and large-scale genetic studies, including burden testing, in >1,500 patients with TOF. We studied gene-targeted mice and conducted cell-based assays to explore the role of KDR genetic variation in the etiology of TOF. Results: Exome sequencing in a family with two siblings affected by TOF revealed biallelic missense variants in KDR. Studies in knock-in mice and in HEK 293T cells identified embryonic lethality for one variant when occurring in the homozygous state, and a significantly reduced VEGFR2 phosphorylation for both variants. Rare variant burden analysis conducted in a set of 1,569 patients of European descent with TOF identified a 46-fold enrichment of protein-truncating variants (PTVs) in TOF cases compared to controls (P = 7 × 10-11). Conclusion: Rare KDR variants, in particular PTVs, strongly associate with TOF, likely in the setting of different inheritance patterns. Supported by genetic and in vivo and in vitro functional analysis, we propose loss-of-function of VEGFR2 as one of the mechanisms involved in the pathogenesis of TOF.

Heart valve development proceeds through coordinated steps by which endocardial cushions (ECs) form thin, elongated, and stratified valves. Wnt signaling and its canonical effector β-catenin are proposed to contribute to endocardial-to-mesenchymal transformation (EMT) through postnatal steps of valvulogenesis. However, genetic redundancy and lethality have made it challenging to define specific roles of the canonical Wnt pathway at different stages of valve formation. We developed a transgenic mouse system that provides spatiotemporal inhibition of Wnt/β-catenin signaling by chemically-inducible overexpression of Dkk1. Unexpectedly, this approach indicates canonical Wnt signaling is required for EMT in the proximal outflow tract (pOFT) but not atrioventricular canal (AVC) cushions. Further, Wnt indirectly promotes pOFT EMT through its earlier activity in neighboring myocardial cells or their progenitors. Subsequently, Wnt/β-catenin signaling is activated in cushion mesenchymal cells where it supports FGF-driven expansion of ECs and then AVC valve extracellular matrix patterning. Mice lacking Axin2, a negative Wnt regulator, have larger valves, suggesting that accumulating Axin2 in maturing valves represents negative feedback that restrains tissue overgrowth rather than simply reporting Wnt activity. Disruption of these Wnt/β-catenin signaling roles that enable developmental transitions during valvulogenesis could account for common congenital valve defects.