Evolutionary theory predicts that sexually antagonistic mutations accumulate differentially on the X chromosome and autosomes in species with an XY sex-determination system, with effects (masculinization or feminization of the X) depending on the dominance of mutations. Organisms with alternative modes of inheritance of sex chromosomes offer interesting opportunities for studying sexual conflicts and their resolution, because expectations for the preferred genomic location of sexually antagonistic alleles may differ from standard systems. Aphids display an XX/X0 system and combine an unusual inheritance of the X chromosome with the alternation of sexual and asexual reproduction. In this study, we first investigated theoretically the accumulation of sexually antagonistic mutations on the aphid X chromosome. Our results show that i) the X is always more favourable to the spread of male-beneficial alleles than autosomes, and should thus be enriched in sexually antagonistic alleles beneficial for males, ii) sexually antagonistic mutations beneficial for asexual females accumulate preferentially on autosomes, iii) in contrast to predictions for standard systems, these qualitative results are not affected by the dominance of mutations. Under the assumption that sex-biased gene expression evolves to solve conflicts raised by the spread of sexually antagonistic alleles, one expects that male-biased genes should be enriched on the X while asexual female-biased genes should be enriched on autosomes. Using gene expression data (RNA-Seq) in males, sexual females and asexual females of the pea aphid, we confirm these theoretical predictions. Although other mechanisms than the resolution of sexual antagonism may lead to sex-biased gene expression, we argue that they could hardly explain the observed difference between X and autosomes. On top of reporting a strong masculinization of the aphid X chromosome, our study highlights the relevance of organisms displaying an alternative mode of sex chromosome inheritance to understanding the forces shaping chromosome evolution.

The peach-potato aphid, Myzus persicae (Sulzer), has a complex reproductive mode in which some lineages reproduce by continuous parthenogenesis, whereas others reproduce sexually once a year. The climate is thought to act directly on the reproductive mode, because sexual eggs are the only form that can resist frost in cold regions. Sexual reproduction necessitates an obligatory host alternation that may result in long-distance dispersal. Here, we examined the genetic variability at seven microsatellite loci of populations of M. persicae in France, where both reproductive modes occur. We provide clear genetic evidence that the breeding system affects genotypic variability, as cyclically parthenogenetic aphids are far more variable than their obligately parthenogenetic counterparts. A temporal decrease in genetic variability and a temporal genetic differentiation effect suggest the existence of selective factors that play an important role in shaping the genetic structure of M. persicae populations. Lastly, differences in the population structure between reproductive modes suggest that the migration associated with the change of host during sexual reproduction lowers the level of population differentiation.

Cyclically parthenogenetic organisms may have facultative asexual counterparts. Such organisms, including aphids, are therefore interesting models for the study of ecological and genetic interactions between lineages differing in reproductive mode. Earlier studies on aphids have revealed major differences in the genetic outcomes of populations that are possibly resulting mostly either from sexual or from asexual reproduction. Besides, notable gene flow between sexual and asexual derivatives has been suspected, which could lead to the emergence of new asexual lineages. The present study examines the interplay between these lineages and is based on analyses of population structure of individuals that may contribute to the pool of sexual reproductive forms in the host alternating aphid Rhopalosiphum padi. Using a Bayesian assignment method, we first show that the sexual forms of R. padi on mating sites encompass two genetically distinct clusters of individuals in the western part of France. The first cluster included unique genotypes of sexual lineages, while the second cluster included facultatively asexual lineages in numerous copies, the reproductive mode of the two clusters being confirmed by reference clones. Sexual reproductive forms produced by sexual and facultatively asexual lineages are thus admixed at mating sites which gives a large opportunity for the two clusters to mate with each other. Nevertheless, this study also highlights, as previously demonstrated, that the two clusters retained high genetic differentiation. Possible explanations for the inferred limited genetic exchanges are advanced in the discussion, but further dedicated investigations are required to solve this paradox.

