Substrate composition directs slime molds behavior

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.


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 100 mM), 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 (200 mM glucose), slime molds started to explore later (P < 0.001 when Figure 1. Fraction of area covered by slime mold, mucus and unexplored substrate -Homogeneous environment. The solid lines correspond to the average index calculated over the 20 replicates, while the shaded areas correspond to the first and third quartiles of the data, the dashed line correspond to the ratio between the average mucus area over the average slime mold area over time. compared to the control). Slime molds placed in an adverse environment (100 mM NaCl) were lagging the most and only started exploring after 3 hours (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 35 hours), 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 18 hours, 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 35 hours, suggesting that the slime molds did not reach its maximum surface area (Fig. 1).
Refinement i.e. appearance of mucus, was observed after 5 hours 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 10 hours, 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 Figure 2. Cumulative areas covered by primary growth, refinement and secondary growth -Homogeneous experiments. The solid line corresponds to the average index calculated over the 20 replicates, while the shaded areas correspond to the first and third quartiles of the data. 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 5 hours 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 "contemplative 50 " 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 www.nature.com/scientificreports www.nature.com/scientificreports/ exploration, and that the presence of repellents triggered a highly pulsatile behavior with small exploration fronts, which were sometimes not detected as separate clusters.  www.nature.com/scientificreports www.nature.com/scientificreports/ 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   www.nature.com/scientificreports www.nature.com/scientificreports/ 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 7 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 12 hours, 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 movement 27 . 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 chemicals 29,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 10 hours. 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 authors 27,52-55 . Slime molds developed multiple pseudopodia and did not exhibit any preferential exploration orientation. At the end of the 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  www.nature.com/scientificreports www.nature.com/scientificreports/ later in nutritive environments. As glucose is only aversive when only above 300 mM 54 , our results suggest that nutritive media depressed migration due to feeding behavior. This allows the organism to remain at a site until nutrients are exhausted [56][57][58] . In previous studies, it was shown that the surface area of substrate covered increases when slime mold responds to nutrient dilution 58,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 surface 44,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 conditions 61 . In slime molds, the extracellular mucus excreted by the slime molds can have different roles: hydrophilic shield to prevent water loss 62 , a lubricating surface over which the slime molds can easily crawl 63 , a defensive coat to protect against invasion by foreign materials and organisms 62,64 , an aid to phagocytosis 65 , a surface that promotes ion-exchange 62 and has externalized spatial memory that helps navigation in unknown environments 26,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. www.nature.com/scientificreports www.nature.com/scientificreports/ 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 sugars 47 . The authors in 47 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 optimization 39,51,67 . It is also viewed by some researchers as an inspiration for bio-computing devices [68][69][70][71][72][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 species 74 . The cluster identification and shape recognition program could be used to differentiate modes of gene expression or to recognize objects 75 . 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 recognition 76 .

Species.
Physarum polycephalum, also known as the true slime mold, belongs to the Amoebozoa, the sister group to fungi and animals 77 . Slime molds are found on organic substrates where they feed on microorganisms such as bacteria or fungi 77 . 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 hour 29,30 . When hygrometry and food availability decrease, the plasmodium turns into an encysted resting stage made of desiccated spherules called sclerotium 29 .
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 35 h. 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.5 cm 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.5 cm in diameter) was punched and www.nature.com/scientificreports www.nature.com/scientificreports/ replaced with a circular slime mold of the same size sitting on oat. In the first and second treatments (nutritive environments) we added glucose (100 mM or 200 mM) to the medium. In the third treatment (adverse environment), we added a known repellent (NaCl 100 mM 55 ) 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.5 cm 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.5 cm). A circular slime mold (2.5 cm in diameter) was placed diametrically opposite to a glucose (attractant) spot of same size, placed 4.5 cm away. In some of the treatments, a salt (NaCl 200 mM, 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 100 mM and 200 mM respectively, the third and fourth keeping the glucose spots with the same concentrations and adding a NaCl 200 mM spot. Previous experiments show no disruption of slime mold behavior due top punching alone 45 .
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 35 hours taking time-lapse photographs every 5 minutes. 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 5 minutes 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-means 78,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 areas 36,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 experiments 58 . 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 (5 min). 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). www.nature.com/scientificreports www.nature.com/scientificreports/ 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 coxme 80 . For the remaining dependent variables, we did linear analyses using the R packages lme4 81 and lmerTest 82 . 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 MuMIn 83 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.