Introduction

The midline is an essential embryonic structure during vertebrate embryogenesis. It provides structural support of vertebrates, and plays a role in tissue patterning and in the establishment of the left right asymmetry (L-R) during internal organogenesis1,2,3. The development of the midline/body axis occurs from the onset of gastrulation. The process is initially marked by mesoderm and endoderm progenitor cell compaction at the dorsal side of the gastrula (dorsal organizer shield) which then undergoes simultaneous morphological changes, cellular intercalation and migration (convergent extension, CE) that is governed by specific signalling pathways to form the body axis4,5. Failure in this process and the associated signalling pathways alter midline formation and result in developmental abnormalities including impaired elongation of the body length6,7,8, establishment of L-R asymmetry9 and tissue patterning10. Many key signalling molecules and signalling pathways regulate midline formation during embryogenesis, including Platelet-Derived Growth Factor (PDGF)11,12, the non-canonical Wnt/Planar Cell Polarity (PCP) pathway13, Bone Morphogenetic Protein (BMPs)14 and Fibroblast Growth Factor (FGFs)15. However, the role of phospholipid signalling in this context remains poorly understood.

Lysophosphatidic acid (LPA) is a unique bioactive lysophospholipid that induces pleiotropic effects through binding to its specific G-protein coupled receptors lpa1–6 (reviewed in refs 16 and 17) in various cell types. LPA signalling affects proliferation, cell survival, motility, morphological rearrangements and differentiation. Dysregulation of LPA signalling is associated with early and late embryonic abnormalities, suggesting critical roles for LPA during development18,19,20,21. Studies in LPA receptor knockout mouse models have shown a prominent role of LPA signalling in particular in neural and vascular development19,22. Lpa1(−/−) mice die postnatally due to an impaired suckling behavior likely associated with a defective olfaction and impaired CNS development19. Lpa2(−/−) mice do not display severe phenotypes although some intracellular signalling pathways such as PLC activation, Ca2+ mobilization, and stress fiber formation are altered23. Lpa3(−/−) female mice display a reproductive impairment, with delayed embryo implantation and embryo spacing alterations24 while Lpa4(−/−) mice have defects in blood vessel formation leading to haemorrhages in many organs at different embryonic stages22. Lpa5(−/−) mice show deficit in response to neuropathic pain (reviewed in ref. 25) and the loss of lpa6 in Xenopus embryo disrupts forebrain development26.

The ectonucleotide pyrophosphatase/phosphodiesterase 2 (enpp2), also known as autotaxin (atx) or lysophospholipase D, is the main enzyme responsible for the synthesis of extracellular LPA27,28. The loss of enpp2 in mice is lethal between E9.5–E10.5 because of severe vascular defects in the yolk sac and embryo25,29,30, hampering the characterization of the LPA function in the early mouse embryo. In the early zebrafish embryo, the enpp2/lpa3 axis regulates the L-R internal organ asymmetry but the molecular mechanisms responsible are unknown31. Furthermore, overexpression of enpp2 in zebrafish induces cardia bifida32. Enpp2/LPA signalling also regulates oligodendrocyte production and differentiation of the zebrafish hindbrain33. Although these studies suggest essential roles for LPA signalling during embryogenesis, the understanding of LPA functions remains limited.

To determine the function of LPA during embryogenesis, we transiently overexpressed the LPA-producing enzyme enpp2 in developing zebrafish. Our gain of function study of enpp2 in the zebrafish embryo shows for the first time the unique role of enpp2 in regulating the cell migration during gastrulation and subsequent formation of the axial midline, as well as the establishment of L-R asymmetry. Furthermore, rescue experiments show that these effects are lpa1–3-dependent and mediated by the Rho/ROCK pathway.

Results

Overexpression of enpp2 alters axis formation in the zebrafish embryo and modulates expression of the midline axis genes shha and ntl

