In situ X-ray-assisted electron microscopy staining for large biological samples

  1. Sebastian Ströh
  2. Eric W Hammerschmith
  3. David W Tank
  4. H Sebastian Seung
  5. Adrian Andreas Wanner  Is a corresponding author
  1. Princeton Neuroscience Institute, Princeton University, United States
  2. Paul Scherrer Institute, Switzerland

Abstract

Electron microscopy of biological tissue has recently seen an unprecedented increase in imaging throughput moving the ultrastructural analysis of large tissue blocks such as whole brains into the realm of the feasible. However, homogeneous, high-quality electron microscopy staining of large biological samples is still a major challenge. To date, assessing the staining quality in electron microscopy requires running a sample through the entire staining protocol end-to-end, which can take weeks or even months for large samples, rendering protocol optimization for such samples to be inefficient. Here, we present an in situ time-lapsed X-ray-assisted staining procedure that opens the ‘black box’ of electron microscopy staining and allows observation of individual staining steps in real time. Using this novel method, we measured the accumulation of heavy metals in large tissue samples immersed in different staining solutions. We show that the measured accumulation of osmium in fixed tissue obeys empirically a quadratic dependence between the incubation time and sample size. We found that potassium ferrocyanide, a classic reducing agent for osmium tetroxide, clears the tissue after osmium staining and that the tissue expands in osmium tetroxide solution, but shrinks in potassium ferrocyanide reduced osmium solution. X-ray-assisted staining gave access to the in situ staining kinetics and allowed us to develop a diffusion-reaction-advection model that accurately simulates the measured accumulation of osmium in tissue. These are first steps towards in silico staining experiments and simulation-guided optimization of staining protocols for large samples. Hence, X-ray-assisted staining will be a useful tool for the development of reliable staining procedures for large samples such as entire brains of mice, monkeys, or humans.

Editor's evaluation

This important study explores the kinetics of heavy metal staining of tissue using time-lapse imaging with X-ray micro computed tomography (CT). Introducing a compelling approach to investigate staining in situ, this work will be of interest to the wide community of scientists preparing biological samples in particular for large-volume electron microscopy. It will become a reference for the field in establishing a quantitative tool for assessing and developing staining protocols.

https://doi.org/10.7554/eLife.72147.sa0

Introduction

In the past decade the image acquisition rates of biological electron microscopy facilities have been scaled to 107–109 pixels per second through parallelization and automation (Hayworth et al., 2014; Schalek et al., 2011; Ren and Kruit, 2016; Denk and Horstmann, 2004; Hayworth et al., 2020; Eberle et al., 2015; Graham et al., 2019; Xu et al., 2017). The increase in imaging throughput for biological samples has been mainly driven by the emerging Neuroscience field of connectomics which aims to densely reconstruct neuronal circuits with synaptic resolution (Briggman et al., 2011; Helmstaedter et al., 2013; Bock et al., 2011; Zheng et al., 2018; Scheffer et al., 2020; Kornfeld et al., 2017; Wanner and Friedrich, 2020; Lee et al., 2016; Kasthuri et al., 2015; Wilson et al., 2019; Schmidt et al., 2017; Svara et al., 2018; Vishwanathan et al., 2017). Also, the analysis of the terabyte-sized electron microscopy datasets produced by these studies is becoming increasingly automated using machine-learning and computer vision (Januszewski et al., 2018; Dorkenwald et al., 2020; Schubert et al., 2019; Dorkenwald et al., 2017; Berning et al., 2015; Staffler et al., 2017; Jain et al., 2010; Buhmann et al., 2020; Vergara et al., 2020; Turner et al., 2020). Both the image acquisition as well as the image analysis for these types of datasets are being scaled by parallelization to large samples on the order of several cubic millimeters or even entire brains. However, a remaining obstacle is the lack of reliable tissue processing and staining protocols for large samples where the smallest dimension is greater than 1 mm.

Existing en bloc electron microscopy staining protocols have been optimized for staining small samples with dimensions of less than 1 mm (Genoud et al., 2018; Hua et al., 2015; Deerinck et al., 2010; Tapia et al., 2012). Using aldehyde-stabilized cryopreservation (McIntyre and Fahy, 2015) the cellular ultrastructure can be preserved even in large tissue blocks with dimensions exceeding 1 mm and entire brains. Despite pioneering work on en bloc staining protocols for whole mouse brains (Mikula and Denk, 2015; Mikula et al., 2012), large en bloc stained samples still suffer from artifacts such as inhomogeneous staining and membrane or tissue cracks.

En bloc sample preparation for electron microscopy generally requires tissue fixation, staining and embedding in resin. First, the macromolecules in the tissue are stabilized via crosslinking by diffusing or perfusing buffered solutions of fixatives such as formaldehyde (Claude and Fullam, 1945; Fox et al., 1985) and glutaraldehyde (Sabatini et al., 1964). Because biological tissue is composed mostly of carbon and other low atomic number elements, the tissue is stained with heavy metals to increase the contrast in electron micrographs (Bahr, 1954; Porter et al., 1945). These heavy metal stains are composed of electron dense atoms such as osmium, lead or uranium. In addition, some of these heavy metals (e.g. osmium) also act as fixatives (Bahr, 1955). Finally, the stained tissue gets dehydrated and embedded in resin. The most commonly used resins are epoxy resins (Maaløe and Andersen, 1956), such as araldite (Glauert et al., 1956) or Epon (Finck, 1960) because of their thermal stability and electron transparency. To inspect the staining quality in an electron microscope, ultrathin sections (<100 nm) are typically collected from the embedded tissue using an ultramicrotome.

Many steps in the classic en bloc electron microscopy staining protocols are based on passive diffusion of chemicals into the tissue. Passive diffusion is one of the main bottlenecks for preparing large samples (smallest dimension >1 mm) for electron microscopy (Burkl and Schiechl, 1968; Medawar, 1941). For small samples (smallest dimension <1 mm) a typical en bloc staining protocol takes about 10–15 days including sectioning (Tapia et al., 2012). However, for large samples or whole brains, a diffusion-based staining protocol takes several weeks or even months (Mikula et al., 2012; Mikula and Denk, 2015; Masís et al., 2018).

To date, electron microscopy staining is a ‘black box’ and changes to staining protocols can only be assessed using an electron microscope, which in turn requires a sample to be run through the entire staining protocol. Conventional approaches for optimizing the parameters of staining protocols rely on sequential screening of hundreds of samples that have been processed end-to-end. However, sequential screening is very inefficient for months-long protocols.

Recently, X-ray based computed micro-tomography (μCT) has been introduced for a relatively fast, macroscopic assessment of staining quality and tissue integrity of resin-embedded whole mouse brains (Mikula and Denk, 2015; Kuan et al., 2020; Dyer et al., 2017). (Mikula and Denk, 2015) showed that the pixel intensities of serial-section EM images and the intensity of the corresponding reslice are similar. Building on this pioneering work, we developed in situ time-lapsed X-ray-assisted staining to observe the staining process while the samples are in the staining solution (Figure 1a). We used X-ray-assisted staining to explore the micro-scale tissue mechanics and the kinetics of the heavy metal diffusion and accumulation in large aldehyde-fixed brain tissue blocks, resulting in new insights on how different staining agents affect the tissue. X-ray-assisted staining opens the ‘black box’ of electron microscopy staining protocols. Each staining step can be monitored and assessed in real time. This enables in silico optimization of electron microscopy staining protocols, which will be particularly useful for the development of staining procedures for large biological samples such as whole brains (Figure 1g).

Figure 1 with 3 supplements see all
X-ray assisted staining.

(a) Experimental setup. The sample was glued into a thin layer of Sylgard and sits in staining solution in a sealed glass vial on a motorized stage. X-rays are emitted from an X-ray source, pass through the sample, and the spatial distribution of the transmitted X-rays is detected and used to form a projection image of the X-ray absorption in the sample. Over time, the heavy metals diffuse into the tissue and the accumulation of heavy metals change the X-ray absorption properties of the sample. For the quantitative results presented in this study we used cylindrical tissue samples extracted with a 4 mm diameter biopsy punch from the cortex of mice. (b) Projection images of a 4 mm punch of mouse cortex incubated in 2% buffered osmium tetroxide (OsO4) solution after 10 min, 1 hr, and 10 hr. The osmium diffuses passively from the cortical surface and the brain ventricles into the tissue and forms a staining front of tissue-bound osmium (large arrowhead) that moves towards the center of the sample as time progresses. The accumulation of heavy metals results in an increase in the pixel intensity (a.u.) of the X-ray projection image, which was measured along the radial axis of cortical depth from the cortical surface towards the white matter. The small arrowhead indicates the upper edge of the 2 mm thin layer of Sylgard. (c) Measured pixel intensities of the X-ray projection images for different concentrations of buffered OsO4 solutions. The higher/darker the intensity, the more X-rays are absorbed. The pixel intensity/X-ray absorption scales linearly with the concentration/density of osmium (black line, Pearson correlation coefficient r=0.99, <10–6; n=3 for each concentration). (d) Propagation of the staining front in 4 mm brain punches immersed in 2% buffered OsO4 solution (n=13) measured along the radial axis of cortical depth (mean ±s.d.). The propagation can be fitted by a quadratic model (dashed black lines) in which the penetration depth x of the staining front is proportional to the square root of the incubation time t (residual standard error: SEres = 14.68 μm). (e) Intensity profiles of the spatial heavy metal accumulation after t=20 hr of incubation in 2% buffered OsO4 (blue, n=13), Ferro-redOs (green, n=6). Ferro-redOs with 5.56% (orange, n=4) and 11.25% formamide (red, n=4) solutions (mean ±s.d.). In Ferro-redOs the heavy metal staining is not homogenous and there is a more densely stained tissue band at a depth of 300–800 μm. (f) Tissue expansion quantified as change in sample height of aldehyde-fixed brain punches immersed in OsO4 (blue, n=13), Ferro-redOs (green, n=6), Ferro-redOs with 5.56% (orange, n=4) and 11.25% formamide (red, n=4), respectively, as a function of the immersion duration (mean ±s.d.). All solutions were buffered by cacodylate. A monomolecular growth model was fitted to the average tissue expansion in OsO4 solution (black dashed line, residual standard error SEres = 0.06%). (g) In situ X-ray-assisted staining of large tissue samples. The images show projections of a whole mouse brain incubated in 2% buffered OsO4 at different points in time across 2.5 days. The arrowheads indicate the location of the staining front of accumulated osmium. Scale bar 2 mm.

Figure 1—source data 1

Experimentally measured accumulation of heavy metals (blue traces) in buffered 2% OsO4 solution in n=13 different samples.

