The transcription factor EBF1 non-cell-autonomously regulates cardiac growth and differentiation.

Reciprocal interactions between non-myocytes and cardiomyocytes regulate cardiac growth and differentiation. Here, we report that the transcription factor Ebf1 is highly expressed in non-myocytes and potently regulates heart development. Ebf1-deficient hearts display myocardial hypercellularity and reduced cardiomyocyte size, ventricular conduction system hypoplasia, and conduction system disease. Growth abnormalities in Ebf1 knockout hearts are observed as early as embryonic day 13.5. Transcriptional profiling of Ebf1-deficient embryonic cardiac non-myocytes demonstrates dysregulation of Polycomb repressive complex 2 targets, and ATAC-Seq reveals altered chromatin accessibility near many of these same genes. Gene set enrichment analysis of differentially expressed genes in cardiomyocytes isolated from E13.5 hearts of wild-type and mutant mice reveals significant enrichment of MYC targets and, consistent with this finding, we observe increased abundance of MYC in mutant hearts. EBF1-deficient non-myocytes, but not wild-type non-myocytes, are sufficient to induce excessive accumulation of MYC in co-cultured wild-type cardiomyocytes. Finally, we demonstrate that BMP signaling induces Ebf1 expression in embryonic heart cultures and controls a gene program enriched in EBF1 targets. These data reveal a previously unreported non-cell-autonomous pathway controlling cardiac growth and differentiation.

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.

Activin A directly impairs human cardiomyocyte contractile function indicating a potential role in heart failure development.

Activin A has been linked to cardiac dysfunction in aging and disease, with elevated circulating levels found in patients with hypertension, atherosclerosis, and heart failure. Here, we investigated whether Activin A directly impairs cardiomyocyte (CM) contractile function and kinetics utilizing cell, tissue, and animal models. Hydrodynamic gene delivery-mediated overexpression of Activin A in wild-type mice was sufficient to impair cardiac function, and resulted in increased cardiac stress markers (N-terminal pro-atrial natriuretic peptide) and cardiac atrophy. In human-induced pluripotent stem cell-derived (hiPSC) CMs, Activin A caused increased phosphorylation of SMAD2/3 and significantly upregulated SERPINE1 and FSTL3 (markers of SMAD2/3 activation and activin signaling, respectively). Activin A signaling in hiPSC-CMs resulted in impaired contractility, prolonged relaxation kinetics, and spontaneous beating in a dose-dependent manner. To identify the cardiac cellular source of Activin A, inflammatory cytokines were applied to human cardiac fibroblasts. Interleukin -1ß induced a strong upregulation of Activin A. Mechanistically, we observed that Activin A-treated hiPSC-CMs exhibited impaired diastolic calcium handling with reduced expression of calcium regulatory genes (SERCA2, RYR2, CACNB2). Importantly, when Activin A was inhibited with an anti-Activin A antibody, maladaptive calcium handling and CM contractile dysfunction were abrogated. Therefore, inflammatory cytokines may play a key role by acting on cardiac fibroblasts, causing local upregulation of Activin A that directly acts on CMs to impair contractility. These findings demonstrate that Activin A acts directly on CMs, which may contribute to the cardiac dysfunction seen in aging populations and in patients with heart failure.

Early onset epilepsy and sudden unexpected death in epilepsy with cardiac arrhythmia in mice carrying the early infantile epileptic encephalopathy 47 gain-of-function FHF1(FGF12) missense mutation.

