Genomic insights into Spiroplasma endosymbionts that induce male-killing and protective phenotypes in the pea aphid.

The endosymbiotic bacteria Spiroplasma (Mollicutes) infect diverse plants and arthropods, and some of which induce male killing, where male hosts are killed during development. Male-killing Spiroplasma strains belong to either the phylogenetically distant Citri-Poulsonii or Ixodetis groups. In Drosophila flies, Spiroplasma poulsonii induces male killing via the Spaid toxin. While Spiroplasma ixodetis infects a wide range of insects and arachnids, little is known about the genetic basis of S. ixodetis-induced male killing. Here, we analyzed the genome of S. ixodetis strains in the pea aphid Acyrthosiphon pisum (Aphididae, Hemiptera). Genome sequencing constructed a complete genome of a male-killing strain, sAp269, consisting of a 1.5 Mb circular chromosome and an 80 Kb plasmid. sAp269 encoded putative virulence factors containing either ankyrin repeat, ovarian tumor-like deubiquitinase, or ribosome inactivating protein domains, but lacked the Spaid toxin. Further comparative genomics of Spiroplasma strains in A. pisum biotypes adapted to different host plants revealed their phylogenetic associations and the diversity of putative virulence factors. Although the mechanisms of S. ixodetis-induced male killing in pea aphids remain elusive, this study underlines the dynamic genome evolution of S. ixodetis and proposes independent acquisition events of male-killing mechanisms in insects.

Rapid turnover of a pea aphid superclone mediated by thermal endurance in central Chile.

Global change drivers are imposing novel conditions on Earth's ecosystems at an unprecedented rate. Among them, biological invasions and climate change are of critical concern. It is generally thought that strictly asexual populations will be more susceptible to rapid environmental alterations due to their lack of genetic variability and, thus, of adaptive responses. In this study, we evaluated the persistence of a widely distributed asexual lineage of the alfalfa race of the pea aphid, Acyrthosiphon pisum, along a latitudinal transect of approximately 600 km in central Chile after facing environmental change for a decade. Based on microsatellite markers, we found an almost total replacement of the original aphid superclone by a new variant. Considering the unprecedented warming that this region has experienced in recent years, we experimentally evaluated the reproductive performance of these two A. pisum lineages at different thermal regimes. The new variant exhibits higher rates of population increase at warmer temperatures, and computer simulations employing a representative temperature dataset suggest that it might competitively displace the original superclone. These results support the idea of a superclone turnover mediated by differential reproductive performance under changing temperatures.

Contrasting Evolutionary Patterns Between Sexual and Asexual Lineages in a Genomic Region Linked to Reproductive Mode Variation in the pea aphid.

Although asexual lineages 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. An SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from 6 populations) revealed that an X-linked region is associated with the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2,000 SNPs or short Indels showing a high degree of association with the phenotypic trait. In OP populations specifically, this region also shows reduced diversity and Tajima's D, consistent with the OP phenotype being a derived trait in aphids. Interestingly, the low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.5%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. These genetic exchanges, combined with the selection of the OP and CP reproductive modes under different climates, probably contribute to the long-term persistence of the cp and op alleles.

Evolutionary consequences of loss of sexual reproduction on male-related traits in parthenogenetic lineages of the pea aphid.

Transition from sexual reproduction to parthenogenesis constitutes a major life-history change with deep evolutionary consequences for sex-related traits, which are expected to decay. The pea aphid Acyrthosiphon pisum shows intraspecific reproductive polymorphism, with cold-resistant cyclically parthenogenetic (CP) lineages that alternate sexual and asexual generations and cold-sensitive obligately parthenogenetic (OP) lineages that produce only asexual females but still males. Here, the genotyping of 219 pea aphid lineages collected in cold-winter and mild-winter regions revealed contrasting population structures. Samples from cold-winter regions consisted mostly of distinct multilocus genotypes (MLGs) usually represented by a single sample (101 different MLGs for 111 samples) and were all phenotyped as CP. In contrast, fewer MLGs were found in mild-winter regions (28 MLGs for 108 samples), all but one being OP. Since the males produced by OP lineages are unlikely to pass on their genes (sexual females being rare in mild-winter regions), we tested the hypothesis that their traits could degenerate due to lack of selection by comparing male production and male reproductive success between OP and CP lineages. Male production was indeed reduced in OP lineages, but a less clear pattern was observed for male reproductive success: females mated with OP males laid fewer eggs (fertilized or not) but OP and CP males fertilized the same proportion of eggs. These differences may stem from the type of selective forces: male production may be counter-selected whereas male performances may evolve under the slower process of relaxed selection. The overall effective reproductive capacity of OP males could result from recent sex loss in OP lineages or underestimated reproductive opportunities.

