Sources of variation in dietary requirements in an obligate nutritional symbiosis

The nutritional symbiosis between aphids and their obligate symbiont, Buchnera aphidicola, is often characterized as a highly functional partnership in which the symbiont provides the host with essential nutrients. Despite this, some aphid lineages exhibit dietary requirements for nutrients typically synthesized by Buchnera, suggesting that some aspect of the symbiosis is disrupted. To examine this phenomenon in the pea aphid, Acyrthosiphon pisum, populations were assayed using defined artificial diet to determine dietary requirements for essential amino acids (EAAs). Six clones exhibiting dependence on EAAs in their diet were investigated further. In one aphid clone, a mutation in a Buchnera amino acid biosynthesis gene could account for the clone's requirement for dietary arginine. Analysis of aphid F1 hybrids allowed separation of effects of the host and symbiont genomes, and revealed that both affect the requirement for dietary EAAs in the clones tested. Amino acid requirements were minimally affected by secondary symbiont infection. Our results indicate that variation among pea aphids in dependence on dietary amino acids can result from Buchnera mutation as well as variation in the host genotype.

Dynamics of a recurrent Buchnera mutation that affects thermal tolerance of pea aphid hosts

Mutations in maternally transmitted symbionts can affect host fitness. In this study we investigate a mutation in an obligate bacterial symbiont (Buchnera), which has dramatic effects on the heat tolerance of pea aphid hosts (Acyrthosiphon pisum). The heat-sensitive allele arises through a single base deletion in a homopolymer within the promoter of ibpA, which encodes a universal small heat-shock protein. In laboratory cultures reared under cool conditions (20°), the rate of fixation (1.4 × 10(-3) substitutions per Buchnera replication) is much higher than the previously estimated mutation rate for single base deletions in homopolymers in the Buchnera genome, implying a strong selective benefit. This mutation recurs in natural populations, but seldom reaches high frequencies, implying that it is only rarely favored by selection. Another potential source of physiological stress in pea aphids is infection by other microorganisms, including facultative bacterial symbionts, which occur in a majority of pea aphids in field populations. Frequency of the heat-sensitive Buchnera allele is negatively correlated with presence of facultative symbionts in both laboratory colonies and field populations, suggesting that these infections impose stress that is ameliorated by ibpA expression. This single base polymorphism in Buchnera has the potential to allow aphid populations to adapt quickly to prevailing conditions.

A stochastic model for a single click of Muller's ratchet

This work presents a new approach to Muller's ratchet, where Haigh's model is approximately mapped into a simpler model that describes the behaviour of a population after a click of the ratchet, i.e., after loss of what was the fittest class. This new model predicts the distribution of times to the next click of the ratchet and is equivalent to a Wright-Fisher model for a population of haploid asexual individuals with one locus and two alleles. Within this model, the fittest members of a population correspond to carriers of one allele, while all other individuals have suboptimal fitness and are represented as carriers of the other allele. In this way, all suboptimal fitness individuals are amalgamated into a single "mutant" class. The approach presented here has some limitations and the potential for improvement. However, it does lead to results for the rate of the ratchet that, over a wide range of parameters, are accurate within one order of magnitude of simulation results. This contrasts with existing approaches, which are designed for only one or other of the two different parameter regimes known for the ratchet and are more accurate only in the parameter regime they were designed for. Numerical results are presented for the mean time between clicks of the ratchet for (i) the Wright-Fisher model, (ii) a diffusion approximation of this model and (iii) individually based simulations of a full model. The diffusion approximation is validated over a wide range of parameters by its close agreement with the Wright-Fisher model. The present work predicts that: (a) the time between clicks of the ratchet is insensitive to the value of the selection coefficient when the genomic mutation rate is large compared with the selection coefficient against a deleterious mutation, (b) the time interval between clicks of the ratchet has, approximately, an exponential distribution (or its discrete analogue). It is thus possible to determine the variance in times between clicks, given the expected time between clicks. Evidence for both (a) and (b) is seen in simulations.

Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum

Buchnera aphidicola is the primary obligate intracellular symbiont of most aphid species. B. aphidicola and aphids have been evolving in parallel since their association started, about 150 Myr ago. Both partners have lost their autonomy, and aphid diversification has been confined to smaller ecological niches by this co-evolution. B. aphidicola has undergone major genomic and biochemical changes as a result of adapting to intracellular life. Several genomes of B. aphidicola from different aphid species have been sequenced in the last decade, making it possible to carry out analyses and comparative studies using system-level in silico methods. This review attempts to provide a systemic description of the symbiotic function of aphid endosymbionts, particularly of B. aphidicola from the pea aphid Acyrthosiphon pisum, by analyzing their structural genomic properties, as well as their genetic and metabolic networks.