OBJECTIVE: Fibroblast growth factor homologous factors (FHFs) are brain and cardiac sodium channel-binding proteins that modulate channel density and inactivation gating. A recurrent de novo gain-of-function missense mutation in the FHF1(FGF12) gene (p.Arg52His) is associated with early infantile epileptic encephalopathy 47 (EIEE47; Online Mendelian Inheritance in Man database 617166). To determine whether the FHF1 missense mutation is sufficient to cause EIEE and to establish an animal model for EIEE47, we sought to engineer this mutation into mice. METHODS: The Arg52His mutation was introduced into fertilized eggs by CRISPR (clustered regularly interspaced short palindromic repeats) editing to generate Fhf1R52H/F+ mice. Spontaneous epileptiform events in Fhf1R52H/+ mice were assessed by cortical electroencephalography (EEG) and video monitoring. Basal heart rhythm and seizure-induced arrhythmia were recorded by electrocardiography. Modulation of cardiac sodium channel inactivation by FHF1BR52H protein was assayed by voltage-clamp recordings of FHF-deficient mouse cardiomyocytes infected with adenoviruses expressing wild-type FHF1B or FHF1BR52H protein. RESULTS: All Fhf1R52H/+ mice experienced seizure or seizurelike episodes with lethal ending between 12 and 26 days of age. EEG recordings in 19-20-day-old mice confirmed sudden unexpected death in epilepsy (SUDEP) as severe tonic seizures immediately preceding loss of brain activity and death. Within 2-53 s after lethal seizure onset, heart rate abruptly declined from 572 ± 16 bpm to 108 ± 15 bpm, suggesting a parasympathetic surge accompanying seizures that may have contributed to SUDEP. Although ectopic overexpression of FHF1BR52H in cardiomyocytes induced a 15-mV depolarizing shift in voltage of steady-state sodium channel inactivation and slowed the rate of channel inactivation, heart rhythm was normal in Fhf1R52H/+ mice prior to seizure. SIGNIFICANCE: The Fhf1 missense mutation p.Arg52His induces epileptic encephalopathy with full penetrance in mice. Both Fhf1 (p.Arg52His) and Scn8a (p.Asn1768Asp) missense mutations enhance sodium channel Nav 1.6 currents and induce SUDEP with bradycardia in mice, suggesting an FHF1/Nav 1.6 functional axis underlying altered brain sodium channel gating in epileptic encephalopathy.

Effect of left atrial ligation-driven altered inflow hemodynamics on embryonic heart development: clues for prenatal progression of hypoplastic left heart syndrome.

Congenital heart defects (CHDs) are abnormalities in the heart structure present at birth. One important condition is hypoplastic left heart syndrome (HLHS) where severely underdeveloped left ventricle (LV) cannot support systemic circulation. HLHS usually initiates as localized tissue malformations with no underlying genetic cause, suggesting that disturbed hemodynamics contribute to the embryonic development of these defects. Left atrial ligation (LAL) is a surgical procedure on embryonic chick resulting in a phenotype resembling clinical HLHS. In this study, we investigated disturbed hemodynamics and deteriorated cardiac growth following LAL to investigate possible mechanobiological mechanisms for the embryonic development of HLHS. We integrated techniques such as echocardiography, micro-CT and computational fluid dynamics (CFD) for these analyses. Specifically, LAL procedure causes an immediate flow disturbance over atrioventricular (AV) cushions. At later stages after the heart septation, it causes hemodynamic disturbances in LV. As a consequence of the LAL procedure, the left-AV canal and LV volume decrease in size, and in the opposite way, the right-AV canal and right ventricle volume increase. According to our CFD analysis, LAL results in an immediate decrease in the left AV canal WSS levels for 3.5-day (HH21) pre-septated hearts. For 7-day post-septated hearts (HH30), LAL leads to further reduction in WSS levels in the left AV canal, and relatively increased WSS levels in the right AV canal. This study demonstrates the critical importance of the disturbed hemodynamics during the heart valve and ventricle development.

Characterization of Vortex Flow in a Mouse Model of Ventricular Dyssynchrony by Plane-Wave Ultrasound Using Hexplex Processing.

The rodent heart is frequently used to study human cardiovascular disease (CVD). Although advanced cardiovascular ultrasound imaging methods are available for human clinical practice, application of these techniques to small animals remains limited due to the temporal and spatial-resolution demands. Here, an ultrasound vector-flow workflow is demonstrated that enables visualization and quantification of the complex hemodynamics within the mouse heart. Wild type (WT) and fibroblast growth factor homologous factor 2 (FHF2)-deficient mice (Fhf2 KO/Y ), which present with hyperthermia-induced ECG abnormalities highly reminiscent of Brugada syndrome, were used as a mouse model of human CVD. An 18-MHz linear array was used to acquire high-speed (30 kHz), plane-wave data of the left ventricle (LV) while increasing core body temperature up to 41.5 °C. Hexplex (i.e., six output) processing of the raw data sets produced the output of vector-flow estimates (magnitude and phase); B-mode and color-Doppler images; Doppler spectrograms; and local time histories of vorticity and pericardium motion. Fhf2 WT/Y mice had repeatable beat-to-beat cardiac function, including vortex formation during diastole, at all temperatures. In contrast, Fhf2 KO/Y mice displayed dyssynchronous contractile motion that disrupted normal inflow vortex formation and impaired LV filling as temperature rose. The hexplex processing approach demonstrates the ability to visualize and quantify the interplay between hemodynamic and mechanical function in a mouse model of human CVD.

