Nr2f1a maintains atrial nkx2.5 expression to repress pacemaker identity within venous atrial cardiomyocytes of zebrafish
Abstract
Maintenance of cardiomyocyte identity is vital for normal heart development and function. However, our understanding of cardiomyocyte plasticity remains incomplete. Here, we show that sustained expression of the zebrafish transcription factor Nr2f1a prevents the progressive acquisition of ventricular cardiomyocyte (VC) and pacemaker cardiomyocyte (PC) identities within distinct regions of the atrium. Transcriptomic analysis of flow-sorted atrial cardiomyocytes (ACs) from nr2f1a mutant zebrafish embryos showed increased VC marker gene expression and altered expression of core PC regulatory genes, including decreased expression of nkx2.5, a critical repressor of PC differentiation. At the arterial (outflow) pole of the atrium in nr2f1a mutants, cardiomyocytes resolve to VC identity within the expanded atrioventricular canal. However, at the venous (inflow) pole of the atrium, there is a progressive wave of AC transdifferentiation into PCs across the atrium toward the arterial pole. Restoring Nkx2.5 is sufficient to repress PC marker identity in nr2f1a mutant atria and analysis of chromatin accessibility identified an Nr2f1a-dependent nkx2.5 enhancer expressed in the atrial myocardium directly adjacent to PCs. CRISPR/Cas9-mediated deletion of the putative nkx2.5 enhancer leads to a loss of Nkx2.5-expressing ACs and expansion of a PC reporter, supporting that Nr2f1a limits PC differentiation within venous ACs via maintaining nkx2.5 expression. The Nr2f-dependent maintenance of AC identity within discrete atrial compartments may provide insights into the molecular etiology of concurrent structural congenital heart defects and associated arrhythmias.
Editor's evaluation
Addressing the role of NR2F transcription factors in the fish heart, this important work provides novel insight into atrial chamber patterning. Solid evidence supports the identification of a new mechanism of pacemaker cell restriction, involving nkx2.5 maintenance in the atrium by nr2f1a. Whether nr2f1 also restricts ventricular fate as shown in the mouse model will require more investigation. This manuscript will be of interest to developmental biologists and pediatric cardiologists.
https://doi.org/10.7554/eLife.77408.sa0Introduction
The vertebrate heart relies on the coordination of multiple specialized cardiomyocyte populations. Cardiomyocytes of the atrium, ventricle, atrioventricular canal (AVC), and pacemaker are each characterized by distinct mechanical, morphological, and electrophysiological properties (Bootman et al., 2006; Brandenburg et al., 2016; Christoffels et al., 2010; Keith and Flack, 1907; Ng et al., 2010; Smyrnias et al., 2010). Pacemaker cardiomyocytes (PCs) of the sinoatrial node (SAN), which is located at the venous pole of the single atrial chamber in zebrafish and the base of the right atrial chamber in mammals, initiate the synchronized wave of contraction that passes through the working atrial cardiomyocytes (ACs), slows at the AVC, and then rapidly activates ventricular cardiomyocytes (VCs) (Arrenberg et al., 2010; Burkhard et al., 2017; Christoffels et al., 2010; Moorman and Christoffels, 2003). Specific gene regulatory programs confer the differentiation and maintenance of these different cardiomyocyte populations within the heart during development (Barth et al., 2005; Ng et al., 2010; Pradhan et al., 2017; Tabibiazar et al., 2003; Targoff et al., 2013; van Weerd and Christoffels, 2016; Xin et al., 2007). Importantly, mutations in genes associated with promoting cardiomyocyte differentiation as well as the maintenance of identity, such as the transcription factors (TFs) NKX2.5, TBX5, and NR2F2, are associated with human congenital heart defects (CHDs), the most common type of congenital malformations found in newborn children (Al Turki et al., 2014; Benjamin et al., 2018; Benson et al., 1999; Cheng et al., 2011; Hoffman and Kaplan, 2002; Loffredo, 2000; Nakamura et al., 2011; Schott et al., 1998; van der Linde et al., 2011). Moreover, structural CHDs are often accompanied by arrhythmias that cause additional complications, morbidity, and death (Bruneau et al., 1999; Ellesøe et al., 2016; Williams and Perry, 2018). In adults, variants in NKX2.5 and TBX5 are also associated with isolated arrhythmias, such as atrial fibrillation, and improper cardiomyocyte gene expression (Benson et al., 1999; Bruneau et al., 1999; Guo et al., 2016; Jhaveri et al., 2018; Ma et al., 2016; Nakashima et al., 2014; Yu et al., 2014). Thus, the acquisition and maintenance of an appropriate number and identity of each cardiomyocyte population within the developing heart is vital for its normal function throughout life.
The Nr2f (Coup-tf) family of nuclear hormone receptor TFs are conserved regulators of atrial development. In humans, mutations in NR2F2 have been associated with multiple types of CHDs, most commonly atrioventricular septal defects (AVSDs) (Al Turki et al., 2014; Nakamura et al., 2011; Qiao et al., 2018; Upadia et al., 2018). Consistent with a role in AC differentiation, NR2F1 and NR2F2 are expressed in ACs of the developing hearts of mice and humans (Al Turki et al., 2014; Li et al., 2016). In vitro studies have shown that both NR2F1 and NR2F2 promote AC differentiation in human embryonic stem cell (ESC)-derived cardiomyocytes (Churko et al., 2018; Devalla et al., 2015). Global Nr2f2 knockout (KO) mice are embryonic lethal with prominent heart defects, including severely hypomorphic atria that lack septa, which is associated with reduced AC differentiation (Pereira et al., 1999). However, at subsequent stages, Nr2f2 appears to maintain atrial identity by directly repressing TFs including Irx4 and Hey2, which promote and maintain the differentiation of VCs. Hence, cardiac-specific conditional KO of Nr2f2 in mice results in ventricularized atria, while its ectopic expression in VCs is sufficient to promote AC identity (Wu et al., 2013). Although Nr2f2 appears to be the primary regulator of AC differentiation and identity in vivo in mammals, we have found that zebrafish Nr2f1a is the functional homolog of mammalian Nr2f2 with respect to heart development (Dohn et al., 2019; Duong et al., 2018). Nr2f1a mutant zebrafish develop smaller atria due to reduced AC differentiation at the venous pole and fail to restrict the size of the AVC, which results in an expansion of AVC markers into the atrium (Duong et al., 2018). Despite the established requirements for Nr2f TFs in promoting and maintaining atrial differentiation, the mechanisms by which Nr2f TFs function within vertebrate ACs are still not well understood.
Interactions between different regulatory programs determine the number and proportions of cardiomyocytes within the ventricle, AVC, atrium, and pacemaker that initially differentiate and are continuously required to reinforce their identity, which maintains compartmental boundaries within the heart that allow it to function properly. At the venous pole of the heart, a conserved regulatory network of TFs including Isl1, Shox/Shox2, and Tbx3 drives PC differentiation and concurrently represses AC identity. Loss of these TFs results in ectopic expression of working AC markers within the SAN and, in turn, a hypoplastic SAN; conversely, ectopic Tbx3 in working murine cardiomyocytes is sufficient to confer PC identity, and overexpression of Shox2 in murine ESCs during differentiation leads to an upregulation of PC genes (Espinoza-Lewis et al., 2009; Espinoza-Lewis et al., 2011; Hoogaars et al., 2007; Ionta et al., 2015; Liang et al., 2015; Liu et al., 2011; Tessadori et al., 2012; van Weerd and Christoffels, 2016). Nkx2.5 performs multiple functions within cardiomyocytes during heart development, consistent with pleiotropic structural and conduction defects found in patients with NKX2.5 mutations (Benson et al., 1999; Ellesøe et al., 2016; Elliott et al., 2003; Jhaveri et al., 2018; McElhinney et al., 2003; Schott et al., 1998; Xie et al., 2013; Yu et al., 2014). At the venous pole of the heart, it is a critical repressor of the PC differentiation regulatory network within working ACs, as loss of Nkx2.5 results in an expansion of PC marker identity within the atrium of mice and zebrafish (Colombo et al., 2018; Espinoza-Lewis et al., 2011; Nakashima et al., 2014). Within VCs at the arterial pole of the heart, Nkx2.5, as well as the TFs Irx4 and Hey2, maintain ventricular identity and repress the atrial gene program. Loss of these genes in mice and zebrafish results in ectopic expression of AC markers within the ventricle (Bao et al., 1999; Bruneau et al., 2000; Koibuchi and Chin, 2007; Pradhan et al., 2017; Targoff et al., 2013). Moreover, in embryonic zebrafish hearts, AC identity is repressed in VCs by FGF signaling, which functions upstream of Nkx2.5 and by inhibition of BMP (de Pater et al., 2012; Pradhan et al., 2017). However, in mice, FGF signaling is also necessary for proper ventricular conduction through regulation of connexin 43 (Sakurai et al., 2013). While the repression of AC identity is central to normal development of the ventricle and SAN, we still lack an understanding of mechanisms that maintain AC identity and plasticity within the embryonic vertebrate heart.
Here, we demonstrate that Nr2f1a is required to inhibit the progressive acquisition of VC and PC identity within different regions of the embryonic zebrafish atrium. At the arterial (outflow) pole of the atrium, Nr2f1a represses VC identity, similar to what has been shown in conditional Nr2f2 KO mice. However, we find it is only sensitized cardiomyocytes within the expanded AVC of nr2f1a mutants, which co-express AC and VC differentiation markers, that progressively resolve to express only VC markers. Importantly, at the venous (inflow) pole of the atrium, we identify a requirement for Nr2f1a in repressing PC identity. While the initial PC population differentiates normally in nr2f1a mutants despite reduced AC differentiation, there is a subsequent expansion of PC markers and complementary loss of Nkx2.5 from the venous pole into the atrium. Electrophysiological studies demonstrate that the transdifferentiated ACs in nr2f1a mutants are adopting central conduction system identity. Genetic epistasis using a heat-shock inducible nkx2.5 transgene showed that overexpressing Nkx2.5 is sufficient to inhibit the transdifferentiation of ACs that occurs in nr2f1a mutant hearts, indicating that within ACs Nr2f1a functions upstream of Nkx2.5 in the conserved genetic hierarchy that limits PC differentiation. Chromatin accessibility integrated with our transcriptomic analysis of sorted ACs identified a putative Nr2f1a-dependent nkx2.5 enhancer that is expressed in a specific subset of venous ACs bordering the SAN. Deletion of this enhancer results in a loss of Nkx2.5 expression and an expansion of a PC reporter within ACs at the venous pole of the atrium. Overall, our results demonstrate that maintenance of Nr2f expression is critical for normal vertebrate heart development through sustaining atrial identity in different AC subpopulations, which may provide insights into the etiology of concomitant congenital structural cardiac and conduction defects found in humans.
Results
Cardiomyocytes within the AVC progressively resolve to VC identity in nr2f1a mutants
Although previous studies have indicated that Nr2f TFs are required to promote atrial differentiation (Duong et al., 2018; Pereira et al., 1999; Wu et al., 2013), we still do not fully understand the effects of Nr2f loss on AC identity. To investigate the consequences of Nr2f loss within the atrium, we performed RNA-seq on flow-sorted cardiomyocytes expressing the AC differentiation marker atrial myosin heavy chain (amhc; also called myh6) from wild-type (WT) and nr2f1a mutant embryos. Embryos carrying the Tg(amhc:EGFP) transgene (Zhang et al., 2013) were dissociated at 48 hr post-fertilization (hpf) and fluorescence-activated cell sorting (FACS) was performed to isolate amhc:EGFP+ cardiomyocytes (Figure 1A). The resulting RNA-seq dataset revealed 2077 upregulated and 1297 downregulated genes with >2-fold change in the nr2f1a mutant amhc:EGFP+ cardiomyocytes relative to WT (Supplementary file 1). Genes associated with AC differentiation, such as amhc itself, were decreased in the amhc:EGFP+ cardiomyocytes captured from nr2f1a mutant embryos relative to the WT embryos, confirming the success of this approach (Figure 1B; Supplementary file 1). However, we also observed two additional trends in the nr2f1a mutant amhc:EGFP+ cardiomyocytes: an increase in the expression of genes associated with VC and outflow tract identity of the heart, such as the irx4 and isl2 paralogs (Bruneau et al., 2000; Bruneau et al., 2001; Targoff et al., 2013; Witzel et al., 2017; Wu et al., 2013); and misexpression of core genes associated with PC differentiation, with genes shown to promote PC differentiation, such as tbx3a and shox, being increased (Hoogaars et al., 2007; Kenyon et al., 2011; Liu et al., 2011; Ribeiro et al., 2007; Sawada et al., 2015), while nkx2.5, a critical repressor of PC differentiation, was decreased (Colombo et al., 2018; Espinoza-Lewis et al., 2011; Nakashima et al., 2014; Figure 1B). Thus, transcriptional analysis implies that nr2f1a mutant atria have increased expression of genes associated with promoting both VC and PC differentiation.