We cloned the full-length enpp2 (NM_200603.1)34 to generate probes for whole mount in situ hybridization (WISH). WISH analysis of the developing zebrafish embryos showed that enpp2 mRNA expression is dynamically regulated during development, which was also confirmed by quantitative RT-PCR analysis (Fig. 1A, Suppl. Fig. 1). Enpp2 was maternally deposited (sphere) and expressed at a very low level in the yolk syncytial layer (YSL) during the early gastrulation stage (50% epiboly and shield) (Fig. 1A, Suppl. Fig. 1A,B). At the end of gastrulation (tail bud), enpp2 was notably expressed in the midline axis and its levels continued to increase during the segmentation period (10 somite stage, 15 somite stage) (Fig. 1A, Suppl. Fig. 1), suggesting a potential role in midline and axis formation. In order to assess the role of enpp2 during embryogenesis, we injected increasing amounts of capped enpp2 mRNA into the zebrafish embryo (1–4 cell stage), determined the morphology (Fig. 1B) and measured the phenotype penetrance (Fig. 1C). The overexpression of enpp2 resulted in significant axis defects and in a kinked notochord in a dose-dependent manner (Fig. 1B–F). This was accompanied by aberrant somite shapes, highlighted by the lack of chevron-shaped somites and shortened body length (Fig. 1B). Embryos injected with enpp2 mRNA at increasing concentrations exhibited dose-dependent penetrance of phenotypes (25 pg: 19.3 ± 6.8%; 50 pg: 25.6 ± 9.6%; 100 pg: 60.7 ± 2.7%; 200 pg: 60.5 ± 4.3%, Fig. 1D). These phenotypes suggest midline and axis defects during embryogenesis. To examine this further, we performed WISH to assess the expression patterns of midline markers shha and ntl to mark the notochord of the zebrafish embryo at 10 somite stage and 24 hours post fertilization, respectively (Fig. 1E,F). The majority of enpp2 -overexpressing embryos displayed abnormal pattern of shha and ntl expression, characterized by kinked, patchy or expanded patterns, or expression in multiple buds (duplicated), which indicated a notochord defect. This finding suggests a role of enpp2 in regulating the midline formation and its impact on the expression of the midline axis genes shha and ntl.

Figure 1
figure 1

Enpp2 overexpression alters midline axis formation in the early embryogenesis of zebrafish.

Developmental series of WISH at designated stages were performed using enpp2 antisense riboprobe. (A) Lateral view and animal/dorsal view of WISH early zebrafish embryos at designated stages. (B) Representative pictures of control, mild (slight delay and smooth somite borders) and severe (developmental delay and midline axial defect) phenotypes in enpp2 injected embryos. (C,D) Quantification of phenotype variant and penetrance (detectable phenotype) following injection with different doses of enpp2 RNA and normalized to control. The sample size (n) is stated as numerical value above each bar. Data are mean ± SEM from at least four independent experiments. Statistical analysis was established by one-way ANOVA; *P < 0.05; ***P < 0.001. Dead embryos were normalised to uninjected embryos (E) Representative WISH pictures of control and enpp2 injected embryos with midline axis gene probes ntl or shha riboprobes (F). >50 embryos used in each experiment. Scale bars: 200 μm.

Enpp2 overexpression induces L-R patterning defect

Establishment of the midline is important for L-R asymmetry and loss of enpp2/LPA signaling has been linked to asymmetry phenotypes31,35,36,37. To establish if increased enpp2/LPA signaling alters the L-R patterning of the zebrafish embryo, we examined the expression of the nodal-related asymmetric genes, lefty1/2 and southpaw (spaw) in the lateral plate mesoderm (LPM) during somite stages in control wild type and enpp2-overexpressing zebrafish. As shown in Fig. 2A–D, lefty1/2 and spaw were mainly expressed on the left side of the LPM during mid-late somitogenesis (15–21ss) in the control wild type zebrafish, which is consistent with previous reports38,39. However, following enpp2 overexpression, the zebrafish displayed a more randomized expression pattern of these genes, either right-sided, absent, or bilaterally expressed in the LPM compared to control (Fig. 2A–D). In addition, the spaw expression domain was shifted posteriorly in some of the enpp2 overexpression embryos suggesting that migration of the lateral plate mesoderm is delayed or perturbed (Fig. 2A–D). Taken together, this demonstrated that enpp2 overexpression modulates L-R patterning and maintenance of the expression of the nodal-related asymmetric genes lefty1/2 and spaw.