The corresponding simulated model was fitted individually to the experimental data of each sample (black traces). The first row of every sample shows the temporal profile of osmium accumulation at different cortical depths, while the second row shows the spatial profile of osmium accumulation at different times after staining onset. The black circles indicate the location of axon tracts that tended to accumulate osmium faster and stronger. The third row shows the measured spatial-temporal accumulation of osmium (left), the corresponding fitted model simulations (middle) and the distribution of osmium in projection images of the sample at the beginning of the experiment and after about 22 hr of incubation (right).

https://cdn.elifesciences.org/articles/72147/elife-72147-fig1-data1-v2.pdf
Figure 1—source data 2

Experimentally measured accumulation of heavy metals in a buffered solution of 2% OsO4 reduced with 2.5% potassium ferrocyanide K4[Fe(CN)6] (Ferro-redOs) in n=6 different samples.

The first row shows the temporal profile of osmium accumulation at different cortical depths, while the second row shows the spatial profile of osmium accumulation at different times after staining onset. The black circles indicate the location of axon tracts that tended to accumulate osmium faster and stronger. The third row shows the measured spatial-temporal accumulation of osmium (left) and the corresponding distribution of osmium in projection images of the sample at the beginning of the experiment and after about 22 hoursr of incubation (right).

https://cdn.elifesciences.org/articles/72147/elife-72147-fig1-data2-v2.pdf
Figure 1—source data 3

Experimentally measured accumulation of heavy metals in buffered 2% Ferro-redOs solution with 5.56% formamide in n=4 different samples.

The first row of every sample shows the temporal profile of osmium accumulation at different cortical depths, while the second row shows the spatial profile of osmium accumulation at different times after staining onset. The black circles indicate the location of axon tracts that tended to accumulate osmium faster and stronger. The third row shows the measured spatial-temporal accumulation of osmium (left) and the corresponding distribution of osmium in projection images of the sample at the beginning of the experiment and after about 22 hoursr of incubation (right).

https://cdn.elifesciences.org/articles/72147/elife-72147-fig1-data3-v2.pdf
Figure 1—source data 4

Experimentally measured accumulation of heavy metals in buffered 2% Ferro-redOs solution with 11.25% formamide in n=4 different samples.

The first row of every sample shows the temporal profile of osmium accumulation at different cortical depths, while the second row shows the spatial profile of osmium accumulation at different times after staining onset. The black circles indicate the location of axon tracts that tended to accumulate osmium faster and stronger. The third row shows the measured spatial-temporal accumulation of osmium (left) and the corresponding distribution of osmium in projection images of the sample at the beginning of the experiment and after about 22 hoursr of incubation (right).

https://cdn.elifesciences.org/articles/72147/elife-72147-fig1-data4-v2.pdf
Figure 1—source data 5

Experimentally measured accumulation of heavy metals in a buffered solution of 2% OsO4 reduced with 2.5% potassium ferricyanide K3[Fe(CN)6] (Ferri-redOs) in n=2 different samples.

The first row shows the temporal profile of osmium accumulation at different cortical depths, while the second row shows the spatial profile of osmium accumulation at different times after staining onset. The black circles indicate the location of axon tracts that tended to accumulate osmium faster and stronger. The third row shows the measured spatial-temporal accumulation of osmium (left) and the corresponding distribution of osmium in projection images of the sample at the beginning of the experiment and after about 22 hoursr of incubation (right).

https://cdn.elifesciences.org/articles/72147/elife-72147-fig1-data5-v2.pdf

Results

We developed in situ time-lapsed X-ray-assisted staining with the goal of facilitating and accelerating the optimization of staining protocols for large biological samples such as whole mouse brains (Abbott et al., 2020). 4mm-brain punches of transcardially perfused mice were immersed in aldehyde fixatives for 36 hr. After washing the sample blocks with cacodylate buffer, the samples were immersed in staining solution and immediately placed in the acquisition chamber of a Zeiss Xradia 520 Versa 3D for X-ray microscopy (Figure 1a). The Xradia can be operated in two different modes:

  • In the μCT mode, projections of the sample are acquired at different rotation angles in order to reconstruct a 3D computed tomograph of the sample. Depending on the required signal to noise and resolution, this mode allowed us to acquire approximately 1–2 computed tomographs per hour.

  • In single-projection mode, a projection of the sample is acquired in a fixed position every few seconds without rotating the sample.

The advantage of the μCT mode is that arbitrary virtual reslices can be extracted to assess the staining progression in various parts of the sample. However, with our X-ray microscope, the temporal resolution in μCT mode was limited to 30–60 min per tomograph. Therefore we performed all quantitative experiments in single-projection mode where the temporal resolution for monitoring of the heavy metal diffusion and accumulation was on the order of a few seconds. Note, in this mode, the pixel intensities of the resulting X-ray projection images are accumulated through the entire sample and therefore the correspondence to the pixel intensities in thin serial-section EM images is less accurate than in μCT mode (Figure 1—figure supplement 1). As soon as the samples got placed in the staining solutions, the heavy metals started accumulating in the immersed tissue. The stained tissue absorbs more X-rays than the unstained tissue, resulting in a noticeable intensity difference that can be used to track the diffusion and accumulation of the heavy metals in the tissue (Figure 1b).

Quadratic scaling of incubation times with sample size

Osmium tetroxide (OsO4) is one of the most commonly used contrast agents for lipid staining in electron microscopy due to its large atomic number and its ability to integrate into cellular membranes (Porter et al., 1945; Palade, 1952; Watson, 1958). The immersion times for OsO4 vary between a few minutes to several days, depending on the sample size and tissue type, and are usually determined empirically (Genoud et al., 2018; Hua et al., 2015; Tapia et al., 2012; Mikula and Denk, 2015). However, the kinetics of OsO4 staining of biological tissue are not well understood and so far had to be determined experimentally by trial and error. We therefore set out to measure the diffusion and accumulation of OsO4 by placing 4 mm punches of aldehyde-fixed mouse brains in 2% buffered OsO4 solution. The osmium diffused into the tissue and accumulated in the sample during the staining process (Figure 1b). The intensity of X-ray absorption in the projection images scales linearly with the local concentration/density of osmium (Figure 1c). The density of OsO4 in the tissue increased beyond the density in the surrounding staining solution, indicating that the density of binding sites for OsO4 in the tissue is higher than the concentration of the OsO4 in the staining solution (Figure 1b). The diffusion and spatio-temporal accumulation of OsO4 results in a staining front that propagates towards the center of the sample (Figure 1b). As expected for diffusive processes (Carnevale et al., 1979), the propagation of the OsO4 staining front can be approximated by a quadratic model (Figure 1d). As in the case of other fixatives (Medawar, 1941), the osmium staining penetration in aldehyde fixed tissue obeys a quadratic scaling law, or rather a “rule of thumb”, for how the necessary incubation time t depends on the staining depth or sample size x:

(1) t x2

This means, for example, that if one would adapt an established OsO4 staining protocol for 3 X larger samples, one would have to prolong the incubation time by 9 X in order to produce comparable staining results. Similarly, the time it takes for any point in the sample to reach a given concentration is proportional to the square of its distance to the sample surface (Carnevale et al., 1979).

Monitoring staining kinetics and tissue deformation

After 20 hr of incubation the staining density of osmium is homogeneous across the first 1000 μm of cortical depth (Figure 1e). In addition, the ultrastructure is well preserved (Figure 1—figure supplement 2). Reduced osmium is another commonly used staining agent that is known to result in higher contrast for electron microscopy images than non-reduced OsO4. Typically, a buffered 2% OsO4 solution is reduced with 2.5% potassium ferrocyanide (K4[Fe(CN)6]) (Hua et al., 2015; Willingham and Rutherford, 1984; Mikula and Denk, 2015), which we will call Ferro-redOs throughout this manuscript. Consistent with previous reports on Ferro-redOs staining in large samples (Hua et al., 2015; Mikula and Denk, 2015) we found that for Ferro-redOs the accumulation of heavy metals peaks at a depth between 300 and 800 μm (Figure 1e), whereas the tissue above or below that depth is stained less. Traditionally, this band of more heavily stained tissue has been associated with precipitated osmium that hinders the diffusion and prevents homogenous staining (Genoud et al., 2018; Hua et al., 2015; Mikula and Denk, 2015), in particular deeper in the tissue. No such band is present in tissue immersed in non-reduced OsO4. (Mikula and Denk, 2015) reported that homogeneous staining of large samples with Ferro-redOs was achieved by adding formamide to the Ferro-redOs solution. However, the underlying mechanisms by which formamide acts are not known (Mikula and Denk, 2015; Genoud et al., 2018). (Mikula and Denk, 2015) hypothesized that formamide might prevent precipitation by generally solubilizing compounds or that it might allow highly charged molecules to cross membranes more easily. Adding formamide to the Ferro-redOs solution indeed resulted in more homogeneous staining (Figure 1e, ), but the heavy metal density was lower than in osmium only. (Mikula and Denk, 2015) reported that the tissue tends to expand at high concentrations of formamide (>50%). But we found that the tissue expands even at lower formamide concentrations, while the amount of expansion depends on both the formamide concentration and the incubation duration (Figure 1f). For a concentration of 11.25% formamide, the sample height increased by about 15% within 20 hr of incubation, whereas for a formamide concentration of 5.56%, the tissue height expanded only by approximately 10% (Figure 1f). Note, in the first 3 hr of incubation in Ferro-redOs with 5.56% formamide the tissue actually shrank, suggesting that there are opposite forces acting on the tissue. Indeed, in Ferro-redOs solution without any formamide, the tissue shrank by about 5% in height (Figure 1f). In contrast, we found that in 2% buffered OsO4 solution the sample height expands by about 5% (Figure 1f). However, we found no sample expansion if the reduced osmium solution was prepared with 2.5% potassium ferricyanide K3[Fe(CN)6]+ 8 a,b), another commonly used reducing agent for OsO4.

