Document Detail

Fibroblast growth factor signalling controls nervous system patterning and pigment cell formation in Ciona intestinalis.
Jump to Full Text
MedLine Citation:
PMID:  25189217     Owner:  NLM     Status:  In-Data-Review    
Abstract/OtherAbstract:
During the development of the central nervous system (CNS), combinations of transcription factors and signalling molecules orchestrate patterning, specification and differentiation of neural cell types. In vertebrates, three types of melanin-containing pigment cells, exert a variety of functional roles including visual perception. Here we analysed the mechanisms underlying pigment cell specification within the CNS of a simple chordate, the ascidian Ciona intestinalis. Ciona tadpole larvae exhibit a basic chordate body plan characterized by a small number of neural cells. We employed lineage-specific transcription profiling to characterize the expression of genes downstream of fibroblast growth factor signalling, which govern pigment cell formation. We demonstrate that FGF signalling sequentially imposes a pigment cell identity at the expense of anterior neural fates. We identify FGF-dependent and pigment cell-specific factors, including the small GTPase, Rab32/38 and demonstrated its requirement for the pigmentation of larval sensory organs.
Authors:
Claudia Racioppi; Ashwani K Kamal; Florian Razy-Krajka; Gennaro Gambardella; Laura Zanetti; Diego di Bernardo; Remo Sanges; Lionel A Christiaen; Filomena Ristoratore
Related Documents :
2483307 - Lymphoid cells, small cell lung cancer cells and epithelial cells share a membrane dete...
3888047 - Characterisation of a cell line from a human lung lymphosarcoma. abundance of vimentin-...
18313807 - 2,3-disubstituted 8-arylamino-3h-imidazo[4,5-g]quinazolines: a novel class of antitumor...
9000807 - Frequency and distribution of chromosomal integration sites of the epstein-barr virus g...
6616447 - Induction of transformation by six classes of chemical carcinogens in adult rat liver e...
1489917 - Locomotion forces generated by a polymorphonuclear leukocyte.
Publication Detail:
Type:  Journal Article     Date:  2014-09-05
Journal Detail:
Title:  Nature communications     Volume:  5     ISSN:  2041-1723     ISO Abbreviation:  Nat Commun     Publication Date:  2014  
Date Detail:
Created Date:  2014-09-05     Completed Date:  -     Revised Date:  -    
Medline Journal Info:
Nlm Unique ID:  101528555     Medline TA:  Nat Commun     Country:  England    
Other Details:
Languages:  eng     Pagination:  4830     Citation Subset:  IM    
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Nat Commun
Journal ID (iso-abbrev): Nat Commun
ISSN: 2041-1723
Publisher: Nature Pub. Group
Article Information
Download PDF
Copyright © 2014, Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.
open-access:
Received Day: 02 Month: 12 Year: 2013
Accepted Day: 25 Month: 07 Year: 2014
Electronic publication date: Day: 05 Month: 09 Year: 2014
Volume: 5E-location ID: 4830
PubMed Id: 25189217
ID: 4164782
Publisher Item Identifier: ncomms5830
DOI: 10.1038/ncomms5830

Fibroblast growth factor signalling controls nervous system patterning and pigment cell formation in Ciona intestinalis
Claudia Racioppi1
Ashwani K. Kamal1
Florian Razy-Krajka2
Gennaro Gambardella3
Laura Zanetti1
Diego di Bernardo3
Remo Sanges4
Lionel A. Christiaena2
Filomena Ristoratoreb1
1Cellular and Developmental Biology Laboratory, Stazione Zoologica Anton Dorhn, Villa Comunale, 80121 Napoli, Italy
2Center for Developmental Genetics, Department of Biology, New York University, New York, New York 1003, USA
3Telethon Institute of Genetics and Medicine (TIGEM), Via P. Castellino 111, 80131 Naples, Italy
4Animal Physiology and Evolution Laboratory, Stazione Zoologica Anton Dorhn, Villa Comunale, 80121 Napoli, Italy
alc121@nyu.edu
bfilomena.ristoratore@szn.it

Distinct combinations of transcription factors (TFs) and signalling molecules define the complex gene regulatory networks that control patterning and commitment of different cell types. In vertebrate embryos, pigment-producing cells originate from the closing neuroepithelium: Retinal pigment epithelium (RPE) of the eye and pigment cells of the pineal organ arising from the neural tube; melanocytes of the inner ear, skin, derive from a population of cells migrating from the neural plate border (neural crest)1. Despite their different embryological origins, all vertebrate pigment cells share some basic characteristics: melanin pigments are synthesized and stored in melanosomes involving tyrosinase (tyr) and tyr-related enzymes. Melanosomes are tissue-specific lysosome-related organelles distributed in the cell by a regulated vesicular trafficking system, which requires several Rab-type small GTPases and their regulators2, 3.

The development of vertebrate RPE and neural retina involves complex interplays among cell-signalling and cell-intrinsic factors: Wnt, bone morphogenetic proteins and fibroblast growth factors (FGFs)4, 5. Among the documented roles for signalling pathways in the formation of the vertebrate visual system, FGF is required to induce neural retina progenitors at the expense of RPE formation4, 6. Nevertheless, understanding the molecular impacts of cell-extrinsic components on cell-specific determinants in the developing of the retina remains a challenging task.

Ascidians belonging to the chordate group represent an attractive model system to investigate the generation of cell diversity in the central nervous system (CNS). Despite their highly divergent adult body plan and genomic architecture, phylogenomic analyses indicated that tunicates are the closest living relatives of the vertebrates7. The reduced genetic redundancy found in the compact tunicate genome facilitates functional analyses of homologous TFs and signalling molecules in the context of the typical chordate body plan of embryos and larvae.

The ascidian Ciona intestinalis emerged as an experimental system of choice due to their simple embryos, rapid development and straightforward transgenesis methods8, 9. Because the ascidian cell lineages and fate maps are well characterized10, the identification of neural precursors starting from early fate restriction facilitates the study of progressive cell fate specification in individual precursors through each division. Thus, it is possible to investigate the cellular and molecular mechanisms underlying fate specification with single-cell resolution during the development of chordate CNS.

The Ciona tadpole larva offers the advantage of a typical chordate body plan with a simple brain11. The anterior sensory vesicle contains two distinct pigment cell containing sensory organs: the geotactic otolith, a single cell containing a melanin granule12 and the photosensitive ocellus, composed of 30 photoreceptor cells, three lens cells and one cup-shaped pigment cell13 Ascidian pigment cells express the melanogenic enzyme-coding genes tyrosinase (tyr) and tyr-related protein14, in a fashion similar to the vertebrate pigment cells15. Moreover, Ciona cephalic melanocyte lineage expresses neural plate border and neural crest specification genes, like Id, Snail, Ets and FoxD and the signalling events and regulatory inputs leading to ascidian ocellus formation are similar to those involved in vertebrate neural crest lineage specification, including canonical Wnt and FoxD-mediated repression of Mitf16.

We recently showed that the FGF/MAPK/Ets pathway is responsible for pigmented cell precursor specification, starting from the late gastrula stage, by rendering their progenies competent to respond to a Wnt signal17. However, the FGF signalling downstream targets involved in the regulatory architecture of pigment cell lineage formation remain uncharacterized. Here, we employed a transcriptomic approach by exploiting targeted interference of the FGF signalling coupled with the enrichment of cells of interest by fluorescence-activated cell sorting (FACS) to obtain pigment cell lineage-specific transcription profiles. By using ‘reverse engineering’ methods, we created a pigment lineage-specific co-regulated network to identify key effectors for pigmentation. By characterizing the transcriptional response to FGF signalling in the pigment cell lineage, we identified molecular mechanisms underlying anterior–posterior patterning of Ciona CNS and identified novel regulators of pigment cell specification and differentiation.


Results
FGF dependent Pigment cell precursors transcription profiles

In ascidians, the pigment cells of the otolith and ocellus arise from the paired a8.25 blastomeres10, 18. During gastrulation, these blastomeres divide and form a9.49 and a9.50 cell pairs: the a9.50s (located in the row IV of the neural plate) are progenitor cells of the anterior brain, while the a9.49s (located in the row III) become fate restricted as pigment cell precursors (PCPs) (Fig. 1a)17. At the mid-neurula stage (stage 15; Four-dimensional Ascidian Body Atlas19), bilateral PCPs divide, forming two cell pairs (a10.98s and a10.97s, Fig. 1b,c). Taking advantage of the stability and brightness of 2xGFP reporter driven by tyrp1/2a enhancer (tyrp1/2a>2xGFP), we observed that a10.97/a10.98s, previously described as a postmitotic pigment cell precursor20, undergo an additional division forming eight cells (four on each side) (Fig. 1d). An extra a10.97 cell division was also suggested in a note added in proof of Abitua et al.16 and recently confirmed in Haupaix et al. (see note added in proof.).