Figure 2
figure 2

Enpp2 modulates nodal-related asymmetry gene expressions, KV formation and ciliogenesis.

Representative images of WISH of the control and enpp2-injected zebrafish embryos at 15–21ss with lefty2 (A) and spaw (C) antisense riboprobes respectively. The midline/notochord is highlighted by ntl expression. The determination of left, bilateral, absence, and right groups was based on the location of lefty2 or spaw expression compared to ntl (midline). The quantification of lefty2 and spaw expression distribution embryos is shown in (B) and (D) respectively. (B) The proportion of lefty2 expression on the left, bilateral, absence and right side are 67.3%, 4.5%, 28.2%, and 0% respectively in the control embryos and 54%, 12.6%, 28.7%, and 4.6% respectively in the enpp2 overexpressed embryos. (D) The proportion of spaw expression on the left, bilateral, absence and right side are 86%, 4%, 8%, and 2% respectively in the control embryos and 61.3%, 21.8%, 12.7%, 4.2% respectively in the enpp2-overexpressed embryos. (E) Representative images of the KV morphology in control- and in enpp2 overexpressing embryos (5–8ss). (F) Quantification KV size in enpp2 overexpressed embryos compared to control. (G) Whole-mount antibody staining of zona occludens 1 (ZO-1) that labels cell junctions and outlines the KV in control enpp2 injected embryos. (H) Representative confocal Z-stack images of cilia in the KV, marked with arl13b GFP, highlighting cilia morphology and distribution in control- and enpp2 overexpressing embryos at 5–8ss. (I) Quantification of the cilia length between control- and enpp2 overexpressing embryos. The sample size (n) is stated as numerical value above each bar. Data are mean ± SEM Statistical analysis was established by t-test; ***P < 0.001. n indicated above bars. White circles highlight the KV lumen. (J) Enpp2 overexpression altered the expression pattern of foxj1, a master regulator of ciliogenesis, observed at 90% epiboly embryo.

Enpp2 modulates the organogenesis and ciliogenesis of the KV

In the zebrafish, the Kupffer’s vesicle (KV) plays an important role in the establishment of L-R patterning by maintaining the flow of nodal morphogens along the body axis40. Ciliated cells in the KV drive gradients of signalling molecules that shape the L-R asymmetry. We next investigated if the formation and morphology of the KV were altered by enpp2 overexpression. Firstly, we examined the lumen shape of the KV at the 6–8-somite stage, when the KV is formed and nodal morphogen flow occurs40. In control embryos, the KV appeared as a flattened sphere with a single round lumen, while in the overexpressed-enpp2 fish, the lumen was smaller, misshapen or absent (Fig. 2E–G). Quantification of KV size showed that the lumen size was dramatically reduced in the enpp2 overexpressed embryos (45.00 ± 7.1%) compared to control (18.46 ± 6.80%) (Fig. 2F). To examine KV cilia, we performed live imaging using the axonemal marker arl13b-GFP41. Live imaging showed that the cilia lengths in the KV were significantly reduced following enpp2 overexpression (2.00 ± 0.04 units), compared to control (2.80 ± 1.44, p < 0.001, n > 200, Fig. 2H,I). This data indicates that enpp2 regulates cilia formation. Foxj1 is a master regulator of ciliogenesis and labels cells in the KV42,43. Analysis of foxj1 expression in enpp2 overexpressing embryos showed de-clustered pattern of foxj1 expression pattern (Fig. 2J). In control embryos foxj1 expression was clustered and compacted in an ovoid-like shape; in enpp2 -overexpressing embryos foxj1 was significantly de-clustered to a linear domain (arrow, Fig. 2J).

