Abstract

Introduction:

Volatile anesthetics, such as isoflurane, generate a state of unconsciousness and analgesia across the animal kingdom and are widely used in clinical settings. Yet, anesthetic mechanisms are poorly understood: the volatile anesthetics are so profligate in their potential effects that it has proven difficult to determine which actions are most causal at the systems level.

Methods:

To test if specific cellular targets mediate the anesthetic effect across a complete, intact nervous system, we imaged neuron activity in the Caenorhabditis elegans head ganglia at cellular resolution. We measured the effect of increasing anesthetic concentrations across a range of identified neurons within the C. elegans nervous system.

Results:

However, rather than dramatic effects on any particular neuronal class, we measured uniform suppression of both neuron activity and connectivity with increasing isoflurane across the nervous system. We find the degree of activity suppression to be proportional to the baseline activity of the neuron in the awake state. Within this context, highly connected neurons, specifically neurons with high in-degree connectivity, are inherently active and display large activity suppression. These include hub interneurons within the C. elegans command locomotory circuit that control behavioral crawling states and contribute to system-wide coherence of neuron dynamics. By analyzing the effect of isoflurane on the activity of two specific hub interneuron classes, AVA and AVE, we show that the large degree of suppression observed in these neurons corresponds to high baseline activity.

Discussion:

Exploiting the small size, simplicity and optical accessibility of C. elegans, our results demonstrate that isoflurane anesthesia globally suppresses activity and connectivity across a wide range of neuron types, and suggest a model of anesthesia in which proportional suppression of activity results in disruption of highly connected, highly active, hub loci that are critical to nervous system coordination and state dynamics.

Introduction

Volatile anesthetics are inhaled drugs that produce a reversible state of unconsciousness and analgesia in organisms spanning the animal kingdom. Despite their ubiquitous clinical usage, no consensus exists on the molecular and cellular action of these drugs. Volatile anesthetics are profligate, with evidence implicating many putative postsynaptic and presynaptic targets. These include agonism of inhibitory GABAA and glycine receptors (Mascia et al., 2000; Mihic et al., 1997; Krasowski and Harrison, 2000), inhibition of excitatory NMDA glutamate (Martin et al., 1995; Solt et al., 2006) and nicotinic acetylcholine receptors (Flood et al., 1997; Xu et al., 2000) as well as disruption of presynaptic neurotransmitter release (Maclver et al., 1996; Perouansky et al., 1995; van Swinderen et al., 1999; Troup et al., 2019; Zalucki et al., 2015). Defects in mitochondrial complex I produce volatile anesthetic hypersensitivity in nematodes, mice and humans (Kayser et al., 2001; Quintana et al., 2012; Zimin et al., 2016), while halogenated ethers have been shown to modulate inhibitory leak current through K+ channels (Heurteaux et al., 2004; Patel et al., 1999). But what are the causal actions of volatile anesthetics across an intact nervous system and, more specifically, do certain neuron-types mediate the effects? These questions remain critical to understanding the effects of volatile anesthetics on nervous system function at the network level.

The mechanisms of volatile anesthetics can be studied through their functional effects on neuronal dynamics and network connectivity. While medical EEG monitoring provides a quantifiable measurement of brain activity associated with anesthesia, the EEG has limited ability to reveal functional changes at the level of individual neurons and simple circuits. By contrast, modern neuro-imaging techniques such as fMRI and fluorescence microscopy, as well as higher resolution and multi-neuron electrical recordings, coupled with quantitative network analyses, have begun to provide a system-wide understanding of the functional effects of volatile anesthetics. Broadly, the anesthetized state displays an overall suppression in network efficiency and a breakdown in information transfer between network nodes, or “hubs” (Mashour, 2017). Intracortical recordings in mice under desflurane anesthesia reveal a breakdown of global functional connectivity resulting from accumulated disruption in local circuits (Hudetz et al., 2020). Our lab has performed similar analyses in C. elegans using single-neuron resolution calcium imaging, showing that the anesthetized state is characterized by a collapse of system-wide functional organization as measured by neuron pair correlations (Awal et al., 2020; Awal et al., 2018).

