Guppies Prefer to Follow Large (Robot) Leaders Irrespective of Own Size

Body size is often assumed to determine how successfully an individual can lead others with larger individuals being better leaders than smaller ones. But even if larger individuals are more readily followed, body size often correlates with specific behavioral patterns and it is thus unclear whether larger individuals are more often followed than smaller ones because of their size or because they behave in a certain way. To control for behavioral differences among differentially-sized leaders, we used biomimetic robotic fish (Robofish) of different sizes. Live guppies (Poecilia reticulata) are known to interact with Robofish in a similar way as with live conspecifics. Consequently, Robofish may serve as a conspecific-like leader that provides standardized behaviors irrespective of its size. We asked whether larger Robofish leaders are preferentially followed and whether the preferences of followers depend on own body size or risk-taking behavior (“boldness”). We found that live female guppies followed larger Robofish leaders in closer proximity than smaller ones and this pattern was independent of the followers’ own body size as well as risk-taking behavior. Our study shows a “bigger is better” pattern in leadership that is independent of behavioral differences among differentially-sized leaders, followers’ own size and risk-taking behavior.


INTRODUCTION
The question of what makes an individual successful in leading others is a long-standing issue King et al., 2009). In shoaling fish, those individuals that occupy front or periphery positions within a shoal are assumed to have the greatest influence on the group's movement direction, hence are capable of leading the other shoal members (Bumann and Krause, 1993;Krause et al., 2000;Lopez et al., 2012;Jolles et al., 2017). Often, occupation of front or peripheral positions is related to motivational or phenotypical differences among individuals Hemelrijk and Kunz, 2005). For example, individuals that take up a front position are often hungrier (Krause et al., 1992;McLean et al., 2018), more risk-taking ("bolder") (Leblond and Reebs, 2006;Nakayama et al., 2012aNakayama et al., ,b, 2016 or simply larger (Krause et al., 1998;Guttridge et al., 2011) than the rest of the group. Mechanistically, those front individuals may move faster (Gueron et al., 1996;Krause et al., 1998;Jolles et al., 2017) or have larger repulsion areas (Hemelrijk and Kunz, 2005;Romenskyy et al., 2017), both resulting in an assortment within the shoal. However, being at the front (i.e., taking the lead) is often not the only factor determining leadership success. Using the golden shiner (Notemigonus crysoleucas), Reebs (2001) showed that a minority of informed large fish was capable of leading a shoal of small fish to a food source, whereas informed small fish had much lower success in leading a shoal of large fish even when occupying the front positions of the shoal. Furthermore, when sticklebacks (Gasterosteus aculeatus) were grouped with two partners of different personalities, they were more likely to follow the partner of similar personality out of refuge (Nakayama et al., 2016). Thus, both body size as well as behavior may determine leadership success in fishes. Moreover, both body size and behavior often covary with each other, for example larger fish can swim faster (Domenici, 2001) than smaller ones or exhibit a certain personality (Polverino et al., 2016). Just recently Romenskyy et al. (2017) concluded that "fish of different sizes cannot be considered simply as particles of different physical size, since their behavior changes with their size." This poses the question whether larger individuals are more often followed than smaller ones simply because they are larger or because they behave in a certain way. Furthermore, we do not know whether following behavior is influenced by the followers' own body size or behavior, or how either attribute may interact with leader size. To answer these questions, we experimentally controlled for the leader's behavior while simultaneously varying its body size through the use of differentially sized biomimetic robotic fish.
Biomimetic robots have become a recent tool to investigate animal behavior Romano et al., 2018). These machines consist of an animal-like part that is either selfpropelled or externally dragged by a robotic unit. Biomimetic robots can be either interactive (closed-loop behavior), which means that they change their behavior in response to the actions of live animals, or static (open-loop behavior), which means that they move and behave in predefined, non-interactive ways (Webb, 2000;Krause et al., 2011;Butail et al., 2015;Romano et al., 2018). Biomimetic robots thus provide the experimenter with a diverse toolset to study social interactions such as the ability to provide completely standardized social cues (e.g., through the use of non-interactive open-loop robots, see Abaid et al., 2012;Phamduy et al., 2014;Bierbach et al., 2018a). Furthermore, the robot's parameters can be set to either resemble those of focal live individuals or show a sharp contrast with them (Butail et al., 2013(Butail et al., , 2014Polverino and Porfiri, 2013). In addition, closedloop-controlled robots allow us to create interactive scenarios that nevertheless follow controlled rules that can be adapted intentionally (Kopman et al., 2013;Landgraf et al., 2013Landgraf et al., , 2014Landgraf et al., , 2016Bonnet et al., 2018;Kim et al., 2018;Datteri, 2020).
To date a wide range of taxa has been shown to accept biomimetic robots as conspecific or heterospecific animals. For example, Halloy et al. (2007) developed a robot that interacted autonomously with live cockroaches and therefore allowed finescaled investigations of their aggregation behavior. Romano et al. (2017a) investigated the lateralization of escape and surveillance responses in locusts during predator-prey interactions with a robot that resembled a predatory bird. Similarly, fish shoals were attacked by a robotic predator fish to investigate their collective predator evasion responses (Swain et al., 2012;Romano et al., 2020). Such interactions can severely impact growth rates and body conditions of the attacked fish, even when experiencing these attacks only for short amounts of time (Polverino et al., 2019). These studies exemplify the extensive use of different fish species in studies with biomimetic robots (but see Romano et al. (2018)) for a more complete list of taxa). More or less natural interaction patterns among live fish and biomimetic robots have now been reported for poeciliids Bierbach et al., 2018a,b;Heathcote et al., 2018), killifish , zebrafish (Kim et al., 2018), golden shiners (Abaid et al., 2013), mormyrids (Donati et al., 2016;Worm et al., 2018), Siamese fighting fish (Romano et al., 2017b) as well as sticklebacks (Faria et al., 2010).
Their success in being accepted as conspecific or heterospecific animals may be due to Nico Tinbergen's (1948) idea of "social releasers, " meaning that only a small subset of perceivable cues are communicative signals. Thus, even minimalistic robot models can exploit species-specific cues that identify conspecifics or heterospecifics (see Landgraf et al., 2016;Datteri, 2020 for discussion). In Halloy et al. (2007), for example, the robot was treated with a cockroach-specific pheromone to facilitate group integration. In poeciliids like the herein used guppies, replicas equipped with realistic glass eyes were found to be followed almost as close as live conspecific partners (Landgraf et al., 2016).
In the current study, we used the so-called Robofish system which is an open-loop controlled (e.g., non-interactive) robot platform that steers an exchangeable 3D-printed fish dummy. In our case, the replica resembled a live female Trinidadian guppy (Poecilia reticulata) (Landgraf et al., 2016). It was recently shown that live guppies interact similarly with Robofish as they do with a live conspecific (Bierbach et al., 2018b). Further, live guppies maintained individual differences in followership patterns exhibited during Robofish trials even when tested consecutively with a live conspecifics (Bierbach et al., 2018b).
We tested differentially-sized live female guppies for their risk-taking behavior (i.e., time to leave a shelter box) and their tendency to follow one of three differentially-sized, robotically steered replicas that all moved almost identically on a predefined trajectory in a large experimental tank. We asked (a) whether larger Robofish leaders are preferentially followed (assuming "bigger is better"), and (b) whether the following tendencies of followers depend on their own body size or their risk taking behavior ("boldness"). As in previous studies (Landgraf et al., 2016;Bierbach et al., 2018a,b), we assumed stronger following tendencies when live fish kept shorter average distances toward to moving Robofish.

