Next Article in Journal
Long-Term Bilateral Neuromuscular Function and Knee Osteoarthritis after Anterior Cruciate Ligament Reconstruction
Previous Article in Journal
An Isolated CNN Architecture for Classification of Finger-Tapping Tasks Using Initial Dip Images: A Functional Near-Infrared Spectroscopy Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Mechanical and Cellular Effects of Rhythmic Vertical Vibrations on Adherent Cell Cultures

by
Dongho Kwak
1,2,*,
Thomas Combriat
2,3,
Alexander Refsum Jensenius
1 and
Petter Angell Olsen
2,4
1
RITMO Centre for Interdisciplinary Studies in Rhythm, Time and Motion, Department of Musicology, University of Oslo, 0371 Oslo, Norway
2
Hybrid Technology Hub, Centre for Organ on a Chip-Technology, Institute of Basic Medical Sciences, University of Oslo, 0372 Oslo, Norway
3
Department of Physics, Njord Center, University of Oslo, 0316 Oslo, Norway
4
Unit for Cell Signaling, Department of Immunology and Transfusion Medicine, Oslo University Hospital, 0372 Oslo, Norway
*
Author to whom correspondence should be addressed.
Submission received: 5 June 2023 / Revised: 23 June 2023 / Accepted: 29 June 2023 / Published: 6 July 2023
(This article belongs to the Section Cellular and Molecular Bioengineering)

Abstract

:
This paper presents an innovative experimental setup that employs the principles of audio technology to subject adherent cells to rhythmic vertical vibrations. We employ a novel approach that combines three-axis acceleration measurements and particle tracking velocimetry to evaluate the setup’s performance. This allows us to estimate crucial parameters such as root mean square acceleration, fluid flow patterns, and shear stress generated within the cell culture wells when subjected to various vibration types. The experimental conditions consisted of four vibrational modes: No Vibration, Continuous Vibration, Regular Pulse, and Variable Pulse. To evaluate the effects on cells, we utilized fluorescence microscopy and a customized feature extraction algorithm to analyze the F-actin filament structures. Our findings indicate a consistent trend across all vibrated cell cultures, revealing a reduction in size and altered orientation (2D angle) of the filaments. Furthermore, we observed cell accumulations in the G1 cell cycle phase in cells treated with Continuous Vibration and Regular Pulse. Our results demonstrate a negative correlation between the magnitude of mechanical stimuli and the size of F-actin filaments, as well as a positive correlation with the accumulations of cells in the G1 phase of the cell cycle. By unraveling these analyses, this study paves the way for future investigations and provides a compelling framework for comprehending the intricate cellular responses to rhythmic mechanical stimulation.

1. Introduction

Mechanical cues within the cellular microenvironment play a crucial role in influencing fundamental processes such as proliferation, migration, and differentiation [1,2,3]. Cells can perceive and convert mechanical stimulation into biochemical signals through mechanosensors, which include ion channels, cell adhesion sites, and the cytoskeleton [4].
In order to mimic a dynamic mechanical environment, numerous strategies for vibrating cell cultures have been developed [2]. Accordingly, vibration has been investigated across diverse cell culture models, such as fibroblast cells [1,5] bone cells [6,7,8], adipose stem cells [9,10], mesenchymal stem cells [11,12] and human induced pluripotent stem cells [13,14]. A common finding in such studies is the notable influence of vibration on the structural arrangement of F-actin filaments within the cytoskeleton.
Actin filaments are one of the components of the cytoskeleton, along with microtubules and intermediate filaments. Among these structures, actin filaments are the thinnest and have a crucial role in organizing the intracellular environment. They are involved in various critical mechanical cellular processes, including providing cellular adhesion, enabling mobility, determining cell shape, contributing to cellular strength, and maintaining cellular structure [15]. The dynamic nature of actin filaments enables them to respond rapidly to the surrounding microenvironment, much like a tensegrity structure, thereby imparting mechanical resilience to cells [16,17]. Additionally, the dynamic properties of actin filaments may play a crucial role in transmitting mechanical stress from the cell surface or membrane to the nucleus [18,19,20].
The cytoskeletal structures show temporal variability and complexity [21,22]. For instance, responses to external mechanical cues of actin filaments within a cell can fluctuate regularly every ∼24 h [23]. However, while most studies involving mechanical stimuli on cells are carried out using continuous stimuli over a specific amount of time, there have only been a handful of experimental studies exploring the temporal regularity of mechanical stimulation of cell cultures using various apparatus. Some of these studies include, for example, rhythmic stretching (a 12 h cycle of stretching at 1 Hz for three consecutive days) of cell cultures using a flexible silicon substrate [24]. Through this method, Rogers et al. [24] demonstrated clock genes (e.g., ARNTL, PER, and SOX) entrainment in adult stem cells, which is important because the entrainment can regulate circadian rhythms in the cells. Additionally, Vágó et al. [25] showed enhanced tissue homeostasis and histogenesis through rhythmic (1 h at the same time of the day for six days) uniaxial compression of chondrogenic cell cultures using biocompatible stainless plates pressing down from above the cells on a Petri dish. These results can have a valuable implication in tissue engineering since rhythmic (dynamic) mechanical stimuli can be used to optimize and improve the integration of tissue constructs within a host body [24]. Such rhythmic “ordered patterns in time” [26] studies could also be relevant in the dynamic cell culture system using a vibration system.
The potential impact of dynamic cell culture systems using vibration setups with controlled temporal variations remains largely unexplored. To address this gap, our study aimed to investigate the effects of rhythmic mechanical stimuli on cell cultures. We devised and carefully characterized an experimental setup that involved subjecting cell cultures to rhythmic vertical vibrations with varying pulse durations and intervals, allowing us to explore different stimulation parameters. In particular, we created three vibration conditions: Continuous Vibration (CV), Regular Pulse (RP), and Variable Pulse (VP). We categorize CV as non-rhythmic and RP and VP as rhythmic conditions, and No Vibration (NV) as a control condition. Departing from conventional non-rhythmic paradigms, we utilized principles from digital audio signal processing to generate rhythmic conditions by manipulating digital signals in various ways.
To measure the performance of the setup, we utilized several techniques. These included three-axis acceleration measurements to characterize the vibration regarding peak and Root Mean Square (RMS) acceleration. This enabled us to further estimate the mechanical movement in terms of RMS velocity and displacement of the vibrations.
In addition to vibration analysis, we measured fluid flow within the chamber of the cell culture plate while it was subjected to vibrations. Such fluid flow can result in shear stress at the cell level that can impact cellular mechanosensing mechanisms, including mechanosensitive ion channels, cell–extracellular matrix adhesions, cell–cell adhesions, and cytoskeleton organization [27]. Consequently, the activation of mechanotransduction pathways through shear stress can significantly influence various cellular behaviors such as alignment, gene expression, differentiation, and migration [27]. Recognizing the significance of fluid dynamics and shear stress and its impact on cellular behavior, we decided to employ particle tracking velocimetry (PTV) to analyze fluid flow. Among several non-contact measurement methods, such as particle image velocimetry (PIV) and laser Doppler velocimetry (LDV), we found that the PTV method is particularly well suited for our study involving a small confined space and microparticles. While PIV is commonly preferred for capturing larger flow fields at a fixed space and time, PTV allows us to track individual particles at any given spatial and temporal points [28]. This provides higher spatial resolution and can be advantageous for unsteady flow and temporal variation analysis, which is assumed under RP and/or VP conditions. The LDV technique can typically cover a range between a few millimeters to several meters per second in terms of velocity. This technique involves scattering of laser light and requires the particles to be sufficiently large for accurate measurement [29]. Therefore, it is not suited to measure microparticles. Consequently, this reinforced our decision to utilize the PTV method. We were able to estimate shear forces in close proximity to the fluid–solid interface where cells were attached, adding a novel insight and introducing an additional mechanical factor in vibration systems. This information is invaluable for comprehensively understanding the vibration apparatus.
To assess the impact of vibration on cellular function, we analyzed F-actin filament structure, focusing on changes in filament length, thickness, and angle. Additionally, we examined the cell cycle distribution as a secondary functional readout on the vibrated cells.
With these methodologies and measurements, we aimed to comprehensively evaluate the experimental setup and shed light on the potential functional implications on cellular responses to dynamic mechanical stimuli.

