Skip to main content
Download PDF

A deficiency in SUMOylation activity disrupts multiple pathways leading to neural tube and heart defects in Xenopus embryos

BMC Genomics volume 20, Article number: 386 (2019) Cite this article

Abstract

Background

Adenovirus protein, Gam1, triggers the proteolytic destruction of the E1 SUMO-activating enzyme. Microinjection of an empirically determined amount of Gam1 mRNA into one-cell Xenopus embryos can reduce SUMOylation activity to undetectable, but nonlethal, levels, enabling an examination of the role of this post-translational modification during early vertebrate development.

Results

We find that SUMOylation-deficient embryos consistently exhibit defects in neural tube and heart development. We have measured differences in gene expression between control and embryos injected with Gam1 mRNA at three developmental stages: early gastrula (immediately following the initiation of zygotic transcription), late gastrula (completion of the formation of the three primary germ layers), and early neurula (appearance of the neural plate). Although changes in gene expression are widespread and can be linked to many biological processes, three pathways, non-canonical Wnt/PCP, snail/twist, and Ets-1, are especially sensitive to the loss of SUMOylation activity and can largely account for the predominant phenotypes of Gam1 embryos. SUMOylation appears to generate different pools of a given transcription factor having different specificities with this post-translational modification involved in the regulation of more complex, as opposed to housekeeping, processes.

Conclusions

We have identified changes in gene expression that underlie the neural tube and heart phenotypes resulting from depressed SUMOylation activity. Notably, these developmental defects correspond to the two most frequently occurring congenital birth defects in humans, strongly suggesting that perturbation of SUMOylation, either globally or of a specific protein, may frequently be the origin of these pathologies.

Background

Post-translational modification by small ubiquitin-like modifier (SUMO) has emerged as a global mechanism for the regulation of protein activity, stability, and localization [1,2,3,4,5]. Its integration into other types of protein modification, such as phosphorylation, acetylation, and ubiquination, allows for coordinated control of diverse biological processes. Although the SUMO (E2) conjugating enzyme is capable of acting directly on target proteins, in the majority of cases a sizable family of SUMO (E3) ligases determine substrate specificity. The SENP family of proteases, likewise exhibiting different specificities, acts to remove SUMO moieties so that the modification is reversible and can be tightly regulated.

The steady-state level of SUMOylated protein in any given instance is typically low and has been referred to as the “SUMO enigma” [6]. This has made detection of SUMO-modified proteins difficult. Nonetheless, a strategy that employed the expression of epitope-tagged SUMO2 in human cells combined with flow cytometry enabled the identification of more than 1600 proteins conjugated to this SUMO isoform [4]. Similar to earlier studies with yeast cells [5, 7], SUMO modification in this case was mostly limited to nuclear proteins or those that are nucleocytoplasmic. This is consistent with the large number of transcription factors, RNA binding proteins, chromatin associated proteins, and cell cycle regulators that are documented targets of SUMOylation [8, 9]. In Xenopus egg extract, a large number of SUMOylated proteins are associated with chromatin [10]. However, a bioinformatic analysis of this proteomic data revealed that an appreciable number (39.5%) of the identified SUMOylated proteins were linked to metabolic processes and translation, indicating the importance of this post-translational modification in cytoplasmic processes as well.

The importance and scope of SUMOylation during early development has been difficult to ascertain, since complete elimination of this activity is embryonic lethal [2, 8, 11,12,13]. From yeast to vertebrates, the essential role of SUMOylation in mitotic processes has defeated approaches that have knocked out this pathway by elimination of the sole E2 conjugating enzyme. Strategies that have eliminated only one of the SUMO isoforms have resulted in less severe phenotypes in the case of SUMO1 [14,15,16], but embryonic lethal phenotypes in the case of SUMO2 [17]. Gene knockout or other methods of gene inactivation that target other components of the SUMOylation machinery (i.e., E3 ligases and SENP proteases) result in a variety of phenotypes [11].

A means to diminish, but not eliminate, SUMOylation activity would enable an assessment of the pathway’s importance and ubiquity during early development. The avian adenovirus protein, Gam1, binds to the SAE1 subunit of the E1 SUMO-activating enzyme and triggers its proteolytic degradation [18, 19], and injection of mRNA encoding Gam1 into Xenopus embryos reduces SUMOylation activity to undetectable levels [20]. Here, we have used Gam1 in order to measure changes in gene expression in SUMOylation-deficient embryos at three developmental stages, early gastrula, late gastrula, and early neurula using microarray technology, which can analyze the activity of approximately 30,000 genes. The data show that SUMOylation impacts many disparate biological processes; nonetheless, certain signaling pathways appear to be particularly sensitive to depletion of this activity and can be correlated with the predominant phenotypes of these embryos that include failure of blastopore and neural tube closure, shortened anterior-posterior (A-P) axis, and defective heart and eye development. Importantly, we show that disruptions in these pathways due to loss of this post-translational modification during the earliest periods of embryogenesis can elicit developmental deficits that correspond to the most frequently occurring human birth defects.

Results and discussion

Gam1 depletion of SUMOylation activity in Xenopus embryos

The adenovirus protein Gam1 triggers the proteolytic destruction of the SUMO E1 activating enzyme [18, 19], providing a strategic alternative to antisense methods, since the effect is immediate and eliminates both existing and de novo accumulation of the activating enzyme. We have reported that injection of mRNA encoding Gam1 into one-cell Xenopus embryos suppresses SUMOylation activity to undetectable levels in extract prepared from early neurula stage embryos [20]. In order to characterize the activity of Gam1 further, we measured the persistence of the protein during development. Fertilized eggs were injected with mRNA (0.5 ng) encoding myc-tagged Gam1; protein extract prepared from embryos at the designated developmental stage was analyzed by western blot. Gam1 protein is present by midblastula stage and reaches a maximum during early gastrula before declining during neurulation (Fig. 1a).

Fig. 1
figure 1

Depletion of SUMOylation activity in embryos injected with Gam1 mRNA. a Expression levels of Gam1 protein during embryogenesis. One-cell embryos were injected with mRNA (0.5 ng) encoding Gam1 with an N-terminal myc tag. Whole cell protein extract was prepared from embryos at the indicated Nieuwkoop–Faber stage and 25 μg taken for western blot analysis using anti-myc antibody. Lanes: 1, midblastula; 2, late blastula; 3, early gastrula; 4, early neurula; 5, late neurula; 6, water injected control. b In vitro SUMOylation assay. All assays contained Ubc9 (E2 enzyme), SUMO1, ATP, and substrate peptide (25 kDa) with (lane 1) or without (lane 2) purified E1 enzyme (500 nM) added. Samples were analyzed by western blot using an antibody specific for a 25 kDa SUMO substrate peptide. c Assays for E1 activity. One-cell embryos were injected with the indicated amount of Gam1 mRNA and allowed to develop to the indicated stage. Whole cell extract was prepared from 20 embryos and an equivalent amount of protein (25 μg) was used as a source of E1 enzyme for each assay

Full size image

An in vitro assay containing E1, E2 (UBC9), SUMO1, and a 25 kDa substrate peptide was used to measure SUMOylation activity by western blot (Fig. 1b). The effect of different amounts of injected Gam1 mRNA on SUMOylation activity was determined at four stages of development using whole cell extract prepared from water or Gam1 injected embryos as a source of E1 enzyme in the in vitro assay (Fig. 1c). All three amounts of injected Gam1 mRNA effectively suppress SUMOylation activity at midblastula. However, at late blastula there is measureable SUMOylation activity in all samples, presumably due to the approximate 2-fold increase during mid- and late blastula in the mRNAs that encode the two subunits (SAE and UBA2) that comprise the E1 enzyme [21]. Suppression of SUMOylation activity is reestablished in early gastrula embryos (0.5 or 1 ng injected mRNA) that then persist until late neurula. The excellent correlation between the levels of Gam1 protein and suppression of E1 activity is in accord with other evidence that expression of the viral protein is an effective method to control the SUMO pathway [18, 19]. The transient appearance of SUMOylation activity at the midblastula transition (MBT) possibly accounts for the survival of these embryos. Injection of higher amounts (5 ng) of Gam1 mRNA is lethal with no survival of injected embryos beyond the gastrula-neurula transition (stage 13). A mutant form of Gam1 that cannot bind to the SAE1 subunit had no effect on viability when injected at the same amount. Based on these results, 0.5 ng of injected Gam1 mRNA was chosen for the subsequent experiments to examine the role of SUMOylation in developing embryos.

Distinct phenotypes in SUMOylation deficient Xenopus embryos

Embryos injected with Gam1 mRNA develop normally up through late blastula compared to water injected or uninjected control embryos. One of the earliest phenotypes of Gam1 embryos is delayed or incomplete closure of the blastopore, which is a hallmark of gastrulation (Fig. 2a). In instances where closure is simply slowed, the embryo usually develops normally after this point. Embryos, in which the blastopore fails to close, do not develop normally and exhibit cascading effects that are apparent at the neurula stage of development. During neurulation, the neural tube forms due to convergent extension of the outermost layer of cells on the embryo [22, 23]. The posterior portion of the neural tube closes at or very near to the point of blastopore closure. Therefore, when the blastopore fails to close properly, the neural tube also fails to close properly and embryos display a spina bifida like phenotype at the late neurula stage. In the most extreme cases, embryos injected with Gam1 present completely open neural tubes with failure of the neural folds to fuse together along the entire A–P axis (Fig. 2b).

Fig. 2
figure 2

Developmental defects of Gam1-injected embryos. a Control (left) and Gam1-injected (right) late gastrula stage embryos showing delayed closure of the blastopore (arrow) in the latter. b Stage 35 (50 hpf) control embryo (left), Gam1 embryo with shortened A–P axis (center), Gam1 embryo with an open neural tube and edema around the heart (right). c An example of an eye defect (cyclopia). d Failure of the heart tube to initiate looping. (Left) A control embryo (96 hpf) in which the two atria (red) and ventricle (black) have fully developed compared to (right) a Gam1 embryo in which the heart tube has failed to undergo looping

Full size image

Embryos that are viable at the late neurula stage usually develop into free-swimming tadpoles (> 48 hpf). Gam1 injected embryos at this stage of development have shortened or bent axes often due to the failure of the neural tube to close properly. Additionally, the majority of embryos display some degree of edema, or fluid collection, in the ventral regions of the embryo specifically around the heart and digestive system (Fig. 2b). Along with edema and a shortened axis, embryos display defects in eye, heart, and gut development, including cyclopia or eyes that are not fully separate (fused eye) (Fig. 2c) and disordered digestive system.

