Skip to main content


Front. Cell Dev. Biol., 08 June 2022
Sec. Morphogenesis and Patterning
This article is part of the Research Topic Correlative Light and Volume Electron Microscopy: Methods and Applications View all 10 articles

Sample Preparation and Warping Accuracy for Correlative Multimodal Imaging in the Mouse Olfactory Bulb Using 2-Photon, Synchrotron X-Ray and Volume Electron Microscopy

  • 1Sensory Circuits and Neurotechnology Lab, The Francis Crick Institute, London, United Kingdom
  • 2Department of Neuroscience, Physiology and Pharmacology, University College London, London, United Kingdom
  • 3ESRF, The European Synchrotron, Grenoble, France
  • 4Department of Physics and Astronomy, University College London, London, United Kingdom
  • 5Diamond Light Source, Harwell Science and Innovation Campus, Didcot, United Kingdom
  • 6School of Physics and Astronomy, University of Southampton, Highfield Campus, Southampton, United Kingdom
  • 7Paul Scherrer Institut, Villigen, Switzerland

Integrating physiology with structural insights of the same neuronal circuit provides a unique approach to understanding how the mammalian brain computes information. However, combining the techniques that provide both streams of data represents an experimental challenge. When studying glomerular column circuits in the mouse olfactory bulb, this approach involves e.g., recording the neuronal activity with in vivo 2-photon (2P) calcium imaging, retrieving the circuit structure with synchrotron X-ray computed tomography with propagation-based phase contrast (SXRT) and/or serial block-face scanning electron microscopy (SBEM) and correlating these datasets. Sample preparation and dataset correlation are two key bottlenecks in this correlative workflow. Here, we first quantify the occurrence of different artefacts when staining tissue slices with heavy metals to generate X-ray or electron contrast. We report improvements in the staining procedure, ultimately achieving perfect staining in ∼67% of the 0.6 mm thick olfactory bulb slices that were previously imaged in vivo with 2P. Secondly, we characterise the accuracy of the spatial correlation between functional and structural datasets. We demonstrate that direct, single-cell precise correlation between in vivo 2P and SXRT tissue volumes is possible and as reliable as correlating between 2P and SBEM. Altogether, these results pave the way for experiments that require retrieving physiology, circuit structure and synaptic signatures in targeted regions. These correlative function-structure studies will bring a more complete understanding of mammalian olfactory processing across spatial scales and time.


Maps are essential tools to accelerate discovery. By providing information in an intuitive structure, they help design strategies that maximise gain while minimising risk. As it has occurred in a myriad of disciplines before, maps play a crucial role in supporting discovery in the life sciences and in neuroscience in particular (Niedworok et al., 2016; Winnubst et al., 2019; Claudi et al., 2020; Ueda et al., 2020; BRAIN Initiative Cell Census Network (BICCN), 2021). Mapping the brain at multiple scales can help gain mechanistic insights into its function from already existing data, and guide the exploration of the uncharted territories that hold the keys for further unsolved functions. Creating those maps requires recording the spatial distribution of relevant brain properties across a volume that encodes functional units at a resolution sufficient to discriminate between individual features. In the mammalian brain, neuronal circuits usually span >1 mm3, while the dimensions of relevant features (nuclei, dendrites, axons, synapses) range between 10 μm and 10 nm (Helmstaedter, 2013; Klinger et al., 2021). Targeted maps carrying more detail (only) where it is needed offer a valuable solution to this vast challenge, and can be generated by means of correlative multimodal imaging (CMI) pipelines (Bock et al., 2011; Briggman et al., 2011; Karreman et al., 2016; Lee et al., 2016; Dyer et al., 2017; Bosch et al., 2022). The reliability of those brain maps will be directly influenced by two main factors: the quality of the images of each individual modality and the quality of their registration into a multimodal map.

The quality of an image is limited by multiple factors including optics, detector sensitivity, exposure time or the structure of the specimen leading to signal contrast. The latter is engineered and optimised by preparing the specimen for a particular imaging technique in a process often referred to as the sample preparation protocol (SPP). An SPP incorporates all transformations necessary to make a specimen compatible with an imaging technique and may involve mechanical procedures (such as slicing), chemical transformations (such as fixing) and signal transformations (such as staining). SPPs can be studied as a sequence of smaller-scale operations sharing a single goal, which we will refer to as “steps”. In that way, a single step may involve slicing a fresh piece of brain tissue to a specific thickness under controlled constraints of temperature, osmolarity and time, while another might involve chemically staining a tissue slice with a buffered solution of osmium tetroxide while controlling for temperature, duration and agitation. Optimising an SPP presents therefore a multivariate challenge that can be approached by making informed decisions to optimise each step individually while monitoring global yields.

The maps guiding our understanding of how neuronal circuits operate in the mammalian brain should convey spatial and functional information across the critical dimensions specified above, recording tissue volumes in the order of 1 mm3. While not a single technique is capable of providing structural and physiological insights at the required scales, two techniques are capable of providing both streams of data separately: serial block-face scanning electron microscopy (SBEM) is the volume electron microscopy technique capable of retrieving the largest continuous tissue volume with minimal specimen distortion (Denk and Horstmann, 2004; Kubota, 2015), whereas multiphoton microscopy offers a window to explore the physiological events taking place up to 1 mm deep inside tissue (Denk et al., 1994; Oheim et al., 2001; Yildirim et al., 2019). Finally, X-ray microscopy is an outstanding resource to reveal subcellular anatomy across large specimens (Kuan et al., 2020; Bosch et al., 2022). Altogether, the previously defined challenge can be addressed by a CMI pipeline that combines in vivo 2-photon (2P) microscopy of intracellular [Ca2+], synchrotron X-ray computed tomography with propagation-based phase contrast (SXRT) and SBEM (Bosch et al., 2022). The SPP for this approach will therefore need to be compatible with all techniques involved and its optimal implementation will be the one that accomplishes two main requisites: 1) being able to deliver an optimal final map, therefore balancing inter-modality trade-offs towards the benefit of the initial research question; and 2) being able to provide the highest yield, making it compatible with more challenging research questions in the future.

The datasets carrying the best spatial resolution in this CMI pipeline will be provided by SBEM (Bosch et al., 2022). In these, the signal is generated by recording the back-scattered electrons arising from a 1–2 keV electron beam scanning the surface of a polished tissue block at an electron dose of typically <20 e/nm2 (Kubota, 2015). For soft tissue to generate contrast under this regime, tissue is stained with heavy metals (Hayat, 2000). Heavy metals bind non-specifically to lipids and proteins, ultimately providing a sample with a topology that faithfully reproduces its ultrastructure (Buckley and Porter, 1967). SBEM adds one challenge to otherwise more conventional transmission electron microscopy imaging: here, part of the electron beam will be absorbed by the sample, locally increasing its charge and potentially influencing the trajectories of neighbouring back-scattered electrons. In order to prevent such charging artefacts, SBEM-specific SPP optimisations have incorporated elements that increase the conductivity of the sample — such as increasing the total metal load employed in staining, in terms of both the number of metallic compounds and their concentrations (Tapia et al., 2012). The diffusion mechanism used to deliver those metals through the tissue remains passive, and therefore incubation time is subject to a quadratic relationship with the distance each compound should travel (Ströh et al., 2021), imposing severe temporal burdens for staining specimens of dimensions beyond 1 mm (Mikula and Denk, 2015). Moreover, abnormalities arising from individual staining steps can interact with each other (Hua et al., 2015; Mikula and Denk, 2015), creating an additional challenge to ensure all compounds optimally stain the tissue across all dimensions. Nevertheless, some SPPs have reported acceptable SBEM results in mouse brain tissue samples of thicknesses ranging from 0.1 mm to a complete mouse brain (∼10 mm) (Briggman et al., 2011; Tapia et al., 2012; Wilke et al., 2013; Hua et al., 2015; Mikula and Denk, 2015; Pallotto et al., 2015; Wanner et al., 2016; Genoud et al., 2018), providing a promising starting point for further optimisations aiming to obtain a multimodal brain map. Interestingly, such SBEM staining also provides good contrast in SXRT (Bushong et al., 2015) because X-ray absorption is proportional to the atomic number. Hence heavy metal-stained samples of thicknesses <1 mm can be imaged efficiently with this technique.

Finally, the registration quality of the individual map layers needs to be accurate enough to correlate any feature of interest across all modalities. For maps aimed at pinpointing the structure of a neural circuit, registration should be accurate at single-cell precision. Since the largest feature to unequivocally identify a single cell within a circuit is its nucleus, the requirement becomes registering at a precision better than the inter-nucleus distance, which is approximately 20 µm in the mitral cell layer of the olfactory bulb (OB) (Panhuber et al., 1985; Frazier and Brunjes, 1988; Parrish-Aungst et al., 2007).

Here we show an optimization of a CMI pipeline aiming to map a mammalian brain circuit with 2P microscopy, SXRT and SBEM employing an automated tissue processor capable of staining up to 30 identified samples simultaneously. Using the mouse OB, we assessed and made improvements on the SPPs and quantified the registration accuracy across imaging modalities. We demonstrate that this pipeline is effective and that single-cell registration accuracy can be readily achieved between 2P and SXRT.

Materials and Methods

All SPPs used in this study are described in detail in Supplementary Table S2 and are available for further analysis (see Data availability statement). The experimental details of samples shown in the figures are listed in Supplementary Table S3. The acronyms used are listed in Supplementary Table S8.


All animals used were 8–24 week-old mice of C57/Bl background and mixed gender. Details of the transgenic mouse lines are listed in Supplementary Table S4. All animal protocols were approved by the Ethics Committee of the board of the Francis Crick Institute and the United Kingdom Home Office under the Animals (Scientific Procedures) Act 1986.

Craniotomy and Durectomy Surgery

In experiments involving in vivo 2-photon (2P) imaging, anaesthesia was induced and maintained with intraperitoneal injection of sleep mix (0.05 mg/kg fentanyl, 5 mg/kg midazolam, 0.5 mg/kg medetomidine) and its depth monitored by toe-pinch reflex tests. Body temperature of mice was measured with a rectal temperature probe and maintained at 36–38°C with a thermoregulator (DC Temperature Controller, FHC, ME United States) heat pad.

