• Photonics Research
  • Vol. 12, Issue 4, 821 (2024)
Wei Shi, Yingchuan He, Jianlin Wang, Lulu Zhou, Jianwei Chen, Liwei Zhou, Zeyu Xi, Zhen Wang, Ke Fang, and Yiming Li*
Author Affiliations
  • Department of Biomedical Engineering, Southern University of Science and Technology, Shenzhen 518055, China
  • show less
    DOI: 10.1364/PRJ.514414 Cite this Article Set citation alerts
    Wei Shi, Yingchuan He, Jianlin Wang, Lulu Zhou, Jianwei Chen, Liwei Zhou, Zeyu Xi, Zhen Wang, Ke Fang, Yiming Li. Aberration correction for deformable-mirror-based remote focusing enables high-accuracy whole-cell super-resolution imaging[J]. Photonics Research, 2024, 12(4): 821 Copy Citation Text show less

    Abstract

    Single-molecule localization microscopy (SMLM) enables three-dimensional (3D) investigation of nanoscale structures in biological samples, offering unique insights into their organization. However, traditional 3D super-resolution microscopy using high numerical aperture (NA) objectives is limited by imaging depth of field (DOF), restricting their practical application to relatively thin biological samples. Here, we developed a unified solution for thick sample super-resolution imaging using a deformable mirror (DM) which served for fast remote focusing, optimized point spread function (PSF) engineering, and accurate aberration correction. By effectively correcting the system aberrations introduced during remote focusing and sample aberrations at different imaging depths, we achieved high-accuracy, large DOF imaging (8 μm) of the whole-cell organelles [i.e., nuclear pore complex (NPC), microtubules, and mitochondria] with a nearly uniform resolution of approximately 35 nm across the entire cellular volume.

    1. INTRODUCTION

    As a powerful tool for imaging biological samples, single-molecule localization-based super-resolution microscopy allows for optical observation of structural details with unprecedented resolution while providing in situ multi-targets imaging capability by multi-color labels [1,2]. By engineering the point spread function (PSF) with different shapes along the axial direction, single-molecule localization microscopy (SMLM) can easily achieve three-dimensional (3D) super-resolution imaging capability. The most widely used 3D super-resolution imaging method for determining the axial position of the single-molecule involves introducing astigmatism into the detection path using a cylindrical lens [3]. However, this technique has limited imaging depth of field (DOF) when a high numerical aperture (NA) objective is used (1  μm). Considering that the typical thickness of a mammalian cell is 5–10 μm [4], this limitation prevents imaging of the entire cellular organelles that can extend over several micrometers throughout the whole cell [5].

    To image samples thicker than the DOF of high-NA objective using SMLM [68], different implementation approaches have been developed, such as z-scan with thin sections imaging [4,911], multi-plane imaging [12,13], and engineered PSF with large DOF imaging [1417]. Z-scan imaging normally requires physically moving the sample or objective to scan at the axial positions of cells and then combining optical sections to obtain the complete cellular structure. It is therefore prone to stitching artefacts due to vibrations caused by stage movement during repeated scans of different sections. Multi-plane detection allows for simultaneous imaging of different focal planes along the axial direction to achieve volumetric imaging. However, compared to single-plane imaging, multi-plane detection often suffers from lower signal-to-noise ratio (SNR) since the total number of photons is distributed across different imaging planes [13]. Thick sample imaging can also be achieved by engineered PSF with large DOF without the need to scan the sample. There are two main ways to generate these engineered PSFs. One is to fabricate a transmission phase mask with a specialized pattern [15], placed in the Fourier plane of the microscope. Another way is to employ a programmable phase modulator, such as a spatial light modulator and deformable mirror (DM). Although both ways can be used to engineer PSFs with desired DOF, programmable phase modulators are often preferred due to their flexibility and ability to correct for both system and sample-induced aberrations. Especially with the development of DM-based optimal PSF engineering, we were able to design an optimal PSF with predefined DOF without compromising the loss of photons [17]. However, there is normally a trade-off between the localization precision and DOF of the PSF optimized. PSFs with longer DOF normally spread in a larger area which leads to a reduction in SNR and loss of localization precision.

    In recent years, adaptive optics (AO) techniques have been widely applied in optical microscopy [18]. In an approach known as remote focusing, fast 3D imaging of thick samples can be achieved by dynamically adjusting the focal plane without physically moving the objective lens or sample [1923]. This process can be executed at rates of a few kilohertz using a DM or other variable optical elements [19,24,25], avoiding the need for physical movement of the sample or objective. In SMLM, this method has demonstrated the capacity to record volumetric data of whole cells up to 10 μm while maintaining axial sample stabilization using a standard widefield microscope equipped with a DM [20]. However, system aberrations often arise when applying the defocus term to the vibrational device [22,26]. Additionally, the inhomogeneous refractive indices of biological tissues can lead to the blurring and distortion of single-molecule emission patterns [11,27]. All these issues can give rise to imaging artifacts and compromise the achievable resolution of SMLM. Hence, accurate correction and modeling of these aberrations are key to acquire high-accuracy super-resolution images across the entire volume.

    Here, we employed a DM to engineer PSFs with optimal 3D localization precision with different DOFs (DMO PSF). Furthermore, the DM was used to fast remote focus to longitudinally record axial stacks of whole cells across an extended range with high accuracy, without the need to physically move the sample. Crucially, we accurately calibrated and corrected system aberrations introduced during refocusing at different imaging planes. Finally, we used in situ single-molecules to estimate the aberrations at each focal plane through global fitting of the pupil function. The globally fitted aberrations were then fed to the DM for sample-induced aberrations correction. In both silicone oil and oil objectives, we achieved consistent high-resolution 3D reconstructions at different depths and obtained a nearly uniform lateral Fourier ring correlation (FRC) [28] resolution (38  nm for silicone oil objective and 35  nm for oil objective).

    2. RESULTS

    A. Remote Focusing Principle and Optical Setup

    As shown in Fig. 1(a), we modified a standard 3D SMLM system by incorporating AO components (for details see Fig. 8 in Appendix B). Remote focusing based on AO is enabled by a DM, which is positioned in conjugate with the objective pupil plane. The DM has emerged as the most prevalent phase modulator for fluorescence detection due to its minimal photon loss and rapid wavefront control capabilities. In this design, the DM serves three purposes: PSF engineering, aberration correction, and focal position shifting for different imaging depths. For precise wavefront control, we meticulously calibrated the experimental DM influence function of each actuator using a Twyman–Green interferometer, which measures the surface deformation of the DM [17,29].

    Schematic diagram illustrating the remote focusing principle. (a) Simplified schematic of a microscope configured for remote focusing using a DM in the imaging path. (b) The DM modulates the wavefront by controlling its defocus phase to shift different refocused focal planes around the nominal focal plane along the optical z-axis. (c) The defocus phase root mean square error (RMSE) for the optical wavefront is indicated when remotely focusing to different axial positions within the imaging optical system.

    Figure 1.Schematic diagram illustrating the remote focusing principle. (a) Simplified schematic of a microscope configured for remote focusing using a DM in the imaging path. (b) The DM modulates the wavefront by controlling its defocus phase to shift different refocused focal planes around the nominal focal plane along the optical z-axis. (c) The defocus phase root mean square error (RMSE) for the optical wavefront is indicated when remotely focusing to different axial positions within the imaging optical system.

    Remote focusing is a high-speed volumetric scanning technique that axially shifts the imaging plane by applying a defocus phase at the pupil of the optical system. This technique avoids the need to move the specimen or objective, thereby enabling accurate volumetric super-resolution imaging. We employed the defocus mode of the DM to refocus the focal plane of the imaging optical system. First, we defocused the beads on the coverslip using the objective’s z-piezo stage. Then, we compensated for this defocus by adjusting the DM’s defocus mode to refocus the beads. We calibrated the relationship between the amplitude of the applied defocus wavefront in the DM and real axial shift of the imaging plane [Fig. 1(b)]. As shown in Fig. 1(c), the amplitude of the applied defocus wavefront does not correlate linearly with the changes in the imaging plane. Especially for refocusing planes far away from the zero axial shift focal plane, the deviation is more pronounced. This is probably due to the limited actuators in the DM which cannot accurately modulate the ideal defocus wavefront. Therefore, it is very important to correct the residual aberrations at each focal plane after applying the defocus wavefront to the DM.

    B. System Aberration Correction of Remote Focus Using Fluorescent Beads

    Astigmatism-based 3D SMLM imaging is the most widely used 3D super-resolution imaging method as it can be achieved by simply introducing a cylindrical lens into the imaging path. The DOF of the astigmatism PSF is approximately at the range of 1.2 μm around the focus [4]. Beyond this depth range, axial sample scanning [9] or optimized extended DOF PSF is necessary [30], as molecules outside the focus cannot be efficiently identified and localized. Furthermore, the 3D imaging capability of astigmatism is suboptimal as the localization precision decreases rapidly at the defocusing positions. In this study, we integrated our previously published DM-optimized PSF (DMO PSF) framework [17], designed for varying imaging depths, with remote focusing to achieve high-quality 3D super-resolution across the entire cell. To obtain the DMO PSF, we employed the influence function of the DM actuators as the solution space to optimize the pupil function of the engineered PSF by minimizing its 3D Cramér–Rao lower bound (CRLB). Subsequently, the DM was used to modulate both the DMO PSF and astigmatism PSF [Figs. 2(a)–2(d)]. Both experimental PSFs were built by axially scanning the beads on the coverslip and interpolating the beads z stack with cubic spline [31]. The localization precision of both experimental PSFs was then compared. As shown in Figs. 2(e)–2(h), the 2 μm DMO PSF exhibited a similar localization precision near the focus compared to the astigmatism PSF. However, the performance of astigmatism deteriorates significantly when the imaging plane is slightly away from the focus. In contrast, the DMO PSF maintains a relatively uniform resolution throughout the optimized axial range. The performance of experimental PSFs in Figs. 2(e)–2(h) agrees well with the theoretical calculations [17].

    Comparison of the localization precision between experimental astigmatism PSF and DMO PSF. (a) Averaged astigmatism PSF of fluorescent beads on the coverslip, with PSFs at axial positions from −800 to 800 nm. (b) The pupil function corresponding to the astigmatism PSF of beads. (c) Averaged DMO PSF of fluorescent beads, identical to the range described in (a). (d) Same as (b) for DMO PSF. (e) x, y, and z localization precisions at different axial positions for astigmatism PSF. For the CRLB calculation, we used 3000 photons and 30 background photons to simulate the typical photon flux of fluorescent dyes. (f) 3D CRLB of astigmatism PSF. (g) x, y, and z localization precisions for DMO PSF. (h) 3D CRLB of DMO PSF.

    Figure 2.Comparison of the localization precision between experimental astigmatism PSF and DMO PSF. (a) Averaged astigmatism PSF of fluorescent beads on the coverslip, with PSFs at axial positions from 800 to 800 nm. (b) The pupil function corresponding to the astigmatism PSF of beads. (c) Averaged DMO PSF of fluorescent beads, identical to the range described in (a). (d) Same as (b) for DMO PSF. (e) x, y, and z localization precisions at different axial positions for astigmatism PSF. For the CRLB calculation, we used 3000 photons and 30 background photons to simulate the typical photon flux of fluorescent dyes. (f) 3D CRLB of astigmatism PSF. (g) x, y, and z localization precisions for DMO PSF. (h) 3D CRLB of DMO PSF.

    We then applied the DM with the defocus wavefront in addition to the optimized pupil function to remotely scan the samples for volumetric super-resolution imaging. To investigate the PSF after adding the defocus wavefront, we imaged the beads on the coverslip with the DM-applied defocus wavefront and refocused the beads using the objective z-stage. As shown in Fig. 3(a), the beads were modulated to the DMO PSF and defocus wavefront of 500  nm (4.62 rad, wavelength 680 nm) amplitude. The beads were then refocused by moving the objective 2 μm above the nominal focal plane. However, the refocused PSF showed significant difference compared to the one at zero-focus position [Figs. 2(c) and 3(a)]. We infer that this is probably due to the mismatch between the theoretical and experimental defocus mode, which leads to residual aberrations. To investigate the residual aberrations introduced besides defocus, we then acquired bead stacks in a range of ±1  μm relative to the focal plane for calibration. We utilized a GPU-based vectorial PSF fitter to fit each bead stack by maximum likelihood estimation (MLE) to retrieve the Zernike-based aberration coefficients [16]. Furthermore, we employed globLoc [32] for the global fitting of images from the bead stacks, enabling high-accuracy estimation of aberrations.

    Correcting aberrations during refocusing with DMO PSF using a DM at various axial positions of the beads on the coverslip. (a) The beads were refocused using a DM to a position 2 μm above the nominal focal plane. (b) The pupil function for the refocused beads. (c) Fitted 21 Zernike coefficients for the beads stack. (d) CRLB for the experimental DMO PSF, calculated using 3000 photons and 30 background photons to simulate the typical photon flux of fluorescent dyes. A comparison of CRLB calculations based on the typical photon flux of fluorescent proteins is shown in Fig. 9 of Appendix C. (e)–(h) The PSF shape, pupil function, Zernike coefficients, and localization precision after aberration correction, respectively. The CRLB exhibits significant improvements, particularly in the z localization precision.

    Figure 3.Correcting aberrations during refocusing with DMO PSF using a DM at various axial positions of the beads on the coverslip. (a) The beads were refocused using a DM to a position 2 μm above the nominal focal plane. (b) The pupil function for the refocused beads. (c) Fitted 21 Zernike coefficients for the beads stack. (d) CRLB for the experimental DMO PSF, calculated using 3000 photons and 30 background photons to simulate the typical photon flux of fluorescent dyes. A comparison of CRLB calculations based on the typical photon flux of fluorescent proteins is shown in Fig. 9 of Appendix C. (e)–(h) The PSF shape, pupil function, Zernike coefficients, and localization precision after aberration correction, respectively. The CRLB exhibits significant improvements, particularly in the z localization precision.

    In this study, we retrieved the coefficients of all 21 tertiary Zernike polynomials from the average experimental bead stack [Fig. 3(c)]. Due to the extra aberrations introduced besides defocus, the pupil function of the refocused DMO PSF [Fig. 3(b)] is quite different from that at the zero focal plane [Fig. 2(d)]. The residual aberrations also led to the degradation of localization precision, particularly in the z-axis direction [Fig. 3(d)]. Therefore, we compensated for the 21 Zernike aberrations using the DM to counteract the additional aberrations introduced after refocusing, making the aberrations at the refocusing position closer to those of the beads on the coverslip [Figs. 3(e) and 3(f)]. As shown in Fig. 3(h), the resulting localization precisions are close to those of the PSF at zero focal plane [Fig. 2(g)].

    We then demonstrated the effectiveness of the residual aberration correction after applying the defocus wavefront using the nucleoporin Nup96 in U2OS cells, a structure often used as a quantitative reference [33]. The top nuclear envelopes of the Nup96-SNAP-AF647 labeled nuclear pore complexes (NPCs) were imaged. The imaging depth is about 4 μm above the coverslip. As shown in Fig. 10(a) of Appendix D, the single-molecule images are very aberrated and the SNR is considerably low due to the additional aberrations introduced when applying the defocus wavefront to the DM. We then applied the calibrated residual aberration correction using the beads on the coverslip as described before. After aberration correction, the single-molecules are clearly visible, and the SNR of the image is improved significantly [Fig. 10(b) in Appendix D]. The number of localized single-molecule events increased sixfold, and the signal background ratio (SBR) of single-molecule events improved by nearly 2.5 times [Fig. 10(c) in Appendix D], with the localization and quantification analysis of single-molecule images implemented in SMAP [34]. The experimental astigmatic PSFs with and without residual aberration correction after applying the defocus wavefront are shown in Fig. 11 of Appendix E. We then applied these two PSFs to image biological samples. As shown in Fig. 12 of Appendix F, the double ring structure of NPCs can be clearly resolved after residual aberration correction [Fig. 12(f) in Appendix F] while it was hardly resolved when only the defocus wavefront was applied to the DM [Fig. 12(c) in Appendix F].

    C. Whole-Cell 3D Imaging with Remote Focusing Using Beads PSF Model

    We first demonstrated the whole-cell imaging capability of our approach using a silicone oil objective where we can match the refractive index of the sample medium with silicone oil [35,36]. The whole-cell nucleus was imaged with 2 μm DMO PSF. Five imaging planes spaced by 1 μm were recorded by remote focusing between different imaging planes using the DM (Appendix A). The five optical section data were then stitched together by redundant cross-correlation to reconstruct the whole nucleus (Appendix A). As shown in Fig. 4, we were able to clearly reconstruct the ring structure in both top and bottom nuclear envelope with FRC resolution of 38 nm [Fig. 4(b)]. In the side-view of the nucleus, we were also able to resolve the double ring structure in the entire nucleus over a depth of 5  μm. Moreover, both the reconstructed diameter and the spacing between the upper and lower rings of the Nup96 structure agree with the reference structure [Figs. 4(h), 4(i) and 4(j), 4(k)]. We then imaged the nucleus with 6 μm DMO PSF (Fig. 13 in Appendix G) for comparison and observed a reduced FRC resolution (50  nm, Fig. 14 in Appendix H). It is also difficult to observe the double ring structure of Nup96 from the side-view images.

    Whole-cell 3D super-resolution imaging of NPC. (a) Overview of the panoramic whole-cell 3D imaging of NPC using DMO PSF, merging five optical sections. (b) FRC analysis of the region boxed in (a). (c) Side-view cross section of the region denoted by the dashed line in (a). (d), (e) Magnified view of the area denoted by the box in (c). (f) View of the bottom surface of the boxed area in (a). (g) View of the top surface of the boxed area in (a). (h), (i) Intensity profile along the white dashed lines in (f) and (g). (j), (k) Intensity profile along the white dashed lines in (d) and (e). The data were acquired from 8000 frames per cycle over 20 cycles from five optical sections, with 120 mW laser power and 20 ms exposure time.

    Figure 4.Whole-cell 3D super-resolution imaging of NPC. (a) Overview of the panoramic whole-cell 3D imaging of NPC using DMO PSF, merging five optical sections. (b) FRC analysis of the region boxed in (a). (c) Side-view cross section of the region denoted by the dashed line in (a). (d), (e) Magnified view of the area denoted by the box in (c). (f) View of the bottom surface of the boxed area in (a). (g) View of the top surface of the boxed area in (a). (h), (i) Intensity profile along the white dashed lines in (f) and (g). (j), (k) Intensity profile along the white dashed lines in (d) and (e). The data were acquired from 8000 frames per cycle over 20 cycles from five optical sections, with 120 mW laser power and 20 ms exposure time.

    To further demonstrate that image quality is maintained throughout the entire cell, we imaged microtubules in COS-7 cells. The 3D reconstruction was assembled from seven optical sections with step sizes of 1 μm, resulting in a whole-cell microtubule super-resolution image with a DOF of 8 μm [Fig. 5(a)]. We were able to clearly observe the microtubule structures throughout the whole cell, with the ability to resolve finer structural details of the microtubules as shown in both the top and bottom layers of the cell [Figs. 5(b) and 5(c)]. Benefitting from the stitching method by redundant cross-correlation, no significant stitching artifacts were observed in the axial profile of the whole-cell microtubule volume, which is aligned of seven optical sections [Fig. 5(d)].

    Whole-cell 3D super-resolution imaging of microtubules. (a) Overview of the panoramic whole-cell 3D imaging of microtubule using DMO PSF, merging seven optical sections. (b) Zoomed bottom surface view of the boxed area denoted in (a). (c) Zoomed top surface view of the boxed area denoted in (a). (d) Side-view cross section of the region denoted by the dashed line in (a). The data were acquired from 7000 frames per cycle over 40 cycles from seven optical sections, with 200 mW laser power and 15 ms exposure time.

    Figure 5.Whole-cell 3D super-resolution imaging of microtubules. (a) Overview of the panoramic whole-cell 3D imaging of microtubule using DMO PSF, merging seven optical sections. (b) Zoomed bottom surface view of the boxed area denoted in (a). (c) Zoomed top surface view of the boxed area denoted in (a). (d) Side-view cross section of the region denoted by the dashed line in (a). The data were acquired from 7000 frames per cycle over 40 cycles from seven optical sections, with 200 mW laser power and 15 ms exposure time.

    D. In situ Aberration Correction of Remote Focus Using Blinking Single-Molecules

    The bead-based PSF model normally works well only when the in situ PSF model within the samples is close to the PSF model on the coverslip. This is often valid when the refractive index of the sample medium is matched to the objective immersion medium. However, the refractive index of the sample medium is often not matched to the objective immersion medium. Especially for a high-NA objective where oil is used as the immersion medium, strong spherical aberrations will be introduced when the focal plane is away from the coverslip. These additional aberrations often result in reduced resolution if they are not corrected. To address this, we employed our recently developed in situ PSF modeling method, uiPSF, to model the PSF directly from the single-molecule blinking events [37]. Subsequently, we carried out in situ aberration correction for different imaging planes to further enhance the imaging resolution across the entire cell.

    Here, we employed a 1.5 NA oil-immersion objective to image the whole cells. Similar to the silicone oil objective imaging, we first calibrated the residual aberrations at different focal planes after applying the defocus wavefront to the DM. As shown in Figs. 6(a) and 6(b), an ideal DMO PSF could be obtained after applying the residual aberration correction. We then applied this PSF to image COS-7 cell at a depth of about 2.5 μm [Fig. 6(c)]. The single-molecule blinking data were used to estimate the in situ PSF model. The pupil function was estimated by globally fitting the single-molecule images via inverse modeling by uiPSF. As shown in Fig. 6(g), the fitted pupil function showed strong spherical and coma aberrations compared to the beads-based PSF on the coverslip. This is probably due to the refractive index mismatch between the sample medium and oil immersion medium.

    In situ aberration correction using blinking single-molecules under remote focusing. (a) The beads on the coverslip were refocused using the DM to 1 μm above the nominal focal plane, with a pre-correction for aberrations applied to the imaging system. (b) Pupil function derived from the PSF of refocused beads. (c) The refocused PSF was obtained from single-molecule blinking data of immunofluorescence-labeled TOM20, recorded at a depth of 2.5 μm in COS-7 cells. (d) Pupil function calculated from the in situ PSF. (e), (f) The PSF shape and its corresponding pupil function after in situ aberration correction. (g) The fitted 21 Zernike coefficients were obtained from the PSF of refocused beads after system aberration correction, the in situ blinking single-molecule, and the in situ PSF after aberration correction using blinking single-molecule.

    Figure 6.In situ aberration correction using blinking single-molecules under remote focusing. (a) The beads on the coverslip were refocused using the DM to 1 μm above the nominal focal plane, with a pre-correction for aberrations applied to the imaging system. (b) Pupil function derived from the PSF of refocused beads. (c) The refocused PSF was obtained from single-molecule blinking data of immunofluorescence-labeled TOM20, recorded at a depth of 2.5 μm in COS-7 cells. (d) Pupil function calculated from the in situ PSF. (e), (f) The PSF shape and its corresponding pupil function after in situ aberration correction. (g) The fitted 21 Zernike coefficients were obtained from the PSF of refocused beads after system aberration correction, the in situ blinking single-molecule, and the in situ PSF after aberration correction using blinking single-molecule.

    Subsequently, we compensated for the aberrations by controlling the DM and then re-collected single-molecule data from the same location for in situ PSF estimation. Through iterative corrections, various aberrations were effectively minimized [Fig. 6(g)]. The resulting single-molecule data displayed a 3D PSF shape [Fig. 6(e)] and pupil function [Fig. 6(f)] similar to those of the beads on the coverslip.

    E. Whole-Cell 3D Imaging with Remote Focusing Using the In situ PSF Model

    We then combined remote focusing with uiPSF-based aberration correction using the DM and an oil objective to image the whole-cell mitochondria. Here, we employed 2 μm DMO PSF and imaged the whole-cell mitochondria with three remote focusing imaging planes spaced by 1 μm. For each refocused position, in situ PSF aberrations were corrected at depths of 0.5, 1.5, and 2.5 μm. Data from these three optical sections were acquired using rapid remote focusing cycles, controlled by the DM. The reconstructed 3D image was assembled from three optical sections, which were then aligned with adjacent optical sections, resulting in a whole-cell mitochondria super-resolution image with a DOF of 3 μm [Fig. 7(a)]. With our approach, the interconnected mitochondrial network was clearly resolved, and the xz cross sections revealed the membrane contour of mitochondria in the axial direction [Fig. 7(b)]. Thanks to the improvement in resolution through in situ aberration correction across the whole cell, we achieved consistently high-resolution 3D reconstructions of organelles at different depths [Figs. 7(c)–7(e)]. Resolution measurements were performed in regions proximate to both the lower [Fig. 7(f)] and upper [Fig. 7(g)] surfaces of the cell using the FRC method, and we obtained a nearly uniform lateral resolution of approximately 35 nm.

    Whole-cell 3D super-resolution imaging of mitochondria. (a) Overview of the panoramic whole-cell 3D imaging of mitochondria using DMO PSF by merging three optical sections. (b) Side-view cross section of the region indicated by the dashed line in (a). (c)–(e) Zoomed-in side-view cross sections from distinct depth areas, as indicated by the dashed lines in (a). (f), (g) FRC analysis of regions enclosed by the boxes in (a), with (f) near the cell’s bottom surface and (g) near the top. The data were acquired from 3000 frames per cycle over 70 cycles from three optical sections, with 200 mW laser power and 15 ms exposure time.

    Figure 7.Whole-cell 3D super-resolution imaging of mitochondria. (a) Overview of the panoramic whole-cell 3D imaging of mitochondria using DMO PSF by merging three optical sections. (b) Side-view cross section of the region indicated by the dashed line in (a). (c)–(e) Zoomed-in side-view cross sections from distinct depth areas, as indicated by the dashed lines in (a). (f), (g) FRC analysis of regions enclosed by the boxes in (a), with (f) near the cell’s bottom surface and (g) near the top. The data were acquired from 3000 frames per cycle over 70 cycles from three optical sections, with 200 mW laser power and 15 ms exposure time.

    3. DISCUSSION

    In this work, we employed a DM to perform super-resolution imaging of the whole cells by remote focusing. We systematically corrected the residual aberrations after applying the defocus wavefront and sample-induced aberrations with uiPSF. Furthermore, we utilized the DMO PSF, which could achieve optimal 3D resolution within a predefined DOF using the DM. These resulted in a significant SNR improvement in single-molecule data post aberration correction. We demonstrated the performance of our approach by conducting whole-cell imaging of nucleoporin Nup96, microtubules, and mitochondria. Compared to the PSF directly optimized by the DM with a large DOF of 6 μm, which exhibited a resolution of approximately 50 nm, our reconstructed super-resolution images exhibited a marked improvement in resolution (35  nm). Moreover, we employed uiPSF to precisely estimate the sample-induced aberrations using the in situ blinking single-molecules when a high-NA oil objective was used.

    In summary, the remote focusing enabled high-resolution whole-cell imaging, particularly in scenarios with high labeling density where large DOF PSF imaging could result in a strong background. With the precise in situ PSF estimation, it can also be further extended to tissue imaging [38,39]. Although we showed superior performance of remote focusing compared to the whole-cell imaging with large DOF PSF, it is worth noting that there is a balance between the number of focus planes imaged and the number of molecules acquired due to photobleaching. A combination of light sheet illumination with remote focusing could minimize the photobleaching while maintaining high-resolution imaging. With careful calibration of the defocus aberration at different wavelengths, multi-color imaging could also be integrated with our method. We believe that our study will promote the use of adaptive optics in SMLM techniques, especially in biological applications that require both a large DOF and high-resolution imaging.

    Acknowledgment

    Acknowledgment. The authors thank Shuang Fu and Mengfan Li for helpful suggestions on DM calibration and PSF optimization.

    APPENDIX A: METHODS AND EXPERIMENTAL DETAILS

    Optical Setups

    In this work, we performed SMLM imaging at room temperature (24°C) on a custom-built microscope (Fig. 8 in Appendix B). We coupled the excitation light (iBEAM-SMART-405-S, 150 mW, iBEAM-SMART-640-S-HP, 200 mW, TOPTICA Photonics) into a single-mode fiber (P3-405BPM-FC-2, Thorlabs) using a dichroic mirror (DMLP 425, Thorlabs) and a fiber coupler (PAF2-A4A, Thorlabs). We adjusted the translation stage to vary the fiber output for different illumination angles. The illumination beam was filtered by a laser clean-up filter (ZET405/488/561/640xv2, Chroma) to remove excess wavelengths and stray light. We used a pair of lenses L1 (75 mm) and L2 (400 mm) with a slit (SP60, Owis) at AP1 for beam collimation and reshaping. The beam was then reflected by the main dichroic mirror (ZT405/488/561/640rpcxt-UF2, Chroma) before entering the objective for sample illumination. The emitted fluorescence was collected by a high-NA objective (NA 1.35, UPLSAPO 100 XS or NA 1.5, UPLSAPO 100 XOHR, Olympus) and imaged by the tube lens TBL (TTL-180-A, Thorlabs) onto the AP2 confined by a slit (SP40, Owis). We utilized two bandpass filters (NF03-405/488/561/635E25 and FF01-676/37-25, Semrock) on the filter wheel to separate the emitted fluorescence from the excitation laser. We set up a 4f system (L3, 125 mm, L4, 75 mm) in the Fourier plane, with a deformable mirror (DM140A-35-P01, Boston Micromachines) for PSF engineering, aberration correction, and focal position shifting for different imaging depths. We could observe the objective’s back focal plane through lens L5 (40 mm) placed in front of the camera. Images were acquired using an sCMOS camera (Dhyana 400BSI V3, Tucsen) with a pixel size of 110 nm in the sample. Furthermore, we integrated a closed-loop focus lock control system into our optical setup. This system introduces the signal from a 785 nm laser (iBEAM-SMART-785-S, 125 mW, TOPTICA Photonics) reflected off the coverslip through a dichroic mirror (FF750-SDi02, Semrock) and detects it with a quadrant photodiode (SD197-23-21-041, Advanced Photonix Inc.), which then provides feedback control to a z-piezo stage (P-726.1CD, Physik Instrumente) to maintain focus stability.

    Detailed layout of the optical setup. M, mirror; DM, dichroic mirror; L, lens; TS, translation stage; FC, fiber coupler; Fiber, single-mode fiber; BFP, back focal plane; FW, filter-wheel; TBL, tube lens; AP, aperture; QPD, quadrant photodiode. The excitation lasers are first reflected by dichroic mirror DM1 and then coupled into a single-mode fiber through the fiber coupler FC. Before being reflected by the main dichroic mirror DM2 to enter the objective for sample illumination, a pair of lenses L1 and L2 along with a slit at AP1 are used for collimating and reshaping the beam. In the imaging path, the fluorescence collected by the objective is transmitted through the dichroic mirror DM2 and filtered by the filters on the filter wheel FW. It is then imaged onto the aperture AP2, which is confined by a slit, using a tube lens TBL. Subsequently, the fluorescence passes through a 4f system composed of lenses L3 and L4, with a deformable mirror installed in the Fourier plane, and ultimately, the fluorescence signal is detected by the camera. Additionally, a beam excited by a 785 nm laser, reflected off the coverslip, is detected by the quadrant photodiode QPD, which provides feedback control to the z-stage for focus locking.

    Figure 8.Detailed layout of the optical setup. M, mirror; DM, dichroic mirror; L, lens; TS, translation stage; FC, fiber coupler; Fiber, single-mode fiber; BFP, back focal plane; FW, filter-wheel; TBL, tube lens; AP, aperture; QPD, quadrant photodiode. The excitation lasers are first reflected by dichroic mirror DM1 and then coupled into a single-mode fiber through the fiber coupler FC. Before being reflected by the main dichroic mirror DM2 to enter the objective for sample illumination, a pair of lenses L1 and L2 along with a slit at AP1 are used for collimating and reshaping the beam. In the imaging path, the fluorescence collected by the objective is transmitted through the dichroic mirror DM2 and filtered by the filters on the filter wheel FW. It is then imaged onto the aperture AP2, which is confined by a slit, using a tube lens TBL. Subsequently, the fluorescence passes through a 4f system composed of lenses L3 and L4, with a deformable mirror installed in the Fourier plane, and ultimately, the fluorescence signal is detected by the camera. Additionally, a beam excited by a 785 nm laser, reflected off the coverslip, is detected by the quadrant photodiode QPD, which provides feedback control to the z-stage for focus locking.

    CRLB calculations before and after aberration correction of the beads on the coverslip during refocusing with DMO PSF using a DM. (a) The CRLB for the experimental DMO PSF was calculated using 800 photons and 10 background photons, simulating the typical photon flux of fluorescent proteins. (b) Localization precision after aberration correction. The CRLB showed significant improvements, particularly in the z localization precision.

    Figure 9.CRLB calculations before and after aberration correction of the beads on the coverslip during refocusing with DMO PSF using a DM. (a) The CRLB for the experimental DMO PSF was calculated using 800 photons and 10 background photons, simulating the typical photon flux of fluorescent proteins. (b) Localization precision after aberration correction. The CRLB showed significant improvements, particularly in the z localization precision.

    Comparison of raw single-molecule images before and after aberration correction, accompanied by the distribution of photons and background photons. (a) Snapshots of an unmodified 2D single-molecule image of Nup96-SNAP-AF647 labeled NPCs in U2OS cell samples refocused to 4 μm above the nominal focal plane. The inset shows a zoomed view of the fluorophore in the boxed area. (b) Snapshots of the 2D single-molecule image after aberration correction, with a significantly improved SNR in the image. The inset clearly shows the successful recovery of the distorted PSF. (c) The photon distribution analysis from 2000 frames of 2D single-molecule images before and after aberration correction reveals that 1644 single-molecule events with an average photon count of 335 were localized before correction, whereas 9680 single-molecule events with an average photon count of 711 were localized after correction. (d) The background photon distribution analysis from the same 2000 frames of 2D single-molecule images before and after aberration correction showed an average background photon count of 29 before correction and an average photon count of 25 after correction. The single-molecule raw images were acquired at 120 mW laser power and 20 ms exposure time.

    Figure 10.Comparison of raw single-molecule images before and after aberration correction, accompanied by the distribution of photons and background photons. (a) Snapshots of an unmodified 2D single-molecule image of Nup96-SNAP-AF647 labeled NPCs in U2OS cell samples refocused to 4 μm above the nominal focal plane. The inset shows a zoomed view of the fluorophore in the boxed area. (b) Snapshots of the 2D single-molecule image after aberration correction, with a significantly improved SNR in the image. The inset clearly shows the successful recovery of the distorted PSF. (c) The photon distribution analysis from 2000 frames of 2D single-molecule images before and after aberration correction reveals that 1644 single-molecule events with an average photon count of 335 were localized before correction, whereas 9680 single-molecule events with an average photon count of 711 were localized after correction. (d) The background photon distribution analysis from the same 2000 frames of 2D single-molecule images before and after aberration correction showed an average background photon count of 29 before correction and an average photon count of 25 after correction. The single-molecule raw images were acquired at 120 mW laser power and 20 ms exposure time.

    Correcting aberrations during refocusing with astigmatism PSF using a DM at various axial positions of the beads on the coverslip. (a) The beads are refocused using a DM to a position 2 μm above the nominal focal plane. (b) The pupil function for the refocused beads. (c) Fitted 21 Zernike coefficients for the beads stack. (d) CRLB for the experimental astigmatism PSF, with 3000 photons and 30 background photons used in the CRLB calculation. (e)–(h) The PSF shape, pupil function, Zernike coefficients, and localization precision after aberration correction, respectively. The CRLB exhibits significant improvements, particularly in the Z-axis localization precision.

    Figure 11.Correcting aberrations during refocusing with astigmatism PSF using a DM at various axial positions of the beads on the coverslip. (a) The beads are refocused using a DM to a position 2 μm above the nominal focal plane. (b) The pupil function for the refocused beads. (c) Fitted 21 Zernike coefficients for the beads stack. (d) CRLB for the experimental astigmatism PSF, with 3000 photons and 30 background photons used in the CRLB calculation. (e)–(h) The PSF shape, pupil function, Zernike coefficients, and localization precision after aberration correction, respectively. The CRLB exhibits significant improvements, particularly in the Z-axis localization precision.

    Experimental astigmatic 3D data of NPCs before and after aberration correction by refocusing. (a) Overview of 3D super-resolution image of Nup96-SNAP-AF647, refocused using a DM at 2 μm above the nominal focal plane. Reconstructed using the Cspline experimental PSF model shown in Fig. 11(a). (b) Magnified view of the area denoted by the box in (a). (c) Side-view cross section of the region denoted by the dashed line in (a). (d)–(f) 3D super-resolution image of NPCs after aberration correction at the refocused position. Significant improvements were observed in the reconstructed images, particularly in the XZ plane, where the double ring structure of NPCs became distinctly visible after aberration correction, as indicated by the white arrows in panels (c) and (f). The data were acquired from 50,000 frames with 120 mW laser power and 20 ms exposure time.

    Figure 12.Experimental astigmatic 3D data of NPCs before and after aberration correction by refocusing. (a) Overview of 3D super-resolution image of Nup96-SNAP-AF647, refocused using a DM at 2 μm above the nominal focal plane. Reconstructed using the Cspline experimental PSF model shown in Fig. 11(a). (b) Magnified view of the area denoted by the box in (a). (c) Side-view cross section of the region denoted by the dashed line in (a). (d)–(f) 3D super-resolution image of NPCs after aberration correction at the refocused position. Significant improvements were observed in the reconstructed images, particularly in the XZ plane, where the double ring structure of NPCs became distinctly visible after aberration correction, as indicated by the white arrows in panels (c) and (f). The data were acquired from 50,000 frames with 120 mW laser power and 20 ms exposure time.

    Optimization of tetrapod PSF using a DM for the beads on the coverslip. (a) The averaged tetrapod PSF for fluorescent beads on the coverslip, with PSFs at axial positions from −2500 to 2500 nm. (b) The pupil function for the beads. (c) Fitted 21 Zernike coefficients for the beads stack. (d) CRLB for the experimental tetrapod PSF, with 3000 photons and 30 background photons used in the CRLB calculation.

    Figure 13.Optimization of tetrapod PSF using a DM for the beads on the coverslip. (a) The averaged tetrapod PSF for fluorescent beads on the coverslip, with PSFs at axial positions from 2500 to 2500 nm. (b) The pupil function for the beads. (c) Fitted 21 Zernike coefficients for the beads stack. (d) CRLB for the experimental tetrapod PSF, with 3000 photons and 30 background photons used in the CRLB calculation.

    Whole-cell 3D super-resolution imaging of NPC using DMO tetrapod PSF. (a) Overview of the panoramic whole-cell 3D imaging of NPC using DMO tetrapod PSF without remote focusing. (b) FRC analysis of the regions indicated in (a). (c) Side-view cross section of the region denoted by the dashed line in (a). (d) Magnified view of the area denoted by the box in (a). (e) View of the bottom surface of the area in (d). (f) View of the top surface of the area in (d). The data were acquired from 100,000 frames with 120 mW laser power and 20 ms exposure time.

    Figure 14.Whole-cell 3D super-resolution imaging of NPC using DMO tetrapod PSF. (a) Overview of the panoramic whole-cell 3D imaging of NPC using DMO tetrapod PSF without remote focusing. (b) FRC analysis of the regions indicated in (a). (c) Side-view cross section of the region denoted by the dashed line in (a). (d) Magnified view of the area denoted by the box in (a). (e) View of the bottom surface of the area in (d). (f) View of the top surface of the area in (d). The data were acquired from 100,000 frames with 120 mW laser power and 20 ms exposure time.

    Comparison of in situ aberration pre-correction and actual aberration correction differences in various cells using single-molecule blinking data. (a) Pre-correction PSF derived from single-molecule blinking data of immunofluorescence-labeled TOM20 in COS7 cells at a depth of 2.5 μm. (b) Pupil function calculated from the in situ PSF of pre-correction. (c), (d) The PSF shape and its corresponding pupil function after in situ aberration correction. (e) The fitted 21 Zernike coefficients obtained from the in-situ single-molecule blinking data of pre-corrected cells, contrasting the differences between pre-corrected and actual corrected cells.

    Figure 15.Comparison of in situ aberration pre-correction and actual aberration correction differences in various cells using single-molecule blinking data. (a) Pre-correction PSF derived from single-molecule blinking data of immunofluorescence-labeled TOM20 in COS7 cells at a depth of 2.5 μm. (b) Pupil function calculated from the in situ PSF of pre-correction. (c), (d) The PSF shape and its corresponding pupil function after in situ aberration correction. (e) The fitted 21 Zernike coefficients obtained from the in-situ single-molecule blinking data of pre-corrected cells, contrasting the differences between pre-corrected and actual corrected cells.

    While using the beads PSF model, we imaged the sample in a refractive index matching buffer. The refractive index of this buffer was matched to that of the silicone immersion oil (1.406), ensuring a constant refractive index axially. When employing the in situ PSF model, we imaged the sample in a refractive index buffer that was not matched to the immersion oil (1.518). The sample was excited using a 642 nm laser with an intensity of 120–200 mW, and the 405 nm laser was adjusted according to the density of blinking single-molecule emitters. Typically, we acquired 100,000–300,000 frames per dataset with an exposure time of 15–20 ms.

    Data Processing and Analysis for Multi-Optical Section

    We have provided an open-source software for analyzing bead stacks and employed a GPU-based vectorial PSF fitter to fit them using MLE. The pupil function of average beads was decomposed into Zernike polynomials and we calculated the coefficients of these polynomials. Additionally, we used an open-source deep learning-based algorithm to analyze the single-molecule data of cell samples encoded with DMO PSF. This open-source software can be accessed in Ref. [42]. Furthermore, for in situ single-molecule blinking data extraction and in situ PSF modeling, as well as the computation of Zernike polynomial coefficients, we utilized our recently published open-source software, which can be accessed in Ref. [43]. We extracted single-molecule images from each optical section, which were acquired in a cyclic manner from the image stacks. These images were then merged and analyzed to generate a localization list for each section.

    Despite having implemented an axial focusing system in the hardware setup (Fig. 8 in Appendix B), the overall lateral drift of the sample and the accuracy of axial refocusing positions can result in slight misalignments of adjacent optical sections. This might lead to a deterioration of the resolution and distortions of the super-resolution images during the merging of sections. When performing drift correction, all sections were corrected using the drift values from the first section to prevent stitching artifacts and distortions in super-resolution images caused by computational offsets across different sections. The drift correction was performed using the redundant cross-correlation method [44] in SMAP [34]. When stitching data from neighboring optical sections, we employed the same redundant cross-correlation method as in drift correction. We scanned the entire cell sample in the axial direction with 1 μm step sizes, maintaining an optical section thickness of 1.2 μm. We ensured overlapping regions between adjacent optical sections to facilitate the use of the redundant cross-correlation method. Initially, we aligned neighboring optical sections in the XY direction, followed by alignment in the XZ direction. This approach enabled precise 3D alignment and stitching of the entire cell data.

    Biological Sample Preparation

    Cell Culture. U2OS cells (Nup96-SNAP No. 300444, Cell Line Services) were grown in DMEM (catalog No. 10569, Gibco) containing 10% (volume fraction) fetal bovine serum (FBS; catalog No. 10099-141C, Gibco), 100 U/mL penicillin and 100 μg/mL streptomycin (PS; catalog No. 15140-122, Gibco), and 1× MEM NEAA (catalog No. 11140-050, Gibco). COS-7 cells (catalog No. 100040, BNCC) were grown in DMEM containing FBS and PS. Cells were cultured in a humidified atmosphere with 5% CO2 at 37°C and passaged every 2 or 3 days. Before cell plating, high-precision 25-mm-round glass coverslips (catalog No. CG15XH, Thorlabs) were cleaned by sequential sonication in 1 M (1 mol/L) potassium hydroxide (KOH), Milli-Q water, and ethanol, and finally irradiated under ultraviolet light for 30 min. For super-resolution imaging, U2OS and COS-7 cells were cultured on the clean coverslips for 2 days with a confluency of 80%.

    SNAP-Tag Labeling of Nup96. To label Nup96, U2OS-Nup96-SNAP cells were prepared as previously reported [33]. Cells were prefixed in 2.4% (mass-to-volume ratio) paraformaldehyde (PFA) for 30 s, permeabilized in 0.4% (volume fraction) Triton X-100 for 3 min, and subsequently fixed in 2.4% PFA for 30 min. The buffer used for fixing and permeabilization was preheated to 37°C before use. Then, cells were quenched in 0.1 M NH4Cl for 5 min and washed twice with PBS. To decrease unspecific binding, cells were blocked for 30 min with Image-iT FX signal enhancer (catalog No. I36933, Invitrogen). For labeling, cells were incubated in dye solution (1 μM SNAP-tag ligand BG-AF647 (catalog No. S9136S, New England Biolabs), 1 mM DTT (catalog No. 1111GR005, BioFroxx), and 0.5% (mass-to-volume ratio) bovine serum albumin (BSA) in PBS for 2 h and washed three times in PBS for 5 min each to remove excess dyes. Lastly, cells were postfixed with 4% PFA for 10 min, washed with PBS three times for 3 min each, and stored at 4°C until imaged.

    Microtubule Labeling. Microtubule samples were prepared as previously described [31]. COS-7 cells were prefixed with 0.3% (volume fraction) glutaraldehyde (GA) and 0.25% Triton X-100 in the cytoskeleton buffer [CB (10 mM MES, 5 mM glucose, 150 mM NaCl, 5 mM MgCl2, 5 mM EGTA, pH 6.1), preheated to 37°C before using] for 1–2 min. Waste liquor was removed and fixed with 2% GA in CB for 10 min (preheated to 37°C before using). We removed waste liquor, quenched the samples in 0.1% (mass-to-volume ratio) NaBH4 (0.01 g in 10 mL PBS) for 7 min, and washed three times for 5 min each with PBS. Then incubate the samples in the permeabilization buffer (0.1% Triton X-100 in 3% BSA) for 10 min. After washing three times for 5 min each with PBS, cells were stained by goat anti-mouse β-tubulin (catalog No. T4026, Sigma, 2.4 mg/mL) with 1:1000 dilution in 3% BSA for 1 h and washed three times for 5 min each with PBS. Cells were then stained with the corresponding secondary antibodies conjugated with AF647 (catalog No. 2289596, Invitrogen, 2 mg/mL) with 1:1000 dilution in 3% BSA for 1 h and washed with PBS. Finally, cells were postfixed with 4% PFA for 10 min, washed with PBS three times for 5 min each, and stored in PBS at 4°C.

    Mitochondria Labeling. Mitochondrial samples were prepared as previously described [45]. COS-7 cells were fixed with 4% PFA (preheated to 37°C before using) in PBS for 12 min, incubated in permeabilization buffer [0.3% CA-630 (catalog No. I8896, Sigma), 0.05% TX-100, 0.1% BSA, and 1× PBS] for 3 min, and then quenched in 0.1 M NH4Cl for 5 min. After being washed three times for 5 min each with PBS, cells were blocked in 3% BSA for 60 min. For labeling, cells were stained by rabbit anti-Tom20 (catalog No. ab78547, Abcam, 1 mg/ mL) with 1:1000 dilution in 3% BSA for 2 h and washed three times for 5 min each with PBS. Cells were then stained with the corresponding secondary antibodies conjugated with AF647 (catalog No. A21245, Invitrogen, 2 mg /mL) with 1:2000 dilution in 3% BSA for 2 h, and washed three times for 5 min each with PBS. Finally, cells were postfixed with 4% PFA for 10 min, washed with PBS three times for 5 min each, and stored in PBS at 4°C.

    Imaging Buffer. Samples were imaged in a refractive index matching buffer, including 50 mM Tris-HCl (pH 8.0), 10 mM NaCl, 10% (mass-to-volume ratio) glucose, 0.5 mg/mL glucose oxidase (G7141, Sigma), 40 μg/mL catalase (C100, Sigma), 35 mM cysteamine, and 28.5% (volume fraction) 2,2’-thiodiethanol (166782, Sigma). The refractive index of the final imaging buffer is 1.406.

    Beads Preparation. We first diluted the 100-nm-diameter crimson beads (custom-designed, Invitrogen) to 1:40,000 in Milli-Q water and vortexed the mixture for 3–5 min. Then, we pipetted 40 μL of 1 M MgCl2 in the center of the 25-mm-diameter coverslip (cleaning protocol refers to the cell culture section) and mixed it with 360 μL of the diluted bead solution on the coverslip. The mixture was incubated for 5 min, washed three times with Milli-Q water, and stored in Milli-Q water at 4°C in the dark.

    APPENDIX B: DETAILED LAYOUT OF THE OPTICAL SETUP

    Figure 8 shows the detailed optical setup of the remote focusing system.

    APPENDIX C: CRLB CALCULATIONS BEFORE AND AFTER ABERRATION CORRECTION OF THE BEADS ON THE COVERSLIP DURING REFOCUSING WITH DMO PSF USING A DM, SIMULATING THE TYPICAL PHOTON FLUX OF FLUORESCENT PROTEINS

    Figure 9 shows the localization precision of the calculated CRLB for beads data before and after aberration correction under the typical photon flux of fluorescent proteins.

    APPENDIX D: COMPARISON OF RAW SINGLE-MOLECULE IMAGES BEFORE AND AFTER ABERRATION CORRECTION, ACCOMPANIED BY THE DISTRIBUTION OF PHOTONS AND BACKGROUND PHOTONS

    Figure 10 shows the distribution of photons and background photons in the raw single-molecule data before and after the residual aberration correction by applying the defocus wavefront with nucleoporin Nup96 in U2OS cells.

    APPENDIX E: CORRECTING ABERRATIONS DURING REFOCUSING WITH ASTIGMATISM PSF USING A DM AT VARIOUS AXIAL POSITIONS OF THE BEADS ON THE COVERSLIP

    Figure 11 shows the astigmatism PSF before and after aberration correction, where beads are refocused using a DM at 2 μm above the nominal focal plane.

    APPENDIX F: EXPERIMENTAL ASTIGMATIC 3D DATA OF NPCS BEFORE AND AFTER ABERRATION CORRECTION BY REFOCUSING

    Figure 12 shows the experimental astigmatic 3D data of NPCs before and after aberration correction, refocused using a DM at 2 μm above the nominal focal plane.

    APPENDIX G: OPTIMIZATION OF TETRAPOD PSF USING A DM FOR THE BEADS ON THE COVERSLIP

    Figure 13 shows the optimization of a 6 μm large DOF tetrapod PSF using a DM for the beads on the coverslip.

    APPENDIX H: WHOLE-CELL 3D SUPER-RESOLUTION IMAGING OF NPC USING DMO TETRAPOD PSF

    Figure 14 shows whole-cell 3D super-resolution imaging of NPC using DMO tetrapod PSF without remote focusing.

    APPENDIX I: COMPARISON OF IN SITU ABERRATION PRE-CORRECTION AND ACTUAL ABERRATION CORRECTION DIFFERENCES IN VARIOUS CELLS USING SINGLE-MOLECULE BLINKING DATA

    Figure 15 shows a comparison of in situ aberrations between actual imaging cells and pre-correction cells, using single-molecule blinking data of immunofluorescence-labeled TOM20 in COS7 cells at a depth of 2.5 μm.

    References

    [1] L. Schermelleh, A. Ferrand, T. Huser. Super-resolution microscopy demystified. Nat. Cell Biol., 21, 72-84(2019).

    [2] M. Lelek, M. T. Gyparaki, G. Beliu. Single-molecule localization microscopy. Nat. Rev. Methods Primers, 1, 39(2021).

    [3] B. Huang, W. Wang, M. Bates. Three-dimensional super-resolution imaging by stochastic optical reconstruction microscopy. Science, 319, 810-813(2008).

    [4] F. Huang, G. Sirinakis, E. S. Allgeyer. Ultra-high resolution 3D imaging of whole cells. Cell, 166, 1028-1040(2016).

    [5] S. Liu, H. Huh, S. H. Lee. Three-dimensional single-molecule localization microscopy in whole-cell and tissue specimens. Annu. Rev. Biomed. Eng., 22, 155-184(2020).

    [6] E. Betzig, G. H. Patterson, R. Sougrat. Imaging intracellular fluorescent proteins at nanometer resolution. Science, 313, 1642-1645(2006).

    [7] S. T. Hess, T. P. Girirajan, M. D. Mason. Ultra-high resolution imaging by fluorescence photoactivation localization microscopy. Biophys. J., 91, 4258-4272(2006).

    [8] M. J. Rust, M. Bates, X. Zhuang. Sub-diffraction-limit imaging by stochastic optical reconstruction microscopy (STORM). Nat. Methods, 3, 793-795(2006).

    [9] B. Huang, S. A. Jones, B. Brandenburg. Whole-cell 3D STORM reveals interactions between cellular structures with nanometer-scale resolution. Nat. Methods, 5, 1047-1052(2008).

    [10] A. M. Rozario, A. Morey, C. Elliott. 3D single molecule super-resolution microscopy of whole nuclear lamina. Front. Chem., 10, 863610(2022).

    [11] F. Xu, D. Ma, K. P. MacPherson. Three-dimensional nanoscopy of whole cells and tissues with in situ point spread function retrieval. Nat. Methods, 17, 531-540(2020).

    [12] M. F. Juette, T. J. Gould, M. D. Lessard. Three-dimensional sub-100 nm resolution fluorescence microscopy of thick samples. Nat. Methods, 5, 527-529(2008).

    [13] B. Hajj, J. Wisniewski, M. El Beheiry. Whole-cell, multicolor superresolution imaging using volumetric multifocus microscopy. Proc. Natl. Acad. Sci. USA, 111, 17480-17485(2014).

    [14] A. Aristov, B. Lelandais, E. Rensen. ZOLA-3D allows flexible 3D localization microscopy over an adjustable axial range. Nat. Commun., 9, 2409(2018).

    [15] E. Nehme, D. Freedman, R. Gordon. DeepSTORM3D: dense 3D localization microscopy and PSF design by deep learning. Nat. Methods, 17, 734-740(2020).

    [16] S. Fu, W. Shi, T. Luo. Field-dependent deep learning enables high-throughput whole-cell 3D super-resolution imaging. Nat. Methods, 20, 459-468(2023).

    [17] S. Fu, M. Li, L. Zhou. Deformable mirror based optimal PSF engineering for 3D super-resolution imaging. Opt. Lett., 47, 3031-3034(2022).

    [18] Q. Zhang, Q. Hu, C. Berlage. Adaptive optics for optical microscopy. Biomed. Opt. Express, 14, 1732-1756(2023).

    [19] S. Y. Kang, M. Duocastella, C. B. Arnold. Variable optical elements for fast focus control. Nat. Photon., 14, 533-542(2020).

    [20] V. Navikas, A. C. Descloux, K. S. Grussmayer. Adaptive optics enables multimode 3D super-resolution microscopy via remote focusing. Nanophotonics, 10, 2451-2458(2021).

    [21] Z. Qin, Z. She, C. Chen. Deep tissue multi-photon imaging using adaptive optics with direct focus sensing and shaping. Nat. Biotechnol., 40, 1663-1671(2022).

    [22] M. Zurauskas, O. Barnstedt, M. Frade-Rodriguez. Rapid adaptive remote focusing microscope for sensing of volumetric neural activity. Biomed. Opt. Express, 8, 4369-4379(2017).

    [23] J. Cui, R. Turcotte, N. J. Emptage. Extended range and aberration-free autofocusing via remote focusing and sequence-dependent learning. Opt. Express, 29, 36660-36674(2021).

    [24] M. Bathe-Peters, P. Annibale, M. J. Lohse. All-optical microscope autofocus based on an electrically tunable lens and a totally internally reflected IR laser. Opt. Express, 26, 2359-2368(2018).

    [25] J. Basumatary, N. Baro, P. Joshi. Scanning single molecule localization microscopy (scanSMLM) for super-resolution volume imaging. Commun. Biol., 6, 1050(2023).

    [26] H. Dibaji, A. K. N. Shotorban, M. Habibi. Axial de-scanning using remote focusing in the detection arm of lighft-sheet microscopy. bioRxiv(2023).

    [27] P. Zhang, D. Ma, X. Cheng. Deep learning-driven adaptive optics for single-molecule localization microscopy. Nat. Methods, 20, 1748-1758(2023).

    [28] R. P. J. Nieuwenhuizen, K. A. Lidke, M. Bates. Measuring image resolution in optical nanoscopy. Nat. Methods, 10, 557-562(2013).

    [29] J. Antonello. Interferometric calibration of a deformable mirror(2020).

    [30] Y. Shechtman, S. J. Sahl, A. S. Backer. Optimal point spread function design for 3D imaging. Phys. Rev. Lett., 113, 133902(2014).

    [31] Y. M. Li, M. Mund, P. Hoess. Real-time 3D single-molecule localization using experimental point spread functions. Nat. Methods, 15, 367-369(2018).

    [32] Y. Li, W. Shi, S. Liu. Global fitting for high-accuracy multi-channel single-molecule localization. Nat. Commun., 13, 3133(2022).

    [33] J. V. Thevathasan, M. Kahnwald, K. Cieslinski. Nuclear pores as versatile reference standards for quantitative superresolution microscopy. Nat. Methods, 16, 1045-1053(2019).

    [34] J. Ries. SMAP: a modular super-resolution microscopy analysis platform for SMLM data. Nat. Methods, 17, 870-872(2020).

    [35] Y. M. Li, E. Buglakova, Y. D. Zhang. Accurate 4Pi single-molecule localization using an experimental PSF model. Opt. Lett., 45, 3765-3768(2020).

    [36] J. Y. Wang, E. S. Allgeyer, G. Sirinakis. Implementation of a 4Pi-SMS super-resolution microscope. Nat. Protoc., 16, 677-727(2021).

    [37] S. Liu, J. Chen, J. Hellgoth. Universal inverse modelling of point spread functions for SMLM localization and microscope characterization. bioRxiv(2023).

    [38] M. E. Siemons, N. A. K. Hanemaaijer, M. H. P. Kole. Robust adaptive optics for localization microscopy deep in complex tissue. Nat. Commun., 12, 3407(2021).

    [39] S. Park, Y. Jo, M. Kang. Label-free adaptive optics single-molecule localization microscopy for whole zebrafish. Nat. Commun., 14, 4185(2023).

    [40] J. Deschamps, J. Ries. EMU: reconfigurable graphical user interfaces for Micro-Manager. BMC Bioinform., 21, 456(2020).

    [41] H. Pinkard, N. Stuurman, I. E. Ivanov. Pycro-Manager: open-source software for customized and reproducible microscope control. Nat. Methods, 18, 226-228(2021).

    [42] https://github.com/Li-Lab-SUSTech/FD-DeepLoc. https://github.com/Li-Lab-SUSTech/FD-DeepLoc

    [43] https://github.com/ries-lab/uiPSF. https://github.com/ries-lab/uiPSF

    [44] Y. N. Wang, J. Schnitzbauer, Z. Hu. Localization events-based sample drift correction for localization microscopy with redundant cross-correlation algorithm. Opt. Express, 22, 15982-15991(2014).

    [45] Y. Zhang, L. K. Schroeder, M. D. Lessard. Nanoscale subcellular architecture revealed by multicolor three-dimensional salvaged fluorescence imaging. Nat. Methods, 17, 225-231(2020).

    [46] https://doi.org/10.5281/zenodo.10580127. https://doi.org/10.5281/zenodo.10580127

    Wei Shi, Yingchuan He, Jianlin Wang, Lulu Zhou, Jianwei Chen, Liwei Zhou, Zeyu Xi, Zhen Wang, Ke Fang, Yiming Li. Aberration correction for deformable-mirror-based remote focusing enables high-accuracy whole-cell super-resolution imaging[J]. Photonics Research, 2024, 12(4): 821
    Download Citation