The functional transfer of genes from the mitochondria to the nucleus: the effects of selection, mutation, population size and rate of self-fertilization

The transfer of mitochondrial genes to the nucleus is a recurrent and consistent feature of eukaryotic genome evolution. Although many theories have been proposed to explain such transfers, little relevant data exist. The observation that clonal and self-fertilizing plants transfer more mitochondrial genes to their nuclei than do outcrossing plants contradicts predictions of major theories based on nuclear recombination and leaves a gap in our conceptual understanding how the observed pattern of gene transfer could arise. Here, with a series of deterministic and stochastic simulations, we show how epistatic selection and relative mutation rates of mitochondrial and nuclear genes influence mitochondrial-to-nuclear gene transfer. Specifically, we show that when there is a benefit to having a mitochondrial gene present in the nucleus, but absent in the mitochondria, self-fertilization dramatically increases both the rate and the probability of gene transfer. However, absent such a benefit, when mitochondrial mutation rates exceed those of the nucleus, self-fertilization decreases the rate and probability of transfer. This latter effect, however, is much weaker than the former. Our results are relevant to understanding the probabilities of fixation when loci in different genomes interact.

Genome-wide functional divergence after the symbiosis of proteobacteria with insects unraveled through a novel computational approach

Symbiosis has been among the most important evolutionary steps to generate biological complexity. The establishment of symbiosis required an intimate metabolic link between biological systems with different complexity levels. The strict endo-cellular symbiotic bacteria of insects are beautiful examples of the metabolic coupling between organisms belonging to different kingdoms, a eukaryote and a prokaryote. The host (eukaryote) provides the endosymbiont (prokaryote) with a stable cellular environment while the endosymbiont supplements the host's diet with essential metabolites. For such communication to take place, endosymbionts' genomes have suffered dramatic modifications and reconfigurations of proteins' functions. Two of the main modifications, loss of genes redundant for endosymbiotic bacteria or the host and bacterial genome streamlining, have been extensively studied. However, no studies have accounted for possible functional shifts in the endosymbiotic proteomes. Here, we develop a simple method to screen genomes for evidence of functional divergence between two species clusters, and we apply it to identify functional shifts in the endosymbiotic proteomes. Despite the strong effects of genetic drift in the endosymbiotic systems, we unexpectedly identified genes to be under stronger selective constraints in endosymbionts of aphids and ants than in their free-living bacterial relatives. These genes are directly involved in supplementing the host's diet with essential metabolites. A test of functional divergence supports a strong relationship between the endosymbiosis and the functional shifts of proteins involved in the metabolic communication with the insect host. The correlation between functional divergence in the endosymbiotic bacterium and the ecological requirements of the host uncovers their intimate biochemical and metabolic communication and provides insights on the role of symbiosis in generating species diversity.

Biological complexity has emerged on earth by the combination of living forms. This combination, called symbiosis, had to overcome the problems caused by the uncoupled metabolisms of the organisms involved. One way to do so was through the loss of genes that were no longer needed for the endosymbiont in the protected cellular environment provided by the host. Another step necessary to adjust both metabolisms was through the change in the function of bacterial proteins to perform new roles in the symbiotic system. In this article, we test such events in symbiotic systems involving an insect and a bacterium by developing a new and simple method to identify proteome-wide functional shifts. Our results show that most of the functional changes occurred at genes involved in metabolic communication with the host and are correlated with the host's ecological traits.

Genomics and evolution of heritable bacterial symbionts

Insect heritable symbionts have proven to be ubiquitous, based on molecular screening of various insect lineages. Recently, molecular and experimental approaches have yielded an immensely richer understanding of their diverse biological roles, resulting in a burgeoning research literature. Increasingly, commonalities and intermediates are being discovered between categories of symbionts once considered distinct: obligate mutualists that provision nutrients, facultative mutualists that provide protection against enemies or stress, and symbionts such as Wolbachia that manipulate reproductive systems. Among the most far-reaching impacts of widespread heritable symbiosis is that it may promote speciation by increasing reproductive and ecological isolation of host populations, and it effectively provides a means for transfer of genetic information among host lineages. In addition, insect symbionts provide some of the extremes of cellular genomes, including the smallest and the fastest evolving, raising new questions about the limits of evolution of life.

Revealing the hidden complexities of mtDNA inheritance

