Fibroblast Growth Factor Signaling Instructs Ensheathing Glia Wrapping of Drosophila Olfactory Glomeruli ======================================================================================================== * Bing Wu * Jiefu Li * Ya-Hui Chou * David Luginbuhl * Liqun Luo ## Abstract The formation of complex yet highly organized neural circuits requires interactions between neurons and glia. During the assembly of the *Drosophila* olfactory circuit, 50 olfactory receptor neuron (ORN) classes and 50 projection neuron (PN) classes form synaptic connections in 50 glomerular compartments in the antennal lobe, each of which represents a discrete olfactory information processing channel. Each compartment is separated from the adjacent compartments by membranous processes from ensheathing glia. Here we show that Thisbe, a fibroblast growth factor (FGF) released from olfactory neurons, particularly local interneurons, instructs ensheathing glia to wrap each glomerulus. The Heartless FGF receptor acts cell-autonomously in ensheathing glia to regulate process extension so as to insulate each neuropil compartment. Overexpressing Thisbe in ORNs or PNs causes over-wrapping of glomeruli to which their axons or dendrites target. Failure to establish the FGF-dependent glia structure disrupts precise ORN axon targeting and discrete glomerular formation. **Significance Statement** This research reports that reciprocal interactions between *Drosophila* olfactory neurons and ensheathing glia mediate the formation of neuronal compartments—groups of synapses that are packed into discrete structures called glomeruli that carry specific olfactory information. Ensheathing glia respond to a neuronal cue, the fibroblast growth factor (FGF) Thisbe, to pattern the boundaries of the nascent compartments. Neural compartments in turn require such glial barriers to separate themselves from neighboring compartments, so as to ensure the correct organization of the olfactory circuit. These findings highlight the importance of glia in the assembly and maintenance of neural circuits and the functions of FGF signaling in these processes. ## Introduction Glia and neurons interact dynamically to coordinate the development and function of neural circuits. For example, *Drosophila* midline glia serve as guideposts that provide spatial cues to direct axonal pathfinding in the embryonic ventral nerve cord (1, 2). Neurons in turn provide signals that weave glia into the fabric of the nervous system. For instance, axonal Neuregulin-1 levels dictate the extent of Schwann cell myelination (3). A myriad of other biological processes in the developing and adult brain, including synapse formation (4-7), elimination (8, 9), and regulation of synaptic activity (10-13), all require extensive communication and cooperation between neurons and glia. The *Drosophila* nervous system contains a diverse array of glial cell types (14), offering opportunities to discover new mechanisms for glia-neuron interactions. Ensheathing glia are a group of cells that insulate neighboring neuropil structures by extending their membranous processes along the outer surface of synaptic neuropil or their sub-compartments, without invading the inner part of the neuropil (15). It remains unknown how ensheathing glia establish such a barrier-like structure, what molecular signals orchestrate this process, and whether the glial barrier is essential for the integrity of the encircled neuropil. The *Drosophila* antennal lobe provides an excellent experimental system for tackling these questions with high resolution. The antennal lobe is organized into 50 discrete neuropil compartments—the glomeruli wherein specific types of olfactory receptor neuron (ORN) axons and projection neuron (PN) dendrites form synapses—each of which is surrounded by ensheathing glia processes (15-17). In addition, the specific and stereotypic projection pattern of neurons in the antennal lobe (18, 19) renders the system convenient for exploring the potential neuronal disorganization caused by malformation of glial structures. Previous studies have utilized this system to find that in adults, ensheathing glia increase their process extension to injured axons (20), help clear degenerating axons (20), and strengthen excitatory interactions between surviving neurons (21). Here we study ensheathing glia in the antennal lobe formation and organization during development. ## Results ### Pupal Development of Antennal Lobe Ensheathing Glia To study the development and function of antennal lobe ensheathing glia, we searched for genetic tools that can specifically label these cells during the pupal stage when the olfactory circuit in the antennal lobe is being assembled. We preselected a number of enhancer-GAL4 lines from the FlyLight GAL4 collection based on the published expression pattern in adult animals (22), and a few *MiMICGAL4* lines (23). We tested whether these GAL4 lines were active at different pupal stages. We found that at 96 hours after puparium formation (h APF), a late stage in antennal lobe morphogenesis, *SPARC-GAL4* (Fig 1D), *GMR56F03-GAL4* (Fig S1) (24), and *GMR10E12-GAL4* (Fig S1) flies expressed GAL4 predominantly in ensheathing glia. As indicated by the *UAS-mCD8GFP* reporter that labeled the cell membranes and processes in the presence of GAL4, these GAL4+ cells were located on the periphery of the antennal lobe and extended their processes to wrap around the antennal lobe as well as individual glomeruli within it, but with minimal invasion into each glomerulus. These are characteristic morphological features of ensheathing glia. We also used *UAS-nuclear-LacZ* to mark the nuclei of these GAL4+ cells (Fig 1E), and found that every LacZ+ nucleus around the antennal lobe (within 10-μm distance to the surface of the antennal lobe) was also positive for Repo, a glia marker in *Drosophila* (25), confirming that these GAL4+ cells are indeed glia. On average ~100 Repo+ cells were detected around the antennal lobe, and ~56 were LacZ+ (Fig 1F). Repo+/LacZ– cells are likely astrocytes and cortex glia present in the vicinity of ensheathing glia. To further distinguish these GAL4+ glia from astrocytes, we used the GAT antibody (26) to mark astrocytes, and found that the majority (~95%) of the GAL4+ glia are negative for GAT (Fig S1). In summary, *SPARC-GAL4*, *GRM56F03-GAL4*, and *GMR10E12-GAL4* mark ensheathing glia around the antennal lobe in the late pupal stage. ![Fig. 1.](http://biorxiv.org/https://www.biorxiv.org/content/biorxiv/early/2017/06/18/151662/F1.medium.gif) [Fig. 1.](http://biorxiv.org/content/early/2017/06/18/151662/F1) Fig. 1. Ensheathing glia morphogenesis during antennal lobe development. *(A–D)* Confocal sections for antennal lobe at 24 h (A), 48 h (B), 72 h (C), and 96 h (D) APF. Ensheathing glia processes are labeled by *SPARC-GAL4* driven (>) *UAS-mCD8GFP*. Ensheathing glia nuclei are marked by *UAS-nuclear-LacZ*. Magenta shows neuropil staining by the N-cadherin (Ncad) antibody. Dashed rectangles are enlarged on the right column. D, dorsal; L, lateral. Scale bars are 10 μm. *(E)* Confocal sections of 96 h antennal lobe. Green shows glia nuclei stained by the Repo antibody, overlaid on the LacZ channel in *E2*. *(F)* Quantification of antennal lobe ensheathing glia cell number at different developmental stages, and the number of total glia around the antennal lobe at the end of pupae stage. Error bars represent standard deviation (SD). N=5 for each time point. *SPARC-GAL4* also drove reporter expression at 24 h APF, which enabled us to investigate the development of antennal lobe ensheathing glia. At 24 h APF, prior to glomerular formation, ensheathing glia processes were restricted to the periphery of the antennal lobe (Fig 1A). From 24 h to 48 h APF, axon terminals of ORNs invade the antennal lobe and target specific subregions to connect with their synaptic partners such as dendrites from PNs (17). We found that around 48h APF when proto-glomeruli first emerged, ensheathing glia processes started to invade the antennal lobe from the lobe surface (Fig 1B). Even at this initial stage of infiltration, ensheathing glia processes displayed a preference to grow along the borders between adjacent proto-glomeruli, instead of growing into them (Fig 1B3). At 72 h APF, as glomerular structures became more clearly separable, ensheathing glia processes had encircled most glomeruli (Fig 1C). These processes remained on glomerular boundaries, with minimal extension into the glomeruli. Ensheathing glia wrapping of olfactory glomeruli was complete by the end of the pupal stage (Fig 1D). We also observed an increase in the number of *SPARC-GAL4*+ ensheathing glia around the antennal lobe from 24 h to 72 h APF (Fig 1F). This increase is at least in part due to glial cell proliferation, since we could generate glia clones via the MARCM technique (27) with heat shock induced mitotic recombination at any point from 24 to 72 h APF (see Fig 3 and Fig S3 below). ### *Heartless* Knockdown in Ensheathing Glia Reduces Processes and Disrupts the Ensheathment Pattern To identify the molecular signals that control the extension of ensheathing glia processes, we tested ~50 candidate adhesion molecules, receptor kinases, and receptor phosphatases by RNA interference (RNAi) and screened for potential defects in ensheathing glia development and antennal lobe organization. We found that expression of two independent RNAi targeting non-overlapping regions of *heartless* (*htl*), a fibroblast growth factor receptor (28-30), caused a 70% reduction of ensheathing glia processes within the antennal lobe (Fig 2A-D). This phenotype is consistent with a previous discovery that FGF signaling is required for astrocytes to infiltrate the larval ventral nerve cord of *Drosophila* (26). However, the antennal lobe provides a unique opportunity to observe the response of neuronal compartmentalization to glia morphogenesis defects. ![Fig. 2.](http://biorxiv.org/https://www.biorxiv.org/content/biorxiv/early/2017/06/18/151662/F2.medium.gif) [Fig. 2.](http://biorxiv.org/content/early/2017/06/18/151662/F2) Fig. 2. Heartless knockdown causes an altered ensheathing glia wrapping pattern and disrupts antennal lobe compartmentalization. *(A–C)* Confocal sections for the antennal lobe at 96 h APF of wild type (A) and two different *UAS-Htl RNAi* constructs (RNAi#1: VDRC6692, RNAi#2: VDRC27180) driven by *SPARC-GAL4 (B, C)*. Ensheathing glia processes are labeled by *SPARC-GAL4> UAS-mCD8GFP*. Neuropil compartments are stained with Ncad antibody. GFP and Ncad intensities along a line between the center of DM6 and DA1 glomeruli (indicated by dashed circles in *A2* and asterisks in *A2–C2*) are plotted in *A3–C3*, with the X-axis indicating the position of each point on the line from DM6 to DA1 (left to right); fluorescence intensity quantified in arbitrary unit. *(D)* Quantification of total GFP intensities normalized by antennal lobe size for wild-type and RNAi expressing flies. *(E)* Quantification of relative standard deviation of Ncad intensities for each antennal lobe of wild-type and RNAi flies. (F) Correlation coefficient of GFP and Ncad intensities as plotted in *A3–C*3. Scale bar is 10 μm. Error bars represent SD. **, p<0.01; \***|, p<0.001, \**\*|\*, p<0.0001. We found that many glomeruli became less clearly separable (Fig 2A2, B2, C2, arrowheads) when *htl* was knocked down. To quantify this effect, we used relative standard deviation (RSD) of the N-cadherin neuropil staining intensity as an index for the degree of compartmentalization, as incomplete antennal lobe compartmentalization would cause obscure glomerular borders and hence smaller variance in the neuropil staining signal. We found a significant reduction of RSD in *htl* knockdown antennal lobes compared with wild type (Fig 2E). The remaining glomerular borders sometimes also lacked ensheathing glia processes that normally would divide the adjacent glomeruli (Fig 2A-C). We also observed ensheathing glia processes that extended within glomeruli, which is rarely observed in wild type antennal lobes (Fig 2B1, arrowhead). This is likely due to poor establishment of glomerular borders. We quantified the localization pattern of the ensheathing glia processes relative to the glomerular compartments by plotting the intensity of the GFP signal derived from glial membranes together with the intensity of Ncadherin signal marking the neuropil (Fig 2A3-C3), and calculated their correlation coefficient (Fig 2F). In control animals, N-cadherin and GFP signals exhibited strong anti-correlation (Fig 2F), since ensheathing glia processes are preferentially located on the borders of neuropil compartments instead of within the glomeruli. However, such anti-correlation was greatly diminished when *htl* was knocked down (Fig 2F), suggesting that loss of *htl* results in a disrupted ensheathment pattern. ### *Heartless* Promotes Ensheathing Glia Survival and Cell-autonomously Regulates Process Elaboration Heartless is involved in cell differentiation, directional migration, and survival in a wide variety of tissues (26, 28, 31-34). Consistent with this notion, we observed that *htl* knockdown caused a 40% reduction of ensheathing glia number near the antennal lobe as assayed by the nuclear-LacZ marker (Fig 3A-D). This decrease is likely due to cell death, since the expression of an apoptosis suppresser, P35 (35), in ensheathing glia largely rescued the reduction of cell number (Fig S2); however, P35 expression did not rescue the reduction of ensheathing glia processes or the defect in glomerular compartmentalization (Fig S2), suggesting that FGF signaling may also regulate morphogenesis of antennal lobe ensheathing glia in addition to its role in supporting glia survival. ![Fig. 3.](http://biorxiv.org/https://www.biorxiv.org/content/biorxiv/early/2017/06/18/151662/F3.medium.gif) [Fig. 3.](http://biorxiv.org/content/early/2017/06/18/151662/F3) Fig. 3. Heartless controls ensheathing glia number and cell-autonomously control their morphology. *(A–C)* Projections of confocal sections along the Z-axis. Red: ensheathing glia cell bodies marked by *SPARC-GAL4> UAS-nuclear-LacZ.* Blue: Ncad antibody stains for antennal lobe neuropil. *(D)* Quantification of ensheathing glia cell number in each antennal lobe. \**\*|\*, p<0.0001. *(E, F)* Wild-type *(E)* and *htlAb42/Ab42 (F)* ensheathing glia single cell MARCM clones in adult antennal lobe. Green shows ensheathing glia processes labeled by *GMR10E12-GAL4> UAS-mCD8GFP*. Blue shows ensheathing glia nuclei marked by *UAS-nuclear-LacZ*. Magenta shows Ncad antibody staining for antennal lobe neuropil. *(G)* Quantification of the volume of processes from each ensheathing glia labeled by MARCM. *(H)* Quantification of the number of borders each MARCM-labeled ensheathing glia contacts. Error bars represent SD. **, p<0.01; \***|, p<0.001. (I–K) Confocal sections of antennal lobe at 24 h *(I)*, 48 h *(J)*, and 72 h *(K)* APF with Htl-GFP signal. Dashed rectangle in *J1* is enlarged and shown as *J3–J5*. The yellow stars mark the center of a proto-glomerulus around which ensheathing glia is wrapping. Scale bars are 10 μm. To dissociate the role of *htl* in controlling glia survival and process extension, we used the MARCM technique (27) to generate single-cell clones of ensheathing glia homozygous for the *htlAb42(null)* (30) allele in an otherwise heterozygous background. In this experiment, *GMR10E12-GAL4* was used to label the homozygous mutant cells. We used *UAS-mCD8GFP* to visualize glial processes (Fig 3E-F), and *UAS-nuclear-LacZ* (Fig 3E1, F2, arrows) to visualize cell bodies of ensheathing glia in these MARCM clones. Compared to wild type, *htlAb42* single-cell clones exhibited a 50% reduction of the volume of glial processes (Fig 3G) and a 50% reduction of total fluorescence intensity (Fig S3). These reductions were accompanied by a decrease in the number of glomeruli each ensheathing glia could access, as quantified by the number of glomerular borders each glia extended to (Fig 3H). Thus, this mosaic experiment demonstrated that *htl* is cell-autonomously required for ensheathing glia process extension. To determine the subcellular localization of the Htl protein, we used a fosmid transgenic line that produces Htl-GFP from an insertion that contains the extended genomic region covering *htl* (36) to analyze GFP signal within and around the antennal lobe during pupal development. At 24 h APF, no Htl-GFP signal was detected inside the antennal lobe (Fig 3I), consistent with the location of ensheathing glia processes (Fig 1A). At 48 h APF, we started to detect Htl-GFP signal in between proto-glomeruli (Fig 3J). By 72 h APF, Htl-GFP signal was detected at most glomerular borders, though a lower level of signal was also detected within glomeruli (Fig 3K), which may originate from cell types other than ensheathing glia, such as astrocytes and potentially neurons. The signal on glomerular borders coincided well with the ensheathing glia wrapping pattern over this developmental course (Fig. 1). These data support the hypothesis that Htl mediates FGF signaling to regulate ensheathing glia process invasion of the antennal lobe and wrapping of individual glomeruli. ### Thisbe, an FGF ligand, Is Required for Ensheathing Glia Wrapping of Glomeruli We next assessed the identity and cellular source of FGF ligands that regulate ensheathing glia development. Htl responds to two FGF ligands, Thisbe (Ths) and Pyramus (Pyr) (37, 38). To test whether *ths* is required for antennal lobe ensheathing glia to form the wrapping pattern, we analyzed ensheathing glia and antennal lobe morphology in animals that carried different *ths* mutant alleles (39). The phenotypes we observed with loss of *htl* were all recapitulated in trans-heterozygous combinations of *ths* alleles (*ths*e02026/*ths*759 and *Df(2R)ED2238*/*ths*759). There was an overall reduction of ensheathing glia processes and cell number in the antennal lobe of trans-heterozygous flies (Fig 4A-B, E-F). Similar to loss of *htl* from ensheathing glia, the antennal lobe of *ths* mutants also showed defective compartmentalization as measured by reduced RSD values of N-cadherin staining (Fig 4G) and a markedly reduced anti-correlation between the intensity of ensheathing glia processes and neuropil signal (Fig 4H). These findings indicate that *ths* is required for ensheathing glia to wrap around glomeruli in the antennal lobe, and that failure to form a correct wrapping pattern could disrupt compartmentalization of antennal lobe neuropil. ![Fig. 4.](http://biorxiv.org/https://www.biorxiv.org/content/biorxiv/early/2017/06/18/151662/F4.medium.gif) [Fig. 4.](http://biorxiv.org/content/early/2017/06/18/151662/F4) Fig. 4. Thisbe is expressed in olfactory neurons and required for ensheathing glia to wrap glomeruli. *(A–D)* Ensheathing glia wrapping pattern in wild type *(A)*, ths mutant *(B; thse02* denotes *thse02026* allele), panneural C155-GAL4 (C) and LN *orb0449-GAL4* (D) driven RNAi against ths. Green shows ensheathing glia processes labeled by *SPARC-QF> QUAS-mtdT*. Magenta shows ensheathing glia nuclei marked with *QUAS-nuclear-LacZ*. *(E–H)* Quantification of ensheathing glia number *(E)*, process intensity (F), relative standard deviation of Ncad staining *(G)*, and the correlation coefficient for the intensities of ensheathing glia process and neuropil staining *(H)*. Error bars represent SD. *, p<0.05; **, p<0.01; \***|, p<0.001; \**\*|\*, p<0.0001; ns, p>0.05. *(I) ths* expression pattern revealed by ths-GAL4> UAS-mCD8GFP (green). Magenta in I1 shows Ncad counterstaining. *(J) ey-FLP* intersects with ths-GAL4 together with *UAS-FRT-stop-FRT-mCD8GFP* to show ths expression pattern in ORNs. *(K) GH146-FLP* intersection shows ths expression pattern in PNs. *(L) ths-GAL4+* LN cell bodies and LN and ORN processes, after *GH146-GAL80* suppressing *ths-GAL4* in most PNs. *(M)* MARCM labels a single LN that is positive for *ths-GAL4*. Green: GFP. Magenta: Ncad. All images are confocal sections of adult antennal lobe except panel (M), which is a projection of Z stacks. Scale bars are 10 μm. ### Thisbe Is Produced by ORNs, PNs, and Antennal Lobe Local Interneurons (LNs) To identify which cell types express Thisbe, we took advantage of a MiMIC insertion (23) located between two coding exons of the *ths* gene. We converted the MiMIC cassette to an artificial exon that contains the coding sequence for *2A-GAL4* (40). Thus, GAL4 can be produced along with the endogenous N-terminus of Ths (encoded by the first two exons), and can drive reporter gene expression in the *ths* pattern. We then used this *Ths-GAL4* to express *UAS-mCD8GFP* to visualize the cell bodies and projections of the Ths-producing cells. In the antennal lobe, neuronal processes from PNs, ORNs, and LNs are all labeled, suggesting that all three major neuronal types produce Ths (Fig 4I). To determine the contributions of each of these cell types, we used an intersectional strategy, in which *Ths-GAL4* was combined with Flp recombinases that are specifically expressed in ORNs (*ey-FLP*) (41) or PNs (*GH146-FLP*) (42). With a FLP-out reporter, *UAS-FRT-stop-FRTmCD8GFP* (42), we found a number of ORN and PN classes were *Ths-GAL4*+ based on glomerular labeling (Fig 4J-K). At least 45 cell bodies over all sections of the antennal lobe were *Ths-GAL4+* without intersection (Fig 4I); however, only 19 of them were *GH146-Flp*+ PNs (which constitutes the majority of PNs innervating ~40 glomeruli), consistent with the small number of glomeruli (~7) labeled by intersection of *GH146-Flp* and *Ths-GAL4* (Fig 4K). This difference is likely due to the contribution from LNs, whose cell bodies are also located around the antennal lobe (43, 44), whereas ORN cell bodies are located in the peripheral sensory organs. We used two approaches to validate that LNs produce Thisbe. First, we suppressed the expression from most PNs by combining *GH146-GAL80* (42) together with *Ths-GAL4*. After PN-derived signal was largely eliminated, around 27 cell bodies around the antennal lobe remained. The projection pattern of these neurons covered the entire antennal lobe (Fig 4L), which is characteristic of the majority of LNs (43). Second, we used *Ths-GAL4* to label single cells by the MARCM technique. We were able to identify MARCM clones for ORNs and PNs (data not shown), as well as LNs (Fig 4M). In summary, Ths is produced by a subset of ORNs, PNs, and LNs. ### LN-derived Thisbe Is Necessary for Ensheathing Glia Wrapping To test in which cell type(s) Thisbe functions to regulate ensheathing glia morphogenesis and antennal lobe compartmentalization, we used RNAi to knock down *ths* in ORNs (*Pebbled-GAL4*), PNs (*GH146-GAL4*), LNs, and all neurons (*C155-GAL4*), respectively, while labeling the ensheathing glia by *SPARC-QF*, which we converted from *SPARC-GAL4* (23). Pan-neuronal knockdown of *ths* recapitulated the phenotypes observed in *ths* mutants (Fig 4C, E-H). Knocking down *ths* specifically in LNs by *orb0449-GAL4* (Fig. S4), a GAL4 line identified from the InSite screen (45), resulted in a mild but significant defect in ensheathing glia wrapping and antennal lobe glomerulus integrity (Fig 4D, E-H). Knocking down *ths* in PNs did not cause significant defects (Fig S5). We have also used MARCM combined with a cell lethal strategy (41) to create a near pan-ORN mutant background for *ths*, and did not find defects in ensheathing glia wrapping (Fig S5). The milder phenotypes in LN knockdown compared to pan-neuronal knockdown could be explained by: 1) pan-neuronal GAL4 may have stronger and/or earlier expression than the LNGAL4, and therefore causes more effective knockdown; 2) LN-GAL4 does not include all LNs; 3) Ths from ORNs and PNs synergize with Ths from LNs. Regardless, our data indicate that LN is an essential cellular source of Ths in the antennal lobe for directing ensheathing glia wrapping. ### Thisbe Can Instruct Glomerular Wrapping with a High Spatial Specificity We have shown that FGF signaling is necessary for ensheathing glia to extend processes into the antennal lobe and demarcate individual glomeruli. To test whether FGF signaling can instruct ensheathing glia to wrap around selected neuropil compartments, we expressed *ths* in only one class of ORNs (VA1v) using a GAL4 under the control of the promoter of the odorant receptor specifically expressed in this class (*Or47b-GAL4*). We observed that overexpression resulted in hyper-wrapped VA1v glomerulus by ensheathing glia (Fig 5B, D). This effect was highly localized, since in addition to VA1v, only the adjacent VA1d glomerulus was slightly hyper-wrapped (Fig 5D), likely due to intensified ensheathing glia processes on its border shared with VA1v. Hyper-wrapping did not extend to the DA1 glomerulus (Fig 5D), which is one glomerulus away from the VA1v glomerulus. There was also an excess of glia cells around the hyper-wrapped VA1v glomerulus (Fig 5F). ![Fig. 5.](http://biorxiv.org/https://www.biorxiv.org/content/biorxiv/early/2017/06/18/151662/F5.medium.gif) [Fig. 5.](http://biorxiv.org/content/early/2017/06/18/151662/F5) Fig. 5. Changes in ensheathing glia processes and number in response to local overexpression of ths. *(A–C)* Confocal sections of adult antennal lobe of wild type *(A)*, overexpression of *UAS-ths* by *Or47b-GAL4 (B)* or by *Mz19-GAL4 (C)*. Green: *SPARC-QF>QUAS-mtdT* labels ensheathing glia processes. Magenta: Ensheathing glia nuclei marked with *QUAS-nuclear-LacZ*. Scale bar is 10 μm. *(D, E)* Quantification of signal intensities of ensheathing glia processes around each glomerulus normalized by the perimeter of the corresponding glomerulus. *(F, G)* Quantification of ensheathing glia numbers within 5-μm distance to the surface of each glomerulus. Error bars represent SD. *, p<0.05; **, p<0.01; \***|, p<0.001; \**\*|\*, p<0.0001; ns, p>0.05. Similarly, overexpressing Ths in *Mz19-GAL4*+ PNs, which send dendrites to DA1 and VA1d, also caused local hyper-wrapping of these glomeruli, as well as a local increase of ensheathing glia cells (Fig 5C, E, G). We have consistently noticed that *Mz19-GAL4* activity is stronger in DA1 PNs than in VA1d PNs (46). Accordingly, glial hyper-wrapping was more pronounced around DA1 than around VA1d. These results suggest that Thisbe acts locally as a spatial cue to instruct ensheathing glia to infiltrate the antennal lobe. ### FGF Signaling Ensures Accurate Neuronal Targeting Lastly, we examined whether glia wrapping defects could disrupt neuronal projections to the glomerular compartments. We used membrane proteins under the control of specific odorant receptor promoters (*Or88a-mtdT* and *Or47b-rCD2*) to label two ORN classes that project their axons to two adjacent glomeruli, VA1v and VA1d (Fig 6A). When *ths* was knocked down with a pan-neuronal driver, *C155-GAL4*, VA1v and VA1d ORN axons still reached their target area in the ventrolateral antennal lobe. However, they failed to establish exclusive territories for their axonal arborization, and instead had partially overlapped axonal terminals (Fig 6B). Likewise, in *ths* mutant animals, VA1v and VA1d axon terminals partially intermingled, rather than confining to their respective glomerular compartments (Fig 6C; quantified in Fig 6D). Thus, FGF signaling between neurons and ensheathing glia is essential for antennal lobe compartmentalization and targeting accuracy for ORN axons. ![Fig. 6.](http://biorxiv.org/https://www.biorxiv.org/content/biorxiv/early/2017/06/18/151662/F6.medium.gif) [Fig. 6.](http://biorxiv.org/content/early/2017/06/18/151662/F6) Fig. 6. Loss of ths results in olfactory receptor neuron targeting defects. *(A–D)* Confocal sections of the adult antennal lobe, with *Or88a-mtdT* and *Or47b-rCD2* labeling ORN axons that target VA1d (magenta) and VA1v (green) glomerulus, respectively. Images within the dashed rectangles in the first column are enlarged in columns 2–4. *(A)* Wild type. *(B) C155-GAL4* driven *ths RNAi*. *(C-D)* ths mutant flies showed severe intermingling of VA1d and VA1v ORN terminals (C) and a mild spillover of VA1d ORN axons (D, arrow). *ths238* denotes the *Df(2R)ED2238* allele. Scale bars are 10 μm. *(E)* Quantification of the accuracy of ORN axon targeting. 1, normal; 2, mild spillover; 3, severe intermingling. Each symbol represents 1 fly with two antennal lobes scored separately and averaged. Error bars represent SD. *, p<0.05; **, p<0.01. ## Discussion The use of discrete neuropil compartments for organizing and signaling information is widespread in invertebrate and vertebrate nervous systems. In both the fly antennal lobe and vertebrate olfactory bulb, axons from different ORN classes are segregated into distinct glomeruli (47). The rodent barrel cortex also uses discrete compartments, the barrels, to represent individual whiskers (48). In this study, we show that FGF signaling between neurons and glia mediate neural compartment formation in the *Drosophila* antennal lobe. Members of the FGF family have diverse functions in a variety of tissues in both vertebrates and invertebrates (49, 50). Vertebrate FGFs regulate not only neural proliferation, differentiation, axon guidance, and synaptogenesis, but also gliogenesis, glial migration, and morphogenesis (51-55). Many of these roles are conserved in invertebrates. For example, Thisbe and Pyramus induce glial wrapping of axonal tracts (32, 33) much like the role of other FGF members played in regulating myelin sheaths in mammals (55). Thisbe and Pyramus also control *Drosophila* astrocyte migration and morphogenesis (26). Likewise, FGF signaling also promotes morphogenesis of mammalian astrocytes (56). Therefore, studying the signaling pathways in *Drosophila* will extend our understanding of the principles of neural development. In ensheathing glia, whose developmental time course and mechanisms have not been well documented prior to this study, we observed a glial response to FGF signaling reminiscent of the paradigm shown before (26, 32, 57); however, the exquisite compartmental structure of the *Drosophila* antennal lobe and genetic access allowed us to further scrutinize changes of neuropil structure and projection patterns that occurred alongside morphological phenotypes in ensheathing glia. We demonstrated a necessity of Ths in local interneurons, although it is possible that ORNs and PNs also contribute. We also tested the function of the other ligand, Pyramus, in antennal lobe development. We did not detect any change in ensheathing glia morphology with *pyramus* RNAi, and double RNAi against *thisbe* and *pyramus* did not enhance the phenotype compared to *thisbe* knockdown alone (data not shown). FGF signaling in glomerular wrapping appears to be highly local. In our overexpression experiments, the hyper-wrapping effect was restricted to the glomerulus where the ligand is excessively produced, and did not spread to nearby non-adjacent glomeruli. These experiments suggest that Thisbe communicates locally to instruct glial ensheathment of the glomeruli, rather than diffusing across several microns to affect nearby glomeruli. Since heparan sulfate proteoglycans are known to act as FGF co-receptors by modulating the activity and spatial distribution of the ligands (50, 58, 59), we speculate that Thisbe in the antennal lobe may be subject to such regulation to limit its diffusion and long-range effect. Our data showed that deficient ensheathment of antennal lobe glomeruli is accompanied by imprecise ORN axon targeting. However, we cannot determine whether these targeting defects reflect initial axon targeting errors, or a failure to stabilize or maintain the discrete targeting pattern. Previous models for the establishment of antennal lobe wiring specificity suggested that the glomerular map is discernable by the time glia processes start to infiltrate the antennal lobe (17). Due to a lack of class-specific ORN markers for early developmental stages, it is still unclear the relative timing between when neighboring ORN classes refine their axonal targeting to discrete compartments and when ensheathing glia barriers are set up. Regardless, our discovery that FGF signaling functions in the formation of discrete neuronal compartments in the antennal lobe highlights an essential role for glia in the precise assembly of neural circuits. ## Methods ### Immunostaining Tissue dissection and immunostaining were performed according to previously described methods (60). Primary antibodies used in this study include rat anti-DNcad (DN-Ex #8; 1:40; DSHB), chicken anti-GFP (1:1000; Aves Labs), rabbit anti-DsRed (1:500; Clontech), mouse anti-rCD2 (OX-34; 1:200; AbD Serotec), rat anti-HA (1μg/ml, Roche), mouse nc82 (1:35; Developmental Studies Hybridoma Bank, [DSHB]), mouse anti-Repo (1:50, DSHB), rabbit anti-β-Galactosidase (MP Biomedicals, 1:125), rabbit anti-GAT (1:3000; a gift from Marc Freeman), and mouse anti-β-Galactosidase (1:1000; Promega). Secondary antibodies were raised in goat or donkey against rabbit, mouse, rat, and chicken antisera (Jackson Immunoresearch), conjugated to Alexa 405, FITC, 568, or 647. Confocal images were collected with a Zeiss LSM 780 and processed with Zen software, ImageJ, and Imaris. ### Mosaic Analysis The hsFlp MARCM analyses were performed as previously described (27, 61) with slight modifications. *GMR10E12-GAL4* was used for labeling ensheathing glia in adult stage. Flies were kept at 18°C and were heat shocked for 30min at 37°C between 0-24h APF. ### Data Analysis Confocal sections of antennal lobes (Fig. 2 and 4) were analyzed by ImageJ software for measuring integrated intensity value, area size, mean, and standard deviation for the region of interest (ROI) manually selected based on N-cadherin counterstaining. Process intensity was calculated using the integrated intensity value of the antennal lobe section crossing DA1 and DM6 glomeruli normalized by the area size of the antennal lobe of that section. The Plot Profile function in ImageJ was used to measure signal intensities for ensheathing glia processes and Ncad along the selected lines. Glia cell number was manually counted based on LacZ staining within 10-μm distance to the surface of the antennal lobe determined by Ncad counterstaining. To quantify the MARCM clones (Fig. 3), confocal sections of marked ensheathing glia were processed by Imaris software through manually thresholding the images by a set of consistent parameters, followed by automatic measurement of the volume and total intensity of glia processes. In the overexpression experiment (Fig. 5), ROI was manually selected by including the inner region and the boundaries of the glomerulus of interest based on Ncad counterstaining. Process intensity was defined as the integrated intensity value divided by the perimeter of the ROI selection. Glia number was manually counted based on LacZ signal within 5-μm distance to the surface of each glomerulus. ORN targeting defects (Fig. 6) were scored by an experimenter blind to the genotype. Graphs were generated using the GraphPad Prism software. Mean ± SD were shown on the graphs. Statistical significance was calculated with Graphpad Prism using two-tailed Student’s t test (Fig. 2-5) or Mann-Whitney test (nonparametric data in Fig. 6). ## Acknowledgments We thank M. Freeman, C. Klämbt, H. Bellen, T. Clandinin, and Vienna and Bloomington Stock Centers for reagents, X. Wang, J. Lui, A. Shuster, J. Ren, H. Li, T. Li, and L. DeNardo for discussion and helpful comments on the manuscript. This work was supported by an NIH grant (R01 DC005982) to L.L., who is an HHMI investigator. * Received June 18, 2017. * Revision received June 18, 2017. * Accepted June 18, 2017. * © 2017, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Mitchell KJ, et al. (1996) Genetic analysis of Netrin genes in Drosophila: Netrins guide CNS commissural axons and peripheral motor axons. Neuron 17(2):203-215. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/S0896-6273(00)80153-1&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=8780645&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=A1996VE84000004&link_type=ISI) 2. 2.Kidd T, Bland KS, & Goodman CS (1999) Slit is the midline repellent for the robo receptor in Drosophila. Cell 96(6):785-794. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(00)80589-9&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=10102267&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000079300100006&link_type=ISI) 3. 3.Michailov GV, et al. (2004) Axonal neuregulin-1 regulates myelin sheath thickness. Science 304(5671):700-703. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMDQvNTY3MS83MDAiO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 4. 4.Pfrieger FW & Barres BA (1997) Synaptic efficacy enhanced by glial cells in vitro. Science 277(5332):1684-1687. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIyNzcvNTMzMi8xNjg0IjtzOjQ6ImF0b20iO3M6Mzc6Ii9iaW9yeGl2L2Vhcmx5LzIwMTcvMDYvMTgvMTUxNjYyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 5. 5.Ullian EM, Sapperstein SK, Christopherson KS, & Barres BA (2001) Control of synapse number by glia. Science 291(5504):657-661. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIyOTEvNTUwNC82NTciO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 6. 6.Mauch DH, et al. (2001) CNS synaptogenesis promoted by glia-derived cholesterol. Science 294(5545):1354-1357. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIyOTQvNTU0NS8xMzU0IjtzOjQ6ImF0b20iO3M6Mzc6Ii9iaW9yeGl2L2Vhcmx5LzIwMTcvMDYvMTgvMTUxNjYyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 7. 7.Christopherson KS, et al. (2005) Thrombospondins are astrocyte-secreted proteins that promote CNS synaptogenesis. Cell 120(3):421-433. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2004.12.020&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=15707899&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000227028900016&link_type=ISI) 8. 8.Paolicelli RC, et al. (2011) Synaptic pruning by microglia is necessary for normal brain development. Science 333(6048):1456-1458. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMzMvNjA0OC8xNDU2IjtzOjQ6ImF0b20iO3M6Mzc6Ii9iaW9yeGl2L2Vhcmx5LzIwMTcvMDYvMTgvMTUxNjYyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 9. 9.Schafer DP, et al. (2012) Microglia sculpt postnatal neural circuits in an activity and complement-dependent manner. Neuron 74(4):691-705. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2012.03.026&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=22632727&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000304747200013&link_type=ISI) 10. 10.Bezzi P, et al. (1998) Prostaglandins stimulate calcium-dependent glutamate release in astrocytes. Nature 391(6664):281-285. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/34651&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=9440691&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000071484400050&link_type=ISI) 11. 11.Araque A, Parpura V, Sanzgiri RP, & Haydon PG (1998) Glutamate-dependent astrocyte modulation of synaptic transmission between cultured hippocampal neurons. The European journal of neuroscience 10(6):2129-2142. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1046/j.1460-9568.1998.00221.x&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=9753099&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000074009000020&link_type=ISI) 12. 12.Araque A, Sanzgiri RP, Parpura V, & Haydon PG (1998) Calcium elevation in astrocytes causes an NMDA receptor-dependent increase in the frequency of miniature synaptic currents in cultured hippocampal neurons. The Journal of neuroscience : the official journal of the Society for Neuroscience 18(17):6822-6829. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIxOC8xNy82ODIyIjtzOjQ6ImF0b20iO3M6Mzc6Ii9iaW9yeGl2L2Vhcmx5LzIwMTcvMDYvMTgvMTUxNjYyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 13. 13.Araque A, et al. (2014) Gliotransmitters travel in time and space. Neuron 81(4):728-739. 14. 14.Stork T, Bernardos R, & Freeman MR (2012) Analysis of glial cell development and function in Drosophila. Cold Spring Harbor protocols 2012(1):1-17. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1101/pdb.top067587&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=22194269&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 15. 15.Awasaki T, Lai SL, Ito K, & Lee T (2008) Organization and postembryonic development of glial cells in the adult central brain of Drosophila. The Journal of neuroscience : the official journal of the Society for Neuroscience 28(51):13742-13753. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIyOC81MS8xMzc0MiI7czo0OiJhdG9tIjtzOjM3OiIvYmlvcnhpdi9lYXJseS8yMDE3LzA2LzE4LzE1MTY2Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 16. 16.Jhaveri D, Sen A, & Rodrigues V (2000) Mechanisms underlying olfactory neuronal connectivity in Drosophila-the atonal lineage organizes the periphery while sensory neurons and glia pattern the olfactory lobe. Developmental biology 226(1):73-87. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1006/dbio.2000.9855&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=10993675&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000089705800006&link_type=ISI) 17. 17.Jefferis GS, et al. (2004) Developmental origin of wiring specificity in the olfactory system of Drosophila. Development 131(1):117-130. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czo5OiIxMzEvMS8xMTciO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 18. 18.Vosshall LB, Wong AM, & Axel R (2000) An olfactory sensory map in the fly brain. Cell 102(2):147-159. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(00)00021-0&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=10943836&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000088366800006&link_type=ISI) 19. 19.Gao Q, Yuan B, & Chess A (2000) Convergent projections of Drosophila olfactory neurons to specific glomeruli in the antennal lobe. Nature neuroscience 3(8):780-785. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/77680&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=10903570&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000167177300012&link_type=ISI) 20. 20.Doherty J, Logan MA, Tasdemir OE, & Freeman MR (2009) Ensheathing glia function as phagocytes in the adult Drosophila brain. The Journal of neuroscience : the official journal of the Society for Neuroscience 29(15):4768-4781. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIyOS8xNS80NzY4IjtzOjQ6ImF0b20iO3M6Mzc6Ii9iaW9yeGl2L2Vhcmx5LzIwMTcvMDYvMTgvMTUxNjYyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 21. 21.Kazama H, Yaksi E, & Wilson RI (2011) Cell death triggers olfactory circuit plasticity via glial signaling in Drosophila. The Journal of neuroscience : the off icial journal of the Society for Neuroscience 31(21):7619-7630. 22. 22.Jenett A, et al. (2012) A GAL4-driver line resource for Drosophila neurobiology. Cell reports 2(4):991-1001. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2012.09.011&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=23063364&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 23. 23.Venken KJ, et al. (2011) MiMIC: a highly versatile transposon insertion resource for engineering Drosophila melanogaster genes. Nature methods 8(9):737-743. 24. 24.Kremer MC, Jung C, Batelli S, Rubin GM, & Gaul U (2017) The glia of the adult Drosophila nervous system. Glia 65(4):606-638. 25. 25.Halter DA, et al. (1995) The homeobox gene repo is required for the differentiation and maintenance of glia function in the embryonic nervous system of Drosophila melanogaster. Development 121(2):317-332. [Abstract](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czo5OiIxMjEvMi8zMTciO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 26. 26.Stork T, Sheehan A, Tasdemir-Yilmaz OE, & Freeman MR (2014) Neuron-glia interactions through the Heartless FGF receptor signaling pathway mediate morphogenesis of Drosophila astrocytes. Neuron 83(2):388-403. 27. 27.Lee T & Luo L (1999) Mosaic analysis with a repressible cell marker for studies of gene function in neuronal morphogenesis. Neuron 22(3):451-461. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/S0896-6273(00)80701-1&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=10197526&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000079483900010&link_type=ISI) 28. 28.Beiman M, Shilo BZ, & Volk T (1996) Heartless, a Drosophila FGF receptor homolog, is essential for cell migration and establishment of several mesodermal lineages. Genes & development 10(23):2993-3002. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXNkZXYiO3M6NToicmVzaWQiO3M6MTA6IjEwLzIzLzI5OTMiO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 29. 29.Shishido E, Higashijima S, Emori Y, & Saigo K (1993) Two FGF-receptor homologues of Drosophila: one is expressed in mesodermal primordium in early embryos. Development 117(2):751-761. [Abstract](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czo5OiIxMTcvMi83NTEiO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 30. 30.Gisselbrecht S, Skeath JB, Doe CQ, & Michelson AM (1996) heartless encodes a fibroblast growth factor receptor (DFR1/DFGF-R2) involved in the directional migration of early mesodermal cells in the Drosophila embryo. Genes & development 10(23):3003-3017. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXNkZXYiO3M6NToicmVzaWQiO3M6MTA6IjEwLzIzLzMwMDMiO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 31. 31.Michelson AM, Gisselbrecht S, Zhou Y, Baek KH, & Buff EM (1998) Dual functions of the heartless fibroblast growth factor receptor in development of the Drosophila embryonic mesoderm. Developmental genetics 22(3):212-229. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1002/(SICI)1520-6408(1998)22:3<212::AID-DVG4>3.0.CO;2-9&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=9621429&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 32. 32.Franzdottir SR, et al. (2009) Switch in FGF signalling initiates glial differentiation in the Drosophila eye. Nature 460(7256):758-761. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/nature08167&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=19597479&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000268670300043&link_type=ISI) 33. 33.Shishido E, Ono N, Kojima T, & Saigo K (1997) Requirements of DFR1/Heartless, a mesoderm-specific Drosophila FGF-receptor, for the formation of heart, visceral and somatic muscles, and ensheathing of longitudinal axon tracts in CNS. Development 124(11):2119-2128. [Abstract](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czoxMToiMTI0LzExLzIxMTkiO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 34. 34.Mandal L, Banerjee U, & Hartenstein V (2004) Evidence for a fruit fly hemangioblast and similarities between lymph-gland hematopoiesis in fruit fly and mammal aorta-gonadal-mesonephros mesoderm. Nature genetics 36(9):1019-1023. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/ng1404&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=15286786&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000223658100032&link_type=ISI) 35. 35.Hay BA, Wolff T, & Rubin GM (1994) Expression of baculovirus P35 prevents cell death in Drosophila. Development 120(8):2121-2129. [Abstract](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czoxMDoiMTIwLzgvMjEyMSI7czo0OiJhdG9tIjtzOjM3OiIvYmlvcnhpdi9lYXJseS8yMDE3LzA2LzE4LzE1MTY2Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 36. 36.Sarov M, et al. (2016) A genome-wide resource for the analysis of protein localisation in Drosophila. eLife 5:el2068. 37. 37.Stathopoulos A, Tam B, Ronshaugen M, Frasch M, & Levine M (2004) pyramus and thisbe: FGF genes that pattern the mesoderm of Drosophila embryos. Genes & development 18(6):687-699. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXNkZXYiO3M6NToicmVzaWQiO3M6ODoiMTgvNi82ODciO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 38. 38.Gryzik T & Muller HA (2004) FGF8-like1 and FGF8-like2 encode putative ligands of the FGF receptor Htl and are required for mesoderm migration in the Drosophila gastrula. Current biology : CB 14(8):659-667. 39. 39.Klingseisen A, Clark IB, Gryzik T, & Muller HA (2009) Differential and overlapping functions of two closely related Drosophila FGF8-like growth factors in mesoderm development. Development 136(14):2393-2402. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czoxMToiMTM2LzE0LzIzOTMiO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 40. 40.Diao F, et al. (2015) Plug-and-play genetic access to drosophila cell types using exchangeable exon cassettes. Cell reports 10(8):1410-1421. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2015.01.059&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=25732830&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 41. 41.Newsome TP, Asling B, & Dickson BJ (2000) Analysis of Drosophila photoreceptor axon guidance in eye-specific mosaics. Development 127(4):851-860. [Abstract](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czo5OiIxMjcvNC84NTEiO3M6NDoiYXRvbSI7czozNzoiL2Jpb3J4aXYvZWFybHkvMjAxNy8wNi8xOC8xNTE2NjIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 42. 42.Hong W, et al. (2009) Leucine-rich repeat transmembrane proteins instruct discrete dendrite targeting in an olfactory map. Nature neuroscience 12(12):1542-1550. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/nn.2442&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=19915565&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 43. 43.Chou YH, et al. (2010) Diversity and wiring variability of olfactory local interneurons in the Drosophila antennal lobe. Nature neuroscience 13(4):439-449. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/nn.2489&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=20139975&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000276073500013&link_type=ISI) 44. 44.Stocker RF, Lienhard MC, Borst A, & Fischbach KF (1990) Neuronal architecture of the antennal lobe in Drosophila melanogaster. Cell and tissue research 262(1):9-34. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1007/BF00327741&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=2124174&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=A1990EA33700002&link_type=ISI) 45. 45.Gohl DM, et al. (2011) A versatile in vivo system for directed dissection of gene expression patterns. Nature methods 8(3):231-237. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/nmeth.1561&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=21473015&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 46. 46.Hong W, Mosca TJ, & Luo L (2012) Teneurins instruct synaptic partner matching in an olfactory map. Nature 484(7393):201-207. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/nature10926&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=22425994&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000303149900027&link_type=ISI) 47. 47.Komiyama T & Luo L (2006) Development of wiring specificity in the olfactory system. Current opinion in neurobiology 16(1):67-73. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/j.conb.2005.12.002&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=16377177&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000236136200010&link_type=ISI) 48. 48.Woolsey TA & Van der Loos H (1970) The structural organization of layer IV in the somatosensory region (SI) of mouse cerebral cortex. The description of a cortical field composed of discrete cytoarchitectonic units. Brain research 17(2):205-242. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/0006-8993(70)90079-X&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=4904874&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=A1970F322900004&link_type=ISI) 49. 49.Ornitz DM & Itoh N (2015) The Fibroblast Growth Factor signaling pathway. Wiley interdisciplinary reviews. Developmental biology 4(3):215-266. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1002/wdev.176&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=25772309&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 50. 50.Muha V & Muller HA (2013) Functions and Mechanisms of Fibroblast Growth Factor (FGF) Signalling in Drosophila melanogaster. International journal of molecular sciences 14(3):5920-5937. 51. 51.Guillemot F & Zimmer C (2011) From cradle to grave: the multiple roles of fibroblast growth factors in neural development. Neuron 71(4):574-588. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2011.08.002&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=21867876&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000294521600005&link_type=ISI) 52. 52.Smith KM, et al. (2006) Midline radial glia translocation and corpus callosum formation require FGF signaling. Nature neuroscience 9(6):787-797. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/nn1705&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=16715082&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000237895200019&link_type=ISI) 53. 53.Perraud F, Labourdette G, Miehe M, Loret C, & Sensenbrenner M (1988) Comparison of the morphological effects of acidic and basic fibroblast growth factors on rat astroblasts in culture. Journal of neuroscience research 20(1):1-11. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1002/jnr.490200102&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=3047413&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=A1988P116700001&link_type=ISI) 54. 54.Reilly JF, Maher PA, & Kumari VG (1998) Regulation of astrocyte GFAP expression by TGF-beta1 and FGF-2. Glia 22(2):202-210. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1002/(SICI)1098-1136(199802)22:2<202::AID-GLIA11>3.0.CO;2-1&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=9537840&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000071771600011&link_type=ISI) 55. 55.Furusho M, Dupree JL, Nave KA, & Bansal R (2012) Fibroblast growth factor receptor signaling in oligodendrocytes regulates myelin sheath thickness. The Journal of neuroscience : the official journal of the Society for Neuroscience 32(19):6631-6641. [Abstract/FREE Full Text](http://biorxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIzMi8xOS82NjMxIjtzOjQ6ImF0b20iO3M6Mzc6Ii9iaW9yeGl2L2Vhcmx5LzIwMTcvMDYvMTgvMTUxNjYyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 56. 56.Kang K, Lee SW, Han JE, Choi JW, & Song MR (2014) The complex morphology of reactive astrocytes controlled by fibroblast growth factor signaling. Glia 62(8):1328-1344. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1002/glia.22684&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=24796693&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 57. 57.Gibson NJ, Tolbert LP, & Oland LA (2012) Activation of glial FGFRs is essential in glial migration, proliferation, and survival and in glia-neuron signaling during olfactory system development. PloS one 7(4):e33828. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0033828&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=22493675&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) 58. 58.Shimokawa K, et al. (2011) Cell surface heparan sulfate chains regulate local reception of FGF signaling in the mouse embryo. Developmental cell 21(2):257-272. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/j.devcel.2011.06.027&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=21839920&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000294387300008&link_type=ISI) 59. 59.Baeg GH & Perrimon N (2000) Functional binding of secreted molecules to heparan sulfate proteoglycans in Drosophila. Current opinion in cell biology 12(5):575-580. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/S0955-0674(00)00134-4&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=10978892&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000089203800007&link_type=ISI) 60. 60.Sweeney LB, et al. (2007) Temporal target restriction of olfactory receptor neurons by Semaphorin-1a/PlexinA-mediated axon-axon interactions. Neuron 53(2):185-200. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2006.12.022&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=17224402&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000245126600006&link_type=ISI) 61. 61.Komiyama T, Carlson JR, & Luo L (2004) Olfactory receptor neuron axon targeting: intrinsic transcriptional control and hierarchical interactions. Nature neuroscience 7(8):819-825. [CrossRef](http://biorxiv.org/lookup/external-ref?access_num=10.1038/nn1284&link_type=DOI) [PubMed](http://biorxiv.org/lookup/external-ref?access_num=15247920&link_type=MED&atom=%2Fbiorxiv%2Fearly%2F2017%2F06%2F18%2F151662.atom) [Web of Science](http://biorxiv.org/lookup/external-ref?access_num=000222930800011&link_type=ISI)