Molecular Mechanisms Underlying Host Plant Specificity in Aphids.

Aphids are serious pests of agricultural and ornamental plants and important model systems for hemipteran-plant interactions. The long evolutionary history of aphids with their host plants has resulted in a variety of systems that provide insight into the different adaptation strategies of aphids to plants and vice versa. In the past, various plant-aphid interactions have been documented, but lack of functional tools has limited molecular studies on the mechanisms of plant-aphid interactions. Recent technological advances have begun to reveal plant-aphid interactions at the molecular level and to increase our knowledge of the mechanisms of aphid adaptation or specialization to different host plants. In this article, we compile and analyze available information on plant-aphid interactions, discuss the limitations of current knowledge, and argue for new research directions. We advocate for more work that takes advantage of natural systems and recently established molecular techniques to obtain a comprehensive view of plant-aphid interaction mechanisms.

Host plants and insecticides shape the evolution of genetic and clonal diversity in a major aphid crop pest.

Understanding the spatiotemporal dynamics of pesticide resistance at the landscape scale is essential to anticipate the evolution and spread of new resistance phenotypes. In crop mosaics, host plant specialization in pest populations is likely to dampen the spread of pesticide resistance between different crops even in mobile pests such as aphids. Here, we assessed the contribution of host-based genetic differentiation to the dynamics of resistance alleles in Myzus persicae, a major aphid pest which displays several insecticide resistance mechanisms. We obtained a representative sample of aphids from a crop mosaic through a suction trap for 7 years and from various crops as a reference collection. We genotyped these aphids at 14 microsatellite markers and four insecticide-resistant loci, analyzed the genetic structure, and assigned host-based genetic groups from field-collected aphids. Four well-defined genetic clusters were found in aerial samples, three of which with strong association with host-plants. The fourth group was exclusive to aerial samples and highly divergent from the others, suggesting mixture with a closely related taxon of M. persicae associated with unsampled plants. We found a sharp differentiation between individuals from peach and herbaceous plants. Individuals from herbaceous hosts were separated into two genetic clusters, one more strongly associated with tobacco. The 4-loci resistance genotypes showed a strong association with the four genetic clusters, indicative of barriers to the spread of resistance. However, we found a small number of clones with resistant alleles on multiple host-plant species, which may spread insecticide resistance between crops. The 7-year survey revealed a rapid turn-over of aphid genotypes as well as the emergence, frequency increase and persistence of clones with resistance to several families of insecticides. This study highlights the importance of considering landscape-scale population structure to identify the risk of emergence and spread of insecticide resistance for a particular crop.

Aphid Resistance in Pisum Affects the Feeding Behavior of Pea-Adapted and Non-Pea-Adapted Biotypes of Acyrthosiphon pisum Differently.