Study Organism and Maintenance
We used female Trinidadian guppies (Poecilia reticulata) that were descendants of wild-caught fish from the Arima River in North Trinidad. Test fish came from large, randomly outbred single-species stocks maintained at the animal care facilities at the Faculty of Life Sciences, Humboldt University of Berlin. We provided a natural 12:12 h light:dark regime and maintained water temperature at 26 • C. Fish were fed twice daily ad libitum with commercially available flake food (TetraMin TM ) and once a week with frozen Artemia shrimps.

The Robofish System
The Robofish is a three-dimensional (3D)-printed guppy-like replica that is attached to a magnetic base. The magnetic base aligns with a wheeled robot that is driving below the actual test tank (88 × 88 cm, coated with white plastic foil) on a transparent second level. Hence, the replica can be moved directly by the robot (Figure 1). The entire system is enclosed in a black, opaque canvas to minimize exposure to external disturbances. The tank is illuminated from above with diffused LED lights. On the floor, a camera is facing upwards to track the robot's movements through the transparent second level. A second camera is fixed above the tank to track both live fish and replicas. Two computers are used for system operation: one PC tracks the robot, computes and sends motion commands to the robot over a wireless channel. The second PC records the video feed of the second camera which is afterward tracked by custom-made software (Mönck et al., 2018). Please see our Supplementary Material for more details on the Robofish system as well as (Landgraf et al., 2016).

