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

Sex-Specific Response to A1BG Loss Results in Female Dilated Cardiomyopathy

Background

Cardiac disease often manifests differently in terms of frequency and pathology between men and women. However, the mechanisms underlying these differences are not fully understood. The glycoprotein A1BG is necessary for proper cardiac function in females but not males. Despite this, the role of A1BG in the female heart remains poorly studied.

Methods

To determine the sex differential function of A1BG, we generated a novel conditional A1bg allele and a novel conditional A1bg Rosa26 knockin allele. Histology, electrocardiography, transcriptional profiling (RNA-seq), transmission electron microscopy, western blot analyses, mass spectrometry, and immunohistochemistry were used to assess cardiac structure and function.

Results

The study reveals that the absence of A1BG results in significant cardiac dysfunction in female but not male mice. Gene expression underscores that A1BG plays a critical role in metabolic processes and the integrity of intercalated discs in female cardiomyocytes. This dysfunction may be related to sex-specific A1BG cardiac interactomes and manifests as structural and functional alterations in the left ventricle indicative of dilated cardiomyopathy, thus suggesting a sex-specific requirement for A1BG in cardiac health.

Conclusion

The loss of A1BG in cardiomyocytes leads to dilated cardiomyopathy in females, not males.

CHD4 and SMYD1 repress common transcriptional programs in the developing heart

Regulation of chromatin states is essential for proper temporal and spatial gene expression. Chromatin states are modulated by remodeling complexes composed of components that have enzymatic activities. CHD4 is the catalytic core of the nucleosome remodeling and deacetylase (NuRD) complex, which represses gene transcription. However, it remains to be determined how CHD4, a ubiquitous enzyme that remodels chromatin structure, functions in cardiomyocytes to maintain heart development. In particular, whether other proteins besides the NuRD components interact with CHD4 in the heart is controversial. Using quantitative proteomics, we identified that CHD4 interacts with SMYD1, a striated muscle-restricted histone methyltransferase that is essential for cardiomyocyte differentiation and cardiac morphogenesis. Comprehensive transcriptomic and chromatin accessibility studies of Smyd1 and Chd4 null embryonic mouse hearts revealed that SMYD1 and CHD4 repress a group of common genes and pathways involved in glycolysis, response to hypoxia, and angiogenesis. Our study reveals a mechanism by which CHD4 functions during heart development, and a previously uncharacterized mechanism regarding how SMYD1 represses cardiac transcription in the developing heart.

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.

Technologies to Study Genetics and Molecular Pathways

Over the last few decades, the study of congenital heart disease (CHD) has benefited from various model systems and the development of molecular biological techniques enabling the analysis of single gene as well as global effects. In this chapter, we first describe different models including CHD patients and their families, animal models ranging from invertebrates to mammals, and various cell culture systems. Moreover, techniques to experimentally manipulate these models are discussed. Second, we introduce cardiac phenotyping technologies comprising the analysis of mouse and cell culture models, live imaging of cardiogenesis, and histological methods for fixed hearts. Finally, the most important and latest molecular biotechniques are described. These include genotyping technologies, different applications of next-generation sequencing, and the analysis of transcriptome, epigenome, proteome, and metabolome. In summary, the models and technologies presented in this chapter are essential to study the function and development of the heart and to understand the molecular pathways underlying CHD.

Conserved mechanisms of NuRD function in hematopoetic gene expression

The Nucleosome Remodeling and Deacetylating Complex (NuRD) is ubiquitously expressed in all metazoans. It combines nucleosome remodeling and histone deacetylating activities to generate inaccessible chromatin structures and to repress gene transcription. NuRD is involved in the generation and maintenance of a wide variety of lineage-specific gene expression programs during differentiation and in differentiated cells. A close cooperation with a large number of lineage-specific transcription factors is key to allow NuRD to function in many distinct differentiation contexts. The molecular nature of this interplay between transcription factors and NuRD is complex and not well understood. This review uses hematopoiesis as a paradigm to highlight recent advances in our understanding of how transcription factors and NuRD cooperate at the molecular level during differentiation. A comparison of vertebrate and invertebrate systems serves to identify the conserved and fundamental concepts guiding functional interactions between transcription factors and NuRD. We also discuss how the transcription factor-NuRD axis constitutes a potential therapeutic target for the treatment of hemoglobinopathies.

Missense Mutation in Human CHD4 Causes Ventricular Noncompaction by Repressing ADAMTS1

BACKGROUND

Left ventricular noncompaction (LVNC) is a prevalent cardiomyopathy associated with excessive trabeculation and thin compact myocardium. Patients with LVNC are vulnerable to cardiac dysfunction and at high risk of sudden death. Although sporadic and inherited mutations in cardiac genes are implicated in LVNC, understanding of the mechanisms responsible for human LVNC is limited.

METHODS

We screened the complete exome sequence database of the Pediatrics Cardiac Genomics Consortium and identified a cohort with a de novo CHD4 (chromodomain helicase DNA-binding protein 4) proband, CHD4M202I, with congenital heart defects. We engineered a humanized mouse model of CHD4M202I (mouse CHD4M195I). Histological analysis, immunohistochemistry, flow cytometry, transmission electron microscopy, and echocardiography were used to analyze cardiac anatomy and function. Ex vivo culture, immunopurification coupled with mass spectrometry, transcriptional profiling, and chromatin immunoprecipitation were performed to deduce the mechanism of CHD4M195I-mediated ventricular wall defects.

RESULTS

CHD4M195I/M195I mice developed biventricular hypertrabeculation and noncompaction and died at birth. Proliferation of cardiomyocytes was significantly increased in CHD4M195I hearts, and the excessive trabeculation was associated with accumulation of ECM (extracellular matrix) proteins and a reduction of ADAMTS1 (ADAM metallopeptidase with thrombospondin type 1 motif 1), an ECM protease. We rescued the hyperproliferation and hypertrabeculation defects in CHD4M195I hearts by administration of ADAMTS1. Mechanistically, the CHD4M195I protein showed augmented affinity to endocardial BRG1 (SWI/SNF-related, matrix-associated, actin-dependent regulator of chromatin, subfamily A, member 4). This enhanced affinity resulted in the failure of derepression of Adamts1 transcription such that ADAMTS1-mediated trabeculation termination was impaired.

CONCLUSIONS

Our study reveals how a single mutation in the chromatin remodeler CHD4, in mice or humans, modulates ventricular chamber maturation and that cardiac defects associated with the missense mutation CHD4M195I can be attenuated by the administration of ADAMTS1.

In Vivo Dissection of Chamber-Selective Enhancers Reveals Estrogen-Related Receptor as a Regulator of Ventricular Cardiomyocyte Identity

BACKGROUND

Cardiac chamber-selective transcriptional programs underpin the structural and functional differences between atrial and ventricular cardiomyocytes (aCMs and vCMs). The mechanisms responsible for these chamber-selective transcriptional programs remain largely undefined.

METHODS

We nominated candidate chamber-selective enhancers (CSEs) by determining the genome-wide occupancy of 7 key cardiac transcription factors (GATA4, MEF2A, MEF2C, NKX2-5, SRF, TBX5, TEAD1) and transcriptional coactivator P300 in atria and ventricles. Candidate enhancers were tested using an adeno-associated virus-mediated massively parallel reporter assay. Chromatin features of CSEs were evaluated by performing assay of transposase accessible chromatin sequencing and acetylation of histone H3 at lysine 27-HiChIP on aCMs and vCMs. CSE sequence requirements were determined by systematic tiling mutagenesis of 29 CSEs at 5 bp resolution. Estrogen-related receptor (ERR) function in cardiomyocytes was evaluated by Cre-loxP-mediated inactivation of ERRα and ERRγ in cardiomyocytes.