Cardiac Pressure Overload Decreases ETV1 Expression in the Left Atrium, Contributing to Atrial Electrical and Structural Remodeling.

BACKGROUND: Elevated intracardiac pressure attributable to heart failure induces electrical and structural remodeling in the left atrium (LA) that begets atrial myopathy and arrhythmias. The underlying molecular pathways that drive atrial remodeling during cardiac pressure overload are poorly defined. The purpose of this study is to characterize the response of the ETV1 (ETS translocation variant 1) signaling axis in the LA during cardiac pressure overload in humans and mouse models and explore the role of ETV1 in atrial electrical and structural remodeling. METHODS: We performed gene expression profiling in 265 left atrial samples from patients who underwent cardiac surgery. Comparative gene expression profiling was performed between 2 murine models of cardiac pressure overload, transverse aortic constriction banding and angiotensin II infusion, and a genetic model of Etv1 cardiomyocyte-selective knockout (Etv1f/fMlc2aCre/+). RESULTS: Using the Cleveland Clinic biobank of human LA specimens, we found that ETV1 expression is decreased in patients with reduced ejection fraction. Consistent with its role as an important mediator of the NRG1 (Neuregulin 1) signaling pathway and activator of rapid conduction gene programming, we identified a direct correlation between ETV1 expression level and NRG1, ERBB4, SCN5A, and GJA5 levels in human LA samples. In a similar fashion to patients with heart failure, we showed that left atrial ETV1 expression is downregulated at the RNA and protein levels in murine pressure overload models. Comparative analysis of LA RNA sequencing datasets from transverse aortic constriction and angiotensin II-treated mice showed a high Pearson correlation, reflecting a highly ordered process by which the LA undergoes electrical and structural remodeling. Cardiac pressure overload produced a consistent downregulation of ErbB4, Etv1, Scn5a, and Gja5 and upregulation of profibrotic gene programming, which includes Tgfbr1/2, Igf1, and numerous collagen genes. Etv1f/fMlc2aCre/+ mice displayed atrial conduction disease and arrhythmias. Correspondingly, the LA from Etv1f/fMlc2aCre/+ mice showed downregulation of rapid conduction genes and upregulation of profibrotic gene programming, whereas analysis of a gain-of-function ETV1 RNA sequencing dataset from neonatal rat ventricular myocytes transduced with Etv1 showed reciprocal changes. CONCLUSIONS: ETV1 is downregulated in the LA during cardiac pressure overload, contributing to both electrical and structural remodeling.

Ionic Mechanisms of Impulse Propagation Failure in the FHF2-Deficient Heart.

RATIONALE: FHFs (fibroblast growth factor homologous factors) are key regulators of sodium channel (NaV) inactivation. Mutations in these critical proteins have been implicated in human diseases including Brugada syndrome, idiopathic ventricular arrhythmias, and epileptic encephalopathy. The underlying ionic mechanisms by which reduced Nav availability in Fhf2 knockout (Fhf2KO) mice predisposes to abnormal excitability at the tissue level are not well defined. OBJECTIVE: Using animal models and theoretical multicellular linear strands, we examined how FHF2 orchestrates the interdependency of sodium, calcium, and gap junctional conductances to safeguard cardiac conduction. METHODS AND RESULTS: Fhf2KO mice were challenged by reducing calcium conductance (gCaV) using verapamil or by reducing gap junctional conductance (Gj) using carbenoxolone or by backcrossing into a cardiomyocyte-specific Cx43 (connexin 43) heterozygous background. All conditions produced conduction block in Fhf2KO mice, with Fhf2 wild-type (Fhf2WT) mice showing normal impulse propagation. To explore the ionic mechanisms of block in Fhf2KO hearts, multicellular linear strand models incorporating FHF2-deficient Nav inactivation properties were constructed and faithfully recapitulated conduction abnormalities seen in mutant hearts. The mechanisms of conduction block in mutant strands with reduced gCaV or diminished Gj are very different. Enhanced Nav inactivation due to FHF2 deficiency shifts dependence onto calcium current (ICa) to sustain electrotonic driving force, axial current flow, and action potential (AP) generation from cell-to-cell. In the setting of diminished Gj, slower charging time from upstream cells conspires with accelerated Nav inactivation in mutant strands to prevent sufficient downstream cell charging for AP propagation. CONCLUSIONS: FHF2-dependent effects on Nav inactivation ensure adequate sodium current (INa) reserve to safeguard against numerous threats to reliable cardiac impulse propagation.

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.

