Host-specific co-evolution likely driven by diet in Buchnera aphidicola

BACKGROUND

Russian wheat aphid (Diuraphis noxia Kurd.) is a severe pest to wheat, and even though resistance varieties are available to curb this pest, they are becoming obsolete with the development of new virulent aphid populations. Unlike many other aphids, D noxia only harbours a single endosymbiont, Buchnera aphidicola. Considering the importance of Buchnera, this study aimed to elucidate commonalities and dissimilarities between various hosts, to better understand its distinctiveness within its symbiotic relationship with D. noxia. To do so, the genome of the D. noxia's Buchnera was assembled and compared to those of other aphid species that feed on diverse host species.

RESULTS

The overall importance of several features such as gene length and percentage GC content was found to be critical for the maintenance of Buchnera genes when compared to their closest free-living relative, Escherichia coli. Buchnera protein coding genes were found to have percentage GC contents that tended towards a mean of ~ 26% which had strong correlation to their identity to their E. coli homologs. Several SNPs were identified between different aphid populations and multiple isolates of Buchnera were confirmed in single aphids.

CONCLUSIONS

Establishing the strong correlation of percentage GC content of protein coding genes and gene identity will allow for identifying which genes will be lost in the continually shrinking Buchnera genome. This is also the first report of a parthenogenically reproducing aphid that hosts multiple Buchnera strains in a single aphid, raising questions regarding the benefits of maintaining multiple strains. We also found preliminary evidence for post-transcriptional regulation of Buchnera genes in the form of polyadenylation.

Co-option of a conserved host glutamine transporter facilitates aphid/Buchnera metabolic integration

Organisms across the tree of life colonize novel environments by partnering with bacterial symbionts. These symbioses are characterized by intimate integration of host/endosymbiont biology at multiple levels, including metabolically. Metabolic integration is particularly important for sap-feeding insects and their symbionts, which supplement nutritionally unbalanced host diets. Many studies reveal parallel evolution of host/endosymbiont metabolic complementarity in amino acid biosynthesis, raising questions about how amino acid metabolism is regulated, how regulatory mechanisms evolve, and the extent to which similar mechanisms evolve in different systems. In the aphid/Buchnera symbiosis, the transporter ApGLNT1 (Acyrthosiphon pisum glutamine transporter 1) supplies glutamine, an amino donor in transamination reactions, to bacteriocytes (where Buchnera reside) and is competitively inhibited by Buchnera-supplied arginine-consistent with a role regulating amino acid metabolism given host demand for Buchnera-produced amino acids. We examined how ApGLNT1 evolved a regulatory role by functionally characterizing orthologs in insects with and without endosymbionts. ApGLNT1 orthologs are functionally similar, and orthology searches coupled with homology modeling revealed that GLNT1 is ancient and structurally conserved across insects. Our results indicate that the ApGLNT1 symbiotic regulatory role is derived from its ancestral role and, in aphids, is likely facilitated by loss of arginine biosynthesis through the urea cycle. Given consistent loss of host arginine biosynthesis and retention of endosymbiont arginine supply, we hypothesize that GLNT1 is a general mechanism regulating amino acid metabolism in sap-feeding insects. This work fills a gap, highlighting the broad importance of co-option of ancestral proteins to novel contexts in the evolution of host/symbiont systems.

Proton gradient mediates hemolymph trehalose influx into aphid bacteriocytes

Aphids harbor proteobacterial endosymbionts such as Buchnera aphidicola housed in specialized bacteriocytes derived from host cells. The endosymbiont Buchnera supplies essential amino acids such as arginine to the host cells and, in turn, obtains sugars needed for its survival from the hemolymph. The mechanism of sugar supply in aphid bacteriocytes has been rarely studied. It also remains unclear how Buchnera acquires its carbon source. The hemolymph sugars in Acyrthosiphon pisum are composed of the disaccharide trehalose containing two glucose molecules. Here, we report for the first time that trehalose is transported and used as a potential carbon source by Buchnera across the bacteriocyte plasma membrane via trehalose transporters. The current study characterized the bacteriocyte trehalose transporter Ap_ST11 (LOC100159441) using the Xenopus oocyte expression system. The Ap_ST11 transporter was found to be proton-dependent with a K_m value ≥700 mM. We re-examined the hemolymph trehalose at 217.8 mM using a fluorescent trehalose sensor. The bacteriocytes did not obtain trehalose by facilitated diffusion along the gradient across cellular membranes. These findings suggest that trehalose influx into the bacteriocytes depends on the extracellular proton-driven secondary electrochemical transporter.

Elucidation of host and symbiont contributions to peptidoglycan metabolism based on comparative genomics of eight aphid subfamilies and their Buchnera