RESULTS

We identified 134 066 and 97 506 regions reproducibly occupied by at least 1 transcription factor or P300, in atria or ventricles, respectively. Enhancer activities of 2639 regions bound by transcription factors or P300 were tested in aCMs and vCMs by adeno-associated virus-mediated massively parallel reporter assay. This identified 1092 active enhancers in aCMs or vCMs. Several overlapped loci associated with cardiovascular disease through genome-wide association studies, and 229 exhibited chamber-selective activity in aCMs or vCMs. Many CSEs exhibited differential chromatin accessibility between aCMs and vCMs, and CSEs were enriched for aCM- or vCM-selective acetylation of histone H3 at lysine 27-anchored loops. Tiling mutagenesis of 29 CSEs identified the binding motif of ERRα/γ as important for ventricular enhancer activity. The requirement of ERRα/γ to activate ventricular CSEs and promote vCM identity was confirmed by loss of the vCM gene profile in ERRα/γ knockout vCMs.

CONCLUSIONS

We identified 229 CSEs that could be useful research tools or direct therapeutic gene expression. We showed that chamber-selective multi-transcription factor, P300 occupancy, open chromatin, and chromatin looping are predictive features of CSEs. We found that ERRα/γ are essential for maintenance of ventricular identity. Finally, our gene expression, epigenetic, 3-dimensional genome, and enhancer activity atlas provide key resources for future studies of chamber-selective gene regulation.

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.

Return of the Tbx5; lineage-tracing reveals ventricular cardiomyocyte-like precursors in the injured adult mammalian heart

The single curative measure for heart failure patients is a heart transplantation, which is limited due to a shortage of donors, the need for immunosuppression and economic costs. Therefore, there is an urgent unmet need for identifying cell populations capable of cardiac regeneration that we will be able to trace and monitor. Injury to the adult mammalian cardiac muscle, often leads to a heart attack through the irreversible loss of a large number of cardiomyocytes, due to an idle regenerative capability. Recent reports in zebrafish indicate that Tbx5a is a vital transcription factor for cardiomyocyte regeneration. Preclinical data underscore the cardioprotective role of Tbx5 upon heart failure. Data from our earlier murine developmental studies have identified a prominent unipotent Tbx5-expressing embryonic cardiac precursor cell population able to form cardiomyocytes, in vivo, in vitro and ex vivo. Using a developmental approach to an adult heart injury model and by employing a lineage-tracing mouse model as well as the use of single-cell RNA-seq technology, we identify a Tbx5-expressing ventricular cardiomyocyte-like precursor population, in the injured adult mammalian heart. The transcriptional profile of that precursor cell population is closer to that of neonatal than embryonic cardiomyocyte precursors. Tbx5, a cardinal cardiac development transcription factor, lies in the center of a ventricular adult precursor cell population, which seems to be affected by neurohormonal spatiotemporal cues. The identification of a Tbx5-specific cardiomyocyte precursor-like cell population, which is capable of dedifferentiating and potentially deploying a cardiomyocyte regenerative program, provides a clear target cell population for translationally-relevant heart interventional studies.

The transcription factor Tbx5 regulates direction-selective retinal ganglion cell development and image stabilization

The diversity of visual input processed by the mammalian visual system requires the generation of many distinct retinal ganglion cell (RGC) types, each tuned to a particular feature. The molecular code needed to generate this cell-type diversity is poorly understood. Here, we focus on the molecules needed to specify one type of retinal cell: the upward-preferring ON direction-selective ganglion cell (up-oDSGC) of the mouse visual system. Single-cell transcriptomic profiling of up- and down-oDSGCs shows that the transcription factor Tbx5 is selectively expressed in up-oDSGCs. The loss of Tbx5 in up-oDSGCs results in a selective defect in the formation of up-oDSGCs and a corresponding inability to detect vertical motion. A downstream effector of Tbx5, Sfrp1, is also critical for vertical motion detection but not up-oDSGC formation. These results advance our understanding of the molecular mechanisms that specify a rare retinal cell type and show how disrupting this specification leads to a corresponding defect in neural circuitry and behavior.

Densely methylated DNA traps Methyl-CpG-binding domain protein 2 but permits free diffusion by Methyl-CpG-binding domain protein 3

The methyl-CpG-binding domain 2 and 3 proteins (MBD2 and MBD3) provide structural and DNA-binding function for the Nucleosome Remodeling and Deacetylase (NuRD) complex. The two proteins form distinct NuRD complexes and show different binding affinity and selectivity for methylated DNA. Previous studies have shown that MBD2 binds with high affinity and selectivity for a single methylated CpG dinucleotide while MBD3 does not. However, the NuRD complex functions in regions of the genome that contain many CpG dinucleotides (CpG islands). Therefore, in this work, we investigate the binding and diffusion of MBD2 and MBD3 on more biologically relevant DNA templates that contain a large CpG island or limited CpG sites. Using a combination of single-molecule and biophysical analyses, we show that both MBD2 and MBD3 diffuse freely and rapidly across unmethylated CpG-rich DNA. In contrast, we found methylation of large CpG islands traps MBD2 leading to stable and apparently static binding on the CpG island while MBD3 continues to diffuse freely. In addition, we demonstrate both proteins bend DNA, which is augmented by methylation. Together, these studies support a model in which MBD2-NuRD strongly localizes to and compacts methylated CpG islands while MBD3-NuRD can freely mobilize nucleosomes independent of methylation status.

Gestational Leucylation Suppresses Embryonic T-Box Transcription Factor 5 Signal and Causes Congenital Heart Disease

Dysregulated maternal nutrition, such as vitamin deficiencies and excessive levels of glucose and fatty acids, increases the risk for congenital heart disease (CHD) in the offspring. However, the association between maternal amino-acid levels and CHD is unclear. Here, it is shown that increased leucine levels in maternal plasma during the first trimester are associated with elevated CHD risk in the offspring. High levels of maternal leucine increase embryonic lysine-leucylation (K-Leu), which is catalyzed by leucyl-tRNA synthetase (LARS). LARS preferentially binds to and catalyzes K-Leu modification of lysine 339 within T-box transcription factor TBX5, whereas SIRT3 removes K-Leu from TBX5. Reversible leucylation retains TBX5 in the cytoplasm and inhibits its transcriptional activity. Increasing embryonic K-Leu levels in high-leucine-diet fed or Sirt3 knockout mice causes CHD in the offspring. Targeting K-Leu using the leucine analogue leucinol can inhibit LARS activity, reverse TBX5 K-Leu modification, and decrease the occurrence of CHD in high-leucine-diet fed mice. This study reveals that increased maternal leucine levels increases CHD risk in the offspring through inhibition of embryonic TBX5 signaling, indicating that leucylation exerts teratogenic effects during heart development and may be an intervening target of CHD.

CHD4 is recruited by GATA4 and NKX2-5 to repress noncardiac gene programs in the developing heart