Defects in the heart include improper looping, disruption of chamber formation, and an absence of blood flow. Normal amphibian hearts contain three chambers, two atria and one ventricle, with looping of the outflow track to the right. Hearts of Gam1-injected embryos often lack discernible chambers and seem to have not developed beyond a heart tube structure due to a failure of proper looping (Fig. 2d). Heart structures that have formed contain little to no blood past the stage in development when control matched embryos have begun producing blood. Surprisingly, these embryos are viable and have heart structures that continue to beat, although contractions are erratic. Defects in the formation of these organs indicate that the temporary disruption of SUMOylation activity during early embryogenesis impacts the developmental program long after this activity is restored.

Inhibition of SUMOylation activity disrupts convergence and extension

The phenotypes of a shortened A–P axis and failure of the blastopore and neural tube to close in the Gam1 embryos indicate that SUMOylation activity is necessary for convergence and extension. Information concerning the involvement of SUMO in cell movement is limited. A morpholino oligonucleotide directed at SUMO1 prevented activin-induced elongation of Xenopus animal caps [15] and SUMOylation of Rac1 GTPase is needed for optimal cell migration in response to hepatocyte growth factor [24]. Since cell movement during gastrulation is primarily internal and cannot be easily observed, we turned to Keller explant assays to determine whether the SUMOylation deficiency in Gam1 embryos impacts this activity [25]. Explants are made by sandwiching together a portion of two embryos dorsal to the blastopore lip, which contain migrating cells. This arrangement causes the mesodermal and endodermal cells to elongate in a plane rather than involuting beneath ectodermal cells, enabling quantification of the degree of convergent extension.

Control embryos (injected with H2O) showed normal convergent extension while those injected with 2.5 ng Gam1 mRNA failed to elongate properly (Fig. 3a). Additionally, when the amount of Gam1 mRNA was doubled to 5.0 ng, explants not only showed no evidence of cell movement, but the two explants failed to adhere to each other. This indicates not only a role of SUMOylation in controlling cell migration, but also in controlling cell adhesion.

Fig. 3
figure 3

Gam1 disrupts convergence and extension. a Examples of Keller sandwich explants prepared from embryos injected with (i) water, (ii) 2.5 ng Gam1 mRNA, (iii) 5 ng Gam1 mRNA. Explants from embryos injected with water (n = 17) or 2.5 ng Gam1 mRNA (n = 18) were measured to quantify effects on b convergence and c extension. The red lines designate the lengths measured in the explants. Center lines in plots show the medians; box limits indicate the 25th and 75th percentiles; whiskers extend 1.5 times the interquartile range from the 25th and 75th percentiles

Full size image

The degree of explant convergent extension was quantified by measuring the boundary of the non-involuting and involuting marginal zones (convergence) and the anterior to posterior length (extension) (Fig. 3b, c). Two-tailed t-tests show that both convergence (p < 0.001) and extension (p < 0.00005) are significantly decreased in Gam1 embryos compared to controls. The disruption of these activities explains, at least in part, many of the observed phenotypes of Gam1 embryos.

Gene expression in SUMOylation deficient Xenopus embryos

Gam1 knockdown of SUMOylation activity reproducibly generated distinct developmental phenotypes, several of which have been previously linked to disruptions of this post-translational modification, including neural tube closure [26], cleft lip/palate [27], heart defects [28,29,30], axial mesodermal defects that lead to midline phenotypes such as cyclopia [31], and disruption of hematopoiesis [32]. In good accord with the Gam1 phenotypes, a morpholino antisense oligonucleotide directed against Xenopus SUMO1 caused shortened and bent A-P axes, incomplete closure of the neural tube, microcephaly, and inhibition of activin-induced elongation of animal caps [15]. The close correspondence between the morphlino and Gam1 phenotypes affirm that the effects resulting from Gam1 expression are chiefly due to interference of SUMOylation activity in these embryos.

The ability to decrease SUMOylation activity to nonlethal levels presents the opportunity to examine the role of this protein modification, which frequently targets transcription factors and several other nuclear proteins, in the regulation of gene expression during early development. Three developmental stages were chosen for transcriptome analysis (30,000 genes) by microarray: early gastrula (stage 10), which shortly follows the initiation of zygotic transcription at the MBT; late gastrula (stage 12), which is at the end of a period of complex cell movement that organizes the three germ layers; and early neurula (stage 14), which marks changes in transcription that underlie organogenesis and patterning of the embryo. Three biological replicates were analyzed in order to account for variation between egg clutches. The R-values of scatter plots of the biological replicates ranged from 0.959 to 0.916, demonstrating excellent correlation between each of the replicates. The number of differentially expressed genes (p < 0.05) between control (water injected) and Gam1 embryos numbers 94 (53 down-regulated, 41 up-regulated) at early gastrula; 447 (263 down-regulated, 221 up-regulated) at late gastrula; and 742 (394 down-regulated, 387 up-regulated) at early neurula. A heatmap of all differentially expressed genes identified from comparisons between control and Gam1 injected embryos at the three developmental stages is presented in Additional file 1: Figure S1. An accompanying data matrix contains the gene names and expression values in the order in which they appear in the heatmap (Additional file 2: Table S1).

In order to validate the microarray data, the same RNA samples were used for qRT-PCR assays of genes selected based upon either the magnitude of the change in expression (dmrta1, rpl8, ccng1), the occurrence of the gene in a pathway involved in patterning the early embryo (xbra, foxc1, wnt8b, gsc, chrd), or a previous microarray analysis of X. laevis development (xpo1, wnt8b, krt) [33]. (A graphical comparison of the microarray and qRT-PCR data is presented in Additional file 3: Figure S2). The expected correlation between microarray and qRT-PCR is 80–87% with respect to the direction of the change (up- versus down-regulation) [33, 34]. The correlation between our microarray and qRT-PCR assays is 93%; these results, along with the scatter plots, indicate that the quality of the microarray data is high. There is also good agreement between the two measurements with regard to quantitation, especially for those genes that show the greatest up- (ccng1 and xpo) and down-regulation (dmrta1).

Biological process ontology of differentially expressed genes

Differentially regulated genes for each time point were assembled into gene lists for analysis using data mining software including MetaCore, BiNGO, DAVID, and the Gene Ontology database [35]. Xenopus gene lists were converted to their human equivalents for analyses using these applications. Volcano plots (Additional file 4: Figure S3) display statistical significance (p value) versus fold change at each time point. These plots reflect not only the increasing number of differentially expressed genes progressing from early gastrula to early neurula, but also the increasing magnitude of these differences. With few exceptions, most changes in expression at early gastrula are only 20–30%; whereas, larger quantitative changes are seen at the later two time points. We have used these lists of genes affected by the loss of SUMOylation activity for various bioinformatics analyses in order to understand the regulatory role of this post-translational modification during early vertebrate development.

The first analysis was to determine the predominant (high level) biological processes affected by loss of SUMOylation activity at the three different developmental time points (Fig. 4). The profiles at these three stages, which span a total of about 7 h, are distinctly different. However, the distribution of affected genes among these biological processes seems to reflect the major cellular activities of embryos at each of these stages. With lengthening of the cell cycle, zygotic transcription begins at the MBT and transcriptome analysis at early gastrula is expected to capture this transition. Indeed, many of the differentially expressed genes at this time point normally show robust induction between Nieuwkoop–Faber stages 8 and 10 [36]. The greatest proportion of affected genes at early gastrula are concerned with Regulation (31%) and Metabolism (18%), primarily as a result of the resumption of transcription and protein synthesis.

Fig. 4
figure 4

Biological processes most affected by loss of SUMOylation activity. Functional enrichment analysis of differentially expressed genes at each time point based on GeneGo categorization of annotated protein function. Only selected groups containing the largest number of genes are presented

Full size image

The expression of 9 out of 94 differentially expressed genes are assigned to Developmental Processes in the early gastrula embryos, of which the most notable are foxc1, foxd1 and twist. All three of these genes, which are highly expressed in mesoderm, can be correlated with morphological phenotypes that later appear in the Gam1 embryos (e.g., deficits in heart and neural tube development). Although GeneGo did not place sox2, a key regulator of neural induction, under Developmental Processes, it is nonetheless notable that its expression is reduced 1.7-fold at early gastrula and remains modestly lower (1.4-fold) at late gastrula.

The most striking change between early and late gastrula embryos is the greatly increased percentage of genes that fall under Developmental Processes (increasing from 9 to 34%) and Multicellular Organism Processes (increasing from 10 to 23%). Indeed, only four high level biological processes predominate at late gastrula. At this time point, several additional members of the Fox family of transcription factors appear in the list of differentially expressed genes ascribed to Developmental Processes as well as genes encoding proteins involved in Wnt, ephrin, and TGF-β signaling, which are all known to play critical roles during gastrulation.

Although the percentage of genes ascribed to Development Processes decreases in early neurula embryos to 16% of the total of differentially expressed genes, the actual number increases, showing that SUMOylation continues to play an important regulatory role in development. While the genes for many of the aforementioned transcription factors and signaling molecules remain mis-regulated, now multiple Hox genes are repressed in the Gam1 embryos at this time point, which is consistent with their normally strong induction during late gastrula stage. Whereas, SUMOylation is often associated with transcriptional repression, several of the Hox genes seemingly require this post-translational modification for activation.

Of special note, expression of the homeobox protein, Pitx2, is down 4-fold at early neurula in Gam1 embryos. This transcription factor is critical for eye development, as well as the formation of asymmetric organs such as the heart and gut, whose perturbation is a predominant phenotype of the SUMOylation deficient embryos. Pitx2 along with foxc1 (down 2.7-fold at early gastrula) are both implicated in Axenfeld-Rieger syndrome. The lists of genes that are categorized into the various biological processes can be found in Additional file 5: Table S2.

Transcription factor network building

The regulation of transcription by SUMOylation can occur on multiple levels that include, for example, changes in chromatin structure (e.g., histone modification and nucleosome positioning) and subnuclear localization (e.g., PML bodies). However, the most common mechanism is likely through the SUMOylation of transcription factors. We used the network building algorithm ‘Transcription Factor Regulation’ to identify common transcription factors that control the maximum number of differentially expressed genes in Gam1 embryos. This algorithm adds a transcription factor to a gene list as a seed node, and builds a network of interactions around it. Gene lists from each time point were analyzed separately and the top 30 transcription factors at each time point were returned (Table 1). The connections made within each network are based on direct control of the target gene by the seed transcription factor and, inversely, target gene control of the seed transcription factor. The lists of transcription factors at each time point are quite similar and include those such as Sp1, HNF4-α, c-myc, and p53, that regulate an especially large number of different genes.

Table 1 Transcription Factor Network Buildinga
Full size table

