Sequential Binding of MEIS1 and NKX2-5 on the Popdc2 Gene: A Mechanism for Spatiotemporal Regulation of Enhancers during Cardiogenesis

Disparity of gene expression in coronary artery disease: insights from MEIS1, HIRA, and Myocardin

INTRODUCTION

Coronary artery disease (CAD) in young adults can have devastating consequences. The cardiac developmental gene MEIS1 plays important roles in vascular networks and heart development. This gene effects on the regeneration capacity of the heart. Considering role of MEIS1 in cardiac tissue development and the progression of myocardial infarction this study investigated the expression levels of the MEIS1, HIRA, and Myocardin genes in premature CAD patients compared to healthy subjects and evaluated the relationships between these genes and possible inflammatory factors.

METHODS AND RESULTS

The study conducted a case-control design involving 35 CAD patients and 35 healthy individuals. Peripheral blood mononuclear cells (PBMCs) were collected, and gene expression analysis was performed using real-time PCR. Compared with control group, the number of PBMCs in the CAD group exhibited greater MEIS1 and HIRA gene expression, with fold changes of 2.45 and 3.6. The expression of MEIS1 exhibited a negative correlation with IL-10 (r= -0.312) expression and positive correlation with Interleukin (IL)-6 (r = 0.415) and tumor necrosis factor (TNF)-α (r = 0.534) gene expression. Moreover, there was an inverse correlation between the gene expression of HIRA and that of IL-10 (r= -0.326), and a positive correlation was revealed between the expression of this gene and that of the IL-6 (r = 0.453) and TNF-α (r = 0.572) genes.

CONCLUSION

This research demonstrated a disparity in expression levels of MEIS1, HIRA, and Myocardin, between CAD and healthy subjects. The results showed that, MEIS1 and HIRA play significant roles in regulating the synthesis of proinflammatory cytokines, namely, TNF-α and IL-6.

TALE transcription factors: Cofactors no more

Exd/PBX, Hth/MEIS and PREP proteins belong to the TALE (three-amino-acid loop extension) superclass of transcription factors (TFs) with an atypical homedomain (HD). Originally discovered as "cofactors" to HOX proteins, revisiting their traditional role in light of genome-wide experiments reveals a strong and reproducible pattern of HOX and TALE co-occupancy across diverse embryonic tissues. While confirming that TALE increases HOX specificity and selectivity in vivo, this wider outlook also reveals novel aspects of HOX:TALE collaboration, namely that HOX TFs generally require pre-bound TALE factors to access their functional binding sites in vivo. In contrast to the restricted expression domains of HOX TFs, TALE factors are largely ubiquitous, and PBX and PREP are expressed at the earliest developmental stages. PBX and MEIS control development of many organs and tissues and their dysregulation is associated with congenital disease and cancer. Accordingly, many instances of TALE cooperation with non HOX TFs have been documented in various systems. The model that emerges from these studies is that TALE TFs create a permissive chromatin platform that is selected by tissue-restricted TFs for binding. In turn, HOX and other tissue-restricted TFs selectively convert a ubiquitous pool of low affinity TALE binding events into high confidence, tissue-restricted binding events associated with transcriptional activation. As a result, TALE:TF complexes are associated with active chromatin and domain/lineage-specific gene activity. TALE ubiquitous expression and broad genomic occupancy, as well as the increasing examples of TALE tissue-specific partners, reveal a universal and obligatory role for TALE in the control of tissue and lineage-specific transcriptional programs, beyond their initial discovery as HOX co-factors.

Atrioventricular node dysfunction in pressure overload-induced heart failure—Involvement of the immune system and transcriptomic remodelling

Heart failure is associated with atrioventricular (AV) node dysfunction, and AV node dysfunction in the setting of heart failure is associated with an increased risk of mortality and heart failure hospitalisation. This study aims to understand the causes of AV node dysfunction in heart failure by studying changes in the whole nodal transcriptome. The mouse transverse aortic constriction model of pressure overload-induced heart failure was studied; functional changes were assessed using electrocardiography and echocardiography and the transcriptome of the AV node was quantified using RNAseq. Heart failure was associated with a significant increase in the PR interval, indicating a slowing of AV node conduction and AV node dysfunction, and significant changes in 3,077 transcripts (5.6% of the transcriptome). Many systems were affected: transcripts supporting AV node conduction were downregulated and there were changes in transcripts identified by GWAS as determinants of the PR interval. In addition, there was evidence of remodelling of the sarcomere, a shift from fatty acid to glucose metabolism, remodelling of the extracellular matrix, and remodelling of the transcription and translation machinery. There was evidence of the causes of this widespread remodelling of the AV node: evidence of dysregulation of multiple intracellular signalling pathways, dysregulation of 109 protein kinases and 148 transcription factors, and an immune response with a proliferation of neutrophils, monocytes, macrophages and B lymphocytes and a dysregulation of 40 cytokines. In conclusion, inflammation and a widespread transcriptional remodelling of the AV node underlies AV node dysfunction in heart failure.

Quantitative Image Processing for Three-Dimensional Episcopic Images of Biological Structures: Current State and Future Directions