The nematode worm C. elegans is an excellent model system for the study of the basic mechanisms of anesthesia given its simple, completely mapped neuro-connectome (Cook et al., 2019), capabilities for comprehensive multi-neuron functional imaging, and behavioral and neurological responses to volatile anesthetics that parallel higher organisms. Its nervous system holds many highly conserved neurotransmitters in common with humans, including glutamate, GABA, and acetylcholine, as well as homologous receptors (Hobert, 2013; Luersen et al., 2016). As with nervous systems across the animal kingdom, its connectome exhibits “rich club” architecture, possessing a cohort of highly-connected hub neurons which are key to maintaining system-wide coherence and function (Towlson et al., 2013; Uzel et al., 2022).

Each of 302 neurons in the hermaphrodite C. elegans has been characterized anatomically, uniquely named, and pan-neuronal gene expression and neurotransmission properties mapped (Taylor et al., 2021; Gendrel et al., 2016; Pereira et al., 2015; Serrano-Saiz et al., 2013). The animal’s small size and transparency allow for volumetric calcium imaging of neuron activity across the entire head of the animal with single cell resolution. The NeuroPAL fluorescent marker system allows for identification of individual neurons within the animals completely mapped neuro-connectome through a combination of differentially expressed nuclear fluorophores and standardized anatomic relationships (Yemini et al., 2021). Thus, identifiable neurons within system-wide multi-neuron activity data sets can be assessed based on innate properties such as neurotransmitter, receptor expression, and synaptic connectivity. These capabilities allow us to systematically investigate how isoflurane anesthesia alters activity patterns of individual neurons within an intact nervous system with a degree of fidelity not possible in higher organisms. In this study, we exploit these advantages of C elegans to measure the functional effects of isoflurane at the neuron and network level and understand how this leads to the alterations in global network dynamics characteristic of the anesthetized state.

Materials and methodsC. elegans strains and maintenance

C. elegans were raised at 20 °C on nematode growth medium agar seeded with a lawn of OP50 E. coli. All animals used in this study were a cross between two transgenic strains: QW1217; OH15263. The first transgenic strain, QW1217 (zfIs124[Prgef-1: GCaMP6s]; otIs355[Prab-3: NLS:tagRFP]), expresses nuclear-localized tagRFP, a fluorophore facilitating neuron tracking, and cytoplasmic GCaMP6s, a fluorescent calcium reporter facilitating capture of neuronal activity (Chen et al., 2013). The second transgenic strain is OH15263 (otIs670), which expresses the NeuroPAL genetic cassette (Yemini et al., 2021). The NeuroPAL cassette induces distinct patterns of expression of three nuclear-localized fluorophores (mNeptune2.5, CyOFP1, and mTagBFP2) across the C. elegans nervous system. The expression of these fluorophores is controlled by a set of promoters, designed to allow neuronal identification to be performed by comparing the relative locations of neurons and relative expression (measured via fluorescence) of the three fluorophores. QW1217 was obtained as a gift from the Alkema Lab (University of Massachusetts Medical School, Worcester, Massachusetts) and OH15263 was obtained from the Caenorhabditis Genetics Center (University of Minnesota, Minneapolis, Minnesota). All C. elegans used in this study were young adult hermaphrodites.

Animal imaging preparation

Animals were paralyzed prior to isoflurane exposure and imaging through immersion in S-basal buffer (100 mM NaCl, 50 mM KPO4 buffer, 5 ug/mL cholesterol) containing 5 mM tetramisole, a selective levamisole receptor agonist. Levamisole receptors are a subtype of nematode-specific nicotinic acetylcholine receptor that, in the C. elegans, are largely expressed at the neuromuscular junction of body wall muscles, although some neuronal expression has been reported (Culetto et al., 2004). Tetramisole exposure induces muscle depolarization, causing spastic contraction and eventually rigid paralysis (Martin et al., 2012). At the concentrations used, paralysis occurs within a few minutes. We have previously shown that at the levels used, C. elegans can be maintained under tetramisole paralysis for at least 2 h while maintaining consistent neuronal dynamics (Chang et al., 2023). After paralysis, animals were embedded in a transparent polyethylene glycol hydrogel before being subjected to anesthetic exposure and imaging procedures (Burnett et al., 2018).