In order to determine if these transcription factors are highly associated with genes controlled by SUMOylation, the same analysis was carried out using the lists of genes not affected by Gam1, which number 11,702 genes from early gastrula, 11,322 genes from late gastrula, and 11,019 genes from early neurula (Table 1). The list of transcription factors connected to changed versus unchanged genes is strikingly similar. At early gastrula, 18 out of 30 transcription factors were common to both sets of genes with the top 5 factors being identical. For the late gastrula, 21 out of 30 factors are common and again the top 5 are identical. At the early neurula, 20 out of 30 factors are common with the top 5 factors remaining identical between the lists. Of the top 5 transcription factors, four are common across all three time points (SP1, HNF4-α, c-myc, p53). The simplest interpretation of this analysis is that these transcription factors emerge solely because they bind to an unusually large number of regulatory sites throughout the genome, a subset of which are regulated by SUMOylation.

Interestingly, the GO biological processes associated with any specific transcription factor are different depending on whether the gene list is built from those changed by Gam1 or those unaffected by Gam1 (Additional file 6: Table S3). For example, the two top transcription factors at early gastrula are Sp1 and HNF4-α. The GO processes most highly associated with Sp1 in Gam1 embryos are eye, sensory organ, and embryo development; whereas, in control embryos they are metabolic processes. It is especially noteworthy that SUMOylation of Sp1 has been directly implicated in Xenopus eye development [37]. For HNF4-α, the enriched GO processes in Gam1 embryos are dendrite morphogenesis, signaling pathways, microtubule polymerization, and cell communication; the GO processes in control embryos are almost entirely metabolic. This analysis indicates that constitutive genes necessary for organism viability, e.g., metabolism and biosynthesis, are mostly controlled in a SUMO-independent manner. Conversely, genes with complex expression patterns are more likely controlled by transcription factors whose activity is modulated by SUMOylation. This GO analysis also suggests that SUMO modification generates distinct pools of a particular transcription factor in order to regulate specific subsets of genes, with SUMOylated forms of these transcription factors biased towards more highly regulated (e.g., cell type, temporal) processes. This view is supported by studies in Drosophila where mutations in the SUMO E1 activating enzyme had no effect on larval growth, but did cause defects in imaginal disk formation [38].

A total of 50 different transcription factors were identified through the Transcription Factor Regulation algorithm as being associated with genes differentially expressed at one or more time points. These factors were analyzed for potential SUMOylation sites and SIMs using the SUMO prediction site algorithm GPS–SUMO [39]. Of the 50 transcription factors, 46 (92%) contain consensus SUMOylation sites, only 1 has no identifiable site, and 3 lack sufficient sequence information to make a prediction (Additional file 7: Table S4). To date, 29 of the 46 have been experimentally verified to be SUMOylated. SIMs appear in 42 of the listed transcription factors. The fact that nearly all of factors identified by network building are targets of SUMOylation and/or contain SIMs indicates that the majority of changes in gene expression induced by Gam1 are due directly to changes in factor activity as a result of this post-translational modification.

The Transcription Factor Regulation algorithm returns the most highly connected factors for a given gene list. A shortcoming of this analysis is that critical transcription factors may be overlooked, if they regulate a small fraction of genes in the entire list. As an example, the gene list at early gastrula was examined manually using a variety of sources to identify candidate transcription factors for each gene. As expected, we found that the transcription factors listed in Table 1 are associated with many of the differentially expressed genes in the list. However, we also found that the transcription factor Nkx2.5, whose role in cardiac differentiation is regulated by SUMOylation [28], has predicted binding sites at the promoters of a subset of these differentially expressed genes (alkbh6, evx1, fam184a, foxc1, foxd1, fzd10a, gapvd1, twist1/2). Of these putative targets of Nkx2.5, foxc1 is particularly significant, since foxc1-depleted Xenopus embryos exhibit similar phenotypes to the Gam1-injected embryos: shortened A–P axis and abnormal heart and gut development [40]. Moreover, a morpholino knockdown of foxc1 protein caused changes in the expression levels of casein kinase 1 epsilon, insulin-like growth factor 2, and keratin 8 mRNAs [40] that we also detect in the early neurula embryos.

The heart defects of Gam1-injected embryos are certainly not due to foxc1 alone. We also measured reductions of other key cardiac regulatory factors such as mespa (1.6-fold), tbx2 (2.1-fold), and tbx3 (2.3-fold) at late gastrula and pitx2 (4.0-fold) at early neurula. Thus, the potential effects of depleted SUMOylation activity on heart development are seemingly propagated through the three stages sampled in this study. Twist, which is down regulated at early gastrula and considerably more at late gastrula, has also been implicated in cardiogenesis (reviewed in [41]). Thus, the predominant phenotypes in heart development can be explained by multiple effects on the cardiac transcription program that results from depletion of SUMOylation activity.

Co-expression analysis

The expression of the differentially expressed genes were compiled across the three time points and compared using the module/cluster detection algorithm implemented in R package WGCNA [42]. The differentially regulated genes were compiled into eight syn-expression clusters based on the similarity of expression changes across all three time points. The gene lists from individual clusters were analyzed using the MetaCore transcriptional regulation algorithm. Each cluster returned a similar set of top transcription factors that are similar to the individual lists returned for both the differentially expressed and unchanged genes (Table 1). This observation further supports the idea that the most highly connected transcription factors are those that are widely used to control the greatest number of genes in any list. This strongly argues that co-expression in this case cannot be traced to a single transcription factor or set of transcription factors acting as master regulators.

Signaling pathways sensitive to diminished SUMOylation activity

We chose a non-lethal level of Gam1 expression in order to identify developmental processes that are potentially regulated by SUMOylation. This strategy will only identify those that are most sensitive to loss of this activity. Because the effects of this post-translational modification are widespread and not limited to gene regulatory proteins, we have used the predominant phenotypes of Gam1-injected embryos as a guide to examine some of the signaling pathways that appear to be the most vulnerable to this perturbation. The list of differentially expressed genes was uploaded into MetaCore and analyzed for their presence in known biochemical pathways. Analysis focused on those pathways with outcomes related to the observed phenotypes of failed blastopore or neural tube closure, shortened A–P axis, and cardiovascular deformations. Three pathways, (i) non-canonical Wnt signaling, (ii) epithelial to mesenchymal transition (EMT), and (iii) the Ets-1 pathway, contain multiple differentially expressed genes and were chosen as examples of cases where disruption of SUMOylation can be related to the observed phenotypes.

Non-canonical Wnt/PCP signaling

Wnt signaling can be divided into two broad groups, canonical and non-canonical, based on outcome. Canonical signaling is dependent on β-catenin and regulates transcription of target genes, while non-canonical is independent of β-catenin and regulates cell polarity and calcium levels [43]. The non-canonical signaling pathway leads to activation of Rho and reorganization of the cytoskeleton and changes in cell polarity. Specific Wnt proteins are classified as canonical or non-canonical based on their ability to induce secondary axes in Xenopus embryos. The latter, which includes Wnt4, Wnt5a, and Wnt11 [44,45,46] have little or no dorsalizing effect and do not activate known Wnt target genes such as nodal3.1 (Xnr3) and sia1 (Siamois), whose expression is not changed in Gam1 embryos.

Several Wnt and Fzd genes are misregulated in SUMOylation deficient embryos (Fig. 5). Three Fzd genes of the non-canonical pathway (fzd2, fzd7 and fzd10) are down regulated at one or more time points. Two Wnt genes (wnt8b and wnt11) are also down regulated. Wnt8b has not been specifically implicated in non-canonical signaling; however, it is often placed in the same group as Wnt11, since both have similar temporal roles, activating canonical pathways in oocytes and pre-MBT embryos, but non-canonical pathways during gastrulation [47, 48]. Both Wnt8 and Wnt11 can bind to Fzd7 [49].

Fig. 5
figure 5

Changes in gene expression in the non-canonical Wnt/PCP signaling pathway. Decreased mRNA levels in Gam1 embryos relative to control embryos is represented by green and increased levels by red. Documented targets of SUMOylation are indicated. Transcriptional regulation is denoted by blue arrows and signaling through protein-protein interactions by black arrows. Log (2) changes in gene expression are listed for the three experimental time points

Full size image

While Fzd2 and Fzd7 have been implicated in planar cell polarity (PCP), the interaction of the latter receptor with Wnt11 ligand appears to be the major determinant of this process [50,51,52,53]. Therefore, the observed decrease in both these receptors and ligands can account for the effect of Gam1 on convergent extension. An early manifestation of convergent extension/PCP is closure of the blastopore, which is generally delayed and often incomplete in Gam1 embryos (Fig. 2a) and is specifically attributed to forces generated by convergent thickening (CT) [54]. Convergent extension underlies the intercalation of mesodermal cells that drives the elongation of the A–P axis, which is also compromised in the SUMOylation deficient embryos. Both Wnt 11 and Fzd7 have been implicated in A–P axis elongation. A constitutively active form of the Rho GTPase, Cdc42, rescues axis elongation effects caused by expression of a truncated form of Fzd7, while a dominant negative form of Cdc42 rescues the inhibition of axis elongation that occurs with Wnt11 and Fzd7 overexpression [50]. In some cases, Fzd7 and Fzd2 act redundantly in convergent extension [55], which likely accounts for the high penetrance of convergent extension phenotypes in Gam1 embryos, since we measure a decrease in expression of both receptors during late gastrula.

Disruption of the non-canonical Wnt pathway can also account for the cardiac and neural tube defects that occur later in embryo development. Wnt 11 and Fzd7 are required at multiple points during Xenopus heart development, not only early in the specification of cardiac progenitor cells, but also in cardiac morphogenesis (septation and outflow tract development) [56,57,58,59,60,61]. The ensuing decreased activity of Dishevelled (Dvl) in Gam1 embryos is fully consistent with the cardiac and neural tube phenotypes of Dvl-deficient mice and Xenopus [62,63,64].

Two other PCP core components, Vangl2 and Prickle3, are modestly reduced at early neurula stage, a period when a complex of these two proteins is required for the radial cell intercalations that sustain neural tube closure [65, 66]. Notably, the polarization of the Prickle3/Vangl2 complex to anterior cell edges has been linked to Wnt11 signaling [67]. Vangl2 also plays a critical role in heart formation that is independent of neural tube closure. In mice, knock down of Vangl2 specifically in undifferentiated cells of the second heart field prevents lengthening of the outflow tract due to a failure of these cells to polarize and differentiate [68]. These early abnormalities lead to later defects in cardiac morphology. Decreased expression of several members of the non-canonical Wnt signaling pathway combined with those of the PCP pathway can account for the most prominent phenotypes of the SUMOylation-deficient Xenopus embryos and strongly suggest that disruptions in this post-translational modification potentially underlie congenital birth defects of the neural tube and heart [30, 69].