Episcopic imaging using techniques such as High Resolution Episcopic Microscopy (HREM) and its variants, allows biological samples to be visualized in three dimensions over a large field of view. Quantitative analysis of episcopic image data is undertaken using a range of methods. In this systematic review, we look at trends in quantitative analysis of episcopic images and discuss avenues for further research. Papers published between 2011 and 2022 were analyzed for details about quantitative analysis approaches, methods of image annotation and choice of image processing software. It is shown that quantitative processing is becoming more common in episcopic microscopy and that manual annotation is the predominant method of image analysis. Our meta-analysis highlights where tools and methods require further development in this field, and we discuss what this means for the future of quantitative episcopic imaging, as well as how annotation and quantification may be automated and standardized across the field.

Genome-wide survey identifies TNNI2 as a target of KLF7 that inhibits chicken adipogenesis via downregulating FABP4

Krüppel-like factor 7 (KLF7) negatively regulates adipocyte differentiation; however, the mechanism underlying its activity in mammals and birds remains poorly understood. To identify genome-wide KLF7-binding motifs in preadipocytes, we conducted a chromatin immunoprecipitation-sequencing analysis of immortalized chicken preadipocytes (ICP2), which revealed 11,063 specific binding sites. Intergenic binding site analysis showed that KLF7 regulates several novel factors whose functions in chicken and mammal adipogenesis are underexplored. We identified a novel regulator, troponin I2 (TNNI2), which is positively regulated by KLF7. TNNI2 is downregulated during preadipocyte differentiation and acts as an adipogenic repressor at least in part by repressing FABP4 promoter activity. In conclusion, we demonstrated that KLF7 functions through cis-regulation of TNNI2, which inhibits adipogenesis. Our findings not only provide the first genome-wide picture of KLF7 associations in preadipocytes but also identify a novel function of TNNI2.

Exploring the growth trait molecular markers in two sheep breeds based on Genome-wide association analysis

Growth traits are quantitative traits controlled by multiple micro-effect genes. we identified molecular markers related to sheep growth traits, which formed the basis of molecular breeding. In this study, we randomly selected 100 Qira Black sheep and 84 German Merino sheep for the blood collection the jugular vein to genotype by using the Illumina Ovine SNP 50K Bead Chip. quality control criteria for statistical analysis were: rejection detection rate < 90% and minimum allele frequency (MAF) < 5%. Then, we performed Genome-wide association studies (GWAS) on sheep body weight, body height, body length, and chest circumference using mixed linear models. After getting 55 SNPs with significant correlation, they were annotated by reference genome of Ovis aries genome (Oar_v4.0) and We obtained a total of 84 candidate genes associated with production traits (BMPR1B, HSD17B3, TMEM63C, etc.). We selected BMPR1B for population validation and found a correlation between the FecB locus and body weight traits. Therefore, this study not only supplements the existing knowledge of molecular markers of sheep growth traits, but also has important theoretical significance and reference value for the mining of functional genes of sheep growth traits.

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.

The Role of POPDC Proteins in Cardiac Pacemaking and Conduction

The Popeye domain-containing (POPDC) gene family, consisting of Popdc1 (also known as Bves), Popdc2, and Popdc3, encodes transmembrane proteins abundantly expressed in striated muscle. POPDC proteins have recently been identified as cAMP effector proteins and have been proposed to be part of the protein network involved in cAMP signaling. However, their exact biochemical activity is presently poorly understood. Loss-of-function mutations in animal models causes abnormalities in skeletal muscle regeneration, conduction, and heart rate adaptation after stress. Likewise, patients carrying missense or nonsense mutations in POPDC genes have been associated with cardiac arrhythmias and limb-girdle muscular dystrophy. In this review, we introduce the POPDC protein family, and describe their structure function, and role in cAMP signaling. Furthermore, the pathological phenotypes observed in zebrafish and mouse models and the clinical and molecular pathologies in patients carrying POPDC mutations are described.

Synergistic roles of homeodomain proteins UNC-62 homothorax and MLS-2 HMX/NKX in the specification of olfactory neurons in Caenorhabditis elegans

General identity of the Caenorhabditis elegans AWC olfactory neuron pair is specified by the OTX/OTD transcription factor CEH-36 and the HMG-box transcription factor SOX-2, followed by asymmetrical differentiation of the pair into two distinct subtypes, default AWCOFF and induced AWCON, through a stochastic signaling event. The HMX/NKX transcription factor MLS-2 regulates the expression of ceh-36 to specify general AWC identity. However, general AWC identity is lost in only one of the two AWC cells in the majority of mls-2 null mutants displaying defective general AWC identity, suggesting that additional transcription factors have a partially overlapping role with MLS-2 in the specification of general AWC identity. Here, we identify a role of unc-62, encoding a homothorax/Meis/TALE homeodomain protein, in the specification of general AWC identity. As in mls-2 null mutants, unc-62 null mutants showed an incomplete penetrance in loss of general AWC identity. However, unc-62; mls-2 double mutants display a nearly complete penetrance of identity loss in both AWC cells. Thus, unc-62 and mls-2 have a partially overlapping function in the specification of general AWC identity. Furthermore, our genetic results suggest that mls-2 and unc-62 act cell autonomously in promoting the AWCON subtype. Together, our findings reveal the sequential roles of the unc-62 and mls-2 pair in AWC development, specification of general AWC identity in early embryogenesis, and asymmetric differentiation of AWC subtypes in late embryogenesis.

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.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

A role of Hey2 transcription factor for right ventricle development through regulation of Tbx2-Mycn pathway during cardiac morphogenesis

