Genome Sequence of the Pea Aphid Acyrthosiphon pisum

Aphids are important agricultural pests and also biological models for studies of insect-plant interactions, symbiosis, virus vectoring, and the developmental causes of extreme phenotypic plasticity. Here we present the 464 Mb draft genome assembly of the pea aphid Acyrthosiphon pisum. This first published whole genome sequence of a basal hemimetabolous insect provides an outgroup to the multiple published genomes of holometabolous insects. Pea aphids are host-plant specialists, they can reproduce both sexually and asexually, and they have coevolved with an obligate bacterial symbiont. Here we highlight findings from whole genome analysis that may be related to these unusual biological features. These findings include discovery of extensive gene duplication in more than 2000 gene families as well as loss of evolutionarily conserved genes. Gene family expansions relative to other published genomes include genes involved in chromatin modification, miRNA synthesis, and sugar transport. Gene losses include genes central to the IMD immune pathway, selenoprotein utilization, purine salvage, and the entire urea cycle. The pea aphid genome reveals that only a limited number of genes have been acquired from bacteria; thus the reduced gene count of Buchnera does not reflect gene transfer to the host genome. The inventory of metabolic genes in the pea aphid genome suggests that there is extensive metabolite exchange between the aphid and Buchnera, including sharing of amino acid biosynthesis between the aphid and Buchnera. The pea aphid genome provides a foundation for post-genomic studies of fundamental biological questions and applied agricultural problems.     

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.     

The cytogenetic architecture of the aphid genome

In recent years aphids, with their well‐defined polyphenism, have become favoured as model organisms for the study of epigenetic processes. The availability of the pea aphid (Acyrthosiphon pisum) genome sequence has engendered much research aimed at elucidating the mechanisms by which the phenotypic plasticity of aphids is inherited and controlled. Yet so far this research effort has paid little attention to the cytogenetic processes that play a vital part in the organisation, expression and inheritance of the aphid genome. Aphids have holocentric chromosomes, which have very different properties from the chromosomes with localised centromeres that are found in most other organisms. Here we review the diverse forms of aphid chromosome behaviour that occur during sex determination and male and female meiosis, often in response to environmental changes and mediated by endocrine factors. Remarkable differences occur, even between related species, that could have significant effects on the inheritance of all or parts of the genome. In relation to this, we review the particular features of the distribution of heterochromatin, rDNA genes and other repetitive DNA in aphid chromosomes, and discuss the part that these may play in the epigenetic modification of chromatin structure and function.     

Diversity in symbiont consortia in the pea aphid complex is associated with large phenotypic variation in the insect host

Virtually all eukaryotes host microbial symbionts that influence their phenotype in many ways. In a host population, individuals may differ in their symbiotic complement in terms of symbiont species and strains. Hence, the combined expression of symbiont and host genotypes may generate a range of phenotypic diversity on which selection can operate and influence host population ecology and evolution. Here, we used the pea aphid to examine how the infection with various symbiotic complements contributes to phenotypic diversity of this insect species. The pea aphid hosts an obligate symbiont (Buchnera aphidicola) and several secondary symbionts among which is Hamiltonella defensa. This secondary symbiont confers a protection against parasitoids but can also reduce the host’s longevity and fecundity. These phenotypic effects of H. defensa infection have been described for a small fraction of the pea aphid complex which encompasses multiple plant-specialized biotypes. In this study, we examined phenotypic differences in four pea aphid biotypes where H. defensa occurs at high frequency and sometimes associated with other secondary symbionts. For each biotype, we measured the fecundity, lifespan and level of parasitoid protection in several aphid lineages differing in their symbiotic complement. Our results showed little variation in longevity and fecundity among lineages but strong differences in their protection level. These differences in protective levels largely resulted from the strain type of H. defensa and the symbiotic consortium in the host. This study highlights the important role of symbiotic complement in the emergence of phenotypic divergence among host populations of the same species.     

Identification of Critical Conditions for Immunostaining in the Pea Aphid Embryos: Increasing Tissue Permeability and Decreasing Background Staining

The pea aphid Acyrthosiphon pisum, with a sequenced genome and abundant phenotypic plasticity, has become an emerging model for genomic and developmental studies. Like other aphids, A. pisum propagate rapidly via parthenogenetic viviparous reproduction, where the embryos develop within egg chambers in an assembly-line fashion in the ovariole. Previously we have established a robust platform of whole-mount in situ hybridization allowing detection of mRNA expression in the aphid embryos. For analyzing the expression of protein, though, established protocols for immunostaining the ovarioles of asexual viviparous aphids did not produce satisfactory results. Here we report conditions optimized for increasing tissue permeability and decreasing background staining, both of which were problems when applying established approaches. Optimizations include: (1) incubation of proteinase K (1 µg/ml, 10 min), which was found essential for antibody penetration in mid- and late-stage aphid embryos; (2) replacement of normal goat serum/bovine serum albumin with a blocking reagent supplied by a Digoxigenin (DIG)-based buffer set and (3) application of methanol rather hydrogen peroxide (H₂O₂) for bleaching endogenous peroxidase; which significantly reduced the background staining in the aphid tissues. These critical conditions optimized for immunostaining will allow effective detection of gene products in the embryos of A. pisum and other aphids.     

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

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

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.     

The pea aphid, Acyrthosiphon pisum: an emerging genomic model system for ecological, developmental and evolutionary studies

