The generation and validation of two NKX2-5 fluorescent reporter human embryonic stem cell lines: UMANe002-A-1 and UMANe002-A-2

The transcription factor NKX2-5 is a highly conserved master regulator of heart development which is widely expressed in cardiac progenitors and cardiomyocytes. Fluorescent reporters of NKX2-5 that minimally perturb normal protein expression can enable the identification, quantification and isolation of NKX2-5-expressing cells in a normal physiological state. Here we report the generation of two new hESC lines with eGFP inserted upstream (5') or downstream (3') of NKX2-5, linked by a cleavable T2A peptide. These complementary reporters produce a robust fluorescent signal in cardiac cells and have wide utility particularly for research on developmental biology and disease modelling.

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.

The multi-lineage transcription factor ISL1 controls cardiomyocyte cell fate through interaction with NKX2.5

Congenital heart disease often arises from perturbations of transcription factors (TFs) that guide cardiac development. ISLET1 (ISL1) is a TF that influences early cardiac cell fate, as well as differentiation of other cell types including motor neuron progenitors (MNPs) and pancreatic islet cells. While lineage specificity of ISL1 function is likely achieved through combinatorial interactions, its essential cardiac interacting partners are unknown. By assaying ISL1 genomic occupancy in human induced pluripotent stem cell-derived cardiac progenitors (CPs) or MNPs and leveraging the deep learning approach BPNet, we identified motifs of other TFs that predicted ISL1 occupancy in each lineage, with NKX2.5 and GATA motifs being most closely associated to ISL1 in CPs. Experimentally, nearly two-thirds of ISL1-bound loci were co-occupied by NKX2.5 and/or GATA4. Removal of NKX2.5 from CPs led to widespread ISL1 redistribution, and overexpression of NKX2.5 in MNPs led to ISL1 occupancy of CP-specific loci. These results reveal how ISL1 guides lineage choices through a combinatorial code that dictates genomic occupancy and transcription.

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.

Nr2f1a maintains atrial nkx2.5 expression to repress pacemaker identity within venous atrial cardiomyocytes of zebrafish

Maintenance of cardiomyocyte identity is vital for normal heart development and function. However, our understanding of cardiomyocyte plasticity remains incomplete. Here, we show that sustained expression of the zebrafish transcription factor Nr2f1a prevents the progressive acquisition of ventricular cardiomyocyte (VC) and pacemaker cardiomyocyte (PC) identities within distinct regions of the atrium. Transcriptomic analysis of flow-sorted atrial cardiomyocytes (ACs) from nr2f1a mutant zebrafish embryos showed increased VC marker gene expression and altered expression of core PC regulatory genes, including decreased expression of nkx2.5, a critical repressor of PC differentiation. At the arterial (outflow) pole of the atrium in nr2f1a mutants, cardiomyocytes resolve to VC identity within the expanded atrioventricular canal. However, at the venous (inflow) pole of the atrium, there is a progressive wave of AC transdifferentiation into PCs across the atrium toward the arterial pole. Restoring Nkx2.5 is sufficient to repress PC marker identity in nr2f1a mutant atria and analysis of chromatin accessibility identified an Nr2f1a-dependent nkx2.5 enhancer expressed in the atrial myocardium directly adjacent to PCs. CRISPR/Cas9-mediated deletion of the putative nkx2.5 enhancer leads to a loss of Nkx2.5-expressing ACs and expansion of a PC reporter, supporting that Nr2f1a limits PC differentiation within venous ACs via maintaining nkx2.5 expression. The Nr2f-dependent maintenance of AC identity within discrete atrial compartments may provide insights into the molecular etiology of concurrent structural congenital heart defects and associated arrhythmias.

Disease-associated non-coding variants alter NKX2-5 DNA-binding affinity

Genome-wide association studies (GWAS) have mapped over 90 % of disease- or trait-associated variants within the non-coding genome, like cis-regulatory elements (CREs). Non-coding single nucleotide polymorphisms (SNPs) are genomic variants that can change how DNA-binding regulatory proteins, like transcription factors (TFs), interact with the genome and regulate gene expression. NKX2-5 is a TF essential for proper heart development, and mutations affecting its function have been associated with congenital heart diseases (CHDs). However, establishing a causal mechanism between non-coding genomic variants and human disease remains challenging. To address this challenge, we identified 8475 SNPs predicted to alter NKX2-5 DNA-binding using a position weight matrix (PWM)-based predictive model. Five variants were prioritized for in vitro validation; four of them are associated with traits and diseases that impact cardiovascular health. The impact of these variants on NKX2-5 binding was evaluated with electrophoretic mobility shift assay (EMSA) using purified recombinant NKX2-5 homeodomain. Binding curves were constructed to determine changes in binding between variant and reference alleles. Variants rs7350789, rs7719885, rs747334, and rs3892630 increased binding affinity, whereas rs61216514 decreased binding by NKX2-5 when compared to the reference genome. Our findings suggest that differential TF-DNA binding affinity can be key in establishing a causal mechanism of pathogenic variants.

Mutational Assessment in NKX2-5 and ACTC1 Genes in Patients with Congenital Cardiac Septal Defect (CCSD) from Ethnic Kashmiri Population

(1) Background globe. The etiology of CHDs is complex and involves both genetic and non-genetic factors. Although, significant progress has been made in deciphering the genetic components involved in CHDs, recent reports have revealed that mutations in Nk2 homeobox5 (NKX2-5) and actin alpha cardiac muscle1 (ACTC1) genes play a key role in CHDs such as atrial and ventricular septum defects. Therefore, the present study evaluates the role of key hotspot mutations in NKX2-5 and ACTC1 genes of congenital cardiac septal defect (CCSD) in ethnic Kashmiri population. (2) Methods: A total of 112 confirmed CHD patients were included in the current study, of which 30 patients were evaluated for mutational analysis for hotspot mutations of NKX2-5 and ACTC1 genes. The total genomic DNA was extracted from the samples (cardiac tissue/blood) and were subjected to amplification for NKX2-5 (exon 1 and 2), and ACTC1 (exon 2) genes by using PCR specific primers to analyze the hotspot mutations in respective exons. The amplified products obtained were sent to Macrogen Korea for sequencing by Sanger’s method. (3) Results: Our results confirmed that not a single mutation was found in either hotspot exon 1 and 2 of NKX2-5 and exon 2 of ACTC1 in the patients included in the current study. Interestingly, a novel synonymous nucleotide variation leading to G > C transversion (GCG > GCC) was found in exon 2 of NKX2-5 gene of CCSD patient. (4) Conclusions: The current findings demonstrated the role of NKX2-5 and ACTC1 in cardiac development. The study will provide an insight in understanding the genetic etiology and highlights the role of newly identified mutations in patients with CDS’s in ethnic Kashmiri population. In silico findings revealed amino acid changes, splice site variation and the creation of new site. Furthermore, the study warrants complete screening of genes involved in CCSDs.

Activation of Nkx2.5 transcriptional program is required for adult myocardial repair

The cardiac developmental network has been associated with myocardial regenerative potential. However, the embryonic signals triggered following injury have yet to be fully elucidated. Nkx2.5 is a key causative transcription factor associated with human congenital heart disease and one of the earliest markers of cardiac progenitors, thus it serves as a promising candidate. Here, we show that cardiac-specific RNA-sequencing studies reveal a disrupted embryonic transcriptional profile in the adult Nkx2.5 loss-of-function myocardium. nkx2.5-/- fish exhibit an impaired ability to recover following ventricular apex amputation with diminished dedifferentiation and proliferation. Complex network analyses illuminate that Nkx2.5 is required to provoke proteolytic pathways necessary for sarcomere disassembly and to mount a proliferative response for cardiomyocyte renewal. Moreover, Nkx2.5 targets embedded in these distinct gene regulatory modules coordinate appropriate, multi-faceted injury responses. Altogether, our findings support a previously unrecognized, Nkx2.5-dependent regenerative circuit that invokes myocardial cell cycle re-entry, proteolysis, and mitochondrial metabolism to ensure effective regeneration in the teleost heart.

