Search

The cardiac conduction system (CCS) orchestrates the electrical impulses that enable coordinated contraction of the cardiac chambers. The T-box transcription factors TBX3 and TBX5 are required for cardiac conduction system development and associated with overlapping and distinct human cardiac conduction system diseases. We evaluated the coordinated role of Tbx3 and Tbx5 in the murine ventricular conduction system (VCS). We engineered a compound Tbx3:Tbx5 conditional knockout allele for both genes located in cis on mouse chromosome 5. Conditional deletion of both T-box transcriptional factors in the ventricular conduction system, using the VCS-specific Mink:Cre, caused loss of VCS function and molecular identity. Combined Tbx3 and Tbx5 deficiency in the adult VCS led to conduction defects, including prolonged PR and QRS intervals and elevated susceptibility to ventricular tachycardia. These electrophysiologic defects occurred prior to detectable alterations in cardiac contractility or histologic morphology, indicative of a primary conduction system defect. Tbx3:Tbx5 double knockout VCS cardiomyocytes revealed a transcriptional shift towards non-CCS-specialized working myocardium, suggesting reprogramming of their cellular identity. Furthermore, optical mapping revealed a loss of VCS-specific conduction system propagation. Collectively, these findings indicate that Tbx3 and Tbx5 coordinate to control VCS molecular fate and function, with implications for understanding cardiac conduction disorders in humans., Competing Interests: Conflict of Interest Statement None declared.

Nuclear envelope (NE) ruptures are emerging observations in Lamin-related dilated cardiomyopathy, an adult-onset disease caused by loss-of-function mutations in Lamin A/C, a nuclear lamina component. Here, we test a prevailing hypothesis that NE ruptures trigger the pathological cGAS-STING cytosolic DNA-sensing pathway using a mouse model of Lamin cardiomyopathy. The reduction of Lamin A/C in cardio-myocyte of adult mice causes pervasive NE ruptures in cardiomyocytes, preceding inflammatory transcription, fibrosis, and fatal dilated cardiomyopathy. NE ruptures are followed by DNA damage accumulation without causing immediate cardiomyocyte death. However, cGAS-STING-dependent inflammatory signaling remains inactive. Deleting cGas or Sting does not rescue cardiomyopathy in the mouse model. The lack of cGAS-STING activation is likely due to the near absence of cGAS expression in adult cardiomyocytes at baseline. Instead, extracellular matrix (ECM) signaling is activated and predicted to initiate pro-inflammatory communication from Lamin-reduced cardiomyocytes to fibroblasts. Our work nominates ECM signaling, not cGAS-STING, as a potential inflammatory contributor in Lamin cardiomyopathy., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2 ) is observed in human ventricles during HF and causes HF in mice., Methods: FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription., Results: FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency., Conclusions: Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling., Competing Interests: Disclosures None.

Background: Conotruncal defects due to developmental abnormalities of the outflow tract (OFT) are an important cause of cyanotic congenital heart disease. Dysregulation of transcriptional programs tuned by NKX2-5 (NK2 homeobox 5), GATA6 (GATA binding protein 6), and TBX1 (T-box transcription factor 1) have been implicated in abnormal OFT morphogenesis. However, there remains no consensus on how these transcriptional programs function in a unified gene regulatory network within the OFT., Methods: We generated mice harboring a 226-nucleotide deletion of a highly conserved cardiac enhancer containing 2 GATA-binding sites located ≈9.4 kb upstream of the transcription start site of Nkx2-5 ( Nkx2-5 ∆enh ) using CRISPR-Cas9 gene editing and assessed phenotypes. Cardiac defects in Nkx2-5 ∆enh/∆enh mice were structurally characterized using histology and scanning electron microscopy, and physiologically assessed using electrocardiography, echocardiography, and optical mapping. Transcriptome analyses were performed using RNA sequencing and single-cell RNA sequencing data sets. Endogenous GATA6 interaction with and activity on the NKX2-5 enhancer was studied using chromatin immunoprecipitation sequencing and transposase-accessible chromatin sequencing in human induced pluripotent stem cell-derived cardiomyocytes., Results: Nkx2-5 ∆enh/∆enh mice recapitulated cyanotic conotruncal defects seen in patients with NKX2-5, GATA6 , and TBX1 mutations. Nkx2-5 ∆enh/∆enh mice also exhibited defects in right Purkinje fiber network formation, resulting in right bundle-branch block. Enhancer deletion reduced embryonic Nkx2-5 expression selectively in the right ventricle and OFT of mutant hearts, indicating that enhancer activity is localized to the anterior second heart field. Transcriptional profiling of the mutant OFT revealed downregulation of important genes involved in OFT rotation and septation, such as Tbx1, Pitx2 , and Sema3c . Endogenous GATA6 interacted with the highly conserved enhancer in human induced pluripotent stem cell-derived cardiomyocytes and in wild-type mouse hearts. We found critical dose dependency of cardiac enhancer accessibility on GATA6 gene dosage in human induced pluripotent stem cell-derived cardiomyocytes., Conclusions: Our results using human and mouse models reveal an essential gene regulatory network of the OFT that requires an anterior second heart field enhancer to link GATA6 with NKX2-5-dependent rotation and septation gene programs., Competing Interests: Disclosures None.

Understanding how the atrial and ventricular heart chambers maintain distinct identities is a prerequisite for treating chamber-specific diseases. Here, we selectively knocked out (KO) the transcription factor Tbx5 in the atrial working myocardium to evaluate its requirement for atrial identity. Atrial Tbx5 inactivation downregulated atrial cardiomyocyte (aCM) selective gene expression. Using concurrent single nucleus transcriptome and open chromatin profiling, genomic accessibility differences were identified between control and Tbx5 KO aCMs, revealing that 69% of the control-enriched ATAC regions were bound by TBX5. Genes associated with these regions were downregulated in KO aCMs, suggesting they function as TBX5-dependent enhancers. Comparing enhancer chromatin looping using H3K27ac HiChIP identified 510 chromatin loops sensitive to TBX5 dosage, and 74.8% of control-enriched loops contained anchors in control-enriched ATAC regions. Together, these data demonstrate TBX5 maintains the atrial gene expression program by binding to and preserving the tissue-specific chromatin architecture of atrial enhancers., Competing Interests: COMPETING INTERESTS The authors have no competing interests to declare.

Genome-wide association studies (GWAS) have linked hundreds of loci to cardiac diseases. However, in most loci the causal variants and their target genes remain unknown. We developed a combined experimental and analytical approach that integrates single cell epigenomics with GWAS to prioritize risk variants and genes. We profiled accessible chromatin in single cells obtained from human hearts and leveraged the data to study genetics of Atrial Fibrillation (AF), the most common cardiac arrhythmia. Enrichment analysis of AF risk variants using cell-type-resolved open chromatin regions (OCRs) implicated cardiomyocytes as the main mediator of AF risk. We then performed statistical fine-mapping, leveraging the information in OCRs, and identified putative causal variants in 122 AF-associated loci. Taking advantage of the fine-mapping results, our novel statistical procedure for gene discovery prioritized 46 high-confidence risk genes, highlighting transcription factors and signal transduction pathways important for heart development. In summary, our analysis provides a comprehensive map of AF risk variants and genes, and a general framework to integrate single-cell genomics with genetic studies of complex traits., (© 2023. Springer Nature Limited.)