Following cell nomenclature previously described in Cole and Meinertzhagen20, we referred to these cells as a11.196 and a11.195 for the two anterior cells as a11.194 and a11.193 for the two posterior ones (Fig. 1d,e). Among the tyrp1/2a+cells, only the two most posterior become otolith and ocellus pigment cells17.

We recently demonstrated that FGF signalling induces PCP’s formation17. We perturbed the endogenous FGF signalling by using a dominant-negative form of the unique Ci-FGF receptor (FGFRDN) and a constitutively active form of Ci-Ets1/2 (Ets:Vp16) (ref. 21), a transcriptional effector of FGF/MAP kinase cascade22 that is predominantly expressed in all the row III cells, including PCPs, at late gastrula stage17. Targeted expression of FGFRDN or Ets:Vp16 in PCPs caused larvae to either lack pigment cells or form extra pigment cells in the sensory vesicle, respectively17 (Fig. 2a–c).

Here, we employed FACS and microarray analyses to obtain PCP lineage-specific transcription profiles and identify FGF downstream genes responsible for Ciona PCPs formation. We used the same tyrp1/2a enhancer14 to express FGFRDN, Ets:VP16 and green fluorescent protein (GFP) and FAC-sort a9.49 cells and their progeny as GFP+ cells (Fig. 2d,e). To avoid a contamination by mesenchymal cells where Ciona enhancers are often ectopically activated, we co-electroporated the MyoD905>YFP (yellow fluorescent protein) construct and counter-selected YFP+ cells as previously described23, 24 (Fig. 2e). To gain insight into the temporal deployment of the FGF signalling regulatory network, we examined the transcription profiles of PCPs isolated at two developmental stages: 8 hpf at ~16 °C (neurula stage) corresponding to early FGF-mediated induction, and 12 hpf at ~18 °C tailbud(TB) stage, where PCPs are already fate restricted as pigment cells (Fig. 2d). The entire data set of 24 microarrays was analysed to identify new candidates involved in the PCP-specific FGF/MAPK/Ets-dependent regulatory network (Supplementary Fig. 1). Correlation analyses of normalized expression values among replicates indicated that microarray experiments were highly reproducible (Supplementary Fig. 1).

To explore the expression data set, we reasoned that pigment cell-specific genes should be downregulated upon misexpression of FGFRDN and/or upregulated by Ets:Vp16 compared with controls (Fig. 2h), given that FGFRDN blocks pigmentation while Ets:Vp16 causes the formation of supernumerary pigment cells in the PCP lineage17 (Fig. 2a–c). The quadruplicate experiments were used to compute adjusted P values, corrected for multiple testing errors.

We first examined the expression values for a training data set consisting of known PCP markers, including Ci-tyr, Ci-tyrp1/2a14, Ci-Tcf/Lef17, Ci-Mitf, Ci-Bmp5/7 and Ci-msxb25 (Fig. 2f). When comparing their expression levels in FGFRDN versus Ets:Vp16 12hpf samples, the melanogenic genes, Ci-tyr, Ci-tyrp1/2a, the direct FGF downstream gene and transcriptional effector of canonical Wnt signalling Ci-Tcf/Lef as well as known factors involved in pigment cell formation, Ci-Mitf and Ci-Bmp5/7, were strongly downregulated (Fig. 2f). Ci-Mitf and Ci-Bmp5/7 expressions decrease between 8 and 12 hpf. This is in agreement with their early expression in the pigment cell lineage (Supplementary Fig. 2a,b; ref. 16) and is probably related to their early regulatory function.

PCP expression of the transcription regulator Ci-msxb decreased over time, while FGFRDN increased its expression at 12 hpf. These observations are consistent with in situ hybridization (ISH) assays revealing that Ci-msxb is expressed in the 8 hpf neural plate, including the PCPs, while it is excluded from most of the a9.49 derivatives in 12 hpf tailbud embryos, except for the most anterior a11.195s and a11.196s, as demonstrated by the co-localization with Ci-tyrp1/2a (Supplementary Fig. 2c,d). According to the extra cell division observed for pigment cell lineage, tyrp1/2a transcripts accumulated preferentially in a10.97s rather than a10.98s retaining a stronger expression in a10.97 progeny (a11.193s/a11.94s) than a10.98s derivatives (a11.195s/a11.96s) that maintains only a residual amount of tyrp1/2a RNAs. Thus, our microarray data recapitulated distinct tissue-specific expression dynamics: Ci-msxb is downregulated within the PCP lineage between 8 and 12 hpf, but it is maintained in surrounding non-pigmented CNS cells and therefore ectopically expressed in 12 hpf PCPs on inhibition of FGF signalling.

Taken together, these observations indicate that our microarray data set efficiently captured known PCP gene activities and probably offers an appropriate representation of pigment cell-specific transcription profiles.

FGF signalling opposes the anterior neural plate fate in PCPs

To consider all the experimental conditions used for microarray analysis, we focused on 912 probe sets that were differentially expressed when comparing the two perturbated conditions to each other and with control samples (FGFRDN versus Control, Ets:Vp16 versus Control and FGFRDN versus Ets:Vp16) at 8 and 12 hpf (Fig. 2g).

Since FGF inhibition reprograms PCPs into non-pigmented cells, we first ought to characterize the transcription profiles underlying this fate change. We reasoned that genes upregulated upon FGF inhibition reflect the CNS fate alternative to PCP. We found 461 probe sets corresponding to 264 unique ‘FGF-inhibited’ genes (Fig. 2h). Most of these genes were also more highly expressed at 8 than 12 hpf (Fig. 2h). Considering that CNS regionalization and cell fate specification take place in early development, we expected that genes involved in these processes would be more highly expressed at early stages. Thus, we reasoned that early ‘FGF-inhibited’ genes would reflect the respecification of a9.49 lineage and illuminate the role of FGF signalling in CNS patterning and early PCP specification.

Among the upregulated genes, that is, more highly expressed in FGFRDN samples and/or decreased expression in Ets:Vp16 in comparison with control samples, we first focused on Ci-Six3/6, a known transcription factor-coding gene expressed in anterior neural domains throughout the animal kingdom26, 27. In Ciona, it is expressed at late gastrula stage in the neural plate IV row of cells, including in daughters of the a9.50 cells (Fig. 3a), the sisters of a9.49 PCP precursors located on neural plate III row (Fig. 1a,b and in Moret et al.26). Double whole-mount ISH (WMISH) demonstrated that Ci-Six3/6 is expressed in a territory that is adjacent but not overlapping with the domain of Ci-tyrp1/2a expression at the neurula to tailbud stages (Fig. 3a,b).

To mirror microarray experiments, we performed double WMISH using Ci-Six3/6 and Ci-tyrp1/2a probes following FGF perturbations. In Ets:Vp16 electroporated embryos and control samples, we did not detect any change in Ci-Six3/6 expression (Fig. 3d), while FGF inhibition caused ectopic Ci-Six3/6 expression in a posterior domain that includes Ci-tyrp1/2a-positive a9.49 derivatives (Fig. 3c). These findings recapitulate the microarray data and indicate that FGF signalling normally inhibits Ci-Six3/6 expression in Ci-tyrp1/2a+ PCPs. To characterize the FGF regulatory network that generate cell identity in the anterior CNS of the tailbud embryos, we analysed the expression pattern of genes encoding TFs that were significantly upregulated when comparing FGFRDN with Ets:Vp16 in 12 hpf PCPs. Among 21 genes, we selected Ci-Pax3/7, Ci-Macho-1, Ci-FoxHa, Ci-msxb and Ci-Pax6. We included additional genes with published expression profiles in the ANISEED database28: Ci-RGR/Nut, Ci-Neurogenin (Ngn) and Ci-COE25 that are widely expressed in the CNS. We then asked whether these genes co-localize with Ci-Six3/6+ anterior neural cells in tailbud embryos. All these genes were expressed in anterior neural domains, overlapping with Ci-Six3/6 (Supplementary Fig. 3a–i).

According to published neural lineages20, this domain located anteriorly to the prospective PCP includes a10.100s, a11.198s and a11.197s, all derived from the a9.50 blastomeres, which already express Ci-Six3/6 in early neurula embryos and are a9.49 sister cells (Supplementary Fig. 3a). Our observations that Ci-Six3/6 and other markers of a9.50 derivatives were upregulated by FGFRDN misexpression in the a9.49 lineage indicate that FGF signalling promotes PCP specification in a9.49 derivatives by inhibiting a9.50 specification including Ci-Six3/6 expression, thus specifying different anterior-to-posterior fates within the neural plate.

