Eusociality is a rare life history strategy that evolved repeatedly in Hymenoptera. At the population genetics level, inbreeding and low effective population size have been suggested as both evolutionary causes and consequences of social life. In this study, we tested these hypotheses by estimating the relative rate of non-synonymous substitution in 169 species to investigate the variation in natural selection efficiency and effective population size throughout the Hymenoptera tree of life. We show that relaxed selection is strongly associated with eusociality. This suggests that the division of reproductive labour decreases effective population size in ants, bees and wasps. Interestingly, we also report a striking and widespread relaxation of selection in both social and non social bees, which indicates that these keystone pollinator species generally feature low effective population sizes. We suggest that a particularly high inbreeding rate in bees might increase the benefits of kin selection, which would explain why most independent origins of eusociality in the tree of life occurred in this taxon. The particularly high load of deleterious mutations we report in the genome of these crucial ecosystem engineer species suggest new concerns about their ongoing population decline.
RESULTS
dN/dS distribution across the Hymenoptera phylogenyWe estimated dN/dS in 3241 gene alignments of 169 Hymenoptera species using the mapNH program (Romiguier et al. 2012, https://github.com/BioPP/testnh) from the testnh program suite (Dutheil and Boussau 2008;Guéguen and Duret 2018). Average genomic dN/dS values are displayed along the Hymenoptera tree on Figure 1. We used this tree obtained by Peters et al. (2017) and its topology thorough all analyses. The largest and smallest mean ratios were inferred for Eucera nigrescens (0.0745) and Macrocentrus marginator (0.0318). As expected for conserved coding regions, the distribution of genomic dN/dS ratios is close to 0 (overall average of 0.0424), indicative of the large prevalence of purifying selection. We observed above average dN/dS ratios in all 5 available eusocial clades: Formicidae (0.0463 ± 0.0019, 3 species), Polistinae/Vespinae wasps (0.0492 ± 0.0023, 3 species), Stenogastrinae (0.0447, 1 species only), the Apis/Bombus/Tetragonula clade (0.0511 ± 0.0057), and Halictidae (0.0613 ± 0.0039; 2 species). These last two clades of bees do not clearly stand out however, as most bees in the dataset (Antophila, species characterized by pollen feeding of larvae: Apidae, Megachilidae, Halictidae, Colettidae, Andrenidae, and Melittidae) show high dN/dS ratios (0.0522 ± 0.0073, 41 species) with no clear dependence on their social organization. Finally, only two purely solitary taxa displayed comparable dN/dS ratios: Siricoidea (0.0466 ± 0.0027, 3 species) and Cynipoidea (0.0447 ± 0.0059, 5 species). We further used simple linear modelling to try and relate variation in dN/dS ratios to life history traits. When needed, phylogenetic independent contrasts where used to account for phylogene...