Main Text

The ability to resolve subcellular features can greatly enhance biological discovery, identification of disease targets, and the development of targeted therapeutics. Super-resolution microscopy methods such as structured illumination microscopy (SIM), stochastic optical reconstruction microscopy (STORM) and stimulated emission depletion microscopy (STED) enable nanoscale imaging of specimens.13 However, these approaches require specialized expertise, costly reagents and microscopes. Expansion microscopy (ExM) is an inexpensive and accessible method that similarly resolves fluorescently labeled structures at subcellular resolution by isotropically enlarging specimens within a swellable polyelectrolyte hydrogel.4 Images of expanded cells show a higher effective resolution due to an increase in the distance between the fluorescent molecules. Super-resolution imaging and ExM are critical tools that are widely used in biological research; however, sample processing and imaging require time consuming manual manipulation, limiting throughput. Although there has been success in increasing the throughput of certain super-resolution microscopy techniques, these methods require highly specialized equipment and expertise which limits the accessibility of these tools.58 To circumvent these limitations, we developed High-Throughput Expansion Microscopy (HiExM) that combines parallel sample processing and automated high-content confocal imaging in a standard 96-well cell culture plate.

Adapting standard ExM to a 96-well plate requires the reproducible delivery of a small volume of gel solution to each well (<1 µL compared to ∼200 µL/sample for slide-based ExM). This small volume is necessary to allow complete gel expansion within the size constraints of wells in a 96-well plate. To achieve this goal, we engineered a device that retains small liquid droplets when dipped into a reservoir of monomeric expansion gel solution. These droplets can then be deposited into the centers of wells across the cell culture plate (Figs. 1a,b). The retention of the gel solution is mediated by a series of reticulated grooves at the end of a cylindrical post which allows for the formation of a pendant droplet (Fig. 1c). We estimate that each post retains ∼230 nL of gel solution as measured by the mass of liquid accumulated on all post tips of a single device after dipping (Supplementary Methods).

HiExM enables gel formation and expansion in a 96-well cell culture plate

a. Schematic representation of HiExM devices showing the key features highlighted in color. b. Example devices used in 96-well cell culture plates. c. Brightfield image of the conical post-tip shows the pattern of grooves that mediate fluid retention. d. Fluid retention at the conical post-tip of the device. Silhouettes taken by an optical comparator of the profile of a single post suspended above a surface (left) and in contact with a surface (right) show a fluid droplet interacting with the device. Upon device insertion, the gel solution fills the space under the conical post tip, forming the toroid gel. e-h. Schematic of HiExM gel deposition and expansion workflow. e. The device is immersed in a shallow reservoir of gel solution. f. Upon removal, the tip of each device post retains a small volume of gel solution. g. Gel solution is deposited by the device into the well centers of the cell culture plate. Brightfield image (right) shows gel geometry and size prior to expansion. Note that gels deposited in HiExM cover ∼1.1 mm2 of the well surface to accommodate the expanded gel, and do not include cells outside the gel footprint. h. Polymerization and expansion are performed with the device in place. Brightfield image (right) shows gel geometry and size after expansion.

The small volume of gel solution delivered to each well by the device forms a toroidal droplet where the inner surface is molded by the conical surface of the post-tip and the outer surface is constrained by surface tension (Figs. 1c,d). The post-tip dimensions allow delivery of a volume of gel solution that expands relative to the diameter of the well (Fig. 1d). The toroidal gel enables isotropic gel expansion within the well constraints while allowing for easy device removal because the inner surface of the gel expands outward away from the conical post-tip, effectively detaching the gel from the device. To achieve this toroidal geometry, our device design incorporates three key features to ensure each post of the device contacts the midpoint at the bottom of each well (Fig. S1). First, two sets of three cantilevers (‘pressure struts’) are located at each end of the device and flank two rows of 6 posts (12 posts per device). These pressure struts apply outward force to the inside of the outer four wells ensuring that the device remains in contact with the well surface throughout the expansion process. Second, the spine of the device incorporates several notches extending through the long side of the device allowing it to bend to accommodate plate inconsistencies between the depths of individual wells. Finally, the post lengths are offset from the center four posts to allow for sequential deformation that ensures each post contacts the well surface. Together, these design features allow for reproducible deposition and expansion of gels in a 96-well plate.