2. Materials and Methods

2.1. Cell Culture

HeLa cells (ATCC, CCL-2) were cultured in a growth medium (Dulbecco’s Modified Eagle Medium; DMEM) containing 10% fetal bovine serum and 1% penicillin–streptomycin in a cell culture flask (25 cm 2 ) in an incubator at 37 C with 5% CO 2 . Before the plating of cells in 8-well glass-bottom slides (µ-slide, ibidi), the wells were coated with poly-L-lysine (Santa Cruz Biotechnology, Dallas, TX, USA). After the addition of poly-L-lysine (0.1% (w/v)), the slides were incubated at 37 C for five minutes, washed with phosphate-buffered saline (PBS) twice, and dried at room temperature for three hours, followed by disinfection under ultraviolet (UV) light for 15 min. Subsequently, trypsinized cells were added at 100,000 cells per well (growth area per chamber: 1 cm 3 ) and were incubated for 24 h before treatment with vibration.

2.2. Experimental Design

2.2.1. Vertical Vibration Setup

For vibration treatment, the ibidi glass-bottom slide (Figure 1a) was taken out of the incubator and placed and secured on top of a round (Ø = 100 mm) baseplate constructed by laser cutting a 4 mm thick PMMA sheet (Figure 1b) that is coupled with the vertical vibration generator (U56001, 3B Scientific) through a stand rod. The rod was connected to the speaker inside the device, and the frequency response of the generator ranged from 0 to 20 kHz. The generator was controlled by a digital signal generator developed in the Max/MSP environment (a visual computer programming language for music and multimedia, Cycling ‘74), where various experimental conditions (e.g., frequency, pulse durations, intervals, and waveform) were controlled and manipulated (a standalone application can be found in Supplementary Code S1). Through digital-to-analog conversion (DAC), an analog audio signal produced by the PC is sent to an arbitrary amplifier and eventually to the vibration generator, as illustrated in Figure 1c.

2.2.2. Experimental Conditions

The cells were treated with vertical vibration using 50 Hz sine tone for 15 min under different conditions that varied in temporal regularity. The conditions were:
  • No Vibration (NV): control (static) condition;
  • Continuous Vibration (CV): continuous (uninterrupted) sound vibration;
  • Regular Pulse (RP): rigidly regular pulses with a fixed pulse duration and interval at 1000 ms;
  • Variable Pulse (VP): irregular pulses with pulse durations and intervals varying unpredictably between 1000 and 1500 ms
We categorized CV as non-rhythmic and RP and VP as rhythmic conditions. The VP condition was included to contrast as much as possible with the precise regularity observed in the CV and RP conditions. Three separate experiments were performed for each condition, using a new batch of cells each time.

2.2.3. The Generation of Numerical Sequences for the Regular (RP) and Variable Pulse (VP) Conditions

The signals used to generate RP and VP were controlled and manipulated in the Max/MSP environment. For RP (Figure 2), a fixed input argument at 1000 ms was used for the “metronome object” to control signals. Accordingly, the pulse duration and interval were set at 1000 ms. For VP, a sequence of pseudorandom numbers was used (Figure 2), which was generated through the “random object” in Max/MSP. A specific sequence of random numbers can be recalled and regenerated using a seed message into the “random object”. When the seed message is not specified, the “random object” uses the system clock as the seed message. The generated numbers alter between 0 and 1, less than any integer argument. To create a sequence with a lower limit of 1000 and an upper limit of 1500, 500 was used as an integer argument for the “random object”. Then, the “plus object” was used to add 1000 to the output numbers of the “random object”, and the output of the “plus object” was used as an input argument into the metronome object. Consequently, the pulse duration and interval varied between 1000 and 1500 ms.
The generated sequences used for VP are stored separately and can be accessed and reused as necessary (Raw Data S1). One of the sequences is visualized in Figure 3 and was tested for its randomness using autocorrelation analysis (Figure 4) and was calculated using the runs test and entropy (Table 1). The blue shaded area of the autocorrelation in Figure 4 shows the confidence interval ( α = 0.05), in which any value not significantly correlated to the previous value is shown. The runs test and entropy results are described and presented in Table 1.

2.3. Vibration Measurement

The vertical vibration setup was characterized using a 3-axis accelerometer (AX3, Axivity). An accelerometer is a type of sensor commonly used to capture various types of motion, including vibrations. The USB rechargeable accelerometer we employed was lightweight (16 g) and is widely used for research purposes. To ensure the vertical nature of the vibration, we opted to utilize a 3-axis accelerometer capable of measuring vibrations in three dimensions: X-axis, Y-axis, and Z-axis. This provides an advantage over traditional approaches that rely on single-axis sensors for measurement [30]. The accelerometer was secured on top of the cell culture plate, and the vibration was measured at a sampling rate of 3200 Hz.

2.4. Particle Tracking Velocimetry (PTV) Measurement

To further characterize the setup and to understand and estimate the fluid dynamics within the chambers of the cell culture plate, the fluid flow under different vibration conditions was characterized using the PTV setup as shown in Figure 5. It was possible to image the plane passing through the middle of the fluid in the glass cuvette that had the same 2D dimensions (10 × 10 mm) as the ibidi slides through the laser sheet created using an especially 3D-printed slit holder. We kept the volume of the fluid (DMEM cell media) the same (350 µL) as in the actual experimental setup. The camera (Grasshopper3 USB3 CCD, Teledyne FLIR) settings were:
  • Shutter: 5 ms;
  • Gain: 12 dB;
  • Framerate: 25 FPS.
The “auto” settings of the camera were disabled.
In fluid mechanics, shear stress, denoted by τ , is a measure of the force per unit area that acts parallel to a surface due to the fluid’s velocity gradient. It is defined as:
τ = μ u y
where μ , u, and y are dynamic viscosity, velocity, and height, respectively. When a fluid flows over a solid surface or boundary, such as at the bottom of the dish, it is commonly assumed that the fluid adheres to the surface without slipping, known as the no-slip condition [31]. This implies that the fluid velocity at the solid surface is assumed to be the same as the velocity of the solid surface. Assuming such a condition at the bottom of the dish, the shear stress experienced by the cells can be approximated as:
τ μ u ( y = y 0 ) y 0
where y 0 represents the first minimal height above the bottom surface where the fluid velocity can be reliably measured using the PTV technique.

2.5. Staining of F-Actin Filaments

