The cardiac outflow tract (OFT) transiently links the ventricles to the aortic sac and forms the arterial valves. Abnormalities in these valves, such as bicuspid aortic valve (BAV), are common congenital anomalies. GATA6-inactivating variants cause cardiac OFT defects and BAV, but their mechanisms are unclear. We generated Gata6STOP/+ mice using CRISPR-Cas9, which show highly penetrant BAV (70%) and membranous ventricular septal defects (43%). These mice exhibited decreased proliferation and increased ISL1-positive progenitor cells in the OFT, indicating abnormal cardiovascular differentiation. Gata6 deletion with the Mef2cCre driver line recapitulated Gata6STOP/+ phenotypes, indicating a cell-autonomous role for Gata6 in the second heart field. Gata6STOP/+ mice showed reduced OFT length and caliber, associated with deficient cardiac neural crest cell contribution, which may cause valvulo-septal defects. RNA-sequencing analysis showed depletion in pathways related to cell proliferation and migration, highlighting Cxcr7 (also known as Ackr3) as a candidate gene. Reduced mesenchymal cell migration and invasion were observed in Gata6STOP/+ OFT tissue. CXCR7 agonists reduced mesenchymal cell migration and increased invasion in wild-type but not in Gata6STOP/+ explants, indicating the GATA6-dependent role of CXCR7 in OFT development and its potential link to BAV.

Congenital heart defects are the commonest birth defects (Hoffman and Kaplan, 2002; van der Linde et al., 2011). Up to 30% of all congenital heart defects involve the outflow tract (OFT) (Srivastava and Olson, 2000; Thom et al., 2006), which briefly connects the embryonic ventricles to the aortic sac at the arterial pole of the developing heart (Rochais et al., 2009). OFT formation relies on the timely contribution of progenitor cell populations from within and outside the cardiac developmental field. Defective OFT development leads to the misalignment of the aorta and pulmonary trunks, and failure to separate systemic and pulmonary circulations (Conway et al., 2003). The endocardial cushions in the OFT are formed via the epithelial–mesenchymal transition of the endocardium and the invasion of extra-cardiac neural crest cells (cNCCs) (Snarr et al., 2008). The latter migrate into the caudal pharyngeal arches to septate the OFT and contribute to valve leaflet patterning, giving rise to the arterial (aortic and pulmonary) valves (Lin et al., 2012; Phillips et al., 2013; Snarr et al., 2008).

Bicuspid aortic valve (BAV) is the commonest congenital heart defect, affecting up to 0.5 to 2% of the population (Masri et al., 2017). BAV is characterized by the presence of two asymmetrical aortic valve leaflets instead of three symmetrical ones (Henderson et al., 2022). BAV disrupts normal aortic valve flow patterns, resulting in abnormal turbulence and increased tissue stress (Ladich et al., 2011). Patients with BAV often develop premature calcific aortic valve disease, leading to valve stenosis and aortic valve insufficiency that requires valve replacement (Della Corte et al., 2007; Prakash et al., 2014; Roth et al., 2020; Tastet et al., 2017; Yuan and Jing, 2010).

The genetic architecture of BAV is complex, with low-penetrance variants, variable clinical expressivity and locus heterogeneity (Giusti et al., 2017), suggesting the involvement of multiple gene–gene and gene–environment interactions. Genome-wide association studies have identified several chromosomal regions, such as 18q, which harbors the gene encoding the GATA6 transcription factor, among several genes associated with BAV (Martin et al., 2007). In mice, Gata6 is expressed in the myocardium, endocardium, cardiac neural crest and vascular smooth muscle (Morrisey et al., 1996), suggesting diverse functions in heart morphogenesis that remain to be fully elucidated.

GATA6 variants lead to a variety of cardiac phenotypes ranging from structural malformations to conduction defects (Lepore et al., 2006; Škorić-Milosavljević et al., 2019). GATA6 variants have been associated with familial BAV (Xu et al., 2018), but have been mostly found in sporadic cases (Alonso-Montes et al., 2018; Ma et al., 2021). The phenotypes of mice lacking Gata6 highlight essential roles in organismal development. Gata6-null mice are early [embryonic day (E) 5.5] embryonic lethal due to extra-embryonic defects (Morrisey et al., 1998), and haploinsufficiency results in 60% right and left leaflet fusion (RL)-type BAV (Gharibeh et al., 2018). GATA6 is required for smooth muscle cell (SMC) differentiation and plays a crucial role in the patterning of the aortic arch arteries (Lepore et al., 2006). However, the underlying cellular and molecular mechanisms leading to BAV are still poorly defined.

In this study, we sought to elucidate the cellular and molecular pathways that underlie the function of GATA6 in aortic valve formation. We identified cellular defects that contribute to aberrant aortic valve morphogenesis and show that GATA6 regulates molecular programs directing cell migratory and invasive processes required for endocardial cushion development. Specifically, we demonstrate that GATA6 is necessary for the pro-migratory and invasive functions of C-X-C chemokine receptor type 7 (CXCR7; also known as atypical chemokine receptor 3 or ACKR3), a key regulator of chemokine signaling.

Gata6-null heterozygous mice show BAV and severe aortic insufficiency

We targeted exon 2 of the Gata6 gene, encoding the large GATA-N domain of the protein (Fig. S1A). Indels introduced by CRISPR-Cas9 editing led to a frameshift mutation giving rise to a premature termination codon at amino acid position 291 (V291X), designated henceforth as Gata6STOP/+ (Fig. S1A,B). Histological analysis of E16.5 Gata6STOP/+ mice revealed a right-non-coronary (RN)-type BAV morphology, based on coronary ostia positions relative to the aortic valve leaflets (Fig. S1C), with nearly 70% (20 of 30) penetrance, and a peri-membranous ventricular septal defect (VSD) with 43% (13 of 30) penetrance (Fig. S1D).

To determine the effect of Gata6 haploinsufficiency on cardiac function, we examined 30-week-old adult mice by spectral Doppler echocardiography of the ascending and descending aorta (Fig. 1A). Severe aortic regurgitation was observed in the ascending and especially the descending aorta in 80% (4 of 5) of Gata6STOP/+ mice (Fig. 1A,B), indicating aortic insufficiency due to valve dysfunction. Moreover, ejection fraction and fractional shortening were also reduced in Gata6STOP/+ mice (Fig. 1C), suggesting systolic dysfunction. Other parameters, i.e. left ventricular (LV) diastolic volume, LV systolic volume and LV mass, as well as flow velocity across the mitral valve and strain (percentage of aortic deformation), were comparable between the Gata6STOP/+ mutant and control mice (Fig. S2). Movat's staining (Garvey et al., 1986) in 52-week-old mice aortic valve sections did not show any obvious differences in collagen deposition, elastin fibers or muscle mucin or fibrin composition (Fig. 1D).

BAV associates with aberrant OFT development in Gata6 haploinsufficient mice

We performed isolectin B4 (IsoB4) whole-mount immunostaining of the ventricular endocardium and blood vessel endothelial cells, followed by IMARIS three-dimensional (3D) modeling and volume rendering to determine OFT morphology at E11.5 (Fig. 2A, Movies 1 and 2). E11.5 Gata6STOP/+ mice displayed a shorter OFT compared to that in littermate controls (Fig. 2B). Moreover, we measured tortuosity as a parameter of OFT curvature and found it to be unchanged (Fig. 2B). We also examined OFT caliber as it relates to the separation of the aortic and pulmonary valves by the aorticopulmonary septum (APS), which occurs between E11.5 and E12.5 (Fernández et al., 2009) (Fig. 2C). OFT dimensions in E12.5 Gata6STOP/+ mice were 20% below normal and led to narrowing of the OFT caliber (Fig. 2D). The shortened distance between primitive pulmonary and aortic valves (Fig. 2D) reduced the APS and gave rise to a more circular OFT (Fig. 2D). The length of the major (ventral-dorsal) axis of the OFT was also significantly reduced (Fig. 2D).

