Introduction

Neuromodulation provides the foundation for neural circuit flexibility, which in turn facilitates adaptive behavioral and physiological responses to ever-changing internal demands and environmental conditions³. Neuromodulators, such as neuro-peptides and biogenic amines, regulate neural circuits through their diverse effects on excitability and synaptic transmission^4-6. Many neural circuits are targeted by multiple neuro-modulators, which adds another layer of flexibility while also ensuring stability. Due to a limited number of in vivo model systems to investigate neuromodulation, most studies to date have used ex vivo, or reduced preparations, often investigating one modulatory system at a time^3,7. Here, we investigate in vivo the key modulatory inputs to Insulin-Producing Cells (IPCs), a group of approximately 14 neurosecretory cells in the pars intercerebralis of the Drosophila brain^8,9. This is particularly timely as recent advances of the Drosophila brain connectome provide an ideal platform to disentangle synaptic connections in an organism with a genetically accessible and numerically simple nervous system^10-12. However, the connectome is based on synaptic connectivity and cannot be used to uncover paracrine connections, which are a hallmark of neuromodulation.

IPCs are an ideal system to investigate principles of neuromodulation in vivo for several reasons: 1) IPCs release Drosophila insulin-like peptides (DILPs), which are analogous to human insulin. The insulin-signaling pathway and its functions are highly conserved across the animal kingdom, from worms over flies to humans^9,13. For instance, insulin signaling is essential for many vital processes including metabolic homeostasis, feeding and foraging, ageing, and growth^8,14-19. Hence, the IPC activity and resultant insulin release need to be carefully regulated on different timescales^20-23. 2) IPCs need to be modulated by different neuronal networks to adjust DILP release according to internal states, metabolic demands, and external inputs^19,23-25. 3) IPCs have been characterized in the brain connectome², facilitating the identification of their synaptic inputs. 4) The genetic toolkit available in Drosophila makes the IPCs not only accessible for in vivo recordings, but also allows for simultaneous manipulation of other neuromodulatory populations via optogenetics. The activation of these neuron populations leads to an intrinsic release of neuromodulators in a controlled manner, which we refer to as ‘intrinsic pharmacology’. 5) Their superficial location in the brain, in combination with genetic tools, enables IPC activity quantification in vivo with approaches like calcium imaging and whole cell patch-clamp recordings. Here, we combined these approaches to investigate how the activity of IPCs is modulated by functional inputs from a set of core aminergic and peptidergic neuron populations. Our results provide a framework for understanding the regulation and possibly dysregulation of an insulinergic system. In humans, dysregulation of the insulinergic system underlies diabetes, obesity, and other metabolic disorders^26,27. In Drosophila, all adult IPCs express DILP2, 3, and 5 and have traditionally been considered a homogeneous population²⁰. However, emerging evidence has indicated that the IPCs could, in fact, represent a heterogeneous population. Specifically, only a subset of IPCs is intrinsically mechanosensitive via Piezo channels, which are required to prevent food overconsumption²⁸. In addition, recording IPC responses to the stimulation of clock neurons or glucose revealed two to three subpopulations of IPCs^29,30. However, how the IPC population and individual IPCs are modulated by different neuronal pathways remains unclear.

In this study, we employed an exhaustive, multi-level approach to leverage Drosophila IPCs as a model for investigating neuromodulatory dynamics and principles of metabolic homeostasis in vivo. Using a combination of single-nucleus transcriptome analyses, anatomical receptor mapping, optogenetics, whole cell patch-clamp recordings, calcium imaging, and connectomics, we unravel how individual IPCs and the IPC population activity are affected by a broad range of modulatory inputs. We demonstrate that certain modulators shift the overall activity of the IPC population towards excitation or inhibition. This shift could in turn alter the responsiveness of the IPCs, possibly to prime the system for higher or lower metabolic demands. In addition to the effects on the overall populational activity, we unveil that IPCs consist of at least three subpopulations regarding their modulatory input. At the molecular level, these heterogeneous responses of individual IPCs can partially be explained by differences in their receptor expression profile. We hypothesize that this heterogeneity helps to buffer the system to allow insulin release to be tailored to the ever-changing demands of the animal. Taken together, our work demonstrates that IPCs are part of a complex neuromodulatory network and integrate various regulatory inputs via different pathways. Furthermore, the regulatory motifs identified in our study regarding the orchestration by multiple modulatory systems, are likely also present in other modulatory systems, allowing for the generalization of our findings.

Results

Single-nucleus transcriptomic analysis identifies novel modulators of IPCs

Previous studies using a combination of molecular techniques, functional studies, and ex vivo recordings have shown that several neuropeptides, biogenic amines, and classical neurotransmitters modulate IPCs^9,24,31. To identify additional modulators of IPCs and confirm previously published findings, we analyzed single-nucleus transcriptomes of IPCs, which were recently sequenced as part of the Fly Cell Atlas project³². Our analysis was based on 392 IPC transcriptomes, out of which 194 were derived from males and 198 from females. All these cells expressed Ilp2, 3 and 5 (Fig. S1A), the best-known marker genes for IPCs. Drosulfakinin (Dsk) neuropeptide, which has previously been shown to be expressed in a subset of IPCs^33,34, was very sparsely expressed in IPCs in this dataset (Fig. S1A). This sparse expression was independently confirmed using a Dsk-T2A-GAL4 knock-in line (Fig. S1B). Nonetheless, high expression of all three Ilp transcripts in these cells confirmed that these transcriptomes were derived from IPCs.

To catalogue all modulators potentially involved in regulating IPC activity, we examined the expression of receptors for all biogenic amines and neuropeptides (Fig. 1A) in IPC transcriptomes. Using a threshold of 5% expression, we found 22 neuropeptide and 17 biogenic amine receptors to be expressed in IPCs. This included receptors previously shown to be expressed in IPCs like leucokinin receptor (Lkr), allatostatinA receptor 1 (AstA-R1), and the serotonin receptor 5-Hy-droxytryptamine 1A (5-HT1A)^24,31,35,36. In addition, our analysis revealed several novel receptors, including 5-Hydroxytryptamine 7 (5-HT7) and 5-Hydroxytryptamine 1B (5-HT1B) receptors, rickets (rk), Trissin receptor (TrissinR), and others (Fig. 1A). Since there were little differences in receptor expression between the two sexes (Fig. S1C), we used the transcriptomes from males and females for all subsequent analyses.

Figure 1:

Interestingly, the insulin receptor (InR) was the highest expressed neuropeptide receptor in IPCs, suggesting that IPCs are modulated in an autocrine manner or by DILPs produced outside the IPCs. Moreover, IPCs express receptors for all classes of biogenic amines (Fig. 1A). The Drosophila genome encodes multiple cognate receptors for all biogenic amines and some neuropeptides. Interestingly, multiple receptors for some of these ligands (e.g., all biogenic amines, Allatostatin-A (AstA) and Diuretic hormone 44 (DH44)) are expressed in IPCs. Since there are functional differences between different receptors for a given ligand^37-39, co-expression of multiple receptors for a specific ligand in IPCs adds complexity to the modulatory input. In addition, IPCs can receive input via fastacting, small molecule neurotransmitters since they express receptors for gamma-aminobutyric acid (GABA), glutamate (Glu) and acetylcholine (ACh; Fig. S1D). Taken together, the gene expression data indicate that IPCs are intricately regulated by synaptic (via neurotransmitters) and paracrine/endocrine (via neuropeptides and biogenic amines) inputs.

IPCs are heterogeneous in their receptor expression profiles

IPCs have generally been regarded as a homogeneous population of approximately 14 cells that can be anatomically identified using antibodies against DILP2, 3 or 5. However, recent evidence suggests that IPCs can be separated into at least two subgroups based on the expression of the mechanosensitive channel Piezo in only a subset of IPCs²⁸ and heterogeneous responses to clock-neuron stimulation and indirectly to glucose^29,30. Our analysis confirmed that Piezo is only expressed in a subset of IPCs (Fig. S1A). Hence, we probed further to see if the receptors expressed in IPCs (Fig. 1A) were differentially expressed across IPCs, lending further support for the heterogeneity of the IPC population. For this, we performed unsupervised clustering of all IPC transcriptomes based on highly variable genes and visualized it using t-SNE analysis (Fig. 1B). This analysis retrieved four clusters that were not very distinct, suggesting a lack of drastic variability in gene expression across IPCs. Nonetheless, examining the expression of 10 neuromodulator receptors that were selected to represent neuropeptides and amines, as well as low and high receptor expression, revealed cluster-specific differences (Fig. 1C, 1D and S2). For instance, 5-HT1A was expressed in all four clusters, AstA-R2 was expressed in two clusters (Fig. 1C and 1D), while short neuro-peptide F receptor (sNPF-R) was only expressed in one cluster (Fig. 1D and S2). Other receptor genes, such as Piezo, were expressed across all clusters but at drastically different levels and in in different fractions of IPCs. Hence, differences in receptor expression profiles, based on transcriptome sequencing, support the conclusion that the IPCs are a heterogeneous population.