Hua et al., 2015 suggested an alternative approach to achieve homogeneous staining in 1 mm brain punches with Ferro-redOs without formamide, in which osmium and the reducing agent are applied separately. First, the samples are immersed in buffered 2% OsO4 solution for 90 min. Subsequently, the samples are placed for 90 min in buffered 2.5% K4[Fe(CN)6] without any washing step in-between. We repeated this procedure for 4 mm brain punches: First, the samples were incubated in 2% buffered OsO4 for 22 hours resulting in homogeneous staining throughout the sample. Next, we placed the sample directly without any washing step in 2.5% buffered K4[Fe(CN)6] solution. (Hua et al., 2015) hypothesized that the main effect of reducing agents such as K4[Fe(CN)6] is to convert VIII-oxidized water-soluble osmium into an VI-oxidized water-soluble form which is thought to generate additional, non-polar OsO2 to be deposited in the membrane increasing the heavy metal content and the contrast in electron microscopy. However, we found that the heavy metal density in the K4[Fe(CN)6] immersed samples did not increase, but rather decrease from the sample surface towards the center as time progressed (Figure 2a). This suggests that K4[Fe(CN)6] removes or “washes out” osmium from the sample (Litman and Barrnett, 1972) resulting in an inverse staining gradient. This “washing effect” is stronger for longer incubation times as well as close to the sample surface and around blood vessels (anecdotal observation in electron microscopy images, data not shown). For incubation times longer than 12 hr, K4[Fe(CN)6] readily dissolves and disintegrates the tissue (Figure 2—figure supplement 1). Interestingly, no ‘washing effect’ was observed if the same procedure was repeated with 2.5% potassium ferricyanide K3[Fe(CN)6], but the sample height shrank by about 2% (Figure 1—figure supplement 3,f, ). Similarly, only a small reduction in heavy metal density could be observed when the sample was washed with double-distilled H2O for 20 hr (Figure 2—figure supplements 2 and 3).

Figure 2 with 5 supplements see all
Kinetics and tissue mechanics of electron microscopy staining.

(a) Spatio-temporal washout of osmium in samples incubated in 2.5% potassium ferrocyanide K4[Fe(CN)6] (n=6). The samples have been stained with 2% OsO4 for 22 hrs prior to placing them in K4[Fe(CN)6]. The osmium gets washed out from the sample surface towards the center as the incubation time increases. (b) Average experimentally measured spatio-temporal osmium density (nmol/mm3) accumulation in samples (n=13) immersed in 2% buffered OsO4 solution. Because of the surface curvature of the cortical samples, the tissue thickness decreases towards the surface (x=0) and results in a lower intensity in the projection images (see tissue geometry model in f). (c) Simulated spatio-temporal osmium density (nmol/mm3) accumulation fitted to the average experimental data in b. (d) Density profiles (nmol/mm3) of the spatial accumulation of osmium (n=13, mean ±s.d.) after 1 hr (orange), 2 hrs (red), and 10 hr (blue) of incubation in 2% buffered OsO4 overlaid with the simulation of the diffusion-reaction-advection model (dashed lines). (e) Density profiles (nmol/mm3) of the temporal accumulation of osmium (n=13, mean ±s.d.) at a depth of 200 μm (blue), 600 μm (red), and 1000 μm (orange) during incubation in 2% buffered OsO4, overlaid with the simulation of the diffusion-reaction-advection model (dashed lines). (f) The geometry of the 4 mm brain punches was modeled as a cylinder with a curved surface. The curvature of the sample surface is approximated by a circle with radius 12HH2+14R2 , where R is diameter of the cylinder and H is the height of the curved surface. The projected tissue thickness as a function of the depth x from the sample surface is then given by d(x): d(x)=2xHH2+14R2-x2 , for 0xH and d(x)=R, for x>H (g) Diffusion-reaction-advection model for osmium staining. The model combines four coupled processes: (1) Free OsO4 diffuses passively into the tissue and (2) binds to the available binding sites. (3) In the presence of freely diffusing OsO4 previously masked binding sites are slowly turned into additional binding sites for OsO4. (4) The sample expands by about 5% in sample height within 20 hr of incubation (see Figure 1f).

Figure 2—source data 1

Experimentally measured washout of heavy metals in buffered 2.5% potassium ferrocyanide K4[Fe(CN)6] solution in n=6 different samples that have previously been stained with buffered OsO4 for 22 hr.

The first row of every sample shows the temporal profile of osmium washout at different cortical depths, while the second row shows the spatial profile of osmium washout at different times after staining onset. The black circles indicate the location of axon tracts. The third row shows the measured spatial-temporal washout of osmium (left) and the corresponding distribution of osmium in projection images of the sample at the beginning of the experiment and after about 22 hr of incubation (right).

https://cdn.elifesciences.org/articles/72147/elife-72147-fig2-data1-v2.pdf
Figure 2—source data 2

Experimentally measured washout of heavy metals in buffered 2.5% potassium ferricyanide K3[Fe(CN)6] solution in n=4 different samples that have previously been stained with buffered OsO4 for 22 hrs.

The first row of every sample shows the temporal profile of osmium washout at different cortical depths, while the second row shows the spatial profile of osmium washout at different times after staining onset. The black circles indicate the location of axon tracts. The third row shows the measured spatial-temporal washout of osmium (left) and the corresponding distribution of osmium in projection images of the sample at the beginning of the experiment and after about 22 hoursr of incubation (right).

https://cdn.elifesciences.org/articles/72147/elife-72147-fig2-data2-v2.pdf

The kinetics of osmium tetroxide staining

The chemistry and diffusion-reaction-advection kinetics of OsO4 staining in aldehyde-fixed tissue are not well understood. In brain tissue, OsO4 is thought to mainly react with carbon-carbon double bonds, in particular of unsaturated lipids, but it has also been reported to react with some amino acids and nucleotides (Khan et al., 1961; Korn, 1967; White et al., 1976; Schroeder, 1980). At the beginning of the OsO4 incubation, the density of OsO4 steeply increases and peaks close to the surface of the sample. As the incubation time progresses, the OsO4 continues to diffuse and accumulates deeper in the tissue (Figure 2b, d and e). Due to the curved surface of the cortical tissue sample, the projected tissue thickness is lower (Figure 2f) and there are fewer binding sites for osmium, which results in a lower density of osmium closer to the sample surface (Figure 2b). At any given time, the accumulation seems to plateau from the surface towards the center of the sample, indicating that the binding of OsO4 saturates (Figure 2d and e). However, we found that the saturated density in the samples slowly increased over time from about 129 nmol/mm3 after 1 hr of incubation to 160 nmol/mm3 after 10 hr of incubation (Figure 2b, d and e), while the density of the solution remained stable at about 82 nmol/mm3. This suggests that additional binding sites for OsO4 become available slowly over time. The reactions underlying the ‘unmasking’ of binding sites are not known, but we hypothesize that it depends at least partially on the availability and the local concentration of unbound, freely diffusing osmium. Hence, at the macroscale level, the accumulation of osmium in biological tissue can be described by four coupled processes (Figure 2g): As the osmium diffuses passively into the tissue, it binds to the available binding sites in the tissue and additional binding sites are slowly being unmasked, while the tissue expands by around 5% within 20 hr of incubation (Figure 1f). The kinetics of these four processes can be described with the following coupled (non-linear) partial differential equations:

(2) f(x,t)t=De2f(x,t)x2kons(x,t)f(x,t)x[f(x,t)v(x,t)]
(3) b(x,t)t=kons(x,t)f(x,t)x[b(x,t)v(x,t)]
(4) s(x,t)t=kons(x,t)f(x,t)+kunmaskm(x,t)f(x,t)x[s(x,t)v(x,t)]
(5) m(x,t)t=kunmaskm(x,t)f(x,t)x[m(x,t)v(x,t)]

f(x,t) is the density of freely diffusing OsO4, b(x,t) is the density of osmium bound to the tissue, s(x,t) is the density of available binding sites and m(x,t) corresponds to the density of masked binding sites at time t and a radial depth x from the cortical surface. De is the effective diffusion coefficient for OsO4, kon is the reaction rate constant for the binding of osmium to the tissue and kunmask is the reaction rate constant at which masked binding sites get unmasked. In addition, the local change of the tissue geometry, e.g. the expansion as observed for OsO4, results in a velocity field v(x,t) that describes the local spatio-temporal flow of the density of each component c(x,t):

(6) xc(x,t)v(x,t)=c(x,t)xv(x,t)+c(x,t)v(x,t)x

The flow consists of an advection term c(x,t)xv(x,t), corresponding to elemental volumes moving with the local flow, and a dilution term c(x,t)v(x,t)x that describes the local volume change. While we assume that the tissue density is homogeneous across the cortical thickness, we consider differences in the sample thickness due to the curvature of the cortical surface. Additional information on the model assumptions can be found in the methods section.

We fitted the parameters of the system of reaction-diffusion-advection equations for a tissue sample with the projection geometry shown in Figure 2f to the average measured spatio-temporal accumulation of OsO4 (Figure 2b–e) as well as to the individual samples (n=13, ). The fitted parameters of the diffusion-reaction-advection model (Figure 2—figure supplement 4) capture the measured diffusion-reaction-advection kinetics of osmium staining well (Figure 2c; residual standard error SEres = 0.026 ± 0.005, mean ±s.d.).

Discussion

In situ X-ray-assisted staining is a new tool for monitoring the different steps of en bloc staining protocols for electron microscopy. By opening the ‘black box’ of staining protocols, X-ray-assisted staining enables, for the first time, to observe directly the diffusion of heavy metals and their microscopic effects on the tissue. While the technique will be useful to develop novel staining protocols for small to intermediate sized samples with dimensions less than 1 mm, it will particularly accelerate and facilitate the development and optimization of new staining protocols for large biological samples such as whole brains (Figure 1g).

We monitored the accumulation of osmium in different commonly used staining solutions and found significant differences in the accumulation densities. In addition, we found that the tissue expands and shrinks significantly depending on the composition of the staining solution. In particular, we found that the tissue height increased by 5% in buffered osmium solution, while the tissue height decreased by 5% in buffered Ferro-redOs solution. We hypothesize that this shrinkage could provide a mechanistic explanation for the long standing issue of diffusion barrier formation in large samples during Ferro-redOs incubation.

The method also revealed that K4[Fe(CN)6], a reducing agent that was believed to enhance the accumulation of osmium in membranes, in contrast can wash out osmium from stained tissue if applied separately. If K4[Fe(CN)6] removes osmium from the sample, why does the protocol modification of Hua et al., 2015 result in enhanced electron microscopy contrast for samples with dimensions of up to 1 mm? A potential explanation could be that K4[Fe(CN)6] clears the cytosol more efficiently than the cell membranes. As a consequence, the osmium content in the cytosol gets washed out more/faster, which might result in an overall electron microscopic contrast increase between the membranes and the cytosol. However, the chemical mechanisms by which K4[Fe(CN)6] interact with the tissue remain poorly understood and shining a light on those processes will require additional work.

