Cardioids reveal self-organizing principles of human cardiogenesis.

Organoids capable of forming tissue-like structures have transformed our ability to model human development and disease. With the notable exception of the human heart, lineage-specific self-organizing organoids have been reported for all major organs. Here, we established self-organizing cardioids from human pluripotent stem cells that intrinsically specify, pattern, and morph into chamber-like structures containing a cavity. Cardioid complexity can be controlled by signaling that instructs the separation of cardiomyocyte and endothelial layers and by directing epicardial spreading, inward migration, and differentiation. We find that cavity morphogenesis is governed by a mesodermal WNT-BMP signaling axis and requires its target HAND1, a transcription factor linked to developmental heart chamber defects. Upon cryoinjury, cardioids initiated a cell-type-dependent accumulation of extracellular matrix, an early hallmark of both regeneration and heart disease. Thus, human cardioids represent a powerful platform to mechanistically dissect self-organization, congenital heart defects and serve as a foundation for future translational research.

Complex Interdependence Regulates Heterotypic Transcription Factor Distribution and Coordinates Cardiogenesis.

Transcription factors (TFs) are thought to function with partners to achieve specificity and precise quantitative outputs. In the developing heart, heterotypic TF interactions, such as between the T-box TF TBX5 and the homeodomain TF NKX2-5, have been proposed as a mechanism for human congenital heart defects. We report extensive and complex interdependent genomic occupancy of TBX5, NKX2-5, and the zinc finger TF GATA4 coordinately controlling cardiac gene expression, differentiation, and morphogenesis. Interdependent binding serves not only to co-regulate gene expression but also to prevent TFs from distributing to ectopic loci and activate lineage-inappropriate genes. We define preferential motif arrangements for TBX5 and NKX2-5 cooperative binding sites, supported at the atomic level by their co-crystal structure bound to DNA, revealing a direct interaction between the two factors and induced DNA bending. Complex interdependent binding mechanisms reveal tightly regulated TF genomic distribution and define a combinatorial logic for heterotypic TF regulation of differentiation.

Development of a 3D atlas of the embryonic pancreas for topological and quantitative analysis of heterologous cell interactions.

Generating comprehensive image maps, while preserving spatial three-dimensional (3D) context, is essential in order to locate and assess quantitatively specific cellular features and cell-cell interactions during organ development. Despite recent advances in 3D imaging approaches, our current knowledge of the spatial organization of distinct cell types in the embryonic pancreatic tissue is still largely based on two-dimensional histological sections. Here, we present a light-sheet fluorescence microscopy approach to image the pancreas in three dimensions and map tissue interactions at key time points in the mouse embryo. We demonstrate the utility of the approach by providing volumetric data, 3D distribution of three main cellular components (epithelial, mesenchymal and endothelial cells) within the developing pancreas, and quantification of their relative cellular abundance within the tissue. Interestingly, our 3D images show that endocrine cells are constantly and increasingly in contact with endothelial cells forming small vessels, whereas the interactions with mesenchymal cells decrease over time. These findings suggest distinct cell-cell interaction requirements for early endocrine cell specification and late differentiation. Lastly, we combine our image data in an open-source online repository (referred to as the Pancreas Embryonic Cell Atlas).

Prioritizing cardiovascular disease-associated variants altering NKX2-5 and TBX5 binding through an integrative computational approach.

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies have mapped >90% of CVD-associated variants within the noncoding genome, which can alter the function of regulatory proteins, such as transcription factors (TFs). However, due to the overwhelming number of single-nucleotide polymorphisms (SNPs) (>500,000) in genome-wide association studies, prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1535 CVD-associated SNPs within TF footprints and putative cardiac enhancers plus 14,218 variants in linkage disequilibrium with genotype-dependent gene expression in cardiac tissues. Using ChIP-seq data from two cardiac TFs (NKX2-5 and TBX5) in human-induced pluripotent stem cell-derived cardiomyocytes, we trained a large-scale gapped k-mer SVM model to identify CVD-associated SNPs that altered NKX2-5 and TBX5 binding. The model was tested by scoring human heart TF genomic footprints within putative enhancers and measuring in vitro binding through electrophoretic mobility shift assay. Five variants predicted to alter NKX2-5 (rs59310144, rs6715570, and rs61872084) and TBX5 (rs7612445 and rs7790964) binding were prioritized for in vitro validation based on the magnitude of the predicted change in binding and are in cardiac tissue eQTLs. All five variants altered NKX2-5 and TBX5 DNA binding. We present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro analysis.