The nucleosome remodeling and deacetylase (NuRD) complex is one of the central chromatin remodeling complexes that mediates gene repression. NuRD is essential for numerous developmental events, including heart development. Clinical and genetic studies have provided direct evidence for the role of chromodomain helicase DNA-binding protein 4 (CHD4), the catalytic component of NuRD, in congenital heart disease (CHD), including atrial and ventricular septal defects. Furthermore, it has been demonstrated that CHD4 is essential for mammalian cardiomyocyte formation and function. A key unresolved question is how CHD4/NuRD is localized to specific cardiac target genes, as neither CHD4 nor NuRD can directly bind DNA. Here, we coupled a bioinformatics-based approach with mass spectrometry analyses to demonstrate that CHD4 interacts with the core cardiac transcription factors GATA4, NKX2-5, and TBX5 during embryonic heart development. Using transcriptomics and genome-wide occupancy data, we characterized the genomic landscape of GATA4, NKX2-5, and TBX5 repression and defined the direct cardiac gene targets of the GATA4-CHD4, NKX2-5-CHD4, and TBX5-CHD4 complexes. These data were used to identify putative cis-regulatory elements controlled by these complexes. We genetically interrogated two of these silencers in vivo: Acta1 and Myh11 We show that deletion of these silencers leads to inappropriate skeletal and smooth muscle gene misexpression, respectively, in the embryonic heart. These results delineate how CHD4/NuRD is localized to specific cardiac loci and explicates how mutations in the broadly expressed CHD4 protein lead to cardiac-specific disease states.

Transcription factor protein interactomes reveal genetic determinants in heart disease

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

DNA Methylation Levels of the TBX5 Gene Promoter Are Associated with Congenital Septal Defects in Mexican Paediatric Patients

The TBX5 gene regulates morphological changes during heart development, and it has been associated with epigenetic abnormalities observed in congenital heart defects (CHD). The aim of this research was to evaluate the association between DNA methylation levels of the TBX5 gene promoter and congenital septal defects. DNA methylation levels of six CpG sites in the TBX5 gene promoter were evaluated using pyrosequencing analysis in 35 patients with congenital septal defects and 48 controls. Average methylation levels were higher in individuals with congenital septal defects than in the controls (p < 0.004). In five CpG sites, we also found higher methylation levels in patients than in the controls (p < 0.05). High methylation levels were associated with congenital septal defects (OR = 3.91; 95% CI = 1.02–14.8; p = 0.045). The analysis of Receiver Operating Characteristic (ROC) showed that the methylation levels of the TBX5 gene could be used as a risk marker for congenital septal defects (AUC = 0.68, 95% CI = 0.56–0.80; p = 0.004). Finally, an analysis of environmental factors indicated that maternal infections increased the risk (OR = 2.90; 95% CI = 1.01–8.33; p = 0.048) of congenital septal defects. Our data suggest that a high DNA methylation of the TBX5 gene could be associated with congenital septal defects.

Proteomic analysis identifies ZMYM2 as endogenous binding partner of TBX18 protein in 293 and A549 cells

The TBX18 transcription factor regulates patterning and differentiation programs in the primordia of many organs yet the molecular complexes in which TBX18 resides to exert its crucial transcriptional function in these embryonic contexts have remained elusive. Here, we used 293 and A549 cells as an accessible cell source to search for endogenous protein interaction partners of TBX18 by an unbiased proteomic approach. We tagged endogenous TBX18 by CRISPR/Cas9 targeted genome editing with a triple FLAG peptide, and identified by anti-FLAG affinity purification and subsequent LC-MS analysis the ZMYM2 protein to be statistically enriched together with TBX18 in both 293 and A549 nuclear extracts. Using a variety of assays, we confirmed binding of TBX18 to ZMYM2, a component of the CoREST transcriptional corepressor complex. Tbx18 is coexpressed with Zmym2 in the mesenchymal compartment of the developing ureter of the mouse, and mutations in TBX18and in ZMYM2 were recently linked to congenital anomalies in the kidney and urinary tract (CAKUT) in line with a possible in vivo relevance of TBX18-ZMYM2 protein interaction in ureter development.

Genetic counseling for congenital heart disease - Practice resource of the national society of genetic counselors

Congenital heart disease (CHD) is an indication which spans multiple specialties across various genetic counseling practices. This practice resource aims to provide guidance on key considerations when approaching counseling for this particular indication while recognizing the rapidly changing landscape of knowledge within this domain. This resource was developed with consensus from a diverse group of certified genetic counselors utilizing literature relevant for CHD genetic counseling practice and is aimed at supporting genetic counselors who encounter this indication in their practice both pre- and postnatally.

Expression Profiling and Functional Analysis of Candidate Col10a1 Regulators Identified by the TRAP Program

Hypertrophic chondrocytes and their specific marker, the type X collagen gene (Col10a1), are critical components of endochondral bone formation during skeletal development. We previously found that Runx2 is an indispensable mouse Col10a1 gene regulator and identified many other transcription factors (TFs) that potentially interact with the 150-bp Col10a1 cis-enhancer. However, the roles of these candidate TFs in Col10a1 expression and chondrocyte hypertrophy have not been elucidated. Here, we focus on 32 candidate TFs recently identified by analyzing the 150-bp Col10a1 enhancer using the transcription factor affinity prediction (TRAP) program. We found that 12 TFs (Hoxa3, Lsx, Evx2, Dlx5, S8, Pax2, Egr2, Mef2a, Barhl2, GKlf, Sox17, and Crx) were significantly upregulated and four TFs (Lhx4, Tbx5, Mef2c, and Hb9) were significantly downregulated in hypertrophic MCT cells, which show upregulation of Col10a1 expression. Most of the differential expression pattern of these TFs conformed with the results obtained from ATDC5 cell model and primary mouse chondrocytes. Notably, Tbx5 was downregulated upon Col10a1 upregulation, overexpression of Tbx5 decreased Col10a1 expression, and knock-down of Tbx5 increased Col10a1 expression in hypertrophic chondrocytes, suggesting that Tbx5 is a negative regulator of Col10a1. We further generated a stable Tbx5-overexpressing ATDC5 cell line and ColX-Tbx5 transgenic mice driven by Col10a1-specific enhancers and promoters. Tbx5 overexpression decreased Col10a1 expression in ATDC5 cells cultured as early as day 7 and in limb tissue on post-natal day 1. Slightly weaker alkaline phosphatase staining was also observed in cell culture on day 7 and in limb digits on embryonic day 17.5, indicating mildly delayed ossification. Further characterization of these candidate Col10a1 transcriptional regulators could help identify novel therapeutic targets for skeletal diseases associated with abnormal chondrocyte hypertrophy.

Hedgehog signalling controls sinoatrial node development and atrioventricular cushion formation

Smoothened is a key receptor of the hedgehog pathway, but the roles of Smoothened in cardiac development remain incompletely understood. In this study, we found that the conditional knockout of Smoothened from the mesoderm impaired the development of the venous pole of the heart and resulted in hypoplasia of the atrium/inflow tract (IFT) and a low heart rate. The blockage of Smoothened led to reduced expression of genes critical for sinoatrial node (SAN) development in the IFT. In a cardiac cell culture model, we identified a Gli2–Tbx5–Hcn4 pathway that controls SAN development. In the mutant embryos, the endocardial-to-mesenchymal transition (EndMT) in the atrioventricular cushion failed, and Bmp signalling was downregulated. The addition of Bmp2 rescued the EndMT in mutant explant cultures. Furthermore, we analysed Gli2⁺ scRNAseq and Tbx5^−/− RNAseq data and explored the potential genes downstream of hedgehog signalling in posterior second heart field derivatives. In conclusion, our study reveals that Smoothened-mediated hedgehog signalling controls posterior cardiac progenitor commitment, which suggests that the mutation of Smoothened might be involved in the aetiology of congenital heart diseases related to the cardiac conduction system and heart valves.