These results indicate that the most posterior Ci-Six3/6+ cells express a unique combination of genes, which are inhibited by FGF signalling in the adjacent PCPs. As further support, we found that anterior CNS markers, as Ci-Pax3/7, Ci-Pax6, Ci-COE, are ectopically expressed posteriorly in PCPs domain when FGF signalling is specifically inhibited in these cells (Supplementary Fig. 4).

Progressive PCP induction by sequential FGF/MAPK signalling

Starting from early gastrulation, the FGF signal is released from the A-nerve cord precursors and induces the a8.25 blastomeres, subsequentially their progenies (a9.49s), as future pigment cells17, 18, 29. To determine the timing of MAPK activation during PCP divisions in a10.97/a10.98s and in their progeny, we performed immunostaining assays of the dual-phosphorylated (dp) form of ERK1/2 (ref. 29). At late gastrula stage, dpERK appears in the a-line neural plate in the row III cells, including a9.49s, but not on the row IV that includes Ci-Six3/6-expressing cells, as shown in Hudson et al.29 Later in development, we detected dpERK staining in the a10.97s but not in a10.98s (Fig. 4a–c). After the second round of PCP division, dpERK staining remains in a10.97 progenies, appearing stronger in the precursors of otolith and ocellus pigmented cells (a11.193s) than in a11.194s (Fig. 4d). These data suggest that the MAPK activation is progressively restricted to the posterior PCPs, the sole to become pigmented at later stage. We next used timed inhibition of MAPK signalling to further test whether sequential induction of FGF signalling corresponds to successive steps in PCP specification. To this aim, we followed PCP division by assaying GFP reporter expression and placed embryos in the pharmacological agent U0126 to block the MEK kinase for 30 min at gastrula (a9.49s), neurula (a10.97s/a10.98s) and early tailbud (a11.193s/a11.194s and a11.195s/a11.196s) stages. We then observed the presence of pigmented cells in the larval sensory vesicle. As already demonstrated, the resulting larvae completely lacked pigment cells when embryos were treated at gastrula stage17. We found that when embryos were treated at neurula and early tailbud stages, larvae deprived of pigment cells decreased (Fig. 5a) while a little effect was observed when embryos were placed in UO126 after the final cell division (a11.193s) of pigment cell lineage. These data are in according to the differential activation of ERK1/2 in PCPs that appeared stronger in the a10.97s and their progenies than the other PCP cells.

We then examined Ci-Six3/6+ and Ci-tyrp1/2a+ cells at tailbud stage following Ci-Six3/6 and Ci-tyrp1/2a transcripts, together with GFP expression driven by tyrp1/2a enhancer (ptyrp1/2a>2xGFP) as pigment cell lineage marker (Fig. 5b). Upon treatment at 5h30′ (a9.49s), Ci-tyrp1/2a expression is lost in PCPs (Fig. 5c) whereas it is strongly reduced when embryos were treated at the neurula stage (10.97s/a10.98s) (Fig. 5d) and appeared almost normal when the drug was applied at early tailbud stage (a11.193s/a11.194s) (Fig. 5e); by contrast, Ci-Six3/6 detection is always expanded posteriorly in the CNS, although to a lesser extent upon later treatments. The results obtained with successive U0126 treatments show that MAPK activity is necessary at each sequential PCP division (a9.49 and a10.97/a10.98 division) to direct PCPs through their final fate from gastrula to tailbud stage and became less effective when PCP became post mitotic (a11.193s/a11.194s).

Expression of new marker genes for pigment cell lineage

We next ought to identify novel genes expressed in the PCP downstream of FGF signalling. Since FGF signalling is implicated in Ciona pigment cell formation17, we reasoned candidate PCP-specific FGF targets should be either downregulated by FGFRDN and/or upregulated by Ets:Vp16 at 8 or 12 hpf and identified 451 probe sets matched to 207 unique annotated genes (Fig. 2h). Gene ontology (GO) terms analyses indicated that the pool of candidate PCP-specific FGF targets are enriched for genes involved in melanization with overrepresented processes such as ‘pigmentation’, ‘pigment binding’, ‘melanocyte differentiation’ and ‘eye pigment biosynthetic process’ thus confirming the validity of our approach and its potential for the discovery of novel PCP-specific candidates involved in melanization that are under the control of FGF signalling (Supplementary Table 1).

To gain insights into the FGF-dependent regulatory network controlling pigment cell formation, we analysed the expression values of candidate PCP genes encoding TFs and signalling molecules. We identified 20 regulatory genes that might contribute to PCP fate specification (Supplementary Fig. 5). We identified the conserved FGF/MAPK pathway components MAPKK7, Dusp1.2.4.5, Dusp6.9 and ELK1/3 (ref. 31), the BMP-family ligands Bmp5/7, Bmp2/4 (ref. 31) and Wnt signalling molecules such as the downstream effector Tcf/Lef, Wnt14/15, Frizzled receptor (Orphan Fz-a), the Wnt-co-receptors Lrp4 and Lrp1b and the secreted Wnt antagonist Dickkopf1 (Dkk1) (ref. 32). Interestingly, genes encoding the secreted BMP antagonists, Ci-noggin33 as well as Ci-BMPER, a BMP-binding antagonist34, were upregulated in FGFRDN compared with Ets:Vp16 samples. Taken together, these data suggest that the FGF/MAPK pathway potentiates BMP signalling during PCP fate specification and provides further insight into the interaction between FGF and Wnt that we previously uncovered17.

To find new PCP-specific markers, we selected 30 genes among the downregulated genes identified from the pairwise comparisons (FGFRDN versus Ets:Vp16 and FGFRDN versus control) for subsequent ISH analyses. We obtained a specific CNS staining for 14 of the PCP candidates. Among these, six genes (Ci-Lrp4, Ci-Bzrap1, Ci-Casz1, Ci-Mad/mnt, Ci-Slc45a2, Ci-Piwi) showing defined PCP expression were selected for further analysis (Fig. 6a–g). Double WMISH using the pigment cell marker, Ci-tyrp1/2a14, was performed to ascertain PCP-specific expression (Fig. 6). Detailed analysis of the expression domain showed that some genes are expressed in all the eight a9.49-derived cells (Ci-Mad/Mnt, Ci-Piwi), while Ci-Casz1 is localized in the posterior six cells; others are exclusively expressed in the posterior two presumptive pigment cells (Ci-Bzrap1, Ci-Slc45a2, Ci-Lrp4) (Fig. 6). Furthermore, in situ assays on embryos where FGF signalling was perturbed confirmed that the expression of some candidates (Ci-Casz1, Ci-Slc45a2 and Ci-Bzrap1) depends upon proper FGF signalling in PCPs (Fig. 7). Thus, our results reveal that each PCP lineage cell at the dorsal midline of the tailbud stage exhibits a unique combination of transcripts. These combinations constitute specific molecular signatures that define individual cell identities in the dorsal CNS.

Ci-Rab32/38 as ‘hub’ gene in Ciona pigment cell formation

To identify key effector genes necessary for Ciona PCPs formation and controlled by FGF signalling, we used a ‘reverse engineering’ approach. This method identifies co-regulated genes, that is, genes that are co-expressed across a set of experimental conditions. Co-regulated genes may be a part of the same pathway, regulated by the same transcription factor, and/or part of the same protein complex35, 36. The reverse engineering approach we used is based on computing correlations across pairs of genes from gene expression profiles35, 36. Specifically, we computed the Spearman correlation coefficient (SCC) values for each pair of probe sets across the 24 microarrays to build a network of co-regulated genes for Ciona PCP (Methods). The resulting gene network consists of 4,755,086 significant (P value <0.01) predicted connections among 30,969 probe sets. To identify groups of candidate FGF target genes in early PCPs, we built a smaller subnetwork (Fig. 8a) by selecting only the top 50 differentially expressed genes (when comparing FGFRDN to control samples at 8 hpf) together with their co-regulated genes from the main network (gene neighbours). In a co-regulatory network, most of the genes are connected to few other genes. Nevertheless, it is possible to identify a relatively small number of genes with a large number of neighbours, called hubs37, 38. In the subnetwork comprising early FGF targets, we identified a set of ‘hub genes’ and ranked them according to the number of neighbours (Supplementary Table 2). The top two hub genes were annotated as unknown protein in the Ciona genome, even though for the second one we found a blast homology with mouse C3a anaphylatoxin chemotactic receptor, part of the rhodopsin family of G protein-coupled receptor39. The remaining hub genes were the transcription factor Ci-NF001, the Ras-related protein Ci-Rab32/38, the Notch co-activator mastermind, Ci-mam and the Melanophilin, Ci-Mlph, a Rab effector involved in melanosome transport along microtubules (Supplementary Table 2).