Mitochondrial DNA (mtDNA) is a pivotal tool in molecular ecology, evolutionary and population genetics. The power of mtDNA analyses derives from a relatively high mutation rate and the apparent simplicity of mitochondrial inheritance (maternal, without recombination), which has simplified modelling population history compared to the analysis of nuclear DNA. However, in biology things are seldom simple, and advances in DNA sequencing and polymorphism detection technology have documented a growing list of exceptions to the central tenets of mitochondrial inheritance, with paternal leakage, heteroplasmy and recombination now all documented in multiple systems. The presence of paternal leakage, recombination and heteroplasmy can have substantial impact on analyses based on mtDNA, affecting phylogenetic and population genetic analyses, estimates of the coalescent and the myriad of other parameters that are dependent on such estimates. Here, we review our understanding of mtDNA inheritance, discuss how recent findings mean that established ideas may need to be re-evaluated, and we assess the implications of these new-found complications for molecular ecologists who have relied for decades on the assumption of a simpler mode of inheritance. We show how it is possible to account for recombination and heteroplasmy in evolutionary and population analyses, but that accurate estimates of the frequencies of biparental inheritance and recombination are needed. We also suggest how nonclonal inheritance of mtDNA could be exploited, to increase the ways in which mtDNA can be used in analyses.

The causes of mutation accumulation in mitochondrial genomes

A fundamental observation across eukaryotic taxa is that mitochondrial genomes have a higher load of deleterious mutations than nuclear genomes. Identifying the evolutionary forces that drive this difference is important to understanding the rates and patterns of sequence evolution, the efficacy of natural selection, the maintenance of sex and recombination and the mechanisms underlying human ageing and many diseases. Recent studies have implicated the presumed asexuality of mitochondrial genomes as responsible for their high mutational load. We review the current body of knowledge on mitochondrial mutation accumulation and recombination, and conclude that asexuality, per se, may not be the primary determinant of the high mutation load in mitochondrial DNA (mtDNA). Very little recombination is required to counter mutation accumulation, and recent evidence suggests that mitochondrial genomes do experience occasional recombination. Instead, a high rate of accumulation of mildly deleterious mutations in mtDNA may result from the small effective population size associated with effectively haploid inheritance. This type of transmission is nearly ubiquitous among mitochondrial genomes. We also describe an experimental framework using variation in mating system between closely related species to disentangle the root causes of mutation accumulation in mitochondrial genomes.

Selection for translational robustness in Buchnera aphidicola, endosymbiotic bacteria of aphids

Its strong intergenerational bottlenecks and effectively asexual reproduction have led Buchnera aphidicola, the endocellular symbiotic bacterium of aphids, to spectacular evolutionary and genomic changes in comparison with its free-living bacterial cousins. These changes summarize into high fixation rates of mildly deleterious destabilizing mutations. This predicts a sharp decline of its fitness and the consequent early demise of this endosymbiotic bacterium. Its survival for hundreds of millions of years casts doubt on genetic drift as the sole evolutionary force and seeks further explanation. We identify in Buchnera selection to increase the robustness of proteins to misfolding translation errors. Translational robustness varies between Buchnera lineages and protein functional categories. Metabolic proteins have been under selection for translational efficiency, whereas evolutionary rates of proteins involved in fundamental cellular processes have been largely determined by selection for translational robustness. We detect the strongest signal of translational robustness in B. aphidicola Cinara cedri with a very similar pattern to that inferred for the most common symbiotic ancestor of Buchnera lineages. This indicates that B. aphidicola Cinara cedri lineage may have probably reached the minimum evolutionary stable gene composition for endosymbiotic lifestyle. The evolutionary patterns from the comparative genomic analyses of these endosymbionts support a paradoxically complex dynamic for apparently simple genomes.

Salinity and temperature effects on physiological responses of Vibrio fischeri from diverse ecological niches