Genome-wide analysis of primary microRNA expression using H3K36me3 ChIP-seq data

MicroRNAs are key players in gene regulatory networks controlling cell homeostasis. Their altered expression has been previously linked to disease outcomes and microRNAs thus serve as biomarkers for disease diagnostics. However, their synthesis and its transcriptional regulation have been challenging to investigate. In this study, we validated the use of H3K36me3 histone modification for the quantification of microRNA transcription levels using data from the ENCODE Consortium and then applied this approach to provide new insight into the cell-type-specific regulation in tissues, cell line models and cardiac disease. In cardiomyocytes derived from patients suffering from septal defects, carrying a G296S mutation in the transcription factor GATA4, we show that microRNA gene transcription is altered in cardiomyocytes carrying this mutation and coincides with novel super-enhancers formed within regulatory domains defined using chromatin interaction profiles. The most prominently elevated primary transcript encodes for let-7a and miR-100 that may target genes in the Hippo signaling pathway. Collectively, our work presents a methodology to quantify microRNA gene expression using histone marker data and paves the way for functional studies of cell-type-specific transcriptional regulation occurring in disease pathology.

Functional and Adaptive Significance of Promoter Mutations That Affect Divergent Myocardial Expressions of TRIM72 in Primates

Cis-regulatory elements play important roles in tissue-specific gene expression and in the evolution of various phenotypes, and mutations in promoters and enhancers may be responsible for adaptations of species to environments. TRIM72 is a highly conserved protein that is involved in energy metabolism. Its expression in the heart varies considerably in primates, with high levels of expression in Old World monkeys and near absence in hominids. Here, we combine phylogenetic hypothesis testing and experimentation to demonstrate that mutations in promoter are responsible for the differences among primate species in the heart-specific expression of TRIM72. Maximum likelihood estimates of lineage-specific substitution rates under local-clock models show that relative to the evolutionary rate of introns, the rate of promoter was accelerated by 78% in the common ancestor of Old World monkeys, suggesting a role for positive selection in the evolution of the TRIM72 promoter, possibly driven by selective pressure due to changes in cardiac physiology after species divergence. We demonstrate that mutations in the TRIM72 promoter account for the differential myocardial TRIM72 expression of the human and the rhesus macaque. Furthermore, changes in TRIM72 expression alter the expression of genes involved in oxidative phosphorylation, which in turn affects mitochondrial respiration and cardiac energy capacity. On a broader timescale, phylogenetic regression analyses of data from 29 mammalian species show that mammals with high cardiac expression of TRIM72 have high heart rate, suggesting that the expression changes of TRIM72 may be related to differences in the heart physiology of those species.

Modeling Human TBX5 Haploinsufficiency Predicts Regulatory Networks for Congenital Heart Disease

Haploinsufficiency of transcriptional regulators causes human congenital heart disease (CHD); however, the underlying CHD gene regulatory network (GRN) imbalances are unknown. Here, we define transcriptional consequences of reduced dosage of the CHD transcription factor, TBX5, in individual cells during cardiomyocyte differentiation from human induced pluripotent stem cells (iPSCs). We discovered highly sensitive dysregulation of TBX5-dependent pathways-including lineage decisions and genes associated with heart development, cardiomyocyte function, and CHD genetics-in discrete subpopulations of cardiomyocytes. Spatial transcriptomic mapping revealed chamber-restricted expression for many TBX5-sensitive transcripts. GRN analysis indicated that cardiac network stability, including vulnerable CHD-linked nodes, is sensitive to TBX5 dosage. A GRN-predicted genetic interaction between Tbx5 and Mef2c, manifesting as ventricular septation defects, was validated in mice. These results demonstrate exquisite and diverse sensitivity to TBX5 dosage in heterogeneous subsets of iPSC-derived cardiomyocytes and predicts candidate GRNs for human CHDs, with implications for quantitative transcriptional regulation in disease.

Cholecystokinin peptide signaling is regulated by a TBX5-MEF2 axis in the heart

The procholecystokinin (proCCK) gene encodes a secreted peptide known to regulate the digestive, endocrine, and nervous systems. Though recently proposed as a biomarker for heart dysfunction, its physiological role in both the embryonic and adult heart is poorly understood, and there are no reports of tissue-specific regulators of cholecystokinin signaling in the heart or other tissues. In the present study, mRNA of proCCK was observed in cardiac tissues during mouse embryonic development, establishing proCCK as an early marker of differentiated cardiomyocytes which is later restricted to anatomical subdomains of the neonatal heart. Three-dimensional analysis of the expression of proCCK and CCKAR/CCKBR receptors was performed using in situ hybridization and optical projection tomography, illustrating chamber-specific expression patterns in the postnatal heart. Transcription factor motif analyses indicated developmental cardiac transcription factors TBX5 and MEF2C as upstream regulators of proCCK, and this regulatory activity was confirmed in reporter gene assays. proCCK mRNA levels were also measured in the infarcted heart and in response to cyclic mechanical stretch and endothelin-1, indicating dynamic transcriptional regulation which might be leveraged for improved biomarker development. Functional analyses of exogenous cholecystokinin octapeptide (CCK-8) administration were performed in differentiating mouse embryonic stem cells (mESCs), and the results suggest that CCK-8 does not act as a differentiation modulator of cardiomyocyte subtypes. Collectively, these findings indicate that proCCK is regulated at the transcriptional level by TBX5-MEF2 and neurohormonal signaling, informing use of proCCK as a biomarker and future strategies for upstream manipulation of cholecystokinin signaling in the heart and other tissues.

Anterior lateral plate mesoderm gives rise to multiple tissues and requires tbx5a function in left-right asymmetry, migration dynamics, and cell specification of late-addition cardiac cells

In this study, we elucidate a single cell resolution fate map in the zebrafish in a sub-section of the anterior Lateral Plate Mesoderm (aLPM) at 18 hpf. Our results show that this tissue is not organized into segregated regions but gives rise to intermingled pericardial sac, peritoneum, pharyngeal arch and cardiac precursors. We further report upon asymmetrical contributions of lateral aLPM-derived heart precursors-specifically that twice as many heart precursors arise from the right side versus the left side of the embryo. Cell tracking analyses and large-scale cell labeling of the lateral aLPM corroborate these differences and show that the observed asymmetries are dependent upon Tbx5a expression. Previously, it was shown that cardiac looping was affected in Tbx5a knock-down and knock-out zebrafish (Garrity et al., 2002; Parrie et al., 2013); our present data also implicate tbx5a function in cell specification, establishment and maintenance of cardiac left-right asymmetry.

Mass spectrometry-based protein-protein interaction networks for the study of human diseases

A better understanding of the molecular mechanisms underlying disease is key for expediting the development of novel therapeutic interventions. Disease mechanisms are often mediated by interactions between proteins. Insights into the physical rewiring of protein-protein interactions in response to mutations, pathological conditions, or pathogen infection can advance our understanding of disease etiology, progression, and pathogenesis and can lead to the identification of potential druggable targets. Advances in quantitative mass spectrometry (MS)-based approaches have allowed unbiased mapping of these disease-mediated changes in protein-protein interactions on a global scale. Here, we review MS techniques that have been instrumental for the identification of protein-protein interactions at a system-level, and we discuss the challenges associated with these methodologies as well as novel MS advancements that aim to address these challenges. An overview of examples from diverse disease contexts illustrates the potential of MS-based protein-protein interaction mapping approaches for revealing disease mechanisms, pinpointing new therapeutic targets, and eventually moving toward personalized applications.

