Axial plane single-molecule super-resolution microscopy of whole cells

Fluorescence nanoscopy has become an indispensable tool for studying organelle structures, protein dynamics, and interactions in biological sciences. Single-molecule localization microscopy can now routinely achieve 10–50 nm resolution through fluorescently labeled specimens in lateral optical sections. However, visualizing structures organized along the axial direction demands scanning and imaging each of the lateral imaging planes with fine intervals throughout the whole cell. This iterative process suffers from photobleaching of tagged probes, is susceptible to alignment artifacts and also limits the imaging speed. Here, we focused on the axial plane super-resolution imaging which integrated the single-objective light-sheet illumination and axial plane optical imaging with single-molecule localization technique to resolve nanoscale cellular architectures along the axial (or depth) dimension without scanning. We demonstrated that this method is compatible with DNA points accumulation for imaging in nanoscale topography (DNA-PAINT) and exchange-PAINT by virtue of its light-sheet illumination, allowing multiplexed super-resolution imaging throughout the depth of whole cells. We further demonstrated this proposed system by resolving the axial distributions of intracellular organelles such as microtubules, mitochondria, and nuclear pore complexes in both COS-7 cells and glioblastoma patient-derived tumor cells. © 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

While 2D SMLM allows capturing super-resolution images in the lateral direction, cellular structures organized along the axial direction (e.g. cytokinetic apparatus or nuclear envelope) mandate the lateral imaging plane to be scanned throughout the cell or tissue specimens, a time-consuming process. In addition, alignment procedures of the reconstructed super-resolution volumes are prone to misalignment and image artifacts especially in presence of specimen-induced aberrations [19]. Although the conventional light-sheet fluorescence microscopy (LSFM) [20][21][22] can directly obtain the axial plane information by using two orthogonal objectives, the system limits the use of objectives with high numerical apertures (NA) which is important for SMLM, and restricts the sample mounting strategy, making it difficult to study cells prepared on regular coverslips.
Single high-NA objective lens (e.g. NA = 1.4, as used in this study) is capable of generating high-quality light sheet in the axial plane. However, only the lateral plane can be captured in conventional SMLM experiments and therefore unable to take direct advantage of the selfgenerated light sheet [23,24]. By designing a specialized microfluidic sample mounting system with a mirror surface, soSPIM [25] and SO-LSM [26] converted the axial light sheet to the lateral direction allowing SMLM with single objective lens. In addition, by using a prism reflected and tilted light sheet generated by a separate objective lens, TILT3D [27] enabled light-sheet illumination in SMLM without the specialized microfluidic-mirror system. However, these developments coupled the light sheet position with one or more sample translation directions resulting in complicated optical readjustments when translating a specimen from one location to another. It is worth mentioning that, Li et al. provided a novel technique, termed axial plane optical microscopy (APOM) [28], which used a remote objective lens and a 45°tilted mirror to convert the axial information to the lateral plane and re-image onto a camera. In this case, the axial plane information could be directly imaged by single shot without scanning. This fast and high-contrast imaging approach is particularly suitable for studying thick samples and live cells.
In this work, we developed an axial plane SMLM integrating the single-objective light-sheet illumination and axial plane optical imaging with single-molecule super-resolution microscopy to resolve nanoscale cellular architecture along the axial (or depth) dimension. This system has a similar optical design as the recently published oblique-plane single-molecule localization microscopy (obSTORM) [29] but focusing on the axial plane super-resolution imaging utilizing two high-NA objective lenses (NA = 1.4). As one of the unique advantages of the axial plane SMLM system, the imaging plane overlaps with the illumination plane formed by a highly confined light sheet generated by the same objective lens for collecting fluorescence emission. This makes the system compatible with DNA points accumulation for imaging in nanoscale topography (DNA-PAINT) as well as its multiplexed version exchange-PAINT [30][31][32][33][34][35] allowing multiplexed super-resolution imaging throughout the depth of a specimen. In this case, utilizing the same type of fluorophore conjugated to orthogonal DNA strands for sequential axial plane multiplexed imaging requires a single laser source for excitation and is free of chromatic aberrations and distortions using a single emission channel. We demonstrated the applications of the presented system by resolving the axial distributions of cellular organelles and macromolecular assemblies such as microtubules, mitochondria, and nuclear pore complexes in both COS-7 cells and glioblastoma patient-derived tumor cells (GBM10).

