A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network.

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.

LATS1/2 control TGFB-directed epithelial-to-mesenchymal transition in the murine dorsal cranial neuroepithelium through YAP regulation.

Hippo signaling, an evolutionarily conserved kinase cascade involved in organ size control, plays key roles in various tissue developmental processes, but its role in craniofacial development remains poorly understood. Using the transgenic Wnt1-Cre2 driver, we inactivated the Hippo signaling components Lats1 and Lats2 in the cranial neuroepithelium of mouse embryos and found that the double conditional knockout (DCKO) of Lats1/2 resulted in neural tube and craniofacial defects. Lats1/2 DCKO mutant embryos had microcephaly with delayed and defective neural tube closure. Furthermore, neuroepithelial cell shape and architecture were disrupted within the cranial neural tube in Lats1/2 DCKO mutants. RNA sequencing of embryonic neural tubes revealed increased TGFB signaling in Lats1/2 DCKO mutants. Moreover, markers of epithelial-to-mesenchymal transition (EMT) were upregulated in the cranial neural tube. Inactivation of Hippo signaling downstream effectors, Yap and Taz, suppressed neuroepithelial defects, aberrant EMT and TGFB upregulation in Lats1/2 DCKO embryos, indicating that LATS1/2 function via YAP and TAZ. Our findings reveal important roles for Hippo signaling in modulating TGFB signaling during neural crest EMT.

Transcriptional Dysregulation Underlies Both Monogenic Arrhythmia Syndrome and Common Modifiers of Cardiac Repolarization.

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 (INa), 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 INa, and that reduced PDGF receptor (PDGFRA [platelet-derived growth factor receptor α]) expression and blunted signal transduction to PI3K was implicated in enhanced late INa. 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.

Transcriptional Patterning of the Ventricular Cardiac Conduction System.

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.

Evolutionarily conserved Tbx5-Wnt2/2b pathway orchestrates cardiopulmonary development.

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.

Correction: Foxf Genes Integrate Tbx5 and Hedgehog Pathways in the Second Heart Field for Cardiac Septation.

[This corrects the article DOI: 10.1371/journal.pgen.1004604.].

De Novo and Rare Variants at Multiple Loci Support the Oligogenic Origins of Atrioventricular Septal Heart Defects.

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.

Cilia gene mutations cause atrioventricular septal defects by multiple mechanisms.

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. Dnah11avc4 did not disrupt SHF Hh signaling and caused AVSDs only concurrently with heterotaxy, a left/right axis abnormality. In contrast, Mks1avc6 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.

Foxf genes integrate tbx5 and hedgehog pathways in the second heart field for cardiac septation.

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.

The pattern of congenital heart defects arising from reduced Tbx5 expression is altered in a Down syndrome mouse model.

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.

Endothelial cells adopt a pro-reparative immune responsive signature during cardiac injury.

Modulation of the heart's immune microenvironment is crucial for recovery after ischemic events such as myocardial infarction (MI). Endothelial cells (ECs) can have immune regulatory functions; however, interactions between ECs and the immune environment in the heart after MI remain poorly understood. We identified an EC-specific IFN responsive and immune regulatory gene signature in adult and pediatric heart failure (HF) tissues. Single-cell transcriptomic analysis of murine hearts subjected to MI uncovered an EC population (IFN-ECs) with immunologic gene signatures similar to those in human HF. IFN-ECs were enriched in regenerative-stage mouse hearts and expressed genes encoding immune responsive transcription factors (Irf7, Batf2, and Stat1). Single-cell chromatin accessibility studies revealed an enrichment of these TF motifs at IFN-EC signature genes. Expression of immune regulatory ligand genes by IFN-ECs suggests bidirectional signaling between IFN-ECs and macrophages in regenerative-stage hearts. Our data suggest that ECs may adopt immune regulatory signatures after cardiac injury to accompany the reparative response. The presence of these signatures in human HF and murine MI models suggests a potential role for EC-mediated immune regulation in responding to stress induced by acute injury in MI and chronic adverse remodeling in HF.

Regrowing the heart, one TREE at a time.

While many animals can completely repair injured tissues, the mammalian heart possesses limited regenerative capabilities. Yan and Cigliola et al. show that AAV-mediated, zebrafish-derived tissue regeneration enhancer elements (TREEs) can direct pro-regenerative gene expression in injured cardiac tissue of mice and pigs that turn off following repair.

ETV2 primes hematoendothelial gene enhancers prior to hematoendothelial fate commitment.

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.

Decoding the PITX2-controlled genetic network in atrial fibrillation.

Atrial fibrillation (AF), the most common sustained cardiac arrhythmia and a major risk factor for stroke, often arises through ectopic electrical impulses derived from the pulmonary veins (PVs). Sequence variants in enhancers controlling expression of the transcription factor PITX2, which is expressed in the cardiomyocytes (CMs) of the PV and left atrium (LA), have been implicated in AF predisposition. Single nuclei multiomic profiling of RNA and analysis of chromatin accessibility combined with spectral clustering uncovered distinct PV- and LA-enriched CM cell states. Pitx2-mutant PV and LA CMs exhibited gene expression changes consistent with cardiac dysfunction through cell type-distinct, PITX2-directed, cis-regulatory grammars controlling target gene expression. The perturbed network targets in each CM were enriched in distinct human AF predisposition genes, suggesting combinatorial risk for AF genesis. Our data further reveal that PV and LA Pitx2-mutant CMs signal to endothelial and endocardial cells through BMP10 signaling with pathogenic potential. This work provides a multiomic framework for interrogating the basis of AF predisposition in the PVs of humans.

Hedgehog signaling activates a mammalian heterochronic gene regulatory network controlling differentiation timing across lineages.

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.

Tbx5 drives Aldh1a2 expression to regulate a RA-Hedgehog-Wnt gene regulatory network coordinating cardiopulmonary development.

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.

The Cardiac TBX5 Interactome Reveals a Chromatin Remodeling Network Essential for Cardiac Septation.

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.