TBX5-encoded T-box transcription factor 5 variant T223M is associated with long QT syndrome and pediatric sudden cardiac death

Long QT syndrome (LQTS) is a genetic disease resulting in a prolonged QT interval on a resting electrocardiogram, predisposing affected individuals to polymorphic ventricular tachycardia and sudden death. Although a number of genes have been implicated in this disease, nearly one in four individuals exhibiting the LQTS phenotype are genotype-negative. Whole-exome sequencing identified a missense T223M variant in TBX5 that cosegregates with prolonged QT interval in a family with otherwise genotype-negative LQTS and sudden death. The TBX5-T223M variant was absent among large ostensibly healthy populations (gnomAD) and predicted to be pathogenic by in silico modeling based on Panther, PolyPhen-2, Provean, SIFT, SNAP2, and PredictSNP prediction tools. The variant was located in a highly conserved region of TBX5 predicted to be part of the DNA-binding interface. A luciferase assay identified a 57.5% reduction in the ability of TBX5-T223M to drive expression at the atrial natriuretic factor promotor compared to wildtype TBX5 in vitro. We conclude that the variant is pathogenic in this family, and we put TBX5 forward as a disease susceptibility allele for genotype-negative LQTS. The identification of this familial variant may serve as a basis for the identification of previously unknown mechanisms of LQTS with broader implications for cardiac electrophysiology.

miR-182-5p is an evolutionarily conserved Tbx5 effector that impacts cardiac development and electrical activity in zebrafish

To dissect the TBX5 regulatory circuit, we focused on microRNAs (miRNAs) that collectively contribute to make TBX5 a pivotal cardiac regulator. We profiled miRNAs in hearts isolated from wild-type, CRE, Tbx5lox/+and Tbx5del/+ mice using a Next Generation Sequencing (NGS) approach. TBX5 deficiency in cardiomyocytes increased the expression of the miR-183 cluster family that is controlled by Kruppel-like factor 4, a transcription factor repressed by TBX5. MiR-182-5p, the most highly expressed miRNA of this family, was functionally analyzed in zebrafish. Transient overexpression of miR-182-5p affected heart morphology, calcium handling and the onset of arrhythmias as detected by ECG tracings. Accordingly, several calcium channel proteins identified as putative miR-182-5p targets were downregulated in miR-182-5p overexpressing hearts. In stable zebrafish transgenic lines, we demonstrated that selective miRNA-182-5p upregulation contributes to arrhythmias. Moreover, cardiac-specific down-regulation of miR-182-5p rescued cardiac defects in a zebrafish model of Holt-Oram syndrome. In conclusion, miR-182-5p exerts an evolutionarily conserved role as a TBX5 effector in the onset of cardiac propensity for arrhythmia, and constitutes a relevant target for mediating the relationship between TBX5, arrhythmia and heart development.

Genetic and Epigenetic Control of Heart Development

A transcriptional program implemented by transcription factors and epigenetic regulators governs cardiac development and disease. Mutations in these factors are important causes of congenital heart disease. Here, we review selected recent advances in our understanding of the transcriptional and epigenetic control of heart development, including determinants of cardiac transcription factor chromatin occupancy, the gene regulatory network that regulates atrial septation, the chromatin landscape and cardiac gene regulation, and the role of Brg/Brahma-associated factor (BAF), nucleosome remodeling and histone deacetylation (NuRD), and Polycomb epigenetic regulatory complexes in heart development.

Hypermethylation-mediated down-regulation of lncRNA TBX5-AS1:2 in Tetralogy of Fallot inhibits cell proliferation by reducing TBX5 expression

Tetralogy of Fallot (TOF) is the most common complex congenital heart disease (CHD) with uncertain cause. Although long non‐coding RNAs (lncRNAs) have been implicated in heart development and several CHDs, their role in TOF is not well understood. This study aimed to investigate how dysregulated lncRNAs contribute to TOF. Using Gene Expression Omnibus data mining, bioinformatics analysis and clinical heart tissue sample detecting, we identified a novel antisense lncRNA TBX5‐AS1:2 with unknown function that was significantly down‐regulated in injured cardiac tissues from TOF patients. LncRNA TBX5‐AS1:2 was mainly located in the nucleus of the human embryonic kidney 293 (HEK293T) cells and formed an RNA‐RNA double‐stranded structure in the overlapping region with its sense mRNA T‐box transcription factor 5 (TBX5), which is an important regulator in heart development. Knock‐down of lncRNA TBX5‐AS1:2 via promoter hypermethylation reduced TBX5 expression at both the mRNA and protein levels by affecting its mRNA stability through RNA‐RNA interaction. Moreover, lncRNA TBX5‐AS1:2 knock‐down inhibited the proliferation of HEK293T cells. In conclusion, these results indicated that lncRNA TBX5‐AS1:2 may be involved in TOF by affecting cell proliferation by targeting TBX5.

TBX5 R264K acts as a modifier to develop dilated cardiomyopathy in mice independently of T-box pathway

Background

TBX5 is a transcription factor that has an important role in development of heart. TBX5 variants in the region encoding the T-box domain have been shown to cause cardiac defects, such as atrial septal defect or ventricular septal defect, while TBX5 variants have also been identified in a few cardiomyopathy patients and considered causative. We identified a TBX5 variant (c.791G>A, p.Arg264Lys), that is over-represented in cardiomyopathy patients. This variant is located outside of the T-box domain, and its pathogenicity has not been confirmed by functional analyses.

Objective

To investigate whether the TBX5 R264K is deleterious and could contribute to the pathogenesis of cardiomyopathy.

Methods and results

We developed mice expressing Tbx5 R264K. Mice homozygous for this variant displayed compensated dilated cardiomyopathy; mild decreased fractional shortening, dilatation of the left ventricle, left ventricular wall thinning and increased heart weight without major heart structural disorders. There was no difference in activation of the ANF promotor, a transcriptional target of Tbx5, compared to wild-type. However, analysis of RNA isolated from left ventricular samples showed significant increases in the expression of Acta1 in left ventricle with concomitant increases in the protein level of ACTA1.

Conclusions

Mice homozygous for Tbx5 R264K showed compensated dilated cardiomyopathy. Thus, TBX5 R264K may have a significant pathogenic role in some cardiomyopathy patients independently of T-box domain pathway.

Eomes and Brachyury control pluripotency exit and germ-layer segregation by changing the chromatin state

The first lineage specification of pluripotent mouse epiblast segregates neuroectoderm (NE) from mesoderm and definitive endoderm (ME) by mechanisms that are not well understood. Here we demonstrate that the induction of ME gene programs critically relies on the T-box transcription factors Eomesodermin (also known as Eomes) and Brachyury, which concomitantly repress pluripotency and NE gene programs. Cells deficient in these T-box transcription factors retain pluripotency and differentiate to NE lineages despite the presence of ME-inducing signals transforming growth factor β (TGF-β)/Nodal and Wnt. Pluripotency and NE gene networks are additionally repressed by ME factors downstream of T-box factor induction, demonstrating a redundancy in program regulation to safeguard mutually exclusive lineage specification. Analyses of chromatin revealed that accessibility of ME enhancers depends on T-box factor binding, whereas NE enhancers are accessible and already activation primed at pluripotency. This asymmetry of the chromatin landscape thus explains the default differentiation of pluripotent cells to NE in the absence of ME induction that depends on activating and repressive functions of Eomes and Brachyury.

