key: cord-0748102-we7p824s
authors: Hong, Lingjuan; Li, Na; Gasque, Victor; Mehta, Sameet; Ye, Lupeng; Wu, Yinyu; Li, Jinyu; Gewies, Andreas; Ruland, Jürgen; Hirschi, Karen K.; Eichmann, Anne; Hendry, Caroline; van Dijk, David; Mani, Arya
title: Prdm6 controls heart development by regulating neural crest cell differentiation and migration
date: 2022-02-02
journal: JCI insight
DOI: 10.1172/jci.insight.156046
sha: d09a3b56e0a112026aef8c1d9380e6a9934c9840
doc_id: 748102
cord_uid: we7p824s

The molecular mechanisms that drive the acquisition of distinct neural crest cell (NCC) fates is still poorly understood. Here, we identified Prdm6 as an epigenetic modifier that temporally and spatially regulates the expression of NCC specifiers and determines the fate of a subset of migrating cardiac NCCs (CNCCs). Using transcriptomic analysis and genetic and fate mapping approaches in transgenic mice, we showed that disruption of Prdm6 was associated with impaired CNCC differentiation, delamination, and migration and led to patent ductus arteriosus (DA) and ventricular noncompaction. Bulk and single-cell RNA-Seq analyses of the DA and CNCCs identified Prdm6 as a regulator of a network of CNCC specification genes, including Wnt1, Tfap2b, and Sox9. Loss of Prdm6 in CNCCs diminished its expression in the pre-epithelial–mesenchymal transition (pre-EMT) cluster, resulting in the retention of NCCs in the dorsal neural tube. This defect was associated with diminished H4K20 monomethylation and G1-S progression and augmented Wnt1 transcript levels in pre-EMT and neural tube clusters, which we showed was the major driver of the impaired CNCC migration. Altogether, these findings revealed Prdm6 as a key regulator of CNCC differentiation and migration and identified Prdm6 and its regulated network as potential targets for the treatment of congenital heart diseases.

Neural crest cells (NCCs) are a transient population of multipotent cells unique to vertebrates that arise from neural folds during embryonic development (1) and migrate throughout the body, giving rise to diverse cell lineages. In mice, cardiac NCCs (CNCCs) -precursors to the heart -are generated at 8-8.5 dpc and arrive in the distal outflow at E9-9.5 (2) . Targeted ablation of Wnt1-expressing CNCCs at 8.5 dpc in mice results in a complex phenotype of craniofacial and cardiovascular outflow tract defects (3) . Similarly, the impaired function of CNCCs in humans underlies the pathogenesis of various complex human congenital disorders collectively known as cardio craniofacial syndromes (4) . Cell lineage analyses have helped uncover the spatial and temporal diversification of NCC lineages (5, 6) . However, the mechanisms underlying the pleiotropic effects of CNCCs are not understood owing to a lack of insight into how subsets of CNCCs are specified and how their subsequent fate is determined.

PRDM6 is a smooth muscle cell-specific (SMC-specific) histone methyltransferase and a member of the PRDM family of transcriptional repressors. It is expressed in the cardiac outflow tract and the ductus arteriosus (DA), a small artery that connects the aorta and pulmonary artery (7) . Mice globally deficient for Prdm6 are embryonically lethal because of abnormal vascular patterning (8) . In humans, loss-of-function mutations in the PRDM6 gene have been the first, and to this date, the only known genetic cause of the familial nonsyndromic patent DA (PDA) (9) . PDA is an extreme example of an isolated cardiac disease that results from failure of closure of the DA, a derivative of the sixth pharyngeal arch, which is largely derived from premigratory CNCCs (10) . The tunica media of the DA and other pharyngeal arch-remodeled arteries

The molecular mechanisms that drive the acquisition of distinct neural crest cell (NCC) fates is still poorly understood. Here, we identified Prdm6 as an epigenetic modifier that temporally and spatially regulates the expression of NCC specifiers and determines the fate of a subset of migrating cardiac NCCs (CNCCs). Using transcriptomic analysis and genetic and fate mapping approaches in transgenic mice, we showed that disruption of Prdm6 was associated with impaired CNCC differentiation, delamination, and migration and led to patent ductus arteriosus (DA) and ventricular noncompaction. Bulk and single-cell RNA-Seq analyses of the DA and CNCCs identified Prdm6 as a regulator of a network of CNCC specification genes, including Wnt1, Tfap2b, and Sox9. Loss of Prdm6 in CNCCs diminished its expression in the pre-epithelial-mesenchymal transition (pre-EMT) cluster, resulting in the retention of NCCs in the dorsal neural tube. This defect was associated with diminished H4K20 monomethylation and G1-S progression and augmented Wnt1 transcript levels in pre-EMT and neural tube clusters, which we showed was the major driver of the impaired CNCC migration. Altogether, these findings revealed Prdm6 as a key regulator of CNCC differentiation and migration and identified Prdm6 and its regulated network as potential targets for the treatment of congenital heart diseases.

JCI Insight 2022;7(4):e156046 https://doi.org/10.1172/jci.insight.156046 to a lesser degree ZsGreen1-positive ductal SMCs in Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice were Ki67 positive (Supplemental Figure 2E ). Taken together, these findings indicated that the patency of the DA in Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice was not caused by impaired proliferation of NCC-derived SMCs. In addition, increased proliferation of non-NCC-derived SMCs appeared as insufficient to compensate for reduced NCC-derived SMCs in Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice. This finding is consistent with a prior report that NCC-derived and not mesodermal-derived SMCs are responsible for ductal closure (20) .

