Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Reduced dosage of β-catenin provides significant rescue of cardiac outflow tract anomalies in a Tbx1 conditional null mouse model of 22q11.2 deletion syndrome

Abstract

The 22q11.2 deletion syndrome (22q11.2DS; velo-cardio-facial syndrome; DiGeorge syndrome) is a congenital anomaly disorder in which haploinsufficiency of TBX1, encoding a T-box transcription factor, is the major candidate for cardiac outflow tract (OFT) malformations. Inactivation of Tbx1 in the anterior heart field (AHF) mesoderm in the mouse results in premature expression of pro-differentiation genes and a persistent truncus arteriosus (PTA) in which septation does not form between the aorta and pulmonary trunk. Canonical Wnt/β-catenin has major roles in cardiac OFT development that may act upstream of Tbx1. Consistent with an antagonistic relationship, we found the opposite gene expression changes occurred in the AHF in β-catenin loss of function embryos compared to Tbx1 loss of function embryos, providing an opportunity to test for genetic rescue. When both alleles of Tbx1 and one allele of β-catenin were inactivated in the Mef2c-AHF-Cre domain, 61% of them (n = 34) showed partial or complete rescue of the PTA defect. Upregulated genes that were oppositely changed in expression in individual mutant embryos were normalized in significantly rescued embryos. Further, β-catenin was increased in expression when Tbx1 was inactivated, suggesting that there may be a negative feedback loop between canonical Wnt and Tbx1 in the AHF to allow the formation of the OFT. We suggest that alteration of this balance may contribute to variable expressivity in 22q11.2DS.

Author summary

To understand the genetic relationship between Tbx1 and canonical Wnt/β-catenin, we performed gene expression profiling and genetic rescue experiments. We found that Tbx1 and β-catenin may provide a negative feedback loop to restrict premature differentiation in the anterior heart field. This is relevant to understanding the basis of variable expressivity of 22q11.2DS, caused by haploinsufficiency of TBX1.

Citation: Racedo SE, Hasten E, Lin M, Devakanmalai GS, Guo T, Ozbudak EM, et al. (2017) Reduced dosage of β-catenin provides significant rescue of cardiac outflow tract anomalies in a Tbx1 conditional null mouse model of 22q11.2 deletion syndrome. PLoS Genet 13(3): e1006687. https://doi.org/10.1371/journal.pgen.1006687

Editor: Gregory S. Barsh, Stanford University School of Medicine, UNITED STATES

Received: January 31, 2017; Accepted: March 13, 2017; Published: March 27, 2017

Copyright: © 2017 Racedo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Microarray data have been deposited to the GEO database (accession number: GSE78125).

