Analysis of Statoliths Displacement in Chara Rhizoids for Validating the Microgravity-Simulation Quality of Clinorotation Modes

In single-celled rhizoids of the green algae Chara , positively gravitropic growth is governed by statoliths kept in a dynamically stable position 10–25 μ m above the cell tip by a complex interaction of gravity and actomyosin forces. Any deviation of the tube-like cells from the tip-downward orientation causes statoliths to sediment onto the gravisensitive subapical cell flank which initiates a gravitropic curvature response. Microgravity experiments have shown that abolishing the net tip-directed gravity force results in an actomyosin-mediated axial displacement of statoliths away from the cell tip. The present study was performed to critically assess the quality of microgravity simulation provided by different operational modes of a Random Positioning Machine (RPM) running with one axis (2D mode) or two axes (3D mode) and different rotational speeds (2D), speed ranges and directions (3D). The effects of 2D and 3D rotation were compared with data from experiments in real microgravity conditions (MAXUS sounding rocket missions). Rotational speeds in the range of 60–85 rpm in 2D and 3D modes resulted in a similar kinetics of statolith displacement as compared to real microgravity data, while slower clinorotation (2–11 rpm) caused a reduced axial displacement and a more dispersed arrangement of statoliths closer to the cell tip. Increasing the complexity of rotation by adding a second rotation axis in case of 3D clinorotation did not increase the quality of microgravity simulation, however, increased side effects such as the level of vibrations resulting in a more dispersed arrangement of statoliths. In conclusion, fast 2D clinorotation provides the most appropriate microgravity simulation for investigating the graviperception mechanism in Chara rhizoids, whereas slower clinorotation speeds and rotating samples around two axes do not improve the quality of microgravity simulation.


Introduction
The International Space Station (ISS), satellites, sounding rockets, parabolic plane flights and drop towers are research facilities that provide years, month, minutes or only seconds of a real-microgravity (μg) environment that has been intensively used in the last decades for gravitational biology research. However, experiments in real μg are scarce, costly and time consuming. Therefore, different μg simulation facilities like clinostats, Random Positioning Machine (RPM), rotating wall vessels and magnetic levitation are frequently used as attractive alternatives for complementing and preparing real-microgravity experiments in the near-Earth orbit. Furthermore, they enable with some limitations, such as sample diameter, biological studies of the impact of altered gravity conditions on Earth (Herranz et al. 2013;Brungs et al. 2016).
Clinostat capabilities include operation modes with one axis (2D clinostats) and two axes (3D clinostats). The classical 2D clinostat usually rotates constantly and slowly at 1-10 revolutions per minute (rpm), whereas fast rotating clinostats (Briegleb 1992;Herranz et al. 2013) run at a selected speed between 60-90 rpm. Due to the physical constraints only small samples placed precisely along the rotation axis will experience acceleration forces which are thought to be below the threshold of the biological system. Specimen larger than a few millimeters or samples not placed very close to the center of the rotation axis will experience centrifugal forces in the periphery.
3D-clinostats and RPMs operate based on the same principle with two independently rotating frames mounted in a gimbal manner allowing complex movements to average the influence of the gravity vector (Hoson et al. 1997;van Loon 2007;Borst and van Loon 2009). In 3D real random mode, the RPM operates with a range of velocities and also the direction of rotation is changed in a real random manner determined by dedicated algorithms. The movement of the samples is randomly generated to avoid a repetitive pattern to which the test samples may adapt to.
The various modes and angular rotation velocities create specific sets of forces and respective side effects like stress and shear forces with different impacts as well as positive and negative accelerations on cells or biological organisms. Due to the fact that reports on clinostat and random positioning experimentation frequently do not adequately state the operational modes and precise physical parameters of the treatment (Hammond and Allen 2011), a comparison of experiment results, calculation of forces and determination of threshold values is in many cases very difficult hampering generalization of conclusions (Herranz et al. 2013).
A diversity of model organisms such as Arabidopsis thaliana seedlings and plant cell cultures, characean rhizoids, human immune cells and protists have been investigated by several European science teams on various microgravity simulation facilities and the results were compared with real μg data (Herranz et al. 2013;Brungs et al. 2016). It became evident that a thoroughly performed assessment of side effects is required to determine the most suitable simulation facility for each biological object chosen. The quality of a simulation method is highly dependent on parameters like organism size, gravity perception time, threshold values and type and dense of culture media (Briegleb 1992;Hoson et al. 1997;Klaus et al. 1998;Hammond and Allen 2011).
In the past decades, 3D clinostats and the RPM operated in the 3D real random mode have been favored by investigators assuming that two rotation axes provide options for more complex ways to average the sensing of the gravity vector, thus, simulating μg more perfectly. To critically assess this assumption, the tip-growing unicellular rhizoids of the green algae Chara were chosen to demonstrate and analyze in detail the effects of the different modes and different rotational velocities on the dynamics of the actin-dependent basipetal movements of statoliths after abolishing the gravity force.
The Chara rhizoid is one of the best studied cell model for research on plant gravity sensing and the processes of gravity-oriented growth (Kiss 1997;Braun and Limbach 2006;Häder et al. 2017). In tip-downward growing rhizoids, statoliths are kept in a dynamically stable position 10-25 μm above the cell tip by two counteracting forces. The gravitational force (1 g) pulling the statoliths toward the cell tip is exactly counteracted by cytoskeletal forces generated by the complexly organized actomyosin system (Braun et al. 2002). Any deviation from the vertical growth direction, however, causes statoliths to sediment onto a specific gravisensitive area of the subapical cell flank, thereby initiating the gravitropic curvature response. Abolishing the gravity force results in an axial displacement of statoliths away from the cell tip (Volkmann et al. 1991;Buchen et al. 1993;Braun et al. 2002). Data on the dynamics of this movement of statoliths is available from various real μg conditions like parabolic plane flights, sounding rocket missions (TEXUS, MAXUS) and several space shuttle missions (Bunchen et al. 1991(Bunchen et al. , 1993Volkmann et al. 1991;Braun et al. 2002;Limbach et al. 2005). Chara rhizoids which have already been studied on 2Dand 3D clinostats (Cai et al. 1997;Hoson et al. 1997;Braun et al. 2002) showed similar effects like in real microgravity.
The present study aims to carefully evaluate the microgravity simulation quality of the different RPM modes by analyzing the basipetal statoliths transport in rhizoids particularly focusing on the relevance of a second rotation axes on the gravity-sensitive statoliths-cytoskeleton based system in this well-understood Chara rhizoid model system.

Rhizoid Culture
Chara globularis Thuill. was collected from a pond in the botanical garden in Bonn. For rhizoid cultivation, thallus segments containing two nodes and one internode were cut and the side branches were removed to induce rhizoid development. Thallus segments were embedded in 1.2% agar (Sigma-Aldrich, Missouri, USA) on microscopic slides, covered with long cover glasses and stored in staining jars under a 60 W tungsten light bulb at a photon fluence rate of approx. 60 μmol * m −2 * s −1 . After 5 days, the rhizoids developed and were used at a length of 5-10 mm for clinorotation experiments and for the MAXUS microgravity experiments.

Random Positioning Machine (RPM)
The "Desktop RPM" from Dutch Space (Leiden, Netherlands) was used to subject characean rhizoids to simulated microgravity conditions (Fig. 1). Its two frames are positioned perpendicular to each other and can be rotated separately. In 2D mode, the software-supported motorcontrolling runs only one of the two frames constantly around one horizontal axis (Fig. 1a). In 3D real random mode, the sample mounted on the inner frame is rotated around two axes with randomly changing rotational speed and directions. In this mode only the speed range is specified by the operator. For the present study, the RPM was operated in the 2D clinorotation mode with 3 rpm and 60 rpm as well as in the 3D real random mode with 2-11 rpm, 30-60 rpm and 60-85 rpm. In 2D clinorotation mode either the inner or the outer frame rotates constantly around one horizontal axis (b). In 3D real random mode both frames rotate with random speeds and directions (c) Chara rhizoids were mounted on the stage of one of the two platforms in such way that the cell tip was precisely aligned to the rotations axis in case of the 2D mode. In the 3D mode, the rhizoid tip was precisely positioned in the center of the two axes. All experiments were performed at 22 For videorecording, a CMOS camera module (RS-OV7949-1818, Conrad Electronic, Hirschau Germany), a 32x objective (Zeiss LD ACHROPLAN Ph2, 32×\0.40; ∞/1.5, Carl Zeiss AG, Oberkochen, Germany) and a white LED (75 • , 14000 mcd, Conrad Electronic, Hirschau, Germany) for illumination were integrated into optomechanical system components (Microbench system, Qioptiq, Paris, France).
Statolith tracking was done by labeling the lower and upper border of the statoliths complex in a sequence of still images prepared from the AVI-video sequences. The distances of the lower and upper border of the statoliths complex to the cell tip was measured with ImageJ (version1.44p, Wayne Rasband, National Institutes of Health, Bethesda, USA).

MAXUS-Sounding Rocket Flights
Characean rhizoids were launched on ESA-funded MAXUS sounding rockets from Esrange (European Space and Sounding Rocket Range) near Kiruna 200 km North off the Arctic Circle in northern Sweden (MAXUS-8 mission). The MAXUS rocket is a non-spinning rocket that uses a rate control system to stabilize the payload. With peak accelerations of up to 13 g, the single-stage solid-fuel rocket motor propels the payload above the Earth atmosphere and 86 s after burn-out and motor separation, the parabolic flight maneuver provides microgravity conditions in a quality of <10 −4 g for approximately 840 s reaching an apogee of approximately 705 km. Thereafter, reentry of the payload terminates the microgravity phase, the parachute opens and the payload lands usually 50-75 km north of Esrange in an uninhabited area.
The microscope experiment module flown on MAXUS was accommodated within a 640 mm diameter cylindrical payload element on a platform structure housing the power supply and electronics in the lower part and the microscope and the late-access unit in the upper part. The biological samples were embedded in 1.2% agar in vacuum-tight growth chambers and mounted in a late access unit, which allowed the insertion of the samples into the payload 90-120 min before launch.
A telemetry system multiplexes and transmits experiment and housekeeping data from the experiment modules and the in-flight service systems to the ground. Microscope operations during flight were executed form a ground terminal in the control room via telecommand channels.

Video Analysis
Statolith tracking and recording the position of the statoliths complex during clinorotation and the MAXUS flight was done by outlining the statoliths complex in a sequence of images prepared from AVI-video sequences. The distance between the upper and the lower border of the statoliths complex from the cell tip was measured with ImageJ, Wayne Rasband, National Institutes of Health, Bethesda, USA.

Vibration Measurements
Vibrations were measured with a plane piezo electrical transducer mounted on the inner or outer frame to detect vibrations in the range of 0-25 kHz. Detected signals were recorded via a PC Soundcard and WaveLab 4.0 (Steinberg) and analyzed with Matlab to visualize the frequency spectrum of the measured vibrations.

Results
Under 1g conditions, statoliths in tip-downward growing rhizoids form a disk-like dense complex located in a dynamically stable position above the cell tip. During the 840 s period of high quality microgravity (typically 10 −4 g) provided by the MAXUS-8 sounding rocket flight, the statoliths became more loosely arranged and were displaced away from the tip toward the subapical area of the cell. The distance of the lower border of the statoliths complex increased up to 23.22 μm, the distance of the upper border increased up to 46.13 μm from the tip at the end of the 840s-period of microgravity, before the microgravity period was terminated by vibrations generated by the reentry of the payload module into the Earth's atmosphere.
The kinetic data of statoliths displacement of four Chara rhizoids during the 840 s period of real-microgravity conditions of a MAXUS-flight was taken as reference for evaluating the quality of microgravity simulation provided by a Random-Positioning Machine running with one axis (2D mode) or two axes (3D mode) and different rotational speeds (2D), speed ranges and directions (3D). All rhizoids which were clinorotated for 20 min in the different modes (n = 11 each) continued to grow straight without showing any curvature response which indicates that sedimentation of statoliths triggering a gravitropic growth correction did not occur during any of the clinostat treatments. Like in real microgravity, the statoliths complex was rearranged from a flat disk-like shape into an axially elongated shape further away from the cell tip.

2D Clinorotation
When rhizoids were clinorotated on the RPM (Fig. 1) in 2D mode (using one axis), only the inner frame was used. Control experiments on the outer frame did not reveal any significant differences in the kinetics of statoliths displacement (data not shown). Slow clinorotation with 3 rpm in 2D mode resulted in the increase of the distance of the lower and the upper border of the statoliths complex from 8.12 μm and 18.95 μm up to 22.53 µm and 39.27 μm after 840 s and 26.21 μm and 42.94 μm after 20 min of clinorotation, respectively. The statoliths at the upper border of the complex reached only 84.9% of the distance compared to real microgravity (MAXUS experiments). The statoliths were considerably more dispersed over the rhizoid diameter and formed a more loosely arranged complex under slow clinorotation than in real-microgravity conditions (Fig. 3a).
Increasing the rotational speed of the single axis to 60 rpm strongly promoted the basipetal displacement of statoliths (Fig. 3b). The distances of the lower border and the upper border of the statoliths complex increased up to 25.27 μm and 46.76 μm after 840 s and up to 29.92 μm and 51.68 μm after 20 min of clinorotation, respectively. Thus, the position of the uppermost statoliths was almost identical (101.3%) to the results from the MAXUS real-microgravity experiment and the shape of the complex was also similar to the relatively dense assembly of statoliths in rhizoids und real microgravity conditions.

3D Clinorotation
In the real random mode, both frames rotated with rotational speeds in a range of 2-11 rpm, 30-60 rpm or 60-85 rpm. Statoliths in rhizoids clinorotated at 2-11 rpm were displaced up to 22.42 μm (lower border) and 38.02 (upper border) after 840 s and up to 25.12 μm and 41.30 μm after 20 min, respectively (Fig. 2c). In comparison with the MAXUS data, the statoliths reached a maximal distance of only about 82.4% of the real-microgravity reference. The statoliths were distributed over a relatively large area forming a rather irregularly outlined complex (Fig. 3).
Increasing the range of rotational speed to 30-60 rpm resulted in a basipetal displacement of statoliths up to 20.98 μm (lower border) and 40.13 μm (upper border) after 840 s and 24.18 μm (lower border) and 42.84 μm (upper border) after 20 min, respectively, which represents only 9.0% of the real-microgravity reference value (Fig. 2d).
At a range of 60-85 rpm rotational speed, statoliths were displaced up to 25.92 μm and 42.45 μm after 13 min and 29.59 μm and 45.16 μm after 20 min,  (Fig. 2e). Thus, statoliths reached approx. 86.9% of the real-microgravity reference value. Statoliths always show trembling movements in all directions also in their resting position at 1g as well as during net-basipetal displacement in real-microgravity and in 2D-clinorotated rhizoids. However, trembling motion was considerably Fig. 3 Outlines of statolith complexes in representative rhizoids at 1g (black area) and after 13 min of clinorotation in the 2D and 3D modes in comparison to real microgravity (μg, MAXUS-8 flight) more pronounced during 3D-clinorotation (data not shown). Furthermore, the statoliths complex was considerably less densely packed during 3D-clinorotation than during 2D clinorotation and in real-microgravity (Fig. 3).

Vibration Analysis
For measuring vibration pattern the samples experienced during 2D and 3D operation, the sensor was fixed next to the sample on the inner frame of the RPM. Fast 2D clinorotation at 60 rpm performed with the inner frame of the RPM caused only moderate vibrations in the range of 0-10 kHz (Fig. 4a). In contrast to the 2D mode, 3D clinorotation in real random mode with 60-85 rpm caused much stronger and unregular vibration patterns and considerably higher relative amplitudes over the whole measured frequency spectrum between 0-25 kHz (Fig. 4b).

Discussion
Characean rhizoids were chosen to evaluate the microgravity-simulation quality of the different 2D-and 3D-operation modes of the RPM due to the fact that the position of the numerous statoliths in the subapical region of rhizoids is precisely controlled by the interaction of actomyosin and gravitational forces (Braun et al. 2002). Only in high-quality real microgravity (weightlessness or free fall situation), the basipetal transport of statoliths toward a new equilibrium position further away from the cell tip is exclusively controlled by the actomyosin system. Any remaining accelerations or vibrations inevitably interfere with the transport process and result in impaired transport kinetics and a more dispersed arrangement of statoliths over the cell diameter. This effect has been reported in slow-rotating clinostat studies e.g. by Cai et al. (1997) and in 3D-clinostat studies by Hoson et al. (1997).
In order to evaluate the μg-simulation quality of the different RPM modes, in this study, the kinetics of the basipetal transport of statoliths was analyzed during the initial 20 min of RPM-operation in the different 2D and 3D modes and compared to the kinetics of the statoliths transport during the 13-min real-microgravity phase of MAXUS-sounding rocket flights. The results show that the different velocities of rotation as well as the 2D and 3D mode have considerable and a very distinctive impact on the kinetics of the statoliths transport.
There is convincing evidence that the μg-simulation quality positively correlates with the rotational speed of the RPM. Fast 2D and 3D clinorotation with rotational velocities of 60 rpm and higher caused a displacement of statoliths whose kinetics was very similar to those observed under real-μg conditions of MAXUS sounding rocket flights. Slower clinorotation impaired the basipetal transport and lead to more dispersed statoliths. Earlier studies have shown that more dispersed statoliths sediment more slowly upon gravistimulation, e.g. when the rhizoid is tilted into a horizontal position at 1g leading to a delayed graviresponse in the form of gravitropic curvature which brings the cell's tip back into the tip-downward orientation within a few Fig. 4 Frequency spectrum during fast 2D clinorotation with 60 rpm performed with the inner frame (a) and 3D clinorotation in real random mode with 60-85 rpm (c). Relative amplitudes of vibrations were measured in a range of 0-25 kHz. The intensity of the relative amplitude between −60 (dark blue) and 30 dB (dark red) is color coded (color bar) hours (Herranz et al. 2013). A similar effect was observed after centrifugation of rhizoids which also led to a disturbed sedimentation process most like caused by an adaptation of the actomyosin system .
The best results, namely, a situation closely resembling the basipetal transport of statoliths in real-μg, were achieved with rhizoids rotating precisely along a single axis with a rotational speed of 60 rpm. Using the second axis did not improve the statoliths transport kinetics, however, resulted in generally more loosely arranged statoliths complexes in slower and faster rotated rhizoids. The vibration sensor measurements clearly suggest an overall higher level and irregularity of the vibration amplitudes to be the main causes for the negative influence of a second rotating axis on the statoliths transport process. The erratic vibration pattern of fast 3D-clinorotation appears to result from the constantly changing rotation direction as well as the constantly varying rotational speed.
These results are consistent with reports of RPM experiments with organisms like free-swimming protists like Euglena und Paramecium as well as Drosophila, fish and some human cell types (Herranz et al. 2013). When the RPM is operated in the real-random mode, microgravity simulation is characterized by a constantly changing rotation velocity and in addition also a constantly changing direction of rotation. Clear evidence for the negative impact of the real-random mode on cellular activities comes from experiments with the fast and sensitive shear stress reporter system of dinoflagellates (Hauslage et al. 2017). Shear stress and hydrodynamic gradients trigger a luciferin/luciferase reaction in dinoflagellates probably via a calcium mediated signaling cascade (Mallipattu et al. 2002;Jin et al. 2013;Von Dassow and Latz 2002). Bioluminescence is part of the defense mechanism of these dinoflagellates. The combination of random speed and direction of the PRM had a stronger impact on bioluminescence of the dinoflagellate Dinocystis noctiluca than only changing the rotational speed or using a fast and constantly rotating clinostat at 60 rpm (Hauslage et al. 2017) clearly indicating a higher stress level caused by the 3D operation mode.
The advantage of the fast 2D mode, the lower level of vibrations and the generally better microgravity quality, however, certainly applies only for small objects like dinoflagellates and characean rhizoids as long as they are perfectly aligned along the rotation axis. Good quality microgravity is provided only in a range of a few millimeters along the axis. With increasing speed, the radius around the rotation axis strongly decreases in which good-quality microgravity is provided. However, in contrast to slower rotational speeds (2-11 rpm) which randomize the gravity vector without effectively preventing sedimentation, from a purely technical perspective, higher rotation speeds (> 60 rpm) provide near weightlessness due to the almost complete prevention of any sedimentation process (Briegleb 1992).
Depending on the sensitivity (namely perception and response time of a specific process or organism), slower 2D and 3D clinorotation were applied for investigations on physiological processes like proliferation/growth and development in bigger organisms like Arabidopsis plants (Herranz et al. 2013). At higher rotational speeds, residual g-forces or even centrifugal forces with increasing distance from a fast rotating axis, induce shear forces and might cause stress symptoms in the organisms. The experiments presented in this paper demonstrate that ground-based facilities have to be carefully evaluated for their nongravitational effects on specific organisms and, also, that high-quality microgravity simulation can be achieved supplementing rare and expensive real-microgravity studies in low Earth orbit.
For small objects like the tip of a Chara rhizoids (embedded in agar), where all processes of the gravitropic signaling pathways occur in a very small area, fast clinorotation in 2D and 3D mode provides a good simulated microgravity environment as long as the tip is perfectly oriented in the center of rotation. However, the lower stress and vibration levels make 2D clinorotation the better choice for investigating gravitropic signaling and the gravitropic response in characean rhizoids.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.