In addition, numerous ZsGreen1-positive cells were detected in the ventricular septum and to a lesser degree in the innermost layer of the left and right ventricles of WT pups at P0.5. The number of ZsGreen1-positive cells in the myocardium of Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice was also significantly diminished compared with WT mice (Supplemental Figure 2 , F and G), a finding that correlated with the ventricular noncompaction in these mice and suggested impaired NCC migration to the heart. Together, these results suggested that loss of Prdm6 expression resulted in impaired NCC migration and Loss of Prdm6 impairs proper migration and differentiation of CNCCs. The reduced number of ZsGreen1-positive cells in the DA of Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice at E17.5 and P0.5 ( Figure 3A and Supplemental Figure 2 , A and C) suggested impaired migration and/or differentiation of CNCCs. Considering the low levels of Prdm6 transcripts, we embarked on using scRNA-Seq analysis of CNCCs at E9.5 (21) . Cardiac NCCs were extracted by dissecting out Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 embryos and their WT littermates (Prdm6 +/+ Wnt1-Cre2 ZsGreen1) from mid otic placode to the posterior edge of the fourth somite ( Figure 4A ). ScRNA-Seq was carried out using the 10x Genomics Chromium system. The quality control (QC) threshold was made using a modified Seurat pipeline based on 2 QC covariates, count depth and the number of genes per barcode (see Methods). There were an average of 2742 and 3102 counts per cell or an average of 1294 and 1432 unique genes detected in WT and KO conditions, respectively. Only samples with RNA integrity number (RIN) greater than 6.5 were used (see Methods for more details). Graph-based clustering was performed using Seurat, and the PHATE (Potential of Heat-diffusion for Affinity-based Trajectory Embedding) technique developed by Moon and van Dijk et al. (22) was used to visualize trajectory structures. Cluster annotation was made manually based on specific cell markers (Supplemental Table 1 ).

PHATE-based analyses revealed 31 distinct clusters ( Figure 4B and Supplemental Figure 3 ). There was no obvious difference in the evolution of these clusters from pre-epithelial mesenchymal transition (preEMT)/EMT transition to migrating CNCCs and their terminal differentiation into cardiac and SMC lineages between WT and KO mice clusters. Although the total number of cells was similar between Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice and their WT littermates, the number of cells in each cluster varied between the 2 genotypes. Specifically, there were more cells in the neural tube and fewer cells in the EMT and migrating NCC clusters of Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice compared with their WT littermates (Prdm6 +/+ Wnt1-Cre2 ZsGreen1) ( Figure 4C ), consistent with our earlier findings from imaging ( Figure 2 , B and C). The single-cell transcript analysis showed that Prdm6 transcript levels were reduced in 2 clusters of the Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice compared with WT mice: 1 of 2 cardiac muscle lineages (cluster 8) and pre-EMT (cluster 10) ( Figure 4D ). Strikingly, in cluster 8, there was reduced expression of SMC markers Tagln and Actg2 and increased expression of endothelial genes and axonal guidance, such as Sema3c, Sema6d, Chd7, and CD34 (http://betsholtzlab.org/VascularSingleCells/database.html), in Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 versus WT littermates (Supplemental Figure 4) , suggesting impaired differentiation of CNCCs. This was consistent with our earlier observation of ZsGreen1 and VE-cadherin coexpression in the DA of Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice at P0.5 (Supplemental Figure 2 , A and B). Our findings collectively indicated that loss of Prdm6 in CNCCs resulted in the presence of fewer differentiated SMCs in the DA due to impaired migration and differentiation of CNCCs.

Smooth muscle-specific Prdm6-KO mice die postnatally from isolated PDA with altered differentiation of SMCs. To investigate the effect of Prdm6 disruption on vascular smooth muscle cell (VSMC) differentiation, we crossbred SM22-Cre mice with homozygote Prdm6 fl/fl mice to generate Prdm6 fl/+ SM22-Cre mice and then further crossbred to obtain Prdm6 fl/fl SM22-Cre mice. Homozygous Prdm6 fl/fl SM22-Cre mice showed no obvious growth retardation or syndromic features (Supplemental Figure 5 ), but all died at P1.5, slightly later compared with Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice. Gross anatomical examination of the cardiovascular system showed a PDA, albeit with a smaller diameter compared with Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice ( Figure 5A) . The cross-section examination of the DA showed complete patency of DA ( Figure 5B ), a completely penetrant trait in both male and female mice ( Table 2 ). We assessed cell proliferation in the DA at E17.5 using Ki67, when ductal cells were still viable in WT mice. Once again, we detected no significant Figure 5 , C and D). This was a striking finding given the patency phenotype and suggested that the loss of Prdm6 in VSMCs impairs their differentiation but not proliferation. Consistent with this, bulk RNA-Seq analysis of E17.5 DA from Prdm-6 fl/fl SM22-Cre mice and WT littermates identified contractile proteins, such as Tagln and Myh11, as the most downregulated genes in the DA of Prdm6 fl/fl SM22-Cre mice versus WT littermates ( Figure 5E and Supplemental Figure 6 ) from 280 differentially expressed genes (Supplemental Table 2 ). Strikingly, the expression levels of Tfap2b and Sox9, which also function as NCC specifiers, were significantly reduced in the DA of Prdm6 fl/fl SM22-Cre versus WT littermates ( Figure 5E and Supplemental Figure 6 ). Real-time quantitative PCR (qRT-PCR) further confirmed the reduced expression of Prdm6 ( Figure 6A ) in both the DA and aorta. However, the transcript levels of contractile proteins Myh11 and Tagln were reduced in the DA and not the descending aorta of Prdm6 fl/fl SM22-Cre mice at E17.5 ( Figure 6 , B and C). The expression levels of Tfap2b and Sox9 were also decreased, whereas those of the endothelial marker Kdr (also known as Flk1, Vegfr2) To identify the Prdm6-regulated gene network, we applied gene set enrichment analysis by using Gene Ontology (GO). Two separate analyses using genes significantly upregulated and downregulated were carried out. Vasculature development, blood vessel development, smooth muscle tissue development, collagen fibril organization, extracellular matrix assembly, and cell-matrix adhesion were the most significantly downregulated, and cell-cell adhesion and negative regulation of cell differentiation were among the most significantly upregulated pathways (Supplemental Table 3 ). We assessed the protein levels of selected differentially expressed genes by immunofluorescent staining of the DA of Prdm6 fl/fl SM22-Cre versus WT littermates at E17.5. We found that the SMC contractile protein Myh11 and fibronectin, a protein synthesized and secreted from SMCs, were significantly reduced in the DA of Prdm6 fl/fl SM22-Cre versus WT littermates (Supplemental Figure 7 , A and B). Fibronectin has been shown to be synthesized to a significantly greater degree in the DA compared with the aorta at early gestation and contributes to the enhanced migration of DA SMCs (23) .

