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

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.

Transcriptome Profiling Revealed Potentially Critical Roles for Digestion and Defense-Related Genes in Insects' Use of Resistant Host Plants: A Case Study with Sitobion Avenae

Using host plant resistance (HPR) in management of insect pests is often environmentally friendly and suitable for sustainable development of agricultural industries. However, this strategy can be limited by rapid evolution of insect populations that overcome HPR, for which the underlying molecular factors and mechanisms are not well understood. To address this issue, we analyzed transcriptomes of two distinct biotypes of the grain aphid, Sitobion avenae (Fabricius), on wheat and barley. This analysis revealed a large number of differentially expressed genes (DEGs) between biotypes 1 and 3 on wheat and barley. The majority of them were common DEGs occurring on both wheat and barley. GO and KEGG enrichment analyses for these common DEGs demonstrated significant expression divergence between both biotypes in genes associated with digestion and defense. Top defense-related common DEGs with the most significant expression changes included three peroxidases, two UGTs (UDP-glycosyltransferase), two cuticle proteins, one glutathione S-transferases (GST), one superoxide dismutase, and one esterase, suggesting their potentially critical roles in the divergence of S. avenae biotypes. A relatively high number of specific DEGs on wheat were identified for peroxidases (9) and P450s (8), indicating that phenolic compounds and hydroxamic acids may play key roles in resistance of wheat against S. avenae. Enrichment of specific DEGs on barley for P450s and ABC transporters suggested their key roles in this aphid's detoxification against secondary metabolites (e.g., alkaloids) in barley. Our results can provide insights into the molecular factors and functions that explain biotype adaptation in insects and their use of resistant plants. This study also has significant implications for developing new resistant cultivars, developing strategies that limit rapid development of insect biotypes, and extending resistant crop cultivars' durability and sustainability in integrated management programs.

The reproductive system of the male and oviparous female of a model organism-the pea aphid, Acyrthosiphon pisum (Hemiptera, Aphididae)

The structure of the reproductive system of the sexual generation-males and oviparous females-of the pea aphid Acyrthosiphon pisum (Harris) (Hemiptera, Aphididae), a serious pest of cultivated plants of Fabaceae, was investigated. For the first time we describe the morphology, histology and ultrastructure of the reproductive system in both morphs of the sexual generation of aphids within one species, using light and fluorescent microscopy, as well as transmission and scanning electron microscopy. The results revealed that males have testes composed of three follicles fused by the upper ends of the vasa efferentia, the vasa deferentia run independently, the accessory glands are asymmetric and the ejaculatory duct shortened. Oviparous females have ovaries composed of seven ovarioles each. The lateral oviducts join to a short common oviduct connected with the unpaired spermatheca and paired accessory glands. Yolky eggs with an aggregation of symbiotic bacteria at the posterior pole are produced. Histologically, the components of genital tracts are broadly similar: the epithelial cells of the walls of the vasa deferentia and accessory glands of the male and oviparous female have secretory functions which correlate with the age of the studied morphs. We also found symbiotic bacteria within the vasa deferentia epithelial cells in males and within the cells of the lateral oviducts of females. Because the pea aphid is listed among the 14 species that are of the greatest economic importance, our results will be useful for managing aphid populations, protecting plants and ensuring global food security.

An integrated protocol for targeted mutagenesis with CRISPR-Cas9 system in the pea aphid

CRISPR-Cas9 technology is a very efficient functional analysis tool and has been developed in several insects to edit their genome through injection of eggs with guide RNAs targeting coding sequences of genes of interest. However, its implementation in aphids is more challenging. Aphids are major pests of crops worldwide that alternate during their life cycle between clonality and sexual reproduction. The production of eggs after mating of sexual individuals is a single yearly event and is necessarily triggered by a photoperiod decrease. Fertilized eggs then experience an obligate 3-month diapause period before hatching as new clonal colonies. Taking into consideration these particularities, we developed in the pea aphid Acyrthosiphon pisum a step-by-step protocol of targeted mutagenesis based on the microinjection within fertilized eggs of CRISPR-Cas9 components designed for the editing of a cuticular protein gene (stylin-01). This protocol includes the following steps: i) the photoperiod-triggered induction of sexual morphs (2 months), ii) the mating and egg collection step (2 weeks), iii) egg microinjection and melanization, iv) the 3-month obligate diapause, v) the hatching of new lineages from injected eggs (2 weeks) and vi) the maintenance of stable lineages (2 weeks). Overall, this 7-month long procedure was applied to three different crosses in order to estimate the impact of the choice of the genetic combination on egg production dynamics by females as well as hatching rates after diapause. Mutation rates within eggs before diapause were estimated at 70-80%. The hatching rate of injected eggs following diapause ranged from 1 to 11% depending on the cross and finally a total of 17 stable lineages were obtained and maintained clonally. Out of these, 6 lineages were mutated at the defined sgRNAs target sites within stylin-01 coding sequence, either at the two alleles (2 lineages) or at one allele (4 lineages). The final germline transmission rate of the mutations was thus around 35%. Our protocol of an efficient targeted mutagenesis opens the avenue for functional studies through genome editing in aphids.

Long-Term Population Studies Uncover the Genome Structure and Genetic Basis of Xenobiotic and Host Plant Adaptation in the Herbivore Tetranychus urticae

