# The Amino Acid Transporter JhI-21 Coevolves with Glutamate Receptors, Impacts NMJ Physiology, and Influences Locomotor Activity in [i]Drosophila[i] Larvae

Anna Ziegler, Hrvoje Augustin, Nathan L. Clark, Martine Berthelot-Grosjean, Mégane Simonnet, Joern R. Steinert, Flore Geillon, Gérard Maniere, David E. Featherstone, Yaël Grosjean

► To cite this version:

Anna Ziegler, Hrvoje Augustin, Nathan L. Clark, Martine Berthelot-Grosjean, Mégane Simonnet, et al. The Amino Acid Transporter JhI-21 Coevolves with Glutamate Receptors, Impacts NMJ Physiology, and Influences Locomotor Activity in [i]Drosophila[i] Larvae. Scientific Reports, 2016, 6, pp.19692. 10.1038/srep19692. hal-01277627

HAL Id: hal-01277627

https://hal.science/hal-01277627v1

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# The Amino Acid Transporter Jhl-21 Coevolves with Glutamate Receptors, Impacts NMJ Physiology, and Influences Locomotor Activity in *Drosophila* Larvae

Received: 24 July 2015

Accepted: 16 December 2015

Published: 25 January 2016

Anna B. Ziegler<fnref n="1" />,²,³, Hrvoje Augustin⁴,†, Nathan L. Clark⁵, Martine Berthelot-Grosjean<fnref n="1" />,²,³, Mégane M. Simonnet<fnref n="1" />,²,³, Joern R. Steinert⁶, Flore Geillon<fnref n="1" />,²,³, Gérard Manière<fnref n="1" />,²,³, David E. Featherstone⁴ & Yael Grosjean<fnref n="1" />,²,³

Changes in synaptic physiology underlie neuronal network plasticity and behavioral phenomena, which are adjusted during development. The *Drosophila* larval glutamatergic neuromuscular junction (NMJ) represents a powerful synaptic model to investigate factors impacting these processes. Amino acids such as glutamate have been shown to regulate *Drosophila* NMJ physiology by modulating the clustering of postsynaptic glutamate receptors and thereby regulating the strength of signal transmission from the motor neuron to the muscle cell. To identify amino acid transporters impacting glutamatergic signal transmission, we used Evolutionary Rate Covariation (ERC), a recently developed bioinformatic tool. Our screen identified ten proteins co-evolving with NMJ glutamate receptors. We selected one candidate transporter, the SLC7 (Solute Carrier) transporter family member Jhl-21 (Juvenile hormone Inducible-21), which is expressed in *Drosophila* larval motor neurons. We show that Jhl-21 suppresses postsynaptic muscle glutamate receptor abundance, and that Jhl-21 expression in motor neurons regulates larval crawling behavior in a developmental stage-specific manner.

The glutamatergic *Drosophila melanogaster* larval neuromuscular junction (NMJ) is a powerful well-established model for the study of synaptic development and function. During the three larval stages the morphology of the NMJ changes dramatically<fnref n="1" />,². From the hatching of a *Drosophila* larva up to the last larval instar, the muscle surface increases faster than the growth of the nerve terminals that innervate it. Despite this, the strength of these synapses is maintained at the same level³. This means that during larval development either the amount of released neurotransmitter or the receptivity of the muscle cell have to be adjusted. This could be achieved via a variety of mechanisms, including addition of new synapses to each junction, and changes in strength of individual synapses.

NMJ strength can also be tuned in previously unsuspected ways. In a previous study, for example, we identified a glial amino acid exchanger, Genderblind (GB), which is capable of tuning synaptic strength by regulating the amount of extracellular glutamate. This glutamate constitutively desensitizes ionotropic glutamate receptors (iGluRs), inhibits their clustering, and thereby suppresses synaptic transmission⁴,⁵. We also showed that the

<footnote n="1">¹CNRS, UMR6265 Centre des Sciences du Goût et de l'Alimentation, F-21000 Dijon, France. ²INRA, UMR1324 Centre des Sciences du Goût et de l'Alimentation, F-21000 Dijon, France. ³Université de Bourgogne Franche-Comté, UMR Centre des Sciences du Goût et de l'Alimentation, F-21000 Dijon, France. ⁴Biological Sciences, University of Illinois at Chicago, Illinois 60607, Chicago, USA. ⁵Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh 15260, Pennsylvania, USA. ⁶MRC Toxicology Unit, University of Leicester, LE1 9HN Leicester, UK. ⁷Present address: Institute of Healthy Ageing, and GEE, University College London, WC1E 6BT London, UK. Correspondence and requests for materials should be addressed to Y.G. (email: yael.grosjean@u-bourgogne.fr)</footnote>
<figure><img src="image_3.png" /></figure>

b
Flow chart:
screening for transporters involved in gutamateric signaling using ERC

<figure><img src="image_4.png" /><figcaption>**Figure 1. Evolutionary rate covariation.** (a) The rates of evolution used in this study describe changes in protein sequences over time. To study if proteins are co-evolving, species trees were generated. In our study those trees were made by using homologue proteins of the following species: *D.melanogaster*, *D. simulans*, *D. sechellia*, *D. yakuba*, *D. erecta*, *D. ananassae*, *D. pseudoobscura*, *D. persimillis*, *D. wilistoni*, *D. mojavensis*, *D. virillis* and, *D. grimshavi*. The rates of a pair of proteins are plotted against each other and a factor is calculated describing the correlation of these rates. The rates of co-evolving proteins (protein A and B) are positive while proteins, which are not co-evolving (protein B and C) have correlation coefficients close to zero or negative. (b) As a positive control we calculated the global ERC value for the six glutamate receptors (mGluRA and iGluRA-E), which are known to act together at the *Drosophila* NMJ. For statistical analysis a mean correlation was calculated and compared to the mean correlation of random sets of six proteins. Next ERC values were calculated between the 39 transporter candidates and the GluRs mentioned above. Ten transporter candidates were showing robust ERC values (p < 0.05) with each glutamate receptor. Next we compared the expression pattern of the six GluRs with the expression pattern of the putative transporters. mGluR is predominantly expressed in the CNS, iGluRs in the carcass. Twelve putative transporters were showing their highest level of expression in either the CNS or the body wall (www.flybase.org), JhI-21 showed robust ERC with GluRs and was expressed highest in the same tissue than GluRs, and was chosen for further investigation. For statistical analysis the mean correlation of JhI-21 and the six GluRs was compared to a random groups of 7 proteins.</figcaption></figure>

presynaptic neuron is capable of secreting non-vesicular glutamate through an unknown transporter to regulate NMJ strength by modulating iGluR clustering⁶.

In order to identify amino acid transporters that might regulate synaptic physiology during development, we used Evolutionary Rate Covariation (ERC). ERC is a recently established bioinformatic method that identifies functional relationships between proteins based on their evolutionary histories. The hypothesis of ERC is that functionally related proteins experience similar evolutionary selective pressures and hence have rates of evolution that correlate across species. ERC values are calculated by generating phylogenetic trees using full protein sequences and computing the correlation between the rates of change of two proteins across the branches of a phylogeny. The resulting values could range from 1 in case of positive correlation to -1 in case of negative correlation (Fig. 1a)⁷. ERC has previously been used to study proteins that are physically interacting or present in the same protein complex⁸⁻¹². However, recent studies demonstrated that functionally related and coexpressed genes reveal positive and significant ERC values as well⁷. In this study, we screened for transporters showing
<table><caption>**Table 1.** Glutamate receptors show overall positive ERC values and highest expression in the CNS or the body wall. The table shows the grouped ERC value of all six glutamate receptors together and its corresponding *p*-value. The positive ERC value indicates the strength of evolutionary covariation between all 6 GluRs. The larval expression pattern of the glutamate receptors according to Flyatlas (www.flybase.org) is summarized as followed: X = low expression (10 – 99), XX = moderate expression (100 – 499). Tissues are abbreviated as followed: central nervous system (CNS), digestive system (DS), malpighian tubules (MT).</caption><thead><tr><th rowspan="2">Gene</th><th colspan="2">ERC Data</th><th colspan="6">GluR expression pattern</th></tr><tr><th>ERC value</th><th>p-value</th><th>CNS</th><th>Body wall</th><th>DS</th><th>MT</th><th>Fat body</th><th>Trachea</th></tr></thead><tbody><tr><td>mGluR</td><td rowspan="6">0.556</td><td rowspan="6">0.00021</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>GluRIIA</td><td></td><td>XX</td><td>X</td><td>X</td><td></td><td></td></tr><tr><td>GluRIIB</td><td>X</td><td>X</td><td>X</td><td>X</td><td></td><td></td></tr><tr><td>GluRIIC</td><td></td><td>XX</td><td>X</td><td></td><td></td><td></td></tr><tr><td>GluRIID</td><td></td><td>XX</td><td>X</td><td></td><td></td><td></td></tr><tr><td>GluRIIE</td><td></td><td>XX</td><td>X</td><td></td><td></td><td></td></tr></tbody></table>

evolutionary covariation with six glutamate receptor subunits. We hypothesized that the co-evolution of amino acid transporters and glutamate receptors might lead to the identification of unknown genes involved in glutamatergic signaling. We subsequently tested the functional relationship between those six GluRs and the transporters by taking advantage of the very well studied *Drosophila* larval NMJ physiology. Six GluRs have been shown to impact synaptic strength at this synapse: a metabotropic glutamate receptor subunit (mGluRA) expressed in motor neurons, and five iGluR subunits (GluRIIA, GluRIIB, GluRIIC, GluRIID, and GluRIIE) forming the ionotropic A- and B-type receptors expressed by the post-synaptic muscle cell¹³,¹⁴ We were particularly interested in transporters co-expressed with either mGluRA in motor-neurons or with iGluR subunits in muscle cells.