The optical setup of the axial plane SMLM
The complete optical diagram of the axial plane SMLM is shown in Fig. 1(A).
In the excitation part, three lasers of 642 nm (2RU-VFL-P-2000-642-B1R, MPB Communications Inc.), 488 nm (Cobolt 06-MLD 488 nm 200 mW, Cobolt Inc.) and 405 nm (DL405-100, CrystaLaser) wavelengths were combined by two dichroic mirrors (Di1, LM01-503-25; Di2, FF01-430/LP-25, Semrock), and then modulated by an acousto-optic tunable filter (AOTF, AOTFnC-400.650-TN, AA OPTO-ELECTRONIC). Thereafter, all the lasers were coupled into a single-mode polarization-maintaining fiber (P1-405BPM-FC, Thorlabs). An achromatic lens (L5, f = 20 mm, AC080-020-A-ML, Thorlabs) was used for collimating the output laser beam. A light sheet was generated by a pair of cylindrical lenses (CL1, f = 200 nm and CL2, f = 20 nm), and its thickness was further confined by placing a slit in the focal plane of CL2. After passing through L6 (f = 175 nm) and an oil-immersion objective lens (Obj1, UPLSAPO 100XO, NA = 1.4, OLYMPUS), hereafter referred to as imaging objective, a quad-band dichroic mirror (Di3, Di03-R405/488/561/635-t3-25X36, Semrock) was placed between these two components to separate the excitation and emission light, a vertical light sheet with ∼ 1.19 µm thickness was created on the sample position for illuminating a certain axial slice. Epi-illumination mode can be switched from light-sheet illumination mode by removing the light-sheet creation part marked with the blue dashed rectangle. Samples were mounted on a custom-made stage with four degrees of freedom in x, y, z and rotation about z. The translations were carried out by stepper motors and piezo stages (MAX341/M, Thorlabs), while a rotation stage (U-651, PI) was assembled on top of the translation stages for controlling the rotation of the sample to allow imaging a particular vertical section of a specimen using a stationary light sheet beam.
In the emission part, the fluorescence signal was collected by the imaging objective Obj1, transmitted through Di3, and relayed to a polarizing beam splitter (PBS252, Thorlabs) by a 4f system with two achromatic lenses L7 and L8 (f = 175 nm). The polarizing beam splitter divided the emission light into two parts. The reflected part was used to image the lateral plane information on a complementary metal-oxide-semiconductor (CMOS) camera (UI-3060CPM-GL, IDS Imaging Development Systems GmbH). The transmitted part passed through a quarter-wave plate (QWP, AHWP05M-600, Thorlabs) and formed an intermediate image at the focal point of another objective lens (Obj2, identical with Obj1), hereafter referred to as remote objective. Then, a 45°tilted mirror was used to convert the axial plane image into the lateral plane with its surface coinciding on the focal point of Obj2. Since the working distance of the high-NA objective was as small as 0.13 mm, only the edge of the 45°tilted mirror was used. Therefore, the mirror needs a high-quality reflecting surface as well as a straight edge. Here, the mirror was made using a cleaved silicon wafer coated by a thin layer of aluminum (150 nm), resulting in a straight and sharp edge. This wafer was glued on a right-angle prism for setting the 45°angle. Then, we fixed the 45°tilted mirror on a six-axis kinematic mount (K6XS, Thorlabs), which was placed on a translation stage (M-460A, Newport) controlled by a picomotor actuator (8303, Newport). With this configuration, the backward emission light, reflected by the 45°tilted mirror, was collected by the remote objective, and passed through the QWP again. The polarizing beam splitter and the QWP were combined to allow the p-polarized fluorescence signal to be detected by a scientific complementary metal-oxide-semiconductor (sCMOS) camera with high sensitivity and low noise (Orca-Flash4.0LT, Hamamatsu). Two achromatic lenses L11 (f = 200 mm) and L12 (f = 100 mm) were used for adjusting the magnification to map an effective pixel size of 130 nm in the specimen [36]. An excitation cleanup filter (CF, FF01-390/482/563/640-25, Semrock) and a matched emission filter (EF, FF01-446/523/600/677-25, Semrock) were used in this system. In addition, a white LED was placed on top of the Obj1 for imaging in the transmission mode.

