Dual pathways in social evolution: Population genetic structure of group-living and solitary species of kleptoparasitic spiders (Argyrodinae: Theridiidae)

Kleptoparasites of social spider colonies do not track hosts' subdivided population structure

Organisms with lower dispersal abilities tend to have more genetically dissimilar populations. The same is true for parasites, whose transmission frequency may depend on the population structure of the host. This should be especially true when hosts and parasites face similar barriers to dispersal. Here, we considered the similarities between host and parasite population structure in a social spider system. In this system, host colonies are typified by rapid growth via internal recruitment followed by budding or fission events when colonies grow too large, with each colony representing a distinct population. Host colonies provide habitat for kleptoparasitic spiders, which steal prey from and may also feed directly on host individuals. We asked whether kleptoparasites exhibit a similar degree of population subdivision as their host. Under the free-mixing hypothesis (i.e., horizontal transmission), kleptoparasites would have well-mixed populations across broader regions than a single host nest, whereas host populations would be strongly genetically structured. Under the host-tracking hypothesis (i.e., vertical transmission), kleptoparasites would have a population structure that parallels that of the host. We conducted a genotype-by-sequencing study to assess the population structure of both hosts and kleptoparasites within three nearby regions in eastern Ecuador. We found strong signatures of population differentiation and bottlenecks in the host species, which is congruent with past studies. However, we found that kleptoparasite populations were well mixed across host nests, with no evidence of recent bottlenecks. These results support our free-mixing hypothesis, suggesting that kleptoparasites follow patterns of horizontal transmission in this social spider system.

Foraging Payoffs Change With Group Size in Kin and Non-kin Groups of an Argyrodinae Kleptoparasitic Spider, Argyrodes miniaceus

Evolutionary transitions from solitary to group-living are ubiquitous in animal systems. While the fitness consequences of group size changes are often investigated, the long-standing debate on whether kinship is a prerequisite of sociality is still ongoing. In the current study, we used kleptoparasitic spiders Argyrodes miniaceus (subfamily Argyrodinae, Theridiidae) as a model system to assess the role of group size on the foraging payoffs of kin and non-kin groups. We set up laboratory-manipulated kin and non-kin foraging groups and used feeding occurrence and duration as proxies for foraging benefits and feeding latency and the number of host attacks as estimates of foraging costs. Compared to solitary individuals, feeding durations of successfully fed individuals in groups was not significantly different from that of solitary foragers in both kin and non-kin groups. The occurrences of feeding decreased significantly in group sizes two and above, in non-kin groups, and in group sizes three and above, in kin groups. In kin groups, groups size two had significantly shorter feeding latencies compared to other group sizes, even though feeding duration did not change systematically with group size. Similarly, the number of attacks from the hosts were highest in non-kin groups with more than two individuals and in kin groups with more than three individuals. The juxtaposition of kin and non-kin group showed that A. miniaceus enjoyed the highest foraging payoffs when being solitary or in small groups (group size two). However, host attacks appeared to hamper feeding occurrences in kin groups, which was not observed in non-kin groups. Our results contrast sharply with the feeding benefits of kinship recorded in kin-based groups of sub-social species present in related subfamilies in the Theridiidae.

Northward geographic diversification of a kleptoparasitic spider Argyrodes lanyuensis (Araneae, Theridiidae) from the Philippine Archipelago to Orchid Island

Oceanic islands are unique geographic systems that promote local adaptations and allopatric speciation in many of their highly endemic taxa. This is a common case in the Philippine Archipelago, where numerous unrelated taxa on islands have been inferred to have diversified in isolation. However, few cases have been reported in invertebrates especially among parasitic organisms. Here, we tested for biogeographical structure in novel populations of the "generalist" kleptoparasitic spider, Argyrodes lanyuensis Yoshida, Tso & Severinghaus, 1998 in the Philippines. Results showed that, in addition to Orchid/Lanyu Island, this species has a wide geographic distribution in the Philippine Archipelago. The estimated divergence time of this lineage using the mitochondrial cytochrome oxidase 1 (mt-CO1) suggests that this species diverged ca 3.12 MYA, during the Pliocene. Two reciprocal monophyletic clades were elucidated in A. lanyuensis, but with limited differentiation across Pleistocene Aggregate Island Complex (PAIC) boundaries and modern-day islands. However, in our analyses of morphological variation, we identified two phenotypically differentiated units in males (Orchid Island, Taiwan + Luzon, Philippine PAIC populations vs. Palawan + West Visayan + Mindanao PAIC populations). We infer that this species diverged in the southern portion of the Philippine Archipelago and only recently colonized Orchid Island. Our study provides new information on the extensive distribution of A. lanyuensis outside Orchid Island, Taiwan, but we documented a very limited geographically associated genetic variation. Our study points to behavioral phenomena such as foraging behavior as essential contributor to the evolutionary process of species diversification, in contrast to the traditionally invoked geographic drivers of divergence.

Machine Learning Data Imputation and Prediction of Foraging Group Size in a Kleptoparasitic Spider

Cost–benefit analysis is widely used to elucidate the association between foraging group size and resource size. Despite advances in the development of theoretical frameworks, however, the empirical systems used for testing are hindered by the vagaries of field surveys and incomplete data. This study developed the three approaches to data imputation based on machine learning (ML) algorithms with the aim of rescuing valuable field data. Using 163 host spider webs (132 complete data and 31 incomplete data), our results indicated that the data imputation based on random forest algorithm outperformed classification and regression trees, the k-nearest neighbor, and other conventional approaches (Wilcoxon signed-rank test and correlation difference have p-value from < 0.001–0.030). We then used rescued data based on a natural system involving kleptoparasitic spiders from Taiwan and Vietnam (Argyrodes miniaceus, Theridiidae) to test the occurrence and group size of kleptoparasites in natural populations. Our partial least-squares path modelling (PLS-PM) results demonstrated that the size of the host web (T = 6.890, p = 0.000) is a significant feature affecting group size. The resource size (T = 2.590, p = 0.010) and the microclimate (T = 3.230, p = 0.001) are significant features affecting the presence of kleptoparasites. The test of conformation of group size distribution to the ideal free distribution (IFD) model revealed that predictions pertaining to per-capita resource size were underestimated (bootstrap resampling mean slopes <IFD predicted slopes, p < 0.001). These findings highlight the importance of applying appropriate ML methods to the handling of missing field data.