Tissue-specific multi-omics analysis of atrial fibrillation

Genome-wide association studies (GWAS) for atrial fibrillation (AF) have uncovered numerous disease-associated variants. Their underlying molecular mechanisms, especially consequences for mRNA and protein expression remain largely elusive. Thus, refined multi-omics approaches are needed for deciphering the underlying molecular networks. Here, we integrate genomics, transcriptomics, and proteomics of human atrial tissue in a cross-sectional study to identify widespread effects of genetic variants on both transcript (cis-eQTL) and protein (cis-pQTL) abundance. We further establish a novel targeted trans-QTL approach based on polygenic risk scores to determine candidates for AF core genes. Using this approach, we identify two trans-eQTLs and five trans-pQTLs for AF GWAS hits, and elucidate the role of the transcription factor NKX2-5 as a link between the GWAS SNP rs9481842 and AF. Altogether, we present an integrative multi-omics method to uncover trans-acting networks in small datasets and provide a rich resource of atrial tissue-specific regulatory variants for transcript and protein levels for cardiovascular disease gene prioritization.

Stage-specific regulation of DNA methylation by TET enzymes during human cardiac differentiation

Changes in DNA methylation are associated with normal cardiogenesis, whereas altered methylation patterns can occur in congenital heart disease. Ten-eleven translocation (TET) enzymes oxidize 5-methylcytosine (5mC) and promote locus-specific DNA demethylation. Here, we characterize stage-specific methylation dynamics and the function of TETs during human cardiomyocyte differentiation. Human embryonic stem cells (hESCs) in which all three TET genes are inactivated fail to generate cardiomyocytes (CMs), with altered mesoderm patterning and defective cardiac progenitor specification. Genome-wide methylation analysis shows TET knockout causes promoter hypermethylation of genes encoding WNT inhibitors, leading to hyperactivated WNT signaling and defects in cardiac mesoderm patterning. TET activity is also needed to maintain hypomethylated status and expression of NKX2-5 for subsequent cardiac progenitor specification. Finally, loss of TETs causes a set of cardiac structural genes to fail to be demethylated at the cardiac progenitor stage. Our data demonstrate key roles for TET proteins in controlling methylation dynamics at sequential steps during human cardiac development.

The Genetic Landscape of Patent Foramen Ovale: A Systematic Review

Patent Foramen Ovale (PFO) is a common postnatal defect of cardiac atrial septation. A certain degree of familial aggregation has been reported. Animal studies suggest the involvement of the Notch pathway and other cardiac transcription factors (GATA4, TBX20, NKX2-5) in Foramen Ovale closure. This review evaluates the contribution of genetic alterations in PFO development. We systematically reviewed studies that assessed rare and common variants in subjects with PFO. The protocol was registered with PROSPERO and followed MOOSE guidelines. We systematically searched English studies reporting rates of variants in PFO subjects until the 30th of June 2021. Among 1231 studies, we included four studies: two of them assessed the NKX2-5 gene, the remaining reported variants of chromosome 4q25 and the GATA4 S377G variant, respectively. We did not find any variant associated with PFO, except for the rs2200733 variant of chromosome 4q25 in atrial fibrillation patients. Despite the scarceness of evidence so far, animal studies and other studies that did not fulfil the criteria to be included in the review indicate a robust genetic background in PFO. More research is needed on the genetic determinants of PFO.

Electrical stimulation applied during differentiation drives the hiPSC-CMs towards a mature cardiac conduction-like cells

Human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) resemble fetal cardiomyocytes and electrical stimulation (ES) has been explored to mature the differentiated cells. Here, we hypothesize that ES applied at the beginning of the differentiation process, triggers both differentiation of the hiPSC-CMs into a specialized conduction system (CS) phenotype and cell maturation. We applied ES for 15 days starting on day 0 of the differentiation process and found an increased expression of transcription factors and proteins associated with the development and function of CS including Irx3, Nkx2.5 and contactin 2, Hcn4 and Scn5a, respectively. We also found activation of intercalated disc proteins (Nrap and β-catenin). We detected ES-induced CM maturation as indicated by increased Tnni1 and Tnni3 expression. Confocal micrographs showed a shift towards expression of the gap junction protein connexin 40 in ES hiPSC-CM compared to the more dominant expression of connexin 43 in controls. Finally, analysis of functional parameters revealed that ES hiPSC-CMs exhibited faster action potential (AP) depolarization, longer intracellular Ca2+ transients, and slower AP duration at 90% of repolarization, resembling fast conducting fibers. Altogether, we provided evidence that ES during the differentiation of hiPSC to cardiomyocytes lead to development of cardiac conduction-like cells with more mature cytoarchitecture. Thus, hiPSC-CMs exposed to ES during differentiation can be instrumental to develop CS cells for cardiac disease modelling, screening individual drugs on a precison medicine type platform and support the development of novel therapeutics for arrhythmias.

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.

Protective Effects of Spermidine and Melatonin on Deltamethrin-Induced Cardiotoxicity and Neurotoxicity in Zebrafish

Increased application of the pyrethroid insecticide deltamethrin has adverse effects on the cardiac system and neurobehavior on the non-target organisms, which has raised the public's attention. Because of spermidine and melatonin considered to have cardioprotective and neuroprotective characteristics, zebrafish were utilized as the model organism to explore the protective effects of spermidine and melatonin against deltamethrin-induced toxicity. We tested the neurobehavior of zebrafish larvae through a rest/wake behavior assay, and evaluated the levels of the expression of Scn5lab, gata4, nkx2.5, hcrt, hcrtr, and aanat2 by qRT-PCR. Besides that cmlc2 was evaluated by whole-mount in situ hybridization. Results have shown that compared with control group, 0.025 mg/L deltamethrin could significantly disturb the cardiac development, downregulating the expression of Scn5lab and transcriptional factors gata4 and nkx2.5, disturbing cardiac looping, resulting in defects in cardiac morphology and function. Moreover, deltamethrin could alter the expression levels of rest/wake genes and cause hyperactivity in zebrafish larvae. Besides, compared with deltamethrin group, the exogenous 0.01 mg/L spermidine and 0.232 mg/L melatonin could significantly rescue the adverse effects of deltamethrin on the cardiac system and neurobehavior in zebrafish. This indicated that spermidine and melatonin have neuroprotective and cardioprotective effects against deltamethrin-induced adverse effects in zebrafish.

Association of Maternal Diabetes Mellitus and Polymorphisms of the NKX2.5 Gene in Children with Congenital Heart Disease: A Single Centre-Based Case-Control Study

BACKGROUND

Congenital heart disease (CHD) is one of the most common birth defects among newborns, accounting for a large proportion of infant mortality worldwide. However, the mechanisms remain largely undefinable. This study aimed to investigate the association of CHD in offspring of mothers with diabetes mellitus (DM) and single nucleotide polymorphisms (SNPs) of NKX2.5.

METHODS AND RESULTS