Pea aphids (Acyrthosiphon pisum) are insects containing genes of bacterial origin with putative functions in peptidoglycan (PGN) metabolism. Of these, rlpA1-5, amiD, and ldcA are highly expressed in bacteriocytes, specialized aphid cells that harbor the obligate bacterial symbiont Buchnera aphidicola, required for amino acid supplementation of the host’s nutrient-poor diet. Despite genome reduction associated with endosymbiosis, pea aphid Buchnera retains genes for the synthesis of PGN while Buchnera of many other aphid species partially or completely lack these genes. To explore the evolution of aphid horizontally-transferred genes (HTGs) and to elucidate how host and symbiont genes contribute to PGN production, we sequenced genomes from four deeply branching lineages, such that paired aphid and Buchnera genomes are now available for 17 species representing eight subfamilies. We identified all host and symbiont genes putatively involved in PGN metabolism. Phylogenetic analyses indicate that each HTG family was present in the aphid shared ancestor, but that each underwent a unique pattern of gene loss or duplication in descendant lineages. While four aphid rlpA gene subfamilies show no relation to symbiont PGN gene repertoire, the loss of aphid amiD and ldcA HTGs coincides with the loss of symbiont PGN metabolism genes. In particular, the coincident loss of host amiD and symbiont murCEF in tribe Aphidini, in contrast to tribe Macrosiphini, suggests either 1) functional linkage between these host and symbiont genes, or 2) Aphidini has lost functional PGN synthesis and other retained PGN pathway genes are non-functional. To test these hypotheses experimentally, we used cell-wall labeling methods involving a d-alanine probe and found that both Macrosiphini and Aphidini retain Buchnera PGN synthesis. Our results imply that compensatory adaptations can preserve PGN synthesis despite the loss of some genes considered essential for this pathway, highlighting the importance of the cell wall in these symbioses.

Throughout evolution, animals have sometimes gained novel abilities by acquiring bacterial genes through horizontal gene transfer. For some insects harboring bacterial symbionts, horizontally-transferred genes may enable hosts to regulate symbiosis by influencing symbiont cell wall metabolism. While mealybug horizontally-transferred genes work collectively to synthesize the symbiont cell wall, the role of aphid horizontally-transferred genes in symbiont cell wall metabolism is unclear. We examined whether different aphid horizontally-transferred genes co-occur with symbiont genes underlying cell wall metabolism across different aphid lineages, indicative of linked function. We included 17 aphid species representing eight distantly related lineages, four of which we sequenced for this study. We found that two of the three horizontally-acquired gene families are present only when symbionts possess cell wall pathway genes, while the third shows no correlation. These results reveal that despite their putative involvement in symbiont cell wall synthesis, aphid horizontally-acquired genes operate independently from one another and likely have lineage-specific functions. Furthermore, we observed that symbiont cell wall synthesis is maintained in one aphid lineage despite loss of genes considered essential for producing the cell wall, implying that other adaptations preserve the cell wall in aphid species with incomplete cell wall synthesis pathways.

B-vitamin nutrition in the pea aphid-Buchnera symbiosis

Various insects that utilize vitamin-deficient diets derive a supplementary supply of these micronutrients from their symbiotic microorganisms. Here, we tested the inference from genome annotation that the symbiotic bacterium Buchnera aphidicola in the pea aphid Acyrthosiphon pisum provides the insect with vitamins B₂ and B₅ but no other B-vitamins. Contrary to expectation, aphid survival over five days of larval development on artificial diets individually lacking each B-vitamin not synthesized by Buchnera was not significantly reduced, despite significantly lower carcass B₁, B₃, B₆ and B₇ concentrations in the aphids on diets lacking each of these B-vitamins than on the vitamin-complete diet. Aphid survival was, however, significantly reduced on diet containing low concentrations (≤0.2 mM) or no pantothenate (B₅). Complementary transcriptome analysis revealed low abundance of the sense-transcript, but high abundance of the antisense transcript, of the Buchnera gene panC encoding the enzyme mediating the terminal reaction in pantothenate synthesis. We hypothesize that metabolic constraints or antisense transcripts may reduce Buchnera-mediated production of pantothenate, resulting in poor aphid performance on pantothenate-free diets. The discrepancy between predictions from genome data and empirical data illustrates the need for physiological study to test functional inferences made from genome annotations.

Widespread expression of conserved small RNAs in small symbiont genomes

Genome architecture of a microbe markedly changes when it transitions from a free-living lifestyle to an obligate symbiotic association within eukaryotic cells. These symbiont genomes experience numerous rearrangements and massive gene loss, which is expected to radically alter gene regulatory networks compared with those of free-living relatives. As such, it remains unclear whether and how these small symbiont genomes regulate gene expression. Here, using a label-free mass-spec quantification approach we found that differential protein regulation occurs in Buchnera, a model symbiont with a reduced genome, when it transitions between two distinct life stages. However, differential mRNA expression could not be detected between Buchnera life stages, despite the presence of a small number of putative transcriptional regulators. Instead a comparative analysis of small RNA expression profiles among five divergent Buchnera lineages, spanning a variety of Buchnera life stages, reveals 140 novel intergenic and antisense small RNAs and 517 untranslated regions that were significantly expressed, some of which have been conserved for ∼65 million years. In addition, the majority of these small RNAs exhibit both sequence covariation and thermodynamic stability, indicators of a potential structural RNA role. Together, these data suggest that gene regulation at the post-transcriptional level may be important in Buchnera. This is the first study to empirically identify Buchnera small RNAs, and we propose that these novel small RNAs may facilitate post-transcriptional regulation through translational inhibition/activation, and/or transcript stability. Ultimately, post-transcriptional regulation may shape metabolic complementation between Buchnera and its aphid host, thus impacting the animal's ecology and evolution.

The central role of the host cell in symbiotic nitrogen metabolism

