- Open Access
Immunity and other defenses in pea aphids, Acyrthosiphon pisum
Genome Biologyvolume 11, Article number: R21 (2010)
Recent genomic analyses of arthropod defense mechanisms suggest conservation of key elements underlying responses to pathogens, parasites and stresses. At the center of pathogen-induced immune responses are signaling pathways triggered by the recognition of fungal, bacterial and viral signatures. These pathways result in the production of response molecules, such as antimicrobial peptides and lysozymes, which degrade or destroy invaders. Using the recently sequenced genome of the pea aphid (Acyrthosiphon pisum), we conducted the first extensive annotation of the immune and stress gene repertoire of a hemipterous insect, which is phylogenetically distantly related to previously characterized insects models.
Strikingly, pea aphids appear to be missing genes present in insect genomes characterized to date and thought critical for recognition, signaling and killing of microbes. In line with results of gene annotation, experimental analyses designed to characterize immune response through the isolation of RNA transcripts and proteins from immune-challenged pea aphids uncovered few immune-related products. Gene expression studies, however, indicated some expression of immune and stress-related genes.
The absence of genes suspected to be essential for the insect immune response suggests that the traditional view of insect immunity may not be as broadly applicable as once thought. The limitations of the aphid immune system may be representative of a broad range of insects, or may be aphid specific. We suggest that several aspects of the aphid life style, such as their association with microbial symbionts, could facilitate survival without strong immune protection.
Aphids face numerous environmental challenges, including infection by diverse pathogens and parasites. These pressures include parasitoid wasps, which consume their hosts as they develop inside, and a variety of viral, bacterial and fungal pathogens. Both parasitoid wasp and fungal pathogens cause significant decline of natural aphid populations [1, 2], and have been suggested as potential agents for biocontrol of these agriculturally destructive pests. While facing such challenges, aphids also cope with predators and abiotic stresses, such as extreme temperature fluctuations. Thus, like most insects, aphids must attempt to survive in a harsh, complex environment.
Insects have a number of defense mechanisms. First, many insects, including aphids, behaviorally avoid predators, pathogens, and environmental stressors [3–6]. When stressors cannot be avoided, insects have a protective cuticle and gut pH inhospitable to many foreign organisms. If these barriers fail, immunological defense mechanisms recognize the invader, triggering a signaling cascade and response. While insects do not have adaptive, antigen-based responses typical of vertebrates, insects do have innate immune responses, which include clotting, phagocytosis, encapsulation, and production of antimicrobial substances [7, 8]. Phagocytosis and encapsulation are referred to as cellular responses as they are mediated by blood cells . Reponses vary depending on the invader, with antimicrobial peptides being central to combating microbes and encapsulation being central to combating larger invaders, such as parasitoids. Until recently, it was presumed that insects were limited to these non-specific innate immune responses and had no specific immunity (for example, the antigen-based immune response of humans). There is, however, increasing evidence for the ability of insects to mount specific immune responses .
Here we focus on the identification of aphid genes that are known to play a role in the recognition and degradation of microbial pathogens in other insects, as these are the invertebrate defense processes that are best understood. In the fruit fly Drosophila melanogaster, recognition of an invasive microbe leads to signal production via four pathways (Toll, immunodeficiency (IMD), c-Jun N-terminal kinase (JNK), and Janus kinase/Signal transducers and activators of transcription (JAK/STAT)) . Each pathway is activated in response to particular pathogens . Signaling triggers the production of a multitude of effectors, including, most notably, antimicrobial peptides (AMPs). Insect AMPs may be 1,000-fold induced in microbe-challenged insects compared to basal levels. In insect genomes annotated to date, these pathways appear well conserved, with most of the key components found across flies (Drosophila spp.), mosquitoes (Aedes aegypti, Anopheles gambiae), bees (Apis mellifera) and beetles (Tribolium castaneum) [13–17].
Because aphids and other insects face diverse challenges, we propose models for several genes critical to other elements of insect stress responses. These include genes encoding heat shock proteins (HSPs), which are synthesized in almost all living organisms when exposed to high temperatures or stress . We also suggest models for genes involved in the synthesis of the alarm pheromone (E)-β farnesene, which aphids release in the presence of predators . While there are undoubtedly many other genes involved in stress and immunological responses, our selection of genes for exploration provides a broad survey of the known insect immune and stress repertoire and will serve as a basis for future exploration of more specific responses.
The pea aphid genome provides novel insights into arthropod immunity for two reasons. First, most of our understanding of insect immune and stress responses comes from holometabolous insects, the group of insects with complete metamorphisis, such as flies, butterflies, beetles and bees. The genome of the hemimetabolous pea aphid, Acyrthosiphon pisum, may thus provide novel insight into immunity and defense in more basal, non-holometabolous insects, which have incomplete metamorphisis. Second, aphids are unique amongst the arthropods sequenced to date in that they are intimately dependent on both obligate and facultative bacterial symbionts for their survival. The aphid symbiont community includes Buchnera aphidicola, obligate and intracellular Gram-negative bacteria that have the ability to synthesize required amino acids not readily available in the aphid diet. Beyond this obligate symbiosis, aphids frequently host one or more additional Gram-negative bacterial symbionts, including most notably Hamiltonella defensa, Serratia symbiotica and Regiella insecticola [20, 21]. Unlike Buchnera, which is present in all aphids and is thus considered a primary symbiont, these bacteria are considered to be facultative, secondary symbionts, because their presence varies within an aphid species . Secondary symbiotic bacteria have been shown to influence several aspects of aphid ecology, including heat tolerance and resistance to parasites and pathogens [23–26]. Specifically, both H. defensa and S. symbiotica confer protection against parasitoid wasp development [27, 28], and R. insecticola decreases A. pisum mortality after exposure to the fungal pathogen Pandora neoaphidis . These are some of the best-studied examples of symbiont-conferred protection .
Aphids thus provide an excellent opportunity to study the immune system of an organism that is dependent on microbial symbionts but is hampered by parasites and pathogens. Despite this, little work has been done to characterize the aphid immune response. Altincicek et al.  found that compared to other insects, stabbing a pea aphid with bacteria elicits reduced lysozyme-like (muramidase) activity, and no detectable activity against live bacteria in hemolymph assays. Furthermore, suppression subtraction hybridization (SSH) of bacterial-challenged aphids uncovered no antimicrobial peptides and few genes of known immune function . These results are surprising given that similar studies in other insects demonstrate that antimicrobial peptide production and upregulation of immune-related genes is a common feature of the insect immune response that can be captured in functional assays such as SSH [32–35]. This suggests that aphids have a significantly reduced or altered immune repertoire.
Using the recently sequenced genome of the pea aphid clone LSR1, in this study, we take two approaches to study immunity and stress in pea aphids. First, we assay presence/absence of a subset of known immune and stress-related genes. Second, we combine functional assays targeting the production of RNA and proteins to gain insight into how pea aphids respond to various challenges. Overall, our results suggest that pea aphids are missing many genes central to immune function in other insects, and that, although pea aphids do mount some response to challenges, the overall immune-response of pea aphids is more limited than that of other insects studied to date.
Results and discussion
Overview of annotation
We focused our manual annotation efforts on a subset of genes involved in the innate, humoral immune response contributing to recognition, signaling and response to bacteria and fungi in arthropods. We also manually annotated some genes involved in more general stress responses (for example, HSPs). All annotations are based on the recently completed sequencing of pea aphid clone LSR1 . All genes manually annotated, as well as those genes that we found to be missing in the pea aphid genome, are listed in Table S1 in Additional file 1. Also in this table, BLAST-based searches revealed that another aphid, Myzus persicae (green peach aphid), has putative homologs for many immune and stress related genes identified in the pea aphid.
Annotation of microbial recognition genes
Peptidoglycan receptor proteins
Upon microbial invasion, Drosophila utilize several pathogen recognition receptors (PRRs) to detect pathogen-specific molecular patterns (for example, cell-surface motifs) . PRRs include peptidoglycan receptor proteins (PGRPs), which recognize peptidoglycans present in cell walls of Gram-positive and Gram-negative bacteria. PGRP-based recognition activates both the Toll and IMD/JNK pathways. PGRPs are highly conserved, with mammals and insect PGRPs sharing a 160 amino acid domain [38, 39]. Thus, it is surprising that pea aphids, in contrast to all other sequenced insects, appear to have no PGRPs. One other sequenced arthropod, the crustacean Daphia pulex, is also missing PGRPs .
Gram-negative binding proteins
GNBPs (Gram-negative binding proteins, a historical misnomer) are thought to detect Gram-positive bacteria . GNBPs and PGRPs are suspected to form a complex. GNBPs then hydrolyze Gram-positive peptidoglycans into small fragments, which are detected by PGRPs [41, 42]. Aphids have two GNBP paralogs, GNBP1 and GNBP2 (see Figure S1a in Additional file 1). Because GNBPs are thought to form a complex with PGRPs, the presence of GNBPs without PGRPs in aphids, as well as in the crustacean D. pulex , calls into question whether GNBPs play a role in bacterial detection in these organisms. Some GNBPs and similar proteins are known to function in fungal recognition , which may be the primary function of these molecules in aphids.
Lectins are a diverse group of sugar binding proteins. Many lectins function in insect immune recognition by binding to polysaccharide chains on the surface of pathogens . Drosophila c-type lectins also appear to facilitate encapsulation of parasitoid invaders, by marking surfaces for hemocyte recruitment . Aphids have five c-type lectin paralogs.
Galectins are another widely-distributed group of lectins . In mosquitoes, galectins are upregulated in response to both bacterial and malaria parasite infection [46, 47]. Insect galectins are thought to be involved in either pathogen recognition, via recognition of β-galactoside, or in phagocytosis . Aphids have two galectin paralogs.
Class C scavenger receptors
In Drosophila, Scavenger receptors exhibit broad affinity towards both Gram-positive and Gram-negative bacteria, but not yeast . Pathogen recognition by class C scavenger receptors in Drosophila facilitates phagocytosis, and natural genetic variation of Drosophila scavenger receptors is correlated with variation in the ability to suppress bacterial infection . While D. melanogaster has four class C scavenger receptor homologs, A. gambiae and A. mellifera have only one. Pea aphids appear to have no class C scavenger receptors.
The Nimrod superfamily and Dscam
Several members of the Nimrod superfamily appear to function as receptors in phagocytosis and bacterial binding [50, 51]. Such insect genes include eater and nimrod. Many of these genes are characterized by a specific EGF (epidermal growth factor) repeat, and are duplicated in the genomes of D. melanogaster, T. castaneum and A. mellifera . We were unable to identify any EGF motif genes in the pea aphid genome.
Complex alternative splicing of Dscam (Down syndrome cell adhesion molecule) generates diverse surface receptors sometimes employed in arthropod innate immune defenses [53–55]. Though we did not manually annotate this complex gene as a part of this initial aphid immune gene project, we did identify multiple predicted protein sequences coded by the aphid genome with strong similarity to Dscam in other insects [GenBank: XP_001951010, XP_001949262, XP_001945921, XP_001951684, XP_001942542]. Further investigations will be necessary to determine the activity and hypervariability of these genes and their transcripts in aphids.
Annotation of signaling pathways
The Toll signaling pathway
The Toll pathway is a signaling cascade involved in both development and innate immunity. In Drosophila, deletion of many of the component genes leads to increased susceptibility to many Gram-positive bacteria and fungal pathogens , and some Gram-negative bacteria and viruses . In addition, upregulation of many components of the Toll pathway is observed following parasitoid wasp invasion . The Toll pathway appears to be intact in pea aphids. We found convincing matches for genes encoding the extracellular cytokine spätzle, the transmembrane receptor Toll, the tube and MyD88 adaptors, the kinase pelle, the inhibitor molecule cactus (a homolog of IkB), cactin, Pellino, Traf, and the transactivator dorsal (Figure 1). The latter two genes are duplicated.
As in other insects, there are several gene families associated with the Toll pathway that are represented in aphids. First, aphids seem to have multiple spätzles that segregate with Drosophila spätzles 1, 2, 3, 4 and 6 in phylogenetic analyses (Figure S1b in Additional file 1). Second, aphids also have a suite of serine proteases and serine protease inhibitors (serpins). Though we did not manually annotate serine proteases and serpins as a part of this initial aphid immune gene project, we did identify multiple predicted protein sequences in the aphid genome with strong similarity to serine proteases and serpins in other insects. In insects, these molecules function in digestion, embryonic development and defense responses towards both microbial and parasitoid wasp invaders [57–59]. In the absence of microbial challenge, the serpin necrotic prevents activation of the Toll pathway, but upon immunological challenge, the Toll pathway is triggered by a cascade of serine proteases, including persephone, which is thought to be specific to fungal challenge . Though it is not clear which of the many aphid serine proteases is homologous to persephone, it is likely that pea aphids have serine proteases capable of triggering the Toll pathway. Finally, aphids also have multiple genes encoding Toll receptors, which function as transmembrane receptors in both mammals and insects. While nine single-copy Toll genes have been identified in D. melanogaster (Toll1 to Toll9), it seems that pea aphids, like other insects, lack some of these genes, but have multiple copies of others (Figure S1c in Additional file 1). In other organisms, some, but not all, Tolls serve a role in immune function, while others function in developmental processes [60–62]. For aphids, it is not yet clear what role each Toll serves.
The JAK/STAT signaling pathway
Like the Toll pathway, in Drosophila, the JAK/STAT pathway is involved in both development and immunity. The JAK/STAT pathway is the least understood of the core insect immune pathways. JAK/STAT pathway induction appears to lead to overproliferation of hemocytes, upregulation of thiolester-containing proteins (TEPs), and an antiviral response . Changes in gene expression following parasitoid wasp invasion of Drosophila larvae suggest a role for the JAK/STAT pathway in parasitoid response . Pea aphids have homologs of all core JAK/STAT genes, including genes encoding the cytokine receptor domeless, JAK tyrosine kinase (aka Hopscotch), and the STAT92E transcription factor (Figure 1). STAT92E appears to be duplicated. No homologs were found for upd (unpaired), considered a key ligand in Drosophila JAK/STAT induction. This ligand is also missing in other insects (for example, A. mellifera) .
IMD and JNK signaling pathways
Surprisingly, pea aphids appear to be missing many crucial components of the IMD signaling pathway. This pathway is critical for fighting Gram-negative bacteria in Drosophila [11, 64], and IMD pathway member knockouts influence susceptibility to some Gram-positive bacteria and fungi as well . IMD-associated genes missing in pea aphids include PGRPs (see above), IMD, dFADD, Dredd and Relish (Rel) (Figure 1). In contrast, conserved one-to-one orthologs of these same genes are found across Drosophila, Apis, Aedes, Anopheles and Tribolium . Cursory BLAST-based searches for these genes in other arthropods suggest that some may be missing (Figure 2). Pea aphids do have homologs for a few pathway members (TAB, TAK, kenny, Iap2 and IRD5; Figure 1).
While missing IMD-associated genes, pea aphids have plausible orthologs for most components of the JNK pathway (Figure 1). In Drosophila, the JNK pathway regulates many developmental processes, as well as wound healing , and has been proposed to play a role in antimicrobial peptide gene expression and cellular immune responses [11, 66]. Genes present include hep, basket, and JRA. Searches for homologs to the Drosophila kayak (kay) gene found an apparently similar transcription factor encoding gene in the A. pisum genome [GenBank: XP_001949014], but this match was largely restricted to the leucine zipper region, and failed tests of reciprocity.
The absence of IMD but presence of JNK in pea aphids is surprising as, in Drosophila, the IMD signaling pathway leads to activation of components of the JNK signaling pathway . Specifically, when TAK, a protein kinase of the IMD pathway, is activated, it triggers the JNK pathway. Whether TAK can be activated without the rest of the IMD pathway is unknown. An alternative IMD-independent activation of JNK, via the inducer Eiger , has been proposed in Drosophila . As Eiger is present in the pea aphid, this mode of activation may serve a critical role in any aphid JNK-based immune response.
Annotation of recognition genes
Introduction of microbes into most insects leads to the production of AMPs by the fat body, an insect immune-response tissue, and occasionally by hemocytes and other tissues [68–71]. These peptides are secreted into the hemolymph, where they exhibit a broad range of activities against fungi and bacteria. The mechanisms of AMP action are poorly understood, but at least in some cases (for example, drosomycin in Drosophila), AMPs destroy invading microbes by disrupting microbial cell membranes, leading to cell lysis .
Antimicrobial peptides are diverse and ubiquitous. They tend to be small molecules (<30 kDa) specialized at attacking particular microbial classes (that is, Gram-positive bacteria, fungi, and so on) [68, 69]. While some antimicrobial peptides are found in only a single insect group (for example, metchnikowin is found only in Drosophila), others are widely dispersed across eukaryotes (for example, defensins are present in fungi, plants and animals). Genomics, coupled with proteomics, has revealed that all sequenced insects, and many other insects, have multiple types of antimicrobial peptides (Figure 2). Pea aphids, surprisingly, are missing many of the antimicrobial peptides common to other insects. For example, while all insect genomes annotated thus far have genes encoding defensins , homology-based searches, phylogenetic-based analyses, transcriptomics (see below), and proteomics (see below) failed to find any signatures of defensins in the pea aphid genome. The presence of defensins in the human louse Pediculus humanus (Figure 2), and in the ancient apterygote insect, the fire brat Thermobia domestica , suggests that defensins have been lost during aphid evolution.
Extensive searches for genes encoding insect cecropins, drosocin (and other proline-rich arthropod AMPs), diptericin (and other glycine-rich AMPs), drosomycin, metchnikowin, formicin, moricin, spingerin, gomesin, tachyplesin, polyphemusin, andropin, gambicin, and virescein also revealed no hits. Weak hits were found for genes that encode for two antimicrobial peptides in other invertebrates: megourin [UniProtKB: P83417], originally isolated from another aphid species, the vetch aphid Megoura viciae (P Bulet et al., unpublished data) and penaeidin [UniProtKB: P81058], originally isolated from the shrimp Penaeus vannamei. The putative pea aphid Megourin (scaffold EQ11086, positions 45,752 to 45,892), however, is highly diverged from that of M. viciae (31% identity) and, compared to its M. viciae counterparts, seems to have a shorter carboxy-terminal region containing a stop-codon (Figure S2 in Additional file 1). Using three different primer pairs, we were unable to amplify products of this putative Megourin from cDNA generated for expression analyses (see below). The highly divergent Penaeidin [GenBank: ACYPI37769] (Figure S2 in Additional file 1) also did not amplify from cDNA.
We found six Thaumatin homologs in the A. pisum genome that show overall sequence and predicted structure similarities to plant thaumatins (Figure 3a, b). Thaumatin-like proteins are disulfide-bridged polypeptides of about 200 residues. Some thaumatins possess antifungal activity in plant tissues after infection . Recently, a thaumatin found in the beetle T. castaneum was shown to inhibit spore germination of the filamentous fungi Beauveria bassiana and Fusarium culmorum . Phylogenetic analyses revealed that A. pisum thaumatins form a monophyletic group closely related to beetle thaumatins (Figure 3c). Since thaumatin-like genes are conspicuously absent from the genomes of Drosophila, Apis, Anopheles, Pediculus and Ixodes (Figure 2), our findings indicate that thaumatins may represent ancient defense molecules that have been lost in several insect species, or have been independently acquired in aphids and beetles. The monophyly of aphid and beetle thaumatins provides no indication of an origin of novel acquisition (Figure 3c).
Lysozymes represent a family of enzymes that degrade bacterial cell walls by hydrolyzing the 1,4-beta-linkages between N-acetyl-D-glucosamine and N-acetylmuramic acid in peptidoglycan heteropolymers . They are ubiquitously distributed among living organisms and are believed to be essential for defense against bacterial infection. Lysozymes are classified into several types (that is, c (chicken), g (goose), i (invertebrate), plant, bacteria and phage types). C-type lysozymes are the most common for metazoa, being found in all vertebrates examined thus far and many invertebrates, including all the previously sequenced insects. For example, D. melanogaster and A. gambiae have at least seven and nine loci for c-type lysozymes, respectively [74, 75]. Insects also have i-type homologs, but their bacteriolytic activities are unclear .
Unlike other insects sequenced thus far, similarity searches demonstrated that A. pisum lacks genes for c-type lysozymes. The analysis further verified that the genome also lacks genes for g-type, plant-type, and phage-type lysozymes. Only three genes for i-type homologs were detected in the genome (Figure S1d in Additional file 1). One of them, Lys1, is highly expressed in the bacteriocyte . Two others, Lys2 and Lys3, are located adjacent to Lys1.
Notably, two genes that appear to have been transferred from bacterial genomes to the A. pisum genome encode bacteriolytic enzymes . One is for a chimeric protein that consists of a eukaryotic carboxypeptidase and a bacterial lysozyme. The other (AmiD) encodes N-acetylmuramoyl-L-alanine amidase, which is not a true lysozyme (1,4-beta-N-acetylmuramidase) but similarly degrades bacterial cell walls. While some of these bacteriolytic-related genes are highly expressed in the bacteriocyte, and lysozymes appear to be upregulated in response to some challenges (see gene expression study, below), assays of bacterioltyic activity of hemolymph from immune-challenged aphids suggest that aphid hemolymph has weak to no lysozyme-like activity . Further studies will determine the role of these gene products.
Chitinases are enzymes that degrade chitin (a long-chain polymer of N-acetyl-D-glucosamine), hydrolyzing 1,4-beta-linkages between N-acetyl-D-glucosamines. Chitinases and lysozymes represent a superfamily of hydrolases, and their catalytic activities are similar. Indeed, some chitinases show lysozyme activity and vice versa . In insects, chitinases are used to degrade the chitin in the exoskeleton and peritrophic membrane during molting, and some are suspected to have antifungal activity, as fungal cell walls also consist of chitin . Similarity searches followed by phylogenetic analyses demonstrated that the genome of A. pisum encodes seven genes for putative chitinase-like proteins . Further studies are required to determine the biochemical properties and substrate specificity of these chitinase-like proteins.
TEPs and Tots
Some TEPs can covalently attach to pathogens and parasites in order to 'mark' them for phagocytosis . Like other insects, aphids have multiple Tep paralogs. Both are homologous to TepIII (Figure S1e in Additional file 1). Homologs of TepI, TepII and TepIV were not found. In contrast, no Turandot (Tot) genes, which encode small peptides induced by severe stress and septic injury in Drosophila [81–83], have been found in aphids or in other insects other than Drosophila spp. Both TEPs and Tots are thought to be regulated by the JAK/STAT pathway.
Phenoloxidase-mediated melanin formation characteristically accompanies wound clotting, phagocytosis and encapsulation of pathogens and parasites . In insects, the inactive enzyme prophenoloxidase (ProPO) is activated by serine proteases to yield phenoloxidase . Aphids appear to have two prophenoloxidase homologs (ProPO1, ProPO2; Figure S1f in Additional file 1), which are homologous to D. melanogaster Diphenol oxidase A3 [Flybase: CG2952].
Nitric oxide synthase
Production of nitric oxide is mediated by the enzyme nitric oxide synthase. Nitric oxide is a highly unstable free radical gas that has been shown to be toxic to both parasites and pathogens. In insects, Nos is upregulated after both parasite and Gram-negative bacterial infection [86, 87]. Like other insects, pea aphids have one Nos homolog.
Heat shock proteins
Though called HSPs, these proteins are produced in response to a range of stresses in both eukaryotic and prokaryotic organisms . They serve as chaperones, facilitating protein folding and stabilization, and as proteases, mediating the degradation of damaged proteins. HSPs may also serve as signaling proteins during immune responses [18, 88]. In many insects, including aphids, HSPs have been shown to be upregulated after septic injury and microbial infection [31, 89–92]. We identified 15 HSPs of varying molecular weight in pea aphids (Figure S1g in Additional file 1).
Gluthione-S-tranferases comprise a diverse class of enzymes that detoxify stress-causing agents, including toxic oxygen free radical species. They are upregulated in some arthropods upon oxidative stress  and microbial challenge [89, 94]. Pea aphids have at least 18 genes encoding gluthione-S-tranferases and many other detoxification enzymes that likely play a role in stress responses . Ramsey et al.  identified many of the genes encoding detoxification enzymes in A. pisum and in Myzus persicae.
Alarm pheromone production
In response to predators, aphids release an alarm pheromone that causes neighboring aphids to become more mobile and to produce more winged than unwinged offspring [19, 96]. These winged offspring have the ability to disperse to enemy-free space. While many insects produce a suite of chemicals that constitute an alarm signal, the aphid alarm pheromone is dominated by a single compound, (E)-β farnesene . While the genes underlying alarm pheromone production have not been fully characterized, we have identified a Farnesyl diphosphate synthase (FPPS) and an Isoprenyl diphosphate synthase (IPPS), which may underlie alarm pheromone production .
We utilized real-time quantitative PCR to conduct a preliminary investigation of the expression of 23 recognition, signaling and response genes in aphids subjected to a number of infection and stress treatments (see Supplementary materials and Table S2 in Additional file 1). While future studies with more biological replicates will be necessary to fully survey gene regulation in the face of stress and infection, this initial survey indicates that aphids do express these genes under both control and infection/stress conditions (Tables S4 and S5 in Additional file 1). This suggests that these genes are functional even in the absence of many other missing immune-related genes.
One expression pattern seen in this initial survey is of particular note. Unlike other insect immune expression studies, we found no strong upregulation of antimicrobial peptides, which frequently exhibit ten-fold or greater upregulation in the face of infection. For example, while Altincicek et al.  observed 20-fold upregulation of Thaumatins in tribolium beetles after stabbing with lipopolysaccaride endotoxin derived from Escherichia coli, we saw modest upregulation (approximately 2-fold) of only one Thaumatin (Thm2) after stabbing aphids (Table S5 in Additional file 1). Furthermore, despite the fact that they are known to suppress fungal germination in beetles, the Thaumatin homologs were not upregulated after fungal infection at the time point included in this study, and were only approximately two-fold upregulated at two additional time points and in a follow-up fungal infection experiment (data not shown) . The role of thaumatins in fighting microbial infections, however, should not be discounted, as they may function in the absence of significant upregulation (that is, they may be constitutively expressed).
Exploration of ESTs from infected and uninfected aphids
In the first of two EST-based experiments, we compared a cDNA library synthesized from the guts of A. pisum that had been fed a Gram-negative pathogen, Dickeya dadantii, to a cDNA library synthesized from uninfected guts. Strikingly, no standard immune-related genes, such as antimicrobial peptides, were identified in the infected sample. The main functional classes differentially expressed were the 'biopolymer metabolism' class, many members of which were down-regulated in infected guts, and 'transport' or 'establishment of localization' classes, whose genes were upregulated in infected guts (Table S6 in Additional file 1). The 'immune response' class, in contrast, was only represented by five genes. Four of these five genes were in the uninfected library, while only one, encoding a leucyl-aminopeptidase, was identified from the infected library; the immune function of leucyl-aminopeptidases is not well understood. Moreover, the 'response to stress/external stimulus/biotic stimulus' classes were not overrepresented in the infected gut library.
In a separate experiment, to further identify aphid immune-relevant genes, we utilized SSH to compare cDNA from E. coli-infected aphids and cDNA from unchallenged aphids. To obtain genes expressed at different phases of the immune response, three RNA samples were extracted 3, 6 and 12 hours after E. coli infection and mixed prior to cDNA synthesis.
Among the 480 ESTs that were sequenced from the subtracted library [GenBank: GD185911 to GD186390], we found some genes with similarity to proteases and protease inhibitors but few other immune-related proteins. Interestingly, SSH-based EST analysis failed to identify any PRRs, such as PGRPs or GNBPs, or any antimicrobial peptides (Table S7 in Additional file 1). It is noteworthy that this aphid experiment was conducted in parallel to a similar Sitophilus weevil experiment, where many immune-related genes (more than 18% of ESTs) were identified, including antibacterial peptides and PRRs . This suggests that the paucity of immune genes identified in A. pisum is not a technical issue but may be a specific feature of aphids . In addition, dot blot analysis demonstrated that only a few genes (less than 5%) were differentially expressed between E. coli-stabbed and unstabbed aphids. These findings indicate that, in contrast to other insects, either aphids respond only weakly to challenge with E. coli or aphid genes and pathways directed against these bacteria are expressed only constitutively.
High performance liquid chromatography
HPLC peptide analyses targeting production of small peptides (for example, antimicrobial peptides) were run on hemolymph samples from pea aphids challenged by three microorganisms: E. coli (Gram-negative bacteria), Micrococcus luteus (Gram-positive bacteria) and Aspergillus fumigatus (fungi). Profiles were compared between control, infected and sterile-stabbed aphids at 6, 12 and 18 hours after challenge. When identified, the production of small peptides was maximal at 18 hours. In E. coli-treated samples, no upregulation could be identified (Figure 4a), in M. luteus-treated samples, there was modest upregulation (data not shown), and in A. fumigatus-treated samples, there was a significant response, though few peaks (Figure 4b). In contrast, a response profile to E. coli from another obligate symbiotic insect (the weevil, Sitophilus oryzae) exhibited at least five well-distinguishable upregulated peaks (Figure 4c). Response being restricted to Gram-positive bacteria and fungi is consistent with previous identification of megourin, an antimicrobial peptide in the aphid Megoura viciae, which appears to have activity against Gram-positive bacteria and fungi, but not against Gram-negative bacteria (P Bulet, unpublished). Because so few distinguishable peaks were present in the aphid samples, we did not choose to identify the associated products, but overall the presence of few inducible peptides suggests a peculiar scarcity of antimicrobial peptides in aphids.
Aphids are one of only a few genomic models for hemimetabolous insects, yet until recently, virtually nothing was known about aphid immune and stress response systems. Here, by coupling gene annotation with functional assays, we see evidence that aphids have some defense systems common to other arthropods (for example, the Toll and JAK/STAT signaling pathways, HSPs, ProPO). Surprisingly, however, several of the genes thought central to arthropod innate immunity are missing in aphids (for example, PGRPs, the IMD signaling pathway, defensins, c-type lysozymes). This calls into question the generality of the current model of insect immunity, and it remains to be determined how aphids protect themselves from the diverse pathogens and parasites that they face.
The fact that we cannot find aphid homologs to many insect immune genes could be a consequence of the large evolutionary distance between aphids and the taxa (in most cases, flies, mosquitoes and bees) from which these genes are known (that is, the split between the ancestors of aphids and these taxa occurred approximately 350 million years ago ), making it challenging to find divergent genes via homology-based searches, even when using highly sensitive methods as done here. Though we cannot preclude this possibility in all cases, in some cases, similar homology-based methods are able to recover homologs in even more distantly related taxa. For example, querying genome databases with Drosophila genes via BLAST recovers putative homologs of PGRPs and defensins in P. humanus (human body louse) and in Ixodes scapularis (deer tick) (Figure 2). The divergence time between Drosophila and these taxa is equal to or greater than that between Drosophila and aphids. Moreover, for some cases, we could identify genomic regions similar to functional genes in other species, but these regions contain large insertions or stop codons (for example, the putative antimicrobial peptide Megourin), indicating they are the result of pseudogenization.
One potential explanation for the lack of known immune-related genes in pea aphids is that aphids mount an alternative, but equal, immune response. Our functional analyses, as well as those of Altincicek et al. , found little evidence for an alternative response. In EST and HPLC analyses, few novel ESTs or peptide signals were recovered from immune-challenge aphids relative to their unchallenged controls. It should be noted, however, that these challenges were primarily limited to exposure to E. coli bacteria. When testing for expression of a few immune genes in response to a wider array of challenges, we do see some evidence of an aphid immune and stress response. Future expression studies, including large-scale transcriptional and proteomic studies, will extend this work and allow for more comprehensive characterization of the full complementation of aphid immune responses.
While we have focused mainly on the humoral component of the innate immune response, it is interesting to note that there is some evidence that the cellular component of pea aphids' innate immune response may also be different to that seen in other insects. While many insects encapsulate parasitoid wasp larvae, smothering them to death with hemocytes (insect immune cells), aphids appear not to have this layer of protection [101, 102]. Aphids, however, appear to recruit some hemocytes to parasitoid eggs, suggesting that cellular immunity may play an alternative, though possibly more limited, role . Better insights into the capacity of the aphid immune system will require further investigation of both the humoral and cellular components of aphid immunity.
The lack of genomic and molecular data regarding immune systems of aphid relatives makes it difficult to establish whether the pea aphid immune system is unique. There are, however, a number of aspects of aphid ecology that could facilitate ecological success without a strong immune defense. Altincicek et al.  proposed three hypotheses to explain the apparent lack of antimicrobial defenses. First, they suggested that contrary to Drosophila, whose natural environment consists of decaying fruit that is colonized by many microbes, aphids exploit phloem sap, which only occasionally contains bacteria and rarely contains entomopathogens. Thus, the risk of encountering pathogens while feeding is more limited. This assumption, however, is only partly true. While probing plants, aphids are capable of acquiring pathogenic bacteria from the surface of their host plants' leaves , and aphids become host to a diverse assemblage of bacteria and fungi under stressful conditions , some of which are pathogenic (NM Gerardo, unpublished data). Furthermore, Sitophilus weevils, which when challenged with E. coli significantly upregulate immune genes , spend their entire larval and nymph stages within sterile cereal grains, indicating that a sterile diet is not likely to explain the absence of antibacterial defenses in aphids.
Altincicek et al.  also suggest that aphids may invest in terminal reproduction in response to an immune challenge, rather than in a costly immune response. In their study, stabbed aphids produced significantly more offspring than untreated aphids within 24 hours of injury. Such an increase in reproduction upon challenge is not uncommon for invertebrates. Biomphalaria snails [105, 106], Acheta crickets , Daphnia waterfleas , and Drosophila flies  have all been shown to increase their investment in reproduction in response to infection. Yet, Drosophila still mount a complex immune response. Furthermore, aphids do not increase their reproductive effort in the face of all immune challenges: fungal infection reduces the number of offspring A. pisum produce within 24 hours of inoculation , and response to stabbing with bacteria seems to be specific to the aphid genotype and to the location of the stab (Barribeau, unpublished data). Therefore, though aphids have the capacity to reproduce many offspring prior to succumbing to some pathogens, it seems that immune competence would still provide increased fitness.
Even without increased reproduction following infection, the prolific reproductive capacity of aphids suggests these insects, in general, may invest most resources towards rapid, early onset reproduction rather than towards fewer, though better-protected offspring (aka, in terms of classical ecological theory, aphids may be r-selected rather k-selection organisms ). Recent theory of the evolution of immunity suggests that such organisms may specifically invest less in costly immune responses [112, 113]. Many characteristics of aphids, including their rapid generation time, short life span and small body size all fit a model of r-selection . Drosophila spp., however, also exhibit many of these characteristics and still invest in a strong defense repertoire.
The third hypothesis proposed by Altincicek et al.  concerning the evolution and maintenance of aphid defense relies on the presence of secondary symbionts that can be found extracellularly in aphids . A. pisum is protected against fungal pathogens by one of these secondary symbionts, Regiella insecticola , and also against the parasitoid wasp Aphidius ervi by another secondary symbiont, Hamiltonella defensa . Such symbiont-mediated host protection may explain why aphids have a reduced (or specialized) antimicrobial defense. This hypothesis seems plausible with regard to the cost of immune gene expression versus the benefit of protection by the secondary endosymbionts. However, it does not explain how the secondary endosymbionts (as Gram-negative bacteria), often present in aphid hemolymph, are themselves perceived and controlled by the aphid immune system. Thus, it is challenging to say whether the presence of secondary symbionts is a cause or a consequence of reduced antimicrobial activity.
Potentially, all of these forces could shape the evolution of aphid stress and immune responses. In order to test these hypotheses (for example, reproductive investment, symbiont-mediated host protection), we need more studies characterizing the global aphid response under more conditions, and in more aphid species. Potential insight from aphid relatives with different lifestyles (for example, those not associated with secondary symbionts, or those that live in soil or other microbe-rich habitats) may be particularly helpful. More broadly, as the pea aphid is the first published genome of a hemimetabolous insect, future analyses of the immune and stress related genes of more insects in this group will facilitate the reconstruction of the evolutionary history of innate immunity and other defenses.
Materials and methods
Bioinformatic screening of the pea aphid genome
Immune and stress gene candidates from other insects (for example, D. melanogaster, A. aegypti, A. gambiae, A. mellifera) were used to query the pea aphid genome. Most searches utilized the blastp search function to search for hits against the predicted A. pisum proteome . For some gene families and putative paralogs, protein sequences were aligned to sequences from other insects and outgroups using ClustalW . These alignments, as well as available EST and full length cDNA sequences, served to refine aphid gene models (exon/intron boundaries, and so on), and to facilitate phylogenetic analyses. In addition, a comprehensive database of all available EST sequences from the green peach aphid, Myzus persicae, was screened using tblastn to search for potential homologs to all immune and stress genes annotated in the pea aphid.
For genes that could not be found in the proteome, we also conducted a tblastn search against all contigs and unassembled reads. Then, a final, more sensitive profile-based search was performed for those immune defense proteins that produced no hits with BLAST searches. For this analysis, insect and other species protein sequences belonging to the family of interest were retrieved from NCBI and aligned with MUSCLE . A hidden Markov model for the alignment was built and calibrated using HMMER . This was used to perform a profile-based search (hmmsearch) against the six-frame translated sequences of the assembled pea aphid genome and the unassembled reads. Additionally, a similar search with PFAM profiles  was also performed for those families encoding PFAM domains in their sequences. Whenever a significant hit was found, the genomic region was analyzed to discard the possibility that it encoded a pseudogene (presence of stop codons, absence of relevant domains, and so on).
Phylogenetic analyses of selected protein families were performed using their corresponding maximum likelihood phylogenetic trees from the pea aphid phylome , deposited in PhylomeDB . When necessary, additional sequences were added to the original PhylomeDB alignment, realigned with MUSCLE and used to reconstruct a maximum likelihood phylogenetic tree, using the JTT (Jones-Taylor-Thornton) model as implemented in PhyML v2.4.4 , assuming a discrete gamma-distribution model with four rate categories and invariant sites, and estimating the gamma shape parameter and the fraction of invariant sites. Cladograms were edited using Dendrogram .
Exploration of ESTs from infected and uninfected aphids
In the first experiment, two EST libraries (one control, one infected) were generated by standard procedures using a SMART cDNA kit (Clontech, Mountain View, California, USA), starting from approximately 1,000 dissected A. pisum midguts for each library. The aphids were clonal, young, reproducing asexuals, which were either fed on control diet or infected by feeding on artificial diet with the Gram-negative aphid pathogen Dickeya dadantii at 106 bacteria per milliliter . Twenty-four hours after infection, control and treated aphids were dissected, and complete guts were transferred immediately into RNeasy solution (Qiagen Valencia, California, USA). ESTs were sequenced according to procedures in Sabater-Munoz et al. .
In another EST-based experiment utilizing SSH and dot-blot technology, we treated aphids (clone LL01) with rifampicin as described in Rahbé et al.  to reduce symbiont load. We challenged wingless fourth-instar aposymbiotic aphids by stabbing them with needles previously dipped into a pellet of overnight cultures of E. coli (TOP10, Invitrogen Carlsbad, California, USA), and then maintained them on fava plants. At 3, 6, and 12 hours post-treatment, we stored surviving aphids at -80°C. To identify genes that are differentially expressed in response to septic injury, we performed SSH using RNAs from immune challenged (3, 6 and 12 hours post-treatment) and untreated aposymbiotic aphids, using the SMART PCR cDNA Synthesis Kit and the PCR-Select cDNA Subtraction Kit (Clontech) according to the manufacturer's instructions and as described in Anselme et al. . After transformation by electroporation, we recovered approximately 1,500 colonies from LB agar plates. We plasmid extracted and sequenced 500 randomly picked colonies (NucleoSpin® Plasmid Kit, Macherey-Nagel, Düren, Germany) utilizing the sequencing center at the University of Valencia (Spain). We compared all sequences against UniProt using blastx. Immune-related gene sequences (Table S7 in Additional file 1) were then compared to the aphid genome using blastn.
To analyze the differential expression status of each EST, we conducted a dot-blot experiment. Briefly, we amplified 344 ESTs from the SSH library by colony PCR with nested PCR primers 1 and 2R from the PCR-Select cDNA Subtraction Kit. We then spotted 10 μl from each PCR product onto two different membranes (Amersham Hybond™-N, GE Healthcare Life Sciences, Piscataway, New Jersey, USA) using a Bio-Dot Microfiltration System (Biorad, Hercules, California, USA). We hybridized membranes with radiolabeled cDNA probes generated by reverse-transcription from RNA extracted from either aposymbiotic aphids stabbed with E. coli or unstabbed aposymbiotic aphids. We synthesized these probes using the Super Script™ First Strand Synthesis system (Invitrogen) for RT-PCR and [α-32P]dCTP, and purified them using Quick Spin Column (Roche Molecular Biochemicals, Indianapolis, Indiana, USA). After exposing blots for up to 24 hours to a Storm PhosphorImager imaging plate (GE Healthcare Life Sciences), we analyzed differential expression by comparison of band intensities between the two membranes. We did not, however, normalize the data, as we failed to see any signal from the Gapdh gene, though the same amount of each PCR product was loaded on both membranes.
Aphids were challenged by abdominal puncture with triple-0 needles dipped in a solution of Gram-negative bacteria (E. coli strain Top10), Gram-positive bacteria (M. luteus) or fungal spores (A. fumigatus). For each microbial treatment, five hemolymph samples from 50 aphids each were collected at four times points (t = 0, 6, 12 and 18 hours).
Hemolymph was flash-extracted by centrifuging (1 minute, 10,000 g, 4°C) live aphids through a 1 ml pipette tip and directly into 40 μl 0.1% trifluoractetic acid contaning 10 μl of saturated phenylthiourea (PTU) for phenoloxidase inhibition. Resulting samples were highly similar to pure hemolymph samples obtained by leg bleeding (>95% band identity by silver-stained SDS-PAGE).
After initial collection, tips were removed and the samples were centrifuged for 5 minutes at 15,000 g. Following addition of 70 μl trifluoractetic acid 0.1%, the supernatant sat for 1 hour at 4°C to allow for protein precipitation prior to a final 10-minute centrifugation at 15,000 g to recover peptides. Samples were evaporated and stored at -20°C until use in HPLC. Chromatography was performed on standard peptide C18-300Å reverse phase columns using water acetonitrile gradients . For retention time standardization, PTU served as an internal standard, and samples were analyzed by area-normalization to unchallenged sample peaks (retention time = 14 minutes, preceding PTU).
expressed sequence tag
Gram-negative binding protein
high performance liquid chromatography
heat shock protein
Janus kinase/Signal transducers and activators of transcription
c-Jun N-terminal kinase
petidoglycan receptor protein
pathogen recognition receptor
suppression subtractive hybridization
Snyder WE, Ives AR: Interactions between specialist and generalist natural enemies: Parasitoids, predators, and pea aphid biocontrol. Ecology. 2003, 84: 91-107. 10.1890/0012-9658(2003)084[0091:IBSAGN]2.0.CO;2.
Hufbauer RA: Aphid population dynamics: does resistance to parasitism influence population size?. Ecol Entomol. 2002, 27: 25-32. 10.1046/j.1365-2311.2002.0379a.x.
Hatano E, Kunert G, Bartram S, Boland W, Gershenzon J, Weisser WW: Do aphid colonies amplify their emission of alarm pheromone?. J Chem Ecol. 2008, 34: 1149-1152. 10.1007/s10886-008-9527-y.
Tarpy DR: Genetic diversity within honeybee colonies prevents severe infections and promotes colony growth. Proc Biol Sci. 2003, 270: 99-103. 10.1098/rspb.2002.2199.
Ha EM, Oh CT, Ryu JH, Bae YS, Kang SW, Jang IH, Brey PT, Lee WJ: An antioxidant system required for host protection against gut infection in/Drosophila. Dev Cell. 2005, 8: 125-132. 10.1016/j.devcel.2004.11.007.
Francke DL, Harmon JP, Harvey CT, Ives AR: Pea aphid dropping behavior diminishes foraging efficiency of a predatory ladybeetle. Entomologia Experimentalis Et Applicata. 2008, 127: 118-124. 10.1111/j.1570-7458.2008.00678.x.
Gagneux S, DeRiemer K, Van T, Kato-Maeda M, de Jong BC, Narayanan S, Nicol M, Niemann S, Kremer K, Gutierrez MC, Hilty M, Hopewell PC, Small PM: Variable host-pathogen compatibility in Mycobacterium tuberculosis. Proc Natl Acad Sci USA. 2006, 103: 2869-2873. 10.1073/pnas.0511240103.
Govind S: Innate immunity in Drosophila: pathogens and pathways. Insect Sci. 2008, 15: 29-43.
Strand MR: The insect cellular immune response. Insect Sci. 2008, 15: 1-14.
Schulenburg H, Boehnisch C, Michiels NK: How do invertebrates generate a highly specific innate immune response?. Mol Immunol. 2007, 44: 3338-3344. 10.1016/j.molimm.2007.02.019.
Boutros M, Agaisse H, Perrimon N: Sequential activation of signaling pathways during innate immune responses in Drosophila. Dev Cell. 2002, 3: 711-722. 10.1016/S1534-5807(02)00325-8.
Dionne MS, Schneider DS: Models of infectious diseases in the fruit fly Drosophila melanogaster. Dis Model Mech. 2008, 1: 43-49. 10.1242/dmm.000307.
Zou Z, Evans JD, Lu ZQ, Zhao PC, Williams M, Sumathipala N, Hetru C, Hultmark D, Jiang HB: Comparative genomic analysis of the Tribolium immune system. Genome Biol. 2007, 8: R177-10.1186/gb-2007-8-8-r177.
Evans JD, Aronstein K, Chen YP, Hetru C, Imler JL, Jiang H, Kanost M, Thompson GJ, Zou Z, Hultmark D: Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Mol Biol. 2006, 15: 645-656. 10.1111/j.1365-2583.2006.00682.x.
Waterhouse RM, Kriventseva EV, Meister S, Xi ZY, Alvarez KS, Bartholomay LC, Barillas-Mury C, Bian GW, Blandin S, Christensen BM, Dong YM, Jiang HB, Kanost MR, Koutsos AC, Levashina EA, Li JY, Ligoxygakis P, MacCallum RM, Mayhew GF, Mendes A, Michel K, Osta MA, Paskewitz S, Shin SW, Vlachou D, Wang LH, Wei WQ, Zheng LB, Zou Z, Severson DW, et al: Evolutionary dynamics of immune-related genes and pathways in disease-vector mosquitoes. Science. 2007, 316: 1738-1743. 10.1126/science.1139862.
Christophides GK, Zdobnov E, Barillas-Mury C, Birney E, Blandin S, Blass C, Brey PT, Collins FH, Danielli A, Dimopoulos G, Hetru C, Hoa N, Hoffmann JA, Kanzok SM, Letunic I, Levashina EA, Loukeris TG, Lycett G, Meister S, Michel K, Muller HM, Osta MA, Paskewitz SM, Reichhart JM, Rzhetsky A, Troxler L, Vernick KD, Vlachou D, Volz J, von Mering C, et al: Immunity-related genes and gene families in Anopheles gambiae. Science. 2002, 298: 159-165. 10.1126/science.1077136.
Sackton TB, Lazzaro BP, Schlenke TA, Evans JD, Hultmark D, Clark AG: Dynamic evolution of the innate immune system in Drosophila. Nat Genet. 2007, 39: 1461-1468. 10.1038/ng.2007.60.
Pockley AG: Heat shock proteins as regulators of the immune response. Lancet. 2003, 362: 469-476. 10.1016/S0140-6736(03)14075-5.
Kunert G, Otto S, Rose USR, Gershenzon J, Weisser WW: Alarm pheromone mediates production of winged dispersal morphs in aphids. Ecol Lett. 2005, 8: 596-603. 10.1111/j.1461-0248.2005.00754.x.
Moran NA, Russell JA, Koga R, Fukatsu T: Evolutionary relationships of three new species of Enterobacteriaceae living as symbionts of aphids and other insects. Appl Environ Microbiol. 2005, 71: 3302-3310. 10.1128/AEM.71.6.3302-3310.2005.
Moran NA, Telang A: Bacteriocyte-associated symbionts of insects - a variety of insect groups harbor ancient prokaryotic endosymbionts. Bioscience. 1998, 48: 295-304. 10.2307/1313356.
Sandstrom JP, Russell JA, White JP, Moran NA: Independent origins and horizontal transfer of bacterial symbionts of aphids. Mol Ecol. 2001, 10: 217-228. 10.1046/j.1365-294X.2001.01189.x.
Tsuchida T, Koga R, Fukatsu T: Host plant specialization governed by facultative symbiont. Science. 2004, 303: 1989-1989. 10.1126/science.1094611.
Russell JA, Moran NA: Costs and benefits of symbiont infection in aphids: variation among symbionts and across temperatures. Proc Biol Sci. 2006, 273: 603-610. 10.1098/rspb.2005.3348.
Montllor CB, Maxmen A, Purcell AH: Facultative bacterial endosymbionts benefit pea aphids Acyrthosiphon pisum under heat stress. Ecol Entomol. 2002, 27: 189-195. 10.1046/j.1365-2311.2002.00393.x.
Chen DQ, Montllor CB, Purcell AH: Fitness effects of two facultative endosymbiotic bacteria on the pea aphid, Acyrthosiphon pisum, and the blue alfalfa aphid, A. kondoi. Entomologia Experimentalis Et Applicata. 2000, 95: 315-323. 10.1023/A:1004083324807.
Oliver KM, Moran NA, Hunter MS: Variation in resistance to parasitism in aphids is due to symbionts not host genotype. Proc Natl Acad Sci USA. 2005, 102: 12795-12800. 10.1073/pnas.0506131102.
Oliver KM, Russell JA, Moran NA, Hunter MS: Facultative bacterial symbionts in aphids confer resistance to parasitic wasps. Proc Natl Acad Sci USA. 2003, 100: 1803-1807. 10.1073/pnas.0335320100.
Scarborough CL, Ferrari J, Godfray HCJ: Aphid protected from pathogen by endosymbiont. Science. 2005, 310: 1781-1781. 10.1126/science.1120180.
Haine ER: Symbiont-mediated protection. Proc Biol Sci. 2008, 275: 353-361. 10.1098/rspb.2007.1211.
Altincicek B, Gross J, Vilcinskas A: Wounding-mediated gene expression and accelerated viviparous reproduction of the pea aphid Acyrthosiphon pisum. Insect Mol Biol. 2008, 17: 711-716. 10.1111/j.1365-2583.2008.00835.x.
Altincicek B, Knorr E, Vilcinskas A: Beetle immunity: Identification of immune-inducible genes from the model insect Tribolium castaneum. Dev Comp Immunol. 2008, 32: 585-595. 10.1016/j.dci.2007.09.005.
Altincicek B, Vilcinskas A: Analysis of the immune-inducible transcriptome from microbial stress resistant, rat-tailed maggots of the drone fly Eristalis tenax. BMC Genomics. 2007, 8: 326-10.1186/1471-2164-8-326.
Altincicek B, Vilcinskas A: Identification of immune-related genes from an apterygote insect, the firebrat Thermobia domestica. Insect Biochem Mol Biol. 2007, 37: 726-731. 10.1016/j.ibmb.2007.03.012.
Anselme C, Perez-Brocal V, Vallier A, Vincent-Monegat C, Charif D, Latorre A, Moya A, Heddi A: Identification of the weevil immune genes and their expression in the bacteriome tissue. BMC Biol. 2008, 6: 43-10.1186/1741-7007-6-43.
The International Aphid Genomics Consortium: Genome sequence of the pea aphid Acyrthosiphon pisum. PloS Biol. 2010, 8: e1000313-10.1371/journal.pbio.1000313.
Kaneko T, Silverman N: Bacterial recognition and signalling by the Drosophila IMD pathway. Cell Microbiol. 2005, 7: 461-469. 10.1111/j.1462-5822.2005.00504.x.
Steiner H: Peptidoglycan recognition proteins: on and off switches for innate immunity. Immunol Rev. 2004, 198: 83-96. 10.1111/j.0105-2896.2004.0120.x.
Werner T, Liu G, Kang D, Ekengren S, Steiner H, Hultmark D: A family of peptidoglycan recognition proteins in the fruit fly Drosophila melanogaster. Proc Natl Acad Sci USA. 2000, 97: 13772-13777. 10.1073/pnas.97.25.13772.
McTaggert SJ, Conlon C, Colbourne JK, Blaxter ML, Little TJ: The components of the Daphnia pulex immune system as revealed by complete genome sequencing. BMC Genomics. 2009, 10: 175-10.1186/1471-2164-10-175.
Lemaitre B, Hoffmann J: The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007, 25: 697-743. 10.1146/annurev.immunol.25.022106.141615.
Goftar M, Gobert V, Matskevich AA, Reichhart JM, Wang CS, Buft TM, BeIvin M, Hoffmann JA, Ferrandon D: Dual detection of fungal infections in Drosophila via recognition of glucans and sensing of virulence factors. Cell. 2006, 127: 1425-1437. 10.1016/j.cell.2006.10.046.
Tanji T, Ohashi-Kobayashi A, Natori S: Participation of a galactose-specific C-type lectin in Drosophila immunity. Biochem J. 2006, 396: 127-138. 10.1042/BJ20051921.
Ao JQ, Ling EJ, Yu XQ: Drosophila C-type lectins enhance cellular encapsulation. Mol Immunol. 2007, 44: 2541-2548. 10.1016/j.molimm.2006.12.024.
Pace KE, Baum LG: Insect galectins: Roles in immunity and development. Glycoconj J. 2004, 19: 607-614. 10.1023/B:GLYC.0000014092.86763.2f.
Dimopoulos G, Seeley D, Wolf A, Kafatos FC: Malaria infection of the mosquito Anopheles gambiae activates immune-responsive genes during critical transition stages of the parasite life cycle. EMBO J. 1998, 17: 6115-6123. 10.1093/emboj/17.21.6115.
Dimopoulos G, Richman A, dellaTorre A, Kafatos FC, Louis C: Identification and characterization of differentially expressed cDNAs of the vector mosquito, Anopheles gambiae. Proc Natl Acad Sci USA. 1996, 93: 13066-13071. 10.1073/pnas.93.23.13066.
Ramet M, Pearson A, Manfruelli P, Li XH, Koziel H, Gobel V, Chung E, Krieger M, Ezekowitz RAB: Drosophila scavenger receptor Cl is a pattern recognition receptor for bacteria. Immunity. 2001, 15: 1027-1038. 10.1016/S1074-7613(01)00249-7.
Lazzaro BP: Elevated polymorphism and divergence in the class C scavenger receptors of Drosophila melanogaster and D. simulans. Genetics. 2005, 169: 2023-2034. 10.1534/genetics.104.034249.
Ju JS, Cho MH, Brade L, Kim JH, Park JW, Ha NC, Soderhall I, Soderhall K, Brade H, Lee BL: A novel 40-kDa protein containing six repeats of an epidermal growth factor-like domain functions as a pattern recognition protein for lipopolysaccharide. J Immunol. 2006, 177: 1838-1845.
Kurucz E, Markus R, Zsamboki J, Folkl-Medzihradszky K, Darula Z, Vilmos P, Udvardy A, Krausz I, Lukacsovich T, Gateff E, Zettervall CJ, Hultmark D, Ando I: Nimrod, a putative phagocytosis receptor with EGF repeats in Drosophila plasmatocytes. Curr Biol. 2007, 17: 649-654. 10.1016/j.cub.2007.02.041.
Somogyi K, Sipos B, Penzes Z, Kurucz E, Zsamboki J, Hultmark D, Ando I: Evolution of genes and repeats in the Nimrod superfamily. Mol Biol Evol. 2008, 25: 2337-2347. 10.1093/molbev/msn180.
Dong YM, Taylor HE, Dimopoulos G: AgDscam, a hypervariable immunoglobulin domain-containing receptor of the Anopheles gambiae innate immune system. PLoS Biol. 2006, 4: e229-10.1371/journal.pbio.0040229.
Brites D, McTaggart S, Morris K, Anderson J, Thomas K, Colson I, Fabbro T, Little TJ, Ebert D, Du Pasquier L: The Dscam homologue of the crustacean Daphnia is diversified by alternative splicing like in insects. Mol Biol Evol. 2008, 25: 1429-1439. 10.1093/molbev/msn087.
Watson FL, Puttmann-Holgado R, Thomas F, Lamar DL, Hughes M, Kondo M, Rebel VI, Schmucker D: Extensive diversity of Ig-superfamily proteins in the immune system of insects. Science. 2005, 309: 1874-1878. 10.1126/science.1116887.
Schlenke TA, Morales J, Govind S, Clark AG: Contrasting infection strategies in generalist and specialist wasp parasitoids of Drosophila melanogaster. PloS Pathogens. 2007, 3: 1486-1501. 10.1371/journal.ppat.0030158.
Kanost MR: Serine proteinase inhibitors in arthropod immunity. Dev Comp Immunol. 1999, 23: 291-301. 10.1016/S0145-305X(99)00012-9.
Krem MM, Di Cera E: Evolution of enzyme cascades from embryonic development to blood coagulation. Trends Biochem Sci. 2002, 27: 67-74. 10.1016/S0968-0004(01)02007-2.
Rawlings ND, Barrett AJ: Evolutionary families of peptidases. Biochem J. 1993, 290: 205-218.
Imler JL, Hoffmann JA: Toll receptors in Drosophila: a family of molecules regulating development and immunity. Toll-Like Receptor Family Members and Their Ligands. Edited by: Beutler B, Wagner H. 2002, Berlin: Springer, 270: 63-79. [Current Topics in Microbiology and Immunology, volume 270]
Leulier F, Lemaitre B: Toll-like receptors - taking an evolutionary approach. Nat Rev Genet. 2008, 9: 165-178. 10.1038/nrg2303.
Tauszig S, Jouanguy E, Hoffmann JA, Imler JL: Toll-related receptors and the control of antimicrobial peptide expression in Drosophila. Proc Natl Acad Sci USA. 2000, 97: 10520-10525. 10.1073/pnas.180130797.
Agaisse H, Perrimon N: The roles of JAK/STAT signaling in Drosophila immune responses. Immunol Rev. 2004, 198: 72-82. 10.1111/j.0105-2896.2004.0133.x.
Lemaitre B, Kromermetzger E, Michaut L, Nicolas E, Meister M, Georgel P, Reichhart JM, Hoffmann JA: A recessive mutation, immume-deficiency (IMD), defines 2 distinct control pathways in the Drosophila host-defense. Proc Natl Acad Sci USA. 1995, 92: 9465-9469. 10.1073/pnas.92.21.9465.
Delaney JR, Stoven S, Uvell H, Anderson KV, Engstrom Y, Mlodzik M: Cooperative control of Drosophila immune responses by the JNK and NF-kappa B signaling pathways. EMBO J. 2006, 25: 3068-3077. 10.1038/sj.emboj.7601182.
Bidla G, Dushay MS, Theopold U: Crystal cell rupture after injury in Drosophila requires the JNK pathway, small GTPases and the TNF homolog Eiger. J Cell Sci. 2007, 120: 1209-1215. 10.1242/jcs.03420.
Igaki T, Kanda H, Yamamoto-Goto Y, Kanuka H, Kuranaga E, Aigaki T, Miura M: Eiger, a TNF superfamily ligand that triggers the Drosophila JNK pathway. EMBO J. 2002, 21: 3009-3018. 10.1093/emboj/cdf306.
Bulet P, Stocklin R: Insect antimicrobial peptides: Structures, properties and gene regulation. Protein Pept Lett. 2005, 12: 3-11. 10.2174/0929866053406011.
Bulet P, Stocklin R, Menin L: Anti-microbial peptides: from invertebrates to vertebrates. Immunol Rev. 2004, 198: 169-184. 10.1111/j.0105-2896.2004.0124.x.
Lemaitre B, Reichhart JM, Hoffmann JA: Drosophila host defense: Differential induction of antimicrobial peptide genes after infection by various classes of microorganisms. Proc Natl Acad Sci USA. 1997, 94: 14614-14619. 10.1073/pnas.94.26.14614.
Hoffmann JA, Reichhart JM: Drosophila innate immunity: an evolutionary perspective. Nat Immunol. 2002, 3: 121-126. 10.1038/ni0202-121.
Shatters RG, Boykin LM, Lapointe SL, Hunter WB, Weathersbee AA: Phylogenetic and structural relationships of the PR5 gene family reveal an ancient multigene family conserved in plants and select animal taxa. J Mol Evol. 2006, 63: 12-29. 10.1007/s00239-005-0053-z.
Jolles P: Lysozymes: Model Enzymes in Biochemistry and Biology. 1996, Basel, Switzerland: Birkhauser
Regel R, Matioli SR, Terra WR: Molecular adaptation of Drosophila melanogaster lysozymes to a digestive function. Insect Biochem Mol Biol. 1998, 28: 309-319. 10.1016/S0965-1748(97)00108-2.
Li B, Calvo E, Marinotti O, James AA, Paskewitz SM: Characterization of the c-type lysozyme gene family in Anopheles gambiae. Gene. 2005, 360: 131-139. 10.1016/j.gene.2005.07.001.
Paskewitz SM, Li B, Kajla MK: Cloning and molecular characterization of two invertebrate-type lysozymes from Anopheles gambiae. Insect Mol Biol. 2008, 17: 217-225. 10.1111/j.1365-2583.2008.00797.x.
Nakabachi A, Shigenobu S, Sakazume N, Shiraki T, Hayashizaki Y, Carninci P, Ishikawa H, Kudo T, Fukatsu T: Transcriptome analysis of the aphid bacteriocyte, the symbiotic host cell that harbors an endocellular mutualistic bacterium, Buchnera. Proc Natl Acad Sci USA. 2005, 102: 5477-5482. 10.1073/pnas.0409034102.
Daimon T, Hamada K, Mita K, Okano K, Suzuki MG, Kobayashi M, Shimada T: A Bombyx mori gene, BmChi-h, encodes a protein homologous to bacterial and baculovirus chitinases. Insect Biochem Mol Biol. 2003, 33: 749-759. 10.1016/S0965-1748(03)00084-5.
Nakabachi A, Shigenobu S, Miyagishima S: Chitinase-like proteins encoded in the genome of the pea aphid, Acyrthosiphon pisum. Insect Mol Biol. 2010, 19: 175-185. 10.1111/j.1365-2583.2009.00985.x.
Lagueux M, Perrodou E, Levashina EA, Capovilla M, Hoffmann JA: Constitutive expression of a complement-like protein in Toll and JAK gain-of-function mutants of Drosophila. Proc Natl Acad Sci USA. 2000, 97: 11427-11432. 10.1073/pnas.97.21.11427.
Agaisse H, Petersen UM, Boutros M, Mathey-Prevot B, Perrimon N: Signaling role of hemocytes in Drosophila JAK/STAT-dependent response to septic injury. Dev Cell. 2003, 5: 441-450. 10.1016/S1534-5807(03)00244-2.
Ekengren S, Hultmark D: A family of Turandot-related genes in the humoral stress response of Drosophila. Biochem Biophys Res Comm. 2001, 284: 998-1003. 10.1006/bbrc.2001.5067.
Ekengren S, Tryselius Y, Dushay MS, Liu G, Steiner H, Hultmark D: A humoral stress response in Drosophila. Curr Biol. 2001, 11: 714-718. 10.1016/S0960-9822(01)00203-2.
Nappi AJ, Christensen BM: Melanogenesis and associated cytotoxic reactions: applications to insect innate immunity. Insect Biochem Mol Biol. 2005, 35: 443-459. 10.1016/j.ibmb.2005.01.014.
Soderhall K, Cerenius L: Role of the prophenoloxidase-activating system in invertebrate immunity. Curr Opin Immunol. 1998, 10: 23-28. 10.1016/S0952-7915(98)80026-5.
Foley E, O'Farrell PH: Nitric oxide contributes to induction of innate immune responses to Gram-negative bacteria in Drosophila. Genes Dev. 2003, 17: 115-125. 10.1101/gad.1018503.
Rivero A: Nitric oxide: an antiparasitic molecule of invertebrates. Trends Parasitol. 2006, 22: 219-225. 10.1016/j.pt.2006.02.014.
Henderson B, Allan E, CoateS ARM: Stress wars: the direct role of host and bacterial molecular chaperones in bacterial infection. Infect Immun. 2006, 74: 3693-3706. 10.1128/IAI.01882-05.
Aguilar R, Jedlicka AE, Mintz M, Mahairaki V, Scott AL, Dimopoulos G: Global gene expression analysis of Anopheles gambiae responses to microbial challenge. Insect Biochem Mol Biol. 2005, 35: 709-719. 10.1016/j.ibmb.2005.02.019.
Dimopoulos G, Christophides GK, Meister S, Schultz J, White KP, Barillas-Mury C, Kafatos FC: Genome expression analysis of Anopheles gambiae: Responses to injury, bacterial challenge, and malaria infection. Proc Natl Acad Sci USA. 2002, 99: 8814-8819. 10.1073/pnas.092274999.
Guedes SD, Vitorino R, Domingues R, Tomer K, Correia AJF, Amado F, Domingues P: Proteomics of immune-challenged Drosophila melanogaster larvae hemolymph. Biochem Biophys Res Comm. 2005, 328: 106-115. 10.1016/j.bbrc.2004.12.135.
Ursic-Bedoya RJ, Lowenberger CA: Rhodnius prolixus: Identification of immune-related genes up-regulated in response to pathogens and parasites using suppressive subtractive hybridization. Dev Comp Immunol. 2007, 31: 109-120. 10.1016/j.dci.2006.05.008.
Girardot F, Monnier V, Tricoire H: Genome wide analysis of common and specific stress responses in adult Drosophila melanogaster. BMC Genomics. 2004, 5: 74-10.1186/1471-2164-5-74.
Rudenko N, Golovchenko M, Edwards MJ, Grubhoffer L: Differential expression of Ixodes ricinus tick genes induced by blood feeding or Borrelia burgdorferi infection. J Med Entomol. 2005, 42: 36-41. 10.1603/0022-2585(2005)042[0036:DEOIRT]2.0.CO;2.
Ramsey JS, Rider DS, Walsh T, de Vos M, Gordon K, Ponnala L, Roe BA, Jander G: Comparative analysis of detoxification enzymes in Acrythrosiphon pisum and Myzus persicae. Insect Mol Biol. 2010, 19: 155-164. 10.1111/j.1365-2583.2009.00973.x.
Kislow CJ, Edwards LJ: Repellent odor in aphids. Nature. 1972, 235: 108-109. 10.1038/235108a0.
Xiangyu JG, Zhang F, Fang YL, Kan W, Zhang GX, Zhang ZN: Behavioural response of aphids to the alarm pheromone component (E)-beta-farnesene in the field. Physiol Entomol. 2002, 27: 307-311. 10.1046/j.1365-3032.2002.00302.x.
Lewis MJ, Prosser IM, Mohib A, Field LM: Cloning and characterisation of a prenyltransferase from the aphid Myzus persicae with potential involvement in alarm pheromone biosynthesis. Insect Mol Biol. 2008, 17: 437-443. 10.1111/j.1365-2583.2008.00815.x.
Grenier AM, Duport G, Pages S, Condemine G, Rahbe Y: The phytopathogen Dickeya dadantii (Erwinia chrysanthemi 3937) is a pathogen of the pea aphid. Appl Environ Microbiol. 2006, 72: 1956-1965. 10.1128/AEM.72.3.1956-1965.2006.
Gaunt MW, Miles MA: An insect molecular clock dates the origin of the insects and accords with palaeontological and biogeographic landmarks. Mol Biol Evol. 2002, 19: 748-761.
Bensadia F, Boudreault S, Guay JF, Michaud D, Cloutier C: Aphid clonal resistance to a parasitoid fails under heat stress. J Insect Physiol. 2006, 52: 146-157. 10.1016/j.jinsphys.2005.09.011.
Carver M, Sullivan DJ: Encapsulative defence reactions of aphids (Hemiptera: Aphididae) to insect parasitoids (Hymenoptera: Aphidiidae and Aphelinidae). Ecology and Effectiveness of Aphidophaga. Edited by: Niemczyk E, Dixon AFG. 1988, The Hague: SPB Academic Publishing, 299-303.
Stavrinides J, McCloskey JK, Ochman H: The pea aphid as both host and vector for the phytopathogenic bacterium, Pseudomonas syringae. Appl Environ Microbiol. 2009, 75: 2230-2235. 10.1128/AEM.02860-08.
Nakabachi A, Ishikawa H, Kudo T: Extraordinary proliferation of microorganisms in aposymbiotic pea aphids, Acyrthosiphon pisum. J Invertebr Pathol. 2003, 82: 152-161. 10.1016/S0022-2011(03)00020-X.
Minchella DJ, Loverde PT: A cost of increased early reproductive effort in the snail Biomphalaria glabarata. Am Nat. 1981, 118: 876-881. 10.1086/283879.
Minchella DJ, Leathers BK, Brown KM, McNair JN: Host and parasite counteradaptations - an example from a fresh-water snail. Am Nat. 1985, 126: 843-854. 10.1086/284456.
Adamo SA: Evidence for adaptive changes in egg laying in crickets exposed to bacteria and parasites. Anim Behav. 1999, 57: 117-124. 10.1006/anbe.1998.0999.
Chadwick W, Little TJ: A parasite-mediated life-history shift in Daphnia magna. Proc Biol Sci. 2005, 272: 505-509. 10.1098/rspb.2004.2959.
Polak M, Starmer WT: Parasite-induced risk of mortality elevates reproductive effort in male Drosophila. Proc Biol Sci. 1998, 265: 2197-2201. 10.1098/rspb.1998.0559.
Baverstock J, Roy HE, Clark SJ, Alderson PG, Pell JK: Effect of fungal infection on the reproductive potential of aphids and their progeny. J Invertebr Pathol. 2006, 91: 136-139. 10.1016/j.jip.2005.11.005.
MacArthur RH, Wilson EO: The Theory of Island Biogeography. 1967, Princeton, NJ: Princeton University Press
Zuk M, Stoehr AM: Immune defense and host life history. Am Nat. 2002, 160: S9-S22. 10.1086/342131.
Miller MR, White A, Boots M: Host life span and the evolution of resistance characteristics. Evolution. 2007, 61: 2-14. 10.1111/j.1558-5646.2007.00001.x.
Pianka ER: On r- and K-selection. Am Nat. 1970, 104: 592-597. 10.1086/282697.
Fukatsu T, Nikoh N, Kawai R, Koga R: The secondary endosymbiotic bacterium of the pea aphid Acyrthosiphon pisum (Insecta: Homoptera). Appl Environ Microbiol. 2000, 66: 2748-2758. 10.1128/AEM.66.7.2748-2758.2000.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic Local Alignment Search Tool. J Mol Biol. 1990, 215: 403-410.
Thompson JD, Higgins DG, Gibson TJ: Clustal-W - Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Eddy SR: Profile hidden Markov models. Bioinformatics. 1998, 14: 755-763. 10.1093/bioinformatics/14.9.755.
Bateman A, Birney E, Cerruti L, Durbin R, Etwiller L, Eddy SR, Griffiths-Jones S, Howe KL, Marshall M, Sonnhammer ELL: The Pfam Protein Families Database. Nucleic Acids Res. 2002, 30: 276-280. 10.1093/nar/30.1.276.
Huerta-Cepas J, Bueno A, Dopazo JQ, Gabaldon T: PhylomeDB: a database for genome-wide collections of gene phylogenies. Nucleic Acids Res. 2008, 36: D491-D496. 10.1093/nar/gkm899.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Huson DH, Richter DC, Rausch C, Dezulian T, Franz M, Rupp R: Dendroscope: An interactive viewer for large phylogenetic trees. BMC Bioinformatics. 2007, 8: 460-10.1186/1471-2105-8-460.
Sabater-Munoz B, Legeai F, Rispe C, Bonhomme J, Dearden P, Dossat C, Duclert A, Gauthier JP, Ducray DG, Hunter W, Dang P, Kambhampati S, Martinez-Torres D, Cortes T, Moya A, Nakabachi A, Philippe C, Prunier-Leterme N, Rahbe Y, Simon JC, Stern DL, Wincker P, Tagu D: Large-scale gene discovery in the pea aphid Acyrthosiphon pisum (Hemiptera). Genome Biol. 2006, 7: R21-10.1186/gb-2006-7-3-r21.
Rahbe Y, Delobel B, Febvay G, Chantegrel B: Aphid-specific triglycerides in symbiotic and aposymbiotic Acyrthosiphon pisum. Insect Biochem Mol Biol. 1994, 24: 95-101. 10.1016/0965-1748(94)90127-9.
Chernysh S, Cociancich S, Briand JP, Hetru C, Bulet P: The inducible antibacterial peptides of the hemipteran insect Palomena prasina: Identification of a unique family of proline-rich peptides and of a novel insect defensin. J Insect Physiol. 1996, 42: 81-89. 10.1016/0022-1910(95)00085-2.
Sonenshine DE, Hynes WL: Molecular characterization and related aspects of the innate immune response in ticks. Front Biosci. 2008, 13: 7046-7063. 10.2741/3209.
Dall'Antonia Y, Pavkov T, Fuchs H, Breiteneder H, Keller W: Crystallization and preliminary structure determination of the plant food allergen Pru av 2. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2005, 61: 186-188. 10.1107/S1744309104033822.
Koga R, Tsuchida T, Sakurai M, Fukatsu T: Selective elimination of aphid endosymbionts: effects of antibiotic dose and host genotype, and fitness consequences. FEMS Microbiol Ecol. 2007, 60: 229-239. 10.1111/j.1574-6941.2007.00284.x.
Harada H, Ishikawa H: Experimental pathogenicity of Erwinia aphidicola to pea aphid, Acyrthosiphon pisum. J Gen Appl Microbiol. 1997, 43: 363-367. 10.2323/jgam.43.363.
Rahbe Y, Febvay G: Protein toxicity to aphids - an invitro test on Acyrthosiphon pisum. Entomologia Experimentalis Et Applicata. 1993, 67: 149-160. 10.1007/BF02386520.
Mietkiewski R, Soper RS, Balazy S: Notes on Zoophthora occidentalis (Thaxter) Batko (entomophthorales, entomophthoraceae). Mycotaxon. 1981, 13: 41-49.
Williamson C, Rybicki EP, Kasdorf GGF, Vonwechmar MB: Characterization of a new picorna-like virus isolated from aphids. J Gen Virol. 1988, 69: 787-795. 10.1099/0022-1317-69-4-787.
van Munster M, Dullemans AM, Verbeek M, Heuvel van den J, Clerivet A, Wilk van der F: Sequence analysis and genomic organization of Aphid lethal paralysis virus: a new member of the family Dicistroviridae. J Gen Virol. 2002, 83: 3131-3138.
Wilson ACC, Dunbar HE, Davis GK, Hunter WB, Stern DL, Moran NA: A dual-genome microarray for the pea aphid, Acyrthosiphon pisum, and its obligate bacterial symbiont, Buchnera aphidicola. BMC Genomics. 2006, 7: 50-10.1186/1471-2164-7-50.
Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, NJ: Humana Press, 365-386.
Pea Aphid Gbrowse. [http://genoweb1.irisa.fr/cgi-bin/gbrowse/gbrowse/aphidbase/]
AphidBase Downloads. [http://www.aphidbase.com/aphidbase/downloads]
Acyrthosiphon pisum EST Database. [http://aphidests.org/]
We thank Angela Douglas, Nancy Moran, Tom Little and members of the International Aphid Genomics Consortium for insightful discussion, and Charles Godfray and two anonymous reviewers for comments that enhanced this manuscript. Cultures of Z. occidentalis were provided by the USDA ARS Collection of Entomopathogenic Fungal Cultures. Samples of ALPV and virus-infection protocols were provided kindly by Bryony Bonning and Liljana Georgievska. Comparison with M. persicae was supported by USDA grant 2005-35604-15446 to Georg Jander.
NMG, SMB, and MG were group leaders for the project. NMG, BA, HA, SMB, MDV, EJD, JDE, AM, MG, IK, AN, BJP, MP, JSR, JT, DT, and CT designed and performed manual gene annotation. TG and SMB conducted phylogenetic analyses. BA and AV conceived of and conducted analyses of Thaumatin. SMB, NMG, CS and BJP performed experiments and analyses for the gene expression study. CA, AH, VPB, AM, and AL conceived of and conducted the SSH study, and CVM constructed the aphid gut libraries. YR conducted the HPLC study. The manuscript was prepared by NMG, SMB, CA, TG and YR with input from MDV, BA, AN, AV and AH. All authors have read and approved the final version of the manuscript.