A basic helix-loop-helix transcription factor Hey2 is expressed in the ventricular myocardium and endocardium of mouse embryos, and Hey2 null mice die perinatally showing ventricular septal defect, dysplastic tricuspid valve and hypoplastic right ventricle. In order to understand region-specific roles of Hey2 during cardiac morphogenesis, we generated Hey2 conditional knockout (cKO) mice using Mef2c-AHF-Cre, which was active in the anterior part of the second heart field and the right ventricle and outflow tract of the heart. Hey2 cKO neonates reproduced three anomalies commonly observed in Hey2 null mice. An earliest morphological defect was the lack of right ventricular extension along the apico-basal axis at midgestational stages. Underdevelopment of the right ventricle was present in all cKO neonates including those without apparent atresia of right-sided atrioventricular connection. RNA sequencing analysis of cKO embryos identified that the gene expression of a non-chamber T-box factor Tbx2 was ectopically induced in the chamber myocardium of the right ventricle. Consistently, mRNA expression of the Mycn transcription factor, which was a cell cycle regulator transcriptionally repressed by Tbx2, was down regulated, and the number of S-phase cells was significantly decreased in the right ventricle of cKO heart. These results suggest that Hey2 plays an important role in right ventricle development during cardiac morphogenesis, at least in part, through mitigating Tbx2-dependent inhibition of Mycn expression.

The M310T mutation in the GATA4 gene is a novel pathogenic target of the familial atrial septal defect

BACKGROUND

Although most cases of atrial septal defect (ASD) are sporadic, familial cases have been reported, which may be caused by mutation of transcription factor GATA binding protein 4 (GATA4). Herein we combined whole-exome sequencing and bioinformatics strategies to identify a novel mutation in GATA4 accounting for the etiology in a Chinese family with ASD.

METHODS

We identified kindred spanning 3 generations in which 3 of 12 (25.0%) individuals had ASD. Punctilious records for the subjects included complete physical examination, transthoracic echocardiography, electrocardiograph and surgical confirming. Whole-exome capture and high-throughput sequencing were performed on the proband III.1. Sanger sequencing was used to validate the candidate variants, and segregation analyses were performed in the family members.

RESULTS

Direct sequencing of GATA4 from the genomic DNA of family members identified a T-to-C transition at nucleotide 929 in exon 5 that predicted a methionine to threonine substitution at codon 310 (M310T) in the nuclear localization signal (NLS) region. Two affected members (II.2 and III.3) and the proband (III.1) who was recognized as a carrier exhibited this mutation, whereas the other unaffected family members or control individuals did not. More importantly, the mutation GATA4 (c.T929C: p.M310T) has not been reported previously in either familial or sporadic cases of congenital heart defects (CHD).

CONCLUSIONS

We identified for the first time a novel M310T mutation in the GATA4 gene that is located in the NLS region and leads to family ASD with arrhythmias. However, the mechanism by which this pathogenic mutation contributes to the development of heart defect and tachyarrhythmias remains to be ascertained.

An update on genetic variants of the NKX2-5

NKX2-5 is a homeodomain transcription factor that plays a crucial role in heart development. It is the first gene where a single genetic variant (GV) was found to be associated with congenital heart diseases in humans. In this study, we carried out a comprehensive survey of NKX2-5 GVs to build a unified, curated, and updated compilation of all available GVs. We retrieved a total of 1,380 unique GVs. From these, 970 had information on their frequency in the general population and 143 have been linked to pathogenic phenotypes in humans. In vitro effect was ascertained for 38 GVs. The homeodomain had the biggest cluster of pathogenic variants in the protein: 49 GVs in 60 residues, 23 in its third α-helix, where 11 missense variants may affect protein-DNA interaction or the hydrophobic core. We also pinpointed the likely location of pathogenic GVs in four linear motifs. These analyses allowed us to assign a putative explanation for the effect of 90 GVs. This study pointed to reliable pathogenicity for GVs in helix 3 of the homeodomain and may broaden the scope of functional and structural studies that can be done to better understand the effect of GVs in NKX2-5 function.

Development of Small Molecule MEIS Inhibitors that modulate HSC activity

Meis1, which belongs to TALE-type class of homeobox gene family, appeared as one of the key regulators of hematopoietic stem cell (HSC) self-renewal and a potential therapeutical target. However, small molecule inhibitors of MEIS1 remained unknown. This led us to develop inhibitors of MEIS1 that could modulate HSC activity. To this end, we have established a library of relevant homeobox family inhibitors and developed a high-throughput in silico screening strategy against homeodomain of MEIS proteins using the AutoDock Vina and PaDEL-ADV platform. We have screened over a million druggable small molecules in silico and selected putative MEIS inhibitors (MEISi) with no predicted cytotoxicity or cardiotoxicity. This was followed by in vitro validation of putative MEIS inhibitors using MEIS dependent luciferase reporter assays and analysis in the ex vivo HSC assays. We have shown that small molecules named MEISi-1 and MEISi-2 significantly inhibit MEIS-luciferase reporters in vitro and induce murine (LSKCD34^l°^w cells) and human (CD34⁺, CD133⁺, and ALDH^hi cells) HSC self-renewal ex vivo. In addition, inhibition of MEIS proteins results in downregulation of Meis1 and MEIS1 target gene expression including Hif-1α, Hif-2α and HSC quiescence modulators. MEIS inhibitors are effective in vivo as evident by induced HSC content in the murine bone marrow and downregulation of expression of MEIS target genes. These studies warrant identification of first-in-class MEIS inhibitors as potential pharmaceuticals to be utilized in modulation of HSC activity and bone marrow transplantation studies.

