Gon4l regulates notochord boundary formation and cell polarity underlying axis extension by repressing adhesion genes

Williams, Margot L. K.; Sawada, Atsushi; Budine, Terin; Yin, Chunyue; Gontarz, Paul; Solnica-Krezel, Lilianna

doi:10.1038/s41467-018-03715-w

Download PDF

Article
Open access
Published: 03 April 2018

Gon4l regulates notochord boundary formation and cell polarity underlying axis extension by repressing adhesion genes

Margot L. K. Williams¹,
Atsushi Sawada^1,2,
Terin Budine¹,
Chunyue Yin²^nAff3,
Paul Gontarz¹ &
…
Lilianna Solnica-Krezel^1,2

Nature Communications volume 9, Article number: 1319 (2018) Cite this article

3139 Accesses
19 Citations
3 Altmetric
Metrics details

Subjects

Abstract

Anteroposterior (AP) axis extension during gastrulation requires embryonic patterning and morphogenesis to be spatiotemporally coordinated, but the underlying genetic mechanisms remain poorly understood. Here we define a role for the conserved chromatin factor Gon4l, encoded by ugly duckling (udu), in coordinating tissue patterning and axis extension during zebrafish gastrulation through direct positive and negative regulation of gene expression. Although identified as a recessive enhancer of impaired axis extension in planar cell polarity (PCP) mutants, udu functions in a genetically independent, partially overlapping fashion with PCP signaling to regulate mediolateral cell polarity underlying axis extension in part by promoting notochord boundary formation. Gon4l limits expression of the cell–cell and cell–matrix adhesion molecules EpCAM and Integrinα3b, excesses of which perturb the notochord boundary via tension-dependent and -independent mechanisms, respectively. By promoting formation of this AP-aligned boundary and associated cell polarity, Gon4l cooperates with PCP signaling to coordinate morphogenesis along the AP embryonic axis.

Simultaneous single-cell three-dimensional genome and gene expression profiling uncovers dynamic enhancer connectivity underlying olfactory receptor choice

Article Open access 15 April 2024

Evolution of tissue-specific expression of ancestral genes across vertebrates and insects

Article 15 April 2024

Spatially organized cellular communities form the developing human heart

Article Open access 13 March 2024

Introduction

Gastrulation is a critical period of animal development during which the three primordial germ layers—ectoderm, mesoderm, and endoderm—are specified, patterned, and shaped into a rudimentary body plan. During vertebrate gastrulation, an elongated anteroposterior (AP) axis emerges as the result of convergence and extension (C&E), a conserved set of gastrulation movements characterized by the concomitant AP elongation and mediolateral (ML) narrowing of the germ layers^1,2. C&E is accomplished by a combination of polarized cell behaviors, including directed migration and ML intercalation behavior (MIB)^3,4,5. During MIB, cells elongate and align their bodies and protrusions in the ML dimension and intercalate preferentially between their anterior and posterior neighbors⁵. This polarization of cell behaviors with respect to the AP axis is regulated by planar cell polarity (PCP) and other signaling pathways^6,7,8,9,10. Because these pathways are essential for MIB and C&E but do not affect cell fates^7,10,11, other mechanisms must spatiotemporally coordinate morphogenesis with embryonic patterning to ensure normal development. BMP, for example, coordinates dorsal–ventral axis patterning with morphogenetic movements by limiting expression of PCP signaling components and C&E to the embryo’s dorsal side¹². In general, though, molecular mechanisms that coordinate gastrulation cell behaviors with axial patterning are poorly understood, and remain one of the key outstanding questions in developmental biology.

Epigenetic regulators offer a potential mechanism by which broad networks of embryonic patterning and morphogenesis genes can be co-regulated. Epigenetic modifiers form protein complexes with chromatin factors that are thought to regulate their binding at specific genomic regions in context-specific ways¹³. The identities, functions, and specificities of chromatin factors with roles during embryogenesis are only now being elucidated, and some have described roles in cell fate specification and embryonic patterning¹⁴. One such chromatin factor is Gon4l, whose homologs have conserved roles in cell cycle regulation and/or embryonic patterning in plants, worms, flies, mice, and fish^{15,16,17,18,19}. However, the contribution of Gon4l or any other chromatin factor to morphogenesis is particularly poorly understood.

Here, we demonstrate a role for zebrafish Gon4l, encoded by ugly duckling (udu), as a regulator of embryonic axis extension during gastrulation. udu was identified in a forward genetic screen for enhancers of short axis phenotypes in PCP mutants, but we find it functions in parallel to PCP signaling. Instead, complete maternal and zygotic udu (MZudu) deficiency produces a distinct set of morphogenetic and cell polarity phenotypes that implicate the notochord boundary in ML cell polarity and cell intercalation during C&E. Extension defects in MZudu mutants are remarkably specific, as internalization, epiboly, and convergence gastrulation movements occur normally. Gene expression profiling reveals that Gon4l regulates expression of a large portion of the zebrafish genome, including genes with roles in housekeeping, patterning, and morphogenesis. Furthermore, Gon4l-associated genomic loci are identified by DNA adenine methyltransferase (Dam) identification^20,21 paired with high-throughput sequencing (DamID-seq), revealing both positive and negative regulation of putative direct targets by Gon4l. Mechanistically, we find that increased expression of epcam and itga3b, direct targets of Gon4l-dependent repression during gastrulation, each contribute to notochord boundary defects in MZudu mutants via a distinct molecular mechanism. This report thereby defines a critical role for a chromatin factor in the regulation of gastrulation cell behaviors in vertebrate embryos: by ensuring proper formation of the AP-aligned notochord boundary and associated ML cell polarity, Gon4l cooperates with PCP signaling to coordinate morphogenesis that extends the AP embryonic axis.

Results

Gon4l regulates axis extension during zebrafish gastrulation

To identify previously unknown regulators of C&E, we performed a three generation synthetic mutant screen^22,23 using zebrafish carrying a hypomorphic allele of the PCP gene knypek (kny)/glypican 4, kny^m818⁸. F0 wild-type (WT) males were mutagenized with N-ethyl-N-nitrosourea (ENU) and outcrossed to WT females. The resulting F1 fish were outcrossed to kny^m818/818 males (rescued by kny RNA injection) to generate F2 families whose F3 offspring were screened at 12 and 24 h post-fertilization (hpf) for short axis phenotypes (Fig. 1a). Screening nearly 100 F2 families yielded eight recessive mutations that enhanced axis extension defects in kny^m818/818 embryos (Fig. 1b–e). vu68 was found to be a L227P kny allele, and vu64 was a Y219* nonsense allele of the core PCP gene trilobite(tri)/vangl2, mutations in which exacerbate kny mutant phenotypes²⁴, demonstrating effectiveness of our screening strategy. We focused on vu66/vu66 mutants, which displayed a pleiotropic phenotype at 24 hpf, including shortened AP axes, reduced tail fins, and heart edema (Fig. 1d, Supplementary Fig.1b).

Employing the simple sequence repeat mapping strategy²⁵, we mapped the vu66 mutation to a small region on Chromosome 16 that contains the udu gene (Fig. 1f). udu encodes a conserved chromatin factor homologous to gon-4 in Caenorhabditis elegans¹⁶ and the closely related Gon4l in mammals^18,19. The previously described udu mutant phenotypes resemble those of homozygous vu66 embryos, including a shorter body axis and abnormal tail fins^18,26, making udu an excellent candidate for this PCP enhancer. Sequencing cDNA of the udu coding region from 24 hpf vu66/vu66 embryos revealed a T to A transversion at position 2261 predicted to change 753Y to a premature STOP codon (Fig. 1g). Furthermore, vu66 failed to complement a known udu^sq1 allele¹⁸ (Supplementary Fig. 1c). Together, these data establish vu66 as an udu allele and identify it as a recessive enhancer of axis extension defects in kny PCP mutant gastrulae.

Complete loss of udu function impairs axial extension

To assess the full role of Gon4l during early development, we eliminated maternal expression of udu¹⁸ using germline replacement²⁷. The resulting WT females carrying udu^vu66/vu66 (udu−/−) germline were crossed to udu^vu66/vu66 germline males to produce 100% embryos lacking both maternal (M) and zygotic (Z) udu function, hereafter referred to as MZudu mutants. These mutants appeared outwardly to develop normally until mid-gastrula stages (Fig. 1h, i), but exhibited clear abnormalities at the onset of segmentation (Fig. 1j, k). MZudu−/− gastrulae specified the three germ layers (Fig. 1l–m, Supplementary Fig. 2), formed an embryonic shield marked by gsc expression (Supplementary Fig. 2d), and completed epiboly on schedule (Fig. 1h–i), but somites were largely absent in mutants (Fig. 1j, k, n, o). Although myoD expression was observed within adaxial cells by whole mount in situ hybridization (WISH), its expression was not detected in nascent somites (Fig. 1o), similar to reported descriptions of Zudu mutants²⁸. Formation of adaxial cells is consistent with normal expression of their inducer shh in the axial mesoderm²⁹ (Supplementary Fig. 2h), and ntla/brachyury expression in the axial mesoderm was also largely intact (Supplementary Fig. 2j). Importantly, MZudu mutants were markedly shorter than age-matched WT controls throughout segmentation (Fig. 1j, k) and at 24 hpf (Supplementary Fig.3e, h). Although increased cell death was observed in MZudu mutants (as in Zudu mutants²⁸), inhibiting apoptosis via injection of RNA encoding the anti-apoptotic mitochondrial protein Bcl-xL³⁰ did not suppress their short axis phenotype (Supplementary Fig. 3a-f) as udu-gfp RNA did (Supplementary Fig. 3g-i). These results demonstrate a specific role for Gon4l in axial extension during gastrulation.