Impaired CNCC migration in Prdm6 fl/fl mice results from deregulated Wnt1 signaling. The induction and subsequent specification of NCCs involves the combinatorial input of multiple signaling pathways. WNT and BMP signaling are key inducers of the neural crest, but their role in CNCC fate and migration has not been well studied. We found that expression levels of Wnt1 were highest in the pre-EMT cluster (cluster 10), followed by 1 of the 4 neural tube clusters (cluster 14) , and dramatically reduced after NCC delamination from the neural tube (Supplemental Figure 8A ). Wnt1 transcript levels were higher in Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 versus WT littermates (Wnt1-Cre2 ZsGreen1) in several clusters, including the pre-EMT cluster ( Figure 7A ), which suggested that Prdm6 may function to suppress Wnt1 during emergence of CNCCs. Pathway analysis of the pre-EMT cluster showed the "Development positive regulation of Wnt/β-catenin signaling at the receptor level" as one of the most significantly upregulated pathways (Supplemental Table 4 ).

Bmp signaling has been implicated in maintaining NCC identity by counteracting the effect of Wnt signaling (24) . Accordingly, we found reduced levels of Bmp4 in the neural tube, cardiac muscle lineages, and SMC clusters upon loss of Prdm6 compared with littermate controls ( Figure 7B , black arrowheads; Supplemental Figure 8B ), concomitant with increased Wnt1 levels. This indicated that Prdm6 controlled a balance between BMP/Wnt signaling, which was impaired in CNCCs of Prdm6 fl/fl mice and impeded their delamination and migration.

To establish the role of Prdm6 in NCC migration, we cultured excised neural tubes of the E9.5 WT and Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 embryos ex vivo to assess migration of CNCCs. We observed a marked reduction of NCC migration after 24 hours in Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 versus WT littermates (Figure 8, A and B) . Postulating that the increased Wnt1 expression underlies the CNCC migration defect, we examined whether the treatment of WT neural tube tissue ex vivo with Wnt1 could phenocopy the Prdm6 fl/fl condition. We therefore repeated these experiments using WT neural tubes and treated them with WNT1(100 ng/mL) for 24 hours. In response, we observed reduced migration of the delaminated NCCs compared with the PBS control, supporting our hypothesis (Figure 8 , C and D). These findings are consistent with a prior report that constitutive activation of β-catenin in the NCC-derived SMC precursor results in a reduced NCC-derived population of SMCs in the DA (20) , leading to its patency. Our findings also identified excess Wnt activation as a mechanism for PDA.

Loss of Prdm6 impairs H4K20 methylation and G1-S phase progression. We next focused on the molecular mechanisms that could link loss of Prdm6 to altered gene expression. We and others had previously shown that PRDM6 regulates H4K20 methylation (9, 25) . We compared H4K20 monomethylation of Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 mice CNCCs at E9.5 with those of the WT littermates ex vivo, and consistent with earlier findings, we observed reduced H4K20 monomethylation in the solitary leading edge of Prdm6 fl/fl and Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 (bottom row) mice at E17.5 stained for Ki67 (red) and αSMA (gray) show the localization of ZsGreen1 (green) and Ki67 (white arrowheads) and Ki67 in ZsGreen1-negative SMCs, indicated by yellow arrowheads. Scale bar: 50 μm. (E) Percentage of Ki67-positive cells in NCC-derived SMCs in WT mice and NCC-derived and non-NCC-derived SMCs in Prdm6 fl/fl Wnt1-Cre2 ZsGreen1 (KO) DA. All SMCs in WT mice were NCC derived. There was no significant difference in the percentage of Ki67-positive NCC-derived SMCs between WT and KO mice. The image intensities were quantified by ImageJ, and the thresholds for positive color detection were kept constant between different images. Each dot represents a biological replicate. The comparison between different groups was done by a 2-tailed unpaired t test, and data are shown as mean ± SEM. The comparisons between multiple groups (E) were done by 1-way ANOVA. A Mann-Whitney test was conducted for non-normally distributed data. ****P < 0.0001. All controls were the corresponding littermates (n = 5-10 per group).

JCI Insight 2022;7(4):e156046 https://doi.org/10.1172/jci.insight.156046

Wnt1-Cre2 ZsGreen1 CNCCs (Figure 8 , E-G). There was no significant difference of H4K20 monomethylation in early migrating NCCs (Supplemental Figure 9, A and B) . Together, these results suggest that loss of Prdm6 reduced H4K20 methylation in CNCCs.

Interestingly, H4K20 monomethylation (H4K20me1) has been shown to promote G1 to S cell cycle progression (26) , which is an essential step for the delamination of premigratory NCCs (27) . H4K20me1 at E-and N-cadherin promoters has also been shown to regulate their transcription and mediate EMT (28) . We therefore hypothesized that H4K20me1 by Prdm6 promotes cell cycle progression and alters the expression of a gene network that is necessary for delamination and migration of the CNCCs. To test this, we used CRISPR/Cas9-based gene editing to inactivate Prdm6 in exon2 of O9 mouse NCCs in vitro. We first differentiated these cells into SMCs ( Figure 9A ) and verified their competency in expressing Prdm6, contractile proteins Myh11 and Tagln, and NCC specifiers Sox9 and Tfap2b at days 0, 1, 3, and 6 ( Figure  9B ). Subsequently, we infected the cells with lentivirus containing CRISPR/Cas9 and transfected them with either control vector or 3 different sgRNAs targeting Prdm6 exon2 ( Figure 9C) and generated 3 clonal cell lines with 8, 55, and 56 bp deletions ( Figure 9D ). qRT-PCR expression analysis confirmed the loss of Prdm6 expression and subsequent downregulation of contractile proteins Myh11 and Tagln and NCC specifiers Tfap2b and Sox9, as well as upregulation of the endothelial maker Kdr in all Prdm6-deficient lines (Figure 9 , E-J), consistent with our in vivo gene expression analyses. We then subjected the Prdm6-deficient CNCCs to FACS analysis to assess the cell cycle activity and observed a marked reduction in G1-S cell cycle progression compared with WT O9 cells ( Figure 10A ). Our analyses revealed a greater number of Prdm6-deficient CNCCs in the G1 phase and fewer in the S and G2/M phase compared with WT CNCCs, indicating G1-S phase arrest. This finding was associated with increased Wnt1 expression in Prdm6-KO versus WT O9 cells ( Figure 10B ).

