Predominance of sexual reproduction in Romanian populations of the aphid Sitobion avenae inferred from phenotypic and genetic structure

Models of coexistence of sexual and asexual lineages in aphids assume that obligate parthenogenetic lineages predominate in areas with mild winter climate because of their high reproductive output, while sexual lineages predominate in areas with severe winter because they produce eggs resistant to frost. To validate this hypothesis in natural conditions, the reproductive mode of populations of the aphid Sitobion avenae was assessed in two very contrasting climatic situations, Romania (severe winter) and Western France (mild winter). To achieve this, reproductive modes were inferred from both (1) the population composition in sexual and asexual forms in autumn, and (2) the genetic structure of Romanian and French populations of S. avenae using microsatellite markers. Romanian populations encompassed a high proportion of sexual forms and were characterised by a very high genotypic diversity and low linkage disequilibrium. In constrast, the French population showed frequent linkage disequilibria, low genetic diversity, and high level of clonal amplification with two asexual genotypes representing over 60% of the sample. In agreement with the model's predictions, these results clearly indicate that sexual reproduction in S. avenae is predominant under the continental climate of Romania, while asexual lineages prevail under the oceanic climate of Western France.     

Seasonal and annual genotypic variation and the effect of climate on population genetic structure of the cereal aphid Sitobion avenae in northern France

Changes in the genetic structure and genotypic variation of the aphid Sitobion avenae collected from cereal crops in northern France were examined by analysing variation at five microsatellite loci across several years and seasons. Little regional and temporal differentiation was detected, as shown by very low FST among populations. Repeated genotypes, significant heterozygote deficits, positive FIS values and frequent linkage disequilibria were found in nearly all samples, suggesting an overall pattern of reproductive mode variation in S. avenae populations. In addition, samples from Brittany (Bretagne) showed greater signs of asexual reproduction than those from the north of France, indicating a trend toward increasing sexuality northward. These patterns of reproductive variation in S. avenae are consistent with theoretical models of selection of aphid reproductive modes by climate. Contrasting with little changes in allelic frequencies, genotypic composition varied substantially in time and, to a lesser extent, in space. An important part of changes in genotypic arrays was due to the variation in frequency distribution of common genotypes, i.e. those that were found at several instances in the samples. Genotypic composition was also shown to vary according to climate, as genotypic diversity in spring was significantly correlated with the severity of the previous winter and autumn. We propose that the genetic homogeneity among S. avenae populations shown here across large temporal and spatial scales is the result of two forces: (i) migration conferred by high dispersal capabilities, and (ii) selection over millions of hectares of cereals (mostly wheat) bred from a narrow genetic base.     

Genetic structure of an aphid studied using microsatellites: cyclic parthenogenesis, differentiated lineages and host specialization

In a previous study, samples of the grain aphid Sitobion avenae (F.) were collected from wheat and adjacent cocksfoot hosts in a population thought to be primarily parthenogenetic, and DNA from individual aphids was analysed with a multilocus technique. Here we have applied single-locus microsatellites and a mitochondrial DNA marker to a subset of the same DNA extracts, and have made several additional inferences about important genetic and population processes in S. avenae . Microsatellite analysis indicated very high levels of genic and genotypic variation. S. avenae fell into three genotypic groups inferred to be almost noninterbreeding, while analysis of linkage and Hardy-Weinberg equilibria suggested high levels of sexual recombination within each genotypic group. Host specialization was evident: one lineage was found only on wheat, and one (bearing many alleles inferred to be introgressed from the blackberry-grass aphid S. fragariae (Walker)) was found only on cocksfoot. The third group of interrelated genotypes was found commonly on both hosts. Although most genotypes were found only once, some were much more numerous in the sample than expected from the frequency of the alleles they contained. This, and rapid temporal changes in genotypic composition of samples, indicates strong selective differences between genotypes and lineages. In the major genotypic group, the commonest genotypes were significantly more homozygous than were rare ones: thus these data may help to explain the frequent observation of homozygous excess in aphid allozymes. The genotype group showing S. avenae -like as well as S. fragariae -like alleles also carried S. fragariae -like mitochondrial DNA in at least 25/31 cases, indicating gender-asymmetrical hybridization.     

