The Glide/Gcm fate determinant controls initiation of collective cell migration by regulating Frazzled
Abstract
Collective migration is a complex process that contributes to build precise tissue and organ architecture. Several molecules implicated in cell interactions also control collective migration, but their precise role and the finely tuned expression that orchestrates this complex developmental process are poorly understood. Here, we show that the timely and threshold expression of the Netrin receptor Frazzled triggers the initiation of glia migration in the developing Drosophila wing. Frazzled expression is induced by the transcription factor Glide/Gcm in a dose-dependent manner. Thus, the glial determinant also regulates the efficiency of collective migration. NetrinB but not NetrinA serves as a chemoattractant and Unc5 contributes as a repellant Netrin receptor for glia migration. Our model includes strict spatial localization of a ligand, a cell autonomously acting receptor and a fate determinant that act coordinately to direct glia toward their final destination.
https://doi.org/10.7554/eLife.15983.001Introduction
Neurons and glia show mutual reliance in many functional and developmental aspects of biology. Glia migrate collectively and over long distances to establish an intricate relationship with neurons. Defective glia migration is associated with several human diseases including glial brain tumors and defective regeneration following injury in the nervous system (Klämbt, 2009; Kocsis and Waxman, 2007; Oudega and Xu, 2006). Hence, a thorough understanding of the molecules involved in the process of glia migration may contribute to the development of therapeutics for these pathologies. Research progress in recent years has revealed the involvement of chemotropic cues in glia migration (von Hilchen et al., 2010; Chen et al., 2010; Spassky et al., 2002; Kinrade et al., 2001; Liu et al., 2012). Netrins, a class of secreted laminin-related extracellular proteins, have been described as chemotropic guidance cues for axons and migrating cells during neural development (Lai Wing Sun et al., 2011; Harris et al., 1996; Ishii et al., 1992; Kennedy et al., 1994; Mitchell et al., 1996; Serafini et al., 1994). In vertebrates, Netrins secreted by the floor plate cells act differentially on the migrating oligodendrocyte precursor cells and this differential outcome is dependant upon the type of receptors expressed in the migrating cells (Spassky et al., 2002; Wolf and Friedl, 2009; Jarjour, 2003; Tsai et al., 2003). Netrins act as chemoattractants through the DCC/Frazzled family of receptors (von Hilchen et al., 2010; Lai Wing Sun et al., 2011; Chan et al., 1996; Keino-Masu et al., 1996; Kolodziej et al., 1996; Timofeev et al., 2012) and chemorepellants through the Unc5 receptor family (Keleman and Dickson, 2001; Labrador et al., 2005). In Drosophila, the two Netrins (NetA and NetB [Harris et al., 1996; Mitchell et al., 1996; Keleman and Dickson, 2001]) and their receptor Frazzled (Fra) mediate the attraction of embryonic longitudinal glia toward the midline (von Hilchen et al., 2010). Despite the extensive knowledge on these ligands and receptors, the transcriptional control underlying chemoattraction and the impact of this pathway remain largely unknown. For example, what is the specific role of the receptors and how do they contribute to the different steps of collective migration, such as initiation, maintenance and arrest? Related to this issue, which transcription factors regulate the cell-specific and timely expression of the receptors? We here investigate the chemoattraction cascade that controls cell migration using the chain of glial cells across the L1 nerve in the developing Drosophila wing (Aigouy et al., 2004, 2008; Berzsenyi et al., 2011; Kumar et al., 2015).
We show that only one of the two Drosophila Netrins, NetB, serves as a chemoattractant for collective glia migration. The role of Fra is to control the time of initiation of glia migration in a dosage-dependent manner, whereas Unc5 acts as the repellant receptor that controls glial arrest. Finally, we identify the transcription factor that controls the expression of Fra at the appropriate time and levels: Glial cell deficient/Glial cell missing (Glide/Gcm or Gcm, for the sake of simplicity), the fate determinant that is expressed early and transiently in the glial lineages (Hosoya et al., 1995; Vincent et al., 1996; Jones et al., 1995). Thus, we find that an early gene, which regulates the expression of transcription factors that execute a specific cell fate, also regulates effector genes that controls late developmental events. To our knowledge, this is the first report showing that a fate determinant directly controls collective cell migration, prompting us to revisit the role and mode of action of these types of molecules during development.
Results
Frazzled expression in the glia of the developing Drosophila wing
Fly wings are innervated by two major sensory nerves that navigate along the so-called L1 vein located at the anterior margin (L1 nerve) and along the L3 vein (L3 nerve) (Figure 1a–d). Glial cells originating from the sensory organ precursors (SOPs) present on the anterior margin migrate proximally, i.e. toward the central nervous system (CNS) following the axon bundle and ensheathing it throughout its length. L1 glia initiate migration at around 18 hr After Puparium Formation (hAPF), reach the level of the Costal nerve at around 22–24 hAPF and join the glial cells on the Radius by 28 hAPF. The migratory process has been accordingly subdivided into three steps: ‘Initiation’, ‘Costa reach’ and ‘Complete migration’ (Figure 1a–c).
To gain insight into the molecular pathway that triggers collective glia migration, we first examined the expression of the Fra chemoattractant receptor by using the pan glial lines repo-Gal4 UAS-PH-GFP (henceforth repo>GFP) or gcm-Gal4 UAS-CD8-GFP (gcm>GFP), which label the glial membranes. Fra is detected in glia at the time these cells begin to move as well as in the underlying axons (Figure 1e–g’’, Figure 1—figure supplement 1a–c, g). The Fra protein seems to be evenly distributed along the L1 glial chain (Figure 1—figure supplement 1d–g). We confirmed these data using the CoinFLP technique (Bosch et al., 2015) to generate WT and fra knock down clones (KD), obtained by means of the UAS-fra-RNAi line, within the same wing (Figure 1h, i). The WT clones covering glia along the L1 nerve and surrounding cells show expression of Fra (GFP-expressing cells in Figure 1k–k’’’), whereas Fra levels are considerably reduced in the glia and in the surrounding cells within the KD clones (RFP expressing cells in Figure 1j–j’’’). Glia can be identified by the expression of the pan glial marker Repo (WT glia are Repo/GFP positive, whereas fra KD glia are Repo/RFP positive).

Expression of Fra in wing glia.
(a–c) Immunolabeled repo>GFP wing (glia in green: anti-GFP; neurons in red: anti-22c10) at different stages. (a) Initiation of migration, (b) reaching the level of the costa, and (c) migration completion. (d) Schematic of a 18 hAPF developing Drosophila wing, insets indicate the regions shown in panels (e–f’’’’’) and panels (i–k’’’). L1 and L3 indicate L1 and L3 nerves. (e–e’’’’’) gcm>GFP/+ 18 hAPF wing immunolabeled with anti-22c10 (red), anti-Fra (gray) and anti-GFP (green). mCD8-GFP was used to label the membrane. (f–g’’) 18 hAPF repo>GFP/+ wing, immunolabeled with anti-22c10 (neurons in red), anti-Fra (gray) and anti-GFP (glia in green). (e–g’’) The presence of Fra in the glial soma (white arrowheads) at the front of migration. The position of the high-magnification panels (g–g”) is highlighted by the white rectangle in (f). Maximum confocal projections are shown in all figures, unless otherwise specified. White arrowheads indicate the glial cells that are expressing Fra. (h) Schematic representations of the coinFLP technique (modified from Bosch et al. (2015); and the phenotypes of the different cells. (i–k’’’) Immunolabeled Fra KD/WT-coinFLP wing at 18 hAPF. The WT clones display GFP labeling at the membranes (anti-GFP), the fra KD clones display RFP labeling at the membranes (anti-RFP); glial nuclei are labelled with anti-Repo in red and anti-Fra is in gray. (i) A projection of 10 confocal sections from a 18 hAPF wing. The arrowheads indicate glial cells and the white lines outline the L1 nerve. (j–k’’’) Individual sections: (j–j’’’) represents section 2/10; (k–k’’’) represents section 5/10, which corresponds to a deeper layer than section 2/10. (j, k) The overlay of the three channels (anti-RFP/Repo, anti-GFP and anti-Fra), (j’, k’) show anti-GFP alone, (j’’,k’’) anti-Fra and (j’’’, k’’’) the overlay of anti-RFP/Repo and anti-Fra. Glial cells are indicated by white arrowheads, the dashed lines indicate the fra KD clones and the continuous lines indicate the WT clones. For technical reasons, RFP (membrane labeling) and Repo (nuclear labeling) are shown in the same channel. Note the decrease in Fra levels in the fra KD clones. The scale bar in (a–c) represents 80 µm, in (e–f) 10 µm and in (i–k) 5 µm.
Thus, migrating glial cells of the peripheral nervous system (PNS) express Fra.
Next, we studied the impact of Fra on glia migration using the semiquantitative approach described by Kumar et al.(2015). In short, we assessed the percentage of wings displaying complete glia migration at 28 hAPF (migratory index: MI), as an estimation of migration efficiency (Kumar et al., 2015). For each genotype, at least 30 wings were analyzed. We first focused on the most characterized loss-of-function (LOF) allele fra3 (Kolodziej et al., 1996). As this mutation is embryonic lethal in homozygous conditions, we analyzed fra3 heterozygous wings and found incomplete L1 glia migration in a significant fraction of samples, as shown by the position of the glial nuclei (Figure 2a–c). The number of glial cells is not affected and hence cannot be the cause of the migratory defect (Figure 2d). We reasoned that nuclei may not migrate properly, but that glial processes may still reach the final destination. To test this possibility, we assessed the migratory index of glial cells labelled by the >GFP transgene, which allowed us to visualize glial processes in flies that were heterozygous for the .. These wings also show incomplete glia migration as shown by the position of glial GFP labeling (Figure 2e–g) suggesting that glial cells require Fra to complete their migration.

Role of Fra in wing glia.
(a, b) 28 hAPF wings labeled with anti-Repo (glial nuclei): (a) WT wing showing complete migration;(b) fra3/+ wing showing incomplete migration (dashed box). (c) Histogram representing the migratory index (MI) of the indicated genotypes, calculated using nuclear labeling (anti-Repo). The MI indicates the percentage of wings displaying complete migration (i.e. in which the glial chain reaches the proximally located glia on the Radius nerve) and was assessed at 28 hAPF unless otherwise specified. (d) Histogram representing the number of glial nuclei in the indicated genotypes. (e) Histogram representing the MI of the indicated genotypes. The MI was calculated using the membrane GFP transgenic line (UAS-mCD8-GFP). (f, g) 28 hAPF wings labeled with anti-GFP (glial processes): (f) repo>GFP wing showing complete migration; (g) repo>GFP/fra3/+ wing showing incomplete migration (dashed box). In this and in the following figures, stars indicate p values: ***p<0.0001; **p<0.001; *p<0.05. Bars indicate the s.e.m. In this and in the following graphs on fixed wings, n≥30. Scale bars: 80 μm.
-
Figure 2—source data 1
Migratory index and repo count of of fra3 wings in WT background.
This file contains underlying raw data for Figure 2.
- https://doi.org/10.7554/eLife.15983.005
In summary, the Fra receptor is expressed in glial cells and is seemingly necessary for their efficient migration.
fra plays an instructive role in L1 glia migration
The lethality of fra3homozygous mutation and the expression of Fra in glia as well as in neurons prompted us to assess the role of glial fra expression in migration specifically. The knock down of fra using the gcm>GFP driver, which is the earliest glial driver, reveals a significant decrease in migration efficiency as compared to that observed in the control wings (Figure 3a, compare blank and light blue columns). To exclude the possibility of off-target effects, we analyzed wings that express the UAS-fra-RNAi together with the UAS-fra transgene, and found complete rescue of the migratory phenotype induced by the fra KD (Figure 3a, patterned light and dark blue column). This strongly suggests that the RNAi line induces a specific phenotype and that fra acts in a cell autonomous manner. Finally, as a complementary approach, we reintroduced fra expression only in the glial cells of fra3/+ animals (Figure 3—figure supplement 1a), and found that this rescues the migratory phenotype, albeit partially. Given the high levels of Fra expression in the gcm>fra GOF wings (Figure 3g–g’’’’), it is unlikely that this partial rescue is due to suboptimal levels of Fra. Rather, the lack of total rescue may be ascribed to the neuronal requirement of Fra. Figure 1e–f’’’’’ shows that Fra is expressed in glia as well as in neurons, and indeed around one third of fra3/+ wings show an axonal navigation phenotype (Figure 3—figure supplement 1b–d). In these wings, axonal navigation is delayed, which may indirectly affect glial migration. To further check the role of Fra neuronal expression in glia migration, we used a driver that is expressed specifically in the L1 neurons: nsyb-Gal4 (West et al., 2015; Pauli et al., 2008; Riabinina et al., 2015). nsyb-driven expression starts to be detected around 18 hAPF; but the expression is more prominent at 22 hAPF and onwards (Figure 3—figure supplement 1f–h”). We then knocked down fra specifically in neurons using nsyb-Gal4 and found no delayed glia migration (Figure 3—figure supplement 1i), except in one third of the wings, which present an axonal navigation defect (not shown). This is in line with the above-mentioned fra3/+ data and with previous data showing that axon navigation defects indirectly affect glia migration (Giangrande et al., 1993).

An instructive role of the chemoattractant receptor Fra in collective glia migration.
(a–b) Histogram representing the MI upon fra knock down (fra KD) or overexpression (fra GOF) using the gcm>GFP/+ line. In (b), the MI was calculated at 24 hAPF. The MI was calculated using the membrane GFP transgenic line. c–c’”’, e–e’”’ and g–g’”’ Expression profiles of Fra in gcm>GFP/+, gcm>fra KD and gcm>fra GOF animals at 24 hAPF. For the sake of consistency, in all genotypes, we show the glial cells that are at the front of the chain. Note the reduced protein levels in gcm>fra KD (e–e’”’) and enhanced levels in gcm>fra GOF animals (g–g’”’) as compared to those found in gcm>GFP/+ animals (c–c’”’). (d, f and h) Wing schematics, boxes indicate the regions shown in the above, high-magnification panels. (i) Graphical representation of the migratory behavior of gcm>GFP/+, gcm>fra KD and gcm>fra GOF wings during three highlighted phases: initiation, costa reach and complete migration (n=10). (j, k) Snapshots from a 21:48 hAPF time-lapse analysis of gcm>GFP/+ and gcm>fra GOF wings. This corresponds to the time by which most control L1 glia (gcm>GFP/+) have reached the level of the Costa (white line) (j), whereas L1 glia overexpressing Gcm (gcm>fra GOF) have already completed migration by that time (white arrow) (k). The two panels show representative samples. (l–o”) Expression profiles of Fra in gcm>GFP/+ and gcm>fra GOF animals at 15 hAPF. See the enhanced protein levels in gcm>fra GOF animals (n–o”) as compared to those found in gcm>GFP/+ animals (l–m”). The position of the high-magnification panels (m–m”) is highlighted by the dashed white rectangle in (l), whereas that of panels (o–o”) is highlighted in (n). Please note that (l–o”) are comprised of few a sections rather than maximum confocal projections. Scale bars: (c–g’”’), (j–o”), 10 μm.
-
Figure 3—source data 1
Migratory index and time-lapse analysis of fra conditional mutants in the gcm>GFP/+ background.
This file contains underlying raw data for Figure 3.
- https://doi.org/10.7554/eLife.15983.007
We then asked whether Fra has an instructive role in glia migration and assessed whether migration is more efficient upon overexpressing fra in glia using the gcm-Gal4 driver (in fra GOF lines). We first checked the MI of fra GOF wings at 28 hAPF and found that the percentage of wings that shows complete migration is higher as compared to that of control wings (Figure 3a, compare white and dark blue columns). Since most control wings show complete migration by 28 hAPF (90%), we also analyzed an earlier stage, when migration has been achieved in only a few control wings (24 hAPF; 12.5%). We found that many more fra-overexpressing wings show complete migration (68.5%, Figure 3b). This strongly suggests that high doses of Fra significantly increase the efficiency of glia migration. Fra levels in fra KD and fra GOF conditions are indeed reduced and increased, respectively, as compared to those observed in control wings (Figure 3c–h).
To clarify why the migratory efficiency decreases in fra KD animals, we performed time-lapse analysis and found that reducing the levels of Fra affects the first step of migration (initiation), as glia start to migrate later than the control chain (Figure 3i, compare the initiation step of white and light blue columns and Figure 3—figure supplement 1e). Accordingly, fra-overexpressing glia start migrating earlier than control glia, indicating that the phenotype is due to precocious initiation (Figure 3i–k, compare the initiation step of white and dark blue columns), and that this is associated with precocious Fra accumulation in fra GOF glial cells compared to that in control glia (Figure 3l–o”). Finally, the live imaging data demonstrates that the speed of migration at the time of initiation is not higher in fra GOF glia than in control glia, showing that the observed phenotype is due to precocious initiation rather than to acceleration (Figure 3—figure supplement 2a).
The cytoplasmic tail of Fra is known to play a major role in mediating Fra-dependent attractive responses in vivo and in cell culture studies (Bashaw and Goodman, 1999; Hong et al., 1999; Ming et al., 1997). We therefore asked whether this region is important in mediating glia migration. A transgenic construct that lacks the Fra intracellular cytoplasmic domain was previously described as a dominant negative mutation (Garbe et al., 2007) and indeed the expression of this reporter significantly reduces migration efficiency (Figure 3a, compare white and green-gray columns).
In summary, we find that acts in a cell autonomous manner in glial cells to mediate migration. Furthermore, our data show that that levels of Fra are critical for the initiation stage of glial cell migration and that this is mediated through the cytoplasmic domain of Fra.
The efficiency of glia migration depends on the dose of Gcm
Interestingly, the migratory phenotype of fra3/+ glia that also carry the gcm>GFP driver is much stronger than that of fra3/+ glia (compare Figure 2e patterned black an gray column, MI = 67% and Figure 4a patterned light gray and white column, MI = 14%). The phenotype is further enhanced in glia that express both the gcm>GFP and a gcm RNAi line (Figure 4a patterned brown and pink column). Since the gcm-Gal4 driver is a hypomorphic gcm allele that results from the insertion of a Gal4-containing transposon into the gcm promoter (Jacques et al., 2009), the above result raised the possibility that Gcm and Fra interact genetically. To explore this possibility, we analyzed the glia migration phenotype in double heterozygous conditions for fra and two other gcm hypomorphic alleles, including the gcmrA87 enhancer trap carrying a LacZ transposon inserted into the gcm promoter and the imprecise excision line gcm (Riabinina et al., 2015; Vincent et al., 1996; Jacques et al., 2009). In addition, we also used a gcm null mutation, gcmN7-4 (Bernardoni et al., 1997). This confirms that reducing the dose of Gcm enhances the fra3-mediated phenotype (Figure 3—figure supplement 2b). Furthermore, we crossed the fra3 mutation with a gcm driver that does not affect the gcm locus, a transgenic line carrying 6Kb of the gcm promoter fused to the Gal4 gene, which is inserted on the third chromosome (Flici et al., 2014). In these wings, we did not observe the enhanced migratory phenotype present in the fra3 , gcm-Gal4 wings (Figure 3—figure supplement 2c). These findings suggest that the two genes act in the same genetic pathway that impinges on glial cell migration.

Gcm affects collective glia migration in a dose-dependent manner.
(a) MI in fra mutants, gcm LOF and rescues of the indicated genotypes: gcm>GFP/+, gcm/fra LOF, gcm>gcm LOF, gcm>gcm KD and gcm>gcm GOF wings. (b) MI calculated at 24 hAPF in gcm>gcm GOF wings. The MI was calculated using the membrane GFP transgenic line. (c) Graphical representation of the migratory behavior of gcm>GFP/+, gcm>gcm LOF and gcm>gcm GOF wings at three migratory phases. (d–f) Snapshots of a 22:48 hAPF time-lapse analysis on gcm>GFP/+, gcm>gcm GOF and gcm>gcm LOF wings. (d) 22:48 hAPF corresponds to the time at which L1 glia have surpassed the level of the Costa (white line) in the control gcm>GFP/+ wing, which is why this time point was chosen to compare the position of L1 glia in the different genetic backgrounds. By this time, (e) glia migration is already complete in gcm>gcm GOF wing (white arrow), while in (f) gcm>gcm LOF wing, L1 glia are still at the level of Costa (white line). Scale bar: (e–g), 10 μm.
-
Figure 4—source data 1
Summary of the role of gcm in glia migration.
This file contains underlying raw data for Figure 4.
- https://doi.org/10.7554/eLife.15983.012
Gcm is a transiently expressed transcription factor that acts very early in glial differentiation (Hosoya et al., 1995; Vincent et al., 1996; Jones et al., 1995). In situ hybridization on wild-type (WT) wings has shown that the gcm RNA becomes detectable by 8–9 hAPF (Van de Bor and Giangrande, unpublished data) and fades in glial cells by 24 hAPF (Popkova et al., 2012). To clarify the role of Gcm on glia migration, we analyzed wings that are only mutant for gcm and used hypomorphic alleles that allow bypassing the lethality of the null mutation. Three allelic conditions were tested: gcm-Gal4/gcmrA87(gcm LOF) (Figure 4a, light brown column, MI = 63%, vs. the MI = 90% of the control line gcm>GFP/+ (white column)), gcm-Gal4 homozygous (Figure 4a, orange column, MI = 41%) and gcmrA87/gcmN7-4 transheterozygous (Figure 3—figure supplement 2b, last column, MI = 28%) animals. Migration is indeed less efficient when the amount of Gcm is reduced and the MI is restored to normal values upon reintroduction of Gcm expression (Figure 4a, patterned pink columns). Finally, we used the UAS-gcm-RNAi line to reduce the amount of Gcm (gcm KD) and also observed a migratory defect (Figure 4a, dark brown column, MI = 25%). The rescue obtained upon co-expressing the UAS-gcm and the UAS-gcm-RNAi transgenes indicates that gcm plays a regulatory role in migration and that the RNAi effects are specific (Figure 4a, patterned pink-dark brown column). Furthermore, overexpressing Gcm using the UAS-gcm transgene (gcm GOF) is sufficient to increase the migration efficiency of glial cells, as the percentage of wings showing complete migration increases compared to that in control animals (Figure 4a, pink column and Figure 4b).
To determine which migratory step is affected by Gcm, we analyzed the gcm LOF (gcm>GFP/gcmrA87) and the gcm GOF (gcm>GFP/+;UAS-gcm/+) wings by confocal time-lapse microscopy and found that migration starts later in Gcm LOF and earlier in gcm GOF wings compared to that observed in control wings (Figure 4c–f). Thus, Gcm affects the initiation of glial cell migration, like Fra, and it does so in a dose-dependent manner: gcm GOF enhances the efficiency of this step and gcm LOF lowers it.
In summary, the levels of Gcm and Fra are crucial for the initiation of glia migration.
Gcm affects migration independently of its role as a fate determinant
The glial to neuron conversion described in gcm mutant flies prompted us to ask whether this defect could impact the glial migratory process indirectly (Hosoya et al., 1995; Vincent et al., 1996; Jones et al., 1995; Bernardoni et al., 1997). We inspected the rate of glia to neuron conversion in gcm KD wings by using the anti-Elav antibody, which specifically recognizes neurons at the analyzed stages (Figure 5—figure supplement 1a–b’’ gcm KD). Only a minor fraction of the wings contains converted cells and only a few cells are converted into neurons, strongly suggesting that fate conversion is not the cause of the altered MI (on average, 14% of the gcm KD wings show up to 10 Repo/Elav double-positive cells along the L1 nerve, >15 wings were analyzed per genotype; none were observed in wild-type wings). In addition, we restrained our analysis to gcm KD wings that do not show Repo/Elav positive cells and still found a strong delay in migration (MI = 45%).
These data strongly suggest that the glial migration phenotype observed in gcm mutant animals is not due to its early requirement in glial cell determination, but specific to cell migration.
Gcm affects migration independently of glial cell number
Another cause for the migratory phenotype observed in the gcm mutant wings might be the control exerted by Gcm on the total number of glial cells. It is indeed plausible that the number of cells in the collective somehow affects the mechanical forces that control migration efficiency, for example through the amount of chemoattractant receptor. The number of Repo-positive cells in gcm LOF and KD backgrounds is indeed lower than that in the wild-type glial chain and, accordingly, gcm GOF wings contain supernumerary glial cells (Figure 5a, light brown, dark brown and pink columns, respectively).

The effects of Gcm on collective glia migration are independent of the number of Repo-positive nuclei.
(a) Bar chart showing the numbers of L1 glial nuclei in the various genotypes. (b) MI in the various genotypes, calculated using the membrane GFP transgenic line.
-
Figure 5—source data 1
Effects of gcm on glia migration are independent of any change in the number of repo+ cells.
This file contains underlying raw data for Figure 5.
- https://doi.org/10.7554/eLife.15983.014
To assess the impact of glial cell number in migration efficiency more directly, we analyzed wings that overexpress proteins promoting or repressing cell division. The exit from the cell cycle results from the timely inactivation of the Cyclin-dependent protein kinase (Cdk) and Cyclin E (Cyc E) complexes. On the one hand, String/Cdc25 encodes a phosphatase that triggers mitosis by activating the Cdc2 kinase, hence enabling cell proliferation (Edgar et al., 1994; Edgar and O'Farrell, 1989; Lasko, 2013). On the other hand, Dacapo functions as an inhibitor of the Cdk–Cyc E complex both in vivo and in vitro, ultimately leading to cell cycle arrest (de Nooij et al., 1996; Lane et al., 1996). First, we produced animals overexpressing String or Dacapo in glia (gcm>GFP) and verified that this induces a significant change in glial cell number compared to that present in control wings (Figure 5a, compare white column with columns with vertical and horizontal lines). Then, we analyzed the migration efficiency in both backgrounds and found that it is not affected (Figure 5b), even though the glial number increase induced by String overexpression and the decrease induced by Dacapo overexpression are comparable to those observed in gcm GOF and LOF, respectively. Finally, we found defective glia migration even in gcm KD wings containing a wild-type number of glial cells, in agreement with the above data (Figure 5—figure supplement 1c).
Thus, the absolute number of glia does not affect migration efficiency, further corroborating the hypothesis that Gcm specifically affects this process.
gcm expression in blood cells does not affect glial migration
Another round of evidence supports the hypothesis that gcm expression in glia, as opposed to effects originating in other cells, is absolutely necessary for their migration . It is already known that gcm is also expressed in hemocytes (Bernardoni et al., 1997; Bataillé et al., 2005; Waltzer et al., 2010). So, we checked whether the severe migratory delay is specific to reducing the amount of gcm in glia or in hemocytes. To clarify this, we used an independent hemocyte driver, collagen-Gal4, which is expressed in the embryonic hemocytes (Asha et al., 2003). We knocked downgcm in hemocytes (collagen-Gal4 crossed with the UAS-gcm-RNAi) and found no defect in glia migration (Figure 5—figure supplement 1d). Furthermore, as Gcm is only present embryonically in blood cells,we used a gcm-Gal4, tubulin-Gal80thermosensitive(ts) line to specifically knock down fra in glia after embryogenesis and confirmed that, in these conditions too, glia migration is affected (Figure 5—figure supplement 1e).
Taken together, these data suggest that Fra and Gcm act cell autonomously in glial cells to control migration on the developing wing disc.
fra is a direct Gcm target
The gene-expression profile of Fra and the observed genetic interaction between Fra and Gcm prompted us to ask whether Gcm acts on glia migration by inducing Fra expression directly. A DNA adenine methyltranferase identification (DAM ID) screen aimed at identifying the direct targets of Gcm does indeed suggest that this potent transcription factor may directly control the expression of genes involved in glia migration, including fra (Cattenoz et al., 2016). There are three canonical Gcm-binding sites (GBS) in the fra locus, two of which are located at the position of a strong DAM ID peak, which is indicative of Gcm binding (Figure 6a). To confirm this data, we assessed whether fra expression is activated by Gcm in Drosophila S2 cells. To do so, we built a GFP reporter under the control of the fra locus containing two GBSs (Figure 6a–a’). qRT PCR assays clearly show an increase in the GFP levels upon co-transfection of the reporter vector with the Gcm expression vector and this effect is dose dependent (Figure 6b, columns with a red color gradient). To further demonstrate that the effect of Gcm on fra is direct (Figure 6a’’), we showed that Gcm-dependent activation of the reporter is completely abolished upon mutagenesis of the two GBSs (Figure 6b, columns with a yellow color gradient). The levels of transfected Gcm were verified by qRT PCR and those of GFP were confirmed by Western blot assays (Figure 6—figure supplement 1a–c).

fra is a direct Gcm target.
(a) Schematic representation of the fra gene indicated in blue. Lines indicate the introns, thin regions indicate the untranslated sequences, thick regions indicate the coding exons, pale blue arrowheads indicate the direction of transcription. GBSs are indicated in red, the black histogram above shows the gcm DAM ID peak at the indicated genomic coordinates. (a’) Represents the two GBSs that were amplified and put in front of the GFP reporter to generate a fra WT plasmid. (a”) The two GBSs were mutated in order to generate a fra Mut plasmid. (b) qRT-PCR analysis of the GFP levels using a WT (red) or a mutant reporter plasmid (yellow) upon co-transfection with increasing doses of the Gcm expression vector (reported as μg) (n=3). (c–e”’) Immunolabeling of gcm>GFP/+, gcm>gcm KD and gcm>gcm GOF wings at 28 hAPF using anti-22c10 (red), anti-Fra (gray) and anti-GFP (green). Note the reduced Fra levels in the gcm>gcm KD and the increased levels in the gcm>gcm GOF wings as compared to those in the gcm>GFP/+ wings (arrows). (f) Wing schematics: the blue rectangle indicates the region shown in (c–c”’) and (e–e”’), the black one, the region shown in (d–d”’) (g) MI of the indicated genotypes calculated using the membrane GFP transgenic line. (h–j) 28 hAPF wings from the indicated genotypes: control (gcm>GFP/+), gcm>gcm KD and rescue (gcm>gcm KD/fra GOF). (h) gcm>GFP/+ wing shows complete migration. (i) gcm>gcm KD wing shows incomplete migration (white bracket). Note that this is a composite image. (j) gcm>gcm KD/fra GOF wing showing complete migration. Scale bar: (c–e”), 10 μm; (h, j), 80 μm.
Finally, we complemented the in vitro data with two in vivo assays. First, we showed that Fra levels are affected in opposite direction in gcm LOF and GOF wings (Figure 6c–f). In gcm GOF wings, the levels of the Fra protein already increase by 15 hAPF, in agreement with the observed precocious initiation of migration (Figure 6—figure supplement 1d–f, compare with Figure 3l–o”). Second, we hypothesized that Fra may constitute an important target of Gcm in L1 glia migration and showed that overexpressing Fra in gcm KD wings is sufficient to completely reverse the migratory phenotype that results from gcm KD (Figure 6g, patterned dark blue and brown columns and Figure 6h–j).
In summary, these data suggest that the role of the transcription factor in glial cell migration is through the direct activation of the of chemoattractant receptor, Fra; this in turn implies that an early fate determinant is also capable of directly controlling late developmental events through inducing the expression of cell-surface molecules in a dose dependent manner..
Role and expression of Netrins
The Fra receptor is known to signal through the two Netrin ligands to mediate cell signaling. To understand which of these ligands were involved in the molecular pathway controlling glial cell migration, we assayed migration in animals that were null mutants for the two netrins: NetAΔ or NetBΔ (Harris et al., 1996; Mitchell et al., 1996; Brankatschk and Dickson, 2006). In this experiments, we find that while NetAΔ mutant glia do not display a defective migratory index (Figure 7a, compare gray and blue columns), migration is significantly affected in NetBΔ mutant glia (Figure 7a, compare gray and red columns).

NetB serves as a chemoattractant in collective glia migration.
(a) MI of the indicated genotypes. Histogram shows the MIs quantified for WT, for NetAΔ or NetBΔ single mutant wings and for NetAΔNetBTM or NetAΔNetBmyc wings. (b) NP4151-Gal4 driven GFP expression of NetB in a 19 hAPF wing. Proximal NetB expression as revealed by the GFP labeling (green). Anti-22C10 is in red. (c, d) Histograms representing the MI of the indicated genotypes. (e) NP4012-Gal4 driven GFP expression of NetA in a 20 hAPF wing. NetA is expressed in the wing epithelium as revealed by the profile of GFP (green). Anti-22C10 is in red. (f, g) Histograms representing the MI of the indicated genotypes. For this figure the MI was calculated by nuclear labeling. Scale bars: (b, e), 80 μm.
-
Figure 7—source data 1
Summary of the role of Netrins in glia migration.
This file contains underlying raw data for Figure 7.
- https://doi.org/10.7554/eLife.15983.019
As it has long been known that Netrins can elicit short-range attraction at the Drosophila embryonic midline (Harris et al., 1996; Mitchell et al., 1996; Brankatschk and Dickson, 2006), we next checked whether NetB could act in a similar manner in the migrating L1 glia. For this purpose, we used a transgenic line that does not express NetA and that only expresses the wild type or the membrane-tethered form of NetB, which is incapable of signal transduction (Brankatschk and Dickson, 2006). The lines NetAΔ NetBTM and NetAΔ NetB had been obtained through homologous recombination and hence express the modified or the wild-type NetB protein at near endogenous levels. Glia migration was found to be comparable to that of wild-type glia in NetAΔ NetB/Y wings, whereas it is delayed in NetAΔ NetBTM/Y wings. This suggests that NetB must be secreted for proper glia migration to occur (Figure 7a, compare light and dark purple columns) implying that the NetB ligand functions as a long-range guiding cue for the receptor Fra that is expressed in glial cells.
We next sought to determine the source of the ligand NetB. To do this, we used a transgenic line that is routinely employed as a reporter of NetB expression (Timofeev et al., 2012; Hayashi et al., 2002). NP4151>UAS GFP reports NetB expression in the proximal region of the wing (19 hAPF) (Figure 7b). This profile of expression fits well with the distal to proximal migration of L1 glia, and we reasoned that if NetB were to act as a chemoattractant, its loss should cause glial migratory defects similar to those induced by the loss of Fra. We therefore knocked down NetB by crossing NP4151-Gal4 to UAS-NetB-RNAi (NetB KD) flies and uncovered a severe migration defect, which could be rescued by simultaneously expressing UAS-NetB and UAS-NetB-RNAi (Figure 7c, first three columns). Moreover, overexpressing NetB in its territory of expression enhances the efficiency of glial cell migration, as revealed by the MI at 24 hAPF (Figure 7d). By contrast, ectopic expression of NetB in the posterior wing compartment using engrailed-Gal4 driver (en>) (Hidalgo, 1994; Lawrence and Morata, 1976) or in the distal part of the wing using GMR 29F05-Gal4 (GMR 29F05>)(Pfeiffer et al., 2008) does not affect the migration efficiency of glial cells (Figure 7—figure supplement 1a).
Finally, NetA is almost ubiquitously expressed in the epithelium, as revealed by the use of the Gal4 transgenic line NP4012 (Figure 7e) (Timofeev et al., 2012; Hayashi et al., 2002) and knocking down or overexpressing NetA with that driver has no impact on glia migration at early or at late stages (Figure 7f,g). Moreover, NetA overexpression in the NetB expression territory fails to enhance migration efficiency (Figure 7—figure supplement 1b) or to rescue the NetB KD phenotype confirming our finding that NetA is not involved in the process of glial cell migration (Figure 7—figure supplement 1c).
Altogether, our data strongly support the hypothesis that secreted NetB in the proximal wing provides a crucial chemoattractant cue for Gcm-mediated Fra expression, hence controlling the efficiency of glia migration. Our data also suggest that NetB may contribute but is not sufficient to trigger directionality.
Unc5 controls the late phase of L1 glia migration
unc5 is a repellant receptor for Netrins and has been previously shown to be transiently expressed and required in the embryonic exit and peripheral glia (PG) associated with both the segmental and intersegmental nerves (von Hilchen et al., 2010; Keleman and Dickson, 2001; Freeman et al., 2003). We therefore asked whether unc5 might also be involved in the migration of glia in the developing wing. To do this, we first analyzed its spatio-temporal expression profile during development. Unc5 is undetectable in the wing disc at 15 hAPF (Figure 8—figure supplement 1a–a”). Its expression is first seen at 18 hAPF at low levels and this expression progressively decreases to completely fade away by 29 hAPF (Figure 8a–d’’’).
If unc5 were to act as a repulsive receptor, the efficiency of glia migration would increase if its expression decreases, but neither RNAi-mediated KD of unc5, nor the null unc58 mutation (Labrador et al., 2005) affect glia migration efficiency at the early and late stages of development (Figure 8e, first three columns; Figure 8—figure supplement 1b). Thus, the loss of Unc5 does not seem to enhance the migration efficiency of L1 glia in the developing wing. We then asked whether Unc5 expression must be tightly regulated, and found that overexpressing Unc5 affects the efficiency of glia migration by delaying it, a phenotype that is opposite to the Fra overexpression phenotype (Figure 8e, compare blank and green columns). This delayed migration phenotype was rescued by reducing the levels of unc5 using the unc5 KD construct or the mutant allele unc58, in this overexpression background. This rescue demonstrates a direct effect of unc5 on glia migration (Figure 8e, patterned green and purple column, Figure 8—figure supplement 1c, third column).

Unc5 may act as a repellant in glia migration.
(a–d”’) Wing immunolabeled with anti-22c10 (red), anti-Unc5 (gray) and anti-GFP (green) in the transgenic line gcm>GFP/+ at different migratory stages. Unc5 starts being expressed in glia and neurons at around 18 hAPF and fades away by 29 hAPF. Arrowheads show glia that express Unc5 at weak levels (a–c’’’) or that do not express Unc5 (d–d’’’). Note that (a–d”) are comprised of a few sections rather than maximum confocal projections. (e–g) Histograms representing the MI of the indicated genotypes, which in this figure was calculated using the membrane GFP transgenic line. Scale bars: (a–d”), 10 μm.
-
Figure 8—source data 1
Summary of role of unc5 in glia migration.
This file contains underlying raw data for Figure 8.
- https://doi.org/10.7554/eLife.15983.023
Thus, fra and unc5 serve opposite roles in glia migration, with fra being necessary to trigger migration and unc5 delaying the migratory process. The two molecules seem to work in the same signaling pathway because the unc5 GOF phenotype is further enhanced by lowering the levels of Fra (Figure 8e, patterned green and light blue column). Also, the migratory phenotype induced by unc5 overexpression is rescued by simultaneously overexpressing Fra (the rescue was analyzed at an early stage for a better quantification; Figure 8f, patterned green and dark blue column). Finally, knocking down unc5 rescues the fra KD phenotype (Figure 8—figure supplement 1c, last two columns).
Altogether, the data described above strongly suggest that Unc5 can act as a repellant but that its expression is not sufficient to affect migration efficiency, which is mostly controlled by Net–Fra interaction. A likely explanation for these results is that Unc5 contributes at late migratory stages, and time-lapse movies indeed show that unc5 overexpression delays migration but does not affect the initiation step (Figure 9—figure supplement 1a). To further test this hypothesis, we assessed the effects of knocking down Unc5 using repo-Gal4, which is a late glial driver compared to gcm-Gal4. In these conditions, we did see an acceleration in migration completion (Figure 8g). Finally, since Unc5 has been identified in microarrays obtained upon overexpressing Gcm (Freeman et al. 2003) and in the DAM ID screen performed to identify the direct targets of Gcm (Cattenoz et al., 2016), we asked whether Unc5 expression is directly induced by Gcm. Co-transfection assays in S2 cells similar to those performed on fra and using an Unc5–GFP reporter containing two GBS close to the unc5 transcription start site, suggest that Unc5 is a rather weak target of Gcm (Figure 8—figure supplement 1d). This was further confirmed in gcm gain-of-function in vivo assays (data not shown) similar to those mentioned above on fra.

Repo regulates Fra at late stages.
(a) Histogram showing the endogenous expression of fra upon S2 cell transfection with a Repo expression vector. The y-axis represents the relative expression levels in cells transfected with Repo compared to that observed in cells not transfected with the Repo expression vector (n=3). (b) Schematic summarizing the regulatory network. (c, d) Graphs representing the MI found at early and late stages in fra GOF animals crossed with the repo-Gal4 line. In this figure, the MI was calculated using the membrane GFP transgenic line.
-
Figure 9—source data 1
Summary of role of repo in glia migration.
This file contains underlying raw data for Figure 9.
- https://doi.org/10.7554/eLife.15983.026
In summary, these data suggest that apart from a role for chemoattraction mediated through Fra and NetB, which determines the initiation of glial cell migration, a later chemorepulsion mediated through Fra-Unc5 ensures the termination of migration. Furthermore, the expression of Unc5 is also directly regulated by the transcription factor gcm suggesting that the same early determinant controls the timely expression of effector genes that initiate and terminate glial cell migration.
Discussion
Collective migration is a complex biological process that allows cells to leave their place of birth and reach their destination in a coordinated and timely manner. Here we dissect the signaling pathway underlying the collective migration of glia on the developing Drosophila wing. The careful coordination between extrinsic cues and intrinsic, cell autonomous mechanisms ensure that timely migration of the glia occurs on the wing disc. The chemoattractant receptor Fra controls glia migration in response to a long-distance signal, thechemoattractant NetB. Fra is expressed before glia start to move and triggers migration initiation in a dose-dependent manner. Such tight control of Fra expression depends upon the Gcm transcription factor, which directly induces the expression of Fra at threshold levels. Thus, the glial determinant factor also affects a specific step of collective migration by regulating the key effector gene.
Collective migration comes in different flavors: streams, chains, sheets and clusters, which all imply tight coordination and cell–cell interactions (Klämbt, 2009; Gilmour et al., 2002; Gupta and Giangrande, 2014; Marín et al., 2010; Rørth, 2003). The small cluster of Drosophila border cells migrates through nurse cells towards the oocyte in response to growth factors (Montell, 2003; Rørth, 2009). The stream of hundreds of proliferating cells of the fish lateral line migrate directionally upon expressing different chemokine receptors within the collective (Dambly-Chaudiere et al., 2007; David et al., 2002; Ghysen and Dambly-Chaudière, 2004; Haas and Gilmour, 2006). The current challenge is to dissect the role of the signaling pathways in the different steps (initiation, maintenance and arrest) and features of collective migration (adhesion, overall velocity and timing) and to analyze this in vivo. The second challenge is to understand how those signaling pathways are regulated, an issue that has been addressed at the post-transcriptional level (Yu and Bargmann, 2001) but much less so at the transcriptional level.
By focusing on the Netrin signaling pathway, which has been extensively studied in the context of axonal navigation and cell proliferation, we have addressed the above questions using the migrating wing glia of Drosophila. The analysis of the mutant phenotypes and the time-lapse approach show that the Fra receptor controls glia migration along the L1 nerve. Reducing the amount of Fra delays migration, whereas excessive Fra in the glial cells triggers their precocious migration. Thus, the Fra receptor plays an instructive role in the first step of collective migration, initiation. These data highlight the importance of quantitative regulation: large cohorts of cells are likely to need strong forces to switch from an immotile to a motile phenotype; therefore only the strong expression of the receptor allows migration toward the chemoattractant. Other migratory collectives also depend on quantitative regulation. Typically, epithelial cells migrating in groups exert much stronger forces than an individual cell before and after epithelial-mesenchymal transition (du Roure et al., 2005). The cell adhesion molecule E-cadherin plays an important role in supporting the directional and efficient collective migration of epithelial cells by mediating adhesion and force generation (Li et al., 2012). Additionally, the collectively migrating border cell cluster cannot exert enough adhesive/pulling force to move between the nurse cells in the absence of E-cadherin (Niewiadomska et al., 1999).
The timely regulation of the Fra receptor and its early role in collective migration allowed us to show that the Gcm fate determinant affects migration through this receptor. Thus, early genes not only trigger the expression of transcription factors that in turn implement a specific developmental program, but also directly contribute to the acquisition of specific phenotypes, such as the migratory potential. Similarly, the Lim-homeodomain transcription factor Islet was shown to specify the electrical properties of motor neurons by repressing the expression of the ion channel Shaker during development, suggesting that the regulation of late genes by early transcription factors might also be a common phenomenon (Wolfram et al., 2012). These observations prompt us to revisit the role of the so-called ‘master regulators’, which indeed control processes beyond fate choices.
Netrins are thought to act either in a gradient at long range, as secreted molecules, or at short range, as membrane-tethered molecules (Lai Wing Sun et al., 2011; Brankatschk and Dickson, 2006). In contrast with studies in the Drosophila embryo and visual system (von Hilchen et al., 2010; Timofeev et al., 2012; Brankatschk and Dickson, 2006), we observed that secreted NetB in the developing wing acts at a long range: glia migration is delayed when solely membrane-tethered NetB is available at near-endogenous levels. While we still do not know whether the distinction between long- and short-range signaling depends on the complement of surface receptors and the associated transduction pathways, it seems that different strategies are used in specific processes. Future studies will reveal whether long-range signaling is specifically dedicated to migration over many cell diameters and/or to large collectives.
NetA and NetB have been suggested to act on the embryonic longitudinal glia through Fra (von Hilchen et al., 2010) and both Netrins are reported to act as chemoattractants, even though in one case, only NetB was proposed to mediate dendritic targeting via Fra (Matthews and Grueber, 2011). Clearly, downregulating or overexpressing NetA has no impact on the migrating glial cells of the wing and NetA cannot rescue the NetB phenotype. While it is possible that NetA may affect other aspects of wing glia biology, our data indicate that NetB alone serves as a chemoattractant ligand to guide L1 glia migration upon signaling to Fra. Studies in other systems will assess whether it is only NetB that works as a long-range chemoattractant. If this were the case, it would be interesting to assess whether the different behavior of NetA and NetB relies on the intrinsic potentials of the two ligands (e.g. the affinity with which NetA and NetB bind to their receptors) or on extrinsic cues (e.g. cell-specific cofactors modulating NetA and/or NetB activity).
The chemoattractant receptor DCC/Fra and the chemorepellant receptor Unc5 have been mostly studied in different cell types, so it is unclear whether the same process and cells require the counteracting activities of these molecules. While mammalian migrating oligodendrocyte precursor cells express both DCC and Unc5, the precise function of these receptors has not been assessed (Tsai et al., 2003). We found that the two receptors Fra and Unc5 are expressed by the same set of PNS glial cells, with Fra appearing earlier than Unc5. The two receptors seem to have a different impact on wing glia migration and their stoichiometry is an important factor, as reducing Unc5 levels does have an effect on migration when fra is downregulated (Figure 8—figure supplement 1b last columns). The fact that overexpressing Unc5 and Fra have opposite effects suggests that Unc5 converts the Netrin-mediated signal from attraction to repulsion in the Drosophila wing glia, in agreement with published data on chimeric molecules. For example, when fused to the cytoplasmic domain of Unc5, the extracellular domain of DCC is as effective in mediating the chemorepellant response to Netrin as a wild type Unc5 (Hong et al., 1999). Additionally, a study performed on Drosophila showed that fusing the cytoplasmic domain of Fra with the extracellular domain of Unc5 signals chemoattraction (Keleman and Dickson, 2001). This may explain how the same cues such as Netrins can work differently on various aspects of collective migratory behavior.
The key element in this study is time: Gcm accumulates and triggers the onset of glial migration through fra but subsequently fades away. The fact that the expression of fra stays on till the end of migration suggests that another player, possibly a direct target of Gcm, maintains its expression. Being a major Gcm target, repo represents a potential candidate. Our in vivo and cell transfection data using a Repo expression vector are in line with this hypothesis and suggest that Repo may play a regulatory role in the maintenance of fra at late migratory stages (Figure 9a). It is known that Repo expression first depends on Gcm but subsequently becomes independent of it through autoregulation (Flici et al., 2014), indicating a transition from early to late events in gliogenesis (see model in Figure 9b). We propose that the initial Fra expression in wing glia strictly depends on Gcm and triggers the initiation of migration (see the graph in Figure 3—figure supplement 1e). Accordingly, the overexpression of Fra using the repo promoter does not induce migratory defects (Figure 9c,d) and repo>fra KD wings show only modest defects in migration initiation (Figure 9—figure supplement 1b). Later on, however, when Gcm fades away, Repo takes over in controlling Fra expression, hence affecting later phases of migration (Figure 9—figure supplement 1b, see model in Figure 9b). Conversely, the unc5 repellant is a weak target of Gcm thatmay require Repo and controls migration completion.
With respect to what has previously been described in embryonic glia migration, we propose a different role for Netrins and their receptors in the migration of L1 glia. The migration of embryonic longitudinal glial cells seems to rely on very early expression of Fra, which is detected in the longitudinal glial precursor, but no longer in the longitudinal glial cells (von Hilchen et al., 2010; Kolodziej et al., 1996). In our work, we do detect Fra expression in the wing glia. However, while in the case of the longitudinal glia both Netrins guide the glia toward the midline, the L1 glia seem to require only NetB. Moreover, ectopic expression of Netrin affects the migration of the longitudinal glia, whereas it does not do so in wing glia, hence suggesting that Netrin has a role in migration efficiency rather than on directionality. And finally, we propose that Fra and Unc5 act in different migratory phases of the same glial population, the first phase being under the direct control of the Gcm glial determinant. Overall, this paper displays a model in which the ligand NetB is expressed early during wing development (Figure 7—figure supplement 2), whereas its receptor Fra is expressed right before migration onset and cell autonomously decides the time of migration initiation of the glial chain. Unc5 appears later and contributes to migration termination. The acquisition of the Fra levels that trigger migration is part of the cell-specification program dictated by the Gcm glial determinant, which not only induces the expression of downstream transcription factors but also directly implements a specific developmental program by triggering the expression of effector genes.
Together, these findings in the Drosophila wing suggest that the dynamic and coordinated action of fate determinants and chemotropic cues contribute to the timely and efficient migration of the glial cells. A similar molecular mechanism, relying on Netrins or other localized attractive cues and their receptors, may be used in other cases of collective migration.
Materials and methods
Fly stocks and genetics
Request a detailed protocolFly stocks were raised at 25°C in standard medium. repo-Gal4 (indicated as repo>) was used to detect glial-specific expression of UAS-PH-GFP (fusion protein between the Pleckstrin homology domain of PLC-d and the GFP coding sequence). gcm-Gal4 UAS-mCD8-GFP (membrane localization) (gcm>GFP) was used as an early glial-specific driver (Jacques et al., 2009). gcm-Gal4,tubGal80ts (Soustelle et al., 2007). The other strains used were gcm 6Kb> 42; fra3 (von Hilchen et al., 2010); UAS-fra (von Hilchen et al., 2010); UAS-fra-RNAi (Manhire-Heath et al., 2013); UAS-fraΔC (Garbe et al., 2007); gcmrA87 (Vincent et al., 1996); gcmN7-4 (Vincent et al., 1996); UAS-gcm (F18A) (Bernardoni et al., 1998); UAS-gcm-RNAi; gcm>GFP/gcm>GFP (used as a homozygous mutant of gcm) (Popkova et al., 2012); UAS-string (Inaba et al., 2011); UAS-dacapo (Lane et al., 1996); gcm34 (Vincent et al., 1996); NetAΔ (von Hilchen et al., 2010; Newquist et al., 2013a, 2013b); NetBΔ (von Hilchen et al., 2010; Newquist et al., 2013a, 2013b); NetAΔ NetBTM (Brankatschk and Dickson, 2006); NetAΔ NetBmycBrankatschk and Dickson, 2006) (note that both NetAΔ NetBTM and NetAΔ NetBmyc encode the c-myc epitope tags); UAS-NetB-RNAi (Manhire-Heath et al., 2013); UAS-NetB (Timofeev et al., 2012); NP4151-Gal4 and NP4012-Gal4 (DGRC, Kyoto) (Timofeev et al., 2012; Hayashi et al., 2002); UAS-NetA-RNAi (Manhire-Heath et al., 2013); UAS-NetA (Newquist et al., 2013a, 2013b); unc58 (Labrador et al., 2005); collagen-Gal4 (Asha et al., 2003); engrailed-Gal4 driver (Hidalgo, 1994; Lawrence and Morata, 1976); GMR 29F05-Gal4 (Pfeiffer et al., 2008); UAS-unc5-RNAi; UAS-unc5 (von Hilchen et al., 2010); nsyb-Gal4>GFP-LAMP (B# 42714); elav-Gal4. The RNAi lines were obtained from Bloomington and/or the VDRC stock center.
To generate the coinFLP clones, the hsFLP;UAS-fra-RNAi built with hsFLP (B# 6) and UAS-fra-RNAi (Manhire-Heath et al., 2013), was crossed with UAS-mCD8-RFP,LexAop-mCD8-GFP;CoinFLP-LexA.GAL4 (B# 58754) (Bosch et al., 2015). Heat shock was carried out in wandering third instar larvae at 37°C for 30 min, wings from female pupae were dissected at 18 hAPF.
The summary table here below provides the genotypes for each figure reporting MI, for the sake of simplicity.
Genotypes | Abbreviated as |
---|---|
Figure 2 | |
WT | Control |
fra3/+ | LOF |
repo>GFP/+ | Control |
repo>GFP/fra3 | fra LOF |
Figure 3 | |
gcm>GFP/+ | Control |
gcm>GFP/+;UAS fra RNAi/+ | fra KD |
gcm>GFP/+;UAS fra RNAi/UAS fra | fra KD, Rescue |
gcm>GFP/+;UAS fra/+ | fra GOF |
gcm>GFP/+UAS fraΔC | |
Figure 4 | |
gcm>GFP/+ | Control |
gcm>GFP/fra3 | fra LOF |
gcm>GFP/fra3;UAS gcm RNAi/+ | |
gcm>GFP/gcmrA87 | gcm LOF |
gcm>GFP/gcmrA87;UAS gcm/+ | gcm LOF, Rescue |
gcm>GFP/gcm>GFP | |
gcm>GFP/gcm>GFP;UAS gcm/+ | Rescue |
gcm>GFP/+;UAS gcm RNAi/+ | gcm KD |
gcm>GFP/+;UAS gcm RNAi/UAS gcm | gcm KD, Rescue |
gcm>GFP/+;UAS gcm/+ | gcm GOF |
Figure 5 | |
gcm>GFP/+ | Control |
gcm>GFP/gcmrA87 | gcm LOF |
gcm>GFP/+;UAS gcm RNAi/+ | gcm KD |
gcm>GFP/+;UAS gcm/+ | gcm GOF |
gcm>GFP/UAS string | string GOF |
gcm>GFP/UAS dacapo | dacapo GOF |
Figure 6 | |
gcm>GFP/+ | Control |
gcm>GFP/+;UAS gcm RNAi/+ | gcm KD |
gcm>GFP/+;UAS gcm RNAi/UAS fra | Rescue |
Figure 7 | |
WT | Control |
NetAΔ/Y | |
NetBΔ/Y | |
NetAΔNetBTM/Y | |
NetAΔNetB/Y | |
NetB>/Y | Control |
NetB>/Y;+/+;UAS NetB RNAi/+ | NetB KD |
NetB>/Y;+/+;UAS NetB RNAi/UAS NetB | NetB KD, Rescue |
NetB>/Y;+/+;UAS NetB/+ | NetB GOF |
NetA>/Y | Control |
NetA>/Y;+/+;UAS NetA RNAi/+ | NetA KD |
NetA>/Y;UAS NetA/+ | NetA GOF |
Figure 8 | |
gcm>GFP/+ | Control |
gcm>GFP/+;UAS unc5 RNAi/+ | unc5 KD |
gcm>GFP/unc58 | unc5 LOF |
gcm>GFP/+;UAS unc5/+ | unc5 GOF |
gcm>GFP/+;UAS unc5; UAS unc5 RNAi | unc5 GOF, Rescue |
gcm>GFP/+;UAS unc5/UAS fra RNAi | |
gcm>GFP/+;UAS unc5/UAS fra | |
gcm>GFP/+;UAS fra/+ | fra GOF |
repo>GFP/+ | Control |
repo>GFP/UAS unc5 RNAi | unc5 KD |
Figure 9 | |
repo>GFP/+ | Control |
repo>GFP/+;UAS fra/+ | fra GOF |
Molecular cloning
Request a detailed protocolFor the fra gene, oligonucleotides (see primers sequence below) surrounding the GBSs were designed with flanking restriction sites for KpnI at the 5’ extremity and NHeI at the 3’ extremity. Each pair of oligonucleotides was used to amplify the genomic region encompassing the GBSs using the Expand High fidelity polymerase (Roche). The amplicons were digested with 20 U of KpnI (NEB # R3142S) and 20 U of NheI (NEB # R3131S) in Cutsmart buffer (NEB # B7204S) for 2 hr min at 37°C. The digested amplicons were then cleaned using the PCR clean-up kit (Macherey-Nagel (MN) # 740609) according to manufacturer’s instructions.
For ligation, 50 ng of the digested probe were used and cloned into the pGreen Pelican vector overnight at 18°C. 1 μL of the ligated product was used for transformation of electro competent DH5α bacterial cells. Bacteria were then kept for 1 hr at 37°C and plated on ampicillin-containing medium. After overnight incubation at 37°C, several colonies were picked up in separate tubes containing LB and incubated overnight at 37°C. The following day, plasmid DNA was extracted using the DNA Purification Kit (MN #740410) according to the manufacturer’s instructions; DNA from positive colonies was identified upon gel electrophoresis and sent for sequencing for final confirmation.
Same procedure was conducted to build the mutated fra reporter plasmid.
Following oligonucleotides were used:
fra WT forward:
5’ GAGAGGTACCGTGTCCAAAAATGCGGGTCTGTTTCTCG 3’
fra WT reverse:
5’GAGAGCTAGCGTTAAGACAAACATGCAGGCATAAAGACATG 3’
fra Mutant forward:
5’GAGAGGTACCGTGTCCAAAAAAAAAAACTGTTTCTCGAAATTGAGTT 3’
fra Mutant reverse:
5’GAGAGCTAGCGTTAAGACAAACAAAAAAAAATAAAGACATGAAATGGATG 3’
Similarly, we built unc5 WT and mutant plasmids with flanking restriction sites for EcoR1 at the 5’ extremity and Kpn1 at the 3’ extremity.
Following oligonucleotides were used:
unc5 WT forward:
5’ GAGAGAATTCTCGTTTTCCCGTTTAGGGCA 3’
unc5 WT reverse:
5’ GAGAGGTACCACTAGCGCTCACCACAGTTC 3’
unc5 Mutant forward 1:
5’ GTGTGAACAGTGATATAAAGTGCACCGTGTAAAAAAAATAGAGATACCT 3’
unc5 Mutant reverse 1:
5’ TTCGTGTGGCACTAGGTTAGGTATCTCTATTTTTTTTACACGGTGCACT 3’
unc5 Mutant forward 2:
5’ ATAAAAACAAGCCGCACACACAGTAGCACAAAAAAAAAAAGGGGCGCAC 3’
unc5 Mutant reverse 2:
5’ CATCGGACGACCACTGCAGTGCGCCCCTTTTTTTTTTTGTGCTACTGTG 3’
Co-transfection and Western blot assays
Request a detailed protocolCo-transfections in S2 cells were carried out using Lipofectamine (Invitrogen). 6 × 106 cells were cultured in 6-well plates containing Schneider medium. In each well, cells were transfected with 1 μg of fra WT or
mutant reporter plasmid, 1 μg of pPAC-lacZ (Flici et al., 2014; Cattenoz et al., 2016) as a transfection control, 0.5 μg or 1 μg or 2 μg of pPac gcm (Cattenoz et al., 2016) expression vector and pPac ‘empty’ (Flici et al., 2014) to make up the volume up to 4 μg. Cells were collected 48 hr after transfection, washed in cold PBS and resuspended in lysis buffer (25 mM Tris-phosphate pH7.8, 2 mM EDTA, 1 mM DTT, 10% glycerol, 1% Triton X-100). Total protein extract was obtained by four freezing-thawing cycles in liquid nitrogen and centrifugation at 4°C at 13,000 g. Protein expression was detected as per standard Western blot procedures. Primary antibodies used were as follows: mouse anti-β-Gal (1/2000, Sigma), rabbit anti-GFP (1/5000, Molecular Probes); and rabbit anti-HRP (1/5000, Jackson ImmunoResearch) were used as secondary antibodies. Each experiment was performed in triplicate.
β-Gal assays were performed to measure the levels of LacZ for each replicate. 20 μL of protein extract mixed with 50 μL of β-Gal assay buffer (60 mM Na2PO4, 40 mM NaH2PO4, 10 mM KCl, 1 mM MgCl2, and 50 mM β-mercaptoethanol) containing ONPG was incubated at 37°C. The reaction was stopped by adding 50 μL of 1M Na2CO3 once the solution turned yellow, DO was analyzed at 415 nm. The levels of GFP were normalized to the LacZ value in each blot and were quantified using ImageJ software. The background was subtracted from each band value and the average was calculated.
Same protocol was used for unc5 S2 cell co-transfection assays.
Reverse transcription and qRT-PCR
Request a detailed protocolTotal RNA was extracted from S2 cells using tri-reagent, 1 μg of purified RNA was reverse transcribed by SuperScript II. qPCR was performed with the Roche LightCycler 480 and Sybr Green Master mix (Roche) using the following oligonucleotides:
GFP forward: ACATGAAGCAGCACGACTTCT
GFP reverse: TTCAGCTCGATGCGGTTCA
Gcm WT forward: 5’GAGAGATCTTATCCCGATCCCCTAGC3’
Gcm WT reverse: 5’CTACTACTACAGCAATACGGG3’
LacZ forward: TGTGCCGAAATGGTCCATCA
LacZ reverse: GTATCGCCAAAATCACCGCC
For each gene, the expression levels were automatically calculated (LightCycler480 Software, release 1.5.0) by calibration to gene-specific standard curves generated on input cDNAs. Collected values, derived from three amplification reactions, each performed in three independent experiments, were normalized to β-gal mRNA amounts.
Immunolabeling and antibodies
Request a detailed protocolPupae of desired stage were collected and fixed in 4% PFA PBS (paraformaldehyde in phosphate buffer saline) overnight at 4°C. They were dissected in PBT (PBS Triton-X100, 0.3%) and wings were given four quick washes of 10 min in PBT and were incubated in the blocking reagent PBT-NGS (5% normal goat serum in PBT) for 60 min at room temperature on a planar shaker. Samples were then incubated overnight in primary antibodies (diluted in PBT-NGS): mouse-anti-Repo labels glia (1:800) and mouse-anti-22c10 labels neurons (1:1000) (DSHB), rabbit-anti-HRP labels neurons (1:1000) (DSHB), chicken-anti-GFP (1:1000) (Abcam), rat anti-Elav labels neurons (1:1000) (DSHB), rabbit-anti-Unc5 and rabbit-anti-Fra (1:500) were gifts from Benjamin Altenhein. After four washes in PBT, wings were incubated for 2 hr at room temperature in secondary antibodies (1:500) raised in mouse, rat, rabbit or chicken and coupled to Cy3, Cy5 or FITC fluorescent dyes diluted in PBT-NGS. Following a final wash in PBT, wings were mounted on slides in Aqua- Poly/Mount medium (Polysciences Inc.).
In vivo imaging
Request a detailed protocolTime-lapse analyses were performed using the standard procedure as described by Aigouy et al. (2004, 2008), Kumar et al. (2015), and Soustelle et al. (2008). Photo bleaching was avoided by using low magnification and reduced exposure time. Maximum projections for time-lapse and confocal images were obtained by using the ImageJ software. Images were annotated by using Adobe Photoshop and Illustrator.
Quantification methods were defined previously (Berzsenyi et al., 2011; Kumar et al., 2015). Briefly, based on the time-lapse movies, the migratory process was subdivided in three phases: the earliest one describes migration initiation around 18 hAPF; this time point defines the movement of the soma of the first cell at the front of the chain. The intermediate phase identifies the time at which the glial chain reaches the level of the nerve on the costa, which is around 22 hAPF; the latest phase refers to migration completion, upon connection of the chain with the proximal glia located on the radius nerve, at around 28 hAPF, see panels shown in Figure 1a–c. The time points (hAPF) that were used to define these three migratory phases were calculated on the basis of time-lapse movies performed in at least 10 control animals (n≥10). In order to quantify the migratory defects, the differences in time between the control and mutant wings were compared and the graphs were plotted using the standard Student's t test method.
Statistical analysis
Request a detailed protocolThe number of wings dissected for each experiment and genotype were more than or equal to 30. The Migratory Index (Kumar et al., 2015) defines the percentage of wings in which glial cells have completed migration at a given time point (28 hAPF in most cases). Graphs were made using Prism software, and the Student's t test method was used for the comparison between two different experimental sets. Error bars indicate the standard error mean (s.e.m.). p values: ***p<0.0001; **p <0.001; *p<0.05.
References
-
Glial chain migration requires pioneer cellsJournal of Neuroscience 28:11635–11641.https://doi.org/10.1523/JNEUROSCI.1361-08.2008
-
Analysis of Ras-induced overproliferation in Drosophila hemocytesGenetics 163:203–215.
-
Glial differentiation does not require a neural ground stateDevelopment 125:3189–3200.
-
glide/gcm is expressed and required in the scavenger cell lineageDevelopmental Biology 191:118–130.https://doi.org/10.1006/dbio.1997.8702
-
Netrins guide Drosophila commissural axons at short rangeNature Neuroscience 9:188–194.https://doi.org/10.1038/nn1625
-
Transcriptional regulation of string (cdc25): a link between developmental programming and the cell cycleDevelopment 120:3131–3143.
-
Development of the zebrafish lateral lineCurrent Opinion in Neurobiology 14:67–73.https://doi.org/10.1016/j.conb.2004.01.012
-
Development and organization of glial cells in the peripheral nervous system of Drosophila melanogasterDevelopment 117:895–904.
-
Collective cell migration: "all for one and one for all"Journal of Neurogenetics 28:190–198.https://doi.org/10.3109/01677063.2014.896911
-
A novel role of the glial fate determinant glial cells missing in hematopoiesisThe International Journal of Developmental Biology 53:1013–1022.https://doi.org/10.1387/ijdb.082726cj
-
Netrin-1 is a chemorepellent for oligodendrocyte precursor cells in the embryonic spinal cordJournal of Neuroscience 23:3735–3744.
-
Roundabout signalling, cell contact and trophic support confine longitudinal glia and axons in the Drosophila CNSDevelopment 128:207–216.
-
Modes and regulation of glial migration in vertebrates and invertebratesNature Reviews Neuroscience 10:769–779.https://doi.org/10.1038/nrn2720
-
Netrins: versatile extracellular cues with diverse functionsDevelopment 138:2153–2169.https://doi.org/10.1242/dev.044529
-
Development: new wrinkles on genetic control of the MBTCurrent Biology 23:R65–67.https://doi.org/10.1016/j.cub.2012.12.012
-
Compartments in the wing of Drosophila: a study of the engrailed geneDevelopmental Biology 50:321–337.https://doi.org/10.1016/0012-1606(76)90155-X
-
E-cadherin plays an essential role in collective directional migration of large epithelial sheetsCellular and Molecular Life Sciences 69:2779–2789.https://doi.org/10.1007/s00018-012-0951-3
-
Slit2 regulates the dispersal of oligodendrocyte precursor cells via Fyn/RhoA signalingJournal of Biological Chemistry 287:17503–17516.https://doi.org/10.1074/jbc.M111.317610
-
Guiding neuronal cell migrationsCold Spring Harbor Perspectives in Biology 2:a001834.https://doi.org/10.1101/cshperspect.a001834
-
Border-cell migration: the race is onNature Reviews Molecular Cell Biology 4:13–24.https://doi.org/10.1038/nrm1006
-
DE-Cadherin is required for intercellular motility during Drosophila oogenesisThe Journal of Cell Biology 144:533–547.https://doi.org/10.1083/jcb.144.3.533
-
Schwann cell transplantation for repair of the adult spinal cordJournal of Neurotrauma 23:453–467.https://doi.org/10.1089/neu.2006.23.453
-
Improved and expanded Q-system reagents for genetic manipulationsNature Methods 12:219–222.https://doi.org/10.1038/nmeth.3250
-
Collective cell migrationAnnual Review of Cell and Developmental Biology 25:407–429.https://doi.org/10.1146/annurev.cellbio.042308.113231
-
Directional guidance of oligodendroglial migration by class 3 semaphorins and netrin-1Journal of Neuroscience 22:5992–6004.
-
Netrin 1 mediates spinal cord oligodendrocyte precursor dispersalDevelopment 130:2095–2105.https://doi.org/10.1242/dev.00424
-
Glide directs glial fate commitment and cell fate switch between neurones and gliaDevelopment 122:131–139.
-
Transcription factor interplay during Drosophila haematopoiesisThe International Journal of Developmental Biology 54:1107–1115.https://doi.org/10.1387/ijdb.093054lw
-
Rab8, POSH, and TAK1 regulate synaptic growth in a Drosophila model of frontotemporal dementiaJournal of Cell Biology 208:931–947.https://doi.org/10.1083/jcb.201404066
-
Mapping proteolytic cancer cell-extracellular matrix interfacesClinical & Experimental Metastasis 26:289–298.https://doi.org/10.1007/s10585-008-9190-2
-
Dynamic regulation of axon guidanceNature Neuroscience 4 Suppl:1169–1176.https://doi.org/10.1038/nn748
-
WASp is required for the correct temporal morphogenesis of rhabdomere microvilliJournal of Cell Biology 164:417–426.https://doi.org/10.1083/jcb.200307048
Decision letter
-
Hugo J BellenReviewing Editor; Baylor College of Medicine, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for submitting your article "The Glide/Gcm fate determinant controls initiation of collective cell migration by regulating Frazzled" for consideration by eLife. Your article has been reviewed by three reviewers who are really familiar with the topic. They all three are positive about your work and we would be glad to consider a revised manuscript for publication in eLife. The evaluation has been overseen by a Reviewing Editor (Hugo Bellen) and a Senior Editor.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
Gupta et al. investigate here the function of netrins and their frazzled (Fra) and Unc5 receptors in glial cell chain migration in the Drosophila wing. Their findings show that frazzled is the major Netrin receptor facilitating initiation of glial migration, and that an early glial fate determinant TF Gcm directly activates frazzled expression. Furthermore, Netrin B but not Netrin A, apparently forms a long-range gradient that guides glia in their migration toward the CNS along peripheral nerves., and Unc5 serves as a repulsive receptor for migration of these wing glia. Transcriptional regulation of axon guidance cues and their receptors remains to be studied in detail. Therefore, this manuscript provides evidence extending our understanding of how these cues and receptors are utilized by glia and how they are transcriptionally regulated. In general, this paper appears to be one that could be of interest to the field and appropriate for eLife following revisions that address the concerns raised below.
1) The lab has previously demonstrated that Pioneer cells trigger chain migration. It would greatly improve this initial part of the paper, if the data on expression of Frazzled on migrating glial cells would be shown more convincingly. Instead of showing one confocal image in Figure 1E-G show nerves of wild type, and nerves of animals with glial or neuronal knockdown of Frazzled. There is a MiMIC insertion line available which would allow easy imaging of Frazzled protein expression (Nagarkar-Jaiswal et al., eLife, 2015).
Fra and Unc5 staining is not easy to see and compare across different conditions; in Figure 2D'-H', Figure 5C'-E', Figure 7A'-D', etc. GFP does not seem to match between different panels, such as in Figure 5C compared to C', or Figure 5D compared to D'. The white used here makes it difficult to separate individual RGB colors when these images are analyzed. It would be far easier for the reader to appreciate these results if there were separate panels for each staining in these figures, especially for Fra and UNC5 ICC. It would be even better to quantify the change in Fra staining in glia vs. neurons in the Fra-KD and GOF experiments using confocal optical sections.
2) It would be quite interesting and relevant in the light that collective migration is triggered by a few pioneers if a clonal analysis could be performed using the existing frazzled null allele. Do all cells have to express Fra? Provided that the wild type twin spot is allowed to stay, do these glial cells overtake mutant cell in front of them?
3) The observed genetic interaction between GCM and FRA is interesting. The fact that the fra3 phenotype is much worse in the background of a heterozygous rA87 insertion – which is a viable enhancer trap in the gcm promoter region and thus a very weak hypomorph – makes me wonder how the fra3 phenotype would look in trans to a gcm null allele? Is it possible to quantify the reduction of gcm mRNA in the different genetic backgrounds?
The authors use a migratory index to characterize the different genotypes. This is clear and straightforward but I wonder given the dramatic differences between 28APF and 24APF whether general developmental problems caused by the expression of the different transgenes also in many cells outside the nervous system (e.g. gcm is active also in the hemocyte lineages?) contribute to the phenotype. This could be controlled by for example comparing the time it takes to form the puparium and to eclose from the pupal cases.
4) Fra is also expressed at high levels (Figure 1F) in the axons that glia migrate along, raising the question of whether neuronal Fra influences glia migration, either directly or indirectly. What happens if Fra is knocked down specifically in neurons?
The authors claim that changes in Gcm expression lead to changes of Fra in glia. However, changes of Gcm expression in glia (using gcm-GAL4) seem to have profound effects on neuronal Fra expression. For example, in Figure 5D'-D', knockdown of Gcm using gcm-GAL4 leads to a decrease of Fra in L1, and even in L3 which is not wrapped by GFP+ processes in these images. Again, in Fra GOF (Figure 5E',E'), the Fra signal co-localizes more with red 22C10 staining than with green GFP. Is this due to leaky expression of gcm-GAL4 in neurons at earlier stages? Or, does this suggest non-cell autonomous functions of glial Gcm on neuronal Fra expression?
5) Pioneer glia vs. follower glia: during glia chain migration, pioneer glia at the tip of the chain are suggested to actively lead followers. Fra, however, appears to be detected all along the chain (Figure 1G and Figure 1—figure supplement 1), raising the question as to whether follower glia also use Fra to initiate and direct their migration independent of pioneer glia. It is technically challenging to use GAL4 drivers to differentially target pioneer or follower glia. However, have the authors considered using MARCM (Aigouy et al., 2008) to randomly target different glia in the wing and monitor their migratory behavior in either WT or Fra mutant clones? This could also address the cell autonomous vs cell-type autonomous (RNAi experiment) requirement for Fra in wing glia migration.
6) It is surprising that the NetB mutant (NetBΔ/Y, Figure 6A) has migration defects of similar severity to the Fra/+ heterozygote, and has milder effects than many other conditions, such as gcm>fra-KD and gcm>gcm-KD. This could be explained by NetB signaling through Unc5 to inhibit glial migration in Fra LOF. If this is true, the NetB mutant should enhance the Fra/+ heterozygote phenotype. Is this the case? Also, does ectopic expression of NetB in a different wing sub-region using drivers such as en-GAL4 or ap-Gal4 mis-direct glial migration?
7) In the subsection “fra plays an instructive role in L1 glia migration”, the authors claim 'Finally, we measured the migration speed at the time of initiation and it was not higher in Fra GOF glia (10 μm/h) than in control glia (40 μm/h)'. This is confusing since the numbers do appear to be significantly different. What is the actual speed of glial migration in Fra GOF vs. control? The authors should also include this important information in graphs using either the speed or duration of glial migration from initiation to completion of migration.
8) The authors try to rescue RNAi phenotypes by co-expression of RNAi and UAS-GOF constructs, stating that this 'excludes the possibility of off target effects'. This is not accurate since this rescue effect could be due to titration of siRNA away from the endogenous mRNA by overexpression of the exogenous mRNA. To directly rule out off target effects the rescue experiments should be done by using either a UAS construct bearing mutants that are insensitive to RNAi or perhaps an RNAi line that targets non coding regions.
9) Throughout the study, there are no details or explanations on the methods of quantification. For instance, Results: subsection “fra plays an instructive role in L1 glia migration”, third paragraph – this is very unclear! how exactly the authors define "initiation? and how was the quantification done? Figure 2J, how does author calculate the hAPF for initiation, costa reach, and complete migration? what is the n number? how early or late in migration do the authors consider a defect in initiation? These all have to be explained fully in the text (the Methods section or the legends). Without these knowledge, the readers will have a hard time to follow and do not understand how the conclusion was reached. Also, for the time-lapse figure, it is not clear what does 22:48 hAPF mean since there are different genotypes shown, and I cannot find any information on this in the legends nor methods.
10) Panel organization: In Figures 2, 3, 6, 7, in which the authors characterize glia migration phenotypes (MI), genotypes are assigned in different panels with the same or different time points. These arrangements make it difficult to compare different genotypes. For example, Figure 2A and B could be combined, and the same is true for Figures 3A and B and Figures 7E and F. In some panels, more genotypes should be included. For example, gcm-GAL4>gcm (GOF) is needed in Figure 3B, and NetB-GAL4> NetB needs to be added to Figure 6C. In Figure 6F and G, all three genotypes should be compared in both time points.
11) Since Gcm is a very important factor in controlling neuronal and glial cell fates, it is crucial to address clearly whether this conversion function is in any part related to the migration phenotype the authors are trying to assessing in the study. However, there is only one experiment provided which wasn't able to give a clear answer. For instance, Figure 4—figure supplement 1: it is not common to analyze the co-positive cells of gcm>GFP and Elav as Elav stained for nucleus and GFP labels the membrane. Can't the author do the costaining of Repo and Elav? this might be more clear and definitive. Also, what is the percentage of gcm KD wings not having extra Elav positive cells? Is it high or low? Isn't it plausible that a high percentage of the gcm KD wings should have extra Elav cells? The gcm KD wings with normal Elav positive cells might as well mean that gcm expression wasn't knock-down and in that case why is the MI still defective? How do the authors reconcile the discrepancies? This section has to be explained clearly since glia-neuron switch is a very important function for Gcm and whether the migration function depends on this needs to be clearly addressed with clear, precise experimental results.
[Editors' note: further revisions were requested prior to acceptance, as described below.]
Thank you for resubmitting your work entitled "The Glide/Gcm fate determinant controls initiation of collective cell migration by regulating Frazzled" for further consideration at eLife. Your revised article has been evaluated by a Senior editor, a Reviewing editor, Hugo Bellen, and three reviewers.
The manuscript has been improved but there are some remaining issues that need to be addressed. The reviewers discussed their reviews on line and came to a consensus that they all agreed with reviewer 3. The reviewers feel that some of the points that they raised where not addressed and Reviewer #3 summarized these concerns.
Reviewer #3:
I read the revised paper and cannot support publication at this stage due to the poor demonstration that fra is expressed by the glia and the lack of a clear demonstration that the phenotype is due to loss of glial fra function. I admit there are some hints supporting this suggestion but I do not think the data are conclusive.
Expression of fra in glia. The authors added a new Figure 1H-J which shows fra expression upon glial knockdown. For me this is not showing that fra is expressed in glia. There is as much background as in the wild type and the control experiment, neuronal knockdown, is not shown based on the argument that elav is also expressed by glia. Please, use nsybGal4 which is specific – and/or generate clones! The quality of the antibody stainings is not consistent and for example in Figure 2E glial knockdown is also removing the neuronal Fra expression. For me the quality of the data is still simply not good enough.
The phenotypes are based on somewhat odd genetic backgrounds. May be there is a fra haploinsufficiency and heterozygous fra/+ animals show a phenotype. But to really determine the Fra phenotype, a true loss of function analysis is required. The argument that MARCM clones are hard to get is weird as the same lab has generated these clones before in another study. Glia-specific fra knockdown animals show roughly the same phenotype as heterozygous fra3/+ animals (MI 58 vs.61). Still only the RNAi induced phenotype can be completely rescued. The authors conclude that this may be due to the reduction of neuronal Fra expression (which was not tested). However, one could also argue that the authors are looking at a genetic background effect. Thus, a MARCM analysis is mandatory.
We had asked for the effect of gcm on the migratory phenotype of fra3/+ heterozygous animals. Although the authors provide arguments for an allelic series of the different gcm alleles they do not show the effect of a true gcm loss of function allele on the migratory index. Provided that the presence of one very weak, homozygous viable, gcm allele has such a profound effect on the migratory index, what would a true null allele do? The allelic series shown in the rebuttal letter indicates that homozygous gcm>GFP animals have the same effect as gcm>GFP gcmN7-4. Since gcmN7-4 is a null this implies that gcm>GFP is a null as well – but why do they survive then? Something is odd here.
https://doi.org/10.7554/eLife.15983.033Author response
[…]
1) The lab has previously demonstrated that Pioneer cells trigger chain migration. It would greatly improve this initial part of the paper, if the data on expression of Frazzled on migrating glial cells would be shown more convincingly. Instead of showing one confocal image in Figure 1E-G show nerves of wild type, and nerves of animals with glial or neuronal knockdown of Frazzled. There is a MiMIC insertion line available which would allow easy imaging of Frazzled protein expression (Nagarkar-Jaiswal et al., eLife, 2015).
Fra and Unc5 staining is not easy to see and compare across different conditions; in Figure 2D'-H', Figure 5C'-E', Figure 7A'-D', etc. GFP does not seem to match between different panels, such as in Figure 5C compared to C', or Figure 5D compared to D'. The white used here makes it difficult to separate individual RGB colors when these images are analyzed. It would be far easier for the reader to appreciate these results if there were separate panels for each staining in these figures, especially for Fra and UNC5 ICC. It would be even better to quantify the change in Fra staining in glia vs. neurons in the Fra-KD and GOF experiments using confocal optical sections.
As suggested by the referees we made sure that the labeling was indeed less prominent in the KD and more intense in the overexpressing animals. We have now included the data that also shows the glial knockdown of fra in a 18 hAPF old wing in Figure 1H-J and we have included panels that only shows Fra or Unc5 labeling in Figure 2C'-G'”, Figure 5C‐E”' and Figure 7A‐D”'.
The difference in GFP intensity in the merged panels is due to the strong red signal. Specifically, at late stages these differences are due to the fact that glia wrap neurons (which are in red) completely. In addition, we have quantified Fra staining in Figure 1E-J between control and fra KD L1 glia by including heat maps. Note that in the KD the neuronal (but not the glial) Fra signal is as high if not higher than in the wt.
Concerning the MiMIC insertion, we tried and stay consistent, using the same antibody that was used in previous papers, in particular in the paper on glial migration from Altenhein and collaborators (References on the used antibody: Von Hilchen et al., 2010; Timofeeev et al., 2012; Joly et al., 2007; Kolodziej et al., 1996; Orr et al., 2014; Gong et al., 1999; Pert et al., 2015).
2) It would be quite interesting and relevant in the light that collective migration is triggered by a few pioneers if a clonal analysis could be performed using the existing frazzled null allele. Do all cells have to express Fra? Provided that the wild type twin spot is allowed to stay, do these glial cells overtake mutant cell in front of them?
We see Fra expression in all L1 glia, which is shown in Figure 1—figure supplement 1(subsection “Frazzled expression in the glia of the developing Drosophila wing”, second paragraph). The clonal analysis is a valuable point, however we think that the data from this experiment will not be easy to interpret. Fra is also expressed in neurons, therefore getting clones that only affect glia will be difficult (low division rate) and even more to get clones that specifically affect the pioneer cells. We did in the past get MARCM clones in pioneer cells, but this is a very rare event. This was useful to analyze the morphology of pioneers, but the analysis of migration phenotypes requires high number of events in order to perform the statistics. More importantly, we have previously shown that several pioneers are at work, therefore mutating one or two of them may not be sufficient to see a phenotype.
3) The observed genetic interaction between GCM and FRA is interesting. The fact that the fra3 phenotype is much worse in the background of a heterozygous rA87 insertion – which is a viable enhancer trap in the gcm promoter region and thus a very weak hypomorph – makes me wonder how the fra3 phenotype would look in trans to a gcm null allele? Is it possible to quantify the reduction of gcm mRNA in the different genetic backgrounds?
The authors use a migratory index to characterize the different genotypes. This is clear and straightforward but I wonder given the dramatic differences between 28APF and 24APF whether general developmental problems caused by the expression of the different transgenes also in many cells outside the nervous system (e.g. gcm is active also in the hemocyte lineages?) contribute to the phenotype. This could be controlled by for example comparing the time it takes to form the puparium and to eclose from the pupal cases.
Concerning the interaction, gcmrA87is a weaker allele than the gcm‐Gal4 line and, accordingly gcmrA87/fra3 has a weaker migratory phenotype that gcm‐Gal4/ fra3. These data are now shown together in Figure 3, to comply with reviewers’ suggestion. For the benefit of the reviewers, we also show here below the rate of conversion and/or the number of repo positive cells in different allelic combinations (gcm-Gal4/ gcmrA87, gcm-Gal4/ gcm-Gal4, gcmrA87/ gcmN7-4 and gcm‐Gal4>gcm RNAi). The inverse correlation between the two phenotypes in the different backgrounds further sustains our hypothesis. We believe that the phenotypic analysis of the allelic series is more sensitive than RNA quantification, given the low number of gcm expressing cells. For the sake of simplicity, we do not think that the data here below should be included. Should the referee feel it necessary, however, we are eager to include them in Figure 5—figure supplement 1.
We did not observe any general developmental defect, such as delay in eclosion from the pupal case.