O-linked β-N-acetylglucosamine transferase plays an essential role in heart development through regulating angiopoietin-1

O-linked N-acetylglucosamine (GlcNAc) transferase (OGT) is the only enzyme catalyzing O-GlcNAcylation. Although it has been shown that OGT plays an essential role in maintaining postnatal heart function, its role in heart development remains unknown. Here we showed that loss of OGT in early fetal cardiomyocytes led to multiple heart developmental defects including hypertrabeculation, biventricular dilation, atrial septal defects, ventricular septal defects, and defects in coronary vessel development. In addition, RNA sequencing revealed that Angiopoietin-1, required within cardiomyocytes for both myocardial and coronary vessel development, was dramatically downregulated in cardiomyocyte-specific OGT knockout mouse hearts. In conclusion, our data demonstrated that OGT plays an essential role in regulating heart development through activating expression of cardiomyocyte Angiopoietin-1.

Pbx4 limits heart size and fosters arch artery formation by partitioning second heart field progenitors and restricting proliferation

Vertebrate heart development requires the integration of temporally distinct differentiating progenitors. However, few signals are understood that restrict the size of the later-differentiating outflow tract (OFT). We show that improper specification and proliferation of second heart field (SHF) progenitors in zebrafish lazarus (lzr) mutants, which lack the transcription factor Pbx4, produces enlarged hearts owing to an increase in ventricular and smooth muscle cells. Specifically, Pbx4 initially promotes the partitioning of the SHF into anterior progenitors, which contribute to the OFT, and adjacent endothelial cell progenitors, which contribute to posterior pharyngeal arches. Subsequently, Pbx4 limits SHF progenitor (SHFP) proliferation. Single cell RNA sequencing of nkx2.5+ cells revealed previously unappreciated distinct differentiation states and progenitor subpopulations that normally reside within the SHF and arterial pole of the heart. Specifically, the transcriptional profiles of Pbx4-deficient nkx2.5+ SHFPs are less distinct and display characteristics of normally discrete proliferative progenitor and anterior, differentiated cardiomyocyte populations. Therefore, our data indicate that the generation of proper OFT size and arch arteries requires Pbx-dependent stratification of unique differentiation states to facilitate both homeotic-like transformations and limit progenitor production within the SHF.

Homeobox gene Meis1 modulates cardiovascular regeneration

Regeneration of cardiomyocytes, endothelial cells and vascular smooth muscle cells (three major lineages of cardiac tissues) following myocardial infarction is the critical step to recover the function of the damaged heart. Myeloid ecotropic viral integration site 1 (Meis1) was first discovered in leukemic mice in 1995 and its biological function has been extensively studied in leukemia, hematopoiesis, the embryonic pattering of body axis, eye development and various genetic diseases, such as restless leg syndrome. It was found that Meis1 is highly associated with Hox genes and their cofactors to exert its regulatory effects on multiple intracellular signaling pathways. Recently with the advent of bioinformatics, biochemical methods and advanced genetic engineering tools, new function of Meis1 has been found to be involved in the cell cycle regulation of cardiomyocytes and endothelial cells. For example, inhibition of Meis1 expression increases the proliferative capacity of neonatal mouse cardiomyocytes, whereas overexpression of Meis1 results in the reduction in the length of cardiomyocyte proliferative window. Interestingly, downregulation of one of the circular RNAs, which acts downstream of Meis1 in the cardiomyocytes, promotes angiogenesis and restores the myocardial blood supply, thus reinforcing better regeneration of the damaged heart. It appears that Meis1 may play double roles in modulating proliferation and regeneration of cardiomyocytes and endothelial cells post-myocardial infarction. In this review, we propose to summarize the major findings of Meis1 in modulating fetal development and adult abnormalities, especially focusing on the recent discoveries of Meis1 in controlling the fate of cardiomyocytes and endothelial cells.

Allele-specific NKX2-5 binding underlies multiple genetic associations with human electrocardiographic traits

The cardiac transcription factor (TF) gene NKX2-5 has been associated with electrocardiographic (EKG) traits through genome-wide association studies (GWASs), but the extent to which differential binding of NKX2-5 at common regulatory variants contributes to these traits has not yet been studied. We analyzed transcriptomic and epigenomic data from induced pluripotent stem cell-derived cardiomyocytes from seven related individuals, and identified ~2,000 single-nucleotide variants associated with allele-specific effects (ASE-SNVs) on NKX2-5 binding. NKX2-5 ASE-SNVs were enriched for altered TF motifs, for heart-specific expression quantitative trait loci and for EKG GWAS signals. Using fine-mapping combined with epigenomic data from induced pluripotent stem cell-derived cardiomyocytes, we prioritized candidate causal variants for EKG traits, many of which were NKX2-5 ASE-SNVs. Experimentally characterizing two NKX2-5 ASE-SNVs (rs3807989 and rs590041) showed that they modulate the expression of target genes via differential protein binding in cardiac cells, indicating that they are functional variants underlying EKG GWAS signals. Our results show that differential NKX2-5 binding at numerous regulatory variants across the genome contributes to EKG phenotypes.

Variational infinite heterogeneous mixture model for semi-supervised clustering of heart enhancers

Motivation

Mammalian genomes can contain thousands of enhancers but only a subset are actively driving gene expression in a given cellular context. Integrated genomic datasets can be harnessed to predict active enhancers. One challenge in integration of large genomic datasets is the increasing heterogeneity: continuous, binary and discrete features may all be relevant. Coupled with the typically small numbers of training examples, semi-supervised approaches for heterogeneous data are needed; however, current enhancer prediction methods are not designed to handle heterogeneous data in the semi-supervised paradigm.

Results

We implemented a Dirichlet Process Heterogeneous Mixture model that infers Gaussian, Bernoulli and Poisson distributions over features. We derived a novel variational inference algorithm to handle semi-supervised learning tasks where certain observations are forced to cluster together. We applied this model to enhancer candidates in mouse heart tissues based on heterogeneous features. We constrained a small number of known active enhancers to appear in the same cluster, and 47 additional regions clustered with them. Many of these are located near heart-specific genes. The model also predicted 1176 active promoters, suggesting that it can discover new enhancers and promoters.

Availability and implementation

We created the ‘dphmix’ Python package: https://pypi.org/project/dphmix/.

Supplementary information

Supplementary data are available at Bioinformatics online.

High-Resolution Episcopic Microscopy (HREM): Looking Back on 13 Years of Successful Generation of Digital Volume Data of Organic Material for 3D Visualisation and 3D Display

High-resolution episcopic microscopy (HREM) is an imaging technique that permits the simple and rapid generation of three-dimensional (3D) digital volume data of histologically embedded and physically sectioned specimens. The data can be immediately used for high-detail 3D analysis of a broad variety of organic materials with all modern methods of 3D visualisation and display. Since its first description in 2006, HREM has been adopted as a method for exploring organic specimens in many fields of science, and it has recruited a slowly but steadily growing user community. This review aims to briefly introduce the basic principles of HREM data generation and to provide an overview of scientific publications that have been published in the last 13 years involving HREM imaging. The studies to which we refer describe technical details and specimen-specific protocols, and provide examples of the successful use of HREM in biological, biomedical and medical research. Finally, the limitations, potentials and anticipated further improvements are briefly outlined.

Furin, a transcriptional target of NKX2-5, has an essential role in heart development and function

The homeodomain transcription factor NKX2-5 is known to be essential for both normal heart development and for heart function. But little is yet known about the identities of its downstream effectors or their function during differentiation of cardiac progenitor cells (CPCs). We have used transgenic analysis and CRISPR-mediated ablation to identify a cardiac enhancer of the Furin gene. The Furin gene, encoding a proprotein convertase, is directly repressed by NKX2-5. Deletion of Furin in CPCs is embryonic lethal, with mutant hearts showing a range of abnormalities in the outflow tract. Those defects are associated with a reduction in proliferation and premature differentiation of the CPCs. Deletion of Furin in differentiated cardiomyocytes results in viable adult mutant mice showing an elongation of the PR interval, a phenotype that is consistent with the phenotype of mice and human mutant for Nkx2-5. Our results show that Furin mediate some aspects of Nkx2-5 function in the heart.

Dynamics of cardiomyocyte transcriptome and chromatin landscape demarcates key events of heart development

Organogenesis involves dynamic regulation of gene transcription and complex multipathway interactions. Despite our knowledge of key factors regulating various steps of heart morphogenesis, considerable challenges in understanding its mechanism still exist because little is known about their downstream targets and interactive regulatory network. To better understand transcriptional regulatory mechanism driving heart development and the consequences of its disruption in vivo, we performed time-series analyses of the transcriptome and genome-wide chromatin accessibility in isolated cardiomyocytes (CMs) from wild-type zebrafish embryos at developmental stages corresponding to heart tube morphogenesis, looping, and maturation. We identified genetic regulatory modules driving crucial events of heart development that contained key cardiac TFs and are associated with open chromatin regions enriched for DNA sequence motifs belonging to the family of the corresponding TFs. Loss of function of cardiac TFs Gata5, Tbx5a, and Hand2 affected the cardiac regulatory networks and caused global changes in chromatin accessibility profile, indicating their role in heart development. Among regions with differential chromatin accessibility in mutants were highly conserved noncoding elements that represent putative enhancers driving heart development. The most prominent gene expression changes, which correlated with chromatin accessibility modifications within their proximal promoter regions, occurred between heart tube morphogenesis and looping, and were associated with metabolic shift and hematopoietic/cardiac fate switch during CM maturation. Our results revealed the dynamic regulatory landscape throughout heart development and identified interactive molecular networks driving key events of heart morphogenesis.

Homeobox Genes and Homeodomain Proteins: New Insights into Cardiac Development, Degeneration and Regeneration

Cardiovascular diseases are the most common cause of human death in the developing world. Extensive evidence indicates that various toxic environmental factors and unhealthy lifestyle choices contribute to the risk, incidence and severity of cardiovascular diseases. Alterations in the genetic level of myocardium affects normal heart development and initiates pathological processes leading to various types of cardiac diseases. Homeobox genes are a large and highly specialized family of closely related genes that direct the formation of body structure, including cardiac development. Homeobox genes encode homeodomain proteins that function as transcription factors with characteristic structures that allow them to bind to DNA, regulate gene expression and subsequently control the proper physiological function of cells, tissues and organs. Mutations in homeobox genes are rare and usually lethal with evident alterations in cardiac function at or soon after the birth. Our understanding of homeobox gene family expression and function has expanded significantly during the recent years. However, the involvement of homeobox genes in the development of human and animal cardiac tissue requires further investigation. The phenotype of human congenital heart defects unveils only some aspects of human heart development. Therefore, mouse models are often used to gain a better understanding of human heart function, pathology and regeneration. In this review, we have focused on the role of homeobox genes in the development and pathology of human heart as potential tools for the future development of targeted regenerative strategies for various heart malfunctions.

