
==== Front
bioRxiv
BIORXIV
bioRxiv
2692-8205
Cold Spring Harbor Laboratory

39229236
10.1101/2024.08.16.607981
preprint
1
Article
Whole-brain mapping in adult zebrafish and identification of a novel tank test functional connectome
http://orcid.org/0000-0003-0565-0945
Rajput Neha
Parikh Kush
Squires Ada
http://orcid.org/0000-0001-5080-390X
Fields Kailyn K.
Wong Matheu
Kanani Dea
http://orcid.org/0000-0001-8790-5184
Kenney Justin W.
Department of Biological Sciences, Wayne State University, Detroit, MI 48202
19 8 2024
2024.08.16.607981https://creativecommons.org/licenses/by/4.0/ This work is licensed under a Creative Commons Attribution 4.0 International License, which allows reusers to distribute, remix, adapt, and build upon the material in any medium or format, so long as attribution is given to the creator. The license allows for commercial use.
nihpp-2024.08.16.607981.pdf
Identifying general principles of brain function requires the study of structure-function relationships in a variety of species. Zebrafish have recently gained prominence as a model organism in neuroscience, yielding important insights into vertebrate brain function. Although methods have been developed for mapping neural activity in larval animals, we lack similar techniques for adult zebrafish that have the advantage of a fully developed neuroanatomy and larger behavioral repertoire. Here, we describe a pipeline built around open-source tools for whole-brain activity mapping in freely swimming adult zebrafish. Our pipeline combines recent advances in histology, microscopy, and machine learning to capture cfos activity across the entirety of the adult brain. Images captured using light-sheet microscopy are registered to the recently created adult zebrafish brain atlas (AZBA) for automated segmentation using advanced normalization tools (ANTs). We used our pipeline to measure brain activity after zebrafish were subject to the novel tank test. We found that cfos levels peaked 15 minutes following behavior and that several regions containing serotoninergic, dopaminergic, noradrenergic, and cholinergic neurons were active during exploration. Finally, we generated a novel tank test functional connectome. Functional network analysis revealed that several regions of the medial ventral telencephalon form a cohesive sub-network during exploration. We also found that the anterior portion of the parvocellular preoptic nucleus (PPa) serves as a key connection between the ventral telencephalon and many other parts of the brain. Taken together, our work enables whole-brain activity mapping in adult zebrafish for the first time while providing insight into neural basis for the novel tank test.
==== Body
pmcIntroduction

A fundamental goal of neuroscience is to understand how patterns of brain activity give rise to behavior. Identifying general principles of brain function is facilitated by cross species comparisons. Over the past two decades, zebrafish have started contributing to our understanding of the brain, a trend that promises to continue due to their low cost, ease of genetic manipulation, and sophisticated behavioral repertoire (Gerlai, 2023; Kenney, 2020; Loring et al., 2020). Although several methods have been developed for whole-brain activity mapping in larval zebrafish (Ahrens et al., 2012; Portugues et al., 2014; Randlett et al., 2015; Shainer et al., 2023), equivalent approaches have yet to be developed for adult stage animals.

Adult and larval zebrafish each have distinct advantages and disadvantages in the study of brain-behavior relationships. Whereas larval animals are amenable to high throughput work due to their small size and transparency, adults have the advantage of mature neuroanatomy and more extensive behavioral repertoire. This behavioral repertoire includes a wide variety of social behaviors (Gerlai, 2014; Jones and Norton, 2015; Kareklas et al., 2023), short and long-term associative, non-associative, and spatial memories (Gerlai, 2020; Kenney, 2020), and different types of exploratory behaviors (Cachat et al., 2010; Rajput et al., 2022; Toms and Echevarria, 2014). Thus, to fully realize the utility of zebrafish as a model organism in neuroscience, methods for whole-brain mapping are also required for adult zebrafish.

Whole-brain activity mapping can yield unexpected insights into brain function that may be lost using more targeted methods. Measuring neural activity across the entire brain also facilitates the use of powerful analytic tools, like network analysis, that captures complex interactions and improves predictions of brain-behavior relationships (Vetere et al., 2017; Wheeler et al., 2013). However, mapping whole-brain activity presents several technical challenges. One roadblock is that the brain of adult animals is not transparent, and thus requires the use of tissue clearing (Richardson et al., 2021). Imaging intact organs presents another technical hurdle due to the increased volume, a challenge met by the recent development of light-sheet microscopy (Hillman et al., 2019). Finally, whole-brain mapping results in large amounts of data that cannot be analyzed via traditional approaches like manual counting and segmentation. We tackled this challenge by combining advances in machine learning to automate cell detection (Tyson et al., 2021) and image registration (Gholipour et al., 2007) with the recently created digital adult zebrafish brain atlas (AZBA) (Kenney et al., 2021). Here, we describe how we have assembled these tools into a pipeline that enables whole-brain activity mapping in adult zebrafish for the first time.

Results

Overview of strategy

We begin by giving an overview of our strategy for whole-brain activity mapping (figure 1) before describing the results of each step in more detail. Following behavior, animals are euthanized and heads fixed in 4% paraformaldehyde overnight. Following careful dissection, brains are rendered optically transparent using iDISCO+ (Renier et al., 2016), which we modified to make it compatible with in situ hybridization chain reaction (HCR) for the detection of cfos mRNA (Choi et al., 2018; Kramer et al., 2018; Kumar et al., 2021). Imaging intact cleared brain tissue was done using light-sheet microscopy. To automatically identify cfos positive cells in the brain, we used the open source CellFinder package (Tyson et al., 2021) that is part of the BrainGlobe suite of Python-based software tools (Claudi et al., 2020). Finally, to automatically parcellate the brain into individual regions, we used advanced normalization tools (ANTs; Avants et al., 2009)) to register images to AZBA (Kenney et al., 2021). The final output of our pipeline is a list of cfos positive cell counts for each brain region and each animal. This enables the use of a variety of downstream analytic tools, one example that we demonstrate here is functional network analysis.

Automated cell detection

After in situ HCR, tissue was cleared using iDISCO+, which allowed us to use light-sheet microscopy to capture whole-brain images in both the cfos (Figure 2A, top) and autofluorescence channels (Figure 2A, bottom). Detection of cfos positive cells was done using CellFinder (Tyson et al., 2021), an artificial neural net-based supervised machine learning algorithm. The first step in the cell detection process uses image filtering to detect cell shaped objects in the cfos image. We found parameters that captured cfos positive cells throughout the entire brain (described in the methods section), including areas with cells of different sizes and densities like the telencephalon (Figure 2B) and cerebellum (Figure 2C). Because the cell detection algorithm generated a lot of overlapping cells, we used a custom written Python script to remove cell candidates that were within 9 μm of one another. We then trained the CellFinder artificial neural network by manually labelling 10,597 cells and 7,303 non-cells across five brains. Non-cells were unambiguously identified by the presence of a signal in both the cfos and autofluorescence channels, suggesting the presence of background bleeding into the cfos channel. Cells only appeared in the cfos channel. The resulting network achieved over 95% accuracy where the cells and non-cells were clearly differentiated across several different brain regions (Figure 2B & C).

During imaging, we noticed that we had sufficient resolution to differentiate cytoplasmic and nuclear localization of cfos. Nuclear staining was characterized by the presence of puncta whereas cytoplasmic staining had a conspicuous dark spot surrounded by more diffuse fluorescence (Figure 2D). This localization of cfos is an indication of how long ago the cell was active as the mRNA is first transcribed in the nucleus before being shuttled to the cytoplasm for translation. To capture this distinct cellular localization, we created and trained an artificial neural net on 2,448 examples of nuclear puncta and 1,916 examples of cytoplasmic staining to differentiate these different patterns of cfos staining. This network also achieved greater than 95% accuracy.

Registration to the adult zebrafish brain atlas

The adult zebrafish brain contains over 200 regions, making manual segmentation implausible. To automate parcellation of brains into individual regions, we used ANTs (Avants et al., 2009) to register brains to AZBA using common autofluorescence images. Initially, we attempted to register the autofluorescence image in AZBA directly to individual autofluorescence images, but the results were inconsistent (data not shown). We had more success by first making an average template by registering together 10 autofluorescence images from present study (Figure 3A). The autofluorescence image from AZBA was then successfully registered to this template brain (Figure 3B). A handful of small anomalies arose from this registration process that we manually fixed using ITK-SNAP (Yushkevich et al., 2019). These arose in parts of the image that tend to be highly variable between individuals, such as where mounting occurs at the ventral hypothalamus and the dorsal sac that extends from the dorsal diencephalon. To segment individual brains, we used the transforms from registering the template autofluorescence brain to individual images (Figure 3C). Using inverse transformations from the registration process, we were also able to bring cfos images into the space of AZBA (Figure 3D).