An Anterior Second Heart Field Enhancer Regulates the Gene Regulatory Network of the Cardiac Outflow Tract.

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

Class I and II Histone Deacetylase Inhibitors as Therapeutic Modulators of Dilated Cardiac Tissue-Derived Mesenchymal Stem/Stromal Cells.

The prevalence of dilated cardiomyopathy (DCM) is increasing globally, highlighting the need for innovative therapeutic approaches to prevent its onset. In this study, we examined the energetic and epigenetic distinctions between dilated and non-dilated human myocardium-derived mesenchymal stem/stromal cells (hmMSCs) and assessed the effects of class I and II HDAC inhibitors (HDACi) on these cells and their cardiomyogenic differentiation. Cells were isolated from myocardium biopsies using explant outgrowth methods. Mitochondrial and histone deacetylase activities, ATP levels, cardiac transcription factors, and structural proteins were assessed using flow cytometry, PCR, chemiluminescence, Western blotting, and immunohistochemistry. The data suggest that the tested HDAC inhibitors improved acetylation and enhanced the energetic status of both types of cells, with significant effects observed in dilated myocardium-derived hmMSCs. Additionally, the HDAC inhibitors activated the cardiac transcription factors Nkx2-5, HOPX, GATA4, and Mef2C, and upregulated structural proteins such as cardiac troponin T and alpha cardiac actin at both the protein and gene levels. In conclusion, our findings suggest that HDACi may serve as potential modulators of the energetic status and cardiomyogenic differentiation of human heart hmMSCs. This avenue of exploration could broaden the search for novel therapeutic interventions for dilated cardiomyopathy, ultimately leading to improvements in heart function.

The G4 resolvase RHAU regulates ventricular trabeculation and compaction through transcriptional and post-transcriptional mechanisms.

The G-quadruplex (G4) resolvase RNA helicase associated with AU-rich element (RHAU) possesses the ability to unwind G4 structures in both DNA and RNA molecules. Previously, we revealed that RHAU plays a critical role in embryonic heart development and postnatal heart function through modulating mRNA translation and stability. However, whether RHAU functions to resolve DNA G4 in the regulation of cardiac physiology is still elusive. Here, we identified a phenotype of noncompaction cardiomyopathy in cardiomyocyte-specific Rhau deletion mice, including such symptoms as spongiform cardiomyopathy, heart dilation, and death at young ages. We also observed reduced cardiomyocyte proliferation and advanced sarcomere maturation in Rhau mutant mice. Further studies demonstrated that RHAU regulates the expression levels of several genes associated with ventricular trabeculation and compaction, including the Nkx2-5 and Hey2 that encode cardiac transcription factors of NKX2-5 and Hey2, and the myosin heavy chain 7 (Myh7) whose protein product is MYH7. While RHAU modulates Nkx2-5 mRNA and Hey2 mRNA at the post-transcriptional level, we uncovered that RHAU facilitates the transcription of Myh7 through unwinding of the G4 structures in its promoter. These findings demonstrated that RHAU regulates ventricular chamber development through both transcriptional and post-transcriptional mechanisms. These results contribute to a knowledge base that will help to understand the pathogenesis of diseases such as noncompaction cardiomyopathy.

Application of 4-D ultrasound-derived regional strain and proteomics analysis in Nkx2-5-deficient male mice.

The comprehensive characterization of cardiac structure and function is critical to better understanding various murine models of cardiac disease. We demonstrate here a multimodal analysis approach using high-frequency four-dimensional ultrasound (4DUS) imaging and proteomics to explore the relationship between regional function and tissue composition in a murine model of metabolic cardiomyopathy (Nkx2-5183P/+). The presented 4DUS analysis outlines a novel approach to mapping both circumferential and longitudinal strain profiles through a standardized framework. We then demonstrate how this approach allows for spatiotemporal comparisons of cardiac function and improved localization of regional left ventricular dysfunction. Guided by observed trends in regional dysfunction, our targeted Ingenuity Pathway Analysis (IPA) results highlight metabolic dysregulation in the Nkx2-5183P/+ model, including altered mitochondrial function and energy metabolism (i.e., oxidative phosphorylation and fatty acid/lipid handling). Finally, we present a combined 4DUS-proteomics z-score-based analysis that highlights IPA canonical pathways showing strong linear relationships with 4DUS biomarkers of regional cardiac dysfunction. The presented multimodal analysis methods aim to help future studies more comprehensively assess regional structure-function relationships in other preclinical models of cardiomyopathy.NEW & NOTEWORTHY A multimodal approach using both four-dimensional ultrasound (4DUS) and regional proteomics can help enhance our investigations of murine cardiomyopathy models. We present unique 4DUS-derived strain maps that provide a framework for both cross-sectional and longitudinal analysis of spatiotemporal cardiac function. We further detail and demonstrate an innovative 4DUS-proteomics z-score-based linear regression method, aimed at characterizing relationships between regional cardiac dysfunction and underlying mechanisms of disease.

