Three-dimensional spatiotemporal focusing of holographic patterns

Abstract

Two-photon excitation with temporally focused pulses can be combined with phase-modulation approaches, such as computer-generated holography and generalized phase contrast, to efficiently distribute light into two-dimensional, axially confined, user-defined shapes. Adding lens-phase modulations to 2D-phase holograms enables remote axial pattern displacement as well as simultaneous pattern generation in multiple distinct planes. However, the axial confinement linearly degrades with lateral shape area in previous reports where axially shifted holographic shapes were not temporally focused. Here we report an optical system using two spatial light modulators to independently control transverse- and axial-target light distribution. This approach enables simultaneous axial translation of single or multiple spatiotemporally focused patterns across the sample volume while achieving the axial confinement of temporal focusing. We use the system's capability to photoconvert tens of Kaede-expressing neurons with single-cell resolution in live zebrafish larvae.

ARTICLE

Received 10 Nov 2015 | Accepted 12 May 2016 | Published 16 Jun 2016

Two-photon excitation with temporally focused pulses can be combined with phase-modulation approaches, such as computer-generated holography and generalized phase contrast, to efciently distribute light into two-dimensional, axially conned, user-dened shapes. Adding lens-phase modulations to 2D-phase holograms enables remote axial pattern displacement as well as simultaneous pattern generation in multiple distinct planes. However, the axial connement linearly degrades with lateral shape area in previous reports where axially shifted holographic shapes were not temporally focused. Here we report an optical system using two spatial light modulators to independently control transverse- and axial-target light distribution. This approach enables simultaneous axial translation of single or multiple spatiotemporally focused patterns across the sample volume while achieving the axial connement of temporal focusing. We use the systems capability to photoconvert tens of Kaede-expressing neurons with single-cell resolution in live zebrash larvae.

DOI: 10.1038/ncomms11928 OPEN

Three-dimensional spatiotemporal focusing of holographic patterns

Oscar Hernandez1,w, Eirini Papagiakoumou1,2, Dimitrii Tanese1, Kevin Fidelin3, Claire Wyart3 & Valentina Emiliani1

1 Wavefront-Engineering Microscopy Group, Neurophotonics Laboratory, CNRS UMR 8250, Paris Descartes University, UFR Biomdicale, 45 rue des Saints-Pres, 75270 Paris Cedex 06, France. 2 Institut national de la sant et de la recherche mdicale (Inserm), France. 3 Institut du Cerveau et de la Moellepinire, UPMC, Inserm UMR S975, CNRS UMR 7225, Campus Hospitalier Piti Salptrire, 47 building de lHpital, 75013 Paris, France. w Present address: CNC

Program, Stanford University, Stanford, California 94305, USA. Correspondence and requests for materials should be addressed to V.E. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928

Since its rst demonstration1, nonlinear two-photon (2P) microscopy has revolutionized diverse research elds including neuronal structural and functional imaging,

photostimulation, laser processing and lithography. Original congurations were based on scanning a tightly focused pulsed laser beam across a dened sample region with galvanometric mirrors, acousto-optic deectors24 or resonant scanners5. Recently, applications including lithography, uncaging, optogenetics and fast functional imaging have motivated scanless 2P-excitation method development, enabling simultaneous illumination of laterally extended regions of interest while preserving micrometre axial connement. This was achieved by axially conning large spots generated by low-numerical aperture (NA) Gaussian beams with temporal focusing (TF)6,7. In TF, a diffraction grating conjugated to the sample plane diffracts the different spectral frequencies comprising an ultrashort excitation pulse towards different directions. The various frequencies thus propagate towards the objective focal plane at different angles, such that spatiotemporal coupling effects8 lead to temporal pulse broadening above and below the focal plane, which remains the only region irradiated at peak powers efcient for 2P excitation. In this way, TF enables axial connement equivalent to that of line-scanning microscopy that is independent of the excitation beams lateral extent. Investigators have applied TF to imaging913, functional imaging14,15, super-resolution imaging16, lithography17 and neuronal photostimulation18,19 by generating circular-spot shapes only.

As a more exible way to elaborate light patterning, phase-modulation approaches using liquid-crystal spatial light modulators, such as computer-generated holography (CGH)20,21 and generalized phase contrast (GPC)22 have demonstrated efcient light sculpting in two-dimensional (2D) user-dened shapes. Combined with TF to axially conne the shapes, CGH and GPC can sculpt 2P-excitation volumes with the micrometre precision needed to illuminate small structures, such as neuronal dendrites22,23. These methods also enable 2P optogenetic stimulation with high temporal resolution and robustness to scattering23,24.

CGH can also sculpt light in three dimensions (3D)25, a feature used to generate multidiffraction-limited traps for optical tweezers26,27 and 3D glutamate uncaging28,29. Although not yet demonstrated with 2P excitation, adding lens-phase modulations to 2D-phase holograms also enables remote axial displacement and 3D positioning of laterally shaped targets27,30,31.

Combining CGHs 3D light-shaping capability with 2P excitation can signicantly broaden the range of possible applications. For example, simultaneous multiplane holographic pattern generation can enable fast 2P multiplane imaging, photostimulation and photo-polymerization methods. In addition, remote axial displacement of holographic targets can couple holographic illumination with a second imaging or stimulation channel, providing independent control of their respective focal planes32, as well as remote volume scanning. However, previous 3D CGH optical congurations could not be implemented with TF because axially shifted excitation planes cannot be simultaneously imaged on the TF grating. This shortcoming restricted CGH-TF to 2D patterns focused at the objective focal plane21.

Here we demonstrate a unique optical system overcoming this limitation. The system achieves remote axial displacement of temporally focused holographic beams as well as multiple temporally focused planes by shaping the incoming wavefront in two steps using two spatial light modulators (SLMs). A rst SLM laterally shapes the target light distribution that is focused on the TF grating, while a second SLM, positioned after the

grating, controls the axial position(s) of the spatiotemporal focal plane(s) in the sample volume. By integrating phase proles that minimize optical aberrations and intensity compensation protocols, we generate spatiotemporally focused patterns with uniform light distribution throughout the entire accessible volume, and demonstrate generation of laterally shaped targets across an unprecedented axial range. We apply axially conned multiplane light illumination to photoconvert in vivo Kaede-protein-expressing neurons in zebrash larvae. The system enables cellular resolution photoconversion of tens of spinal cord neurons occupying separate axial planes.

ResultsOptical system. In the original optical system we designed (Fig. 1, see Methods), the output of a Ti:Sapphire laser was expanded to illuminate a rst liquid crystal on silicon SLM (LCOS-SLM; SLM1) used to generate 2D holographic light patterns at the focal plane of a rst lens where a blazed grating (G) was placed for TF. TF enables 2D CGH pattern generation with the axial connement of a line-scanning microscope, practically decoupling axial connement from lateral extent21 (Supplementary Fig. 1). Two telescopes (L2, L3; L4, OBJ1) relayed the holographic pattern on the sample plane. The phase holograms addressed on SLM1 were calculated using a Gerchberg and Saxton (GS)-based algorithm33. A second LCOS-SLM (SLM2) was placed at the confocal plane of the rst telescope and imaged at the back focal plane of the microscope objective through the second telescope. By addressing SLM2 with single or multiple lens phases, we axially displaced single or multiple spatiotemporally focused patterns across the sample volume.

