For the following experiments, a microfluidic compression device15 was used to model tumor-immune interactions within a physiologically relevant 3D microenvironment, and at the same time, can apply static and dynamic compressions to tumor spheroids (Figure 1). The device components (Figure 1A1) and their assembly order (Figure 1A2) are shown in Figure 1A. The metal frame serves as the base of the device, with the cell chamber layer (L1) positioned within the frame. The pressure control unit, consisting of L2 and L3 plasma bonded together, is placed above L1, followed by the manifold (Figure 1A2). Together, L1 and the pressure control unit form the functional microfluidic device and are sandwiched between the metal frame and the manifold with screws. The manifold contains inlet ports for fluidic tubing to supply pressure to the device, as well as media reservoirs that facilitate long-term live-cell imaging experiments.

Figure 1: Microfluidic compression device for studies of interactions between immune cells and tumor spheroids under mechanical compression. (A1) Photograph of the components required for device assembly, shown from left to right: manifold, pressure control unit (L2 + L3), cell chamber layer (L1), and metal frame. (A2) 3D schematic of the device components, shown from top to bottom: manifold, pressure control unit (L2 + L3), cell chamber layer (L1), and metal frame. PDMS layers L2 and L3 are plasma-bonded to form the pressure control unit, which is positioned above the cell chamber. The manifold and metal frame secure the assembled device. (B1) Cross-sectional schematic of a compression unit before and after compression. ECM-embedded tumor spheroids and NK cells are placed in the cell chamber. The initial distance between the bottom of L1 and the bottom of L2 before compression is denoted as h. Upon application of pressure to the pressure control unit, L2 deflects downward, causing the PDMS piston to compress the tumor spheroid. The downward displacement of L2 following compression is denoted as Δh. (B2) Merged brightfield and GFP images of MCF-7 tumor spheroid with NK cells before (left) and after (right) compression. MCF-7 cells are GFP-labeled, while NK cells are not labeled. Representative NK cells are circled in blue, while yellow outlines mark the peripheries of the tumor spheroids. Scale bar = 50 µm. Please click here to view a larger version of this figure.
Once assembled, the microfluidic compression device can be used to compress tumor spheroids as illustrated in Figure 1B. Here, the collagen-embedded tumor spheroids with NK cells are introduced into the cell chamber (L1). Upon application of pressure to the pressure chamber (L3), the deformable membrane layer (L2) deflects downward, causing the Polydimethylsiloxane (PDMS) piston to transmit compressive forces to the spheroid. The pressure chamber is controlled by an external pressure source. Additional details regarding the design, fabrication, and characterization of the microfluidic compression device have previously been published in another study15.
We define the distance between the bottom of the cell chamber (L1) and the bottom of the deformable membrane (L2) as h (Figure 1B1, left panel). Before compression, this distance is 200 µm. Upon application of fluid pressure to the pressure control unit, the deformable membrane deflects downward, reducing the distance h by Δh (Δh; Figure 1B1, right panel) and thereby compressing the collagen-embedded tumor spheroid. In this system, compressive strain is defined as Δh/h, where h is the initial distance between L1 and L2 before compression, and Δh is the downward displacement of L2 after compression. Membrane displacement was previously calibrated by visualizing fluorescent beads adhered to the bottom of L2. A detailed description of the calibration procedure and compressive strain calculations can be found in a previous study17. An example calibration curve can be seen in Supplementary Figure 2.
Figure 1B2 shows a representative merged brightfield and GFP images of an MCF-7 spheroid with NK cells embedded in 1.5 mg/mL collagen before (left panel) and after compression (right panel). Tumor spheroid boundaries were outlined in yellow. Following compression, the spheroid exhibited an increase in projected area, indicating deformation in response to the successful application of mechanical load. For all experiments presented in this manuscript, we used a compression strain of Δh/h = 50%. We note that the spheroid diameter ranges from 140 µm to 200 µm. This means that the actual compression strain on the spheroid is slightly smaller than 50%. Supplementary Video 1 and Supplementary Video 2 show representative X-Y z-stacks acquired in the GFP fluorescence channel of an MCF-7 spheroid embedded in 1.5 mg/mL collagen before and after application of the 50% compressive strain, respectively.
The microfluidic compression device was used to investigate the effects of NK cells, compression, and their combination on tumor spheroid death (Figure 2). This is made possible by the design of the device, which contains 12 culture wells. Of these wells, six can be subjected to compression, while the remaining six serve as uncompressed controls. This configuration enables four experimental conditions (Figure 2A): spheroids alone with no compression (sp), spheroids with NK cells but no compression (+NK), spheroids subjected to compression (+comp), and spheroids exposed to both compression and NK cells (+comp+NK).

Figure 2: Mechanical compression modulates NK cell–mediated killing of tumor spheroids. (A) Time-lapse fluorescence images of MCF-7 spheroids embedded in 1.5 mg/mL collagen under four conditions: spheroid-only control (sp), spheroids with NK cells (+NK), spheroids with compression (+comp), and spheroids with both compression and NK cells (+comp+NK). GFP fluorescence labels viable MCF-7 tumor cells. The yellow outline indicates the region of interest (ROI), which was defined using the 18-h image and applied to all time points for fluorescence quantification. Total GFP fluorescence intensity within the ROI was measured and used for subsequent analysis. Scale bar = 50 µm. (B) Total GFP fluorescence intensity of one spheroid over time under 4 conditions. Each colored line represents the total fluorescence intensity of a single spheroid, normalized to its initial intensity at t = 0. The black line is the second-order polynomial fitting of the fluorescence intensity data pooled from 7 spheroids cultured under the control condition and was used as the reference for calculating tumor cell loss in (C). (C) Normalized GFP loss for the +comp, +NK, and +comp+NK conditions was calculated by subtracting the GFP intensity at t = 18 h of each spheroid from the spheroid-only reference line (sp). Therefore, each data point represents the normalized GFP loss of a single spheroid. Error bars represent SEM. Please click here to view a larger version of this figure.
The effects of NK cells, compression, and their combination on tumor spheroids can be visualized in the representative images of MCF-7 spheroids in 1.5 mg/mL collagen (Figure 2A). Because the MCF-7 cells used to generate the tumor spheroids were GFP-labeled, tumor viability could be monitored through changes in GFP fluorescence. GFP and brightfield images were acquired every 10 min for 18 h using an inverted epifluorescence microscope. Notably, spheroids subjected to compression exhibited a larger projected area than uncompressed spheroids, confirming successful transmission of mechanical load from the deformable membrane to the tumor spheroid. Over time, a progressive loss of GFP fluorescence was observed in the +NK and +comp+NK conditions (Figure 2A), indicating NK cell-mediated tumor cell killing.
We use the loss of GFP fluorescence as a measure of tumor cell death. The GFP-labeled MCF-7 cells used in this study were generated through lentiviral transduction. Briefly, a lentiviral vector was used to introduce an enhanced green fluorescent protein (eGFP) expression cassette driven by the elongation factor-1α (EF1α) promoter, along with a puromycin resistance gene driven by the Rous sarcoma virus (RSV) promoter. Following infection, the lentiviral expression cassettes were integrated into the genomic DNA of the MCF-7 cells, resulting in stable GFP expression. Puromycin selection was subsequently used to eliminate non-transduced cells. The GFP labeling of MCF-7 cells enabled tumor cell viability to be quantified by monitoring changes in GFP fluorescence intensity over time (Figure 2B). A strong correlation between GFP fluorescence loss and established markers of cell death, including the live/dead marker propidium iodide (PI) and the apoptotic marker annexin-V, has previously been established18, validating GFP fluorescence loss as a surrogate measure of tumor cell death.
To quantify the effects of the four experimental conditions on tumor viability, the total GFP fluorescence intensity of each tumor spheroid was measured at every time point, background-subtracted, and normalized to its initial value (Figure 2B). Total spheroid fluorescence was calculated within a predefined region of interest (ROI). This ROI was obtained by manually outlining the boundary of the tumor spheroid in the 18-h image (yellow outline on 18 h panel, Fig. 2A) using an in-house MATLAB script. The total fluorescence intensity within this ROI was then calculated for each time point and corrected by subtracting the average background fluorescence measured from a separate ROI selected elsewhere in the image. The resulting background-corrected fluorescence values were subsequently normalized to the initial fluorescence intensity of the spheroid.
Each colored curve represents a normalized GFP intensity of an individual spheroid under a specific experimental condition in Figure 2B. The curve is obtained using a second-order polynomial fitting to the experimental data. The purple curves represent spheroids subjected to compression alone (+comp, Figure 2B), the red curves represent spheroids with NK cells (+NK, Figure 2B), and the blue curves represent spheroids exposed to both compression and NK cells (+comp+NK, Figure 2B). In contrast, the black curve is the second-order polynomial fit to the data of seven spheroids cultured under the spheroid-only condition serving as a control (sp, Figure 2B). In Figure 2B, only the polynomial fitted data, but not the experimental data points, are presented for clarity. Notably, the GFP fluorescence intensity of the spheroid-only (sp) and spheroid + compression (+comp) conditions exhibited a general upward trend over time (Figure 2B), consistent with continued tumor growth or maintenance of viability. In contrast, the +NK and +comp+NK conditions demonstrated a progressive decrease in GFP fluorescence intensity (Figure 2B), indicating tumor cell death.
Cumulative tumor cell loss was quantified by subtracting the fitted normalized GFP fluorescence intensity of each spheroid under a given condition (+comp, +NK, or +comp+NK) from the control condition at 18 h (Figure 2C). As a result, each data point represents the normalized GFP loss of an individual spheroid relative to the spheroid-only control (Figure 2C). Accordingly, the +comp condition reflects GFP loss resulting from compression alone (purple, Figure 2C), the +NK condition reflects GFP loss resulting from NK cells alone (red, Figure 2C), and the +comp+NK condition reflects GFP loss resulting from the combined effects of compression and NK cells (blue, Figure 2C). The +comp condition exhibited the lowest degree of GFP loss (Figure 2C), suggesting that compression alone had a minimal effect on spheroid viability. In contrast, both the +NK and +comp+NK conditions showed significantly greater GFP loss than the +comp condition (Figure 2C), indicating that NK cells are the primary drivers of tumor cell death. Although the +comp+NK condition exhibited a greater average GFP loss than the +NK condition, this difference did not reach statistical significance (Figure 2C). These findings suggest that compression does not impair NK cell-mediated killing of tumor cells. Taken together, these findings demonstrate that the microfluidic compression device can be successfully used to quantify the independent and combined effects of compression and NK cells on tumor spheroid viability within a physiologically relevant 3D environment.
In contrast, suboptimal outcomes may arise when (1) spheroids are located outside the piston boundary, (2) spheroids are positioned too close to the piston boundary, and (3) multiple spheroids are located within one spheroid diameter of each other (Figure 3). In these cases, the corresponding fields should be excluded from analysis.

Figure 3: Representative example of a sequence of images of spheroids excluded from analysis. Time-lapse brightfield (BF), GFP fluorescence, and merged images of MCF-7 spheroids with NK cells in 1.5 mg/mL collagen before (left column) and after compression (right two columns) are shown. Pre-compression images were acquired 20 min before pressure application and 10 min before the start of time-lapse imaging. Images were subsequently collected every 10 min for 18 h, with 0 h corresponding to the first image acquired immediately after compression and 18 h representing the end of the experiment. Blue arrow indicates the edge of the PDMS piston (L2). This imaging field was excluded from analysis because two spheroids were present within the compression region and were located too close to the piston boundary. As a result, spheroid deformation and GFP fluorescence intensity could not be reliably quantified during compression. Scale bar = 50 µm. Please click here to view a larger version of this figure.
(1) Spheroids located outside the piston boundary do not experience compression because pressure is applied only within the piston region. These spheroids should therefore be excluded from both compressed and control conditions. (2) Spheroids positioned on or near the piston boundary (Figure 3, top spheroid) should be excluded, as the spheroid is experiencing nonuniform mechanical loading. (3) Spheroids located within one spheroid diameter of neighboring spheroids should be excluded (Figure 3). Close proximity between spheroids complicates the quantification of GFP fluorescence intensity because the region of interest (ROI) is defined using the final time point. Over time, spheroid boundaries may become indistinguishable, making it challenging to delineate an ROI that accurately encompasses only a single spheroid. These issues can be minimized by ensuring that spheroids are centered within the compression region during device loading. Careful attention to spheroid distribution, positioning, and device setup is essential for achieving uniform compression and generating reliable, reproducible experimental results.
| Components | Concentration |
| α-MEM media | |
| Fetal bovine serum (FBS) | 10% |
| Penicillin/streptomycin | 1% |
Table 1: Complete MCF-7 medium.
| Components | Concentration |
| α-MEM media | |
| Myo-inositol | 0.2 mM |
| 2-mercaptoethanol | 0.1 mM |
| Fetal bovine serum (FBS) | 12.50% |
| Folic acid | 0.02 mM |
| Horse serum (HS) | 12.50% |
| Penicillin/streptomycin | 1% |
Table 2: NK complete medium.
Supplementary Figure 1: Schematic of the pressure control system. An air compressor supplies compressed air (up to 20 psi) to the pressure controller. The pressure controller regulates the air pressure delivered to the pressure reservoir based on user-defined settings. The controller has a maximum operating pressure of 200 mbar, a pressure stability of 0.005% full scale (FS), a pressure resolution of 0.003% FS, and a response time of 9 ms, enabling accurate and reproducible pressure control. The controller is connected to a computer via universal serial bus (USB), allowing the user to program the magnitude and profile of the applied pressure using the accompanying software. The software supports multiple compression waveforms, including constant, sinusoidal, triangular, and periodic compression. In this schematic, the pressure is set to 120 mbar, and compressed air at this pressure is delivered to the pressurized reservoir through the pneumatic tube. The pressurized reservoir contains 50 mL of deionized (DI) water and serves to convert pneumatic pressure into fluidic pressure. Once pressurized, the reservoir drives fluid through fluid delivery tubing to the pressure control unit of the microfluidic compression device (see Figure 1 for details of the device), where the fluidic pressure is used to apply compression to the tumor spheroid. Note that fluid pressure is used to avoid air bubbles in the PDMS device. Please click here to download this file.
Supplementary Figure 2: Calibration curve for the microfluidic compression device showing the relationship between piston displacement (, µm) and the pressure applied to the pressure control chamber (KPa). Experimental measurements are represented by dots, while the solid line corresponds to the COMSOL simulation. Error bars indicate the standard error of the mean (SEM) from eight independent experimental measurements but are too small to be visible. This figure is reproduced from Suh et al.15 with permission from the Royal Society of Chemistry. Please click here to download this file.
Supplementary Video 1: X-Y z-stack of a GFP-labeled MCF-7 spheroid before compression. The spheroid was embedded in 1.5 mg/mL collagen within the microfluidic compression device and imaged using the GFP fluorescence channel with a 1 µm z-step size. Then, ImageJ was used to enhance contrast. Scale bar = 50 µm. Please click here to download this file.
Supplementary Video 2: A compressive strain of 50% (Δh/h = 50%) was then applied to the GFP-labeled MCF-7 spheroid embedded in 1.5 mg/mL collagen. Following compression, an X–Y z-stack was acquired in the GFP fluorescence channel using a 1 µm step size and processed in ImageJ to enhance contrast. Scale bar = 50 µm. Please click here to download this file.