Association of maternal hypertension and diabetes with variants of the NKX2-5, LEFTY1 and LEFTY2 genes in children with congenital heart defects: a case-control study from Pakistani Population.

BACKGROUND: Globally, congenital heart defect (CHD) is the most common congenital malformation, responsible for higher morbidity and mortality in the pediatric population. It is a complex multifactorial disease influenced by gene-environment and gene-gene interactions. The current study was the first attempt to study these polymorphisms in common clinical phenotypes of CHD in Pakistan and the association between maternal hypertension and diabetes with single nucleotide polymorphisms (SNPs) in children. METHODS: A total of 376 subjects were recruited in this current case-control study. Six variants from three genes were analyzed by cost-effective multiplex PCR and genotyped by minisequencing. Statistical analysis was done by GraphPad prism and Haploview. The association of SNPs and CHD was determined using logistic regression. RESULTS: The risk allele frequency was higher in cases as compared to healthy subjects, but the results were not significant for rs703752. However, stratification analysis suggested that rs703752 was significantly associated with the tetralogy of Fallot. The rs2295418 was significantly associated with maternal hypertension (OR = 16.41, p = 0.003), while a weak association was present between maternal diabetes and rs360057 (p = 0.08). CONCLUSION: In conclusion, variants in transcriptional and signaling genes were associated with Pakistani pediatric CHD patients that showed varied susceptibility between different clinical phenotypes of CHD. In addition, this study was the first report regarding the significant association between maternal hypertension and the LEFTY2 gene variant.

Knockout of the Cardiac Transcription Factor NKX2-5 Results in Stem Cell-Derived Cardiac Cells with Typical Purkinje Cell-like Signal Transduction and Extracellular Matrix Formation.

The human heart controls blood flow, and therewith enables the adequate supply of oxygen and nutrients to the body. The correct function of the heart is coordinated by the interplay of different cardiac cell types. Thereby, one can distinguish between cells of the working myocardium, the pace-making cells in the sinoatrial node (SAN) and the conduction system cells in the AV-node, the His-bundle or the Purkinje fibres. Tissue-engineering approaches aim to generate hiPSC-derived cardiac tissues for disease modelling and therapeutic usage with a significant improvement in the differentiation quality of myocardium and pace-making cells. The differentiation of cells with cardiac conduction system properties is still challenging, and the produced cell mass and quality is poor. Here, we describe the generation of cardiac cells with properties of the cardiac conduction system, called conduction system-like cells (CSLC). As a primary approach, we introduced a CrispR-Cas9-directed knockout of the NKX2-5 gene in hiPSC. NKX2-5-deficient hiPSC showed altered connexin expression patterns characteristic for the cardiac conduction system with strong connexin 40 and connexin 43 expression and suppressed connexin 45 expression. Application of differentiation protocols for ventricular- or SAN-like cells could not reverse this connexin expression pattern, indicating a stable regulation by NKX2-5 on connexin expression. The contraction behaviour of the hiPSC-derived CSLCs was compared to hiPSC-derived ventricular- and SAN-like cells. We found that the contraction speed of CSLCs resembled the expected contraction rate of human conduction system cells. Overall contraction was reduced in differentiated cells derived from NKX2-5 knockout hiPSC. Comparative transcriptomic data suggest a specification of the cardiac subtype of CSLC that is distinctly different from ventricular or pacemaker-like cells with reduced myocardial gene expression and enhanced extracellular matrix formation for improved electrical insulation. In summary, knockout of NKX2-5 in hiPSC leads to enhanced differentiation of cells with cardiac conduction system features, including connexin expression and contraction behaviour.

Mef2c factors are required for early but not late addition of cardiomyocytes to the ventricle.

