Host-based divergence in populations of the pea aphid: insights from nuclear markers and the prevalence of facultative symbionts

In North America, the pea aphid Acyrthosiphon pisum encompasses ecologically and genetically distinct host races that offer an ideal biological system for studies on sympatric speciation. In addition to its obligate symbiont Buchnera, pea aphids harbour several facultative and phylogenetically distant symbionts. We explored the relationships between host races of A. pisum and their symbiotic microbiota to gain insights into the historical process of ecological specialization and symbiotic acquisition in this aphid. We used allozyme and microsatellite markers to analyse the extent of genetic differentiation between populations of A. pisum on pea, alfalfa and clover in France. In parallel, we examined: (i) the distribution of four facultative symbionts; and (ii) the genetic variation in the Buchnera genome across host-associated populations of A. pisum. Our study clearly demonstrates that populations of A. pisum on pea, clover and alfalfa in France are genetically divergent, which indicates that they constitute distinct host races. We also found a very strong association between host races of A. pisum and their symbiotic microbiota. We stress the need for phylogeographic studies to shed light on the process of host-race formation and acquisition of facultative symbionts in A. pisum. We also question the effects of these symbionts on aphid host fitness, including their role in adaptation to a host plant.

[Identification of the horizontal gene transfer based on phylogenetic data]

We suggest a new procedure to search for the genes with horizontal transfer events in their evolutionary history. The search is based on analysis of topology difference between the phylogenetic trees of gene (protein) groups and the corresponding phylogenetic species trees. Numeric values are introduced to measure the discrepancy between the trees. This approach was applied to analyze 40 prokaryotic genomes classified into 132 classes of orthologs. This resulted in a list of the candidate genes for which the hypothesis of horizontal transfer in evolution looks true.

Consequences of reductive evolution for gene expression in an obligate endosymbiont

The smallest cellular genomes are found in obligate symbiotic and pathogenic bacteria living within eukaryotic hosts. In comparison with large genomes of free-living relatives, these reduced genomes are rearranged and have lost most regulatory elements. To test whether reduced bacterial genomes incur reduced regulatory capacities, we used full-genome microarrays to evaluate transcriptional response to environmental stress in Buchnera aphidicola, the obligate endosymbiont of aphids. The 580 genes of the B. aphidicola genome represent a subset of the 4500 genes known from the related organism, Escherichia coli. Although over 20 orthologues of E. coli heat stress (HS) genes are retained by B. aphidicola, only five were differentially expressed after near-lethal heat stress treatments, and only modest shifts were observed. Analyses of upstream regulatory regions revealed loss or degradation of most HS (sigma32) promoters. Genomic rearrangements downstream of an intact HS promoter yielded upregulation of a functionally unrelated and an inactivated gene. Reanalyses of comparable experimental array data for E. coli and Bacillus subtilis revealed that genome-wide differential expression was significantly lower in B. aphidicola. Our demonstration of a diminished stress response validates reports of temperature sensitivity in B. aphidicola and suggests that this reduced bacterial genome exhibits transcriptional inflexibility.

Low and homogeneous copy number of plasmid-borne symbiont genes affecting host nutrition in Buchnera aphidicola of the aphid Uroleucon ambrosiae

The bacterial endosymbiont of aphids, Buchnera aphidicola, often provides amino acids to its hosts. Plasmid amplification of leucine (leuABCD) and tryptophan (trpEG) biosynthesis genes may be a mechanism by which some Buchnera over-produce these nutrients. We used quantitative polymerase chain reaction to assess the leuABCD/trpEG copy variability within Uroleucon ambrosiae, an aphid with a wide diet breadth and range. Both leuABCD and trpEG abundances are: (i) similar for aphids across 15 populations, and (ii) low compared to Buchnera from other aphid species (particularly trpEG). Consequently, the plasmid location of trpEG combined with Buchnera's chromosomal polyploidy may functionally limit, rather than increase, tryptophan production within Uroleucon ambrosiae.

Why are the genomes of endosymbiotic bacteria so stable?

The comparative analysis of three strains of the endosymbiotic bacterium Buchnera aphidicola has revealed high genome stability associated with an almost complete absence of chromosomal rearrangements and horizontal gene transfer events during the past 150 million years. The loss of genes involved in DNA uptake and recombination in the initial stages of endosymbiosis probably underlies this stability. Gene loss, which was extensive during the initial steps of Buchnera evolution, has continued in the different Buchnera lineages since their divergence.

Extraordinary proliferation of microorganisms in aposymbiotic pea aphids, Acyrthosiphon pisum