To determine the direct role of H4K20me1 in regulating CNCC cell cycle progression, we inactivated retinoblastoma-associated protein Rb1 in O9 cells using Rb1-specific siRNA. RB1 protein plays a critical role in recruiting H4K20 methyltransferases, such as Suv4-20h1/KMT5B and Suv4-20h12/ KMT5C29, and mediates mono-, di-, and tri-methylation of H4K20 (29) . Knockdown of Rb1 in O9 cells resulted in modest but persistently reduced G1-S progression compared with control siRNA, as determined by FACS analysis ( Figure 10C ). There was an inverse relationship between G1-S progression and Wnt1 expression (R = 0.83; Figure 10C ).

Our data suggested that loss of Prdm6 reduced H4K20me1 and increased Wnt1 levels, impeding CNCC migration, but the order of the events remained unknown. We examined the effect of Wnt1 on histone modifications in these cells and found that it increased H4K20me1 levels compared with PBS, as predicted ( Figure 10 , D and E). This finding suggests that reduced H4K20me1 after Prdm6 knockdown is not caused by increased Wnt1. Accordingly, treatment of Prdm6-deficient O9 cells with Wnt inhibitor (IWP-4, 5 μM, 04-0036) rescued the migration-defective phenotype (Figure 10 , F and G); however, it did not change H4K20me1 ( Figure 10H ). In conclusion, we suggest a model whereby Prdm6-H4K20me1 synchronizes the cell cycle, which is necessary for the regulation of Wnt1 transcript levels and CNCC migration. In contrast, the loss of Prdm6 results in reduced H4K20 monomethylation and G1-S cell cycle progression, resulting in increased Wnt1 expression and impaired CNCC migration.

NCCs give rise to a diverse cell lineage and therefore are increasingly gaining attention for their potential role in tissue regeneration. The mouse phenotyping and fate mapping in our study showed that Prdm6 regulates differentiation and migration of CNCCs, which are the precursor of SMCs in the tunica media of the DA and pharyngeal arches (16) . Consequently, loss of Prdm6 results in early postnatal death from PDA (30) (31) (32) . The importance of CNCCs in DA closure has been previously shown by the ablation of myocardia in early postmigratory NCC derivatives, albeit the exact molecular mechanisms remain to be elaborated (32) . Mutations in histone modifiers KMT2D (histone-lysine N-methyltransferase 2D) and KDM6A (lysine-specific demethylase 6A) (33) and EHMT1 (euchromatic histone-lysine N-methyltransferase 1) (34) have also been associated with syndromic congenital heart disease due to defects in NCC development (35) . However, in contrast to PRDM6 mutations, defects in these genes appear to affect both cardiac and cranial NCCs.

(red) and costained with antibodies against αSMA (green), CD31 (gray), and DAPI (blue). Scale bar: 50 μm. (D) Quantification of Ki67-positive cells in the smooth muscle cell layer. (E) Volcano plot of bulk RNA-Seq data of Prdm6 fl/fl SM22-Cre mice DA compared with control mice at E17.5. The log fold change and the log P values are shown on x and y axes, respectively. The comparison between different groups was done by a 2-tailed unpaired t test, and data are shown as mean ± SEM. The significantly changed transcripts in Prdm6 fl/fl SM22-Cre versus control mice are shown in red and the locations of specific genes of interest are shown. The significantly altered transcripts were determined using the following thresholds: adjusted P < 0.05, unpaired t test, FDR < 0.05, Benjamini-Hochberg corrected. All controls were the corresponding littermates (n = 4 per group). Our investigations ex vivo demonstrated that increased Wnt1 underlies impaired CNCC migration in Prdm6-deficient mice, which was rescued by Wnt inhibition. These findings are consistent with a prior report that constitutive activation of β-catenin (ctnnb1Δex3) results in a reduced NCC-derived population of SMCs in the DA and early postnatal death (20) . Although NCC migration in ctnnb1Δex3 was not explored, our fate mapping and functional studies suggest that excessive Wnt activation can result in impaired CNCC migration.

Our comprehensive genomic and epigenomic approach identified a Prdm6-dependent gene regulatory network that underlies the regulation of a subset of CNCCs during development. The neural crest gene regulatory networks consist of multiple complex subnetworks that are sequentially activated to control the key steps of neural crest formation, delamination, migration, and differentiation. Of particular importance is the mutually inhibitory crosstalk between Wnt and Bmp signaling in that Bmp counteracts proliferation promoted by Wnt while Wnt antagonizes BMP-dependent neuronal differentiation (36) . The coordinated actions of these 2 pathways result in activation of transcriptional cascades that drive EMT and migration and cooperate with signaling systems that induce cell differentiation (37) . In our models, we showed that the absence of Prdm6 tips the balance in favor of Wnt1 by reducing Bmp4 and increasing Wnt1 transcripts in several clusters, leading to impaired CNCC migration.

