Multisite Assembly of Gateway Induced Clones (MAGIC): a flexible cloning toolbox with diverse applications in vertebrate model systems.

Here we present the Multisite Assembly of Gateway Induced Clones (MAGIC) system, which harnesses site-specific recombination-based cloning via Gateway technology for rapid, modular assembly of between 1 and 3 "Entry" vector components, all into a fourth, standard high copy "Destination" plasmid backbone. The MAGIC toolkit spans a range of in vitro and in vivo uses, from directing tunable gene expression, to driving simultaneous expression of microRNAs and fluorescent reporters, to enabling site-specific recombinase-dependent gene expression. All MAGIC system components are directly compatible with existing multisite gateway Tol2 systems currently used in zebrafish, as well as existing eukaryotic cell culture expression Destination plasmids, and available mammalian lentiviral and adenoviral Destination vectors, allowing rapid cross-species experimentation. Moreover, herein we describe novel vectors with flanking piggyBac transposon elements for stable genomic integration in vitro or in vivo when used with piggyBac transposase. Collectively, the MAGIC system facilitates transgenesis in cultured mammalian cells, electroporated mouse and chick embryos, as well as in injected zebrafish embryos, enabling the rapid generation of innovative DNA constructs for biological research due to a shared, common plasmid platform.

MEK signaling represents a viable therapeutic vulnerability of KRAS-driven somatic brain arteriovenous malformations.

Brain arteriovenous malformations (bAVMs) are direct connections between arteries and veins that remodel into a complex nidus susceptible to rupture and hemorrhage. Most sporadic bAVMs feature somatic activating mutations within KRAS, and endothelial-specific expression of the constitutively active variant KRASG12D models sporadic bAVM in mice. By leveraging 3D-based micro-CT imaging, we demonstrate that KRASG12D-driven bAVMs arise in stereotypical anatomical locations within the murine brain, which coincide with high endogenous Kras expression. We extend these analyses to show that a distinct variant, KRASG12C, also generates bAVMs in predictable locations. Analysis of 15,000 human patients revealed that, similar to murine models, bAVMs preferentially occur in distinct regions of the adult brain. Furthermore, bAVM location correlates with hemorrhagic frequency. Quantification of 3D imaging revealed that G12D and G12C alter vessel density, tortuosity, and diameter within the mouse brain. Notably, aged G12D mice feature increased lethality, as well as impaired cognition and motor function. Critically, we show that pharmacological blockade of the downstream kinase, MEK, after lesion formation ameliorates KRASG12D-driven changes in the murine cerebrovasculature and may also impede bAVM progression in human pediatric patients. Collectively, these data show that distinct KRAS variants drive bAVMs in similar patterns and suggest MEK inhibition represents a non-surgical alternative therapy for sporadic bAVM.

Protocol for optical, aqueous-based clearing of murine tissues using EZ Clear.

Tissue clearing is an essential prerequisite for 3D volumetric imaging of larger tissues or organs. Here, we present a detailed protocol for optical, aqueous-based clearing of adult murine tissues using EZ Clear. We describe steps to ensure successful perfusion and fixation of organs from the adult mouse and supply guidelines for optimal lipid removal, refractive index matching, and tissue clearing. Finally, we provide imaging parameters for visualizing both exogenous perfused fluorescent dyes and endogenous fluorescence reporters in the adult mouse. For complete details on the use and execution of this protocol, please refer to Hsu et al.1.

GPU-Accelerated RSF Level Set Evolution for Large-Scale Microvascular Segmentation.

Microvascular networks are challenging to model because these structures are currently near the diffraction limit for most advanced three-dimensional imaging modalities, including confocal and light sheet microscopy. This makes semantic segmentation difficult, because individual components of these networks fluctuate within the confines of individual pixels. Level set methods are ideally suited to solve this problem by providing surface and topological constraints on the resulting model, however these active contour techniques are extremely time intensive and impractical for terabyte-scale images. We propose a reformulation and implementation of the region-scalable fitting (RSF) level set model that makes it amenable to three-dimensional evaluation using both single-instruction multiple data (SIMD) and single-program multiple-data (SPMD) parallel processing. This enables evaluation of the level set equation on independent regions of the data set using graphics processing units (GPUs), making large-scale segmentation of high-resolution networks practical and inexpensive. We tested this 3D parallel RSF approach on multiple data sets acquired using state-of-the-art imaging techniques to acquire microvascular data, including micro-CT, light sheet fluorescence microscopy (LSFM) and milling microscopy. To assess the performance and accuracy of the RSF model, we conducted a Monte-Carlo-based validation technique to compare results to other segmentation methods. We also provide a rigorous profiling to show the gains in processing speed leveraging parallel hardware. This study showcases the practical application of the RSF model, emphasizing its utility in the challenging domain of segmenting large-scale high-topology network structures with a particular focus on building microvascular models.