Symbiotic nitrogen recycling enables animals to thrive on nitrogen-poor diets and environments. It traditionally refers to the utilization of animal waste nitrogen by symbiotic micro-organisms to synthesize essential amino acids (EAAs), which are translocated back to the animal host. We applied metabolic modelling and complementary metabolite profiling to investigate nitrogen recycling in the symbiosis between the pea aphid and the intracellular bacterium Buchnera, which synthesizes EAAs. The results differ from traditional notions of nitrogen recycling in two important respects. First, aphid waste ammonia is recycled predominantly by the host cell (bacteriocyte) and not Buchnera. Host cell recycling is mediated by shared biosynthetic pathways for four EAAs, in which aphid transaminases incorporate ammonia-derived nitrogen into carbon skeletons synthesized by Buchnera to generate EAAs. Second, the ammonia substrate for nitrogen recycling is derived from bacteriocyte metabolism, such that the symbiosis is not a sink for nitrogenous waste from other aphid organs. Host cell-mediated nitrogen recycling may be general among insect symbioses with shared EAA biosynthetic pathways generated by the loss of symbiont genes mediating terminal reactions in EAA synthesis.

Serratia symbiotica from the aphid Cinara cedri: a missing link from facultative to obligate insect endosymbiont

The genome sequencing of Buchnera aphidicola BCc from the aphid Cinara cedri, which is the smallest known Buchnera genome, revealed that this bacterium had lost its symbiotic role, as it was not able to synthesize tryptophan and riboflavin. Moreover, the biosynthesis of tryptophan is shared with the endosymbiont Serratia symbiotica SCc, which coexists with B. aphidicola in this aphid. The whole-genome sequencing of S. symbiotica SCc reveals an endosymbiont in a stage of genome reduction that is closer to an obligate endosymbiont, such as B. aphidicola from Acyrthosiphon pisum, than to another S. symbiotica, which is a facultative endosymbiont in this aphid, and presents much less gene decay. The comparison between both S. symbiotica enables us to propose an evolutionary scenario of the transition from facultative to obligate endosymbiont. Metabolic inferences of B. aphidicola BCc and S. symbiotica SCc reveal that most of the functions carried out by B. aphidicola in A. pisum are now either conserved in B. aphidicola BCc or taken over by S. symbiotica. In addition, there are several cases of metabolic complementation giving functional stability to the whole consortium and evolutionary preservation of the actors involved.

Large-scale label-free quantitative proteomics of the pea aphid-Buchnera symbiosis

Many insects are nutritionally dependent on symbiotic microorganisms that have tiny genomes and are housed in specialized host cells called bacteriocytes. The obligate symbiosis between the pea aphid Acyrthosiphon pisum and the γ-proteobacterium Buchnera aphidicola (only 584 predicted proteins) is particularly amenable for molecular analysis because the genomes of both partners have been sequenced. To better define the symbiotic relationship between this aphid and Buchnera, we used large-scale, high accuracy tandem mass spectrometry (nanoLC-LTQ-Orbtrap) to identify aphid and Buchnera proteins in the whole aphid body, purified bacteriocytes, isolated Buchnera cells and the residual bacteriocyte fraction. More than 1900 aphid and 400 Buchnera proteins were identified. All enzymes in amino acid metabolism annotated in the Buchnera genome were detected, reflecting the high (68%) coverage of the proteome and supporting the core function of Buchnera in the aphid symbiosis. Transporters mediating the transport of predicted metabolites were present in the bacteriocyte. Label-free spectral counting combined with hierarchical clustering, allowed to define the quantitative distribution of a subset of these proteins across both symbiotic partners, yielding no evidence for the selective transfer of protein among the partners in either direction. This is the first quantitative proteome analysis of bacteriocyte symbiosis, providing a wealth of information about molecular function of both the host cell and bacterial symbiont.

Age- and morph-dependent activation of the lysosomal system and Buchnera degradation in aphid endosymbiosis

Endosymbiosis in aphids is maintained through a mutualistic association between the host and a symbiotic bacterium, Buchnera, which is harbored in specialized host cells called bacteriocytes. Here, we examined the changes in the Buchnera density in bacteriocytes in relation to the development and polyphenism of the host aphid. Buchnera density in the winged morph aphids, alatae, decreased drastically around the final ecdysis, whereas in the wingless morph aphids, apterae, Buchnera density decreased after the final ecdysis. Thereafter, in both apterae and alatae, Buchnera density was maintained at a constant level until 10 days and then again decreased gradually until 18 days after the final ecdysis. Cytochemical analysis with LysoTracker reagent and quantitative RT-PCR analysis revealed that the number of lysosome-like acidic organelles and the amount of lysosome-related gene (lysozyme and cathepsin L) transcripts increased drastically in the bacteriocytes of alatae around the final ecdysis. Electron microscopy of alatae bacteriocytes around the final ecdysis revealed many Buchnera with irregular electron-dense areas in their cytoplasm that were enclosed by a distended symbiosome membrane. These findings indicated that age- and morph-dependent decreases in Buchnera density coincided with activation of the host lysosomal system and the increased degradation of Buchnera.

The striking case of tryptophan provision in the cedar aphid Cinara cedri

Buchnera aphidicola BCc has lost its symbiotic role as the tryptophan supplier to the aphid Cinara cedri. We report the presence of a plasmid in this endosymbiont that contains the trpEG genes. The remaining genes for the pathway (trpDCBA) are located on the chromosome of the secondary endosymbiont "Candidatus Serratia symbiotica." Thus, we propose that a symbiotic consortium is necessary to provide tryptophan.

Biosynthetic potentials from species-specific metabolic networks