Using these criteria, we found a strong ERC value between these glutamate receptor subunits and the amino acid transporter JhI-21. Consistent with the idea that ERC predicts functional relationships, we found that JhI-21 negatively regulates iGluR clustering at NMJs and plays a role in locomotion control during late larval development. Investigating the reason for these effects, we discovered differential expression of JhI-21 in the central nervous system neurons and at the NMJ during larval feeding and wandering stages. Taken together, our results demonstrate how ERC can be used to find novel previously unsuspected roles for proteins, and reveal for the first time a role for JhI-21 in glutamatergic synapse function and behavior.

## Results

**JhI-21 is identified as a component of the glutamatergic signaling pathway by Evolutionary Rate Covariation (ERC).** Six glutamate receptor subunits (mGluRA, GluRIIA-E) are known to be key members for signal transmission at the glutamatergic 3rd-instar larval NMJ. We identified their orthologs in 12 *Drosophila* species and calculated their ERC values (Fig. 1b). The six glutamate receptors showed overall positive scores indicating robust rate covariation. The mean ERC value between all possible pairs of those 6 proteins is 0.556 and is strongly significant (*p* = 0.00021, permutation test; Table 1). Positive ERC values are typically found for proteins acting as subunits in the same complex, as in the case for the ionotropic glutamate receptor subunits GluRIIA-E. Glutamate receptors in larval NMJs localize to postsynaptic densities on the surface of the muscle cell and always contain GluRIIC, GluRIID, and GluRIIE subunits. In addition, “A-Type” receptors contain GluRIIA while “B-Type” receptors contain GluRIIB. In contrast, the metabotropic mGluRA receptor subunit is expressed in the innervating motor neuron, and acts in a glutamatergic feedback loop regulating the amount of transmitter released¹⁵. Even though metabotropic and ionotropic glutamate receptors are not co-localized within the same cell at the NMJ of 3rd-instar larvae, they do act in the same intercellular signaling pathway and therefore the ERC values between the two different types of receptors are very high. These conditions demonstrate that ERC can detect proteins that function together despite being expressed in different cells. This finding enabled us to use this technique to screen for previously unrecognized candidate amino acid transporters associated with glutamatergic signal transmission.

In total, we first screened 39 confirmed and putative amino acid transporter homologs¹⁶,¹⁷ for rate covariation with the six glutamate receptors mentioned above (Table 2 and Supplemental Table1). Ten putative transporters showed positive and significant (*p* < 0.05) ERC values with the six glutamate receptor subunits (Table 2, Fig. 2). To prove the functional relationship between these GluRs and our positive hits we decided to focus next on the glutamatergic larval NMJ. The cell bodies of motor-neurons expressing mGluRA are located in the ventral nerve cord of the larval CNS. Those neurons project their axons towards the body wall where they contact muscle cells, which express iGluRs. Therefore, in parallel we used a second screening procedure to determine which transporters are most highly expressed in either the CNS (where motor neuron cell bodies are housed) or the body wall (location of NMJs) (Tables 1 and 2). We found eleven putative transporter genes (Table 2, Fig. 2). We selected one significant hit from our ERC analysis, JhI-21, which also matched very well our second criterion concerning the expression pattern in the neuromuscular system (Fig. 2).

The single ERC values of JhI-21 compared to the different glutamatergic receptors range from 0.298 to 0.59 (Supplemental Table 1). The global ERC value for JhI-21 and the six glutamate receptor subunits that we examined is 0.3788 (*p* = 0.0019, permutation test) (Table 2). This very high co-evolutionary statistic suggested that JhI-21 and the glutamate receptor subunits co-evolved and therefore may be functionally related.
<table><caption>**Table 2. Co-evolution of putative transporter genes with glutamate receptors and summary of their larval expression pattern.** The table shows group ERC values of 39 putative amino acid transporters and glutamate receptor subunits (mGluRA, GluRIIA, GluRIIB, GluRIIC, GluRIID, and GluRIIE) with their corresponding *p*-values. The magnitude of positive ERC values indicate the strength of evolutionary covariation between the amino acid transporter and the 6 GluRs. The expression pattern according to Fly Atlas (www.flybase.org) is summarized as followed: X = low expression (10 – 99), XX = moderate expression (100 – 499), XXX = high level expression (500 – 999), XXXX = very high expression (>999). Tissues are abbreviated as followed: central nervous system (CNS), digestive system (DS), malpighian tubules (MT).</caption><thead><tr><th rowspan="2">Transporter</th><th colspan="2">ERC Data</th><th colspan="6">Transporter expression pattern</th></tr><tr><th>ERC value</th><th>p-value</th><th>CNS</th><th>Body wall</th><th>DS</th><th>MT</th><th>Fat body</th><th>Trachea</th></tr></thead><tbody><tr><td>bdg</td><td>0.465</td><td>0.0001</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>CG8785</td><td>0.3828</td><td>0.0007</td><td></td><td></td><td>XXXX</td><td>XX</td><td></td><td>X</td></tr><tr><td>CG13384</td><td>0.4405</td><td>0.0007</td><td>XX</td><td>XX</td><td>XX</td><td>XX</td><td>XXX</td><td>XX</td></tr><tr><td>Jhl-21</td><td>0.3788</td><td>0.0019</td><td>XXX</td><td>XX</td><td>XX</td><td>X</td><td>X</td><td>XX</td></tr><tr><td>CG16700</td><td>0.4073</td><td>0.0023</td><td>X</td><td>XX</td><td>XX</td><td>X</td><td></td><td>X</td></tr><tr><td>CG9413</td><td>0.305</td><td>0.0028</td><td>XX</td><td>XXX</td><td>XXXX</td><td>X</td><td></td><td>X</td></tr><tr><td>CG7255</td><td>0.4127</td><td>0.0045</td><td></td><td>X</td><td>X</td><td>X</td><td>X</td><td></td></tr><tr><td>DAT</td><td>0.3223</td><td>0.0075</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>mnd</td><td>0.2603</td><td>0.0122</td><td>XX</td><td>XX</td><td>X</td><td>XX</td><td>X</td><td>X</td></tr><tr><td>CG13795</td><td>0.23</td><td>0.0467</td><td></td><td>X</td><td>X</td><td>X</td><td>XXXX</td><td></td></tr><tr><td>blot</td><td>0.183</td><td>0.0624</td><td>XX</td><td>XXX</td><td>XX</td><td>XX</td><td>X</td><td>XXX</td></tr><tr><td>ine</td><td>0.1678</td><td>0.0737</td><td>X</td><td>X</td><td>XXXX</td><td>XX</td><td></td><td></td></tr><tr><td>gb</td><td>0.1493</td><td>0.1334</td><td>XXX</td><td>X</td><td>XX</td><td>X</td><td>XXX</td><td>X</td></tr><tr><td>CG1698</td><td>0.1077</td><td>0.2114</td><td></td><td>X</td><td>XXX</td><td>X</td><td></td><td>X</td></tr><tr><td>CG32079</td><td>0.107</td><td>0.2472</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VGlut</td><td>0.0648</td><td>0.2617</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CG43066</td><td>0.0193</td><td>0.4055</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SerT</td><td>0.0298</td><td>0.4108</td><td>X</td><td></td><td>X</td><td></td><td></td><td></td></tr><tr><td>CG5535</td><td>0.0737</td><td>0.4433</td><td>XX</td><td>XX</td><td>XX</td><td>X</td><td>X</td><td>XX</td></tr><tr><td>CG4476</td><td>-0.0002</td><td>0.597</td><td>XX</td><td></td><td>X</td><td></td><td></td><td></td></tr><tr><td>CG15279</td><td>-0.0185</td><td>0.6095</td><td></td><td>XX</td><td>X</td><td>XXXX</td><td></td><td>XX</td></tr><tr><td>CG13796</td><td>-0.0432</td><td>0.6126</td><td>X</td><td>XX</td><td>XX</td><td>X</td><td>XXXX</td><td>XX</td></tr><tr><td>CG4991</td><td>-0.063</td><td>0.7208</td><td></td><td>X</td><td>XX</td><td></td><td></td><td>X</td></tr><tr><td>CG12531</td><td>-0.072</td><td>0.600</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>NAAT1</td><td>-0.0828</td><td>0.739</td><td></td><td></td><td>XX</td><td>XX</td><td></td><td></td></tr><tr><td>Eaat1</td><td>-0.087</td><td>0.74</td><td>XXX</td><td></td><td>XX</td><td>X</td><td>XX</td><td></td></tr><tr><td>path</td><td>-0.0598</td><td>0.797</td><td>XXX</td><td>XX</td><td>XX</td><td>X</td><td>XX</td><td>X</td></tr><tr><td>CG5549</td><td>-0.1158</td><td>0.8007</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>slif</td><td>-0.1537</td><td>0.9087</td><td></td><td>XX</td><td>XXX</td><td>XXXX</td><td></td><td>X</td></tr><tr><td>CG7888</td><td>-0.1735</td><td>0.929</td><td>XX</td><td></td><td>X</td><td>XXX</td><td></td><td></td></tr><tr><td>Vmat</td><td>-0.1715</td><td>0.9422</td><td>XXX</td><td></td><td>XX</td><td></td><td></td><td></td></tr><tr><td>CG8850</td><td>-0.227</td><td>0.9604</td><td></td><td>X</td><td></td><td>XXX</td><td>XXX</td><td></td></tr><tr><td>CG15088</td><td>-0.2275</td><td>0.9631</td><td>X</td><td></td><td>X</td><td>XX</td><td></td><td></td></tr><tr><td>CG1139</td><td>-0.2303</td><td>0.9736</td><td></td><td>XXX</td><td>X</td><td>XXX</td><td></td><td>X</td></tr><tr><td>VGAT</td><td>-0.2725</td><td>0.9929</td><td>X</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CG1607</td><td>-0.3837</td><td>0.9988</td><td>XXX</td><td>XX</td><td>XX</td><td>XX</td><td>XXX</td><td>XX</td></tr><tr><td>CG13248</td><td>-0.434</td><td>0.9998</td><td>XX</td><td></td><td>X</td><td></td><td></td><td></td></tr><tr><td>CG10804</td><td>-0.435</td><td>0.9999</td><td>XX</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Eaat2</td><td>-0.5347</td><td>1.00</td><td>XXX</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Jhl-21 was initially discovered in a screen for juvenile hormone (JH) inducible genes¹⁸,¹⁹. It was classified as a SLC7A5-11 family member based on sequence analysis⁴, and amino acid uptake into cultured cells²⁰. The *Drosophila* genome encodes four more putative SLC7A5-11 paralogs: Genderblind (GB), Minidiscs (Mnd), CG9413 and CG1607⁴,²⁰. GB controls extracellular glutamate levels, which in turn regulates the number of iGluRs in the glutamatergic NMJ⁴,²¹. This unexpected function -regulation of iGluRs in synapses- was recently shown to be also conserved in mice²².