Vibrio fischeri is a bioluminescent bacterial symbiont of sepiolid squids (Cephalopoda: Sepiolidae) and monocentrid fishes (Actinopterygii: Monocentridae). V. fischeri exhibit competitive dominance within the allopatrically distributed squid genus Euprymna, which have led to the evolution of V. fischeri host specialists. In contrast, the host genus Sepiola contains sympatric species that is thought to have given rise to V. fischeri that have evolved as host generalists. Given that these ecological lifestyles may have a direct effect upon the growth spectrum and survival limits in contrasting environments, optimal growth ranges were obtained for numerous V. fischeri isolates from both free-living and host environments. Upper and lower limits of growth were observed in sodium chloride concentrations ranging from 0.0% to 9.0%. Sepiola symbiotic isolates possessed the least variation in growth throughout the entire salinity gradient, whereas isolates from Euprymna were the least uniform at <2.0% NaCl. V. fischeri fish symbionts (CG101 and MJ101) and all free-living strains were the most dissimilar at >5.0% NaCl. Growth kinetics of symbiotic V. fischeri strains were also measured under a range of salinity and temperature combinations. Symbiotic V. fischeri ES114 and ET101 exhibited a synergistic effect for salinity and temperature, where significant differences in growth rates due to salinity existed only at low temperatures. Thus, abiotic factors such as temperature and salinity have differential effects between free-living and symbiotic strains of V. fischeri, which may alter colonization efficiency prior to infection.

Accelerated evolutionary rates in tropical and oceanic parmelioid lichens (Ascomycota)

Background

The rate of nucleotide substitutions is not constant across the Tree of Life, and departures from a molecular clock have been commonly reported. Within parmelioid lichens, the largest group of macrolichens, large discrepancies in branch lengths between clades were found in previous studies. Using an extended taxon sampling, we test for presence of significant rate discrepancies within and between these clades and test our a priori hypothesis that such rate discrepancies may be explained by shifts in moisture regime or other environmental conditions.

Results

In this paper, the first statistical evidence for accelerated evolutionary rate in lichenized ascomycetes is presented. Our results give clear evidence for a faster rate of evolution in two Hypotrachyna clades that includes species occurring in tropical and oceanic habitats in comparison with clades consisting of species occurring in semi-arid and temperate habitats. Further we explore potential links between evolutionary rates and shifts in habitat by comparing alternative Ornstein-Uhlenbeck models.

Conclusion

Although there was only weak support for a shift at the base of a second tropical clade, where the observed nucleotide substitution rate is high, overall support for a shift in environmental conditions at cladogenesis is very strong. This suggests that speciation in some lichen clades has proceeded by dispersal into a novel environment, followed by radiation within that environment. We found moderate support for a shift in moisture regime at the base of one tropical clade and a clade occurring in semi-arid regions and a shift in minimum temperature at the base of a boreal-temperate clade.

An effect of 16S rRNA intercistronic variability on coevolutionary analysis in symbiotic bacteria: molecular phylogeny of Arsenophonus triatominarum

The genes of ribosomal RNA are the most popular and frequently used markers for bacterial phylogeny and reconstruction of insect-symbiont coevolution. In primary symbionts, such as Buchnera and Wigglesworthia, genome economization leads to the establishment of a single copy of these sequences. In phylogenetic studies, they provide sufficient information and yield phylogenetic trees congruent with host evolution. In contrast, other symbiotic lineages (e.g., the genus Arsenophonus) carry a higher number of rRNA copies in their genomes, which may have serious consequences for phylogenetic inference. In this study, we show that in Arsenophonus triatominarum the degree of heterogeneity can affect reconstruction of phylogenetic relationships and mask possible coevolution between the symbiont and its host. Phylogenetic arrangement of individual rRNA copies was used, together with a calculation of their divergence time, to demonstrate that the incongruent 16S rDNA trees and low nucleotide diversity in the secondary symbiont could be reconciled with the coevolutionary scenario.

Population structure, levels of selection, and the evolution of intracellular symbionts

Many obligately intracellular symbionts exhibit a characteristic set of genetic changes that include an increase in substitution rates, loss of many genes, and apparent destabilization of many proteins and structural RNAs. Authors have suggested that these changes are due to increased mutation rates, or, more commonly, decreased effective population size due to population bottlenecks at the symbiont or, perhaps, host level. I propose that the increase in substitution rates and accumulation of deleterious mutations is a consequence of the population structure imposed on the endosymbionts by strict host association, loss of horizontal transmission and potentially conflicting levels of selection. I analyze a population genetic model of endosymbiont evolution, and demonstrate that substitution rates will increase, and the effect of those substitutions on endosymbiont fitness will become more deleterious as horizontal transmission among hosts decreases. Additionally, I find that there is a critical level of horizontal transmission below which natural selection cannot effectively purge deleterious mutations, leading to an expected loss of fitness over time. This critical level varies across loci with the degree of correlation between host and endosymbiont fitness, and may help explain differential retention and loss of certain genes.

Mutational interference and the progression of Muller's ratchet when mutations have a broad range of deleterious effects