Anatomical mapping of neuromodulator receptors complements single-nucleus transcriptome profiling of IPCs

Single-nucleus transcriptomes can prove extremely powerful in examining global gene expression differences at a cellular resolution. However, it has limitations, including contamination by ambient RNA, which can give rise to false positives⁴⁰. To assess the reliability of our expression analysis, we verified the expression of 14 neuromodulator receptors in IPCs using anatomical techniques. We primarily used T2A-GAL4 knock-in lines^41,42 to drive expression of nuclear mCherry under the control of endogenous neuromodulator receptor promoters and then quantified the percentage of IPCs expressing mCherry for each receptor line using immunohistochemistry (Fig. 1E, 1F and S3-S5). In accordance with previous anatomical investigations^17,43-46 and our transcriptome analysis, receptors for AstA, LK, sNPF, and 5-HT were all expressed in the IPCs (Fig. 1E, 1F, S3 and S4). Importantly, the expression of several novel neuromodulator receptors identified by our sequencing analysis was also confirmed. For instance, all five serotonin receptor transcripts were detected in the IPC transcriptomes, one of which (5-HT7) was further confirmed using anatomical mapping (Fig. 1A, 1F, S3 and S4). These receptor mapping experiments also indicated that not all IPCs are identical, as some receptors were expressed in only a subset of IPCs (Fig. 1F and S4). For example, AstA-R2 was expressed in 86%, Lkr in 63%, and sNPF-R and myosuppressin (MS) receptor 1 (MsR1) in about 20% of IPCs. Thus, each receptor type is expressed by a different fraction of IPCs, implying that IPCs consist of more than two subgroups. Across individuals, however, the fractions of IPCs expressing a given receptor were consistent, suggesting that receptor expression is tightly regulated (Fig. 1F). We used percent expression similar and lower than corazonin receptor (CrzR) (Fig. 1F) as the upper bound to determine a receptor as not expressed in IPCs based on anatomical mapping (Fig. S5). This was justified based on the fact that the average CrzR expression (13.3%) was lower than one pair of IPCs (2 out of 14 or 14.3%) per brain, and previous studies reported a lack of CrzR expression in IPCs^29,47. Our anatomical mapping largely correlated with our single-nucleus sequencing analysis. Exceptions to this included the two MS receptors, MsR1 and MsR2. T2A-GAL4 lines for both receptors drove expression a small subset of IPCs but the expression of their transcripts in the IPC transcriptomes was below the 5% threshold used here (Fig. 1D-F, S3 and S4). Hence, neuromodulator receptors with low expression levels could be expressed in the IPCs in addition to the ones reported here based on the single-nucleus transcriptomic analysis. Other differences between the gene expression analysis and anatomical mapping were exemplified by the dopamine receptor Dop2R and the octopamine receptor Octα2R, which were both highly expressed in the IPC transcriptomes but were not detected using anatomical mapping (Fig. 1A, 1F and S5). To determine if these discrepancies were due to the presence of additional receptor transcript variants with different C-termini that may not be under the control of the T2A-GAL4 knock-in lines, we utilized a Trojan-GAL4 which should represent expression of other Octα2R transcript variants. Surprisingly, this general GAL4 line for Octa2R also did not drive expression in IPCs (Fig. S5). Given that G-protein coupled receptors (GPCRs) are usually expressed in low amounts, it is highly unlikely that these receptors represent false positives due to ambient RNA contamination during sample preparation for transcriptome sequencing⁴⁸. Hence, inadequate transgenes are the likely cause for the discrepancy between transcriptome analysis and anatomical mapping.

In summary, anatomical receptor mapping complemented the single-nucleus transcriptome analysis and confirmed the heterogeneous expression of several novel neuromodulator receptors in IPCs. Using two different methods to assess expression of receptors in IPCs also demonstrated the necessity of parallel approaches to obtain accurate insights into receptor expression in complex modulatory systems.

IPC activity is modulated by aminergic and peptidergic neurons

Previous studies and our receptor expression data suggested that the IPC activity could be modulated by several different modulatory systems acting in parallel, the combined activity of which could ultimately determine insulin release. However, the functional significance of these pathways remains unclear. Therefore, we next sought to investigate how the activation of specific modulatory populations affects IPC activity in vivo. To this end, we recorded the activity of individual IPCs via whole cell patch-clamp while optogenetically activating different populations of modulatory neurons optogenetically via CsChrimson (Fig. 2A and 1B). We used nine different driver lines targeting different populations of modulatory neurons (ModNs) for the activation experiments. Receptors for most of the neuromodulators found in these neurons were expressed in IPCs according to our expression analyses. We further included a driver line for neurons expressing adipokinetic hormone (AKH), since it was previously shown to influence DILP3 release in larvae⁴⁹. Hence, all nine of these ModNs have been suggested to affect IPC activity or play a role in metabolic homeostasis and feeding^{1,9,17,24,29,31,46,50-54}. These specific ModN driver lines targeted neurons containing the biogenic amines dopamine (DANs), octopamine (OANs), or serotonin (5-HTNs), as well as the neuropeptides AstA (AstANs), AKH (AKHNs), diuretic hormone 31 (DH31Ns), LK (LKNs), MS (MSNs), or tachykinin (TKNs, Fig. S6). We also included a line we expected to target sNPF-expressing neurons, but anatomical labeling revealed later that the line was in fact empty and, therefore, lent itself as an ideal negative control (Fig. S6, S7B and S7C). While some of these modulators are expressed and released by large populations of ModNs forming distinct sub-populations (e.g., AstA, DA, and 5-HT), others are expressed in very few, individually identifiable neurons (e.g., LK and AKH). Accordingly, the number of targeted neurons varied between driver lines (Fig. S6).

Figure 2:

During these intrinsic pharmacology experiments, we first recorded the baseline activity of IPCs, which had a membrane potential of about -40 mV (without subtraction of the 13mV liquid junction potential⁵⁵). The spike frequency was 0.8 Hz on average, in flies fed ad libitum. Activating individual ModN populations showed strong and immediate effects on IPC activity in several cases (Fig. 2C). To classify IPC responses using a data-driven approach, we implemented threshold criteria for the membrane potential to cluster the IPCs depending on the effects of ModN activation into excited, inhibited, or unaffected (Fig. 2C, for details see Material and Methods). Excitation was characterized by a fast depolarization of the membrane potential, often accompanied by an increase in spike frequency shortly after the onset of activation (Fig. 2C, magenta). Opposite to this, inhibition was characterized by a hyperpolarization of the membrane potential and a reduction in spike frequency (Fig. 2C, blue). Unaffected IPCs displayed no significant changes after activation compared to the baseline (Fig. 2C, black). For each recorded IPC, we activated one ModN population repeatedly for ten trials. Individual IPCs showed consistent responses to repeated activation (Fig. 2C and S7A). Using this approach, we quantified and categorized the effect of each ModN population on the IPC activity. We then analyzed the global effect of each ModN population on the IPC population activity. To this end, we calculated a general effect coefficient based on the mean area under the curve (AUC) of the membrane potential of all recorded IPCs after ModN activation (see Material and Methods for calculation details), describing a global trend for the shift in IPC activity (Fig. 2D). This coefficient was highest for DAN activation (0.51), meaning the overall activity of the IPCs was strongly shifted towards excitation. OAN activation also led to a positive coefficient (0.2), indicating a weaker excitatory shift. Activation of TKNs (0.01) and MSNs (−0.04) both resulted in a coefficient near zero and therefore did not seem to shift the overall activity. AKHNs (−0.05), 5-HTNs (−0.08), and DH31Ns (−0.09) shifted the activity of the IPC activity slightly towards inhibition. Activating AstANs (−0.39) had a stronger net inhibitory effect, and the activation of LKNs (−0.58) led to the strongest inhibition of the IPC population.