Sample preparation
Single layer beads sample preparation: A cleaned 25-mm-diameter round coverslip was incubated with 200 µL Poly-L-Lysine (P4707, Sigma-Aldrich) for 30 min. Then the coverslip was rinsed thrice with deionized water, dried and incubated with 200 µL 100 nm-diameter crimson beads (custom-designed, Invitrogen) diluted to 1:10 5 in Phosphate Buffered Saline (PBS) (10010023, Gibco) for 30 min. After rinsing thrice with deionized water, the coverslip was mounted on a custom-made sample holder, and 20 µL 97% TDE (166782, Sigma-Aldrich) was added on top of the coverslip. Another pre-cleaned coverslip was carefully placed on top of it to extend the buffer over the whole coverslip. The excessive buffer was removed with tissue papers. Finally, the sample was sealed with the two-component silicon dental glue (picodent twinsil speed 22, Dental-Produktions und Vertriebs GmbH, Wipperfürth).
Agarose embedded beads sample preparation: 8 mg of agarose (A9045-5G, Sigma-Aldrich) was added to 400 µL beads diluted to 1:10 5 in PBS. The mixture was heated to 69°C for 5-10 min and vortexed to blend evenly. A cleaned coverslip was mounted on a custom-made sample holder. Then, 400 µL prepared mixture was added on top of the coverslip. Another pre-cleaned coverslip was immediately placed on top of it to extend the mixture over the whole coverslip. Then the sample was put in refrigerator to solidify. Finally, the sample was sealed with the two-component silicon dental glue.
Immunofluorescence labeling: Mammalian cells (COS-7 or GBM10) were cultured at 37°C with 5% CO 2 on autoclaved coverslips. The culture media contained 10% FBS (30-2020, ATCC), 1% PS (15140122, Gibco) in DMEM (30-2002, ATCC). Cells were fixed and immune-labeled when the confluence reached ∼80%. Mitochondria or microtubules labeling: COS-7 cells were fixed for 15 min at room temperature (RT) with 37°C pre-warmed 3% Paraformaldehyde (PFA, 15710, Electron Microscopy Sciences) and 0.5% Glutaraldehyde (GA, 16019, Electron Microscopy Sciences) in PBS, then rinsed twice with PBS and treated for 7 min with freshly prepared fluorescence quenching buffer containing 0.1% NaBH 4 (452882, Sigma-Aldrich) in PBS. After rinsing thrice with PBS, cells were permeabilized and blocked for 30 min with 3% BSA (001-000-162, Jackson ImmunoResearch) and 0.2% Triton X-100 (X100, Sigma-Aldrich) in PBS, gently rocked at RT. After blocking, cells were incubated with the primary antibody diluted to 1:500 in 1% BSA and 0.2% Triton X-100 in PBS, overnight at 4°C. The primary antibodies for mitochondria and microtubules were anti-Tom20 (Rabbit polyclonal, sc-11415, Santa Cruz Biotechnology) and anti-α-tubulin (T5168, Sigma-Aldrich), respectively. Cells were subsequently rinsed thrice for 5 min each time with wash buffer (0.05% Triton X-100 in PBS), and then incubated with the secondary antibody diluted to 1:500 in 1% BSA and 0.2% Triton X-100 in PBS for 4 h. The secondary antibodies for mitochondria and microtubules were AF647-anti-rabbit (A21245, Molecular Probes) and AF647-anti-mouse (A21236, Molecular Probes), respectively. Then cells were rinsed thrice for 5 min each time with wash buffer and post-fixed with 4% PFA in PBS for 10 min at RT. After rinsing thrice with PBS, cells were stored in PBS at 4°C. DNA-PAINT labeling: Only the secondary antibody was replaced to DNA-conjugated anti-mouse P1 and/or anti-rabbit P4 according to the primary antibody, and the other labeling operations followed the above-mentioned protocol. Nuclear pore complex labeling: Cells were rinsed for 20 s with 37°C pre-warmed 2.4% PFA in PBS, then extracted with 37°C pre-warmed 0.4% Triton X-100 in PBS for 3 min. After extraction, cells were fixed with 37°C pre-warmed 2.4% PFA in PBS for 30 min at RT, then rinsed thrice with PBS and treated for 7 min with freshly prepared fluorescence quenching buffer (0.1% NaBH 4 in PBS). After rinsing thrice with PBS, cells were permeabilized and blocked for 1 h with 5% BSA in PBS, gently rocked at RT. Then, after blocking, cells were sequentially incubated with the primary antibody Nup98 (2598S, Cell Signaling Technology) and secondary antibody AF647-anti-Rabbit (A21245, Molecular Probes), the labeling operations followed the above-mentioned protocol.