Time course for cfos expression

To effectively map whole-brain activity we need to know at what point after behavior cfos expression peaks. We exposed fish to a commonly used behavioral task, the novel tank test, and euthanized animals 5, 15, 30, 60, or 120 minutes following the behavior (Figure 4). We also had two control groups: (1) fish that were euthanized immediately after removal from their housing racks, and (2) fish that were brought into the behavioral room and euthanized an hour later, mimicking the habituation to the behavioral room we use for fish that were exposed to the novel tank (i.e., time = 0). A sex × time ANOVA found a large effect of time (P < 0.001, η2 = 0.54), a trend towards a small effect of sex (P = 0.07, η2 = 0.059), and no interaction (P = 0.46). Using a Dunnet’s t-test to compare all groups to the home tank (HT) control group, we found a large increase in cfos cell density at 15 minutes (P = 0.00067, d = 2.07) when cfos activity peaked (Figure 4A & B).

We also examined how the proportion of nuclear and cytoplasmic stained cells changed across time (Figure 4C). A cell type × time ANOVA found a large main effect of time (P < 0.001, η2 = 0.40) and no overall effect of cell type (P = 0.95). There was also a large interaction between cell type and time (P = 0.0082, η2 = 0.13). FDR corrected paired t-tests at each time point found that there were more nuclear than cytoplasmic stained cells at 5 minutes (P = 0.048). This trend switched to more cytoplasmic than nuclear stained cells at 15 and 30 minutes, although the differences at these time points were not statistically significant (P’s = 0.16 & 0.22, respectively).

Cell types active during the novel tank test

AZBA contains several stains that can be used to identify different cell types across brain regions such as 5-hydroxytryptamine (5-HT), tyrosine hydroxylase (TH), and choline acetyltransferase (ChAT) (Kenney et al., 2021). To determine if exposure to a novel tank results in the activation of regions containing these neuronal cell types, we brought home tank and 15-minute cfos brains into the same space as AZBA, averaged the images together, and looked for overlap between the stains in AZBA and elevated cfos (Figure 5). For regions expressing 5-HT (Figure 5A), we saw an increase in cfos in the paraventricular organ (PVO), intermediate nucleus (IN), and caudal zone of the periventricular hypothalamus (Hc). For TH, which labels dopaminergic and noradrenergic cells, we saw overlap in the ventromedial thalamic nucleus (VM), the posterior part of the parvocellular preoptic nucleus (PPp), paracommissural nucleus (PCN), and Hc (Figure 5B). Finally, for ChAT, we saw overlap in the paraventricular gray zone of the optic tectum (PGZ; Figure 5C). Although we can see overlap at the regional level, our findings are only tentative because the cfos and antibody stained images come from separate brains, so we cannot make claims at the cellular level. Nonetheless, this demonstrates how our approach can be used to generate hypotheses about roles different neurotransmitters may play in the underlying a behavior.

Brain network analysis

We used functional network analysis to gain insight into the organization of brain activity that underlies exploration of a novel tank (Pinho et al., 2023; Vetere et al., 2017; Wheeler et al., 2013). Using cfos counts from the 15-minute time point, we computed the correlated activity between all 143 gray matter regions across animals (Figure 6). To filter the correlation matrix to generate a network, we used efficiency cost optimization where the network density is chosen such that it balances the inclusion of edges to increase global and local efficiency against the putative cost of including additional connections (Fallani et al., 2017). We found a density of 2.5% maximized the efficiency cost optimization quality function (Figure 7A). This resulted in a network with 256 edges and an average degree of 3.6, which is consistent with other functional brain networks generated using different imaging modalities (Fallani et al., 2017). This network also exhibited small world properties: its average shortest path length between nodes was 5.6, which is similar to the average path length of the average from equivalently dense random networks (3.9) with much higher clustering (0.38 versus 0.024). This yielded a small world coefficient greater than 1 (11.0) indicating the expected small world property (Humphries and Gurney, 2008). We also computed degree and eigenvector centrality for each node to uncover brain regions that may play outsized roles in the network (Figure 7C). This uncovered four regions that were in the top 10 for each of these centrality measures: the ventral nucleus of the ventral telencephalon (Vv), the dorsal zone of the ventral telencephalon (Vd-dd), the dorsal most zone of the ventral telencephalon (Vdd), and the anterior part of the parvocellular preoptic nucleus (PPa).

Next, we used the Louvain algorithm (Blondel et al., 2008) to identify 10 distinct communities in the network (Figure 7B). Using the network and community structure, we categorized the roles that different nodes play in interconnecting different parts of the network (Guimerà and Amaral, 2005): provincial hubs (highly connected within its community, but not between communities), connector hubs (highly connected both within and between communities), peripheral nodes (low connectivity within and between communities), and non-hub connectors (low connectivity within a community, but high between communities). Interestingly, the PPa, which was identified as important based on centrality measures, arises as a connector hub. The PPa interconnects a module dominated by regions of the ventral telencephalon with other parts of the preoptic area (SC and PPp), thalamus (VM, CP, and ZL) and hypothalamus (ATN, Hv, Hc, and Hd). Thus, our network analysis points to the PPa and ventral telencephalon as likely playing an important role in regulating behavior during exploration of a novel tank.

Discussion

In the present study, we introduce a pipeline for performing whole-brain activity mapping in adult zebrafish. Our pipeline combines several recently developed tools: a digital brain atlas for adult zebrafish (Kenney et al., 2021), registration using ANTs (Avants et al., 2011), machine learning tools for automated cell detection (Tyson et al., 2021), tissue clearing (Renier et al., 2014), light-sheet microscopy (Reynaud et al., 2014), and in situ HCR (Choi et al., 2018) for detecting cfos. Importantly, all the computational tools are open access and free to use. Furthermore, to aid in the implementation of this pipeline, we have included a bench protocol (Supplemental file 1). The primary stumbling blocks for implementing this pipeline are likely to be access to a light-sheet microscope for whole-brain imaging and sufficient computational power for training and applying the registration and CellFinder machine learning algorithms. The former issue is partly mitigated by the increased availability of light-sheet microscopes, particularly in core facilities. Access to computational resources can be addressed by using tools like Google Colaboratory (Bisong, 2019) or high performance computing facilities available at many institutions.

Cfos to capture whole-brain activity

We captured neural activity using in situ HCR to detect cfos mRNA. We chose this approach for several reasons: (1) there are a paucity of antibodies for detecting cfos protein in zebrafish, none of which are known to work in whole-mount tissue-cleared samples, (2) in situ HCR probes are small (~150 bp), which easily penetrates chunks of intact tissue like the adult zebrafish brain, and (3) cfos is one of the most widely used markers of neural activity due to autoinhibition of transcription that results in low background, high signal-to-noise, and good temporal resolution (Chung, 2015; Lucibello et al., 1989). The findings in the present study further support these rationales: we saw even penetration of cfos staining throughout the brain (Figures 2 and 4B) and the levels of background cfos staining were low, with an approximately 3.5 fold increase in cfos density 15 minutes following behavior compared to quiescent animals removed directly from their housing racks (Figure 4A). The increase in cfos was also tightly coupled to the behavior, peaking 15 minutes after exposure to the novel tank before decreasing to baseline levels by 60 minutes. Interestingly, if we look at only cells that have nuclear staining, we see the increase begins as soon as 5 minutes after behavior. The higher cfos density at 15 minutes is likely due to the opportunity for increased transcription which would be expected to create a brighter signal resulting a larger number of detectable cells. The time to maximal cfos we observed is faster than is seen in rodents, where it is often found to peak at 30 minutes post-stimulation (Ding et al., 1994; Guzowski et al., 2001; Kovács, 1998; Zangenehpour and Chaudhuri, 2002). The reasons for this time difference between zebrafish and rodents is unclear. Nonetheless, it emphasizes the importance of performing time course analysis when establishing new methods for brain mapping in different species.

