- Research article
- Open Access
A pancreatic exocrine-like cell regulatory circuit operating in the upper stomach of the sea urchin Strongylocentrotus purpuratus larva
BMC Evolutionary Biologyvolume 16, Article number: 117 (2016)
Digestive cells are present in all metazoans and provide the energy necessary for the whole organism. Pancreatic exocrine cells are a unique vertebrate cell type involved in extracellular digestion of a wide range of nutrients. Although the organization and regulation of this cell type is intensively studied in vertebrates, its evolutionary history is still unknown. In order to understand which are the elements that define the pancreatic exocrine phenotype, we have analyzed the expression of genes that contribute to specification and function of this cell-type in an early branching deuterostome, the sea urchin Strongylocentrotus purpuratus.
We defined the spatial and temporal expression of sea urchin orthologs of pancreatic exocrine genes and described a unique population of cells clustered in the upper stomach of the sea urchin embryo where exocrine markers are co-expressed. We used a combination of perturbation analysis, drug and feeding experiments and found that in these cells of the sea urchin embryo gene expression and gene regulatory interactions resemble that of bona fide pancreatic exocrine cells. We show that the sea urchin Ptf1a, a key transcriptional activator of digestive enzymes in pancreatic exocrine cells, can substitute for its vertebrate ortholog in activating downstream genes.
Collectively, our study is the first to show with molecular tools that defining features of a vertebrate cell-type, the pancreatic exocrine cell, are shared by a non-vertebrate deuterostome. Our results indicate that the functional cell-type unit of the vertebrate pancreas may evolutionarily predate the emergence of the pancreas as a discrete organ. From an evolutionary perspective, these results encourage to further explore the homologs of other vertebrate cell-types in traditional or newly emerging deuterostome systems.
The ability of a species to survive in a new environment is strictly related to its capacity to successfully digest and assimilate different food sources. A characteristic that is shared by all the living forms is the ability to digest large molecules available in their habitat to produce energy. Digestion relies on the function of enzymes that are able to break down proteins, lipids and sugars in the diet. The number of specialized enzyme-producing cell types involved in food digestion increases with the complexity of the animal body and the variety of the meal. Metazoans typically have specialized exocrine cell types for the secretion of digestive enzymes. For example, gland cells that produce and release zymogens in the gastric cavity are scattered along the hydra gastroepidermis [1, 2] and along the digestive filaments of corals . Cells rich in granular inclusions involved in extracellular digestion are localized along the foregut and in the gut caeca of different species of flat worms [4–6]. In different regions of the Drosophila midgut there are enterocytes that secrete a vast array of digestive enzymes, with expression levels matching organismal requirements . A specialized organ geared toward the production of digestive enzymes and hormones, the pancreas, appeared uniquely in the vertebrate lineage. The pancreas is a complex organ evolved mainly for food digestion (exocrine cells) and maintaining blood sugar levels (endocrine cells). The exocrine pancreas constitutes the majority of the mass of the pancreas and it is composed by exocrine cells grouped into acini that secrete enzymes into the intestine . The pancreas originates early in development from progenitor cells located within the endodermal epithelium. Two members of the basic helix-loop-helix (bHLH) family of protein are known to be critical for the exocrine pancreas differentiation and identity, Ptf1a (pancreatic transcriptional factor 1) and Mist1 [9, 10]. Ptf1a has a double role, first in early pancreas specification, and later in acinar cell differentiation and maintenance [11–13], while Mist1 is necessary for full maturation of the acinar cell phenotype . In the mouse, it has been found that in the absence of Hnf1β, Ptf1a expression is not induced, leading to defective specification and reduction of the ventral and dorsal pancreas . Furthermore, several signaling factors are involved in balancing exocrine cell proliferation and differentiation. Among all, Notch appears to prevent pancreatic exocrine development, while FGF signaling mediates growth, morphogenesis and differentiation of exocrine cells [16–18]. In both mouse and zebrafish, Notch and its target genes inhibit the activity of the Ptf1 complex in the exocrine cells, and loss of Notch function results in accelerated development of exocrine pancreas .
One of the most abundant and evolutionary conserved miRNA involved in pancreas development is miR-375. MiR-375 negatively regulates glucose-stimulated insulin exocytosis by targeting myothrophin mRNA, a protein involved in insulin secrection in pancreatic β-cells . A loss of function approach in zebrafish revealed that miR-375 is essential for the formation of insulin producing β-cells because its knockdown mainly results in malformation of the endocrine pancreas . In addition, studies on mice lacking miR-375 demonstrated that it controls β- and α-cells mass by regulating genes involved in cellular growth and differentiation .
Echinoderms—a group of animals that comprise sea urchins, sea stars, sea lilies, brittle stars and sea cucumbers- belong to the deuterostome clade. The relatively simple development of sea urchin embryos, coupled with the availability of molecular tools for gene perturbation, enable comparative studies on cell specification and developmental mechanisms with other deuterostomes . In particular, echinoderm systems present, in the comparison with vertebrate deuterostomes, the great advantage of a simplified gene toolkit, since this group split before the duplication events occurred at the emergence of Craniata. Importantly for this study, Strongylocentrotus purpuratus has an indirect development that generates a bilateral free-swimming larva whose behavior is centered on feeding.
The sea urchin larval gut is a simple tripartite tract composed of a muscular esophagus that exhibits strong contractions, a large spherical stomach with sphincters at both openings, and a small tubular intestine, which exits through the anus. In the stomach of sea urchin larvae, there are cells with the morphological features of zymogen secreting cells that are able to synthesize digestive enzymes, such as β-glucanase and α-amylase [24–27]. Gastric exocrine cells have been described also in the stomach of the adult form of S. purpuratus. These cells are highly concentrated in the upper stomach and electron microscopy analyses reveal that they are strikingly similar to mammalian zymogen cells, such as cells of the gastric mucosa and pancreatic acinar cells .
Although pancreatic cell types have been extensively characterized in the main vertebrate model systems, most of the studies in non-vertebrate animals have examined gut cells at the morphological level only and molecular information is still very scarce. Therefore, a few data are currently available on the genes/genetic characterization of the pancreatic exocrine-like cell types in non-vertebrate metazoans. In this study, we use a candidate gene approach to identify and characterize the pancreatic exocrine cell type in the sea urchin embryo and larva. We defined a unique population of cells clustered in the sea urchin larva upper stomach that responds to food and expresses the homologs of the pancreatic genes Ptf1a and Mist1 together with at least three digestive enzymes. Our findings are a first step to understand the evolution of pancreatic exocrine cells. We report that the components of the specification pathway and expression for zymogen genes are active in an early branching deuterostome, which suggests that the exocrine cell type may predate the evolution of the vertebrate pancreas.
Animal husbandry, embryo and larva cultures, feeding experiments
Adult S. purpuratus were obtained from Patrick Leahy (Kerchoff Marine Laboratory, California Institute of Technology, Pasadena, CA, USA) and housed in circulating seawater aquaria at the Stazione Zoologica Anton Dohrn of Naples. Adult sea urchin maintenance, embryo and larva cultures and feeding experiments were performed as previously described .
RNA whole mount in situ hybridization
For fluorescent whole mount in situ hybridization (FISH), we followed the protocol outlined in Cole et al. 2009 with the modification described in . Signal was developed with fluorophore-conjugated tyramide (1:400 reagent diluents, Perkin Elmer) following instructions. For all the genes, labeled probes were transcribed from linearized DNA as described in . SpmiR-375 probe has been synthesized and DIG labelled from Exiqon, and the sequence is: 5’/DigN/TGACGCGAGCCGAACGAACAAA/3’DigN/. The double FISH procedure for SpCpa2L and SpmiR-375 was performed as cited with the only exception that the miRNA probe concentration was 0025 pmol/μl and the samples were hybridized 5 days at 42 °C. mir375 probe has already been used by Christodoulou et al. 2010. Primers used to amplify the riboprobes are summarized in Additional file 1. Templates of all the probes were sequenced prior to probe generation and cloned in the pGEM®-T Easy Vector (Promega, Madison, WI, USA). Sense probes were synthesized to test the specificity of the antisense probe signal. FISH was imaged with a Zeiss 510 Meta confocal microscope.
The γ-secretase inhibitor DAPT (N-[N-(3,5-difluorophenacetyl)-L-alanyl]-S-phenylglycine t-butyl ester, Sigma–Aldrich, St. Louis, MO)  was dissolved in DMSO and added at 17 h to a final concentration of 8 μM, to avoid the early and toxic effects of the drug, as described in Materna et al. 2012. A corresponding volume of DMSO was added as control.
Perturbation experiments with MO injection
For each experiment and for each morpholino oligonucleotide (MO), 200 to 400 eggs were injected with approximately 2-4pl of oligonucleotide injection solution and each experiment was repeated three times. As a negative control, fertilized eggs were injected with the standard control morpholino (GeneTools, Philomath, OR) and compared side-by-side with knockdown embryos. The injection of the standard control morpholino did not have any effect on the development of embryos. In a pilot experiment, different concentrations of the MASOs were injected and the morphant’s phenotype was observed in order to test for the right concentration and avoid secondary effects or MASO toxicity. The highest concentration that was not toxic for the embryos was used in each experiment. The injected embryos were cultured at 15 °C until 50 h and 67 h of development, when they were fixed and used to analyze the expression of target genes in FISH or qPCR experiments. A translation-blocking antisense MO against SpHnf1a was used at a final concentration in the injection solution of 200 μM (5’-CTAGTTCGTCACCCGAATGCAGCAT-3’, first tested by Peterson and Davidson 2011). A translation-blocking antisense MO against SpPtf1a was newly designed and used at a final concentration of 300 μM (SpPtf1a_MO1: 5’-ATATTTTCCATAGTGATCTCTGAGT-3’). A second translation-blocking antisense MO against SpPtf1a was injected at a final concentration of 150 μM and used to confirm the specific effect of SpPtf1a MO_1 (SpPtf1a_MO2: 5’-GCGCAGGTGGATTATCAAATTGTTC-3’).
RNA extraction and quantitative real-time PCR (qPCR)
For qPCR analysis of temporal expression profiles, total RNA was isolated from cultures derived from three mixed batches and extracted with Eurozol (EuroClone, Celbio, Milan, Italy). For SpPtf1a morphants and relative controls and for feeding experiments, RNA was extracted using the RNAquous kit (Ambion). For all the experiments, samples were treated with DNase I (Ambion) to remove DNA contamination as described by the manufacturer. First-strand cDNA was synthesized from total RNA using the SuperScript VILO™ cDNA Synthesis Kit (Invitrogen) following the manufacturer’s instructions. Reactions were performed as described in  using the ViiA7 REAL TIME PCR with SYBR Green chemistry (Applied Biosystems, Foster City, CA). For all qPCR experiments, ubiquitin mRNA was used as a reference because remains relatively constant during development . For knockdown and feeding experiments, three biological replicates were analyzed, ddCt values were calculated between experiment and control embryos and converted to fold differences using ubiquitin as a reference (FD = 2 ΔΔCt). A threshold of 2-fold difference was chosen as a significant change. Primer sets were chosen to amplify products 100 to 200 bp in length and sequences are reported in Additional file 1. qPCR experiments in different feeding conditions and the relative statistics were performed as previously described in Perillo and Arnone 2014.
The coding sequences of E47, Rbpl, rnPtf1a and spPtf1a were cloned into the pCDNA3.1 backbone. Ptf1a Luciferase assay was conducted as previous described . Briefly, 293 t or Hela cells were plated at a density of 1x105 cells/well in a 24-well plate and transfected with 512.5 ng total DNA/well using Fugene 6 transfection reagent (Promega). For individual constructs, the following concentrations were used (per well): 12.5 ng renilla (pRL-TK, Promega), 125 ng 4x ptf1a:luciferase reporter (contains 4x Ptf1a binding motif cloned from rat Chymotrypsinogen B promoter, gift of Dr. Masashi Kawaichi), 125 ng human E47, 125 ng Rbpl and 125 ng of either rat rnPtf1a or sea urchin spPtf1a. Different combinations of plasmids were cotransfected and pcDNA3.1 (Invitrogen) was used as filler DNA to make the total amount of DNA per well uniform across different conditions. Two days post-transfection, cells were harvested and analyzed for firefly and renilla luciferase activities using the Dual-luciferase reporter assay system (Promega) per manufactor’s protocol.
To identify putative Ptf1a binding sites on sea urchin genes of interest, we used the list of Ptf1a binding sites of the mouse acinar digestive enzyme gene promoters published by Masui et. al. (2007). From this list we built base frequency tables for the E-box and TC-box separately, i.e. counts the number of times each base occurs at each position. The base frequency table has four rows (one row for each letter: A, C, G and T) and the number of columns is equal to the motif length. Then, we calculated the weight matrices for the E-box and TC-box, the matrices have the same number of rows and columns as the corresponding frequency tables, with the value at each position being the natural logarithm of the value from the frequency table divided by the number of sequences, i.e. the weight matrices contain the estimates of the log-probabilities of each base occurring at each position in true binding sites, based on the sample of known sites. Since the gap between the two boxes is about one or two helical DNA-turns, we introduced a penalty function based on a cosine function, whose period and phase shift are 10.5 and 21 base pairs separately. Thus for a given sequence, we can calculate a matching score based on the sum of the weights that correspond to the sequence, which should be equal to the log-probability of seeing that sequence given that it is a binding site, and also the gap penalty function for the distance between the boxes.
Characterization of regulatory and differentiation signatures for sea urchin pancreatic exocrine-like cells
We identified sea urchin orthologs of regulatory genes with a known conserved role in pancreatic exocrine cell development. The embryonic temporal expression of the transcriptional factors (TFs) analyzed in this study agrees with data already published from other authors [36, 37], but our analysis includes also late larval stages.
The hepatocyte nuclear factor (HNF) gene family is part of the HNF homeobox class and includes HNF1α and HNF1β, two paralogs with known interchangeable functions in several contexts [38, 39]. The sea urchin genome has only one homolog named SpHnf1. Its expression begins in the presumptive endoderm at blastula stage and by mid-gastrula appears confined to the gut . To clarify its later functions in development, we defined the spatial expression of SpHnf1 in post-gastrula stage embryos. During gastrulation, SpHnf1 is expressed throughout the hindgut and the midgut, but particularly enriched in a group of cells just below the cardiac sphincters (Fig. 1a). This expression fades away in the early pluteus (Fig. 1b) and it is not detectable in older stages, consistently with the low quantitative level of expression exhibited by this gene at larval stages (Fig. 1g).
The sea urchin ortholog of Ptf1a  has two distinct domains of expression. At gastrula stage (48hpf, Fig. 1c), transcripts start to accumulate in cells situated right-lateral toward the dorsal side of the midgut, next to the cardiac sphincter, and some hours later it equally expands on both sides. The second domain of expression is in the ectoderm, particularly at the apical organ and in scattered cells of the ciliary band. Throughout the larval stages, the stomach cells increase in number, and we also observe an increase in SpPtf1a positive cells (Fig. 1d). The expression of another bHLH gene, SpMist1 , starts at late gastrula stage in the upper stomach and in scattered cells of the apical organ and ciliary band (Fig. 1e, f). SpMist1 gene expression increases during gastrulation (2 days), although at early pluteus stage (3 days) SpPtf1a decreases while SpMist1 reaches its maximum expression (Fig. 1g). Similarly, the expression of both genes decreases at 7 days after fertilization and increases again after a few days.
We next identified the genes that most likely encode for enzymes used in food digestion. We compared the predicted ORFs of digestive enzyme gene families with the ESTs of an annotated sea urchin late larva cDNA library (MPMGp691, http://owww.molgen.mpg.de/ag_seaurchin/, ). The carboxypeptidase family was the first to be analyzed. Among 7 sea urchin orthologs, only SpCpa2L showed a significant representation in the above mentioned larval cDNA library. As shown in Fig. 1h , SpCpa2L transcripts are expressed at gastrula stage in two groups of cells in the upper midgut, next to the cardiac sphincter. In a temporal progression, SpCpa2L expression begins at the upper lateral stomach cells and later it extends ventrally from both sides until coalescing at the ventral most cells, while leaving a gap at the dorsal most part of the stomach (Fig. 1i and k). Based on the EST abundance at larval stages, SpPnlp2 and SpPnlp5 were selected among the lipases and blast analysis revealed that they encode the same protein, hereby named SpPnlp2/5. Likewise, SpAmy3 was selected among the amylases. For both genes, transcripts are detectable one day after SpCpa2L starts to be expressed, and they are produced in the same above described cells of the upper stomach (Fig. 1j-m). Consistently, the levels of expression of these three genes steadily increase during development and it reaches a peak at 10 days (Fig. 1n).
Pancreatic exocrine-like cells are clustered in the sea urchin larva upper stomach
Our examination of the expression of individual pancreatic regulatory and terminal differentiation genes suggests that they are expressed within a common exocrine-like cell population. To further test this premise, we tested co-expression of these factors by double fluorescent in situ. The endodermal expression of SpPtf1a is initially confined to a group of ~3 cells in the midgut within the broad SpHnf1 domain (Fig. 2a). We next asked if other homologs of pancreatic exocrine genes were co-expressed in those cells. The first digestive enzyme that appears in the SpPtf1a and SpHnf1 positive domain is SpCpa2L (Fig. 2b). At gastrula stage, a group of cells in the dorsal-right midgut co-express SpCpa2L and SpPtf1a. Some hours later, at the pluteus stage, cells expressing these two markers appear ventrally, forming a belt-like structure (Fig. 2c; see also Fig. 2I in  for co-expression of these two genes at the pluteus stage). Remarkably, only a few cells of the dorsal stomach do not express these markers. Considering the high expression of SpCpa2L transcripts, all confined to the Ptf1a positive cells (SpPtf1a+), we used this gene expression as a marker for the pancreatic exocrine-like cells. Thus, we analyzed SpMist, SpPnlp2/5 and SpAmy3 spatial expressions relative to SpCpa2L at larval stages. As shown in Fig. 2c-e, all the orthologs of exocrine genes that we considered in this study are co-expressed in the pancreatic exocrine-like cells. Remarkably, the same cells express the evolutionary conserved pancreatic microRNA miRNA-375. Figure 2f shows that miRNA-375 is strongly expressed throughout the ciliary band and in the stomach, where it is localized exclusively in the pancreatic exocrine-like cells. In addition, this gene is expressed in the coelomic pouches and the probe signal is stronger in the right one.
Notch signaling inhibits differentiation towards pancreatic exocrine cell fate
Given its known role in pancreatic differentiation, we next tested the role of the Notch signaling pathway in sea urchin pancreatic exocrine-like cell differentiation, using the Notch inhibitor DAPT. The first striking effect of this treatment is the expansion of the SpPtf1a + cells in the gastrula and larva upper midgut/stomach (Fig. 3). Notably, the SpPtf1a expansion in the treated embryos never overcomes the cardiac sphincter border. Despite the strong effect on the endodermal cells, ectodermal SpPtf1a expression is not affected. Moreover, 35 % of the analyzed plutei showed an ectopic expression of SpPtf1a in a small group of cells in the anus (Fig. 3d, red arrow). Conversely, the expression of the early endodermal gene SpHnf1 is not affected by the DAPT treatment at gastrula stage (Fig. 3a, b), in agreement with the previously reported data at blastula stage . When Notch signaling is perturbed, the expression of SpMist1 exhibits opposite effects than that of SpPtf1a. At the pluteus stage, the gene expression in the stomach is unaltered, while in the esophagus and in the apical organ SpMist1 is turned on also in cells that are adjacent to the wild type SpMist1+ cells(Fig. 3e and f and arrows in Fig. 3f).
Many modulators of Notch signaling are known, among which there is the O-fucose β-1,3 N-acetylglucosaminyltransferase Fringe, whose role is to potentiate Notch activation by Delta . Since a Notch signaling requiring Fringe is involved in specification of endoderm and boundary formation during sea urchin embryonic development , we asked if Notch inhibition affects Fringe expression also in the pluteus larva. The sea urchin ortholog of Fringe, SpFng, is widely expressed in the pluteus mesoderm, ectoderm and endoderm (Fig. 3g). Our analysis demonstrated that perturbation of Notch signaling causes a consistent reduction of SpFng expression in all the tissues where it is normally expressed (Fig. 3h), even in the pluteus.
Transcriptional control of sea urchin pancreatic exocrine-like cell development
We next asked if the TFs expressed early in development are involved in a gene regulatory module that specifies pancreatic exocrine cells. To this aim, specific morpholino oligonucleotides were injected in the sea urchin fertilized eggs to block either SpHnf1 or SpPtf1a translation. SpHnf1 MO sequence specificity has already been demonstrated by Peterson and Davidson (2001), and we used that same morpholino in our experiments. In SpHnf1 morphants, SpPtf1a expression in the midgut is downregulated, while its expression in the ectoderm remains unaltered (Fig. 4a, b). Fluorescent in situ hybridization (FISH) or quantitative PCR (qPCR) analyses were performed on larvae to test the effects of SpPtf1a MO perturbation on putative target genes. We observed that the spatial expression of SpCpa2L, SpPnlp2/5 and SpAmy3 is strongly reduced (Fig. 4e-h). For SpCpa2L, which among the enzymes shows the highest transcript expression at 3 days, changes at the transcript levels were also analyzed by qPCR analysis, confirming that its expression is strongly downregulated (Fig. 4k). Conversely, in SpPtf1a morphants additional groups of SpPtf1a + cells appear in the upper stomach, both ventrally and dorsally (Fig. 4j) and qPCR data confirm SpPtf1a upregulation (Fig. 4k). Additionally, in SpPtf1 morphants there is a modest increase of SpMist1 expression (Fig. 4k), and a slight increase of SpFng transcripts (Fig. 4i, j and k). Particularly, SpFng and SpPtf1 seem to be expressed by adjacent and not overlapping cells (Fig. 4j, insert). To test if the effects of SpPtf1a MO were specific, we injected a second MO to block SpPtf1a translation (SpPtf1a MO_2) and we checked the expression of SpCpa2L, since this gene is the first TD marker of this cell type to appear and also the most abundant at pluteus stage. Also in this case, we found that the SpCpa2L mRNA is downregulated in the morphants (Fig. 4c and d), confirming that the two morpholinos produce comparable effects.
SpPtf1a function in mammalian cells
Ptf1a is part of the PTF1 heterotrimeric complex that also includes one of the common E proteins (E47) and RBPL . To test the ability of the sea urchin Ptf1a protein to substitute for mammalian Ptf1a in this complex, we performed luciferase essays in which either SpPtf1a or rat Ptf1a were co-transfected with mammalian E47 and RBPL, and evaluated for an ability to transactivate a PTF1-responsive luciferase reporter containing four tandem repeats of the PTF1 binding site from the rat chymotrypsinogen promoter (Obata et al., 2001). Figure 5a and b show that, in both 293 t and HeLa cells, SpPtf1a can effectively partner with mammalian PTF1 components to activate the rat chymotrypsinogen element, albeit with lower activity than rat Ptf1a (Fig. 5b). As a correlate of this conserved function, we hypothesized that the sea urchin orthologs of mammalian PTF1a target genes would also display similar PTF1a regulatory motifs. In the promoter of genes encoding mammalian pancreatic digestive enzymes, PTF1 has been found to bind to an extended bipartite DNA sequence that consists of an E-box and a TC-box, spaced one or two helical DNA-turns apart. Masui and coauthors have reported a list of PTF1-binding sites of the acinar digestive enzyme gene promoters for mouse . We used the list of binding sites from that paper and built a position weight matrix for the E-box and TC-box, taking into consideration of the distance between these two elements. This matrix was used to look for Ptf1a binding sites in the promoters of the genes that were affected in SpPtf1a morpholino perturbations. Interestingly, we identified canonical Ptf1a binding sites in sea urchin genes displaying in vivo spPtf1a dependence, including SpFng, SpCpa2L, SpPnlp2/5 and SpAmy3. An alignment of the documented Ptf1a binding site in the rat Ctrb1 gene with the putative Ptf1a binding sites in each of the S. purpuratus genes is shown in Fig. 5c.
Pancreatic exocrine-like gene expression response to different feeding regimes
We next asked if pancreatic exocrine-like cells respond to different feeding regimes. To this aim, S. purpuratus embryos were grown and fed on algae until day 10 to reach the normal larval morphology, then food was removed from the culture media. The culture was split into two batches, of which one was fed for 24 h and the other was fasted for the same time interval. First, we observed larval morphology and found that the different feeding regimes show a clear effect on gut sizes. The stomach and intestine of larvae that were fed is bigger than the ones of the fasted larvae, suggesting that the whole gut epithelium responds to the presence of food (Fig. 6a, b). Second, we compared by qPCR the expression levels of SpPtf1a and the previously described digestive enzymes in feeding vs. starving conditions (Fig. 6c). The analysis shows that the expression of the digestive enzymes increases upon feeding. Differently from its putative target genes, SpPtf1a expression is higher in fasted larvae than in fed ones. To understand if the increased expression of those genes was due to the same pancreatic exocrine-like cells or to other cells, we tested gene expression through FISH analysis, confirming that the spatial expression of the above genes does not change in feeding or starving conditions (data not shown).
In the present study, we show that, in the sea urchin larva, cells specialized in food digestion are clustered in the upper stomach. These cells display a differentiation pathway similar to zymogen cells, where regulative genes (such as Ptf1a, Hnf1 and Mist1) and terminal differentiation genes (those encoding digestive enzymes) are co-expressed. Within this upper stomach region, only a few cells on the ventral side do not express the pancreatic markers, denoting a different gene regulatory module and a different function. Figure 7a presents a schematic view of the ontogenetic process leading to the formation of pancreatic exocrine-like clustered cells in the sea urchin embryo and larva. From a temporal point of view, the increase of SpPtf1a expression is necessary for the following expression of at least three digestive enzymes that have been analyzed in this study. The expression of these terminal differentiation genes remains constantly high throughout the larval stage, underlying their importance for larval development. SpPtf1a and SpMist1 expression oscillate, and after initial activation, decrease, suggesting that these early TFs may be necessary to activate digestive enzyme transcription but also other factors are involved in maintaining the right level of enzymes during feeding. Interestingly, the pancreatic exocrine-like cells also express the islet-specific microRNA-375. A former study showed that also miR-100, miR-125 and let-7, markers of neurosecretory endocrine cells, are localized in the same area of the sea urchin stomach . These cells seem to be of the pancreatic exocrine-like type, although further analysis will be necessary to test co-expression.
In the developing pancreas, Notch signaling represses endocrine development and exocrine cell differentiation [16, 19, 47]. The early role of Notch signaling in the sea urchin embryos has been extensively studied [48–51]. In this study, we identified a putative recruitment of Notch signaling at the late gastrula stage in blocking the differentiation of stomach cells towards an exocrine-like fate in the pluteus larva. We found that Notch signaling represses SpPtf1a expression in those cells that are adjacent to the SpPtf1a + cells. Moreover, we observed an ectopic expression of SpPtf1a in a group of cells of the anus, suggesting that the Notch signaling could be responsible for the direct or indirect repression of SpPtf1a expression in these intestinal cells. Nonetheless, we observed that when Notch signaling is blocked, SpMist1 ectopically accumulates in the cells of the apical organ that are adjacent to SpMist1+ cells. One possibility is that a secretory fate is repressed by a Notch lateral-inhibition mechanism in ectodermal cells that are committed to a not-secretory fate. Alternatively, this ectodermal effect may be due to an extra round of cell division. In a similar way, our observations suggest that Notch signaling prevents the exocrine fate of a group of cells of the esophagus. In some cases, we observed that when Notch signaling is blocked, pancreatic genes that are expressed in both the ectoderm and the endoderm lose expression in a specific territory, and not in the other. For instance, the transcription of SpPtf1a in the ectoderm and that of SpMist1 in the stomach do not change in DAPT treated plutei, similarly to what Materna and Davidson observed for SpShr2 and SpDelta in the early embryo . Our data suggest that Notch signaling could have different effects not only during early development, but it can also finely regulate cell differentiation during development in the sea urchin larva. Nonetheless, to deeply analyze the fine mechanisms that underlie regulation of midgut fates, further knockdown experiments of the Notch signaling machinery will be important to enhance our understanding in the future.
The Ptf1a transcriptional complex controls maturation of pancreatic exocrine cells and many efforts has been made to identify its cofactors and targets [52, 53]. In the pancreas, Hnf1a activates the transcription of Ptf1a, which is part of a transcriptional complex that controls pancreatic exocrine cell maturation and regulates the production of digestive enzymes in acinar cells [13, 54]. Here we found that this gene regulatory module is conserved in an early-branching deuterostome. Based on the homology with mammalian binding sequences, we bioinformatically identified several putative Ptf1a binding sites on the promoters of the sea urchin digestive enzymes which transcription is strongly downregulated in Ptf1a perturbed embryos. We also showed that SpPtf1a can bind mammalian target genes to activate transcription. However, further in vivo experiments are necessary to test the functionality of these binding sites predicted in silico. In this study, we also demonstrated that SpPtf1a represses its own gene expression in some stomach cells. We reported that in SpPtf1a morphants, more SpPtf1a + cells appear as single pouches at the ventral and dorsal stomach in cells where normally the gene is never expressed. These cells appear as a distinct group of cells that alternate with cells where SpPtf1a is absent, like a chessboard (Fig. 4i). Conversely, blocking the Notch pathway results in a continuous ring of Ptf1a + cells (Fig. 4l). Hence, the presence of different mechanisms supports the role of a Notch lateral-inhibition model in the stomach domain. Another finding that emerged from our data is that SpPtf1a may repress the expression of SpFng. Due to its role, the Fringe ortholog is widely expressed in the embryos, and this could explain the borderline qPCR values in SpPtf1a morphants. However, our bioinformatic analysis found binding sites for Ptf1a in the SpFng promoter region. To understand if the effect of DAPT treatment on SpFng is specific for an exocrine cell fate, further experiments should define whether Ptf1a interacts with corepressor proteins, as in the case of Nkx6 in mammals , to block SpFng transcription and consequentially Notch activation, in a negative feedback loop manner.
A summary of the above described putative gene interactions and the proposed gene model leading to specification and differentiation of sea urchin pancreatic exocrine-like cells are schematized in Fig. 7b.
The sea urchin larva is a very plastic organism, able to adapt its metabolic and growth programs to changes in the environmental energy context, e.g. arms and ciliary band growth is influenced by different food availability [56–61]. Feeding has a central role in the life of sea urchin larvae, since they may live up to several months in the plankton, therefore they have to be considered as fully developed animals, even if their body axis will be rearranged after metamorphosis. In a previous study, we found that in the sea urchin larva, a member of the insulin family is expressed in a feeding-dependent fashion, opening new questions on whether other genes are influenced by food assumption or deprivation . As expected during feeding, when more digestive enzymes are necessary, we found that the transcription of a carboxipeptidase, a pancreatic lipase and an amylase are significantly increased, whereas SpPtf1a expression showed an opposite behavior. It will be important to further investigate which other factors are involved in the binding of the Ptf1a complex to its targets. Previous transcriptomic analysis reported that SpCpa2L and SpPnlp2/5 transcripts are strongly reduced when SpLox is knocked down  and it is known that its vertebrate ortholog, Pdx1, regulates the development of digestive enzyme-producing acinar cells and insulin expression [63, 64]. In the sea urchin, SpLox is expressed in a group of cells between the lower stomach and the intestine (the pyloric sphincter) , in a gut domain and in cells that are different from the herein characterized pancreatic exocrine-like cells. Thus, a possible scenario that could be worth testing in further experiments is that SpLox could act through a hypothetic signaling pathway to eventually regulate the expression of digestive enzymes, accordingly with food availability.
The pancreatic exocrine-like cells that we describe are localized in the same area where there are also endocrine-like cells producing a member of the insulin family, ILP1 . We cannot rule out the possibility that these two cell types are just in close proximity and do not represent a unique cell-type sharing both functions. Generation of antibodies for the digestive enzymes will be necessary to test if these exocrine markers are coexpressed in the ILP+ cells, but we can at least conclude that both cell types are localized in the same area of the stomach. In the cephalocordate amphioxus, endocrine and exocrine cells do not exhibit any gut specialization [8, 66–68]. In the hagfish, an islet-like organ is independent of exocrine acini, and, in the teleost, endocrine islets are scattered among the more widespread exocrine parenchyma . Thus, two possible evolutionary scenarios could be put forward. First, the pancreas and the exocrine and endocrine-like cells (closely localized) in the sea urchin gut are the result of convergent evolution, thus implying that this feature was lost in early branching chordates. Second, the clustering of the sea urchin pancreatic exocrine-like cells could be the result of an independent evolution from an ancestral pancreatic cell-type. To test these hypotheses and to infer the evolutionary origin of the pancreas, homologs of pancreatic genes should be studied in other echinoderms as well as in other non-chordate deuterostomes, such as the hemichordates.
Here, we report that a gene regulatory module required for vertebrate pancreatic exocrine cell development is shared by a non-chordate deuterostome. Our study, together with previous studies [65, 70], supports the idea that the sea urchin gut is a multifunctional organ and it is composed of highly differentiated cells, which in many aspects are similar to vertebrates, therefore providing an amenable system in which to study development and function of gut cells.
Chera S, de Rosa R, Miljkovic-Licina M, Dobretz K, Ghila L, Kaloulis K, Galliot B. Silencing of the hydra serine protease inhibitor Kazal1 gene mimics the human SPINK1 pancreatic phenotype. J Cell Sci. 2006;119:846–57.
Lentz TL. Intramitochondrial glycogen granules in digestive cells of Hydra. J Cell Biol. 1966;29:162–7.
Goldberg WM. Gastrodermal structure and feeding responses in the scleractinian Mycetophyllia reesi, a coral with novel digestive filaments. Tissue Cell. 2002;34:246–61.
Brennan GP, Ramasamy P. Ultrastructure of the gut caecal epithelium of Pricea multae (Monogenea: Polyopisthocotylea). Parasitol Res. 1996;82:312–8.
Jones MK, Hughes-Stamm SR, East RM, Cribb TH. Ultrastructure of the digestive tract of Gyliauchen nahaensis (Platyhelminthes, Digenea), an inhabitant of the hindgut of herbivorous fishes. J Morphol. 2000;246:198–211.
Brennan GP, Ramasamy P. Ultrastructure of the surface structures and electron immunogold labeling of peptide immunoreactivity in the nervous system of Pseudothoracocotyla indica (Polyopisthocotylea: Monogenea). Parasitol Res. 1996;82:638–46.
Lemaitre B, Miguel-Aliaga I. The digestive tract of Drosophila melanogaster. Annu Rev Genet. 2013;47:377–404.
Slack JM. Developmental biology of the pancreas. Development. 1995;121:1569–80.
Krapp A, Knöfler M, Ledermann B, Bürki K, Berney C, Zoerkler N, Hagenbüchle O, Wellauer PK. The bHLH protein PTF1-p48 is essential for the formation of the exocrine and the correct spatial organization of the endocrine pancreas. Genes Dev. 1998;12:3752–63.
Lemercier C, To RQ, Swanson BJ, Lyons GE, Konieczny SF. Mist1: a novel basic helix-loop-helix transcription factor exhibits a developmentally regulated expression pattern. Dev Biol. 1997;182:101–13.
Cockell M, Stevenson BJ, Strubin M, Hagenbüchle O, Wellauer PK. Identification of a cell-specific DNA-binding activity that interacts with a transcriptional activator of genes expressed in the acinar pancreas. Mol Cell Biol. 1989;9:2464–76.
Kawaguchi Y, Cooper B, Gannon M, Ray M, MacDonald RJ, Wright CVE. The role of the transcriptional regulator Ptf1a in converting intestinal to pancreatic progenitors. Nat Genet. 2002;32:128–34.
Krapp A, Knöfler M, Frutiger S, Hughes GJ, Hagenbüchle O, Wellauer PK. The p48 DNA-binding subunit of transcription factor PTF1 is a new exocrine pancreas-specific basic helix-loop-helix protein. EMBO J. 1996;15:4317–29.
Pin CL, Rukstalis JM, Johnson C, Konieczny SF. The bHLH transcription factor Mist1 is required to maintain exocrine pancreas cell organization and acinar cell identity. J Cell Biol. 2001;155:519–30.
Haumaitre C, Barbacci E, Jenny M, Ott MO, Gradwohl G, Cereghini S. Lack of TCF2/vHNF1 in mice leads to pancreas agenesis. Proc Natl Acad Sci U S A. 2005;102:1490–5.
Hald J, Hjorth JP, German MS, Madsen OD, Serup P, Jensen J. Activated Notch1 prevents differentiation of pancreatic acinar cells and attenuate endocrine development. Dev Biol. 2003;260:426–37.
Kim SK, Hebrok M. Intercellular signals regulating pancreas development and function. Genes Dev. 2001;15:111–27.
Norgaard GA, Jensen JN, Jensen J. FGF10 signaling maintains the pancreatic progenitor cell state revealing a novel role of Notch in organ development. Dev Biol. 2003;264:323–38.
Esni F, Ghosh B, Biankin AV, Lin JW, Albert MA, Yu X, MacDonald RJ, Civin CI, Real FX, Pack MA, Ball DW, Leach SD. Notch inhibits Ptf1 function and acinar cell differentiation in developing mouse and zebrafish pancreas. Development. 2004;131:4213–24.
Poy MN, Eliasson L, Krutzfeldt J, Kuwajima S, Ma X, Macdonald PE, Pfeffer S, Tuschl T, Rajewsky N, Rorsman P, Stoffel M. A pancreatic islet-specific microRNA regulates insulin secretion. Nature. 2004;432:226–30.
Kloosterman WP, Lagendijk AK, Ketting RF, Moulton JD, Plasterk RHA. Targeted inhibition of miRNA maturation with morpholinos reveals a role for miR-375 in pancreatic islet development. PLoS Biol. 2007;5:e203.
Poy MN, Hausser J, Trajkovski M, Braun M, Collins S, Rorsman P, Zavolan M, Stoffel M. miR-375 maintains normal pancreatic alpha- and beta-cell mass. Proc Natl Acad Sci U S A. 2009;106:5813–8.
McClay DR. Evolutionary crossroads in developmental biology: sea urchins. Development. 2011;138:2639–48.
Burke RD. Structure of the digestive tract of the pluteus larva of Dendraster excentricus (Echinodermata: Echinoida). Zoomorphology. 1981;98:209–25.
Vacquier VD. The appearance of -1,3-glucanohydrolase activity during the differentiation of the gut of sand dollar plutei. Dev Biol. 1971;26:1–10.
Vacquier VD, Korn LJ, Epel D. The appearance of -amylase activity during gut differentiation in sand dollar plutei. Dev Biol. 1971;26:393–9.
Vacquier VD. The effects of glucose and lithium chloride on the appearance of -1,3-glucanohydrolase activity in sand dollar plutei. Dev Biol. 1971;26:11–6.
Holland ND, Lauritis JA: The fine structure of the gastric exocrine cells of the purple sea urchin, Strongylocentrotus purpuratus. Trans Amer Microscop Soc. 1968;87:201-209.
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.
Andrikou C, Iovene E, Rizzo F, Oliveri P, Arnone MI. Myogenesis in the sea urchin embryo: the molecular fingerprint of the myoblast precursors. Evodevo. 2013;4:33.
Andrikou C, Pai C-Y, Su Y-H, Arnone MI: Logics and properties of a genetic regulatory program that drives embryonic muscle development in an echinoderm. Elife. 2015;4.
Hughes JN, Dodge N, Rathjen PD, Rathjen J. A novel role for gamma-secretase in the formation of primitive streak-like intermediates from ES cells in culture. Stem Cells. 2009;27:2941–51.
Rast JP, Amore G, Calestani C, Livi CB, Ransick A, Davidson EH. Recovery of developmentally defined gene sets from high-density cDNA macroarrays. Dev Biol. 2000;228:270–86.
Nemer M, Rondinelli E, Infante D, Infante AA. Polyubiquitin RNA characteristics and conditional induction in sea urchin embryos. Dev Biol. 1991;145:255–65.
Dong PDS, Provost E, Leach SD, Stainier DYR. Graded levels of Ptf1a differentially regulate endocrine and exocrine fates in the developing pancreas. Genes Dev. 2008;22:1445–50.
Howard-Ashby M, Materna SC, Brown CT, Chen L, Cameron RA, Davidson EH. Gene families encoding transcription factors expressed in early development of Strongylocentrotus purpuratus. Dev Biol. 2006;300:90–107.
Howard-Ashby M, Materna SC, Brown CT, Chen L, Cameron RA, Davidson EH. Identification and characterization of homeobox transcription factor genes in Strongylocentrotus purpuratus, and their expression in embryonic development. Dev Biol. 2006;300:74–89.
Bach I, Mattei MG, Cereghini S, Yaniv M. Two members of an HNF1 homeoprotein family are expressed in human liver. Nucleic Acids Res. 1991;19:3553–9.
Booth H, Holland P: Annotation, nomenclature and evolution of four novel homeobox genes expressed in the human germ line. Gene. 2007;387:7–14.
Poustka AJ, Groth D, Hennig S, Thamm S, Cameron A, Beck A, Reinhardt R, Herwig R, Panopoulou G, Lehrach H. Generation, annotation, evolutionary analysis, and database integration of 20,000 unique sea urchin EST clusters. Genome Res. 2003;13:2736–46.
Annunziata R, Perillo M, Andrikou C, Cole AG, Martinez P, Arnone MI. Pattern and process during sea urchin gut morphogenesis: the regulatory landscape. Genesis. 2014;52:251–68.
Materna SC, Davidson EH. A comprehensive analysis of Delta signaling in pre-gastrular sea urchin embryos. Dev Biol. 2012;364:77–87.
Moloney DJ, Panin VM, Johnston SH, Chen J, Shao L: Fringe is a glycosyltransferase that modifies Notch. Nature. 2000;406:369–375.
Peterson RE, McClay DR. A Fringe-modified Notch signal affects specification of mesoderm and endoderm in the sea urchin embryo. Dev Biol. 2005;282:126–37.
Masui T, Long Q, Beres TM, Magnuson MA, MacDonald RJ. Early pancreatic development requires the vertebrate Suppressor of Hairless (RBPJ) in the PTF1 bHLH complex. Genes Dev. 2007;21:2629–43.
Christodoulou F, Raible F, Tomer R, Simakov O, Trachana K, Klaus S, Snyman H, Hannon GJ, Bork P, Arendt D. Ancient animal microRNAs and the evolution of tissue identity. Nature. 2010;463:1084–8.
Murtaugh LC, Stanger BZ, Kwan KM, Melton DA. Notch signaling controls multiple steps of pancreatic differentiation. Proc Natl Acad Sci U S A. 2003;100:14920–5.
Croce JC, McClay DR. Dynamics of Delta/Notch signaling on endomesoderm segregation in the sea urchin embryo. Development. 2010;137:83–91.
Sherwood DR, McClay DR: Identification and localization of a sea urchin Notch homologue: insights into vegetal plate regionalization and Notch receptor regulation. Development. 1997;124:3363–3374.
Sherwood DR, McClay DR: LvNotch signaling plays a dual role in regulating the position of the ectoderm-endoderm boundary in the sea urchin embryo. Development. 2001;128:2221–32.
Materna SC, Swartz SZ, Smith J. Notch and Nodal control forkhead factor expression in the specification of multipotent progenitors in sea urchin. Development. 2013;140:1796–806.
Thompson N, Gésina E, Scheinert P, Bucher P, Grapin-Botton A. RNA profiling and chromatin immunoprecipitation-sequencing reveal that PTF1a stabilizes pancreas progenitor identity via the control of MNX1/HLXB9 and a network of other transcription factors. Mol Cell Biol. 2012;32:1189–99.
Arda HE, Benitez CM, Kim SK. Gene regulatory networks governing pancreas development. Dev Cell. 2013;25:5–13.
Petrucco S, Wellauer PK, Hagenbüchle O. The DNA-binding activity of transcription factor PTF1 parallels the synthesis of pancreas-specific mRNAs during mouse development. Mol Cell Biol. 1990;10:254–64.
Schaffer AE, Freude KK, Nelson SB, Sander M. Nkx6 transcription factors and Ptf1a function as antagonistic lineage determinants in multipotent pancreatic progenitors. Dev Cell. 2010;18:1022–9.
Adams DK, Sewell MA, Angerer RC, Angerer LM. Rapid adaptation to food availability by a dopamine-mediated morphogenetic response. Nat Commun. 2011;2:592.
Heyland A, Reitzel AM, Hodin J. Thyroid hormones determine developmental mode in sand dollars (Echinodermata: Echinoidea). Evol Dev. 2004;6:382–92.
Reitzel AM, Heyland A: Reduction in morphological plasticity in echinoid larvae: relationship of plasticity with maternal investment and food availability. Evol Ecol Res. 2007;9:109–121.
Strathmann RR, Fenaux L, Strathmann MF: Heterochronic developmental plasticity in larval sea urchins and its implications for evolution of nonfeeding larvae. Evolution. 1992;46:972–986.
Strathmann RR, Foley GP, Hysert AN. Loss and gain of the juvenile rudiment and metamorphic competence during starvation and feeding of bryozoan larvae. Evol Dev. 2008;10:731–6.
Carrier TJ, King BL, Coffman JA. Gene Expression Changes Associated With the Developmental Plasticity of Sea Urchin Larvae in Response to Food Availability. Biol Bull. 2015;228:171–80.
Annunziata R, Martinez P, Arnone MI. Intact cluster and chordate-like expression of ParaHox genes in a sea star. BMC Biol. 2013;11:68.
Ohlsson H, Karlsson K, Edlund T. IPF1, a homeodomain-containing transactivator of the insulin gene. EMBO J. 1993;12:4251–9.
Zhou Q, Law AC, Rajagopal J, Anderson WJ, Gray PA, Melton DA. A multipotent progenitor domain guides pancreatic organogenesis. Dev Cell. 2007;13:103–14.
Cole AG, Rizzo F, Martinez P, Fernandez-Serra M, Arnone MI. Two ParaHox genes, SpLox and SpCdx, interact to partition the posterior endoderm in the formation of a functional gut. Development. 2009;136:541–9.
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.
Lecroisey C, Le Pétillon Y, Escriva H, Lammert E, Laudet V. Identification, evolution and expression of an insulin-like peptide in the cephalochordate Branchiostoma lanceolatum. PLoS One. 2015;10:e0119461.
Pieler T, Chen Y. Forgotten and novel aspects in pancreas development. Biol Cell. 2006;98:79–88.
Youson JH, Al-Mahrouki AA. Ontogenetic and phylogenetic development of the endocrine pancreas (islet organ) in fish. Gen Comp Endocrinol. 1999;116:303–35.
Annunziata R, Arnone MI. A dynamic regulatory network explains ParaHox gene control of gut patterning in the sea urchin. Development. 2014;141:2462–72.
The authors wish to thank Giovanna Benvenuto for confocal microscopy, Marco Borra for quantitative Real Time PCR and Xiaogang Zhong for bioinformatics assistance. We would also like to thank Pat Leahy (KML, Caltech, Pasadena, CA) for providing adult S. purpuratus and Davide Caramiello (SZN) for animal maintenance. MP would like to thank Steven Z. Swartz for stimulating discussions and paper revisions. This work was supported by MIUR (premiale PANTRAC to MIA) and NIH grant DK61215 (to SDL). M.P. has been supported by a SZN PhD fellowship and by a fellowship from POR Campania FSE 2007-2013 Project-MODO, Model Organism.
MIA and MP conceived the project, SDL and YJW designed the bioinformatic analysis and the experiments in mammalian cells. MP performed the experimental work, YJW performed the luciferase essays. The authors contributed equally to the interpretation of the results and the writing of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Supplementary Table 1. (PDF 55 kb)