To test our device for high-throughput gel expansion, we first performed ProExM using ammonium persulfate (APS) and tetramethylethylenediamine (TEMED) to initiate gel polymerization (Fig. S2a).9 Since APS rapidly accelerates polymerization of polyacrylamide hydrogels, TEMED-containing expansion gel solution is first delivered by the device and removed leaving a ≤1µl droplet in each well of the plate on ice. Another device is then used to deliver APS-containing expansion gel solution at the same volume and the polymerization reaction is initiated by placing the plate on a 50°C surface. Given the large air-liquid-interface of the gel, droplet delivery and polymerization are performed in a nitrogen filled glove bag to minimize oxygen inhibition of the reaction. Although we observed gel formation across wells, polymerization and expansion were inconsistent, likely due to inconsistent radical formation and rapid evaporation of the small volume relative to the rate of polymerization (Movie S1). To overcome this limitation, we tested photochemical initiators including Lithium phenyl(2,4,6-trimethylbenzoyl)phosphinate (LAP), diphenyl(2,4,6-trimethylbenzoyl)phosphine oxide (TPO), and 1-[4-(2-Hydroxyethoxy)-phenyl]-2-hydroxy-2-methyl-1-propane-1-one (Irgacure 2959) which control polymerization of polyacrylamide gels when exposed to UV light.10,11 LAP has shown success in slide-based ExM using PEG-based hydrogels.12,13 We found that Irgacure 2959 allowed reproducible gel formation and expansion when compared to either APS/TEMED or other photoinitiators (Movie S2).

We next tested the ability of the HiExM platform to achieve nanoscale image resolution of biological specimens. Briefly, A549 cells were cultured, fixed, permeabilized and immunostained with alpha-tubulin antibodies to visualize microtubules as a reference structure before and after expansion. Stained samples were then incubated with Acryloyl-X (AcX) overnight at 4°C to allow anchoring of native proteins to the polymer matrix during gel formation. Device posts were immersed in the expansion gel solution containing 0.1% Irgacure 2959. Devices were inserted into the well plate followed by exposure to UV light (365nm) in an anoxic environment at room temperature to initiate the polymerization reaction. Cells embedded in the resulting gels were digested with Proteinase K then expanded in deionized water overnight (Fig. S2b, Fig. S3). Visual inspection using an epifluorescence microscope showed both robust gel formation and fluorescence signal retention across wells. Although we observed bleaching of AlexaFluor dyes in our Irgacure 2959 photopolymerized gels, Cyanine-based Fluorescent (CF) dyes are robust against bleaching in HiExM and result in reproducible signal retention (Fig. S4).14 We also titrated AcX and ProteinaseK to optimize signal retention in HiExM and found the 50 µg/mL AcX and 1U/mL ProteinaseK was most consistent for A549 cells (Fig. S5). This optimization step is likely critical for different cell types. In some cases, we observed residual Hoechst signal underneath expanded cells due to digested cells that were incompletely removed during wash steps during the expansion process (Fig. S6). This signal does not impact interpretation of results.

A major workflow bottleneck associated with super-resolution imaging, including current ExM protocols, is the significant time and expertise required for image acquisition. For example, the vertical component of expansion increases the depth of the sample, necessitating more optical z-slices to image the full depth of the sample. Moreover, expansion effectively decreases the field-of-view due to increasing sample size, necessitating acquisition of more images to analyze the same number of cells at lower resolution. To enable rapid 3D image capture of expanded cells in 96-well plates, we used a high content confocal microscope (Opera Phenix system, Perkin Elmer) (Fig. 2a). This approach allowed autonomous capture of images across the entire plate at a resolution of ∼115 nm compared to ∼463 nm for unexpanded samples using the same microscope, objective (63X, 1.15NA water immersion), and imaging wavelength (Fig. S7). Image resolution was calculated using a decorrelation analysis.15 We used the PreciScan plug-in in the Harmony analysis software to first accurately target regions of interest based on Hoechst staining to reduce acquisition times at higher magnification. Briefly, we imaged plates at 5x across entire wells and stored the coordinates of nuclei for each well. These coordinates were then used to direct the 63x magnification objective to the relevant fields of view.