In order to characterize the optical properties of the holographic patterns, we used a microscope conguration using two opposite-facing objectives20,21. By doing so, holographic patterns illuminated a thin uorescent lm and 2P-excited uorescence was collected by a second objective (OBJ2) placed opposite to OBJ1 and imaged on a CCD (charged couple device) camera. For 3D reconstruction of illumination volumes, OBJ2 was xed and focused on the uorescent layer, while OBJ1 was

BE

SLM1

CL L1

L2

Ti:Sapphire laser

G

CCD L5 OBJ2 FFP OBJ1 L4 L3 SLM2

Figure 1 | Experimental set-up for 3D-CGH-TF. The output beam of a Ti:Sapphire laser is magnied using a beam expander (BE) and projected on a rst SLM (SLM1). SLM1 modulates the beams phase so that light forms a user-dened intensity pattern on the diffraction grating (G) after passing through lens L1. The rst diffraction order is collimated by the lens L2 and directed on a second SLM (SLM2). SLM2 imprints a lens-phase modulation that enables precise axial positioning of the spatiotemporal focal plane. The laser beam is then relayed and scaled by lenses L3 and L4 to the excitation objective (OBJ1) pupil size. OBJ1 is mounted on a piezo positioner so that it focuses and axially scans the excitation beam across a thin uorescent layer. A second objective (OBJ2), always focused on the uorescent layer, collects emitted uorescence and forms an image on a CCD camera. Two cross-oriented cylindrical lenses (CL), with focal lengths of equal power and opposite sign, are used to suppress the zero-order spot of the rst SLM64.

2 NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928 ARTICLE

moved along the axial direction with a piezo-scanning stage20,21. For experiments on characterizing SLMs diffraction efciency, the holographic patterns were directly imaged on the CCD camera (image of the transmitted laser light). For the experiments of photoconversion and optogenetic stimulation in zebrash larvae, the system was coupled to a HiLo imaging set-up (see following sections).

The optical system enabled remote axial displacement of temporally focused holographic patterns as well as generation of multiple spatiotemporally focused holographic targets in distinct axial planes. In a simplied conguration where mirrors replaced the SLM2 and the grating, the system also generated 2P multiplane arbitrarily shaped non-temporally focused patterns.

Remote axial shift of spatiotemporally focused patterns. First, we demonstrated remote axial displacement of a temporally focused holographic pattern by separating input-beam-phase modulation in two steps, controlling rst the target lateral light distribution, and, second, its axial position. In contrast to conventional CGH, axial displacement of temporally focused holographic patterns cannot be achieved by simply adding a Fresnel lens to the phase hologram34 because the lens effect, by axially displacing the excitation plane with respect to the TF grating, would forbid the spatial and temporal focal planes to coincide35. We resolved this issue by implementing a novel, two-SLM strategy. SLM1 generated 2D CGH illumination patterns focused on the grating G, which dispersed the spectral components of the illumination pattern on SLM2. SLM2, conjugated to the objective back focal plane, was addressed with a Fresnel lens-phase prole to control the targets axial position in the sample volume. This design enables the spatial and temporal focal planes to coincide at the grating, and to be jointly translated by SLM2 across the sample volume axial extent.

The effect of spherical aberrations on targets generated out of OBJ1s nominal focal plane (Supplementary Fig. 2) was minimized by describing the objective focal sphere36 within the approximation of small defocus and high NA. We thus addressed SLM2 with lens-phase proles featuring spherical phase28,37:

f r

kDz

n2 NA2r2

q 1 where k is the free-space wavenumber, n is the refractive index of the immersion medium, NA is the numerical aperture of the objective, r

rNA feq (with feq being the equivalent focal length: feq f3f4 fobj and r

x2 y2 z2

p ) is the normalized pupil at the

SLM2 plane and Dz the axial displacement of the holographic light pattern in the sample volume.

With the above strategy, we demonstrated remote axial translation of temporally focused holographic patterns by displacing a 20-mm-diameter holographic spot throughout a 130-mm axial range (Fig. 2ac), the accessible eld of excitation (FOEz) without hologram aliasing (Supplementary Note 1). Within this range, we could remotely displace holographic patterns while conserving target shape sharpness (Fig. 2a, bottom) and axial connement (510 mm; Fig. 2b,c). The integrated 2P uorescence intensity decreased by approximately vefold at the edges because of the diffraction efciency of SLM2 (Supplementary Note 1 and Supplementary Fig. 3ac).

Furthermore, we generated targets with micrometre axial connement (1025 mm) and sharp shapes across a broader axial range (300 mm, that is, approximately two times greater than the FOEz). However, in this extended range, the full-width at half-maximum (FWHM) broadened (B25 mm at 300 mm; Fig. 2b,c), uorescence intensity decreased (B10 times at 300 mm) and spot quality deteriorated (Supplementary Fig. 4)

as a consequence of optical aberrations, SLM2 diffraction efciency and the nite dimensions of mirrors and lenses placed after SLM2 that can vignette beam edges for highly diverging/ converging wavefronts, among other factors.

For applications requiring constant illumination across the axial range, attenuating illumination intensity at axial positions near the centre of the FOE can compensate the position-dependent intensity decrease (Supplementary Fig. 3g). This can be accomplished either by reducing the incident laser power or by redirecting light into an extra spot similar to the procedure described in Supplementary Fig. 3df.

3D-CGH-TF spatiotemporally focused patterns. Next, we demonstrated that decoupling lateral and axial light shaping into two separate steps also enables generation of multiplane spatio-temporally focused patterns. In this conguration, SLM1 is tiled into n vertical regions (that is, parallel to the orientation of the grating lines and orthogonal to the grating linear dispersion), with n equal to the number of axial planes populated by the nal illumination pattern. Each tile generates a 2D target shape focused through L1 on the TF grating (Fig. 3a, top, Supplementary Fig. 5). SLM2, vertically tiled into n independent lens-phase proles, then axially displaces each shape into its target axial plane (Fig. 3a, bottom). In this conguration, for each of the n targeted planes, both spatial and temporal foci coincide at the grating, while the n phase-lens proles addressed on SLM2 enable independent remote displacement of each shaped target in the sample volume (Fig. 3b,c).

In agreement with previous observations38, reducing the hologram width, Dy, along the direction perpendicular to the gratings dispersion did not compromise spot axial connement, as the reduction of the hologram width, Dx, in the orthogonal direction would do (Fig. 3d). Moreover, tiling holograms along the x axis would create crosstalk between axially separated planes because of the dispersion at the grating that generates at SLM2 a lateral spatial overlapping of the holograms generated by SLM1, making impossible to imprint independent axial shifts to the corresponding patterns.

However, vertical hologram resizing introduced a spatial lter that vertically elongated speckles (Fig. 3e). In particular, the vertical mean speckle size, estimated by calculating the vertical autocorrelation width39, sy, increased by a factor of two with respect to the original size (sy 2symin) when the hologram

vertical dimension was reduced to one-fourth of the full SLM aperture (Fig. 3f). In contrast, the horizontal autocorrelation width sx was unaffected by vertical hologram resizing (Fig. 3f).

The limited size of the SLM can constitute the limiting factor dening the maximum number of axially separated planes in 3D-CGH-TF. Theoretically, the maximum number of planes is determined by the number of pixel rows on the SLM, but practically there is a tradeoff between the number of axial planes and the lateral resolution requirements for each application. As far as the spot shape is not distorted and spots contour remains well dened, the smoothing of the excitation holographic spot to some extent can be even benecial21,38,40.

