Atlas pipeline workflow

Last modified by drodarie on 2023/11/13 12:31

Creating a mouse cell atlas model

From annotation displaying to neuron density estimations

workflow_complete_plot.jpg

Content

This collab is meant for modelists interested in using Blue Brain atlas suite [1] and the Brain Scaffold Builder (BSB) [2] to build circuits of the mouse brain.

This collab contains notebooks to present the mouse brain cell atlas pipeline: its steps and results.

Through these notebooks, you should learn:

  • How to manipulate the Allen Brain reference atlases and brain region hierarchy, using the voxcell library.
  • How to realign ISH datasets to the Nissl reference atlas (introduction to Deep-Atlas [3]).
  • How to compute the orientation field followed by the neurons inside a region (introduction to atlas-direction-vectors).
  • How to compute the depth field following the orientation field.
  • How to compute cell, neuron and glia densities in each region of the mouse brain (first introduction to atlas-densities).
  • How to estimate neuron sub-type densities based on literature findings (second introduction to atlas-densities).
  • How to refine neuron densities based on linear densities.
  • How to scale, rotate and bend neuron morphologies, based on the orientation and depth fields.
  • How to construct a circuit based on the resulting cell atlas model in BSB.

Each item in the list correspond to a notebook in the lab. Items in grey are notebooks still in preparation and will be integrated in the HBP collab, as soon as possible.

Jupyter notebooks

Inside the lab, you should find notebooks numbered so that they follow the natural steps of the pipeline.

Some parts of the notebooks might require more than the 2GB RAM memory available on the HBP lab. These steps will be highlighted in the notebook.

About the figure.

Figure 1. Mouse declive cell placement workflow: 

Each image represents a sagittal slice of the result of the step of our process to reconstruct the mouse declive. Black arrows show the order of the workflow steps. 

A. Our corrected annotation atlas of declive is shown in colors over the Nissl volume in levels of grey. The granular layer appears in orange, the molecular layer in blue and Purkinje layer in green. 

B. To each voxel of the declive, a 3D direction normalized vector is computed corresponding to the main axis of the axons in the region. Colors represent the orientation vectors norm on their respective plane, black arrows their projected vector. 

C. Distance of each voxel of declive to the outside border of the molecular layer, following the orientation field, expressed in micrometers. 

D. E. Respectively Neuron and inhibitory neuron density in logarithmic scale. 

F. Soma position of the different neuron types of the declive, displayed over the annotation atlas. Each cell type appears in a different color and size corresponding to its radius. The annotation atlas’ colors correspond to A. 

G. Projection of the Purkinje cell morphology displayed in colors over the annotation atlas. Each morphology has been rotated, scaled, and bended following the orientation and depth fields. The annotation atlas’ colors correspond to A. 

References.

  1. D. Rodarie et al., “A method to estimate the cellular composition of the mouse brain from heterogeneous datasets,” PLoS Comput Biol, vol. 18, no. 12, p. e1010739, Dec. 2022, doi: 10.1371/journal.pcbi.1010739.
  2. R. De Schepper et al., “Model simulations unveil the structure-function-dynamics relationship of the cerebellar cortical microcircuit,” Commun Biol, vol. 5, no. 1, Art. no. 1, Nov. 2022, doi: 10.1038/s42003-022-04213-y.

  3. J. Krepl et al., “Supervised Learning With Perceptual Similarity for Multimodal Gene Expression Registration of a Mouse Brain Atlas,” Front. Neuroinform., vol. 15, p. 691918, Jul. 2021, doi: 10.3389/fninf.2021.691918.