Data acquisition
Point spread function (PSF) acquisition using fluorescence beads: A prepared single layer fluorescence beads sample was placed on the sample stage controlled by piezo stages. After aligning the 45°tilted mirror, the axial plane image of a bead was identified on the camera. Then the PSF data was taken by scanning the sample along the y-axis with 100 nm step size. The power density of the illumination light sheet was ∼0.3-0.6 kW/cm 2 .
Single-color SMLM imaging data acquisition: An agarose embedded fluorescence beads sample was used for aligning the 45°tilted mirror. The beads sample was excited by the 642 nm laser with an intensity of 0.3-0.6 kW/cm 2 . We optimized the imaging system by adjusting the tilted mirror such that the illumination light sheet plane overlapped with the imaged axial plane on the camera. After optimization, the beads sample was replaced with the cell specimen. Before SMLM data acquisition, the specimen was excited at a high laser power to push most fluorophores into the dark state therefore reduce the background. The SMLM imaging was performed at 20-50 frames/s. The total number of frames per dataset ranged from 10,000 to 30,000 frames. The power density of the light-sheet illumination was ∼2-5 kW/cm 2 .
Two-channel SMLM imaging data acquisition: After aligning the 45°tilted mirror following the above protocol, cell specimen was placed on the sample stage. A wide-field reference image under transmission light was recorded for two-channel alignment later. Then, the blinking dataset of the first channel was recorded at 20 frames/s for 20,000-30,000 frames. After data acquisition of the first channel was finished, the first imaging buffer was removed, the sample was washed thrice with PBS, and then the second imaging buffer was infused. The buffer exchange was accomplished by using a custom-made fluid exchange tool. This tool included two 1-meter-long thin plastic tubes connected with needles on all ends. The needle on one end of the tubing was inserted into the imaging medium, while the other needle was inserted into a syringe. One of the syringes/tubing combinations was used to draw buffer from the sample chamber and the other was used to infuse new buffer/imager strands to the system. Then, the blinking dataset of the second channel was recorded at 20 frames/s for 20,000-30,000 frames.