Aposymbiotic pea aphids, which were deprived of their intracellular symbiotic bacterium, Buchnera, exhibit growth retardation and no fecundity. High performance liquid chromatographic (HPLC) analysis revealed that these aposymbiotic aphids, when reared on broad bean plants, accumulated a large amount of histamine. To assess the possibility of extraordinary proliferation of microorganisms other than Buchnera, we enumerated eubacteria and fungi in aphids using the real-time quantitative PCR method that targets genes encoding small-subunit rRNAs. The result showed that these microorganisms were extremely abundant in the aposymbiotic aphids reared on plants. Microbial communities in aposymbiotic aphids were further profiled by phylogenetic analysis of small-subunit rDNAs. Of 172 nonchimeric sequences of fungal 18S rDNAs, 138 (80.2%) belonged to the phylum Ascomycota. Among them, 21 clustered within a monophyletic group consisting of insect-pathogenic fungi and yeast-like symbionts of homopteran insects. Thirty-one (18.0%), two (1.2%), and one (0.6%) clones were clustered within the Basidiomycota, Zygomycota, and Oomycota, respectively. Of 167 nonchimeric sequences of eubacterial 16S rDNAs, 84 (50.3%) belonged to the gamma-subdivision of Proteobacteria to which most primary endosymbionts of insects and prolific histamine producers belong. Forty (24.0%), 25 (15.0%), 10 (6.0%), and five (3.0%) clones were clustered within alpha-Proteobacteria, Cytophaga-Flavobacterium-Bacteroides (CFB) group, Actinobacteria, and beta-Proteobacteria, respectively. Three had no phylogenetic association with known taxonomic divisions. None of the sequences studied in this study coincided exactly with those deposited in GenBank.

Reductive genome evolution in Buchnera aphidicola

We have sequenced the genome of the intracellular symbiont Buchnera aphidicola from the aphid Baizongia pistacea. This strain diverged 80-150 million years ago from the common ancestor of two previously sequenced Buchnera strains. Here, a field-collected, nonclonal sample of insects was used as source material for laboratory procedures. As a consequence, the genome assembly unveiled intrapopulational variation, consisting of approximately 1,200 polymorphic sites. Comparison of the 618-kb (kbp) genome with the two other Buchnera genomes revealed a nearly perfect gene-order conservation, indicating that the onset of genomic stasis coincided closely with establishment of the symbiosis with aphids, approximately 200 million years ago. Extensive genome reduction also predates the synchronous diversification of Buchnera and its host; but, at a slower rate, gene loss continues among the extant lineages. A computational study of protein folding predicts that proteins in Buchnera, as well as proteins of other intracellular bacteria, are generally characterized by smaller folding efficiency compared with proteins of free living bacteria. These and other degenerative genomic features are discussed in light of compensatory processes and theoretical predictions on the long-term evolutionary fate of symbionts like Buchnera.

Extremely low levels of genetic polymorphism in endosymbionts (Buchnera) of aphids (Pemphigus)

Molecular evolutionary studies have suggested that vertically transmitted endosymbionts are subject to accumulation of deleterious mutations through genetic drift. Predictions of this hypothesis for patterns of intraspecific polymorphism were borne out in the single relevant study available, on the symbiont Buchnera aphidicola of Uroleucon ambrosiae. In order to examine the generality of this result, we surveyed DNA sequence variation in Buchnera of the distantly related aphid, Pemphigus obesinymphae. In contrast to Uroleucon species, Pemphigus species have complex life cycles with few dispersal stages. Despite these differences, P. obesinymphae showed patterns of variation at two Buchnera loci and one mitochondrial locus that were remarkably similar to those reported previously for Buchnera of U. ambrosiae. In the western US, Buchnera was nearly monomorphic, and in the eastern US, synonymous divergence ranged from 0.08 to 0.16%. Most polymorphisms involved rare alleles, consistent with a recent range of ancestral polymorphism, probably due to demographic fluctuations in aphid populations. These results support the generality of small effective population size in Buchnera and their aphid hosts.

Acceleration of genomic evolution caused by enhanced mutation rate in endocellular symbionts

Endosymbionts, which are widely observed in nature, have undergone reductive genome evolution because of their long-term intracellular lifestyle. Here we compared the complete genome sequences of two different endosymbionts, Buchnera and a protist mitochondrion, with their close relatives to study the evolutionary rates of functional genes in endosymbionts. The results indicate that the rate of amino acid substitution is two times higher in symbionts than in their relatives. This rate increase was observed uniformly among different functional classes of genes, although strong purifying selection may have counterbalanced the rate increase in a few cases. Our data suggest that, contrary to current views, neither the Muller's ratchet effect nor the slightly deleterious mutation theory sufficiently accounts for the elevated evolutionary rate. Rather, the elevated evolutionary rate appears to be mainly due to enhanced mutation rate, although the possibility of relaxation of purifying selection cannot be ruled out.