Resistant genotypes of crops have emerged as an alternative and sustainable solution to pesticide use against pest insects. The resistance depends on the genetic diversity of the host plant and the pest species and can cause an alteration of the insect behavior. The aim of this work was to characterize the resistance level of different Pisum genotypes (one P. fulvum and five P. sativum genotypes) to two biotypes of the aphid Acyrthosiphon pisum, respectively adapted to pea and alfalfa, by measuring the individual aphid weight and analyzing aphid feeding behavior by electropenetrography (EPG). Aphid body mass was influenced by Pisum genotypes reflecting variation in their resistance level. P. fulvum was the most resistant to the A. pisum pea biotype (ArPo28 clone) and showed intermediate resistance to the A. pisum alfalfa biotype (LSR1 clone). The resistance levels of the five P. sativum genotypes to the two aphid biotypes were variable and more pronounced for the alfalfa biotype. EPG data showed that ArPo28 on P. fulvum and LSR1 on all the Pisum genotypes spent shorter time phloem feeding compared to ArPo28 on P. sativum genotypes, indicating that the resistance of Pisum genotypes to non-adapted A. pisum resides in mesophyll and phloem cells. In the meantime, ArPo28 on P. sativum genotypes with a different level of resistance spent a similar length of time phloem feeding, indicating that the quality of phloem sap of the resistance genotypes may not be optimal for the aphid. The study indicated that the resistance of Pisum genotypes to the two A. pisum biotypes involves different genetic factors and mechanisms that affect the aphid differently.

A major-effect genetic locus, ApRVII, controlling resistance against both adapted and non-adapted aphid biotypes in pea.

KEY MESSAGE: A genome-wide association study for pea resistance against a pea-adapted biotype and a non-adapted biotype of the aphid, Acyrthosiphon pisum, identified a genomic region conferring resistance to both biotypes. In a context of reduced insecticide use, the development of cultivars resistant to insect pests is crucial for an integrated pest management. Pea (Pisum sativum) is a crop of major importance among cultivated legumes, for the supply of dietary proteins and nitrogen in low-input cropping systems. However, yields of the pea crop have become unstable due to plant parasites. The pea aphid (Acyrthosiphon pisum) is an insect pest species forming a complex of biotypes, each one adapted to feed on one or a few related legume species. This study aimed to identify resistance to A. pisum and the underlying genetic determinism by examining a collection of 240 pea genotypes. The collection was screened against a pea-adapted biotype and a non-adapted biotype of A. pisum to characterize their resistant phenotype. Partial resistance was observed in some pea genotypes exposed to the pea-adapted biotype. Many pea genotypes were completely resistant to non-adapted biotype, but some exhibited partial susceptibility. A genome-wide association study, using pea exome-capture sequencing data, enabled the identification of the major-effect quantitative trait locus ApRVII on the chromosome 7. ApRVII includes linkage disequilibrium blocks significantly associated with resistance to one or both of the two aphid biotypes studied. Finally, we identified candidate genes underlying ApRVII that are potentially involved in plant-aphid interactions and marker haplotypes linked with aphid resistance. This study sets the ground for the functional characterization of molecular pathways involved in pea defence to the aphids but also is a step forward for breeding aphid-resistant cultivars.

Influence of "protective" symbionts throughout the different steps of an aphid-parasitoid interaction.

Microbial associates are widespread in insects, some conferring a protection to their hosts against natural enemies like parasitoids. These protective symbionts may affect the infection success of the parasitoid by modifying behavioral defenses of their hosts, the development success of the parasitoid by conferring a resistance against it or by altering life-history traits of the emerging parasitoids. Here, we assessed the effects of different protective bacterial symbionts on the entire sequence of the host-parasitoid interaction (i.e., from parasitoid attack to offspring emergence) between the pea aphid, Acyrthosiphon pisum, and its main parasitoid, Aphidius ervi and their impacts on the life-history traits of the emerging parasitoids. To test whether symbiont-mediated phenotypes were general or specific to particular aphid-symbiont associations, we considered several aphid lineages, each harboring a different strain of either Hamiltonella defensa or Regiella insecticola, two protective symbionts commonly found in aphids. We found that symbiont species and strains had a weak effect on the ability of aphids to defend themselves against the parasitic wasps during the attack and a strong effect on aphid resistance against parasitoid development. While parasitism resistance was mainly determined by symbionts, their effects on host defensive behaviors varied largely from one aphid-symbiont association to another. Also, the symbiotic status of the aphid individuals had no impact on the attack rate of the parasitic wasps, the parasitoid emergence rate from parasitized aphids nor the life-history traits of the emerging parasitoids. Overall, no correlations between symbiont effects on the different stages of the host-parasitoid interaction was observed, suggesting no trade-offs or positive associations between symbiont-mediated phenotypes. Our study highlights the need to consider various sequences of the host-parasitoid interaction to better assess the outcomes of protective symbioses and understand the ecological and evolutionary dynamics of insect-symbiont associations.

Global patterns in genomic diversity underpinning the evolution of insecticide resistance in the aphid crop pest Myzus persicae.

The aphid Myzus persicae is a destructive agricultural pest that displays an exceptional ability to develop resistance to both natural and synthetic insecticides. To investigate the evolution of resistance in this species we generated a chromosome-scale genome assembly and living panel of >110 fully sequenced globally sampled clonal lines. Our analyses reveal a remarkable diversity of resistance mutations segregating in global populations of M. persicae. We show that the emergence and spread of these mechanisms is influenced by host-plant associations, uncovering the widespread co-option of a host-plant adaptation that also offers resistance against synthetic insecticides. We identify both the repeated evolution of independent resistance mutations at the same locus, and multiple instances of the evolution of novel resistance mechanisms against key insecticides. Our findings provide fundamental insights into the genomic responses of global insect populations to strong selective forces, and hold practical relevance for the control of pests and parasites.

Secondary Symbionts Affect Foraging Capacities of Plant-Specialized Genotypes of the Pea Aphid.

Ecological specialization is widespread in animals, especially in phytophagous insects, which have often a limited range of host plant species. This host plant specialization results from divergent selection on insect populations, which differ consequently in traits like behaviors involved in plant use. Although recent studies highlighted the influence of symbionts on dietary breadth of their insect hosts, whether these microbial partners influence the foraging capacities of plant-specialized insects has received little attention. In this study, we used the pea aphid Acyrthosiphon pisum, which presents distinct plant-specialized lineages and several secondary bacterial symbionts, to examine the possible effects of symbionts on the different foraging steps from plant searching to host plant selection. In particular, we tested the effect of secondary symbionts on the aphid capacity (1) to explore habitat at long distance (estimated through the production of winged offspring), (2) to explore habitat at short distance, and (3) to select its host plant. We found that secondary symbionts had a variable influence on the production of winged offspring in some genotypes, with potential consequences on dispersal and survival. By contrast, symbionts influenced both short-distance exploration and host plant selection only marginally. The implication of symbionts' influence on insect foraging capacities is discussed.

Alarm Pheromone Responses Depend on Genotype, but Not on the Presence of Facultative Endosymbionts in the Pea Aphid Acyrthosiphon pisum.

Aphids use an alarm pheromone, E-ß farnesene (EBF), to warn conspecifics of potential danger. The antennal sensitivity and behavioural escape responses to EBF can be influenced by different factors. In the pea aphid, Acyrthosiphon pisum, different biotypes are adapted to different legume species, and within each biotype, different genotypes exist, which can carry or not Hamiltonella defensa, a bacterial symbiont that can confer protection against natural enemies. We investigate here the influence of the aphid genotype and symbiotic status on the escape behaviour using a four-way olfactometer and antennal sensitivity for EBF using electroantennograms (EAGs). Whereas the investigated three genotypes from two biotypes showed significantly different escape and locomotor behaviours in the presence of certain EBF doses, the infection with H. defensa did not significantly modify the escape behaviour and only marginally influenced the locomotor behaviour at high doses of EBF. Dose-response curves of EAG amplitudes after stimulation with EBF differed significantly between aphid genotypes in correlation with behavioural differences, whereas antennal sensitivity to EBF did not change significantly as a function of the symbiotic status. The protective symbiont H. defensa does thus not modify the olfactory sensitivity to the alarm pheromone. How EBF sensitivity is modified between genotypes or biotypes remains to be investigated.

Involvement of chemosensory proteins in host plant searching in the bird cherry-oat aphid.

Chemosensory systems are considered to play an important role in host plant selection in herbivorous insects. However, few studies have focused on chemosensory proteins (CSPs) for aphid host-location mechanisms. The roles of CSPs in searching for different Poaceae species (wheat, barley, triticale, maize and sorghum) were tested in Rhopalosiphum padi, an important cereal pest. The olfactometer assays showed that R. padi responds to plant odors. Seven R. padi CSP genes were identified. Influence of aphid morph, tissue and starvation state on expression patterns of CSPs was evaluated. Expression levels of CSP1, CSP4, CSP5 and CSP6 in winged aphids were significantly higher than those in wingless ones. Transcription levels of four genes (CSP1, CSP4, CSP5 and CSP6) were relatively higher in the head with antennae, and the four genes tended to be upregulated following starvation. Silencing of three CSPs (CSP4, CSP5 and CSP6) altered aphid host-location behavior in response to the five different host plants tested. Three volatile compounds of host plants (octanal, [E]-2-hexenol and linalool) have significant attraction to winged R. padi according to the four-arm olfactometer tests. Molecular docking predicted hydrogen bonding sites which played key roles in the binding of CSP4, CSP5 and CSP6 with volatile compounds. Knockdown of CSP4 or CSP5 significantly decreased the staying time of R. padi in the arms with octanal. However, knockdown of CSP6 could not affect the response of R. padi to octanal. These results bring evidence for the involvement of three CSPs in R. padi host-location behavior.

Isolation of Insect Bacteriocytes as a Platform for Transcriptomic Analyses.

Over the past few decades, various techniques have been developed and optimized for the accurate measurement of RNA abundance in cells or tissues. These methods have been instrumental in gaining insight in complex systems such as host-symbiont associations. The pea aphid model has recently emerged as a powerful and experimentally tractable system for studying symbiotic relationships and it is the subject of a growing number of molecular studies. Nevertheless, the lack of standardized protocols for the collection of bacteriocytes, the specialized host cells harboring the symbionts, has limited its use. This chapter provides a simple, step-by-step dissection protocol for the rapid isolation of aphid bacteriocytes. We then describe an adapted protocol for efficient extraction and purification of bacteriocyte RNA that can be used for most downstream transcriptomic analyses.

Hosting certain facultative symbionts modulates the phenoloxidase activity and immune response of the pea aphid Acyrthosiphon pisum.

The pea aphid Acyrthosiphon pisum hosts different facultative symbionts (FS) which provide it with various benefits, such as tolerance to heat or protection against natural enemies (e.g., fungi, parasitoid wasps). Here, we investigated whether and how the presence of certain FS could affect phenoloxidase (PO) activity, a key component of insect innate immunity, under normal and stressed conditions. For this, we used clones of A. pisum of different genetic backgrounds (LL01, YR2 and T3-8V1) lacking FS or harboring one or two (Regiella insecticola, Hamiltonella defensa, Serratia symbiotica + Rickettsiella viridis). Gene expression and proteomics analyses of the aphid hemolymph indicated that the two A. pisum POs, PPO1 and PPO2, are expressed and translated into proteins. The level of PPO genes expression as well as the amount of PPO proteins and phenoloxidase activity in the hemolymph depended on both the aphid genotype and FS species. In particular, H. defensa and R. insecticola, but not S. symbiotica + R. viridis, caused a sharp decrease in PO activity by interfering with both transcription and translation. The microinjection of different types of stressors (yeast, Escherichia coli, latex beads) in the YR2 lines hosting different symbionts affected the survival rate of aphids and, in most cases, also decreased the expression of PPO genes after 24 h. The amount and activity of PPO proteins varied according to the type of FS and stressor, without clear corresponding changes in gene expression. These data demonstrate that the presence of certain FS influences an important component of pea aphid immunity.

Aphid infestation differently affects the defences of nitrate-fed and nitrogen-fixing Medicago truncatula and alters symbiotic nitrogen fixation.

Legumes can meet their nitrogen requirements through root nodule symbiosis, which could also trigger plant systemic resistance against pests. The pea aphid Acyrthosiphon pisum, a legume pest, can harbour different facultative symbionts (FS) influencing various traits of their hosts. It is therefore worth determining if and how the symbionts of the plant and the aphid modulate their interaction. We used different pea aphid lines without FS or with a single one (Hamiltonella defensa, Regiella insecticola, Serratia symbiotica) to infest Medicago truncatula plants inoculated with Sinorhizobium meliloti (symbiotic nitrogen fixation, SNF) or supplemented with nitrate (non-inoculated, NI). The growth of SNF and NI plants was reduced by aphid infestation, while aphid weight (but not survival) was lowered on SNF compared to NI plants. Aphids strongly affected the plant nitrogen fixation depending on their symbiotic status, suggesting indirect relationships between aphid- and plant-associated microbes. Finally, all aphid lines triggered expression of Pathogenesis-Related Protein 1 (PR1) and Proteinase Inhibitor (PI), respective markers for salicylic and jasmonic pathways, in SNF plants, compared to only PR1 in NI plants. We demonstrate that the plant symbiotic status influences plant-aphid interactions while that of the aphid can modulate the amplitude of the plant's defence response.

Latitudinal trend in the reproductive mode of the pea aphid Acyrthosiphon pisum invading a wide climatic range.

The maintenance of sexuality is a puzzling phenomenon in evolutionary biology. Many universal hypotheses have been proposed to explain the prevalence of sex despite its costs, but it has been hypothesized that sex could be also retained by lineage-specific mechanisms that would confer some short-term advantage. Aphids are good models to study the maintenance of sex because they exhibit coexistence of both sexual and asexual populations within the same species and because they invade a large variety of ecosystems. Sex in aphids is thought to be maintained because only sexually produced eggs can persist in cold climates, but whether sex is obligate or facultative depending on climatic conditions remains to be elucidated. In this study, we have inferred the reproductive mode of introduced populations of the pea aphid Acyrthosiphon pisum in Chile along a climatic gradient using phenotypic assays and genetic-based criteria to test the ecological short-term advantage of sex in cold environments. Our results showed a latitudinal trend in the reproductive mode of Chilean pea aphid population from obligate parthenogenesis in the north to an intermediate life cycle producing both parthenogenetic and sexual progeny in the southernmost locality, where harsh winters are usual. These findings are congruent with the hypothesis of the ecological short-term advantage of sex in aphids.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum.

BACKGROUND: Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts. RESULTS: We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes. CONCLUSIONS: These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Intraspecific difference among herbivore lineages and their host-plant specialization drive the strength of trophic cascades.

Trophic cascades - the indirect effect of predators on non-adjacent lower trophic levels - are important drivers of the structure and dynamics of ecological communities. However, the influence of intraspecific trait variation on the strength of trophic cascade remains largely unexplored, which limits our understanding of the mechanisms underlying ecological networks. Here we experimentally investigated how intraspecific difference among herbivore lineages specialized on different host plants influences trophic cascade strength in a terrestrial tri-trophic system. We found that the occurrence and strength of the trophic cascade are strongly influenced by herbivores' lineage and host-plant specialization but are not associated with density-dependent effects mediated by the growth rate of herbivore populations. Our findings stress the importance of intraspecific heterogeneities and evolutionary specialization as drivers of trophic cascade strength and underline that intraspecific variation should not be overlooked to decipher the joint influence of evolutionary and ecological factors on the functioning of multi-trophic interactions.

A large genomic insertion containing a duplicated follistatin gene is linked to the pea aphid male wing dimorphism.

Wing dimorphisms have long served as models for examining the ecological and evolutionary tradeoffs associated with alternative phenotypes. Here, we investigated the genetic cause of the pea aphid (Acyrthosiphon pisum) male wing dimorphism, wherein males exhibit one of two morphologies that differ in correlated traits that include the presence or absence of wings. We mapped this trait difference to a single genomic region and, using third generation, long-read sequencing, we identified a 120 kb insertion in the wingless allele. This insertion includes a duplicated follistatin gene, which is a strong candidate gene in the minimal mapped interval to cause the dimorphism. We found that both alleles were present prior to pea aphid biotype lineage diversification, we estimated that the insertion occurred millions of years ago, and we propose that both alleles have been maintained in the species, likely due to balancing selection.