Because several Wnt and Fzd genes are down regulated in the Gam1 embryos, it is reasonable to speculate that they share one or more common transcriptional regulators that can coordinate their expression; however, it has not been possible to discern a single master regulator. Transcription factor network building indicated p53 activity as the most likely source for dysregulation of non-canonical Wnt signaling. At all three time points, p53 is in the top five transcription factors most highly connected to differentially expressed genes. Surprisingly, p53 shows a two-fold increase in expression at late gastrula in Gam1 embryos. However, regulation of p53 activity by SUMOylation is complex with the modification having a positive, negative, or no effect depending on the target gene in question [70]. Moreover, ANKRD11, a coactivator of p53 [71], shows a five-fold decrease in expression by early neurula stage. Our results indicate that SUMOylation of p53 and/or one of its partners plays a critical role in the regulation of non-canonical Wnt pathway genes during Xenopus gastrulation and neurulation (Fig. 5). In addition, decreased expression of transcription factors GATA6 [72], Pitx2 [73], and T2 (Xbra3) [74] can also contribute specifically to the decrease in wnt11 mRNA expression (Fig. 5).

Epithelial to mesenchymal transition

The epithelial to mesenchymal transition (EMT) is a hallmark of early vertebrate development and patterning of the embryo [75, 76]. In Xenopus, EMT changes allow for the presumptive mesoderm to migrate inward during gastrulation, and disruptions in EMT signaling lead to developmental defects in gastrulation and later patterning of the embryo [77]. One important pathway involves control of EMT genes by the transcription factors snail 1 (snai1) and twist (twist1) [75, 76, 78]. Although the expression of snai1 and twist1 genes requires high mobility group A-T hook 2 (HMGA2), the two genes appear to be under distinct regulatory control [79]. HMGA2 binding to the snai1 promoter involves formation of a complex with phosphorylated SMADs activated by TGFβ signaling [80]. Conversely, HMGA2 is able to bind directly to the twist1 promoter in the absence of TGFβ signaling [79, 80]. HMGA2 is a known target of SUMOylation [81, 82] that has been directly implicated in Xenopus EMT and migration of neural crest cells [83] as well as Xenopus cardiogenesis [84]. The different behavior of HMGA2 at the snai1 and twist1 promoters is reflected in the different effect of Gam1 on the expression of these two genes. Twist expression was reduced at all three time points. In contrast, no misregulation of snai1 was observed.

The microarray analysis detected disruptions in the twist/snail pathway that potentially increase epithelial and decrease mesenchymal characteristics that contribute to the observed gastrulation defects and failure of the blastopore to close (Fig. 6). Twist and Snail1 form a complex to control expression of the EMT-specific transcription factor, Zeb-2, which is down regulated 1.3-fold at early neurula. Zeb-2 is responsible both for activation of genes that promote a mesenchymal phenotype and repression of genes that promote an epithelial phenotype [85]. Moreover, Zeb-2 is a target of SUMOylation, which controls its activity in a promoter-specific fashion [86]. ZO-2 (tjp2), an integral component of tight junctions normally down-regulated by Zeb-2 in mesenchymal cells, allowing for increased cell movement, is up-regulated at late gastrula and early neurula in Gam1 embryos. Vimentin (vim), claudin 4 (cldn4), and connexin26 (gjb2) normally up regulated in mesenchymal cells, are down regulated in Gam1 embryos [87, 88].

Fig. 6
figure 6

Changes in gene expression in the Twist/Snail pathway. Decreased mRNA levels in Gam1 embryos relative to control embryos is represented by green and increased levels by red. Documented targets of SUMOylation are indicated. Transcriptional regulation is denoted by black arrows. Log (2) changes in gene expression are listed for the three experimental time points

Full size image

SUMO possibly impacts this pathway at a second point. NF-κB/p65 an activator of Twist and Zeb2 expression [89], is regulated by a negative feedback loop involving the SUMOylation of the p65 (RELA) subunit [90] that can account for the decreased expression these two transcription factors. In combination, the altered expression of the proteins in this pathway is expected to contribute to an epithelial-like character that impedes cell migration during gastrulation and neurulation, resulting in phenotypes such as the observed incomplete closure of the neural tube and shortened A–P axis. On the other hand, the constant level of snail1 expression can account for the unperturbed levels that we measure of other regulators of the EMT (e.g., E-cadherin, N-cadherin, fibronectin, and occludin).

Ets-1 pathway

Pathway analysis of the microarray data revealed that perturbation of Ets-1 activity is another likely contributor to the phenotypes exhibited by Gam1 embryos [91,92,93,94]. The regulation of Ets-1 during development is complex and occurs through the integration of multiple types of post-translational modification: phosphorylation, acetylation, ubiquitination, and SUMOylation [92, 94, 95]. Furthermore, Ets factors can either enhance or repress transcription [96]. Ets-1 can be SUMOylated at two lysine residues (K15 and K227), resulting in the loss of transcriptional activation activity [92]. Conversely, Ets-1 phosphorylated at threonine 38 by MAPK, as a result of Ras activation, increases its transcriptional activation activity (Fig. 7). The SUMOylation of two targets in the Ras/Raf-MEK-MAPK phosphorylation pathway indirectly influence Ets-1. Ras-1 is a target for SUMO modification, which is required for Ras-dependent activation of MAPK [97]. In contrast, SUMOylation of the downstream target MEK blocks its interaction and activation of MAPK [98]. Down regulation of several genes controlled by Ets-1 in Gam1 embryos indicates that inactivation of the Ras pathway and the resulting dephosphorylation of Ets-1 has the greatest impact on the activity of this transcription factor.

Fig. 7
figure 7

Changes in gene expression in the Ets-1 pathway. Decreased mRNA levels in Gam1 embryos relative to control embryos is represented by green and increased levels by red. Documented targets of SUMOylation are indicated. Connections by arrows represent activation and T-bars represent inhibition. Log (2) changes in gene expression are listed for the three experimental time points

Full size image

Knockdown of Ets-1 in Xenopus neural crest cells results in defects in delamination of the neural tube and disruption of cell migration that lead to malformed outflow tracts of the heart; whereas, depletion of Ets-1 activity in heart mesoderm results in distinct phenotypes that include loss of the entire endocardium and failure to form discrete chambers [93]. Two Ets-1 target genes, platelet derived growth factor α (pdgfa) and PDGFα receptor (pdgfra), are down regulated. There is abundant evidence that PDGFα is necessary very early in heart development for the proper migration of cardiomyocytes to the midline and formation of the heart tube [99,100,101]. Ets-1 works in concert with Sp1 to activate transcription of the PDGFα gene [102]. Like Ets-1, Sp1 is repressed by SUMOylation [103] and activated by MAPK [104], potentially amplifying the effect of Gam1 on the expression of pdgfa and pdgfra.

Gastrulation defects in developing Xenopus embryos are also observed when PDGFα signaling is disrupted [105,106,107]. PDGFα is expressed in the blastocoel roof of the developing embryo and its receptor is expressed in the migrating mesoderm. This complementary expression of ligand and receptor orients and directs the movement of mesodermal cells inward during gastrulation. Interfering with PDGFα signaling randomizes internal movement of these cells and prevents proper convergent extension [106]. This disruption leads to improper patterning and phenotypes that include incomplete blastopore closure and shortened axis elongation. Mice with mutations in the PDGFαR display severe skeletal abnormalities, cleft face, and spina bifida [108, 109], with the latter being a frequent phenotype in Gam1 embryos.

While the disruption of Ets-1 (and Sp1) activity can account for the immediate effects on PDGFα and its receptor, we note that there are also several downstream changes in gene expression that can contribute to the phenotypes of Gam1 embryos, including DUSP1 (cell migration) [110] and IGF2 (neural and eye development) [111, 112]. Of particular note, is the increased expression of the transcription factor Egr1 (late gastrula) which is a positive regulator of MyoD and a negative regulator of Xbra [113]. Consistent with the change in Egr1, we measure increased expression of MyoD (late neurula) and decreased Xbra expression (late gastrula/early neurula). Both MyoD and Xbra are critical factors in several developmental pathways.

Conclusions

It is now appreciated that SUMO impinges on a substantial number of biological processes, explaining why deletion of E2 conjugating enzyme is embryonic lethal in mice and zebrafish [2, 8, 11,12,13]. For this reason, it has been difficult to assess the role of this post-translational modification in live animals. While the magnitude of many of the changes in gene expression in Gam1 embryos is modest, the cumulative effect of disrupting multiple pathways that impinge upon a particular developmental event accounts for the frequency of the commonly observed phenotypes. The experiments here make evident that regulation by SUMO during the earliest period of embryogenesis is critical to proper development at later stages. Injection of Gam1 mRNA suppressed SUMOylation activity up to late neurula (26 hpf), yet generated phenotypes that appear later in development (e.g., heart looping, 48 hpf), demonstrating the early embryonic origins that give rise to the two most common human birth defects: heart and neural tube. Transcriptome analysis has revealed those pathways that are most susceptible to diminished SUMOylation activity. The frequent occurrence of heart defects in Gam1 embryos is in accord with the prominent role SUMO plays in the regulation of cardiac gene expression [114, 115]. Indeed, in a cohort of 87 babies exhibiting atrial septal defects, five carried mutations that reduced the promoter activity of the SUMO1 gene [30], demonstrating a direct connection between deficits in this posttranslational modification and congenital heart defects. It is also clear from the experiments described here that the PCP pathway is exceptionally sensitive to depressed SUMOylation activity and this most certainly accounts for the frequent neural tube defects of the Gam1 embryos [116, 117] and supports a proposed relationship between the PCP pathway and risk for spina bifida [69]. Our transcriptomics analysis has identified disruptions in pathways that have especially great potential to contribute to these pathologies and that warrant closer clinical examination. Targeted inhibition of SUMOylation activity by Gam1 will be useful for further investigation of this post-translational modification and its causative role in birth defects.

Methods

Embryo preparation

Female X. laevis were injected with 500 units of human chorionic gonadotropin (HCG) at least 12 h prior to spawning. Testes were isolated and stored at 4 °C in high salt MBS (0.7 mM CaCl2, 108 mM NaCl, 1 mM KCl, 1 mM MgSO4, 5 mM HEPES, pH 7.8, 2.5 mM NaHCO3) for up to a week. Minced testis and eggs were mixed together with 3 mL 1/3 MMR (0.1 M NaCl, 2.0 mM KCl, 1 mM MgSO4, 2 mM CaCl2, 5 mM HEPES, pH 7.8) and left to fertilize at room temperature. After 10 min, eggs were flooded with additional 1/3 MMR and after an additional 30 min eggs were washed with 2% cysteine (< 5 min) in order to remove the jelly coating, upon which eggs were washed 8 times with 1/3 MMR. Embryos were injected in the animal hemisphere (volumes ranged between 5 to 30 nl) and allowed to develop at room temperature in 1/3 MMR. At least 30 min prior to harvest of the testes, the male frog is submerged in an ice water bath using a mesh basket to keep the ice away from the skin surface. The frog is decapitated, and pithed prior to removal of the testes. After removal of the testes tissue, the heart is cut and the frog pithed to ensure death prior to disposal. Female frogs are rested at least 3 months between spawnings.