HiExM is compatible with nanoscale expansion and automated image acquisition of human cells with minimal distortion

a. Expanded gels in a 96 well plate imaged at 5x with an inset showing an individual gel shown on the right. b. Imaging fields shown for pre-expansion and post-expansion samples. Prior to expansion, 61 fields were imaged at 63x around the center of the gel. After expansion, 12 fields were imaged at 20x magnification. Max projections for each field were then stitched together and manually inspected to identify fields in the pre-expansion and post-expansion image sets containing matching cells. Arrows show an example of a field representing the same cells before and after expansion. c. Representative registered fluorescence image of microtubules before (left) and after (right) expansion showing the increase in resolution conferred by expansion. In this case, both pre-expansion and post-expansion images were taken at 63x magnification. d. Representative error curve calculated using non-rigid registration of 43 independent fields of view. The shaded region denotes one standard error of the mean. e. Error is not dependent on the location of a given field of view within the gel. Each data point represents the average percent error (up to 40 μm measurement length post-expansion) for a given field of view as determined by an NRR analysis. Seven to ten fields of view were analyzed across six wells. The distance from the center of each field of view to the center of the well was measured manually in FIJI using the stitched pre-expansion image set. f. Box plots showing the distribution of average error among gels across three plates (average = 3.63% +/− 1.39). Average error percentages for each well plate: 3.73 +/− 0.879, 3.89 +/− 1.65, and 3.43 +/− 1.44. g. Expansion factor is not dependent on the location of a given field of view within the gel. Images taken before and after expansion of the same cell or group of cells were used to measure the ratio of length scales in the two images. h. Box plots showing the distribution of expansion factors of 58 gels across three plates. Expansion factors ranged from 3.5-5.1x across gels with an average of 4.16 +/− 0.394. Average expansion factors for each well plate showed tight distribution with some variability across plates: 3.78 +/− 0.171, 4.26 +/− 0.171, and 4.53 +/− 0.321. Comparison of the three well plates was performed using one-way ANOVA and post-hoc Tukey HSD test (*p-value<<0.05).

We next benchmarked the performance of our platform relative to slide-based proExM.9 Image distortion due to expansion was measured using a non-rigid registration (NRR) algorithm.4 Prior to the application of HiExM, images of microtubule stained A549 cells were acquired using the Opera Phenix from 61 fields around the center of each well (Fig. 2a). Cells were then expanded as above and 12 fields were imaged per well. Images from pre- and post-expansion were stitched together and the images for each well were then compared side-by-side to identify matched cells. Although we observed that gels did not maintain a fixed position pre- and post-expansion, close inspection of the images allowed us to align their orientations. Distortion was calculated using the root mean squared error (RMSE) of a 2-dimensional deformation vector field by comparing microtubule morphology in the same cell before and after expansion in six to ten fields in each well across six arbitrarily selected wells (Figs. 2c,d). The average percent error (up to 40 μm measurement length post-expansion) was plotted against the distance from the center of the gel (Fig. 2e). Our analysis found no correlation (R2 = 0.003) between expansion error and distance from the gel center. In fact, we observed only a small number of outliers (>3 s.d.) in any given well (<2 instances per gel across 6 gels). These outliers are not expected to affect biological interpretation. The outliers, which are closer to the center of the gel, likely occur when the device post is not in full contact with the well surface, resulting in a higher volume of gel at the center that resists gel swelling under the device (Fig. S8). Expansion error was also highly consistent when comparing wells across multiple plates (average = 3.68%), comparable to slide-based ExM (Fig. 2f).4,9 We observed only a small number of image fields (7 of 58) with average percent error above 6% (maximum 6.78%). Moreover, calculated error showed that photoinitiation resulted in markedly lower spatial errors compared to APS and TEMED polymerization chemistry (Fig. S9). Analyzing cell seeding densities also showed that cell confluency did not affect image distortion as measured by NRR (Fig. S10).

