Our study system (3,347 nests across 22 ha) is the largest known European polydomous system of F. exsecta Nylander, 1846, located in a semi-wet meadow in Central Romania in the Eastern Carpathians [54]. This system also hosts the myrmecopathogenic fungus P. formicae (Entomophthoromycota, Fig. 1a and Fig. 1b) formerly reported by us as P. myrmecophaga [see 26], but renamed following the description of Małagočka et al. [27].
In a first experiment, we randomly selected 80 F. exsecta nests from the central part of the supercolony, maintaining a minimum distance of 3 m between neighboring nests. All nests were checked for infected ants before and during the experiments by carefully inspecting vegetation within 0.5 m of the nest margin. As nest mound size is a fair indicator of the number of F. exsecta ants residing within [55] we calculated the above-ground volume of each experimental nest mound [56]: V = 1/2 × π × r1 × r2 × h, where r1 is the largest radius at the bottom, r2 is perpendicular to r1, and h is the height aboveground.
Ant corpses were obtained by collecting live, uninjured individuals from the nest surface two days before the experiments. Ants were placed in small vials, where they died within hours. To ensure that no individuals were initially infected by the fungus, fresh corpses were placed on moist cotton in vials and kept in a cool, dark room for two days until the experiment. This method promotes fungal growth [25], allowing us to exclude ants that were already infected. None of the corpses showed infection. To avoid potential nest-specific effects of chemical cues, experimental corpses were returned to their nest of origin. The experiments were conducted on the 14th and 17th of August 2012, with 40 nests tested per day. Additionally, 432 nests were checked for fungal prevalence during the study period. All dead ants attached to grass blades were checked for fungal infection as described above.
To mimic the appearance of fungus-killed ants, we fixed a single experimental corpse (hereafter ‘nest corpse’) to a Festuca pratensis grass blade using a minutia pin through the thorax, at ~ 8 cm height and ~ 1 cm from the nest mound edge—consistent with the typical position of infected individuals [25–26]. The exact location of the grass blade was chosen randomly. Formica exsecta ants are active across the nest surface and nearby vegetation. To test the effect of distance, another corpse (‘distant corpse’) was pinned in the same way 0.5 m away from the nest mound along the same axis. In addition, to test the ants’ reaction to a non-biological object, a white ant-sized dummy (constructed from three small polystyrene balls) was pinned to a grass blade ~ 1 cm from the nest mound edge, on the opposite side of the nest. We used only a ‘nest dummy’ (and not a ‘distant dummy’), as our aim was to assess whether ants would treat a non-corpse object as foreign material to be removed from the immediate nest surroundings, rather than to compare it with corpses across distances. Corpses and dummies were placed out 10 min before the first observation in the afternoon, when F. exsecta activity is higher (between 4 and 6 PM). At each nest, we performed eight 1-min scan observations, separated by 15-min intervals. Four observers simultaneously monitored 10 nests per day. We recorded the number of ants active on the vegetation surrounding each nest. Discovery was defined as the first physical contact recorded during these scans. Because observations were conducted every 15 minutes, some first contacts may have occurred between scans. However, this approach provides a conservative and standardized measure of detection probability across treatments, ensuring comparability between nest-corpses, distant corpses, and dummies. Disposal was defined as the complete removal of the corpse from the pin, and the time to disposal was recorded. After the observation session, any remaining corpse or dummy was left in place. Their status was checked the following day at 4 PM, and again on day 3 (for the first group of nests) and day 4 (for the second group).
In a follow-up experiment conducted in August 2025, we compared discovery and removal rates between corpses/dummies placed on grass blades versus directly on the ground to test whether corpse removal from elevated positions (grass blades) reflects a specific adaptive response to summit disease rather than general necrophoric behavior. We randomly selected 20 F. exsecta nests from the same population as in the 1st experiment and assigned them to four treatments (n = 5 nests per treatment): (1) corpse on a grass blade ~ 1 cm from the nest edge, (2) corpse on the ground ~ 1 cm from the nest edge, (3) dummy on a grass blade ~ 1 cm from the nest edge, and (4) dummy on the ground ~ 1 cm from the nest edge. All corpses were prepared as described above. Observation protocols followed those of the first experiment, except that in this case, we performed six 1-minute scans over 90 minutes for each nest. Two observers monitored 10 nests simultaneously.
Statistical analyses
The relationship between worker activity on grass blades and nest size was assessed using a Linear Mixed Model (LMM, N = 640 obs.) approach. Nest size was introduced as an independent variable, while observation series and nest code were applied as random factors. Differences in corpse/dummy discovery rates were analyzed with the help of a Cox regression model (proportional hazard approach, N = 240 corpses). The type of corpse, the mound size and the average number of ants active on grass blades were applied as independent variables. Nest code was included as a random factor in order to handle dependencies. A similar approach was applied for the analysis of corpse removal rates. In the follow-up experiment, differences in the probability of corpse discovery and removal among treatments (corpse on the ground, corpse on the grass, dummy on the ground and dunny on the grass, N = 80 corpses/dummies) were assessed using Fisher’s exact tests due to small expected counts in some categories. Then, pairwise comparisons of proportions were performed using pairwise.prop.test with Benjamini–Hochberg correction to account for multiple testing [57]. Proportions of discovery and removal were calculated for each treatment and are reported as percentages. All statistical analyses were carried out using R 4.1.0 (R Development Core Team 2021 [58]). LMM was performed using the lmer function in the lme4 package [59]. Cox regression analysis was performed with the coxme package [60]. Tukey’s HSD test was used to calculate post-hoc comparisons for each factor using the glht function in the multcomp package [61]. ggplot2 package [62] was used for graphs.