Following shaving, skin disinfection with 1% chlorhexidine and skin incision over the skull, a custom head-fix implant with a central hole was fixed to the skull with medical super glue (Vetbond, 3 M, Maplewood MN, United States) and dental cement (Paladur, Heraeus Kulzer GmbH, Hanau, Germany; Simplex Rapid Liquid, 437 Associated Dental Products Ltd., Swindon, United Kingdom) so that the nasal bone above olfactory bulbs were exposed through the headplate hole. A ∼2 mm-diameter craniotomy over the left olfactory bulb was made with a dental drill (Success 40, Osada, Tokyo, Japan). ACSF (10 mM HEPES, 10 mM glucose, 125 mM NaCl, 5 mM KCl, 2 mM MgSO4.7H2O, 2 mM CaCl2.2H2O, pH = 7.4) was used to moisten the brain throughout the surgery. Exposed dura was removed with forceps and needles. Then a drop of warm 2% low-melt agarose diluted in ACSF was applied, quickly followed by the placement of a glass window cut from a coverslip [borosilicate glass #1 thickness (150 μm)]. The edge of the window was glued to the nasal bones with medical super glue (Vetbond, 3 M, Maplewood MN, United States).

Some SPPs also include measures to prevent the formation of edema, including the use of ice-cold ACSF during surgery to cool down the skull from drilling and the use of dexamethasone (20 μl/mouse intraperitoneal injection) 2–4 h prior to surgery or dissection (Supplementary Table S2).

In vivo 2-Photon Imaging

Following surgery, animals were transferred to a 2P microscope (Scientifica Multiphoton VivoScope coupled with a MaiTai DeepSee Laser tuned to 910–940 nm, <30 mW underneath objective) and held by a custom head-fix implant holder. Respiration was measured at the right nostril with a flow sensor (A3100, Honeywell, NC, United States) and the signal was digitised with a Power 1401 ADC board (CED, Cambridge, United Kingdom). Odours were delivered at inhalation onset through custom-made olfactometers (Ackels et al., 2021; Bosch et al., 2022) whose output port was <1 cm away from the left nostril. Images (512 × 512 pixels, field of view 550 × 550 μm) were taken with a resonant-galvo scanner at 30 Hz frame rate. Six or twelve planes covering ∼300 μm in z, including the mitral cell layer, were acquired with a 16 × 0.8 NA water-immersion objective (Nikon) coupled to a piezo motor (PI Instruments, United Kingdom). The fluorescently labelled glomerulus was used to position the field-of-view.

Following functional imaging, blood vessels were fluorescently labelled by injecting the tracer sulforhodamine 101 intraperitoneally (50 µl of a 0.1 M solution) (Sigma Aldrich). A dataset with voxels sized 0.358 µm in x, y and 5 µm in z containing the entire functionally-recorded volume was then acquired. Red and green channels were used to capture the vasculature and the resting GCaMP signals, respectively.

Dissection, Slicing and Fixation

For immersion fixation, animals were culled by cervical dislocation or decapitation after intraperitoneal injection of an overdose of pentobarbital (600 mg/kg pentobarbitone with mepivacaine) or sleep mix (0.05 mg/kg fentanyl, 5 mg/kg midazolam, 0.5 mg/kg medetomidine). Ice-cold dissection solution (4.6% sucrose, 0.563 mM CaCl2, 65.2 mM NaH2PO4.H2O, with 0.02% sodium azide, bubbled with 95% O2/5% CO2, osmolarity 280–320 mOsm/l) was used both to keep the tissue cold during dissection and as the bathing solution during vibratome slicing (Leica VT1200S). Dorsal OB slices with nominal thicknesses of 0.4–1 mm were taken by near-horizontal slicing, setting the vibratome blade tangent to the OB surface. Half and whole OB samples were cut with a scalpel. Specimens were then transferred to cold fixative (see below) and post-fixed overnight at 4°C. The time between sacrifice and fixation was kept to a minimum (achieving 10 min on average for the first samples dissected and sliced from an animal).

For perfusion, animals were anaesthetized with an overdose of pentobarbital (600 mg/kg pentobarbitone with mepivacaine) administered intraperitoneally. The skin over the chest was removed and incisions made medially-laterally across the abdomen muscle and across the diaphragm. Before cardiac arrest took place, room temperature fixative (see below) was delivered by a cannula inserted into the left atrium of the heart driven by a peristaltic pump (4.5–7.2 ml/min, 80–100 ml). The brain was dissected out and post-fixed overnight at 4°C and OBs vibratome-sliced on the following day in 0.15 M NCB [sodium cacodylate (aq.), pH = 7.4, 280–320 mOsm/l].

The fixatives used were:

a) 1% or 2% glutaraldehyde (GA) with 2 mM CaCl2 and 0.02% sodium azide in 0.15 M sodium cacodylate buffer (NCB), based on Pallotto et al. (2015).

b) 1.25% GA and 2.5% paraformaldehyde (PFA) or 2% GA and 2.5% PFA with 2 mM CaCl2 in 0.08 M NCB, based on Hua et al. (2015).

c) 0.5% GA or 1.25% GA or 2% GA or 1% PFA or 2.5% PFA or 4% PFA with 2 mM CaCl2 in 0.08 M NCB.

Epifluorescence and ex vivo 2P Imaging

For samples containing a fluorescently-labelled glomerulus, after a brief contact with fixative, epifluorescence images tiling the entire specimen were recorded with a widefield fluorescence microscope (Olympus Axioplan2 with LEJ “ebq 100 isolated-z” lamp) before overnight postfixation.

After postfixation, two volumes covering the fluorescently-labelled glomerulus (located using the epifluorescence images as guides) and its surrounding regions were acquired and stitched together (voxels sized 1.074 µm in x, y and 5 µm in z) with the aforementioned 2P system (laser wavelength 800 nm). This dataset featured glomeruli and surface blood vessels from the autofluorescence signal and was used to aid the relocation of functionally imaged volume in later acquired SXRT/EM datasets as described in Bosch et al. (2022).

Staining, Dehydration and Embedding

Fixed OB slices were washed in 0.15 M NCB [sodium cacodylate (aq.) with or without 0.02% sodium azide, pH = 7.4, 280–320 mOsm/l, 3 × 1 h], then stained using an automated tissue processor (Leica EMTP). The full descriptions can be found in Supplementary Table S2. In brief, specimens were processed with one of the following staining protocols:

a) 2% OsO4, 3% potassium ferrocyanide, 2 mM CaCl2 in 0.15 M NCB (1 h, 4°C); 1% thiocarbohydrazide (TCH, aq., 1 h, 50°C); 2% OsO4 (aq., 1 h, 20°C); 2% uranyl acetate (UA, aq., overnight, 4°C); lead aspartate (aq., pH = 5.5, 2 h, 60°C) with water washes between each step, protocol referred to as “1sRO half time”.

b) same as a) but doubled incubation times of all steps except UA, protocol referred to as “1sRO”.

c) 2% OsO4 in 0.15 M NCB (1.5 h, 20°C) followed by 2.5% potassium ferrocyanide in 0.15 M NCB (1.5 h, 20°C) without water wash; 1% TCH (aq., 45 min, 40°C); 2% OsO4 (aq., 1.5 h, 20°C); 1% UA (aq.) warmed to 50°C for 2 h after overnight at 4°C; lead aspartate (aq., pH = 5.0, 2 h, 50°C) with water washes between each step, protocol referred to as “2sRO”.

d) same as c) but TCH (aq.) incubated for 1 h at 50°C and lead aspartate (pH = 5.5) at 60°C, used in exp. A and exp. B in Figure 3B1.

e) same as c) but all warm incubations happening at 60°C and all steps performed manually. There were 6 groups of samples in this SPP, corresponding to terminating the staining at each of the 6 reagents (Figure 3D1). Early finishing samples were pooled in water at room temperature while others finished their course of staining.

Samples were then dehydrated with 75%, 90%, 100%, 100% ethanol series, transferred to propylene oxide and infiltrated with 25%, 50%, 75%, 100%, 100% hard epon diluted in propylene oxide, and finally polymerized for 72 h at 60–70°C. For experiments in Figures 3D1,2, samples were dehydrated with 100% ethanol at 4°C overnight before being transferred to propylene oxide.

The use of an automated tissue processor enabled processing a large amount of identified samples simultaneously (theoretically up to 48, in this study the maximum processed in parallel was 30). This implementation allowed us to systematically compare different sources of staining variability at scale with samples processed in parallel. In turn, the optimised protocols obtained with this implementation will have a predictable success rate for reasonably large sample sizes and will provide a valuable starting point for estimating success rates of scaled-up approaches.

We considered a “sample” one piece of tissue, whether it comes from the same or different animals. We consider one “batch” the group of samples processed simultaneously in the same run of the automated tissue processor (Supplementary Figure S2). “Experiments” may contain samples arising from the same batch or from different batches. The relationship between samples, batches and experiments is described in detail in Supplementary Table S7.

LXRT Imaging

After staining and embedding, the samples were imaged with a laboratory-source X-ray microscope (LXRT; Zeiss Xradia 510 Versa, 40 kV, 3 W, LE2-4 filter, 1601 projections over 180° of sample rotation, 3–20 s exposure time per frame, one tile, 3–6 μm pixel size). The reconstructed tomogram was used to screen for artefacts.

SXRT Imaging

The SXRT dataset analysed in this work is C525, described previously in (Bosch et al., 2022). In brief, samples were imaged with synchrotron X-ray computed tomography exploiting propagation-based phase contrast at the I13-2 beamline of the Diamond Lightsource (Didcot, United Kingdom; 22 keV, 3001 projections over 180° of sample rotation, 0.4 s exposure time per frame, 2–24 tiles, 325 nm isotropic voxel size). The reconstructed tomogram is referred to as the SXRT dataset.

Trimming and SBEM Imaging

To meet the sample size and conductance requirements of our electron microscopes, the resin embedded samples were further processed prior to electron microscopy. Firstly, trims were planned from rendered LXRT or SXRT datasets in Amira (Thermo Fischer Scientific) or Imaris (Bitplane). 2P datasets were warped into LXRT or SXRT space and overlaid on top to help target and preserve the functionally imaged volume in the trimmed block, as described in Bosch et al. (2022). With the help of this trimming plan, embedded samples were trimmed to a suitable shape and size with a razor blade and a glass and/or diamond knife (Diatome) at an ultramicrotome (Leica UC7). Subsequently, the sides of the sample were painted with silver colloidal suspension (SPI supplies) and sputter-coated with a layer of 10 nm-thick platinum (Quorum Q150R S).

2D images used for ultrastructure checks (Supplementary Figure S4A) were obtained in a 3View2—Zeiss Merlin SBEM under high vacuum (1.5 kV, 0.2 nA, 2 μs/pixel, 10 nm or 50 nm pixel size) using an OnPoint (Gatan) back-scattered electron detector.

The volume SBEM dataset analysed in this study comes from sample C525a, described previously in Bosch et al. (2022). This volume SBEM dataset was obtained with a 3View2—variable pressure Zeiss Sigma SBEM using focal charge compensation (2.5 kV, 30 μm aperture, (80 nm)3 voxel size).

Artefact Categorisation

