Skip to main content

Pancreatic and intestinal endocrine cells in zebrafish share common transcriptomic signatures and regulatory programmes



Endocrine cells of the zebrafish digestive system play an important role in regulating metabolism and include pancreatic endocrine cells (PECs) clustered in the islets of Langerhans and the enteroendocrine cells (EECs) scattered in the intestinal epithelium. Despite EECs and PECs are being located in distinct organs, their differentiation involves shared molecular mechanisms and transcription factors. However, their degree of relatedness remains unexplored. In this study, we investigated comprehensively the similarity of EECs and PECs by defining their transcriptomic landscape and comparing the regulatory programmes controlled by Pax6b, a key player in both EEC and PEC differentiations.


RNA sequencing was performed on EECs and PECs isolated from wild-type and pax6b mutant zebrafish. Data mining of wild-type zebrafish EEC data confirmed the expression of orthologues for most known mammalian EEC hormones, but also revealed the expression of three additional neuropeptide hormones (Proenkephalin-a, Calcitonin-a and Adcyap1a) not previously reported to be expressed by EECs in any species. Comparison of transcriptomes from EECs, PECs and other zebrafish tissues highlights a very close similarity between EECs and PECs, with more than 70% of genes being expressed in both endocrine cell types. Comparison of Pax6b-regulated genes in EECs and PECs revealed a significant overlap. pax6b loss-of-function does not affect the total number of EECs and PECs but instead disrupts the balance between endocrine cell subtypes, leading to an increase of ghrelin- and motilin-like-expressing cells in both the intestine and pancreas at the expense of other endocrine cells such as beta and delta cells in the pancreas and pyyb-expressing cells in the intestine. Finally, we show that the homeodomain of Pax6b is dispensable for its action in both EECs and PECs.


We have analysed the transcriptomic landscape of wild-type and pax6b mutant zebrafish EECs and PECs. Our study highlights the close relatedness of EECs and PECs at the transcriptomic and regulatory levels, supporting the hypothesis of a common phylogenetic origin and underscoring the potential implication of EECs in metabolic diseases such as type 2 diabetes.


The enteroendocrine cells (EECs) and pancreatic endocrine cells (PECs) play a crucial role in the control of metabolism of all vertebrates. The PECs, notably via the antagonistic action of glucagon and insulin, secreted respectively by the alpha and beta endocrine cells, regulate glycaemia. The EECs also participate in this control through the release of some hormones, like glucagon-like peptide-1 (GLP-1) or glucose-dependent insulinotropic polypeptide (GIP), which display incretin effect, enhancing the glucose-stimulated insulin secretion. The EECs secrete more than 15 hormones that modulate many physiological processes including the release of digestive enzymes and bile from the exocrine pancreas and gallbladder, control of intestinal motility, sensing of nutrients and microbial metabolites and stimulation or suppression of appetite [1]. The EECs were initially classified into distinct subtypes according to the hormone they secrete; however, subsequent analyses revealed that there is an extensive co-expression of several hormones in individual EEC [2]. Recent single-cell RNA-seq studies have indeed confirmed this co-expression in murine EECs thereby identifying at least 9 EEC subtypes expressing various combinations of neuropeptides or hormones [3, 4].

EECs display similarities with pancreatic endocrine cells (PECs). Some hormones are produced from both pancreatic and intestinal endocrine cells, like somatostatin (Sst), ghrelin (Grhl) or peptide YY (Pyy) [5, 6]. Also, many transcription factors are involved in the differentiation of both EECs and PECs such as NeuroD, Pax6, Isl1, Nkx2.2, Insm1 and Arx [7,8,9,10]. For example, Pax6 is crucial for the differentiation of PECs, and its inactivation disturbs the adequate proportion of the different endocrine cell subtypes. Indeed, in Pax6 KO mice [11, 12] and in pax6b−/− zebrafish [13], the number of alpha, beta and delta pancreatic cells is significantly reduced while epsilon cells are increased. In the small intestine, Pax6 has been reported to be required for the differentiation of some EEC subtypes [14, 15]; however, it is still not known which set of genes is regulated by Pax6 during EEC and PEC differentiation and whether the Pax6-regulated network is shared by these two tissues.

EECs have been detected in the gut of many vertebrate and invertebrate species, and their formation appears early during animal evolution dating back to the cnidarian-bilaterian ancestor [16,17,18]. In contrast, the pancreas and the PECs appeared much later with the first vertebrate species [19, 20]. The similarity between EECs and PECs could be due to several reasons. The first possibility could be a common phylogenetic origin; the PECs could have indeed derived from some EECs and moved out from the gut to form clusters and to secrete other hormones like insulin and glucagon. This scenario is compatible with the “sister cell-type model” [21] where two cell types derived from an ancestral cell type by segregation of function (e.g. secretion of distinct hormones) due to the gain or loss of transcription factors. In this model, the two sister cell types still share many characteristics and can display similarities in their transcriptomes [22]. Another reason of the PEC and EEC similarities could be the co-option of the EEC (or neuronal) regulatory pathways by the pancreatic cells [23]. To approach this matter, we have measured in the present study the degree of similarity between EECs and PECs at the genome-wide expression level and determined whether they share regulatory programmes. We previously determined the transcriptomic landscape of zebrafish PECs [24]. In contrast, the EECs from zebrafish have not yet been fully characterized, and only a few hormones have been shown be expressed in these cells [25,26,27,28,29]. The first objective of this study was to perform a detailed characterization of zebrafish EECs by determining their transcriptomic profiles and identifying all neuropeptide hormones and regulatory genes expressed in these cells. Secondly, the transcriptomic landscape of zebrafish EECs was compared with those of zebrafish PECs to measure the degree of similarity. Finally, we analysed the consequences of pax6b inactivation on the transcriptome of both pancreatic and intestinal endocrine cells and show that Pax6b controls a large set of genes identical in both EECs and PECs further supporting their close relatedness. Finally, we discuss our findings in light of the evolution of EECs and pancreatic cells, supporting the model of a common phylogenetic origin and bringing some implications for deciphering the defects involved in type 2 diabetes.


Characterization of the enteroendocrine cells in zebrafish

In mice, at least 15 different hormones have been described to be expressed in the gut [1, 8]. Out of these, we identified in the zebrafish genome at least one orthologue for 11 of them, 9 showing syntenic conservation at their genomic loci with the murine/human genes, thereby validating orthology relationship (Table 1). We did not find any orthologous genes coding for gastrin (gast), secretin (sct), the pancreatic polypeptide (pp) nor motilin (mln); however, a motilin-like (mlnl) gene, proposed to be the functional equivalent of motilin, is present in the zebrafish genome [31]. In order to identify all peptide hormones expressed in the gut of zebrafish larva, we performed RNA-seq on isolated EECs. This was achieved by microdissection of the intestines from 4-dpf zebrafish transgenic Tg(pax6b:GFP)ulg515tg larvae, which express GFP in the EECs [32], followed by cell dissociation and selection of GFP+ EECs by FACS. The restricted expression of Pax6b:GFP in the EEC and not in the enteric neurons (EN) was confirmed by immunohistochemistry showing a total absence of colocalisation between the Hu enteric marker and pax6b:GFP (Additional file 1: Fig. S1). RNA-seq was performed on four independent EEC preparations, and about 40 million paired-end reads were obtained from each Illumina library, 70–80% of which mapped to the zebrafish genome. We provide the expression level of all genes in counts per million (CPM) and in fragments per kilobase of exon model per million reads mapped (FPKM) in additional file 2: Table S1. This expression profiling allowed us to determine that, amongst the 12 zebrafish genes corresponding to EEC hormones, 11 are expressed in the zebrafish EECs at significant levels (above 100 CPM or 100 FPKM), neurotensin (nts) being the only hormone not expressed in EECs (Table 1). The expression of many of these hormones was further validated by whole-mount in situ hybridization (WISH) thereby locating the region of the gut displaying the strongest expression and defining the onset of their expression (Fig. 1 and Table 1). The number of EECs detected in the gut of 4-dpf zebrafish larvae by WISH varies for each hormone; for example, more than 20 cells expressing peptide YY-b (pyyb), cholecystokinin-a (ccka), preproglucagon-a (gcga), galanin (galn) or insulin-like 5a (insl5a) transcripts were detected (Fig. 1a, b, g–i; Table 1 and Additional file 3: Fig. S2) while only few cells express ghrelin (ghrl), glucose-dependent insulinotropic polypeptide (gip) and peptide YY-a (pyya) (Fig. 1d–f and Table 1). These differences of expression are also reflected in the RNA-seq data: gcga, ppyb, ccka, and insl5a reaching more than 2000 CPM while ghrl and gip range around 100 CPM (Table 1).