A strong effect of AT mutational bias on amino acid usage in Buchnera is mitigated at high-expression genes

The advent of full genome sequences provides exceptionally rich data sets to explore molecular and evolutionary mechanisms that shape divergence among and within genomes. In this study, we use multivariate analysis to determine the processes driving genome-wide patterns of amino usage in the obligate endosymbiont Buchnera and its close free-living relative Escherichia coli. In the AT-rich Buchnera genome, the primary source of variation in amino acid usage differentiates high- and low-expression genes. Amino acids of high-expression Buchnera genes are generally less aromatic and use relatively GC-rich codons, suggesting that selection against aromatic amino acids and against amino acids with AT-rich codons is stronger in high-expression genes. Selection to maintain hydrophobic amino acids in integral membrane proteins is a primary factor driving protein evolution in E. coli but is a secondary factor in Buchnera. In E. coli, gene expression is a secondary force driving amino acid usage, and a correlation with tRNA abundance suggests that translational selection contributes to this effect. Although this and previous studies demonstrate that AT mutational bias and genetic drift influence amino acid usage in Buchnera, this genome-wide analysis argues that selection is sufficient to affect the amino acid content of proteins with different expression and hydropathy levels.

Comparative molecular evolution of primary (Buchnera) and secondary symbionts of aphids based on two protein-coding genes

A+T content, phylogenetic relationships, codon usage, evolutionary rates, and ratio of synonymous versus non-synonymous substitutions have been studied in partial sequences of the atpD and aroQ/pheA genes of primary ( Buchnera) and secondary symbionts of aphids and a set of selected non-symbiotic bacteria, belonging to the five subdivisions of the Proteobacteria. Compared to the homologous genes of the last group, both genes belonging to Buchnera behave in a similar way, showing a higher A+T content, forming a monophyletic group, a loss in codon bias, especially in third base position, an evolutionary acceleration and an increase in the number of non-synonymous substitutions, confirming previous results reported elsewhere for other genes. When available, these properties have been partly observed with the secondary symbionts, but with values that are intermediate between Buchnera and free living Proteobacteria. They show high A+T content, but not as high as Buchnera, a non-solved phylogenetic position between Buchnera, and the other gamma-Proteobacteria, a loss in codon bias, again not as high as in Buchnera and a significant evolutionary acceleration in the case of the three atpD genes, but not when considering aroQ/pheA genes. These results give support to the hypothesis that they are symbionts at different stages of the symbiotic accommodation to the host.

50 million years of genomic stasis in endosymbiotic bacteria

Comparison of two fully sequenced genomes of Buchnera aphidicola, the obligate endosymbionts of aphids, reveals the most extreme genome stability to date: no chromosome rearrangements or gene acquisitions have occurred in the past 50 to 70 million years, despite substantial sequence evolution and the inactivation and loss of individual genes. In contrast, the genomes of their closest free-living relatives, Escherichia coli and Salmonella spp., are more than 2000-fold more labile in content and gene order. The genomic stasis of B. aphidicola, likely attributable to the loss of phages, repeated sequences, and recA, indicates that B. aphidicola is no longer a source of ecological innovation for its hosts.

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.

Extreme genome reduction in Buchnera spp.: toward the minimal genome needed for symbiotic life

Buchnera is a mutualistic intracellular symbiont of aphids. Their association began about 200 million years ago, with host and symbiont lineages evolving in parallel since that time. During this coevolutionary process, Buchnera has experienced a dramatic decrease of genome size, retaining only essential genes for its specialized lifestyle. Previous studies reported that genome size in Buchnera spp. is very uniform, suggesting that genome shrinkage occurred early in evolution, and that modern lineages retain the genome size of a common ancestor. Our physical mapping of Buchnera genomes obtained from five aphid lineages shows that the genome size is not conserved among them, but has been reduced down to 450 kb in some species. Here we show evidence of six species with a genome size smaller than Mycoplasma genitalium, the smallest bacterial genome reported thus far (580 kb). Our findings strongly suggest that the Buchnera genome is still experiencing a reductive process toward a minimum set of genes necessary for its symbiotic lifestyle.

An experimental validation of orphan genes of Buchnera, a symbiont of aphids