Single-molecule localization analysis
Here, the axial plane PSF was fitted with a 2D elliptical Gaussian function, where u k (x, z) is the expected photon count for a given pixel k, I 0 is the peak intensity, (x 0 , z 0 ) is the location of the fluorophore in the axial plane, w x and w z are the widths of the PSF in the x and z directions respectively, and b is the fluorescence background.
Maximum likelihood estimator (MLE) was used to fit data with Poisson statistics through maximizing the likelihood function, where the likelihood of the parameter θ given the data D is modeled as Poisson distribution given by the model u k and the observed counts d k at pixel k. Drift correction was performed using a 2D redundant cross-correlation method [37,38]. The entire localization dataset was split into n segments (each segment contains 2,000 frames). For each segment, a 2D histogram was generated where the value in each pixel represents the number of localizations falling within that pixel. Cross-correlation was calculated to find the shift distance between each pair of segments, where the correlation peak with sub-pixel resolution was identified after Fourier interpolation [37]. In this process, for each pair of segments (i and j), the shift distance can be obtained as (x shift )}, forming an overdetermined system for (n − 1) independent shifts between adjacent segments (r 1 , r 2 , · · · , r n−1 ). A temporary solution was found to minimize the overall error of the system using least square. This temporary solution was substituted back into all equations and residual errors were calculated. An equation from the system was removed if its error was larger than an empirical threshold (60 nm). This process was repeated until no equation gave an error larger than the threshold or the system of linear equations was no longer at full rank. In the case of two-channel imaging using exchange-PAINT, single molecule data were acquired first for Tom20 and then for α-tubulin. For the first channel, we corrected the sample drift by 2D redundant cross-correlation method with the reference being the last data segment. For the second channel, the first segment was used as the reference to correct the drift. The sample was tightly fixed on the custom-made sample stage and we assumed that the stage shift during the brief buffer changing period was negligible. Future applications requiring long term stability of the system will benefit from active stage stabilization techniques, for example using transmission images [39] or large polystyrene beads for high precision feedback [40].  between (x , y , z ) and (x, y, z) can be described as, In the axial plane SMLM system, α = 90 • , where (x , y , z ) = (x, z, y). Therefore, the x-z plane of the rotated intermediate image will be detected on the final image plane.
The detected image from a point source of unpolarized light, i.e. the 3D intensity point spread function I(x, y, z), can be derived as: where ì I A and ì I B are two vectors initially defined in Kim et al. [41,Eq. (9)] with the spherical coordinate system. Here, we modified the 2D intensity PSF into 3D together with the Spherical-to-Cartesian Coordinate transformation (k x , k y , k z ) ≡ n λ sinθcosφ, n λ sinθsinφ, − n λ cosθ , then Eq. (4) can be derived as: where j = A or B, k x k y k z -axes are parallel with x y z -axes respectively, P(k x , k y ) is the pupil function in k space with a length ρ = n/λ, n is the refractive index of the immersion medium, and λ is the emission wavelength in vacuum [42,43]. The detailed description of the pupil function P(k x , k y ) and electric field ì E j , j = A or B are defined in Kim et al. [41,].

Working principle and point spread function of the axial plane SMLM
Our design began with the axial plane optical microscopy (APOM) method [28]. A vertical slice of fluorescently stained biological specimen was excited by a light sheet generated using a high-NA objective lens (NA = 1.4). The emitted fluorescence signal was collected by the same objective lens Obj1, making the system compatible with regular coverslip-mounted samples ( Fig. 1(B)). A remote high-NA objective lens (NA = 1.4) was used for both forming an intermediate image with its focal point coinciding on a mirror surface oriented 45°with respect to its optical axis, and collecting the backward reflected fluorescence signal from the 45°tilted mirror. This design effectively rotated the axial plane image to the lateral focal plane of the remote objective ( Fig. 1(C) and (D)). To double the efficiency of the photon collection from single-molecule emissions, compared with previous APOM design [28], we combined a polarizing beam splitter with a QWP such that the p-polarized single-molecule fluorescence signal was transmitted, double passed through the QWP before being detected by an sCMOS camera. This design also allowed imaging both lateral and axial planes simultaneously ( Fig. 1(A)). Comparing to conventional lateral plane SMLM, the presented axial plane SMLM decreased the light efficiency, ∼7% fluorescence emission was detected on account of 50% loss at the polarizing beam splitter, 36% at the remote objective and 78.36% because of the asymmetric pupil ( Fig. 1(E)). Despite the significant loss of emitted photons, a reasonable signal-to-noise ratio (SNR) of single molecule emissions was achieved because of the compatibility of light-sheet illumination. This allowed us to reconstruct an axial plane super-resolution image without scanning through localizing the detected stochastically switching but isolated single-molecule emission patterns in the x-z plane.
To characterize the PSF in the axial plane SMLM system, we extended a previous 2D analysis of APOM response [41] to 3D and compared the theoretical response function with the experimentally obtained PSF using fluorescence beads by scanning the sample stage in the y direction, perpendicular to the imaged x-z plane. In comparison with an ideal PSF, we observed that our experimental PSF showed high similarity with the theoretical one, however, with a small amount of spherical aberration (127 mλ, Fig. 1(F), and Appendix, Fig. 5(A)). Due to the asymmetric loss of fluorescence within the collection cone of the remote objective, induced by the 45°angle of the tilted mirror, the detected axial PSF evolved asymmetrically along the y-axis (Appendix, Fig. 5(A)). When the designed focus of the remote objective coincided  on the mirror surface, we observed the elliptical shaped emission pattern with its long axis aligned with the z (axial) direction and the short axis aligned with the x direction (directions referencing to the sample plane). This elliptical PSF shape maintained within a y-range of 800 nm around the focal plane. The emission pattern started to differ significantly from the elliptical shape at ± 500 nm away from the focus. Considering that our axial plane SMLM system utilized light-sheet illumination with full width half maximum (FWHM) of 1.19 ± 0.04 µm (Appendix, Fig. 5(E)-(H)), we reasoned that the detected PSF from this narrowly illuminated sheet could be modeled as elliptical Gaussian with its long and short axes along the z and x directions respectively (Appendix, Fig. 5(I)). In addition, our axial plane SMLM system used two oil-immersion objectives with large NAs (1.4) for imaging specimens in aqueous environment which would result in the refractive index mismatch aberrations [44]. Under our present imaging conditions, we observed the elliptical PSF shape could be maintained in z direction over 20 µm which was currently sufficient for single-molecule super-resolution experiments demonstrated here and for future live cell imaging tasks.