X-ray assisted staining opens the door for in silico staining protocol optimization, by enabling to measure important heavy metal staining parameters such as the effective diffusion coefficients and reaction rate constants directly in situ. We derived a diffusion-reaction-advection model that describes the accumulation of osmium in homogeneous tissue by four distinct, but coupled, processes. Osmium diffuses passively into the tissue and binds there quickly to available binding sites, for example unsaturated lipids. Thereby the tissue expands by about 5% in height. Additional binding sites are being unmasked over the course of hours, which enables the tissue to accumulate additional osmium. While the binding of osmium appears to be nearly instantaneous, the passive diffusion and unmasking reactions are slow and appear to be rate-limiting for the staining. The mechanisms and reactions underlying the unmasking of additional binding sites for osmium are not known. Unsaturated lipids, one of the prominent binding sites for osmium, can be masked by lipid-protein complexes. In the presence of osmium, these complexes could dissociate and make the previously masked unsaturated lipid available for binding osmium (Clayton, 1959; Lauder and Beynon, 1993; Ashhurst, 1961; Wigglesworth, 1975). A recent study used time-of-flight secondary ion mass spectrometry to analyze the colocalization of different fatty acid species with osmium at the tissue surface (Belazi et al., 2009). As expected, the unsaturated fatty acids colocalize well with the tissue-bound osmium. However, they also report a remarkable complementary clustering of saturated and unsaturated fatty acids in osmium stained tissue that appears not to be present in unstained tissue. Hence, it could be that osmium can convert saturated lipids to unsaturated lipids through hydrogenation/oxidation. Alternatively, the masked binding sites could also be interpreted as a different binding partner that reacts much slower with osmium.

Biological samples such as brains are non-homogeneous media. They are composed of different tissue and cell types and therefore the diffusion of staining agents is expected to be heterogeneous. Nevertheless, the fitted model parameters reproduce the observed kinetics of osmium accumulation relatively well up to a depth of about 700–800 μm from the surface. Deeper in the tissue osmium that diffused from secondary sources such as ventricles or axon tracts appear to contribute to the accumulation of osmium. For example, we found that white matter and myelinated axon tracks tend to accumulate osmium faster and stronger compared to gray matter (, ). Such secondary sources and tissue inhomogeneities are not considered in our model. Still, the fitted model parameters appear consistent with comparable values from the literature, although we also find some discrepancies. For example, in mouse brains we expect to find about 117.975 nmol/mm3 double-bonds on unsaturated lipids (see methods for the derivation). Because the dominant binding reaction for osmium is assumed to be the formation of mono- and di-ester bonds with unsaturated lipids which has been described to result in the deposition of OsO2, the density of double bonds should roughly correspond to the density of available binding sites. But the proposed model predicts 2.999 times more binding sites, even without unmasking. One interpretation would be that unsaturated lipids are not the dominant reaction partner and that other reaction partners such particular amino acids or nucleotides account for about two third of the available binding sites. Alternatively, this could indicate that each unsaturated lipid double bond gets associated with on average 3 osmium atoms, for example through Os3O5- and Os3O5H- isotopes that have previously reported to exist in osmium-stained adipose mouse tissue (Belazi et al., 2009). Similarly, it has been reported that if the mono- and di-ester products are left standing in solution, a highly insoluble trimer osmium product is formed (Schroeder, 1980).

We expect that X-ray-assisted staining will be very useful in establishing a more refined picture of the heterogeneity of the heavy metal diffusion processes in different tissue types and compartments as well as in different conditions. For example, we also measured and modeled the accumulation of osmium in liver tissue (Figure 2—figure supplement 5). It will also be interesting to examine how changes in the extracellular space and different fixation protocols or microwave-assisted incubation affect the diffusion dynamics (Korogod et al., 2015; Fulton and Briggman, 2020; Jensen and Harris, 1989; Login and Dvorak, 1988).

Furthermore, X-ray-assisted staining can be used at higher resolution and in combination with computed tomography for producing a more detailed picture of the tissue mechanics involved in different staining steps. We and several other labs have noticed that brain tissue expands significantly during incubation in ‘osmium amplifiers’ such as the osmiophilic thiocarbohydrazide (TCH) (Seligman et al., 1966). (Mikula and Denk, 2015) reported that nitrogen bubbles are formed during TCH incubation that could be responsible for the formation of tissue cracks in large samples. X-ray-assisted staining could be used to detect and monitor the formation of tissue cracks and to explore ways to reduce these artifacts in large samples.

Aside from accelerating the development of novel staining protocols this new method will also allow for improved closed-loop quality control of the staining process. Large samples tend to be biologically more variable and more difficult to precisely fine-tune all the parameters in the staining protocol ahead of time to account for these individual differences. Hence, X-ray-assisted staining can be used to adjust the parameters in real time to compensate for these individual differences. This is particularly important for experiments in which the sample yield is low or the individual sample is precious, as for example in the case of human brains or functionally imaged brains of behaving animals that require several months of behavioral training. However, this will require calibrating what is measured with X-ray microscopy to the corresponding appearance at an electron microscopic level, because the image formation process and the resulting contrast, signal-to-noise, etc. for both imaging modalities depend not only on the staining protocol, but also on the imaging modality, energy, flux, etc. (Du and Jacobsen, 2018).

Methods

Sample preparation