A case-control study of 620 mothers of CHD patients and 620 mothers of healthy children admitted to Hunan Children's Hospital from November 2017 to December 2019 was conducted. We collected the mothers' information by questionnaire and detected children's NKX2.5 variants with a MassARRAY system. The interaction coefficient (γ) was used to quantify the estimated gene-environment interactions. Univariate and multivariate analyses both showed that the infants had a higher risk of CHD if their mothers had a history of DM, including gestational DM (GDM) during this pregnancy (adjusted odds ratio [aOR = 4.98]), GDM in previous pregnancies (aOR = 4.30), and pregestational DM (PGDM) in the 3 months before this pregnancy (aOR = 6.78). Polymorphisms of the NKX2.5 gene at rs11802669 (C/C vs. T/T: aOR = 4.97; C/T vs. T/T: aOR = 2.15) and rs2277923 (T/T vs. C/C, aOR = 1.74; T/C vs. C/C, aOR = 1.61) were significantly associated with the risk of CHD in offspring. In addition, significant interactions between maternal DM and NKX2.5 genetic variants at rs11802669 (aOR = 8.12) and rs2277923 (aOR = 17.72) affecting the development of CHD were found.

CONCLUSIONS

These results suggest that maternal DM, NKX2.5 genetic variants, and their interactions are significantly associated with the risk of CHD in offspring.

Differential chamber-specific expression and regulation of long non-coding RNAs during cardiac development

Cardiovascular development is governed by a complex interplay between inducting signals such as Bmps and Fgfs leading to activation of cardiac specific transcription factors such as Nkx2.5, Mef2c and Srf that orchestrate the initial steps of cardiogenesis. Over the last decade we have witnessed the discovery of novel layers of gene regulation, i.e. post-transcriptional regulation exerted by non-coding RNAs. The function role of small non coding RNAs has been widely demonstrated, e.g. miR-1 knockout display several cardiovascular abnormalities during embryogenesis. More recently long non-coding RNAs have been also reported to modulate gene expression and function in the developing heart, as exemplified by the embryonic lethal phenotypes of Fendrr and Braveheart knock out mice, respectively. In this study, we investigated the differential expression profile during cardiogenesis of previously reported lncRNAs in heart development. Our data revealed that Braveheart, Fendrr, Carmen display a preferential adult expression while Miat, Alien, H19 preferentially display chamber-specific expression at embryonic stages. We also demonstrated that these lncRNAs are differentially regulated by Nkx2.5, Srf and Mef2c, Pitx2 > Wnt > miRNA signaling pathway and angiotensin II and thyroid hormone administration. Importantly isoform-specific expression and distinct nuclear vs cytoplasmic localization of Braveheart, Carmen and Fendrr during chamber morphogenesis is observed, suggesting distinct functional roles of these lncRNAs in atrial and ventricular chambers. Furthermore, we demonstrate by in situ hybridization a dynamic epicardial, myocardial and endocardial expression of H19 during cardiac development. Overall our data support novel roles of these lncRNAs in different temporal and tissue-restricted fashion during cardiogenesis.

CRISPR-Knockout Screen Identifies Dmap1 as a Regulator of Chemically Induced Reprogramming and Differentiation of Cardiac Progenitors

Direct in vivo reprogramming of cardiac fibroblasts into myocytes is an attractive therapeutic intervention in resolving myogenic deterioration. Current transgene-dependent approaches can restore cardiac function, but dependence on retroviral delivery and persistent retention of transgenic sequences are significant therapeutic hurdles. Chemical reprogramming has been established as a legitimate method to generate functional cell types, including those of the cardiac lineage. Here, we have extended this approach to generate progenitor cells that can differentiate into endothelial cells and cardiomyocytes using a single inhibitor protocol. Depletion of terminally differentiated cells and enrichment for proliferative cells result in a second expandable progenitor population that can robustly give rise to myofibroblasts and smooth muscle. Deployment of a genome-wide knockout screen with clustered regularly interspaced short palindromic repeats-guide RNA library to identify novel mediators that regulate the reprogramming revealed the involvement of DNA methyltransferase 1-associated protein 1 (Dmap1). Loss of Dmap1 reduced promoter methylation, increased the expression of Nkx2-5, and enhanced the retention of self-renewal, although further differentiation is inhibited because of the sustained expression of Cdh1. Our results hence establish Dmap1 as a modulator of cardiac reprogramming and myocytic induction. Stem Cells 2019;37:958-972.

Nkx-2.5 Regulates MDR1 Expression via Its Upstream Promoter in Breast Cancer Cells

BACKGROUND

Increased expression of MDR1 gene is one of the major mechanisms responsible for multidrug resistance in cancer cells. Two alternative promoters, upstream and downstream, are responsible for transcription of MDR1 gene in the human. However, the molecular mechanism regarding the transactivation of MDR1 upstream promoter (USP) has not been determined.

METHODS

Dual-luciferase reporter gene assays were used to assess the effect of Nkx-2.5 on MDR1 USP activity using reporter plasmids for human MDR1 USP and its mutants. MDR1 mRNA level was examined by quantitative real-time PCR. The direct binding of Nkx-2.5 to the USP of MDR1 was evaluated by promoter enzyme immunoassays and chromatin immunoprecipitation assays.

RESULTS

Nkx-2.5 significantly stimulates the transactivation of MDR1 USP and increases MDR1 mRNA expression in MCF7 breast cancer cells. Reporter gene assays with deleted MDR1 USPs showed that the Nkx-2.5-binding site is located between positions -71 and +12. Mutation of the Nkx-2.5-binding site at nucleotide +4 to +10 markedly reduced the Nkx-2.5-mediated activation of MDR1 USP activity. A promoter binding immunoassay and a chromatin immunoprecipitation assay revealed that Nkx-2.5 binds directly to the region +4/+10 of human MDR1 USP.

CONCLUSION

The results in the present study show Nkx-2.5 is a positive regulator for the transactivation of MDR1 USP in MCF7 breast cancer cells. Our findings will help elucidate the regulatory mechanism responsible for the multidrug resistant cancer phenotype.

Single cell RNA-seq and ATAC-seq analysis of cardiac progenitor cell transition states and lineage settlement

Formation and segregation of cell lineages forming the heart have been studied extensively but the underlying gene regulatory networks and epigenetic changes driving cell fate transitions during early cardiogenesis are still only partially understood. Here, we comprehensively characterize mouse cardiac progenitor cells (CPCs) marked by Nkx2-5 and Isl1 expression from E7.5 to E9.5 using single-cell RNA sequencing and transposase-accessible chromatin profiling (ATAC-seq). By leveraging on cell-to-cell transcriptome and chromatin accessibility heterogeneity, we identify different previously unknown cardiac subpopulations. Reconstruction of developmental trajectories reveal that multipotent Isl1+ CPC pass through an attractor state before separating into different developmental branches, whereas extended expression of Nkx2-5 commits CPC to an unidirectional cardiomyocyte fate. Furthermore, we show that CPC fate transitions are associated with distinct open chromatin states critically depending on Isl1 and Nkx2-5. Our data provide a model of transcriptional and epigenetic regulations during cardiac progenitor cell fate decisions at single-cell resolution.

Targeted next-generation sequencing of thirteen causative genes in Chinese patients with congenital hypothyroidism

To identify the spectrum and prevalence of thirteen causative genes mutations in congenital hypothyroidism (CH) patients, we collected blood samples and extracted genomic DNA of 106 CH patients, and designed a customized targeted next-generation sequencing panel containing 13 CH-causing genes to detect mutations. A total of 132 mutations were identified in 65.09% of patients (69/106) on the following nine genes: DUOX2, TG, TPO, TSHR, TTF1, TTF2, NKX2-5, PAX8 and GNAS. 69.70% (92/132) mutations related to thyroid dyshormonogenesis genes, including DUOX2 (n = 49), TG (n = 35), and TPO (n = 8). 21.21% (28/132) mutations related to thyroid dysgenesis genes, including TSHR (n = 19), TTF1 (n = 5), TTF2 (n = 1), PAX8 (n = 2), and NKX2-5 (n = 1). 9.09% (12/132) mutations related to GNAS, which was associated with thyrotropin resistance. No mutation of THRA, TSHB, IYD or SLC5A5 was detected. Among 69 mutations detected patients, 41 (59.42%) patients were two or more mutations detected, and mutations of 30 (43.48%) patients related to two or three genes. According to the pathomechanism of the mutant genes, 57.97% CH patients were classified as thyroid dyshormonogenesis. Overall, DUOX2, TG and TSHR mutations were the most common genetic defects in Chinese CH patients, and thyroid dyshormonogenesis could be the first genetic etiology of CH in Chinese. Besides, multiple mutations accounts for a part of genetic pathogenesis.

