Comparative analysis of the Acyrthosiphon pisum genome and expressed sequence tag-based gene sets from other aphid species

Genomic Identification and Functional Analysis of JHAMTs in the Pond Wolf Spider, Pardosa pseudoannulata

Juvenile hormone (JH) plays a critical role in many physiological activities of Arthropoda. Juvenile hormone acid methyltransferase (JHAMT) is involved in the last steps of JH biosynthesis as an important rate-limiting enzyme. In recent studies, an increasing number of JHAMTs were identified in arthropods, but no JHAMT was reported in spiders. Herein, eight JHAMTs were identified in the pond wolf spider, Pardosa pseudoannulata, all containing the well conserved S-adenosyl-L-methionine binding motif. JHAMT-1 and the other seven JHAMTs were located at chromosome 13 and chromosome 1, respectively. Multiple alignment and phylogenetic analysis showed that JHAMT-1 was grouped together with insect JHAMTs independently and shared high similarities with insect JHAMTs compared to the other seven JHAMTs. In addition, JHAMT-1, JHAMT-2, and JHAMT-3 were highly expressed in the abdomen of spiderlings and could respond to the stimulation of exogenous farnesoic acid. Meanwhile, knockdown of these three JHAMTs caused the overweight and accelerated molting of spiderlings. These results demonstrated the cooperation of multi-JHAMTs in spider development and provided a new evolutionary perspective of the expansion of JHAMT in Arachnida.

Research advances in plant-insect molecular interaction

Acute and precise signal perception and transduction are essential for plant defense against insects. Insect elicitors-that is, the biologically active molecules from insects' oral secretion (which contains regurgitant and saliva), frass, ovipositional fluids, and the endosymbionts-are recognized by plants and subsequently induce a local or systematic defense response. On the other hand, insects secrete various types of effectors to interfere with plant defense at multiple levels for better adaptation. Jasmonate is a main regulator involved in plant defense against insects and integrates with multiple pathways to make up the intricate defense network. Jasmonate signaling is strictly regulated in plants to avoid the hypersensitive defense response and seems to be vulnerable to assault by insect effectors at the same time. Here, we summarize recently identified elicitors, effectors, and their target proteins in plants and discuss their underlying molecular mechanisms.

Transcriptome Analysis of Green Peach Aphid (Myzus persicae): Insight into Developmental Regulation and Inter-Species Divergence

Green peach aphid (Myzus persicae) and pea aphid (Acyrthosiphon pisum) are two phylogenetically closely related agricultural pests. While pea aphid is restricted to Fabaceae, green peach aphid feeds on hundreds of plant species from more than 40 families. Transcriptome comparison could shed light on the genetic factors underlying the difference in host range between the two species. Furthermore, a large scale study contrasting gene expression between immature nymphs and fully developed adult aphids would fill a previous knowledge gap. Here, we obtained transcriptomic sequences of green peach aphid nymphs and adults, respectively, using Illumina sequencing technology. A total of 2244 genes were found to be differentially expressed between the two developmental stages, many of which were associated with detoxification, hormone production, cuticle formation, metabolism, food digestion, and absorption. When searched against publically available pea aphid mRNA sequences, 13,752 unigenes were found to have no homologous counterparts. Interestingly, many of these unigenes that could be annotated in other databases were involved in the "xenobiotics biodegradation and metabolism" pathway, suggesting the two aphids differ in their adaptation to secondary metabolites of host plants. Conversely, 3989 orthologous gene pairs between the two species were subjected to calculations of synonymous and nonsynonymous substitutions, and 148 of the genes potentially evolved in response to positive selection. Some of these genes were predicted to be associated with insect-plant interactions. Our study has revealed certain molecular events related to aphid development, and provided some insight into biological variations in two aphid species, possibly as a result of host plant adaptation.

Insect-induced effects on plants and possible effectors used by galling and leaf-mining insects to manipulate their host-plant

Gall-inducing insects are iconic examples in the manipulation and reprogramming of plant development, inducing spectacular morphological and physiological changes of host-plant tissues within which the insect feeds and grows. Despite decades of research, effectors involved in gall induction and basic mechanisms of gall formation remain unknown. Recent research suggests that some aspects of the plant manipulation shown by gall-inducers may be shared with other insect herbivorous life histories. Here, we illustrate similarities and contrasts by reviewing current knowledge of metabolic and morphological effects induced on plants by gall-inducing and leaf-mining insects, and ask whether leaf-miners can also be considered to be plant reprogrammers. We review key plant functions targeted by various plant reprogrammers, including plant-manipulating insects and nematodes, and functionally characterize insect herbivore-derived effectors to provide a broader understanding of possible mechanisms used in host-plant manipulation. Consequences of plant reprogramming in terms of ecology, coevolution and diversification of plant-manipulating insects are also discussed.