Next, we analyzed the effects of each ModN population on individual IPCs in detail. Activation of DANs and OANs either excited or had no effect on individual IPCs (Fig. 2E). DAN activation resulted in a strong excitation in the majority of IPCs (7/8), and only 1/8 IPCs remained unaffected. Activation of OANs also strongly excited most IPCs (7/10), while 3/10 IPCs were unaffected (Fig. 2E). In contrast, activation of MSNs and 5-HTNs either inhibited or had no effect on individual IPCs (Fig. 2F). However, in one IPC, there was reliably a single spike upon activation of 5-HTNs but the average membrane potential did not meet our threshold and therefore it was not labeled excited. Activation of MSNs led to the inhibition of 2/7 IPCs, while 5/7 remained unaffected. Activating 5-HTNs inhibited 4/9 IPCs and had no effect on 5/9 (Fig. 2F). Surprisingly, activating AstANs, AKHNs, DH31Ns, LKNs, or TKNs, had different and even opposite effects on different subsets of IPCs. Here, we found IPCs exhibiting all three response types (excited, inhibited, and unaffected) during activation of the same ModNs (Fig. 2G). When activating TKNs, the IPC response was distributed in 2/11 excited IPCs, 1/11 inhibited, and the majority with 8/11 remained unaffected. With the activation of AKHNs, only 1/8 IPCs was excited and 2/8 inhibited, while 5/8 remained unaffected. Activating DH31Ns had a stronger effect than AKHNs, with 2/8 excited IPCs, 3/8 inhibited, and only 3/8 remained unaffected. Interestingly, AstAN activation resulted in strong excitation and inhibition. Hence, we decided to analyze the effects of AstANs on IPCs in more detail by increasing the sample size. We found that 10/44 recorded IPCs were excited upon activation, while the largest subset with 32/44 were inhibited, and only 2/44 remained unaffected. LKNs showed an even stronger inhibitory effect, since the activation led to only 1/12 excited IPCs, while 11/12 were inhibited. We also recorded from IPCs while activating the line that was expected to express CsChrimson in sNPFNs but was in fact empty (Fig. S6 and 7D). Accordingly ‘activating’ this line while recording from IPCs had no effect on the IPC activity, except for one cell where the membrane potential increased slightly (Fig. S7D). To analyze the effect strength for all lines, we subtracted the baseline AUC, calculated before the activation onset, from the AUC after activation (ΔAUC). To determine if the effects were significant for the respective line, we calculated the p-values using the Wilcoxon signed-rank test for all recorded IPCs and for individual clusters with a minimum number of 5 recorded IPCs. Our analysis revealed that only the effects of DANs, OANs, AstANs, and LKNs were significant (Fig. 2H). For DANs and OANs, we found a significant effect for the population of all recorded IPCs and in the excited cluster. For LKNs, on the other hand, the IPC population and the inhibited cluster showed a significant effect. AstAN activation led to a significant effect in the IPC population, the excited and the inhibited cluster. TKNs, AKHNs, MSNs, 5-HTNs, and DH31Ns also had effects on the IPC activity, especially in the clusters we found for the respective line, but they were not statistically significant. Taken together, these results show that different populations of ModNs have distinct effects on the IPC activity. The heterogeneous effects the activation of individual ModN populations had on IPCs could arise from the heterogeneous receptor profiles exhibited by individual neurons in the IPC population. In addition, our intrinsic pharmacology approach does not only drive neuromodulator release, but also the release of classical, fast-acting transmitters contained in the ModNs. Therefore, IPC responses could be driven both by the respective neuromodulator and the co-release of fast-acting neurotransmitters possibly present in those neurons. This co-transmitter release occurs under naturalistic circumstances⁵⁶, and our IPC transcriptome analysis showed that IPCs express receptors for the classical transmitters ACh, Glu, and GABA (Fig. S1).

The high temporal resolution of our in vivo patch-clamp recordings enabled us to analyze the detailed temporal dynamics of IPC modulation. Here, we focused on DANs, OANs, AstANs and LKNs, which had strong, significant effects on the IPC activity (Fig. 2H). We quantified the latency from activation onset to IPC de- or hyperpolarization for each ModN population (Fig. 3A). All ModNs had rapid effects on the IPC activity, with DANs driving the fastest excitatory effects at a median latency of 10 ms, while OANs had the slowest with a median latency of 122 ms. Inhibitory effects of AstAN and LKN activation occurred with a similar median latency of 49 ms and 43 ms, respectively. Classical monosynaptic connections have a response latency of a few milliseconds⁵⁷, which is much faster than the latencies we observed. We therefore assume that the majority of the effects we observed were not due to direct, monosynaptic input to the IPCs. Instead, they are likely driven by polysynaptic or modulatory connections. Only the excitatory effects of DANs and the single excited IPC during LKN activation were on a timescale that could indicate monosynaptic connections.

Figure 3:

Heterogeneous IPC responses occur within the same animal

The heterogeneity we observed in the functional connectivity between ModNs and individual IPCs during patch-clamp recordings could have three different origins: 1) It could reflect differences in synaptic connectivity between different animals. 2) Functional connectivity between ModN populations and IPCs could change over time. For example, the heterogeneity could arise from changes in the animal’s internal state, including starvation or changes in behavior, such as increased leg flailing, which readily occurs during in vivo recordings. We controlled for some, but not all these aspects. For example, we provided food to all animals ad libitum until we prepared them for experiments, and we continuously perfusing the brain with the same, sugar-containing saline. 3) The heterogeneity could be a persistent feature of IPC modulation, allowing for a more nuanced shift of IPC population activity through the action of modulatory neurons on individual IPCs. This nuanced modulation might be beneficial considering that insulin release is critical for metabolic homeostasis and needs to be flexibly adjusted to various contexts with a relatively small population of IPCs.

To determine if the heterogeneity was a result of inter-individual differences in connectivity and/or internal states, we tested whether IPC responses to ModN activation were less heterogeneous when recording from multiple IPCs in a single animal, rather than inferring the population activity from recordings across animals. We therefore subsequently recorded from 11 IPCs in the same fly while activating AstANs. Remarkably, responses across the population of IPCs recorded successively in the same fly were just as heterogeneous as responses recorded across flies: 2/11 IPCs were excited, 7/11 inhibited, and 2/11 remained unaffected in the same fly (Fig. 3B and 3C). These heterogeneous responses and their distribution in a single animal mirrored the responses of IPCs recorded across multiple animals (Fig. 3B and 3C). Recording the activity of multiple IPCs in four additional animals showed similar heterogeneous response profiles (Fig. S8). Hence, heterogeneous IPC responses are not accounted for by inter-individual differences. Instead, the heterogeneity occurred within single animals and appears to be a persistent feature of IPC modulation.

Modulation of individual IPCs underlies shifts in population activity

Our recordings from multiple IPCs in the same animal revealed that the heterogeneous responses are not based on functional differences between individuals. In principle, however, the het-erogeneity could still arise from changes in connectivity over time, for example due to internal state changes. Therefore, we next recorded from multiple IPCs simultaneously via in vivo widefield calcium imaging, to address whether excitatory and inhibitory responses were elicited at the same time in the same animal and to gain a better understanding of the IPC population dynamics (Fig. 4A). To this end, we recorded the calcium dynamics in IPCs during the activation of DANs, OANs, or MSNs. We chose these lines since they generated the full spectrum of responses from primarily excitatory to inhibitory in patch-clamp recordings. In addition, we selected lines targeting AstANs and LKNs, which elicited heterogeneous effects in our patch-clamp recordings. We tested different activation durations to account for the slower dynamics of the calcium indicator (Fig. S9A). A 5 s long activation seemed most appropriate, since it elicited consistent IPC responses while ensuring that IPCs retained naturalistic activity patterns throughout the duration of the experiment (Fig. 4A-C, and S9B). Activation of the selected driver lines evoked strong and immediate fluorescence changes across the IPC population (Fig. 4 C). When comparing IPC responses during two consecutive activations, we found that the same subset of IPCs showed a fluorescence increase during both activations while other IPCs consistently showed a decrease (Fig. 4D). Hence, individual IPC responses were consistent across trials, just like in our electrophysiology experiments. Therefore, we performed a threshold-based analysis to cluster the responses of individual IPCs into excited, inhibited, and unaffected, analogously to the approach used for electro-physiology data (Fig. 4E, see for details Material and Methods). Imaging multiple IPCs simultaneously revealed that the three activity clusters occurred within the same animal, and at the same time (Fig. 4E). We visualized the ΔF/F traces for each recorded IPC for all animals in the respective ModN line as heatmaps (Fig. 4F). These high-resolution heatmaps depict the effects of the respective ModN activation on all recorded IPCs, across all animals. For example, the activation of AstANs led to a fluorescence increase in the majority of IPCs across all animals, while the activation of MSNs decreased the fluorescence in most IPCs in all tested animals (Fig. 4F). However, we also found differences between individual animals. For example, during the activation of OANs, IPCs in flies #2, 5, and 6 were primarily excited, while IPCs in fly #4 were primarily inhibited. Interestingly, we found heterogeneous IPC responses across all driver lines, independent of the strength or valence of the global trend of the IPC population (Fig. 4F). Recording from multiple IPCs at the same time uncovered that within one cluster, the response seemed homogeneous, e.g., in fly #4 during the activation of DANs (Fig. 4F).