Table 1 Identification and expression level of zebrafish enteroendocrine hormones
Fig. 1.

Expression of different enteroendocrine hormones in the zebrafish gut. Whole-mount in situ hybridization (WISH) obtained with different hormonal probes at 3 dpf (ah, k) or 4 dpf (i, k). Ventral views of embryos with anterior on the left. The dotted lines represent the location of the gut. The probes used are as follows: adcyap1a, adenylate cyclase-activating polypeptide 1a; ccka, cholecystokinin a; gcga, Glucagon a; pyya/b, Peptide YYa/b; ghrl, Ghrelin; gip, glucose-dependent insulinotropic polypeptide; insl5a, insulin-like peptide 5a; calca, calcitonin. Scale bars, 50 μm. il Fluorescent in situ hybridization (FISH) (DAPI staining shown in grey)

Data mining of the zebrafish EEC RNA-seq with the GO term “hormone” also revealed a strong expression of three neuropeptide transcripts not described so far as expressed in mammalian EECs: adenylate cyclase-activating polypeptide 1a (adcyap1a) (previously named PACAP), proenkephalin-a (penka) and calcitonin (calca) (Table 1). Fluorescent in situ hybridization experiments (FISH) confirmed the expression of these transcripts in the zebrafish EECs (Fig. 1j–l). An average of 68 penka+ and 12 calcitonin+ EECs were detected in the bulb and the mid-gut and more than 80 cells expressing adcyap1a were found all over the gut (Table 1 and Additional file 3: Fig. S2). As enteric neurons (ENs) surround the gut at that stage and express a variety of neuropeptides, some being expressed also by EEC, we next determined whether some of the hormone expressed by EEC would be also detected in ENs. For that purpose, we retrieve the transcriptomic data of zebrafish ENs, recently determined by RNA-seq after sorting the ENs by FACS from Tg(phox2b::GFP) larvae at 7 dpf [30]. Comparison of the EN and EEC RNA-seq data indicates that most EEC hormones display a much higher expression in EECs compared to ENs, except for sst1.2, sst2 and vipb (see Table 1). FISH using the EN-specific marker Phox2b confirms the strong expression all over the gut of vipb in ENs while the vipb+ EECs are mostly located in the posterior gut (Additional file 4: Fig. S3A and B). Interestingly, RNA-seq data indicate that adcyap1a is exclusively expressed by EECs while adcyap1b is restricted to ENs. FISH confirms the strong expression of adcyap1a in EECs and not in phox2b+ ENs (Additional file 4: Fig. S3 C). However, IHC using an anti-ADCYAP1 antibody reveals staining in pax6b:GFP+ EECs as well as at the level of EN axon fibres in larvae, as previously shown [33] (Fig. 2b). Similar results were observed in the adult intestine (Fig. 2d). Taken together, all these data confirm the expression of the Adcyap1a in EECs and of Adcyap1b in ENs. As for the two other new hormones Calca and Penka, RNA-seq comparison indicates that they are exclusively expressed by EECs, and this was confirmed for Penka for which we observe strong enkephalin immunostainings in EECs but not in ENs at the larva stage as well as in adults (Fig. 2a, c).

Fig. 2.

Immunostaining of enkephalin and Adcyap1 in the zebrafish intestine. a, b Confocal images of the intestines from 5-dpf Tg(pax6b:GFP) zebrafish larva stained with antibodies directed against GFP (green), Hu (white), ENKEPHALIN (a; red) or ADCYAP1 (b; red). c, d Confocal images of the intestine of adult Tg(8.5nkx2.2:GFP) stained with antibodies directed against GFP (green), ENKEPHALIN (c; red) or ADCYAP1 (d; red). Arrows show the location of some EECs expressing Adcyap1 (b, d) or enkephalin (c). Adcyap1 immunostaining is also observed at the level of EN axon fibres in gut larva (b) and in the muscularis layer of the adult intestine (shown by an asterisk in d). The zebrafish Adcyap1a and Adcyap1b display respectively 6 and 3 mismatches with human ADCYAP1 at the level of the epitope reacting with the Adcyap1 antibody, explaining the stronger signal of Adcyap1b in EN fibres. (The nkx2.2:GFP transgene is expressed specifically in EECs like Pax6b). Scale bar, 20 μm

Close relatedness between the enteroendocrine and pancreatic endocrine cells

Although EECs and PECs are known to share common features, we wanted to evaluate the degree of similarity at a more global level. Thus, we compared the RNA-seq data of the zebrafish PECs and EECs as well as of other zebrafish cell types and organs used as control (Fig. 3a). The heatmap shows that all endocrine cells cluster together apart from other pancreatic or intestinal cells and other tissues. Indeed, the EEC cluster is much closer to the PEC cluster than to the ductal or the acinar pancreatic cells or even to the whole intestine tissue which is mainly composed of enterocytes and of only 1–2% of EECs. In agreement with previous reports showing shared features between pancreatic cells and neurons [23, 34], this clustering analysis reveals also transcriptomic similarities of brain tissue not only with pancreatic but also with intestinal endocrine cells.

Fig. 3.

Similar gene expression profiles between zebrafish EECs and PECs. Global comparison of transcriptomic profiles (RNA-seq) from distinct zebrafish tissues/organs. a Clustered heatmap displaying the Euclidean distance matrix between every pair of zebrafish RNA-seq datasets. The order of the different tissues is identical for the two axes (rows and columns) as shown by the different colours on the upper part and left part of the matrix. Each zebrafish tissue (analysed in triplicate) was compared with the other samples as presented in the matrix. Darker colour indicates closer distance (i.e. more similar transcriptomes). EEC transcriptome is mostly similar to the transcriptome of embryonic pancreatic endocrine cells (PECs) or to the adult alpha, beta or delta pancreatic cells. The accession number of all used RNA-seq datasets is given in the “Methods” section. b Venn diagram showing the overlap (74%) of genes expressed in both EECs and PECs above the threshold of 10 CPM. A total of 8264 genes and 8326 genes were expressed above 10 CPM (counts per million) respectively in EECs and in PECs

The close relatedness of EECs and PECs is also illustrated by the high percentage of genes expressed in both tissues: amongst all genes expressed in pancreatic or intestinal endocrine cells above the threshold of 10 normalized CPM (8326 and 8264 genes, respectively), 74% of them are expressed in both tissues (Fig. 3b), while if the same analysis is done for genes expressed in EECs and liver, only 45% of them are expressed in both tissues (data not shown). Similar percentages are obtained when considering genes coding for transcription factors (TF): 72% of all expressed TF-coding genes are detected in both PECs and EECs, while only 45% are expressed in both pancreatic endocrine cells and liver. The list of 481 TF expressed in both PECs and EECs (given in Additional file 5: Table S2) includes notably most of the transcription factors reported to be important for zebrafish PEC differentiation (Table 2). For example, arxa and pax4, known to be an important determinant of cell fate of pancreatic endocrine cells [35, 36], are also expressed in the zebrafish gut, and like in pancreas, their expression is mainly non-overlapping (Fig. 4a). Fev, known to be highly expressed in pancreatic endocrine cells [37], is also amongst of the most highly expressed transcription factors in the zebrafish EECs (Table 2) and is detected by FISH in many cells all over the gut (Fig. 4b). Similarly, insm1b [9] is also highly expressed in zebrafish EECs (Fig. 4c, d and Table 2) in addition to its expression in enteric neurons (labelled with asterisk in Fig. 4d) as described in mice. Expression of pdx1, visualized using the tg(BAC pdx1:GFP) [38], is also detected in scattered EECs of the mid- and posterior gut (Fig. 4e, g), in addition to a widespread expression in cells of the rostral part of the intestine as previously described [39] (Fig. 4e, f). The expression of pdx1 in EECs of the mid- and posterior gut was confirmed by WISH (Fig. 4h) and by immunohistochemistry (Fig. 4i). As expected, ascl1a was barely detected in the transcriptome of 4-dpf EECs (Table 2) as we have previously shown that its expression in the secretory precursor cells turns off as soon as the cells pursue their differentiation process [27]. In the same way, sox4b, coexpressed with ascl1a in the secretory precursor cells, is hardly detected in the transcriptome of the EECs at 4 dpf (Table 2), indicating that sox4b is also switched off with the maturation of the EECs. In contrast, Nkx6.1, Nkx6.2 and Mnx1 are the only 3 transcription factors involved in zebrafish PEC differentiation that are not expressed in EECs. Many of these EEC/PEC transcription factors are expressed at lower levels in enteric neurons as revealed by the EN RNA-seq data (Table 2) with the exception of insm1b, neurod1 and ascl1a which are expressed above 50 CPM.