Mechanisms underlying distinct specification, commitment, and differentiation phases of cell fate determination remain undefined due to difficulties capturing these processes. Here, we interrogate the activity of ETV2, a transcription factor necessary and sufficient for hematoendothelial differentiation, within isolated fate intermediates. We observe transcriptional upregulation of Etv2 and opening of ETV2-binding sites, indicating new ETV2 binding, in a common cardiac-hematoendothelial progenitor population. Accessible ETV2-binding sites are active at the Etv2 locus but not at other hematoendothelial regulator genes. Hematoendothelial commitment coincides with the activation of a small repertoire of previously accessible ETV2-binding sites at hematoendothelial regulators. Hematoendothelial differentiation accompanies activation of a large repertoire of new ETV2-binding sites and upregulation of hematopoietic and endothelial gene regulatory networks. This work distinguishes specification, commitment, and sublineage differentiation phases of ETV2-dependent transcription and suggests that the shift from ETV2 binding to ETV2-bound enhancer activation, not ETV2 binding to target enhancers, drives hematoendothelial fate commitment., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Frem1 has been linked to human face shape variation, dysmorphology, and malformation, but little is known about its regulation and biological role in facial development., Results: During midfacial morphogenesis in mice, we observed Frem1 expression in the embryonic growth centers that form the median upper lip, nose, and palate. Expansive spatial gradients of Frem1 expression in the cranial neural crest cell (cNCC) mesenchyme of these tissues suggested transcriptional regulation by a secreted morphogen. Accordingly, Frem1 expression paralleled that of the conserved Sonic Hedgehog (Shh) target gene Gli1 in the cNCC mesenchyme. Suggesting direct transcriptional regulation by Shh signaling, we found that Frem1 expression is induced by SHH ligand stimulation or downstream pathway activation in cNCCs and observed GLI transcription factor binding at the Frem1 transcriptional start site during midfacial morphogenesis. Finally, we found that FREM1 is sufficient to induce cNCC proliferation in a concentration-dependent manner and that Shh pathway antagonism reduces Frem1 expression during pathogenesis of midfacial hypoplasia., Conclusions: By demonstrating that the Shh signaling pathway regulates Frem1 expression in cNCCs, these findings provide novel insight into the mechanisms underlying variation in midfacial morphogenesis., (© 2022 The Authors. Developmental Dynamics published by Wiley Periodicals LLC on behalf of American Association for Anatomy.)

Background: Brugada syndrome (BrS) is an inherited arrhythmia syndrome caused by loss-of-function variants in the cardiac sodium channel gene SCN5A (sodium voltage-gated channel alpha subunit 5) in ≈20% of subjects. We identified a family with 4 individuals diagnosed with BrS harboring the rare G145R missense variant in the cardiac transcription factor TBX5 (T-box transcription factor 5) and no SCN5A variant., Methods: We generated induced pluripotent stem cells (iPSCs) from 2 members of a family carrying TBX5-G145R and diagnosed with Brugada syndrome. After differentiation to iPSC-derived cardiomyocytes (iPSC-CMs), electrophysiologic characteristics were assessed by voltage- and current-clamp experiments (n=9 to 21 cells per group) and transcriptional differences by RNA sequencing (n=3 samples per group), and compared with iPSC-CMs in which G145R was corrected by CRISPR/Cas9 approaches. The role of platelet-derived growth factor (PDGF)/phosphoinositide 3-kinase (PI3K) pathway was elucidated by small molecule perturbation. The rate-corrected QT (QTc) interval association with serum PDGF was tested in the Framingham Heart Study cohort (n=1893 individuals)., Results: TBX5-G145R reduced transcriptional activity and caused multiple electrophysiologic abnormalities, including decreased peak and enhanced "late" cardiac sodium current (I Na ), which were entirely corrected by editing G145R to wild-type. Transcriptional profiling and functional assays in genome-unedited and -edited iPSC-CMs showed direct SCN5A down-regulation caused decreased peak I Na , and that reduced PDGF receptor ( PDGFRA [platelet-derived growth factor receptor α]) expression and blunted signal transduction to PI3K was implicated in enhanced late I Na . Tbx5 regulation of the PDGF axis increased arrhythmia risk due to disruption of PDGF signaling and was conserved in murine model systems. PDGF receptor blockade markedly prolonged normal iPSC-CM action potentials and plasma levels of PDGF in the Framingham Heart Study were inversely correlated with the QTc interval ( P <0.001)., Conclusions: These results not only establish decreased SCN5A transcription by the TBX5 variant as a cause of BrS, but also reveal a new general transcriptional mechanism of arrhythmogenesis of enhanced late sodium current caused by reduced PDGF receptor-mediated PI3K signaling.

Many developmental signaling pathways have been implicated in lineage-specific differentiation; however, mechanisms that explicitly control differentiation timing remain poorly defined in mammals. We report that murine Hedgehog signaling is a heterochronic pathway that determines the timing of progenitor differentiation. Hedgehog activity was necessary to prevent premature differentiation of second heart field (SHF) cardiac progenitors in mouse embryos, and the Hedgehog transcription factor GLI1 was sufficient to delay differentiation of cardiac progenitors in vitro. GLI1 directly activated a de novo progenitor-specific network in vitro, akin to that of SHF progenitors in vivo, which prevented the onset of the cardiac differentiation program. A Hedgehog signaling-dependent active-to-repressive GLI transition functioned as a differentiation timer, restricting the progenitor network to the SHF. GLI1 expression was associated with progenitor status across germ layers, and it delayed the differentiation of neural progenitors in vitro, suggesting a broad role for Hedgehog signaling as a heterochronic pathway., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

Analyzing single-cell transcriptomes promises to decipher the plasticity, heterogeneity, and rapid switches in developmental cellular state transitions. Such analyses require the identification of gene markers for semi-stable transition states. However, there are nontrivial challenges such as unexplainable stochasticity, variable population sizes, and alternative trajectory constructions. By advancing current tipping-point theory-based models with feature selection, network decomposition, accurate estimation of correlations, and optimization, we developed BioTIP to overcome these challenges. BioTIP identifies a small group of genes, called critical transition signal (CTS), to characterize regulated stochasticity during semi-stable transitions. Although methods rooted in different theories converged at the same transition events in two benchmark datasets, BioTIP is unique in inferring lineage-determining transcription factors governing critical transition. Applying BioTIP to mouse gastrulation data, we identify multiple CTSs from one dataset and validated their significance in another independent dataset. We detect the established regulator Etv2 whose expression change drives the haemato-endothelial bifurcation, and its targets together in CTS across three datasets. After comparing to three current methods using six datasets, we show that BioTIP is accurate, user-friendly, independent of pseudo-temporal trajectory, and captures significantly interconnected and reproducible CTSs. We expect BioTIP to provide great insight into dynamic regulations of lineage-determining factors., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Endothelial and erythropoietic lineages arise from a common developmental progenitor. Etv2 is a master transcriptional regulator required for the development of both lineages. However, the mechanisms through which Etv2 initiates the gene-regulatory networks (GRNs) for endothelial and erythropoietic specification and how the two GRNs diverge downstream of Etv2 remain incompletely understood. Here, by analyzing a hypomorphic Etv2 mutant, we demonstrate different threshold requirements for initiation of the downstream GRNs for endothelial and erythropoietic development. We show that Etv2 functions directly in a coherent feedforward transcriptional network for vascular endothelial development, and a low level of Etv2 expression is sufficient to induce and sustain the endothelial GRN. In contrast, Etv2 induces the erythropoietic GRN indirectly via activation of Tal1, which requires a significantly higher threshold of Etv2 to initiate and sustain erythropoietic development. These results provide important mechanistic insight into the divergence of the endothelial and erythropoietic lineages., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

