Belly roll, a GPI-anchored Ly6 protein, regulates Drosophila melanogaster escape behaviors by modulating the excitability of nociceptive peptidergic interneurons
Abstract
Appropriate modulation of escape behaviors in response to potentially damaging stimuli is essential for survival. Although nociceptive circuitry has been studied, it is poorly understood how genetic contexts affect relevant escape responses. Using an unbiased genome-wide association analysis, we identified an Ly6/α-neurotoxin family protein, Belly roll (Bero), which negatively regulates Drosophila nociceptive escape behavior. We show that Bero is expressed in abdominal leucokinin-producing neurons (ABLK neurons) and bero knockdown in ABLK neurons resulted in enhanced escape behavior. Furthermore, we demonstrated that ABLK neurons responded to activation of nociceptors and initiated the behavior. Notably, bero knockdown reduced persistent neuronal activity and increased evoked nociceptive responses in ABLK neurons. Our findings reveal that Bero modulates an escape response by regulating distinct neuronal activities in ABLK neurons.
Editor's evaluation
This convincing and fundamental study uses unbiased genome-wide association analysis to identify a gene, called Bero, encoding Ly6/α-neurotoxin family protein, that affects the way larval Drosophila respond to nociceptive stimuli. This discovery is followed up by the identification of neurons in which Bero function is relevant for the modulation of nociceptive behaviour, namely the abdominal leucokinin-producing neurons. These neurons are activated by nociceptive sensory neurons and can initiate escape behavior. In these neurons Bero modulates both persistent and evoked activities. This elegant work will be of interest to neurobiologists working on genes, neural circuits, and behaviour.
https://doi.org/10.7554/eLife.83856.sa0Introduction
Appropriate escape behavior in response to threatening stimuli is essential for organismal survival (Branco and Redgrave, 2020; Burrell, 2017; Chin and Tracey, 2017; Im and Galko, 2012; Peirs and Seal, 2016). In both vertebrates and invertebrates, the modulation of escape behavior depends on environmental contexts as well as intensity and modality of stimuli, which are sensed and processed in a multilayered neuronal network, subject to genetic variation, that integrates external and internal sensory stimuli (Burnett et al., 2016; Jennings et al., 2014; Mu et al., 2012; Ohyama et al., 2015). Within the network, peptidergic neuromodulation plays a critical role in neuronal information processing (Nässel and Winther, 2010; Taghert and Nitabach, 2012).
When Drosophila melanogaster larvae encounter noxious stimuli, such as harsh mechanical stimulations intrinsic to attacks by parasitoid wasps, they show stereotypic escape behaviors composed of an initial abrupt bending and subsequent corkscrew-like rolling (Hwang et al., 2007; Onodera et al., 2017; Tracey et al., 2003). In larvae, the noxious stimuli are mainly sensed by class IV dendritic arborization neurons (Class IV neurons) underneath the body wall (Tracey et al., 2003). In addition, a type of tracheal dendrite neuron (v’td2 neurons) and class III dendritic arborization neurons (Class III neurons) are known to detect noxious light stimuli and noxious cold stimuli, respectively (Imambocus et al., 2022; Turner et al., 2016). Recent studies have identified several important downstream components of the nociceptive neural circuitry (Ohyama et al., 2015; Takagi et al., 2017; Yoshino et al., 2017). In parallel, several genes and neurotransmitters have been characterized as context-dependent modulators of the behaviors (Burgos et al., 2018; Dason et al., 2020; Hu et al., 2017; Kaneko et al., 2017). While most of these studies identified the nociceptive neural circuitry by connectomic approaches, we used a genetic strategy to provide additional insights into the nociceptive modulation. Importantly, we discovered that D. melanogaster larvae displayed an extensive natural variation in nociceptive escape behavior. This broad behavioral variation allowed us to perform an unbiased genome-wide association (GWA) analysis to explore genetic variants affecting the nociceptive escape behaviors.
In this study, we have demonstrated that an Ly6/α-neurotoxin family protein, Belly roll (Bero), modulates the neuronal excitability of peptidergic interneurons that initiates and facilitates Drosophila nociceptive escape behaviors. We first identified the bero gene through GWA analysis and validated its role in regulating nociception using pan-neuronal bero knockdown and null mutant animals. By searching for bero-expressing neurons, we identified abdominal leucokinin-producing neurons (ABLK neurons) that have been reported to be involved in controlling larval escape behaviors (Hu et al., 2020; Imambocus et al., 2022). We then found that bero knockdown in ABLK neurons resulted in enhanced escape behaviors, similar to bero mutants. Furthermore, we demonstrated that ABLK neurons responded to the activation of nociceptors and initiated nociceptive escape behavior. Notably, bero knockdown in ABLK neurons resulted in reduced persistent activities and increased evoked nociceptive responses in ABLK neurons, implying an essential role of bero in regulating neuronal activities. Based on our results regarding the ABLK neurons, we propose that these neurons may integrate external and internal sensory stimuli, thereby orchestrating nociceptive escape behaviors.
Results
Natural variation in nociceptive responses in wild-type strains of the Drosophila melanogaster Genetic Reference Panel
Upon noxious thermal stimulation, Drosophila larvae show stereotypical nociceptive responses: bending and rolling behavior (Figure 1A; Tracey et al., 2003). Using the Heat Probe Assay (HPA; Tracey et al., 2003), we noticed that there is a large difference in the latency of rolling behavior (rolling latency) between the two commonly used strains, w1118 and Canton-S (Figure 1B). In order to evaluate natural variation in nociceptive responses, we took advantage of the Drosophila melanogaster Genetic Reference Panel (DGRP; Huang et al., 2014; Mackay et al., 2012), a collection of wild-caught Drosophila lines. We quantified the rolling behavior in 38 representative inbred wild-type strains of DGRP and observed a considerable variation in the behavior: the strain-specific median rolling latency ranged from 2.83 s to 10 s and the occurrence rates of rolling behavior (rolling probability) within 10 s ranged from around 30% to almost 100% among these strains (Figure 1C and D). We categorized the rolling probability into three response classes (rolling probability in 2, 5, and 10 s) based on previous studies (Onodera et al., 2017; Terada et al., 2016; Tracey et al., 2003) and found that the three response classes are highly correlated with one another (Figure 1—figure supplement 1A). To identify genetic variants associated with the behavioral variation, we performed multiple trials of GWA mapping analyses using the publicly available web-based analysis tool, DGRP 2.0 (Huang et al., 2014), with four statistical metrics that characterize larval rolling behavior and are mutually dependent to some extent (Figure 1E, Figure 1—source data 1, ‘Materials and methods’). We then identified in total 37 single-nucleotide polymorphisms (SNPs), two deletions, and three insertions in or near 36 candidate genes (Figure 1—source data 2), at a threshold p-value of 1.0 × 10–5, a commonly used threshold p-value in previous studies with DGRP (Figure 1F, Figure 1—figure supplement 1B; Mackay and Huang, 2018).
Identification of CG9336/belly roll (bero) through an RNA interference knockdown screen of candidate genes
To explore impacts of these candidate genes on nociceptive responses, we carried out a secondary functional screen by pan-neuronal RNA interference (RNAi)-mediated gene knockdown. By evaluating 17 RNAi strains that were currently available, we found that pan-neuronal knockdown of the CG9336 gene resulted in significant enhancement of the rolling probability (Figure 2A). To further test whether the phenotype was specifically caused by suppression of CG9336, we generated two distinct short-hairpin RNAs (shRNAs) against CG9336 (shRNA#1 and #2; Figure 2B; ‘Materials and methods’), validated the efficiency of shRNA#2 (Figure 2—figure supplement 1A and B; ‘Materials and methods’), and found that both shRNAs caused a comparable enhancement (Figure 2E, Figure 2—figure supplement 1C). We also performed three sets of effector-control tests for UAS-bero shRNA#2 with different genomic backgrounds. In each group, both the UAS-bero shRNA#2 control animal and the corresponding controls exhibited comparable rolling behavior, which ensures that any observed phenotypes were not due to leaky expression from the UAS transgene (Figure 2E). We then designated the gene as belly roll (bero), after a high-jump style in athletics. bero encodes a member of the Ly6/α-neurotoxin protein superfamily, which is characterized by a three-finger structure, including various membrane-tethered and secreted polypeptides (Loughner et al., 2016). Using several bioinformatics tools, we predicted a signal peptide sequence at the amino terminus, two possible glycosylation sites, a GPI-anchorage site, a putative transmembrane region at the carboxyl terminus, and five putative disulfide bridges, which are typical structural features of the Ly6/α-neurotoxin protein superfamily (Figure 2C, Figure 2—figure supplement 1D; Almagro Armenteros et al., 2019; Baek et al., 2021; Eisenhaber et al., 1999; Gupta and Brunak, 2001; Krogh et al., 2001). A three-dimensional model of Bero was generated by protein structure prediction using AlphaFold2 (Jumper et al., 2021; Figure 2D).
The bero-associated SNP, 2L_20859945_SNP, is located in the second intron of CG9336 isoform A (FBtr0081427) and comprises two nucleotide variations in the DGRP core 40 lines: a minor allele, cytosine, and a major allele, thymine (Figure 2B). We found that most of the 38 DGRP core lines with the bero minor allele (C) showed higher responsiveness than those with the major allele (T) upon noxious thermal stimulations (Figure 2—figure supplement 2A). Based on these observations and our bero knockdown results, we speculated that animals with the bero minor allele (C) might have lower expression levels of bero compared with those possessing the major allele (T). Thus, using reverse transcriptase PCR, we compared the expression levels of bero in the larval central nervous system (CNS) among four DGPR strains that showed extremely quick or slow behavioral responses. We found that animals with the minor allele have lower expression levels of bero (Figure 2—figure supplement 2B). Using online sequence analysis based on JASPAR 2022 (Castro-Mondragon et al., 2022), we predicted that the neighboring nucleotide sequences of the bero-associated SNP might be a target site of certain transcription factors that control the expression level of bero in CNS (Supplementary file 2). Furthermore, we generated a bero null mutant line using CRISPR-mediated mutagenesis, which deleted the entire protein-coding sequence (Figure 2—figure supplement 1E; ‘Materials and methods’), and observed in HPA experiments that both bero heterozygous (beroKO/+) and homozygous (beroKO/beroKO) mutant animals displayed an enhanced rolling escape behavior (Figure 2F). Taken together, these results strongly suggest that bero expression in neurons negatively regulates nociceptive escape behavior.
bero is expressed in several subgroups of peptidergic neurons in the larval CNS
We next sought to find which neurons express bero in the larval CNS. To do this, we examined the overlap between bero-expressing neurons labeled in a protein-trap strain bero-YFP (CPTI-001654; Lowe et al., 2014), in which the endogenous Bero protein is fused with a yellow fluorescent protein variant mVenus, and candidate neurons labeled by a CD4:tdTomato fusion protein driven by distinct GAL4-drivers. We found no bero-YFP expression in the Class IV neurons (nociceptors) or in the downstream neurons that have been reported to be required for the nociceptive escape behavior (Figure 3—figure supplement 1A; Ohyama et al., 2015). Consistent with this, bero knockdown specifically in the corresponding neurons did not alter the rolling behavior (Figure 3—figure supplement 1B). Instead, we observed bero-YFP expression in several groups of peptidergic neurons specified by the bHLH transcription factor dimmed (dimm), including insulin-producing cells (IPCs) and a subset of Eclosion hormone-producing neurons (EH neurons) in the larval brain, and abdominal leucokinin-producing neurons (ABLK neurons) in the lateral region of abdominal ganglions in the larval ventral nerve cord (VNC; Figure 3A and B). In support of our observations, a recent single-cell transcriptomic atlas showed that bero/CG9336 is selectively enriched in peptidergic neurons (Corrales et al., 2022). Bero is also expressed in certain non-neuronal tissues, such as midline glia (Figure 3A). However, we did not study Bero function in the glia in this study.
bero expression in ABLK neurons negatively regulates nociceptive rolling escape behavior
To identify the particular subgroups of neurons that require bero expression for negative regulation of nociceptive escape behavior, we examined rolling escape behavior in larvae with cell type-specific RNAi knockdowns of bero. We found that knockdown of bero using Lk-GAL4, a leucokinin-producing neuron-specific driver, caused enhanced rolling escape behavior (Figure 3E). In contrast, knockdown of bero using either Eh-GAL4, a EH neuron-specific driver, or Ilp2-GAL4, an IPC-specific driver, did not result in any changes in rolling escape behavior (Figure 3E). LK neurons can be classified into four anatomically distinct subgroups in the larval CNS: anterior LK neurons (ALKs), lateral horn LK neurons (LHLKs), subesophageal LK neurons (SELKs), and abdominal LK neurons (ABLKs; de Haro et al., 2010). bero is expressed in all seven pairs of ABLK neurons and is delivered to their neurites, but not expressed in any other LK neurons (Figure 3A and C). Consistent with this, the seven pairs of bero-YFP labeled neurons were immunoreactive for Leucokinin (Figure 3C). High-resolution imaging of individual optical sections showed that Bero localizes on both the plasma membrane and internal membrane of ABLK neurons (Figure 3C, middle). We expressed the dendrite and axon terminal markers, namely DenMark and synaptotagmin-eGFP (syt:eGFP), in ABLK neurons (Figure 3D), demonstrating that Bero distributes along both dendric and axonal regions of the neurons. To investigate the role of bero expression in ABLK neurons in modulating rolling behavior, we tried to induce ABLK-specific bero overexpression using an intersectional flip-out approach (Bohm et al., 2010; Simpson, 2016; Figure 3—figure supplement 1C). We observed that bero overexpression suppressed rolling behavior (Figure 3F), indicating that bero expression in ABLK neurons plays an essential role in the negative regulation of nociceptive behavior. In contrast, in the effector-control test, both the UAS-Bero:FLAG control animal and the corresponding controls showed comparable levels of rolling behavior (Figure 3F). Hence, we focused our analysis on the physiological properties of ABLK neurons and the functional roles of bero in these neurons.
bero promotes persistent fluctuating activities and inhibits evoked nociceptive responses in ABLK neurons
We next addressed whether ABLK neurons respond to stimulation of nociceptors. To test this, we expressed a blue light-gated cation channel ChR2.T159C (a variant of channelrhodopsin-2; Berndt et al., 2011) in nociceptors, and a red fluorescence Ca2+ indicator jRCaMP1b (Dana et al., 2016) in LK neurons. We monitored calcium responses in ABLK neurons in fillet preparations and observed that cell bodies and neurite regions of ABLK neurons displayed persistent fluctuations of intracellular Ca2+ concentration (Figure 4A and B, Figure 4—figure supplement 1A), which was consistent with a previous study (Okusawa et al., 2014). Notably, optogenetic activation of the nociceptors evoked weak but reliable evoked neuronal activities in ABLK neurons, suggesting that these neurons respond to nociceptive inputs (Figure 4C and D; Hu et al., 2020). We used larvae reared on a medium with no additional all trans-retinal (ATR) as controls since optogenetic activation requires extra ATR in the food. A Ca2+ rise was also detected in the control larvae reared on a medium with no additional ATR, but the rising phase preceded the onset of the blue light illumination (Figure 4E). This indicates that the Ca2+ rise could be attributed to the persistent fluctuation in ABLK neurons.
To further investigate the physiological roles of bero in ABLK neurons, we tested whether knockdown of bero alters persistent activities and/or nociceptive responses of the neurons. We found that bero knockdown almost completely abolished persistent fluctuating activities (Figure 4A and B), whereas it increased the amplitude of evoked nociceptive responses (Figure 4C, D and F). Given that bero knockdown enhanced rolling escape behavior in the HPA experiments (Figure 3E), we speculated that the evoked nociceptive responses in ABLK neurons are important for facilitating rolling escape behavior. Given the phenotype of the bero knockdown on spontaneous persistent activities in ABLK neurons, we performed a behavioral analysis of undisturbed free locomotion and found that bero knockdown did not affect larval free locomotion (Figure 4—figure supplement 2A–C).
Remarkably, under ATR (+) condition, persistent fluctuations in ABLK neurons of bero knockdown animals appeared to be enhanced compared to controls (Figure 4D, Figure 4—figure supplement 1B). We postulate that the 561 nm laser used for jRCaMP1b excitation weakly activates ChR2.T159C in nociceptors, inducing weak nociceptive activity. Furthermore, ABLK neurons with bero knockdown displayed increased sensitivity to these nociceptive inputs than controls, which may provide a possible explanation. Therefore, the baseline neural activity under ATR (+) conditions may not represent proper spontaneous activity. Thus, we only compared spontaneous activity under the ATR (–) condition.
Collectively, we have revealed that bero plays a critical role in not only generating a relatively higher level of persistent fluctuating activities, but also inhibiting the evoked nociceptive responses of the ABLK neurons, thereby suppressing the nociceptive escape behavior (Figure 6C).
Proper dynamics of acute evoked activity in ABLK neurons is necessary for the control of rolling escape behavior
To test whether direct manipulations of ABLK activity lead to escape behaviors, we optogenetically activated the neurons using a red light-gated cation channel CsChrimson (Klapoetke et al., 2014). We first activated pan-LK neurons and found that red light illumination elicited a robust rolling behavior that was never observed in control larvae (Figure 5A and B). We then employed the intersectional flip-out approach to selectively express CsChrimson in ABLK or SELK neurons (Figure 5B). Surprisingly, optogenetic activation of ABLK neurons almost exclusively resulted in not rolling but bending behavior, which is considered as a separable initial phase of the entire sequence of escape motions (Figure 5A, Figure 5—figure supplement 1A and B). In contrast, activation of SELK neurons induced robust rolling behavior in all tested larvae.
Based on this observation and bero knockdown results, we hypothesize that acute evoked activity of ABLK neurons may enhance nociceptive escape behavior. To test this, we silenced ABLK neurons by expressing the inward rectifying potassium channel Kir2.1, which mildly enhanced rolling escape behavior (Figure 5D). This result may be attributed to the silencing of persistent activity in ABLK neurons, which may disinhibit acute evoked nociceptive activity. Alternatively, chronic silencing of ABLK neurons may lead to modulation of nociceptive circuitry during development. Thus, in order to manipulate the activity of the neurons transiently, we combined the heat stimulation with optogenetic manipulation. We optically silenced or activated ABLK neurons during the first second of noxious thermal stimulation, such that the heat stimulation and optogenetic manipulation of the neurons were performed almost simultaneously (Figure 5C; ‘Materials and methods’). Inhibition of ABLK neurons using Guillardia theta anion channelrhodopsin-1 (GtACR1; Mohammad et al., 2017) significantly impaired rolling escape behavior (Figure 5E), suggesting that acute evoked activity of ABLK neurons is essential for nociceptive escape behavior. Unexpectedly, however, concurrent activation of ABLK neurons also resulted in significantly decreased rolling escape behavior (Figure 5F). These seemingly inconsistent results suggest that the neuronal activity induced by optogenetic activation is not comparable with the naturally evoked nociceptive responses in ABLK neurons, and proper dynamics of acute activity in ABLK neurons is essential for the facilitation of the rolling escape behavior. The rolling response of animals expressing myr:GFP in ABLK neurons was comparable in both ATR (+) and ATR (–) groups under light illumination, indicating that ATR does not influence the behavior (Figure 5G).
DH44 and octopamine are functional neurotransmitters of ABLK neurons
We finally explored the functional neurotransmitters from ABLK neurons that enhance nociceptive behavior. Previous studies have reported that ABLK neurons produce leucokinin (LK; de Haro et al., 2010) and corticotropin-releasing factor (CRF)-like diuretic hormone (diuretic hormone 44, DH44; Zandawala et al., 2018). Our immunohistochemistry tests additionally identified that ABLK neurons express tyrosine decarboxylase 2 (Tdc2), a key synthetic enzyme for tyramine (TA) and octopamine (OA; Figure 6A). In contrast, ABLK neurons are neither GABAergic, glutamatergic, nor cholinergic (Figure 6A).
To test whether these neurotransmitters are required for the facilitation of nociceptive escape behavior, we used the intersectional flip-out approach to specifically suppress the expression of Lk, Dh44 or Tyramine β hydroxylase (Tbh; a key synthetic enzyme converting tyramine to octopamine) in ABLK neurons by RNAi and examined the rolling escape behavior of these larvae in HPA experiments. We found that suppression of Dh44 or Tbh reduced the rolling probability, indicating functional roles of DH44 and octopamine in the facilitation of nociceptive escape behavior (Figure 6B). However, in the case of Lk suppression, despite the reduction in rolling probability, the corresponding effector-control also showed significantly reduced rolling response, suggesting that the observed phenotype may be attributed to UAS transgene leaky expression (Figure 6—figure supplement 1A).
A previous study showed that knockdown of the ionotropic GABAA receptor Resistance to dieldrin (Rdl) in all LK neurons increased the probability of escape behaviors only at low intensity of nociceptor stimulation (Hu et al., 2020). Consistently, our result showed that ABLK-neuron-specific knockdown of Rdl enhanced rolling escape behavior upon lower thermal stimulations (temperature of the heat probe: 44°C; Figure 6—figure supplement 1B). Another previous study suggested that the 5-HT1B receptor regulates larval turning behavior by changing the periodic activity level of ABLK neurons (Okusawa et al., 2014). We examined the effects of ABLK-neuron-specific knockdown of 5-HT1B on rolling escape behavior in HPA experiments and found that the suppression of 5-HT1B did not show any significant effects (Figure 6—figure supplement 1C).
Discussion
Natural variation in escape behaviors
In this study, we revealed that D. melanogaster larvae display an extensive natural variation in nociceptive escape behaviors. The behavioral variations in a Drosophila population may provide it with adaptive advantages in ever-changing environments, especially when coexisting with geographically and seasonally variable parasitoid wasp species (Bergland et al., 2014; Fleury et al., 2004). In an environment with a high abundance of wasps, it might be adaptive for the larvae to be sensitized to noxious stimuli, whereas in an environment with fewer wasps, a high nociceptive sensitivity might provide only limited protection or even interfere with efficient execution of other essential behaviors such as foraging.
We identified nociception-associated genetic variants through an unbiased GWA analysis (Figure 1—source data 1), which has been employed to discover causative genes corresponding to a given phenotype in a straightforward manner (Long et al., 2014). A strongly associated SNP in the bero locus is located in an intronic region, suggesting that it affects gene expression. Indeed, we discovered a significant correlation between the SNP and the expression level of bero: animals with a homozygous cytosine allele have lower expression levels of bero and are more sensitive to noxious heat stimulations compared to those with a homozygous thymine allele. Consistently, pan-neuronal bero knockdown, bero heterozygous mutant (beroKO/bero+), and homozygous mutant (beroKO/beroKO) animals showed enhanced nociceptive rolling responses, implying that the reduction of the bero expression level has a profound impact on escape behavior. Notably, most homozygous mutant animals with the Canton-S background did not survive to adulthood, suggesting that bero plays an essential role in animal development as well as nociception.
In the validation experiments testing 17 candidate genes, we found that only one pan-neuronal RNAi knockdown produced a significant effect on nociceptive rolling escape behavior (Figure 2A). One interpretation of the low yield is that most of the variants identified by GWA analysis were false positives since the p-values associating these variants exceeded the Bonferroni-corrected 10% significance threshold, 6.95 × 10–8; the alternative is that some polymorphisms in the DGRP strains are gain-of-function mutations. A third possibility is that some candidates are required only in non-neuronal tissues for the behavior. We also noticed that the variants in genes essential for nociception, such as Transient receptor potential cation channel A1 and painless hardly affected the behaviors. One reason is that these genes are essential for survival in natural environments and hence have no functional variants.
How does bero control the activities of ABLK neurons?
Bero has been identified as a member of the Ly6/α-neurotoxin protein superfamily (Hijazi et al., 2009), whereas little is known about its biological functions (Khan et al., 2017). Notably, mammalian Ly6 proteins affect nicotinic acetylcholine receptors (nAChR) to modulate neuronal excitability (Tsetlin, 2015). A Drosophila Ly6-like protein, Quiver (Qvr)/Sleepless (Sss), plays dual roles in sleep regulation, not only suppressing neuronal excitability by enhancing the expression levels and open probability of Shaker potassium channels but also reducing cholinergic synaptic transmission by antagonizing nAChR (Wu et al., 2014; Wu et al., 2010). In our study, the knockdown of bero abolished persistent activity and increased nociceptive responses in ABLK neurons. At least two mechanisms could account for the role of bero in the two distinct neuronal activities. One interpretation could be that bero can regulate the two neuronal activities by modulating the functions of distinct transmembrane proteins separately. Alternatively, bero may target a single protein required for generating the persistent activities, which then constrains the evoked nociceptive responses in ABLK neurons and thereby downregulates the nociceptive rolling behavior (Figure 6—figure supplement 2). In the latter scenario, the bero-targeting protein may elevate intracellular calcium concentration, which interferes with the opening of voltage-gated cation channels, thereby suppressing nociceptive responses (Budde et al., 2002; Oh et al., 2021). Although our current model is based on preliminary observations, it may serve as a conceptual framework for future studies on the roles of bero and ABLK neurons in controlling nociceptive escape behaviors.
In both scenarios, a key issue is the identity of Bero-interacting membrane proteins in ABLK neurons. ABLK neurons express 5-hydroxytryptamine receptor 1B (5-HT1B; Okusawa et al., 2014), leucokinin receptor (Okusawa et al., 2014), insulin-like receptor (Luo et al., 2013), RYamide receptor (Veenstra and Khammassi, 2017), leucine-rich repeat-containing G protein-coupled receptor 4 (Imambocus et al., 2022), and Ecdysis-triggering hormone receptor (Kim et al., 2006). Among them, 5-HT1B has been reported to be required for the maintenance of persistent activity (Okusawa et al., 2014); however, our preliminary results suggest that it is unlikely to be the effectors of Bero (Figure 6—figure supplement 1C). As such, the identity of the relevant Bero-interacting partners still needs to be established.
ABLK neurons initiate the nociceptive escape behaviors
In this study, we identified ABLK neurons as a group of bero+ neurons regulating nociceptive escape behaviors. Activation of all LK neurons was sufficient to induce robust rolling behavior (Figure 5A), which is consistent with previous studies (Hu et al., 2020; Imambocus et al., 2022). Furthermore, we succeeded in the specific activation of two distinct subgroups of LK neurons: ABLK or SELK neurons, both of which can be activated upon nociceptive stimulations (Figure 4D; Hu et al., 2020). A previous study hypothesized that LK neurons decode the decision-related activity, and the activation of LK neurons produces the nociceptive escape behavior (Hu et al., 2020). Here, we showed that the specific activation of ABLK neurons does not result in complete rolling, but rather robust bending behavior, which is considered to be an initial posture of rolling (Burgos et al., 2018; Figure 5A). Based on this observation, we hypothesize that the activation of ABLK neurons induces the initiation of rolling. Consistent with this idea, we found that bero knockdown in LK neurons resulted in elevated levels of nociceptive responses in ABLK neurons (Figure 4D and F) and concomitantly enhanced nociceptive rolling escape behavior (Figure 3E). Notably, activation of SELK neurons leads to the robust rolling behavior (Figure 5A); although this is not the major concern of this study, the pre- and postsynaptic circuitry of SELK neurons should be investigated as an important command circuitry for nociceptive rolling.
We also showed that either activation or inhibition of ABLK neurons during noxious thermal stimulation impairs rolling behavior (Figure 5E and F), suggesting that maintaining ABLK neuronal activity within an appropriate physiological range and/or dynamics is essential for promoting the behavior, supported by previous findings that a proper level of ABLK neuronal activity is crucial for controlling noxious light-evoked turning behavior (Okusawa et al., 2014). More specifically, we propose a scenario that is based on the Ca2+-dependent inhibition of transient excitation in ABLK neurons, as discussed above. In this scenario, transient activation of CsChrimson results in membrane depolarization and a large influx of Ca2+ ions into the cytoplasm of ABLK neurons. This surge of Ca2+ ions can suppress voltage-gated cation channels, thereby inhibiting subsequent nociceptive responses for a substantial period of time (Budde et al., 2002; Oh et al., 2021). Conversely, transient activation of GtACR1 suppresses nociceptive responses in the neurons, resulting in the nociceptive insensitivity. Another scenario proposes that activation of ABLK neurons may not only elicit immediate bending behavior but also exert a delayed suppressive effect on subsequent nociceptive responses to thermal stimulation through modulation of downstream circuits.
Previous research has showed that Kir2.1-mediated silencing of LK neurons increased the noxious touch response (Imambocus et al., 2022). Consistent with these findings, our results showed that Kir2.1-mediated silencing of ABLK neurons mildly enhanced rolling escape behavior in the HPA experiments (Figure 5D), which differs from effect of GtACR1-mediated inhibition (Figure 5E). We propose two possible scenarios. In the first scenario, Kir2.1-mediated hyperpolarization attenuates the persistent activities, thereby facilitating nociceptive responses. An alternative scenario suggests that continuous silencing may modulate the nociceptive circuitry during development, leading to increased nociceptive sensitivity. In contrast, GtACR1-mediated transient inhibition acutely interrupts nociceptive responses without affecting the network excitability.
ABLK neurons as a potential neuromodulatory hub for nociceptive escape behaviors
A previous study reported that ABLK neurons showed evoked neuronal responses upon UV/blue light stimulation, and the responses were required for light-avoidance behavior (Imambocus et al., 2022). Of note, the roles of ABLK neurons in promoting noxious light-avoidance behavior and in initiating nociceptive rolling escape behavior do not necessarily conflict with each other; indeed, combined with nociceptor activation, the blue light stimulation can enhance nociceptive rolling escape behavior (Wietek et al., 2017). These findings imply that ABLK neurons receive converging inputs and affect differential outputs when animals encounter different modalities of noxious stimuli. Upon UV/blue light stimulation, for example, the evoked responses of ABLK neurons are employed for light-avoidance behavior; in contrast, upon noxious heat stimulation, the transient activities with appropriate dynamics facilitate the rolling. It will be interesting to test whether bero knockdown in ABLK neurons increases the UV light-induced activities in the neurons and enhances light-avoidance behavior.
ABLK neurons have known functions in the regulation of water balance and food intake in adults (Liu et al., 2015; Zandawala et al., 2018), suggesting that the persistent activities of the neurons may reflect stresses from desiccation and starvation. Therefore, if persistent activities suppress evoked nociceptive response in the neurons, bero would help integrate stresses with nociceptive inputs in the neurons. In other words, bero expressed in ABLK neurons could modulate nociceptive escape behavior based on the extent of stress in the animal.
Our study revealed that DH44 and octopamine released from ABLK neurons are required for the facilitation of nociceptive rolling behavior. DH44 is an ortholog of mammalian CRF. In mice, CRF+ neurons in the central amygdala mediate conditioned flight (jump escape behavior; Fadok et al., 2017), where CRF release in the central amygdala increases mechanosensitivity and nocifensive reflexes (Bourbia et al., 2010; Ji et al., 2013; Johnson et al., 2015). Moreover, CRF mediates stress-induced thermal hyperalgesia in rats (Itoga et al., 2016), suggesting that the role of DH44 in facilitating nociceptive response might be evolutionarily conserved between insects and mammals.
In summary, our analyses revealed that bero modulates two distinct neuronal activities of a particular subset of bero-expressing neurons, ABLK neurons, which initiate and facilitate the nociceptive escape behavior in D. melanogaster larvae (Figure 6C). Our findings thus raise the possibility that bero has an essential role in monitoring stresses, which may modulate the nociceptive escape behavior (Figure 6—figure supplement 2).
Materials and methods
Resource availability
Lead contact
Request a detailed protocolFurther information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Tadao Usui (usui.tadao.3c@kyoto-u.ac.jp).
Materials availability
Request a detailed protocolLines generated and described in this study are available on request from the lead contact.
Data and code availability
Request a detailed protocolAll data reported in this paper are available from source data tables for each figure.
All original code and scripts are available at this address: https://github.com/TadaoUsui/0_Bero (copy archived at Usui and Li, 2023).
Experimental model and subject details
Fly stocks
Request a detailed protocolD. melanogaster larvae were reared at 25°C and 75–80% humidity with a 12 hr light/dark cycle on standard fly food. Transgenic lines were mainly maintained in either white (w1118) or yellow vermillion (y1, v1) mutant backgrounds unless stated otherwise. No sex-specific effects were part of this study. For fly line details, see Key Resources Table. Lines were obtained from the Bloomington Drosophila Stock Center or the Vienna Drosophila Stock Center unless stated otherwise. For behavioral analysis, wandering third-instar larvae of both sexes were used in this study (120 hr ± 12 hr AEL unless stated otherwise). For calcium imaging, female wandering third-instar larvae were used in this study. Detailed genotypes of experimental and control animals are listed in source data table of each figure.
Method details
Generation of plasmids and transgenic strains
Request a detailed protocol20XUAS-BeroFLAG
Request a detailed protocolBero cDNA was obtained from the Drosophila Genetics Resource Center (DGRC) and amplified from clone FI02856 (DGRC Stock 1621396; https://dgrc.bio.indiana.edu//stock/1621396; RRID:DGRC_1621396) by PCR. FLAG-tagged Bero was generated by inserting the FLAG sequence after the signal peptide sequence at position 26 of the Bero cDNA using overlap-PCR. Primers containing the FLAG-tag sequence were used for amplification and cloning into the pJFRC7-20XUAS-IVS-mCD8::GFP vector via NotI and XhaI restriction sites (the primer sequences are listed in Supplementary file 1). Transgene was made by φC31-mediated genomic integration into the attP2 landing site (WellGenetics Inc, Taipei, Taiwan). pJFRC7-20XUAS-IVS-mCD8::GFP was a gift from Gerald Rubin (Addgene plasmid # 26220; http://n2t.net/addgene:26220; RRID:Addgene_26220).
VALIUM20-berosh#1 and VALIUM20-berosh#2
Request a detailed protocolThe bero RNAi lines were generated as described (Ni et al., 2008). Briefly, the annealed oligo DNA hybrid containing one of the 21-nt sequences targeting the coding region of bero (sh#1: GCACCAAGGACGAGTGCAACG; sh#2: CCTCTATGCCGTTCGTTAAGC) was cloned into a pVALIUM20 vector (the entire oligo DNA sequences are listed in Supplementary file 1). Transgene cassette insertion was achieved by φC31-mediated genomic integration into the attP2 landing site on the 3rd chromosome (WellGenetics Inc).
Generation of bero knockout strains
Request a detailed protocolCRISPR-mediated mutagenesis was performed by WellGenetics Inc using modified methods of Kondo and Ueda, 2013. In brief, the upstream gRNA sequence CAGACTGATCATAACGGCCA[CGG] and downstream gRNA sequence CATCCTGCTCTTCTTCGGCG[TGG] were cloned into U6 promoter plasmids separately. A 3xP3-RFP cassette containing two loxP sites and two homology arms was cloned into pUC57 Kan as the donor template for repair (the primers used are listed in Supplementary file 1). CG9336/bero targeting gRNAs and hs-Cas9 were supplied in DNA plasmids, together with the donor plasmid for microinjection into embryos of control strain w1118. F1 flies carrying selection marker of 3xP3-RFP were further validated by genomic PCR and sequencing (the primers used are listed in Supplementary file 1). This CRISPR procedure generated a 1629 bp deletion allele in the CG9336/bero gene, replacing the entire CDS of the gene with the 3xP3-RFP cassette (Figure 2—figure supplement 1E). This bero KO strain was then outcrossed to Canton-S for 11 generations and established as an isogenic line using 3xP3-RFP as a selection marker.
Immunohistochemistry and confocal imaging
Request a detailed protocolWandering third-instar larvae were dissected in PBS and fixed with 4% formaldehyde in PBS for 20 min at room temperature. After fixation, the larvae were washed five times with PBST (PBS containing 0.3% Triton X-100) and blocked in PBST containing 2% BSA (filtrated with 0.22 μm filter) for 30 min at room temperature. The samples were then incubated with primary antibodies at 4°C for 2 d. After five washes with PBST, the samples were incubated with the corresponding secondary antibodies for 1 hr at room temperature. After further washes, the samples were mounted with ProLong Glass Antifade Mountant (Thermo Fisher, Carlsbad, CA). Images were taken with a Nikon C1Si confocal microscope and processed with Fiji (ImageJ, NIH, Bethesda, MD).
The following primary antibodies were used: chicken polyclonal anti-GFP (1:1000; ab13970, Abcam); rat monoclonal anti-DN-cadherin (1:100; DN-Ex #8, Developmental Studies Hybridoma Bank [DSHB], Iowa City, IA); rabbit polyclonal anti-Lk (1:100; Ohashi and Sakai, 2018); mouse monoclonal anti-FLAG (M2) (1:500; F-3165, Sigma-Aldrich, Burlington, MA); rat monoclonal anti-DYKDDDDK Epitope Tag (L5) (1:500; NBP1-06712, Novus Biologicals, Littleton, CO); mouse monoclonal anti-ChAT (1:50; ChAT4B1, DSHB); rabbit polyclonal anti-GABA (1:100; A2052, Sigma-Aldrich); rabbit polyclonal anti-Tdc2 (1:1000; ab128225, Abcam, Cambridge, UK). Secondary antibodies were the following: Alexa Fluor 488-conjugated AffiniPure donkey polyclonal anti-chicken IgY (IgG) (H+L) (1:500; 703-545-155, Jackson ImmunoResearch Labs, West Grove, PA); Alexa Fluor 405-conjugated goat polyclonal anti-rat IgG (H+L) (1:500; ab175673, Abcam); Alexa Fluor 405-conjugated goat polyclonal anti-rabbit IgG (H+L) (1:500; A31556, Thermo Fisher Scientific); Alexa Fluor 546-conjugated goat polyclonal anti-mouse IgG (H+L) (1:500; A11030, Molecular Probes, Eugene, OR); Alexa Fluor 546-conjugated goat polyclonal anti-rat IgG (H+L) (1:500; A-11081, Molecular Probes); Alexa Fluor 555-conjugated goat polyclonal anti-mouse IgG (H+L) (1:500; A-21424, Thermo Fisher Scientific); Alexa Fluor 546-conjugated goat polyclonal anti-rabbit IgG (H+L) (1:500; A-11035, Thermo Fisher Scientific).
Validation of bero RNAishRNA#2
Request a detailed protocolWe validated the expression level of Bero (Bero-YFP, anti-GFP) in ABLK neurons in control and pan-neuronal bero knockdown larvae (nSyb>bero RNAishRNA#2). bero-expressing neurons were labeled by bero-YFP (CPTI-001654) and immunohistochemistry (chicken anti-GFP, 1:1000; Alexa Fluor 488 anti-chicken, 1:500). The Immunohistochemistry was performed as described above. z stacks were obtained using a Nikon C1Si confocal microscope and processed with Fiji (ImageJ; NIH). Optical setting such as laser power and detector gain were identical for all samples. The cell body in individual ABLK neurons was manually selected as a region of interest (ROI). The normalized fluorescence was defined as FBero-YFP/Fbackground, where FBero-YFP and Fbackground are defined as the mean fluorescence intensity of the ROI and background region, respectively.
Genome-wide association analysis
Request a detailed protocolRolling phenotypes with four distinct statistical metrics (Figure 1—source data 1) covering the 38 representative DGRP lines were submitted to the DGRP2 web-based analysis pipeline (http://dgrp2.gnets.ncsu.edu; Huang et al., 2014; Mackay et al., 2012). Utilizing a linear mixed model, the analysis first adjusts the phenotype for the effects of five major chromosomal inversions (In(2L)t, In(2R)NS, In(3R)P, In(3R)K, and In(3R)Mo) and Wolbachia infection (Mackay et al., 2012). Then the analysis pipeline performed association tests for individual genetic variants with minor allele frequencies of 0.05 or greater, running a linear mixed model using the FaST-LMM algorithm (Lippert et al., 2011). We defined GWA significantly associated genetic variants (Figure 1—source data 2) by a p-value of 1.0 × 10–5, which is a nominal threshold commonly used in previous DGRP studies (Mackay and Huang, 2018). The data were visualized with a custom program written in MATLAB (The MathWorks, Inc, Natick, MA). All genes within 1 kb upstream or downstream of significant associated genetic variants were taken to be associated genes (Figure 1—source data 2).
Gene expression analysis
Request a detailed protocolTo quantitate gene expression of bero in the central nervous system (CNS) across distinct strains, we examined expression levels in the pooled CNS from two wandering third-instar larvae for the following strains: Canton-S (control); four DGRP lines (RAL_208, RAL_315, RAL_391, RAL_705); and the bero KO strain which had been outcrossed to Canton-S for 11 generations. Total RNA from the respective CNS was extracted using Sepasol-RNA I (Nacalai tesque, Kyoto, Japan). Extracted RNA samples were further purified using RNeasy Mini Kit (QIAGEN N.V., Venlo, Netherlands). cDNA was synthesized from purified total RNA using ReverTra Ace qPCR RT Master Mix with gDNA Remover (Toyobo Co., Ltd., Osaka, Japan) following the manufacturer’s instructions. RT-PCR was performed on Mastercycler gradient thermocycler (Eppendorf, Hamburg, Germany) using KOD-Plus-Neo (Toyobo). Expression levels of the target mRNA were normalized against those of αTub84B in the same samples. The primers used are listed in Supplementary file 1. The original files of the full raw unedited gels and figures with the uncropped gels are available in Figure 2—figure supplement 2—source data 2.
Heat Probe Assay (thermal nociception behavioral assay)
Request a detailed protocolThe thermal nociception behavioral assays were performed as previously described (Onodera et al., 2017; Tracey et al., 2003) with slight modifications. Animals were raised at 25°C in an incubator with 12 hr light/dark cycles on standard fly food. Humidity was manually controlled (75–80%). Wandering third-instar larvae of both sexes were picked up from the vial, rinsed twice with deionized water, and transferred to a 140 × 100 mm Petri dish with a bed of fresh 2% agarose gel. Larvae were touched laterally in abdominal segments (A4–A6) with a noxious heat probe (a soldering iron) at 46°C, unless otherwise stated.
The behavioral responses were recorded and later analyzed. The response latency was defined as the time elapsed between the delivery of the heat stimulation and the animals’ initiation of the first 360° rotation. The response latency was computed as 10.05 s if the latency was more than 10 s. Sample sizes were estimated based on previous publications in the field (Jovanic et al., 2016; Ohyama et al., 2015). To reduce possible variation, each genotype was tested multiple times on different days, and the behavioral responses were analyzed in a blinded fashion. The data were visualized with a custom program written in MATLAB (The MathWorks, Inc).
Optogenetic behavioral assays
Request a detailed protocolAnimals for optogenetic manipulation experiments were grown in the dark at 25°C for 5 d on fly food containing 0.5 mM all trans-retinal (R2500; Sigma-Aldrich) unless otherwise stated, and humidity was manually controlled (75–80%).
The optogenetic system consists of a dark chamber with a circuit that controls two arrays of LEDs placed in such a way as to uniformly cover the entire arena. To conduct the experiments, individual wandering third-instar larvae were placed on a 100 × 100 mm matte coated copper plate covered with 10 mL of deionized water. The larva was covered with the water film to be able to crawl on the plate. A 30-s-long pulse of 590 nm orange light (29.3 μW/mm2) was delivered (Amber LUXEON Rebel LED; Quadica Developments Inc, Alberta, Canada). The intensity at the water surface was measured using a HIOKI optical sensor (HIOKI E.E. CORPORATION, Nagano, Japan). Two infrared LED lights (LDR2-90IR2-850; CCS Inc, Kyoto, Japan) surrounding the plate allowed recording of the larval response in darkness. A custom software program written in LabVIEW (National Instruments, Austin, TX) and a multifunction DAQ device (NI USB-6210; National Instruments) were used to control the delivery of LED light pulses and record behavioral responses using a GE60 CCD imager (Library Co., Ltd, Tokyo, Japan). Videos were then examined manually and further analyzed automatically using FIMTrack (Version 2.1; https://github.com/kostasl/FIMTrack; Risse et al., 2017; Lagogiannis, 2016). We then performed the decoding analysis using a support vector machine (SVM)-based custom program written in MATLAB (The MathWorks, Inc) to determine whether each larva displayed rolling behavior in a frame-by-frame manner. Rolling was defined as at least one complete 360° rotation along the body axis. Bending was defined as a C-shape-like twitching without rolling, typically observed as incomplete rolling behavior. Both manually and automatically examined data were visualized with a custom program written in MATLAB (The MathWorks, Inc).
To combine the optogenetic manipulation with the Heat Probe Assay, individual wandering third-instar larvae were placed in a 140 × 100 mm Petri dish as described. The heat stimulation and light illumination were delivered almost simultaneously: a 1-s-long light pulse of 590 nm orange LED (~30 μW/mm2; M590L3; ThorLabs, Newton, NJ) or 565 nm lime LED (~80 μW/mm2; M565L3; ThorLabs) was delivered using a TTL-controlled LED driver (LEDD1B; ThorLabs), respectively, for optogenetic activation or inhibition. The onset of the TTL pulse for the LED driver was synchronized semi-automatically with the heat stimulation using a microcontroller board (Arduino Nano Every; Arduino.cc, Lugano, Switzerland) with a custom-written program (Arduino IDE 2.0; Arduino). The behavioral responses were recorded and then analyzed offline as described.
Locomotion assay
Request a detailed protocolAnimals for locomotion analysis were reared in the dark at 25°C for 4 d on fly food, and humidity was manually controlled (75–80%). For experiments, individual wandering larvae were placed on a 100 × 100 mm matte-coated aluminum plate covered with a bed of 2% agarose gel. Larval locomotion was recorded using a CMOS imager (EG130-B; Shodensha, Osaka, Japan) at a rate of 30 frames per second for 1 min. Larvae that exceeded the imaging area were excluded from the analysis. The locomotor characteristics were automatically analyzed using FIMTrack. The total number of ‘left-bended’ and ‘right-bended’ frames was used to calculate the number of frames for head casting. Average velocity was calculated over 15-frame moving windows.
Bioinformatic analysis of Bero
Request a detailed protocolSignalP-5.0 (https://services.healthtech.dtu.dk/service.php?SignalP-5.0) identified the amino terminal signal peptide (Almagro Armenteros et al., 2019). NetNGlyc-1.0 (https://services.healthtech.dtu.dk/service.php?NetNGlyc-1.0) identified two possible glycosylation sites (Gupta and Brunak, 2001). Big-PI Predictor (https://mendel.imp.ac.at/gpi/gpi_server.html) predicted the GPI-anchor site (Eisenhaber et al., 1999). The putative transmembrane region was predicted using TMHMM-2.0 (https://services.healthtech.dtu.dk/services/TMHMM-2.0/) (Krogh et al., 2001). A three-dimensional model of Bero, along with putative disulfide bridges, was generated using AlphaFold2 and visualized with PyMOL (Jumper et al., 2021).
Calcium imaging in larvae
Request a detailed protocolThird wandering larvae were dissected in calcium-free external saline solution (120 mM NaCl, 3 mM KCl, 4 mM MgCl2, 10 mM NaHCO3, 5 mM TES, 10 mM HEPES, 10 mM glucose, 10 mM sucrose, 10 mM trehalose, 1 mM sodium L-glutamate; osmolality was adjusted to 305 mOsm/kg with water, and the final pH was adjusted to 7.25 with NaOH; modified from Xiang et al., 2010). Buffer was exchanged with oxygenated external saline solution (1.5 mM CaCl2 added to calcium-free external saline solution). The cell bodies were imaged at one frame per second.
470 nm LED blue light (M470L5; ThorLabs) was applied 250 s after the start of imaging and delivered using a TTL-controlled LED driver (LEDD1B; ThorLabs) for optogenetic activation or inhibition. The onset of the TTL pulse for the LED driver was operated using a microcontroller board (Arduino Uno R3; Arduino.cc) with a custom-written program (Arduino IDE 1.0; Arduino.cc). Blue light was applied at 0.37 mW/mm2 for 2.5 s. jRCaMP1b was excited with a 561 nm diode laser (Sapphire 561 LP; Coherent). Images ware acquired at 512 × 512 pixels, with a 250 ms exposure time. The fluorescence signal was captured with images at 1 fps through 610/60 nm bandpass filters (Chroma Technology Corp., Bellows Falls, VT) using an inverted microscope (IX-71; Olympus, Tokyo, Japan) equipped with an EMCCD camera (iXon X3; Andor Technology Ltd., Belfast, UK) with a ×20 objective lens (UCPlanFLN ×20/0.45; Olympus). The data were analyzed with a custom program written in MATLAB (The MathWorks, Inc). The cell body in individual neurons was manually selected as an ROI from the confocal time series.
For the analysis of the persistent activity, the fluorescence change was defined as
where F(t) and B(t) are the fluorescence intensity of the ROI and background region at time point t, respectively; F0(t) and B0(t) are defined as the median of the five minimal values of fluorescence intensity of the ROI and background region, respectively, from t = −10 s to t = 10 s. The sum of ∆FPersistent/F0(t) values was used as the area under the curve (AUC).
For the analysis of the evoked activities, the fluorescence change was defined as
where F(t) and B(t) are the fluorescence intensity of the ROI and background region at time point t, respectively; F0 and B0 are defined as the median of the five minimal values of fluorescence intensity of the ROI and background region, respectively, during a 50 s period just before blue light stimulation was initiated. F0 was used as the baseline fluorescence. The maximum intensity of ∆F/F0(t) from t = 0 s to t = 5 s was used as the ∆Fmax/F0.
Quantification and statistical analysis
Statistics
Sample sizes were similar to those of previous studies. For comparisons of two groups, an unpaired Welch’s t-test was used, or a nonparametric Wilcoxon rank-sum test, as noted in figure legends. For multiple comparison, a Kruskal–Wallis test followed by Dunn’s multiple-comparison test was used. For comparisons of optogenetic activation-induced behaviors, a Chi-squared test was used. Differences were considered significant for p<0.05 (*p<0.05, **p<0.01, ***p<0.001, ****p<0.0001). Statistical testing was performed using MATLAB (The MathWorks, Inc) or R (version 3.6.1, R Development Core Team). Exact statistical tests and p-values for all quantitative data comparisons are listed in the source data tables.
Appendix 1
Data availability
All of the source data are submitted as part of our manuscript.
References
-
Signalp 5.0 improves signal peptide predictions using deep neural networksNature Biotechnology 37:420–423.https://doi.org/10.1038/s41587-019-0036-z
-
The neural basis of escape behavior in vertebratesAnnual Review of Neuroscience 43:417–439.https://doi.org/10.1146/annurev-neuro-100219-122527
-
Calcium-dependent inactivation of neuronal calcium channelsNature Reviews. Neuroscience 3:873–883.https://doi.org/10.1038/nrn959
-
Comparative biology of pain: what Invertebrates can tell us about how nociception worksJournal of Neurophysiology 117:1461–1473.https://doi.org/10.1152/jn.00600.2016
-
JASPAR 2022: the 9th release of the open-access database of transcription factor binding profilesNucleic Acids Research 50:D165–D173.https://doi.org/10.1093/nar/gkab1113
-
Nociceptive circuits: can’t escape detectionCurrent Biology 27:R796–R798.https://doi.org/10.1016/j.cub.2017.07.031
-
Prediction of potential GPI-modification sites in proprotein sequencesJournal of Molecular Biology 292:741–758.https://doi.org/10.1006/jmbi.1999.3069
-
ConferencePrediction of glycosylation across the human proteome and the correlation to protein functionProceedings of the Pacific Symposium. pp. 310–322.https://doi.org/10.1142/9789812799623_0029
-
Nociceptive neurons protect Drosophila larvae from parasitoid waspsCurrent Biology 17:2105–2116.https://doi.org/10.1016/j.cub.2007.11.029
-
Pokes, sunburn, and hot sauce: Drosophila as an emerging model for the biology of nociceptionDevelopmental Dynamics 241:16–26.https://doi.org/10.1002/dvdy.22737
-
Stress-induced hyperalgesiaProgress in Neurobiology 121:1–18.https://doi.org/10.1016/j.pneurobio.2014.06.003
-
Independent optical excitation of distinct neural populationsNature Methods 11:338–346.https://doi.org/10.1038/nmeth0914-972
-
Predicting transmembrane protein topology with a hidden Markov model: application to complete genomesJournal of Molecular Biology 305:567–580.https://doi.org/10.1006/jmbi.2000.4315
-
Fast linear mixed models for genome-wide association studiesNature Methods 8:833–835.https://doi.org/10.1038/nmeth.1681
-
Serotonin and insulin-like peptides modulate leucokinin-producing neurons that affect feeding and water homeostasis in DrosophilaThe Journal of Comparative Neurology 523:1840–1863.https://doi.org/10.1002/cne.23768
-
Charting the genotype-phenotype map: lessons from the Drosophila melanogaster genetic reference panelWiley Interdisciplinary Reviews. Developmental Biology 7:e289.https://doi.org/10.1002/wdev.289
-
Optogenetic inhibition of behavior with anion channelrhodopsinsNature Methods 14:271–274.https://doi.org/10.1038/nmeth.4148
-
Drosophila neuropeptides in regulation of physiology and behaviorProgress in Neurobiology 92:42–104.https://doi.org/10.1016/j.pneurobio.2010.04.010
-
Leucokinin signaling regulates hunger–driven reduction of behavioral responses to noxious heat in DrosophilaBiochemical and Biophysical Research Communications 499:221–226.https://doi.org/10.1016/j.bbrc.2018.03.132
-
Serotonin and downstream leucokinin neurons modulate larval turning behavior in DrosophilaThe Journal of Neuroscience 34:2544–2558.https://doi.org/10.1523/JNEUROSCI.3500-13.2014
-
Fimtrack: an open source tracking and locomotion analysis software for small animalsPLOS Computational Biology 13:e1005530.https://doi.org/10.1371/journal.pcbi.1005530
-
Rationally subdividing the fly nervous system with versatile expression reagentsJournal of Neurogenetics 30:185–194.https://doi.org/10.1080/01677063.2016.1248761
-
SoftwareCode and scripts, version swh:1:rev:0cc10eb6d3607a48c1522a52b4f1217d51d9f3eeSoftware Heritage.
-
Rudimentary expression of Ryamide in Drosophila melanogaster relative to other Drosophila species points to a functional decline of this neuropeptide GeneInsect Biochemistry and Molecular Biology 83:68–79.https://doi.org/10.1016/j.ibmb.2017.03.001
-
Characterization of a set of abdominal neuroendocrine cells that regulate stress physiology using colocalized diuretic peptides in DrosophilaCellular and Molecular Life Sciences 75:1099–1115.https://doi.org/10.1007/s00018-017-2682-y
Article and author information
Author details
Funding
Japan Society for the Promotion of Science (21K06264)
- Tadao Usui
Toray Science Foundation (Toray Science and Technology Grant)
- Tadao Usui
Japan Society for the Promotion of Science (15H02400)
- Tadashi Uemura
Japan Society for the Promotion of Science (JSPS Research Fellowship for Young Scientist)
- Yuma Tsukasa
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank M Futamata, S Oki, and H Imai for excellent technical assistance; members of the Uemura lab for experimental advice and discussions; J Hejna for feedback on the manuscript; T Sakai (Tokyo Metropolitan Univ.) for the anti-Lk primary antibody; J Simpson (UCSB) for tsh-LexA and Scr-LexA strains; A Claridge-Chang (Duke-NUS Medical School) and K Emoto (U Tokyo) for UAS-GtACR1 strain; Drosophila Stocks of Ehime University, Kyoto Drosophila Stock Center, and Bloomington Drosophila Stock Center for providing fly stocks. We also thank FlyBase and Drosophila Genomics Resource Center. This study was supported by KAKENHI from the Japan Society for the Promotion of Science (JSPS) grants 21K06264 to TUs and 15H02400 to TUe, by Toray Science and Technology Grant from Toray Science Foundation to TUs and by a JSPS Research Fellowship for Young Scientists to YT.
Copyright
© 2023, Li et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,353
- views
-
- 112
- downloads
-
- 2
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Genetics and Genomics
Paternal obesity has been implicated in adult-onset metabolic disease in offspring. However, the molecular mechanisms driving these paternal effects and the developmental processes involved remain poorly understood. One underexplored possibility is the role of paternally-induced effects on placenta development and function. To address this, we investigated paternal high-fat diet-induced obesity in relation to sperm histone H3 lysine 4 tri-methylation signatures, the placenta transcriptome and cellular composition. C57BL6/J male mice were fed either a control or high-fat diet for 10 weeks beginning at 6 weeks of age. Males were timed-mated with control-fed C57BL6/J females to generate pregnancies, followed by collection of sperm, and placentas at embryonic day (E)14.5. Chromatin immunoprecipitation targeting histone H3 lysine 4 tri-methylation (H3K4me3) followed by sequencing (ChIP-seq) was performed on sperm to define obesity-associated changes in enrichment. Paternal obesity corresponded with altered sperm H3K4me3 at promoters of genes involved in metabolism and development. Notably, sperm altered H3K4me3 was also localized at placental enhancers. Bulk RNA-sequencing on placentas revealed paternal obesity-associated sex-specific changes in expression of genes involved in hypoxic processes such as angiogenesis, nutrient transport, and imprinted genes, with a subset of deregulated genes showing changes in H3K4me3 in sperm at corresponding promoters. Paternal obesity was also linked to impaired placenta development; specifically, a deconvolution analysis revealed altered trophoblast cell lineage specification. These findings implicate paternal obesity-effects on placenta development and function as one potential developmental route to offspring metabolic disease.
-
- Developmental Biology
- Genetics and Genomics
We present evidence implicating the BAF (BRG1/BRM Associated Factor) chromatin remodeler in meiotic sex chromosome inactivation (MSCI). By immunofluorescence (IF), the putative BAF DNA binding subunit, ARID1A (AT-rich Interaction Domain 1 a), appeared enriched on the male sex chromosomes during diplonema of meiosis I. Germ cells showing a Cre-induced loss of ARID1A arrested in pachynema and failed to repress sex-linked genes, indicating a defective MSCI. Mutant sex chromosomes displayed an abnormal presence of elongating RNA polymerase II coupled with an overall increase in chromatin accessibility detectable by ATAC-seq. We identified a role for ARID1A in promoting the preferential enrichment of the histone variant, H3.3, on the sex chromosomes, a known hallmark of MSCI. Without ARID1A, the sex chromosomes appeared depleted of H3.3 at levels resembling autosomes. Higher resolution analyses by CUT&RUN revealed shifts in sex-linked H3.3 associations from discrete intergenic sites and broader gene-body domains to promoters in response to the loss of ARID1A. Several sex-linked sites displayed ectopic H3.3 occupancy that did not co-localize with DMC1 (DNA meiotic recombinase 1). This observation suggests a requirement for ARID1A in DMC1 localization to the asynapsed sex chromatids. We conclude that ARID1A-directed H3.3 localization influences meiotic sex chromosome gene regulation and DNA repair.