Pesticide resistance arises rapidly in arthropod herbivores, as can host plant adaptation, and both are significant problems in agriculture. These traits have been challenging to study as both are often polygenic and many arthropods are genetically intractable. Here, we examined the genetic architecture of pesticide resistance and host plant adaptation in the two-spotted spider mite, Tetranychus urticae, a global agricultural pest. We show that the short generation time and high fecundity of T. urticae can be readily exploited in experimental evolution designs for high-resolution mapping of quantitative traits. As revealed by selection with spirodiclofen, an acetyl-CoA carboxylase inhibitor, in populations from a cross between a spirodiclofen-resistant and a spirodiclofen-susceptible strain, and which also differed in performance on tomato, we found that a limited number of loci could explain quantitative resistance to this compound. These were resolved to narrow genomic intervals, suggesting specific candidate genes, including acetyl-CoA carboxylase itself, clustered and copy variable cytochrome P450 genes, and NADPH cytochrome P450 reductase, which encodes a redox partner for cytochrome P450s. For performance on tomato, candidate genomic regions for response to selection were distinct from those responding to the synthetic compound and were consistent with a more polygenic architecture. In accomplishing this work, we exploited the continuous nature of allele frequency changes across experimental populations to resolve the existing fragmented T. urticae draft genome to pseudochromosomes. This improved assembly was indispensable for our analyses, as it will be for future research with this model herbivore that is exceptionally amenable to genetic studies.

Genomic transitions during host race and species formation

Darwin recognized species as discontinuous, yet considered them to be formed by an incremental process of natural selection. Recent theoretical work on 'genome-wide congealing' is bridging this gap between the gradualism of divergent selection and rapid genome-wide divergence, particularly during ecological speciation-with-gene-flow. Host races and species of phytophagous insects, displaying a spectrum of divergence and gene flow among member taxa, provide model systems for testing predicted non-linear transitions from 'genic' divergence at a few uncoupled loci to 'genomic' divergence with genome-wide coupling of selected loci and strong reproductive isolation. Integrating across natural history, genomics, and evolutionary theory, emerging research suggests a tipping point from 'genic' to 'genomic' divergence between host races and species, during both sympatric speciation and secondary contact.

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.

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.

Genomic Scans across Three Eucalypts Suggest that Adaptation to Aridity is a Genome-Wide Phenomenon

Widespread species spanning strong environmental (e.g., climatic) gradients frequently display morphological and physiological adaptations to local conditions. Some adaptations are common to different species that occupy similar environments. However, the genomic architecture underlying such convergent traits may not be the same between species. Using genomic data from previous studies of three widespread eucalypt species that grow along rainfall gradients in southern Australia, our probabilistic approach provides evidence that adaptation to aridity is a genome-wide phenomenon, likely to involve multiple and diverse genes, gene families and regulatory regions that affect a multitude of complex genetic and biochemical processes.

Targeted re-sequencing confirms the importance of chemosensory genes in aphid host race differentiation

Host‐associated races of phytophagous insects provide a model for understanding how adaptation to a new environment can lead to reproductive isolation and speciation, ultimately enabling us to connect barriers to gene flow to adaptive causes of divergence. The pea aphid (Acyrthosiphon pisum) comprises host races specializing on legume species and provides a unique system for examining the early stages of diversification along a gradient of genetic and associated adaptive divergence. As host choice produces assortative mating, understanding the underlying mechanisms of choice will contribute directly to understanding of speciation. As host choice in the pea aphid is likely mediated by smell and taste, we use capture sequencing and SNP genotyping to test for the role of chemosensory genes in the divergence between eight host plant species across the continuum of differentiation and sampled at multiple locations across western Europe. We show high differentiation of chemosensory loci relative to control loci in a broad set of pea aphid races and localities, using a model‐free approach based on principal component analysis. Olfactory and gustatory receptors form the majority of highly differentiated genes and include loci that were already identified as outliers in a previous study focusing on the three most closely related host races. Consistent indications that chemosensory genes may be good candidates for local adaptation and barriers to gene flow in the pea aphid open the way to further investigations aiming to understand their impact on gene flow and to determine their precise functions in response to host plant metabolites.

The genome of Diuraphis noxia, a global aphid pest of small grains

BACKGROUND

The Russian wheat aphid, Diuraphis noxia Kurdjumov, is one of the most important pests of small grains throughout the temperate regions of the world. This phytotoxic aphid causes severe systemic damage symptoms in wheat, barley, and other small grains as a direct result of the salivary proteins it injects into the plant while feeding.

RESULTS

We sequenced and de novo assembled the genome of D. noxia Biotype 2, the strain most virulent to resistance genes in wheat. The assembled genomic scaffolds span 393 MB, equivalent to 93% of its 421 MB genome, and contains 19,097 genes. D. noxia has the most AT-rich insect genome sequenced to date (70.9%), with a bimodal CpG(O/E) distribution and a complete set of methylation related genes. The D. noxia genome displays a widespread, extensive reduction in the number of genes per ortholog group, including defensive, detoxification, chemosensory, and sugar transporter groups in comparison to the Acyrthosiphon pisum genome, including a 65% reduction in chemoreceptor genes. Thirty of 34 known D. noxia salivary genes were found in this assembly. These genes exhibited less homology with those salivary genes commonly expressed in insect saliva, such as glucose dehydrogenase and trehalase, yet greater conservation among genes that are expressed in D. noxia saliva but not detected in the saliva of other insects. Genes involved in insecticide activity and endosymbiont-derived genes were also found, as well as genes involved in virus transmission, although D. noxia is not a viral vector.

CONCLUSIONS

This genome is the second sequenced aphid genome, and the first of a phytotoxic insect. D. noxia's reduced gene content of may reflect the influence of phytotoxic feeding in shaping the D. noxia genome, and in turn in broadening its host range. The presence of methylation-related genes, including cytosine methylation, is consistent with other parthenogenetic and polyphenic insects. The D. noxia genome will provide an important contrast to the A. pisum genome and advance functional and comparative genomics of insects and other organisms.

Bacterial communities associated with host-adapted populations of pea aphids revealed by deep sequencing of 16S ribosomal DNA

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

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.