Super-resolution imaging of whole cells with the axial plane SMLM
To demonstrate the super-resolution capability of the light-sheet illuminated axial plane SMLM system, we imaged immuno-labeled α-tubulin with Alexa Fluor 647 in fixed COS-7 cells. In comparison to the diffraction-limited axial plane image, the super-resolution image displayed distinct and isolated individual tubular structures and the overall interlaced microtubule network surrounding the nucleus -spherical shaped void showing little fluorescent staining, in the x-z plane (Fig. 2(A) and (B)). Given the outer diameter of a microtubule structure as ∼25 nm, and the size of primary and secondary labeling as ∼17.5 nm, the expected width (FWHM) of a cross-section profile of a microtubule can range from 70 to 92 nm depending on the achieved localization precision [45]. We quantified these resolved individual tubular structures by their intensity profiles in both lateral and axial directions (Fig. 2(D), and Appendix, Fig. 6(A) and (B)) and the average FWHMs in the lateral and axial directions are 88 ± 13 nm and 103 ± 12 nm (mean ± std., n = 30) respectively.
We further imaged the nucleoporin protein Nup98, part of the nuclear pore complex (NPC), localized near the center channel of the NPC scaffold, immuno-labeled with Alexa Fluor 647 in fixed COS-7 cells. We observed that the axial plane SMLM system was capable of resolving individual labeled NPCs as isolated clusters in both lateral and axial directions whereas the conventional axial plane image resulted in diffraction-limited resolution where only the NPC-decorated membrane contour was discernible (Fig. 2(E) and (F)). From the quantitative measurements in both lateral and axial directions of thirty distinctly resolved NPCs (Appendix, Fig. 6(C) and (D)), we found the FWHMs in the lateral and axial directions as 79 ± 6 nm and 103 ± 16 nm (mean ± std., n = 30) respectively. There was no significant difference in these measurements with respect to the imaging depth (∼8 µm) (additional super-resolved images are shown in Appendix, Fig. 7).

