#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Fast quantitative time lapse displacement imaging of endothelial cell invasion


Authors: Christian Steuwe aff001;  Marie-Mo Vaeyens aff002;  Alvaro Jorge-Peñas aff002;  Célie Cokelaere aff001;  Johan Hofkens aff003;  Maarten B. J. Roeffaers aff001;  Hans Van Oosterwyck aff002
Authors place of work: Centre for Membrane Separations, Adsorption, Catalysis and Spectroscopy for Sustainable Solutions (cMACS), Department of Microbial and Molecular Systems (MS), KU Leuven, Leuven, Belgium aff001;  Biomechanics Section (BMe), Department of Mechanical Engineering, KU Leuven, Leuven, Belgium aff002;  Molecular Imaging and Photonics, Department of Chemistry, KU Leuven, Leuven, Belgium aff003;  Prometheus, Division of Skeletal Tissue Engineering, KU Leuven, Leuven, Belgium aff004
Published in the journal: PLoS ONE 15(1)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0227286

Summary

In order to unravel rapid mechano-chemical feedback mechanisms in sprouting angiogenesis, we combine selective plane illumination microscopy (SPIM) and tailored image registration algorithms — further referred to as SPIM-based displacement microscopy — with an in vitro model of angiogenesis. SPIM successfully tackles the problem of imaging large volumes while upholding the spatial resolution required for the analysis of matrix displacements at a subcellular level. Applied to in vitro angiogenic sprouts, this unique methodological combination relates subcellular activity — minute to second time scale growing and retracting of protrusions — of a multicellular systems to the surrounding matrix deformations with an exceptional temporal resolution of 1 minute for a stack with multiple sprouts simultaneously or every 4 seconds for a single sprout, which is 20 times faster than with a conventional confocal setup. Our study reveals collective but non-synchronised, non-continuous activity of adjacent sprouting cells along with correlations between matrix deformations and protrusion dynamics.

Keywords:

Light – Collagens – Fluorescence microscopy – Fluorescence imaging – Deformation – Optical lenses – Gels – Angiogenesis

Introduction

Traction Force Microscopy (TFM) is a powerful tool to study cell-matrix mechanical interactions, as it allows to quantify unknown cellular forces (tractions) from the displacements (deformations) they induce in the extracellular matrix, provided the mechanical (elastic) properties of the matrix are known. While the technique is well established for 2D cell culture systems (i.e. for cells cultured on relatively soft, elastic substrates) [14], resolving tractions for 3D cultures (with cells being embedded in a matrix-mimicking hydrogel) is computationally much more challenging. It requires the use of numerical methods (such as the finite element method) to solve the inverse mechanical problem that underlies TFM and that relates unknown tractions to measured displacements [57]. The fact that the measured displacements contain errors (e.g. due to a limited number of sampling points and due to microscopy limitations), renders the inverse problem ill-posed, thereby requiring proper regularization to avoid error propagation from displacements to tractions [8]. One can avoid the computationally challenging step of 3D traction recovery (or 4D, when also resolving in time) by only solving for matrix deformations (by means of image registration) and expressing cell-matrix mechanical interactions in terms of matrix deformation-based metrics [9]. This so called displacement microscopy–combining live cell optical microscopy with image registration—has the additional advantage that mechanical characterisation of the matrix is not required (provided all experiments are performed in one and the same matrix).

Traction force and displacement microscopy techniques typically rely on point scanning techniques. Much faster imaging can be provided with Selective Plane Illumination Microscopy (SPIM). SPIM [1012] uses a sheet of light to illuminate samples, as proposed by Voie et al. [13] and further improved by the Stelzer group [14,15]. While SPIM is generally regarded as a fast scanning technique compared to point scanning techniques [16,17], recent studies demonstrate even faster setups by using for example multiple sheets [18], and many extensions have been developed such as Limited-View LSM [19] and Tiling Lightsheet SPIM [20]. Fast and efficient illumination restricted to the observed focal volume additionally minimizes photobleaching and phototoxicity [21,22]. At the same time, live cell imaging with SPIM imposes design challenges [2325] on the sample geometry that needs to be tailored to the SPIM setup. SPIM has been typically applied for imaging large biological samples such as embryos [14,23,26,27], large tissues such as mouse brain samples [28] and tissue mimics such as spheroids, organoids or cellular aggregates [24,25,2931], with only a few studies reporting on single-cell SPIM [3234].

So far SPIM has not been combined with any traction force or displacement microscopy, nor has it been used to study angiogenic sprout dynamics. In the context of angiogenesis—the growth of blood vessels from the pre-existing vasculature–only a few studies exist that have reported on displacements around angiogenic sprouts, again using point scanning techniques for live cell imaging [35,36]. However, fast interactions between matrix and sprout—with minute to second scale protrusions and retractions [37,38]—in larger multicellular volumes remain difficult to monitor with point scanning techniques because of the relatively slow scanning, but are vital for understanding sprout progression [39].

To fill this methodological gap, we developed a SPIM-compatible in vitro model of angiogenesis for applying SPIM-based displacement microscopy. While the signal leading to sprout initiation is a pro-angiogenic chemical signal—with Vascular Endothelial Growth Factor (VEGF) being the most prominent one [40]—it is well known that endothelial cells respond to mechanical cues (such as matrix stiffness) as well and need to apply forces to their environment in order to move forward [4147]. The importance of these cell-matrix mechanical interactions for sprout progression explains the rationale for studying them by means of displacement microscopy. By applying SPIM-based displacement microscopy analysis to the in vitro model of angiogenesis, we demonstrate that sprouts do neither continuously nor synchronously exert forces on their surrounding extracellular matrix. Finally, we focus on short term dynamics of growing protrusions and validate that fast volumetric imaging is vital for following the sprout–matrix interaction during these processes at the minute to second time scale.

Methods and materials

The study was approved by the Medical Ethics Commission of the KU Leuven University Hospitals (file number S54744).

First, to gain imaging speed and to increase the monitored volume, we replace point scanning imaging techniques [5,48] by fast plane scanning using orthogonal lenses as offered by SPIM (see Fig 1A and S1 and S2 Figs for the used SPIM setup). Next, we develop a robust, SPIM-compatible sample design for in situ 3D endothelial cells invading and forming sprouts into a collagen hydrogel with embedded fluorescent microspheres (beads) as fiducial markers (Fig 1B as well as S3 Fig). In this in vitro model sprouting is induced by sphingosine-1-phosphate (S1P), a pro-angiogenic factor that is added to the collagen matrix (as well as other pro-angiogenic growth factors that are present in the EGM2 culture medium). Sprouting tip cells then laterally inhibit neighbouring cells to become tip cells as well. The amount of sprouting (i.e. the amount of endothelial cells that turn from a quiescent to a migratory phenotype) is controlled by a number of signaling mechanisms that regulate lateral inhibition, one of the most studied ones being Dll4-Notch signaling [40,49,50]. Finally, we tailor displacement microscopy algorithms to calculate both absolute and incremental displacements (Fig 1C and 1D) of the collagen matrix. Due to its capability to deal with large distortions in complex materials[51], Free Form Deformation-based non-rigid image registration is used for both approaches to compare images of the fluorescent beads embedded in the collagen and extract the corresponding matrix displacement fields. Absolute displacements (tref = tfix) result from the registration and comparison of the image stack of fluorescent beads acquired under cell activity at each time point ti against the (fixed reference) image stack acquired in a stress-free state of the hydrogel at timepoint tfix. The latter was obtained at the end of the experiment by adding Cytochalasin D [6,52] (CytoD), a cell-permeable mycotoxin, disrupting the cells’ actin cytoskeleton and relaxing the sample (see S4 Fig and S1 Video). At the same time, incremental displacements (tref = ti−τ) can be computed by registering image stacks of fluorescent beads acquired at sequential time points (non-fixed reference) spaced by a time interval τ. (See S5 Fig for notations and comparison of absolute versus incremental displacements). We highlight that both approaches are complementary: while absolute displacements are indicative for the direction and intensity with which cells are mechanically interacting with the matrix, the incremental displacements better reveal the temporal dynamics with which these interactions change. All methods are elaborated below in more detail.

Workflow of 4D displacement microscopy around <i>in vitro</i> angiogenic sprouts.
Fig. 1. Workflow of 4D displacement microscopy around in vitro angiogenic sprouts.
(a) SPIM setup with mounted sample (see S1 and S2 Figs for specifications). (b) SPIM-compatible in vitro model of angiogenesis showing endothelial cells invading and sprouting into a collagen gel with embedded fluorescent beads. (c) Live fluorescence cell imaging of the sprouts and beads captures cell-matrix mechanical interactions and allows mapping sprout-induced matrix deformations by (d) registering the fluorescent beads and sprout segmentation in time.

SPIM unit

The SPIM (or light sheet) imaging setup is based on a light sheet illuminating a thin (3 μm) section of the sample at the focal plane of the imaging objective (S1 and S2 Figs).

Light sheet generation of the SPIM microscope

In its simplest form, SPIM setups combine a single cylindrical lens together with an objective in a telescope configuration to generate a sheet/plane of light along the direction of the optical path. Here we decided to generate a virtual light sheet (S1 and S2 Figs) to yield a better-defined illuminated volume inside the sample. The virtual light sheet is formed by the beam waist of a focussed Gaussian beam that is being swept over a plane perpendicular to the imaging objective this way illuminating a planar volume inside the sample. If the sweeping frequency is synchronized and sufficiently high in comparison to the exposure time of the camera, the swept gaussian beam will appear as a plane of light. To deflect the beam angularly, we use a galvanometric mirror (62xxH series, Cambridge Technology, USA) with a suitable driver. A scan lens (f = 80 mm) converts the angular movement of the galvanometric mirror into a vertical, parallel translation. A tube lens (f = 200 mm) in combination with the illumination objective (Nikon Plan Flour 10x, NA = 0.3, working distance: 16 mm) acts as a telescope. Consequently, the objective back focal plane is effectively conjugated to the scanner. An angular beam tilt in the scanner plane therefore results in an angular beam tilt in the illumination objective back focal plane and is translated into a vertical translation behind the objective. The generated virtual Gaussian light sheet is adjusted to overlap with the focal plane of the imaging objective. In the vertical direction, this is done by adjusting an angular offset of the galvo mirror. In z-direction (perpendicular to the plane of the sheet), the adjustment is performed by changing the angle of a mirror in the conjugate plane of the illumination objective. The width of the Gaussian focus with NA = 0.3 is a nominal 820 nm with a length (confocal parameter) of just 4 μm. As light sheet imaging is a trade-off between field of view and light sheet thickness, we purposefully further reduced the numerical aperture with an iris to have a focus length of > 50 μm with a width of about 3 μm.

The waveform for driving the galvanometric mirror is customised and generated by a programme written in C#. The digital to analogue conversion is performed by a data acquisition card (NI USB-6341).

Two laser sources have been used throughout this study. For two-colour displacement imaging we employed a 488 nm Argon ion laser (Model 163-C1210, Spectra Physics, Mountain View, CA, USA). Note that the red fluorescent beads could also be sufficiently excited using this wavelength. For DAPI diffusion experiments we coupled a 375 nm UV diode laser into the system (Excelsior, 375 nm, Spectra Physics, Mountain View, CA, USA).

Detection subsystem

The imaging unit of the SPIM platform (S1 and S2 Figs) obeys the principles of a widefield fluorescence microscope. The sample is in first instance imaged by an objective lens (Nikon NIR APO, 40x, NA = 0.8, working distance: 3.5 mm). Our samples are completely immersed in water which demands an infinity corrected water dipping lens. Furthermore, for the live cell experiments we heat up the whole medium to 37 degrees. The contrast between cells and fiducial markers is generated by introducing two different fluophores, with green emission (LifeAct/GFP, 520 nm, F-actin of the cells) and red emission (605 nm, fluorescent beads). A tube lens (f = 200 mm) forms the primary image on two charge coupled device (CCD) chips of two separate cameras (Hamamatsu ImageEM CCD digital camera, resolution: 512×512 pixels, pixel pitch: 16 μm, frame rate: up to 33 frames per second (fps) in full frame mode, Hamamatsu Photonics K.K., Japan). Each of the cameras is synchronised to the sweeping frequency of the sheet generating signal to avoid random ghost images. A dichroic filter splits red and green emission light for dual colour imaging (600 nm, 25.2 x 35.6 mm, Techspec–Edmund optics). Residual excitation light is blocked by adequate filters for both emission channels. For the fluorescent proteins LifAct/GFP, we use a 500 nm long pass filter (Chroma, USA). For imaging the fluorescent fiducial marker beads, we use a 600 nm bandpass filter (FB600–40, Thorlabs). In order to overlay the two planes of both cameras, we used a sample containing beads only and removed the green bandpass filter (F1, Edmund #67–030, 520 nm FWHM:40 nm) and replaced it with a 500 nm long pass filter. The second camera was then overlaid with the first image by means of a micrometer stage, with an error limited to the optical resolution of the setup. The effective 20 nm broad spectral detection band for the fiducial markers provides strong signals and sufficient contrast. For DAPI we change the filter set of the green detection channel to a blue-shifted (440–20 nm) filter.

SPIM movement unit, sample containment and heating

We opted for sample movement, hence unlike on a confocal or wide field microscope, the sample is moved for acquiring three-dimensional image stacks. The sample is translated along the optical axis of the imaging unit (z-axis) and hence through the cross-section of the imaging focal and illumination plane. Since the z-movement is critical for reliable image acquisition we use a computer controlled piezo stage (Q-545.20, Physik Instrumente, Germany). Using this piezo and by carefully adapting its driver, we care for sample inertia. For X, Y and coarse Z movements, a compact 4 axis motor assembly USB-4D-Stage was used (Picard-Industries, USA). The objectives are mounted in a container of acrylic glass and sealed using vacuum grease. For live-cell experiments, cell medium is filled into the chamber and warmed up using a heating element. A proportional-integral-derivative (PID) controller was implemented on an Arduino Uno microcontroller board (S1 and S2 Figs). It gets feedback from a thermistor electronic element which is changing its resistance upon temperature variations. The Arduino output signal is amplified and controls the heating element. After a warming-up period and a short overshoot, the temperature regulation of the PID-controller is stable enough for live-cell experiments (S1 and S2 Figs).

SPIM acquisition and control software

We use Micro-Manager, a software package for the control and synchronisation of the SPIM platform parts. For this software package, a plugin was written by the developers of OpenSPIM. This plugin controls the cameras and sample stages (for xyz movement of the sample) to acquire three-dimensional stacks. The software was adapted for stable dual colour imaging. We also supplied an additional customised driver for the Physik Intrumente piezo stage. The waveform generation is not yet included in the Micro-Manager software but has its own graphical user interface. Similarly, the Arduino controller is read out simply by USB (RS232 protocol) and a terminal programme.

In vitro angiogenesis model preparation

Cell culture and transduction

Pooled wild-type human vein umbilical cord (HUVEC) cells were purchased from Angio-Proteomie (Boston, USA). The cells were subcultured at approximately 90% confluency or used for displacement microscopy experiments up to passage 6. The cells were maintained in complete endothelial growth medium (EGM-2 BulletKit, Lonza, Switzerland), which was exchanged 3 times per week. Before live cell fluorescence microscopy, 70% confluent wild type HUVEC cells were transduced with adenoviral LifeAct-GFP2 (Ibidi, Martinsried, Germany) at multiplicity of infection (MOI) of 10 and incubated for 16 to 24 h. For diffusion experiments, non-fluorescent wild-type cells were stained with DAPI (Invitrogen, Belgium) in situ.

Rat tail collagen type I hydrogel preparation

Collagen mixes with a total volume of 1 ml before casting hydrogels in the tube containers were prepared on ice in a 50 ml conical tube by mixing collagen type I (Merck Millipore, Darmstadt, Germany) with 150 μl sodium bicarbonate buffer (15.6 mg/mL) and the appropriate variable amount of EGM2 containing (5 μl of) fluorescent bead solution (2% solids, 0.2 μm diameter, carboxylated, ex/em 580/605, Invitrogen, Belgium) to a final collagen concentration as indicated in the diffusion results, between 0.5 and 3 mg/ml. To initiate the polymerization reaction, a pH switch of the acidic collagen mixture was induced by using 1 M sodium hydroxide (Sodium hydroxide, VWR Chemicals, USA), after which the collagen was supplemented with pro-angiogenic factor Sphingosine-1-Phosphate (S1P, Sigma-Aldrich, Belgium) to a final concentration of 1 μM. Using a Pasteur pipette, the collagen was then introduced into the container to fill the 15 mm long part of the tube (S3 Fig). The hydrogels were allowed to polymerize and equilibrate in an incubator with 5% CO2 at 37°C for a minimum of 30 minutes, before proceeding with the cell seeding.

In vitro angiogenesis model tuned for SPIM

The in vitro angiogenesis model (S3 Fig) was an endothelial invasion assay where HUVEC cells, starting from a densely seeded cell layer, invade and form (~50–250 μm long) sprouts overnight within a 1.5 mg/ml collagen hydrogel containing pro-angiogenic factor S1P [5357]. The physical sample setup was designed to meet the requirements for fast imaging with SPIM, which dictates that the sample has to be transparent from two sides. Cylindrical imaging chambers with 2 mm inner diameter were fabricated from laboratory grade tubing (Bola, S 1815–07). This fluoroplastic polymer (FEP) tube exhibits a refractive index which is similar to the one of water. Hence imaging aberrations are minimized when imaging through the plastic using a water dipping objective. The tube was sealed at the bottom employing two component glue and was sterilized by UV light (30 min) before hydrogel casting. An amount of about 10 μl of the collagen mixture was pipetted into the imaging tube. After polymerization as described previously, 60000 cells from a 70–90% confluent HUVEC culture in 10 μl of EGM2 were seeded through the opening and the tube was placed for 15 minutes in the incubator at 37°C, 5% CO2. This allowed the cells to sink and form a cellular layer on top of the collagen. Afterwards, abundant medium was added to the sample. The cells were then allowed to initiate invasion and form sprouts overnight before live imaging sprout dynamics.

Imaging

Imaging single sprouts

Samples were transferred to the SPIM setup and live imaged at 37°C in M199 medium supplemented with 125 mM HEPES to buffer the medium to a neutral pH in ambient air. The fluorescence signals of the LifeAct-GFP2-cells and the beads were recorded using the setup described above. Spectral filtering ensured that beads and cells were imaged separately and simultaneously. Large volume timelapses consist of 200 x 200 x 200 μm3 volumes (512 pixels * 512 pixels * 201 planes), which were acquired for approximately 4h with temporal resolutions of 1 min (Res = 1 min). For local (thin volume) fast scanning, volumes of 200 x 200 x 20 μm3 (512 pixels * 512 pixels * 21 planes) with temporal resolutions up to 4 seconds (Res = 4 s) were recorded. In both scan modes, data was thus captured with voxel sizes of 0.39 x 0.39 x 1 μm3. To obtain the relaxed end state required to compute absolute displacements, cell relaxation was chemically induced by the actin cytoskeleton inhibitor cytochalasin D, and the gradual relaxation of cells was monitored during up to 4 hours of relaxation (S1 Video and S4 Fig).

Imaging DAPI diffusion experiments

To estimate the diffusion time preceding the chemically induced relaxing of the sprouts, experiments with 4’,6-diamidino-2-phenylindole (DAPI) were performed. DAPI is a cell-membrane permeable molecule (MW = 350 g/mol) that interacts with double stranded DNA which results in the formation of a fluorescent complex [58,59]. DAPI was added before imaging at final concentrations between 10 to 30 μM depending on the experiment. Subsequently, the samples were imaged for up to 3 hours. Fluorescence intensities were measured in ImageJ and plotted in Matlab (The MathWorks Inc., Natick, MA USA) fitting a diffusion model. The average inflection point was then estimated based on the distribution of inflection points (S4 Fig).

Imaging acellular control collagen samples

Cell-free control experiments were performed by recording timelapses of independent samples of representative collagen gels containing fluorescent beads but no cells. Two series of timings were chosen to mimic the cellular experiments: (1) 3D z-stacks of > 20 micron were acquired every 10–15 seconds for 2–3 hours. (2) Similar stacks were scanned every 4 seconds for 20 minutes. Other imaging parameters were similar to those from the cellular experiments.

Computational analysis

Enhancement of bead images

The quality of the acquired z-stacks of the fluorescent beads has a great impact on the calculation of the cell-induced displacement fields. To boost blob-like structures and reduce noise, these images were processed by a difference of Gaussian filter and combined with a contrast stretching operation [51].

Cell mask calculation

To obtain a mask of the cellular structures, z—stacks containing sprouts were first enhanced by a two-step process consisting of a penalized least squares denoising [60] followed by a contrast stretching operation to smoothen and highlight the cell layer and sprouting cells from the background. Then, processed images were used to segment cellular structures by Otsu’s thresholding algorithm [61].

Sample drift compensation

Enhanced images of the fluorescent beads were used to compensate for the translational shifts in the experimental setup, including the drift of the microscope stage and/or global sample displacements caused by the activity of the cell layer present during the timelapse acquisitions. Phase-correlation-based rigid image registration was used to efficiently correct for the average global shift between images of fluorescent beads. Drift compensation is usually performed by rigidly registering (i.e. aligning) the bead images at every time point to the bead image of a pre-selected fixed time-point that acts as a global reference, which is commonly selected as the bead image after mechanical relaxation of cells. This strategy works fine for small hydrogel deformations showing a fast decay and thus, drift compensation algorithms benefit from having large regions with negligible bead displacements in the acquired field of view. However, this approach hinders the calculations under large deformations or when significant bead displacements are present within the whole field of view, as expected for images of multiple sprouts growing into a collagen gel. Instead, assuming that both experimental drifts and cell-induced hydrogel displacements will be smaller between consecutive time points than between each time point and the global reference, images of fluorescent beads on consecutive time points were rigidly registered to exploit the temporal coherency over long time lapse acquisitions. Finally, the position of rigidly registered images was expressed in a global frame of reference.

SPIM-based 4D displacement microscopy

Drift-compensated images of fluorescent beads were used to compute the local matrix displacements induced by the cell layer and the sprouts. The calculation of the 4D matrix displacements was formulated as a multiscale Free Form Deformation (FFD) non-rigid image registration process, which has been proven suitable for dealing with full-field deformations under large displacement and strain regimes[51]. Briefly, in FFD-based displacement field calculation, two bead images are compared, where one bead image is warped to match another fixed one (reference image). The algorithm overlays a regular grid over the reference image, of which the nodes are defined as the control points of multivariate B-splines curves. Then, the position of the control points is tuned iteratively during an optimization process, which warps the moving image until it matches the reference image. This yields as output a full displacement field, i.e. a vector value at each voxel of the registered images. Although the method warps the moving image to fit the reference image, both the underlying B-spline-based transformation model and the resulting displacement field are defined from the reference to moving image. To cope with different levels of matrix deformations while providing smooth displacement fields, the described FFD-based registration scheme was performed following a coarse-to-fine multiscale strategy. This implies that bead images were first registered using a relatively large grid size (i.e. few control points) and then, the output was used as the starting point for the next scale where the grid is refined, repeating the process iteratively for multiple refinement scales. Specifically, a three-level multiscale approach was used here with 180 x 180 x 44 voxels, 90 x 90 x 22 voxels, and 45 x 45 x 11 voxels grid size for the coarsest, intermediate and finest scales, respectively. During the optimization process that drives the registration by warping the moving image towards the reference image, a normalized correlation coefficient was used as the metric that iteratively calculates the similarity between both images. Furthermore, a stochastic gradient descent method with adaptive estimation of the step size was selected as optimizer due to its efficiency for registering very large image volumes: full-field displacement fields of 3D images containing ~512 x 512 x 200 voxels were obtained within 30 to 45 seconds in a standard Intel Core i7 PC.

Removal of engulfed beads from the calculation of the displacement fields

Beads engulfed by the cells move independently of the cell-matrix mechanical interactions. Thus, to remove these beads from further data analysis and minimize errors in the calculation of the matrix displacements, the binary mask of the segmented cellular structures (sprouts and cell layer) was supplied to the FFD-based displacement field calculation algorithm. This mask ensured that engulfed beads were not corrupting the matrix displacements in the surroundings to the cell-matrix interface.

Absolute versus incremental displacement fields

To obtain absolute displacement fields (S5B Fig), the images of fluorescent beads in the stressed hydrogel under the effects of cellular tractions were registered to the image of fluorescent beads in the relaxed hydrogel acquired after chemically inducing mechanical relaxation of the cells. Hence, the image of fluorescent beads at each time point showing a stressed hydrogel state acted as the moving image during the calculation of FFD-based displacements, and the one for the relaxed hydrogel as the (fixed) reference image (tref = tfix). Absolute displacements provide quantitative information on the magnitude, direction and global distribution of cell-matrix mechanical interactions and it is the standard approach in Traction Force Microscopy experiments. However, it might conceal the instantaneous cell-matrix interactions occurring between stressed hydrogel states along the timelapse. Therefore, to unravel the mechanical interactions among sprouts within a given time-lag τ, the incremental displacements were computed (S5B Fig), where the image of fluorescent beads at each time point ti was used as the moving image and the one at time point ti−τ as the reference image in FFD-based registration (tref = ti−τ). In this way for every timepoint incremental displacement values are obtained compared to a time-lag τ earlier. Note that, while the magnitude of the incremental displacements reflects the instantaneous local changes in the displacement field within a time-lag τ, their direction should not be interpreted as the push-pull behavior of the sprouts (S5C and S5D Fig). Similar to absolute displacements, incremental displacements are computed at the same temporal resolution Res used during the acquisition, with Res τ. Therefore, the temporal dynamics of the cell-matrix mechanical interactions can be analyzed at different time scales (specified by τ) with a temporal resolution Res as good as the one provided by the imaging system.

Full-field and local rms displacements

To quantify the cell-matrix mechanical interaction over time, the root-mean-square (rms) displacement was computed from the displacement field for each time point in a selected region of interest as:


with ‖Un2 being the L2-norm of the displacement vector Un at location n within a region of interest defined by a total of N spatial points.

Two types of regions have been considered. A rectangular box around the sprout(s) of interest that provides the rms displacement of the full-field, and a mask covering the displacement field within a close vicinity (~ 6 micron) of the sprout surface, which provides a local rms displacement value. The local mask was obtained by a binary dilation operation applied to the cell mask of the selected sprout.

The full-field and local rms displacements have been computed for both absolute and incremental displacement fields.

Shape morphometry

The dynamics of the cell morphology were assessed qualitatively by contour evolution plots and quantitatively by the length of individual protrusions and the normalized volume difference over time.

Contour evolution plots display the changes of the cell shape over time by overlapping the color-coded contours of the Z-projected cell mask at different time points.

To calculate the length of a selected protrusion, the morphological skeleton of its Z-projected cell mask was first extracted. Then, the protrusion length at each time point was defined as the length of the shortest path within the corresponding skeleton from its base to its protrusion tip.

Finally, the normalized volume difference at each time point was obtained as:


where Mtref and Mti are the segmented cell masks at specified reference time point tref and at current time point ti, respectively. Analogously to the calculation of absolute and incremental displacement fields, the reference time point can be selected to be a fixed (tref = tfix) or relative (tref = ti—τ) time point.

Implementation and visualization

The computational workflow was implemented using custom-made code written in Matlab, except for FFD-based displacement field estimation, which was computed using the open-source software package elastix [62]. The integration of the elastix -based FFD image registration on the main computational workflow was done as previously described [51,63]. The software program Paraview 5.2 (Kitware Inc., NY USA) was used for the 3D rendering of segmented sprouts and the visualization of 3D vector fields.

Validation

Stability of the experimental setup and noise level of the displacements

The SPIM-based displacement microscopy platform was extensively validated by control experiments with cell-free collagen samples. To confirm the stability of the experimental setup and its suitability for calculation of cell-induced matrix displacements, acellular control hydrogels with embedded fluorescent beads were prepared and imaged as previously described to check the existence of unwanted motion of beads. Acquired time-series of bead image volumes were first processed and compensated for possible stage drifts and then, analyzed by spatio-temporal image correlation spectroscopy (STICS) in 3D (S6 Fig). Briefly, the three-dimensional spatial correlation map between every bead image at a given time-lag Δt apart from each other was computed and averaged in time. Then, the obtained averaged correlation peak was compared with the one obtained for Δt = 0 (time averaged autocorrelation). Advecting beads would in general shift the averaged correlation peak over different values of Δt, while diffusing beads would broaden it. To quantify the peak location and width at each Δt, the mean and standard deviation of a fitted multidimensional Gaussian function was used. Finally, note that bead advection was a priori expected to be negligible for drift-compensated images of fluorescent beads.

Additionally, the stability of the bead locations was qualitatively assessed by building a kymographic representation, showing spatial bead distributions in time of randomly selected z-planes of the drift-compensated images of fluorescent beads (S7 Fig). If no bead motion occurred, the beads were stable over time, resulting in straight lines along the time axis of the kymograph. On the other hand, diffusing beads would show an oscillating line along the time axis, while advecting beads would show a drifted line along the time axis.

Finally, the local non-rigid displacements of the drift-compensated images of fluorescent beads were computed for the control experiments to detect the existence of spurious displacements and get the background noise level of the FFD-based image registration for the current experimental setup.

Correlations and sample sizes

Auto- and cross-correlations in S1 Table between the temporal evolution of rms displacements and morphological metrics were calculated in Matlab (The MathWorks Inc., Natick, MA USA). For assessing the diffusion rate of the system, measurements from five independent samples were considered. For the STICS analysis, based on the collagen control experiments, nine stacks gathered from three independent samples were analyzed, for each of the three collagen concentrations measured.

Results and discussion

While the assessment of overall sprout progression in the in vitro model of angiogenesis (that typically runs over 24 hours) does not require any high temporal resolutions (as sprout progression is typically expressed as growth/hour or growth/day), capturing the fast morphological changes of the growing and retracting protrusions (and the displacements of the surrounding matrix) requires a minute to second time scale resolution (see S2 Video). Applying our fast 4D displacement microscopy methodology to our in vitro model of sprouting angiogenesis allows monitoring displacement fields not only at sufficiently high temporal resolution (1 minute) but also for multiple sprouts simultaneously (Fig 2 and S3S5 Videos) captured within a 200 x 200 x 200 μm3 volume with voxel sizes of 0.39 x 0.39 x 1 μm3. For comparable spatial resolution and image quality, this is 20 times faster than with a conventional confocal microscope (measured here for a Fluoview FV1000, Olympus, Tokyo, Japan). The few other groups investigating displacements around in vitro sprouts have reported temporal resolutions of 20 minutes for approximately half of the stack size [36] or 15 seconds for only a single plane [35]. Absolute displacements (Fig 2A) unveil that adjacent sprouts invading from a dense layer of cells display qualitatively similar matrix displacement fields. The general pulling activity at the sprout protrusions is oriented towards the sprout base, although with variable magnitude (between 1 μm and 6 μm, see S3 Video) among sprouts. At the same time, the cell layer induces large displacements at the base of the sprouts directed towards the sprout tip. While absolute displacements (tref = tfix) inform us on the overall magnitude and direction of sprout-matrix mechanical interactions, key short-term dynamics are hidden by the combined activity of multiple sprouts. Short-term dynamics of cell-matrix mechanical interactions were unmasked by computing incremental displacements (tref = tiτ) (Methods and S5 Fig) for every timepoint ti with a time-lag τ of 5 minutes. Fast clustered fluctuations of displacements now vary spatially, temporally and from sprout to sprout (Fig 2B and S6 Video). Quantifying the incremental matrix deformations in close vicinity of each segmented sprout confirmed these observations. At timepoint t = 40 min., sprout (viii) is particularly displacing the matrix (3 μm) while others are stalled. (Fig 2F and S8 Fig). All together these results reveal that neighbouring sprouts invade collectively (within 24 hours) but in a non-synchronous, highly dynamic fashion; which is in line with current studies reporting that asynchronous selection of tip cells–regulated by asynchronous individual cell protein dynamics—is required for normal network formation [6466].

4D SPIM-based displacement microscopy enables matrix displacement quantification at high temporal resolution, simultaneously around multiple sprouts and in relation to protrusion dynamics.
Fig. 2. 4D SPIM-based displacement microscopy enables matrix displacement quantification at high temporal resolution, simultaneously around multiple sprouts and in relation to protrusion dynamics.
(a) Absolute (tref = tfix at t = 255 min) and (b) incremental (tref = ti−τ,with τ = 5 min) displacements of indicated timepoints (see S5 Fig). For visualization purposes the smallest displacements (< 2 μm for (a) and < 0.75 μm for (b)) are not shown. For the full range of displacement magnitudes, see S6 Video. (c) Segmented and colour-coded sprouts. One selected sprout is marked by a dashed box. (d) Contour evolution plots of sprout (ii) marked in (c) and projected in z-direction and time, showing the dynamics of its protrusions during the final 55 minutes of the chemically-induced relaxation. (e) Absolute displacements near the surface of sprout c-(ii). (f) Incremental rms displacements as a function of time in close vicinity of 4 selected sprouts from (c) (see S8 Fig for additional graphs). (g) Activity profile during relaxation of sprout c-(ii) as assessed by its length (blue line; see also white dashed line in (d)) and velocity (red line) as a function of time. Protrusion velocity is positive for protrusion extension and negative for retraction. Absolute rms displacements in close vicinity of the selected protrusion are plotted as well (black shaded line). Absolute displacements are obtained from a stress-free state at t = 255 min after addition of 4 μM CytoD at t = 135 min. Acquired volume of ~200 x 200 x 200 micron3 with 1 minute temporal resolution.

The computational methods in the SPIM-based displacement microscopy were further extended to capture changes of sprout morphological features (Methods) and relate them to matrix displacements here illustrated with selected sprout (ii) from Fig 2C for timepoints t = 200 min to t = 255 min (Fig 2D, 2E and 2G). Colour-coding the contour of the Z-projected cell mask over time (Fig 2D) allows the inspection of the sprout morphological evolution at a glance. Note that the high temporal resolution provided by SPIM together with the region-based local analysis allows capturing subtle morphological changes such as the time-varying retraction (length and velocity) of a single protrusion upon the addition of CytoD (at t = 135 min.) to obtain a stress-free state (at t = 255 min.). Matrix displacements around the same protrusion (Fig 2E) gradually decrease alongside with the protrusion retraction (Fig 2G).

While SPIM allows for mapping displacements of multiple adjacent sprouts with subcellular detail—displacements within voxel sizes of 0.39 x 0.39 x 1 μm3 were calculated around sprout protrusions (subcellular sized cellular extensions)—the system also provides the flexibility to scan more local, smaller volumes (~200 x 200 x 30 μm3) and with improved temporal resolution of every 10 (Fig 3A) or every 4 seconds (Fig 3B–3J). This fast-local scanning is particularly relevant for investigating fast cell morphological dynamics. From a region with adjacent sprouts (Fig 3A), the most active (i) and stalled sprouts (ii) were identified based on the incremental displacement fields (τ = 5 min) (see contour plot in Fig 3A and S6 Video), and defined as the sprouts with highest and lowest local average displacements maximum (in 5 minutes), which was up to 0.44 μm for the active one, and up to 0.25 μm for the stalled one. Note that co-existing active and stalled behaviour exposes the non-synchronicity of neighbouring sprouts described above. Sprout shape dynamics (contour evolution plots, Fig 3A) were quantified by monitoring the normalized volume difference over time (see Methods) and compared to the local incremental displacements, showing moderate to high correlations (see S1 Table) between sprout morphological changes and surrounding matrix deformations. Note that these correlations are measured for a single sprout, with multiple complex-shaped protrusions contributing to the displacement field.

Fast 4D displacement microscopy around dynamic sprout protrusions.
Fig. 3. Fast 4D displacement microscopy around dynamic sprout protrusions.
(a) 3D rendering of sprouts imaged with fast local SPIM scanning (every 10 seconds), with the magnitude of local incremental matrix displacements (tref = ti− 5 min) mapped on the surface. Contour evolution plot and activity profile of an active (i) and stalled sprout (ii) (see text for selection criteria) illustrate the positive correlation (S1 Table) between the normalized sprout volume difference and the incremental rms displacements around the selected sprouts. (b) 3D rendering of a retracting sprout imaged every 4 seconds with the local displacement magnitudes mapped on the surface. (c, d) Contour evolution plots showing shape dynamics for temporal resolutions of 1 minute and 4 seconds, respectively. (e-h) Local incremental displacement fields (tref = ti− 5 min) at indicated timepoints. (i, j) Activity profiles of protrusions (i) and (ii) of the sprout in (b). Protrusion (i) shows the highest inverse correlation of its length with the local incremental rms displacements (S1 Table).

For refined metrics, we hence fast-scan isolated sprouts with more sharply defined protrusions and perform local displacement microscopy analysis (Fig 3B–3H and S7 Video). Contour evolution plots illustrate the image acquisition improvement when switching from scanning with 1 minute resolution used for large volumes (Fig 2D and Fig 3C) to almost continuously capturing a single sprout by locally scanning every 4 seconds (Fig 3D). This is also illustrated in S8 Video. In absence of a CytoD induced stress-free state, we focus on incremental matrix displacements (τ = 5 min) to assess protrusion dynamics, while keeping a temporal resolution of 4 seconds (see Methods).

Quantification of the protrusion length and the velocity [67] of its change confirms a time-varying retraction process along with an increase of the local incremental displacements around the protrusion (Fig 3I and 3J). This (inverse) relation between protrusion length and local incremental matrix displacement is particularly strong for protrusion (i) (Fig 3B and 3I and S1 Table). In Fig 3E–3H, the sprout is extending its protrusions while the fluorescent beads move away from the protrusion, followed by the main retraction event (~20 micron between t = 8 min and t = 20 min, see Fig 3I and 3J) during which the fluorescent beads (S9 Video left) move towards the protrusion. Together with the displacement field (S9 Video right), these observations lead to the conclusion that the observed main retraction event is associated with an increase of contractile forces in the protrusion, which are pulling on the matrix. As such, it is different from retraction events observed after CytoD application (Fig 2D and 2G), which leads to protrusion shortening accompanied by a gradual reduction of both absolute and incremental displacements (see displacement pattern comparison in S9 Fig).

Conclusion and perspectives

In summary, the high temporal resolution of SPIM-based displacement microscopy allows correlating cell-matrix mechanical interaction to morphological changes at minute to second time scales. Extending displacement microscopy with incremental displacement calculations revealed subtle cell-matrix interactions that would have been overlooked by mere absolute displacement calculations. The possibility to simultaneously follow several large multicellular systems at subcellular level and with high temporal resolution opens new perspectives in mechanobiology, such as for monitoring cell-matrix interactions during angiogenesis [56,67] as well as other multicellular systems such as organoids [68]. Studying other fast subcellular sprouting processes, such as adhesion remodeling, would also be possible at the current imaging parameters, with average adhesion lifetimes reported to be of the order of ~15 min [69]. We believe that even faster processes—such as the protein turnover within apparently stable focal adhesions at the second to millisecond time scale [70,71]—could also be examined with SPIM by reducing the imaged volume. In conclusion, our results demonstrate the rich potential of SPIM-based dynamic displacement imaging for understanding sprout development in angiogenesis.

Supporting information

S1 Fig [a]
SPIM microscope.

S2 Fig [a]
SPIM setup and microscopy chamber used for the presented experiments.

S3 Fig [a]
Schematic of the 3D angiogenesis sample design, compatible with SPIM-based displacement microscopy.

S4 Fig [a]
Chemically-induced sprout relaxation.

S5 Fig [a]
Absolute versus incremental displacement calculations for different hypothetical scenarios of cellular pushing and pulling activity.

S6 Fig [a]
STICS analysis consecutively in three dimensions confirms sample and system stability, a prerequisite to recover reliable displacements in cellular samples.

S7 Fig [a]
Kymographic representation of the 3D STICS analysis as a visual control for fluorescent bead attachment to the collagen.

S8 Fig [top]
Local data processing around individual sprouts.

S9 Fig [a]
Comparison of protrusion retractions: sprout pulling versus sprout relaxing.

S1 Table [docx]
Correlations between displacements and morphology.

S1 Video [spim]
Timelapse of the chemically-induced sprout relaxation.

S2 Video [left]
Increasing the temporal resolution from 5 to 1 minute for imaging fast morphological sprout dynamics.

S3 Video [avi]
Timelapse of the full field absolute displacements around angiogenic sprouts.

S4 Video [avi]
Timelapse of the local absolute displacement field around angiogenic sprouts.

S5 Video [avi]
Timelapse of the magnitude of the absolute displacements closest to the sprout mapped on the segmented sprout surface.

S6 Video [left]

S7 Video [left]

S8 Video [left]
Increased temporal resolution from 1 minute to 4 seconds for imaging fast morphological protrusion dynamics.

S9 Video [left]
Timelapse of retracting sprout protrusions.


Zdroje

1. Dembo M, Oliver T. Imaging the Traction Stresses Exerted by Locomoting Cells with the Elastic Substratum Method. 1996;70: 2008–2022.

2. Dembo M, Wang Y. Stresses at the Cell-to-Substrate Interface during Locomotion of Fibroblasts. 1999;76: 2307–2316.

3. Butler J, Fredberg JJ. Traction fields, moments, and strain energy that cells exert on their surroundings. 2002. doi: 10.1152/ajpcell.00270.2001 11832345

4. del Alamo JC, Meili R, Alonso-latorre B, Rodrı J, Juan CA, Aliseda A, et al. Spatio-temporal analysis of eukaryotic cell motility by improved force cytometry. 2007.

5. Legant WR, Miller JS, Blakely BL, Cohen DM, Genin GM, Chen CS. Measurement of mechanical tractions exerted by cells in three-dimensional matrices. Nat Methods. 2010;7: 969–971. doi: 10.1038/nmeth.1531 21076420

6. Steinwachs J, Metzner C, Skodzek K, Lang N, Thievessen I, Mark C, et al. Three-dimensional force microscopy of cells in biopolymer networks. Nat Methods. 2016;13: 171–176. doi: 10.1038/nmeth.3685 26641311

7. Condor M, Steinwachs J, Mark C, Garcia-Aznar JM, Fabry B. Traction Force Microscopy in 3-Dimensional Extracellular Matrix. Curr Protoc Cell Biol. 2017;75: 1–20. doi: 10.1002/cpcb.24 28627753

8. Schwarz US, Soiné JRD. Traction force microscopy on soft elastic substrates: A guide to recent computational advances ☆. BBA—Mol Cell Res. 2015;1853: 3095–3104. doi: 10.1016/j.bbamcr.2015.05.028 26026889

9. Stout DA, Bar-Kochba E, Estrada JB, Toyjanova J, Kesari H, Reichner JS, et al. Mean deformation metrics for quantifying 3D cell–matrix interactions without requiring information about matrix material properties. Proc Natl Acad Sci. 2016;113: 2898–2903. doi: 10.1073/pnas.1510935113 26929377

10. Santi PA. Light sheet fluorescence microscopy: A review. J Histochem Cytochem. 2011;59: 129–138. doi: 10.1369/0022155410394857 21339178

11. Reynaud EG, Peychl J, Huisken J, Tomancak P. Guide to light-sheet microscopy for adventurous biologists. Nat Publ Gr. 2015;12: 30–34. doi: 10.1038/nmeth.3222 25549268

12. Gao L, Shao L, Chen B-C, Betzig E. 3D live fluorescence imaging of cellular dynamics using Bessel beam plane illumination microscopy. Nat Protoc. 2014;9: 1083–1101. doi: 10.1038/nprot.2014.087 24722406

13. Voie A, Burns D, Spelman F. Orthogonal-plane fluorescence optical sectioning: three- dimensional imaging of macroscopic biological specimens. 1993;170: 8371260. doi: 10.1111/j.1365-2818.1993.tb03346.x 8371260

14. Huisken J, Swoger J, Bene F Del, Wittbrodt J, Stelzer EHK. Live Embryos by Selective Plane Illumination Microscopy. 2004;305: 1007–1010.

15. Engelbrecht CJ, Stelzer EHK. Resolution enhancement in a light-sheet-based microscope (SPIM). 2006;31: 1477–1479.

16. Stelzer EHK. Light-sheet fluorescence microscopy for quantitative biology. Nat Methods. 2014;12: 23–26. doi: 10.1038/nmeth.3219 25549266

17. Thorn K, Kellogg D. A quick guide to light microscopy in cell biology. Mol Biol Cell. 2016;27. doi: 10.1091/mbc.E15-02-0088 26768859

18. Mohan K, Purnapatra SB, Mondal PP. Three dimensional fluorescence imaging using multiple light-sheet microscopy. PLoS One. 2014;9: 1–8. doi: 10.1371/journal.pone.0096551 24911061

19. Rasmi CK, Mohan K, Madhangi M, Rajan K, Nongthomba U, Mondal PP. Limited-view light sheet fluorescence microscopy for three dimensional volume imaging. 2017;im. doi: 10.1063/1.4938536

20. Fu Q, Martin BL, Matus DQ, Gao L. Imaging multicellular specimens with real-time optimized tiling light-sheet selective plane illumination microscopy. Nat Commun. 2016;7. doi: 10.1038/ncomms11088 27004937

21. Verveer PJ, Swoger J, Pampaloni F, Greger K, Marcello M, Stelzer EHK. High-resolution three-dimensional imaging of large specimens with light sheet-based microscopy. Nat Methods. 2007;4: 311–313. doi: 10.1038/nmeth1017 17339847

22. Magidson V, Khodjakov A. Circumventing photodamage in live-cell microscopy. Methods Cell Biol. 2013;114: 1–15. doi: 10.1016/B978-0-12-407761-4.00001-4

23. Udan RS, Piazza VG, Hsu CC -w. C, Hadjantonakis A a.-K, Dickinson ME. Quantitative imaging of cell dynamics in mouse embryos using light-sheet microscopy. Development. 2014;141: 4406–4414. doi: 10.1242/dev.111021 25344073

24. Patra B, Peng Y-S, Peng C-C, Liao W-H, Chen Y-A, Lin K-H, et al. Migration and vascular lumen formation of endothelial cells in cancer cell spheroids of various sizes. Biomicrofluidics. 2014;8: 052109. doi: 10.1063/1.4895568 25332736

25. Gualda EJ, Simão D, Pinto C, Alves PM, Brito C, Ries J, et al. Imaging of human differentiated 3D neural aggregates using light sheet fluorescence microscopy. Front Cell Neurosci. 2014;8: 1–10. doi: 10.3389/fncel.2014.00001

26. Rasmi CK, Madhangi M, Nongthomba U, Pratim Mondal P. Curtailed light sheet microscopy for rapid imaging of macroscopic biological specimens. Microsc Res Tech. 2016;79: 455–458. doi: 10.1002/jemt.22665 27059099

27. Duncan LH, Moyle MW, Shao L, Sengupta T, Ikegami R, Kumar A, et al. Isotropic Light-Sheet Microscopy and Automated Cell Lineage Analyses to Catalogue Caenorhabditis elegans Embryogenesis with Subcellular Resolution. J Vis Exp. 2019; 1–15. doi: 10.3791/59533 31233035

28. Funane T, Hou SS, Zoltowska KM, Van Veluw SJ, Berezovska O, Kumar ATN, et al. Selective plane illumination microscopy (SPIM) with time-domain fluorescence lifetime imaging microscopy (FLIM) for volumetric measurement of cleared mouse brain samples. Rev Sci Instrum. 2018;89. doi: 10.1063/1.5018846 29864842

29. Andilla J, Jorand R, Olarte OE, Dufour AC, Cazales M, Montagner YLE, et al. Imaging tissue-mimic with light sheet microscopy: A comparative guideline. Sci Rep. 2017;7: 1–14. doi: 10.1038/s41598-016-0028-x

30. Pampaloni F., Richa R., Ansari N., Stelzer E.H.K. (2015) Live Spheroid Formation Recorded with Light Sheet-Based Fluorescence Microscopy. In: Verveer P. (eds) Advanced Fluorescence Microscopy. Methods in Molecular Biology (Methods and Protocols), vol 1251. Humana Press, New York, NY

31. Schmitz A, Fischer SC, Mattheyer C, Pampaloni F, Stelzer EHK. Multiscale image analysis reveals structural heterogeneity of the cell microenvironment in homotypic spheroids. Sci Rep. 2017;7: 1–13. doi: 10.1038/s41598-016-0028-x

32. Galland R, Grenci G, Aravind A, Viasnoff V, Studer V, Sibarita J. 3D high- and super- resolution imaging using single-objective SPIM. 2015;12. doi: 10.1038/nmeth.3402 25961414

33. Planchon TA, Liang G, Milkie DE, Davidson MW, Galbraith JA, Galbraith CG, et al. Bessel beam plane illumination. Nat Methods. 2011;8: 417–423. doi: 10.1038/nmeth.1586 21378978

34. Gebhardt JCM. Single Molecule Imaging of Transcription Factor Binding to DNA in Live Mammalian Cells. 2013;10: 421–426. doi: 10.1038/nmeth.2411.Single

35. Kniazeva E, Weidling JW, Singh R, Botvinick EL, Digman MA, Gratton E, et al. Quantification of local matrix deformations and mechanical properties during capillary morphogenesis in 3D. Integr Biol. 2012;4: 431–439. doi: 10.1039/c2ib00120a 22281872

36. Du Y, Herath SCB, Wang Q, Wang D, Asada HH, Chen PCY. Three-Dimensional Characterization of Mechanical Interactions between Endothelial Cells and Extracellular Matrix during Angiogenic Sprouting. Sci Rep. 2016; 21362.

37. Fischer RS, Gardel M, Ma X, Adelstein R, Clare M. Myosin II mediates local cortical tension to guide endothelial cell branching morphogenesis and migration in 3D. Curr Biol. 2010;19: 260–265. doi: 10.1016/j.cub.2008.12.045.Myosin

38. Xue N, Bertulli C, Sadok A, Yan Y, Huang S. Dynamics of filopodium-like protrusion and endothelial cellular motility on one- dimensional extracellular matrix fibrils. 2014.

39. Santos-Oliveira P, Correia A, Rodrigues T, Ribeiro-Rodrigues TM, Matafome P, Rodríguez-Manzaneque JC, et al. The Force at the Tip—Modelling Tension and Proliferation in Sprouting Angiogenesis. PLoS Comput Biol. 2015;11: 1–20. doi: 10.1371/journal.pcbi.1004436 26248210

40. Carmeliet P, De Smet F, Loges S, Mazzone M. Branching morphogenesis and antiangiogenesis candidates: Tip cells lead the way. Nat Rev Clin Oncol. 2009;6: 315–326. doi: 10.1038/nrclinonc.2009.64 19483738

41. Geudens I, Gerhardt O. Coordinating cell behaviour during blood vessel formation. Development. 2011;138: 4569–4583. doi: 10.1242/dev.062323 21965610

42. Bordeleau F, Mason BN, Lollis EM, Mazzola M, Zanotelli MR, Somasegar S, et al. Matrix stiffening promotes a tumor vasculature phenotype. Proc Natl Acad Sci U S A. 2017;114: 492–497. doi: 10.1073/pnas.1613855114 28034921

43. Edgar LT, Underwood CJ, Guilkey JE, Hoying JB, Weiss JA. Extracellular matrix density regulates the rate of neovessel growth and branching in sprouting angiogenesis. PLoS One. 2014;9: 1–10. doi: 10.1371/journal.pone.0085178 24465500

44. Mason BN, Starchenko A, Williams RM, Bonassar LJ, Reinhart-King CA. Tuning three-dimensional collagen matrix stiffness independently of collagen concentration modulates endothelial cell behavior. Acta Biomater. 2013;9: 4635–4644. doi: 10.1016/j.actbio.2012.08.007 22902816

45. Miron-Mendoza M, Seemann J, Grinnell F. The differential regulation of cell motile activity through matrix stiffness and porosity in three dimensional collagen matrices. Biomaterials. 2010;31: 6425–6435. doi: 10.1016/j.biomaterials.2010.04.064 20537378

46. Shamloo A, Heilshorn SC. Matrix density mediates polarization and lumen formation of endothelial sprouts in VEGF gradients. Lab Chip. 2010;10: 3061–3068. doi: 10.1039/c005069e 20820484

47. Lee PF, Bai Y, Smith RL, Bayless KJ, Yeh AT. Angiogenic responses are enhanced in mechanically and microscopically characterized, microbial transglutaminase crosslinked collagen matrices with increased stiffness. Acta Biomater. 2013;9: 7178–7190. doi: 10.1016/j.actbio.2013.04.001 23571003

48. Hur SS, Zhao Y, Li Y-S, Botvinick E, Chien S, SUNGSIKH Ur, et al. Live Cells Exert 3-Dimensional Traction Forces on Their Substrata. Cell Mol Bioeng. 2009;2: 425–436. doi: 10.1007/s12195-009-0082-6 19779633

49. Eilken HM, Adams RH. Dynamics of endothelial cell behavior in sprouting angiogenesis. Curr Opin Cell Biol. 2010;22: 617–625. doi: 10.1016/j.ceb.2010.08.010 20817428

50. Potente M, Gerhardt H, Carmeliet P. Basic and therapeutic aspects of angiogenesis. Cell. 2011;146: 873–887. doi: 10.1016/j.cell.2011.08.039 21925313

51. Jorge-Peñas A, Bové H, Sanen K, Vaeyens MM, Steuwe C, Roeffaers M, et al. 3D full-field quantification of cell-induced large deformations in fibrillar biomaterials by combining non-rigid image registration with label-free second harmonic generation. Biomaterials. 2017;136: 86–97. doi: 10.1016/j.biomaterials.2017.05.015 28521203

52. Hall MS, Long R, Feng X, Huang YL, Hui CY, Wu M. Toward single cell traction microscopy within 3D collagen matrices. Exp Cell Res. 2013;319: 2396–2408. doi: 10.1016/j.yexcr.2013.06.009 23806281

53. Bayless KJ, Kwak H Il, Su SC. Investigating endothelial invasion and sprouting behavior in three-dimensional collagen matrices. Nat Protoc. 2009;4: 1888–1898. doi: 10.1038/nprot.2009.221 20010936

54. Nguyen D-HTD, Stapleton SC, Yang MT, Cha SS, Choi CK, Galie PA, et al. Biomimetic model to reconstitute angiogenic sprouting morphogenesis in vitro. Proc Natl Acad Sci U S A. 2013;110: 6712–7. doi: 10.1073/pnas.1221526110 23569284

55. Title B, Title S, Title C, Year C, Holdername C, Author C, et al. Mechanical and Chemical Signaling in Angiogenesis. 2013. doi: 10.1007/s10456-012-9300-2

56. Trappmann B, Baker BM, Polacheck WJ, Choi CK, Burdick JA, Chen CS. Matrix degradability controls multicellularity of 3D cell migration. Nat Commun. 2017;8: 1–8. doi: 10.1038/s41467-016-0009-6

57. Koh W, Stratman AN, Sacharidou A, Davis GE. Chapter 5 In Vitro Three Dimensional Collagen Matrix Models of Endothelial Lumen Formation During Vasculogenesis and Angiogenesis. Methods Enzymol. 2008;443: 83–101. doi: 10.1016/S0076-6879(08)02005-3 18772012

58. Barcellona ML, Cardiel G, Gratton E, Gardiel G, Gratton E. Time-resolved fluorescence of DAPI in solution and bound to polydeoxynucleotides. Biochem Biophys Res Commun. 1990;170: 270–280. doi: 10.1016/0006-291x(90)91270-3 2372293

59. Favilla R, Stecconi G, Cavatorta P, Sartor G, Mazzini A. The interaction of DAPI with phospholipid vesicles and micelles. Biophys Chem. 1993;46: 217–226. doi: 10.1016/0301-4622(93)80015-b 8343569

60. Garcia D. Robust smoothing of gridded data in one and higher dimensions with missing values. Comput Stat Data Anal. 2010;54: 1167–1178. doi: 10.1016/j.csda.2009.09.020 24795488

61. Otsu N. A Threshold Selection Method from Gray-Level Histograms. IEEE Trans Syst Man Cybern. 1979;9: 62–66. doi: 10.1109/TSMC.1979.4310076

62. Klein S, Staring M, Murphy K, Viergever MA, Pluim JPW. Elastix: A toolbox for intensity-based medical image registration. IEEE Trans Med Imaging. 2010;29: 196–205. doi: 10.1109/TMI.2009.2035616 19923044

63. Jorge-Peñas A, Izquierdo-Alvarez A, Aguilar-Cuenca R, Vicente-Manzanares M, Garcia-Aznar JM, Van Oosterwyck H, et al. Free form deformation-based image registration improves accuracy of traction force microscopy. PLoS One. 2015;10: 1–22. doi: 10.1371/journal.pone.0144184 26641883

64. Ubezio B, Blanco RA, Geudens I, Stanchi F, Mathivet T, Jones ML, et al. Synchronization of endothelial Dll4-Notch dynamics switch blood vessels from branching to expansion. Elife. 2016;5: 1–32. doi: 10.7554/eLife.12167 27074663

65. Bentley K, Chakravartula S. The temporal basis of angiogenesis. Philos Trans R Soc B Biol Sci. 2017;372. doi: 10.1098/rstb.2015.0522 28348255

66. Cruys B, Wong BW, Kuchnio A, Verdegem D, Cantelmo AR, Conradi LC, et al. Glycolytic regulation of cell rearrangement in angiogenesis. Nat Commun. 2016;7: 1–15. doi: 10.1038/ncomms12240 27436424

67. Fischer RS, Gardel M, Ma X, Adelstein RS, Waterman CM. Local Cortical Tension by Myosin II Guides 3D Endothelial Cell Branching. Curr Biol. 2009;19: 260–265. doi: 10.1016/j.cub.2008.12.045 19185493

68. Mammoto T, Mammoto A, Ingber DE. Mechanobiology and Developmental Control. Annu Rev Cell Dev Biol. 2013;29: 27–61. doi: 10.1146/annurev-cellbio-101512-122340 24099083

69. Zinn A, Goicoechea SM, Kreider-Letterman G, Maity D, Awadia S, Cedeno-Rosario L, et al. The small GTPase RhoG regulates microtubule-mediated focal adhesion disassembly. Sci Rep. 2019;9: 1–15. doi: 10.1038/s41598-018-37186-2

70. Stehbens S, Wittmann T. Analysis of focal adhesion turnover: A quantitative live cell imaging example. 2015; 335–346. doi: 10.1016/B978-0-12-420138-5.00018-5.Analysis

71. Wolfenson H, Lavelin I, Geiger B. Dynamic Regulation of the Structure and Functions of Integrin Adhesions. 2014;24. doi: 10.1016/j.devcel.2013.02.012.Dynamic


Článok vyšiel v časopise

PLOS One


2020 Číslo 1
Najčítanejšie tento týždeň
Najčítanejšie v tomto čísle
Kurzy

Zvýšte si kvalifikáciu online z pohodlia domova

Získaná hemofilie - Povědomí o nemoci a její diagnostika
nový kurz

Eozinofilní granulomatóza s polyangiitidou
Autori: doc. MUDr. Martina Doubková, Ph.D.

Všetky kurzy
Prihlásenie
Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.

Prihlásenie

Nemáte účet?  Registrujte sa

#ADS_BOTTOM_SCRIPTS#