Post-Pleistocene radiation of the pea aphid complex revealed by rapidly evolving endosymbionts

Adaptation to different resources has the potential to cause rapid species diversification, but few studies have been able to quantify the time scale of recent adaptive radiations. The pea aphid, Acyrthosiphon pisum, a model of speciation for host-specialized parasites, consists of several biotypes (races or species) living on distinct legume hosts. To document this radiation, we used rapidly evolving sequences from Buchnera, the maternally transmitted bacterial endosymbiont of aphids. Analyses of Buchnera pseudogene sequences revealed that 11 host-associated biotypes sort mostly into distinct matrilines despite low sequence divergence. A calibration based on divergence times of 7 sequenced genomes of Buchnera allowed us to date the last maternal ancestor of these biotypes between 8,000 and 16,000 years, with a burst of diversification at an estimated 3,600-9,500 years. The recency of this diversification, which is supported by microsatellite data, implies that the pea aphid complex ranks among the most rapid adaptive radiations yet documented. This diversification coincides with post-Pleistocene warming and with the domestication and anthropogenic range expansion of several of the legume hosts of pea aphids. Thus, we hypothesize that the new availability or abundance of resources triggered a cascade of divergence events in this newly formed complex.

Two species of symbiotic bacteria present in the soybean aphid (Hemiptera: Aphididae)

Aphids, which feed solely on plant phloem sap, have developed symbiotic associations with bacteria that provide them with the amino acids that are lacking in phloem. Three soybean aphid (Aphis glycines Mat samura) populations were screened for the presence of Buchnera aphidicola and three common species of secondary aphid symbionts (Serratia symbiotica, Hamiltonella defensa, and Regiella insecticola). Diagnostic polymerase chain reaction and subsequent DNA sequencing showed the presence of two species of symbiotic bacteria present in all three soybean aphid populations tested: B. aphidicola and Arsenophonus sp. Although Buchnera is commonly found in aphids, Arsenophonus is most commonly found in whiteflies (Hemiptera: Aleyrodidae), making the soybean aphid unique among aphids that have been tested for the presence of Arsenophonus.

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.

Evolution of the secondary symbiont "Candidatus serratia symbiotica" in aphid species of the subfamily lachninae

Buchnera aphidicola BCc, the primary endosymbiont of the aphid Cinara cedri (subfamily Lachninae), is losing its symbiotic capacity and might be replaced by the coresident "Candidatus Serratia symbiotica." Phylogenetic and morphological analyses within the subfamily Lachninae indicate two different "Ca. Serratia symbiotica" lineages and support the longtime coevolution of both symbionts in C. cedri.

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.

Plasmids derived from Gifsy-1/Gifsy-2, lambdoid prophages contributing to the virulence of Salmonella enterica serovar Typhimurium: implications for the evolution of replication initiation proteins of lambdoid phages and enterobacteria

Gifsy-1 and Gifsy-2 are lambdoid prophages which contribute to the virulence of Salmonella enterica serovar Typhimurium. The nucleotide sequence of the replication region of both prophages is identical, and similar in organization to the replication region of bacteriophage lambda. To investigate the replication of the Gifsy phages and the relationship between Gifsy and host chromosome replication, a plasmid which contained all the genes and regulatory sequences required for autonomous replication in bacterial cells was constructed. This plasmid, pGifsy, was stably maintained in Escherichia coli cells. The helicase loader of the Gifsy phages is very similar to the DnaC protein of the host, a feature characteristic of a large group of prophages common in the sequenced genomes of pathogenic enterobacteria. This DnaC-like protein showed no similarity to the helicase loader of bacteriophage lambda and closely related phages. Interestingly, unlike plasmids derived from bacteriophage lambda (lambda plasmids), pGifsy did not require a gene encoding the putative helicase loader for replication, although deletion of this gene resulted in a decrease in plasmid copy number. Under these conditions, it was shown that the plasmid utilized the helicase loader coded by the host. On the other hand, the viral protein could not substitute for DnaC in bacterial chromosome replication. The results of the current study support the hypothesis that the enterobacterial helicase loader is of viral origin. This hypothesis explains why the gene for DnaC, the protein central to both replication initiation and replication restart in E. coli, is present in the genomes of Escherichia, Shigella, Salmonella and Buchnera, but not in the genomes of related enterobacteria.

Aphid thermal tolerance is governed by a point mutation in bacterial symbionts