Enpp2 overexpression induces defects in the midline axis and in the KV organogenesis

The KV formation is dependent on the correct migration and clustering of distinct precursor cells, the dorsal fore runner cells (DFC), during mid-gastrulation (75% epiboly stage)44,45. The DFC express sox17 during early embryogenesis45. In enpp2 overexpressing embryos, the expression of sox17 is readily observed, indicating that the specification of the DFC is not altered. However, the cluster formation of the DFC is significantly impaired in enpp2 overexpressing embryos (51.2 ± 6.2%) compared to control (13.1 ± 1.1%). Enpp2 overexpressing displayed an altered pattern of sox17 expression (Fig. 3A,B), similar to the one observed with foxj1 expression (Fig. 2F). The disorganisation of the DFC could lead to reduced or absent of KV formation. By measuring and observing the KV at stage 6–8ss, we detected a smaller or absent KV in enpp2 overexpressing embryos (Fig. 2E). The de-clustered expression of sox17 and foxj1 expressions suggest the observed phenotypes could be caused by the impaired cell migration during gastrulation.

Figure 3
figure 3

Enpp2 overexpression induces alteration in cellular migration, the expression pattern of KV precursors, dorsal organizer and midline genes (gsc and shha) during early embryogenesis.

(A) Representative images of WISH with sox17 riboprobe to mark the integrity of the DFC cluster at 75% epiboly. (B) Quantification of KV defect and/alteration of DFC cluster, as highlighted by sox17 WISH. (C) The gsc expression is observed in the dorsal organizer of enpp2 overexpressing embryos. However, the expression is also detected in the axial mesoderm precursors suggesting that convergence and cell migration during midline formation is impaired. (D) Representative images of WISH of shha at 90% epiboly showing a truncated axis and laterally expanded expression of shha at the blastomderm margin. (E) Quantification of the broadened shha expression domain at the blastoderm margin in control and enpp2 overexpressing embryos at 90% epiboly. (F,G) Representative images of time lapse series showing cell migration during midline formation (10 hpf- early somitogenesis) in wild type and enpp2-injected embryos. Image stacks were taken every 15 minutes. Time course for each image is indicated. The enpp2 overexpressed embryo failed to undergo a proper CE (G) which results in bend and broadened midline compared to wild type (F). (B,E) The sample size (n) is stated as numerical value above each bar. Data are mean ± SEM. Statistical analysis was established by t-test; *P < 0.05; **P < 0.01. n indicated above bars.

Since we found the convergence and extension movements (CE) during midline formation were impaired, we next assessed if the dorsal organizer is formed correctly following enpp2 overexpression by performing WISH using goosecoid (gsc), a dorsal organizer marker46, during the mid epiboly stage (Fig. 3C). We observed gsc expression in the dorsal organizer but gsc expression was also detected in the axial mesoderm, suggesting that precursors do not properly converge into the dorsal organizer during midline formation (Fig. 3C). Consistent with this result, the expression of the axial mesoderm marker shha showed expanded expression at the blastoderm margin and a shortened anterior to posterior expression pattern in enpp2 injected embryos (61.96 ± 6.8% in the enpp2 overexpressed embryo compared to 1.58 ± 1.6% in control, Fig. 3D,E) confirming an impaired convergence extension during gastrulation. We next performed time-lapse analysis of cell movement from the end of gastrulation until early somitogenesis to observe the midline formation. Cell migration defects were detected during epiboly in enpp2 overexpressing embryos (Fig. 3F,G, Suppl. Movie 12). In wild type embryos, cells intercalated and initiated the elongation during mid-line formation, the cells in the enpp2 overexpressing embryos failed to converge and extend, resulting in a broader midline area and a bent midline formation (Fig. 3F,G, Suppl. Movie 12). The detected cell migration phenotypes were in agreement with the altered gene expression patterns of gsc, ntl, spaw and shha in enpp2 overexpressing embryos.