After looking through all LXRT datasets, 12 types of artefacts were defined according to their appearance in the LXRT datasets (Supplementary Table S1). For each sample, a score of 0 (absent) or 1 (present) was given to each artefact type and the following variables compiled: “experiment ID,” “experiment start date,” “fixation method,” “fixative formula,” “staining method,” “anaesthetics,” “culling method,” “use of dexamethasone,” “use of cold buffer,” “mouse line,” “tissue type,” “slice index,” “thickness” and “2P history”. For overall occurrence rate quantification, all samples from all experiments were included. Other analyses were done using samples from relevant experiments (Supplementary Table S5). Analyses were performed with Python3 (Van Rossum and Drake, 2009) in a Jupyter notebook (Kluyver et al., 2016). The artefact scoring tables and analysis scripts are available (see Data availability statement).

Sample Volume Estimation

In total, the final volume of 80 samples was measured from LXRT datasets of the resin-embedded samples (Supplementary Table S6). Sample volume was measured as the number of voxels containing tissue times the voxel size. The voxels containing tissue data were masked in Fiji (Schindelin et al., 2012) by thresholding the signal and discarding low-signal background voxels. For samples with “crack” artefact, the volume occupied by “crack” was excluded from the tissue voxel count.

To estimate the volume of the fresh tissue samples, we utilised landmark pairs from the set of “warping landmarks” used to calculate the warp field (see results and below). Paired landmarks mark the same features of a sample in both datasets. Changes in tissue volume after sample preparation will therefore be reflected in changes in distances between landmarks. For one sample, 20 landmarks from the 2P dataset were randomly selected and their distances to every other landmark calculated. The top 400 distances were taken (Df1_400). Distances between these landmarks were also calculated for the structural dataset of the same sample (Ds1_400). The amount of sample volume change (r) was taken to be the averaged ratio between Ds1_400 and Df1_400 cubed. The fresh sample volume is therefore the CT volume divided by r (Supplementary Figure S5).

Tracing, Warping and Warping Accuracy Estimation

Dataset annotations, including somata seeds and blood vessel tracings, were generated with the data storage, display and annotation tool webKnossos (Boergens et al., 2017).

Warping landmarks were seeded manually using the Fiji plugin BigWarp (Bogovic et al., 2016) and then exported. Landmarks were any distinguishable features shared by the datasets involved, such as blood vessel branching points. The warping function was then calculated using the exported landmark file and a function from the webKnossos skeleton analysis toolbox built in the Matlab (Mathworks) environment. This toolbox (Bosch et al., 2022) is available online (see Data availability statement).

For the C525a dataset, the warping function between 2P and SBEM was calculated using 171 landmarks seeded throughout the dataset (Figure 6A). Similarly, 64 and 51 landmarks were used for 2P to SXRT warping and SXRT to SBEM warping, respectively. These sets of paired landmarks were referred to as “warping landmarks” in the results section.

For warping accuracy quantification, blood vessels in the upper EPL within a (160 μm3) bounding box were traced in the 2P, SXRT and SBEM datasets of sample C525a. To trace the same sets of blood vessels in all three datasets, the bounding boxes used were centred on a functionally imaged soma in 2P and its warped location in SXRT and SBEM. Due to the resolution limit, thin blood vessels were not traceable in the SXRT dataset and the corresponding blood vessel tracings in the 2P and SBEM datasets were deleted for consistency. After this, all blood vessels branching points (junctions) were named and matched across the three datasets using blood vessel shape as a guide. Lastly, freely ending blood vessel tracings were trimmed up until a junction. Therefore, the same blood vessels were covered to the same extent in all three datasets and hence could be compared later. Processed tracings containing those annotations in the different datasets are available online (see Data availability statement).

The blood vessel tracings were then warped into a common framework (e.g., SBEM space, Figure 6C). Each tracing was then split into segments at blood vessel junctions. For each segment, the average segment length (among 2P, SBEM and SXRT) was calculated. This length was then divided by an interpolation step size of 500 nm, giving the number of steps (n) to interpolate for this segment. Subsequently, each segment from each dataset was interpolated with n steps, giving n points on the blood vessel that are matched across all three datasets. The tracing accuracy at that location between any two datasets was then estimated to be the distance between the matched points (Figure 6D). The overall warping accuracy between any two datasets was taken to be the mean of all matched-point distances.


Variables Affecting Performance in Staining Soft Tissue

The olfactory bulb (OB) contains the first neuronal circuit in the olfactory sensory pathway, the glomerular column (Mori et al., 1999). Thick (>0.4 mm) dorsal slices of the bulb contain all the histological layers embedding that circuit (Figure 1A). We prepared 415 0.4-1 mm-thick vibratome-cut slices of mouse olfactory bulb for electron microscopy (Figure 1A1), of which 328 were first dorsal slices (Figure 1A2). However, first slices are asymmetrical, with one side containing biological features potentially capable of acting as a barrier for the diffusion of the staining solutions (e.g., meninges, tightly packed axon bundles), which could in turn translate into asymmetrical staining artefacts. To facilitate the identification of specific artefact patterns, we also prepared 87 second to fourth dorsal slices (Figure 1A3) which, while also containing all histological layers, present a symmetrical structure that allows free diffusion of staining solutions through their shortest dimension. Additionally, to account for possible effects only attributable to larger specimen sizes, we prepared 8 and 4 large samples consisting of complete and half olfactory bulbs, respectively, adding up to a total of 427 samples.


FIGURE 1. Categorisation of staining artefacts from mouse olfactory bulb samples. (A1) Positions of the first and second dorsal slices taken from the mouse olfactory bulb. (A2) LXRT volume of a first dorsal slice and LXRT image of the dashed region in the volume. (A3) same as (A2) but for a second dorsal slice. (B) Structure of an SPP. Staining is split into key steps. Steps that vary between SPP tests are shown in coloured boxes. (C) 12 types of artefacts: name, schematic of appearance, overall occurrence rate in OB samples and LXRT image of a representative sample. (D1) Examples of samples with co-occurring artefacts. They respectively contain “sideways” (white asterisk) and “bubble” (orange); “crack” (yellow) and “smoky” (grey); and “undefined” (purple), “patchy” (blue) and “bubble” (orange). (D2) Number and percentage of OB samples with zero, one or multiple artefacts. Abbreviations: onl: olfactory nerve layer, gl: glomerular layer, epl: external plexiform layer, mcl: mitral cell layer, ipl: internal plexiform layer, gcl: granule cell layer, img: imaging. Os: osmium tetroxide, R: reducing agent e.g., potassium ferrocyanide, TCH: thiocarbohydrazide, U: uranyl acetate, Pb: lead aspartate, arte: artefact. Scale bars are shared across all images in the same panel.

We processed the samples following a standard protocol structure (Briggman et al., 2011; Tapia et al., 2012; Wilke et al., 2013; Hua et al., 2015; Mikula and Denk, 2015; Pallotto et al., 2015; Wanner et al., 2016; Genoud et al., 2018) (Figure 1B). After all in vivo studies had finished on the animal, the brain tissue was fixed, stained with osmium tetroxide, reduced with potassium ferrocyanide, incubated with the bridging agent thiocarbohydrazide, and further stained with osmium tetroxide, uranyl acetate and lead aspartate. Finally, samples were dehydrated in a series of ethanol solutions, embedded in hard epon resin and cured in an oven for 72 h at 60°C. A number of variations were implemented to this protocol backbone with the aim to optimise the yield of samples holding optimal staining quality. Staining quality was routinely assessed with X-ray tomography using a laboratory-based micro-CT (LXRT), a technique well-suited to resolve the different histological layers (Bosch et al., 2022) as well as the dynamics of heavy metal diffusion when staining brain tissue (Ströh et al., 2021). Staining, dehydration and embedding steps were conducted on an automated tissue processor. Thereby staining variability would appear systematic, and large sample groups could be processed simultaneously to address specific staining variables.

Categories of Staining Artefacts

We observed a wide variety of SPP-derived artefacts (Figure 1C). Twelve artefact categories were defined: After examining all stained samples, we described all artefacts into eleven distinct types (Supplementary Table S1). An additional category “undefined” collected the samples (16%, 66/427) that showed poor staining but could not be categorised into any of the defined types. Each sample was assessed independently for the presence of any of the artefacts. Notably, more than one artefacts can co-occur in one sample (Figure 1D). Finally, around 20% of the samples do not contain any artefact and were referred to as “perfect” (Figure 1D2). In our samples, the most frequently observed artefacts were “sideways” and “crack”, which collectively occurred in 60% of all samples. Other insights could be extracted from the prevalence (% of experiments affected) and penetrance (% of affected samples within an experiment) of the artefacts (Supplementary Figures S1, S2), yet we focussed our attention on understanding the aetiology of “sideways” and “crack”, since even modest improvements in the occurrence of these two major artefacts would convey a meaningful improvement in the overall experimental yield.

Diffusion Kinetics and Reduction of Osmium Tetroxide Underly the Most Common “Sideways” Artefact

The “sideways” artefact consists of a poorly stained region which is fully enclosed within the sample and displaced closer towards the brain surface compared to the sliced surface in first slices (Figure 2A). The transition from well stained to poorly stained region may be abrupt, gradual, or both, giving “clear”, “fuzzy”, and “compound boundary sideways”, respectively (Figure 2B). In second slices, where two sliced surfaces make up most of the surface area, the “sideways” artefact occurs at a different anatomical location yet the boundaries remain closer to the brain surface than to the sliced surfaces (Figure 2A). Moreover, in first slices where the brain surface was locally perturbed (because of e.g., some broken meninges exposing the underlying tissue), the distance from that point in the brain surface to the artefact border was similar to the distance between the artefact border and the sliced surface (Figure 2A). Finally, the overall pattern resembles that of a staining artefact previously described (Hua et al., 2015; Mikula and Denk, 2015) (Supplementary Figure S3).


FIGURE 2. Occurrence and proposed mechanism of “sideways” and “central” artefacts. (A) LXRT images of four samples with “sideways” and schematics of hypothesised causes. In the last two samples, the positions of thin GL and broken meninges are highlighted with arrowheads, and part of the “sideways” boundary is marked with a dashed line. In schematics, the size of triangles represent the diffusion rate. (B) Subdivision of “sideways” by boundary type and by complexity. (C) Summary of three SPPs. Staining is split into key steps. Length and colour of staining steps represent incubation time and temperature, respectively. Double dash on box represents overnight incubation. Thin lines connecting staining reagents represent water washes. (D1) Occurrence rate of “sideways” with “clear” and “fuzzy” boundaries in the three SPPs. (D2) Occurrence rate of “simple” or “compound” “sideways” in the three SPPs. (D3) Penetrance of “sideways” in experiments processed with the three SPPs. Each dot represents one experiment. (E) Samples with (2nd and 4th, 1 mm-thick) and without (1st and 3rd, 0.6 mm-thick) the “central” artefact. (F) Occurrence of “central” depends on slice nominal thickness. Each data point represents the % of samples with “central” artefact within an experimental batch (all samples stained with 2sRO protocols). Error bar shows mean ± standard error of mean. Abbreviations: GL: glomerular layer, imm: immersion, per: perfusion, dehy: dehydration, emb: embedding, Os: osmium tetroxide, R: potassium ferrocyanide, TCH: thiocarbohydrazide, U: uranyl acetate, Pb: lead aspartate. Scale bars are shared across all images in the same panel.