HDAC4 and 5 repression of TBX5 is relieved by protein kinase D1

TBX5 is a T-box family transcription factor that regulates heart and forelimb development in vertebrates and functional deficiencies in this protein result in Holt-Oram syndrome. Recently, we have shown that acetylation of TBX5 potentiates its activity and is important for heart and limb development. Here we report that class II histone deacetylases HDAC4 and HDAC5 associate with TBX5 and repress its role in cardiac gene transcription. Both HDAC4 and HDAC5 deacetylate TBX5, which promotes its relocation to the cytoplasm and HDAC4 antagonizes the physical association and functional cooperation between TBX5 and MEF2C. We also show that protein kinase D1 (PRKD1) relieves the HDAC4/5-mediated repression of TBX5. Thus, this study reveals a novel interaction of HDAC4/5 and PRKD1 in the regulation of TBX5 transcriptional activity.

A reference map of murine cardiac transcription factor chromatin occupancy identifies dynamic and conserved enhancers

Mapping the chromatin occupancy of transcription factors (TFs) is a key step in deciphering developmental transcriptional programs. Here we use biotinylated knockin alleles of seven key cardiac TFs (GATA4, NKX2-5, MEF2A, MEF2C, SRF, TBX5, TEAD1) to sensitively and reproducibly map their genome-wide occupancy in the fetal and adult mouse heart. These maps show that TF occupancy is dynamic between developmental stages and that multiple TFs often collaboratively occupy the same chromatin region through indirect cooperativity. Multi-TF regions exhibit features of functional regulatory elements, including evolutionary conservation, chromatin accessibility, and activity in transcriptional enhancer assays. H3K27ac, a feature of many enhancers, incompletely overlaps multi-TF regions, and multi-TF regions lacking H3K27ac retain conservation and enhancer activity. TEAD1 is a core component of the cardiac transcriptional network, co-occupying cardiac regulatory regions and controlling cardiomyocyte-specific gene functions. Our study provides a resource for deciphering the cardiac transcriptional regulatory network and gaining insights into the molecular mechanisms governing heart development.

The Methyl-CpG-Binding Domain 2 and 3 Proteins and Formation of the Nucleosome Remodeling and Deacetylase Complex

The Nucleosome Remodeling and Deacetylase (NuRD) complex uniquely combines both deacetylase and remodeling enzymatic activities in a single macromolecular complex. The methyl-CpG-binding domain 2 and 3 (MBD2 and MBD3) proteins provide a critical structural link between the deacetylase and remodeling components, while MBD2 endows the complex with the ability to selectively recognize methylated DNA. Hence, NuRD combines three major arms of epigenetic gene regulation. Research over the past few decades has revealed much of the structural basis driving formation of this complex and started to uncover the functional roles of NuRD in epigenetic gene regulation. However, we have yet to fully understand the molecular and biophysical basis for methylation-dependent chromatin remodeling and transcription regulation by NuRD. In this review, we discuss the structural information currently available for the complex, the role MBD2 and MBD3 play in forming and recruiting the complex to methylated DNA, and the biological functions of NuRD.

Contribution of Mass Spectrometry-Based Proteomics to Discoveries in Developmental Biology

Understanding multicellular organism development from a molecular perspective is no small feat, yet this level of comprehension affords clinician-scientists the ability to identify root causes and mechanisms of congenital diseases. Inarguably, the maturation of molecular biology tools has significantly contributed to the identification of genetic loci that underlie normal and aberrant developmental programs. In combination with cell biology approaches, these tools have begun to elucidate the spatiotemporal expression and function of developmentally-regulated proteins. The emergence of quantitative mass spectrometry (MS) for biological applications has accelerated the pace at which these proteins can be functionally characterized, driving the construction of an increasingly detailed systems biology picture of developmental processes. Here, we review the quantitative MS-based proteomic technologies that have contributed significantly to understanding the role of proteome regulation in developmental processes. We provide a brief overview of these methodologies, focusing on their ability to provide precise and accurate proteome measurements. We then highlight the use of discovery-based and targeted mass spectrometry approaches in model systems to study cellular differentiation states, tissue phenotypes, and spatiotemporal subcellular organization. We also discuss the current application and future perspectives of MS proteomics to study PTM coordination and the role of protein complexes during development.

Conservation and divergence of protein pathways in the vertebrate heart

Heart disease is the leading cause of death in the western world. Attaining a mechanistic understanding of human heart development and homeostasis and the molecular basis of associated disease states relies on the use of animal models. Here, we present the cardiac proteomes of 4 model vertebrates with dual circulatory systems: the pig (Sus scrofa), the mouse (Mus musculus), and 2 frogs (Xenopus laevis and Xenopus tropicalis). Determination of which proteins and protein pathways are conserved and which have diverged within these species will aid in our ability to choose the appropriate models for determining protein function and to model human disease. We uncover mammalian- and amphibian-specific, as well as species-specific, enriched proteins and protein pathways. Among these, we find and validate an enrichment in cell-cycle-associated proteins within Xenopus laevis. To further investigate functional units within cardiac proteomes, we develop a computational approach to profile the abundance of protein complexes across species. Finally, we demonstrate the utility of these data sets for predicting appropriate model systems for studying given cardiac conditions by testing the role of Kielin/chordin-like protein (Kcp), a protein found as enriched in frog hearts compared to mammals. We establish that germ-line mutations in Kcp in Xenopus lead to valve defects and, ultimately, cardiac failure and death. Thus, integrating these findings with data on proteins responsible for cardiac disease should lead to the development of refined, species-specific models for protein function and disease states.

Spatiotemporal Gene Coexpression and Regulation in Mouse Cardiomyocytes of Early Cardiac Morphogenesis

Background Heart tube looping to form a 4-chambered heart is a critical stage of embryonic heart development, but the gene drivers and their regulatory targets have not been extensively characterized at the cell-type level. Methods and Results To study the interaction of signaling pathways, transcription factors (TFs), and genetic networks in the process, we constructed gene co-expression networks and identified gene modules highly activated in individual cardiomyocytes at multiple anatomical regions and developmental stages using previously published single-cell RNA-seq data. Function analyses of the modules uncovered major pathways important for spatiotemporal cardiomyocyte differentiation. Interestingly, about half of the pathways were highly active in cardiomyocytes at the outflow tract (OFT) and atrioventricular canal, including well-known pathways for cardiac development and many newly identified ones. We predicted that these OFT-atrioventricular canal pathways were regulated by a large number of TFs actively expressed at the OFT-atrioventricular canal cardiomyocytes, with the prediction supported by motif enrichment analysis, including 10 TFs that have not been previously associated with cardiac development (eg, Etv5, Rbpms, and Baz2b). Furthermore, we found that TF targets in the OFT-atrioventricular canal modules were most significantly enriched with genes associated with mouse heart developmental abnormalities and human congenital heart defects, in comparison with TF targets in other modules, consistent with the critical developmental roles of OFT. Conclusions By analyzing gene co-expression at single cardiomyocytes, our systematic study has uncovered many known and additional new important TFs and their regulated molecular signaling pathways that are spatiotemporally active during heart looping.

Single-Cell Transcriptomics Reveals Chemotaxis-Mediated Intraorgan Crosstalk During Cardiogenesis

RATIONALE

We hypothesized that the differentiation processes of cardiac progenitor cell (CP) from first and second heart fields (FHF and SHF) may undergo the unique instructive gene regulatory networks or signaling pathways, and the precise SHF progression is contingent on the FHF signaling developmental cues.

OBJECTIVE

We investigated how the intraorgan communications control sequential building of discrete anatomic regions of the heart at single-cell resolution.

METHODS AND RESULTS

By single-cell transcriptomic analysis of Nkx2-5 (NK2 homeobox 5) and Isl1 (ISL LIM homeobox 1) lineages at embryonic day 7.75, embryonic day 8.25, embryonic day 8.75, and embryonic day 9.25, we present a panoramic view of distinct CP differentiation hierarchies. Computational identifications of FHF- and SHF-CP descendants revealed that SHF differentiation toward cardiomyocytes underwent numerous step-like transitions, whereas earlier FHF progressed toward cardiomyocytes in a wave-like manner. Importantly, single-cell pairing analysis demonstrated that SHF-CPs were attracted to and expanded FHF-populated heart tube region through interlineage communications mediated by the chemotactic guidance (MIF [macrophage migration inhibitory factor]-CXCR2 [C-X-C motif chemokine receptor 2]). This finding was verified by pharmacological blockade of this chemotaxis in embryos manifesting limited SHF cell migration and contribution to the growth of the outflow tract and right ventricle but undetectable effects on the left ventricle or heart tube initiation. Genetic loss-of-function assay of Cxcr2 showed that the expression domain of CXCR4 was expanded predominantly at SHF. Furthermore, double knockout of Cxcr2/Cxcr4 exhibited defective SHF development, corroborating the redundant function. Mechanistically, NKX2-5 directly bound the Cxcr2 and Cxcr4 genomic loci and activated their transcription in SHF.

CONCLUSIONS

Collectively, we propose a model in which the chemotaxis-mediated intraorgan crosstalk spatiotemporally guides the successive process of positioning SHF-CP and promoting primary heart expansion and patterning upon FHF-derived heart tube initiation.

Single-Cell Transcriptomic Analyses of Cell Fate Transitions during Human Cardiac Reprogramming

Direct cellular reprogramming provides a powerful platform to study cell plasticity and dissect mechanisms underlying cell fate determination. Here, we report a single-cell transcriptomic study of human cardiac (hiCM) reprogramming that utilizes an analysis pipeline incorporating current data normalization methods, multiple trajectory prediction algorithms, and a cell fate index calculation we developed to measure reprogramming progression. These analyses revealed hiCM reprogramming-specific features and a decision point at which cells either embark on reprogramming or regress toward their original fibroblast state. In combination with functional screening, we found that immune-response-associated DNA methylation is required for hiCM induction and validated several downstream targets of reprogramming factors as necessary for productive hiCM reprograming. Collectively, this single-cell transcriptomics study provides detailed datasets that reveal molecular features underlying hiCM determination and rigorous analytical pipelines for predicting cell fate conversion.

Development and evolution of the metazoan heart

The mechanisms of the evolution and development of the heart in metazoans are highlighted, starting with the evolutionary origin of the contractile cell, supposedly the precursor of cardiomyocytes. The last eukaryotic common ancestor is likely a combination of several cellular organisms containing their specific metabolic pathways and genetic signaling networks. During evolution, these tool kits diversified. Shared parts of these conserved tool kits act in the development and functioning of pumping hearts and open or closed circulations in such diverse species as arthropods, mollusks, and chordates. The genetic tool kits became more complex by gene duplications, addition of epigenetic modifications, influence of environmental factors, incorporation of viral genomes, cardiac changes necessitated by air‐breathing, and many others. We evaluate mechanisms involved in mollusks in the formation of three separate hearts and in arthropods in the formation of a tubular heart. A tubular heart is also present in embryonic stages of chordates, providing the septated four‐chambered heart, in birds and mammals passing through stages with first and second heart fields. The four‐chambered heart permits the formation of high‐pressure systemic and low‐pressure pulmonary circulation in birds and mammals, allowing for high metabolic rates and maintenance of body temperature. Crocodiles also have a (nearly) separated circulation, but their resting temperature conforms with the environment. We argue that endothermic ancestors lost the capacity to elevate their body temperature during evolution, resulting in ectothermic modern crocodilians. Finally, a clinically relevant paragraph reviews the occurrence of congenital cardiac malformations in humans as derailments of signaling pathways during embryonic development.

The cardiac regulatory toolkit contains many factors including epigenetic, genetic, viral, hemodynamic, and environmental factors, but also transcriptional activators, repressors, duplicated genes, redundancies and dose‐dependancies.

Numerous toolkits regulate mechanisms including cell‐cell interactions, EMT, mitosis patterns, cell migration and differentiation and left/right sidedness involved in the development of endocardial cushions, looping, septum complexes, pharyngeal arch arteries, chamber and valve formation and conduction system.

Evolutionary development of the yolk sac circulation likely preceded the advent of endothermy in amniotes.

Parallel evolutionary traits regulate the development of contractile pumps in various taxa often in conjunction with the gut, lungs and excretory organs.

Navigating the non-coding genome in heart development and Congenital Heart Disease

Congenital Heart Disease (CHD) is characterised by a wide range of cardiac defects, from mild to life-threatening, which occur in babies worldwide. To date, there is no cure to CHD, however, progress in surgery has reduced its mortality allowing children affected by CHD to reach adulthood. In an effort to understand its genetic basis, several studies involving whole-genome sequencing (WGS) of patients with CHD have been undertaken and generated a great wealth of information. The majority of putative causative mutations identified in WGS studies fall into the non-coding part of the genome. Unfortunately, due to the lack of understanding of the function of these non-coding mutations, it is challenging to establish a causal link between the non-coding mutation and the disease. Thus, here we review the state-of-the-art approaches to interpret non-coding mutations in the context of CHD and address the following questions: What are the non-coding sequences important for cardiac function? Which technologies are used to identify them? Which resources are available to analyse them? What mutations are expected in these non-coding sequences? Learning from developmental process, what is their expected role in CHD?

T-box genes and retinoic acid signaling regulate the segregation of arterial and venous pole progenitor cells in the murine second heart field

The arterial and venous poles of the mammalian heart are hotspots of congenital heart defects (CHD) such as those observed in 22q11.2 deletion (or DiGeorge) and Holt-Oram syndromes. These regions of the heart are derived from late differentiating cardiac progenitor cells of the Second Heart Field (SHF) located in pharyngeal mesoderm contiguous with the elongating heart tube. The T-box transcription factor Tbx1, encoded by the major 22q11.2 deletion syndrome gene, regulates SHF addition to both cardiac poles from a common progenitor population. Despite the significance of this cellular addition the mechanisms regulating the deployment of common progenitor cells to alternate cardiac poles remain poorly understood. Here we demonstrate that Tbx5, mutated in Holt-Oram syndrome and essential for venous pole development, is activated in Tbx1 expressing cells in the posterior region of the SHF at early stages of heart tube elongation. A subset of the SHF transcriptional program, including Tbx1 expression, is subsequently downregulated in Tbx5 expressing cells, generating a transcriptional boundary between Tbx1-positive arterial pole and Tbx5-positive venous pole progenitor cell populations. We show that normal downregulation of the definitive arterial pole progenitor cell program in the posterior SHF is dependent on both Tbx1 and Tbx5. Furthermore, retinoic acid (RA) signaling is required for Tbx5 activation in Tbx1-positive cells and blocking RA signaling at the time of Tbx5 activation results in atrioventricular septal defects at fetal stages. Our results reveal sequential steps of cardiac progenitor cell patterning and provide mechanistic insights into the origin of common forms of CHD.