enpp2 modulates cell migration via lpar1–3 Rho/ROCK

In order to dissect the signalling mechanisms underlying LPA’s effect in development, we first examined the expression profile of its receptors (lpa1–3) in the early embryo by WISH. We detected expression of lpa1–3 in the early embryo (Suppl. Fig. 2 and data not shown). We found that lpa1–3 were ubiquitously expressed at the animal portion and at the margin between cells and yolk portion (yolk syncytial layer) during blastula until onset of gastrulation. During the late gastrulation and early segmentation stages, lpa1 was expressed in the mesoderm tissue adjacent to the body axis, while lpa2–3 were expressed along the midline during zebrafish (Suppl. Fig. 2). Taken together, the expression of lpa13 suggests that these receptors may play an important role during gastrulation and the midline formation.

We next addressed the role played by LPA receptors in the phenotypes observed after enpp2-overexpression. To this end we performed rescue experiments using the well-established LPA receptor (Lpa1–3) antagonist Ki1642547. Ki16425 is a specific lpa receptor antagonist that efficiently antagonize lpa1–3 but not lpa4, lpa6a, and lpa6b and it has been used to antagonise LPA signalling in zebrafish previously34. We incubated enpp2-overexpressing embryos immediately after RNA injection with increasing concentrations of Ki16425 and assessed the impact on phenotypes by using shha and foxj1 expression as a read-out at the end of gastrulation (Fig. 4). enpp2 injected embryos treated with Ki16425 showed a significant reduction of ennp2-induced CE and ciliogenesis phenotypes in a dose-dependent manner (Fig. 4A–C). Following enpp2 overexpression, the CE phenotype characterised by a short axis and a broadened shha expression was detected in 66.5 ± 7.8% of the injected embryos (Fig. 4E). However, following and Ki16425 treatment, the number of embryos displaying this phenotype decreased to 41.3 ± 4.9% (Fig. 4E). Similarly, declustered foxj1 expression was observed in 53 ± 3.8% of the vehicle control embryos but decreased to 28.3 ± 6.2% in Ki16425-treated embryos (Fig. 4E). The rescue effect was maintained until later stages of development as we observed significantly reduced morphological phenotypes and almost normal expression levels of ntl, shha and its downstream gene gli2 following Ki16425 treatments at 24 hours post fertilisation (hpf, Fig. 2, Suppl. Fig. 3).

Figure 4
figure 4

Pharmacological blocking of lpa1–3 or Rho/ROCK signalling rescue enpp2 overexpression induced midline phenotype.

Representative WISH images of shha (A,B) and foxj1 expression (C-D) at 90% epiboly of enpp2 overexpressing embryos treated with the lpa1–3 antagonist Ki16425(A,C) and Rho/ROCK signalling inhibitor Y27632 (B,D). Enpp2 overexpressing embryos display expanded expression of shha (A,B) and foxj1 (C,D) illustrated by pointed in arrows. Quantifications of the of shha expression domain in enpp2 overexpressing embryo following vehicle or Ki16425 (E) or Y27632 treatment (F). The phenotype penetrance following enpp2 injection and Ki16425 treatment and Y27632 were measured at 90% epiboly (E,F). (E,F) The sample size (n) is stated as numerical value above each bar from at least two independent experiments. Data are mean ± SEM. Statistical analysis was established by t-test; *P < 0.05; **P < 0.01. n indicated above bars.

