The organization of serotonergic fibers in the Pacific angelshark brain: Neuroanatomical and supercomputing analyses
Data files
Jul 14, 2025 version files 1.96 GB
-
PacificAngelshark_CoronalSectionSet_Nissl.zip
734.44 MB
-
PacificAngelshark_CoronalSectionSet_UnstainedWet.zip
634.23 MB
-
PacificAngelshark_SimulatedFiberDensities.zip
592.04 MB
-
Reader_SimulatedFiberDensities.nb
13.29 KB
-
README.md
4.83 KB
Abstract
Serotonergic axons (fibers) are a universal feature of all vertebrate brains. They form meshworks, typically quantified with regional density measurements, and appear to support neuroplasticity. The self-organization of this system remains poorly understood, partly because of the strong stochasticity of individual fiber trajectories. In an extension to our previous analyses of the mouse brain, serotonergic fibers were investigated in the brain of the Pacific angelshark (Squatina californica), a representative of a unique (ray-like) lineage of the squalomorph sharks. The fundamental cytoarchitecture of the angelshark brain was examined, including the expression of ionized calcium binding adaptor molecule 1 (Iba1, AIF-1), and serotonergic fibers were visualized with immunohistochemistry. In the forebrain, serotonergic fibers had the tendency to advance toward the dorsal pallium and also accumulated at significantly higher densities at pial borders. Next, a population of serotonergic fibers was modeled inside a digital model of the angelshark brain by using a supercomputing simulation. The simulated fibers were defined as sample paths of fractional Brownian motion (FBM), a continuous-time stochastic process. The results reproduced key features of serotonergic fiber densities in the telencephalon, a brain division with a considerable physical uniformity and no major “obstacles” (dense axon tracts). The study provides further evidence that serotonergic fibers can be successfully modeled as paths of a rigorously defined stochastic process, and that a rich repertoire of self-organizing behaviors can be produced by axons that are inherently stochastic but also respond to external forces.
Authors: Skirmantas Janusonis, Ralf Metzler, Thomas Vojta
FILE: PacificAngelshark_CoronalSectionSet_Nissl.zip
This set of images is a series of coronal sections through a Pacific angelshark (Squatina californica) adult female brain. The specimen was collected in Santa Barbara County (California, USA) in September 2017.
The brain was embedded in gelatin (in two blocks) and sectioned at 40-micrometer thickness on a freezing microtome. The sections have been stained with a Nissl dye (thionine).
The sections are in the sequential rostro-caudal order, separated by 160 micrometers. The first section (“PacificAngelshark_Sp24_Thionine_1s1”) contains a scale bar that applies to all sections with no magnification in the name. These images have been imaged in bright field with a 1X objective. Some images have two partially overlapping parts (“p1” and “p2”).
Some details are shown at higher magnifications. These images have been captured in bright field and contain a scale bar (their names also include the magnification; e.g., “_5X”).
Each physical slide had 6 mounted sections. The file name “PacificAngelshark_Sp24_Thionine_AsB” refers to the shark specimen number (“Sp24”), the slide number (A), and the section number (B) on this slide. Some sections in the middle of the set are incomplete (e.g., “14s1”) because they were at the beginning of the second block.
The imaging system: Zeiss AxioImager Z1, with the following objectives: 1X (Plan-Neofluar, NA = 0.025), 5X (EC Plan-Neofluar, NA = 0.16), 10X (Plan-Apochromat, NA = 0.45), and 20X (Plan-Apochromat, NA = 0.80).
All sections have been processed and imaged by Dr. Skirmantas Janusonis (University of California, Santa Barbara).
FILE: PacificAngelshark_CoronalSectionSet_UnstainedWet.zip
This set of images is a series of coronal sections through a Pacific angelshark (Squatina californica) adult female brain. The specimen was collected in Santa Barbara County (California, USA) in September 2017.
The brain was embedded in gelatin (in two blocks) and sectioned at 40-micrometer thickness on a freezing microtome. The sections have been imaged on glass slides, unstained and uncoverslipped, after a quick dip in water.
The sections are in the sequential rostro-caudal order, separated by 160 micrometers. The first section (“PacificAngelshark_Sp24_Thionine_1s1”) contains a scale bar that applies to all sections with no magnification in the name. These images have been imaged in bright field with a 1X objective. Some images have two partially overlapping parts (“p1” and “p2”).
Each physical slide had 6 mounted sections. The file name “PacificAngelshark_Sp24_Thionine_AsB” refers to the shark specimen number (“Sp24”), the slide number (A), and the section number (B) on this slide. Some sections in the middle of the set are incomplete (e.g., “14s1”) because they were at the beginning of the second block.
The imaging system: Zeiss AxioImager Z1, with a 1X objective (Plan-Neofluar, NA = 0.025).
All sections have been processed and imaged by Dr. Skirmantas Janusonis (University of California, Santa Barbara).
FILE: PacificAngelshark_SimulatedFiberDensities.zip
This dataset contains the simulated densities of 4800 serotonergic fibers modeled as fractional Brownian motion paths with Hurst index H = 0.8. Each path consists of 2^25 (approximately 33.6 million) steps. In the simulations, distances were measured in “pixel units” whose size corresponds to the size of the pixels in the images used to create the 3D-brain model. One pixel unit corresponds to 12.5 micrometers in the physical angelshark brain.
The fiber densities were recorded using a cubic grid with cells of a linear size of 25.0 micrometers (i.e., 2 pixel-units). The densities are stored in the files 011.dat through 999.dat. Each file contains the densities at one rostrocaudal level (one coronal section). The filename indicates the rostrocaudal position in pixel-units.
Each data file is organized in columns. The first three columns (igx, igy, igz) contain the labels (counters) of the grid cells in the three dimensions. The next three columns (x, y, x) contain the position of the center of each cell in pixel-units. The seventh column (distrib) contains the fiber density, and the eighth column (log(distrib)) contains the natural log of distrib (or -100 if the density is zero). The densities are normalized as probability densities (i.e., the integral of the density over the entire shape is normalized to unity).
FILE: Reader_SimulatedFiberDensities.nb
This script (written in Wolfram Mathematica 14) reads simulated fiber densities (in the *.dat format) and outputs them as TIFF images.