Following treatment with vertical vibration, the cells were incubated for 15 min in the CO 2 incubator at 37 C. Afterward, the cells were rinsed with cold PBS once and incubated with cold 4% paraformaldehyde (PFA) at room temperature for 15 min. Following this, the cells were washed twice with PBS and stored at +4 C. For staining of F-actin, fixed cells were permeabilized by incubation with 0.1% Triton X-100 in PBS at room temperature for 15 min. Next, the cells were incubated with AlexaFluor-488 Phalloidin (ThermoFisher, diluted 1:50 in 4% BSA/PBS) for one hour at room temperature and then washed three times with PBS to remove any unbound dye. Stained cells were stored at +4 C until microscopy analysis.

2.6. Microscopy and Image Acquisition

The ZEISS Elyra PS1 system was used with standard filter sets and laser lines with a Plan-APOCHROMAT, 63x/1.4 NA oil objective using the structured illumination microscopy (SIM) mode of the system. SIM imaging was performed using five grid rotations (0.51 mm grid) for 15 Z-planes with a 0.110 µm spacing. For SIM image reconstruction, the ZEN black software (MicroImaging, Carl Zeiss) was used with the following “method” parameters:
  • Processing: manual;
  • Noise filter: −4;
  • SR frequency weighting: 1;
  • Baseline cut, sectioning: 100/83/83;
  • PSF: theoretical;
  • Output: SR-SIM.
For image analysis, SIM images displayed as maximum intensity projections rendered from all Z-planes were used.
Nine different locations across the entire growth area of a well were imaged. Each vibration treatment was independently repeated three times, resulting in a total of 27 images per condition for analysis. All the acquired images can be found in Supplementary Raw Data S2.

2.7. Image Analysis of F-Actin Filament Structures

The feature extraction of the acquired fluorescence images was performed using the Scikit and Scipy image packages in Python for further analysis. An “Actin Analyzer” algorithm was developed for this particular analysis and is available in Supplementary Code S2.
The algorithm depicted in Figure 6 was developed for this analysis to detect and extract the F-actin filament structural features. Filaments represented in the images were first segmented using a thresholding method on the output of a Laplacian filter (an edge-detecting method; step 2 of Figure 6) applied to the images. For the detection to be reliable even with images slightly out of focus, an adaptive threshold was used (step 3 of Figure 6), and pixels p ( x , y ) belonging to filaments were considered to follow the criteria:
Δ p ( x , y ) Δ p std x , y ( Δ p )
where Δ p ( x , y ) , Δ p , and std x , y ( Δ p ) denote the value of the laplacian at position ( x , y ) , the average of the laplacian over the whole image, and its standard deviation, respectively.
The found contours were then skeletonized (i.e., made one pixel thick), and extremely short contours were considered as noise and filtered out (steps 4 and 5 of Figure 6). Finally, slight gaps in the filaments were filled by connecting filaments whose ends were both at a distance shorter than a set value and of similar relative orientation (i.e., angle; step 6 of Figure 6). Additionally, the extreme outliers (over 100 µm) were filtered out in the post-data analysis.

2.8. Analysis of Cell Cycle Distribution

Cell cycle distribution was analyzed using the Click-iT Plus EdU Flow Cytometry Assay Kit from ThermoFisher, following the manufacturer’s instructions. Briefly, 24 h after the vibration treatment (in 8-well glass-bottom slides), cells were exposed to 1 mM of Click-iT-EdU for 2 h. After trypsinization, cells from two “8-wells” were collected in the same tube and fixed, permeabilized, and stained with Alexa 488 Click-it-EdU label reagent to detect newly synthesized DNA. Total DNA was stained by incubating cells in FxCycle PI/RNase staining solution (ThermoFisher, Waltham, MA, USA) for 30 min. The BD Accuri C6 Plus flow cytometer (Becton Dickinson, Franklin Lakes, NJ, USA) was used to quantify the percentage of cells in the G1, S, and G2/M cell cycle phases. Dot plots showing total DNA (PI stained) were analyzed using gatings shown in Supplementary Figure S1. Four parallels were analyzed for each treatment.

2.9. Statistical Analysis

A t test was applied to compare the means of two independent data sets with normal distribution. For non-parametric data sets, Kruskal–Wallis (non-parametric ANOVA) was applied to test the significance level between experimental groups. As a post hoc test, the Wilcoxon rank-sum test was applied to compare the medians of the data sets with skewed distribution. All tests were performed using SciPy.stats package in Python. A p value < 0.05 was considered statistically significant. All the plots were generated in Python.

3. Results

3.1. Characterization of Vertical Vibration in Terms of Acceleration Parameters

The vertical vibration setup was characterized using a three-axis accelerometer (see Section 2). The measured acceleration in three dimensions (X, Y, and Z) showed acceleration only in the Z (vertical) direction in Figure 7, confirming that the direction of the vibration is vertical or perpendicular to the surface the setup rests on. A visual comparison in Figure 8 presents the differences between the experimental conditions. The fast Fourier transform (FFT) analysis revealed a few odd harmonics present above the fundamental frequency in the measured acceleration (Figure 9). The reason for the odd harmonics and the fundamental frequency being slightly higher (∼51 Hz) than the input frequency (50 Hz) could have been due to the mechanical resonances between the Ibidi cell culture slide and its lid. The total amount of time of vibration and silences (pulse intervals), the highest acceleration peaks, root-mean-square (RMS) accelerations, crest factors (the ratio between the highest peaks and the RMS values), RMS velocity (m/s), and RMS displacement (m) are presented in Figure 10. Additionally, the distribution of measured acceleration between the different experimental conditions is illustrated in Figure 11. The summary of the data sets is presented in Table 2.

3.2. Fluid Flow Patterns and Shear Stress Estimation in the Wells of the Cell Culture Plate through the Particle Tracking Velocimetry (PTV) Method

The PTV tracks revealed two vortexes in the plane where the direction of the flow was divided into left and right, resulting in negative and positive velocity and shear stress values at the bottom of the dish, respectively (Figure 12). Due to the symmetry of the glass cuvette, the flow pattern in that region is expected to be co-planar. The shear stress at the boundary layer (close to the bottom surface of the glass cuvette) was calculated using the velocimetry data measured from the images and is presented in Figure 12. Although the shear stress in non-constant flow observed under RP and VP is only an approximation and is not as reliable as in constant flow under CV, it was possible to visualize patterns with different regularities under RP and VP.
Under CV, due to the direction of the fluid flow that resulted in distinct vortexes, the estimated shear stress was the highest in the middle of each vortex close to the bottom surface of the cuvette. The shear stress varied more regularly under RP than VP. Furthermore, we analyzed the motion of the particles using an analyzer developed to extract motion over time [32]. The motiongram is a great tool to evaluate movement between successive frames in time, as presented in Figure 13. According to these estimations, the cells were exposed to more regular mechanical stress under RP than VP.

3.3. Impact of Vertical Vibration on the F-Actin Filaments in HeLa Cells