Proteomic-based approaches to cardiac development and disease

Congenital malformations, or structural birth defects, are now the leading cause of infant mortality in the United States and Europe (Dolk et al., 2010; Heron et al., 2009). Of the congenital malformations, congenital heart disease (CHD) is the most common (Dolk et al., 2010; Heron et al., 2009). Thus, a molecular understanding of heart development is an essential goal for improving clinical approaches to CHD. However, CHDs are commonly a result of genetic defects that manifest themselves in a spatial and temporal manner during the early stages of embryogenesis, leaving them mostly intractable to mass spectrometry-based analysis. Here, we describe the technologies and advancements in the field of mass spectrometry over the past few years that have begun to provide insights into the molecular and cellular basis of CHD and prospects for these types of approaches in the future.

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.

Holt-Oram Syndrome With Multiple Cardiac Abnormalities

Holt-Oram syndrome (HOS) is a rare monogenic disorder characterized by upper limb abnormalities, congenital heart defects and/or conduction abnormalities. It is determined by mutations of TBX5 gene and is inherited in an autosomal dominant manner. Penetrance is complete, but variable expressivity is present, which gives sometimes diagnostic difficulties. Our case is a young adult with a personal history of preaxial polydactyly operated in infancy, multiple cardiac malformations (atrial septal defect, bicuspid aortic valve, left ventricular non-compaction) and radiologic findings consistent with HOS. Family history is negative for HOS. In conclusion, we present a case of HOS diagnosed in the adult period to highlight the diagnostic problems for the proband and the family and the importance of an early diagnostic.

CHD4 and the NuRD complex directly control cardiac sarcomere formation

Cardiac development relies on proper cardiomyocyte differentiation, including expression and assembly of cell-type-specific actomyosin subunits into a functional cardiac sarcomere. Control of this process involves not only promoting expression of cardiac sarcomere subunits but also repressing expression of noncardiac myofibril paralogs. This level of transcriptional control requires broadly expressed multiprotein machines that modify and remodel the chromatin landscape to restrict transcription machinery access. Prominent among these is the nucleosome remodeling and deacetylase (NuRD) complex, which includes the catalytic core subunit CHD4. Here, we demonstrate that direct CHD4-mediated repression of skeletal and smooth muscle myofibril isoforms is required for normal cardiac sarcomere formation, function, and embryonic survival early in gestation. Through transcriptomic and genome-wide analyses of CHD4 localization, we identified unique CHD4 binding sites in smooth muscle myosin heavy chain, fast skeletal α-actin, and the fast skeletal troponin complex genes. We further demonstrate that in the absence of CHD4, cardiomyocytes in the developing heart form a hybrid muscle cell that contains cardiac, skeletal, and smooth muscle myofibril components. These misexpressed paralogs intercalate into the nascent cardiac sarcomere to disrupt sarcomere formation and cause impaired cardiac function in utero. These results demonstrate the genomic and physiological requirements for CHD4 in mammalian cardiac development.

Single-cell transcriptomics reconstructs fate conversion from fibroblast to cardiomyocyte

Direct lineage conversion offers a new strategy for tissue regeneration and disease modelling. Despite recent success in directly reprogramming fibroblasts into various cell types, the precise changes that occur as fibroblasts progressively convert to the target cell fates remain unclear. The inherent heterogeneity and asynchronous nature of the reprogramming process renders it difficult to study this process using bulk genomic techniques. Here we used single-cell RNA sequencing to overcome this limitation and analysed global transcriptome changes at early stages during the reprogramming of mouse fibroblasts into induced cardiomyocytes (iCMs). Using unsupervised dimensionality reduction and clustering algorithms, we identified molecularly distinct subpopulations of cells during reprogramming. We also constructed routes of iCM formation, and delineated the relationship between cell proliferation and iCM induction. Further analysis of global gene expression changes during reprogramming revealed unexpected downregulation of factors involved in mRNA processing and splicing. Detailed functional analysis of the top candidate splicing factor, Ptbp1, revealed that it is a critical barrier for the acquisition of cardiomyocyte-specific splicing patterns in fibroblasts. Concomitantly, Ptbp1 depletion promoted cardiac transcriptome acquisition and increased iCM reprogramming efficiency. Additional quantitative analysis of our dataset revealed a strong correlation between the expression of each reprogramming factor and the progress of individual cells through the reprogramming process, and led to the discovery of new surface markers for the enrichment of iCMs. In summary, our single-cell transcriptomics approaches enabled us to reconstruct the reprogramming trajectory and to uncover intermediate cell populations, gene pathways and regulators involved in iCM induction.

Formation of a TBX20-CASZ1 protein complex is protective against dilated cardiomyopathy and critical for cardiac homeostasis

By the age of 40, one in five adults without symptoms of cardiovascular disease are at risk for developing congestive heart failure. Within this population, dilated cardiomyopathy (DCM) remains one of the leading causes of disease and death, with nearly half of cases genetically determined. Though genetic and high throughput sequencing-based approaches have identified sporadic and inherited mutations in a multitude of genes implicated in cardiomyopathy, how combinations of asymptomatic mutations lead to cardiac failure remains a mystery. Since a number of studies have implicated mutations of the transcription factor TBX20 in congenital heart diseases, we investigated the underlying mechanisms, using an unbiased systems-based screen to identify novel, cardiac-specific binding partners. We demonstrated that TBX20 physically and genetically interacts with the essential transcription factor CASZ1. This interaction is required for survival, as mice heterozygous for both Tbx20 and Casz1 die post-natally as a result of DCM. A Tbx20 mutation associated with human familial DCM sterically interferes with the TBX20-CASZ1 interaction and provides a physical basis for how this human mutation disrupts normal cardiac function. Finally, we employed quantitative proteomic analyses to define the molecular pathways mis-regulated upon disruption of this novel complex. Collectively, our proteomic, biochemical, genetic, and structural studies suggest that the physical interaction between TBX20 and CASZ1 is required for cardiac homeostasis, and further, that reduction or loss of this critical interaction leads to DCM. This work provides strong evidence that DCM can be inherited through a digenic mechanism.

A molecular understanding of cardiomyocyte development is an essential goal for improving clinical approaches to CHD. While TBX20 is an essential transcription factor for heart development and its disease relevance is well established, many fundamental questions remain about the mechanism of TBX20 function. Principle among these is how TBX20 mutations associated with adult dilated cardiomyopathy circumvent (DCM) the essential embryonic requirement for TBX20 in heart development. Here we report using an integrated approach that TBX20 complexes with the cardiac transcription factor CASZ1 in vivo. We confirmed TBX20 and CASZ1 interact biochemically and genetically, and show mice heterozygous for both Tbx20 and Casz1 die, beginning at 4 to 8 weeks post birth, exhibiting hallmarks of DCM. Interestingly, the human mutant TBX20^F256I bypasses the early essential requirement for TBX20 but leads to DCM. We report here that TBX20^F256I disrupts the TBX20-CASZ1 interaction, ascribing clinical relevance to this protein complex. Further, by using quantitative proteomics we have identified the molecular pathways altered in TBX20-CASZ1-mediated DCM. Together, these results identify a novel interaction between TBX20 and CASZ1 that is essential for maintaining cardiac homeostasis and imply that DCM can be inherited through a digenic mechanism.