To further define the cellular mechanism causing OFT maldevelopment in Gata6STOP/+ mice, we examined cell proliferation in E9.5 hearts, specifically in the S phase of the cell cycle by bromodeoxyuridine (BrdU) immunostaining (Fig. 2E,E′). Proliferation was reduced in Gata6STOP/+ OFT (Fig. 2E,E′), in both the myocardium and endocardium (Fig. 2F). These findings indicate that the shortened and narrowed OFT in Gata6STOP/+ mice can be ascribed to decreased myocardial and endocardial cell proliferation.

Cell-autonomous requirement for GATA6 in the second heart field

We next sought to establish GATA6 requirement in cardiovascular development (Kisanuki et al., 2001). Endothelial-specific Gata6 deletion using the Tie2Cre driver line resulted in overriding aorta with 45% (5 of 11) penetrance (Fig. S3A,B). Gata6 deletion in the myocardium, endocardium (Stanley et al., 2002) and epicardium (Zhou et al., 2008) using the Nkx2.5Cre driver line caused fully penetrant VSD (Fig. S3C,D) and was lethal at birth (Table S1) (van Berlo et al., 2010). In addition, Gata6flox/flox; Nkx2.5Cre mutants showed extensive hypertrabeculation and a thinner compact ventricular wall, consistent with the key roles of GATA6 in the endocardium and myocardium during ventricular chamber development (van Berlo et al., 2010). None of these endocardial- or endothelial-specific models recapitulated the BAV phenotype (Fig. S3B,D).

To determine GATA6 requirement in early cardiac progenitors, we used the Mef2cCre driver line, active in anterior or secondary heart field (SHF)-derived endothelial and myocardial progenitors that give rise to the OFT, right ventricle and ventricular septum (Verzi et al., 2005). Gata6flox/flox; Mef2cCre mutants fully recapitulated the Gata6STOP/+ BAV and VSD phenotypes (Fig. 3A), with 75% (9 of 12) and 58% (7 of 12) penetrance, respectively (Fig. 3B). These data suggest that GATA6 is required cell-autonomously in SHF-derived progenitors.

To gain insight, we examined the SHF marker islet-1 (ISL1) (Quaranta et al., 2018) in Gata6STOP/+ mice. ISL1 marks progenitors for cardiomyocytes, SMCs and endothelial cell lineages, and is gradually switched off as progenitors incorporate into the arterial pole of the heart between E8.5 and E10.5 (Bu et al., 2009). ISL1-positive immunostaining was significantly increased in the distal OFT myocardium, and marginally in the endocardium of E9.5 Gata6STOP/+ hearts (Fig. 3C,D), suggesting that SHF-progenitor differentiation is deficient, whereas no differences were detected in pharyngeal mesoderm (Fig. S3E). These data indicate that defective OFT development in Gata6STOP/+ mice can be ascribed to impaired SHF-derived progenitor differentiation.

Deficient cNCC contribution and SMC differentiation in Gata6STOP/+ OFT

cNCCs constitute a major reservoir of mesenchyme progenitors that are required for patterning of the endocardial cushions in the OFT and aortic arch arteries (Erhardt et al., 2021; Jain et al., 2011). We examined the expression of Sema3c (Fig. 4A), a marker of post-migratory cells and transcriptionally regulated by GATA6 (Brown et al., 2001; Kodo et al., 2009; Lepore et al., 2006). In situ hybridization on E12.5 Gata6STOP/+ OFT frontal sections revealed markedly reduced Sema3c staining, suggesting reduced cNCC presence (Fig. 4A). Moreover, GATA6 is required for neural crest cell-derived SMC differentiation (Lepore et al., 2006; Losa et al., 2017). Analysis of smooth muscle differentiation along distal, medial and proximal segments of E12.5 Gata6STOP/+ mice OFT (Fig. 4B) revealed deficient smooth muscle actin (SMA, encoded by Acta2)-positive immunostaining in the APS, consistent with reduced cNCC presence and/or impaired SMC differentiation (Fig. 4C).

Gata6 regulates a pro-migratory gene signature in OFT development

For further insight, we performed RNA sequencing (RNA-seq) on E11.5 Gata6STOP/+ and control OFTs (Fig. 5). EdgeR identified 113 differentially expressed genes, of which 53 were upregulated and 60 downregulated (Fig. 5A; Table S2). As expected, Gata6 expression was decreased by RNA-seq [log2(fold change or FC)=−0.42331; adjusted P-value=0.022647]. Quantitative real-time PCR confirmed that Gata6 was expressed at 60% of the wild-type level (Fig. S3F). Ingenuity Pathway Analysis uncovered only a few enriched categories, i.e. ‘congenital anomaly of cardiovascular system’ (Hand1, Foxf1) and ‘differentiation of bone cells’ (Postn) (Fig. 5B,D). To expand the gene categories, we performed gene set enrichment analysis (GSEA) against ‘Hallmark’ sets (Liberzon et al., 2015) (Fig. 5C). Enriched pathways were ‘G2M checkpoint’ and ‘KRAS signaling UP’, potentially highlighting a proliferative defect.

In contrast, Gata6STOP/+ OFT gene expression was depleted for functional categories such as ‘cell movement of neurons’ (Cxcr7, Erbb3), ‘chemotaxis’ (Sema3a, Robo2, Ntn1, Nrp2, Ntrk2), ‘invasion of cells’ (Diras2, Itga4) and ‘proliferation of endothelial cells’ (Anxa8) (Fig. 5B,D). Related categories such as ‘angiogenesis’ (Pdgfra, Sfrp2) and ‘development of epithelial tissue’ (Col12a1, Lamc2) were also depleted (Fig. 5B,D). Consistent with this, GSEA revealed depleted gene sets such as ‘epithelial-mesenchymal transition’ (Fig. 5C) (von Gise and Pu, 2012) and those related to cellular stress pathways such as ‘DNA repair’ and ‘P53 pathway’, which is potentially linked to ‘oxidative phosphorylation’ (Ndor1), ‘glycolysis’ (Gpc1) and lipid metabolism dysfunction [i.e. ‘cholesterol homeostasis’ (Apoe) and ‘adipogenesis’] (Fig. 5C,D). Moreover, depleted ‘myogenesis’ (Mstn, Mybpc1) is consistent with GATA6 requirement for myocardial differentiation (Sharma et al., 2020).

Altogether, the bioinformatic analysis indicates that GATA6 regulates a set of processes associated with cellular proliferation and especially migration, potentially involved in development and patterning of the OFT.

CXCR7 mediates GATA6-dependent migration and invasion

Endocardial cushion explant assays can recapitulate morphogenetic processes occurring in vivo, including the epithelial-mesenchymal transition and proliferation (Runyan and Markwald, 1983; Timmerman et al., 2004). Therefore, to corroborate the potential link between OFT underdevelopment and defective cell migration in Gata6STOP/+ mutants, we performed explant assays (Fig. S4A). SMA+ mesenchymal cells migrating out from OFT explants are derived from both the endocardial-mesenchymal transition and neural crest cells that have colonized the endocardial cushions (MacGrogan et al., 2016; Ridge et al., 2021). Consistent with the assessment in the RNA-seq analysis, mesenchymal migratory and invasive behaviors in E11.5 Gata6STOP/+ explants were below normal, whereas proliferation was unchanged (Fig. S4B,C).

We then focused on Cxcr7 as a potential mediator of Gata6 function in the OFT. Cxcr7 gene expression was found to be decreased in Gata6STOP/+ OFT (FC=−1.27, adjusted P-value=0.0196) (Fig. 5A; Table S2) and, importantly, Cxcr7 mutant mice recapitulate valvulo-septal phenotypes found in the Gata6STOP/+ mutants (Sierro et al., 2007). CXCR7, together with the chemokine receptor CXCR4 and its ligand CXCL12 (also known as SDF1), regulates cellular migration and proliferation in many developmental and pathological settings (Döring et al., 2014).