Circulative, "nonpropagative" virus transmission: an orchestra of virus-, insect-, and plant-derived instruments

Species of plant viruses within the Luteoviridae, Geminiviridae, and Nanoviridae are transmitted by phloem-feeding insects in a circulative, nonpropagative manner. The precise route of virus movement through the vector can differ across and within virus families, but these viruses all share many biological, biochemical, and ecological features. All share temporal and spatial constraints with respect to transmission efficiency. The viruses also induce physiological changes in their plant hosts resulting in behavioral changes in the insects that optimize the transmission of virus to new hosts. Virus proteins interact with insect, endosymbiont, and plant proteins to orchestrate, directly and indirectly, virus movement in insects and plants to facilitate transmission. Knowledge of these complex interactions allows for the development of new tools to reduce or prevent transmission, to quickly identify important vector populations, and to improve the management of these economically important viruses affecting agricultural and natural plant populations.

Comparative transcriptomic analyses revealed divergences of two agriculturally important aphid species

BACKGROUND

Grain aphid (Sitobion avenae F) and pea aphid (Acyrthosiphon pisum) are two agriculturally important pest species, which cause significant yield losses to crop plants each year by inflicting damage both through the direct effects of feeding and by vectoring debilitating plant viruses. Although a close phylogenetic relationship between grain aphid and pea aphid was proposed, the biological variations between these two aphid species are obvious. While the host ranges of grain aphid is restricted to cereal crops and in particular wheat, that of pea aphid is wider, mainly colonizing leguminous plant species. Until now, the genetic factors underlying the divergence between grain aphid and pea aphid still remain unclear due to the limited genomic data of grain aphid available in public databases.

RESULTS

Based on a set of transcriptome data of grain aphid generated by using Roche 454 GS-FLX pyrosequencing, comparative analysis between this set of transcriptome data of grain aphid and mRNA sequences of pea aphid available in the public databases was performed. Compared with mRNA sequences of pea aphid, 4,857 unigenes were found to be specifically presented in the transcriptome of grain aphid under the rearing conditions described in this study. Furthermore, 3,368 orthologous pairs which could be calculated with both nonsynonymous (Ka) and synonymous (Ks) substitutions were used to infer their sequence divergences. The average differences in the coding, 5' and 3' untranslated regions of these orthologs were 10.53%, 21.29% and 18.96%, respectively. Moreover, of 340 orthologs which were identified to have evolved in response to positive selection based on the rates of Ka and Ks substitutions, 186 were predicted to be involved in secondary metabolism and xenobiotic metabolisms which might contribute to the divergence of these two aphid species.

CONCLUSIONS

The comprehensive transcriptome divergent sequence analysis between grain aphid and pea aphid provides an invaluable resource for the investigation of genes involved in host plant adaptation and evolution. Moreover, the demonstration of divergent transcriptome sequences between grain aphid and pea aphid pave the way for the investigation of the molecular mechanisms underpinning the biological variations of these two agriculturally important aphid species.

Accelerated evolution of morph-biased genes in pea aphids

Phenotypic plasticity, the production of alternative phenotypes (or morphs) from the same genotype due to environmental factors, results in some genes being expressed in a morph-biased manner. Theoretically, these morph-biased genes experience relaxed selection, the consequence of which is the buildup of slightly deleterious mutations at these genes. Over time, this is expected to result in increased protein divergence at these genes between species and a signature of relaxed purifying selection within species. Here we test these theoretical expectations using morph-biased genes in the pea aphid, a species that produces multiple morphs via polyphenism. We find that morph-biased genes exhibit faster rates of evolution (in terms of dN/dS) relative to unbiased genes and that divergence generally increases with increasing morph bias. Further, genes with expression biased toward rarer morphs (sexual females and males) show faster rates of evolution than genes expressed in the more common morph (asexual females), demonstrating that the amount of time a gene spends being expressed in a morph is associated with its rate of evolution. And finally, we show that genes expressed in the rarer morphs experience decreased purifying selection relative to unbiased genes, suggesting that it is a relaxation of purifying selection that contributes to their faster rates of evolution. Our results provide an important empirical look at the impact of phenotypic plasticity on gene evolution.

