Quantitative mapping of dynamic 3D transport in growing cells via volumetric spatio-temporal image correlation spectroscopy (vSTICS)

Abstract: Quantitatively mapping three-dimensional (3D) flow, diffusion, and particle density in crowded living cells remains challenging because most dynamic optical microscopy measurements are effectively planar and existing analysis methods struggle with dense, noisy volumetric data. We introduce volumetric spatio-temporal image correlation spectroscopy (vSTICS), a framework that recovers voxel-resolved flow, diffusion coefficients, and particle densities from 3D fluorescence time series. Growing Camellia japonica pollen tubes were imaged with field-synthesis lattice light-sheet microscopy, and localized 3D spatio-temporal correlation analysis was applied to overlapping volumetric samples to generate maps of velocity, diffusion, and density. Validation with synthetic flow-diffusion simulations showed accurate recovery of seeded transport parameters, including velocities near $3$ $μ$m s$<sup>{-1}$ and diffusion near $10<sup>{-3}$ $μ$m$<sup>2$ s$<sup>{-1}$. Fluorescent microsphere experiments verified particle number and point spread function readouts and measured diffusion coefficients of $0.3 \pm 0.1$ $μ$m$<sup>2$ s$<sup>{-1}$ in gel, consistent with imaging-FCS measurements of $0.5 \pm 0.2$ $μ$m$<sup>2$ s$<sup>{-1}$. Applied to mitochondria in pollen tubes, vSTICS resolved a bidirectional reverse-fountain pattern with slower anterograde transport ($0.1$-$1$ $μ$m s$<sup>{-1}$) and faster retrograde motion peaking near $3$ $μ$m s$<sup>{-1}$, plus a retrograde corridor about $2$ $μ$m wide. Density and diffusion maps indicated a denser, more advective core and higher peripheral diffusion. High-density sub-diffraction vesicle mapping produced similar velocity landscapes with about ten-fold higher particle densities. These results establish vSTICS as a practical method for quantitative 3D mapping of intracellular transport and refines the reverse-fountain model by revealing asymmetric, predominantly transverse circulation.

• The paper introduces vSTICS, a volumetric image-correlation method that extracts voxel-level velocity, diffusion, and particle density from 3D fluorescence movies, avoiding trajectory tracking in crowded samples.
• The method recovered simulated transport values closely, matched microsphere and imFCS diffusion benchmarks, and mapped pollen-tube mitochondrial flow with retrograde speeds near 3 μm/s—about three times anterograde speeds.
• The study shows that vesicles at densities near 30 particles/μm³ remain analyzable with vSTICS, while PSF anisotropy, fluctuation sampling, and stationarity constrain axial velocity and diffusion measurements.

Overview and motivation

Quantitative characterization of intracellular transport has historically been constrained by two limitations: most optical measurements of flow are effectively planar, and existing analysis methods degrade in dense, noisy volumetric data. Single-particle tracking (SPT) fails when particle density is high and trajectories overlap; particle image velocimetry (PIV) provides ensemble velocity fields but is insensitive to diffusion and fragile at low signal-to-noise ratio. Spatio-temporal image correlation spectroscopy (STICS) addresses some of these shortcomings but had been applied almost exclusively to quasi-2D systems (2603.27484).

The paper introduces volumetric spatio-temporal image correlation spectroscopy (vSTICS), a framework that recovers voxel-resolved velocity vectors, diffusion coefficients, and particle densities directly from 3D fluorescence time series. The method is paired with field-synthesis lattice light-sheet microscopy (LLSM), chosen for its high volume rates and reduced photobleaching. Validation proceeds through synthetic flow–diffusion simulations and fluorescent microsphere calibrations, followed by application to mitochondria and sub-diffraction cargo vesicles in growing Camellia japonica pollen tubes—a bidirectional transport system that serves as a stringent biological benchmark.

Mathematical framework

vSTICS computes the normalized spatial intensity fluctuation autocorrelation function on overlapping volumetric regions of interest (ROIs) across a sequence of temporal lags $\tau$. The molecular detection efficiency (MDE), modeled as a 3D Gaussian with lateral radius $\omega_{xy}$ and axial radius $\omega_z$, transforms concentration distributions into fluorescence signals. The non-normalized autocorrelation is evaluated efficiently via the Wiener–Khinchin theorem as the inverse Fourier transform of the image power spectrum.

At each lag, the spatio-temporal correlation function is fit with an asymmetric 3D Gaussian whose parameters encode all three transport quantities simultaneously:

• Velocity: the fitted peak positions $(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$ drift linearly with $\tau$, and linear fitting yields local velocity vectors.
• Density: the zero-lag amplitude $g(\vec{0}) = 1/\langle N \rangle$ is inversely proportional to the mean number of particles in the effective observation volume $V_{eff} = \pi^{3/2}\omega_{xy}^{2}\omega_z$.
• Diffusion: decay of the normalized correlation amplitude with $\tau$, fit with a 3D diffusion model, yields $D$.

A notable methodological choice is estimating $D$ from amplitude decorrelation rather than time-dependent peak broadening (as in iMSD), which the authors report is more robust in strong-flow regimes where PSF-scale decorrelation remains accessible at LLSM volume rates.

Validation with simulations and calibration standards

Simulations seeded particles undergoing pure flow, combined flow and diffusion, and combined dynamics with additive noise. With ground-truth values of $\omega_{xy}$0, $\omega_{xy}$1, and density $\omega_{xy}$2 particles/$\omega_{xy}$3m³, vSTICS recovered $\omega_{xy}$4, $\omega_{xy}$5, $\omega_{xy}$6, and density $\omega_{xy}$7 particles/$\omega_{xy}$8m³. Axial recovery shows somewhat larger error bars than transverse components, foreshadowing the PSF-anisotropy limitation discussed below.

Experimental calibrations used agar-immobilized 50 nm fluorospheres and fluorospheres suspended in 0.3% Gelrite (viscosity raised from 1 cP to 8 cP, reducing the theoretical Stokes–Einstein diffusion coefficient to 0.5 $\omega_{xy}$9m²/s). Key results:

Measurement	vSTICS	Independent reference
Particle count (stationary beads)	$\omega_z$0	$\omega_z$1 (segmentation, full beads only)
Structure factor $\omega_z$2	$\omega_z$3	$\omega_z$4 (segmentation)
Diffusion (beads in gel)	$\omega_z$5	$\omega_z$6 (imFCS)

The discrepancy between vICS and segmentation particle counts arises because segmentation counts only fully contained microspheres while correlation analysis reports partial objects within ROI crops; the authors treat this as consistent rather than contradictory. The agreement between vTICS and imFCS diffusion values confirms quantitative accuracy within the imaging bandwidth, though the vTICS estimate sits slightly below the theoretical value.

Mitochondrial transport: an asymmetric reverse fountain

Applying vSTICS to MitoTracker Green-labeled mitochondria in growing pollen tubes (55 × 55 × 10 µm volumes at 0.5 s per volume, 125 time points) produced simultaneous maps of density, diffusion, and flow. Three findings stand out:

Asymmetric bidirectional flow. Anterograde motion ranges from 0.1 to 1 $\omega_z$7m s⁻¹ while retrograde transport peaks near 3 $\omega_z$8m s⁻¹, yielding a retrograde-to-anterograde velocity ratio of approximately 3—contradicting prior symmetric-flow assumptions in the literature. Transverse velocity histograms are bimodal, whereas the axial component is narrowly centered at zero ($\omega_z$9m s⁻¹), indicating circulation is predominantly confined to transverse planes. This refines the classic "reverse-fountain" model: the circulation is largely transverse rather than axially symmetric.

Geometric compartmentalization. Skewed-Gaussian fits to transverse speed profiles yield a retrograde corridor FWHM of $(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$0m ($(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$1). For tubes of diameter $(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$2m, the anterograde corridor occupies roughly 20% of the tube diameter. These geometric constraints can be incorporated directly into mechanochemical models of tip growth.

Coupled density–diffusion gradients. Mitochondrial densities average $(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$3 particles/$(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$4m³ overall, with a denser core (~6 particles/$(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$5m³ center vs. ~4 at the periphery). Diffusion coefficients span 0.1–1 $(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$6m²/s, systematically higher at the periphery. Combined with the velocity field, this pattern implies a Péclet number $(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$7 near the core (advection-dominated) and more diffusive behavior toward the edges—consistent with cytoskeletal architecture imposing transverse compartmentalization of advective versus stochastic motion.

Vesicle transport at high density

Cargo vesicles (~100 nm diameter, sub-diffraction-limited) labeled with FM 1-43 present densities roughly tenfold higher than mitochondria (~30 vs. ~3 particles/$(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$8m³), placing them beyond the reach of SPT or PIV. After segmenting away the plasma membrane with Ilastik, vSTICS recovered density maps showing enrichment near the growing tip and along retrograde corridors, together with bimodal velocity landscapes similar to the mitochondrial case. The apex-region vesicle density of ≈30 $(\mu_{\xi}, \mu_{\eta}, \mu_{\phi})$9m⁻³ is consistent with earlier FRAP-based estimates in lily pollen tubes, providing independent cross-validation. Notably, the correlation amplitude decay was insignificant for vesicles, so no diffusion coefficient is reported—an honest limitation of the fluctuation-sampling bandwidth for these conditions.

Practical parameterization

The paper supplies concrete guidance for adapting vSTICS: ROI and TOI tiling must balance spatial resolution against correlation signal-to-noise. The working configuration used $\tau$0 px, $\tau$1 px, shifts of 5 and 2 px, and a TOI of 5 frames with unit shift. The authors emphasize that recovered parameters depend on detector noise, background normalization, axial sampling, temporal windowing, and rolling z-stack acquisition, and recommend simulation-based characterization as a mandatory first step on any new platform or specimen. All code is publicly available on GitHub.

Limitations and open questions

The most significant limitation is PSF anisotropy: large structure factors $\tau$2 reduce sensitivity to $\tau$3, inflate uncertainty in axial peak localization, and bias velocities toward the imaging plane. Field synthesis mitigates anisotropy in calibration samples, but in situ refractive-index mismatch and scattering can enlarge $\tau$4. Proposed remedies include tighter light-sheet confinement, deconvolution with measured PSFs, and smaller $\tau$5 with increased axial sampling.

Interpretation also rests on assumptions that may fail in practice: locally stationary dynamics within each time-of-interest window and adequate fluctuation sampling. Rapidly varying growth or flow demands shorter TOIs and higher volume rates; slow diffusion requires longer TOIs, which in turn require mechanical stability over the sampling interval. An open question left by the work is how well the amplitude-decay approach to diffusion estimation performs for very high-density species like vesicles, where decorrelation was not resolvable here.

Conclusion

vSTICS extends image correlation spectroscopy from single-plane measurements into full volumetric transport mapping, recovering density, diffusion, and directed flow from a single-pass analysis of 3D fluorescence time series. Validated against simulations, microsphere calibrations, and imFCS benchmarks, it was applied to living pollen tubes to reveal an asymmetric, predominantly transverse reverse-fountain circulation with a retrograde-to-anterograde velocity ratio near 3 and a quantified retrograde corridor of about 2 µm. The framework, together with its parameterization guidelines and open-source implementation, constitutes a transferable template for quantitative 3D transport mapping in crowded biological systems such as axons, cilia, ducts, and microvascular networks, provided that PSF anisotropy and stationarity assumptions are carefully managed.