In order to determine the impact of vertical vibration on the cells, the F-actin filaments in cultures of adherent HeLa cells were examined. After subjecting the cells to a 15 min vibration treatment followed by 1 -min incubation without vibration, the cells were fixed, and the F-actin filaments were stained. We then used our custom algorithm (see Section 2) to measure the angle, length, and thickness of the F-actin filaments.
We report median values for the filament lengths and thicknesses here since the data sets show a non-normal (asymmetric) distribution (see Section 2). A detailed statistical summary of the data sets is available in Supplementary Table S1. In general, the data show a pattern of reduced F-actin filament size in all vibrated cell cultures, as presented in Figure 14a. In terms of filament lengths, treatment of all three vibration conditions produced shorter median filament lengths (CV: 0.60 µm, p < 0.001; RP: 0.64 µm, p < 0.001; VP: 0.64 µm, p < 0.001) than the control (NV: 0.68 µm). When it comes to filament thicknesses, compared to the control (NV: 0.100 µm), CV (0.095 µm, p < 0.001) and RP (0.097 µm, p < 0.001) resulted in thinner median F-actin filaments, while VP (0.101 µm, p < 0.001) produced slightly thicker filaments (Figure 14a). The impact of the vibration patterns on filament angle was also noticeable in the data. The filament angle data set showed a normal distribution; thus, t tests were employed to compare the differences in the means. Compared to the control (NV: 6.61 degrees), the mean angle of the filaments of the cell cultures was altered under CV (−1.67 degrees, p < 0.001), RP (1.05 degrees, p < 0.05), and VP (−0.47 degrees, n.s.) (Figure 14b).

3.4. Effect of Vertical Vibration on Cell Cycle Distribution

As a follow-up functional readout, we next analyzed if vertical vibration had an effect on the cell cycle distribution of the cells. The cell cycle analysis was performed 24 h after vibrational treatment using the EdU incorporation assay. As seen in Figure 15, the control (NV) exhibited a cell cycle distribution consisting of 44.88% in G1 phase, 45.02% in S phase, and 5.95% in G2/M phase. In CV-treated cells, a significant increase in G1-phase cells (47.23%, p < 0.001) and a significant decrease in G2/M-phase cells (5.08%, p < 0.01) was observed. RP-treated cells likewise displayed a significant increase in G1-phase cells (45.88%, p < 0.05), while VP treatment did not result in significant changes in the cell cycle distribution (Figure 15).

3.5. Correlations in Relation to the Magnitude of the Mechanical Stimuli

According to our data, the filament size (length and thickness) and accumulation of cells in the G1 phase exhibited a strong correlation with the mechanical parameters, while the filament angle appeared to be less/not influenced. The F-actin filaments of cells subjected to vibration exhibited a reduction in both length and thickness compared to the control (NV) cells. This effect was more pronounced in cells subjected to the non-rhythmic condition (CV), which had a higher RMS acceleration (Figure 10b) and shear stress compared to the rhythmic conditions (RP and VP) (Figure 12d), indicating a negative correlation between the magnitude of mechanical stimuli and the size of F-actin filament structure (Figure 16). Additionally, our data indicate that vibrated cells displayed a trend of G1-phase accumulation, with CV having a greater effect on cell cycle distribution than RP and VP. As a result, the observed increase in the level of G1-phase cells was positively correlated with the magnitude of the mechanical stimuli (Figure 16).

4. Discussion

The aim of this study was to develop a methodological protocol where cell cultures can be vertically vibrated rhythmically with varied temporal regularity (e.g., Continuous Vibration (CV), Regular Pulse (RP), and Variable Pulse (VP)). This was achieved by using a vertical vibration system and concept borrowed from music technology (i.e., the generation of rhythmic sound vibrations through audio signal processing in Max/MSP). The experimental setup was characterized in terms of acceleration (g), displacement (m), velocity (m/s), motiongram (quantity of motion), and estimated shear stress (Pa) in the cell culture chambers.
In our study, the rhythmic vibration conditions (RP and VP) translated into more dynamic movement with a varying degree of temporal regularity (rigid pulse duration and interval at 1000 ms or varied pulse duration and interval between 1000 and 1500 ms), lower RMS acceleration (RP: 0.97 g, VP: 1.03 g), and lower shear stress accumulation range (RP: 1.12–2.13 Pa, VP: 0.37–1.07 Pa) than the non-rhythmic condition (CV; RMS acceleration: 1.68 g; shear stress range: 3.18–3.53 Pa). Consequently, RP and VP produced mechanical stress that varied dynamically over the surface of the cell culture plate. CV also resulted in varied mechanical stress over the surface; however, the stress was constant and stable in comparison (Figure 12 and Figure 13).
To quantify the effects of vibration on the F-actin structure, an automated batch image analysis and feature extraction algorithm was developed to analyze the fluorescent F-actin images (“Actin Analyzer”; Supplementary Code S2). Batch processing was used to conduct an automated analysis of 108 images, and multiple tests were employed to verify the accuracy of the data produced by the algorithm (Supplementary Figures S2–S4). The effect of the three vibration conditions on the structure of F-actin filaments was evaluated by quantifying filament length, thickness, and angle.
The analysis revealed significant alterations in the structure of F-actin filaments in the vibrated HeLa cells compared to the non-vibrated cells and variations in the extent of changes observed in the results from CV, RP, and VP treatments. The main differentiating factor among the CV, RP, and VP conditions is the temporal regularity of the vibrations, which consequently affects the overall magnitude of the mechanical parameters. The significant negative correlation observed between shear stress and F-actin filament length, thickness, and angle (Figure 16) implies that, as the total accumulated vibration time and magnitude of mechanical stress increased, the filaments of the vibrated cells underwent a greater reduction in size. Actin filaments are known to be highly dynamic structures, regulated by both cellular mechanisms and mechanical forces [33], and several studies have reported that mechanical stress can cause actin filaments to shorten [34,35]. In addition, actin filaments can form a denser network under compressive loads, resulting in a shorter mean filament length [36]. However, when cells are subjected to shear stress, actin filaments can also be stretched [36,37]. Further characterization of the mechanisms underlying our observed changes in F-actin structure was beyond the scope of this study.
Mechanical forces can affect the progression of the cell cycle, as evidenced by previous studies that have reported contradictory findings. While some studies have demonstrated that mechanical stress can result in G1 cell cycle arrest [38,39], others have reported increased rates of progression through the G1 and G2 phases [40]. In our study, we observed an increase in the number of cells in the G1 phase of the cell cycle 24 h after CV and RP treatment (2.4% and 1.0% increases, respectively) (Figure 15).
To summarize, the main objective of this work was to introduce an experimental setup for exposing adherent cells to rhythmic vertical vibrations. While our focus was not solely on characterizing biological responses and underlying mechanisms, we sought to lay a solid foundation for future investigations. Various techniques were employed to characterize and quantify the mechanical forces exerted on the cells by such vibrations. By analyzing changes in F-actin structure and cell cycle distribution, we have provided examples of biological readouts that can be utilized to measure the functional effects of different vibration regimes. The described experimental system and characterization methodology can serve as a model for future studies investigating the impact of rhythmic mechanical stimulation on cells. We encourage subsequent investigations to delve deeper into the multifaceted nature of this phenomenon by considering critical variables such as amplitude, duration, frequency, and pattern of vibration to better comprehend the effect of rhythmic mechanical stimulation on important cellular features. Such features include cell proliferation, morphology, migration, gene and protein expression, signaling pathways, metabolic activity, and more. By exploring these variables, a deeper understanding of the specific cellular responses and underlying mechanisms in the context of rhythmic mechanical stimulation can be attained. Through our research, we aspired to not only illuminate the vast potential in rhythmic mechanical stimulation but also to catalyze further explorations that will unlock more profound insights into the intriguing interplay between mechanical forces and cellular behavior. This contribution to the field promises to pave the way for transformative discoveries in understanding cellular responses to mechanical stimulation.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/bioengineering10070811/s1, Code S1: Audio signal generation and manipulation Max/MSP patch—the standalone application can be downloaded at https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.7747414 4 June 2023; Code S2: Actin Analyzer, an image processing and feature extraction algorithm developed in Python, can be downloaded at https://github.com/donghodk/actin-analyzer 4 June 2023; Figure S1: Cell proliferation analysis using EdU proliferation kit; Figure S2: Algorithm test using Gaussian blur filter of filament lengths and thicknesses; Figure S3: Algorithm test using Gaussian blur filter of filament angle; Figure S4: Algorithm test comparing several algorithm parameters with extracted features overlaid on the original microscopic image; Raw Data S1: Pseudonumber sequences can be found at https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.7923984 4 June 2023; Raw Data S2: Raw microscopic images can be found at https://0-doi-org.brum.beds.ac.uk/10.18710/ALOBQK 4 June 2023; Table S1: A summary of the feature extraction data sets.