Studies of genome-scale metabolic networks allow for qualitative and quantitative descriptions of an organism's capability to convert nutrients into products. The set of synthesizable products strongly depends on the provided nutrients as well as on the structure of the metabolic network. Here, we apply the method of network expansion and the concept of scopes, describing the synthesizing capacities of an organism when certain nutrients are provided. We analyze the biosynthetic properties of four species: Arabidopsis thaliana, Saccharomyces cerevisiae, Buchnera aphidicola, and Escherichia coli. Matthäus et al. have recently developed a method to identify clusters of scopes, reflecting specific biological functions and exhibiting a hierarchical arrangement, using the network comprising all reactions in KEGG. We extend this method by considering random sets of nutrients on well-curated networks of the investigated species from BioCyc. We identify structural properties of the networks that allow to differentiate their biosynthetic capabilities. Furthermore, we evaluate the quality of the clustering of scopes applied to the species-specific networks. Our study provides a novel assessment of the biosynthetic properties of different species.

A small microbial genome: the end of a long symbiotic relationship?

Intracellular bacteria are characterized by genome reduction. The 422,434-base pair genome of Buchnera aphidicola BCc, primary endosymbiont of the aphid Cinara cedri, is approximately 200 kilobases smaller than the previously sequenced B. aphidicola genomes. B. aphidicola BCc has lost most metabolic functions, including the ability to synthesize the essential amino acid tryptophan and riboflavin. In addition, most retained genes are evolving rapidly. Possibly, B. aphidicola BCc is losing its symbiotic capacity and is being complemented (and might be replaced) by the highly abundant coexisting secondary symbiont.

Different levels of transcriptional regulation due to trophic constraints in the reduced genome of Buchnera aphidicola APS

Symbiotic associations involving intracellular microorganisms and animals are widespread, especially for species feeding on poor or unbalanced diets. Buchnera aphidicola, the obligate intracellular bacterium associated with most aphid species, provides its hosts with essential amino acids (EAAs), nutrients in short supply in the plant phloem sap. The Buchnera genome has undergone severe reductions during intracellular evolution. Genes for EAA biosynthesis are conserved, but most of the transcriptional regulatory elements are lost. This work addresses two main questions: is transcription in Buchnera (i) regulated and (ii) scaled to aphid EAA demand? Two microarray experiments were designed for profiling the gene expression in Buchnera. The first one was characterized by a specific depletion of tyrosine and phenylalanine in the aphid diet, and the second experiment combined a global diminution of EAAs in the aphid diet with a sucrose concentration increase to manipulate the aphid growth rate. Aphid biological performance and budget analysis (the balance between EAAs provided by the diet and those synthesized by Buchnera) were performed to quantify the nutritional demand from the aphids toward their symbiotic bacteria. Despite the absence of known regulatory elements, a significant transcriptional regulation was observed at different levels of organization in the Buchnera genome: between genes, within putative transcription units, and within specific metabolic pathways. However, unambiguous evidence for transcriptional changes underpinning the scaling of EAA biosynthesis to aphid demand was not obtained. The phenotypic relevance of the transcriptional response from the reduced genome of Buchnera is addressed.

Codon usage bias and tRNA over-expression in Buchnera aphidicola after aromatic amino acid nutritional stress on its host Acyrthosiphon pisum

Codon usage bias and relative abundances of tRNA isoacceptors were analysed in the obligate intracellular symbiotic bacterium, Buchnera aphidicola from the aphid Acyrthosiphon pisum, using a dedicated 35mer oligonucleotide microarray. Buchnera is archetypal of organisms living with minimal metabolic requirements and presents a reduced genome with high-evolutionary rate. Codonusage in Buchnera has been overcome by the high mutational bias towards AT bases. However, several lines of evidence for codon usage selection are given here. A significant correlation was found between tRNA relative abundances and codon composition of Buchnera genes. A significant codon usage bias was found for the choice of rare codons in Buchnera: C-ending codons are preferred in highly expressed genes, whereas G-ending codons are avoided. This bias is not explained by GC skew in the bacteria and might correspond to a selection for perfect matching between codon-anticodon pairs for some essential amino acids in Buchnera proteins. Nutritional stress applied to the aphid host induced a significant overexpression of most of the tRNA isoacceptors in bacteria. Although, molecular regulation of the tRNA operons in Buchnera was not investigated, a correlation between relative expression levels and organization in transcription unit was found in the genome of Buchnera.

A dual-genome microarray for the pea aphid, Acyrthosiphon pisum, and its obligate bacterial symbiont, Buchnera aphidicola

Background

The best studied insect-symbiont system is that of aphids and their primary bacterial endosymbiont Buchnera aphidicola. Buchnera inhabits specialized host cells called bacteriocytes, provides nutrients to the aphid and has co-speciated with its aphid hosts for the past 150 million years. We have used a single microarray to examine gene expression in the pea aphid, Acyrthosiphon pisum, and its resident Buchnera. Very little is known of gene expression in aphids, few studies have examined gene expression in Buchnera, and no study has examined simultaneously the expression profiles of a host and its symbiont. Expression profiling of aphids, in studies such as this, will be critical for assigning newly discovered A. pisum genes to functional roles. In particular, because aphids possess many genes that are absent from Drosophila and other holometabolous insect taxa, aphid genome annotation efforts cannot rely entirely on homology to the best-studied insect systems. Development of this dual-genome array represents a first attempt to characterize gene expression in this emerging model system.

Results

We chose to examine heat shock response because it has been well characterized both in Buchnera and in other insect species. Our results from the Buchnera of A. pisum show responses for the same gene set as an earlier study of heat shock response in Buchnera for the host aphid Schizaphis graminum. Additionally, analyses of aphid transcripts showed the expected response for homologs of known heat shock genes as well as responses for several genes with unknown functional roles.

Conclusion