Other markers of neural activity have gained traction in recent years in zebrafish, such as the phosphorylated forms of ribosomal protein S6 (pS6) and extracellular regulated kinase 1/2 (pERK1/2). Our data suggests that cfos as an activity marker compares favorably to these options. For pS6, the signal-to-noise ratio is comparable to what we see for cfos, with an approximately 2–4 fold increase over baseline both in vivo in zebrafish (Butler et al., 2018; Parada et al., 2024; Scaia et al., 2022) and in vitro neuronal cell culture (Kenney et al., 2015). However, the time course of elevated pS6 is notably slower, taking an hour or more to peak (Kenney et al., 2015; Parada et al., 2024) compared to 15 minutes for cfos (Figure 4). In contrast, pERK1/2 activity peaks quickly, within 2–5 minutes, but the signal-to-noise ratio is ~0.5–1, considerably lower than cfos (Randlett et al., 2015; Venincasa et al., 2021). This low signal-to-noise ratio likely arises from higher background levels of pERK due to the wide variety of cellular processes that it regulates (Cargnello and Roux, 2011). Thus, the best choice of stain depends on the behavioral paradigm. Large, rapid responses to brief behavioral stimuli are best captured by pERK. However, more subtle responses may be missed due to the low signal-to-noise ratio. S6 phosphorylation excels at capturing long lasting steady-state neural activity, as suggested by Maruska et al (2020) and would excel for behaviors lasting 30 minutes or more. Cfos represents a solid middle ground that is ideal for capturing neural activity from behaviors lasting on the order of 5–10 minutes, like the novel tank test used in the present study.

Registration to AZBA to identify cell types

We were able to successfully register our brains to AZBA using ANTs (Avants et al., 2009). To do so, we first used ANTs to make an average template from our images by registering 9 brains to a single brain and averaging them together. The autofluorescence image in ABZA was then registered to this average template, yielding good results (Figure 3). We chose this method because we found that registering the autofluorescence image from AZBA to individual brains gave inconsistent results. This is likely because the autofluorescence image in AZBA is also an average of many brains (Kenney et al., 2021). We chose ANTs because the non-linear symmetric diffeomorphic image registration it employs has been consistently found to be one of the best algorithms for 3D image registration (Klein et al., 2009; Murphy et al., 2011). The tool is also well documented and straightforward to use. Finally, ANTs has recently grown in popularity for image registration in larval zebrafish (Marquart et al., 2017; Shainer et al., 2023), which provided a starting point for identifying the best parameters for registration in our samples.

Following registration to AZBA, we were able to identify potential neuronal cell types relevant to the novel tank test (Figure 5). We found that several regions containing high levels of 5-HT were active during behavior, such as the PVO, IN, and Hc. Consistent with this, several papers have implicated 5-HT as contributing to exploration of a novel tank using pharmacological approaches (Beigloo et al., 2024; Maximino et al., 2013; Nowicki et al., 2014; Wong et al., 2010). Similarly, there was overlap in cfos activity in several regions that express tyrosine hydroxylase (VM, PPp, PCN, and Hc), implicating these populations of dopaminergic or noradrenergic neurons in novel tank behavior (Kacprzak et al., 2017; Nabinger et al., 2023). Of the cfos positive cells that overlap with TH, our network analysis suggests that the PPp may be of particular importance in regulating exploratory behavior, as it is one of the few non-hub connectors (Figure 7D). The PPp also has a direct connection to the PPa region, which ranks highly in both eigenvector and degree centrality (Figure 7B), and connects to the thalamic VM region, another area high in TH expression. This suggests that the PPp and VM may may act in concert to mediate the effects of the dopaminergic system on exploration. However, one important caveat to these interpretations is that we are comparing averaged cfos images to averaged neurotransmitter-related stains in AZBA, and thus we cannot definitively identify the specific cell types that are active. This would require co-staining of brains with both cfos and various neuronal cell-type markers to determine if the activity of these specific cell types changes.

Novel tank functional connectome

Using our whole brain mapping data, we generated the first novel tank functional connectome. The novel tank test is one of the most widely used behavioral tests in adult zebrafish, often used to study exploratory and anxiety-related behaviors (Blaser et al., 2010; Kalueff et al., 2013; Luca and Gerlai, 2012; Rajput et al., 2022; Spence et al., 2006). Our functional network analysis identifies several key regions that are engaged during exploration of a novel tank for the first time (Figures 6 and 7). In particular, the medial portion of the ventral telencephalon stands out, where several subregions (the Vv, Vd-dd, Vc, Vd-vd, and Vp) rank highly on at least one measure of centrality (Figure 7C). These regions are also highly interconnected, a fact that is clear from both the correlation matrix (Figure 6) and the community they form in the network (dark orange in Figure 7B). Based on molecular markers, these regions of the ventral telencephalon are thought to correspond to the mammalian subpallial amygdala (i.e., the central and medial amygdala) and basal ganglia (Mueller, 2022; Porter and Mueller, 2020). In mammals, these brain regions have been found to be important for a wide range of behaviors, from defensive, anxiety-related, and social behaviors to motor control (Fadok et al., 2018; Grillner and Robertson, 2016; Raam and Hong, 2021). Our findings that the ventral telencephalon appears to be engaged during the novel tank test is reasonable given that novelty and exploration would be expected to engage circuits involved in decision making, emotional regulation, and muscle coordination.

In examining how the regions of the ventral telencephalon interact with the rest of the brain, a few interesting trends emerge. Notably, the interaction of ventral telencephalic regions with many other communities is anti-correlated (i.e., the dark green, light orange, and grey communities in Figure 7B). This suggests the presence of strong inhibitory connections between the medial ventral telencephalon and other parts of the brain. Consistent with this interpretation, the ventral telencephalon has been found to contain a substantial number of inhibitory GABAergic neurons (Porter and Mueller, 2020). Our network analysis suggests that these inhibitory connections are most likely present between the ventral telencephalon and the Vmn (mesencephalic nucleus of the trigeminal nerve), End (entopeduncular nucleus in the lateral portion of the ventral telencephalon), and from the BSTm (bed nucleus of the stria terminalis, medial portion in the dorsal telencephalon) to R (rostrolateral nucleus in the thalamus). However, given that our findings are correlational in nature, techniques like tract tracing and direct manipulation would be needed to confirm these interactions.

Our network analysis also identified the PPa as a region of high importance. The PPa was high in both eigenvector and degree centrality (Figure 7C) and was one of the few connector hub nodes (Figure 7D). In examining its place in the network (Figure 7B), the PPa interconnects with several regions of the ventral telencephalon and, working in concert with the PPp, mediates their interactions with parts of the network that contain several thalamic and hypothalamic regions (magenta cluster in Figure 7B). To our knowledge, the correspondence between the PPa and PPp in teleosts and tetrapods has not been determined. Based on the expression of neuropeptides like oxytocin and arginine vasopressin, parts of the PPp are thought to be equivalent to the supraoptic nucleus in mammals (Herget et al., 2014). In larval zebrafish, the preoptic area has recently been implicated in behaviors such as navigation, thermoregulation, and stress reactivity (Corradi et al., 2022; Palieri et al., 2024). However, the preoptic area in larval zebrafish cannot be differentiated into subregions like the PPa and PPp due to a lack of cytoarchitectural boundaries (Herget et al., 2014). This makes it unclear as to what specific regions in the adult would subsume the functions identified in larval animals. Future work should determine the role that these different subregions might play in different aspects of exploration and anxiety-like behavior in adult zebrafish.

Summary

The present study provides an open-source framework for performing whole-brain mapping in adult zebrafish. This work also yielded the first description of brain activity that underlies the novel tank test, suggesting the medial ventral telencephalon may play an important role in one of the most widely used behavioral tasks in adult zebrafish. Taken together, we anticipate that our pipeline will help generate insights into the principles of brain function by enhancing the utility of adult zebrafish as a model organism.

Methods

Animals

Zebrafish

Subjects were 8–10 month old zebrafish of the TU strain from both sexes. Fish were bred and raised at Wayne State University and within two generations of animals obtained from the Zebrafish International Resource Center (ZIRC, catalog ID: ZL84) at the University of Oregon. Fish were maintained in high-density racks under standard conditions: water temperature of 27.5 ± 0.5 °C, salinity of 500 ± 10 μS, and pH of 7.4 ± 0.2. Lighting followed a 14:10 light:dark cycle, with lights on at 8:00 AM. Fish were fed twice daily with dry feed (Gemma 300, Skretting, Westbrook, ME, USA) in the morning and brine shrimp (Artemia salina, Brine Shrimp Direct, Ogden, UT, USA) in the afternoon.