Funding: This work was supported by the Foundation Leducq. It was also supported by NIH grants, P01HD070454 (BEM), R56HL129807 (CC, DZ and BEM) and R01 GM111987 01 (EMO) and an AHA grant 12POST9100003 (SER). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The 22q11.2 deletion syndrome (22q11.2DS), also known as velo-cardio-facial syndrome (MIM# 192430) or DiGeorge syndrome (MIM# 188400) is a congenital malformation disorder that is caused by a hemizygous 1.5–3 million base pair (Mb) deletion of chromosome 22q11.2. It occurs with a frequency of 1:1,000 fetuses [1] and 1:4,000 live births [2]. Approximately 60–70% of affected 22q11.2DS individuals have congenital heart disease (CHD) due to malformations of the aortic arch and/or cardiac outflow tract [3]. There are over 46 known coding genes in the 3 Mb region, including TBX1 (T-box 1; MIM# 602054), encoding a T-box containing transcription factor [4]. TBX1 has been considered the strongest candidate gene for CHD, based upon studies of mouse models [5–7] and discovery of mutations in some non-deleted patients [8, 9]. The basis of variable phenotypic expression is under intense investigation. Understanding responsible genetic factors upstream and downstream of TBX1 is necessary to test for relevancy as modifiers in human 22q11.2DS patients. We are taking mouse genetic approaches to identify genes and networks that may act as modifiers.

Tbx1 heterozygous mice have mild aortic arch anomalies or ventricular septal defects, at reduced penetrance, while all homozygous null mutant mice die at birth and have a persistent truncus arteriosus (PTA), which is the most serious heart defect that occurs in 22q11.2DS patients [5–7]. In mammals, Tbx1 is expressed strongly in the embryonic pharyngeal apparatus, but not the heart tube itself suggesting that its critical functions are in this tissue [4].

In the early vertebrate embryo, the heart forms as a bilateral cardiac crescent of mesodermal cells, termed the first heart field that fuses to form the primitive heart tube [10, 11]. Additional mesodermal cells derived from the pharyngeal apparatus, referred to as the second heart field (SHF) migrates and helps to expand the heart tube in both directions [12] [13] [13–16]. These cells remain in a progenitor state, allowing them to migrate and build the length of the heart tube, where they differentiate into smooth and cardiac muscle and endothelial cells [17, 18]. The SHF itself, can be further subdivided to the anterior heart field (AHF or anterior SHF) forming the cardiac OFT and right ventricle as well as the posterior SHF forming the inflow tract, respectively, based upon gene expression and cell lineage studies [19–21]. Of interest, Tbx1 is strongly expressed in the pharyngeal mesoderm, including the AHF, but it is not noticeably expressed in the posterior SHF or heart tube [22–24]. One of the key functions of AHF cells is to maintain a progenitor cell state and to prevent premature differentiation. [25] Gene expression profiling of the AHF, within pharyngeal arches two to six, in Tbx1-/- embryos versus wild type littermates [24] and embryonic stem cell lineage studies [22], suggest that Tbx1 serves to restrict premature differentiation of the pharyngeal mesoderm, so as to allow the OFT to elongate properly [25]. However, the tissue specificity and key molecular mechanisms are not well defined.

The basis for premature differentiation in the AHF in Tbx1 mutant embryos is unknown. Major signaling pathways likely have a role in this process. The canonical Wnt signaling pathway is mediated by β-catenin, which has critical functions in most aspects of embryonic development. There are multiphasic functions of Wnt/β-catenin in the pharyngeal mesoderm required for heart development [26]. Several years ago, it was shown that canonical Wnt/β-catenin has a major role in the AHF in forming the cardiac OFT [27]. Further, one study showed that increased or decreased Wnt/β-catenin in the pharyngeal mesenchyme (DermoCre) resulted in a decrease or increase in Tbx1 expression, implicating antagonistic functions upstream of Tbx1 [28]. However, genetic interaction studies were not explored nor were gene expression profiling performed to understand possible molecular connections. Such studies would provide possible modifier genes to investigate in human 22q11.2DS to understand its variable expressivity. In this report we performed genetic rescue experiments between Tbx1 and β-catenin in the AHF, using mouse models.

Results

Constitutive β-catenin expression in the AHF promotes differentiation

Wnt/β-catenin and Tbx1 are expressed in the opposite domains of the SHF, with Tbx1 higher in the AHF and Wnt/β-catenin higher in the posterior SHF, as denoted by Wnt2 and Mef2c-AHF-Cre [18] lineage compared to canonical Wnt signaling (Fig 1A–1E). We were interested in further exploring the function of β-catenin when completely diminished (Mef2c-AHF-Cre/+;β-cateninflox/flox, referred to as β-cat LOF [29]) or constitutively active (Mef2c-AHF-Cre/+;β-cateninEx3/+ [30], referred to as β-cat GOF) in the AHF. To identify downstream genes affected by these changes, gene expression profiling was performed on the distal pharyngeal apparatus containing the AHF micro-dissected from β-cat LOF and β-cat GOF embryos at E9.5 (Fig 1F–1H). Note that the dissection of the AHF did not include the heart tube. In order to highlight the genes with the greatest fold change, we created a dot plot of log2 fold changes (Fig 1I). Loss of both β-catenin alleles in the Mef2c-AHF-Cre domain resulted in strongly reduced expression of muscle structural genes in the AHF, while constitutive activation of β-catenin in this domain, had the opposite effect and caused a strong increase in expression of the same genes in the AHF (Fig 1I). This increase was strikingly similar to that in the AHF of global Tbx1 null mutant embryos that were previously reported [22, 24, 31]. We then examined cardiac phenotypes upon inactivation of Tbx1 in the Mef2c-AHF-Cre lineage to determine whether Tbx1 had a specific role in the AHF.

thumbnail
Download:
Fig 1. Constitutive β-catenin expression in the AHF promotes differentiation.

(A) Whole mount in situ hybridization (WMISH) of Wnt2 in a wild type (Tbx1+/+) embryo at embryonic day 9.5 (E9.5) and a sagittal section from the heart region of the embryo is shown in B. (C) Mef2c-AHF-Cre lineage in a wild type embryo in a sagittal section. (D) Visualization of β-catenin signaling using a TCF/Lef:H2B-GFP reporter allele in a sagittal section of a wild type embryo. (E) WMISH of Tbx1 in a wild type embryo. (F) Mef2c-AHF-Cre lineage shown by green fluorescence in a wild type embryo (Mef2c-AHF-Cre/+;ROSA26-GFP f/+) at E9.5. Numbers denote the pharyngeal arches 1, 2 and 3. (G) Whole mount embryo showing the region that has been dissected for the experiments (rectangle in dorsal to the heart). (H) The micro-dissected region containing the Mef2c-AHF-Cre lineage is shown from a frontal view. (I) Comparison of global gene expression changes in micro-dissected tissues of β-catenin GOF and β-catenin LOF embryos at E9.5. Differentially expressed genes (p < 0.05 and FC > 1.5) in at least one of the two comparisons, β-catenin LOF vs control (x-axis) or β-catenin GOF vs control (y-axis), were plotted. Red dots denote cardiac differentiation genes. Abbreviations: heart (H), pharyngeal arch (PA), outflow tract (OFT), ventricle (V), atrium (A), FC = fold change.

https://doi.org/10.1371/journal.pgen.1006687.g001

Persistent truncus arteriosus in Tbx1 LOF embryos

To determine a specific role of Tbx1 in the Mef2c-AHF-Cre domain, we generated two different genotypes, Mef2c-AHF-Cre/+;Tbx1f/- and Mef2c-AHF-Cre/+;Tbx1f/f. Embryos at E14.5 with both genotypes had a persistent truncus arteriosus (PTA) with complete penetrance (n = 50; Fig 2A). Most but not all had an accompanying ventricular septal defect (VSD; n = 30; Fig 2B and 2C), in contrast to Tbx1-/- mutant embryos, which all has a PTA with a VSD. The PTA was observed as early as E12.5 (S1 Fig). Due to the similarity in phenotype (Fig 2A), the two genotypes were combined and further referred to as Tbx1 LOF. In Tbx1 LOF embryos at E9.5, the pharyngeal apparatus and individual arches within appeared grossly normal (Fig 2F–2G) compared to control littermates (Fig 2D and 2E). This is distinctly different as compared to the global Tbx1-/- null mutant embryos or mesoderm specific Tbx1 conditional loss of function embryos at this stage [32] that have a severely hypoplastic distal pharyngeal apparatus. This rules out extreme morphology defects, such as absence of neural crest cell populations, as being responsible for the presence of a PTA in Tbx1 LOF embryos. We also performed lineage tracing (Fig 2D–2G) and observed that the Mef2c-AHF-Cre lineage in the AHF was only slightly reduced in Tbx1 LOF embryos versus control littermates at E9.5 (Fig 2H). By in situ hybridization analysis, Tbx1 expression was greatly reduced in Tbx1 LOF embryos (S1 Fig) and this was confirmed by qRT-PCR (Fig 2I). Cell proliferation and apoptosis in the Tbx1 LOF versus control embryos did not show any significant difference in the Mef2c-AHF-Cre lineage in the AHF region between groups at E9.5 (S2 and S3 Figs). This is different than what was previously found for Tbx1-/- [33] or Nkx2-5Cre [22] conditional mutant embryos, which have changes in proliferation and apoptosis. We suggest the improved appearance of the distal pharyngeal apparatus in Tbx1 LOF embryos is due to differences in the Mef2c-AHF-Cre recombination domain. In relation to β-catenin, we noted a slight decrease in Tbx1 expression in β-cat GOF embryos (S4 Fig). This was consistent, although not as dramatic, as was found previously using a broader mesenchymal Cre driver (DermoCre) [28]. We found β-catenin mRNA is significantly increased in expression in the AHF of Tbx1 LOF embryos by qRT-PCR (Fig 2I).

thumbnail
Download:
Fig 2. Congenital heart defects in Tbx1 LOF embryos.

(A) Heart phenotype analysis of Tbx1 LOF embryos at E14.5 generated from two different crosses. PTA-VSD refers to a PTA associated with a ventricular septal defect (VSD) while PTA refers to hearts that did not show a VSD. Partial septation in Tbx1 LOF embryos means a PTA with presence of a short septum at some level of the OFT and complete septation between the ventricles. N = total number of hearts observed per group. Significance between Tbx1 LOF and controls calculated by Fisher’s exact test (p < 0.001). Note that Mef2c-AHF-Cre;Tbx1flox/flox embryos had additional phenotypes (three with double outlet right ventricle, DORV and one with tetralogy of Fallot, TOF, as indicated); (B) H&E histological sections of the heart of a control embryo at E14.5, with a typical ventricular septum is shown in the inset on the lower right part of the image. (C) Tbx1 LOF embryo with a PTA-VSD. (D) Mef2c-AHF-Cre lineage tracing by using a GFP reporter allele in a control embryo at E9.5 and a representative sagittal section is shown in E. (F) Mef2c-AHF-Cre lineage tracing in a Tbx1 LOF embryo and a representative sagittal section of the embryo is shown in G. DAPI fluorescent stain to visualize nuclei and identify the tissue is shown in blue. (H) Mef2c-AHF-Cre lineage quantification from the area shown in the inset in G for control and Tbx1 LOF embryos. (I) Detection of Tbx1 and β-catenin by qRT-PCR in control and Tbx1 LOF embryos. Statistical significance of the difference in gene expression was estimated using two-tailed t-test, FC = fold change, p values < 0.05. Error bars = standard deviation (SD). Abbreviations: aorta (Ao), pulmonary trunk (PT), left atrium (LA), right atrium (RA), left ventricle (LV), right ventricle (RV), pharyngeal arch (PA), outflow tract (OFT), 1, 2 and 3 indicate the first, second and third pharyngeal arches (PA), respectively.

https://doi.org/10.1371/journal.pgen.1006687.g002

β-catenin and Tbx1 conditional mutants have opposite roles in AHF muscle cell differentiation

As for β-catenin gain of function in the AHF, we were interested in determining the function of loss and gain of Tbx1 in the AHF. We previously generated a tissue specific constitutively expressing Tbx1 gene [34]. Homozygous mice were crossed with Mef2c-AHF-Cre mice to overexpress Tbx1 in the same domain as other alleles, and the embryos are referred to as Tbx1 GOF. Gene expression profiling of Tbx1 LOF and GOF embryos was performed of AHF tissue at E9.5, to test whether loss or gain of Tbx1 would have opposing effects on muscle structural protein differentiation genes and to compare with findings of β-catenin loss and gain mutant embryos (Fig 3A–3C). The dot plots of global gene expression changes in the AHF between β-cat GOF versus Tbx1 LOF embryos showed increase in gene expression in the same direction (Fig 3A). The genes with the largest increase were the muscle structural protein genes. Similarly, β-cat LOF versus Tbx1 GOF showed the same strong decrease of expression of muscle differentiation genes. The genes with the largest decrease were the muscle structural protein genes. These results provide functional genetic insight as to the previously implicated antagonistic relationship between Tbx1 and β-catenin in the pharyngeal mesenchyme [28], that they perhaps are needed to balance cell differentiation. Some of the genes were tested by qRT-PCR for Tbx1 LOF and β-cat LOF embryos (Fig 3C). We found top genes decreased in expression in Tbx1 LOF embryos were not generally decreased in β-cat LOF embryos as top genes that were increased in expression. Tbx1 was slightly increased in expression in β-cat LOF embryos (Fig 3C).

thumbnail
Download:
Fig 3. Opposing β-catenin and Tbx1 conditional mutants have same effect on expression of pro-differentiation genes in the AHF.

(A) Comparison of global gene expression changes in the micro-dissected AHF of β-catenin GOF and Tbx1 LOF embryos at E9.5. Plotted are differentially expressed genes (p < 0.05 and FC > 1.5) in at least one of the two comparisons, Tbx1 LOF vs controls (x-axis) or β-catenin GOF vs controls (y-axis). (B) Comparison of global gene expression changes in the micro-dissected AHF from Tbx1 GOF and β-catenin LOF embryos at E9.5. Plotted are differentially expressed genes (p < 0.05 and FC > 1.5) in at least one of the two comparisons, Tbx1 GOF vs controls (x-axis) or β-catenin LOF vs controls (y-axis). Red dots denote cardiac differentiation genes. (C) qRT-PCR analysis of Tbx1 LOF versus β-catenin LOF (FC = Fold Change) of selected genes demonstrating the opposite gene expression changes in these mutant embryos.

https://doi.org/10.1371/journal.pgen.1006687.g003

Reduced β-catenin dosage significantly rescues heart defects in Tbx1 LOF embryos

Based upon the opposing gene expression changes between Tbx1 and β-catenin in the AHF, that also included β-catenin mRNA itself in Tbx1 LOF embryos (Fig 2I), we tested whether we could rescue heart defects in the Tbx1 conditional loss of function mutant embryos by inactivating one allele of β-catenin in the Mef2c-AHF-Cre domain (Mef2c-AHF-Cre/+;Tbx1f/- and Mef2c-AHF-Cre/+;Tbx1f/f). Details of the background and crosses are provided in the Methods section and details of the control genotypes are provided in S1 Table. Inactivation of one allele of β-catenin in the Mef2c-AHF-Cre domain did not result in any cardiovascular defects (S1 Table and [27]). Significant rescue (p < 0.001, Fisher’s exact test) was obtained in both sets of double mutant rescue genotype embryos (Fig 4). Upon combining all double mutant embryos together (n = 56), a total of 61% (n = 34/56) showed some rescue of the PTA phenotype (Fig 4). Specifically, complete distal OFT and partial proximal OFT septation and/or complete septation between the ventricles were present in these hearts (Fig 5A–5D). Ten percent showed complete rescue. Additional and more posterior sections can be found in S5 Fig Lineage tracing of the double mutant rescue genotyped embryos showed no significant difference in the number of cells in the AHF, of Mef2c-AHF-Cre lineage compared to the control or Tbx1 LOF embryos at E9.5 (Fig 5E). Finally, qRT-PCR was performed and Tbx1 and β-catenin mRNA expression were reduced in the AHF of these embryos compared to the control (Fig 5F).

thumbnail
Download:
Fig 4. Phenotype analysis of significantly rescued Tbx1 LOF embryos.

Phenotypes in embryos in which one β-catenin loss of function allele to either Mef2c-AHF-Cre/+;Tbx1f/- embryos (A) or to Mef2c-AHF-Cre/+;Tbx1f/f embryos (B) was done to lower the dosage of β-catenin within the AHF. Significant rescue (p < 0.001, Fisher’s exact test) was obtained in both sets of double mutant embryos. Middle panels show the percent within the groups with PTA (C and C’), DORVor TOF associated with a VSD or PTA without VSD (dark grey) and those ones with partial septation or a normal heart (light grey). Partial septation in Tbx1 LOF embryos means a PTA with presence of a short septum at some level of the OFT (D and D’) and complete septation between the ventricles (8/50) while partial septation in the rescue genotype embryos (34/56) means a range of noticeably less severe phenotypes including: a longer partial OFT septation (E and E’) and either complete septation between ventricles (15/56) or VSD (10/56) or normal OFT (F and F’) with a VSD (4/56) or normal heart [5/56] (F and F’). Abbreviations: aorta (Ao), pulmonary trunk (PT), right ventricle (RV), outflow tract (OFT), PTA (persistent truncus arteriosus).

https://doi.org/10.1371/journal.pgen.1006687.g004

thumbnail
Download:
Fig 5. Histology analysis of representative embryos with the “rescue” genotype.

Transverse H&E histological sections of hearts from E14.5 significantly rescued embryos (Tbx1 LOF with loss of one allele of β-catenin in the Mef2c-AHF-Cre domain). (A) An embryonic heart with a thin septum formed that separates the Ao and the PT (top black arrow), and a small VSD compared to the usual PTA-VSD in Tbx1 LOF embryos (Fig 2C). (B) Example showing a DORV (middle black arrow), with a separate Ao and PT; normal ventricular septation is present. (C) Example of a heart showing a partially rescued septation between the Ao and PT (middle black arrow) and normal septation between the two ventricles (right black arrow). (D) Rescued septation between the Ao and PT. (E) Mef2c-AHF-Cre lineage quantification from the same area shown in Fig 2G for control, Tbx1 LOF and rescued embryos. (F) Detection of Tbx1 and β-catenin by qRT-PCR in control, Tbx1 LOF and rescued embryos. Statistical significance of the difference in gene expression was estimated using two-tailed t-test, FC = fold change, p values < 0.05. Error bars = standard deviation (SD). Abbreviations: RA = right atrium, RV = right ventricle, LA = left atrium, LV = left ventricle, Ao = Aorta, PT = pulmonary trunk, VS = ventricular septum, PTA = persistent truncus arteriosus, VSD = ventricular septal defect, DORV = double outlet right ventricle, OFT = outflow tract.

https://doi.org/10.1371/journal.pgen.1006687.g005

Genes for cardiomyocyte differentiation are normalized in rescued embryos

Since we identified the greatest increase of expression in Tbx1 LOF and decrease in β-cat LOF embryos pertaining to muscle differentiation genes, we tested if there is global normalization of expression in the embryos of the double mutant, rescue genotype. For this test, gene expression profiling was performed on these embryos, in the same way for the individual mutant embryos and we found this to be the case. Expression of genes with greatest increase in Tbx1 LOF embryos (>1.3 fold), primarily the differentiation genes and greatest decrease in β-cat LOF embryos were largely normalized in rescued embryos (Fig 6A). However, we did not observe this same strong finding for genes increased in β-cat LOF embryos.

thumbnail
Download:
Fig 6. Gene expression of cardiac muscle differentiation genes are normalized in rescued embryos.

(A) Comparison of global expression changes. Differentially expressed genes (p < 0.05; n = 3,636) in either Tbx1 LOF or β-catenin LOF were sorted by their FC between Tbx1 LOF and controls (x-axis). The sorted genes were grouped (50 genes per group) and then the average expression changes for each group in the Tbx1 LOF, β-catenin LOF or rescue embryos (vs controls) were plotted in the y-axis. (B) Heatmap showing the expression of selected key genes increased in expression in Tbx1 LOF embryos but decreased in expression in β-catenin LOF embryos.

https://doi.org/10.1371/journal.pgen.1006687.g006

As mentioned, most of the genes with the strongest increase of expression in Tbx1 LOF and decrease in β-cat LOF embryos that were normalized (p<0.01) in rescued embryos, were genes that encode smooth or cardiac muscle genes (Figs 3C and 6B). This also included major transcription factors such as Pitx2, Tbx5, Gata4 and Gata6, that are required for cardiac muscle differentiation [35–37]. The canonical Wnt gene, Wnt2, showed a similar pattern (Figs 3C and 6B). A full heatmap of the experiment is shown in S6 Fig Some additional genes of note, increased in expression in Tbx1 LOF embryos by gene profiling include Myocd, Bmp2, Bmp10, Erbb4 and Sfrp5, which were oppositely affected in β-cat LOF embryos, and normalized in rescued embryos (Fig 6B and S6 Fig). This data supports the idea that pro-differentiation by canonical Wnt/β-catenin might be modulated by Tbx1 in the Mef2c-AHF-Cre lineage.

Not all genes with increase in expression in Tbx1 LOF embryos and decrease in β-cat LOF showed normalization in rescued embryos (Hand1, Zfpm2, Smarcd3 and Tbx20; Fig 6B and S7 Fig). Further, genes reduced in expression in both types of LOF mutants (S7 Fig) might not be relevant for the observed rescued phenotype since they were not normalized in rescue genotyped embryos. This suggests that other pathways are required for OFT development, and explains, in part, why complete rescue did not occur. Nonetheless it provides insights as to the nature of the relationship of the two genes, Tbx1 and β-catenin as well as their independent functions.

Discussion

Tbx1 and β-catenin in the AHF

Loss of β-catenin using various pharyngeal mesoderm engineered Cre drivers, including Mesp1-Cre [38], Nkx2-5-Cre [39], Isl1-Cre [40, 41] and Mef2c-AHF-Cre [27] results in embryonic lethality due to the presence of cardiac outflow tract defects. The mechanisms mediating these abnormalities, in particularly within the pharyngeal mesoderm of the AHF, have not been well defined. This is especially important because there are many divergent and distinctive functions of β-catenin during cardiac development [42–45] [45] [38, 46, 47]. Our interest was to follow up on a previous study in which Tbx1 expression was affected oppositely by loss or constitutively active β-catenin in the pharyngeal mesenchyme using a mesenchymal Cre, termed DermoCre [28]. Based upon this finding, we investigated the two genes in the AHF tissue at stage E9.5, when the heart tube is elongating. We found that β-catenin promotes muscle differentiation in the AHF. We also found that Tbx1 and Wnt/β-catenin act antagonistically to provide a balance of expression of pro-differentiation genes in the AHF that may be required for cardiac outflow tract development. This sheds new light onto the importance of the two genes in heart development as outlined in the model shown in Fig 7.

thumbnail
Download:
Fig 7. Model for Tbx1 and β-catenin function in the SHF.

In the model, the triangle represents the Mef2c-AHF-Cre lineage and the Tbx1 expression pattern, before migrating into the heart tube at E9.5. Left panel: Tbx1 expression is strongest in the AHF and weakest in the posterior SHF (pSHF) while Wnt/β-catenin expression is opposite. Left panel depicts a possible double negative feedback loop between the two genes in the SHF, required for normal cardiac outflow tract (OFT) formation. Middle panel shows increased differentiation in the AHF when Tbx1 is inactivated or β-catenin is constitutively active in the AHF. This results in premature differentiation within this tissue. Right panel depicts the rescue genotype in which both alleles of Tbx1 and one allele of β-catenin was inactivated in the AHF. Significant rescue of heart defects was obtained. Abbreviations: A = anterior, P = posterior, OFT = outflow tract.

https://doi.org/10.1371/journal.pgen.1006687.g007

In the model in Fig 7, we illustrate the Tbx1 expression domain in the SHF as a triangle, with the strongest expression anteriorly, in the AHF tissue and weakest in the posterior SHF. On the other hand, Wnt/β-catenin expression and function is strongest in the posterior SHF and weakest in the AHF, at E9.5 during mouse embryogenesis. In the panel on the left, we created a simple negative feedback loop, which is consistent with our findings in this study and previous findings using DermoCre [28]. The center panel of the model illustrates the situation when Tbx1 is inactivated or β-catenin is constitutively active in the AHF. Here, in these embryos differentiation occurs prematurely in the AHF, prior to reaching the heart tube (Fig 7, middle panel). This results in impaired cardiac outflow tract development.

In our study, we found that loss of Tbx1 in the Mef2c-AHF-Cre domain, along with loss of one allele of β-catenin provided significant rescue of heart defects (Fig 7, right panel). This supports the importance of their interaction. However, rescue is not complete. An explanation for this is that β-catenin is only partially diminished when one allele is inactivated, such that complete normalization is not possible. Another explanation is that both Tbx1 and Wnt/β-catenin act in many complex pathways in the AHF at this time point (E9.5), for which only the overlapping functions were normalized in rescued embryos [48–51]. Further, genes changing in expression at E9.5 may only be partially reflected in PTA defects observed at E12.5. We also note that the defects in the cardiac outflow tract between Tbx1 loss and β-catenin gain in the Mef2c-AHF-Cre domain are different [27], supporting this idea. In particular, neonatal lethality occurs in Tbx1 LOF embryos due to the presence of a PTA, while gain of β-catenin results in mid-gestational lethality with a short, hypercellular outflow tract.

Molecular mechanisms for Tbx1 function

One of the main functions of Tbx1 in the AHF is to maintain a progenitor state and restrict premature differentiation prior to reaching the elongating heart tube [22, 24]. Supporting this idea, constitutive overexpression of Tbx1 in the Mef2c-AHF-Cre domain results in a decrease in expression of muscle differentiation genes [22]. Based upon our gene expression profiling data, we suggest that Tbx1 may directly or indirectly, repress expression of key transcription factors that regulate this process i.e., it maintains the AHF cell fate. We suggest that there is a small decrease in AHF cell numbers but more importantly, a change in cell fate. Lack of observable morphology defects in the distal pharyngeal apparatus and lack of significant change in proliferation or apoptosis of the AHF progeny at E9.5 support this. Interestingly, in Tbx1 LOF mutant embryos, we found an increase in expression of genes required for cardiomyocyte specification, such as Gata4, Tbx5 and Smarcd3 (Baf60c) [52–55] [56–61] [35, 62–65]. Tbx5 and Gata4 proteins are co-expressed and physically interact to regulate expression of downstream muscle structural protein genes [52] [53]. The combination of Tbx5, Gata4 and Smarcd3 are sufficient to differentiate mouse embryonic mesoderm to beating cardiomyocytes [54]. Another intermediate protein is Serum Response Factor (SRF), which directly promotes expression of genes encoding muscle structural proteins that were found increased in Tbx1 null mutant embryos [22]. Inactivation of Tbx1 resulted in expansion of expression of SRF protein but not mRNA [22]. Similarly, we did not find Srf expression levels altered in Tbx1 LOF embryos. Of interest, the above transcription factors may interact with SRF protein to induce differentiation [66], supporting a continued role of SRF in Tbx1 biology [22] [36, 67]. Of interest, the Wnt2, Tbx5, Tbx20, Gata4 and Gata6 genes are expressed and have function in the posterior SHF for formation of the inflow tract. It is not yet known if any are directly or indirectly regulated by Tbx1.

An additional role of Tbx1 may be to restrict posterior SHF fate in the AHF so as to maintain the appropriate sub-populations within the SHF for proper heart development. We previously found expression of these posterior SHF genes were greatly expanded in the AHF tissue in Tbx1 global null mutant embryos [24] and were increased in the same tissue in the conditional mutant embryos by qRT-PCR. It was previously found that Wnt2 and Gata6 act in the same genetic pathway in the posterior SHF during heart development and when inactivated cause atrial septal defects among other anomalies [68, 69]. We observed an increase in expression, but did not identify atrial septal defects. Since we did not observe a severe morphological defect in Tbx1 LOF embryos at E9.5, we suggest that some of these molecular changes will then affect later development.

Relation of mouse developmental biology to 22q11.2 deletion syndrome

One of the challenges in human genetics is to identify risk factors of complex traits, such as congenital heart disease [70, 71]. The 22q11.2DS, although rare in the general population, offers a relatively homogenous cohort to investigate the basis of variable phenotypic heterogeneity among affected individuals. Rare deleterious DNA variants altered in expression in Tbx1 LOF embryos and acting antagonistically to canonical Wnt/β-catenin, might act as genetic modifiers of CHD in 22q11.2DS. Examination of whole exome sequence of 22q11.2DS subjects [72] is underway with a larger cohort, to identify such variants connected to Tbx1 and Wnt/β-catenin gene networks or pathways needed to provide proper balance of critical cell fate choices. The work here provides a basis in the near future to translate efforts to studies of human subjects.

Summary

In this study, we showed that inactivation of Tbx1 in the AHF using Mef2c-AHF-Cre allele, results in a PTA that is also observed in the most seriously affected 22q11.2DS (velo-cardio-facial/DiGeorge syndrome) patients. The PTA defect in Tbx1 conditional loss of function mutant embryos, was partially, but significantly rescued by decreasing one allele of the β-catenin gene in this domain, and this also resulted in a normalization of gene expression changes specifically for muscle differentiation but not necessarily for other classes of genes. Thus, we conclude that Tbx1 in the Mef2c-AHF-Cre domain acts antagonistically with Wnt/β-catenin in the SHF to modulate differentiation prior to entering the heart tube.

Methods

Mouse mutant alleles

The following mouse mutant alleles used in this study have been previously described:

Tbx1+/- [7], Tbx1f/+ (flox = f) [73], Tbx1-GFP [34], β-cateninf/+ and β-cateninE3/+ [29, 30], Mef2c-AHF-Cre/+ [18], ROSA26-GFPf/+ (RCE:loxP)[74] and Wnt/β-catenin signaling reporter mice (Tg(TCF/Lef1-HIST1H2BB/EGFP)61Hadj/J; TCF/Lef:H2B-GFP [75]. To generate Mef2c-AHF-Cre/+;Tbx1f/- mutant embryos (Tbx1 LOF), Mef2c-AHF-Cre/+ transgenic male mice were crossed to Tbx1+/- mice to obtain male Mef2c-AHF-Cre/+;Tbx1+/- mice that were then crossed with Tbx1f/f females. Alternatively, to generate Mef2c-AHF-Cre/+;Tbx1f/f mutant embryos, Mef2c-AHF-Cre/+ transgenic male mice were crossed to Tbx1f/f mice to obtain male Mef2c-AHF-Cre/+;Tbx1f/+ mice, and these were then crossed with Tbx1f/f females. Wild type and Me2fc-AHF-Cre/+;Tbx1f/+ littermates were used as controls for the experiments (First Tbx1 LOF and rescue crosses, S1 Table). Tbx1 gain of function embryos (Tbx1 GOF) were generated by crossing male Mef2c-AHF-Cre/+ mice with Tbx1-GFPf/f females. To generate Mef2c-AHF-Cre/+;β-cateninf/f mutant embryos (β-cat LOF), male Mef2c-AHF-Cre/+ transgenic mice were crossed to β-cateninf/f females to obtain male Mef2c-AHF-Cre/+;β-cateninf/+ mice that were then crossed with β-cateninf/f females. β-catenin gain of function (β-cat GOF) embryos i.e, Mef2c-AHF-Cre/+;β-cateninE3/+ or male Mef2c-AHF-Cre/+ transgenic mice were crossed to β-cateninE3/E3 females. Double mutant embryos were generated by addition of one copy of the β-cateninflox allele to Tbx1 LOF embryos resulting in what we denote as rescue genotyped embryos. In this case, the females used for the experimental crosses were of the Tbx1f/f;β-catenin f/f genotype, which have been maintained as an inbred line deriving from a mixed C57Bl/6; Swiss Webster background. The reporter ROSA26-GFP f/+ allele was added to the Tbx1f/f and Tbx1f/f;β-catenin f/f lines when visualizing Mef2c-AHF-Cre lineage. To evaluate Wnt/β-catenin signaling in wild type embryos, TCF/Lef:H2B-GFP/+ reporter mice were used.

The Mef2c-AHF-Cre/+;Tbx1f/+ and the Mef2c-AHF-Cre/+;Tbx1+/- mice are congenic in Swiss Webster. The Tbx1f/f;β-cateninf/f mice are in an inbred line, as above. The Mef2c-AHF-Cre/+;Tbx1+/- x Tbx1f/f; β-cateninf/f crosses were performed 2 years before the Mef2c-AHF-Cre/+;Tbx1f/+ x Tbx1f/f; β-cateninf/f crosses. The Tbx1f/f and ROSA26-GFPf/+ lines are congenic in Swiss Webster. The β-cateninE3/+ and β-cateninf/f mice were in a mixed C57Bl/6; Swiss Webster background. To exclude the possibility that a strain background might affect the possible rescue by β-catenin LOF allele, half of both Tbx1 LOF and the rescue genotyped embryos were generated by using Tbx1f/+;β-cateninf/+ females (second Tbx1 LOF and rescue crosses; S1 Table). Here, both Tbx1 LOF and the rescue genotyped embryos were littermates. The PCR strategies for mouse genotyping have been described in the original reports and are available upon request. All experiments including mice were carried out according to regulatory standards defined by the NIH and the Institute for Animal Studies, Albert Einstein College of Medicine (https://www.einstein.yu.edu/administration/animal-studies/), IACUC protocol # 2013–0405.

Ethics statement

Institutional Animal Care and Use Committee (IACUC) approved this research. The IACUC number is 20160507.

Whole mount in situ hybridization

Whole-mount RNA in situ hybridization with non-radioactive probes was performed as previously described [76, 77], using PCR-based probes, Tbx1 [78], Wnt2 forward primer: 5’ TGGCTCTGGCTCCCTCTGCT 3’ and reverse primer: 5’ CAGGGAGCCTGCCTCTCGGT 3’ and Wnt4 forward primer: 5’ CCGCGAGCAATTGGCTGTACC 3’ and reverse primer: 5’ TGGAACCTGCAGCCACAGCG 3’. Following whole-mount protocol, the embryos were fixed overnight in 4% paraformaldehyde (PFA) and then dehydrated through a graded ethanol series, embedded in paraffin and sectioned at 10 μm. Minimum of 5 embryos from 3 independent litters were analyzed per embryonic stage.

Proliferation and apoptosis on tissue sections

After fixation as described above, frozen sections were obtained at a thickness of 10 μm and then permeabilized in 0.5% Triton X-100 for 5 min. Blocking was performed with 5% serum (goat or donkey) in PBS/0.1% Triton X-100 (PBT) for 1 hour. Primary antibody was diluted in blocking solution (1:500) and incubated for 1 hour. Proliferation of cells was assessed by immunofluorescence using the primary antibody anti-phospho Histone H3 (Ser10), a mitosis marker (06–570 Millipore). Sections were washed in PBT and incubated with a secondary antibody for 1 hour. Secondary antibody was Alexa Fluor 568 goat a-rabbit IgG (A11011 Invitrogen) at 1:500. Slides were mounted in hard-set mounting medium with DAPI (Vector Labs H-1500). Images were captured using a Zeiss Axio Observer microscope. To perform statistical analysis of cell proliferation, we first counted the Mef2c-AHF-Cre, GFP positive cells in the pharyngeal apparatus located behind the heart in embryo sections and then calculated the average cell counts per tissue section for each embryo. Then we counted all proliferating cells in each section and calculated the ratio of proliferating cells within the Mef2c-AHF-Cre lineage. Then, we estimated the mean and standard error of the average cell counts for controls, Tbx1 LOF and rescued embryos and compared them using the t-test. Apoptosis was assessed on 10 μm thick frozen sections by using TMR Red In situ Cell Death kit (2156792 Roche) following the manufacturer’s instructions. Natural GFP from the reporter or an antibody for GFP (Abcam 6290) was used to distinguish the AHF cells in both assays described above. Representations of the complete AHF region from at least 4 embryos per genotype from at least 3 independent litters were used in each assay.

Wnt/β-catenin signaling in mouse wild type embryos

Wnt/β-catenin signaling reporter mice, TCF/Lef:H2B-GFP [75] were used to observe Wnt/β-catenin signaling by direct fluorescence of green fluorescent protein (GFP) in wild type embryos at embryonic day E9.5 (19–21 somite pairs). Mouse embryos were fixed and cryosectioned at 10 μm. Slides were mounted in hard-set mounting medium with DAPI to visualize DNA (Vector Labs H-1500). Images were then captured using a Zeiss Axio Observer microscope. Nuclear Wnt/β-catenin signaling was counted as the GFP positive signal that co-localized with the DNA. A minimum of 5 embryos from 3 independent litters was analyzed.

Mouse embryo heart histology and phenotypic analysis

Mouse embryos were isolated in phosphate-buffered saline (PBS) and fixed in 10% neutral buffered formalin (Sigma Corp.) overnight. Following fixation, the embryos were dehydrated through a graded ethanol series, embedded in paraffin and sectioned at 5 μm. All histological sections were stained with hematoxylin and eosin using standard protocols. Staining was performed in the Einstein Histopathology Core Facility (http://www.einstein.yu.edu/histopathology/page.aspx). For Tbx1 LOF mutants, a total of 70 hearts at E14.5 were obtained from more than 50 independent crosses and analyzed morphologically using light microscopy. For the rescue crosses, 56 hearts at E14.5 were obtained and the Fisher’s exact test was performed to compare the proportion of rescued phenotypes observed between rescued genotype hearts and the Tbx1 LOF mutants.

Direct fluorescence: Mef2c-AHF-Cre lineage tracing

Images were generated from GFP expressing embryos by direct fluorescence immediately following dissection. For tissue sections, embryos were fixed for 2 hours with embryos stage ≤ E10.5 (30–32 somite pairs). Fixation was carried out in 4% PFA in PBS at 4°C. After fixation, tissue was washed in PBS and then cryoprotected in 30% sucrose in PBS overnight at 4°C. Embryos were embedded in OCT and cryosectioned at 10 μm. Images were then captured using a Zeiss Axio Observer microscope.

Gene expression profiling on microarrays

Embryos at E9.5 (19–21 somites pairs) were used for global gene expression studies. To obtain enough RNA for microarray hybridization experiments, microdissected AHFs (defined here as: pharyngeal arches 2–6) from 27 of each of the following genotypes: Tbx1 LOF and its control (Tbx1f/+), Tbx1 GOF and its control (Tbx1-GFP/+), β-cat LOF and its control (β-cateninf/+), β-cat GOF and its control (β-catE3/+), rescue and its control (Tbx1f/+;β-cateninf/+), were pooled in groups of three or six according to the genotype. For this experiment we used controls that did not have Cre. Between 4–6 microarrays were performed per genotype in 2–3 batches. The tissue was homogenized in Buffer RLT (QIAGEN). Total RNA was isolated with the RNeasy Micro Kit according to the manufacturer’s protocol. Quality and quantity of total RNA were determined using an Agilent 2100 Bioanalyzer (Agilent) and an ND-1000 Spectrophotometer (NanoDrop), respectively. Biotinylated single-stranded cDNA targets were amplified from 100 nanograms (ng) starting total RNA using the Ovation RNA Amplification System V2 and FLOvation cDNA Biotin Module V2 (NuGEN). A total of 3.75 mg of cDNA was hybridized to the GeneChip Test3 array (Affymetrix) to test the quality of the labeled target. Nucleic acid samples that passed quality control were then hybridized to the Affymetrix Mouse GeneST 1.0 chip. Hybridization, washing, staining and scanning were performed in the Genomics Core at Einstein (https://www.einstein.yu.edu/research/shared-facilities/cores/46/genomics/) according to the Affymetrix manual.

Microarray data analysis

Data analysis was performed in the R statistical package. GeneChip data were pre-processed by the ‘oligo’ package [79], which implements Robust Multichip Average (RMA) algorithm with background correction, quantile normalization and gene level summarization [80]. Afterwards, for convenience of comparison, only probe-set assigned to genes were kept for subsequent analysis. Multiple probe-sets for the same genes were collapsed by “average” to obtain a single measurement per gene [81]. As some arrays were assayed in different batches, we performed UPGMA (unweighted pair group method with arithmetic mean) clustering of samples by transcriptomic profile similarities based on the Spearman correlation coefficients. This analysis indicated clear batch effects, especially for β-cat LOF and Tbx1 LOF data (data not shown). Hence, we applied ComBat, an efficient batch effect removal approach, to remove batch effects [82]. This analysis detected some individual arrays of poor quality that were then excluded. In the end, to keep a balance between controls and mutants, we analyzed 4 arrays per genotype. The ‘Limma’ package was used for determining differential expression [83]. To address the issue that adjustment of batch effect by any linear model based approach (including ComBat) can introduce a systematic correlation structure in the data, which may lead to exaggerated confidence in differential expression analysis [84], we accounted for this correlation in Limma by adding ‘blocking for batch’ in the model. In the end, genes with p-values < 0.05 were further explored. The microarray data has been deposited to the GEO database (accession number: GSE78125).

Quantitative RT-PCR

Embryos at E9.5 (19–21 somites pairs) were used for quantitative gene expression studies of microdissected AHFs from each of the following genotypes: Tbx1 LOF and its control (Tbx1f/+), Tbx1 GOF, β-cat LOF, β-cat GOF and rescue were pooled in groups according to genotype. Tbx1f/+ was used as control. To obtain enough total RNA and minimize the variability of gene expression in individual embryos, each biological replicate of RNA contained microdissected AHFs from six embryos of the same genotype at E9.5 collected from at least 3 independent litters. Three biological replicates were performed per genotype. The tissue was immediately frozen, samples were homogenized and total RNA was isolated with the RNeasy Micro Kit (Qiagen). Quality and quantity of total RNA was determined using an Agilent 2100 Bioanalyzer (Agilent) and a ND-1000 Spectrophotometer (NanoDrop), respectively. Single-stranded cDNA targets were amplified from 100 nanograms (ng) starting total RNA using the Ovation RNA Amplification System V2 and FL- Ovation cDNA Biotin Module V2 (NuGEN). The mRNA levels were measured using TaqMan Gene Expression assays (Applied Biosystems) for each gene and were carried out in triplicate using 18S (RNA, 18S ribosomal 1), Actb (Actin, beta) and B2m (Beta-2-microglobulin) genes as normalization controls. TaqMan probes and primer sets were obtained from the Applied Biosystems Gene Expression Assay database (http://allgenes.com). Samples were processed in standard 96-well plates (20 ul final volume per reaction and each reaction in triplicate containing 25 ng of cDNA) on an ABI 7900HT Q-PCR apparatus. The SDS 2.2 software platform (Applied Biosystems) was used for the computer interface with the ABI 7900HT PCR System to generate normalized data, compare samples, and calculate the relative quantity. Statistical significance of the difference in gene expression was estimated using ANOVA and the two-tailed t-test independently when type of comparison allowed it.

Supporting information

S1 Table. Embryo genotypes used as controls for Tbx1 LOF and rescue experiments (Figs 2 and 4).

The control embryo genotypes are listed in the first column. N is the number of embryos examined for heart and aortic arch anomalies. The % that is of a normal phenotype is indicated. Those with defects are identical to what was previously published as indicated.

https://doi.org/10.1371/journal.pgen.1006687.s001

(PDF)

S1 Fig. Tbx1 inactivation in Tbx1 LOF embryos.

Whole mount in situ hybridization (WMISH) of Tbx1 in a control embryo at embryonic day 9.5 (E9.5) and the corresponding sagittal section is shown in B. (C) WMISH of Tbx1 in a AHF conditional mutant and the corresponding sagittal section is shown in D. 1 and 2 indicate the first and second pharyngeal arches, respectively. H&E histological sections of embryos at E12.5, control (E and F) and Tbx1 LOF (G and H). Abbreviations: heart (H), aorta (Ao), pulmonary trunk (PT) and persistent truncus arteriosus (PTA).

https://doi.org/10.1371/journal.pgen.1006687.s002

(PDF)

S2 Fig. Proliferation of AHF in control, Tbx1 LOF and rescue embryos.

(A) Immunofluorescence images of sagittal sections to visualize the AHF lineage (GFP, green) and cell proliferation (anti-phospho Histone H3 (Ser10); red); in control, Tbx1 LOF and rescue embryos are shown. DAPI fluorescent stain to visualize nuclei and identify the tissue is shown in blue. Statistical analysis was performed to determine whether cell proliferation was the same or different between groups of embryos by two-tailed t-test, p value <0.05. Error bars = standard deviation (SD). Abbreviations: outflow tract (OFT), right ventricle (RV).

https://doi.org/10.1371/journal.pgen.1006687.s003

(PDF)

S3 Fig. Apoptosis of AHF in control, Tbx1 LOF and rescue embryos.

(A) Immunofluorescence images of sagittal sections to visualize the AHF lineage (GFP, green) and apoptosis (TUNEL, red); in control, Tbx1 LOF and rescue embryos are shown. DAPI fluorescent stain to visualize nuclei and identify the tissue is shown in blue. Statistical analysis was performed to determine whether the number of dead cells was the same or different between groups of embryos by two-tailed t-test, p value <0.05. Error bars = standard deviation (SD. Abbreviations: outflow tract (OFT), right ventricle (RV).

https://doi.org/10.1371/journal.pgen.1006687.s004

(PDF)

S4 Fig. Detection of Tbx1 and β-catenin by qRT-PCR in control, Tbx1 LOF and β-catenin GOF embryos.

Statistical significance of the difference in gene expression was estimated using two-tailed t-test, FC = fold change, p values < 0.05. Note that β-catenin mRNA is not expected to have a significant change in the AHF of β-catenin GOF embryos due to be a constitutive protein activation.

https://doi.org/10.1371/journal.pgen.1006687.s005

(PDF)

S5 Fig. Posterior histology analysis of representative embryos with the “rescue” genotype.

Additional transverse H&E histological sections of rescued hearts (Tbx1 LOF with loss of one allele of β-catenin in the Mef2c-AHF-Cre domain) at E14.5. (A—F) and (G–L) sections show two embryonic hearts with rescued septation between the two ventricles throughout the heart. Abbreviations: left atrium (LA), right atrium (RA), left ventricle (LV), right ventricle (RV).

https://doi.org/10.1371/journal.pgen.1006687.s006

(PDF)

S6 Fig. Full heatmap of gene expression changes in Tbx1 LOF, β-catenin and rescue genotypes (Fig 6B).

Heatmap showing the expression changes for all genes differentially expressed (p < 0.01) between Tbx1 LOF, β-catenin LOF or rescue embryos vs their respective controls.

https://doi.org/10.1371/journal.pgen.1006687.s007

(PDF)

S7 Fig. Gene expression of cardiac morphogenesis genes.

(A) Quantitative PCR was performed on micro-dissected AHF of Tbx1 LOF, β-catenin LOF and rescue embryos to detect the expression levels of selected genes known important for cardiac morphogenesis. Tbx1f/+ was used as control for expression plates due to number of genes tested and array design. Statistical significance of the difference in gene expression was estimated using ANOVA; p values < 0.05. Asterisks note those genes which expressions were, at least, significantly different between control and rescue embryos. FC = fold change. Error bars = standard deviation (SD).

https://doi.org/10.1371/journal.pgen.1006687.s008

(PDF)

Acknowledgments

We thank Dr Deborah Yelon for cmlc2 and nkx2.5 probes. We thank Dr Tao Wang for statistical analysis. We thank the Genomics Core at Einstein for performing Affymetrix gene expression arrays. We also want to thank the Histopathology Facility at Einstein.

Author Contributions

  1. Conceptualization: SER EMO CLC BEM.
  2. Data curation: ML TG DZ SER.
  3. Formal analysis: TG DZ ML.
  4. Funding acquisition: CLC DZ BEM.
  5. Investigation: SER EH GSD.
  6. Methodology: SER DZ.
  7. Project administration: BEM.
  8. Resources: BEM.
  9. Supervision: SER EMO CLC DZ BEM.
  10. Validation: SER GSD EH DZ BEM.
  11. Visualization: SER EH ML GSD EMO CLC DZ BEM.
  12. Writing – original draft: SER BEM.
  13. Writing – review & editing: EH TG EMO CLC DZ BEM.

References

  1. 1. Gross SJ, Stosic M, McDonald-McGinn DM, Bassett AS, Norvez A, Dhamankar R, et al. Clinical Experience with Single-Nucleotide Polymorphism-Based Noninvasive Prenatal Screening for 22q11.2 Deletion Syndrome. Ultrasound Obstet Gynecol. 2015.
  2. 2. Burn J, Goodship J. Developmental genetics of the heart. Curr Opin Genet Dev. 1996;6(3):322–5. pmid:8791511
  3. 3. Goldmuntz E, Driscoll D, Budarf ML, Zackai EH, McDonald-McGinn DM, Biegel JA, et al. Microdeletions of chromosomal region 22q11 in patients with congenital conotruncal cardiac defects. J Med Genet. 1993;30(10):807–12. PubMed Central PMCID: PMCPMC1016559. pmid:7901419
  4. 4. Chapman DL, Garvey N, Hancock S, Alexiou M, Agulnik SI, Gibson-Brown JJ, et al. Expression of the T-box family genes, Tbx1-Tbx5, during early mouse development. Dev Dyn. 1996;206(4):379–90. pmid:8853987
  5. 5. Jerome LA, Papaioannou VE. DiGeorge syndrome phenotype in mice mutant for the T-box gene, Tbx1. Nat Genet. 2001;27(3):286–91. pmid:11242110
  6. 6. Lindsay EA, Vitelli F, Su H, Morishima M, Huynh T, Pramparo T, et al. Tbx1 haploinsufficieny in the DiGeorge syndrome region causes aortic arch defects in mice. Nature. 2001;410(6824):97–101. pmid:11242049
  7. 7. Merscher S, Funke B, Epstein JA, Heyer J, Puech A, Lu MM, et al. TBX1 is responsible for cardiovascular defects in velo-cardio-facial/DiGeorge syndrome. Cell. 2001;104(4):619–29. pmid:11239417
  8. 8. Yagi H, Furutani Y, Hamada H, Sasaki T, Asakawa S, Minoshima S, et al. Role of TBX1 in human del22q11.2 syndrome. Lancet. 2003;362(9393):1366–73. pmid:14585638
  9. 9. Stoller JZ, Epstein JA. Identification of a novel nuclear localization signal in Tbx1 that is deleted in DiGeorge syndrome patients harboring the 1223delC mutation. Hum Mol Genet. 2005;14(7):885–92. pmid:15703190
  10. 10. Saga Y, Miyagawa-Tomita S, Takagi A, Kitajima S, Miyazaki J, Inoue T. MesP1 is expressed in the heart precursor cells and required for the formation of a single heart tube. Development. 1999;126(15):3437–47. pmid:10393122
  11. 11. Stalsberg H, DeHaan RL. The precardiac areas and formation of the tubular heart in the chick embryo. Dev Biol. 1969;19(2):128–59. pmid:4180350
  12. 12. Mjaatvedt CH, Nakaoka T, Moreno-Rodriguez R, Norris RA, Kern MJ, Eisenberg CA, et al. The outflow tract of the heart is recruited from a novel heart-forming field. Dev Biol. 2001;238(1):97–109. pmid:11783996
  13. 13. Cai CL, Liang X, Shi Y, Chu PH, Pfaff SL, Chen J, et al. Isl1 identifies a cardiac progenitor population that proliferates prior to differentiation and contributes a majority of cells to the heart. Dev Cell. 2003;5(6):877–89. pmid:14667410
  14. 14. Kelly RG, Brown NA, Buckingham ME. The arterial pole of the mouse heart forms from Fgf10-expressing cells in pharyngeal mesoderm. Dev Cell. 2001;1(3):435–40. pmid:11702954
  15. 15. Buckingham M, Meilhac S, Zaffran S. Building the mammalian heart from two sources of myocardial cells. Nat Rev Genet. 2005;6(11):826–35. pmid:16304598
  16. 16. Waldo KL, Kumiski DH, Wallis KT, Stadt HA, Hutson MR, Platt DH, et al. Conotruncal myocardium arises from a secondary heart field. Development. 2001;128(16):3179–88. pmid:11688566
  17. 17. Zaffran S, Kelly RG, Meilhac SM, Buckingham ME, Brown NA. Right ventricular myocardium derives from the anterior heart field. Circ Res. 2004;95(3):261–8. pmid:15217909
  18. 18. Verzi MP, McCulley DJ, De Val S, Dodou E, Black BL. The right ventricle, outflow tract, and ventricular septum comprise a restricted expression domain within the secondary/anterior heart field. Dev Biol. 2005;287(1):134–45. pmid:16188249
  19. 19. Bajolle F, Zaffran S, Meilhac SM, Dandonneau M, Chang T, Kelly RG, et al. Myocardium at the base of the aorta and pulmonary trunk is prefigured in the outflow tract of the heart and in subdomains of the second heart field. Dev Biol. 2008;313(1):25–34. pmid:18005956
  20. 20. Dominguez JN, Meilhac SM, Bland YS, Buckingham ME, Brown NA. Asymmetric fate of the posterior part of the second heart field results in unexpected left/right contributions to both poles of the heart. Circ Res. 2012;111(10):1323–35. pmid:22955731
  21. 21. Vincent SD, Buckingham ME. How to make a heart: the origin and regulation of cardiac progenitor cells. Curr Top Dev Biol. 2010;90:1–41. pmid:20691846
  22. 22. Chen L, Fulcoli FG, Tang S, Baldini A. Tbx1 regulates proliferation and differentiation of multipotent heart progenitors. Circ Res. 2009;105(9):842–51. PubMed Central PMCID: PMCPMC2796444. pmid:19745164
  23. 23. Watanabe Y, Zaffran S, Kuroiwa A, Higuchi H, Ogura T, Harvey RP, et al. Fibroblast growth factor 10 gene regulation in the second heart field by Tbx1, Nkx2-5, and Islet1 reveals a genetic switch for down-regulation in the myocardium. Proc Natl Acad Sci U S A. 2012;109(45):18273–80. PubMed Central PMCID: PMCPMC3494960. pmid:23093675
  24. 24. Liao J, Aggarwal VS, Nowotschin S, Bondarev A, Lipner S, Morrow BE. Identification of downstream genetic pathways of Tbx1 in the second heart field. Dev Biol. 2008;316(2):524–37. PubMed Central PMCID: PMCPMC2494702. pmid:18328475
  25. 25. Nevis K, Obregon P, Walsh C, Guner-Ataman B, Burns CG, Burns CE. Tbx1 is required for second heart field proliferation in zebrafish. Dev Dyn. 2013;242(5):550–9. PubMed Central PMCID: PMCPMC3676967. pmid:23335360
  26. 26. Gessert S, Kuhl M. The multiple phases and faces of wnt signaling during cardiac differentiation and development. Circ Res. 2010;107(2):186–99. pmid:20651295
  27. 27. Ai D, Fu X, Wang J, Lu MF, Chen L, Baldini A, et al. Canonical Wnt signaling functions in second heart field to promote right ventricular growth. Proc Natl Acad Sci U S A. 2007;104(22):9319–24. PubMed Central PMCID: PMCPMC1890492. pmid:17519332
  28. 28. Huh SH, Ornitz DM. Beta-catenin deficiency causes DiGeorge syndrome-like phenotypes through regulation of Tbx1. Development. 2010;137(7):1137–47. PubMed Central PMCID: PMCPMC2835329. pmid:20215350
  29. 29. Brault V, Moore R, Kutsch S, Ishibashi M, Rowitch DH, McMahon AP, et al. Inactivation of the beta-catenin gene by Wnt1-Cre-mediated deletion results in dramatic brain malformation and failure of craniofacial development. Development. 2001;128(8):1253–64. pmid:11262227
  30. 30. Harada N, Tamai Y, Ishikawa T, Sauer B, Takaku K, Oshima M, et al. Intestinal polyposis in mice with a dominant stable mutation of the beta-catenin gene. EMBO J. 1999;18(21):5931–42. PubMed Central PMCID: PMCPMC1171659. pmid:10545105
  31. 31. Ivins S, Lammerts van Beuren K, Roberts C, James C, Lindsay E, Baldini A, et al. Microarray analysis detects differentially expressed genes in the pharyngeal region of mice lacking Tbx1. Dev Biol. 2005;285(2):554–69. pmid:16109395
  32. 32. Zhang Z, Cerrato F, Xu H, Vitelli F, Morishima M, Vincentz J, et al. Tbx1 expression in pharyngeal epithelia is necessary for pharyngeal arch artery development. Development. 2005;132(23):5307–15. pmid:16284121
  33. 33. Xu H, Morishima M, Wylie JN, Schwartz RJ, Bruneau BG, Lindsay EA, et al. Tbx1 has a dual role in the morphogenesis of the cardiac outflow tract. Development. 2004;131(13):3217–27. pmid:15175244
  34. 34. Freyer L, Nowotschin S, Pirity MK, Baldini A, Morrow BE. Conditional and constitutive expression of a Tbx1-GFP fusion protein in mice. BMC Dev Biol. 2013;13:33. PubMed Central PMCID: PMCPMC3765320. pmid:23971992
  35. 35. Hiroi Y, Kudoh S, Monzen K, Ikeda Y, Yazaki Y, Nagai R, et al. Tbx5 associates with Nkx2-5 and synergistically promotes cardiomyocyte differentiation. Nat Genet. 2001;28(3):276–80. pmid:11431700
  36. 36. Maitra M, Schluterman MK, Nichols HA, Richardson JA, Lo CW, Srivastava D, et al. Interaction of Gata4 and Gata6 with Tbx5 is critical for normal cardiac development. Dev Biol. 2009;326(2):368–77. PubMed Central PMCID: PMCPMC2651674. pmid:19084512
  37. 37. Wang C, Cao D, Wang Q, Wang DZ. Synergistic activation of cardiac genes by myocardin and Tbx5. PLoS One. 2011;6(8):e24242. PubMed Central PMCID: PMCPMC3163680. pmid:21897873
  38. 38. Klaus A, Saga Y, Taketo MM, Tzahor E, Birchmeier W. Distinct roles of Wnt/beta-catenin and Bmp signaling during early cardiogenesis. Proc Natl Acad Sci U S A. 2007;104(47):18531–6. PubMed Central PMCID: PMCPMC2141811. pmid:18000065
  39. 39. Kwon C, Arnold J, Hsiao EC, Taketo MM, Conklin BR, Srivastava D. Canonical Wnt signaling is a positive regulator of mammalian cardiac progenitors. Proc Natl Acad Sci U S A. 2007;104(26):10894–9. PubMed Central PMCID: PMCPMC1904134. pmid:17576928
  40. 40. Lin L, Cui L, Zhou W, Dufort D, Zhang X, Cai CL, et al. Beta-catenin directly regulates Islet1 expression in cardiovascular progenitors and is required for multiple aspects of cardiogenesis. Proc Natl Acad Sci U S A. 2007;104(22):9313–8. PubMed Central PMCID: PMCPMC1890491. pmid:17519333
  41. 41. Kwon C, Qian L, Cheng P, Nigam V, Arnold J, Srivastava D. A regulatory pathway involving Notch1/beta-catenin/Isl1 determines cardiac progenitor cell fate. Nat Cell Biol. 2009;11(8):951–7. PubMed Central PMCID: PMCPMC2748816. pmid:19620969
  42. 42. Marvin MJ, Di Rocco G, Gardiner A, Bush SM, Lassar AB. Inhibition of Wnt activity induces heart formation from posterior mesoderm. Genes Dev. 2001;15(3):316–27. PubMed Central PMCID: PMCPMC312622. pmid:11159912
  43. 43. Schneider VA, Mercola M. Wnt antagonism initiates cardiogenesis in Xenopus laevis. Genes Dev. 2001;15(3):304–15. PubMed Central PMCID: PMCPMC312618. pmid:11159911
  44. 44. Tzahor E, Lassar AB. Wnt signals from the neural tube block ectopic cardiogenesis. Genes Dev. 2001;15(3):255–60. PubMed Central PMCID: PMCPMC312627. pmid:11159906
  45. 45. Ueno S, Weidinger G, Osugi T, Kohn AD, Golob JL, Pabon L, et al. Biphasic role for Wnt/beta-catenin signaling in cardiac specification in zebrafish and embryonic stem cells. Proc Natl Acad Sci U S A. 2007;104(23):9685–90. PubMed Central PMCID: PMCPMC1876428. pmid:17522258
  46. 46. Mishina Y, Hanks MC, Miura S, Tallquist MD, Behringer RR. Generation of Bmpr/Alk3 conditional knockout mice. Genesis. 2002;32(2):69–72. pmid:11857780
  47. 47. Bondue A, Blanpain C. Mesp1: a key regulator of cardiovascular lineage commitment. Circ Res. 2010;107(12):1414–27. pmid:21148448
  48. 48. Aggarwal VS, Liao J, Bondarev A, Schimmang T, Lewandoski M, Locker J, et al. Dissection of Tbx1 and Fgf interactions in mouse models of 22q11DS suggests functional redundancy. Hum Mol Genet. 2006;15(21):3219–28. pmid:17000704
  49. 49. Theveniau-Ruissy M, Perez-Pomares JM, Parisot P, Baldini A, Miquerol L, Kelly RG. Coronary stem development in wildtype and Tbx1 null mouse hearts. Dev Dyn. 2015.
  50. 50. Sinha T, Li D, Theveniau-Ruissy M, Hutson MR, Kelly RG, Wang J. Loss of Wnt5a disrupts second heart field cell deployment and may contribute to OFT malformations in DiGeorge syndrome. Hum Mol Genet. 2015;24(6):1704–16. PubMed Central PMCID: PMCPMC4381755. pmid:25410658
  51. 51. Okubo T, Kawamura A, Takahashi J, Yagi H, Morishima M, Matsuoka R, et al. Ripply3, a Tbx1 repressor, is required for development of the pharyngeal apparatus and its derivatives in mice. Development. 2011;138(2):339–48. pmid:21177346
  52. 52. Garg V, Kathiriya IS, Barnes R, Schluterman MK, King IN, Butler CA, et al. GATA4 mutations cause human congenital heart defects and reveal an interaction with TBX5. Nature. 2003;424(6947):443–7. pmid:12845333
  53. 53. Lickert H, Takeuchi JK, Von Both I, Walls JR, McAuliffe F, Adamson SL, et al. Baf60c is essential for function of BAF chromatin remodelling complexes in heart development. Nature. 2004;432(7013):107–12. pmid:15525990
  54. 54. Takeuchi JK, Bruneau BG. Directed transdifferentiation of mouse mesoderm to heart tissue by defined factors. Nature. 2009;459(7247):708–11. PubMed Central PMCID: PMCPMC2728356. pmid:19396158
  55. 55. Li Q, Guo ZK, Chang YQ, Yu X, Li CX, Li H. Gata4, Tbx5 and Baf60c induce differentiation of adipose tissue-derived mesenchymal stem cells into beating cardiomyocytes. Int J Biochem Cell Biol. 2015;66:30–6. pmid:26071180
  56. 56. Grepin C, Robitaille L, Antakly T, Nemer M. Inhibition of transcription factor GATA-4 expression blocks in vitro cardiac muscle differentiation. Mol Cell Biol. 1995;15(8):4095–102. PubMed Central PMCID: PMCPMC230648. pmid:7623805
  57. 57. Ip HS, Wilson DB, Heikinheimo M, Tang Z, Ting CN, Simon MC, et al. The GATA-4 transcription factor transactivates the cardiac muscle-specific troponin C promoter-enhancer in nonmuscle cells. Mol Cell Biol. 1994;14(11):7517–26. PubMed Central PMCID: PMCPMC359288. pmid:7935467
  58. 58. Molkentin JD, Kalvakolanu DV, Markham BE. Transcription factor GATA-4 regulates cardiac muscle-specific expression of the alpha-myosin heavy-chain gene. Mol Cell Biol. 1994;14(7):4947–57. PubMed Central PMCID: PMCPMC358867. pmid:8007990
  59. 59. Jiang Y, Evans T. The Xenopus GATA-4/5/6 genes are associated with cardiac specification and can regulate cardiac-specific transcription during embryogenesis. Dev Biol. 1996;174(2):258–70. pmid:8631498
  60. 60. Opdecamp K, Nakayama A, Nguyen MT, Hodgkinson CA, Pavan WJ, Arnheiter H. Melanocyte development in vivo and in neural crest cell cultures: crucial dependence on the Mitf basic-helix-loop-helix-zipper transcription factor. Development. 1997;124(12):2377–86. pmid:9199364
  61. 61. Zeisberg EM, Ma Q, Juraszek AL, Moses K, Schwartz RJ, Izumo S, et al. Morphogenesis of the right ventricle requires myocardial expression of Gata4. J Clin Invest. 2005;115(6):1522–31. PubMed Central PMCID: PMCPMC1090473. pmid:15902305
  62. 62. Li QY, Newbury-Ecob RA, Terrett JA, Wilson DI, Curtis AR, Yi CH, et al. Holt-Oram syndrome is caused by mutations in TBX5, a member of the Brachyury (T) gene family. Nat Genet. 1997;15(1):21–9. pmid:8988164
  63. 63. Basson CT, Bachinsky DR, Lin RC, Levi T, Elkins JA, Soults J, et al. Mutations in human TBX5 [corrected] cause limb and cardiac malformation in Holt-Oram syndrome. Nat Genet. 1997;15(1):30–5. pmid:8988165
  64. 64. Liberatore CM, Searcy-Schrick RD, Yutzey KE. Ventricular expression of tbx5 inhibits normal heart chamber development. Dev Biol. 2000;223(1):169–80. pmid:10864469
  65. 65. Garrity DM, Childs S, Fishman MC. The heartstrings mutation in zebrafish causes heart/fin Tbx5 deficiency syndrome. Development. 2002;129(19):4635–45. pmid:12223419
  66. 66. Christoforou N, Chellappan M, Adler AF, Kirkton RD, Wu T, Addis RC, et al. Transcription factors MYOCD, SRF, Mesp1 and SMARCD3 enhance the cardio-inducing effect of GATA4, TBX5, and MEF2C during direct cellular reprogramming. PLoS One. 2013;8(5):e63577. PubMed Central PMCID: PMCPMC3660533. pmid:23704920
  67. 67. Turbendian HK, Gordillo M, Tsai SY, Lu J, Kang G, Liu TC, et al. GATA factors efficiently direct cardiac fate from embryonic stem cells. Development. 2013;140(8):1639–44. PubMed Central PMCID: PMCPMC3621482. pmid:23487308
  68. 68. Alexandrovich A, Arno M, Patient RK, Shah AM, Pizzey JA, Brewer AC. Wnt2 is a direct downstream target of GATA6 during early cardiogenesis. Mech Dev. 2006;123(4):297–311. pmid:16621466
  69. 69. Tian Y, Yuan L, Goss AM, Wang T, Yang J, Lepore JJ, et al. Characterization and in vivo pharmacological rescue of a Wnt2-Gata6 pathway required for cardiac inflow tract development. Dev Cell. 2010;18(2):275–87. PubMed Central PMCID: PMCPMC2846539. pmid:20159597
  70. 70. Edwards JJ, Gelb BD. Genetics of congenital heart disease. Current opinion in cardiology. 2016;31(3):235–41. PubMed Central PMCID: PMC4868504. pmid:26872209
  71. 71. Azhar M, Ware SM. Genetic and Developmental Basis of Cardiovascular Malformations. Clinics in perinatology. 2016;43(1):39–53. PubMed Central PMCID: PMC5215107. pmid:26876120
  72. 72. Guo T, Chung JH, Wang T, McDonald-McGinn DM, Kates WR, Hawula W, et al. Histone Modifier Genes Alter Conotruncal Heart Phenotypes in 22q11.2 Deletion Syndrome. Am J Hum Genet. 2015;97(6):869–77. PubMed Central PMCID: PMCPMC4678435. pmid:26608785
  73. 73. Arnold JS, Werling U, Braunstein EM, Liao J, Nowotschin S, Edelmann W, et al. Inactivation of Tbx1 in the pharyngeal endoderm results in 22q11DS malformations. Development. 2006;133(5):977–87. pmid:16452092
  74. 74. Sousa VH, Miyoshi G, Hjerling-Leffler J, Karayannis T, Fishell G. Characterization of Nkx6-2-derived neocortical interneuron lineages. Cereb Cortex. 2009;19 Suppl 1:i1–10. PubMed Central PMCID: PMCPMC2693535.
  75. 75. Ferrer-Vaquer A, Piliszek A, Tian G, Aho RJ, Dufort D, Hadjantonakis AK. A sensitive and bright single-cell resolution live imaging reporter of Wnt/ss-catenin signaling in the mouse. BMC Dev Biol. 2010;10:121. PubMed Central PMCID: PMCPMC3017038. pmid:21176145
  76. 76. Hidai H, Bardales R, Goodwin R, Quertermous T, Quertermous EE. Cloning of capsulin, a basic helix-loop-helix factor expressed in progenitor cells of the pericardium and the coronary arteries. Mech Dev. 1998;73(1):33–43. pmid:9545526
  77. 77. Alappat SR, Zhang Z, Suzuki K, Zhang X, Liu H, Jiang R, et al. The cellular and molecular etiology of the cleft secondary palate in Fgf10 mutant mice. Dev Biol. 2005;277(1):102–13. pmid:15572143
  78. 78. Funke B, Epstein JA, Kochilas LK, Lu MM, Pandita RK, Liao J, et al. Mice overexpressing genes from the 22q11 region deleted in velo-cardio-facial syndrome/DiGeorge syndrome have middle and inner ear defects. Hum Mol Genet. 2001;10(22):2549–56. pmid:11709542
  79. 79. Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010;26(19):2363–7. PubMed Central PMCID: PMCPMC2944196. pmid:20688976
  80. 80. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4(2):249–64. pmid:12925520
  81. 81. Miller JA, Cai C, Langfelder P, Geschwind DH, Kurian SM, Salomon DR, et al. Strategies for aggregating gene expression data: the collapseRows R function. BMC Bioinformatics. 2011;12:322. PubMed Central PMCID: PMCPMC3166942. pmid:21816037
  82. 82. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27. pmid:16632515
  83. 83. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47. PubMed Central PMCID: PMCPMC4402510. pmid:25605792
  84. 84. Nygaard V, Rodland EA, Hovig E. Methods that remove batch effects while retaining group differences may lead to exaggerated confidence in downstream analyses. Biostatistics. 2016;17(1):29–39. PubMed Central PMCID: PMCPMC4679072. pmid:26272994