Experimental Setup
To provide live female guppies with differently sized Robofish leaders, we used three replicas that differed only in body size [replica 1 = 20 mm standard length (SL); replica 2 = 25 mm SL, replica 3 = 30 mm SL, see Supplementary Figure S1]. All replicas were equipped with 3 mm glass eyes. Thus, the relative eye size decreased from the smallest to the largest replica (0.15-0.10). Such a decrease in relative eye size is common to all vertebrates including fishes (see Richardson et al., 2015) and thus our replicas mirrored the naturally found change in relative eye size during ontogeny. As we used transparent screws to attach the replica to its magnet foot, all replicas regardless of size kept the same distance to the water surface (1 cm, at 10 cm water level). We used only females to avoid sex differences in attraction toward Robofish, which resembles a female guppy. Test fish were randomly chosen to span the natural body size variation of this species (ranging from 18.0 to 32.0 mm, mean = 25.6 mm, SD = 4.1 mm, n = 88 fish used in this study). To measure body size, fish were transferred into a water-filled petri dish placed upon millimeter paper after the behavioral testing. We took a picture from centrally above using a SLR camera (Canon EOS 400D) and measured standard length of the fish from these pictures using ImageJ software (Schindelin et al., 2012). We tested 90 fish but had to remove 2 fish from the analysis as the video recording failed due to technical issues.
To initiate a trial, we transferred individual test fish into an opaque PVC cylinder located at the lower left corner of the test tank. The PVC cylinder had an opening (diameter 3 cm), which was closed with a sponge. We removed the sponge after 1 min of acclimation and noted the time each fish took to leave the cylinder as a proxy for its risk-taking tendency (i.e., "boldness"), which is thought to correlate with following tendencies (Nakayama et al., 2012a(Nakayama et al., ,b, 2016Jolles et al., 2015). We initiated the Robofish's movement sequence when the live fish left the cylinder (i.e., one body length away from the cylinder's border). Robofish moved along a zig-zag pattern with a maximum speed of 15 cm/s and reduced its speed at the turning points to almost 0 cm/s, before accelerating again to the predefined maximum speed. This stop-and-go motion pattern led to an average speed of 10 cm/s. A zigzag movement was found to increase the likelihood of the Robofish to be followed (Landgraf et al., 2016) and the differently-sized replicas did not differ in exhibited velocities (see Supplementary Material). During the trials, Robofish moved to the opposite corner and then counter-clockwise to its start position. This round was repeated for a second time and a trial took about 60 s in total (see Figure 2 for an example track as well as Supplementary Video S1). Each trial was videotaped for subsequent tracking and the test fish was transferred back to its holding tank after size measurement was completed. Videos were recorded at 30 fps and also tracking was performed at the same sampling frame rate via Biotracker (Mönck et al., 2018). We analyzed the first 50 s after the fish left the shelter box, resulting in 1,500 frames analyzed per trial. We calculated the inter-individual distance between focal fish and Robofish as the average distance between subjects for all 1,500 frames. IID has been shown to reflect a live fish's tendency to follow the moving Robofish (Landgraf et al., 2016;Bierbach et al., 2018a,b).

Statistical Analysis
We initially log10 transformed both recorded continuous response variables (IID, time to leave start box) to match a Gaussian distribution. We then used the IID as dependent variable in an ANCOVA (unianova package in SPSS 25) with "leader size" as a fixed factor. We included "live fish body size (SL in mm)" and "time to leave start box [log10(s)]" as well as their interactions with "leader size" as covariates in the model. Non-significant interaction terms were removed from the final model. In order to test whether differently-sized live fish differ in their risk aversion tendency, we further correlated live fish body size (SL in mm) with time to leave shelter [log10 (s)] using Pearson's correlation.

RESULTS
Regardless of own size [non-significant covariate "live fish body size" F (1, 83) = 3.36; p = 0.071], live guppies tested with larger Robofish replicas followed significantly closer than those tested with smaller ones [significant effect of factor "leader size" F (2, 78) = 4.49; p = 0.009, Figure 2]. There was no size assortative pattern detectable (i.e., smaller live fish did not follow smaller replicas closer than larger ones and vice versa) as suggested by a non-significant interaction term "leader size × live fish body size" [F (4, 79) = 0.49; p = 0.62]. Also, the time each fish took to leave the start box had no significant influence on its following behavior [F (1, 83) = 0.98; p = 0.33] and there was no significant interaction with leader size detectable [F (2, 79) = 1.51; p = 0.23]. We found no significant correlation between live fish's body size and their tendency to leave the start box (Pearsons's r = 0.190, p = 0.073).

