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.

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.

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.

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.

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 Ca2+ 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) Ca2+ 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 Ca2+ sequestration, enhanced trans-sarcolemmal calcium fluxes, and AF, establishing a mechanism underlying the genetic basis for a Ca2+-dependent pathway for AF risk.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

Mutations in the nuclear Lamin A/C gene (LMNA) cause diverse degenerative disorders, including malignant dilated cardiomyopathy in adults. A prevailing hypothesis postulates that LMNA mutations cause nuclear envelope ruptures that trigger pathogenic inflammatory signaling via the cGAS-STING cytosolic DNA-sensing pathway. Here, we provide evidence against this hypothesis, using a mouse model of LMNA-related cardiomyopathy that mimics Lamin A/C protein reduction observed in patient cardiomyocytes. We observed that pervasive nuclear envelope ruptures preceded the onset of cardiac transcriptional modulation and dilated cardiomyopathy. Nuclear ruptures activated DNA damage response without causing immediate cardiomyocyte death. However, cGAS-STING downstream cytokine genes remained inactive in the mutant cardiomyocytes. Deleting cGas or Sting did not alleviate cardiomyopathy. Instead, extracellular matrix signaling was predicted to emanate from Lamin A/C-reduced cardiomyocytes to communicate with fibroblasts in the heart. These findings suggest that cGAS-STING is not a major pathogenetic contributor to LMNA-related dilated cardiomyopathy in adult humans.

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.