PRDM6 has been previously shown by our group and others to mediate H4K20 methylation in vitro (9, 25) . H4K20 monomethylation has been shown to be critical for cell cycle progression (38, 39) . In the current study, we showed that Prdm6 regulates H4K20 methylation in migrating NCCs and that this event is critical for cell cycle progression and subsequent induction of EMT and migration of specific CNCCs. We showed that H4K20 methylation and cell cycle activation are inversely related to Wnt1 expression. The causal relationship between H4K20 methylation, cell cycle activation, and Wnt1 expression was supported by knocking down the retinoblastoma-associated protein Rb1, which has been shown to interact with H4-K20 methyltransferases Suv4-20h1/KMT5B and Suv4-20h12/ KMT5C (40) and to mediate their methylation of H4K20 (41) . Accordingly, the knockdown of Rb1 led to G1-S phase arrest and was associated with higher Wnt1 expression. These findings linked reduced H4K20 methylation in Wnt1-Cre2 Prdm6 fl/fl mice to increased Wnt1 expression in multiple CNCC clusters. One of the most important implications of our study is the role of temporal and spatial epigenetic regulation in the specification of NCCs. It is likely in this context that elevated Wnt1 impairs a subset of CNCCs that contributes to the DA, whereas elevated Wnt1 is a necessary component of trunk NCCs, including their entry into the S-phase and their facilitated migration (42) . Strikingly, environmental factors, such as teratogens and infectious disorders during the first trimester of pregnancy, have been associated with isolated persistent PDA and altered DNA or histone modifications (43) (44) (45) (46) . The characterization of Prdm6 function in NCCs has, therefore, a broad implication for the pathogenesis of inherited and acquired congenital heart diseases and may lead to identification of novel targets for drug development. 

Generation of SM22-Cre KO mice. The day of plug detection was considered as E0.5. Prdm6-flox mice were generated in-house and were intercrossed to generate homozygous Prdm6-flox mice (Prdm6 fl/fl ). Prdm6 fl/fl mice were viable, fertile, and indistinguishable from control littermates. Prdm6 fl/fl mice were intercrossed with SM22-Cre transgenic mice that expressed high levels of Cre-recombinase in SMCs (Supplemental Figure 10) .

Neural crest-restricted inactivation of the Prdm6 gene. To determine the function of Prdm6 in cardiac neural crest-derived cells, homozygous Prdm6 fl/fl mice were crossed with transgenic mice expressing Cre-recombinase under the transcriptional control of the neural crest-restricted Wnt1 promoter. Prdm6 fl/+ Wnt1-Cre2 offspring were then interbred with Prdm6 fl/+ mice to generate Prdm6 fl/fl Wnt1-Cre2 mutants and control littermates Prdm6 +/+ Wnt1-Cre2 (Supplemental Figure 11) .

For lineage tracing, Prdm6 fl/fl mice were crossed with ZsGreen1 to generate Prdm6 fl/+ ZsGreen1 mice and then interbred with ZsGreen1 to generate Prdm6 fl/+ ZsGreen1/ZsGreen1 mice. To knock out Prdm6 in Wnt1-positive NCCs, the Prdm6 fl/+ ZsGreen1/ZsGreen1 mice were intercrossed with Prdm6 fl/+ Wnt1-Cre2 mice (Supplemental Figure 12) .

Wnt1-Cre2 mice were purchased from The Jackson Laboratory, stock number 022501 (full strain name B6.Cg-E2f1Tg(Wnt1-cre)2Sor/J). ZsGreen1 mice were received from Daniel Greif (Yale University) and are also available from The Jackson Laboratory (Rosa-CAG-LSL-ZsGreen1-WPRE, Jax 007906). All mice were fed ad libitum and housed at constant ambient temperature in a 12-hour light/12-hour dark cycle. Both male and female mice were studied. All studies in animals were conducted in accordance with the NIH Guide for the Care and Use of Laboratory Animals (National Academies Press, 2011).

Mouse embryo and tissue preparation and IHC. Mouse embryos and tissues were collected and rinsed in ice-cold Dulbecco's PBS and fixed by incubation for 15-20 minutes in PBS containing 4% PFA (Santa Cruz, sc-281692) at 4°C. They were rinsed twice in ice-cold PBS and incubated for 24 hours in 30% sucrose/PBS at 4°C, followed by Tissue-Tek OCT compound (catalog 4583) in plastic molds and allowed to freeze on dry-ice and stored at -80°C. The tissue submerged in OCT compound was sectioned onto 5 μm-10 μm slices in a cryostat and dried on Fisherbrand Superfrost Plus slides (Fisher Scientific, catalog 22-037-246) 12 to 16 hours at room temperature. P0.5 mice tissue were collected after anesthesia. Hearts were then removed and treated as described above.

Immunofluorescence. All immunofluorescence staining was performed on 5 μm-10 μm frozen sections. For immunostaining, the slides were washed 3 times in 1× PBS, followed by 3 washes in PBS plus 0.1% Triton X-100 and incubated in the blocking solution containing 5% BSA plus 1× PBS plus 0.1% Triton X-100 for 2-3 hours. Subsequently, the slides were incubated in the primary antibody overnight, followed by 3 washes in PBS plus 0.1% Triton X-100 for 10 minutes each. After 2 hours of secondary antibody incubation, slides were washed 3 times in PBS plus 0.1% Triton X-100 for 10 minutes each. The specimens were overlaid with ProLong-Gold Antifade solution and a coverslip. All primary antibodies (Supplemental Table 5 ) were diluted 1:300 and secondary antibodies at 1:300. As a control, we used no primary as well as PRDM6 KO mice tissues. Otherwise, only validated and previously published antibodies were used. Fluorescence images were obtained by a Zeiss 4 laser confocal microscope or Leica SP8 confocal microscope, and the intensity was measured using the same laser output, gain, and offset for each set of antibodies tested. Images were quantified with ImageJ (NIH) and adjusted for the area.