DISCUSSION
Live guppies followed larger Robofish leaders closer than smaller ones and this pattern was independent of the followers' own body size as well as risk-taking behavior. This result is consistent with a study using golden shiners that found large individuals are more readily followed than small ones (Reebs, 2001). However, this is the first study that showed a "bigger is better" pattern in leadership in shoaling fish while controlling for the potentially confounding effects of behavioral differences (in small and large fish) by using a biomimetic robot.
Body size in fishes is often inevitably linked to specific behavioral patterns (Domenici, 2001;Polverino et al., 2016) and it is thus experimentally difficult to disentangle which cue (body size or linked behavior) is used by individuals that have to choose whom to follow among conspecifics of different sizes. While researchers from the field of sexual selection make use of video animations in binary choice tests to decouple behavior from body size and keep either one constant while varying the other (see Fisher and Rosenthal, 2007;Fisher et al., 2009;Gierszewski et al., 2017;Sommer-Trembo et al., 2017), the study of collective movement has largely relied on the use of live stimuli (but see Stowers et al., 2017 for a working Virtual Reality set-up). We addressed this issue by using a biomimetic robot toward live guppies show natural interaction patterns (Landgraf et al., 2016). Still, open-loop controlled robots (e.g., those that do not adjust their behavior in real-time to the actions of a partner) may have some short-comings in terms of reproducing natural behavioral feedbacks (see discussion in Bierbach et al., 2018b), they provide a powerful tool to present live animals with a standardized and almost identical set of social cues in a full contact design (i.e., not only visual cues available). It is thus possible to control for differences in leaders' behavior that may affect leadership success like swimming speed (Jolles et al., 2017) and movement patterns (Ioannou et al., 2015;Nakayama et al., 2016), while simultaneously varying the parameter of interest (e.g., body size of the leader).
Our results show that live guppies followed larger Robofish closer than smaller ones and there are several (although not mutually exclusive) explanations for this result: (a) Larger individuals are often focused on by predators (Pocklington and Dill, 1995) and, in case of poeciliid females, experience more harassment by males (Herdman et al., 2004;Agrillo et al., 2006). Thus, small female guppies can benefit from associating with larger ones, as doing so may take away the attention of predators and/or harassing males. Larger individuals, however, can benefit from associating with other large individuals to minimize the oddity effect during predation (Peuhkuri, 1997;Krause and Ruxton, 2002). (b) Larger individuals in species with continuous growth throughout their lives are older and thus more experienced than smaller, younger individuals. As most teleost fishes show indeterminate growth, larger and thus older ones may have accumulated experience that provides them with fitness benefits, for example in terms of predator recognition (Brown and Smith, 1998;Holmes and McCormick, 2010) and navigation (Odling-Smee and Braithwaite, 2003). It might thus be beneficial for followers (regardless of own size) to associate FIGURE 2 | Following behavior of live guppies toward differently sized Robofish replicas. (A) Example track of a trial with Robofish. Fish were introduced into the start box (gray circle, lower left corner) and released into the tank after 1 min. Robofish then moved on a predefined zig-zag trajectory through the tank until it reached its start position. This movement was repeated a second time and a trial lasted about 60 s in total. (B) Inter-individual distance (Distance to Robofish) between live fish and differentially-sized Robofish replicas. Shown are medians and interquartile ranges (whiskers at 10 and 90% percentiles, with outliers as dots).
with those experienced, larger phenotypes. This concept is also assumed to be important in the context of social learning (see Laland, 2004). For example, younger thus smaller guppy females copy the mate choice of older thus larger ones more readily (Dugatkin and Godin, 1993;Amlacher and Dugatkin, 2005). (c) Larger individuals are more conspicuous than smaller ones (both visual and non-visual, see Faria et al., 2010). As staying together as a group is assumed beneficial (Krause and Ruxton, 2002) and environmental conditions such as water turbidity (Borner et al., 2015) often hamper prober sensing of conspecifics, following the larger, more conspicuous ones might be under positive selection.
Live fish's own body size did not affect following tendencies in our study, and although we found no evidence for size-assortative following, sorting by size is common in fish shoals in the wild . While fish might have a ubiquitous preference to follow larger individuals, assortative patterns might simply arise mechanistically because those front-occupying individuals may swim faster  leading to an assortment within the shoal . Also other factors may contribute to observed size assortative patterns, most likely mediated through competition among group members Croft et al., 2009). Body size in many species correlates with an individual's fighting ability and, as a result, its dominance status (Arnott and Elwood, 2009). In turn, smaller, inferior individuals despite their preference to associate with larger, dominant conspecifics, may suffer from increased aggressive behaviors, which, ultimately, can lead to exclusion from the group (Magurran and Seghers, 1991;Hoare et al., 2000). However, in Atlantic mollies (Poecilia mexicana), a species closely related to the guppy, size-assortative patterns in both sexes have been found to arise only gradually over time, e.g., a size-assortative pattern was only detectable after a few days of familiarization . Thus, despite a possible preference to associate and/or follower larger individuals, there are many other factors that can lead to observable size-assortative patterns in the wild Croft et al., 2009).
We found no evidence that follower's risk-taking behavior affected their tendencies to follow Robofish leaders of different size. This result is in contrast to studies in sticklebacks, where shyer individuals are better followers and are less likely to initiate leadership behavior themselves (Leblond and Reebs, 2006;Nakayama et al., 2012aNakayama et al., ,b, 2016. Besides possible species-specific differences, reinforcing feedbacks due to mutual influences among leaders and followers may have led to the observed personality-dependent following behavior in sticklebacks (Harcourt et al., 2009). Future comparative studies will thus help to pinpoint species-specific leadership behaviors as well as general patterns observable across taxa. However, this would need experimenters to take the same experimental approach while investigating a variety of different species -at least for small fish a platform like Robofish or similar tools that allow experimenters to adapt their systems quickly to different species (i.e., by exchanging the replicas) seem to be a promising tool for this.
A leader's behavior is clearly influencing its leadership success (Ioannou et al., 2015;Nakayama et al., 2016), but independent of size-specific behavioral differences, body size seems to be the sole determinant of success in leading others in a "bigger is better" way at least in the guppy. When followers have ubiquitous preferences for large-bodied leaders, theoretical and practical considerations of collective behavior will strongly benefit from accounting for these size-specific leadership patterns.
Biomimetic robots allow experimenters to gain control of the animal stimulus, which is not possible using live animals. This may have several potential benefits for the study of animal behavior: First, standardized test assays become available that minimize variation of provided stimuli which might come into play when using live animals as social stimulus. Second, decoupling of behavior and morphology becomes available as we demonstrated in the current study. Third, interactive robots can validate theoretically proposed mechanisms of social interactions among animals as their rules of interaction can be systematically manipulated. Forth, the use of artificial stimuli can reduce the number of individuals tested during experimentation (though same or better data quality is achieved through highly standardized test procedures) and thus help to promote the 3R principle in behavioral research. In fact, future developments may further increase the usage of biomimetic robots for example through use of soft robotic techniques (Li et al., 2017;Katzschmann et al., 2018;Gao et al., 2019) or a better integration of acoustic or olfactory cues which are important communication channels in many species, including fish (Ward et al., 2002;Bass and McKibben, 2003;Halloy et al., 2007). Biomimetic robots, together with video animations or virtual reality platforms, are thus able to strongly assist biologists answering a wide variety of question that could not be answered through the use of classical experimental assays.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material.