The transcription factor NFATC2 induces β cell proliferation in mouse and human islets. However, the genomic targets that mediate these effects have not been identified. We expressed active forms of Nfatc2 and Nfatc1 in human islets. By integrating changes in gene expression with genomic binding sites for NFATC2, we identified approximately 2200 transcriptional targets of NFATC2. Genes induced by NFATC2 were enriched for transcripts that regulate the cell cycle and for DNA motifs associated with the transcription factor FOXP. Islets from an endocrine-specific Foxp1, Foxp2, and Foxp4 triple-knockout mouse were less responsive to NFATC2-induced β cell proliferation, suggesting the FOXP family works to regulate β cell proliferation in concert with NFATC2. NFATC2 induced β cell proliferation in both mouse and human islets, whereas NFATC1 did so only in human islets. Exploiting this species difference, we identified approximately 250 direct transcriptional targets of NFAT in human islets. This gene set enriches for cell cycle-associated transcripts and includes Nr4a1. Deletion of Nr4a1 reduced the capacity of NFATC2 to induce β cell proliferation, suggesting that much of the effect of NFATC2 occurs through its induction of Nr4a1. Integration of noncoding RNA expression, chromatin accessibility, and NFATC2 binding sites enabled us to identify NFATC2-dependent enhancer loci that mediate β cell proliferation.

The gene regulatory networks that coordinate the development of the cardiac and pulmonary systems are essential for terrestrial life but poorly understood. The T-box transcription factor Tbx5 is critical for both pulmonary specification and heart development, but how these activities are mechanistically integrated remains unclear. Here using Xenopus and mouse embryos, we establish molecular links between Tbx5 and retinoic acid (RA) signaling in the mesoderm and between RA signaling and sonic hedgehog expression in the endoderm to unveil a conserved RA-Hedgehog-Wnt signaling cascade coordinating cardiopulmonary (CP) development. We demonstrate that Tbx5 directly maintains expression of aldh1a2, the RA-synthesizing enzyme, in the foregut lateral plate mesoderm via an evolutionarily conserved intronic enhancer. Tbx5 promotes posterior second heart field identity in a positive feedback loop with RA, antagonizing a Fgf8-Cyp regulatory module to restrict FGF activity to the anterior. We find that Tbx5/Aldh1a2-dependent RA signaling directly activates shh transcription in the adjacent foregut endoderm through a conserved MACS1 enhancer. Hedgehog signaling coordinates with Tbx5 in the mesoderm to activate expression of wnt2/2b, which induces pulmonary fate in the foregut endoderm. These results provide mechanistic insight into the interrelationship between heart and lung development informing CP evolution and birth defects., Competing Interests: SR, JS, XY, AR, KA, PC, KI, MH, IM, AZ No competing interests declared, (© 2021, Rankin et al.)

Congenital Heart Disease (CHD), malformations of the heart present at birth, is the most common class of life-threatening birth defect (Hoffman (1995) [1], Gelb (2004) [2], Gelb (2014) [3]). A major research challenge is to elucidate the genetic determinants of CHD and mechanistically link CHD ontogeny to a molecular understanding of heart development. Although the embryonic origins of CHD are unclear in most cases, dysregulation of cardiovascular lineage specification, patterning, proliferation, migration or differentiation have been described (Olson (2004) [4], Olson (2006) [5], Srivastava (2006) [6], Dunwoodie (2007) [7], Bruneau (2008) [8]). Cardiac differentiation is the process whereby cells become progressively more dedicated in a trajectory through the cardiac lineage towards mature cardiomyocytes. Defects in cardiac differentiation have been linked to CHD, although how the complex control of cardiac differentiation prevents CHD is just beginning to be understood. The stages of cardiac differentiation are highly stereotyped and have been well-characterized (Kattman et al. (2011) [9], Wamstad et al. (2012) [10], Luna-Zurita et al. (2016) [11], Loh et al. (2016) [12], DeLaughter et al. (2016) [13]); however, the developmental and molecular mechanisms that promote or delay the transition of a cell through these stages have not been as deeply investigated. Tight temporal control of progenitor differentiation is critically important for normal organ size, spatial organization, and cellular physiology and homeostasis of all organ systems (Raff et al. (1985) [14], Amthor et al. (1998) [15], Kopan et al. (2014) [16]). This review will focus on the action of signaling pathways in the control of cardiomyocyte differentiation timing. Numerous signaling pathways, including the Wnt, Fibroblast Growth Factor, Hedgehog, Bone Morphogenetic Protein, Insulin-like Growth Factor, Thyroid Hormone and Hippo pathways, have all been implicated in promoting or inhibiting transitions along the cardiac differentiation trajectory. Gaining a deeper understanding of the mechanisms controlling cardiac differentiation timing promises to yield insights into the etiology of CHD and to inform approaches to restore function to damaged hearts., (Copyright © 2021. Published by Elsevier Ltd.)

Atrioventricular septal defects (AVSD) are a severe congenital heart defect present in individuals with Down syndrome (DS) at a > 2000-fold increased prevalence compared to the general population. This study aimed to identify risk-associated genes and pathways and to examine a potential polygenic contribution to AVSD in DS. We analyzed a total cohort of 702 individuals with DS with or without AVSD, with genomic data from whole exome sequencing, whole genome sequencing, and/or array-based imputation. We utilized sequence kernel association testing and polygenic risk score (PRS) methods to examine rare and common variants. Our findings suggest that the Notch pathway, particularly NOTCH4, as well as genes involved in the ciliome including CEP290 may play a role in AVSD in DS. These pathways have also been implicated in DS-associated AVSD in prior studies. A polygenic component for AVSD in DS has not been examined previously. Using weights based on the largest genome-wide association study of congenital heart defects available (2594 cases and 5159 controls; all general population samples), we found PRS to be associated with AVSD with odds ratios ranging from 1.2 to 1.3 per standard deviation increase in PRS and corresponding liability r 2 values of approximately 1%, suggesting at least a small polygenic contribution to DS-associated AVSD. Future studies with larger sample sizes will improve identification and quantification of genetic contributions to AVSD in DS.

Atrial fibrillation (AF) is the most common cardiac arrhythmia, yet the molecular signature of the vulnerable atrial substrate is not well understood. Here, we delineated a distinct transcriptional signature in right versus left atrial cardiomyocytes (CMs) at baseline and identified chamber-specific gene expression changes in patients with a history of AF in the setting of end-stage heart failure (AF+HF) that are not present in heart failure alone (HF). We observed that human left atrial (LA) CMs exhibited Notch pathway activation and increased ploidy in AF+HF but not in HF alone. Transient activation of Notch signaling within adult CMs in a murine genetic model is sufficient to increase ploidy in both atrial chambers. Notch activation within LA CMs generated a transcriptomic fingerprint resembling AF, with dysregulation of transcription factor and ion channel genes, including Pitx2, Tbx5, Kcnh2, Kcnq1, and Kcnip2. Notch activation also produced distinct cellular electrophysiologic responses in LA versus right atrial CMs, prolonging the action potential duration (APD) without altering the upstroke velocity in the left atrium and reducing the maximal upstroke velocity without altering the APD in the right atrium. Our results support a shared human/murine model of increased Notch pathway activity predisposing to AF.