Endothelial cells adopt a pro-reparative immune responsive signature during cardiac injury.

Modulation of the heart's immune microenvironment is crucial for recovery after ischemic events such as myocardial infarction (MI). Endothelial cells (ECs) can have immune regulatory functions; however, interactions between ECs and the immune environment in the heart after MI remain poorly understood. We identified an EC-specific IFN responsive and immune regulatory gene signature in adult and pediatric heart failure (HF) tissues. Single-cell transcriptomic analysis of murine hearts subjected to MI uncovered an EC population (IFN-ECs) with immunologic gene signatures similar to those in human HF. IFN-ECs were enriched in regenerative-stage mouse hearts and expressed genes encoding immune responsive transcription factors (Irf7, Batf2, and Stat1). Single-cell chromatin accessibility studies revealed an enrichment of these TF motifs at IFN-EC signature genes. Expression of immune regulatory ligand genes by IFN-ECs suggests bidirectional signaling between IFN-ECs and macrophages in regenerative-stage hearts. Our data suggest that ECs may adopt immune regulatory signatures after cardiac injury to accompany the reparative response. The presence of these signatures in human HF and murine MI models suggests a potential role for EC-mediated immune regulation in responding to stress induced by acute injury in MI and chronic adverse remodeling in HF.

The Role and Therapeutic Implications of Inflammation in the Pathogenesis of Brain Arteriovenous Malformations.

Brain arteriovenous malformations (bAVMs) are focal vascular lesions composed of abnormal vascular channels without an intervening capillary network. As a result, high-pressure arterial blood shunts directly into the venous outflow system. These high-flow, low-resistance shunts are composed of dilated, tortuous, and fragile vessels, which are prone to rupture. BAVMs are a leading cause of hemorrhagic stroke in children and young adults. Current treatments for bAVMs are limited to surgery, embolization, and radiosurgery, although even these options are not viable for ~20% of AVM patients due to excessive risk. Critically, inflammation has been suggested to contribute to lesion progression. Here we summarize the current literature discussing the role of the immune system in bAVM pathogenesis and lesion progression, as well as the potential for targeting inflammation to prevent bAVM rupture and intracranial hemorrhage. We conclude by proposing that a dysfunctional endothelium, which harbors the somatic mutations that have been shown to give rise to sporadic bAVMs, may drive disease development and progression by altering the immune status of the brain.

EZ Clear for simple, rapid, and robust mouse whole organ clearing.

Tissue clearing for whole organ cell profiling has revolutionized biology and imaging for exploration of organs in three-dimensional space without compromising tissue architecture. But complicated, laborious procedures, or expensive equipment, as well as the use of hazardous, organic solvents prevent the widespread adoption of these methods. Here, we report a simple and rapid tissue clearing method, EZ Clear, that can clear whole adult mouse organs in 48 hr in just three simple steps. Samples stay at room temperature and remain hydrated throughout the clearing process, preserving endogenous and synthetic fluorescence, without altering sample size. After wholemount clearing and imaging, samples processed with EZ Clear can be subjected to downstream applications, such as tissue embedding and cryosectioning followed by standard histology or immunofluorescent staining without loss of fluorescence signal from endogenous or synthetic reporters. Furthermore, we demonstrate that wholemount adult mouse brains processed with EZ Clear can be successfully immunolabeled for fluorescent imaging while still retaining signal from endogenous fluorescent reporters. Overall, the simplicity, speed, and flexibility of EZ Clear make it easy to adapt and implement in diverse imaging modalities in biomedical research.

A human model of arteriovenous malformation (AVM)-on-a-chip reproduces key disease hallmarks and enables drug testing in perfused human vessel networks.