Among the selected hub genes, we focused on the small GTPase Ci-Rab32/38. In vertebrate melanocytes, the two homologous proteins, Rab32 and Rab38, localize to secretory vesicles of the trans-Golgi network where they function redundantly to transport Tyrosinase- and Tyrp1-containing vesicles to the forming melanosomes40. The C. intestinalis genome encodes a unique ortholog, Ci-Rab32/38/rab-rp1/ltd (Ci-Rab32/38), showing a significant differential expression when comparing FGFRDN and control at 8 hpf, as confirmed by in situ assays on embryos electroporated with FGFRDN and Ets:Vp16 (Supplementary Fig. 6a–c).

Starting at the early neurula stage, Ci-Rab32/38 is specifically expressed in two cells that correspond to a10.97s and subsequently in a11.193s and a11.194s at tailbud stage as confirmed by the co-localization with Ci-tyrp1/2a (Fig. 8b,c).

The specific and early expression of Ci-Rab32/38 in the pigment cell lineage suggested a role in Ciona melanization by regulating vesicle dynamics during pigment cell formation. In keeping with this hypothesis, we observed that a group of genes also implicated in vesicular trafficking and melanosome biogenesis were differentially expressed when FGF signalling is blocked in PCPs between 8 and 12 hpf. This group includes the RabGAP Wdr67, the Rab effector MyRip/Exophilin and other genes involved in exocytosis and vesicle docking (RABEPK, Rab6-interacting kinesin-like protein, MPHOSPH1, ARF3, Doc2a/Rabphilin41, which we found expressed in a11.195s and a11.196s PCPs (Fig. 6g).

These data indicate that FGF signalling regulates Rab GTPases and associated cellular effectors controlling vesicular trafficking during Ciona pigment cell formation.

Role of Rab32/38 in Ciona pigment cell formation

In vertebrates, mutations in several Rab GTPases and their regulators result in pigmentation defects due to altered platelet and melanosome formation42. The chocolate (cht) mouse, characterized by a recessive Gly19 to Val mutation, displayed dramatically reduced numbers of melanosomes in adult RPE, in addition to a mild phenotype in skin melanocytes43, 44. Depletion of Rab32 in mildly hypopigmented cht melanocytes severely impairs pigmentation and melanosome morphology indicating that the hypopigmentation defects in cht melanocytes was due to a functional compensation of Rab32 (ref. 40). To investigate a possible function of the single Ci-Rab32/38 ortholog gene during Ciona pigmentation process, we first studied the effects of a dominant-negative Ci-Rab32/38 mutant (Rab32/38G19V). PCP-targeted expression of Rab32/38G19V in PCPs using the tyr enhancer14 altered both otolith and ocellus pigmentation. In more than 70% of larvae, pigmented sensory organs appeared either not or partially melanized (Fig. 9a). Considering that mouse Rab38G19V is unstable because it lacks posttranslational lipid modification that results in premature degradation40, we introduced another point mutation in Ci-Rab32/38G19V altering a second residue in the highly conserved GTPase fold domain (Rab32/38G19V+T23N) creating a dominant-negative mutated form (T23N)45. Targeted expression of Ci-Rab32/38G19V+T23N in pigment cell lineage caused a drastic reduction of the proportion of larvae (5%) showing two pigmented cells within the sensory vesicle, with more than 90% of perturbed larvae had alteration in at least one of the two pigmented cells (Fig. 9a). To further confirm the role of Ci-Rab32/38, we used short hairpin RNA (shRNA)-mediated RNA interference approach46 to knock down Ci-Rab32/38 activity and analyze the effect on pigmented cell formation. We used combinations of shRNA constructs targeting two sites in the Ci-Rab32/38 coding region. Double WMISH assays on embryos electroporated with either the U6>shRab32/38-C and -D or U6>shRab32/38-E and -F combination, showed that Ci-Rab32/38-targeting shRNAs markedly reduced its endogenous expression in PCPs (Fig. 9c,d). Expression of two pairs of constructs (U6>shRab32/38-C, -D and U6>shRab32/38-E, -F) resulted in up to 90% larvae showing defective pigment cell melanization (Fig. 9b). Taken together, these data indicate that the complete melanization of the pigment cell is thus significantly affected by Ci-Rab32/38G19V+T23N and Ci-Rab32/38 knockdown experiments, suggesting that proper Ci-Rab32/38 activity has an essential role in pigment cell terminal differentiation.


Discussion

The combination of cell sorting and transcriptome analyses identified pigment cell lineage candidate genes and provided new insights into the mechanisms of CNS patterning in Ciona. We have shown that FGF signalling functions as successive inputs rather than continuously during induction of PCPs. Sequential FGF inputs influence defined regulatory networks successively controlling early fate specification, the selection of definitive pigment cells within the PCP lineage and the processes that underlie pigment cell differentiation at later stages (summarized in Fig. 10).

Our data elucidate the roles of FGF signalling pathway during differential fate specification and anteroposterior patterning of the Ciona CNS. We show that, in the absence of FGF signalling, PCPs change their fate into anterior CNS cells, as judged by ectopic expression in the a9.49-derived PCPs of genes, which are normally expressed in 9.50s derivatives. According to these observations, we propose that FGF signalling imposes a pigment cell identity at the expense of the anterior CNS cell fate, distinguishing between anterior and posterior identities during the formation of Ciona nervous system. This is reminiscent of the situation in vertebrates, where FGF, together with Wnt signalling act to induce a posterior identity on neural tissue during neurulation47, 48.

In chick, FGF derived from the lens ectoderm is necessary to maintain the retina domain in the distal optic vesicle49. Moreover in FGF9 mouse mutants, the RPE expands into the presumptive retina suggesting that FGF9 is required to define the boundary between retina and RPE50. Our analyses suggest that in Ciona FGF signalling plays an initial role during the induction of pigmented versus anterior CNS cells. This could be related to the position of photoreceptor and pigmented cell precursors on the III row cells of the neural plate, both in close proximity to FGF signalling molecules, Ci-FGF8/17/18 and Ci-FGF9/16/20 (refs 11, 18) (Fig. 10). Although all PCPs have the potential to complete pigment cell differentiation, only the two posterior a11.113s will form the otolith and melanocyte of the ocellus. It is likely that specific transcription profiles impose their cell identity to uniquely defined PCP derivatives. Microarray data identified new candidates involved in Ciona pigment cell formation and showing a distinct gene expression profiles among PCPs (summarized in Fig. 10). Although some factors, such as Rab32/38 and Ci-Slc45a2, are associated with melanogenesis in vertebrates3, 51, to the best of our knowledge, the others (Ci-Mad/Mnt, Ci-Piwi, Ci-Lrp4, Ci-Bzrap1, Ci-Doc2a) have not been implicated in pigmentation before. Interestingly, Ci-tyr, Ci-tyrp1/2a and Ci-tyrp1/2b are expressed in PCPs long before melanization begins14 (Fig. 10). This precocious expression of tyr in the pigment cell lineage is a feature shared with the vertebrates15. However, it is not clear why only the two posterior PCPs, a11.193 cell pairs, are finally melanized although all the cells express melanogenic genes. It has been suggested that this could be related to the differential tyr activity in the presumptive pigment cells, where it is finally lost in the most anterior PCP cells that do not become pigmented52. Additional melanocyte-specific constituents influence melanosome biogenesis, such as the availability of the tyr substrate, L-tyrosine, as well as the accumulation of the L-DOPA cofactor that governs its enzymatic activity. Perhaps the permeability of PCP cell membrane controls the time and rate of access of the tyr substrate and cofactor to the enzyme. Here, we showed that the expression of a subset of transcripts encoding membrane transporters (that is, Ci-Slc2a6, Ci-Slc12a9, Ci-Slc7a6, Ci-Slc25a23) was decreased when FGF signalling was blocked in PCPs. In particular, Ci-Slc45a2, the homologue of SLC45A2, which is responsible for oculocutaneous albinism type 4 in humans51, is selectively expressed in the posterior a11.193 cell pairs. These data suggest that FGF signalling could regulate melanization in the two posterior PCP cells by increasing the permeability and regulating the ionic composition of these cells.

Here, we modelled a pigment cell lineage-specific transcriptional network to identify key effector genes involved in the early PCP induction and identified PCP-specific Ci-Rab32/38 as a hub gene. In mice, Rab32 and Rab38 control the trafficking of melanogenic enzymes (Tyr and Tyrp1) and are critical for melanosome maturation40. Moreover, it has been demonstrated that loss of Rab38 function alone causes a dramatic reduction in melanosome number within adult mouse RPE cells43. Also in Drosophila, mutation of Rab32/38 ortholog, Rab-RP1, causes eye pigmentation defects53.

These evidences parallel our finding in C. intestinalis where targeted perturbations of Ci-Rab32/38 activity interfere with a terminal melanization of the sensory organs. Considering the evolutionary relationship between the vertebrate small Rab GTPase members, Rab32 and Rab38, and the unique Ciona homologue, Ci-Rab32/38, ascidian pigment cells provide a simple and suitable system to identify the interacting partners for Rab32/38 and to study the functions of their regulators as well as the functional relationships with other factors during melanosome maturation.

We showed that FGF signalling is involved in the early pigment cell type-specific regulation of several components of melanosome vesicular trafficking such as Ci-Rab32/38 and the melanophilin, Ci-Mlph, which are differential expressed in all pairwise comparisons at 8 hpf. Further studies are required to investigate the role of FGF signalling in the regulation of Ci-Rab32/38 and other melanosome-specific factors, thus offering insights on the molecular events underlying pigment cell biology. Furthermore, our microarray data revealed that FGF regulates the expression of known target-associated albinism or reduced pigmentation, such as HPS1 and LYST (Hermansky–Pudlak and Chédiak–Higashi syndromes54, 55). Therefore, it appears that Ciona pigment cells are endowed with the basic molecular machinery that underlies vertebrate melanosome biogenesis and that Ciona melanosomes may thus follow a biogenesis process similar to that observed in more complex organisms. In summary, the importance of Ci-Rab32/38 function during pigment cell formation together with the identification of Ci-Rab32/38 and Ci-Mlph as central hub genes in the inferred 8 hpf network, strongly suggest that FGF signalling regulates the expression of key effectors required to prime the molecular machinery for melanization since the early FGF-mediated induction on PCPs.

Although detailed dissection of the different steps of Ciona melanosome maturation will require more in-depth analyses, Ciona pigmented cells provide an opportunity to analyse the conserved molecular basis of melanosome formation in vivo with high spatiotemporal resolution avoiding redundancy and compensations of vesicular trafficking components as well as the endocytic network that make difficult to dissect these cellular processes in a more complex system like vertebrates.

Note added in proof: During correction of proof we found this paper related to our work. Nicolas Haupaix, Philip B. Abitua, Cathy Sirour, Hitoyoshi Yasuo,Michael Levine and Clare Hudson. (2014). Ephrin-mediated restrictionof ERK1/2 activity delimits the number of pigment cells in the /Ciona/CNS. Developmental Biology. In Press. PII: S0012-1606(14)00346-7; DOI: 10.1016/j.ydbio.2014.07.010.


Methods
Animals and embryo electroporation

Adult C. intestinalis were collected from the Gulf of Naples and, for microarray experiments, obtained from M-REP (San Diego, CA, USA). Ripe oocytes and sperm were collected surgically and kept separately until in vitro fertilization. We used a chemical dechorionation to eliminate chorion and follicular cells surrounding and protecting the eggs. Dechorionated eggs were then in vitro fertilized from different individuals. Fertilized eggs were washed in Millipore-filtered sea water (MFSW) transferred in a solution containing 0.77 M Mannitol and with the amount of DNA used for electroporation that was typically 50–80 μg of each plasmid. The electroporation was carried out in Bio-Rad Gene Pulser 0.4 cm cuvettes, using Gene Pulser II (Bio-Rad)56, 57. Each experiment was performed at least three times. Embryos were staged according to the developmental timeline established in Hotta et al.19 Embryo imaging capture was performed with a Zeiss Axio Imager M1 and a Zeiss LSM 510 META confocal microscope.

Cell dissociation and FACS

Embryonic cell suspensions for FACS were obtained as follows: stage-selected embryos were rinsed three times in artificial calcium and magnesium-free sea water (CMF-ASW: 449 mM NaCl, 33 mM Na2SO4, 9 mM KCl, 2.15 mM NaHCO3, 10 mM Tris-Cl pH 8.2, 2.5 mM EGTA) at room temperature. Embryos were dissociated by thorough pipetting for 2–3 min in 1 ml CMF-ASW supplemented with 0.1% trypsin (w/v, Sigma, T- 4799). Trypsin activity was blocked by adding ice-cold CMF-ASW supplemented with 0.1% bovine serum albumin (BSA, Sigma, A-3311). Cell suspensions were then rinsed three times with ice-cold CMF-ASW+BSA (cell were pelleted at 4 °C, by centrifugation at 3,000 r.p.m. (~500–700 g), for 2 min). Cell suspensions were maintained on ice and filtered into 5 ml polystyrene round bottom tubes equipped with a 40 μm cell-strainer cap (Falcon) before cell sorting23. Fluorescent-FACSof ptyrp1/2a>2xGFP and ptyrp1/2a>GFP-labelled PCP founder lineage cells was performed using, positive selection of GFP and negative selection of a MyoD905-driven Venus YFP on the FACS Aria (BD Bioscienes) sorter (New York University, Center for Genomics & Systems Biology). GFP and YFP fluorescence distributions were gated on the population of cells. Special attention was devoted to exclude debris and cell aggregates by gating. GFP fluorescence was collected using a 500 nm Dichroic long pass filter and a 510 nm discriminating optical filter. YFP fluorescence was collected using 550 nm Dichroic long pass filter and a 550 nm discriminating optical filter. For every sort, a non-fluorescent cell suspension as well as GFP-only and YFP-only samples were analysed to determine the optimal settings for sorting. Spectral overlap was adjusted using these samples to compensate for the GFP fluorescence component detected by the YFP detector and for the YFP fluorescence component detected by the GFP detector. Collimated GFP and YFP fluorescence was separated using 500 nm dichroic long pass filter. GFP fluorescence (fluorescence lower than 500 nm) was reflected towards a photomultiplier detector (PMT) with a 510 nm discriminating filter in front of it. YFP fluorescence was detected from the transmitted 500 nm dichroic long pass filter using a 550 nm dichroic long pass optical filter. A 550 nm discriminating optical filter was placed in the light path towards the designated YFP fluorescence PMT detector23.

RNA Extraction and microarray hybridization

The RNAqueous-Micro Kit (Ambion) was used for RNA extraction. RNA quality was assayed on Agilent Bioanalyzer 2100 with RNA 6000 Pico Chips and the associated kit (Agilent). RNA amplification and probe synthesis were performed using The Ovation Pico WTA System following manufacturer‘s recommendations (NuGen Technologies). The fragmented and labelled single-stranded cDNA targets were generated using the Encore Biotin Module (NuGen Techlonogies) following manufacturer’s recommendations (#4200-12). Hybridization to the C. intestinalis custom-design GeneChip was conducted using standard Affymetrix procedures (Affymetrix GeneChip Expression Analysis Technical Manual (P/N 702232 Rev. 2) at the New York University, Center for Genomics & Systems Biology.

Microarray data analysis

Affymetrix DAT files were converted into CEL files using proprietary software of the Affymetrix scanner. CEL files were loaded into the R statistical environment58 (version 2.13.1) using the libraries Affy59 and makecdfenv from the Bioconductor60 collection of packages (version 2.8). The microarrays used in the present study represent an Affymetrix custom chip designed at the University of California, Berkeley and already used in gene expression studies of C. intestinalis23. Background correction, normalization and probe set summarization were performed using the Robust Multiarray Averaging method61 with default parameters producing a single normalized expression value for each probe set in each sample. We first applied independent filtering of the expression measures to remove all the probe sets presenting low expression levels in all the samples or generally no expression variations. This increases the efficiency of the analysis by mitigating the effects of P value correction during the selection of differentially expressed genes62. Cut-offs for the independent filters were defined for each time point to maximize the number of significantly differentially expressed genes and were executed by performing the following comparisons at each time point: FGFRDN versus control, Ets:Vp16 versus control and FGFRDN versus Ets:Vp16. Expression values were analysed for differential expression using linear models and the empirical Bayes method from the Limma package63 with default parameters. P values were adjusted using the Benjamini–Hochberg method64. We selected, as differentially expressed genes, only those genes associated to probe sets showing a corrected P value smaller than 0.1, which indicate a false discovery rate smaller than 10%.

Microarray probe annotations and GO enrichment analyses

Nucleotide sequences used to design the probes for each probe sets were obtained from the Affymetrix library file CINT06a520380F.sif. The next transcript models and annotations were downloaded from the ANISEED database28 on March 2011: (1) JGI version 1, (2) KYOTOGRAIL2005, (3) KH and (4) ENSEMBL. ANISEED provides, for each transcript model sequence as well functional (homologues, GO, and so on) annotations. Transcript models and annotations were parsed and joined together to form a unique collection. Probe set sequences were then aligned against the collection of transcript models using BLAST+ (ref. 65) with the megablast algorithm option. Results were filtered to select exclusive hits with coverage >75% and identity >90%. Based on the selected hits, we associated each probe set to the corresponding transcript models and annotations. To perform such analysis, we developed custom Perl scripts that make extensive use of the BioPerl toolkit66. Finally, a MySQL database containing all the information downloaded and generated was used to collect and manage the data. GO functional enrichment analyses were performed using the collected GO annotations associated to the probe sets. We compared the GO classifications of specific subsets of probe sets (representing significantly differentially expressed transcripts and transcripts connected to hub genes) to the GO classifications of all the probe sets annotated on the whole microarray. For each GO class, we used the Fisher exact text and the P values were corrected using the Benjamini–Hochberg method. A summary table with the collected annotations and the code used in the analyses are publicly accessible at http://bit.ly/nPepdH.

Reverse engineering of Ciona PCP co-regulatory network

We computed the SCC for each pair of probes across the 24 microarrays, obtaining a correlation matrix of dimension 30,969 × 30,969. We estimated the SCC significance for each pair of probes by computing t statistics of each SCC value and then using a Student’s t-test distribution to estimate the P value. To control the number of false positives owing to the multiple hypotheses test problem, we estimated the degrees of freedom of the t-test distribution from the data by fitting the parameters of a Student’s t-location scale distribution to the t statistics computed for all the probe pairs. We estimated the parameters by minimizing the squared error between the theoretical and the empirical distribution.

We also built a smaller early responsive to FGF signalling pathway subnetwork by starting from the top 50 differentially expressed transcripts (when comparing FGFRDN versus control at 8 hpf) and the transcripts connected to in the main network (gene neighbours). To identify ‘hub genes’ in this subnetwork, we ranked all the probe sets in the subnetwork according to their degree (that is, the number of connections a probe set has in the network).

In situ hybridization

The following ISH probes were obtained from plasmids contained in the C. intestinalis gene collection release I: tyrp1/2a (GC31h05), Bmp5/7 (GC46m03), Bzrap1(GC05h08), Doc2a (GC21m06), Lrp4 (GC42f11), Mad/Mnt (GC02d23), Rab32/38 (GC40p15), Six3/6 (GC11m13), FoxHa (GC32f03), Neurogenin (GC07n07), Pax3/7 (GC42e20), COE (GC02i14), Macho-1 (GC18g07), Msxb (GC42h24), Piwi (GC23f14) and full ORF Gateway-compatible clone: Casz1 (VES66_L10), RGR/Nut2 (VES94_E20), Pax6 (VES87_B20).

PCR amplification of transcription templates was with the following oligos: Slc45a2 fw (5′-TCCCAGTGGGTATGATGTTG-3′) and Slc45a2 rev (5′-TTATTGGTTGCAGCACGAAG-3′). DIG- and FLUO-labelled probes were transcribed with T7 RNA polymerase (Roche) using DIG and FLUO labelling mix (Roche) and purified with the RNeasy Mini Kit (Qiagen).

Single and double fluorescent ISH was performed as follows: Embryos were fixed at the different developmental stages for 2 h in 4% MEM-PFA (4% paraformaldehyde, 0.1 M MOPS pH 7.4, 0.5 M NaCl, 1 mM EGTA, 2 mM MgSO4, 0.05% Tween 20), then rinsed in cold phosphate-buffered saline (PBS) and stored in 75% ethanol at −20 °C. Before rehydration, embryos were treated with 2% H2O2 in methanol for 30 min at room temperature, and then rehydrated and manipulated as described through the post-hybridization washes. Embryos were then equilibrated in TNT (100 mM Tris pH 7.5, 150 mM NaCl, 0.1% Tween 20) and blocked in TNB (0.5% Roche blocking reagent in 100 mM Tris pH 7.5, 150 mM NaCl) for 1–2 h. Anti-FLUO-POD and anti-DIG-POD (Roche) antibodies were diluted 1:1,000 in TNB and incubated with embryos overnight at 4 °C. Embryos were washed extensively in TNT. For tyramide signal amplification, Cy3- and Cy5-coupled TSA Plus reagent (Perkin Elmer) was diluted 1:100 in amplification diluent and added to embryos for 5–10 min at room temperature. For 4,6-diamidino-2-phenylindole staining, DAPI (Sigma) was diluted 1:10,000 in PBT, and embryos were incubated overnight at 4 °C, then washed in PBT67.

Embryo manipulation

For antibody staining, embryos were fixed in 4% MEM-PFA for 30 min, rinsed several times in PBT (PBS with 0.1% Tween 20) and incubated with anti-dpERK antibody (Rabbit mAb, Cell Signalling technology, 1:100) and anti-GFP (Mouse mAb, Roche, 1:300) with 2% normal goat serum in PBT at 4 °C overnight. Embryos were washed in PBT and then incubated with donkey anti-mouse secondary antibody (1:500) coupled to Alexa Fluor 488 (Life Technologies) and anti-rabbit secondary antibody (1:500) coupled to Alexa Fluor 647 (Life Technologies) in PBT with 2% normal goat serum for 2 h at room temperature, then washed in PBT.

For direct GFP and mCherry visualization, embryos were fixed in MEM-FA (3.7% methanol-free formaldehyde, 0.1 M MOPS pH 7.4, 0.5 M NaCl, 2 mM MgSO4, 1 mM EGTA) for 30 min, rinsed several times in PBS with 0.15% Triton X-100, 0.05% Tween 20 for 1 h.

Drug treatments

Embryos were obtained by in vitro fertilization, then chemically dechorionated and grown in MFSW at room temperature, until the desired stage. Experimental embryos were treated with 4 μM MEK inhibitor compound U0126 (Sigma) in MFSW (U0126 stock solution: 2 mM in dimethylsulphoxide) (Sigma) for 30 min and then washed with MFSW to remove the compound. Treated embryos were allowed to develop to the desired stage and then fixed for whole-mount ISH or immunohistochemistry.

Construct preparation

For cell sorting, PCPs were GFP labelled using ptyrp1/2a>2Xgfp (for 8 hpf) and ptyrp1/2a>GFP (for 12 hpf); the ptyrp1/2a minimal enhancer region is the same described in Esposito and coauthors15. For ptyrp1/2a>GFP, the GFP coding sequence was digested with XhoI and SmaI restriction enzymes and fused downstream from ptyrp1/2a promoter, replacing mCherry in the ptyrp1/2a>mCherry vector, previously digested with XhoI-EcoRV to eliminate mCherry.

For ptyrp1/2a>2xGFP and ptyrp1/2a>LacZ, the tyrp1/2a enhancer was cloned upstream of the 2xGFP coding sequence and in the in the pBlueScript 1230 vector using NotI and EcoRI restriction sites.

For ptyr>Rab32/38, the cDNA of Rab38/32 was amplified by PCR using as template the cDNA clone (N. Satoh Gene Collection 1 ID: CiGC40p15), with the oligonucleotides Rab32/38_Not FW (5′-AAGGAAAAAAGCGGCCGCTTTAATGACCGAACAAAGAAAGGA-3′) and Rab32/38_EcoRI REV (5′-CCGGAATTCCGGTTATCTGCAACAACCGCTCTTTTC-3′). Ci-Rab32/38 cDNA, digested with NotI/EcoRI was inserted in place of mCherry in the construct ptyrp1a>mCherry, previously digested with NotI-EcoRI to eliminate mCherry.

Two different site-directed mutations in the GTPase domain of ptyr>Rab32/38 construct were induced to test Ci-Rab32/38 function. The first (ptyr>Rab32/38G19V) was prepared by a Gly to Val substitution in position 19 in the highly conserved phosphate/Mg2+ domain, where GTP binding pocket is located. The second point mutation was a threonine to asparagine substitution in position 23 (ptyr>Rab32/38G19V+T23N), that disrupts the Mg2+ binding site, resulting in greatly reduced affinity for guanine nucleotides. The point mutations were performed using (Quick Change Site-Directed Mutagenesis Kit from Stratagene/Agilent) using the following oligonucleotides:

Rab32/38G19VF (5′-gttcttgtgatcggcgagttagttgttgggaagacgagcgtcatt-3′),

Rab32/38G19VR (5′-aatgacgctcgtcttcccaacaactaactcgccgatcacaagaac-3′),

Rab32/38T23N F (5′-cgagttagttgttgggaagaatagcgtcattaaaagatat-3′),

Rab32/38T23N R (5′-atatcttttaatgacgctattcttcccaacaactaactcg-3′),

shRNA-mediated RNA interference knockdown of Rab32/38 function

The coding sequence of Ci- Rab32/38 was used to identify siRNA target sites with the DSIR68 and Public TRC Portal ( http://www.broadinstitute.org/rnai/public/seq/search) algorithms. Candidate 21-mer target sites were trimmed to 20 mer, and those that begin with A or G and end with A were further considered. Candidate 20 mer were BLASTed against the C. intestinalis genome on the Ghost database and selected to avoid polymorphic region as well as 20 mer with potential off-targets (that is, 16 nucleotide between positions 2 and 17 on the antisense). Target sense sequences were reverse complemented, and sense and antisense sequences were combined to mimic the hairpin structure corresponding to the microRNA miR-2213 (ref. 69). A vector containing the Ciona U6 promoter70 was used to clone complementary oligos encoding four shRNA constructs:

shRAB32/38-C-F 5′-agatGAAGAACGcgTTGGATAAttTGTAATTaGtAACATGTTATCCAAATCGTTCTTttt-3′

shRAB32/38-C-R 5′-aattaaaAAGAACGATTTGGATAACATGTTaCtAATTACAaaTTATCCAAcgCGTTCTTC-3′

shRAB32/38-D-F 5′-agatGAAGAACGcgTTGGATAAttTGTAATTaGtAACATGTTATCCAAATCGTTCTTttt-3′

shRAB32/38-D-R 5′-aattaaaAAGAACGATTTGGATAACATGTTaCtAATTACAaaTTATCCAAcgCGTTCTTC-3′

shRAB32/38-E-F 5′-agatGAAGAACGcgTTGGATAAttTGTAATTaGtAACATGTTATCCAAATCGTTCTTttt-3′

shRAB32/38-E-R 5′-aattaaaAAGAACGATTTGGATAACATGTTaCtAATTACAaaTTATCCAAcgCGTTCTTC-3′

shRAB32/38-F-F 5′-agatGAAGAACGcgTTGGATAAttTGTAATTaGtAACATGTTATCCAAATCGTTCTTttt-3′

shRAB32/38-F-R 5′-aattaaaAAGAACGATTTGGATAACATGTTaCtAATTACAaaTTATCCAAcgCGTTCTTC-3′

The data showed used 50 μg each of U6>shRab32/38-C, -D, -E, -F per electroporation.


Author contributions

C.R., L.A.C., F.R conceived and designed the experiments. C.R. and A.K.K. performed the experiments. C.R., A.K.K., R.S., L.A.C., F.R. analysed the data. R.S. did bioinformatical analysis of microarray data. G.G., D.d.B. did transcriptional networking analysis. F.R.-K., L.Z., L.A.C. contributed reagents/materials/analysis tools. C.R. L.A.C. F.R. wrote the paper.


Additional information

Accession codes: Microarray data are available in the ArrayExpress database ( www.ebi.ac.uk/arrayexpress) under accession number E-MTAB-2759.

How to cite this article: Racioppi, C. et al. Fibroblast growth factor signalling controls nervous system patterning and pigment cell formation in Ciona intestinalis. Nat. Commun. 5:4830 doi: 10.1038/ncomms5830 (2014).


Supplementary Material Supplementary Information

Supplementary Figures 1-6 and Supplementary Tables 1-2



We thank Maria Francone for her helpful technical assistance. We are grateful to Pui Leng Ip for support with the FACS and Farhana Salek for technical assistance in preparing constructs. C.R. and A.K.K. have been supported by SZN PhD fellowship. C.R. has been supported by a short-term fellowship from EMBO. This work is supported by the NIGMS/NIH (grant R01GM096032).


References
Quevedo W. C. Jr, & Fleischmann R. D.,Developmental biology of mammalian melanocytes. J. Invest. Dermatol.75, 116–120 (Year: 1980).6771335
Raposo G., & Marks M. S.,Melanosomes—dark organelles enlighten endosomal membrane transport. Nat. Rev. Mol. Cell Biol.8, 786–797 (Year: 2007).17878918
Ohbayashi N., & Fukuda M.,Role of Rab family GTPases and their effectors in melanosomal logistics. J. Biochem.151, 343–351 (Year: 2012).22323658
Yaar M., & Park H. Y.,Melanocytes: a window into the nervous system. J. Invest. Dermatol.132, 835–845 (Year: 2012).22158549
Fuhrmann S.,Eye morphogenesis and patterning of the optic vesicle. Curr. Top. Dev. Biol.93, 61–84 (Year: 2010).20959163
Hyer J.,, Mima T., & Mikawa T.,FGF1 patterns the optic vesicle by directing the placement of the neural retina domain. Development125, 869–877 (Year: 1998).9449669
Delsuc F.,, Brinkmann H.,, Chourrout D., & Philippe H.,Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature439, 965–968 (Year: 2006).16495997
Stolfi A., & Christiaen L.,Genetic and genomic toolbox of the chordate Ciona intestinalis. Genetics192, 55–66 (Year: 2012).22964837
Lemaire P.,Evolutionary crossroads in developmental biology: the tunicates. Development138, 2143–2152 (Year: 2011).21558365
Nishida H.,Cell lineage analysis in ascidian embryos by intracellular injection of a tracer enzyme. III. Up to the tissue restricted stage. Dev. Biol.121, 526–541 (Year: 1987).3582738
Meinertzhagen I. A.,, Lemaire P., & Okamura Y.,The neurobiology of the ascidian tadpole larva: recent developments in an ancient chordate. Annu. Rev. Neurosci.27, 453–485 (Year: 2004).15217340
Tsuda M.,, Sakurai D., & Goda M.,Direct evidence for the role of pigment cells in the brain of ascidian larvae by laser ablation. J. Exp. Biol.206, 1409–1417 (Year: 2003).12624175
Horie T.,, Orii H., & Nakagawa M.,Structure of ocellus photoreceptors in the ascidian Ciona intestinalis larva as revealed by an anti-arrestin antibody. J. Neurobiol.65, 241–250 (Year: 2005).16118796
Esposito R.,et al.New insights into the evolution of metazoan tyrosinase gene family. PLoS ONE7, e35731 (Year: 2012).22536431
Tief K.,, Hahne M.,, Schmidt A., & Beermann F.,Tyrosinase, the key enzyme in melanin synthesis, is expressed in murine brain. Eur. J. Biochem.241, 12–16 (Year: 1996).8898882
Abitua P. B.,, Wagner E.,, Navarrete I. A., & Levine M.,Identification of a rudimentary neural crest in a non-vertebrate chordate. Nature492, 104–107 (Year: 2012).23135395
Squarzoni P.,, Parveen F.,, Zanetti L.,, Ristoratore F., & Spagnuolo A.,FGF/MAPK/Ets signalling renders pigment cell precursors competent to respond to Wnt signal by directly controlling Ci-Tcf transcription. Development138, 1421–1432 (Year: 2011).21385767
Nishida H.,Induction of brain and sensory pigment cells in the ascidian embryo analyzed by experiments with isolated blastomeres. Development112, 389–395 (Year: 1991).
Hotta K.,et al.A web-based interactive developmental table for the ascidian Ciona intestinalis, including 3D real-image embryo reconstructions: I. From fertilized egg to hatching larva. Dev. Dyn.236, 1790–1805 (Year: 2007).17557317
Cole A. G., & Meinertzhagen I. A.,The central nervous system of the ascidian larva: mitotic history of cells forming the neural tube in late embryonic Ciona intestinalis. Dev. Biol.271, 239–262 (Year: 2004).15223332
Davidson B.,, Shi W.,, Beh J.,, Christiaen L., & Levine M.,FGF signalling delineates the cardiac progenitor field in the simple chordate, Ciona intestinalis. Genes Dev.20, 2728–2738 (Year: 2006).17015434
Wasylyk B.,, Hagman J., & Gutierrez-Hartmann A.,Ets transcription factors: nuclear effectors of the Ras-MAP-kinase signalling pathway. Trends Biochem. Sci.23, 213–216 (Year: 1998).9644975
Christiaen L.,et al.The transcription/migration interface in heart precursors of Ciona intestinalis. Science320, 1349–1352 (Year: 2008).18535245
Razy-Krajka F.,et al.Collier/OLF/EBF-dependent transcriptional dynamics control pharyngeal muscle specification from primed cardiopharyngeal progenitors. Dev. Cell29, 263–276 (Year: 2014).24794633
Imai K. S.,, Hino K.,, Yagi K.,, Satoh N., & Satou Y.,Gene expression profiles of transcription factors and signalling molecules in the ascidian embryo: towards a comprehensive understanding of gene networks. Development131, 4047–4058 (Year: 2004).15269171
Moret F.,et al.The dopamine-synthesizing cells in the swimming larva of the tunicate Ciona intestinalis are located only in the hypothalamus-related domain of the sensory vesicle. Eur. J. Neurosci.21, 3043–3055 (Year: 2005).15978015
Oliver G.,et al.Six3, a murine homologue of the sine oculis gene, demarcates the most anterior border of the developing neural plate and is expressed during eye development. Development121, 4045–4055 (Year: 1995).8575305
Tassy O.,et al.The ANISEED database: digital representation, formalization, and elucidation of a chordate developmental program. Genome Res.20, 1459–1468 (Year: 2010).20647237
Hudson C.,, Lotito S., & Yasuo H.,Sequential and combinatorial inputs from Nodal, Delta2/Notch and FGF/MEK/ERK signalling pathways establish a grid-like organisation of distinct cell identities in the ascidian neural plate. Development134, 3527–3537 (Year: 2007).17728350
Lanner F., & Rossant J.,The role of FGF/Erk signalling in pluripotent cells. Development137, 3351–3360 (Year: 2010).20876656
Pregizer S., & Mortlock D. P.,Control of BMP gene expression by long-range regulatory elements. Cytokine Growth Factor Rev.20, 509–515 (Year: 2009).19900834
Zilberberg A.,, Yaniv A., & Gazit A.,The low density lipoprotein receptor-1, LRP1, interacts with the human frizzled-1 (HFz1) and down-regulates the canonical Wnt signalling pathway. J. Biol. Chem.279, 17535–17542 (Year: 2004).14739301
Smith W. C., & Harland R. M.,Expression cloning of noggin, a new dorsalizing factor localized to the Spemann organizer in Xenopus embryos. Cell70, 829–840 (Year: 1992).1339313
Moser M.,et al.BMPER, a novel endothelial cell precursor-derived protein, antagonizes bone morphogenetic protein signalling and endothelial cell differentiation. Mol. Cell. Biol.23, 5664–5679 (Year: 2003).12897139
Belcastro V.,et al.Transcriptional gene network inference from a massive dataset elucidates transcriptome organization and gene function. Nucleic Acids Res.39, 8677–8688 (Year: 2011).21785136
Gambardella G.,et al.Differential network analysis for the identification of condition-specific pathway activity and regulation. Bioinformatics29, 1776–1785 (Year: 2013).23749957
Basso K.,et al.Reverse engineering of regulatory networks in human B cells. Nat. Genet.37, 382–390 (Year: 2005).15778709
Carro M. S.,et al.The transcriptional network for mesenchymal transformation of brain tumours. Nature463, 318–325 (Year: 2010).20032975
Nilsson G.,et al.C3a and C5a are chemotaxins for human mast cells and act through distinct receptors via a pertussis toxin-sensitive signal transduction pathway. J. Immunol.157, 1693–1698 (Year: 1996).8759757
Wasmeier C.,et al.Rab38 and Rab32 control post-Golgi trafficking of melanogenic enzymes. J. Cell Biol.175, 271–281 (Year: 2006).17043139
Hutagalung A. H., & Novick P. J.,Role of Rab GTPases in membrane traffic and cell physiology. Physiol. Rev.91, 119–149 (Year: 2011).21248164
Bultema J. J.,, Ambrosio A. L.,, Burek C. L., & Di Pietro S. M.,BLOC-2, AP-3, and AP-1 proteins function in concert with Rab38 and Rab32 proteins to mediate protein trafficking to lysosome-related organelles. J. Biol. Chem.287, 19550–19563 (Year: 2012).22511774
Lopes V. S.,, Wasmeier C.,, Seabra M. C., & Futter C. E.,Melanosome maturation defect in Rab38-deficient retinal pigment epithelium results in instability of immature melanosomes during transient melanogenesis. Mol. Biol. Cell.18, 3914–3927 (Year: 2007).17671165
Loftus S. K.,et al.Mutation of melanosome protein RAB38 in chocolate mice. Proc. Natl Acad. Sci. USA99, 4471–4476 (Year: 2002).11917121
Alto N. M.,Rab32 is an A-kinase anchoring protein and participates in mitochondrial dynamics. J. Cell. Biol.158, 659–668 (Year: 2002).12186851
Wang W.,, Razy-Krajka F.,, Siu E.,, Ketcham A., & Christiaen L.,NK4 antagonizes Tbx1/10 to promote cardiac versus pharyngeal muscle fate in the ascidian second heart field. PLoS Biol.11, e1001725 (Year: 2013).24311985
Dono R.,Fibroblast growth factors as regulators of central nervous system development and function. Am. J. Physiol. Regul. Integr. Comp. Physiol.284, R867–R881 (Year: 2003).12626354
Fletcher R. B.,, Baker J. C., & Harland R. M.,FGF8 spliceforms mediate early mesoderm and posterior neural tissue formation in Xenopus. Development133, 1703–1714 (Year: 2006).16554360
Dellinger R. P.,et al.Effects of inhaled nitric oxide in patients with acute respiratory distress syndrome: results of a randomized phase II trial. Inhaled Nitric Oxide in ARDS Study Group. Crit. Care Med.26, 15–23 (Year: 1998).9428538
Abashian A.,et al.Measurement of the CP violation parameter sin2 phi(1) in B(0)(d) meson decays. Phys. Rev. Lett.86, 2509–2514 (Year: 2001).11289969
Newton J. M.,et al.Mutations in the human orthologue of the mouse underwhite gene (uw) underlie a new form of oculocutaneous albinism, OCA4. Am. J. Hum. Genet.69, 981–988 (Year: 2001).11574907
Whittaker J. R.,Tyrosinase in the presumptive pigment cells of ascidian embryos: tyrosine accessibility may initiate melanin synthesis. Dev. Biol.30, 441–454 (Year: 1973).4634082
Ma J.,, Plesken H.,, Treisman J. E.,, Edelman-Novemsky I., & Ren M.,Lightoid and Claret: a rab GTPase and its putative guanine nucleotide exchange factor in biogenesis of Drosophila eye pigment granules. Proc. Natl Acad. Sci. USA101, 11652–11657 (Year: 2004).15289618
Wei M. L.,Hermansky-Pudlak syndrome: a disease of protein trafficking and organelle function. Pigment Cell Res.19, 19–42 (Year: 2006).16420244
Zhao H.,et al.On the analysis of the pathophysiology of Chediak-Higashi syndrome. Defects expressed by cultured melanocytes. Lab. Invest.71, 25–34 (Year: 1994).8041115
Christiaen L.,, Wagner E.,, Shi W., & Levine M.. Electroporation of transgenic DNAs in the sea squirt Ciona. Cold Spring Harb. Protoc doi:10.1101/pdb.prot5345 (Year: 2009).
Christiaen L.,, Wagner E.,, Shi W., & Levine M.. Isolation of sea squirt (Ciona) gametes, fertilization, dechorionation, and development. Cold Spring Harb. Protoc doi:10.1101/pdb.prot5344 (Year: 2009).
R Development Core Team. R: A language and environment for statistical computing R Foundation for Statistical Computing (Year: 2008).
Gautier L.,, Cope L.,, Bolstad B. M., & Irizarry R. A.,affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics.20, 307–315 (Year: 2004).14960456
Gentleman R. C.,et al.Bioconductor: open software development for computational biology and bioinformatics. Genome. Biol.5, R80 (Year: 2004).15461798
Irizarry R. A.,et al.Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics.4, 249–264 (Year: 2003).12925520
Bourgon R.,, Gentleman R., & Huber W.,Independent filtering increases detection power for high-throughput experiments. Proc. Natl Acad. Sci. USA107, 9546–9551 (Year: 2010).20460310
Smyth G. K.,Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol.3, Article3 (Year: 2004).16646809
Benjamini Y., & Hochberg Y.,Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B57, 289–300 (Year: 1995).
Camacho C.,et al.BLAST+: architecture and applications. BMC Bioinformatics10, 421 (Year: 2009).20003500
Stajich J. E.,et al.The Bioperl toolkit: perl modules for the life sciences. Genome Res.12, 1611–1618 (Year: 2002).12368254
Christiaen L.,, Wagner E.,, Shi W., & Levine M.. Whole-mount in situ hybridization on sea squirt (Ciona intestinalis) embryos. Cold Spring Harb. Protoc doi:10.1101/pdb.prot5348 (Year: 2009).
Auger H.,et al.Similar regulatory logic in Ciona intestinalis for two Wnt pathway modulators, ROR and SFRP-1/5. Dev. Biol.329, 364–373 (Year: 2009).19248777
Hendrix D.,, Levine M., & Shi W.,miRTRAP, a computational method for the systematic identification of miRNAs from high throughput sequencing data. Genome Biol.11, R39 (Year: 2010).20370911
Nishiyama A., & Fujiwara S.,RNA interference by expressing short hairpin RNA in the Ciona intestinalis embryo. Dev. Growth Differ.50, 521–529 (Year: 2008).18510713

Article Categories:
  • Article


Previous Document:  The Challenges of Success: Maintaining Access to High Quality PCI in the Face of Declining Procedura...
Next Document:  Same-Sex Domestic Violence: Prevalence, Unique Aspects, and Clinical Implications.