Heavy metals passively diffuse through soft tissue following a quadratic kinetic (Mikula and Denk, 2015; Ströh et al., 2021), i.e., staining three times the depth takes nine times longer. From our observations, we hypothesised that the brain surface is less permeable than the sliced surface to the diffusion of heavy metals. If so, reagents would reach locations through the sliced (permeable) border more easily. This combined with insufficient incubation times could lead to the “sideways” artefacts we observed.

In order to assess whether diffusion time affected the presence of this artefact, we modified the protocol so the duration of all staining steps was doubled (Figure 2C). This led to a slight change in the pattern of the “sideways” artefact — with it now displaying only clear boundaries (Figures 2D1,D2).

A second mechanism that could contribute to the “sideways” artefact is the reduction of osmium VIII to osmium IV inside the tissue during staining (Hua et al., 2015). A prior study showed mm-thick mouse brain tissue samples could be evenly stained for SBEM imaging if the reducing agent was presented after the tissue had been already stained with osmium tetroxide (instead of presenting the fresh tissue with reduced osmium tetroxide). We refer to this protocol as “2-step reduced osmium” (2sRO) (Figure 2C). 2sRO also differed in that the incubation of the bridging agent TCH was milder (in time and in temperature), the uranyl acetate incubation (typically overnight at 4°C) was extended with a short (2 h) hot (50°C) period and the incubation of lead aspartate was milder (in temperature). Finally, this protocol employed a distinct fixation method — perfusion — to the one used before — immersion. We implemented the complete 2sRO to replicate its reported performance and compared it across our benchmarked protocol (1sRO). Indeed, experiments stained with 2sRO displayed a distinct penetrance of the “sideways” artefact, now more variable and displaying extreme values – such as being entirely absent in all samples of 38% of the experiments (6/16) (Figure 2D3).

The diffusion of staining solutions through the specimen was also found associated with another artefact type: “central”. “Central” displays a poorly stained region fully enclosed inside the sample, but in this case the unstained region presents virtually no metal staining — with no histological structure revealed by LXRT, and its location appears equally distant to all sample borders (Figure 2E). It is similar to previously reported artefacts associated with insufficient staining times (Hua et al., 2015; Mikula and Denk, 2015; Ströh et al., 2021), and it was completely absent in slices with thickness smaller or equal than 0.6 mm (Figure 2F).

Altogether, we identified two possible mechanisms behind the most prevalent artefact “sideways”: diffusion kinetics of the staining solutions and reduction of osmium tetroxide during the first staining step, and the implemented changes removed this artefact from all samples in 38% of the batches.

“Cracks” can be Reduced With Glutaraldehyde Fixation and Gentle TCH Reactions

The second most prevalent artefact was “crack” (Figure 1C), whose presence precludes any further analysis on the affected sample: if this artefact occurs in the volume of interest, samples become unusable as cellular features cannot be annotated across even small cracks. In fact, even if “crack” artefacts were found outside of the volume of interest, the sample should also be treated as unusable because smaller cracks were likely to be found across the rest of the volume. We further classified the observed cracks into two non-mutually-exclusive subtypes, based on whether the crack did not follow any histological feature (“crack I”) or whether it followed histological borders (“crack II”) (Figure 3A).


FIGURE 3. Occurrence and proposed mechanisms of “crack” subtypes. (A1) LXRT images of samples with “crack” subtypes and schematics of their hypothesised causes. In schematics, “layer a” and “layer b” represent two tissue layers and the dashed box represents the starting tissue volume i.e., when fresh. Dots represent gas bubbles. Arrows represent expansive forces experienced by the tissue. Top, 0.6 mm-thick representative samples from experiments A, B and C referred to in (B1). Bottom, samples showing “crack II” at GL-EPL boundary, MCL, and GCL, respectively. (A2) Number and percentage of “crack” samples with only “crack I” or “crack II” or with both “crack” subtypes. (A3) Percentage of “crack II” samples with the artefact occurring at GL-EPL boundary, MCL or GCL. (B1) “Crack” occurrence in samples from four experiments with different TCH incubation conditions and GA concentrations. (B2) “Crack” occurrence in 0.6 mm-thick samples fixed with different fixatives. Empty squares indicate non-evaluated cases. (C) Volume of samples that do and do not contain “crack”. Each dot represents one sample. The boxes extend from the first to the third quartile, with the line representing the median. The whiskers extend to the most extreme, non-outlier data points. (*: p < 0.5, **: p < 0.005, ***: p < 0.0005, two-tailed paired t-test). (D1) Experimental design of a step-by-step investigation of artefact occurrence. Double dash on box represents overnight incubation. Length and colour of staining steps represent incubation time and temperature, respectively. (D2) LXRT images of representative samples from each experimental group described in (D1). (E) Occurrence of “crack I” (grey line), “crack II” (blue line) and “crack” (blue bars) in each experimental group described in (D1). (F1) Sample volumes of embedded samples of all experimental groups described in (D1), measured in LXRT datasets. Each light dot represents one sample. Each dark dot represents the average of all samples from an experimental group. For the “fresh” group (triangle), the sample volume was estimated from warping data and the data point is connected to the next staining step with a dashed line. (F2) Ratio of the mean sample volume of each experimental group to that of the previous group. A dashed line marks ratio = 1 i.e., no volume change. Abbreviations: GL: glomerular layer, EPL: external plexiform layer, MCL: mitral cell layer, GCL: granule cell layer, w: with, GA: glutaraldehyde, TCH: thiocarbohydrazide, PFA: paraformaldehyde, imm: immersion, Os: osmium tetroxide, R: potassium ferrocyanide, TCH: thiocarbohydrazide, U: uranyl acetate, Pb: lead aspartate, dehy: dehydration, EtOH: ethanol. Scale bars are shared across all images in the same panel.

“Crack” might be generated by expansive forces overwhelming the forces that keep the tissue together (Figure 3A1). Therefore, we expected chemical fixation to play a role in preventing the appearance of cracks — by making the tissue less fragile. Paraformaldehyde and glutaraldehyde are the two most commonly used fixatives for electron microscopy (EM), employed either alone or in combination (Morris, 1965). Glutaraldehyde provides a stronger fixation, derived from its two aldehyde groups that enable reactions cross-linking protein molecules (Kiernan, 2000). We compared sample batches fixed with commonly used fixative solutions employing distinct concentrations of glutaraldehyde. On the other hand, it has been reported that the bridging agent TCH generates gaseous nitrogen as a subproduct when binding to the fixed osmium (Guha and De, 1924; Seligman et al., 1966; Mikula and Denk, 2015). We assumed the stress generated would be related to the amount of gas created, which in turn should correlate with the amount of TCH-osmium reactions. We thereby devised two stress scenarios, differing in the duration and temperature of the TCH step. We observed that fewer samples carried “crack” artefact when processed with the milder TCH incubation, even when fixed with the milder fixative solution (Figure 3B1), and that crack severity decreased at the same time (Figure 3A1). Moreover, glutaraldehyde was required in the fixative solution in order to prevent the appearance of cracks (Figure 3B2). Finally, cracks of all types were more frequently present in larger samples (Figure 3C).

We then aimed to identify at which point in the staining protocol “crack” appeared by stopping the staining process at every single step (Figure 3D). As expected, “crack” started appearing during the TCH step, and their presence increased progressively throughout the consecutive staining steps (Figure 3E). Moreover, this experiment showed that samples expand during the staining protocol (Figure 3D2), that expansion starts during the TCH step and that samples continue to expand during the following osmium tetroxide step (Figure 3F).

Altogether, “crack” artefact can be prevented by including glutaraldehyde in the fixative and by minimising the TCH reaction.

The Effect of Prior in vivo 2P Imaging on Staining Success

Next, we evaluated the effect of prior in vivo 2P imaging on the experimental yield of the SPPs (Figure 4). The throughput (the number of samples processed per batch) was severely reduced by the addition of the in vivo imaging step (Figure 4A). Nevertheless, on average the success rate of the experiment remained unaltered by the incorporation of this imaging modality (Figure 4B1), despite introducing some 2P-specific artefacts (Figures 4B3,B4). Edema is considered a source of surgery-derived histology-scale damage (Betz and Coester, 1990; Judkewitz et al., 2009). To minimise this factor, we introduced a series of edema-preventing measures before and during surgery: mice were administered with the anti-inflammatory drug dexamethasone before the experiment, ice-cold buffer was frequently applied to cool down the drilling of the skull and laser exposure during imaging was minimised whenever possible. Nevertheless, we did not observe a significant effect of these additional measures (Figure 4C).


FIGURE 4. The effect of 2P and dexamethasone on artefact occurrence. (A) Number of samples with (“2P”) and without in vivo functional imaging (“no 2P”) from each experiment. Each light dot represents “2P” samples or “no 2P” samples of one experiment. Solid lines connect “2P” samples to “no 2P”samples from the same experiment. Black dots represent the mean. Dashed lines connect the mean values. Asterisks represent statistical significance (***: p < 0.0005, two-tailed paired t-test). (B1) Percentage of perfect samples from each experiment. (B2–B4) Percentage of samples with “sideways” (B2), “necrotic” (B3) and “sandy” (B4). (C) Artefacts in samples administered with and without dexamethasone within the same experiment (n = 13 and 14, respectively).

In sum, we optimised the SPP to obtain good quality images from all modalities with sufficient yield. At this point we turned our attention towards the second main factor affecting the throughput of the CMI workflow: the quality of the registration of the different modalities into a multimodal map.

Warping Accuracy Across Imaging Modalities

In order to reliably identify neurons functionally-recorded in vivo in the anatomical map acquired ex vivo, the datasets arising from each imaging modality need to be spatially correlated at cellular precision (Helmstaedter, 2013). By identifying conserved pairs of landmarks in two modalities reporting the same space, a warp field relating the two spaces can be computed (Figure 5A1) (Bogovic et al., 2016). This set of landmarks was termed “warping landmarks”. This warp field can later be used to translate any spatial annotation across the two coordinate systems (Figure 5A2). Blood vessels can be used to set warping landmarks for correlating in vivo 2P and SBEM datasets, since their structures can be reliably retrieved by both modalities. They often describe chiral spatial patterns that occur at a density potentially sufficient to map tissue at single cell body precision (Bock et al., 2011; Briggman et al., 2011; Lee et al., 2016; Bosch et al., 2020; Miettinen et al., 2021; Bosch et al., 2022). As blood vessels are also resolved by SXRT, they can be used as a source for landmarks for warping across the three dataset modalities (Figure 5B).