Brain arteriovenous malformations (AVMs) are a disorder wherein abnormal, enlarged blood vessels connect arteries directly to veins, without an intervening capillary bed. AVMs are one of the leading causes of hemorrhagic stroke in children and young adults. Most human sporadic brain AVMs are associated with genetic activating mutations in the KRAS gene. Our goal was to develop an in vitro model that would allow for simultaneous morphological and functional phenotypic data capture in real time during AVM disease progression. By generating human endothelial cells harboring a clinically relevant mutation found in most human patients (activating mutations within the small GTPase KRAS) and seeding them in a dynamic microfluidic cell culture system that enables vessel formation and perfusion, we demonstrate that vessels formed by KRAS4AG12V mutant endothelial cells (ECs) were significantly wider and more leaky than vascular beds formed by wild-type ECs, recapitulating key structural and functional hallmarks of human AVM pathogenesis. Immunofluorescence staining revealed a breakdown of adherens junctions in mutant KRAS vessels, leading to increased vascular permeability, a hallmark of hemorrhagic stroke. Finally, pharmacological blockade of MEK kinase activity, but not PI3K inhibition, improved endothelial barrier function (decreased permeability) without affecting vessel diameter. Collectively, our studies describe the creation of human KRAS-dependent AVM-like vessels in vitro in a self-assembling microvessel platform that is amenable to phenotypic observation and drug delivery.

Bone marrow endothelial dysfunction promotes myeloid cell expansion in cardiovascular disease.

Abnormal hematopoiesis advances cardiovascular disease by generating excess inflammatory leukocytes that attack the arteries and the heart. The bone marrow niche regulates hematopoietic stem cell proliferation and hence the systemic leukocyte pool, but whether cardiovascular disease affects the hematopoietic organ's microvasculature is unknown. Here we show that hypertension, atherosclerosis and myocardial infarction (MI) instigate endothelial dysfunction, leakage, vascular fibrosis and angiogenesis in the bone marrow, altogether leading to overproduction of inflammatory myeloid cells and systemic leukocytosis. Limiting angiogenesis with endothelial deletion of Vegfr2 (encoding vascular endothelial growth factor (VEGF) receptor 2) curbed emergency hematopoiesis after MI. We noted that bone marrow endothelial cells assumed inflammatory transcriptional phenotypes in all examined stages of cardiovascular disease. Endothelial deletion of Il6 or Vcan (encoding versican), genes shown to be highly expressed in mice with atherosclerosis or MI, reduced hematopoiesis and systemic myeloid cell numbers in these conditions. Our findings establish that cardiovascular disease remodels the vascular bone marrow niche, stimulating hematopoiesis and production of inflammatory leukocytes.

FOXO1 represses sprouty 2 and sprouty 4 expression to promote arterial specification and vascular remodeling in the mouse yolk sac.

Establishing a functional circulatory system is required for post-implantation development during murine embryogenesis. Previous studies in loss-of-function mouse models showed that FOXO1, a Forkhead family transcription factor, is required for yolk sac (YS) vascular remodeling and survival beyond embryonic day (E) 11. Here, we demonstrate that at E8.25, loss of Foxo1 in Tie2-cre expressing cells resulted in increased sprouty 2 (Spry2) and Spry4 expression, reduced arterial gene expression and reduced Kdr (also known as Vegfr2 and Flk1) transcripts without affecting overall endothelial cell identity, survival or proliferation. Using a Dll4-BAC-nlacZ reporter line, we found that one of the earliest expressed arterial genes, delta like 4, is significantly reduced in Foxo1 mutant YS without being substantially affected in the embryo proper. We show that FOXO1 binds directly to previously identified Spry2 gene regulatory elements (GREs) and newly identified, evolutionarily conserved Spry4 GREs to repress their expression. Furthermore, overexpression of Spry4 in transient transgenic embryos largely recapitulates the reduced expression of arterial genes seen in conditional Foxo1 mutants. Together, these data reveal a novel role for FOXO1 as a key transcriptional repressor regulating both pre-flow arterial specification and subsequent vessel remodeling within the murine YS.

RBFOX2 is required for establishing RNA regulatory networks essential for heart development.

