The level of HAND1 controls the specification of multipotent cardiac and extraembryonic progenitors

The level of HAND1 controls the specification of multipotent cardiac and extraembryonic progenitors | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results The level of HAND1 controls the specification of multipotent cardiac and extraembryonic progenitors Adam T Lynch , Naomi Phillips , Megan Douglas , Marta Dorgnach , I-Hsuan Lin , Antony D Adamson , Zoulfia Darieva , Jessica Whittle , Neil A Hanley , Nicoletta Bobola , Matthew J Birket doi: https://doi.org/10.1101/2024.08.15.607916 Adam T Lynch 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Naomi Phillips 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Megan Douglas 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Marta Dorgnach 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site I-Hsuan Lin 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Antony D Adamson 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Zoulfia Darieva 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Jessica Whittle 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Neil A Hanley 1 Faculty of Biology, Medicine and Health, University of Manchester , UK 2 College of Medicine & Health, University of Birmingham , Edgbaston, Birmingham, B15 2TT, UK 3 University Hospitals Birmingham NHS Foundation Trust , Birmingham, B15 2GW, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Nicoletta Bobola 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Matthew J Birket 1 Faculty of Biology, Medicine and Health, University of Manchester , UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: matthew.birket{at}manchester.ac.uk Abstract Full Text Info/History Metrics Supplementary material Data/Code Preview PDF Summary Diverse sets of progenitors contribute to the development of the embryonic heart, but the mechanisms of their specification have remained elusive. Here, using a human pluripotent stem cell (hPSC) model, we deciphered cardiac and non-cardiac lineage trajectories in differentiation and identified transcription factors underpinning cell specification, identity and function. We discovered a concentration-dependent, fate determining function in mesodermal progenitors for the basic helix-loop-helix transcription factor HAND1 and uncovered its gene regulatory network. At low level, HAND1 directs differentiation towards multipotent juxta-cardiac field progenitors able to make cardiomyocytes and epicardial cells, whereas at high level it promotes the development of extraembryonic mesoderm. Importantly, HAND1-low progenitors can be propagated in their multipotent state. This detailed mechanistic insight into human development will accelerate the delivery of effective disease modelling, including for congenital heart disease, and cell therapy-based regenerative medicine. Introduction Development of the mammalian four-chamber heart is highly complex and depends on the specification and exquisitely regulated, phased differentiation of multiple populations of progenitor cells to a broad final array of muscular (cardiomyocytes) and non-muscular lineages. 1 Errors in any part of this process have been catalogued in human congenital heart disease and have illustrated necessary aspects of genetic regulation. 2 , 3 However, how a number of the critical early human embryonic fate decisions are made has remained elusive. During gastrulation, cell diversity along the primitive streak is largely enabled by differential intensity and duration of BMP and Nodal extracellular signaling. 4 Cardiogenic mesodermal progenitors arise under the combined actions of BMP, Nodal/Activin, and WNT signaling and migrate from the primitive streak in two broad waves. 5 – 7 The initial wave forms the first heart field (FHF) and leads to the linear heart tube. Subsequently, second heart field (SHF) progenitors contribute progressively to the heart’s reshaping into its final four-chamber structure. Within this framework, transcription factors (TFs) controlling the identity of terminally differentiated cells such as cardiomyocytes have been well described. However, how an overlapping set of TFs specify progenitors, including those retaining multipotency, is poorly understood. Recently, single cell genomics, including directly in human embryos, has reinforced the idea that progenitors contributing to the FHF and SHF are actually made up of subpopulations with discrete identities and fate assignments, 8 – 10 such as the juxta-cardiac field (JCF) which overlaps the FHF and is marked by the basic helix-loop-helix (bHLH) TF, Hand1. The JCF gives rise to cardiomyocytes and epicardial cells in the heart but also contributes cells to the extraembryonic mesoderm. 11 , 12 Hand1 is functionally important in each cell type in this array. 13 – 18 However, the underlying mechanisms that determine lineage commitment and control multipotency within progenitor cells remain unresolved, hampering effective regenerative medicine aimed at accurately programming specific human cardiac cell types. In this study we aimed to address this knowledge gap by using human pluripotent stem cells (hPSCs) to model a wide repertoire of cardiac cell development and explore whether a multipotent embryonic progenitor state could be propagated in vitro. Specifically, by using lineage tracking and single cell genomics, we aimed to discover the gene regulatory networks (GRNs) that were operational in specifying the different cardiac cell fates, and mechanistically support advances in the field where 3D models have been used to promote the development of multiple heart cell types. 19 – 21 The outcome established an upstream role for HAND1 in programming the specification and GRNs of multipotent progenitors within the JCF and related extraembryonic mesoderm. The data supports a TF concentration-dependent mechanism of cell fate determination. Results Multi-lineage cell diversity can be achieved in 3D by modulating the duration of Activin and BMP signaling We hypothesized that mimicking the patterning of the primitive streak would achieve diverse differentiation including endoderm and non-cardiac embryonic and extraembryonic mesoderm alongside cardiogenic mesoderm. We adapted our previously characterized cardiomyocyte differentiation model in hPSCs 22 and modulated the duration of Activin and BMP signaling with small molecule inhibitors ( Figure 1A ). In the continued presence of Activin A and BMP4, the ALK5 inhibitor SB431542 (SB) and the ALK1-3 inhibitor DMH1 terminated receptor-regulated SMAD2/3 and SMAD1/5/8 phosphorylation respectively ( Figure 1B ). R-SMADs compete for common partner SMAD4, therefore we hypothesized these inhibitions would enhance the downstream activity of the uninhibited pathway. In line with this, SB promoted development of posterior primitive streak-like mesoderm with increased TBXT , MESP1, and HOXB1 expression compared to the control, whereas DMH1 suppressed MESP1 and promoted anterior primitive streak-like mesoderm and endoderm development with increased MIXL1 and FOXA2 ( Figure 1C ). To deconvolute this developmental diversity, we created a dual cardiac-endoderm fluorescent reporter hESC line by taking an NKX2-5-GFP reporter hESC line and introducing a knock-in SOX17 - T2A-dTomato ( Figure S1A ). With DMH1, SOX17-Tom expression became prominent at day 4 of differentiation and persisted over several days ( Figures 1D–G ). By day 8, NKX2-5-GFP was also detectable with DMH1 whereas it was abolished by a saturating dose of SB. We titrated the small molecules and measured the reporters in concert with surface markers PDGFRα (broad mesoderm) and CD235 (anterior-like mesoderm) ( Figure 1E ). 21 At day 4, PDGFRα and CD235 showed an inverse pattern through the signaling gradient ( Figures 1E and 1F ). At day 7.5, SOX17-Tom could be detected simultaneously with NKX2-5-GFP ( Figures 1E and 1G ) and while mostly mutually exclusive, a population of SOX17+ NKX2-5+ cells was also detected ( Figure 1E ). To help clarify the identity of this dual population we introduced a knock-in FOXA2-T2A-mTagBFP2 reporter ( Figure S1B ), an endoderm marker also potentially expressed in a population of cardiac progenitors. 9 , 23 In this triple reporter, FOXA2-BFP was entirely co-expressed with SOX17-Tom at day 7.5 when NKX2-5-GFP+ cells emerged ( Figure S1C ), consistent with data in early human embryos, 10 and when isolated at day 5, SOX17+ FOXA2+ cells generated very few cardiomyocytes suggesting these were largely not cardiac progenitors ( Figures S1D and S1E ). If FOXA2 marks cardiac mesoderm in humans as reported in mice, 9 , 23 it may be too transient to detect in our system. Download figure Open in new tab Figure 1. Multi-lineage cell diversity in 3D can be achieved by modulating the duration of Activin and BMP signaling. A) Schematic of hESC 3D EB differentiation protocol with pathway inhibitors introduced at day 2 to promote cellular diversity. All factors were removed at day 3. B) Western blots of receptor-SMADs 6 hours after small molecule addition following the protocol in A. C) RT-qPCR analysis of marker gene expression in a differentiation time-course following the protocol in A. D) Live imaging of SOX17-Tom and NKX2-5-GFP at 3 timepoints of differentiation. E) Representative flow cytometric analyses of surface markers PDGFRα and CD235, together with SOX17-Tom and NKX2-5-GFP at 2 timepoints of differentiation with the indicated concentrations of SB and DMH1. F) Positive cell proportions for day 4 and G) for day 7.5. H) Immunostaining of whole mount EBs stained for FOXA2 (red) at day 4, and α-actinin (red) and WT1 (green) at day 12. I) RNA-seq quantifications of NKX2-5 and FOXA2 after the sorting of 8 populations at day 7.5 based on SOX17-Tom and NKX2-5-GFP. The boxplots follow standard Tukey representations and are colored by the lineage markers. Typical population sorting gates are indicated in E. J) Hierarchical clustering of cell type and lineage markers based on RNA-seq data from the 8 populations. Expression levels are represented as Z-score normalized log2 FC differences. Data are represented as mean ± SD (n=2 independent biological experiments for C, F and G, and n=3 for the RNA-seq data of I and J). Scale bars represent 300 μm. A.U., arbitrary units; EB, embryoid body; Exp., experiment; FC, fold change; SD, standard deviation; STM, septum transversum. See also Figure S1 . Establishing anterior-posterior primitive streak-like patterning allowed us to select cardiac-permissive concentrations of SB and/or DMH1 (both 0.6 μM). At these concentrations we noted that WT1+ mesothelial / epicardial cells developed in control conditions and following the addition of SB, but not following DMH1 treatment, whereas cardiomyocytes formed in all conditions ( Figure 1H ). This led us to question what molecular mechanisms limit cardiac differentiation under suboptimal (high BMP or high Activin signaling) conditions. And what markers discriminate cardiac progenitors with different developmental potentials, including those with or without the ability to make epicardial cells? Based on NKX2-5-GFP and SOX17-Tom expression at day 7.5, a timepoint when cardiac progenitors may still not be fully differentiated, we sorted 8 populations for RNA-seq which we hypothesized would be developmentally distinct including cardiac and closely related lineages and thus important to characterize ( Figures 1E and 1I ). Developmental marker analysis and hierarchical clustering showed that the GFP+ cells had a uniformly myocardial identity except when co-expressed with Tomato, which had an endoderm signature specific of thymic epithelial cells ( Figure 1J ). 24 Many ‘cardiac progenitor’ TFs including GATA4/5/6 , TBX5 , and MEIS1/2 were expressed at similar levels across control and SB-treated cells independently of GFP expression. However, whereas GFP+ cells in these conditions had a myocardial-commitment signature of high SRF , MEF2C , and MYOCD expression, the GFP-neg cells had higher HAND1 and other markers described to demarcate the JCF, proepicardium, septum transversum mesenchyme, and extraembryonic mesoderm 11 , 12 . In contrast, with DMH1, emerging GFP+ Tom-neg cells displayed an anterior SHF-like profile with high FGF8 and FGF10 expression and an absence of posterior HOXB1 . Expression values for some key markers are shown in Figure S1F . HAND1 and HOX transcription factors regulate the cardiac inhibitory effect of high BMP signaling We next performed chromatin accessibility analysis on the isolated populations. To obtain epigenetic signatures, we collated all regions which were unique to population classes based on NKX2-5-GFP and SOX17-Tom status ( Figure S2A ). Motif enrichment analysis provided insight into the TFs controlling the identify of each population ( Figure 2A ). Recognition motifs of GATA, CTCF, MEIS/PBX, MEF2, NKX2-5, and TBX20 were significantly enriched in the cardiac GFP+ population. In contrast, motifs of TEAD, WT1, HOX, and HAND were enriched in the GFP-neg Tom-neg population. This was consistent with the expression of these TFs and their local chromatin state ( Figures S2B–D ). We attribute the HAND motif enrichment to the expression of HAND1 , and the HOX motif enrichment primarily to the expression of the HOXB cluster, principally HOXB1–6 ( Figures S2B–D ). The WT1 motif enrichment may be explained by the corresponding presence of epicardial cells, but there is also evidence for WT1 playing a role in cardiac progenitors. 25 , 26 The SOX17+ population was dominated by the enrichment of motifs and associated co-expression of endoderm TFs FOXA2, SOX17, and HNF4A ( Figures 2A , S2B, and S2E ). The NKX2-5+ SOX17+ population additionally showed a significant enrichment of the motif for NKX2-5, confirming its functional role in these cells ( Figure 2A ). Download figure Open in new tab Figure 2. HAND1 and HOX transcription factors regulate the cardiac inhibitory effect of high BMP signaling. A) Transcription factor motif enrichment analysis of ATAC-seq data from sorted cells at day 7.5. The results are from an input of all regions of open chromatin uniquely identified in each of the 4 population classes as indicated by the colored boxes and based only on SOX17-Tom and NKX2-5-GFP status. For each motif, the average background value is indicated by the dotted line. Significance above background for any given population *p<1E-6. B) Typical EBs at day 10 under control or SB-treated (day 2–3) conditions, from wild-type, YAP1-null, WT1-null, HOXB1–3-null, and HAND1-null hESCs. The brightfield image is overlayed by the NXK2-5-GFP signal, with quantification by flow cytometry shown in C). Each point displays the result of an independent biological experiment. Data are represented as mean ± SD (n=3–7). Significance between any population was assessed by One-Way ANOVA with Tukey’s multiple comparison test *p<0.05. D) Typical flow cytometric analysis of GFP fluorescence in wild-type and HAND1-null EBs at day 10 with and without SB treatment. E) Immunostaining of whole mount day 12 EBs stained for α-actinin (red) and WT1 (green). Scale bars represent 300 μm. Con, control/vehicle; EB, embryoid bodies. See also Figure S2 and Figure S3 . Considering these data, we hypothesized that some of the TFs active in the NKX2-5-neg mesodermal cells in the high BMP signaling environment, i.e. HAND1, TEAD, HOX, and WT1, might have an important role in the specification of JCF or extraembryonic mesodermal progenitors, and so potentially determine cell fate in this high BMP environment and direct cells away from becoming cardiomyocytes (NXK2-5-GFP+). To test this, we made a series of gene knockout hESC lines ( Figures S3A–F ). To overcome the genetic redundancy of TEAD and HOX, we created YAP1-null hESCs to inhibit the TEAD-YAP1 complex, and generated a triple knockout of HOXB1-2-3 (HOXB1–3-null) to restrict HOX activity ( Figures S3A–C ). The impact of each TF type on cardiac differentiation and NKX2-5-GFP expression was assessed ( Figures 2B-E ). Although YAP1-null cells could express PDGFRα at day 4 ( Figure S3D ), they then ceased to proliferate and failed to express any GFP even under control conditions. In contrast, the WT1-null, HOXB1–3-null and HAND1-null hESCs could all differentiate to cardiomyocytes under control conditions, with GFP levels equivalent to the wild-type ( Figures 2B and 2C ). The WT1-null behaved like the wild-type in terms of cardiac repression by SB, whereas the HOXB1–3-null and the HAND1-null showed significant resistance to SB, which in the case of the HAND1-null was essentially complete ( Figure 2D ). WT1+ epicardial cells still formed from the HOXB1–3-null (data not shown) but were completely absent in the HAND1-null where the EBs instead contained an increased number of cardiomyocytes consistent with the increased NKX2-5-GFP expression ( Figure 2E ). These data position HAND1 and the HOXB family as fate determining factors in a high BMP signaling environment and suggest that HAND1 may play a major role in the development of JCF-like multipotent progenitors. At higher levels, HAND1 and the HOXB family may block cardiac fate and potentially drive an extraembryonic fate. Single-cell RNA-sequencing supports a functional role for HAND1 in the development of a multipotent ‘JCF’-like population To gain insight into the role of HAND1 in cardiac and non-cardiac progenitors, we characterized the cellular identities in the wild-type and HAND1-null cell models during differentiation by single-cell RNA-sequencing. To promote the development of a broad range of cell types, we combined populations of vehicle, SB-treated and DMH1-treated cells (as described above) at a ratio of 0.5:0.25:0.25 and collected cells at 7 timepoints. Using UMAP with mutual nearest neighbor (MNN) correction, unsupervised clustering and detailed marker analysis on the combined transcriptomes of >100k wild-type and HAND1-null cells, we could annotate 8 broad population types: endoderm (END), mesoderm (MES), mesodermal progenitors (MP), cardiac progenitors (CP), endothelial (EC), cardiomyocytes (CM), epicardial (EPI), and fibroblast-like cells (FLC) and calculated assignment frequencies by day for wild-type ( Figures 3A and 3A ’ ) and for HAND1-null ( Figure 3B and 3B ’ ). Figure 3C shows the levels of prominent markers of each cell population in the wild-type. Whereas the populations in the upper part of the MNN-UMAP plots are largely overlapping ( Figures 3A and 3B ), there was a clear absence of epicardial and fibroblast-like cell development in the HAND1-null. Thus, HAND1 promotes epicardial and fibroblast differentiation, consistent with the previous data ( Figure 2E ). Download figure Open in new tab Figure 3. Single-cell RNA-sequencing shows that HAND1 is important for epicardial and fibroblast-like cell development. MNN-UMAP plots of A) wild-type and B) HAND1-null cells from day 3–10 of differentiation colored by cell type. The dashed red lasso highlights the loss of epicardial and fibroblast-like cells in the HAND1-null. A’) The frequency of each cell type by day in wild-type and B’) HAND1-null. C) The levels of prominent markers identifying each population in wild-type cells. To gain further biological insight into cardiac cell lineage development, we performed line-specific analyses using only hyper variable genes (HVGs) implicated in cardiac development (∼1000 genes), while retaining the previous cell type annotations. Monocle 3 was then used to assign pseudotime and lineage trajectory data from root mesoderm. 27 In the wild-type, the major mesoderm derivatives still partitioned into four terminal lineages ( Figure 4A ). Lineage trajectory analysis enabled the inference of their developmental origins ( Figure 4A ’ ). From differential gene expression analysis of the two lineages, we deduced that the early bifurcation point in the Monocle trajectory may represent an FHF- and SHF-like split in the mesoderm ( Figure S4A and A’ ). We noted HAND1 , TBX2 and HCN4 expression on the FHF-like lineage ( Figures 4B ). TBX2 lineage tracing showed it marks cells which make a similar contribution to the left ventricle as HAND1 , 28 and HCN4 is recognized as an FHF progenitor marker. 29 The SHF-like lineage showed prominent ISL1 , FGF10, FOXC1/2 and JAG1 expression ( Figure 4B ), which are well-established SHF markers. 8 , 30 – 32 Whereas the SHF-like lineage appeared to progress solely to cardiomyocytes, the FHF-like lineage – which at its early stage included a strong JCF-like signature ( HAND1 , MAB21L2 , BNC2 and MSX1/2 ) 8 , 11 , 12 – displayed multilineage fate also giving rise to epicardial and fibroblast-like cells. Additionally, we detected a side branch from the FHF/JCF-like progenitors with high transient CDX2, GATA2 , and TBX3 expression and low TBX5 , which may represent an extraembryonic mesoderm-like lineage ( Figures 4A ’ and 4B ). 33 Download figure Open in new tab Figure 4. Single-cell RNA-sequencing supports a functional role for HAND1 in the development of a multipotent ‘JCF’-like population. A-D) UMAP plots from Monocle3: A) colored by cell type in the wild-type and C) in the HAND1-null line; A’) colored by pseudotime with predicted trajectories overlayed, in the wild-type and C’) in the HAND1-null line. Endpoint cell types and FHF/JCF- and SHF-like populations have been annotated and their marker frequencies shown. B) Gene expression of lineage and cell type markers for the wild-type and D) HAND1-null lines. E) SCENIC GRN analyses of cardiac progenitors extracted from the two lineage trajectories of the wild-type differentiation and the single lineage of the HAND1-null trajectory. F) Activity of selected regulons enriched in cardiac progenitor lineages mapped to UMAP plots of wild-type and HAND1-null cells. White indicates the regulon activity is below threshold. The number of gene targets in each regulon is indicated following the name of the TF. G) Heatmap of relative regulon activity by cell line, cell type, and subpopulation, grouped by lineage and stage of enrichment. The cardiac progenitor data represent cells extracted from each lineage. *Identified in HAND1-null population. # Identified in CPC lineage analysis. DEG, differentially expressed genes; FHF, first heart field; JCF, juxta-cardiac field; SHF, second heart field; GRN, gene regulatory network; SCENIC, single-cell regulatory network interference and clustering. See also Figures S4 and S5 . In contrast to the wild-type, the HAND1-null showed a single dominant trajectory to cardiomyocytes ( Figures 4C and 4C ’ ). A small number of WT1+ (EPI) and DCN+ (FLC) cells clustered together, possibly explained by partial compensation from HAND2 which was ectopically activated in 10% of the HAND1-null mesoderm, but the development of FHF/JCF-like and extraembryonic-like progenitors and epicardial and fibroblast-like cell derivatives, was almost entirely lost. The main lineage of progenitors in the HAND1-null had a SHF-like signature with an equivalent or higher frequency of FGF10 + (51% vs 46%) ISL1 + (43% vs 43%), and TBX1 + (18% vs 10%) cells compared to the wild-type SHF-like lineage ( Figures 4A ’ and 4C’ ). Following the assumption that the main differences between the wild-type and HAND1-null differentiations can be explained by a change in cardiac lineage selection and, therefore, by the gene regulatory networks (GRNs) coordinating these lineages, we performed GRN analysis using SCENIC to identify active TFs and their target genes (regulons). 34 , 35 To maximize regulon detection power but also ensure robustness, we performed analyses on 1) all wild-type cells 2) wild-type CP lineages-only (shown in Figure S4A ) and 3) all HAND1-null cells, and ran the algorithm multiple times retaining only regulons identified in every run. A total of 307 regulons met these criteria ( Figure S4B ). Using our lineage analysis, we first focused on regulons with FHF/JCF- or SHF-bias in activity and examined these in the wild-type and HAND1-null populations ( Figures 4F and 4G ). The SHF-like lineage was characterized by increased activity of FOXC1/2, SIX1, PRRX2 and MEIS1 regulons, and regulons controlled by regulators of cell proliferation including E2F family members ( Figure 4F and 4G ). In contrast, the early FHF/JCF-like lineage was characterized by increased activity of HOX, JDP2, MSX1/2, FOXF1 and LEF1 regulons and many other regulons which were later enriched in the epicardial or fibroblast-like cells such as CREB3L1 ( Figure 4F and 4G ). The HAND1-null progenitor lineage most resembled the SHF-like lineage of the wild-type and had lower activity of the FHF/JCF-associated regulons ( Figure 4F and 4G ). A HAND1 regulon was identified in the wild-type cardiac progenitor analysis ( Figure 4G ) but not in analyses of the full datasets where the regulon was lost due to lack of motif enrichment detected using cisTarget. 34 MESP1, another bHLH TF with an established role in cardiac mesoderm, was similarly lost, suggesting an insensitivity of the method to correctly prune / retain the regulons of some bHLH TFs using the standard gene regulatory space. Our SCENIC analysis allowed us to rank the most specific regulons for each cell type ( Figures S5 ). Comparing the cardiomyocytes in the wild-type to the HAND1-null revealed that the dominant regulons were controlled largely by the same factors; this included classical cardiac TFs MEF2, SRF, GATA4, and TBX5; and metabolic regulators including estrogen-related receptors (ESRRA and ESRRG) and the transcriptional co-activator PGC-1α (PPARGC1A). These data support the conclusion that the major function of HAND1 in this system is not in determining the gross identity of cardiomyocytes but rather in controlling the type of mesodermal progenitors that are made. HAND1 controls chromatin and enhancer landscapes to pattern mesoderm and program cell fate To investigate how HAND1 functions in mesodermal progenitors, we first epitope-tagged endogenous HAND1 in wild-type cells to map its DNA binding by ChIP-seq ( Figure S6A ). HAND1 could be detected already by day 3 and was increased by SB treatment ( Figure 5A ). At day 3, we identified 26,002 and 34,808 HAND1 binding sites in control and SB conditions, respectively, largely in distal intergenic and intronic locations ( Figure S6B ). HAND1-bound regions were enriched for the previously identified long non-classical E-box HAND1 motif, but also for the recognition motifs of other factors known to be important in mesoderm including GATA6, EOMES, ZIC3 and TBX6 ( Figure 5B ). Next, to investigate how HAND1 activity influences chromatin structure and gene expression we employed a tightly controlled system in which HAND1 could be transiently expressed at the mesoderm stage of HAND1-null cells. To do this we introduced a doxycycline (dox)-inducible HAND1-BFP transgene via lentivirus into ∼30% of the cells. Dox was added at day 2.5 of differentiation, and 12 hours later we isolated HAND1-BFP+ and BFP-negative cells for ATAC-seq and RNA-seq ( Figure 5C ). HAND1 induced a dramatic restructuring of the chromatin, largely in distal intergenic and intronic locations, with ∼17k regions displaying significantly increased accessibility and ∼9k regions displaying significantly decreased accessibility ( Figures 5C ’ and S6B ). A total of 1481 genes were significantly upregulated, and 1525 genes were significantly downregulated by the HAND1 transgene. Integrating the ATAC-seq and ChIP-seq data revealed that HAND1 binding was enriched in the regions of increased accessibility ( Figure 5C ’ ). Examining these genomic subsets allowed us to gain further insight into the TFs impacting the chromatin environment as a result of HAND1 expression. This showed that the above HAND1 motif was uniquely enriched in the regions of HAND1 binding with increased accessibility and not in the regions with decreased accessibility even where HAND1 was bound ( Figure S6C ). HAND1 binding and chromatin opening was also specifically associated with the binding motifs of SMAD, FOXO1, FOXF1, and PBX/MEIS. Conversely, regions of decreased accessibility displayed an exclusive enrichment of the classical E-box motifs CAGCTG and CACCTG, particularly in regions bound by HAND1. This suggests that the binding partners of HAND1 may determine its occupancy at different motifs. Download figure Open in new tab Figure 5. HAND1 controls chromatin and enhancer landscapes to pattern mesoderm and program cell fate. A) Western blot of HAND1-AM-Tag (Active Motif) at day 3 and 4 of differentiation in different conditions compared to undifferentiated hESCs. B) Motif enrichment analysis of HAND1 ChIP-seq peaks with control (low HAND1) and SB (high HAND1) sets merged. Significance is relative to the background control. C) Doxycycline-inducible HAND1-BFP transgene expressed in HAND1-null cells for cell sorting by HAND1 expression (12-hour induction) and molecular analysis by ATAC- and RNA-seq (both performed on 3 biological replicates). C’) Euler plots for the ATAC-seq analysis representing the number of differentially accessible chromatin regions in HAND1+ vs HAND1-negative cells and the subset with detected HAND1 binding. D) Schematic illustration of the motifs enriched in chromatin with increased or decreased accessibility as a result of HAND1 expression. E) Identification of HAND1[+] (activating) and HAND1[-] (repressing) regulons. F) Chord plot showing significantly enriched gene ontology terms for HAND1 target genes. The red and blue text coloring indicates terms enriched in activated targets and repressed targets respectively. G) UMAP plots showing the activity of constant and HAND1-sensitive regulons in endoderm and mesoderm of wild-type and HAND1-null cells. The activity of the HAND1[+] and HAND1[-] regulons are shown with the domains of strong activity highlighted. H–J) ATAC-seq in HAND1-null and HAND1+ (12 h induction) populations, and HAND1 binding at low-HAND1 (CTRL) and high-HAND1 (SB) levels at the H) FOXF1 locus, I) HOXB cluster and J) CDX2 locus. K) Gene regulatory network model of HAND1. * Direct early target of HAND1. The dashed line indicates a link at high HAND1 levels. EB, embryoid bodies; CTRL, control / vehicle-only; DEG, differentially expressed genes; EMT, epithelial-mesenchymal transition; ExEm, extraembryonic mesoderm; GO, gene ontology; GRN, gene regulatory network; JCF, juxta-cardiac field; WT, wild-type. See also Figure S6 . Regions of decreased accessibility were also enriched in motifs for GSC and SOX17, consistent with the down regulation of these genes ( Figure S6C ). Supportive findings were observed in the analysis of published HAND1 ChIP-seq data from a later time-point in mesoderm differentiation (data not shown). 36 Figure 5D shows a schematic representation of these changes. Next, we aimed to understand the HAND1 GRN in mesodermal progenitors. First, we identified high confidence target genes (regulons) using the strategy illustrated in Figure 5E . The regulons comprised 335 activated targets and 350 repressed targets. 51% of the activated and 58% of the repressed genes were concordantly changed by the 12 hours of HAND1 transgene expression, suggesting those as probable early targets. Gene Ontology analysis provided an overview of the biological processes impacted by HAND1 ( Figure 5F ). Activated targets were enriched for TFs and coactivators (including MESP1/2 , FOXF1 and TEAD1 ), genes involved in different aspects of heart morphogenesis, anterior-posterior pattern formation (including HOXB2–9 ) and regulators of cell migration and epithelial-mesenchymal transition (EMT) (including SNAI1/2 and CRB2 ). Conversely, genes potentially repressed by HAND1 were associated with processes including neural development, WNT signaling and endoderm development (including SOX17 and NODAL ). We used the gene-set analysis tool AUCell to quantify the activity of the activating HAND1[+] (335g) and repressive HAND1[-] (350g) regulons in our single cell data and compared the results to other relevant regulons ( Figure 5G ). We noted that FOXF1 and HOXB6 were early positive targets of HAND1, and the activity of their regulons overlapped, whereas the regulon of SOX17 overlapped with the HAND1[-] regulon. These patterns were disrupted in the HAND1-null, while the patterns of endoderm-restricted FOXA2 and HNF4A activity were unchanged. To overcome its absence, a regulon controlled by the cardiac mesoderm master regulator MESP1 was defined using our co-expression data from SCENIC and published RNA-seq data from Mesp1 overexpression. 37 HAND1 was confirmed as a MESP1 target in our system. The activity of the MESP1 regulon was not significantly changed by HAND1, suggesting that MESP1 mostly acts upstream and independently of HAND1. HAND1 binding and increased chromatin accessibility was evident at the FOXF1 locus and HOXB cluster ( Figures 5H and I ). Accessibility at the FOXF1 promoter was robustly increased later, as seen in the high-HAND1 populations in our day 7.5 ATAC-seq data ( Figure S2F ). A total of 48 genes were common targets of HAND1, FOXF1, and HOXB6, including the genes encoding these upstream regulators as well as other HOXB TFs, suggesting positive feedback and cross-activation. Common targets also included the TFs LEF1 , ID2 , and TEAD1 . Targets specific to FOXF1 and HOXB6 included JDP2 , MSX1/2 , PRDM6 , and PRRX1 , the activities of which were all increased in the FHF/JCF-like lineage, fibroblasts or epicardial cells ( Figures 4G and S4B ). In line with the potential of HAND1 to promote extraembryonic mesoderm specification, we identified HAND1 binding at the CDX2 locus, but only in the high-HAND1 condition ( Figure 5J ). CDX2 was also a predicted target of FOXF1. The above data allowed us to construct a GRN model of how HAND1 could regulate FHF/JCF and extraembryonic lineage assignment ( Figure 5K ). HAND1-expressing progenitors are expandable and unique in their multipotency for cardiomyocytes and epicardial cells Markers indicating the developmental potential of cardiac progenitors derived from hPSCs are badly needed to reduce heterogeneity in disease modelling and to develop cell therapies. Because the level of HAND1 expression was seemingly indicative of developmental lineage and fate propensity, as illustrated in Figure 6A , we addressed whether it could be used to distinguish progenitor types and identify and isolate multipotent JCF-like progenitors. To test this, we introduced a HAND1-T2A-Tomato knock-in reporter into the NKX2-5-GFP hESC line to create a double reporter ( Figure S7A ). HAND1-Tom expression became evident by day 4 of differentiation and was increased by SB and decreased by DMH1 treatment, in agreement with the level of HAND1 ( Figures 6B and 5A ). Cardiac progenitors are a transient population in this model and differentiate after little self-renewal. Aiming to circumvent uncontrolled differentiation we examined our GRN data for correlated regulators of progenitor self-renewal and identified MYC- and N-MYC-controlled regulons in mesoderm and early cardiac progenitors ( Figure 6C ). The activity of these two regulons declined simultaneously with differentiation towards cardiomyocytes, which displayed an increased activity of PGC-1α ( PPARGC1A ) and SRF ( Figure 6D ). Both MYC and N-MYC have been shown to play a role in cardiac progenitor expansion in development and in vitro. 22 , 38 – 40 A dox-inducible MYC transgene activated from day 4.75 was able to halt differentiation and maintain the activity of the MYC and MYCN regulons ( Figure 6E ). This enabled the clonal expansion of these cells supported by IGF-1 and a Hedgehog pathway activator. BMP4, FGF8 and the WNT pathway activators CHIR and WNT3A all promoted HAND1-Tom expression ( Figures S7B–F) . FGF and BMP signals also acted cooperatively to promote the induction of NKX2-5 expression, whereas WNT signals did not. To maintain progenitors with HAND1-Tom-high, Tom-low, and Tom-neg statuses, generated in SB-, vehicle- and DMH1-treated EBs respectively, we optimized three maintenance conditions. Download figure Open in new tab Figure 6. HAND1-expressing progenitors are expandable and unique in their multipotency for cardiomyocytes and epicardial cells. A) Schematic illustration of HAND1 level with primitive streak-like patterning of hESCs into cardiac and extraembryonic mesoderm. B) Live images (brightfield and fluorescence) showing HAND1-Tomato in EBs at day 5 of differentiation. C) UMAP plots showing the target genes AUC score (regulon activity) for MYC and N-MYC with differentiation in wild-type cells. White indicates the score is below threshold. D) Target genes AUC score for MYC, N-MYC, PGC-1α (PPARGC1A) and SRF through pseudotime of mesoderm differentiation. E) Scaled target gene AUC scores in bulk populations by differentiation day and impact of induction of a dox-inducible MYC transgene at day 5–6. F) Expansion of HAND1-high, HAND1-low, and HAND1-neg progenitors from SB, control and DMH1-treated conditions respectively. HAND1-high progenitors were maintained in FGF8, BMP4, and CHIR; HAND1-low progenitors were maintained in FGF8, BMP4, and WNT3A, and HAND1-negative progenitors were maintained in FGF8-only (see Methods). Flow cytometric analysis of HAND1-Tom and NKX2-5-GFP after 6 days of expansion. G) Differentiation of progenitors to epicardial and endothelial cells assessed by immunostaining for WT1 and CD31, respectively and additionally by the flow cytometric analysis of CD31. H) Differentiation of progenitors to cardiomyocytes assessed by immunostaining for cTroponin T. I) Quantification of differentiated cell types by population. Data in I represent mean ± SD, n = 3 independent biological experiments. Scale bars represent 150 μm. AUC, area under the curve. EB, embryoid body; Epi, epicardial cells; Endo, endothelial cells; CM, cardiomyocyte. See also Figure S7 . Following expansion for 6 days, fluorescence was measured ( Figure 6F ) . HAND1-Tom level correlated with a resistance to NKX2-5-GFP expression, which emerged largely from HAND1-Tom-low or Tom-neg cells. The expanded populations were further differentiated by removing dox and providing specification signals. Differentiation to endothelial and WT1+ mesothelial / epicardial cells was promoted by low-level FGF plus VEGF ( Figure 6G ), whereas cardiomyocyte differentiation was promoted by high-level FGF plus TGFβ ( Figure 6H ). All populations could generate CD31+ endothelial cells, whereas WT1+ cells were generated exclusively from HAND1-Tom-high and Tom-low cells. HAND1-Tom-high progenitors were unable to make cardiomyocytes, whereas the other populations did so efficiently. HAND1-Tom-low progenitors were multipotent and able to make all three cell types with good efficiency suggesting that they represent a JCF-like population ( Figure 6I ). Discussion Through a systematic study of cell diversity in a human model of early cardiac development, we have deciphered lineage trajectories and identified GRNs controlling cell identity and fate. Focusing on the decision point in mesoderm when progenitors may follow either a cardiac FHF/JCF, SHF, or an extraembryonic path, we have uncovered a fate-determining role for the bHLH factor HAND1. The recent seminal discovery of the JCF, a progenitor population marked by Hand1 , which contributes to both cardiac and extraembryonic tissues including to the epicardial lineage, underpins the developmental significance of these results. 11 , 12 We found that HAND1 is a master regulator of JCF progenitor specification in mesoderm; in our 3D differentiation model its knockout caused a switch in fate towards a SHF-like lineage. HAND1 restructures chromatin and initiates GRNs including those controlled by HOX and FOXF1. At high levels HAND1 drives an extraembryonic differentiation trajectory. Importantly, HAND1-low progenitors can be expanded and demonstrate a unique multipotency to give rise to both cardiomyocytes and epicardial cells. In line with previous research supporting the importance of signal duration in early development, 41 , 42 we confirmed that the duration of BMP and Activin signaling was a major factor in programming and patterning endoderm and mesoderm in differentiation. We further queried why cardiomyocyte development is inhibited or delayed in mesodermal progenitors exposed to a high BMP signaling environment, which programs progenitors into making more epicardial and fibroblast-like cells. The close proximity of the FHF/JCF and the extraembryonic mesoderm, and overlap in their molecular profiles, makes this question non-trivial. 11 In our system, the expression of GATA4 , TBX5 , MEIS1 , and PDGFRA failed to distinguish closely related progenitors fated to, or not to, become cardiomyocytes. However, by analyzing differentially accessible chromatin we deduced that, whereas GATA was universally important, the TFs MEIS, MEF2, NKX2-5, TBX20, and CTCF, presumably acting with some cooperativity, were distinctive in shaping the chromatin structure of cells becoming cardiomyocytes. Conversely, the TFs TEAD, HOX, WT1, and HAND1 were predicted to play the corresponding role in cells not progressing to cardiomyocytes. Focusing on the latter, the function of the Hippo pathway TEAD-YAP1 was found to be important in cardiac progenitor proliferation and survival, consistent with the pathway’s known importance in controlling proliferation and organ size in development and survival in PSC differentiation. 43 , 44 WT1 knockout did not affect cardiomyocyte differentiation, or the inhibitory effect of high BMP signaling levels, thus not supporting a major function for WT1 in human cardiac progenitor formation or differentiation as reported in mouse and zebrafish. 25 , 26 The enrichment of motifs for HOX TFs was unexpected as their role in cardiac progenitors has been little studied. However, Hoxb1 is required to maintain proliferation in the posterior SHF in mice 45 and was observed in the JCF. 11 The HOXB1–3-null showed that HOX TFs may play a significant role in restricting cardiac differentiation in a high BMP signaling environment. Finally, HAND1 was found to be an activator of HOX and the major factor in interpreting the response to increased BMP signaling, and promoting epicardial and fibroblast differentiation over cardiomyocytes. Single cell RNA-seq showed that, in response to modulation of only BMP and Activin signal duration, wild-type cells followed multiple developmental trajectories into FHF/JCF-, extraembryonic mesoderm-, and SHF-like lineages. Interestingly, the fibroblast-like cells, which were generated in large number in the high BMP condition, did not come from the epicardial cells but instead directly from earlier progenitors; some from HAND1 + TBX5 -neg CDX2 -neg cells and others from HAND1+ CDX2+ cells which may represent an extraembryonic mesoderm-like population. Despite a degree of compensation by HAND2, as observed in Hand1-null mice, 46 the loss of HAND1 resulted in the cells preferentially committing to a SHF-like lineage, suggesting that HAND1 assigns cells to the other lineages. The validity of these trajectories was supported by our comprehensive GRN analysis. The recent seminal discovery of a shared cardiac-extraembryonic-progenitor marked by Hand1 has raised the question of what function this bHLH TF has in the formation or differentiation capacity of multipotent progenitors. 11 , 12 Understanding the developmental functions of Hand1 has been complicated by its important early role in placental development. 13 , 15 However, using conditional mutants, a function in cardiac cell proliferation has been described, 14 , 16 , 17 and it has been linked to congenital heart disease. 18 Using a hESC model, free of complications related to extraembryonic dependence, we found that HAND1 could induce a rapid change in chromatin accessibility in mesoderm. Working with putative binding partners including GATA6, ZIC2/3, FOXO1, FOXF1, PBX/MEIS, and SMAD, our results suggest it plays a major role in activating GRNs defining the FHF/JCF and extraembryonic mesoderm, while inhibiting endoderm differentiation. HOXB3 , HOXB6 and FOXF1 were found to be early HAND1 target genes and inferred to be central components of the GRN directed by HAND1, along with other TFs including MSX1. More work is needed to fully understand the functions of these downstream factors. In our partner manuscript, we show that HOX can inhibit cardiomyocyte differentiation by partnering with MEIS1 and preventing it from binding with GATA to drive differentiation. 47 At high levels, HAND1 bound and activated CDX2 expression. CDX2 is known to promote the development of extraembryonic mesoderm from hPSCs and inhibit cardiac differentiation, as we observed in this condition. 33 , 48 Positioning HAND1 in the developmental hierarchy, our data confirmed that HAND1 lies immediately downstream of the bHLH TF MESP1, a master regulator of cardiac mesoderm development. 37 They reinforce each other’s expression and have many shared target genes including SNAI2 to promote EMT, and the ability to repress SOX17 and endoderm development. However, according to the MESP1 ChIP-seq data of Liang et al., they appear to share a limited number of binding sites (14% at low HAND1 and 20% at high HAND1 levels) and their co-binding does not correlate with HAND1-mediated chromatin activation. 49 This suggests that MESP1 functions more widely to program mesoderm and HAND1 specifically patterns it to promote FHF/JCF and extraembryonic fates. The additional promotion of EMT and migration provided by HAND1 is consistent with JCF progenitors being the most highly migratory. 8 Finally, with regenerative medicine in mind, and building on our previous work, 22 , 50 we took advantage of the ability of MYC and N-MYC to promote cell self-renewal in cardiac development to study progenitor fate as a function of HAND1 level. Each progenitor type required a carefully tuned signaling environment to maintain its identity. WNT3A, FGF8, and BMP4 could support the self-renewal of HAND1-low progenitors, which proved to be multipotent for cardiomyocytes and epicardial cells. This provides further evidence of shared lineage at this stage of development. With the added potential to make endothelial cells, these progenitors are of significant interest for cardiac regeneration where all three cell types have been shown important for tissue repair. 51 – 53 In summary, this study provides a comprehensive map of human cardiovascular development and reveals important regulatory mechanisms, most notably a lineage-determining function for HAND1. Controlling these events is crucial for achieving accurate differentiation and maximizing the opportunities of hPSCs for translational applications. These findings will particularly support the modelling of congenital heart disease and the development of cell-based therapies. Limitations of the study This study has not extended to an in-depth analysis of the putative partners of HAND1 including FOXF1, SMADs, and PXB/MEIS, to know if they function together on chromatin and whether partner choice impacts the developmental program. Another limitation is the measurable compensatory activity of intact HAND2 in HAND1-null mesoderm and whether this influences the phenotype. The generation of a HAND1/2-null line will be required to investigate this. The inference of regulons in our SCENIC analysis was limited to the default regulatory sequences, but represents an important foundation on which to build. Author contributions M.B. conceived the study; A.L., N.P., M.D., A.A., N.H., and M.B. designed and conducted the experiments; A.L., N.P., I.L., and M.B. analyzed data; Z.D., J.W. and N.B., supported the data analysis; M.B. and A.L. wrote the paper. Declaration of interests The authors declare no competing interests. STAR Methods Maintenance of human ESC lines HES3 NKX2-5 eGFP/w hESCs 54 and gene-edited derivatives were maintained on a layer of mitotically-inactivated mouse embryonic fibroblasts (MEFs) in DMEM/F12-based medium containing 20% (v/v) KnockOut™ Serum Replacement (Gibco), 100 mM non-essential amino acids (Gibco), 2 mM GlutaMAX™ (Gibco), 0.1 mM β-mercaptoethanol (Gibco) and 10 ng/ml bFGF (Miltenyi Biotech). hPSCs were passaged using TrypLE Select (Thermo Fisher Scientific) every 3-4 days. Gene editing by CRISPR-Cas9 Fluorescent reporters for SOX17 (dTomato), FOXA2 (mTagBFP2) and HAND1 (dTomato) were generated using a 3’ gene targeting strategy to incorporate the fluorophore downstream of a T2A peptide. C-terminal epitope-tagged HAND1 was generated using a targeting strategy to incorporate an AM-Tag (Active Motif). Targeting plasmids incorporated 700-800 bp homology arms and a loxP-flanked neomycin resistance cassette and were assembled using NEBuilder HiFi Assembly (NEB). The plasmids have been deposited as follows: pSOX17-T2A-dTom-PGK-Neo (Addgene ID 223190), pFOXA2-T2A-mTagBFP2-PGK-Neo (Addgene ID 223191), pHAND1-T2A-dTom-PGK-Neo (Addgene ID 223190), pHAND1-AM-Tag-PGK-Neo (Addgene ID 223216). We recently described the generation of a HAND1-T2A-Tomato reporter using the same method. 55 In brief, ribonucleoprotein (RNP) complexes of CRISPR-Cas9 were made by combining 120 pmol of crRNA:tracrRNA mix (IDT) with 67 pmol of Cas9 (Alt-R™ S.p. Cas9 Nuclease V3, IDT), final concentrations 1.3 and 1.1 µM respectively. gRNA sequences are listed in Table S1 . 0.5x10 6 cells were transfected with the RNP mix, 120 pmol Alt-R® Cas9 Electroporation Enhancer (1.3 µM final) and 1 µg of donor plasmid (pSOX17-T2A-dTom-PGK-Neo, pFOXA2-T2A-mTagBFP2-PGK-Neo pHAND1-T2A-dTom-PGK-Neo or pHAND1-T2A-Am-Tag). Transfection was performed by electroporation using an Amaxa Nucleofector II combined with Human Stem Cell Nucleofector Solution 2 (Lonza, ref #: VPH-5022) and programme B-16. Cells were selected with 100 μg/ml neomycin for 5 days and then electroporated with 1 µg of pCAG_CRE_IRES_GFP plasmid (Addgene 26646) expressing Cre recombinase and GFP to excise the neomycin selection cassette. 48 hours later, single GFP+ cells were sorted by FACS (BD Influx) into MEF-coated 96-well plates. DNA was isolated from clones using a Qiagen DNA Mini kit for PCR screening using primers spanning the limits of the homology arms. PCRs were performed using Q5 High-Fidelity DNA polymerase (NEB) or MyFi™ DNA Polymerase (Bioline), with the supplier’s standard protocol. Clones with correct left and right arm integration were further examined for the presence of a targeted or untargeted allele to verify zygosity. Accurate targeting was confirmed by Sanger sequencing. Gene knockouts (KO) of YAP1 , WT1 , HOXB1-2-3 and HAND1 were generated by targeting the earliest coding exon of each gene with a pair of gRNAs predicted to generate a frameshift deletion. The HOXB1-2-3 triple KO line was made by sequential mutation of HOXB1 , HOXB2 and HOXB3 . RNP complexes were generated and delivered as above. gRNA sequences are listed in Table S2 . 3 days after targeting, cells were cloned by FACS into MEF-coated 96-well plates. DNA was isolated from clones using a Qiagen DNA Mini kit for PCR screening using primers spanning the predicted deletion. Homozygous mutation was confirmed by Sanger sequencing. TO-HAND1 lentiviral vector cloning The HAND1 coding sequence (CDS) was obtained commercially (Horizon Discovery; Clone 3162118) and cloned into pENTR4-FLAG (Addgene ID 17423). A T2A-mTagBFP2 fragment was inserted downstream of the CDS to generate an entry vector containing FLAG-HAND1-T2A-mTagBFP2 (deposited as Addgene ID 223192). The entry vector was recombined with the destination vector FU-tetO-Gateway (Addgene #43914) using Gateway™ LR Clonase™ II Enzyme Mix (Thermo Fisher Scientific) and transformed into One Shot™ Stbl3™ Chemically Competent cells (Thermo Fisher Scientific). Lentivirus production and hPSC transduction The lentiviral packaging vectors (pMDLgp.pRRE [Addgene ID 12251]; pMD.G [Addgene ID 12259] and pRSV-Rev [Addgene ID 12253]), together with the transgene-containing destination vector (TO-HAND1 as described above or TO-MYC [Addgene ID 19775]), were transfected into Lenti-X 293T cells (Takara Bio) using Lipofectamine 2000 (Thermo Fisher Scientific). The supernatants were collected daily for 2 days, and the lentiviruses concentrated by ultracentrifugation and stored at -80°C. The TO-HAND1-T2A-BFP system was introduced by transducing HAND1-null hPSCs with two lentiviruses, one carrying the HAND1-T2A-mTagBFP2 transgene under the control of a TRE-CMV promoter (described above), and the second carrying the tetracycline transactivator (Addgene ID 19780). A TO-MYC system was introduced by transducing the hPSCs with two lentiviruses, one carrying a MYC T58A transgene under the control of a TRE-CMV promoter (Addgene ID 19775), and the second carrying the tetracycline transactivator (Addgene ID 19780). Polyclonal populations were used throughout. hPSC differentiation Differentiations were performed in serum-free BPEL medium 56 containing 1 μg/ml insulin using an embryoid body (EB) system. EBs were formed by depositing 2500 cells in 50 μl differentiation medium per well in V-bottomed 96-well plates (Greiner). The following growth factors were present for the first 3 days of differentiation: 25 ng/ml BMP4 (R&D Systems), 25 ng/ml Activin A (Miltenyi Biotech) and 1.5-1.75 μM CHIR99021 (Selleckchem). Where indicated, the ALK5 inhibitor SB431542 (Generon), ALK2/3 inhibitor DMH1 (Merck), or vehicle-only (DMSO), were added on day 2 at the specified concentration as a 2 μL bolus and mixed. On day 3, wells were refreshed with 100 μl BPEL. On day 6, wells were refreshed with 100 μl BPEL supplemented with 100 pg/ml bFGF (Miltenyi Biotech). On day 9, wells were refreshed with 100 μL BPEL. Where indicated, on day 7.5, cells were dissociated with 1X TrypLE Select and sorted by Tomato and GFP fluorescence using FACS (BD Aria) for RNA- and ATAC-seq. DAPI was used as a viability dye. Doxycycline-inducible HAND1 ( TO-HAND1-T2A-BFP) hESCs were differentiated in monolayers in BPEL plus growth factors, as described above. On day 2.5, 1 µg/ml doxycycline was added for 12 hours. On day 3, cells were dissociated with 1X TrypLE Select and sorted by BFP level (negative, medium and high) using FACS (BD Influx) for RNA- and ATAC-seq. 3 µM DR™ was used as a viability dye. Progenitor cell maintenance and differentiation For progenitor culture of TO-MYC cells. On day 4.75, doxycycline was added to EBs at 1 μg/ml for 24 hours. After dissociation with TrypLE Select, progenitors were seeded on a thin layer of undiluted Matrigel (14 μl/cm2) (growth factor reduced; Corning 354230) in organ culture dishes (BD Falcon) at a cell density of 8x10 3 /cm 2 in BPEL medium containing 1 μg/ml doxycycline, 5 μM SB431542, 100 ng/ml LONG R3 IGF-1 (hereafter, IGF-1) and 1 μM SAG (Selleckchem) plus additional supplements as indicated: 5-100 ng/ml FGF8b (Proteintech), 0.5-25 ng/ml BMP4 (R&D Systems), 200 ng/ml WNT3A (Proteintech) and 1 µM CHIR99021 (Selleckchem). Cultures were refreshed every 2 days. Cardiomyocyte differentiation was performed by seeding cells at a density of 2x10 5 cells/cm 2 on Matrigel-coated (1:100 dilution) 96-well plates (BD Falcon) in BPEL medium supplemented with 100 ng/ml IGF-1, 1 μM SAG, 1 ng/ml bFGF and 2.5ng/ml TGFβ1 (Peprotech). On day 2, medium was changed to BPEL supplemented with 100 ng/ml IGF-1, 100 pg/ml bFGF and 2.5ng/ml TGFβ1. On day 4, medium was changed to BPEL supplemented with 100 ng/ml IGF-1 and 10 pg/ml bFGF. On day 6, medium was changed to BPEL supplemented with 100 ng/ml IGF-1, 100 nM T3 (Merck) and 1 μM Dexamethasone (Merck). Cells were fixed for immunocytochemistry on day 8. Endothelial and epicardial differentiation was performed by seeding cells at a density of 6x10 4 cells/cm 2 on Matrigel-coated (1:100 dilution) 96-well plates (BD Falcon) in BPEL medium supplemented with 5 μM SB431542, 100 ng/ml IGF-1, 50 ng/ml VEGF165 (Proteintech) and 50 pg/ml bFGF. This medium was refreshed at day 2, 4 and 6. On day 8, cells were dissociated with TrypLE for flow cytometric analysis or fixed for immunocytochemistry. Flow cytometry and sorting EBs were dissociated into single cells with TrypLE (1X or 10X depending on the stage). Staining was performed in FACS buffer (0.5% BSA/1 mM EDTA in PBS) with the following antibodies ( Table S3 ): anti-PDGFRα-BV421 (BD Biosciences 562799), anti-TRA181-FITC (BD Biosciences cat no. 560883), and CD235a-APC (BD Biosciences cat no. 561775). Staining was performed at room temperature for 15 minutes followed by washing with FACS buffer. Measurements were performed using LSR Fortessa (BD Biosciences). A viability dye was added to exclude dead cells (DAPI/DR/7-AAD, depending on the panel). Cell sorting was performed on a BD Aria or a BD Influx (BD Biosciences), with samples maintained at 4°C throughout. Flow cytometric data was analyzed with FlowJo (v10). RNA isolation and quantitative reverse transcription PCR (RT-qPCR) RNA was extracted using TRIzol followed by purification using the Purelink RNA Mini Kit (Thermo Fisher Scientific), according to the manufacturer’s instructions. Genomic DNA was removed by on-column DNAse digestion. cDNA was generated from 1 µg of total RNA using the Tetro cDNA Synthesis Kit (Bioline) according to the manufacturer’s instructions. qPCR was performed using Power SYBR™ Green PCR Master Mix (Thermo Fisher Scientific). Relative gene expression was obtained by normalising to the housekeeping genes RPLPO and GUSB. Thermocycling was performed on a 7900HT qPCR system (Applied Biosystems) or QuantStudio 1 (Thermo Fisher Scientific). The ΔΔCt method was used to calculate relative gene expression. Primers are listed in Table S4 . Bulk RNA-seq Quality and integrity of the RNA samples were assessed using a 4200 TapeStation (Agilent Technologies) and then libraries generated using the Illumina® Stranded mRNA Prep Ligation kit (Illumina, Inc.) according to the manufacturer’s protocol. Briefly, total RNA (typically 0.025-1 µg) was used as input material from which polyadenylated mRNA was purified using poly-T, oligo-attached, magnetic beads. Next, the mRNA was fragmented under elevated temperature and then reverse transcribed into first strand cDNA using random hexamer primers and in the presence of Actinomycin D (thus improving strand specificity whilst mitigating spurious DNA-dependent synthesis). Following removal of the template RNA, second strand cDNA was then synthesized to yield blunt-ended, double-stranded cDNA fragments. Strand specificity was maintained by the incorporation of deoxyuridine triphosphate (dUTP) in place of dTTP to quench the second strand during subsequent amplification. Following a single adenine (A) base addition, adapters with a corresponding, complementary thymine (T) overhang were ligated to the cDNA fragments. Pre-index anchors were then ligated to the ends of the double-stranded cDNA fragments to prepare them for dual indexing. A subsequent PCR amplification step was then used to add the index adapter sequences to create the final cDNA library. The adapter indices enabled the multiplexing of the libraries, which were pooled prior to cluster generation using a cBot instrument. The loaded flow-cell was then paired-end sequenced (76 + 76 cycles, plus indices) on an Illumina HiSeq4000 instrument. Finally, the output data was demultiplexed and BCL-to-Fastq conversion performed using Illumina’s bcl2fastq software (v2.20.0.422). Input sequence was mapped using STAR (v2.7.7a). Gene annotation used was Human Genocode 37 for hg38/GRCh38 Ensembl (v103). Read counts were generated by STAR. Ensembl gene identifiers and gene names were assigned to each gene. For day 3 analysis, DESeq2 (v3.19) was used to identify genes differentially expressed in either HAND1-BFP high or HAND1-BFP med fractions compared to the HAND1-BFP neg fraction (FDR < 0.05). Genes responding to dox treatment independently of HAND1 were identified and subtracted by performing control experiments in a line not carrying the TO-HAND1 lentivirus. All visualization of RNA-seq analysis was carried out in R Studio (v4.0.5). ATAC-Seq After cell sorting, samples were prepared for ATAC-seq analysis using an ATAC-Seq Kit (Active Motif) as per the manufacturer’s instructions. Briefly, cells were gently washed with ice cold PBS and lysed with ATAC lysis buffer. Tagmentation master mix was prepared (Tagmentation buffer, PBS, 0.5% Digitonin, 10% Tween 20, Assembled Transposomes) and tagmentation carried out at 37°C for 30 minutes at 800 rpm. DNA was purified and then PCR amplified using unique combinations of i7 and i5 indexing primers per sample. After PCR amplification, samples were purified using SPRI beads and quality determined using Agilent Tapestation 2200. Barcoding primers are listed in Table S5 . The adapter indices enabled the multiplexing of the libraries, which were pooled prior to cluster generation using a cBot instrument. The loaded flow-cell was then paired-end sequenced (76 + 76 cycles, plus indices) on an Illumina HiSeq4000 instrument. Finally, the output data was demultiplexed and BCL-to-Fastq conversion performed using Illumina’s bcl2fastq software, (v.2.20.0.422). Briefly, peaks were called using MACS2 (day 7.5 samples) or MACS3 (day 3 samples) as 200 bp regions centered on the MACS peak summit. Overlapping regions were combined to create longer regions. Data were normalized based on sequencing depth. Peak annotation was performed using UpsetR (v1.4.0) and ChIPSeeker (v3.19). For the day 3 samples (HAND1-BFP+and HAND1-BFPneg), differential accessibility analysis was carried out using DiffBind (v3) and DESeq2 (v3.19) with a default threshold of false discovery rate (FDR) <0.05. Single cell RNA-seq: Sample and library preparation EBs were dissociated to single cells as described above. Cells from different conditions were combined as follows: 25% SB-treated cells, 25% DMH1-treated cells and 50% control (vehicle-only) cells. Samples were collected at day 3, 4, 5, 6, 7, 8 and 10. Additional wild-type replicates at days 7 and 8 were included. A total of 1-2 x 10 6 cells were used for each sample. The Single Cell Fixed RNA Sample Preparation Kit (10X Genomics) was used to fix and cryopreserve the cells according to the manufacturer’s instructions. In brief, cells were fixed with Fixation Buffer containing 4% formaldehyde (ThermoFisher Scientific) overnight at 4°C. 24 hours later, cells were quenched with freshly prepared quench buffer and stored at -80°C. Frozen cells were thawed and washed in 0.5X PBS/0.02% BSA (Miltenyi)/0.2 U/ml RNase Inhibitor (Roche). Barcoding was performed using the Chromium Next GEM Single Cell Fixed RNA Human Transcriptome Probe Kit (10X Genomics) according to the manufacturer’s instructions. Up to 2 x 10 6 cells per timepoint and genotype (wild-type and HAND1-null) were re-suspended in Hybridisation Mix with a single Human WTA barcode (10x Genomics, Inc. Pleasanton, USA) according to the manufacturer’s protocol (CG000527 Rev C). Briefly, samples were incubated overnight with probe pairs designed to hybridize to mRNA, samples were pooled and washed to remove unbound probes. Nanoliter-scale Gel Beads-in-emulsion (GEMs) were generated by loading barcoded Gel Beads, a master mix containing pooled, probe-hybridized cells, and partitioning oil onto a Chromium Next GEM chip. Cells were delivered at a limiting dilution, such that the majority (90-99%) of generated GEMs contain no cell, while the remainder largely contain a single cell. The Gel Beads were then dissolved, primers, including an Illumina TruSeq Read 1 sequence, a 16-nucleotide 10x Barcode, a 12-nucleotide unique molecular identifier (UMI) and partial Capture Sequence 1, were released, and any co-partitioned cells lysed. Following GEM generation, a ligation step sealed the nick between the left- and right-hand probe and subsequently the Gel Bead primer was hybridised to the capture sequence on the ligated probe pair and extended by a polymerase adding the UMI, 10x Barcode and partial Read 1 sequence before a heat-denaturation step inactivated the enzymes. Post-incubation, GEMs were broken and the recovered ligated product pre-amplified before being cleaned up by SPRIselect. Illumina-compatible sequencing libraries were constructed by adding P5, P7, i5 and i7 sample indexes, and Illumina Small Read 2 sequences to the 10x barcoded, ligated probe products via Sample Index PCR followed by a SPRIselect size-selection. scRNA-seq analysis: Library sequencing The resulting sequencing library comprised standard Illumina paired-end constructs flanked with P5 and P7 sequences. The 16 bp 10x Barcode and 12 bp UMI were encoded in Read 1, while Read 2 sequenced the ligated probe insert, constant sequence and the 8bp Probe Barcode. Sample indexes were incorporated as the i5 and i7 index reads. Paired-end sequencing (28:90) was performed on the Illumina NovaSeq6000 platform. The .bcl sequence data were processed for QC purposes using bcl2fastq software (v. 2.20.0.422) and the resulting .fastq files assessed using FastQC (v. 0.11.3), FastqScreen (v. 0.14.0) and FastqStrand (v. 1.13.0) prior to pre-processing with the CellRanger pipeline. scRNA-seq analysis: Read mapping to genome and cell filtering Raw sequencing data were processed using the 10x Genomics Cell Ranger pipeline (v7.1.0). Base call (BCL) files generated by the sequencer were converted to FASTQ files using “cellranger mkfastq”. The FASTQ files were then mapped against the pre-built human reference package from 10X Genomics (GRCh38-2020-A) using “cellranger multi” to demultiplex and produce the gene-cell barcode matrix for individual samples. The single-cell data were processed in R environment (v4.1) following the workflow documented in Orchestrating Single-Cell Analysis with Bioconductor. 57 Briefly, for each sample, the HDF5 file generated by Cell Ranger was imported into R to create a SingleCellExperiment object. A combination of median absolute deviation (MAD), as implemented by the “isOutlier” function in the scuttle R package (v1.4.0) and exact thresholds was used to identify and subsequently remove low quality cells before data integration. scRNA-seq: Data integration, visualization, cell clustering and annotation The log-normalized expression values of the combined data were computed using the “multiBatchNorm” function from the batchelor R package (v1.10.0). The per-gene variance of the log-expression profile was modelled using the “modelGeneVarByPoisson” function from the scran R package (v1.22.1) and top 5000 highly variable genes were selected. The mutual nearest neighbors (MNN) approach implemented by the “fastMNN” function from the batchelor R package was used to perform batch correction. For the combined analysis, the first 50 dimensions of the MNN low-dimensional corrected coordinates for all cells were used as input to produce the uniform manifold approximation and projection (UMAP) using the “runUMAP” function from the scater R package (v1.22.0). Putative cell clusters were identified using the Leiden algorithm from the igraph R package (v1.5.1). Cluster-specific markers were identified using the “findMarkers” function from the scran R package, which performs pairwise t-tests between clusters. scRNA-seq: Lineage trajectory analysis and marker gene expression Separate wild-type and HAND1-null analyses were performed using Monocle3 (v.1.3.4). The raw counts were pre-processed using the “preprocess_cds” function. The HVGs were as identified from the combined analysis but filtered to include only genes implicated in heart development, 58 excluding cell cycle-related genes (from the Seurat package dataset ‘cc.genes.updated.2019’). Dimensionality reduction by UMAP and in the examples was performed using a minimum distance of 1 and sqrt(n cells) nearest neighbours. Functions “learn_graph” and ”order_cells” were performed using default settings. To find consensus trajectories the UMAP settings were varied, and the dominant paths were taken. For gene expression plots a minimum expression threshold of 0.7 was used. HAND1 target gene identification For the HAND1-high against HAND1-low cell comparison in the wild-type mesoderm, the HAND1-high population was defined as cells with HAND1 expression in the top quantile, and the HAND1-low population as cells with expression in the bottom quantile. For the wild-type against HAND1-null mesoderm comparison, the mesoderm clusters of both were selected. Differential expression analysis of HAND1-null against wild-type and HAND1-high against HAND1-low in the mesoderm was performed using the “nbTestSH” function (FDR < 0.05) from the sSeq package. 59 Genes affected concordantly were considered putative HAND1 targets. To generate the final regulons, these targets were further filtered using the HAND1 ChIP-seq data and only genes with a HAND1 binding site within 500kb of the gene’s transcription start site were retained. Gene ontology analysis Gene ontology (GO) enrichment analysis of differentially expressed genes was performed using the enrichR R package (v3.2). Results were ranked by the combined scores (calculated as log(P-value) biormultiplied by z-score) and terms with an adjusted P-value < 0.05 were considered significant. The GOplot package was used to generate a circularly composited overview of HAND1 target genes and their assigned Biological Processes (2023) GO terms. Gene regulatory network analysis by SCENIC After QC filtering and pre-processing as above, the PySCENIC pipeline was followed as described. 35 Network inference was performed using the GRNBoost2 algorithm. For motif enrichment analysis and TF-regulon prediction by cisTarget, two ranking databases were used: ‘hg38 refseq-r80 10kb_up_and_down_tss.genes_vs_motifs.rankings.feather’ and ‘encode_20190621 ChIP_seq_transcription_factor.hg38 refseq-r80 10kb_up_and_down_tss.max.genes_vs_tracks.rankings.feather’. The wild-type, cardiac progenitors-only and HAND1-null cells were analysed separately. GRNBoost2 and cisTarget were run 5 times using different seed values and only regulators identified in all 5 runs were retained. Regulon target gene lists were aggregated and filtered to include all genes identified in at least 2 runs. Where unique regulators were identified in the cardiac progenitors or the HAND1-null cells they were added to the final list of regulons. Custom regulons were defined for HAND1 and MESP1 as described in the results. The MESP1 regulon was defined using published RNA-seq data from a Mesp1 overexpression system. 37 Genes upregulated at 12 hours (dox vs control; ‘batch 1’) were identified using DESeq2 (v3.19) (FDR 1.5) and intersected with SCENIC adjacencies in our data to enrich for direct targets. AUCell was performed in R to score the activity of these regulons across all cells. The automatically selected regulon activity thresholds were used for plotting. Regulon specificity scores (RSS) were calculated and those above the 95 th percentile marked. Constructing a gene regulatory network model for HAND1 A hierarchical gene regulatory network aiming to explain early HAND1-dependent regulatory events was constructed using the data from SCENIC together with the HAND1 target genes, identified as above, and HAND1 ChIP-seq data. The model was drawn using BioTapestry. Immunoblotting EBs were lysed in RIPA buffer (Sigma) supplemented with Halt™ Protease and Phosphatase Inhibitor Cocktail (Thermo Fisher Scientific), DNAseI and 5 mM EDTA. Protein lysate was quantified using the BCA assay (Biorad). 10-20 µg of protein lysate was run through MiniPROTEAN TGX Stain-free gels 4-12% gels (Biorad). Protein was transferred to nitrocellulose membranes using the Turbo transfer system (Biorad) then blocked in blocking buffer (5% milk/0.1% Tween-20/PBS). Membranes were incubated overnight at 4°C in blocking buffer with the appropriate antibodies ( Table S3 ). Appropriate HRP-conjugated secondary antibodies were added for 1 hour at room temperature in blocking buffer. Blots were developed using ECL™ Prime Western Blotting Detection Reagents (Amersham™) and imaged using ChemiDoc™ Imaging System (Biorad). Immunocytochemistry and microscopy EBs/cells grown on coverslips were fixed with 10% neutral buffered formalin (Sigma-Aldrich) for 20/15 minutes, respectively, followed by 3X PBS washes. EBs were permeabilised with 0.5% TritonX-100/PBS for 1 hour at room temperature, whereas cells on coverslips were permeabilised for 8 minutes. Cells were blocked with with 4% donkey serum/PBS for 1 hour at room temperature. Antibodies ( Table S3 ) were incubated overnight at 4°C in blocking buffer. The next day, EBs were washed 3X with PBS-tween, followed by incubation with AlexaFluor secondary antibodies ( Table S3 ) for 1 hour at room temperature. EBs were stained with 1 µg/ml Hoechst for 15 minutes at room temperature, then washed with PBS/0.1% BSA. EBs were mounted onto coverslips using ProLong™ Gold Antifade Mountant (Invitrogen). Images were taken with an Evos™ 5000™ Imaging System (Thermo Fisher Scientific) or on a Leica SP8 confocal microscope. Chromatin Immunoprecipitation (ChIP) and ChIP-seq ChIP was performed largely as previously described. 60 In brief, cells were differentiated as EBs as described above to D3. 5 µM of SB (or an equivalent volume of DMSO for vehicle-only controls) was added on D2-D3. EBs were dissociated into single cells with TryPLE and fixed with 1% formaldehyde (Pierce) for 10 minutes at room temperature, then quenched with 0.125 M Glycine for 5 minutes at room temperature. Cells were centrifuged and washed three times with cold PBS, then the pellets snap-frozen in a dry ice/ethanol bath and stored at -80°C until required. Pellets were thawed on ice and re-suspended in Sonication/IP buffer (140 mM NaCl, 50 mM HEPES pH 7.5, 1 mM EDTA pH 8.0, 1%, 0.1% SDS, 0.1% sodium deoxcycholate, 1X protease inhibitor). 3 x 10 6 cells in 150 µl per TPX 1.5 ml tube (Diagenode) was sonicated for 40 cycles (30 seconds on/30 seconds off on High Power) using a Bioruptor (Diagenode). 10 µg of anti-AM-Tag antibody (Active Motif) was incubated with 40 µl Protein G Dynabeads (Invitrogen) for 3 hours at room temperature on a rotator. Sonicated chromatin was diluted with 1.5 ml of Chromatin Dilution Buffer (25 mM Tris pH 7.5, 5 mM EDTA, 1% Triton X-100, 0.1% SDS) plus protease inhibitors and centrifuged at 13,600 x g for 30 minutes at 4°C. The supernatant was transferred to a protein low-bind tube (Eppendorf) and incubated with antibody/bead complex overnight on a rotator in the cold room. Wash steps and reverse cross-linking was performed as previously described. 60 Libraries were prepared using MicroPlex Library Preparation Kit v3 (Diagenode) according to the Manufacturer’s instructions. Library quality was assessed using a 4200 TapeStation (Agilent Technologies) and DNA quantified using the Qubit. The adapter indices enabled the multiplexing of the libraries, which were pooled prior to loading on to the appropriate flow-cell. This was then paired-end sequenced (59 + 59 cycles, plus indices) on an Illumina NovaSeq6000 instrument. Finally, the output data was demultiplexed and BCL-to-Fastq conversion performed using Illumina’s bcl2fastq software (v2.20.0.422). ChIP-seq analysis Reads were trimmed using Trim Galore! to remove sequencing adaptors and poor-quality bases. Reads were aligned to the human genome (hg38) using Bowtie2 (v2.5.3). Peaks were called using MACS2 (v2.2.9.1) using default settings. Overlaps with ATAC-seq regions were identified using the R IRanges package (v.2.36). Peak annotation was performed using ChIPSeeker (v3.19). Transcription factor motif enrichment analysis From peaksets obtained using MACS, we used Bedtools multiinter to combine the replicates of each sample and return the replicated peaks. For the day 7.5 samples, regions unique to each population defined only by SOX17-Tom and NKX2-5-GFP were used for the motif enrichment analysis. In the final analysis, to obtain comparable significance values, all peaksets were reduced to the same number of regions (4562) by random sampling. Multiple samples were tested to ensure results were representative. Motif enrichment analysis was performed using Homer (v4.4) ‘findMotifsGenome.pl’. Options ‘-mask’ and ‘-size 200’ and default background regions were used. A p-value of < 1E-6 was selected as significant. Statistical Analysis Statistical details of experiments can be found in the Figure legends, including p-values and numbers of samples analyzed. Unless noted otherwise, adjusted p-values of below 0.05 were treated as significant. Plots were generated with ggplot2 (v3.5.1) in R or GraphPad Prism 9. Data availability Sequencing data has been deposited in the ArrayExpress database at EMBL-EBI under the following accession numbers: RNA-seq data (day 7.5): E-MTAB-14276 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14276?key=ed2cc2f0-ed66-4bdd-85eb-76593189c9e7 ATAC-seq data (day 7.5): E-MTAB-14277 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14277?key=40259d46-fd3f-410a-ab11-ec601106dcf7 RNA-seq data (day 3): E-MTAB-14278 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14278?key=f7129751-c31b-454c-bafe-0e27f5eed0ff ATAC-seq data (day 3): E-MTAB-14279 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14279?key=ce4726ec-91a6-4419-baf6-0b40c0750748 ChIP-seq data (day 3): E-MTAB-14280 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14280?key=a545b95f-db45-4b89-b179-a562e9ad7fd3 Single cell RNA-seq data: E-MTAB-14285 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14285?key=31df2b45-a691-4a51-85d7-75da9401a202 Acknowledgements This work was supported and funded by the UKRI Future Leaders Fellowship (MR/T041668/1 to MJB); the British Heart Foundation (FS/PhD/23/29429 to MJB and NB); an EPSRC/MRC PhD studentship (EP/L014904/1 to NP); and the Biotechnology and Biological Sciences Research Council (BB/T007761 to NB and BB/X016684 to NB and MJB). We thank all the core facilities staff who supported this research. We thank David Elliot and Christine Mummery for providing the NKX2-5-GFP hESC reporter line. We thank Andy Sharrocks for reviewing the manuscript. Footnotes https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14276?key=ed2cc2f0-ed66-4bdd-85eb-76593189c9e7 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14277?key=40259d46-fd3f-410a-ab11-ec601106dcf7 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14278?key=f7129751-c31b-454c-bafe-0e27f5eed0ff https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14279?key=ce4726ec-91a6-4419-baf6-0b40c0750748 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14280?key=a545b95f-db45-4b89-b179-a562e9ad7fd3 https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-14285?key=31df2b45-a691-4a51-85d7-75da9401a202 References 1. ↵ Meilhac , S.M. , and Buckingham , M.E . ( 2018 ). The deployment of cell lineages that form the mammalian heart . Nat. Rev. Cardiol . 15 , 705 – 724 . doi: 10.1038/s41569-018-0086-9 . OpenUrl CrossRef PubMed 2. ↵ Houyel , L. , and Meilhac , S.M . ( 2021 ). Heart Development and Congenital Structural Heart Defects . Annu. Rev. Genomics Hum. Genet . 22 , 257 – 284 . doi: 10.1146/annurev-genom-083118-015012 . OpenUrl CrossRef 3. ↵ Morton , S.U. , Quiat , D. , Seidman , J.G. , and Seidman , C.E . ( 2022 ). Genomic frontiers in congenital heart disease . Nat. Rev. Cardiol . 19 , 26 – 42 . doi: 10.1038/s41569-021-00587-4 . OpenUrl CrossRef PubMed 4. ↵ Bardot , E.S. , and Hadjantonakis , A.-K . ( 2020 ). Mouse gastrulation: Coordination of tissue patterning, specification and diversification of cell fate . Mech. Dev . 163 , 103617 . doi: 10.1016/j.mod.2020.103617 . OpenUrl CrossRef PubMed 5. ↵ Devine , W.P. , Wythe , J.D. , George , M. , Koshiba-Takeuchi , K. , and Bruneau , B.G . ( 2014 ). Early patterning and specification of cardiac progenitors in gastrulating mesoderm . eLife 3 , e03848 . doi: 10.7554/eLife.03848 . OpenUrl CrossRef PubMed 6. ↵ Lescroart , F. , Chabab , S. , Lin , X. , Rulands , S. , Paulissen , C. , Rodolosse , A. , Auer , H. , Achouri , Y. , Dubois , C. , Bondue , A. , et al. ( 2014 ). Early lineage restriction in temporally distinct populations of Mesp1 progenitors during mammalian heart development . Nat. Cell Biol . 16 , 829 – 840 . doi: 10.1038/ncb3024 . OpenUrl CrossRef PubMed Web of Science 7. ↵ Meilhac , S.M. , Esner , M. , Kelly , R.G. , Nicolas , J.-F. , and Buckingham , M.E . ( 2004 ). The clonal origin of myocardial cells in different regions of the embryonic mouse heart . Dev. Cell 6 , 685 – 698 . OpenUrl CrossRef PubMed Web of Science 8. ↵ Dominguez , M.H. , Krup , A.L. , Muncie , J.M. , and Bruneau , B.G . ( 2023 ). Graded mesoderm assembly governs cell fate and morphogenesis of the early mammalian heart . Cell 186 , 479 – 496 .e23. doi: 10.1016/j.cell.2023.01.001 . OpenUrl CrossRef PubMed 9. ↵ Ivanovitch , K. , Soro-Barrio , P. , Chakravarty , P. , Jones , R.A. , Bell , D.M. , Gharavy , S.N.M. , Stamataki , D. , Delile , J. , Smith , J.C. , and Briscoe , J . ( 2021 ). Ventricular, atrial, and outflow tract heart progenitors arise from spatially and molecularly distinct regions of the primitive streak . PLOS Biol . 19 , e3001200 . doi: 10.1371/journal.pbio.3001200 . OpenUrl CrossRef PubMed 10. ↵ Tyser , R.C.V. , Mahammadov , E. , Nakanoh , S. , Vallier , L. , Scialdone , A. , and Srinivas , S . ( 2021 ). Single-cell transcriptomic characterization of a gastrulating human embryo . Nature 600 , 285 – 289 . doi: 10.1038/s41586-021-04158-y . OpenUrl CrossRef 11. ↵ Tyser , R.C.V. , Ibarra-Soria , X. , McDole , K. , Arcot Jayaram , S. , Godwin , J. , van den Brand , T.A.H. , Miranda , A.M.A. , Scialdone , A. , Keller , P.J. , Marioni , J.C. , et al. ( 2021 ). Characterization of a common progenitor pool of the epicardium and myocardium . Science 371 , eabb2986 . doi: 10.1126/science.abb2986 . OpenUrl Abstract / FREE Full Text 12. ↵ Zhang , Q. , Carlin , D. , Zhu , F. , Cattaneo , P. , Ideker , T. , Evans , S.M. , Bloomekatz , J. , and Chi , N.C . ( 2021 ). Unveiling Complexity and Multipotentiality of Early Heart Fields . Circ. Res . 129 , 474 – 487 . doi: 10.1161/CIRCRESAHA.121.318943 . OpenUrl CrossRef PubMed 13. ↵ Firulli , A.B. , McFadden , D.G. , Lin , Q. , Srivastava , D. , and Olson , E.N . ( 1998 ). Heart and extra-embryonic mesodermal defects in mouse embryos lacking the bHLH transcription factor Hand1 . Nat. Genet . 18 , 266 – 270 . doi: 10.1038/ng0398-266 . OpenUrl CrossRef PubMed Web of Science 14. ↵ Firulli , B.A. , George , R.M. , Harkin , J. , Toolan , K.P. , Gao , H. , Liu , Y. , Zhang , W. , Field , L.J. , Liu , Y. , Shou , W. , et al. ( 2020 ). HAND1 loss-of-function within the embryonic myocardium reveals survivable congenital cardiac defects and adult heart failure . Cardiovasc. Res . 116 , 605 – 618 . doi: 10.1093/cvr/cvz182 . OpenUrl CrossRef 15. ↵ Riley , P. , Anaon-Cartwight , L. , and Cross , J.C . ( 1998 ). The Hand1 bHLH transcription factor is essential for placentation and cardiac morphogenesis . Nat. Genet . 18 , 271 – 275 . doi: 10.1038/ng0398-271 . OpenUrl CrossRef PubMed Web of Science 16. ↵ Risebro , C.A. , Smart , N. , Dupays , L. , Breckenridge , R. , Mohun , T.J. , and Riley , P.R . ( 2006 ). Hand1 regulates cardiomyocyte proliferation versus differentiation in the developing heart . Development 133 , 4595 – 4606 . doi: 10.1242/dev.02625 . OpenUrl Abstract / FREE Full Text 17. ↵ Vincentz , J.W. , Toolan , K.P. , Zhang , W. , and Firulli , A.B . ( 2017 ). Hand factor ablation causes defective left ventricular chamber development and compromised adult cardiac function . PLOS Genet . 13 , e1006922 . doi: 10.1371/journal.pgen.1006922 . OpenUrl CrossRef 18. ↵ Li , L. , Wang , J. , Liu , X.-Y. , Liu , H. , Shi , H.-Y. , Yang , X.-X. , Li , N. , Li , Y.-J. , Huang , R.-T. , Xue , S. , et al. ( 2017 ). HAND1 loss-of-function mutation contributes to congenital double outlet right ventricle . Int. J. Mol. Med . 39 , 711 – 718 . doi: 10.3892/ijmm.2017.2865 . OpenUrl CrossRef 19. ↵ Meier , A.B. , Zawada , D. , De Angelis , M.T. , Martens , L.D. , Santamaria , G. , Zengerle , S. , Nowak-Imialek , M. , Kornherr , J. , Zhang , F. , Tian , Q. , et al. ( 2023 ). Epicardioid single-cell genomics uncovers principles of human epicardium biology in heart development and disease . Nat. Biotechnol . 41 , 1787 – 1800 . doi: 10.1038/s41587-023-01718-7 . OpenUrl CrossRef 20. Schmidt , C. , Deyett , A. , Ilmer , T. , Haendeler , S. , Caballero , A.T. , Novatchkova , M. , Netzer , M.A. , Ginistrelli , L.C. , Juncosa , E.M. , Bhattacharya , T. , et al. ( 2023 ). Multi-chamber cardioids unravel human heart development and cardiac defects . Cell 186 , 5587 – 5605 .e27. doi: 10.1016/j.cell.2023.10.030 . OpenUrl CrossRef PubMed 21. ↵ Lee , J.H. , Protze , S.I. , Laksman , Z. , Backx , P.H. , and Keller , G.M . ( 2017 ). Human Pluripotent Stem Cell-Derived Atrial and Ventricular Cardiomyocytes Develop from Distinct Mesoderm Populations . Cell Stem Cell 21 , 179 – 194 .e4. doi: 10.1016/j.stem.2017.07.003 . OpenUrl CrossRef PubMed 22. ↵ Birket , M.J. , Ribeiro , M.C. , Verkerk , A.O. , Ward , D. , Leitoguinho , A.R. , den Hartogh , S.C. , Orlova , V.V. , Devalla , H.D. , Schwach , V. , Bellin , M. , et al. ( 2015 ). Expansion and patterning of cardiovascular progenitors derived from human pluripotent stem cells . Nat. Biotechnol . 33 , 970 – 979 . doi: 10.1038/nbt.3271 . OpenUrl CrossRef PubMed 23. ↵ Bardot , E. , Calderon , D. , Santoriello , F. , Han , S. , Cheung , K. , Jadhav , B. , Burtscher , I. , Artap , S. , Jain , R. , Epstein , J. , et al. ( 2017 ). Foxa2 identifies a cardiac progenitor population with ventricular differentiation potential . Nat. Commun . 8 , 14428 . doi: 10.1038/ncomms14428 . OpenUrl CrossRef PubMed 24. ↵ Wei , Q. , and Condie , B.G . ( 2011 ). A Focused In Situ Hybridization Screen Identifies Candidate Transcriptional Regulators of Thymic Epithelial Cell Development and Function . PLOS ONE 6 , e26795 . doi: 10.1371/journal.pone.0026795 . OpenUrl CrossRef PubMed 25. ↵ Marques , I.J. , Ernst , A. , Arora , P. , Vianin , A. , Hetke , T. , Sanz-Morejón , A. , Naumann , U. , Odriozola , A. , Langa , X. , Andrés-Delgado , L. , et al. ( 2022 ). Wt1 transcription factor impairs cardiomyocyte specification and drives a phenotypic switch from myocardium to epicardium . Development 149 , dev200375 . doi: 10.1242/dev.200375 . OpenUrl CrossRef 26. ↵ Martínez-Estrada , O.M. , Lettice , L.A. , Essafi , A. , Guadix , J.A. , Slight , J. , Velecela , V. , Hall , E. , Reichmann , J. , Devenney , P.S. , Hohenstein , P. , et al. ( 2010 ). Wt1 is required for cardiovascular progenitor cell formation through transcriptional control of Snail and E-cadherin . Nat. Genet . 42 , 89 – 93 . doi: 10.1038/ng.494 . OpenUrl CrossRef PubMed Web of Science 27. ↵ Trapnell , C. , Cacchiarelli , D. , Grimsby , J. , Pokharel , P. , Li , S. , Morse , M. , Lennon , N.J. , Livak , K.J. , Mikkelsen , T.S. , and Rinn , J.L . ( 2014 ). The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells . Nat. Biotechnol . 32 , 381 – 386 . doi: 10.1038/nbt.2859 . OpenUrl CrossRef PubMed 28. ↵ Aanhaanen , W.T.J. , Brons , J.F. , Domínguez , J.N. , Rana , M.S. , Norden , J. , Airik , R. , Wakker , V. , de Gier-de Vries , C. , Brown , N.A. , Kispert , A. , et al. ( 2009 ). The Tbx2+ Primary Myocardium of the Atrioventricular Canal Forms the Atrioventricular Node and the Base of the Left Ventricle . Circ. Res . 104 , 1267 – 1274 . doi: 10.1161/CIRCRESAHA.108.192450 . OpenUrl Abstract / FREE Full Text 29. ↵ Liang , X. , Wang , G. , Lin , L. , Lowe , J. , Zhang , Q. , Bu , L. , Chen , Y. , Chen , J. , Sun , Y. , and Evans , S.M . ( 2013 ). HCN4 Dynamically Marks the First Heart Field and Conduction System Precursors . Circ. Res . 113 , 399 – 407 . doi: 10.1161/CIRCRESAHA.113.301588 . OpenUrl Abstract / FREE Full Text 30. ↵ Cai , C.-L. , Liang , X. , Shi , Y. , Chu , P.-H. , Pfaff , S.L. , Chen , J. , and Evans , S . ( 2003 ). Isl1 identifies a cardiac progenitor population that proliferates prior to differentiation and contributes a majority of cells to the heart . Dev. Cell 5 , 877 – 889 . OpenUrl CrossRef PubMed Web of Science 31. High , F.A. , Jain , R. , Stoller , J.Z. , Antonucci , N.B. , Lu , M.M. , Loomes , K.M. , Kaestner , K.H. , Pear , W.S. , and Epstein , J.A . ( 2009 ). Murine Jagged1/Notch signaling in the second heart field orchestrates Fgf8 expression and tissue-tissue interactions during outflow tract development . J. Clin. Invest . 119 , 1986 – 1996 . doi: 10.1172/JCI38922 . OpenUrl CrossRef PubMed Web of Science 32. ↵ Kelly , R.G. , Brown , N.A. , and Buckingham , M.E . ( 2001 ). The Arterial Pole of the Mouse Heart Forms from Fgf10-Expressing Cells in Pharyngeal Mesoderm . Dev. Cell 1 , 435 – 440 . doi: 10.1016/S1534-5807(01)00040-5 . OpenUrl CrossRef PubMed Web of Science 33. ↵ Bulger , E.A. , McDevitt , T.C. , and Bruneau , B.G . ( 2024 ). CDX2 dose-dependently influences the gene regulatory network underlying human extraembryonic mesoderm development . Biol. Open 13 , bio060323 . doi: 10.1242/bio.060323 . OpenUrl CrossRef 34. ↵ Aibar , S. , González-Blas , C.B. , Moerman , T. , Huynh-Thu , V.A. , Imrichova , H. , Hulselmans , G. , Rambow , F. , Marine , J.-C. , Geurts , P. , Aerts , J. , et al. ( 2017 ). SCENIC: single-cell regulatory network inference and clustering . Nat. Methods 14 , 1083 – 1086 . doi: 10.1038/nmeth.4463 . OpenUrl CrossRef PubMed 35. ↵ Van de Sande , B. , Flerin , C. , Davie , K. , De Waegeneer , M. , Hulselmans , G. , Aibar , S. , Seurinck , R. , Saelens , W. , Cannoodt , R. , Rouchon , Q. , et al. ( 2020 ). A scalable SCENIC workflow for single-cell gene regulatory network analysis - Nature Protocols . Nat. Protoc . 15 , 2247 – 2276 . doi: 10.1038/s41596-020-0336-2 . OpenUrl CrossRef PubMed 36. ↵ Tsankov , A.M. , Gu , H. , Akopian , V. , Ziller , M.J. , Donaghey , J. , Amit , I. , Gnirke , A. , and Meissner , A . ( 2015 ). Transcription factor binding dynamics during human ES cell differentiation . Nature 518 , 344 – 349 . doi: 10.1038/nature14233 . OpenUrl CrossRef PubMed 37. ↵ Lin , X. , Swedlund , B. , Ton , M.-L.N. , Ghazanfar , S. , Guibentif , C. , Paulissen , C. , Baudelet , E. , Plaindoux , E. , Achouri , Y. , Calonne , E. , et al. ( 2022 ). Mesp1 controls the chromatin and enhancer landscapes essential for spatiotemporal patterning of early cardiovascular progenitors . Nat. Cell Biol . 24 , 1114 – 1128 . doi: 10.1038/s41556-022-00947-3 . OpenUrl CrossRef 38. ↵ Charron , J. , Malynn , B.A. , Fisher , P. , Stewart , V. , Jeannotte , L. , Goff , S.P. , Robertson , E.J. , and Alt , F.W . ( 1992 ). Embryonic lethality in mice homozygous for a targeted disruption of the N-myc gene . Genes Dev . 6 , 2248 – 2257 . doi: 10.1101/gad.6.12a.2248 . OpenUrl Abstract / FREE Full Text 39. Davis , A.C. , Wims , M. , Spotts , G.D. , Hann , S.R. , and Bradley , A . ( 1993 ). A null c-myc mutation causes lethality before 10.5 days of gestation in homozygotes and reduced fertility in heterozygous female mice . Genes Dev . 7 , 671 – 682 . doi: 10.1101/gad.7.4.671 . OpenUrl Abstract / FREE Full Text 40. ↵ Muñoz-Martín , N. , Sierra , R. , Schimmang , T. , Campo , C.V. del, and Torres , M. ( 2019 ). Myc is dispensable for cardiomyocyte development but rescues Mycn-deficient hearts through functional replacement and cell competition . Development 146 . doi: 10.1242/dev.170753 . OpenUrl Abstract / FREE Full Text 41. ↵ Loh , K.M. , Ang , L.T. , Zhang , J. , Kumar , V. , Ang , J. , Auyeong , J.Q. , Lee , K.L. , Choo , S.H. , Lim , C.Y.Y. , Nichane , M. , et al. ( 2014 ). Efficient Endoderm Induction from Human Pluripotent Stem Cells by Logically Directing Signals Controlling Lineage Bifurcations . Cell Stem Cell 14 , 237 – 252 . doi: 10.1016/j.stem.2013.12.007 . OpenUrl CrossRef PubMed Web of Science 42. ↵ Teague , S. , Primavera , G. , Chen , B. , Liu , Z.-Y. , Yao , L. , Freeburne , E. , Khan , H. , Jo , K. , Johnson , C. , and Heemskerk , I . ( 2024 ). Time-integrated BMP signaling determines fate in a stem cell model for early human development . Nat. Commun . 15 , 1471 . doi: 10.1038/s41467-024-45719-9 . OpenUrl CrossRef 43. ↵ Croci , O. , Fazio , S.D. , Biagioni , F. , Donato , E. , Caganova , M. , Curti , L. , Doni , M. , Sberna , S. , Aldeghi , D. , Biancotto , C. , et al. ( 2017 ). Transcriptional integration of mitogenic and mechanical signals by Myc and YAP . Genes Dev . 31 , 2017 – 2022 . doi: 10.1101/gad.301184.117 . OpenUrl Abstract / FREE Full Text 44. ↵ LeBlanc , L. , Lee , B.-K. , Yu , A.C. , Kim , M. , Kambhampati , A.V. , Dupont , S.M. , Seruggia , D. , Ryu , B.U. , Orkin , S.H. , and Kim , J . ( 2018 ). Yap1 safeguards mouse embryonic stem cells from excessive apoptosis during differentiation . eLife 7 , e40167 . doi: 10.7554/eLife.40167 . OpenUrl CrossRef 45. ↵ Stefanovic , S. , Laforest , B. , Desvignes , J.-P. , Lescroart , F. , Argiro , L. , Maurel-Zaffran , C. , Salgado , D. , Plaindoux , E. , De Bono , C. , Pazur , K. , et al. ( 2020 ). Hox-dependent coordination of mouse cardiac progenitor cell patterning and differentiation . eLife 9 , e55124 . doi: 10.7554/eLife.55124 . OpenUrl CrossRef 46. ↵ Firulli , B.A. , McConville , D.P. , Byers , J.S. , Vincentz , J.W. , Barnes , R.M. , and Firulli , A.B . ( 2010 ). Analysis of a Hand1 hypomorphic allele reveals a critical threshold for embryonic viability . Dev. Dyn. Off. Publ. Am. Assoc. Anat . 239 , 2748 – 2760 . doi: 10.1002/dvdy.22402 . OpenUrl CrossRef PubMed 47. ↵ Darieva , Z. , Zarrineh , P. , Phillips , N. , Mallen , J. , Mora , A.G. , Donaldson , I. , Bridoux , L. , Douglas , M. , Henriques , S.F.D. , Schulte , D. , et al. ( 2024 ). The evolutionary ancient MEIS transcription factors actuate lineage-specific transcription to establish cardiac fate . Preprint at bioRxiv , doi: 10.1101/2024.07.18.604075 10.1101/2024.07.18.604075. OpenUrl Abstract / FREE Full Text 48. ↵ Rao , J. , Pfeiffer , M.J. , Frank , S. , Adachi , K. , Piccini , I. , Quaranta , R. , Araúzo-Bravo , M. , Schwarz , J. , Schade , D. , Leidel , S. , et al. ( 2016 ). Stepwise Clearance of Repressive Roadblocks Drives Cardiac Induction in Human ESCs . Cell Stem Cell 18 , 341 – 353 . doi: 10.1016/j.stem.2015.11.019 . OpenUrl CrossRef 49. ↵ Liang , Q. , Wang , S. , Zhou , X. , Li , Y. , Xing , S. , Sha , Y. , Yang , F. , Huang , W. , Liu , N. , Li , Z. , et al. ( 2022 ). Essential role of MESP1-RING1A complex in cardiac differentiation . Dev. Cell 57 , 2533 – 2549 .e7. doi: 10.1016/j.devcel.2022.10.009 . OpenUrl CrossRef 50. ↵ Schwach , V. , Gomes Fernandes , M. , Maas , S. , Gerhardt , S. , Tsonaka , R. , van der Weerd , L. , Passier , R. , Mummery , C.L. , Birket , M.J. , and Salvatori , D.C.F. ( 2020 ). Expandable human cardiovascular progenitors from stem cells for regenerating mouse heart after myocardial infarction . Cardiovasc. Res . 116 , 545 – 553 . doi: 10.1093/cvr/cvz181 . OpenUrl CrossRef 51. ↵ Bargehr , J. , Ong , L.P. , Colzani , M. , Davaapil , H. , Hofsteen , P. , Bhandari , S. , Gambardella , L. , Le Novère , N. , Iyer , D. , Sampaziotis , F. , et al. ( 2019 ). Epicardial cells derived from human embryonic stem cells augment cardiomyocyte-driven heart regeneration . Nat. Biotechnol . 37 , 895 – 906 . doi: 10.1038/s41587-019-0197-9 . OpenUrl CrossRef 52. Cheng , Y.-C. , Hsieh , M.L. , Lin , C.-J. , Chang , C.M.C. , Huang , C.-Y. , Puntney , R. , Wu Moy , A. , Ting , C.-Y. , Herr Chan , D.Z. , Nicholson , M.W. , et al. ( 2023 ). Combined Treatment of Human Induced Pluripotent Stem Cell–Derived Cardiomyocytes and Endothelial Cells Regenerate the Infarcted Heart in Mice and Non-Human Primates . Circulation 148 , 1395 – 1409 . doi: 10.1161/CIRCULATIONAHA.122.061736 . OpenUrl CrossRef 53. ↵ Liu , Y.-W. , Chen , B. , Yang , X. , Fugate , J.A. , Kalucki , F.A. , Futakuchi-Tsuchida , A. , Couture , L. , Vogel , K.W. , Astley , C.A. , Baldessari , A. , et al. ( 2018 ). Human embryonic stem cell– derived cardiomyocytes restore function in infarcted hearts of non-human primates . Nat. Biotechnol . 36 , 597 – 605 . doi: 10.1038/nbt.4162 . OpenUrl CrossRef PubMed 54. ↵ Elliott , D.A. , Braam , S.R. , Koutsis , K. , Ng , E.S. , Jenny , R. , Lagerqvist , E.L. , Biben , C. , Hatzistavrou , T. , Hirst , C.E. , Yu , Q.C. , et al. ( 2011 ). NKX2-5eGFP/w hESCs for isolation of human cardiac progenitors and cardiomyocytes . Nat. Methods 8 , 1037 – 1040 . doi: 10.1038/nmeth.1740 . OpenUrl CrossRef PubMed Web of Science 55. ↵ Lynch , A.T. , Douglas , M. , Kimber , S.J. , and Birket , M.J . ( 2024 ). The generation and validation of a dual cardiac HAND1-Tomato NKX2-5-GFP human embryonic stem cell line UMANe002-A-3 . Stem Cell Res . 77 , 103342 . doi: 10.1016/j.scr.2024.103342 . OpenUrl CrossRef 56. ↵ Ng , E.S. , Davis , R. , Stanley , E.G. , and Elefanty , A.G . ( 2008 ). A protocol describing the use of a recombinant protein-based, animal product-free medium (APEL) for human embryonic stem cell differentiation as spin embryoid bodies . Nat Protoc . 3 , 768 – 776 . doi: 10.1038/nprot.2008.42 . OpenUrl CrossRef PubMed 57. ↵ Amezquita , R.A. , Lun , A.T.L. , Becht , E. , Carey , V.J. , Carpp , L.N. , Geistlinger , L. , Marini , F. , Rue-Albrecht , K. , Risso , D. , Soneson , C. , et al. ( 2020 ). Orchestrating single-cell analysis with Bioconductor . Nat. Methods 17 , 137 – 145 . doi: 10.1038/s41592-019-0654-x . OpenUrl CrossRef PubMed 58. ↵ Zhou , W.-Z. , Li , W. , Shen , H. , Wang , R.W. , Chen , W. , Zhang , Y. , Zeng , Q. , Wang , H. , Yuan , M. , Zeng , Z. , et al. ( 2022 ). CHDbase: A comprehensive knowledgebase for congenital heart disease-related genes and clinical manifestations . Genomics Proteomics Bioinformatics . doi: 10.1016/j.gpb.2022.08.001 . OpenUrl CrossRef 59. ↵ Yu , D. , Huber , W. , and Vitek , O . ( 2013 ). Shrinkage estimation of dispersion in Negative Binomial models for RNA-seq experiments with small sample size . Bioinformatics 29 , 1275 – 1282 . doi: 10.1093/bioinformatics/btt143 . OpenUrl CrossRef PubMed Web of Science 60. ↵ Sullivan , A.E. , and Santos , S.D.M . ( 2020 ). An Optimized Protocol for ChIP-Seq from Human Embryonic Stem Cell Cultures . STAR Protoc . 1 , 100062 . doi: 10.1016/j.xpro.2020.100062 . OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted August 17, 2024. Download PDF Supplementary Material Data/Code Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following The level of HAND1 controls the specification of multipotent cardiac and extraembryonic progenitors Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share The level of HAND1 controls the specification of multipotent cardiac and extraembryonic progenitors Adam T Lynch , Naomi Phillips , Megan Douglas , Marta Dorgnach , I-Hsuan Lin , Antony D Adamson , Zoulfia Darieva , Jessica Whittle , Neil A Hanley , Nicoletta Bobola , Matthew J Birket bioRxiv 2024.08.15.607916; doi: https://doi.org/10.1101/2024.08.15.607916 Share This Article: Copy Citation Tools The level of HAND1 controls the specification of multipotent cardiac and extraembryonic progenitors Adam T Lynch , Naomi Phillips , Megan Douglas , Marta Dorgnach , I-Hsuan Lin , Antony D Adamson , Zoulfia Darieva , Jessica Whittle , Neil A Hanley , Nicoletta Bobola , Matthew J Birket bioRxiv 2024.08.15.607916; doi: https://doi.org/10.1101/2024.08.15.607916 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Developmental Biology Subject Areas All Articles Animal Behavior and Cognition (7651) Biochemistry (17746) Bioengineering (13928) Bioinformatics (42066) Biophysics (21499) Cancer Biology (18650) Cell Biology (25579) Clinical Trials (138) Developmental Biology (13409) Ecology (19947) Epidemiology (2067) Evolutionary Biology (24374) Genetics (15633) Genomics (22557) Immunology (17775) Microbiology (40505) Molecular Biology (17217) Neuroscience (88796) Paleontology (667) Pathology (2845) Pharmacology and Toxicology (4836) Physiology (7664) Plant Biology (15179) Scientific Communication and Education (2047) Synthetic Biology (4304) Systems Biology (9839) Zoology (2272)