We examined gene expression under heat shock of an insect and its bacterial symbiont in a single assay using a dual-genome microarray. Further, our results indicate that microarrays are a useful tool for inferring functional roles of genes in A. pisum and other insects and suggest that the pea aphid genome may contain many gene paralogs that are differentially regulated.

Plasmids in the aphid endosymbiont Buchnera aphidicola with the smallest genomes. A puzzling evolutionary story

Buchnera aphidicola, the primary endosymbiont of aphids, has undergone important genomic and biochemical changes as an adaptation to intracellular life. The most important structural changes include a drastic genome reduction and the amplification of genes encoding key enzymes for the biosynthesis of amino acids by their translocation to plasmids. Molecular characterization through different aphid subfamilies has revealed that the genes involved in leucine and tryptophan biosynthesis show a variable fate, since they can be located on plasmids or on the chromosome in different lineages. This versatility contrasts with the genomic stasis found in three distantly related B. aphidicola strains already sequenced. We present the analysis of three B. aphidicola strains (BTg, BCt and BCc) belonging to aphids from different tribes of the subfamily Lachninae, that was estimated to harbour the bacteria with the smallest genomes. The presence of both leucine and tryptophan plasmids in BTg, a chimerical leucine-tryptophan plasmid in BCt, and only a leucine plasmid in BCc, indicates the existence of many recombination events in a recA minus bacterium. In addition, these B. aphidicola plasmids are the simplest described in this species, indicating that plasmids are also involved in the genome shrinkage process.

Phloem-sap feeding by animals: problems and solutions

The incidence of phloem sap feeding by animals appears paradoxical. Although phloem sap is nutrient-rich compared with many other plant products and generally lacking in toxins and feeding deterrents, it is consumed as the dominant or sole diet by a very restricted range of animals, exclusively insects of the order Hemiptera. These insects display two sets of adaptations. First, linked to the high ratio of non-essential:essential amino acids in phloem sap, these insects contain symbiotic micro-organisms which provide them with essential amino acids. For example, bacteria of the genus Buchnera contribute up to 90% of the essential amino acids required by the pea aphid Acyrthosiphon pisum feeding on Vicia faba. Second, the insect tolerance of the very high sugar content and osmotic pressure of phloem sap is promoted by their possession in the gut of sucrase-transglucosidase activity, which transforms excess ingested sugar into long-chain oligosaccharides voided via honeydew. Various other animals consume phloem sap by proxy, through feeding on the honeydew of phloem-feeding hemipterans. Honeydew is physiologically less extreme than phloem sap, with a higher essential:non-essential amino acid ratio and lower osmotic pressure. Even so, ant species strongly dependent on honeydew as food may benefit from nutrients derived from their symbiotic bacteria Blochmannia.

Metabolic interdependence of obligate intracellular bacteria and their insect hosts

Mutualistic associations of obligate intracellular bacteria and insects have attracted much interest in the past few years due to the evolutionary consequences for their genome structure. However, much less attention has been paid to the metabolic ramifications for these endosymbiotic microorganisms, which have to compete with but also to adapt to another metabolism--that of the host cell. This review attempts to provide insights into the complex physiological interactions and the evolution of metabolic pathways of several mutualistic bacteria of aphids, ants, and tsetse flies and their insect hosts.

Mutational and selective pressures on codon and amino acid usage in Buchnera, endosymbiotic bacteria of aphids

We have explored compositional variation at synonymous (codon usage) and nonsynonymous (amino acid usage) positions in three complete genomes of Buchnera, endosymbiotic bacteria of aphids, and also in their orthologs in Escherichia coli, a close free-living relative. We sought to discriminate genes of variable expression levels in order to weigh the relative contributions of mutational bias and selection in the genomic changes following symbiosis. We identified clear strand asymmetries, distribution biases (putative high-expression genes were found more often on the leading strand), and a residual slight codon bias within each strand. Amino acid usage was strongly biased in putative high-expression genes, characterized by avoidance of aromatic amino acids, but above all by greater conservation and resistance to AT enrichment. Despite the almost complete loss of codon bias and heavy mutational pressure, selective forces are still strong at nonsynonymous sites of a fraction of the genome. However, Buchnera from Baizongia pistaciae appears to have suffered a stronger symbiotic syndrome than the two other species.

Evolution of minimal-gene-sets in host-dependent bacteria

Several attempts have been made to identify the minimal set of genes that is required for life using computational approaches or studies of deletion mutants. These experiments resemble those already performed by nature; a few hundred million years ago an ancestor of Escherichia coli was domesticated by aphids, which resulted in the elimination of 70-75% of the original bacterial genome. Amazingly, the small genomes of these imprisoned bacteria are more stable than those of their free-living relatives. Minimal-gene-sets that have evolved naturally are largely species-specific, with the exception of a small set of core genes that are required for information processing. Comparative genomics of host-dependent bacteria have shown that minimal-gene-sets can persist in nature for tens of millions of years provided that the environment is rich in nutrients, that the host population size is large and that there is a strong host-level selection for bacterial gene functions.

Developmental origin and evolution of bacteriocytes in the aphid-Buchnera symbiosis