Gon4l regulates formation of the notochord boundary

Time-lapse Nomarski (Fig. 2a, b) and confocal (Fig. 2c, d) microscopy of dorsal mesoderm in MZudu−/− gastrulae revealed reduced definition and regularity of the boundary between axial and paraxial mesoderm, hereafter referred to as the notochord boundary, compared to WT (Fig. 2a, b, arrowheads). The ratio of the total/net length of notochord boundaries was significantly higher in MZudu−/− than in WT gastrulae at all-time points (Fig. 2e, two-way ANOVA p < 0.0001), indicative of decreased straightness. Interestingly, laminin was detected by immunostaining at the notochord boundary of both WT and MZudu−/− embryos (Fig. 2f, g), indicating that MZudu mutants form a bona fide boundary, albeit an irregular one. These results demonstrate that Gon4l is necessary for proper formation of the notochord boundary during gastrulation.

Cell polarity and intercalation are reduced in MZudu mutants

In vertebrate gastrulae, C&E is achieved chiefly through ML intercalation of polarized cells that elongate and align their cell bodies with the ML embryonic axis^3,5,7,9. To determine whether cell polarity defects underlie reduced axis extension in MZudu−/− gastrulae, we measured cell orientation (the angle of a cell’s long axis with respect to the embryo’s ML axis) and cell body elongation (length-to-width or aspect ratio (AR)) in confocal time-lapse series of fluorescent membrane-labeled WT (Fig. 3a–d) and MZudu−/− (Fig. 3e–h) gastrulae. Given the irregular notochord boundaries observed in MZudu−/− gastrulae (Fig. 2d, e), we examined the time course of cell polarization according to a cell’s position with respect to the boundary, i.e., boundary-adjacent “edge” cells versus those one or two cell diameters away (hereafter −1 and −2, respectively), and so on (see Fig. 3). Most WT axial mesoderm cells at midgastrulation (80% epiboly) were largely ML oriented and somewhat elongated, with boundary-adjacent “edge” cells being the least well oriented (median angle = 24.6°) (Fig. 3a–d). However, at the end of the gastrula period 90 min later, edge cells became highly aligned and elongated (median angle = 13.6°) similar to internal cell rows (Fig. 3a–d). All MZudu−/− axial mesoderm cells exhibited significantly reduced ML alignment (Fig. 3g, Kolmogorov–Smirnov test p < 0.0001), and elongation (Fig. 3h, Mann–Whitney test p < 0.0001) at 80% epiboly, but 90 min later only the edge cells remained less aligned than WT (median angle = 17.0°) (Fig. 3g), although AR of the edge and −1 cells remained reduced (Fig. 3h). These results indicate significantly reduced ML orientation of axial mesoderm cells in MZudu−/− gastrulae, a defect that persisted only in boundary-adjacent cells at late gastrulation. Importantly, this reduction in ML cell polarity was accompanied by significantly fewer cell intercalation events within the axial mesoderm of MZudu mutants compared to WT (Fig. 3i–k, T-test p < 0.05). As ML intercalation is the key cellular behavior required for vertebrate C&E^3,5, we conclude that this is likely the primary cause of axial mesoderm extension defects in MZudu mutants. We also observed significantly fewer mitoses in MZudu−/− gastrulae (Fig. 3l–n, T-tests), consistent with reports of Zudu mutants²⁸. Because decreased cell proliferation and the resulting reduced number of axial mesoderm cells were demonstrated to impair extension in zebrafish³¹, this could also be a contributing factor. Together, reduction of ML cell polarity, ML cell intercalations, and cell proliferation provide a suite of mechanisms resulting in impaired axial extension in MZudu−/− gastrulae. Combined with irregular notochord boundaries observed in MZudu mutants, we hypothesize that this boundary provides a ML orientation cue that contributes to ML polarization of axial mesoderm cells, and that this cue is absent or reduced in MZudu−/− gastrulae.

Gon4l functions independently of PCP signaling

Molecular control of ML cell polarity underlying C&E movements in vertebrate embryos is chiefly attributed to PCP and Gα12/13 signaling^6,7,8,9,10. While zygotic loss of udu enhanced axial extension defects in kny^m818/m818 PCP mutants (Fig. 1), it was unclear whether udu functions within or parallel to the PCP network. To address this, we generated compound MZudu;Zkny^fr6/fr6 (a nonsense/null kny allele⁸) mutants utilizing germline replacement as described above. Strikingly, these compound mutant embryos were substantially shorter than single MZudu or kny^fr6/fr6 mutants (also strict maternal udu mutants, which have no obvious phenotype) (Fig. 4a–d). Likewise, interference with vangl2/tri function in MZudu mutants by injection of MO1-vangl2 antisense morpholino oligonucleotide (MO)³² also exacerbated axis extension defects of MZudu mutants (Supplementary Fig. 4a-d). That reduced levels of PCP components kny or tri enhanced short axis phenotypes resulting from complete udu deficiency provides evidence that Gon4l affects axial extension via a parallel pathway. Furthermore, expression domains of genes encoding Wnt/PCP signaling components kny, tri, and wnt5 in MZudu−/− gastrulae were comparable to WT (Supplementary Fig. 4e-l). Finally, we found that the asymmetric intracellular localization of Prickle (Pk)-GFP, a core PCP component, to anterior cell membranes in WT gastrulae during C&E movements^33,34 was not affected in MZudu−/− gastrulae (Fig. 4e, f). This is consistent with intact PCP signaling in MZudu mutant gastrulae, and provides further evidence that Gon4l functions largely in parallel to PCP to regulate ML cell polarity and axis extension during gastrulation.

Gon4l-dependent boundary cue cooperates with PCP signaling

In addition to PCP signaling, notochord boundaries are required for proper C&E of the axial mesoderm in ascidian embryos³⁵ and are involved in the polarization of intercalating cells during Xenopus gastrulation¹. Boundary defects observed in MZudu−/− gastrulae are not a common feature of mutants with reduced C&E, however, as notochord boundaries in kny^fr6/fr6 gastrulae were straighter than in WT (Fig. 4g). Consistent with cell polarity defects previously reported in kny mutants⁸, all axial mesoderm cells failed to align ML within kny−/− embryos at 80% epiboly regardless of their position relative to the notochord boundary (Fig. 4h–k). After 90 min, however, kny−/− cells in the edge (and to a lesser extent, −1) position attained distinct ML orientation (median angle = 22.1°)(Fig. 4j). Indeed, the nearer a kny−/− cell was to the notochord boundary, the more ML aligned it was likely to be. This suggests that the notochord boundary provides a ML orientation cue that is independent of PCP signaling and functions across approximately two cell diameters. Furthermore, this boundary-associated cue appears to operate later in gastrulation, whereas PCP-dependent cell polarization is evident by 80% epiboly. Importantly, distinct cell polarity phenotypes observed in MZudu−/− and kny−/− gastrulae provide further evidence that Gon4l functions in parallel to PCP signaling.

To assess how PCP signaling interacts with the proposed Gon4l-dependent boundary cue, we examined the polarity of axial mesoderm cells in compound zygotic kny−/−;udu−/− mutant gastrulae (Fig. 4l–o). As observed in single kny−/− mutant gastrulae⁸ (Fig. 4h–k), both ML orientation and elongation of all axial mesoderm cells were reduced in double mutants at 80% epiboly, regardless of position with respect to the notochord boundary (Fig. 4n, o). By late gastrulation, however, edge cells in double mutant gastrulae failed to attain the ML alignment observed in kny−/− mutants, but instead remained largely randomly oriented (median angle = 38.0°) (Fig. 4n). This exacerbation of cell orientation defects correlated with stronger axis extension defects in compound kny−/−;udu−/− compared to single kny−/− mutants (Fig. 1e). This supports our hypothesis that a Gon4l-dependent boundary cue regulates ML alignment of axial mesoderm cells independent of PCP signaling, and that these two mechanisms partially overlap in time and space to cooperatively polarize all axial mesoderm cells and promote axial extension.

Loss of Gon4l results in large-scale gene expression changes