FIGURE 5. Warping using common features e.g., blood vessels and cell correlation by annotation warping. (A1) Common features found throughout the volume of two datasets (i.e., paired landmarks) are used to calculate a warping function (“W”). (A2) The warping function can be used to warp any annotations from one dataset to another. The inputs needed are the x, y, z coordinates of the annotations in the native dataset space and the outputs are the x, y, z coordinates of those annotations in the target dataset space. Annotation warping is used to correlate ROIs from two datasets. ROIs can be correlated as long as the warped ROI falls to a position that is closer to the target ROI than any other ROI. For (A1,A2), cubes represent datasets (e.g., 2P, LXRT, SXRT, SBEM), acquired from the same piece of tissue. The datasets may differ in the dimensions, voxel sizes and the exact region of the sample, differences here represented by cube dimensions, landmark sizes and landmark positions, respectively. Other differences may involve a reversed order in one of the spatial dimensions (such as an horizontal flip in all images), making it necessary to use chiral structures as a source for warping landmarks. (B) Example images showing different blood vessels (arrowheads) resolved by 2P, SXRT and EM datasets.

We aimed to quantify the precision at which one could warp datasets obtained by those three modalities in a previously established correlative workflow (Figures 6A,B) (Bosch et al., 2022). For this purpose, we manually traced blood vessels in the same region spanning (160 µm3) in each original dataset (Figure 6A). Then, these tracings were warped into a common reference space using the warp field calculated with the previously described set of “warping landmarks” (Figure 6C). Blood vessel branching points are unique junctions and many of them (n = 39) could be identified across all imaging modalities, defining blood vessel segments (n = 35). We interpolated a defined number of points within each blood vessel segment in all modalities and paired the resulting nodes (Figure 6D). This allowed for calculating the distance across modalities between paired nodes of the blood vessel (Figure 6E). Warping accuracy depends on the accuracy at which landmarks are resolved across the dataset and therefore its value can vary throughout the dataset’s topology. We observed that the range of warping accuracies recorded depended on the pair of imaging modalities being warped: it ranged between 0.6 and 27.3 µm when warping 2P to EM or to SXRT and between 0.2 and 9.9 µm when warping SBEM to SXRT (Figures 6E,F). On average, any point in the blood vessel could be mapped at ∼14 µm accuracy between the 2P and the SBEM dataset and at ∼9 µm accuracy between the 2P and the SXRT dataset, and regions between SXRT and SBEM could be mapped at an accuracy of ∼2 µm (Figures 6F,G), implying that any region can be mapped across any of the three imaging modalities at single cell precision (Figure 6H).


FIGURE 6. Accurate warping allows for cell correlation among 2P, SXRT and SBEM. (A) Position of landmarks used to generate warping functions between two datasets and position of blood vessel tracings used to quantify warping accuracy in 2P, SXRT and SBEM datasets of a sample. Blue, black and red dots represent landmarks used between 2P and SXRT datasets, between 2P and SBEM datasets and between SXRT and SBEM datasets, respectively. (B) Fields of view of 2P (red, warped), SBEM (yellow, warped) and SXRT (blue, native space) datasets. (C) Same blood vessels traced in the 2P (red, warped), SXRT (blue, warped) and SBEM (yellow, native space) datasets. (D) Quantification of warping accuracy. The nodes in the blood vessel traced in one imaging modality (red) are paired to the equivalent nodes in the other modality (yellow). The distance between two paired nodes (e.g., “d1”) reflects the warping accuracy at that location. Large dots represent nodes in the original tracing. Lines represent paths of the original tracing. Small dots represent nodes interpolated along the tracing path. (E1,E2) Native SBEM blood vessel tracing coloured by the mapped 2P-SXRT and 2P-SBEM warping accuracy, respectively. (F) Distribution of warping accuracies (i.e., distances between paired nodes). Purple, orange and green represent 2P-SXRT, 2P-SBEM and SXRT-SBEM accuracies, respectively. Mean warping accuracies are marked by triangles. (G) Mean warping accuracy between each pair of imaging modalities. (H) Example mitral cell and surrounding blood vessels warped from 2P (red) and SXRT (blue) datasets into SBEM (yellow) space. The diameter of the sphere is 20 µm, similar to the diameter of a mitral cell soma.


We have identified and addressed two major bottlenecks in the performance of a correlative multimodal imaging pipeline capable of resolving the structure and function of a neural circuit in the mouse olfactory bulb: the yield of sample preparation and the accuracy of dataset correlation across 2P, SXRT and SBEM modalities.

Our study focuses on mouse olfactory bulb (OB) dorsal slices. The OB contains the first neuronal circuit in the olfactory sensory pathway (Mori et al., 1999), the glomerular column. The architecture of these columns is compact (<1 mm3) and modular (>1000 columns in each OB hemisphere) (Schwarz et al., 2018) and their anatomical location makes them accessible for in vivo studies involving minimally invasive surgeries (e.g., imaging intracellular [Ca2+] through a cranial window). Combined, these properties make the OB a good model for studying systems neuroscience in the mammalian brain (Schaefer and Margrie, 2007). For in vivo experiments aimed at recording the circuit’s function, 2P datasets can map a dorsal field of view of ∼mm2, but depth resolution of 2P is poorer than the lateral one (Chaigneau et al., 2011). Moreover, its signal-to-noise ratio decreases with the depth of the region being imaged, making it increasingly challenging to retrieve a signal beyond a depth of 0.5 mm deep inside the brain (Wachowiak et al., 2013; Ackels et al., 2020).

Our horizontal slices provide a maximal overlap between the 2P-imaged region and the stained tissue that will later be imaged at higher structural resolution with SXRT and SBEM. However, these dorsal sections present anatomical asymmetries: one of their sides is a vibratome-cut surface, whereas the other side is the anatomically-preserved surface of the brain. We generated a set of samples that would present a more symmetrical disposition to diffusing solutions: second slices from the same brain region.

Two main artefacts dominated the sample preparation protocol (SPP). An asymmetrical unstained region here denoted “sideways” was the most frequent artefact type. Its pattern was reminiscent of an artefact defined previously (Hua et al., 2015; Mikula and Denk, 2015) although in our samples it displayed an asymmetrical morphology (Supplementary Figure S3). We found it to be related to two independent parameters: the diffusion of heavy metals through the tissue and the reduction of the first osmium tetroxide staining step. This is consistent with its symmetrical appearance in previous reports that analysed punch-dissected tissue cylinders presenting a symmetrical structure to freely diffusing solutions (Hua et al., 2015). Addressing those two factors allowed for removing its presence in all samples in a subset of experimental batches.

The next most common artefact was the presence of cracks in the stained tissue, possibly generated by an imbalance between expansive and cohesive forces in the tissue during sample preparation. It is known that osmium tetroxide can provoke both the expansion and the shrinkage of brain tissue samples during staining (Ströh et al., 2021) and that the bridging agent TCH generates gaseous nitrogen as a by-product of binding to fixed osmium (Guha and De, 1924; Seligman et al., 1966; Mikula and Denk, 2015). Consistent with that, we found that a milder TCH reaction can minimise the generation of cracks in samples extending multiple mm3 in volume, while not compromising the overall staining quality. On the other hand, strong sample fixation with glutaraldehyde was necessary to prevent tissue from cracking during the staining process. Finally, larger samples displayed cracks more frequently. While dehydration can lead to severe shrinkage of fresh soft tissues (Buytaert et al., 2014) we observed tissue had overall expanded after dehydration. Altogether, these findings allowed extracting simple specific requirements for future protocols to minimise the appearance of cracks: 1) employing fixatives containing at least 1.25% glutaraldehyde, 2) reducing the strength of the TCH reaction, and 3) minimising the smallest dimension of the sample being processed whenever possible. Further, albeit more difficult to implement improvements are: 1) finding optimal formulations for fixative solutions for a specific sample, 2) replacing TCH for a reagent that does not create bubbles — such as pyrogallol (Mikula and Denk, 2015; Song et al., 2022), and 3) minimising the sample expansion processes by e.g., counteracting osmium-derived effects with suitable additives (Hua et al., 2015; Ströh et al., 2021).

A subset of the experiments reported were obtained by chemically fixing the tissue by immersion in fixative (Supplementary Table S2). All artefacts in immersion-fixed samples could be diagnosed by observing an LXRT dataset covering the entire sample. Such 3D datasets can be obtained at a rate of multiple samples per day with benchtop LXRTs (Metscher, 2009; Bushong et al., 2015; Gutiérrez et al., 2018) and of nearly one sample per minute at specialised hard X-ray microtomography beamlines at synchrotron X-ray facilities (Walker et al., 2014; Rau, 2017; Harkiolaki et al., 2018; Walsh et al., 2021; Bosch et al., 2022). A distinct approach could involve fixing the samples via intracardiac perfusion of the fixative (Gage et al., 2012) (Supplementary Figure S4). Doing so could prevent some artefacts that might be specific to immersion-fixation (Supplementary Figure S4A) and would remove the requirement of a fast dissection of the tissue of interest, making this approach compatible with processing some regions that might be out of reach or processing multiple different regions in the same animal. On the other hand, maintaining ultrastructural integrity is generally more challenging with perfusion (Palay et al., 1962; Crang et al., 1988; Maunsbach and Afzelius, 1999; Tao-Cheng et al., 2007; Glauert and Lewis, 2014) (Supplementary Figure S4B). It is worth noting that evaluating sample quality by its ultrastructure has a lower throughput than relying on LXRT for that purpose. Therefore, choosing perfusion as a fixation method could create a quality control bottleneck in the experimental pipeline (Supplementary Figure S4C).

