$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
The protocol described above enables the preparation of high-quality Arabidopsis thaliana nuclei suitable for SMLM. Several rounds of optimization were performed to improve nuclear integrity and reduce cytoplasmic and cellular debris that compromise image quality. After each optimization step, the quality of the nuclear preparations was evaluated by confocal imaging using an Olympus IX81 spinning-disk microscope before proceeding to super-resolution acquisition.
Confocal imaging
The procedure described above is indeed compatible with confocal imaging. For this purpose, the LMA solution used in step 3.1 can be based on PBS (or any neutral buffer or medium). MEA and GLOX are only necessary for SMLM. On the other hand, the JF-based Hoechst used for SMLM (step 3.3) is too diluted to be detectable by conventional imaging methods. DNA labeling for confocal or epifluorescence imaging was thus achieved using Hoechst H33258 at 1/500 dilution (stock at 1 mg/mL).
Using the present optimized protocol, nuclei were isolated from 15-day-old Arabidopsis wild-type seedlings. The actively elongating form of RNA Polymerase II, marked by phosphorylation of Ser2 of the heptapeptide tandem repeats on the C-terminal domain (CTD), was immunostained with a specific primary antibody29 and an AF647-conjugated secondary antibody26. Confocal inspection illustrated in Figure 1 shows various isolated and well-preserved nuclei with minimal debris, most of which did not show Hoechst, AF647, or chlorophyll autofluorescence signals. Chromocenters (revealed by Hoechst) and the nucleolus were clearly distinguishable. Notably, RNA Polymerase II signal was absent from both chromocenters and the nucleolus (arrows), consistent with the expected localization of the enzyme in its active, transcriptionally engaged form.