Human genetic studies identified a strong association between loss of function mutations in RBFOX2 and hypoplastic left heart syndrome (HLHS). There are currently no Rbfox2 mouse models that recapitulate HLHS. Therefore, it is still unknown how RBFOX2 as an RNA binding protein contributes to heart development. To address this, we conditionally deleted Rbfox2 in embryonic mouse hearts and found profound defects in cardiac chamber and yolk sac vasculature formation. Importantly, our Rbfox2 conditional knockout mouse model recapitulated several molecular and phenotypic features of HLHS. To determine the molecular drivers of these cardiac defects, we performed RNA-sequencing in Rbfox2 mutant hearts and identified dysregulated alternative splicing (AS) networks that affect cell adhesion to extracellular matrix (ECM) mediated by Rho GTPases. We identified two Rho GTPase cycling genes as targets of RBFOX2. Modulating AS of these two genes using antisense oligos led to cell cycle and cell-ECM adhesion defects. Consistently, Rbfox2 mutant hearts displayed cell cycle defects and inability to undergo endocardial-mesenchymal transition, processes dependent on cell-ECM adhesion and that are seen in HLHS. Overall, our work not only revealed that loss of Rbfox2 leads to heart development defects resembling HLHS, but also identified RBFOX2-regulated AS networks that influence cell-ECM communication vital for heart development.

Computational insights on coronary artery function.

Collateral arteries may act as natural bypasses that reduce hypoperfusion after a coronary blockage. 3D imaging of neonatal and adult mouse hearts, plus human fetal and diseased adult hearts, is now used to computationally predict flow within the heart, and understand the cardioprotective role of collateral arteries in vivo.

Identification and characterization of conserved noncoding cis-regulatory elements that impact Mecp2 expression and neurological functions.

While changes in MeCP2 dosage cause Rett syndrome (RTT) and MECP2 duplication syndrome (MDS), its transcriptional regulation is poorly understood. Here, we identified six putative noncoding regulatory elements of Mecp2, two of which are conserved in humans. Upon deletion in mice and human iPSC-derived neurons, these elements altered RNA and protein levels in opposite directions and resulted in a subset of RTT- and MDS-like behavioral deficits in mice. Our discovery provides insight into transcriptional regulation of Mecp2/MECP2 and highlights genomic sites that could serve as diagnostic and therapeutic targets in RTT or MDS.

Identification of diverse tumor endothelial cell populations in malignant glioma.

BACKGROUND: Glioblastoma is the most common and aggressive type of primary brain tumor, as most patients succumb to the disease less than two years after diagnosis. Critically, studies demonstrate that glioma recruits surrounding blood vessels, while some work suggests that tumor stem cells themselves directly differentiate into endothelial cells, yet the molecular and cellular dynamics of the endothelium in glioma are poorly characterized. The goal of this study was to establish molecular and morphological benchmarks for tumor associated vessels (TAVs) and tumor derived endothelial cells (TDECs) during glioblastoma progression. METHODS: Using In-Utero Electroporation and CRISPR/Cas9 genome engineering to generate a native, immunocompetent mouse model of glioma, we characterized vascular-tumor dynamics in three dimensions during tumor progression. We employed bulk and single-cell RNA-Sequencing to elucidate the relationship between TAVs and TDECs. We confirmed our findings in a patient derived orthotopic xenograft (PDOX) model. RESULTS: Using a mouse model of glioma, we identified progressive alteration of vessel function and morphogenesis over time. We also showed in our mouse model that TDECs are a rare subpopulation that contributes to vessels within the tumor, albeit to a limited degree. Furthermore, transcriptional profiling demonstrates that both TAVs and TDECs are molecularly distinct, and both populations feature extensive molecular heterogeneity. Finally, the distinct molecular signatures of these heterogeneous populations are also present in human glioma. CONCLUSIONS: Our findings show extensive endothelial heterogeneity within the tumor and tumor microenvironment and provide insights into the diverse cellular and molecular mechanisms that drive glioma vascularization and angiogenesis during tumorigenesis.

Somatic Gain of KRAS Function in the Endothelium Is Sufficient to Cause Vascular Malformations That Require MEK but Not PI3K Signaling.