Spatial and temporal genetic variability in French populations of the peach-potato aphid,Myzus persicae

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.     

Climate and agricultural context shape reproductive mode variation in an aphid crop pest

In aphids, reproductive mode is generally assumed to be selected for by winter climate. Sexual lineages produce frost-resistant eggs, conferring an advantage in regions with cold winters, while asexual lineages predominate in regions with mild winters. However, habitat and resource heterogeneities are known to exert a strong influence on sex maintenance and might modulate the effect of climate on aphid reproductive strategies. We carried out a hierarchical sampling in northern France to investigate whether reproductive mode variation of the aphid Rhopalosiphum padi is driven by winter climate conditions, by habitat and resource heterogeneities represented by a range of host plants or by both factors. We confirmed the coexistence in R. padi populations of two genetic clusters associated with distinct reproductive strategies. Asexual lineages predominated, whatever the surveyed year and location. However, we detected a between-year variation in the local contribution of both clusters, presumably associated with preceding winter severity. No evidence for host-driven niche differentiation was found in the field on six Poaceae among sexual and asexual lineages. Two dominant multilocus genotypes (∼70% of the sample), having persisted over a 10-year period, were equally abundant on different plant species and locations, indicating their large ecological tolerance. Our results fit theoretical predictions of the influence of winter climate on the balance between sexual and asexual lineages. They also highlight the importance of current agricultural practices which seem to favour a small number of asexual generalist genotypes and their migration across large areas of monotonous environments.     

Grain aphid population structure: no effect of fungal infections in a 2-year field study in Denmark

1 Sitobion avenae (F.) is a serious pest in Danish cereal crops. To understand the population genetic structure, aphids were sampled in seven different winter wheat (Triticum sativum Lamarck) fields throughout Denmark. The aphids were genotyped with seven microsatellite markers. In total, 2075 aphids were collected and 1203 of these were genotyped. 2 The Danish S. avenae populations displayed very high genotypic diversity, high percentages of unique genotypes and low linkage disequilibria; this is likely to be a result of genetic recombination encompassed by their holocyclic lifestyle. The populations showed very limited differentiation and no sign of isolation by distance. Almost all the genetic variation was ascribed within the populations rather than between populations, probably due to a high migration rate at approximate 10% per generation. 3 Seasonal changes in clonal diversity and distribution of asexual summer generations of S. avenae within the infestation period in a single winter wheat field were followed over two consecutive years by weekly sampling from 60 plots each of 20 × 20 m. Clonal diversity was high in all samples with no dominant clonal lineages and no significant difference in the genotypic diversity between weeks or between years. However, a temporal genetic differentiation effect, throughout the infestation, suggests that selective factors or high temporal migration play an important role in shaping the genetic structure S. avenae. 4 Analyses of fungal infected and uninfected aphids were performed to test whether some clonal linage were more often infected by fungi from the Entomophthorales under field conditions. In total, 54 progeny from aphids with Entomophthorales were genotyped and compared with 422 uninfected aphid genotypes. The Entomophthorales‐infected aphid genotypes did not cluster out together, suggesting that these fungal pathogens did not affect the population differentiation or clonal distribution of S. avenae in a Danish agroecosystem. 5 Our findings indicate that S. avenae populations can be controlled using conservation biological control     

Genetic interactions influence host preference and performance in a plant-insect system