Semi-automatic detection of myocardial trabeculation using cardiovascular magnetic resonance: correlation with histology and reproducibility in a mouse model of non-compaction

BACKGROUND

The definition of left ventricular (LV) non-compaction is controversial, and discriminating between normal and excessive LV trabeculation remains challenging. Our goal was to quantify LV trabeculation on cardiovascular magnetic resonance (CMR) images in a genetic mouse model of non-compaction using a dedicated semi-automatic software package and to compare our results to the histology used as a gold standard.

METHODS

Adult mice with ventricular non-compaction were generated by conditional trabecular deletion of Nkx2-5. Thirteen mice (5 controls, 8 Nkx2-5 mutants) were included in the study. Cine CMR series were acquired in the mid LV short axis plane (resolution 0.086 × 0.086x1mm3) (11.75 T). In a sub set of 6 mice, 5 to 7 cine CMR were acquired in LV short axis to cover the whole LV with a lower resolution (0.172 × 0.172x1mm3). We used semi-automatic software to quantify the compacted mass (Mc), the trabeculated mass (Mt) and the percentage of trabeculation (Mt/Mc) on all cine acquisitions. After CMR all hearts were sliced along the short axis and stained with eosin, and histological LV contouring was performed manually, blinded from the CMR results, and Mt, Mc and Mt/Mc were quantified. Intra and interobserver reproducibility was evaluated by computing the intra class correlation coefficient (ICC).

RESULTS

Whole heart acquisition showed no statistical significant difference between trabeculation measured at the basal, midventricular and apical parts of the LV. On the mid-LV cine CMR slice, the median Mt was 0.92 mg (range 0.07-2.56 mg), Mc was 12.24 mg (9.58-17.51 mg), Mt/Mc was 6.74% (0.66-17.33%). There was a strong correlation between CMR and the histology for Mt, Mc and Mt/ Mc with respectively: r2 = 0.94 (p < 0.001), r2 = 0.91 (p < 0.001), r2 = 0.83 (p < 0.001). Intra- and interobserver reproducibility was 0.97 and 0.8 for Mt; 0.98 and 0.97 for Mc; 0.96 and 0.72 for Mt/Mc, respectively and significantly more trabeculation was observed in the Mc Mutant mice than the controls.

CONCLUSION

The proposed semi-automatic quantification software is accurate in comparison to the histology and reproducible in evaluating Mc, Mt and Mt/ Mc on cine CMR.

Heterogeneity of adult masseter muscle satellite cells with cardiomyocyte differentiation potential

Although resident cardiac stem cells have been reported, regeneration of functional cardiomyocytes (CMs) remains a challenge. The present study identifies an alternative progenitor source for CM regeneration without the need for genetic manipulation or invasive heart biopsy procedures. Unlike limb skeletal muscles, masseter muscles (MM) in the mouse head are developed from Nkx2-5 mesodermal progenitors. Adult masseter muscle satellite cells (MMSCs) display heterogeneity in developmental origin and cell phenotypes. The heterogeneous MMSCs that can be characterized by cell sorting based on stem cell antigen-1 (Sca1) show different lineage potential. While cardiogenic potential is preserved in Sca1+ MMSCs as shown by expression of cardiac progenitor genes (including Nkx2-5), skeletal myogenic capacity is maintained in Sca1- MMSCs with Pax7 expression. Sca1+ MMSC-derived beating cells express cardiac genes and exhibit CM-like morphology. Electrophysiological properties of MMSC-derived CMs are demonstrated by calcium transients and action potentials. These findings show that MMSCs could serve as a novel cell source for cardiomyocyte replacement.

Fenbuconazole exposure impacts the development of zebrafish embryos

Fenbuconazole (FBZ), a triazole-containing fungicide, is widely used in agriculture and horticulture. In the present study, the development and cardiac functioning were observed and determined in zebrafish embryos exposed to FBZ at 5, 50 and 500 ng/L nominal concentrations for 72 h. The results showed that 500 ng/L FBZ significantly increased pericardial edema rate, spine curvature rate, disturbed cardiac function, and led a shortened lower jaw. The transcription of genes such as tbx5, nkx2.5, tnnt2, gata4, bmp2b, myl7 was altered, which might be responsible for the cardiac developmental and functioning defects in the larvae. The deformation in bone development might be related with the impaired transcription levels of shh and bmp2b. The transcription of cyp26a1 (encoding retinoic acid metabolism enzyme) was significantly up-regulated in the 500 ng/L group, which might be a reason causing the teratogenic effect of FBZ. These results suggest that FBZ could have toxic effects on embryonic development, which should be considered in the risk evaluation of FBZ application.

ETV1 activates a rapid conduction transcriptional program in rodent and human cardiomyocytes

Rapid impulse propagation is a defining attribute of the pectinated atrial myocardium and His-Purkinje system (HPS) that safeguards against atrial and ventricular arrhythmias, conduction block, and myocardial dyssynchrony. The complex transcriptional circuitry that dictates rapid conduction remains incompletely understood. Here, we demonstrate that ETV1 (ER81)-dependent gene networks dictate the unique electrophysiological characteristics of atrial and His-Purkinje myocytes. Cardiomyocyte-specific deletion of ETV1 results in cardiac conduction abnormalities, decreased expression of rapid conduction genes (Nkx2-5, Gja5, and Scn5a), HPS hypoplasia, and ventricularization of the unique sodium channel properties that define Purkinje and atrial myocytes in the adult heart. Forced expression of ETV1 in postnatal ventricular myocytes (VMs) reveals that ETV1 promotes a HPS gene signature while diminishing ventricular and nodal gene networks. Remarkably, ETV1 induction in human induced pluripotent stem cell-derived cardiomyocytes increases rapid conduction gene expression and inward sodium currents, converting them towards a HPS phenotype. Our data identify a cardiomyocyte-autonomous, ETV1-dependent pathway that is responsible for specification of rapid conduction zones in the heart and demonstrate that ETV1 is sufficient to promote a HPS transcriptional and functional program upon VMs.

Ronin Governs Early Heart Development by Controlling Core Gene Expression Programs

Ronin (THAP11), a DNA-binding protein that evolved from a primordial DNA transposon by molecular domestication, recognizes a hyperconserved promoter sequence to control developmentally and metabolically essential genes in pluripotent stem cells. However, it remains unclear whether Ronin or related THAP proteins perform similar functions in development. Here, we present evidence that Ronin functions within the nascent heart as it arises from the mesoderm and forms a four-chambered organ. We show that Ronin is vital for cardiogenesis during midgestation by controlling a set of critical genes. The activity of Ronin coincided with the recruitment of its cofactor, Hcf-1, and the elevation of H3K4me₃ levels at specific target genes, suggesting the involvement of an epigenetic mechanism. On the strength of these findings, we propose that Ronin activity during cardiogenesis offers a template to understand how important gene programs are sustained across different cell types within a developing organ such as the heart.

•

Ronin displays complex expression patterns during embryogenesis

•

Ronin is critical for heart growth

•

Ronin regulates genetic growth programs

•

Ronin binding influences H3K4me₃ levels at target genes

Novel Point Mutations in the NKX2.5 Gene in Pediatric Patients with Non-Familial Congenital Heart Disease

