http://orcid.org/0000-0002-7906-7173
Figures
Abstract
Synaptic structure and activity are sensitive to environmental alterations. Modulation of synaptic morphology and function is often induced by signals from glia. However, the process by which glia mediate synaptic responses to environmental perturbations such as hypoxia remains unknown. Here, we report that, in the mutant for Trachealess (Trh), the Drosophila homolog for NPAS1 and NPAS3, smaller synaptic boutons form clusters named satellite boutons appear at larval neuromuscular junctions (NMJs), which is induced by the reduction of internal oxygen levels due to defective tracheal branches. Thus, the satellite bouton phenotype in the trh mutant is suppressed by hyperoxia, and recapitulated in wild-type larvae raised under hypoxia. We further show that hypoxia-inducible factor (HIF)-1α/Similar (Sima) is critical in mediating hypoxia-induced satellite bouton formation. Sima upregulates the level of the Wnt/Wingless (Wg) signal in glia, leading to reorganized microtubule structures within presynaptic sites. Finally, hypoxia-induced satellite boutons maintain normal synaptic transmission at the NMJs, which is crucial for coordinated larval locomotion.
Author summary
Oxygen is essential for animals to maintain their life such as growth, metabolism, responsiveness, and movement. It is therefore important to understand how animal cells trigger hypoxia response and adapt to hypoxia thereafter. Both mammalian vascular and insect tracheal branches are induced to enhance the oxygen delivery. However, the study of hypoxia response in the nervous system remains limited. In this study, we assess the morphology of Drosophila neuromuscular junctions (NMJs), a model system to study development and function of synapses, in two hypoxia conditions, one with raising wild-type larvae in hypoxia, and the other in the trachealess (trh) mutant in which the trachea is defective, causing insufficient oxygen supply. Interestingly, the glial processes, normally wrapping around the axons, invade into synaptic boutons of NMJs under hypoxia. Also, the hypoxia-induced factor Sima activates the Wg signal in glia and the secreted Wg signal reorganizes the synaptic boutons as a response to hypoxia. This synaptic bouton reorganization might maintain normal synaptic transmission and locomotion of larvae.
Citation: Chen P-Y, Tsai Y-W, Cheng Y-J, Giangrande A, Chien C-T (2019) Glial response to hypoxia in mutants of NPAS1/3 homolog Trachealess through Wg signaling to modulate synaptic bouton organization. PLoS Genet 15(8): e1007980. https://doi.org/10.1371/journal.pgen.1007980
Editor: Hongyan Wang, Duke-NUS Medical School, SINGAPORE
Received: January 18, 2019; Accepted: July 15, 2019; Published: August 5, 2019
Copyright: © 2019 Chen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This study was supported by an internal grant from the Institute of Molecular Biology, Academia Sinica to C-TC. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Animals need oxygen and food, not only to sustain life, but also for motility. In vertebrates, oxygen and nutrients are delivered through the vascular systems to organs and tissues throughout the body. To maintain proper nutrient and oxygen supply, and thus physiological functions, the vascular system is also highly coordinated with the nervous system during development. Indeed, the vascular and nervous systems resemble each other in terms of their anatomical structures and developmental processes [1]. In the brain, nerves and vessels, form close associations and are in physical contact through the third player astrocytes to form neurovascular units (NVU) [2]. Such organization is essential for controlling oxygen and glucose delivery through the blood vessels by neuronal activity, and this regulatory process is mediated through the coupled astrocytes [3]. However, some invertebrates lack the complex vascular systems [4]. In nematodes, oxygen is supplied simply by ambient diffusion to inner cells [5]. Insects such as Drosophila have evolved a prototype of the tracheal system to deliver oxygen and a primitive vascular system, the dorsal vessel, to facilitate nutrient delivery [6]. However, the physical association of nerves, trachea, and glial processes has also been demonstrated at the NMJs of adult Drosophila flight muscles [7].
Animals respond to changing oxygen levels by altering their oxygen delivery system. Insufficient oxygen levels (hypoxia) activate a broad range of genes to re-establish body homeostasis. One crucial regulator of these hypoxia-responsive genes is the sequence-specific DNA-binding transcription factor hypoxia inducible factor 1 (HIF-1) [8]. HIF-1 consists of α and β subunits that form heterodimers [9]. Whereas HIF-1β is expressed constitutively, HIF-1α protein levels are modulated by oxygen levels [10]. Under normal oxygen conditions (normoxia), oxygen-dependent prolyl hydroxylases (PHDs) catalyze hydroxylation of a conserved prolyl residue in the central oxygen-dependent degradation (ODD) domain of HIF-1α [11]. Hydroxylation of HIF-1α promotes interaction with Von Hippel Lindau (VHL), which is the substrate recognition subunit of the cullin2-based E3 ubiquitin ligase, leading to HIF-1α ubiquitination and proteasomal degradation. Under hypoxia, prolyl hydroxylation does not occur, HIF-1α proteins are stabilized and are translocated from the cytoplasm to the nucleus where they form heterodimers with HIF-1β to activate transcription of target genes. One major class of target genes encoding the Fibroblast Growth Factor (FGF) is involved in inducing angiogenesis in mammals. In Drosophila, the FGF member encoded by Branchless (Bnl) induces tracheal branching [12]. When oxygen levels are reduced, oxygen-starved cells express Bnl as a chemo-attractant to guide the growth tracheal terminal branches toward them [13].
In addition to adaptations of the respiratory system, the nervous system also responds to hypoxia. Oxygen levels modulate the survival, proliferation, and differentiation of radial glial cells (RGCs) in the human cerebral cortex. Interestingly, physiological hypoxia (3% O2) induces neurogenesis and differentiation of RGCs into glutamatergic neurons [14]. Hypoxia induces neurite outgrowth in PC12 cells through activation of A2A receptor [15]. Brief exposure to anoxia and hypoglycemia caused axonal remodeling in hippocampal neurons, including presynaptic protrusion of filopodia and formation of multi-innervated spines [16]. Under hypoxia or upon depletion of PHD2, upregulation of the actin cross-linker Filamin A (FLNA) induces generation of more immature spines [17]. Astrocytes have been shown to play a crucial role in ischemic tolerance via the activation of P2X7 receptors, which trigger upregulation of HIF-1α [18].
Neuronal PAS (NPAS) proteins containing a DNA-binding Per-Arnt-Sim domain function in vascular and nervous system development. In mice, NPAS1 is responsible for cortical interneuron generation [19], whereas NPAS3 is required for adult neurogenesis [20]. NPAS1 and NPAS3 also play key roles in lung development [21, 22]. The homolog of NPAS1/3 in Drosophila, Trachealess (Trh), has been well studied for its involvement in formation of the respiratory tracheal system. Trh is a master regulator of tracheal cell fates, activating gene expression to induce tracheal development [23, 24]. However, the role of Trh in the development of other tissues, particularly the nervous system, is unknown. In this study, we found altered synaptic bouton morphology at the NMJs of trh1/trh2 mutants. By performing trh-RNAi knockdown and UAS-trh transgene rescue experiments, we show that trh is required in tracheal cells for normal bouton formation. Defective tracheal branching in the trh1/trh2 mutant mimics the effect of hypoxic conditions during larval development, and supplying higher than normal oxygen levels restored normal bouton morphology. We further show that glial cells respond to hypoxia by elevating Wnt/Wg expression to mediate synaptic bouton reorganization through HIF1-α/Sima in Drosophila. Finally, we reveal that this morphological change may be linked to normal synaptic transmission and locomotion of larvae.
Results
trh modulates synaptic bouton formation non-cell autonomously
To understand the possible role of Trh in synapse formation, we examined NMJ morphology in the trh mutant. Since both trh1 and trh2 loss-of-function alleles are homozygous lethal [23, 25, 26], we examined the trans-heterozygous trh1/trh2 mutant that survived to adult stages and compared it to wild-type (w1118) and heterozygous trh1/+ controls. Synaptic boutons of w1118 and trh1/+ NMJs were evenly spaced along the axonal terminals, displaying the typical “beads-on-a-string” pattern (Fig 1A, upper and middle panels, enlarged images at right). Strikingly, the trh1/trh2 mutant larvae exhibited aberrant NMJ morphology. Small synaptic boutons formed clusters or surrounded a normal large bouton (Fig 1A, bottom panel); a phenotype described as “satellite boutons” [27]. This satellite bouton phenotype in the trh1/trh2 mutant was detected at a high frequency; 20.7 ± 4.1% (n = 10) of total boutons were satellite ones, more than three-fold increases compared to 4.3 ± 1.7% (n = 10) for trh1/+ and 6.3 ± 1.7% (n = 10) for w1118 (Fig 1B). Although the percentage of satellite boutons in the trh1/trh2 mutant was greatly increased, the total bouton number was slightly higher than that observed in controls (74.1 ± 4.7, n = 10 in w1118; 89.6 ± 5.2, n = 10 in trh1/+; and 95.5 ± 9.1, n = 11 in trh1/trh2, bottom panel in Fig 1B). Also, the muscle areas were not significantly different from each other (S1A Fig). Given the small size and clustering of satellite boutons in the trh1/trh2 mutant, we examined whether these boutons express synaptic proteins normally. We found that the synaptic vesicle protein Synapsin (Syn in Fig 1A) was normally distributed relative to control, but the active zone protein Bruchpilot (Brp) was expressed at higher levels in the trh1/trh2 mutant (S1B Fig). The postsynaptic glutamate receptor, as revealed by GluRIII (S1B Fig) and GluRIIA (S1C Fig) signals, as well as dPAK (S1C Fig) were also localized in satellite boutons, which were surrounded by the subsynaptic reticulum protein Dlg (S1D Fig). Thus, although the Brp signal intensity in the trh1/trh2 mutant was stronger than in controls, the composition of synaptic proteins in satellite boutons was largely similar to that of normal-sized boutons.
(A, C, E) Images show NMJs of muscle 6/7 immunostaining for Syn (green) and HRP (red) in w1118, trh1/+, and trh1/trh2 (A), btl-GAL4 and btl-GAL4/trh-RNAi (C), and UAS-trh trh1/trh2, btl-GAL4 trh1/trh2, and btl-GAL4/UAS-trh trh1/trh2 (E). Dashed squares in (A) are enlarged in right panels with white arrows indicating clusters of satellite boutons. Scale bars are 10 μm. (B, D, F) Bar graphs show mean ± standard error of mean for percentages of satellite boutons to total boutons (satellite + normal boutons)(B, upper graph, D, F) or total bouton number (B, bottom graph) for each genotype. Dots in each bar represent the distribution and number of data. Statistical significance by Mann-Whitney test is shown (n.s., no significance; *, p < 0.5; **, p < 0.01; ***, p < 0.001).
As trh is expressed in both tracheal and nervous systems in embryonic stages [28], altered bouton morphology in the trh1/trh2 mutant could be due to a lack of trh in neurons, tracheal cells or other cells/tissues. Therefore, we performed trh-RNAi knockdown by using tissue-specific GAL4 drivers for trachea (btl-GAL4), neurons (elav-GAL4), glia (repo-GAL4), and muscles (MHC-GAL4). We observed a dramatic increase in satellite boutons upon tracheal trh-RNAi knockdown using btl-GAL4 (10.2 ± 1.7%, n = 10) compared to those in the trh-RNAi (2.7 ± 1.3%, n = 10) or btl-GAL4 (3.3 ± 0.9%, n = 12) control (Fig 1C and 1D. In contrast, trh-RNAi knockdown by other tissue-specific GAL4 drivers failed to significantly increase satellite boutons, displaying similar amounts of satellite boutons to respective GAL4 drivers (2.0 ± 0.7%, n = 10 in elav-GAL4 v.s. 1.4 ± 0.6%, n = 10 in elav-GAL4/trh-RNAi; 2.8 ± 1.0%, n = 10 in repo-GAL4 v.s. 3.4 ± 0.9%, n = 10 in repo-GAL4/trh-RNAi; and 1.3 ± 0.8%, n = 10 in MHC-GAL4 v.s. 2.1 ± 1.0%, n = 10 in MHC-GAL4/trh-RNAi, Fig 1D).
To further confirm the necessity of tracheal trh for normal bouton formation, we performed a rescue experiment for the trh1/trh2 phenotype. When UAS-trh expression was driven by tracheal btl-GAL4 in the trh1/trh2 mutant, the satellite bouton phenotype was suppressed (1.8 ± 1.4%, n = 13, Fig 1E and 1F). Controls bearing only btl-GAL4 (14.3 ± 3.3%, n = 10) or UAS-trh (21.1 ± 3.3%, n = 10) still contained large amounts of satellite boutons. These results indicate that trh is required in the trachea for normal bouton formation.
Hypoxia induces satellite bouton formation
Apart from specifying the tracheal cell fate, Trh is also involved in the branching of tubular structures during post-embryonic stages [24]. Therefore, we examined the tracheal phenotypes in the trh1/trh2 larvae (S2A Fig) and observed an increase in the number of terminal branches in the dorsal branch of the third segment (S2B Fig, 5.7 ± 0.15, n = 10 for trh1/+, and 7.5 ± 0.28, n = 11 for trh1/trh2). Furthermore, we identified morphological defects such as tracheal breaks and tangles, suggesting structural defects in the trh1/trh2 larvae (arrows in S2A Fig). Tracheal branching activity is enhanced under hypoxia [12]. Thus, the increased terminal branches in trh1/trh2 could be a compensatory mechanism for defective trachea formation.
To understand whether trh1/trh2 mutant cells are under hypoxia, we used the hypoxia biosensor GFP-ODD, in which the GFP is fused to the ODD domain of Sima, under the control of the ubiquitin-69E (ubi) promoter [29]. We first confirmed that GFP-ODD signal was low under normoxia (21% O2) and enhanced under hypoxia (5% O2) in wild-type late-stage embryos when tracheal tubules are already formed and functioning [29]. Indeed, enhanced GFP signal was ubiquitous under hypoxia in wild-type embryos with some pronounced focal GFP signals (Fig 2A, upper row, and S2C Fig). The signal of mRFP-nls, also under the control of the ubi promoter as an internal control, remained constant under hypoxia (Fig 2A, bottom row, and S2D Fig). Quantification of the GFP/RFP ratio revealed a significant difference between normoxia and hypoxia (0.18 ± 0.05, n = 6 at 21% O2 and 0.84 ± 0.15, n = 6 at 5% O2, Fig 2B). We then examined whether oxygen supply is deficient in the trh1/trh2 mutant by measuring the GFP/RFP ratios. We detected a higher GFP/RFP ratio (0.82 ± 0.15, n = 6) in the mutant compared to that in heterozygous trh1/+ (0.09 ± 0.02, n = 6) in the 21% O2 condition, supporting that the trh1/trh2 mutant senses reduced oxygen levels internally (Fig 2A and 2B).
(A) Images of late-stage embryos of w1118, trh1/+, and trh1/trh2 carrying the GFP-ODD reporter (top row) and mRFP-nls control (bottom row) under normoxia (21% O2) or hypoxia (5% O2). Scale bar represents 40 μm. (C) Images of NMJs from muscle 6/7 immunostained for Syn (green) and HRP (red) from w1118 in 21% or 5% O2 (left panels) and trh1/trh2 in 21% or 50% O2 (right panels). White arrows indicate clusters of satellite boutons. Scale bar represents 10 μm. (B, D) Bar graphs show averages of the GFP/RFP immunofluorescence intensity ratios (B), or percentages of satellite boutons to total boutons (D). Statistical significance assayed by Mann-Whitney test is shown (**, p < 0.01; ***, p < 0.001).
Thus, formation of satellite boutons in the trh1/trh2 mutant could be caused by hypoxia. To test this hypothesis, we reared wild-type larvae under hypoxia (5%) and assessed synaptic bouton morphology at the third instar stage (Fig 2C). Consistently, small clustered boutons increased, mirroring the satellite bouton phenotype (Fig 2D, 4.6 ± 1.5%, n = 10 in 21% O2 v.s. 19.4 ± 1.2%, n = 11 in 5% O2). Furthermore, when the trh1/trh2 mutant was raised in a high oxygen level (50% O2), normal bouton morphology was restored (Fig 2C); the satellite bouton phenotype was almost completely suppressed (Fig 2D, 15.2 ± 2.5%, n = 10 in 21% O2 v.s. 0.74 ± 0.49%, n = 10 in 50% O2). These results suggest that hypoxia due to the defective tracheal system in the trh1/trh2 mutant induces the satellite bouton phenotype, and that this phenotype can be suppressed by extra oxygen supply.
Glial HIF-1α/Sima mediates satellite bouton formation
HIF-1α/Sima mediates the response to low oxygen supply [30, 31]. The protein levels of Sima are increased in wild-type Drosophila embryos subjected to hypoxia [32], leading to transcriptional activation of downstream target genes and the induction of tracheal branching [12]. We overexpressed Sima in tracheal cells, neurons, glia, or muscle cells by tissue-specific drivers to investigate which types of cells may play a role in modulating synaptic bouton formation. Overexpressing Sima in trachea caused embryonic lethality, preventing us from observing NMJ phenotypes. Larvae in which Sima was overexpressed in muscles, neurons, or glia could survive to the third instar stage, allowing us to examine the bouton phenotypes. We found that Sima overexpression in glia gave the highest number of satellite boutons (Fig 3A, and 11.9 ± 2.3%, n = 10 for repo-GAL4/UAS-sima in Fig 3B, v.s. 2.8 ± 1.0%, n = 10 for repo-GAL4 in Fig 1D). However, Sima overexpression in neurons by elav-GAL4 failed to induce satellite boutons (2.6 ± 0.8%, n = 8 in elav-GAL4/UAS-sima in Fig 3B v.s. 2.0 ± 0.7%, n = 10 for elav-GAL4 in Fig 1D). Also, muscle expression of Sima also maintained a basal level of satellite boutons (1.7 ± 1.0%, n = 10 for MHC-GAL4/UAS-sima in Fig 3B v.s. 1.3 ± 0.8%, n = 10 for MHC-GAL4 in Fig 1D). This result shows that the hypoxia-responding factor Sima is capable of inducing satellite bouton formation when overexpressed in glia.
(A, C, E) Images showing NMJs from muscle 6/7 immunostained for Syn (green) and HRP (red) in repo-GAL4, and repo-GAL4/UAS-sima (A), sima-RNAi trh1/trh2, repo-GAL4 trh1/trh2, and repo-GAL4/sima-RNAi trh1/trh2 (C), sima-RNAi, repo-GAL4, and repo-GAL4/sima-RNAi under hypoxia at 5% O2 (E). White arrows indicate clusters of satellite boutons. Scale bars represent 10 μm. (B, D, F) Bar graphs show percentages (mean ± SEM) of satellite boutons. Statistical significance was assayed by one-way ANOVA with Tukey's Multiple Comparison post-test (B) or Mann-Whitney test (D, F) with n.s., no significance; *, p <0.05; **, p < 0.01; ***, p < 0.001.
The Drosophila PHD protein Fatiga (Fga) promotes prolyl hydroxylation and degradation of Sima in normaxia [33]. Consistently, glial knockdown by expressing fga-RNAi led to a significant increase of satellite boutons (S3A and S3B Fig, 8.3 ± 1.3%, n = 10 for repo-GAL4/fga-RNAi) from the control (4.3 ± 0.7%, n = 10 for fga-RNAi). Instead, tracheal knockdown by btl-GAL4 (3.8 ± 0.7%, n = 9 for btl-GAL4/fga-RNAi) showed no apparent difference to the control. This result is consistent with that Sima upregulation in glia induces satellite bouton formation.
If glial Sima is the factor responding to hypoxia in the trh1/trh2 mutant, reducing Sima level in glia would suppress satellite bouton formation. Accordingly, we expressed the sima-RNAi transgene, which could reduce sima expression in both transcript and protein levels (S3C and S3D Fig), using repo-GAL4 in the trh1/trh2 mutant. As our prediction, the satellite bouton phenotype was suppressed upon glial sima knockdown (3.6 ± 1.0%, n = 10 in repo-GAL4/sima-RNAi trh1/trh2), as compared to controls of the trh1/trh2 mutant carrying either the UAS-sima-RNAi transgene (9.2 ± 1.5%, n = 10 in sima-RNAi trh1/trh2) or the repo-GAL4 driver (12.2 ± 2.2%, n = 10 in repo-GAL4 trh1/trh2) that displayed high percentages of satellite boutons (Fig 3C and 3D). We also tested whether low oxygen level-induced satellite bouton formation is mediated through Sima in glia. Satellite bouton phenotypes were detected in controls carrying either UAS-sima-RNAi (12.5 ± 2.7%, n = 11) or repo-GAL4 (16.4 ± 2.2%, n = 10) when raised in 5% O2. However, almost no satellite boutons were detected in larvae carrying both repo-GAL4 and UAS-sima-RNAi (0.45 ± 0.45%, n = 11 in repo-GAL4/sima-RNAi) when raised in the same condition (Fig 3E and 3F).
We also examined whether the Sima protein levels are changed in hypoxia or in the trh mutant. We found ubiquitous increases in the Sima levels in the wild-type control under the 5% O2 condition or in the trh1/trh2 mutant (S3E and S3F Fig). The increases could be identified in glial processes along the peripheral nerves and in different subtypes of glia. Thus, glial Sima could play the role to mediate hypoxia in the trh1/trh2 mutant and in the low O2 condition to modulate synaptic bouton formation.
Wg signals mediate glial Sima activity to modulate bouton morphology
Next, we explored possible signals transduced from glia to neurons in response to hypoxia. The glia-secreted Wingless (Wg) signaling molecule regulates synaptic growth at Drosophila NMJs [34, 35]. Therefore, we examined whether Wg can be induced under hypoxia in the trh1/trh2 mutant. Wg signals were enriched around the synaptic boutons of wild-type NMJs (Fig 4A). Whereas the pattern of Wg signals at trh1/+ NMJs was similar to that of w1118, we detected much higher levels of Wg signals at the trh1/trh2 NMJ (Fig 4A). Quantification of Wg immunofluorescence intensities normalized to co-stained HRP in trh1/trh2 (Wg/HRP: 0.49 ± 0.07, n = 9, Fig 4B) revealed ~3-fold increases relative to w1118 (0.18 ± 0.03, n = 8) and trh1/+ (0.15 ± 0.02, n = 9). We then examined whether glial Sima is required for the enhanced Wg expression in the trh1/trh2 mutant (Fig 4C). The Wg level relative to the HRP level in the trh1/trh2 mutant carrying repo-GAL4 (0.68 ± 0.16, n = 8 for repo-GAL4 trh1/trh2) was also about 3 folds to the repo-GAL4 control (0.23 ± 0.04, n = 11, Fig 4D). When we reduced sima levels in the trh1/trh2 mutant by repo-GAL4-driven UAS-sima-RNAi, Wg signals were suppressed to a level equivalent to that in the repo-GAL4 control (0.25 ± 0.05, n = 9 for repo-GAL4/sima-RNAi trh1/trh2). Interestingly, sima-RNAi knockdown in glia of the repo-GAL4 control had no effect on the Wg level (0.25 ± 0.02, n = 10 for repo-GAL4/sima-RNAi), suggesting that Sima is induced in the trh1/trh2 mutant to upregulate Wg expression but has no role in basal Wg expression in the wild-type. Taken together, we suggest that glial Sima is required for Wg upregulation at the NMJs of the trh1/trh2 mutant.
(A, C, E) Images showing NMJs of muscle 6/7 immunostained for Wg (green) and HRP (red) in w1118, trh1/+, and trh1/trh2 (A); repo-GAL4, repo-GAL4/sima-RNAi, repo-GAL4 trh1/trh2, and repo-GAL4/sima-RNAi trh1/trh2 (C); and for Syn (green) and HRP (red) in wg-RNAi trh1/trh2 and repo-GAL4/wg-RNAi trh1/trh2 (E). White arrows indicate clusters of satellite boutons. Scale bars represent 5 μm. (B, D, F) Bar graphs show averages (mean ± SEM) of normalized Wg to HRP intensities (B, D), and percentages (mean ± SEM) of satellite boutons (F). Statistical significance by Mann-Whitney test is shown (n.s., no significance; *, p < 0.05; ***, p < 0.001).
If glia-secreted Wg is responsible for satellite bouton induction in the trh1/trh2 mutant, then glia-specific wg knockdown should phenocopy sima knockdown in glia to suppress satellite bouton formation. In control, satellite boutons were still prominent in the trh1/trh2 mutant bearing only UAS-wg-RNAi (Fig 4E and 4F, 17.0 ± 2.8%, n = 8 for wg-RNAi trh1/trh2). The satellite bouton phenotype was suppressed in the trh1/trh2 mutant bearing both repo-GAL4 and UAS-wg-RNAi (4.4 ± 1.4%, n = 10 for repo-GAL4/wg-RNAi trh1/trh2). Also, overexpression of Sima by repo-GAL4 induced satellite bouton formation (Fig 3B, 11.9 ± 2.3%, n = 10 for repo-GAL4/UAS-sima), which was suppressed by co-expressing the wg-RNAi transgene (Fig 3B, 2.3 ± 1.3%, n = 7 for repo-GAL4/UAS-sima UAS-wg-RNAi). Taken together, these results strongly suggest that glia-secreted Wg mediates Sima activity in promoting satellite bouton formation in the trh mutant.
As Wg signals are secreted from both glia and presynaptic neurons [34, 36], we found that reduction of Wg signals from neurons also suppressed satellite bouton formation in the trh1/trh2 mutant (Fig 5A and 5B, 10.1 ± 1.5%, n = 10 for elav-GAL4 trh1/trh2; 6.0 ± 1.0%, n = 8 for elav-GAL4/wg-RNAi trh1/trh2). Also, neuronal and glial but not muscle overexpression of Wg induced satellite bouton formation (S4A and S4B Fig). Glial and neuronal knockdown of Wg reduced total Wg levels at the trh1/trh2 NMJ, although the reduction by neuronal knockdown was not significant (S4C and S4D Fig). Thus, neuronal Wg also contributes to satellite bouton formation. However, only glial overexpression of Sima induced higher levels of Wg, but not neuronal or muscle overexpression of Sima (S4E and S4F Fig). Thus, while neuronal expression of Wg contributes to the overall level at NMJs in the trh mutant, Sima-induced Wg expression is likely glial-specific.
(A) Images showing NMJ 6/7 immunostained for Syn (green) and HRP (red) in elav-GAL4 trh1/trh2, elav-GAL4/wg-RNAi trh1/trh2, elav-GAL4/fz2-RNAi trh1/trh2, fz2-RNAi trh1/trh2, MHC-GAL4 trh1/trh2, and MHC-GAL4/fz2-RNAi trh1/trh2. Scale bars represent 10 μm. (B) Bar graph shows averages (mean ± SEM) of percentages of satellite boutons. (C) Images showing NMJs of muscle 6/7 immunostained for HRP (red) and Futsch (green) in trh1/+ and trh1/ trh2. White arrows indicate microtubule loops, and images within dotted squares are amplified in right panels. Scale bars represent 10 μm. (D) Bar graphs show average numbers (mean ± SEM) of Futsch loops. (E) Images of live tissues showing NMJs of muscle 6/7 carrying repo-cyto-GFP (green) labeled for HRP (red) in trh1/+ and trh1/trh2. Synaptic overlap by glia is shown as yellow color. Sections of overlapping boutons with dashed lines are shown in bottom panels. Scale bar represents 10 μm. (F) Bar graph shows percentages (mean ± SEM) of synaptic areas overlapping with glial processes. Statistical significance by Mann-Whitney test is shown (n.s., no significance; *, p < 0.05; **, p < 0.01; ***, p < 0.001).
Presynaptic Wg signaling modulates bouton morphology in the trh mutant
At larval NMJs, Wg signaling in pre-synaptic and post-synaptic sites through distinct signaling pathways [35–37]. However, the receptor Frizzled2 (Fz2) is involved in transducing the signaling pathway activities in both sites. We tested whether pre- or post-synaptic Wg signaling is involved in satellite bouton formation by fz2 knockdown in the trh mutant. Neuronal expression of fz2-RNAi by elav-GAL4 suppressed satellite bouton formation in the trh1/trh2 mutant (Fig 5A and 5B, 11.8 ± 2.7%, n = 9 for fz2-RNAi trh1/trh2; 10.1 ± 1.5%, n = 10 for elav-GAL4 trh1/trh2; and 2.5 ± 1.1%, n = 9 for elav-GAL4/fz2-RNAi trh1/trh2). In contrast, muscle expression of fz2-RNAi by MHC-GAL4 had no significant difference to the controls (6.4 ± 1.6%, n = 10 for MHC-GAL4 trh1/trh2 and 7.3 ± 1.6%, n = 8 for MHC-GAL4/fz2-RNAi trh1/trh2). Thus, the elevated Wg level in the trh mutant is mainly transduced in presynaptic sites to modulate satellite bouton formation.
Inactivation of Wg signaling leads to a reduction of the more stabilized microtubule loops within synaptic boutons, which could be visualized by immunostaining for the microtubule-binding protein Futsch [35, 36]. We examined whether elevated Wg signaling alters microtubule loops in presynaptic boutons of the trh1/trh2 mutant by Futsch immunostaining and found significantly more Futsch-positive loops within the boutons of the trh mutant (Fig 5C and 5D, 2.92 ± 0.29, n = 12 for trh1/+; and 9.56 ± 1.11, n = 9 for trh1/trh2), supporting the elevation of presynaptic Wg signaling.
Glial processes invade synaptic boutons to match the growth of NMJs [38], which intrigued us to assess whether glia in the trh1/trh2 mutant exhibits morphological change. In a live fillet preparation for imaging NMJs, we found that glial processes labeled by GFP invaded the area of synaptic boutons in the trh1/trh2 mutant, whereas glial processes were relatively restrained from the bouton areas in the control (Fig 5E). Quantification of the glial process overlaying the synaptic bouton area revealed significantly greater area of overlap in the trh1/trh2 mutant relative to control (Fig 5F, 6.2 ± 1.0%, n = 10 for trh1/+; and 16.7 ± 3.7%, n = 10 for trh1/trh2). This increased extent of glial processes in the synaptic area may facilitate signal transduction from glia to synaptic boutons for structural reorganization. Taken together, these results suggest that Wg plays a prominent role in the trh1/trh2 mutant to transduce the hypoxia signal from glia to modify presynaptic bouton structure.
Impaired crawling behavior in the trh mutant
Given the evident morphological changes at trh1/trh2 NMJs, we wondered if locomotion is affected in mutant larvae. We observed larvae crawling under free-movement conditions and found that wild-type control and trh1/+ heterozygous larvae presented smooth crawling paths (Fig 6A), with an average speed of 0.64 ± 0.09 mm/s (n = 13) in w1118 and 0.45 ± 0.08 mm/s (n = 12) in trh1/+ (Fig 6B). However, the trh1/trh2 larvae had much shorter paths and a slower speed of 0.14 ± 0.02 mm/s (n = 14). The head turning angle of the trh1/trh2 mutant was comparable to both controls, not contributing to the slow movement (Fig 6C, 7.8 ± 0.6 degree/s, n = 13 for w1118; 9.7 ± 1.2 degree/s, n = 12 for trh1/+; and 11.3 ± 1.9 degree/s, n = 14 for trh1/trh2). Larval crawling is a rhythmic behavior involving a series of periodic strides (S1 Movie) [39]. We noticed uncoordinated crawling in the trh1/trh2 larvae, with their posterior body segments failing to follow the rhythmic movement (S2 Movie). We recorded larval forward crawling and constructed kymographs to represent the stride cycle. In wild-type larvae, normal and consistent periodic strides were apparent with regular head and tail displacements (Fig 6D, left panel). Similar to the wild-type, head movements of trh1/trh2 larvae were smooth and periodic, albeit slower. However, tail movements of trh1/trh2 larvae were abrupt (Fig 6D, right panel). While wild-type larvae crawled completely normal (n = 10), 69.2 ± 9.6% (n = 10) of the strides of trh1/trh2 larvae were uncoordinated (Fig 6E), which might contribute to the slower crawling of the trh1/trh2 mutant.
(A) Representative 5-minute crawling traces for w1118, trh1/+, and trh1/trh2. (B, C) Bar graphs show (B) average speeds (mean ± SEM in mm/s), and (C) average rotational speeds (mean ± SEM in degree/s). (D) Kymographs show four forward strides from tail to head (shown by arrow) for w1118 and trh1/trh2. The green lines show smooth transitions between successive strides, whereas orange lines show abrupt transitions. Scale bars for 5s and 1mm are shown. (E) Bar graph shows uncoordinated stride percentages (mean ± SEM). Statistical significance by Mann-Whitney test is shown (n.s., no significance; ***, p < 0.001).
This uncoordinated stride cycle prompted us to examine the bouton morphology in A2-A6 segments of the trh1/trh2 larvae (Fig 7A and 7B). Strikingly, large numbers of satellite boutons were detected in segments A2 (18.0 ± 4.3%, n = 10) and A3 (20.7 ± 4.1%, n = 10) and an intermediate level of satellite boutons was detected in segment A4 (9.2 ± 2.6%, n = 9). The numbers of satellite boutons were low in more posterior segments A5 (4.9 ± 2.3%, n = 10) and A6 (3.2 ± 1.7%, n = 9). In wild-type larvae, all segments had relatively normal bouton morphology except the A3 segment (0.6 ± 0.4%, n = 10 in A2, 6.3 ± 1.7%, n = 10 in A3, 0.7 ± 0.4%, n = 10 in A4, 0.5 ± 0.5%, n = 10 in A5, and 0.9 ± 0.6%, n = 9 in A6). While high numbers of satellite boutons appeared in the trh mutant, the total numbers of boutons in individual segments were comparable between wild-type and trh1/trh2 (Fig 7C, w1118: A2, 176.6 ± 15.0, n = 10; A3, 76.1 ± 4.8, n = 9; A4, 85.8 ± 4.1, n = 10; A5, 66.1 ± 5.4, n = 10; A6, 40.0 ± 3.7, n = 9; trh1/trh2: A2, 173.2 ± 8.8, n = 10; A3, 106.6 ± 15.3, n = 10; A4, 74.1 ± 9.9, n = 9; A5, 60.7 ± 9.1, n = 10; A6, 33.1 ± 3.6, n = 9). Therefore, these analyses suggest that satellite boutons are prone to appear in more anterior than posterior segments.
(A) Images show NMJ 6/7 from segments A2-A6 immunostained for Syn (green) and HRP (red) in w1118 and trh1/ trh2. White arrows indicate clusters of satellite boutons. Scale bar represents 10 μm. (B, C) Line graphs show average percentages (mean ± SEM) of satellite boutons (B) and total bouton numbers (C) for w1118 (magenta lines) and trh1/ trh2 (blue lines). Statistical significance by Mann-Whitney test is shown (n.s., no significance; *, p < 0.05; **, p < 0.01; ***, p < 0.001).
We further examined whether glial Sima and Wg have any role on modifying bouton morphology in the posterior A6 segment of the trh1/trh2 larvae. With the satellite boutons at a basal level in the A6 segment of trh1/trh2 (Fig 7B), we tested whether overexpression of Sima or Wg could induce satellite boutons in the trh1/trh2 mutant. Overexpression of Sima by repo-GAL4 in trh1/trh2 induced some satellite boutons (S5A and S5B Fig, 6.8 ± 2.3%, n = 10), which showed no significant difference to the trh1/trh2 mutant carrying repo-GAL4 (2.8 ± 1.5%, n = 8). Overexpression of Wg by repo-GAL4 in trh1/trh2 displayed a basal level of satellite boutons (2.6 ± 1.5%, n = 9). Also, glial wg-RNAi knockdown suppressed satellite bouton formation in the A3 segment of trh1/trh2 (Fig 4F), but had no effect on the morphological phenotype in the A6 segment (S5C and S5D Fig). Thus, the analysis of these data suggests that Wg and Sima might have relatively specific roles to induce satellite bouton formation in the anterior A3 segment.
Given the satellite bouton phenotype in the trh1/trh2 larvae, we assessed basal synaptic transmission properties, firstly at muscle 6 of the A3 segment. The amplitude of spontaneous release, or the miniature evoked junctional potential (mEJP), was slightly but non-significantly reduced (Fig 8A), from 1.4 ± 0.1 mV (n = 10) in w1118 to 1.2 ± 0.1 mV (n = 13) in trh1/trh2 (Fig 8C). Comparable frequencies were detected between w1118 (1.9 ± 0.1 Hz, n = 10) and trh1/trh2 (2.4 ± 0.3 Hz, n = 13)(Fig 8D). The amplitudes of EJP were also comparable between w1118 (54.6 ± 2.1 mV, n = 10) and trh1/trh2 (50.9 ± 4.2 mV, n = 13)(Fig 8B and 8E). The quantal content, calculated by dividing the EJP amplitude with that of mEJP, was slightly but non-significantly increased, from 40.2 ± 3.3 (n = 10) in w1118 to 46.9 ± 4.5 (n = 13) in trh1/trh2 (Fig 8F). We then evaluated the synaptic transmission properties of muscle 6 for the A6 segment. Between the w1118 control and the trh1/trh2 mutant, the mEJP amplitudes (1.2 ± 0.1, n = 9 v.s. 1.5 ± 0.2 mV, n = 8), the mEJP frequencies (2.0 ± 0.3 Hz, n = 9, v.s. 1.4 ± 0.2 Hz, n = 8), and the EJP amplitudes (46.2 ± 3.7 mV, n = 9 v.s. 41.6 ± 5.7 mV, n = 8) remained similar without significant difference (Fig 8A–8E). However, the slight increase in mEJP and the slight decrease in EJP in trh1/trh2 lead to a significant reduction in the quantal content, from 38.9 ± 3.7 (n = 9) in w1118 to 27.9 ± 1.5 (n = 8) in trh1/trh2 (Fig 8F). Thus, the impaired synaptic activity of the A6 segment in the mutant larvae may underlie the defective stride cycles of the posterior segments.
(A, B) Representative traces of mEJPs (A) and EJPs (B) recorded from NMJs of muscle 6 of segment A3 (left panels) and A6 (right panels) for w1118 and trh1/trh2. The scales for both mEJP and EJP traces are depicted at right. (C-F) Bar graphs show averages (mean ± SEM) for (C) mEJP amplitudes, (D) mEJP frequencies, (E) EJP amplitudes and (F) quantal contents (QC) calculated as EJP/mEJP. Statistical significance by Mann-Whitney test is shown (n.s., no significance and *, p < 0.05).
Discussion
Here, we demonstrate that Trh, a member of the NPAS protein family, non-cell autonomously regulates synaptic bouton formation at NMJs through a hypoxic response from glia. We observed small-sized and clustered satellite boutons at the NMJs of the trh mutant larvae or larvae reared at low oxygen levels. The abnormal bouton morphology at the trh NMJs could be suppressed by reducing the level of the hypoxia-inducible factor Sima in glia. We further show that Sima enhanced the Wg signal from glia to cause satellite bouton formation. Although normal synaptic transmission was detected at NMJs located in an anterior segment of larvae bearing satellite boutons, reduced synaptic transmission was found in a posterior segment lacking satellite boutons of the trh mutant, suggesting that glia-induced satellite bouton formation might be a homeostatic response in restoring normal synaptic transmission. Imbalanced synaptic activities at the anterior and posterior NMJs of the trh mutant might contribute to the uncoordinated stride cycles detected in the trh mutant, slowing larval crawl speed (S6 Fig). Thus, we provide a model for studying the glial responses that modulate synaptic bouton reorganization and activities during hypoxia.
Glia play a critical role in satellite bouton formation in the trh mutant
Our results suggest that Trh has a late developmental role in tracheal morphogenesis, in addition to its well-characterized role in early tracheal cell fate specification [23, 24]. We observed defective tracheal structure in the trh mutant (S2A Fig), which may result in hypoxic conditions inside the larval body. The increase in terminal branch number (S2A and S2B Fig) may be a response to oxygen supply deficiency [12]. Moreover, the increases in Sima protein levels and ODD-GFP reporter expression indicate reduced internal oxygen levels (Fig 2A and 2B and S3E and S3F Fig). Finally, the satellite bouton phenotype in the trh mutant was recapitulated by rearing larvae under hypoxia, and it was suppressed by rearing larvae under hyperoxia. Taken together, these observations suggest that cells in the trh mutant sense low oxygen levels caused by the defective tracheal system and respond by elevating Sima protein levels. It is not clear how profound this effect is for other types of larval cells. Based on our ODD-GFP and Sima immunostaining patterns (Fig 2A and 2B and S3E and S3F Fig), many types of cells are likely to be affected [40].
We suggest that glia is the major cell type mediating satellite bouton formation in the trh mutant under hypoxia. While Sima was increased ubiquitously, manipulating the levels of Sima or Ftg, the negative regulator of Sima, in glia modulates satellite bouton formation (Fig 3 and S3A and S3B Fig). Elevated Sima levels induce tracheal sprouting in tracheal cells, as well as protrusions in non-tracheal cells [12]. Interestingly, we also observed protrusion of glial processes into synaptic area in the trh mutant, indicative of a glial response (Fig 5E and 5F). Several types of cells in Drosophila have been shown to respond to hypoxia [32, 41]. For instance, under hypoxia, elevated Sima levels induce the expression of Breathless (Btl, the FGF receptor) in tracheal cells that branch out seeking cells that express Branchless (Bnl)/FGF, with this latter process also being partially dependent on Sima [12, 13]. In an alternative pathway, atypical soluble guanylyl cyclases can mediate graded and immediate hypoxia responses mainly in neurons [42, 43]. Drosophila glia have not been reported to sense and respond to hypoxia, but mammalian astrocytes in the central nervous system have been shown to be involved in these processes. In a mouse model for middle cerebral artery occlusion, astrocyte activation was shown to play a crucial role in ischemic tolerance, which is mediated through P2X7 receptor-activated HIF-1α upregulation [18]. Under physiological hypoxia, reduced mitochondrial respiration leads to the release of intracellular calcium and exocytosis of ATP-containing vesicles, thereby signaling the brainstem to modulate animal breathing [44]. Our results reveal a role for Drosophila larval glia in sensing hypoxia via the conventional HIF-1α/Sima pathway. We also demonstrate that under hypoxia, glia modulate the formation of synaptic boutons (Fig 3E and 3F). These results clearly place the glia-modulated morphology of synaptic boutons in the context of hypoxia responses.
Sima elevates Wg signal expression in glia in satellite bouton formation
Our study further establishes that in response to hypoxia, Wg is a glial signal that modulates synaptic bouton formation. Two sources of Wg, presynaptic motor neurons and glia, are involved in synaptic growth and remodeling [34, 36]. Our results suggest that Sima upregulates the level of Wg secreted from glia to modulate synapse formation in the trh mutant or in control larvae grown under hypoxia. In hypoxic macrophages, HIF-1α ediates the induction of Wnt11, which is a mammalian homolog of Wg [45]. It is likely wg is a direct target of Sima in Drosophila. We found the HIF1-α binding motif (CGTG) at the -269 nucleotide sequence in the wg promoter. Also, direct binding of Sima to this consensus site was reported in a systematic ChIP-seq experiment [46]. We further show that the level of Wg is controlled by glial Sima in the trh mutant (Fig 4C and 4D) and that Wg mediates Sima-induced satellite bouton formation at trh NMJs (Fig 3B). Also, glial overexpression of Sima upregulates Wg levels at NMJs (S4E and S4F Fig). Wg is also expressed from presynapses [35]. Thus, neuronal wg knockdown partially suppressed satellite bouton phenotypes in the trh mutant (Fig 5A and 5B), while neuronal wg overexpression in wild-type induced the phenotypes (S4A and S4B Fig). These results are consistent with the idea that presynaptic Wg contributes to the overall pool of Wg at NMJs of the trh mutant. As a secreted morphogen, Wg functions in both pre-synaptic and post-synaptic sites [35–37]. At presynaptic terminals, the canonical Wg pathway induces microtubule loop formation to regulate synaptogenesis. We also detected an increase in microtubule loops in the trh mutant (Fig 5C and 5D), consistent with a role for Wg signaling in modulating synaptic reorganization. Postsynaptic Wg signaling leads to subsynaptic reticulum differentiation [35], which was not apparent in the trh mutant (S1D Fig), suggesting that Wg might be a component of the complex hypoxia response that induces synaptic bouton reorganization. Brief exposure to hypoxia induces immature spines and impaired synaptic function in hippocampal neurons [17]. The morphological change to satellite boutons at the A3 segment of the trh NMJs was not accompanied by altered synaptic transmission (Fig 8), which may be compensated during long-term hypoxia.
The satellite boutons, also named as bunch boutons, have been described in spastin mutants [47, 48]. As an AAA ATPase, Spastin severs microtubules to facilitate transport to distal axon segments [49]. Accordingly, the spastin mutant also exhibits a lack of microtubules at terminal boutons [48]. In contrast, the trh mutant presented an increase of the more stabilized microtubule loops (Fig 5C and 5D). Microtubule loops have been linked to synaptic bouton stabilization, and an excess of microtubule loops has been associated with increased synaptic bouton formation [50, 51]. The altered morphology of satellite boutons may be part of the structural changes necessary to maintain normal synaptic transmission under hypoxia. The trh and spastin mutants also exhibit differences in synaptic function, with loss of spastin function slightly enhancing spontaneous synaptic transmission release but reducing evoked synaptic transmission [48]. Thus, although the morphology of synaptic boutons at trh NMJs resembles that of spastin mutants, satellite boutons at trh NMJs retain synaptic functions, unlike the impaired synaptic transmission of spastin mutant boutons.
Differential morphological and physiological changes of anterior and posterior segments in the trh NMJs
The size of NMJs in muscles 6/7 decreases from the anterior to posterior segments, which could represent a coupling with muscle growth [52, 53], thereby maintaining similar electrophysiological efficacy at anterior and posterior NMJs (Fig 8). Interestingly, our findings show that synaptic responses in the trh mutant differ, with satellite boutons only appearing in anterior segments (Fig 7A and 7B). Furthermore, synaptic transmission at trh NMJs remained normal in the anterior A3 segment but was impaired in the posterior A6 segment (Fig 8). These observations are consistent with the idea that satellite bouton formation is a part of a homeostatic response to restore synaptic activity. Why synapses are not reorganized in the posterior segments remains elusive. We failed to detect an upregulation of Wg in the A6 segment (S4C and S4D Fig), and glial Sima overexpression even in the A6 segment of the trh mutant failed to increase satellite boutons significantly (S5A and S5B Fig). Thus, the upregulation of Wg by Sima may be segment-dependent, which awaits further study. Motor neurons in the ventral nerve cord project much longer axons to muscles in posterior segments compared to anterior ones. It has been shown that axonal transport to posterior segments is more vulnerable to inefficient transport conditions. For example, mutation of long-chain Acyl-CoA synthetase impairs the balance between anterograde and retrograde transport, causing distally-biased axonal aggregations and affecting the growth and functioning of synapses [54]. It is possible that glia-derived Wg signals may not be efficiently transported to posterior segments during hypoxia. This polar difference in synaptic activity and bouton morphology may contribute to the uncoordinated movements of the trh mutant larvae. Alternatively, defective locomotion in posterior segments of the trh mutant is independent of the glial modulation of bouton morphological changes. Larval forward locomotion, propelled by peristaltic contraction, is controlled by different circuits targeting anterior and posterior segments. The GABAergic SEZ-LN1 neurons specifically control posterior A6 and A7 segmental muscle contraction by inhibiting A27h premotor neurons, which promotes longitudinal muscle contraction during larval forward crawling [55]. Specific alteration of the circuit in the posterior segments may lead to the locomotion defect in the trh mutant.
Materials and methods
Fly stocks
All flies were reared at 25°C. w1118 was used as wild-type control and to backcross with trh1 or trh2. The sources of fly strains are as follows: trh2, elav-GAL4, MHC-GAL4, repo-GAL4, UAS-trh, UAS-sima, UAS-sima-RNAi, Nrv2-GAL4, and His2Av-mRFP were obtained from Bloomington Drosophila Stock Center (BDSC); trh1, NP6293-GAL4, and UAS-wg-RNAi from Kyoto Stock Center; and UAS-trh-RNAi, UAS-fga-RNAi, and UAS-fz2-RNAi from Vienna Drosophila Resource Center (VDRC). Also used were btl-GAL4[56], moody-GAL4 [57], alrm-GAL4[58], UAS-wg, and GFP-ODD [29]. The repo-cyto-GFP line was generated with the sequence for cytoplasmic GFP under the control of the 4.3 kb repo promoter, which recapitulates the full repo expression pattern.
Hypoxia or hyperoxia rearing conditions
Larvae in a food vial were transferred at 1 day after egg laying (AEL) to a ProOx (model 110, BioSpherix, Lacona, NY) oxygen-controlled chamber. Oxygen or nitrogen was infused into the chamber to a desired concentration (5% or 50%), which was maintained until assay.
Immunostaining
The NMJ phenotypes were analyzed as previously described [59]. For live tissue preparation to detect repo-cyto-GFP expression, non-fixed larvae were dissected and the larval fillets were incubated with anti-horseradish peroxidase (HRP, 1:10) in phosphate buffered saline (PBS) for 10 minutes. Primary antibodies used were against Synapsin (3C11, mouse, 1:100; Developmental Studies Hybridoma Bank, DSHB), HRP-Cy5 (rabbit, 1:100; Jackson ImmunoResearch), Dlg (mouse, 1:100, DSHB), GluRIIA (mouse, 1:100, DSHB), dPAK (rabbit, 1:1000), GluRIII (rabbit, 1:1000), Brp (nc82, mouse, 1:100, DHSB), Sima (guinea pig, 1:1000)[60], Repo (mouse, 1:1000, DSHB), Wg (4D4, mouse, 1:10, DSHB), and Futsch (22C10, mouse, 1:100, DSHB). Secondary antibodies used were anti-rabbit or -mouse 488, Cy3, or Cy5 (1:1000, Jackson ImmunoResearch). For Wg immunostaining, larval preparations of different genotypes were marked and immunostaining was performed in the same test tube.
Image acquisition and processing
NMJs in muscle 6/7 of A3 segments (or A2-A6 in Fig 7 and A6 in S4 and S5 Figs) of wandering third-instar larvae were analyzed. Confocal images were acquired via LSM510 confocal microscopy (Carl Zeiss) using 40x water, 40x water immersion (for live tissue in Fig 5E), or 100x oil objectives. All presented images are projections of confocal z-stacks. Numbers of satellite boutons, total boutons, and microtubule loops were counted manually. The percentage of satellite boutons was calculated as the number of satellite boutons divided by the total boutons (satellite + normal ones) for each NMJ, and the average percentage is calculated from about 6–13 NMJs for each genotype. The immunofluorescence intensities of Wg and HRP were analyzed by ImageJ. HRP-positive regions were chosen to measure mean intensities of Wg and HRP. After subtracting the intensity with the background one, the ratio of Wg levels to HRP levels was presented as the normalized Wg intensity. The overlapping area of GFP and HRP projections was chosen by the “AND” operator in ImageJ, which was divided by HRP area for the percentage. Each dot in the bar graph represents the data from a single NMJ of a larva, and 6–13 NMJs from 2–5 independent experiments were pooled for quantification. Embryos were acquired by means of LSM510 confocal microscopy (Carl Zeiss) using a 20x objective, and were analyzed as previously described [29]. For Fig 2B, each dot in the bar graph represents data from a single embryo in which fluorescence was measured in at least 35 cells.
Electrophysiological recordings
Basal transmission properties were analyzed at NMJs of muscle 6/7 in specified segments of wandering third-instar larvae as previously described [61], with some modifications. The larval body wall was dissected in cold calcium-free HL3 solution and recorded in HL3 solution containing 0.4 mM CaCl2 at room temperature.
Crawling behavior
Mid third instar larvae (feeding stage) were placed on black agar plates (2% agar with black food coloring in 25 × 20 cm2 dishes) at room temperature for filming. Video recording by a Sony Xperia Z1 camera started after 1 min habituation and lasted for 5 min, and it was analyzed using Ctrax software [62]. The (x, y) positions were used to calculate the crawling distance between two successive frames, and crawling speed was derived by dividing total distance travelled by time. The change in angle of larvae between two frames was divided by time to represent rotational angles. The forward crawling assay was a modification of a previous study [39]. Larvae were transferred into a tunnel (~1 mm width) made in 2% black agar. Specimens were video-recorded for 3–10 minutes using a Leica S8 APO microscope. Kymographs were constructed using the MultipleKymograph plug-in for ImageJ (NIH). Only forward crawling was counted, and 7–10 steps for each of ten larvae were analyzed for each genotype.
Supporting information
S1 Fig. Presynaptic and postsynaptic structures in trh mutant.
(A) Bar graphs show averages of muscle areas (w1118, 61.2 ± 1.6 x103 μm2, n = 12; trh1/+, 70.0 ± 2.8 x103 μm2, n = 11; and trh1/trh2, 60.5 ± 4.9 x103 μm2, n = 9). (B, C, D) Images showing NMJs of muscle 6/7 immunostained for Brp, GluRIII, and HRP (B), GluRIIA, dPAK, and HRP (C), and Dlg and HRP (D) in w1118, trh1/+, and trh1/trh2. Scale bars are 10 μm. Statistical significance in (A) was assayed by Mann-Whitney test is shown (n.s., no significance).
https://doi.org/10.1371/journal.pgen.1007980.s001
(TIF)
S2 Fig. Tracheal defects in the trh mutant.
(A) A bright-field view of tracheal dorsal branches of trh1/+ (left top panel) and trh1/trh2 (right top panel), and tracheal dorsal trunks of trh1/trh2 (bottom panels). Numbers denote terminal branches, and arrows indicate a tracheal break (bottom left) and a tracheal tangle (bottom right). Scale bar represents 50 μm. (B) Bar graph shows averages of dorsal terminal branches. (C) Bar graph shows averages of the GFP immunofluorescence intensity (w1118 in 21% O2, 20.3 ± 7.5, n = 6; w1118 in 5% O2, 79.3 ± 18.0, n = 6; trh1/+ in 21% O2, 13.3 ± 3.6, n = 6; trh1/trh2 at 21% O2, 108.2 ± 14.7, n = 6). (D) Bar graph shows averages of RFP immunofluorescence intensity (w1118 in 21% O2, 102.0 ± 15.5, n = 6; w1118 at 5% O2, 109.6 ± 13.1, n = 6; trh1/+ at 21% O2, 143.0 ± 15.7, n = 6; trh1/trh2 at 21% O2, 151.5 ± 16.8, n = 6). Statistical significance by Mann-Whitney test is shown (n.s., no significance; *, p < 0.05; **, p < 0.01; ***, p < 0.001).
https://doi.org/10.1371/journal.pgen.1007980.s002
(TIF)
S3 Fig. Fga in satellite bouton formation, efficiency of sima-RNAi knockdown and Sima upregulation in hypoxia and the trh mutant.
(A) Images showing NMJs from muscle 6/7 immunostained for Syn (green) and HRP (red) in fga-RNAi, and repo-GAL4/fga-RNAi and btl-GAL4/fga-RNAi. (B) Bar graph shows percentages (mean ± SEM) of satellite boutons. Statistical significance was assayed by Mann-Whitney test (n.s., no significance; *, p <0.05). (C, D) RT-PCR for sima mRNA (C) or Western blot for Sima protein (D) expressions in da-GAL4 (control) or da-GAL4/sima-RNAi. The controls Rpl19 mRNA (C) and GAPDH protein (D) are comparable in both genotypes. (E) Images showing central nerve cord (left 3 panels) and peripheral nerve (right 2 panels) immunostained for Sima (green) and Repo (red) for w1118 in 21% O2, w1118 in 5% O2 for 4hrs, and trh1/trh2 in 21% O2. Arrows indicate nuclei of glia with high Sima protein levels. Enlarged images of a nucleus in each condition are shown in left bottom site. (F) Images showing Sima (green) localization in cells of His2Av-mRFP-expressing subperineurial (moody-GAL4), wrapping (Nrv2-GAL4), astrocyte-like (alrm-GAL4) and perineurial (NP6293-GAL4) glia in 5% O2 for 4hrs. Lateral views of the Sima-enriched subtype glia are shown in the left and bottom panels. Scale bar represents 20 μm.
https://doi.org/10.1371/journal.pgen.1007980.s003
(TIF)
S4 Fig. Wg in satellite bouton formation and induction by Sima in glia.
(A) Images showing NMJ 6/7 immunostained for Syn (green) and HRP (red). (B) Bar graph shows percentages (mean ± SEM) of satellite boutons: UAS-wg, 3.1 ± 1.0%, n = 10; repo-GAL4/UAS-wg, 12.7 ± 1.9%, n = 10; elav-GAL4/UAS-wg, 9.6 ± 2.3%, n = 6, and MHC-GAL4/UAS-wg, 4.4 ± 1.4%, n = 10. (C, E) Images showing NMJ 6/7 immunostained for Wg (green) and HRP (red). (D, F) Bar graph shows averages (mean ± SEM) of normalized Wg to HRP intensities. (D) Wg/HRP ratios: trh1/trh2 (A3), 0.35 ± 0.05, n = 7; trh1/trh2 (A6), 0.18 ± 0.04, n = 8; repo-GAL4/wg-RNAi trh1/trh2 (A3), 0.16 ± 0.04, n = 7, elav-GAL4/wg-RNAi trh1/trh2 (A3), 0.25 ± 0.04, n = 6. (F) Wg/HRP ratios: UAS-sima, 0.24 ± 0.06, n = 8; repo-GAL4/UAS-sima, 0.54 ± 0.12, n = 9; elav-GAL4/UAS-sima, 0.20 ± 0.04, n = 7, and MHC-GAL4/UAS-sima, 0.17 ± 0.03, n = 7. Statistical significance was assayed by Mann-Whitney test (n.s., no significance; *, p < 0.05, ***, p < 0.001).
https://doi.org/10.1371/journal.pgen.1007980.s004
(TIF)
S5 Fig. Wg and Sima have minor effects on satellite bouton formation in A6 segment of the trh mutant.
(A, C) Images showing NMJs of muscle 6/7 immunostained for Syn (green) and HRP (red) in the A6 segment. (B, D) Bar graphs show percentages (mean ± SEM) of satellite boutons in the A6 segment of (B) repo-GAL4 trh1/trh2, repo-GAL4/UAS-sima trh1/trh2, and repo-GAL4/UAS-wg trh1/trh2; and (D) wg-RNAi trh1/trh2, 2.9 ± 1.7%, n = 10 and repo-GAL4/wg-RNAi trh1/trh2, 1.5 ± 1.0%, n = 10. Statistical significance by Mann-Whitney test is shown (n.s., no significance).
https://doi.org/10.1371/journal.pgen.1007980.s005
(TIF)
S6 Fig. Schematic model for satellite bouton induction and crawling defects in the trh mutant.
(A) In wild-type larvae with normal oxygen supply, the low basal Sima level in glia induces basal Wg expression and secretion to regulate normal synaptic bouton formation. In the trh mutant with a defect tracheal system and limited oxygen supply, Sima is upregulated to induce Wg expression and secretion, thus induce satellite bouton formation in the anterior segments of larvae. (B) Wild-type larvae crawl normally with normal sets of synaptic boutons at NMJs of A3 and A6 segments. In the trh mutant, the appearance of satellite boutons in the anterior segments (e.g. A3) may be a part of homeostasis for restoring normal synaptic activity while the posterior segments (e.g. A6) were unable to compensate for the reduction of synaptic activity, leading to uncoordinated peristaltic movement.
https://doi.org/10.1371/journal.pgen.1007980.s006
(TIF)
S1 Movie. Forward crawling of w1118.
The movie shows freely crawling of w1118 larva in mid-third instar stage for 8.6 seconds with 1x time speed.
https://doi.org/10.1371/journal.pgen.1007980.s007
(AVI)
S2 Movie. Forward crawling of trh1/trh2.
The movie shows freely crawling of trh1/trh2 larva in mid-third instar stage for 7.3 seconds with 1x time speed.
https://doi.org/10.1371/journal.pgen.1007980.s008
(AVI)
S1 Dataset. Raw data for Figs 1–8 and S1–S5 Figs.
https://doi.org/10.1371/journal.pgen.1007980.s009
(XLSX)
Acknowledgments
We thank Y. H. Sun (Academia Sinica), S. Luschnig (Universität Münster), T Leung (National University of Singapore), A. DiAntonio (Washington University), BDSC, Kyoto Stock Center, VDRC, DSHB, and Taiwan Fly Core, for providing reagents; IMB Bioinformatics-Biology Service Core and NPAS Neuro-electrophysiology Core for technical support; as well as H. Li, V. Nithianandam and all members of C.-T. Chien’s laboratory for discussion and comments.
References
- 1. Tam SJ, Watts RJ. Connecting vascular and nervous system development: angiogenesis and the blood-brain barrier. Annu Rev Neurosci. 2010;33:379–408. pmid:20367445
- 2. Banerjee S, Bhat MA. Neuron-glial interactions in blood-brain barrier formation. Annu Rev Neurosci. 2007;30:235–258. pmid:17506642
- 3. Iadecola C. Neurovascular regulation in the normal brain and in Alzheimer's disease. Nat Rev Neurosci. 2004;5(5):347. pmid:15100718
- 4. Zacchigna S, Lambrechts D, Carmeliet P. Neurovascular signalling defects in neurodegeneration. Nat Rev Neurosci. 2008;9(3):169. pmid:18253131
- 5.
Lee DL, Atkinson HJ. Physiology of nematodes: Macmillan Press; 1976.
- 6.
Manning G, Krasnow MA. Development of the Drosophila tracheal system. The development of Drosophila melanogaster. 1. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press; 1993. p. 609–685.
- 7. Danjo R, Kawasaki F, Ordway RW. A tripartite synapse model in Drosophila. PloS one. 2011;6(2):e17131. pmid:21359186
- 8. Semenza GL, Nejfelt MK, Chi SM, Antonarakis SE. Hypoxia-inducible nuclear factors bind to an enhancer element located 3'to the human erythropoietin gene. Proc Natl Acad Sci U S A. 1991;88(13):5680–5684. pmid:2062846
- 9. Wang GL, Semenza GL. Purification and characterization of hypoxia-inducible factor 1. J Biol Chem. 1995;270(3):1230–1237. pmid:7836384
- 10. Stroka DM, Burkhardt T, Desbaillets I, Wenger RH, Neil DA, Bauer C, et al. HIF-1 is expressed in normoxic tissue and displays an organ-specific regulation under systemic hypoxia. FASEB J. 2001;15(13):2445–2453. pmid:11689469
- 11. Jaakkola P, Mole DR, Tian Y-M, Wilson MI, Gielbert J, Gaskell SJ, et al. Targeting of HIF-α to the von Hippel-Lindau ubiquitylation complex by O2-regulated prolyl hydroxylation. Science. 2001;292(5516):468–472. pmid:11292861
- 12. Centanin L, Dekanty A, Romero N, Irisarri M, Gorr TA, Wappner P. Cell autonomy of HIF effects in Drosophila: tracheal cells sense hypoxia and induce terminal branch sprouting. Dev Cell. 2008;14(4):547–558. pmid:18410730
- 13. Jarecki J, Johnson E, Krasnow MA. Oxygen regulation of airway branching in Drosophila is mediated by branchless FGF. Cell. 1999;99(2):211–220. pmid:10535739
- 14. Ortega JA, Sirois CL, Memi F, Glidden N, Zecevic N. Oxygen levels regulate the development of human cortical radial glia cells. Cereb Cortex. 2016;27(7):3736–3751.
- 15. O’driscoll C, Gorman A. Hypoxia induces neurite outgrowth in PC12 cells that is mediated through adenosine A2A receptors. Neuroscience. 2005;131(2):321–329. pmid:15708476
- 16. Nikonenko I, Jourdain P, Muller D. Presynaptic remodeling contributes to activity-dependent synaptogenesis. J Neurosci. 2003;23(24):8498–8505. pmid:13679418
- 17. Segura I, Lange C, Knevels E, Moskalyuk A, Pulizzi R, Eelen G, et al. The oxygen sensor PHD2 controls dendritic spines and synapses via modification of filamin A. Cell Rep. 2016;14(11):2653–2667. pmid:26972007
- 18. Hirayama Y, Ikeda-Matsuo Y, Notomi S, Enaida H, Kinouchi H, Koizumi S. Astrocyte-mediated ischemic tolerance. J Neurosci. 2015;35(9):3794–3805. pmid:25740510
- 19. Stanco A, Pla R, Vogt D, Chen Y, Mandal S, Walker J, et al. NPAS1 represses the generation of specific subtypes of cortical interneurons. Neuron. 2014;84(5):940–953. pmid:25467980
- 20. Pieper AA, Wu X, Han TW, Estill SJ, Dang Q, Wu LC, et al. The neuronal PAS domain protein 3 transcription factor controls FGF-mediated adult hippocampal neurogenesis in mice. Proc Natl Acad Sci U S A. 2005;102(39):14052–14057. pmid:16172381
- 21. Zhou S, Degan S, Potts EN, Foster WM, Sunday ME. NPAS3 is a trachealess homolog critical for lung development and homeostasis. Proc Natl Acad Sci U S A. 2009;106(28):11691–11696. pmid:19581591
- 22. Levesque BM, Zhou S, Shan L, Johnston P, Kong Y, Degan S, et al. NPAS1 regulates branching morphogenesis in embryonic lung. Am J Respir Cell Mol Biol. 2007;36(4):427–434. pmid:17110583
- 23. Chung S, Chavez C, Andrew DJ. Trachealess (Trh) regulates all tracheal genes during Drosophila embryogenesis. Dev Biol. 2011;360(1):160–172. pmid:21963537
- 24. Wilk R, Weizman I, Shilo B-Z. trachealess encodes a bHLH-PAS protein that is an inducer of tracheal cell fates in Drosophila. Genes Dev. 1996;10(1):93–102. pmid:8557198
- 25. Jurgens G, Wieschaus E, Nusslein-Volhard C, Kluding H. Mutations affecting the pattern of the larval cuticle in Drosophila melanogaster. II. Zygotic loci on the third chromosome. Wilehm Roux Arch Dev Biol. 1984.
- 26. Younossi-Hartenstein A, Hartenstein V. The role of the tracheae and musculature during pathfinding of Drosophila embryonic sensory axons. Dev Biol. 1993;158(2):430–447. pmid:8344461
- 27. Menon KP, Carrillo RA, Zinn K. Development and plasticity of the Drosophila larval neuromuscular junction. Wiley Interdiscip Rev Dev Biol. 2013;2(5):647–670. pmid:24014452
- 28. Isaac DD, Andrew DJ. Tubulogenesis in Drosophila: a requirement for the trachealess gene product. Genes Dev. 1996;10(1):103–117. pmid:8557189
- 29. Misra T, Baccino-Calace M, Meyenhofer F, Rodriguez-Crespo D, Akarsu H, Armenta-Calderón R, et al. A genetically encoded biosensor for visualising hypoxia responses in vivo. Biol Open. 2017;6(2):296–304. pmid:28011628
- 30. Centanin L, Gorr TA, Wappner P. Tracheal remodelling in response to hypoxia. J Insect Physiol. 2010;56(5):447–454. pmid:19482033
- 31. Wang GL, Jiang B-H, Rue EA, Semenza GL. Hypoxia-inducible factor 1 is a basic-helix-loop-helix-PAS heterodimer regulated by cellular O2 tension. Proc Natl Acad Sci U S A. 1995;92(12):5510–5514. pmid:7539918
- 32. Lavista-Llanos S, Centanin L, Irisarri M, Russo DM, Gleadle JM, Bocca SN, et al. Control of the hypoxic response in Drosophila melanogaster by the basic helix-loop-helix PAS protein similar. Mol Cell Biol. 2002;22(19):6842–6853. pmid:12215541
- 33. Bruick RK, McKnight SL. A conserved family of prolyl-4-hydroxylases that modify HIF. Science. 2001;294(5545):1337–1340. pmid:11598268
- 34. Kerr KS, Fuentes-Medel Y, Brewer C, Barria R, Ashley J, Abruzzi KC, et al. Glial wingless/Wnt regulates glutamate receptor clustering and synaptic physiology at the Drosophila neuromuscular junction. J Neurosci. 2014;34(8):2910–2920. pmid:24553932
- 35. Packard M, Koo ES, Gorczyca M, Sharpe J, Cumberledge S, Budnik V. The Drosophila Wnt, wingless, provides an essential signal for pre-and postsynaptic differentiation. Cell. 2002;111(3):319–330. pmid:12419243
- 36. Miech C, Pauer H-U, He X, Schwarz TL. Presynaptic local signaling by a canonical wingless pathway regulates development of the Drosophila neuromuscular junction. J Neurosci. 2008;28(43):10875–10884. pmid:18945895
- 37. Mathew D, Ataman B, Chen J, Zhang Y, Cumberledge S, Budnik V. Wingless signaling at synapses is through cleavage and nuclear import of receptor DFrizzled2. Science. 2005;310(5752):1344–1347. pmid:16311339
- 38. Brink DL, Gilbert M, Xie X, Petley-Ragan L, Auld VJ. Glial processes at the Drosophila larval neuromuscular junction match synaptic growth. PloS one. 2012;7(5):e37876. pmid:22666403
- 39. Heckscher ES, Lockery SR, Doe CQ. Characterization of Drosophila larval crawling at the level of organism, segment, and somatic body wall musculature. J Neurosci. 2012;32(36):12460–12471. pmid:22956837
- 40. Ke Q, Costa M. Hypoxia-inducible factor-1 (HIF-1). Mol Pharmacol. 2006;70(5):1469–1480. pmid:16887934
- 41. Arquier N, Vigne P, Duplan E, Hsu T, Therond PP, Frelin C, et al. Analysis of the hypoxia-sensing pathway in Drosophila melanogaster. Biochem J. 2006;393(2):471–480.
- 42. Morton DB, Stewart JA, Langlais KK, Clemens-Grisham RA, Vermehren A. Synaptic transmission in neurons that express the Drosophila atypical soluble guanylyl cyclases, Gyc-89Da and Gyc-89Db, is necessary for the successful completion of larval and adult ecdysis. J Exp Biol. 2008;211(10):1645–1656.
- 43. Vermehren-Schmaedick A, Ainsley JA, Johnson WA, Davies S-A, Morton DB. Behavioral responses to hypoxia in Drosophila larvae are mediated by atypical soluble guanylyl cyclases. Genetics. 2010;186(1):183–196. pmid:20592263
- 44. Angelova PR, Kasymov V, Christie I, Sheikhbahaei S, Turovsky E, Marina N, et al. Functional oxygen sensitivity of astrocytes. J Neurosci. 2015;35(29):10460–10473. pmid:26203141
- 45. Mori H, Yao Y, Learman BS, Kurozumi K, Ishida J, Ramakrishnan SK, et al. Induction of WNT11 by hypoxia and hypoxia-inducible factor-1α regulates cell proliferation, migration and invasion. Sci Rep. 2016;6:21520. pmid:26861754
- 46. Davis CA, Hitz BC, Sloan CA, Chan ET, Davidson JM, Gabdank I, et al. The Encyclopedia of DNA elements (ENCODE): data portal update. Nucleic Acids Res. 2017;46(D1):D794–D801.
- 47. Ozdowski EF, Gayle S, Bao H, Zhang B, Sherwood NT. Loss of Drosophila melanogaster p21-activated kinase 3 suppresses defects in synapse structure and function caused by spastin mutations. Genetics. 2011;189(1):123–135. pmid:21705760
- 48. Sherwood NT, Sun Q, Xue M, Zhang B, Zinn K. Drosophila spastin regulates synaptic microtubule networks and is required for normal motor function. PLoS Biol. 2004;2(12):e429. pmid:15562320
- 49. Karabay A, Yu W, Solowska JM, Baird DH, Baas PW. Axonal growth is sensitive to the levels of katanin, a protein that severs microtubules. J Neurosci. 2004;24(25):5778–5788. pmid:15215300
- 50. Roos J, Hummel T, Ng N, Klämbt C, Davis GW. Drosophila Futsch regulates synaptic microtubule organization and is necessary for synaptic growth. Neuron. 2000;26(2):371–382. pmid:10839356
- 51. Franco B, Bogdanik L, Bobinnec Y, Debec A, Bockaert J, Parmentier M-L, et al. Shaggy, the homolog of glycogen synthase kinase 3, controls neuromuscular junction growth in Drosophila. J Neurosci. 2004;24(29):6573–6577. pmid:15269269
- 52. Koch I, Schwarz H, Beuchle D, Goellner B, Langegger M, Aberle H. Drosophila ankyrin 2 is required for synaptic stability. Neuron. 2008;58(2):210–222. pmid:18439406
- 53. Nijhof B, Castells-Nobau A, Wolf L, Scheffer-de Gooyert JM, Monedero I, Torroja L, et al. A new Fiji-based algorithm that systematically quantifies nine synaptic parameters provides insights into Drosophila NMJ morphometry. PLoS Comput Biol. 2016;12(3):e1004823. pmid:26998933
- 54. Liu Z, Huang Y, Zhang Y, Chen D, Zhang YQ. Drosophila Acyl-CoA synthetase long-chain family member 4 regulates axonal transport of synaptic vesicles and is required for synaptic development and transmission. J Neurosci. 2011;31(6):2052–2063. pmid:21307243
- 55. Tastekin I, Khandelwal A, Tadres D, Fessner ND, Truman JW, Zlatic M, et al. Sensorimotor pathway controlling stopping behavior during chemotaxis in the Drosophila melanogaster larva. Elife. 2018;7:e38740. pmid:30465650
- 56. Shiga Y, Tanaka‐Matakatsu M, Hayashi S. A nuclear GFP/β‐galactosidase fusion protein as a marker for morphogenesis in living Drosophila. Dev Growth Differ. 1996;38(1):99–106.
- 57. Bainton RJ, Tsai LT-Y, Schwabe T, DeSalvo M, Gaul U, Heberlein U. moody encodes two GPCRs that regulate cocaine behaviors and blood-brain barrier permeability in Drosophila. Cell. 2005;123(1):145–156. pmid:16213219
- 58. Stork T, Sheehan A, Tasdemir-Yilmaz OE, Freeman MR. Neuron-glia interactions through the Heartless FGF receptor signaling pathway mediate morphogenesis of Drosophila astrocytes. Neuron. 2014;83(2):388–403. pmid:25033182
- 59. Wang M, Chen P-Y, Wang C-H, Lai T-T, Tsai P-I, Cheng Y-J, et al. Dbo/Henji modulates synaptic dPAK to gate glutamate receptor abundance and postsynaptic response. PLoS Genet. 2016;12(10):e1006362. pmid:27736876
- 60. Wang C-W, Purkayastha A, Jones KT, Thaker SK, Banerjee U. In vivo genetic dissection of tumor growth and the Warburg effect. Elife. 2016;5:e18126. pmid:27585295
- 61. Tsai P-I, Kao H-H, Grabbe C, Lee Y-T, Ghose A, Lai T-T, et al. Fak56 functions downstream of integrin alphaPS3betanu and suppresses MAPK activation in neuromuscular junction growth. Neural Dev. 2008;3(1):26.
- 62. Branson K, Robie AA, Bender J, Perona P, Dickinson MH. High-throughput ethomics in large groups of Drosophila. Nat Methods. 2009;6(6):451. pmid:19412169