Although Buchnera sp. APS, an intracellular symbiont of pea aphids, is a close relative of Escherichia coli, its genome has been extensively modified because of its prolonged intracellular life. In our previous studies on the Buchnera genome, computer analysis predicted three "orphan" genes, yba2, yba3, and yba4, which are open reading frames (ORFs) with no homologs in the database. In this paper, we successfully validated all these orphan genes by RT-PCR and Northern hybridization. The present study also revealed that yba3 and yba4 formed an operon, suggesting that they function in concert. Sequences around transcriptional start sites suggests that these genes are under the control of sigma 70. In view of codon usage and AT bias observed in these genes, it is likely that Buchnera have maintained them for an evolutionarily long time.

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.

Genome size reduction through multiple events of gene disintegration in Buchnera APS

The evolution of the endosymbiont Buchnera during its adaptation to intracellular life involved a massive reduction in its genome. By comparing the orthologous genes of Buchnera, Escherichia coli and Vibrio cholerae, we show that the minimal genome size of Buchnera arose from multiple events of gene disintegration dispersed over the whole genome. The elimination of the genes was a continuous process that began with gene inactivation and progressed until the DNA corresponding to the pseudogenes were completely deleted.

Accumulation of species-specific amino acid replacements that cause loss of particular protein functions in Buchnera, an endocellular bacterial symbiont

Endosymbiotic bacteria live in animal cells and are transmitted vertically at the time of the host's reproduction. In view of their small and asexual populations with infrequent chances of recombination, these endocellular bacteria are expected to accumulate mildly deleterious mutations. Previous studies showed that the DNA sequences of these bacteria evolved faster than those of free-living bacteria. In this study, we compared all the ORFs of Buchnera, an endocellular bacterial symbiont of aphids, with those of 34 other prokaryotic organisms and estimated the effect of the accelerated evolution of Buchnera on the functions of its proteins. It was revealed that Buchnera proteins contain many mutations at the sites where sequences are conserved in their orthologues in many other organisms. In addition, amino acid replacements at the conserved sites are mostly changes to physicochemically different amino acids. These results suggest that functions and conformations of Buchnera proteins have been seriously impaired or strongly modified. Indeed, extensive loss of functional motifs was observed in some Buchnera proteins. In many Buchnera proteins mutations were not detected evenly throughout each molecule but tended to accumulate in some functional units, possibly leading to loss of specific functions. As Buchnera has an unusual and limited gene repertory, it is conceivable that the manner of interactions among its proteins has been changed, and thus, functional constraints over their amino acid residues have also been changed during evolution. This may account for the loss of some functional units only in the Buchnera proteins. We obtained evidence that amino acid replacements in Buchnera were not always deleterious, but neutral or, in some cases, even positively selected.

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.

Molecular systematics of aphids and their primary endosymbionts

Aphids constitute a monophyletic group within the order Homoptera (i.e., superfamily Aphidoidea). The Aphidoidea originated in the Jurassic about 150 my ago from some aphidiform ancestor whose origin can be traced back to about 250 my ago. They exhibit a mutualistic association with intracellular bacteria (Buchnera sp.) related to Escherichia coli. Buchnera is usually considered the aphids' primary endosymbiont. The association is obligate for both partners. The 16S rDNA-based phylogeny of Buchnera from four aphid families showed complete concordance with the morphology-based phylogeny of their aphid hosts, which pointed to a single original infection in a common ancestor of aphids some 100-250 my ago followed by cospeciation of aphids and Buchnera. This study concentrated on the molecular phylogeny of both the aphids and their primary endosymbionts of five aphid families including for the first time representatives of the family Lachnidae. We discuss results based on two Buchnera genes (16S rDNA and the beta subunit of the F-ATPase complex) and on one host mitochondrial gene (the subunit 6 of the F-ATPase complex). Although our data do not allow definitive evolutionary relationships to be established among the different aphid families, some traditionally accepted groupings are put into question from both bacterial and insect data. In particular, the Lachnidae and the Aphididae, which from morphological data are considered recently evolved sister groups, do not seem to be as closely related as is usually accepted. Finally, we discuss our results in the light of the proposed parallel evolution of aphids and their endosymbionts.

Mutualists and parasites: how to paint yourself into a (metabolic) corner

Eukaryotes have developed an elaborate series of interactions with bacteria that enter their bodies and/or cells. Genome evolution of symbiotic and parasitic bacteria multiplying inside eukaryotic cells results in both convergent and divergent changes. The genome sequences of the symbiotic bacteria of aphids, Buchnera aphidicola, and the parasitic bacteria of body louse and humans, Rickettsia prowazekii, provide insights into these processes. Convergent genome characteristics include reduction in genome sizes and lowered G+C content values. Divergent evolution was recorded for amino acid and cell wall biosynthetic genes. The presence of pseudogenes in both genomes provides examples of recent gene inactivation events and offers clues to the process of genome deterioration and host-cell adaptation.