The yield of perfectly stained samples in an experiment was not altered when extending the pipeline with prior in vivo 2P imaging. However, those experiments have a much lower throughput (Figure 4A). A low throughput combined with a moderate yield can make it challenging to predict future stock of samples with both optimal 2P imaging and optimal staining. Importantly, sample preparation does not end with staining: trimming the specimen will be necessary to accommodate its dimensions to some acquisition techniques (Supplementary Figure S4C) (Karreman et al., 2016; Ronchi et al., 2021; Bosch et al., 2022), hard X-ray tomography imaging may inflict damage from radiation dose and temperature changes (Henderson, 1995; Howells et al., 2009; Du and Jacobsen, 2018), volume electron microscopy techniques imply slicing or destruction of the sample and therefore add their own risk of a specific sample becoming irreversibly damaged, and optimal targeting of the features of interest to be imaged with the distinct modalities precludes extracting the information sought by the initial experimental question (Bosch et al., 2020; Walter et al., 2020). In light of these successive challenges, it is imperative for a successful CMI experiment to prepare a sufficient number of samples at the first stage. Given the low throughput nature of in vivo 2P experiments (in particular in e.g., awake or even behaving animals), this can currently partially be achieved by e.g., obtaining samples from multiple brain regions in a given experiment or by coordinated efforts of multiple researchers or even research labs. Notably, any improvement to the SPP will reduce this initial experimental burden and allow for increasingly complex physiological analysis to be combined with structural investigations. The approach described in this study highlights the major bottlenecks in a protocol, but it does not provide an exhaustive readout of all sources of variability, nor does it describe the spectrum of ultrastructural defects one can encounter in optimally-stained specimens as assessed by LXRT. Therefore, it may be useful in the future to both reassess the performance of novel protocols with higher throughputs and efficiencies — to find the causes for rarer staining deficits — and to conduct focussed analyses on the quality of ultrastructure. This iterative process can guide the development of more efficient and robust protocols, able to stain larger specimens reliably.

Finally, we demonstrated that tissue locations imaged in vivo can be recalled at single-cell precision in an SXRT dataset. Because SXRT datasets can be obtained on multi-mm3 specimens, this strategic improvement provides an early readout of circuit structure and can inform on how to precisely trim precious specimens, preserving virtually all 2P-imaged cells in a trimmed specimen for follow-up higher resolution structural imaging techniques that impose limitations on sample dimensions and require extraordinary sample quality control standards. Notably, this approach may also be useful for experiments requiring a correlation accuracy better than 10 µm: in those cases, it would narrow the landscape containing the feature of interest from multiple mm3 down to a volume of ∼(10 µm)3. This would relax the requirements imposed to a further disambiguation algorithm, which would operate in the range between (10 µm)3 and the required accuracy. Direct 2P-SXRT correlation therefore generates experimental compatibility of in vivo 2P with not only the different implementations of volume EM but also with other synchrotron hard X-ray modalities such as nanoholotomography (Kuan et al., 2020), laminography (Helfen et al., 2013; Witte et al., 2020) or ptychography (Rodenburg et al., 2007; Holler et al., 2017; Shahmoradian et al., 2017). Adding non-destructive techniques to the workflow not only enables more diverse insights to be recalled from the same sample, but also makes it more resilient, ensuring insights are obtained from a larger yield of samples. Additionally, combining multimodal insights from a single sample may allow developing more robust image processing algorithms, by providing insights free from technique-specific artefacts (e.g., using blood vessels from the X-ray CT datasets, with isotropic resolution, to guide the fine alignment of volume EM images). Altogether, we are confident that improved staining and reliable warping will enable exhaustive studies retrieving physiology, circuit structure and synaptic signatures in targeted regions of slices covering multiple mm3.

Our optimisation was aimed at obtaining a single sample preparation protocol that is compatible with both SXRT and SBEM, since staining for the latter did provide a sufficient signal-to-noise ratio in the former. However, it is possible that better SXRT datasets could be obtained if the samples were prepared following specially SXRT-optimised preparation protocols. Those could, for example, aim to maximise the phase contrast in the detected X-ray beam while minimising the absorption and thereby bring benefits in terms of both image contrast and specimen size. Doing so could be envisaged by using weaker metal stainings. This could modify, in turn, the landscape of compatible techniques for future multimodal imaging approaches in systems neuroscience and biological soft tissue research more generally, e.g., enable to extract insights on molecular composition and transcriptomic expression in the tissue alongside its function and structure (Ståhl Patrik et al., 2016; Domart et al., 2020; Close et al., 2021).

Altogether, we report improvements in two limiting steps of a correlative multimodal imaging pipeline for systems neuroscience in the mouse olfactory bulb in sample preparation and in dataset correlation accuracy. These insights will allow us to make informed decisions to troubleshoot and improve the yield of the existing protocols and to design novel pipelines relying solely on hard X-rays to retrieve the location of the cell bodies of interest. These frameworks will pave the way for the design of more resilient workflows to interrogate both structure and function of neural circuits in the mammalian brain.

Data Availability Statement

The datasets presented in this study can be found in the online repositories:,

Ethics Statement

The animal study was reviewed and approved by Ethics Committee of the board of the Francis Crick Institute and the United Kingdom Home Office under the Animals (Scientific Procedures) Act 1986.

Author Contributions

YZ, ATS and CB conceptualised the research project and designed experiments; YZ, TA, AP, M-CZ, AB, and CB performed experiments and acquired data; YZ, ATS and CB analysed and interpreted the data; YZ, ATS and CB wrote the manuscript draft; YZ and CB generated the figures, and all authors reviewed and provided input to the manuscript.


This research was funded in whole, or in part, by the Wellcome Trust (FC001153). For the purpose of Open Access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission. This work was carried out with the support of Diamond Light Source, beamline I13-2 (proposal 20274) and the TOMCAT beamline of the Swiss Light Source at the Paul Scherrer Institut (proposal 20190417). This work was supported by the Francis Crick Institute, which receives its core funding from Cancer Research UK (FC001153), the UK Medical Research Council (FC001153), and the Wellcome Trust (FC001153); by the UK Medical Research Council (grant reference MC_UP_1202/5) to ATS, by a Physics of Life grant to ATS and AP (grant reference EP/W024292/1) and a DFG postdoctoral fellowship to TA. ATS is a Wellcome Trust Investigator (110174/Z/15/Z). AP acknowledges funding from the European Research Council under the European Union’s Horizon 2020 Research and Innovation Programme (grant No. 852455).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.


We are grateful to the biological research and electron microscopy science technology platforms of the Francis Crick Institute. We thank Marta Pallotto, Kara Fulton and Kevin Briggman for support and critical discussions on the initial configuration of the staining protocol, to Manuel Berning for support in developing the warpAnnotations toolbox, to Norman Rzepka for support in the webknossos infrastructure, and to the members of the sample preparation working group of the volumeEM community for insightful discussions on the topic.

Supplementary Material

The Supplementary Material for this article can be found online at:


Ackels, T., Erskine, A., Dasgupta, D., Marin, A. C., Warner, T. P. A., Tootoonian, S., et al. (2021). Fast Odour Dynamics Are Encoded in the Olfactory System and Guide Behaviour. Nature 593 (7860), 558–563. doi:10.1038/s41586-021-03514-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ackels, T., Jordan, R., Schaefer, A. T., and Fukunaga, I. (2020). Respiration-Locking of Olfactory Receptor and Projection Neurons in the Mouse Olfactory Bulb and its Modulation by Brain State. Front. Cell. Neurosci. 14, 220. doi:10.3389/fncel.2020.00220

PubMed Abstract | CrossRef Full Text | Google Scholar

Betz, A. L., and Coester, H. C. (1990). Effect of Steroids on Edema and Sodium Uptake of the Brain During Focal Ischemia in Rats. Stroke 21 (8), 1199–1204. doi:10.1161/01.STR.21.8.1199

PubMed Abstract | CrossRef Full Text | Google Scholar

Bock, D. D., Lee, W.-C. A., Kerlin, A. M., Andermann, M. L., Hood, G., Wetzel, A. W., et al. (2011). Network Anatomy and In Vivo Physiology of Visual Cortical Neurons. Nature 471 (7337), 177–182. doi:10.1038/nature09802

PubMed Abstract | CrossRef Full Text | Google Scholar

Boergens, K. M., Berning, M., Bocklisch, T., Bräunlein, D., Drawitsch, F., Frohnhofen, J., et al. (2017). webKnossos: Efficient Online 3D Data Annotation for Connectomics. Nat. Methods 14 (7), 691–694. doi:10.1038/nmeth.4331

PubMed Abstract | CrossRef Full Text | Google Scholar

Bogovic, J. A., Hanslovsky, P., Wong, A., and Saalfeld, S. (2016). “Robust Registration of Calcium Images by Learned Contrast Synthesis,” in 2016 IEEE 13th International Symposium on Biomedical Imaging (ISBI) (IEEE), 1123–1126. doi:10.1109/isbi.2016.7493463

CrossRef Full Text | Google Scholar

Bosch, C., Ackels, T., Pacureanu, A., Zhang, Y., Peddie, C. J., Berning, M., et al. (2022). Functional and Multiscale 3D Structural Investigation of Brain Tissue through Correlative In Vivo Physiology, Synchrotron Micro-Tomography and Volume Electron Microscopy. Nat. Commun. 13 (1), 2923. doi:10.1038/s41467-022-30199-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Bosch, C., Schaefer, A. T., Pacureanu, A., Walter, A., Karreman, M. A., and Dykas, M. (2020). “Correlated Multimodal Imaging Across Scales,” in Imaging Modalities in Biological and Preclinical Research: A Compendium. Editor A. Walter (Bristol, United Kingdom: IoP-IPEM), III.3-1–III.3-11.

Google Scholar

BRAIN Initiative Cell Census Network (BICCN) (2021). A Multimodal Cell Census and Atlas of the Mammalian Primary Motor Cortex. Nature 598 (7879), 86–102. doi:10.1038/s41586-021-03950-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Briggman, K. L., Helmstaedter, M., and Denk, W. (2011). Wiring Specificity in the Direction-Selectivity Circuit of the Retina. Nature 471 (7337), 183–188. doi:10.1038/nature09818

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckley, I. K., and Porter, K. R. (1967). Cytoplasmic Fibrils in Living Cultured Cells. Protoplasma 64 (4), 349–380. doi:10.1007/BF01666538

PubMed Abstract | CrossRef Full Text | Google Scholar

Bushong, E. A., Johnson, D. D., Kim, K.-Y., Terada, M., Hatori, M., Peltier, S. T., et al. (2015). X-Ray Microscopy as an Approach to Increasing Accuracy and Efficiency of Serial Block-Face Imaging for Correlated Light and Electron Microscopy of Biological Specimens. Microsc. Microanal. 21 (1), 231–238. doi:10.1017/S1431927614013579

PubMed Abstract | CrossRef Full Text | Google Scholar

Buytaert, J., Goyens, J., De Greef, D., Aerts, P., and Dirckx, J. (2014). Volume Shrinkage of Bone, Brain and Muscle Tissue in Sample Preparation for Micro-CT and Light Sheet Fluorescence Microscopy (LSFM). Microsc. Microanal. 20 (4), 1208–1217. doi:10.1017/S1431927614001329

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaigneau, E., Wright, A. J., Poland, S. P., Girkin, J. M., and Silver, R. A. (2011). Impact of Wavefront Distortion and Scattering on 2-photon Microscopy in Mammalian Brain Tissue. Opt. Express 19 (23), 22755–22774. doi:10.1364/OE.19.022755