Table 2 Expression level in EECs of transcription factors involved in zebrafish PEC differentiation
Fig. 4.

Expression profiles of pax4, arx, pdx1 and fev in the zebrafish gut. Fluorescent in situ hybridization (FISH) performed on 66 hpf (a), 72 hpf (bd) or 96 hpf (ei) embryos. Views of the gut (delimited by dotted lines; anterior to the left) showing the non-overlapping expression of pax4 (green) and arxa (red) (a), numerous EECs expressing fev (b) and insm1b (d). Higher magnification showing insm1+ EECs located within the gut epithelium and insm1+ enteric neurons (shown by asterisks) outside but juxtaposed to the epithelium. e View of the gut from the tg(BAC pdx1:GFP) larvae (GFP immunostaining) with higher magnification of the anterior (f) and posterior (g) parts. FISH using pdx1 probe (h) and immunostaining using Pdx1 antibody. P, pancreas. Scale bars, 20 μM

The number of genes which are selectively expressed in PECs or in EECs is relatively low (gene lists available in Additional file 6: Table S3). Gene Ontology (GO) enrichment analyses of the PEC-specific genes highlight terms such as “cell adhesion” due to the enriched expression of many cadherins and other cell adhesion molecules. This enrichment makes senses as the PECs are clustered into islets while EECs are scattered in the gut. The EEC-specific genes display an enrichment of GO terms including “drug transmembrane transport”, “ABC transporter” or “oxidation process” due to the enriched expression of many transmembrane transporters and cytochrome P450 type molecules (Additional file 6: Table S3).

In conclusion, all these data indicate a close relatedness between PECs and EECs at the level of their transcriptomes and of their regulatory factors, suggesting the involvement of similar regulatory programmes in these two cell types.

Shared Pax6-dependent regulatory programmes between pancreatic and intestinal endocrine cells

To investigate whether pancreatic and intestinal endocrine cells are controlled by similar regulatory programmes, we analysed the effect of Pax6 inactivation on the transcriptome of both cell types. Zebrafish has two pax6 paralogs; however, only pax6b is expressed in PECs and EECs. To isolate pax6b mutant EECs and PECs, pax6bsa0086 heterozygous fish harbouring the Tg(pax6b:GFP) transgene were inbred and pax6bsa0086 homozygous embryos were selected based on lens abnormalities [13]. The expression of the transgene Tg(pax6b:GFP) was not perturbed in the homozygous mutants (data not shown), allowing us to isolate by FACS the PECs or EECs from pax6b−/− mutant larvae after microdissection of the dorsal pancreatic bud or intestine, respectively. RNA-seq was performed on three independent preparations of EECs or PECs purified from pax6bsa0086 mutant and from wild-type zebrafish (see the “Methods” section). The principal component analysis shows a tight clustering of replicates underscoring a good reproducibility of the data and showing distinct clusters for mutant and wild-type cells for both PECs and EECs (Fig. 5a). Differential expression analyses identified 2824 and 1634 Pax6b-regulated genes in PECs and in EECs, respectively (FDR < 0.1) (Fig. 5b, d) (expression of all genes and of Pax6b-regulated genes are given respectively in Additional files 7 and 8: Table S4 and S5). Comparison of these two sets of genes reveals that there is a large set of 517 genes which are altered by pax6b inactivation in both organs (Fig. 5c); this overlap is statistically significant according to Fisher’s exact test (p value < 2.2e−16). Furthermore, most of these genes (424 genes) are regulated in the same way in the pancreas and in the intestine (215 upregulated and 209 downregulated genes) (Fig. 5c), indicating that a part of the Pax6-dependent gene regulatory network is thus identical in EECs and PECs. To explore the pathways and biological processes regulated by Pax6b, we performed a Gene Ontology analysis on these 3 sets of Pax6b-regulated genes: the common set for PECs and EECs (517 genes) or the specific one for PEC (2307 genes) or EECs (1117 genes) (see Additional file 8: Table S5). GO terms associated to “cell-cell signalling”, “regulation of appetite” or “hormone activity” were enriched significantly in the common set as several hormones expressed in PECs and EECs are affected in the pax6 mutant (e.g. pyya, sst2, pyyb, see below). Such an effect results probably from a change in the expression of several transcription factors (e.g. pdx1, neurod1, rfx6, nkx2.2, see below) as revealed in the GO lists “endocrine pancreas development” and “sequence-specific DNA binding”. Also, the GO analysis highlights the role of Pax6b in the control of ion transport (K+, Ca++ and Cl−) as well as calcium signalling. Interestingly, “ion transport” and “calcium channel/binding activity” were also GO-enriched terms amongst the genes regulated by Pax6b only in PECs or in EECs, supporting an important role of Pax6b for controlling these pathways.

Fig. 5.

RNA-seq analysis of EECs and PECs from wild-type and pax6b mutants identifies a set of genes displaying similar Pax6b regulation in both the pancreas and intestine. a, Principal component analysis (PCA) on all EEC and PEC RNA-seq data obtained from wild-type and pax6bsa0086 mutants. The close clustering of the triplicate wild-type and mutant samples demonstrates the high reproducibility of the data. b, d MA plots showing all the upregulated (in blue) and downregulated (in red) genes in the pax6bsa0086 mutant PECs (b) and EECs (d). c Venn diagrams showing the overlap of genes regulated in both EECs and PECs; amongst the 517 regulated genes, 215 and 209 genes are respectively up- and downregulated in both cell types as shown in the blue and red Venn diagrams. e, f Expression ratio (in log2 of fold change from pax6b mutant versus wild-type) of hormones expressed in PECs (e) and in EECs (f); downregulated hormones are in red, upregulated hormones are in blue and not statistically are affected in grey; RNA-seq values are shown in Tables 3 and S4

Drastic increase of ghrelin/motilin-like-expressing cells in pax6b−/− pancreatic and intestinal endocrine cells with a concomitant decrease of other endocrine cell subtypes

The above analyses indicate a regulation by Pax6 of the expression of several hormones in PECs and EECs. Amongst them, ghrelin (ghrl) and motilin-like (mlnl) genes are strongly upregulated in both the pancreas and intestine of pax6b mutants (Fig. 5e, f and Table 3). FISH experiments confirmed an increase in the number of ghrl+ cells and mlnl+ cells in both the pancreas and intestine (Fig. 6A–H). Furthermore, double fluorescent staining demonstrated that these two hormones are often co-expressed in the same endocrine cells (Fig. 6I, J). In both tissues, the increase of ghrl and mlnl gene expression is concomitant to a decrease of other hormones: galn, pyyb and calca genes are indeed significantly downregulated in EECs of pax6bsa0086 mutants, while in the pancreas, the expression of notably the insulin (ins) and somatostatin-1.1, 1.2 and 2 genes is strongly decreased (Fig. 5e and Table 3). To determine whether these modifications in hormone expression result from a change in cell subtype proportion in pax6b mutant PECs, we analysed the expression of markers of these different pancreatic cell subtypes that we and others have previously identified (Additional file 9: Table S6) [24, 40]. We found that about 40% of genes showing an enriched expression in beta cell (with a fold enrichment > 4-fold) were downregulated in the islet of pax6b−/− embryos, while only 2% were upregulated (Fig. 7, upper panel), strongly suggesting a loss of beta cells. This loss was further confirmed by FISH experiments: firstly, ndufa4l2, a gene showing a restricted expression in zebrafish pancreatic beta cells, was not detected anymore in the pax6b mutants (Fig. 7a–d). Similarly, the expression of the two beta cell-enriched transcription factors pdx1 and nkx6.2 were strongly decreased in PECs of pax6b mutants while could still be detected at low levels in the ventral pancreatic progenitors (Fig. 7e–h, Table 3). All these data indicate the loss of beta cells in pax6b−/−sa0086 embryos. Concerning the pancreatic delta cells, 35% of delta-enriched genes were downregulated in the pax6bsa0086 mutants and 12% were upregulated (Fig. 7, upper panel). The decreased number of delta cells in pax6b mutants was further supported by FISH through the strong reduction of cells expressing somatostatin2 (sst2) as well as the loss of the delta-specific markers laminin C2 (lamc2) and the hormone peptide Neuromedin Ba (nmba) (Fig. 7i–l, Table 3). The two delta-specific transcription factors Cdx4 and Hhex were also downregulated in PECs of pax6b mutants (Table 3). Inversely, 40% of genes with enriched expression in pancreatic epsilon cells are upregulated in pax6b mutants while only 2% are decreased (Fig. 7, upper panel). Like ghrelin and motilin-like, mboat4, another epsilon-specific marker, is strongly stimulated following pax6b inactivation as revealed by FISH (Fig. 7m, n) and by the RNA-seq data (Table 3). Thus, these data also indicate an increase of epsilon cells in pax6bsa0086 mutant embryos. Finally, the alpha cells in pax6b mutants display an “ambiguous” phenotype with a respective decrease and increase of 26% and 12% of alpha enriched genes (Fig. 7, upper panel). The RNA-seq data indeed reveal a strong decrease of the alpha cell markers gcgb and pnoca while others like gcga, arx and scinlb are even increased (Table 3). FISH confirmed these results. Indeed, although gcga and gcgb are detected in the same alpha cells in wild-type embryos, gcgb expression was severely reduced in the pax6b mutants in contrast to gcga that remained highly expressed (Fig. 7i, j). Like gcgb, the alpha cell marker pnoca (prepronociceptin a) was drastically downregulated (Fig. 7q, r). On the other hand, the alpha cell marker scindl (Fig. 7u, v) and the alpha cell determinant arxa (Fig. 7s, t) were still expressed. Double labelling of arxa and scinlb with ghrelin revealed that these two markers are actually expressed in both alpha and epsilon cells in zebrafish (Fig. 7s–x). Regarding the pan-endocrine markers, which are not enriched in a specific endocrine cell subtype, only a minority (about 15%) are either down- or upregulated in the pax6b−/− mutants (Fig. 7, upper panel). Altogether, our data indicate that there is a loss of beta cells and a strong reduction of delta cells in the pax6b−/− mutants with a concomitant increase in the number of epsilon cells. Alpha cells are still generated in the pax6bsa0086 mutants, but these cells are misdifferentiated.

