A quarter-million-year-old polymorphism drives reproductive mode variation in the pea aphid

Although asexual linages evolved from sexual lineages in many different taxa, the genetics of sex loss remains poorly understood. We addressed this issue in the pea aphid Acyrthosiphon pisum, whose natural populations encompass lineages performing cyclical parthenogenesis (CP) and producing one sexual generation per year, as well as obligate parthenogenetic (OP) lineages that can no longer produce sexual females but can still produce males. A SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from six populations) showed that a single X-linked region controls the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2000 SNPs or short Indels showing a high degree of association with the phenotypic trait. Comparison of de novo genome assemblies built from long reads did not reveal large structural rearrangements between CP and OP lineages within the candidate region. This reproductive polymorphism still appears relatively ancient, as we estimated its age at ~0.25 million years from the divergence between cp and op alleles. The low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.4%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. This contagious asexuality, combined with environment-induced selection (each reproductive mode being favored under different climates) probably contributes to the long-term persistence of the cp and op alleles. ### Competing Interest Statement The authors have declared no competing interest.