Figure 4:

The threshold-based cluster analysis revealed three clusters based on the responses of the IPCs to the activation of the respective driver line for all ModNs tested (Fig. 4G). To further quantify the effects of ModNs on IPCs and to assess whether they were statistically significant, we compared the activity of each IPC before activation to the activity after activation based on the AUC of individual ΔF/F traces (Fig. 4H). First, we performed this analysis across the entire IPC population for each driver line. Second, we analyzed the IPCs by cluster, i.e., the excited, inhibited, and unaffected IPCs separately. AstANs, OANs, and DANs had significant overall effects on the population activity, driving a net increase in IPC activity (Fig. 4H, gray). For these ModNs, the effects were also significant when considering the inhibited and excited cluster separately (Fig. 4H, blue and magenta). MSN activation significantly inhibited the IPC population activity, an effect that was carried solely by a significant effect on the inhibited IPC cluster (Fig. 4H, blue). LKN activation had significant effects in the excited and inhibited clusters, but those canceled each other out in the pooled data, resulting in no net excitation or inhibition (Fig. 4H, gray). The global effect of each ModN population on the IPCs is summarized in Fig. 4I. This analysis demonstrated the effect strength during the activation of the respective ModNs, with AstANs, OANs, and DANs shifting the population towards an excited state, while MSNs shifted the population activity towards inhibition (Fig. 4I). These shifts likely change the responsiveness of the IPCs to additional inputs, with AstAN, OAN, and DAN input increasing the responsiveness of the system, while MSN dampens it.

In addition to the population-level analysis, calcium imaging allows the investigation of receptor-dependent calcium dynamics that are not coupled to electrical activity. For example, intracellular calcium mediates changes in gene expression in addition to transmitter release⁵⁸, which might underlie longterm changes in IPC activity. Our intrinsic pharmacology experiments during calcium imaging revealed that activating AstANs and LKNs caused calcium dynamics that were not necessarily predicted from electrophysiological data. For example, in both cases, we observed strong increases in calcium in the IPCs, which occurred in addition to the inhibition recorded via patch-clamp. Especially AstAN activation resulted in a fast and robust calcium influx across the majority of recorded IPCs. This calcium increase could lead to an increase in insulin release, but also to alterations in gene transcription, for example. To investigate how the electrophysiological activity recorded in patch-clamp translates into calcium signals in IPCs, we recorded the spontaneous baseline activity from one IPC via patch-clamp while performing calcium imaging of 13 IPC cell bodies (ROIs) without activating any ModNs. When comparing the membrane potential with the fluorescence traces of the same IPC, fluorescence increases matched the occurrence of individual action potentials faithfully (Fig. 4J, gray boxes and purple traces). Accordingly, bursts of action potentials led to a prolonged and stronger calcium influx (Fig. 4J asterisks).

By superimposing the action potentials and the associated fluorescence traces, the close correlation became even more apparent (Fig. 4K), with the expected delay for calcium imaging rooted in the kinetics of GCaMP6m⁵⁹. However, when inspecting the membrane potential and the fluorescence trace of the patched IPC more closely, we noticed that in some instances calcium influxes were not driven by action potentials, but by subthreshold depolarization (Fig. 4J, arrowheads). Hence, while there was a strong correlation between IPC spiking and calcium influx, calcium dynamics were not solely coupled directly to action potentials. Interestingly, when comparing the fluorescence traces of all 13 IPCs with each other, we found that the fluorescence dynamics were synchronized between some, but not all IPCs (Fig. 4J). This suggests that some IPCs either receive the same input or form reciprocal connections, while others receive different inputs and are not coupled to the same sub-population. This further demonstrates that IPCs form distinct clusters.

IPCs receive heterogeneous synaptic inputs

After investigating the functional inputs to IPCs, we next analyzed direct synaptic connections. To this end, we used data from the FlyWire whole-brain connectome^10,12,60-62 to identify presynaptic neurons that connect to IPCs² via at least 5 synapses (see Materials & Methods for more details). These neurons were predicted to contain either ACh, Glu, GABA, or unknown transmitters (Fig. 5 A – D). ‘Unknown’ refers to neurons for which transmitter prediction scores were low, and therefore inconclusive⁶¹, which are primarily modulatory neurons. We grouped the presynaptic neurons according to their neurotransmitter identity and calculated the overall synaptic input each IPC receives from neurons of each group (Fig. 5E). This analysis revealed that the IPCs receive heterogeneous inputs from three classical transmitters, which reflects their heterogenous modulatory inputs. This heterogeneity is most apparent in GABAergic connections, since only eight IPCs receive inputs from neurons expressing GABA. Presynaptic neurons containing ACh or Glu provide input to all IPCs. However, the synaptic strength of the individual connections was heterogeneous, ranging from 7 synapses to 63 for cholinergic and 8 to 36 synapses for glutamatergic inputs. Even though we cannot define the transmitters released by neuron in the ‘unknown’ category, they also form strongly heterogeneous connections to IPCs. The heterogeneity in synaptic inputs likely contributes to the heterogeneous responses of the IPC population we observed in functional connectivity experiments.

Figure 5:

Taken together, our intrinsic pharmacology experiments showed that the different modulatory populations had immediate effects on the IPC activity. In calcium imaging experiments, the overall trends of the IPC population activity ranged from a strong calcium increase during the activation of AstANs to mostly inhibited during the activation of MSNs (Fig. 6A). For DANs, OANs, and MSNs the activation led to comparable global shifts in the population activity between calcium imaging and patch-clamp. DANs and OANs had a strong excitatory effect, while MSNs had a slightly inhibitory effect. For AstANs and LKNs, calcium imaging revealed additional activity dynamics compared to patch-clamp recordings. In patch-clamp recordings, AstANs elicited a strong inhibitory effect on the majority of IPCs and a strong excitation in a smaller subset. In calcium imaging, AstAN activation led to a strong calcium increase in the majority of IPCs. LKN activation also had a stronger inhibitory effect in patch-clamp recordings compared to calcium imaging experiments, which is apparent when comparing the global trends of both lines between methods (Fig. 6A). Those similarities and differences in the overall effect were not only driven by the effect strength, but also by the number of affected IPCs. For DAN, OAN, and MSN activation, the number of cells per cluster was similar in patch-clamp and imaging, while for AstANs, and LKNs the affected fraction of IPCs per cluster differed considerably (Fig. 6B). On a single cell level, we found heterogeneous responses while imaging multiple cells at the same time during activation of the same ModN population. This is a strong indicator that the response clusters we observed in patch-clamp recordings are not due to state-dependent shifts between recordings. Instead, the calcium imaging recordings indicated that the heterogeneous responses are a persistent feature of the IPC population. Within any given cluster, the responses appeared highly homogeneous, indicating that within one response cluster the modulatory input to the associated IPCs was similar.

Figure 6:

To link our receptor expression analyses with the valence of the connectivity observed in electrophysiological recordings and calcium imaging, we analyzed the G-protein coupling scores for the receptors of interest, including multiple splice isoforms where relevant (see Table S2 for details; http://athina.biol.uoa.gr/bio-informatics/PRED-COUPLE2/). To simplify comparisons, we predicted that a neuromodulator could be excitatory if at least one of its receptors (including transcript variants) coupled to Gs or Gq proteins. Similarly, a neuromodulator was predicted to be inhibitory if at least one of its receptor variants was predicted to couple to Gi protein. When comparing the predictions and anatomical mapping results with the results from patch-clamp recordings and calcium imaging, it became apparent that the majority of ModN activation effects were correctly predicted by the receptor profile (Fig. 6C). For example, AstA receptors were predicted to signal via both excitatory and inhibitory secondary messenger pathways, which could explain the heterogeneity and the differences observed between patch-clamp experiments and calcium imaging. On the other hand, we expected that LKN activation would primarily result in IPC excitation, but this contrasts with our patch-clamp and calcium imaging experiments. Hence, the effects of LKN activation on IPCs are likely due to other factors, e.g., indirect connections between LKNs and IPCs, and potentially other neurotransmitters co-released by LKNs upon activation. However, for the majority of investigated ModNs, their global effect on the IPC activity can be predicted based on receptor expression and G-protein coupling, which can also account for the heterogeneous effects.

Discussion

Our results show that the activity of IPCs in Drosophila is affected by multiple neuromodulatory inputs. We analyzed the receptor profile of IPCs and the effects of nine key ModN populations on the IPCs on the level of single-cell responses and population activity shifts (Fig. 6C and 6D). Hence, our study provides a multi-level quantification of the complex regulation of IPCs, which underlies their ability to respond to changing metabolic demands and sensory inputs.

General effects of ModNs on IPC activity

Neuromodulation enables multiplexing in neuronal circuits and lends flexibility to seemingly hard-wired connections. In contrast to classical, fast-acting neurotransmitters, neuromodulators can act far away from their release sites, and thus change neuronal circuit dynamics across the entire nervous system. Due to the complexity of neuromodulatory systems, studies examining their functions or modes of action have often been limited to ex vivo preparations or in vitro physiological assays⁶³. Here, we leveraged the IPCs of Drosophila as a model system to investigate neuromodulation by combining multiple approaches, including in vivo patch-clamp recordings and calcium imaging, as well as receptor expression mapping and connectomics. Our results revealed that the IPC population activity can be shifted towards excited or inhibited states on fast timescales through the action of different ModN populations. Hence, some ModNs have strong, immediate effects on the global IPC activity, which could prime the system for pending metabolic demands by altering the excitability of IPCs and ultimately affect insulin release. Prolonged activation protocols during patch-clamp recordings and calcium imaging experiments demonstrated that these effects persist during continued activation (Fig. S7A, S7B and S9A). This indicates that the input of the respective ModNs could also affect the IPC activity on longer timescales. These effects on long timescales could in turn affect the IPCs on two levels: 1) sustained excitation could lead to a direct increase in insulin release, and 2) continuous activation of IPCs through modulatory inputs could lead to altered gene expression, e.g., resulting in increased expression of ilps^64,65. Prolonged inhibition could have the opposite effects, with a decrease of insulin release and ultimately a decrease in ilp expression. In both scenarios, the prolonged modulation of IPC activity and therefore insulin release could alter nutritional state-dependent behaviors, e.g., foraging behavior³⁵.

Although most ModNs investigated in our study coexpress other signaling molecules, e.g., additional neuropeptides and fast-acting neurotransmitters⁵⁶, we focus our discussion on functions of the neuromodulators that define the populations, especially in relation to insulin signaling. For example, activating DANs and OANs shifted the IPC population activity towards an excited state in patch-clamp recordings and calcium imaging. In Drosophila, DA is implicated in modulating locomotion, learning, and courtship behavior^66-68. In female flies, DA directly acts on the IPCs via DopR1, which affects ovarian dormancy via effects on DILP signaling⁶⁹. Our results demonstrate that DANs have a fast and strong excitatory effect on the IPC activity in vivo, on short time scales. However, ovarian dormancy is accompanied by diminished insulin signaling and is likely controlled over longer time scales. Hence, there might be an additional long-term effect of dopamine on IPCs, which we have not explored here.

Our finding that OANs had strong, fast, excitatory effects on IPCs in patch-clamp recordings and calcium imaging, was somewhat unexpected. This is because the actions of OA and insulin are antagonistic in several contexts. For example, OA drives increases in locomotor activity while insulin decreases locomotion^25,70-72. Accordingly, OA is released during increased physical activity, such as flight and walking, while insulin release is inhibited^23,73. Moreover, knockdown of the OAMB receptor in IPCs increases Dilp2 expression^51,74. However, our receptor expression analysis revealed that IPCs express a variety of OA receptors besides OAMB. In addition, OA has been reported to have varying effects in foraging behavior. While it has been demonstrated that OA is in general mediating hyperactivity in starved flies^25,70, activating a specific octopaminergic neuron type suppresses odor tracking during foraging⁷⁵. Similarly, reduced insulin levels increase odor-driven food-searching behavior^76,77. Therefore, OA is not always acting antagonistically to insulin, which means that context-dependent differences in the effect of OA on IPCs are conceivable. Our G-protein coupling predictions indicate that the expressed OA receptors could be excitatory and inhibitory, hence activation of OA receptors could lead to heterogeneous responses in IPCs. In addition, neurons containing OA co-express other transmitters, in particular tyramine, Glu, GABA, and ACh⁵⁶. Therefore, activating OANs in our experiments could lead to the release of several transmitters, which could bind to a variety of receptors resulting in an overall excitatory effect. Importantly, since we activated OANs intrinsically, these cotransmitters are likely also released under natural circumstances.

In contrast, MSN and LKN activation led to an inhibition of the IPC activity in patch-clamp and calcium imaging experiments (Fig. 6B). MS is known to have relaxant properties on insect muscles and leads to crop enlargement in female flies during food consumption after mating⁷⁸. Thus, MS enables the consumption of larger meals, which is essential for mated females to account for their higher energy demands. This could be supported by the concurrent IPC inhibition we observed, since this would result in lowered insulin levels and therefore a higher feeding drive³⁵. Insulin signaling also contributes to the modulation of mating in flies. Reduced insulin levels lead to a reduction of re-mating in mated females. Hence, the inhibitory effect of MSNs on IPCs could also act as a feedback signal to reduce premature re-mating⁷⁹. In addition, AKH-producing cells also express receptors for MS⁸⁰. This suggests that MS could be regulating metabolic homeostasis via the modulation of insulin and AKH in the context of reproduction, where the inhibition of IPCs is one of several parallel effects mediated by MS to modulate post-mating behaviors in females.

The overall inhibitory effect of LKNs on the IPCs was more pronounced in patch-clamp recordings than in calcium imaging, which is partially due to the higher sensitivity of electrophysiological recordings to membrane hyperpolarization. In calcium imaging, the IPC population activity coefficient was slightly positive, while it was strongly negative in patch-clamp recordings (Fig. 6A). However, the cluster distribution shows that the inhibited IPC clusters contained the highest numbers of IPCs during LKN activation in both methods, meaning that the majority of IPCs were inhibited (Fig. 6B). IPCs have been shown to be anatomically and functionally downstream of LKNs^17,46. For example, Lkr knockdown in IPCs affects starvation tolerance and sleep. Specifically, a subset of LKNs increases their activity during starvation and are required for starvation-induced suppression of sleep⁴⁶. Insulin is also implicated in the modulation of sleep, since a loss of DILPs, especially DILP2 and DILP3, reduces sleep depth and duration^81,82. Our study supports those findings by showing that LKN activation, which would occur during starvation under normal conditions, has a strong inhibitory effect on IPCs. The resulting decreased DILP release is expected to cause a reduction of sleep, mimicking the effects of starvation mediated by LKNs.

Interestingly, with AstANs we observed different effects when comparing patch-clamp recordings with calcium imaging. In patch-clamp recordings, AstAN activation had a strong inhibitory effect on the majority of IPCs. AstA is known to inhibit feeding and promote sleep^24,44,83. Therefore, the inhibitory effect of AstAN activation on IPCs could lead to an immediate decrease in DILP release, which is then followed by a long-term reduction in IPC activity resulting from the AstA-driven cessation of food intake. However, we also recorded a fast and strong excitation in a small subset of IPCs during patch-clamp experiments. In addition, we observed a strong calcium influx in the majority of IPCs upon AstAN activation in calcium imaging experiments. Our G-protein coupling predictions can explain these heterogeneous effects, but also the differences between patch-clamp and calcium imaging. AstAR1 and AstAR2 are both predicted to couple to Gi proteins, supporting the prediction of an inhibitory effect of AstA. However, AstAR2, which is almost exclusively expressed in IPCs within the brain, is also predicted to couple to Gq subunits. Gq subunits activate cellular pathways that lead to an intracellular calcium release that subsequently causes calcium channels in the membrane to open and the influx of extracellular calcium⁸⁴. Therefore, the presence of AstAR2 receptors opens the possibility of a calcium increase in parallel to the overall inhibition of the IPCs (see Table S2 for more details). The heterogeneity resulting from AstAN activation in both patch-clamp and calcium imaging can therefore be explained by heterogeneous expression of receptors across individual IPCs since AstAR1 and AstAR2 are not expressed in all IPCs. In addition, AstAN activation could also lead to excitatory input via co-expressed transmitters. Previous single cell transcriptome analyses indicated that at least some AstANs co-express ACh and glutamate^56,85 and our expression analyses revealed that IPCs express receptors for these transmitters (Fig. S1C). In addition, our connectome analysis revealed that all IPCs can receive input from putative cholinergic and glutamatergic neurons (Fig. 5). Moreover, indirect pathways from the AstANs to IPCs via interneurons could also drive excitatory responses, especially considering the long average latency of the excitatory effect (Fig. 3A).