As a nuclear-localized chromatin factor^18,36 (Supplementary Fig. 3j), Gon4l is unlikely to influence morphogenesis directly. To identify genes regulated by Gon4l with potential roles in morphogenesis, we performed RNA sequencing (RNA-seq) in MZudu−/− and WT tailbud-stage embryos. Analysis of relative transcript levels revealed that more than 11% of the genome was differentially expressed (p_adj < 0.05, ≥twofold change) in MZudu mutants (Fig. 5a). Of these ~2950 differentially expressed genes, 1692 exhibited increased expression in MZudu mutants compared to 1259 with decreased expression (Supplementary data 1). Functional annotation analysis revealed that genes downregulated in MZudu−/− gastrulae were enriched for ontology terms related to chromatin structure, transcription, and translation (Fig. 5b, c), while upregulated genes were enriched for terms related to biosynthesis, metabolism, and protein modifications (Fig. 5b, d). Notably, many of these misregulated genes are considered to have “housekeeping” functions, in that they are essential for cell survival. We further examined genes with plausible roles in morphogenesis, including those encoding signaling and adhesion molecules, and found the majority of genes in both classes were expressed at higher levels in MZudu−/− gastrulae (Fig. 5e, f). This putative increase in adhesion was of particular interest given the tissue boundary defects observed in MZudu mutants.

DamID-seq identifies putative direct targets of Gon4l

To determine genomic loci with which Gon4l protein associates, and thereby distinguish direct from indirect targets of Gon4l regulation, we employed DNA adenine methyltransferase (Dam) identification paired with next generation sequencing (DamID-seq)^20,21. To this end, we generated Gon4l fused to Escherichia coli Dam, which methylates adenine residues within genomic regions in its close proximity²⁰. Small equimolar amounts of RNA encoding a Myc-tagged Gon4l-Dam fusion or a Myc-tagged GFP-Dam control were injected into one-celled embryos, and genomic DNA was collected at tailbud stage (see Methods section). In support of this Gon4l-Dam fusion being functional, it localized to the nuclei of zebrafish embryos and partially rescued MZudu−/− embryonic phenotypes (Supplementary Fig. 5a-k). Methylated, and therefore Gon4l-proximal, genomic regions were then selected using methylation specific restriction enzymes and adaptors, amplified to produce libraries, and sequenced. Because Dam is highly active, even an untethered version methylates DNA within open chromatin regions²⁰, and so libraries generated from embryos expressing GFP-Dam served as controls.

Unbiased genome-wide analysis of DamID reads revealed a significant enrichment of Gon4l-Dam over GFP-Dam in promoter regions, and a significant underrepresentation of Gon4l within intergenic regions (Fig. 6a, T-test p < 0.05, Fig. 6b). Although no global difference was detected within gene bodies (Fig. 6a), examination of individual loci revealed ~4500 genes and over 2300 promoters in which at least one region was highly Gon4l-enriched (P_adj ≤ 0.01, ≥fourfold enrichment over GFP) (Fig. 6c–g, Supplementary Data 2, 3). Of these, ~1000 genes were co-enriched for Gon4l in both the promoter and gene body (Fig. 6c, g). Levels of Gon4l association across a gene and its promoter were significantly correlated (Supplementary Fig. 5l, Spearmann correlation p < 0.0001), indicating co-enrichment or co-depletion for Gon4l at both gene features of many loci. Within gene bodies, we found robust enrichment specifically within 5′ untranslated regions (UTRs) (Fig. 6b), consistent with association of Gon4l at or near transcription start sites (Fig. 6d). Of the ~2950 genes differentially expressed in MZudu−/− gastrulae, approximately 28% (812) were also enriched for Gon4l at the gene body (492), promoter (170), or both (150) (Fig. 6g), and will hereafter be described as putative direct Gon4l targets. histh1, for example, was among the most downregulated genes by RNA-seq in MZudu mutants (Supplementary Data 1) and was highly enriched for Gon4l at both its promoter and gene body (Fig. 6c). By contrast, tbx6 expression was also reduced in MZudu mutants (Supplementary Data 1), but exhibited no enrichment of Gon4l over GFP controls (Fig. 6f). Approximately 35% and 53% of genes enriched for Gon4l in only the gene body or only the promoter, respectively, were positively regulated (i.e., downregulated in MZudu mutants), as were 50% of genes co-enriched at both features (Fig. 6h). Furthermore, among these positively regulated genes, differential expression levels (the degree to which WT expression exceeded MZudu−/−) correlated positively and significantly with Gon4l-enrichment levels across both gene bodies and promoters (Fig. 6h, Spearmann correlation p < 0.01). A similar correlation was not observed for negatively regulated genes, hence, the highest levels of Gon4l enrichment were associated with positive regulation of gene expression. These results implicate Gon4l as both a positive and negative regulator of gene expression during zebrafish gastrulation.

Gon4l limits itga3b to promote boundary straightness

We next examined our list of putative direct Gon4l target genes for those with potential roles in tissue boundary formation and/or cell polarity. epcam, which encodes Epithelial cell adhesion molecule (EpCAM), stood out because it was not only enriched for Gon4l by DamID (Fig. 7a) and upregulated in MZudu−/− gastrulae (Fig. 5f, Supplementary Data 1), but was also identified in a Xenopus overexpression screen for molecules that disrupt tissue boundaries³⁷. Furthermore, EpCAM negatively regulates non-muscle myosin activity³⁸, making it a compelling candidate. We also chose to examine itga3b, which encodes Integrinα3b, because as a component of a laminin receptor³⁹, it is an obvious candidate for a molecule involved in formation of a tissue boundary at which laminin is highly enriched (Fig. 2f). DamID revealed a region within the itga3b promoter at which Gon4l was highly enriched (Fig. 7b, c), and itga3b expression was increased in MZudu−/− gastrulae by RNA-seq (Supplementary Data 1), which was also validated (along with epcam) by qRT-PCR (Fig. 7d, e). We found that overexpression of either epcam or itga3b by RNA injection into WT embryos recapitulated the irregular notochord boundaries (Fig. 7f, g) and reduced ML orientation and elongation of axial mesoderm cells (Supplementary Figs. 6d-g, 7c-f) seen in MZudu−/− gastrulae. Conversely, while injection of an epcam MO (MO2-epcam⁴⁰, Supplementary Fig. 7b) had no effect on MZudu−/− boundaries (Fig. 7h), injection of a translation-blocking itga3b MO (MO1-itga3b⁴¹, Supplementary Fig. 6b) significantly improved boundary straightness in MZudu mutants compared with control-injected siblings (Fig. 7i, two-way ANOVA p < 0.0001). This indicates that excess itga3b is largely responsible for reduced straightness of notochord boundaries in MZudu−/− gastrulae. Furthermore, disruption of WT notochord boundaries by excess itga3b was abrogated by co-injection with a MO against lama5 (MO3-lama5⁴², Supplementary Fig. 6c) (Fig. 7j), which encodes its ligand Lamininα5⁴¹, indicating that this effect is ligand dependent. Both itga3b and lama5 are expressed within the axial mesoderm of WT zebrafish gastrulae^43,44. Moreover, although an Itgα3b-GFP fusion was not enriched at WT notochord boundaries, it became increasingly localized to the plasma membrane of axial mesoderm cells (but not presomitic mesoderm or neuroectodermal cells) as gastrulation proceeded (Fig. 7k–m). This implies that intracellular localization (and likely function) of Integrinα3b is regulated in a stage- and tissue-dependent manner, as was reported for other Integrins during Xenopus gastrulation⁴⁵. Finally, to determine whether excess itga3b was sufficient to phenocopy the kny−/−;udu−/− double mutant phenotype identified in our synthetic screen (Fig. 1), we overexpressed itga3b in kny−/− embryos and found that it exacerbated the short axis phenotype of these PCP mutants (Fig. 7n–p), an effect not produced by injection of control GFP RNA. Together, these results indicate that negative regulation of itga3b expression by Gon4l is essential for proper notochord boundary formation and axis extension in zebrafish gastrulae. Interestingly, ML cell polarity defects were not suppressed in MZudu−/− embryos injected with itga3b (or epcam) MO compared to control-injected mutant siblings (Supplementary Figs. 6h-o, 7g-n), despite improved boundary straightness. This implies that other Gon4l-dependent boundary properties are involved, and/or that loss of udu affects ML polarity cell-autonomously, for example, by making cells unable to respond to the boundary-associated polarity cue.

Loss of Gon4l and excess epcam reduce boundary tension

In WT zebrafish and Xenopus embryos, notochord boundaries straighten over time (Fig. 2e) and accumulate myosin⁴⁶, implying that they are under tension. Because boundaries of MZudu mutants are less straight than those of WT gastrulae, we hypothesized that tension at these boundaries is reduced. To test this, we laser-ablated interfaces between axial mesoderm cells in live WT and MZudu−/− gastrulae and recorded recoil of adjacent cell vertices as a measure of tissue tension⁴⁷ (Fig. 8a, b). Cell interfaces were classified according to an established convention as V junctions (actively shrinking to promote cell intercalation), T junctions (not shrinking)^47,48 or Edge junctions (falling on and thus comprising the notochord boundary). We found that recoil distances at WT Edge junctions were significantly greater than that of WT V (Kruskal–Wallis test p = 0.0001) or T junctions (Kruskal–Wallis test p = 0.0026), demonstrating that the notochord boundary is under greater tension than the rest of the tissue (Fig. 8c–e). Consistent with our hypothesis, recoil distance was significantly smaller in MZudu−/− than in WT gastrulae for all classes of junctions (Fig. 8c–e, two-way ANOVA p < 0.01), and this decrease was largest and most significant in Edge junctions (Fig. 8c, two-way ANOVA p < 0.0001).