RNA synthesis

Linearized plasmid (2.5 μg) was added to a mixture containing transcription buffer (Promega), 10 mM DTT, 2.5 mM cap analog, 50 units RNasin, 200 units RNA polymerase, and 0.5 mM each UTP, ATP, CTP, GTP in a total volume of 50 uL. The reaction was incubated at 37 °C for 90 min, followed by the addition of another 200 units of polymerase and incubation for 90 additional minutes. The reaction mixture was extracted twice with phenol pH 4.5 and twice with 24:1 chloroform: isoamylalcohol. The mRNA was precipitated with ethanol, suspended in 50 ul of H2O and run though a NucAway spin column (Ambion) to remove unincorporated nucleotides. mRNA concentrations were determined using a NanoDrop spectrometer (Thermo Scientific), and samples were stored at − 80 °C until use.

SUMOylation assays

Gam1 mRNA was injected into one-cell embryos, which were allowed to develop to the indicated stage (midblastula, late blastula, early gastrula, and midneurula). Twenty embryos per time point were homogenized in 22 uL of SUMO reaction buffer (Boston Biochem) and spun at 16,000 rcf for 5 min. Supernatant was removed and the protein concentration was determined by Bradford assay. Each assay conatined 25 μg of whole cell extract, 60 uM SUMO1, 5 uM Ubc9, 5 uM E2-25K (substrate peptide), 25 mM Mg-ATP in a total reaction volume of 25 uL. Control reactions contained 500 nM SAE1/SAE2 (E1 enzyme) in place of cell extract. Reactions were incubated at 37 °C for four hours. SDS loading dye was mixed with the reactions and vortexed briefly in order to stop the reaction. Assays were analyzed by SDS polyacrylamide gel electrophoresis followed by western blot.

Western blots

Gam1 carries a single (10 amino acid) N-terminal myc tag that was detected using c-myc antibody (sc-40, Santa Cruz Biotechnology) at 1:500 dilution. For SUMOylation assays an antibody (A-603, Boston Biochem) directed against the SUMO1 substrate peptide (E-25 K) was used at 1:500 dilution. All blots were visualized using a goat anti-rabbit IgG:alkaline phosphatase fusion protein at 1:3000 dilution in conjunction with BCIP/NBT color development (Bio-Rad).

Keller sandwich explant assays

One-cell embryos injected with either Gam1 mRNA (2.5 ng or 5.0 ng) or an equivalent volume of H2O and allowed to develop in 1/3 MMR until early gastrula stage (~ 9 hpf) when the blastopore lip was initially beginning to form. Embryos were maintained in 1× MBS buffer (88 mM NaCl, 1 mM KCl, 0.7 mM CaCl2, 1 mM MgSO4, 5 mM HEPES pH 7.8, 2.5 mM NaHCO3) during the dissection, which followed the procedure described, by Keller and Danilchik (1998).

Once the sandwich was fabricated, it was cultured in the dissection dish for one hour allowing the explants to adhere. Explant sandwiches were then moved to an agarose-coated Petri dish containing 1X MBS supplemented with 50 μg/mL gentamicin and left to culture for an additional 14 h. Samples were photographed immediately after moving them to the Petri dish (10 hpf) and again at 24 hpf using an Olympus SZX16 microscope. Scale bars on the CellSens software were used to measure the length and width of the sandwich at the mentioned time points. Explants from water-injected or Gam1-injected embryos were compared to each other to determine the level of convergence and extension. Sample sets were analyzed for significant differences using a Student’s T-Test.

Microarray analysis

