- 1Graduate School of Life Science, Hokkaido University, Sapporo, Japan
- 2Research Institute for Electronic Science, Hokkaido University, Sapporo, Japan
- 3Field Science Center for Northern Biosphere, Tomakomai Experimental Forest, Hokkaido University, Tomakomai, Japan
Protists ubiquitously live in nature and play key roles in the food web chain. Their habitats consist of various geometrical structures, such as porous media and rigid surfaces, affecting their motilities. A kind of protist, Stentor coeruleus, exhibits free swimming and adhering for feeding. Under environmental and culture conditions, these organisms are often found in sediments with complex geometries. The determination of anchoring location is essential for their lives. However, the factors that induce the behavioral transition from swimming to adhering are still unknown. In this study, we quantitatively characterized the behavioral transitions in S. coeruleus and observed the behavior in a chamber with dead ends made by a simple structure mimicking the environmental structures. As a result, the cell adheres and feeds in narrow spaces between the structure and the chamber wall. It may be reasonable for the organism to hide itself from predators and capture prey in these spaces. The behavioral strategy for the exploration and exploitation of spaces with a wide variety of geometries in their habitats is discussed.
1 Introduction
Animal navigation is widely observed in nature. For example, birds perform homing behaviors (Landreth and Ferguson, 1967) that have been actively studied in neuroscience (Capaldi et al., 1999), and ants explore the locations of food, which can be explained by simple algorithms such as cellular automata (Deneubourg et al., 1990). Even in bacteria and protists, a navigation ability is required to survive, and their habitats consist of diverse microscopic structures that affect their behaviors. Indeed, the behavioral responses of these organisms have been studied in various microspace geometries: porous media (Bhattacharjee and Datta, 2019), rigid surfaces (Kantsler et al., 2013; Ohmura et al., 2018, 2021), dead ends (Kunita et al., 2014), corners (Théry et al., 2021), interstices between inclined plates (Ishikawa and Kikuchi, 2018) and confined areas (Beppu et al., 2017, 2021; Ostapenko et al., 2018; Bentley et al., 2021). From these studies, it was determined that cell behaviors change in response to the geometrical structures of living spaces.
A kind of protist, namely, ciliates, live in various natural environments, such as fresh water and even the deep sea (Hausmann et al., 2002; Lynn, 2017). These organisms play key roles in the food web chain in the hydrosphere (Sanders and Wickham, 1993; Lynn, 2008). Many kinds of ciliates Vorticella (Sleigh and Barlow, 1976), Opercularia (Sládecek, 1981; Hartmann et al., 2007) and Stentor (Tartar, 1961) anchor to underwater surfaces and consume bacteria, small ciliates and small detritus. Alternatively, they are preyed upon by large multicellular organisms. Thus, their behaviors play an important role in ecosystems in their habitats (Sherr and Sherr, 2002).
Stentor is a large single-celled organism (0.5–1 mm) belonging to the phylum Ciliophora in the order Heterotrichida (Adl et al., 2019); this organism has been well studied for wound healing and regeneration (Tartar, 1961; Blauch et al., 2017; Marshall, 2021; Zhang et al., 2021). The entire body surface of Stentor is covered in cilia for swimming. On the edge of the oral region (peristome), the cilia are also tightly packed (membranelle), i.e., two or three rows of cilia (20–25 cilia in each row (Randall and Jackson, 1958)). The membranelles are components of a ring-like oral structure (membranellar band). The membranelles collectively beat and stir the fluid to generate a feeding vortex (Maupas, 1888), which brings food to the oral apparatus.
Stentor exhibits free swimming. They become sessile once attached to the structures such as dead leaves and twig duckweed, Spirogyra mats and dead cattail leaves (Tartar, 1961). Switching its motility may depend on the intra- and extracellular conditions. Light affects their behaviors; that is, a negative photoresponse can be caused by the colored pigment stentorin under the near cell surface (Tartar, 1961; Song et al., 1980a, 1980b; Iwatsuki and Song, 1989; Iwatsuki, 1992). Recently, Trinh et al. reported that “complex decision-making” in Stentor (Dexter et al., 2019) relates to intracellular conditions (Trinh et al., 2019).
In their behavior, Stentor mainly exhibits three shape states: a droplet, a cone and a trumpet (Tartar, 1961). In free swimming, the cell forms a cone shape; on the other hand, trumpet-shaped cells slowly swim or attach to substances for feeding. Mechanical and electrical stimuli produce a quick contraction from cone or trumpet shapes to droplet shapes as an avoidance behavior. During the adherence process, the Stentor anchors to a surface using an adherent structure (holdfast) at the posterior end, in which mucus is secreted (Andrews, 1945). The determination of anchoring location is essential for capturing prey and hiding from predators. However, it is still unclear what external cues induce the behavioral change from swimming to adhering. In this study, we quantitatively evaluated the effect of extracellular structures on Stentor’s behavior.
2 Materials and methods
2.1 Cell cultures
S. coeruleus was collected from the pond at Tomakomai located Hokkaido, Japan (42.67°N, 141.59°W) and was maintained in a 96-well cell culture plate (TR5003, TrueLine, United States; 6.4 mm diameter) containing filtered pond water using a Millex 0.22 µm fileter (Merck Millipore, Burlington, MA, United States) at room temperature in the dark. We placed Cryptomonas paramecium into the culture plate as the Stentor’s feed, which was axenically grown in a CHM medium [0.1% CH3COONa·3H2O, 0.1% “Lab-Lemco” powder (Oxoid L29) in Milli-Q water; The pH 6.9 ± 0.1 with 1 N NaOH] was adjusted under light-dark conditions (12 h:12 h) at room temperature every three or 4 days. We transferred C. paramecium to a fresh CHM medium every week.
In the experiments using chambers with/without structure, cells were collected from the Shiribetsu River located in Hokkaido, Japan (42.81°N, 140.69°W) and cultured at 25°C in a 6 cm diameter plastic dish containing modified Peters’ solution (0.55 mM CaCl2, 0.15 mM MgSO4, 0.15 mM K2CO3 and 0.75 mM Na2CO3, adjusted pH 7.4 ± 0.1 with 1 N HCl) in the dark (de Terra, 1966) with one wheat grain (MRP-706, Marukan, Osaka, Japan). The cells were transferred into a fresh medium to maintain the cultures every 2–3 weeks. The difference in the strain was not significant in our experiment (detailed in Supplementary Figure S1).
2.2 Observation of the shape transitions of the cell
Before the observations, we washed the cells in fresh pond water twice and starved them for 3 days. Prior to transferring a cell into an observation chamber, we washed the cell with fresh pond water again and equilibrated it for over 20 min. The observation chamber was made of a silicon sheet (244-6012-03, HAGITEC Co., Ltd., Chiba, Japan; 0.2 mm depth, 5 mm inner diameter) and covered by a cover slip (Supplementary Figure S2). Cell behavior was observed by an inverted microscope IX73 (Olympus, Tokyo, Japan) equipped with a PLANPO X1.25 objective (Olympus, Tokyo, Japan). Images were recorded by using a CMOS camera (ORCA-Flash4.0 v2, Hamamatsu Photonics, Hamamatsu, Japan) with an exposure time of 5.0 ms and a frame rate of 10 frames per second (fps) under a brightness field using weak 700 nm light (M700L4, Thorlabs, New Jersey, United States) to reduce the photoresponse effects (Song et al., 1980a, 1980b; Iwatsuki, 1992; Song, 1999).
We used 14 cells to measure the swimming speed and the cell length over 25 min after placing them in the chamber. We recorded each cell for 24.7 ± 3.4 min (SD, 14 cells, maximum 32.0 min, minimum 18.6 min, total 345.1 min: Supplementary Figure S3). To measure the duration time of each shape transition and the diameter of the swimming trajectory in the trumpet shape, we recorded the cell behaviors and analyzed the data (Supplementary Figure S4). Because of the very low frequency at the droplet state, we recorded the behaviors of the cells in the droplet state by adding an external mechanical stimulus to obtain a transition vector and conducted computational classification of the cell state (Supplementary Figure S4).
Due to the speed of the contraction from the trumpet shape to the droplet shape, we recorded the contractions of eight cells by the high-speed camera FASTCAM Mini AX50 (Photron, Tokyo, Japan) at 2,000 fps and an exposure time of 1/2,000 s. Images were collected by using the IX73 microscope with a UNLPlan X10 (Olympus, Tokyo, Japan) objective under 700 nm light.
When we observed the shape transitions from the trumpet shape to the droplet shape, we added an external stimulus by striking the chamber using a tweezer (Meister 0-SA, Rubis, Switzerland).
2.3 Image analyses and measurements of the center of the cell, swimming speed, cell length and swimming trajectory
First, we analyzed the data containing no external stimulus. We cutoff sequential images to several sections (2,000–5,000 frames) and obtained each background image by calculating the local median intensity. Then, the sequential images subtracted from each background were subjected to Gaussian blurring (radius = 1.1 pixels), and binary images were acquired using the MaxEntropy threshold algorithm (Kapur et al., 1985). After that, we obtained the centers Ci: ri = (xi,yi) (mm) of the cells at the ith frame. To measure the cell length, we obtained the length li (θ) from the center of the cell to the edge by converting images to polar coordinates by 2° each (Supplementary Figure S5). The cell length Li was defined by Li = li (θ1) + li (θ2), where li (θ1) and li (θ2) are the top two local maxima of li (θ), and θ1 and θ2 denote the directions of the anterior and posterior directions, respectively (Supplementary Figure S5).
Instantaneous swimming speeds si = |vi| were calculated from the center of the cell Ci, defined by vi = 0.5 [(ri+1−ri-1) · fps] (mm/s). Because of the difficulties of obtaining the shape and the center of the cell due to the thin posterior region and the fusing to the air bubble, we obtained a smoothed time series of the cell length L’i and speed s’i from the median values for ± 5 s (50 frames) for displaying the overview of behaviors; ± 3 s (30 frames) for obtaining the time duration from cone to trumpet; and ± 1 s (10 frames) for obtaining the time duration from droplet to cone and from droplet to trumpet.
In the contraction process from trumpet to droplet by adding an external mechanical stimulus, we manually extracted three points of the posterior end (P) and two edges of the oral apparatus, OA1 and OA2, due to defocusing in striking. Then, we calculated the middle point (M) between OA1 and OA2. After that, the cell length L was determined as the distance of the MP.
2.4 The estimation of the frequencies of each cell state
To obtain an overview of the behaviors and compare them with the computational classification of the cell states (see Section 2.6), we divided the behaviors into three states (droplet, cone and trumpet) with reference to Tartar’s descriptive classification (Tartar, 1961). Then, we obtained the frequency of the cell states and the distributions of the cell length and swimming speed at each state. The normalized length
All image analyses were conducted by using Python 3.7.1 (https://www.python.org/) and ImageJ/Fiji (Schindelin et al., 2012).
2.5 Obtaining the vector field of the behavioral transition of the cell
We described the behaviors of cells in a 2D state field. The field consists of the normalized length
Then, we derived the vector field of the transition of the cell state, and we divided the state field into a grid pattern at regular interval 0.1 (normalized). We calculated the median vector in a grid. We eliminated the vector in the grid less than 0.05%.
2.6 The computational classification of cellular states in the state field
We eliminated the transition periods (±15 s of droplet state, ±60 s of cone state, ±30 s of trumpet state) to classify the steady states of the cells computationally. Hierarchical cluster analysis was carried out to divide the cell states into three clusters using the “AgglomerativeClustering” function included in the “sklearn.cluster” library in Python (Pedregosa et al., 2011). The parameters used in the clustering method were the number of clusters (“nclusters”): 3, “affinity”: Euclidean, “linkage”: complete. Downsampling was applied to the data points (±5 s median values) by 10-increment sampling due to the limitation of the computational memory.
2.7 Measurements of the swimming trajectory and evaluation of the diameter of the local swimming trajectory in the trumpet shape
The swimming trajectories were obtained from the centers of the cell Ci. To compare the trajectories of a cell in the cone state (swimming) and in the trumpet state (rotating), we extracted 400 frames (40 s) of sequential images from each swimming state (20 sections, 14 cells) and sequential images from each trumpet state (2502 ± 901 frames (SD, 20 sections); maximum 4435 frames; minimum 1111 frames; 14 cells).
The evaluation of the diameter of the local swimming trajectory in trumpet shape was conducted from each lap trajectory extracted from the trajectory in the rotating state (133 sections, 11 cells). The diameters of each trajectory were calculated by the maximum distance of the position (Supplementary Figure S6).
2.8 Estimations of the characteristic time of the shape transitions
From the obtained time series of cell length Li, we estimated the characteristic time of duration for the shape transitions. First, we obtained a smoothed time series of cell length L’ i from the median (following the above). The starting points of each shape transition were decided at the cross point of the line before the transition and after the transition (from cone to trumpet), the point of the minimum length from the contraction to after 20 s (from droplet to cone or trumpet) and the point at which external stimuli were added (from trumpet to droplet).
The duration times of each cell deformation were estimated from the deformation rate var ε(t) by fitting the following exponential functions. The deformation rate was defined as ε(t) = [L(t)−L0]/L0, where L0 is the cell length at the starting time of shape transitions in elongation or at the final state in contraction (t = 40 ms). In elongation, the rate was determined by calculating the averaged length for 10 s before each transition. To determine the duration times of each transition τ, the least-squares method (in Python using the package scipy. optimize.curve_fit) was conducted using the following function:
In the contraction process, we used the mechanical model of the damped harmonic oscillator to model the contraction process in Vorticella convallaria (Moriyama et al., 1998). The model included the viscous resistance around the cell. The least-squares method was conducted by using the following function:
The standard deviation of the time series of the strain was calculated for ±1 s in the elongations or ±0.5 ms in the contraction.
In the derivation of the correlation of the cell length and swimming speed, we calculated the normalized values by the maximum value of the speed in the transition is 1, and we defined the length at that time as 1. We fitted these values by power laws in two sections (1–1.25 and 1.25–1.45 about normalized length).
2.9 Preparation of the chambers made by agar to compare of two shapes
The observation chambers were obtained by mold casting in two processes: the preparation of prime molds made by polydimethylsiloxane (PDMS) and the formation of a gel chamber by pouring the agarose gel into PDMS prime molds. We made two kinds of chambers. One is the chamber (diameter 5 mm, depth 0.3 mm) with a disk-like structure (diameter 2.5 mm), and the another is the chamber (diameter 5 mm, depth 0.3 mm ) without any structure (Supplementary Figure S7).
The prime mold data were designed by Fusion 360 software (Autodesk, California, United States). The data were sent to a CAMM-3 Model PNC-3200 milling machine (Roland DG Corporation, Hamamatsu, Japan) by MODELA Player Version 3.7 (Roland DG Corporation, Hamamatsu, Japan), and prime molds were made by cutting a ZW-100 (Roland DG Corporation, Hamamatsu, Japan). Cutting involved three processes: surface cutting, rough cutting and finish cutting. Each cutting parameter is shown in Supplementary Table S1. The parameters in surface cutting and rough cutting were the same. For rough cutting and finish cutting of the chamber with structure, we used the micro end mill, MHRH230 φ0.2 × 2 (NS TOOL, Tokyo, Japan). In other processes, the end mill 2CEM-G-2.5D-2 (MonotaRO, Hyogo Japan) was used.
PDMS was prepared by mixing the prepolymer and the curing agent of the SYLGARD™184 silicone elastomer kit (Dow Corporate, Michigan, United States) in a 10:1 (weight:weight) ratio. The mixture was mixed for 20–30 min by plastic stick by hand. After that, it was vacuum-degassed for approximately 1 h. The mixture was poured into molds and vacuum-degassed again for approximately 1 h. Then, PDMS prime chambers covering a glass slide were horizontally kept at 25°C. After 48 h, the PDMS prime chambers were removed from the molds.
Prior to obtaining the observation chambers, the surface of the PDMS chamber was washed with ethanol and Milli-Q water and sonicated in reverse osmosis water. These processes were repeated three times or more. Then, the surface was rinsed with Milli-Q water. After drying, the surface was coated with plasma (PC-400T, STREX, Osaka, Japan) for 10–20 s for hydrophilization.
The mixture of the 2% (weight per volume) agarose (PrimeGel™ Agarose LE 1–20K, Takara Bio, Shiga, Japan) medium that melted the fresh modified Peters’ solution was poured into PDMS molds and covered by a cover slip to flatten the agar chambers. After 2–3 min, the cover slip was vertically removed, and the agar chambers were withdrawn from the PDMS prime mold. If we could not remove the agar chambers by the above methods, then the chamber was removed by sonication using the ultrasonic cleaner US-5KS (SND, Nagano, Japan) containing reverse osmosis water.
2.10 The observation method in the chamber with/without a structure
Prior to the observations, we washed the cells in a fresh medium twice and starved them for several hours. Then, we placed the gel chambers on a cover slip, and the surface of the chambers was wetted by the medium and covered on a glass slip with a plasma coating (20 s). After that, a glass slide on the chamber was slid until a swimming arena appeared, and the cell was transferred into the chamber from the edge of the glass slide by using a pipette tip. Its inside was coated with the 2% 2-methacryloyloxyethyl phosphorylcholine (MPC) polymer Lipidure® CM5206 (NOF Corporation, Tokyo, Japan) in ethanol (weight per volume) to prevent the cell surface from adhering to the inside of a tip. In this operation, most of the cell quickly contracted due to external stimuli or hydrophilicity of the gel chambers; however, the cell swam after several seconds. Last, we slowly slid the cover slip to cover the observation chamber.
The behaviors were observed by using a stereomicroscope SZX16 (Olympus, Tokyo, Japan) equipped with a C-Mount Adaptor 0.5X, U-TV0.5XC-3 (Olympus, Tokyo, Japan) and a CMOS camera, ORCA-spark C11440-36U (Hamamatsu Photonics, Hamamatsu, Japan). The exposure time was 50 ms at 4 fps. A 2 × 2 binning scheme under a dim white brightness field, filtered by a red right pass plate for approximately 40 min, was utilized. The quantum sensor MQ-610 (Apogee Instruments, Utah, United States) could not detect the observation light due to weakness [<1 μmol s−1 m−2 ∼ 0.2 W/m2 (650 nm)]. Less than 0.2 W/m2 red light has little influence on the phototactic behavior (Song et al., 1980a).
2.11 The analyses of behavioral changes in the chamber with/without a structure
We obtained the behavioral data (center of the cell, swimming speed and cell length) according to Section 2.3. The different points were the obtained background images (using sequential images for 15 min by calculating local maximum intensity), the parameter of the Gaussian blur [radius = 0.5 (pixels)] and the threshold algorithm [triangle methods (Zack et al., 1977)]. When we could not acquire the exact cell position automatically, we manually excluded the noise and dust and prepared the threshold of the circularity or 4 pixel opening process (due to the high circularity of the cell compared with dust). After that, we obtained the position by the automatic method according to Section 2.3. To estimate the cell length, we converted the polar coordinate by 1° each. These acquired time series of length and speed are shown in Supplementary Figures S8, S9.
To calculate the local staying frequency in chambers, we divided the quasi-2D physical fields by 36 mm2 (6 mm × 6 mm) at regular intervals of 56 × 56 grids. Then, the probability density p(xi, yj) in each grid was obtained from the number of observed centers of the cells in each grid. p(xi, yj) was normalized as
Adhesive spots of the cells, duration times and the directions of anterior end adherence were manually extracted from the sequential images. Crescent areas are defined as the region within 1 mm from the tip of the crescent. We did not count the cells staying for less than 10 s as adhesion. All statistical analyses were conducted in Python 3.7.1 using the “scipy.stat” package version 1.1.0.
3 Results
3.1 Cell behaviors in free swimming
After placing them in the chamber, we recorded the behavior of the cells in a quasi-2D disk chamber to perform quantitative measurements of the swimming speed and the cell length over 25 min. The cell exhibits three kinds of shape and behavioral transitions between these states (droplet, cone and trumpet) (Figure 1A). These shapes and transitions are consistent with the previous description (Tartar, 1961). In our experiments, some cells change their states frequently (Figure 1B, Cell A), while others take the steady cone state (Figure 1B, Cell B). There is a variation in the frequency of changing the states among those cells in this situation (Supplementary Figure S3). In the following analyses, we used the behavioral data containing these states. In our observation, cone-shaped cells and trumpet-shaped cells are observed in 72.3% and 25.7% of cells, respectively. Droplet-shaped cells are rarely observed in 2.0% (Figure 1C, Supplementary Figures S13A,B). The distributions of normalized length in each state (cone = 1) are distributed to three regions (Figure 1C). Cone-shaped cells perform free swimming, and droplets and trumpet-shaped cells mostly adhere or swim slowly. Then, we plotted the behavioral data in Figure 1C on the 2D state field. Figure 1D shows that the plots mainly form three clusters corresponding to each state. The vectors in the 2D field represent the velocity of the behavioral change. The graph allows the cell states and transitions in Figure 1A to be evaluated quantitatively and describes the behavioral flows between the cell states. Because two transitions from droplet to trumpet and from trumpet to droplet go through the same path toward mutually opposite directions in the state field, the velocities with respect to the behavioral changes are canceled and underestimated. Hierarchal clustering analysis, which connects the nearest clusters or dots repeatedly and makes several clusters, is conducted on the dots in the steady cell states (Figure 1E). When the computational clustering divides the cell state into three, the result corresponds to the classification in manual clustering (Figures 1D,E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021469/fcell-10-1021469-HTML/image_m/fcell-10-1021469-g001.jpg)
FIGURE 1. Three modes and transitions of S. coeruleus and the variety of behavioral modes in the quasi-2D disk chamber. (A) Schematic images of the cell states and the transitions. The cell takes three shapes (droplet, cone and trumpet). The cell elongates or contracts between these states, corresponding to each arrow. In the shape transition from cone to trumpet, the switching behavior occurs from swimming to adhering to a substrate. (B) The time series of cell length (blue line) and swimming speed (orange line). The graph displays the behaviors of two cells (Cell A, B). The cell states are shown on the top of each graph. Cell A often experiences a shape transition; on the other hand, Cell B remains in a swimming cone state for a longer observation time. (C) The frequency of the cell state and the distributions of the cell length and speed in each state. Each state was manually divided (droplet, 2.0%; cone, 72.3%; trumpet, 25.7%). The distributions of the length and the speed were normalized by calculating averaged values in the cone state. The length in the droplet is shorter than that in the cone. On the other hand, the length in the trumpet shape is higher than that in the cone. Swimming speed is separated into two regions, the swimming state (cone) and the adhering or lower swimming state (droplet and trumpet). (D) The transition vector of the cells. The dots represent the cell states in the 2D field containing the cell length and speed at time points. The speed and length are normalized by cone = 1. The dots were manually classified into three by cell state Tartar’s classification (droplet, purple; cone, orange; trumpet, light blue). (E) Hierarchical clustering of the steady cell states. The dots in the 2D state field are divided into three clusters. The computational clustering is similar to manual classification (D).
3.2 Shape transitions of the cells
To characterize the behavioral transitions accompanied by the shape deformations in Figures 1A,D, we measured the time series of the cell length and estimated the duration times in each shape deformation (Figure 2). The duration time of each transition was estimated by exponential fitting to the averaged time series of the deformation rate. Because the contraction process from cone to droplet was observed only three times in total recording time (over 6 h), the process was eliminated from the analysis. Therefore, the Stentor almost takes the droplet shape caused by the contraction of the trumpet shape. Then, the contracted droplet cell becomes cone shaped and swims away (Figure 2A) or becomes trumpet shaped again (Figure 2C). Additionally, the duration times are 16.2 ± 0.2 s (SD, n = 13, 10 cells, Figure 2B) and 40.1 ± 0.51 s (SD, n = 10, eight cells, Figure 2D). Figure 2E shows another elongation process from cone to trumpet, and the duration time was 107.9 ± 1.47 s (SD, n = 20, 14 cells, Figure 2F). The contractile process from trumpet to droplet is very quick, more than ten thousand times faster than other transitions (Figure 2G). In the contraction process, the two-mode exponential function (Moriyama et al., 1998) fits well [2.1 ± 0.51 ms and 15.6 ± 2.3 ms (SD, n = 8, eight cells, Figure 2H)] with the averaged deformation rate rather than the function with a single time scale (details in Supplementary Figure S10).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021469/fcell-10-1021469-HTML/image_m/fcell-10-1021469-g002.jpg)
FIGURE 2. The shape transitions between the cell states and the estimation of each duration time of the transitions. (A,C,E,G) Sequential images of the shape transitions from droplet to cone (A), from droplet to trumpet (C), from cone to trumpet (E) and from trumpet to droplet (G). The time is shown at the upper left in each image. The beginning time of the shape transition is adjusted at 0 s. In the transitions from droplet to cone (A) and cone to trumpet (E), the cell moves actively. The black curved line in (A) and (E) represents the boundary of the chamber. All bars are 0.2 mm. (B,D,F,H) The time series of shape deformation rates in each transition from droplet to cone (B), from droplet to trumpet (D), from cone to trumpet (F) and from trumpet to droplet (H). The gray regions represent the sections before each transition. The colored solid lines are a time series of the averaged deformation rate. Individual data and the standard deviation of the rate are shown by gray lines and transparent colored regions, respectively. The duration times are estimated by exponential fitting to the averaged deformation rate represented as black dashed lines. These characteristic times are 16.2 ± 0.2 s (SD, n = 13, 10 cells) from droplet to cone (A,B), 40.1 ± 0.51 s (SD, n = 10, eight cells) from droplet to trumpet (C,D), and 107.9 ± 1.47 s (SD, n = 20, 14 cells) from cone to trumpet (E,F). The duration time from trumpet to droplet (G,H) is estimated at 2.1 ± 0.51 ms and 15.6 ± 2.3 ms (SD, n = 8, eight cells).
3.3 Behavioral switching from swimming to adhering
In this section, we focus on the transition from cone to trumpet (Figures 2E,F) related to the determination of the adhering location. In the transition, the swimming speed of the cell decreases, accompanied by increasing cell length (Figures 3A,B). The correlation between the length and speed of the cell shows a negative correlation, and two kinds of relationships are observed in the transition by scaling methods (s∼l−6.3 and s∼l0.2, Figure 3B). In the behavioral transition, the swimming trajectory changed from straight forward to rotating (Figure 3C, Supplementary Video S1). The trajectories in the swimming and trumpet states are shown in Figures 3D,E. The true diameter of the trajectory in the free-swimming cone state could not be measured because the cone-shaped cell usually swims in a straight trajectory and the swimming is affected by the chamber’s wall due to the restriction of the observation area. The diameter of the rotating trajectory of the trumpet cell (Figure 3E) is distributed at 1.4 ± 0.3 mm (SD, n = 133, 11 cells, Figure 3F). The parameters of the beating forms of the membranellar band are different between the cone state and trumpet state (Supplementary Figure S11).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021469/fcell-10-1021469-HTML/image_m/fcell-10-1021469-g003.jpg)
FIGURE 3. The processes of behavioral switching from swimming (cone shaped) to adhering (trumpet shaped). (A) Time series of the cell length (blue line) and swimming speed (orange line) of the cell in the transition from cone to trumpet (Figures 2E,F). Averaged lines are represented as solid bold lines, and individual data are represented as solid thin lines. Standard deviations of values for each cell at a time point are the transparent regions. These values are normalized by Cone = 1. It shows the negative correlation between the cell length and swimming speed. (B) The correlation between cell length and speed. The green solid line represents the average speed. Gray lines are individual transition data, and averaged data are fitted by power laws. The characteristic factors of the law are −6.3 in the range of 1–1.25 (normalized length) and 0.2 in the range of 1.25–1.45 (normalized length). (C) The trajectory of a cell in the behavioral transition from swimming to adhering. As the cell elongates, the cell becomes locally rotating. The color represents elapsed time from the beginning of the transition. The dashed line represents the wall of the chamber. (D,E) The comparison between the swimming trajectories at swimming steady state and steady rotating state. (F) The distribution of the swimming pass diameter in the steady rotating state. The diameter was 1.4 ± 0.3 mm (SD, n = 133, 11 cells).
3.4 The influence of a structure on behaviors in an arena
In nature, Stentor elongates and adheres to a better location for better living conditions. In these processes, the cell obtains some information about the outer environment. To determine where they tend to adhere, we focused on the geometry in the outer swimming environment and observed their behaviors in the quasi-2D gel chambers to encourage cell adhesion with/without the structure (Figure 4A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1021469/fcell-10-1021469-HTML/image_m/fcell-10-1021469-g004.jpg)
FIGURE 4. The behavioral change in swimming to adhesion in the crescent areas. (A) Image of the geometry of the gel chamber with the structure in the experiments. The dashed line represents the wall of the chamber. Crescent areas are indicated as two orange regions. The depth of the chamber is 0.3 mm. (B) The frequency of the cell states in the gel chamber with/without the structure. In the chamber without the structure (control), the cell mainly takes the cone state in 74.1% (orange), and trumpet (light blue) and droplet (purple) states are observed in 22.4% and 3.5%, respectively. In the chamber with the structure, the cell mainly takes a trumpet shape in 72.4% and a droplet shape in 10.7%. Cone-shaped cells are observed in only 16.9% of cells. (C) The difference in the probability of three states in the simple chamber and in the chamber with structure. The colors correspond with each state indicated in (B). This bar plot exhibits the probability of cell states in the simple chamber and the chamber with structure. The cone shape was dominant in the simple chambers, but the trumpet type was dominant in the chambers containing structures. The probability of the cell states was tested by Welch’s t-test, revealing significant differences between the chamber with/without structure (p = 0.0026 in droplet; p = 6.3 × 10−7 in cone; p = 2.0 × 10−6 in trumpet). The bars are standard errors (12 cells). (D,E) 2D maps of the probability density observed in the cells in each chamber. (F) The comparison of presence probability within the crescent areas and other regions. Except for adhering, presence probabilities are 614.9 ± 136.6 s/mm2 (SE, n = 12, 12 cells) in crescent and 43.4 ± 11.1 s/mm2 (SE, n = 12, 12 cells) in other regions. The difference is statistically significant in Welch’s t-test. Black bars are standard errors. (G) The locations and duration times of adhesion in the cell. The color and size of the dots represent the duration of the adhesion. The color bar represents the logarithmic adherence time scale. (H) The directions of the cell in the adhesive state. The color and length of the vectors represent the adhering time in one direction. (I) The comparison of adhering time prior to the detachment within the crescent areas and other regions. Ten cells adhered in crescent areas (18 events) and 7 cells adhered in other regions (13 events) in observed 12 cells (31 events). Averaged adhering times prior to detachment are 603.0 ± 209.2 s (SE, n = 18, 10 cells) in the crescent areas and 411.2 ± 187.1 s (SE, n = 13, seven cells) in other regions. Welch’s t-test revealed that they did not have a significant difference (p = 0.50). Black bars are standard errors. (J–L) Typical sequential images of the process from swimming to adhesion in the crescent area [in the rectangle in (A)]. (J) Before the cell adheres, the cell enters the crescent area. (K) Adhesive cells become trumpet-shaped from cone-shaped cells. (L) The cell adheres and feeds in the crescent area. The direction of the anterior end directs toward the opposite side to the dead end. All bars are 0.5 mm. (M) Time series of the cell length and swimming speed in the process from swimming to adhesion at the crescent area. Images in (J)–(L) represent the sections in the graph.
The structure in the chamber causes a change in the frequency of the cell states. In the chamber without the structure (control), the cell mainly takes the cone state in 74.1%, and trumpet and droplet states are observed for 22.4% and 3.5%, respectively (Figure 4B). However, in the chamber with the structure, the cell mainly takes a trumpet shape in 72.4% and a droplet shape in 10.7%. Cone-shaped cells were observed in only 16.9% (Figure 4B). In the chamber with the structure, trumpet- and droplet-shaped cells mostly adhere, and swimming cells often collide with the structure. Thus, a distribution of lower swimming speeds is often observed (Supplementary Figure S12A,C). The averaged each state have significant differences in Welch’s t-test (p = 0.0026 in droplet, p = 6.3 × 10−7 in cone, p = 2.0 × 10−6 in trumpet, Figure 4C). The geometry in the outer environment affects the modulation of their states.
The 2D color maps in Figures 4D,E represent the averaged probability densities of the cell’s position in each chamber. In the quasi-2D chamber without a structure, the probability density is uniformly distributed along the wall (Figure 4D). On the other hand, the cell frequently stays in the locations between the wall and the structure, which are termed the crescent areas (in the orange area in Figure 4A) (Figure 4E). Except for adhering, existence frequencies were 614.9 ± 136.6 s/mm2 (SE, n = 12, 12 cells) in crescent and 43.4 ± 11.1 s/mm2 (SE, n = 12, 12 cells) in other regions. Welch’s t-test was applied to the values, revealing a significant difference between them (p = 0.0015, Figure 4F). It shows the increase of the presence in the crescent areas.
Our experiment suggests that the cell tends to stay in and adhere to the crescent areas. Figure 4G displays the positions and duration times of adhesion in 12 cells. All 12 cells adhere to substrate in the chamber with the structure. Though two cells did not adhere to the crescent areas, the others (10 cells) exhibited switching behavior from swimming to adhesion (from cone to trumpet) in crescent areas. In the 10 cells, 9 cells maintained their feeding (trumpet state) with the oral structure directed to the opposite side of the dead end (Figures 4G,H,L). Although the adhering frequency in crescent areas (28.0 events/mm2) increases compared to the one in other regions (0.9 events/mm2), the difference between adhering times per cell prior to the detachment were not significant in Welch’s t-test [603.0 ± 209.2 s (SE, n = 18, 10 cells) in crescent areas and 411.2 ± 187.1 s (SE, n = 13, seven cells) in other regions; p = 0.50; Figure 4I].
Figures 4J–L shows typical sequential images of the adhesion processes in the crescent area (Supplementary Video S2). First, the cell encounters the crescent area and goes back and forth (Figure 4J). Then, the posterior shortly elongates and contacts the wall. The cell sometimes swims backward. Then, the posterior region adheres to the substrate and elongates over 1–2 min (Figure 4K). In this process, we can sometimes observe the bending behavior (sequential images from 96 to 120 s in Figure 4K). Although contraction and elongation between the trumpet and droplet occur, the cell typically remains adherent at the same location for more than 35 min (Figure 4M, Supplementary Figure S9).
Modulation of their states is affected by the structure in an arena. In the crescent areas, the cell not only stays frequently but also switches its behavior and adheres for a long time with its oral structure directed to the opposite side of the dead end.
4 Discussion
Protists perform some behavioral modulations according to their surrounding environments. The ciliate Stentor exhibits changing shape and motility, and several studies have shown the complexity of these behavioral transitions (Jennings, 1906; Wood, 1969b; Dexter et al., 2019; Trinh et al., 2019). The environmental information around the cell affects the behavioral transitions. Here, we quantitatively evaluated the behaviors of S. coeruleus and then measured the characteristic durations of their behavioral transitions. Considering the swimming speed and cell length, the cell states are computationally classified into three states, consistent with the previous descriptive classification of the three states (Tartar, 1961). Among those three states, we revealed that S. coeruleus has 4 transition processes accompanied by each characteristic duration time. Previously, the durations time from droplet to trumpet and trumpet to droplet have been estimated by qualitative observations, and our results are consistent with their reports (Wood, 1969a; Kristensen et al., 1974). The contraction process in S. coeruleus from trumpet to droplet contributes to the contractile filaments, myoneme (Randall and Jackson, 1958; Huang and Pitelka, 1973). The contraction time in Vorticella convallaria (Moriyama et al., 1998), which also uses a similar contractile filament, the spasmoneme (Lynn, 2008) is almost the same as that in S. coeruleus. In this paper, we also characterized the duration of the elongation processes from droplet to cone and from cone to trumpet. The elongation processes are associated with the sliding of microtubule ribbons in the km-fibers (Randall and Jackson, 1958; Huang and Pitelka, 1973). In our study, the duration times of each elongation were different. This means that the question of whether all three elongation processes relate to sliding in the km-fibers is still open.
In the transition from cone to trumpet, we found a negative correlation between the cell length l and swimming speed s with the power laws s∼l−6.3 and s∼l0.2 (Figure 3B). Hydrodynamically, in a moving sphere (radius r) with the same energy, the swimming speed is proportional to r−1 by Stokes law. We found that the exponent between observed length and speed is different from that derived from the stokes law. It is considered that the difference in the relationship between our experiment and the theory is derived from the difference in the ciliary waveforms between the cone state and the trumpet state (Supplementary Figure S11, Supplementary Video S3).
The changes in swimming speed and trajectory in S. coeruleus play crucial roles in exploring feeding location. In microorganisms, the modulations induced by taxis-caused chemical substances, light and flows have been well studied experimentally and mathematically (Wan and Jékely, 2021). In Stentor, photoresponse (Song et al., 1980a; Song et al., 1980b; Iwatsuki, 1992), external mechanical stimuli (Wood, 1969b) and chemical avoidance responses (Jennings, 1906; Dexter et al., 2019; Trinh et al., 2019) have been reported previously. The behaviors of microorganisms are strongly influenced by extracellular geometries (Kantsler et al., 2013; Kunita et al., 2014; Beppu et al., 2017; Beppu et al.,2021; Ishikawa and Kikuchi, 2018; Nishigami et al., 2018; Ohmura et al., 2018; Ohmura et al., 2021; Ostapenko et al., 2018; Bhattacharjee and Datta, 2019; Bentley et al., 2021; Okuyama et al., 2021; Théry et al., 2021) corresponding to mud, dead leaves or other varieties of sediments in their habitats. In our study, we evaluated the effect of a structure on behavioral transition. In the chamber with the crescent areas, the probability density of the cell’s position at the crescent areas is higher, and the cell becomes more adherent than the cells at other locations. Although the frequency of presence except for adhering had a significant difference, adhering time per cell prior to detachment did not have a significant difference. These statistical analyses indicated that the bias to adhere to the crescent was dominated by an increase in adhesion than a decrease in detachment.
Previously, a directional swimming organism such a Chlamydomonas accumulates the corner in the geometrical chamber, and the phenomenon has been described as a theoretical model (Théry et al., 2021). Although the cell accumulates in the corner without changing the behavioral state (Mondal et al., 2021), we revealed that S. coeruleus, which is frequently found in a free-swimming form (Reynierse and Walsh, 1967), tended to adhere to the chamber with the structure, and the behavioral change promoted accumulation in the crescent areas. The anterior end of the cell faces the opposite side of the dead end. Hydrodynamically, S. coeruleus, which belongs to the puller-like swimmer in the squirmer model, tends to swim away from the wall (Lighthill, 1952; Blake, 1971; Lauga and Powers, 2009). When the posterior region of the cell, which secretes a sticky substance, attaches to a substrate (Andrews, 1945), the cell directs the anterior toward the opposite side of the dead end; this plays a role in positive feeding efficiency.
Moreover, it is known that the Stentor causes two feeding vortices at the left and right sides of the oral apparatus on a focal plane (Maupas, 1888). Proximity to the attachment surface or substrate causes vortex flow near the sessile filter feeders, Vorticella and Opercularia, generating feeding currents around the oral apparatus similar to Stentor (Liron and Blake, 1981; Hartmann et al., 2007; Pepper et al., 2010, 2013). The feeding efficiency influenced by the feeding vortex depends on the environmental situations and theoretical definitions. The eddies contribute to decreasing the long-term feeding efficiency defined by the clearance rate (the volume of nutrient water crossed per time) due to recirculation when we consider that there is no nutrient in the flow processed the oral region once (Liron and Blake, 1981; Hartmann et al., 2007; Pepper et al., 2013, 2021). On the other hand, when we modeled that the organisms probabilistically capture prey, the vortex array contributes to the positive efficiency of their feedings because the recirculation allows an increase in feeding chances of starfish larvae (Gilpin et al., 2017). Additionally, vortex flow enhances the mixing and transport of substances, i.e., nutrients and wastes, around the cell (Mondal et al., 2021). Combining our results and these reports, we determined that S. coeruleus adheres to narrow areas not only by avoiding being the prey of potential predators but also by capturing prey and mixing nutrients around the cell.
How organisms perceive the narrowness and shape of space is an interesting future topic. Two possibilities are discussed here: 1) the interaction between the flow in the narrow space and the ciliary beating that causes it, 2) the concentration of some chemical secretions (due to suppression of diffusion by narrow space). The Stentor swims by moving their cilia to move the water around the body. Ciliary beating is modulated in response to the reaction force cilia receive from water. This modulation is the interaction of the cilia with the surrounding flow of water (Brumley et al., 2015; Wan et al., 2020). In a narrow space, the surrounding flow changes, so it may react to the change in the flow field around the body (Liron and Blake, 1981; Pepper et al., 2010). Sensing a geometric shape may also be possible in principle since the flow field pattern changes with the spatial shape. In other words, it is a mechanism that can be called hydrodynamic sensing for spatial shape. On the other hand, the second possibility is chemical diffusion sensing. The Stentor secretes various chemical substances, such as signal substances and excretory substances (Andrews, 1945, 1946; Miyake et al., 2001). These chemicals move away from the body by diffusion and advection if the space is sufficiently open. However, in a narrow space, these outflows are suppressed, and the concentration becomes high. This change in a chemical concentration may be a cue for the organism to switch the swimming behavior.
The above two effects do not have to be mutually exclusive, and both may work as possible. Alternatively, in a narrow space where the flow tends to stagnate, bacteria and other microorganisms are likely to grow, and it is possible that the Stentor will react to it. In any case, the ability of Stentor to perceive the narrowness and shape of space simply by moving around in space, in other words, without perceiving the overall shape from a bird’s-eye view, is surprising in terms of the information processing ability of living things.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
SE, KS, TN, and YN contributed to conceptualizing the research and designing the experiments. SE, YN, and OK performed the experiments. SE acquired and analyzed the data. SE wrote the draft of the manuscript. SE, KS, OK, TN, and YN contributed to manuscript revision and read and approved the submitted version.
Funding
This research was supported by JSPS KAKENHI Grant Numbers, JP16H06280 (TN), JP21H05308 (YN), JP21H05310 (TN), JP19KK0180 (YN), and JP20K03871 (KS). A part of this study was supported by the “Dynamic Alliance for Open Innovation Bridging Human, Environment and Materials” of the Ministry of Education, Culture, Sports, Science and Technology of Japan (TN, KS, and YN), the JSPS Core-to-Core Program, A. Advanced Research Networks (YN), the Project of Junior Scientist Promotion in Hokkaido University (YN), Photo-Excitonix Project in Hokkaido University (TN), Hokkaido University Support Program for Frontier Research (YN), Global Station for Soft Matter at Hokkaido University (KS, TN and YN), the Sasakawa Scientific Research Grant from The Japan Science Society Grant No. 2021-6029 (SE), and the establishment of university fellowships towards the creation of science technology innovation, Grant Number JPMJFS2101 (SE).
Acknowledgments
We thank Takuji Ishikawa and Kenta Ishimoto for helpful discussions and comments from theoretical aspects. We also thank Taku Morimoto and the Division of Technical Staffs in the Research Institute for Electronic Science, Hokkaido University for supporting the fabrication technique.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.1021469/full#supplementary-material
References
Adl, S. M., Bass, D., Lane, C. E., Lukeš, J., Schoch, C. L., Smirnov, A., et al. (2019). Revisions to the classification, nomenclature, and diversity of eukaryotes. J. Eukaryot. Microbiol. 66, 4–119. doi:10.1111/jeu.12691
Andrews, E. A. (1946). Ingestion organs in folliculinids and in stentors. J. Morphol. 79, 419–444. doi:10.1002/jmor.1050790304
Andrews, E. A. (1945). Stentor’s anchoring organs. J. Morphol. 77, 219–232. doi:10.1002/jmor.1050770206
Bentley, S. A., Anagnostidis, V., Schlogelhofer, H. L., Gielen, F., and Wan, K. Y. (2021). Phenotyping single-cell motility in microfluidic confinement. bioRxiv[preprint]. doi:10.1101/2021.12.24.474109
Beppu, K., Izri, Z., Gohya, J., Eto, K., Ichikawa, M., and Maeda, Y. T. (2017). Geometry-driven collective ordering of bacterial vortices. Soft Matter 13, 5038–5043. doi:10.1039/c7sm00999b
Beppu, K., Izri, Z., Sato, T., Yamanishi, Y., Sumino, Y., and Maeda, Y. T. (2021). Edge current and pairing order transition in chiral bacterial vortices. Proc. Natl. Acad. Sci. U. S. A. 118, e2107461118. doi:10.1073/PNAS.2107461118
Bhattacharjee, T., and Datta, S. S. (2019). Bacterial hopping and trapping in porous media. Nat. Commun. 10, 2075–2110. doi:10.1038/s41467-019-10115-1
Blake, J. R. (1971). A spherical envelope approach to ciliary propulsion. J. Fluid Mech. 46, 199–208. doi:10.1017/S002211207100048X
Blauch, L. R., Gai, Y., Khor, J. W., Sood, P., Marshall, W. F., and Tang, S. K. Y. (2017). Microfluidic guillotine for single-cell wound repair studies. Proc. Natl. Acad. Sci. U. S. A. 114, 7283–7288. doi:10.1073/pnas.1705059114
Brumley, D. R., Polin, M., Pedley, T. J., and Goldstein, R. E. (2015). Metachronal waves in the flagellar beating of Volvox and their hydrodynamic origin. J. R. Soc. Interface 12, 20141358. doi:10.1098/rsif.2014.1358
Capaldi, E. A., Robinson, G. E., and Fahrbach, S. E. (1999). Neuroethology of spatial learning: The birds and the bees. Annu. Rev. Psychol. 50, 651–682. doi:10.1146/annurev.psych.50.1.651
de Terra, N. (1966). Culture of stentor coeruleus on colpidium campylum. J. Protozool. 13, 491–492. doi:10.1111/j.1550-7408.1966.tb01947.x
Deneubourg, J. L., Aron, S., Goss, S., and Pasteels, J. M. (1990). The self-organizing exploratory pattern of the argentine ant. J. Insect Behav. 3, 159–168. doi:10.1007/BF01417909
Dexter, J. P., Prabakaran, S., and Gunawardena, J. (2019). A complex hierarchy of avoidance behaviors in a single-cell eukaryote. Curr. Biol. 29, 4323–4329.e2. doi:10.1016/j.cub.2019.10.059
Gilpin, W., Prakash, V. N., and Prakash, M. (2017). Vortex arrays and ciliary tangles underlie the feeding-swimming trade-off in starfish larvae. Nat. Phys. 13, 380–386. doi:10.1038/nphys3981
Hartmann, C., Özmutlu, Ö., Petermeier, H., Fried, J., and Delgado, A. (2007). Analysis of the flow field induced by the sessile peritrichous ciliate Opercularia asymmetrica. J. Biomech. 40, 137–148. doi:10.1016/j.jbiomech.2005.11.006
Hausmann, K., Hülsmann, N., Polianski, I., Schade, S., and Weitere, M. (2002). Composition of benthic protozoan communities along a depth transect in the eastern Mediterranean Sea. Deep Sea Res. Part I Oceanogr. Res. Pap. 49, 1959–1970. doi:10.1016/S0967-0637(02)00095-X
Huang, B., and Pitelka, D. R. (1973). The contractile process in the ciliate, Stentor coeruleus: I. The role of microtubules and filaments. J. Cell. Biol. 57, 704–728. doi:10.1083/jcb.57.3.704
Ishikawa, T., and Kikuchi, K. (2018). Biomechanics of Tetrahymena escaping from a dead end. Proc. Biol. Sci. 285, 20172368. doi:10.1098/rspb.2017.2368
Iwatsuki, K., and Song, P. (1989). The ratio of extracellular Ca2+ to K+ ions affects the photoresponses in Stentor coeruleus. Comp. Biochem. Physiol. A Comp. Physiol. 92, 101–106. doi:10.1016/0300-9629(89)90749-4
Iwatsuki, K. (1992). Stentor coeruleus shows positive photokinesis. Photochem. Photobiol. 55, 469–471. doi:10.1111/j.1751-1097.1992.tb04264.x
Jennings, H. S. (1906). Behavior of the lower organisms. New York: Columbia University Press. doi:10.1037/10817-000
Kantsler, V., Dunkel, J., Polin, M., and Goldstein, R. E. (2013). Ciliary contact interactions dominate surface scattering of swimming eukaryotes. Proc. Natl. Acad. Sci. U. S. A. 110, 1187–1192. doi:10.1073/pnas.1210548110
Kapur, J. N., Sahoo, P. K., and Wong, A. K. C. (1985). A new method for gray-level picture thresholding using the entropy of the histogram. Comput. Vis. Graph. Image Process. 29, 273–285. doi:10.1016/0734-189X(85)90125-2
Kristensen, B. I., Nielsen, L. E., and Rostgaard, J. (1974). Variations in myoneme birefringence in relation to length changes in Stentor coeruleus. Exp. Cell. Res. 85, 127–135. doi:10.1016/0014-4827(74)90222-5
Kunita, I., Kuroda, S., Ohki, K., and Nakagaki, T. (2014). Attempts to retreat from a dead-ended long capillary by backward swimming in Paramecium. Front. Microbiol. 5, 270. doi:10.3389/fmicb.2014.00270
Landreth, H. F., and Ferguson, D. E. (1967). Newt orientation by sun-compass. Nature 215, 516–518. doi:10.1038/215516a0
Lauga, E., and Powers, T. R. (2009). The hydrodynamics of swimming microorganisms. Rep. Prog. Phys. 72, 096601. doi:10.1088/0034-4885/72/9/096601
Lighthill, M. J. (1952). On the squirming motion of nearly spherical deformable bodies through liquids at very small Reynolds numbers. Commun. Pure Appl. Math. 5, 109–118. doi:10.1002/cpa.3160050201
Liron, N., and Blake, J. R. (1981). Existence of viscous eddies near boundaries. J. Fluid Mech. 107, 109–129. doi:10.1017/S0022112081001699
Lynn, D. H. (2017). “Ciliophora,” in Handbook of the protists (Switzerland: Springer International Publishing AG), 679–730. doi:10.1007/978-3-319-28149-0
Lynn, D. H. (2008). The ciliated Protozoa. 3rd ed. Dordrecht: Springer Netherlands. doi:10.1007/978-1-4020-8239-9
Marshall, W. F. (2021). Regeneration in Stentor coeruleus. Front. Cell. Dev. Biol. 9, 753625–753713. doi:10.3389/fcell.2021.753625
Maupas, E. (1888). Recherches experimentales sur la multiplication des infusoires cilies. Arch. Zool. Exp. Gen. 46, 498–516. Available at: https://www.biodiversitylibrary.org/bibliography/79165.
Miyake, A., Harumoto, T., and Iio, H. (2001). Defence function of pigment granules in Stentor coeruleus. Eur. J. Protistol. 37, 77–88. doi:10.1078/0932-4739-00809
Mondal, D., Prabhune, A. G., Ramaswamy, S., and Sharma, P. (2021). Strong confinement of active microalgae leads to inversion of vortex flow and enhanced mixing. eLife 10, 676633–e67723. doi:10.7554/eLife.67663
Moriyama, Y., Hiyama, S., and Asai, H. (1998). High-speed video cinematographic demonstration of stalk and zooid contraction of Vorticella convallaria. Biophys. J. 74, 487–491. doi:10.1016/S0006-3495(98)77806-3
Nishigami, Y., Ohmura, T., Taniguchi, A., Nonaka, S., Manabe, J., Ishikawa, T., et al. (2018). Influence of cellular shape on sliding behavior of ciliates. Commun. Integr. Biol. 11, 15066666–e1506675. doi:10.1080/19420889.2018.1506666
Ohmura, T., Nishigami, Y., Taniguchi, A., Nonaka, S., Ishikawa, T., and Ichikawa, M. (2021). Near-wall rheotaxis of the ciliate Tetrahymena induced by the kinesthetic sensing of cilia. Sci. Adv. 7, eabi5878. doi:10.1126/sciadv.abi5878
Ohmura, T., Nishigami, Y., Taniguchi, A., Nonaka, S., Manabe, J., Ishikawa, T., et al. (2018). Simple mechanosense and response of cilia motion reveal the intrinsic habits of ciliates. Proc. Natl. Acad. Sci. U. S. A. 115, 3231–3236. doi:10.1073/pnas.1718294115
Okuyama, K., Nishigami, Y., Ohmura, T., and Ichikawa, M. (2021). Accumulation of Tetrahymena pyriformis on interfaces. Micromachines 12, 1339. doi:10.3390/mi12111339
Ostapenko, T., Schwarzendahl, F. J., Böddeker, T. J., Kreis, C. T., Cammann, J., Mazza, M. G., et al. (2018). Curvature-guided motility of microalgae in geometric confinement. Phys. Rev. Lett. 120, 068002. doi:10.1103/PhysRevLett.120.068002
Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., et al. (2011). Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 12, 2825–2830. Available at: https://scikit-learn.org/stable/about.html#citing-scikit-learn.
Pepper, R. E., Riley, E. E., Baron, M., Hurot, T., Nielsen, L. T., Koehl, M. A. R., et al. (2021). The effect of external flow on the feeding currents of sessile microorganisms. J. R. Soc. Interface 18, 20200953. doi:10.1098/rsif.2020.0953
Pepper, R. E., Roper, M., Ryu, S., Matsudaira, P., and Stone, H. A. (2010). Nearby boundaries create eddies near microscopic filter feeders. J. R. Soc. Interface 7, 851–862. doi:10.1098/rsif.2009.0419
Pepper, R. E., Roper, M., Ryu, S., Matsumoto, N., Nagai, M., and Stone, H. A. (2013). A new angle on microscopic suspension feeders near boundaries. Biophys. J. 105, 1796–1804. doi:10.1016/j.bpj.2013.08.029
Randall, J. T., and Jackson, S. F. (1958). Fine structure and function in Stentor polymorphus. J. Biophys. Biochem. Cytol. 4, 807–830. doi:10.1083/jcb.4.6.807
Reynierse, J. H., and Walsh, G. L. (1967). Behavior modification in the protozoan Stentor re-examined. Psychol. Rec. 17, 161–165. doi:10.1007/BF03393700
Sanders, R. W., and Wickham, S. A. (1993). Planktonic protozoa and metazoa: Predation, food quality and population control. Mar. Microb. Food Webs 7, 197–223. Available at: http://hdl.handle.net/11858/00-001M-0000-000F-E43A-E.
Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682. doi:10.1038/nmeth.2019
Sherr, E. B., and Sherr, B. F. (2002). Significance of predation by protists in aquatic microbial food webs. Ant. Van Leeuwenhoek 81, 293–308. doi:10.1023/A:1020591307260
Sládecek, V. (1981). Indicator value of the genus Opercularia (ciliata). Hydrobiologia 79, 229–232. doi:10.1007/BF00006317
Sleigh, M. A., and Barlow, D. (1976). Collection of food by Vorticella. Trans. Am. Microsc. Soc. 95, 482. doi:10.2307/3225140
Song, P.-S., Häder, D.-P., and Poff, K. L. (1980a). Phototactic orientation by the ciliate, Stentor coeruleus. Photochem. Photobiol. 32, 781–786. doi:10.1111/j.1751-1097.1980.tb04055.x
Song, P.-S., Häder, D.-P., and Poff, K. L. (1980b). Step-up photophobic response in the ciliate, Stentor coeruleus. Arch. Microbiol. 126, 181–186. doi:10.1007/BF00511225
Song, P.-S. (1999). Photosensory signal transduction in Stentor coeruleus. J. Photoscience 6, 37–45.
Théry, A., Wang, Y., Dvoriashyna, M., Eloy, C., Elias, F., and Lauga, E. (2021). Rebound and scattering of motile: Chlamydomonas algae in confined chambers. Soft Matter 17, 4857–4873. doi:10.1039/d0sm02207a
Trinh, M. K., Wayland, M. T., and Prabakaran, S. (2019). Behavioural analysis of single-cell aneural ciliate, Stentor roeseli, using machine learning approaches. J. R. Soc. Interface 16, 20190410. doi:10.1098/rsif.2019.0410
Wan, K. Y., Hürlimann, S. K., Fenix, A. M., McGillivary, R. M., Makushok, T., Burns, E., et al. (2020). Reorganization of complex ciliary flows around regenerating Stentor coeruleus. Philos. Trans. R. Soc. Lond. B Biol. Sci. 375, 20190167. doi:10.1098/rstb.2019.0167
Wan, K. Y., and Jékely, G. (2021). Origins of eukaryotic excitability. Philos. Trans. R. Soc. Lond. B Biol. Sci. 376, 20190758. doi:10.1098/rstb.2019.0758
Wood, D. C. (1969a). Electrophysiological studies of the protozoan, Stentor coeruleus. J. Neurobiol. 1, 363–377. doi:10.1002/neu.480010402
Wood, D. C. (1969b). Parametric studies of the response decrement produced by mechanical stimuli in the protozoan, Stentor coeruleus. J. Neurobiol. 1, 345–360. doi:10.1002/neu.480010309
Zack, G. W., Rogers, W. E., and Latt, S. A. (1977). Automatic measurement of sister chromatid exchange frequency. J. Histochem. Cytochem. 25, 741–753. doi:10.1177/25.7.70454
Keywords: geometrical cue, ciliates, Stentor, spatial exploration, behavioral transition, protists
Citation: Echigoya S, Sato K, Kishida O, Nakagaki T and Nishigami Y (2022) Switching of behavioral modes and their modulation by a geometrical cue in the ciliate Stentor coeruleus. Front. Cell Dev. Biol. 10:1021469. doi: 10.3389/fcell.2022.1021469
Received: 17 August 2022; Accepted: 17 October 2022;
Published: 01 November 2022.
Edited by:
Helena Soares, Fac. Ciências Universidade de Lisboa, PortugalReviewed by:
Wallace Marshall, University of California, San Francisco, United StatesAnne-Marie Tassin, Centre National de la Recherche Scientifique (CNRS), France
Judith Van Houten, University of Vermont, United States
Copyright © 2022 Echigoya, Sato, Kishida, Nakagaki and Nishigami. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yukinori Nishigami, nishigami@es.hokudai.ac.jp