Animal use procedures were approved by the Princeton University Institutional Animal Care and Use Committee (protocol number 2000) and carried out in accordance with National Institutes of Health standards (AAALAC International Institutional Number: Unit #1001, PHS assurance ID D16-00273). The experiments were carried out with C57BL/6 J mice (Jackson Laboratories, Bar Harbor, Maine, USA) of both sexes, age ranging from 13 to 58 weeks (mean ± sd: 27±9 weeks). Prior to tissue sample extraction, the mice were anesthetized by isoflurane inhalation (4%) and euthanized with an intraperitoneal injection of a ketamine (100 mg/kg) and xylazine (10 mg/kg) overdose. Next, the animals were perfused transcardially with ≥50 ml fixative solution (1.3–1.5% glutaraldehyde (GA) and 2.5–2.6% paraformaldehyde (PFA) in 0.14–0.15 M cacodylate buffer with 2.0–2.1 mM CaCl2 at pH 7.4; (GA, PFA) and cacodylate buffer Electron Microscopy Sciences, Hatfield, PA; CaCl2, Sigma Aldrich, St. Louis, MO). After perfusion, the animal was left for 15–60 min on ice. Subsequently, the brain was removed in order to extract tissue samples at the center of the cortical surface of each hemisphere using biopsy punches of 4 mm diameter. The liver tissue was extracted with a 4 mm biopsy punch from the center of the liver after perfusion. The samples were kept in fixative solution at 4 °C for 36 hr (4 mm biopsy punches) or for 12 hr (whole mouse brain). Next, the samples were washed 6–7 times for 30 min in 0.15 M cacodylate buffer with 2 mM CaCl2.

Glass vials were prepared prior to the experiments with a 1–2 mm layer of Sylgard 184 (Electron Microscopy Sciences, Hatfield, PA) at the bottom. The whole mouse brain was glued with a small drop of Vetbond tissue adhesive (3 M, Hanna Pharmaceutical Supply Co., Inc, Wilmington, DE) to the bottom of the vial. For the 4 mm biopsy punch samples, a hole with a diameter of 4 mm was cut out of the Sylgard layer. Next, the samples were placed into the sylgard hole and were glued to the bottom of the glass vial with a small drop of Vetbond tissue adhesive. Subsequently, the samples were immersed in fresh cacodylate buffer solution for the transfer to the X-ray microscope. At the microscope, presets (see X-ray microscopy section below) were loaded for image acquisition. After that the washing solution was replaced with the staining solution and the sample was placed into the X-ray microscope. The acquisition of X-ray images was typically started within 5.66±2.35 min. (mean ±s.d.) after adding the staining solution. The following staining solutions have been assessed at room temperature for this study (all buffered with 0.15 M cacodylate buffer and 2 mM CaCl2):

  • Osmium tetroxide: 1%, 2% and 3% OsO4 (prepared from 2 mL 4% osmium tetroxide solution, Electron Microscopy Sciences, Hatfield, PA)

  • Potassium ferrocyanide: 2.5% potassium ferrocyanide (prepared from crystalline 25 g potassium hexacyanoferrate(II) trihydrate powder, Sigma Aldrich, St. Louis, MO)

  • Ferro-redOs: 2% OsO4, 2.5% potassium ferrocyanide

  • Ferro-redOs in 11.25% formamide: 2% OsO4, 2.5% potassium ferrocyanide, 11.25% formamide (prepared from 25 mL formamide solution, Sigma Aldrich, St. Louis, MO)

  • Ferro-redOs in 5.56% formamide: 2% OsO4, 2.5% potassium ferrocyanide, 5.56% formamide

  • Potassium ferricyanide: 2.5% potassium ferricyanide (EMS #13746-66-2)

  • Potassium ferricyanide reduced osmium: 2% OsO4, 2.5% potassium ferricyanide

Sample embedding, sectioning, and electron microscopy

After staining, the samples were washed with double-distilled H2O 3 times for 30 min. Next, the tissue was dehydrated with iteratively increasing concentrations of ethanol (EtOH, 200 proof, anhydrous,≥99.5%, Sigma Aldrich, St. Louis, MO) in double-distilled H2O for 30 min per step: 25% EtOH, 50% EtOH, 75% EtOH, 90% EtOH, 95% EtOH, 100% EtOH, 100% EtOH. Subsequently, the samples were transferred to 100% acetonitrile (ACS grade, Electron Microscopy Sciences, Hatfield, PA) for two 30-min washes. The tissue was embedded in hard Epon in an iterative concentration series: 1:2 Epon: acetonitrile for 12 hr, 1:1 Epon: acetonitrile for 24 hr, 2:1 Epon: acetonitrile for 2 days with fresh resin every 8–12 hr. Finally, the sample was left on the bench at room temperature for 6 days, with fresh resin every 8–12 hr, before it was placed in a 60 °C oven for curing for 5–7 days.

Serial sections of ~40 nm thickness were collected using an UC7 ultramicrotome from Leica. To assess the ultrastructural preservation of the tissue, the sections were imaged with a Zeiss Sigma SEM (Zeiss Microscopy, Oberkochen, Germany) with a backscattered electron detector (BSD) and the following acquisition parameters: 38 μs dwell time, 4 nm pixel size, probe current 200 nA, EHT voltage 6 kV, BSD Gain = High, all sectors of the detector ON.

X-ray microscopy

All X-ray microscopy experiments were performed on a Zeiss Xradia 520 Versa 3D X-ray microscope (Zeiss, Thornwood, NY). The samples were immersed in heavy metal solution in glass containers tightly sealed with Parafilm (Bemis Company, Inc, Neenah, WI) just before being placed in the recording chamber of the Xradia. Test images were acquired to check imaging parameters. For the biopsy punch samples, projections at a fixed view angle (rotation angle step size = 0°) were taken every 22–60 s resulting in about 1300–1350 images per sample over 22 hr of continuous acquisition. For the whole mouse brain, projections were taken every 22–60 s and after each projection the sample was rotated by 1°. All images were acquired using the image acquisition software Scout-and-Scan Control System (Zeiss, Thornwood, NY) with the following parameters:

  • Objective: 0.4 X

  • Tube Voltage = 100 kV

  • Output power = 7 W

  • Source filter = air (no filter)

  • Exposure time = 20 s.

  • Effective pixel size = 2–3 μm (biopsy punches) and 13–14 μm (whole mouse brain)

Image preprocessing

First, the individual images and metadata of each projection of each sample were extracted from the proprietary *.txrm and *.xrm files as saved by the Xradia image acquisition software. In addition, a time stamp of the acquisition time point was extracted from the image metadata for each projection. A custom cross-correlation based image alignment procedure (Wanner et al., 2016) was used to correct for translational offsets between images of subsequent time points. The projection images were flat-field corrected to account for radial intensity variations due to the cone-shaped X-ray beam as follows:

  1. For each sample, a median image immedian was calculated from the 2nd to 5th projection images right after acquisition onset.

  2. The following quadratic 2d function was fitted to log(immedian) using the built-in function lsqcurvefit in MATLAB:

    • Z(x,y)=b5x2+b4y2+b3x+b2y+b1

  3. At each time t, the projection image was flat-field corrected using Z(x,y) with the fitted parameters:

    • im^t(x,y)=imt(x,y)exp(Z(x,y))max(exp(Z))

Measuring the X-ray absorption and local heavy metal density

The pixel intensities of the detected absorbed X-rays were used as a proxy for the density of heavy metal accumulation in the tissue. In each sample, we quantified the pixel intensities in each projection taken at time t along a cross-section that pointed from the cortical surface radially towards the white matter (Figure 1b). The cross-sections were manually selected such that they started at a well detectable surface boundary and they were located at least 1 mm away from the vertical sample borders. For each sample i, the intensity wi(x,t) of the transmitted X-rays was quantified by averaging the pixel intensity of the cross-section at a depth x at time t within a width of ±25 pixels perpendicular to the axis of the cross-section. x=0 corresponds to the cortical surface of the brain samples or the surface of the liver for the liver samples. The intensity of X-ray absorption was defined as ui(x,t)=10,000-wi(x,t). For each time point t, the baseline intensity bi(t) was calculated by averaging the pixels outside of the sample within –50 μm and –150 μm from the sample surface along the cross-section within a width of ±25 pixels perpendicular to the axis of the cross-section. To account for baseline fluctuation in the staining solution and the X-ray illumination the intensity of X-ray absorption was baseline-corrected by uib(x,t)=ui(x,t)-bi(t)-median(bi) .

To standardize uib(x,t), the following quantities where calculated for each sample:

  1. bi_=mean(uib(x,t)) was calculated for all pixels along the cross-section outside of the sample within 0 μm and –200 μm from the sample surface.

  2. uimin=bi_10.90.9(mbi_)

  3. uiref=uimina(muimin)

  4. uinorm=(1+a)(muimin)

Finally, uib(x,t) was standardized to u^i(x,t)=uib(x,t)-uiref/uinorm.

This procedure ensures that the origin is at 0 for the linear relationship between the X-ray absorption intensities in the projection images and the different OsO4 concentrations as shown in Figure 1c and that a pixel intensity of 0.9 roughly corresponds to the saturated staining in 2% OsO4. For this, the factors a=0.61749 and m=6420 were determined empirically.

Measuring the linear relationship between X-ray absorption and OsO4 concentration

The X-ray absorption/pixel intensity of cacodylate-buffered solutions with 1%, 2% and 3% OsO4 was measured the same way as for the X-ray absorption in stained tissue.

Monitoring of tissue expansion and shrinkage

Intensity based heuristics were used in a semi-automated procedure to detect the cortical surface in each image of each sample. The glass vial bottom was used as a reference point to measure the distance to the cortical surface H(t) at time t. Right after placing the samples in the X-ray the height of the samples was H0=H(t=0) and the change in sample height (expansion or shrinkage) was quantified as H(t)-H0H0100%.

Tissue expansion modeling

The expansion of the tissue in buffered OsO4 solution was modeled by the monomolecular, saturating growth function g(t)=a1-a2exp(-tτ). The built-in function lsqcurvefit of MATLAB was used for parameter fitting. The fitted parameters for the average expansion in Figure 1f were a1=5.8697%, a2=5.6924%, τ=588.23 min. Under the assumption that the expansion of the tissue domain in buffered osmium solution is isotropic, the corresponding velocity field v(x,t) of the flow in Equation 6 can be characterized by:

v(x,t)=x(1+g(t)/100)t=xa2100τexp(tτ).

Tissue projection geometry model

While we assumed that the tissue density is homogeneous across the brain punch, we took into account the variation in the tissue thickness of the 4 mm brain punch. Figure 2f shows a sketch of the modeled tissue geometry. The geometry of the 4 mm brain punches is modeled as a cylinder with a curved surface. The diameter of the cylinder is R and the curvature of the sample surface is approximated by a circle with radius 12HH2+14R2 . The thickness of the sample as a function of the depth x from the sample surface is then given by d(x):

d(x)=2xHH2+14R2-x2 , for 0xH and d(x)=R, for x>H.

Diffusion-reaction-advection modeling

The diffusion-reaction-advection model described by Equations 2–5 makes the following assumptions:

  • the cortical surface is in contact with an infinite source of staining agents at constant concentration,

  • the composition of the 4 mm brain punches can be approximated by a semi-infinite, homogeneous medium,

  • the tissue thickness is given by d(x)=2xHH2+14R2-x2 , for 0xH and d(x)=R, for x>H

  • the kinetics of OsO4 diffusion can be described by Fick’s second law,

  • the binding of OsO4 to the tissue is irreversible and can be described by first-order reaction kinetics

  • the density of available OsO4 binding sites is limited,

  • the unmasking of additional binding sites is irreversible, follows first-order reaction kinetics and depends on the local concentration of free OsO4

  • the tissue expands isotropically and can be modeled by a monomolecular, saturating growth function, which induces a local velocity field v(x,t)=-x(1+g(t)/100)t=-xa2100τexp(-tτ).

Because the projected intensity also includes any X-ray absorbing components surrounding the sample, we have to consider the change in free osmium density outside of the sample o(x,t) for fitting the model to the experimental data:

(7) o(x,t)t=-xo(x,t)v(x,t)

The initial conditions (t=0) for the system of coupled partial differential Equations 2–5; 7 were chosen as follows:

  • concentration of free OsO4 in the sample: f(x0,t=0)=C0R and f(x>0,t=0)=0

  • density of bound OsO4: b(x,t=0)=0

  • density of available binding sites: s(x0,t=0)=0 and s(x>0,t=0)=S0d(x)

  • density of masked binding sites: m(x0,t=0)=0 and m(x>0,t=0)=M0d(x)

  • concentration of free OsO4 surrounding the sample: o(x0,t=0)=C0(D-R) and o(x>0,t=0)=C0(Dd(x)), where D=25mm is the vial diameter.

The boundary conditions of the system of coupled partial differential Equations 2–5; 7 for x0,L with maximal depth L=3mm were chosen as follows:

  • concentration of free OsO4: f(x=0,t)=C0R and f(x=L,t)=0

  • density of bound OsO4: b(x=0,t)=0 and b(x=L,t)=B0d(L)

  • density of available binding sites: s(x=0,t)=0 and s(x=L,t)=S0d(L)

  • density of masked binding sites: m(x=0,t)=0 and m(x=L,t)=M0d(L)

  • concentration of free OsO4 surrounding the sample: o(x=0,t)=C0(D-R) and o(x=L,t)=C0(D-d(L))

  • No flux beyond x=L: f(x=L,t)x=s(x=L,t)x=m(x=L,t)x=o(x=L,t)x=0

The MATLAB function pdepe was used to solve/simulate the initial-boundary value problem described above.

Diffusion-reaction-advection model fitting

C0 was calculated by averaging uib(x,t) over x[-150,-50] μm from the sample surface (i.e. outside of the sample) and t[0,1200] min and dividing it by the vial diameter D=25mm and by the intensity-to-concentration scaling factor FI2C=214.6308mm3mmolmm . FI2C was calculated from the linear fit in Figure 1c. In each fitting iteration, simulated data was produced according to the diffusion-reaction-advection model above. Within the fitting procedure, FI2C(o(x,t)+f(x,t)+b(x,t)) was fitted to uib(x,t) in a least square displacement sense within the domain x[100,700] μm from the sample surface and t[0,1200] min by using the trust-region-reflective algorithm of the nonlinear least-squares solver lsqnonlin in MATLAB to optimize the N=6 model parameters De , kon , kunmask , S0 , M0 and H. These parameters were fitted to the average experimentally measured osmium accumulation from n=13 samples within the domains x[100,700] μm and t[0,1200] min. To evaluate the quality of the model fits, the residual standard error was calculated as SEres=(x,t)uib(x,t)-(FI2C(o(x,t)+f(x,t)+b(x,t)))2df , were df corresponds to the degree of freedom.

Estimation of the average double-bond/ osmium binding site density in the mouse brain

The average phospholipid density in the mouse brain is 56,241±4000 nmol/g (Barceló-Coblijn et al., 2007) and for an average mouse brain volume of 508.91±23.42 mm3 (Badea et al., 2007) with a weight of 0.427±0.024 g (The Jackson Laboratory, 2021), a cubic millimeter brain tissue contains about 47.19 nmol phospholipids. Lipids in the brain carry on average 2.5 double bonds (Fitzner et al., 2020) which results in 117.975 nmol/mm3 double bonds to which osmium tetroxide could bind through osmium mono- or di-ester formation. Our proposed model predicts 2.999 times more binding sites (353.774±27.460 nmol/mm3), even without unmasking (Figure 2—figure supplement 4).

Data availability

The code to analyze the X-ray projection images and to model the accumulation of heavy metals can be found on https://github.com/adwanner/XrayAssistedStaining, (copy archived at swh:1:rev:62c9930d08bbdd1602614d515cd20b1cb1f89d15). All X-ray data is available for download on https://www.ebi.ac.uk/empiar/ with dataset ID EMPIAR-10782.

The following data sets were generated
    1. Ströh S
    2. Hammerschmith EW
    3. Tank DW
    4. Seung HS
    5. Wanner AA
    (2022) EMPIAR
    ID EMPIAR-10782. In situ X-ray assisted electron microscopy staining for large biological samples.

References

    1. Burkl W
    2. Schiechl H
    (1968) A study of osmium tetroxide fixation
    The Journal of Histochemistry and Cytochemistry 16:157–161.
    https://doi.org/10.1177/16.3.157
    1. Finck H
    (1960) Epoxy resins in electron microscopy
    The Journal of Biophysical and Biochemical Cytology 7:27–30.
    https://doi.org/10.1083/jcb.7.1.27
    1. Fox CH
    2. Johnson FB
    3. Whiting J
    4. Roller PP
    (1985) Formaldehyde fixation
    The Journal of Histochemistry and Cytochemistry 33:845–853.
    https://doi.org/10.1177/33.8.3894502
    1. Maaløe O
    2. Andersen B
    (1956)
    Symposium of the Society for General Microbiology
    Bacterial anatomy, Symposium of the Society for General Microbiology, Cambridge, University Press.
    1. Ren Y
    2. Kruit P
    (2016) Transmission electron imaging in the delft multibeam scanning electron microscope 1
    Journal of Vacuum Science & Technology B, Nanotechnology and Microelectronics 34:06KF02.
    https://doi.org/10.1116/1.4966216

Decision letter

  1. Andreas T Schaefer
    Reviewing Editor; The Francis Crick Institute, United Kingdom
  2. Catherine Dulac
    Senior Editor; Harvard University, United States

Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "In situ X-ray assisted electron microscopy staining for large biological samples" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Catherine Dulac as the Senior Editor. The reviewers have opted to remain anonymous.

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

As you can see from the detailed reviews, all reviewers (and myself as well) felt that this work is potentially an important contribution.

However, as the central claim of the work is to use CT as a benchmark for the quality of staining for EM, it will be critical to directly relate those two measures, i.e. to perform EM imaging on previously CT imaged samples and relate the two measures of contrast to each other. This will be the critical experimental dataset for revision. There are other valuable comments below. Please edit the manuscript accordingly.

Reviewer #1 (Recommendations for the authors):

The authors are in a very good position to disambiguate current divergences in the initial steps of osmium staining of soft tissues for electron microscopy imaging purposes.

First, it would be very informative to assess the effect of a different iron salt commonly used to reduce osmium: potassium ferricyanide (e.g. by adding the experimental group 2% OsO4, 2.5% potassium ferricyanide to figure panels 1e-f).

Second, ref. Hua et al., 2015 implemented the osmium reduction posterior to the osmium incubation, by adding potassium ferrocyanide without intermediate washes. The authors have already explored this staining protocol in Supp Figure 5: samples being stained with OsO4 for 22h were transferred to a potassium ferrocyanide solution. It would be very informative to view those results in summary (ie averaging the 20h timepoint datasets across samples) along with the osmium and reduced osmium baselines shown in Figure 1e, and to complement the insight on sample expansion for that protocol, showing it against the same baselines, as in Figure 1f. An additional control group could be beneficial here: transferring the osmium-stained samples to a water vial, to account for decrease of staining by dilution and diffusion of unbound osmium. In this case it would not be expected to see sample degradation regardless of the incubation time.

Third, after staining with osmium or reduced osmium, samples are meant to be washed with water before proceeding to the next staining step, aiming to remove any unbound osmium from the media. The goal of optimizing the first staining step would therefore be to maximize the bound osmium remaining in the sample after the water wash(es). To that end, it would be interesting to see the evolution of intensity vs depth for samples placed in staining solution for 20h and then placed in water, for at least the most typical stains explored (osmium and reduced osmium).

Finally, it would be very useful if the authors released the code used for pre-processing single X-ray projections. Given the popularity of benchtop X-ray uCT devices, doing so could simplify and accelerate the adoption of this monitoring method of sample preparation protocols by the user community.

These additions could slightly improve the scope and completeness of an already robust and solid study.

Reviewer #2 (Recommendations for the authors):

Reference, page 4: Claude and Fullam (1945) used formaldehydes + osmium fixation. Only Sabatini introduced glutaraldehyde. It might be good to slightly modify the text.

Page 9: the authors refer to figure 1b to comment the staining after 20hours. Yet, fig 1b shows example at 10 hours max.

Page 9: "we found that for reduced osmium the accumulation of heavy metals peaks at a depth between 300- 800 μm (Figure 1e, Supplementary Figure 2), whereas the tissue above or below that depth is stained less." Could this be due to the anatomy/heterogeneity of the various cortical layers?

Page 9: about the shrinking/expansion of the tissue upon fixation with osmium or reduced osmium. This is very interesting as it could explain a few artifacts observed with several sample types. Yet in this experiment, the cortex biopsy is fitted inside a hole, of the same shape, made inside the sylgard support layer. It is thus possible that the tissue is constrained and can't expand in the lateral directions. Would this enhance the vertical expansion? And therefore such expansion would be less than 5%? Please comment on this.

In line with this, it might be good to adapt the drawings of figure 1a. The way the biopsy is displayed, it seems as if it was just placed onto the vial bottom, when in fact is is embedded in sylgard.

Page 12: about the osmium density: if I am not wrong these numbers are estimates, extrapolated from many assumptions re the density of unsaturated lipids, double bonds etc ... this should be made clear in the main text and not only in the methods. Even though the authors made their best possible efforts to describe the reactions within the tissues and to estimate the amount of binding sites, I doubt we have a solid proof of the actual concentration of bound or free osmium in this experiment.

Page 12: the authors write "While we assume that the tissue density is homogeneous across the cortical thickness ..." I don't think this is a fair assumption. The cortex is organized in layers of different cellular densities.

Page 14: discussion. what is not addressed in this paper is the secondary osmium blackening during dehydration. This has been described previously and attributed to unbound osmium and later reduction during the process. IT might be good to elaborate a bit on this. Of course, this leads to the comments made above, which highlight the need to check the samples in EM, i.e. to go, for some of them, through the full sample preparation.

Page 16: methods.

- The concentration ranges are surprising. Can you explain: "fixative solution 1.3 - 1.5% glutaraldehyde (GA) and 2.5 - 2.6% paraformaldehyde (PFA) in 0.14 - 0.15M cacodylate buffer with 2.0 - 2.1mM CaCl2"

- What is the rationale behind these incubation times. Intuitively, one would expect to keep whole brains longer in the fixative? "The samples were kept in fixative solution at 4 °C for 36 hours (4-mm biopsy punches) or for 12 hours (whole mouse brain)."

- About the sylgard molds: can you comment on the osmium infiltration. The sides of the biopsy might not be as well exposed as the top.

Page 17: X-ray microscopy: the working temperature is not indicate, though it is a very important parameter for the osmium fixation.

Reviewer #3 (Recommendations for the authors):

The manuscript would be strengthened if the authors compared the actual membrane contrast observed in cross-sections from their samples to demonstrate if these intensities can be quantitatively interpreted (beyond the description of staining gradients). I suggest to cut thin sections from a reduced osmium and osmium tissue cylinder and quantify the staining intensity/contrast across the 4mm diameter to compare with the x-ray profiles.

[Editors’ note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "In situ X-ray assisted electron microscopy staining for large biological samples" for further consideration by eLife. Your revised article has been evaluated by Catherine Dulac (Senior Editor) and a Reviewing Editor.

The manuscript has been improved but there are some remaining issues that need to be addressed by small textual adjustments, as outlined below:

Please take note of the comments made by reviewer 3 regarding the strict correlation of X-ray and EM measure. Please adjust your textual description accordingly. While e.g. Mikula and Denk do provide some indication that this is the case (as does your work here), I would not want to give the impression that this is a question that has experimentally conclusively and quantitatively been settled. There is still room for direct experimental comparison for X-rays with EM under different conditions (although not necessarily in this manuscript).

In short, please tune down the corresponding comments and add a sentence or two in the discussion to mention the caveats and limitations discussed below and maybe point towards further such experiments. In this context it could be useful to expand a bit on the physics background and mention how X-ray contrast will depend e.g. on photon energy and that e.g. at the absorption edges the relationship of X-ray absorption and signal in the EM can become more complicated.

Reviewer #2 (Recommendations for the authors):

The revised version of the paper is efficiently addressing the points that I have raised and I am satisfied with the rebuttal answers as well.

In my opinion, the paper is ready for publication and expect the scientific community, especially the EM community to greatly benefit from this contribution.

Reviewer #3 (Recommendations for the authors):

My recommendation to the authors was to correlate their observed X-ray intensities to what is actually observed by EM in thin sections. In response, they cited Figure 4b from Mikula and Denk, 2015 as evidence that the X-ray reconstructed pixel intensities are correlated with what is observed in EM sections and state "Mikula and Denk (2015) showed that the pixel intensities of serial-section EM images and the intensity of the corresponding reslice are highly correlated."

I think this is an overstatement for several reasons:

1) As far as I can tell, Mikula and Denk Figure 4b is a n=1 measurement

2) The figure in question doesn't report a correlation coefficient, indeed their figure shows a peak SEM normalized intensity at the core of their sample and the X-ray normalized intensity has a peak at the edge of the sample

3) Mikula and Denk themselves are far more cautious in discussing this: "…X-ray microcomputed tomography, which exhibits an image contrast similar to that of SEM (Figure 4b). This allowed us to quickly test an embedded brain for defects and distortion before deciding whether to finish preparing it for electron microscopy." Mikula and Denk were clearing using X-ray reconstructions to check for gross defects/cracks, not to make quantitative comparisons.

