Identification and Expression Analysis of Zebrafish Glypicans during Embryonic Development

Abstract

Heparan sulfate Proteoglycans (HSPG) are ubiquitous molecules with indispensable functions in various biological processes. Glypicans are a family of HSPG’s, characterized by a Gpi-anchor which directs them to the cell surface and/or extracellular matrix where they regulate growth factor signaling during development and disease. We report the identification and expression pattern of glypican genes from zebrafish. The zebrafish genome contains 10 glypican homologs, as opposed to six in mammals, which are highly conserved and are phylogenetically related to the mammalian genes. Some of the fish glypicans like Gpc1a, Gpc3, Gpc4, Gpc6a and Gpc6b show conserved synteny with their mammalian cognate genes. Many glypicans are expressed during the gastrulation stage, but their expression becomes more tissue specific and defined during somitogenesis stages, particularly in the developing central nervous system. Existence of multiple glypican orthologs in fish with diverse expression pattern suggests highly specialized and/or redundant function of these genes during embryonic development.

About the Authors:

Mansi Gupta

Affiliation: DFG- Center for Regenerative Therapies Dresden/ Cluster of Excellence (CRTD), Biotechnology Center, Technische Universität Dresden, Dresden, Saxony, Germany

Michael Brand

* E-mail: [email protected]

Affiliation: DFG- Center for Regenerative Therapies Dresden/ Cluster of Excellence (CRTD), Biotechnology Center, Technische Universität Dresden, Dresden, Saxony, Germany

Introduction

The members of the glypican family are extracellular matrix components which play essential roles in various biological processes. They are Heparan Sulfate Proteoglycans, composed of a cysteine-rich protein core to which heparan sulfate (HS) sugar chains are covalently attached at the C-terminal end. They usually mediate their function anchored on the exoplasmic cell membrane via a GPI-linkage, but can also be cleaved off the membrane and in some cases enter the circulation [1].

The HS sugar chains, being highly sulfated, attract a number of growth factors, due to their negative surface charge. Consequently, glypicans are known to modulate the activity of various growth factors like Wnt, Hedgehog (HH), Fibroblast growth factors (Fgfs) and BMP [2-5]. However not all functions of glypicans are mediated by their HS chains [6]. There are 6 glypican genes in humans (GPC1-6) and two in Drosophila (dally and dally-like), which code for the protein core [7]. Among glypicans, dally was the first mutant to be isolated in Drosophila displaying developmental defects in the eye, brain and wing [1]. This was attributed to a reduction in DPP signaling. Since then, various genetic and biochemical studies have associated glypicans to different developmental processes, growth and disease progression. They can function as low affinity co-receptors for growth factors or assist in their transport across cells [8]. In rodents, Gpc1 regulates brain size via the modulation of Fgf signaling [4]. It is over-expressed in tumors from different organs [9] and was recently associated with the pathogenesis of a liver disorder, biliary artresia [10]. GPC3 mutations in humans and mice result in Simpson-Golabi Behmel syndrome, characterized by pre- and postnatal skeletal anomalies and craniofacial malformations [11]. There is also a very high occurrence of GPC3 over-expression in hepatocellular carcinoma where it is considered as a serum marker and potential therapeutic target [12]. Gpc5 is a very specific enhancer of HH signaling and stabilizes the interaction between HH and its receptor Patched1 [3]. Gpc5 knockdown was seen in patients with neural tube defects [13] and this gene has also been implicated in nephrotic syndrome [14].

Among all glypicans, Gpc4 influences extremely diverse development processes. In Xenopus and zebrafish, Gpc4 regulates the convergent extension movements during gastrulation [2,15]. Its interaction with Fgf2 is also important for proper forebrain patterning in Xenopus [16]. Astrocytes secrete Gpc4 and Gpc6 which guide the formation of excitatory synapses [17]. Gpc4 also regulates insulin signaling via its interaction with the insulin receptor, and its levels in circulation correlate with an increased BMI [1]. Finally mutations in Gpc6 result in omodysplasia, characterized by shortened limbs and facial dysmorphism [18].

In contrast to humans and mice, the study of glypican function in zebrafish has so far been limited to gpc3 and gpc4. A gpc4 mutant in fish (knypek) displays defective convergent extension movements during gastrulation resulting in a reduced body length [2]. Gpc4 is directly involved in mediating non-canonical Wnt signaling in the embryos. Rescued knypek mutants also display defects in craniofacial cartilage development in larval and adult stages [19]. Gpc3 inhibits canonical Wnt β-catenin signaling after getting cleaved from the surface by Notum hydrolase and this regulation is also necessary for proper gastrulation [20].