The atrioventricular canal (AVC) resolves to Vmhc-expressing cardiomyocytes in nr2f1a mutant hearts.
(A) Schematic for isolation of atrial cardiomyocytes (ACs) using Tg(amhc:EGFP) transgene for RNA-seq and assay for transpose-accessible chromatin sequencing (ATAC-seq) at 48 hr post-fertilization (hpf). (B) Differential expression of genes associated with ventricular cardiomyocyte (VC)/arterial pole, AC, and pacemaker cardiomyocyte (PC) differentiation in nr2f1a mutants compared to wild-type (WT). (C–H) IHC for Amhc (blue) and Vmhc (red) in nr2f1a and WT hearts at 48, 72, and 96 hpf. White arrows indicate region of overlapping Amhc+ and Vmhc+ cardiomyocytes in the AVC. Blue arrows indicate overt arterial position of the AVC of nr2f1a mutant atria. Number of embryos examined - 48 hpf: WT (n=3), nr2f1a-/- (n=4); 72 hpf: WT (n=4), nr2f1a-/- (n=6); 96 hpf: WT (n=7), nr2f1a-/- (n=10). (I) Schematic for cell quantification: Amhc+/Vmhc- (blue) cardiomyocytes (CMs) mark the atria; Amhc+/Vmhc+ (green/blue) cardiomyocytes mark the AVC; Amhc-/Vmhc+ (green) cardiomyocytes mark the ventricles; DsRed2-NLS (pan-cardiac) marks the cardiomyocyte nuclei. (J–L) Quantification of Amhc+/Vmhc- cardiomyocytes, Amhc+/Vmhc+ (AVC) cardiomyocytes, and Amhc-/Vmhc+ cardiomyocytes. Error bars indicate s.e.m. 48 hpf: WT (n=4), nr2f1a-/- (n=5); 72 hpf: WT (n=5), nr2f1a-/- (n=8); 96 hpf: WT (n=5), nr2f1a-/- (n=12). Scale bar indicates 50 μm. Differences between WT and nr2f1a-/- were analyzed using ANOVA with multiple comparisons. *p=0.05–0.001, **p<0.001.
Conditional KO of Nr2f2 in murine ACs results in their acquisition of VC identity (Wu et al., 2013). Furthermore, while our transcriptomic analysis here and previous analysis of the ventricular differentiation marker ventricular myosin heavy chain (vmhc; also called myh7) (Duong et al., 2018) also suggests there may be an expansion of VC/arterial pole identity into the atrium of nr2f1a mutants, it was still not clear if a fate switch of ACs to VCs occurred within nr2f1a mutant atria at 48 hpf comparable to what was reported in the conditional Nr2f2 KO mice (Wu et al., 2013). To examine cardiomyocyte identities more closely within nr2f1a mutant hearts, we performed immunohistochemistry (IHC) for the cardiomyocyte differentiation markers Amhc and Vmhc (Figure 1C–H) and quantified cardiomyocyte number by using the pan-cardiac Tg(myl7:DsRed2-NLS) transgene (Figure 1l; Mably et al., 2003). Although we found a significant overlap of Amhc and Vmhc in the heart tubes at 30 hpf, which was indistinguishable between WT and nr2f1a mutant hearts (Figure 1—figure supplement 1A–A”, E–E”), the number of Amhc+/Vmhc+ cardiomyocytes, which we classify as cardiomyocytes of the AVC, becomes significantly restricted by 48 hpf in the hearts of WT embryos (Figure 1C and J–L; Figure 1—figure supplement 1B–B”). Subsequently, from 48 to 96 hpf there remained a relatively constant and small number of Amhc+/Vmhc+ cardiomyocytes in the hearts of WT embryos (Figure 1C, E, G and J–L; Figure 1—figure supplement 1B–D”; Figure 1—figure supplement 2B). However, by 48 hpf in nr2f1a mutant hearts, there were still a significant number of cardiomyocytes co-expressing Amhc and Vmhc (Figure 1D and J–L; Figure 1—figure supplement 1F–F’’; Figure 1—figure supplement 2B), consistent with their expanded AVC that we reported previously (Duong et al., 2018). Interestingly, the number of Amhc+/Vmhc+ cardiomyocytes within this region of nr2f1a mutant hearts was progressively reduced through 96 hpf (Figure 1F, H and J–L; Figure 1—figure supplement 1G–G”, H-H”; Figure 1—figure supplement 2B). While Amhc+/Vmhc+ cardiomyocytes remained constant over these stages, there was a modest increase of Amhc+/Vmhc- cardiomyocytes in the atria of WT hearts. However, total Amhc+ cardiomyocytes (Amhc+/Vmhc- and Amhc+/Vmhc+) were found to diminish from 48 to 96 hpf in the nr2f1a mutant hearts, which was primarily due to the decrease in Amhc+/Vmhc+ cardiomyocytes as the number of Amhc+/Vmhc- cardiomyocytes, while dramatically reduced compared to WT hearts, did not change significantly over these stages (Figure 1J; Figure 1—figure supplement 1F–H” and Figure 1—figure supplement 2A and B). Conversely, there were more total Vmhc+ (Amhc-/Vmhc+ and Amhc+/Vmhc+) cardiomyocytes in nr2f1a mutant hearts compared to WT hearts across all timepoints, which primarily reflects the expansion of Vmhc expression within the enlarged AVC of the nr2f1a mutant hearts (Figure 1K; Figure 1—figure supplement 1F–H” and Figure 1—figure supplement 2B and C). Furthermore, the number of Vmhc+-only (Amhc-/Vmhc+) cardiomyocytes progressively increased from 48 to 96 hpf in nr2f1a mutant hearts, complementing the loss of Amhc+/Vmhc+ cardiomyocytes within the AVC (Figure 1K; Figure 1—figure supplement 1F–H” and Figure 1—figure supplement 2C and D). Interestingly, the total number of cardiomyocytes in nr2f1a mutant hearts remained constant from 48 to 96 hpf, stages that showed a progressive increase in the total number of cardiomyocytes in WT hearts (Figure 1L), implying that Nr2f1a loss within the atrium affects overall heart growth. Cell death, as assayed with IHC for active Caspase 3, and proliferation, as assayed with IHC for phosphohistone H3 (pHH3), were not affected at these stages in nr2f1a mutant hearts (Figure 1—figure supplement 3A–D; Figure 1—figure supplement 4A–F). Similarly, EdU pulses at 48 or 72 hpf showed that Amhc+ cardiomyocytes in nr2f1a mutants hearts also do not have reduced proliferation (Figure 1—figure supplement 4G–M). However, Amhc- cardiomyocytes, which are predominantly VCs, did show reduced proliferation (Figure 1—figure supplement 4G–M), suggesting the lack of increase in heart size in nr2f1a mutants primarily reflects a lack of ventricular growth. Thus, our data indicate that changes in the different cardiomyocyte populations of nr2f1a mutants hearts are not due to increased cell death or proliferation. However, the lack of observed heart growth is potentially due to Nr2f1a loss cell non-autonomously affecting proliferation of VCs, which will be investigated in the future. Collectively, these results show that the expansion and acquisition of VC identity within nr2f1a mutant hearts is regionally restricted, with cardiomyocytes comprising the expanded AVC that initially expresses both AC and VC differentiation markers resolving to VC identity.
PC identity progressively expands throughout nr2f1a mutant atria
Our results show that in the absence of Nr2f TFs not all Amhc+ cardiomyocytes within the atrium acquire VC identity. Moreover, our transcriptomic data indicated there was also misexpression of core regulatory genes that control PC differentiation in the nr2f1a mutant Amhc+ cardiomyocytes (Figure 1A and B). Orthologs of genes associated with promoting PC differentiation and function, including tbx3a, shox, and fgf13a (Burkhard and Bakkers, 2018; Colombo et al., 2018; Espinoza-Lewis et al., 2009; Hoogaars et al., 2007; Kenyon et al., 2011; Liu et al., 2011; Poon et al., 2016; Ribeiro et al., 2007; Sawada et al., 2015; Wang et al., 2011; Wang et al., 2017; Yang et al., 2016), were predominantly increased, while nkx2.5, a critical repressor of PC identity in mice and zebrafish (Colombo et al., 2018; Hoogaars et al., 2007; Nakashima et al., 2014), was decreased. To determine if there are effects on the PCs in nr2f1a mutants, we used the SqET33-mi59B enhancer trap line, which for clarity we refer to as Et(fgf13a:EGFP), as it reports fgf13a expression and is expressed in PCs (Poon et al., 2016). At 48 hpf, we found no overt difference in Et(fgf13a:EGFP) expression between the WT and nr2f1a mutant embryos (Figure 2A and B; Figure 2—figure supplement 1A–A’’’, D–D’’’), despite the smaller atria in the nr2f1a mutants. However, while Et(fgf13a:EGFP) expression is normally restricted to the venous pole in WT hearts, we found a progressive expansion of Et(fgf13a:EGFP) expression from the venous pole throughout the Amhc+ cardiomyocytes in the atrium of the nr2f1a mutant hearts through 96 hpf (Figure 2C–F; Figure 2—figure supplement 1B–C’’’, E–F’’’). Quantification of Amhc+ cardiomyocytes co-expressing the Et(fgf13a:EGFP) transgene via incorporation of the Tg(myl7:DsRed2-NLS) transgene confirmed the concurrent increase in fgf13a:EGFP+/Amhc+ cardiomyocytes and decrease in fgf13a:EGFP-/Amhc+ cardiomyocytes within the nr2f1a mutant atria (Figure 2G and H), even as the total number of Amhc+ cardiomyocytes decreases (Figure 2I). Equivalent results were found with IHC for the SAN marker Isl1 (Liang et al., 2015; Sun et al., 2007; Tessadori et al., 2012; Figure 2—figure supplement 2). Additionally, real-time quantitative PCR (RT-qPCR) on flow-sorted cardiomyocytes from 96 hpf embryos using the pan-cardiac Tg(myl7:EGFP) transgene (Huang et al., 2003) showed an increase in PC differentiation marker expression, as well as vmhc expression, and decrease in amhc at these stages in nr2f1a mutant hearts (Figure 2J and K), reflecting the progressive expansion observed with expression of Et(fgf13a:EGFP) and Isl1 with IHC analysis.

Pacemaker cardiomyocyte (PC) identity expands from the venous pole across the atrium in nr2f1a mutant hearts.
(A–F) IHC for Amhc (blue), Vmhc (red), and fgf13a:EGFP (green). White arrows indicate boundaries of Et(fgf13a:EGFP) expression within the atrium. Number of embryos examined - 48 hr post-fertilization (hpf): wild-type (WT) (n=10), nr2f1a-/- (n=12); 72 hpf: WT (n=10), nr2f1a-/- (n=10); 96 hpf: WT (n=7), nr2f1a-/- (n=10). (G–I) Quantification of Amhc+/fgf13a:EGFP-, Amhc+/fgf13a:EGFP+, and total Amhc+ cardiomyocytes (CMs) using the Tg(myl7:DsRed2-NLS) transgene within the hearts of WT and nr2f1a mutants. Error bars indicate s.e.m.; 48 hpf: WT (n=16), nr2f1a-/- (n=11); 72 hpf: WT (n=14), nr2f1a-/- (n=11); 96 hpf: WT (n=15), nr2f1a-/- (n=12). (J) Schematic for the isolation of cardiomyocytes at 96 hpf using the Tg(myl7:EGFP) transgene. (K) Fold change of marker genes relative to β-actin from real-time quantitative PCR (RT-qPCR) on isolated cardiomyocytes of WT and nr2f1a mutants at 96 hpf. Scale bar indicates 50 μm. Differences between WT and nr2f1a-/- were analyzed using ANOVA with multiple comparisons. *p=0.05–0.001, **p<0.001.
Because the total number of cardiomyocytes in the hearts of the nr2f1a mutants does not change over the stages examined (Figure 1L), the aforementioned analysis of cardiomyocyte death and proliferation (Figure 1—figure supplement 3 and Figure 1—figure supplement 4) would suggest that the arterial-directed expansion of Et(fgf13a:EGFP) expression in Amhc+ cardiomyocytes reflects an acquisition of PC identity in these cardiomyocytes and is not due to the incorporation of newly differentiating cardiomyocytes at the venous pole of the heart. We sought to further test this hypothesis via lineage tracing of cardiomyocytes using a temporal differentiation assay with the transgene Tg(myl7:NLS-KikGR), which expresses the photoconvertible NLS-KikGR in all cardiomyocytes (Lazic and Scott, 2011). The cardiomyocyte-specific NLS-KikGR was photoconverted from green to red in the hearts of embryos also carrying the Et(fgf13a:EGFP) transgene at 72 hpf. The hearts were then analyzed at 96 hpf for the presence of red and green fgf13a:EGFP+ cardiomyocytes (Figure 2—figure supplement 3A). Importantly, the green, nuclear KikGR was visible within the cardiomyocytes even with expression of the cytosolic fgf13a:EGFP. We found that in both WT and nr2f1a mutant hearts all fgf13a:EGFP+ cardiomyocytes within the atria co-expressed both photoconverted (red) NLS-KikGR and newly expressed (green) NLS-KikGR, including those in the expanded fgf13a:EGFP+ expression domain and venous pole of nr2f1a mutant hearts (Figure 2—figure supplement 3B–E”). The absence of green-only, newly differentiated NLS-KikGR+ cardiomyocytes in the nr2f1a mutant atria at 96 hpf indicates that all of the fgf13a:EGFP+ cardiomyocytes were present in the atrium at 72 hpf. Thus, our data support that venous Amhc+ ACs are progressively transdifferentiating into PCs in the absence of Nr2f1a.
Next, we wanted to determine if the ACs that are acquiring PC marker gene expression in the nr2f1a mutants are functionally adopting PC characteristics. PCs have different electrophysiological properties than working cardiomyocytes, including slower depolarization and the ability to spontaneously generate action potentials (Christoffels et al., 2010; Lakatta et al., 2010; Mesirca et al., 2015; Schram et al., 2002). Moreover, given that the intrinsic automaticity of coupled cells is a function of their total electrical capacitance (Honjo et al., 1996; Joyner and van Capelle, 1986; Mangoni and Nargeot, 2008; Watanabe et al., 1995), one might predict that if there are more ACs adopting PC or central conduction system-like characteristics then the heart rates of nr2f1a mutants may be slower than WT embryos. We assessed heart rate using high-speed time lapse imaging (Figure 3A and B; Figure 3—videos 1–6). At 48 hpf, there was no difference in heart rate between WT and nr2f1a mutant hearts (Figure 3C). However, as the heart rate of WT embryos increased through 96 hpf, the heart rate of nr2f1a mutant hearts failed to increase (Figure 3C). To examine the electrophysiological properties of the ACs that are acquiring PC marker expression in nr2f1a mutant hearts, we analyzed myocardial electrophysiology using optical voltage mapping (Colombo et al., 2018; Mosimann et al., 2015; Panáková et al., 2010). The intrinsic bradycardia we observed appeared to be related to effects on the duration and extent of repolarization with evidence of modest hyperpolarization in both heterozygous and homozygous nr2f1a mutant hearts, as well as substantial prolongation of repolarization in the homozygous nr2f1a mutant hearts, as conservatively assessed based on the action potential duration at 20% repolarization (APD20) (Figure 3D). Measuring the APD20 within the center of the atria (Figure 3A), which acquires PC marker expression in the nr2f1a mutants based on the IHC analysis, showed that even at 48 hpf the ACs in nr2f1a mutant hearts depolarize more slowly compared to WT and heterozygous nr2f1a atria (Figure 3E; Figure 3—figure supplement 1A). These effects were not associated with any changes in Phase 4 depolarization (data not shown). There was evidence of significantly slower conduction across the nr2f1a mutant atria by 96 hpf (Figure 3F–L; Figure 3—figure supplement 1B), but with no effect on Vmax (Figure 3—figure supplement 1C), consistent with less chamber-like and more central conduction system-like intercellular coupling characteristics. Remarkably, there was no evidence of the emergence of the AVC slow conduction zone in the nr2f1a mutant hearts that is characteristic of the AV junction and the resultant sequential chamber contraction that is the hallmark of the vertebrate heart (Figure 3G–L), supporting that the cardiomyocytes in the AVC of nr2f1a mutants also have acquired intermediate physiological properties with the VCs. Collectively, our data support that more venous ACs in nr2f1a mutant hearts are adopting a functional identity that exhibits delayed repolarization and reduced intercellular coupling, closer to the physiology of the cardiomyocytes of the central conduction system and paralleling the PC-associated gene expression changes we observe within the nr2f1a mutant atria.

Atrial cardiomyocytes (ACs) in nr2f1a mutants function as pacemaker cardiomyocytes (PCs).
(A) Schematic of point placements to measure heart rate (cyan) and action potential duration at 20% repolarization (APD20) (magenta) in a 48 hr post-fertilization (hpf) heart. (B) Representative kymograph used to analyze heart rate in wild-type (WT) and nr2f1a mutant embryos. (C) Quantification of heart rate in WT and nr2f1a mutant embryos; 48 hpf: WT (n=6), nr2f1a-/- (n=6); 72 hpf: WT (n=6), nr2f1a-/- (n=6); 96 hpf: WT (n=6), nr2f1a-/- (n=6). (D) Representative atrial action potentials of WT, nr2f1a+/-, and nr2f1a-/- embryos at 96 hpf. (E) APD20 in WT and nr2f1a-/- hearts; 48hpf: WT (n=21), nr2f1a-/- (n=4); 72 hpf: WT (n=13), nr2f1a-/- (n=9); 96 hpf: WT (n=9), nr2f1a-/- (n=14). (F) Atrial velocity in WT and nr2f1a-/- hearts; 72 hpf: WT (n=13), nr2f1a-/- (n=9); 96 hpf: WT (n=9), nr2f1a-/- (n=14). (G–L) Representative isochronal maps illustrating the positions of the depolarizing wave front in 5 ms intervals for WT and nr2f1a-/- embryos at 48, 72, and 96 hpf; 48 hpf: WT (n=20), nr2f1a-/- (n=6) 72 hpf: WT (n=12), nr2f1a-/- (n=7) 96 hpf: WT (n=11), nr2f1a-/- (n=14). Differences between WT and nr2f1a-/- analyzed using ANOVA with multiple comparisons. Error bars in all graphs indicate s.e.m. *p=0.05–0.001, **p<0.001.
Nr2f1a represses SAN identity upstream of Nkx2.5
We next wanted to understand how Nr2f1a may be regulating the PC gene regulatory network within the atrium. We prioritized examining Nkx2.5 in nr2f1a mutant hearts, as Nkx2.5 loss in mice and zebrafish results in an expansion of PCs within the atrium (Colombo et al., 2018; Espinoza-Lewis et al., 2011; Nakashima et al., 2014) and its expression was reduced in the bulk RNA-seq data of nr2f1a mutant amhc:EGFP+ cardiomyocytes (Figure 1B). IHC for Nkx2.5 and either the Et(fgf13a:EGFP) reporter or Isl1 both showed that in WT hearts at 48 through 96 hpf Nkx2.5 is expressed throughout the atrial myocardium but is predominantly excluded from the PCs at the venous pole (Figure 4A–C”; Figure 4—figure supplement 1), consistent with what has been reported in zebrafish and mice (Colombo et al., 2018; Mommersteeg et al., 2007). At 48 hpf in nr2f1a mutants, Nkx2.5 shows a similar complementary pattern with the Et(fgf13a:EGFP) reporter despite the smaller atria (Figure 4D). As the Et(fgf13a:EGFP) reporter expands toward the AVC into the Amhc+ ACs through 96 hpf in nr2f1a mutants, Nkx2.5 expression recedes from the venous pole and predominantly maintains the complementary expression pattern (Figure 4E–F”). Equivalent results were obtained using the TgBAC(–36Nkx2.5:ZsYellow) transgene (Zhou et al., 2011; Figure 4—figure supplement 2). Thus, a progressive loss of Nkx2.5 within ACs correlates with and abuts the expansion of PC identity in nr2f1a mutant atria.