Stroh et al., do provide a comparison of X-ray vs. EM intensities (new Figure 1 —figure supplement 1) with the caveat that the X-ray intensities are measured from a projection instead of a computed slice reconstruction. Perhaps due to this caveat the intensities do not appear to be very well correlated (although it would be preferable to actually show the correlation scatter plot and the measured correlation coeff). In light of this (and the above comments regarding the Mikula and Denk data), I remain unsure to what degree the measured X-ray intensities that form the foundation of the manuscript are informative about actual contrast in EM sections.

My original point was related to the author’' Figure 1e: the figure shows that reduced osmium stains tissue with a maximum intensity of ~1.1 (a.u.) compared to osmium alone at ~0.9 (a.u.). Because these intensities are not calibrated against the appearance of the staining in EM sections, their interpretation is limited. The authors state they were not able to test whether this particular finding is consistent with the intensity/contrast of EM sections due to embedding issues of reduced osmium samples.

https://doi.org/10.7554/eLife.72147.sa1

Author response

Essential revisions:

As you can see from the detailed reviews, all reviewers (and myself as well) felt that this work is potentially an important contribution.

However, as the central claim of the work is to use CT as a benchmark for the quality of staining for EM, it will be critical to directly relate those two measures, i.e. to perform EM imaging on previously CT imaged samples and relate the two measures of contrast to each other. This will be the critical experimental dataset for revision. There are other valuable comments below. Please edit the manuscript accordingly.

We added Figure 1 —figure supplement 3 showing the ultrastructural preservation after 20 hours of incubation in OsO4 and Figure 1 —figure supplement 1 that shows the covariation of the EM and X-ray pixel intensities. Note, in the presented study we acquired X-ray projection images that represent the cumulative tissue and heavy metal density along the direction of projection through a 4 mm thick tissue punch. Therefore our X-ray projection pixel intensities cannot directly be compared to the EM pixel intensity of a thin section. As has been shown previously (see for example Figure 4b in Mikula and Denk 2015), in computed X-ray tomograms the intensity scales linearly with EM intensity, if the pixel intensity in an EM section is compared to the corresponding reslice of a computed X-ray tomogram.