We therefore turned our attention to testing explicitly whether Jhl-21 and iGluRs in the *Drosophila* are functionally related, as predicted by ERC.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Candidate amino acid transporters involved in glutamatergic transmission. The *Drosophila* genome encodes 39 genes showing homology to amino acid transporters⁴,⁴⁴. 10 of those were predicted to be co-evolving with GluRs by usage of ERC (blue) (ERC-values >0.23 and *p* < 0.05). Out of those 39 putative transporter genes 11 show the highest level of expression in the same tissue as either mGluRA (highest expression in the CNS) or iGluRs (highest expression in the body wall) (green). We picked one candidate gene, JhI-21, for further analysis.</figcaption></figure>

**JhI-21 is expressed in motor neurons at the glutamatergic NMJ of 3rd instar larvae.** We first investigated whether JhI-21 is expressed at the larval neuromuscular junction close to the glutamate receptor subunits with which ERC revealed it to coevolve. The *Drosophila* larval NMJ is a tripartite synapse containing the presynaptic motor neuron, the postsynaptic muscle cell, and adjacent glia. All three cell-types are known to be involved in regulating synaptic development and physiology²³. First, we designed a polyclonal anti-JhI-21 antibody and confirmed its specificity in JhI-21 null mutant embryos (Figure S1). As shown in Fig. 3a, we could detect the anti-JhI-21 labeling at the NMJ. This anti-JhI-21 signal co-localized with the anti-HRP, which marks presynaptic motor terminals. These data therefore strongly suggest that JhI-21 is expressed in motor neuron terminals.

However, the close proximity of different cell types (neurons, muscle, or glia) at the NMJ makes it difficult to conclude cell-type expression based on light microscopy alone. To address this problem, we generated transgenic flies in which Gal4 is expressed under the control of a putative *JhI-21* regulatory region, and then examined the expression pattern of UAS-nSyb::GFP using the Gal4/UAS system²⁴ (Fig. 3b). Synaptobrevin (Syb) is a synaptic terminal protein, and therefore nSyb::GFP will be enriched at synaptic endings when expressed in neurons²⁵. Consistent with the conclusion that JhI-21 is expressed in motor neurons, nSyb::GFP expressed under control of *JhI-21*-Gal4 was localized to synaptic terminals. To ensure that nSyb::GFP did not localize similarly when expressed in glia or muscles, we also expressed nSyb::GFP under the control of well-characterized glial and muscle cell Gal4 drivers (repo-Gal4 and 24B-Gal4, respectively; Fig. 3c,d). As expected, the nSyb::GFP fluorescence pattern in these cases was drastically different from when nSyb::GFP is expressed in motor neurons. We therefore conclude that JhI-21 is normally expressed in motor neurons and (based on antibody staining) at least in part localized to motor neuron terminals.

To determine whether JhI-21 was expressed in other parts of the motor neurons, we examined the expression of JhI-21 protein in the ventral nerve cord (VNC), and compared it to the location of glutamatergic motor neurons cell bodies marked using the OK6-Gal4, which is expressed in motor neurons, or OK371-Gal4, which is expressed in glutamatergic neurons, using membrane-bound GFP (mCD8::GFP) as a reporter transgene. As shown in Fig. 3e, anti-JhI-21 labeling was detected in cell bodies at the VNC. While most neurons appear to express low levels of JhI-21, some neurons identified as a subset of glutamatergic motor neurons by the OK6-Gal4 or OK371-Gal4 driven expression of mCD8::GFP, express JhI-21 at a higher level (solid arrow in Fig. 3e; white bordered arrow in Fig. 3f). According to the results obtained with our anti-JhI-21 antibody and our *JhI-21*-Gal4 transgene, JhI-21 seems to be expressed in neurons of 3rd-instar larvae, with highest expression in motor neurons (Fig. 3, Figure S2). However, our *JhI-21*-Gal4 transgene recapitulates just partially the endogenous expression of JhI-21 expression in the brain (Figure S2), suggesting that JhI-21 expression in motor neurons may be variable or actively regulated.

**Loss of JhI-21 expression increases iGluR clustering in NMJs.** The expression pattern of JhI-21, along with its co-evolution with glutamate receptors, raised the possibility that JhI-21 regulates NMJ physiology. To test this hypothesis, we performed two-electrode voltage-clamp electrophysiology to measure spontaneous excitatory junction currents (sEJCs) in control and *JhI-21* mutant larval NMJs (muscle 6, segment A3). We used the morphic *JhI-21* allele P{SUPor-P}JhI-21[KG00977], (hereafter referred to as ‘*JhI-21* KG’), which was generated by the BDGP Gene Disruption Project and carries a P{SUPor-P} transposable element in the first exon of the *JhI-21* gene²⁶. The *JhI-21* KG allele was used homozygous or in combination with a deficiency (*Df1* is *Df(2L)esc-P3-0*), in which the *JhI-21* gene is deleted²⁷. Another strain, P*JhI-21* [EP1187] (*JhI-21* EP), contains UAS-binding sites and was used to over-express JhI-21²⁸. As controls, we used *w*¹¹¹⁸ (control 1) and *P*\*82 (control 2). *P*\*82 is a clean excision allele of the *JhI-21* P-element [KG00977] (Fig. 4a)²⁹. As shown in Fig. 4d, sEJC amplitudes are strongly correlated with *JhI-21* expression as measured by q-PCR. sEJC amplitude distributions measured in the two control strains are nearly indistinguishable from each other (P > 0.05, n.s., Kolmogorov-Smirnov test). The *JhI-21* morphs, however, display significantly larger spontaneous postsynaptic currents (P < 0.05). In contrast, larvae overexpressing *JhI-21* show smaller postsynaptic currents (P < 0.01 compared to control 1; Fig. 4b,c; Figure S3).