Moreover, vertical SLM tiling induced a lateral tilt on the intensity propagation due to the asymmetric illumination of the objective back focal plane. Addressing SLM1 and SLM2 symmetrically along the vertical axis of the SLM eliminated this effect (Supplementary Fig. 6). Hologram resizing also affected the intensity of the illumination pattern at the sample plane. Specically, holograms projected at the centre of the objective back aperture generated brighter targets than holograms projected on the side (Supplementary Fig. 7 and Supplementary Note 2). This effect adds to the position-dependent diffraction

NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928

a b

100 m

z = 50 m 0 m 50 m 100 m

1

0 1.00.8

0.6

0.4

0.2

0.0

30

25

20

15

10

FOEz

Axial position (m)

Fluorescence intensity

(a.u.)

Axial position (m)

c

FOEz

Axial confinement

FWHM (m)

z y

x y

CGG 70

60 50 40 30 20 10 0

200 300

5

300 200 100 0 100 200 300

0 300 200 100 0 100

CGH-TF

Figure 2 | Axial displacement of spatiotemporally focused patterns. (a) Axial displacement of a 20-mm-diameter temporally focused holographic spot. Top, orthogonal maximum uorescence intensity projection of the spot, for axial displacements of 50 and100 mm from the focal plane (0 mm). Bottom, corresponding xy uorescence intensity cross-sections. Scale bar, 20 mm. The colour bar refers to normalized intensity. (b) Axial prole of the integrated uorescence intensity of the 20-mm-diameter holographic spot for different axial displacements. (c) Axial connement (FWHM) of the proles shown in b (CGH-TF; red data) compared with the axial connement of a 20-mm-diameter holographic spot without TF (CGH; blue data). Data were tted with a parabolic function (dashed lines) in both cases. White area in b,c represents the eld of excitation (FOEz).

efciency determined by the SLM pixel size (Supplementary Note 1 and Supplementary Fig. 3). Thus, generation of homogeneous 3D-CGH-TF patterns required compensating both lateral and axial position-dependent intensity variations, which in both cases was higher at the FOE centre and falling off towards the periphery. For lateral shape generation with SLM1, we weighted target intensity input to the GS algorithm such that targets occupying low-efciency regions were brighter than targets position at the efcient FOEx,y central zone (Supplementary

Fig. 3). To achieve uniformity in the axial direction, we scaled the vertical tile size, allocating greater area to holograms projected at the objective pupil periphery with respect to the centre (Supplementary Note 2 and Supplementary Fig. 7). This enabled generation of homogeneous light patterns (Fig. 3g) across the whole excitation volume. Holographic patterns generated in distinct planes yielded uorescence intensity distributions showing that axial connement was also well conserved (Fig. 3g). However, an unavoidable background in the intermediate planes appeared when targets were laterally aligned (Fig. 3g, right bottom panel, and Supplementary Fig. 8).

3D computer-generated holography. We demonstrate a simplied, single-SLM system for non-temporally focused multiplane pattern generation, replacing the grating and second SLM with mirrors (Fig. 4a); we patterned multiple targets occupying separate axial planes (Fig. 4b,c, Supplementary Movie 1) with a modied GS algorithm for multiplane pattern projection41 (Supplementary Fig. 9). By using the spherical phase expression described in equation (1) to translate the various targets into different axial planes, with r

rNA feq1 where feq1 f1f3f2f4 fobj,

we minimized the effects of spherical aberrations on targets generated out of the focal plane, improving agreement between experimental and expected axial positions for the different targets compared with 3D-CGH implemented with parabolic-lens phases42 (Supplementary Fig. 2). We compensated diffraction

efciency-induced intensity variations by weighting target intensity inputs to the multiplane GS algorithm, as previously described (Supplementary Fig. 3), achieving uniform intensity shape generation throughout the accessible B240 240 260 mm3 volume (Supplementary Note 1).

In vivo 3D-patterned selective photoconversion of neurons. We leveraged the axial specicity achievable with 3D-CGH-TF to photoconvert Kaede protein-expressing neurons in live Tg(HuC:Gal4; UAS:Kaede) double transgenic zebrash larvae43. Kaede is a green photoactivable uorescent protein that, when exposed to ultraviolet light, undergoes a photo-induced protein cleavage, red shifting its uorescence emission spectrum44. With 2P excitation, efcient Kaede photoconversion can be obtained for wavelengths in the range of 760800 nm45. We rst performed Kaede photoconversion in spinal neurons (Fig. 5a, blue and red areas) to demonstrate the ability for single-cell resolution. Next, we demonstrated the robustness of 3D-CGH-TF to scattering by moving to the brain (Fig. 5a, dark yellow area) that is a more scattering region (Supplementary Fig. 10).

To monitor photoconversion in the spinal cord, we combined multiplane patterned photostimulation with a two-colour HiLo (high/low-frequency sequential acquisition) imaging system46 (see Methods and Supplementary Fig. 11). First, a HiLo z-stack exciting green uorescence (excitation wavelength lexc 473 nm)

was acquired to localize Kaede-expressing neurons from which we selected photoconversion targets. We then photoconverted selected cells with 2P holographically patterned illumination (stimulation time 0.20100 s; excitation power density:0.034.0 mW mm 2; Methods. We note that power densities are always given relatively to the area of the spots surface). A second HiLo z-stack exciting red uorescence (lexc 561 nm)

was acquired after the photostimulation protocol in order to assess the photoconversion efciency. For photoconversion experiments in the brain, scattering quickly deteriorates the quality of HiLo images. We therefore monitored the extent of

4 NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928 ARTICLE

a b

Plane A

[afii9843]

[afii9843]

0

1

0

SLM 1 SLM 2

Plane A

Plane B

Plane APlane B

x

z

y

x

Plane B

x y

z y

c

1.0

0.2

z

x

d e f

12 mm 10 mm 8 mm

x

y

Axial confinement

FWHM (m)

8 2.0

1

0 1.5

1.0

0.5 2 4 6 8 10 12 Hologram width (mm)

x y

7

6

1.0

0.5

0.0100 50 50 100

Axial position (m)

0 100 50 50 100Axial position (m)

0 100 50 50 100Axial position (m)

0

Autocorrelation width

FWHM (m)

6 mm 2 mm

4 mm

5 0 2 4 6 8 10 12

Hologram width (mm)

y

x

g

Fluorescence

intensity (a.u.)

Fluorescence

intensity (a.u.)

1.0

0.5

0.0

Fluorescence

intensity (a.u.)

1.0

0.5

0.0

Fluorescence

intensity (a.u.)

1.0

0.5

0.0

Fluorescence

intensity (a.u.)

1.0

0.5

0.0

100 50 50 100 Axial position (m)

0

100 50 50 100 Axial position (m)

0

Figure 3 | Multiplane spatiotemporally focused pattern generation. (a) Top, tiled phase proles addressed to SLM1 for encoding the words neuro (planeA) and photonics (plane B). Bottom, Fresnel lens-phase proles addressed to SLM2 to axially displace each holographic pattern generated by SLM1 on separated planes at 20 mm (plane A) and 20 mm (plane B). (b) Left, xy 2P uorescence intensity cross-sections at planes A and B generated by the

holograms in a. Right, orthogonal maximum 2P intensity projection along x (top) and y (bottom). (c) Orthogonal uorescence intensity projection of spatiotemporally focused patterns created by projecting 20-mm-diameter holographic spots in one, two, three and four planes at positions laterally shifted and in four planes at positions axially aligned (from left to right). Scale bars, 20 mm. (d) Axial connement (FWHM) of a 20-mm-diameter holographic spot for different hologram widths, tested in both directions, Dx (blue data) and Dy (red data) parallel and perpendicular to gratings dispersion, respectively.

(e) xy 2P uorescence intensity cross-sections of a 20-mm-diameter holographic spot generated with different hologram widths Dy. (f) Autocorrelation width at the sample plane of a 20-mm-diameter holographic spot as a function of the hologram width along the x- (blue data; sx) and y directions (red data; sy). (g) Axial prole of the integrated uorescence intensity of the holographic spots shown on c. Colour bars refer to normalized intensity.

NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928

a c

BE SLM1

Ti:Sapphire laser

CL

L1

L2

150

100

50

50

100

150

z x

M

z = 130 m

CCD L5 FFP

M

OBJ2

OBJ1

L4

L3

b

m

0

50

50

z = 130 m

z = 65 m

z = 0 m

z = 65 m

0

m

0

y

m

50

Figure 4 | 3D-CGH. (a) Schematic of the optical set-up for 3D-CGH. In this case the diffraction grating G and SLM2 were replaced by mirrors. (b) 2P uorescence images of 3D-CGH patterns, depicting the letters a, b and c at three different axial positions, z 130 mm, 0 and 130 mm,

respectively. The phase prole used to project these patterns was calculated using a multiplane GS algorithm (Supplementary Fig. 9). Weighting of the input patterns according to their lateral and axial position (Supplementary Fig. 3) enabled diffraction efciency correction and generated equal-intensity light patterns. (c) Volumetric reconstruction of three-dimensional distribution of 5-mm-diameter holographic spots (see also Supplementary Movie 1).

photoconversion by moving the sample to a 2P galvo-based scanning imaging system in another set-up (see Methods). Green and red uorescence in the 2P imaging system were excited at 780 nm, and scanning was performed at 0.74 Hz.

In the spinal cord, we rst used 2D-CGH-TF to demonstrate simultaneous photoconversion (lphot 800 nm) of multiple

neurons in a single plane with cellular precision (Fig. 5a,b). Then, we used 3D-CGH-TF to photoconvert isolated neurons (Fig. 5c) or groups of multiple neurons (Fig. 5d) in distinct axial planes. The axial connement achievable with 3D-CGH-TF enabled single-layer selectivity that would otherwise not be possible using CGH deprived of TF (Supplementary Fig. 12a). In all experiments performed, photoconversion induced a B10-fold increase in the ratio of pre- and post-photoconversion red uorescence.

To demonstrate in-depth photoconversion, we performed a second set of experiments in the brain of zebrash larvae by photoconverting neuronal ensembles in two different axial planes using 35-mm-diameter temporally focused holographic spots separated by B80 mm (Fig. 5e), with the rst plane positioned at B90 mm from the surface. Although the axial connement deteriorates with depth because of scattering, the red uorescence intensity proles along the axial direction showed selective and efcient photoconversion (Fig. 5d,e), which was not achievable with 3D-CGH alone (Supplementary Fig. 12b). Strikingly, axial resolution and overall shape were also preserved despite light propagation through the zebrash brain (Supplementary Fig. 10). In agreement with previous results in rats23,24, these results demonstrate that TF enables robust in-depth propagation of shaped patterns through scattering media.

It is important to note that 2P Kaede photoconversion is a low-efciency process and requires relatively high illumination doses (B0.044 mW mm 2) and long (from 200 ms to a few hundreds of seconds) stimulation protocols. Other type of applications such as activation of neurons via 2P-mediated optogenetic stimulation can work with drastically reduced

photostimulation time (5 pulses of 50 ms), as well as lower illumination doses (0.040.60 mW mm 2; Supplementary Note 3).

A still open question for multiplane multispot photo-stimulation concerns the possible temperature rise induced by laser illumination. For the illumination doses we used in the optogenetics experiments (Supplementary Note 3), we estimated that this does not exceed few degrees Celsius. Specically, by solving the heat diffusion equation47 with typical thermal optical parameters of tissue found in literature4850, we estimated the mean temperature rise to be of the order of 0.060.25 C, using illumination pulses of 50 ms in the range of 0.040.6 mW mm 2 and a 10-mm diameter spot (Picot et al., private communication). Using 10 spots of the same size distributed randomly over a FOE of 100 100 mm2 would generate a mean temperature rise of

0.93.7 C. Diffusion over the tissue volume enables the temperature to return to its equilibrium value in less than 150 ms, which was the interval between illumination pulses (Picot et al., private communication).

DiscussionWe have demonstrated a unique optical system enabling remote axial displacement of temporally focused holographic patterns, as well as generation of multiple temporally focused holographic targets occupying separate axial planes. In our two-step system, the rst SLM is addressed with phase holograms controlling the transverse target light distribution. The second SLM, positioned after the TF grating, is addressed with Fresnel lens-phase functions and controls target axial position. We demonstrated that this conguration can jointly translate single or multiple spatiotemporally focused patterns across the sample volume.

We demonstrated axial displacement of a single temporally focused holographic pattern across an axial range (300 mm), roughly two times greater than the nominal accessible axial FOE, FOEz (130 mm). For axial shifts of Dz

j j FOEZ=2, spot

shape and axial resolution (510 mm) were well conserved, while uorescence intensity decreased approximately vefold

6 NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928 ARTICLE

a

d

3D-CGH-TF

z = 18 m

Fluorescence intensity

1.0

0.5

0.050 50 Axial position (m)

b

2D-CGH-TF

25 25

0

z = 3 m

Fluorescence intensity

1.0

0.5

0.050 50 Axial position (m)

y

x

25 25

0

Fluorescence intensity

c

1.0

0.5

0.050 50 Axial position (m)

84 (m)

z = 12 m

3D-CGH-TF

z = 1 m z = 9 m

25 25

0

z

x

z

x

y

x

e

3D-CGH-TF

z = 34 m

z = 50 m

Fluorescence intensity

1.0

0.5

0.0100 50 50

Axial position (m)

z = 84 m

z = 18 m

y

x

z

x

0

100

Figure 5 | 3D simultaneous 2P photoconversion of Kaede in vivo. (a) Merged brighteld and wideeld uorescence images of a double transgenic Tg(HuC:gal4; UAS:kaede) zebrash larvae. Red and blue squares represent the approximate areas where we performed photoconversion. Scale bar, 400 mm.

(b) Left, overlaid green and red HiLo uorescence images before and after photoconversion, respectively. Right, orthogonal maximum red uorescence intensity projection showing 14 photoconverted neurons on a single axial plane (illumination density 0.4 mW mm 2, 200 pulses of 50 ms). Scale bars, 60 mm. (c) Orthogonal maximum uorescence intensity projection of overlaid HiLo pre- and post-photoconversion images (green and red uorescence, respectively). Three single cells were photoconverted on separated axial planes (4.0 mW mm 2, one pulse of 200 ms). Scale bar, 60 mm. (d) Simultaneous 3D photoconversion of neural ensembles in the spinal cord. Left, overlaid HiLo pre- and post-photoconversion uorescence images, where three 35-mm-diameter holographic spots projected at z 18 mm, 3 and 18 mm were used for photoconversion (0.03 mW mm 2, 2,000 pulses of 50 ms).

Right, axial distributions of green pre- and red post-photoconversion integrated uorescence intensity over z for the spots projected at the three different planes. Scale bars, 60 mm. (e) Simultaneous 3D photoconversion of neural ensembles in the zebrash brain. Left, overlaid 2P-excited green- and red post-photoconversion uorescence images, where two 35-mm-diameter holographic spots projected at z 38 mm and 50 mm were used for

photoconversion (0.11 mW mm 2, 9,000 pulses of 50 ms). Scale bar, 20 mm. Middle, Orthogonal maximum 2P-excited uorescence intensity projection of overlaid green and red-post-photoconversion images. Scale bar, 20 mm. Right, axial distributions of green pre- and red post-photoconversion integrated 2P uorescence intensity over z for the spot at z 38 mm (solid lines) and the one at z 50 mm (dotted lines). z-values in all cases are given as distances

from the focal plane of the objective, which for the spinal cord experiments was at B60 mm and for the brain at B90 mm from the sh surface (where green uorescence was starting). Positive z-values are closer to the surface. lphot 800 nm.

(consistent with the roughly twofold decrease in illumination intensity, Supplementary Fig. 3). Displacements up to twice the accessible axial FOE were also possible, although requiring compensation mechanisms to correct for broadening of axial connement (B25 mm for Dz 300 mm) and intensity losses.

With our approach, applications requiring continuous fast scanning of a single spatiotemporally focused target will be limited by the SLM refresh rate (60200 Hz). Replacing SLM2

with a tunable lens51 would overcome this limitation. However, it would not enable generation of multiple spatiotemporally focused targets at distinct axial planes.

Alternative solutions for axial displacement of a temporally focused shape used either variable group velocity dispersion (GVD)38,5255 or mechanical axial translation of the TF grating35. Indeed, introducing variable GVD to the input laser pulse enabled efcient axial displacement of a Gaussian beams temporal focal

NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928

plane. However, GVD-induced axial displacement strongly depends on the autocorrelation width of the illumination patterns, which for holographic spots is small due to speckles. Therefore, in order to achieve axial displacement of temporally focused holographic patterns, spatial ltering is required to increase the correlation width. However, this degrades spatial resolution and decreases illumination intensity38. Furthermore, the GVD approach achieved axial displacements of only a few microns38. Greater scanning ranges could be obtained by combining 2D-CGH with the optical design proposed by Dana et al.35 for axial scanning of a temporally focused line. The system implemented on-axis-light propagation and mechanical axial translation of the TF grating. However, for high-magnication(M) objectives, the axial shift, d, was limited by the long required grating translation, Dpd M2. Thus, the maximum shift

achieved with a 40 objective was inferior to 30 mm.

In addition to limited range, neither GVD shift nor grating translation enable the generation of multiple temporally focused targets at distinct axial planes. Here we overcome this limitation with a two-step, two-SLM system, addressing a rst SLM with multiple vertically tiled phase holograms, each encoding the light distribution of a single plane, and a second SLM with an equal number of Fresnel lens phases, which individually control the axial position for each plane. In this way, all targets are projected to the TF plane, which is kept at a xed position, while the second SLM imposes the spatial wavefront curvature needed to displace each plane axially.

The number of spatiotemporal focal planes that can be generated with this design is a tradeoff between the number of available pixels and spatial resolution. Here we demonstrated generation of up to four spatiotemporal focal planes. In this case, the vertical spatial resolution (that is, the average speckle size in the y axis), where tiling occurs, decreases roughly twofold because of the reduced size of the hologram at the objective back aperture. LCOS-SLM devices with increased number of pixels should increase the number of achievable planes and reduce hologram resolution deterioration.

Moreover, we implemented a simpler single-SLM system to generate multiplane CGH light patterns without TF. Although this approach has a poor axial connement with respect to 3D-CGH-TF, it could still be a powerful method for applications using sparsely expressing samples. Previous schemes for multi-plane generation of laterally extended shapes through high-NA objectives reported patterns spanning a limited (tens of microns) axial range27,30,31 and used linear excitation. Here we exploited the whole excitation volume reachable with the LCOS-SLM using high-NA objectives. This was achieved by minimizing spherical aberration of targets generated out of the nominal objective focal plane by accounting for the high-NA objective in the expression of the Fresnel lens-phase prole. In addition, throughout the accessible volume, we obtained uniform light distribution among the target shapes. Weighting the amplitude of each target input to the GS algorithm according to position enabled compensation of intensity loss because of diffraction efciency, which decreases with increasing distance from the excitation volume centre.

We demonstrated the capabilities of the 3D-CGH-TF system by performing multiplane 2P photoconversion of Kaede-expressing neurons. Previous in vivo experiments photo-converting Kaede in single neurons of zebrash larvae used one-photon visible light illumination43,56,57. Here we demonstrated 2P-mediated single-cell resolution simultaneous photoconversion of multiple Kaede-expressing neurons. The unprecedented spatial precision achieved by our approach can, in the future, be extended to track the morphology of single neurons across the entire nervous system of the zebrash larva.

For multiplane photoconversion of multiple isolated targets (for instance, in sparsely expressing samples), 3D-CGH holds an advantage over 3D-CGH-TF by enabling light shaping in a greater number of planes. However, 3D-CGH-TF enables a better axial connement, a key parameter for experiments requiring illumination of spatially nearby multiple targets or large areas. Moreover, generation of multiplane illumination patterns with 3D-CGH also produces spurious light in the intermediate planes absent in illumination geometries produced with 3D-CGH-TF. In addition, we have previously demonstrated that TF is particularly robust to scattering making 3D-CGH-TF more suitable than 3D-CGH for applications requiring in-depth illumination through scattering media23,24.

The same optical design could be used to photoswitch other proteins such as photactivatable green uorescent protein58 or kindling uorescent protein (KFP1)59 to precisely track the 3D position of specic cells in vivo during embryo development. Combined with optogenetics (Supplementary Note 3) or uncaging, multiplane generation of spatiotemporally focused patterns will enable simultaneous control of neurons and substructures in different planes, as well as provide a exible mean to stimulate locations lying above or below the imaging plane. Combined with extended depth-of-eld imaging60, multiplane light patterning could also improve the spatial specicity of functional voltage or calcium physiology by shaping light on cells or structures of interest61.

In this study, the two-SLM system is combined with HiLo imaging. Combining this approach with 2P imaging would also enable deep photostimulation and imaging in vivo. Decoupling of lateral and axial wavefront shaping could also be adopted in optical designs different from the one presented here. For example, placing a second SLM at the Fourier plane of a fast switchable array should enable 3D-encoded multisite 2P microscopy62 or high-speed 3D holographic light patterning63.

Methods

Two-SLM optical set-up. The optical system, schematically depicted in Fig. 1 and Supplementary Fig. 11, was built around a commercial Olympus IX71 inverted microscope, modied in order to accommodate two opposite-facing objectives, OBJ1 and OBJ2, for excitation and uorescence collection, respectively. To this end, the condenser lens of the microscope was substituted with a dielectric mirror and an Olympus LUMPLFL60xW/IR2, NA 0.90 objective (OBJ1). The expanded ( 10) beam of a Ti:Sapphire laser (MaiTai Deep-See, Spectra-Physics) covered the

active area of a rst LCOS-SLM (X10468-07, Hamamatsu Photonics; SLM1), which modulated the phase of the incoming beam to create a rst image of the desired intensity pattern on the diffraction grating (830 l/mm, 53004ZD02-035R, Richardson Gratings; G) for TF through the lens L1 (f1 500 mm). 2D-phase

holograms were calculated using a standard GS algorithm20,33. The rst diffraction order was subsequently collimated by lens L2 (f2 500 mm) and impinged on a

second SLM (X10468-07, Hamamatsu Photonics; SLM2), which was imaged at the back focal plane of the excitation objective, OBJ1, via a 2:1 telescope (lenses L3, f3 1,000 mm and L4, f4 500 mm). Suppression of the zero-order spot arising

from SLM1 was achieved by using two cylindrical lenses (fL1 1,000 mm and

fL2 1,000 mm) oriented at 45 and 45 with respect to the grating

lines64.

Holographic light patterns generated at the sample volume illuminated a thin spin-coated uorescent layer of rhodamine-6G in polymethyl methacrylate 2% w/v in chloroform and the induced uorescence was imaged on a CCD camera (CoolSNAP HQ2, Roper Scientic) through OBJ2 (Olympus UPLSAPO60XW, NA 1.2). For 3D reconstruction of illumination volumes, OBJ2 was xed and focused on the uorescent layer, while OBJ1 was moved along the axial direction with a piezo positioner of 1 mm range when working in closed loop (PI N-725.2A PIFOC). The two SLMs, the CCD camera, the piezo positioner, lasers and other electronic components of the set-up were controlled by a custom-developed interface in LabVIEW. GS-based algorithms were run in MATLAB.

When the set-up was used for generation of multiplane holographic patterns (not temporally focused) a mirror replaced the diffraction grating, and SLM2 was either used in reectance mode by applying only the atness correction phase mask of the device or was replaced by a mirror. Both the grating and SLM2 were mounted on magnetic bases enabling fast switching between the different congurations. The multiplane GS algorithm used in this case was run in a custom-designed C software interface, Wavefront Designer20.

8 NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928 ARTICLE

Two-colour HiLo imaging system. High-resolution multiplane uorescence imaging of zebrash larvae was achieved by coupling an optical set-up for two-colour HiLo microscopy46 to the Olympus IX71 microscope. Two continuous wave 473-nm (Laser Quantum, Ciel 350 mW) and 561-nm (CNI laser, MGL-N-561-500 mW) lasers were co-aligned in the same optical path with the dichroic mirror D2 (Semrock, Di02-R514) and collimated with lenses L6(f6 125 mm), L7 (f7 150 mm) and L10 (f10 35 mm) to illuminate an oscillating

diffuser plate (Optotune LSR-3005-10) that was imaged through lenses L8(f8 75 mm) and L9 (f9 200 mm) at the back aperture of the excitation objective,

OBJ1 (Supplementary Fig. 11). The D1 dichroic mirror reected the collected uorescence to a CMOS camera (Hamamatsu Photonics, Orca Flash 4.0-V2) through the appropriate lter cube (FC2) for green (dichroic mirror Semrock FF495-Di02, emission lter Semrock FF01-520/35-25) or red uorescence (dichroic mirror Semrock Di02-R561, emission lter Semrock FF595-Di02).

The sectioned image was computed with custom scripts written in MATLAB65. The cutoff frequency used to merge the low- and high-frequency components was chosen such as kcE0.1 klow, where klow is the frequency of the low-pass lter applied to the uniform illumination image. With those parameters, we measured an axial resolution of 3.2-mm FWHM for the emitted uorescence (Supplementary

Fig. 11b). The axial resolution was measured using a Rhodamine-6G thin layer. First, we recorded a z-stack with uniform illumination (oscillating diffuser on) and then with speckle illumination (oscillating diffuser off). The two stacks were then processed in MATLAB to generate an axially resolved HiLo z-stack using algorithms previously described34,65. The axial resolution shown in Supplementary Fig. 11 is the axial resolution measured on the processed HiLo z-stack.

2P galvo-based scanning imaging system. 2P imaging of photoconverted zebrash larvae in the brain performed by a mode-locked Ti-Sapphire laser source (Coherent Chameleon Vision II, pulse width 140 fs, tuning range 6801,080 nm). The femtosecond pulsed beam was raster-scanned on the sample via a pair of xy galvanometric mirrors (3 mm aperture, 6215H series, Cambridge Technology) imaged at the back aperture of the microscope objective ( 40 W APO NIR,

Nikon) through an afocal telescope (scan lens: f 100 mm, tube lens: f 300 mm).

