Super-resolution fluorescence microscopy is an important tool in biomedical research for its ability to discern features smaller than the diffraction limit. However, due to its difficult implementation and high cost, the universal application of super-resolution microscopy is not feasible. In this paper, we propose and demonstrate a new kind of super-resolution fluorescence microscopy that can be easily implemented and requires neither additional hardware nor complex post-processing. The microscopy is based on the principle of stepwise optical saturation (SOS), where steps of raw fluorescence images are linearly combined to generate an image with a -fold increase in resolution compared with conventional diffraction-limited images. For example, linearly combining (scaling and subtracting) two images obtained at regular powers extends resolution by a factor of beyond the diffraction limit. The resolution improvement in SOS microscopy is theoretically infinite but practically is limited by the signal-to-noise ratio. We perform simulations and experimentally demonstrate super-resolution microscopy with both one-photon (confocal) and multiphoton excitation fluorescence. We show that with the multiphoton modality, the SOS microscopy can provide super-resolution imaging deep in scattering samples.
Super-resolution fluorescence microscopy by stepwise optical saturation
Yide Zhang, Prakash D. Nallathamby, Genevieve D. Vigil, Aamir A. Khan, Devon E. Mason, Joel D. Boerckel, Ryan K. Roeder, and Scott S. Howard
Email: firstname.lastname@example.org, email@example.com
Super-resolution fluorescence microscopy techniques, such as stimulated emission depletion (STED) microscopy [1, 2] and its related reversible saturable optical fluorescence transitions (RESOLFT) microscopy , photoactivated localization microscopy (PALM) , stochastic optical reconstruction microscopy (STORM) , and structured illumination microscopy (SIM) [6, 7], have enabled a dramatic development in modern biology by being able to discern fluorescent molecules or features that are closer together than the diffraction limit [8, 9]. Many super-resolution techniques, however, only work well on thin and nearly transparent samples . Super-resolution imaging deep in scattering samples is still a challenge, since optical aberrations and scattering severely degrade resolution and signal-to-noise ratio (SNR), so obtaining even diffraction-limited performance is difficult . Meanwhile, the implementation of most super-resolution techniques is difficult and expensive, which hinders the universal application of super-resolution microscopy in many labs.
To increase the imaging depth, some super-resolution methods have been combined with multiphoton microscopy (MPM) [12, 11, 13, 14], an imaging technique widely used in biomedical research for its inherent 3D resolution, deep penetration, and minimal phototoxicity [15, 16, 17]. On the other hand, to reduce the complexity and cost of super-resolution microscopy, novel techniques such as saturated excitation (SAX) microscopy and its variants [18, 19, 13, 20, 21] and fluorescence emission difference (FED) microscopy  have been developed, whose implementations are easier and cheaper compared to traditional super-resolution methods. A super-resolution technique with both deep penetration and easy implementation would be desirable and widely utilized.
In this paper, we propose and demonstrate a new kind of super-resolution fluorescence microscopy using the principle of stepwise optical saturation (SOS), which is designed to work with both confocal and MPM modalities and to permit deep penetration. In its simplest form, two-step SOS, two conventional fluorescence images are linearly combined to produce a super-resolution image with a -fold increase in spatial resolution. In general, an -step microscopy uses the linear combination of raw images to generate a super-resolved SOS image, which has a -fold increase in spatial resolution compared to diffraction-limited images. The improvement in resolution is theoretically infinite, but practically, it is limited by the SNR performance of the resulting SOS image. We perform simulations and experiments with both one-photon excitation fluorescence (1PEF) and two-photon excitation fluorescence (2PEF) to validate the resolution improvement. Similar to SAX microscopy, fluorophore saturation is used in SOS microscopy but the excitation power used to cause the saturation is relatively weak, so effects like photobleaching and photodamage are reduced. We show that the implementation of SOS microscopy is straightforward and requires neither additional hardware nor complex post-processing, which are usually required in SAX microscopy. With the SOS methods, it is easy to produce super-resolution images using a conventional fluorescence microscope.
2 Principle of Stepwise Optical Saturation
The principle of SOS microscopy is based on a two-level fluorophore model depicted in a Jablonski diagram in Fig. 1(a), similar to our previous work in . We denote the populations of the ground and excited singlet states as and , respectively, as functions of time. Note that , where is the concentration of the fluorophore. To make this model applicable to both one-photon and multiphoton excitation, the excitation rate is written as , where is the incident photon flux, is the number of excitation photons needed for a fluorophore to emit one photon ( for 1PEF, for 2PEF, etc.), is the cross-section for -photon excitation, and is the pulse gain factor which accounts for the temporal pulse profile of the excitation [24, 25]. The fluorescence rate is , where is the fluorescence lifetime. This two-level model is an empirical abstraction of the dynamics of the singlet and triplet states and the intersystem crossing [26, 27]. A rate equation describing this model can be written as
We can convert and to experimentally measurable excitation irradiance, , and fluorescence intensity, , by and , respectively, where is Planck’s constant, is the velocity of light, is the excitation wavelength, is the fluorescence detection efficiency, and is the observation time. Denoting and , from Eq. (1), we have
The SOS technique only requires the linear combination of conventional fluorescence intensity images. Therefore it is valid to solely consider steady-state solutions of the fluorophore model. Dropping the temporal dependency, the steady-state solution of Eq. 2 is
Although higher excited singlet and triplet states and photobleaching are not considered in this two-level system, the model can precisely describe the fluorophore saturation behavior. Fig. 1(b) plots the simulated fluorescence-excitation relation described by this model for the 1PEF case. In the range of low excitation irradiance, the fluorescence-excitation relation is linear and unsaturated, while in the range of high excitation irradiance, the relation becomes nonlinear and saturated. Note that the saturation phenomenon can happen without a high excitation intensity. As shown in the inset of Fig. 1(b), even for a moderate excitation irradiance, the simulated curve can deviate from the unsaturated linear curve. We call this phenomenon “weak saturation” and it is the basis for the SOS technique. Because of the “weak saturation” phenomenon, SOS microscopy is able to produce super-resolution images at a relatively low excitation intensity. Under the “weak saturation” condition, denoting , Eq. 3 can be Taylor expanded to
An -step SOS microscopy needs fluorescence images to be collected. Here we consider one-dimensional spatial dependence of excitation and fluorescence. The extension to a multi-dimensional case is trivial. A Gaussian excitation profile is assumed in the focus, with the focal irradiance and radius . In the case of -photon excitation, this profile is effectively powered to the -th. A subscript is added to the excitation and fluorescence intensities of the -th step image among the images. We can separate the spatial dependency by denoting . Then for the -th step, we have and Eq. 4 can be written as
In Eq. 5, high powers of , such as , , etc., represent components with higher spatial frequency, while is the diffraction-limited component. With the definition of , an -th order power component, , has a -fold increase in spatial resolution. Due to the magnitude difference among each components, the spatial resolution of is dominated by the lowest power of .
The idea of an -step SOS microscopy is to eliminate the lowest powers of by the linear combination of steps of conventional images obtained at different excitation intensities. We assume the excitation intensities of each step follow and they stay in the “weak saturation” region. The resulting SOS image is , where the coefficients are chosen such that the lowest power of in is . In other words, an -step SOS image will have a -fold increase in spatial resolution.
The coefficients of the linear combination, , , , , of the images with excitation intensities, , , , , are calculated by solving variables from equations, where the variables are the coefficients and the equations are the constraints that the lowest powers of in the resulting -step SOS image, , are eliminated. Since this is an under-constrained problem, we set without loss of generality. Table 1 presents the coefficients for the linear combination of SOS methods from two-step to four-step.
For instance, in two-step SOS, we need two images, and , obtained at excitation intensities, and , respectively. With the power series in Eq. 5 and the two-step coefficients, and , in Tab. 1, the resulting two-step SOS image is , which is dominated by the component . Compared to the diffraction-limited component , has a -fold increase in spatial resolution. In general, for an -step SOS, a super-resolved image with a -fold resolution improvement compared to the diffraction limit can be achieved by the linear combination with the coefficients in Tab. 1. Although Tab. 1 only provides SOS coefficients from two-step to four-step due to limited space for this article, the coefficients for any number of steps can be calculated using the methods described above. Alternatively, by inspecting the mathematical expressions in Tab. 1, a rule about the formation of SOS coefficients can be observed. Specifically, for two adjacent steps and , where , if , the coefficient in -step SOS is the product of in -step SOS and a factor ; if , then in -step SOS is
Figure 2 graphically illustrates how super-resolution images are obtained in SOS microscopy. Generally, in an -step microscopy, steps of diffraction-limited raw images [here we only plot point spread functions (PSFs) for illustration] obtained at different excitation intensities, , , , , are linearly combined to produce a super-resolution SOS image with a -fold increase in spatial resolution (a PSF -fold narrower than a diffraction-limited PSF). The linear combination coefficients, , , , , are functions of the intensities and are presented in Tab. 1. For instance, in two-step SOS, we need two images, and , obtained at excitation intensities, and , respectively. With the two-step coefficients, and , in Tab. 1, the resulting two-step SOS image is dominated by the component . Compared to the diffraction-limited component , has a -fold increase in spatial resolution. Similarly, three-step SOS needs three raw images and the resulting SOS image, , is dominated by the components , which has a -fold increase in spatial resolution. The same principle applies to -step SOS, which is able to generate a super-resolution image with a -fold increase in spatial resolution.
3 Numerical Results
The super-resolution capability of SOS microscopy is first investigated using numerical simulations. The simulation is based on the two-level model described above and the parameters are given in Section 5. Here several diffraction-limited PSFs obtained at different excitation intensities are used to generate SOS images from two-step to six-step. The resulting PSFs of the SOS microscopy for both 1PEF and 2PEF are shown in Fig. 3. For the 1PEF case, the conventional diffraction-limited PSF has a full-width at half-maximum (FWHM) of 228.9 nm, while the two-step to six-step SOS’s FWHMs are 162.1 nm, 132.5 nm, 114.9 nm, 103.1 nm, and 94.5 nm, corresponding to a 1.41, 1.73, 2.00, 2.22, and 2.42 folds increase in resolution, respectively. For the 2PEF case, the conventional diffraction-limited FWHM is 265.3 nm, while the FWHMs of the two-step to six-step SOS are 187.7 nm, 153.5 nm, 133.3 nm, 119.5 nm, and 109.5 nm, which, identically, correspond to a 1.41, 1.73, 2.00, 2.22, and 2.42 folds increase in resolution, respectively. It is clearly shown that the resolution improvements for both 1PEF and 2PEF are equivalent, and for an -step SOS, whether it is 1PEF or 2PEF, the resolution improvement is exactly -fold.
Although the principle of SOS microscopy is derived with one-dimensional objects. The extension to two- or multi-dimensional cases is trivial. Here we apply the SOS methods to a two-dimensional simulation target, as shown in Fig. 4. The target is a star-like object consisting of several lines, where the closet distance of the lines is 140 nm. The PSFs in Fig. 3 are used as kernels to simulate the diffraction-limited and SOS images for both 1PEF and 2PEF cases. The diffraction-limited images for both cases are blurred and the center ends of the lines cannot be differentiated. With the help of the SOS microscopy, super-resolved images of the object can be obtained. In the simplest two-step SOS, an obvious resolution improvement can be seen for both 1PEF and 2PEF. As the number of SOS steps increases, the resolution improvement is more significant. The performance of the images in Fig. 4 matches well with the FWHM results in Fig. 3, as the center ends of the lines in the object are successfully differentiated in three-step to six-step 1PEF SOS [Figs. 4(c-f)] and four-step to six-step 2PEF SOS [Figs. 4(j-l)], whose FWHMs are smaller than the closet distance of the lines (140 nm) in the object.
4 Experimental Results
We now present the experimental validation of the proposed SOS methods. Although SOS microscopy can provide a -fold increase in spatial resolution with the linear combination of steps of conventional images, and theoretically this resolution improvement can be infinite, practically, however, an experimental SOS image generated from more than two steps of raw images is generally unacceptable due to its extremely poor SNR. For example, for two steps of 1PEF images obtained with irradiances and , the two-step coefficients, based on Tab. 1, are and ; for three-step 1PEF images with irradiances , , and , the coefficients are , , and , respectively. Since the excitation irradiances for each image are approximate, if we assume the fluorescence intensities of raw images are the same and equal to , we can estimate the intensity of the resulting SOS image after the linear combination, . Hence, the intensities of two-step and three-step SOS images will be and , respectively, which are one order and two orders of magnitude lower compared with the raw intensity . Considering the shot noise and the propagation of errors in SOS algorithms , the SNR of a two-step SOS image will be an order of magnitude lower than a conventional image, while for SOS images generated from more than two steps, the SNR will be at least two orders of magnitude lower. Therefore in the experiments we only present confocal (1PEF) and multiphoton (2PEF) two-step SOS images for their satisfactory SNR performance and easy implementation. Figure 5 shows a representative example of how two raw images of Alexa Fluor 488 phalloidin labeled F-actin in fixed cells are obtained and processed to generate the super-resolved two-step SOS image. The raw images in Figs. 5(a) and 5(b) were obtained with Nikon’s commercial software. These two images were then imported into Matlab and converted to matrices and . Based on the knowledge of the excitation powers of these two images, the corresponding two-step SOS image was computed using the coefficients shown in Tab. 1. In the 2PEF two-step SOS case of Fig. 5, we have , mW, and mW, so the SOS coefficients are and and the resulting SOS image is . The resolution improvement of the two-step SOS image is appreciable. The materials and methods used in the experiments in this Section are described in Section 5.
Next, we imaged a phantom sample consisting of subdiffractive fluorescent nanoparticles (NPs) embedded in a hydrogel matrix (0.25%, ~55 nm RITC--SiEDTA NPs dispersed in 1% agarose gel). The conventional and two-step SOS images of the sample obtained with 1PEF and 2PEF are presented in Fig. 6. Compared to conventional, diffraction-limited 1PEF and 2PEF images in Figs. 6(a) and 6(d), the 1PEF and 2PEF two-step SOS images in Figs. 6(b) and 6(e) are super-resolved where NPs at different sites are better differentiated. The intensity profiles of a cluster of NPs obtained with 1PEF and 2PEF are fitted to Gaussian curves and plotted in Figs. 6(c) and 6(f), respectively. For the 1PEF case, the diffraction-limited Gaussian curve has a FWHM of 442.2 nm, while the two-step SOS’s FWHM is 325.1 nm, corresponding to a 1.36-fold increase in resolution. For the 2PEF case, the diffraction-limited FWHM is 636.9 nm, while the FWHM of the two-step SOS is 474.9 nm, which corresponds to a 1.34-fold increase in resolution. The resolution improvements in both 1PEF and 2PEF two-step SOS microscopy are close to the theoretical value, . The discrepancy, which is more severe for the 2PEF case, is due to the limited SNR of the experimentally obtained images.
We further studied the resolution improvement enabled by two-step SOS microscopy using a standard biological test slide [F-actin labeled with Alexa Fluor 488 phalloidin in fixed BPAE cells, FluoCells prepared slide #1, F36924] . Figure 7 shows the conventional and SOS images of the F-actin filaments with both 1PEF and 2PEF. By comparing the conventional [Figs. 7(a) and 7(d)] with the two-step SOS [Figs. 7(b) and 7(e)] images, one can see the significant sharpening of the actin structures enabled by the super-resolution SOS method in both 1PEF and 2PEF imaging. The intensity profiles of a few closely located filaments shown in Figs. 7(c) and 7(f) demonstrate the finer apparent width and better resolution of the filaments in two-step SOS images. In conventional 1PEF and 2PEF images, the differentiation of these filaments is not possible.
Having verified the resolution improvement with a standard sample, we then validated the SOS methods using a biological sample [F-actin labeled with Alexa Fluor 488 phalloidin in fixed endothelial colony forming cells (ECFCs)] prepared in our labs. The conventional and two-step imaging results for 1PEF and 2PEF are shown in Fig. 8. Since the fluorophore concentration used to stain this sample is lower than the one used in the standard test slide (Fig. 7), the fluorophores, though successfully bind to F-actin, look like subdiffractive beads when imaged. The improvement in resolution and imaging quality of two-step SOS microscopy [Figs. 8(b) and 8(e)] can be easily seen compared with the conventional diffraction-limited 1PEF and 2PEF images [Figs. 8(a) and 8(d)]. The intensity profiles in Figs. 8(c) and 8(f) show that, in both 1PEF and 2PEF imaging, the neighboring conventionally unresolvable fluorophores can be clearly resolved using two-step SOS microscopy.
Since two-step SOS microscopy was able to generate super-resolution multiphoton (2PEF) images, we then investigated the advantages of 2PEF, such as deep penetration and confined photodamage, in addition to super-resolution, in 2PEF two-step SOS microscopy. We used the same subdiffractive fluorescent NPs used in Fig. 6 (0.25%, ~55 nm RITC--SiEDTA NPs) but dispersed in a scattering phantom also comprising polystyrene beads within the agarose gel to mimic the optical scattering inside a brain tissue . The experimental 2PEF conventional and two-step SOS images at various depths from the coverslip surface are shown in Fig. 9. By virtue of the deep penetration of 2PEF, for a diffraction-limited 2PEF image obtained at each depth (from 0 m to 75 m) in a scattering phantom, we are able to extract a corresponding super-resolved two-step SOS image where conventionally unresolvable NPs can be clearly resolved. In principle, as long as a raw fluorescence image with an acceptable SNR can be obtained, a corresponding super-resolution two-step SOS image can be generated, regardless of the experimental conditions such as imaging depth.
5 Materials and Experimental Methods
Alexa Fluor-488 conjugated phalloidin (Invitrogen), Ammonium hydroxide (, 20%, VWR), (3-aminopropyl) triethoxysilane (APTES, , ¿98%, Sigma, St. Louis, MO), Aqua-Mount (Lerner laboratories), endothelial growth medium (EBM-2, Lonza), ethanol (200 proof, VWR), fetal bovine serum (FBS, GE life sciences), 10% w/v 1.1 micron Latex beads, polystyrene (PS beads, Sigma), 1000U/mL penicillin-streptomyocin (Corning), rhodamine B isothiocyanate (RITC, , Sigma), tetraethylorthosilicate (TEOS, , ¿99%, Sigma), triton-x-100 (, Sigma-Aldrich) and N-(trimethoxysillylpropyl)ethylenediaminetriacete, sodium salt, 30% in water (SiEDTA, , GELEST) were all used as-received.
5.2 Dispersion of 55 nm Fluorescent Silica Nanoparticles in Agarose Gel and Brain Optical Phantoms
A stock solution of silanized RITC was obtained by reacting RITC with APTES in a 1:1 mole ratio, overnight (RITC-APTES). Fluorescent silica nanoparticles (NPs) were synthesized by magnetically stirring 10 ml of ethanol and 2 ml of water in a 20 mL glass vial. 0.3 mL of TEOS and 0.075 mL of RITC-APTES were added sequentially and the solution was left stirring. 0.6 mL of 20% w/v was then added rapidly to the ethanolic solution of TEOS/RITC-APTES. The mixture was stirred for 12 h at room temperature, in the dark. The fluorescent silica particles (RITC- NPs) were recovered by centrifugation at 12000 relative centrifugal force (rcf) for 45 min, discarding the supernatant and resuspending the pellet in ethanol via ultrasonication. The as-prepared RITC- NPs were amphiphilic and could be dispersed in either ethanol or water. 10 mL of a 1% w/v solution of RITC- NPs in 80/20 ethanol/water solution was vigorously stirred with 0.1 mL of 30% w/v SiEDTA solution for 24 h in the dark, modifying previously established methods [31, 32, 33]. The silanized RITC--SiEDTA NPs were rinsed thrice in 50/50 ethanol/water by centrifugation at 12000 rcf for 45 min and discarding the supernatant. After the final rinse, the RITC--SiEDTA NPs were resuspended in de-ionized water via ultrasonication. RITC--SiEDTA NPs were prepared at 0.5%, 0.25% and 0.125% w/v in 1% agarose by mixing 500 L, 250 L and 125 L of 1% w/v RITC--SiEDTA NPs with 500 L, 750 L and 875 L of a warm 2% agarose solution, respectively. 100 L of the RITC--SiEDTA NPs in 1% agarose were deposited on glass slides and the agarose was allowed to spread and set at room temperature. Additionally, 37 L of 10% w/v PS beads were added to each mixture. 100 L of the mixture of PS beads and RITC--SiEDTA NPs in 1% agarose were then spread on glass slides and the agarose was allowed to set at room temperature. The mixture of fluorescent silica particles, agarose and polystyrene beads mimicked the optical scattering inside a brain tissue .
5.3 Simulation Parameters
The simulation parameters in this research are extracted from previous works [24, 23]. Unless otherwise noted, all parameters are set as follows: , ns, s, , , , , . For the simulations with 1PEF, we use , nm, , . For the simulations with 2PEF, we use , nm, , . In Figs. 3 and 4, the excitation irradiances of each step for 1PEF are , , , , , , and the excitation irradiances for 2PEF are , , , , , , where the -step SOS uses irradiances from to .
5.4 Optical System
Both confocal and MPM imaging were performed on a Nikon A1R-MP confocal system. All samples were imaged using a 100x, 1.45 NA, oil-immersion objective (Nikon Plan Apo) except for the NPs in brain phantom sample, which was imaged using a 40x, 1.15 NA, water-immersion objective (Nikon Apo LWD). The excitation in the confocal imaging was generated by a LU4/LU4A laser unit. The samples with 1PEF were excited at a wavelength of 488 nm (biological test slide and fixed cells sample) and 561 nm (NPs in gel sample). The pinhole size in all confocal imaging was set to 1.2 AU. The excitation in the MPM was generated by a Spectra-Physics Mai Tai DeepSee femtosecond laser. All samples with 2PEF were excited at a wavelength of 800 nm. The laser power levels for both confocal and MPM systems were measured at the sample plane with a digital optical power meter (Thorlabs PM100D) and a S120C detector (Thorlabs S120C). A careful power calibration for both confocal and MPM lasers was performed to correspond laser power settings (%) with actual focal power measurements (mW). This calibration was required in order to correctly apply the two-step SOS method because the linear combination coefficients in SOS microscopy are based on excitation powers.
5.5 Image Acquisition and Processing
All images were 512x512 pixels and were line-averaged for 16 times for the best SNR. As two-step SOS microscopy requires a linear combination of two images, for both confocal and MPM imaging, each sample was imaged twice at two different laser powers to generate two raw images. For the NPs in gel sample (Fig. 6), the confocal images were obtained with laser powers of 1.0% (3.7 W) and 1.3% (4.6 W), the PMT HV gain of 110, and the pixel time of 12.1 s; the MPM imaging was performed with laser powers of 1.0% (3.45 mW) and 1.2% (3.6 mW), the PMT HV gain of 140, and the pixel time of 12.1 s. For the biological test slide (Fig. 7), the confocal images were obtained with laser powers of 1.0% (2.8 W) and 1.2% (3.2 W), the PMT HV gain of 125, and the pixel time of 4.8 s; the MPM imaging was performed with laser powers of 1.0% (3.45 mW) and 1.2% (3.6 mW), the PMT HV gain of 125, and the pixel time of 12.1 s. For the fixed cells sample (Fig. 8), the confocal images were obtained with laser powers of 1.0% (2.8 W) and 1.2% (3.2 W), the PMT HV gain of 125, and the pixel time of 12.1 s; the MPM imaging was performed with laser powers of 1.0% (3.45 mW) and 1.2% (3.6 mW), the PMT HV gain of 140, and the pixel time of 12.1 s. For the NPs in brain phantom sample (Fig. 9), the MPM imaging at various depth was performed, with laser powers of 1.0% (3.45 mW) and 1.2% (3.6 mW), the PMT HV gain of 130, and the pixel time of 12.1 s.
The resolution improvement in SOS microscopy is observed from the narrowing of its PSF. The improvement is fundamentally resulted from the broadening, or the increased FWHM, of its optical transfer function (OTF), which is the Fourier transform of the PSF and the indicator of how much spatial frequency content the image can support. To demonstrate the broadening of spatial frequency content, we plot PSFs and corresponding OTFs of two steps of images obtained at different excitation intensities and their processed two-step SOS image in Fig. 10. From the first step image to the second, the excitation intensity is increased; therefore the PSF is broadened due to the saturation behavior. Nevertheless, the processed two-step SOS image has a PSF narrower than both of them, as shown in Fig. 10(a). Correspondingly, the OTF of the first step image is a Gaussian function with a limited frequency bandwidth, while the OTF of the second step brings about an evident side lobe around the main lobe, as shown in Fig. 10(b). In combination with the main lobe, the second step’s side lobe contributes to an “increased” frequency bandwidth compared to the first step, although the “increase” of bandwidth here is not regular as there are zeros in the spatial frequency band. With the linear combination algorithm in two-step SOS microscopy, the OTF of the processed two-step SOS image is a single-lobed Gaussian function, which has the largest FWHM and the highest spatial resolution. Note that the zeros in the second step image’s OTF is eliminated in the two-step SOS image due to the linear combination such that the “increase” of spatial frequency bandwidth here in SOS microscopy is in the regular sense. This strategy is different than deconvolution because deconvolution cannot increase the spatial frequency bandwidth of an image.
Compared with established super-resolution methods including STED/RESOLFT, PALM, STORM, and SIM, the implementation of SOS microscopy is very easy, fast and requires no new instrumentation. Only two raw images are required to generate a super-resolution SOS image with a -fold increase in spatial resolution, and theoretically, the resolution improvement can be infinite. However, the super-resolution capability of SOS microscopy is limited by a sacrifice in SNR. In practice, an experimental SOS image generated from more than two steps of raw images is generally unacceptable for its SNR is at least two orders of magnitude lower than a conventional image. Therefore, it is very important to accommodate both resolution and SNR in SOS imaging. Here we show that the SNR of a SOS microscopy can be calculated analytically, which can be used as a guideline for researchers to optimize the resolution and SNR performance in SOS imaging. For an -step SOS, , we consider a single pixel at and denote the experimentally obtained values of and as random variables and , respectively. The mean of is , which can be calculated using Eq. 3 if the excitation intensities and optical parameters are known. Due to the high detection efficiency of photomultiplier tubes (PMTs) in our confocal and two-photon microscopes, it is valid to consider shot noise as the only noise source in modeling the system. Therefore, the variance of , which equals to , is . Since are mutually independent, the mean and variance of an SOS image are , , respectively. Consequently, the SNR of measuring is .
Although SOS shares a similar physical principle with SAX microscopy, the implementation of SOS is not only simpler but also superior for its ”weak saturation” and versatility. SOS works with ”weak saturation” where the excitation power is relatively low (e.g., mW in 2PEF SOS) compared with the reported excitation intensities used in SAX microscopy (e.g., mW in 2PEF SAX) ; it permits super-resolution imaging deep in scattering samples within the limit of multiphoton depth when combined with multiphoton microscopy; it also enables super-resolution 3D imaging when combined with depth resolved modalities. In principle, as long as a conventional fluorescence image with an acceptable SNR can be obtained, regardless of the fluorophores used, the excitation scheme utilized, and the image depth, a corresponding super-resolution SOS image can always be generated.
Potential limitations of SOS microscopy include photobleaching and image mismatch. Photobleaching is usually accompanied by high excitation intensity and illumination dose . Although SOS microscopy works in the “weak saturation” region with low excitation intensities, the prolonged integration time aiming for better SNR performance may cause photobleaching, which degrades imaging quality and resolution. This problem can be mitigated using anti-bleaching agents , or special fluorophores like quantum dots  and nanodiamonds . Image mismatch is an issue arising in methods that require subtraction of two images . Since SOS microscopy requires the linear combination of several raw images of the same field of view, image mismatch could happen due to sample drift, which can be eliminated using image registration algorithms .
In this paper, we have proposed and demonstrated a new kind of super-resolution fluorescence microscopy using the principle of SOS. SOS microscopy employs the “weak saturation” phenomenon, which does not require high excitation intensities and therefore the damage to sample can be minimized. With a set of properly chosen linear combination coefficients, steps of raw fluorescence images can be linearly combined to generate an -step SOS image, which provides a -fold increase in spatial resolution compared with conventional diffraction-limited images. For example, linearly combining two fluorescence images obtained at regular powers extends resolution by a factor of beyond the diffraction limit. The resolution improvement of an -step SOS microscopy has been validated with simulations. With the imaging experiments performed on various samples including fixed cells, we have demonstrated the super-resolution capability of the two-step SOS microscopy with both confocal (1PEF) and multiphoton (2PEF) modalities. Owing to the deep penetration of 2PEF, the multiphoton two-step SOS has been verified to provide super-resolution imaging deep in scattering samples. The implementation of SOS microscopy is straightforward and requires neither additional hardware nor complex post-processing that are usually required in super-resolution microscopes. A conventional fluorescence microscope can be easily converted into an SOS microscope with super-resolution imaging capability.
National Science Foundation (NSF) (CBET-1554516, DMR-1309587); Defense Advanced Research Projects Agency (DARPA-14-56-A2P-PA-055); American Heart Association (16SDG31230034); Indiana CTSI (NIH/NCRR UL1TR001108).
Zhang’s research was supported by the Berry Family Foundation Graduate Fellowship of Advanced Diagnostics & Therapeutics (AD&T), University of Notre Dame. Roeder and Nallathamby’s research was supported by a grant from the Walther Cancer Foundation. The authors further acknowledge the Notre Dame Integrated Imaging Facility (NDIIF) for the use of the A1R-MP confocal microscopy in NDIIF’s Optical Microscopy Core, the Notre Dame Center for Nanoscience and Nanotechnology (NDnano), and the Notre Dame Center for Environmental Science and Technology (CEST).
-  S. W. Hell and J. Wichmann, “Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy,” Opt. Lett. 19, 780–782 (1994).
-  T. A. Klar, S. Jakobs, M. Dyba, A. Egner, and S. W. Hell, “Fluorescence microscopy with diffraction resolution barrier broken by stimulated emission,” Proc. Natl. Acad. Sci. USA 97, 8206–8210 (2000).
-  M. Hofmann, C. Eggeling, S. Jakobs, and S. W. Hell, “Breaking the diffraction barrier in fluorescence microscopy at low light intensities by using reversibly photoswitchable proteins,” Proc. Natl. Acad. Sci. USA 102, 17,565–17,569 (2005).
-  E. Betzig, G. H. Patterson, R. Sougrat, O. W. Lindwasser, S. Olenych, J. S. Bonifacino, M. W. Davidson, J. Lippincott-Schwartz, and H. F. Hess, ‘‘Imaging intracellular fluorescent proteins at nanometer resolution,” Science 313, 1642–1645 (2006).
-  M. J. Rust, M. Bates, and X. Zhuang, “Sub-diffraction-limit imaging by stochastic optical reconstruction microscopy (STORM),” Nat. Methods 3, 793–796 (2006).
-  M. G. L. Gustafsson, “Surpassing the lateral resolution limit by a factor of two using structured illumination microscopy,” J. Microsc. 198, 82–87 (2000).
-  M. G. L. Gustafsson, “Nonlinear structured-illumination microscopy: Wide-field fluorescence imaging with theoretically unlimited resolution,” Proc. Natl. Acad. Sci. USA 102, 13,081–13,086 (2005).
-  S. W. Hell, “Far-field optical nanoscopy,” Science 316, 1153–1158 (2007).
-  F. Balzarotti, Y. Eilers, K. C. Gwosch, A. H. Gynnå, V. Westphal, F. D. Stefani, J. Elf, and S. W. Hell, “Nanometer resolution imaging and tracking of fluorescent molecules with minimal photon fluxes,” Science 355, 606–612 (2017).
-  S. Gigan, “Optical microscopy aims deep,” Nat. Photonics 11, 14–16 (2017).
-  P. W. Winter, A. G. York, D. D. Nogare, M. Ingaramo, R. Christensen, A. Chitnis, G. H. Patterson, and H. Shroff, “Two-photon instant structured illumination microscopy improves the depth penetration of super-resolution imaging in thick scattering samples,” Optica 1, 181–191 (2014).
-  A. Egner, S. Jakobs, and S. W. Hell, “Fast 100-nm resolution three-dimensional microscope reveals structural plasticity of mitochondria in live yeast,” Proc. Natl. Acad. Sci. USA 99, 3370–3375 (2002).
-  A. D. Nguyen, F. Duport, A. Bouwens, F. Vanholsbeeck, D. Egrise, G. Van Simaeys, P. Emplit, S. Goldman, and S.-P. Gorza, “3D super-resolved in vitro multiphoton microscopy by saturation of excitation,” Opt. Express 23, 22,667–22,675 (2015).
-  J. J. Field, K. A. Wernsing, S. R. Domingue, A. M. Allende Motz, K. F. DeLuca, D. H. Levi, J. G. DeLuca, M. D. Young, J. A. Squier, and R. A. Bartels, “Superresolved multiphoton microscopy with spatial frequency-modulated imaging,” Proc. Natl. Acad. Sci. USA 113, 6605–6610 (2016).
-  W. Denk, J. H. Strickler, and W. W. Webb, “Two-photon laser scanning fluorescence microscopy,” Science 248, 73–76 (1990).
-  W. R. Zipfel, R. M. Williams, and W. W. Webb, “Nonlinear magic: multiphoton microscopy in the biosciences,” Nat. Biotechnol. 21, 1369–1377 (2003).
-  E. E. Hoover and J. A. Squier, “Advances in multiphoton microscopy technology,” Nat. Photonics 7, 93–101 (2013).
-  K. Fujita, M. Kobayashi, S. Kawano, M. Yamanaka, and S. Kawata, “High-resolution confocal microscopy by saturated excitation of fluorescence,” Phys. Rev. Lett. 99, 228,105 (2007).
-  J. Humpolíčková, A. Benda, and J. Enderlein, “Optical saturation as a versatile tool to enhance resolution in confocal microscopy,” Biophys. J. 97, 2623–2629 (2009).
-  G. Vigil, Y. Zhang, A. Khan, and S. Howard, “Description of deep saturated excitation multiphoton microscopy for super-resolution imaging,” J. Opt. Soc. Am. A 34, 1217–1223 (2017).
-  R. Oketani, A. Doi, N. I. Smith, Y. Nawa, S. Kawata, and K. Fujita, “Saturated two-photon excitation fluorescence microscopy with core-ring illumination,” Opt. Lett. 42, 571 (2017).
-  C. Kuang, S. Li, W. Liu, X. Hao, Z. Gu, Y. Wang, J. Ge, H. Li, and X. Liu, “Breaking the diffraction barrier using fluorescence emission difference microscopy,” Sci. Rep. 3, 1441 (2013).
-  Y. Zhang, G. D. Vigil, L. Cao, A. A. Khan, D. Benirschke, T. Ahmed, P. Fay, and S. S. Howard, “Saturation-compensated measurements for fluorescence lifetime imaging microscopy,” Opt. Lett. 42, 155–158 (2017).
-  C. Eggeling, A. Volkmer, and C. A. M. Seidel, “Molecular photobleaching kinetics of rhodamine 6G by one- and two-photon induced confocal fluorescence microscopy,” ChemPhysChem 6, 791–804 (2005).
-  C. Xu and W. W. Webb, “Measurement of two-photon excitation cross sections of molecular fluorophores with data from 690 to 1050 nm,” J. Opt. Soc. Am. B 13, 481–491 (1996).
-  E. Gatzogiannis, X. Zhu, Y.-T. Kao, and W. Min, “Observation of frequency-domain fluorescence anomalous phase advance due to dark-state hysteresis,” J. Phys. Chem. Lett. 2, 461–466 (2011).
-  Y. Yonemaru, M. Yamanaka, N. I. Smith, S. Kawata, and K. Fujita, “Saturated excitation microscopy with optimized excitation modulation,” ChemPhysChem 15, 743–749 (2014).
-  Y. Zhang, A. A. Khan, G. D. Vigil, and S. S. Howard, “Investigation of signal-to-noise ratio in frequency-domain multiphoton fluorescence lifetime imaging microscopy,” J. Opt. Soc. Am. A 33, B1–B11 (2016).
-  R. M. Zucker, “Quality assessment of confocal microscopy slide based systems: Performance,” Cytometry Part A 69A, 659–676 (2006).
-  K. Yoshida, I. Nishidate, T. Ishizuka, S. Kawauchi, S. Sato, and M. Sato, “Multispectral imaging of absorption and scattering properties of in vivo exposed rat brain using a digital red-green-blue camera,” J. Biomed. Opt. 20, 051,026 (2015).
-  P. D. Nallathamby, J. Hopf, L. E. Irimata, T. L. McGinnity, and R. K. Roeder, “Preparation of fluorescent Au-SiO core-shell nanoparticles and nanorods with tunable silica shell thickness and surface modification for immunotargeting,” J. Mater. Chem. B 4, 5418–5428 (2016).
-  P. D. Nallathamby, N. P. Mortensen, H. A. Palko, M. Malfatti, C. Smith, J. Sonnett, M. J. Doktycz, B. Gu, R. K. Roeder, W. Wang, and S. T. Retterer, “New surface radiolabeling schemes of super paramagnetic iron oxide nanoparticles (SPIONs) for biodistribution studies,” Nanoscale 7, 6545–6555 (2015).
-  W. Wang, P. D. Nallathamby, C. M. Foster, J. L. Morrell-Falvey, N. P. Mortensen, M. J. Doktycz, B. Gu, and S. T. Retterer, ‘‘Volume labeling with Alexa Fluor dyes and surface functionalization of highly sensitive fluorescent silica (SiO) nanoparticles,” Nanoscale 5, 10,369–10,375 (2013).
-  T. Bernas, J. P. Robinson, E. K. Asem, and B. Rajwa, “Loss of image quality in photobleaching during microscopic imaging of fluorescent probes bound to chromatin,” J. Biomed. Opt. 10, 064,015 (2005).
-  J. Vogelsang, R. Kasper, C. Steinhauer, B. Person, M. Heilemann, M. Sauer, and P. Tinnefeld, “A reducing and oxidizing system minimizes photobleaching and blinking of fluorescent dyes,” Angew. Chem. Int. Edit. 47, 5465–5469 (2008).
-  D. R. Larson, “Water-soluble quantum dots for multiphoton fluorescence imaging in vivo,” Science 300, 1434–1436 (2003).
-  E. Rittweger, K. Y. Han, S. E. Irvine, C. Eggeling, and S. W. Hell, “STED microscopy reveals crystal colour centres with nanometric resolution,” Nat. Photonics 3, 144–147 (2009).
-  N. S. Alexander, G. Palczewska, P. Stremplewski, M. Wojtkowski, T. S. Kern, and K. Palczewski, “Image registration and averaging of low laser power two-photon fluorescence images of mouse retina,” Biomed. Opt. Express 7, 2671–2691 (2016).