The 21st International Conference on Bioinformatics (InCoB2022): neuroscience
- Research
- Open access
- Published:
Computational modeling of light processing in the habenula and dorsal raphe based on laser ablation of functionally-defined cells
BMC Neuroscience volume 25, Article number: 22 (2024)
Abstract
Background
The habenula is a major regulator of serotonergic neurons in the dorsal raphe, and thus of brain state. The functional connectivity between these regions is incompletely characterized. Here, we use the ability of changes in irradiance to trigger reproducible changes in activity in the habenula and dorsal raphe of zebrafish larvae, combined with two-photon laser ablation of specific neurons, to establish causal relationships.
Results
Neurons in the habenula can show an excitatory response to the onset or offset of light, while neurons in the anterior dorsal raphe display an inhibitory response to light, as assessed by calcium imaging. The raphe response changed in a complex way following ablations in the dorsal habenula (dHb) and ventral habenula (vHb). After ablation of the ON cells in the vHb (V-ON), the raphe displayed no response to light. After ablation of the OFF cells in the vHb (V-OFF), the raphe displayed an excitatory response to darkness. After ablation of the ON cells in the dHb (D-ON), the raphe displayed an excitatory response to light. We sought to develop in silico models that could recapitulate the response of raphe neurons as a function of the ON and OFF cells of the habenula. Early attempts at mechanistic modeling using ordinary differential equation (ODE) failed to capture observed raphe responses accurately. However, a simple two-layer fully connected neural network (NN) model was successful at recapitulating the diversity of observed phenotypes with root-mean-squared error values ranging from 0.012 to 0.043. The NN model also estimated the raphe response to ablation of D-off cells, which can be verified via future experiments.
Conclusion
Lesioning specific cells in different regions of habenula led to qualitatively different responses to light in the dorsal raphe. A simple neural network is capable of mimicking experimental observations. This work illustrates the ability of computational modeling to integrate complex observations into a simple compact formalism for generating testable hypotheses, and for guiding the design of biological experiments.
Background
Light enables vision and also influences non-visual systems such as mood and alertness [1, 2]. The impact of light on mood and alertness is important for human well-being, but the underlying neurobiology is only partially understood [3]. Zebrafish are simple vertebrates whose transparent tissues are amenable to optical imaging. Their complex behaviour (e.g., swimming, hiding, freezing) provide a way for studying mood, anxiety, alertness, and other phenomena. The habenula is an evolutionarily conserved region of the brain that is of central importance in mood disorders, and is involved in processing responses to light. The habenula plays a pivotal role in motor and cognitive behaviours, by influencing the release of important neuromodulators such as serotonin, dopamine, epinephrine and histamine [4,5,6,7]. In clinical studies, deep brain stimulation of the habenula has been successful at triggering remission of therapy-resistant depression [8]. The role of the habenula in light processing is not fully understood, but light-evoked activity has been detected in the habenula of several different species, including the rat [9], pigeon [10] and zebrafish [11, 12].
Although the zebrafish brain is far simpler than mammalian brains, the zebrafish habenula has strong homology to the mammalian habenula. Just as the mammalian habenula contains two subdomains (medial and lateral), the zebrafish habenula can be divided into two subdomains - the dorsal (dHb) and ventral (vHb) habenula. In rat [9] and pigeon [10], electrophysiological recordings demonstrated that light triggers both excitation and inhibition in the habenula. Evoked activity, either phasic or sustained, was detected in both medial and lateral subdomains. Similar results have also been shown in zebrafish [12].
The dorsal raphe (DR) nucleus, a major structure downstream of the habenula, also participates in regulating brain state in response to lighting conditions as shown by recent studies in which larval zebrafish prefer different lighting conditions after optogenetic manipulation of the raphe [13] or habenula [12]. Tracing experiments have suggested that the medial habenula can modulate DR activity via the interpeduncular nucleus [14, 15], while the lateral habenula has direct and indirect projections [16]. The roles of these different pathways in mediating the effects of light is unknown. One way to confirm causal relationships between anatomical regions of the brain is to perform localized perturbations, such as via two-photon laser ablation to lesion specific cells [17, 18]. Two-photon laser ablation can provide a causal relationship because it allows ablation of functionally specific cells, which cannot be easily done with a genetically encoded tool. Nevertheless, establishing neuronal connectivity is rarely simple, and the results of localized perturbations can be extremely non-intuitive to interpret directly.
Mathematical modeling is used in many scientific studies where experimental evidence provides non-obvious clues about underlying phenomena: the model provides an interpretation of the evidence (a candidate mechanism), while computation provides unbiased evaluation (across millions or billions of candidates) for whether the candidate interpretation is indeed consistent with the observations. Here, we investigated raphe response to light (i.e., calcium response of neurons in the dorsal raphe), immediately after functionally specific (ON or OFF) cells of the habenula were ablated. Specifically, we developed a two-layer feed forward neural network model, that used the habenula response to light as input to simulate the raphe response. The goal of this modeling is to provide a compact computational replica of the experimental findings, because summarizing and recapitulating a complex system can give insight into causal patterns and hypotheses for future experimental testing.
Results
The habenula displays a broad and dynamic response to irradiance change
The zebrafish habenula consists of neurons surrounding neuropils that are innervated by afferent neurons [19,20,21,22]. We first characterized habenula activity evoked by pulses of light. Two-photon imaging was performed on a transgenic zebrafish line expressing the calcium indicator GCaMP3 throughout the habenula [23] (Fig. 1A). Resonant-scanning, combined with piezo-driven focusing, was used to record the activity of cells at multiple focal planes (Fig. 1B, C). With a step size of 10 μm, so that each cell would be sampled only once, most of habenula could be covered with 5 planes at a rate of 1 Hz. Habenula activity was monitored as the larva was exposed to 20- second pulses of blue light. We used relatively long pulses, rather than brief flashes, to allow longer characterization of responses. Steps involved in the analysis of habenula response to light are shown in Fig. 1D. In total, 2974 cells were identified in the habenula, and an analysis of their responses to light pulses showed the presence of multiple neuronal subtypes (Fig. 1E). We found neurons that showed excitatory response to light onset (the ON response) and also neurons that showed excitatory response to light offset (the OFF response). The finding here is consistent with previous findings using a short (1 s) red light pulse [11] and a 10-second blue light pulse [12].
We used the temporal light responses to group dorsal habenula (dHb) and ventral habenula (vHb) cells into individual clusters. The 2974 cells were first classified into dHb (1747 cells) or vHb (1227) region, depending on physical location. Then, k-means clustering was performed for each region (see Methods). Resulting clusters consisted of cells that were excited by light (ON cells), cells that were excited by subsequent darkness followed by sharply spiked response (phasic OFF cells), or cells excited by darkness followed by gradual decrease (tonic OFF cells). Both ON and OFF responses were seen in both dHB and vHb (Fig. 1F). This suggests that light OFF (or darkness) causes specific events in specific habenula cells, distinct from the ON cells. Thus, depending on their light response, the six habenula clusters identified were named D-ON-Tonic, D-OFF-Tonic, D-OFF-Phasic, V-ON-Tonic, V-OFF-Tonic, and V-OFF-Phasic. We refer to these six as the habenula subtypes.
Raphe response is modulated differently by ablating different regions of the habenula
Given the dynamic response in the habenula as shown in Fig. 1 and the persistence of dorsal raphe (DR) inhibitory response by light ON, as shown previously [13], we further tested if there is any functional link between the habenula light response and the DR inhibitory response by light ON. To test this, we used the technique of laser ablation or “lesioning” [12, 24], which allows localized disruption. In general, lesioning led to a sustained increase in GCaMP6f fluorescence throughout the targeted cell (Fig. 2A). Damage to fibers of passage was controlled for by imaging: namely, we only used lesions that did not lead to transient or sustained GCaMP6f fluorescence increase in regions away from the targeted site (e.g. in the contra-lateral habenula or pallium). There was no obvious change in the dHb response when cells in the ventral regions of the habenula were targeted, indicating that fibers of passage, which innervate the dorsal habenula, are not damaged. The workflow for our analysis of raphe response to light can be found in Fig. 2C.
In concordance to results from our earlier studies, we found again that with an intact habenula, the DR showed an inhibitory response during light exposure (black curves, Fig. 2D-F). When the laser was targeted to the left dorsal habenula neuropil regions that are activated by light ON, there was now light-evoked excitation rather than inhibition in the DR (Fig. 2D). Lesioning ON cells in the ventral region of the habenula led to a reduction in tonic activity in the DR, and few raphe cells showed any response to illumination change (Fig. 2E). When OFF cells were targeted in the ventral region, there was again a reduction in tonic activity of serotonergic neurons, accompanied by strong excitation upon transition to darkness (Fig. 2F). These lesion studies show that the distinct functional populations of the habenula, in response to light onset and offset, dynamically regulate response in DR serotonergic neurons. However, the nature of the habenula influence on the raphe was not clear from our experiments, and so we next developed a computational model.
Model construction
The goal of the modeling is to recapitulate the behavior of raphe (the output) in response to the habenula (the input) when different parts of the habenula are ablated or not. We curated the input to the modeling as follows. For each of the six habenula subtypes (Fig. 1F), we first created characteristic responses over a 100 s time frame that included a single 20s light exposure window (see Methods and Fig. 3A). After smoothing these curves, we used linear interpolation to obtain values at a resolution of 0.01s (the experimental measurements were taken at 1 s intervals). This resulted in a 9901-long vector for the neuronal response for each of the six habenula subtypes, for each of the four experimental conditions (3 ablations and 1 non-ablation condition). When a particular habenula subtype was ablated, its neuronal response was set to zero for that experimental condition. The 9901-long vectors of neuronal response were gathered for all six habenula subtypes and all four ablation or non-ablation conditions into a 39,604⨯6 matrix. Each column corresponds to one habenula subtype, and each row corresponds to the habenula response from time t = 1s to t = 100s in increments of 0.01s, concatenated across four experimental conditions. This habenula matrix is the input for our efforts to model how the raphe responds to the habenula (Fig. 3B).
For the raphe region of the brain, the same processing steps (averaging, smoothing, interpolation) were performed to obtain a raphe response profile under each of the four experimental conditions: unablated, D-ON ablated, V-ON ablated and V-OFF ablated. This resulted in a vector of size 39,604⨯1 for the raphe output dataset (the target output for our modeling, Fig. 3C).
We first tried to obtain mechanistic models of habenula-raphe interactions by developing multiple Ordinary Differential Equation (ODE) models that employed incoherent feedforward loops. The ODE models used biologically-informed networks of neuronal activation and inhibition, similar to those in Additional File 2. However, our ODE modelling efforts failed to recapitulate raphe responses from habenula inputs.
We then proceeded to implement simple multi-layer perceptron (MLP) feed forward neural networks (NN). Such networks would have an input layer with six nodes (the six habenula subtypes) and an output layer with one node (raphe response). Between these layers would be an intermediate hidden layer that would assimilate the habenula inputs and pass a signal to the raphe. In the simplest form, this hidden layer would have a single node and this design of the neural network also failed to recapitulate raphe responses. As the next simplest design, we tried a model with two nodes in the hidden layer, and this simple model was able to fit the observed raphe responses to a satisfactory degree, as estimated by the root mean square error (RMSE) between the experimentally observed and model-estimated Raphe behaviors. The weights and biases of the best fit model can be seen as arrow labels on Fig. 4A. The fits of model-estimated vs. observed raphe responses and the RMSE values are shown in Fig. 4B. The RMSE values for the unablated case, the V-ON ablated case and the V-OFF ablated case were all low and comparable (0.014, 0.015 and 0.012 respectively). The RMSE for the D-ON ablated experiment was comparatively higher at 0.0435, but still included a neuronal activity peak during the light ON period and a second peak once light was switched OFF, as observed in the experiments. The overall RMSE for all four ablation experiments was found to be 0.0246.
Because this is a simple network, the edge weights can also be used to draw inferences about neuronal influences over other neurons in the network, hinting towards the beginnings of mechanistic understanding. The weights suggest that the D-ON-Tonic cells in the habenula have an inhibitory effect on both intermediate nodes, while the V-ON-Tonic cells have an activating effect. It is important to note that since the corresponding neurons in dorsal and ventral regions (ON-Tonic, OFF-Tonic and OFF-Phasic) have mostly similar (but not identical) behaviors to light (Fig. 2), it is possible to get similar (but not identical) raphe behavior in the unablated experiments by assigning the weights of the Dorsal nodes to the ventral and vice versa. However, this would not fit the observations from the D-ON and V-ON ablation experiments. It is interesting that the two ON-Tonic nodes have opposing effects on the intermediate nodes. In contrast both D-OFF-Phasic and V-OFF-Phasic nodes have only activating effects on both intermediate nodes. Also, while the V-OFF-Tonic node has an inhibitory effect on both intermediate nodes, the D-OFF-Tonic node has an activating effect on one intermediate node and an inhibitory effect on the other intermediate node.
Finally, our model was also used to estimate the potential raphe response to a hypothetical D-OFF ablation experiment (Fig. 4C). Such an experiment in the future may be used to validate/ fine tune the parameters of the model. To test the robustness of our model we also performed ten-fold cross validation on our dataset. Briefly, for each ablation experiment, we split the raphe datasets into ten non-overlapping sets of neurons. Each set was taken to be the validation set for one iteration of the cross validation, while the other nine were taken as the training set. The average raphe responses were computed for the training and validation datasets. An MLP was fit to the training set and the resulting model (the “best-fit MLP”) was evaluated using the RMSE relative to the training set and relative to the validation set. This was repeated for ten iterations. The range of training and test errors over the ten-fold cross validation can be found in Additional File 3.
Interpreting the NN model parameters
To seek mechanistic understanding from our best fit neural network (NN) model, we computed the contribution of each node in our model, to nodes in the successive layer (Fig. 5A-D). According to our model, in the unablated case, the inhibitory phenotype shown in the raphe response to light activation is mostly driven by an inhibitory effect of the H1 node on the output raphe node. This can be seen by the negative values of H1→R in the leftmost column of Fig. 5A between 20 and 40s (period when light is on). It is also important to note that the weight of the edge between H1 and R in our model is negative, and so strong activators of the H1 node will cause strong inhibition of the raphe node. Looking at the contribution of individual habenula subtype nodes to H1, we see that H1 activation during light is mostly due to strong activation by the V-ON-Tonic cells (middle column, Fig. 5A). In contrast, the recovery of this inhibitory phenotype after light is switched off, is due to activation of the raphe node by the H2 node. This can be seen from positive values of H2→R between time 40–50 s in the leftmost column of Fig. 5A.
In the D-ON ablated case, we see that raphe activation during light activation occurs as a competition between inhibiting H1 and activating H2 nodes, resulting in overall activation of the raphe node (leftmost column Fig. 5B). As earlier, this is a consequence of V-ON-T activation of the H1 and H2 nodes. In addition, the ablation of D-ON results in the removal of D-ON inhibition of H2, which also contributes to activation of the R node by H2 (as opposed to no effect in the unablated case). In the V-ON ablated case, the lack of strong activation of H1 and H2 nodes (seen earlier) is now absent, and hence is responsible for an overall loss of signal from H1 and H2 to the R node (Fig. 5C). Finally in the V-OFF ablated case, the strong phasic response after light offset is driven mostly by H2 activation of the R node, as evidenced by the positive values of H2→R between time 40–50 s in the leftmost column of Fig. 5D. In addition, ablation of the V-OFF-Phasic node results in removal of the inhibitory effect due to V-OFF phasic activation of H1 node in the unablated case.
Model sensitivity analysis
To obtain the sensitivity of model estimation to changes in each parameter of the NN model, we individually scaled each NN parameter by a factor between 0.1 and 10 (0.1x– 10x) and compared the model simulations against experimental observations of raphe behavior. For each ablation experiment, the error is compared as root mean square deviation after the simulated raphe response has been translated to have a median response of 1 over the first 10s (out of the 100s window). As seen in Fig. 4C, the NN model was more sensitive to changes in some parameters and less sensitive to others. For example, the model is relatively less sensitive to changes in the weight of the edges from D-ON-Tonic to the H1 and H2 nodes. In contrast, the model is much more sensitive to changes in the edges from V-ON-Tonic to H2, which is unsurprising, considering the magnitude of the contributions of V-ON-Tonic to raphe responses, seen in Fig. 5A-D. Also, as a general trend, higher error values are found when parameters are scaled up rather than scaled down. Furthermore, the model is less sensitive to changes in edges going into the H1 node than into the H2 node.
Discussion
Typically, ablation of upstream neuronal regions causes downstream regions to lose responsiveness to external stimuli. Remarkably, in our experiments, lesioning different types of functionally-specific neurons in the habenula caused the downstream neuronal responses in the dorsal raphe to change in several unexpected ways– from inhibitory to excitatory, from light-excitatory to dark-excitatory, and from tonic to phasic. Our work used brief exposures to light and dark in larval zebrafish, because light and dark are strong stimuli for these brain regions and therefore aid dissection of neural connectivity. Another reason we used light and dark is because they are evolutionarily connected to predators and danger response. In humans, danger response is linked to psychiatric conditions such as anxiety and depression, and the human habenula is known to have an influential role in anxiety and depression.
Previous work showed that neural activity in the dorsal raphe is generally inhibited by light [13], and we also found that dorsal raphe neural activity (measured by calcium fluorescence intensity change) was negatively correlated with light exposure. We then identified different populations of habenula cells (ON and OFF cells) that were positively or negatively correlated with light. We used laser ablation to disable specific types of neurons in the habenula, specifically the ON-cells or specifically the OFF- cells. The raphe response to light and dark, after each type of ablation, was altered in remarkable ways which contain intricate information about how the dorsal raphe is affected by different types of neurons in the habenula.
To characterize and analyze this unexpected diversity of responses, we analyzed neuronal imaging and constructed a simple multilayer perceptron model that was able to recapitulate raphe responses from habenula input. We chose an MLP model for our estimations because compared to other statistical algorithms for prediction, a simple MLP model allowed both high regression accuracy and the ability to draw insights from the model parameters. Our computational simulations using the MLP model mimic the biological observations by producing four different states of the DR: (1) inhibited by light; (2) spontaneously active regardless of irradiance; (3) excited by light and (4) excited by darkness. While we do not claim that the model in Fig. 4A is the only model that is capable of recapitulating raphe responses from habenula inputs, it could be one among a set of similar equally-likely models. Such a set of equally-likely models would include uncertainty about the parameter values as well as mechanistic uncertainty. These variations cannot be disambiguated with the available data, but we can confidently assert that a small number of connections (as in Fig. 4A) is sufficient to explain the flipped responses seen in the ablation experiments seen in Fig. 2. The problem of choosing among possible models depends on whether the goal is fitting or prediction. Our primary goal is to describe the data rather than to generalize to different cases, but our model has parameters that might predict raphe behavior in unknown cases, so it illustrates one plausible scenario with a logical explanation for interpretation of the experimental results. When additional experiments are performed in the future, the resulting data can further improve the modeling. While we did perform ten-fold cross validation with our model, there were challenges with such validation. Since the input and output data are representations of average behavior of habenula and raphe, there is no strict one-to-one correspondence between the input and output instances. For true one-to-one correspondence, the habenula and the raphe responses would need to be measured in the same fish, which is not feasible with our ablation experiments. Furthermore, since each individual row or “instance” in the input and output corresponds to habenula/raphe behavior at time intervals separated by just 0.01s, it is very likely that the input and output values of any randomly chosen validation set would have very close neighbors (in value) in the training set. This defeats the purpose of a random validation set, as these values are not independent of each other. Hence, we chose to perform 10⨯ cross-validation by computing average behavior from a subset of raphe neurons, while using the average of the left-out neurons as the validation (Additional File 3).
Upon deeper inspection of the weights and the contributions of each node in our neural network to the final raphe response, we observed that most of the interesting phenotypes seen with the raphe response to light are dominated by the activity of V-ON-Tonic nodes and its activation of the H1 and H2 nodes, with smaller contributions from the other habenula subtypes. If our modelling estimations are accurate, that would also imply the presence of biological equivalents of the H1 and H2 nodes, which act as intermediates between the habenula and raphe. While previous studies have suggested the rostromedial tegmental nucleus as one such mediator [16], future experiments can help identify the anatomical region and the precise set of neurons involved in this mediation.
In summary, this study used computational modeling to delineate specific hypotheses about the interplay between the dorsal and ventral habenula during raphe response to light and dark, with systems-level correlations emerging from analysis of the best fit model. Computational modeling is particularly valuable in the field of neuroscience because even a single lesion in the brain can yield rich and non-intuitive phenotypes. Our model reveals that a simple network can create a variety of complex behaviors, and can explain the seemingly paradoxical observations observed with ablation experiments. Other aspects of the model, such as the hypothesized nature of each functional connection and the relative strength of each effect are more speculative, and their usefulness is more likely to lie in the holistic illustration of a plausible scenario rather than a literal specification of calcium regulation in the neurons. Qualitative modeling can shed light on confusing phenomena and guide the choice of validation experiments, while providing testable hypotheses for future experiments.
Conclusion
Because the habenula has a causal influence on anxiety/depression disorders, there is urgent need to map functional connectivity from the habenula to the raphe at the level of neurons or neuron clusters. When neurons in specific regions of the habenula were ablated, neurons in the dorsal raphe showed reproducible patterns of change in how they responded to light and dark. Response patterns exhibited a surprising variety of qualitative changes, depending on which habenula site was ablated. Computational modeling was undertaken because simulating this variety could aid the design of future experiments and the disambiguation of potential mechanisms. An artificial neural network was trained by statistical machine learning and was successful at recapitulating the observed range of response profiles. In other words, the surprising variety we observed could be captured using a simple, compact formalism. Future work should continue to upgrade models of the habenula until both the model mimics the microanatomy in not just its outputs but also its internal formalisms. Inevitable future growth in data availability and computational resources will expand modelling research toward comprehensive integration of neuroscience experimentation and simulation. Such integration will be necessary for the overarching biomedical objective, which is mapping the functional connectivity of zebrafish neuronal structures that are homologous to structures implicated in human psychiatric disorders.
Methods
Experimental methods
Fish lines
Experiments were performed in accordance with guidelines issued by the Institutional Animal Care and Use Committee (IACUC) of the Biological Resource Centre at Biopolis, Singapore. Zebrafish (Danio rerio) lines used for this study were: Tg(tph2:GAL4, UAS: Kaede)y228 [25], GAL4s1011t [26] and Tg(elavl3:GCaMP6f)a12200.
Two-photon calcium imaging
Zebrafish larvae (aged 5–10 days-post-fertilization, dpf) were immobilized in mivacurium (1.5 mg/ml) and embedded in low-melting temperature agarose (2.0% in E3) in a glass-bottom dish (Mat Tek). They were imaged on a Nikon two-photon microscope (A1RMP), using a 25x water immersion objective (NA = 1.1). The femtosecond laser (Coherent Vision II) was tuned to 920 nm for GCaMP6f imaging. Stacks were collected in resonant-scanning mode with 2x pixel averaging. The sample size was based on [11]. Blue light stimulus was generated by 5 mm blue LEDs (458 nm peak emission), which was powered by a 5 V TTL signal from a control computer and synchronized with image capture using a National Instruments DAQ board, controlled by the Nikon Elements software. Each light pulse was 20 s long and followed by 20 s dark with a total of 4 pulses of light. Light intensity at the sample was 0.13 mW/cm2.
Laser ablation
Larval zebrafish (6–8 days post fertilization, dpf) were subjected to laser ablation for ON or OFF cells in the habenula, as described in [27]. Briefly, Tg(Elavl3:GCaMP6f) larvae were anaesthetized and then mounted in 2% low melting temperature agarose. After identifying habenula ON and OFF cells, lesions were conducted via several pulses (100–500 msec) of the femto-second laser (960 nm). Lesioning was monitored by time-lapse imaging before and after each pulse, and was terminated when there was a localized increase in GCaMP6f fluorescence. Animals were discarded if lesioning caused bursting of blood vessels in the habenula.
Image analysis
Initial data pre-processing
All image processing steps were performed following the same protocol as our earlier work [27]. In short, raw images obtained were first registered to correct for any vertical/horizontal movement artefacts using cross-correlation. Then, a median filter of size 3 was applied to remove noise. A darker region outside the region of interest was chosen as the background and subtracted from the image to remove background noise. Non-linear trends in the data were detrended using polynomials of order 2–5. Data was then normalized into Z-scores by subtracting the overall mean and dividing by the standard deviation. A rolling window average was then used to smooth noisy traces where necessary. Where possible, cells were segmented (see below) or images were directly analysed as pixels (using the Thunder platform [28] for fast pixel-based clustering and factorization).
Cell segmentation
Each stack was scaled 2x in imageJ (RRID: SCR_003070), then maximally projected to a single image, which was then subjected to a minimum filter and unsharp mask to sharpen the boundary of cells. ROIs were identified using the “find maxima…” command, as a way to localize regional darkest point as the center of each ROI. The boundary of the ROI was outlined by “analyze particle…” that connects bright pixels into mosaic-like tessellated plane, encircling each darkest point. Each ROI was then numbered sequentially using the ImageJ ROI Manager and mapped back to the original despeckled image stack. Manual segmentation was done here to delete extraneous ROIs outside the habenula or dorsal raphe and to encircle cells that were not detected by the algorithm (< 10% of total ROIs). In the last step, “Set measurements…” and “measure” in ImageJ provided the mean fluorescence value of all pixels within each ROI across the entire image stack and the x-y coordinates of each ROI. Time-lapse series in which z drifting occurred were excluded, as in this case ROIs could not be defined.
Statistical and modeling methods
Classification of the habenula cells
The habenula was segmented and then, to identify evoked responses to light, the spatiotemporal calcium dataset was subjected to k-means clustering as described in [27]. Habenula cells were classified into being dorsal or ventral depending on their location (position in the Z-stack), and k-means clustering was performed individually for cells in each region, as in our earlier studies [13, 27], to find cells that responded similarly to the light pulses. The clustering method was k-means and the distance measure was correlation, and this configuration ran quickly and was able to recapitulate well-established classifications of neuronal phenotype such as tonic and phasic. Not knowing a priori the ideal number of neuronal clusters, we ran k-means clustering with different k values from 2 to 10 and compared results. For each region (dorsal/ventral), we manually inspected the clustering results and judged the optimal k to be three (see Additional File 1 for clustering results for k values 2–5).
Defining model input and output
For each of the six habenula subtypes identified via k-means clustering, we first obtained a mean response for a 40 s window that comprised of a 20 s light period followed by a 20s dark period. This response was obtained by averaging across four subsequent light/dark cycles for each cell, and then computing the mean response of all cells in each cluster for this window. Although at the end of the 40 s window a new light cycle was initiated in the experiments, we observed that some phenotypes had not yet decayed to the baseline. Hence to simulate complete response for each habenula subtype, we first padded the response on each side as follows. For the time before light exposure, we appended a 20s period of baseline activity, simulating unexcited neurons. For padding after the 40s window, we used PCHIP (Piecewise Cubic Hermite Interpolating Polynomial) interpolation to extend the response for 40 more seconds, bringing the total to 100s for each habenula cluster (20s unexcited + 20 s light exposure + 60 s decay to baseline). These 100s habenula responses were then smoothed to remove noise and were used as input to our model.
Similarly for the dorsal raphe, we constructed a mean 100s behavior for raphe response to each ablation experiment as follows. Since the raphe response across subsequent light/dark cycles were not consistent for each cell unlike the habenula, we used only the response to the first light exposure to obtain mean 40s response for each experiment. These were then padded and smoothed similar to the habenula data to obtain 100s responses for the raphe. This was used as the “true” output against which our model to fit its simulated responses.
Neural network model
To simulate raphe responses from habenula input, we used a simple multilayer perceptron (MLP) neural network (NN). The network had 6 input nodes (corresponding to each habenula subtype), a hidden layer with two nodes and a third layer with one output node. All nodes in the hidden and output layer used a Rectified Linear Unit (ReLU) activation function. Any other hyperparameters were either set to default values (e.g., learning rate), or were not included (e.g., dropout). Earlier attempts with ordinary differential equation models employing incoherent feedforward loops, or an MLP with a one-node hidden layer proved insufficient to recapitulate raphe behavior from habenula input. The final NN model was implemented and fitted using the Keras (v2.2.4-tf) and TensorFlow (v2.1.0) APIs in Python 3.7.6.
Data availability
The images obtained from Zebrafish experiments can be found at the Figshare database (DOI: https://doi.org/10.6084/m9.figshare.21721058). The codes used in analysis can be found on github at the URL https://github.com/nsuhasj/habenulaRaphe/. Any other datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- DAQ:
-
Data Acquisition
- dHb:
-
dorsal habenula
- D–OFF:
-
Dorsal OFF cells
- D–ON:
-
Dorsal ON cells
- DR:
-
Dorsal Raphe
- LED:
-
Light emitting diode
- MLP:
-
Multi–layer perceptron
- NA:
-
Numerical aperture
- NN:
-
Neural Network
- ODE:
-
Ordinary Differential Equation
- PCHIP:
-
Piecewise Cubic Hermite Interpolating Polynomial
- ReLU:
-
Rectified Linear Unit
- RMSE:
-
Root mean square error
- ROI:
-
Region of Interest
- TTL:
-
Transistor–Transistor Logic
- vHb:
-
ventral habenula
- V:
-
OFF–Ventral OFF cells
- V:
-
ON–Ventral ON cells
References
Vandewalle G, Schwartz S, Grandjean D, Wuillaume C, Balteau E, Degueldre C, et al. Spectral quality of light modulates emotional brain responses in humans. Proc Natl Acad Sci U S A. 2010;107(45):19549–54.
Badia P, Myers B, Boecker M, Culpepper J, Harsh JR. Bright light effects on body temperature, alertness, EEG and behavior. Physiol Behav. 1991;50(3):583–8.
Chellappa SL, Gordijn MCM, Cajochen C. Can light make us bright? Effects of light on cognition and sleep. Prog Brain Res. 2011;190:119–33.
Jhou TC, Fields HL, Baxter MG, Saper CB, Holland PC. The rostromedial tegmental nucleus (RMTg), a GABAergic afferent to midbrain dopamine neurons, encodes aversive stimuli and inhibits motor responses. Neuron. 2009;61(5):786–800.
Quina LA, Tempest L, Ng L, Harris JA, Ferguson S, Jhou TC, et al. Efferent pathways of the mouse lateral habenula. J Comp Neurol. 2015;523(1):32–60.
Morley BJ, Spangler KM, Javel E. The development of somatostatin immunoreactivity in the interpeduncular nucleus of the cat. Brain Res. 1985;352(2):241–8.
Wang RY, Aghajanian GK. Physiological evidence for habenula as major link between forebrain and midbrain raphe. Science. 1977;197(4298):89–91.
Sartorius A, Kiening KL, Kirsch P, von Gall CC, Haberkorn U, Unterberg AW, et al. Remission of major depression under deep brain stimulation of the lateral habenula in a therapy-refractory patient. Biol Psychiatry. 2010;67(2):e9–11.
Zhao H, Rusak B. Circadian firing-rate rhythms and light responses of rat habenular nucleus neurons in vivo and in vitro. Neuroscience. 2005;132(2):519–28.
Semm P, Demaine C. Electrophysiology of the pigeon’s habenular nuclei: evidence for pineal connections and input from the visual system. Brain Res Bull. 1984;12(1):115–21.
Dreosti E, Vendrell Llopis N, Carl M, Yaksi E, Wilson SW. Left-right asymmetry is required for the habenulae to respond to both visual and olfactory stimuli. Curr Biol. 2014;24(4):440–5.
Zhang BB, Yao YY, Zhang HF, Kawakami K, Du JL. Left Habenula mediates light-preference behavior in zebrafish via an asymmetrical visual pathway. Neuron. 2017;93(4):914–e9284.
Cheng RK, Krishnan S, Jesuthasan S. Activation and inhibition of tph2 serotonergic neurons operate in tandem to influence larval zebrafish preference for light over darkness. Sci Rep. 2016;6:20788.
Lima LB, Bueno D, Leite F, Souza S, Gonçalves L, Furigo IC, et al. Afferent and efferent connections of the interpeduncular nucleus with special reference to circuits involving the habenula and raphe nuclei. J Comp Neurol. 2017;525(10):2411–42.
Agetsuma M, Aizawa H, Aoki T, Nakayama R, Takahoko M, Goto M, et al. The habenula is crucial for experience-dependent modification of fear responses in zebrafish. Nat Neurosci. 2010;13(11):1354–6.
Sego C, Gonçalves L, Lima L, Furigo IC, Donato J, Metzger M. Lateral habenula and the rostromedial tegmental nucleus innervate neurochemically distinct subdivisions of the dorsal raphe nucleus in the rat. J Comp Neurol. 2014;522(7):1454–84.
Cheng Z, Han Y, Wei B, Lai B, Li B, Cui M, et al. Probing neuronal functions with precise and targeted laser ablation in the living cortex. Optica. 2021;8(12):1559.
Muto A, Kawakami K. Ablation of a neuronal Population using a two-photon laser and its Assessment using calcium imaging and behavioral Recording in zebrafish larvae. J Vis Exp. 2018;(136).
Hendricks M, Jesuthasan S. Asymmetric innervation of the habenula in zebrafish. J Comp Neurol. 2007;502(4):611–9.
Miyasaka N, Morimoto K, Tsubokawa T, Higashijima S, ichi, Okamoto H, Yoshihara Y. From the olfactory bulb to higher brain centers: genetic visualization of secondary olfactory pathways in zebrafish. J Neurosci. 2009;29(15):4756–67.
Amo R, Fredes F, Kinoshita M, Aoki R, Aizawa H, Agetsuma M, et al. The habenulo-raphe serotonergic circuit encodes an aversive expectation value essential for adaptive active avoidance of danger. Neuron. 2014;84(5):1034–48.
Turner KJ, Hawkins TA, Yáñez J, Anadón R, Wilson SW, Folgueira M. Afferent connectivity of the zebrafish Habenulae. Front Neural Circuits. 2016;10:30.
Krishnan S, Mathuru AS, Kibat C, Rahman M, Lupton CE, Stewart J, et al. The right dorsal habenula limits attraction to an odor in zebrafish. Curr Biol. 2014;24(11):1167–75.
Chung SH, Mazur E. Surgical applications of femtosecond lasers. J Biophotonics. 2009;2(10):557–72.
Yokogawa T, Hannan MC, Burgess HA. The dorsal raphe modulates sensory responsiveness during arousal in zebrafish. J Neurosci. 2012;32(43):15205–15.
Scott EK, Baier H. The cellular architecture of the larval zebrafish tectum, as revealed by gal4 enhancer trap lines. Front Neural Circuits. 2009;3:13.
Cheng RK, Krishnan S, Lin Q, Kibat C, Jesuthasan S. Characterization of a thalamic nucleus mediating habenula responses to changes in ambient illumination. BMC Biol. 2017;15(1):104.
Freeman J, Vladimirov N, Kawashima T, Mu Y, Sofroniew NJ, Bennett DV, et al. Mapping brain activity at scale with cluster computing. Nat Methods. 2014;11(9):941–50.
Acknowledgements
The authors gratefully acknowledge Guo Yuge and Seetha Lakshmi Krishnan for data analysis ideas and helpful discussions.
Funding
This research is supported by the Singapore Ministry of Education Academic Research Fund Tier 2 grant MOE2019-T2-1-138 and by Duke-NUS SRP Phase 2 Research Block Grant to LTK; core funding from the Institute of Cell and Molecular Biology to SJ; and an NTU undergraduate student fellowship to AIK. The Duke-NUS block grant paid the publication fees. Any opinions, findings and conclusions or recommendations expressed in this material are those of the author(s) and do not reflect the views of the funders. Funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Author information
Authors and Affiliations
Contributions
RKC and SJ designed the zebrafish experiments. RKC and AIK conducted the experiments. NSJ, AIK and LTK designed and performed the statistical analysis and modeling. AIK, NSJ, RKC, LTK and SJ wrote the manuscript. All authors have read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Experiments were carried out in accordance with protocols approved by the Institutional Animal Care and Use Committee.
Consent for publication
Not applicable.
Conflict of interest
The authors claim no conflicts of interests.
About this supplement
This article has been published as part of BMC Neuroscience Volume 25 Supplement 1, 2024: The 21st International Conference on Bioinformatics (InCoB2022): neuroscience. The full contents of the supplement are available online at https://bmcneurosci.biomedcentral.com/articles/supplements/volume-25-supplement-1.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Cheng, RK., Jagannathan, N.S., Kathrada, A.I. et al. Computational modeling of light processing in the habenula and dorsal raphe based on laser ablation of functionally-defined cells. BMC Neurosci 25 (Suppl 1), 22 (2024). https://doi.org/10.1186/s12868-024-00866-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12868-024-00866-z