During heart formation, the heart grows and undergoes dramatic morphogenesis to achieve efficient embryonic function. Both in fish and amniotes, much of the growth occurring after initial heart tube formation arises from second heart field (SHF)-derived progenitor cell addition to the arterial pole, allowing chamber formation. In zebrafish, this process has been extensively studied during embryonic life, but it is unclear how larval cardiac growth occurs beyond 3 days post-fertilisation (dpf). By quantifying zebrafish myocardial growth using live imaging of GFP-labelled myocardium we show that the heart grows extensively between 3 and 5 dpf. Using methods to assess cell division, cellular development timing assay and Kaede photoconversion, we demonstrate that proliferation, CM addition, and hypertrophy contribute to ventricle growth. Mechanistically, we show that reduction in Mef2c activity (mef2ca+/-;mef2cb-/-), downstream or in parallel with Nkx2.5 and upstream of Ltbp3, prevents some CM addition and differentiation, resulting in a significantly smaller ventricle by 3 dpf. After 3 dpf, however, CM addition in mef2ca+/-;mef2cb-/- mutants recovers to a normal pace, and the heart size gap between mutants and their siblings diminishes into adulthood. Thus, as in mice, there is an early time window when SHF contribution to the myocardium is particularly sensitive to loss of Mef2c activity.

Role of Cdx factors in early mesodermal fate decisions.

Murine cardiac and hematopoietic progenitors are derived from Mesp1+ mesoderm. Cdx function impacts both yolk sac hematopoiesis and cardiogenesis in zebrafish, suggesting that Cdx family members regulate early mesoderm cell fate decisions. We found that Cdx2 occupies a number of transcription factor loci during embryogenesis, including key regulators of both cardiac and blood development, and that Cdx function is required for normal expression of the cardiogenic transcription factors Nkx2-5 and Tbx5 Furthermore, Cdx and Brg1, an ATPase subunit of the SWI/SNF chromatin remodeling complex, co-occupy a number of loci, suggesting that Cdx family members regulate target gene expression through alterations in chromatin architecture. Consistent with this, we demonstrate loss of Brg1 occupancy and altered chromatin structure at several cardiogenic genes in Cdx-null mutants. Finally, we provide evidence for an onset of Cdx2 expression at E6.5 coinciding with egression of cardiac progenitors from the primitive streak. Together, these findings suggest that Cdx functions in multi-potential mesoderm to direct early cell fate decisions through transcriptional regulation of several novel target genes, and provide further insight into a potential epigenetic mechanism by which Cdx influences target gene expression.

Gene-environment interaction impacts on heart development and embryo survival.

Congenital heart disease (CHD) is the most common type of birth defect. In recent years, research has focussed on identifying the genetic causes of CHD. However, only a minority of CHD cases can be attributed to single gene mutations. In addition, studies have identified different environmental stressors that promote CHD, but the additive effect of genetic susceptibility and environmental factors is poorly understood. In this context, we have investigated the effects of short-term gestational hypoxia on mouse embryos genetically predisposed to heart defects. Exposure of mouse embryos heterozygous for Tbx1 or Fgfr1/Fgfr2 to hypoxia in utero increased the incidence and severity of heart defects while Nkx2-5+/- embryos died within 2 days of hypoxic exposure. We identified the molecular consequences of the interaction between Nkx2-5 and short-term gestational hypoxia, which suggest that reduced Nkx2-5 expression and a prolonged hypoxia-inducible factor 1α response together precipitate embryo death. Our study provides insight into the causes of embryo loss and variable penetrance of monogenic CHD, and raises the possibility that cases of foetal death and CHD in humans could be caused by similar gene-environment interactions.

Nkx2-5 defines a subpopulation of pacemaker cells and is essential for the physiological function of the sinoatrial node in mice.

The sinoatrial node (SAN), the primary cardiac pacemaker, consists of a head domain and a junction/tail domain that exhibit different functional properties. However, the underlying molecular mechanism defining these two pacemaker domains remains elusive. Nkx2-5 is a key transcription factor essential for the formation of the working myocardium, but it was generally thought to be detrimental to SAN development. However, Nkx2-5 is expressed in the developing SAN junction, suggesting a role for Nkx2-5 in SAN junction development and function. In this study, we present unambiguous evidence that SAN junction cells exhibit unique action potential configurations intermediate to those manifested by the SAN head and the surrounding atrial cells, suggesting a specific role for the junction cells in impulse generation and in SAN-atrial exit conduction. Single-cell RNA-seq analyses support this concept. Although Nkx2-5 inactivation in the SAN junction did not cause a malformed SAN at birth, the mutant mice manifested sinus node dysfunction. Thus, Nkx2-5 defines a population of pacemaker cells in the transitional zone. Despite Nkx2-5 being dispensable for SAN morphogenesis during embryogenesis, its deletion hampers atrial activation by the pacemaker.