Symbiosis is a ubiquitous phenomenon generating biological complexity, affecting adaptation, and expanding ecological capabilities. However, symbionts, which can be subject to genetic limitations such as clonality and genomic degradation, also impose constraints on hosts. A model of obligate symbiosis is that between aphids and the bacterium Buchnera aphidicola, which supplies essential nutrients. We report a mutation in Buchnera of the aphid Acyrthosiphon pisum that recurs in laboratory lines and occurs in field populations. This single nucleotide deletion affects a homopolymeric run within the heat-shock transcriptional promoter for ibpA, encoding a small heat-shock protein. This Buchnera mutation virtually eliminates the transcriptional response of ibpA to heat stress and lowers its expression even at cool or moderate temperatures. Furthermore, this symbiont mutation dramatically affects host fitness in a manner dependent on thermal environment. Following a short heat exposure as juveniles, aphids bearing short-allele symbionts produced few or no progeny and contained almost no Buchnera, in contrast to aphids bearing symbionts without the deletion. Conversely, under constant cool conditions, aphids containing symbionts with the short allele reproduced earlier and maintained higher reproductive rates. The short allele has appreciable frequencies in field populations (up to 20%), further supporting the view that lowering of ibpA expression improves host fitness under some conditions. This recurring Buchnera mutation governs thermal tolerance of aphid hosts. Other cases in which symbiont microevolution has a major effect on host ecological tolerance are likely to be widespread because of the high mutation rates of symbiotic bacteria and their crucial roles in host metabolism and development.

Aphids are sap-feeding insects that depend on obligate bacterial symbionts of the genus Buchnera for biosynthesis of needed nutrients. Studying Buchnera gene expression in pea aphids, we identified a recurring mutation (a single-base deletion) in the transcriptional promoter of the small heat-shock protein, ibpA. This mutation arose and was fixed twice in sublines derived from a single female aphid in the lab and kept at constant 20 °C. Experiments using aphid lines that differed only in the presence of this Buchnera mutation revealed that it eliminates the ibpA transcriptional response to heat shock and affects ibpA expression at low temperatures. In aphids containing Buchnera with the mutation, a short heat treatment as juveniles leads to elimination of most or all symbionts and to reproductive failure; the same treatment has little effect on aphids containing Buchnera without the mutation. Conversely, at constant lower temperatures, aphids with Buchnera bearing this mutation enjoy a reproductive advantage. Pea aphid populations are polymorphic for the Buchnera mutation, suggesting that it is maintained at substantial frequencies by selection. This study indicates that mutations in obligate symbionts can have major consequences for host fitness and geographic distributions.

A recurring single-nucleotide deletion in a heat-shock transcriptional promoter in the bacterialBuchnera symbiont of the pea aphid(Acyrthosiphon pisum) dramatically affects host fitness in a thermal-dependent manner.

Optimization and clinical validation of a pathogen detection microarray

New design and optimization of pathogen detection microarrays is shown to allow robust and accurate detection of a range of pathogens. The customized microarray platform includes a method for reducing PCR bias during DNA amplification.

DNA microarrays used as 'genomic sensors' have great potential in clinical diagnostics. Biases inherent in random PCR-amplification, cross-hybridization effects, and inadequate microarray analysis, however, limit detection sensitivity and specificity. Here, we have studied the relationships between viral amplification efficiency, hybridization signal, and target-probe annealing specificity using a customized microarray platform. Novel features of this platform include the development of a robust algorithm that accurately predicts PCR bias during DNA amplification and can be used to improve PCR primer design, as well as a powerful statistical concept for inferring pathogen identity from probe recognition signatures. Compared to real-time PCR, the microarray platform identified pathogens with 94% accuracy (76% sensitivity and 100% specificity) in a panel of 36 patient specimens. Our findings show that microarrays can be used for the robust and accurate diagnosis of pathogens, and further substantiate the use of microarray technology in clinical diagnostics.

Temporal fragmentation of speciation in bacteria

Because bacterial recombination involves the occasional transfer of small DNA fragments between strains, different sets of niche-specific genes may be maintained in populations that freely recombine at other loci. Therefore, genetic isolation may be established at different times for different chromosomal regions during speciation as recombination at niche-specific genes is curtailed. To test this model, we separated sequence divergence into rate and time components, revealing that different regions of the Escherichia coli and Salmonella enterica chromosomes diverged over a approximately 70-million-year period. Genetic isolation first occurred at regions carrying species-specific genes, indicating that physiological distinctiveness between the nascent Escherichia and Salmonella lineages was maintained for tens of millions of years before the complete genetic isolation of their chromosomes.

