- Research article
- Open Access
Aphids acquired symbiotic genes via lateral gene transfer
BMC Biology volume 7, Article number: 12 (2009)
Aphids possess bacteriocytes, which are cells specifically differentiated to harbour the obligate mutualist Buchnera aphidicola (γ-Proteobacteria). Buchnera has lost many of the genes that appear to be essential for bacterial life. From the bacteriocyte of the pea aphid Acyrthosiphon pisum, we previously identified two clusters of expressed sequence tags that display similarity only to bacterial genes. Southern blot analysis demonstrated that they are encoded in the aphid genome. In this study, in order to assess the possibility of lateral gene transfer, we determined the full-length sequences of these transcripts, and performed detailed structural and phylogenetic analyses. We further examined their expression levels in the bacteriocyte using real-time quantitative RT-PCR.
Sequence similarity searches demonstrated that these fully sequenced transcripts are significantly similar to the bacterial genes ldcA (product, LD-carboxypeptidase) and rlpA (product, rare lipoprotein A), respectively. Buchnera lacks these genes, whereas many other bacteria, including Escherichia coli, a close relative of Buchnera, possess both ldcA and rlpA. Molecular phylogenetic analysis clearly demonstrated that the aphid ldcA was derived from a rickettsial bacterium closely related to the extant Wolbachia spp. (α-Proteobacteria, Rickettsiales), which are intracellular symbionts of various lineages of arthropods. The evolutionary origin of rlpA was not fully resolved, but it was clearly demonstrated that its double-ψ β-barrel domain is of bacterial origin. Real-time quantitative RT-PCR demonstrated that ldcA and rlpA are expressed 11.6 and 154-fold higher in the bacteriocyte than in the whole body, respectively. LdcA is an enzyme required for recycling murein (peptidoglycan), which is a component of the bacterial cell wall. As Buchnera possesses a cell wall composed of murein but lacks ldcA, a high level of expression of the aphid ldcA in the bacteriocyte may be essential to maintain Buchnera. Although the function of RlpA is not well known, conspicuous up-regulation of the aphid rlpA in the bacteriocyte implies that this gene is also essential for Buchnera.
In this study, we obtained several lines of evidence indicating that aphids acquired genes from bacteria via lateral gene transfer and that these genes are used to maintain the obligately mutualistic bacterium, Buchnera.
Aphids are hemipteran insects that have close associations with various lineages of microorganisms. Most aphid species harbour the obligate mutualist (usually called primary symbiont), Buchnera aphidicola (γ-Proteobacteria), within the cytoplasm of specialized cells called bacteriocytes [1–4]. Since the initial infection more than 100 million years ago , Buchnera have been subjected to strict vertical transmission through host generations, and the mutualism between Buchnera and their host has evolved to the point that neither can reproduce in the absence of the other. Buchnera cannot proliferate outside bacteriocytes and, when deprived of Buchnera, the host insects suffer retarded growth and sterility, as they are obligately dependent on Buchnera for the supply of essential nutrients they cannot synthesize, and which are scarce in their diet of phloem sap [4, 6–9]. During the process of co-evolution with the host, Buchnera has lost a number of genes that appear to be essential for bacterial existence (the genomes of Buchnera range from 420 to 650 kb in size and contain 400 to 600 genes; [10–13]); this raises the question of how Buchnera survive within the host bacteriocyte. One of several possible explanations for the absence of these genes is that some genes were transferred from the genome of an ancestor of Buchnera to the genome of an aphid ancestor and are now expressed under the control of the host nucleus. Such lateral gene transfer (LGT) should take place in the germ line for the transferred gene to be inherited through the generations of the recipient. During most of their life stages, Buchnera are confined within bacteriocytes, which are segregated from germ cells; however, the symbionts are freed from the maternal bacteriocytes before being transmitted to the next generation. In cases of parthenogenetic reproduction, Buchnera cells are transferred into the parthenogenetic blastoderm-stage embryos; Buchnera are localized proximal to the host germ cells during early development of the host. Moreover, in cases of sexual reproduction, Buchnera enter sexual eggs at the pre-cellularization stage; at this stage, there are no membranous barriers between Buchnera and the germ lines [14, 15]. Such localization of Buchnera cells proximal to host germ lines might provide opportunities for the LGT from Buchnera into the germ lines.
In addition to Buchnera, a number of aphid strains harbour other maternally transmitted intracellular bacteria, such as Rickettsia (α-Proteobacteria), Spiroplasma (Mollicutes), and various γ-proteobacterial microbes, including Hamiltonella defensa, Regiella insecticola, Serratia symbiotica, and Arsenophonus species [16–28]. These 'secondary symbionts' are often shared between divergent insect lineages. For example, Hamiltonella and Arsenophonus are observed in scattered strains and species of aphids, psyllids, whiteflies and planthoppers [22, 24, 29]. Wolbachia lineages (α-Proteobacteria, Rickettsiales) are observed in a wide variety of arthropods [30–32], though only one case of infection has been reported in aphids . These suggest that secondary symbionts undergo horizontal transfer among matrilines within and between species. They are also transmitted vertically [16–28], but this appears to be achieved in a less tightly controlled manner in comparison to the case of Buchnera . Whereas Buchnera exist as passive symbionts within their hosts, which in turn have evolved mechanisms to maintain and transmit the Buchnera [14, 15, 33], secondary symbionts overcome host immune responses and invade various types of host cells, including germ cells [1, 18, 21, 23, 25, 27, 34]. Thus, there are likely to have been frequent opportunities for aphids to acquire genomic fragments from these symbiotic bacteria during evolution.
We previously performed transcriptome analysis of the bacteriocyte of the pea aphid Acyrthosiphon pisum, to elucidate the host mechanisms required to maintain Buchnera . This study identified a number of aphid genes that are highly expressed in the bacteriocyte. Among them, two genes (corresponding to the cDNA clusters R2C00193 and R2C00214, which consist of 10 and four expressed sequence tags, respectively) exhibited similarity only to prokaryotic genes, and not to those of extant Buchnera lineages. Southern blot analysis confirmed that they are encoded in the aphid genome.
In the present study, we show the detailed analysis of the phylogenetic positions, domain structures, and expression profiles of these genes, thus revealing their evolutionary history and functional roles.
Full-length sequencing of cDNA clones
In the previous study, the sequences of the transcripts corresponding to the cDNA clusters (unigenes) R2C00193 and R2C00214 were not fully determined, as the cap-trapper cDNA clones were sequenced only from the 5' end. In the present study, all the cap-trapper clone inserts relevant to these unigenes were amplified by PCR using vector primers (T3 and T7) and sequenced from both ends to obtain full-length sequences. In the case of R2C00214, all of the four clones had an identical sequence of 1312 bp encoding a polypeptide of 226 amino acid residues (Figure 1). Full-length sequences for R2C00193 were approximately 1 kb in length, with slight variations mainly in the putative untranslated regions (UTRs). They encoded a polypeptide of 220 amino acid residues (Figure 2). These full-length unigenes are hereafter referred to as R2C00214F (DDBJ: AB435382) and R2C00193F (DDBJ: AB435384 and AB435385), respectively.
Basic Local Alignment Search Tool (BLAST) search demonstrated that the product of R2C00214F has significant sequence similarity to the bacterial enzyme LD-carboxypeptidase (LdcA), and the microcin C7 self-immunity protein (MccF) that are produced by Gram-negative bacteria (Figure 1A). The top BLAST hit for the R2C00214F product was the hypothetical protein WD1015 [Wolbachia endosymbiont of Drosophila melanogaster (Wolbachia wMel; α-Proteobacteria, Rickettsiales)] (RefSeq: NP_966741) (E = 1 × 10-23), which has not been fully annotated; however, the analysis of the conserved domains of the gene product performed using the CD-search at the National Center for Biotechnology Information (NCBI) website indicated that the gene encodes the bacterial LdcA belonging to the peptidase S66 family (pfam02016, E = 1 × 10-47). The subordinate hits were either LdcA or MccF, the latter of which mediates resistance against microcin C7, an antimicrobial peptide that is secreted by enterobacteria and inhibits the growth of bacterial species phylogenetically related to the producing strains . The mechanism through which MccF mediates resistance against microcin C7 is uncertain; however, MccF belongs to the peptidase S66 family, and all the residues required for LdcA activity are conserved in it . Thus, in this paper, we collectively refer to these proteins belonging to the S66 family as 'LdcAs'. Putative orthologs of R2C00214F are found in a variety of bacteria, but not in eukaryotes, except for the fungus Gibberella zeae (RefSeq: XP_383840), implying that the two distantly related organisms, namely, the aphid and the plant pathogenic fungi, independently acquired ldcA from a bacterium via LGT. We discuss this possible mode of inheritance via LGT below.
R2C00214F appeared to lack the sequences required to encode the middle region of canonical LdcAs (Figure 1). To check the corresponding genomic sequences, the preliminary genome assembly of the pea aphid (Acyr_1.0) was screened using R2C00214F as the query sequence. The entire coding sequence (CDS) of ldcA was located in a single scaffold; however, the genomic DNA had a sequence corresponding to the middle region of the LdcAs that was missing from the R2C00214F gene product. This suggests that the sequence fragment represents an intron of the R2C00214F gene (Figure 1). In order to search for splice variants, we further amplified cDNAs for the aphid LdcA by RT-PCR using specific primers and determined their sequences. Unexpectedly, the amplified cDNAs essentially consisted of a single type of sequence variant that contained a sequence (402 bp = 134 amino acids) corresponding to the middle region of LdcA. This long form of the transcript (DDBJ: AB435383) encoded a 360-amino acid-long polypeptide sequence, while the short form (R2C00214F, DDBJ: AB435382) encoded a 226-amino acid-long polypeptide sequence (Figure 1). The long form and the short form appeared to be splice variants as cap-trapper libraries rarely contain inappropriate artifacts that do not reflect the mRNA structures in vivo [37, 38]. The terminal dinucleotides of the insertion sequence were GT-CG, which is similar to the canonical splicing signal GT-AG. It has previously been verified that the GT-CG set can also be used as a splicing signal . The short form of the transcript was not detected by RT-PCR; this might, at least in part, reflect the low level of expression of this truncated form of the transcript in the sample used in the study. The inconsistency might also be due to some bias in constructing the cDNA library in the previous study and/or in the RT-PCR employed in the present study.
Figure 1B shows the intron-exon boundaries of the aphid ldcA gene. The long-form transcript consists of two exons and a single intron. Only the second exon encodes the open reading frame (ORF) of the protein. The short-form transcript consists of three exons and two introns; the middle region of the second exon of the long-form transcript is spliced out as the second intron. These exon/intron organizations were verified by PCR cloning.
The long-form transcript was also characterized by BLAST similarity search. Once again, the top BLAST hit was the hypothetical protein WD1015 (Wolbachia wMel) (RefSeq: NP_966741) (E = 5 × 10-81). The subordinate hits were similar to those obtained with the short-form transcripts, but with much smaller E-values. The amino acid sequence of the long-form transcript exhibited 45% and 24% identity to the LdcA proteins of Wolbachia wMel and Escherichia coli, respectively (Figure 1A). Three catalytically active sites identified in Pseudomonas LdcA (Ser-126, His-304, and Glu-227)  were conserved in the aphid LdcA (Figure 1A). No other domain structure was observed in the protein.
Putative rare lipoprotein A
The BLAST search revealed that the R2C00193F gene product is significantly similar to a bacterial protein, rare lipoprotein A (RlpA) (Figure 2A, 2B). The top BLAST hit was a putative RlpA family protein [Bradyrhizobium sp. BTAi1 (α-Proteobacteria, Rhizobiales)] (RefSeq: YP_001239851) (E = 9 × 10-15), and essentially all of the subordinate hits were thus annotated rare lipoprotein A. Homologous sequences of the pea aphid putative rlpA gene were observed in various bacteria, but not in eukaryotes, except for two other aphid species, Aphis gossypii (GenBank: DR391796) and Toxoptera citricida (GenBank: CD450666). Domain analysis revealed that the region detected by the similarity search corresponds to the double-ψ β-barrel (DPBB) fold, which is the domain conserved in RlpA proteins. Although the function of RlpA is not well understood , the DPBB fold is suspected to be an enzymatic domain .
Using RT-PCR cloning, two types of sequences were identified. As expected, these sequences corresponded to the transcripts originally found in the sequence cluster of R2C00193F (DDBJ: AB435384 and AB435385). These contained putative full CDSs encoding 220-amino acid polypeptide sequences (Figure 2A, 2B). These sequences appeared to be from distinct alleles, with two nucleotide discrepancies in their CDSs resulting in a single amino acid difference (denoted with an asterisk in Figure 2A).
Three other domain structures were observed in the pea aphid putative RlpA (Figure 2A, 2B). At the N-terminal region, a eukaryotic signal peptide motif was identified. BLAST search of the remaining sequences revealed that two regions adjacent to the DPBB domain are similar to the inhibitor cysteine-knot (ICK) motif of three antimicrobial peptides – Alo-1, Alo-2, and Alo-3 – of the harlequin beetle Acrocinus longimanus (Swiss-Prot: P83651, P83652, P83653)  (Figure 2C). The ICK motif presents a unique knotted topology of three disulphide bridges, with one disulphide penetrating through a macrocycle formed by the other two disulphides and interconnecting the peptide backbones . The ICK family proteins are relatively small (typically less than 40 residues in length), and are found in various lineages of eukaryotes including plants, molluscs, arachnids and insects, exhibiting various biological activities such as toxic, antimicrobial, and insecticidal activities [42, 43]. This motif was observed also in the putative ORFs of two other aphid transcripts. However, the domain has never been found in bacterial proteins, including RlpA.
To reveal the exon/intron structure of the pea aphid putative rlpA, a preliminary genome assembly of the pea aphid (Acyr_1.0) was screened using R2C00193F as the query sequence. The pea aphid rlpA locus was split into two distinct scaffolds (Figure 2B). The rlpA gene consists of three exons and two introns. The first exon contains the eukaryotic signal peptide, the second contains one of the cysteine-rich domains, and the third contains the DPBB domain and another cysteine-rich domain. The boundaries of the protein domains were consistent with the locations of introns. The chimeric structure of the aphid RlpAs might have come into being as the result of exon-shuffling  involving prokaryotic and eukaryotic elements.
Pea aphid ancestor acquired ldcA via LGT from a Wolbachia-like bacterium
The amino acid sequence of the aphid putative LdcA was subjected to molecular phylogenetic analysis (Figure 3). The phylogenetic tree demonstrated with robust statistical support (98% in NJ; 97% in ML; 1.0 in Bayesian) that the aphid gene is most closely related to the clade of LdcAs of rickettsial bacteria, especially Wolbachia (RefSeq: NP_966741) and Orientia tsutsugamushi (RefSeq: YP_001248242). This branching pattern can be most simply explained by the hypothesis that the aphid acquired ldcA via LGT from Wolbachia or some other rickettsial bacteria, many of which are known to be intracellular symbionts of insects. The putative orthologous gene detected in the plant pathogenic fungus G. zeae was distantly related to the aphid ldcA (Figure 3); this further suggests that the ancestors of A. pisum and G. zeae independently acquired the genes from different lineages of bacteria.
Common ancestor of three species of aphids acquired rlpAvia LGT
The amino acid sequences of putative RlpAs of A. pisum, A. gossypii, and T. citricida were subjected to molecular phylogenetic analysis with RlpAs of various bacterial lineages (Figure 4). The highly conserved DPBB domains were aligned and used for this analysis. The phylogenetic positions of aphid RlpAs were not clearly resolved with a high level of statistical support. However, to date, no rlpA genes have been observed in any eukaryotes, except aphids. Moreover, the phylogenetic tree showed that the aphid rlpAs are monophyletic and that the phylogenetic relationships were congruent with the species tree of aphids . This suggests that the common ancestor of these three aphid species acquired the rlpA gene from a bacterium via LGT. The relatively low resolution of the phylogenetic positions of the aphid rlpA may be partly due to the high evolutionary rate of the aphid lineages (see below).
Aphid rlpAencodes a functional protein
To test the integrity of the functionality of the aphid RlpAs, K A /K S ratios of the DPBB-encoding sequences of aphid rlpAs were calculated. The ratios between A. pisum and A. gossypii, and between A. pisum and T. citricida were 0.45 and 0.30, respectively. Both of the ratios were significantly smaller than 1 (p-values were 0.02 and 0.001, respectively). This indicates that the aphid genes are not pseudogenes, but are functional and contribute to the fitness of the insects. However, the K A /K S values were somewhat higher than those of other bacteria [for example, 0.03, Nitrosomonas eutropha (RefSeq: NC_008344) vs. Nitrosomonas europaea (RefSeq: NC_004757); 0.04, Yersinia mollaretii (RefSeq: NZ_AALD01000017) vs. Yersinia bercovieri (RefSeq: NZ_AALC01000036); the K S values between other pairs were saturated]. This indicates that the selective constraints on amino acid substitutions in the DPBB domains are 8 to 15 times more relaxed in aphids than bacteria, under the assumption that synonymous sites evolve neutrally.
ldcA and rlpAare highly expressed specifically in the bacteriocyte
To examine the expression profiles of ldcA and rlpA, we quantified their transcripts in the bacteriocyte and in the whole body using real-time quantitative RT-PCR (Figure 5). The results clearly demonstrated that ldcA and rlpA are actively transcribed in the bacteriocyte. Transcripts for ldcA and rlpA were 11.6 and 154-fold more abundant in the bacteriocyte than in the whole body, respectively. It is also notable that the copy numbers of their transcripts in the bacteriocyte were comparable to those of the control transcript encoding ribosomal protein L7 (RpL7), indicating that their expression levels are relatively high. High levels of expression of these genes in the bacteriocyte strongly suggest that they are not only functional, but they play important roles in maintaining the symbiotic relationship with the obligate mutualist, Buchnera.
Aphids have recruited genes from bacterial genomes via LGT
We have obtained strong evidence for two cases of bacteria-to-insect LGT. The genes encoded in the aphid genome that are expressed in the bacteriocyte were demonstrated to be significantly similar only to the bacterial genes, ldcA and rlpA. Quantitative RT-PCR further verified that these genes are highly expressed in the bacteriocyte. The orthologs of such genes are absent in Buchnera, the obligate mutualistic bacteria that are harboured in the bacteriocytes. These findings imply that the aphid ldcA and rlpA have compensational functions to support the survival of Buchnera.
Although until recently it was believed that LGT plays an important role exclusively in the evolution of unicellular organisms, especially prokaryotes [46–48], the accumulating genomic data is now revealing that LGT also affects the genomic content of multicellular eukaryotes with segregated germ cells. DNA sequences with significant similarity to genes of Wolbachia, an endocellular rickettsial bacterium, have been observed in the genomes of a wide range of arthropods and filarial nematodes [49–52]. Wolbachia is a maternally transmitted endosymbiont that can enter the germ line of host animals , which facilitates bacterial DNA transfer to the host nucleus. However, many of the transferred Wolbachia genes appear to be in the process of pseudogenization, and even intact "genes" are not expressed at a significant level [51, 52], implying that these transferred genes do not confer novel functions on the host organisms.
The cases discovered in the present study are especially interesting in that these transferred bacterial genes not only retain their functionality, but are highly expressed in the bacteriocyte, which is the cell that harbours Buchnera. The molecular phylogenetic analysis clearly indicated that the aphid ldcA is closely related to that of Wolbachia, and of other rickettsial bacteria. Although infections of Wolbachia and Rickettsia are sporadically observed among the aphid species [17, 23, 25], the ISO strain that was used in the present study lacks such symbionts ; this suggests that the previous infection left only a transferred gene as a footprint, while the source bacterium disappeared. With regard to rlpA, it was clearly demonstrated that this gene also was of bacterial origin, but its phylogenetic position has not been fully resolved.
Eukaryote-type structures of the genes and transcripts
Recent studies have revealed that LGTs from bacteria can occur in metazoa [49–52]. However, these transferred genes cannot function unless they obtain eukaryotic promoters, since the gene expression systems of prokaryotes and eukaryotes differ. The likelihood of promoter acquisition seems very low, as suggested by the previously reported lack of expression of laterally transferred genes [51, 52]. The aphid ldcA and rlpA are highly expressed in the bacteriocytes, clearly indicating that these genes have acquired eukaryotic promoters, although the mechanism of promoter acquisition has yet to be determined. The cDNAs for the aphid ldcA and rlpA were originally observed in a cDNA library constructed by the cap-trapper method that targets the 5' cap structure and 3' poly-A tails of eukaryotic mRNAs [37, 38]. This suggests that both the mRNAs for the aphid ldcA and rlpA have the 5' cap structure and 3' poly-A tail that bacterial mRNAs lack. Indeed, polyadenylation signals (AAUAAA) are observed in the 3'-UTR of the transcripts. Screening of the genome scaffold followed by PCR cloning revealed that the genes have spliceosomal-type introns. This type of intron has not been observed in bacterial genes, suggesting that these genes acquired introns after they were transferred into the aphid genome.
LdcA may be used to control Buchnera
LdcA is an enzyme required for recycling murein (peptidoglycan), a component of the bacterial cell wall. LdcA releases the terminal D-alanine from L-alanyl-D-glutamyl-meso-diaminopimelyl-D-alanine, which contains turnover products of murein. The disruption of E. coli ldcA results in bacteriolysis during the stationary phase, indicating that the reaction is essential for bacterial survival . In the mutant, due to a defect in murein recycling, the unusual murein precursor uridine 5'-pyrophosphoryl N-acetylmuramyl-tetrapeptide accumulates, and the overall cross-linkage of murein decreases dramatically. This is interpreted as a reflection of the increased incorporation of tetrapeptide precursors that can only function as acceptors and not as donors in the cross-linking reaction.
Buchnera has cell walls composed of murein , but it lacks ldcA . Although the evolutionary origin of the aphid ldcA seems to be from rickettsial bacteria and not from Buchnera, it is intriguing to note that this gene is highly expressed in the bacteriocyte. Aphids may control the proliferation of Buchnera using ldcA, which was recruited from another symbiotic bacterium that previously had resided in aphids.
Chimeric structure of putative RlpA
The molecular phylogenetic tree indicated that the LGT of rlpAs occurred before the divergence of the three aphid species. On the basis of fossil records, this divergence is inferred to date back to more than 50 million years ago [45, 55]. Even if the transferred genes successfully acquire sequence elements that allow their expression, contribution of the genes to the host fitness, or strategies enabling the selfish propagation of the genes, would be required for the maintenance of the transferred genes in the host genome for such a long period of time. The functional role of the rlpA in any bacteria is not well understood; however, RlpA suppresses the E. coli mutant of Prc that cleaves the C terminus of FtsI , suggesting that rlpA plays an important role in bacteria. Domain analyses revealed that, in addition to the conserved DPBB domain, the aphid RlpA has three other domains that are not found among bacterial orthologs. This implies that RlpA might have gained novel functions that are yet to be determined. Although the function of RlpA is not well understood, the high level of expression of the aphid rlpA in the bacteriocyte implies that this gene is also essential for Buchnera.
In this study, several lines of evidence indicated that aphids acquired genes from bacteria via LGT, and are using such genes to maintain the obligately mutualistic bacteria, Buchnera. Phylogenetic analysis clearly demonstrated that one of the genes was derived from a rickettsial bacterium that is closely related to the extant Wolbachia. This is the first report of functional genes that were laterally transferred from symbiotic bacteria to metazoa. The cases presented here are of special interest in that these transferred bacterial genes not only retain their functionality, but are highly expressed in the bacteriocyte that is differentiated so as to harbour Buchnera, which lack such genes.
Strain ISO, a parthenogenetic clone of the pea aphid Acyrthosiphon pisum which is free from secondary symbionts including Wolbachia, was used in this study. The insects were reared on Vicia faba at 15°C in a long-day regime of 16 hr light and 8 hr dark. Parthenogenetic apterous adults (12 to 15 days old) were used for the experiments.
Cloning of the pea aphid genes
Genomic DNA was extracted from the whole body of the pea aphid using a DNeasy Kit (Qiagen). Total RNA was extracted from bacteriocytes and first strand cDNA was prepared as described previously . PCR was performed using various sets of gene specific primers (Additional file 1). PCR products were either purified and sequenced directly or cloned using the pGEM-T easy vector system (Promega).
Characterization of gene products by similarity searches
Homologous protein sequences and conserved domains were detected by BLASTP similarity searches at the website of the NCBI using deduced amino acid sequences as queries . The presence and location of signal peptides were predicted using the program SignalP 3.0 . Statistical tests of homology between two amino acid sequences were conducted with bl2seq. Default parameters were used except for the matrix (-M), set to BLOSUM80, the gap existence penalties (-G), set to 11, and the theoretical database size (-d), set to 127,836,513, the size of Swiss-Prot release 55.0.
Molecular phylogenetic analysis
Multiple protein sequences were aligned using the program package MAFFT 5.8 , followed by manual refinement. Amino acid sites corresponding to alignment gap(s) were omitted from the data set. Only unambiguously aligned amino acid sequences were used for the phylogenetic analysis. The aligned sequence data are shown in Additional file 2.
Phylogenetic trees were inferred by the neighbour joining, the maximum likelihood and the Bayesian methods. Neighbour-joining trees were constructed using the program package Xced . The distance matrix was estimated by the maximum likelihood distance method assuming the JTT model with among-site rate heterogeneity. The bootstrap probability for each node was calculated by generating 1000 bootstrap replicates. Maximum likelihood trees were estimated using the program package RAxML . In the analysis, the JTT model was used as a substitution model for amino acids. To incorporate the effect of among-site rate heterogeneity, a mixed model (one invariable rate plus Γ distributed rates) was used. The support values for the internal nodes were inferred by 1000 bootstrap replicates. In the Bayesian inference, we used the program MrBayes 3.1.2 . The JTT +Γ +Inv model was used as a substitution model. In total, 4100 trees were obtained (ngen 410,000, samplefreq 100), and the first 2000 of these were considered as the 'burn in' and discarded. We checked that the potential scale reduction factor was approximately 1.00 for all parameters and that the average standard deviation of split frequencies converged towards zero. K S and K A values were calculated as described previously . Statistical significance of the obtained K A /K S value was tested against a bootstrap distribution of K A /K S values, which was generated by 10,000 bootstrap resamplings of codons from the original alignment.
Real-time quantitative RT-PCR
RNA was isolated from whole bodies and bacteriocytes of 12 to 15-day-old parthenogenetic apterous adults using TRIzol reagent, followed by RNase-free DNase I treatment. Each whole body sample and bacteriocyte sample was derived from one individual and a batch of bacteriocytes that were collected from about 10 individuals, respectively. First-strand cDNAs were synthesized using pd(N)6 primer and PrimeScript reverse transcriptase (Takara). Quantification was performed with the LightCycler instrument and FastStart DNA MasterPLUS SYBR Green I kit (Roche), as described previously . The primers used were: ldcA-677F (CAACCTGACGCTAGTCGAGAACT), ldcA-758R (CACGTCCTCCAAGAACACGAT), rlpA-449F (CGGCGGACGGTAAGGTAAT), and rlpA-529R (ACTGTACCGGGCCTGTGTTC). The running parameters were: 95°C for 10 min, followed by 45 cycles of 95°C for 10 s, 55°C for 5 s, and 72°C for 4 s. Results were analyzed using the LightCycler software version 3.5 (Roche), and relative expression levels were normalized to mRNA for the ribosomal protein RpL7. Statistical analyses were performed using the Mann-Whitney U-test.
Source of the genomic data of the pea aphid
The preliminary genome assembly of A. pisum was obtained from the Human Genome Sequencing Center at Baylor College of Medicine through the web site at http://www.hgsc.bcm.tmc.edu/projects/aphid/.
Buchner P: Endosymbiosis of animals with plant microorganisms. 1965, New York: Interscience
Munson MA, Baumann P, Kinsey MG: Buchnera gen. nov. and Buchnera aphidicola sp. nov., a taxon consisting of the mycetocyte-associated, primary endosymbionts of aphids. Int J Syst Bacteriol. 1991, 41: 566-568.
Baumann P, Baumann L, Lai CY, Rouhbakhsh D, Moran NA, Clark MA: Genetics, physiology, and evolutionary relationships of the genus Buchnera: intracellular symbionts of aphids. Annu Rev Microbiol. 1995, 49: 55-94. 10.1146/annurev.mi.49.100195.000415.
Douglas AE: Nutritional interactions in insect-microbial symbioses: aphids and their symbiotic bacteria Buchnera. Annu Rev Entomol. 1998, 43: 17-37. 10.1146/annurev.ento.43.1.17.
Moran NA, Munson MA, Baumann P, Ishikawa H: A molecular clock in endosymbiotic bacteria is calibrated using the insect hosts. P Roy Soc Lond B Bio. 1993, 253: 167-171. 10.1098/rspb.1993.0098.
Febvay G, Liadouze I, Guillaud J, Bonnot G: Analysis of energetic amino acid metabolism in Acyrthosiphon pisum: a multidimensional approach to amino acid metabolism in aphids. Arch Insect Biochem. 1995, 29: 45-69. 10.1002/arch.940290106.
Sasaki T, Ishikawa H: Production of essential amino acids from glutamate by mycetocyte symbionts of the pea aphid, Acyrthosiphon pisum. J Insect Physiol. 1995, 41: 41-46. 10.1016/0022-1910(94)00080-Z.
Nakabachi A, Ishikawa H: Differential display of mRNAs related to amino acid metabolism in the endosymbiotic system of aphids. Insect Biochem Mol Biol. 1997, 27: 1057-1062. 10.1016/S0965-1748(97)00092-1.
Nakabachi A, Ishikawa H: Provision of riboflavin to the host aphid, Acyrthosiphon pisum, by endosymbiotic bacteria, Buchnera. J Insect Physiol. 1999, 45: 1-6. 10.1016/S0022-1910(98)00104-8.
Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H: Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature. 2000, 407: 81-86. 10.1038/35024074.
Tamas I, Klasson L, Canback B, Naslund AK, Eriksson AS, Wernegreen JJ, Sandstrom JP, Moran NA, Andersson SG: 50 million years of genomic stasis in endosymbiotic bacteria. Science. 2002, 296: 2376-2379. 10.1126/science.1071278.
van Ham RC, Kamerbeek J, Palacios C, Rausell C, Abascal F, Bastolla U, Fernandez JM, Jimenez L, Postigo M, Silva FJ, Tamames J, Viguera E, Latorre A, Valencia A, Morán F, Moya A: Reductive genome evolution in Buchnera aphidicola. Proc Natl Acad Sci USA. 2003, 100: 581-586. 10.1073/pnas.0235981100.
Perez-Brocal V, Gil R, Ramos S, Lamelas A, Postigo M, Michelena JM, Silva FJ, Moya A, Latorre A: A small microbial genome: the end of a long symbiotic relationship?. Science. 2006, 314: 312-313. 10.1126/science.1130441.
Braendle C, Miura T, Bickel R, Shingleton AW, Kambhampati S, Stern DL: Developmental origin and evolution of bacteriocytes in the aphid-Buchnera symbiosis. PLoS Biol. 2003, 1: E21-10.1371/journal.pbio.0000021.
Miura T, Braendle C, Shingleton A, Sisk G, Kambhampati S, Stern DL: A comparison of parthenogenetic and sexual embryogenesis of the pea aphid Acyrthosiphon pisum (Hemiptera: Aphidoidea). J Exp Zoolog B Mol Dev Evol. 2003, 295: 59-81.
Unterman BM, Baumann P, McLean DL: Pea aphid symbiont relationships established by analysis of 16S rRNAs. J Bacteriol. 1989, 171: 2970-2974.
Chen DQ, Campbell BC, Purcell AH: A new rickettsia from a herbivorous insect, the pea aphid Acyrthosiphon pisum (Harris). Curr Microbiol. 1996, 33: 123-128. 10.1007/s002849900086.
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.
Fukatsu T, Tsuchida T, Nikoh N, Koga R: Spiroplasma symbiont of the pea aphid, Acyrthosiphon pisum (Insecta: Homoptera). Appl Environ Microbiol. 2001, 67: 1284-1291. 10.1128/AEM.67.3.1284-1291.2001.
Darby AC, Douglas AE: Elucidation of the transmission patterns of an insect-borne bacterium. Appl Environ Microbiol. 2003, 69: 4403-4407. 10.1128/AEM.69.8.4403-4407.2003.
Koga R, Tsuchida T, Fukatsu T: Changing partners in an obligate symbiosis: a facultative endosymbiont can compensate for loss of the essential endosymbiont Buchnera in an aphid. Proc Biol Sci. 2003, 270 (1533): 2543-2550. 10.1098/rspb.2003.2537.
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.
Gomez-Valero L, Soriano-Navarro M, Perez-Brocal V, Heddi A, Moya A, Garcia-Verdugo JM, Latorre A: Coexistence of Wolbachia with Buchnera aphidicola and a secondary symbiont in the aphid Cinara cedri. J Bacteriol. 2004, 186: 6626-6633. 10.1128/JB.186.19.6626-6633.2004.
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.
Sakurai M, Koga R, Tsuchida T, Meng XY, Fukatsu T: Rickettsia symbiont in the pea aphid Acyrthosiphon pisum: novel cellular tropism, effect on host fitness, and interaction with the essential symbiont Buchnera. Appl Environ Microbiol. 2005, 71: 4069-4075. 10.1128/AEM.71.7.4069-4075.2005.
Tsuchida T, Koga R, Fukatsu T: Host plant specialization governed by facultative symbiont. Science. 2004, 303: 1989-10.1126/science.1094611.
Moran NA, Dunbar HE: Sexual acquisition of beneficial symbionts in aphids. Proc Natl Acad Sci USA. 2006, 103: 12803-12806. 10.1073/pnas.0605772103.
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.
Moran NA, McCutcheon JP, Nakabachi A: Genomics and evolution of heritable bacterial symbionts. Annu Rev Genet. 2008, 42: 165-190. 10.1146/annurev.genet.41.110306.130119.
Jeyaprakash A, Hoy MA: Long PCR improves Wolbachia DNA amplification: wsp sequences found in 76% of sixty-three arthropod species. Insect Mol Biol. 2000, 9: 393-405. 10.1046/j.1365-2583.2000.00203.x.
Werren JH, Windsor DM: Wolbachia infection frequencies in insects: evidence of a global equilibrium?. Proc Biol Sci. 2000, 267: 1277-1285. 10.1098/rspb.2000.1139.
Knight J: Meet the Herod bug. Nature. 2001, 412: 12-14. 10.1038/35083744.
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.
Serbus LR, Sullivan W: A cellular basis for Wolbachia recruitment to the host germline. PLoS Pathog. 2007, 3: e190-10.1371/journal.ppat.0030190.
Duquesne S, Petit V, Peduzzi J, Rebuffat S: Structural and functional diversity of microcins, gene-encoded antibacterial peptides from enterobacteria. J Mol Microbiol Biotechnol. 2007, 13: 200-209. 10.1159/000104748.
Korza HJ, Bochtler M: Pseudomonas aeruginosa LD-carboxypeptidase, a serine peptidase with a Ser-His-Glu triad and a nucleophilic elbow. J Biol Chem. 2005, 280: 40802-40812. 10.1074/jbc.M506328200.
Carninci P, Kvam C, Kitamura A, Ohsumi T, Okazaki Y, Itoh M, Kamiya M, Shibata K, Sasaki N, Izawa M, Muramatsu M, Hayashizaki Y, Schneider C: High-efficiency full-length cDNA cloning by biotinylated CAP trapper. Genomics. 1996, 37: 327-336. 10.1006/geno.1996.0567.
Carninci P, Westover A, Nishiyama Y, Ohsumi T, Itoh M, Nagaoka S, Sasaki N, Okazaki Y, Muramatsu M, Schneider C, Hayashizaki Y: High efficiency selection of full-length cDNA by improved biotinylated cap trapper. DNA Res. 1997, 4: 61-66. 10.1093/dnares/4.1.61.
Burset M, Seledtsov IA, Solovyev VV: Analysis of canonical and non-canonical splice sites in mammalian genomes. Nucleic Acids Res. 2000, 28: 4364-4375. 10.1093/nar/28.21.4364.
Bass S, Gu Q, Christen A: Multicopy suppressors of prc mutant Escherichia coli include two HtrA (DegP) protease homologs (HhoAB), DksA, and a truncated RlpA. J Bacteriol. 1996, 178: 1154-1161.
Castillo RM, Mizuguchi K, Dhanaraj V, Albert A, Blundell TL, Murzin AG: A six-stranded double-psi beta barrel is shared by several protein superfamilies. Structure. 1999, 7: 227-236. 10.1016/S0969-2126(99)80028-8.
Barbault F, Landon C, Guenneugues M, Meyer JP, Schott V, Dimarcq JL, Vovelle F: Solution structure of Alo-3: a new knottin-type antifungal peptide from the insect Acrocinus longimanus. Biochemistry. 2003, 42: 14434-14442. 10.1021/bi035400o.
Gracy J, Le-Nguyen D, Gelly JC, Kaas Q, Heitz A, Chiche L: KNOTTIN: the knottin or inhibitor cystine knot scaffold in 2007. Nucleic Acids Res. 2008, 36: D314-319. 10.1093/nar/gkm939.
Gilbert W: Why genes in pieces?. Nature. 1978, 271: 501-10.1038/271501a0.
Heie O: Paleontology and phylogeny. Aphids Their Biology, Natural Enemies and Control. Edited by: Harrewijn P. 1987, Amsterdam: Elsevier, 2A: 367-391.
Jain R, Rivera MC, Lake JA: Horizontal gene transfer among genomes: the complexity hypothesis. Proc Natl Acad Sci USA. 1999, 96: 3801-3806. 10.1073/pnas.96.7.3801.
Rivera MC, Lake JA: The ring of life provides evidence for a genome fusion origin of eukaryotes. Nature. 2004, 431: 152-155. 10.1038/nature02848.
Pallen MJ, Wren BW: Bacterial pathogenomics. Nature. 2007, 449: 835-842. 10.1038/nature06248.
Kondo N, Nikoh N, Ijichi N, Shimada M, Fukatsu T: Genome fragment of Wolbachia endosymbiont transferred to X chromosome of host insect. Proc Natl Acad Sci USA. 2002, 99: 14280-14285. 10.1073/pnas.222228199.
Fenn K, Conlon C, Jones M, Quail MA, Holroyd NE, Parkhill J, Blaxter M: Phylogenetic relationships of the Wolbachia of nematodes and arthropods. PLoS Pathog. 2006, 2: e94-10.1371/journal.ppat.0020094.
Dunning Hotopp JC, Clark ME, Oliveira DC, Foster JM, Fischer P, Torres MC, Giebel JD, Kumar N, Ishmael N, Wang S, Ingram J, Nene RV, Shepard J, Tomkins J, Richards S, Spiro DJ, Ghedin E, Slatko BE, Tettelin H, Werren JH: Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science. 2007, 317: 1753-1756. 10.1126/science.1142490.
Nikoh N, Tanaka K, Shibata F, Kondo N, Hizume M, Shimada M, Fukatsu T: Wolbachia genome integrated in an insect chromosome: evolution and fate of laterally transferred endosymbiont genes. Genome Res. 2008, 18: 272-280. 10.1101/gr.7144908.
Templin MF, Ursinus A, Holtje JV: A defect in cell wall recycling triggers autolysis during the stationary growth phase of Escherichia coli. Embo J. 1999, 18: 4108-4117. 10.1093/emboj/18.15.4108.
Houk EJ, Griffiths GW, Hadjokas NE, Beck SD: Peptidoglycan in the cell wall of the primary intracellular symbiote of the pea aphid. Science. 1977, 198: 401-403. 10.1126/science.198.4315.401.
Clark MA, Moran NA, Baumann P: Sequence evolution in bacterial endosymbionts having extreme base compositions. Mol Biol Evol. 1999, 16: 1586-1598.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.
Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33: 511-518. 10.1093/nar/gki198.
Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30: 3059-3066. 10.1093/nar/gkf436.
Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Miyata T, Yasunaga T: Molecular evolution of mRNA: a method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application. J Mol Evol. 1980, 16: 23-36. 10.1007/BF01732067.
The preliminary genomic sequence data of Acyrthosiphon pisum, which was accomplished with support from National Human Genome Research Institute, were obtained from Human Genome Sequencing Center at Baylor College of Medicine through the web site at http://www.hgsc.bcm.tmc.edu/projects/aphid/. This study was financially supported by the Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science and Technology of Japan (NN) and Research Fellowship of the Japan Society for the Promotion of Science for Young Scientists (AN).
AN conceived the study, performed the molecular work, and contributed to the structural analysis. NN performed the structural and phylogenetic analyses. NN and AN collaboratively designed the study and prepared the manuscript. Both authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Alignments of the amino acid sequences of (A) LD- carboxypeptidase (the product of ldcA) and (B) rare lipoprotein A (the product of rlpA). Asterisks (*) and circles (O) represent amino acid positions that are occupied by identical and chemically similar amino acids for all the sequences compared, respectively. Gaps (-) were inserted to increase sequence similarity. Dots (.) represent residues identical to those of the pea aphid proteins. Ambiguously aligned regions are trimmed and the numbers of the amino acid residues trimmed are shown in parenthesis. The GenBank accession numbers of sequences used are included in the additional file itself. (PDF 51 KB)
About this article
- Basic Local Alignment Search Tool
- Lateral Gene Transfer
- Aphid Species
- Molecular Phylogenetic Analysis
- Basic Local Alignment Search Tool Search