Taken together, the differences we observed for AstAN activation between patch-clamp recordings and calcium imaging emphasize the importance of utilizing parallel approaches to unravel different layers of functional effects caused by modulatory inputs to IPCs. Our receptor expression analysis complements the electrophysiological data by providing valuable information about the foundation for the heterogeneous effects we observed in the functional connectivity.

IPCs form a heterogeneous population that is differentially modulated by sets of neuromodulators acting in parallel

Since the IPCs form a cluster of morphologically indistinguishable neurons and all of them express DILPs 2, 3, and 5, they have traditionally been regarded as a homogeneous population. Recently, this dogma has been challenged by a study that found the mechanosensitive channel Piezo in only a subset of IPCs²⁸. We found additional proof of heterogeneity within the IPC population based on receptor expression analyses and IPC responses to ModN activation. Intriguingly, cellular and functional heterogeneity also occurs in human pancreatic β-cells. This heterogeneity is assumed to be a mechanism that allows for dynamic transitions between states of high insulin production and cellular stress recovery⁸⁶, which could be a general feature of insulin releasing cells.

In our study, single-nucleus transcriptomes of IPCs and anatomical receptor mapping experiments revealed that modulatory and other receptors are not uniformly expressed across the IPC population. Instead, most receptors are only expressed in subsets of IPCs. Activating ModN populations while recording the IPC responses via patch-clamp revealed that the heterogeneous receptor profiles of individual IPCs have functional significance. Calcium imaging experiments for selected ModNs revealed that the same modulatory input led to simultaneous, heterogeneous responses across the IPC population, within the same animal. Therefore, the heterogeneity is not based on internal state changes. Finally, analyzing connectomics data revealed that IPCs also receive heterogeneous inputs via direct synapses, both regarding input profile and input strength. The question that arises from our results is whether we investigate stochastic heterogeneity, or whether the heterogeneity is a control mechanism for insulin release. Since we demonstrated similar heterogeneous receptor expression as well as similar responses to a specific modulatory input across animals, stochastic heterogeneity, as for example demonstrated in the development of photoreceptors of Drosophila⁸⁷, is unlikely. For example, over 88% of IPCs expressed at least one receptor for dopamine, which matches our recordings in patch-clamp and calcium imaging closely, where the activation of DANs had an effect on 87% of IPCs. In addition, the G-protein coupling analysis accurately predicted a strong excitation of the IPC population accompanied by a weak inhibition by DA.

Since the heterogeneity occurs on the receptor level and in physiological responses, it is likely to be beneficial for the system. We hypothesize that heterogeneous activation of individual IPCs facilitates a more nuanced insulin release and dynamic transitions between states of high insulin demand and recovery, comparable to the heterogeneity found in human pancreatic β-cells. A fine-tuned release from the IPCs is especially crucial when considering that insulin does not only regulate metabolic homeostasis, but also acts as a neuromodulator in the brain^77,88,89. In the olfactory systems of flies and rats, insulin reduces food-searching behavior in satiated animals⁸⁸ through the inhibition of a specific subset of olfactory sensory neurons⁷⁷. Therefore, insulin release needs to be adjusted to metabolic demands, and for example be reduced in foraging animals⁹⁰. We have three hypotheses for how the modulation of IPCs could alter their activity dynamics and ultimately insulin release: 1) Heterogeneous excitation and inhibition of IPCs could lead to a direct adjustment of the net output by increasing or decreasing insulin release by individual IPCs. 2) We revealed that on top of the heterogeneity, the population activity was shifted towards excitation or inhibition. Shifting the IPCs into a more excited state could increase the responsiveness of the system to enable a fast response to a secondary excitatory input (Fig. 6E), leading to a stronger release of insulin, for example in the context of feeding. Receiving an additional inhibitory input, on the other hand, would only lead to a moderate insulin release (Fig. 6E). In the same fashion, shifting IPCs towards inhibition could make the system less sensitive to small changes on the input level in a generally low insulin demand state. This would enable a moderate insulin release upon arrival of a secondary excitatory input (Fig. 6E). On the other hand, an additional inhibitory input would lead to a drastic shut down of IPC activity (Fig. 6E) as observed during locomotion²³. 3) The heterogeneous input could regulate the ratio of DILP2, 3, and 5 release. For example, the secretion of either DILP2 or DILP3 is selectively increased depending on the dietary conditions⁴⁹. Moreover, DILP3 expression, but not DILP2 or 5, has been tightly connected to the activity of AstANs⁴⁴.

Generally, heterogeneity adds flexibility and stability to neuromodulatory systems. In the case of IPCs, this could allow a relatively small number of neurons to integrate a large variety of inputs to achieve flexible insulin release tailored to the current internal demands.

Competing interest statement

The authors declare no competing interests, financial or otherwise.

Author contributions

Writing – original draft, M.H., M.Z., and J.M.A.; writing – review & editing, M.H., R.S.B., M.Z., A.S.C., I.S.B., S.H., S.L., F.C.-M., and J.M.A.; patch-clamp recordings, R.S.B., A.S.C., S.L., F.C.-M.; calcium imaging, M.H. and I.S.B.; receptor expression, M.Z. and S.H.; anatomy, M.H., R.S.B., and S.H.; connectomics analysis, M.H.; connectomics IPCS identification, M.Z.; software, M.H., M.Z., F.C.-M., and J.M.A.; formal analysis, M.H., R.S.B., M.Z., I.S.B., F.C.-M., and J.M.A; supervision, M.Z. and J.M.A.; project conceptualization, M.H., M.Z., and J.M.A; project administration, J.M.A; funding acquisition, M.Z. and J.M.A.

Acknowledgements

We thank Haluk Lacin (University of Missouri-Kansas City) for providing fly-lines for activation experiments, Jan A. Veenstra (University of Bordeaux) for providing the DILP-2 antibody, and Susan Morton for providing the RFP antibody. We are grateful to Konrad Öchsner for technical assistance and to Charlotte Helfrich-Förster and Wolfgang Rössler (all Julius-Maximilians-Universi-tät of Würzburg) for sharing resources; We also thank Nicole Enslinger for pre-liminary calcium imaging experiments and Emilia Derksen (both JMU) for technical support. We are grateful to Dick R. Nässel (Stockholm University), Tanja A. Godenschwege (Florida Atlantic University), Mattias Alenius (Umeå University), and Chris J. Dallmann (JMU) for valuable feedback on the manuscript. We thank the Princeton FlyWire team and members of the Murthy and Seung labs, as well as members of the Allen Institute for Brain Science, for development and maintenance of FlyWire (supported by BRAIN Initiative grants MH117815 and NS126935 to Murthy and Seung). We also acknowledge members of the Princeton FlyWire team and the FlyWire consortium for neuron proofreading and annotation.This work was supported by a grant from the Deutsche Forschungsgemeinschaft (DFG) to J.M.A. via the Emmy Noether program (DFG AC 371/1-1), by a grant from the DFG to J.M.A. as part of the NSF/CIHR/DFG/FRQ/UKRI-MRC Next Generation Networks for Neuroscience (Neuronex) Program (DFG AC 371/2-1), by a grant to J.M.A. from the Bavaria California Technology Center (BaCaTec), and by a grant from the DFG to M.Z. (DFG; ZA1296/1-1).

Materials and Methods

Fly husbandry

All flies were reared on standard fly food (cornmeal-agar-molasses medium) and kept in a climate chamber at 60% relative humidity, 25°C, and a 12h:12h light/dark cycle. Experiments were performed on mated females 3 – 6 days after eclosion. For experiments with optogenetic activation of CsChrimson, 300 µM all-trans-retinal (R2500, Sigma-Aldrich, Steinheim, Germany) was added to standard fly food and experimental flies were kept on that food after eclosion in darkness until usage.

Fly genotypes

Table 1 shows all genotypes used for anatomical receptor mapping and Table 2 shows the genotypes used for activation experiments as well as brain and VNC anatomical labeling.

Table 1:

Table 2:

Immunohistochemistry

Anatomical receptor mapping