Multiplexed super-resolution imaging with the axial plane SMLM and exchange-PAINT
DNA-PAINT and exchange-PAINT offer a robust and versatile labeling approach for SMLM, owning to the transient binding of "imager" strands to their complementary "docking" strands [30,31]. Comparing to conventional switchable SMLM probes, DNA-PAINT exhibits effectively nonexistent photobleaching with the near infinite supply of fresh "imager" strands from the buffer, and allows straightforward sequential imaging of multiple targets, namely exchange-PAINT [31]. However, these free-diffusing "imager" strands significantly contribute to fluorescence background during single-molecule imaging, limiting the majority of its applications to biological structures near the coverslip surface. The axial plane SMLM system allows super-resolution imaging of the axial plane throughout a whole cell and is compatible with DNA-PAINT and exchange-PAINT by virtue of its light-sheet illumination providing a highly confined excitation (light sheet FWHM 1.19 ± 0.04 µm, Appendix, Fig. 5(E)-(I)). The axial plane SMLM excites single molecules only at the axial optical section that is in focus on the detector but not the freely diffusing probes nor the out-of-focus docked probes and therefore minimizes the background.
To demonstrate DNA-PAINT imaging using the axial plane SMLM, we imaged immuno-labeled α-tubulin with DNA-PAINT "docking" strands in fixed COS-7 cells, and with corresponding "imager" conjugated with Atto 655. From the acquired single-molecule blinking dataset, we observed a nearly 11-fold reduction of fluorescence background when comparing between lightsheet and epi illumination ( Fig. 1(G) and (H)). While the significant fluorescence background, resulting from an ocean of free diffusing probes and out-of-focus docked probes, made singlemolecule emissions difficult to observe by using epi illumination, it could be clearly visualized under light-sheet illumination in this system. In the reconstructed super-resolution image of microtubules network, we quantified the resolved individual tubulin profiles in both lateral and axial directions as 75 nm and 106 nm respectively (Appendix, Fig. 8).
The demonstrated compatibility of the axial plane SMLM with oligo-based probes allowed us to perform multiplexed super-resolution imaging using exchange-PAINT. Here we demonstrated multiplexed super-resolution imaging in the axial plane resolving both microtubules and mitochondria labeled with orthogonal DNA "docking" strands P1 and P4 in fixed COS-7 cells, respectively. We obtained multiplexed axial plane super-resolution images through whole cells in the axial cross-section, displaying super-resolved α-tubulin and Tom20 through the same axial plane ( Fig. 3 and Appendix, Fig. 9). We found that the axial distribution of these two cellular targets was organized spatially exclusively with individual microtubules extending alongside the Tom20-decorated mitochondria surface. In addition, by resolving multiple protein species with the same fluorescent dye, exchange-PAINT allowed multiplexed super-resolution imaging with a single laser source in excitation and without chromatic aberrations in emission.

Resolving the nucleus of glioblastoma multiforme cells with the axial plane SMLM
Our presented SMLM system allows efficient super-resolution profiling of thick specimens in the axial plane. We demonstrated this unique compatibility of the axial plane SMLM with thick specimens by imaging immuno-labeled Nup98 in glioblastoma multiforme cells (GBM10 cells), whose cell height ranges from 10 to 30 µm usually consisting a large nucleus (>10 µm). Conventional imaging of a GBM10 cell throughout its extended height requires scanning in the axial direction with a step size usually 400 nm and subsequently ∼50 optical sections are required to reveal its axial profile. Using the axial plane SMLM, we were able to resolve Nup98-labeled nuclear envelope throughout the entire depth of the tumor cell without axial scanning. At various imaging depths, we found individual pores were distinctly resolved whereas their diffraction-limited counterparts blurred the individual NPCs into a continuous contour of the nucleus (Fig. 4). We observed consistent achieved resolution with different imaging depths as shown in Nup98 axial plane SMLM images with total depths of 14 µm and 17 µm (Fig. 4(A) and (D)). We measured localized Nup98 cluster profiles with a single reconstruction in the lateral direction at different imaging depths and found FWHMs of 66 nm and 69 nm at 1.8 µm and 16.0 µm respectively. The corresponding axial profiles displayed FWHMs of 93 nm and 95 nm respectively (Fig. 4(F)). Further studies of GBM10 immune response using multiplexed axial plane SMLM could allow axial profiling of multiple proteins (e.g. antigens) on these treatment-refractive cells to infer mechanisms driving immune response and to guide the development of more effective and safer immunotherapies.

Conclusion
In this work, we developed an axial plane super-resolution system by integrating single-objective light-sheet illumination and axial plane optical imaging with single-molecule detection and localization. We demonstrated the proposed system by resolving the axial distribution of different subcellular targets such as microtubules, mitochondria and NPCs over the whole-cell depth. The presented axial plane SMLM is compatible with DNA-PAINT and exchange-PAINT due to its native light-sheet illumination and axial detection allowing multiplexed axial plane super-resolution imaging resolving multiple subcellular targets in a sequential manner. Due to the ability to achieve consistent resolution through an imaging depth, the system allows us to resolve individual nuclear pore complexes decorating the entire nuclear envelope in large glioblastoma patient-derived tumor cells. As a unique technique for efficient super-resolution profiling in the depth dimension, we hope the axial plane SMLM will be a useful tool to study intra-and extra-cellular structures and processes along the axial direction in both fixed and living specimens.