We evaluated CXCR7 requirement for GATA6-mediated cellular motility. We tested the effect of VUF11207, a specific CXCR7 agonist that induces recruitment of the cytosolic scaffolding protein β-arrestin, followed by its internalization and downstream pathway activation (Wijtmans et al., 2012). Thus, whereas supplementing wild-type OFT explants with 100 nM VUF11207 had no effect (Fig. S5A,A′,B,B′,D), supplementing with 1 µM VUF11207 led to a significant dose-dependent decline in migration, defined as the furthest distance migrated (Fig. S5A,A′,C,Cʹ,D). In contrast, a dose-dependent increase was found with 1 µM of VUF11207 for the invasion parameter, defined as the furthest distance invaded, relative to that seen in the control (Fig. S5A″,C″,D), whereas supplementing with 100 nM VUF11207 had no effect (Fig. S5A″,B″,D). However, providing the explants with either 100 nM or 1 µM of the agonist did not result in any proliferative changes (Fig. S5D).

Next, consistent with our pilot experiment, supplementing control explants with 1 µM VUF11207 resulted in decreased migration (Fig. 6A,A′,C) and increased invasion (Fig. 6A,A″,C), confirming that migratory processes taking place in endocardial cushion development are regulated by CXCR7. However, Gata6STOP/+ explants did not respond to 1 µM VUF11207 supplementation (Fig. 6B,B″,C), suggesting that migration and invasion mediated by CXCR7 are dependent on Gata6. No effects on proliferation were observed in either Gata6STOP/+ or littermate control explants (Fig. 6C).

CXCR7 modulates CXCR4 function by sequestering CXCL12 (Boldajipour et al., 2008; Chen et al., 2022; Naumann et al., 2010). To determine whether CXCR7 modulation of CXCR4 is dependent on Gata6, we tested the effects of AMD3100, a dual-function chemical compound known to simultaneously antagonize CXCR4 (Debnath et al., 2013) and allosterically activate CXCR7 (Kalatskaya et al., 2009). Having previously defined an optimal dose-dependent AMD3100 concentration (Ferreira Oliveira, 2020), we found that supplementing wild-type OFT explants with 1 µM AMD3100 did not affect migration (Fig. 6D,D′,F) but significantly increased invasion (Fig. 6D,D″,F) and proliferation (Fig. 6D,D′,F). This suggests that CXCR4 regulates mesenchymal cell proliferation in OFT explants, whereas CXCR7 likely regulates migration and invasion. Importantly, supplementing Gata6STOP/+ OFT explants with 1 µM AMD3100 did not alter mesenchymal cell migration, invasion or proliferation relative to that in Gata6STOP/+ OFT explants in the absence of AMD3100 (Fig. 6E,E″,F). Thus, Gata6 is required for the pro-invasive effects mediated by simultaneous CXCR4 inhibition and CXCR7 activation via AMD3100, but not for pro-proliferative effects.

In this study, we show that Gata6 haploinsufficiency results in highly penetrant BAV and characterize the phenotype to OFT alignment anomalies, including a membranous VSD, a common component of congenital heart disease in humans. We elucidate the deficient cellular differentiation and proliferation mechanisms associated with congenital heart defects in Gata6STOP/+ mice. We find a requirement for functional hierarchy between GATA6 and CXCR7 in mesenchymal cell migration and invasion in endocardial cushion development.

The failure to elongate the OFT likely underlies the ventricular-arterial alignment defects observed in Gata6STOP/+ embryos. The OFT length in Gata6STOP/+ embryos may be insufficient for normal ventricular septal alignment and atrioventricular canal endocardial cushion fusion (Gittenberger-de Groot et al., 2014). Normally, SHF progenitors are added progressively to the arterial pole, which subsequently becomes incorporated into the ventricular outlets at the base of the great arteries (Rochais et al., 2009). Deficient SHF deployment is the primary cause of septal defects at the ventricular-arterial connection with the OFT (Gittenberger-de Groot et al., 2014). The delayed differentiation of ISL1-positive progenitors and decreased proliferation of distal OFT endocardial/myocardial lineages observed in Gata6STOP/+ mice are consistent with dysregulation in SHF-derived progenitors.

A cell-autonomous GATA6 requirement in the SHF is further supported by conditional Gata6 inactivation, which reproduces the Gata6STOP/+ mutant phenotypes. Moreover, interactions between the SHF and cNCCs are crucial for normal OFT development, and perturbed SHF signaling indirectly impacts cNCC distribution in the endocardial cushions (Ward et al., 2005), causing mispatterning of the aortic valve leaflets (Waldo et al., 2005). GATA6 is also required cell-autonomously in post-migratory cNCCs for Sema3c expression, smooth muscle differentiation and aortic arch artery patterning (Lepore et al., 2005). Medial arterial defects have been previously linked to BAV aortopathy (Yassine et al., 2017), and aortic structural integrity is likely compromised in the adult Gata6STOP/+ mice with BAV, based on our spectral Doppler echocardiography experiments revealing retrograde flow (ascending and descending) and severe aortic regurgitation.

Cell movement and invasion are necessary for organogenesis (Aman and Piotrowski, 2010). GATA factors control processes that repress epithelial characteristics and confer migratory capabilities to epithelial cells involved in mesoderm development, neural crest formation and wound healing (Campbell et al., 2011). Consistent with this, our transcriptional profiling suggests that pro-migratory and invasive processes are globally impaired in Gata6STOP/+ OFT development. These pathways principally involve canonical cell guidance molecules and receptors, including semaphorins, ephrins and netrins, and oncogenic pathways such as ERBB3 and NTRK. Mis-regulation of these pathways has been associated with developmental patterning defects and congenital heart disease (Claro and Ferro, 2020; Erickson et al., 1997; Hang et al., 2021; Mommersteeg et al., 2015).

We focused on Cxcr7 as potentially acting downstream of GATA6. CXCR7 is a decoy receptor for CXCL12–CXCR4 signaling with broad effects on chemotaxis, cell proliferation and migration (Chédotal et al., 2005; Shi et al., 2020). Moreover, this chemokine pathway plays key role(s) in aortic valve development (Ivins et al., 2015; Ridge et al., 2021; Sierro et al., 2007; Yu et al., 2011) and its disruption recapitulates Gata6STOP/+ phenotypes (Duval et al., 2022; Ivins et al., 2015; Sierro et al., 2007; Yu et al., 2011). As expected, we found that agonist-mediated activation of CXCR7 inhibited migration and increased invasion, likely by blocking the pro-migratory effects of CXCL12–CXCR4 signaling (Wijtmans et al., 2012). We also expected that blocking CXCL12 binding to CXCR4 might impair migration (Ridge et al., 2021; Yu et al., 2011) but, given that the antagonist simultaneously activated CXCR7 (Kalatskaya et al., 2009), only effects on invasion and proliferation were found. This anti-invasive and anti-proliferative response is in agreement with previous findings indicating that CXCL12–CXCR4 signaling promotes migratory processes and inhibits proliferation in the valve mesenchyme (Ridge et al., 2021). Regardless, even in the presence of highly potent CXCR7 activity modulators, no effects were found in Gata6STOP/+ mutants, suggesting that GATA6 is required for CXCR7-mediated migratory or invasive processes in endocardial cushion formation. Our findings are summarized in Fig. S5E,F. Proliferation was not affected in the Gata6STOP/+ explants or by the modulators of CXCR7 activity, suggesting that proliferation and migration/invasion may be independently regulated in this context. We have not found any bioinformatic evidence of GATA-binding sequences in proximity to the Cxcr7 gene supportive of direct transcriptional regulation, suggesting that GATA6 control of Cxcr7 could be indirect. Further research is required to understand the mechanisms by which GATA6 regulates CXCR7 and CXCL12–CXCR4 signaling during endocardial cushion formation and in congenital heart disease.

Animal studies and generation of the Gata6STOP/+ mouse line

Animal experiments adhered to Animal Research: Reporting of In Vivo Experiments (ARRIVE) guidelines and complied with ethical standards outlined in the UK Animal (Scientific Procedures) Act 1986, Royal Decree 53/2013, EU Directive 2010/63/EU and Recommendation 2007/526/EC on the protection of animals used for experimental and other scientific purposes, enacted in Spanish law under Real Decreto 1201/2005. Approval was obtained from the CNIC Animal Experimentation Ethics Committee, Complutense University of Madrid, and Community of Madrid (reference PROEX 155.7/20).