We next measured expansion factors within individual gels, between gels, and across plates by manually measuring the distance between matched pairs of cell landmarks in FIJI before and after expansion. We found that the expansion factor was highly reproducible across different regions of the gel confirming robust expansion using our method (Fig. 2g). The average expansion factor of gels sampled across three plates was 4.16 +/− 0.394 (Fig. 2h), consistent with proExM.9 We next compared the distributions of expansion factors across 3 well plates by ANOVA and found statistically significant differences between expansion factor distributions. Thus, while uniform within gels, the expansion factor is somewhat variable between gels and plates. We attribute these differences primarily to the small size of the gels, making them vulnerable to the effects of evaporation. This variability should be taken into consideration for studies where absolute length measurements between plates are important for biological interpretation. Collectively, these data show that the HiExM method results in reproducible nanoscale and isotropic expansion of cells within gels, across wells, and plates.

The ability to image cells at increased resolution using HiExM provides an opportunity to analyze nanoscale cellular features for a range of applications including cell phenotyping as well as drug and toxicology studies. As proof of concept, we analyzed the effects of the cardiotoxic chemotherapy agent doxorubicin (dox) on human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) cultured in 96 well plates (Fig. 3a).16 Commonly used dox concentrations in in vitro cardiotoxicity studies range between 0.1 μM to 10 μM, depending on experimental goals1720. Although studies largely focus on electrophysiology, structural changes at the cellular level also provide critical information for assessing cardiotoxicity. For example, changes in nuclear morphology as well as DNA damage are associated with cellular stress and apoptosis, which are commonly observed at high dox concentrations (>200 nM) in in vitro studies using conventional microscopy.1721 In contrast, heart tissue is potentially exposed to lower dox concentrations in clinical treatments.16,22,23 Because HiExM has the sensitivity to detect nanoscale cellular features, we compared nuclei using Hoechst signal intensity in CMs treated with 1 nM, 10 nM, 100 nM, and 1 µM dox after 24 hours compared with DMSO (Fig. 3b). Although we did not observe a significant change in overall nuclear morphology (Fig. S11), images showed increasing Hoechst density at the nuclear periphery in a dox concentration dependent manner.

HiExM shows altered nuclear morphology of cardiomyocytes treated with low doses of doxorubicin

a. hiPSC-CMs imaged post-expansion in HiExM. b. Example images of hiPSC-CM nuclei after dox treatment for 24 hrs before (top) and after (bottom) expansion. c. Example images of hiPSC-CM nuclei taken before (top) and after (bottom) expansion. d. Heatmaps of nuclei in c showing pixel groupings for subsequent analysis. Heatmaps were generated from the image mask after dilation, such that the outer edge represents a contour just beyond the nuclear periphery. e. Relative Hoechst intensity plotted as a function of pixel position relative to the edge of the dilated mask in pre-(top) and post-(bottom) expansion nuclei. Colored curves represent dox concentrations or DMSO (ctrl). Shaded regions represent SEM for n = 56, 71, 64, 92, and 62 nuclei in the pre-expansion images and n = 4 replicates in the expanded case (118, 111, 110, 113, and 77 nuclei analyzed for DMSO control, 1 nM, 10 nM, 100 nM, and 1 µM Dox, respectively). f. Insets of plots in e highlighting the nuclear periphery. g. Rate-of-change analysis of curves for pre-(left) and post-(right) expansion images. Outer edge values were obtained by determining the maximum values of the derivatives for each condition, and inner edge values were obtained by determining the minimum values of the derivatives within the domain 0 - 0.4. h. Curvature analysis for pre-(left) and post-(right) expansion data. Values represent the minima of the second derivative for curves within the domain 0 – 0.4. i. Average Hoechst intensities plotted as gradients for each dox concentration. Scale bar represents 10 percent distance from the edge of the nucleus to the center. Error bars represent SEM as in e. * denotes significance from an independent two-sample t-test (p < 0.05). *Scale bar represents biological scale assuming 4x expansion.

To quantify this effect, we measured the gradient slope of Hoechst intensity around the nuclear periphery. To this end, we masked 529 nuclei from background in dox treated and control samples (Fig. 3c). We selected a single optical slice that intersected the nucleus perpendicularly along the z-axis for our analysis to avoid out-of-focus signal at the nuclear edge (Fig. S12). To account for irregular nuclear shapes, we iteratively grouped pixels into rings for each nucleus starting with the edge of the mask (Fig 3d). The mask is dilated outward from the edge of the nucleus to include a background reading as reference. The average pixel intensity value for each ring represents the average Hoechst intensity as a function of relative distance from the edge of the nucleus to its center (Fig. 3e, Movie S3). Notably, HiExM revealed a peak of Hoechst intensity at the edge of the nucleus that is not resolved in pre-expansion images (Fig. 3f). A line-scan-based analysis showed the same trend (Fig. S13). To further evaluate these differences, we next calculated the first derivative of each curve and found the maximum and minimum values, which correspond to the gradient slopes on either side of the peak (Fig. 3g, Fig. S14). In this case, a higher magnitude of gradient slope reflects a ‘sharper’ transition from the nuclear edge to the background (outside) or from the nuclear edge to the nucleoplasm (inside). Post-expansion, the gradient slopes on both sides of the peak showed a significant increase at dox concentrations as low as 10 nM compared to DMSO, whereas no trend is observed in pre-expansion images (Fig. 3g).

In addition to gradient slopes at the nuclear edge, we also quantified the pattern of Hoechst signal at the nuclear edge by calculating the curvature of the peak using the second derivative (Fig. 3h). These data also showed a significant increase in Hoechst intensity at 10 nM dox compared to DMSO, whereas no significant difference is observed using conventional confocal microscopy (Fig. 3i). We note that, for both the outer gradient slope and peak curvature, values decrease from 100 nM to 1 µM dox, consistent with increased chromatin condensation associated with apoptosis typically observed at higher dox concentrations19,24. In contrast, CMs treated with increasing concentrations of H2O2, another cellular stressor, did not show similar differences in Hoechst patterns (Fig. S15). Although direct evidence specifically linking dox to increased DNA condensation at the nuclear periphery is limited, the known pro-apoptotic effects of dox strongly suggest that our observations correlate with these changes.17,19,25 Thus, we suggest that HiExM can detect early alterations in response to dox treatment using Hoechst signal, setting the stage for further investigation to understand the consequences across a range of drugs and conditions. Overall, our analysis demonstrates the ability of HiExM to detect subtle cellular changes compared to standard confocal microscopy that could improve the resolution of drug screening approaches.

Our work demonstrates that HiExM is a simple and inexpensive method that can increase the throughput of expansion microscopy. Another recent method also reports a scalable ExM-based method with comparable throughput, albeit this approach requires the use of a custom cell culture microplate and manual imaging.8 Our device is fabricated with CNC milling, a common fabrication technique available in most academic settings, making our method widely accessible. The HiExM device can also be produced by injection molding which can further extend its accessibility and distribution in research settings. We have also adapted our device to fit cell culture plates of varying dimensions including a 24 well plate which can accommodate greater expansion and significantly higher resolution (Fig. S16).26 Moreover, the use of a photo-active gel solution enabled consistent expansion in our HiExM workflow. In this work, we found that the use of photoinitiation with irgacure caused signal loss for to varying degrees depending on the fluorophore; coupled with the signal dilution inherent to the expansion process, signal intensity is an important consideration when using HiExM. Our platform is amenable to other photo-active expansion chemistries such as the Photo-ExM chemistry which does not require CF dyes or an anoxic environment (Fig. S17)12,13.

Another advantage of HiExM is that it is fully compatible with autonomous imaging, allowing data acquisition of thousands of cells within a day. HiExM achieved a resolution of ∼115 nm in a 96-well culture plate using a 63x water immersion objective with a 1.15 numerical aperture. By employing microscopes with a higher numerical aperture, we expect our platform can achieve significantly greater resolution. We also show that HiExM error measurements are similar to proExM and suggest that expansion distortions should be considered for robust data interpretation. Recent innovative work that developed a microplate-based ExM protocol maintains hydrogel geometry and orientation improves coordinate-based imaging to accurately assess distortions with ease.27 Another recent work used micropatterning to apply expandable protein grids to coverslips for local expansion factor measurement. This method could be adapted to HiExM to enable expansion measurements in large datasets where non-rigid registration is impractical.28 Taken together, HiExM is an accessible nanoscale imaging platform designed for scalable sample preparation and autonomous imaging, broadening the use of ExM in biological research and pharmaceutical applications.

Materials and methods

Cell culture

Human A549 cells were maintained in culture at a density between 2×103 and 1×104 cells/cm2 in DMEM modified with 10% fetal bovine serum and 1% penicillin-streptomycin. Glass and plastic 96-well plates (Cellvis, P96-1.5H-N or P96-1.5P) were seeded at ∼1×104 and 7×104 cells/cm2 for expansion experiments.

Induced pluripotent stem cell-derived cardiomyocytes were purchased from Cellular Dynamics and cultured on glass-bottom 96-well plates at 30,000 cells/well. For drug treatment, CMs were plated on gelatin for 7 days (Cellular Dynamics) prior to drug treatment and incubated with varying concentrations of Doxorubicin (Tocris 2252) for 24 hours prior to fixation and subsequent sample preparation.

Immunofluorescence

Briefly, cells were fixed with 4% PFA for 15 minutes at room temperature and then blocked (1% bovine serum albumin and 0.1% Triton X100 in PBS) for one hour. Cells were then incubated for one hour at 37°C with primary antibodies in a blocking solution followed by incubation with secondary antibodies in blocking solution for one hour at 37 °C. Cells were washed between each stage in PBS three times at 0-, 5-, and 10-minute intervals.

HiExM Device Fabrication

HiExM devices were fabricated from polysulfone (McMaster-Carr 86735K74) using a Trak DPM2 CNC milling machine using the design provided in the supplementary files.

HiExM Workflow

Cells were fixed with 4% paraformaldehyde (PFA) and immunostained as described above and treated with 50 µM (A549 cells) or 100 µM (cardiomyocytes) Acryloyl-X, SE (6-((Acryloyl)amino)hexanoic acid, succinimidyl ester abbreviated as AcX (Invitrogen A20770)) overnight at 4°C. Cells were then washed twice with PBS, once with DI water, and the remaining liquid was gently aspirated from the wells. Gel solution composed of 4.04M sodium acrylate (AK Scientific R624), 0.506M acrylamide (MilliporeSigma A4058), 0.00934M n, n’-methylenebisacrylamide (MilliporeSigma M1533), 1.84M sodium chloride (Invitrogen AM9759) and 1.45X PBS 1mL (Invitrogen AM9625) was prepared in a light-protected tube. 1-[4-(2-Hydroxyethoxy)-phenyl]-2-hydroxy-2-methyl-1-propane-1-one (Irgacure 2959, MilliporeSigma 410896) was dissolved in DMSO at a concentration of 1.64M and this solution was mixed with the gel solution to a final Irgacure 2959 concentration of 0.0328M at room temperature. The HiExM protocol was performed in a glove bag purged two times with dry nitrogen (MilliporeSigma Z530112). The gel solution was poured onto an absorbent pad at room temperature to form a fluid reservoir. Devices were dipped into the gel solution and immediately inserted into the well plate. The plate was then irradiated for 60 seconds with UVA light using a Uvitron Portaray curing lamp placed 20 inches beneath the plate (Uvitron International UV2270). Approximately 100 µL of Proteinase K (NEB P8107S) digestion solution was added to each gel-containing well and incubated for ∼6 hours at room temperature. After digestion, Hoechst 33324 (Fisher scientific #62249) in DI water (1:5000 dilution) was added to wells for 5 minutes, then the 96-well plate was submerged in ∼4 liters of DI water overnight under constant agitation while the gels were held in place by the device. Finally, the plate was removed from the water bath and devices were removed from the well plate for imaging. For each well, excess water was aspirated and two drops of mineral oil was added to secure the expanded gel within the well.

Image Acquisition

Confocal imaging was performed on an Opera Phenix imaging system (PerkinElmer). For NRR analysis, pre-expansion images were acquired with a 63x water objective (NA 1.15) in stacks of 6 planes spaced 0.5 microns apart and post-expansion images of the same plate were acquired using a 20x water objective (NA 1.0) in stacks of 50 planes spaced one micron apart. Prior to expansion, 61 tiled fields at the center of each well were imaged at 63x to capture cells that would be incorporated into the gel. After expansion, 12 fields in a cross pattern were imaged at 20x to capture a representative set of images for each gel across multiple wells. Based on our analysis, the centers of gels in HiExM should be avoided in imaging as these regions can include cells that are torn and disfigured due to detaching the gels from the device. Movies S1 and S2 were obtained on an EVOS M5000 epifluorescence microscope (ThermoFisher Scientific).

For drug treatment experiments, expanded images were acquired at 63x magnification using the PreciScan plugin. First, each well was imaged for Hoechst at 5x in 4 quadrants to cover the full footprint of the well. An analysis in Harmony determined the coordinates of each nucleus in three dimensions and selected a subset of those coordinates to image again at 63x. The sample sizes for each condition in expanded samples were as follows: 118, 111, 110, 113, and 77 cells for DMSO control, 1 nM, 10 nM, 100 nM, and 1 µM Doxorubuicn (dox), respectively across the four independent experiments. Pre-expansion images were acquired at 63x magnification without the PreciScan plugin. Sample sizes for each condition in pre-expansion samples were as follows: 56, 71, 64, 92, and 62 for DMSO control, 1 nM, 10 nM, 100 nM, and 1 µM Doxorubuicn (dox), respectively.

HiExM samples were imaged on the Opera Phenix at 63x with the following parameters: 100% laser power for all channels; 200 ms exposure for Hoechst, 500-1000+ ms exposure for other channels depending on the strength of the stain and the laser; 60 optical sections with 1 micron spacing; and 4-20 fields of view per well depending on the cell density and sample size requirements. Therefore, imaging one full 96-well plate (60 wells total to avoid outer wells) requires 3 hr to 64 hr depending on the combination of parameters used.

Distortion and Expansion Factor Analysis

Pre-Expansion/Post-Expansion Image Curation

Prior to applying HiExM, 61 fields of view were imaged at 63x centered in each well. After expansion, gels were imaged at 20x in 12 fields of view centered in each well in a cross pattern as in Figure 2b. All images were converted to maximum projections. Images from pre-expansion and post expansion were stitched and examined to find matching cells. For fields of view with matching cells, the corresponding maximum projections for the individual fields of view in pre- and post-expansion were used for subsequent analyses. This image curation is the rate limiting step in our NRR and expansion factor data analyses due to the time required to identify matched cells.

Non-rigid registration

Non-rigid registration was performed to measure isotropic expansion within individual gels and across gels. For analysis represented in Fig. 2d, all fields of view were used in which a cell or group of cells could be identified in both the pre-expansion and the post-expansion image sets. The distance from the center of the analyzed field of view to the center of the well was measured in the pre-expansion image set. For Fig. 2f, a single field of view near the periphery of the gel was arbitrarily chosen to represent the whole gel. Non-rigid registration was performed using a custom MATLAB analysis pipeline as described.4 Briefly, pre and post expansion images are coarsely aligned using the TurboReg plugin for ImageJ.14 This pair of outputs is then histogram normalized and masks are generated to exclude regions with no features by applying a gaussian blur. B-spline-based registration package in MATLAB was used to do non-rigid registration between the images. This procedure was performed on image fields measuring 400px x 400px for all analyses.

The resulting data from a non-rigid registration analysis is a curve that represents measurement error as a function of measurement length. To synthesize these curves into single points (as represented in Figs. 2c-f), The average percent error up to 40 µm measurement length (as measured after expansion) was used. The average percent error is represented as the mean of the ratios of error to measurement length for each individual data point in the error curve up to 40 µm measurement length.

Expansion Factor

To calculate expansion factors, images taken before and after expansion of the same cell or group of cells were compared. In both images, the distance between two easily identifiable features was measured manually in FIJI, and the ratio of those measurements was used to calculate the expansion factor.

Nuclear Periphery Analysis

Image stacks of Hoechst-stained nuclei were compiled and the midplane of each nucleus was manually identified as the plane with the largest nuclear cross-sectional area. A threshold was applied to midplane images to generate masks. These masks as well as their respective raw midplane images were further analyzed using a custom Python script as follows: First, the mask was dilated to encompass background signal outside of the nucleus. Second, the pixels at the edge of the mask were identified and their position and intensity information were stored. Third, pixels in the 4-neighborhood of each pixel from the outer-edge pixel group were identified and their position and intensity information stored (note that only pixels inside the boundary defined by the outer-edge pixel group were stored in this way). This third step was iterated for each new ‘ring’ until all pixels in the dilated mask were grouped. Finally, the mean pixel intensity values of each group were used and the list of mean values was run through linear interpolation to arrive at a list of length 500. To determine the inflection-point slopes of the resulting curves, the first derivative of each curve was calculated using the finite difference approximation. To determine the curvature at the peak intensities of each curve, the second derivative was calculated in the same way. Significance was determined using a two-sample t-test (p < 0.05).

Data Sharing Plans

Solidworks files for our devices as well as code for our analysis are available at the following link: https://github.com/lboyerlab/hiExM_Supplementary_Files.

Acknowledgements

We thank Paul Tillberg and members of the Boyer and Boyden labs, especially Shiwei Wang, Chi Zhang and Asmamaw Wassie for helpful discussions. We also thank the Edgerton Student Shop and Mark Belanger for assisting in the fabrication of devices. We thank the Center for the Development of Therapeutics at the Broad Institute and acknowledge funding through the S10 grant NIH OD-026839-01. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. 2141064 and 1122374 to JHD and TS, respectively. This work is supported by NIH EB024261, NIH 1R01AG070831, HHMI, NIH 1R01MH123403, NIH R01MH124606, and NIH 1R56AG069192 to ESB. This research is supported in part by the National Research Foundation, Prime Minister’s Office, Singapore under its Campus for Research Excellence and Technological Enterprise (CREATE) programme, through Singapore MIT Alliance for Research and Technology (SMART): Critical Analytics for Manufacturing Personalised-Medicine (CAMP) Inter-Disciplinary Research Group. ESB acknowledges funding from Lisa Yang, Lore McGovern and John Doerr. This work is also supported by funding from NIH R01HL140471, the Leila and Harold G. Mathers Foundation, John J. Jarve Seed Fund and the Deshpande Center to LAB.

Disclosures

ESB is an inventor on multiple patents related to expansion microscopy, and co-founder of a company exploring commercial applications of expansion microscopy. JHD, PM, CMDS, ALA, ESB, and LAB are inventors on a provisional patent filed for the HiExM device.