Rationale: The heartbeat is organized by the cardiac conduction system (CCS), a specialized network of cardiomyocytes. Patterning of the CCS into atrial node versus ventricular conduction system (VCS) components with distinct physiology is essential for the normal heartbeat. Distinct node versus VCS physiology has been recognized for more than a century, but the molecular basis of this regional patterning is not well understood., Objective: To study the genetic and genomic mechanisms underlying node versus VCS distinction and investigate rhythm consequences of failed VCS patterning., Methods and Results: Using mouse genetics, we found that the balance between T-box transcriptional activator, Tbx5 , and T-box transcriptional repressor, Tbx3 , determined the molecular and functional output of VCS myocytes. Adult VCS-specific removal of Tbx5 or overexpression of Tbx3 re-patterned the fast VCS into slow, nodal-like cells based on molecular and functional criteria. In these cases, gene expression profiling showed diminished expression of genes required for VCS-specific fast conduction but maintenance of expression of genes required for nodal slow conduction physiology. Action potentials of Tbx5 -deficient VCS myocytes adopted nodal-specific characteristics, including increased action potential duration and cellular automaticity. Removal of Tbx5 in vivo precipitated inappropriate depolarizations in the atrioventricular (His)-bundle associated with lethal ventricular arrhythmias. TBX5 bound and directly activated cis -regulatory elements at fast conduction channel genes required for fast physiological characteristics of the VCS action potential, defining the identity of the adult VCS., Conclusions: The CCS is patterned entirely as a slow, nodal ground state, with a T-box dependent, physiologically dominant, fast conduction network driven specifically in the VCS. Disruption of the fast VCS gene regulatory network allowed nodal physiology to emerge, providing a plausible molecular mechanism for some lethal ventricular arrhythmias.

The mechanisms used by embryos to pattern tissues across their axes has fascinated developmental biologists since the founding of embryology. Here, using single-cell technology, we interrogate complex patterning defects and define a Hedgehog (Hh)-fibroblast growth factor (FGF) signaling axis required for anterior mesoderm lineage development during gastrulation. Single-cell transcriptome analysis of Hh-deficient mesoderm revealed selective deficits in anterior mesoderm populations, culminating in defects to anterior embryonic structures, including the pharyngeal arches, heart, and anterior somites. Transcriptional profiling of Hh-deficient mesoderm during gastrulation revealed disruptions to both transcriptional patterning of the mesoderm and FGF signaling for mesoderm migration. Mesoderm-specific Fgf4 / Fgf8 double-mutants recapitulated anterior mesoderm defects and Hh-dependent GLI transcription factors modulated enhancers at FGF gene loci. Cellular migration defects during gastrulation induced by Hh pathway antagonism were mitigated by the addition of FGF4 protein. These findings implicate a multicomponent signaling hierarchy activated by Hh ligands from the embryonic node and executed by FGF signals in nascent mesoderm to control anterior mesoderm patterning., Competing Interests: The authors declare no competing interest.

Rationale: ZO-1 (Zona occludens 1), encoded by the tight junction protein 1 ( TJP1 ) gene, is a regulator of paracellular permeability in epithelia and endothelia. ZO-1 interacts with the actin cytoskeleton, gap, and adherens junction proteins and localizes to intercalated discs in cardiomyocytes. However, the contribution of ZO-1 to cardiac physiology remains poorly defined., Objective: We aim to determine the role of ZO-1 in cardiac function., Methods and Results: Inducible cardiomyocyte-specific Tjp1 deletion mice ( Tjp1 fl/fl ; Myh6 Cre/Esr1* ) were generated by crossing the Tjp1 floxed mice and Myh6 Cre/Esr1* transgenic mice. Tamoxifen-induced loss of ZO-1 led to atrioventricular (AV) block without changes in heart rate, as measured by ECG and ex vivo optical mapping. Mice with tamoxifen-induced conduction system-specific deletion of Tjp1 ( Tjp1 fl/fl ; Hcn4 CreERt2 ) developed AV block while tamoxifen-induced conduction system deletion of Tjp1 distal to the AV node ( Tjp1 fl/fl ; Kcne1 CreERt2 ) did not demonstrate conduction defects. Western blot and immunostaining analyses of AV nodes showed that ZO-1 loss decreased Cx (connexin) 40 expression and intercalated disc localization. Consistent with the mouse model study, immunohistochemical staining showed that ZO-1 is abundantly expressed in the human AV node and colocalizes with Cx40. Ventricular conduction was not altered despite decreased localization of ZO-1 and Cx43 at the ventricular intercalated disc and modestly decreased left ventricular ejection fraction, suggesting ZO-1 is differentially required for AV node and ventricular conduction., Conclusions: ZO-1 is a key protein responsible for maintaining appropriate AV node conduction through maintaining gap junction protein localization.

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.

LMNA encodes nuclear Lamin A/C that tethers lamina-associated domains (LADs) to the nuclear periphery. Mutations in LMNA cause degenerative disorders including the premature aging disorder Hutchinson-Gilford progeria, but the mechanisms are unknown. We report that Ser22-phosphorylated (pS22) Lamin A/C was localized to the nuclear interior in human fibroblasts throughout the cell cycle. pS22-Lamin A/C interacted with a subset of putative active enhancers, not LADs, at locations co-bound by the transcriptional activator c-Jun. In progeria-patient fibroblasts, a subset of pS22-Lamin A/C-binding sites were lost, whereas new pS22-Lamin A/C-binding sites emerged in normally quiescent loci. New pS22-Lamin A/C binding was accompanied by increased histone acetylation, increased c-Jun binding, and upregulation of nearby genes implicated in progeria pathophysiology. These results suggest that Lamin A/C regulates gene expression by enhancer binding. Disruption of the gene regulatory rather than LAD tethering function of Lamin A/C may underlie the pathogenesis of disorders caused by LMNA mutations., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Long noncoding RNAs (lncRNAs) localize in the cell nucleus and influence gene expression through a variety of molecular mechanisms. Chromatin-enriched RNAs (cheRNAs) are a unique class of lncRNAs that are tightly bound to chromatin and putatively function to locally cis-activate gene transcription. CheRNAs can be identified by biochemical fractionation of nuclear RNA followed by RNA sequencing, but until now, a rigorous analytic pipeline for nuclear RNA-seq has been lacking. In this study, we survey four computational strategies for nuclear RNA-seq data analysis and develop a new pipeline, Tuxedo-ch, which outperforms other approaches. Tuxedo-ch assembles a more complete transcriptome and identifies cheRNA with higher accuracy than other approaches. We used Tuxedo-ch to analyze benchmark datasets of K562 cells and further characterize the genomic features of intergenic cheRNA (icheRNA) and their similarity to enhancer RNAs (eRNAs). We quantify the transcriptional correlation of icheRNA and adjacent genes and show that icheRNA is more positively associated with neighboring gene expression than eRNA or cap analysis of gene expression (CAGE) signals. We also explore two novel genomic associations of cheRNA, which indicate that cheRNAs may function to promote or repress gene expression in a context-dependent manner. IcheRNA loci with significant levels of H3K9me3 modifications are associated with active enhancers, consistent with the hypothesis that enhancers are derived from ancient mobile elements. In contrast, antisense cheRNA (as-cheRNA) may play a role in local gene repression, possibly through local RNA:DNA:DNA triple-helix formation., Competing Interests: The authors declare no competing interests.