Isoflurane anesthesia and imaging procedure

Hydrogel embedded animals were equilibrated to specific atmospheric concentrations of isoflurane to produce stable depths of anesthesia. To induce anesthesia to a given atmospheric percentage of isoflurane, we first equilibrate 50 mL of S-basal buffer to that percentage of isoflurane by exposing the buffer to a surface atmosphere of isoflurane at that percentage in a sealed petri dish for 30 min. The hydrogel embedded animal is then submerged in the pre-equilibrated buffer pool, and the pool and animal then re-exposed to the specified percentage of atmospheric isoflurane for an additional 30 min. All S-basal used in this procedure includes 5 mM tetramisole to maintain muscular paralysis.

We have previously shown that 50% of animals no longer respond to sharp touch when anesthetized to 3% atmospheric isoflurane using this procedure, a level which is analogous to a surgical plane of anesthesia in humans (Awal et al., 2018). Therefore, atmospheric concentrations of 2 and 4% isoflurane were chosen to, respectively, represent light and moderate depths of anesthesia. All animals in this study were subjected to a progressive anesthetization and imaging procedure. After paralysis and hydrogel encapsulation, animals were rested for 30 min, and then imaged for 10 min, generating a baseline “0%” recording. Animals were then exposed to 2% isoflurane for 30 min, followed by 10 min of imaging, then exposed to 4% isoflurane for 30 min, followed by a final 10 min of imaging.

Light sheet microscopy

Volumetric fluorescent microscopy was performed using a Dual Inverted Selective Plane Illumination (diSPIM) fluorescence microscope (Applied Scientific Instrumentation, Eugene, OR) with 0.8 NA 40x water immersion objectives (Nikon USA, Melville, NY; Awal et al., 2020). To capture neuronal activity, tagRFP and GCaMP6s fluorescence were imaged in the head region of each animal (n = 15) for three sessions of 10 min each, at 0, 2, and 4% isoflurane concentrations. Volumes were captured at a rate of 2 Hz with a voxel size of 0.1625 × 0.1625 × 1 μm. Before and after each activity imaging session, a single volume capturing tagRFP and the three NeuroPAL fluorophores (mNeptune2.5, CyOFP1, and mTagBFP2) was obtained to facilitate neuronal identification. An example image displaying the maximum intensity projection of a C. elegans head with NeuroPAL fluorophores is displayed in Figure 1A.

Capture of neuronal dynamics in C. elegans at varying depths of isoflurane anesthesia. (A) Maximum intensity projection image of the head of a C. elegans showing false-color representation of the NeuroPAL fluorescence markers (mTagBFP2 is shown in blue, CyOFP1 is shown in green, and mNeptune2.5 is shown in red) used for neuron identification. Identified neurons are as labeled. The enclosures delineate ganglia including the anterior pharyngeal bulb (yellow); the anterior ganglion (magenta); the dorsal ganglion (turquoise); the lateral ganglia (scarlet). The ventral ganglion, which is located more caudally is not shown. (B) Neuronal dynamics as measured by GCaMP6s fluorescence (ΔF/F normalized) captured throughout head ganglia in an individual exemplar animal during three distinct states: awake, mildly anesthetized (equilibrated to 2% atmospheric isoflurane), and moderately anesthetized (equilibrated to 4% atmospheric isoflurane). In between each 10-min recording period, animals were equilibrated to the next greater level of anesthesia for 30 min. Only neurons that were able to be identified via NeuroPAL in at least one recording were shown. Successful identification in one recording did not guarantee identification in any other recording, leading to gaps in data as shown. Citation: Chang et al (2025). License: CC BY. Source: PLoS One.

Signal extraction and neuron identification

RFP-labeled neuronal nucleus tracking and extraction of GCaMP6s fluorescence signal from neuronal soma was performed as described in previous work from our group (Awal et al., 2020). 150 neuronal activity traces were extracted from each 10-min volumetric recording using custom Python and MATLAB (Mathworks, USA) scripts. Neuron ID was performed using pre and post-imaging NeuroPAL volumes. Neurons in the pre and post imaging volumes corresponding to each 10-min activity imaging session were identified using the NeuroPAL method (Yemini et al., 2021), considering relative anatomical position of neurons and relative brightness of the three NeuroPAL fluorophores in each neuronal nucleus. Neurons in which putative IDs in both pre and post imaging NeuroPAL volumes matched and in which the specified nuclei in both pre and post volumes could be verifiably linked to the same neuron in the corresponding activity imaging session were considered positively IDed.