Rho/ROCK is one of the common downstream pathways of enpp2/LPA signalling and it plays an essential role in regulating cell shape and migration through actin stabilization of the cytoskeleton16,48. Furthermore, the Rho/ROCK pathway is an important regulator of CE cell movement during midline formation49,50. We inhibited the Rho/ROCK pathway by using a selective inhibitor of the Rho-associated protein kinase p160 ROCK, Y2763251. Inhibiting Rho/ROCK after enpp2 overexpression significantly reversed the observed phenotypes. The penetrance and severity of CE and clustering of DFC were significantly reduced after Y27632 treatment (Fig. 4B,D,F). Y27632 treatment reduced the expanded shh expression phenotype from 62.4 ± 8.5% in the vehicle control to 33.9 ± 4.9% in treated embryos (Fig. 4F). The foxj1 expression pattern was 53 ± 3.8% in the vehicle control and decreased to 29.2 ± 4.3% in Y27632-treated embryos (Fig. 4F). The rescue effect of Y27632 in enpp2 injected embryos was maintained until later stages of development as we observed significantly reduced morphological phenotypes and almost normal expression patterns of ntl and spaw at 24 hpf (Suppl. Fig. 3). Taken together this data indicates that enpp2 induces the midline axis phenotype through activation of the Rho/ROCK pathway in a lpa1–3 mediated manner.

Discussion

The role of LPA signalling in early embryonic development is limited19,22,23,24,31. Here we found a role of enpp2/LPA signalling in cell migration during gastrulation, which interfered with the midline axis formation and the establishment of the L-R asymmetry.

The establishment of the axial/midline mesoderm and the L-R asymmetry are crucial to the formation of the body plan in the zebrafish embryo. The body axis arises from the dorsal organizer shield6,52,53. Midline formation involves a complex morphogenetic movement, called the convergence and extension movement. During this process, the axial mesoderm cells migrate together and intercalate (converge) towards the dorsal organizer and then migrate towards the anterior and posterior poles (elongation) to give rise to the future midline4,5,54. Enpp2 is expressed at the blastula stage in the zebrafish and it is highly expressed in the margin and midline during gastrulation and somitogenesis respectively. During epiboly, enpp2 is expressed in the YSL and it is later expressed in the dorsal organizer suggesting a role in tissue patterning. lpa13 are also expressed in the embryo during gastrulation. During early gastrulation lpa1 is expressed at the dorsal organizer shield, where the CE occurs to form the body axis54. At the end of gastrulation lpa1 expression is found in the mesoderm adjacent to the midline. In contrast, lpa2b3 are expressed at the midline and at the end of gastrulation. These expression patterns suggest that enpp2/LPA signalling is involved in epiboly and gastrulation. To assess the role of enpp2 in the zebrafish embryogenesis, we overexpressed enpp2 mRNA and assessed the embryonic phenotype at different stages of the zebrafish development. We observed an alteration of cell migration towards the midline during gastrulation. Cells failed to complete the CE e.g. to migrate towards the dorsal organizers and to elongate to form the midline, resulting in perturbed expression of midline and lateral plate mesoderm genes (shha ntl and spaw), bent body axis and shortened body length. The impaired cell migration during CE is likely responsible for the expanded expression of midline genes, bend/kinked midline in the enpp2 -overexpressing embryo and the shortened body plan at later stages of development. The Rho/ROCK pathway is a common and a likely downstream pathway of the enpp2/LPA signalling axis, and particularly through multiple LPA G-protein coupled receptors, including the lpa1–616,48. Rho/ROCK is a key regulator of cell shape and migration through its action on actin stabilization of cytoskeleton49. Rho/ROCK pathway is an important regulator of CE process during midline formation and regulates left-right asymmetry49,50. Indeed ROCK has been previously shown to play an obligatory role in the morphogenetic movement during embryonic organogenesis. Inhibiting ROCK with Y27632 in the early chick embryo blocked cell migration and fusion of the bilateral heart primordia, axis formation, movement of Hensen’s node (similar to KV), and establishment of L-R asymmetry50. We demonstrated that the enpp2 overexpression midline defect phenotype could be rescued by blocking ROCK, suggesting that the phenotype is mediated by Rho/ROCK. In agreement with our findings, the Rho/ROCK pathway regulates CE cell movement during gastrulation49,55. rock2b in zebrafish is specifically expressed in DFC and KV and later in midline56. Knock-down of rock2b function using morpholinos leads to altered anterior posterior patterning, random distribution of ciliated cells and loss of L-R symmetry in the KV56. Similarly, loss of Rock2 function in mice alters organ laterality and asymmetry gene expression including the heart, hypochord and notochord50. Interestingly, over activation of sphingosine-1-phosphate (S1P) signalling via its cognate G protein-coupled receptor S1pr2 results in a very similar endoderm convergence phenotype as enpp257. This phenotype is mediated via a RhoGEF-dependent pathway57 and implies shared downstream mechanisms. Rho/ROCK is also downstream of the planar cell polarity (PCP) pathway that is a well-known regulator of CE during gastrulation. The vertebrate PCP pathway is a non-canonical Wnt signalling cascade that can modulate the actin cytoskeleton via small GTPases including Rho58 to control polarized cell behaviours such as CE movements and orientation of stereocilia in the inner ear49,59,60,61,62. An intriguing possibility is that LPA signalling interacts with the PCP pathway via the extracellular matrix. Enpp2 has several positively charged surface residues and it may interact with heparin sulphate proteoglycans like glypicans that are well known regulators of the PCP pathway63,64.