Identification of cell type-specific cis- regulatory elements (CREs) is crucial for understanding development and disease, although identification of functional regulatory elements remains challenging. We hypothesized that context-specific CREs could be identified by context-specific non-coding RNA (ncRNA) profiling, based on the observation that active CREs produce ncRNAs. We applied ncRNA profiling to identify rod and cone photoreceptor CREs from wild-type and mutant mouse retinas, defined by presence or absence, respectively, of the rod-specific transcription factor (TF) Nrl Nrl -dependent ncRNA expression strongly correlated with epigenetic profiles of rod and cone photoreceptors, identified thousands of candidate rod- and cone-specific CREs, and identified motifs for rod- and cone-specific TFs. Colocalization of NRL and the retinal TF CRX correlated with rod-specific ncRNA expression, whereas CRX alone favored cone-specific ncRNA expression, providing quantitative evidence that heterotypic TF interactions distinguish cell type-specific CRE activity. We validated the activity of novel Nrl -dependent ncRNA-defined CREs in developing cones. This work supports differential ncRNA profiling as a platform for the identification of cell type-specific CREs and the discovery of molecular mechanisms underlying TF-dependent CRE activity., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2020. Published by The Company of Biologists Ltd.)

Atrial fibrillation (AF), defined by disorganized atrial cardiac rhythm, is the most prevalent cardiac arrhythmia worldwide. Recent genetic studies have highlighted a major heritable component and identified numerous loci associated with AF risk, including the cardiogenic transcription factor genes TBX5, GATA4, and NKX2-5. We report that Tbx5 and Gata4 interact with opposite signs for atrial rhythm controls compared with cardiac development. Using mouse genetics, we found that AF pathophysiology caused by Tbx5 haploinsufficiency, including atrial arrhythmia susceptibility, prolonged action potential duration, and ectopic cardiomyocyte depolarizations, were all rescued by Gata4 haploinsufficiency. In contrast, Nkx2-5 haploinsufficiency showed no combinatorial effect. The molecular basis of the TBX5/GATA4 interaction included normalization of intra-cardiomyocyte calcium flux and expression of calcium channel genes Atp2a2 and Ryr2. Furthermore, GATA4 and TBX5 showed antagonistic interactions on an Ryr2 enhancer. Atrial rhythm instability caused by Tbx5 haploinsufficiency was rescued by a decreased dose of phospholamban, a sarco/endoplasmic reticulum Ca2+-ATPase inhibitor, consistent with a role for decreased sarcoplasmic reticulum calcium flux in Tbx5-dependent AF susceptibility. This work defines a link between Tbx5 dose, sarcoplasmic reticulum calcium flux, and AF propensity. The unexpected interactions between Tbx5 and Gata4 in atrial rhythm control suggest that evaluating specific interactions between genetic risk loci will be necessary for ascertaining personalized risk from genetic association data.

Dominant mutations of Gata4, an essential cardiogenic transcription factor (TF), were known to cause outflow tract (OFT) defects in both human and mouse, but the underlying molecular mechanism was not clear. In this study, Gata4 haploinsufficiency in mice was found to result in OFT defects including double outlet right ventricle (DORV) and ventricular septum defects (VSDs). Gata4 was shown to be required for Hedgehog (Hh)-receiving progenitors within the second heart field (SHF) for normal OFT alignment. Restored cell proliferation in the SHF by knocking-down Pten failed to rescue OFT defects, suggesting that additional cell events under Gata4 regulation is important. SHF Hh-receiving cells failed to migrate properly into the proximal OFT cushion, which is associated with abnormal EMT and cell proliferation in Gata4 haploinsufficiency. The genetic interaction of Hh signaling and Gata4 is further demonstrated to be important for OFT development. Gata4 and Smo double heterozygotes displayed more severe OFT abnormalities including persistent truncus arteriosus (PTA). Restoration of Hedgehog signaling renormalized SHF cell proliferation and migration, and rescued OFT defects in Gata4 haploinsufficiency. In addition, there was enhanced Gata6 expression in the SHF of the Gata4 heterozygotes. The Gata4-responsive repressive sites were identified within 1kbp upstream of the transcription start site of Gata6 by both ChIP-qPCR and luciferase reporter assay. These results suggested a SHF regulatory network comprising of Gata4, Gata6 and Hh-signaling for OFT development., Competing Interests: The authors have declared that no competing interests exist.

Risk for Atrial Fibrillation (AF), the most common human arrhythmia, has a major genetic component. The T-box transcription factor TBX5 influences human AF risk, and adult-specific Tbx5 -mutant mice demonstrate spontaneous AF. We report that TBX5 is critical for cellular Ca 2+ homeostasis, providing a molecular mechanism underlying the genetic implication of TBX5 in AF. We show that cardiomyocyte action potential (AP) abnormalities in Tbx5 -deficient atrial cardiomyocytes are caused by a decreased sarcoplasmic reticulum (SR) Ca 2+ ATPase (SERCA2)-mediated SR calcium uptake which was balanced by enhanced trans-sarcolemmal calcium fluxes (calcium current and sodium/calcium exchanger), providing mechanisms for triggered activity. The AP defects, cardiomyocyte ectopy, and AF caused by TBX5 deficiency were rescued by phospholamban removal, which normalized SERCA function. These results directly link transcriptional control of SERCA2 activity, depressed SR Ca 2+ sequestration, enhanced trans-sarcolemmal calcium fluxes, and AF, establishing a mechanism underlying the genetic basis for a Ca 2+ -dependent pathway for AF risk., Competing Interests: WD, BL, LT, KS, RN, FA, SM, SL, MG, YW, YL, HV, LS, MB, IM, CW No competing interests declared, (© 2019, Dai et al.)

Biomechanical cues dynamically control major cellular processes, but whether genetic variants actively participate in mechanosensing mechanisms remains unexplored. Vascular homeostasis is tightly regulated by hemodynamics. Exposure to disturbed blood flow at arterial sites of branching and bifurcation causes constitutive activation of vascular endothelium contributing to atherosclerosis, the major cause of coronary artery disease (CAD) and ischemic stroke (IS). Conversely, unidirectional flow promotes quiescent endothelium. Genome-wide association studies (GWAS) have identified chromosome 1p32.2 as strongly associated with CAD/IS; however, the causal mechanism related to this locus remains unknown. Using statistical analyses, assay of transposase accessible chromatin with whole-genome sequencing (ATAC-seq), H3K27ac/H3K4me2 ChIP with whole-genome sequencing (ChIP-seq), and CRISPR interference in human aortic endothelial cells (HAECs), our results demonstrate that rs17114036, a common noncoding polymorphism at 1p32.2, is located in an endothelial enhancer dynamically regulated by hemodynamics. CRISPR-Cas9-based genome editing shows that rs17114036-containing region promotes endothelial quiescence under unidirectional shear stress by regulating phospholipid phosphatase 3 (PLPP3). Chromatin accessibility quantitative trait locus (caQTL) mapping using HAECs from 56 donors, allelic imbalance assay from 7 donors, and luciferase assays demonstrate that CAD/IS-protective allele at rs17114036 in PLPP3 intron 5 confers increased endothelial enhancer activity. ChIP-PCR and luciferase assays show that CAD/IS-protective allele at rs17114036 creates a binding site for transcription factor Krüppel-like factor 2 (KLF2), which increases the enhancer activity under unidirectional flow. These results demonstrate that a human SNP contributes to critical endothelial mechanotransduction mechanisms and suggest that human haplotypes and related cis -regulatory elements provide a previously unappreciated layer of regulatory control in cellular mechanosensing mechanisms., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