sEJC frequencies were not significantly different between genotypes (control 1 = 2.6 ± 0.5 Hz; control 2 = 2.3 ± 0.3 Hz; *JhI-21* KG = 1.9 ± 0.2 Hz; *JhI-21* KG/*Df1* = 2.5 ± 0.3 Hz; tub-Gal4/*JhI-21* EP = 1.7 ± 0.3 Hz; n.s.) ruling out the possibility of changes in presynaptic release.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** *Jhl-21 is expressed in presynaptic motor nerve terminals.* (a) Confocal image of third-instar larval NMJ formed on ventral longitudinal muscles 6 and 7 marked by anti-HRP (magenta) co-stained with antibodies against Jhl-21 (green). (b–d), Representative confocal projections of third-instar larval NMJs formed on ventral longitudinal muscles 6 and 7, stained with antibodies against HRP (magenta) and against GFP (green). Only *Jhl-21*-Gal4 driven transgenic synapse-tethered GFP (nSyb::GFP) expression co-localizes with anti-HRP labeling (B). Glial (*repo*)-Gal4 or muscle (24-B)-Gal4 driven nSyb::GFP does not show co-localization with anti-HRP labeled motor-neurons. (e,f) Ventral nerve cord (VNC) of third instar larva. (E) Cell bodies of motor-neurons marked by OK6-Gal4 driving the expression of membrane bound mCD8GFP (green), and are co-labeled by anti-Jhl-21 (magenta). (F) Cell bodies of glutamatergic neurons marked by OK371-Gal4 driving the expression of membrane bound mCD8GFP (green), and are co-labeled by anti-Jhl-21 (magenta). Two arrowheads are pointing two examples of cells expressing both the GFP marker and Jhl-21. Scale Bar NMJ = 20 mm. Scale Bar VNC = 40 µm.</figcaption></figure>

As expected based on changes in sEJC size, overexpression of Jhl-21 in the *Jhl-21 EP* strain led to a reduction in evoked excitatory junction currents (eEJCs; Figure S4). This suggests that increased Jhl-21 leads to less muscle excitation. However, we observed no significant differences between the controls and loss-of-function mutant genotypes (Figure S4).

The number of NMJ branches, number, and bouton area were not statistically different in between controls and *Jhl-21* mutant alleles (Fig. 4e).

The dramatic changes in sEJC amplitude without apparent changes in NMJ morphology or frequency of neurotransmitter release suggest that *Jhl-21* mutant NMJs might cause alterations in the number of
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 4.** The electrophysiological response at the muscle 6/7 NMJ is controlled by Jhl-21. (a) Schematic representation of the *Jhl-21* genomic locus. Exons are indicated by boxes, translated exons of the *Jhl-21* gene by black boxes, 5' and 3' untranslated regions of the *Jhl-21* gene as gray boxes. Triangles represent the inserting region of P-elements used to generate hypo-(stripped triangle) or hypermorph (grey triangle) *Jhl-21* alleles. (b) Representative traces for two-electrode voltage clamp experiments from the larval muscle 6 (LIII). As wildtype we used *w*¹¹¹⁸ (control 1) and a clean excision of P(KG00977) (control 2). (c) Relative cumulative frequency histogram of sEJC ("mini") amplitudes from different genotypes in third-instar *Drosophila* larvae. Rightward shift (*Jhl-21* KG and *Jhl-21* KG/Df1, blue) indicates increase in the abundance of current-conducting postsynaptic receptors, i.e larger synaptic currents. Note the shift to the left (i.e. decreased receptor number) in Jhl-21 overexpression mutants (gray) N = 6-10 animals, 800-3.400 events measured. Kolmogorov-Smirnov test was used to compare the cumulative distributions two by two (*P* < 0.05; **P** < 0.01). (d) Negative correlation between the *Jhl-21* mRNA levels and the number of sEJC amplitudes at the NMJ. (e) Motor neurons were stained using anti-HRP and confocal images were taken at the 6/7 NMJ. Neither *Jhl-21* KG nor *Jhl-21* KG/Df1 showed alterations of synaptic morphology in terms of NMJ branches, bouton number, or bouton area using a Kruskal-Wallis test. Error bars represent SEM in d and e.</figcaption></figure>

postsynaptic glutamate receptors. To test this, we measured postsynaptic glutamate receptor protein abundance immunohistochemically.

*Drosophila* muscle cells express two different subtypes of heterotetrameric iGluRs, called 'A-type' and 'B-type', which can be distinguished immunocytochemically using antibodies against the subunits unique to each receptor type: GluRIIA, and GluRIIB³⁰⁻³³.
<figure><img src="image_12.png" /><figcaption>**Figure 5.** Postsynaptic glutamate receptor immunoreactivity is inversely proportional to *Jhl-21* expression levels. Genotypes used are the same than the ones shown in Fig. 2. (a) *Left*, representative confocal images showing the accumulation of GluRIIA receptor subunits at the NMJ in *Jhl-21* hypomorphs (*Jhl-21* KG/Df1) compared to a control genotype control 2; *middle*, quantification of the NMJ GluRIIA and GluRIIB abundance in various genotypes; *right*, negative correlation between the *Jhl-21* mRNA levels and the abundance of type-A glutamate receptor subunits at the NMJ. (b) Immunohistochemistry of type B glutamate receptor subunits (GluRIIB) shows similar (negative) correlation with *Jhl-21* transcript levels (N = 4–13 animals). Error bars represent SEM; statistical test: Kruskal-Wallis test. (*p* < 0.05; **p** < 0.01; ***p** < 0.001).</figcaption></figure>

Using anti-GluRIIA and anti-GluRIIB antibodies, we compared the immunoreactivity for both A- and B-type receptors in wild-type controls and in the *Jhl-21* alleles. The strongest *Jhl-21* hypomorph allele (*Jhl-21* KG/Df1) causes a ~5-fold increase in the postsynaptic receptor abundance for A-type receptors and a ~3.5-fold increase in B-type receptors compared to control. In contrast, overexpression of *Jhl-21* (*tub-Gal4/Jhl-21* EP) caused a significant decrease in A- and B-type receptors compared to control genotypes at the NMJ (Fig. 5). We obtained the same effect when looking at the GluRIIC subunit, which is shared by A- and B-type receptors (Fig. 6). Overall, as observed for sEJC amplitudes, the iGluR immunoreactivity was negatively correlated to the amount of larval *Jhl-21* mRNA detected by q-PCR (Fig. 5). Taken together, these results show that *Jhl-21* negatively regulates glutamatergic transmission and the abundance of iGluR protein at the larval NMJ.

We next wondered if the action of *Jhl-21* on iGluR clustering could be through the action of glutamate, based on previous work showing that the related transporter GB controls iGluR abundance via regulation of extra-cellular glutamate<sup>4,5</sup>. Consistent with this idea, we could fully rescue the phenotype of our strongest mutant by bathing *Jhl-21* KG/Df1 NMJs with 2 mM glutamate during 24 h when measuring GluRIIC staining (Fig. 6). Unfortunately our intense efforts to test if *Jhl-21* could transport (or not) glutamate, by glutamate quantification in the hemolymph or by using the S2 cell model system, failed to give conclusive results (Figure S7).

**Jhl-21 regulates locomotor behavior.** Our results reveal *Jhl-21* as an unexpected regulator of NMJ iGluR abundance and spontaneous synaptic transmission strength. But what role does this novel form of regulation play? To test whether changes in *Jhl-21* expression at the NMJ and/or in the CNS could have behavioral consequences we measured speed and meandering (turning rate) at two different physiological stages of 3rd-instar larvae: feeding and wandering (shortly before pupation).

In the absence of food, wildtype feeding-stage larvae moved significantly faster than wildtype wandering-stage larvae (0.093 ± 0.003 cm/sec, and 0.070 ± 0.003 cm/sec respectively; *p* < 0.0001). Also, feeding-stage larvae had
<figure><img src="image_13.png" /><figcaption>**Figure 6.** Postsynaptic GluRIIC immunoreactivity is dependent on *JhI-21* expression and ambient glutamate. Bar graphs showing postsynaptic GluRIIC abundance in a control genotype (precise excision of *JhI-21* P[KG00977]) and a *JhI-21* hypomorph allele *JhI-21* KG/Df1, and showing postsynaptic reduced GluRIIC abundance in *JhI-21* KG/Df1 mutants incubated with 2 mM ambient glutamate or without for 0h or 24 hours. (N = 3-6 animals per genotype and condition). \*p < 0.05 (control compared to the other incubation times 2 by 2); \*\*p < 0.01; ns: not significant (Mann-Whitney tests). Error bars represent SEM.</figcaption></figure>

a lower turning rate compared to the wandering-stage larvae (1651 ± 151 deg/cm for feeding-stage larvae compared to 2918 ± 335 deg/cm for wandering ones; *p* < 0.01; Fig. 7a,b).

To test whether JhI-21 regulates locomotor behavior, we used the hypomorphic *JhI-21* KG allele. Wandering homozygous *JhI-21* KG hypomorphs still showed a significant decrease in speed (from 0.078 ± 0.003 cm/sec in feeding stage to 0.066 ± 0.003 cm/sec in wandering stage; *p* < 0.05), but no change in meandering (2263 ± 210 deg/cm in feeding stage versus 2677 ± 269 deg/cm in wandering stage; n.s.). In the strongest viable mutants, *JhI-21* KG/Df1 hypomorphs, neither speed (0.083 ± 0.003 cm/sec in feeding stage and 0.081 ± 0.003 cm/sec in wandering stage) nor meandering (2685 ± 331 deg/cm in feeding stage and 2126 ± 195 deg/cm in wandering stage) differed between feeding and wandering larvae (Fig. 7a,b; Figure S5). Therefore, *JhI-21* mutants do not exhibit normal differences between feeding and wandering larvae in locomotor characteristics (speed and meandering) that are displayed by control animals.

To confirm the association between *JhI-21* expression level and locomotor behavior, we divided within each genotype the mean value scored for feeding animals by the mean value scored for wandering animals, and plotted this ratio against the average *JhI-21* mRNA levels. The ratio for speed negatively correlates with *JhI-21* mRNA levels, while the ratio for meandering positively correlates with the expression levels of *JhI-21* mRNA (Fig. 7a,b).

We also looked at other locomotor phenotypes such as the number of stops and go, and the peristaltic waves along the length of the larval body axis. No difference between feeding and wandering larvae were found (Figure S6). Since all genotypes tested in this assay showed also no difference in developmental time or lethality (Fig. 7c), *JhI-21* expression at the NMJ and/or in the CNS seems to specifically regulate locomotor behavior shifts (speed and meandering) that normally occur in late 3rd-instar larvae.

**Synapse physiology differs in feeding and wandering 3<sup>rd</sup>-instar larvae.** If JhI-21 regulates NMJ physiology, and this regulation affects feeding and wandering behavior, then there should be differences in NMJ physiology between feeding and wandering larvae. To test this, we compared sEJC amplitude distributions in feeding and wandering animals. In our control genotype we observed overall smaller amplitudes of sEJCs when larvae were in the wandering stage compared to feeding stage larvae (*p*<0.0001). This change could be explained by JhI-21-mediated inhibition of iGluR clustering at the postsynaptic muscle cell during the wandering stage (Fig. 4). To test this hypothesis explicitly, we compared sEJCs in our strongest JhI-21 hypomorph (*JhI-21* KG/Df1). The JhI-21 hypomorphs exhibit overall larger spontaneous miniature postsynaptic currents in wandering stage compared to feeding 3<sup>rd</sup>-instar larvae, consistent with the hypothesis (*p*<0.001; Fig. 8).

**JhI-21 is differentially expressed in feeding and wandering 3rd-instar larvae.** How does JhI-21 regulate feeding and wandering behavior? One possibility is that JhI-21 expression at the NMJ, and thus the strength of NMJ regulation, differs between feeding and wandering stages. To test this we stained the nervous system of feeding and wandering 3rd instar larvae with the anti-JhI-21 antibody and analyzed the expression levels. In the VNC, where the cell bodies of motor neurons are located, expression of JhI-21 is significantly higher in feeding animals than in wandering ones (feeding stage 1.275 ± 0.266, wandering stage 0.531 ± 0.117; *p* < 0.05) (Fig. 9a-c). The anti-Bruchpilot antibody nc82 labels the neuropil³⁴ and was used as an internal control for this experiment. As shown in Fig. 9d, the signal obtained with nc82 was constant between feeding-stage (94 ± 6) and wandering-stage larvae (107 ± 10). In contrast JhI-21 staining could never be detected at the NMJ of feeding-stage animals, but wandering stage animals showed strong anti-JhI-21 labeling in the NMJ (Fig. 9e,f). Thus, JhI-21 subcellular localization appears to shift between feeding and wandering stages. During feeding stage,
<figure><img src="image_14.png" /><figcaption>**Figure 7.** *Jhl-21 mutants lack age-dependent shift in locomotor behaviors. (a,b) Parameters of locomotion were analyzed in wildtype larvae (control 1, *w*¹¹¹⁸) and Jhl-21 mutants (*Jhl-21* KG and *Jhl-21* KG/*Df*1). (a) right, Speed is significantly decreased in wandering wildtype larvae compared to feeding individuals. *****p* < 0.0001. This difference is still present, although reduced in *Jhl-21* KG mutants. *p* < 0.05. Difference in speed between feeding and wandering larvae is not present in *Jhl-21* KG/*Df*1 mutants. Left, the ratio of average speed of feeding larvae/average speed of wandering larvae correlates with *Jhl-21* mRNA levels. (b) right, the meander describes amount of turnings in degree per cm and was used to characterize the turning behavior. 3rd instar wildtype larvae increase their average turning in wandering stage. **p* < 0.01. Neither *Jhl-21* KG nor *Jhl-21* KG/*Df*1 mutants show a significant change in turning behavior between both 3rd-instar larval stages. Left, the ratio of average meander of feeding larvae/average speed of wandering larvae shows negative correlation with *Jhl-21* mRNA levels. (c) Jhl-21 mutants show no difference in lethality or developmental time as compared to the control genotype. Statistical test in (a,b): 2-Way-ANOVA, N = 30, Error bars represent SEM. Statistical test in c Mantel-Cox test: N = 131-334 per genotype.</figcaption></figure>

Jhl-21 is predominantly in the cell bodies within the VNC and not in the motor terminals. During wandering stage, Jhl-21 is strongly localized in the motor terminals at the NMJ, and less abundant in the cell bodies within the VNC.

## Discussion

**Evolutionary rate co-variation identified Jhl-21 as a member of glutamatergic signaling.** In this study we used the “covariation of protein evolutionary rates” to screen for proteins that might play unsuspected roles in glutamatergic synapse physiology. It has been previously demonstrated that ERC signatures provide a powerful method to reveal functionally related proteins or proteins acting in the same complex⁸,³⁵,³⁶. As expected, physically interacting ionotropic glutamate receptor (iGluR) subunits GluRIIA, GluRIIB, GluRIIC,
Frequency distribution histogram

<figure><img src="image_15.png" /><figcaption>**Figure 8.** The age-dependent shift in electrophysiological response at the muscle 6/7 NMJ is controlled by Jhl-21. Cumulative frequency histograms of sEJC ("mini") amplitudes from third-instar feeding and wandering *Drosophila* larvae. Leftward shift in wandering control animals indicates an age-dependent decrease in the abundance of current-conducting postsynaptic glutamate receptors, i.e smaller synaptic currents. This shift is inverted to rightward in the *Jhl-21* KG/Df1 allele. N = 4–5 animals (representing 10–17 NMJs), 2400–3600 events measured per phenotype; Statistical test: Kolmogorov-Smirnov test (***P*<*0.001*, ****P*<*0.0001*).</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 9.** Jhl-21 is differentially expressed in the nervous system of feeding and wandering 3rd-instar control larvae. (a,b) Confocal projection of *w*<sup>1118</sup> 3rd-instar larval VNC stained with anti-Jhl-21 (magenta) and NC82 antibody (green). (c,d) Quantification of signal intensity reveals higher expression of Jhl-21 in the ventral nerve cord of feeding larvae compared to wandering ones (\*\*p*<*0.01*, N = 12, Mann Whitney test). Intensity of nc82 does not change between these larval stages (N = 12, unpaired t-test). (e,f) Confocal images of third-instar larval (LIII) NMJs of muscles 6 and 7, stained with antibodies against HRP (stains all neuronal membrane, green) and Jhl-21 (magenta). At the NMJ anti-Jhl-21 labeling is detectable at wandering stage, not at feeding stage. Scale Bar (brains) =40 mm. Scale Bar (NMJs) =20 mm. Error bars represent SEM in c and d.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 10.** JhI-21 co-evolved with glutamate receptors, and regulates strength of glutamatergic signaling at the larval NMJ. JhI-21 is expressed in motor neurons together with mGluRA, and regulates synaptic strength by inhibiting the clustering of iGluRs at the postsynaptic muscle cell specifically during the wandering stage in third-instar larvae. The schematic representation illustrates two possible functions of JhI-21 on iGluR clustering in wandering larvae. On the left is shown the direct export of glutamate close to the synapse, which inhibits iGluR clustering. On the right, an indirect action of JhI-21 is proposed through the transport of leucine. Leucine acts on the activity of the Glutamate dehydrogenase (Gdh), which can reciprocally catalyze the production of a-ketoglutarate (a-KG) to glutamate within the motoneuron. Glutamate then could be released to the hemolymph surrounding the NMJ via an unidentified transporter, and act on iGluR clustering.</figcaption></figure>

GluRIID, and GluRIIE showed overall positive ERC values. Even if the metabotropic glutamate receptor mGluRA is neither physically interacting nor co-expressed in the same cell than iGluR subunits, it also showed positive ERC values when compared to iGluR subunits, likely due to its action in the same neurophysiological pathway. We therefore considered whether the function of two proteins in the same synapse can be sufficient to gain positive ERC values even though both the proteins are not expressed in the same cell.

We next used ERC to determine whether any of 39 putative amino acid or biogenic amine transporters genes were co-evolving with GluRs, and got 10 significant hits. Based on ERC scores and tissue expression, we selected the amino acid antiporter JhI-21 for further analysis. Specifically, we sought to determine whether JhI-21 did indeed function in the glutamatergic signaling pathway as suggested by ERC.

Using two independent strategies we showed that JhI-21 is localized in glutamatergic motor neurons in 3^rd^-instar larvae. JhI-21 showed the strongest correlation in terms of co-evolution with mGluRA, which is also expressed in motor neurons^15^. Although co-expressed proteins tend to have higher ERC values in general^7^, we want to highlight that JhI-21 also showed positive ERC score with the post-synaptically expressed iGluRs, indicating a previously unsuspected but important role for JhI-21 in adjusting the strength of glutamatergic neuromuscular transmission.

Interestingly, other significant hits are not even expressed in the CNS or the body wall. For example CG8785 is expressed only in the digestive system and malpighian tubules. Both tissues also contain iGluRs. What might be the role of iGluRs in the digestive and *Drosophila* renal system and what could be the link to the co-evolution with CG8785 is an interesting subject for further studies.

**Possible direct action of JhI-21 in adjusting synaptic strength at the larval NMJ.** The fact that JhI-21 is expressed at the glutamatergic synapse allowed us to further elucidate the functioning of the NMJ in *Drosophila*. Specifically, we found that JhI-21 expression in motor neurons leads to inhibition of synaptic transmission by reducing the clustering abundance of iGluRs in neuromuscular junctions of late stage 3^rd^-instar larvae. This could be linked to the extracellular glutamate concentration since additional application of glutamate can compensate a lack of JhI-21 activity in *JhI-21* hypomorph mutants when measuring the amount of iGluR expression at the NMJ (Fig. 6). It is therefore reasonable to ask whether JhI-21 has a direct impact on extracellular glutamate levels at the NMJ.

Indeed JhI-21 was classified as a SLC7A5-11 family member based on sequence analysis⁴. These transporters can be found in many species across the animal kingdom and other members of this transporter family have been proven to regulate glutamatergic signaling directly by adjusting extracellular glutamate levels in mammals and *Drosophila*⁴,²⁰. This extracellular glutamate is mostly independent of synaptic vesicular release⁶,³⁷,³⁸, and is partly attributable to glial expressed SLC7A11 transporters⁴,²¹,³⁹. Those transporters in mammals act as hererodimers consisting of a 'light chain' SLC7A5-11 core catalytic subunit in combination with a 'heavy chain' subunit (such as 4F2hc protein) that is required for trafficking to the plasma membrane⁴⁰⁻⁴³. GB, which is the homologue of JhI-21, previously revealed to be expressed in a particular subset of glia at the larval NMJ, was the first *Drosophila* SLC7A5-11 member shown to directly impact the physiology of this glutamatergic synapse by exporting glutamate into the hemolymph⁴,⁴⁴. This led us speculate that JhI-21 might also export glutamate directly. Nevertheless, when compared to GB, JhI-21 does not affect the overall hemolymph glutamate concentration²⁹. Presumably, JhI-21 modulates postsynaptic iGluR levels via some other mechanism, or regulates glutamate only very locally near the synapse, in contrast to GB, which controls glutamate levels more globally (Fig. 10). Using JhI-21 expressing
S2 cell we could not demonstrate that this transporter might directly or not transport glutamate (Figure S7), although we found that incubating larval NMJs in glutamate was able to rescue the JhI-21 phenotype (Fig. 6).

We also considered whether JhI-21 might be involved in loading glutamate into neurotransmitter vesicles and might thereby have a direct impact in glutamatergic transmission. No SLC7A5-11 family member is known to be a component of synaptic vesicle membrane. Also it has been shown that at the larval NMJ the vesicular glutamate transporter (VGlut, belonging to the SCL17 family) is necessary and sufficient to fill vesicles with glutamate. Absence of VGlut leads to empty vesicles at the NMJ⁴⁵. This indicates that there is no other transporter involved in the filing of glutamate into these neurotransmitter vesicles.

**Possible indirect action of JhI-21 in adjusting synaptic strength at the larval NMJ.** It has been previously shown that non-vesicular glutamate release is partially dependent on the amount of glutamate present in the motor neuron⁶. This intracellular glutamate pool is regulated via glutamate metabolizing enzymes. Glutamate oxaloacetate transaminase (GOT) produces glutamate from aspartate. Glutamate decarboxylase (GAD) catalyzes the decarboxylation of glutamate to γ-aminobutyric acid (GABA), and glutamate synthase (GS) converts glutamate to glutamine⁴,⁴⁴. In addition glutamate dehydrogenase (GDH) catalyzes the reversible formation of glutamate to α-ketoglutarate⁴⁶. It has been shown that misexpression of glutamate metabolizing enzymes alters postsynaptic iGluR clustering by regulating presynaptic intracellular glutamate concentrations⁶. Instead of directly transporting glutamate, JhI-21 could therefore function to transport metabolites, activators, or suppressors associated with regulation of glutamate levels.

For example, it has been shown that JhI-21 in S2 cells transports leucine, which acts as an allosteric activator of the glutamate dehydrogenase (Gdh) and could thereby modify intracellular glutamate levels (Fig. 10)²⁰,⁴⁷.

**Specific role of JhI-21 during late larval development.** JhI-21 expression was also previously shown to be dependent on Juvenile Hormone (JH)¹⁸,¹⁹. The role of this hormone in many insects is to maintain juvenile morphological characteristics⁴⁸. In *Drosophila*, several studies suggested a role for JH in adult behaviors like foraging and sexual maturation¹⁹,⁴⁹. The role of JH in larvae, however, remains unclear. In *Drosophila* larvae, the titer of JH III, the active form of JH in Diptera, is high during the first and second larval stages. In feeding stage 3rd-instar larvae, JH III levels decrease drastically before increasing again during wandering prior to pupation⁵⁰,⁵¹. These differences in JH III levels might be responsible for the differences of JhI-21 expression and subcellular localization that we observed in this study. This shift in JH III-dependent expression of JhI-21, and the resulting shift in postsynaptic sensitivity occurs at the same time as the behavioral switch in late 3rd-instar larvae. When placed on agar plate, which is low in nutrients, feeding larvae moved at a higher speed compared to wandering individuals. This hypermobility in feeding larvae could reflect aggressive searching for food. Wandering larvae, on the other hand, move at a lower pace but turn twice as much compared to feeding larvae (Fig. 7). This behavior could reflect a search for the best place to start pupation. Both behaviors are altered in *JhI-21* hypermorphs, where feeding and wandering animals move at nearly the same speed and display similar turning rates. Although we could demonstrate that JhI-21 is involved in the age dependent shift of synaptic strength between feeding and wandering larvae, we can not rule out that other JhI-21 expressing cells in the brain are presumably involved in those behavioral changes. In any case, it is nonetheless clear that JhI-21 is important for regulating the change of behavior during late 3rd-instar life.

In summary, we show that ERC values can be used to screen for proteins that work in the same neurotransmitter pathway especially when combined with a second screening method such as comparison of expression patterns. Using this method, we identified JhI-21 as a novel regulator of synaptic glutamate signaling. We found that JhI-21 is expressed in motor neurons where it regulates the developmental specificity of synaptic strength by inhibiting the clustering of post-synaptic iGluRs (Fig. 10). We were also able to highlight the role of JhI-21 on late larval locomotor activity, and show that changes in the subcellular distribution of JhI-21 within motor neurons might explain the differential regulation of JhI-21 on NMJ strength and behavior.

Glutamate is the predominant excitatory neurotransmitter in the central nervous system in mammals. Thus glutamate is involved in the control of a wide range of brain functions. It is also a key player in many neurological diseases⁵²,⁵³. In this context glutamate transporters play a key role in regulating extracellular glutamate levels to maintain dynamic synaptic signaling processes. Therefore our identification of JhI-21 as an important actor on glutamatergic and locomotor physiology in *Drosophila* is suggesting that its possible ortholog (LAT-1³⁶) might also have an impact on such pathways in mammals. Beside the known glutamate transporters (e.g. vGluTs loading glutamate into synaptic vesicles, and EAATs removing the excess of glutamate by surrounding glial cells)⁵⁴, other transporters such as LAT-1, which are possibly not transporting glutamate could also have a major impact on glutamate receptor physiology in the mammalian nervous system. This would be a major clue to explore new strategies to cure neurological diseases.

Our data are also highlighting the potential of a new emerging technique in finding genes that co-evolve: ERC. This technology proved to be extremely powerful in the recent past to find molecular partners that interact in the same cells³⁶. Here we show that it has also the power to decipher specific pathways such as the glutamatergic physiology and the age specific control of locomotor activity even if the products of these genes are expressed in different cell types (JhI-21 and glutamate receptors). Therefore the ERC technique would be particularly useful to reveal some pluricellular molecular networks such as the interaction between glial cells and neurons, or the plasticity of interacting neurons during development or during a learning task by identifying genes that are co-evolving. Such applications are facilitated by the ERC analysis webserver that will perform custom analysis genome-wide for user-chosen genes⁵⁵.
## Materials and Methods

**Calculation of evolutionary rate covariation statistics.** ERC values were calculated as previously described³⁶. Briefly, orthologous genes sequences from 12 *Drosophila* genomes were aligned and used to estimate gene-specific branch lengths across the species phylogeny. These branch lengths were transformed into relative rates of evolution based on the average genome-wide amount of divergence for each branch¹⁰. Rate covariation (ERC) values were calculated from the relative rates as a correlation coefficient for each gene pair. Statistical tests on groups of genes, such as the glutamate receptors, were performed by comparing the mean ERC value to 10,000 random sets of the same size using a permutation test.

**Drosophila stocks and genetics.** Flies were grown and maintained on regular corn medium, at 25 °C, in a 12 h/12 h light/dark cycle. Control *Drosophila melanogaster* used in this study were *Oregon-R*, P\*82 (precise excision of P{SUPor-P}*JhI-21*KG00977*; electrophysiology and iGluR quantification), and *w*¹¹¹⁸ (electrophysiology, iGluR quantification, and behavior). *JhI-21* mutants P{SUPor-P}*JhI-21*KG00977 (BL12970), P*JhI-21*EP1187 (not available any longer at Bloomington, but still available in our laboratory), and Df(2L)esc-P3-0 (BL3131) were previously characterized²⁹ and obtained from the Bloomington Stock Center. *JhI-21* alleles were re-balanced over *CyO-GFP*. P{TubP-Gal4}LL7/TM3, Sb was obtained from Bloomington (BL5138) and re-balanced over TM3GFP, Ser. UAS-*nSyb-GFP*⁵⁶ provided by B. Hovemann, (Ruhr-University Bochum), *repo-Gal4* (BL7415, obtained from Bloomington), and 24-B-Gal4 (BL1767, obtained from Bloomington), OK6-Gal4 and OK-371-Gal4 were provided by H. Abele (Heinrich-Heine University Düsseldorf). All genotypes used for behavioral experiments were backcrossed 5 times to an isogenized *w*¹¹¹⁸ fly stock.

**JhI-21-Gal4.** UAS-Sequences were removed from pUAST attB⁵⁷ using EcoRI and HindIII. Vector was treated with Klenow fragment and religated to make attB. Heat-shock (HS) minimal promotor-, Gal4-, and hsp7αolyA-sequence were removed from pCHS-Gal4⁵⁸ using NotI and cloned into attB to make Gal4 attB. JhI-21 promoter fragment was PCR amplified from genomic DNA of Canton-S flies by using primers, which introduced KpnI sites for cloning (underlined): *GGTACCGGGATTCTTCTGCTTACCCTCT* and *GGTACCGCACCGATAGGAGGATGTATTC*. The amplicon was subcloned to pGEM®-T Easy (Promega), re-excised and cloned to Gal4 attB using KpnI. The resulting transgenic flies were generated by site-directed integration using attB44 (provided by J. Bischof, University of Zürich). Injections were performed by Genetic Services Inc. (Sudbury, MA, USA).

**qPCR for quantification of JhI-21 expression.** For reverse transcription RT-PCR, total RNA was isolated using Trizol (Invitrogen, Carlsbad, CA) extraction⁵⁹. For the *JhI-21* expression experiments, total RNA was isolated from whole larvae. RNA (500 ng, quantified spectrophotometrically) was reverse transcribed using oligo-dT primers and standard methods. 10% of the cDNA product was used to amplify *JhI-21* and *actin5C* cDNA fragments by PCR. A 122 fragment of *JhI-21* was amplified using the following primer pair: TTGTTTACCACGGCGAAATAG and CTTTGTGACGGAGGAGCTACA; a 200 bp fragment of *actin5C* was simultaneously amplified using CAAGCCTCCATTCCCAAGAAC and CGTGAAATCGTCCGTGACATC primer pair.

Real-time PCR was performed using an MJResearch Opticon2 real-time thermocycler and quantitative fluorescent detection of SYBR green-labeled PCR product. Relative mRNA abundance was calculated using the "ΔΔCT method"⁶⁰.

**Immunocytochemistry and Confocal microscopy.** 2 Rabbit polyclonal anti-JhI-21 antibodies were raised against a synthesized peptide (DGEEKIVLKRRKLLINGVA) by Thermo Scientific/Open Biosystems, using standard methods and used at a dilution of 1:250. Both gave similar results, and we kept the one giving less background staining. The JhI-21 peptide epitope represents amino acids 30-48 of the predicted *Drosophila* JhI-21 protein. Specificity was verified on JhI-21 null mutant embryos (Figure S2).

NMJs were dissected under *Drosophila* saline supplied with 2 mM Glutamate. Glutamate present in the buffer prevents retraction of glia from the NMJ [4]. Imaging was performed on larval ventral longitudinal muscles (VLM) 6 and 7. For GluRIIA or GluRIIB stainings, NMJs were fixed 30 min in Bouin's fixative (Sigma). Immunostaining measurements of postsynaptic glutamate receptor abundance were performed as previously described⁴,⁶,³¹,³³,⁶¹. For JhI-21 localization experiments, NMJs and larval brains were dissected in PBS + 2 mM glutamate (NMJs) or directly in fixative (brains) followed by a 30 min fixation in 4% PFA in PBS. Primary antibodies were incubated overnight at 4 °C. Mouse monoclonal anti-GFP (Sigma Aldrich, Saint Quentin Fallavier, France) was used at 1:200. Mouse monoclonal anti-ElaV (9F8A9) was obtained from the University of Iowa Developmental Studies Hybridoma Bank (Iowa City, US) and used either at 1:200 or at 1:500. A 1-5 h washing step was performed with at least 3 solution changes, before the incubation of secondary antibodies for 2 h at RT (for embryos) or overnight at 4 °C (for the rest). DylightTM594-conjugated goat anti-Horseradish Peroxidase (HRP) antibody was obtained from Jackson ImmunoReserach Laboratories Inc. (West Grove, US) and used at 1:500. TRITC-conjugated anti-HRP antibody was obtained from Jackson ImmunoReserach Laboratories Inc. (West Grove, US) and used at 1:100. Secondary antibodies (Alexa Fluor 488 or 594 Dye) were obtained from Molecular Probes and were diluted at 1:400. Preparations were mounted in Vectashield (H-1000, Vector Labs) before imaging using a Leica TCS SP2 AOBS or an Olympus Fluoview FV500 laser scanning confocal system. Confocal projections were scanned at 1 µm section intervals, and were orientated and cropped with LeicaLight (also used to obtain Z-projections). Measurements of postsynaptic glutamate receptor density in confocal stacks were made by quantifying mean postsynaptic immunofluorescence intensity relative to fluorescence in surrounding muscle tissue: F synapse /F background membrane⁴.
**Electrophysiology.** All electrophysiological recordings were obtained at 19°C from third instar (110-120 hr after egg laying) larval ventral longitudinal muscle 6 (A3-A4) at -60 mV holding potential using two-electrode voltage clamp technique (TEVC), as previously described⁴,⁶,²³. Dissections and electrophysiology were performed under glutamate free *Drosophila* HL-3 saline⁴. Electrodes for TEVC were filled with 3M KCl yielding a resistance of 30–40 MOhm. Spontaneous excitatory junctional currents (sEJCs) or “minis” were detected and analyzed using Clampfit9/10 template-matching method that identifies synaptic events based on shape matching to a data-based ideal template³³. Axon GeneClamp 500B/AxoClamp 900A and Digidata 1550 (Molecular Devices, Sunnyvale, CA) were used for data acquisition.

The bathing experiment was performed with 2 mM glutamate during 24 h on NMJ cultures as previously described⁴.

**Staging of feeding and wandering LIII larvae.** 3rd-instar larvae were picked from the food (to obtain feeding stage larvae) or the wall (to obtain wandering stage larvae) of a stock vial and then transferred to a dish containing a piece of food for stage confirmation. Animals that stayed in the food were considered as “feeding” larvae. Larvae escaping the food were determined to be “wandering” 3rd-instar.

**Behavior.** Staged larvae were washed in ddH₂O before conducting behavioral analysis. Movement analysis was performed on a 9 cm-circular 2%-agar-agar medium. This agar medium was used upside down to have a perfect flat surface and was placed in the middle of a bigger Pétri dish (13 cm diameter). The space between the edge of the Pétri dish and the agar was filled with ddH₂O to avoid escaping larvae. A single larva was placed on the middle of the agar. Recordings and analysis were performed automatically using EthoVision XT software (Noldus information technology, Wageningen, Netherlands). Dynamic subtraction function was chosen to distinguish the larva from the background. Data collection was started two minutes after the larva was placed on the agar plate. A ten-minute movie was recorded for each individual. Speed and meander behaviors were calculated using the same software, since these parameters reflect most of the locomotor activity characteristics of larvae in our behavioral set up.

**Lethality and Developmental time.** Flies were allowed to lay eggs for 20 h. Approximately 3 × 10⁶ embryos were transferred to fresh vials. 40 h after egg-laying, non-hatched embryos were counted. The amount of pupae was determined each day. After 10 days hatched adults were counted and the amount of lethality calculated.

**Statistics.** All data were transferred to Prism 5.0d (Graphpad) for statistical analysis and tested for normal distribution using the D'Agostino and Pearson omnibus normality test. Two pairs of normally distributed data were analyzed using the paired-t-test. Pairs of data, which did not pass the normality test, were analyzed using the Mann-Whitney test. Three sets of not normally distributed data were analyzed by Kruskal-Wallis test; for quantification of GluRIIA and GluRIIB we compared both controls with each 'test' genotype (i.e. we did 3 comparisons for all 5 genotypes). Normally distributed data with two nominal variables were analyzed using the two-way ANOVA followed by Bonferroni post test. Survival curves were analyzed using the Chi2 test and cumulative frequency histograms by usage of the Kolmogorov-Smirnov test.

**References**

1. Schmid, A. et al. Non-NMDA-type glutamate receptors are essential for maturation but not for initial assembly of synapses at *Drosophila* neuromuscular junctions. *J Neurosci* **26**, 11267–11277 (2006).

2. Schmid, A. et al. Activity-dependent site-specific changes of glutamate receptor composition *in vivo*. *Nature Neuroscience* **11**, 659–666 (2008).

3. Li, H., Peng, X. & Cooper, R. L. Development of Drosophila larval neuromuscular junctions: maintaining synaptic strength. *Neuroscience* **115**, 505–513 (2002).

4. Augustin, H., Grosjean, Y., Chen, K., Sheng, Q. & Featherstone, D. E. Nonvesicular release of glutamate by glial xCT transporters suppresses glutamate receptor clustering *in vivo*. *J Neurosci* **27**, 111–123 (2007).

5. Chen, K., Augustin, H. & Featherstone, D. E. Effect of ambient extracellular glutamate on Drosophila glutamate receptor trafficking and function. *Journal of Comparative Physiology* **195**, 21–29 (2009).

6. Featherstone, D. E., Rushton, E. & Broadie, K. Developmental regulation of glutamate receptor field size by nonvesicular glutamate release. *Nature neuroscience* **5**, 141–146 (2002).

7. Clark, N. L., Alani, E. & Aquadro, C. F. Evolutionary rate covariation reveals shared functionality and coexpression of genes. *Genome research* **22**, 714–720 (2012).

8. Pazos, F. & Valencia, A. Similarity of phylogenetic trees as indicator of protein-protein interaction. *Protein engineering* **14**, 609–614 (2001).

9. Goh, C. S. & Cohen, F. E. Co-evolutionary analysis reveals insights into protein-protein interactions. *Journal of molecular biology* **324**, 177–192 (2002).

10. Sato, T., Yamanishi, Y., Kanehisa, M. & Toh, H. The inference of protein-protein interactions by co-evolutionary analysis is improved by excluding the information about the phylogenetic relationships. *Bioinformatics* (Oxford, England) **21**, 3482–3489 (2005).

11. Lovell, S. C. & Robertson, D. L. An integrated view of molecular coevolution in protein-protein interactions. *Molecular biology and evolution* **27**, 2567–2575 (2010).

12. Clark, N. L., Alani, E. & Aquadro, C. F. Evolutionary rate covariation in meiotic proteins results from fluctuating evolutionary pressure in yeasts and mammals. *Genetics* **193**, 529–538 (2013).

13. Parmentier, M. L., Pin, J. P., Bockaert, J. & Grau, Y. Cloning and functional expression of a Drosophila metabotropic glutamate receptor expressed in the embryonic CNS. *J Neurosci* **16**, 6687–6694 (1996).

14. DiAntonio, A. Glutamate receptors at the Drosophila neuromuscular junction. *International review of neurobiology* **75**, 165–179 (2006).

15. Bogdanik, L. et al. The Drosophila metabotropic glutamate receptor DmGluRA regulates activity-dependent synaptic facilitation and fine synaptic morphology. *J Neurosci* **24**, 9105–9116 (2004).

16. Romero-Calderon, R. et al. A screen for neurotransmitter transporters expressed in the visual system of Drosophila melanogaster identifies three novel genes. *Developmental neurobiology* **67**, 550–569 (2007).

17. Thimgan, M. S., Berg, J. S. & Stuart, A. E. Comparative sequence analysis and tissue localization of members of the SLC6 family of transporters in adult Drosophila melanogaster. *The Journal of experimental biology* **209**, 3383–3404 (2006).
18. Dubrovsky, E. B., Dubrovskaya, V. A., Bilderback, A. L. & Berger, E. M. The isolation of two juvenile hormone-inducible genes in Drosophila melanogaster. *Developmental biology* **224**, 486–495 (2000).

19. Dubrovsky, E. B., Dubrovskaya, V. A. & Berger, E. M. Juvenile hormone signaling during oogenesis in Drosophila melanogaster. *Insect biochemistry and molecular biology* **32**, 1555–1565 (2002).

20. Reynolds, B. *et al*. Drosophila expresses a CD98 transporter with an evolutionarily conserved structure and amino acid-transport properties. *The Biochemical journal* **420**, 363–372 (2009).

21. De Bundel, D. *et al*. Loss of system x(c)- does not induce oxidative stress but decreases extracellular glutamate in hippocampus and influences spatial working memory and limbic seizure susceptibility. *J Neurosci* **31**, 5792–5803 (2011).

22. Williams, L. E. & Featherstone, D. E. Regulation of hippocampal synaptic strength by glial xCT. *J Neurosci* **34**, 16093–16102, doi: 10.1523/JNEUROSCI.1267-14.2014 (2014).

23. Danjo, R., Kawasaki, F. & Ordway, R. W. A tripartite synapse model in Drosophila. *PloS one* **6**, e17131 (2011).

24. Brand, A. H. & Perrimon, N. Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. *Development* (Cambridge, England) **118**, 401–415 (1993).

25. Sudhof, T. C., Baumert, M., Perin, M. S. & Jahn, R. A synaptic vesicle membrane protein is conserved from mammals to Drosophila. *Neuron* **2**, 1475–1481 (1989).

26. Bellen, H. J. *et al*. The BDGP gene disruption project: single transposon insertions associated with 40% of Drosophila genes. *Genetics* **167**, 761–781 (2004).

27. Parks, A. L. *et al*. Systematic generation of high-resolution deletion coverage of the Drosophila melanogaster genome. *Nature genetics* **36**, 288–292 (2004).

28. Rorth, P. A modular misexpression screen in Drosophila detecting tissue-specific phenotypes. *Proceedings of the National Academy of Sciences of the United States of America* **93**, 12418–12422 (1996).

29. Piyankarage, S. C., Augustin, H., Featherstone, D. E. & Shippy, S. A. Hemolymph amino acid variations following behavioral and genetic changes in individual Drosophila larvae. *Amino acids* **38**, 779–788 (2010).

30. Marrus, S. B., Portman, S. L., Allen, M. J., Moffat, K. G. & DiAntonio, A. Differential localization of glutamate receptor subunits at the Drosophila neuromuscular junction. *J Neurosci* **24**, 1406–1415 (2004).

31. Chen, K. & Featherstone, D. E. Discs-large (DLG) is clustered by presynaptic innervation and regulates postsynaptic glutamate receptor subunit composition in Drosophila. *BMC biology* **3**, 1 (2005).

32. Chen, K., Merino, C., Sigrist, S. J. & Featherstone, D. E. The 4.1 protein coracle mediates subunit-selective anchoring of Drosophila glutamate receptors to the postsynaptic actin cytoskeleton. *J Neurosci* **25**, 6667–6675 (2005).

33. Featherstone, D. E. *et al*. An essential Drosophila glutamate receptor subunit that functions in both central neuropil and neuromuscular junction. *J Neurosci* **25**, 3199–3208 (2005).

34. Wagh, D. A. *et al*. Bruchpilot, a protein with homology to ELKS/CAST, is required for structural integrity and function of synaptic active zones in Drosophila. *Neuron* **49**, 833–844 (2006).

35. Clark, N. L. *et al*. Coevolution of interacting fertilization proteins. *PLoS genetics* **5**, e1000570 (2009).

36. Findlay, G. D. *et al*. Evolutionary rate covariation identifies new members of a protein network required for Drosophila melanogaster female post-mating responses. *PLoS genetics* **10**, e1004108 (2014).

37. Miele, M., Boutelle, M. G. & Fillenz, M. The source of physiologically stimulated glutamate efflux from the striatum of conscious rats. *The Journal of physiology* **497** (Pt 3), 745–751 (1996).

38. Bradford, H. F., Young, A. M. & Crowder, J. M. Continuous glutamate leakage from brain cells is balanced by compensatory high-affinity reuptake transport. *Neuroscience letters* **81**, 296–302 (1987).

39. Baker, D. A., Xi, Z. X., Shen, H., Swanson, C. J. & Kalivas, P. W. The origin and neuronal function of *in vivo* nonsynaptic glutamate. *J Neurosci* **22**, 9134–9141 (2002).

40. Sato, H., Tamba, M., Ishii, T. & Bannai, S. Cloning and expression of a plasma membrane cystine/glutamate exchange transporter composed of two distinct proteins. *The Journal of biological chemistry* **274**, 11455–11458 (1999).

41. Chillaron, J., Roca, R., Valencia, A., Zorzano, A. & Palacin, M. Heteromeric amino acid transporters: biochemistry, genetics, and physiology. *American journal of physiology* **281**, F995–1018 (2001).

42. Wagner, C. A., Lang, F. & Broer, S. Function and structure of heterodimeric amino acid transporters. *Am J Physiol Cell Physiol* **281**, C1077–1093 (2001).

43. Verrey, F. *et al*. CATs and HATs: the SLC7 family of amino acid transporters. *Pflugers Arch* **447**, 532–542 (2004).

44. Grosjean, Y., Grillet, M., Augustin, H., Ferveur, J. F. & Featherstone, D. E. A glial amino-acid transporter controls synapse strength and courtship in Drosophila. *Nature neuroscience* **11**, 54–61 (2008).

45. Daniels, R. W. *et al*. A single vesicular glutamate transporter is sufficient to fill a synaptic vesicle. *Neuron* **49**, 11–16 (2006).

46. Papadopoulou, D. & Louis, C. The gene coding for glutamate dehydrogenase in Drosophila melanogaster. *Biochemical genetics* **28**, 337–346 (1990).

47. Meier, C., Ristic, Z., Klauser, S. & Verrey, F. Activation of system L heterodimeric amino acid exchangers by intracellular substrates. *The EMBO journal* **21**, 580–589 (2002).

48. Nijhout, H. F. Physiological Control of Molting in Insects *Amer. Zool.* **21** 631–640 (1981).

49. Meunier, N., Belgacem, Y. H. & Martin, J. R. Regulation of feeding behaviour and locomotor activity by takeout in Drosophila. *The Journal of experimental biology* **210**, 1424–1434 (2007).

50. Bownes, M. & Rembold, H. The titre of juvenile hormone during the pupal and adult stages of the life cycle of Drosophila melanogaster. *European journal of biochemistry/FEBS* **164**, 709–712 (1987).

51. Sliter, T. J., Sedlak, B. J., Baker, F. C. & Schooley, D. A. Juvenile hormone in Drosophila melanogaster-Identification and titer determination during development. *Insect Biochemistry* **17**, 161–165 (1987).

52. Meldrum, B. S. Glutamate as a neurotransmitter in the brain: review of physiology and pathology. *J Nutr* **130**, 1007S–1015S (2000).

53. Sundaram Sundaram, R., Gowtham L. & Nayak, B. S. The role of excitatory neurotransmitter glutamate in brain physiology and pathology. *Asian Journal of Pharmaceutical & Clinical Research* **5**, 1, doi: 82151880 (2012).

54. Vandenberg, R. J. & Ryan, R. M. Mechanisms of glutamate transport. *Physiol Rev* **93**, 1621–1657, doi:10.1152/physrev.00007.2013 (2013).

55. Wolfe, N. W. & Clark, N. L. ERC analysis: web-based inference of gene function via evolutionary rate covariation. *Bioinformatics*, doi: 10.1093/bioinformatics/btv454 (2015).

56. Ito, K. *et al*. The organization of extrinsic neurons and their implications in the functional roles of the mushroom bodies in Drosophila melanogaster Meigen. *Learning & memory* (Cold Spring Harbor, N.Y 5, 52–77 (1998).

57. Bischof, J., Maeda, R. K., Hediger, M., Karch, F. & Basler, K. An optimized transgenesis system for Drosophila using germ-line-specific phiC31 integrases. *Proceedings of the National Academy of Sciences of the United States of America* **104**, 3312–3317 (2007).

58. Apitz, H. pChs-Gal4, a vector for the generation of Drosophila Gal4 lines driven by identified enhancer elements. *Dros. Inf. Serv* **85**, 118–120. (2002).

59. Roberts, D. B. *Drosophila: A Practical Approach*. (Ed2. Oxford: Oxford UP, 1998).

60. Horz, H. P., Barbrook, A., Field, C. B. & Bohannan, B. J. Ammonia-oxidizing bacteria respond to multifactorial global change. *Proceedings of the National Academy of Sciences of the United States of America* **101**, 15136–15141 (2004).

61. Liebl, F. L. & Featherstone, D. E. Genes involved in Drosophila glutamate receptor expression and localization. *BMC neuroscience* **6**, 44 (2005).
62. Featherstone, D. E. et al. Presynaptic glutamic acid decarboxylase is required for induction of the postsynaptic receptor field at a glutamatergic synapse. *Neuron* **27**, 71–84 (2000).

63. Robinson, S. W., Nugent, M. L., Dinsdale, D. & Steinert, J. R. Prion protein facilitates synaptic vesicle release by enhancing release probability. *Human molecular genetics* **23**, 4581–4596 (2014).

## Acknowledgements

We thank Pr. Bernhard Hovemann, Dr. Johannes Bischof, and Dr. Hermann Abele for providing fly stocks. We thank Serge Loquin and José Solonot for their technical help. We also thank Ingrid Adjovi and Isabelle Chauvel for conducting preliminary experiments concerning larval locomotor activity and immunostaining in Y.G.'s laboratory. Research in Y.G.'s laboratory is supported by the Centre National de la Recherche Scientifique, the European Union (ERC Starting Grant, GliSFCo-311403), the Agence Nationale de la Recherche (ANR-JCJC: GGCB-2010), the Conseil Régional de Bourgogne (Faber), and the Université de Bourgogne Franche-Comté. Research in D.E.F.'s laboratory was supported by the Muscular Dystrophy Association, and by the National Institutes of Health, National Institute of Neurological Disorders and Stroke (Grant Nr. R01NS045628).

## Author Contributions

A.Z. performed most of the cloning necessary to obtain the *JhI-21-Gal4* allele, JhI-21 immunohistology, behavioral experiments, analyzed ERC data, and contributed to the biochemical experiments. H.A. did the q-PCR experiments, the GluR immunohistology, and part of the electrophysiology related to sEJCs. N.C. produced and analyzed the ERC data. MBG measured JhI-21 developmental time, performed the FRT based JhI-21 deletion, and helped with the behavioral experiments. M.S. contributed to the biochemical experiments, and assisted with the JhI-21-Gal4 allele cloning. J.S. performed the electrophysiological experiments related to the difference between feeding and wandering larvae, and quantified the EJCs on NMJs. FG performed the biochemical experiments and cell culture to test JhI-21 glutamate transport ability. G.M. helped with the behavioral experiments. A.Z., D.F. and Y.G. designed the general scientific strategy and wrote the article with input from all the other authors.

## Additional Information

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing financial interests: The authors declare no competing financial interests.

**How to cite this article:** Ziegler, A. B. et al. The Amino Acid Transporter JhI-21 Coevolves with Glutamate Receptors, Impacts NMJ Physiology, and Influences Locomotor Activity in *Drosophila* Larvae. *Sci. Rep.* **6**, 19692; doi: 10.1038/srep19692 (2016).

<figure><img src="image_18.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/