Sex determination was based on secondary sex characteristics such as shape, color, and the presence of pectoral fin tubercles (McMillan et al., 2015). Confirmation was conducted post-experimentation by euthanizing the animals and observing the presence or absence of eggs. All experimental procedures were conducted under the ethical approval of the Wayne State University Institutional Animal Care and Use Committee (Protocol ID: 21-02-3238).

Behavioral stimuli and tissue collection

The novel tank test was used as the behavioral stimulus, using tanks that were distinct from housing tanks. Behavioral tanks were open top five-sided (15 × 15 × 15 cm) and made from frosted acrylic (TAP Plastics, Stockton, CA, USA). Each tank was filled to a height of 12 cm with 2.5 L of fish facility water and housed within a white corrugated plastic enclosure to minimize external disturbances and diffuse light.

One week before the novel tank test, animals were housed in 2-liter tanks divided into two chambers with transparent dividers. Male and female pairs were kept in each chamber to enable identification of individuals without social isolation or tagging. A day prior to the experiment, animals were acclimatized to the behavior room for one hour before being placed back on the housing racks. On the day of the experiment, animals were removed from the housing rack and allowed to acclimate in the behavioral room for one hour. After acclimation, animals were individually transferred to a novel tank and allowed to explore the tank for 6 minutes. Water was replaced between animals. After six minutes, fish were removed and placed back in their home tank for a designated periods of time (5, 15, 30, 60, or 120 minutes) prior to euthanization. A subset of animals was euthanized one hour after acclimation to the room (home tank control) and another set of animals were euthanized immediately after removal from the housing racks (rack control).

Animals were euthanized by immersion in ice cold water for 5 minutes to induce anesthesia and then decapitated using a sharp blade. Heads were then washed in ice-cold phosphate buffered saline (PBS) for 60 seconds to allow for blood drainage, and then fixed in 4% paraformaldehyde in PBS overnight. Brains were then dissected in ice cold PBS and subject to iDISCO+ and in situ HCR.

Histology

Tissue pre-treatment

We adapted the iDISCO+ protocol (Renier et al., 2016) for zebrafish brain tissue staining. Following dissection, brain samples were washed for 30 minutes, three times, in PBS at room temperature. This was followed by dehydration using a methanol concentration gradient (20, 40, 60, 80, and 100%) for 30 min each. Samples were further washed in 100% methanol, chilled on ice, and then incubated in chilled 5% hydrogen peroxide in methanol overnight at 4°C. The next day, the samples were rehydrated through a reverse methanol series (80%, 60%, 40%, 20%) at room temperature, followed by a 1 h PBS wash, two 1 h PBS-T washes (1x PBS, 0.1% Tween 20), and a 3 h PBS-T wash. Samples were then equilibrated overnight in 5× SSCT (sodium chloride sodium citrate/0.1% Tween-20) buffer.

In-situ HCR

We modified the original HCR method described by Choi and colleagues (2018) and informed by the work of Kumar et al (2021). Samples were first prepared by acetylation in 0.25% v/v acetic anhydride solution in ultrapure water for 30 min. Samples were then washed in ultrapure water three times for 5 mins and then equilibrated in probe hybridization buffer (30% formamide, 5x SSC, 9 mM citric acid, 0.1% Tween-20, 50 μg/mL heparin, 1x Denhardt’s solution, 10% Dextran sulfate) for 15 min at room temperature. Samples were then incubated in probe hybridization buffer for 1h at 37 °C with shaking and then incubated with 1 μM of cfos probes in hybridization buffer at 37 °C with shaking for 48–60 hours. Samples were then washed with probe wash buffer (30% formamide, 5x SSCT, 9 mM citric acid, 50 μg/mL heparin) three times at 37°C, then twice with 5× SSCT for 1 h each with shaking. The tissue was then equilibrated in amplification buffer (5x SSC, 0.1% Tween-20, 10% Dextran sulfate) at room temperature for 1h with shaking. Alexa647 labeled hairpins (B1) were prepared by heating to 95 °C for 90 seconds prior to cooling at room temperature in the dark. We diluted 7.5 pmol of each hairpin into 125 μL of amplification buffer for each sample. Samples were incubated for 48–60 hours in the dark at room temperature. Finally, tissue was washed in 5× SSCT, 5 times for 1 hour each before being washed overnight in 5X SSCT.

Tissue clearing

Samples were first dehydrated in a series of methanol-water mixtures (20%, 40%, 60%, 80%, 100%) at room temperature for 1 hour each and then left in 100% methanol overnight. The next day, samples were incubated at room temperature in a mixture of 66% dichloromethane and 33% methanol for 3 hours followed by two 15-minute washes in dichloromethane. After removing the dichloromethane, samples were incubated and stored in dibenzyl ether at room temperature for at least 24 hours until imaging.

Brain imaging and processing

Cleared samples were imaged on a LaVision BioTec UltraMicroscope II (Miltenyi Biotec, Auburn, CA) using Imspector software for image acquisition. The microscope setup included a 4.2 Megapixel sCMOS camera and a 2x objective lens with a dipping cap with spherical aberration correction. Images were taken at a magnification of 6.4x. Samples were mounted on the sample holder using an ultraviolet cured resin (NOA 61, Norland Products, Jamesburg, NJ) with a refractive index (1.56) that matched DBE. Imaging was conducted from the right laser sheet with a 4 μm step size using dynamic horizontal focus. Both 480 nm autofluorescence and 640 nm signal channels were used. The imaging settings used were: 90% laser power, 200-ms exposure time, 50% sheet width, sheet numerical aperture of XX. Acquired images were stitched using Terastitcher (Bria and Iannello, 2012).

Computational analysis

Automated cell detection

For the automated detection and quantification of cfos positive cells, we utilized the Python-based software, CellFinder (Tyson et al., 2021). It comprises two steps: cell candidate detection and cell classification. The initial step of cell detection identifies cell-like objects in the image. We optimized parameters to capture as many cell-like objects in our images as possible. Running from the Linux terminal, we used the following command for cell detection:

cellfinder -s path/to/folder/signal/channel/cfos -b /path/to/folder/background/channel/AF -o path/to/output1 -v 3.990 0.943 0.943 --orientation sal --no-register --no-classification -- soma-diameter 5 --threshold 3 --ball-xy-size 2 --ball-z-size 7 --ball- overlap-fraction 0.77 --log-sigma-size 0.1 --save-csv --batch-size 64 --epochs 100

After detecting cell candidates, a customized python script was used to remove cell candidates that were within 9 μm of one another.

Napari was utilized for visualization and labelling. We manually annotated 10,597 cells and 7,303 non-cells across five brains for training the artificial neural network. Cellfinder was trained using the following command:

Cellfinder_train -y path/to/brain1_labels.yml path/to/brain2_labels.yml path/to/brain3_labels.yml path/to/brain4_labels.yml path/to/brain5_labels.yml -o /trained_network --batch-size 64 --epochs 100 --no-save-checkpoints -- save-progress

The trained network achieved 96.1% accuracy. Finally, the trained network was applied to all the experimental brains to classify the detected cell candidates into cells and non-cells. This was achieved by utilizing the following command:

cellfinder -s /path/to/folder/signal/channel/cfos/ -b /path/to/folder/background/channel/AF/ -o path/to/output -v 3.990 0.943 0.943 --orientation sal --no-register --no-detection --soma- diameter 5 --threshold 3 --ball-xy-size 2 --ball-z-size 7 --ball- overlap-fraction 0.77 --log-sigma-size 0.1 --save-csv --trained-model /trained_network/model.h5

Differentiating nuclear and cytoplasmic stained cells

To differentiate between cytoplasmic and nuclear puncta, we developed a convolutional neural network (CNN) built in Python using the TensorFlow library. The architecture of the CNN is outlined in Table S2. Cfos Images from 10 brains were labelled, totaling 2,448 puncta and 1,916 cytoplasmic labels. A training dataset was created by isolating 11ˣ11ˣ11 pixel cubes centered around each of the labeled cells. The dataset was split 80/20 into a training set and a testing set. The input data was augmented through a series of horizontal and vertical flips, 90° rotations, and 2-pixel horizontal translations to create a total training dataset of 13,706 puncta and 10,724 cytoplasmic labels. No data augmentation was performed on the testing set. The model was trained using an NVIDIA GeForce 3090 GPU for 500 epochs. The batch size was 32, the weight decay rate was 0.0005, and the learning rate was 0.0001. The model achieved an accuracy of 95.3% on the testing set.

Brain registration