Out of 45 imaging sessions (n = 15 animals at three anesthesia conditions each), an average of 44.6 neurons were positively IDed per session, with a minimum of 21 and a maximum of 69 (excluding one session each at 2% isoflurane 4% isoflurane in which no positive identifications could be made, and one session at 4% in which only one positive identification could be made due to insufficient image quality). A total of n = 991 neurons were positively identified at any of the three anesthesia conditions, with n = 760/678/448 neurons being identified at the 0%/2%/4% isoflurane conditions. 256 neurons were positively IDed across all three conditions. Figure 1A shows the IDed GCaMP activity traces obtained from a single animal across the three conditions. Because only the head of each animal was imaged, we only considered neurons of the head ganglia and the most anterior neurons of the ventral nerve cord to be possible targets for identification. Figure 2 summarizes the number of times each neuron type in the C. elegans head (segregated by ganglion) was positively IDed across all 15 animals. Supplementary Figure S1 may be referenced as a comprehensive survey of neuron identification for the dataset analyzed in this manuscript, providing specification and quantification of which neurons classes were identified in which animals and recording conditions.

Survey of neurons identified across recording conditions. Color indicates the number of times each neuron class was positively identified in each of the three recording conditions (0, 2, and 4% atmospheric isoflurane). Each neuron class could have been identified as many as 15 times across the sample of 15 animals. A gray entry signifies that a neuron class was never positively identified. Red text labels indicate neurons classified as interneurons within the C. elegans connectome.

Data analysis

IDed GCaMP6s signal traces were individually normalized before analysis to ΔF/F0, with F0 calculated as the mean value of the lowest 1% of fluorescence observed in each trace. Because fluorophores tend to fluoresce particularly brightly when first illuminated, the first 30 s of signal recorded from each imaging session were excised prior to normalization and were therefore omitted from all subsequent analysis. To assess levels of overall neuronal activity within each IDed signal trace, “signal variance” (SV), or the standard deviation of the normalized signal, were calculated for each trace. Neuronal traces were assorted based on features derived from neuron ID: neurotransmitter identity (Pereira et al., 2015; i.e., “ergic” type) and rich club membership (Towlson et al., 2013; Uzel et al., 2022). C. elegans possesses 11 neurons that satisfy the statistical properties of a rich club: they are a set of high-degree neurons which are more highly interconnected than would be expected based on random chance alone. To assess isoflurane induced changes in neuronal activity levels, we assessed the difference in signal variance (Δ signal SD), or “activity shift” for each neuron, between the anesthetized states (2 and 4% isoflurane) and the awake state (0% isoflurane).

To assess isoflurane induced changes in synchrony between synaptically linked neuron-pairs, we generated a list of all pairs of IDed neurons putatively linked by chemical synapses using the hermaphrodite connectome generated by Cook et al. (2019) Because we were primarily interested in drug exposure induced changes, we generated subset populations of neuron-pairs found at both 0 and 2% isoflurane conditions (n = 1,418) and at both 0 and 4% isoflurane conditions (n = 672). These neuron pairs were then assorted further based on neurotransmitter identity of the presynaptic partner as derived from neuron ID (Pereira et al., 2015). Synchrony between neuron pairs was calculated as the Pearson correlation coefficient (PCC) between the differentiated activity signals of pair members.

To investigate the interaction between neuron connectivity and isoflurane-induced neuronal activity change, we determined the in-degree and out-degree of each neuron type using the hermaphrodite chemical connectome (Cook et al., 2019). In-degree was defined as a neuron’s number of presynaptic partners and out-degree was defined as a neuron’s number of post-synaptic partners.

Statistical methods

For dataset validation, to quantify how our identified neuron sample differs from the population of all tracked neurons, we applied two-way mixed ANOVA with signal variance as outcome. One-way ANOVA was applied to assess for differences in mean signal variance across animals within each recording condition.

