Abstract
Free full text
Substrate composition directs slime molds behavior
Associated Data
Abstract
Cells, including unicellulars, are highly sensitive to external constraints from their environment. Amoeboid cells change their cell shape during locomotion and in response to external stimuli. Physarum polycephalum is a large multinucleated amoeboid cell that extends and develops pseudopods. In this paper, changes in cell behavior and shape were measured during the exploration of homogenous and non-homogenous environments that presented neutral, and nutritive and/or adverse substances. In the first place, we developed a fully automated image analysis method to measure quantitatively changes in both migration and shape. Then we measured various metrics that describe the area covered, the exploration dynamics, the migration rate and the slime mold shape. Our results show that: (1) Not only the nature, but also the spatial distribution of chemical substances affect the exploration behavior of slime molds; (2) Nutritive and adverse substances both slow down the exploration and prevent the formation of pseudopods; and (3) Slime mold placed in an adverse environment preferentially occupies previously explored areas rather than unexplored areas using mucus secretion as a buffer. Our results also show that slime molds migrate at a rate governed by the substrate up until they get within a critical distance to chemical substances.
Introduction
Large-scale spatial patterns in biology are common and knowing how these patterns evolve and what are their functional role, enables us to understand the evolution of biocomplexity (see e.g.1–4). Morphogenesis has been studied in length at the cell level (see e.g.5–8); cells are highly sensitive to geometrical and mechanical constraints from their microenvironment and respond to these conditions by changing shape (see e.g.9,10); these transformations impact cell migration and growth (see e.g.7,11–13). Cellular migration is a fundamental property of every cell and it is crucial for the development and morphogenesis of animal body plans and organ systems (see e.g.14–16). Cell migration is either in a random direction or directed towards localized cues17–20. Mechanisms of cellular movement have been mostly studied in chemotactic cells, such as neutrophils17, bacteria21, Ciliata22, fungi23 and cellular slime molds19.
Due to its extremely fast migration rate and highly irregular shape, the acellular slime mold Physarum Polycephalum represents a prime example of differentiated growth and thus offers an attractive model for the analysis of morphogenesis dynamics underlying cellular migration and exploration24–28. P polycephalum is a giant single-celled organism that can grow to cover several square meters. Its morphology includes search fronts that are connected to a system of intersecting veins, in which oscillatory flows of the protoplasm “shuttle streaming” take place. This vein network allows (1) an efficient distribution of chemical signals, oxygen, nutrients over large distances and (2) cell migration at a speed of few centimeters per hour29,30. The driving force for this protoplasm streaming is a periodic, peristaltic contraction and relaxation of the veins due to the actin-myosin interaction, which is regulated by oscillations of intracellular chemicals such as calcium31–33. As it explores its environment, the slime mold extends temporary arm-like projections named pseudopods. It also secretes continuously a thick extracellular slime34. The glycoprotein nature of the extracellular slime coat endows P polycephalum with unique protective and structural properties that favor survival of the migrating, naked slime mold35. As the slime mold is foraging, it avoids areas covered with this mucus, which marks previously explored areas36,37.
In the presence of chemical substances in the environment, P. Polycephalum shows directional movements towards or away from the stimuli (i.e. chemotaxis). Physarum morphology, evolution and behaviors are strongly affected by the availability, location and concentration of nutrients. When the slime mold senses attractants (e.g. food cues) using specific receptors located on the membrane, the oscillation frequency in the pseudopod closest to the attractant increases, causing cytoplasm to flow towards the attractant38. On the contrary, when repellents such as salts are sensed, the oscillation frequency decreases and the slime mold moves away from the repellent38. Although slime molds lack the complex hardware of animals with brains, they live in environments that are as complex and they face the same decision-making challenges39. Hence, acellular slime molds have been the subject of a wide range of studies showing that they can solve complex biological and computational problems without any specialized nervous tissue24,36,37,40–46.
In this paper the objectives are to characterize the morphology and dynamics of Physarum exploring various environments. First, we investigate how movement is affected by homogeneous environmental conditions: adverse environment (using salt as a repellent47; nutritive environments (using glucose as a chemo-attractant48,49 with 2 different concentrations) and a neutral environment (using plain agar). Second, we analyze the geometrical evolution of slime molds placed at a distance from a nutritive spot (glucose), with and without a repelling spot (salt) in between. We characterize slime molds’ movement both temporally and spatially, to capture the full dynamics. To this aim, we develop a program that automatically analyzes sequences of images to track the areas covered and explored by the slime mold, the slime mold shape, the refinement and secondary growth cycles, as well as the distance to the nutritive spot.
Results
Homogeneous environment
In order to study the influence of the environment on slime mold expansion rate, we analyzed the areas covered by slime mold, unexplored substrate and mucus over time, as shown in Fig. 1.
In a neutral (control) and slightly nutritive environment (glucose at 100mM), the slime molds started to spread from the very beginning of the experiment (Fig. 1; Table 1 and Fig. 5 in S1 Appendix, P>0.05). By contrast, in a highly nutritive environment (200mM glucose), slime molds started to explore later (P<0.001 when compared to the control). Slime molds placed in an adverse environment (100mM NaCl) were lagging the most and only started exploring after 3hours (P<0.001 when compared to the control). Once the slime molds started to explore, they all grew at the same rate (Fig. 1; Table 2, Fig. 6 in S1 Appendix, P>0.05 when compared to the control) except the ones placed in a highly nutritive environment which were slowed down (P<0.001 when compared to the control).
At the end of the experiment (after 35hours), the slime molds reached a similar surface area in a control environment and in an adverse environment (Table 3 and Fig. 7 in S1 Appendix, P>0.05 when compared to the control). Interestingly, after reaching a plateau at 18hours, the area covered by the slime molds in an adverse environment oscillated with seemingly cyclic fluctuations (Fig. 1). In both a slightly and a highly nutritive environment, the slime molds reached a higher final surface area than the slime molds placed in a control environment (P<0.001 in both comparisons) and covered approximately 30% of arena at the end experiment. It is worth noting that in a highly nutritive environment, the area of the surface covered by the slime molds never reached a plateau after 35hours, suggesting that the slime molds did not reach its maximum surface area (Fig. 1).
Refinement i.e. appearance of mucus, was observed after 5hours in the control environment. In all other environments, mucus appeared later (Table 4 and Fig. 8 in S1 Appendix: P<0.001 for all treatments when compared to the control). In a highly nutritive environment, mucus was only observed after 10hours, which marked the strongest delay in the refinement process. Once the mucus started to be apparent, its surface area grew quicker in the control environment than in the other three treatments (Table 5 and Fig. 9 in S1 Appendix; P<0.001 for all treatments when compared to the control). Thus the area of the surface covered by mucus at the end of the experiment was the largest in the control environment where it reached 75% of the arena against 55%, 40% and 35% for the slightly nutritive, the adverse and the highly nutritive environments respectively (Table 6 and Fig. 10 in S1 Appendix; P<0.001 for all treatments when compared to the control).
Hence, slime molds placed in a control environment explored almost all the arena leaving only 5% of the arena unexplored while in the other treatments the area of the surface unexplored were significant: 15%, 35% and 38% for the slightly nutritive, the highly nutritive and the adverse environments respectively. Interestingly, although the growth rate dynamics differed between highly nutritive and adverse environments, the final unexplored surface areas were similar. In a highly nutritive environment the slime molds grew slowly and steadily while in an adverse environment slime molds grew rather quickly but after a long delay.
Next, we analyzed the evolution of the cumulative areas covered by primary growth, refinement and secondary growth (Fig. 2). The cumulative area covered by secondary growth, which reveals the cyclic nature of the exploration process, was the highest in the adverse environment (480% coverage) followed by the control environment (380%), the slightly nutritive environment (250%) and the highly nutritive environment (180%). All comparisons lead to significant differences P<0.05, except control vs. adverse environment (Table 7 and Fig. 11 in S1 Appendix). This observation confirms that exploration was slowed down by the presence of nutrients, and that the pulsatile behavior (i.e the exploration of previously explored area) was stimulated by repellents.
In accordance with the previous results, Fig. 3 (average migration rate for the treatments) shows a higher migration rate for the control treatment than for the other treatments (Table 8 and Fig. 12 in S1 Appendix: P<0.001 for each pairwise comparison). While slime molds exploring the highly nutritive environment were slower than slime molds exploring the slightly nutritive or the adverse environment (P<0.001 each), these two showed no significant differences (P>0.05).
The slime molds exploring the adverse environment showed the highest probability to explore a previously explored substrate than the other treatments as shown in Fig. 4 and Supplementary Materials (Tendency for secondary growth: Table 9 and Fig. 13 in S1 Appendix: P<0.001 for each pairwise comparison). When exploring a highly nutritive environment, slime molds also displayed a significant positive tendency for secondary growth (P<0.001 for each pairwise comparison) but significantly less strong than on the adverse environment (P<0.001). For the other treatments, the measured proportion of secondary growth was not different from the expected proportion of secondary growth, indicating that slime molds did not avoid previously explored substrate and explored randomly (Fig. 13 in S1 Appendix). The peaks observed within the first 5hours of the experiment correspond to an isotropic extension immediately followed by a refinement process that occurred before the slime mold started to explore continuously its environment. This behavior is often observed when a slime mold is introduced in a new environment and is referred as “contemplative50” i.e. the slime mold migrates, retracts and moves again. The peak was larger in an adverse environment.
We then analyzed the evolution of the shape of the slime molds contour. Note that the experimental set ups in which slime molds were placed exhibited radial symmetry. Hence, no preferential expansion direction was expected. We thus focused on contour shape, not orientation.
As expected, circularity was initially one in all tests (circular slime mold spot), and increased over time as the contour shape departed from a circle (Fig. 2 in S1 Appendix). In the control and nutritive environments, circularity remained between 1.05 and 1.10, whereas it fluctuated between 1.05 and 1.30 in the adverse environment. This observation suggests that, in an adverse environment, slime molds explored the petri dish by spreading and thinning over larger areas than in the other environments, which led to shape changes and a decrease of slime mold circularity. However fluctuations among the 20 replicates were too high to identify any trend in the evolution of slime mold circularity.
The eccentricity index was initially close to zero (circular cell) and increased up to almost 0.8 over time, with important fluctuations in all the treatments (Fig. 3 in S1 Appendix). Eccentricity is an indicator of the number of pseudopodia. But a non-eccentric convex hull can enclose non-circular contours of slime mold, since pseudopodia can develop in a symmetric fashion. That is why no major difference was noted between the treatments. This result highlights the absence of preferred expansion direction in symmetric, homogeneous environments.
Solidity decreased with the emergence of pseudopodia, since slime mold branching disrupted the initially convex shape of the slime mold (Fig. 5). In the control environment, in which the exploration rate was the highest, the decrease of solidity of the slime mold area was the highest (and the fastest) decreasing from 1 to 0.3 and then becoming relatively stable, with fluctuations of +/−0.05 (Table 10 and Fig. 14 in S1 Appendix; P<0.01 for all paired comparisons except adverse environment vs. slightly nutritive environment, where P>0.05). Slower exploration resulted in a slower and steadier loss of solidity as observed in the nutritive and adverse environments. The highly nutritive environment yielded the highest solidity at the end of the experiment (0.6), which confirmed that the presence of glucose slowed down exploration.
We next focused on the number of clusters, corresponding to the number of pseudopodia (Fig. 6). Initially the slime molds stretched as a single cluster. Once mucus started to be apparent, slime molds usually divided up into several clusters, and started the active exploration phase. In highly nutritive and adverse environments, the number of clusters over time was lower than in the other two treatments (new pseudopod number: Table 11 and Fig. 15 in S1 Appendix, and P<0.01 for all paired comparisons except for adverse environment vs. highly nutritive environment). This observation confirmed that the presence of concentrated nutrients slowed down the exploration, and that the presence of repellents triggered a highly pulsatile behavior with small exploration fronts, which were sometimes not detected as separate clusters.
Spot experiments
In the spot experiments, we studied the influence of discrete distributions of nutrients and repellents on exploration dynamics. When looking at the evolution of slime mold, mucus and unexplored substrate over time (Fig. 7), we only observed marginal difference among the treatments, which all exhibited similar patterns of exploration, e.g. similar percentage of non-explored area and similar mucus accumulation. The presence of an adverse spot only delayed the appearance of the first pseudopod (first movement: Table 12 and Fig. 15 in S1 Appendix, P<0.05) but not the first appearance of mucus (first appearance of mucus: Table 15 in S1 Appendix, not significant). The only noticeable differences lie in the surface area reached at the end of the experiment: slime molds that were offered a highly nutritive spot grew larger (final surface area: Table 14 and Fig. 17 in S1 Appendix; P<0.01). By contrast, the area of the surface covered by mucus was lower (mucus final surface: Table 17 and Fig. 18 in S1 Appendix; P<0.01) than slime molds that were offered a slightly nutritive spot. In comparison with the experiments conducted in homogeneous environments, we did not observe any expansion/refinement cycles in the spot experiments, meaning that slime mold spread steadily towards the food source despite the presence of an obstacle on the way.
The exploration behavior in the spot experiments was similar to that observed in the control environment in homogeneous experiments as observed in Fig. 8, which includes the average percentage of unexplored area for the homogeneous and spot experiments shown in Figs 1 and and77 respectively. This suggests that the spatial exploration of slime mold depended mostly on the substrate and not on the geometric distribution of the nutritive and adverse stimuli.
The cumulative areas covered by secondary growth for the spot experiments (Fig. 9) were also similar for all treatments (Table 18 in S1 Appendix, P>0.05), suggesting again, that isolated spots with nutritive or adverse stimuli did not alter the overall exploration of slime molds when growing on the same, control, substrate.
The trends of migration rate, as shown in Fig. 10, show that slime molds were not affected by the difference between treatments, showing only a slight effect of the concentration of the food spot (P<0.05), as shown on Table 19 in S1 Appendix. This effect showed that the migration rate was slightly superior when slime molds were offered a higher than a lower nutritive spot.
Similarly, looking at the predilection of slime molds to grow towards mucus (secondary growth), as shown in Fig. 11, no significant differences were observed between treatments (Table 19 in S1 Appendix), which suggests that the growth type is not influenced by the existence or concentration of discrete attracting/repelling spots.
The results obtained for the four different shape indexes (Fig. 4 in S1 Appendix) for the spot experiments support the hypothesis that discrete spots of nutrients or repellents did not affect the overall expansion dynamics and exploration cycles. This interpretation is confirmed by the average number of pseudopods, which was found to be correlated to the formation of mucus during the exploration phase: in all the spot experiments, the number of clusters increases from 1 to 2.5 within around 12hours, to reach a plateau afterwards. In other words, less exploration cycles were observed in non-homogenous environments, yielding less pseudopodia.
The evolution of the shortest distance from the slime mold cell to the glucose spot is shown in Fig. 12, which can be viewed as a “survival” plot, displaying the proportion of the replicate (P) in which slime mold has not reached the nutritive spot at a given time. Both increasing the concentration of nutrient in the spot (from slightly nutritive to highly nutritive) and adding an aversive spot increased the time to reach the food patch by decreasing the probability to reach it (time to food patch: Table 22 and Fig. 18 in S1 Appendix, nutritive: P<0.01; aversive: P<0.05).
Discussion
Exploration in slime molds involves two different processes: area extension and movement27. Slime molds locomotion and morphogenesis depend on the response of the organism to the environmental conditions such as light, hygrometry, pH, or the presence of chemicals29,51. In the present paper, we show that distributions of nutritive and aversive cues affected drastically the exploration pattern of slime molds.
The typical exploration behavior of a slime mold (in control conditions) started by a stretching period where the slime molds grew uniformly in all directions for 10hours. Then, the contour of slime mold lost circularity when the first pseudopodia appeared, which also corresponds to the first occurrence of mucus. This phenomenon is typical of the directed digitated growth, or branching phase, described by numerous authors27,52–55. Slime molds developed multiple pseudopodia and did not exhibit any preferential exploration orientation. At the end of the experiment, almost all the arena was explored by the slime mold. Thus the exploration was characterized by three phases: (i) Primary growth only, in the quasi-absence of mucus (5–10hours); (ii) Combination of primary and secondary growth; (iii) Secondary growth only, when the slime mold stops exploring new substrate areas.
Regarding the result obtained for the homogeneous distribution of nutritive cues, we noted first, an environment containing a uniform distribution of nutrients slowed down the exploration of slime mold, mainly by delaying secondary growth and increasing the period of the pulsatile exploration/refinement movements. The area not explored by the slime molds was 3 times larger (respectively 7 times larger) than in the environment deprived of nutrient (control case) for a slightly nutritive (respectively highly nutritive) environment. The exploration rate was almost linear for highly nutritive environments, while for other treatments, the area covered by the slime molds reached a plateau after a period of stretching, which indicates secondary growth and slime mold displacement. This means that slime molds that explored nutritive environments never exhibited a Phase (iii) in their exploration pattern.
Second, on substrates with higher nutrient concentrations, the slime molds grew in a more compact fashion, i.e. slime molds presented the highest solidity index and the lowest number of pseudopodia (clusters). Additionally, the appearance of mucus, which indicates that the slime mold was withdrawing, occurred much later in nutritive environments. As glucose is only aversive when only above 300mM54, our results suggest that nutritive media depressed migration due to feeding behavior. This allows the organism to remain at a site until nutrients are exhausted56–58. In previous studies, it was shown that the surface area of substrate covered increases when slime mold responds to nutrient dilution58,59. Here, we confirmed these observations and noted that slime mold tended to migrate and grow faster on substrates with the lowest concentration of nutrients, thus maximizing nutrient intake and optimizing the trade-off between nutrient foraging and nutrient intake.
Similarly, for a homogeneous distribution of aversive cues we noted the aversion of slime mold to salt manifested itself through longer contemplative behavior, delayed primary growth and a higher probability to crawl on previously explored surface. In addition, the slime molds grew more compactly and with less pseudopods. This suggests that slime molds were actively avoiding contact with the aversive surface44,50,60. In the absence of cell walls, slime mold has no other protection from the environment than mucus. Indeed, in bacteria for instance, mucus is used as protective barrier for the cells against harsh external conditions61. In slime molds, the extracellular mucus excreted by the slime molds can have different roles: hydrophilic shield to prevent water loss62, a lubricating surface over which the slime molds can easily crawl63, a defensive coat to protect against invasion by foreign materials and organisms62,64, an aid to phagocytosis65, a surface that promotes ion-exchange62 and has externalized spatial memory that helps navigation in unknown environments26,36,66. Here, we can add a new function for the mucus i.e. a buffer to move in adverse environment. In other words we demonstrate for the first time that the mucus may be used as a safe haven to avoid prolonged and repeated contacts with hazardous substances encountered in the environment.
Finally, for non-homogeneous distribution of aversive and non-aversive cues, e.g. spot experiments, pulsatile movements were limited and slime molds responded in the same way regardless of the concentration of glucose used as attractant and regardless of the presence of a salt spot on the way to the glucose spot. Slime molds grew in a more compact fashion, i.e. slime molds presented the highest solidity index and the lowest number of pseudopodia (clusters). Additionally, the evolution of the areas covered by slime mold and mucus over time indicates that the response of slime molds to heterogeneous environments was similar to that in the control case. This result suggests that in the spot experiments, the exploration behavior of slime mold is mostly controlled by the substrate. Our observation confirms that salt reception can be affected by the presence of sugars47. The authors in47 showed that the “apparent” enthalpy change accompanying salt perception decreases with increase of sugar concentration.
The proposed image analysis program allows extracting information on expansion rates, geometric changes and probability of occupancy. Ongoing developments aim to acquire high quality images of slime mold exploration tests and to expand the code’s capabilities for extracting topological information on the networks formed by slime molds. Slime molds have been shown to be capable of distributed sensing, parallel information processing, and decentralized optimization39,51,67. It is also viewed by some researchers as an inspiration for bio-computing devices68–73. Our image analysis methods will help such research avenues by improving data collection.
Trajectories of individuals within a flock could also be described to understand whether or not members of the flock can inform and influence the travel direction of other individuals. Such a finding would allow understanding how group decisions are made among gregarious species74. The cluster identification and shape recognition program could be used to differentiate modes of gene expression or to recognize objects75. Object identification is an important pillar to explain associative memory or to track species interactions in an ecosystem. The tools and approach presented here are thus applicable to any problem of network dynamics or pattern recognition76.
Methods
Species
Physarum polycephalum, also known as the true slime mold, belongs to the Amoebozoa, the sister group to fungi and animals77. Slime molds are found on organic substrates where they feed on microorganisms such as bacteria or fungi77. The vegetative morph of P. polycephalum, the plasmodium, is a vast multinucleate cell that can grow to cover up to a few square meters and crawl at speeds from 0.1 to few centimeters per hour29,30. When hygrometry and food availability decrease, the plasmodium turns into an encysted resting stage made of desiccated spherules called sclerotium29.
Rearing conditions
Experiments were initiated with a total of 10 sclerotia per strain (Southern Biological, Victoria, Australia). We cultivated slime molds on a 1% agar medium with rolled oat flakes, slime molds were fed every day and the medium was replaced daily. Slime molds were 2 weeks old when the experiment started. All experiments were carried out in the dark at 25°C temperature and 70% humidity, and ran for 35h. Pictures were taken with a Canon 70D digital camera.
Experimental setup
Initially we monitored the exploration movement evoked in slime molds in a homogeneous environment. Each slime mold was placed in the center of a circular arena (14.5cm in diameter) with a layer of agar (1% in water) mixed with non-nutritive cellulose (5%). Adding cellulose to the agar mix proved to be useful to obtain a homogeneous pigmentation and to enhance the color contrast between the substrate and slime mold, therefore improving the identification process. A circular hole (2.5cm in diameter) was punched and replaced with a circular slime mold of the same size sitting on oat. In the first and second treatments (nutritive environments) we added glucose (100mM or 200mM) to the medium. In the third treatment (adverse environment), we added a known repellent (NaCl 100mM55) to the medium. Lastly, in the fourth treatment, the medium remained unchanged (neutral environment i.e. control treatment).
Subsequently, to investigate how chemotaxis modified the exploration behavior, we introduced discrete spots of attractants/repellants within a neutral substrate made of plain agar. In these so-called “spot experiments”, we followed a procedure similar to that for the homogeneous environments. The arena consisted of a 14.5cm diameter petri dish filled with plain 1% agar mixed with non-nutritive cellulose (5%). Once the agar had set, we punched two or three holes in a line configuration (diameter of each hole 2.5cm). A circular slime mold (2.5cm in diameter) was placed diametrically opposite to a glucose (attractant) spot of same size, placed 4.5cm away. In some of the treatments, a salt (NaCl 200mM, repellent) spot was added between the slime mold and the attractant spot. A total of 4 different treatments were tested: the first and the second with a single spot of glucose at concentrations of 100mM and 200mM respectively, the third and fourth keeping the glucose spots with the same concentrations and adding a NaCl 200mM spot. Previous experiments show no disruption of slime mold behavior due top punching alone45.
All slime molds were fed just before the experiment so we assumed that they were in the same physiological state. The experiment consisted of a total of 8 different treatments. We replicated the experiment 20 times for each treatment and monitored each arena for 35hours taking time-lapse photographs every 5minutes. Figure 13 shows the experimental set-up for homogeneous environments (left) and discrete distributions of attractants/repellants (right).
Image processing
Time-lapse images were taken every 5minutes for a total of 420 pictures for each replicate. First, the outside of the petri dish was masked by manually identifying the circular contour; then, the images were converted into the ab* color space (which is the CLAB space without the L* lighting component) since this color space yielded better segmentation results. Then, the clustering algorithm k-means78,79 was used to classify every pixel into one of two categories or clusters: slime mold or not-slime mold. This last category was further refined based on pixel history, becoming either unexplored surface or mucus. Mucus is the substance left by slime after refinement, which acts as an external memory on explored areas36,66.
The change of class from unexplored substrate to slime mold, defined as primary growth, means that slime mold reached a point it had never explored. Similarly, secondary growth is the change from mucus to slime mold, meaning that the cell is revisiting a location. Lastly, if the slime mold recedes, e.g. a pixel goes from slime mold to non-slime mold, it becomes mucus, and is defined as refinement. Note that once a pixel is identified as mucus, it can never be classified as unexplored substrate in the following time frames. Two videos are provided as Supplementary Material, S2 and S3 show one replicate as original images and trinarized images respectively.
Image analysis in space and time
We computed indexes to characterize slime mold geometry dynamics and averaged them over the 20 replicates to obtain statistically representative measures. In order to quantify the differences on distinct substrates, we first calculated the fraction of the petri dish area covered by slime mold, mucus and unexplored substrate over time. The total area, the lighting conditions and the test duration were the same for all treatments, both in the homogeneous and spot experiments. Note that glucose only provides energy to slime mold, which is not gaining significant mass during the experiments58. In other words, slime mold is changing its area by mostly by stretching and contracting, therefore changing its area density.
We then computed the cumulative area of primary growth, refinement, and secondary growth over the full period of the experiments comparing two consecutive images at the time. The cumulative area covered by primary growth is indicative of the total area of exploration, therefore it is always smaller or equal to the total area of the dish. The cumulative area covered by secondary growth indicates whether slime mold expansion is monotonic (dominated by primary growth) which results in a smaller magnitude, or cyclic (secondary growth dominated, with pulsatile movements) which results in a larger magnitude. The cumulative area covered by refinement indicates slime mold density changes. Within a given time interval, if the area covered by primary plus secondary growth equals that covered by refinement, then slime mold keeps the same density, whereas if it is superior, the slime mold stretches (e.g. density decreases). If secondary growth is negligible and if the area covered by primary growth equals the area covered by refinement, then slime mold displaces mass.
We calculated the extent of growth for each pair of consecutive images as the distance from each pixel where growth occurred (both primary and secondary) to the closest pixel classified as slime mold in the previous image. We calculated the migration rate as the ratio between the maximum extent of growth and the time interval between images (5min). Then we delineated the region explored by slime mold within the interval as the contour of the slime mold with an offset distance corresponding to the migration rate (see the Supplementary Material Fig. 1 in S1 Appendix for more details).
We estimated the fraction of secondary growth relative to the total number of pixels in the region of expansion. We then calculated the fraction of “expected secondary growth”, which would have occurred if secondary growth had happened randomly. If the measured secondary growth fraction is higher (respectively, lower) than the expected one, this means that slime mold has a bias towards mucus (respectively, unexplored substrate).
Additionally, we computed four shape parameters indicative of the contour of slime mold: circularity, eccentricity, solidity and number of clusters. Circularity (C) is defined as:
where P and A are the perimeter and area of the shape of slime mold at a given time. Eccentricity (E) is calculated as the ratio between the distance between the foci and the major axis length, as follows:
In which a and b are the lengths of the major and minor axes, respectively. Solidity (S) is the ratio between the area of the slime mold contour and the area of its convex hull. Lastly, we measured the number of clusters by performing an erosion operation along the contour of slime mold, after which only the clusters of high concentration of slime mold remain, which provided the number of pseudopodia or clusters of slime mold.
For the spot experiments, we also determined the distance from the slime mold to the glucose spot at every time, as the minimum distance between the contour of slime mold and the glucose spot. The evolution of the distance to glucose over time was analyzed in a way similar to a survival analysis.
Statistics
The full description of the statistics is provided as part of the Supplementary Information; Appendix S1 includes the statistical procedure and results; while Appendix S4 is an R markdown allowing analysis reproduction. When dependent variables lasted until the occurrence of certain event, we conducted survival analyses using the R package coxme80. For the remaining dependent variables, we did linear analyses using the R packages lme481 and lmerTest82. We tested as fixed factors the four treatments in the homogeneous experiment and both the attracting and the repelling spots and their interaction. The date of the experiment was considered as a random factor. Finally, we performed a nested model comparison using the R package MuMIn83 by ordering models according to their Akaike criterion and represented the selected model by plotting each estimator, their 95%CI and p-values.
Author summary
Physarum polycephalum, also called slime mold, is a giant single-celled organism that can grow to cover several square meters, forming search fronts that are connected to a system of intersecting veins. An original experimental protocol allowed tracking the shape of slime mold placed in homogenous substrates containing an attractant (glucose) or a repellent (salt), or in homogeneous substrates that contained an attractive spot (glucose), an eccentric slime mold and a repulsive spot (salt) in between. For the first time, the rate of exploration of unexplored areas (primary growth) and the rate of extension in previously explored areas (secondary growth) were rigorously measured, by means of a sophisticated image analysis program. This paper shows that the chemical composition of the substrate has more influence on the morphology and growth dynamics of slime mold than that of concentrated spots of chemicals. It was also found that on a repulsive substrate, slime mold exhibits a bias towards secondary growth, which suggests that the mucus produced during slime mold migration acts as a protective shell in adverse environments.
Supplementary information
Acknowledgements
This work was supported by the U.S National Science Foundation, under grant CMMI#1552368: “CAREER: Multiphysics Damage and Healing of Rocks for Performance Enhancement of Geo-Storage Systems - A Bottom-Up Research and Education Approach.” A.D. and A.B. were supported by a grant from the Agence Nationale de la Recherche (Reference Number: ANR-17-CE02-0019-01 -SMART- CELL).
Author contributions
Author Fernando Patino-Ramirez (F.P.) carried out the image analysis, wrote the corresponding sections about methods and results created the figures in the manuscript and was responsible of the image analysis section on S1 and the videos S2 and S3. A.B. developed the statistical analysis, including the corresponding sections in Supplementary Materials S1 and S4. Authors Aurèle Boussard (A.B.) and Audrey Dussutour (A.D.) carried out the experiments in the laboratory and the corresponding image acquisition. Chloé Arson (C.A.) drafted the manuscript and contributed in the planning, analysis and interpretation of the experiments and the data together with A.D., F.P. and A.B. All the authors have approved the submitted version of this manuscript and have agreed both to be accountable for the author’s own contributions.
Competing interests
The authors declare no competing interests.
Footnotes
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
These authors contributed equally: Chloé Arson and Audrey Dussutour.
Contributor Information
Fernando Patino-Ramirez, Email: ude.hcetag@pf.
Audrey Dussutour, Email: rf.3eslt-vinu@ruotussud.yerdua.
Supplementary information
is available for this paper at 10.1038/s41598-019-50872-z.
References
Articles from Scientific Reports are provided here courtesy of Nature Publishing Group
Full text links
Read article at publisher's site: https://doi.org/10.1038/s41598-019-50872-z
Read article for free, from open access legal sources, via Unpaywall: https://www.nature.com/articles/s41598-019-50872-z.pdf
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
Smart citations by scite.ai
Explore citation contexts and check if this article has been
supported or disputed.
https://scite.ai/reports/10.1038/s41598-019-50872-z
Article citations
Hydrophobic Barriers for Directing Physarum polycephalum Propulsion and Navigation.
ACS Omega, 8(44):41649-41654, 25 Oct 2023
Cited by: 0 articles | PMID: 37970039 | PMCID: PMC10634242
Behavioural changes in slime moulds over time.
Philos Trans R Soc Lond B Biol Sci, 378(1874):20220063, 20 Feb 2023
Cited by: 1 article | PMID: 36802777 | PMCID: PMC9939273
An effective solution to numerical and multi-disciplinary design optimization problems using chaotic slime mold algorithm.
Eng Comput, 38(suppl 4):2739-2777, 30 May 2021
Cited by: 6 articles | PMID: 34092833 | PMCID: PMC8164690
Substrate and cell fusion influence on slime mold network dynamics.
Sci Rep, 11(1):1498, 15 Jan 2021
Cited by: 4 articles | PMID: 33452314 | PMCID: PMC7810851
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
BioStudies: supplemental material and supporting data
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Substrate and cell fusion influence on slime mold network dynamics.
Sci Rep, 11(1):1498, 15 Jan 2021
Cited by: 4 articles | PMID: 33452314 | PMCID: PMC7810851
Slime molds response to carbon nanotubes exposure: from internalization to behavior.
Nanotoxicology, 15(4):511-526, 11 Mar 2021
Cited by: 0 articles | PMID: 33705250
Toxicity of heavy metals and insecticides on slime mold Physarum polycephalum.
J Toxicol Sci, 3(4):293-303, 01 Nov 1978
Cited by: 0 articles | PMID: 739566
Intelligent behaviors of amoeboid movement based on complex dynamics of soft matter.
Soft Matter, 4(1):57-67, 01 Dec 2007
Cited by: 7 articles | PMID: 32907084
Review
Funding
Funders who supported this work.
Agence Nationale de la Recherche (1)
Grant ID: ANR-17-CE02-0019-01 -SMART- 644 CELL
National Science Foundation (1)
Grant ID: CMMI#1552368