Background and objective: Congenital heart disease (CHD) is the most common birth abnormality in the structure or function of the heart that affects approximately 1% of all newborns. Despite its prevalence and clinical importance, the etiology of CHD remains mainly unknown. Somatic and germline mutations in cardiac specific transcription factor genes have been identified as the factors responsible for various forms of CHD, particularly ventricular septal defects (VSDs), tetralogy of Fallot (TOF), and atrial septal defects (ASDs). p. NKX2.5 is a homeodomain protein that controls many of the physiological processes in cardiac development including specification and proliferation of cardiac precursors. The aim of our study was to evaluate the NKX2.5 gene mutations in sporadic pediatric patients with clinical diagnosis of congenital heart malformations. Materials and methods: In this study, we investigated mutations of the NKX2.5 gene’s coding region in 105 Iranian pediatric patients with non-familial CHD by polymerase chain reaction-single stranded conformation polymorphism (PCR-SSCP) and direct sequencing. Results: We observed a total of four mutations, of which, two were novel DNA sequence variants in the coding region of exon 1 (c. 95 A > T and c. 93 A > T) and two others were previously reported as single-nucleotide polymorphisms (SNPs), namely rs72554028 (c. 2357 G > A) and rs3729753 (c. 606 G > C) in exon 2. Further, observed mutations are completely absent in normal healthy individuals (n = 92). Conclusion: These results suggest that NKX2.5 mutations are highly rare in CHD patients. However, in silico analysis proves that c.95 A > T missense mutation in NKX2.5 gene is probably pathogenic and may be contributing to the risk of sporadic CHD in the Iranian population.

[Association of single nucleotide polymorphisms of transcription factors with congenital heart diseases in the Chinese population: a Meta analysis]

OBJECTIVE

To study the association of single nucleotide polymorphisms (SNPs) of transcription factors (NKX2.5, GATA4, TBX5, and FOG2) with congenital heart disease (CHD) in the Chinese population.

METHODS

PubMed, Google Scholar, CNKI, Wanfang Data, and Weipu Data were searched for articles on the association of SNPs of target genes with CHD in the Chinese population. If one locus was mentioned in at least two articles, the random or fixed effect model was used to perform a pooled analysis of study results and to calculate the pooled OR and its 95%CI. If a locus was mentioned in only one article, related data were extracted from this article to analyze the association between the SNPs of this locus and CHD.

RESULTS

Twenty-three articles were included. The Meta analysis showed that there were significant differences between the CHD and control groups in the genotype and allele frequencies of GATA4 rs1139244 and rs867858 and the genotype frequency of GATA4 rs904018, while there were no significant differences in the SNPs of the other genetic loci between the two groups. The single-article analysis showed that there were significant differences between the two groups in the allele frequencies of NKX2.5 rs118026695/rs703752, GATA4 rs884662/rs12825/rs12458/rs3203358/rs4841588, and TBX5 rs6489956. There were no significant differences in the SNPs of FOG2 locus between the two groups.

CONCLUSIONS

The SNPs of some loci in NKX2.5, GATA4, and TBX5 are associated with CHD in the Chinese population, but the association between the SNPs of FOG2 locus and the development of CHD has not been found yet.

[Association of single nucleotide polymorphisms of transcription factors with congenital heart diseases in the Chinese population: a Meta analysis]

OBJECTIVE

To study the association of single nucleotide polymorphisms (SNPs) of transcription factors (NKX2.5, GATA4, TBX5, and FOG2) with congenital heart disease (CHD) in the Chinese population.

METHODS

PubMed, Google Scholar, CNKI, Wanfang Data, and Weipu Data were searched for articles on the association of SNPs of target genes with CHD in the Chinese population. If one locus was mentioned in at least two articles, the random or fixed effect model was used to perform a pooled analysis of study results and to calculate the pooled OR and its 95%CI. If a locus was mentioned in only one article, related data were extracted from this article to analyze the association between the SNPs of this locus and CHD.

RESULTS

Twenty-three articles were included. The Meta analysis showed that there were significant differences between the CHD and control groups in the genotype and allele frequencies of GATA4 rs1139244 and rs867858 and the genotype frequency of GATA4 rs904018, while there were no significant differences in the SNPs of the other genetic loci between the two groups. The single-article analysis showed that there were significant differences between the two groups in the allele frequencies of NKX2.5 rs118026695/rs703752, GATA4 rs884662/rs12825/rs12458/rs3203358/rs4841588, and TBX5 rs6489956. There were no significant differences in the SNPs of FOG2 locus between the two groups.

CONCLUSIONS

The SNPs of some loci in NKX2.5, GATA4, and TBX5 are associated with CHD in the Chinese population, but the association between the SNPs of FOG2 locus and the development of CHD has not been found yet.

Isolation and characterization of ventricular-like cells derived from NKX2-5eGFP/w and MLC2vmCherry/w double knock-in human pluripotent stem cells

Human pluripotent stem cell (hPSC)-derived cardiomyocytes (CMs) are a promising source for cell transplantation into the damaged heart, which has limited regenerative ability. Many methods have been developed to obtain large amounts of functional CMs from hPSCs for therapeutic applications. However, during the differentiation process, a mixed population of various cardiac cells, including ventricular, atrial, and pacemaker cells, is generated, which hampers the proper functional analysis and evaluation of cell properties. Here, we established NKX2-5^eGFP/w and MLC2v^mCherry/w hPSC double knock-ins that allow for labeling, tracing, purification, and analysis of the development of ventricular cells from early to late stages. As with the endogenous transcriptional activities of these genes, MLC2v-mCherry expression following NKX2-5-eGFP expression was observed under previously established culture conditions, which mimic the in vivo cardiac developmental process. Patch-clamp and microelectrode array electrophysiological analyses showed that the NKX2-5 and MLC2v double-positive cells possess ventricular-like properties. The results demonstrate that the NKX2-5^eGFP/w and MLC2v^mCherry/w hPSCs provide a powerful model system to capture region-specific cardiac differentiation from early to late stages. Our study would facilitate subtype-specific cardiac development and functional analysis using the hPSC-derived sources.

In vivo cardiovascular toxicity induced by acetochlor in zebrafish larvae

The risk of acetochlor to human health is still unclear, prompting concern over its risk, especially to pesticide suicides population, occupational population (farmers, retailers and pharmaceutical workers), and special population (young children and infants, pregnant women, older people, and those with compromised immune systems). This study was to explore the toxic effect and the possible mechanism of toxic action of acetochlor using zebrafish larvae whose toxicity profiles have been confirmed to be strikingly similar with mammalian. The result indicated that the toxic target organ of acetochlor was cardiovascular system. Thus, cardiovascular toxicity evaluation was investigated systematically. The main phenotypes of cardiovascular toxicity induced by acetochlor were bradycardia, pericardial edema, circulation defect, and thrombosis; Malformed heart was confirmed by histopathological examination. Thrombosis which maybe triggered by bradycardia was further studied using o-dianisidine for erythrocyte staining; Substantial thrombus in the caudal vein and significantly reduced heart red blood cells (RBCs) intensity which can reflect the thrombosis degree were observed in zebrafish in a concentration-dependent manner. Additionally, the mRNA expression level of Nkx2.5 and Gata4 related to induction of cardiac program were down-regulated significantly by quantitative real-time polymerase chain reaction (qRT-PCR), which could cause defects in the cardiovascular system. For the first time, our results demonstrated that acetochlor induced cardiovascular toxicity, and down-regulation of Nkx2.5 and Gata4 might be its possible molecular basis. Our data generated here might provide novel insights into cardiovascular disease risk following acetochlor exposure to human, especially to pesticide suicides population, occupational population and special population.

Morpho-functional characterization of the systemic venous pole of the reptile heart