Because excess epcam and itga3b were sufficient to reduce boundary straightness (Fig. 7), we next tested whether epcam and itga3b overexpression could likewise affect boundary tension. Injection of WT embryos with epcam RNA was sufficient to lower tissue tension at the notochord boundary and throughout the axial mesoderm (Fig. 8c–e). Curiously, although we found itga3b to be largely causative of reduced boundary straightness in MZudu mutants (Fig. 7), its overexpression did not produce a similar reduction in tension (Fig. 8c–e), implying a tension-independent role for Itgα3b in boundary straightness. Accordingly, treatment of itga3b-overexpressing WT embryos with Calyculin A, a myosin phosphatase inhibitor that increases myosin contractility⁴⁹ and notochord boundary tension (Supplementary Fig. 8a), did not restore their boundary straightness despite increased tension (Fig. 8f, Supplementary Fig. 8c). By contrast, Calyculin A treatment significantly increased boundary tension and straightness in epcam-overexpressing WT embryos (Fig. 8g, Supplementary Fig. 8d, two-way ANOVA), indicating a tension-dependent role for EpCAM in boundary formation. Consistent with its previously described inhibition of myosin activity³⁸, these data support a role for EpCAM as a negative regulator of myosin-dependent tissue tension at the notochord boundary. Notably, treatment of MZudu mutants with Calyculin A did not restore boundary straightness or axial mesoderm cell polarity (Supplementary Fig. 8b, e-h), demonstrating that increasing tension alone was not sufficient to improve these phenotypes. Given that the effect of EpCAM on the boundary is largely tension dependent, this likely explains why the epcam MO failed to improve boundary straightness in MZudu mutants (Fig. 7h). Together these results implicate excess EpCAM and Integrinα3b as key molecular defects underlying reduced boundary tension and reduced boundary straightness, respectively, observed in MZudu−/− gastrulae. We propose a model whereby Gon4l limits expression of these adhesion molecules to ensure proper formation of the notochord boundary, which together with additional boundary-independent roles of Gon4l cooperates with PCP signaling to promote ML cell polarity underlying C&E gastrulation movements (Fig. 8h).

Discussion

Substantial advances have been made in defining signaling pathways that regulate gastrulation cell behaviors and shape the vertebrate body plan^6,7,9,10, but epigenetic control of these morphogenetic processes remains largely unexplored. Here we have described the conserved chromatin factor Gon4l as a regulator of polarized cell behaviors underlying axis extension during zebrafish gastrulation. We identified a large number of Gon4l target genes, including many with known or predicted roles in morphogenesis, and linked misregulation of a subset of these genes to specific morphogenetic defects. Because Gon4l does not bind DNA directly, we predict that Gon4l-enriched genomic loci are direct targets of chromatin modifying protein complexes with which Gon4l associates³⁶. Only a fraction of the thousands of Gon4l-enriched genes and promoters exhibited corresponding changes in gene expression during gastrulation, which may reflect one or more possible scenarios. Either Gon4l does not alter expression of all loci with which it associates, loci recently occupied by Gon4l do not yet reflect changes in transcript levels, or because DamID provides a “history” of Gon4l association, both formerly and currently occupied loci are represented in our data. We also identified loci at which Gon4l was depleted compared to GFP-Dam controls (Supplementary Fig. 5), and speculate they represent open chromatin regions with which Gon4l does not associate. Surprisingly, although a larger number of putative Gon4l direct target genes were negatively regulated (i.e., upregulated in MZudu mutants), the highest levels of Gon4l enrichment were correlated with positive regulation of gene expression (Fig. 6). This argues against Gon4l acting strictly as a negative regulator of gene expression, a role assigned to it based on in vitro evidence and thought to be mediated by its interactions with Histone deacetylases³⁶. Our data instead indicate that Gon4l acts as both a positive and negative regulator of gene expression during zebrafish gastrulation, implying context-specific interactions with multiple epigenetic regulatory complexes.

Phenotypes caused by complete udu deficiency are conspicuously pleiotropic, but our studies point to remarkably specific roles for udu in gastrulation morphogenesis. Loss of udu function reduced tissue extension without affecting mesendoderm internalization, epiboly, prechordal plate migration, or convergence (see below). Moreover, the dorsal gastrula organizer and all three germ layers were specified (Fig. 1, Supplementary Fig. 2), indicating that MZudu mutants do not suffer from a general delay or arrest of development; rather Gon4l regulates a specific subset of gastrulation cell behaviors, including ML cell polarity and intercalation in the axial mesoderm. Importantly, the role of Gon4l in these processes is independent of PCP signaling as supported by several lines of evidence (Fig. 4, Supplementary Fig. 4). Additionally, whereas PCP mutants exhibited reduced convergence^7,8, evidenced by a larger number of cell rows in kny−/− axial mesoderm (Fig. 4), MZudu mutants contained a normal number of axial mesoderm cell rows (Fig. 3), indicating no obvious convergence defect. Despite these apparently parallel functions, RNA-seq and DamID-seq experiments revealed that Gon4l regulates expression of some PCP genes, including wnt11, wnt11r, prickle1a, prickle1b, celsr1b, celsr2, and fzd2 (Supplementary Data 1-3). However, functional redundancies within the PCP network^50,51 likely allow for intact PCP signaling in MZudu mutants despite misregulation of some PCP components. Notably, gastrula morphology and cell polarity defects in MZudu mutants were also distinct from ventralized and dorsalized patterning mutants with impaired C&E^11,12.

Modulation of adhesion at tissue boundaries has been implicated as a driving force of cell intercalation^46,52, and indeed genes annotated as encoding adhesion molecules tended to be expressed at higher levels in MZudu−/− gastrulae (Fig. 5). Two of these, itga3b and epcam, exhibited increased expression in MZudu−/− gastrulae by RNA-seq and qRT-PCR, were identified as putative direct Gon4l targets by DamID (Fig. 7), and were each sufficient to reduce notochord boundary straightness and ML cell polarity when overexpressed in WT embryos (Fig. 7, Supplementary Figs. 6, 7). Furthermore, decreasing levels of Integrinα3b (but not EpCAM) improved boundary straightness in MZudu mutants (Fig. 7), while overexpression of epcam (but not itga3b) was sufficient to reduce tissue tension at the notochord boundary of WT embryos (Fig. 8). In addition, increasing tension with Calyculin A in epcam-overexpressing (but not itga3b overexpressing) WT embryos normalized their notochord boundary straightness (Fig. 8). These observations support a model in which misregulation of each of these molecules contributes to distinct cellular defects in MZudu mutants: excess itga3b disrupts boundary straightness in a tension-independent manner via interactions with its laminin ligand, while EpCAM negatively regulates myosin-dependent boundary tension. EpCAM was similarly shown to negatively regulate non-muscle myosin contractility in Xenopus gastrulae³⁸, where experimental perturbation of myosin activity disrupts notochord boundary formation⁴⁶.

Links between myosin-dependent tissue tension, boundary formation, and cell intercalation are well established. In addition to Xenopus notochord boundaries⁴⁶, myosin accumulation increases tension at compartment boundaries in Drosophila⁵³, which was shown to bias cell intercalations⁵⁴. We found that decreased boundary straightness in epcam-overexpressing WT gastrulae could be rescued by pharmacologically increasing tissue tension, demonstrating a causal relationship between tension and boundary straightness (Fig. 8). However, increasing tissue tension was not sufficient to improve boundary straightness in MZudu−/− gastrulae (Supplementary Fig. 8), implying additional molecular bases for this phenotype, namely excess itga3b. Furthermore, restoration of boundary straightness (with itga3b MO) in MZudu−/− gastrulae was not sufficient to normalize ML cell orientation, raising questions as to the functional relationship between these elements in this context. We speculate that increased boundary straightness and tension together may be required to promote ML polarity of boundary-adjacent cells. Alternatively, improved boundary straightness in MZudu−/− gastrulae may have restored the polarity cue, but additional gene expression changes rendered mutant cells unable to respond to it. Although it is not yet clear whether and how boundary straightness and tension affect the polarity of adjacent cells, evidence from this study and others strongly implicate the notochord boundary in ML polarization of cells undergoing C&E^1,35. Namely, within kny−/− PCP mutants, only boundary-adjacent edge cells became ML aligned during late gastrulation. Together with the additive cell polarity defects observed in compound kny−/−;udu−/− mutants, this is strong genetic evidence that the boundary provides a Gon4l-dependent, PCP-independent cell polarity cue (Fig. 4).

In this study of zebrafish Gon4l, we have begun to dissect the logic of epigenetic regulation of gastrulation cell behaviors by revealing a key role for this chromatin factor in limiting expression of specific genes with specific morphogenetic consequences. This has important developmental implications, as C&E gastrulation movements are sensitive to both gain and loss of gene function^7,10. We propose that by negatively regulating Integrinα3b and EpCAM levels, Gon4l promotes development of the anteroposteriorly aligned notochord boundary and influences ML polarity and intercalation of axial mesoderm cells. In cooperation with PCP signaling, this cue coordinates ML cell polarity with embryonic patterning to drive AP embryonic axis extension (Fig. 8h).

Methods

Zebrafish strains and embryo staging

