Temporal stability of sex ratio distorter prevalence in natural populations of the isopod Armadillidium vulgare

In the terrestrial isopod Armadillidium vulgare, many females produce progenies with female-biased sex ratios due to two feminizing sex ratio distorters (SRD): Wolbachia endosymbionts and a nuclear non-mendelian locus called the f element. To investigate the potential impact of these SRD on the evolution of host sex determination, we analyzed their temporal distribution in six A. vulgare populations sampled between 2003 and 2017, for a total of 29 time points. SRD distribution was heterogeneous among populations despite their close geographic locations, so that when one SRD was frequent in a population, the other SRD was rare. In contrast with spatial heterogeneity, our results overall did not reveal substantial temporal variability in SRD prevalence within populations, suggesting equilibria in SRD evolutionary dynamics may have been reached or nearly so. Temporal stability was also generally reflected in mitochondrial and nuclear variation. Nevertheless, in a population, a Wolbachia strain replacement coincided with changes in mitochondrial composition but no change in nuclear composition, thus constituting a typical example of mitochondrial sweep caused by endosymbiont rise in frequency. Rare incongruence between Wolbachia strains and mitochondrial haplotypes suggested the occurrence of intraspecific horizontal transmission, making it a biologically relevant parameter for Wolbachia evolutionary dynamics in A. vulgare. Overall, our results provide an empirical basis for future studies on SRD evolutionary dynamics in the context of multiple sex determination factors co-existing within a single species, to ultimately evaluate the impact of SRD on the evolution of host sex determination mechanisms and sex chromosomes.