Mammals evolved from reptile-like ancestors, and while the mammalian heart is driven by a distinct sinus node, a sinus node is not apparent in reptiles. We characterized the myocardial systemic venous pole, the sinus venosus, in reptiles to identify the dominant pacemaker and to assess whether the sinus venosus remodels and adopts an atrium-like phenotype as observed in mammals. Anolis lizards had an extensive sinus venosus of myocardium expressing Tbx18. A small sub-population of cells encircling the sinuatrial junction expressed Isl1, Bmp2, Tbx3, and Hcn4, homologues of genes marking the mammalian sinus node. Electrical mapping showed that hearts of Anolis lizards and Python snakes were driven from the sinuatrial junction. The electrical impulse was delayed between the sinus venosus and the right atrium, allowing the sinus venosus to contract and aid right atrial filling. In proximity of the systemic veins, the Anolis sinus venosus expressed markers of the atrial phenotype Nkx2-5 and Gja5. In conclusion, the reptile heart is driven by a pacemaker region with an expression signature similar to that of the immature sinus node of mammals. Unlike mammals, reptiles maintain a sinuatrial delay of the impulse, allowing the partly atrialized sinus venosus to function as a chamber.

[NKX2.5 and TBX5 gene mutations in in vitro fertilization children with congenital heart disease]

OBJECTIVE

To explore the differences of NKX2.5 and TBX5 gene mutations between in vitro fertilization (IVF) children with congenital heart disease (CHD) and naturally conceived children with CHD.

METHODS

Blood samples from 68 IVF children with CHD and 98 naturally conceived children with CHD were collected. The mutations in coding regions 1 and 2 of the NKX2.5 gene, and coding regions 4, 5, and 8 of the TBX5 gene were examined by polymerase chain reaction (PCR) and DNA sequencing.

RESULTS

An A-to-G mutation at nucleotide 63 (c.63A>G) in coding region 1 of the NKX2.5 gene was found in both IVF and naturally conceived children with CHD. There were no significant differences in genotype and allele frequencies at c.63A>G locus of the NKX2.5 gene between the two groups. No mutations were detected in coding region 2 of the NKX2.5 gene and coding regions 4, 5 and 8 of the TBX5 gene.

CONCLUSIONS

There is no difference in NKX2.5 and TBX5 gene mutations between IVF and naturally conceived children with CHD. Therefore, it is presumed that assisted reproductive technology may not lead to mutations in the NKX2.5 and TBX5 genes.

Developmental toxicity in flounder embryos exposed to crude oils derived from different geographical regions

Crude oils from distinct geographical regions have distinct chemical compositions, and, as a result, their toxicity may be different. However, developmental toxicity of crude oils derived from different geographical regions has not been extensively characterized. In this study, flounder embryos were separately exposed to effluents contaminated by three crude oils including: Basrah Light (BLO), Pyrenees (PCO), and Sakhalin Vityaz (SVO), in addition to a processed fuel oil (MFO-380), to measure developmental toxicity and for gene expressions. Each oil possessed a distinct chemical composition. Edema defect was highest in embryos exposed to PCO and MFO-380 that both have a greater fraction of three-ring PAHs (33% and 22%, respectively) compared to BLO and SVO. Observed caudal fin defects were higher in embryos exposed to SVO and MFO-380, which are both dominated by naphthalenes (81% and 52%, respectively). CYP1A gene expressions were also highest in embryos exposed to SVO and MFO-380. Higher incidence of cardiotoxicity and lower nkx 2.5 expression were detected in embryos exposed to PCO. Unique gene expression profiles were observed in embryos exposed to crude oils with distinct compositions. This study demonstrates that crude oils of different geographical origins with different compositional characteristics induce developmental toxicity to different degrees.

Next-generation sequencing of NKX2.1, FOXE1, PAX8, NKX2.5, and TSHR in 100 Chinese patients with congenital hypothyroidism and athyreosis

BACKGROUND

The abnormal expression of certain transcription factors (NKX2.1, FOXE1, NKX2.5, and PAX8) and thyroid stimulating hormone receptor (TSHR) genes has been associated with athyreosis, which is a form of thyroid dysgenesis (TD). We aimed to identify candidate gene mutations in CH patients with athyreosis and to establish the genotype-phenotype correlations in a Chinese population.

METHODS

The exons and flanking sequences of NKX2.1, FOXE1, NKX2.5, PAX8, and TSHR were screened by next-generation sequencing and further confirmed by direct Sanger sequencing. The mutation frequencies were calculated and compared against databases. The relationship between genotype and phenotype was also determined.

RESULTS

Seven variants were detected in TSHR-p.P52T, p.G132R, p.M164K, p.R450H, p.C700E, p.A522V, and p.R528S. The p. G132R, p. M164K and p. R528S variants were first identified in public databases. Five variants (p.G44D, p.G360V, p.R401Q, p.L418I, and p.E453Q) were found in NKX2.1 and one variant (p.P243T) was detected in FOXE1. In addition, one variant (p.N291I) was found in NKX2.5 and two variants (p.A355V and c.-26G>A) were detected in PAX8.

CONCLUSIONS

Our study indicated that TSHR mutations have phenotypic variability and has further expanded the mutation spectrum of TSHR. We also revealed that the rate of NKX2.1, FOXE1, NKX2.5, and PAX8 mutations were low in patients with CH and athyreosis, in contrast to the higher rate of TSHR mutations.

Disruption of spatiotemporal hypoxic signaling causes congenital heart disease in mice

Congenital heart disease (CHD) represents the most prevalent inborn anomaly. Only a minority of CHD cases are attributed to genetic causes, suggesting a major role of environmental factors. Nonphysiological hypoxia during early pregnancy induces CHD, but the underlying reasons are unknown. Here, we have demonstrated that cells in the mouse heart tube are hypoxic, while cardiac progenitor cells (CPCs) expressing islet 1 (ISL1) in the secondary heart field (SHF) are normoxic. In ISL1+ CPCs, induction of hypoxic responses caused CHD by repressing Isl1 and activating NK2 homeobox 5 (Nkx2.5), resulting in decreased cell proliferation and enhanced cardiomyocyte specification. We found that HIF1α formed a complex with the Notch effector hes family bHLH transcription factor 1 (HES1) and the protein deacetylase sirtuin 1 (SIRT1) at the Isl1 gene. This complex repressed Isl1 in the hypoxic heart tube or following induction of ectopic hypoxic responses. Subsequently, reduced Isl1 expression abrogated ISL1-dependent recruitment of histone deacetylases HDAC1/5, inhibiting Nkx2.5 expression. Inactivation of Sirt1 in ISL1+ CPCs blocked Isl1 suppression via the HIF1α/HES1/SIRT1 complex and prevented CHDs induced by pathological hypoxia. Our results indicate that spatial differences in oxygenation of the developing heart serve as signals to control CPC expansion and cardiac morphogenesis. We propose that physiological hypoxia coordinates homeostasis of CPCs, providing mechanistic explanations for some nongenetic causes of CHD.

Transcription factor-induced activation of cardiac gene expression in human c-kit+ cardiac progenitor cells

Although transplantation of c-kit+ cardiac progenitor cells (CPCs) significantly alleviates post-myocardial infarction left ventricular dysfunction, generation of cardiomyocytes by exogenous CPCs in the recipient heart has often been limited. Inducing robust differentiation would be necessary for improving the efficacy of the regenerative cardiac cell therapy. We assessed the hypothesis that differentiation of human c-kit+ CPCs can be enhanced by priming them with cardiac transcription factors (TFs). We introduced five different TFs (Gata4, MEF2C, NKX2.5, TBX5, and BAF60C) into CPCs, either alone or in combination, and then examined the expression of marker genes associated with the major cardiac cell types using quantitative RT-PCR. When introduced individually, Gata4 and TBX5 induced a subset of myocyte markers. Moreover, Gata4 alone significantly induced smooth muscle cell and fibroblast markers. Interestingly, these gene expression changes brought by Gata4 were also accompanied by morphological changes. In contrast, MEF2C and NKX2.5 were largely ineffective in initiating cardiac gene expression in CPCs. Surprisingly, introduction of multiple TFs in different combinations mostly failed to act synergistically. Likewise, addition of BAF60C to Gata4 and/or TBX5 did not further potentiate their effects on cardiac gene expression. Based on our results, it appears that GATA4 is able to potentiate gene expression programs associated with multiple cardiovascular lineages in CPCs, suggesting that GATA4 may be effective in priming CPCs for enhanced differentiation in the setting of stem cell therapy.