Aphids display an abundance of adaptations that are not easily studied in existing model systems. Here we review the biology of a new genomic model system, the pea aphid, Acyrthosiphon pisum. We then discuss several phenomena that are particularly accessible to study in the pea aphid: the developmental genetic basis of polyphenisms, aphid-bacterial symbioses, the genetics of adaptation and mechanisms of virus transmission. The pea aphid can be maintained in the laboratory and natural populations can be studied in the field. These properties allow controlled experiments to be performed on problems of direct relevance to natural aphid populations. Combined with new genomic approaches, the pea aphid is poised to become an important model system for understanding the molecular and developmental basis of many ecologically and evolutionarily relevant problems.     

Characterisation of immune responses in the pea aphid, Acyrthosiphon pisum

The innate immune system of insects provides effective defence against a range of parasites and pathogens. The pea aphid, Acyrthosiphon pisum, is a novel study system for investigating host-parasite interactions due to its complex associations with both well-characterised bacterial symbionts and a diversity of pathogens and parasites, including several important biological control agents. However, little is known about the cellular and humoral immune responses of aphids. Here we identify three morphologically distinct types of haemocytes in circulation that we name prohemocytes, granulocytes and oenocytoids. Granulocytes avidly phagocytose Gram negative Escherechia coli and Gram positive Micrococcus luteus while oenocytoids exhibit melanotic activity. Prohaemocytes increase in abundance immediately following an immune challenge, irrespective of the source of stimulus. Pea aphids form melanotic capsules around Sephadex beads but do not form cellular capsules. We also did not detect any antimicrobial peptide activity in the haemolymph using zone of inhibition assays. We discuss these results in relation to recent findings from the pea aphid genome annotation project that suggest that aphids have a reduced immune gene repertoire compared to other insects.     

Early progress in aphid genomics and consequences for plant-aphid interactions studies

Aphids occupy a niche comprising two conceptual realms: a micron-scale feeding site beneath the plant surface, in which a syringe-like appendage mediates chemical exchange with a specific plant cell type; and the larger realm of a metazoan with sensory organs, a nervous system, and behavior, all responsive to the condition of the host plant and the broader environment. The biology that connects these realms is not well understood, but new details are emerging with the help of genomic tools. The power of these tools is set to increase substantially now that the first genome of an aphid is being sequenced and annotated. This has been possible because a community of aphid researchers focused their efforts to develop and share genomic resources through an international consortium. This complete genome sequence, along with other resources, should permit major advances in understanding the complex and peculiar biological traits responsible for aphids' evolutionary success and their damaging effects on agriculture. This review highlights early progress in the application of aphid genomics and identifies key issues of plant-aphid interactions likely to benefit as molecular tools are further developed. Use of this new knowledge could make significant contributions to crop protection against these and other phloem-feeding insects.     

Hitching a lift on the road to speciation

Understanding how speciation can take place in the presence of homogenizing gene flow remains a major challenge in evolutionary biology. In the early stages of ecological speciation, reproductive isolation between populations occupying different habitats is expected to be concentrated around genes for local adaptation. These genomic regions will show high divergence while gene exchange in other regions of the genome should continue relatively unimpaired, resulting in low levels of differentiation. The problem is to explain how speciation progresses from this point towards complete reproductive isolation, allowing genome‐wide divergence. A new study by Via and West (2008) on speciation between host races of the pea aphid, Acyrthosiphon pisum, introduces the mechanism of ‘divergence hitchhiking’ which can generate large ‘islands of differentiation’ and facilitate the build‐up of linkage disequilibrium, favouring increased reproductive isolation. This idea potentially removes a major stumbling block to speciation under continuous gene flow.     

Chromosome‐level genome assembly for the horned‐gall aphid provides insights into interactions between gall‐making insect and its host plant

The aphid Schlechtendalia chinensis is an economically important insect that can induce horned galls, which are valuable for the medicinal and chemical industries. Up to now, more than twenty aphid genomes have been reported. Most of the sequenced genomes are derived from free‐living aphids. Here, we generated a high‐quality genome assembly from a galling aphid. The final genome assembly is 271.52 Mb, representing one of the smallest sequenced genomes of aphids. The genome assembly is based on contig and scaffold N50 values of the genome sequence are 3.77 Mb and 20.41 Mb, respectively. Nine‐seven percent of the assembled sequences was anchored onto 13 chromosomes. Based on BUSCO analysis, the assembly involved 96.9% of conserved arthropod and 98.5% of the conserved Hemiptera single‐copy orthologous genes. A total of 14,089 protein‐coding genes were predicted. Phylogenetic analysis revealed that S. chinensis diverged from the common ancestor of Eriosoma lanigerum approximately 57 million years ago (MYA). In addition, 35 genes encoding salivary gland proteins showed differentially when S. chinensis forms a gall, suggesting they have potential roles in gall formation and plant defense suppression. Taken together, this high‐quality S. chinensis genome assembly and annotation provide a solid genetic foundation for future research to reveal the mechanism of gall formation and to explore the interaction between aphids and their host plants.     

Featured organism: reductive evolution in bacteria: Buchnera sp., Rickettsia prowazekii and Mycobacterium leprae

Obligate intracellular bacteria commonly have much reduced genome sizes compared to their nearest free-living relatives. One reason for this is reductive evolution: the loss of genes rendered non-essential due to the intracellular habitat. This can occur because of the presence of orthologous genes in the host, combined with the ability of the bacteria to import the protein or metabolite products of the host genes. In this article we take a look at three such bacteria whose genomes have been fully sequenced. Buchnera is an endosymbiont of the pea aphid, Acyrthosiphon pisum, the relationship between these two organisms being so essential that neither can reproduce in the absence of the other. Rickettsia prowazekii is the causative agent of louse-borne typhus in humans and Mycobacterium leprae infection of humans leads to leprosy. Both of these human pathogens have fastidious growth requirements, which has made them very difficult to study.