Figure 1: Examples of three different ROIs with extracted nuclei (rows), displaying the active RNA Polymerase II immunolabeled with AF647-conjugated antibody (red), and DNA counterstained with Hoechst H33258 (cyan). Confocal imaging was used to assess sample quality prior to Single-Molecule Localization Microscopy (SMLM). The preparation contains multiple intact, well-preserved nuclei clearly separated from debris, with chromocenters and nucleoles readily visible (arrows), indicating suitability for subsequent SMLM. DIC: Differential Interference Contrast. Scale bar: 3 μm (all panels). Please click here to view a larger version of this figure.
STORM imaging
For super-resolution imaging, samples were stained for RNA Polymerase II and counterstained for DNA with JF549-Hoechst. Among a few DNA-binding SMLM dyes (e.g., SiR-DNA, JF646-Hoechst), JF549-Hoechst was chosen as a second label in combination with AF647 for two-color imaging. Ten to thirty thousand frames were sequentially acquired for AF647 (RNA Pol II) and JF549 (DNA). Representative results are shown in Figure 2. To simulate diffraction-limited images, the maximum intensity projections (MIP) of all detected blinking events in the raw acquisitions were generated using ImageJ2/Fiji30,31 (Figure 2, Diffraction-limited MIP). To achieve sub-diffraction-limited localizations of single-molecule signals, the raw image stacks were processed with the ThunderSTORM plugin in Fiji32 (Figure 2, SMLM). The typical analysis workflow is detailed below. The corresponding images were generated using ThunderSTORM average shifted histogram representation.
SMLM data processing by ThunderSTORM was obtained using Gaussian fitting with Weighted Least-Squares (WLSQ) estimation (Figure 3A, top). To further improve the quality and resolution of the reconstructed images, additional filtering over sigma and uncertainty was applied. An appropriate sigma value corresponds to the expected width of a true single-molecule emission, whereas excessively small or large values typically represent noise and overlapping or slightly defocused emitters, leading to reconstruction artifacts, and therefore were filtered out (Figure 3A and histograms in Figure 3B). On the other hand, localization uncertainty values can be used to select the best localizations with respect to the calculated fit accuracy33. Uncertainty reflects how precisely ThunderSTORM can localize a single molecule based on photon statistics and background noise. Localizations were therefore filtered on their uncertainty for values below 30 nm (Figure 3A and histograms in Figure 3C; Supplementary File 1). Following this second filtering step, detections with poor positional precision were removed, resulting in sharper molecular signals and reduced background noise.
Sample drift correction was carried out using ThunderSTORM cross-correlation or the COMET (Cost-function Optimized Maximal overlap drift EsTimation) algorithm (https://github.com/gpufit/comet). Note that in Figure 3A, the drift was minimal and does not appear on the non-corrected images.
After reconstruction, the nucleolus and the overall distribution of DNA within the nucleus became clearly resolved. Active RNA Pol II formed small clusters throughout the nucleoplasm and was distinctly excluded from compact DNA regions, consistent with its transcriptional activity. JF549-Hoechst produced a strong signal under SMLM conditions; Surprisingly, reconstructed images showed reduced detection of chromocenters compared with confocal images of Hoechst H33258-labeled DNA, although occasional dense DNA foci remained distinguishable (Figure 2 or Figure 3, close-up views).
To quantify the improvement in image resolution during processing, we performed Fourier Ring Correlation (FRC) analysis. FRC provides an objective estimate of spatial resolution by comparing two independent reconstructions of the same dataset and measuring their correlation across spatial frequencies34. This method, therefore, enables direct assessment of resolution gains from the fitted localizations to the fully processed data. In the example shown in Figure 4, the resolution of the AF647 channel improved from 90.1 nm to 47.7 nm after filtering (Figure 4A–C). Similarly, the resolution of the JF549 channel increased from 56.6 nm to 27.9 nm following processing.
The representative results presented here illustrate the effectiveness of the optimized nuclear preparation protocol in preserving nuclear integrity and enabling high-quality SMLM imaging. Confocal inspection confirms that isolated nuclei are largely free of cytoplasmic debris, with chromocenters and nucleoli clearly distinguishable, demonstrating that the protocol maintains structural features critical for downstream super-resolution analysis. The subsequent SMLM reconstructions further highlight the technique’s capacity to resolve sub-diffraction molecular distributions, as exemplified by the clustered localization of active RNA Polymerase II within the nucleoplasm and its exclusion from compact DNA regions. Image processing in ThunderSTORM, including filtering based on sigma and localization uncertainty, as well as drift correction, provides a robust framework to refine single-molecule localizations and reduce background noise. Fourier Ring Correlation analysis quantitatively confirms the resolution enhancement achieved through these steps, offering a reliable metric for assessing the quality of reconstructed images. Together, these results validate the protocol and provide a clear workflow for interpreting SMLM data, allowing users to identify both global nuclear organization and fine-scale molecular patterns.

Figure 2: Two-dimensional STORM imaging of Arabidopsis thaliana Col-0 wild-type nucleus immunostained with an antibody against the active RNA Polymerase II, and counterstained for DNA (JF549-Hoechst). Diffraction-limited maximum intensity projection (MIP) of the acquired streams for AF647 (left), JF549 (middle), and merged channels (right) is shown for an isolated nucleus. Corresponding ThunderSTORM reconstructions are displayed below. A 1.2 µm × 1.2 µm region of interest was selected and enlarged to visualize finer structural details. Scale bars: 1 µm. Please click here to view a larger version of this figure.

Figure 3: SMLM localization filtering workflow. (A) Localizations obtained after applying the weighted least-squares fitting method in ThunderSTORM (top, detected localizations). A 1.2 µm × 1.2 µm region of interest is enlarged to highlight structural details. The image appears blurred and partially saturated due to the presence of poorly fitted or out-of-focus fluorophore localizations. Localizations remaining after sigma and uncertainty filtering (middle). Both AF647 (RNA Pol II) and JF549 (DNA) channels show substantially improved definition compared to the top panel. The drift was corrected using the COMET algorithm (bottom). (B) Corresponding distributions of sigma values for the raw localizations with the selected filtering ranges highlighted in red. Sigma filtering was applied by selecting the dominant peak of the distribution and removing events with excessively small or large sigma values. (C) Corresponding distributions of uncertainty values for the raw localizations with the selected filtering ranges highlighted in red. Localizations with high uncertainty were eliminated. Scale bars: 1 µm. Please click here to view a larger version of this figure.

Figure 4: Fourier Ring Correlation (FRC) analysis of the nucleus presented in Figure 3A. FRC curves estimate the effective spatial resolution for the AF647 and JF549 channels following sequential application of (A) WLSQ analysis, (B) sigma filtering, and (C) uncertainty filtering. After applying all filters, the resolution improved from 90.1 nm to 47.7 nm for the AF647 channel and from 56.6 nm to 27.9 nm for the JF549 channel. The decrease in FRC-derived resolution values indicates enhanced ability to resolve finer structural details. Please click here to view a larger version of this figure.
Supplementary File 1: ThunderSTORM fitting parameters. The raw ThunderSTORM protocol file corresponding to the fitted image includes the plugin version used in the present work, specific camera settings, and detailed settings for the WLSQ fitting method. Further post-processing filtering parameters include the sigma and uncertainty ranges.Please click here to download this file.