Image registration was performed using ANTs (Avants et al., 2009). For the non-linear diffeomorphic step, four parameters were optimized: cross-correlation, gradient step, update field variance in voxel space, and total field variance in voxel space to achieve the best alignment. Using the optimized parameters, brain registration was carried out in two steps: first, an average brain template was created, and second, AZBA was registered to this average template.

Before registration, images were downsampled to 4 μm isotropic using brainreg from the BrainGlobe suite of tools (Tyson et al., 2021):

brainreg /path/to/raw/data /path/to/output/directory -v 3.990 0.943 0.943 --orientation sal --atlas azba_zfish_4um –debug

The average template was generated using 10 autofluorescence images. Initially, nine autofluorescence images were individually brought into the space a single image (template) using the following ANTs command:

antsRegistration --dimensionality 3 --float 1 -o [${AF_sample_1_for_avg_},${ AF_sample_1_for_avg-warped}] -- interpolation WelchWindowedSinc -u 0 -r [${AF_template.nii},${AF_sample_1.nii},1] -t Rigid[0.1] -m MI[${AF_template.nii},${AF_sample_1.nii},1,32,Regular,0.25] -c [200 x 200 × 200 × 0,1e-8,10] --shrink-factors 12×8×4×2 --smoothing-sigmas 4×3×2×1vox -t Affine[0.1] -m MI[${AF_template.nii},${AF_sample_1.nii}, 1,32,Regular,0.25] -c [200 × 200 × 200 × 0,1e-8,10] --shrink-factors 12×8×4×2 --smoothing-sigmas 4×3×2×1vox -t SyN[0.3,4,0] -m CC[${AF_template.nii},${AF_sample_1.nii}, 1,3] -c [200 × 200 × 200 x 200, 1e-6,10] --shrink-factors 12×8×4×2 --smoothing-sigmas 4×3×2×1vox --verbose 1

These outputs were then used to create an average image using the ‘AverageImages’ command in ANTs. Next, the autofluorescence image from AZBA was registered to the average template using the following command:

antsRegistration --dimensionality 3 --float 1 -o [${AZBA_to_avg_temp_},${AZBA_to_avg_temp-warped}] --interpolation WelchWindowedSinc -u 0 -r [${avg_template.nii.gz},${AZBA/20180628_AF_average.nii.gz},1] -t Rigid[0.1] -m MI[${avg_template.nii.gz},${AZBA/20180628_AF_average.nii.gz},1,32,Regu lar,0.25] -c [200 × 200 × 200 × 0,1e-8,10] --shrink-factors 12×8×4×2 - -smoothing-sigmas 4×3×2×1vox -t Affine[0.1] -m MI[${avg_template.nii.gz},${AZBA/20180628_AF_average.nii.gz}, 1,32,Regular,0.25] -c [200 × 200 × 200 × 0,1e-8,10] --shrink-factors 12×8×4×2 --smoothing-sigmas 4×3×2×1vox -t SyN[0.3,4,0] -m CC[${avg_template.nii.gz},${AZBA/20180628_AF_average.nii.gz}, 1,3] -c [200 × 200 × 200 × 200, 1e-6,10] --shrink-factors 12×8×4×2 -- smoothing-sigmas 4×3×2×1vox --verbose 1

To bring the segmentation from AZBA into the space of the template we used the following command:

antsApplyTransforms -d 3 --float -n NearestNeighbor -i /AZBA/2021–08- 22_AZBA_segmentation.nii.gz -r avg_template.nii.gz -o AZBA_to_avg_temp_transformed.nii.gz -t AZBA_to_avg_temp_1Warp.nii.gz - t AZBA_to_avg_temp_0GenericAffine.mat

Finally, the newly generated average template image was used as a reference image and was registered onto individual autofluorescence images:

antsRegistration --dimensionality 3 --float 1 -o [${AF_sample_},${AF_sample-warped}] --interpolation WelchWindowedSinc -u 0 -r [${AF_sample.nii},${avg_template.nii.gz },1] -t Rigid[0.1] -m MI[${AF_sample.nii},${ avg_template.nii.gz },1,32,Regular,0.25] -c [200 × 200 × 200 × 0,1e-8,10] --shrink-factors 12×8×4×2 --smoothing- sigmas 4×3×2×1vox -t Affine[0.1] -m MI[${AF_sample.nii},${ avg_template.nii.gz },1,32,Regular,0.25] -c [200 × 200 × 200 × 0,1e- 8,10] --shrink-factors 12×8×4×2 --smoothing-sigmas 4×3×2×1vox -t SyN[0.3,4,0] -m CC[${AF_sample.nii},${ avg_template.nii.gz },1,3] -c [200 × 200 × 200 × 200, 1e-6,10] --shrink-factors 12×8×4×2 -- smoothing-sigmas 4×3×2×1vox --verbose 1

Finally, segmentation of individual brains was done using the same command as above but applied to the segmentation file as the floating image.

Cfos cell counts and network analysis

R (version 4.1.1) was used for network analysis and to integrate the output from Cellfinder with the brain segmentation using the RNifti package (Clayden et al., 2021) to read in the segmentation files. The number cfos positive cells in each brain were summed excluding white matter and clear labelled regions yielding 143 gray matter regions for analysis.

Network analysis was performed using the igraph (version 2.0.2) package (Csardi and Nepusz, 2006). The network was generated by treating the correlation matrix (Figure 6) as an adjacency matrix. For thresholding we chose the network density using efficiency cost optimization to maximize the quality function (Fallani et al., 2017): j=Eg+Elρ

Where Eg is the global efficiency, El is the average of the local efficiency, and ρ is the network density. For the calculations of global and local efficiency we used a binarized network based on the absolute value of the correlations.

For identifying node roles, we calculated the within module degree z-score: zi=κi+κ¯siσκsi

Where ki is the number of connections between node i and other nodes in the same community and κ¯si is the average of over all nodes in a community; σκsi is the standard deviation of the number of connections in a community. We also calculated the participation coefficient: Pi=1−∑s=1Nc(KisKi)2

Where Nc is the number of communities, Kis is the number of connections between node i and all other nodes in community s, and Ki is the degree of node i. The definitions of the above equations and the boundaries for the different types of nodes were based on the Guimerà and Amaral (2005).

The small worldness parameter was calculated as described in (Humphries and Gurney, 2008): σ=LgLrandCgCrand

Where Lg is the average shortest path length between all nodes of the network, Lrand is the average shortest path length between all nodes in an equivalent random network, Cg is the clustering coefficient of the network, and Crand is the clustering coefficient of an equivalent random network. For random network parameters, we took the average from 1,000 instances of Edros-Renyi random networks (Erdös and Rényi, 2011) with an equivalent number of nodes and edges as the target network.

Statistical analysis

Statistical analysis was done using R. Data were analyzed using 2 × 2 ANOVAs as indicated in the results. For the overall time course cfos data, Dunnet’s t-tests were used to compare all other groups to the home tank control group (time = 0). False discovery rate (FDR; Benjamini and Hochberg, 1995) corrected paired t-tests at each time point were used for cytoplasmic versus nuclear data.

Supplementary Material

Supplement 1

Supplement 2

Acknowledgments and Funding

We thank Jacob Hudock and Dinh Luong for excellent care of the zebrafish and facility maintenance. This work was supported by the National Institutes for General Medical Sciences (R35GM142566, JWK) and the Richard Barber Interdisciplinary Research Program (KKF, JWK).

Data and code availability

Data and code are available at github:

https://github.com/KenneyLab/RajputEtAl_2024_Whole_brain_mapping

Figure 1. Overview of method for mapping neural activity in adult zebrafish. Following behavior, zebrafish are euthanized and brains carefully removed. In situ HCR is then used to label cfos. Brians are then cleared using iDISCO and imaged using light-sheet microscopy. Cells are then detected using CellFinder and brains are registered to AZBA. Regional cfos counts are then used to generate brain networks for further analysis.

Figure 2. Staining for cfos and identifying cfos positive cells. A) Adult zebrafish brain stained for cfos (top) and the corresponding autofluorescence image (bottom). Scale bar is 0.5 mm. B & C) Zoomed in sections of the brain corresponding to red squares in part A showing cfos staining and autofluorescence with labelling of cells (yellow arrows) and non-cells (pink triangles). Scale bars are 0.1 mm. D) Examples of cfos staining in the cell nucleus and cytoplasm. Scale bar is 10 μm.