PubMed Abstract | CrossRef Full Text | Google Scholar

Claudi, F., Petrucco, L., Tyson, A., Branco, T., Margrie, T., and Portugues, R. (2020). BrainGlobe Atlas API: A Common Interface for Neuroanatomical Atlases. Joss 5 (54), 2668. doi:10.21105/joss.02668

CrossRef Full Text | Google Scholar

Close, J. L., Long, B. R., and Zeng, H. (2021). Spatially Resolved Transcriptomics in Neuroscience. Nat. Methods 18 (1), 23–25. doi:10.1038/s41592-020-01040-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Crang, R. F. E., Crang, R. F. E., and Klomparens, K. L. (1988). Artifacts in Biological Electron Microscopy. Springer US.

Google Scholar

Denk, W., Delaney, K. R., Gelperin, A., Kleinfeld, D., Strowbridge, B. W., Tank, D. W., et al. (1994). Anatomical and Functional Imaging of Neurons Using 2-Photon Laser Scanning Microscopy. J. Neurosci. Methods 54 (2), 151–162. doi:10.1016/0165-0270(94)90189-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Denk, W., and Horstmann, H. (2004). Serial Block-Face Scanning Electron Microscopy to Reconstruct Three-Dimensional Tissue Nanostructure. PLoS Biol. 2 (11), e329. doi:10.1371/journal.pbio.0020329

PubMed Abstract | CrossRef Full Text | Google Scholar

Domart, F., Cloetens, P., Roudeau, S., Carmona, A., Verdier, E., Choquet, D., et al. (2020). Correlating STED and Synchrotron XRF Nano-Imaging Unveils Cosegregation of Metals and Cytoskeleton Proteins in Dendrites. eLife 9, e62334. doi:10.7554/eLife.62334

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, M., and Jacobsen, C. (2018). Relative Merits and Limiting Factors for X-Ray and Electron Microscopy of Thick, Hydrated Organic Materials. Ultramicroscopy 184, 293–309. doi:10.1016/j.ultramic.2017.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Dyer, E. L., Gray Roncal, W., Prasad, J. A., Fernandes, H. L., Gürsoy, D., De Andrade, V., et al. (2017). Quantifying Mesoscale Neuroanatomy Using X-Ray Microtomography. Eneuro 4 (5), 1–18. doi:10.1523/ENEURO.0195-17.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Frazier, L. L., and Brunjes, P. C. (1988). Unilateral Odor Deprivation: Early Postnatal Changes in Olfactory Bulb Cell Density and Number. J. Comp. Neurol. 269 (3), 355–370. doi:10.1002/cne.902690304

PubMed Abstract | CrossRef Full Text | Google Scholar

Gage, G. J., Kipke, D. R., and Shain, W. (2012). Whole Animal Perfusion Fixation for Rodents. JoVE (65), e3564. doi:10.3791/3564

PubMed Abstract | CrossRef Full Text | Google Scholar

Genoud, C., Titze, B., Graff-Meyer, A., and Friedrich, R. W. (2018). Fast Homogeneous En Bloc Staining of Large Tissue Samples for Volume Electron Microscopy. Front. Neuroanat. 12, 76. doi:10.3389/fnana.2018.00076

PubMed Abstract | CrossRef Full Text | Google Scholar

Glauert, A. M., and Lewis, P. R. (2014). Biological Specimen Preparation for Transmission Electron Microscopy. Princeton University Press.

Google Scholar

Guha, P. C., and De, S. C. (1924). CXLVII.-A New Method for the Preparation of Thiocarbohydrazide: Mono- and Di-thio-p-urazine. J. Chem. Soc. Trans. 125 (0), 1215–1218. doi:10.1039/CT9242501215

CrossRef Full Text | Google Scholar

Gutiérrez, Y., Ott, D., Töpperwien, M., Salditt, T., and Scherber, C. (2018). X‐Ray Computed Tomography and its Potential in Ecological Research: A Review of Studies and Optimization of Specimen Preparation. Ecol. Evol. 8 (15), 7717–7732. doi:10.1002/ece3.4149

PubMed Abstract | CrossRef Full Text | Google Scholar

Harkiolaki, M., Darrow, M. C., Spink, M. C., Kosior, E., Dent, K., and Duke, E. (2018). Cryo-Soft X-Ray Tomography: Using Soft X-Rays to Explore the Ultrastructure of Whole Cells. Emerg. Top. life Sci. 2 (1), 81–92. doi:10.1042/etls20170086

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayat, M. A. (2000). In Principles and Techniques of Electron Microscopy : Biological Applications. Editor M. A. Hayat (Cambridge: Cambridge University Press).

Google Scholar

Helfen, L., Xu, F., Suhonen, H., Urbanelli, L., Cloetens, P., and Baumbach, T. (2013). Nano-Laminography for Three-Dimensional High-Resolution Imaging of Flat Specimens. J. Inst. 8 (05), C05006. doi:10.1088/1748-0221/8/05/c05006

CrossRef Full Text | Google Scholar

Helmstaedter, M. (2013). Cellular-Resolution Connectomics: Challenges of Dense Neural Circuit Reconstruction. Nat. Methods 10 (6), 501–507. doi:10.1038/nmeth.2476

PubMed Abstract | CrossRef Full Text | Google Scholar

Henderson, R. (1995). The Potential and Limitations of Neutrons, Electrons and X-Rays for Atomic Resolution Microscopy of Unstained Biological Molecules. Quart. Rev. Biophys. 28 (2), 171–193. doi:10.1017/S003358350000305X

CrossRef Full Text | Google Scholar

Holler, M., Guizar-Sicairos, M., Tsai, E. H. R., Dinapoli, R., Müller, E., Bunk, O., et al. (2017). High-Resolution Non-Destructive Three-Dimensional Imaging of Integrated Circuits. Nature 543 (7645), 402–406. doi:10.1038/nature21698

PubMed Abstract | CrossRef Full Text | Google Scholar

Howells, M. R., Beetz, T., Chapman, H. N., Cui, C., Holton, J. M., Jacobsen, C. J., et al. (2009). An Assessment of the Resolution Limitation Due to Radiation-Damage in X-Ray Diffraction Microscopy. J. Electron. Spectrosc. Relat. Phenom. 170 (1-3), 4–12. doi:10.1016/j.elspec.2008.10.008

CrossRef Full Text | Google Scholar

Hua, Y., Laserstein, P., and Helmstaedter, M. (2015). Large-Volume En-Bloc Staining for Electron Microscopy-Based Connectomics. Nat. Commun. 6 (1), 7923–7927. doi:10.1038/ncomms8923

PubMed Abstract | CrossRef Full Text | Google Scholar

Judkewitz, B., Rizzi, M., Kitamura, K., and Häusser, M. (2009). Targeted Single-Cell Electroporation of Mammalian Neurons In Vivo. Nat. Protoc. 4 (6), 862–869. doi:10.1038/nprot.2009.56

PubMed Abstract | CrossRef Full Text | Google Scholar

Karreman, M. A., Mercier, L., Schieber, N. L., Solecki, G., Allio, G., Winkler, F., et al. (2016). Fast and Precise Targeting of Single Tumor Cells In Vivo by Multimodal Correlative Microscopy. J. Cell Sci. 129 (2), 444–456. doi:10.1242/jcs.181842

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiernan, J. A. (2000). Formaldehyde, Formalin, Paraformaldehyde and Glutaraldehyde: What They Are and What They Do. Micros. Today 8 (1), 8–13. doi:10.1017/S1551929500057060

CrossRef Full Text | Google Scholar

Klinger, E., Motta, A., Marr, C., Theis, F. J., and Helmstaedter, M. (2021). Cellular Connectomes as Arbiters of Local Circuit Models in the Cerebral Cortex. Nat. Commun. 12 (1), 2785. doi:10.1038/s41467-021-22856-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Kluyver, T., Ragan-Kelley, B., Pérez, F., Granger, B. E., Bussonnier, M., Frederic, J., et al. (2016). Jupyter Notebooks-A Publishing Format for Reproducible Computational Workflows.

Google Scholar

Kuan, A. T., Phelps, J. S., Thomas, L. A., Nguyen, T. M., Han, J., Chen, C.-L., et al. (2020). Dense Neuronal Reconstruction Through X-Ray Holographic Nano-Tomography. Nat. Neurosci. 23 (12), 1637–1643. doi:10.1038/s41593-020-0704-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubota, Y. (2015). New Developments in Electron Microscopy for Serial Image Acquisition of Neuronal Profiles. Microscopy 64 (1), 27–36. doi:10.1093/jmicro/dfu111

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, W.-C. A., Bonin, V., Reed, M., Graham, B. J., Hood, G., Glattfelder, K., et al. (2016). Anatomy and Function of an Excitatory Network in the Visual Cortex. Nature 532 (7599), 370–374. doi:10.1038/nature17192

PubMed Abstract | CrossRef Full Text | Google Scholar

Maunsbach, A. B., and Afzelius, B. A. (1999). “Fixative Application,” in Biomedical Electron Microscopy. Editors A. B. Maunsbach,, and B. A. Afzelius (San Diego: Academic Press), 79–102. doi:10.1016/b978-012480610-8/50006-8

CrossRef Full Text | Google Scholar

Metscher, B. D. (2009). MicroCT for Comparative Morphology: Simple Staining Methods Allow High-Contrast 3D Imaging of Diverse Non-Mineralized Animal Tissues. BMC Physiol. 9 (1), 11. doi:10.1186/1472-6793-9-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Miettinen, A., Zippo, A. G., Patera, A., Bonnin, A., Shahmoradian, S. H., Biella, G. E. M., et al. (2021). Micrometer-Resolution Reconstruction and Analysis of Whole Mouse Brain Vasculature by Synchrotron-Based Phase-Contrast Tomographic Microscopy. bioRxiv, 2021.2003.2016.435616. doi:10.1101/2021.03.16.435616

CrossRef Full Text | Google Scholar

Mikula, S., and Denk, W. (2015). High-Resolution Whole-Brain Staining for Electron Microscopic Circuit Reconstruction. Nat. Methods 12 (6), 541–546. doi:10.1038/nmeth.3361

PubMed Abstract | CrossRef Full Text | Google Scholar

Mori, K., Nagao, H., and Yoshihara, Y. (1999). The Olfactory Bulb: Coding and Processing of Odor Molecule Information. Science 286 (5440), 711–715. doi:10.1126/science.286.5440.711

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, J. K. (1965). A Formaldehyde Glutaraldehyde Fixative of High Osmolality for Use in Electron Microscopy. J. Cell. Biol. 27, 1A–149A.

Google Scholar