Table 3 Expression level of some selected markers in EEC and PEC isolated from wild-type and pax6bsa0086 mutant larvae
Fig. 6.

Increase in the number of cells co-expressing ghrl and mlnl in the pax6b−/− PECs and EECs. FISH using ghrelin (ghrl, labelled in red) and motilin-like (mlnl, labelled in green) probes on 2.5 dpf (A, B, E, F and I) and 4.5 dpf (C, D, G, H, J) zebrafish larvae. Confocal views of the pancreatic islet (A, B, E, F) and of the gut (C, D, G, H). I, J Overlays of ghrl (I′ and J′) and mlnl (I″ and J″) stainings demonstrating co-labelling with the two probes. A, E, I, I′, and I″ are views of the pancreatic islet from a wild-type larva while B and F are the pancreatic islet from pax6b−/−. C, G Views of the intestines from wild-type larvae; D, H, J, J′, J″ Views of the gut from a pax6b−/− larvae. Scale bars: 20 microns for A-F and I panels, 50 microns for G, H and J panels

Fig. 7.

pax6b inactivation leads to a loss of pancreatic beta cells, reduction of delta cells, increase of epsilon cells and abnormal alpha cells. (Upper panel) percentage of markers from beta-, delta-, alpha-, epsilon- and pan-endocrine cells which are downregulated (in red) or upregulated (in blue) in the pax6b−/− PECs. The list of markers for each cell subtype is shown in Additional file 9: Table S6. The up- or downregulation was determined from RNA-seq data from pax6bsa0086 mutant versus wild-type. (Lower panel) Expression analysis of cell subtype markers by WISH (a, b) and FISH (cx) on 2 dpf zebrafish pax6b mutant and wild-type embryos as noted on the top of the panel; probes are indicated on the left of each picture. Scale bars: 100 µm for panels a and b and 20µm for panels c to x

To verify that the increase of ghrl+/mlnl+ endocrine cells in the intestine of pax6b−/− larvae (Fig. 6) is also concomitant to a decrease of other EECs, we analysed by FISH the expression of pyyb, a gene highly expressed in the intestine and strongly downregulated in pax6b mutants (Table 3). The number of EECs expressing pyyb was indeed drastically reduced in the intestine of pax6bsa0086 mutants (Fig. 8d, e). Thus, all these FISH confirms that pax6b has a similar role in the control of the proportion of the various endocrine cell types in both intestine and pancreas.

Fig. 8.

The effect of Pax6b on the number of EEC expressing pyyb and mlnl does not depend on its homeodomain. FISH showing the respective increase of mlnl+ EECs and decrease of pyyb+ EECs in the pax6bsa0086 null mutants (b, e), while the hypomorphic pax6bsunrise mutants harbouring a mutation in pax6b homeobox do not display modifications in the expression of EEC hormones (c, f)

Like in the pancreas, the homeodomain of Pax6b is not crucial for the differentiation of the endocrine cells of the intestine

Previous studies have shown that a missense mutation in the homeobox of pax6b (sunrise allele) gene causes eye developmental defects but does not perturb the development of pancreatic endocrine cells [13, 41], indicating that the homeodomain of Pax6b is not crucial for endocrine cell differentiation in the pancreas in contrast to eye development. To get a clue whether the mode of action of Pax6b is also similar between EECs and PECs, we next investigated whether the EEC differentiation was affected in the sunrise mutants. FISH revealed no effect of the sunrise mutation on the expression of pyyb and of mlnl in contrast to the null pax6bsa0086 mutation (Fig. 8c, f). This indicates that the homeodomain of Pax6b is not crucial for the differentiation of endocrine cells both in the intestine and in the pancreas.


By RNA sequencing of FACS-purified cells, we report in this study the transcriptome of enteroendocrine cells (EECs) from zebrafish highlighting the repertoire of peptide hormones and regulatory factors expressed by these cells. The comparison of these RNA-seq data with other zebrafish tissues confirms the close similarity of EECs with pancreatic endocrine cells (PECs). Furthermore, we show that Pax6b is required for the proper expression of many hormones in both the pancreas and intestine, acting via a large set of common gene targets. Inactivation of zebrafish Pax6b does not affect the total number of endocrine cells but affects the balance between endocrine cell subtypes, leading to an increase of ghrelin+/motilin-like+ cells and a reduction of other endocrine cell subtypes in both tissues. The Pax6b homeodomain is also dispensable in the two organs. All these data highlight the high similarity between EECs and PECs.

The present study reveals also that the zebrafish EECs express most of the hormones known to be secreted by mammalian EECs. Only 4 EEC hormones were not detected in zebrafish. Secretin (Sct), pancreatic polypeptide (PP) and gastrin could not be detected as no orthologous gene is present in the zebrafish genome. The absence of sct and pp genes is in agreement with previous studies indicating the lack of these two genes in teleost [42,43,44]. The detection of PP by immunofluorescence at the level of zebrafish EECs [29] and zebrafish PECs [45] is probably due to cross-reactions of the antibodies with Pyya and/or Pyyb as significant sequence similarities are observed between Pyy and PP mature peptides (see Additional file 10: Fig. S4). Although the neurotensin gene is present in the zebrafish genome, no expression was detected in EECs, at least in 4 dpf larvae. On the other hand, the present transcriptomic profiling of zebrafish EECs reveals the identification of three neuropeptide transcripts not described so far in EECs: proenkephalin-a (penka), calcitonin (calca) and adenylate cyclase-activating polyeptide 1a (adcyap1a). We could confirm at the protein level the expression of two of them, enkephalin and Adcyap1 (Fig. 2). Further studies should be performed to confirm the expression of Calca at the protein level and determine the role of these 3 neuropeptides on the digestive tract. To get a hint whether these hormones could be also expressed in the mammalian intestine, we searched into the NCBI GEO repository and into the Human Protein Atlas database [46] for the presence of these transcripts in the gut. We found that Penk is expressed in the colon of 8-week-old male mice [47] and at a low level in the human colon. Similarly, human CALCA and more specially CALCB (calcitonin-related polypeptide B) transcripts can be detected in bulk RNA-seq of the gastrointestinal tract mainly at the level of the colon. Finally, ADCYAP1 expression is detected in the human gastrointestinal tract mostly in the colon as well as in beta pancreatic cells [48]. Furthermore, the recent single-cell RNA-seq data from the small murine intestine reveal that some enterochromaffin cells express Adcyap1 at low levels [3]. Thus, all these data suggest that the three hormones Penk, Calca(/b) and Adcyap1 could also be expressed in some mammalian EECs, and their expression should be verified by immunostaining notably at the level of the colon.