Whole-mount staining. Whole-mount microscopy was carried out as previously described (47) . Briefly, E9.5 mice were harvested and kept in ice-cold PBS for 5 minutes to drain blood from the umbilical cord and then fixed in 4% PFA at 4°C and washed at room temperature. Samples were dehydrated in methanol/PBS series (20%, 40%, 60%, 80%) for 30 minutes at room temperature and then 100% for the entire day. Samples were chilled over ice and bleached in 5% H 2 O 2 in methanol overnight at 4°C, washed with 80% methanol/PBS for and KO (Prdm6 fl/fl Wnt1-Cre2 ZsGreen1) mice at E9.5. The NCCs are shown as ZsGreen1-positive (green) and the neural tube is stained with DAPI (blue) and phalloidin (red). Scale bar: 500 μm. (B) Quantification of migration distances; each dot represents a biological replicate. Migration is measured as distance from the periphery at the junction with 2 perpendicular radiuses. (C) Migration of NCCs from ex vivo cultured neural tube of WT mice at E9.5 treated with WNT1 for 24 hours and stained with DAPI (blue), Tfap2b antibody (red), phalloidin (green). Scale bars: 500 μm. (D) Quantification data of C. Each dot represents a biological replicate. (E) H4K20 monomethylation (red) of NCCs in the leading edge of ex vivo cultured neural tube of WT (Prdm6 +/+ Wnt1-Cre2 ZsGreen1) and KO (Prdm6 fl/fl Wnt1-Cre2 ZsGreen1) mice at E9.5 stained with DAPI (blue) and phalloidin (gray). Scale bar: 50 μm. (F) Schematic depiction of ZsGreen1-positive (green) migrating NCCs, divided into 2 populations of migrating (clustered) and leading edge (single cells). (G) Quantification of H4K20 monomethylation from E; each dot represents a cell. The comparison between different groups was done by a 2-tailed unpaired t test, and data are shown as mean ± SEM. *P < 0.05, ****P < 0.0001. All controls were the corresponding littermates (n = 3 per group). 30 minutes at room temperature, and rehydrated with the methanol series in PBS/0.2% Triton X-100 (fresh prepared): 80%, 60%, 40%, 20%, 0%; 30 minutes each at room temperature. For immunolabeling, pretreated samples were incubated in 1× PBS/0.2% Triton X-100/20% DMSO/0.3 M glycine at 37°C for 24 hours and blocked in 1× PBS/0.2% Triton X-100/10% DMSO/6% goat serum at 37°C for 24 hours and subsequently washed with 1× PBS/0.2% Tween 20 with 10 μg/mL heparin and incubated with primary antibody in 1× PBS/0.2% Tween 20 with 10 μg/mL heparin/5% DMSO/3% donkey serum at 37°C for 24 hours. After several washes, samples were incubated with secondary antibody in 1× PBS/0.2% Tween 20 with 10 μg/mL heparin/3% donkey serum at 37°C for 24 hours. For clearing, samples were incubated overnight in 50%, 80%, and 100% tetrahydrofuran/H 2 O (THF, Sigma-Aldrich, 186562-12X100ML) for 2 hours each time. Samples were then incubated in dichloromethane (DCM, Sigma-Aldrich, 270997-12X100ML) until they sunk to the bottom for 5 minutes and in dibenzyl ether (Sigma-Aldrich, 108014-1KG) until the samples were clear (20 minutes to 2 hours). We then examined samples using light sheet microscopy.

Isolation and culture of the neural tube. Isolation of the neural tube was carried out as described by Elise R. Pfaltzgraff, Nathan A. Mundell, and Patricia A. Labosky (48) . For timed pregnancies, dams were checked in the morning for a vaginal plug and the embryo was considered as E0.5 at noon of the day the plug was observed. Mice were euthanized, and the embryos were removed at E9.5 and placed on ice-cold sterile DPBS. Using insulin needles, neural tubes were cut between the mid otic placode and the posterior edge of the fourth somite. The yolk sacs were kept for genotyping and the remaining embryonic tissues were discarded. The cardiac neural tube-containing segments were placed in collagenase/Dispase (Sigma-Aldrich, catalog 269-638) at room temperature for 8 minutes and subsequently washed in ice-cold sterile DPBS. The non-neural ectoderm was gently removed from the tissue, and the somites were separated away from the cardiac neural tube. The remaining parts of the mesoderm were gently removed using a cut-down 1 mL pipette tip. The isolated neural tube was washed for 30 seconds with culture medium (ES101-B) with fibroblast growth factor (GF003) and placed in the center of Matrigel-coated (Corning, catalog 354234) 24-well glass-bottom cell culture plate (Nest Scientific, catalog 801006) and incubated at 37°C.

RNA extraction, cDNA synthesis, and qPCR. The RNA was extracted using Qiagen miRNeasy Micro kit (catalog 1071023). Next, 100-500 ng of RNA was used for cDNA synthesis (iScript cDNA synthesis, Bio-Rad, 1708890). No-RT and no-RNA controls were used to ensure DNase digestion and no residual contamination. The cDNA was diluted to different concentrations and eventually were diluted by 5-fold to obtain qPCR values in the linear range. qPCR primers (Supplemental Table 6 ) were from PrimerBank. A dilution curve that included primers and qPCR Master Mix were utilized to ensure that PCR primers did not self-anneal and to reduce background amplification. The qPCR was performed with 2× SYBR Green Master Mix (Bio-Rad, 170-8880), 0.5 μM primer (forward and reverse), and 5 μL of 1:5 diluted cDNA and rest water. All qPCR reactions had 3 technical replicates and at least 3 biological replicates per genotype. All values were normalized to the reference gene Gapdh.

