Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT

User menu

Search

  • Advanced search
eNeuro

eNeuro

Advanced Search

 

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Blog
    • Collections
    • Podcast
  • TOPICS
    • Cognition and Behavior
    • Development
    • Disorders of the Nervous System
    • History, Teaching and Public Awareness
    • Integrative Systems
    • Neuronal Excitability
    • Novel Tools and Methods
    • Sensory and Motor Systems
  • ALERTS
  • FOR AUTHORS
  • ABOUT
    • Overview
    • Editorial Board
    • For the Media
    • Privacy Policy
    • Contact Us
    • Feedback
  • SUBMIT
PreviousNext
Research ArticleNew Research, Cognition and Behavior

Local Versus Global Effects of Isoflurane Anesthesia on Visual Processing in the Fly Brain

Dror Cohen, Oressia H. Zalucki, Bruno van Swinderen and Naotsugu Tsuchiya
eNeuro 13 July 2016, 3 (4) ENEURO.0116-16.2016; DOI: https://doi.org/10.1523/ENEURO.0116-16.2016
Dror Cohen
1School of Psychological Sciences, Monash University, Melbourne 3168, Victoria, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Dror Cohen
Oressia H. Zalucki
2Queensland Brain Institute, The University of Queensland, Brisbane 4072, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Oressia H. Zalucki
Bruno van Swinderen
2Queensland Brain Institute, The University of Queensland, Brisbane 4072, Queensland, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Naotsugu Tsuchiya
1School of Psychological Sciences, Monash University, Melbourne 3168, Victoria, Australia
3Monash Institute of Cognitive and Clinical Neuroscience, Monash University, Melbourne 3168, Victoria, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Naotsugu Tsuchiya
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

What characteristics of neural activity distinguish the awake and anesthetized brain? Drugs such as isoflurane abolish behavioral responsiveness in all animals, implying evolutionarily conserved mechanisms. However, it is unclear whether this conservation is reflected at the level of neural activity. Studies in humans have shown that anesthesia is characterized by spatially distinct spectral and coherence signatures that have also been implicated in the global impairment of cortical communication. We questioned whether anesthesia has similar effects on global and local neural processing in one of the smallest brains, that of the fruit fly (Drosophila melanogaster). Using a recently developed multielectrode technique, we recorded local field potentials from different areas of the fly brain simultaneously, while manipulating the concentration of isoflurane. Flickering visual stimuli (‘frequency tags’) allowed us to track evoked responses in the frequency domain and measure the effects of isoflurane throughout the brain. We found that isoflurane reduced power and coherence at the tagging frequency (13 or 17 Hz) in central brain regions. Unexpectedly, isoflurane increased power and coherence at twice the tag frequency (26 or 34 Hz) in the optic lobes of the fly, but only for specific stimulus configurations. By modeling the periodic responses, we show that the increase in power in peripheral areas can be attributed to local neuroanatomy. We further show that the effects on coherence can be explained by impacted signal-to-noise ratios. Together, our results show that general anesthesia has distinct local and global effects on neuronal processing in the fruit fly brain.

  • anesthesia
  • consciousness
  • drosophila
  • frequency tagging
  • isoflurane
  • SSVEP

Significance Statement

Understanding the neural basis of general anesthesia is important for both clinical and consciousness research. Studies in humans show that general anesthesia has distinct local and global effects. Here, we show homologous findings in the fruit fly brains, taking us a step closer to understanding how the loss of consciousness under general anesthesia is evolutionarily conserved across different neuroanatomies. Our unique combination of methods demonstrates that (1) frequency tagging can be used to dissect the neural mechanisms of general anesthesia, (2) anesthesia manipulations deepen our mechanistic understanding of neural processing, and (3) simple modeling can help to clarify unexpected results.

Introduction

Volatile general anesthetics, such as isoflurane, abolish behavioral responsiveness in all animals, but the neural underpinnings of this phenomenon remain unclear (van Swinderen and Kottler, 2014). Although the cellular and molecular mechanisms through which general anesthetics work have been quite well characterized (Franks, 2008; Garcia et al., 2010; Brown et al., 2011), it is unclear what aspect of neural activity is at the core of the profound disconnection from the environment that is induced by all general anesthetics. The difficulty in understanding the mechanisms of general anesthesia may be attributed in part to these drugs targeting multiple processes, from sleep circuits to the synaptic release machinery (van Swinderen and Kottler, 2014).

In humans, general anesthetics have several stereotypical effects on neural activity, as measured by the electroencephalogram (EEG), the best known of which is the increase in delta (0.5-4 Hz) power that is associated with alternation between highly coordinated UP (depolarized) and DOWN (hyperpolarized) states, which is also observed during human non-rapid eye movement (REM) sleep (Murphy et al., 2011; Lewis et al., 2012). Human EEG studies using propofol anesthesia also show an increase in coherent frontal oscillations in the alpha band (8-12 Hz), a potential mechanism for impaired cortical communication (Cimenser et al., 2011; Supp et al., 2011). Human studies combining transcranial magnetic stimulation (TMS) with EEG show that midazolam, propofol, and xenon dramatically disrupt corticocortical communication in response to a TMS pulse (Ferrarelli et al., 2010; Sarasso et al., 2015). These are in agreement with the theoretical suggestion that anesthetics cause the loss of consciousness by interrupting the global integration of cortical activity (Alkire et al., 2008).

The fly model offers a unique opportunity for studying anesthetic action, because it offers the smallest brain (∼100,000 neurons) that is potentially affected by general anesthetics in the same way as the human brain. Isoflurane anesthesia abolishes behavioral responsiveness in fruit flies (Kottler et al., 2013; Zalucki et al., 2015a, Zalucki et al., 2015b), and this is associated with decreased brain activity (van Swinderen, 2006). Genetic manipulations in Drosophila melanogaster are shedding new light on anesthetic action, suggesting that general anesthesia might also involve presynaptic mechanisms as well as the potentiation of sleep circuits (van Swinderen and Kottler, 2014). However, it is currently unclear whether the effects of general anesthesia on neural processing are conserved across all brains, regardless of specific neuroanatomy. To investigate this, we recorded neural activity from multiple regions of the fly brain simultaneously during wakefulness and isoflurane anesthesia, while also measuring brain and behavioral responses to exogenous stimuli.

We used a recently developed multielectrode preparation (Paulk et al., 2013) to record evoked local field potentials (LFPs) across the fly brain in response to flickering visual stimuli. The flickering stimuli produced a periodic response, known as steady-state visually evoked potentials (SSVEPs; Norcia et al., 2015), that allowed us to accurately track the responses in the frequency domain across brain structures, from the optic lobes to the central brain. We hypothesized that isoflurane would globally reduce the power of the SSVEP throughout the brain, but that impaired signal transmission would have a greater effect in the central brain compared with the optic lobes. We found that isoflurane indeed reduced SSVEP power and coherence in central brain areas, but that, surprisingly, responses in the periphery actually increased under isoflurane exposure. We explain these results using a simple model based on known fly neuroanatomy, which suggests a possibility that isoflurane induces an imbalance of the On and Off visual pathways. We further show that the relationship between SSVEP power and coherence can be explained by explicitly considering the relationship between evoked responses and spontaneous brain activity. These results suggest that volatile anesthetics have distinct local and global level effects in all brains, regardless of their specific neuroanatomy, but also that local neuroanatomy is key to understanding anesthetic effects.

Materials and Methods

Animals

Female laboratory-reared D. melanogaster (Canton S wild type) flies (3–7 d post eclosion) were collected under cold anesthesia and positioned for tethering. Flies were dorsally glued to a tungsten rod using dental cement (Synergy D6 FLOW A3.5/B3, Coltène Whaledent), which was cured with blue light. Dental cement was applied to the neck to stabilize the head. The wings of the flies were glued to the tungsten rod to prevent wingbeats or attempted flight during recording. Tethered flies were positioned above a 45.5 mg air-supported Styrofoam ball (Fig. 1a,b), similar to that described by Paulk et al. (2013).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Experimental procedure and paradigm. a, Experimental setup. Flies were dorsally fixed to a tungsten rod and placed on an air-supported ball, where they could freely walk. Flickering stimuli at 13 or 17 Hz were presented through two LED screens to the left and right. Isoflurane in different volumetric concentrations was delivered through a rubber hose. An air puff was used as a startle stimulus to gauge the responsiveness of the flies. A 16-contact electrode probe mounted on an electrode holder was inserted laterally from the left. Only the electrode holder is visible at the depicted scale. b, A close-up view contralateral to the insertion site showing the fly, isoflurane delivery hose, and probe base. c, Example of spontaneous (no presentation of visual stimuli), bipolar rereferenced data before anesthesia (0% isoflurane) from a half-brain probe recording (see Electrode probe insertion). A standardized fly brain is shown for comparison (Paulk et al., 2013, 2015). The electrode contacts are indicated by white dots (not to scale). Channels are grouped as peripheral, estimated to correspond to the optic lobe, and central, estimated to correspond to the central brain.

Electrode probe insertion

Linear silicon probes with 16 electrodes (Neuronexus Technologies) were inserted laterally to the eye of the fly, and perpendicularly to the curvature of the eye. Insertion was performed with the aid of a micromanipulator (Merzhauser), with the electrode recording sites facing posteriorly. For the majority of experiments (14 flies), probes with electrode site separation of 25 μm (3mm-25-177) and 375 μm from base to tip (Fig. 1c) were used. This probe covers approximately half of the brain and henceforth is referred to as the “half-brain probe.” In two additional flies, a probe of 3mm-50-177, with electrode site separation of 50 μm and measuring 703 μm from base to tip was used. This probe covers approximately the whole brain and is referred to as the “whole-brain probe.” Probe tip width (33 μm), base width (123 μm), thickness (15 μm), and electrode site area (177 μm2) are identical for both probes.

A sharpened fine tungsten wire (0.01 inch diameter, A-M Systems) acted as the reference electrode and was placed superficially in the thorax. Recordings were made using a Tucker-Davis Technologies multichannel data acquisition system with a sampling rate of 25 kHz.

The probes were fully inserted until all electrode sites were recording neural activity, which was confirmed by the presentation of visually flickering stimuli (1 and 13 Hz; see Visual stimuli) and observing SSVEPs at the most peripheral electrode site (furthest from the probe tip). The probe was then gently retracted until the most peripheral site showed little to no neural activity. We assumed that this indicated that the most peripheral site was placed just outside the eye. This ensured consistent probe insertion depth among flies.

Visual stimuli

Flickering blue lights (spectral peak at 470 nm with a 30 nm half-peak width) were presented through two LED panels (Fig. 1a). The panels were flickered on and off (square wave, 50% duty cycle) at 13.4 Hz (hereafter cited as 13 Hz) or 16.6 Hz (hereafter cited as 17 Hz), and either to the left or to the right of the fly. There were thus eight possible flicker configurations. These are [off 13], [off 17], [13 off], [17 off], [13 13], [13 17], [17 13], and [17 17], where the number represents the flicker frequency and the location represents the left or right LED panel. An “off” signifies that we turned off the blue LED lights in the respective panel. Visual flickers were presented in sets of 80 trials, consisting of 10 presentation of each flicker configuration. A trial lasted 2.3 s, and the intertrial interval was 0.8 s, taking 248 s to complete the 80 trials. The flicker configuration order was randomly generated with the added restriction that consecutive trials consisted of different flicker configurations. Panel voltage levels were recorded at 25 kHz in the same recording system as the electrophysiological signals. The LED lights were turned off except during the period of visual stimulation (see Fig. 4a).

In three of the flies (two with the whole-brain probe and one with the half-brain probe, all with the graded anesthesia manipulation; see Isoflurane delivery), we also included an additional flicker configuration of [1 1] (1 Hz flicker in both panels). This stimulus was presented once for 20 s before the start of the 80 trials described above. These three flies are used for evaluating the modeled SSVEPs (see Modeling the SSVEPs).