Fhf2 gene deletion causes temperature-sensitive cardiac conduction failure.

Fever is a highly conserved systemic response to infection dating back over 600 million years. Although conferring a survival benefit, fever can negatively impact the function of excitable tissues, such as the heart, producing cardiac arrhythmias. Here we show that mice lacking fibroblast growth factor homologous factor 2 (FHF2) have normal cardiac rhythm at baseline, but increasing core body temperature by as little as 3 °C causes coved-type ST elevations and progressive conduction failure that is fully reversible upon return to normothermia. FHF2-deficient cardiomyocytes generate action potentials upon current injection at 25 °C but are unexcitable at 40 °C. The absence of FHF2 accelerates the rate of closed-state and open-state sodium channel inactivation, which synergizes with temperature-dependent enhancement of inactivation rate to severely suppress cardiac sodium currents at elevated temperatures. Our experimental and computational results identify an essential role for FHF2 in dictating myocardial excitability and conduction that safeguards against temperature-sensitive conduction failure.

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 Murine Myh6MerCreMer Knock-In Allele Specifically Mediates Temporal Genetic Deletion in Cardiomyocytes after Tamoxifen Induction.

A mouse model that mediates temporal, specific, and efficient myocardial deletion with Cre-LoxP technology will be a valuable tool to determine the function of genes during heart formation. Mhy6 encodes a cardiac muscle specific protein: alpha-myosin heavy chain. Here, we generated a new Myh6-MerCreMer (Myh6(MerCreMer/+)) inducible Cre knock-in mouse by inserting a MerCreMer cassette into the Myh6 start codon. By crossing knock-in mice with Rosa26 reporter lines, we found the Myh6(MerCreMer/+) mice mediate complete Cre-LoxP recombination in cardiomyocytes after tamoxifen induction. X-gal staining and immunohistochemistry analysis revealed that Myh6-driven Cre recombinase was specifically activated in cardiomyocytes at embryonic and adult stages. Furthermore, echocardiography showed that Myh6(MerCreMer/+) mice maintained normal cardiac structure and function before and after tamoxifen administration. These results suggest that the new Myh6(MerCreMer/+) mouse can serve as a robust tool to dissect the roles of genes in heart development and function. Additionally, myocardial progeny during heart development and after cardiac injury can be traced using this mouse line.

Generation of a tamoxifen inducible Tnnt2MerCreMer knock-in mouse model for cardiac studies.

Tnnt2, encoding thin-filament sarcomeric protein cardiac troponin T, plays critical roles in heart development and function in mammals. To develop an inducible genetic deletion strategy in myocardial cells, we generated a new Tnnt2:MerCreMer (Tnnt2(MerCreMer/+)) knock-in mouse. Rosa26 reporter lines were used to examine the specificity and efficiency of the inducible Cre recombinase. We found that Cre was specifically and robustly expressed in the cardiomyocytes at embryonic and adult stages following tamoxifen induction. The knock-in allele on Tnnt2 locus does not impact cardiac function. These results suggest that this new Tnnt2(MerCreMer/+) mouse could be applied towards the temporal genetic deletion of genes of interests in cardiomyocytes with Cre-LoxP technology. The Tnnt2(MerCreMer/+) mouse model also provides a useful tool to trace myocardial lineage during development and repair after cardiac injury.

Scn1b deletion leads to increased tetrodotoxin-sensitive sodium current, altered intracellular calcium homeostasis and arrhythmias in murine hearts.