Reviewer #1 (Recommendations for the authors):

The authors are in a very good position to disambiguate current divergences in the initial steps of osmium staining of soft tissues for electron microscopy imaging purposes.

First, it would be very informative to assess the effect of a different iron salt commonly used to reduce osmium: potassium ferricyanide (e.g. by adding the experimental group 2% OsO4, 2.5% potassium ferricyanide to figure panels 1e-f).

We added Figure 1 —figure supplement 7 and 8 as well as Figure 2 —figure supplement 3 that show the spatio-temporal accumulation of osmium reduced with 2.5% potassium ferricyanide K3[Fe(CN)6]. The accumulation of heavy metals for the buffered solution with 2% OsO4 and 2.5% potassium ferricyanide is not as homogeneous as for 2% OsO4 (Figure 1 —figure supplement 8a). But in contrast to the experiments where osmium was reduced with 2.5% potassium ferrocyanide, the sample does not shrink if incubated in buffered solution of 2% osmium reduced with 2.5% potassium ferricyanide (Figure 1 —figure supplement 8b). In contrast to potassium ferrocyanide, no ‘washout’ effect was observed if the sample was placed in buffered solution of 2.5% potassium ferricyanide after 20hrs of incubation in buffered 2% OsO4 (Figure 1 —figure supplement 8e).

Second, ref. Hua et al., 2015 implemented the osmium reduction posterior to the osmium incubation, by adding potassium ferrocyanide without intermediate washes. The authors have already explored this staining protocol in Supp Figure 5: samples being stained with OsO4 for 22h were transferred to a potassium ferrocyanide solution. It would be very informative to view those results in summary (ie averaging the 20h timepoint datasets across samples) along with the osmium and reduced osmium baselines shown in Figure 1e, and to complement the insight on sample expansion for that protocol, showing it against the same baselines, as in Figure 1f.

We added a corresponding Figure 1 —figure supplement 8.

An additional control group could be beneficial here: transferring the osmium-stained samples to a water vial, to account for decrease of staining by dilution and diffusion of unbound osmium. In this case it would not be expected to see sample degradation regardless of the incubation time.

Third, after staining with osmium or reduced osmium, samples are meant to be washed with water before proceeding to the next staining step, aiming to remove any unbound osmium from the media. The goal of optimizing the first staining step would therefore be to maximize the bound osmium remaining in the sample after the water wash(es). To that end, it would be interesting to see the evolution of intensity vs depth for samples placed in staining solution for 20h and then placed in water, for at least the most typical stains explored (osmium and reduced osmium).

We added Figure 2 —figure supplement 4 and 5 that show the dynamics of heavy metal washout in double distilled H2O for 22 hours after 22 hrs of incubation in 2% buffered OsO4. First, the sample was quickly flushed with double distilled H2O four times to completely remove any buffered OsO4 solution. Subsequently, the sample was immersed in double distilled H2O for 22 hrs. The accumulation of heavy metals decreases by about 6% in a depth of 100-1200μm.

Our hypothesis is that this can be explained by unbound osmium that diffuses out of the sample and by the slight expansion of the sample (Figure 2 —figure supplement 5b). But in contrast to the washout effect in K3[Fe(CN)6], the reduction in heavy metal density appears to be small.

Finally, it would be very useful if the authors released the code used for pre-processing single X-ray projections. Given the popularity of benchtop X-ray uCT devices, doing so could simplify and accelerate the adoption of this monitoring method of sample preparation protocols by the user community.

The python and MATLAB code for processing the X-ray data is now available on https://github.com/adwanner/XrayAssistedStaining.

Reviewer #2 (Recommendations for the authors):

Reference, page 4: Claude and Fullam (1945) used formaldehydes + osmium fixation. Only Sabatini introduced glutaraldehyde. It might be good to slightly modify the text.

We corrected the text and citations accordingly.

Page 9: the authors refer to figure 1b to comment the staining after 20hours. Yet, fig 1b shows example at 10 hours max.

We corrected it accordingly.

Page 9: "we found that for reduced osmium the accumulation of heavy metals peaks at a depth between 300- 800 μm (Figure 1e, Supplementary Figure 2), whereas the tissue above or below that depth is stained less." Could this be due to the anatomy/heterogeneity of the various cortical layers?

Indeed, the effects of reduced osmium accumulation could depend at least in parts of the tissue composition, i.e. the anatomy/heterogeneity of the cortical layers in the tissue block. It is, however, important to note that this effect is not present in non-reduced osmium immersed brain tissue samples that have the same composition. We added the sentences that ‘No such band is present in tissue immersed in non-reduced OsO4.’

Page 9: about the shrinking/expansion of the tissue upon fixation with osmium or reduced osmium. This is very interesting as it could explain a few artifacts observed with several sample types. Yet in this experiment, the cortex biopsy is fitted inside a hole, of the same shape, made inside the sylgard support layer. It is thus possible that the tissue is constrained and can't expand in the lateral directions. Would this enhance the vertical expansion? And therefore such expansion would be less than 5%? Please comment on this.

In line with this, it might be good to adapt the drawings of figure 1a. The way the biopsy is displayed, it seems as if it was just placed onto the vial bottom, when in fact is is embedded in sylgard.

The contribution to the vertical expansion from the embedding in sylgard is expected to be minimal, because the sylgard layer was typically only 2 mm thick, whereas the height of the tissue samples was typically about 6.75 mm tall. The upper end of the sylgard layer can be seen in Figure 1b at the very bottom of the image (small arrowhead). Hence, only a small fraction of the sample sits in sylgard. We adapted Figure 1 accordingly.

That the contribution from the restriction in sylgard is small can also be seen from the following considerations:

2mm×π(a×124mm)22mm×π(124mm)2=2mm×(a1)2×π(124mm)2=2mm×(11/a)2×π(a×124mm)2

For buffered 2% OsO4 we measured an expansion of 5%, i.e. b=1.05 and according to 3. above the ‘true’ homogeneous expansion factor would have been a=1.049.

Page 12: about the osmium density: if I am not wrong these numbers are estimates, extrapolated from many assumptions re the density of unsaturated lipids, double bonds etc ... this should be made clear in the main text and not only in the methods. Even though the authors made their best possible efforts to describe the reactions within the tissues and to estimate the amount of binding sites, I doubt we have a solid proof of the actual concentration of bound or free osmium in this experiment.

The osmium density (as shown in Figure 2b) was measured using the pixel intensities which were calibrated using the concentration series shown in Figure 1c. No further assumptions are being made. In the discussion section highlighted by the reviewer we compare this measured osmium density to the density of double bonds of unsaturated lipids in the brain as estimated from numbers reported in the literature. The double bonds of unsaturated lipids are thought to be the main binding site for osmium. Comparing the measured density of osmium to the estimated density of binding sites we found that the measured osmium density is 2.999 times higher than the estimated density of double bonds in the brain. We therefore hypothesize that with long osmium incubation, 3 osmium atoms get associated with 1 double bond.

Page 12: the authors write "While we assume that the tissue density is homogeneous across the cortical thickness ..." I don't think this is a fair assumption. The cortex is organized in layers of different cellular densities.

We agree that the presented model and the underlying assumptions significantly simplify the actual bio-chemical environment. But despite these oversimplifications, the resulting model presented in Figure 2 captures the dynamics of osmium accumulation well. While modeling the tissue densities of different cortical layers in detail and calculating individual diffusion constants is beyond the presented study, we do agree that it would be very interesting to investigate how the diffusion constants differ in different tissue types. To illustrate this, we measured and modeled the staining kinetics of 2% buffered OsO4 in 4mm punches of liver tissue (Figure 2 - Figure Supplement 7). Interestingly, the effective diffusion coefficient appears to be >4X larger in liver tissue compared to brain tissue.

Page 14: discussion. what is not addressed in this paper is the secondary osmium blackening during dehydration. This has been described previously and attributed to unbound osmium and later reduction during the process. IT might be good to elaborate a bit on this. Of course, this leads to the comments made above, which highlight the need to check the samples in EM, i.e. to go, for some of them, through the full sample preparation.

We added Figure 1 - Figure Supplement 3that shows the ultrastructural preservation after 20 hours of osmium incubation. In addition, we performed additional experiments in which we measured how the osmium density changes during the H2O washing step after the OsO4 incubation. These results of the wash-out experiments are shown in Figure 2 - Figure Supplement 4 and 5. We do find a slight decrease in heavy metal density that can likely be explained on one hand by the wash out of unbound osmium from the sample as well as by expansion of the tissue which reduces the local heavy metal density.

Page 16: methods.

- The concentration ranges are surprising. Can you explain: "fixative solution 1.3 - 1.5% glutaraldehyde (GA) and 2.5 - 2.6% paraformaldehyde (PFA) in 0.14 - 0.15M cacodylate buffer with 2.0 - 2.1mM CaCl2"

Some of the mice have been perfused by different experimenters that used slightly different effective fixative formulations.

- What is the rationale behind these incubation times. Intuitively, one would expect to keep whole brains longer in the fixative? "The samples were kept in fixative solution at 4 °C for 36 hours (4-mm biopsy punches) or for 12 hours (whole mouse brain)."

The 36 hours of post-perfusion incubation in fixative corresponds to the duration that we and other labs typically use for fixing brain tissue. The 12 hours post-perfusion incubation in fixative of the whole mouse brain might indeed be on the lower end for what might be necessary for thoroughly fixing a whole mouse brain after perfusion.

- About the sylgard molds: can you comment on the osmium infiltration. The sides of the biopsy might not be as well exposed as the top.

The sylgard layer was typically only 2 mm thick, whereas the height of the tissue samples was typically about 6.75 mm tall. The upper end of the sylgard layer can be seen in Figure 1b at the very bottom of the image (small arrowhead). Hence, only a small fraction of the sample sits in sylgard and the measurements have been taken several millimeters away from the sylgard, therefore we do not expect any major impact of the sylgard onto the osmium infiltration measurements. We adapted Figure 1 accordingly.

Page 17: X-ray microscopy: the working temperature is not indicate, though it is a very important parameter for the osmium fixation.

All staining and washing steps have indeed been done at room temperature. We added that information accordingly in the sample preparation section of the methods.

Reviewer #3 (Recommendations for the authors):

The manuscript would be strengthened if the authors compared the actual membrane contrast observed in cross-sections from their samples to demonstrate if these intensities can be quantitatively interpreted (beyond the description of staining gradients). I suggest to cut thin sections from a reduced osmium and osmium tissue cylinder and quantify the staining intensity/contrast across the 4mm diameter to compare with the x-ray profiles.