Adult zebrafish were raised and maintained according to established methods⁵⁵ in compliance with standards established by the Washington University Animal Care and Use Committee. Embryos were obtained from natural mating and staged according to morphology as described⁵⁶. All studies on WT were carried out in AB* or AB*/Tübingen backgrounds. Additional lines used include kny^m818, kny^fr6⁸, udu^sq1¹⁸, and udu^vu66 (this work). Embryos of these strains generated from heterozygous intercrosses were genotyped by PCR after completion of each experiment. Germline-replaced fish were generated by the method described in ref. ²⁷. Briefly, donor embryos from udu^vu66/+ intercrosses or females with udu^vu66/vu66 germline were injected with synthetic RNA encoding GFP-nos1-3′UTR⁵⁷, and WT host embryos were injected with a MO against dead end1 (MO1-dnd1)²⁷ to eliminate host germ cells. Cells were transplanted from the embryonic margin of donor blastulae to the embryonic margin of hosts at sphere stage, and both hosts and donors were cultured in agarose-coated plates. Host embryos were screened for GFP + germ cells at 36–48 hpf, and the genotype of corresponding donors was determined by phenotype. All putative udu^vu66/vu66 germline hosts were raised to adulthood and confirmed by crossing to udu^vu66/+ animals prior to use in experiments. Fish were chosen from their home tank to be crossed at random. The resulting embryos were also chosen from the dish at random for injection and inclusion in experiments. No fewer than six (but often many more) embryos of each condition or genotype were analyzed in at least three independent trials for each experiment in which measurements were made, except where noted. These rare cases reflect difficulty in obtaining sufficient MZudu−/− embryos due to limited availability and productivity of germline-replaced females.

Synthetic mutant screening

WT male fish were mutagenized by the chemical mutagen ENU as described in ref. ²². Briefly, adult male fish were incubated four to six times, for 1 h each, in a solution of 3.5 mM ENU at 21 °C. After treatment, fish were allowed to recover for 2–3 h in 10 mg/L 3-aminobenzoic acid ethyl ester (MESAB) at 17–19 °C and an additional day without MESAB before returning to the circulating water system. Approximately 3 weeks after mutagenesis, these males were outcrossed to WT females to produce F1 families. F2 families were obtained by crossing F1 fish with fish homozygous for the hypomorphic knypek allele kny^m818 (rescued by injection with synthetic kny WT RNA). F3 embryos obtained from F2 cross were screened by morphology at 12 and 24 hpf to identify recessive enhancers of the kny^m818/m818 short axis mutant phenotype⁸.

Positional cloning

We employed the positional cloning approach using a panel of CA simple sequence length polymorphism markers representing 25 linkage groups²⁵ to map the vu66 mutation to chromosome 16 between the markers of z17403 and z15431. Given phenotypic similarities between vu66/vu66 and the udu^sq1/sq1 mutant phenotype¹⁸, we sequenced udu cDNA from 24 hpf vu66/vu66 mutant embryos revealing a T2261A transversion that is predicted to create Y753STOP nonsense mutation. We designed a dCAPS⁵⁸ marker for the vu66 mutation and confirmed that no recombination occurred in 810 vu66/vu66 homozygous embryos.

Microinjection

One-celled embryos were aligned within agarose troughs generated using custom-made plastic molds and injected with 1–3 pL volumes using pulled glass needles. Synthetic mRNAs for injection were made by in vitro transcription from linearized plasmid DNA templates using Invitrogen mMessage mMachine kits. Doses of RNA per embryo were as follows: 100 pg membrane Cherry, 50 pg membrane GFP, 25 pg udu-gfp, 200 pg epcam, 200 pg itga3b, 50 pg itga3b-GFP, 1 pg gfp-dam-myc, 3 pg udu-dam-myc for DamID experiments, 20 pg gfp-dam-myc for MZudu−/− rescue, and 150 pg GFP-nos1-3′UTR for germline transplantation. To assess Pk-GFP localization, embryos were injected at one-cell stage with membrane Cherry RNA, then injected with 15 pg Drosophila prickle-GFP and 20 pg H2B-RFP RNAs into a single blastomere at 16-cell stage as described^33,34. Injections of MOs were carried out as for synthetic RNA. Doses of MOs per embryo were as follows: 3 ng MO1-dnd1²⁷, 4 ng MO1-tri/vangl2³², 1 ng MO2-epcam⁴⁰, 2 ng MO1-itga3b⁴¹, and 1 ng MO3-lama5⁴².

WISH