Galvanometric mirrors were driven by two servo drivers (MicroMax series 671, Cambridge Technology) controlled by a Digital/Analog converter board (PCI-6110, National Instrument). Emitted uorescence was collected by a bre-coupled detection scheme66. The bre exit was imaged on two photomultiplier tubes GaAsP (H10770-40 SEL, Hamamatsu Photonics, active area 5 mm) by a set of three matching asphere lenses (f 23.5 mm, Melles Griot #LAG-32.5-23.5-C).

Following the bre exit, uorescence light was ltered with an infrared-light-blocking lter (FF01-750sp, Semrock), split into two channels by a dichroic mirror (FF555-Di03, Semrock) and detected through two emission lters (FF01-510/84 and FF02-617/73, Semrock). The whole system was built around a commercial upright microscope (SliceScope, Scientica). 2P imaging laser power was tuned by combining an electrically controlled liquid crystal variable phase retarder (LRC-200-IR1, Meadowlark Optics) and a polarizer cube (BB-050-IR1, Meadowlark Optics) at the exit of the laser source.

Green and red uorescence z-stacks of photoconverted Kaede in the zebrash brain were acquired by scanning the excitation beam (780 nm) at 0.74 Hz (full frame) and averaging 1020 frames for each plane.

Photoconversion protocol. First, a HiLo z-stack in the green channel(200 200 100 mm3) was recorded to map the location of neuronal cells for

photoconversion. On the basis of these images, we calculated phase holograms that produced the corresponding 2D or 3D illumination patterns. We typically used 5-mm-diameter holographic spots to target single cells and 3035-mm-diameter holographic spots to target sets of neurons. In order to quantify the efciency of photoconversion, we also recorded the corresponding z-stack in the red channel before photoconversion.

Simultaneous 2P photoconversion (lphot 800 nm) of all targets was performed

while monitoring the uorescence in the red channel. We typically observed a tenfold increase of red uorescence in the targeted cells. Photoconversion during uorescence imaging was minimized by keeping the total acquisition time below 2 min and laser power at the sample plane below 20 mW. To minimize thermal damage during photoconversion, we delivered trains of 50-ms pulses, low laser intensity B0.044.0 mW mm 2 (power densities are always given relatively to the area of the spots surface) for periods of time that ranged from 200 ms to a few hundred seconds depending on the laser intensity.

Transgenic lines. Experiments were performed on Danio rerio larvae between 2 and 6 days post fertilization following procedures approved by the Institutional Ethics Committee Darwin in the Institut du Cerveau et de la Moelle pinire (ICM). AB and TL strains of wild-type (WT) larvae were obtained from laboratorys stock of adults. Embryos and larvae were raised in an incubator at28.5 C until shortly before recordings were performed. For photoconversion experiments, we used Tg(HuC:gal4; UAS:kaede)43 where the HuC promoter drives pan-neuronal expression of Gal4 and Kaede at the larval stage. Tg(pkd2l1:gal4; UAS:ChR2-H134R-mCherry; UAS:GCaMP5G)65,67 were used for combination of optogenetics and calcium imaging in Supplementary Data. Before performing

image acquisitions, embryos were dechorionated and screened for uorescence at 1 days post fertilization. Larvae screened for Kaede uorescence were later embedded laterally in 1.5% agarose. Larvae were anaesthetized in 0.02% tricain (MS-222, Sigma-Aldrich, USA).

Data availability. The data that support the ndings of this study are available from the corresponding author upon request.

References

1. Denk, W., Strickler, J. H. & Webb, W. W. Two-photon laser scanning uorescence microscopy. Science 248, 7376 (1990).

2. Salom, R. et al. Ultrafast random-access scanning in two-photon microscopy using acousto-optic deectors. J. Neurosci. Methods 154, 161174 (2006).3. Reddy, G. D., Kelleher, K., Fink, R. & Saggau, P. Three-dimensional random access multiphoton microscopy for functional imaging of neuronal activity. Nat. Neurosci. 11, 713720 (2008).

4. Grewe, B. F., Langer, D., Kasper, H., Kampa, B. M. & Helmchen, F. High-speed in vivo calcium imaging reveals neuronal network activity with near-millisecond precision. Nat. Methods 7, 399405 (2010).

5. Nguyen, Q. T., Callamaras, N., Hsieh, C. & Parker, I. Construction of a two-photon microscope for video-rate Ca2 imaging. Cell Calcium 30, 383393 (2001).

6. Oron, D., Tal, E. & Silberberg, Y. Scanningless depth-resolved microscopy. Opt. Express 13, 14681476 (2005).

7. Zhu, G., van Howe, J., Durst, M., Zipfel, W. & Xu, C. Simultaneous spatial and temporal focusing of femtosecond pulses. Opt. Express 13, 21532159 (2005).

8. Akturk, S., Gu, X., Bowlan, P. & Trebino, R. Spatio-temporal couplings in ultrashort laser pulses. J. Opt. 12, 093001 (2010).

9. Tal, E., Oron, D. & Silberberg, Y. Improved depth resolution in video-rate line-scanning multiphoton microscopy using temporal focusing. Opt. Lett. 30, 16861688 (2005).

10. Vaziri, A. & Shank, C. V. Ultrafast wideeld optical sectioning microscopy by multifocal temporal focusing. Opt. Express 18, 1964519655 (2010).11. Isobe, K. et al. Enhancement of lateral resolution and optical sectioning capability of two-photon uorescence microscopy by combining temporal-focusing with structured illumination. Biomed. Opt. Express 4, 23962410 (2013).

12. Choi, H. et al. Improvement of axial resolution and contrast in temporally focused wideeld two-photon microscopy with structured light illumination. Biomed. Opt. Express 4, 9951005 (2013).

13. Song, Q. et al. Two-dimensional spatiotemporal focusing of femtosecond pulses and its applications in microscopy. Rev. Sci. Instrum. 86, 083701 (2015).

14. Losonczy, A., Zemelman, B. V, Vaziri, A. & Magee, J. C. Network mechanisms of theta related neuronal activity in hippocampal CA1 pyramidal neurons. Nat. Neurosci. 13, 967972 (2010).

15. Dana, H. et al. Hybrid multiphoton volumetric functional imaging of large-scale bioengineered neuronal networks. Nat. Commun. 5, 3997 (2014).

16. Vaziri, A., Tang, J., Shroff, H. & Shank, C. V. Multilayer three-dimensional super resolution imaging of thick biological samples. Proc. Natl Acad. Sci. USA 105, 2022120226 (2008).

17. Kim, D. & So, P. T. C. High-throughput three-dimensional lithographic microfabrication. Opt. Lett. 35, 16021604 (2010).

18. Andrasfalvy, B. K., Zemelman, B. V, Tang, J. & Vaziri, A. Two-photon single-cell optogenetic control of neuronal activity by sculpted light. Proc. Natl. Acad. Sci. USA 107, 1198111986 (2010).

19. Rickgauer, J. P., Deisseroth, K. & Tank, D. W. Simultaneous cellular-resolution optical perturbation and imaging of place cell ring elds. Nat. Neurosci. 17, 18161824 (2014).

20. Lutz, C. et al. Holographic photolysis of caged neurotransmitters. Nat. Methods 5, 821827 (2008).

21. Papagiakoumou, E., de Sars, V., Oron, D. & Emiliani, V. Patterned two-photon illumination by spatiotemporal shaping of ultrashort pulses. Opt. Express 16, 2203922047 (2008).

22. Papagiakoumou, E. et al. Scanless two-photon excitation of channelrhodopsin-2. Nat. Methods 7, 848854 (2010).

23. Bgue, A. et al. Two-photon excitation in scattering media by spatiotemporally shaped beams and their application in optogenetic stimulation. Biomed. Opt. Express 4, 28692879 (2013).

24. Papagiakoumou, E. et al. Functional patterned multiphoton excitation deep inside scattering tissue. Nat. Photon. 7, 274278 (2013).

25. Piestun, R., Spekor, B. & Shamir, J. Wave elds in three dimensions: analysis and synthesis. J. Opt. Soc. Am. A 13, 18371848 (1996).

26. Curtis, J. E., Koss, B. A. & Grier, D. G. Dynamic holographic optical tweezers. Opt. Commun. 207, 169175 (2002).

27. Sinclair, G. et al. Interactive application in holographic optical tweezers of a multi-plane Gerchberg-Saxton algorithm for three-dimensional light shaping. Opt. Express 12, 16651670 (2004).

NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms11928

28. Yang, S. et al. Three-dimensional holographic photostimulation of the dendritic arbor. J. Neural Eng. 8, 46002 (2011).

29. Anselmi, F., Ventalon, C., Bgue, A., Ogden, D. & Emiliani, V. Three-dimensional imaging and photostimulation by remote-focusing and holographic light patterning. Proc. Natl Acad. Sci. USA 108, 1950419509 (2011).

30. Hilario, P. L., Villangca, M. J. & Tapang, G. Independent light elds generated using a phase-only spatial light modulator. Opt. Lett. 39, 20362039 (2014).

31. Ma, B. et al. Generation of three-dimensional optical structures by dynamic holograms displayed on a twisted nematic liquid crystal display. Appl. Phys. B 110, 531537 (2013).

32. Emiliani, V. et al. Wave front engineering for living cells microscopy. Opt. Express 13, 13951405 (2005).

33. Gerchberg, R. W. & Saxton, W. O. A pratical algorithm for the determination of the phase from image and diffraction pictures. Optik (Stuttg) 35, 237246 (1972).

34. Zahid, M. et al. Holographic photolysis for multiple cell stimulation in mouse hippocampal slices. PLoS ONE 5, e9431 (2010).

35. Dana, H. & Shoham, S. Remotely scanned multiphoton temporal focusing by axial grism scanning. Opt. Lett. 37, 29132915 (2012).

36. Richards, B. & Wolf, E. Electromagnetic diffraction in optical systems. II. Structure of the image eld in an aplanatic system. Proc. R Soc. Lond. A 253, 358379 (1959).

37. Botcherby, E. J., Juskaitis, R., Booth, M. J. & Wilson, T. Aberration-free optical refocusing in high numerical aperture microscopy. Opt. Lett. 32, 20072009 (2007).

38. Leshem, B., Hernandez, O., Papagiakoumou, E., Emiliani, V. & Oron, D. When can temporally focused excitation be axially shifted by dispersion? Opt. Express 22, 70877098 (2014).

39. Dainty, J. C. Laser Speckle and Related Phenomena. Topics in Applied Physics 9 (Springer-Verlag, 1975).

40. Matar, S., Golan, L. & Shoham, S. Reduction of two-photon holographic speckle using shift-averaging. Opt. Express 19, 2589125899 (2011).

41. Haist, T., Schnleber, M. & Tiziani, H. Computer-generated holograms from 3D-objects written on twisted-nematic liquid crystal displays. Opt. Commun. 140, 299308 (1997).

42. Leach, J. et al. Interactive approach to optical tweezers control. Appl. Opt. 45, 897903 (2006).

43. Scott, E. K. et al. Targeting neural circuitry in zebrash using GAL4 enhancer trapping. Nat. Methods 4, 323326 (2007).

44. Ando, R., Hama, H., Yamamoto-Hino, M., Mizuno, H. & Miyawaki, A. An optical marker based on the UV-induced green-to-red photoconversion of a uorescent protein. Proc. Natl Acad. Sci. USA 99, 1265112656 (2002).

45. Isobe, K. et al. Measurement of two-photon excitation spectrum used to photoconvert a uorescent protein (Kaede) by nonlinear Fourier-transform spectroscopy. Biomed. Opt. Express 1, 687693 (2010).

46. Lim, D., Chu, K. K. & Mertz, J. Wide-eld uorescence sectioning with hybrid speckle and uniform-illumination microscopy. Opt. Lett. 33, 18191821 (2008).

47. Carslaw, H. S. & Jaeger, J. C. Conduction of Heat in Solids (Clarendon Press, 1959).

48. Elwassif, M. M., Kong, Q., Vazquez, M. & Bikson, M. Bio-heat transfer model of deep brain stimulation induced temperature changes. Conf. Proc. IEEE Eng. Med. Biol. Soc. 1, 35803583 (2006).

49. Johansson, J. D. Spectroscopic method for determination of the absorption coefcient in brain tissue. J. Biomed. Opt. 15, 057005 (2010).

50. Yizhar, O., Fenno, L. E., Davidson, T. J., Mogri, M. & Deisseroth, K. Optogenetics in neural systems. Neuron 71, 934 (2011).

51. Grewe, B. F., Voigt, F. F., van t Hoff, M. & Helmchen, F. Fast two-layer two-photon imaging of neuronal cell populations using an electrically tunable lens. Biomed. Opt. Express 2, 20352046 (2011).

52. Durst, M. E., Zhu, G. & Xu, C. Simultaneous spatial and temporal focusing for axial scanning. Opt. Express 14, 1224312254 (2006).

53. Durst, M. E., Zhu, G. & Xu, C. Simultaneous spatial and temporal focusing in nonlinear microscopy. Opt. Commun. 281, 17961805 (2008).

54. Suchowski, H., Oron, D., Silberberg, Y., Suchowski Oron, D. & Silberberg, Y.,H. Generation of a dark nonlinear focus by spatio-temporal coherent control. Opt. Commun. 264, 482487 (2006).55. Du, R. et al. Analysis of fast axial scanning scheme using temporal focusing with acousto-optic deectors. J. Mod. Opt. 56, 8184 (2009).

56. Hatta, K., Tsujii, H. & Omura, T. Cell tracking using a photoconvertible uorescent protein. Nat. Protoc. 1, 960967 (2006).

57. Warp, E. et al. Emergence of patterned activity in the developing zebrash spinal cord. Curr. Biol. 22, 93102 (2012).

58. Patterson, G. H. A photoactivatable GFP for selective photolabeling of proteins and cells. Science 297, 18731877 (2002).

59. Chudakov, D. M. et al. Kindling uorescent proteins for precise in vivo photolabeling. Nat. Biotechnol. 21, 191194 (2003).

60. Quirin, S., Peterka, D. S. & Yuste, R. Instantaneous three-dimensional sensing using spatial light modulator illumination with extended depth of eld imaging. Opt. Express 21, 16007 (2013).

61. Foust, A. J., Zampini, V., Tanese, D., Papagiakoumou, E. & Emiliani, V. Computer-generated holography enhances voltage dye uorescence discrimination in adjacent neuronal structures. Neurophotonics 2, 021007 (2015).

62. Ducros, M., Goulam Houssen, Y., Bradley, J., de Sars, V. & Charpak, S. Encoded multisite two-photon microscopy. Proc. Natl Acad. Sci. USA 110, 1313813143 (2013).

63. Go, M. A. et al. Four-dimensional multi-site photolysis of caged neurotransmitters. Front. Cell. Neurosci. 7, 231 (2013).

64. Hernandez, O., Guillon, M., Papagiakoumou, E. & Emiliani, V. Zero-order suppression for two-photon holographic excitation. Opt. Lett. 39, 59535956 (2014).

65. Lauterbach, M. A., Ronzitti, E., Sternberg, J. R., Wyart, C. & Emiliani, V. Fast calcium imaging with optical sectioning via HiLo microscopy. PLoS ONE 10, e0143681 (2015).

66. Ducros, M. et al. Efcient large core ber-based detection for multi-channel two-photon uorescence microscopy and spectral unmixing. J. Neurosci. Methods 198, 172180 (2011).

67. Fidelin, K. et al. State-dependent modulation of locomotion by GABAergic spinal sensory neurons. Curr. Biol. 25, 30353047 (2015).

Acknowledgements

We thank Jenna Sternberg for help with the preparation of the biological samples, Marc Guillon for helpful discussion on the implementation of the multiplane GS algorithm, Marcel Lauterbach and Emiliano Ronzitti for discussions on the design of the HiLo set-up, Emiliano Ronzitti and Amanda Foust for building the 2P imaging system, Benot Forget and Alexis Picot for useful discussions on heating effects,Vincent de Sars for the continued development of Wavefront Designer softwareand Amanda Foust for critical reading of the manuscript. V.E. and E.P. acknowledge the Agence Nationale de la Recherche (grants ANR-12-BSV5-0011-01, Neurholog and ANR-15-CE19-0001-01, 3DHoloPAc), VE acknowledges the National Institutes of Health (NIH 1-U01-NS090501-01). O.H. acknowledges the programme Nanotechnologies France-Israel for nancial support. C.W. and K.F. acknowledge the City of Paris Emergence programme, the Human Frontier Science Program (HFSP) Research Grant #RGP0063/2014 and the European Research Council (ERC) starter grant OptoLoco #311673.

Author contributions

O.H., E.P. and V.E. designed the experiments. O.H. built up the optical system, performed the experiments and analysed the data. E.P. participated in the development of the optical system and experiments. D.T. participated in 2P imaging experiments of photoconverted Kaede. K.F. prepared zebrash samples. C.W. helped in designing the zebrash experiments. E.P. and V.E. wrote the manuscript, with contributions from O.H. and C.W. V.E. supervised the project.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions

Web End =reprintsandpermissions .

How to cite this article: Hernandez, O. et al. Three-dimensional spatiotemporal focusing of holographic patterns. Nat. Commun. 7:11928 doi: 10.1038/ncomms11928 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

10 NATURE COMMUNICATIONS | 7:11928 | DOI: 10.1038/ncomms11928 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications