Quantitative mapping and minimization of super-resolution optical imaging artifacts (original) (raw)

Super-resolution image simulation with SuReSim.

In order to simulate disappearance of a filament from a realistic microtubule network, a real super-resolution image of microtubules (Fig. 1c) was used as a support for SuReSim data simulation. Raw data of blinking Alexa 647–labeled microtubules imaged using TIRF were reconstructed using ThunderSTORM maximum likelihood multiemitter fitting and then loaded into the SuReSim software and all filaments were traced using the editor function and the WIMP file saved. SuReSim was used to generate a simulated super-resolution reconstruction of all filaments, which was then convolved by a Gaussian PSF to generate a simulated reference image. The object in the WIMP file corresponding to the filament highlighted in Figure 1d,e was deleted, and SuReSim was used again to render a simulated super-resolution reconstruction, except this time missing a filament. SuReSim was also used for the simulations in Supplementary Note 1.

Cell lines and primary cells.

HeLa cells (Figs. 1 and 2) were kindly provided by M. Marsh, MRC LMCB, University College London, and cultured in phenol-red free DMEM (Gibco) supplemented with 2 mM GlutaMAX (Gibco), 50 U/ml penicillin, 50 μg/ml streptomycin and 10% FBS (FBS; Gibco). CHO cells (Fig. 3) were cultured in phenol red-free Minimum Essential Medium Alpha (MEMα; Gibco) supplemented with 10% FBS (Gibco) and 1% penicillin/streptomycin (Sigma).

Rat hippocampal neurons and glial cells (Supplementary Notes 10 and 11) were harvested from embryonic day 18 pups, following established guidelines of the European Animal Care and Use Committee (86/609/CEE) and approval of the local ethics committee (agreement D13-055-8), and cultured in Neurobasal medium (Gibco) supplemented with 2 mM GlutaMAX-I (Gibco) and B27 supplement (Gibco). All cells were grown at 37 °C in a 5% CO2 humidified incubator.

Sample preparation for fixed microtubules.

For TIRF-SMLM imaging of microtubules (Figs. 1 and 2), 13 mm diameter, thickness #1.5 coverslips were washed overnight in 1:1 HCl:methanol and washed 5 times in ddH2O and twice in 90% isopropyl alcohol. Coverslips were then incubated overnight in poly-L-lysine (0.01%) (Sigma Aldrich) and rinsed twice in PBS. HeLa cells were seeded on these coverslips and grown overnight in 12-well plates. Cells were fixed with 4% PFA in cytoskeleton buffer (10 mM MES, pH 6.1, 150 mM NaCl, 5 mM EGTA, 5 mM glucose, 5 mM MgCl2) for 15 min at 37 °C, washed 3× with PBS, then permeabilized with 0.1% Triton X-100 in PBS for 10 min and blocked in 2.5% BSA in PBS for a further 30 min. Samples were then labeled with 2 μg/ml anti-α-tubulin (DM1A mouse monoclonal, Sigma Aldrich) in 2.5% BSA in PBS for 1 h, followed by 3× washes with PBS and labeling with Alexa Fluor 647–labeled goat anti-mouse secondary antibody (Invitrogen) (2 μg/ml in 2.5% BSA in PBS) for 1 h. Samples were washed 3× with PBS and fixed again in 4% PFA in cytoskeleton buffer for 10 min, before being washed 3× with PBS. Samples were mounted on a Parafilm-formed gasket3 in STORM buffer (150 mM Tris, pH 8, 10 mM NaCl, 1% glycerol, 1% glucose, 1% BME), sealed with clear nail varnish (Maybelline) and imaged within 3 h of mounting.

For imaging in different focal volumes (Supplementary Note 5), COS cells were fixed with glutaraldehyde and labeled with two monoclonal mouse anti-α-tubulin antibodies (DM1A and B-5-1-2, both from Sigma) and a goat anti-mouse Alexa Fluor 647–labeled secondary antibody (Thermo Fisher Scientific). Samples were mounted in Smart Buffer (Abbelight) for imaging.

For widefield super-resolution imaging of microtubules (Fig. 3), CHO cells were seeded on ultra-clean3 8 mm diameter thickness #1.5 coverslips (Zeiss) at a density of 1 × 105 per 35 mm dish. Fixation was performed with 4% PFA in a modified version of cytoskeleton stabilizing buffer (CSB) (5 mM KCl, 0.1 mM NaCl, 4 mM NaHCO3, 11 mM Na2HPO4, 2 mM MgCl2, 5 mM PIPES, 2 mM EGTA, pH 6.9) for 15 min at 37 °C, followed by washing with the same CSB (without PFA). Additional permeabilization was performed (0.05% Triton X-100 in CSB) for 5 min followed by three washing steps using 0.05% Tween-20 in the modified version of CSB and blocking in 5% BSA (Sigma) for 40 min. Microtubules were stained and submitted to a secondary fixation step as described above. 100 nm TetraSpeck beads (Life Technologies) were added at a dilution of 1:1000 in PBS for 10 min to each coverslip. Coverslips were mounted on clean microscope slides3 in 100 mM mercaptoethylamine (Sigma) at pH 7.3 and all imaging was performed within 3 h of mounting.

Fixed microtubule imaging.

Fixed microtubule samples were imaged by TIRF-SMLM (Figs. 1 and 2) on a N-STORM microscope (Nikon Instruments), using a 100× TIRF objective (Plan-APOCHROMAT 100×/1.49 Oil, Nikon) with additional 1.5× magnification. A reference TIRF image was acquired with 5% power 647 nm laser illumination and 100 ms exposure time, before SMLM data acquisition of 40 000 frames at 100% power 647 nm illumination with 405 nm stimulation and an exposure time of 30 ms per frame.

Imaging with different illumination regimes (Supplementary Note 5) was performed on an N-STORM microscope using a 100×, 1.49 NA objective as above, but with no additional magnification and an exposure time of 15 ms. Prior to dSTORM imaging a reference image was acquired using a high-pressure mercury lamp (Intensilight, Nikon) with a Cy5 filter cube (Nikon); the filter cube was then switched and the laser illumination set to either vertical (i.e. widefield), HILO, or TIRF. A second reference image was then acquired, this time with laser illumination. A cylindrical lens was inserted into the detection path and 60,000 frame dSTORM data set acquired at this angle.

Widefield and super-resolution imaging for fusion (Fig. 3) was carried out on a Zeiss Elyra PS.1 inverted microscope system, using a 100× TIRF objective (PlanAPOCHROMAT 100×/1.46 Oil, Zeiss) and additional 1.6× magnification. The sample was illuminated with a 642 nm laser operating at 100% laser power. 45000 frames were acquired with a 20 ms exposure time per frame.

Reconstruction algorithms for dSTORM data.

The freely available software packages ThunderSTORM17 (Figs. 1,2,3, Supplementary Notes 5, 7, and 10), SRRF (Fig. 3; Supplementary Notes 7 and 10) and QuickPALM18 (Fig. 3 and Supplementary Note 10) were used for super-resolution image reconstruction. Images labeled 'MLE' were reconstructed with ThunderSTORM with the integrated PSF method with maximum likelihood fitting and multiemitter fitting enabled. Drift correction was performed post-localization and images were rendered using a normalized 20 nm Gaussian. Images labeled 'SRRF' were analyzed with the most appropriate parameters for each individual data set and drift corrected during analysis. Images labeled 'CoM' were reconstructed using QuickPALM with the default parameters, following drift correction of the raw data using the NanoJ-SRRF package. The particle tables from QuickPALM were then loaded into ThunderSTORM for rendering using a normalized 20 nm Gaussian.

SIM imaging.

For SQUIRREL analysis of SIM images (Supplementary Note 4), FluoCells prepared slide #2 (Invitrogen) with BPAE cells stained with Texas Red-X phalloidin and Alexa Fluor 488-tubulin was imaged on a Zeiss Elyra PS.1 system, using a 63× NA 1.4 objective with additional 1.6× magnification for SIM and widefield acquisition. For actin imaging, 'Low SNR' images were acquired with a 561 nm laser at 0.05% laser power, using 100 ms exposure time, and 5 grid rotations. 'High SNR' images were acquired with a 561 nm laser at 5% laser power, 100 ms exposure, 5 grid rotations. Widefield images were acquired with a 561 nm laser at 0.2% laser power, 100 ms exposure time. SIM reconstructions were generated with the Zeiss Elyra Zen software using automatic settings. For microtubule imaging, raw SIM data were acquired with a 488 nm laser at 10% laser power using 100 ms exposure time and three grid rotations. The SIM reconstruction was generated using FairSIM19.

Generation and analysis of synthetic data at different z-positions.

The bead images used in Supplementary Note 5 were obtained from the open source data set 'z-stack-Bead-2D-Exp-as-stack' available to download from the SMLMS Challenge 2016 website,

http://bigwww.epfl.ch/smlm/challenge2016/datasets/Beads/Data/data.html (data used here were downloaded on 4th September 2017). This data set comprises 151 slices of an image of six fluorescent beads covering the z-range −750 nm to 750 nm (step size = 10 nm). The central x,y location of each of the six beads in this was determined at the central plane of the z-stack, and this was used to define the center of a 3.3 μm × 3.3 μm region about each bead. For generation of the data set containing PSFs from all z-positions in each frame, these regions were pasted into an image where the x,y coordinates mapped to a specific z-position from the bead image stack. The target x,y coordinates for pasting the images were spaced such that there was 5 μm between adjacent bead centers, and regions were randomly from the six original bead images. This was repeated 1000 times to generate a 1000-frame data set. Gaussian-Poisson noise was then added to the image stack to mask the edges of the pasted bead images. This data set was then analyzed with SRRF and ThunderSTORM (default software settings in both cases) to produce a single super-resolution image for each algorithm. The reference image was an average projection of all 1000 frames. For generation of the data set containing constant z-positions in each frame, regions from the bead z-stack were again selected and pasted but this time z was varied between slices as opposed to the x,y position within each frame. 10 frames were produced for each z-position, and noise was added again as above. ThunderSTORM and SRRF analyses (default settings) were then run on this image stack 10 slices at a time to generate a single super-resolution image for each z-position. The reference was the average of the 10 frames corresponding to z = 0 nm.

For assessing the impact of out-of-focus fluorescence on defect detection (Supplementary Note 5), a test structure was simulated consisting of three semicircles of radius 500 nm and axial tilt ranging from −50 nm to +750 nm. A widefield representation of this structures was produced via convolution with a 3D PSF generated using the ImageJ PSF Generator plugin20 with the following settings: Born and Wolf optical model, numerical aperture 1.4, wavelength 640 nm, z-range 1500 nm, z-step size 2 nm. Single molecule blinking data sets were generated with custom-written simulation software with the same PSF used for rendering molecule appearances, and were binned into 100 nm 'camera' pixels with Gaussian-Poisson noise. This was performed for both the defect-free structure and an artifactual equivalent where 100 nm stretches of the structure had been deleted. These data sets were analyzed using weighted-least-squares fitting in ThunderSTORM.

VACV sample preparation and imaging.

2.5 × 106 VACV particles (WR strain, EGFP-F18 in tk locus21) were diluted in 100 μl 1 mM TRIS, pH 8, sonicated for 3 × 30 s and incubated on gridded #1.5 glass-bottom petri dishes (Zell-Kontakt GmbH) for 1 h at room temperature and fixed with 4% PFA in PBS for 15 min. Samples were quenched with 50 mM NH4Cl in PBS for 10 min, washed in PBS, and incubated in permeabilization/blocking buffer (1% Triton X-100, 5% BSA, 1% FBS in PBS) for 30 min. Samples were labeled in blocking/permeabilization buffer overnight at 4 °C or 2 h at room temperature with anti-GFP nanobodies (Chromotek), labeled in-house with Alexa Fluor 647 NHS-ester (Life Technologies) with a dye-to- protein ratio of approximately 1, as previously described22. Samples were then washed 3× with PBS, fixed in 4% PFA in PBS for 10 min, quenched with 50 mM NH4Cl in PBS for 10 min and washed in PBS.

VACV samples were imaged in STORM buffer on a Zeiss Elyra PS.1 system, using a 100× TIRF objective with additional 1.6× magnification (as above) for SIM, SRRF and SMLM acquisition (Supplementary Note 7). Buffer was exchanged to PBS and STED images were acquired on a Leica SP8, re-localizing the same region of interest based on the grid. SMLM data acquisition parameters were 30,000 frames at 100% laser power 647 nm illumination with 405 nm stimulation and an exposure time of 33 ms per frame.

Clathrin coated pits sample preparation and imaging.