HIRA directly targets the enhancers of selected cardiac transcription factors during in vitro differentiation of mouse embryonic stem cells

HIRA is a histone chaperone known to modulate gene expression through the deposition of H3.3. Conditional knockout of Hira in embryonic mouse hearts leads to cardiac septal defects. Loss of function mutation in HIRA, together with other chromatin modifiers, was found in patients with congenital heart diseases. However, the effects of HIRA on gene expression at earlier stages of cardiogenic mesoderm differentiation have not yet been studied. Differentiation of mouse embryonic stem cells (mESCs) towards cardiomyocytes mimics some of these early events and is an accepted model of these early stages. We performed RNA-Seq and H3.3-HA ChIP-seq on both WT and Hira-null mESCs and early cardiomyocyte progenitors of both genotypes. Analysis of RNA-seq data showed differential down regulation of cardiovascular development-related genes in Hira-null cardiomyocytes compared to WT cardiomyocytes. We found HIRA-dependent H3.3 deposition at these genes. In particular, we observed that HIRA influenced directly the expression of the transcription factors Gata6, Meis1 and Tbx2, essential for cardiac septation, through H3.3 deposition. We therefore identified new direct targets of HIRA during cardiac differentiation.

Epigenome-wide association study reveals differential DNA methylation in individuals with a history of myocardial infarction

Cardiovascular diseases (CVDs) are the leading causes of death worldwide and represent a substantial economic burden on public health care systems. Epigenetic markers have potential as diagnostic markers before clinical symptoms have emerged, and as prognostic markers to inform the choice of clinical intervention. In this study, we performed an epigenome-wide association study (EWAS) for CVDs, to identify disease-specific alterations in DNA methylation. CpG methylation in blood samples from the northern Sweden population health study (NSPHS) (n = 729) was assayed on the Illumina Infinium HumanMethylation450 BeadChip. Individuals with a history of a CVD were identified in the cohort. It included individuals with hypertension (N = 147), myocardial infarction (MI) (N = 48), stroke (N = 27), thrombosis (N = 22) and cardiac arrhythmia (N = 5). Differential DNA methylation was observed at 211 CpG-sites in individuals with a history of MI (q <0.05). These sites represent 196 genes, of which 42 have been described in the scientific literature to be related to cardiac function, cardiovascular disease, cardiogenesis and recovery after ischemic injury. We have shown that individuals with a history of MI have a deviating pattern of DNA methylation at many genomic loci of which a large fraction has previously been linked to CVD. Our results highlight genes that might be important in the pathogenesis of MI or in recovery. In addition, the sites pointed out in this study can serve as candidates for further evaluation as potential biomarkers for MI.

Revised roles of ISL1 in a hES cell-based model of human heart chamber specification

The transcription factor ISL1 is thought to be key for conveying the multipotent and proliferative properties of cardiac precursor cells. Here, we investigate its function upon cardiac induction of human embryonic stem cells. We find that ISL1 does not stabilize the transient cardiac precursor cell state but rather serves to accelerate cardiomyocyte differentiation. Conversely, ISL1 depletion delays cardiac differentiation and respecifies nascent cardiomyocytes from a ventricular to an atrial identity. Mechanistic analyses integrate this unrecognized anti-atrial function of ISL1 with known and newly identified atrial inducers. In this revised view, ISL1 is antagonized by retinoic acid signaling via a novel player, MEIS2. Conversely, ISL1 competes with the retinoic acid pathway for prospective cardiomyocyte fate, which converges on the atrial specifier NR2F1. This study reveals a core regulatory network putatively controlling human heart chamber formation and also bears implications for the subtype-specific production of human cardiomyocytes with enhanced functional properties.

HIC2 regulates isoform switching during maturation of the cardiovascular system

Physiological changes during embryonic development are associated with changes in the isoform expression of both myocyte sarcomeric proteins and of erythrocyte haemoglobins. Cell type-specific isoform expression of these genes also occurs. Although these changes appear to be coordinated, it is unclear how changes in these disparate cell types may be linked. The transcription factor Hic2 is required for normal cardiac development and the mutant is embryonic lethal. Hic2 embryos exhibit precocious expression of the definitive-lineage haemoglobin Hbb-bt in circulating primitive erythrocytes and of foetal isoforms of cardiomyocyte genes (creatine kinase, Ckm, and eukaryotic elongation factor Eef1a2) as well as ectopic cardiac expression of fast-twitch skeletal muscle troponin isoforms. We propose that HIC2 regulates a switching event within both the contractile machinery of cardiomyocytes and the oxygen carrying systems during the developmental period where demands on cardiac loading change rapidly.

New kids on the block: The Popeye domain containing (POPDC) protein family acting as a novel class of cAMP effector proteins in striated muscle