Codevelopment of the lungs and heart underlies key evolutionary innovations in the transition to terrestrial life. Cardiac specializations that support pulmonary circulation, including the atrial septum, are generated by second heart field (SHF) cardiopulmonary progenitors (CPPs). It has been presumed that transcription factors required in the SHF for cardiac septation, e.g., Tbx5 , directly drive a cardiac morphogenesis gene-regulatory network. Here, we report instead that TBX5 directly drives Wnt ligands to initiate a bidirectional signaling loop between cardiopulmonary mesoderm and the foregut endoderm for endodermal pulmonary specification and, subsequently, atrial septation. We show that Tbx5 is required for pulmonary specification in mice and amphibians but not for swim bladder development in zebrafish. TBX5 is non-cell-autonomously required for pulmonary endoderm specification by directly driving Wnt2 and Wnt2b expression in cardiopulmonary mesoderm. TBX5 ChIP-sequencing identified cis -regulatory elements at Wnt2 sufficient for endogenous Wnt2 expression domains in vivo and required for Wnt2 expression in precardiac mesoderm in vitro. Tbx5 cooperated with Shh signaling to drive Wnt2b expression for lung morphogenesis. Tbx5 haploinsufficiency in mice, a model of Holt-Oram syndrome, caused a quantitative decrement of mesodermal-to-endodermal Wnt signaling and subsequent endodermal-to-mesodermal Shh signaling required for cardiac morphogenesis. Thus, Tbx5 initiates a mesoderm-endoderm-mesoderm signaling loop in lunged vertebrates that provides a molecular basis for the coevolution of pulmonary and cardiac structures required for terrestrial life., Competing Interests: The authors declare no conflict of interest.

The mammalian heart has a limited regenerative capacity and typically progresses to heart failure following injury. Here, we defined a hedgehog (HH)-Gli1-Mycn network for cardiomyocyte proliferation and heart regeneration from amphibians to mammals. Using a genome-wide screen, we verified that HH signaling was essential for heart regeneration in the injured newt. Next, pharmacological and genetic loss- and gain-of-function of HH signaling demonstrated the essential requirement for HH signaling in the neonatal, adolescent, and adult mouse heart regeneration, and in the proliferation of hiPSC-derived cardiomyocytes. Fate-mapping and molecular biological studies revealed that HH signaling, via a HH-Gli1-Mycn network, contributed to heart regeneration by inducing proliferation of pre-existing cardiomyocytes and not by de novo cardiomyogenesis. Further, Mycn mRNA transfection experiments recapitulated the effects of HH signaling and promoted adult cardiomyocyte proliferation. These studies defined an evolutionarily conserved function of HH signaling that may serve as a platform for human regenerative therapies.

Over 500 genetic loci have been associated with risk of cardiovascular diseases (CVDs); however, most loci are located in gene-distal non-coding regions and their target genes are not known. Here, we generated high-resolution promoter capture Hi-C (PCHi-C) maps in human induced pluripotent stem cells (iPSCs) and iPSC-derived cardiomyocytes (CMs) to provide a resource for identifying and prioritizing the functional targets of CVD associations. We validate these maps by demonstrating that promoters preferentially contact distal sequences enriched for tissue-specific transcription factor motifs and are enriched for chromatin marks that correlate with dynamic changes in gene expression. Using the CM PCHi-C map, we linked 1999 CVD-associated SNPs to 347 target genes. Remarkably, more than 90% of SNP-target gene interactions did not involve the nearest gene, while 40% of SNPs interacted with at least two genes, demonstrating the importance of considering long-range chromatin interactions when interpreting functional targets of disease loci., Competing Interests: LM, DS, NS, IA, AJ, GH, GB, IM, EM, MN No competing interests declared, (© 2018, Montefiori et al.)

The National Institutes of Health have made substantial investments in genomic studies and technologies to identify DNA sequence variants associated with human disease phenotypes. The National Heart, Lung, and Blood Institute has been at the forefront of these commitments to ascertain genetic variation associated with heart, lung, blood, and sleep diseases and related clinical traits. Genome-wide association studies, exome- and genome-sequencing studies, and exome-genotyping studies of the National Heart, Lung, and Blood Institute-funded epidemiological and clinical case-control studies are identifying large numbers of genetic variants associated with heart, lung, blood, and sleep phenotypes. However, investigators face challenges in identification of genomic variants that are functionally disruptive among the myriad of computationally implicated variants. Studies to define mechanisms of genetic disruption encoded by computationally identified genomic variants require reproducible, adaptable, and inexpensive methods to screen candidate variant and gene function. High-throughput strategies will permit a tiered variant discovery and genetic mechanism approach that begins with rapid functional screening of a large number of computationally implicated variants and genes for discovery of those that merit mechanistic investigation. As such, improved variant-to-gene and gene-to-function screens-and adequate support for such studies-are critical to accelerating the translation of genomic findings. In this White Paper, we outline the variety of novel technologies, assays, and model systems that are making such screens faster, cheaper, and more accurate, referencing published work and ongoing work supported by the National Heart, Lung, and Blood Institute's R21/R33 Functional Assays to Screen Genomic Hits program. We discuss priorities that can accelerate the impressive but incomplete progress represented by big data genomic research., (© 2018 American Heart Association, Inc.)

The noncoding genome is pervasively transcribed. Noncoding RNAs (ncRNAs) generated from enhancers have been proposed as a general facet of enhancer function and some have been shown to be required for enhancer activity. Here we examine the transcription-factor-(TF)-dependence of ncRNA expression to define enhancers and enhancer-associated ncRNAs that are involved in a TF-dependent regulatory network. TBX5, a cardiac TF, regulates a network of cardiac channel genes to maintain cardiac rhythm. We deep sequenced wildtype and Tbx5 -mutant mouse atria, identifying ~2600 novel Tbx5 -dependent ncRNAs. Tbx5-dependent ncRNAs were enriched for tissue-specific marks of active enhancers genome-wide. Tbx5-dependent ncRNAs emanated from regions that are enriched for TBX5-binding and that demonstrated Tbx5-dependent enhancer activity. Tbx5 -dependent ncRNA transcription provided a quantitative metric of Tbx5 -dependent enhancer activity, correlating with target gene expression. We identified RACER , a novel Tbx5 -dependent long noncoding RNA (lncRNA) required for the expression of the calcium-handling gene Ryr2 . We illustrate that TF-dependent enhancer transcription can illuminate components of TF-dependent gene regulatory networks.

We recently described a new class of long noncoding RNAs (lncRNAs) that are distinguished by especially tight chromatin association and whose presence is strongly correlated to expression of nearby genes. Here, we examine the cis-enhancer mechanism of this class of chromatin-enriched RNA (cheRNA) across multiple human cell lines. cheRNAs are largely cell type specific and provide the most reliable chromatin signature to predict cis-gene transcription in every human cell type examined. Targeted depletion of three cheRNAs decreases expression of their neighboring genes, indicating potential co-activator function, and single-molecule fluorescence in situ hybridization (smFISH) of one cheRNA-distal target gene pair suggests a spatial overlap consistent with a role in chromosome looping. Additionally, the cheRNA HIDALGO stimulates the fetal hemoglobin subunit gamma 1 (HBG1) gene during erythroid differentiation by promoting contacts to a downstream enhancer. Our results suggest that multiple cheRNAs activate proximal lineage-specific gene transcription.