The acceleration of cardiomyogenesis in embryonic stem cells in vitro by serum depletion does not increase the number of developed cardiomyocytes

The differentiation of pluripotent embryonic stem (ES) cells into various lineages in vitro represents an important tool for studying the mechanisms underlying mammalian embryogenesis. It is a key technique in studies evaluating the molecular mechanisms of cardiomyogenesis and heart development and also in embryotoxicology. Herein, modest modifications of the basic protocol for ES cell differentiation into cardiomyocytes were evaluated in order to increase the yield and differentiation status of developed cardiomyocytes. Primarily, the data show that ES cell cultivation in the form of non-adherent embryoid bodies (EBs) for 5 days compared to 8 days significantly improved cardiomyogenic differentiation. This is illustrated by the appearance of beating foci in the adherent EBs layer at earlier phases of differentiation from day 10 up to day 16 and by the significantly higher expression of genes characteristic of cardiomyogenic differentiation (sarcomeric alpha actinin, myosin heavy chain alpha and beta, myosin light chain 2 and 7, and transcriptional factor Nkx2.5) in EBs cultivated under non-adherent conditions for 5 days. The ratio of cardiomyocytes per other cells was also potentiated in EBs cultivated in non-adherent conditions for only 5 days followed by cultivation in adherent serum-free culture conditions. Nevertheless, the alteration in the percentage of beating foci among these two tested cultivation conditions vanished at later phases and also did not affect the total number of cardiomyocytes determined as myosin heavy chain positive cells at the end of the differentiation process on day 20. Thus, although these modifications of the conditions of ES cells differentiation may intensify cardiomyocyte differentiation, the final count of cardiomyocytes might not change. Thus, serum depletion was identified as a key factor that intensified cardiomyogenesis. Further, the treatment of EBs with N-acetylcysteine, a reactive oxygen species scavenger, did not affect the observed increase in cardiomyogenesis under serum depleted conditions. Interestingly, a mild induction of the ventricular-like phenotype of cardiomyocytes was observed in 5-day-old EBs compared to 8-day-old EBs. Overall, these findings bring crucial information on the mechanisms of ES cells differentiation into cardiomyocytes and on the establishment of efficient protocols for the cardiomyogenic differentiation of ES cells. Further, the importance of determining the absolute number of formed cardiomyocyte-like cells per seeded pluripotent cells in contrast to the simple quantification of the ratios of cells is highlighted.

Maternal exposure to ambient PM2.5 exaggerates fetal cardiovascular maldevelopment induced by homocysteine in rats

Maternal exposure to airborne particulate matter with aerodynamic diameter <2.5 µm (PM_2.5 ) during pregnancy and lactation periods is associated with filial congenital cardiovascular diseases. This study aimed to investigate the toxic effects of maternal exposure to ambient levels of PM_2.5 on filial cardiovascular maldevelopment induced by homocysteine. Using a 2 × 2 factorial design, rats were randomized into four groups and were exposed to ambient PM_2.5 or filtered air (FA) throughout the pregnancy and lactation periods coupled with the administration of either homocysteine (HCY) or normal saline (NS) daily from gestation days 8-10. Morphological changes in the heart, myocardial apoptosis, expressions of cardiac progenitor transcriptional factors, and levels of cytokines were investigated in the offspring. The apoptosis-like changes in the myocardium were seen in the FA plus HCY-treated group and more obviously in the PM_2.5 plus HCY-treated group, which was in accordance with an increased myocardial apoptosis rate in the two groups. PM_2.5 exposure resulted in significantly decreased Nkx2-5 protein level and GATA4 and Nkx2-5 mRNA expressions, and significantly increased TNF-α and IL-1β levels. There were significant interactions between PM_2.5 exposure and HCY-treatment that PM_2.5 exposure reduced Nkx2-5 protein levels and GATA4 and Nkx2-5 mRNA expressions in the HCY-treated groups. These results suggest that maternal exposure to PM_2.5 , even at the ambient levels in urban regions in China, exaggerates filial cardiovascular maldevelopment induced by HCY in a murine model, exacerbating structural abnormalities in the filial cardiac tissue, which is possibly associated with oxidative stress and reduced GATA4 and Nkx2-5 transcription factor expressions. © 2016 Wiley Periodicals, Inc. Environ Toxicol 32: 877-889, 2017.

Excess Imidacloprid Exposure Causes the Heart Tube Malformation of Chick Embryos

As a neonicotinoid pesticide, imidacloprid is widely used to control sucking insects on agricultural planting and fleas on domestic animals. However, the extent to which imidacloprid exposure has an influence on cardiogensis in early embryogenesis is still poorly understood. In vertebrates, the heart is the first organ to be formed. In this study, to address whether imidacloprid exposure affects early heart development, the early chick embryo has been used as an experimental model because of its accessibility at its early developmental stage. The results demonstrate that exposure of the early chick embryo to imidacloprid caused malformation of heart tube. Furthermore, the data reveal that down-regulation of GATA4, NKX2.5, and BMP4 and up-regulation of Wnt3a led to aberrant cardiomyocyte differentiation. In addition, imidacloprid exposure interfered with basement membrane breakdown, E-cadherin/laminin expression, and mesoderm formation during the epithelial-mesenchymal transition (EMT) in gastrula chick embryos. Finally, the DiI-labeled cell migration trajectory indicated that imidacloprid restricted the cell migration of cardiac progenitors to primary heart field in gastrula chick embryos. A similar observation was also obtained from the cell migration assay of scratch wounds in vitro. Additionally, imidacloprid exposure negatively affected the cytoskeleton structure and expression of corresponding adhesion molecules. Taken together, these results reveal that the improper EMT, cardiac progenitor migration, and differentiation are responsible for imidacloprid exposure-induced malformation of heart tube during chick embryo development.

Disruption of cardiogenesis in human embryonic stem cells exposed to trichloroethylene

Trichloroethylene (TCE) is ubiquitous in our living environment, and prenatal exposure to TCE is reported to cause congenital heart disease in humans. Although multiple studies have been performed using animal models, they have limited value in predicting effects on humans due to the unknown species-specific toxicological effects. To test whether exposure to low doses of TCE induces developmental toxicity in humans, we investigated the effect of TCE on human embryonic stem cells (hESCs) and cardiomyocytes (derived from the hESCs). In the current study, hESCs cardiac differentiation was achieved by using differentiation medium consisting of StemPro-34. We examined the effects of TCE on cell viability by cell growth assay and cardiac inhibition by analysis of spontaneously beating cluster. The expression levels of genes associated with cardiac differentiation and Ca²⁺ channel pathways were measured by immunofluorescence and qPCR. The overall data indicated the following: (1) significant cardiac inhibition, which was characterized by decreased beating clusters and beating rates, following treatment with low doses of TCE; (2) significant up-regulation of the Nkx2.5/Hand1 gene in cardiac progenitors and down regulation of the Mhc-7/cTnT gene in cardiac cells; and (3) significant interference with Ca²⁺ channel pathways in cardiomyocytes, which contributes to the adverse effect of TCE on cardiac differentiation during early embryo development. Our results confirmed the involvement of Ca²⁺ turnover network in TCE cardiotoxicity as reported in animal models, while the inhibition effect of TCE on the transition of cardiac progenitors to cardiomyocytes is unique to hESCs, indicating a species-specific effect of TCE on heart development. This study provides new insight into TCE biology in humans, which may help explain the development of congenital heart defects after TCE exposure. © 2015 Wiley Periodicals, Inc. Environ Toxicol 31: 1372-1380, 2016.