The cyclic 3′,5′-adenosine monophosphate (cAMP) signalling pathway constitutes an ancient signal transduction pathway present in prokaryotes and eukaryotes. Previously, it was thought that in eukaryotes three effector proteins mediate cAMP signalling, namely protein kinase A (PKA), exchange factor directly activated by cAMP (EPAC) and the cyclic-nucleotide gated channels. However, recently a novel family of cAMP effector proteins emerged and was termed the Popeye domain containing (POPDC) family, which consists of three members POPDC1, POPDC2 and POPDC3. POPDC proteins are transmembrane proteins, which are abundantly present in striated and smooth muscle cells. POPDC proteins bind cAMP with high affinity comparable to PKA. Presently, their biochemical activity is poorly understood. However, mutational analysis in animal models as well as the disease phenotype observed in patients carrying missense mutations suggests that POPDC proteins are acting by modulating membrane trafficking of interacting proteins. In this review, we will describe the current knowledge about this gene family and also outline the apparent gaps in our understanding of their role in cAMP signalling and beyond.

•

Popeye domain containing (POPDC) proteins are novel class of cAMP effector proteins.

•

POPDC proteins control membrane trafficking of interacting proteins.

•

POPDC proteins play a role in cardiac pacemaking and atrioventricular conduction.

•

Mutations of POPDC genes are causing muscular dystrophy.

Gene Expression Networks in the Murine Pulmonary Myocardium Provide Insight into the Pathobiology of Atrial Fibrillation

The pulmonary myocardium is a muscular coat surrounding the pulmonary and caval veins. Although its definitive physiological function is unknown, it may have a pathological role as the source of ectopic beats initiating atrial fibrillation. How the pulmonary myocardium gains pacemaker function is not clearly defined, although recent evidence indicates that changed transcriptional gene expression networks are at fault. The gene expression profile of this distinct cell type in situ was examined to investigate underlying molecular events that might contribute to atrial fibrillation. Via systems genetics, a whole-lung transcriptome data set from the BXD recombinant inbred mouse resource was analyzed, uncovering a pulmonary cardiomyocyte gene network of 24 transcripts, coordinately regulated by chromosome 1 and 2 loci. Promoter enrichment analysis and interrogation of publicly available ChIP-seq data suggested that transcription of this gene network may be regulated by the concerted activity of NKX2-5, serum response factor, myocyte enhancer factor 2, and also, at a post-transcriptional level, by RNA binding protein motif 20. Gene ontology terms indicate that this gene network overlaps with molecular markers of the stressed heart. Therefore, we propose that perturbed regulation of this gene network might lead to altered calcium handling, myocyte growth, and contractile force contributing to the aberrant electrophysiological properties observed in atrial fibrillation. We reveal novel molecular interactions and pathways representing possible therapeutic targets for atrial fibrillation. In addition, we highlight the utility of recombinant inbred mouse resources in detecting and characterizing gene expression networks of relatively small populations of cells that have a pathological significance.

Genome-Wide Temporal Profiling of Transcriptome and Open Chromatin of Early Cardiomyocyte Differentiation Derived From hiPSCs and hESCs

RATIONALE

Recent advances have improved our ability to generate cardiomyocytes from human induced pluripotent stem cells (hiPSCs) and human embryonic stem cells (hESCs). However, our understanding of the transcriptional regulatory networks underlying early stages (ie, from mesoderm to cardiac mesoderm) of cardiomyocyte differentiation remains limited.

OBJECTIVE

To characterize transcriptome and chromatin accessibility during early cardiomyocyte differentiation from hiPSCs and hESCs.

METHODS AND RESULTS

We profiled the temporal changes in transcriptome and chromatin accessibility at genome-wide levels during cardiomyocyte differentiation derived from 2 hiPSC lines and 2 hESC lines at 4 stages: pluripotent stem cells, mesoderm, cardiac mesoderm, and differentiated cardiomyocytes. Overall, RNA sequencing analysis revealed that transcriptomes during early cardiomyocyte differentiation were highly concordant between hiPSCs and hESCs, and clustering of 4 cell lines within each time point demonstrated that changes in genome-wide chromatin accessibility were similar across hiPSC and hESC cell lines. Weighted gene co-expression network analysis (WGCNA) identified several modules that were strongly correlated with different stages of cardiomyocyte differentiation. Several novel genes were identified with high weighted connectivity within modules and exhibited coexpression patterns with other genes, including noncoding RNA LINC01124 and uncharacterized RNA AK127400 in the module related to the mesoderm stage; E-box-binding homeobox 1 (ZEB1) in the module correlated with postcardiac mesoderm. We further demonstrated that ZEB1 is required for early cardiomyocyte differentiation. In addition, based on integrative analysis of both WGCNA and transcription factor motif enrichment analysis, we determined numerous transcription factors likely to play important roles at different stages during cardiomyocyte differentiation, such as T and eomesodermin (EOMES; mesoderm), lymphoid enhancer-binding factor 1 (LEF1) and mesoderm posterior BHLH transcription factor 1 (MESP1; from mesoderm to cardiac mesoderm), meis homeobox 1 (MEIS1) and GATA-binding protein 4 (GATA4) (postcardiac mesoderm), JUN and FOS families, and MEIS2 (cardiomyocyte).

CONCLUSIONS

Both hiPSCs and hESCs share similar transcriptional regulatory mechanisms underlying early cardiac differentiation, and our results have revealed transcriptional regulatory networks and new factors (eg, ZEB1) controlling early stages of cardiomyocyte differentiation.

Homeodomain proteins in action: similar DNA binding preferences, highly variable connectivity