Grouping flicker configurations as ipsilateral or contralateral

The eight flicker configurations were chosen to isolate the effects of flicker frequency (13 vs 17 Hz), flicker interaction (e.g., [13 off] vs [13 17]), and flicker location (e.g., [13 off] vs [off 13]). However, we found that grouping trials as either ipsilateral or contralateral simplified the results and was sufficient to substantiate all our claims (see Results). Under this classification scheme, trials in which a single flicker was presented at the left panel were labeled ipsilateral (see Fig. 3a,b, [13 off] and [17 off]). Trials in which a single flicker was presented at the right panel were labeled contralateral (see Fig. 3a,b, [off 13] and [off 17]). Because the ipsilateral flicker dominated the response, we classified trials in which the same flicker was presented in both panels as ipsilateral (see Fig. 3a,b), and trials in which different flickers were presented in both panels as ipsilateral at the frequency of the ipsilateral panel and contralateral at the frequency of the contralateral panel (see Fig. 3a,b).

Isoflurane delivery

Isoflurane was delivered onto the fly through a rubber hose connected to an evaporator (Mediquip; Fig. 1a,b). The isoflurane was blown onto the fly at a constant flow of 2l/min and continuously vacuumed from the opposite side of the fly. Following the gas chromatography procedure described by Kottler et al. (2013) for measuring isoflurane concentration, we found that the actual concentration near the fly body was 0.3% (vol) when the concentration at the evaporator was set to 1%. Throughout this article, we report isoflurane concentration as the linearly estimated concentration at the fly body, not at the evaporator.

Isoflurane concentrations were manipulated in either a graded or a binary manner over the blocks. In the graded manipulation (N = 3 with the half-brain probe and N = 2 with the whole-brain probe), concentrations were incrementally and sequentially increased over five levels and then reduced to 0%; 0% (air) → 0.06% → 0.18% → 0.3% → 0.45% → 0.6% → 0% (recovery). In the binary manipulation (N = 10 with the half-brain probe), the isoflurane concentration was manipulated over three blocks [0% (air) → 0.6% → 0% (recovery)]. Throughout the article, we distinguish the two periods of 0% isoflurane as 0% (air) and 0% (recovery), before and after drug exposure, respectively. In one fly in which we used the half-brain probe, we administered the graded manipulation up to 0.45%. In a subset of the flies (N = 8 of 14 with the half-brain probe), an additional recovery block [0% (recovery 2)] was performed.

Air puff stimuli and behavioral responsiveness

An olfactory stimulus controller (custom built) was used to deliver six air puffs to gauge the behavioral responsiveness of the flies in each concentration of isoflurane. The inter-air puff duration was ∼1.5 s. Air puffs were delivered before and after the presentation of visual stimuli (see Fig. 4a). Fly movement activity was recorded with a 602f-2 Firewire Camera (Basler) and a 1-6010 12× zoom lens (Navitar) at 30 frames/s, time locked to the onset of the air puff. We used the video data to assess the behavioral responsiveness of the flies under anesthesia (see Movement analysis).

Experimental protocol

After inserting a probe and confirming the visible responses of the flies to an air puff, we initiated our experimental protocol. An experiment consisted of several blocks, each at a different concentration of isoflurane (see Fig. 4a). Each block started with the delivery of a series of air puffs, used to gauge the responses of the fly and to establish the depth of anesthesia (see Movement analysis). Thirty seconds after the startle stimulus, 80 trials of visual flickers were presented (see Visual stimuli). After the completion of 80 trials, the flies were left for an additional 30 s, and then a second episode of air puffs was delivered. After the last air puff was delivered, the isoflurane concentration was immediately changed. Flies were left for 180 s to adjust to the new isoflurane concentration before the next block commenced.

Movement analysis

To confirm the depth of anesthesia, fly movements were analyzed in response to the air puffs. The recorded movies were analyzed to extract the amount of overall movement using custom software written in MATLAB (MathWorks). First, movies were downsampled to 5 frames/s and converted to grayscale. Second, individual images were annotated with the corresponding isoflurane concentration (k = [0 (air), 0.06, 0.18, 0.3, 0.45, 0.6, 0 (recovery), 0 (recovery 2)]%) and saved. Third, images were cropped to include only the body of the fly and were tailored for each fly. Fourth, the mean square error (MSE) between consecutive images across all pixels was calculated, giving one MSE value for each pair of consecutive frames at each isoflurane concentration, k%;Embedded Image where N is the total number of pixels in each image, Embedded Image and Embedded Image represent the grayscale value of pixel (x,y) in frame i and i + 1, respectively, and the sum is taken over all pixels in the image. Finally, the resulting values were averaged over the two episodes of six air puffs in each block with isoflurane concentration k (∼90 frames in total; see Fig. 4a) to obtain Embedded Image . For comparison across flies, we further normalized the values for each fly by dividing the value in k% isoflurane (= Embedded Image ) by the value in 0% isoflurane (air; = Embedded Image ). We refer to the resulting quantity as the movement index (MIk). MI values above and below 1.0 indicate increased and decreased movement compared with 0% isoflurane, respectively. When computing the MI for the recovery period, we used the images from the last experimental block of each fly, which was 0% (recovery) for eight flies and 0% (recovery 2) for five flies.

Local field potential analysis

Electrophysiological data was recorded at 25 kHz and downsampled to 1000 Hz for all subsequent analyses. The most peripheral electrode site was removed from the analysis as it was outside the brain (see Electrode probe insertion). The remaining 15 electrodes sites were bipolar rereferenced by subtracting neighboring electrodes to obtain a set of 14 differential signals, which we refer to as “channels” hereafter (Fig. 1c).