To rule out a possible hemocyte contribution in glia migration we downregulated gcm using a hemocyte driver collagen‐Gal4 and found no effect on glia migration. Further, we also knocked down fra using a temperature sensitive Gal80 construct of gcm driver and specifically knocked down fra post embryonically. Since gcm is only expressed in blood cells embryonically, this eliminates any possible contribution of hemocyte Gcm expression to the observed phenotypes. These analyses confirm that gcm and fra directly affect glia migration as in these conditions too glia migration is significantly delayed (Figure 5—figure supplement 1D and E).
4) Fra is also expressed at a high levels (Figure 1F) in the axons that glia migrate along, raising the question of whether neuronal Fra influences glia migration, either directly or indirectly. What happens if Fra is knocked down specifically in neurons?
The authors claim that changes in Gcm expression lead to changes of Fra in glia. However, changes of Gcm expression in glia (using gcm-GAL4) seem to have profound effects on neuronal Fra expression. For example, in Figure 5D'-D', knockdown of Gcm using gcm-GAL4 leads to a decrease of Fra in L1, and even in L3 which is not wrapped by GFP+ processes in these images. Again, in Fra GOF (Figure 5E', E'), the Fra signal co-localizes more with red 22C10 staining than with green GFP. Is this due to leaky expression of gcm-GAL4 in neurons at earlier stages? Or, does this suggest non-cell autonomous functions of glial Gcm on neuronal Fra expression?
The neuronal expression of Fra is an important issue and we thank the reviewer(s) for the comment. We analyzed L1 axons of 17 hAPF fra3 wings and found that about 30% of them show a transient phenotype of delayed axonal growth (see now Figure 3—figure supplement 2B-D). We also analyzed wings where fra is knocked down using an Elav driver and got a very strong migratory phenotype (Data not shown in the manuscript). However, the available Elav drivers (we tested two of them) are also expressed in wing glia and epithelial cells (Author response image 2). Other known drivers are either expressed in neuronal subsets (motoneurons) or too late for our purpose. Thus, in the wing, we were unable to use a driver that specifically activates or knocks down gene expression in neurons. These analyses suggested that downregulating fra in neurons may indirectly affect glia migration and prompted us to use a driver that specifically affects glia.