Some similarity between EECs and PECs has been previously recognized notably thanks to the shared expression of markers including some hormones and transcription factors. Using the zebrafish, the present study demonstrates the similitude at the transcriptomic level and by identifying a common Pax6b-dependant gene regulatory network. When considering the ontogeny and evolutionary aspects of EECs and PECs, it is tempting to speculate that these two endocrine cell types not only share developmental pathways but also have a common ancestor. Indeed, on the ontogeny aspect, both EECs and PECs derive from the endoderm layer and their specification is controlled by lateral inhibition via the Delta/Notch signalling pathway, leading to the formation of scattered endocrine cells within the epithelium. Inhibition of the Notch pathway allows the activation of proneuronal bHLH factors in endocrine precursors that drive the formation of EECs and PECs [49, 50]. As proposed by others [16], the EECs have probably co-opted for these neuronal gene programmes to trigger their differentiation early during the evolution of metazoa. While the EECs are present in the gut of invertebrates, the first signs of the pancreatic organ appear much later with the first vertebrates: the agnathan species, such as lampreys or hagfishes, possess at the adult stage one or two aggregates of insulin-expressing cells just next to the intestine and near the bile duct without associated exocrine tissue while other endocrine cells expressing glucagon, somatostatin or PYY-like peptides are present within the intestinal epithelium [19, 20, 51]. Furthermore, at the larva stage, some of these insulin+ cells are still detected within the intestinal epithelium of lampreys [52]. Thus, these observations suggest that, when the first vertebrate species appeared, some enteroendocrine cells could move out from the gut to form clusters, possibly through the upregulation of cell adhesion molecules that have led to the formation of pancreatic endocrine islets. This hypothesis is supported by studies on sea urchin larvae showing the presence of cells expressing insulin-like peptide in the gut region located at the stomach-intestine boundary and expressing the pdx1 orthologous gene [18, 53]. Furthermore, the generation of the first PECs in the zebrafish embryos is also consistent with this idea: indeed, the first insulin-expressing cells appear scattered within the anterior endodermal gut epithelium at an early stage (14 hpf) and subsequently migrate dorsally from the embryonic gut to form a cluster and generate the principal pancreatic islet (named the dorsal pancreatic bud) [45]. Thus, all these observations entail the notion that PECs are in fact EEC-like cells, explaining why PECs share so many features with EECs. This concept is also supported by previous reports showing the conversion of EECs into insulin-expressing beta-like cells upon targeted expression or inactivation of some transcription factors [54,55,56]. The close relatedness between EECs and PECs underlined in our study has also important implications for the interpretation of GWAS for type 2 diabetes (T2D). Indeed, as many of the genomic variants linked to T2D are near (or in) genes expressed in PECs [57, 58], it is widely considered that at least part of the defects causing T2D stems from PEC dysfunctions. However, as the transcriptomes of PECs and EECs are highly similar, the defects causing T2D may also have their origin from EECs.

The present study confirms that Pax6 is required for the proper expression of many hormones from EECs and PECs and changes the proportion of the endocrine cell subtype without affecting the total number of endocrine cells in both the pancreas and intestine. The pancreatic phenotype of zebrafish pax6b-null mutants is comparable to the Pax6 mutant mice as the same major pancreatic hormones are up- and downregulated in the two species, although the effect on beta cells looks more drastic in zebrafish than in mice (Pax6 KO mice displaying only 70% reduction) [11, 59]. The almost complete loss of beta cells in zebrafish may derive from a stronger decrease in the expression of the beta-specific factor Pdx1 (as well as Mnx1 and Nkx6.2) in zebrafish pax6b mutants. On the other hand, in the intestine, there is a striking contrast between zebrafish and mice Pax6 mutants: while Pax6 is required for the expression of several EEC hormones in both species, the identity of these hormones is strikingly divergent. Indeed, amongst the 8 EEC hormones analysed in Pax6 mutant mice, Gcg (GLP1 and 2), Sst and Gip genes were strongly downregulated while Pyy was not significantly affected (ghrelin and motilin were not analysed in Pax6 KO murine intestine). In zebrafish, pax6b inactivation causes a drastic downregulation of pyyb and does not affect significantly gcga, sst2 and gip gene expression. The reason for such differences between mice and zebrafish is unclear and probably rely on changes during vertebrate evolution either at the level of genomic Pax6 binding sites or of downstream regulatory programmes. Further studies including ChIP-seq experiments will be required to answer that question.


By determining the transcriptomic landscape of zebrafish enteroendocrine cells, we report in this study all peptide hormones and regulatory factors expressed in these cells. The close relatedness of EECs and PECs observed at the transcriptomic level was further supported by analysing the effect of pax6b loss-of-function revealing that a large set of genes is regulated similarly by Pax6b in both cell types. These results support the hypothesis of a common phylogenetic origin of EECs and PECs; they also underscore the potential implication of EECs in metabolic diseases such as type 2 diabetes.


Zebrafish maintenance, transgenic and mutant lines and isolation of EECs and PECs by FACS

Zebrafish (Danio rerio) were raised according to standard protocols and staged according to Kimmel et al. [60]. The following zebrafish transgenic and mutant lines were used: Tg(pax6b:GFP)ulg515 [61], TgBAC(pdx1:EGFP)bns13 [38], Tg(− 8.5nkx2.2a:GFP)ia2 [62], pax6bsa0086 and pax6bsunrise [13]. Enteroendocrine cells (EECs) were isolated by dissecting the gut from about 200 Tg(pax6b:GFP)ulg515 larvae at 4 dpf, taking care of not including pancreatic tissue. Cell dissociation was next performed by incubation in HBSS 1× supplemented with 100 U/ml collagenase IV and 0.3 U/ml Dispase (Life Technologies) for 10 min. Cells were washed in HBSS (Mg2+ and Ca2+ free) containing 1% BSA, and GFP-expressing EECs were selected by two consecutive FACS purifications, the first in the “yield mode” and the second in the “purity mode”, using FACS Aria II. Four replicates of EEC containing about 3000 cells were prepared. Pancreatic endocrine cells (PECs) were also obtained from the Tg(pax6b:GFP)ulg515 line [61] by dissecting the dorsal pancreatic bud from about 200 27-hpf transgenic embryos. FACS selection was performed as described for EECs except that cell dissociation was performed in Tryple Select 1X (Gibco) supplemented with 100 U/ml collagenase IV (Life Technologies) for 5 min. For the preparations of EECs and PECs from pax6b null mutant embryos, the pax6sa0086 line [13] was first crossed with the Tg(pax6b:GFP)ulg515 line; heterozygous pax6bsa0086 fish harbouring the transgene (pax6b:GFP) were inbred to generate homozygous pax6bsa0086 transgenic embryos which were selected based on the absence or reduction of lens. The isolation of EECs and PECs from pax6bsa0086 homozygotes were performed in triplicates following the same procedure than for the wild-type larvae. The accuracy of pax6bsa0086 homozygote selection was verified after the RNA-seq by checking the presence of the null sa0086 allele in 100% of pax6b reads in the mutant samples.

cDNA synthesis, library preparation and sequencing

Each EEC or PEC sample obtained after FACS was directly pelleted by centrifugation and resuspended in 3.5 μl of reaction buffer, lysed by freezing in liquid nitrogen and stored at − 80 °C according to the Smart-seq2 protocol [63]. cDNA was synthesized and amplified by a 13-cycle PCR reaction. The quality of cDNA was verified by 2100 High Sensitivity DNA assay (Agilent Technologies); 1-ng cDNA was used for preparing each cDNA library using the Nextera-XT kit (Illumina) and sequenced on Hi-seq 2000 to obtain around 40–60 millions of reads (100 base paired-ends).

RNA-seq data analyses