Niedworok, C. J., Brown, A. P. Y., Jorge Cardoso, M., Osten, P., Ourselin, S., Modat, M., et al. (2016). aMAP Is a Validated Pipeline for Registration and Segmentation of High-Resolution Mouse Brain Data. Nat. Commun. 7 (1), 11879. doi:10.1038/ncomms11879

PubMed Abstract | CrossRef Full Text | Google Scholar

Oheim, M., Beaurepaire, E., Chaigneau, E., Mertz, J., and Charpak, S. (2001). Two-Photon Microscopy in Brain Tissue: Parameters Influencing the Imaging Depth. J. Neurosci. Methods 111 (1), 29–37. doi:10.1016/S0165-0270(01)00438-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Palay, S. L., McGee-Russell, S. M., Gordon, S., and Grillo, M. A. (1962). Fixation of Neural Tissues for Electron Microscopy by Perfusion with Solutions of Osmium Tetroxide. J. Cell Biol. 12 (2), 385–410. doi:10.1083/jcb.12.2.385

PubMed Abstract | CrossRef Full Text | Google Scholar

Pallotto, M., Watkins, P. V., Fubara, B., Singer, J. H., and Briggman, K. L. (2015). Extracellular Space Preservation Aids the Connectomic Analysis of Neural Circuits. Elife 4, e08206. doi:10.7554/eLife.08206

PubMed Abstract | CrossRef Full Text | Google Scholar

Panhuber, H., Laing, D. G., Willcox, M. E., Eagleson, G. K., and Pittman, E. A. (1985). The Distribution of the Size and Number of Mitral Cells in the Olfactory Bulb of the Rat. J. Anat. 140 ( Pt 2) (Pt 2Pt 2), 297–308.

PubMed Abstract | Google Scholar

Parrish-Aungst, S., Shipley, M. T., Erdelyi, F., Szabo, G., and Puche, A. C. (2007). Quantitative Analysis of Neuronal Diversity in the Mouse Olfactory Bulb. J. Comp. Neurol. 501 (6), 825–836. doi:10.1002/cne.21205

PubMed Abstract | CrossRef Full Text | Google Scholar

Rau, C. (2017). Imaging with Coherent Synchrotron Radiation: X-Ray Imaging and Coherence Beamline (I13) at Diamond Light Source. Synchrotron Radiat. News 30 (5), 19–25. doi:10.1080/08940886.2017.1364530

CrossRef Full Text | Google Scholar

Rodenburg, J. M., Hurst, A. C., Cullis, A. G., Dobson, B. R., Pfeiffer, F., Bunk, O., et al. (2007). Hard-X-Ray Lensless Imaging of Extended Objects. Phys. Rev. Lett. 98 (3), 034801. doi:10.1103/PhysRevLett.98.034801

PubMed Abstract | CrossRef Full Text | Google Scholar

Ronchi, P., Mizzon, G., Machado, P., D’Imprima, E., Best, B. T., Cassella, L., et al. (2021). High-Precision Targeting Workflow for Volume Electron Microscopy. J. Cell Biol. 220 (9), e202104069. doi:10.1083/jcb.202104069

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaefer, A. T., and Margrie, T. W. (2007). Spatiotemporal Representations in the Olfactory System. Trends Neurosci. 30 (3), 92–100. doi:10.1016/j.tins.2007.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: An Open-Source Platform for Biological-Image Analysis. Nat. Methods 9 (7), 676–682. doi:10.1038/nmeth.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwarz, D., Kollo, M., Bosch, C., Feinauer, C., Whiteley, I., Margrie, T. W., et al. (2018). Architecture of a Mammalian Glomerular Domain Revealed by Novel Volume Electroporation Using Nanoengineered Microelectrodes. Nat. Commun. 9 (1), 183. doi:10.1038/s41467-017-02560-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Seligman, A. M., Wasserkrug, H. L., and Hanker, J. S. (1966). A New Staining Method (Oto) for Enhancing Contrast of Lipid-Containing Membranes and Droplets in Osmium Tetroxide-Fixed Tissue with Osmiophilic Thiocarbohydrazide (Tch). J. Cell. Biol. 30 (2), 424–432. doi:10.1083/jcb.30.2.424

PubMed Abstract | CrossRef Full Text | Google Scholar

Shahmoradian, S. H., Tsai, E. H. R., Diaz, A., Guizar-Sicairos, M., Raabe, J., Spycher, L., et al. (2017). Three-Dimensional Imaging of Biological Tissue by Cryo X-Ray Ptychography. Sci. Rep. 7 (1), 6291. doi:10.1038/s41598-017-05587-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, K., Feng, Z., and Helmstaedter, M. (2022). High-Contrast En-Bloc Staining of Mouse Whole-Brain Samples for EM-Based Connectomics. bioRxiv [Preprint]. doi:10.1101/2022.03.30.486341

CrossRef Full Text | Google Scholar

Ståhl, P. L., Salmén, F., Vickovic, S., Lundmark, A., Navarro, J. F., Magnusson, J., et al. (2016). Visualization and Analysis of Gene Expression in Tissue Sections by Spatial Transcriptomics. Science 353 (6294), 78–82. doi:10.1126/science.aaf2403

PubMed Abstract | CrossRef Full Text | Google Scholar

Ströh, S., Hammerschmith, E. W., Tank, D. W., Seung, H. S., and Wanner, A. A. (2021). In Situ X-Ray Assisted Electron Microscopy Staining for Large Biological Samples. bioRxiv.

Google Scholar

Tao-Cheng, J.-H., Gallant, P. E., Brightman, M. W., Dosemeci, A., and Reese, T. S. (2007). Structural Changes at Synapses after Delayed Perfusion Fixation in Different Regions of the Mouse Brain. J. Comp. Neurol. 501 (5), 731–740. doi:10.1002/cne.21276

PubMed Abstract | CrossRef Full Text | Google Scholar

Tapia, J. C., Kasthuri, N., Hayworth, K. J., Schalek, R., Lichtman, J. W., Smith, S. J., et al. (2012). High-Contrast En Bloc Staining of Neuronal Tissue for Field Emission Scanning Electron Microscopy. Nat. Protoc. 7 (2), 193–206. doi:10.1038/nprot.2011.439

PubMed Abstract | CrossRef Full Text | Google Scholar

Ueda, H. R., Dodt, H.-U., Osten, P., Economo, M. N., Chandrashekar, J., and Keller, P. J. (2020). Whole-Brain Profiling of Cells and Circuits in Mammals by Tissue Clearing and Light-Sheet Microscopy. Neuron 106 (3), 369–387. doi:10.1016/j.neuron.2020.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Rossum, G., and Drake, F. (2009). 10. CA: Scotts Valley, 1593511. CreateSpace. Python 3 Reference Manual.

Google Scholar

Wachowiak, M., Economo, M. N., Diaz-Quesada, M., Brunert, D., Wesson, D. W., White, J. A., et al. (2013). Optical Dissection of Odor Information Processing In Vivo Using GCaMPs Expressed in Specified Cell Types of the Olfactory Bulb. J. Neurosci. 33 (12), 5285–5300. doi:10.1523/JNEUROSCI.4824-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, S. M., Schwyn, D. A., Mokso, R., Wicklein, M., Müller, T., Doube, M., et al. (2014). In Vivo Time-Resolved Microtomography Reveals the Mechanics of the Blowfly Flight Motor. PLoS Biol. 12 (3), e1001823. doi:10.1371/journal.pbio.1001823

PubMed Abstract | CrossRef Full Text | Google Scholar

Walsh, C. L., Tafforeau, P., Wagner, W. L., Jafree, D. J., Bellier, A., Werlein, C., et al. (2021). Imaging Intact Human Organs with Local Resolution of Cellular Structures Using Hierarchical Phase-Contrast Tomography. Nat. Methods 18 (12), 1532–1541. doi:10.1038/s41592-021-01317-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Walter, A., Paul-Gilloteaux, P., Plochberger, B., Sefc, L., Verkade, P., Mannheim, J. G., et al. (2020). Correlated Multimodal Imaging in Life Sciences: Expanding the Biomedical Horizon. Front. Phys. 8. doi:10.3389/fphy.2020.00047

CrossRef Full Text | Google Scholar

Wanner, A. A., Genoud, C., Masudi, T., Siksou, L., and Friedrich, R. W. (2016). Dense EM-Based Reconstruction of the Interglomerular Projectome in the Zebrafish Olfactory Bulb. Nat. Neurosci. 19 (6), 816–825. doi:10.1038/nn.4290

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilke, S. A., Antonios, J. K., Bushong, E. A., Badkoobehi, A., Malek, E., Hwang, M., et al. (2013). Deconstructing Complexity: Serial Block-Face Electron Microscopic Analysis of the Hippocampal Mossy Fiber Synapse. J. Neurosci. 33 (2), 507–522. doi:10.1523/JNEUROSCI.1600-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Winnubst, J., Bas, E., Ferreira, T. A., Wu, Z., Economo, M. N., Edson, P., et al. (2019). Reconstruction of 1,000 Projection Neurons Reveals New Cell Types and Organization of Long-Range Connectivity in the Mouse Brain. Cell 179 (1), 268–281 e213. doi:10.1016/j.cell.2019.07.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Witte, K., Späth, A., Finizio, S., Donnelly, C., Watts, B., Sarafimov, B., et al. (2020). From 2D STXM to 3D Imaging: Soft X-Ray Laminography of Thin Specimens. Nano LettNano, doi:10.1021/acs.nanolett.9b04782

CrossRef Full Text | Google Scholar

Yildirim, M., Sugihara, H., So, P. T. C., and Sur, M. (2019). Functional Imaging of Visual Cortical Layers and Subplate in Awake Mice with Optimized Three-Photon Microscopy. Nat. Commun. 10 (1), 177. doi:10.1038/s41467-018-08179-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: staining, warping, olfactory bulb, 2-photon calcium imaging, synchrotron X-ray, volume EM, correlative multimodal imaging

Citation: Zhang Y, Ackels T, Pacureanu A, Zdora M-C, Bonnin A, Schaefer AT and Bosch C (2022) Sample Preparation and Warping Accuracy for Correlative Multimodal Imaging in the Mouse Olfactory Bulb Using 2-Photon, Synchrotron X-Ray and Volume Electron Microscopy. Front. Cell Dev. Biol. 10:880696. doi: 10.3389/fcell.2022.880696

Received: 21 February 2022; Accepted: 22 April 2022;
Published: 08 June 2022.

Edited by:

Christopher Guerin, Vlaams Instituut voor Biotechnologie, Belgium

Reviewed by:

Yunfeng Hua, Shanghai Jiao Tong University, China
Kara Fulton, Harvard Medical School, United States

Copyright © 2022 Zhang, Ackels, Pacureanu, Zdora, Bonnin, Schaefer and Bosch. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Andreas T. Schaefer,; Carles Bosch,

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.