Nkx2.5 expression recedes from venous pole in nr2f1a mutant atria.
(A–F’’’) IHC for Nkx2.5 (purple), Amhc (blue), and fgf13a:EGFP (SAN - green) in wild-type (WT) and nr2f1a mutant embryos from 48 to 96 hr post-fertilization (hpf). White arrows indicate border of Nkx2.5+ and fgf13a:EGFP+ cardiomyocytes. Number of embryos examined - 48 hpf: WT (n=7), nr2f1a-/- (n=6); 72 hpf: WT (n=7), nr2f1a-/- (n=10); 96 hpf: WT (n=14), nr2f1a-/- (n=23). Scale bar indicates 50 μm.
To determine if restoring Nkx2.5 in nr2f1a mutants can repress the transdifferentiation of ACs to PCs, we used the Tg(hsp70l:nkx2.5-EGFP) transgene (George et al., 2015). Nr2f1a+/- fish that carry the Tg(myl7:DsRed2-NLS) and Et(fgf13a:EGFP) transgenes were crossed to nr2f1a+/- fish carrying the Tg(hsp70l:nkx2.5-EGFP) transgene (Figure 5A). The resulting embryos were sorted for the Et(fgf13a:EGFP) transgene and then heat-shocked at 20 hpf, as a single heat-shock at ~20 hpf using the Tg(hsp70l:nkx2.5-EGFP) transgene is sufficient to rescue nkx2.5 mutants to adulthood (George et al., 2015). Following heat-shock, the embryos were sorted based on the presence of Nkx2.5-EGFP expression and allowed to develop until 48 and 96 hpf, at which point the Nkx2.5-EGFP is no longer visible, and the number of fgf13a:EGFP+ cardiomyocytes was quantified (Figure 5A). At 48 hpf, there was a modest decrease in the number of fgf13a:EGFP+ cardiomyocytes in Nkx2.5-EGFP+ WT hearts compared to the Nkx2.5-EGFP- WT siblings (Figure 5B, C, J and K). However, Nkx2.5-EGFP+ nr2f1a mutants had a complete absence of fgf13a:EGFP+ cardiomyocytes at the venous pole (Figure 5F, G, J and K), suggesting that Nkx2.5 is sufficient to repress PC differentiation in nr2f1a mutant hearts and, further, that nr2f1a mutant ACs may be sensitized to increased Nkx2.5 expression. At 96 hpf, both Nkx2.5-EGFP+ WT and nr2f1a mutant hearts had an increased number of fgf13a:EGFP+ cardiomyocytes compared to those at 48 hpf, yet Nkx2.5-EGFP+ nr2f1a mutants still had significantly fewer fgf13a:EGFP+ cardiomyocytes than mutants lacking the transgene (Figure 5D, E and H–K). While we see a repression of the fgf13a:EGFP+ cardiomyocyte expansion at 48 hpf, we found that the overexpression of Nkx2.5 did not rescue the heart rate of nr2f1a mutants at 96 hpf when we observe the resumption of the PC reporter (Figure 5—figure supplement 1), suggesting that despite the inhibition of the reporter expansion, additional Nkx2.5-independent mechanisms downstream of Nr2f1a also contribute to the inhibition of PC identity. We found that heat-shock induction of Nkx2.5-EGFP at 40 hpf could also repress the expansion of fgf13a:EGFP+ cardiomyocytes in nr2f1a mutant atria (Figure 5—figure supplement 2A–C, F, G, J and K). However, the repression of fgf13a:EGFP+ cardiomyocytes within the mutant atria was not as effective or prolonged as the heat-shock induction of Nkx2.5-EGFP at 20 hpf (Figure 5—figure supplement 2D, E and H–K). Overall, these data suggest that Nr2f1a prevents the expansion of PC identity within the venous pole of the atrium at least in part through promoting or maintaining Nkx2.5 in ACs.

Nkx2.5 induction at 20 hr post-fertilization (hpf) represses pacemaker cardiomyocyte (PC) marker expansion in nr2f1a mutant atria.
(A) Schematic of timeline for heat-shock experiments. (B–I) IHC of representative hearts for Amhc (blue), fgf13a:EGFP (green), and myl7:DsRed2-NLS (red) used to quantify Amhc+/fgf13a:EGFP+ cardiomyocytes (CMs) marking PCs from wild-type (WT) and nr2f1a mutant embryos with and without Nkx2.5-EGFP. White arrows indicate boundaries of Tg(fgf13a:EGFP) expression. Scale bar indicates 50 μm. (J–K) Quantification of Amhc+/fgf13a:EGFP+/myl7:DsRed2-NLS+ (PCs) and Amhc+/fgf13a:EGFP-/myl7:DsRed2-NLS+ cardiomyocytes in the atria with and without Nkx2.5-EGFP induction; 48 hpf WT: Nkx2.5-EGFP- (n=6), Nkx2.5-EGFP+ (n=4); 48 hpf nr2f1a-/-: Nkx2.5-EGFP- (n=7), Nkx2.5-EGFP+ (n=5); 96 hpf WT: Nkx2.5-EGFP- (n=4), Nkx2.5-EGFP+ (n=4); 96 hpf nr2f1a-/-: Nkx2.5-EGFP- (n=7), Nkx2.5-EGFP+ (n=4). Differences between WT and nr2f1a-/- at the different time points were analyzed using Student’s t-test. Error bars indicate s.e.m. *p=0.05–0.001, **p<0.001.
A putative nkx2.5 enhancer is expressed within the venous pole of the atrium
To identify cis-regulatory elements that may control nkx2.5 expression and be affected by Nr2f1a loss, we assayed chromatin accessibility by performing assay for transpose-accessible chromatin sequencing (ATAC-seq) on flow-sorted amhc:EGFP+ cardiomyocytes from 48 hpf embryos, similar to what is described above (Figure 1A; Figure 6—figure supplement 1). Changes in called peaks of WT and nr2f1a mutant amhc:EGFP+ cardiomyocytes were scored relative to the nearest differentially expressed genes (Supplementary file 1). HOMER was used to determine enrichment for canonical Nr2f TF binding sites on the putative enhancer elements. One of the peaks that was revealed to close in nr2f1a mutant atria and contain an Nr2f binding site was ~55 kb upstream of the nkx2.5 transcriptional start site (Figure 6A). Although this region did not show significant sequence conservation, enhancer conservation can vary between tissue types, with enhancers that drive expression in the heart often being weakly conserved (Blow et al., 2010). This putative enhancer region was placed into a reporter vector upstream of the basal E1b promoter to drive GFP expression (Figure 6B; Birnbaum et al., 2012), and transgenic lines were generated using Tol2 (Kwan et al., 2007). Remarkably, stable transgenic lines showed that this putative enhancer is expressed in a group of cardiomyocytes at the venous pole of the atrium (Figure 6C), which are directly adjacent to Isl1+ PCs (Figure 6D and D’). The Tg(–55nkx2.5:EGFP) transgene was then crossed into nr2f1a mutant allele carriers to determine if the expression of this putative nkx2.5 enhancer is affected in nr2f1a mutants. We found the number of atria with Tg(–55nkx2.5:EGFP) expression was significantly decreased in nr2f1a mutants compared to WT, suggesting that Nr2f1a is required for enhancer expression or for those cardiomyocytes to develop (Figure 6C, E and F). Together, these data suggest that a putative Nr2f1a-dependent nkx2.5 enhancer promotes expression in ACs adjacent to PCs.

A putative nkx2.5 enhancer is expressed in atrial cardiomyocytes (ACs) adjacent to the sinoatrial node (SAN).
(A) Comparison of assay for transpose-accessible chromatin sequencing (ATAC-seq) peaks from wild-type (WT) and nr2f1a mutant amhc:EGFP+ cardiomyocytes at 48 hr post-fertilization (hpf), ~55 kb upstream of nkx2.5. Asterisk denotes canonical Nr2f binding site. (B) Schematic of –55 kb nkx2.5 enhancer reporter construct. (C) IHC for Tg(–55nkx2.5:EGFP) enhancer (green), Amhc (blue), and Vmhc (red) in hearts at 72 hpf. White arrow indicates –55nkx2.5:EGFP expression. (D–D’) IHC for –55nkx2.5:EGFP (green), Isl1 (purple), Amhc (blue), and sarcomeric myosin heavy chain (Mhc: pan-cardiac - red). The transgenic enhancer reporter is expressed in the ACs directly adjacent to the Isl1+ pacemaker cardiomyocytes (PCs). (E) IHC for Tg(–55nkx2.5:EGFP) enhancer (green), Amhc (blue), and Vmhc (red) in an nr2f1a mutant heart at 72 hpf. Expression of transgenic enhancer reporter is lost in nr2f1a mutants. Yellow arrow indicates venous pole devoid of –55nkx2.5:EGFP expression. (F) Comparison of percentage of WT and nr2f1a mutant embryos with –55nkx2.5:EGFP expression in the heart at 72 hpf WT (n=16); nr2f1a-/- (n=16), using Fisher’s exact test. (G,H) IHC for Nkx2.5 (green), Amhc (purple), and Mhc (red) in uninjected control and –55nkx2.5 crispant embryo hearts at 96 hpf. White arrow in (G) indicates the border of Nkx2.5+ cardiomyocytes extends close to the venous pole. Yellow arrow in (H) indicates that the Nkx2.5+ cardiomyocyte border is located farther from the venous pole of the atrium in crispant embryos. (I) Quantification of Nkx2.5+/Amhc+ cardiomyocytes in uninjected control (n=11) and –55nkx2.5 crispant (n=11) embryos at 96 hpf. (J,K) IHC for fgf13a:EGFP (green), Amhc (blue), and myl7:DsRed2-NLS (red) in uninjected control and –55nkx2.5 crispant Tg(fgf13a:EGFP)+ embryo hearts at 96 hpf. White arrows in (J) indicate the region of fgf13a:EGFP+ cardiomyocytes at the venous pole. Yellow arrows in (K) indicates that the region of fgf13a:EGFP+ cardiomyocytes is expanded from the venous pole in the atrium of crispant embryos. (L) Quantification of fgf13a:EGFP+/Amhc+ cardiomyocytes in uninjected control (n=10) and –55nkx2.5 crispant (n=11) embryo hearts at 96 hpf. Differences between uninjected control and –55nkx2.5 crispants were analyzed using Student’s t-test. Error bars indicate s.e.m. *p=0.05–0.001, **p<0.001. Scale bars in all images indicate 50 μm.
To determine if the putative nkx2.5 enhancer is required to promote Nkx2.5 expression and restrict PC identity within the venous pole of the atrium, we employed an established, highly efficient CRISPR/Cas9 system to remove the enhancer in embryos with guide RNAs (gRNAs) that flank the nkx2.5 enhancer sequence (Hoshijima et al., 2019). We found that Tg(–55nkx2.5:EGFP) crispant (CRISPR/Cas9-injected) embryos completely lacked or had significantly reduced EGFP expression within their hearts and other tissues where the enhancer is expressed (Figure 6—figure supplement 2A–C). PCR for the transgenic locus in individual Tg(–55nkx2.5:EGFP) embryos and the endogenous locus in individual WT embryos showed the gRNA pair used significantly abrogated both the transgenic and WT enhancer loci, respectively (Figure 6—figure supplement 2D and E). IHC for Nkx2.5 in –55nkx2.5 crispants showed that Nkx2.5 did not extend as far toward to the venous pole of their atria and that they had a significant reduction in the number of Nkx2.5+/Amhc+ cardiomyocytes within the atria compared to control embryos (Figure 6G–I). Conversely, –55nkx2.5 crispant embryos concurrently carrying the Et(fgf13a:EGFP) and Tg(myl7:DsRed2-NLS) transgenes showed an arterial-directed expansion of Et(fgf13a:EGFP) expression within their atria and an increase in the number of fgf13a:EGFP+ cardiomyocytes at the venous pole compared to control embryos (Figure 6J–L). Altogether, these data support that this putative nkx2.5 enhancer is required to promote or maintain Nkx2.5 expression within a subpopulation of ACs that in turn may limit the extent of PCs within the venous pole of the atrium.
Discussion
Vertebrate embryonic cardiomyocytes maintain a high degree of plasticity, with their identity being continually reinforced even after they have overtly differentiated (Barth et al., 2005; Ng et al., 2010; Pradhan et al., 2017; Tabibiazar et al., 2003; Targoff et al., 2013; van Weerd and Christoffels, 2016; Xin et al., 2007). Here, our data have advanced our understanding of the requirements of Nr2f TFs within the vertebrate heart and show that they are required to maintain AC identity while concurrently repressing both VC and PC identity in different regions of the atrium (Figure 7). In mice, Nr2f2 maintains AC identity by repressing the expression of TFs that promote a ventricular identity program, including Irx4 and Hey2 (Wu et al., 2013). Moreover, murine Nr2f2 appears to be required over a relatively long developmental period after ACs express differentiation markers to prevent the acquisition of VC identity within the atrium. Specifically, loss of Nr2f2 by e9.5 using the Myh6-Cre or by e12.5 using the inducible Myh6-MerCreMer both result in the ectopic expression of VC differentiation genes within ACs (Wu et al., 2013). Conversely, ectopic Nr2f2 expression in VCs is sufficient to induce AC gene expression (Wu et al., 2013). Although our previous work had shown that zebrafish Nr2f1a promotes atrial differentiation (Duong et al., 2018), it was still unclear if ACs in zebrafish nr2f1a mutants acquired VC identity, similar to what was shown for Nr2f2 loss in mice. We now show that there is acquisition of VC identity within the nr2f1a mutant atria. However, with the murine conditional Nr2f2 KO approaches, the two atria were formed, but enlarged and appeared to express VC markers throughout the atria (Wu et al., 2013). In contrast to what was shown in the conditional KO mice, the acquisition of VC identity in zebrafish nr2f1a mutants was limited to the population of less-differentiated cardiomyocytes that comprise the expanded AVC, which extends into the smaller single atria. Moreover, our data show that initially within the early zebrafish heart tube there is significant overlap of Vmhc and Amhc expression. Normally, the population of cardiomyocytes expressing both Vmhc and Amhc becomes significantly refined to the AVC. However, this population of cardiomyocytes that comprise the AVC is not restricted within the nr2f1a mutants and they then progressively resolve to VC identity in an arterial-to-venous direction via losing expression of AC differentiation markers. We presently cannot rule out that these overt differences in effects on the ACs are due to slightly different species-specific requirements of Nr2f TFs within mice and zebrafish or technical differences from the global zebrafish nr2f1a mutants vs the conditional Nr2f2 approaches. Nevertheless, taken together, the trends reinforce that there are conserved requirements for vertebrate Nr2f TFs in repressing VC identity within ACs. However, our results with global loss of Nr2f1a support that not all ACs have the potential to transdifferentiate into VCs and that this transformation occurs in a sensitized AC population.