Author Contributions

Conceptualization, D.K. and A.R.J.; methodology, D.K., T.C. and P.A.O.; software, D.K. and T.C.; validation, D.K. and P.A.O.; formal analysis, D.K., T.C. and P.A.O.; investigation, D.K. and P.A.O.; resources, A.R.J. and P.A.O.; data curation, D.K.; writing—original draft preparation, D.K.; writing—review and editing, D.K., T.C., A.R.J. and P.A.O.; visualization, D.K. and T.C.; supervision, A.R.J. and P.A.O.; project administration, D.K. All authors have read and agreed to the published version of the manuscript.

Funding

Financial support from UiO: Life Science and the Research Council of Norway through its Centers of Excellence (CoE) funding scheme, project numbers 262613 (ABINO), 262762 (RITMO CoE) and 262613 (HTH CoE). The funders had no role in the study design, data collection, analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data can be found in the supplementary files.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Mohammed, T.; Murphy, M.F.; Lilley, F.; Burton, D.R.; Bezombes, F. The effects of acoustic vibration on fibroblast cell migration. Mater. Sci. Eng. C 2016, 69, 1256–1262. [Google Scholar] [CrossRef]
  2. Yuan, H.; Xing, K.; Hsu, H.Y. Trinity of three-dimensional (3D) scaffold, vibration, and 3D printing on cell culture application: A systematic review and indicating future direction. Bioengineering 2018, 5, 57. [Google Scholar] [CrossRef] [Green Version]
  3. Wall, M.; Butler, D.; Haj, A.E.; Bodle, J.C.; Loboa, E.G.; Banes, A.J. Key developments that impacted the field of mechanobiology and mechanotransduction: Mechanobiology and signal transduction. J. Orthop. Res. 2017, 36, 605–619. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Matthews, B.D.; Overby, D.R.; Mannix, R.; Ingber, D.E. Cellular adaptation to mechanical stress: Role of integrins, Rho, cytoskeletal tension and mechanosensitive ion channels. J. Cell Sci. 2006, 119, 508–518. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Enomoto, U.; Imashiro, C.; Takemura, K. Collective cell migration of fibroblasts is affected by horizontal vibration of the cell culture dish. Eng. Life Sci. 2020, 20, 402–411. [Google Scholar] [CrossRef]
  6. Kulkarni, R.N.; Voglewede, P.A.; Liu, D. Mechanical vibration inhibits osteoclast formation by reducing DC-STAMP receptor expression in osteoclast precursor cells. Bone 2013, 57, 493–498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Sakamoto, M.; Fukunaga, T.; Sasaki, K.; Seiryu, M.; Yoshizawa, M.; Takeshita, N.; Takano-Yamamoto, T. Vibration enhances osteoclastogenesis by inducing RANKL expression via NF-kB signaling in osteocytes. Bone 2019, 123, 56–66. [Google Scholar] [CrossRef] [PubMed]
  8. Pravitharangul, A.; Suttapreyasri, S.; Leethanakul, C. Iliac and mandible osteoblasts exhibit varied responses to LMHF vibration: Sites affect bone cell responses to vibration. Cell Biol. Int. 2018, 42, 1349–1357. [Google Scholar] [CrossRef]
  9. Halonen, H.T.; Ihalainen, T.O.; Hyväri, L.; Miettinen, S.; Hyttinen, J.A.K. Cell adhesion and culture medium dependent changes in the high frequency mechanical vibration induced proliferation, osteogenesis, and intracellular organization of human adipose stem cells. J. Mech. Behav. Biomed. Mater. 2020, 101, 103419. [Google Scholar] [CrossRef]
  10. Tirkkonen, L.; Halonen, H.; Hyttinen, J.; Kuokkanen, H.; Sievänen, H.; Koivisto, A.M.; Mannerström, B.; Sándor, G.K.B.; Suuronen, R.; Miettinen, S.; et al. The effects of vibration loading on adipose stem cell number, viability and differentiation towards bone-forming cells. J. R. Soc. Interface 2011, 8, 1736–1747. [Google Scholar] [CrossRef] [Green Version]
  11. Cho, H.; Seo, Y.K.; Jeon, S.; Yoon, H.H.; Choi, Y.K.; Park, J.K. Neural differentiation of umbilical cord mesenchymal stem cells by sub-sonic vibration. Life Sci. 2012, 90, 591–599. [Google Scholar] [CrossRef] [PubMed]
  12. Kim, I.; Song, Y.; Lee, B.; Hwang, S. Human mesenchymal stromal cells are mechanosensitive to vibration stimuli. J. Dent. Res. 2012, 91, 1135–1140. [Google Scholar] [CrossRef] [PubMed]
  13. Kanie, K.; Sakai, T.; Imai, Y.; Yoshida, K.; Sugimoto, A.; Makino, H.; Kubo, H.; Kato, R. Effect of mechanical vibration stress in cell culture on human induced pluripotent stem cells. Regen. Ther. 2019, 12, 27–35. [Google Scholar] [CrossRef] [PubMed]
  14. Terao, Y.; Kurashina, Y.; Tohyama, S.; Fukuma, Y.; Fukuda, K.; Fujita, J.; Takemura, K. An effective detachment system for human induced pluripotent stem cells cultured on multilayered cultivation substrates using resonance vibrations. Sci. Rep. 2019, 9, 15655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Fletcher, D.A.; Mullins, R.D. Cell mechanics and the cytoskeleton. Nature 2010, 463, 485–492. [Google Scholar] [CrossRef] [Green Version]
  16. Ingber, D.E.; Wang, N.; Stamenović, D. Tensegrity, cellular biophysics, and the mechanics of living systems. Rep. Prog. Phys. 2014, 77, 046603. [Google Scholar] [CrossRef] [Green Version]
  17. Reilly, C.B.; Ingber, D.E. Multi-scale modeling reveals use of hierarchical tensegrity principles at the molecular, multi-molecular, and cellular levels. Extrem. Mech. Lett. 2018, 20, 21–28. [Google Scholar] [CrossRef]
  18. Burridge, K.; Guilluy, C. Focal adhesions, stress fibers and mechanical tension. Exp. Cell Res. 2016, 343, 14–20. [Google Scholar] [CrossRef] [Green Version]
  19. Haque, F.; Lloyd, D.J.; Smallwood, D.T.; Dent, C.L.; Shanahan, C.M.; Fry, A.M.; Trembath, R.C.; Shackleton, S. SUN1 interacts with Nuclear Lamin A and Cytoplasmic Nesprins to provide a physical connection between the nuclear lamina and the cytoskeleton. Mol. Cell. Biol. 2006, 26, 3738–3751. [Google Scholar] [CrossRef] [Green Version]
  20. Wang, N.; Tytell, J.D.; Ingber, D.E. Mechanotransduction at a distance: Mechanically coupling the extracellular matrix with the nucleus. Nat. Rev. Mol. Cell Biol. 2009, 10, 75–82. [Google Scholar] [CrossRef]
  21. Goldbeter, A.; Gérard, C.; Gonze, D.; Leloup, J.C.; Dupont, G. Systems biology of cellular rhythms. FEBS Lett. 2012, 586, 2955–2965. [Google Scholar] [CrossRef] [PubMed]
  22. Kruse, K.; Jülicher, F. Oscillations in cell biology. Curr. Opin. Cell Biol. 2005, 17, 20–26. [Google Scholar] [CrossRef] [PubMed]
  23. Hoyle, N.P.; Seinkmane, E.; Putker, M.; Feeney, K.A.; Krogager, T.P.; Chesham, J.E.; Bray, L.K.; Thomas, J.M.; Dunn, K.; Blaikley, J.; et al. Circadian actin dynamics drive rhythmic fibroblast mobilization during wound healing. Sci. Transl. Med. 2017, 9, eaal2774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Rogers, E.H.; Fawcett, S.A.; Pekovic-Vaughan, V.; Hunt, J.A. Comparing circadian dynamics in primary derived stem cells from different sources of human adult tissue. Stem Cells Int. 2017, 2017, 2057168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Vágó, J.; Katona, E.; Takács, R.; Dócs, K.; Hajdú, T.; Kovács, P.; Zákány, R.; van der Veen, D.R.; Matta, C. Cyclic uniaxial mechanical load enhances chondrogenesis through entraining the molecular circadian clock. J. Pineal Res. 2022, 73, e12827. [Google Scholar] [CrossRef]
  26. Kwak, D.; Olsen, P.A.; Danielsen, A.; Jensenius, A.R. A trio of biological rhythms and their relevance in rhythmic mechanical stimulation of cell cultures. Front. Psychol. 2022, 13, 867191. [Google Scholar] [CrossRef]
  27. De Belly, H.; Paluch, E.K.; Chalut, K.J. Interplay between mechanics and signalling in regulating cell fate. Nat. Rev. Mol. Cell Biol. 2022, 23, 465–480. [Google Scholar] [CrossRef]
  28. Jolley, M.J.; Russell, A.J.; Quinn, P.F.; Perks, M.T. Considerations When Applying Large-Scale PIV and PTV for Determining River Flow Velocity. Front. Water 2021, 3, 709269. [Google Scholar] [CrossRef]
  29. Scislo, L. Single-Point and Surface Quality Assessment Algorithm in Continuous Production with the Use of 3D Laser Doppler Scanning Vibrometry System. Sensors 2023, 23, 1263. [Google Scholar] [CrossRef]
  30. Liu, J.; Wang, P.; Tian, X. Vibration displacement measurement based on three axes accelerometer. In Proceedings of the 2017 Chinese Automation Congress (CAC), Jinan, China, 20–22 October 2017; pp. 2374–2378. [Google Scholar] [CrossRef]
  31. Day, M.A. The no-slip condition of fluid dynamics. Erkenntnis 1990, 33, 285–296. [Google Scholar] [CrossRef]
  32. Jensenius, A. Using motiongrams in the study of musical gestures. In Proceedings of the International Computer Music Conference, New Orlenas, LA, USA, 6–11 November 2006; pp. 499–502. [Google Scholar]
  33. Kueh, H.Y.; Brieher, W.M.; Mitchison, T.J. Dynamic stabilization of actin filaments. Proc. Natl. Acad. Sci. USA 2008, 105, 16531–16536. [Google Scholar] [CrossRef]
  34. Butt, T.; Mufti, T.; Humayun, A.; Rosenthal, P.B.; Khan, S.; Khan, S.; Molloy, J.E. Myosin motors drive long range alignment of actin filaments. J. Biol. Chem. 2010, 285, 4964–4974. [Google Scholar] [CrossRef] [Green Version]
  35. Gao, J.; He, L.; Zhou, L.; Jing, Y.; Wang, F.; Shi, Y.; Cai, M.; Sun, J.; Xu, H.; Jiang, J.; et al. Mechanical force regulation of YAP by F-actin and GPCR revealed by super-resolution imaging. Nanoscale 2020, 12, 2703–2714. [Google Scholar] [CrossRef] [PubMed]
  36. De La Cruz, E.M.; Gardel, M.L. Actin mechanics and fragmentation. J. Biol. Chem. 2015, 290, 17137–17144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Yoshigi, M.; Hoffman, L.M.; Jensen, C.C.; Yost, H.J.; Beckerle, M.C. Mechanical force mobilizes zyxin from focal adhesions to actin filaments and regulates cytoskeletal reinforcement. J. Cell Biol. 2005, 171, 209–215. [Google Scholar] [CrossRef] [Green Version]
  38. Taïeb, H.M.; Garske, D.S.; Contzen, J.; Gossen, M.; Bertinetti, L.; Robinson, T.; Cipitria, A. Osmotic pressure modulates single cell cycle dynamics inducing reversible growth arrest and reactivation of human metastatic cells. Sci. Rep. 2021, 11, 13455. [Google Scholar] [CrossRef] [PubMed]
  39. Gupta, V.K.; Chaudhuri, O. Mechanical regulation of cell-cycle progression and division. Trends Cell Biol. 2022, 32, 773–785. [Google Scholar] [CrossRef]
  40. Tarannum, N.; Singh, R.; Woolner, S. Sculpting an Embryo: The Interplay between Mechanical Force and Cell Division. J. Dev. Biol. 2022, 10, 37. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A schematic diagram of the vertical vibration setup and its components. (a) Side view of the thin 8-well glass-bottom slide. (b) Exploded view of the PMMA baseplate. The black paper is included for better visualization, and the paper is not expected to impact the vibration propagation. (c) The vertical vibration experimental setup, including a PC, amplifier, and vibration generator.