Symbiotic relationships between bacteria and insect hosts are common. Although the bacterial endosymbionts have been subjected to intense investigation, little is known of the host cells in which they reside, the bacteriocytes. We have studied the development and evolution of aphid bacteriocytes, the host cells that contain the endosymbiotic bacteria Buchnera aphidicola. We show that bacteriocytes of Acyrthosiphon pisum express several gene products (or their paralogues): Distal-less, Ultrabithorax/Abdominal-A, and Engrailed. Using these markers, we find that a subpopulation of the bacteriocytes is specified prior to the transmission of maternal bacteria to the embryo. In addition, we discovered that a second population of cells is recruited to the bacteriocyte fate later in development. We experimentally demonstrate that bacteriocyte induction and proliferation occur independently of B. aphidicola. Major features of bacteriocyte development, including the two-step recruitment of bacteriocytes, have been conserved in aphids for 80–150 million years. Furthermore, we have investigated two cases of evolutionary loss of bacterial symbionts: in one case, where novel extracellular, eukaryotic symbionts replaced the bacteria, the bacteriocyte is maintained; in another case, where symbionts are absent, the bacteriocytes are initiated but not maintained. The bacteriocyte represents an evolutionarily novel cell fate, which is developmentally determined independently of the bacteria. Three of five transcription factors we examined show novel expression patterns in bacteriocytes, suggesting that bacteriocytes may have evolved to express many additional transcription factors. The evolutionary transition to a symbiosis in which bacteria and an aphid cell form a functional unit, similar to the origin of plastids, has apparently involved extensive molecular adaptations on the part of the host cell.

Molecular markers show that bacteriocytes, the aphid cells that house the bacterial endosymbionts, are specified in a conserved two-step process that does not depend on the presence of the bacteria

The evolution of the heat-shock protein GroEL from Buchnera, the primary endosymbiont of aphids, is governed by positive selection

The heat-shock protein GroEL is a double-ring-structured chaperonin that assists the folding of many newly synthesized proteins in Escherichia coli and the refolding in vitro, with the cochaperonin GroES, of conformationally damaged proteins. This protein is constitutively overexpressed in the primary symbiotic bacteria of many insects, constituting approximately 10% of the total protein in Buchnera, the primary endosymbiont of aphids. In the present study, we perform a maximum likelihood (ML) analysis to unveil the selective constraints in GroEL. In addition, we apply a new statistical approach to determine the patterns of evolution in this highly interesting protein. The main conclusion derived from our analysis is that GroEL has suffered an accelerated rate of amino acid substitution upon the symbiotic integration of Buchnera into the aphids. It is most interesting that the ML analysis of codon substitutions in the different branches of the phylogenetic tree strongly supports the action of positive selection in the different lineages of BUCHNERA: Additionally, the new sliding window analysis of the complete groEL sequence reveals different regions of the molecule under the action of positive selection, mainly located in the apical domain, that are important for both peptide and GroES binding.

Reduction of bacterial genome size and expansion resulting from obligate intracellular lifestyle and adaptation to soil habitat

Prokaryotic organisms are exposed in the course of evolution to various impacts, resulting often in drastic changes of their genome size. Depending on circumstances, the same lineage may diverge into species having substantially reduced genomes, or such whose genomes have undergone considerable enlargement. Genome reduction is a consequence of obligate intracellular lifestyle rendering numerous genes expendable. Another consequence of intracellular lifestyle is reduction of effective population size and limited possibility of gene acquirement via lateral transfer. This causes a state of relaxed selection resulting in accumulation of mildly deleterious mutations that can not be corrected by recombination with the wild type copy. Thus, gene loss is usually irreversible. Additionally, constant environment of the eukaryotic cell renders that some bacterial genes involved in DNA repair are expandable. The loss of these genes is a probable cause of mutational bias resulting in a high A+T content. While causes of genome reduction are rather indisputable, those resulting in genome expansion seem to be less obvious. Presumably, the genome enlargement is an indirect consequence of adaptation to changing environmental conditions and requires the acquisition and integration of numerous genes. It seems that the need for a great number of capabilities is common among soil bacteria irrespective of their phylogenetic relationship. However, this would not be possible if soil bacteria lacked indigenous abilities to exchange and accumulate genetic information. The latter are considerably facilitated when housekeeping genes are physically separated from adaptive loci which are useful only in certain circumstances.

Molecular characterization of the leucine cluster in Buchnera sp. strain PSY, a primary endosymbiont of the aphid Pemphigus spyrothecae

Buchnera strains from most aphid subfamilies studied to date have been found to carry the leucine gene cluster (leuA, -B, -C, and -D) on a plasmid, an organization unique among bacteria. Here, however, we demonstrate a classical chromosomal location of the cluster in Buchnera sp. strain PSY from the aphid Pemphigus spyrothecae (subfamily Pemphiginae). The genes that flank leuABCD in Buchnera sp. strain PSY appear to be adjacent in the genome of Buchnera sp. strain APS, a strain carrying a leucine plasmid. We propose that the presence of a leucine plasmid predates the diversification of symbiotic Buchnera and that the chromosomal location observed in Buchnera sp. strain PSY arose by a transfer of the leucine genes from a plasmid to the chromosome.

Compositional correlation and codon usage studies in Buchnera aphidicola

Compositional distributions in three different codon positions as well as codon usage biases of all available DNA sequences of Buchnera aphidicola genome have been analyzed. It was observed that GC levels among the three codon positions is I>II>III as observed in other extremely high AT rich organisms. B. aphidicola being an AT rich organism is expected to have A and/or T at the third positions of codons. Overall codon usage analyses indicate that A and/or T ending codons are predominant in this organism and some particular amino acids are abundant in the coding region of genes. However, multivariate statistical analysis indicates two major trends in the codon usage variation among the genes; one being strongly correlated with the GC contents at the third synonymous positions of codons, and the other being associated with the expression level of genes. Moreover, codon usage biases of the highly expressed genes are almost identical with the overall codon usage biases of all the genes of this organism. These observations suggest that mutational bias is the main factor in determining the codon usage variation among the genes in B. aphidicola.

The process of genome shrinkage in the obligate symbiont Buchnera aphidicola

BACKGROUND

Very small genomes have evolved repeatedly in eubacterial lineages that have adopted obligate associations with eukaryotic hosts. Complete genome sequences have revealed that small genomes retain very different gene sets, raising the question of how final genome content is determined. To examine the process of genome reduction, the tiny genome of the endosymbiont Buchnera aphidicola was compared to the larger ancestral genome, reconstructed on the basis of the phylogenetic distribution of gene orthologs among fully sequenced relatives of Escherichia coli and Buchnera.

RESULTS

The reconstructed ancestral genome contained 2,425 open reading frames (ORFs). The Buchnera genome, containing 564 ORFs, consists of 153 fragments of 1-34 genes that are syntenic with reconstructed ancestral regions. On the basis of this reconstruction, 503 genes were eliminated within syntenic fragments, and 1,403 genes were lost from the gaps between syntenic fragments, probably in connection with genome rearrangements. Lost regions are sometimes large, and often span functionally unrelated genes. In addition, individual genes and regulatory regions have been lost or eroded. For the categories of DNA repair genes and rRNA genes, most lost loci fall in regions between syntenic fragments. This history of gene loss is reflected in the sequences of intergenic spacers at positions where genes were once present.

CONCLUSIONS

The most plausible interpretation of this reconstruction is that Buchnera lost many genes through the fixation of large deletions soon after the acquisition of an obligate endosymbiotic lifestyle. An implication is that final genome composition may be partly the chance outcome of initial deletions and that neighboring genes influence the likelihood of loss of particular genes and pathways.

The phylogenetic distribution of frataxin indicates a role in iron-sulfur cluster protein assembly

Much has been learned about the cellular pathology of Friedreich's ataxia, a recessive neurodegenerative disease resulting from insufficient expression of the mitochondrial protein frataxin. However, the biochemical function of frataxin has remained obscure, hampering attempts at therapeutic intervention. To predict functional interactions of frataxin with other proteins we investigated whether its gene specifically co-occurs with any other genes in sequenced genomes. In 56 available genomes we identified two genes with identical phylogenetic distributions to the frataxin/cyaY gene: hscA and hscB/JAC1. These genes have not only emerged in the same evolutionary lineage as the frataxin gene, they have also been lost at least twice with it, and they have been horizontally transferred with it in the evolution of the mitochondria. The proteins encoded by hscA and hscB, the chaperone HSP66 and the co-chaperone HSP20, have been shown to be required for the synthesis of 2Fe-2S clusters on ferredoxin in proteobacteria. JAC1, an ortholog of hscB, and SSQ1, a paralog of hscA, have been shown to be required for iron-sulfur cluster assembly in mitochondria of Saccharomyces cerevisiae. Combining data on the co-occurrence of genes in genomes with experimental and predicted cellular localization data of their proteins supports the hypothesis that frataxin is directly involved in iron-sulfur cluster protein assembly. They indicate that frataxin is specifically involved in the same sub-process as HSP20/Jac1p.

The impact of host plant on the abundance and function of symbiotic bacteria in an aphid

The black-bean aphid Aphis fabae bears populations of coccoid symbiotic bacteria Buchnera spp. at 2.0–3.2×107cellsmg−1aphidmass and rod-shaped secondary symbionts of uncertain taxonomic affiliation at 0.1–0.6×107cellsmg−1aphidmass. Buchnera provides essential amino acids, supplementing the poor supply in the aphid diet of plant phloem sap. Comparison of the performance of A. fabae containing and experimentally deprived of their bacteria showed that the bacteria caused increased larval mass of aphids reared on Chenopodium album and Papaver dubium plants, but not when reared on Lamium purpureum. In the aphids reared on L. purpureum, the density of the bacteria, especially the secondary symbionts, was significantly elevated, and bacterial-mediated production of the essential amino acid threonine was reduced, even though the essential amino acid content of phloem exudates from L. purpureum had a low threonine content. It is proposed that the shortfall in threonine, possibly compounded by the high density of secondary symbionts, may contribute to the poor performance of the aphids on L. purpureum. This study offers the first evidence to suggest plant-mediated interference with the nutritional function of symbiotic bacteria in any phytophagous insect.

Parallel acceleration of evolutionary rates in symbiont genes underlying host nutrition

The overproduction of essential amino acids by Buchnera aphidicola, the primary bacterial mutualist of aphids, is considered an adaptation for increased production of nutrients that are lacking in aphids' diet of plant sap. Given their shared role in host nutrition, amino acid biosynthetic genes of Buchnera are expected to experience parallel changes in selection that depend on host diet quality, growth rate, and population structure. This study evaluates the hypothesis of parallel selection across biosynthetic pathways by testing for correlated changes in evolutionary rates at biosynthetic genes of Buchnera. Previous studies show fast evolutionary rates at tryptophan biosynthetic genes among Buchnera associated with the aphid genus Uroleucon and suggest reduced purifying selection on symbiont nutritional functions in this aphid group. Here, we test for parallel rate acceleration at other amino acid biosynthetic genes of Buchnera-Uroleucon, including those for leucine (leuABC) and isoleucine/valine biosynthesis (ilvC). Ratios of nonsynonymous to synonymous substitutions (d(N)/d(S)) were estimated using codon-based maximum-likelihood methods that account for the extreme AT compositional bias of Buchnera sequences. A significant elevation in d(N)/d(S) at biosynthetic loci but not at two housekeeping genes sampled (dnaN and tuf) suggests reduced host-level selection on biosynthetic capabilities of Buchnera-Uroleucon. In addition, the discovery of trpEG pseudogenes in Buchnera-U. obscurum further supports reduced selection on amino acid biosynthesis.

Expression of host S-adenosylmethionine decarboxylase gene and polyamine composition in aphid bacteriocytes

Differential cDNA display and quantitative RT-PCR revealed that mRNA of host S-adenosylmethionine decarboxylase (SAMDC) was abundant only in the aphid endosymbiotic system well organized in young hosts, suggesting that SAMDC plays some important roles in the system. SAMDC is a key enzyme to synthesize polyamines that are known to be involved in a large array of biological events including protein synthesis, DNA stabilization, DNA replication, and cell proliferation. As the first step to investigate roles of polyamines in the endosymbiotic system, polyamine composition in bacteriocytes was determined by high performance liquid chromatography. As a result, we found that bacteriocytes contained virtually an only single polyamine, spermidine. The spermidine content of bacteriocytes fluctuated with time in the course of development and aging of the host aphid. This is the first report of polyamine assessment in a prokaryote-eukaryote endocellular symbiotic system, which demonstrated a unique polyamine composition.

Quantifying nutrient production by the microbial symbionts in an aphid

The symbiotic bacteria Buchnera sp. provide aphids with essential amino acids, nutrients in short supply in the aphid diet of plant phloem sap. The contribution of Buchnera-derived amino acids to net protein growth of the aphid Aphis fabae was quantified from the protein growth of aphids reared on chemically defined diets lacking individual amino acids. The amino acid production rates varied among the nine essential amino acids over the range 8–156 pmol microg(−1)protein day(−1) (for tryptophan and leucine, respectively), equivalent to 0.02-0.33 fmol Buchnera(−1)day(−1). In a complementary metabolic analysis, the aphids incorporated radioactivity from dietary [(14)C]glutamic acid into the essential amino acids isoleucine, lysine and threonine. Incorporation into isoleucine was significantly elevated by the omission of dietary isoleucine, indicating that dietary supply may affect the biosynthetic rates of certain amino acids by Buchnera. Aphids experimentally deprived of Buchnera did not synthesize essential amino acids from dietary glutamic acid. The mortality of aposymbionts was high over 7 days on the phenylalanine-free diet, and their assimilation of dietary leucine was depressed on the complete diet, suggesting that both the absence of bacteria-derived amino acids and the low rates of assimilation of certain dietary amino acids may contribute to the poor growth of these insects.

Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS

Almost all aphid species (Homoptera, Insecta) have 60-80 huge cells called bacteriocytes, within which are round-shaped bacteria that are designated Buchnera. These bacteria are maternally transmitted to eggs and embryos through host generations, and the mutualism between the host and the bacteria is so obligate that neither can reproduce independently. Buchnera is a close relative of Escherichia coli, but it contains more than 100 genomic copies per cell, and its genome size is only a seventh of that of E. coli. Here we report the complete genome sequence of Buchnera sp. strain APS, which is composed of one 640,681-base-pair chromosome and two small plasmids. There are genes for the biosyntheses of amino acids essential for the hosts in the genome, but those for non-essential amino acids are missing, indicating complementarity and syntrophy between the host and the symbiont. In addition, Buchnera lacks genes for the biosynthesis of cell-surface components, including lipopolysaccharides and phospholipids, regulator genes and genes involved in defence of the cell. These results indicate that Buchnera is completely symbiotic and viable only in its limited niche, the bacteriocyte.

Polyamine composition and expression of genes related to polyamine biosynthesis in an aphid endosymbiont, Buchnera

Polyamine composition in an aphid endosymbiotic bacterium, Buchnera sp., was determined by high-performance liquid chromatographic analysis. We found that Buchnera contained virtually only a single polyamine, spermidine. The spermidine content of Buchnera was considerably higher in young aphids and tended to decrease with the age of the host. Expression of speD and speE, whose gene products are key enzymes in the synthesis of spermidine, was analyzed by real-time quantitative reverse transcription-PCR. It was shown that the levels of their mRNAs fluctuated in line with the spermidine content.

Identifying the determinants in the equatorial domain of Buchnera GroEL implicated in binding Potato leafroll virus

Luteoviruses avoid degradation in the hemolymph of their aphid vector by interacting with a GroEL homolog from the aphid's primary endosymbiotic bacterium (Buchnera sp.). Mutational analysis of GroEL from the primary endosymbiont of Myzus persicae (MpB GroEL) revealed that the amino acids mediating binding of Potato leafroll virus (PLRV; Luteoviridae) are located within residues 9 to 19 and 427 to 457 of the N-terminal and C-terminal regions, respectively, of the discontinuous equatorial domain. Virus overlay assays with a series of overlapping synthetic decameric peptides and their derivatives demonstrated that R13, K15, L17, and R18 of the N-terminal region and R441 and R445 of the C-terminal region of the equatorial domain of GroEL are critical for PLRV binding. Replacement of R441 and R445 by alanine in full-length MpB GroEL and in MpB GroEL deletion mutants reduced but did not abolish PLRV binding. Alanine substitution of either R13 or K15 eliminated the PLRV-binding capacity of the other and those of L17 and R18. In the predicted tertiary structure of GroEL, the determinants mediating virus binding are juxtaposed in the equatorial plain.