KEY POINTS: Na(+) current (INa) results from the integrated function of a molecular aggregate (the voltage-gated Na(+) channel complex) that includes the ß subunit family. Mutations or rare variants in Scn1b (encoding the ß1 and ß1B subunits) have been associated with various inherited arrhythmogenic syndromes, including Brugada syndrome and sudden unexpected death in patients with epilepsy. We used Scn1b null mice to understand better the relation between Scn1b expression, and cardiac electrical function. Loss of Scn1b caused, among other effects, increased amplitude of tetrodotoxin-sensitive INa, delayed after-depolarizations, triggered beats, delayed Ca(2+) transients, frequent spontaneous calcium release events and increased susceptibility to polymorphic ventricular arrhythmias. Most alterations in Ca(2+) homeostasis were prevented by 100 nM tetrodotoxin. We propose that life-threatening arrhythmias in patients with mutations in Scn1b, a gene classically defined as ancillary to the Na(+) channel α subunit, can be partly consequent to disrupted intracellular Ca(2+) homeostasis. ABSTRACT: Na(+) current (INa) is determined not only by the properties of the pore-forming voltage-gated Na(+) channel (VGSC) α subunit, but also by the integrated function of a molecular aggregate (the VGSC complex) that includes the VGSC ß subunit family. Mutations or rare variants in Scn1b (encoding the ß1 and ß1B subunits) have been associated with various inherited arrhythmogenic syndromes, including cases of Brugada syndrome and sudden unexpected death in patients with epilepsy. Here, we have used Scn1b null mouse models to understand better the relation between Scn1b expression, and cardiac electrical function. Using a combination of macropatch and scanning ion conductance microscopy we show that loss of Scn1b in juvenile null animals resulted in increased tetrodotoxin-sensitive INa but only in the cell midsection, even before full T-tubule formation; the latter occurred concurrent with increased message abundance for the neuronal Scn3a mRNA, suggesting increased abundance of tetrodotoxin-sensitive NaV 1.3 protein and yet its exclusion from the region of the intercalated disc. Ventricular myocytes from cardiac-specific adult Scn1b null animals showed increased Scn3a message, prolonged action potential repolarization, presence of delayed after-depolarizations and triggered beats, delayed Ca(2+) transients and frequent spontaneous Ca(2+) release events and at the whole heart level, increased susceptibility to polymorphic ventricular arrhythmias. Most alterations in Ca(2+) homeostasis were prevented by 100 nM tetrodotoxin. Our results suggest that life-threatening arrhythmias in patients with mutations in Scn1b, a gene classically defined as ancillary to the Na(+) channel α subunit, can be partly consequent to disrupted intracellular Ca(2+) homeostasis in ventricular myocytes.

Growth and hemodynamics after early embryonic aortic arch occlusion.

The majority of severe clinically significant forms of congenital heart disease (CHD) are associated with great artery lesions, including hypoplastic, double, right or interrupted aortic arch morphologies. While fetal and neonatal interventions are advancing, their potential ability to restore cardiac function, optimal timing, location, and intensity required for intervention remain largely unknown. Here, we combine computational fluid dynamics (CFD) simulations with in vivo experiments to test how individual pharyngeal arch artery hemodynamics alter as a result of local interventions obstructing individual arch artery flow. Simulated isolated occlusions within each pharyngeal arch artery were created with image-derived three-dimensional (3D) reconstructions of normal chick pharyngeal arch anatomy at Hamburger-Hamilton (HH) developmental stages HH18 and HH24. Acute flow redistributions were then computed using in vivo measured subject-specific aortic sinus inflow velocity profiles. A kinematic vascular growth-rendering algorithm was then developed and implemented to test the role of changing local wall shear stress patterns in downstream 3D morphogenesis of arch arteries. CFD simulations predicted that altered pressure gradients and flow redistributions were most sensitive to occlusion of the IVth arches. To evaluate these simulations experimentally, a novel in vivo experimental model of pharyngeal arch occlusion was developed and implemented using two-photon microscopy-guided femtosecond laser-based photodisruption surgery. The right IVth arch was occluded at HH18, and resulting diameter changes were followed for up to 24 h. Pharyngeal arch diameter responses to acute hemodynamic changes were predicted qualitatively but poorly quantitatively. Chronic growth and adaptation to hemodynamic changes, however, were predicted in a subset of arches. Our findings suggest that this complex biodynamic process is governed through more complex forms of mechanobiological vascular growth rules. Other factors in addition to wall shear stress or more complex WSS rules are likely important in the long-term arterial growth and patterning. Combination in silico/experimental platforms are essential for accelerating our understanding and prediction of consequences from embryonic/fetal cardiovascular occlusions and lay the foundation for noninvasive methods to guide CHD diagnosis and fetal intervention.

Isolation and characterization of embryonic stem cell-derived cardiac Purkinje cells.

The cardiac Purkinje fiber network is composed of highly specialized cardiomyocytes responsible for the synchronous excitation and contraction of the ventricles. Computational modeling, experimental animal studies, and intracardiac electrical recordings from patients with heritable and acquired forms of heart disease suggest that Purkinje cells (PCs) may also serve as critical triggers of life-threatening arrhythmias. Nonetheless, owing to the difficulty in isolating and studying this rare population of cells, the precise role of PC in arrhythmogenesis and the underlying molecular mechanisms responsible for their proarrhythmic behavior are not fully characterized. Conceptually, a stem cell-based model system might facilitate studies of PC-dependent arrhythmia mechanisms and serve as a platform to test novel therapeutics. Here, we describe the generation of murine embryonic stem cells (ESC) harboring pan-cardiomyocyte and PC-specific reporter genes. We demonstrate that the dual reporter gene strategy may be used to identify and isolate the rare ESC-derived PC (ESC-PC) from a mixed population of cardiogenic cells. ESC-PC display transcriptional signatures and functional properties, including action potentials, intracellular calcium cycling, and chronotropic behavior comparable to endogenous PC. Our results suggest that stem-cell derived PC are a feasible new platform for studies of developmental biology, disease pathogenesis, and screening for novel antiarrhythmic therapies.