RATIONALE: We previously identified somatic activating mutations in the KRAS (Kirsten rat sarcoma viral oncogene homologue) gene in the endothelium of the majority of human sporadic brain arteriovenous malformations; a disorder characterized by direct connections between arteries and veins. However, whether this genetic abnormality alone is sufficient for lesion formation, as well as how active KRAS signaling contributes to arteriovenous malformations, remains unknown. OBJECTIVE: To establish the first in vivo models of somatic KRAS gain of function in the endothelium in both mice and zebrafish to directly observe the phenotypic consequences of constitutive KRAS activity at a cellular level in vivo, and to test potential therapeutic interventions for arteriovenous malformations. METHODS AND RESULTS: Using both postnatal and adult mice, as well as embryonic zebrafish, we demonstrate that endothelial-specific gain of function mutations in Kras (G12D or G12V) are sufficient to induce brain arteriovenous malformations. Active KRAS signaling leads to altered endothelial cell morphogenesis and increased cell size, ectopic sprouting, expanded vessel lumen diameter, and direct connections between arteries and veins. Furthermore, we show that these lesions are not associated with altered endothelial growth dynamics or a lack of proper arteriovenous identity but instead seem to feature exuberant angiogenic signaling. Finally, we demonstrate that KRAS-dependent arteriovenous malformations in zebrafish are refractory to inhibition of the downstream effector PI3K but instead require active MEK (mitogen-activated protein kinase kinase 1) signaling. CONCLUSIONS: We demonstrate that active KRAS expression in the endothelium is sufficient for brain arteriovenous malformations, even in the setting of uninjured adult vasculature. Furthermore, the finding that KRAS-dependent lesions are reversible in zebrafish suggests that MEK inhibition may represent a promising therapeutic treatment for arteriovenous malformation patients. Graphical Abstract: A graphical abstract is available for this article.

S1PR1 regulates the quiescence of lymphatic vessels by inhibiting laminar shear stress-dependent VEGF-C signaling.

During the growth of lymphatic vessels (lymphangiogenesis), lymphatic endothelial cells (LECs) at the growing front sprout by forming filopodia. Those tip cells are not exposed to circulating lymph, as they are not lumenized. In contrast, LECs that trail the growing front are exposed to shear stress, become quiescent, and remodel into stable vessels. The mechanisms that coordinate the opposed activities of lymphatic sprouting and maturation remain poorly understood. Here, we show that the canonical tip cell marker Delta-like 4 (DLL4) promotes sprouting lymphangiogenesis by enhancing VEGF-C/VEGF receptor 3 (VEGFR3) signaling. However, in lumenized lymphatic vessels, laminar shear stress (LSS) inhibits the expression of DLL4, as well as additional tip cell markers. Paradoxically, LSS also upregulates VEGF-C/VEGFR3 signaling in LECs, but sphingosine 1-phosphate receptor 1 (S1PR1) activity antagonizes LSS-mediated VEGF-C signaling to promote lymphatic vascular quiescence. Correspondingly, S1pr1 loss in LECs induced lymphatic vascular hypersprouting and hyperbranching, which could be rescued by reducing Vegfr3 gene dosage in vivo. In addition, S1PR1 regulates lymphatic vessel maturation by inhibiting RhoA activity to promote membrane localization of the tight junction molecule claudin-5. Our findings suggest a potentially new paradigm in which LSS induces quiescence and promotes the survival of LECs by downregulating DLL4 and enhancing VEGF-C signaling, respectively. S1PR1 dampens LSS/VEGF-C signaling, thereby preventing sprouting from quiescent lymphatic vessels. These results also highlight the distinct roles that S1PR1 and DLL4 play in LECs when compared with their known roles in the blood vasculature.

Epithelial Vegfa Specifies a Distinct Endothelial Population in the Mouse Lung.

The lung microvasculature is essential for gas exchange and commonly considered homogeneous. We show that VEGFA from the epithelium is required for a distinct endothelial cell (EC) population in the mouse lung. Vegfa is predominantly expressed by alveolar type 1 (AT1) cells and locally required to specify a subset of ECs. Single-cell RNA sequencing (scRNA-seq) reveals that ∼15% of lung ECs are transcriptionally distinct-marked by Carbonic anhydrase 4 (Car4)-and arise from bulk ECs, as suggested by trajectory analysis. Car4 ECs have extensive cellular projections and are separated from AT1 cells by a limited basement membrane without intervening pericytes. Car4 ECs are specifically lost upon epithelial Vegfa deletion; without Car4 ECs, the alveolar space is aberrantly enlarged despite the normal appearance of myofibroblasts. Lung Car4 ECs and retina tip ECs have common and distinct features. These findings support a signaling role of AT1 cells and shed light on alveologenesis.