One cell embryos were injected with either 0.5 ng Gam1 mRNA in a total volume of 9.2 nL or with H2O as a control. Embryos were allowed to develop in 1/3 MMR at room temperature until the desired time points. Samples were taken at early gastrula (9 hpf), late gastrula (13.5 hpf) and early neurula (16.5 hpf). Twenty embryos from each time point were homogenized in 500 uL proteinase K buffer (100 mM Tris pH 8, 150 mM NaCl, 12.5 mM EDTA pH 8.0, 1% SDS, 200 μg proteinase K) and incubated at 37 °C for one hour. Samples were then extracted twice with an equal volume of phenol, pH 4.5, and twice with 24:1 chloroform:isoamylalcohol. RNA was purified using an RNAEasy spin column (Qiagen, #74104) and precipitated with 2.5 volumes of ice-cold ethanol. The concentration and purity of the samples was determined using a NanoDrop spectrometer. Prior to microarray data collection, the RNA integrity number (RIN) for each sample was determined by the Notre Dame Genomics and Bioinformatics Core Facility (GBCF) using an Agilent 2100 BioAnalyzer.

Total RNA (50–500 ng) was used for cDNA synthesis according to the GeneChip 3′ IVT Express Kit for Affymetrix microarray chips. The cDNA was then converted to double-stranded DNA and used as a template for transcription of amplified RNA (aRNA) which was labeled with a biotin-conjugated nucleotide. aRNA samples (15 μg) were subsequently purified and fragmented to produce probes for hybridization onto Xenopus laevis 2.0 GeneChip 3’expression arrays. Three biological replicates were carried out for Gam1 and H2O injected samples at each time point utilizing a total of 18 microarray chips. Analysis of the microarray chips was carried out with the Affymetrix GeneChip System. Affymetrix .cel files were analyzed using the Bioconductor software package (http://www.bioconductor.org/). The mean fluorescence intensity was derived from a log2 transformation of the data and normalized using the quantile normalization method. A Student’s t test was used to determine if there was a significant difference between the Gam1 and H2O embryo gene expression values. A p value less than 0.05 was used as a cut off for differential expression. The fold change, between control and Gam1 samples, for each expression value was also calculated. Gene lists for those genes differentially expressed (p < 0.05) were compiled for bioinformatics analysis. Since the Xenopus laevis genome is not fully annotated, some probes on the 2.0 chip, which do not have identified names or functions; only annotated genes were used for bioinformatics analysis.

Microarray data analysis: MetaCore

Gene lists were uploaded to the MetaCore software suite of programs; only annotated genes were used for each analysis and genes were grouped differently depending on the analysis. Genes were grouped into lists based on differential expression at each time point (EG, LG, EN), not-differential expression at each time point (EG, LG, EN), or common expression patterns across all three time points (Cluster 1–8). MetaCore was used to analyze the gene lists for (i) transcription factor networks, (ii) functional enrichment, and (iii) pathway analysis each using an algorithm with specific well-defined instructions for creating the network based on established interactions.

In order to remove the redundancies created by the MetaCore enrichment algorithms, GeneGo IDs identified in the initial analysis were further categorized using GO Term classification software CateGOrizer (v. 3.218). Gene ontologies are hierarchically organized and the CateGOrizer software allows the specific sub-terms to be grouped into their broader categories providing a more accurate enrichment picture.

Pathway maps in MetaCore are created by compiling known biochemical processes or signaling cascades into graphical images. These maps can be used to determine the interconnectedness of uploaded genes in a biological context. Additionally, pathway maps can be used to determine downstream gene targets not present on the microarray. Gene lists were uploaded into MetaCore while pathways relevant to development (heart and nervous system), cell adhesion, cell cycle, cytoskeleton remodeling, and epithelial to mesenchymal transition were manually searched for any matching genes. Pathways containing differentially regulated genes were analyzed further for their connection to the observed developmental phenotypes.

Identification of SUMO consensus sequences and SIMs

Potential sites of SUMOylation and SUMO interaction motifs (SIM) were identifed using GPS–SUMO (http://sumosp.biocuckoo.org/) using the low threshold setting for the former. X. laevis reference sequences were accessed through XenBase (http://www.xenbase.org).

Real time polymerase chain reaction

RNA samples (4 μg) were mixed with random hexamers (Promega, #C118A) and heated at 70 °C for 5 min. Samples were then placed on ice for at least 5 min. Reverse transcription reactions containing 1 uL GoScript reverse transcriptase (Promega, #A501C) in a final volume of 20 uL were placed in a thermal cycler with program sequence: 25 °C for 10 min, 42 °C for 50 min, 70 °C for 15 min. Dilution calculations of the reverse transcription (RT) reaction were based on a proportional amount of cDNA being created from the total RNA added to the reaction thereby producing a solution with a cDNA concentration of 200 ng/uL. RT samples were diluted with H2O to give a final cDNA concentration of 5 ng/uL. An aliquot of 50 uL of cDNA was mixed with an additional 200 uL of H2O to give a final working concentration of 1 ng/uL which was used in SYBR green PCR reactions. An additional 50 uL aliquot of each 5 μg/uL RT sample was taken in order to make standard curve dilutions for testing primer efficiency. A four point standard curve was created with 10 fold dilutions from the 5 ng/uL cDNA arbitrarily termed: 1000 (undiluted), 100, 10, and 1. A standard curve was generated for each primer pair in order to test the primer efficiency. A master mix containing 25 uM of both forward and reverse primers and 12.5 uL SYBR Green Master Mix (Applied Biosystems, #4309155) was prepared for each sample. The master mix (15 uL) and cDNA (10 ul) was added to each well of a 96 well plate. Samples were run with the following program parameters: 50 °C for 2 min, 95 °C for 15 min, 95 °C for 15 s, 60 °C for 30 s, 72 °C for 30 s (last three steps repeated for 45 cycles), 95 °C for 15 s, 60 °C for 1 min, ramping temperature for 20 min, 95 °C for 15 s. The ramping temperature function was used to produce melt curves for calculation of primer efficiency.

Abbreviations

Dvl:

Dishevelled

GO:

Gene ontology

MBT:

Midblastula transition

PCP:

Planar cell polarity

SIM:

SUMO interaction motif

References

  1. Chymkowitch P, Nguea PA, Enserink JM. SUMO-regulated transcription: challenging the dogma. BioEssays. 2015;37:1095–105.

    Article  CAS  PubMed  Google Scholar 

  2. Eifler K, Vertegaal AC. Mapping the SUMOylated landscape. FEBS J. 2015;282:3669–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Flotho A, Melchior F. Sumoylation: a regulatory protein modification in health and disease. Annu Rev Biochem. 2013;82:357–85.

    Article  CAS  PubMed  Google Scholar 

  4. Hendriks IA, D'Souza RC, Yang B, Verlaan-de Vries M, Mann M, Vertegaal AC. Uncovering global SUMOylation signaling networks in a site-specific manner. Nat Struct Mol Biol. 2014;21:927–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Makhnevych T, Sydorskyy Y, Xin X, Srikumar T, Vizeacoumar FJ, Jeram SM, Li Z, Bahr S, Andrews BJ, Boone C, et al. Global map of SUMO function revealed by protein-protein interaction and genetic networks. Mol Cell. 2009;33:124–35.

    Article  CAS  PubMed  Google Scholar 

  6. Hay RT. SUMO: a history of modification. Mol Cell. 2005;18:1–12.

    Article  CAS  PubMed  Google Scholar 

  7. Wohlschlegel JA, Johnson ES, Reed SI, Yates JR 3rd. Global analysis of protein sumoylation in Saccharomyces cerevisiae. J Biol Chem. 2004;279:45662–8.

    Article  CAS  PubMed  Google Scholar 

  8. Eifler K, Vertegaal AC. SUMOylation-mediated regulation of cell cycle progression and Cancer. Trends Biochem Sci. 2015;40:779–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Iñiguez-Lluhi JA. SUMO Modification and Transcriptional Regulation. In: Wilson VG, editor. SUMO Regulation of Cellular Processes. Netherlands: Springer; 2009. p. 13–40.

    Chapter  Google Scholar 

  10. Ma L, Aslanian A, Sun H, Jin M, Shi Y, Yates JR 3rd, Hunter T. Identification of small ubiquitin-like modifier substrates with diverse functions using the Xenopus egg extract system. Mol Cell Proteomics. 2014;13:1659–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Lomeli H, Vazquez M. Emerging roles of the SUMO pathway in development. Cell Mol Life Sci. 2011;68:4045–64.

    Article  CAS  PubMed  Google Scholar 

  12. Nacerddine K, Lehembre F, Bhaumik M, Artus J, Cohen-Tannoudji M, Babinet C, Pandolfi PP, Dejean A. The SUMO pathway is essential for nuclear integrity and chromosome segregation in mice. Dev Cell. 2005;9:769–79.

    Article  CAS  PubMed  Google Scholar 

  13. Nowak M, Hammerschmidt M. Ubc9 regulates mitosis and cell survival during zebrafish development. Mol Biol Cell. 2006;17:5324–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Evdokimov E, Sharma P, Lockett SJ, Lualdi M, Kuehn MR. Loss of SUMO1 in mice affects RanGAP1 localization and formation of PML nuclear bodies, but is not lethal as it can be compensated by SUMO2 or SUMO3. J Cell Sci. 2008;121(Pt 24):4106–13.

    Article  CAS  PubMed  Google Scholar 

  15. Yukita A, Michiue T, Danno H, Asashima M. XSUMO-1 is required for normal mesoderm induction and axis elongation during early Xenopus development. Dev Dyn. 2007;236:2757–66.

    Article  CAS  PubMed  Google Scholar 

  16. Zhang FP, Mikkonen L, Toppari J, Palvimo JJ, Thesleff I, Janne OA. Sumo-1 function is dispensable in normal mouse development. Mol Cell Biol. 2008;28:5381–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wang L, Wansleeben C, Zhao S, Miao P, Paschen W, Yang W. SUMO2 is essential while SUMO3 is dispensable for mouse embryonic development. EMBO Rep. 2014;15:878–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Boggio R, Colombo R, Hay RT, Draetta GF, Chiocca S. A mechanism for inhibiting the SUMO pathway. Mol Cell. 2004;16:549–61.

    Article  CAS  PubMed  Google Scholar 

  19. Boggio R, Passafaro A, Chiocca S. Targeting SUMO E1 to ubiquitin ligases: a viral strategy to counteract sumoylation. J Biol Chem. 2007;282:15376–82.

    Article  CAS  PubMed  Google Scholar 

  20. Malik MQ, Bertke MM, Huber PW. Small ubiquitin-like modifier (SUMO)-mediated repression of the Xenopus oocyte 5 S rRNA genes. J Biol Chem. 2014;289:35468–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Session AM, Uno Y, Kwon T, Chapman JA, Toyoda A, Takahashi S, Fukui A, Hikosaka A, Suzuki A, Kondo M, et al. Genome evolution in the allotetraploid frog Xenopus laevis. Nature. 2016;538:336–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Delaune E, Lemaire P, Kodjabachian L. Neural induction in Xenopus requires early FGF signalling in addition to BMP inhibition. Development. 2005;132:299–310.

    Article  CAS  PubMed  Google Scholar 

  23. Keller RE. The cellular basis of epiboly: an SEM study of deep-cell rearrangement during gastrulation in Xenopus laevis. J Embryol Exp Morphol. 1980;60:201–34.

    CAS  PubMed  Google Scholar 

  24. Castillo-Lluva S, Tatham MH, Jones RC, Jaffray EG, Edmondson RD, Hay RT, Malliri A. SUMOylation of the GTPase Rac1 is required for optimal cell migration. Nat Cell Biol. 2010;12:1078–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Keller R, Danilchik M. Regional expression, pattern and timing of convergence and extension during gastrulation of Xenopus laevis. Development. 1988;103:193–209.

    Article  CAS  PubMed  Google Scholar 

  26. Beaudin AE, Stover PJ. Insights into metabolic mechanisms underlying folate-responsive neural tube defects: a Minireview. Birth Defects Res A Clin Mol Teratol. 2009;85:274–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Alkuraya FS, Saadi I, Lund JJ, Turbe-Doan A, Morton CC, Maas RL. SUMO1 haploinsufficiency leads to cleft lip and palate. Science. 2006;313:1751.

    Article  PubMed  Google Scholar 

  28. Kim EY, Chen L, Ma Y, Yu W, Chang J, Moskowitz IP, Wang J. Expression of sumoylation deficient Nkx2.5 mutant in Nkx2.5 haploinsufficient mice leads to congenital heart defects. PLoS One. 2011;6:e20803.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kim EY, Chen L, Ma Y, Yu W, Chang J, Moskowitz IP, Wang J. Enhanced desumoylation in murine hearts by overexpressed SENP2 leads to congenital heart defects and cardiac dysfunction. J Mol Cell Cardiol. 2012;52(3):638–49.

    Article  CAS  PubMed  Google Scholar 

  30. Wang J, Chen L, Wen S, Zhu H, Yu W, Moskowitz IP, Shaw GM, Finnell RH, Schwartz RJ. Defective sumoylation pathway directs congenital heart disease. Birth Defects Res A Clin Mol Teratol. 2011;91(6):468–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Moreno-Ayala R, Schnabel D, Salas-Vidal E, Lomeli H. PIAS-like protein Zimp7 is required for the restriction of the zebrafish organizer and mesoderm development. Dev Biol. 2015;403:89–100.

    Article  CAS  PubMed  Google Scholar 

  32. Lee H-Y, Johnson KD, Fujiwara T, Boyer ME, Kim S-I, Bresnick EH. Controlling hematopoiesis through Sumoylation-dependent regulation of a GATA factor. Mol Cell. 2009;36:984–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Altmann CR, Bell E, Sczyrba A, Pun J, Bekiranov S, Gaasterland T, Brivanlou AH. Microarray-based analysis of early development in Xenopus laevis. Dev Bio. 2001;236:64–75.

    Article  CAS  Google Scholar 

  34. Dallas PB, Gottardo NG, Firth MJ, Beesley AH, Hoffmann K, Terry PA, Freitas JR, Boag JM, Cummings AJ, Kees UR. Gene expression levels assessed by oligonucleotide microarray analysis and quantitative real-time RT-PCR -- how well do they correlate? BMC Genomics. 2005;6:59.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25:25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Nieuwkoop PD, Faber J. Normal tables of Xenopus laevis Daudin. Amsterdam: Elsevier-North Holland; 1956.

    Google Scholar 

  37. Terada K, Furukawa T. Sumoylation controls retinal progenitor proliferation by repressing cell cycle exit in Xenopus laevis. Dev Biol. 2010;347:180–94.

    Article  CAS  PubMed  Google Scholar 

  38. Kanakousaki K, Gibson MC. A differential requirement for SUMOylation in proliferating and non-proliferating cells during Drosophila development. Development. 2012;139:2751–62.

    Article  CAS  PubMed  Google Scholar 

  39. Zhao Q, Xie Y, Zheng Y, Jiang S, Liu W, Mu W, Liu Z, Zhao Y, Xue Y, Ren J. GPS-SUMO: a tool for the prediction of sumoylation sites and SUMO-interaction motifs. Nucleic Acids Res. 2014;42(W1):W325–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Cha JY, Birsoy B, Kofron M, Mahoney E, Lang S, Wylie C, Heasman J. The role of FoxC1 in early Xenopus development. Dev Dyn. 2007;236:2731–41.

    Article  CAS  PubMed  Google Scholar 

  41. VanDusen NJ, Firulli AB. Twist factor regulation of non-cardiomyocyte cell lineages in the developing heart. Differentiation. 2012;84:79–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9(1):559.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Rao TP, Kuhl M. An updated overview on Wnt signaling pathways: a prelude for more. Circ Res. 2010;106:1798–806.

    Article  CAS  PubMed  Google Scholar 

  44. Du SJ, Purcell SM, Christian JL, McGrew LL, Moon RT. Identification of distinct classes and functional domains of Wnts through expression of wild-type and chimeric proteins in Xenopus embryos. Mol Cell Biol. 1995;15:2625–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Torres MA, Yang-Snyder JA, Purcell SM, DeMarais AA, McGrew LL, Moon RT. Activities of the Wnt-1 class of secreted signaling factors are antagonized by the Wnt-5A class and by a dominant negative cadherin in early Xenopus development. J Cell Biol. 1996;133:1123–37.

    Article  CAS  PubMed  Google Scholar 

  46. Ungar AR, Kelly GM, Moon RT. Wnt4 affects morphogenesis when misexpressed in the zebrafish embryo. Mech Dev. 1995;52:153–64.

    Article  CAS  PubMed  Google Scholar 

  47. Cui Y, Brown JD, Moon RT, Christian JL. Xwnt-8b: a maternally expressed Xenopus Wnt gene with a potential role in establishing the dorsoventral axis. Development. 1995;121:2177–86.

    Article  CAS  PubMed  Google Scholar 

  48. Tao Q, Yokota C, Puck H, Kofron M, Birsoy B, Yan D, Asashima M, Wylie CC, Lin X, Heasman J. Maternal wnt11 activates the canonical wnt signaling pathway required for axis formation in Xenopus embryos. Cell. 2005;120:857–71.

    Article  CAS  PubMed  Google Scholar 

  49. Medina A, Reintsch W, Steinbeisser H. Xenopus frizzled 7 can act in canonical and non-canonical Wnt signaling pathways: implications on early patterning and morphogenesis. Mech Dev. 2000;92:227–37.

    Article  CAS  PubMed  Google Scholar 

  50. Djiane A, Riou J, Umbhauer M, Boucaut J, Shi D. Role of frizzled 7 in the regulation of convergent extension movements during gastrulation in Xenopus laevis. Development. 2000;127:3091–100.

    Article  CAS  PubMed  Google Scholar 

  51. Penzo-Mendez A, Umbhauer M, Djiane A, Boucaut JC, Riou JF. Activation of Gβγ signaling downstream of Wnt-11/Xfz7 regulates Cdc42 activity during Xenopus gastrulation. Dev Biol. 2003;257:302–14.

    Article  CAS  PubMed  Google Scholar 

  52. Tanegashima K, Zhao H, Dawid IB. WGEF activates rho in the Wnt-PCP pathway and controls convergent extension in Xenopus gastrulation. EMBO J. 2008;27:606–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Zallen JA. Planar polarity and tissue morphogenesis. Cell. 2007;129:1051–63.

    Article  CAS  PubMed  Google Scholar 

  54. Shook DR, Kasprowicz EM, Davidson LA, Keller R. Large, long range tensile forces drive convergence during Xenopus blastopore closure and body axis elongation. eLife. 2018;7.

  55. Yu H, Ye X, Guo N, Nathans J. Frizzled 2 and frizzled 7 function redundantly in convergent extension and closure of the ventricular septum and palate: evidence for a network of interacting genes. Development. 2012;139:4383–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Abu-Elmagd M, Mulvaney J, Wheeler GN. Frizzled-7 is required for Xenopus heart development. Biol Open. 2017;6:1861–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Garriock RJ, D'Agostino SL, Pilcher KC, Krieg PA. Wnt11-R, a protein closely related to mammalian Wnt11, is required for heart morphogenesis in Xenopus. Dev Biol. 2005;279:179–92.

    Article  CAS  PubMed  Google Scholar 

  58. Gessert S, Kuhl M. The multiple phases and faces of wnt signaling during cardiac differentiation and development. Circ Res. 2010;107:186–99.

    Article  CAS  PubMed  Google Scholar 

  59. Pandur P, Lasche M, Eisenberg LM, Kuhl M. Wnt-11 activation of a non-canonical Wnt signalling pathway is required for cardiogenesis. Nature. 2002;418:636–41.

    Article  CAS  PubMed  Google Scholar 

  60. Zhang Z, Rankin SA, Zorn AM. Different thresholds of Wnt-frizzled 7 signaling coordinate proliferation, morphogenesis and fate of endoderm progenitor cells. Dev Biol. 2013;378:1–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. van Vliet PP, Lin L, Boogerd CJ, Martin JF, Andelfinger G, Grossfeld PD, Evans SM. Tissue specific requirements for WNT11 in developing outflow tract and dorsal mesenchymal protrusion. Dev Biol. 2017;429:249–59.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. Gray RS, Bayly RD, Green SA, Agarwala S, Lowe CJ, Wallingford JB. Diversification of the expression patterns and developmental functions of the Dishevelled gene family during chordate evolution. Dev Dyn. 2009;238:2044–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Hamblet NS, Lijam N, Ruiz-Lozano P, Wang J, Yang Y, Luo Z, Mei L, Chien KR, Sussman DJ, Wynshaw-Boris A. Dishevelled 2 is essential for cardiac outflow tract development, somite segmentation and neural tube closure. Development. 2002;129:5827–38.

    Article  CAS  PubMed  Google Scholar 

  64. Wallingford JB, Harland RM. Neural tube closure requires Dishevelled-dependent convergent extension of the midline. Development. 2002;129:5815–25.

    Article  CAS  PubMed  Google Scholar 

  65. Ossipova O, Chu CW, Fillatre J, Brott BK, Itoh K, Sokol SY. The involvement of PCP proteins in radial cell intercalations during Xenopus embryonic development. Dev Biol. 2015;408:316–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Ossipova O, Kim K, Sokol SY. Planar polarization of Vangl2 in the vertebrate neural plate is controlled by Wnt and myosin II signaling. Biol Open. 2015;4:722–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Chu CW, Sokol SY. Wnt proteins can direct planar cell polarity in vertebrate ectoderm. eLife. 2016;5.

  68. Ramsbottom SA, Sharma V, Rhee HJ, Eley L, Phillips HM, Rigby HF, Dean C, Chaudhry B, Henderson DJ. Vangl2-regulated polarisation of second heart field-derived cells is required for outflow tract lengthening during cardiac development. PLoS Genet. 2014;10:e1004871.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Wen S, Zhu H, Lu W, Mitchell LE, Shaw GM, Lammer EJ, Finnell RH. Planar cell polarity pathway genes and risk for spina bifida. Am J Med Genet A. 2010;152A:299–304.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Hock A, Vousden KH. Regulation of the p53 pathway by ubiquitin and related proteins. Int J Biochem Cell Biol. 2010;42:1618–21.

    Article  CAS  PubMed  Google Scholar 

  71. Neilsen PM, Cheney KM, Li CW, Chen JD, Cawrse JE, Schulz RB, Powell JA, Kumar R, Callen DF. Identification of ANKRD11 as a p53 coactivator. J Cell Sci. 2008;121:3541–52.

    Article  CAS  PubMed  Google Scholar 

  72. Afouda BA, Martin J, Liu F, Ciau-Uitz A, Patient R, Hoppler S. GATA transcription factors integrate Wnt signalling during heart development. Development. 2008;135:3185–90.

    Article  CAS  PubMed  Google Scholar 

  73. Zhou W, Lin L, Majumdar A, Li X, Zhang X, Liu W, Etheridge L, Shi Y, Martin J, Van de Ven W, et al. Modulation of morphogenesis by noncanonical Wnt signaling requires ATF/CREB family-mediated transcriptional activation of TGFβ2. Nat Genet. 2007;39:1225–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Smith JC, Conlon FL, Saka Y, Tada M. Xwnt11 and the regulation of gastrulation in Xenopus. Philos Trans R Soc Lond Ser B Biol Sci. 2000;355(1399):923–30.

    Article  CAS  Google Scholar 

  75. Kang Y, Massague J. Epithelial-mesenchymal transitions: twist in development and metastasis. Cell. 2004;118:277–9.

    Article  CAS  PubMed  Google Scholar 

  76. Wang Y, Shi J, Chai K, Ying X, Zhou BP. The role of snail in EMT and tumorigenesis. Curr Cancer Drug Targets. 2013;13:963–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Keller R, Shook D. Dynamic determinations: patterning the cell behaviours that close the amphibian blastopore. Philos Trans R Soc Lond Ser B Biol Sci. 2008;363:1317–32.

    Article  Google Scholar 

  78. Hotz B, Arndt M, Dullat S, Bhargava S, Buhr HJ, Hotz HG. Epithelial to mesenchymal transition: expression of the regulators snail, slug, and twist in pancreatic cancer. Clin Cancer Res. 2007;13:4769–76.

    Article  CAS  PubMed  Google Scholar 

  79. Tan EJ, Thuault S, Caja L, Carletti T, Heldin CH, Moustakas A. Regulation of transcription factor twist expression by the DNA architectural protein high mobility group A2 during epithelial-to-mesenchymal transition. J Biol Chem. 2012;287:7134–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Thuault S, Tan EJ, Peinado H, Cano A, Heldin CH, Moustakas A. HMGA2 and Smads co-regulate SNAIL1 expression during induction of epithelial-to-mesenchymal transition. J Biol Chem. 2008;283:33437–346.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Cao X, Clavijo C, Li X, Lin HH, Chen Y, Shih HM, Ann DK. SUMOylation of HMGA2: selective destabilization of promyelocytic leukemia protein via proteasome. Mol Cancer Ther. 2008;7:923–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Wu J, Kim S, Kwak MS, Jeong JB, Min HJ, Yoon H-G, Ahn J-H, Shin J-S. High mobility group Nucleosomal binding domain 2 (HMGN2) SUMOylation by the SUMO E3 ligase PIAS1 decreases the binding affinity to nucleosome Core particles. J Biol Chem. 2014;289:20000–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Macri S, Simula L, Pellarin I, Pegoraro S, Onorati M, Sgarra R, Manfioletti G, Vignali R. Hmga2 is required for neural crest cell specification in Xenopus laevis. Dev Biol. 2016;411:25–37.

    Article  CAS  PubMed  Google Scholar 

  84. Monzen K, Ito Y, Naito AT, Kasai H, Hiroi Y, Hayashi D, Shiojima I, Yamazaki T, Miyazono K, Asashima M, et al. A crucial role of a high mobility group protein HMGA2 in cardiogenesis. Nature Cell Biol. 2008;10:567–74.

    Article  CAS  PubMed  Google Scholar 

  85. Vandewalle C, Comijn J, De Craene B, Vermassen P, Bruyneel E, Andersen H, Tulchinsky E, Van Roy F, Berx G. SIP1/ZEB2 induces EMT by repressing genes of different epithelial cell-cell junctions. Nucleic Acids Res. 2005;33:6566–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Long J, Zuo D, Park M. Pc2-mediated sumoylation of Smad-interacting protein 1 attenuates transcriptional repression of E-cadherin. J Biol Chem. 2005;280(42):35477–89.

    Article  CAS  PubMed  Google Scholar 

  87. Ma X, Miao H, Jing B, Pan Q, Zhang H, Chen Y, Zhang D, Liang Z, Wen Z, Li M. Claudin-4 controls the proliferation, apoptosis, migration and in vivo growth of MCF-7 breast cancer cells. Oncol Rep. 2015;34:681–90.

    Article  CAS  PubMed  Google Scholar 

  88. Polusani SR, Kalmykov EA, Chandrasekhar A, Zucker SN, Nicholson BJ. Cell coupling mediated by connexin 26 selectively contributes to reduced adhesivity and increased migration. J Cell Sci. 2016;129:4399–410.

    CAS  PubMed  PubMed Central  Google Scholar 

  89. Pires BRB, Mencalha AL, Ferreira GM, de Souza WF, Morgado D. íaz JA, Maia AM, Corrêa S, Abdelhay ESFW. NF-kappaB is involved in the regulation of EMT genes in breast Cancer cells. PLoS One. 2017;12(1):e0169622.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  90. Liu Y, Bridges R, Wortham A, Kulesz-Martin M. NF-k B repression by PIAS3 mediated RelA SUMOylation. PLoS One. 2012;7(5):e37636.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Hollenhorst PC, Chandler KJ, Poulsen RL, Johnson WE, Speck NA, Graves BJ. DNA specificity determinants associate with distinct transcription factor functions. PLoS Genet. 2009;5:e1000778.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  92. Ji Z, Degerny C, Vintonenko N, Deheuninck J, Foveau B, Leroy C, Coll J, Tulasne D, Baert JL, Fafeur V. Regulation of the Ets-1 transcription factor by sumoylation and ubiquitinylation. Oncogene. 2007;26:395–406.

    Article  CAS  PubMed  Google Scholar 

  93. Nie S, Bronner ME. Dual developmental role of transcriptional regulator Ets1 in Xenopus cardiac neural crest vs. heart mesoderm. Cardiovasc Res. 2015;106:67–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Tootle TL, Rebay I. Post-translational modifications influence transcription factor activity: a view from the ETS superfamily. BioEssays. 2005;27:285–98.

    Article  CAS  PubMed  Google Scholar 

  95. Chen J, Fu Y, Day DS, Sun Y, Wang S, Liang X, Gu F, Zhang F, Stevens SM, Zhou P, et al. VEGF amplifies transcription through ETS1 acetylation to enable angiogenesis. Nat Commun. 2017;8(1):383.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Sharrocks AD. The ETS-domain transcription factor family. Nat Rev Mole Cell Biol. 2001;2(11):827–37.

    Article  CAS  Google Scholar 

  97. Nie M, Xie Y, Loo JA, Courey AJ. Genetic and proteomic evidence for roles of Drosophila SUMO in cell cycle control, Ras signaling, and early pattern formation. PLoS One. 2009;4:e5905.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  98. Kubota Y, O'Grady P, Saito H, Takekawa M. Oncogenic Ras abrogates MEK SUMOylation that suppresses the ERK pathway and cell transformation. Nat Cell Biol. 2011;13:282–91.

    Article  CAS  PubMed  Google Scholar 

  99. Bax NAM, Bleyl SB, Gallini R, Wisse LJ, Hunter J, Van Oorschot AAM, Mahtab EAF, Lie-Venema H, Goumans M-J, Betsholtz C, et al. Cardiac malformations in Pdgfralpha mutant embryos are associated with increased expression of WT1 and Nkx2.5 in the second heart field. Dev Dyn. 2010;239(8):2307–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Bloomekatz J, Singh R, Prall OWJ, Dunn AC, Vaughan M, Loo C-S, Harvey RP, Yelon D. Platelet-derived growth factor (PDGF) signaling directs cardiomyocyte movement toward the midline during heart tube assembly. eLife. 2017;6:e21172.

    Article  PubMed  PubMed Central  Google Scholar 

  101. El-Rass S, Eisa-Beygi S, Khong E, Brand-Arzamendi K, Mauro A, Zhang H, Clark KJ, Ekker SC, Wen X-Y. Disruption of pdgfra alters endocardial and myocardial fusion during zebrafish cardiac assembly. Biol Open. 2017;6(3):348–57.

    CAS  PubMed  PubMed Central  Google Scholar 

  102. Santiago FS, Khachigian LM. Ets-1 stimulates platelet-derived growth factor A-chain gene transcription and vascular smooth muscle cell growth via cooperative interactions with Sp1. Circ Res. 2004;95:479–87.

    Article  CAS  PubMed  Google Scholar 

  103. Spengler ML, Brattain MG. Sumoylation inhibits cleavage of Sp1 N-terminal negative regulatory domain and inhibits Sp1-dependent transcription. J Biol Chem. 2006;281:5567–74.

    Article  CAS  PubMed  Google Scholar 

  104. Milanini-Mongiat J, Pouysségur J, Pages G. Identification of two Sp1 phosphorylation sites for p42/p44 mitogen-activated protein kinases. J Biol Chem. 2002;277:20631–9.

    Article  CAS  PubMed  Google Scholar 

  105. Ataliotis P, Symes K, Chou MM, Ho L, Mercola M. PDGF signalling is required for gastrulation of Xenopus laevis. Development. 1995;121:3099–110.

    Article  CAS  PubMed  Google Scholar 

  106. Damm EW, Winklbauer R. PDGF-A controls mesoderm cell orientation and radial intercalation during Xenopus gastrulation. Development. 2011;138:565–75.

    Article  CAS  PubMed  Google Scholar 

  107. Nagel M, Tahinci E, Symes K, Winklbauer R. Guidance of mesoderm cell migration in the Xenopus gastrula requires PDGF signaling. Development. 2004;131:2727–36.

    Article  CAS  PubMed  Google Scholar 

  108. Andrae J, Gallini R, Betsholtz C. Role of platelet-derived growth factors in physiology and medicine. Genes Dev. 2008;22:1276–312.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  109. Soriano P. The PDGF alpha receptor is required for neural crest cell development and for normal patterning of the somites. Development. 1997;124:2691–700.

    Article  CAS  PubMed  Google Scholar 

  110. Moncho-Amor V, Ibanez de Caceres I, Bandres E, Martinez-Poveda B, Orgaz JL, Sanchez-Perez I, Zazo S, Rovira A, Albanell J, Jimenez B, et al. DUSP1/MKP1 promotes angiogenesis, invasion and metastasis in non-small-cell lung cancer. Oncogene. 2011;30:668–78.

    Article  CAS  PubMed  Google Scholar 

  111. Luehders K, Sasai N, Davaapil H, Kurosawa-Yoshida M, Hiura H, Brah T, Ohnuma S. The small leucine-rich repeat secreted protein Asporin induces eyes in Xenopus embryos through the IGF signalling pathway. Development. 2015;142:3351–61.

    Article  CAS  PubMed  Google Scholar 

  112. Pera EM, Wessely O, Li SY, De Robertis EM. Neural and head induction by insulin-like growth factor signals. Dev Cell. 2001;1:655–65.

    Article  CAS  PubMed  Google Scholar 

  113. Nentwich O, Dingwell KS, Nordheim A, Smith JC. Downstream of FGF during mesoderm formation in Xenopus: the roles of Elk-1 and Egr-1. Dev Biol. 2009;336:313–26.

    Article  CAS  PubMed  Google Scholar 

  114. Lee A, Oh JG, Gorski PA, Hajjar RJ, Kho C. Post-translational modifications in heart failure: small changes, big impact. Heart Lung Circ. 2016;25:319–24.

    Article  PubMed  Google Scholar 

  115. Mendler L, Braun T, Muller S. The ubiquitin-like SUMO system and heart function: from development to disease. Circ Research. 2016;118:132–44.

    Article  CAS  Google Scholar 

  116. Ueno N, Greene ND. Planar cell polarity genes and neural tube closure. Birth Defects Res C. 2003;69:318–24.

    Article  CAS  Google Scholar 

  117. Wallingford JB. Planar cell polarity, ciliogenesis and neural tube defects. Hum Mol Genet. 2006;15:R227–34.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

We thank Dr. Olivia Cox for assistance with the bioinformatics analysis and members of the Notre Dame Genomics & Bioinformatics Core Facility. We are grateful to Dr. Susanna Chiocca for providing a clone of Gam1. We acknowledge the expertise acquired from the organizers and instructors of the Cold Spring Harbor Laboratory Xenopus course. Constructive comments on the manuscript came from Olivia Cox, Joseph O’Tousa, and Holly Goodson.

Funding

This work was supported by the Faculty Scholarship Award Program (Office of Research, University of Notre Dame), providing funds for collection, analysis, and interpretation of data. The Genomics Pilot Projects Program (Eck Institute for Global Health, University of Notre Dame) provided funds for microarray experiments. These sources of support were not involved in the study design, analysis, interpretation of the results, or preparation of the manuscript.

Availability of data and materials

All data files have been uploaded to the Gene Expression Omnibus repository, including, experiment protocols, sample specifications, original CEL files, and processed Log2 GC-RMA signals under the following accession number GSE116164. Plasmids used in this study are available upon request from the corresponding author.

Author information

Authors and Affiliations

  1. Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana, 46556, USA

    Michelle M. Bertke, Kyle M. Dubiak & Paul W. Huber

  2. Department of Computer Science and Engineering, University of Notre Dame, Notre Dame, Indiana, 46556, USA

    Laura Cronin & Erliang Zeng

  3. Harper Cancer Research Institute, University of Notre Dame, Notre Dame, Indiana, USA

    Paul W. Huber

  4. Center for Stem Cells and Regenerative Medicine, University of Notre Dame, Notre Dame, Indiana, 46556, USA

    Paul W. Huber

  5. Present Address: College of Computer, Mathematical, and Natural Sciences, University of Maryland, College Park, MD, 20742, USA

    Michelle M. Bertke

  6. Present Address: Division of Biostatistics and Computational Biology, Iowa Institute for Oral Health Research, University of Iowa, Iowa City, IA, 52242, USA

    Erliang Zeng

  7. Present Address: Department of Preventive & Community Dentistry, College of Dentistry, University of Iowa, Iowa City, IA, 52242, USA

    Erliang Zeng

  8. Present Address: Department of Biostatistics, University of Iowa, Iowa City, IA, 52242, USA

    Erliang Zeng

  9. Present Address: Department of Biomedical Engineering, University of Iowa, Iowa City, IA, 52242, USA

    Erliang Zeng

Authors
  1. Michelle M. Bertke
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Kyle M. Dubiak
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Laura Cronin
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Erliang Zeng
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Paul W. Huber
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

MMB and PWH conceived and designed the project. MMB performed the experiments. KMD performed the Gam1 western blot. MMB, LC, EZ, and PWH carried out the bioinformatics analysis. MMB, EZ, and PWH wrote the manuscript with input from all authors. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Paul W. Huber.

Ethics declarations

Ethics approval

All animal experiments were performed in accordance with protocols approved by the University of Notre Dame Institutional Animal Care and Use Committee (Protocol Number 18-02-4408).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Heatmap showing the hierarchical clustering of all differentially expressed genes (DEGs) identified from comparisons between control (water injected) and Gam1 injected embryos at the three developmental stages (early gastrula, late gastrula, and early neurula, respectively). Numbers refer to the individual biological triplicates. (XLSX 323 kb)

Additional file 2:

Table S1. Gene expression data matrix associated with the heatmap as shown in Additional file 1: Figure S1. Genes and expression values in the matrix are listed in the same order as they appear in heatmap. (XLSX 34 kb)

Additional file 3:

Figure S2. Comparison of microarray data with measurement by qRT-PCR. The fold change between Gam1 and control embryos for the indicated gene as measured by microarray and qRT-PCR are compared with each pair of bars corresponding to early gastrula, late gastrula, and early neurula, respectively. Asterisks indicate samples in which the RNA level was too low to measure by qRT-PCR. (DOCX 157 kb)

Additional file 4:

Figure S3. Volcano plots of differentially expressed genes between Gam1 and control embryos. Fold change (log2) is plotted versus statistical significance at (A) early gastrula, (B) late gastrula, and (C) early neurula. Significance values of p < 0.05 (black line) indicates the cut off for accepted differential expression. (DOCX 22 kb)

Additional file 5:

Table S2. Gene Lists Constituting the Biological Processes Most Affected by Loss of SUMOylation Activity. (TIF 66126 kb)

Additional file 6:

Table S3. Biological Processes Associated with Transcription Factors Identified in Network Building. (TIF 21375 kb)

Additional file 7:

Table S4. SUMO Targets Sites and SUMO Interaction Motifs (SIM) in Top Transcription Factors from Network Building. (TIF 21809 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bertke, M.M., Dubiak, K.M., Cronin, L. et al. A deficiency in SUMOylation activity disrupts multiple pathways leading to neural tube and heart defects in Xenopus embryos. BMC Genomics 20, 386 (2019). https://doi.org/10.1186/s12864-019-5773-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-019-5773-3

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us