One-way ANOVA was used to compare difference in mean signal variance across anesthesia conditions for different neuronal subgroups. Post-hoc pairwise comparisons were performed using the Tukey–Kramer method. Differences between pre- and post-anesthesia PCC distributions were performed using the Kolmogorov–Smirnov test, and separation of means was quantified using Cohen’s d.

Correlation between neuronal in-degree or out-degree and activity shift was determined using Spearman’s rank correlation coefficient. Comparison of activity shift between rich club neurons and all other neurons was performed using the Wilcoxson rank-sum test. A non-parametric test was chosen because of a large discrepancy between sample sizes of compared populations.

Pre- and post-anesthesia signal variance within individual neuron classes (AVA and AVE) was compared using the paired-sample t-test. Correlation between baseline signal variance and isoflurane-induced activity shift was assessed using linear regression.

Code accessibility

Custom image and data analysis code as well as original imaging data sets are available from the authors upon request.

ResultsNeuron identification and imaging the response to isoflurane anesthesia

There is evidence that volatile anesthetics interact with a large array of neurotransmitter systems, both postsynaptically and presynaptically (Chau, 2010; Weinrich and Worcester, 2018). However, it remains unclear if there is a selective effect of isoflurane on any distinct class of neurons that drives the functional effects of anesthesia. Employing fluorescence microscopy in C. elegans, we can simultaneously image individual neuron activity across many neurons in the animal’s head and identify them within the complete connectome using the NeuroPAL marker system. Figure 1A illustrates an example image of the NeuroPAL fluorescence markers and identified neuron labels within the head of C. elegans. Figure 1B illustrates the calcium fluorescence (ΔF/F normalized GCaMP6s) of identified neurons from an individual example animal during 10 min trials of spontaneous activity at progressively higher levels of anesthesia. In order to effectively measure the functional elements that initiate and drive anesthesia (i.e., those that are affected first, most dramatically), we focused our imaging assays over a range of mild to moderate levels of isoflurane [i.e., 0, 2, 4% atmospheric isoflurane, as we determined earlier (Awal et al., 2018)]. Figure 2 illustrates the rate at which each neuron type was successfully identified within our imaging assays, using the NeuroPAL fluorescence markers (see methods) over 15 separate animals. Using this system we successfully identified a majority of neurons in the animals head (63.7% of neurons were IDed at least one time) and identified >30% of neurons in at least five imaging trials. Thus, while identification rates were not comprehensive within individual trials, identified neurons within the complete data set sufficiently cover a broad census of neurons and synaptic subtypes (see Supplementary Figure S1). In particular, we routinely captured numerous interneurons, indicated in red in Figure 2.

The average calcium fluorescence across all imaged neurons within each animal (n = 15) is shown in Figure 3A, with a marked pattern of decreasing fluorescence indicating an increase in neuronal quiescence at greater depths of anesthesia. To further assess the population of neurons we successfully IDed via NeuroPAL compared to the ensemble effects of anesthesia, we compared the distribution of neuronal activity in the set of all tracked neurons to that of the IDed neurons. Neuron activity is quantified by the standard deviation of the ΔF/F normalized GCaMP6s fluorescence signal for each neuron over a given recording block (i.e., the “signal variance”). The distribution and mean across all imaged neurons and the subset of IDed neurons at each condition are shown in Figure 3B. There is a small observable difference in median signal variance (SV), with a lower median SV in the IDed subpopulation compared to the tracked superset in all three anesthesia conditions (p values between the mean of all imaged neurons and the IDed subset were 0%: 5 × 10−6, 2%: 4 × 10−4, 4%: 2 × 10−11, via one sample t-test). However, SV interquartile ranges are grossly similar when comparing all tracked neurons and IDed neurons across conditions.

Neuronal activity signal quantified by signal variance across the set of identified and unidentified neurons, and across animals. (A) Mean neuronal GCaMP fluorescence signal as recorded from each of n = 15 hermaphrodite C. elegans across specified depths of anesthesia. Neuronal activity has been averaged across 150 neurons independently tracked from each animal at each depth. (B) Violin plots and histograms summarizing distribution of neuronal activity signal as quantified by signal variance (standard deviation of GCaMP fluorescence signal). The distribution of the population of all neurons recorded at each anesthesia recording condition (0, 2, and 4% atmospheric isoflurane) are shown, as well as the distribution of the subset of neurons successfully identified at each condition. White circles show the population median and shaded bars show the interquartile range. (C) Comparison of signal variance distribution and means between the population of identified neurons and non-identified neurons at each recording condition. Means are shown by line plots. Boxplots show median and interquartile range (IQR), with whiskers showing non-outlier minima and maxima, with outliers being values > 1.5 times the IQR away from the IQR. Individual outliers have not been plotted. (D) Distribution of signal variance across recording conditions and animals. Boxplots show median and interquartile range (IQR), with whiskers showing non-outlier minima and maxima. Individual outlier data points have been plotted. See main text for statistical analysis.

As illustrated in Figure 3C, we further compared activity measurements and fluorescence SV, between the subsets of IDed neurons and non-IDed neurons. Although SV distributions are qualitatively largely similar, and while both IDed and non-IDed populations exhibited a dose-dependent depression of SV in response to isoflurane, the IDed population means were consistently lower than the non-IDed population mean (IDed vs. non-IDed, 0%: 0.107 SV vs. 0.135 SV, 2%: 0.083 SV vs. 0.098 SV, 4%: 0.062 SV vs. 0.083 SV). To quantify these differences, we performed two-way mixed ANOVA assessing IDed vs. not IDed and isoflurane condition (0%/2%/4%) as predictive factors for neuronal SV. As expected, depth of anesthesia was significantly predictive for SV, with p = 3 × 10−56. The IDed vs. non-IDed population was also shown to be predictive for SV with p = 2 × 10−17, suggesting that the IDed population does differ significantly from the non-IDed population in terms of SV. However, interaction between the factors of being IDed and anesthesia depth was not significant at p = 0.085, which crucially suggests that the rate of SV change in response to anesthetic is consistent across the IDed and non-IDed populations. This can be visually observed in Figure 3C, in which the lines depicting the change in SV from 0 to 2% isoflurane and 2 to 4% isoflurane are largely parallel between the IDed and non-IDed populations. Thus while the sample of IDed neurons is not akin to an ideal random sampling from the overall set of tracked neurons, the population of IDed and non-IDed neurons react to isoflurane in a similar way in terms of SV. Therefore, analysis of the metric of SV in response to isoflurane exposure in the IDed subpopulation of neurons should be representative of SV response to isoflurane in neurons in C. elegans overall.

To further compare the response across individual animals, we plotted separately the distribution of SV in IDed neurons of all 15 animals and 3 recording conditions in Figure 3D. While the neuron populations at 0% isoflurane exhibit comparatively large numbers of highly active outlier neurons, this behavior is consistent across all 15 animals. Indeed, one-way ANOVA comparing the mean SV across all 15 animals at the 0% isoflurane recording condition does not reject the null hypothesis that means are equal, with p = 0.097. The mean SV across animals does appear to differ at the 2% (p = 3 × 10−11) and 4% (p = 3 × 10−4) isoflurane conditions. However, some variation between animals in response to isoflurane is expected and consistent with previous observations. While some animals clearly differ from others in terms of distribution, given the largely overlapping interquartile ranges and similar counts of IDed neurons across animals, there do not appear to be gross outlier effects from any specific animal or few animals.

The response of neurotransmitter subpopulations to isoflurane anesthesia

Leveraging the cell-specific information of the C. elegans connectome, we can categorize the identified neurons in our imaging data sets by either their presynaptic or postsynaptic properties and assess how neurotransmitter identity corresponds to the degree to which that neuron’s activity is altered in the anesthetized state. As in higher organisms, the majority of neurons in C. elegans only express the molecular machinery required to produce one neurotransmitter type (Pereira et al., 2015). This establishes a clean division of the captured neurons into largely mutually exclusive groups. Figure 4A and Table 1 show how different subpopulations of neurons defined by neurotransmitter identity behave under isoflurane anesthesia. As shown in Figure 4A, left, we observe that the dynamic range of neuronal activity, measured by signal variance (SV), becomes suppressed with progressively deeper planes of anesthesia when measured across all identified neurons. These results are similar to that of ensemble measurements for all tracked neurons in the data set (SV for all neurons including those tracked and measured but not identified: 0.125 SV at 0% isoflurane, 0.093 SV at 2% isoflurane, and 0.079 at 4% isoflurane) suggesting that our IDed neuron data set accurately represents the effects of isoflurane on activity across the entire nervous system. These results also recapitulate our previous findings measured across all tracked neurons in the C. elegans head (Awal et al., 2020). Similar observations have been documented in mammals in which isoflurane anesthesia is accompanied by overall reduction in activity (Ying et al., 2009). Comparing neuronal subpopulations based on neurotransmitter type, we find that acetylcholinergic neurons exhibit potent isoflurane-induced suppression of activity, with a 28% drop in signal variance (SV) at 2% (0.103 SV) isoflurane and a 55% drop at 4% (0.064 SV) isoflurane compared to controls (0% 0.142 SV). Glutamatergic neurons exhibit a 20% drop in SV at 2% isoflurane (0.089 to 0.071 SV) and a 31% drop SV at 4% isoflurane (to 0.061 SV). GABAergic neurons exhibit a 22% drop in SV from 0% (0.152) to 2% (0.118) isoflurane (although this is not statistically significant by ANOVA analysis) and a 40% drop at 4% isoflurane (0.091 SV). Strikingly, dopaminergic neurons, which are sensory and show low baseline activity levels at 0% isoflurane, exhibits no significant change in SV across all depths of anesthesia.

Suppression of neuronal activity and synchrony of distinct neuron populations under isoflurane anesthesia. (A) Neuronal signal variance measured as standard deviation of ΔF/F normalized calcium signal. Included are all positively identified neurons recorded from n = 15 animals, sorted by neurotransmitter identity and level of anesthesia (n = 760/678/448 neurons at 0%/2%/4% isoflurane). Red bars mark population means. Difference of means across anesthesia conditions for each subgroup of neurons were compared using one-way ANOVA, and post-hoc pairwise comparisons were performed using the Tukey–Kramer method. See Table 1 for statistical detail. * p < 0.05; ** p < 0.01; *** p < 0.001 (B) Distribution of signal synchrony of neuron pairs connected by chemical synapses as measured by Pearson correlation coefficient (PCC). Only neuron pairs that were positively identified in both the awake condition and in the respective comparison condition (either 2% or 4% isoflurane) were included. Correlation distributions from neuron pairs drawn from awake animals (grey histograms) were compared against those drawn from animals equilibrated 2 and 4% isoflurane (purple histograms) using the Kolmogorov–Smirnov test, and separation between distribution means was measured using Cohen’s d. (C) As in (B), but neuron pairs have additionally been sorted by the neurotransmitter identity of the presynaptic partner.

One-way ANOVATukey–Kramer post-hocAll Neurons0%2%p = 3 × 10−7***p = 1 × 10−17***
n = 760/678/4480%4%p = 5 × 10−18***2%4%p = 3 × 10−4***Acetylcholinergic0%2%p = 2 × 10−6***p = 2 × 10−17***
n = 335/315/2230%4%p = 9 × 10−19***2%4%p = 5 × 10−5***Glutamatergic0%2%p = 0.012*p = 3 × 10−4***
n = 235/213/1260%4%p = 4 × 10−4***2%4%p = 0.375GABAergic0%2%p = 0.131p = 0.006**
n = 46/36/260%4%p = 0.005**2%4%p = 0.336Dopaminergicp = 0.939
n = 52/43/28

ANOVA results, neuronal signal variance across depths of anesthesia for neuronal populations defined by neurotransmitter expression.

Difference of means across anesthesia conditions for each subgroup of neurons were compared using one-way ANOVA, and post-hoc pairwise comparisons were performed using the Tukey–Kramer method. See Figure 4A for visualization. *p < 0.05; **p < 0.01; ***p < 0.001.

Signal coherence of acetylcholinegeric and glutamaterigic synaptic pairs are most sensitive to isoflurane

While the suppression of neuronal activity appears to be a key feature of anesthesia, the anesthetized state is also characterized by a disruption in neuronal connectivity. Our previous work captured isoflurane induced general neuronal discoherence at the cellular level in C. elegans. We measured a reduction in the mean activity correlation between all pairs of imaged neurons within the nervous system as a global network metric (Awal et al., 2020), as well as suppression of other metrics such as mutual information and transfer entropy (Awal et al., 2018; Chang et al., 2023). Within our identified neuron data set we first measure activity correlations calculated using the Pearson correlation coefficient (PCC) across all identified synaptic pairs (i.e., the activity correlation between all pairs of identified neurons that are known to be connected via a direct synapse within the connectome). As might be expected, these synaptic pairs show greater correlation than the mean correlation across all pairs of imaged neurons without regard to synaptic connectivity: the mean correlation across all possible pairs of tracked neurons at 0% Isoflurane is 0.0766, which is substantially lower than that of identified synaptic pairs at 0.216. As shown in Figure 4B, with this correlation metric, we measure a mild decrease in the synaptic correlations between 0 and 2% (mean PCC, 0.216 to 0.192) and a larger decrease between 0 and 4% (mean PCC, 0.214 to 0.135) across all identified synapses. Shifts in correlation distributions are quantified by the mean correlation and tested for significance using the Kolmogorov–Smirnov test. These effects recapitulate our previous findings of decreased correlation across the ensemble of tracked but non-identified neurons in the animal’s head (Awal et al., 2020).

The identified dataset allows for further refinement based on individual synapse type. Applying the C. elegans chemical connectome to our identified neurons, we categorized sets of neuron pairs that are known to be linked via particular chemical synapses. Thus, the PCC can be calculated for any given subpopulation of synapses to determine how that select distribution is altered in the anesthetized state. For example, we were able to extract n = 559 neuron pairs that satisfied the following conditions: (1) putatively connected by a chemical synapse, (2) acetylcholinergic presynaptic partner, (3) both members of the neuron pair were identified in both the 0 and 2% isoflurane trials for that animal. Figure 4C displays the probability distribution histograms of neuron correlations for all ID synaptic pairs of the type and condition indicated. Acetylcholinergic and glutamatergic synaptic pairs at 2% isoflurane exhibit a significant suppression of PCC compared to baseline, with acetylcholinergic pairs decreasing 19% (mean PCC, 0.281 to 0.227) and glutamatergic synaptic pairs decreasing by 17% (mean PCC, 0.215 to 0.178). The effects become larger at 4% with an overall 45% decrease (mean PCC, 0.277 to 0.151) in correlation of acetylcholinergic pairs and a 46% decrease (mean PCC, 0.213 to 0.115) in glutamatergic synaptic pairs. We measured no such significance in the GABAergic pairs at 2%, although it should be noted the limited number of physical GABAergic synapses in the C. elegans nervous system (i.e., resulting in low n in the subclass) restrict the statistical power of this measurement. At 4%, GABAergic pairs show a significant increase in correlation (of 51%, mean PCC, 0.128 to 0.262) consistent with a loss of inhibitory/anti-correlated signaling. Dopaminergic pairs show no detectable change at 2% isoflurane and a smaller but significant decrease of 22% (mean PCC, 0.182 to 0.143) at 4%. These results suggest that acetylcholinergic and glutamatergic synaptic pairs experience disruption in synchrony at the earlier 2% isoflurane plane of anesthesia, and that they may drive the overall initial global dyssynchrony we observe across the entire system. By the time the moderate plane of 4% isoflurane is reached, the C. elegans nervous system experiences a more uniform disruption in signal coherence across all synapse types.

Neuron in-degree strongly correlates with isoflurane-induced neuronal activity suppression

To what extent do the functional effects of isoflurane depend on the total connectivity of individual neurons, i.e., their overall connectedness, or in- and out-degree? The out-degree of a neuron refers to the number of postsynaptic partners it has: i.e., how many neurons it sends information to. Conversely, the in-degree of a neuron refers to the number of presynaptic partners that neuron has: i.e., how many other neurons it receives information from. Figure 5 visualizes how out- and in