For SSVEP analysis, we segmented the data into 2.3 s epochs according to the flicker configuration and isoflurane concentration. We removed line noise at 50 Hz using the rmlinesmovingwinc.m function from the Chronux toolbox (http://chronux.org/; Mitra and Bokil, 2007) with three tapers, a window size of 0.7 s, and a step size of 0.35 s.

Analyzing power

For each fly, we denote the power of the LFP during visual stimulation at frequency f, in channel i (1–14), flicker configuration l (1–8), and isoflurane concentration k% ([0 (air), 0.06, 0.18, 0.3, 0.45, 0.6, 0 (recovery)]) as Embedded Image , with subscript E meaning “evoked.” Embedded Image is in units of 10log10 (μV2), averaged (in the log scale) over the 10 repetitions of the flicker configuration (see Visual stimuli; see Fig. 4a). Embedded Image was calculated over the 2.3 s trial period using the multitaper method (mtspectrumc.m, http://chronux.org/; Mitra and Bokil, 2007) with three tapers, giving a half bandwidth of ∼0.87 Hz (Mitra and Pesaran, 1999), which is sufficiently fine for our claims in this article. We denote spontaneous power at frequency f in channel i. and isoflurane k% as Embedded Image , with subscript S meaning “spontaneous.” Embedded Image is the power averaged across four 2.3-s-long segments before the start of the visual flicker presentation in units of 10log10 (μV2; see Fig. 4a).

When presenting the results for k = 0% (air), we corrected for baseline levels by subtracting the spontaneous power from SSVEP power, as follows:Embedded Image (1.1)

(subscript B for baseline correction). Embedded Image is reported in decibels, emphasizing that the subtraction is performed after conversion to the log scale.

We use the symbols f1 and f2 to refer to the tag frequency or twice the tag frequency, respectively. The frequencies corresponding to f1 and f2 are flicker configuration dependent (e.g., when the flicker configuration was [13 13], f1 = 13 Hz and f2 = 26 Hz). We refer to power at frequency fn as the average power from −0.5 to +0.5 around the frequency of interest, as follows:Embedded Image (1.2)where N = 4 is the number of frequency bins over which the sum is evaluated. The baseline-corrected SSVEP power Embedded Image was obtained by the substitution of Embedded Image for Embedded Image in Equation 1.2.

When reporting SSVEP power for ipsilateral and contralateral flicker configurations, we separately averaged the flicker configurations for each grouping and the corresponding tags (see Grouping flicker configurations as ipsilateral or contralateral; see Fig. 3a,b). For example, Embedded Image refers to the average power at f1 across the six flicker configurations where the flicker was presented ipsilateral to the probe insertion site (13 Hz for [13 off], [13 13], and [13 17], and 17 Hz for [17 off], [17 17], and [17 13]). The baseline-corrected SSVEP power Embedded Image was obtained by substituting Embedded Image for Embedded Image and repeating the derivation.

The effect of k% isoflurane on SSVEP power is denoted by the symbol Δ and obtained by subtracting respective values in 0% (air) isoflurane, as follows:Embedded Image (1.3)

When reporting the effect on power at the tagged frequency (f = f1 or f2), we averaged the power around the tagged frequency, as in Equation 1.2. To obtain the average for ipsilateral/contralateral flicker configurations (Embedded Image ), we repeated the derivation above with Embedded Image substituted for Embedded Image .

Analyzing coherence

We analyzed coherence between channel pairs using the function coherencyc.m in the Chronux toolbox (Mitra and Bokil, 2007) with five tapers, giving a half bandwidth of 1.40 Hz (Mitra and Pesaran, 1999), which is sufficient for our claims. Our notation and terminology for coherence parallel those used for reporting power, as described below.

SSVEP coherence for channel pair (i,j), flicker configuration l, and isoflurane concentration k%, Embedded Image , is calculated for the SSVEPs over the 2.3 s trials and averaged over the 10 repetitions of the flicker configuration. As spontaneous coherence Embedded Image , we report coherence averaged across four 2.3-s-long segments before the start of visual flicker presentation (see Fig. 4a).

Baseline-corrected SSVEP coherence is used when presenting results in 0% (air) isoflurane and defined as follows:Embedded Image (2.1)

As for power, we refer to coherence at frequency fn (Embedded Image ) as the average coherence from −0.5 to +0.5 Hz around the tagged frequency.

We calculated SSVEP and spontaneous coherence between all channel pairs, resulting in 91 (14*13/2) unique values at every frequency. To summarize these data in a concise way, we grouped channel pairs into periphery (channels 1–6) and center (channels 9–14; Fig. 1c). We report periphery (P), center-periphery (CP), and center (C) coherence as averaged across all pairs of the electrodes within periphery, between center and periphery, and within center, respectively;Embedded Image (2.2)where the superscripts P, CP, and C replace the channel pair superscript. The upper and lower limits of the sums reflect the grouping into peripheral (1–6) and center (9–14) channels, and take into account the fact that coherence is invariant with respect to channel order (C ij = C ji ), while excluding coherence between a channel and itself (C ii = 1). Our results were not sensitive to the exact grouping, such that other schemes, for example, periphery = channels 2–5, center = channels 10–13, gave similar results. We obtained the analogous quantities for baseline-corrected SSVEP coherence at fn, Embedded Image , Embedded Image , and Embedded Image by substitution of Embedded Image for Embedded Image in Equation 2.2.

Paralleling the power analysis, we report SSVEP coherence for ipsilateral and contralateral configurations at f1 and f2 [e.g., Embedded Image refers to the average coherence at f1 across six flicker conditions, where the flicker was presented ipsilateral to the probe insertion site (see Fig. 3a,b)]. The baseline-corrected SSVEP coherence [e.g., Embedded Image ] and the effect of k% isoflurane [e.g.,Embedded Image ] were defined similarly to the analogous quantities for power [e.g., Embedded Image and Embedded Image ].

Modeling the SSVEPs

We modeled the SSVEPs as the sum of two separate linear responses corresponding to the On and Off pathways (see Fig. 5a). The input (depicted as a square wave) is differentiated to extract points of luminance increments and decrements before splitting into two streams corresponding to the On and Off pathways. The responses of the On and Off pathways are summed to give the modeled SSVEP. Mathematically, the output of the model is given by the following:Embedded Image (3.1)where Embedded Image is the output voltage, Embedded Image and Embedded Image are the impulse responses of the On and Off pathways, respectively, and Embedded Image denotes convolution in the time domain. uon and uoff are the half-wave rectified inputs to the On and Off pathways, such that uon(t) = 1 and uoff(t) = 1 signify an increase and a decrease in luminance at time t, respectively. Note that beyond the rectification nonlinearity the model is a linear multiple-input/single-output model (Bendat and Piersol, 2000). Our model assumes that the effect of isoflurane on the SSVEPs can be explained by changes to the impulse responses Embedded Image and Embedded Image alone.

We note that the nonlinearity is cancelled if the impulse response of the On pathway is identical and opposite to the impulse response of the Off pathway (Regan and Regan, 1988). Subbing Embedded Image into Equation 3.1 and using Embedded Image , we obtain the following:Embedded Image (3.2)

Because Embedded Image is simply the (linearly) differentiated input (Embedded Image ), Equation 3.2 shows that the model reduces to a single linear operation of the (linearly differentiated) input when Embedded Image .

The frequency response of the model is given by the following:Embedded Image where the convolution in Equation 3.1 is replaced by multiplication, and capital letters represent the Fourier transforms of their respective variables. The power spectrum of the response of the model is given by the following:Embedded Image (3.3)where we used Embedded Image and Embedded Image for the responses of the On and Off pathways to their respective inputs. The symbol * represents conjugation, and Re() denotes taking the real part.

We note two things about Equation 3.3. First, the response at frequency f is only a function of the responses of the On and Off pathways at frequency f; there is no contribution from other frequencies. In the context of our experiments, this means that the prediction of the model for SSVEP power at f1 and f2 depends only on the stimulus, and the properties of the transfer functions (Embedded Image ) at f1 and f2 .

Second, the prediction of the model for SSVEP power depends on the SSVEP power of the On (Embedded Image and Off pathways (Embedded Image , but also on the cross-spectrum between the responses (Embedded Image ).

The On and Off impulse responses at each channel i and isoflurane concentration k Embedded Image were estimated by averaging the LFP over 20 on–off cycles of the [1 1] flicker configuration, which we presented to three flies (see Visual stimuli; see Fig. 5b). Because the input is a square wave, the half-wave rectification effectively transforms the input into two pulse trains (see Fig. 5a). We then used the estimated impulse responses at each channel and the isoflurane concentration together with Equation 3.1 to predict SSVEPs for the [13 13] and [17 17] flicker configurations by setting the input to a 50% duty cycle square wave with periods (1/13 Hz) and (1/17 Hz), respectively. By computing the Fourier transform of the modeled SSVEPs at each channel i, the two flicker configurations l ([13 13] and [17 17]) and isoflurane concentration k, we obtained the prediction of the model for SSVEP power Embedded Image . The prediction of the model for the effect of k% isoflurane Embedded Image is obtained by substituting Embedded Image for the measured SSVEP power (Embedded Image ) in Equation 1.3. The quantities Embedded Image and Embedded Image were obtained by averaging from −0.5 to +0.5 Hz around the frequency of interest. Note that the model does not have any degrees of freedom for fitting, once the impulse responses are determined by the simple averaging of the data over 20 on–off cycles of the [1 1] flicker configuration. Nothing further is estimated or fit from the data.

Evaluating the SSVEP model

To investigate the relationship between the model and the data, we performed linear regression between the model-predicted and observed effects of isoflurane on SSVEP power, as follows:Embedded Image where b and c were estimated in R (https://www.r-project.org/; R Core Development Team, 2015) using the lm function. We report the Pearson’s correlation coefficient Embedded Image between the prediction of the model and the observed data, and 95% confidence intervals on the slope (b) and intercept (c) obtained by the confint function. Note that a perfect fit between model and data is given by Embedded Image and the line (b = 1, c = 0).

We performed the regression over three flies, all channels (1–14), two flicker configurations ([13 13] and [17 17]), and both f1 and f2, giving 168 paired data points in total. We used the highest concentration of isoflurane [represented as the superscript H in ΔSHE (see Fig. 5g)] presented to each fly, where k = 0.6% for two flies and k = 0.45% for one fly (see Isoflurane delivery).

Signal-to-noise ratio-based estimation of coherence

To investigate the relationship between evoked power and coherence, we assumed a linear framework in which the SSVEPs for each channel pair are related through a linear transfer function in the presence of noise (see Fig. 6b). This framework is conceptually related to the SSVEP model but is completely independent in its implementation and evaluation.

In this framework, the SSVEP at channel i, given by vi (t) passes through the linear transfer function Hij (f) to give the SSVEP at channel j, vj (t). Independent noise enters at each channel separately [ni (t) and nj (t)] to give the recorded responses yi (t) and yj (t). Under these assumptions, squared coherence between channel pairs has an analytical description (for the detailed derivation, see Bendat and Piersol, 2000), as follows:Embedded Image where Embedded Image , and Embedded Image are the power spectrums of ni , nj , vi , and vj , respectively. If we define the signal-to-noise ratios (SNRs) as Embedded Image and Embedded Image , then:Embedded Image (4.1)

Thus, in this simplified setting coherence is totally determined by the SNRs at the respective channels.

To evaluate the SNR-based coherence estimate, we quantified the SNR at each channel and used Equation 4.1 to obtain the prediction of the model of SSVEP coherence. First, we recalculated Embedded Image using the same number of tapers used for the coherence analysis (i.e., five tapers; see Local field potential analysis). We then fitted power law noise to the observed SSVEP power Embedded Image , as follows:Embedded Image (4.2)where we excluded f values from −1.4 to +1.4 Hz around f1 and f2 for each flicker configuration (1.4 Hz corresponds to the half-bandwidth for the coherence measurement, see Local field potential analysis). The purpose of the fit is to estimate the level of neural activity that is not directly tagged by visual flickers. Following the convention in the SSVEP literature (Norcia et al., 2015), we considered the nontagged activity representing the level of noise (n(t) or SN(f)). The parameters Embedded Image and Embedded Image were estimated by linear regression in the log–log scale in the range 1–50 Hz and used to define the noise spectrum.

We define the SNR at frequency f as the observed SSVEP power at f divided by the interpolated noise spectrum at f as follows:Embedded Image (4.3)

An example of the estimated noise spectrums and resulting SNRs for two exemplar channels is shown in Figure 6, c and d. Figure 6e shows the resulting coherence estimate.

The predicted effect of k% isoflurane on coherence for channel pair (i, j) is obtained by subtracting the predicted coherence in 0% (air) isoflurane from the predicted coherence in k% isoflurane, as follows:Embedded Image

For the further analyses, we grouped the electrode pairs into P, CP, and C, and for ipsilateral and contralateral flicker configurations separately as described before (see Fig. 6f).

To provide an overall measure of fit for the SNR-based estimation of coherence, we calculated the MSE between the prediction of the model and the observed effect of 0.6% isoflurane on SSVEP coherence in each fly, across all channel pairs (91), flicker configurations (8), and f1 and f2, as follows:Embedded Image (4.4)where N is the total number of terms in the sum N = 91*8*2 = 1456 (for each of 13 flies).

Separating the contribution of “noise” and “signal” to the SNR-based estimation of coherence

The SNR-based estimation of coherence is completely determined by the SNR (Eq. 4.1), which in turn is a function of the estimated noise levels from nontagged frequency as well as the observed SSVEP power at the tagged frequency (Eq. 4.3). To isolate the contribution of the noise and SSVEP power to the coherence estimate, we defined two additional variants of SNR. In the first, for SNRFN, (where subscript FN is fixed noise), we fixed the noise spectrum to that fitted in 0% (air) isoflurane, as follows:Embedded Image thereby removing the influence of isoflurane on noise levels. In the second, for SNRFE (where subscript FE is fixed SSVEP power), we fixed the SSVEP power to that observed in 0% (air) isoflurane, as follows:Embedded Image thereby removing the influence of isoflurane on SSVEP power. Embedded Image and Embedded Image are used together with Equation 4.1 to obtain two additional estimates of coherence (Embedded Image and Embedded Image ), the effects of k% isoflurane (Embedded Image and Embedded Image ), as well as the grouped coherence over electrode pairs and flicker configurations as described before (Embedded Image , Embedded Image , and Embedded Image , presented in Fig. 6g,h). Paired t tests between the MSEs (Eq. 4.4; obtained for each fly separately), and between the observed effects of isoflurane on coherence and those predicted by each SNR variant were used for assessing statistical significance.

Statistical analysis

We used R (https://www.r-project.org/, R Core Development Team, 2015) and lme4 (Bates et al., 2015) to perform linear mixed-effect analysis of the data. Throughout, the response variable is either power Embedded Image or coherence Embedded Image with four factors corresponding to isoflurane, channel location, harmonic, and the fly.

As to the factors flicker configuration and response frequency, only a subset of all combinations is relevant for our claims. Specifically, when the flicker configuration was [13 13], [13 off], or [off 13] we analyzed power at f1 = 13 Hz and f2 = 26 Hz. When the flicker configuration was [17 17], [17 off], or [off 17] we analyzed power at f1 = 17 Hz and f2 = 34 Hz. When the flicker configuration was [13 17] or [17 13], we analyzed power at f1 = 13 and 17 Hz, and f2 = 26 and 34 Hz. Thus, replacing the response frequency factor, we included the factor flicker location that corresponds to the division into ipsilateral and contralateral flicker configurations (see Fig. 3a,b; see Grouping flicker configurations as ipsilateral and contralateral).

Among those factors, we focused on the crucial isoflurane-specific effects by including interactions between the isoflurane and flicker location factors, isoflurane and harmonic factors, and isoflurane and channel location factors, as well as the triple interaction among isoflurane, channel location, and harmonic factors. In addition, our results in 0% (air) isoflurane imply harmonic-dependent effects for flicker location and channel location, so we included interactions between flicker location and harmonic, and channel location and harmonic. We included random intercepts for all random effects to account for possible offsets between the levels of each factor.

To test for the effect of a given factor or interaction, we performed likelihood ratio tests between the full model described above and a reduced model without the factor or interaction in question (Bates et al., 2015). When applicable, we adjusted p values using the false discovery rate (Yekutieli and Benjamini, 1999).

Results

Evoked responses vary across the fly brain

We presented flickering visual stimuli to awake and anesthetized flies while recording LFPs from different areas of the fly brain (Fig. 1a–c). Before characterizing the effects of anesthesia on the LFPs, we investigated whether different brain areas showed different responses to the visual flicker in 0% isoflurane (air). We first confirmed that in all our recordings (N = 16) we were able to detect the SSVEPs, whose magnitude depended on brain region and the location of the visual flicker (Fig. 2). When we presented 13 Hz flicker ipsilateral (on the same side) to the probe insertion site (Fig. 2a), we saw clear periodic waveforms in the time domain (Fig. 2b) as well as a clear peak at 13 Hz and its multiples (harmonics) in the frequency domain (Fig. 2c), reflecting robust SSVEPs. Figure 2d summarizes the average SSVEP power at 13 Hz (blue) and at its harmonic (26 Hz, red), for each channel. SSVEP power at both 13 Hz (f1) and 26 Hz (f2) was highest around channels 3–6, roughly corresponding to the medulla of the optic lobe (Fig. 1c, cyan structure), and lower responses in channels 8–14, corresponding to the higher-order central structures of the fly brain, as expected and consistent with other SSVEP studies in the fly (Paulk et al., 2013, 2015).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

SSVEP recordings before anesthesia (0% isoflurane). a, Schema of an experiment showing the electrode inserted laterally from the left. The left LED panel is shown flickering at 13 Hz, corresponding to the [13 off] flicker configuration. b, Exemplar mean bipolar rereferenced SSVEP, averaged over 10 trials in the [13 off] condition. The same data from one fly are presented in c and d in different formats. c, Exemplar baseline-corrected SSVEP power spectrum, averaged over the same 10 trials in b (SEB(f); see Local field potential analysis). The blue and the red lines mark the first (f1 = 13 Hz) and second (f2 = 26 Hz) harmonic, respectively. d, Baseline-corrected SSVEP power at f1 (blue) and f2 (red) for the 10 trials of the [13 off] condition (SEB(f1/2)). Note that the narrow shaded areas represent the SD across 10 trials, showing the robust and repeatable nature of the SSVEP paradigm. The grouping into peripheral and central channels is depicted at the bottom. The channels are consistently aligned on the x-axis, b–d.

We observed clear and spatially specific responses for both frequency tags (13 and 17 Hz; main effect of channel locationa, χ2 = 494.0, p < 10−16 See statistical table, row j–l for further detail.). However, we found that a subset of the eight flicker configurations could explain much of the variance in the SSVEP spatial response profiles we observed. Specifically, classifying flicker configurations as either ipsilateral or contralateral to the insertion site (Fig. 3a,b) allowed us to concisely summarize the SSVEP spatial response profiles in the first harmonic (f1 = 13 and f1 = 17; Fig. 3c,d, respectively) and second harmonic (f2 = 26 and f2 = 34; Fig. 3e,f, respectively) across all flicker configurations, into a single figure depicting the four combinations of harmonic and flicker location [(f1,f2) × (ipsilateral, contralateral); Fig. 3g]. SSVEP power for ipsilateral flicker configurations was much higher than contralateral ones (main effect of flicker locationb: χ2 = 211.0, p < 10−16), which is in line with the visual information traveling from the optic lobes to the center of the brain (Fig. 3g; N = 13 flies, 0% isoflurane). For both ipsilateral and contralateral flickers, the responses were stronger at the first harmonic (f1 = 13 or 17 Hz) than at the second harmonic (f2 = 26 or 34 Hz; main effect of harmonic c : χ2 = 718.0, p < 10−16; Fig. 3g). In particular, SSVEP power for contralateral flickers at f2 was weakest and did not show increased responsiveness in peripheral channels (Fig. 3g; confirmed as the significant interaction between flicker location and harmonicd: χ2 = 190.0, p < 10−16). This observation suggests that the second harmonic, f2, reflects more local processing, which is evoked mostly when the flicker is presented to the ipsilateral side. In contrast, the first harmonic, f1, may reflect more global processing, showing a large response even when the flicker is presented to the opposite side of the insertion site.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Baseline-corrected SSVEP power (SEB) and coherence (CEB) before anesthesia (0% isoflurane). a, b, Grouping trials according to ipsilateral and contralateral flicker configurations. A trial consisted of a presentation of one of the following eight flicker configurations: [off 13], [off 17], [13 off], [17 off], [13 13], [13 17], [17 13], and [17 17]. Trials were classified as either ipsilateral or contralateral according to the location of the flicker, with respect to electrode insertion site. a, When analyzing SSVEP power at f1 = 13 Hz or f2 = 26 Hz, trials in which a single flicker was shown ([13 0] and [0 13], row 1) were classified according to the location of the flicker. Trials in which different flickers are shown at each panel ([13 17], [17 13], row 2) are classified according to the location of the 13 Hz flicker. Trials in which both panels show the same flicker ([13 13], row 3) are classified as ipsilateral, as this component dominated the response (see c–f). Trials in which only a 17 Hz flicker is shown ([17 0], [0 17], and [17 17]) are excluded from the grouping. b, Corresponding classification scheme when analyzing SSVEP power at f1 = 17 Hz or f2 = 34 Hz. c–f, Group average (N = 13) baseline-corrected SSVEP power (SEB; see Local field potential analysis) for f1 = 13 (c), f1 = 17 (d), f2 = 26 (e), and f2 = 34 (f) Hz for each of the eight flicker configurations. Grouping flicker configurations as ipsilateral or contralateral accounted for much of the variance, as indicated by the color code (see legend). Error bars represent the SEM across flies (N = 13). g, Group average (N = 13) baseline-corrected SSVEP power at f1 and f2 for ipsilateral and contralateral flicker configurations. Shaded area represents the SEM across flies. h, Group average (N = 13) baseline-corrected SSVEP coherence (CEB) for P, CP, and C channel pairs. Schematics of the fly brain with superimposed examples of channel pairs from each grouping are shown at the bottom. SSVEP coherence followed a similar trend to SSVEP power: higher coherence at f1 than f2 and a decrease toward the center. Contralateral flickers evoked coherence predominantly at f1. Error bars represent SEM across flies (N = 13).

In contrast to these large effects, separating the response profiles for ipsilateral and contralateral flicker configurations into the individual flicker configurations had a minimal but significant effect (main effect of flicker configuratione: χ2 = 17.6, p = 0.015), demonstrating that the grouping into ipsilateral and contralateral flicker configurations accounts for much of the variance.

Evoked coherence varies across the fly brain

Next, we assessed whether SSVEP coherence showed brain region-specific patterns. Coherence measures the strength of linear dependency between two variables in the frequency domain (Bendat and Piersol, 2000), and was used in a similar preparation to investigate closed and open loop behavior in flies (Paulk et al., 2015). Paulk et al. (2015) observed increased SSVEP coherence when flies were engaged in closed-loop behavior, compared to open-loop behavior, where flies were not in control. Notably, SSVEP power alone did not distinguish between the two conditions.

To summarize the coherence data we grouped channels into periphery and center (Fig. 3h), and averaged channel pairs across the periphery, periphery-center, and center. Example pairs from each grouping are shown at the bottom of Figure 3h.

Overall, SSVEP coherence showed a pattern similar to that for SSVEP power. Higher coherence was observed at f1 (=13 or 17 Hz) than at f2 (=26 or 34 Hz), which is clearly seen by comparing blue (f1) and red (f2) bars in Figure 3h (main effect of harmonicf: χ2 = 179.0, p < 10−16), and at ipsilateral than at contralateral flickers, which are seen as lighter bars (ipsilateral) and darker bars (contralateral) in Figure 3h (main effect of flicker locationg: χ2 = 29.5, p < 10−6). The effect of channel location is also strong, with the highest coherence observed between peripheral pairs, followed by peripheral-central pairs, and weakest for central pairs (Fig. 3h; main effect of channel locationh: χ2 = 62.3, p < 10−10). Similarly to SSVEP power, there was an interaction between flicker location and harmonici (χ2 = 6.2, p < 0.02). While contralateral flicker configurations barely evoked coherent SSVEP activity throughout the brain at f2, they evoked location-dependent coherence at f1. Ipsilateral flicker configurations, however, evoked similar location-dependent coherence at both f1 and f2 (Fig. 3h). Together, these similar patterns of results for SSVEP power and SSVEP coherence (Fig. 3g,h) suggest a strong relationship.

The results so far have shown that SSVEP power and coherence have a characteristic spatial response profile, with higher values in the peripheral optic lobe than in the central regions, and that responses at f2 may reflect more local processing, as observed in the limited response to contralateral flicker configurations (Fig. 3g,h). In what follows, we investigated how a volatile general anesthetic, isoflurane, affects these distinct visual responses in the fly brain.

Isoflurane reduces behavioral responses

In flies, like other animals, the behavioral effects of general anesthesia are investigated through behavioral responses to noxious stimuli, such as mechanical vibrations (Kottler et al., 2013; Zalucki et al., 2015a, Zalucki et al., 2015b). In our paradigm, we delivered a series of startling air puffs to the tethered fly exposed to different concentrations of isoflurane (Fig. 4a, blue rectangles). To quantify the responses to the air puffs, we analyzed video recordings of the experiments (see Movement analysis). Before any anesthesia (0% isoflurane), flies responded to the air puffs by moving their legs and abdomen, and this was visible as differences in pixel intensities between consecutive frames of the video recording (Fig. 4b, left column). The 0.6% isoflurane administration rendered flies completely inert, as was evident in the small differences between consecutive frames (Fig. 4b, right column). After the isoflurane concentration was reset to 0%, flies regained pre-anesthesia responsiveness (Fig. 4c): the MI (see Movement analysis) was significantly <1 at 0.6% isofluranej (p < 0.008, paired two-tailed t test; df = 12) and was not different from 1 at the end of the recovery periodk (p = 0.130). MI was significantly lower during 0.6% isoflurane administration than after the recovery periodl (p < 0.003). This analysis confirms that isoflurane abolishes behavioral responsiveness in fruit flies, as demonstrated in previous studies (van Swinderen, 2006; Kottler et al., 2013; Zalucki et al., 2015a, Zalucki et al., 2015b), and also that flies can recover from isoflurane in this preparation.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Isoflurane anesthesia has region- and harmonic-dependent effects on SSVEP power. a, Experimental protocol. An experiment consisted of multiple blocks, each at a different concentration of isoflurane (top black line). Each block proceeded with (1) air puffs (light blue rectangles); (2) 30 s of rest; (3) 80 trials of flicker presentation, corresponding to 10 presentations for each of the eight flicker configurations (gray rectangles); (4) 30 s of rest; (5) air puff; (6) isoflurane concentration change; and (7) 180 s of rest for adjustment to the new isoflurane concentration. b, Isoflurane abolishes behavioral responses. Consecutive video frames (first and second row) in response to an air puff before any anesthesia was administered (0%, left column) and before 0.6% isoflurane exposure (right column). In 0% isoflurane, flies respond to the air puff by moving, which is seen as the large difference in pixel intensity between consecutive frames (left, third row). After 0.6% isoflurane exposure, flies do not respond to the air puff, and there are only small differences between consecutive frames (right, third row). c, Quantifying behavioral responses. Group average (N = 13) movement index (see Movement analysis) was reduced during exposure to 0.6% isoflurane and rebounded after isoflurane levels were reset to 0%. Error bars represent the SEM across flies. d, Isoflurane reduces spontaneous brain activity (ΔSs), measured over four segments of 2.3 s before the start of the presentation of the visual stimuli. Group average (N = 13) effect of 0.6% isoflurane on spontaneous power (ΔS0.6S; see Local field potential analysis). Power is averaged across all channels. The average power for 20–30 and 80–90 Hz is significantly reduced. Error bars represent the SEM across flies. e, Isoflurane reduces SSVEP power (ΔSE) at f1 but increases power at f2 in a concentration-dependent manner. SSVEP power at f1 = 13 Hz (blue) and f2 = 26 Hz (red) for the [13 off] flicker configuration (indicated by the schematic above), at increasing concentrations of isoflurane. For each fly, the SSVEP is first averaged over peripheral channels (triangles, channels 1–6) or central channels (circles, channels 9–14). The channel average is further averaged across flies (N = 3). Error bars reflect the SEM across flies. f, Isoflurane increases SSVEP power at f2 for ipsilateral but not for contralateral flicker configurations. Spatial profile of SSVEP power at f1 (blue) and f2 (red) for contralateral (dark) and ipsilateral (light) flicker configurations in 0.6% isoflurane (ΔS0.6E). SSVEP power is averaged across ipsilateral or contralateral flicker configurations (Fig. 3a,b) first, then is averaged across flies (N = 13). Shaded areas represent the SEM across flies. The SSVEP power at f1 is reduced in central channels for all flicker configurations, indicating an effect on global neural processing. In contrast, SSVEP power at f2 is increased at the periphery, but only for ipsilateral flicker configurations, indicating an effect on local neural processing. The peripheral and central channels over which the average was taken in e are depicted at the bottom of f. ***p < 0.001 and **p < 0.01 in c and d.

View this table:
  • View inline
  • View popup
Table 1:

Statistical table

Isoflurane attenuates spontaneous brain activity

Previous work has shown that the attenuated motor behavior in fruit flies is accompanied by attenuated spontaneous brain activity, quantified as a reduction in mean power in the 20–30 and 80–90 Hz frequency bands (van Swinderen, 2006). We replicated the same effect here using the multiple electrode preparation, by averaging spontaneous power (= Embedded Image ; for the definition, see Local field potential analysis) across each frequency band (f) and across all channels (i) at k = 0.6% isoflurane concentration. Figure 4d shows the group average (N = 13) effect of 0.6% isoflurane on spontaneous power and confirms a significant reduction due to anesthesia in the 20–30m and 80–90 Hzn frequency bands (paired two-tailed t test: df = 12; p < 0.00005 and p < 0.002, respectively).

Isoflurane has opposite effects on SSVEP power at f1 and f2

Previous studies in humans show that general anesthetics have spatially distinct effects on spectral power and coherence (Cimenser et al., 2011). In these studies, the primary sensory areas tend to remain reliably responsive, but higher-order areas show markedly reduce responsivity (Supp et al., 2011; Liu et al., 2012; Mashour, 2013). Isoflurane is known to potentiate GABAergic neurons, resulting in increased inhibition (Alkire et al., 2008; Garcia et al., 2010), and this is consistent with the attenuated spontaneous brain activity observed above and in previous reports (van Swinderen, 2006). Thus, we hypothesized that isoflurane would generally reduce neural responses, but that this reduction would be brain region dependent.

To assess the effects of anesthesia on SSVEP power, we presented visual flickers during exposure to increasing concentrations of isoflurane (Fig. 4a, gray rectangles).

In line with our expectations, we observed a concentration-dependent reduction in SSVEP power at f1. Furthermore, this reduction was more pronounced in central than peripheral areas (Fig. 4e, blue circles vs blue triangles; mean of channels 9–14, N = 3). Surprisingly, responses at f2 increased under anesthesia, but only in peripheral areas (Fig. 4e, red triangles vs circles). We confirmed a strong effect for isoflurane concentration (main effect of isofluraneo: N = 3, χ2 = 631.36, p < 10−16) as well as an interaction between harmonic and isofluranep (N = 3, χ2 = 434.7, p < 10−16).

To better understand the dissociation between the responses at f1 and f2, we collected data from 10 additional flies in which we manipulated isoflurane concentration in a binary manner [0% (air) → 0.6% → 0% (recovery)]. We found that isoflurane reduced SSVEP power at f1 (i.e., Embedded Image < 0) for both ipsilateral flicker configurations (Fig. 4f, light blue) and contralateral flicker configurations (Fig. 4f, dark blue). In contrast, isoflurane increased SSVEP power at f2 (i.e., Embedded Image > 0) for peripheral areas but only in response to ipsilateral flicker configurations (Fig. 4f, light and dark red). This region- and flicker-specific dissociation was confirmed by a strong interaction between isoflurane and channel locationq (χ2 = 187, p < 10−16) and isoflurane and flicker locationr (χ2 = 31.26, p < 0.01), as well as the triple interaction among isoflurane, channel location and harmonics (χ2 = 23.09, p < 0.048).

The reduction in SSVEP power at f1 was observed for both ipsilateral and contralateral conditions (Fig. 4f), suggesting a reduction in global levels of neuronal processing. Further, that this reduction is more pronounced in the central brain is consistent with isoflurane modulating sleep/wake pathways in the central brain (Kottler et al., 2013), in addition to possibly also impairing signal transmission from the periphery to the center. At the same time, the increase in SSVEP power at f2 in the periphery was not observed for contralateral flickers, suggesting that this increase may be attributed to isoflurane acting on some local circuit in the periphery. In the following, we provide a potential explanation with simple, yet biologically plausible, modeling of the SSVEPs.

A minimal model explains the opposing effects of isoflurane on SSVEP power at f1 and f2

A global reduction in SSVEP power at f1 is in line with the reduced neural responsiveness described previously (van Swinderen, 2006) and with global impairment of neural communication across the brain (Alkire et al., 2008). However, the local increase in SSVEP power at f2 in the periphery is not consistent with these. Here, we propose a minimal model that explains these results in a quantitative manner.

First, we considered what type of processing of the input can result in a response at f2. Linear models of SSVEPs from human EEGs have demonstrated a reasonable fit to the observed data (Capilla et al., 2011). However, we can immediately reject purely linear models because our flickering stimuli consisted of a square wave, whose Fourier decomposition consists of only odd harmonics (f1, f3, f5,…). A linear transformation of the input signal cannot result in power at frequencies that are not present at the input in the first place (Norcia et al., 2015). This suggests that a nonlinear process is involved in the generation of the power at f2. The fly visual system is known to contain nonlinear processing (Egelhaaf and Borst, 1989; Reiff et al., 2010; Clark et al., 2011; Behnia et al., 2014)—could a physiologically based, well established nonlinearity account for the unexpected increase in power at f2 that we observed?

A prominent property in visual processing in animals, including fruit flies, is the segregation of the input pathway into luminance increment-responsive (On) and luminance decrement-responsive (Off) pathways (Joesch et al., 2010). The splitting of processing into these two pathways is captured by a nonlinearity in the form of half-wave rectification (Regan and Regan, 1988). Half-wave rectification is also implemented in the fly visual system (Reiff et al., 2010) and represents a biologically plausible, yet simple, nonlinearity.

Figure 5, a and b, summarizes our model, which is based on previous models of nonlinear SSVEP generation (Regan and Regan, 1988). First, the input is linearly differentiated to extract points of luminance change before passing through two opposite half-wave rectifiers, corresponding to segregation into the On and Off pathways. The result is two pulse trains with the same period as the input stimulus and a time delay of half of the stimulus period. The two pulse trains are separately linearly processed by the On and Off pathways, and are finally summed to give the recorded response. We estimated the impulse responses of the On and Off pathways for each channel from the response to a 20 s, 1 Hz flicker that was obtained before the main 13/17 Hz flicker blocks (Fig. 5b; see Modeling the SSVEPs).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

A minimal model explains the unexpected increase in SSVEP power at f2 due to isoflurane. a, Modeling the SSVEPs. The input (depicted as a blue square wave) is (linearly) differentiated to extract points of luminance increments and decrements before splitting into two streams corresponding to the On (pink) and Off (orange) pathways. Each pathway is modeled as a linear operation determined by the impulse response of the pathway. The responses of the two pathways are summed to give the recorded SSVEP. b, The impulse responses of the On and Off pathways were estimated from the response to a 1 Hz flickering stimulus (blue). An example from one channel is shown. No other parameters are fitted from the data. c, Exemplar On (pink) and Off (orange) impulse responses obtained from the 1 Hz flicker presented in both panels [1 1] in 0% isoflurane (air). Note that the negative of the On impulse response (gray) is not identical to the Off impulse response. d, The power spectra of the model output. When the negative of the On impulse response is the same as the Off impulse response, there is no power at f2 (gray). When the empirical On and Off impulse responses are used, the power spectrum has a sharp peak at f2 (black). e, Comparison between the model output (black) and the recorded SSVEP (blue, average across 10 trials) to a [13 0] stimulus in 0% (left) and 0.6% isoflurane (right) in the time domain. An example from one channel is shown. f, Corresponding comparison to e in the frequency domain. Spectra of model output (left) and recorded data (right, averaged across 10 trials of the [13 0] flicker configuration) in 0% (green) and 0.6% (blue) isoflurane show that the model correctly predicts that isoflurane increases SSVEP power at f2. g, The SSVEP model predictions are in excellent agreement with the observed effects of isoflurane. The model correctly predicts the reduction in power at f1 (blue) and the increase in power at f2 (red) for each of three flies (marked by a cross, square, or asterisk), across all channels (14) and both flicker configuration ([13 13] or [17 17]; n = 168, ρ = 0.76). The empirical line of best fit (dashed black) closely resembles the line of perfect fit (solid black).

The model predicts that if the Off pathway impulse response is the exact negative of the On pathway impulse response (Fig. 5c, gray line), there will be no power at the second harmonic (Fig. 5d, gray line). The symmetry between the On and Off responses cancels the nonlinearity (see Eq. 3.2 in Materials and Methods). When the impulse responses for the On and Off pathways are asymmetric, as is expected and consistent with known fly neurophysiology (Behnia et al., 2014), the half-wave rectification is in effect and a prominent peak at f2 is observed (Fig. 5d, black line).

Our minimal model is effective in explaining the opposing effects of anesthesia at f1 and f2 in the time (Fig. 5e) and frequency (Fig. 5f) domain representations of the SSVEP, explaining that isoflurane anesthesia increases the power at f2 by changing the impulse responses of the On and Off pathways. We emphasize that the model is completely determined by the response to the 1 Hz stimulus, which is then used to predict responses for the [13 13] and [17 17] flicker configurations. No parameters are fitted after computing the impulse responses.

To evaluate the model, we computed the correlation coefficient (ρ) and the line of best fit between the model-predicted and observed SSVEP power at f1 and f2 (see Evaluating the SSVEP model). We found excellent agreement between the model prediction and the actual data in 0% (air)t isoflurane (n = 168; ρ2 = 0.9; 95% confidence interval: for slope, 1.0–1.21; for intercept, −6.01 to −2.43) and in the highest concentration of isoflurane delivered to each flyu (n = 168; ρ2 = 0.95; 95% confidence interval: for slope, 1.05–1.15; for intercept, −1.70 to 0.75). Most importantly, the model accurately predicts the effects of isoflurane on SSVEP power. Figure 5g shows the observed effects versus the predicted effects of isoflurane on SSVEP power, and demonstrates that the model captures both the increase at f2 (red) and the decrease at f1 (blue) for each of three flies (marked by cross, asterisk, and square). The predicted and observed effects of isoflurane show a strong linear relationshipv (dashed black line; ρ2 = 0.76; df = 167; 95% confidence interval: for slope, 0.722–0.94; for intercept, 2.50–4.48) that closely resembles a perfect fit (solid black line).

Thus, assuming a minimal, yet biologically plausible, nonlinearity, our model provides a possible explanation for why isoflurane anesthesia unexpectedly increased SSVEP power at f2. Isoflurane is known to potentiate GABAergic neurons (Alkire et al., 2008; Garcia et al., 2010) so a “global” reduction in neural responses is expected. Consistent with this, we observed that SSVEP power at f1 was globally reduced, while SSVEP power at f2 increased, but only locally at the periphery. Our model, however, predicts that if isoflurane caused an imbalance between the On and Off pathways, the nonlinearity can be enhanced, causing increased SSVEP power at f2. The imbalance between the On and Off pathways is likely to emerge only when the local circuit is strongly engaged. Such strong engagement is much more likely for ipsilateral than for contralateral flicker configurations, which is consistent with the observation that the increase at f2 was observed only for ipsilateral flicker configurations. While our model cannot pinpoint the cellular/molecular mechanisms underlying this change, one potential cause is a widespread impairment in synaptic efficacy, independent of sleep circuits, that here results in affected On and Off responses (van Swinderen and Kottler, 2014; Zalucki et al., 2015).

Isoflurane has opposite effects on SSVEP coherence at f1 and f2

We next investigated how isoflurane anesthesia affected the observed SSVEP coherence. Generally, the effects were closely related to the changes observed for SSVEP power. Following the same procedure for 0% isoflurane (Fig. 3h), we summarized the results by averaging coherence among pairs of recording sites within periphery, between periphery and center, and within center (see Analyzing SSVEP coherence). As expected, 0.6% isoflurane significantly modulated SSVEP coherencew (main effect of isoflurane, χ2 = 185, p < 10−16).

The effects of 0.6% isoflurane on SSVEP coherence (ΔCE , N = 13 flies) are qualitatively similar to those on power, in terms of channel pair location, flicker location, and harmonic, as shown in Figure 6a. Isoflurane reduced coherence at f1 but increased coherence at f2 (Fig. 6a, red vs blue bars; interaction between harmonic and isofluranex: χ2 = 146, p < 10−16). The reduction at f1 was greater in the center (Fig. 6a, right column), while the increase at f2 was predominantly observed at the periphery (Fig. 6a, left column; interaction between isoflurane and channel locationy: χ2 = 29.25, p < 0.00001). The reduction at f1 was observed for all flicker configurations (Fig. 6a, light and dark blue), but the increase at f2 was only observed for ipsilateral flicker configurationsz (Fig. 6a, light red vs dark red; interaction between isoflurane and flicker location: χ2 = 12.35, p < 0.009). The triple interaction among isoflurane, harmonic, and channel was not significantaa (χ2 = 0.66, p = 0.72). The results imply that in our paradigm there is a strong connection between SSVEP power and SSVEP coherence. In the following section, we dissect this by assuming a linear framework that provides an estimate of coherence based on the signal-to-noise ratios of tagged power in the frequency domain.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

A minimal model based on the SNRs of SSVEP power explains the effects of isoflurane on SSVEP coherence. a, Group average (N = 13) of the observed effects of anesthesia on SSVEP coherence (ΔCE). Isoflurane decreases SSVEP coherence at f1 (blue) and increases coherence at f2 (red). Isoflurane decreases coherence at f1 for all flicker configurations throughout the brain. Isoflurane increases SSVEP coherence at f2 at the periphery, but only for ipsilateral flicker configurations (light red). Schematics of the fly brain with superimposed examples of channel pairs from each grouping are shown at the bottom. Error bars represent the SEM across flies. b, Linear framework for SNR-based coherence estimates. The SSVEP in channel vi (t) is related to the SSVEP in channel vj (t) through the transfer function H(f). Independent noise ni (t) and nj (t) enters at each channel separately to give the recorded SSVEPs yi (t) and yj (t). Under this scheme, SSVEP coherence has an analytic expression based on the SNR at each channel, given by Equation 4.1 (see SNR-based estimation of coherence). c–e, Estimation of SNR and coherence prediction from the data. c, Noise levels were estimated from nontagged frequencies at each channel, isoflurane concentration, and flicker configurations by fitting power-law noise to the SSVEP spectrum (see SNR-based estimation of coherence). Exemplar average SSVEP spectra and power-law fits for the [13 off] flicker configuration in 0% isoflurane for two channels, indexed by i and j, are shown. A schematic of the fly brain and channel locations is shown at the top. d, The SNR of channels i and j, obtained by dividing the spectrum of the SSVEP by the power-law fit in the linear scale. e, Example of SSVEP coherence prediction for channel i and j in d based on the SNR through Equation 4.1. f, The SNR-based model correctly predicts the effects of isoflurane on coherence (Embedded Image ). Group average (N = 13) SNR-based prediction of the effects of 0.6% isoflurane. The format and color scheme is the same as in a. g, h, Coherence predictions using different definitions of the SNR. g, Prediction based on SNRs using noise levels in 0% isoflurane (Embedded Image ; see SNR-based estimation of coherence). h, Prediction based on SNR using SSVEP power levels in 0% isoflurane (Embedded Image ). i, Quality of coherence prediction from each model. MSE between the observed (a) and each of the three predictions (f–h) averaged across all flies, channels, flicker configurations, and f1 and f2. This demonstrates that the effects of isoflurane on SSVEP coherence are largely attributed to the effects of isoflurane on SSVEP power, not on noise. Error bars represent the SEM across flies (N = 13). ***p < 0.001 and *p < 0.05.

A minimal model explains the opposing effects of isoflurane on SSVEP coherence at f1 and f2

The observation that isoflurane affected coherence and power in a similar way suggests that in our data the two measures are linked. We explain the opposing effects of isoflurane on SSVEP coherence at f1 and f2 with another simple model. The model assumes that in each pair of channels, one channel (Fig. 6b, vi (t)) receives an input from the initial sensory processing (i.e., Fig. 6b, On/Off response box) and the other (vj (t)) is a linearly filtered version of the first (represented by the transfer function H(f)). Finally, independent noise (ni (t) and nj (t)) enters at each channel, giving the two output voltages (yi (t) and yj (t)). This simple framework allows us to apply an analytic derivation of coherence based on the SNRs at each channel (see Signal-to-noise ratio-based estimation of coherence; Bendat and Piersol, 2000).

The model involves two main assumptions about the nature of the LFP. First, the relationship between the SSVEPs in each pair of channels is linear. Second, spontaneous activity is independent of evoked activity. There is evidence that purely linear processing can provide a good description of the neural responses in the Drosophila brain (Bialek et al., 1991; Behnia et al., 2014), and independence between evoked and spontaneous neural activity is often assumed when assessing neural connectivity in human LFP studies (Truccolo et al., 2002; Wang et al., 2008). Thus, both assumptions are physiologically plausible. However, given the present lack of understanding regarding the physiological underpinning and inter-area properties of the LFP, we cannot propose a mechanistic, physiology-based perspective for this model.

To quantify the SNRs, we first estimated the noise level by fitting power-law noise to the power spectrum at the nontagged frequencies during visual stimulation for each channel, flicker configuration, and isoflurane concentration (Fig. 6c; see SNR-based estimation of coherence). Note that the SSVEP paradigm allows us to operationally regard power at the tagged frequency (f1 and f2) as signal and power at non-tagged frequencies as noise (Norcia et al., 2015). Dividing the measured SSVEP power by the estimated noise levels (in the linear scale) provides our estimation of the SNR (Fig. 6d; see SNR-based estimation of coherence). The SNR estimates together with Equation 4.1 provide a coherence estimate (Fig. 6e). Finally, we separately obtained estimates of the SSVEP coherence in 0% and 0.6% isoflurane concentrations to predict the effects of isoflurane on coherence (Embedded Image ; Fig. 6f–h).

The predicted effects of isoflurane on SSVEP coherence is in excellent agreement with the observed data (Fig. 6a, observed coherence, f, for the model prediction). The model captures the general decrease of coherence at f1 as well as the increase of coherence in the periphery at f2 for ipsilateral flicker configurations.

In this framework, the effects of isoflurane on noise level as well as on SSVEP power both contribute to the SNR-based prediction of coherence. But what is the relative contribution of nontagged noise and tagged signal to our successful prediction of SSVEP coherence? To isolate the relative contribution, we recalculated the SNR by fixing either noise or signal to 0% isoflurane levels, which we call SNRFN and SNRFE (see Separating the contribution of “noise” and “signal” to the SNR-based estimation of coherence). The results (Fig. 6g,h), clearly show that the contribution of the signal (or evoked response) is much more important for the model prediction.

We formally confirmed the above observation by computing the mean squared error between each SNR-based prediction (Embedded Image , Embedded Image , and Embedded Image ) and the observed data (Embedded Image ) across all channel pairs and flicker configurations (Fig. 6i; see Separating the contribution of noise and signal to the SNR-based estimation of coherence). Disregarding the effect of isoflurane on power (Embedded Image ) resulted in considerably worse predictions thanab Embedded Image (p < 0.0001, df = 12) and ac Embedded Image (p < 0.0001, df = 12). However, disregarding the effects of isoflurane on noise resulted in only slightly (but significantly) worse predictionsad (Embedded Image vs Embedded Image : p < 0.040, df = 12). This means that the observed effects of isoflurane on SSVEP coherence, which is a global decrease of coherence at f1 and a local (peripheral) increase of coherence at f2, is largely attributed to the effect of isoflurane on SSVEP power at the tagging frequency of the stimulus, rather than to general effects on nontagged frequencies.

Discussion

In this study, we showed that isoflurane has distinct local and global effects on the fruit fly brain. This was made possible by our approach that combines pharmacological manipulation of the states of the brain through anesthetics, perturbation of the neural circuits through periodic visual stimuli, and analysis and modeling of behavior and neural data. Together, these components synergistically provide a fuller picture of the effects of isoflurane anesthesia on visual processing, which may generalize to the brains of animals other than flies.

As to the mechanisms of anesthesia, recent studies (Alkire et al., 2008; Mashour, 2013) suggest that reduced cortical communication is at the core of the anesthetic-induced loss of consciousness. In particular, increased synchronous activity induced by anesthesia has been suggested to adversely interfere with the communication between brain areas (Supp et al., 2011; Lewis et al., 2012; Sarasso et al., 2014), which may explain the failure of the propagation of evoked responses from primary to higher-order areas (Supp et al., 2011; Liu et al., 2012; Mashour, 2013). The volatile general anesthetic isoflurane also abolishes behavioral responses in fruit flies at concentrations similar to those required for human anesthesia (van Swinderen, 2006; Kottler et al., 2013; Zalucki et al., 2015a, Zalucki et al., 2015b; present study), suggesting that the neural mechanisms through which this anesthetic works may be conserved in most animals. Here we questioned whether isoflurane has distinct effects on local and global processing in the fruit fly brain, and thereby investigate whether an entirely different brain neuroanatomy might reflect similar fundamental effects on neural processing under general anesthesia.

Using a multielectrode preparation allowed us to record from different brain areas simultaneously, and to assess brain region-dependent effects. By presenting flickering visual stimuli, we could isolate the neural response in the frequency domain. The frequency decomposition revealed specific effects of anesthesia on the first harmonic (f1, 13 or 17 Hz) and second harmonic (f2, 26 or 34 Hz), which reflected global and local visual processing. Our results show that the reduction in behavioral responses is accompanied by attenuated spontaneous brain activity, and this was also true for the SSVEPs in the central brain, which were reduced for all stimulus configurations, indicating an effect on global neuronal processing at f1. In contrast, and to our surprise, local responses at f2 in peripheral areas increased, but only for ipsilateral flicker configurations. Modeling the SSVEPs was crucial to understand this unexpected effect, explaining that the f2 power increase in the periphery can be attributed to isoflurane-induced changes of the On and Off response pathways in the optic lobes. We further showed that the analogous effects of isoflurane on coherence can be explained by explicitly considering how isoflurane affects the tagged brain activity (both f1 and f2). Overall, the reduction in SSVEP power and coherence in the central brain fits with the view that general anesthetics target interarea neural communication, impairing the transmission of the visually evoked responses from the optic lobes to central brain structures.

One possibility to describe our finding is to separate the effects of isoflurane into the central and the peripheral, rather than the local and global, as we presented the effects in the article. We think that using the term “central effect” is not as well suited because it fails to capture the fact that the effect of isoflurane on SSVEP power and coherence at f1 generalized over flicker configurations (both ipsilateral and contralateral; Fig. 4f). In addition, we argue that this effect involves the transmission of the signal between the periphery and the center. As such, we think that the response at f1 reflects global processing, rather than central processing, and f2 reflects local (which can be also considered as peripheral) processing.

Evoked and spontaneous activity

The characterization of evoked responses, as opposed to spontaneous activity, through the delivery of a controlled input can reveal additional information about the system. In our experiment, the SSVEPs increased in peripheral areas at f2 for specific flicker configurations, revealing a clear difference between the effects of isoflurane on the periphery and center of the fly brain. The use of evoked activity in studying general anesthesia may be particularly important because it allows tracking a stimulus-related neural process across the brain, potentially making it easier to identify impaired interarea communication. In SSVEP paradigms, the signal is operationally defined as activity at the tag and its harmonics, and this assumption makes it straightforward to define SNRs of evoked activity. This is more difficult with spontaneous activity, where signal and noise cannot be easily separated. Our operational definition of signal and noise, following the tradition of SSVEP studies (Norcia et al., 2015), allowed us to explicitly consider how SSVEPs at the tag frequency combine with nontagged activity (through the quantification of the SNR) to influence coherence. In our data, the effects of isoflurane on SSVEP coherence could be largely attributed to the effects of isoflurane on SSVEP power, as opposed to effects on surrounding, non-stimulus-related activity (Fig. 6f–i).

Focusing on neural activity at predefined frequencies, however, is also a limitation of SSVEP paradigms as this only probes the behavior of the system in a narrow range: the tag and its harmonics. This is particularly important in the context of nonlinear systems whose frequency response can be highly input dependent. We expect that both our modeling of the SSVEPs and the SNR-based estimation of coherence will need to be expanded when the system is evaluated over a broader dynamic range.

Neural substrate of the SSVEPs

Our modeling of the SSVEPs concisely yet plausibly accounts for the unexpected increase in power at f2 observed in the periphery. Given the vast literature on elementary motion detection circuitry in flies (Egelhaaf and Borst, 1989; Reisenman et al., 2003; Borst and Euler, 2011), it may be possible to provide more comprehensive modeling. However, for the purpose of explaining the unexpected effects of anesthesia, our minimal modeling was sufficient and provided a physiologically plausible explanation: isoflurane most likely affected the responses of local On and Off pathways, which, combined with the presentation of a periodic stimulus, resulted in increased power at f2. We note that, in principle, the model can be used to predict the response to arbitrary flicker configuration, but in this work we focused on the observation that required further explanation: the unexpected increase in power at f2. We are currently investigating whether isoflurane-resistant mutants, which have recently been identified (Kottler et al., 2013; Zalucki et al., 2015), can be used to further clarify the mechanisms involved in the global decrease and local increase in responsiveness that we observed.

Even for our simple model, it is not straightforward to assign a fine neural substrate to the SSVEP because there are many connections between the fly optic lobes, such that stimulation in one lobe causes activation in the other (Haag and Borst, 2008). While our recordings (and those by Paulk et al., 2015) clearly show that SSVEP power is much smaller when the flicker is presented to the opposite eye, the broad-field flicker prevents us from precisely disentangling the relative contribution of each optic lobe to the LFP. Another factor is the aggregate nature of the LFP; while the first On and Off responsive cells may be observed as early as the lamina (Reiff et al., 2010), we cannot tell how much these cells contribute to the LFP, compared to other downstream neurons. Future studies separating the contributions of the On or Off pathways to the LFP via genetic manipulations and the use of stimuli that target each pathway separately will help to clarify the neural substrate of the SSVEP.

Slow-wave and interarea neural communication

Sleep and general anesthesia are defined by similar criteria, and there is evidence for some shared mechanisms (Franks, 2008). The involvement of sleep mechanisms in the impairment of cortical communication observed in general anesthesia (Ferrarelli et al., 2010; Sarasso et al., 2015) is not established, but one possibility is that the stereotypical DOWN states that manifest as the human EEG slow wave and are observed in both general anesthesia and non-REM sleep states, may prevent long-range coordinated activation (Sarasso et al., 2014).

Recent findings extend the proposed relationship between sleep and anesthesia to fruit flies, in which genetic manipulations of sleep circuits can confer both resistance and hypersensitivity to isoflurane (Kottler et al., 2013). However, to date, there has been no evidence of a slow wave in sleep or anesthesia in flies (van Swinderen, 2006; Kirszenblat and van Swinderen, 2015), and we found no evidence of it here. Thus, anesthetics may target sleep circuits in all brains but only produce a slow wave in some. Instead, the mechanism for the reduced responsiveness in the central brain that we observed under isoflurane may be a combination of potentiated sleep circuits and compromised synaptic efficacy, which has been demonstrated in flies (van Swinderen and Kottler, 2014; Zalucki et al., 2015). While sleep circuits seem unlikely to modulate the responses of the peripheral On and Off pathways, the globally compromised synaptic efficacy could cause an imbalance in the responses of the On and Off pathways, resulting in the unexpected increase in power at f2.

Outlook

Bottom-up approaches that focus on molecular mechanisms have considerably improved our understanding of anesthetic drugs and have identified a promising set of potential target sites (Franks, 2008; Garcia et al., 2010; Brown et al., 2011). On the other hand, it remains unclear how effects at the molecular level affect large-scale neuronal circuits. Instead, top-down approaches that focus on global effects are providing evidence that general anesthetics share a common end point in the reduction of interarea communication (Lee et al., 2013; Mashour, 2014; Sarasso et al., 2015). Using the metrics developed for characterizing these global effects (Casali et al., 2013; Lee et al., 2015) in conjunction with the genetic manipulations available in Drosophila is a promising direction. Studies that manipulate the state of the brain and external perturbations can be combined with signal-processing techniques and modeling to help us understand how anesthetic effects at the molecular level change the global state of the brain.

Acknowledgments

Acknowledgments: We thank Angelique C. Paulk for helpful discussion and assistance using the multielectrode preparation.

Footnotes

  • The authors declare no competing financial interests.

  • N.T. was funded by an Australian Research Council Future Fellowship (FT120100619) and the Discovery Project (Grant DP130100194). B.V.S. was funded by National Health and Medical Research Council Project APP1103923.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    Alkire MT, Hudetz AG, Tononi G (2008) Consciousness and anesthesia. Science 322:876–880. doi:10.1126/science.1149213 pmid:18988836
    OpenUrlAbstract/FREE Full Text
  2. ↵
    Bates D, Mächler M, Bolker B, Walker S (2015) Fitting linear mixed-effects models using lme4. J Stat Softw 67. doi:10.18637/jss.v067.i01
    OpenUrlCrossRef
  3. ↵
    Behnia R, Clark DA, Carter AG, Clandinin TR, Desplan C (2014) Processing properties of ON and OFF pathways for Drosophila motion detection. Nature 512:427–430. doi:10.1038/nature13427 pmid:25043016
    OpenUrlCrossRefPubMed
  4. ↵
    Bendat JS, Piersol AG (2000) Random data: analysis and measurement procedures, Ed 3. New York: Wiley.
  5. ↵
    Bialek W, Rieke F, de van Steveninck RR, Warland D (1991) Reading a neural code. Science 252:1854–1857. pmid:2063199
    OpenUrlAbstract/FREE Full Text
  6. ↵
    Borst A, Euler T (2011) Seeing things in motion: models, circuits, and mechanisms. Neuron 71:974–994. doi:10.1016/j.neuron.2011.08.031
    OpenUrlCrossRefPubMed
  7. ↵
    Brown EN, Purdon PL, Van Dort CJ (2011) General anesthesia and altered states of arousal: a systems neuroscience analysis. Annu Rev Neurosci 34:601–628. doi:10.1146/annurev-neuro-060909-153200 pmid:21513454
    OpenUrlCrossRefPubMed
  8. ↵
    Capilla A, Pazo-Alvarez P, Darriba A, Campo P, Gross J (2011) Steady-state visual evoked potentials can be explained by temporal superposition of transient event-related responses. PLoS One 6:e14543. doi:10.1371/journal.pone.0014543
    OpenUrlCrossRefPubMed
  9. ↵
    Casali AG, Gosseries O, Rosanova M, Boly M, Sarasso S, Casali KR, Casarotto S, Bruno AM, Laureys S, Tononi G, Massimini M (2013) A theoretically based index of consciousness independent of sensory processing and behavior. Sci Transl Med 5:198ra105. doi:10.1126/scitranslmed.3006294
    OpenUrlAbstract/FREE Full Text
  10. ↵
    Cimenser A, Purdon PL, Pierce ET, Walsh JL, Salazar-Gomez AF, Harrell PG, Tavares-Stoeckel C, Habeeb K, Brown EN (2011) Tracking brain states under general anesthesia by using global coherence analysis. Proc Natl Acad Sci U S A 108:8832–8837. doi:10.1073/pnas.1017041108
    OpenUrlAbstract/FREE Full Text
  11. ↵
    Clark DA, Bursztyn L, Horowitz MA, Schnitzer MJ, Clandinin TR (2011) Defining the computational structure of the motion detector in Drosophila. Neuron 70:1165–1177. doi:10.1016/j.neuron.2011.05.023
    OpenUrlCrossRefPubMed
  12. ↵
    Egelhaaf M, Borst A (1989) Transient and steady-state response properties of movement detectors. J Opt Soc Am A 6:116–127. pmid:2921651
    OpenUrlCrossRefPubMed
  13. ↵
    Ferrarelli F, Massimini M, Sarasso S, Casali A, Riedner BA, Angelini G, Tononi G, Pearce RA (2010) Breakdown in cortical effective connectivity during midazolam-induced loss of consciousness. Proc Natl Acad Sci U S A 107:2681–2686. doi:10.1073/pnas.0913008107 pmid:20133802
    OpenUrlAbstract/FREE Full Text
  14. ↵
    Franks NP (2008) General anaesthesia: from molecular targets to neuronal pathways of sleep and arousal. Nat Rev Neurosci 9:370–386. doi:10.1038/nrn2372 pmid:18425091
    OpenUrlCrossRefPubMed
  15. ↵
    Garcia PS, Kolesky SE, Jenkins A (2010) General anesthetic actions on GABA(A) receptors. Curr Neuropharmacol 8:2–9. doi:10.2174/157015910790909502 pmid:20808541
    OpenUrlCrossRefPubMed
  16. ↵
    Haag J, Borst A (2008) Electrical coupling of lobula plate tangential cells to a heterolateral motion-sensitive neuron in the fly. J Neurosci 28:14435–14442. doi:10.1523/JNEUROSCI.3603-08.2008 pmid:19118177
    OpenUrlAbstract/FREE Full Text
  17. ↵
    Joesch M, Schnell B, Raghu S, Reiff DF, Borst A (2010) ON and OFF pathways in Drosophila motion vision. Nature 468:300–304. doi:10.1038/nature09545 pmid:21068841
    OpenUrlCrossRefPubMed
  18. ↵
    Kirszenblat L, van Swinderen B (2015) The yin and yang of sleep and attention. Trends Neurosci 38:776–786. doi:10.1016/j.tins.2015.10.001 pmid:26602764
    OpenUrlCrossRefPubMed
  19. ↵
    Kottler B, Bao H, Zalucki O, Imlach W, Troup M, van Alphen B, Paulk A, Zhang B, van Swinderen B (2013) A sleep/wake circuit controls isoflurane sensitivity in Drosophila. Curr Biol 23:594–598. doi:10.1016/j.cub.2013.02.021
    OpenUrlCrossRefPubMed
  20. ↵
    Lee U, Ku S, Noh G, Baek S, Choi B, Mashour GA (2013) Disruption of frontal-parietal communication by ketamine, propofol, and sevoflurane. Anesthesiology 118:1264–1275. doi:10.1097/ALN.0b013e31829103f5 pmid:23695090
    OpenUrlCrossRefPubMed
  21. ↵
    Lee U, Blain-Moraes S, Mashour GA (2015) Assessing levels of consciousness with symbolic analysis. Philos Trans A Math Phys Eng Sci 373:20140117 doi:10.1098/rsta.2014.0117
    OpenUrlAbstract/FREE Full Text
  22. ↵
    Lewis LD, Weiner VS, Mukamel EA, Donoghue JA, Eskandar EN, Madsen JR, Anderson WS, Hochberg LR, Cash SS, Brown EN, Purdon PL (2012) Rapid fragmentation of neuronal networks at the onset of propofol-induced unconsciousness. Proc Natl Acad Sci U S A 109:E3377–E3386. doi:10.1073/pnas.1210907109
    OpenUrlAbstract/FREE Full Text
  23. ↵
    Liu X, Lauer KK, Ward BD, Rao SM, Li SJ, Hudetz AG (2012) Propofol disrupts functional interactions between sensory and high-order processing of auditory verbal memory. Hum Brain Mapp 33:2487–2498. doi:10.1002/hbm.21385 pmid:21932265
    OpenUrlCrossRefPubMed
  24. ↵
    Mashour GA (2013) Cognitive unbinding: a neuroscientific paradigm of general anesthesia and related states of unconsciousness. Neurosci Biobehav Rev 37:2751–2759. doi:10.1016/j.neubiorev.2013.09.009 pmid:24076246
    OpenUrlCrossRefPubMed
  25. ↵
    Mashour GA (2014) Top-down mechanisms of anesthetic-induced unconsciousness. Front Syst Neurosci 8:115. doi:10.3389/fnsys.2014.00115 pmid:25002838
    OpenUrlCrossRefPubMed
  26. ↵
    Mitra P, Bokil H (2007) Observed brain dynamics. New York: Oxford UP.
  27. ↵
    Mitra PP, Pesaran B (1999) Analysis of dynamic brain imaging data. Biophys J 76:691–708. doi:10.1016/S0006-3495(99)77236-X pmid:9929474
    OpenUrlCrossRefPubMed
  28. ↵
    Murphy M, Bruno M-A, Riedner BA, Boveroux P, Noirhomme Q, Landsness EC, Brichant J-F, Phillips C, Massimini M, Laureys S, Tononi G, Boly M (2011) Propofol anesthesia and sleep: a high-density EEG study. Sleep 34:283. pmid:21358845
    OpenUrlPubMed
  29. ↵
    Norcia AM, Appelbaum LG, Ales JM (2015) The steady-state visual evoked potential in vision research: a review. J Vis 15(6):4 1–46. doi:10.1167/15.6.4
    OpenUrlCrossRef
  30. ↵
    Paulk AC, Zhou Y, Stratton P, Liu L, van Swinderen B (2013) Multichannel brain recordings in behaving Drosophila reveal oscillatory activity and local coherence in response to sensory stimulation and circuit activation. J Neurophysiol 110:1703–1721. doi:10.1152/jn.00414.2013 pmid:23864378
    OpenUrlAbstract/FREE Full Text
  31. ↵
    Paulk AC, Kirszenblat L, Zhou Y, van Swinderen B (2015) Closed-loop behavioral control increases coherence in the fly brain. J Neurosci 35:10304–10315. doi:10.1523/JNEUROSCI.0691-15.2015 pmid:26180205
    OpenUrlAbstract/FREE Full Text
  32. ↵
    Regan MP, Regan D (1988) A frequency domain technique for characterizing nonlinearities in biological systems. J Theor Biol 133:293–317. doi:10.1016/S0022-5193(88)80323-0
    OpenUrlCrossRef
  33. ↵
    Reiff DF, Plett J, Mank M, Griesbeck O, Borst A (2010) Visualizing retinotopic half-wave rectified input to the motion detection circuitry of Drosophila. Nat Neurosci 13:973–978. doi:10.1038/nn.2595 pmid:20622873
    OpenUrlCrossRefPubMed
  34. ↵
    Reisenman C, Haag J, Borst A (2003) Adaptation of response transients in fly motion vision. I: experiments. Vision Res 43:1293–1309. doi:10.1016/S0042-6989(03)00091-9
    OpenUrlCrossRef
  35. ↵
    Sarasso S, Rosanova M, Casali AG, Casarotto S, Fecchio M, Boly M, Gosseries O, Tononi G, Laureys S, Massimini M (2014) Quantifying cortical EEG responses to TMS in (un)consciousness. Clin EEG Neurosci 45:40–49. doi:10.1177/1550059413513723 pmid:24403317
    OpenUrlAbstract/FREE Full Text
  36. ↵
    Sarasso S, Boly M, Napolitani M, Gosseries O, Charland-Verville V, Casarotto S, Rosanova M, Casali A, Brichant J-F, Boveroux P, Rex S, Tononi G, Laureys S, Massimini M (2015) Consciousness and complexity during unresponsiveness induced by propofol, xenon, and ketamine. Curr Biol 25:3099–3105. doi:10.1016/j.cub.2015.10.014 pmid:26752078
    OpenUrlCrossRefPubMed
  37. ↵
    Supp GG, Siegel M, Hipp JF, Engel AK (2011) Cortical hypersynchrony predicts breakdown of sensory processing during loss of consciousness. Curr Biol 21:1988–1993. doi:10.1016/j.cub.2011.10.017 pmid:22100063
    OpenUrlCrossRefPubMed
  38. ↵
    Truccolo WA, Ding M, Knuth KH, Nakamura R, Bressler SL (2002) Trial-to-trial variability of cortical evoked responses: implications for the analysis of functional connectivity. Clin Neurophysiol 113:206–226. pmid:11856626
    OpenUrlCrossRefPubMed
  39. ↵
    van Swinderen B (2006) A succession of anesthetic endpoints in the Drosophila brain. J Neurobiol 66:1195–1211. doi:10.1002/neu.20300 pmid:16858694
    OpenUrlCrossRefPubMed
  40. ↵
    van Swinderen B, Kottler B (2014) Explaining general anesthesia: a two step hypothesis linking sleep circuits and the synaptic release machinery. Bioessays 36:372–381. doi:10.1002/bies.201300154 pmid:24449137
    OpenUrlCrossRefPubMed
  41. ↵
    Wang X, Chen Y, Ding M (2008) Estimating Granger causality after stimulus onset: a cautionary note. Neuroimage 41:767–776. doi:10.1016/j.neuroimage.2008.03.025 pmid:18455441
    OpenUrlCrossRefPubMed
  42. ↵
    Yekutieli D, Benjamini Y (1999) Resampling-based false discovery rate controlling multiple test procedures for correlated test statistics. J Stat Plan Inference 82:171–196. doi:10.1016/S0378-3758(99)00041-5
    OpenUrlCrossRef
  43. ↵
    Zalucki OH, Menon H, Kottler B, Faville R, Day R (2015a) Syntaxin1A-mediated resistance and hypersensitivity to isoflurane in Drosophila melanogaster. Anesthesiology 122:1060–1074.
    OpenUrl
  44. ↵
    Zalucki OH, Day R, Kottler B, Karunanithi S, van Swinderen B (2015b) Behavioral and electrophysiological analysis of general anesthesia in 3 background strains of Drosophila melanogaster. Fly 9:7–15.
    OpenUrl

Synthesis

The decision was a result of the Reviewing Editor Gustavo Deco and the peer reviewers coming together and discussing their recommendations until a consensus was reached. A fact-based synthesis statement explaining their decision and outlining what is needed to prepare a revision is listed below. The following reviewers agreed to reveal their identity: Karin Nordström, Miguel Dasilva

Review 1:

The first para in the Introduction is too vague to be useful. For example, the sentence on lines 41-44, does not really say anything, but is then continued with "..to reconcile these varied observations..." Which varied observations are you referring to?

In general in the introduction it is very unclear what species you are referring to. I assume that the second paragraph refers to humans, but this is never specified. Since you are comparing flies with humans, it is paramount that you specify who you are talking about.

I am not convinced about your grouping of experiments into ipsilateral and contralateral, described in Figure 3a, as the stimuli are quite different (no stimuli on half the screen, vs different stimuli). The grouping may be perfectly sound, but needs to be justified, with data.

The main point (based on the title) of the paper is described on page 24, i.e. after many pages of reading. I think it is a bit of a shame that it takes so long to get to this, and that many busy readers may thus miss it. Is there any way to simplify figures 1-2 and the associated text (as these are mostly methods figures), so that readers do not give up before getting to the results?

Lines 552, 553 and 554 have superscript letters after 'period'. What do these mean?

The paragraph on lines 558-575 seems to belong under the next heading.

On line 624 you introduce biological non-linearities. All fly vision people are drilled in the knowledge that the visual system is full of non-linearities, so you might want to rephrase this not to annoy this audience.

Again on line 647, it is well established that the ON and OFF pathways are asymmetric (See any recent paper from the Reiser, Clandinin or Borst lab). If you do not want to annoy people from the fly vision community, rephrase.

Lines 671-674, your reasoning here is not very clear. Please clarify.

Additionally I do not understand why the f1 harmonic is so large in the periphery for contralateral stimulation (Figure 3b and 3c)? Could that be an effect of your definition into contralateral and ipsilateral, or how do you think this signal has traveled from the contralateral visual field to the medulla on the other side?

Figure 2a: blue on black does not print well at all. Can you modify?

Figure 2c, please include a scale bar for the y-axis.

Can you make figure 4b bigger, as it is an important control (you can make figure 4f smaller).

Review 2:

In my opining there are some concerns that should be explained by the authors:

1. The classification of flicker configurations as ipsilateral or contralateral to the insertion site seems controversial and would need further explanation. It would be most helpful to see that ipsilateral and bilateral stimulations at the same frequency do not show significant differences. In that regard, further explanation would be required on how trials where different tagging frequencies were used at both sites were classified as ipsi or contralateral. Proof of their similarity with those trials where one single frequency was presented at both sites would also be needed.

2. The authors claim that the second harmonic reflects more local processing, because f2 is stronger when the flicker is shown to the ipsilateral site. However, according to their classification of the different types of flickers, both ipsilateral-only and bilateral-at-the-same-frequency are included in this category. This, at least, should be discussed. Also, why do f2-ipsilateral and f1-contralateral (Figure 3b) have similar power?

3. Is the Chi-square test the best way to show significant differences in all cases? What I mean is that it is difficult to follow whether statistical significance accounts for one specific frequency or both. For example, in figure 3c, why is an ANOVA not used to show the differences between central and peripheral locations? Please, explain.

4. In terms of the minimal model to explain the effects of isoflurane on SSVEP power, all presented data refer to ipsilateral stimulation, could the model account for contralateral stimulation too? Did the authors consider testing experimentally the prediction that isoflurane changed the balance of the ON, OFF responses? Why not? This would need some explanation.

5. What is the reason for applying the minimal model based on SNRs? Can the authors justify from a physiology-related perspective the necessity of using this model in the same way they did for the previous one related to SSVEP power?

6. The authors claim that they distinguish between local and global effects of isoflurane. I wonder whether it would be more precise talking about central and peripheral effects. Could the authors, please justify the use of the former term?

7. In the discussion (line 824), the authors admit that "the broad-field flicker prevents us from precisely disentangling the relative contribution of each optic lobe to the LFP" Could the authors explain whether this represents a concern when establishing the difference between ipsi and contralateral stimulation in terms of the scope of the paper?

Back to top

In this issue

eneuro: 3 (4)
eNeuro
Vol. 3, Issue 4
July/August 2016
  • Table of Contents
  • Index by author
Email

Thank you for sharing this eNeuro article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Local Versus Global Effects of Isoflurane Anesthesia on Visual Processing in the Fly Brain
(Your Name) has forwarded a page to you from eNeuro
(Your Name) thought you would be interested in this article in eNeuro.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Local Versus Global Effects of Isoflurane Anesthesia on Visual Processing in the Fly Brain
Dror Cohen, Oressia H. Zalucki, Bruno van Swinderen, Naotsugu Tsuchiya
eNeuro 13 July 2016, 3 (4) ENEURO.0116-16.2016; DOI: 10.1523/ENEURO.0116-16.2016

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Share
Local Versus Global Effects of Isoflurane Anesthesia on Visual Processing in the Fly Brain
Dror Cohen, Oressia H. Zalucki, Bruno van Swinderen, Naotsugu Tsuchiya
eNeuro 13 July 2016, 3 (4) ENEURO.0116-16.2016; DOI: 10.1523/ENEURO.0116-16.2016
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Significance Statement
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Footnotes
    • References
    • Synthesis
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • anesthesia
  • consciousness
  • Drosophila
  • frequency tagging
  • isoflurane
  • SSVEP

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

New Research

  • Deciding while acting - Mid-movement decisions are more strongly affected by action probability than reward amount
  • CaMKIIα promoter-controlled circuit manipulations target both pyramidal cells and inhibitory interneurons in cortical networks
  • Gas7 is a novel dendritic spine initiation factor
Show more New Research

Cognition and Behavior

  • Deciding while acting - Mid-movement decisions are more strongly affected by action probability than reward amount
  • Environment Enrichment Facilitates Long-Term Memory Consolidation Through Behavioral Tagging
  • Effects of cortical FoxP1 knockdowns on learned song preference in female zebra finches
Show more Cognition and Behavior

Subjects

  • Cognition and Behavior

  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Latest Articles
  • Issue Archive
  • Blog
  • Browse by Topic

Information

  • For Authors
  • For the Media

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
  • Feedback
(eNeuro logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
eNeuro eISSN: 2373-2822

The ideas and opinions expressed in eNeuro do not necessarily reflect those of SfN or the eNeuro Editorial Board. Publication of an advertisement or other product mention in eNeuro should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in eNeuro.