Juvenile hormone biosynthesis gene expression in the corpora allata of honey bee (Apis mellifera L.) female castes

Juvenile hormone (JH) controls key events in the honey bee life cycle, viz. caste development and age polyethism. We quantified transcript abundance of 24 genes involved in the JH biosynthetic pathway in the corpora allata-corpora cardiaca (CA-CC) complex. The expression of six of these genes showing relatively high transcript abundance was contrasted with CA size, hemolymph JH titer, as well as JH degradation rates and JH esterase (jhe) transcript levels. Gene expression did not match the contrasting JH titers in queen and worker fourth instar larvae, but jhe transcript abundance and JH degradation rates were significantly lower in queen larvae. Consequently, transcriptional control of JHE is of importance in regulating larval JH titers and caste development. In contrast, the same analyses applied to adult worker bees allowed us inferring that the high JH levels in foragers are due to increased JH synthesis. Upon RNAi-mediated silencing of the methyl farnesoate epoxidase gene (mfe) encoding the enzyme that catalyzes methyl farnesoate-to-JH conversion, the JH titer was decreased, thus corroborating that JH titer regulation in adult honey bees depends on this final JH biosynthesis step. The molecular pathway differences underlying JH titer regulation in larval caste development versus adult age polyethism lead us to propose that mfe and jhe genes be assayed when addressing questions on the role(s) of JH in social evolution.

Aphid protein effectors promote aphid colonization in a plant species-specific manner

Microbial pathogens and pests produce effectors to modulate host processes. Aphids are phloem-feeding insects, which introduce effectors via saliva into plant cells. However, it is not known if aphid effectors have adapted to modulate processes in specific plant species. Myzus persicae is a polyphagous insect that colonizes Arabidopsis thaliana and Nicotiana benthamiana, while the pea aphid Acyrthosiphon pisum specializes on colonizing plant species of the family Fabaceae. We found that M. persicae reproduction increased on transgenic Arabidopsis, producing the M. persicae effectors C002, PIntO1 (Mp1), and PIntO2 (Mp2), whereas reproduction of M. persicae did not increase on Arabidopsis producing the A. pisum orthologs of these three proteins. Plant-mediated RNA interference experiments showed that c002- and PIntO2-silenced M. persicae produce less progeny on Arabidopsis and N. benthamiana than nonsilenced aphids. Orthologs of c002, PIntO1, and PIntO2 were identified in multiple aphid species with dissimilar plant host ranges. We revealed high nonsynonymous versus synonymous nucleotide substitution rates within the effector orthologs, indicating that the effectors are fast evolving. Application of maximum likelihood methods identified specific sites with high probabilities of being under positive selection in PIntO1, whereas those of C002 and PIntO2 may be located in alignment gaps. In support of the latter, a M. persicae c002 mutant without the NDNQGEE repeat region, which overlaps with an alignment gap in C002, does not promote M. persicae colonization on Arabidopsis. Taken together, these results provide evidence that aphid effectors are under positive selection to promote aphid colonization on specific plant species.

Evolutionary study of duplications of the miRNA machinery in aphids associated with striking rate acceleration and changes in expression profiles

Background

The sequencing of the genome of the pea aphid Acyrthosiphon pisum revealed an unusual expansion of the miRNA machinery, with two argonaute-1, two dicer-1 and four pasha gene copies. In this report, we have undertaken a deeper evolutionary analysis of the phylogenetic timing of these gene duplications and of the associated selective pressures by sequencing the two copies of ago-1 and dcr-1 in different aphid species of the subfamily Aphidinae. We have also carried out an analysis of the expression of both copies of ago-1 and dcr-1 by semi-quantitative PCR in different morphs of the pea aphid life cycle.

Results

The analysis has shown that the duplication of ago-1 occurred in an ancestor of the subfamily Aphidinae while the duplication of dcr-1 appears to be more recent. Besides, it has confirmed a pattern of one conserved copy and one accelerated copy for both genes, and has revealed the action of positive selection on several regions of the fast-evolving ago-1b. On the other hand, the semi-quantitative PCR experiments have revealed a differential expression of these genes between the morphs of the parthenogenetic and the sexual phases of Acyrthosiphon pisum.

Conclusions

The discovery of these gene duplications in the miRNA machinery of aphids opens new perspectives of research about the regulation of gene expression in these insects. Accelerated evolution, positive selection and differential expression affecting some of the copies of these genes suggests the possibility of a neofunctionalization of these duplicates, which might play a role in the display of the striking phenotypic plasticity of aphids.

Deep sequencing of the transcriptomes of soybean aphid and associated endosymbionts

BACKGROUND

The soybean aphid has significantly impacted soybean production in the U.S. Transcriptomic analyses were conducted for further insight into leads for potential novel management strategies.

METHODOLOGY/PRINCIPAL FINDINGS

Transcriptomic data were generated from whole aphids and from 2,000 aphid guts using an Illumina GAII sequencer. The sequence data were assembled de novo using the Velvet assembler. In addition to providing a general overview, we demonstrate (i) the use of the Multiple-k/Multiple-C method for de novo assembly of short read sequences, followed by BLAST annotation of contigs for increased transcript identification: From 400,000 contigs analyzed, 16,257 non-redundant BLAST hits were identified; (ii) analysis of species distributions of top non-redundant hits: 80% of BLAST hits (minimum e-value of 1.0-E3) were to the pea aphid or other aphid species, representing about half of the pea aphid genes; (iii) comparison of relative depth of sequence coverage to relative transcript abundance for genes with high (membrane alanyl aminopeptidase N) or low transcript abundance; (iv) analysis of the Buchnera transcriptome: Transcripts from 57.6% of the genes from Buchnera aphidicola were identified; (v) identification of Arsenophonus and Wolbachia as potential secondary endosymbionts; (vi) alignment of full length sequences from RNA-seq data for the putative salivary gland protein C002, the silencing of which has potential for aphid management, and the putative Bacillus thuringiensis Cry toxin receptors, aminopeptidase N and alkaline phosphatase.

CONCLUSIONS/SIGNIFICANCE

THIS STUDY PROVIDES THE MOST COMPREHENSIVE DATA SET TO DATE FOR SOYBEAN APHID GENE EXPRESSION: This work also illustrates the utility of short-read transcriptome sequencing and the Multiple-k/Multiple-C method followed by BLAST annotation for rapid identification of target genes for organisms for which reference genome sequences are not available, and extends the utility to include the transcriptomes of endosymbionts.

Aphids: a model for polyphenism and epigenetics

Environmental conditions can alter the form, function, and behavior of organisms over short and long timescales, and even over generations. Aphid females respond to specific environmental cues by transmitting signals that have the effect of altering the development of their offspring. These epigenetic phenomena have positioned aphids as a model for the study of phenotypic plasticity. The molecular basis for this epigenetic inheritance in aphids and how this type of inheritance system could have evolved are still unanswered questions. With the availability of the pea aphid genome sequence, new genomics technologies, and ongoing genomics projects in aphids, these questions can now be addressed. Here, we review epigenetic phenomena in aphids and recent progress toward elucidating the molecular basis of epigenetics in aphids. The discovery of a functional DNA methylation system, functional small RNA system, and expanded set of chromatin modifying genes provides a platform for analyzing these pathways in the context of aphid plasticity. With these tools and further research, aphids are an emerging model system for studying the molecular epigenetics of polyphenisms.

Juvenile hormone titre and related gene expression during the change of reproductive modes in the pea aphid

Most aphids show reproductive polyphenism, i.e. they alternate their reproductive modes from parthenogenesis to sexual reproduction in response to short photoperiods. Although juvenile hormone (JH) has been considered a likely candidate for regulating the transition from asexual to sexual reproduction after photoperiod sensing, there are few studies investigating the direct relationship between JH titres and the reproductive-mode change. In addition, the sequencing of the pea aphid genome has allowed identification of the genes involved in the JH pathway, which in turn allows us to examine their expression levels in relation to the reproductive-mode change. Using liquid chromatography-mass spectrometry in the pea aphid, JHIII titre was shown to be lower in aphids producing sexual morphs under short-day conditions than in aphids producing parthenogenetic morphs under long-day conditions. The expression levels of genes upstream and downstream of JH action were quantified by real-time quantitative reverse-transcription-PCR across the reproductive-mode change. The expression level of JH esterase, which is responsible for JH degradation, was significantly higher in aphids reared under short-day conditions. This suggests that the upregulation of the JH degradation pathway may be responsible for the lower JHIII titre in aphids exposed to short-days, leading to the production of sexual morphs.

Comparison of gene repertoires and patterns of evolutionary rates in eight aphid species that differ by reproductive mode

In theory, the loss of sexual reproduction is expected to result in the accumulation of deleterious mutations. In aphids, two main types of life cycle, cyclic and obligate parthenogenesis, represent respectively “sexual” and “asexual” reproductive modes. We used the complete pea aphid genome and previously published expressed sequence tags (ESTs) from two other aphid species. In addition, we obtained 100,000 new ESTs from five more species. The final set comprised four sexual and four asexual aphid species and served to test the influence of the reproductive mode on the evolutionary rates of genes. We reconstructed coding sequences from ESTs and annotated these genes, discovering a novel peptide gene family that appears to be among the most highly expressed transcripts from several aphid species. From 203 genes found to be 1:1 orthologs among the eight species considered, we established a species tree that partly conflicted with taxonomy (for Myzus ascalonicus). We then used this topology to evaluate the dynamics of evolutionary rates and mutation accumulation in the four sexual and four asexual taxa. No significant increase of the nonsynonymous to synonymous ratio or of nonsynonymous mutation numbers was found in any of the four branches for asexual taxa. We however found a significant increase of the synonymous rate in the branch leading to the asexual species Rhopalosiphum maidis, which could be due to a change in the mutation rate or to an increased number of generations implied by its change of life cycle.

Genomics of environmentally induced phenotypes in 2 extremely plastic arthropods

Understanding how genes and the environment interact to shape phenotypes is of fundamental importance for resolving important issues in adaptive evolution. Yet, for most model species with mature genetics and accessible genomic resources, we know little about the natural environmental factors that shape their evolution. By contrast, animal species with deeply understood ecologies and well characterized responses to environmental cues are rarely subjects of genomic investigations. Here, we preview advances in genomics in aphids and waterfleas that may help transform research on the regulatory mechanisms of phenotypic plasticity. This insect and crustacean duo has the capacity to produce extremely divergent phenotypes in response to environmental stimuli. Sexual fate and reproductive mode are condition-dependent in both groups, which are also capable of altering morphology, physiology and behavior in response to biotic and abiotic cues. Recently, the genome sequences for the pea aphid Acyrthosiphon pisum and the waterflea Daphnia pulex were described by their respective research communities. We propose that an integrative study of genome biology focused on the condition-dependent transcriptional basis of their shared plastic traits and specialized mode of reproduction will provide broad insight into adaptive plasticity and genome by environment interactions. We highlight recent advances in understanding the genome regulation of alternative phenotypes and environmental cue processing, and we propose future research avenues to discover gene networks and epigenetic mechanisms underlying phenotypic plasticity.

Diversification of genes for carotenoid biosynthesis in aphids following an ancient transfer from a fungus

The pea aphid genome was recently found to harbor genes for carotenoid biosynthesis, reflecting an ancestral transfer from a fungus. To explore the evolution of the carotene desaturase gene family within aphids, sequences were retrieved from a set of 34 aphid species representing numerous deeply diverging lineages of aphids and analyzed together with fungal sequences retrieved from databases. All aphids have at least one copy of this gene and some aphid species have up to seven, whereas fungal genomes consistently have a single copy. The closest relatives of aphids, adelgids, also have carotene desaturase; these sequences are most closely related to those from aphids, supporting a shared origin from a fungal to insect transfer predating the divergence of adelgids and aphids. Likewise, all aphids, and adelgids, have carotenoid profiles that are consistent with their biosynthesis using the acquired genes of fungal origin rather than derivation from food plants. The carotene desaturase was acquired from a fungal species outside of Ascomycota or Basidiomycota and closest to Mucoromycotina among sequences available in databases. In aphids, an ongoing pattern of gene duplication is indicated by the presence of both anciently and recently diverged paralogs within genomes and by the presence of a high frequency of pseudogenes that appear to be recently inactivated. Recombination among paralogs is evident, making analyses of patterns of selection difficult, but tests of selection for a nonrecombining region indicates that duplications tend to be followed by bouts of positive selection. Species of Macrosiphini, which often show color polymorphisms, typically have a larger number of desaturase copies relative to other species sampled in the study. These results indicate that aphid evolution has been accompanied by ongoing evolution of carotenogenic genes, which have undergone duplication, recombination, and occasional positive selection to yield a wide variety of carotenoid profiles in different aphid species.

Effector proteins that modulate plant--insect interactions

Insect herbivores have highly diverse life cycles and feeding behaviors. They establish close interactions with their plant hosts and suppress plant defenses. Chewing herbivores evoke characteristic defense responses distinguishable from general mechanical damage. In addition, piercing-sucking hemipteran insects display typical feeding behavior that suggests active suppression of plant defense responses. Effectors that modulate plant defenses have been identified in the saliva of these insects. Tools for high-throughput effector identification and functional characterization have been developed. In addition, in some insect species it is possible to silence gene expression by RNAi. Together, this technological progress has enabled the identification of insect herbivore effectors and their targets that will lead to the development of novel strategies for pest resistances in plants.

Genomic identification of a putative circadian system in the cladoceran crustacean Daphnia pulex

Essentially nothing is known about the molecular underpinnings of crustacean circadian clocks. The genome of Daphnia pulex, the only crustacean genome available for public use, provides a unique resource for identifying putative circadian proteins in this species. Here, the Daphnia genome was mined for putative circadian protein genes using Drosophila melanogaster queries. The sequences of core clock (e.g. CLOCK, CYCLE, PERIOD, TIMELESS and CRYPTOCHROME 2), clock input (CRYPTOCHROME 1) and clock output (PIGMENT DISPERSING HORMONE RECEPTOR) proteins were deduced. Structural analyses and alignment of the Daphnia proteins with their Drosophila counterparts revealed extensive sequence conservation, particularly in functional domains. Comparisons of the Daphnia proteins with other sequences showed that they are, in most cases, more similar to homologs from other species, including vertebrates, than they are to those of Drosophila. The presence of both CRYPTOCHROME 1 and 2 in Daphnia suggests the organization of its clock may be more similar to that of the butterfly Danaus plexippus than to that of Drosophila (which possesses CRYPTOCHROME 1 but not CRYPTOCHROME 2). These data represent the first description of a putative circadian system from any crustacean, and provide a foundation for future molecular, anatomical and physiological investigations of circadian signaling in Daphnia.

Tangible benefits of the aphid Acyrthosiphon pisum genome sequencing for aphid proteomics: Enhancements in protein identification and data validation for homology-based proteomics

Homology-driven proteomics promises to reveal functional biology in insects with sparse genome sequence information. A proteomics study comparing plant virus transmission competent and refractive genotypes of the aphid Schizaphis graminum isolated numerous candidate proteins involved in virus transmission, but limited genome sequence information hampered their identification. The complete genome of the pea aphid, Acyrthosiphon pisum, released in 2008, enabled us to double the number of protein identifications beyond what was possible using available EST libraries and other insect sequences. This was concomitant with a dramatic increase of the number of MS and MS/MS peptide spectra matching the genome-derived protein sequence. LC-MS/MS proved to be the most robust method of peptide detection. Cross-matching spectral data to multiple EST sequences and error tolerant searching to identify amino acid substitutions enhanced the percent coverage of the Schizaphis graminum proteins. 2-D electrophoresis provided the protein pI and MW which enabled the refinement of the candidate protein selection and provided a measure of protein abundance when coupled to the spectral data. Thus, the homology-based proteomics pipeline for insects should include efforts to maximize the number of peptide matches to the protein to increase certainty in protein identification and relative protein abundance.

Combining next-generation sequencing strategies for rapid molecular resource development from an invasive aphid species, Aphis glycines

Background

Aphids are one of the most important insect taxa in terms of ecology, evolutionary biology, genetics and genomics, and interactions with endosymbionts. Additionally, many aphids are serious pest species of agricultural and horticultural plants. Recent genetic and genomic research has expanded molecular resources for many aphid species, including the whole genome sequencing of the pea aphid, Acrythosiphon pisum. However, the invasive soybean aphid, Aphis glycines, lacks in any significant molecular resources.

Methodology/Principal Findings

Two next-generation sequencing technologies (Roche-454 and Illumina GA-II) were used in a combined approach to develop both transcriptomic and genomic resources, including expressed genes and molecular markers. Over 278 million bp were sequenced among the two methods, resulting in 19,293 transcripts and 56,688 genomic sequences. From this data set, 635 SNPs and 1,382 microsatellite markers were identified. For each sequencing method, different soybean aphid biotypes were used which revealed potential biotype specific markers. In addition, we uncovered 39,822 bp of sequence that were related to the obligatory endosymbiont, Buchnera aphidicola, as well as sequences that suggest the presence of Hamiltonella defensa, a facultative endosymbiont.

Conclusions and Significance

Molecular resources for an invasive, non-model aphid species were generated. Additionally, the power of next-generation sequencing to uncover endosymbionts was demonstrated. The resources presented here will complement ongoing molecular studies within the Aphididae, including the pea aphid whole genome, lead to better understanding of aphid adaptation and evolution, and help provide novel targets for soybean aphid control.

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.