We have previously found by single molecule analysis that extracellular matrix composition influences the mobility of Fgf8 forming a morphogen gradient in the extracellular matrix [21,22] Due to a lack of detailed information on fish glypicans as extracellular matrix components, we began to systematically characterize these genes. Firstly we identified and isolated 10 glypican genes from the fish. We studied their phylogeny with respect to humans and analyzed their expression pattern at various stages of embryonic development. Our findings indicate that glypicans are overall conserved between mammals and zebrafish, and may serve both generalized and highly tissue-specific functions in developing tissues.

Materials and Methods

Ethics Statement

All animal experiments were carried out in strict accordance with European Union and German laws (Tierschutzgesetz). All experimental procedures were approved by the animal ethics committee of the TU Dresden and the Landesdirektion Sachsen (approval number: AZ 24D-9168.11-1/2008-4). This institutional review board specifically approved this study.

Zebrafish husbandry

Zebrafish were raised and maintained as described previously [23]. The wild-type line used was TL. Zebrafish embryos were obtained by natural spawning of adult fish and staged according to hours post fertilization (hpf) or standard criteria [24]

Bioinformatics analysis

Ensembl Zv9, GenBank and DFCI EST databases were used to identify zebrafish glypican sequences. Human glypican sequences were blasted against all three databases and the obtained fish sequences were confirmed for the presence of the glypican domain. Mega5.1 software was used for phylogenetic analysis. ClustalW and the Jalview software were used for generating and viewing the multiple sequence alignment, respectively. Protein sequence similarity and identities were obtained using NCBI BLASTp. Cinteny server (http://cinteny.cchmc.org/) and genomicus server were used for syntenic analysis.

Molecular Cloning

Zebrafish mRNA was isolated from 24 and 48 hpf old embryos using the Trizol/Phenol-Chloroform method. cDNA was prepared from RNA using SuperScriptIII First-Strand Synthesis system (Invitrogen). The open reading frames of glypican genes were cloned into Topo vector using the primers given in Table 1. gpc2 full length sequence was obtained by performing RACE using SMARTer RACE cDNA amplification kit (Clontech) with primer: CAGCCCTGAAACACCTTAGCAGAGA for 5’ RACE and primer: AGACGCGCGGCAGGTACCTGCCAGCAG for 3’ RACE.

Gene Name Forward Primer Reverse Primer

Primers for cloning

gpc1a ATGGATCTGACAGCGGTCGC TGATCTAGATTATCGTCTGAGCAATAGACTC

gpc1b ATGGGTTTTGTCTCGCTGGT GCATCTAGATTATCGCCTAAGCAAGACTGT

gpc2 ATGAAGATGATGAAGGTGGTGATGAAGAT TTACAGACACAGACAGAGAGTAAAGC

gpc3 ATGATGCCTGGACTGAAGTTG AGCTCGAGATCACTGAAGACCCAGTGTTATG

gpc4 CCTGGATCCATGAAGATGATCGTTGTGTT GTGCTCGAGTTATCTTGTTTGGAGAGTGA

gpc5a ACGGGATCCATGTCTCTTCACCTAAAATC GTCTCGAGTTCATTGTGTCTCCTTACTGG

gpc5b AATGCTCCGCGGGACAGCA TACTCGAGATCAAGGCCACAGGAGGAGCT

gpc5c CGAATTCATGTCACGCGTGAATGTCAGCT TCTCGAGTTCGACCCTTAGAGCTAAGGTATG

gpc6a TAGGATCCATGGTGAAGACACCTGTCGT CACTCGAGTTATCTCCAGCAGAGCAGCA

gpc6b AAGGATCCATGTGGGCTGTGTGCGCGCT GTCTCGAGCTATCTGCAGGTCAGTGCCA

Primers for RT-PCR

gpc2 CGA AGC TGC GCG GAG TCC CG ACCAGTTGGTCCAGGTCCGTCC

gpc3 AGTGCGCTCATCCTTTCAGT CATACGCAGAGCTCCCTTTC

gpc5a GGATACTGCCTGAACGTGGT TTTTCAGCCTGCGCTTTAAT

gpc5b TCTACCTTCGAGGCGATGAT CTCGCATCACGTTGAGACAT

gpc6a CACCAAAGGCTTCAGTCTCC CCATAAGTCCCTGGACGAAG

gpc6b ACGTCAACCTGGAGGAGATG ATGGGCTCCATCACAGACTC

Table 1. Primer list.

CSV

Download CSV

RT-PCR

cDNA was prepared from different developmental stages as mentioned above and used for RT-PCR. Full length primers (Table 1) were used to amplify gpc1a, gpc1b, gpc4 and gpc5c. The primers used for gpc2, gpc3, gpc5a, gpc5b, gpc6a and gpc6b are also mentioned in Table 1. The PCR reaction was carried out using DreamTaq DNA polymerase (Fermentas) for 30-35 cycles.

In Situ Hybridization

Embryos at the desired stage were fixed in 4% PFA and in situ hybridization was carried out according to the protocol described previously [25]. Full-length antisense probes were synthesized using T7/SP6 polymerase, Digoxigenin label (Roche) and linear Topo vectors. The color was developed using BM purple and used for imaging. Flat mounts were prepared after removing yolk sac and mounted in glycerol.

Results

Identification of Fish Glypicans

There are 6 glypican genes in the human genome (GPC1-6). These sequences were blasted against Ensembl Zv9, GenBank and DFCI zebrafish EST databases to obtain the corresponding zebrafish glypicans. A previous study also reported 6 fish glypicans [10] but we identified 10 genes and named them according to their similarity to the human orthologs. The accession numbers of all genes from the different databases are shown in Table 2. New GenBank accession numbers were obtained for most sequences. Also indicated are the corresponding human orthologs and their peptide length.

Gene Name Chromosome No. Ensembl Gene ID GenBank Acc No. (old) DFCI Identifier GenBank Acc. No. (new) Peptide length (a.a.) Human homolog/Peptide length Comments

gpc1a 22 ENSDARG00000019341^* BC109411 TC366647 KC836776 554 HsGPC1/558 Gpc1^**

gpc1b 2 ENSDARG00000090456 BC053161^* TC368439 KC836777 541 HsGPC1/558

gpc2 14 ENSDARG00000037150 BC133103^* TC411825 KC791423 576 HsGPC2/579 Identified by 5’ and 3’ RACE

gpc3 14 ENSDARG00000032199 XM_682922^* TC371608 - 590 HsGPC3/603 [31]

gpc4 14 ENSDARG00000015472^* NM_131860 TC419247 - 557 HsGPC4/556 [2]

gpc5a Zv9_NA ENSDARG00000088858 XM_001920514 TC414869^* KC791424 582 HsGPC5/572

gpc5b 22 ENSDARG00000024588^* -- TC376562 KC791425 610 HsGPC5/572 Gpc5^**

gpc5c 2 ENSDARG00000074082^* -- -- KC999393 523 HsGPC5/572

gpc6a

1. 12

2.

1. ENSDARG00000091739

2.

1. BC151902^*

2.

3.

TC370142 KC999394 562 HsGPC6/555

gpc6b 9 ENSDARG00000036468 XM_003199260^* NP13316404 KC999395 543 HsGPC6/555 Gpc6^**

Table 2. Accession numbers of zebrafish glypicans.

* cloned sequences; HsGPC: Human glypican

** Previously reported names [10]

CSV

Download CSV

Glypican gene family members can be grouped into two sub-families, as reported previously [7]: GPC1/2/4/6 and GPC3/5 family (Figure 1). Phylogenetic comparison of fish, human and Drosophila sequences revealed that all fish genes cluster with their human orthologs. The two genes in Drosophila, dally and dally-like protein each belong to a different family. Multiple orthologs of human genes are commonly found in fish genome due to a duplication event which occurred before the radiation of teleosts [26]. Corresponding to HsGPC2 we identified one new ortholog in the fish genome (Drgpc2). Two new orthologs were identified corresponding to HsGPC1 (Drgpc1a and Drgpc1b) and HsGPC6 (Drgpc6a and Drgpc6b) each, and 3 corresponding to HsGPC5 (Drgpc5a, Drgpc5b and Drgpc5c). All glypican proteins exhibit more than 56% sequence similarity with their corresponding human orthologs (Table 3).

[Figure omitted. See PDF.]

Figure 1. Phylogenetic tree of glypican genes from Danio rerio, Homo sapiens and Drosophila melanogastor.

Phylogenetic comparison was carried out for glypicans genes sequences from zebrafish (Dr in black), humans and Drosophila (Hs and Dm in grey). The evolutionary tree was constructed using Neighbor-Joining method in the MEGA5.1 software using the Poisson model and partially deleted dataset. 500 bootstrap replications were used as a test of phylogeny and the values are indicated next to the branch. Branch length corresponds to evolutionary distances which denote the number of amino acid substitutions per site. Scale bar: 0.1

https://doi.org/10.1371/journal.pone.0080824.g001

Identity % Similarty %

DrGpc1a vs HsGPC1 54 70

DrGpc1b vs HsGPC1 52 68

DrGpc2 vs HsGPC2 46 63

DrGpc3 vs HsGPC3 43 58

DrGpc4 vs HsGPC4 58 72

DrGpc5a vs HsGPC5 55 73

DrGpc5b vs HsGPC5 39 56

DrGpc5c vs HsGPC5 39 59

DrGpc6a vs HsGPC6 68 81

DrGpc6b vs HsGPC6 60 77

Table 3. Protein sequence similarity.

CSV

Download CSV

In order to identify true orthologs of HsGPC1, HsGPC5 and HsGPC6, we investigated the chromosomal syntenic relationship between fish and human glypican genes by comparing common markers flanking these genes. The gene cluster containing GPC3 and GPC4 is well conserved in Eumetazoa and is also seen in zebrafish (Figure 2B) [7]. Of the two gpc1 in fish, gpc1a shares one adjacent marker with the HsGPC1 (Figure 2A), although both Drgpc1a and Drgpc1b share an almost similar amino acid sequence identity with HsGPC1 (Table 3). This indicates that Drgpc1a is likely to be a functional ortholog of HsGPC1. Fish has two orthologs for HsGPC6, namely Drgpc6a and Drgpc6b. Both genes share some synteny with the human form, such that certain markers are common between HsGPC6 and Drgpc6a and others are common between HsGPC6 and Drgpc6b (Figure 2C and 2D). This suggests that both zebrafish genes are derived from the same ancestral locus but have undergone rearrangements resulting in the final gene arrangements. Currently, Drgpc6a is annotated on two chromosomes, Chr1 and Chr12 (Ensembl Zv9), but only the Chr12 version is syntenic with HsGPC6.

[Figure omitted. See PDF.]

Figure 2. Syntenic conservation between zebrafish and humans glypicans.

The figure depicts arrangement of genes surrounding glypicans on corresponding chromosomes in humans and zebrafish. Chromosome numbers from humans (Hs) and zebrafish (Dr) are indicated on the left. Glypican genes are represented in red and conserved markers around glypicans are in green or yellow. Non-conserved genes are in grey. The diagram is not to scale.

https://doi.org/10.1371/journal.pone.0080824.g002

In the human genome, GPC5 and GPC6 are present adjacently, but this contiguity is lost in the fish genome. There are 3 copies of gpc5 in fish and Drgpc5a has the maximum amino acid sequence identity to HsGPC5 (Table 3) but none are syntenic with humans. On fish Chr22, Drgpc5b and Drgpc1a are present contiguously, an arrangement not seen in mammals. Absence of synteny implies that along with whole genome duplication, individual gene duplications for fish gpc5 might have resulted in the observed multiple paralogs.

Glypican conserved domain structure

The glypican family of proteins is characterized by a large and highly conserved N-terminal glypican domain and a C-terminal region containing the residues for heparan sulfate and Gpi-anchor attachment. Multiple sequence alignment confirmed the presence of 14 conserved cysteine residues in all zebrafish glypicans (Figure 3). All these residues are involved in disulfide linkages and are necessary to maintain the folded structure [27]. The HS attachment site is always present within 50 residues of the C-terminal placing them very close to the plasma membrane [28]. It consists of repeating Ser-Gly (SG) cluster (n≥2) flanked by acidic residues on both sides [29]. There are 4 SG repeats in Gpc1a, Gpc1b, Gpc6a and Gpc6b; 3 repeats in Gpc4; 2 repeats in Gpc2, Gpc3, Gpc5a and 1 repeat in Gpc5b. Gpc5c and also Gpc5b contain a separate SG dipeptide followed by acidic amino acids, a potential chondroitin sulfate priming motif. Since HsGpc5 is known to carry both heparan sulfate and chondroitin sulfate chains [30], this feature may be conserved even in zebrafish Gpc5.

[Figure omitted. See PDF.]

Figure 3. Multiple sequence alignment of glypican protein sequences from zebrafish.

All protein sequences were aligned to compare conserved features between glypicans. Residues highlighted in blue are well conserved. Black highlights indicate the 14 conserved cysteine residues; yellow highlights are the multiple Ser-Gly repeats flanked by acidic residues for sugar attachment; in red is the potential chondroitin sulfate attachment site.

https://doi.org/10.1371/journal.pone.0080824.g003

Spatiotemporal expression of glypicans during early development

The composition of extracellular matrix differs considerably between specific tissues and stages of development. A previous study has focused on the expression of fish glypicans at 5 day post fertilization only in the developing liver [10]. We now analyze the temporal and spatial expression of all glypicans from cleavage to pharyngula stages using Reverse transcriptase PCR (RT-PCR) and in situ hybridization (ISH) (Figures 4-7). As reported before, gpc4 mRNA is present at all stages of embryonic development [2]. A constant expression (both maternal and zygotic) was also seen for gpc2, gpc3, gpc5a, gpc5b, gpc6a and gpc6b. Maternal expression was detected for gpc1b, but zygotic expression started only during segmentation.. Maternal expression was not seen for gpc1a and gpc5c. gpc1a transcripts were first detected at gastrulation stage and gradually increased during segmentation period and gpc5c was detected from late segmentation stages.

[Figure omitted. See PDF.]

Figure 4. Temporal mRNA expression of glypicans genes during early embryonic development.

RT-PCR was performed for all 10 genes using cDNA from specified stages of development. β-actin was used as a control. Ep, epiboly; som, somites; hpf, hours post fertilization.

https://doi.org/10.1371/journal.pone.0080824.g004

[Figure omitted. See PDF.]

Figure 5. Spatial expression of zebrafish glypicans during cleavage to early segmentation stages.

In situ hybridization was performed to detect glypican mRNA during early developmental stages. Lateral views with anterior to the top (A, B, C) and dorsal to the right (B, C). Developmental stages and gene names are indicated. tb, tailbud; hb, hindbrain; en, endoderm; me, mesencephalon; n.d: not detected.

https://doi.org/10.1371/journal.pone.0080824.g005

[Figure omitted. See PDF.]

Figure 6. Spatial expression pattern of zebrafish glypicans during segmentation and pharyngula stages.

In situ hybridization was carried out to detect glypican mRNA at 20 somite stage: 19 hpf and pharyngula stages: 24-26 hpf. Lateral view of whole mount embryos, anterior to the left. Expression pattern of gpc1a (A, A’), gpc1b (B, B’); subfigure in B is an optical section indicating expression in yolk cells (arrowheads) and branchial arch primordium, gpc2 (C, C’), gpc3 (D, D’), gpc4 (E, E’), gpc5a (F, F’), gpc5b (G, G’), gpc5c (H, H’), gpc6a (I, I’), gpc6b (J, J’). di, diencephalon; te, telencephalon; hb, hindbrain; tg, tegmentum; nc, neural crest; pm, paraxial mesoderm; som, somites; ba, branchial arch; y, yolk; mb, midbrain; ff, fin fold; ov, otic vesicle; fp, floor plate; psom, posterior somites; ep, epiphysis; tri, trigeminal placode; cl, cloaca; pllp, posterior lateral line primordium; cg, cranial ganglion.

https://doi.org/10.1371/journal.pone.0080824.g006

[Figure omitted. See PDF.]

Figure 7. Specific expression domains for gpc1a, gpc1b and gpc5b.

Transverse cross sections revealed the expression details of gpc1a, gpc1b and gpc5b. Dotted line indicates the plane of section. pllp: posterior lateral line primordium; pba: posterior branchial arch; my: myotome; fp: floor plate; nc: notochord.

https://doi.org/10.1371/journal.pone.0080824.g007

These findings were confirmed and extended by whole mount in situ hybridization. As seen by RT-PCR, maternal contribution was observed for gpc2, gpc4, gpc5a, gpc5b, gpc6a and gpc6b (Figure 5A). gpc1b was also detected by ISH at cleavage stages but not later during gastrulation. gpc3 transcripts could not be detected either during blastula or gastrula stages by ISH probably due to low levels of transcripts.

Gastrulation marks the process when cell movements result in the formation of the germ layers and embryonic axis. Gpc3 and Gpc4 influence this process in a mutually exclusive manner [20]. Along with gpc4, we also detected gpc2, gpc5a, gpc5b, gpc6a and gpc6b transcripts at 60% epiboly (Figure 5B). Further functional studies will uncover whether all glypican genes function independently or exhibit redundancy with each other during early stages.

We next examined the expression of glypicans during the segmentation period. Beginning from the 2-somite stage, the expression domains of glypicans became more tissue-restricted (Figure 5C). Whereas gpc2, gpc5b and gpc6b were seen ubiquitously on the dorsal side of the embryo, gpc1a, gpc1b and gpc5a showed restricted expression domains in the developing nervous system. gpc1a expression was seen throughout the brain primordium, gpc1b was present near the hindbrain and tailbud region and gpc6a in the mesencephalic region. gpc5a was selectively expressed in the ventral endoderm. gpc3 or gpc5c were not detected.

By the 20 somite stage, all glypicans were detected by ISH. Most genes were expressed in a variety of neural tissues. gpc1a was detected in the primordial telencephalon, diencephalon and hindbrain (Figure 6A). By the 24 hpf stage, its expression became more confined to these regions and was also detected in the posterior lateral line primordium (Figure 6A’, Figure 7A). gpc1b was initially seen in the trigeminal placode and the yolk cells during segmentation (Figure 6B) and at 24 hpf, in somites (Figure 6B’) and the branchial arch (Figure 7B).

gpc3 was expressed very specifically in the midbrain, coinciding with the time of active patterning processes in this tissue (Figure 6D, 6D’). This is in contrast to previous observation where gpc3 was suggested to be restricted to the prospective hindbrain [31]. gpc4 was expressed broadly throughout the embryo but excluded from most dorsal regions, including the telencephalon (Figure 6E. E’). gpc5a was detected in the midbrain, otic vesicle, primordial fin fold and cloaca (Figure 6F) and this pattern persisted at 24 hpf (Figure 6F’). gpc5b was restricted to the floor plate from the 20 somite stage (Figure 6G), but at 24 hpf, its expression was additionally detected throughout the nervous system (Figures 6G’, 7C)

gpc6a was initially very strongly expressed in the epiphysis and the trigeminal placode (Figure 6I) and at 24 hpf, it was also seen in the telencephalon, tegmentum, cranial ganglia and the hindbrain (Figure 6I’). gpc6b was present uniformly in the nervous system (Figure 6J, 6J’).

gpc2 had the most widespread expression domain, seen ubiquitously in the embryo (Figure 6C, 6C’). In contrast, gpc5c was expressed outside the nervous tissue, in posterior somites (Figure 6H, 6H’).

Discussion

A previous study of glypican expression pattern in mice revealed that these genes are expressed predominantly in the embryonic brain [32]. In support of this, we also observed widespread expression of zebrafish glypicans in the developing nervous system. Apart from gpc5c, which is found in the posterior somites, all zebrafish genes are present either in very specific domains or more generally in the brain. Hence glypicans might have a conserved function during the development and patterning of the nervous system. Several studies in mammals have revealed the importance of glypicans in brain development [4,13] . In Xenopus also, Gpc4 is required for proper patterning of the forebrain [16]. Interestingly, the zebrafish gpc4 mutant, knypek, displays craniofacial skeletal defects during the larval to adult stages due to improper chondrocyte and cartilage growth ([33].

A direct comparison of expression domains between mouse versus zebrafish glypicans reveals interesting analogous features. During rat and mouse embryogenesis, Gpc1 is predominantly expressed in the developing nervous system and skeletal system. It is specifically present in the proliferating neural progenitors of forebrain, midbrain and hindbrain at E14/ E18 pharyngula stages [34]. The expression of zebrafish gpc1a seen at 24 hpf corresponds well with that of rat/mouse Gpc1. Moreover we observed partial synteny between human GPC1 and zebrafish gpc1a (Figure 2A). Hence based on comparable expression pattern and partial synteny, we conclude that zebrafish gpc1a is more likely to be the true ortholog of mammalian Gpc1.

In a previous study on gpc3, it was reported that it is ubiquitously present during early zebrafish development [31]. We detected transcripts of gpc3 by RT-PCR at all stages, but these could not be detected by ISH, possibly due to very low levels. Interestingly, gpc3 expression in zebrafish does not correlate with mammalian expression. GPC3 is the causative gene for Simpson-Golabi Behmel Syndrome and is strongly up-regulated in hepatocellular carcinoma. In humans and mice, it is expressed in hepatocytes, several kidney structures, mesenchyme of limb buds, vertebrae and the liver [35]. In contrast, we observed a very specific expression of zebrafish gpc3 only in the midbrain. It is likely that this function of Gpc3 is preserved only in mammals, although later developmental stages in zebrafish remain to be analyzed.

Another parallel can be drawn between the Gpc5 expression patterns. In Xenopus, gpc5 is expressed strongly during neurulation and becomes restricted to the floor plate, somites and optic vesicle in the late neural stages [13]. It has been implicated in neural tube closure defects [13]. Zebrafish gpc5b is also specific to the floor plate and is diffusely present in the nervous system. This indicates that gpc5b might be the cognate zebrafish gene and also suggests for a potential conserved role of gpc5b in neural tube defects. In mice and humans, Gpc5 is prominent in kidney, limb and brain both in developing and adult tissues [30]. But so far, we did not examine later stages of zebrafish development for gpc5 expression.

Taken together, zebrafish glypicans exhibit quite dynamic expression pattern during early embryonic development. The 10 genes identified in this study are selectively turned on in different tissues and hence have the potential to mediate important morphogenetic functions either independently or in conjunction with each other. Multiple signaling pathways are known to be modulated by glypicans in a context dependent manner [2-5], and our data will therefore contribute to elucidating the mechanisms by which glypicans mediate these functions.

Acknowledgments

We thank Matthias Nowak and Anja Machate for support during the early stages of this work, Kalaimathy Singaravelu and Ramu Chenna for help with the bioinformatics analysis, Marika Fischer and Jitka Michling for fish care and Stefan Hans for very useful comments on the manuscript.

Author Contributions

Conceived and designed the experiments: MG MB. Performed the experiments: MG. Analyzed the data: MG MB. Contributed reagents/materials/analysis tools: MG. Wrote the manuscript: MG.

Citation: Gupta M, Brand M (2013) Identification and Expression Analysis of Zebrafish Glypicans during Embryonic Development. PLoS ONE 8(11): e80824. https://doi.org/10.1371/journal.pone.0080824

References

1. Nakato H, Futch TA, Selleck SB (1995) The division abnormally delayed (dally) gene: a putative integral membrane proteoglycan required for cell division patterning during postembryonic development of the nervous system in Drosophila. Development 121: 3687-3702. PubMed: 8582281.

2. Topczewski J, Sepich DS, Myers DC, Walker C, Amores A et al. (2001) The zebrafish glypican knypek controls cell polarity during gastrulation movements of convergent extension. Dev Cell 1: 251-264. doi:10.1016/S1534-5807(01)00005-3. PubMed: 11702784.

3. Kitagawa H, Nakamura H, Kawaguchi Y, Tsuji H, Satone T et al. (1995) Magnetic-evoked compound muscle action potential neuromonitoring in spine surgery. Spine (Phila Pa 1976) 20: 2233-2239. PubMed: 8545718.

4. Jen YH, Musacchio M, Lander AD (2009) Glypican-1 controls brain size through regulation of fibroblast growth factor signaling in early neurogenesis. Neural Dev 4: 33. doi:10.1186/1749-8104-4-33. PubMed: 19732411.

5. Jackson SM, Nakato H, Sugiura M, Jannuzi A, Oakes R et al. (1997) dally, a Drosophila glypican, controls cellular responses to the TGF-beta-related morphogen, Dpp. Development 124: 4113-4120. PubMed: 9374407.

6. Lin X, Perrimon N (2002) Developmental roles of heparan sulfate proteoglycans in Drosophila. Glycoconj J 19: 363-368. doi:10.1023/A:1025329323438. PubMed: 12975617.

7. Filmus J, Capurro M, Rast J (2008) Glypicans. Genome Biol 9: 224. doi:10.1186/gb-2008-9-5-224. PubMed: 18505598.

8. Häcker U, Nybakken K, Perrimon N (2005) Heparan sulphate proteoglycans: the sweet side of development. Nat Rev Mol Cell Biol 6: 530-541. doi:10.1038/nrm1681. PubMed: 16072037.

9. Aikawa T, Whipple CA, Lopez ME, Gunn J, Young A et al. (2008) Glypican-1 modulates the angiogenic and metastatic potential of human and mouse cancer cells. J Clin Invest 118: 89-99. doi:10.1172/JCI32412. PubMed: 18064304.

10. Cui S, Leyva-Vega M, Tsai EA, Eauclaire SF, Glessner JT et al. (2013) Evidence From Humans and Zebrafish That GPC1 Is a Biliary Atresia Susceptibility Gene. Gastroenterology.

11. Pilia G, Hughes-Benzie RM, MacKenzie A, Baybayan P, Chen EY et al. (1996) Mutations in GPC3, a glypican gene, cause the Simpson-Golabi-Behmel overgrowth syndrome. Nat Genet 12: 241-247. doi:10.1038/ng0396-241. PubMed: 8589713.

12. Filmus J, Capurro M (2013) Glypican-3: a marker and a therapeutic target in hepatocellular carcinoma. FEBS J.

13. Bassuk AG, Muthuswamy LB, Boland R, Smith TL, Hulstrand AM et al. (2013) Copy number variation analysis implicates the cell polarity gene glypican 5 as a human spina bifida candidate gene. Hum Mol Genet 22: 1097-1111. doi:10.1093/hmg/dds515. PubMed: 23223018.

14. Kobayashi H, Namba K, Mori T, Nakato Y (1995) Energy distribution of interface states in the band gap of GaAs determined from x-ray photoelectron spectra under biases. Phys Rev B Condens Matter 52: 5781-5788. doi:10.1103/PhysRevB.52.5781. PubMed: 9981767.

15. Aida Y, Kusumoto K, Nakatomi K, Takada H, Pabst MJ et al. (1995) An analogue of lipid A and LPS from Rhodobacter sphaeroides inhibits neutrophil responses to LPS by blocking receptor recognition of LPS and by depleting LPS-binding protein in plasma. J Leukoc Biol 58: 675-682. PubMed: 7499965.

16. Galli A, Roure A, Zeller R, Dono R (2003) Glypican 4 modulates FGF signalling and regulates dorsoventral forebrain patterning in Xenopus embryos. Development 130: 4919-4929. doi:10.1242/dev.00706. PubMed: 12930779.

17. Nakamura H, Kitagawa H, Kawaguchi Y, Tsuji H, Takano H et al. (1995) Intracortical facilitation and inhibition after paired magnetic stimulation in humans under anesthesia. Neurosci Lett 199: 155-157. doi:10.1016/0304-3940(95)12031-X. PubMed: 8584247.

18. Campos-Xavier AB, Martinet D, Bateman J, Belluoccio D, Rowley L et al. (2009) Mutations in the heparan-sulfate proteoglycan glypican 6 (GPC6) impair endochondral ossification and cause recessive omodysplasia. Am J Hum Genet 84: 760-770. doi:10.1016/j.ajhg.2009.05.002. PubMed: 19481194.

19. Fukushima M, Tsuchiyama Y, Nakato T, Yokoi T, Ikeda H et al. (1995) A female heterozygous patient with Fabry's disease with renal accumulation of trihexosylceramide detected with a monoclonal antibody. Am J Kidney Dis 26: 952-955. doi:10.1016/0272-6386(95)90061-6. PubMed: 7503071.

20. Flowers GP, Topczewska JM, Topczewski J (2012) A zebrafish Notum homolog specifically blocks the Wnt/beta-catenin signaling pathway. Development 139: 2416-2425. doi:10.1242/dev.063206. PubMed: 22669824.

21. Bökel C, Brand M (2013) Generation and interpretation of FGF morphogen gradients in vertebrates. Curr Opin Genet Dev, 23: 415–22. PubMed: 23669552.

22. Yu SR, Burkhardt M, Nowak M, Ries J, Petrásek Z et al. (2009) Fgf8 morphogen gradient forms by a source-sink mechanism with freely diffusing molecules. Nature 461: 533-536. doi:10.1038/nature08391. PubMed: 19741606.

23. Brand M, Granato M, Nüsslein-Volhard C (2002) Keeping and raising zebrafish. In: C. Nüsslein-VolhardR. Dahm. Zebrafish : A Practical Approach. Oxford: Oxford University Press. pp. 7-37.

24. Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF (1995) Stages of embryonic development of the zebrafish. Dev Dyn 203: 253-310. doi:10.1002/aja.1002030302. PubMed: 8589427.

25. Reifers F, Böhli H, Walsh EC, Crossley PH, Stainier DY et al. (1998) Fgf8 is mutated in zebrafish acerebellar (ace) mutants and is required for maintenance of midbrain-hindbrain boundary development and somitogenesis. Development 125: 2381-2395. PubMed: 9609821.

26. Wittbrodt J, Meyer A, Schartl M (1998) More genes in fish? Bioessays 20: 511-515. doi:10.1002/(SICI)1521-1878(199806)20:6.

27. Svensson G, Awad W, Håkansson M, Mani K, Logan DT (2012) Crystal structure of N-glycosylated human glypican-1 core protein: structure of two loops evolutionarily conserved in vertebrate glypican-1. J Biol Chem 287: 14040-14051. doi:10.1074/jbc.M111.322487. PubMed: 22351761.

28. Filmus J, Selleck SB (2001) Glypicans: proteoglycans with a surprise. J Clin Invest 108: 497-501. doi:10.1172/JCI13712. PubMed: 11518720.

29. Zhang L, David G, Esko JD (1995) Repetitive Ser-Gly sequences enhance heparan sulfate assembly in proteoglycans. J Biol Chem 270: 27127-27135. doi:10.1074/jbc.270.45.27127. PubMed: 7592967.

30. Saunders S, Paine-Saunders S, Lander AD (1997) Expression of the cell surface proteoglycan glypican-5 is developmentally regulated in kidney, limb, and brain. Dev Biol 190: 78-93. doi:10.1006/dbio.1997.8690. PubMed: 9331333.

31. De Cat B, Muyldermans SY, Coomans C, Degeest G, Vanderschueren B et al. (2003) Processing by proprotein convertases is required for glypican-3 modulation of cell survival, Wnt signaling, and gastrulation movements. J Cell Biol 163: 625-635. doi:10.1083/jcb.200302152. PubMed: 14610063.

32. Luxardi G, Galli A, Forlani S, Lawson K, Maina F et al. (2007) Glypicans are differentially expressed during patterning and neurogenesis of early mouse brain. Biochem Biophys Res Commun 352: 55-60. doi:10.1016/j.bbrc.2006.10.185. PubMed: 17107664.

33. LeClair EE, Mui SR, Huang A, Topczewska JM, Topczewski J (2009) Craniofacial skeletal defects of adult zebrafish Glypican 4 (knypek) mutants. Dev Dyn 238: 2550-2563. doi:10.1002/dvdy.22086. PubMed: 19777561.

34. Litwack ED, Ivins JK, Kumbasar A, Paine-Saunders S, Stipp CS et al. (1998) Expression of the heparan sulfate proteoglycan glypican-1 in the developing rodent. Dev Dyn 211: 72-87. doi:10.1002/(SICI)1097-0177(199801)211:1. PubMed: 9438425.

35. Iglesias BV, Centeno G, Pascuccelli H, Ward F, Peters MG et al. (2008) Expression pattern of glypican-3 (GPC3) during human embryonic and fetal development. Histol Histopathol 23: 1333-1340. PubMed: 18785116.