Phytophagous insects generally feed on a restricted range of host plants, using a number of different sensory and behavioural mechanisms to locate and recognize their host plants. Phloem-feeding aphids have been shown to exhibit genetic variation for host preference of different plant species and genetic variation within a plant species can also have an effect on aphid preference and acceptance. It is known that genotypic interactions between barley genotypes and Sitobion avenae aphid genotypes influence aphid fitness, but it is unknown if these different aphid genotypes exhibit active host choice (preference) for the different barley genotypes. Active host choice by aphid genotypes for particular plant genotypes would lead to assortative association (non-random association) between the different aphid and plant genotypes. The performance of each aphid genotype on the plant genotypes also has the ability to enhance these interactions, especially if the aphid genotypes choose the plant genotype that also infers the greatest fitness. In this study, we demonstrate that different aphid genotypes exhibit differential preference and performance for different barley genotypes. Three out of four aphid genotypes exhibited preference for (or against) particular barley genotypes that were not concordant with differences in their reproductive rate on the specific barley genotype. This suggests active host choice of aphids is the primary mechanism for the observed pattern of non-random associations between aphid and barley genotypes. In a community context, such genetic associations between the aphids and barley can lead to population-level changes within the aphid species. These interactions may also have evolutionary effects on the surrounding interacting community, especially in ecosystems of limited species and genetic diversity.     

Admixed sexual and facultatively asexual aphid lineages at mating sites

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.     

Evolutionary and functional insights into reproductive strategies of aphids

Aphids are among the few organisms capable of reproducing either sexually or asexually. This plasticity in reproductive mode is viewed as an adaptive response to cope with seasonal changes. Clonal reproduction occurs during the growing season allowing rapid population increase, while sexual reproduction occurs during late summer and leads to frost-resistant eggs that can survive winter conditions. This shift between these two extreme reproductive modes is achieved by using the same genotype, i.e. within the same genetic clone, and is triggered by photoperiodic changes perceived by the aphid brain or visual system. Advances have been made recently to depict genetic programs that relate to the regulation of reproductive modes in aphids. These studies have benefited from the rapid development of genomic and post-genomic resources obtained through the International Aphid Genomics Consortium. Here, we underline the importance of several candidate genes in the switch from clonal to sexual reproduction in aphids and whose roles await full validation. Besides reproductive mode variation expressed at the genotypic level, aphid species also frequently encompass lineages which have lost the sexual phase and hence the alternating clonal and sexual reproductive phases of the life cycle. This coexistence of sex and asexual reproduction within the same species raises questions on its evolutionary and ecological significance. We summarize the knowledge accumulated to date on the maintenance of sex as well as on the origin and evolution of asexuality in aphids. By combining functional genomics, genetic and ecological approaches on reproductive plasticity and polymorphism, we hope to obtain an integrative view of the evolutionary forces shaping aphid reproductive strategies, from gene to population and species levels.     

Evidence for predominant clones in a cyclically parthenogenetic organism provided by combined demographic and genetic analyses

Aphids are particularly interesting models in the study of genetic and demographic components of plant adaptation because of their breeding system which combines parthenogenesis and sexual reproduction (i.e. cyclical parthenogenesis), and the frequent emergence of host-adapted races reported in this group. In this paper, patterns of host adaptation were assessed on local populations of the aphid Sitobion avenae by following their demographic and genetic structure in a maize field for two consecutive years. The existence of putative generalist (polyphagous) or specialized (host-adapted) genotypes was also investigated by comparing the genotypic distribution of this aphid on maize and other cultivated host plants, using five microsatellite loci. Although population dynamics revealed strong variation in aphid abundance during the colonization period on maize, two genotypes identified at seven additional microsatellite loci were predominant and exhibited stable frequencies over cropping season and between years. Based on present and earlier studies, these two prevalent genotypes were shown to survive on different host plants other than maize, to colonize large geographical zones and to persist parthenogenetically for several years. All these data strongly suggest that these two genotypes are asexual generalist clones that could have been favoured by agricultural practices encountered in western Europe. Besides these two clones, a continual replacement of rare genotypes was observed on maize in both years. Hypotheses involving selection via aphid-plant interactions and natural enemies were proposed for explaining the disappearance of these genotypes on maize.     

Genetic architecture of sexual and asexual populations of the aphid Rhopalosiphum padi based on allozyme and microsatellite markers

Cyclical parthenogens, including aphids, are attractive models for comparing the genetic outcomes of sexual and asexual reproduction, which determine their respective evolutionary advantages. In this study, we examined how reproductive mode shapes genetic structure of sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) populations of the aphid Rhopalosiphum padi by comparing microsatellite and allozyme data sets. Allozymes showed little polymorphism, confirming earlier studies with these markers. In contrast, microsatellite loci were highly polymorphic and showed patterns very discordant from allozyme loci. In particular, microsatellites revealed strong heterozygote excess in asexual populations, whereas allozymes showed heterozygote deficits. Various hypotheses are explored that could account for the conflicting results of these two types of genetic markers. A strong differentiation between reproductive modes was found with both types of markers. Microsatellites indicated that sexual populations have high allelic polymorphism and heterozygote deficits (possibly because of population subdivision, inbreeding or selection). Little geographical differentiation was found among sexual populations confirming the large dispersal ability of this aphid. In contrast, asexual populations showed less allelic polymorphism but high heterozygosity at most loci. Two alternative hypotheses are proposed to explain this heterozygosity excess: allele sequence divergence during long-term asexuality or hybrid origin of asexual lineages. Clonal diversity of asexual lineages of R. padi was substantial suggesting that they could have frozen genetic diversity from the pool of sexual lineages. Several widespread asexual genotypes were found to persist through time, as already seen in other aphid species, a feature seemingly consistent with the general-purpose genotype hypothesis.     

Genetic structure and clonal diversity of an introduced pest in Chile, the cereal aphid Sitobion avenae

In Chile, the aphid Sitobion avenae is of recent introduction, lives on cultivated and wild Poaceae, and is thought to reproduce by permanent parthenogenesis. In order to study the genetic variability and population structure of this species, five microsatellite loci were typed from individual aphids collected from different cultivated and wild host plants, from different geographical zones, and years. Chilean populations showed a high degree of heterozygosity and a low genetic variability across regions and years, with four predominant genotypes representing nearly 90% of the sample. This pattern of low clonal diversity and high heterozygosity was interpreted as the result of recent founder events from a few asexually reproducing genotypes. Most geographical and temporal variation observed in the genetic composition resulted from fluctuations of a few predominant clones. In addition, comparisons of the genotypes found in Chile with those described in earlier surveys of S. avenae populations in Western Europe led us to identify 'superclones' with large geographical distribution and high ecological success, and to make a preliminary exploration of the putative origin(s) of S. avenae individuals introduced to Chile.     

Gene flow between sexual and facultatively asexual lineages of an aphid species and the maintenance of reproductive mode variation

Many organisms considered as strictly clonal may in fact experience some rare events of sexual reproduction with their sexual relatives. However, the rate of sexual–asexual gene flow has rarely been assessed mainly because its evaluation is difficult to achieve in the field. In the cyclically parthenogenetic aphid Rhopalosiphum padi , two main sets of lineages, differing in their investment in sexual reproduction and in their genetic attributes, co-exist even at a very fine scale: the 'sexual' lineages which have a full commitment to the sexual reproduction, and the 'facultatively asexual' lineages, which allocate investment in the sexual and parthenogenetic reproduction. This system offers a unique opportunity to tackle the genetic interactions between two contrasting reproductive modes. Here, we provide evidence that gene flow occurred between sexual and facultatively asexual lineages of R. padi. We carefully examined the shuffling in phenotypic and genotypic variation following a sexual reproduction event that took place in the field. Combining genotypic data and phenotypic measurements showed that this gene mixing led to the production of a wide array of reproductive modes, including strictly asexual lineages. Finally, we discuss the central role played by facultatively asexual lineages on the maintenance of reproductive mode variation.     

The genetic outcomes of sex and recombination in long-term functionally parthenogenetic lineages of Australian Sitobion aphids

The typical life cycle of an aphid is cyclical parthenogenesis which involves the alternation of sexual and asexual reproduction. However, aphid life cycles, even within a species, can encompass everything on a continuum from obligate sexuality, through facultative sexuality to obligate asexuality. Loss of the sexual cycle in aphids is frequently associated with the introduction of a new pest and can occur for a number of environmental and genetic reasons. Here we investigate loss of sexual function in Sitobion aphids in Australia. Specifically, we aimed to determine whether an absence of sexual reproduction in Australian Sitobion results from genetic loss of sexual function or environmental constraints in the introduced range. We addressed our aims by performing a series of breeding experiments. We found that some lineages have genetically lost sexual function while others retain sexual function and appear environmentally constrained to asexuality. Further, in our crosses, using autosomal and X-linked microsatellite markers, we identified processes deviating from normal Mendelian segregation. We observed strong deviations in X chromosome transmission through the sexual cycle. Additionally, when progeny genotypes were examined across multiple loci simultaneously we found that some multilocus genotypes are significantly over-represented in the sample and that levels of heterozygosity were much higher than expected at almost all loci. This study demonstrates that strong biases in the transmission of X chromosomes through the sexual cycle are likely to be widespread in aphids. The mechanisms underlying these patterns are not clear. We discuss several possible alternatives, including mutation accumulation during periods of functional asexuality and genetic imprinting.     

Evidence for gene flow and local clonal selection in field populations of the grain aphid (Sitobion avenae) in Britain revealed using microsatellites

Samples of the grain aphid, Sitobion avenae (F.), a major European pest of cereals, were collected in June and July 1997 from fields sown with winter wheat in a rough transect south-west of Rothamsted, UK. These aphids were genotyped at four microsatellite loci known from previous studies to be highly polymorphic. Allelic frequencies were similar between samples collected in the fields and in the 12.2 m high suction trap at Rothamsted, and there were many widespread genotypes (clones), providing evidence that the species is highly migratory. However, field samples were found to display a high level of genotypic heterogeneity (= variable clonal composition), most probably the result of clonal selection. The suction trap genotypes sample were slightly different from the field samples, indicative of the inclusion of genotypes from plant hosts (cereals and grasses, Poaceae) other than winter wheat and/or genotype-biased emigration from the field. The relevance of these data to modelling of aphid outbreaks is briefly discussed.     

Genetic Control of Contagious Asexuality in the Pea Aphid

Although evolutionary transitions from sexual to asexual reproduction are frequent in eukaryotes, the genetic bases of such shifts toward asexuality remain largely unknown. We addressed this issue in an aphid species where both sexual and obligate asexual lineages coexist in natural populations. These sexual and asexual lineages may occasionally interbreed because some asexual lineages maintain a residual production of males potentially able to mate with the females produced by sexual lineages. Hence, this species is an ideal model to study the genetic basis of the loss of sexual reproduction with quantitative genetic and population genomic approaches. Our analysis of the co-segregation of ∼300 molecular markers and reproductive phenotype in experimental crosses pinpointed an X-linked region controlling obligate asexuality, this state of character being recessive. A population genetic analysis (>400-marker genome scan) on wild sexual and asexual genotypes from geographically distant populations under divergent selection for reproductive strategies detected a strong signature of divergent selection in the genomic region identified by the experimental crosses. These population genetic data confirm the implication of the candidate region in the control of reproductive mode in wild populations originating from 700 km apart. Patterns of genetic differentiation along chromosomes suggest bidirectional gene flow between populations with distinct reproductive modes, supporting contagious asexuality as a prevailing route to permanent parthenogenesis in pea aphids. This genetic system provides new insights into the mechanisms of coexistence of sexual and asexual aphid lineages.     

Host range expansion of an introduced insect pest through multiple colonizations of specialized clones

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 effects of reproductive specialization on energy costs and fitness genetic variances in cyclical and obligate parthenogenetic aphids

Organisms with coexisting sexual and asexual populations are ideal models for studying the consequences of either reproductive mode on the quantitative genetic architecture of life-history traits. In the aphid Rhopalosiphum padi, lineages differing in their sex investment coexist but all share a common parthenogenetic phase. Here, we studied multiple genotypes of R. padi specialized either for sexual and asexual reproduction and compared their genetic variation in fitness during the parthenogenetic phase. Specifically, we estimated maintenance costs as standard metabolic rate (SMR), together with fitness (measured as the intrinsic rate of increase and the net reproductive rate). We found that genetic variation (in terms of broad-sense heritability) in fitness was higher in asexual genotypes compared with sexual genotypes. Also, we found that asexual genotypes exhibited several positive genetic correlations indicating that body mass, whole-animal SMR, and apterous individuals production are contributing to fitness. Hence, it appears that in asexual genotypes, energy is fully allocated to maximize the production of parthenogenetic individuals, the simplest possible form of aphid repertoire of life-histories strategies.     

Molecular markers linked to breeding system differences in segregating and natural populations of the cereal aphid Rhopalosiphum padi L.

The aphid Rhopalosiphum padi shows coexistence of sexual and asexual populations, providing an opportunity to study the evolution of breeding system variation in the context of theories on the origin and maintenance of sex. However, assessments of the distribution of sexual and asexual lineages of this aphid are complicated by the difficulties in rapidly characterizing their breeding system. To facilitate this task and to gain insight into the genetic relatedness between sexual and asexual genotypes, molecular markers linked to breeding system differences were recently developed. In this study, we have successfully converted a random amplified polymorphic DNA (RAPD) marker associated with life-cycle variation in R. padi into a codominant sequences-characterized amplified region (SCAR). Segregating and natural populations of known breeding systems were examined to evaluate the life cycle-SCAR marker association. Complete linkage was found in segregating populations while the association averaged 94% in field populations. Detailed analysis of allelic distribution and sequence divergence of the SCAR locus among sexual and asexual populations provides further evidence for a unique and apparently ancient loss of sexuality in this aphid. It also suggests that occasional gene flow occurs between populations differing in their breeding system, mediated by males produced by 'asexual' lineages. This system provides the possibility for the recurrent emergence of new asexual lineages, ensuring the longer persistence of asexuality, and would have important implications for the assessment of costs and benefits of sex in aphids.     

Facultative symbiont infections affect aphid reproduction

Some bacterial symbionts alter their hosts reproduction through various mechanisms that enhance their transmission in the host population. In addition to its obligatory symbiont Buchnera aphidicola, the pea aphid Acyrthosiphon pisum harbors several facultative symbionts influencing several aspects of host ecology. Aphids reproduce by cyclical parthenogenesis whereby clonal and sexual reproduction alternate within the annual life cycle. Many species, including the pea aphid, also show variation in their reproductive mode at the population level, with some lineages reproducing by cyclical parthenogenesis and others by permanent parthenogenesis. While the role of facultative symbionts has been well studied during the parthenogenetic phase of their aphid hosts, very little is known on their possible influence during the sexual phase. Here we investigated whether facultative symbionts modulate the capacity to produce sexual forms in various genetic backgrounds of the pea aphid with controlled symbiont composition and also in different aphid genotypes from natural populations with previously characterized infection status and reproductive mode. We found that most facultative symbionts exhibited detrimental effects on their hosts fitness under sex-inducing conditions in comparison with the reference lines. We also showed that the loss of sexual phase in permanently parthenogenetic lineages of A. pisum was not explained by facultative symbionts. Finally, we demonstrated that Spiroplasma infection annihilated the production of males in the host progeny by inducing a male-killing phenotype, an unexpected result for organisms such as aphids that reproduce primarily through clonal reproduction.