GATA4 , an essential cardiogenic transcription factor, provides a model for dominant transcription factor mutations in human disease. Dominant GATA4 mutations cause congenital heart disease (CHD), specifically atrial and atrioventricular septal defects (ASDs and AVSDs). We found that second heart field (SHF)-specific Gata4 heterozygote embryos recapitulated the AVSDs observed in germline Gata4 heterozygote embryos. A proliferation defect of SHF atrial septum progenitors and hypoplasia of the dorsal mesenchymal protrusion, rather than anlage of the atrioventricular septum, were observed in this model. Knockdown of the cell-cycle repressor phosphatase and tensin homolog ( Pten ) restored cell-cycle progression and rescued the AVSDs. Gata4 mutants also demonstrated Hedgehog (Hh) signaling defects. Gata4 acts directly upstream of Hh components: Gata4 activated a cis -regulatory element at Gli1 in vitro and occupied the element in vivo. Remarkably, SHF-specific constitutive Hh signaling activation rescued AVSDs in Gata4 SHF-specific heterozygous knockout embryos. Pten expression was unchanged in Smoothened mutants, and Hh pathway genes were unchanged in Pten mutants, suggesting pathway independence. Thus, both the cell-cycle and Hh-signaling defects caused by dominant Gata4 mutations were required for CHD pathogenesis, suggesting a combinatorial model of disease causation by transcription factor haploinsufficiency.

TBX5 is a member of the T-box transcription factor family and is primarily known for its role in cardiac and forelimb development. Human patients with dominant mutations in TBX5 are characterized by Holt-Oram syndrome, and show defects of the cardiac septa, cardiac conduction system, and the anterior forelimb. The range of cardiac defects associated with TBX5 mutations in humans suggests multiple roles for the transcription factor in cardiac development and function. Animal models demonstrate similar defects and have provided a useful platform for investigating the roles of TBX5 during embryonic development. During early cardiac development, TBX5 appears to act primarily as a transcriptional activator of genes associated with cardiomyocyte maturation and upstream of morphological signals for septation. During later cardiac development, TBX5 is required for patterning of the cardiac conduction system and maintenance of mature cardiomyocyte function. A comprehensive understanding of the integral roles of TBX5 throughout cardiac development and adult life will be critical for understanding human cardiac morphology and function., (© 2017 Elsevier Inc. All rights reserved.)

Cardiac rhythm is extremely robust, generating 2 billion contraction cycles during the average human life span. Transcriptional control of cardiac rhythm is poorly understood. We found that removal of the transcription factor gene Tbx5 from the adult mouse caused primary spontaneous and sustained atrial fibrillation (AF). Atrial cardiomyocytes from the Tbx5-mutant mice exhibited action potential abnormalities, including spontaneous depolarizations, which were rescued by chelating free calcium. We identified a multitiered transcriptional network that linked seven previously defined AF risk loci: TBX5 directly activated PITX2, and TBX5 and PITX2 antagonistically regulated membrane effector genes Scn5a, Gja1, Ryr2, Dsp, and Atp2a2 In addition, reduced Tbx5 dose by adult-specific haploinsufficiency caused decreased target gene expression, myocardial automaticity, and AF inducibility, which were all rescued by Pitx2 haploinsufficiency in mice. These results defined a transcriptional architecture for atrial rhythm control organized as an incoherent feed-forward loop, driven by TBX5 and modulated by PITX2. TBX5/PITX2 interplay provides tight control of atrial rhythm effector gene expression, and perturbation of the co-regulated network caused AF susceptibility. This work provides a model for the molecular mechanisms underpinning the genetic implication of multiple AF genome-wide association studies loci and will contribute to future efforts to stratify patients for AF risk by genotype., Competing Interests: The authors declare that they have no competing interests., (Copyright © 2016, American Association for the Advancement of Science.)