25 hAPF immunolabeled wing of the elav-Gal4,UASmCD8GFP transgenic line (anti-GFP in green and anti-Repo in red).
Note the presence of the mCD8GFP in the axons, in the glia (Repo-positive cells in red, arrows) and in the epithelium (asterisk).
We think there is no change in the neuronal expression of Fra upon gcm KD or GOF. Figure 5D‐D’” is only one example and the decrease might be due to the fact that these images are few sections and not the whole projection. Concerning the role of Gcm on neuronal expression of Fra. We never observed the Gcm driver being expressed in neurons. For all the above reasons, we do not think there is a non-cell autonomous effect.
5) Pioneer glia vs. follower glia: during glia chain migration, pioneer glia at the tip of the chain are suggested to actively lead followers. Fra, however, appears to be detected all along the chain (Figure 1G and Figure 1—figure supplement 1), raising the question as to whether follower glia also use Fra to initiate and direct their migration independent of pioneer glia. It is technically challenging to use GAL4 drivers to differentially target pioneer or follower glia. However, have the authors considered using MARCM (Aigouy et al., 2008) to randomly target different glia in the wing and monitor their migratory behavior in either WT or Fra mutant clones? This could also address the cell autonomous vs cell-type autonomous (RNAi experiment) requirement for Fra in wing glia migration.
The reviewers raise an interesting point, i.e. what is the role of Fra, since it is all along the L1 chain. We believe that the suggested MARCM clonal analysis is beyond the scope of this manuscript. We indeed used MARCM to follow the morphology of single pioneers but mutant clones on single pioneer glial cells are very rare, which makes it hard to draw conclusions that are statistically significant. If we go for big clones, they will also affect neurons, which will not help assessing the cell autonomous effect. See also answer to Q2.
To comply with the comment, we asked whether the ligand Net B accumulates on specific glia, based on the recent paper from Pert et al., (Biology Open 2015) showing Net B accumulation in Fra expressing cells. Immunolabeling assays on young wings indeed show preferential Net B accumulation in some cells at the front of the chain. We speculate that the receptor is all over the glia but pioneer glia are closer to the ligand than the other glia and some pioneers accumulate it, perhaps due to receptor‐mediated endocytosis. This may in turn triggers the ‘promigratory’ downstream pathway in those cells.

6) It is surprising that the NetB mutant (NetBΔ/Y, Figure 6A) has migration defects of similar severity to the Fra/+ heterozygote, and has milder effects than many other conditions, such as gcm>fra-KD and gcm>gcm-KD. This could be explained by NetB signaling through Unc5 to inhibit glial migration in Fra LOF. If this is true, the NetB mutant should enhance the Fra/+ heterozygote phenotype. Is this the case? Also, does ectopic expression of NetB in a different wing sub-region using drivers such as en-GAL4 or ap-Gal4 mis-direct glial migration?
We thank the reviewer(s) for the thoughtful comments. The reviewer is correct, the NetB phenotype might be milder than that observed in fra3/+ because of the effects of NetB on unc5. Along this line, Figure 8—figure supplement 1C shows that UAS‐RNAi for both unc5 and fra has a weaker phenotype than UAS‐fra‐RNAi alone. As per the double mutant NetBΔ/Y; fra3/+, this genotype would not inform on the role of Unc5, because in a hemizygous NetBΔ/Y animal both Fra and Unc5 should be absent anyway, therefore removing the fra receptor would not make any difference. We thought that this experiment may nevertheless be informative and tell whether fra is also submitted to other signals, as suggested by recent studies (Neuhaus-Follini et al; 2015). We hence performed the suggested experiment and it does seem that fra also depends on other signals, see data on the right graph. We would not include this data in the main text, however, as they go beyond the scope of the manuscript, which is to show the impact of a fate determinant on late aspects of cell specification. As requested, we show the data on the ectopic expression of NetB using en‐Gal4 in the posterior wing compartment and using GMR 29F05-Gal4 that is expressed in the distal part of the wing. These experiments show that ectopic NetB does not affect glia migration (Now Figure 8—figure supplement 1A). Thus, NetB does not act as a guidance molecule triggering directionality. This is also discussed in the manuscript on page 23.

7) In the subsection “fra plays an instructive role in L1 glia migration”, the authors claim 'Finally, we measured the migration speed at the time of initiation and it was not higher in Fra GOF glia (10 μm/h) than in control glia (40 μm/h)'. This is confusing since the numbers do appear to be significantly different. What is the actual speed of glial migration in Fra GOF vs. control? The authors should also include this important information in graphs using either the speed or duration of glial migration from initiation to completion of migration.
To comply with the comments and clarify our statements, we performed further analyses that are now incorporated in Figure 3—figure supplement 2. The time-lapse analysis shows that fra GOF glia start moving earlier than control glia (Figure 2I, Figure 3—figure supplement 1E). The gcm>fra GOF glia phenotype is solely due to altered initiation rather than to acceleration and we have added a speed graph of control and fra GOF glia in Figure 3—figure supplement 2A. The speed analysis was done on the most proximal glial cell of the chain (μm/h, y-axis; hAPF, x‐axis). The distance covered by the front cell was measured by analyzing the position of the glial soma every hour, as already shown in Kumar et al., 2015. We followed the cell starting from the point of initiation till complete migration. The migratory speed increases both in control and fra GOF glia until they reach the level of costa and then it decreases, however, fra GOF glia display a more uniform speed than that of control glia: they accelerate less in the early phases and slow down less in the late phases.
8) The authors try to rescue RNAi phenotypes by co-expression of RNAi and UAS-GOF constructs, stating that this 'excludes the possibility of off target effects'. This is not accurate since this rescue effect could be due to titration of siRNA away from the endogenous mRNA by overexpression of the exogenous mRNA. To directly rule out off target effects the rescue experiments should be done by using either a UAS construct bearing mutants that are insensitive to RNAi or perhaps an RNAi line that targets non coding regions.
We take the point made by the reviewer. However, we would like to emphasize that the RNAi data are further supported by the opposite phenotype obtained in gain of function experiments and by the fact that we could rescue the fra3mutant phenotype by reintroducing a UAS‐fra transgene in glia (Figure 3—figure supplement 1A). We nevertheless played down the interpretation of the RNAi data in the main text, to comply with the comment.
9) Throughout the study, there are no details or explanations on the methods of quantification. For instance, Results: subsection “fra plays an instructive role in L1 glia migration”, third paragraph – this is very unclear! how exactly the authors define "initiation? and how was the quantification done? Figure 2J, how does author calculate the hAPF for initiation, costa reach, and complete migration? what is the n number? how early or late in migration do the authors consider a defect in initiation? These all have to be explained fully in the text (the Methods section or the legends). Without these knowledge, the readers will have a hard time to follow and do not understand how the conclusion was reached. Also, for the time-lapse figure, it is not clear what does 22:48 hAPF mean since there are different genotypes shown, and I cannot find any information on this in the legends nor methods.
As the referee suggests, we included a detailed description of how the time-laps analysis was performed and how the quantification was made in the Methods section of the manuscript (in vivoimaging paragraph). These methods were already defined in our previous publications (Berzseny et al., 2011, Kumar et al., 2015).
“Based on the time-lapse videos, the migratory process has been subdivided in three phases: the earliest one describes migration initiation around 18 hAPF; this time point defines the movement of the soma of the first cell at the front of the chain. […] In order to quantify the migratory defects, the differences in time between the control and mutant wings were compared and the graphs were plotted using the standard Student's t test method.”
We have used the above phases to select the images from the time-lapsevideos. Typically, Figure 2J, Kshows a snap shot at 21h48, which corresponds to the time by which wild type glia have reached the Costa, in average. Figure 3D-F,shows a snap shot at 22h48, which corresponds to the time by which GOF glia have completed migration, in average. We have clarified this in the text.
10) Panel organization: In Figures 2, 3, 6, 7, in which the authors characterize glia migration phenotypes (MI), genotypes are assigned in different panels with the same or different time points. These arrangements make it difficult to compare different genotypes. For example, Figure 2A and B could be combined, and the same is true for Figures 3A and B and Figures 7E and F. In some panels, more genotypes should be included. For example, gcm-GAL4>gcm (GOF) is needed in Figure 3B, and NetB-GAL4> NetB needs to be added to Figure 6C. In Figure 6F and G, all three genotypes should be compared in both time points.
As per referee’s suggestion we have reorganized Figure 2A, B; Figure 3A, B and Figure 7E, Finto one panel. We have also included the migratory index of gcm>gcm GOF in Figure 3Aand the migratory index of NetB>NetB GOF in Figure 6C. Finally, as requested, we have also included all the genotypes in Figure 6F and G.
11) Since Gcm is a very important factor in controlling neuronal and glial cell fates, it is crucial to address clearly whether this conversion function is in any part related to the migration phenotype the authors are trying to assessing in the study. However, there is only one experiment provided which wasn't able to give a clear answer. For instance, Figure 4—figure supplement 1: it is not common to analyze the co-positive cells of gcm>GFP and Elav as Elav stained for nucleus and GFP labels the membrane. Can't the author do the costaining of Repo and Elav? this might be more clear and definitive. Also, what is the percentage of gcm KD wings not having extra Elav positive cells? Is it high or low? Isn't it plausible that a high percentage of the gcm KD wings should have extra Elav cells? The gcm KD wings with normal Elav positive cells might as well mean that gcm expression wasn't knock-down and in that case why is the MI still defective? How do the authors reconcile the discrepancies? This section has to be explained clearly since glia-neuron switch is a very important function for Gcm and whether the migration function depends on this needs to be clearly addressed with clear, precise experimental results.

