Eliana Battistella1, Jan Schniete1, Katrina Wesencraft1, Juan F Quintana2, Gail McConnell1. 1. Department of Physics, SUPA, University of Strathclyde, 107 Rottenrow East, Glasgow G4 0NG, UK. 2. Wellcome Centre for Integrative Parasitology, Institute of Biodiversity Animal Health and Comparative Medicine, University of Glasgow, Glasgow, UK.
Abstract
Rapid imaging of large biological tissue specimens such as ultrathick sections of mouse brain cannot easily be performed with a standard microscope. Optical mesoscopy offers a solution, but thus far imaging has been too slow to be useful for routine use. We have developed two different illuminators for light-sheet mesoscopy with the Mesolens and we demonstrate their use in high-speed optical mesoscale imaging of large tissue specimens. The first light-sheet approach uses Gaussian optics and is straightforward to implement. It provides excellent lateral resolution and high-speed imaging, but the axial resolution is poor. The second light-sheet is a more complex Airy light-sheet that provides sub-cellular resolution in three dimensions that is comparable in quality to point-scanning confocal mesoscopy, but the light-sheet method of illuminating the specimen reduces the imaging time by a factor of 14. This creates new possibilities for high-content, higher-throughput optical bioimaging at the mesoscale.
Rapid imaging of large biological tissue specimens such as ultrathick sections of mouse brain cannot easily be performed with a standard microscope. Optical mesoscopy offers a solution, but thus far imaging has been too slow to be useful for routine use. We have developed two different illuminators for light-sheet mesoscopy with the Mesolens and we demonstrate their use in high-speed optical mesoscale imaging of large tissue specimens. The first light-sheet approach uses Gaussian optics and is straightforward to implement. It provides excellent lateral resolution and high-speed imaging, but the axial resolution is poor. The second light-sheet is a more complex Airy light-sheet that provides sub-cellular resolution in three dimensions that is comparable in quality to point-scanning confocal mesoscopy, but the light-sheet method of illuminating the specimen reduces the imaging time by a factor of 14. This creates new possibilities for high-content, higher-throughput optical bioimaging at the mesoscale.
Light-sheet fluorescence microscopy (LSFM) has gained popularity since its introduction in 1993 (Voie et al., 1993). In LSFM, the sample is illuminated with a thin sheet of light at a 90° angle to the detection path such that the illuminated plane coincides with the focal plane of the detection objective. This method of illuminating the sample limits fluorescence excitation to a thin layer in the specimen and avoids the creation of out-of-focus background fluorescence. LSFM has become the tool of choice in developmental biology for imaging live specimens and sub-millimeter scale fixed and optically cleared specimens (Krzic et al., 2012; Bassi et al., 2015), out-performing confocal laser scanning microscopy (CLSM) in speed (Schmid et al., 2013) and is a gentler imaging technique in that it causes less photodamage and less phototoxicity (Jemielita et al., 2013).The light-sheet is commonly generated in one of the two ways. The first method is to pass light through a cylindrical lens which focuses the beam in one dimension, thus creating a thin sheet of light at the focal plane of the detection objective lens (Engelbrecht et al., 2010; Yang et al., 2014). The second method uses a narrow excitation beam in the same plane but of a symmetrical cross-section that is brought to a focus in the center of the field of view (FOV) (Mertz and Kim, 2010; Fahrbach et al., 2013; Tomer et al., 2015). This focused beam is scanned across the FOV and a camera with a global or rolling shutter is used to acquire the image only from the vicinity of the illumination beam at any instant. This approach is called digitally scanned laser light-sheet microscopy (DSLM) (Keller et al., 2008, 2010). An evolution of this method is the confocal line detection light-sheet in which only one line of pixels is recorded by synchronizing detection and illumination (Fahrbach et al., 2013). These methods are known to reduce stripe artifacts and increase the axial resolution if used in combination with non-diffractive beams (Kafian et al., 2020; Ricci et al., 2022).LSFM has been successful with small specimens but there are currently few options for the three-dimensional imaging of large specimens with high spatial resolution. The mesoSPIM initiative (Voigt et al., 2019) allows the imaging of specimens up to 1 cm3 with an isotropic resolution of 6.5 μm. The mesoSPIM system provides a variable FOV of between 2 and 21 mm via a zoom macroscope, thus large specimens such as a whole mouse brain can be imaged in toto at low resolution. There are also two-photon mesoscale but non-light-sheet imaging methods (Stirman et al., 2016; Sofroniew et al., 2016; Demas et al., 2021; Yu et al., 2021) for in-vivo microscopy. These techniques are intended for the functional imaging of sub-regions of large tissue volumes that are either pre-defined by the user or determined stochastically. In these methods, the acquisition speed is usually prioritized over the acquisition volume but the working distances are short and so the imaging volumes are modest.The Mesolens (McConnell et al., 2016) is a complex objective with a magnification of 4× and a numerical aperture (NA) of 0.47. This gives a large improvement in image information content compared with commercial objective lenses. It is possible to image the internal structure of every cell within a 6 mm long, 3 mm thick mouse embryo or adult Drosophila if the tissue is cleared appropriately (McConnell et al., 2016; McConnell and Amos, 2018). To achieve volumetric imaging confocal point-scanning has been used, but it is very slow to acquire data. For example, confocal imaging with the Mesolens within a reduced FOV to study a 4.4 mm x 3 mm x 3 mm mouse embryo at Nyquist resolution (4 pixels/μm in xy, and 1000 images in z) with a pixel dwell time of 1 μs and a frame average = 2 takes over 117 h.It may seem obvious to combine the now common LSFM method with the high information-gathering capacity of the Mesolens. The NA of the Mesolens is such that a light-sheet of only a few microns thick is ideally preferred. However, this is not straightforward to produce because of fundamental limits in Gaussian optics. There is always a trade-off between the FOV and light-sheet thickness owing to the nature of Gaussian beam propagation (Buytaert et al., 2012):where is the beam waist radius at a given distance z from the focal point, is the beam waist radius at the focal point, and is the Rayleigh range, given byand is the distance either side of the focus along the direction of propagation where the beam waist is no larger than √2 times the waist at the focus and λ is the wavelength of the light. According to Equations (1) and (2), a light-sheet 6 mm wide cannot have a thickness of less than approximately 50 μm (Born and Wolf, 2019). We cannot reduce the Rayleigh length of the beam to reduce the beam diameter and use a scanned light-sheet method because the sensor-shifting cameras required to image the Mesolens FOV with Nyquist sampled resolution are global shutter only (Schniete et al., 2018; Shaw et al., 2021). Although the Mesolens has a usable FOV of 6 mm in diameter, our sensor-shifting camera restricts the FOV to 4.4 × 3 mm. Because of a long processing time, it takes this camera approximately 25 s to record and transfer each frame to the controlling computer.The limits of Gaussian optics in forming a light-sheet have been circumvented using non-diffractive beams such as Bessel beams and Airy beams ((Fahrbach et al., 2010; Vettenburg et al., 2014; Kafian et al., 2020). In particular, the Airy electric field envelope can be described as the function (Siviloglou and Christodoulides, 2007):where describes the Airy electric field envelope, which depends on and s, respectively, a normalized propagation distance and a transverse coordinate; is a positive parameter, typically <<1. The Airy function is one of the two independent solutions of the Airy equation . The propagation distance and the light-sheet resolution are related to the parameter (Vettenburg et al., 2014) which is experimentally determined. The theoretical FOV for an Airy light-sheet is given by (Vettenburg et al., 2014).where λ is the light source wavelength, n is the refractive index, NA is the effective numerical aperture of the illumination system, and the parameter α is experimentally determined based on the measured beam profile.The axial resolution r is given by Vettenburg et al. (2014).To date, no sufficiently thin non-diffracting beam has been described with a propagation distance matched to the wide FOV of the Mesolens.Here we describe light-sheets which we have devised, both Gaussian and Airy, that are compatible with the Mesolens for the imaging of large specimens. In spite of the slow framing rate of our high-resolution sensor-shifting camera, this has allowed us to image more than 14 times faster than point-scanning confocal mesoscopy. Our Gaussian illuminator was designed using a conventional combination of cylindrical lenses and gave a minimum light sheet thickness of 30 μm across the 3 mm FOV set by the camera detector, and is suitable for fast, low-resolution scanning. The Airy illuminator was developed from the starting point of the design of Pende et al. (2018) but we implemented an off-axis micro-translation of a Powell lens to produce an Airy beam. We obtained an Airy beam 7.8 μm thick which filled the full FOV of the Mesolens. This thickness is commensurate with the axial resolution of the Mesolens, and the axial resolution produced by the Airy illuminator was comparable to that obtained with point-scanning confocal illumination. It was not necessary to deconvolve or otherwise process the data to remove the subsidiary maxima inherent in the Airy beam because the short focal depth of the Mesolens allows these signals to be excluded optically.
Results
Measurement of light-sheet thickness
The thicknesses of Gaussian and Airy light-sheets showed in Figure 1 were measured as described under STAR Methods where it is also described the custom-made specimen holder shown in Figure 2. Figure 3 shows the light-sheet thickness measurements, with experimentally measured data and theoretical fits.
Figure 1
Schematic diagram of both light-sheet illuminators for use with the Mesolens
In both cases, the laser source is a Coherent Sapphire 488-10 CDRH laser emitting light at 488 nm with a maximum of 10 mW average power at the sample plane, with a single Gaussian mode of 0.8 mm diameter and M2 value ∼1.1.
(A) The laser beam was expanded to a diameter of 12 mm (not shown) and a flip mirror (not shown) was used to select the light path used for imaging. The light-sheet coming from the front and represented with a cross is the Airy light-sheet, the light-sheet on the right is the Gaussian illuminator. The chromatic reflector (Chroma) is placed at a 30° angle of incidence to introduce the incoherent light source for widefield epi-fluorescence imaging, while the compensation plate (Chroma) is oriented at a −30° angle to compensate for the astigmatism in the fluorescence signal. The emission filter after the Mesolens is a custom-designed 100 mm diameter filter (Chroma), with 535 ± 20 nm collection. A thermoelectric Peltier cooled camera (VNP-29MC, Vieworks) with a chip-shifting mechanism was used for image acquisition.
(B) Schematic diagram of the optical path created to produce the Gaussian light-sheet.
(C) Schematic diagram of the optical path used to create the Airy light-sheet. The red arrow shows the direction of offset for the Powell lens to create the Airy beam.
Figure 2
A computer-aided design drawing of the disassembled specimen holder and multipart chamber to support the long-term use of immersion fluid for light-sheet mesoscopy with the Mesolens
Specimens were mounted in their clearing solutions in a quartz cuvette cut to provide an inner wall height of 3 mm. The cut cuvette was then placed into a bottom plate fabricated from black acetal resin plastic. A 3D printed spacer was used to surround the cuvette on three sides for stability, leaving the fourth side open to permit light-sheet illumination, as indicated by the blue arrow. A 60 mm diameter borosilicate coverslip was placed on top of the cuvette and spacer. A top plate comprised a 30 mm diameter “O” ring (not shown) was brought into contact with the coverslip and four screws bring the sections together to provide a stable mount. Immersion fluid is added to the bath created by the contact of the top plate and coverslip of the specimen slide for long-term imaging.
Figure 3
Measurements of Airy and Gaussian light-sheet thickness across the large field of view (FOV)
(A–F) The parameter x is the length perpendicular to the light-sheet, and y is the position of measurement relative to the edge of the FOV, with experimental measured data and theoretical fits shown. The measured FWHM for the Airy light-sheet was (A) 7.8 μm ± 1.3 μm at y = 0 mm position, (B) 10.1 μm ± 1.5 μm at the central position (y = 1.5 mm), and (C) 11.1 μm ± 0.9 μm at y = 3 mm. The measured FWHM for the Gaussian light-sheet was (D) 36.8 μm ± 4.1 μm at y = 0 mm position, (E) 30.1 μm ± 5.3 μm at the central position (y = 1.5 mm), and (F) 31.3 μm ± 2.1 μm at y = 3 mm. The light-sheet profile was fitted with the Airy function (Siviloglou and Christodoulides, 2007). This measurement confirms that we generated a static Airy light-sheet with quasi-constant thickness across the entire Mesolens FOV that provides a 4-fold axial resolution improvement compared to a static Gaussian light-sheet.
Schematic diagram of both light-sheet illuminators for use with the MesolensIn both cases, the laser source is a Coherent Sapphire 488-10 CDRH laser emitting light at 488 nm with a maximum of 10 mW average power at the sample plane, with a single Gaussian mode of 0.8 mm diameter and M2 value ∼1.1.(A) The laser beam was expanded to a diameter of 12 mm (not shown) and a flip mirror (not shown) was used to select the light path used for imaging. The light-sheet coming from the front and represented with a cross is the Airy light-sheet, the light-sheet on the right is the Gaussian illuminator. The chromatic reflector (Chroma) is placed at a 30° angle of incidence to introduce the incoherent light source for widefield epi-fluorescence imaging, while the compensation plate (Chroma) is oriented at a −30° angle to compensate for the astigmatism in the fluorescence signal. The emission filter after the Mesolens is a custom-designed 100 mm diameter filter (Chroma), with 535 ± 20 nm collection. A thermoelectric Peltier cooled camera (VNP-29MC, Vieworks) with a chip-shifting mechanism was used for image acquisition.(B) Schematic diagram of the optical path created to produce the Gaussian light-sheet.(C) Schematic diagram of the optical path used to create the Airy light-sheet. The red arrow shows the direction of offset for the Powell lens to create the Airy beam.A computer-aided design drawing of the disassembled specimen holder and multipart chamber to support the long-term use of immersion fluid for light-sheet mesoscopy with the MesolensSpecimens were mounted in their clearing solutions in a quartz cuvette cut to provide an inner wall height of 3 mm. The cut cuvette was then placed into a bottom plate fabricated from black acetal resin plastic. A 3D printed spacer was used to surround the cuvette on three sides for stability, leaving the fourth side open to permit light-sheet illumination, as indicated by the blue arrow. A 60 mm diameter borosilicate coverslip was placed on top of the cuvette and spacer. A top plate comprised a 30 mm diameter “O” ring (not shown) was brought into contact with the coverslip and four screws bring the sections together to provide a stable mount. Immersion fluid is added to the bath created by the contact of the top plate and coverslip of the specimen slide for long-term imaging.Measurements of Airy and Gaussian light-sheet thickness across the large field of view (FOV)(A–F) The parameter x is the length perpendicular to the light-sheet, and y is the position of measurement relative to the edge of the FOV, with experimental measured data and theoretical fits shown. The measured FWHM for the Airy light-sheet was (A) 7.8 μm ± 1.3 μm at y = 0 mm position, (B) 10.1 μm ± 1.5 μm at the central position (y = 1.5 mm), and (C) 11.1 μm ± 0.9 μm at y = 3 mm. The measured FWHM for the Gaussian light-sheet was (D) 36.8 μm ± 4.1 μm at y = 0 mm position, (E) 30.1 μm ± 5.3 μm at the central position (y = 1.5 mm), and (F) 31.3 μm ± 2.1 μm at y = 3 mm. The light-sheet profile was fitted with the Airy function (Siviloglou and Christodoulides, 2007). This measurement confirms that we generated a static Airy light-sheet with quasi-constant thickness across the entire Mesolens FOV that provides a 4-fold axial resolution improvement compared to a static Gaussian light-sheet.The length y corresponds to the distance across the FOV from the edge where the beam enters the photometric area.First, we present the Gaussian results. The intensity profile along y allowed FWHM measurement of Gaussian light-sheet thickness as 36.8 μm ± 4.1 μm at y = 0 mm position, (e) 30.1 μm ± 5.3 μm at the central position (y = 1.5 mm), and (f) 31.3 μm ± 2.1 μm at y = 3 mm.Secondly, the Airy light-sheet intensity profile along y was measured by taking the FWHM of the first maximum only. The thickness of the Airy light-sheet was measured as 7.8 μm ± 1.3 μm at y = 0 mm position, 10.1 μm ± 1.5 μm at the central position (y = 1.5 mm), and 11.1 μm ± 0.9 μm at y = 3 mm. The Airy light-sheet profile was fitted with the Airy function as shown in Equation (3). This measurement confirmed we generated a static Airy light-sheet with quasi-constant thickness across the entire Mesolens FOV that provides an almost 4-fold axial resolution improvement compared to a static Gaussian light-sheet. The experimental measurements compared well with theory. From Equation (4), and with λ = 488 nm, n = 1.515, effective NA of the final lens element before the specimen = 0.043, and α = 0.7, we calculate a theoretical FOV of 3306 μm. From Equation (5) we expected the theoretical axial resolution from the Airy light-sheet to be 7.6 μm, which is close to our measured values of FWHM.
Comparison of point-scanning confocal, Gaussian light-sheet and Airy light-sheet imaging with the Mesolens
Figure 4 shows images of mouse pancreas cleared in a 1:2 mixture of benzyl alcohol:benzyl benzoate (BABB), with nuclei stained with propidium iodide (PI) with xy, xz, and yz cross-sectional views obtained with Gaussian and Airy light-sheets and confocal point-scanning illumination. In the confocal images, individual cell nuclei are clearly resolved in three dimensions throughout the tissue volume. However, with Gaussian light-sheet illumination cell nuclei are barely resolved in the xy plane and cannot be resolved at all in the yz or xz planes. The relatively thick illuminating beam directly accounts for the poor resolving power in the z-direction. We attribute the unexpected decrease in lateral resolution performance to the Gaussian light-sheet exciting fluorescence with a thicker slice than both confocal and Airy illumination, the result of which is analogous to the capture of out-of-focus fluorescence in widefield imaging.
Figure 4
A comparison of point-scanning confocal, Gaussian light-sheet, and Airy light-sheet imaging with the Mesolens using optically cleared mouse pancreas specimens prepared a fluorescent nuclear marker
(A–L) Confocal images are shown in (A) xy, (B) xz and (C) yz, with a region of interest (ROI) indicated with a white box. Gaussian images are shown in (D) xy, (E) xz, and (F) yz, with a second ROI indicated with a white box. Airy images are shown in (G) xy, (H) xz, and (I) yz, with a third ROI highlighted with a white box. The ROIs are shown with a digital zoom in (J) confocal, (K) Gaussian, and (L) Airy. All images have been converted to 8-bit depth and are shown with the false-colour “mpl-inferno” look-up table. The axial resolution of the images obtained with the Gaussian illuminator is markedly poorer than those obtained with the Airy illuminator, which performs similarly to the confocal method. The striping artifacts visible in (D) are directed from right to left while in (G) are from bottom to top: this is owing to the orientation of the sample with respect to the Gaussian and Airy light-sheets, respectively. Scale bar: (A–I): 500 μm, (J–L): 100 μm.
A comparison of point-scanning confocal, Gaussian light-sheet, and Airy light-sheet imaging with the Mesolens using optically cleared mouse pancreas specimens prepared a fluorescent nuclear marker(A–L) Confocal images are shown in (A) xy, (B) xz and (C) yz, with a region of interest (ROI) indicated with a white box. Gaussian images are shown in (D) xy, (E) xz, and (F) yz, with a second ROI indicated with a white box. Airy images are shown in (G) xy, (H) xz, and (I) yz, with a third ROI highlighted with a white box. The ROIs are shown with a digital zoom in (J) confocal, (K) Gaussian, and (L) Airy. All images have been converted to 8-bit depth and are shown with the false-colour “mpl-inferno” look-up table. The axial resolution of the images obtained with the Gaussian illuminator is markedly poorer than those obtained with the Airy illuminator, which performs similarly to the confocal method. The striping artifacts visible in (D) are directed from right to left while in (G) are from bottom to top: this is owing to the orientation of the sample with respect to the Gaussian and Airy light-sheets, respectively. Scale bar: (A–I): 500 μm, (J–L): 100 μm.The Airy light-sheet resolution performance is almost identical to that obtained with confocal point scanning. It was not necessary to deconvolve or otherwise process the data to remove the subsidiary maxima inherent in the Airy beam because the short focal depth of the Mesolens (McConnell et al., 2016) allows these signals to be excluded optically. The stripe patterns characteristic of light-sheet illumination resulting from tissue scatter (Mayer et al., 2018) were more pronounced with the Airy light-sheet because of the thinness of the illumination source.
Comparison of Gaussian and Airy light-sheet imaging with deconvolution
Figure 5 shows a maximum intensity projection image of 100 optical sections of mouse intestine imaged with both Gaussian and Airy light-sheets, with and without post-acquisition deconvolution. For visualization, the Contrast Limited Adaptive Histogram Equalization (CLAHE) algorithm (Zuiderfeld, 1994) was applied to (c) and (d). The resolution of both the raw and deconvolved images obtained with the Airy illuminator provides the highest spatial resolution, with the deconvolved image providing the best signal-to-background ratio. We note that the raw image obtained with the Airy illuminator provides substantially better spatial resolution than both the raw and deconvolved images acquired with Gaussian illumination. There is some improvement in the image quality when the Airy light-sheet image data are deconvolved; however, this adds additional time and requires more processing power. Three-dimensional animation of the imaging volume acquired with both Gaussian and Airy light-sheets are shown in Videos S1 and S2.
Figure 5
A comparison of raw and deconvolved 8-bit images of optically cleared mouse intestine prepared with a nuclear marker
(A) Raw and (C) deconvolved xy images obtained with the Gaussian illuminator.
(B) Raw and (D) deconvolved xy images obtained with the Airy illuminator.
(E) Normalised line intensity profile as indicated by the region of interest shown in (A) that was also applied to images (B–D). These are shown as (i) raw data from Gaussian beam illumination, (ii) raw data from Airy beam illuminator, (iii) deconvolved data from Gaussian beam illumination, and (iv) deconvolved data from Airy beam illumination. The resolution of both the raw and deconvolved images obtained with the Airy illuminator provides highest spatial resolution, with the deconvolved image providing the best signal-to-background ratio. The individual peaks in the raw Airy intensity profile correspond to individual nuclei. We note that the raw image obtained with the Airy illuminator provides better spatial resolution than the raw or deconvolved images acquired with Gaussian illumination. For visualization purpose, CLAHE was applied to (C) and (D). Scale bar = 500 μm.
A comparison of raw and deconvolved 8-bit images of optically cleared mouse intestine prepared with a nuclear marker(A) Raw and (C) deconvolved xy images obtained with the Gaussian illuminator.(B) Raw and (D) deconvolved xy images obtained with the Airy illuminator.(E) Normalised line intensity profile as indicated by the region of interest shown in (A) that was also applied to images (B–D). These are shown as (i) raw data from Gaussian beam illumination, (ii) raw data from Airy beam illuminator, (iii) deconvolved data from Gaussian beam illumination, and (iv) deconvolved data from Airy beam illumination. The resolution of both the raw and deconvolved images obtained with the Airy illuminator provides highest spatial resolution, with the deconvolved image providing the best signal-to-background ratio. The individual peaks in the raw Airy intensity profile correspond to individual nuclei. We note that the raw image obtained with the Airy illuminator provides better spatial resolution than the raw or deconvolved images acquired with Gaussian illumination. For visualization purpose, CLAHE was applied to (C) and (D). Scale bar = 500 μm.
Video S1. Three-dimensional visualisation of the optically cleared mouse intestine related to Figure 5a
The movie was produced using the proprietary visualisation and analysis software Imaris with the ‘Fire’ look-up table. The volume was acquired with the Gaussian light-sheet illuminator and the visualised sub-volume is composed of 100 optical sections.
Video S2. Three-dimensional visualisation of the optically cleared mouse intestine related to Figure 5b
The movie was produced using the proprietary visualisation and analysis software Imaris with the ‘Fire’ look-up table. The volume is acquired with the Airy light-sheet illuminator and the visualised sub-volume is composed of 100 optical sections.
Airy light-sheet Mesolens imaging of cleared mouse brain is shown in Figure 6 to demonstrate the application of the method with the same specimen mounting protocol but with antibody labeling and an alternative clearing strategy (Murakami et al., 2018). We imaged the forebrain region, including the dorsal and ventral 3rd ventricle, the hippocampus, the corpus callosum, and the external capsule of the caudoputamen as previous work have shown astrocytic and potential ependymal reactivity in these brain structures following infection (Amin et al., 2009; Cespuglio et al., 2019; Quintana et al., 2022; Tesoriero et al., 2018). After deconvolution, dendrites of astrocytes and ependymal cells are clearly resolved based on the reactivity of the putative marker glial fibrillary acidic protein (GFAP), including in the corpus callosum (Figure 6C), circumventricular organs (Figures 6B–6E), and thalamus (Figure 6D), thus confirming that sub-cellular resolution is retained over this large FOV.
Figure 6
Airy light-sheet Mesolens imaging of cleared coronal section of the mouse forebrain
(A) A maximum intensity projection of 100 Airy light-sheet images of mouse brain stained with anti-GFAP mouse monoclonal antibody coupled with Alexa Fluor 488 and cleared using CUBIC-R, shown after deconvolution. A z-step size of 1.5 μm was used; therefore, spanning a total specimen thickness of 150 μm. Digitally zoomed regions (B) and (E) show both sub-cellular resolution of astrocytes and ependymal cells, and (C) and (D) show GFAP reactive astrocytes only. The regions shown in the image correspond to the commissure of the cingulate cortex (B), external capsule and caudoputamen (C and D), and 3rd ventricle (E). The look-up table for (A–E) is mpl-inferno, and the images were converted to 8-bit values (0-255). Scale bar: (A) 1000 μm, (B–E) 50 μm. The vertical line is not indicative of stitching: this is an artifact of the imaging sensor.
Airy light-sheet Mesolens imaging of cleared coronal section of the mouse forebrain(A) A maximum intensity projection of 100 Airy light-sheet images of mouse brain stained with anti-GFAP mouse monoclonal antibody coupled with Alexa Fluor 488 and cleared using CUBIC-R, shown after deconvolution. A z-step size of 1.5 μm was used; therefore, spanning a total specimen thickness of 150 μm. Digitally zoomed regions (B) and (E) show both sub-cellular resolution of astrocytes and ependymal cells, and (C) and (D) show GFAP reactive astrocytes only. The regions shown in the image correspond to the commissure of the cingulate cortex (B), external capsule and caudoputamen (C and D), and 3rd ventricle (E). The look-up table for (A–E) is mpl-inferno, and the images were converted to 8-bit values (0-255). Scale bar: (A) 1000 μm, (B–E) 50 μm. The vertical line is not indicative of stitching: this is an artifact of the imaging sensor.
Discussion
The intended vast increase in speed relative to confocal point-scanning resulting from the combination of light-sheet illumination and Mesolens imaging with a camera has been demonstrated in this work, without a compromise in spatial resolution over the large FOV. For example, the 117 h mentioned above as needed for complete imaging of a 4.4 mm x 3 mm x 3 mm mouse embryo can be reduced to under 8.5 h.We show that the Mesolens FOV can be filled either with a Gaussian light-sheet or a thinner Airy light-sheet which provides better axial resolution. Either approach can be followed without any need for moving parts or adaptive optical elements and use only off-the-shelf optics. We regard the Gaussian and Airy light-sheets as having different advantages. Although the Airy light-sheet outperforms the Gaussian illuminator in optical sectioning, the Gaussian light-sheet illuminator is valuable for fast scanning, particularly when fluorescent objects are sparse and punctate or filamentous so out-of-focus glare is minimal.In a triangular diagram, in which the apices represent the maxima in information content, spatial resolution, and time resolution, the Airy light-sheet and Mesolens combination plainly scores highly along the information content and spatial resolution axes. We note that the resolution in the large images was far higher than can be reproduced here because of limitations of space, hence why only example regions of interest (ROIs) are shown as indicators of resolution performance in Figure 4, Figure 5, Figure 6. However, the time resolution is limited by the framing speed of available cameras with pixel numbers greater than 150 Megapixels. With a static light-sheet, the high time resolution needed for recording fast electrophysiological events such as action potentials lasting for less than 1 ms may, eventually, be achieved at the mesoscale if the total pixel number of high-speed cameras can be improved. Where low signal levels are encountered, the high optical throughput of the Mesolens may prove advantageous (McConnell et al., 2016).We have shown that our method is compatible with a range of biological specimens that have been prepared with different fluorescent dyes or which express photoproteins, and that have been optically cleared using various protocols. Although in this work we have used a laser with a fixed wavelength of 488 nm for illumination, we anticipate that our system could be upgraded to include additional laser wavelengths using achromatic optical elements for multiplexed mesoscale light-sheet imaging.We also note that our method is compatible with different specimen mounting media, with no measurable change in resolution performance. In addition, the use of commercially available quartz cuvettes which can be readily cut to size with a diamond saw allows the creation of specimen holders optimized for specimens of widely different sizes.In this study, no attempt was made to rotate the specimen or the light-sheet, so the well-known light-sheet streaking artifact was sometimes seen, as in Figures 4G and 5A–5D. Indeed, in both the mouse pancreas and mouse intestine datasets the streaking artifacts were more prominent in the Gaussian light-sheet images. This may be a consequence of the thicker Gaussian light-sheet profile experiencing greater scattering effects from the tissue, or that the self-healing properties of the non-diffractive Airy beam are minimizing streaking artifacts deeper into the tissue. It may be possible to reduce the streaking artifacts by scanning the beam in the plane of illumination very quickly relative to the camera exposure time (de Medeiros et al., 2015).
Limitations of the study
The greatest limitation of this approach is the slow speed of the high-resolution sensor-shifting camera. New CMOS devices with the same resolution and FOV are becoming commercially available which are appreciably faster. A perceived drawback of the Airy light-sheet illuminator itself might be the small alignment tolerance, but in practice, only one micrometer screw needs routine adjustment and sufficient accuracy can be achieved by simple manual control. A more exacting limitation of our new technology, and in light-sheet mesoscopy, is the need for the specimen to be optically transparent to depths of a few millimeters. There are many optical clearing techniques that can render specimens of this scale transparent, but all require the tissue to be fixed. At present, therefore, light-sheet optical mesoscopy cannot easily be applied to most living tissues. It may be possible to overcome this somewhat by using multiphoton excitation in combination with the new optical systems we have reported, but the peak power of the laser sources required would exceed 20 W and this may prove damaging to living tissue, even across a large FOV. A final limitation is that the working distance of the Mesolens is 3 mm, and this restricts the thickness of specimens to the same distance.
STAR★Methods
Key resources table
Resource availability
Lead contact
Further information and requests for resources and should be directed to and will be fulfilled by the lead contact, Eliana Battistella (eliana.battistella@strath.ac.uk).
Materials availability
This study did not generate new unique reagents.
Experimental model and subject details
Trypanosoma brucei brucei strains
Pleomorphic bloodstream form parasites AnTat 1.1E were used for infecting mice.
Mouse strains
Nine-week-old wildtype C57BL/6 male mice were used. Mice were housed under 12:12 light-dark cycles, at constant room temperature and humidity (19–23C–45–65% humidity). Eight-week-old wildtype C57BL/6J female mice were used. Animals were infected with T. brucei brucei Antat 1.1 and naive controls. Animals were only used for the purpose of these experiments. Mice were housed under 12:12 light-dark cycles, at constant room temperature and humidity (21-22C - 50–60% humidity). All experimental procedures using animals were conducted in strict accordance with the United Kingdom Animals (Scientific Procedures) Act, 1986 and approved by the Home Office (UK).
Method details
A schematic diagram of both light-sheet illuminator configurations is shown in Figure 1. We illuminated the specimens along the shorter dimension of the camera FOV (3 mm) using the Airy light-sheet and along the longer dimension (4.4 mm) using the Gaussian light-sheet.The static Gaussian light-sheet and the static Airy light-sheet were generated through two distinct optical paths and used a flip mirror to direct the required light-sheet. The light source used for both illuminators was a Coherent Sapphire 488-10 CDRH laser emitting light at 488 nm with a maximum of 4 mW average power at the sample plane, with a single Gaussian mode of 0.8 mm diameter and M2 value ∼1.1.
Gaussian light-sheet illuminator
The 0.8 mm diameter laser beam was expanded to a diameter of 8 mm using a home-made telescope (f = 25 mm LA1951-A Thorlabs, f = 250 mm LA1461-A Thorlabs) and then directed through a cylindrical lens (f = 300 mm LJ1558RM-A, Thorlabs) to form a static Gaussian light-sheet at the specimen plane of the Mesolens approximately 300 mm from the cylindrical element. The cylindrical lens was mounted on a vertical breadboard by means of an optical cage system to precisely position it at the height of the Mesolens focal plane. To generate a 3 mm-long light-sheet (which corresponds to the FOV of the Mesolens) and considering a 488 nm wavelength, from Equations (1) and (2) the theoretically calculated minimum beam thickness is 30 μm, with a thickness of 42 μm at the edges of the FOV.
Airy light-sheet illuminator
The laser beam was expanded to a diameter of 1.2 mm using two home-made telescopes (f = 25 mm LA1951-A Thorlabs, f = 150 mm LA1433-A Thorlabs; and f = 500 mm LA1908-A Thorlabs, f = 125 mm LA1986-A Thorlabs) and guided through the following optical elements. Firstly, the beam propagated through two identical aspheric lenses (AL, f = 20 mm. AL2520-A, Thorlabs) separated by 80 mm. A Powell lens (PL, fan angle = 10°, Laserline Optics Canada) was positioned at 40 mm distance from the first aspheric lens. Its axial position was regulated via an XY translation mount (CXY1QA, Thorlabs), translation being achieved by a screw of pitch 250 μm/revolution. This allowed the axial position of the Powell lens to be readily shifted by a few hundred nanometres from the main optical axis to generate the Airy light-sheet. Two cylindrical lenses were used to collimate and refocus the beam into the specimen. The first cylindrical lens (CL, f = 75 mm LJ1703RM-A, Thorlabs) was positioned 95 mm from the second aspheric lens. The second identical cylindrical lens was positioned at 75 mm from the previous cylindrical lens. The power direction of both cylindrical lenses was oriented perpendicularly to the beam propagation direction. The aspheric optical elements were mounted on a vertical breadboard using an optical cage system, which was precisely positioned in front of the Mesolens specimen stage (Prior V31XYZE) by means of four micrometre-adjustable translation stages, two for the lateral position and two for the vertical position. The second cylindrical lens was a distance of approximately 110 mm from the edge of the FOV closest to the laser source.
Practical adjustments during use
The initial setup of the Airy light-sheet illuminator was achieved by standard use of pinholes. In practical use, the only adjustment requiring frequent checking, i.e. before each work session, is the translation of the Powell lens in the direction shown by the red arrow in Figure 1C. The thinness of the light-sheet was monitored by observing the beam directly while turning the fine pitch translation control screw.
Measurement of light-sheet thickness
We used a reflective surface to directly image the light-sheet illumination and measure its thickness across the entire field of view by evaluating the full width at half maximum (FWHM) of a section of the reflected image (Voie et al., 1993; McConnell et al., 2011). The reflective element was a type 1.5 borosilicate microscope coverslip (631-0150P, VWR) coated with a 7 nm-thin layer of aluminum. An evaporator (Edwards, coating system E306A) was used for the aluminum coating. The aluminium-coated coverslip was placed in a quartz cuvette (CV10Q3500F, Thorlabs). The aluminised coverslip was placed at an angle of 45° with respect to the detection objective, facing the incoming light-sheet optical path. The coverslip was held in place by means of a 1.1% agarose (A9414-10G, Sigma-Aldrich) plug and the cuvette was filled with immersion oil type LDF (16,241, Cargille). For this measurement, the laser was set at minimum power to avoid any camera damage and the fluorescence emission filters were removed from the detection path to directly image the reflected laser light. In addition, we used this target-sample consisting of the aluminium-coated coverslip to precisely set the Powell lens position, the cylindrical lenses power axis, and the translation stage positions. This was accomplished through an iterative process that consisted in observing the light-sheet reflection through the Mesolens optics and minimising its profile thickness by either translating the Powell lens or the micrometre-adjustable stages. For all measurements, the correction collars of the Mesolens were set to be compatible with oil immersion, and type LDF oil was used throughout.
Specimen preparation
Specimens of mouse tissue were prepared for imaging with the light-sheet illuminators.
Mouse pancreas
A 9-week-old male C57BL/6 mouse was euthanized by cervical dislocation. The mouse was immobilised in a supine position on a dissection block and sprayed with 70% ethanol. A vertical incision was made along the abdominal midline, the skin was retracted and the pancreas was carefully dissected with forceps. The tissue was washed in a 50 mL Falcon tube of cold 1X PBS (10010023, ThermoFisher Scientific, Renfrew, UK) using a rotating mixer (3x, 30 min). The tissue was then fixed in 4% PFA, rotating at 4°C overnight. The specimen was then washed in 1X PBS, rotating at room temperature (RT) for 30 min, 3x.The pancreas was dehydrated, stepwise in methanol (MeOH) (322415-1L, Sigma-Aldrich, Gillingham, UK)/H2O (20%, 40%, 60%, 80%, 100%, 100%) rotating at RT (30 min per stage). The specimen was stored in 100% MeOH overnight. The specimen was then rehydrated, stepwise in MeOH/H2O (80%, 60%, 40%, 20%) rotating at RT (30 min per stage). The specimen was then washed in 1X PBS, rotating at RT for 30 min, 3x. The specimen was treated with 100 μM RNAse (EN0531, ThermoFisher Scientific, Renfrew, UK) in 1X PBS for 3 h, rotating at RT. Without washing, the specimen was transferred to an aluminum foil-covered Falcon tube (to reduce bleaching by ambient light) and stained with 100 μM PI (P4864-10 mL, ThermoFisher Scientific, Renfrew, UK), rotating at RT overnight. Next, the fluorescently-stained specimen was washed in 1X PBS (3x, 30 min) in the dark.The specimen was then dehydrated stepwise in MeOH (as above, but in a foil-covered Falcon tube). To begin clearing, the pancreas was transferred to a foil-covered glass vial containing a 1:1 solution of MeOH (100%) and BABB (1:2 (v/v) benzyl alcohol (402834-500 mL, Sigma-Aldrich, Gillingham, UK) to benzyl benzoate (B6630-1L, Sigma-Aldrich, Gillingham, UK)), gently rocking overnight. Finally, tissue was placed in 100% BABB, gently rocking until cleared.
Mouse intestine
Specimens of mouse intestine were prepared, using the same mouse type and euthanasia method as for the mouse pancreas. After excision, tissue was cut into 6 mm long sections and rinsed in PBS (3x, 5 min). The specimens were placed in a Falcon tube of fixative (4% PFA) and were placed on a tube rotator overnight at RT. Following fixation, the specimens were washed three times in 1X PBS (10010023, ThermoFisher Scientific, Renfrew, UK), 30 min each. To remove absorbing pigment from hemoglobin, the specimens were placed in 35% H2O2 (349887-500 mL, Sigma-Aldrich, Gillingham, UK) for 18 h, again with gentle rotation, then washed in 1X PBS (3x, 5 min). Next, the bleached specimens were treated with 100 μM RNAse in 1X PBS for 1.5 h at RT with gentle agitation. Without washing, the specimens were then added to 10 μM PI and gently rocked for 4 h. After the dye loading and for subsequent steps, the specimens were covered with aluminum foil.The fluorescently stained specimens were washed in 1X PBS (3x, 5 min) and then dehydrated through a methanol series (50% MeOH, 75% MeOH, 100% anhydrous MeOH, 100% anhydrous MeOH, each step 1 h) with gentle agitation. BABB was introduced by placing the specimens in a 1:1 by volume mixture of anhydrous MeOH and BABB, the latter being a 1:2 mixture of benzyl alcohol to benzyl benzoate and rocked for 1 h. The specimens were removed from the MeOH/BABB mix and placed in 100% BABB and gently rocked for at least 24 h before imaging.
Mouse brain
Eight-week-old female C57BL/6J mice (Jackson Labs, UK) were inoculated by intraperitoneal injection with 104 T b. brucei AnTat 1.1E. Parasitaemia was monitored every two days from the superficial tail vein, examined using phase microscopy, and scored using the rapid “matching” method (Herbert and Lumsden, 1976). At 21 days post-infection, mice were culled by cervical dislocation under terminal anesthesia with isoflurane and decapitated to remove intact brain. After overnight fixation in 10% neutral buffered formalin, samples were cleared using CUBIC (TCI Chemicals) [Tainaka et al., 2018]. Briefly, specimens were first delipidated using CUBIC-L for a period of ∼6–7 days at 37°C. After clearing was completed, samples were washed three times in 1X PBS (including an overnight wash at 4°C). Blocking was performed for 2 h at room temperature with 1X PBS supplemented with 5% Fetal Bovine Serum (Invitrogen) and 0.25% Tween-20 (Invitrogen) and stained with anti-GFAP mouse monoclonal antibody coupled to Alexa Fluor 488 (Clone GA5, 1/100, Thermo) for 3 days at 37°C and with gentle agitation. Specimens were then washed three times with 1X PBS to remove excess antibody and incubated in CUBIC-R+ (TCI chemicals) for 3 days for refractive index matching.
Specimen mounting
A schematic diagram of the custom-designed specimen mount for light-sheet imaging with the Mesolens is shown in Figure 2. For light-sheet imaging, the cleared specimens were mounted in their respective clearing agents in a quartz cuvette (CV10Q3500F, Thorlabs) that was cut with a Beuhler diamond saw to produce walls of 3 mm height. This small chamber was then placed within a custom-designed mount fabricated from black acetal resin plastic. This mount comprised an optical window compatible with brightfield transmission imaging (not performed in this study) that when embedded in the resin plastic formed a cylindrical well of 10 mm diameter and 3 mm height. A 3D printed square-C shape was fabricated using PLA plastic to provide a spacer to give access to the incoming light-sheet illuminator. A 60 mm diameter borosilicate coverslip (0107999098, Marienfeld) was placed on top of the cuvette. A top plate, fabricated in the same black acetal resin plastic as the bottom mount, was used to sandwich the specimen held in the cuvette and coverslip, and this plate was held in place with corner screws. The design of the top plate was such that a 30 mm diameter well was formed between the acetal plastic and the coverslip, allowing oil immersion with type LDF immersion oil (16241, Cargille). The Mesolens objective is equipped with correction collars to minimise spherical aberrations and in this work, the collars were set to match the refractive index of the clearing solutions used.For comparisons with confocal imaging, specimens were placed in the black acetal resin plastic plate and surrounded by the clearing media. A 60 mm diameter borosilicate coverslip (0107999098, Marienfeld) was placed on top of the specimen. This was held in place using the same black acetal resin plastic plate and screws used for light-sheet imaging, and type LDF immersion oil (16241, Cargille) was added to the bath created between the coverslip and the top plate.The intestine specimen was oriented with the longitudinal section facing the incident light-sheet, while the pancreas specimen was not mounted with any specific directionality. The axial step size between images in z for the intestine specimen was 1.5 μm for the Airy light-sheet illumination and 3.0 μm for the Gaussian light-sheet illumination. For the pancreas specimen, axial step size was 2.0 μm for the Airy light-sheet illumination, 3.0 μm for the Gaussian light-sheet illumination, and 4.0 μm for confocal laser scanning.The brain specimen was oriented so that an optical section corresponded to a coronal section, with the retrosplenial area being the first part of the specimen illuminated by the light-sheet. The region considered in the first instance was visibly identified by the ventral part of the retrosplenial area and by part of the corpus callosum. The specimen was then translated along the y direction to image both the thalamus and hypothalamus in the infected brain specimen. The axial displacement z was chosen to satisfy the Nyquist sampling criterion (Shannon, 1949) of 1.5 μm for the Airy light-sheet illumination, and 3.0 μm for the Gaussian light-sheet illumination.Before mounting the sample, each cuvette was carefully cleaned by removing large impurities with isopropanol and wiping with lens tissue paper. The cuvette was subsequently immersed overnight in a bath of 2% detergent (Hellmanex III, Hellma) in distilled water, and dried with lens tissue paper before use.
Image acquisition parameters
Light-sheet image datasets were acquired using the custom ‘MesoCam’ software, which controlled a sensor-shifting camera and acquired images of 12-bit depth in the TIFF format. Nyquist sampling for the Mesolens FOV in widefield (and therefore light-sheet) modality corresponds to a 224 nm pixel size, for a total of 260MP. There are no commercially available cameras with said resolution, but the Mesolens camera detector (VNP-29MC, Vieworks) is equipped with a chip-shifting mechanism that by shifting the image by 1/3 of a pixel in a 3 × 3 square allows Nyquist sampling (Schniete et al., 2018; Shaw et al., 2021).For light-sheet imaging of mouse pancreas, images were acquired with 100 ms exposure time for each of the 9 sensor positions and with a gain setting equal to 1. The entire volume was composed of 615 slices, which was acquired in approximately 5 h. For light-sheet imaging of mouse intestine, all images were acquired with 200 ms exposure with a gain setting equal to 1. The entire volume was composed of 1013 slices and was acquired in approximately 8 h for Airy illumination, and 4 h for Gaussian illumination. The increase in time for Airy illumination was due to the thinner light-sheet requiring approximately double the number of images for the same specimen thickness.For imaging of mouse brain, all images were acquired with 1000 ms exposure time and with a gain setting equal to 50. The acquisition of a single optical section took 25 s as the overall acquisition speed was limited by the time employed by the MesoCam software to transfer and save the acquired image from the camera to the PC. The acquisition volume for the Airy light-sheet modality was composed of 933 slices, at the full 3.0 × 4.4 mm FOV. The entire volume was acquired in approximately 8 h.Confocal Mesolens imaging of mouse pancreas specimens were performed for comparison of the image quality with the light-sheet illuminators. Details of the Mesolens are already reported (McConnell et al., 2016), so only the imaging parameters used in these experiments are described here. For fluorescence excitation of PI, laser powers of no more than 5 mW (Laserbank, Cairn Research, Faversham, UK) at a wavelength of 488 nm was used, with less than 200 μW of total laser power incident on the specimen during scanned imaging. Fluorescence from the PI stain was spectrally separated from the 488 nm excitation using a 550 nm dichroic filter (DMLP550R, Thorlabs, Newton, USA) with a 600 nm long-wave pass filter (FEL0600, Thorlabs, Newton, USA), and was detected using a photomultiplier tube (P30-09, Senstech, Egham, UK). A galvo mirror scan speed of 40 Hz was used, and each 12,000 pixel × 12,000 pixel image took around 4.8 min to acquire. The full 6 × 6 mm FOV of the Mesolens in confocal point-scanning mode was scanned. Images were saved as 11-bit OME TIFF files.
Image post-processing
The images were visualised in two dimensions and processed with Fiji (Schindelin et al., 2012) using virtual stacks. The outliers due to camera noise were removed using the batch processor for virtual stacks with the ‘RemoveOutliers’ function with a radius equal to 3.The maximum intensity projection for each dataset was computed with the ‘Z Project’ function. The three-dimensional data handling and visualisation was possible thanks to a Python script using the Dask library and the visualisation framework napari. Imaris (Oxford Instruments, 2021) was used for the three-dimensional visualisation and rendering.Light-sheet image datasets were deconvolved using proprietary software (Huygens, Scientific Volume Imaging). The deconvolution operation was performed using the classical maximum likelihood estimation (CMLE) method on a sub-volume comprising 100 optical slices, which corresponds to an approximately 50 GB file size. This sub-volume file size was chosen since our computational power was limited by the available server RAM size. Furthermore, the process was time-consuming, exceeding 12 h computational time per dataset. For the Airy light-sheet, deconvolution parameters were chosen to approximate a Gaussian light-sheet with a 12 μm waist, directed from the bottom of the image, and with no focus offset. For the Gaussian light-sheet, the waist was set at 25 μm. The theoretical point spread function of the system was generated automatically by the software, using the previously described Gaussian parameters. The other parameters were set to describe the excitation and peak emission wavelengths (488 and 520 nm), the sampling intervals (x and y: 224.4 nm, z: 1500 nm for the Airy light-sheet, z: 3000 nm for the Gaussian light-sheet), and the optical parameters (numerical aperture 0.46, LDF oil immersion refractive index: 1.515).Image post-processing was performed subsequently to the image acquisition and therefore did not affect the acquisition time. The removal of outliers took approximately 3 h, and the maximum intensity projection was performed in a few minutes.The images presented in this manuscript were converted to 8-bit for the purpose of visualisation.
Authors: Zhengyi Yang; Martynas Prokopas; Jonathan Nylk; Clara Coll-Lladó; Frank J Gunn-Moore; David E K Ferrier; Tom Vettenburg; Kishan Dholakia Journal: Biomed Opt Express Date: 2014-09-05 Impact factor: 3.732
Authors: Matthew Jemielita; Michael J Taormina; April Delaurier; Charles B Kimmel; Raghuveer Parthasarathy Journal: J Biophotonics Date: 2012-12-14 Impact factor: 3.207
Authors: Florian O Fahrbach; Vasily Gurchenkov; Kevin Alessandri; Pierre Nassoy; Alexander Rohrbach Journal: Opt Express Date: 2013-06-03 Impact factor: 3.894
Authors: Raju Tomer; Matthew Lovett-Barron; Isaac Kauvar; Aaron Andalman; Vanessa M Burns; Sethuraman Sankaran; Logan Grosenick; Michael Broxton; Samuel Yang; Karl Deisseroth Journal: Cell Date: 2015-12-17 Impact factor: 41.582
Authors: Fabian F Voigt; Daniel Kirschenbaum; Evgenia Platonova; Stéphane Pagès; Robert A A Campbell; Rahel Kastli; Martina Schaettin; Ladan Egolf; Alexander van der Bourg; Philipp Bethge; Karen Haenraets; Noémie Frézel; Thomas Topilko; Paola Perin; Daniel Hillier; Sven Hildebrand; Anna Schueth; Alard Roebroeck; Botond Roska; Esther T Stoeckli; Roberto Pizzala; Nicolas Renier; Hanns Ulrich Zeilhofer; Theofanis Karayannis; Urs Ziegler; Laura Batti; Anthony Holtmaat; Christian Lüscher; Adriano Aguzzi; Fritjof Helmchen Journal: Nat Methods Date: 2019-09-16 Impact factor: 28.547