Associations between microbes and animals are ubiquitous and hosts may benefit from harbouring microbial communities through improved resource exploitation or resistance to environmental stress. The pea aphid, Acyrthosiphon pisum, is the host of heritable bacterial symbionts, including the obligate endosymbiont Buchnera aphidicola and several facultative symbionts. While obligate symbionts supply aphids with key nutrients, facultative symbionts influence their hosts in many ways such as protection against natural enemies, heat tolerance, color change and reproduction alteration. The pea aphid also encompasses multiple plant-specialized biotypes, each adapted to one or a few legume species. Facultative symbiont communities differ strongly between biotypes, although bacterial involvement in plant specialization is uncertain. Here, we analyse the diversity of bacterial communities associated with nine biotypes of the pea aphid complex using amplicon pyrosequencing of 16S rRNA genes. Combined clustering and phylogenetic analyses of 16S sequences allowed identifying 21 bacterial OTUs (Operational Taxonomic Unit). More than 98% of the sequencing reads were assigned to known pea aphid symbionts. The presence of Wolbachia was confirmed in A. pisum while Erwinia and Pantoea, two gut associates, were detected in multiple samples. The diversity of bacterial communities harboured by pea aphid biotypes was very low, ranging from 3 to 11 OTUs across samples. Bacterial communities differed more between than within biotypes but this difference did not correlate with the genetic divergence between biotypes. Altogether, these results confirm that the aphid microbiota is dominated by a few heritable symbionts and that plant specialization is an important structuring factor of bacterial communities associated with the pea aphid complex. However, since we examined the microbiota of aphid samples kept a few generations in controlled conditions, it may be that bacterial diversity was underestimated due to the possible loss of environmental or transient taxa.

Asexuality confers demographic advantages to invasive taxa, but generally limits adaptive potential for colonizing of new habitats. Therefore, pre-existing adaptations and habitat tolerance are essential in the success of asexual invaders. We investigated these key factors of invasiveness by assessing reproductive modes and host-plant adaptations in the pea aphid, Acyrthosiphon pisum, a pest recently introduced into Chile. The pea aphid encompasses lineages differing in their reproductive mode, ranging from obligatory cyclical parthenogenesis to fully asexual reproduction. This species also shows variation in host use, with distinct biotypes specialized on different species of legumes as well as more polyphagous populations. In central Chile, microsatellite genotyping of pea aphids sampled on five crops and wild legumes revealed three main clonal genotypes, which showed striking associations with particular host plants rather than sampling locations. Phenotypic analyses confirmed their strong host specialization and demonstrated parthenogenesis as their sole reproductive mode. The genetic relatedness of these clonal genotypes with corresponding host-specialized populations from the Old World indicated that each clone descended from a particular Eurasian biotype, which involved at least three successful introduction events followed by spread on different crops. This study illustrates that multiple introductions of highly specialized clones, rather than local evolution in resource use and/or selection of generalist genotypes, can explain the demographic success of a strictly asexual invader.

The pea aphid, Acyrthosiphon pisum, encompasses distinct host races specialized on various Fabaceae species, but the extent of genetic divergence associated with ecological specialization varies greatly depending on plant and geographic origins of aphid populations. Here, we studied the genetic structure of French sympatric pea aphid populations collected on perennial (pea and faba bean) and annual (alfalfa and red clover) hosts using 14 microsatellite loci. Classical and Bayesian population genetics analyses consistently identified genetic clusters mostly related to plant origin: the pea/faba bean cluster was highly divergent from the red clover and the alfalfa ones, indicating they represent different stages along the continuum of genetic differentiation. Some genotypes were assigned to a cluster differing from the one expected from their plant origin while others exhibited intermediate genetic characteristics. These results suggest incomplete barriers to gene flow. However, this limited gene flow seems insufficient to prevent ecological specialization and genetic differentiation in sympatry.

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome-wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host-adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.

Variation in traits affecting preference for, and performance on, new habitats is a key factor in the initiation of ecological specialisation and adaptive speciation. However, habitat and resource use also involves other traits whose influence on ecological and genetic divergence remains poorly understood. In the present study, we investigated the extent of variation of life-history traits among sympatric populations of the pea aphid Acyrthosiphon pisum, which shows several host races that are specialised on various plants of the family Fabaceae plants and is an established model for ecological speciation. First, we assessed the community structure of microbial partners within host populations of the pea aphid. The effect of these microbes on host fitness is uncertain, although there is growing evidence that they may modulate various important adaptive traits of their host such as plant utilisation and resistance against natural enemies. Second, we performed a multivariate analysis on several ecologically relevant features of host populations recorded in the present and previous studies (including microbial composition, colour morph, reproductive mode, and male dispersal phenotype), enabling the identification of correlations between phenotypic traits. We discuss the ecological significance of these associations of traits in relation to the habitat characteristics of pea aphid populations, and their consequences for the evolution of ecological specialisation and sympatric speciatio