Confocal multispot microscope for fast and deep imaging in semicleared tissues

Although perfectly transparent specimens are imaged faster with light-sheet microscopy, less transparent samples are often imaged with two-photon microscopy leveraging its robustness to scattering; however, at the price of increased acquisition times. Clearing methods that are capable of rendering strongly scattering samples such as brain tissue perfectly transparent specimens are often complex, costly, and time intensive, even though for many applications a slightly lower level of tissue transparency is sufficient and easily achieved with simpler and faster methods. Here, we present a microscope type that has been geared toward the imaging of semicleared tissue by combining multispot two-photon excitation with rolling shutter wide-field detection to image deep and fast inside semicleared mouse brain. We present a theoretical and experimental evaluation of the point spread function and contrast as a function of shutter size. Finally, we demonstrate microscope performance in fixed brain slices by imaging dendritic spines up to 400-μm deep. © The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI. [DOI: 10.1117/1.

Brain function and structure are intrinsically correlated and therefore best elucidated on an organ-wide scale in threedimensional (3-D) at cellular resolution.Observing such large volumes, as, for example, intact mammalian brains, requires imaging techniques, which can penetrate deep into the tissue and ideally also afford fast image acquisition to gather statistically relevant datasets.Both issues, such as depth and speed, have been addressed by recent advances in tissue clearing methodologies, [1][2][3] which are capable of rendering the brain transparent and thus provide optical access.The reduced scattering of cleared samples not only allows for deeper tissue imaging but also enables wide-field detection schemes that suffer from scatter-induced optical cross talk.Light-sheet microscopy (LSM), 4,5 a method taking advantage of wide-field detection to allow for fast acquisition speeds at cellular resolution, requires the sample to be perfectly transparent and consequently, a variety of samples that can be studied with LSM strongly depends on the effectiveness of the tissue clearing method.The most successful methodologies to render a brain transparent are also the most complex, costly, and time-consuming ones, whereas simpler and cheaper methods are now emerging, which render the brain not perfectly clear but instead "semitransparent." Two-photon microscopy (TPM), on the other hand, is a wellestablished technique 6 for imaging even in highly scattering media 7 owing to decreased scattering at long wavelengths and the signal's quadratic dependence on the excitation intensity.Being a point-scanning technique, however, TPM suffers from slow scan rates.To increase imaging speed, several strategies have been explored, for example, multispot multiphoton microscopy (MMM) using either single spot detectors arrays 8 or wide-field sensors. 9An increase in the speed gained by multiplexing the excitation; however, comes with a penalty in optical cross talk, which prohibits the application of MMM deep within scattering samples.][12] Here, we present a mixed microscope, which aims to capitalize on the advantages of MMM while preserving the optical sectioning capabilities of confocal systems.We combined fast and robust multispot two-photon fluorescence excitation with wide-field confocal line detection to reject cross talk.Our method benefits from the increased penetration depth of twophoton excitation in scattering media and the increased imaging rates of multiplexed excitation while exploiting the rolling shutter of a fast sCMOS camera.We present a full theoretical and experimental evaluation of the microscope performance and further validate our method by imaging fine neuronal structures 400-μm deep inside partially cleared mouse brain.
The microscope was built using a commercial inverted microscope (DMLFSA, Leica) with few modifications to the excitation and detection path [Fig.1(a)].Two-photon excitation was achieved with a high-power-pulsed Ti:Sa laser (Chameleon Ultra II, Coherent, 150-fs pulse width, 80-MHz repetition rate).The power at the back aperture of the microscope objective was kept below 740 mW.A line of 15 equidistant laser spots was generated at the focal plane of the microscope objective using a diffractive optic element (DOE, MS-454-920-Y-A, Holo/or).The DOE was reimaged onto the first of a pair of galvanometric mirrors (G1222HS-C355, Galvoline) using a 4f-system of unity magnification, and the second galvo was placed into a common conjugate plane using two spherical mirrors again of unity magnification.A second telescope in a 4f configuration reimaged both galvos and the DOE through the entry port of the inverted microscope onto the back aperture of the microscope objective (Zeiss, LD LCI Plan-Apochromat 25×∕0.8Imm Corr).
In the sample, the line of excitation spots was scanned rapidly along its axis to generate a homogeneous excitation line, which in turn was scanned perpendicularly to obtain a sheet of excitation of 90 μm × 96 μm.This field of view (FOV) resulted from the full DOE diffraction angle (5.27 deg) and total effective magnification.Fluorescence was collected in epidetection and separated from the excitation light by a longpass dichroic mirror (FF665-Di02-25, Semrock) and bandpass filters (FF01-520/35-25 and FF01-680/SP-25, Semrock).A tube lens (450 mm) formed the primary image onto the chip of a fast sCMOS camera (PCO.edge5.5, PCO) such that 1 pixel corresponded to 96 nm in the object plane.By sweeping the excitation line synchronously to the line readout of the camera [Fig.1(b)], a virtual confocal slit, adjustable in height and speed, partially rejected out-of-focus fluorescence and reduced cross talk induced by tissue scattering.
3-D stacks were acquired by scanning the microscope objective along the optical axis.A custom-written LabWindows code [National Instrument (NI)] in conjunction with a data acquisition card (PCI-6221, NI) controlled and synchronized the microscope.To demonstrate illumination homogeneity, as well as shutter synchronization, the same pollen grain was imaged in three different configurations: conventional single beam scanning and multispot scanning with and without rolling shutter [Fig.1(c)].The slight inhomogeneity in fluorescent signal, from left to right, that can be observed in the middle and right image, is due to an inhomogeneity in diffraction efficiency of our DOE that we estimated around 0.8 between the first and last diffraction orders.
To evaluate the optical performances of the microscope, we model the image formation process similarly to what has been done for light-sheet microscopy with confocal slit detection. 13e define x as the axis parallel to slit direction, and z as the optical axis of the microscope.We neglect all aberrations.We consider the intensity at ðx; y; zÞ of light emitted from a point object at (0, 0, 0) illuminated by a single excitation spot centered at ðx 0 ; y 0 ; z 0 Þ.Assuming that the laser focus and camera are perfectly conjugated and synchronized, the camera sensor lies at z ¼ z 0 , and the slit aperture is centered at y 0 .At a specific position of the scanning beam, since we are imaging on a pixelated array detector, the partial point spread function (PSF) of the microscope is simply the product of the PSF of a wide-field imaging system, denoted as H WF , the laser intensity I laser , and the aperture function of the slit AF E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 3 2 6 ; 7 5 2 d 2 H TOT ðx; y; z; x 0 ; y 0 ; z 0 ¼ zÞ ¼ H WF ðx; y; zÞ • I 2 laser ðx 0 ; y 0 ; zÞ • AFðy; y 0 Þdx 0 dy 0 ; (1) with E Q -T A R G E T ; t e m p : i n t r a l i n k -; ; 3 2 6 ; 7 0 0 0 elsewhere ; where J 0 is the zero order of the Bessel function, NA is the numerical aperture, k 0 is the wave vector modulus in vacuum, n is the refractive index of the medium between the object and the objective lens, ω 0 and ωðzÞ are the beam radii at the positions 0 and z, respectively, and s is the size of the slit.ωðzÞ can be written as ωðzÞ , where z R is the Rayleigh range.H WF is given by the Airy formula for Fraunhofer diffraction at a circular aperture 14 and I laser as the intensity of a focused Gaussian beam since we are slightly underfilling the objective back aperture.
As the laser is scanned over the FOV, H TOT corresponds to the integral sum of the intensity contributions d 2 H TOT of each laser position within the slit.Assuming that the camera extent along the x-axis is much larger than the slit width and the lateral extent of the partial PSF contributions, we have E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 3 2 6 ; 4 2 3 H TOT was numerically computed with Mathematica (Wolfram Research), using: λ 0 ¼ 510 nm, n ¼ 1.45, NA ¼ 0.8, and ω 0 ¼ 580 nm.The refractive index, corresponding to a solution of 63% of TDE in PBS, was measured with a refractometer.ω 0 was extracted by the corresponding beam waist measured at the back focal plane of the microscope objective with the knife-edge method.The full width at half maximum (FWHM) was calculated as a function of slit size from the PSF profile H TOT for the three axes, respectively.The resolution of the microscope was evaluated as a function of slit size both theoretically and experimentally in terms of its PSF and their FWHM and is summarized in Fig. 2(a).The PSF was measured experimentally by imaging fluorescent microspheres of 0.1-μm diameter in z-stacks (200-nm step size) for five different slits (3, 6, 12, 30, and 60 pixels).The beads (505-mm to 515-nm emission wavelength, Invitrogen) were diluted in a 63% 2,2′-thiodiethanol (TDE)/distilled water gel containing 4% of agarose.Twenty PSFs were manually extracted and fitted with Gaussian curves to obtain the FWHM along the 3-D [Fig.1(a)].The FWHMs were 368 AE 3 nm in x, 366 AE 3 nm in y, and 2.08 AE 0.02 μm in z for a slit size of 60 pixels [n ¼ 20, error is standard deviation (SD)].As the slit size decreased, the FWHMs in y and z decreased and reached 283 AE 9 nm and 1.63 AE 0.04 μm, respectively, for a slit of three pixels.As expected, resolution along the slit (x) was independent of the slit size.According to theoretical calculations, with the slit fully open, the lateral and axial FWHMs are equal to 328 nm and 1.69 μm, respectively.As the silt size decreases, the FWHMs along the y and z axes decrease to 274 nm and 1.56 μm, respectively, for a slit of three pixels.The experimental values for the FWHMs are consistent with the theoretical predictions (solid lines) albeit slightly bigger, especially for larger slits (> 10 pixels) probably due to spherical aberrations.Since the rolling shutter rejects photons from out-of-focus planes, its contribution to increase image contrast was evaluated by analyzing the ratio between signal and background.The same bead was imaged successively with six slit sizes [3, 6, 12, 30, 60, and 1000 pixels, Fig. 1(b)] and its xy profile was fitted by a two-dimensional Gaussian function [Fig.2(c)].The images were corrected by subtracting the dark noise of the camera for each slit size.The normalized contrast, averaged over 10 beads, is shown in Fig. 2(d).With the slit corresponding to a global shutter (1000 pixels), a normalized contrast of 0.26 AE 0.03 was measured.By gradually closing down the slit, more out-of-focus photons were rejected and the contrast increased up to its maximum at a slit size of 12 pixels, corresponding to approximately three PSFs.At this optimal shutter size, contrast was improved fourfold compared with the global shutter.When the slit size was further reduced, the rolling shutter not only rejected background photons but also signal photons and therefore contrast began to drop.The theoretical contrast, calculated from the normalized ratio of the theoretical signal extracted from the PSF calculation and the experimental background extracted from the bead images, was in agreement with the experimentally measured data points [Fig.1(d), solid line].
After a comprehensive evaluation of resolution and contrast of the microscope on fluorescent beads, the imaging performance was tested on semicleared brain samples, which in addition to increased scattering, also exhibit macroscopic mismatches in refractive index.Adult mice (Thy1-GFPM 15 ) were prepared and their brains optically cleared according to the protocol described in Costantini et al. 16 With this clearing methodology, we achieved transmittances of 30% to 40% at 900 nm resulting in a refractive index compatible with our microscope objective.After overnight fixation in paraformaldehyde (PFA), brains were cut in 2-mm-thick slices and cleared by two successive incubations.The slices were put in a 20-mL PSB solution containing 30% and 63% of TDE for 1 and 3 h while rotating at room temperature. 16This procedure reduces tissue scattering in samples, which were previously fixed in PFA; however, without rendering them entirely transparent.All experimental protocols were designed in accordance with Italian laws and were approved by the Italian Minister of Health (authorization no.790/2016-PR) (Fig. 3).
Imaging performance in tissue was evaluated from z stacks of 400-μm depth (step size 2 μm) corresponding to the working distance of the microscope objective.For each depth, six images corresponding to 3, 6, 12, 30, 60, and 1000 pixels slits were successively acquired.Resolution in tissue was estimated by fitting a line profile through the neck of dendritic spines located throughout the stack [Figs.2(a) and 2(b)] with a Gaussian curve [Fig.2(c)] and comparing the FWHM [n ¼ 11, error bars are standard error of the mean (SEM), Fig. 2(d)].The reported FWHMs are slightly elevated compared with values obtained with beads due to increased scattering in the semicleared sample and spine necks not necessarily being structures under the resolution limit of the microscope.As a general trend, resolution decreased with increasing shutter size.
Image contrast was quantified at each depth by the normalized Shannon entropy of the discrete cosine transform. 17or a comparative study, six regions of interest were selected over the entire FOV [Fig.2(d)] and analyzed through depth  for all shutters.The dark stripe visible in the full FOV imaging in Fig. 2(e) is due to a suboptimal DOE, which has a lower intensity in the zero-order spot.The contrast values of the six ROIs were averaged at each depth and then binned over 30 μm.We normalized each average contrast value by the global shutter contrast (corresponding to a slit size of 1000 pixels) and reported its evolution with depth for all shutter sizes as shown in Fig. 2(g).The general trend of increased scattering and therefore reduced contrast with increasing shutter size and increasing depth was also evident in tissue.Similar to the imaging of beads, the best contrast was achieved with a slit size of 12 pixels corresponding to three PSFs throughout the entire depths of the stack.
In summary, we have presented a mixed microscope, which operates with a robust and fast multispot two-photon excitation and a cross talk rejecting wide-field detection.The microscope was fully characterized by imaging subresolution sized fluorescent beads and semicleared mouse brains.However, we point out that the proposed method is poorly compatible with uncleared samples and in vivo measurements.
The frame rate of the microscope is currently limited by the galvo speed to ∼100-μs per pixel line (corresponding to 10 Hz for a 1000 × 1000 pixels image).A future implementation would benefit from a resonant scan system, which would shift the bottleneck to the camera line readout time (on the order of 10 μs), therefore corresponding to 100 Hz for a 1000 × 1000 pixels image.The gain in speed scales essentially with the number of multiplexed excitation spots as long as the laser source is capable of providing the same energy per beamlet.With current Ti:Sa lasers, up to 20 beamlets of enough power (75 mW) can be generated, resulting in a 20× increase in imaging speed.By comparison, a typical two-photon microscope with a pixel dwell time of 1 μs achieves a frame rate of 1 Hz for the same size image.Considering that the 3-D reconstruction of large volumes such as 2 mm × 6 mm × 0.5 mm of brain tissue requires 5 days with a standard TP system, a 20× increase in acquisition speed is highly desirable.Whereas a conventional MMM struggles with optical cross talk present in scattering tissues even at shallow depths, in the microscope presented here suppression of cross talk with a confocal rolling shutter was sufficient to distinguish dendritic spines down to a depth of 400 μm in semicleared mouse brains.

Fig. 1
Fig. 1 (a) Schematic of the microscope setup, (b) scheme of the rolling shutter as a synchronized confocal slit on the sCMOS camera sensor, and (c) fluorescent pollen grain for three different microscope configurations.From left to right: a single beam is scanned, 15 beams are scanned, 15 beams are scanned with a synchronized 10 pixels slit.

Fig. 2
Fig. 2 (a) FWHM of the PSF as a function of slit size in x , y, and z (n ¼ 20 beads, errors are SD, solid line: theoretical prediction).(b) Fluorescent bead for six slit sizes.Laser power: 740 mW, line exposure time: 1.6 ms.(c) y intensity profiles and Gaussian fit of bead shown in (d).The amplitude of the Gaussian fits S and the corresponding background B was used to calculate contrast in (d).(d) Normalized contrast as a function of slit size.(n ¼ 10 beads, error is SD, solid line: theoretical prediction.)

Fig. 3
Fig. 3 (a) and (b) Exemplary line profiles through the necks of dendritic spines.Scale bar 5 μm.(c) Intensity profiles for the line profile shown in (a) for all shutter sizes.(d) FWHM of Gaussian fits.(n ¼ 11, error bars are SEM.)(e) Maximum intensity projection of the full FOV through the entire depth of the stack (shutter 12 pixels).Indicated are the six ROIs (25 μm × 25 μm) used in the contrast analysis.Scale bar 15 μm.(f) Exemplary images at different depths and different shutter sizes.Brightness has been adjusted for each image.Scale bar 5 μm.(g) Evolution of relative contrast, calculated as the normalized Shannon entropy of the discrete cosine transform with depths for all shutter sizes binned over 30 μm, and normalized by the global shutter contrast (corresponding to a slit of 1000 pixels).Error is SEM.