Furthermore, overexpression of enpp2 altered the migration and clustering of dorsal fore runner cells. This leads to perturbed formation of the KV and ciliogenesis, and subsequently to a failure to establish a correct L-R asymmetry in the embryo. Overexpression of enpp2 could result in effects independent of LPA signalling. However, the efficient pharmacological rescue with the specific LPA receptor (Lpa1–3) antagonist Ki16425 suggests that the phenotypes observed are mediated by LPA receptors. Furthermore, our findings are consistent with the previously reported roles of enpp2/LPA signalling during L-R patterning in mice and zebrafish. Loss of enpp2 in mice leads to phenotypes such as: axial mesoderm/notochord defects, and abnormal organ positioning29,65, and are associated with abnormal L-R patterning66 Furthermore, knock-down of enpp2 or lpa3 in zebrafish using morpholinos results in altered KV formation, ciliogenesis and L-R asymmetry31. In addition, over activation of enpp2/LPA signalling induces cardiac bifida in zebrafish32. Cardiac bifida is caused by the failure of bilateral mesoderm to migrate and coalesce into a single central heart tube67. Similar to enpp2/LPA signaling, over activation of sphingosine-1-phosphate (S1P) signalling via its cognate G protein-coupled receptor S1pr2 impairs endoderm convergence and results in cardia bifida57. Taken together, available data suggests that LPA signaling needs to be tightly regulated during multiple stages of embryonic development.

Conclusion

Midline formation and establishment of L-R symmetry is crucial in the embryonic development. By overexpressing enpp2 in the zebrafish embryo, our study for the first time shows the unique role of enpp2/LPA in regulating mesendoderm cell migration in the establishment of axial midline and L-R asymmetry in the embryogenesis through Rho/ROCK pathway via lpa1–3.

Material and Methods

All experimental work performed in this study was approved by the Monash Animal Ethics Committee (MARP/2013/096), in accordance with the requirements of the National Health & Medical Research Council of Australia.

Zebrafish

Zebrafish embryos were collected immediately after fertilization, maintained at 28.5 °C, and staged by developmental time (hours post-fertilization, hpf) using morphological criteria68. The Tübingen was used as a wild type strain. To perform in situ, the embryos at the designated stages were fixed in 4% PFA (2 hrs RT (<24 hpf) or overnight 4 °C)), washed with PBS-0.1% Tween (PBS-T, 5 min, twice), dechorionated, dehydrated with MeOH (5 min RT, twice) and stored in MeOH at −20 °C until processed for in situ.

Microinjection and phenotypic analysis