Nkx2-5(+)islet1(+) mesenchymal precursors generate distinct spleen stromal cell subsets and participate in restoring stromal network integrity.

Secondary lymphoid organ stromal cells comprise different subsets whose origins remain unknown. Herein, we exploit a genetic lineage-tracing approach to show that splenic fibroblastic reticular cells (FRCs), follicular dendritic cells (FDCs), marginal reticular cells (MRCs), and mural cells, but not endothelial cells, originate from embryonic mesenchymal progenitors of the Nkx2-5(+)Islet1(+) lineage. This lineage include embryonic mesenchymal cells with lymphoid tissue organizer (LTo) activity capable also of supporting ectopic lymphoid-like structures and a subset of resident spleen stromal cells that proliferate and regenerate the splenic stromal microenvironment following resolution of a viral infection. These findings identify progenitor cells that generate stromal diversity in spleen development and repair and suggest the existence of multipotent stromal progenitors in the adult spleen with regenerative capacity.

Epigenomic and Transcriptomic Dynamics During Human Heart Organogenesis.

RATIONALE: There is growing evidence that common variants and rare sequence alterations in regulatory sequences can result in birth defects or predisposition to disease. Congenital heart defects are the most common birth defect and have a clear genetic component, yet only a third of cases can be attributed to structural variation in the genome or a mutation in a gene. The remaining unknown cases could be caused by alterations in regulatory sequences. OBJECTIVE: Identify regulatory sequences and gene expression networks that are active during organogenesis of the human heart. Determine whether these sites and networks are enriched for disease-relevant genes and associated genetic variation. METHODS AND RESULTS: We characterized ChromHMM (chromatin state) and gene expression dynamics during human heart organogenesis. We profiled 7 histone modifications in embryonic hearts from each of 9 distinct Carnegie stages (13-14, 16-21, and 23), annotated chromatin states, and compared these maps to over 100 human tissues and cell types. We also generated RNA-sequencing data, performed differential expression, and constructed weighted gene coexpression networks. We identified 177 412 heart enhancers; 12 395 had not been previously annotated as strong enhancers. We identified 92% of all functionally validated heart-positive enhancers (n=281; 7.5× enrichment; P<2.2×10-16). Integration of these data demonstrated novel heart enhancers are enriched near genes expressed more strongly in cardiac tissue and are enriched for variants associated with ECG measures and atrial fibrillation. Our gene expression network analysis identified gene modules strongly enriched for heart-related functions, regulatory control by heart-specific enhancers, and putative disease genes. CONCLUSIONS: Well-connected hub genes with heart-specific expression targeted by embryonic heart-specific enhancers are likely disease candidates. Our functional annotations will allow for better interpretation of whole genome sequencing data in the large number of patients affected by congenital heart defects.

NKX2-5 Variant in Two Siblings with Thyroid Hemiagenesis.

Thyroid hemiagenesis (THA) is an inborn absence of one thyroid lobe of largely unknown etiopathogenesis. The aim of the study was to reveal genetic factors responsible for thyroid maldevelopment in two siblings with THA. None of the family members presented with congenital heart defect. The samples were subjected to whole-exome sequencing (WES) (Illumina, TruSeq Exome Enrichment Kit, San Diego, CA 92121, USA). An ultra-rare variant c.839C>T (p.Pro280Leu) in NKX2-5 gene (NM_004387.4) was identified in both affected children and an unaffected father. In the mother, the variant was not present. This variant is reported in population databases with 0.0000655 MAF (GnomAD v3, dbSNP rs761596254). The affected amino acid position is moderately conserved (positive scores in PhyloP: 1.364 and phastCons: 0.398). Functional prediction algorithms showed deleterious impact (dbNSFP v4.1, FATHMM, SIFT) or benign (CADD, PolyPhen-2, Mutation Assessor). According to ACMG criteria, variant is classified as having uncertain clinical significance. For the first time, NKX2-5 gene variants were found in two siblings with THA, providing evidence for its potential contribution to the pathogenesis of this type of thyroid dysgenesis. The presence of the variant in an unaffected parent, carrier of p.Pro280Leu variant, suggests potential contribution of yet unidentified additional factors determining the final penetrance and expression.

Conjugated activation of myocardial-specific transcription of Gja5 by a pair of Nkx2-5-Shox2 co-responsive elements.