Development of barium-based low viscosity contrast agents for micro CT vascular casting: Application to 3D visualization of the adult mouse cerebrovasculature.

Recent advances in three-dimensional (3D) fluorescence microscopy offer the ability to image the entire vascular network in entire organs, or even whole animals. However, these imaging modalities rely on either endogenous fluorescent reporters or involved immunohistochemistry protocols, as well as optical clearing of the tissue and refractive index matching. Conversely, X-ray-based 3D imaging modalities, such as micro CT, can image non-transparent samples, at high resolution, without requiring complicated or expensive immunolabeling and clearing protocols, or fluorescent reporters. Here, we compared two "homemade" barium-based contrast agents to the field standard, lead-containing Microfil, for micro-computed tomography (micro CT) imaging of the adult mouse cerebrovasculature. The perfusion pressure required for uniform vessel filling was significantly lower with the barium-based contrast agents compared to the polymer-based Microfil. Accordingly, the barium agents showed no evidence of vascular distension or rupture, common problems associated with Microfil. Compellingly, perfusion of an aqueous BaCl2 /gelatin mixture yielded equal or superior visualization of the cerebrovasculature by micro CT compared to Microfil. However, phosphate-containing buffers and fixatives were incompatible with BaCl2 due to the formation of unwanted precipitates. X-ray attenuation of the vessels also decreased overtime, as the BaCl2 appeared to gradually diffuse into surrounding tissues. A second, unique formulation composed of BaSO4 microparticles, generated in-house by mixing BaCl2 and MgSO4 , suffered none of these drawbacks. These microparticles, however, were unable to pass small diameter capillary vessels, conveniently labeling only the arterial cerebrovasculature. In summary, we present an affordable, robust, low pressure, non-toxic, and straightforward methodology for 3D visualization of the cerebrovasculature.

FishNET: An automated relational database for zebrafish colony management.

The zebrafish Danio rerio is a powerful model system to study the genetics of development and disease. However, maintenance of zebrafish husbandry records is both time intensive and laborious, and a standardized way to manage and track the large amount of unique lines in a given laboratory or centralized facility has not been embraced by the field. Here, we present FishNET, an intuitive, open-source, relational database for managing data and information related to zebrafish husbandry and maintenance. By creating a "virtual facility," FishNET enables users to remotely inspect the rooms, racks, tanks, and lines within a given facility. Importantly, FishNET scales from one laboratory to an entire facility with several laboratories to multiple facilities, generating a cohesive laboratory and community-based platform. Automated data entry eliminates confusion regarding line nomenclature and streamlines maintenance of individual lines, while flexible query forms allow researchers to retrieve database records based on user-defined criteria. FishNET also links associated embryonic and adult biological samples with data, such as genotyping results or confocal images, to enable robust and efficient colony management and storage of laboratory information. A shared calendar function with email notifications and automated reminders for line turnover, automated tank counts, and census reports promote communication with both end users and administrators. The expected benefits of FishNET are improved vivaria efficiency, increased quality control for experimental numbers, and flexible data reporting and retrieval. FishNET's easy, intuitive record management and open-source, end-user-modifiable architecture provides an efficient solution to real-time zebrafish colony management for users throughout a facility and institution and, in some cases, across entire research hubs.

A cellular atlas of Pitx2-dependent cardiac development.

The Pitx2 gene encodes a homeobox transcription factor that is required for mammalian development. Disruption of PITX2 expression in humans causes congenital heart diseases and is associated with atrial fibrillation; however, the cellular and molecular processes dictated by Pitx2 during cardiac ontogeny remain unclear. To characterize the role of Pitx2 during murine heart development we sequenced over 75,000 single cardiac cell transcriptomes between two key developmental timepoints in control and Pitx2 null embryos. We found that cardiac cell composition was dramatically altered in mutants at both E10.5 and E13.5. Interestingly, the differentiation dynamics of both anterior and posterior second heart field-derived progenitor cells were disrupted in Pitx2 mutants. We also uncovered evidence for defects in left-right asymmetry within atrial cardiomyocyte populations. Furthermore, we were able to detail defects in cardiac outflow tract and valve development associated with Pitx2 Our findings offer insight into Pitx2 function and provide a compilation of gene expression signatures for further detailing the complexities of heart development that will serve as the foundation for future studies of cardiac morphogenesis, congenital heart disease and arrhythmogenesis.