PCP4 regulates Purkinje cell excitability and cardiac rhythmicity.

Cardiac Purkinje cells are important triggers of ventricular arrhythmias associated with heritable and acquired syndromes; however, the mechanisms responsible for this proarrhythmic behavior are incompletely understood. Here, through transcriptional profiling of genetically labeled cardiomyocytes, we identified expression of Purkinje cell protein-4 (Pcp4), a putative regulator of calmodulin and Ca2+/calmodulin-dependent kinase II (CaMKII) signaling, exclusively within the His-Purkinje network. Using Pcp4-null mice and acquired cardiomyopathy models, we determined that reduced expression of PCP4 is associated with CaMKII activation, abnormal electrophysiology, dysregulated intracellular calcium handling, and proarrhythmic behavior in isolated Purkinje cells. Pcp4-null mice also displayed profound autonomic dysregulation and arrhythmic behavior in vivo. Together, these results demonstrate that PCP4 regulates cardiac excitability through both Purkinje cell-autonomous and central mechanisms and identify this modulator of CaMKII signaling as a potential arrhythmia-susceptibility candidate.

Computational fluid dynamics of developing avian outflow tract heart valves.

Hemodynamic forces play an important role in sculpting the embryonic heart and its valves. Alteration of blood flow patterns through the hearts of embryonic animal models lead to malformations that resemble some clinical congenital heart defects, but the precise mechanisms are poorly understood. Quantitative understanding of the local fluid forces acting in the heart has been elusive because of the extremely small and rapidly changing anatomy. In this study, we combine multiple imaging modalities with computational simulation to rigorously quantify the hemodynamic environment within the developing outflow tract (OFT) and its eventual aortic and pulmonary valves. In vivo Doppler ultrasound generated velocity profiles were applied to Micro-Computed Tomography generated 3D OFT lumen geometries from Hamburger-Hamilton (HH) stage 16-30 chick embryos. Computational fluid dynamics simulation initial conditions were iterated until local flow profiles converged with in vivo Doppler flow measurements. Results suggested that flow in the early tubular OFT (HH16 and HH23) was best approximated by Poiseuille flow, while later embryonic OFT septation (HH27, HH30) was mimicked by plug flow conditions. Peak wall shear stress (WSS) values increased from 18.16 dynes/cm(2) at HH16 to 671.24 dynes/cm(2) at HH30. Spatiotemporally averaged WSS values also showed a monotonic increase from 3.03 dynes/cm(2) at HH16 to 136.50 dynes/cm(2) at HH30. Simulated velocity streamlines in the early heart suggest a lack of mixing, which differed from classical ink injections. Changes in local flow patterns preceded and correlated with key morphogenetic events such as OFT septation and valve formation. This novel method to quantify local dynamic hemodynamics parameters affords insight into sculpting role of blood flow in the embryonic heart and provides a quantitative baseline dataset for future research.

Hemodynamic patterning of the avian atrioventricular valve.

In this study, we develop an innovative approach to rigorously quantify the evolving hemodynamic environment of the atrioventricular (AV) canal of avian embryos. Ultrasound generated velocity profiles were imported into Micro-Computed Tomography generated anatomically precise cardiac geometries between Hamburger-Hamilton (HH) stages 17 and 30. Computational fluid dynamic simulations were then conducted and iterated until results mimicked in vivo observations. Blood flow in tubular hearts (HH17) was laminar with parallel streamlines, but strong vortices developed simultaneous with expansion of the cushions and septal walls. For all investigated stages, highest wall shear stresses (WSS) are localized to AV canal valve-forming regions. Peak WSS increased from 19.34 dynes/cm(2) at HH17 to 287.18 dynes/cm(2) at HH30, but spatiotemporally averaged WSS became 3.62 dynes/cm(2) for HH17 to 9.11 dynes/cm(2) for HH30. Hemodynamic changes often preceded and correlated with morphological changes. These results establish a quantitative baseline supporting future hemodynamic analyses and interpretations.