ETHICS STATEMENT
The experiments reported here comply with the current German laws approved by LaGeSo Berlin (Reg. 0117/16 to DB).

AUTHOR CONTRIBUTIONS
DB and JK designed the study. DB, HM, JK, and TL built the Robofish system. DB, HM, MH, and JL performed the experiments. DB and PR analyzed the data. DB wrote the first draft of the manuscript. All authors contributed to the final version of the manuscript and approved the submission.

ACKNOWLEDGMENTS
We would like to thank David Lewis for his valuable help in raising our test fish. Furthermore, we like to thank Angelika Szengel and Hai Nguyen for their help in developing the Robofish and Christopher Schutz and Apollonia Landes for their help during experimentation. This manuscript has been released as a pre-print at bioRxiv (Bierbach et al., 2018c).

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe. 2020.00441/full#supplementary-material FIGURE S1 | Photograph of differently-sized replicas. Left (20 mm SL) and middle (25 mm SL) replicas are unprocessed 3D printed blanks that were later on equipped with glass eyes and color-painted as shown for the 30 mm replica on the right. FIGURE S2 | Speed distributions during a trial (a) as well as speed profiles (b) of the three differently-sized replicas. In (a) fraction of speeds with a 1 cm/s bin size are shown for a trial run (60 s). There was no significant difference among replicas (Kolmogorov-Smirnoff test: P > 0.41). In (b) the speed profiles of Robofish with all three replicas are shown. There was no significant difference in accelerating/ decelerating among replicas (Kolmogorov-Smirnoff test: P > 0.99). Please note that this pattern of acceleration and deceleration is typically found between two turning points along the zig-zag movement. Maximum speed was set to 15 cm/s (dashed line).
VIDEO S1 | Example track of a live guppy following the Robofish replica.