Figure 3. Registration of brain images to AZBA. A) Image of 10 brains registered and averaged. B) Segmentation from AZBA applied to average brain in A. C) Segmentation from AZBA applied to an individual zebrafish brain. D) An individual cfos brain brought into the space of AZBA. Scale bars are 0.5 mm.

Figure 4. Time course for cfos expression following exploration of a novel tank. A) Cfos cell density across the entire brain in animals taken off the rack, that remained in their home tank (HT), or a different times after exploration (5, 15, 30, 60, or 120 minutes). * - p < 0.05 compared to the HT group. B) Cfos stained brains from each time point were brought into the space of AZBA and averaged and displayed in the coronal plane. The numbers on the left of image are the distance (in mm) from the anterior most portion of the brain. Scale bar is 0.5 mm. C) Number of cfos cells classified as nuclear or cytoplasmic at each time point. * - p < 0.05 difference between the number of nuclear and cytoplasmic cells at that time point. Sample sizers were as follows: rack: female: n=4, male: n=4; HT: female: n=5, male: n=4; 5 min: female: n = 6, male: n = 5; 15 min: female: n=7, male: n=6; 30 min: female: n=5, male: n=5; 60 min: female: n=4, male: n=4; 120 min: female: n=5, male: n=3.

Figure 5. Overlap between cfos expression and neurotransmitter-related stains in AZBA. Regional overlap for A) 5-HT, B) TH, and C) ChAT. Scale bar is 0.5mm. Numbers on left are distance from anterior most portion of the brain in mm.

Figure 6. Correlation matrix of cfos activity across the zebrafish brain. Entries in the matrix are Pearson correlations between brain regions across animals euthanized 15 minutes after the novel tank test. Regions are organized based on common ontological levels. Regional abbreviations and ontological levels can be found in Table S1.

Figure 7. Analysis of the brain network active during the novel tank test. A) Efficiency-cost optimization for different network densities. J: quality function (see methods), El: local efficiency, EG: global efficiency. B) Network filtered at a density of 2.5%. Connections between nodes represent suprathreshold correlations from Figure 6. Color of connections represents the strength (darker means higher absolute value) and direction (red: positive, blue: negative) of the correlation. Node colors correspond to communities. Regions not in the giant component are not shown. C) Degree and eigenvector centrality for the top 10 brain regions. Gray bars are those regions that are in the top 10 for both degree and eigenvector centrality. D) Identification of the role that each node plays in the network based on within module degree z-score and participation coefficient.
==== Refs
References

Ahrens MB , Li JM , Orger MB , Robson DN , Schier AF , Engert F , Portugues R . 2012. Brain-wide neuronal dynamics during motor adaptation in zebrafish. Nature 485 :471–477. doi:10.1038/nature11057 22622571
Avants B , Tustison NJ , Song G . 2009. Advanced Normalization Tools: V1.0. The Insight Journal. doi:10.54294/uvnhin
Avants BB , Tustison NJ , Song G , Cook PA , Klein A , Gee JC . 2011. A reproducible evaluation of ANTs similarity metric performance in brain image registration. NeuroImage 54 :2033–2044. doi:10.1016/j.neuroimage.2010.09.025 20851191
Beigloo F , Davidson CJ , Gjonaj J , Perrine SA , Kenney JW . 2024. Individual differences in the boldness of female zebrafish are associated with alterations in serotonin function. Journal of Experimental Biology jeb.247483 . doi:10.1242/jeb.247483
Benjamini Y , Hochberg Y . 1995. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society: Series B (Methodological) 57 :289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
Bisong E. 2019. Google Colaboratory In: Bisong E , editor. Building Machine Learning and Deep Learning Models on Google Cloud Platform: A Comprehensive Guide for Beginners. Berkeley, CA: Apress. pp. 59–64. doi:10.1007/978-1-4842-4470-8_7
Blaser RE , Chadwick L , McGinnis GC . 2010. Behavioral measures of anxiety in zebrafish (Danio rerio). Behavioural Brain Research 208 :56–62. doi:10.1016/j.bbr.2009.11.009 19896505
Blondel VD , Guillaume J-L , Lambiotte R , Lefebvre E . 2008. Fast unfolding of communities in large networks. Journal of statistical mechanics: theory and experiment 2008 :P10008.
Bria A , Iannello G . 2012. TeraStitcher - A tool for fast automatic 3D-stitching of teravoxel-sized microscopy images. BMC Bioinformatics 13 :316. doi:10.1186/1471-2105-13-316 23181553
Butler JM , Whitlow SM , Roberts DA , Maruska KP . 2018. Neural and behavioural correlates of repeated social defeat. Sci Rep 8 :6818. doi:10.1038/s41598-018-25160-x 29717159
Cachat J , Stewart A , Grossman L , Gaikwad S , Kadri F , Chung KM , Wu N , Wong K , Roy S , Suciu C , Goodspeed J , Elegante M , Bartels B , Elkhayat S , Tien D , Tan J , Denmark A , Gilder T , Kyzar E , DiLeo J , Frank K , Chang K , Utterback E , Hart P , Kalueff AV . 2010. Measuring behavioral and endocrine responses to novelty stress in adult zebrafish. Nat Protoc 5 :1786–1799. doi:10.1038/nprot.2010.140 21030954
Cargnello M , Roux PP . 2011. Activation and Function of the MAPKs and Their Substrates, the MAPK-Activated Protein Kinases. Microbiology and Molecular Biology Reviews 75 :50–83. doi:10.1128/mmbr.00031-10 21372320
Choi HMT , Schwarzkopf M , Fornace ME , Acharya A , Artavanis G , Stegmaier J , Cunha A , Pierce NA . 2018. Third-generation in situ hybridization chain reaction: multiplexed, quantitative, sensitive, versatile, robust. Development (Cambridge, England) 145 :dev165753. doi:10.1242/dev.165753 29945988
Chung L. 2015. A Brief Introduction to the Transduction of Neural Activity into Fos Signal. Development & Reproduction 19 :61–67. doi:10.12717/DR.2015.19.2.061 27004262
Claudi F , Petrucco L , Tyson A , Branco T , Margrie T , Portugues R . 2020. BrainGlobe Atlas API: a common interface for neuroanatomical atlases. Journal of Open Source Software 5 :2668. doi:10.21105/joss.02668
Clayden J , Cox B , Jenkinson M , Reynolds R , Fissell K , Gailly J , Adler M . 2021. RNifti: Fast R and C++ access to NIfTI images. R package version 170 1 .
Corradi L , Bruzzone M , Maschio M dal , Sawamiphak S , Filosa A . 2022. Hypothalamic Galanin-producing neurons regulate stress in zebrafish through a peptidergic, self-inhibitory loop. Current Biology 32 :1497–1510.e5. doi:10.1016/j.cub.2022.02.011 35219430
Csardi G , Nepusz T . 2006. The igraph software package for complex network research. InterJournal, Complex Systems 1695 :1–9.
Ding JM , Carver WC , Terracio L , Buggy J . 1994. Proto-oncogene c-fos and the regulation of vasopressin gene expression during dehydration. Molecular Brain Research 21 :247–255. doi:10.1016/0169-328X(94)90255-0 8170349
Erdös P , Rényi A . 2011. On the evolution of random graphsOn the Evolution of Random Graphs. Princeton University Press. pp. 38–82. doi:10.1515/9781400841356.38
Fadok JP , Markovic M , Tovote P , Lüthi A . 2018. New perspectives on central amygdala function. Current Opinion in Neurobiology, Neurobiology of Behavior 49 :141–147. doi:10.1016/j.conb.2018.02.009 29522976
Fallani FDV , Latora V , Chavez M . 2017. A Topological Criterion for Filtering Information in Complex Brain Networks. PLOS Computational Biology 13 :e1005305. doi:10.1371/journal.pcbi.1005305 28076353
Gerlai R. 2023. Zebrafish (Danio rerio): A newcomer with great promise in behavioral neuroscience. Neuroscience & Biobehavioral Reviews 144 :104978. doi:10.1016/j.neubiorev.2022.104978 36442644
Gerlai R. 2020. Evolutionary conservation, translational relevance and cognitive function: The future of zebrafish in behavioral neuroscience. Neuroscience & Biobehavioral Reviews 116 :426–435. doi:10.1016/j.neubiorev.2020.07.009 32681940
Gerlai R. 2014. Social behavior of zebrafish: From synthetic images to biological mechanisms of shoaling. Journal of Neuroscience Methods, Measuring Behavior 234 :59–65. doi:10.1016/j.jneumeth.2014.04.028 24793400
Gholipour A , Kehtarnavaz N , Briggs R , Devous M , Gopinath K . 2007. Brain Functional Localization: A Survey of Image Registration Techniques. IEEE Transactions on Medical Imaging 26 :427–451. doi:10.1109/TMI.2007.892508 17427731
Grillner S , Robertson B . 2016. The Basal Ganglia Over 500 Million Years. Current Biology 26 :R1088–R1100. doi:10.1016/j.cub.2016.06.041 27780050
Guimerà R , Amaral LAN . 2005. Cartography of complex networks: Modules and universal roles. Journal of Statistical Mechanics: Theory and Experiment 2005 :1–13. doi:10.1088/1742-5468/2005/02/P02001
Guzowski JF , Setlow B , Wagner EK , McGaugh JL . 2001. Experience-Dependent Gene Expression in the Rat Hippocampus after Spatial Learning: A Comparison of the Immediate-Early GenesArc, c-fos, and zif268. J Neurosci 21 :5089–5098. doi:10.1523/JNEUROSCI.21-14-05089.2001 11438584
Herget U , Wolf A , Wullimann MF , Ryu S . 2014. Molecular neuroanatomy and chemoarchitecture of the neurosecretory preoptic-hypothalamic area in zebrafish larvae. Journal of Comparative Neurology 522 :1542–1564. doi:10.1002/cne.23480 24127437
Hillman EMC , Voleti V , Li W , Yu H . 2019. Light-Sheet Microscopy in Neuroscience. Annual Review of Neuroscience 42 :295–313. doi:10.1146/annurev-neuro-070918-050357
Humphries MD , Gurney K . 2008. Network ‘Small-World-Ness’: A Quantitative Method for Determining Canonical Network Equivalence. PLOS ONE 3 :e0002051. doi:10.1371/journal.pone.0002051 18446219
Jones LJ , Norton WHJ . 2015. Using zebrafish to uncover the genetic and neural basis of aggression, a frequent comorbid symptom of psychiatric disorders. Behavioural Brain Research, SI: Neuropsychiatric Spectra 276 :171–180. doi:10.1016/j.bbr.2014.05.055
Kacprzak V , Patel NA , Riley E , Yu L , Yeh J- RJ , Zhdanova IV . 2017. Dopaminergic control of anxiety in young and aged zebrafish. Pharmacology Biochemistry and Behavior 157 :1–8. doi:10.1016/j.pbb.2017.01.005 28408289
Kalueff AV , Gebhardt M , Stewart AM , Cachat JM , Brimmer M , Chawla JS , Craddock C , Kyzar EJ , Roth A , Landsman S , Gaikwad S , Robinson K , Baatrup E , Tierney K , Shamchuk A , Norton W , Miller N , Nicolson T , Braubach O , Gilman CP , Pittman J , Rosemberg DB , Gerlai R , Echevarria D , Lamb E , Neuhauss SCFF , Weng W , Bally-Cuif L , Schneider H , Schneider , and the Zebrafish Neuros H. 2013. Towards a comprehensive catalog of zebrafish behavior 1.0 and beyond. Zebrafish 10 :70–86. doi:10.1089/zeb.2012.0861 23590400
Kareklas K , Teles MC , Nunes AR , Oliveira RF . 2023. Social zebrafish: Danio rerio as an emerging model in social neuroendocrinology. Journal of Neuroendocrinology 35 :e13280. doi:10.1111/jne.13280 37165563
Kenney JW . 2020. Associative and nonassociative learning in adult zebrafishBehavioral and Neural Genetics of Zebrafish. Elsevier. pp. 187–204. doi:10.1016/b978-0-12-817528-6.00012-7
Kenney JW , Sorokina O , Genheden M , Sorokin A , Armstrong JD , Proud CG . 2015. Dynamics of Elongation Factor 2 Kinase Regulation in Cortical Neurons in Response to Synaptic Activity. J Neurosci 35 :3034–3047. doi:10.1523/JNEUROSCI.2866-14.2015 25698741
Kenney JW , Steadman PE , Young O , Shi MT , Polanco M , Dubaishi S , Covert K , Mueller T , Frankland PW . 2021. A 3D adult zebrafish brain atlas (AZBA) for the digital age. eLife 10 :e69988. doi:10.7554/eLife.69988 34806976
Klein A , Andersson J , Ardekani BA , Ashburner J , Avants B , Chiang M-C , Christensen GE , Collins DL , Gee J , Hellier P , Song JH , Jenkinson M , Lepage C , Rueckert D , Thompson P , Vercauteren T , Woods RP , Mann JJ , Parsey RV . 2009. Evaluation of 14 nonlinear deformation algorithms applied to human brain MRI registration. NeuroImage 46 :786–802. doi:10.1016/j.neuroimage.2008.12.037 19195496
Kovács KJ . 1998. Invited review c-Fos as a transcription factor: a stressful (re)view from a functional map. Neurochemistry International 33 :287–297. doi:10.1016/S0197-0186(98)00023-0 9840219
Kramer EE , Steadman PE , Epp JR , Frankland PW , Josselyn SA . 2018. Assessing Individual Neuronal Activity Across the Intact Brain: Using Hybridization Chain Reaction (HCR) to Detect Arc mRNA Localized to the Nucleus in Volumes of Cleared Brain Tissue. Current Protocols in Neuroscience 84 :e49. doi:10.1002/cpns.49 29944213
Kumar V , Krolewski DM , Hebda-Bauer EK , Parsegian A , Martin B , Foltz M , Akil H , Watson SJ . 2021. Optimization and evaluation of fluorescence in situ hybridization chain reaction in cleared fresh-frozen brain tissues. Brain Struct Funct 226 :481–499. doi:10.1007/s00429-020-02194-4 33386994
Loring MD , Thomson EE , Naumann EA . 2020. Whole-brain interactions underlying zebrafish behavior. Current Opinion in Neurobiology, Whole-brain interactions between neural circuits 65 :88–99. doi:10.1016/j.conb.2020.09.011 33221591
Luca RM , Gerlai R . 2012. In search of optimal fear inducing stimuli: Differential behavioral responses to computer animated images in zebrafish. Behavioural Brain Research 226 :66–76. doi:10.1016/j.bbr.2011.09.001 21920389
Lucibello FC , Lowag C , Neuberg M , Müller R . 1989. trans-repression of the mouse c-fos promoter: a novel mechanism of Fos-mediated trans-regulation. Cell 59 :999–1007.2513130
Marquart GD , Tabor KM , Horstick EJ , Brown M , Geoca AK , Polys NF , Nogare DD , Burgess HA . 2017. High-precision registration between zebrafish brain atlases using symmetric diffeomorphic normalization. GigaScience 6 . doi:10.1093/gigascience/gix056
Maruska KP , Butler JM , Field KE , Forester C , Augustus A . 2020. Neural Activation Patterns Associated with Maternal Mouthbrooding and Energetic State in an African Cichlid Fish. Neuroscience 446 :199–212. doi:10.1016/j.neuroscience.2020.07.025 32707292
Maximino C , Puty B , Benzecry R , Araújo J , Lima MG , de Jesus Oliveira Batista E , Renata de Matos Oliveira K , Crespo-Lopez ME , Herculano AM . 2013. Role of serotonin in zebrafish (Danio rerio) anxiety: Relationship with serotonin levels and effect of buspirone, WAY 100635, SB 224289, fluoxetine and para-chlorophenylalanine (pCPA) in two behavioral models. Neuropharmacology 71 :83–97. doi:10.1016/j.neuropharm.2013.03.006 23541719
Mueller T. 2022. The Everted Amygdala of Ray-Finned Fish: Zebrafish Makes a Case. Brain Behavior and Evolution 97 :321–335. doi:10.1159/000525669 35760049
Murphy K , van Ginneken B , Reinhardt JM , Kabus S , Ding K , Deng X , Cao K , Du K , Christensen GE , Garcia V , Vercauteren T , Ayache N , Commowick O , Malandain G , Glocker B , Paragios N , Navab N , Gorbunova V , Sporring J , de Bruijne M , Han X , Heinrich MP , Schnabel JA , Jenkinson M , Lorenz C , Modat M , McClelland JR , Ourselin S , Muenzing SEA , Viergever MA , De Nigris D , Collins DL , Arbel T , Peroni M , Li R , Sharp GC , Schmidt-Richberg A , Ehrhardt J , Werner R , Smeets D , Loeckx D , Song G , Tustison N , Avants B , Gee JC , Staring M , Klein S , Stoel BC , Urschler M , Werlberger M , Vandemeulebroucke J , Rit S , Sarrut D , Pluim JPW . 2011. Evaluation of Registration Methods on Thoracic CT: The EMPIRE10 Challenge. IEEE Transactions on Medical Imaging 30 :1901–1920. doi:10.1109/TMI.2011.2158349 21632295
Nabinger DD , Altenhofen S , Buatois A , Facciol A , Peixoto JV , da Silva JMK , Chatterjee D , Rübensam G , Gerlai R , Bonan CD . 2023. Acute administration of a dopamine D2/D3 receptor agonist alters behavioral and neural parameters in adult zebrafish. Progress in Neuro-Psychopharmacology and Biological Psychiatry 125 :110753. doi:10.1016/j.pnpbp.2023.110753 36934998
Nowicki M , Tran S , Muraleetharan A , Markovic S , Gerlai R . 2014. Serotonin antagonists induce anxiolytic and anxiogenic-like behavior in zebrafish in a receptor-subtype dependent manner. Pharmacology Biochemistry and Behavior 126 :170–180. doi:10.1016/j.pbb.2014.09.022 25284132
Palieri V , Paoli E , Wu YK , Haesemeyer M , Kadow ICG , Portugues R . 2024. The preoptic area and dorsal habenula jointly support homeostatic navigation in larval zebrafish. Current Biology 34 :489–504.e7. doi:10.1016/j.cub.2023.12.030 38211586
Parada CDC , Mayer U , Chagnaud BP . 2024. The Dorsal Part of the Anterior Tuberal Nucleus Responds to Auditory Stimulation in Zebrafish (Danio rerio). eNeuro 11 . doi:10.1523/ENEURO.0062-24.2024
Pinho JS , Cunliffe V , Kareklas K , Petri G , Oliveira RF . 2023. Social and asocial learning in zebrafish are encoded by a shared brain network that is differentially modulated by local activation. Commun Biol 6 :1–13. doi:10.1038/s42003-023-04999-5 36596887
Porter BA , Mueller T . 2020. The Zebrafish Amygdaloid Complex – Functional Ground Plan, Molecular Delineation, and Everted Topology. Frontiers in Neuroscience 14 :608. doi:10.3389/fnins.2020.00608 32765204
Portugues R , Feierstein CE , Engert F , Orger MB . 2014. Whole-Brain Activity Maps Reveal Stereotyped, Distributed Networks for Visuomotor Behavior. Neuron 81 :1328–1343. doi:10.1016/j.neuron.2014.01.019 24656252
Raam T , Hong W . 2021. Organization of neural circuits underlying social behavior: A consideration of the medial amygdala. Current Opinion in Neurobiology, The Social Brain 68 :124–136. doi:10.1016/j.conb.2021.02.008 33940499
Rajput N , Parikh K , Kenney JW . 2022. Beyond bold versus shy: Zebrafish exploratory behavior falls into several behavioral clusters and is influenced by strain and sex. Biology Open 11 :bio059443. doi:10.1242/bio.059443 36039864
Randlett O , Wee CL , Naumann EA , Nnaemeka O , Schoppik D , Fitzgerald JE , Portugues R , Lacoste AMB , Riegler C , Engert F , Schier AF . 2015. Whole-brain activity mapping onto a zebrafish brain atlas. Nature Methods 12 :1039–1046. doi:10.1038/nmeth.3581 26778924
Renier N , Adams EL , Kirst C , Wu Z , Azevedo R , Kohl J , Autry AE , Kadiri L , Umadevi Venkataraju K , Zhou Y , Wang VX , Tang CY , Olsen O , Dulac C , Osten P , Tessier-Lavigne M . 2016. Mapping of Brain Activity by Automated Volume Analysis of Immediate Early Genes. Cell 165 :1789–1802. doi:10.1016/j.cell.2016.05.007 27238021
Renier N , Wu Z , Simon DJ , Yang J , Ariel P , Tessier-Lavigne M . 2014. IDISCO: A simple, rapid method to immunolabel large tissue samples for volume imaging. Cell 159 :896–910. doi:10.1016/j.cell.2014.10.010 25417164
Reynaud EG , Peychl J , Huisken J , Tomancak P . 2014. Guide to light-sheet microscopy for adventurous biologists, Nature Publishing Group. doi:10.1038/nmeth.3222
Richardson DS , Guan W , Matsumoto K , Pan C , Chung K , Ertürk A , Ueda HR , Lichtman JW . 2021. Tissue clearing. Nat Rev Methods Primers 1 :1–24. doi:10.1038/s43586-021-00080-9
Scaia MF , Akinrinade I , Petri G , Oliveira RF . 2022. Sex Differences in Aggression Are Paralleled by Differential Activation of the Brain Social Decision-Making Network in Zebrafish. Front Behav Neurosci 16 . doi:10.3389/fnbeh.2022.784835
Shainer I , Kuehn E , Laurell E , Al Kassar M , Mokayes N , Sherman S , Larsch J , Kunst M , Baier H . 2023. A single-cell resolution gene expression atlas of the larval zebrafish brain. Science Advances 9 :eade9909. doi:10.1126/sciadv.ade9909 36812331
Spence R , Fatema MK , Reichard M , Huq KA , Wahab MA , Ahmed ZF , Smith C . 2006. The distribution and habitat preferences of the zebrafish in Bangladesh. Journal of Fish Biology 69 :1435–1448. doi:10.1111/j.1095-8649.2006.01206.x
Toms CN , Echevarria DJ . 2014. Back to basics: Searching for a comprehensive framework for exploring individual differences in zebrafish (Danio Rerio) behavior. Zebrafish 11 :325–340. doi:10.1089/zeb.2013.0952 24921670
Tyson AL , Rousseau CV , Niedworok CJ , Keshavarzi S , Tsitoura C , Cossell L , Strom M , Margrie TW . 2021. A deep learning algorithm for 3D cell detection in whole mouse brain image datasets. PLOS Computational Biology 17 :e1009074. doi:10.1371/journal.pcbi.1009074 34048426
Venincasa MJ , Randlett O , Sumathipala SH , Bindernagel R , Stark MJ , Yan Q , Sloan SA , Buglo E , Meng QC , Engert F , Züchner S , Kelz MB , Syed S , Dallman JE . 2021. Elevated preoptic brain activity in zebrafish glial glycine transporter mutants is linked to lethargy-like behaviors and delayed emergence from anesthesia. Sci Rep 11 :3148. doi:10.1038/s41598-021-82342-w 33542258
Vetere G , Kenney JW , Tran LM , Xia F , Steadman PE , Parkinson J , Josselyn SA , Frankland PW . 2017. Chemogenetic Interrogation of a Brain-wide Fear Memory Network in Mice. Neuron 94 :363–374.e4. doi:10.1016/j.neuron.2017.03.037 28426969
Wheeler AL , Teixeira CM , Wang AH , Xiong X , Kovacevic N , Lerch JP , McIntosh AR , Parkinson J , Frankland PW . 2013. Identification of a functional connectome for long-term fear memory in mice. PLoS Comput Biol 9 :e1002853.23300432
Wong K , Elegante M , Bartels B , Elkhayat S , Tien D , Roy S , Goodspeed J , Suciu C , Tan J , Grimes C , Chung A , Rosenberg M , Gaikwad S , Denmark A , Jackson A , Kadri F , Chung KM , Stewart A , Gilder T , Beeson E , Zapolsky I , Wu N , Cachat J , Kalueff AV . 2010. Analyzing habituation responses to novelty in zebrafish (Danio rerio). Behavioural Brain Research 208 :450–457. doi:10.1016/j.bbr.2009.12.023 20035794
Yushkevich PA , Pashchinskiy A , Oguz I , Mohan S , Schmitt JE , Stein JM , Zukić D , Vicory J , McCormick M , Yushkevich N , Schwartz N , Gao Y , Gerig G . 2019. User-Guided Segmentation of Multi-modality Medical Imaging Datasets with ITK-SNAP. Neuroinformatics 17 :83–102. doi:10.1007/s12021-018-9385-x 29946897
Zangenehpour S , Chaudhuri A . 2002. Differential induction and decay curves of c-fos and zif268 revealed through dual activity maps. Molecular Brain Research 109 :221–225. doi:10.1016/S0169-328X(02)00556-9 12531532