Bulk RNA-Seq. The reads were trimmed for quality and length. All the reads were minimum base quality of Q30, and reads with minimum length of 45 were used for further analysis. The reads were aligned to the mouse UCSC reference genome mm10 (49, 50) using TopHat2 (51) . The alignment data were converted to per-gene counts using Cufflinks (52) and further analyzed using Cuffdiff and R (http://www.r-project.org/index.html). The results were visualized using R.

ScRNA-Seq. E9.5 embryos were harvested and briefly placed in ice-cold DPBS and then placed on a plastic petri dish resting on ice. Using insulin needles, the cardiac neural crest was dissected out from the mid otic placode to the posterior edge of the fourth somite. The tissues were removed and placed in a 1.5 mL conical snap-cap plastic tube on ice and 1× enzyme digestion buffer containing 100 μL of Bacillus licheniformis stock (100 mg/mL, Sigma-Aldrich, P5380), 895 μL of DPBS, and 5 μL of 0.5 M EDTA. The tube was quickly shaken thoroughly to mix and resuspend, incubated in ice-water bath, and triturated using a cut-down 1 mL pipette tip for 15 seconds every 2 minutes. Cells were filtered through a 40 μm strainer, rinsed with 10 mL of ice-cold PBS/BSA, and weighed. Subsequently, 300 g of the cells were transferred to 15 mL conical tube, centrifuged at 300g for 5 minutes at 4°C, and the supernatant discarded. The pelleted cells were then resuspended in 100 μL of ice-cold DPBS and counted on a Countess II (Thermo Fisher Scientific). Cell count and viability were determined using trypan blue dye exclusion.

QC. QC was evaluated from the distribution of QC covariates in data set visualizations and clustering after multiple times of analyzing the data and comparison with bulk RNA-Seq data. The analysis showed that use of fraction of counts from mitochondrial genes per barcode was too stringent and resulted in significant JCI Insight 2022;7(4):e156046 https://doi.org/10.1172/jci.insight.156046 loss of signal since cardiovascular progenitor cells and migrating NCCs had a high degree of mitochondrial content. A more permissive QC threshold based on only 2 QC covariates, count depth and the number of genes per barcode, was considered as more appropriate. Capture sites that were either empty or contained multiple cells were filtered. DoubletFinder (53) and high-count depth threshold were also used to exclude doublets. Normalization protocol included count depth scaling (AKA CPM normalization). Gene normalization, which constitutes scaling gene counts to have zero mean and unit variance (z scores), was applied. Normalized data were log-transformed for use with downstream analysis. Data correction for biological effects based on cell cycle signals was avoided because of the presence of a large number of proliferating cells.

Single-cell sequencing pipeline. Ten thousand cells per sample were used for the 10x Genomics Cell Ranger pipeline with a combined human genome and the ZsGreen1 gene to obtain count matrices for the WT and KO conditions. There were an average of 2742 and 3102 counts per cell or an average of 1294 and 1432 unique genes detected in WT and KO conditions, respectively. Only samples with RIN greater than 6.5 were used. Library preparations and sequencing were performed by the Yale Center for Genome Analysis. Briefly, gel beads in emulsion (GEMs) were created using Chromium Next GEM Single Cell 3′ Gel Beads v3.1 kit (10x Genomics) and barcoded. GEMs generated were used for cDNA synthesis and library preparation using the Chromium Single Cell 3′ Library Kit v3.1 (10x Genomics) and sequenced on NovaSeq 6000 system using HiSeq 100 base pair reads and dual indexing.

ScRNA-Seq data alignment. scRNA-Seq samples were demultiplexed using MULTI-Seq and barcoded and pooled on droplet microfluidic emulsion and sequenced with the 10x Genomics Chromium system. Cell Ranger software v1.3.1 was used for library demultiplexing, FASTQ file generation, read alignment, and unique molecular identifier (UMI) quantification. The resulting raw UMI counts in each cell were normalized to their library size. The normalized counts were square-root transformed. We then used a combination of 3 metrics to find differentially expressed genes across conditions: (a) the binary logarithm of fold change between mean counts; (b) the Wasserstein or earth mover's distance (EMD); and (c) adjusted P value from a 2-sided Mann-Whitney U test with continuity and Benjamini-Hochberg correction, as described by Moon et al. (22) . A P value less than or equal to 0.01 was used as the threshold for statistical significance. The top 30 differentially expressed genes (up-or downregulated, ranked by Wasserstein distance, were represented in heatmaps. To identify putative cellular functions changed across conditions, we performed PANTHER-GO gene set enrichment analysis for each cluster (54) using the default human PANTHER-GO reference list as a background.

scRNA-Seq data analysis. We used the standard scRNA-Seq analysis pipeline for clustering (55) . We removed cells with less than that 200 and more than 6000 raw UMIs to account for doublets. Next, we removed genes that were expressed in fewer than 3 cells and removed cells that expressed fewer than 200 genes. We also removed UMIs pertaining to mitochondrial and ribosomal genes. The counts in each cell were normalized to their library size. Then, normalized counts were square-root transformed, similar to a log transform but not requiring addition of a pseudo count. Data preprocessing was performed in Python (v3.7.4) using Scanpy (v1.4.6) (56, 57).

We visually observed batch effects between conditions in 2D cellular embeddings. To remove these batch effects for downstream analysis, we used an approximate batch-balanced k nearest neighbor (kNN) graph for manifold learning (BB-kNN batch-effect correction) using Scanpy's fast approximation implementation (56, 57) . For each cell, the 3 nearest neighboring cells in each condition were identified by Euclidean distance in 100-dimensional principal component analysis (PCA) space. This kNN graph was used as the basis for downstream analysis.

To visualize the scRNA-Seq data, we implemented 2 nonlinear dimension reduction methods and used the BB-kNN batch-corrected connectivity matrix as input for uniform manifold approximation and projection (UMAP) (58) and PHATE (22) . UMAPs were generated using a minimum distance of 0.5. PHATE projections were generated with a gamma parameter of 0.

For cell clustering, we used the Louvain community detection method (59) with the BB-kNN graph. We used high-resolution community detection and merged clusters based on expression of hand-picked markers for the cell types (Supplemental Figure 13) . In order to minimize the effect of dropout for clustering, we imputed counts using Markov affinity-based graph imputation of cells (MAGIC) (60).

To find the differentially expressed genes between the WT and KO conditions in each cell type, we used 3 metrics: the 1-dimensional Wasserstein or EMD, an adjusted P value from a 2-sided Mann-Whitney U test with continuity and Benjamini-Hochberg correction, and the binary logarithm of fold-change between mean counts. Significance was set to an adjusted P value of 0.01. The EMD can be defined as

the minimal cost to transform distribution to another and has been previously used to perform differential gene expression analysis

Single-cell suspensions were prepared and adjusted to contain 2 × 10 6 cells/mL. Single-cell suspensions were stained with Hoechst (10 μg/mL) and pyronin Y for 15 minutes at 37°C. Statistics. In vivo studies included a minimum of 3 mice in each group. In vitro studies were carried out in more than 3 independent experiments. The comparison between different groups was done by a 2-tailed unpaired t test. The normality was tested by Kolmogorov-Smirnov test. An F statistic was calculated to determine whether variances were different between samples, and the P values were then corrected with Welch's 2-tailed t test. The comparisons between multiple groups were done by 1-way ANOVA. A Mann-Whitney test was conducted for non-normally distributed data. Fisher's exact test was carried out for the continuous variables. The preparation of graphs and all statistical analyses, including 2-tailed Student's t tests, 2-way ANOVA (SigmaPlot), and testing for equal variance, were carried out using GraphPad Prism 8.1 Project software. P values less than 0.05 were considered significant. Data are presented as mean ± SEM. Fluorescence images were evaluated using ImageJ (NIH). Code and data availability. All differential gene expression analyses from bulk and scRNA

Neural crest cell plasticity and its limits

In the beginning: generating NCC diversity

Temporal-spatial ablation of neural crest in the mouse results in cardiovascular defects

The neural crest in cardiac congenital anomalies

Spatiotemporal structure of cell fate decisions in murine neural crest

Neural crest lineage analysis: from past to future trajectory

Influence of white-coat hypertension on left ventricular deformation 2-and 3-dimensional speckle tracking study

Prdm6 is essential for cardiovascular development in vivo

Mutations in the histone modifier PRDM6 are associated with isolated nonsyndromic patent ductus arteriosus

Pathogenesis of persistent truncus arteriosus and dextroposed aorta in the chick embryo after neural crest ablation

Fate of the mammalian cardiac neural crest

Novel TFAP2B mutations that cause Char syndrome provide a genotype-phenotype correlation

Syndromic patent ductus arteriosus: evidence for haploinsufficient TFAP2B mutations and identification of a linked sleep disorder

Fate of the mammalian cardiac neural crest

The heart of the neural crest: cardiac neural crest cells in development and regeneration

Neural crest cell contribution to the developing circulatory system: implications for vascular morphology?

Apoptosis and regulation of Bax and Bcl-X proteins during human neonatal vascular remodeling

Role of apoptosis in the closure of neonatal ductus arteriosus

Differentiation, dedifferentiation, and apoptosis of smooth muscle cells during the development of the human ductus arteriosus

A subpopulation of smooth muscle cells, derived from melanocyte-competent precursors, prevents patent ductus arteriosus

Neural crest cells in cardiovascular development

Visualizing structure and transitions in high-dimensional biological data

Fibronectin, hyaluronan, and a hyaluronan binding protein contribute to increased ductus arteriosus smooth muscle cell migration

Neural crest stem cell maintenance by combinatorial Wnt and BMP signaling

PRDM6 is enriched in vascular precursors during development and inhibits endothelial cell proliferation, survival, and differentiation

Suv4-20h1 promotes G1 to S phase transition by downregulating p21 WAF1/CIP1 expression in chronic myeloid leukemia K562 cells

Association between the cell cycle and neural crest delamination through specific regulation of G1/S transition

SET8 promotes epithelial-mesenchymal transition and confers TWIST dual transcriptional activities

Role of the RB1 family in stabilizing histone methylation at constitutive heterochromatin

The prostaglandin receptor EP4 triggers remodelling of the cardiovascular system at birth

Genetic model of selective COX2 inhibition reveals novel heterodimer signaling

Myocardin regulates expression of contractile genes in smooth muscle cells and is required for closure of the ductus arteriosus in mice

Exome sequencing identifies MLL2 mutations as a cause of Kabuki syndrome

Loss-of-function mutations in euchromatin histone methyl transferase 1 (EHMT1) cause the 9q34 subtelomeric deletion syndrome

Epigenetic regulation of cardiac neural crest cells

Wnt/BMP signal integration regulates the balance between proliferation and differentiation of neuroepithelial cells in the dorsal spinal cord

Regulatory logic underlying diversification of the neural crest

SETting the clock for histone H4 monomethylation

Monomethylation of histone H4-lysine 20 is involved in chromosome structure and stability and is essential for mouse development

The retinoblastoma protein regulates pericentric heterochromatin

The SUV4-20 inhibitor A-196 verifies a role for epigenetics in genomic integrity

Canonical Wnt activity regulates trunk neural crest delamination linking BMP/noggin signaling with G1/S transition

Associative DNA methylation changes in children with prenatal alcohol exposure

Alteration of gene expression, DNA methylation, and histone methylation in free radical scavenging networks in adult mouse hippocampus following fetal alcohol exposure

Epigenetic landscape of amphetamine and methamphetamine addiction in rodents

Prenatal antiepileptic exposure associates with neonatal DNA methylation differences

Whole-mount three-dimensional imaging of internally localized immunostained cells within mouse embryos

Isolation and culture of neural crest cells from embryonic murine neural tube

The human genome browser at UCSC

GENCODE reference annotation for the human and mouse genomes

TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions

Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks

DoubletFinder: doublet detection in single-cell RNA sequencing data using artificial nearest neighbors

Pathway enrichment analysis and visualization of omics data using g:Profiler, GSEA, Cytoscape and Enrich-mentMap

Spatial reconstruction of single-cell gene expression data

SCANPY: large-scale single-cell gene expression data analysis

BBKNN: fast batch alignment of single cell transcriptomes

Uniform Manifold Approximation and Projection for Dimension Reduction

Fast unfolding of communities in large networks

Recovering gene interactions from single-cell data using data diffusion

SigEMD: a powerful method for differential gene expression analysis in single-cell RNA sequencing data

Earth mover's distance (EMD): a true metric for comparing biomarker expression levels in cell populations

Single-cell longitudinal analysis of SARS-CoV-2 infection in human airway epithelium identifies target cells, alterations in gene expression, and cell state changes

This work was supported by grants from NIH (5R35HL135767A to AM). The authors acknowledge the use of the shared imaging resources from the Yale Translational Research Imaging Center (Y-TRIC).Address correspondence to: Arya Mani, Yale Cardiovascular Research Center, 300 George Street, New Haven, Connecticut, 06511, USA. Email: arya.mani@yale.edu.

LH designed and conducted experiments, acquired data, analyzed data, and helped with preparing the manuscript. NL generated floxed mice. VG analyzed scRNA-Seq data. SM analyzed bulk RNA-Seq data. LY, YW, and JL analyzed data. AG and JR provided the mice. KKH and AE supervised the data analysis. CH helped with the writing of the manuscript and provided instruction on the preparation of the data. DVD supervised the analysis of scRNA-Seq data. AM designed and supervised the research studies and wrote the manuscript.