The standard knockout Gata6STOP mouse line was generated using CRISPR-Cas9 technology. We designed a specific CRISPR RNA (crRNA) using the Spanish National Biotechnology Centre (CNB)-Spanish National Research Council (CSIC) web tool Breaking-Cas (https://bioinfogp.cnb.csic.es/tools/breakingcas/) and CRISPOR-TEFOR online tool (http://crispor.gi.ucsc.edu/crispor.py). The sequences with higher specificity, fewer off-targets and higher efficiency scores were selected (Table S3). crRNA sequences were used as annealed two-part synthetic crRNA (Alt-RR CRISPR-Cas9 crRNA, 2 nmol, Integrated DNA Technologies or IDT) and trans-activating crRNA (tracrRNA; Alt-RR CRISPR-Cas9 tracrRNA, 5 nmol, IDT, 1072532) molecules, known as the Alt-R CRISPR-Cas9 system. We used the Alt-RR Streptococcus pyogenes Cas9 nuclease (IDT, 1081058). Single-stranded DNA templates were custom synthetic genes (Megamer single-stranded gene fragments, IDT). The complete guide RNA complex (annealed crRNA and tracrRNA) was diluted in microinjection buffer (1 mM Tris HCl, pH 7.5; 0.1 mM EDTA) and incubated with Cas9 protein to obtain active ribonucleoprotein complexes. The final concentrations of components in ribonucleoprotein preparations were 0.15-0.61 pmol/μl of crRNA and tracrRNA, 20-30 ng/μl of Cas9 protein and 10 ng/μl of single-stranded oligodeoxynucleotides (Table S3). The final injection mixes were passed through centrifugal filter units (UFC30VV25, EMD Millipore) and spun at 21,000 g for 5 min at room temperature. Reagents were microinjected into one-cell fertilized C57BL/6 mouse embryos (Harms et al., 2014). Pups were screened for the targeted mutations using external primers by PCR (Table S3). Sanger sequencing was performed to confirm target gene editing. The selected founder was backcrossed to the C57BL/6 background.

Other mouse strains used in this study are as follows: Gata6flox (Sodhi et al., 2006), Tie2Cre (Kisanuki et al., 2001), Nkx2.5Cre (Stanley et al., 2002) and Mef2cCre (Verzi et al., 2005).

Tissue processing for histological procedures

For histological procedures, whole embryos or torsos were fixed in 4% paraformaldehyde (PFA, Electron Microscopy Sciences, 50980487) overnight at 4°C. Paraffin-embedded embryos/torsos were cut into 7 μm sections using a HistoCore AUTOCUT rotary microtome (Leica Biosystems). Hematoxylin and Eosin staining was performed according to standard protocols.

RNA probe synthesis and in situ hybridization

Antisense RNA probes were designed to be exon spanning or complementary to the 3′ untranslated region sequences. Targeted regions were amplified by PCR using 100 ng of cDNA from hearts. The PCR product was ligated to the pGEM-T easy vector (Promega, A1360) and transformed into DH5α competent Escherichia coli cells. Plasmid DNA was isolated from positive colonies, linearized and transcribed by T7 or Sp6 polymerase (RPOLT7-RO and RPOLSP6-RO, Roche). The RNA probes were purified and stored at −80°C. Primer sequences used for probe generation are listed in Table S3. In situ hybridization was performed as described previously (Kanzler et al., 1998).

Immunohistochemistry

Paraffin-embedded 7 μm sections were citrate-unmasked and incubated overnight with the primary antibodies. For antigens that did not require amplification, fluorescent dye-conjugated secondary antibodies were directly incubated for 1 h. For signal amplification, a secondary biotin-conjugated antibody was used, and signal detection was performed using the Tyramide Signal Amplification Plus Fluorescein kit (PerkinElmer, NEL741B001KT). All antibodies used are listed in Table S3.

Whole-mount immunofluorescence

E11.5 embryos were fixed in 4% PFA overnight at 4°C and then kept in PBS at 4°C. Tissue clearing was performed with CUBIC-1 (Susaki et al., 2014) at 37°C for 1.5 days. Embryos were incubated with primary antibody endomucin (V.7C7) (Santa Cruz Biotechnology, sc-65495) at 4°C for 3 days. Incubation with secondary antibody was performed at 4°C for 3 days. Then, incubation with CUBIC-2 (Susaki et al., 2014) was performed at room temperature overnight. Quantification of OFT length and tortuosity and 3D modeling were performed with IMARIS software. Six embryos of each genotype were studied. Statistical analyses were performed using unpaired two-tailed Student's t-test, and differences were considered statistically significant at P<0.05. Data are presented as mean±s.d.

OFT explant assay

Collagen gels were prepared by pouring a solution of 1.5 mg/ml rat-tail collagen type I (Corning, 354236), filtered 1 N NaOH, 10× Dulbecco's modified Eagle medium (DMEM; 11430030, Gibco, Thermo Fisher Scientific) and Biopak water into four-well microculture dishes (500 μl per well). The gels were solidified inside a 37°C and 5% CO2 incubator for 20 min before being washed three times for 30 min with explant medium in a solution containing 1× DMEM, GlutaMAX-I CTS (A12860-01, Gibco, Thermo Fisher Scientific), 10% fetal bovine serum (F7524-1654682, Sigma-Aldrich) and a 1× solution of penicillin (10 U/ml), streptomycin (10 mg/ml; P4333, Sigma-Aldrich) and fungizone (250 µg/ml, SV30078.01, Hyclone). The gel was pre-soaked using explant medium containing the drug of interest. Drugs are listed in Table S3. The following day, OFTs were dissected from E11.5 embryos in fresh 1× DMEM, opened, placed with the endocardium side lying face down on gel and allowed to attach for 24 h at 37°C and 5% CO2. The next day, we added 500 μl of explant medium, containing the required drug and 1% insulin-transferrin-selenium (ITS; 51500056, Gibco, Thermo Fisher Scientific). After 3 days in culture, explants were fixed in 4% PFA and immunofluorescence was performed. Statistical analyses were performed using unpaired two-tailed Student's t-test. Data are presented as mean±s.d.

Microscopy and confocal imaging

Brightfield imaging was performed using an Olympus BX51 Microscope and Olympus cellSense software. Confocal images were acquired using a Nikon A1-R confocal microscope and Zeiss LSM700 confocal microscope. Confocal data were processed using the ZEN 2012 software (black edition) and image analysis was performed using the Fiji image software. Images were processed in Adobe Photoshop Creative Suite 5.1.

Echocardiographic recordings

LV function and wall thickness were analyzed by transthoracic echocardiography in 30-week-old male mice. Mice were mildly anesthetized by inhalation of isoflurane/oxygen (1-2%/98.75%), adjusted to obtain a target heart rate of 450±50 beats per minute, and examined with a 30 MHz transthoracic echocardiography probe. Images were obtained with a Vevo 2100 imaging system (VisualSonics). LV long-axis M-mode views were obtained as described previously (Cruz-Adalia et al., 2010). From these images, LV mass, LV systolic volume and LV diastolic volume were measured, LV systolic function was assessed by estimating the LV shortening fraction and LV ejection fraction, and the percentage of strain was assessed by measuring the percentage of ascendant aorta deformation by estimating the diameter in diastole and systole of the aortic arch (Cruz-Adalia et al., 2010). Statistical analyses were performed using unpaired two-tailed Student's t-test. Data are presented as mean±s.d.

RNA-seq

RNA was isolated from E11.5 Gata6STOP/+ and control OFTs. Samples were distributed in three pools of four pairs of OFT per genotype. Tissue was homogenized with a pestle mechanical homogenizer and RNA was extracted with an Arcturus PicoPure RNA Isolation kit (Thermo Fisher Scientific, KIT0214). RNA libraries were prepared using the NEBNext Ultra II Directional RNA Library Prep Kit (New England Biolabs) and sequenced in a Nextseq 2000 Illumina sequencer using a 60 bp single-end elongation protocol. Sequenced reads were checked for quality control and pre-processed using Cutadapt v1.18 (Kechin et al., 2017) to remove adapter contaminants. Resulting reads were aligned and gene expression quantified using RSEM v1.2.3 (Li and Dewey, 2011) over mouse reference GRCm38 with Ensembl genebuild. Differential gene expression was analyzed with the EdgeR R package (v3.32.1 on R 4.0.3) (Robinson et al., 2010). Genes with one count per million in at least three samples were defined as expressed and retained for later analysis. Counts were normalized by the trimmed mean of M-values (TMM) method. Differential gene expression was tested using a generalized linear model as implemented in the EdgeR package. Genes showing altered expression with an adjusted P-value <0.05 were considered differentially expressed. The set of differentially expressed genes was used for functional analysis with Ingenuity Pathway Analysis Software (Qiagen-IPA) (Krämer et al., 2014), where we used Benjamini–Hochberg adjusted P-value <0.05 for significance. GSEA was performed on the complete set of expressed genes, against the Hallmark term database. A nominal P-value <0.1 was used to select for significantly enriched gene sets.

Statistical analysis

Sample sizes, statistical tests and P-values are specified in the corresponding figure legends and corresponding subsections of the Materials and Methods. For comparisons between two groups, mean±s.d. is represented and unpaired two-tailed Student's t-test was performed. For experiments comparing two groups of categorical variables, mean±s.d. per group is represented and a Fisher's exact test was performed. Differences were considered statistically significant at P<0.05 (Table S4). Statistical analysis and graphical representation were performed using GraphPad Prism v8.

The authors thank A. Galicia and L. Méndez for mouse husbandry; Dra. Anabel Rojas from Universidad Pablo de Olavide, Seville, Spain, for kindly providing the Gata6flox line; the Centro Nacional de Investigaciones Cardiovasculares Carlos III (CNIC) Genomics, Transgenesis, Comparative Medicine and Microscopy (especially V. Labrador) Units for technical support; and M. J. Gómez (Bioinformatics Unit, CNIC) for expert advice with the RNA-seq analysis. The CNIC is supported by the Instituto de Salud Carlos III, the Spanish Ministry of Science and Innovation and the Pro-CNIC Foundation, and is a Severo Ochoa Center of Excellence (grant CEX2020001041-S) financed by MCIN/AEI/10.13039/501100011033.

Author Contributions

Conceptualization: D.M., J.L.d.l.P.; Methodology: R.P.-S.; Validation: D.M., J.L.d.l.P.; Formal analysis: R.P.-S.; Investigation: R.P.-S.; Resources: J.L.d.l.P.; Data curation: R.P.-S.; Writing - original draft: R.P.-S., D.M.; Writing - review & editing: D.M., J.L.d.l.P.; Supervision: D.M., J.L.d.l.P.; Funding acquisition: J.L.d.l.P.

Funding

This study was supported by grants PID2022-104776RB-100 and CB16/11/00399 (CIBER CV) from Ministerio de Ciencia e Innovación and Agencia Estatal de Investigación (MCIN/AEI/10.13039/501100011033), and the 'la Caixa' Foundation (reference HR23-00084) to J.L.d.l.P. Support for this publication also came from the European Regional Development Fund. Open access funding provided by Centro Nacional de Investigaciones Cardiovasculares. Deposited in PMC for immediate release.

Data availability

Data are deposited in the NCBI Gene Expression Omnibus database under the accession number GSE255237.

Alonso-Montes
,
C.
,
Martín
,
M.
,
Martínez-Arias
,
L.
,
Coto
,
E.
,
Naves-Díaz
,
M.
,
Morís
,
C.
,
Cannata-Andía
,
J. B.
and
Rodríguez
,
I.
(
2018
).
Variants in cardiac GATA genes associated with bicuspid aortic valve
.
Eur. J. Clin. Invest.
48
,
e13027
.
Aman
,
A.
and
Piotrowski
,
T.
(
2010
).
Cell migration during morphogenesis
.
Dev. Biol.
341
,
20
-
33
.
Boldajipour
,
B.
,
Mahabaleshwar
,
H.
,
Kardash
,
E.
,
Reichman-Fried
,
M.
,
Blaser
,
H.
,
Minina
,
S.
,
Wilson
,
D.
,
Xu
,
Q.
and
Raz
,
E.
(
2008
).
Control of chemokine-guided cell migration by ligand sequestration
.
Cell
132
,
463
-
473
.
Brown
,
C. B.
,
Feiner
,
L.
,
Lu
,
M. M.
,
Li
,
J.
,
Ma
,
X.
,
Webber
,
A. L.
,
Jia
,
L.
,
Raper
,
J. A.
and
Epstein
,
J. A.
(
2001
).
PlexinA2 and semaphorin signaling during cardiac neural crest development
.
Development
128
,
3071
-
3080
.
Bu
,
L.
,
Jiang
,
X.
,
Martin-Puig
,
S.
,
Caron
,
L.
,
Zhu
,
S.
,
Shao
,
Y.
,
Roberts
,
D. J.
,
Huang
,
P. L.
,
Domian
,
I. J.
and
Chien
,
K. R.
(
2009
).
Human ISL1 heart progenitors generate diverse multipotent cardiovascular cell lineages
.
Nature
460
,
113
-
117
.
Campbell
,
K.
,
Whissell
,
G.
,
Franch-Marro
,
X.
,
Batlle
,
E.
and
Casanova
,
J.
(
2011
).
Specific GATA factors act as conserved inducers of an endodermal-EMT
.
Dev. Cell
21
,
1051
-
1061
.
Chédotal
,
A.
,
Kerjan
,
G.
and
Moreau-Fauvarque
,
C.
(
2005
).
The brain within the tumor: new roles for axon guidance molecules in cancers
.
Cell Death Differ.
12
,
1044
-
1056
.
Chen
,
D.
,
Xia
,
Y.
,
Zuo
,
K.
,
Wang
,
Y.
,
Zhang
,
S.
,
Kuang
,
D.
,
Duan
,
Y.
,
Zhao
,
X.
and
Wang
,
G.
(
2022
).
Author correction: crosstalk between SDF-1/CXCR4 and SDF-1/CXCR7 in cardiac stem cell migration
.
Sci. Rep.
12
,
2709
.
Claro
,
V.
and
Ferro
,
A.
(
2020
).
Netrin-1: focus on its role in cardiovascular physiology and atherosclerosis
.
JRSM Cardiovasc. Dis.
9
,
2048004020959574
.
Conway
,
S. J.
,
Kruzynska-Frejtag
,
A.
,
Kneer
,
P. L.
,
Machnicki
,
M.
and
Koushik
,
S. V.
(
2003
).
What cardiovascular defect does my prenatal mouse mutant have, and why?
Genesis
35
,
1
-
21
.
Cruz-Adalia
,
A.
,
Jiménez-Borreguero
,
L. J.
,
Ramírez-Huesca
,
M.
,
Chico-Calero
,
I.
,
Barreiro
,
O.
,
López-Conesa
,
E.
,
Fresno
,
M.
,
Sánchez-Madrid
,
F.
and
Martín
,
P.
(
2010
).
CD69 limits the severity of cardiomyopathy after autoimmune myocarditis
.
Circulation
122
,
1396
-
1404
.
Debnath
,
B.
,
Xu
,
S.
,
Grande
,
F.
,
Garofalo
,
A.
and
Neamati
,
N.
(
2013
).
Small molecule inhibitors of CXCR4
.
Theranostics
3
,
47
-
75
.
Della Corte
,
A.
,
Bancone
,
C.
,
Quarto
,
C.
,
Dialetto
,
G.
,
Covino
,
F. E.
,
Scardone
,
M.
,
Caianiello
,
G.
and
Cotrufo
,
M.
(
2007
).
Predictors of ascending aortic dilatation with bicuspid aortic valve: a wide spectrum of disease expression
.
Eur. J. Cardiothorac. Surg.
31
,
397
-
404; discussion 404-5
.
Döring
,
Y.
,
Pawig
,
L.
,
Weber
,
C.
and
Noels
,
H.
(
2014
).
The CXCL12/CXCR4 chemokine ligand/receptor axis in cardiovascular disease
.
Front. Physiol.
5
,
212
.
Duval
,
V.
,
Alayrac
,
P.
,
Silvestre
,
J. S.
and
Levoye
,
A.
(
2022
).
Emerging roles of the atypical chemokine receptor 3 (ACKR3) in cardiovascular diseases
.
Front. Endocrinol. (Lausanne)
13
,
906586
.
Erhardt
,
S.
,
Zheng
,
M.
,
Zhao
,
X.
,
Le
,
T. P.
,
Findley
,
T. O.
and
Wang
,
J.
(
2021
).
The cardiac neural crest cells in heart development and congenital heart defects
.
J. Cardiovasc. Dev. Dis.
8
,
89
.
Erickson
,
S. L.
,
O'shea
,
K. S.
,
Ghaboosi
,
N.
,
Loverro
,
L.
,
Frantz
,
G.
,
Bauer
,
M.
,
Lu
,
L. H.
and
Moore
,
M. W.
(
1997
).
ErbB3 is required for normal cerebellar and cardiac development: a comparison with ErbB2-and heregulin-deficient mice
.
Development
124
,
4999
-
5011
.
Fernández
,
B.
,
Durán
,
A. C.
,
Fernández-Gallego
,
T.
,
Fernández
,
M. C.
,
Such
,
M.
,
Arqué
,
J. M.
and
Sans-Coma
,
V.
(
2009
).
Bicuspid aortic valves with different spatial orientations of the leaflets are distinct etiological entities
.
J. Am. Coll. Cardiol.
54
,
2312
-
2318
.
Ferreira Oliveira
,
V. L.
(
2020
).
Unraveling the molecular mechanisms behind the regulation of mesenchymal cell proliferation during cardiac valve remodeling. PhD thesis, Centro Nacional de Investigaciones Cardiovasculares Carlos III, Madrid
. http://hdl.handle.net/10486/694423
Garvey
,
W.
,
Fathi
,
A.
,
Bigelow
,
F.
,
Carpenter
,
B.
and
Jimenez
,
C.
(
1986
).
Improved Movat pentachrome stain
.
Stain Technol.
61
,
60
-
62
.
Gharibeh
,
L.
,
Komati
,
H.
,
Bossé
,
Y.
,
Boodhwani
,
M.
,
Heydarpour
,
M.
,
Fortier
,
M.
,
Hassanzadeh
,
R.
,
Ngu
,
J.
,
Mathieu
,
P.
,
Body
,
S.
et al. 
(
2018
).
GATA6 regulates aortic valve remodeling, and its haploinsufficiency leads to right-left type bicuspid aortic valve
.
Circulation
138
,
1025
-
1038
.
Gittenberger-De Groot
,
A. C.
,
Calkoen
,
E. E.
,
Poelmann
,
R. E.
,
Bartelings
,
M. M.
and
Jongbloed
,
M. R.
(
2014
).
Morphogenesis and molecular considerations on congenital cardiac septal defects
.
Ann. Med.
46
,
640
-
652
.
Giusti
,
B.
,
Sticchi
,
E.
,
De Cario
,
R.
,
Magi
,
A.
,
Nistri
,
S.
and
Pepe
,
G.
(
2017
).
Genetic bases of bicuspid aortic valve: the contribution of traditional and high-throughput sequencing approaches on research and diagnosis
.
Front. Physiol.
8
,
612
.
Hang
,
P. Z.
,
Zhu
,
H.
,
Li
,
P. F.
,
Liu
,
J.
,
Ge
,
F. Q.
,
Zhao
,
J.
and
Du
,
Z. M.
(
2021
).
The emerging role of BDNF/TrkB signaling in cardiovascular diseases
.
Life (Basel)
11
,
70
.
Harms
,
D. W.
,
Quadros
,
R. M.
,
Seruggia
,
D.
,
Ohtsuka
,
M.
,
Takahashi
,
G.
,
Montoliu
,
L.
and
Gurumurthy
,
C. B.
(
2014
).
Mouse genome editing using the CRISPR/Cas system
.
Curr. Protoc. Hum. Genet.
83
,
15.7.1
-
15.727
.
Henderson
,
D. J.
,
Eley
,
L.
,
Turner
,
J. E.
and
Chaudhry
,
B.
(
2022
).
Development of the human arterial valves: understanding bicuspid aortic valve
.
Front. Cardiovasc. Med.
8
,
802930
.
Hoffman
,
J. I.
and
Kaplan
,
S.
(
2002
).
The incidence of congenital heart disease
.
J. Am. Coll. Cardiol.
39
,
1890
-
1900
.
Ivins
,
S.
,
Chappell
,
J.
,
Vernay
,
B.
,
Suntharalingham
,
J.
,
Martineau
,
A.
,
Mohun
,
T. J.
and
Scambler
,
P. J.
(
2015
).
The CXCL12/CXCR4 axis plays a critical role in coronary artery development
.
Dev. Cell
33
,
455
-
468
.
Jain
,
R.
,
Engleka
,
K. A.
,
Rentschler
,
S. L.
,
Manderfield
,
L. J.
,
Li
,
L.
,
Yuan
,
L.
and
Epstein
,
J. A.
(
2011
).
Cardiac neural crest orchestrates remodeling and functional maturation of mouse semilunar valves
.
J. Clin. Invest.
121
,
422
-
430
.
Kalatskaya
,
I.
,
Berchiche
,
Y. A.
,
Gravel
,
S.
,
Limberg
,
B. J.
,
Rosenbaum
,
J. S.
and
Heveker
,
N.
(
2009
).
AMD3100 is a CXCR7 ligand with allosteric agonist properties
.
Mol. Pharmacol.
75
,
1240
-
1247
.
Kanzler
,
B.
,
Kuschert
,
S. J.
,
Liu
,
Y. H.
and
Mallo
,
M.
(
1998
).
Hoxa-2 restricts the chondrogenic domain and inhibits bone formation during development of the branchial area
.
Development
125
,
2587
-
2597
.
Kechin
,
A.
,
Boyarskikh
,
U.
,
Kel
,
A.
and
Filipenko
,
M.
(
2017
).
cutPrimers: a new tool for accurate cutting of primers from reads of targeted next generation sequencing
.
J. Comput. Biol.
24
,
1138
-
1143
.
Kisanuki
,
Y. Y.
,
Hammer
,
R. E.
,
Miyazaki
,
J.
,
Williams
,
S. C.
,
Richardson
,
J. A.
and
Yanagisawa
,
M.
(
2001
).
Tie2-Cre transgenic mice: a new model for endothelial cell-lineage analysis in vivo
.
Dev. Biol.
230
,
230
-
242
.
Kodo
,
K.
,
Nishizawa
,
T.
,
Furutani
,
M.
,
Arai
,
S.
,
Yamamura
,
E.
,
Joo
,
K.
,
Takahashi
,
T.
,
Matsuoka
,
R.
and
Yamagishi
,
H.
(
2009
).
GATA6 mutations cause human cardiac outflow tract defects by disrupting semaphorin-plexin signaling
.
Proc. Natl. Acad. Sci. USA
106
,
13933
-
13938
.
Krämer
,
A.
,
Green
,
J.
,
Pollard
,
J.
, Jr.
and
Tugendreich
,
S.
(
2014
).
Causal analysis approaches in ingenuity pathway analysis
.
Bioinformatics
30
,
523
-
530
.
Ladich
,
E.
,
Nakano
,
M.
,
Carter-Monroe
,
N.
and
Virmani
,
R.
(
2011
).
Pathology of calcific aortic stenosis
.
Future Cardiol.
7
,
629
-
642
.
Lepore
,
J. J.
,
Cappola
,
T. P.
,
Mericko
,
P. A.
,
Morrisey
,
E. E.
and
Parmacek
,
M. S.
(
2005
).
GATA-6 regulates genes promoting synthetic functions in vascular smooth muscle cells
.
Arterioscler. Thromb. Vasc. Biol.
25
,
309
-
314
.
Lepore
,
J. J.
,
Mericko
,
P. A.
,
Cheng
,
L.
,
Lu
,
M. M.
,
Morrisey
,
E. E.
and
Parmacek
,
M. S.
(
2006
).
GATA-6 regulates semaphorin 3C and is required in cardiac neural crest for cardiovascular morphogenesis
.
J. Clin. Invest.
116
,
929
-
939
.
Li
,
B.
and
Dewey
,
C. N.
(
2011
).
RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome
.
BMC Bioinformatics
12
,
323
.
Liberzon
,
A.
,
Birger
,
C.
,
Thorvaldsdóttir
,
H.
,
Ghandi
,
M.
,
Mesirov
,
J. P.
and
Tamayo
,
P.
(
2015
).
The molecular signatures database (MSigDB) hallmark gene set collection
.
Cell Syst.
1
,
417
-
425
.
Lin
,
C. J.
,
Lin
,
C. Y.
,
Chen
,
C. H.
,
Zhou
,
B.
and
Chang
,
C. P.
(
2012
).
Partitioning the heart: mechanisms of cardiac septation and valve development
.
Development
139
,
3277
-
3299
.
Losa
,
M.
,
Latorre
,
V.
,
Andrabi
,
M.
,
Ladam
,
F.
,
Sagerström
,
C.
,
Novoa
,
A.
,
Zarrineh
,
P.
,
Bridoux
,
L.
,
Hanley
,
N. A.
,
Mallo
,
M.
et al. 
(
2017
).
A tissue-specific, Gata6-driven transcriptional program instructs remodeling of the mature arterial tree
.
Elife
6
,
e31362
.
Ma
,
M.
,
Li
,
Z.
,
Mohamed
,
M. A.
,
Liu
,
L.
and
Wei
,
X.
(
2021
).
Aortic root aortopathy in bicuspid aortic valve associated with high genetic risk
.
BMC Cardiovasc. Disord.
21
,
413
.
MacGrogan
,
D.
,
D'amato
,
G.
,
Travisano
,
S.
,
Martinez-Poveda
,
B.
,
Luxán
,
G.
,
Del Monte-Nieto
,
G.
,
Papoutsi
,
T.
,
Sbroggio
,
M.
,
Bou
,
V.
,
Gomez-Del Arco
,
P.
et al. 
(
2016
).
Sequential ligand-dependent notch signaling activation regulates valve primordium formation and morphogenesis
.
Circ. Res.
118
,
1480
-
1497
.
Martin
,
L. J.
,
Ramachandran
,
V.
,
Cripe
,
L. H.
,
Hinton
,
R. B.
,
Andelfinger
,
G.
,
Tabangin
,
M.
,
Shooner
,
K.
,
Keddache
,
M.
and
Benson
,
D. W.
(
2007
).
Evidence in favor of linkage to human chromosomal regions 18q, 5q and 13q for bicuspid aortic valve and associated cardiovascular malformations
.
Hum. Genet.
121
,
275
-
284
.
Masri
,
A.
,
Svensson
,
L. G.
,
Griffin
,
B. P.
and
Desai
,
M. Y.
(
2017
).
Contemporary natural history of bicuspid aortic valve disease: a systematic review
.
Heart
103
,
1323
-
1330
.
Mommersteeg
,
M. T.
,
Yeh
,
M. L.
,
Parnavelas
,
J. G.
and
Andrews
,
W. D.
(
2015
).
Disrupted Slit-Robo signalling results in membranous ventricular septum defects and bicuspid aortic valves
.
Cardiovasc. Res.
106
,
55
-
66
.
Morrisey
,
E. E.
,
Ip
,
H. S.
,
Lu
,
M. M.
and
Parmacek
,
M. S.
(
1996
).
GATA-6: a zinc finger transcription factor that is expressed in multiple cell lineages derived from lateral mesoderm
.
Dev. Biol.
177
,
309
-
322
.
Morrisey
,
E. E.
,
Tang
,
Z.
,
Sigrist
,
K.
,
Lu
,
M. M.
,
Jiang
,
F.
,
Ip
,
H. S.
and
Parmacek
,
M. S.
(
1998
).
GATA6 regulates HNF4 and is required for differentiation of visceral endoderm in the mouse embryo
.
Genes Dev.
12
,
3579
-
3590
.
Naumann
,
U.
,
Cameroni
,
E.
,
Pruenster
,
M.
,
Mahabaleshwar
,
H.
,
Raz
,
E.
,
Zerwes
,
H. G.
,
Rot
,
A.
and
Thelen
,
M.
(
2010
).
CXCR7 functions as a scavenger for CXCL12 and CXCL11
.
PLoS One
5
,
e9175
.
Phillips
,
H. M.
,
Mahendran
,
P.
,
Singh
,
E.
,
Anderson
,
R. H.
,
Chaudhry
,
B.
and
Henderson
,
D. J.
(
2013
).
Neural crest cells are required for correct positioning of the developing outflow cushions and pattern the arterial valve leaflets
.
Cardiovasc. Res.
99
,
452
-
460
.
Prakash
,
S. K.
,
Bossé
,
Y.
,
Muehlschlegel
,
J. D.
,
Michelena
,
H. I.
,
Limongelli
,
G.
,
Della Corte
,
A.
,
Pluchinotta
,
F. R.
,
Russo
,
M. G.
,
Evangelista
,
A.
,
Benson
,
D. W.
et al. 
(
2014
).
A roadmap to investigate the genetic basis of bicuspid aortic valve and its complications: insights from the International BAVCon (Bicuspid Aortic Valve Consortium)
.
J. Am. Coll. Cardiol.
64
,
832
-
839
.
Quaranta
,
R.
,
Fell
,
J.
,
Rühle
,
F.
,
Rao
,
J.
,
Piccini
,
I.
,
Araúzo-Bravo
,
M. J.
,
Verkerk
,
A. O.
,
Stoll
,
M.
and
Greber
,
B.
(
2018
).
Revised roles of ISL1 in a hES cell-based model of human heart chamber specification
.
Elife
7
,
e31706
.
Ridge
,
L. A.
,
Kewbank
,
D.
,
Schütz
,
D.
,
Stumm
,
R.
,
Scambler
,
P. J.
and
Ivins
,
S.
(
2021
).
Dual role for CXCL12 signaling in semilunar valve development
.
Cell Rep.
36
,
109610
.
Robinson
,
M. D.
,
Mccarthy
,
D. J.
and
Smyth
,
G. K.
(
2010
).
edgeR: a Bioconductor package for differential expression analysis of digital gene expression data
.
Bioinformatics
26
,
139
-
140
.
Rochais
,
F.
,
Mesbah
,
K.
and
Kelly
,
R. G.
(
2009
).
Signaling pathways controlling second heart field development
.
Circ. Res.
104
,
933
-
942
.
Roth
,
G. A.
,
Mensah
,
G. A.
,
Johnson
,
C. O.
,
Addolorato
,
G.
,
Ammirati
,
E.
,
Baddour
,
L. M.
,
Barengo
,
N. C.
,
Beaton
,
A. Z.
,
Benjamin
,
E. J.
,
Benziger
,
C. P.
et al. 
(
2020
).
Global burden of cardiovascular diseases and risk factors, 1990-2019: update from the GBD 2019 study
.
J. Am. Coll. Cardiol.
76
,
2982
-
3021
.
Runyan
,
R. B.
and
Markwald
,
R. R.
(
1983
).
Invasion of mesenchyme into three- dimensional collagen gels: a regional and temporal analysis of interaction in embryonic heart tissue
.
Dev. Biol.
95
,
108
-
114
.
Sharma
,
A.
,
Wasson
,
L. K.
,
Willcox
,
J. A.
,
Morton
,
S. U.
,
Gorham
,
J. M.
,
Delaughter
,
D. M.
,
Neyazi
,
M.
,
Schmid
,
M.
,
Agarwal
,
R.
,
Jang
,
M. Y.
et al. 
(
2020
).
GATA6 mutations in hiPSCs inform mechanisms for maldevelopment of the heart, pancreas, and diaphragm
.
Elife
9
,
e53278
.
Shi
,
Y.
,
Riese
,
D. J.
, 2nd
and
Shen
,
J.
(
2020
).
The role of the CXCL12/CXCR4/CXCR7 chemokine axis in cancer
.
Front. Pharmacol.
11
,
574667
.
Sierro
,
F.
,
Biben
,
C.
,
Martínez-Muñoz
,
L.
,
Mellado
,
M.
,
Ransohoff
,
R. M.
,
Li
,
M.
,
Woehl
,
B.
,
Leung
,
H.
,
Groom
,
J.
,
Batten
,
M.
et al. 
(
2007
).
Disrupted cardiac development but normal hematopoiesis in mice deficient in the second CXCL12/SDF-1 receptor, CXCR7
.
Proc. Natl. Acad. Sci. USA
104
,
14759
-
14764
.
Škorić-Milosavljević
,
D.
,
Tjong
,
F. V. Y.
,
Barc
,
J.
,
Backx
,
A.
,
Clur
,
S. B.
,
Van Spaendonck-Zwarts
,
K.
,
Oostra
,
R. J.
,
Lahrouchi
,
N.
,
Beekman
,
L.
,
Bökenkamp
,
R.
et al. 
(
2019
).
GATA6 mutations: characterization of two novel patients and a comprehensive overview of the GATA6 genotypic and phenotypic spectrum
.
Am. J. Med. Genet. A
179
,
1836
-
1845
.
Snarr
,
B. S.
,
Kern
,
C. B.
and
Wessels
,
A.
(
2008
).
Origin and fate of cardiac mesenchyme
.
Dev. Dyn.
237
,
2804
-
2819
.
Sodhi
,
C. P.
,
Li
,
J.
and
Duncan
,
S. A.
(
2006
).
Generation of mice harbouring a conditional loss-of-function allele of Gata6
.
BMC Dev. Biol.
6
,
19
.
Srivastava
,
D.
and
Olson
,
E. N.
(
2000
).
A genetic blueprint for cardiac development
.
Nature
407
,
221
-
226
.
Stanley
,
E. G.
,
Biben
,
C.
,
Elefanty
,
A.
,
Barnett
,
L.
,
Koentgen
,
F.
,
Robb
,
L.
and
Harvey
,
R. P.
(
2002
).
Efficient Cre-mediated deletion in cardiac progenitor cells conferred by a 3'UTR-ires-Cre allele of the homeobox gene Nkx2-5
.
Int. J. Dev. Biol.
46
,
431
-
439
.
Susaki
,
E. A.
,
Tainaka
,
K.
,
Perrin
,
D.
,
Kishino
,
F.
,
Tawara
,
T.
,
Watanabe
,
T. M.
,
Yokoyama
,
C.
,
Onoe
,
H.
,
Eguchi
,
M.
,
Yamaguchi
,
S.
et al. 
(
2014
).
Whole-brain imaging with single-cell resolution using chemical cocktails and computational analysis
.
Cell
157
,
726
-
739
.
Tastet
,
L.
,
Capoulade
,
R.
,
Clavel
,
M. A.
,
Larose
,
É.
,
Shen
,
M.
,
Dahou
,
A.
,
Arsenault
,
M.
,
Mathieu
,
P.
,
Bédard
,
É.
,
Dumesnil
,
J. G.
et al. 
(
2017
).
Systolic hypertension and progression of aortic valve calcification in patients with aortic stenosis: results from the PROGRESSA study
.
Eur. Heart J. Cardiovasc. Imaging
18
,
70
-
78
.
Thom
,
T.
,
Haase
,
N.
,
Rosamond
,
W.
,
Howard
,
V. J.
,
Rumsfeld
,
J.
,
Manolio
,
T.
,
Zheng
,
Z. J.
,
Flegal
,
K.
,
O'donnell
,
C.
,
Kittner
,
S.
et al. 
(
2006
).
Heart disease and stroke statistics--2006 update: a report from the American Heart Association Statistics Committee and Stroke Statistics Subcommittee
.
Circulation
113
,
e85
-
e151
.
Timmerman
,
L. A.
,
Grego-Bessa
,
J.
,
Raya
,
A.
,
Bertrán
,
E.
,
Pérez-Pomares
,
J. M.
,
Díez
,
J.
,
Aranda
,
S.
,
Palomo
,
S.
,
Mccormick
,
F.
,
Izpisúa-Belmonte
,
J. C.
et al. 
(
2004
).
Notch promotes epithelial-mesenchymal transition during cardiac development and oncogenic transformation
.
Genes Dev.
18
,
99
-
115
.
Van Berlo
,
J. H.
,
Elrod
,
J. W.
,
Van Den Hoogenhof
,
M. M.
,
York
,
A. J.
,
Aronow
,
B. J.
,
Duncan
,
S. A.
and
Molkentin
,
J. D.
(
2010
).
The transcription factor GATA-6 regulates pathological cardiac hypertrophy
.
Circ. Res.
107
,
1032
-
1040
.
Van Der Linde
,
D.
,
Konings
,
E. E.
,
Slager
,
M. A.
,
Witsenburg
,
M.
,
Helbing
,
W. A.
,
Takkenberg
,
J. J.
and
Roos-Hesselink
,
J. W.
(
2011
).
Birth prevalence of congenital heart disease worldwide: a systematic review and meta-analysis
.
J. Am. Coll. Cardiol.
58
,
2241
-
2247
.
Verzi
,
M. P.
,
Mcculley
,
D. J.
,
De Val
,
S.
,
Dodou
,
E.
and
Black
,
B. L.
(
2005
).
The right ventricle, outflow tract, and ventricular septum comprise a restricted expression domain within the secondary/anterior heart field
.
Dev. Biol.
287
,
134
-
145
.
Von Gise
,
A.
and
Pu
,
W. T.
(
2012
).
Endocardial and epicardial epithelial to mesenchymal transitions in heart development and disease
.
Circ. Res.
110
,
1628
-
1645
.
Waldo
,
K. L.
,
Hutson
,
M. R.
,
Ward
,
C. C.
,
Zdanowicz
,
M.
,
Stadt
,
H. A.
,
Kumiski
,
D.
,
Abu-Issa
,
R.
and
Kirby
,
M. L.
(
2005
).
Secondary heart field contributes myocardium and smooth muscle to the arterial pole of the developing heart
.
Dev. Biol.
281
,
78
-
90
.
Ward
,
C.
,
Stadt
,
H.
,
Hutson
,
M.
and
Kirby
,
M. L.
(
2005
).
Ablation of the secondary heart field leads to tetralogy of Fallot and pulmonary atresia
.
Dev. Biol.
284
,
72
-
83
.
Wijtmans
,
M.
,
Maussang
,
D.
,
Sirci
,
F.
,
Scholten
,
D. J.
,
Canals
,
M.
,
Mujić-Delić
,
A.
,
Chong
,
M.
,
Chatalic
,
K. L.
,
Custers
,
H.
,
Janssen
,
E.
et al. 
(
2012
).
Synthesis, modeling and functional activity of substituted styrene-amides as small-molecule CXCR7 agonists
.
Eur. J. Med. Chem.
51
,
184
-
192
.
Xu
,
Y. J.
,
Di
,
R. M.
,
Qiao
,
Q.
,
Li
,
X. M.
,
Huang
,
R. T.
,
Xue
,
S.
,
Liu
,
X. Y.
,
Wang
,
J.
and
Yang
,
Y. Q.
(
2018
).
GATA6 loss-of-function mutation contributes to congenital bicuspid aortic valve
.
Gene
663
,
115
-
120
.
Yassine
,
N. M.
,
Shahram
,
J. T.
and
Body
,
S. C.
(
2017
).
Pathogenic mechanisms of bicuspid aortic valve aortopathy
.
Front. Physiol.
8
,
687
.
Yu
,
S.
,
Crawford
,
D.
,
Tsuchihashi
,
T.
,
Behrens
,
T. W.
and
Srivastava
,
D.
(
2011
).
The chemokine receptor CXCR7 functions to regulate cardiac valve remodeling
.
Dev. Dyn.
240
,
384
-
393
.
Yuan
,
S. M.
and
Jing
,
H.
(
2010
).
The bicuspid aortic valve and related disorders
.
Sao Paulo Med. J.
128
,
296
-
301
.
Zhou
,
B.
,
Von Gise
,
A.
,
Ma
,
Q.
,
Rivera-Feliciano
,
J.
and
Pu
,
W. T.
(
2008
).
Nkx2-5- and Isl1-expressing cardiac progenitors contribute to proepicardium
.
Biochem. Biophys. Res. Commun.
375
,
450
-
453
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.