Deleterious mutations can accumulate in asexual haploid genomes through the process known as Muller's ratchet. This process has been described in the literature mostly for the case where all mutations are assumed to have the same effect on fitness. In the more realistic situation, deleterious mutations will affect fitness with a wide range of effects, from almost neutral to lethal. To elucidate the behavior of the ratchet in this more realistic case, simulations were carried out in a number of models, one where all mutations have the same effect on selection [one-dimensional (1D) model], one where the deleterious mutations can be divided into two groups with different selective effects [two-dimensional (2D) model], and finally one where the deleterious effects are distributed. The behavior of these models suggests that deleterious mutations can be classified into three different categories, such that the behavior of each can be described in a straightforward way. This makes it possible to predict the ratchet rate for an arbitrary distribution of fitness effects using the results for the well-studied 1D model with a single selection coefficient. The description was tested and shown to work well in simulations where selection coefficients are derived from an exponential distribution.

Muller’s ratchet in symbiont populations

Muller's ratchet, the inevitable accumulation of deleterious mutations in asexual populations, has been proposed as a major factor in genome degradation of obligate symbiont organisms. Essentially, if left unchecked the ratchet will with certainty cause extinction due to the ever increasing mutational load. This paper examines the evolutionary fate of insect symbionts, using mathematical modelling to simulate the accumulation of deleterious mutations. We investigate the effects of a hierarchical two level population structure. Since each host contains its own subpopulation of symbionts, there will be a large number of small symbiont populations linked indirectly via selection on the host level. We show that although the separate subpopulations will accumulate deleterious mutations quickly, the symbiont population as a whole will be protected from extinction by selection acting on the hosts. As a consequence, the extent of genome degradation observed in present day symbionts is more likely to represent loss of functions that were (near-) neutral to the host, rather than a snap shot of a decline towards complete genetic collapse.

GRAST: a new way of genome reduction analysis using comparative genomics

MOTIVATION

Establishment of intra-cellular life involved a profound re-configuration of the genetic characteristics of bacteria, including genome reduction and rearrangements. Understanding the mechanisms underlying these phenomena will shed light on the genome rearrangements essential for the development of an intra-cellular lifestyle. Comparison of genomes with differences in their sizes poses statistical as well as computational problems. Little efforts have been made to develop flexible computational tools with which to analyse genome reduction and rearrangements.

RESULTS

Investigation of genome reduction and rearrangements in endosymbionts using a novel computational tool (GRAST) identified gathering of genes with similar functions. Conserved clusters of functionally related genes (CGSCs) were detected. Heterogeneous gene and gene cluster non-functionalization/loss are identified between genome regions, functional gene categories and during evolution. Results show that gene non-functionalisation has accelerated during the last 50 MY of Buchnera's evolution while CGSCs have been static.

Biology bacteriocyte-associated endosymbionts of plant sap-sucking insects

Psyllids, whiteflies, aphids, and mealybugs are members of the suborder Sternorrhyncha and share a common property, namely the utilization of plant sap as their food source. Each of these insect groups has an obligatory association with a different prokaryotic endosymbiont, and the association is the result of a single infection followed by maternal, vertical transmission of the endosymbionts. The result of this association is the domestication of the free-living bacterium to serve the purposes of the host, namely the synthesis of essential amino acids. This domestication is probably in all cases accompanied by a major reduction in genome size. The different properties of the genomes and fragments of the genomes of these endosymbionts suggest that there are different constraints on the permissible evolutionary changes that are probably a function of the gene repertoire of the endosymbiont ancestor and the gene losses that occurred during the reduction of genome size.

Inheritance and recombination of mitochondrial genomes in plants, fungi and animals

It is generally assumed that mitochondrial genomes are uniparentally transmitted, homoplasmic and nonrecombining. However, these assumptions draw largely from early studies on animal mitochondrial DNA (mtDNA). In this review, we show that plants, animals and fungi are all characterized by episodes of biparental inheritance, recombination among genetically distinct partners, and selfish elements within the mitochondrial genome, but that the extent of these phenomena may vary substantially across taxa. We argue that occasional biparental mitochondrial transmission may allow organisms to achieve the best of both worlds by facilitating mutational clearance but continuing to restrict the spread of selfish genetic elements. We also show that methodological biases and disproportionately allocated study effort are likely to have influenced current estimates of the extent of biparental inheritance, heteroplasmy and recombination in mitochondrial genomes from different taxa. Despite these complications, there do seem to be discernible similarities and differences in transmission dynamics and likelihood of recombination of mtDNA in plant, animal and fungal taxa that should provide an excellent opportunity for comparative investigation of the evolution of mitochondrial genome dynamics.

The roles of positive and negative selection in the molecular evolution of insect endosymbionts

The evolutionary rate acceleration observed in most endosymbiotic bacteria may be explained by higher mutation rates, changes in selective pressure, and increased fixation of deleterious mutations by genetic drift. Here, we explore the forces influencing molecular evolution in Blochmannia, an obligate endosymbiont of Camponotus and related ant genera. Our goals were to compare rates of sequence evolution in Blochmannia with related bacteria, to explore variation in the strength and efficacy of negative (purifying) selection, and to evaluate the effect of positive selection. For six Blochmannia pairs, plus Buchnera and related enterobacteria, estimates of sequence divergence at four genes confirm faster rates of synonymous evolution in the ant mutualist. This conclusion is based on higher dS between Blochmannia lineages despite their more recent divergence. Likewise, generally higher dN in Blochmannia indicates faster rates of nonsynonymous substitution in this group. One exception is the groEL gene, for which lower dN and dN/dS compared to Buchnera indicate exceptionally strong negative selection in Blochmannia. In addition, we explored evidence for positive selection in Blochmannia using both site-and lineage-based maximum likelihood models. These approaches confirmed heterogeneity of dN/dS among codon sites and revealed significant variation in dN/dS across Blochmannia lineages for three genes. Lineage variation affected genes independently, with no evidence of parallel changes in dN/dS across genes along a given branch. Our data also reveal instances of dN/dS greater than one; however, we do not interpret these large dN/dS ratios as evidence for positive selection. In sum, while drift may contribute to an overall rate acceleration at nonsynonymous sites in Blochmannia, variable selective pressures best explain the apparent gene-specific changes in dN/dS across lineages of this ant mutualist. In the course of this study, we reanalyzed variation at Buchnera groEL and found no evidence of positive selection that was previously reported.

Protein Interactions Limit the Rate of Evolution of Photosynthetic Genes in Cyanobacteria

Using a bioinformatic approach, we analyzed the correspondence in genetic distance matrices between all possible pairwise combinations of 82 photosynthetic genes in 10 species of cyanobacteria. Our analysis reveals significant correlations between proteins linked in a conserved gene order and between structurally identified interacting protein scaffolds that coordinate the binding of cofactors involved in photosynthetic electron transport. Analyses of amino acid substitution rates suggest that the tempo of evolution of genes encoding core metabolic processes in the photosynthetic apparatus is highly constrained by protein-protein, protein-lipid, and protein-cofactor interactions (collectively called "protein interactions"). These interactions are critical for energy transduction, primary charge separation, and electron transport and effectively act as an internal selection pressure governing the conservation of clusters of photosynthetic genes in oxygenic prokaryotic photoautotrophs. Consequently, although several proteins within the photosynthetic apparatus are biophysically and physiologically inefficient, selection has not significantly altered the genes encoding these essential proteins over billions of years of evolution. In effect, these core proteins have become "frozen metabolic accidents."

Regulation of transcription in a reduced bacterial genome: nutrient-provisioning genes of the obligate symbiont Buchnera aphidicola

Buchnera aphidicola, the obligate symbiont of aphids, has an extremely reduced genome, of which about 10% is devoted to the biosynthesis of essential amino acids needed by its hosts. Most regulatory genes for these pathways are absent, raising the question of whether and how transcription of these genes responds to the major shifts in dietary amino acid content encountered by aphids. Using full-genome microarrays for B. aphidicola of the host Schizaphis graminum, we examined transcriptome responses to changes in dietary amino acid content and then verified behavior of individual transcripts using quantitative reverse transcriptase PCR. The only gene showing a consistent and substantial (>twofold) response was metE, which underlies methionine biosynthesis and which is the only amino acid biosynthetic gene retaining its ancestral regulator (metR). In another aphid host, Acyrthosiphon pisum, B. aphidicola has no functional metR and shows no response in metE transcript levels to changes in amino acid concentrations. Thus, the only substantial transcriptional response involves the one gene for which an ancestral regulator is retained. This result parallels that from a previous study on heat stress, in which only the few genes retaining the global heat shock promoter showed responses in transcript abundance. The irreversible losses of transcriptional regulators constrain ability to alter gene expression in the context of environmental fluctuations affecting the symbiotic partners.

Germline bottlenecks, biparental inheritance and selection on mitochondrial variants: a two-level selection model

Selection on mitochondrial mutations potentially occurs at different levels: at the mitochondria, cell, and organism levels. Several factors affect the strength of selection at these different levels; in particular, mitochondrial bottlenecks during germline development and reduced paternal transmission decrease the genetic variance within cells, while they increase the variance between cells and between organisms, thus decreasing the strength of selection within cells and increasing the strength of selection between cells and organisms. However, bottlenecks and paternal transmission also affect the effective mitochondrial population size, thus affecting genetic drift. In this article, we use a simple model of a unicellular life cycle to investigate the effects of bottlenecks and paternal transmission on the probability of fixation of mitochondrial mutants and their frequency at mutation-selection equilibrium. We find that bottlenecks and reduced paternal transmission decrease the mean frequency of alleles with sm>sc (approximately), where sm and sc are the strengths of selection for an allele within and between cells, respectively, and increase the frequency of alleles with sm<sc. Effects on fixation probabilities are different; for example, bottlenecks reduce the fixation probability of mutants with sm>0 (unless sm is very small relative to sc) and increase the fixation probability of mutants with sm<0.

Endosymbiosis: lessons in conflict resolution

Endosymbiotic bacteria live within a host species. There are many and diverse examples of such relationships, the study of which provides important lessons for ecology and evolution

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.

Transmission of symbiotic bacteria Buchnera to parthenogenetic embryos in the aphid Acyrthosiphon pisum (Hemiptera: Aphidoidea)

All phloem-feeding aphids have an absolute requirement for their primary bacterial symbionts Buchnera sp. The bacteria are transmitted vertically to either embryos in the viviparous morph or to eggs in the oviparous morph, with the implication that the symbiont population regularly passes through a genetic 'bottleneck', i.e. only a small proportion of the maternal symbiont population is transmitted to offspring. In this paper, we visualise this process in viviparous aphids using a specific immunolabelling technique for Buchnera. The images show a stream of bacteria originating from a single mycetocyte and entering the embryo, possibly via a membranous conduit, and individual bacterial cells free in the haemocoel of the aphid. Staining within the embryo blastoderm suggests over expression of antigen, perhaps indicative of rapid bacterial division immediately following infection.

Genome evolution in an insect cell: distinct features of an ant-bacterial partnership

Bacteria that live exclusively within eukaryotic host cells include not only well-known pathogens, but also obligate mutualists, many of which occur in diverse insect groups such as aphids, psyllids, tsetse flies, and the ant genus Camponotus (Buchner, 1965; Douglas, 1998; Moran and Telang, 1998; Baumann et al., 2000; Moran and Baumann, 2000). In contrast to intracellular pathogens, these primary (P) endosymbionts of insects are required for the survival and reproduction of the host, exist within specialized host cells called bacteriocytes, and undergo stable maternal transmission through host lineages (Buchner, 1965; McLean and Houk, 1973). Due to their long-term host associations and close phylogenetic relationship with well-characterized enterobacteria (Fig. 1), P-endosymbionts of insects are ideal model systems to examine changes in genome content and architecture that occur in the context of beneficial, intracellular associations. Since these bacteria have not been cultured outside of the host cell, they are difficult to study with traditional genetic or physiological approaches. However, in recent years, molecular and computational approaches have provided important insights into their genetic diversity and ecological significance. This review describes some recent insights into the evolutionary genetics of obligate insect-bacteria symbioses, with a particular focus on an intriguing association between the bacterial endosymbiont Blochmannia and its ant hosts.

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.

Molecular evidence for host-adapted races of the fungal endophyte Epichloë bromicola after presumed host shifts

Host shifts of plant-feeding insects and parasites promote adaptational changes that may result in the formation of host races, an assumed intermediate stage in sympatric speciation. Here, we report on genetically differentiated and host-adapted races of the fungal endophyte Epichloë bromicola, which presumably emerged after a shift from the grass Bromus erectus to other Bromus hosts. Fungi of the genus Epichloë (Ascomycota) and related anamorphs of Neotyphodium are widespread endophytes of cool-season grasses. Sexually reproducing strains sterilize the host by formation of external fruiting structures (stromata), whereas asexual strains are asymptomatic and transmitted via seeds. In E. bromicola, strains infecting B. erectus are sexual, and strains from two woodland species, B. benekenii and B. ramosus, are asexual and seed transmitted. Analyses of amplified fragment length polymorphism fingerprinting and of intron sequences of the tub2 and tef1 genes of 26 isolates from the three Bromus hosts collected at natural sites in Switzerland and nearby France demonstrated that isolates are genetically differentiated according to their host, indicating that E. bromicola does not form a single, randomly mating population. Phylogenetic analyses of sequence data did not unambiguously resolve the exact origin of asexual E. bromicola strains, but it is likely they arose from within sexual populations on B. erectus. Incongruence of trees derived from different genes may have resulted from recombination at some time in the recent history of host strains. Reciprocal inoculations of host plant seedlings showed that asexual isolates from B. benekenii and B. ramosus were incapable of infecting B. erectus, whereas the sexual isolates from B. erectus retained the assumed ancestral trait of broad compatibility with Bromus host seedlings. Because all isolates were interfertile in experimental crosses, asexual strains may not be considered independent biological species. We suggest that isolates infecting B. benekenii and B. ramosus represent long-standing host races or incipient species that emerged after host shifts and that may evolve through host-mediated reproductive isolation toward independent species.

Genome evolution in bacterial endosymbionts of insects

Many insect species rely on intracellular bacterial symbionts for their viability and fecundity. Large-scale DNA-sequence analyses are revealing the forces that shape the evolution of these bacterial associates and the genetic basis of their specialization to an intracellular lifestyle. The full genome sequences of two obligate mutualists, Buchnera aphidicola of aphids and Wigglesworthia glossinidia of tsetse flies, reveal substantial gene loss and an integration of host and symbiont metabolic functions. Further genomic comparisons should reveal the generality of these features among bacterial mutualists and the extent to which they are shared with other intracellular bacteria, including obligate pathogens.

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.

Recent and ancient asexuality in Timema walkingsticks

Determining the evolutionary age of asexual lineages should help in inferring the temporal scale under which asexuality and sex evolve and assessing selective factors involved in the evolution of asexuality. We used 416 bp of the mitochondrial COI gene to infer phylogenetic relationships of virtually all known Timema walkingstick species, including extensive intraspecific sampling for all five of the asexuals and their close sexual relatives. The asexuals T. douglasi and T. shepardii were very closely related to each other and evolutionarily young (less than 0.5 million years old). For the asexuals T. monikensis and T. tahoe, evidence for antiquity was weak since only one population of each was sampled, intraspecific divergences were low, and genetic distances to related sexuals were high: maximum-likelihood molecular-clock age estimates ranged from 0.26 to 2.39 million years in T. monikensis and from 0.29-1.06 million years in T. tahoe. By contrast, T. genevieve was inferred to be an ancient asexual, with an age of 0.81 to 1.42 million years. The main correlate of the age of asexual lineages was their geographic position, with younger asexuals being found further north.

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.

Intraspecific variation in symbiont genomes: bottlenecks and the aphid-buchnera association

Buchnera are maternally transmitted bacterial endosymbionts that synthesize amino acids that are limiting in the diet of their aphid hosts. Previous studies demonstrated accelerated sequence evolution in Buchnera compared to free-living bacteria, especially for nonsynonymous substitutions. Two mechanisms may explain this acceleration: relaxed purifying selection and increased fixation of slightly deleterious alleles under drift. Here, we test the divergent predictions of these hypotheses for intraspecific polymorphism using Buchnera associated with natural populations of the ragweed aphid, Uroleucon ambrosiae. Contrary to expectations under relaxed selection, U. ambrosiae from across the United States yielded strikingly low sequence diversity at three Buchnera loci (dnaN, trpBC, trpEG), revealing polymorphism three orders of magnitude lower than in enteric bacteria. An excess of nonsynonymous polymorphism and of rare alleles was also observed. Local sampling of additional dnaN sequences revealed similar patterns of polymorphism and no evidence of food plant-associated genetic structure. Aphid mitochondrial sequences further suggested that host bottlenecks and large-scale dispersal may contribute to genetic homogenization of aphids and symbionts. Together, our results support reduced N(e) as a primary cause of accelerated sequence evolution in Buchnera. However, our study cannot rule out the possibility that mechanisms other than bottlenecks also contribute to reduced N(e) at aphid and endosymbiont loci.

Vertical transmission of biosynthetic plasmids in aphid endosymbionts (Buchnera)

This study tested for horizontal transfer of plasmids among Buchnera aphidicola strains associated with ecologically and phylogenetically related aphid hosts (Uroleucon species). Phylogenetic congruence of Buchnera plasmid (trpEG and leuABC) and chromosomal (dnaN and trpB) genes supports strictly vertical long-term transmission of plasmids, which persist due to their contributions to host nutrition rather than capacity for infectious transfer. Synonymous divergences indicate elevated mutation on plasmids relative to chromosomal genes.