Model summarizing the consequences of Nr2f1a loss on atrial development.
Following its initial requirement promoting atrial cardiomyocyte (AC) differentiation (Duong et al., 2018), Nr2f1a is required to maintain AC identity within the arterial (outflow) and venous (inflow) regions of the zebrafish atrium. In the absence of Nr2f1a, cardiomyocytes within the enlarged atrioventricular canal (AVC), where there is expression of both AC and ventricular cardiomyocyte (VC) differentiation markers, progressively lose expression of AC differentiation markers and only express VC differentiation markers. Concurrently, more venous ACs progressively gain PC differentiation marker expression.
Our data also support a new paradigm whereby Nr2f TFs are required to limit the acquisition of PC identity within more venous ACs. Work predominantly in mice has revealed a transcriptional regulatory network that promotes the differentiation of PCs within the SAN at the venous pole of the right atrium (Burkhard et al., 2017; Christoffels et al., 2010; van Weerd and Christoffels, 2016). Tbx3, Isl1, and Shox/Shox2 TFs are part of a core regulatory network that promotes PC differentiation within venous ACs, while Nkx2.5 represses PC identity in working ACs (Blaschke et al., 2007; Espinoza-Lewis et al., 2009; Espinoza-Lewis et al., 2011; Hoogaars et al., 2007; Liang et al., 2015; Liu et al., 2011; Mommersteeg et al., 2007; Nakashima et al., 2014; Sun et al., 2007; Weinberger et al., 2012; Wiese et al., 2009). Although the characterization of PC differentiation has been less extensive in zebrafish, functional analysis of Isl1, Shox/Shox2, and Nkx2.5 support that this regulatory network is fundamentally conserved (Blaschke et al., 2007; Colombo et al., 2018; de Pater et al., 2009; Hoffmann et al., 2013; Tessadori et al., 2012). Our data show that in nr2f1a mutants there is a progressive expansion of central conduction system identity from the venous pole toward the arterial pole of the atrium, supporting a role for Nr2f TFs in regulating this conserved network within ACs. The progressive expansion, where the number of PCs is initially unchanged in nr2f1a mutant atria despite being smaller, contrasts some with Nkx2.5 mutant mice and zebrafish (Colombo et al., 2018; Mommersteeg et al., 2007; Nakashima et al., 2014), which show an initial expansion of PC markers and ectopic differentiation at earlier stages. Although the studies of Nr2f2 conditional KO mice did not examine PC and SAN differentiation, manual interrogation of their data, which was from whole hearts, does show they have increased expression of Tbx3 and Fgf13 (Wu et al., 2013), implying that this requirement in repressing PC identity may also be conserved within the vertebrate atrium. However, determining if murine Nr2f2 possesses similar requirements in repressing PC identity at the venous pole of the right atrium in mammals will require additional analysis.
Our data support that one mechanism by which Nr2f1a intersects with the PC regulatory network and represses the progressive acquisition of PC identity within ACs is through maintaining Nkx2.5 expression at the venous pole of the atrium. Interestingly, Nkx2.5 induction at early stages of heart development can initially completely repress PC reporter expression in the nr2f1a mutant atria at later stages. Thus, these results imply that the transcriptional complexes reflecting the complementary changes in Nkx2.5 and PC gene expression observed at later stages must actually be established much earlier during heart development and that Nr2f1a loss actually sensitizes these ACs to Nkx2.5. Furthermore, early Nkx2.5 induction does not produce a permanent repression of PC reporter expansion, as the PC reporter expands by 96 hpf, though less than in nr2f1a mutant atria. While one interpretation could be that Nkx2.5 expression needs to be maintained within the ACs over a longer period than provided by the heat-shock to prevent PC expansion, as the Nkx2.5-GFP is only thought to last ~8 hr (George et al., 2015), a long requirement would seem unlikely as nkx2.5 mutants can also be rescued to adulthood following a single heat-shock at the 20 somite stage. Therefore, we postulate that Nr2f1a must also concurrently affect the AC expression of other TFs within the core PC regulatory network that promote PC differentiation, which warrants further investigation. Furthermore, it is interesting that the putative Tg(–55nkx2.5:EGFP) enhancer begins expression at later stages (~72 hpf), correlating with the PC marker expansion and implying the enhancer may be activated to reinforce nkx2.5 expression at this border. Additionally, Nr2f1a appears to be expressed equivalently throughout the cardiomyocytes of the zebrafish atrium (Duong et al., 2018; Figure 6—figure supplement 3) and we did not observe overlap of VC and PC markers within the nr2f1a mutant atria at the stages examined, suggesting that Nr2f1a at least functionally, if not directly, must interact with other signaling pathways and factors to confer the regional effects within these AC populations. Altogether, our current data support that Nr2f1a must intersect and functionally interact with multiple gene regulatory networks to mediate its regional requirements preventing the progressive acquisition of VC and PC identity in different AC subpopulations. The nature of these factors and signaling will be investigated in the future.
In considering the implications of our data, it is interesting that in humans, congenital malformations, including AVSDs, are associated with increased incidence of arrhythmias (Williams and Perry, 2018). Moreover, the arrhythmias can arise from complications secondary to the structural defects or pleiotropic requirements in the development of the myocardium and the conduction system (Bruneau et al., 1999; Ellesøe et al., 2016; Williams and Perry, 2018). A prime example of the latter is mutations in NKX2.5, which are associated with numerous CHDs, including atrial and ventricular septal defects, and conduction defects, consistent with its requirements in working cardiomyocyte and PC differentiation shown in experimental models (Benson et al., 1999; Ellesøe et al., 2016; Elliott et al., 2003; Jhaveri et al., 2018; McElhinney et al., 2003; Schott et al., 1998; Xie et al., 2013; Yu et al., 2014). Similar to NKX2.5, NR2F2 mutations in humans have now been associated with myriad CHDs, including atrial septal defects and AVSDs (Al Turki et al., 2014; Nakamura et al., 2011; Poot et al., 2007; Qiao et al., 2018; Upadia et al., 2018). While it is presently unclear if there are also accompanying arrhythmias or conduction defects in these patients, given our results, we speculate that further evaluation of patients with NR2F2 mutations will show that they have or are at higher risk to develop arrhythmias.
In conclusion, mutations in genes, such as NR2F2, that are required for the maintenance of cardiomyocyte identity have been associated with multiple types of CHDs in humans (Al Turki et al., 2014; Benjamin et al., 2018; Benson et al., 1999; Cheng et al., 2011; Hoffman and Kaplan, 2002; Loffredo, 2000; Nakamura et al., 2011; Schott et al., 1998; van der Linde et al., 2011). Hence, elucidating the transcriptional mechanisms that maintain cardiomyocyte identity and regulate cardiomyocyte plasticity may provide insights into the etiology of CHDs, as well as help refine stem cell-derived tissue engineering and regenerative strategies. Our results have delineated novel requirements for Nr2f TFs in maintaining AC identity in vertebrates, which may help provide a framework for our understanding of the molecular etiology of concurrent structural malformations and conduction defects in patients.
Materials and methods
Zebrafish lines used
Request a detailed protocolAdult zebrafish were raised and maintained under standard laboratory conditions (Westerfield, 2000). Transgenic lines used: Tg(–5.1myl7:DsRed2-NLS)f2 (Mably et al., 2003), SqET33-mi59B (Poon et al., 2016), TgBAC(–36nkx2.5:ZsYellow)fb7 (Zhou et al., 2011), Tg(myh6:EGFP)s958 (Zhang et al., 2013), Tg(hsp70l:nkx2.5-EGFP)fcu1 (George et al., 2015), Tg(myl7:EGFP)f1 (Huang et al., 2003), Tg(myl7:NLS-KikGR) (Lazic and Scott, 2011). Three nr2f1a mutant alleles were used: nr2f1ael512 and nr2f1aci1009, which were reported previously (Duong et al., 2018), and nr2f1aci1017, which was identified from an ENU screen. Additional details of this allele will be reported elsewhere (manuscript in preparation). The nr2f1aci1017 allele was used for RNA-seq and ATAC-seq. The nr2f1ael512 and nr2f1aci1009 alleles were used for the photoconversion assay. All other experiments were completed using the nr2f1aci1009 allele. WT and nr2f1a+/-embryos were overtly indistinguishable, including with respect to their hearts. Therefore, unless otherwise indicated, WT in figures in which embryos were genotyped refers to either WT or nr2f1a+/-embryos. For all experiments, embryos were grown in 0.003% 1-phenyl-2-thiourea to prevent pigmentation.
Generation of transgenic Tg(-55nkx2.5:EGFP) enhancer line
Request a detailed protocolThe putative nkx2.5 enhancer transgenic line was generated using standard Tol2/Gateway methods (Kwan et al., 2007). The enhancer region at ~55 kb upstream of nkx2.5 was amplified with PCR of genomic DNA and cloned into the pDONR221 middle entry vector. Subsequently, it was transferred into the E1b-GFP-Tol2 Gateway destination vector (Addgene plasmid #37846, Birnbaum et al., 2012). To generate transgenic lines, one cell embryos were co-injected with 25 pg of the –55nkx2.5:EGFP plasmid and 25 pg Tol2 mRNA (Kwan et al., 2007). Injected embryos were raised to adulthood and transgenic founders were identified via outcrossing to WT fish. Multiple founders with overtly equivalent GFP expression within the atrium were identified. The line that was used in this study has been designated Tg(–55nkx2.5:EGFP)ci1015.
Generation of –55nkx2.5 enhancer crispants
Request a detailed protocolTo generate –55nkx2.5 enhancer crispants (CRISPR/Cas9-injected embryos), crRNAs flanking the endogenous enhancer sequence were designed using IDT’s custom gRNA design tool (https://www.idtdna.com/). crRNA:tracrRNA duplexes were prepared as described previously (Hoshijima et al., 2019). Briefly, equal volumes of 100 μM target-specific crRNA and common tracrRNA were mixed together and annealed by heating to 95°C for 5 min, then cooled at 0.1 °C/s to 25°C, followed by rapid cooling to 4°C in a PCR machine to make a 50 μM solution of the crRNA:tracrRNA duplex, referred to as the duplex guide RNA (dgRNA). An equal volume of duplex buffer (IDT) was then added to further dilute the dgRNA to 25 μM. To generate 5 μM dgRNA;Cas9 complexes, we combined: 1 μL of each 25 μM dgRNA, 0.4 μL 61 μM Cas9 stock (Alt-R S.p. Cas9 Nuclease V3, IDT), 1.6 μL H2O, and 1 μL phenol red. Prior to injection, the solution was incubated at 37°C for 5 min and then placed at room temperature. We injected 1 nL of the 5 μM dgRNA:Cas9 complex into one-cell stage embryos. Efficiency of dgRNA:Cas9 complexes was assessed using PCR on individual embryos.
IHC and cell quantification
Request a detailed protocolIHC was performed as previously described (Waxman et al., 2008). Embryos were fixed in 1% formaldehyde in PBS for 1 hr at room temperature and then washed 1× in PBS, 2× in 0.1% saponin/1× PBS followed by blocking in saponin blocking solution (0.1% saponin, 10% sheep serum, 1× PBS, 2 mg/mL BSA) for 1 hr. Primary antibodies were diluted in saponin blocking solution and applied to embryos overnight at 4°C. Embryos were washed with 0.1% saponin/1× PBS multiple times before being incubated with secondary antibodies diluted in saponin blocking solution for 2 hr at room temperature. Embryos were washed multiple times with 0.1% saponin/1× PBS before imaging. Antibody information is provided in Supplementary file 2.
IHC for quantification of cardiomyocytes in the AVC required the use of two rabbit primary antibodies (anti-DsRed and anti-Vmhc; Song et al., 2019). To accommodate this, after blocking, embryos were first incubated with anti-Amhc (anti-Myh6/S46, mouse IgG1) and anti-DsRed (rabbit) overnight at 4°C. Secondary antibodies (anti-mouse IgG1-Dylight 405 and anti-rabbit-TRITC) were then applied for 2 hr at room temperature. After washing multiple times with 0.1% saponin/1× PBS, embryos were again incubated in saponin blocking solution for 1 hr. Anti-Vmhc (rabbit) was then diluted in saponin blocking solution and applied to embryos overnight at 4°C. Embryos were then washed multiple times with 0.1% saponin/1× PBS. Anti-rabbit-Alexa-488 was diluted in saponin blocking solution and embryos were incubated in secondary antibody for 1 hr at room temperature, followed by multiple washes with 0.1% saponin/1× PBS.
IHC for activated Caspase 3 was performed as previously described (Rydeen and Waxman, 2016). Embryos were fixed in 4% paraformaldehyde overnight at 4°C then dehydrated in a methanol series and left in 100% methanol overnight at –20°C. Embryos were then rehydrated in PBST (PBS +0.1% Tween-20), permeabilized for 20 min using PDT (PBST + 1% DMSO) supplemented with 0.3% Triton-X, and then incubated in blocking solution (1× PBS, 10% sheep serum, 0.1% Tween-20) for 1 hr. Primary antibodies were diluted in blocking solution and applied to embryos overnight at 4°C. Embryos were then washed 3× for 20 min with PDT and blocked again in blocking solution for 30 min. Secondary antibodies were diluted in blocking solution and applied to embryos for 2 hr at room temperature. Following secondary antibody incubation, embryos were washed multiple times with PDT.
Following IHC, embryos were imaged on a Nikon A1R inverted confocal microscope. Nikon’s Denoise-AI was employed on images taken using an HD resonance scanner. Cardiomyocytes were counted using Photoshop and ImageJ. For quantification of PCs, DsRed2-NLS+/Amhc+/fgf13a:EGFP+ cardiomyocytes were counted as PCs, while the DsRed2-NLS+/Amhc+/fgf13a:EGFP- cardiomyocytes were counted as atrium. For AVC counting, DsRed2-NLS+/Amhc+/Vmhc- cardiomyocytes were counted as atrium, DsRed2-NLS+/Amhc+/Vmhc+ cardiomyocytes were counted as AVC, and DsRed2-NLS+/Amhc-/Vmhc+ cardiomyocytes were counted as ventricle.
EdU assay for cell proliferation
Request a detailed protocolEdU labeling was performed using a Click-iT EdU Alexa Fluor Imaging Kit (Molecular Probes) according to the manufacturer’s instructions. At 48 or 72 hpf, embryos carrying the Tg(myl7:DsRed2-NLS) transgene were placed in glass vials (30 embryos per vial) and incubated with 10 mM EdU for 30 min on ice. EdU was then washed out and embryos were returned to Petri dishes and incubated at 28.5°C until 72 or 96 hpf, respectively. Embryos were then fixed in 1% formaldehyde in PBS and IHC was performed using saponin protocol listed above. After IHC, embryos were post-fixed in 2% formaldehyde in PBS for 1 hr at room temperature followed by washes with 1× PBS. The Click-iT reaction was then performed according to the manufacturer’s protocol. Embryos were imaged on a Nikon A1R inverted confocal, as described above. In each chamber, the number of EdU+/DsRed2-NLS+ cardiomyocytes were counted and normalized to the total number of DsRed2-NLS+ cardiomyocytes to determine the proliferation index.
Cardiomyocyte temporal differentiation assay
Request a detailed protocolEmbryos from a cross of adult nr2f1a+/-; Tg(myl7:NLS-KikGR) zebrafish and nr2f1a+/-; Et(fgf13a:EGFP) zebrafish were sorted for both the Et(fgf13a:EGFP) and Tg(myl7:NLS-KikGR) transgenes. As the Et(fgf13a:EGFP) is expressed in the PCs and the Tg(myl7:NLS-KikGR) is expressed in all cardiomyocytes, we were able to distinguish between these two transgenes before photoconversion of the NLS-KikGR. At 72 hpf, embryos carrying both transgenes were anesthetized in 0.16 mg/mL tricaine and exposed to UV light using a DAPI filter on Zeiss M2BioV12 Stereo microscope for 30 min to photoconvert the NLS-KikGR. At 96 hpf, photoconverted embryos were fixed in 1% formaldehyde in PBS for 1 hr at room temperature, then washed 1× in PBS followed by 2× in 0.1% saponin/1× PBS. Embryos were imaged the same day they were fixed using a Nikon A1R inverted confocal microscope with an HD resonance scanner. Nikon’s Denoise-AI was employed on acquired images.
Heart rate analysis
Request a detailed protocolTo analyze heart rate, embryos were anesthetized in 0.16 mg/mL tricaine and mounted in 1% low-melt agarose. Hearts were imaged on a Nikon Ti-2 SpectraX Widefield microscope with an Andor Xyla 4.2 megapixel, 16-bit sCMOS monochromatic camera at 200 frames per second (fps) in a controlled temperature chamber (28.5°C). Analysis performed in Nikon Elements software.
Heat-shock experiments
Request a detailed protocolEmbryos resulting from a cross of adult nr2f1a+/-; Et(fgf13a:EGFP); Tg(myl7:DsRed2-NLS); and nr2f1a+/; Tg(hsp70l:nkx2.5-EGFP) zebrafish were first sorted for the Et(fgf13a:EGFP) transgene, as it is also expressed at low levels in the epidermis and in the olfactory pits (Poon et al., 2016). At 20 or 40 hpf, 10 Et(fgf13a:EGFP)+ embryos per tube were placed into 0.5 mL PCR tubes with 200 µL embryo water and heat-shocked for 30 min by raising the temperature to 37°C using a thermocycler. Embryos were then returned to Petri dishes and placed in a 28.5°C incubator. 3 to 4 hr after the heat-shock, when the Nkx2.5-EGFP becomes visible, embryos were sorted based on the presence of Nkx2.5-EGFP, which could be distinguished over the fluorescence from the Et(fgf13a:EGFP) transgene as the Nkx2.5-EGFP signal is much stronger than that of the Et(fgf13a:EGFP). Nkx2.5-EGFP is not visible by 10 hr post heat-shock (George et al., 2015) and, as such, does not interfere with quantification of fgf13a:EGFP+ cardiomyocytes. Embryos were then returned to the 28.5°C incubator and allowed to develop until the desired timepoints when they were harvested, processed for IHC, imaged with confocal microscopy, and the cardiomyocytes counted as described above. Embryos without the Tg(hsp70l:nkx2.5-EGFP) transgene that were heat-shocked served as controls. All embryos were genotyped for the nr2f1a mutant allele as well as the Tg(hsp70l:nkx2.5-EGFP) transgene following analysis, to ensure embryos were sorted correctly after the heat-shock.
Fluorescence-activated cell sorting
Request a detailed protocolEGFP+ cardiomyocytes from WT and nr2f1a mutant with the Tg(amhc:EGFP), Tg(myl7:EGFP), or Tg(amhc:EGFP); Tg(myl7:DsRed2-NLS) transgenes were isolated using FACS at 48 or 96 hpf, as indicated. Briefly, embryos were dissociated as described previously (Holowiecki et al., 2020; Samsa et al., 2016; Stachura and Traver, 2011) with the following modifications. WT and nr2f1a mutant embryos were harvested, rinsed in Hank’s Balanced Salt Solution with Ca2+ and Mg2+ and transferred to Eppendorf tubes in 500 µL of FACS media (0.9× PBS, 2% fetal bovine serum). 60 µL of Liberase (Roche) was added and a pipette was used to triturate the embryos. Embryos were then incubated at 32.5°C for 15 min before being triturated again. This process was repeated until the embryos were completely dissociated. Following dissociation, the cell suspensions were centrifuged at 1300 rpm (100 × g) for 5 min at 4°C. The supernatant was discarded, and the pellets were resuspended in a total volume of 400 µL. Samples were then filtered through a 40 µm strainer and collected into a FACS tube. FACS was performed by the CCHMC Research Flow Cytometry Core using a 70 µm nozzle and cells were collected into Eppendorf tubes for downstream applications.
Real-time qPCR
Request a detailed protocol2229 WT and 704 nr2f1a mutant Tg(myl7:EGFP)+ were captured using flow-sorting at 96 hpf. PolyA RNA was collected from the cardiomyocytes using the Single Cell RNA Purification Kit (Norgen Biotek) and then amplified using the SeqPlex RNA Amplification Kit (Sigma-Aldrich), which generates double-stranded cDNA. RT-qPCR was performed using SYBR Green Mix in a Bio-Rad CFX-96 PCR machine. Expression levels were standardized to β-actin and analyzed using the 2−ΔΔCT Livak method. Primer sequences are provided in Supplementary file 3.
RNA-seq and ATAC-seq analysis
Request a detailed protocolFor RNA-seq, 7797 WT and 4118 nr2f1a mutant Tg(amhc:EGFP)+ were captured from flow-sorting of ~200 WT and nr2f1a mutant embryos at 48 hpf. RNA was collected using the Single Cell RNA Purification Kit (Norgen Biotek). PolyA RNA was then amplified using the MessageAmpII aRNA Amplification Kit (Thermo Fisher) and submitted to the CCHMC Sequencing core for quality control and library generation. Single-end sequencing was performed. RNA-seq data were submitted to the CCHMC Bioinformatics (BMI) core for analysis. To assess the need for trimming of adapter sequences and bad quality segments, the RNA-seq reads in FASTQ format were first subjected to quality control using FastQC v0.11.5, Trim Galore! V0.4.2, and Cutadapt v1.9.1 (Andrews, 2010; Krueger, 2012; Martin, 2011). The trimmed reads were aligned to the reference zebrafish genome version GRCz10/danRer10 with the program STAR v2.5.2 (Dobin et al., 2013). Aligned reads were stripped of duplicate reads using Picard v1.89 (Broad Institute, 2018). Gene-level expression was assessed by counting features for each gene, as defined in the NCBI’s RefSeq database (Li et al., 2009). Read counting was done with the program feature Counts v1.5.3 from the Rsubread package (Liao et al., 2019). Raw counts were normalized as transcripts per million. Differential gene expressions between conditions were assessed with the R v3.4.4 package DESeq2 v1.18.1 (Love et al., 2014).
For ATAC-seq, 8230 WT and 4338 nr2f1a mutant Tg(amhc:EGFP)+; Tg(myl7:DsRed2-NLS)+ were captured from flow-sorting of ~200 WT and nr2f1a mutant embryos at 48 hpf. ATAC-seq libraries were generated from the cells using the Nextera Kit (Illumina), as previously reported (Buenrostro et al., 2015). Libraries were submitted to the CCHMC Sequencing core for paired-end sequencing. ATAC-seq data were then submitted to the CCHMC BMI Core for analysis. Reads in FASTQ format were subjected to quality control, as above. The trimmed reads were aligned to the zebrafish genome (GRCz10/danRer10) using STAR v2.5.2 with splice awareness turned off. Aligned reads were stripped of duplicate reads using Picard v1.89. Peaks were called using MACS v2.1.2 using the broad peaks mode (Zhang et al., 2008). Peaks with q-value <0.01 were selected for further analysis. Common peaks among all samples were obtained by merging called peaks at 50% overlap using BEDtools v2.27.0 (Quinlan and Hall, 2010). The common peaks, originally in BED format, were converted to a Gene Transfer Format to enable fast counting of reads under the peaks with the program feature Counts v1.5.3 (Rsubread package) (Liao et al., 2019). Differential open chromatin relative to control samples was assessed with the R package DESeq2 v1.18.1 (Love et al., 2014). Each differentially open chromatin region was annotated with closest differentially expressed gene from the RNA-seq. Differentially open chromatin regions are selected for motif enrichment analysis using HOMER v4.10 (Heinz et al., 2010) and known zebrafish motifs from CIS-BP v1.92 (Weirauch et al., 2014) are used to detect motif occurrences.
RNA-seq and ATAC-seq data used are deposited in GEO - accession number GSE194054.
Electrophysiological analysis of hearts
Request a detailed protocolOptical mapping of action potential duration and conduction velocities were performed similar to what was previously reported (Colombo et al., 2018; Mosimann et al., 2015; Panáková et al., 2010). Briefly, hearts from zebrafish embryos were isolated at 48, 72, and 96 hpf and stained with the transmembrane potential-sensitive dye FluoVolt (Invitrogen) for 20 min. To measure action potential duration, fluorescence intensities were recorded with a high-speed CCD camera (Redshirt Imaging) with 14-bit resolution at 2000 fps. Action potentials were temporally (800 Hz low-pass cutoff) and spatially (3×3 pixel average) filtered to enhance signal-to-noise ratios. Action potential duration at 20% repolarization (APD20) was defined as the absolute time difference between the decay phase and the upstroke phase of the local action potential at 20% repolarization and was measured at the midpoint of the atria. Conduction velocity vectors across the hearts were derived using an established algorithm (Bayly et al., 1998) with custom scripts in MATLAB (Version R2018B, Mathworks) and estimated from the local action potentials (Panáková et al., 2010). All electrophysiological measurements were averaged across the relevant anatomical segments for each heart. Embryos were genotyped following analysis.
Statistical analysis
Request a detailed protocolComparisons between groups were analyzed using Student’s t-test and one- and two-way ANOVAs with multiple comparisons as appropriate. Fisher’s exact test was used to determine if two proportions were statistically distinct. Statistical analysis used for the different analyses is indicated in the figure legends. Statistical analysis was performed using GraphPad Prism. A p-value of <0.05 was considered statistically significant for all analysis. Error bars indicate s.e.m.
Materials availability
Request a detailed protocolMaterials created that were used in this study are available from the authors without restrictions.
Data availability
Sequencing data have been deposited in GEO - accession number GSE194054.
-
NCBI Gene Expression OmnibusID GSE194054. Nr2f1a maintains atrial nkx2.5 expression to repress pacemaker identity within venous atrial cardiomyocytes of zebrafish.
References
-
Rare variants in NR2F2 cause congenital heart defects in humansAmerican Journal of Human Genetics 94:574–585.https://doi.org/10.1016/j.ajhg.2014.03.007
-
Estimation of conduction velocity vector fields from epicardial mapping dataIEEE Transactions on Bio-Medical Engineering 45:563–571.https://doi.org/10.1109/10.668746
-
Mutations in the cardiac transcription factor Nkx2.5 affect diverse cardiac developmental pathwaysThe Journal of Clinical Investigation 104:1567–1573.https://doi.org/10.1172/JCI8154
-
Coding exons function as tissue-specific enhancers of nearby genesGenome Research 22:1059–1068.https://doi.org/10.1101/gr.133546.111
-
Chip-Seq identification of weakly conserved heart enhancersNature Genetics 42:806–810.https://doi.org/10.1038/ng.650
-
Calcium signalling during excitation-contraction coupling in mammalian atrial myocytesJournal of Cell Science 119:3915–3925.https://doi.org/10.1242/jcs.03223
-
Axial tubule junctions control rapid calcium signaling in atriaThe Journal of Clinical Investigation 126:3999–4015.https://doi.org/10.1172/JCI88241
-
Chamber-Specific cardiac expression of Tbx5 and heart defects in Holt-Oram syndromeDevelopmental Biology 211:100–108.https://doi.org/10.1006/dbio.1999.9298
-
Cardiomyopathy in irx4-deficient mice is preceded by abnormal ventricular gene expressionMolecular and Cellular Biology 21:1730–1736.https://doi.org/10.1128/MCB.21.5.1730-1736.2001
-
Atac-Seq: A method for assaying chromatin accessibility genome-wideCurrent Protocols in Molecular Biology 109:21.https://doi.org/10.1002/0471142727.mb2129s109
-
On the evolution of the cardiac pacemakerJournal of Cardiovascular Development and Disease 4:4.https://doi.org/10.3390/jcdd4020004
-
Development of the pacemaker tissues of the heartCirculation Research 106:240–254.https://doi.org/10.1161/CIRCRESAHA.109.205419
-
Bmp signaling exerts opposite effects on cardiac differentiationCirculation Research 110:578–587.https://doi.org/10.1161/CIRCRESAHA.111.261172
-
Star: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Cardiac homeobox gene Nkx2-5 mutations and congenital heart disease: associations with atrial septal defect and hypoplastic left heart syndromeJournal of the American College of Cardiology 41:2072–2076.https://doi.org/10.1016/s0735-1097(03)00420-0
-
Ectopic expression of Nkx2.5 suppresses the formation of the sinoatrial node in miceDevelopmental Biology 356:359–369.https://doi.org/10.1016/j.ydbio.2011.05.663
-
Tbx5 loss-of-function mutation contributes to atrial fibrillation and atypical Holt-Oram syndromeMolecular Medicine Reports 13:4349–4356.https://doi.org/10.3892/mmr.2016.5043
-
The incidence of congenital heart diseaseJournal of the American College of Cardiology 39:1890–1900.https://doi.org/10.1016/s0735-1097(02)01886-7
-
Correlation between electrical activity and the size of rabbit sino-atrial node cellsThe Journal of Physiology 496 ( Pt 3):795–808.https://doi.org/10.1113/jphysiol.1996.sp021728
-
Tbx3 controls the sinoatrial node gene program and imposes pacemaker function on the atriaGenes & Development 21:1098–1112.https://doi.org/10.1101/gad.416007
-
The form and nature of the muscular connections between the primary divisions of the vertebrate heartJournal of Anatomy and Physiology 41:172–189.
-
The tol2kit: a multisite gateway-based construction kit for Tol2 transposon transgenesis constructsDevelopmental Dynamics 236:3088–3099.https://doi.org/10.1002/dvdy.21343
-
The sequence alignment/map format and samtoolsBioinformatics 25:2078–2079.https://doi.org/10.1093/bioinformatics/btp352
-
Transcription factor ISL1 is essential for pacemaker development and functionThe Journal of Clinical Investigation 125:3256–3268.https://doi.org/10.1172/JCI68257
-
Functional redundancy between human SHOX and mouse Shox2 genes in the regulation of sinoatrial node formation and pacemaking functionThe Journal of Biological Chemistry 286:17029–17038.https://doi.org/10.1074/jbc.M111.234252
-
Epidemiology of cardiovascular malformations: prevalence and risk factorsAmerican Journal of Medical Genetics 97:319–325.https://doi.org/10.1002/1096-8628(200024)97:4<319::aid-ajmg1283>3.0.co;2-e
-
Tbx5 mutations contribute to early-onset atrial fibrillation in Chinese and CaucasiansCardiovascular Research 109:442–450.https://doi.org/10.1093/cvr/cvw003
-
Genesis and regulation of the heart automaticityPhysiological Reviews 88:919–982.https://doi.org/10.1152/physrev.00018.2007
-
Nkx2.5 mutations in patients with congenital heart diseaseJournal of the American College of Cardiology 42:1650–1655.https://doi.org/10.1016/j.jacc.2003.05.004
-
Functional role of voltage gated Ca (2+) channels in heart automaticityFrontiers in Physiology 6:19.https://doi.org/10.3389/fphys.2015.00019
-
Molecular pathway for the localized formation of the sinoatrial nodeCirculation Research 100:354–362.https://doi.org/10.1161/01.RES.0000258019.74591.b3
-
Cardiac chamber formation: Development, genes, and evolutionPhysiological Reviews 83:1223–1267.https://doi.org/10.1152/physrev.00006.2003
-
Chamber identity programs drive early functional partitioning of the heartNature Communications 6:8146.https://doi.org/10.1038/ncomms9146
-
5.78 MB terminal deletion of chromosome 15q in a Girl, evaluation of NR2F2 as candidate gene for congenital heart defectsEuropean Journal of Medical Genetics 54:354–356.https://doi.org/10.1016/j.ejmg.2010.12.004
-
Nkx2-5 suppresses the proliferation of atrial myocytes and conduction systemCirculation Research 114:1103–1113.https://doi.org/10.1161/CIRCRESAHA.114.303219
-
Differential gene expressions in atrial and ventricular myocytes: insights into the road of applying embryonic stem cell-derived cardiomyocytes for future therapiesAmerican Journal of Physiology. Cell Physiology 299:C1234–C1249.https://doi.org/10.1152/ajpcell.00402.2009
-
The orphan nuclear receptor COUP-TFII is required for angiogenesis and heart developmentGenes & Development 13:1037–1049.https://doi.org/10.1101/gad.13.8.1037
-
Development of the cardiac conduction system in zebrafishGene Expression Patterns 21:89–96.https://doi.org/10.1016/j.gep.2016.08.003
-
Proportional growth failure and oculocutaneous albinism in a girl with a 6.87 MB deletion of region 15q26.2→qterEuropean Journal of Medical Genetics 50:432–440.https://doi.org/10.1016/j.ejmg.2007.08.003
-
A novel NR2F2 loss-of-function mutation predisposes to congenital heart defectEuropean Journal of Medical Genetics 61:197–203.https://doi.org/10.1016/j.ejmg.2017.12.003
-
Cardiomyocyte FGF signaling is required for Cx43 phosphorylation and cardiac gap junction maintenanceExperimental Cell Research 319:2152–2165.https://doi.org/10.1016/j.yexcr.2013.05.022
-
Isolation and characterization of single cells from zebrafish embryosJournal of Visualized Experiments 1:53877.https://doi.org/10.3791/53877
-
Cellular dissection of zebrafish hematopoiesisMethods in Cell Biology 101:75–110.https://doi.org/10.1016/B978-0-12-387036-0.00004-9
-
Nkx genes are essential for maintenance of ventricular identityDevelopment 140:4203–4213.https://doi.org/10.1242/dev.095562
-
Novel de novo pathogenic variant in the NR2F2 gene in a boy with congenital heart defect and dysmorphic featuresAmerican Journal of Medical Genetics. Part A 176:1423–1426.https://doi.org/10.1002/ajmg.a.38700
-
Birth prevalence of congenital heart disease worldwideJournal of the American College of Cardiology 58:2241–2247.https://doi.org/10.1016/j.jacc.2011.08.025
-
Conditional knockout of fgf13 in murine hearts increases arrhythmia susceptibility and reveals novel ion channel modulatory rolesJournal of Molecular and Cellular Cardiology 104:63–74.https://doi.org/10.1016/j.yjmcc.2017.01.009
-
Modulation of pacemaker activity of sinoatrial node cells by electrical load imposed by an atrial cell modelThe American Journal of Physiology 269:H1735–H1742.https://doi.org/10.1152/ajpheart.1995.269.5.H1735
-
Localization of islet-1-positive cells in the healthy and infarcted adult murine heartCirculation Research 110:1303–1310.https://doi.org/10.1161/CIRCRESAHA.111.259630
-
BookThe Zebrafish Book. A Guide for the Laboratory Use of Zebrafish Danio rerioUniversity of Oregon Press.
-
Arrhythmias and conduction disorders associated with atrial septal defectsJournal of Thoracic Disease 10:S2940–S2944.https://doi.org/10.21037/jtd.2018.08.27
-
Isl2b regulates anterior second heart field development in zebrafishScientific Reports 7:41043.https://doi.org/10.1038/srep41043
-
Atrial identity is determined by a COUP-TFII regulatory networkDevelopmental Cell 25:417–426.https://doi.org/10.1016/j.devcel.2013.04.017
-
Mutational spectrum of the Nkx2-5 gene in patients with lone atrial fibrillationInternational Journal of Medical Sciences 11:554–563.https://doi.org/10.7150/ijms.8407
-
Model-Based analysis of ChIP-Seq (MACS)Genome Biology 9:R137.https://doi.org/10.1186/gb-2008-9-9-r137
Decision letter
-
Sigolène M MeilhacReviewing Editor; Imagine-Institut Pasteur, France
-
Didier YR StainierSenior Editor; Max Planck Institute for Heart and Lung Research, Germany
-
Sigolène M MeilhacReviewer; Imagine-Institut Pasteur, France
-
Vincent M ChristoffelsReviewer; Amsterdam University Medical Centers, Netherlands
Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.
Decision letter after peer review:
Thank you for submitting your article "Nr2f1a maintains atrial nkx2.5 expression to repress pacemaker identity within venous atrial cardiomyocytes" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Sigolène M Meilhac as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Edward Morrisey as the Senior Editor. The following individuals involved in the review of your submission have agreed to reveal their identity: Vincent M. Christoffels (Reviewer #3).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
As you will see, reviewers refer to a strong paper addressing patterning of the atrial chamber and a striking new mechanism of pacemaker cell restriction, but they also indicate that further work would be necessary for publication. In particular, three major issues have been raised.
1) The claim of transdifferentiation is currently not supported by the data. The definition of AC/VC/PC requires clarification, as well as that of the boundary between chambers. The same definition should be used throughout the paper. The tracking of a limited number of markers (fgf13a:GFP transgene for example) to ascribe cell lineage or inconsistencies between markers (cells that take on ventricular gene expression in mutants have more of an atrial electrophysiological character) raise concerns about the conclusion of transdifferentiation. Inconsistencies between figures (1D versus 1J-J, 2F versus 4F, 1B and 2K, 1J-L versus S1 and 2G-I) should be resolved. An in-depth analysis of the ATAC-seq and RNA-seq data in wild-type and mutant atrial cells and comparison with epigenetic state data from other models would provide better insights into lineage changes. In addition, cell tracking is required to conclude transdifferentiation.
2) Alternative mechanisms through which pacemaker cells are expanded should be addressed experimentally: cell death, cell proliferation, or progenitor recruitment.
3) The identification of nkx2.5 as a target of nr2f1a to repress pacemaker cells needs strengthening. Please explain how nrf2f1a is envisioned to act as a suppressor of pacemaker phenotype (atrial-specific activator of nkx2-5) as the gene is expressed in pacemaker cells and in the posterior second heart field progenitors the atria and sinus venosus (including pacemaker cells) derive from. Maintenance of nkx2.5 expression by nr2f1a at the venous pole is incongruent with the atrial to ventricular shifts observed. More details on nr2f1a expression and co-expression with nkx2.5 and pacemaker markers would be useful. The authors assume the identified enhancer regulates nkx2-5, but do not provide evidence for this. Overexpression of nkx2.5 not completely suppresses but delays the appearance of pacemaker cells. Examination of nkx2.5 rescue should be extended at a functional level (cardiac function and electrophysiology).
Detailed comments of the reviewers are below. Please take into account all comments and provide more explanations on the methods as detailed by the reviewers. If you are not able to address experimentally some reviewer comments within the revision period, please discuss the limitations or alternative hypotheses in the manuscript.
Reviewer #1 (Recommendations for the authors):
– Please provide a clear indication of the biological system under investigation in the title.
– Abbreviations complicate reading.
– How fluorescence of Nkx2.5-EGFP is distinguished over the Et(fgf13a:EGFP) transgene?
– Was cell counting done in 3D or in 2D in a projection image? Can a partial cell count in 2D explain variations between Fig2G+H and 1J? Please clarify the method section.
– Please complete the criteria of differential expression in the RNA-seq (p-value or adjusted p-value ? cut-off ?).
– "ectopic expression" (line 139) is not correct since amhc positive cells normally do express ventricular and pacemaker markers.
– Two mutant lines are described in the method section. Which line is used for each experiment? How have they been validated for the loss of nr2f1a expression?
– How are vmhc and AVC markers in the transcriptomic dataset?
– Provide single-channel images for Figure 1C-H, 2A-F to clarify marker overlap.
Reviewer #2 (Recommendations for the authors):
1. Figure 1: the overlap in Amhc and Vmhc should be more clearly shown. As there appears to be a significant overlap by 48 hpf (1D), what do these markers show at earlier time points? Finally, the image shown in 1D does not appear to be congruent with the data in 1J-L – based on the image, I would expect to see more vmhc + CMs in the mutant atrium at 48 hpf.
2. Line 168: The cell death and proliferation data are important to argue it is a fate conversion vs loss/expansion of populations that is occurring. It should be shown, or some form of fate tracing (photoconversion, etc.) should be used if statements like those in lines 174/5 are to be made.
3. Figure 4: The extent of fgf13a:GFP + cells in the atrium in nr2f1a mutants at 96hpf appears to differ markedly in Figures 2 (most of the cells are +) vs 4 (appears to be a lower proportion of the atrium as shown). How are chamber boundaries being defined in these analyses?
4. Figure 5: Somewhat confusing – how are new PCs made from 48 to 96hpf if Nkx2.5 is still present? ACs are still lost in "rescued" hearts – presumably, these now make VCs? Importantly, does the nkx2.5 transgene rescue/normalize heart rate or conduction? It is important to distinguish an effect on the fgf13a transgene (which is convincing) versus a more important effect on PC and AC fate and FUNCTION is reflected in these results.
Reviewer #3 (Recommendations for the authors):
1) Please provide more details in figure 1 regarding the RNAseq analysis. Why did the authors use a cut-off of >2 fold difference rather than significance? Could they provide an MA plot or so to show distributions of FC and p values of all transcripts detected? How many transcripts were differentially expressed? What kind of biological processes were changed?
2) In figure 1B the behavior of pacemaker markers is inconsistent. Please provide data or references showing the expression of these markers in fish pacemaker cells (e.g. I did not find these data for tbx3a etc.). Furthermore, central pacemaker markers like isl1a (Isl1 homologue) and shox2a (Shox2 homologue) are downregulated in mutants? That is unexpected, and also inconsistent with the data given in figure 2K.
3) Arterial pole of the atrium is confusing, as arterial pole usually refers to the actual arterial pole of the heart (distal/downstream portion of the ventricle).
4) Please explain how the authors envision nrf2f1a/Nrf2f2 acts as a suppressor of pacemaker phenotype (atrial-specific activator of nkx2-5) as the gene is expressed in pacemaker cells and in the posterior second heart field progenitors the atria and sinus venosus (including pacemaker cells) derive from.
5) The cell numbers seem inconsistent between the different figures. In suppl Figure 1, there are about 75 amhc+ cells in wt, and 20 or so in mutants. In figure 1, there are close to 100 amhc+ cells in wt, and about 50 in mutants? In figure 2 there are 60 amhc+ cells in wt, and 40, 20 and 10 amhc+ cells in mutants at subsequent stages. If one adds the amhc+ and fgf13a-gfp+ cells, the numbers are 40, 30, 20, respectively (figure 2c). Which numbers are correct? Is the number of amhc+ cells very low in mutants from stage 48 hpf onwards, without changing, or is there a decrease during development?
A more general issue related to the cell number issues, the data can be explained by transdifferentiation of atrial cells into pacemaker cells, as the authors seem to favor, or by reduced addition of cardiomyocytes fated to become atrial, while the normal atrial-pacemaker demarcation programs are still active, or by reduced atrial differentiation. The latter 2 options are consistent with the data of suppl Figure 1. The authors should resolve this, and provide additional support for the trans-differentiation hypothesis. Thus, does a fraction of cells expressing pacemaker markers derive from atrial marker expressing cells in mutants, or is there an increased population of "de novo" pacemaker cells (not previously atrial marker positive) in the venous pole and a decreased population of atrial cells because of defective recruitment from progenitors or defective atrial differentiation.
6) Line 201: "Thus, we would predict…" Why would the heart rate be predicted to decrease? More pacemaker cells in mutants may also result in increased automaticity and increased spontaneous depolarizations.
7) The relation between hyperpolarization, increased repolarization times and increased numbers of pacemaker cells in mutants is not clear to me. Pacemaker cells have increased, not decreased resting membrane potentials, decreased Vmax (due to lower sodium channel activity), and unstable resting membrane potentials. The authors could analyse the electrophysiological properties of the pacemaker cells in wt and mutants (fgf13a-gfp label is perfect for that purpose), and also investigate the transcriptomes for differential expression of particular ion handling proteins known to be involved in the action potential morphology, membrane potential and intercellular conduction.
8) Line 217 "…slower conduction…" and line 219 "… more conduction system-like…" If anything, the conduction system shows very high conduction velocities. Perhaps refer to pacemaker properties. The remark regarding Vmax here is not clear. Vmax is low in pacemaker cells compared to atrial cells.
9) Line 242, transdifferentiation. See comment 4.
10) The authors assume the identified enhancer regulates nkx2-5, but do not provide evidence for this. The enhancer could be removed from the genome, for example. Did the authors analyse the topologically associated domain of nkx2-5 to see whether the enhancer is part of it? Is any interaction data available?
Please provide more detailed data regarding the ATAC-seq experiment. How many accessible sites, how many differential accessible sites, examples of known atrial accessible sites, UCSC browser views of a few examples to get an impression regarding data quality and signal-to-noise ratios. What is the correlation between RNAseq identified target genes and differential accessibility?
https://doi.org/10.7554/eLife.77408.sa1Author response
Essential revisions:
As you will see, reviewers refer to a strong paper addressing patterning of the atrial chamber and a striking new mechanism of pacemaker cell restriction, but they also indicate that further work would be necessary for publication. In particular, three major issues have been raised.
(1.A) The claim of transdifferentiation is currently not supported by the data.
This first comment is essentially the same as Essential Revision comment #2. Please see the response to Essential Revision comment #2 below for how we specifically addressed issues of transdifferentiation vs progenitor recruitment.
(1.B) The definition of AC/VC/PC requires clarification, as well as that of the boundary between chambers. The same definition should be used throughout the paper.
Thank you for this suggestion to help clarify our references to the specific cardiomyocyte populations. In the revised manuscript, we have now changed the references to cardiomyocyte populations examined, such as Amhc+/Vmhc-, so that they are defined primarily by their expression of the differentiation and transgenic markers used for the specific analysis. We then state our interpretations with respect to what that means regarding the atrial, atrioventricular canal, ventricular, and pacemaker cardiomyocyte populations. The boundaries between chambers and cardiomyocyte populations were determined by the expression of the indicated markers. The references to the cardiomyocyte populations with the designated markers have been changed throughout the text and figure legends, the main figures, and the figure supplements.
(1.C) The tracking of a limited number of markers (fgf13a:GFP transgene for example) to ascribe cell lineage or inconsistencies between markers (cells that take on ventricular gene expression in mutants have more of an atrial electrophysiological character) raise concerns about the conclusion of transdifferentiation.
a. We acknowledge that unfortunately there are a limited number of specific markers to use for some of the cardiomyocyte populations, in particular the pacemaker cardiomyocytes within the zebrafish hearts. This is something we have struggled with in our experiments. We did not use some established markers in that analysis, such as Bmp4, because they are expressed in the atrioventricular canal and pacemaker cells at the venous pole. This created some ambiguity to effects seen at later stages in the mutants due to the expression in both regions of the heart and the smaller atria. With other markers that should be expressed more specifically in the pacemaker cardiomyocytes, such as shox and shox2, there were technical issues in that the probes for both standard and HCR in situ hybridization, while marking the population, were not consistent or robust enough for us to perform quantitative analysis. Therefore, we primarily relied on immunohistochemistry for Isl1 and the fgf13a:EGFP reporter, because we found they were the most robust markers for the pacemaker cardiomyocytes currently available for analysis in the zebrafish hearts.
b. Regarding the electrophysiological characteristics (lines 248-280 of the revised manuscript), for the analysis of the action potentials we focused on regions that are closer to the middle of the morphological atrium where cardiomyocytes would express Amhc+ and fgf13a:EGFP+ and not express ventricular markers within the enlarged atrioventricular canal of nr2f1a mutant hearts based on our immunohistochemical analysis. Our data does show that this atrial region has electrophysiological properties more consistent with the central conduction system and not the WT atrial cardiomyocytes. The atrioventricular canal region also has altered conduction properties, as there is less slowing of conduction in the nr2f1a mutant atrioventricular canals. The lack of slowing of conduction at the atrioventricular canal of nr2f1a mutants is consistent with these cardiomyocytes taking on ventricular-like conduction properties and the gene expression changes we observe. Despite the improper conduction properties, as we showed previously larger endocardial cushions are still formed (Duong et al., 2018), supporting that the Amhc+/Vmhc+ cardiomyocytes still must have atrioventricular cardiomyocyte properties that promote valve formation. Please also see responses to Reviewer #3’s comments 6, 7, and 8.
(1.D) Inconsistencies between figures (1D versus 1J-J, 2F versus 4F, 1B and 2K, 1J-L versus S1 and 2G-I) should be resolved.
Thank you for point these out. We have tried to address each of the inconsistencies for the listed figures.
a. With respect to Figures 1D vs 1J-J, it appears the origin of this comment is Reviewer #2’s comment 1, “Finally, the image shown in 1D does not appear to be congruent with the data in 1J-L – based on the image, I would expect to see more vmhc + CMs in the mutant atrium at 48 hpf.” In considering this comment, it is not clear to us what Reviewer #2 found inconsistent about this data. Our data show there is an increase in total Vmhc+ cardiomyocytes in the nr2f1a mutant hearts compared to WT sibling hearts. However, the increase in total Vmhc+ cardiomyocytes comes from the combined Amhc-/Vmhc+ and Amhc+/Vmhc+ cardiomyocytes populations, which are marking ventricular and atrioventricular canal cardiomyocytes respectively. There is not an increase in the Vmhc+-only (Amhc-/Vmhc+) ventricular cardiomyocytes in the nr2f1a mutants relative to the WT control embryos. This is illustrated in Figure 1 —figure supplement 2 of the revised manuscript (Supplemental Figure 1 of the original manuscript), which shows there is not a statistically significant difference in the number of Vmhc+-only (Amhc-/Vmhc+) cardiomyocytes from these conditions from 48-96 hpf. To address this comment and better convey the expansion of Vmhc expression within the nr2f1a mutant hearts, we have provided the individual image channels in a revised Figure 1 —figure supplement 1 and also show the proportions of Vmhc+ cardiomyocyte populations in the hearts in Figure 1 —figure supplement 2D.
b. With respect to Figures 2F vs 4F, it appears this comment originated from Reviewer #2’s comment 3, “The extent of fgf13a:GFP + cells in the atrium in nr2f1a mutants at 96hpf appears to differ markedly in Figures 2 (most of the cells are +) vs 4 (appears to be a lower proportion of the atrium as shown). How are chamber boundaries being defined in these analyses?” We respectfully disagree that there is a significant difference in what is being shown between the effects in these images. In both Figure 2F and Figure 4F, the fgf13a:EGFP expression extends approximately half way across smaller atrium of the nr2f1a mutant hearts and overlaps significantly with Amhc+ cardiomyocytes. There can be variation in both the morphology of the hearts and the expansion of the fgf13a:EGFP expression in the nr2f1a mutants, with the representative heart in Figure 2F being more distended than Figure 4F. However, we are not basing the analysis and interpretations on cardiac or chamber morphology. We are basing the interpretation on the overlap of Amhc and fgf13a:EGFP expression. Importantly, our quantification of cardiomyocytes shows in Figure 2H that despite their morphology the nr2f1a mutant hearts have consistently more Amhc+/fgf13a:EGFP+ cardiomyocytes than sibling control embryos by 96 hpf. To address this comment and clarify our presentation of this data, we have revised Figure 4 to include the individual panels so that they are also showing Amhc expression within the hearts, which allows for better visualization of Amhc and fgf13a:EGFP expression within the atrium compared to the merged image in Figure 4F.
c. With respect to Figures 1B vs 2K, the RNA-seq in Figure 1B was performed on sorted amhc:EGFP+ cardiomyocytes at 48hpf. While Figure 2K was RT-qPCR for the respective genes from myl7:EGFP+ (pan-cardiac) cardiomyocytes flow sorted at 96hpf. We performed the experiments this way as they are complementary analyses and provide insights into changes observed at the earlier stages of analysis that generated the hypotheses and those found to change at later stages of heart development following our temporal analysis. The sorting using a pan-cardiac marker was chosen because the smaller amhc:EGFP+ population observed at 96 hpf was technically difficult to sort and so we could better analyze changes in ventricular markers within the whole hearts. Figure 2K shows that by 96 hpf, the hearts overall have increased expression of selected ventricular and pacemakers markers, consistent with what we observe with the IHC shown in Figures 1 and 2.
d. With respect to Figures 1J-L vs S1, we do not understand what exactly could be inconsistent with these graphs. As we indicated in the Figure S1 legend of the original manuscript, the data shown in Figure S1 (Figure 1 —figure supplement 2 of the revised manuscript) is the same data shown in Figures 1J-L. These are not separate experiments. We separated out the individual cardiomyocytes populations in Figure 1 —figure supplement 2 to make it easier for the reader to observe the effects on the individual cardiomyocyte populations, as the stacked graphs may be difficult to easily observe some of the differences.
e. With respect to Figure 1J-L vs 2G-I, thank you for pointing this out. We understand why there many have been some confusion. There can be some slight experiment-to-experiment variability in the number of cardiomyocytes quantified within the hearts, which can be affected by a variety of factors, including the embryos, the IHC, and the microscope used. Importantly, despite slight differences in cardiomyocyte numbers for these different experiments, the trends on the cardiomyocyte populations quantified in the graphs between the original Figure 1J-L and Figure 2G-I were the same and did not affect the interpretation that there were decreased Amhc+ cardiomyocytes and increased Vmhc+ (ventricular) and fgf13a:EGFP+ (pacemaker) cardiomyocytes in the nr2f1a mutant hearts. However, to alleviate any confusion, we have replaced the graphs in Figure 2G-I with an experimental replicate showing the same trends, but the cardiomyocytes counts were more similar to those of Figure 1J-L. Again, the trends are the same from all the biological replicates and do not affect the interpretation regarding the changes in cardiomyocyte populations.
(1.E) An in-depth analysis of the ATAC-seq and RNA-seq data in wild-type and mutant atrial cells and comparison with epigenetic state data from other models would provide better insights into lineage changes.
We appreciate this comment and suggestion. Certainly, gaining a greater understanding of the level of conservation for the roles of Nr2fs within cardiomyocytes shown in this study is a goal that we have. However, it is not clear to us from the comment how these proposed comparisons should be done. For these types of comparisons, ATAC-seq and RNA-seq would need to be available from pretty similar cell types and collected in a similar manner, as chromatin state and expression can be very dynamic depending on factors such as developmental stage and cell type. Therefore, it is not clear to us that there is an available dataset we could use for this comparison to look at conserved changes in chromatin or gene expression, and even if there was, that it would immediately provide additional insight into the lineage changes. Another problem is that many cardiac enhancers are not highly conserved in the heart compared to other tissues (Blow et al., 2010). While the putative nkx2.5 enhancer may not be highly conserved at the sequence level, we cannot rule out that is functionally conserved. However, this additional analysis would take significantly more time to perform. We have modified the text to specifically mention the lack of conservation of the putative nkx2.5 enhancer element within the heart and that this may not be unusual, in the revised manuscript (lines 335-337).
(1.F) In addition, cell tracking is required to conclude transdifferentiation.
This comment addresses the same issue as Essential Revision comment 2. We have included additional experimental data in the revised manuscript to support that cardiomyocytes are transdifferentiating in the nr2f1a mutant hearts. Please see our response to Essential Revisions comment 2 below.
2) Alternative mechanisms through which pacemaker cells are expanded should be addressed experimentally: cell death, cell proliferation, or progenitor recruitment.
A major concern from the reviewers was that they felt we had not ruled out the possibility that the additional pacemaker cardiomyocytes in the nr2f1a mutant hearts could be added from another source, such as through progenitor recruitment and newly-differentiated pacemaker cardiomyocytes at the venous pole, and were not from the transdifferentiation of atrial (Amhc+) cardiomyocytes. We did mention in the original manuscript that we did not see effects on cardiomyocyte proliferation or death within the hearts of nr2f1a mutants. However, we did not show this data as it was negative. Furthermore, our data in Figure 1L show that the nr2f1a mutant hearts effectively do not grow from 48 to 96 hpf. Therefore, it is unlikely the expansion of pacemaker cardiomyocyte gene expression within the nr2f1a mutant hearts could be explained by being from newly-added cardiomyocytes. To explain the pacemaker expansion being from newly-differentiating pacemaker cardiomyocytes with a static number of total cardiomyocytes, one would have to posit there is a treadmilling effect with atrial cardiomyocytes dying at a significant rate and being replaced by newly-differentiating pacemaker cardiomyocytes. To specifically address the issue of progenitor recruitment/newly-differentiating cardiomyocytes vs transdifferentiation, in the revised manuscript we have included the data examining cell death with active Caspase 3 and cell proliferation using both pHH3 and EdU. We also performed analysis of the temporal addition of newly-differentiating cardiomyocytes within the nr2f1a mutant hearts with the previous reported Tg(myl7:NLS-KikGR) transgene. Altogether, these data support that the Amhc+ cardiomyocytes in the venous pole of the atrium are acquiring pacemaker gene expression. The text of the revised manuscript now refers to these experiments and the new supplemental figures (lines 189-198 and 229-247, and Figure 1 —figure supplement 3, Figure 1 —figure supplement 4, and Figure 2 —figure supplement 3).
(3.A) The identification of nkx2.5 as a target of nr2f1a to repress pacemaker cells needs strengthening. Please explain how nrf2f1a is envisioned to act as a suppressor of pacemaker phenotype (atrial-specific activator of nkx2-5) as the gene is expressed in pacemaker cells and in the posterior second heart field progenitors the atria and sinus venosus (including pacemaker cells) derive from.
We agree that we do not yet understand how Nr2f1a mechanistically promotes/maintains Nkx2.5 expression in venous atrial cardiomyocytes and represses pacemaker identity in atrial cardiomyocytes. Since Nr2f1a is expressed throughout atrial cardiomyocytes (Figure 6 —figure supplement 3), including pacemaker cardiomyocytes, we hypothesize that Nr2f1a must interact with other factors and signaling to provide these regional effects. The identification of these factors and signaling extends beyond this manuscript. We specifically discuss and state this hypothesis in the Discussion of the revised manuscript (lines 436-443).
(3.B) Maintenance of nkx2.5 expression by nr2f1a at the venous pole is incongruent with the atrial to ventricular shifts observed. More details on nr2f1a expression and co-expression with nkx2.5 and pacemaker markers would be useful.
a. Unfortunately, it is not clear to us how maintenance of nkx2.5 expression is incongruent with the atrial to ventricular shifts. In nkx2.5 zebrafish mutants, the hearts fail to maintain ventricular identity and transition to atrial and pacemaker cardiomyocyte identity (Colombo et al., 2018). As Nkx2.5 expression is not lost in the atrioventricular canal of the nr2f1a mutants, but the more venous atrial cardiomyocytes, we would not predict that the ability to promote or maintain ventricular identity would be affected in the nr2f1a mutant hearts.
b. To address the comment about co-expression of Nr2f1a and Nkx2.5, we have added immunohistochemistry of Nr2f1a and Nkx2.5 with pacemaker markers. However, they cannot be examined together as the antibodies were both generated in rabbits. The data of the co-expression are mentioned in the Results and Discussion and supplemental figures of the revised manuscript (lines 287-291 and 436-443 Figure 4 —figure supplement 4 and Figure 6 —figure supplement 3).
(3.C) The authors assume the identified enhancer regulates nkx2-5, but do not provide evidence for this.
In the revised manuscript, we have examined the function of the putative nkx2.5 enhancer using CRISPR/Cas9 to delete it. We find that deletion of the putative nkx2.5 enhancer in embryos leads to a loss of Nkx2.5+ cardiomyocytes and an expansion of the fgf13a:EGFP reporter in atrial cardiomyocytes. This new experiment has been described in the Results (lines 348-364) and included in a revised Figure 6 and new Figure 6 —figure supplement 2 of the revised manuscript.
(3.D) Overexpression of nkx2.5 not completely suppresses but delays the appearance of pacemaker cells. Examination of nkx2.5 rescue should be extended at a functional level (cardiac function and electrophysiology).
We assessed heart rate following induction of Nkx2.5 at the 20 somite stage in nr2f1a mutants at 96 hpf. Despite the strong effects on the fgf13:EGFP reporter in nr2f1a mutants, this does not restore heart rate at that stage. One reason for the failure to rescue heart rate could be that at this stage the pacemaker cardiomyocytes have already started to expand. Furthermore, this result is consistent with the hypothesis that additional factors besides Nkx2.5 are needed to fully rescue atrial cardiomyocyte and restrict pacemaker cardiomyocyte identities within the venous pole of the atrium. This data is now mentioned in the Results (lines 314-319) and included in a new Figure 5 —figure supplement 1 of the revised manuscript.
(3.E) Detailed comments of the reviewers are below. Please take into account all comments and provide more explanations on the methods as detailed by the reviewers. If you are not able to address experimentally some reviewer comments within the revision period, please discuss the limitations or alternative hypotheses in the manuscript.
We have responded to all the comments and comments not addressed in the Essential Revisions below. We have expanded the Methods as requested by the reviewers, and we have included the Methods used for new experiments added to the revised manuscript.
Reviewer #1 (Recommendations for the authors):
– Please provide a clear indication of the biological system under investigation in the title.
We have revised the title to indicate that the studies use zebrafish.
– Abbreviations complicate reading.
We have tried to limit the use of abbreviations for the cardiomyocyte populations in the revised manuscript. This was in large part achieved through changing how the populations of cardiomyocytes are referenced. When we did use abbreviations, we made sure they are clearly defined.
– How fluorescence of Nkx2.5-EGFP is distinguished over the Et(fgf13a:EGFP) transgene?
As we indicated in the Methods of the original manuscript, prior to the heat-shock, we first sorted embryos for those that were transgenic for the fgf13a:EGFP expression based on its expression in the epidermis. Following the heat-shock, transgenic and non-transgenic hsp70l:nkx2.5-EGFP embryos were sorted based on induced Nkx2.5-EGFP expression, which shortly after the heat-shock is significantly brighter than the transgenic fgf13a:EGFP expression. In the revised manuscript, we have tried to clarify how we distinguished the expression of these transgenes with a statement in the Results (lines 301-303) and in the Methods (lines 584-588). In the initial manuscript, we unintentionally omitted that the genotypes of all embryos were subsequently confirmed with PCR genotyping following imaging for the presence or absence of the hsp70l:nkx2.5-EGFP transgene. This is now also included in the revised Methods (lines 591-593).
– Was cell counting done in 3D or in 2D in a projection image? Can a partial cell count in 2D explain variations between Fig2G+H and 1J? Please clarify the method section.
The cardiomyocytes were counted using MaxIP projections from confocal images of the whole hearts. We do not think partial counts could explain variation between these figures. It is more likely some experiment-to-experiment variation that can occur. Please see our response to Essential Revisions comment 1D for additional details of how we addressed this issue in the revised manuscript.
– Please complete the criteria of differential expression in the RNA-seq (p-value or adjusted p-value ? cut-off ?).
This data should have been included in Supplementary file 1. There is now a sheet with all the RNA-seq data included in a revised Supplementary file 1.
– "ectopic expression" (line 139) is not correct since amhc positive cells normally do express ventricular and pacemaker markers.
Thank you. We have modified this sentence (line 150 of the revised manuscript).
– Two mutant lines are described in the method section. Which line is used for each experiment? How have they been validated for the loss of nr2f1a expression?
The nr2f1ael512 and nr2f1aci1009 mutant alleles lines were reported previously (Duong et al., 2018). The nr2f1aci1017 allele was verified phenotypically, through complementation tests with the previous lines, and Western blots for Nr2f1a. As indicated in the text, additional details of the allele are in another manuscript that will be published shortly. As the nr2f1a mutant cardiac phenotype was published previously, it was not necessary to belabor the point of the additional alleles. The nr2f1aci1017 allele was used for the RNA and ATAC-seq experiments. The nr2f1ael512 and the nr2f1aci1009 alleles were used for the temporal differentiation assay with the myl7:KikGR. All the rest of the experiments in the manuscript were performed with the nr2f1aci1009 allele. There was no experimental rationale for using these lines for different experiments, other than they had the necessary transgenes to perform desired experiments crossed into them.
– How are vmhc and AVC markers in the transcriptomic dataset?
Vmhc (myh7) for some reason was not found in the RNA-seq data sets from sorted amhc:EGFP+ cardiomyocytes. However, other markers of ventricular identity were, as we indicated. We know from previous analysis (Duong et al., 2018) and RT-qPCR shown in Figure 2K that vmhc expression is increased in the whole hearts. Despite the trends we observe, it could reflect that a limited number of these cardiomyocytes were captured due to lower amhc expression in these cardiomyocytes from nr2f1a mutant hearts and the expression didn’t meet the defined cut-offs. Furthermore, it is difficult to distinguish some of the atrioventricular canal markers from pacemaker markers at this stage. Again, as we reported previously, when examining the heart atrioventricular canal markers are also expanded (Duong et al., 2018).
– Provide single-channel images for Figure 1C-H, 2A-F to clarify marker overlap.
Single-channel images have been provided for Figures 1 and 2 in Figure 1 —figure supplement 1 and Figure 2 —figure supplement 1 of the revised manuscript.
Reviewer #2 (Recommendations for the authors):
1. Figure 1: the overlap in Amhc and Vmhc should be more clearly shown. As there appears to be a significant overlap by 48 hpf (1D), what do these markers show at earlier time points? Finally, the image shown in 1D does not appear to be congruent with the data in 1J-L – based on the image, I would expect to see more vmhc + CMs in the mutant atrium at 48 hpf.
a. We have provided images of hearts at an earlier time points (30 hpf). These show that there is actually significant overlap of Vmhc and Amhc at these earlier stages in the heart tube, which to our knowledge has not really been reported previously at the shown stage. Importantly, we did not find a difference in the overlap of Vmhc and Amhc between WT and nr2f1a mutant hearts at this stage. This suggests that Nr2f1a is required to refine the number of Vmhc and Amhc cardiomyocytes in the atrioventricular canal, consistent with our previous report (Duong et al. 2018). This new data and related interpretations are now described in the Results and Discussion of the revised manuscript (lines 160-162 and 385-390) and in new panels included in Figure 1- supplement 1.
b. We have tried to address the issue of inconsistency between Figure 1D and 1J-L in the response to Essential Revisions comment 1D.
2. Line 168: The cell death and proliferation data are important to argue it is a fate conversion vs loss/expansion of populations that is occurring. It should be shown, or some form of fate tracing (photoconversion, etc.) should be used if statements like those in lines 174/5 are to be made.
Thanks you. We have provided the cardiomyocyte death and proliferation analysis, and additional lineage tracing/temporal differentiation assay to examine newly-differentiating cells in the nr2f1a mutant hearts to address the issue of fate conversion/transdifferentiation. Please see our response to Essential Revision comment 2 for additional details of how we addressed this issue.
3. Figure 4: The extent of fgf13a:GFP + cells in the atrium in nr2f1a mutants at 96hpf appears to differ markedly in Figures 2 (most of the cells are +) vs 4 (appears to be a lower proportion of the atrium as shown). How are chamber boundaries being defined in these analyses?
Please see our response to Essential Revision comment 1D for how we addressed this issue.
4. Figure 5: Somewhat confusing – how are new PCs made from 48 to 96hpf if Nkx2.5 is still present? ACs are still lost in "rescued" hearts – presumably, these now make VCs? Importantly, does the nkx2.5 transgene rescue/normalize heart rate or conduction? It is important to distinguish an effect on the fgf13a transgene (which is convincing) versus a more important effect on PC and AC fate and FUNCTION is reflected in these results.
a. It is difficult to address this issue because it is not clear to us from this comment what exactly Reviewer #2 finds confusing. Normally, Nkx2.5 is expressed in virtually all cardiomyocytes within the zebrafish heart except for pacemaker cardiomyocytes at the venous pole of the heart. This is now shown in Figure 4 —figure supplement 1 of the revised manuscript. In the nr2f1a mutants, Nkx2.5 is not completely lost throughout the heart. It just recedes from the venous pole of the heart, which corresponds with the new marker expression for the pacemaker cardiomyocytes (Figure 4). As Nkx2.5 is still expressed throughout the heart in the mutants, we would not expect to observe difference in the effect on ventricular cardiomyocytes.
b. In the “rescued” hearts where Nkx2.5-EGFP is induced, it is not clear that Nkx2.5-EGFP is still visible at the later stages following the heat-shock. The previous characterization of the hsp70l:nkx2.5-EGFP line showed that Nkx2.5-EGFP was no longer detectable by 10 hrs after the heat shock (George et al., 2015). While inducing Nkx2.5-EGFP at the 20 somite stage is sufficient to rescue nkx2.5 mutant embryos to adulthood (George et al., 2015), it does not mean it is perpetually expressed following the heat shock at the stages examined in the analysis. Moreover, we do not see the Nkx2.5-EGFP at 48 through 96 hpf following the heat-shock at 20 somites. Similar to the nkx2.5 mutant rescue, we see maximal repression of the pacemaker reporter (fgf13a:EGFP) when induced at the 20 somite stage. These results suggest that even though we observe the regression of Nkx2.5 at later stages in the mutants, Nkx2.5 and Nr2f1a complexes must be required earlier to establish these boundaries. We include this in the Discussion of the revised manuscript (lines 421-425 and 429).
c. We did assess rescue of heart rate following Nkx2.5-EGFP induction as requested. Please see our response to Essential Revisions comment 3.D.
Reviewer #3 (Recommendations for the authors):
1) Please provide more details in figure 1 regarding the RNAseq analysis. Why did the authors use a cut-off of >2 fold difference rather than significance? Could they provide an MA plot or so to show distributions of FC and p values of all transcripts detected? How many transcripts were differentially expressed? What kind of biological processes were changed?
We did not present the RNA-seq data this way as we did not think it provided additional information that was necessary. It served as the basis for additional in vivo investigation of a hypothesis in the embryos. All of the RNA-seq data are now included in a sheet in Supplement File 1. They were unintentionally omitted from the excel sheet originally submitted. The expression can be represented in different ways. We felt the differences were best illustrated using fold difference changes that were found in this data set.
2) In figure 1B the behavior of pacemaker markers is inconsistent. Please provide data or references showing the expression of these markers in fish pacemaker cells (e.g. I did not find these data for tbx3a etc.). Furthermore, central pacemaker markers like isl1a (Isl1 homologue) and shox2a (Shox2 homologue) are downregulated in mutants? That is unexpected, and also inconsistent with the data given in figure 2K.
a. Please see our response to Essential Revision comment 1D.
b. References to the expression of the pacemaker marker genes in the zebrafish have been added to the revised manuscript (lines 147,148, 207-209).
3) Arterial pole of the atrium is confusing, as arterial pole usually refers to the actual arterial pole of the heart (distal/downstream portion of the ventricle).
We welcome suggestions for what would be the best way to refer to this anatomical position in the context of our study as we would like to be as clear as possible. Given the changes in heart morphology and markers within the different regions of the nr2f1a mutant hearts, it was not clear how to best refer to this portion of the atrium, which is not often specifically referred to as most studies have not had reason to refer to these different regions within the zebrafish atrium. Therefore, it seemed clearest to us to use arterial and venous as relative anatomical position markers with respect to the direction of blood flow in the heart, as is done with respect to the whole hearts in fish and mammals, and contrast these positions with our references to the pacemaker cardiomyocytes at the venous pole of the heart/atrium. In the Abstract and Introduction of the revised manuscript, we have also defined arterial and venous as the outflow and inflow portion of the atrium to help clarify this point (lines 27, 29, 110, and 114).
4) Please explain how the authors envision nrf2f1a/Nrf2f2 acts as a suppressor of pacemaker phenotype (atrial-specific activator of nkx2-5) as the gene is expressed in pacemaker cells and in the posterior second heart field progenitors the atria and sinus venosus (including pacemaker cells) derive from.
Please see our responses to Essential Revisions comments 3.A and 3.B, Reviewer #1’s major comment 3.
5) The cell numbers seem inconsistent between the different figures. In suppl Figure 1, there are about 75 amhc+ cells in wt, and 20 or so in mutants. In figure 1, there are close to 100 amhc+ cells in wt, and about 50 in mutants? In figure 2 there are 60 amhc+ cells in wt, and 40, 20 and 10 amhc+ cells in mutants at subsequent stages. If one adds the amhc+ and fgf13a-gfp+ cells, the numbers are 40, 30, 20, respectively (figure 2c). Which numbers are correct? Is the number of amhc+ cells very low in mutants from stage 48 hpf onwards, without changing, or is there a decrease during development?
A more general issue related to the cell number issues, the data can be explained by transdifferentiation of atrial cells into pacemaker cells, as the authors seem to favor, or by reduced addition of cardiomyocytes fated to become atrial, while the normal atrial-pacemaker demarcation programs are still active, or by reduced atrial differentiation. The latter 2 options are consistent with the data of suppl Figure 1. The authors should resolve this, and provide additional support for the trans-differentiation hypothesis. Thus, does a fraction of cells expressing pacemaker markers derive from atrial marker expressing cells in mutants, or is there an increased population of "de novo" pacemaker cells (not previously atrial marker positive) in the venous pole and a decreased population of atrial cells because of defective recruitment from progenitors or defective atrial differentiation.
a. Please see our response to Essential Revisions comment 1.D for how we addressed concerns of the inconsistencies.
b. Please see our response to Essential Revisions comment 2 for how we experimentally addressed concerns of transdifferentiation vs. new addition of cardiomyocytes.
6) Line 201: "Thus, we would predict…" Why would the heart rate be predicted to decrease? More pacemaker cells in mutants may also result in increased automaticity and increased spontaneous depolarizations.
Thank you for the comment. The reviewer is correct. One might readily imagine ways in which increased numbers of pacemaker cells could lead to any number of heart rate outcomes depending on their cell autonomous features, their capacitance, and other biophysical properties as well as particular patterns of intercellular coupling. In our description of the anticipated heart rate findings, we neglected to fully outline our rationale. Given the typical intrinsic automaticity of a coupled cellular mass is a function of the total electrical capacitance of that mass, we simply hypothesized that the intrinsic heart rate might be slower as the mass of coupled pacemaker cells increased. Clearly, this assumes that cellular characteristics and coupling are unchanged, but as these properties were also measured directly in the work described we did not elaborate on these elements at this stage. We have revised the text to be more explicit on this point (lines 252-255).
7) The relation between hyperpolarization, increased repolarization times and increased numbers of pacemaker cells in mutants is not clear to me. Pacemaker cells have increased, not decreased resting membrane potentials, decreased Vmax (due to lower sodium channel activity), and unstable resting membrane potentials. The authors could analyse the electrophysiological properties of the pacemaker cells in wt and mutants (fgf13a-gfp label is perfect for that purpose), and also investigate the transcriptomes for differential expression of particular ion handling proteins known to be involved in the action potential morphology, membrane potential and intercellular conduction.
We apologize if the relationship between the various cellular parameters, myocardial physiology, and cell identity markers was not clear. Our intent was to convey that the atria in the nr2f1a mutants exhibit intermediate electrophysiologic characteristics consistent with cardiomyocytes of the central conduction system, although we cannot readily classify them as pacemaker, atrial, or ventricular cardiomyocytes. While we cannot perform the analysis in the fgf13a:EGFP transgenic fish, we did analyze the electrophysiology in different regions of the atria that would be consistent with them having fgf13a:EGFP+ and fgf13a:EGFP- cardiomyocytes in our studies and did not identify discrete physiological boundaries. This intermediate state which we identified in the nr2f1a-/- atria exhibits atrial Vmax and ventricular repolarization with hyperpolarization of the resting membrane potential and no significant changes in phase 4 depolarization. These features and the slowed conduction observed suggest cellular electrophysiology and cell-cell interactions most consistent with the central conduction system. This is of importance for several reasons not least the known gradients of function between atrial cardiomyocyte and central conduction system observed from mid-atrium distally in humans. We have clarified the relevant text in the Introduction, Results, and Discussion (lines 118, 252-255, 266, 267, 272, 273, 408).
8) Line 217 "…slower conduction…" and line 219 "… more conduction system-like…" If anything, the conduction system shows very high conduction velocities. Perhaps refer to pacemaker properties. The remark regarding Vmax here is not clear. Vmax is low in pacemaker cells compared to atrial cells.
We apologize that we did not make it clear that we were referring not to the ventricular (His-Purkinje) conduction system but rather to the central conduction system whose biology we believe these experiments are helping to elucidate. We have clarified the relevant text of the revised manuscript (lines 118, 252-255, 266, 267, 272, 273, 408).
9) Line 242, transdifferentiation. See comment 4.
Please see our response to Essential Revisions comment 2 for how we experimentally addressed concerns of transdifferentiation vs. new addition of cardiomyocytes.
10) The authors assume the identified enhancer regulates nkx2-5, but do not provide evidence for this. The enhancer could be removed from the genome, for example. Did the authors analyse the topologically associated domain of nkx2-5 to see whether the enhancer is part of it? Is any interaction data available?
Please provide more detailed data regarding the ATAC-seq experiment. How many accessible sites, how many differential accessible sites, examples of known atrial accessible sites, UCSC browser views of a few examples to get an impression regarding data quality and signal-to-noise ratios. What is the correlation between RNAseq identified target genes and differential accessibility?
a. Thank you for the suggestion. We have provided experimental analysis of the putative nkx2.5 enhancer. Please see response to Essential Revisions comment 3.C.
b. While it would be ideal and we would have loved to analyze the topologically association of chromatin in these zebrafish cardiomyocytes, presently it is not possible to perform the techniques, such as HI-C, to analyze chromatin conformation in these cells. Therefore, we cannot presently directly address this comment experimentally.
c. The ATAC-seq data and their relationship to changes in gene expression from the RNA-seq are summarized in a sheet in Supplementary file 1. Examples of UCSC browser views for select genes have now been added to a new Figure 6 —figure supplement 1 of the revised manuscript as requested to show the quality of the data. Additionally, Bigwig files were submitted to GEO.
https://doi.org/10.7554/eLife.77408.sa2Article and author information
Author details
Funding
National Heart, Lung, and Blood Institute (R01 HL137766)
- Joshua S Waxman
National Heart, Lung, and Blood Institute (R01 HL141186)
- Joshua S Waxman
National Heart, Lung, and Blood Institute (R01 HL154522)
- Joshua S Waxman
National Heart, Lung, and Blood Institute (F31 HL152600)
- Kendall E Martin
NIH Office of the Director (R24 OD017870)
- Calum A MacRae
Leducq Foundation
- Calum A MacRae
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank members of the Waxman lab for fruitful discussions and the CCHMC veterinary staff for care of the zebrafish. This work was supported by funding from the National Institutes of Health: [R01 HL137766, R01 HL141186, and R01 HL154522 to JSW], [F31 HL152600 to KEM], and [R24 OD017870 and the Leducq Foundation to CAM].
Ethics
All zebrafish husbandry and experiments were performed in accordance with protocols approved by the Cincinnati Children's Hospital Medical Center Institutional Animal Care and Use Committee (IACUC). Protocol # 2020-0091.
Senior Editor
- Didier YR Stainier, Max Planck Institute for Heart and Lung Research, Germany
Reviewing Editor
- Sigolène M Meilhac, Imagine-Institut Pasteur, France
Reviewers
- Sigolène M Meilhac, Imagine-Institut Pasteur, France
- Vincent M Christoffels, Amsterdam University Medical Centers, Netherlands
Version history
- Received: January 27, 2022
- Preprint posted: February 24, 2022 (view preprint)
- Accepted: April 28, 2023
- Version of Record published: May 15, 2023 (version 1)
Copyright
© 2023, Martin et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 498
- Page views
-
- 42
- Downloads
-
- 0
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Microbiology and Infectious Disease
Lateral partitioning of proteins and lipids shapes membrane function. In model membranes, partitioning can be influenced both by bilayer-intrinsic factors like molecular composition and by bilayer-extrinsic factors such as interactions with other membranes and solid supports. While cellular membranes can departition in response to bilayer-intrinsic or -extrinsic disruptions, the mechanisms by which they partition de novo are largely unknown. The plasma membrane of Mycobacterium smegmatis spatially and biochemically departitions in response to the fluidizing agent benzyl alcohol, then repartitions upon fluidizer washout. By screening for mutants that are sensitive to benzyl alcohol, we show that the bifunctional cell wall synthase PonA2 promotes membrane partitioning and cell growth during recovery from benzyl alcohol exposure. PonA2’s role in membrane repartitioning and regrowth depends solely on its conserved transglycosylase domain. Active cell wall polymerization promotes de novo membrane partitioning and the completed cell wall polymer helps to maintain membrane partitioning. Our work highlights the complexity of membrane–cell wall interactions and establishes a facile model system for departitioning and repartitioning cellular membranes.
-
- Cell Biology
- Computational and Systems Biology
mTORC1 senses nutrients and growth factors and phosphorylates downstream targets, including the transcription factor TFEB, to coordinate metabolic supply and demand. These functions position mTORC1 as a central controller of cellular homeostasis, but the behavior of this system in individual cells has not been well characterized. Here, we provide measurements necessary to refine quantitative models for mTORC1 as a metabolic controller. We developed a series of fluorescent protein-TFEB fusions and a multiplexed immunofluorescence approach to investigate how combinations of stimuli jointly regulate mTORC1 signaling at the single-cell level. Live imaging of individual MCF10A cells confirmed that mTORC1-TFEB signaling responds continuously to individual, sequential, or simultaneous treatment with amino acids and the growth factor insulin. Under physiologically relevant concentrations of amino acids, we observe correlated fluctuations in TFEB, AMPK, and AKT signaling that indicate continuous activity adjustments to nutrient availability. Using partial least squares regression modeling, we show that these continuous gradations are connected to protein synthesis rate via a distributed network of mTORC1 effectors, providing quantitative support for the qualitative model of mTORC1 as a homeostatic controller and clarifying its functional behavior within individual cells.