To address this issue we included the co-‐labeling of Repo and Elav in (Figure 5—figure supplement 1A-B”), as requested. Concerning the possible indirect effects of glia to neuron conversion on migration, we find that only 14% of the gcm KD wings show conversion, whereas 75% of the KD wings show delayed migration. Further, to make sure that glia to neuron transformation in gcm knockdown wings does not affect glia migration we analyzed two sets of KD wings. The first set comprises of wings with fate transformation and the migratory index of these wings is very low, 25%. The second set includes wings without fate transformation; in this case, the migratory index is 45%, which is still significantly different from the control (Figure 5—figure supplement 1C). Last but not the least, the number of converted cells is very low (10). We also analyzed an allelic combination that is weaker than the KD and that only shows fate conversion in 10% of the wings. We still find a migratory delay in these wings (gcm‐Gal4/gcmrA87animals, See Figure 3A). This analysis confirms that Gcm affects glia migration irrespective of the change in the number of glial cells due to the fate transformation, which altogether fits perfectly with the finding that Gcm directly induces the expression of the Fra receptor.
[Editors' note: further revisions were requested prior to acceptance, as described below.]
[…]
Reviewer #3:
I read the revised paper and cannot support publication at this stage due to the poor demonstration that fra is expressed by the glia and the lack of a clear demonstration that the phenotype is due to loss of glial fra function. I admit there are some hints supporting this suggestion but I do not think the data are conclusive.
Concerning the detailed points described in, we addressed the three issues, as previously agreed:
1) Better evidence for Fra expression/requirement in glia
2) Migratory assay in fra4/gcm null wings
3) Mitotic clones
Expression of fra in glia. The authors added a new Figure 1H-Jwhich shows fra expression upon glial knockdown. For me this is not showing that fra is expressed in glia. There is as much background as in the wild type and the control experiment, neuronal knockdown, is not shown based on the argument that elav is also expressed by glia. Please, use nsybGal4 which is specific – and/or generate clones! The quality of the antibody stainings is not consistent and for example in Figure 2E glial knockdown is also removing the neuronal Fra expression. For me the quality of the data is still simply not good enough.
We followed the suggestion and used the nsybGal4 line as a neuronal driver. We first assessed expression of this driver in wing neurons, next we used the nsybGal4 line to knock down Fra in neurons and analyze glial migration. We found no defects in the migratory index, except in wings that carry an axonal navigation defect. These data are now presented in Figure 3—figure supplement 1. The wings carrying axon navigation defects were of course eliminated from the graph, as it has already been published that axonal defects result in glial migration defects.
We add Fra immunolabelling at 18 hAPF in combination with a second glial reporter (Repo>GFP). See Figure 1.
We removed panels H-J from Figure 1, as requested.
We improved the quality of Figure 2 showing glial Fra expression at later stages, as requested.
The phenotypes are based on somewhat odd genetic backgrounds. May be there is a fra haploinsufficiency and heterozygous fra/+ animals show a phenotype. But to really determine the Fra phenotype, a true loss of function analysis is required. The argument that MARCM clones are hard to get is weird as the same lab has generated these clones before in another study. Glia-specific fra knockdown animals show roughly the same phenotype as heterozygous fra3/+ animals (MI 58 vs.61). Still only the RNAi induced phenotype can be completely rescued. The authors conclude that this may be due to the reduction of neuronal Fra expression (which was not tested). However, one could also argue that the authors are looking at a genetic background effect. Thus, a MARCM analysis is mandatory.
We complied with the request of the referee and analyzed the migratory phenotype in transheterozygous animals fra3/gcm null (N7-4), which further confirms the data obtained with the hypomorphs and the KD line. These data are now presented in Figure 3—figure supplement 2, second last column.
We would like to clarify that the allelic series showing the reduction in the number of Repo+ cells in the last rebuttal letter shows the effects of gcmrA87/gcmN7-4and not gcm>GFP/ gcmN7-4, as claimed by the referee (see the same figure panel). We have already mentioned in the manuscript that gcmrA87 is a weak gcm hypomorphic allele. Therefore, the effect it produces on the number of repo cells when combined with null gcmN7-4 allele is similar to gcm>GFP/gcm>GFP.
We had asked for the effect of gcm on the migratory phenotype of fra3/+ heterozygous animals. Although the authors provide arguments for an allelic series of the different gcm alleles they do not show the effect of a true gcm loss of function allele on the migratory index. Provided that the presence of one very weak, homozygous viable, gcm allele has such a profound effect on the migratory index, what would a true null allele do? The allelic series shown in the rebuttal letter indicates that homozygous gcm>GFP animals have the same effect as gcm>GFP gcmN7-4. Since gcmN7-4 is a null this implies that gcm>GFP is a null as well – but why do they survive then? Something is odd here.
We performed mitotic clone analyses. We ordered the strain to produce MARCM clones in WT or in Fra KD cells. We also ordered strains to perform a more refined clonal analysis using the recently developed coinFLP technique in a wild type and in a Fra KD background. The crosses from this second approach went faster. Given the time frame we had agreed upon and the fact that the PhD student could only come for one month, we pursued the analyses of the coinFLP clones, which was indeed very successful. We could identify wild type and mutant clones in the same wing using two membrane reporters (GFP- green cells: WT cells; RFP – red cells: Fra DK cells). In this way, we could confirm that red (Fra KD) cells contain much less Fra than green (WT) cells. Furthermore, we could specifically identify WT and Fra KD glia using the Repo marker that is detected in glial nuclei. All these data are now shown in Figure 1.
https://doi.org/10.7554/eLife.15983.034Article and author information
Author details
Funding
CEFIPRA-4403-1 (graduate student fellowship)
- K VijayRaghavan
- Angela Giangrande
Agence Nationale de la Recherche (international award)
- Angela Giangrande
Fondation pour la Recherche Médicale (labelisation)
- Angela Giangrande
ARC Centre of Excellence for Coherent X-Ray Science (Projet grant)
- Angela Giangrande
Ligue Contre le Cancer (Grant regional)
- Angela Giangrande
USIAS
- Angela Giangrande
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Benjamin Altenhein, Thomas Kidd, Barry Dickson, Matthias Landgraf, Iris Salecker, Wesley B Grueber, DHSB, VDRC, Kyoto Stock Center and the Bloomington Stock Center for reagents and flies. We thank Claude Delaporte, Celine Diebold and the fly, cell separation and imaging facilities for technical assistance. We thank Yoshi Yuasa, Wael Bazzi as well as the other members of the lab for valuable input and comments on the manuscript. INSERM, CNRS, UDS, Hôpital de Strasbourg, ARC, INCA, Indo-French Center for the Promotion of Advanced Research (CEFIPRA), Fondation pour la Recherche Médicale en France (FRM) and ANR grants supported this work. T Gupta, A Kumar were funded by AFM, ARC, CEFIPRA and FRM fellowships, respectively, and PB Cattenoz is supported by the ANR. The IGBMC was also supported by a French state fund through the ANR labex.
Reviewing Editor
- Hugo J Bellen, Baylor College of Medicine, United States
Publication history
- Received: March 11, 2016
- Accepted: October 12, 2016
- Accepted Manuscript published: October 14, 2016 (version 1)
- Version of Record published: November 8, 2016 (version 2)
Copyright
© 2016, Gupta et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,076
- Page views
-
- 300
- Downloads
-
- 4
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Stem Cells and Regenerative Medicine
Cell-free DNA (cfDNA) present in the bloodstream or other bodily fluids holds potential as a non-invasive diagnostic for early disease detection. However, it remains unclear what cfDNA markers might be produced in response to specific tissue-level events. Organoid systems present a tractable and efficient method for screening cfDNA markers. However, research investigating the release of cfDNA from organoids is limited. Here, we present a scalable method for high-throughput screening of cfDNA from cardiac organoids. We demonstrate that cfDNA is recoverable from cardiac organoids, and that cfDNA release is highest early in differentiation. Intriguingly, we observed that the fraction of cell-free mitochondrial DNA appeared to decrease as the organoids developed, suggesting a possible signature of cardiac organoid maturation, or other cardiac growth-related tissue-level events. We also observe alterations in the prevalence of specific genomic regions in cardiac organoid-derived cfDNA at different timepoints during growth. In addition, we identify cfDNA markers that were increased upon addition of cardiotoxic drugs, prior to the onset of tissue demise. Together, these results indicate that cardiac organoids may be a useful system towards the identification of candidate predictive cfDNA markers of cardiac tissue development and demise.
-
- Developmental Biology
- Neuroscience
Fenestrated and blood-brain barrier (BBB)-forming endothelial cells constitute major brain capillaries, and this vascular heterogeneity is crucial for region-specific neural function and brain homeostasis. How these capillary types emerge in a brain region-specific manner and subsequently establish intra-brain vascular heterogeneity remains unclear. Here, we performed a comparative analysis of vascularization across the zebrafish choroid plexuses (CPs), circumventricular organs (CVOs), and retinal choroid, and show common angiogenic mechanisms critical for fenestrated brain capillary formation. We found that zebrafish deficient for Gpr124, Reck, or Wnt7aa exhibit severely impaired BBB angiogenesis without any apparent defect in fenestrated capillary formation in the CPs, CVOs, and retinal choroid. Conversely, genetic loss of various Vegf combinations caused significant disruptions in Wnt7/Gpr124/Reck signaling-independent vascularization of these organs. The phenotypic variation and specificity revealed heterogeneous endothelial requirements for Vegfs-dependent angiogenesis during CP and CVO vascularization, identifying unexpected interplay of Vegfc/d and Vegfa in this process. Mechanistically, expression analysis and paracrine activity-deficient vegfc mutant characterization suggest that endothelial cells and non-neuronal specialized cell types present in the CPs and CVOs are major sources of Vegfs responsible for regionally restricted angiogenic interplay. Thus, brain region-specific presentations and interplay of Vegfc/d and Vegfa control emergence of fenestrated capillaries, providing insight into the mechanisms driving intra-brain vascular heterogeneity and fenestrated vessel formation in other organs.