The role of mutational dynamics in genome shrinkage

Genome shrinkage occurs after whole genome duplications (WGDs) and in the evolution of parasitic or symbiotic species. The dynamics of this process, whether it occurs by single gene deletions or also by larger deletions are however unknown. In yeast, genome shrinkage has occurred after a WGD. Using a computational model of genome evolution, we show that in a random genome single gene deletions cannot explain the observed pattern of gene loss in yeast. The distribution of genes deleted per event can be very well described by a geometric distribution, with a mean of 1.1 genes per event. In terms of deletions of a stretch of base pairs, we find that a geometric distribution with an average of 500-600 base pairs per event describes the data very well. Moreover, in the model, as in the data, gene pairs that have a small intergenic distance are more likely to be both deleted. This proves that simultaneous deletion of multiple genes causes the observed pattern of gene deletions, rather than deletion of functionally clustered genes by selection. Furthermore, we found that in the bacterium Buchnera aphidicola larger deletions than in yeast are necessary to explain the clustering of deleted genes. We show that the excess clustering of deleted genes in B. aphidicola can be explained by the clustering of genes in operons. Therefore, we show that selection has little effect on the clustering of deleted genes after the WGD in yeast, while it has during genome shrinkage in B. aphidicola.

Genetic evidence from mitochondrial, nuclear, and endosymbiont markers for the evolution of host plant associated species in the aphid genus Hyalopterus (Hemiptera: Aphididae)

Over the past several decades biologists' fascination with plant-herbivore interactions has generated intensive research into the implications of these interactions for insect diversification. The study of closely related phytophagous insect species or populations from an evolutionary perspective can help illuminate ecological and selective forces that drive these interactions. Here we present such an analysis for aphids in the genus Hyalopterus (Hemiptera: Aphididae), a cosmopolitan group that feeds on plants in the genus Prunus (Rosaceae). Hyalopterus currently contains two recognized species associated with different Prunus species, although the taxonomy and evolutionary history of the group is poorly understood. Using mitochondrial COI sequences, 16S rDNA sequences from the aphid endosymbiont Buchnera aphidicola, and nine microsatellite loci we investigated population structure in Hyalopterus from the most commonly used Prunus host species throughout the Mediterranean as well as in California, where the species H. pruni is an invasive pest. We found three deeply divergent lineages structured in large part by specific associations with plum, almond, and peach trees. There was no evidence that geographic or temporal barriers could explain the overall diversity in the genus. Levels of genetic differentiation are consistent with that typically attributed to aphid species and indicate divergence times older than the domestication of Prunus for agriculture. Interestingly, in addition to their typical hosts, aphids from each of the three lineages were frequently found on apricot trees. Apricot also appears to act as a resource mediated hybrid zone for plum and almond associated lineages. Together, results suggest that host plants have played a role in maintaining host-associated differentiation in Hyalopterus for as long as several million years, despite worldwide movement of host plants and the potential for ongoing hybridization.

Conservation of the links between gene transcription and chromosomal organization in the highly reduced genome of Buchnera aphidicola

BACKGROUND

Genomic studies on bacteria have clearly shown the existence of chromosomal organization as regards, for example, to gene localization, order and orientation. Moreover, transcriptomic analyses have demonstrated that, in free-living bacteria, gene transcription levels and chromosomal organization are mutually influenced. We have explored the possible conservation of relationships between mRNA abundances and chromosomal organization in the highly reduced genome of Buchnera aphidicola, the primary endosymbiont of the aphids, and a close relative to Escherichia coli.

RESULTS

Using an oligonucleotide-based microarray, we normalized the transcriptomic data by genomic DNA signals in order to have access to inter-gene comparison data. Our analysis showed that mRNA abundances, gene organization (operon) and gene essentiality are correlated in Buchnera (i.e., the most expressed genes are essential genes organized in operons) whereas no link between mRNA abundances and gene strand bias was found. The effect of Buchnera genome evolution on gene expression levels has also been analysed in order to assess the constraints imposed by the obligate symbiosis with aphids, underlining the importance of some gene sets for the survival of the two partners. Finally, our results show the existence of spatial periodic transcriptional patterns in the genome of Buchnera.

CONCLUSION

Despite an important reduction in its genome size and an apparent decay of its capacity for regulating transcription, this work reveals a significant correlation between mRNA abundances and chromosomal organization of the aphid-symbiont Buchnera.

Strict host-symbiont cospeciation and reductive genome evolution in insect gut bacteria

Host-symbiont cospeciation and reductive genome evolution have been identified in obligate endocellular insect symbionts, but no such example has been identified from extracellular ones. Here we first report such a case in stinkbugs of the family Plataspidae, wherein a specific gut bacterium is vertically transmitted via “symbiont capsule.” In all of the plataspid species, females produced symbiont capsules upon oviposition and their gut exhibited specialized traits for capsule production. Phylogenetic analysis showed that the plataspid symbionts constituted a distinct group in the γ-Proteobacteria, whose sister group was the aphid obligate endocellular symbionts Buchnera. Removal of the symbionts resulted in retarded growth, mortality, and sterility of the insects. The host phylogeny perfectly agreed with the symbiont phylogeny, indicating strict host-symbiont cospeciation despite the extracellular association. The symbionts exhibited AT-biased nucleotide composition, accelerated molecular evolution, and reduced genome size, as has been observed in obligate endocellular insect symbionts. These findings suggest that not the endocellular conditions themselves but the population genetic attributes of the vertically transmitted symbionts are probably responsible for the peculiar genetic traits of these insect symbionts. We proposed the designation “Candidatus Ishikawaella capsulata” for the plataspid symbionts. The plataspid stinkbugs, wherein the host-symbiont associations can be easily manipulated, provide a novel system that enables experimental approaches to previously untouched aspects of the insect-microbe mutualism. Furthermore, comparative analyses of the sister groups, the endocellular Buchnera and the extracellular Ishikawaella, would lead to insights into how the different symbiotic lifestyles have affected their genomic evolution.

Evidence of host-symbiont cospeciation in an insect gut symbiont suggests that long-term vertical transmission and population structure are central forces driving the genomic changes characteristic of insect nutritional symbionts.

Modular organization in the reductive evolution of protein-protein interaction networks

BACKGROUND

The variation in the sizes of the genomes of distinct life forms remains somewhat puzzling. The organization of proteins into domains and the different mechanisms that regulate gene expression are two factors that potentially increase the capacity of genomes to create more complex systems. High-throughput protein interaction data now make it possible to examine the additional complexity generated by the way that protein interactions are organized.

RESULTS

We have studied the reduction in genome size of Buchnera compared to its close relative Escherichia coli. In this well defined evolutionary scenario, we found that among all the properties of the protein interaction networks, it is the organization of networks into modules that seems to be directly related to the evolutionary process of genome reduction.

CONCLUSION

In Buchnera, the apparently non-random reduction of the modular structure of the networks and the retention of essential characteristics of the interaction network indicate that the roles of proteins within the interaction network are important in the reductive process.

Dynamics of reductive genome evolution in mitochondria and obligate intracellular microbes

Reductive evolution in mitochondria and obligate intracellular microbes has led to a significant reduction in their genome size and guanine plus cytosine content (GC). We show that genome shrinkage during reductive evolution in prokaryotes follows an exponential decay pattern and provide a method to predict the extent of this decay on an evolutionary timescale. We validated predictions by comparison with estimated extents of genome reduction known to have occurred in mitochondria and Buchnera aphidicola, through comparative genomics and by drawing on available fossil evidences. The model shows how the mitochondrial ancestor would have quickly shed most of its genome, shortly after its incorporation into the protoeukaryotic cell and prior to codivergence subsequent to the split of eukaryotic lineages. It also predicts that the primary rickettsial parasitic event would have occurred between 180 and 425 million years ago (MYA), an event of relatively recent evolutionary origin considering the fact that Rickettsia and mitochondria evolved from a common alphaproteobacterial ancestor. This suggests that the symbiotic events of Rickettsia and mitochondria originated at different time points. Moreover, our model results predict that the ancestor of Wigglesworthia glossinidia brevipalpis, dated around the time of origin of its symbiotic association with the tsetse fly (50-100 MYA), was likely to have been an endosymbiont itself, thus supporting an earlier proposition that Wigglesworthia, which is currently a maternally inherited primary endosymbiont, evolved from a secondary endosymbiont.

Genome reduction of the aphid endosymbiont Buchnera aphidicola in a recent evolutionary time scale

Genome reduction, a typical feature of symbiotic bacteria, was analyzed in the last stages of evolution of Buchnera aphidicola, the primary aphid endosymbiont, in two neutrally evolving regions: the pseudogene cmk and an intergenic region. These two regions were examined in endosymbionts from several lineages of their aphid host Rhopalosiphum padi, and different species of the same genus, whose divergence times ranged from 0.62 to 19.51 million years. Estimates of nucleotide substitution rates were between 4.3 and 6.7x10(-9) substitution/site/year, with G or C nucleotides being substituted around four times more frequently than A or T. Two different types of indel events were detected, of which many were small (1-10 nt) but one was large (about 200 nucleotides). With respect to the large one and considering the proportion and size of the deletions and insertions, the reduction rate was 1.3x10(-8) lost nucleotides/site/year. We propose a stepwise scenario for the last stages of evolution in B. aphidicola: together with a very slow and gradual degradation, considerable indels would punctually emerge. The only restriction to large deletion fixation is that the lost fragment does not contain essential genes.

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.

Hundreds of flagellar basal bodies cover the cell surface of the endosymbiotic bacterium Buchnera aphidicola sp. strain APS

Buchnera aphidicola is the endosymbiotic bacterium of the pea aphid. Due to its small genome size, Buchnera lacks many essential genes for autogenous life but obtains nutrients from the host. Although the Buchnera cell is nonmotile, it retains clusters of flagellar genes that lack the late genes necessary for motility, including the flagellin gene. In this study, we show that the flagellar genes are actually transcribed and translated and that the Buchnera cell surface is covered with hundreds of hook-basal-body (HBB) complexes. The abundance of HBB complexes suggests a role other than motility. We discuss the possibility that the HBB complex may serve as a protein transporter not only for the flagellar proteins but also for other proteins to maintain the symbiotic system.

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.

Selection acts on DNA secondary structures to decrease transcriptional mutagenesis

Single-stranded DNA is more subject to mutation than double stranded. During transcription, DNA is transiently single stranded and therefore subject to higher mutagenesis. However, if local intra-strand secondary structures are formed, some bases will be paired and therefore less sensitive to mutation than unpaired bases. Using complete genome sequences of Escherichia coli, we show that local intra-strand secondary structures can, as a consequence, be used to define an index of transcription-driven mutability. At gene level, we show that natural selection has favoured a reduced transcription-driven mutagenesis via the higher than expected frequency of occurrence of intra-strand secondary structures. Such selection is stronger in highly expressed genes and suggests a sequence-dependent way to control mutation rates and a novel form of selection affecting the evolution of synonymous mutations.

Genome sequence evolution results from the interplay between mutagenesis and natural selection. Mutations occur as the result of biochemical or physical alteration of DNA and/or from the errors made by polymerases while replicating DNA. As many mutations tend to be detrimental to the organism's fitness, natural selection favours a decrease in mutation rate. Hence, many mechanisms have evolved to control mutation rate. The mechanisms described to date have relied on (i) the existence of enzymes repairing the damaged DNA or correcting mismatched bases, which are mechanisms having an effect on whole genome mutation rate, and (ii) the avoidance in the sequence of repetition that could be misread by the polymerases, which is a sequence-dependent local control of mutation rate. In the present paper, the authors suggest that another sequence-dependent control of mutation exists and shapes the overall evolution of the genome. Using a comparative analysis of Escherichia coli genomes, they show that local secondary structures that are formed during the transcription of genes into RNA can modulate the base-to-base mutation rate. Moreover, the authors show that natural selection seems to have favoured the occurrence of such structures to minimise mutability, especially in the most expressed genes. This paper proposes a new way in which gene sequences can be constrained by natural selection.

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.

Chance and necessity in the evolution of minimal metabolic networks

It is possible to infer aspects of an organism's lifestyle from its gene content. Can the reverse also be done? Here we consider this issue by modelling evolution of the reduced genomes of endosymbiotic bacteria. The diversity of gene content in these bacteria may reflect both variation in selective forces and contingency-dependent loss of alternative pathways. Using an in silico representation of the metabolic network of Escherichia coli, we examine the role of contingency by repeatedly simulating the successive loss of genes while controlling for the environment. The minimal networks that result are variable in both gene content and number. Partially different metabolisms can thus evolve owing to contingency alone. The simulation outcomes do preserve a core metabolism, however, which is over-represented in strict intracellular bacteria. Moreover, differences between minimal networks based on lifestyle are predictable: by simulating their respective environmental conditions, we can model evolution of the gene content in Buchnera aphidicola and Wigglesworthia glossinidia with over 80% accuracy. We conclude that, at least for the particular cases considered here, gene content of an organism can be predicted with knowledge of its distant ancestors and its current lifestyle.