Transcription factor ETV1 is essential for rapid conduction in the heart

Rapid impulse propagation in the heart is a defining property of pectinated atrial myocardium (PAM) and the ventricular conduction system (VCS) and is essential for maintaining normal cardiac rhythm and optimal cardiac output. Conduction defects in these tissues produce a disproportionate burden of arrhythmic disease and are major predictors of mortality in heart failure patients. Despite the clinical importance, little is known about the gene regulatory network that dictates the fast conduction phenotype. Here, we have used signal transduction and transcriptional profiling screens to identify a genetic pathway that converges on the NRG1-responsive transcription factor ETV1 as a critical regulator of fast conduction physiology for PAM and VCS cardiomyocytes. Etv1 was highly expressed in murine PAM and VCS cardiomyocytes, where it regulates expression of Nkx2-5, Gja5, and Scn5a, key cardiac genes required for rapid conduction. Mice deficient in Etv1 exhibited marked cardiac conduction defects coupled with developmental abnormalities of the VCS. Loss of Etv1 resulted in a complete disruption of the normal sodium current heterogeneity that exists between atrial, VCS, and ventricular myocytes. Lastly, a phenome-wide association study identified a link between ETV1 and bundle branch block and heart block in humans. Together, these results identify ETV1 as a critical factor in determining fast conduction physiology in the heart.

A Novel Missense Mutation of GATA4 in a Chinese Family with Congenital Heart Disease

Background

Congenital heart disease (CHD) is the most prevalent type of birth defect in human, with high morbidity in infant. Several genes essential for heart development have been identified. GATA4 is a pivotal transcription factor that can regulate the cardiac development. Many GATA4 mutations have been identified in patients with different types of CHD.

Aims

In this study, the NKX2-5, HAND1 and GATA4 coding regions were sequenced in a family spanning three generations in which seven patients had CHD. Disease-causing potential variation in this family was evaluated by bioinformatics programs and the transcriptional activity of mutant protein was analyzed by the dual luciferase reporter assay.

Results

A novel GATA4 mutation, c.C931T (p.R311W), was identified and co-segregated with the affected patients in this family. The bioinformatics programs predicted this heterozygous mutation to be deleterious and the cross-species alignment of GATA4 sequences showed that the mutation occurred within a highly conserved amino acid. Even though it resided in the nuclear localization signal domain, the mutant protein didn’t alter its intracellular distribution. Nevertheless, further luciferase reporter assay demonstrated that the p.R311W mutation reduced the ability of GATA4 to activate its downstream target gene.

Conclusions

Our study identified a novel mutation in GATA4 that likely contributed to the CHD in this family. This finding expanded the spectrum of GATA4 mutations and underscored the pathogenic correlation between GATA4 mutations and CHD.

Inhibition of Cardiomyogenesis in Embryocarcinoma Cells Induced by Long-Term High Level of Glucose

BACKGROUND/AIMS

Cardiac myocytes constitute the first differentiated cell type during mammalian heart formation with the ability to beat spontaneously and rhythmically. Hyperglycemia is a primary risk factor for cardiovascular disease in pre-gestational diabetes mellitus (PGDM). However, the impact that hyperglycemia has on cardiac progenitors or on precursors differentiation remains poorly understood. The aim of the present study is to investigate whether hyperglycemia affects cardiomyogenesis of embryocarcinoma cells.

METHODS

P19CL6 cells differentiation induced by 1% DMSO was evaluated under either normal glucose (5.6 mmol/L) or high level of glucose concentrations (20 mmol/L or 40 mmol/L). To investigate the effect of long-term high level of glucose on cardiomyocytes differentiation, sarcomeric α-actinin, peroxisome proliferator-activated receptor coactivator-1 (PGC-1α), transcription factor GATA4 and Nkx2.5 were assessed by qRT-PCR analysis, western blot and immunofluorescence.

RESULTS

We observed that long-term high level of glucose markedly reduced P19CL6 cells differentiation into cardiomyocytes. The change in PGC-1α expression was consistent with changes in cardiac muscle myosin expression after exposure to 20 mmol/L or 40 mmol/L of glucose. On the other hand, the high level of glucose concentration profoundly decreased both GATA4 and Nkx2-5 expressions from day 6 to day 12 after differentiation, which was induced by 1% DMSO.

CONCLUSION

Our results elucidate that the effect resulting from the long-term exposure of cardiac progenitors to high level of glucose is associated with decreased expression of GATA4 and Nkx2.5, providing a novel mechanism by which high glucose is able to affect cell differentiation.

Establishment of Relational Model of Congenital Heart Disease Markers and GO Functional Analysis of the Association between Its Serum Markers and Susceptibility Genes

PURPOSE

The purpose of present study was to construct the best screening model of congenital heart disease serum markers and to provide reference for further prevention and treatment of the disease.

METHODS

Documents from 2006 to 2014 were collected and meta-analysis was used for screening susceptibility genes and serum markers closely related to the diagnosis of congenital heart disease. Data of serum markers were extracted from 80 congenital heart disease patients and 80 healthy controls, respectively, and then logistic regression analysis and support vector machine were utilized to establish prediction models of serum markers and Gene Ontology (GO) functional annotation.

RESULTS

Results showed that NKX2.5, GATA4, and FOG2 were susceptibility genes of congenital heart disease. CRP, BNP, and cTnI were risk factors of congenital heart disease (p < 0.05); cTnI, hs-CRP, BNP, and Lp(a) were significantly close to congenital heart disease (p < 0.01). ROC curve indicated that the accuracy rate of Lp(a) and cTnI, Lp(a) and BNP, and BNP and cTnI joint prediction was 93.4%, 87.1%, and 97.2%, respectively. But the detection accuracy rate of the markers' relational model established by support vector machine was only 85%. GO analysis suggested that NKX2.5, GATA4, and FOG2 were functionally related to Lp(a) and BNP.

CONCLUSIONS

The combined markers model of BNP and cTnI had the highest accuracy rate, providing a theoretical basis for the diagnosis of congenital heart disease.

Irx3 is required for postnatal maturation of the mouse ventricular conduction system

The ventricular conduction system (VCS) orchestrates the harmonious contraction in every heartbeat. Defects in the VCS are often associated with life-threatening arrhythmias and also promote adverse remodeling in heart disease. We have previously established that the Irx3 homeobox gene regulates rapid electrical propagation in the VCS by modulating the transcription of gap junction proteins Cx40 and Cx43. However, it is unknown whether other factors contribute to the conduction defects observed in Irx3 knockout (Irx3(-/-)) mice. In this study, we show that during the early postnatal period, Irx3(-/-) mice develop morphological defects in the VCS which are temporally dissociated from changes in gap junction expression. These morphological defects were accompanied with progressive changes in the cardiac electrocardiogram including right bundle branch block. Hypoplastic VCS was not associated with increased apoptosis of VCS cardiomyocytes but with a lack of recruitment and maturation of ventricular cardiomyocytes into the VCS. Computational analysis followed by functional verification revealed that Irx3 promotes VCS-enriched transcripts targeted by Nkx2.5 and/or Tbx5. Altogether, these results indicate that, in addition to ensuring the appropriate expression of gap junctional channels in the VCS, Irx3 is necessary for the postnatal maturation of the VCS, possibly via its interactions with Tbx5 and Nkx2.5.