Antisense riboprobes were transcribed using NEB T7 or T3 RNA polymerase and labeled with digoxygenin (DIG) (Roche). WISH was performed according to ref. ⁵⁹. Briefly, embryos were fixed overnight in 4% paraformaldehyde (PFA) in phosphate buffered saline (PBS), rinsed in PBS+0.1% tween 20 (PBT), and dehydrated into methanol. Embryos were then rehydrated into PBT, incubated for at least 2 h in hybridization solution with 50% formamide (in 0.75 M sodium chloride, 75 mM sodium citrate, 0.1% tween 20, 50 μg/mL heparin (Sigma), and 200 μg/mL tRNA) at 70 °C, then hybridized overnight at 70 °C with antisense probes diluted approximately 1 ng/μl in hybridization solution. Embryos were washed gradually into 2X SSC buffer (0.3 M sodium chloride, 30 mM sodium citrate), and then gradually from SSC to PBT. Embryos were blocked at room temperature for several hours in PBT with 2% goat serum and 2 mg/mL bovine serum albumin (BSA), then incubated overnight at 4 °C with anti-DIG antibody (Roche #11093274910) at 1:5000 in block. Embryos were rinsed extensively in PBT, and then in staining buffer (PBT+100 mM Tris pH 9.5, 50 mM MgCl₂, and 100 mM NaCl) prior to staining with BM Purple solution (Roche).

Immunofluorescent staining

Embryos were fixed in 4% PFA, rinsed in PBT, digested briefly in 10 μg/mL proteinase K, refixed in 4% PFA, rinsed in PBT, and blocked in 2 mg/mL BSA + 2% goat serum in PBT. Embryos were then incubated overnight in rabbit anti-Laminin (Sigma L9393) at 1:200, mouse anti-Myc (Cell Signaling 2276) at 1:1000, or rabbit anti-phospho histone H3 (Upstate 06-570) at 1:500 in blocking solution, rinsed in PBT, and incubated overnight in Alexa Fluor 488 anti-Rabbit IgG, 568 anti-Rabbit, or 568 anti-Mouse (Invitrogen) at 1:1000 in PBT. Embryos were co-stained with 4′,6-Diamidino-2-Phenylindole, Dihydrochloride and rinsed in PBT prior to mounting in agarose for confocal imaging.

TUNEL staining

Terminal deoxynucleotidyl transferase (TdT) dUTP Nick-End Labeling (TUNEL) staining to detect apoptosis was carried out according to the instructions for the ApopTag Peroxidase in situ apoptosis detection kit (Millipore) with modifications. Briefly, embryos were fixed in 4% PFA, digested with 10 μg/mL proteinase K, refixed with 4% PFA, and post-fixed in chilled ethanol:acetic acid 2:1, rinsing in PBT after each step. Embryos were incubated overnight with TdT and rinsed in stop/wash buffer, then blocked, incubated with anti-DIG antibody, and stained in Roche BM Purple staining solution.

Microscopy

Live embryos expressing fluorescent proteins or fixed embryos subjected to immunofluorescent staining were mounted in 0.75% low-melt agarose in glass bottomed 35-mm petri dishes for imaging using a modified Olympus IX81 inverted spinning disc confocal microscope equipped with Voltran and Cobolt steady-state lasers and a Hamamatsu ImagEM EM CCD digital camera. For live time-lapse series, 60 μm z-stacks with a 2 μm step were collected every three to ten minutes (depending on the experiment) for 3 h using a ×40 dry objective lens. Embryo temperature was maintained at 28.5 °C during imaging using a Live Cell Instrument stage heater. When necessary, embryos were extracted from agarose after imaging for genotyping. For immunostained embryos, 200 μm z-stacks with a 1 or 2 μm step were collected using a ×10 or ×20 dry objective lens, depending on the experiment. Bright field and transmitted light images of live embryos and in situ hybridizations were collected using a Nikon AZ100 macroscope.

Laser ablation tension measurements

Embryos were injected at one-cell stage with mCherry mRNA and mounted at 80% epiboly for imaging (as described above) on a Zeiss 880 Airyscan 2-photon inverted confocal microscope. An infrared laser tuned to 710–730 nm was used to ablate fluorescently labeled cell interfaces, immediately followed by a quick time-lapse series of ten images with no interval to record recoil after each ablation event. Images were collected using a ×40 water immersion objective lens, and a Zeiss stage heater was used to maintain embryo temperature at 28.5 °C. Approximately 8–12 cell interfaces were ablated per embryo. Image series were analyzed using ImageJ to determine the inter-vertex distance of each cell interface prior to and immediately after ablation and used to calculate recoil distance.

Calyculin A treatment

Embryos were mounted in agarose as described above, then embryo medium containing 50 nM Calyculin A (Sigma-Aldrich) was added approximately 30 min prior to the start of imaging and remained throughout the imaging period.

Image analysis

ImageJ was used to visualize and manipulate all microscopy datasets. For immunostained embryos, multiple z-planes were projected together to visualize the entire region of interest. For live embryo analysis, a single z-plane through the length of the axial mesoderm was chosen for each time point. When possible, embryo images were analyzed prior to genotyping, and all images were coded during analysis. To measure cell orientation and elongation, the AP axis in all embryo images was aligned prior to manual outlining of cells. A fit ellipse was used to measure orientation of each cell’s major axis and its AR. The TissueAnalyzer ImageJ package⁶⁰ was used to automatically segment time-lapse series of axial mesoderm and detect T1 transitions. Boundary straightness was measured by manually tracing the notochord boundary to determine total length, then dividing it by the length of a straight line connecting the ends of the boundary (net length). To assess Pk-GFP localization, isolated cells expressing Pk-GFP were scored according to subcellular localization of GFP signal.

Quantitative RT-PCR

Total RNA was isolated from tailbud stage WT and MZudu−/− embryos homogenized in Trizol (Life Technologies), 1 μg of which was used to synthesize cDNA using the iScript kit (BioRad) following manufacturer’s protocol. SYBR green (BioRad) qRT-PCR reactions were run in a CFX Connect Real-Time PCR detection system (BioRad) in technical triplicate. Primers used are as follows:

epcam: F-TGAGGACGGGGATTGAGAAC

R-GAGCCTGCCATCCTTGTCAT

itga3b: F-CCGGTGTTGGGAGAAGAGAC

R-CTTGAAGAAACCACACGAAGGG

EF1a: F-CTGGAGGCCAGCTCAAACAT

R-ATCAAGAAGAGTAGTACCGCTAGCATTAC

Rpl13a: F-TCTGGAGGACTGTAAGAGGTATGC

R-AGACGCACAATCTTGAGAGCAG

RNA sequencing and analysis

RNA for sequencing was isolated from 50 WT or MZudu−/− embryos per sample at tailbud stage according to instructions for the Dynabeads mRNA direct kit (Ambion). Embryos from two clutches per genotype (e.g., WT A and B) were collected, then divided in two (e.g., WT A1, A2, B1, and B2) to yield four independently prepared libraries representing two biological and two technical replicates per genotype. Libraries for were prepared according to instructions for the Epicentre ScriptSeq v2 RNA-seq Library preparation kit (Illumina). Briefly, RNA was enzymatically fragmented prior to cDNA synthesis. cDNA was then 3′ tagged, purified using Agencourt AMPure beads, and PCR amplified, at which time sequencing indexes were added. Indexed libraries were then purified and submitted to the Washington University Genome Technology Access Center for sequencing using an Illumina HiSeq 2500 to obtain single-ended 50 bp reads. Raw reads were mapped to the zebrafish GRCz10 reference genome using STAR (2.4.2a)⁶¹ with default parameters. FeatureCounts (v1.4.6) from the Subread package⁶² was used to quantify the number of uniquely mapped reads (phred score ≥10) to gene features based on the Ensembl annotations (v83). Significantly differentially expressed genes were determined by using DESeq2 in the negative binomial distribution model⁶³ with a cutoff of adjusted p-value ≤0.05 and fold change ≥2.0. Heatmaps were built using the heatmaps2 package in R, and other plots were built using the ggplot2 package in R.

DamID-seq

E. coli DNA adenine methyltransferase (EcoDam) was cloned from the pIND(V5) EcoDam plasmid²¹ (a kind gift from Dr. Bas Van Steensel, Netherlands Cancer Institute) by Gibson assembly⁶⁴ into a 3′ Gateway entry vector containing 6 Myc tags and an SV40 polyA signal (p3E-MTpA). The resulting p3E-EcoDam-MTpA vector was then Gateway cloned into PCS2 + downstream of udu cDNA or eGFP to produce C terminal fusions. The resulting udu-dam-myc and gfp-dam-myc plasmids were linearized by KpnI digestion and transcribed using the mMessage mMachine SP6 in vitro transcription kit (Ambion). WT AB* embryos were injected at one-cell stage with 3 pg mRNA encoding Gon4l-Dam-Myc or with 1 pg encoding GFP-Dam-Myc as controls. The remainder of the protocol was carried out largely as in ref. ²¹ with modifications. Genomic (g)DNA was collected at tailbud stage using a Qiagen DNeasy kit with the addition of RNAse A. gDNA was digested with DpnI overnight, followed by ligation of DamID adaptors. Un-methylated regions were destroyed by digestion with DpnII, and then methylated regions were amplified using primers complementary to DamID adaptors. Two identical 50 μl PCR reactions were performed and pooled for each sample to reduce amplification bias, and three independent biological replicates were collected per condition. Amplicons were purified using a Qiagen PCR purification kit, then digested with DpnII to remove DamID adaptors. Finally, samples were purified using Agencourt AMPure XP beads and submitted for library preparation and sequencing at the Washington University GTAC using an Illumina HiSeq 2500 to obtain single-ended 50 bp reads.

Analysis of DamID-seq data

Raw reads were aligned to zebrafish genome GRCz10 by using bwa mem (v0.7.12) with default parameters⁶⁵, then sorted and converted into bam format by using SAMtools (v1.2)⁶⁶. The zebrafish genome was divided into continuous 1000 bp bins, and FeatureCounts (v1.4.6) from the Subread package⁶² was used to quantify the number of uniquely mapped reads (phred score ≥10) in each bin. Significantly differentially Gon4l-associated bins were determined by using DESeq2 in the negative binomial distribution model⁶³ with stringent cutoff: adjusted p-value ≤0.01 and fold change ≥4.0. Promoter regions were defined as 2 kb upstream of the transcription start site of a gene based on Ensembl annotations (v83). Wiggle and bigwig files were created from bam files using IGVtools (v2.3.60) (https://software.broadinstitute.org/software/igv/igvtools) and wigToBigWig (v4)⁶⁷, respectively. BigWig tracks were visualized using IGV (v2.3.52)⁶⁸.

Statistical analysis

Graphpad Prism 6 and 7 software was used to perform statistical analyses and generate graphs of data collected from embryo images. The statistical tests used varied as appropriate for each experiment and are described in the text and figure legends. Data were tested for normal distribution, and nonparametric tests (Mann–Whitney and Kolmogorov–Smirnov) were used for all non-normally distributed data. Normally distributed data with similar variance between groups were analyzed using parametric tests (T-tests and ANOVAs). All tests used were two tailed. Differential expression and differential enrichment analysis of RNA and DamID sequencing data were completed as described above. Panther⁶⁹ was used to classify differentially expressed genes and to produce pie charts; DAVID Bioinformatic Resources⁷⁰ was used for functional annotation analysis.

Subcloning

The full-length udu open reading frame was subcloned from WT cDNA using following primers:

udu cacc F1: CACCATGGGATGGAAACGCAAGTCTTC

udu TGA R: TCAGTCCTGCTCTTCATCAGTGGC

udu R: GTCCTGCTCTTCATCAGTGGCCGAC

udu cDNAs with or without a stop codon were cloned into the pENTR/D-TOPO (Thermo Fisher) vector, which were then Gateway cloned into PCS2 + upstream of polyA signal or eGFP to produce a C terminal fusion.

The full-length epcam open reading frame was subcloned from WT cDNA using the following primers:

epcam F: GGATCCCATCGATTCGATGAAGGTTTTAGTTGCCTTG

epcam R: ACTCGAGAGGCCTTGTTAAGAAATTGTCTCCATCTC

The 5′ portion of the itga3b open reading frame was cloned from WT cDNA, and the 3′ portion was cloned from a partial cDNA clone (GE/Dharmacon) using the following primers:

5′ itga3b F: TTTGCAGGCGCGCCGGATCCCATCGATTCGATGGCCGGAAAGTCTCTG

5′ itga3b R: ATTTGAGTGAGTATGGAATGGAGATGTTGAGCG

3′ itga3b F: TCAACATCTCCATTCCATACTCACTCAAATACTCAGG

3′ itga3b R: GTTCTAGAGGTTTAAACTCGAGAGGCCTTGTCAGAACTCCTCCGTCAG

The resulting amplicons were Gibson cloned⁶⁴ into PJS2 (a derivative of PCS2+) linearized with EcoRI. To create an Itgα3b-GFP fusion, the itga3b open reading frame minus the stop codon was amplified from this plasmid and Gibson cloned upstream of an eGFP open reading frame into PJS2.

Data availability

The authors declare that all data supporting the findings of this study are available within the article and its Supplementary Information files or from the corresponding author upon reasonable request.

Processed RNA-seq and DamID-seq data are available in Supplementary Data files 1-3 of this publication, and raw data have been deposited in the Gene Expression Omnibus (GEO) database under the accession codes:

GSE96575 (for RNA-seq) (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=ipwzmaoqjpajbqd&acc=GSE96575) and GSE96576 (for DamID-seq) (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=apkloooyvxsjxij&acc=GSE96576).

References

Shih, J. & Keller, R. Patterns of cell motility in the organizer and dorsal mesoderm of Xenopus laevis. Development 116, 915–930 (1992).
CAS PubMed Google Scholar
Warga, R. M. & Kimmel, C. B. Cell movements during epiboly and gastrulation in zebrafish. Development 108, 569–580 (1990).
CAS PubMed Google Scholar
Shih, J. & Keller, R. Cell motility driving mediolateral intercalation in explants of Xenopus laevis. Development 116, 901–914 (1992).
CAS PubMed Google Scholar
Sepich, D. S., Calmelet, C., Kiskowski, M. & Solnica-Krezel, L. Initiation of convergence and extension movements of lateral mesoderm during zebrafish gastrulation. Dev. Dyn. 234, 279–292 (2005).
Article CAS PubMed Google Scholar
Keller, R. et al. Mechanisms of convergence and extension by cell intercalation. Philos. Trans. R. Soc. Lond. B Biol. Sci. 355, 897–922 (2000).
Article CAS PubMed Central PubMed Google Scholar
Wallingford, J. B. et al. Dishevelled controls cell polarity during Xenopus gastrulation. Nature 405, 81–85 (2000).
Article ADS CAS PubMed Google Scholar
Jessen, J. R. et al. Zebrafish trilobite identifies new roles for Strabismus in gastrulation and neuronal movements. Nat. Cell Biol. 4, 610–615 (2002).
Article CAS PubMed Central PubMed Google Scholar
Topczewski, J. et al. The zebrafish glypican knypek controls cell polarity during gastrulation movements of convergent extension. Dev. Cell 1, 251–264 (2001).
Article CAS PubMed Google Scholar
Williams, M., Yen, W., Lu, X. & Sutherland, A. Distinct apical and basolateral mechanisms drive planar cell polarity-dependent convergent extension of the mouse neural plate. Dev. Cell 29, 34–46 (2014).
Article CAS PubMed Central PubMed Google Scholar
Lin, F. et al. Essential roles of Gα_12/13 signaling in distinct cell behaviors driving zebrafish convergence and extension gastrulation movements. J. Cell Biol. 169, 777–787 (2005).
Article CAS PubMed Central PubMed Google Scholar
Solnica-Krezel, L. et al. Mutations affecting cell fates and cellular rearrangements during gastrulation in zebrafish. Development 123, 67–80 (1996).
CAS PubMed Google Scholar
Myers, D. C., Sepich, D. S. & Solnica-Krezel, L. Bmp activity gradient regulates convergent extension during zebrafish gastrulation. Dev. Biol. 243, 81–98 (2002).
Article CAS PubMed Google Scholar
Chen, T. & Dent, S. Y. Chromatin modifiers and remodellers: regulators of cellular differentiation. Nat. Rev. Genet. 15, 93–106 (2014).
Article CAS PubMed Google Scholar
Ma, Y. et al. The chromatin remodeling protein Bptf promotes posterior neuroectodermal fate by enhancing Smad2-activated wnt8a expression. J. Neurosci. 35, 8493-8506 (2015).
Article CAS PubMed Google Scholar
Brownfield, L. et al. Arabidopsis DUO POLLEN3 is a key regulator of male germline development and embryogenesis. Plant Cell 21, 1940–1956 (2009).
Article CAS PubMed Central PubMed Google Scholar
Friedman, L., Santa Anna-Arriola, S., Hodgkin, J. & Kimble, J. gon-4, a cell lineage regulator required for gonadogenesis in Caenorhabditis elegans. Dev. Biol. 228, 350–362 (2000).
Article CAS PubMed Google Scholar
Bulchand, S., Menon, S. D., George, S. E. & Chia, W. Muscle wasted: a novel component of the Drosophila histone locus body required for muscle integrity. J. Cell Sci. 123, 2697-2707 (2010).
Article CAS PubMed Google Scholar
Liu, Y. et al. The zebrafish udu gene encodes a novel nuclear factor and is essential for primitive erythroid cell development. Blood 110, 99–106 (2007).
Article CAS PubMed Google Scholar
Lu, P. et al. The Justy mutation identifies Gon4-like as a gene that is essential for B lymphopoiesis. J. Exp. Med. 207, 1359–1367 (2010).
Article CAS PubMed Central PubMed Google Scholar
van Steensel, B. & Henikoff, S. Identification of in vivo DNA targets of chromatin proteins using tethered dam methyltransferase. Nat. Biotechnol. 18, 424–428 (2000).
Article PubMed Google Scholar
Vogel, M. J., Peric-Hupkes, D. & van Steensel, B. Detection of in vivo protein–DNA interactions using DamID in mammalian cells. Nat. Protoc. 2, 1467–1478 (2007).
Article CAS PubMed Google Scholar
Solnica-Krezel, L., Schier, A. F. & Driever, W. Efficient recovery of ENU-induced mutations from the zebrafish germline. Genetics 136, 1401–1420 (1994).
CAS PubMed Central PubMed Google Scholar
Mullins, M. C., Hammerschmidt, M., Haffter, P. & Nüsslein-Volhard, C. Large-scale mutagenesis in the zebrafish: in search of genes controlling development in a vertebrate. Curr. Biol. 4, 189–202 (1994).
Article CAS PubMed Google Scholar
Marlow, F. et al. Functional interactions of genes mediating convergent extension, knypek and trilobite, during the partitioning of the eye primordium in zebrafish. Dev. Biol. 203, 382–399 (1998).
Article CAS PubMed Google Scholar
Knapik, E. W. et al. A microsatellite genetic linkage map for zebrafish (Danio rerio). Nat. Genet. 18, 338–343 (1998).
Article CAS PubMed Google Scholar
Hammerschmidt, M. et al. Mutations affecting morphogenesis during gastrulation and tail formation in the zebrafish, Danio rerio. Development 123, 143–151 (1996).
CAS PubMed Google Scholar
Ciruna, B. et al. Production of maternal-zygotic mutant zebrafish by germ-line replacement. Proc. Natl Acad. Sci. USA 99, 14919-14924 (2002).
Article ADS CAS PubMed Central PubMed Google Scholar
Lim, C. H., Chong, S. W. & Jiang, Y. J. Udu deficiency activates DNA damage checkpoint. Mol. Biol. Cell 20, 4183-4193 (2009).
Article CAS PubMed Central PubMed Google Scholar
Blagden, C. S., Currie, P. D., Ingham, P. W. & Hughes, S. M. Notochord induction of zebrafish slow muscle mediated by Sonic hedgehog. Genes Dev. 11, 2163-2175 (1997).
Article CAS PubMed Central PubMed Google Scholar
Payne, E. M. et al. Ddx18 is essential for cell-cycle progression in zebrafish hematopoietic cells and is mutated in human AML. Blood 118, 903–915 (2011).
Article CAS PubMed Central PubMed Google Scholar
Liu, Y., Sepich, D. S. & Solnica-Krezel, L. Stat3/Cdc25a-dependent cell proliferation promotes embryonic axis extension during zebrafish gastrulation. PLoS Genet. 13, e1006564 (2017).
Article PubMed Central PubMed Google Scholar
Park, M. & Moon, R. T. The planar cell-polarity gene stbm regulates cell behaviour and cell fate in vertebrate embryos. Nat. Cell Biol. 4, 20–25 (2002).
Article CAS PubMed Google Scholar
Yin, C., Kiskowski, M., Pouille, P. A., Farge, E. & Solnica-Krezel, L. Cooperation of polarized cell intercalations drives convergence and extension of presomitic mesoderm during zebrafish gastrulation. J. Cell Biol. 180, 221–232 (2008).
Article CAS PubMed Central PubMed Google Scholar
Ciruna, B., Jenny, A., Lee, D., Mlodzik, M. & Schier, A. F. Planar cell polarity signalling couples cell division and morphogenesis during neurulation. Nature 439, 220–224 (2006).
Article ADS CAS PubMed Central PubMed Google Scholar
Veeman, M. T. et al. Chongmague reveals an essential role for laminin-mediated boundary formation in chordate convergence and extension movements. Development 135, 33–41 (2008).
Article CAS PubMed Google Scholar
Lu, P. et al. The developmental regulator protein Gon4l associates with protein YY1, co-repressor Sin3a, and histone deacetylase 1 and mediates transcriptional repression. J. Biol. Chem. 286, 18311-18319 (2011).
Article CAS PubMed Central PubMed Google Scholar
Maghzal, N., Vogt, E., Reintsch, W., Fraser, J. S. & Fagotto, F. The tumor-associated EpCAM regulates morphogenetic movements through intracellular signaling. J. Cell Biol. 191, 645–659 (2010).
Article CAS PubMed Central PubMed Google Scholar
Maghzal, N., Kayali, H. A., Rohani, N., Kajava, A. V. & Fagotto, F. EpCAM controls actomyosin contractility and cell adhesion by direct inhibition of PKC. Dev. Cell 27, 263–277 (2013).
Article CAS PubMed Google Scholar
Nishiuchi, R. et al. Characterization of the ligand-binding specificities of integrin alpha3beta1 and alpha6beta1 using a panel of purified laminin isoforms containing distinct alpha chains. J. Biochem. 134, 497–504 (2003).
Article CAS PubMed Google Scholar
Slanchev, K. et al. The epithelial cell adhesion molecule EpCAM is required for epithelial morphogenesis and integrity during zebrafish epiboly and skin development. PLoS Genet. 5, e1000563 (2009).
Article PubMed Central PubMed Google Scholar
Carney, T. J. et al. Genetic analysis of fin development in zebrafish identifies furin and hemicentin1 as potential novel fraser syndrome disease genes. PLoS Genet. 6, e1000907 (2010).
Article PubMed Central PubMed Google Scholar
Webb, A. E. et al. Laminin alpha5 is essential for the formation of the zebrafish fins. Dev. Biol. 311, 369–382 (2007).
Article CAS PubMed Google Scholar
Thisse, B. & Thisse, C. ZFIN Direct Data Submission http://zfin.org (2004).
Pollard, S. M. et al. Essential and overlapping roles for laminin alpha chains in notochord and blood vessel formation. Dev. Biol. 289, 64–76 (2006).
Article CAS PubMed Google Scholar
Ramos, J. W., Whittaker, C. A. & DeSimone, D. W. Integrin-dependent adhesive activity is spatially controlled by inductive signals at gastrulation. Development 122, 2873-2883 (1996).
CAS PubMed Google Scholar
Fagotto, F., Rohani, N., Touret, A. S. & Li, R. A molecular base for cell sorting at embryonic boundaries: contact inhibition of cadherin adhesion by ephrin/Eph-dependent contractility. Dev. Cell 27, 72–87 (2013).
Article CAS PubMed Google Scholar
Rauzi, M., Verant, P., Lecuit, T. & Lenne, P. F. Nature and anisotropy of cortical forces orienting Drosophila tissue morphogenesis. Nat. Cell Biol. 10, 1401–1410 (2008).
Article CAS PubMed Google Scholar
Shindo, A. & Wallingford, J. B. PCP and septins compartmentalize cortical actomyosin to direct collective cell movement. Science 343, 649–652 (2014).
Article ADS CAS PubMed Central PubMed Google Scholar
Suzuki, A. & Itoh, T. Effects of calyculin A on tension and myosin phosphorylation in skinned smooth muscle of the rabbit mesenteric artery. Br. J. Pharmacol. 109, 703–712 (1993).
Article CAS PubMed Central PubMed Google Scholar
Kilian, B. et al. The role of Ppt/Wnt5 in regulating cell shape and movement during zebrafish gastrulation. Mech. Dev. 120, 467–476 (2003).
Article CAS PubMed Google Scholar
Nikaido, M., Law, E. W. & Kelsh, R. N. A systematic survey of expression and function of zebrafish frizzled genes. PLoS ONE 8, e54833 (2013).
Article ADS CAS PubMed Central PubMed Google Scholar
Green, J. B., Dominguez, I. & Davidson, L. A. Self-organization of vertebrate mesoderm based on simple boundary conditions. Dev. Dyn. 231, 576–581 (2004).
Article PubMed Google Scholar
Monier, B., Pélissier-Monier, A., Brand, A. H. & Sanson, B. An actomyosin-based barrier inhibits cell mixing at compartmental boundaries in Drosophila embryos. Nat. Cell Biol. 12, 60–65 (2010). sup pp. 61–69.
Article CAS PubMed Google Scholar
Umetsu, D. et al. Local increases in mechanical tension shape compartment boundaries by biasing cell intercalations. Curr. Biol. 24, 1798–1805 (2014).
Article CAS PubMed Google Scholar
Westerfield, M. The Zebrafish Book. A Guide for the Laboratory Use of Zebrafish (Danio rerio) (University of Oregon Press, Eugene, OR, 1993).
Kimmel, C. B., Ballard, W. W., Kimmel, S. R., Ullmann, B. & Schilling, T. F. Stages of embryonic development of the zebrafish. Dev. Dyn. 203, 253–310 (1995).
Article CAS PubMed Google Scholar
Köprunner, M., Thisse, C., Thisse, B. & Raz, E. A zebrafish nanos-related gene is essential for the development of primordial germ cells. Genes Dev. 15, 2877-2885 (2001).
PubMed Central PubMed Google Scholar
Neff, M. M., Neff, J. D., Chory, J. & Pepper, A. E. dCAPS, a simple technique for the genetic analysis of single nucleotide polymorphisms: experimental applications in Arabidopsis thaliana genetics. Plant J. 14, 387–392 (1998).
Article CAS PubMed Google Scholar
Thisse, C. & Thisse, B. High-resolution in situ hybridization to whole-mount zebrafish embryos. Nat. Protoc. 3, 59–69 (2008).
Article CAS PubMed Google Scholar
Aigouy, B., Umetsu, D. & Eaton, S. Segmentation and quantitative analysis of epithelial tissues. Methods Mol. Biol. 1478, 227–239 (2016).
Article CAS PubMed Google Scholar
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
Article CAS PubMed Google Scholar
Liao, Y., Smyth, G. K. & Shi, W. The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 41, e108 (2013).
Article PubMed Central PubMed Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Article PubMed Central PubMed Google Scholar
Gibson, D. G. et al. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat. Methods 6, 343–345 (2009).
Article CAS PubMed Google Scholar
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760 (2009).
Article CAS PubMed Central PubMed Google Scholar
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
Article PubMed Central PubMed Google Scholar
Kent, W. J., Zweig, A. S., Barber, G., Hinrichs, A. S. & Karolchik, D. BigWig and BigBed: enabling browsing of large distributed datasets. Bioinformatics 26, 2204-2207 (2010).
Article CAS PubMed Central PubMed Google Scholar
Thorvaldsdóttir, H., Robinson, J. T. & Mesirov, J. P. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief. Bioinform 14, 178–192 (2013).
Article PubMed Google Scholar
Mi, H., Muruganujan, A., Casagrande, J. T. & Thomas, P. D. Large-scale gene function analysis with the PANTHER classification system. Nat. Protoc. 8, 1551–1566 (2013).
Article PubMed Google Scholar
Huang, da W., Sherman, B. T. & Lempicki, R. A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57 (2009).
Article CAS Google Scholar

Download references

Acknowledgements

We thank Dr. Bas Van Steensel for EcoDam plasmids, Drs. Christine and Bernard Thisse for WISH probes, Dr. Bo Zhang and Dr. Scott Higdon for bioinformatics help, Dr. Matthew Hass for DamID advice, Bisiayo Fashemi for assistance with image analysis, the Washington University Center for Cellular Imaging for assistance with laser ablation experiments, and the Washington University Genome Technology Access Center for library preparation and sequencing services. The National Institutes of Health grants R01GM55101 and R35GM118179 to L.S.-K. and F32GM113396 to M.L.K.W., and a W.M. Keck Foundation Fellowship to M.L.K.W. in part supported this study.

Author information

Chunyue Yin
Present address: Division of Pediatric Gastroenterology, Hepatology, and Nutrition, Cincinnati Children’s Hospital, Cincinnati, OH, 45229, USA

Authors and Affiliations

Department of Developmental Biology, Washington University School of Medicine, Saint Louis, MO, 63110, USA
Margot L. K. Williams, Atsushi Sawada, Terin Budine, Paul Gontarz & Lilianna Solnica-Krezel
Department of Biological Sciences, Vanderbilt University, Nashville, TN, 37235, USA
Atsushi Sawada, Chunyue Yin & Lilianna Solnica-Krezel

Authors

Margot L. K. Williams
View author publications
You can also search for this author in PubMed Google Scholar
Atsushi Sawada
View author publications
You can also search for this author in PubMed Google Scholar
Terin Budine
View author publications
You can also search for this author in PubMed Google Scholar
Chunyue Yin
View author publications
You can also search for this author in PubMed Google Scholar
Paul Gontarz
View author publications
You can also search for this author in PubMed Google Scholar
Lilianna Solnica-Krezel
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

M.L.K.W., A.S., and L.S.-K. designed the study. M.L.K.W., A.S., and T.B. performed experiments. C.Y. participated in the initial forward genetic screen. P.G. performed bioinformatic analysis. M.L.K.W. and L.S.-K. wrote the manuscript.

Corresponding author

Correspondence to Lilianna Solnica-Krezel.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Supplementary Information(PDF 22557 kb)

Peer Review File(PDF 2005 kb)

Description of Additional Supplementary Files(DOCX 15 kb)

Supplementary Data 1(XLSX 323 kb)

Supplementary Data 2(XLSX 180 kb)

Supplementary Data 3(XLSX 115 kb)

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Williams, M.L.K., Sawada, A., Budine, T. et al. Gon4l regulates notochord boundary formation and cell polarity underlying axis extension by repressing adhesion genes. Nat Commun 9, 1319 (2018). https://doi.org/10.1038/s41467-018-03715-w

Download citation

Received: 13 June 2017
Accepted: 06 March 2018
Published: 03 April 2018
DOI: https://doi.org/10.1038/s41467-018-03715-w

This article is cited by

Newly identified Gon4l/Udu-interacting proteins implicate novel functions
- Su-Mei Tsai
- Kuo-Chang Chu
- Yun-Jin Jiang
Scientific Reports (2020)
Pcdh18a regulates endocytosis of E-cadherin during axial mesoderm development in zebrafish
- Bernadett Bosze
- Yosuke Ono
- Steffen Scholpp
Histochemistry and Cell Biology (2020)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.