The sinoatrial node (SAN) is the primary pacemaker in the heart. During cardiogenesis, Shox2 and Nkx2-5 are co-expressed in the junction domain of the SAN and regulate pacemaker cell fate through a Shox2-Nkx2-5 antagonism. Cx40 is a marker of working myocardium and an Nkx2-5 transcriptional output antagonized by Shox2, but the underlying regulatory mechanisms remain elusive. Here we characterized a bona fide myocardial-specific Gja5 (coding gene of Cx40) distal enhancer consisting of a pair of Nkx2-5 and Shox2 co-bound elements in the regulatory region of Gja5. Transgenic reporter assays revealed that neither element alone, but the conjugation of both elements together, drives myocardial-specific transcription. Genetic analyses confirmed that the activation of this enhancer depends on Nkx2-5 but is inhibited by Shox2 in vivo, and its presence is essential for Gja5 expression in the myocardium but not the endothelial cells of the heart. Furthermore, chromatin conformation analysis showed an Nkx2-5-dependent loop formation between these two elements and the Gja5 promoter in vivo, indicating that Nkx2-5 bridges the conjugated activation of this enhancer by pairing the two elements to the Gja5 promoter.

Axenfeld-Rieger syndrome-associated mutants of the transcription factor FOXC1 abnormally regulate NKX2-5 in model zebrafish embryos.

FOXC1 is a member of the forkhead family of transcription factors, and whose function is poorly understood. A variety of FOXC1 mutants have been identified in patients diagnosed with the autosomal dominant disease Axenfeld-Rieger syndrome, which is mainly characterized by abnormal development of the eyes, particularly those who also have accompanying congenital heart defects (CHD). However, the role of FOXC1 in CHD, and how these mutations might impact FOXC1 function, remains elusive. Our previous work provided one clue to possible function, demonstrating that zebrafish foxc1a, an orthologue of human FOXC1 essential for heart development, directly regulates the expression of nkx2.5, encoding a transcriptional regulator of cardiac progenitor cells. Abnormal expression of Nkx2-5 leads to CHD in mice and is also associated with CHD patients. Whether this link extends to the human system, however, requires investigation. In this study, we demonstrate that FOXC1 does regulate human NKX2-5 expression in a dose-dependent manner via direct binding to its proximal promoter. A comparison of FOXC1 mutant function in the rat cardiac cell line H9c2 and zebrafish embryos suggested that the zebrafish embryos might serve as a more representative model system than the H9c2 cells. Finally, we noted that three of the Axenfeld-Rieger syndrome FOXC1 mutations tested increased, whereas a fourth repressed the expression of NKX2-5 These results imply that mutant FOXC1s might play etiological roles in CHD by abnormally regulating NKX2-5 in the patients. And zebrafish embryos can serve as a useful in vivo platform for rapidly evaluating disease-causing roles of mutated genes.

AMPKß1 and AMPKß2 define an isoform-specific gene signature in human pluripotent stem cells, differentially mediating cardiac lineage specification.

AMP-activated protein kinase (AMPK) is a key regulator of energy metabolism that phosphorylates a wide range of proteins to maintain cellular homeostasis. AMPK consists of three subunits: α, ß, and γ. AMPKα and ß are encoded by two genes, the γ subunit by three genes, all of which are expressed in a tissue-specific manner. It is not fully understood, whether individual isoforms have different functions. Using RNA-Seq technology, we provide evidence that the loss of AMPKß1 and AMPKß2 lead to different gene expression profiles in human induced pluripotent stem cells (hiPSCs), indicating isoform-specific function. The knockout of AMPKß2 was associated with a higher number of differentially regulated genes than the deletion of AMPKß1, suggesting that AMPKß2 has a more comprehensive impact on the transcriptome. Bioinformatics analysis identified cell differentiation as one biological function being specifically associated with AMPKß2. Correspondingly, the two isoforms differentially affected lineage decision toward a cardiac cell fate. Although the lack of PRKAB1 impacted differentiation into cardiomyocytes only at late stages of cardiac maturation, the availability of PRKAB2 was indispensable for mesoderm specification as shown by gene expression analysis and histochemical staining for cardiac lineage markers such as cTnT, GATA4, and NKX2.5. Ultimately, the lack of AMPKß1 impairs, whereas deficiency of AMPKß2 abrogates differentiation into cardiomyocytes. Finally, we demonstrate that AMPK affects cellular physiology by engaging in the regulation of hiPSC transcription in an isoform-specific manner, providing the basis for further investigations elucidating the role of dedicated AMPK subunits in the modulation of gene expression.