The embryos were injected with 1 nl of the corresponding RNA at different concentrations. Full-length enpp2/atx constructs (25, 50, 75, 100, 150 and 200 pg) or Arl13b GFP mRNA (50–100 pg) or H2B (75 pg) were injected into the 1–4 cell stage embryos. Phenotype penetrances and mortality rates were quantified within 9–24 hrs post injection. For assessment of Kupffer’s vesicle (KV), embryos were examined at the 5–8 somite stages.

Quantitative Polymerase Chain Reaction (qPCR)

qPCR was carried out using TaqMan Universal master mix (Applied Biosystems, Foster City, CA) and the 7900HT Fast Real-Time PCR system (Applied Biosystems) and SYBR Green (Sigma). Primer set targeting enpp2: forward (TTCTCCTTCATCCTTCCACAC) and reverse (GTAACTCCACATCTCGCAGG) (NM_200603) were used. A total of 10 μl reaction was prepared and repeated in triplicate per test. The relative quantitation was achieved by applying the comparative CT method (ΔΔCT) whereby the mRNA levels were normalized against the level of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) mRNA with enpp2 expression at the shield stage used as the reference.

Whole Mount In Situ Hybridization (WISH)

WISH was performed following the established method with minor modifications69,70. The primer sequence used for PCR amplification and cloning the following genes, lpa13, enpp2/atx, shha, foxj1 (Suppl. Table 1). The resulting DNA fragments were subcloned into pGEMT-easy vector or pCS2+ vector; of which were digested to be transcribed into sense (enpp2) and/or antisense mRNA sequence/probe.

Embryo Immersion/Drug treatment

Following enpp2/atx mRNA injection, embryos were immersed in Ringer’s solution containing 0.5% DMSO and supplemented with the Rho/ROCK inhibitor, Y27632 (5 and 10 μM, Sigma-Aldrich) or lpar1–3 antagonist, Ki16425 at different concentrations (1, 5, 7.5 and 10 μM Sigma-Aldrich). Y27632 was used at 10 μM and Ki16425 was used at 7.5 μM for the rescue experiments. WT control embryos immersed with the drugs from early gastrulation and embryos did not develop detectable phenotype at 24 hpf. The embryos were then collected and fixed with 4% PFA and processed for WISH.

Live Imaging

To visualize the KV, wild type or over-expressed zebrafish embryos were dechorionated at 5–8ss stage and mounted in 0.8–1% of Low Melting-point Agar (LMP, Sigma). The KV was then visualized using an Olympus Dissecting microscope MVX10 using the CellSens Standard software (v8.1, Olympus). In some experiments, cilia in embryos KV were visualized under a Leica SP5 confocal microscope and the Leica application suit advanced software (v3, Leica) using a 20 × 1.0 NA water immersion lens. For imaging of cilia, Z-projections of the entire KV were captured with slow acquisition with some frame averaging using arl13b-GFP, a live axoneme marker of cilia, to investigate cilia in the living embryo as studied in ref. 41. To image cell migration during midline axis formation, the embryos were injected with H2B together with or without enpp2/atx and visualized starting at the tail bud stage (10 hpf) and recorded by time lapse for 8–10 hours. The embryo was visualized every 10–15 minutes.

Statistical analysis

All sets of experiments were performed at least three times in triplicates, unless specified (n refers to the number of independent experiments performed on different cell cultures). Data-sets were expressed as mean ± standard error of the mean (SEM). Significance of the differences was evaluated using the unpaired t-test (two-tailed) or the one and two-way ANOVA followed by Dunnett’s multiple comparisons test in the in vivo experiment. Statistical significance was established at *p < 0.05, **p < 0.01 and ***p < 0.001, ****=<0.0001. α = 0.5, using Prism (v.6, GraphPad).

Additional Information

How to cite this article: Frisca, F. et al. Role of ectonucleotide pyrophosphatase/phosphodiesterase 2 in the midline axis formation of zebrafish. Sci. Rep. 6, 37678; doi: 10.1038/srep37678 (2016).

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