Sequences were trimmed in order to remove adaptors and low-quality bases. Trimmed reads were mapped into the zebrafish genome GRCz11 (Ensembl Release 92, using the STAR software v.2.5.4b [64]. Gene expression was measured from the mapped reads by using built-in STAR module (--quantMode GeneCounts) and are expressed in counts of reads per million (CPM) [65]. The RNA-seq raw data have been deposited in the Gene Expression Omnibus (GEO) under the accession number GSE149081. The comparison of the transcriptome of EECs and PECs with other zebrafish tissues was performed using the DESeq2 R package [66]. This comparison includes RNA-seq data from the heart (ArrayExpress: E-MTAB-460; GEO: GSE71755), brain (ArrayExpress: E-MTAB-460), pancreatic cells [24, 67], liver (GEO: GSE82246) and intestine (GEO: GSE83195). The reproducibility of RNA-seq data from wild-type and pax6b mutants was verified by a principal component analysis (PCA) obtained using the DESeq R package [68]. Differential expression (DE) analysis was performed using the R package DESeq2 [66], to identify all genes displaying significant change in the expression between wild-type and pax6b mutants (with fold discovery rate < 10%). The significance in the overlap of the sets of genes regulated by Pax6b in EECs and in PECs was determined with the R built-in statistical Fisher’s exact test.

Gene Ontology enrichment analysis

Gene Ontology (GO) enrichment analysis was performed on the different gene sets (genes selectively expressed either in EECs or in PECs, pax6b-regulated genes only in PECs, only in EECs and in both PECs and EECs) using the DAVID bioinformatics resources 6.8 [69] taking as background all the zebrafish genes. The enrichment analysis was focused on the GO biological process, molecular function and KEGG pathways with a statistical Fisher exact test p value < 0.05.

In situ hybridization and immunohistochemistry

Antisense RNA probe for the different genes was prepared as described by Thisse and Thisse [70]. Briefly, primers were designed to amplify a part of the transcript that is used as a template to synthesize the probe. The reverse primer at the 5′ contains the minimal promoter sequence for T3 RNA polymerase (5′-AATTAACCCTCACTAAAGGGAG-3′); templates were amplified by RT-PCR. Whole-mount in situ hybridization and fluorescent in situ hybridization (WISH and FISH, respectively) were performed as described previously [71] on wild-type (AB strain), pax6sa0086 null mutant or pax6bsunrise mutant embryos/larvae [13]. Immunohistochemistry (IHC) on whole-mount embryos/larvae was performed as described [72]. The antibodies used are the guinea pig polyclonal anti-PDX1 antibody (ab47308, Abcam) used at a dilution of 200×, the chicken anti-GFP (Aves lab) used at a dilution of 1000×, the rabbit anti-enkephalin (T4294; Peninsula Laboratories International, Inc.) used at a 400× dilution, the rabbit anti-PACAP (Adcyap1) (T-4473.0050; from Peninsula Laboratories International, Inc.) used at a 300× dilution, the monoclonal mouse Hu antibody (16A11: Invitrogen Cat. #A-21271) used at a 1000× dilution and Alexa Fluor secondary antibodies at 1000× dilution. Stained embryos were mounted in Prolong (Invitrogen) with DAPI and imaged using SP5 confocal microsope (Leica). For the immunostainings on adult tissue, the intestine of transgenic Tg(-8.5nkx2.2a:GFP)ia2 adult zebrafish was dissected and fixed for 1 day in PFA 4% rinsed in PBS and embedded in TEK using standard procedures and sectioned using cryostat for generating 10-μm sections which were mounted on glass slides. The primary and secondary antibodies were used as the same dilution as described above for whole-mount larvae.

Availability of data and materials

All data generated or analysed during this study are included in this published article, its supplementary information files and publicly available repositories. The raw RNA-seq datasets of EECs and PECs from wild-type and pax6b mutants have been deposited on Gene Expression Omnibus (GEO) under the accession number GSE149081 ( The EEC and PEC RNA-seq were compared to other published RNA-seq data from the heart (ArrayExpress: E-MTAB-460; GEO: GSE71755), brain (ArrayExpress: E-MTAB-460), liver (GEO: GSE82246) and intestine (GEO: GSE83195).



Enteroendocrine cell


Pancreatic endocrine cell


Enteric neuron


Whole-mount in situ hybridization


Fluorescent in situ hybridization


Count of reads per million


Fragments per kilobase of exon model per million reads mapped


Standard errors








Peptide YY




Glucose-dependent insulinotropic polypeptide


Vasoactive intestinal peptide




Insulin-like peptide 5


Neuromedin B a


Adenylate cyclase-activating polypeptide 1








  1. 1.

    Gribble FM, Reimann F. Enteroendocrine cells: chemosensors in the intestinal epithelium. Annu Rev Physiol. 2016;78:277–99.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Habib AM, Richards P, Cairns LS, Rogers GJ, Bannon CAM, Parker HE, et al. Overlap of endocrine hormone expression in the mouse intestine revealed by transcriptional profiling and flow cytometry. Endocrinology. 2012;153:3054–65.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Haber AL, Biton M, Rogel N, Herbst RH, Shekhar K, Smillie C, et al. A single-cell survey of the small intestinal epithelium. Nature. 2017;551:333–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Glass LL, Calero-Nieto FJ, Jawaid W, Larraufie P, Kay RG, Göttgens B, et al. Single-cell RNA-sequencing reveals a distinct population of proglucagon-expressing cells specific to the mouse upper small intestine. Mol Metab. 2017;6:1296–303.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Böttcher G, Sjöberg J, Ekman R, Håkanson R, Sundler F. Peptide YY in the mammalian pancreas: immunocytochemical localization and immunochemical characterization. Regul Pept. 1993;43:115–30.

    Article  PubMed  Google Scholar 

  6. 6.

    Wierup N, Svensson H, Mulder H, Sundler F. The ghrelin cell: a novel developmentally regulated islet cell in the human pancreas. Regul Pept. 2002;107:63–9.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Naya FJ, Huang HP, Qiu Y, Mutoh H, DeMayo FJ, Leiter AB, et al. Diabetes, defective pancreatic morphogenesis, and abnormal enteroendocrine differentiation in BETA2/neuroD-deficient mice. Genes Dev. 1997;11:2323–34.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    May CL, Kaestner KH. Gut endocrine cell development. Mol Cell Endocrinol. 2010;323:70–5.

    CAS  Article  Google Scholar 

  9. 9.

    Gierl MS, Karoulias N, Wende H, Strehle M, Birchmeier C. The zinc-finger factor Insm1 (IA-1) is essential for the development of pancreatic beta cells and intestinal endocrine cells. Genes Dev. 2006;20:2465–78.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Beucher A, Gjernes E, Collin C, Courtney M, Meunier A, Collombat P, et al. The homeodomain-containing transcription factors Arx and Pax4 control enteroendocrine subtype specification in mice. PLoS One. 2012;7:e36449.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Sander M, Neubuser A, Kalamaras J, Ee HC, Martin GR, German MS. Genetic analysis reveals that PAX6 is required for normal transcription of pancreatic hormone genes and islet development. Genes Dev. 1997;11:1662–73

    CAS  Article  Google Scholar 

  12. 12.

    Heller RS, Jenny M, Collombat P, Mansouri A, Tomasetto C, Madsen OD, et al. Genetic determinants of pancreatic epsilon-cell development. Dev Biol. 2005;286:217–24.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Verbruggen V, Ek O, Georlette D, Delporte F, Von Berg V, Detry N, et al. The Pax6b homeodomain is dispensable for pancreatic endocrine cell differentiation in zebrafish. J Biol Chem. 2010;285:13863–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Larsson LI, St-Onge L, Hougaard DM, Sosa-Pineda B, Gruss P. Pax 4 and 6 regulate gastrointestinal endocrine cell development. Mech Dev. 1998;79:153–9 Accessed 29 Jun 2017.

    CAS  Article  Google Scholar 

  15. 15.

    Hill ME, Asa SL, Drucker DJ. Essential requirement for Pax6 in control of enteroendocrine proglucagon gene transcription. Mol Endocrinol. 1999;13:1474–86

    CAS  Article  Google Scholar 

  16. 16.

    Hartenstein V, Takashima S, Hartenstein P, Asanad S, Asanad K. bHLH proneural genes as cell fate determinants of entero-endocrine cells, an evolutionarily conserved lineage sharing a common root with sensory neurons. Dev Biol. 2017;431:36–47.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Hartenstein V, Martinez P. Structure, development and evolution of the digestive system. Cell Tissue Res. 2019;377:289–92.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Annunziata R, Andrikou C, Perillo M, Cuomo C, Arnone MI. Development and evolution of gut structures: from molecules to function. Cell Tissue Res. 2019;377:445–58.

    Article  PubMed  Google Scholar 

  19. 19.

    Heller RS. The islets of Langerhans 2010;654:21–37. doi:

  20. 20.

    Youson JH, Al-Mahrouki AA. Ontogenetic and phylogenetic development of the endocrine pancreas (islet organ) in fishes. Gen Comp Endocrinol. 1999;116:303–35.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Arendt D. The evolution of cell types in animals: emerging principles from molecular studies. Nat Rev Genet. 2008;9:868–82.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Kin K, Nnamani MC, Lynch VJ, Michaelides E, Wagner GP. Cell-type phylogenetics and the origin of endometrial stromal cells. Cell Rep. 2015;10:1398–409.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Arntfield ME, van der Kooy D. β-Cell evolution: how the pancreas borrowed from the brain: the shared toolbox of genes expressed by neural and pancreatic endocrine cells may reflect their evolutionary relationship. BioEssays. 2011;33:582–7.

    Article  Google Scholar 

  24. 24.

    Tarifeño-Saldivia E, Lavergne A, Bernard A, Padamata K, Bergemann D, Voz MLML, et al. Transcriptome analysis of pancreatic cells across distant species highlights novel important regulator genes. BMC Biol. 2017;15:21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Bates JM, Mittge E, Kuhlman J, Baden KN, Cheesman SE, Guillemin K. Distinct signals from the microbiota promote different aspects of zebrafish gut differentiation 2006;297:374–386. Accessed 21 Aug 2019.

  26. 26.

    Chen YH, Lu YF, Ko TY, Tsai MY, Lin CY, Lin CC, et al. Zebrafish cdx1b regulates differentiation of various intestinal cell lineages. Dev Dyn. 2009;238:1021–32.

    CAS  Article  Google Scholar 

  27. 27.

    Flasse LCLC, Stern DGDG, Pirson JLJL, Manfroid I, Peers B, Voz MLMLML. The bHLH transcription factor Ascl1a is essential for the specification of the intestinal secretory cells and mediates Notch signaling in the zebrafish intestine. Dev Biol. 2013;376:187–97.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Roach G, Heath Wallace R, Cameron A, Emrah Ozel R, Hongay CF, Baral R, et al. Loss of ascl1a prevents secretory cell differentiation within the zebrafish intestinal epithelium resulting in a loss of distal intestinal motility. Dev Biol. 2013;376:171–86.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Wallace KN, Akhter S, Smith EM, Lorent K, Pack M. Intestinal growth and differentiation in zebrafish. Mech Dev. 2005;122:157–73.

    CAS  Article  Google Scholar 

  30. 30.

    Roy-Carson S, Natukunda K, Chou H, Pal N, Farris C, Schneider SQ, et al. Defining the transcriptomic landscape of the developing enteric nervous system and its cellular environment. BMC Genomics. 2017;18:290.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Liu Y, Li S, Huang X, Lu D, Liu X, Ko WH, et al. Identification and characterization of a motilin-like peptide and its receptor in teleost. Gen Comp Endocrinol. 2013;186:85–93.

    CAS  Article  Google Scholar 

  32. 32.

    Delporte FM, Pasque V, Devos N, Manfroid I, Voz ML, Motte P, et al. Expression of zebrafish pax6b in pancreas is regulated by two enhancers containing highly conserved cis-elements bound by PDX1, PBX and PREP factors. BMC Dev Biol. 2008;8:53. PMID: 18485195; PMCID: PMC2409314.

  33. 33.

    Uyttebroek L, Shepherd IT, Harrisson F, Hubens G, Blust R, Timmermans JP, et al. Neurochemical coding of enteric neurons in adult and embryonic zebrafish (Danio rerio). J Comp Neurol. 2010;518:4419–38.

    CAS  Article  Google Scholar 

  34. 34.

    van Arensbergen J, García-Hurtado J, Moran I, Maestro MA, Xu X, Van de Casteele M, et al. Derepression of Polycomb targets during pancreatic organogenesis allows insulin-producing beta-cells to adopt a neural gene activity program. Genome Res. 2010;20:722–32.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Collombat P, Mansouri A, Hecksher-Sorensen J, Serup P, Krull J, Gradwohl G, et al. Opposing actions of Arx and Pax4 in endocrine pancreas development. Genes Dev. 2003;17:2591–603.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Djiotsa J, Verbruggen V, Giacomotto J, Ishibashi M, Manning E, Rinkwitz S, et al. Pax4 is not essential for beta-cell differentiation in zebrafish embryos but modulates alpha-cell generation by repressing arx gene expression. BMC Dev Biol. 2012;12:37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Ohta Y, Kosaka Y, Kishimoto N, Wang J, Smith SB, Honig G, et al. Convergence of the insulin and serotonin programs in the pancreatic β-cell. Diabetes. 2011;60:3208–16.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Helker CSM, Mullapudi S-T, Mueller LM, Preussner J, Tunaru S, Skog O, et al. A whole organism small molecule screen identifies novel regulators of pancreatic endocrine development. Development. 2019;146:dev172569.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Yee NS, Lorent K, Pack M. Exocrine pancreas development in zebrafish. Dev Biol. 2005;284:84–101

    CAS  Article  Google Scholar 

  40. 40.

    Spanjaard B, Hu B, Mitic N, Olivares-Chauvet P, Janjuha S, Ninov N, et al. Simultaneous lineage tracing and cell-type identification using CRISPR–Cas9-induced genetic scars. Nat Biotechnol. 2018;36:469–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Kleinjan DA, Bancewicz RM, Gautier P, Dahm R, Schonthaler HB, Damante G, et al. Subfunctionalization of duplicated zebrafish pax6 genes by cis-regulatory divergence. PLoS Genet. 2008;4:e29.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Cardoso JC, Vieira FA, Gomes AS, Power DM. The serendipitous origin of chordate secretin peptide family members. BMC Evol Biol. 2010;10:135.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Hort Y, Baker E, Sutherland GR, Shine J, Herzog H. Gene duplication of the human peptide YY gene (PYY) generated the pancreatic polypeptide gene (PPY) on chromosome 17q21.1. Genomics. 1995;26:77–83.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Sundström G, Larsson TA, Brenner S, Venkatesh B, Larhammar D. Evolution of the neuropeptide Y family: new genes by chromosome duplications in early vertebrates and in teleost fishes. Gen Comp Endocrinol. 2008;155:705–16.

    Article  PubMed  Google Scholar 

  45. 45.

    Biemar F, Argenton F, Schmidtke R, Epperlein S, Peers B, Driever W. Pancreas development in zebrafish: early dispersed appearance of endocrine hormone expressing cells and their convergence to form the definitive islet. Dev Biol. 2001;230:189–203.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Tissue-based map of the human proteome. Science. 2015;347:1260419.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Mutch DM, Anderle P, Fiaux M, Mansourian R, Vidal K, Wahli W, Williamson G, Roberts MA. Regional variations in ABC transporter expression along the mouse intestinal tract. Physiol Genomics. 2004;17(1):11–20. PMID: 14679303.

  48. 48.

    Segerstolpe Å, Palasantza A, Eliasson P, Andersson E-M, Andréasson A-C, Sun X, et al. Single-cell transcriptome profiling of human pancreatic islets in health and type 2 diabetes. Cell Metab. 2016;24:593–607.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Crosnier C, Vargesson N, Gschmeissner S, Ariza-McNaughton L, Morrison A, Lewis J. Delta-Notch signalling controls commitment to a secretory fate in the zebrafish intestine. Development. 2005;132:1093–104.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Jensen J, Pedersen EE, Galante P, Hald J, Heller RS, Ishibashi M, et al. Control of endodermal endocrine development by Hes-1. Nat Genet. 2000;24:36–44.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Falkmer S, Dafgård E, el-Salhy M, Engström W, Grimelius L, Zetterberg A. Phylogenetical aspects on islet hormone families: a minireview with particular reference to insulin as a growth factor and to the phylogeny of PYY and NPY immunoreactive cells and nerves in the endocrine and exocrine pancreas. Peptides. 1985;6(Suppl 3):315–20.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Youson JH. The agnathan enteropancreatic endocrine system: phylogenetic and ontogenetic histories, structure, and function1. Am Zool. 2000;40:179–99.

    Google Scholar 

  53. 53.

    Perillo M, Arnone MI. Characterization of insulin-like peptides (ILPs) in the sea urchin Strongylocentrotus purpuratus: insights on the evolution of the insulin family. Gen Comp Endocrinol. 2014;205:68–79.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Ariyachet C, Tovaglieri A, Xiang G, Lu J, Shah MS, Richmond CA, et al. Reprogrammed stomach tissue as a renewable source of functional β cells for blood glucose regulation. Cell Stem Cell. 2016;18:410–21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Chen Y-J, Finkbeiner SR, Weinblatt D, Emmett MJ, Tameire F, Yousefi M, et al. De novo formation of insulin-producing “neo-β cell islets” from intestinal crypts. Cell Rep. 2014;6:1046–58.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Talchai C, Xuan S, Kitamura T, DePinho RA, Accili D. Generation of functional insulin-producing cells in the gut by Foxo1 ablation. Nat Genet. 2012;44:406–12.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Thurner M, van de Bunt M, Torres JM, Mahajan A, Nylander V, Bennett AJ, et al. Integration of human pancreatic islet genomic data refines regulatory mechanisms at type 2 diabetes susceptibility loci. Elife. 2018;7.

  58. 58.

    Mahajan A, Taliun D, Thurner M, Robertson NR, Torres JM, Rayner NW, et al. Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps. Nat Genet. 2018;50:1505–13.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Scott Heller R, Stoffers DA, Liu A, Schedl A, Crenshaw EB, Madsen OD, et al. The role of Brn4/Pou3f4 and Pax6 in forming the pancreatic glucagon cell identity. Dev Biol. 2004;268:123–34.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF. Stages of embryonic development of the zebrafish. Dev Dyn. 1995;203:253–310

    CAS  Article  Google Scholar 

  61. 61.

    Delporte FM, Pasque V, Devos N, Manfroid I, Voz ML, Motte P, Biemar F, Martial JA, Peers B. Expression of zebrafish pax6b in pancreas is regulated by two enhancers containing highly conserved cis-elements bound by PDX1, PBX and PREP factors. BMC Dev Biol. 2008 May 16;8:53. PMID: 18485195; PMCID: PMC24093.

  62. 62.

    Pauls S, Zecchin E, Tiso N, Bortolussi M, Argenton F. Function and regulation of zebrafish nkx2.2a during development of pancreatic islet and ducts. Dev Biol. 2007;304:875–90.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Picelli S, Faridani OR, Björklund AK, Winberg G, Sagasser S, Sandberg R. Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc. 2014;9:171–81.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. Epub 2014 Sep 25. PMID: 25260700; PMCID: PMC4287950.

  66. 66.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Ghaye AP, Bergemann D, Tarifeño-Saldivia E, Flasse LC, Von Berg V, Peers B, et al. Progenitor potential of nkx6.1-expressing cells throughout zebrafish life and during beta cell regeneration. BMC Biol. 2015;13:70.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.

    CAS  Article  Google Scholar 

  69. 69.

    Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    CAS  Article  Google Scholar 

  70. 70.

    Thisse C, Thisse B. High-resolution in situ hybridization to whole-mount zebrafish embryos. Nat Protoc. 2008;3:59–69.

    CAS  Article  Google Scholar 

  71. 71.

    Mavropoulos A, Devos N, Biemar F, Zecchin E, Argenton F, Edlund H, et al. sox4b is a key player of pancreatic alpha cell differentiation in zebrafish. Dev Biol. 2005;285:211–23.

    CAS  Article  PubMed  Google Scholar 

  72. 72.

    Flasse LC, Pirson JL, Stern DG, Von Berg V, Manfroid I, Peers B, et al. Ascl1b and Neurod1, instead of Neurog3, control pancreatic endocrine cell fate in zebrafish. BMC Biol. 2013;11:78.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


We are grateful for Christian Helker and Didier Stainier for the zebrafish transgenic line TgBAC(pdx1:EGFP)bns13, M. Takamiya and Uhle Strähle for the sunrise mutant line, Francesco Argenton for the Tg(nkx2.2::GFP) and Kenneth Wallace for the phox2b probe. We thank the following technical platforms: GIGA-Zebrafish (H Pendeville), GIGA-Cell Imaging and Flow Cytometry (S Ormenese and S Raafat), GIGA-histology (C Humblet) and the GIGA-Genotranscriptomic platforms (W Coppieters, and L Karim).


AL, JP, LF and A-SR were supported by FRIA (Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture) and the Léon Fredericq fund. ETS was supported by WBI (Wallonie-Bruxelles International), Becas Chile and the Léon Fredericq fund. This work was funded by the FNRS-FRS, the ZENCODE-ITN european project 643062, and the Fonds Speciaux from the ULiège (University of Liège). BP, IM and MLV are associate researchers from FRS/FNRS (Fonds National pour la Recherche Scientifique).

Author information




ETS and JP carried out the RNA-seq experiments. AL and ETS performed the RNA-seq analyses. AL, ASR, LF and MLV performed the WISH and FISH experiments. AL carried out the bioinformatics analyses concerning the transcriptome comparison, differential gene expression and cell subtype markers. BP performed the IF. BP and MLV conceived the study. BP, AL, MLV and IM participated in the interpretation of the data and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Marianne L. Voz or Bernard Peers.

Ethics declarations

Ethics approval and consent to participate

All animal experiments were conducted according to the national guidelines and were approved by the ethical committee of the University of Liège (protocol numbers 1328 and 1557). Consent to participate is not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Figure S1.

The pax6b:GFP transgene is not expressed in enteric neurones. Immunofluorescence on 5 dpf Tg(pax6b:GFP) larva using antibodies against GFP (green) and against the enteric neurone marker Hu (red). Confocal vues of the bulb intestine (upper panels), mid-intestine (middle panels) and posterior intestine (lower panels) showing no colocalisation of GFP with Hu. Dapi staining is in blue. Scale bar =50μm.

Additional file 2: Table S1.

Expression level of all genes in zebrafish enteroendocrine cells. The expression levels is given for all genes in CPM (counts per million) for the four EEC samples (column C to F) as well as the mean and the standard deviation obtained for each gene (columns G and H). The last column gives the mean expression levels in FPMK.

Additional file 3: Figure S2.

Quantification of EECs expressing different neuropeptide transcripts in 4 dpf zebrafish larvae. The number of EECs expressing each hormones was determined by counting the labelled cells after WISH using the corresponding hormone probes on 4 dpf larvae. Each point in the graph represents the number of labelled cells in one larva. Bars represent the mean values and S.E.

Additional file 4: Figure S3.

Expression of vipb and adcyap1a transcripts in zebrafish EECs. Confocal images of the zebrafish gut from larvae stained by double fluorescent in situ hybridization (FISH) using the enteric neurone marker phox2b (green) and the vipb or adcyap1a probes (red). (A) general vue of the gut showing co-localisation of vipb and phox2b in ENs at the level of anterior intestine (upper left part of the image) and vipb+ EECs in the posterior intestine (bottom right of image). (B) higher magnification of the gut showing three vipb+ phox2b+ ENs and two vipb+ EECs. (C) adcyap1a+ cells are distinct from phox2b+ ENs which are located outside the intestinal epithelium (Dapi staining in blue).

Additional file 5: Table S2.

List of transcription factors expressed in both EECs and PECs, specific for EECs and specific for PECs. The expression level is given in FPKM. The TF expressed in both PECs and EECs were selected based on their expression level above the threshold of 1 FPKM in both cell types. The TF expressed specifically in PECs were obtained by selecting those expressed above 1 FPKM only in PECs and with an expression ratio of PEC/EEC above 5-fold. Inversely, EEC-specific TF were obtained by selecting those expressed above 1 FPKM only in EECs and with an expression ration EEC/PEC above 5-fold.

Additional file 6: Table S3.

Lists of genes selectively expressed in EECs and in PECs with their gene ontology enrichment analyses. The genes specifically expressed in PECs were selected by their expression above 5 CPM in PECs and below 1 CPM in EECs (sheet1). Inversely EEC-specific genes were selected by their expression above 5 CPM in EECs and below 1 in PECs (sheet 2). The GO terms obtained for the PEC-specific genes and the EEC-specific genes are given in sheets 3 and 4 respectively.

Additional file 7: Table S4.

Expression level of all genes EECs and PECs (wild type and pax6bsa0086). The expression level of all genes is given for wild-type and pax6b mutant EECs and PECs in normalized CPM. The expression levels in each sample are given in the excel sheet 1(“CPM samples”) and the means and standard deviations are given in excel sheet 2 (“mean and Sdt dev.”).

Additional file 8: Table S5.

List of Pax6b-regulated genes and GO enrichment analysis. The excel sheets 1, 3 and 5 give respectively the lists of genes regulated by Pax6b in both PECs and EECs (sheet 1), in PECs only (sheet 3) and in EECs only (sheet 5) together with the expression level in normalized CPM in wild type and pax6b mutant EECs and PECs, the expression fold changes and the P-adjusted values. The excel sheets 2, 4 and 6 display respectively the enrichment of GO terms and pathways for the Pax6b-regulated genes in PECs and EECs (sheet 2), in PECs only (sheet 4) and in EECs only (sheet 6) together with the P-values for the Fischer exact test and the lists of genes associated to each GO term.

Additional file 9: Table S6.

list of markers for each pancreatic endocrine cell subtypes. The list of markers for the pancreatic alpha-, beta-, delta- and epsilon- cell subtypes was obtained from previous published RNA-seq studies [24, 40]. This list was used to determine the proportion of cell subtype markers regulated by Pax6b as shown in the upper panel of Fig. 7.

Additional file 10: Figure S4.

Sequence similarity between pancreatic polypeptide and of peptide YY. Alignment of the amino-acid sequence of pancreatic polypeptide (PP) from mice and rat with peptide YY from mice and zebrafish (DrPyy A or b).

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lavergne, A., Tarifeño-Saldivia, E., Pirson, J. et al. Pancreatic and intestinal endocrine cells in zebrafish share common transcriptomic signatures and regulatory programmes. BMC Biol 18, 109 (2020).

Download citation


  • Pancreatic endocrine cells
  • Enteroendocrine
  • Transcriptome
  • RNA-seq
  • Pax6
  • Zebrafish
  • Hormone
  • Transcriptome comparison