Flies of the F1 generation from crosses with GAL4 lines from Table 1 and UAS-NLS-mCherry were used for anatomical receptor mapping. We primarily used the recently generated receptor-T2A-GAL4 knock-in lines ^41,42 to drive expression of nuclear mCherry. Since the T2A-GAL4 knock-in lines were generated by inserting the T2A-GAL4 coding sequence right before the stop codon for the gene, an accurate expression can be assumed since the GAL4 is expressed in the same temporal and spatial patterns, and at similar levels, as the native receptor. Flies were prefixed at room temperature for 2.5 hours in a 2 mL tube with 1.5 mL 4% paraformaldehyde and 20 μL 0.5% PBT (phosphate-buffered saline with 0.5% TritonX-100; Sigma-Aldrich, Steinheim, Germany). The fixed flies were washed 4 × 15 minutes each with 0.5% PBT. Brains were dissected afterwards and collected directly in 400 μL blocking solution (5% NGS in PBT) and incubated for 1.5 hours to block non-specific binding sites. Next, samples were incubated in primary antibodies (Table 3) diluted in blocking solution for 48 hours at 4°C. Afterwards, the brains were washed 4 × 15 minutes each in PBST at room temperature before they were incubated in the secondary anti-bodies (Table 3) for 48 hours at 4°C. Next, the samples were washed for 4 × 15 minutes in PBT, followed by two final washing steps in PBS for 15 minutes each before they were mounted in Fluoro mount-G (Invitrogen, Waltham, MA, USA). Z-stacks were acquired using a Leica TCS SPE confocal microscope (Leica Microsystems, Wetzlar, Germany) equipped with Leica Microsystems LAS AF software, a photo-multiplier tube, and 488, 532, and 635 nm solid-state lasers for excitation. We used a 20× glycerol immersion objective for whole-mount scans and 40× oil immersion objective for detailed images. All images were captured at a resolution of 1024×1024 pixels. To reduce background noise, all frames were captured two to three times and then averaged. Images were analyzed with ImageJ.

Table 3:

Table 4:

Expression analysis of driver lines used for activation experiments

To extract the brain and VNC from flies of the driver lines we used in the opto-genetic activation experiments, flies were anesthetized on ice and then fixed in 4% PFA in 0.1 M PBS for 3 hours followed by three subsequent washes in 0.1 M PBS. They were then dissected in 0.1 M PBS at room temperature. Prior to dissection, each fly was shortly submerged in 70% ethanol for ∼30s. The dissected brains were then washed in PBT for 3 x 15 minutes. To block non-specific binding sites the brains and VNCs were pre-incubated in 10% Normal Goat Serum (NGS) in PBT (blocking buffer) overnight at 4°C. The primary antibodies (see Table 3) were diluted in blocking buffer and this solution was then used to incubate brains and VNCs for 48 hours at 4°C. Next, the tissues were washed 3 × 15 min. in PBT and, then incubated overnight at 4°C (or for 3 hours at RT) with the secondary antibodies diluted in blocking buffer (see Table 3). After washing for 3 × 15 min. with PBT, the solution was replaced with Vectashield Antifade Mounting Medium (H-1000, Vector Laboratories, CA, USA) and the brains and VNCs were mounted on adhesion slides (Polysine slides, Thermo Scientific, Waltham, MA, USA). After covering the samples with Vectashield under a coverslip, the edges were sealed with nail polish and stored at 4°C until images were taken. Fluorescence images were acquired with a confocal laser scanning microscope (CLSM, Leica TCS SP8 WLL) via the Leica Application Suite X (LAS X, Leica Microsystems) using HC PL APO 10×/0.4, HC PL APO 20×/0.75 IMM objectives. Fluorescence signals were detected in serial stacks with a resolution of 1024 × 1024 pixels using 3 channels in sequential scanning mode. The GFP signal, which labeled the respective MN population of the GAL4 lines in Table 2 was acquired using an excitation laser wavelength of 488 nm and a HyD detector tuned to 498 – 540 nm. For Alexa Fluor 555 visualization, which was used to label the IPCs via the DILP2 antibody staining, we excited the fluorophore with the laser set to 555 nm and detected the emission with a PMT tuned to 565 – 620 nm. In the third channel, Alexa Fluor 635 bound to the nc82 antibody was excited with a wavelength of 633 nm and the emitted light was detected by a second PMT tuned to 643 – 693 nm. The image stacks were saved as a lif files from which maximum projections and scale bars were created using the Leica Application Suite X. Images were analyzed with ImageJ.

In vivo whole cell patch-clamp recordings