Rat glial cells (from embryonic day 18 pups) were cultured on 18-mm coverslips at a density of 4000 /cm2, respectively. After 9 days in culture, samples were fixed using 4% PFA in PEM (80 mM PIPES, 2 mM MgCl2, 5 mM EGTA, pH 6.8) for 10 min. For PAINT imaging23 of clathrin coated pits (CCPs) in glial cells, fixed samples were incubated with a rabbit anti-clathrin primary antibody (abCam, catalog #21679) overnight at 4 °C, then with an anti-rabbit DNA–conjugated secondary antibody (Ultivue) for 1 h at room temperature.

DNA-PAINT imaging of CCPs in glial cells (Supplementary Note 10) was performed on a N-STORM microscope using a 100× objective as above. The same glial cell (present in low numbers in hippocampal cultures) was imaged in serial dilutions of Imager-650 (2 mM stock, from lower to higher concentration) in Imaging Buffer (Ultivue). The sample was illuminated at 647 nm (50% laser power) and a sequence of 20,000 images at 33 Hz was acquired for each Imager-650 dilution, before switching to a higher concentration of Imager-650 in Imaging Buffer.

Actin rings sample preparation and imaging.

Rat hippocampal neurons (from embryonic day 18 pups) were cultured on 18-mm coverslips at a density of 10,000 After 9 days in culture, samples were fixed using 4% PFA in PEM (80 mM PIPES, 2 mM MgCl2, 5 mM EGTA, pH 6.8) for 10 min. Preparation of actin-stained neurons for SMLM was performed similarly to the protocol described in24, with minor modifications. After blocking, neurons were incubated with a mouse anti-map2 primary antibody (Sigma Aldrich, catalog #M4403) for 1h 30 min at room temperature (RT), then with a Alexa Fluor 488–labeled donkey anti-mouse secondary antibody (Thermo Fisher) for 45 min at RT, then with 0.5 mM phalloidin–Alexa Fluor 647 (Thermo Fisher) overnight at 4 °C. Neurons were mounted in a modified STORM buffer (50 mM Tris, pH 8, 10 mM NaCl, 10% glucose, 100 mM mercaptoethylamine, 3.5 U/ml pyranose oxidase, 40 μg/mL catalase) complemented with 0.05 mM phalloidin–Alexa Fluor 647, to mitigate phalloidin unbinding during acquisition and imaged immediately.

Neuron samples were imaged on a N-STORM microscope using a 100× objective as above (Supplementary Note 11). The sample was illuminated at 100% laser power at 647 nm. A sequence of 60,000 images at 67 Hz was acquired. Images were rendered with ThunderSTORM using a normalized 20 nm Gaussian from particle tables generated with SMAP, a MATLAB based software package developed by the Ries group at the EMBL, Heidelberg. Localizations were determined using a probability based method after background subtraction by wavelet filtering and lateral drift was corrected by cross-correlation.

Visibility analysis.

To quantify the quality of the super-resolution reconstructions of parallel actin rings (Supplementary Note 11), a normalized visibility similar to that described in Geissbuehler et al.25 was calculated as follows. Average intensity profiles were plotted for a 0.5 × 1 μm stretch of axon containing 5 actin rings for each of the 120 reconstructed images. The MATLAB function findpeaks was used to find the 5 peak positions in the average profile measured from the 60,000 frames reconstruction, and mean pairwise visibility was calculated as follows.

_I_max,i and _I_max,i+1 are the intensities at peak positions i and i + 1, respectively, where i denotes the index of the actin ring in the sampled regions and _I_min,i_→_i + 1 is the intensity at the midpoint of two adjacent peaks. Higher visibilities correspond to a greater ability to differentiate between two structures up to a maximum value of .

Color maps.

Color maps used for displaying images ('NanoJ-Orange'), error maps ('SQUIRREL-errors') and FRC maps ('SQUIRREL-FRC') are provided in the NanoJ-SQUIRREL software package.

Software availability.

The version of NanoJ-SQUIRREL used in this paper is available as Supplementary Software. The most recent version of NanoJ-SQUIRREL can be downloaded and installed in ImageJ and Fiji automatically by following the instructions in the manual, available here: https://bitbucket.org/rhenriqueslab/nanoj-squirrel. Source code is also available at the same website.

Data availability.

The data that support the findings of this study are available from the corresponding authors upon reasonable request. Sample data sets can be downloaded from https://bitbucket.org/rhenriqueslab/nanoj-squirrel.