Homeodomain proteins are evolutionary conserved proteins present in the entire eukaryote kingdom. They execute functions that are essential for life, both in developing and adult organisms. Most homeodomain proteins act as transcription factors and bind DNA to control the activity of other genes. In contrast to their similar DNA binding specificity, homeodomain proteins execute highly diverse and context-dependent functions. Several factors, including genome accessibility, DNA shape, combinatorial binding and the ability to interact with many transcriptional partners, diversify the activity of homeodomain proteins and culminate in the activation of highly dynamic, context-specific transcriptional programs. Clarifying how homeodomain transcription factors work is central to our understanding of development, disease and evolution.

Integrated live imaging and molecular profiling of embryoid bodies reveals a synchronized progression of early differentiation

Embryonic stem cells can spontaneously differentiate into cell types of all germ layers within embryoid bodies (EBs) in a highly variable manner. Whether there exists an intrinsic differentiation program common to all EBs is unknown. Here, we present a novel combination of high-throughput live two-photon imaging and gene expression profiling to study early differentiation dynamics spontaneously occurring within developing EBs. Onset timing of Brachyury-GFP was highly variable across EBs, while the spatial patterns as well as the dynamics of mesendodermal progression following onset were remarkably similar. We therefore defined a 'developmental clock' using the Brachyury-GFP signal onset timing. Mapping snapshot gene expression measurements to this clock revealed their temporal trends, indicating that loss of pluripotency, formation of primitive streak and mesodermal lineage progression are synchronized in EBs. Exogenous activation of Wnt or BMP signaling accelerated the intrinsic clock. CHIR down-regulated Wnt3, allowing insights into dependency mechanisms between canonical Wnt signaling and multiple genes. Our findings reveal a developmental clock characteristic of an early differentiation program common to all EBs, further establishing them as an in vitro developmental model.

HIRA Is Required for Heart Development and Directly Regulates Tnni2 and Tnnt3

Chromatin remodelling is essential for cardiac development. Interestingly, the role of histone chaperones has not been investigated in this regard. HIRA is a member of the HUCA (HIRA/UBN1/CABIN1/ASF1a) complex that deposits the variant histone H3.3 on chromatin independently of replication. Lack of HIRA has general effects on chromatin and gene expression dynamics in embryonic stem cells and mouse oocytes. Here we describe the conditional ablation of Hira in the cardiogenic mesoderm of mice. We observed surface oedema, ventricular and atrial septal defects and embryonic lethality. We identified dysregulation of a subset of cardiac genes, notably upregulation of troponins Tnni2 and Tnnt3, involved in cardiac contractility and decreased expression of Epha3, a gene necessary for the fusion of the muscular ventricular septum and the atrioventricular cushions. We found that HIRA binds GAGA rich DNA loci in the embryonic heart, and in particular a previously described enhancer of Tnni2/Tnnt3 (TTe) bound by the transcription factor NKX2.5. HIRA-dependent H3.3 enrichment was observed at the TTe in embryonic stem cells (ESC) differentiated toward cardiomyocytes in vitro. Thus, we show here that HIRA has locus-specific effects on gene expression and that histone chaperone activity is vital for normal heart development, impinging on pathways regulated by an established cardiac transcription factor.

Spatiotemporal regulation of enhancers during cardiogenesis

With the advance in chromatin immunoprecipitation followed by high-throughput sequencing, there has been a dramatic increase in our understanding of distal enhancer function. In the developing heart, the identification and characterisation of such enhancers have deepened our knowledge of the mechanisms of transcriptional regulation that drives cardiac differentiation. With next-generation sequencing techniques becoming widely accessible, the quantity of data describing the genome-wide distribution of cardiac-specific transcription factor and chromatin modifiers has rapidly increased and it is now becoming clear that the usage of enhancers is highly dynamic and complex, both during the development and in the adult. The identification of those enhancers has revealed new insights into the transcriptional mechanisms of how tissue-specific gene expression patterns are established, maintained, and change dynamically during development and upon physiological stress.

Heartbreak hotel: a convergence in cardiac regeneration

In February 2016, The Company of Biologists hosted an intimate gathering of leading international researchers at the forefront of experimental cardiovascular regeneration, with its emphasis on ‘Transdifferentiation and Tissue Plasticity in Cardiovascular Rejuvenation’. As I review here, participants at the workshop revealed how understanding cardiac growth and lineage decisions at their most fundamental level has transformed the strategies in hand that presently energize the prospects for human heart repair.

Hox Genes in Cardiovascular Development and Diseases

Congenital heart defects (CHD) are the leading cause of death in the first year of life. Over the past 20 years, much effort has been focused on unraveling the genetic bases of CHD. In particular, studies in human genetics coupled with those of model organisms have provided valuable insights into the gene regulatory networks underlying CHD pathogenesis. Hox genes encode transcription factors that are required for the patterning of the anterior–posterior axis in the embryo. In this review, we focus on the emerging role of anteriorly expressed Hox genes (Hoxa1, Hoxb1, and Hoxa3) in cardiac development, specifically their contribution to patterning of cardiac progenitor cells and formation of the great arteries. Recent evidence regarding the cooperative regulation of heart development by Hox proteins with members of the TALE-class of homeodomain proteins such as Pbx and Meis transcription factors is also discussed. These findings are highly relevant to human pathologies as they pinpoint new genes that increase susceptibility to cardiac anomalies and provide novel mechanistic insights into CHD.