Figure 1. A schematic diagram of the vertical vibration setup and its components. (a) Side view of the thin 8-well glass-bottom slide. (b) Exploded view of the PMMA baseplate. The black paper is included for better visualization, and the paper is not expected to impact the vibration propagation. (c) The vertical vibration experimental setup, including a PC, amplifier, and vibration generator.
Bioengineering 10 00811 g001
Figure 2. A section from the Max/MSP patch illustrating how the numerical sequences were generated to manage the digital audio signals regularly or irregularly within the Max/MSP environment. The red arrowheads were added to indicate the message or signal flow. The metronome object and a fixed integer argument ( i.e., 1000) are used to vary the signals regularly for the Regular Pulse (RP) condition. The random object that generates a sequence of pseudorandom numbers from 0 and 1 less than a set argument is used as input arguments into the metronome object to vary the signals irregularly for the Variable Pulse (VP) condition.
Figure 2. A section from the Max/MSP patch illustrating how the numerical sequences were generated to manage the digital audio signals regularly or irregularly within the Max/MSP environment. The red arrowheads were added to indicate the message or signal flow. The metronome object and a fixed integer argument ( i.e., 1000) are used to vary the signals regularly for the Regular Pulse (RP) condition. The random object that generates a sequence of pseudorandom numbers from 0 and 1 less than a set argument is used as input arguments into the metronome object to vary the signals irregularly for the Variable Pulse (VP) condition.
Bioengineering 10 00811 g002
Figure 3. (a) Visual representations of a pseudorandom number sequence used for the Variable Pulse durations and intervals. (b) A density histogram analysis shows the distribution of the durations and intervals. The two colored lines, the kernel density estimations, illustrate an approximately uniform distribution of the series of numbers within the set range between 1000 and 1500.
Figure 3. (a) Visual representations of a pseudorandom number sequence used for the Variable Pulse durations and intervals. (b) A density histogram analysis shows the distribution of the durations and intervals. The two colored lines, the kernel density estimations, illustrate an approximately uniform distribution of the series of numbers within the set range between 1000 and 1500.
Bioengineering 10 00811 g003
Figure 4. The autocorrelation analysis of a pseudorandom number sequence used for the varied (a) pulse durations and (b) pulse intervals. The maximum lag is 360 in this data set. Any data point falling within the blue shaded area ( α = 0.05) indicates that it is not significantly correlated with the previous point. The analysis illustrates that the coefficient values of the autocorrelation are close to 0, and the consecutive numbers are mostly not significantly related in the sequence. Therefore, the sequence is most likely produced randomly.
Figure 4. The autocorrelation analysis of a pseudorandom number sequence used for the varied (a) pulse durations and (b) pulse intervals. The maximum lag is 360 in this data set. Any data point falling within the blue shaded area ( α = 0.05) indicates that it is not significantly correlated with the previous point. The analysis illustrates that the coefficient values of the autocorrelation are close to 0, and the consecutive numbers are mostly not significantly related in the sequence. Therefore, the sequence is most likely produced randomly.
Bioengineering 10 00811 g004
Figure 5. Particle Tracking Velocimetry (PTV) setup with a VWR glass cuvette (10 × 10 mm, the same dimensions as a well of the ibidi glass bottom 8-well slide) secured on the 3D-printed slit holder on the vertical vibration generator. The slit holder and glass cuvette were held down firmly on the baseplate with screws. The laser sheet in the cuvette was created using a CO 2 laser and the slit, which was 0.1 mm in width and 10 mm in height. The glass cuvette was filled with DMEM cell media (350 µL) and microparticles (Ø = 10 µm). A Grasshopper3 USB3 CCD camera was used to capture images. The images were taken using a stroboscopic technique by keeping the framerate of the camera at 25 FPS while the vibration generator produced a 50 Hz sine tone.
Figure 5. Particle Tracking Velocimetry (PTV) setup with a VWR glass cuvette (10 × 10 mm, the same dimensions as a well of the ibidi glass bottom 8-well slide) secured on the 3D-printed slit holder on the vertical vibration generator. The slit holder and glass cuvette were held down firmly on the baseplate with screws. The laser sheet in the cuvette was created using a CO 2 laser and the slit, which was 0.1 mm in width and 10 mm in height. The glass cuvette was filled with DMEM cell media (350 µL) and microparticles (Ø = 10 µm). A Grasshopper3 USB3 CCD camera was used to capture images. The images were taken using a stroboscopic technique by keeping the framerate of the camera at 25 FPS while the vibration generator produced a 50 Hz sine tone.
Bioengineering 10 00811 g005
Figure 6. Schematic diagram of the image feature extraction algorithm to analyze F-actin filament structures (filament angle, length, and thickness) using fluorescence microscopic images. The two shaded areas in image 4 indicate the zoomed-in positions shown in images 5 and 6. (1) Acquired fluorescence microscopic image used as input; (2)–(6) images filtered through specific steps, and extracted features (filament angle, length, and thickness) were exported for further data analysis.
Figure 6. Schematic diagram of the image feature extraction algorithm to analyze F-actin filament structures (filament angle, length, and thickness) using fluorescence microscopic images. The two shaded areas in image 4 indicate the zoomed-in positions shown in images 5 and 6. (1) Acquired fluorescence microscopic image used as input; (2)–(6) images filtered through specific steps, and extracted features (filament angle, length, and thickness) were exported for further data analysis.
Bioengineering 10 00811 g006
Figure 7. A comparison of the acceleration measurement of the vertical vibration setup between X, Y, and Z dimensions shows acceleration only in the Z (vertical) direction.
Figure 7. A comparison of the acceleration measurement of the vertical vibration setup between X, Y, and Z dimensions shows acceleration only in the Z (vertical) direction.
Bioengineering 10 00811 g007
Figure 8. A comparison of measured acceleration between the experimental conditions. Note the continuous (uninterrupted) acceleration measured under Continuous Vibration (CV), rigidly the same pulse duration and interval under Regular Pulse (RP), and varying pulse duration and intervals under Variable Pulse (VP).
Figure 8. A comparison of measured acceleration between the experimental conditions. Note the continuous (uninterrupted) acceleration measured under Continuous Vibration (CV), rigidly the same pulse duration and interval under Regular Pulse (RP), and varying pulse duration and intervals under Variable Pulse (VP).
Bioengineering 10 00811 g008
Figure 9. Fast Fourier Transform (FFT) analysis of the acceleration measurement of each experimental condition. (Top row) The full range of the FFT analysis shows a few minor peaks around the odd harmonics above the fundamental frequency. (Bottom row) The same FFT analysis zoomed in around 50 Hz to show the fundamental frequency (∼51 Hz) in the measured acceleration.
Figure 9. Fast Fourier Transform (FFT) analysis of the acceleration measurement of each experimental condition. (Top row) The full range of the FFT analysis shows a few minor peaks around the odd harmonics above the fundamental frequency. (Bottom row) The same FFT analysis zoomed in around 50 Hz to show the fundamental frequency (∼51 Hz) in the measured acceleration.
Bioengineering 10 00811 g009
Figure 10. Visual representation of the parameters between the different vibration conditions. (a) A comparison of the total vibration and silence times between the conditions. (b) A comparison of the highest acceleration peak and RMS acceleration. (c) The crest factors (the ratio between the peak and RMS) show how “peaky” the signal is. The higher the value, the more the peakiness. CV scores show a lesser crest factor than the other two conditions, meaning the acceleration was more stable, and there was less variation between the acceleration peaks and the RMS throughout the sampled period. This corresponds to the RMS acceleration value, where CV scored the highest RMS value. (d) Comparison of the RMS velocity and (e) RMS displacement between different conditions.
Figure 10. Visual representation of the parameters between the different vibration conditions. (a) A comparison of the total vibration and silence times between the conditions. (b) A comparison of the highest acceleration peak and RMS acceleration. (c) The crest factors (the ratio between the peak and RMS) show how “peaky” the signal is. The higher the value, the more the peakiness. CV scores show a lesser crest factor than the other two conditions, meaning the acceleration was more stable, and there was less variation between the acceleration peaks and the RMS throughout the sampled period. This corresponds to the RMS acceleration value, where CV scored the highest RMS value. (d) Comparison of the RMS velocity and (e) RMS displacement between different conditions.
Bioengineering 10 00811 g010
Figure 11. Distribution of the acceleration of the vertical vibration under each condition. Under RP and VP, the distribution is dominated by the lowest category (0-1) in (Top row) and around 0 g in (Bottom row), which represent the intervals (silences) between pulses. (Top row) Sub-divisions of absolute values of the measured acceleration. (Bottom row) Histograms of the measured acceleration.
Figure 11. Distribution of the acceleration of the vertical vibration under each condition. Under RP and VP, the distribution is dominated by the lowest category (0-1) in (Top row) and around 0 g in (Bottom row), which represent the intervals (silences) between pulses. (Top row) Sub-divisions of absolute values of the measured acceleration. (Bottom row) Histograms of the measured acceleration.
Bioengineering 10 00811 g011
Figure 12. PTV results. The maximum intensity projection images, velocity fields, and shear stress plot at the boundary layer close to the bottom surface of the cuvette are shown according to the experimental conditions (a) CV, (b) RP, and (c) VP. (d) Distribution of estimated shear stress in absolute values, which shows a smoothed out distribution and where the values are concentrated.
Figure 12. PTV results. The maximum intensity projection images, velocity fields, and shear stress plot at the boundary layer close to the bottom surface of the cuvette are shown according to the experimental conditions (a) CV, (b) RP, and (c) VP. (d) Distribution of estimated shear stress in absolute values, which shows a smoothed out distribution and where the values are concentrated.
Bioengineering 10 00811 g012
Figure 13. Motion analysis of the particles from the images obtained using the PTV setup (Figure 5). The motiongrams (in black background) and the quantity of motion (QoM) plotted over time (1000 frames) reveal the contrasting degrees of regularity in the motion of the particles between CV, RP, and VP conditions.
Figure 13. Motion analysis of the particles from the images obtained using the PTV setup (Figure 5). The motiongrams (in black background) and the quantity of motion (QoM) plotted over time (1000 frames) reveal the contrasting degrees of regularity in the motion of the particles between CV, RP, and VP conditions.
Bioengineering 10 00811 g013
Figure 14. Analysis of the extracted features from the microscopic images. The results are statistically compared between the experimental conditions (NV, CV, RP, and VP). Red stars indicate the test that compared one of the vibration conditions to NV (control), gray stars indicate the test that compared the in-between vibration conditions, and ns indicates not significant. (a) The boxplot analysis (top) shows the distribution of the data. The box represents 50% of the data points, the green dotted line in the box represents the mean, and the orange line represents the median. The bar plots (bottom) show the differences in median filament lengths, thicknesses, and mean filament angles. (b) Filament angle distribution is presented with a density histogram and polar plots. Note that the angle does not show directions except for the 2D orientation within the collected microscopic images. The red dotted line in the polar plots indicates the mean value. * p < 0.05, *** p < 0.001.
Figure 14. Analysis of the extracted features from the microscopic images. The results are statistically compared between the experimental conditions (NV, CV, RP, and VP). Red stars indicate the test that compared one of the vibration conditions to NV (control), gray stars indicate the test that compared the in-between vibration conditions, and ns indicates not significant. (a) The boxplot analysis (top) shows the distribution of the data. The box represents 50% of the data points, the green dotted line in the box represents the mean, and the orange line represents the median. The bar plots (bottom) show the differences in median filament lengths, thicknesses, and mean filament angles. (b) Filament angle distribution is presented with a density histogram and polar plots. Note that the angle does not show directions except for the 2D orientation within the collected microscopic images. The red dotted line in the polar plots indicates the mean value. * p < 0.05, *** p < 0.001.
Bioengineering 10 00811 g014
Figure 15. Cell cycle distribution in control (NV) and vibrated (CV, RP, VP) HeLa cells analyzed 24 h after vibration treatment. Red stars indicate a significant difference relative to the control (NV) cells. Flow cytometry dot plots showing the gating strategy are included in Supplementary Figure S1. * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 15. Cell cycle distribution in control (NV) and vibrated (CV, RP, VP) HeLa cells analyzed 24 h after vibration treatment. Red stars indicate a significant difference relative to the control (NV) cells. Flow cytometry dot plots showing the gating strategy are included in Supplementary Figure S1. * p < 0.05, ** p < 0.01, *** p < 0.001.
Bioengineering 10 00811 g015
Figure 16. Correlation between the extracted features from the microscopic images (median filament length and thickness and mean filament angle) and characterized mechanical parameters (RMS acceleration and estimated shear stress). r is the Pearson correlation coefficient. Markers: NV—star, CV—circle, RP—square, VP—triangle.
Figure 16. Correlation between the extracted features from the microscopic images (median filament length and thickness and mean filament angle) and characterized mechanical parameters (RMS acceleration and estimated shear stress). r is the Pearson correlation coefficient. Markers: NV—star, CV—circle, RP—square, VP—triangle.
Bioengineering 10 00811 g016
Table 1. Randomness tests of the pulse durations and intervals used to generate the Variable Pulse (VP) condition. The total amount of numbers (n) was 360 for each sequence. The p value of the runs test indicates that if p < 0.05, then the sequence was not produced randomly. The entropy values are computed using Python’s Scipy.stats.entropy function. The closer the entropy value is to the maximum value, the higher the randomness of the sequence of numbers.
Table 1. Randomness tests of the pulse durations and intervals used to generate the Variable Pulse (VP) condition. The total amount of numbers (n) was 360 for each sequence. The p value of the runs test indicates that if p < 0.05, then the sequence was not produced randomly. The entropy values are computed using Python’s Scipy.stats.entropy function. The closer the entropy value is to the maximum value, the higher the randomness of the sequence of numbers.
nRuns Test (p Value)Entropy Max ValueEntropy Value
Pulse duration3600.45.895.27
Pulse interval3600.1 5.27
Table 2. A comparison of parameters between the different vibration conditions.
Table 2. A comparison of parameters between the different vibration conditions.
Data SetCVRPVP
Input frequency (Hz)505050
Measured fundamental frequency (Hz)∼51∼51∼51
Total vibration time (s)900450448.84
Total silence (intervals) time (s)0450451.16
Acceleration peak (g)3.123.923.96
Acceleration RMS (g)1.690.971.03
Crest factor1.864.033.84
Velocity RMS (m/s)0.0490.0290.031
Displacement RMS (m)0.000150.00010.0001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kwak, D.; Combriat, T.; Jensenius, A.R.; Olsen, P.A. Characterization of Mechanical and Cellular Effects of Rhythmic Vertical Vibrations on Adherent Cell Cultures. Bioengineering 2023, 10, 811. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering10070811

AMA Style

Kwak D, Combriat T, Jensenius AR, Olsen PA. Characterization of Mechanical and Cellular Effects of Rhythmic Vertical Vibrations on Adherent Cell Cultures. Bioengineering. 2023; 10(7):811. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering10070811

Chicago/Turabian Style

Kwak, Dongho, Thomas Combriat, Alexander Refsum Jensenius, and Petter Angell Olsen. 2023. "Characterization of Mechanical and Cellular Effects of Rhythmic Vertical Vibrations on Adherent Cell Cultures" Bioengineering 10, no. 7: 811. https://0-doi-org.brum.beds.ac.uk/10.3390/bioengineering10070811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop