Post-Pleistocene radiation of the pea aphid complex revealed by rapidly evolving endosymbionts

'Drifting' Buchnera genomes track the microevolutionary trajectories of their aphid hosts

Evolution of Buchnera-aphid host symbioses is often studied among species at macroevolutionary scales. Investigations within species offer a different perspective about how eco-evolutionary processes shape patterns of genetic variation at microevolutionary scales. Our study leverages new and publicly available whole-genome sequencing data to study Buchnera-aphid host evolution in Myzus persicae, the peach potato aphid, a globally invasive and polyphagous pest. Across 43 different asexual, clonally reproducing isofemale strains, we examined patterns of genomic covariation between Buchnera and their aphid host and considered the distribution of mutations in protein-coding regions of the Buchnera genome. We found Buchnera polymorphisms within aphid strains, suggesting the presence of genetically different Buchnera strains within the same clonal lineage. Genetic distance between pairs of Buchnera samples was positively correlated to genetic distance between their aphid hosts, indicating shared evolutionary histories. However, there was no segregation of genetic variation for both M. persicae and Buchnera with plant host (Brassicaceae and non-tobacco Solanaceae) and no associations between genetic and geographic distance at global or regional spatial scales. Abundance patterns of non-synonymous mutations were similar to synonymous mutations in the Buchnera genome, and both mutation classes had similar site frequency spectra. We hypothesize that a predominance of neutral processes results in the Buchnera of M. persicae to simply 'drift' with the evolutionary trajectory of their aphid hosts. Our study presents a unique microevolutionary characterization of Buchnera-aphid host genomic covariation across multiple aphid clones. This provides a new perspective on the eco-evolutionary processes generating and maintaining polymorphisms in a major pest aphid species and its obligate primary endosymbiont.

Systematics and phylogeny of the entomopathogenic nematobacterial complexes Steinernema-Xenorhabdus and Heterorhabditis-Photorhabdus

Entomopathogenic nematodes of the genera Steinernema and Heterorhabditis, along with their bacterial symbionts from the genera Xenorhabdus and Photorhabdus, respectively, are important biological control agents against agricultural pests. Rapid progress in the development of genomic tools has catalyzed a transformation of the systematics of these organisms, reshaping our understanding of their phylogenetic and cophlylogenetic relationships. In this review, we discuss the major historical events in the taxonomy and systematics of this group of organisms, highlighting the latest advancements in these fields. Additionally, we synthesize information on nematode-bacteria associations and assess the existing evidence regarding their cophylogenetic relationships.

Comparative genomics of the primary endosymbiont Buchnera aphidicola in aphid hosts and their coevolutionary relationships

BACKGROUND

Coevolution between modern aphids and their primary obligate, bacterial endosymbiont, Buchnera aphidicola, has been previously reported at different classification levels based on molecular phylogenetic analyses. However, the Buchnera genome remains poorly understood within the Rhus gall aphids.

RESULTS

We assembled the complete genome of the endosymbiont Buchnera in 16 aphid samples, representing 13 species in all six genera of Rhus gall aphids by shotgun genome skimming method. We compared the newly assembled genomes with those from GenBank to comprehensively investigate patterns of coevolution between the bacteria Buchnera and their aphid hosts. Buchnera genomes were mostly collinear, and the pan-genome contained 684 genes, in which the core genome contained 256 genes with some lineages having large numbers of tandem gene duplications. There has been substantial gene-loss in each Buchnera lineage. We also reconstructed the phylogeny for Buchnera and their host aphids, respectively, using 72 complete genomes of Buchnera, along with the complete mitochondrial genomes and three nuclear genes of 31 corresponding host aphid accessions. The cophylogenetic test demonstrated significant coevolution between these two partner groups at individual, species, generic, and tribal levels.

CONCLUSIONS

Buchnera exhibits very high levels of genomic sequence divergence but relative stability in gene order. The relationship between the symbionts Buchnera and its aphid hosts shows a significant coevolutionary pattern and supports complexity of the obligate symbiotic relationship.

Mobile elements create strain-level variation in the services conferred by an aphid symbiont

Heritable, facultative symbionts are common in arthropods, often functioning in host defence. Despite moderately reduced genomes, facultative symbionts retain evolutionary potential through mobile genetic elements (MGEs). MGEs form the primary basis of strain-level variation in genome content and architecture, and often correlate with variability in symbiont-mediated phenotypes. In pea aphids (Acyrthosiphon pisum), strain-level variation in the type of toxin-encoding bacteriophages (APSEs) carried by the bacterium Hamiltonella defensa correlates with strength of defence against parasitoids. However, co-inheritance creates difficulties for partitioning their relative contributions to aphid defence. Here we identified isolates of H. defensa that were nearly identical except for APSE type. When holding H. defensa genotype constant, protection levels corresponded to APSE virulence module type. Results further indicated that APSEs move repeatedly within some H. defensa clades providing a mechanism for rapid evolution in anti-parasitoid defences. Strain variation in H. defensa also correlates with the presence of a second symbiont Fukatsuia symbiotica. Predictions that nutritional interactions structured this coinfection were not supported by comparative genomics, but bacteriocin-containing plasmids unique to co-infecting strains may contribute to their common pairing. In conclusion, strain diversity, and joint capacities for horizontal transfer of MGEs and symbionts, are emergent players in the rapid evolution of arthropods.

Pea Aphid (Acyrthosiphon pisum) Host Races Reduce Heat-Induced Forisome Dispersion in Vicia faba and Trifolium pratense

Although phloem-feeding insects such as aphids can cause significant damage to plants, relatively little is known about early plant defenses against these insects. As a first line of defense, legumes can stop the phloem mass flow through a conformational change in phloem proteins known as forisomes in response to Ca²⁺ influx. However, specialized phloem-feeding insects might be able to suppress the conformational change of forisomes and thereby prevent sieve element occlusion. To investigate this possibility, we triggered forisome dispersion through application of a local heat stimulus to the leaf tips of pea (Pisum sativum), clover (Trifolium pratense) and broad bean (Vicia faba) plants infested with different pea aphid (Acyrthosiphon pisum) host races and monitored forisome responses. Pea aphids were able to suppress forisome dispersion, but this depended on the infesting aphid host race, the plant species, and the age of the plant. Differences in the ability of aphids to suppress forisome dispersion may be explained by differences in the composition and quantity of the aphid saliva injected into the plant. Various mechanisms of how pea aphids might suppress forisome dispersion are discussed.

Aphid male wing polymorphisms are transient and have evolved repeatedly

Polymorphic phenotypes have long been used to examine the maintenance of genetic variation within and between species. Most studies have focused on persistent polymorphisms, which are retained across species boundaries, and their positive effects on speciation rates. Far less is known about the macroevolutionary impacts of more transient polymorphisms, which are also common. Here we investigated male wing polymorphisms in aphids. We estimated the phylogenetic history of wing states across species, along with several other traits that could affect wing evolution. We found that male wing polymorphisms are transient: they are found in only ~4% of extant species, but have likely evolved repeatedly across the phylogeny. We reason that the repeated evolution of transient polymorphisms might be facilitated by the existence of the asexual female wing plasticity, which is common across aphids, and would maintain the wing development program even in species with wingless males. We also discovered that male wingedness correlates positively with host plant alternation and host plant breadth, and that winged morphs and wing polymorphisms may be associated with higher speciation rates. Our results provide new evolutionary insights into this well-studied group and suggest that even transient polymorphisms may impact species diversification rates.

Climatic oscillation promoted diversification of spinous assassin bugs during Pleistocene glaciation

Insect speciation is among the most fascinating topics in evolutionary biology; however, its underlying mechanisms remain unclear. Allopatric speciation represents one of the major types of speciation and is believed to have frequently occurred during glaciation periods, when climatic oscillation may have caused suitable habitats to be fragmented repeatedly, creating geographical isolation among populations. However, supporting evidence for allopatric speciation of insects in East Asia during the Pleistocene glaciation remains lacking. We aim to investigate the effect of climatic oscillation during the Pleistocene glaciation on the diversification pattern and evolutionary history of hemipteran insects and to test the hypothesis of Pleistocene species stability using spinous assassin bugs Sclomina (Hemiptera: Reduviidae), a small genus widely distributed in southern China but was later found to have cryptic species diversity. Here, using the whole mitochondrial genome (mitogenome) and nuclear ribosomal RNA genes, we investigated both interspecific and intraspecific diversification patterns of spinous assassin bugs. Approximate Bayesian computation, ecological niche modeling, and demographic history analyses were also applied to understand the diversification process and driven factors. Our data suggest that the five species of Sclomina are highly diverged, despite three of them currently being cryptic. Speciation occurred during the Pleistocene when suitable distribution areas were possibly fragmented. Six phylogeographic groups in the type species S. erinacea were identified, among which two groups underwent expansion during the early Last Glacial Period and after Last Glacier Maximum. Our analyses suggest that this genus may have experienced climate-driven habitat fragmentation and postglacial expansion in the Pleistocene, promoting allopatric speciation and intraspecific diversification. Our results reveal underestimated species diversity in a small insect group and illustrate a remarkable example of allopatric speciation of insects in East Asia promoted by Pleistocene climatic oscillations. These findings provide important insights into the speciation processes and aid the conservation of insect species diversity.

Generalism in nature: a community ecology perspective

Life on Earth is complex and generally abounds in food webs with other living organisms in terms of an ecological community. Besides such complexity, and the fact that populations of most living organisms have never been studied in terms of their molecular ecology, it is best to tread carefully when describing a given species as a ‘generalist’, more especially in terms of dietary and habitat breadth. We very much doubt that population homogeneity ever exists—because populations are always undergoing molecular-genetic changes, sometimes rapid, in response to various ecological challenges (e.g. climate, intra- and interspecific competition). In any case, a population may already have begun to undergo cryptic speciation. Such entities can occupy different habitats or exhibit different dietary breadths as a result of various ecological interactions formed over different spatial scales. These scales include everything from local (including islands) to geographic. The fossil evidence reveals that specialisations have existed over vast swathes of time. Besides, as is well documented, evolution only occurs as a result of adaptations leading to specialisation, and ultimately, specialist entitles, i.e. species and lower levels of ecological-evolutionary divergence. Here, focusing on diet, we posit that the terms mono-, oligo-and polyphagous are more accurate in relation to the dietary breadth of animals, with omnivory adopted in the case of organisms with very different food items. Thus, we strongly urge that the dubious and unscientific term ‘generalism’ be dropped in favour of these more precise and scientifically accurate terms directly relating to levels of phagy.

The phylogenetic and evolutionary analyses of detoxification gene families in Aphidinae species

Detoxification enzymes play significant roles in the interactions between insects and host plants, wherein detoxification-related genes make great contributions. As herbivorous pests, aphids reproduce rapidly due to parthenogenesis. They are good biological materials for studying the mechanisms that allow insect adaptation to host plants. Insect detoxification gene families are associated with insect adaptation to host plants. The Aphidinae is the largest subfamily in the Aphididae with at least 2483 species in 256 genera in 2 tribes: the Macrosiphini (with 3/4 of the species) and the Aphidini. Most aphid pests on crops and ornamental plants are Aphidinae. Members of the Aphidinae occur in nearly every region of the world. The body shape and colour vary significantly. To research the role that detoxification gene families played in the process of aphid adaptation to host evolution, we analyzed the phylogeny and evolution of these detoxification gene families in Aphidinae. In general, the P450/GST/CCE gene families contract, whereas the ABC/UGT families are conserved in Aphidinae species compared to these families in other herbivorous insects. Genus-specific expansions of P450 CYP4, and GST Delta have occurred in the genus Acyrthosiphon. In addition, the evolutionary rates of five detoxification gene families in the evolution process of Aphidinae are different. The comparison of five detoxification gene families among nine Aphidinae species and the estimated relative evolutionary rates provided herein support an understanding of the interaction between and the co-evolution of Aphidinae and plants.

Defensive Symbionts and the Evolution of Parasitoid Host Specialization

Insect host-parasitoid interactions abound in nature and are characterized by a high degree of host specialization. In addition to their behavioral and immune defenses, many host species rely on heritable bacterial endosymbionts for defense against parasitoids. Studies on aphids and flies show that resistance conferred by symbionts can be very strong and highly specific, possibly as a result of variation in symbiont-produced toxins. I argue that defensive symbionts are therefore an important source of diversifying selection, promoting the evolution of host specialization by parasitoids. This is likely to affect the structure of host-parasitoid food webs. I consider potential changes in terms of food web complexity, although the nature of these effects will also be influenced by whether maternally transmitted symbionts have some capacity for lateral transfer. This is discussed in the light of available evidence for horizontal transmission routes. Finally, I propose that defensive mutualisms other than microbial endosymbionts may also exert diversifying selection on insect parasitoids.

Natural selection shapes maintenance of orthologous sRNAs in divergent host-restricted bacterial genomes

Historically it has been difficult to study the evolution of bacterial small RNAs (sRNAs) across distantly related species. For example, identifying homologs of sRNAs is often difficult in genomes that have undergone multiple structural rearrangements. Also, some types of regulatory sRNAs evolve at rapid rates. The high degree of genomic synteny among divergent host-restricted bacterial lineages, including intracellular symbionts, is conducive to sRNA maintenance and homolog identification. In turn, symbiont genomes can provide us with novel insights into sRNA evolution. Here, we examine the sRNA expression profile of the obligate symbiont of psyllids, Carsonella ruddii, which has one of the smallest cellular genomes described. Using RNA-seq, we identified 36 and 32 antisense sRNAs (asRNAs) expressed by Carsonella from the psyllids Bactericera cockerelli (Carsonella-BC) and Diaphorina citri (Carsonella-DC), respectively. The majority of these asRNAs were associated with genes that are involved in essential amino acid biosynthetic pathways. Eleven of the asRNAs were conserved in both Carsonella lineages and the majority were maintained by selection. Notably, five of the corresponding coding sequences are also the targets of conserved asRNAs in a distantly related insect symbiont, Buchnera. We detected differential expression of two asRNAs for genes involved in arginine and leucine biosynthesis occurring between two distinct Carsonella-BC life stages. Using asRNAs identified in Carsonella, Buchnera, and Profftella which are all endosymbionts, and Escherichia coli, we determined that regions upstream of these asRNAs encode unique conserved patterns of AT/GC richness, GC skew, and sequence motifs which may be involved in asRNA regulation.

Wing plasticity and associated gene expression varies across the pea aphid biotype complex

Developmental phenotypic plasticity is a widespread phenomenon that allows organisms to produce different adult phenotypes in response to different environments. Investigating the molecular mechanisms underlying plasticity has the potential to reveal the precise changes that lead to the evolution of plasticity as a phenotype. Here, we study wing plasticity in multiple host-plant adapted populations of pea aphids as a model for understanding adaptation to different environments within a single species. We describe the wing plasticity response of different "biotypes" to a crowded environment and find differences within as well as among biotypes. We then use transcriptome profiling to compare a highly plastic pea aphid genotype to one that shows no plasticity and find that the latter exhibits no gene expression differences between environments. We conclude that the loss of plasticity has been accompanied by a loss of differential gene expression and therefore that genetic assimilation has occurred. Our gene expression results generalize previous studies that have shown a correlation between plasticity in morphology and gene expression.

Comparison of native and non-native predator consumption rates and prey avoidance behavior in North America and Europe

Novel predator-prey interactions can contribute to the invasion success of non-native predators. For example, native prey can fail to recognize and avoid non-native predators due to a lack of co-evolutionary history and cue dissimilarity with native predators. This might result in a competitive advantage for non-native predators. Numerous lady beetle species were globally redistributed as biological control agents against aphids, resulting in novel predator-prey interactions. Here, we investigated the strength of avoidance behavior of the pea aphid (Acyrthosiphon pisum) toward chemical cues of native lady beetles and non-native Asian Harmonia axyridis and European Coccinella septempunctata and Hippodamia variegata in North America, hypothesizing that cues of non-native lady beetles induce weaker avoidance behavior than cues of co-evolved native lady beetles. Additionally, we compared aphid consumption of lady beetles, examining potential predation advantages of non-native lady beetles. Finally, we compared cue avoidance behavior between North American and European pea aphid populations and aphid consumption of native and non-native lady beetles in North America and Europe. In North America, pea aphids avoided chemical cues of all ladybeetle species tested, regardless of their origin. In contrast to pea aphids in North America, European pea aphids did not avoid cues of the non-native H. axyridis. The non-native H. axyridis and C. septempunctata were among the largest and most voracious lady beetle species tested, on both continents. Consequently, in North America non-native lady beetle species might have a competitive advantage on shared food resources due to their relatively large body size, compared to several native American lady beetle species. In Europe, however, non-native H. axyridis might benefit from missing aphid cue avoidance as well as a large body size. The co-evolutionary time gap between the European and North American invasion of H. axyridis likely explains the intercontinental differences in cue avoidance behavior and might indicate evolution in aphids toward non-native predators.

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.

Low Spontaneous Mutation Rate and Pleistocene Radiation of Pea Aphids

Accurate estimates of divergence times are essential to understand the evolutionary history of species. It allows linking evolutionary histories of the diverging lineages with past geological, climatic, and other changes in environment and shed light on the processes involved in speciation. The pea aphid radiation includes multiple host races adapted to different legume host plants. It is thought that diversification in this system occurred very recently, over the past 8,000–16,000 years. This young age estimate was used to link diversification in pea aphids to the onset of human agriculture, and led to the establishment of the pea aphid radiation as a model system in the study of speciation with gene flow. Here, we re-examine the age of the pea aphid radiation, by combining a mutation accumulation experiment with a genome-wide estimate of divergence between distantly related pea aphid host races. We estimate the spontaneous mutation rate for pea aphids as 2.7×10-10 per haploid genome per parthenogenic generation. Using this estimate of mutation rate and the genome-wide genetic differentiation observed between pea aphid host races, we show that the pea aphid radiation is much more ancient than assumed previously, predating Neolithic agriculture by several hundreds of thousands of years. Our results rule out human agriculture as the driver of diversification of the pea aphid radiation, and call for re-assessment of the role of allopatric isolation during Pleistocene climatic oscillations in divergence of the pea aphid complex.

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.

Genome Evolution of the Obligate Endosymbiont Buchnera aphidicola

An evolutionary consequence of uniparentally transmitted symbiosis is degradation of symbiont genomes. We use the system of aphids and their maternally inherited obligate endosymbiont, Buchnera aphidicola, to explore the evolutionary process of genome degradation. We compared complete genome sequences for 39 Buchnera strains, including 23 newly sequenced symbiont genomes from diverse aphid hosts. We reconstructed the genome of the most recent shared Buchnera ancestor, which contained 616 protein-coding genes, and 39 RNA genes. The extent of subsequent gene loss varied across lineages, resulting in modern genomes ranging from 412 to 646 kb and containing 354-587 protein-coding genes. Loss events were highly nonrandom across loci. Genes involved in replication, transcription, translation, and amino acid biosynthesis are largely retained, whereas genes underlying ornithine biosynthesis, stress responses, and transcriptional regulation were lost repeatedly. Aside from losses, gene order is almost completely stable. The main exceptions involve movement between plasmid and chromosome locations of genes underlying tryptophan and leucine biosynthesis and supporting nutrition of aphid hosts. This set of complete genomes enabled tests for signatures of positive diversifying selection. Of 371 Buchnera genes tested, 29 genes show strong support for ongoing positive selection. These include genes encoding outer membrane porins that are expected to be involved in direct interactions with hosts. Collectively, these results indicate that extensive genome reduction occurred in the ancestral Buchnera prior to aphid diversification and that reduction has continued since, with losses greater in some lineages and for some loci.

Bioactivity-guided isolation of rosmarinic acid as the principle bioactive compound from the butanol extract of Isodon rugosus against the pea aphid, Acyrthosiphon pisum

Aphids are agricultural pest insects that transmit viruses and cause feeding damage on a global scale. Current pest control practices involving the excessive use of synthetic insecticides over many years have resulted in aphid resistance to a number of pesticides. In nature, plants produce secondary metabolites during their interaction with insects and these metabolites can act as toxicants, antifeedants, anti-oviposition agents and deterrents towards the insects. In a previous study, we demonstrated that the butanol fraction from a crude methanolic extract of an important plant species, Isodon rugosus showed strong insecticidal activity against the pea aphid, Acyrthosiphon pisum. To further explore this finding, the current study aimed to exploit a bioactivity-guided strategy to isolate and identify the active compound in the butanol fraction of I. rugosus. As such, reversed-phase flash chromatography, acidic extraction and different spectroscopic techniques were used to isolate and identify the new compound, rosmarinic acid, as the bioactive compound in I. rugosus. Insecticidal potential of rosmarinic acid against A. pisum was evaluated using standard protocols and the data obtained was analyzed using qualitative and quantitative statistical approaches. Considering that a very low concentration of this compound (LC90 = 5.4 ppm) causes significant mortality in A. pisum within 24 h, rosmarinic acid could be exploited as a potent insecticide against this important pest insect. Furthermore, I. rugosus is already used for medicinal purposes and rosmarinic acid is known to reduce genotoxic effects induced by chemicals, hence it is expected to be safer compared to the current conventional pesticides. While this study highlights the potential of I. rugosus as a possible biopesticide source against A. pisum, it also provides the basis for further exploration and development of formulations for effective field application.

Resolving the taxonomic status of biocontrol agents belonging to the Aphidius eadyi species group (Hymenoptera: Braconidae: Aphidiinae): an integrative approach

Species that belong to the Aphidius eadyi group have been used as biocontrol agents against Acyrthosiphon pisum worldwide. However, despite their extensive use, there are still gaps in our knowledge about their taxonomy and distribution. In this study, we employed an integrative taxonomic approach by combining genetic analyses (mtDNA COI barcoding) with standard morphological analyses and geometric morphometrics of forewing shape. We identified three species within the A. eadyi species group, viz., A. smithi, A. eadyi and A. banksae. Genetic separation of all three species was confirmed, with mean genetic distances between species ranging from 5 to 7.4%. The following morphological characters were determined as the most important for separating species of the A. eadyi group: number and shape of costulae on the anterolateral part of the petiole, shape of the central areola on the propodeum, and shape and venation of the forewings. The differences in wing shape of all three species were statistically significant, but with some overlapping. We identified A. banksae as a widely distributed pea aphid parasitoid, whose known range covers most of the western Palaearctic (from the UK to Israel). Aphidius banksae is diagnosed and redescribed.

Untargeted Metabolomics Approach Reveals Differences in Host Plant Chemistry Before and After Infestation With Different Pea Aphid Host Races

The pea aphid (Acyrthosiphon pisum), a phloem-sucking insect, has undergone a rapid radiation together with the domestication and anthropogenic range expansion of several of its legume host plants. This insect species is a complex of at least 15 genetically different host races that can all develop on the universal host plant Vicia faba. However, each host race is specialized on a particular plant species, such as Medicago sativa, Trifolium pratense, or Pisum sativum, which makes it an attractive model insect to study ecological speciation. Previous work revealed that pea aphid host plants produce a specific phytohormone profile depending on the host plant - host race combination. Native aphid races induce lower defense hormone levels in their host plant than non-native pea aphid races. Whether these changes in hormone levels also lead to changes in other metabolites is still unknown. We used a mass spectrometry-based untargeted metabolomic approach to identify plant chemical compounds that vary among different host plant-host race combinations and might therefore, be involved in pea aphid host race specialization. We found significant differences among the metabolic fingerprints of the four legume species studied prior to aphid infestation, which correlated with aphid performance. After infestation, the metabolic profiles of M. sativa and T. pratense plants infested with their respective native aphid host race were consistently different from profiles after infestation with non-native host races and from uninfested control plants. The metabolic profiles of P. sativum plants infested with their native aphid host race were also different from plants infested with non-native host races, but not different from uninfested control plants. The compounds responsible for these differences were putatively identified as flavonoids, saponins, non-proteinogenic amino acids and peptides among others. As members of these compound classes are known for their activity against insects and aphids in particular, they may be responsible for the differential performance of host races on native vs. non-native host plants. We conclude that the untargeted metabolomic approach is suitable to identify candidate compounds involved in the specificity of pea aphid - host plant interactions.

Fast Evolution and Lineage-Specific Gene Family Expansions of Aphid Salivary Effectors Driven by Interactions with Host-Plants

Effector proteins play crucial roles in plant-parasite interactions by suppressing plant defenses and hijacking plant physiological responses to facilitate parasite invasion and propagation. Although effector proteins have been characterized in many microbial plant pathogens, their nature and role in adaptation to host plants are largely unknown in insect herbivores. Aphids rely on salivary effector proteins injected into the host plants to promote phloem sap uptake. Therefore, gaining insight into the repertoire and evolution of aphid effectors is key to unveiling the mechanisms responsible for aphid virulence and host plant specialization. With this aim in mind, we assembled catalogues of putative effectors in the legume specialist aphid, Acyrthosiphon pisum, using transcriptomics and proteomics approaches. We identified 3,603 candidate effector genes predicted to be expressed in A. pisum salivary glands (SGs), and 740 of which displayed up-regulated expression in SGs in comparison to the alimentary tract. A search for orthologs in 17 arthropod genomes revealed that SG-up-regulated effector candidates of A. pisum are enriched in aphid-specific genes and tend to evolve faster compared with the whole gene set. We also found that a large fraction of proteins detected in the A. pisum saliva belonged to three gene families, of which certain members show evidence consistent with positive selection. Overall, this comprehensive analysis suggests that the large repertoire of effector candidates in A. pisum constitutes a source of novelties promoting plant adaptation to legumes.

Disentangling the Causes for Faster-X Evolution in Aphids

The faster evolution of X chromosomes has been documented in several species, and results from the increased efficiency of selection on recessive alleles in hemizygous males and/or from increased drift due to the smaller effective population size of X chromosomes. Aphids are excellent models for evaluating the importance of selection in faster-X evolution because their peculiar life cycle and unusual inheritance of sex chromosomes should generally lead to equivalent effective population sizes for X and autosomes. Because we lack a high-density genetic map for the pea aphid, whose complete genome has been sequenced, we first assigned its entire genome to the X or autosomes based on ratios of sequencing depth in males (X0) to females (XX). Then, we computed nonsynonymous to synonymous substitutions ratios (dN/dS) for the pea aphid gene set and found faster evolution of X-linked genes. Our analyses of substitution rates, together with polymorphism and expression data, showed that relaxed selection is likely to be the greatest contributor to faster-X because a large fraction of X-linked genes are expressed at low rates and thus escape selection. Yet, a minor role for positive selection is also suggested by the difference between substitution rates for X and autosomes for male-biased genes (but not for asexual female-biased genes) and by lower Tajima’s D for X-linked compared with autosomal genes with highly male-biased expression patterns. This study highlights the relevance of organisms displaying alternative chromosomal inheritance to the understanding of forces shaping genome evolution.

Multi-scale characterization of symbiont diversity in the pea aphid complex through metagenomic approaches

BACKGROUND

Most metazoans are involved in durable relationships with microbes which can take several forms, from mutualism to parasitism. The advances of NGS technologies and bioinformatics tools have opened opportunities to shed light on the diversity of microbial communities and to give some insights into the functions they perform in a broad array of hosts. The pea aphid is a model system for the study of insect-bacteria symbiosis. It is organized in a complex of biotypes, each adapted to specific host plants. It harbors both an obligatory symbiont supplying key nutrients and several facultative symbionts bringing additional functions to the host, such as protection against biotic and abiotic stresses. However, little is known on how the symbiont genomic diversity is structured at different scales: across host biotypes, among individuals of the same biotype, or within individual aphids, which limits our understanding on how these multi-partner symbioses evolve and interact.

RESULTS

We present a framework well adapted to the study of genomic diversity and evolutionary dynamics of the pea aphid holobiont from metagenomic read sets, based on mapping to reference genomes and whole genome variant calling. Our results revealed that the pea aphid microbiota is dominated by a few heritable bacterial symbionts reported in earlier works, with no discovery of new microbial associates. However, we detected a large and heterogeneous genotypic diversity associated with the different symbionts of the pea aphid. Partitioning analysis showed that this fine resolution diversity is distributed across the three considered scales. Phylogenetic analyses highlighted frequent horizontal transfers of facultative symbionts between host lineages, indicative of flexible associations between the pea aphid and its microbiota. However, the evolutionary dynamics of symbiotic associations strongly varied depending on the symbiont, reflecting different histories and possible constraints. In addition, at the intra-host scale, we showed that different symbiont strains may coexist inside the same aphid host.

CONCLUSIONS

We present a methodological framework for the detailed analysis of NGS data from microbial communities of moderate complexity and gave major insights into the extent of diversity in pea aphid-symbiont associations and the range of evolutionary trajectories they could take.

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long-term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant-specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model-based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population-specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host-plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host-plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.

Rapid evolution of aphid pests in agricultural environments

Aphids constitute a major group of crop pests that inflict serious damages to plants, both directly by ingesting phloem and indirectly as vectors of numerous diseases. In response to intense and repeated human-induced pressures, such as insecticide treatments, the use of resistant plants and biological agents, aphids have developed a series of evolutionary responses relying on adaptation and phenotypic plasticity. In this review, we highlight some remarkable evolutionary responses to anthropogenic pressures in agroecosystems and discuss the mechanisms underlying the ecological and evolutionary success of aphids. We outline the peculiar mode of reproduction, the polyphenism for biologically important traits and the diverse and flexible associations with microbial symbionts as key determinants of adaptive potential and pest status of aphids.

Symbiont interactions with non-native hosts limit the formation of new symbioses

Background

Facultative symbionts are common in eukaryotes and can provide their hosts with significant fitness benefits. Despite the advantage of carrying these microbes, they are typically only found in a fraction of the individuals within a population and are often non-randomly distributed among host populations. It is currently unclear why facultative symbionts are only found in certain host individuals and populations. Here we provide evidence for a mechanism to help explain this phenomenon: that when symbionts interact with non-native host genotypes it can limit the horizontal transfer of symbionts to particular host lineages and populations of related hosts.

Results

Using reciprocal transfections of the facultative symbiont Hamiltonella defensa into different pea aphid clones, we demonstrate that particular symbiont strains can cause high host mortality and inhibit offspring production when injected into aphid clones other than their native host lineage. However, once established, the symbiont’s ability to protect against parasitoids was not influenced by its origin. We then demonstrate that H. defensa is also more likely to establish a symbiotic relationship with aphid clones from a plant-adapted population (biotype) that typically carry H. defensa in nature, compared to clones from a biotype that does not normally carry this symbiont.

Conclusions

These results provide evidence that certain aphid lineages and populations of related hosts are predisposed to establishing a symbiotic relationship with H. defensa. Our results demonstrate that host-symbiont genotype interactions represent a potential barrier to horizontal transmission that can limit the spread of symbionts, and adaptive traits they carry, to certain host lineages.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1143-z) contains supplementary material, which is available to authorized users.

Establishment and maintenance of aphid endosymbionts after horizontal transfer is dependent on host genotype

Animal-associated microbial communities have important effects on host phenotypes. Individuals within and among species differ in the strains and species of microbes that they harbour, but how natural selection shapes the distribution and abundance of symbionts in natural populations is not well understood. Symbionts can be beneficial in certain environments but also impose costs on their hosts. Consequently, individuals that can or cannot associate with symbionts will be favoured under different ecological circumstances. As a result, we predict that individuals within a species vary in terms of how well they accept and maintain symbionts. In pea aphids, the frequency of endosymbionts varies among host-plant-associated populations ('biotypes'). We show that aphid genotypes from different biotypes vary in how well they accept and maintain symbionts after horizontal transfer. We find that aphids from biotypes that frequently harbour symbionts are better able to associate with novel symbionts than those from biotypes that less frequently harbour symbionts. Intraspecific variation in the ability of hosts to interact with symbionts is an understudied factor explaining patterns of host-symbiont association.

Genetic Structure of the Bacterial Endosymbiont Buchnera aphidicola from Its Host Aphid Schlechtendalia chinensis and Evolutionary Implications

Buchnera aphidicola is a primary symbiotic bacterium which provides essential amino acids to aphids. In this study, we sequenced nuclear 16s rDNA and atpAGD genes for 156 individuals of B. aphidicola from eight geographically distant populations to investigate the genetic diversity and structure of B. aphidicola associated to the sumac gall aphid Schlechtendalia chinensis in central and southern China. Our analyses of the combined sequences showed that B. aphidicola from S. chinensis had high haplotype and nucleotide diversity (h = 0.893; π = 0.00164). One of the 16 haplotypes detected had a wide geographic distribution across the central and southern China and was probably the ancestral haplotype of B. aphidicola from S. chinensis. A network and phylogenetic analysis revealed a geographic structure in which the 16 haplotypes of B. aphidicola were divided into the northern and southern clades separated by the Yangtze River. The two clades diverged from each other at 22.1 ± 3.7 Mya according to our divergence time estimations. Therefore, the modern genetic structure in B. aphidicola from S. chinensis has been probably impacted by historical geological events. Combined with the data from GenBank, we also reconstructed the phylogenetic relationships of three aphid subfamilies and their symbiont bacteria. The results indicated significant topological correlations between the aphid and bacterial phylogenies at interspecific levels.

Insect-bacteria parallel evolution in multiple-co-obligate-aphid association: a case in Lachninae (Hemiptera: Aphididae)

Parallel phylogenies between aphid and its obligate symbiont Buchnera are hot topics which always focused on aphid lower taxonomic levels. Symbionts in the subfamily Lachninae are special. Buchnera in many lachnine species has undergone functional and genome size reduction that was replaced by other co-obligate symbionts. In this study, we constructed the phylogenetic relationships of Lachninae with a combined dataset of five genes sequenced from Buchnera to estimate the effects of a dual symbiotic system in the aphid-Buchnera cospeciation association. The phylogeny of Buchnera in Lachninae was well-resolved in the combined dataset. Each of the genera formed strongly supported monophyletic groups, with the exception of the genus Cinara. The phylogeny based on sequences from Buchnera was divided into five tribes according to the clades of the Lachninae hosts tree, with the phylogenies of Buchnera and Lachninae being generally congruent. These results first provided evidence of parallel evolution at the aphid subfamily level comprehensively and supported the view that topological congruence between the phylogenies of Buchnera and Lachninae would not be interfered with the other co-obligate symbionts, such as Sarretia, in aphid-entosymbiont association. These results also provided new insight in understanding host-plant coevolution in lachnine lineages.

Aphid specialism as an example of ecological-evolutionary divergence

Debate still continues around the definition of generalism and specialism in nature. To some, generalism is equated solely with polyphagy, but this cannot be readily divorced from other essential biological factors, such as morphology, behaviour, genetics, biochemistry, chemistry and ecology, including chemical ecology. Viewed in this light, and accepting that when living organisms evolve to fill new ecological-evolutionary niches, this is the primal act of specialisation, then perhaps all living organisms are specialist in the broadest sense. To illustrate the levels of specialisation that may be found in a group of animals, we here provide an overview of those displayed by a subfamily of hemipteran insects, the Aphididae, which comprises some 1600 species/subspecies in Europe alone and whose members are specialised in a variety of lifestyle traits. These include life cycle, host adaptation, dispersal and migration, associations with bacterial symbionts (in turn related to host adaptation and resistance to hymenopterous wasp parasitoids), mutualisms with ants, and resistance to insecticides. As with polyphagy, these traits cannot easily be separated from one another, but rather, are interconnected, often highly so, which makes the Aphididae a fascinating animal group to study, providing an informative, perhaps unique, model to illustrate the complexities of defining generalism versus specialism.

Genomic adaptation to agricultural environments: cabbage white butterflies (Pieris rapae) as a case study

Background

Agricultural environments have long presented an opportunity to study evolution in action, and genomic approaches are opening doors for testing hypotheses about adaptation to crops, pesticides, and fertilizers. Here, we begin to develop the cabbage white butterfly (Pieris rapae) as a system to test questions about adaptation to novel, agricultural environments. We focus on a population in the north central United States as a unique case study: here, canola, a host plant, has been grown during the entire flight period of the butterfly over the last three decades.

Results

First, we show that the agricultural population has diverged phenotypically relative to a nonagricultural population: when reared on a host plant distantly related to canola, the agricultural population is smaller and more likely to go into diapause than the nonagricultural population. Second, drawing from deep sequencing runs from six individuals from the agricultural population, we assembled the gut transcriptome of this population. Then, we sequenced RNA transcripts from the midguts of 96 individuals from this canola agricultural population and the nonagricultural population in order to describe patterns of genomic divergence between the two. While population divergence is low, 235 genes show evidence of significant differentiation between populations. These genes are significantly enriched for cofactor and small molecule metabolic processes, and many genes also have transporter or catalytic activity. Analyses of population structure suggest the agricultural population contains a subset of the genetic variation in the nonagricultural population.

Conclusions

Taken together, our results suggest that adaptation of cabbage whites to an agricultural environment occurred at least in part through selection on standing genetic variation. Both the phenotypic and genetic data are consistent with the idea that this pest has adapted to an abundant and predictable agricultural resource through a narrowing of niche breadth and loss of genetic variants rather than de novo gain of adaptive alleles. The present research develops genomic resources to pave the way for future studies using cabbage whites as a model contributing to our understanding of adaptation to agricultural environments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3787-2) contains supplementary material, which is available to authorized users.

Modulation of Legume Defense Signaling Pathways by Native and Non-native Pea Aphid Clones

The pea aphid (Acyrthosiphon pisum) is a complex of at least 15 genetically different host races that are native to specific legume plants, but can all develop on the universal host plant Vicia faba. Despite much research, it is still unclear why pea aphid host races (biotypes) are able to colonize their native hosts while other host races are not. All aphids penetrate the plant and salivate into plant cells when they test plant suitability. Thus plants might react differently to the various pea aphid host races. To find out whether legume species vary in their defense responses to different pea aphid host races, we measured the amounts of salicylic acid (SA), the jasmonic acid-isoleucine conjugate (JA-Ile), other jasmonate precursors and derivatives, and abscisic acid (ABA) in four different species (Medicago sativa, Trifolium pratense, Pisum sativum, V. faba) after infestation by native and non-native pea aphid clones of various host races. Additionally, we assessed the performance of the clones on the four plant species. On M. sativa and T. pratense, non-native clones that were barely able to survive or reproduce, triggered a strong SA and JA-Ile response, whereas infestation with native clones led to lower levels of both phytohormones. On P. sativum, non-native clones, which survived or reproduced to a certain extent, induced fluctuating SA and JA-Ile levels, whereas the native clone triggered only a weak SA and JA-Ile response. On the universal host V. faba all aphid clones triggered only low SA levels initially, but induced clone-specific patterns of SA and JA-Ile later on. The levels of the active JA-Ile conjugate and of the other JA-pathway metabolites measured showed in many cases similar patterns, suggesting that the reduction in JA signaling was due to an effect upstream of OPDA. ABA levels were downregulated in all aphid clone-plant combinations and were therefore probably not decisive factors for aphid-plant compatibility. Our results suggest that A. pisum clones manipulate plant-defense signaling to their own advantage, and perform better on their native hosts due to their ability to modulate the SA- and JA-defense signaling pathways.

Population genetic structure of cotton pink bollworm, Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae) using mitochondrial cytochrome oxidase I (COI) gene sequences from India

Pink bollworm (PBW), Pectinophora gossypiella is one of the most destructive pest's globally inflicting huge economic losses in cotton even during later stages of crop growth. In the present investigation, the population genetic structure, distribution, and genetic diversity of P. gossypiella in cotton growing zones of India using partial mitochondrial DNA cytochrome oxidase-I (COI) gene was addressed. The overall haplotype (Hd), number of nucleotide differences (K), and nucleotide diversity (π) were 0.3028, 0.327, and 0.00047, respectively which suggest that entire population exhibited low level of genetic diversity. Zone-wise clustering of population revealed that central zone recorded low level of Hd (0.2730) as compared to north (0.3619) and south (0.3028) zones. The most common haplotype (H1) reported in all 19 locations could be proposed as ancestral/original haplotype. This haplotype with one mutational step formed star-like phylogeny connected with 11 other haplotypes. The phylogenetic relationship studies revealed that most haplotypes of populations are closely related to each other. Haplotype 5 was exclusively present in Dharwad (South zone) shared with populations of Hanumangarh and Bathinda (North zone). The result indicated that there is no isolation by distance effect among the Indian populations of PBW. The present study reports a low genetic diversity among PBW populations of India and H1, as ancestral haplotype from which other haplotypes have evolved suggests that the migration and dispersal over long distance and invasiveness are major factors.

Differential gene expression according to race and host plant in the pea aphid

Host-race formation in phytophagous insects is thought to provide the opportunity for local adaptation and subsequent ecological speciation. Studying gene expression differences amongst host races may help to identify phenotypes under (or resulting from) divergent selection and their genetic, molecular and physiological bases. The pea aphid (Acyrthosiphon pisum) comprises host races specializing on numerous plants in the Fabaceae and provides a unique system for examining the early stages of diversification along a gradient of genetic and associated adaptive divergence. In this study, we examine transcriptome-wide gene expression both in response to environment and across pea aphid races selected to cover the range of genetic divergence reported in this species complex. We identify changes in expression in response to host plant, indicating the importance of gene expression in aphid-plant interactions. Races can be distinguished on the basis of gene expression, and higher numbers of differentially expressed genes are apparent between more divergent races; these expression differences between host races may result from genetic drift and reproductive isolation and possibly divergent selection. Expression differences related to plant adaptation include a subset of chemosensory and salivary genes. Genes showing expression changes in response to host plant do not make up a large portion of between-race expression differences, providing confirmation of previous studies' findings that genes involved in expression differences between diverging populations or species are not necessarily those showing initial plasticity in the face of environmental change.

A genomic perspective on the generation and maintenance of genetic diversity in herbivorous insects

Understanding the processes that generate and maintain genetic variation within populations is a central goal in evolutionary biology. Theory predicts that some of this variation is maintained as a consequence of adapting to variable habitats. Studies in herbivorous insects have played a key role in confirming this prediction. Here, we highlight theoretical and conceptual models for the maintenance of genetic diversity in herbivorous insects, empirical genomic studies testing these models, and pressing questions within the realm of evolutionary and functional genomic studies. To address key gaps, we propose an integrative approach combining population genomic scans for adaptation, genome-wide characterization of targets of selection through experimental manipulations, mapping the genetic architecture of traits influencing fitness, and functional studies. We also stress the importance of studying the maintenance of genetic variation across biological scales-from variation within populations to divergence among populations-to form a comprehensive view of adaptation in herbivorous insects.

Performances of survival, feeding behavior, and gene expression in aphids reveal their different fitness to host alteration

Insect populations feeding on different plant species are under selection pressure to adapt to these differences. A study integrating elements of the ecology, behavior, and gene expression of aphids on different host plants has not yet been well-explored. The present study explores the relationship between host fitness and survival, feeding behavior, and salivary gland gene expression of a pea (Pisum sativum) host race of Acyrthosiphon pisum feeding on a common host Vicia faba and on three genetically-related hosts (Vicia villosa, Medicago truncatula, and Medicago sativa). Life table data indicated that aphids on non-favored hosts exhibited small size, low reproduction rate, slow population increase and individual development, and long lifespan. Electrical penetration graph results showed that the aphids spent significantly less time in passive ingestion of phloem sap on all non-preferred host plants before acclimation. After a period of acclimation on M. truncatula and V. villosa, pea host race individuals showed improved feeding behavior. No individuals of the pea host race completed its life history on M. sativa. Interestingly, the number of host-specific differentially-expressed salivary gland genes was negatively correlated with the fitness of aphids on this host plant. This study provided important cues in host plant specialization in aphids.

Enemy-free space promotes maintenance of host races in an aphid species

The enormous biodiversity of herbivorous insects may arise from ecological speciation via continuous host-plant switches. Whether such switches are successful depends on the trade-off between different selection pressures that act on herbivores. Decreased herbivore performance due to suboptimal nutrition might be compensated for by a reduced natural enemy pressure. As a consequence, an "enemy-free space" on a certain plant might facilitate host-plant switches and maintain biotypes. To test this hypothesis, we used the pea aphid (Acyrthosiphon pisum) complex, which consists of at least 11 genetically distinct host races that are native to specific legume host plants but can all develop on the universal host plant Vicia faba. Three A. pisum host races native to Trifolium pratense, Pisum sativum, and Medicago sativa were investigated in experiments on their respective host plants and on the universal host plant V. faba. We found that hoverflies preferred to oviposit on P. sativum and the universal host V. faba. Since feeding by hoverfly larvae suppressed aphid population growth on these host plants, the native hosts M. sativa and T. pratense provided enemy-free space for the respective A. pisum races. Mobile predators, such as ants and ladybird beetles, preferred Pisum race aphids on V. faba over P. sativum. Thus, all three of the native host plants studied supply enemy-free space for A. pisum compared to the universal host V. faba. Reducing encounters between aphid races on V. faba would reduce gene flow among them and could contribute to maintaining the host races.

Selection on a genetic polymorphism counteracts ecological speciation in a stick insect

The interplay between selection and aspects of the genetic architecture of traits (such as linkage, dominance, and epistasis) can either drive or constrain speciation [1-3]. Despite accumulating evidence that speciation can progress to "intermediate" stages-with populations evolving only partial reproductive isolation-studies describing selective mechanisms that impose constraints on speciation are more rare than those describing drivers. The stick insect Timema cristinae provides an example of a system in which partial reproductive isolation has evolved between populations adapted to different host plant environments, in part due to divergent selection acting on a pattern polymorphism [4, 5]. Here, we demonstrate how selection on a green/melanistic color polymorphism counteracts speciation in this system. Specifically, divergent selection between hosts does not occur on color phenotypes because melanistic T. cristinae are cryptic on the stems of both host species, are resistant to a fungal pathogen, and have a mating advantage. Using genetic crosses and genome-wide association mapping, we quantify the genetic architecture of both the pattern and color polymorphism, illustrating their simple genetic control. We use these empirical results to develop an individual-based model that shows how the melanistic phenotype acts as a "genetic bridge" that increases gene flow between populations living on different hosts. Our results demonstrate how variation in the nature of selection acting on traits, and aspects of trait genetic architecture, can impose constraints on both local adaptation and speciation.

Engineering plants for aphid resistance: current status and future perspectives

The current status of development of transgenic plants for improved aphid resistance, and the pros and cons of different strategies are reviewed and future perspectives are proposed. Aphids are major agricultural pests that cause significant yield losses of crop plants each year. Excessive dependence on insecticides for aphid control is undesirable because of the development of insecticide resistance, the potential negative effects on non-target organisms and environmental pollution. Transgenic plants engineered for resistance to aphids via a non-toxic mode of action could be an efficient alternative strategy. In this review, the distribution of major aphid species and their damages on crop plants, the so far isolated aphid-resistance genes and their applications in developments of transgenic plants for improved aphid resistance, and the pros and cons of these strategies are reviewed and future perspectives are proposed. Although the transgenic plants developed through expressing aphid-resistant genes, manipulating plant secondary metabolism and plant-mediated RNAi strategy have been demonstrated to confer improved aphid resistance to some degree. So far, no aphid-resistant transgenic crop plants have ever been commercialized. This commentary is intended to be a helpful insight into the generation and future commercialization of aphid-resistant transgenic crops in a global context.

Widespread host-dependent hybrid unfitness in the pea aphid species complex

Linking adaptive divergence to hybrid unfitness is necessary to understand the ecological factors contributing to reproductive isolation and speciation. To date, this link has been demonstrated in few model systems, most of which encompass ecotypes that occupy relatively early stages in the speciation process. Here we extend these studies by assessing how host-plant adaptation conditions hybrid fitness in the pea aphid, Acyrthosiphon pisum. We made crosses between and within five pea aphid biotypes adapted to different host plants and representing various stages of divergence within the complex. Performance of F1 hybrids and nonhybrids was assessed on a "universal" host that is favorable to all pea aphid biotypes in laboratory conditions. Although hybrids performed equally well as nonhybrids on the universal host, their performance was much lower than nonhybrids on the natural hosts of their parental populations. Hence, hybrids, rather than being intrinsically deficient, are maladapted to their parents' hosts. Interestingly, the impact of this maladaptation was stronger in certain hybrids from crosses involving the most divergent biotype, suggesting that host-dependent postzygotic isolation has continued to evolve late in divergence. Even though host-independent deficiencies are not excluded, hybrid maladaptation to parental hosts supports the hypothesis of ecological speciation in this complex.

Whole-genome re-sequencing of non-model organisms: lessons from unmapped reads

Unmapped reads are often discarded from the analysis of whole-genome re-sequencing, but new biological information and insights can be uncovered through their analysis. In this paper, we investigate unmapped reads from the re-sequencing data of 33 pea aphid genomes from individuals specialized on different host plants. The unmapped reads for each individual were retrieved following mapping to the Acyrthosiphon pisum reference genome and its mitochondrial and symbiont genomes. These sets of unmapped reads were then cross-compared, revealing that a significant number of these unmapped sequences were conserved across individuals. Interestingly, sequences were most commonly shared between individuals adapted to the same host plant, suggesting that these sequences may contribute to the divergence between host plant specialized biotypes. Analysis of the contigs obtained from assembling the unmapped reads pooled by biotype allowed us to recover some divergent genomic regions previously excluded from analysis and to discover putative novel sequences of A. pisum and its symbionts. In conclusion, this study emphasizes the interest of the unmapped component of re-sequencing data sets and the potential loss of important information. We here propose strategies to aid the capture and interpretation of this information.

Dynamics of copy number variation in host races of the pea aphid

Copy number variation (CNV) makes a major contribution to overall genetic variation and is suspected to play an important role in adaptation. However, aside from a few model species, the extent of CNV in natural populations has seldom been investigated. Here, we report on CNV in the pea aphid Acyrthosiphon pisum, a powerful system for studying the genetic architecture of host-plant adaptation and speciation thanks to multiple host races forming a continuum of genetic divergence. Recent studies have highlighted the potential importance of chemosensory genes, including the gustatory and olfactory receptor gene families (Gr and Or, respectively), in the process of host race formation. We used targeted resequencing to achieve a very high depth of coverage, and thereby revealed the extent of CNV of 434 genes, including 150 chemosensory genes, in 104 individuals distributed across eight host races of the pea aphid. We found that CNV was widespread in our global sample, with a significantly higher occurrence in multigene families, especially in Ors. We also observed a decrease in the gene probability of being completely duplicated or deleted (CDD) with increase in coding sequence length. Genes with CDD variants were usually more polymorphic for copy number, especially in the P450 gene family where toxin resistance may be related to gene dosage. We found that Gr were overrepresented among genes discriminating host races, as were CDD genes and pseudogenes. Our observations shed new light on CNV dynamics and are consistent with CNV playing a role in both local adaptation and speciation.

Genomic regions repeatedly involved in divergence among plant-specialized pea aphid biotypes

Understanding the genetic bases of biological diversification is a long-standing goal in evolutionary biology. Here, we investigate whether replicated cases of adaptive divergence involve the same genomic regions in the pea aphid, Acyrthosiphon pisum, a large complex of genetically differentiated biotypes, each specialized on different species of legumes. A previous study identified genomic regions putatively involved in host-plant adaptation and/or reproductive isolation by performing a hierarchical genome scan in three biotypes. This led to the identification of 11 F(ST) outliers among 390 polymorphic microsatellite markers. In this study, the outlier status of these 11 loci was assessed in eight biotypes specialized on other host plants. Four of the 11 previously identified outliers showed greater genetic differentiation among these additional biotypes than expected under the null hypothesis of neutral evolution (α < 0.01). Whether these hotspots of genomic divergence result from adaptive events, intrinsic barriers or reduced recombination is discussed.

Aphid-encoded variability in susceptibility to a parasitoid

BACKGROUND

Many animals exhibit variation in resistance to specific natural enemies. Such variation may be encoded in their genomes or derived from infection with protective symbionts. The pea aphid, Acyrthosiphon pisum, for example, exhibits tremendous variation in susceptibility to a common natural enemy, the parasitic wasp Aphidius ervi. Pea aphids are often infected with the heritable bacterial symbiont, Hamiltonella defensa, which confers partial to complete resistance against this parasitoid depending on bacterial strain and associated bacteriophages. That previous studies found that pea aphids without H. defensa (or other symbionts) were generally susceptible to parasitism, together with observations of a limited encapsulation response, suggested that pea aphids largely rely on infection with H. defensa for protection against parasitoids. However, the limited number of uninfected clones previously examined, and our recent report of two symbiont-free resistant clones, led us to explicitly examine aphid-encoded variability in resistance to parasitoids.

RESULTS

After rigorous screening for known and unknown symbionts, and microsatellite genotyping to confirm clonal identity, we conducted parasitism assays using fifteen clonal pea aphid lines. We recovered significant variability in aphid-encoded resistance, with variation levels comparable to that contributed by H. defensa. Because resistance can be costly, we also measured aphid longevity and cumulative fecundity of the most and least resistant aphid lines under permissive conditions, but found no trade-offs between higher resistance and these fitness parameters.

CONCLUSIONS

These results indicate that pea aphid resistance to A. ervi is more complex than previously appreciated, and that aphids employ multiple tactics to aid in their defense. While we did not detect a tradeoff, these may become apparent under stressful conditions or when resistant and susceptible aphids are in direct competition. Understanding sources and amounts of variation in resistance to natural enemies is necessary to understand the ecological and evolutionary dynamics of antagonistic interactions, such as the potential for coevolution, but also for the successful management of pest populations through biological control.

Inheritance patterns of secondary symbionts during sexual reproduction of pea aphid biotypes

Herbivorous insects frequently harbor bacterial symbionts that affect their ecology and evolution. Aphids host the obligatory endosymbiont Buchnera, which is required for reproduction, together with facultative symbionts whose frequencies vary across aphid populations. These maternally transmitted secondary symbionts have been particularly studied in the pea aphid, Acyrthosiphon pisum, which harbors at least 8 distinct bacterial species (not counting Buchnera) having environmentally dependent effects on host fitness. In particular, these symbiont species are associated with pea aphid populations feeding on specific plants. Although they are maternally inherited, these bacteria are occasionally transferred across insect lineages. One mechanism of such nonmaternal transfer is paternal transmission to the progeny during sexual reproduction. To date, transmission of secondary symbionts during sexual reproduction of aphids has been investigated in only a handful of aphid lineages and 3 symbiont species. To better characterize this process, we investigated inheritance patterns of 7 symbiont species during sexual reproduction of pea aphids through a crossing experiment involving 49 clones belonging to 9 host-specialized biotypes, and 117 crosses. Symbiont species in the progeny were detected with diagnostic qualitative PCR at the fundatrix stage hatching from eggs and in later parthenogenetic generations. We found no confirmed case of paternal transmission of symbionts to the progeny, and we observed that maternal transmission of a particular symbiont species (Serratia symbiotica) was quite inefficient. We discuss these observations in respect to the ecology of the pea aphid.

Evolutionary relationships of Pemphigus and allied genera (Hemiptera: Aphididae: Eriosomatinae) and their primary endosymbiont, Buchnera aphidicola

Aphids harbor primary endosymbionts, Buchnera aphidicola, in specialized cells within their body cavities. Aphids and Buchnera have strict mutualistic relationships in nutrition exchange. This ancient association has received much attention from researchers who are interested in endosymbiotic evolution. Previous studies have found parallel phylogenetic relationships between non-galling aphids and Buchnera at lower taxonomic levels (genus, species). To understand whether relatively isolated habitats such as galls have effect on the parallel relationships between aphids and Buchnera, the present paper investigated the phylogenetic relationships of gall aphids from Pemphigus and allied genera, which induce pseudo-galls or galls on Populus spp. (poplar) and Buchnera. The molecular phylogenies inferred from three aphid genes (COI, COII and EF-1α) and two Buchnera genes (gnd, 16S rRNA gene) indicated significant congruence between aphids and Buchnera at generic as well as interspecific levels. Interestingly, both aphid and Buchnera phylogenies supported three main clades corresponding to the galling locations of aphids, namely leaf, the joint of leaf blade and petiole, and branch of the host plant. The results suggest phylogenetic conservatism of gall characters, which indicates gall characters are more strongly affected by aphid phylogeny, rather than host plants.

Brood ball-mediated transmission of microbiome members in the dung beetle, Onthophagus taurus (Coleoptera: Scarabaeidae)

Insects feeding on plant sap, blood, and other nutritionally incomplete diets are typically associated with mutualistic bacteria that supplement missing nutrients. Herbivorous mammal dung contains more than 86% cellulose and lacks amino acids essential for insect development and reproduction. Yet one of the most ecologically necessary and evolutionarily successful groups of beetles, the dung beetles (Scarabaeinae) feeds primarily, or exclusively, on dung. These associations suggest that dung beetles may benefit from mutualistic bacteria that provide nutrients missing from dung. The nesting behaviors of the female parent and the feeding behaviors of the larvae suggest that a microbiome could be vertically transmitted from the parental female to her offspring through the brood ball. Using sterile rearing and a combination of molecular and culture-based techniques, we examine transmission of the microbiome in the bull-headed dung beetle, Onthophagus taurus. Beetles were reared on autoclaved dung and the microbiome was characterized across development. A ~1425 bp region of the 16S rRNA identified Pseudomonadaceae, Enterobacteriaceae, and Comamonadaceae as the most common bacterial families across all life stages and populations, including cultured isolates from the 3^rd instar digestive system. Finer level phylotyping analyses based on lepA and gyrB amplicons of cultured isolates placed the isolates closest to Enterobacter cloacae, Providencia stuartii, Pusillimonas sp., Pedobacter heparinus, and Lysinibacillus sphaericus. Scanning electron micrographs of brood balls constructed from sterile dung reveals secretions and microbes only in the chamber the female prepares for the egg. The use of autoclaved dung for rearing, the presence of microbes in the brood ball and offspring, and identical 16S rRNA sequences in both parent and offspring suggests that the O. taurus female parent transmits specific microbiome members to her offspring through the brood chamber. The transmission of the dung beetle microbiome highlights the maintenance and likely importance of this newly-characterized bacterial community.