Fly handling and preparation were performed as described in ⁹³. All experiments were performed at room temperature under daylight conditions. Before preparation, flies were cold anesthetized at 4°C for immobilization and the head and thorax were positioned in a custom-made shim plate fly holder and mounted with UV glue (Proformic C1001, VIKO UG, Munich, Germany). After removing the forelegs, the remaining stumps and the proboscis were glued together for maximal stability. Afterwards, a window was cut into the posterior cuticle of the head. Trachea, fat tissue, and the ocellar ganglion were removed to expose the brain and gain access to the IPC cell bodies in the pars intercerebralis (PI). The fly was then transferred with the fly holder into a fixed stage fluorescence upright microscope (Olympus BX51WIF, Olympus Corporation, Japan) or a customized Slicescope (Scientifica, Uckfield, UK). Live images of the fly brain were acquired with a high-resolution camera (SciCam Pro, Scientifica, US) using an image acquisition software (OCULAR, Digital Optics Limited, Auckland, NZ). The brain was continuously perfused with carbonated (95% O2 and 5% CO2) extracellular saline containing 103 mM NaCl, 3 mM KCl, 5 mM N-[Tris(hydroxymethyl)methyl]-2-aminoethanesulfonic acid, 10 mM D-(+)-glucose, 8 mM D-(+)-Trehalose dihydrate, 26 mM NaHCO3, 1 mM NaH2PO4, 1.5 mM CaCl2.2H2O, 4 mM MgCl2.6H2O, and osmolarity adjusted to 273-275 mOsm ⁵⁵ during the preparation and the experiment. For enzymatic dissociation of the neural sheath from the brain above the IPCs, 0.025 % (w/v in extracellular saline) collagenase (Sigma-Aldrich #C5138) was applied locally in the PI region using a blunt thin-walled pipette. After identifying the GFP-expressing IPC cell bodies, whole cell patch-clamp recordings were performed using thick-walled patch pipettes (4 – 8 MΩ resistance) containing intracellular saline (140 mM potassium aspartate, 10 mM HEPES, 1 mM EGTA, 4 mM MgATP, 0.5 mM Na3 GTP, and 1 mM KCl, adjusted to 260–275 mOsm, pH 7.3). Continuous time series of the membrane potential and spike events were captured in current clamp mode with an AxoPatch MultiClamp 200B (Molecular Devices LLC., San Jose, CA, USA). A part of data were also collected with a MultiClamp 700B amplifier (Molecular Devices). All data were recorded with a Digidata 1440A analog-digital converter (Molecular Devices) controlled by the pCLAMP 10 software using a 10 kHz low-pass filter and a 20 kHz sampling rate. All optogenetics experiments were carried out using a 625 nm LED adjusted to an intensity of ∼4.4 mW/cm² at the fly for CsChrimson activation. The LED was operated using TTL triggers controlled by preset protocols in pCLAMP. Each recording consisted of 10 activations with a stimulation length of 100 ms and a 10 s inter-stimulus-intervals. The TTL triggers were recorded simultaneously with the membrane potential and stored in abf format (Axon Binary File).

In vivo calcium imaging

Calcium imaging experiments were preceded by the same preparation protocol as for whole cell patch-clamp recordings. The calcium concentration of the extracellular saline solution was increased to match the saline used for calcium imaging used before ⁹⁴, with a final concentration of 2.5 mM CaCl2.2H2O All calcium imaging data were acquired with a Prime BSI Express Scientific CMOS camera (Teledyne Photometrics, Tucson, AZ, USA) in an in vivo fluorescence microscope (SliceScope, Scientifica) controlled by the ImageJ based software Micro-Manager 2.0 ⁹⁵. Time series were recorded in a single z-plane at a framerate of 10 Hz for 2 min per trial and stored as tiff stacks. The x and y dimensions of the image frames were adjusted individually to the respective position of the IPCs in each fly. The genetically encoded calcium indicator GCaMP6m ⁵⁹ was expressed in IPCs, while CsChrimson was expressed in the respective MN driver line (see Table 2). GCaMP6m was then excited via a tunable LED (pE-4000 universal fluorescence illumination system, CoolLED Ltd., Andover, UK) set to 470 nm. The light beam was directed through a GFP filter set (Chroma, Bellow Falls, VT, USA), including a 495 nm long-pass filter and an excitation filter with a bandwidth of 450-490 nm. An emission filter with a bandwidth of 500-550 nm blocked unspecific background light before emitted light from the specimen was detected by the sCMOS sensor of the camera. CsChrimson was activated in the same manner as in patch-clamp recordings, but the activation protocols were adapted to the 2 min trials. The LED trigger signal was stored together with a trace of TTL pulses for each frame of the camera as abf files. We used different lengths of activation: 5 x 100 ms, 5 x 1 s, 2 x 5 s, or 1 x 15 s per time series, with the 5 s protocol being used for the majority of the analysis and figures.

Data Analysis

Single-nucleus transcriptomic analysis of IPCs

The pre-processed single-nucleus transcriptome data for IPCs was downloaded from https://www.flycellatlas.org/ ³². This data was further processed and analyzed with the Seurat package (version 4.1.1; https://satijalab.org/seurat/arti-cles/pbmc3k_tutorial.html) in R-Studio ⁹⁶. The original data comprised of 658 single cell transcriptomes, derived from 250 males and 250 females, and sequenced using the Smart-seq2 technology. We used a stringent criterion (Ilp2 > 0 & Ilp3 > 0 & Ilp5 > 0) to retrieve 392 IPC transcriptomes (from 194 male and 198 female flies transcriptomes) which express all three Ilps. Standard commands in Seurat were used to perform the t-distributed Stochastic Neighbor Embedding (t-SNE) analysis and create dot plots and feature plots. Heatmaps were generated using the heatmap package (1.0.12).

Patch-clamp

For patch-clamp recordings, all analysis steps and figure creations were carried out using custom-written MATLAB scripts (Mathworks, Natick, USA). To analyze the recorded patch-clamp traces we first detected all spike events using a custom-written function. To remove high-frequency noise from the traces, a median filter was applied using a time window of 250 ms to each of the 10 time series that compose an optogenetic stimulation experiment. We centered all values around zero by subtracting the median of each trace and we calculated the median of all 10 recordings to estimate the global response of each cell across trials. In order to provide a measure of the baseline activity variability, the 10th and 90th percentiles of all 10 filtered traces were computed in a time window of 900 ms starting 100 ms after recording onset. Thus, we established two thresholds to evaluate the effects of MN populations on the IPC activity. The median of the 10 traces was then used to represent the global response of the recorded IPC, and for the clustering into excited, inhibited, and non-responsive neurons. If the median trace was above the 90th percentile in a time window of 1 s starting at activation onset for longer than 50 consecutive milliseconds, the cell was classified as ‘excited’. Conversely, an IPC was classified as ‘inhibited’ if its activity was below the 10th percentile for 50 consecutive milliseconds, and as ‘unaffected’ if neither of the above cases applied. We ran this analysis across all recorded IPCs and across all tested activation experiments. Thus, we quantified the percentage of IPCs that were excited, inhibited, or unaffected by the activation of each MN population.

To further analyze the spiking activity, we computed the firing rate of each IPC across the 10 trials using non-overlapping time bins of 250 ms. To obtain an overall representation of the population response across all cells, we computed the mean of the median traces of all recorded cells and the mean of their individual firing rates, providing an overview of the collective neuronal activity during the activation experiments. The latency analysis of postsynaptic potentials was conducted using a Savitzky-Golay filter of 4th-order polynomial with four distinct time windows (400 ms, 200 ms, 100 ms, 50 ms) to capture fast dynamics of IPC activity while filtering out individual action potentials. After individual trial traces were filtered, we subtracted the median of each time series to center them to zero. We calculated the median of all 10 recordings and determined the 10th and 90th percentiles of baseline activity, following the previously described procedure for each time window. For cells previously classified as ‘excited’ or ‘inhibited’ we determined the latency of the stimulus as the time delay, following stimulus onset, at which the median of all 10 trials crossed the superior threshold (for excitation) or the inferior threshold (for inhibition). To estimate the overall latency, we averaged the values obtained from the four filter window lengths for each cell. Finally, to compare the magnitude of the IPC modulation via different modulatory populations, we computed the area under the curve (AUC) of the median traces with a trapezoidal numerical integration from activation onset to 1 s after onset. In the same manner, we calculated the baseline AUC from 1 s before stimulus onset to onset and then compared both values, p-values were calculated using the Wilcoxon signed-rank test. To compute a simple effect coefficient, we averaged the AUCs across all IPCs for the respective driver line, which indicated the activation shift of the population activity towards excitation or inhibition. By normalizing them to 1, with 1 being the maximum AUC value across all excited cells of the respective driver line, we were then able to compare the effects between activation experiments.

Calcium imaging

The pre-processing of all calcium imaging data was carried out with custom-written Python code based on the fly2p module as described before ²³. To pre-process the imaging data, each frame of a trial was spatially filtered with a two-pixel-wide Gaussian filter and x/y motion corrected by using phase correlation to a reference image that was generated from the mean over 50 frames in midst of a trial. Fluorescence changes were computed as ΔF/F = (F-F0)/F0 for each pixel over all frames, with F0 being defined as fluorescence baseline averaged over the 10% of lowest-intensity frames in each trial. Next, individual IPC cell bodies were manually delimited as regions of interest (ROIs). From these ROIs, mean ΔF/F values were calculated throughout the time series by averaging pixels within the ROI, and stored as csv files. All further analysis steps were conducted using custom-written MATLAB scripts. First, the ΔF/F traces were aligned and visualized with the respective LED trigger signal. The distribution of the ΔF/F values of all ROIs was calculated in a 5 s window before the first stimulus across all 4 trials for each animal as a baseline. The ΔF/F traces for each ROI from 5 s before stimulus onset to 5 s after stimulus cessation were then superimposed and the mean of these traces was quantified the percent of IPCs per cluster based on the total number of IPCs recorded in every animal. To quantify how the activation of the respective modulatory neurons shifted the population activity of the IPCs, we also calculated an effect coefficient. To this end, we first subtracted the median of a 1 s baseline window before the stimulus onset from the ΔF/F traces and then calculated the area under the curve (AUC) as described for patch-clamp experiments for each ROI, from activation onset to 5 s after activation offset. Finally, we compared the AUC values of the traces without baseline subtraction in a 5 s time window before activation onset to the AUC values in a time window spanning from activation onset to 5 s after activation ceased. These two values were statistically compared, p-values were calculated using a Wilcoxon signed-rank test. calculated. The mean ΔF/F traces for each IPC across all animals tested in the respective driver line were then visualized in a heatmap, normalized to the highest ΔF/F as 1 for each animal. To analyze the effect of each MN population on the IPC activity further, we employed a threshold-based clustering analysis. To this end, we calculated the individual 20th and 80th percentiles of the ΔF/F value distributions in the base-line window for each animal and each IPC separately and used them as upper and lower thresholds. If the ΔF/F trace rose above the upper threshold for a continuous frame number equal 25% of all frames in a window after stimulus onset until 5 s after stimulus offset, this IPC was classified as ‘excited’. If the trace decreased beneath the lower threshold for the same time window, the IPC was classified as ‘inhibited’. If an increase or decrease did not meet the time criteria or if the trace stayed between the thresholds, the cell was classified as ‘unaffected’. Afterwards the mean and standard deviation were calculated for each cluster and each tested driver line. Then, we quantified the percent of IPCs per cluster based on the total number of IPCs recorded in every animal. To quantify how the activation of the respective modulatory neurons shifted the population activity of the IPCs, we also calculated an effect coefficient. To this end, we first subtracted the median of a 1 s baseline window before the stimulus onset from the ΔF/F traces and then calculated the area under the curve (AUC) as described for patch-clamp experiments for each ROI, from activation onset to 5 s after activation offset. Finally, we compared the AUC values of the traces without baseline subtraction in a 5 s time window before activation onset to the AUC values in a time window spanning from activation onset to 5 s after activation ceased. These two values were statistically compared, p-values were calculated using a Wilcoxon signed-rank test.

FlyWire whole-brain connectome

To analyze synaptic connections between IPCs and presynaptic partners, we used the FlyWire whole-brain connectome⁹⁷, which is based on the full adult fly brain (FAFB)¹⁰. All IPCs² and their presynaptic partners (synapse threshold ≥ 5) were identified by using the Connectome Data explorer (Codex)⁶². We then exported the pre- and postsynaptic neuron IDs, as well as predicted transmitters in the presynaptic neurons. For the transmitter prediction score, we used a stringent threshold of 0.62⁶¹, for all predictions below that threshold the neurotransmitter type was labeled as ‘unknown’. The presynaptic inputs were then pooled based on the predicted transmitters and visualized as a connectivity map with the individual IPCs and synaptic weights using a custom-written MATLAB Script.

Finally, all figures were finalized with CorelDRAW (Cascade Parent Limited, Ottawa, Canada).