We added Figure 1 —figure supplement 3showing the ultrastructural preservation after 20 hours of incubation in OsO4 and Figure 1 —figure supplement 1 that shows the covariation of the EM and X-ray pixel intensities. Note, in the presented study we acquired X-ray projection images that represent the cumulative tissue and heavy metal density along the direction of projection through a 4 mm thick tissue punch. Therefore our X-ray projection pixel intensities cannot directly be compared to the EM pixel intensity of a thin section. As has been shown previously (see for example Figure 4b in Mikula and Denk 2015), in computed X-ray tomograms the intensity scales linearly with EM intensity, if the pixel intensity in an EM section is compared to the corresponding reslice of a computed X-ray tomogram.

We also tried to acquire EM images of tissue punches stained with reduced osmium, but unfortunately the resin could not infiltrate the samples sufficiently well to allow us to collect serial sections.

[Editors’ note: further revisions were suggested prior to acceptance, as described below.]

The manuscript has been improved but there are some remaining issues that need to be addressed by small textual adjustments, as outlined below:

Please take note of the comments made by reviewer 3 regarding the strict correlation of X-ray and EM measure. Please adjust your textual description accordingly. While e.g. Mikula and Denk do provide some indication that this is the case (as does your work here), I would not want to give the impression that this is a question that has experimentally conclusively and quantitatively been settled. There is still room for direct experimental comparison for X-rays with EM under different conditions (although not necessarily in this manuscript).

In short, please tune down the corresponding comments and add a sentence or two in the discussion to mention the caveats and limitations discussed below and maybe point towards further such experiments. In this context it could be useful to expand a bit on the physics background and mention how X-ray contrast will depend e.g. on photon energy and that e.g. at the absorption edges the relationship of X-ray absorption and signal in the EM can become more complicated.

Reviewer #3 (Recommendations for the authors):

My recommendation to the authors was to correlate their observed X-ray intensities to what is actually observed by EM in thin sections. In response, they cited Figure 4b from Mikula and Denk, 2015 as evidence that the X-ray reconstructed pixel intensities are correlated with what is observed in EM sections and state "Mikula and Denk (2015) showed that the pixel intensities of serial-section EM images and the intensity of the corresponding reslice are highly correlated."

I think this is an overstatement for several reasons:

1) As far as I can tell, Mikula and Denk Figure 4b is a n=1 measurement

2) The figure in question doesn't report a correlation coefficient, indeed their figure shows a peak SEM normalized intensity at the core of their sample and the X-ray normalized intensity has a peak at the edge of the sample

3) Mikula and Denk themselves are far more cautious in discussing this: "…X-ray microcomputed tomography, which exhibits an image contrast similar to that of SEM (Figure 4b). This allowed us to quickly test an embedded brain for defects and distortion before deciding whether to finish preparing it for electron microscopy." Mikula and Denk were clearing using X-ray reconstructions to check for gross defects/cracks, not to make quantitative comparisons.

Stroh et al., do provide a comparison of X-ray vs. EM intensities (new Figure 1 —figure supplement 1) with the caveat that the X-ray intensities are measured from a projection instead of a computed slice reconstruction. Perhaps due to this caveat the intensities do not appear to be very well correlated (although it would be preferable to actually show the correlation scatter plot and the measured correlation coeff). In light of this (and the above comments regarding the Mikula and Denk data), I remain unsure to what degree the measured X-ray intensities that form the foundation of the manuscript are informative about actual contrast in EM sections.

My original point was related to the authors' Figure 1e: the figure shows that reduced osmium stains tissue with a maximum intensity of ~1.1 (a.u.) compared to osmium alone at ~0.9 (a.u.). Because these intensities are not calibrated against the appearance of the staining in EM sections, their interpretation is limited. The authors state they were not able to test whether this particular finding is consistent with the intensity/contrast of EM sections due to embedding issues of reduced osmium samples.

Thank you very much for the constructive feedback on the revised manuscript. We agree with Reviewer 3 that adding some quantitative evidence for the correlation between electron microscopy images and X-ray images would strengthen the revised manuscript. On one hand, we adjusted the wording regarding the similarity/correlation between EM and micro-CT in the introduction. On the other hand, we added a note of

caution to the discussion that for a qualitative comparison between EM and micro-CT appearance the two imaging modalities need to be calibrated and the precise image acquisition parameters need to be considered.

In addition, we added additional data comparing EM and micro-CT appearance similar to Mikula et al., to Figure 1 —figure supplement 1. In two different en bloc stained brain samples we show the correspondence between the EM and micro-CT appearance. As suggested by Reviewer 3 we included the (indeed high) correlation statistics and corresponding scatter plots.

To conclude, we have tuned down the claims regarding the correspondence between EM and X-ray imaging and at the same time added more examples and quantifications supporting the claims. We believe that these changes and additional information have improved the revised manuscript.

https://doi.org/10.7554/eLife.72147.sa2

Article and author information

Author details

  1. Sebastian Ströh

    Princeton Neuroscience Institute, Princeton University, Princeton, United States
    Contribution
    Conceptualization, Data curation, Validation, Investigation, Methodology, Writing – original draft, Writing – review and editing
    Contributed equally with
    Eric W Hammerschmith
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8128-2617
  2. Eric W Hammerschmith

    Princeton Neuroscience Institute, Princeton University, Princeton, United States
    Contribution
    Conceptualization, Data curation, Validation, Investigation, Methodology, Writing – original draft, Writing – review and editing
    Contributed equally with
    Sebastian Ströh
    Competing interests
    No competing interests declared
  3. David W Tank

    Princeton Neuroscience Institute, Princeton University, Princeton, United States
    Contribution
    Formal analysis, Supervision, Funding acquisition, Validation, Methodology, Writing – review and editing
    Competing interests
    is an inventor of US Patent Application 16/681,028
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9423-4267
  4. H Sebastian Seung

    Princeton Neuroscience Institute, Princeton University, Princeton, United States
    Contribution
    Formal analysis, Supervision, Funding acquisition, Validation, Methodology, Writing – review and editing
    Competing interests
    is an inventor of US Patent Application 16/681,028
  5. Adrian Andreas Wanner

    1. Princeton Neuroscience Institute, Princeton University, Princeton, United States
    2. Paul Scherrer Institute, Villigen, Switzerland
    Contribution
    Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    adrian.wanner@psi.ch
    Competing interests
    is an inventor of US Patent Application 16/681,028
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5864-8577

Funding

CV Starr Fellowship in Neuroscience by the Princeton University

  • Adrian Andreas Wanner

National Institutes of Health (NS104648)

  • David W Tank

National Institutes of Health (1R01EY027036)

  • H Sebastian Seung

National Institutes of Health (U01NS090562)

  • H Sebastian Seung

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Michael Bozlar for advice on potassium ferrocyanide and osmium chemistry, Stefan N Oline for excellent assistance with animal perfusions and Alyssa M Wilson and Zhihao Zheng for comments on the manuscript. The authors acknowledge the use of Princeton’s Imaging and Analysis Center, which is partially supported through the Princeton Center for Complex Materials (PCCM), a National Science Foundation (NSF)-MRSEC program (DMR-2011750). AAW acknowledges support by the CV Starr Fellowship in Neuroscience by the Princeton University. This work was supported by the National Institutes of Health (NIH) grants U19 NS104648, NIH/NEI (1R01EY027036 to HSS), NIH/NINDS (1R01NS104926 to HSS), NINDS/NIH (U01NS090562 to HSS).

Ethics

Animal use procedures were approved by the Princeton University Institutional Animal Care and Use Committee (protocol number 2000) and carried out in accordance with National Institutes of Health standards (AAALAC International Institutional Number: Unit #1001, PHS assurance ID D16-00273).

Senior Editor

  1. Catherine Dulac, Harvard University, United States

Reviewing Editor

  1. Andreas T Schaefer, The Francis Crick Institute, United Kingdom

Version history

  1. Preprint posted: June 20, 2021 (view preprint)
  2. Received: July 13, 2021
  3. Accepted: October 18, 2022
  4. Accepted Manuscript published: October 20, 2022 (version 1)
  5. Version of Record published: November 15, 2022 (version 2)

Copyright

© 2022, Ströh, Hammerschmith et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 725
    Page views
  • 114
    Downloads
  • 5
    Citations

Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Sebastian Ströh
  2. Eric W Hammerschmith
  3. David W Tank
  4. H Sebastian Seung
  5. Adrian Andreas Wanner
(2022)
In situ X-ray-assisted electron microscopy staining for large biological samples
eLife 11:e72147.
https://doi.org/10.7554/eLife.72147

Further reading

    1. Cell Biology
    2. Neuroscience
    Elisabeth Jongsma, Anita Goyala ... Collin Yvès Ewald
    Research Article Updated

    The amyloid beta (Aβ) plaques found in Alzheimer’s disease (AD) patients’ brains contain collagens and are embedded extracellularly. Several collagens have been proposed to influence Aβ aggregate formation, yet their role in clearance is unknown. To investigate the potential role of collagens in forming and clearance of extracellular aggregates in vivo, we created a transgenic Caenorhabditis elegans strain that expresses and secretes human Aβ1-42. This secreted Aβ forms aggregates in two distinct places within the extracellular matrix. In a screen for extracellular human Aβ aggregation regulators, we identified different collagens to ameliorate or potentiate Aβ aggregation. We show that a disintegrin and metalloprotease a disintegrin and metalloprotease 2 (ADM-2), an ortholog of ADAM9, reduces the load of extracellular Aβ aggregates. ADM-2 is required and sufficient to remove the extracellular Aβ aggregates. Thus, we provide in vivo evidence of collagens essential for aggregate formation and metalloprotease participating in extracellular Aβ aggregate removal.

    1. Computational and Systems Biology
    2. Neuroscience
    Marjorie Xie, Samuel P Muscinelli ... Ashok Litwin-Kumar
    Research Article Updated

    The cerebellar granule cell layer has inspired numerous theoretical models of neural representations that support learned behaviors, beginning with the work of Marr and Albus. In these models, granule cells form a sparse, combinatorial encoding of diverse sensorimotor inputs. Such sparse representations are optimal for learning to discriminate random stimuli. However, recent observations of dense, low-dimensional activity across granule cells have called into question the role of sparse coding in these neurons. Here, we generalize theories of cerebellar learning to determine the optimal granule cell representation for tasks beyond random stimulus discrimination, including continuous input-output transformations as required for smooth motor control. We show that for such tasks, the optimal granule cell representation is substantially denser than predicted by classical theories. Our results provide a general theory of learning in cerebellum-like systems and suggest that optimal cerebellar representations are task-dependent.