Atrioventricular septal defects (AVSDs) are a common severe form of congenital heart disease (CHD). In this study we identified deleterious non-synonymous mutations in two cilia genes, Dnah11 and Mks1, in independent N-ethyl-N-nitrosourea-induced mouse mutant lines with heritable recessive AVSDs by whole-exome sequencing. Cilia are required for left/right body axis determination and second heart field (SHF) Hedgehog (Hh) signaling, and we find that cilia mutations affect these requirements differentially. Dnah11 avc 4 did not disrupt SHF Hh signaling and caused AVSDs only concurrently with heterotaxy, a left/right axis abnormality. In contrast, Mks1 avc 6 disrupted SHF Hh signaling and caused AVSDs without heterotaxy. We performed unbiased whole-genome SHF transcriptional profiling and found that cilia motility genes were not expressed in the SHF whereas cilia structural and signaling genes were highly expressed. SHF cilia gene expression predicted the phenotypic concordance between AVSDs and heterotaxy in mice and humans with cilia gene mutations. A two-step model of cilia action accurately predicted the AVSD/heterotaxyu phenotypic expression pattern caused by cilia gene mutations. We speculate that cilia gene mutations contribute to both syndromic and non-syndromic AVSDs in humans and provide a model that predicts the phenotypic consequences of specific cilia gene mutations., (© The Author 2016. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Congenital heart disease (CHD) has a complex genetic etiology, and recent studies suggest that high penetrance de novo mutations may account for only a small fraction of disease. In a multi-institutional cohort surveyed by exome sequencing, combining analysis of 987 individuals (discovery cohort of 59 affected trios and 59 control trios, and a replication cohort of 100 affected singletons and 533 unaffected singletons) we observe variation at novel and known loci related to a specific cardiac malformation the atrioventricular septal defect (AVSD). In a primary analysis, by combining developmental coexpression networks with inheritance modeling, we identify a de novo mutation in the DNA binding domain of NR1D2 (p.R175W). We show that p.R175W changes the transcriptional activity of Nr1d2 using an in vitro transactivation model in HUVEC cells. Finally, we demonstrate previously unrecognized cardiovascular malformations in the Nr1d2tm1-Dgen knockout mouse. In secondary analyses we map genetic variation to protein-interaction networks suggesting a role for two collagen genes in AVSD, which we corroborate by burden testing in a second replication cohort of 100 AVSDs and 533 controls (p = 8.37e-08). Finally, we apply a rare-disease inheritance model to identify variation in genes previously associated with CHD (ZFPM2, NSD1, NOTCH1, VCAN, and MYH6), cardiac malformations in mouse models (ADAM17, CHRD, IFT140, PTPRJ, RYR1 and ATE1), and hypomorphic alleles of genes causing syndromic CHD (EHMT1, SRCAP, BBS2, NOTCH2, and KMT2D) in 14 of 59 trios, greatly exceeding variation in control trios without CHD (p = 9.60e-06). In total, 32% of trios carried at least one putatively disease-associated variant across 19 loci,suggesting that inherited and de novo variation across a heterogeneous group of loci may contribute to disease risk.

Human mutations in the cardiac transcription factor gene TBX5 cause congenital heart disease (CHD), although the underlying mechanism is unknown. We report characterization of the endogenous TBX5 cardiac interactome and demonstrate that TBX5, long considered a transcriptional activator, interacts biochemically and genetically with the nucleosome remodeling and deacetylase (NuRD) repressor complex. Incompatible gene programs are repressed by TBX5 in the developing heart. CHD mis-sense mutations that disrupt the TBX5-NuRD interaction cause depression of a subset of repressed genes. Furthermore, the TBX5-NuRD interaction is required for heart development. Phylogenetic analysis showed that the TBX5-NuRD interaction domain evolved during early diversification of vertebrates, simultaneous with the evolution of cardiac septation. Collectively, this work defines a TBX5-NuRD interaction essential to cardiac development and the evolution of the mammalian heart, and when altered may contribute to human CHD., (Copyright © 2016 Elsevier Inc. All rights reserved.)

A mutant transcription factor that has been linked to congenital heart disease has wider effects than previously thought.

Background: Nearly half of all individuals with Down Syndrome (DS) have some type of congenital heart defect (CHD), suggesting that DS sensitizes to CHD but does not cause it. We used a common mouse model of DS, the Ts65Dn mouse, to study the contribution of Tbx5, a known modifier of CHD, to heart defects on a trisomic backgroun. Mice that were heterozygous for a Tbx5 null allele were crossed with Ts65Dn mice. Thoraxes of progeny were fixed in 10% formalin, embedded in paraffin, and sectioned for analysis of CHD. Gene expression in embryonic hearts was examined by quantitative PCR and in situ hybridization. A TBX5 DNA binding site was verified by luciferase assays., Methods: Mice that were heterozygous for a Tbx5 null allele were crossed with Ts65Dn mice. Thoraxes of progeny were fixed in 10% formalin, embedded in paraffin, and sectioned for analysis of CHD. Gene expression in embryonic hearts was examined by quantitative PCR and in situ hybridization. A TBX5 DNA binding site was verified by luciferase assays., Results: We crossed mice that were heterozygous for a Tbx5 null allele with Ts65Dn mice. Mice that were trisomic and carried the Tbx5 mutation (Ts65Dn;Tbx5 (+/-) ) had a significantly increased incidence of overriding aorta compared to their euploid littermates. Ts65Dn;Tbx5 (+/-) mice also showed reduced expression of Pitx2, a molecular marker for the left atrium. Transcript levels of the trisomic Adamts1 gene were decreased in Tbx5 (+/-) mice compared to their euploid littermates. Evidence of a valid binding site for TBX5 upstream of the trisomic Adamts1 locus was also shown., Conclusion: Haploinsufficiency of Tbx5 and trisomy affects alignment of the aorta and this effect may stem from deviations from normal left-right patterning in the heart. We have unveiled a previously unknown interaction between the Tbx5 gene and trisomy, suggesting a connection between Tbx5 and trisomic genes important during heart development.

Zfp57 is a maternal-zygotic effect gene that maintains genomic imprinting. Here we report that Zfp57 mutants exhibited a variety of cardiac defects including atrial septal defect (ASD), ventricular septal defect (VSD), thin myocardium, and reduced trabeculation. Zfp57 maternal-zygotic mutant embryos displayed more severe phenotypes with higher penetrance than the zygotic ones. Cardiac progenitor cells exhibited proliferation and differentiation defects in Zfp57 mutants. ZFP57 is a master regulator of genomic imprinting, so the DNA methylation imprint was lost in embryonic heart without ZFP57. Interestingly, the presence of imprinted DLK1, a target of ZFP57, correlated with NOTCH1 activation in cardiac cells. These results suggest that ZFP57 may modulate NOTCH signaling during cardiac development. Indeed, loss of ZFP57 caused loss of NOTCH1 activation in embryonic heart with more severe loss observed in the maternal-zygotic mutant. Maternal and zygotic functions of Zfp57 appear to play redundant roles in NOTCH1 activation and cardiomyocyte differentiation. This serves as an example of a maternal effect that can influence mammalian organ development. It also links genomic imprinting to NOTCH signaling and particular developmental functions.

Cardiac left ventricular outflow tract (LVOT) defects represent a common but heterogeneous subset of congenital heart disease for which gene identification has been difficult. We describe a 46,XY,t(1;5)(p36.11;q31.2)dn translocation carrier with pervasive developmental delay who also exhibited LVOT defects, including bicuspid aortic valve (BAV), coarctation of the aorta (CoA) and patent ductus arteriosus (PDA). The 1p breakpoint disrupts the 5' UTR of AHDC1, which encodes AT-hook DNA-binding motif containing-1 protein, and AHDC1-truncating mutations have recently been described in a syndrome that includes developmental delay, but not congenital heart disease [Xia, F., Bainbridge, M.N., Tan, T.Y., Wangler, M.F., Scheuerle, A.E., Zackai, E.H., Harr, M.H., Sutton, V.R., Nalam, R.L., Zhu, W. et al. (2014) De Novo truncating mutations in AHDC1 in individuals with syndromic expressive language delay, hypotonia, and sleep apnea. Am. J. Hum. Genet., 94, 784-789]. On the other hand, the 5q translocation breakpoint disrupts the 3' UTR of MATR3, which encodes the nuclear matrix protein Matrin 3, and mouse Matr3 is strongly expressed in neural crest, developing heart and great vessels, whereas Ahdc1 is not. To further establish MATR3 3' UTR disruption as the cause of the proband's LVOT defects, we prepared a mouse Matr3(Gt-ex13) gene trap allele that disrupted the 3' portion of the gene. Matr3(Gt-ex13) homozygotes are early embryo lethal, but Matr3(Gt-ex13) heterozygotes exhibit incompletely penetrant BAV, CoA and PDA phenotypes similar to those in the human proband, as well as ventricular septal defect (VSD) and double-outlet right ventricle (DORV). Both the human MATR3 translocation breakpoint and the mouse Matr3(Gt-ex13) gene trap insertion disturb the polyadenylation of MATR3 transcripts and alter Matrin 3 protein expression, quantitatively or qualitatively. Thus, subtle perturbations in Matrin 3 expression appear to cause similar LVOT defects in human and mouse., (© The Author 2015. Published by Oxford University Press.)

The Second Heart Field (SHF) has been implicated in several forms of congenital heart disease (CHD), including atrioventricular septal defects (AVSDs). Identifying the SHF gene regulatory networks required for atrioventricular septation is therefore an essential goal for understanding the molecular basis of AVSDs. We defined a SHF Hedgehog-dependent gene regulatory network using whole genome transcriptional profiling and GLI-chromatin interaction studies. The Forkhead box transcription factors Foxf1a and Foxf2 were identified as SHF Hedgehog targets. Compound haploinsufficiency for Foxf1a and Foxf2 caused atrioventricular septal defects, demonstrating the biological relevance of this regulatory network. We identified a Foxf1a cis-regulatory element that bound the Hedgehog transcriptional regulators GLI1 and GLI3 and the T-box transcription factor TBX5 in vivo. GLI1 and TBX5 synergistically activated transcription from this cis-regulatory element in vitro. This enhancer drove reproducible expression in vivo in the posterior SHF, the only region where Gli1 and Tbx5 expression overlaps. Our findings implicate Foxf genes in atrioventricular septation, describe the molecular underpinnings of the genetic interaction between Hedgehog signaling and Tbx5, and establish a molecular model for the selection of the SHF gene regulatory network for cardiac septation.