Why are the genomes of endosymbiotic bacteria so stable?

A matter of background: DNA repair pathways as a possible cause for the sparse distribution of CRISPR-Cas systems in bacteria

The absence of CRISPR-Cas systems in more than half of the sequenced bacterial genomes is intriguing, because their role in adaptive immunity and their frequent transfer between species should have made them almost ubiquitous, as is the case in Archaea. Here, we investigate the possibility that the success of CRISPR-Cas acquisition by horizontal gene transfer is affected by the interactions of these systems with the host genetic background and especially with components of double-strand break repair systems (DSB-RS). We first described the distribution of systems specialized in the repair of double-strand breaks in Bacteria: homologous recombination and non-homologous end joining. This allowed us to show that such systems are more often positively or negatively correlated with the frequency of CRISPR-Cas systems than random genes of similar frequency. The detailed analysis of these co-occurrence patterns shows that our method identifies previously known cases of mechanistic interactions between these systems. It also reveals other positive and negative patterns of co-occurrence between DSB-RS and CRISPR-Cas systems. Notably, it shows that the patterns of distribution of CRISPR-Cas systems in Proteobacteria are strongly dependent on the epistatic groups including RecBCD and AddAB. Our results suggest that the genetic background plays an important role in the success of adaptive immunity in different bacterial clades and provide insights to guide further experimental research on the interactions between CRISPR-Cas and DSB-RS. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Random Genetic Drift and Selective Pressures Shaping the Blattabacterium Genome

Estimates suggest that at least half of all extant insect genera harbor obligate bacterial mutualists. Whereas an endosymbiotic relationship imparts many benefits upon host and symbiont alike, the intracellular lifestyle has profound effects on the bacterial genome. The obligate endosymbiont genome is a product of opposing forces: genes important to host survival are maintained through physiological constraint, contrasted by the fixation of deleterious mutations and genome erosion through random genetic drift. The obligate cockroach endosymbiont, Blattabacterium - providing nutritional augmentation to its host in the form of amino acid synthesis - displays radical genome alterations when compared to its most recent free-living relative Flavobacterium. To date, eight Blattabacterium genomes have been published, affording an unparalleled opportunity to examine the direction and magnitude of selective forces acting upon this group of symbionts. Here, we find that the Blattabacterium genome is experiencing a 10-fold increase in selection rate compared to Flavobacteria. Additionally, the proportion of selection events is largely negative in direction, with only a handful of loci exhibiting signatures of positive selection. These findings suggest that the Blattabacterium genome will continue to erode, potentially resulting in an endosymbiont with an even further reduced genome, as seen in other insect groups such as Hemiptera.

Inhibition of NHEJ repair by type II-A CRISPR-Cas systems in bacteria

Type II CRISPR-Cas systems introduce double-strand breaks into DNA of invading genetic material and use DNA fragments to acquire novel spacers during adaptation. These breaks can be the substrate of several DNA repair pathways, paving the way for interactions. We report that non-homologous end-joining (NHEJ) and type II-A CRISPR-Cas systems only co-occur once among 5563 fully sequenced prokaryotic genomes. We investigated experimentally the possible molecular interactions using the NHEJ pathway from Bacillus subtilis and the type II-A CRISPR-Cas systems from Streptococcus thermophilus and Streptococcus pyogenes. Our results suggest that the NHEJ system has no effect on CRISPR immunity. On the other hand, we provide evidence for the inhibition of NHEJ repair by the Csn2 protein. Our findings give insights on the complex interactions between CRISPR-Cas systems and repair mechanisms in bacteria, contributing to explain the scattered distribution of CRISPR-Cas systems in bacterial genome.

The Evolutionary Relationship between Alternative Splicing and Gene Duplication

The protein diversity that exists today has resulted from various evolutionary processes. It is well known that gene duplication (GD) along with the accumulation of mutations are responsible, among other factors, for an increase in the number of different proteins. The gene structure in eukaryotes requires the removal of non-coding sequences, introns, to produce mature mRNAs. This process, known as cis-splicing, referred to here as splicing, is regulated by several factors which can lead to numerous splicing arrangements, commonly designated as alternative splicing (AS). AS, producing several transcripts isoforms form a single gene, also increases the protein diversity. However, the evolution and manner for increasing protein variation differs between AS and GD. An important question is how are patterns of AS affected after a GD event. Here, we review the current knowledge of AS and GD, focusing on their evolutionary relationship. These two processes are now considered the main contributors to the increasing protein diversity and therefore their relationship is a relevant, yet understudied, area of evolutionary study.

Identification and analysis of integrons and cassette arrays in bacterial genomes

Integrons recombine gene arrays and favor the spread of antibiotic resistance. Their broader roles in bacterial adaptation remain mysterious, partly due to lack of computational tools. We made a program – IntegronFinder – to identify integrons with high accuracy and sensitivity. IntegronFinder is available as a standalone program and as a web application. It searches for attC sites using covariance models, for integron-integrases using HMM profiles, and for other features (promoters, attI site) using pattern matching. We searched for integrons, integron-integrases lacking attC sites, and clusters of attC sites lacking a neighboring integron-integrase in bacterial genomes. All these elements are especially frequent in genomes of intermediate size. They are missing in some key phyla, such as α-Proteobacteria, which might reflect selection against cell lineages that acquire integrons. The similarity between attC sites is proportional to the number of cassettes in the integron, and is particularly low in clusters of attC sites lacking integron-integrases. The latter are unexpectedly abundant in genomes lacking integron-integrases or their remains, and have a large novel pool of cassettes lacking homologs in the databases. They might represent an evolutionary step between the acquisition of genes within integrons and their stabilization in the new genome.

Found and Lost: The Fates of Horizontally Acquired Genes in Arthropod-Symbiotic Spiroplasma

Horizontal gene transfer (HGT) is an important mechanism that contributed to biological diversity, particularly in bacteria. Through acquisition of novel genes, the recipient cell may change its ecological preference and the process could promote speciation. In this study, we determined the complete genome sequence of two Spiroplasma species for comparative analyses and inferred the putative gene gains and losses. Although most Spiroplasma species are symbionts of terrestrial insects, Spiroplasma eriocheiris has evolved to be a lethal pathogen of freshwater crustaceans. We found that approximately 7% of the genes in this genome may have originated from HGT and these genes expanded the metabolic capacity of this organism. Through comparison with the closely related Spiroplasma atrichopogonis, as well as other more divergent lineages, our results indicated that these HGT events could be traced back to the most recent common ancestor of these two species. However, most of these horizontally acquired genes have been pseudogenized in S. atrichopogonis, suggesting that they did not contribute to the fitness of this lineage that maintained the association with terrestrial insects. Thus, accumulation of small deletions that disrupted these foreign genes was not countered by natural selection. On the other hand, the long-term survival of these horizontally acquired genes in the S. eriocheiris genome hinted that they might play a role in the ecological shift of this species. Finally, the implications of these findings and the conflicts among gene content, 16S rRNA gene sequencing, and serological typing, are discussed in light of defining bacterial species.

Two host clades, two bacterial arsenals: evolution through gene losses in facultative endosymbionts

Bacterial endosymbiosis is an important evolutionary process in insects, which can harbor both obligate and facultative symbionts. The evolution of these symbionts is driven by evolutionary convergence, and they exhibit among the tiniest genomes in prokaryotes. The large host spectrum of facultative symbionts and the high diversity of strategies they use to infect new hosts probably impact the evolution of their genome and explain why they undergo less severe genomic erosion than obligate symbionts. Candidatus Hamiltonella defensa is suitable for the investigation of the genomic evolution of facultative symbionts because the bacteria are engaged in specific relationships in two clades of insects. In aphids, H. defensa is found in several species with an intermediate prevalence and confers protection against parasitoids. In whiteflies, H. defensa is almost fixed in some species of Bemisia tabaci, which suggests an important role of and a transition toward obligate symbiosis. In this study, comparisons of the genome of H. defensa present in two B. tabaci species (Middle East Asia Minor 1 and Mediterranean) and in the aphid Acyrthosiphon pisum revealed that they belong to two distinct clades and underwent specific gene losses. In aphids, it contains highly virulent factors that could allow protection and horizontal transfers. In whiteflies, the genome lost these factors and seems to have a limited ability to acquire genes. However it contains genes that could be involved in the production of essential nutrients, which is consistent with a primordial role for this symbiont. In conclusion, although both lineages of H. defensa have mutualistic interactions with their hosts, their genomes follow distinct evolutionary trajectories that reflect their phenotype and could have important consequences on their evolvability.

Signs of neutralization in a redundant gene involved in homologous recombination in Wolbachia endosymbionts

Genomic reduction in bacterial endosymbionts occurs through large genomic deletions and long-term accumulation of mutations. The latter process involves successive steps including gene neutralization, pseudogenization, and gradual erosion until complete loss. Although many examples of pseudogenes at various levels of degradation have been reported, neutralization cases are scarce because of the transient nature of the process. Gene neutralization may occur due to relaxation of selection in nonessential genes, for example, those involved in redundant functions. Here, we report an example of gene neutralization in the homologous recombination (HR) pathway of Wolbachia, a bacterial endosymbiont of arthropods and nematodes. The HR pathway is often depleted in endosymbiont genomes, but it is apparently intact in some Wolbachia strains. Analysis of 12 major HR genes showed that they have been globally under strong purifying selection during the evolution of Wolbachia strains hosted by arthropods, supporting the evolutionary importance of the HR pathway for these Wolbachia genomes. However, we detected signs of recent neutralization of the ruvA gene in a subset of Wolbachia strains, which might be related to an ancestral, clade-specific amino acid change that impaired DNA-binding activity. Strikingly, RuvA is part of the RuvAB complex involved in branch migration, whose function overlaps with the RecG helicase. Although ruvA is experiencing neutralization, recG is under strong purifying selection. Thus, our high phylogenetic resolution suggests that we identified a rare example of targeted neutralization of a gene involved in a redundant function in an endosymbiont genome.

Comparative genomics reveals surprising divergence of two closely related strains of uncultivated UCYN-A cyanobacteria

Marine planktonic cyanobacteria capable of fixing molecular nitrogen (termed 'diazotrophs') are key in biogeochemical cycling, and the nitrogen fixed is one of the major external sources of nitrogen to the open ocean. Candidatus Atelocyanobacterium thalassa (UCYN-A) is a diazotrophic cyanobacterium known for its widespread geographic distribution in tropical and subtropical oligotrophic oceans, unusually reduced genome and symbiosis with a single-celled prymnesiophyte alga. Recently a novel strain of this organism was also detected in coastal waters sampled from the Scripps Institute of Oceanography pier. We analyzed the metagenome of this UCYN-A2 population by concentrating cells by flow cytometry. Phylogenomic analysis provided strong bootstrap support for the monophyly of UCYN-A (here called UCYN-A1) and UCYN-A2 within the marine Crocosphaera sp. and Cyanothece sp. clade. UCYN-A2 shares 1159 of the 1200 UCYN-A1 protein-coding genes (96.6%) with high synteny, yet the average amino-acid sequence identity between these orthologs is only 86%. UCYN-A2 lacks the same major pathways and proteins that are absent in UCYN-A1, suggesting that both strains can be grouped at the same functional and ecological level. Our results suggest that UCYN-A1 and UCYN-A2 had a common ancestor and diverged after genome reduction. These two variants may reflect adaptation of the host to different niches, which could be coastal and open ocean habitats.

Genome sequence of Candidatus Riesia pediculischaeffi, endosymbiont of chimpanzee lice, and genomic comparison of recently acquired endosymbionts from human and chimpanzee lice

The obligate-heritable endosymbionts of insects possess some of the smallest known bacterial genomes. This is likely due to loss of genomic material during symbiosis. The mode and rate of this erosion may change over evolutionary time: faster in newly formed associations and slower in long-established ones. The endosymbionts of human and anthropoid primate lice present a unique opportunity to study genome erosion in newly established (or young) symbionts. This is because we have a detailed phylogenetic history of these endosymbionts with divergence dates for closely related species. This allows for genome evolution to be studied in detail and rates of change to be estimated in a phylogenetic framework. Here, we sequenced the genome of the chimpanzee louse endosymbiont (Candidatus Riesia pediculischaeffi) and compared it with the closely related genome of the human body louse endosymbiont. From this comparison, we found evidence for recent genome erosion leading to gene loss in these endosymbionts. Although gene loss was detected, it was not significantly greater than in older endosymbionts from aphids and ants. Additionally, we searched for genes associated with B-vitamin synthesis in the two louse endosymbiont genomes because these endosymbionts are believed to synthesize essential B vitamins absent in the louse’s diet. All of the expected genes were present, except those involved in thiamin synthesis. We failed to find genes encoding for proteins involved in the biosynthesis of thiamin or any complete exogenous means of salvaging thiamin, suggesting there is an undescribed mechanism for the salvage of thiamin. Finally, genes encoding for the pantothenate de novo biosynthesis pathway were located on a plasmid in both taxa along with a heat shock protein. Movement of these genes onto a plasmid may be functionally and evolutionarily significant, potentially increasing production and guarding against the deleterious effects of mutation. These data add to a growing resource of obligate endosymbiont genomes and to our understanding of the rate and mode of genome erosion in obligate animal-associated bacteria. Ultimately sequencing additional louse p-endosymbiont genomes will provide a model system for studying genome evolution in obligate host associated bacteria.

Small but powerful, the primary endosymbiont of moss bugs, Candidatus Evansia muelleri, holds a reduced genome with large biosynthetic capabilities

Moss bugs (Coleorrhyncha: Peloridiidae) are members of the order Hemiptera, and like many hemipterans, they have symbiotic associations with intracellular bacteria to fulfill nutritional requirements resulting from their unbalanced diet. The primary endosymbiont of the moss bugs, Candidatus Evansia muelleri, is phylogenetically related to Candidatus Carsonella ruddii and Candidatus Portiera aleyrodidarum, primary endosymbionts of psyllids and whiteflies, respectively. In this work, we report the genome of Candidatus Evansia muelleri Xc1 from Xenophyes cascus, which is the only obligate endosymbiont present in the association. This endosymbiont possesses an extremely reduced genome similar to Carsonella and Portiera. It has crossed the borderline to be considered as an autonomous cell, requiring the support of the insect host for some housekeeping cell functions. Interestingly, in spite of its small genome size, Evansia maintains enriched amino acid (complete or partial pathways for ten essential and six nonessential amino acids) and sulfur metabolisms, probably related to the poor diet of the insect, based on bryophytes, which contains very low levels of nitrogenous and sulfur compounds. Several facts, including the congruence of host (moss bugs, whiteflies, and psyllids) and endosymbiont phylogenies and the retention of the same ribosomal RNA operon during genome reduction in Evansia, Portiera, and Carsonella, suggest the existence of an ancient endosymbiotic Halomonadaceae clade associated with Hemiptera. Three possible scenarios for the origin of these three primary endosymbiont genera are proposed and discussed.

Bacterial insertion sequences: their genomic impact and diversity

Insertion sequences (ISs), arguably the smallest and most numerous autonomous transposable elements (TEs), are important players in shaping their host genomes. This review focuses on prokaryotic ISs. We discuss IS distribution and impact on genome evolution. We also examine their effects on gene expression, especially their role in activating neighbouring genes, a phenomenon of particular importance in the recent upsurge of bacterial antibiotic resistance. We explain how ISs are identified and classified into families by a combination of characteristics including their transposases (Tpases), their overall genetic organisation and the accessory genes which some ISs carry. We then describe the organisation of autonomous and nonautonomous IS‐related elements. This is used to illustrate the growing recognition that the boundaries between different types of mobile element are becoming increasingly difficult to define as more are being identified. We review the known Tpase types, their different catalytic activities used in cleaving and rejoining DNA strands during transposition, their organisation into functional domains and the role of this in regulation. Finally, we consider examples of prokaryotic IS domestication. In a more speculative section, we discuss the necessity of constructing more quantitative dynamic models to fully appreciate the continuing impact of TEs on prokaryotic populations.

Genome degeneration and adaptation in a nascent stage of symbiosis

Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis-allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.

The evolution of genomic instability in the obligate endosymbionts of whiteflies

Many insects depend on ancient associations with intracellular bacteria to perform essential metabolic functions. These endosymbionts exhibit striking examples of convergence in genome architecture, including a high degree of structural stability that is not typical of their free-living counterparts. However, the recently sequenced genome of the obligate whitefly endosymbiont Portiera revealed features that distinguish it from other ancient insect associates, such as a low gene density and the presence of perfectly duplicated sequences. Here, we report the comparative analysis of Portiera genome sequences both within and between host species. In one whitefly lineage (Bemisia tabaci), we identify large-scale structural polymorphisms in the Portiera genome that exist even within individual insects. This variation is likely mediated by recombination across identical repeats that are maintained by gene conversion. The complete Portiera genome sequence from a distantly related whitefly host (Trialeurodes vaporarium) confirms a history of extensive genome rearrangement in this ancient endosymbiont. Using gene-order-based phylogenetic analysis, we show that the majority of rearrangements have occurred in the B. tabaci lineage, coinciding with an increase in the rate of nucleotide substitutions, a proliferation of short tandem repeats (microsatellites) in intergenic regions, and the loss of many widely conserved genes involved in DNA replication, recombination, and repair. These results indicate that the loss of recombinational machinery is unlikely to be the cause of the extreme structural conservation that is generally observed in obligate endosymbiont genomes and that large, repetitive intergenic regions are an important substrate for genomic rearrangements.

Genome evolution and phylogenomic analysis of Candidatus Kinetoplastibacterium, the betaproteobacterial endosymbionts of Strigomonas and Angomonas

It has been long known that insect-infecting trypanosomatid flagellates from the genera Angomonas and Strigomonas harbor bacterial endosymbionts (Candidatus Kinetoplastibacterium or TPE [trypanosomatid proteobacterial endosymbiont]) that supplement the host metabolism. Based on previous analyses of other bacterial endosymbiont genomes from other lineages, a stereotypical path of genome evolution in such bacteria over the duration of their association with the eukaryotic host has been characterized. In this work, we sequence and analyze the genomes of five TPEs, perform their metabolic reconstruction, do an extensive phylogenomic analyses with all available Betaproteobacteria, and compare the TPEs with their nearest betaproteobacterial relatives. We also identify a number of housekeeping and central metabolism genes that seem to have undergone positive selection. Our genome structure analyses show total synteny among the five TPEs despite millions of years of divergence, and that this lineage follows the common path of genome evolution observed in other endosymbionts of diverse ancestries. As previously suggested by cell biology and biochemistry experiments, Ca. Kinetoplastibacterium spp. preferentially maintain those genes necessary for the biosynthesis of compounds needed by their hosts. We have also shown that metabolic and informational genes related to the cooperation with the host are overrepresented amongst genes shown to be under positive selection. Finally, our phylogenomic analysis shows that, while being in the Alcaligenaceae family of Betaproteobacteria, the closest relatives of these endosymbionts are not in the genus Bordetella as previously reported, but more likely in the Taylorella genus.

Functional and evolutionary analysis of the genome of an obligate fungal symbiont

Nutritional symbionts of insects include some of the most bizarre genomes studied to date, with extremely reduced size, biased base composition, and limited metabolic abilities. A monophyletic group of aphids within the subfamily Cerataphidinae have lost the bacterial symbiont common to all other Aphididae (Buchnera aphidicola), which have been replaced by a eukaryotic one, the yeast-like symbiont (YLS). As symbionts are expected to experience reduced effective population size (Ne) and largely clonal life cycles, we used this system as a model to test the hypothesis that chronically high levels of genetic drift will result in an increase in size of a eukaryotic genome. We sequenced the genome of the YLS of the aphid Cerataphis brasiliensis and observed elevated rates of protein sequence evolution and intron proliferation in YLS orthologs relative to those of its closest-sequenced relative, consistent with predictions. A moderate amount of repetitive DNA was found along with evidence of directed mutation to prevent proliferation of repetitive elements. Despite increased intron numbers, the overall genome structure appears not to have undergone massive expansion and is around 25 Mb in size. Compared with Buchnera, the YLS appears to have a much broader metabolic repertoire, though many gene families have been reduced in the YLS relative to related fungi. The patterns observed in the YLS genome suggest that its symbiotic lifestyle is permissive to intron proliferation and accelerated sequence evolution, though other factors appear to limit its overall genome expansion.

A cost-effective and universal strategy for complete prokaryotic genomic sequencing proposed by computer simulation

BACKGROUND

Pyrosequencing techniques allow scientists to perform prokaryotic genome sequencing to achieve the draft genomic sequences within a few days. However, the assemblies with shotgun sequencing are usually composed of hundreds of contigs. A further multiplex PCR procedure is needed to fill all the gaps and link contigs into complete chromosomal sequence, which is the basis for prokaryotic comparative genomic studies. In this article, we study various pyrosequencing strategies by simulated assembling from 100 prokaryotic genomes.

FINDINGS

Simulation study shows that a single end 454 Jr. run combined with a paired end 454 Jr. run (8 kb library) can produce: 1) ~90% of 100 assemblies with < 10 scaffolds and ~95% of 100 assemblies with < 150 contigs; 2) average contig N50 size is over 331 kb; 3) average single base accuracy is > 99.99%; 4) average false gene duplication rate is < 0.7%; 5) average false gene loss rate is < 0.4%.

CONCLUSIONS

A single end 454 Jr. run combined with a paired end 454 Jr. run (8 kb library) is a cost-effective way for prokaryotic whole genome sequencing. This strategy provides solution to produce high quality draft assemblies for most of prokaryotic organisms within days. Due to the small number of assembled scaffolds, the following multiplex PCR procedure (for gap filling) would be easy. As a result, large scale prokaryotic whole genome sequencing projects may be finished within weeks.

Metabolic networks of Sodalis glossinidius: a systems biology approach to reductive evolution

BACKGROUND

Genome reduction is a common evolutionary process affecting bacterial lineages that establish symbiotic or pathogenic associations with eukaryotic hosts. Such associations yield highly reduced genomes with greatly streamlined metabolic abilities shaped by the type of ecological association with the host. Sodalis glossinidius, the secondary endosymbiont of tsetse flies, represents one of the few complete genomes available of a bacterium at the initial stages of this process. In the present study, genome reduction is studied from a systems biology perspective through the reconstruction and functional analysis of genome-scale metabolic networks of S. glossinidius.

RESULTS

The functional profile of ancestral and extant metabolic networks sheds light on the evolutionary events underlying transition to a host-dependent lifestyle. Meanwhile, reductive evolution simulations on the extant metabolic network can predict possible future evolution of S. glossinidius in the context of genome reduction. Finally, knockout simulations in different metabolic systems reveal a gradual decrease in network robustness to different mutational events for bacterial endosymbionts at different stages of the symbiotic association.

CONCLUSIONS

Stoichiometric analysis reveals few gene inactivation events whose effects on the functionality of S. glossinidius metabolic systems are drastic enough to account for the ecological transition from a free-living to host-dependent lifestyle. The decrease in network robustness across different metabolic systems may be associated with the progressive integration in the more stable environment provided by the insect host. Finally, reductive evolution simulations reveal the strong influence that external conditions exert on the evolvability of metabolic systems.

An interdependent metabolic patchwork in the nested symbiosis of mealybugs

Highly reduced genomes of 144-416 kilobases have been described from nutrient-provisioning bacterial symbionts of several insect lineages [1-5]. Some host insects have formed stable associations with pairs of bacterial symbionts that live in specialized cells and provide them with essential nutrients; genomic data from these systems have revealed remarkable levels of metabolic complementarity between the symbiont pairs [3, 4, 6, 7]. The mealybug Planococcus citri (Hemiptera: Pseudococcidae) contains dual bacterial symbionts existing with an unprecedented organization: an unnamed gammaproteobacteria, for which we propose the name Candidatus Moranella endobia, lives inside the betaproteobacteria Candidatus Tremblaya princeps [8]. Here we describe the complete genomes and metabolic contributions of these unusual nested symbionts. We show that whereas there is little overlap in retained genes involved in nutrient production between symbionts, several essential amino acid pathways in the mealybug assemblage require a patchwork of interspersed gene products from Tremblaya, Moranella, and possibly P. citri. Furthermore, although Tremblaya has the smallest cellular genome yet described, it contains a genomic inversion present in both orientations in individual insects, starkly contrasting with the extreme structural stability typical of highly reduced bacterial genomes [4, 9, 10].

Function-specific accelerations in rates of sequence evolution suggest predictable epistatic responses to reduced effective population size

Changes in effective population size impinge on patterns of molecular evolution. Notably, slightly deleterious mutations are more likely to drift to fixation in smaller populations, which should typically also lead to an overall acceleration in the rates of evolution. This prediction has been validated empirically for several endosymbiont and island taxa. Here, we first show that rate accelerations are also evident in bacterial pathogens whose recent shifts in virulence make them prime candidates for reduced effective population size: Bacillus anthracis, Bordetella parapertussis, Mycobacterium leprae, Salmonella enterica typhi, Shigella spp., and Yersinia pestis. Using closely related genomes to analyze substitution rate dynamics across six phylogenetically independent bacterial clades, we demonstrate that relative rates of coding sequence evolution are biased according to gene functional category. Notably, genes that buffer against slightly deleterious mutations, such as chaperones, experience stronger rate accelerations than other functional classes at both nonsynonymous and synonymous sites. Although theory predicts altered evolutionary dynamics for buffer loci in the face of accumulating deleterious mutations, to observe even stronger rate accelerations is surprising. We suggest that buffer loci experience elevated substitution rates because the accumulation of deleterious mutations in the remainder of the genome favors compensatory substitutions in trans. Critically, the hyper-acceleration is evident across phylogenetically independent clades, supporting the hypothesis that reductions in effective population size predictably induce epistatic responses in genes that buffer against slightly deleterious mutations.

Protein homeostasis and the phenotypic manifestation of genetic diversity: principles and mechanisms

Changing a single nucleotide in a genome can have profound consequences under some conditions, but the same change can have no consequences under others. Indeed, organisms can be surprisingly robust to environmental and genetic perturbations. Yet, the mechanisms underlying such robustness are controversial. Moreover, how they might affect evolutionary change remains enigmatic. Here, we review the recently appreciated central role of protein homeostasis in buffering and potentiating genetic variation and discuss how these processes mediate the critical influence of the environment on the relationship between genotype and phenotype. Deciphering how robustness emerges from biological organization and the mechanisms by which it is overcome in changing environments will lead to a more complete understanding of both fundamental evolutionary processes and diverse human diseases.

Symbionts and pathogens: what is the difference?

The ecological relationships that organisms establish with others can be considered as broad and diverse as the forms of life that inhabit and interact in our planet. Those interactions can be considered as a continuum spectrum, ranging from beneficial to detrimental outcomes. However, this picture has revealed as more complex and dynamic than previously thought, involving not only factors that affect the two or more members that interact, but also external forces, with chance playing a crucial role in this interplay. Thus, defining a particular symbiont as mutualist or pathogen in an exclusive way, based on simple rules of classification is increasingly challenging if not unfeasible, since new methodologies are providing more evidences that depict exceptions, reversions and transitions within either side of this continuum, especially evident at early stages of symbiotic associations. This imposes a wider and more dynamic view of a complex landscape of interactions.

Structure and dynamics of the operon map of Buchnera aphidicola sp. strain APS

BACKGROUND

Gene expression regulation is still poorly documented in bacteria with highly reduced genomes. Understanding the evolution and mechanisms underlying the regulation of gene transcription in Buchnera aphidicola, the primary endosymbiont of aphids, is expected both to enhance our understanding of this nutritionally based association and to provide an intriguing case-study of the evolution of gene expression regulation in a reduced bacterial genome.

RESULTS

A Bayesian predictor was defined to infer the B. aphidicola transcription units, which were further validated using transcriptomic data and RT-PCR experiments. The characteristics of B. aphidicola predicted transcription units (TUs) were analyzed in order to evaluate the impact of operon map organization on the regulation of gene transcription.On average, B. aphidicola TUs contain more genes than those of E. coli. The global layout of B. aphidicola operon map was mainly shaped by the big reduction and the rearrangements events, which occurred at the early stage of the symbiosis. Our analysis suggests that this operon map may evolve further only by small reorganizations around the frontiers of B. aphidicola TUs, through promoter and/or terminator sequence modifications and/or by pseudogenization events. We also found that the need for specific transcription regulation exerts some pressure on gene conservation, but not on gene assembling in the operon map in Buchnera. Our analysis of the TUs spacing pointed out that a selection pressure is maintained on the length of the intergenic regions between divergent adjacent gene pairs.

CONCLUSIONS

B. aphidicola can seemingly only evolve towards a more polycistronic operon map. This implies that gene transcription regulation is probably subject to weak selection pressure in Buchnera conserving operons composed of genes with unrelated functions.

Functional convergence in reduced genomes of bacterial symbionts spanning 200 My of evolution

The main genomic changes in the evolution of host-restricted microbial symbionts are ongoing inactivation and loss of genes combined with rapid sequence evolution and extreme structural stability; these changes reflect high levels of genetic drift due to small population sizes and strict clonality. This genomic erosion includes irreversible loss of genes in many functional categories and can include genes that underlie the nutritional contributions to hosts that are the basis of the symbiotic association. Candidatus Sulcia muelleri is an ancient symbiont of sap-feeding insects and is typically coresident with another bacterial symbiont that varies among host subclades. Previously sequenced Sulcia genomes retain pathways for the same eight essential amino acids, whereas coresident symbionts synthesize the remaining two. Here, we describe a dual symbiotic system consisting of Sulcia and a novel species of Betaproteobacteria, Candidatus Zinderia insecticola, both living in the spittlebug Clastoptera arizonana. This Sulcia has completely lost the pathway for the biosynthesis of tryptophan and, therefore, retains the ability to make only 7 of the 10 essential amino acids. Zinderia has a tiny genome (208 kb) and the most extreme nucleotide base composition (13.5% G + C) reported to date, yet retains the ability to make the remaining three essential amino acids, perfectly complementing capabilities of the coresident Sulcia. Combined with the results from related symbiotic systems with complete genomes, these data demonstrate the critical role that bacterial symbionts play in the host insect’s biology and reveal one outcome following the loss of a critical metabolic activity through genome reduction.

Mobile genetic element proliferation and gene inactivation impact over the genome structure and metabolic capabilities of Sodalis glossinidius, the secondary endosymbiont of tsetse flies

BACKGROUND

Genome reduction is a common evolutionary process in symbiotic and pathogenic bacteria. This process has been extensively characterized in bacterial endosymbionts of insects, where primary mutualistic bacteria represent the most extreme cases of genome reduction consequence of a massive process of gene inactivation and loss during their evolution from free-living ancestors. Sodalis glossinidius, the secondary endosymbiont of tsetse flies, contains one of the few complete genomes of bacteria at the very beginning of the symbiotic association, allowing to evaluate the relative impact of mobile genetic element proliferation and gene inactivation over the structure and functional capabilities of this bacterial endosymbiont during the transition to a host dependent lifestyle.

RESULTS

A detailed characterization of mobile genetic elements and pseudogenes reveals a massive presence of different types of prophage elements together with five different families of IS elements that have proliferated across the genome of Sodalis glossinidius at different levels. In addition, a detailed survey of intergenic regions allowed the characterization of 1501 pseudogenes, a much higher number than the 972 pseudogenes described in the original annotation. Pseudogene structure reveals a minor impact of mobile genetic element proliferation in the process of gene inactivation, with most of pseudogenes originated by multiple frameshift mutations and premature stop codons. The comparison of metabolic profiles of Sodalis glossinidius and tsetse fly primary endosymbiont Wiglesworthia glossinidia based on their whole gene and pseudogene repertoires revealed a novel case of pathway inactivation, the arginine biosynthesis, in Sodalis glossinidius together with a possible case of metabolic complementation with Wigglesworthia glossinidia for thiamine biosynthesis.

CONCLUSIONS

The complete re-analysis of the genome sequence of Sodalis glossinidius reveals novel insights in the evolutionary transition from a free-living ancestor to a host-dependent lifestyle, with a massive proliferation of mobile genetic elements mainly of phage origin although with minor impact in the process of gene inactivation that is taking place in this bacterial genome. The metabolic analysis of the whole endosymbiotic consortia of tsetse flies have revealed a possible phenomenon of metabolic complementation between primary and secondary endosymbionts that can contribute to explain the co-existence of both bacterial endosymbionts in the context of the tsetse host.

Correlations between bacterial ecology and mobile DNA

Several factors can affect the density of mobile DNA in bacterial genomes including rates of exposure to novel gene pools, recombination, and reductive evolution. These traits are difficult to measure across a broad range of bacterial species, but the ecological niches occupied by an organism provide some indication of the relative magnitude of these forces. Here, by analyzing 384 bacterial genomes assigned to three ecological categories (obligate intracellular, facultative intracellular, and extracellular), we address two, related questions: How does the density of mobile DNA vary across the Bacteria? And is there a statistically supported relationship between ecological niche and mobile element gene density? We report three findings. First, the fraction of mobile element genes in bacterial genomes ranges from 0 to 21% and decreases significantly: facultative intracellular > extracellular > obligate intracellular bacteria. Results further show that the obligate intracellular bacteria that host switch have a higher mobile DNA gene density than the obligate intracellular bacteria that are vertically transmitted. Second, while bacteria from the three ecological niches differ in their average mobile DNA contents, the ranges of mobile DNA found in each category overlap a surprising extent, suggesting bacteria with different lifestyles can tolerate similar amounts of mobile DNA. Third, mobile DNA gene densities increase with genome size across the entire dataset, and the significance of this correlation is dependent on the obligate intracellular bacteria. Further, mobile DNA gene densities do not correlate with evolutionary relationships in a 16S rDNA phylogeny. These findings statistically support a compelling link between mobile element evolution and bacterial ecology.

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.

Convergent evolution of metabolic roles in bacterial co-symbionts of insects

A strictly host-dependent lifestyle has profound evolutionary consequences for bacterial genomes. Most prominent is a sometimes-dramatic amount of gene loss and genome reduction. Recently, highly reduced genomes from the co-resident intracellular symbionts of sharpshooters were shown to exhibit a striking level of metabolic interdependence. One symbiont, called Sulcia muelleri (Bacteroidetes), can produce eight of the 10 essential amino acids, despite having a genome of only 245 kb. The other, Baumannia cicadellinicola (gamma-Proteobacteria), can produce the remaining two essential amino acids as well as many vitamins. Cicadas also contain the symbiont Sulcia, but lack Baumannia and instead contain the co-resident symbiont Hodgkinia cicadicola (alpha-Proteobacteria). Here we report that, despite at least 200 million years of divergence, the two Sulcia genomes have nearly identical gene content and gene order. Additionally, we show that despite being phylogenetically distant and drastically different in genome size and architecture, Hodgkinia and Baumannia have converged on gene sets conferring similar capabilities for essential amino acid biosynthesis, in both cases precisely complementary to the pathways conserved in Sulcia. In contrast, they have completely divergent capabilities for vitamin biosynthesis. Despite having the smallest gene set known in bacteria, Hodgkinia devotes at least 7% of its proteome to cobalamin (vitamin B(12)) biosynthesis, a significant metabolic burden. The presence of these genes can be explained by Hodgkinia's retention of the cobalamin-dependent version of methionine synthase instead of the cobalamin-independent version found in Baumannia, a situation that necessitates retention of cobalamin biosynthetic capabilities to make the essential amino acid methionine.

Phylogenetic ubiquity and shuffling of the bacterial RecBCD and AddAB recombination complexes

RecBCD and AddAB are bacterial enzymes that share similar helicase and nuclease activities and initiate repair of DNA double-strand breaks by homologous recombination. Examination of the phylogenetic distribution of AddAB and RecBCD revealed that one or the other complex is present in most sequenced bacteria. In addition, horizontal gene transfer (HGT) events involving addAB and recBCD appear to be common, with the genes encoding one complex frequently replacing those encoding the other. HGT may also explain the unexpected identification of archaeal addAB genes. More than 85% of addAB and recBCD genes are clustered on the genome, suggesting operon structures. A few organisms, including the Mycobacteria, encode multiple copies of these complexes of either the same or mixed classes. The possibility that the enzymatic activities of the AddAB and RecBCD enzymes promote their horizontal transfer is discussed, and the distribution of AddAB/RecBCD is compared to that of the RecU/RuvC resolvases. Finally, it appears that two sequence motifs, the Walker A box involved in ATP binding and an iron-sulfur-cysteine cluster, are present only in subsets of AddB proteins, suggesting the existence of mechanistically distinct classes of AddB.

Gene conversion maintains nonfunctional transposable elements in an obligate mutualistic endosymbiont

Long-term bacterial endosymbionts typically exhibit reduced genomes, lack genes encoding recombination functions and transposable elements, such as insertion sequences (ISs). In sharp contrast, I found that ISs constitute 2.4% of the genome of the obligate mutualistic endosymbiont Wolbachia wBm. Although no IS copy is transpositionally functional, I show that ISs persist in wBm because of frequent recombinational gene conversion (GC) homogenizing homologous IS sequences. These results not only indicate that there exists a functional recombination molecular machinery in wBm, but they also suggest that, by slowing down the rate of IS degradation and loss, GC may represent a major force influencing reductive evolution in wBm.

Evidence for homologous recombination in intracellular chemosynthetic clam symbionts

Homologous recombination is a fundamental mechanism for the genetic diversification of free-living bacteria. However, recombination may be limited in endosymbiotic bacteria, as these taxa are locked into an intracellular niche and may rarely encounter sources of foreign DNA. This study tested the hypothesis that vertically transmitted endosymbionts of deep-sea clams (Bivalvia: Vesicomyidae) show little or no evidence of recombination. Phylogenetic analysis of 13 loci distributed across the genomes of 14 vesicomyid symbionts revealed multiple, well-supported inconsistencies among gene tree topologies, and maximum likelihood-based tests rejected a hypothesis of shared evolutionary history (linkage) among loci. Further, multiple statistical methods confirmed the presence of recombination by detecting intragenic breakpoints in two symbiont loci. Recombination may be confined to a subset of vesicomyid symbionts, as some clades showed high levels of genomic stability, whereas others showed clear patterns of homologous exchange. Notably, a mosaic genome is present in symB, a symbiont lineage shown to have been acquired laterally (i.e., nonvertically) by Vesicomya sp. JdF clams. The majority of loci analyzed here supported a tight sister clustering of symB with the symbiont of a host species from the Mid-Atlantic Ridge, whereas others placed symB in a clade with symA, the dominant phylotype of V. sp. JdF clams. This result raises the hypothesis that lateral symbiont transfer between hosts may facilitate recombination by bringing divergent symbiont lineages into contact. Together, the data show that homologous recombination contributes to the diversification of vesicomyid clam symbionts, despite the intracellular lifestyle of these bacteria.

Symbioses: a key driver of insect physiological processes, ecological interactions, evolutionary diversification, and impacts on humans

Symbiosis is receiving increased attention among all aspects of biology because of the unifying themes it helps construct across ecological, evolutionary, developmental, semiochemical, and pest management theory. Insects show a vast array of symbiotic relationships with a wide diversity of microorganisms. These relationships may confer a variety of benefits to the host (macrosymbiont), such as direct or indirect nutrition, ability to counter the defenses of plant or animal hosts, protection from natural enemies, improved development and reproduction, and communication. Benefits to the microsymbiont (including a broad range of fungi, bacteria, mites, nematodes, etc.) often include transport, protection from antagonists, and protection from environmental extremes. Symbiotic relationships may be mutualistic, commensal, competitive, or parasitic. In many cases, individual relationships may include both beneficial and detrimental effects to each partner during various phases of their life histories or as environmental conditions change. The outcomes of insect-microbial interactions are often strongly mediated by other symbionts and by features of the external and internal environment. These outcomes can also have important effects on human well being and environmental quality, by affecting agriculture, human health, natural resources, and the impacts of invasive species. We argue that, for many systems, our understanding of symbiotic relationships will advance most rapidly where context dependency and multipartite membership are integrated into existing conceptual frameworks. Furthermore, the contribution of entomological studies to overall symbiosis theory will be greatest where preoccupation with strict definitions and artificial boundaries is minimized, and integration of emerging molecular and quantitative techniques is maximized. We highlight symbiotic relations involving bark beetles to illustrate examples of the above trends.

The evolutionary history of symbiotic associations among bacteria and their animal hosts: a model

A model to explain the evolutionary history of animal-bacteria obligatory mutualistic symbiosis is presented. Dispensability of genes and genetic isolation are key factors in the reduction process of these bacterial genomes. Major steps in such genome reductive evolution, leading towards primary endosimbiosis, and the possibility of complementation or replacement by a secondary symbiont are also indicated. Yet, we need to understand what happens at the beginning of the adaptative process towards an obligate mutualistic relationship. For this purpose, we propose to sequence the complete genome of SOPE, the primary endosymbiont of the rice weevil.

Uncovering rate variation of lateral gene transfer during bacterial genome evolution

BACKGROUND

Large scale genome arrangement, such as whole gene insertion/deletion, plays an important role in bacterial genome evolution. Various methods have been employed to study the dynamic process of gene insertions and deletions, such as parsimony methods and maximum likelihood methods. Previous maximum likelihood studies have assumed that the rate of gene insertions/deletions is constant over different genes. This assumption is unrealistic. For instance, it has been shown that informational genes are less likely to be laterally transferred than non-informational genes. However, how much of the variation in gene transfer rates is due to the difference between informational genes and non-informational genes is unclear. In this study, a Gamma-distribution was incorporated in the likelihood estimation by considering rate variation for gene insertions/deletions between genes. This makes it possible to address whether a difference between informational genes and non-informational genes is the main contributor to rate variation of lateral gene transfers.

RESULTS

The results show that models incorporating rate variation fit the data better than do constant rate models in many phylogenetic groups. Even though informational genes are less likely to be laterally transferred than non-informational genes, the degree of rate variation for insertions/deletions did not change dramatically and remained high even when informational genes were excluded from the study. This suggests that the variation in rate of insertions/deletions is not due mainly to the simple difference between informational genes and non-informational genes. Among genes that are not classified as informational and among the informational genes themselves, there are still large differences in the rates that these genes are inserted and deleted.

CONCLUSION

While the difference in informational gene rates contributes to rate variation, it is only a small fraction of the variation present; instead, a substantial amount of rate variation for insertions/deletions remains among both informational genes and among non-informational genes.

Reconstructing the ancestor of Mycobacterium leprae: the dynamics of gene loss and genome reduction

We have reconstructed the gene content and order of the last common ancestor of the human pathogens Mycobacterium leprae and Mycobacterium tuberculosis. During the reductive evolution of M. leprae, 1537 of 2977 ancestral genes were lost, among which we found 177 previously unnoticed pseudogenes. We find evidence that a massive gene inactivation took place very recently in the M. leprae lineage, leading to the loss of hundreds of ancestral genes. A large proportion of their nucleotide content ( approximately 89%) still remains in the genome, which allowed us to characterize and date them. The age of the pseudogenes was computed using a new methodology based on the rates and patterns of substitution in the pseudogenes and functional orthologous genes of closely related genomes. The position of the genes that were lost in the ancestor's genome revealed that the process of function loss and degradation mainly took place through a gene-to-gene inactivation process, followed by the gradual loss of their DNA. This suggests a scenario of massive genome reduction through many nearly simultaneous pseudogenization events, leading to a highly specialized pathogen.

Reductive genome evolution from the mother of Rickettsia

The Rickettsia genus is a group of obligate intracellular α-proteobacteria representing a paradigm of reductive evolution. Here, we investigate the evolutionary processes that shaped the genomes of the genus. The reconstruction of ancestral genomes indicates that their last common ancestor contained more genes, but already possessed most traits associated with cellular parasitism. The differences in gene repertoires across modern Rickettsia are mainly the result of differential gene losses from the ancestor. We demonstrate using computer simulation that the propensity of loss was variable across genes during this process. We also analyzed the ratio of nonsynonymous to synonymous changes (Ka/Ks) calculated as an average over large sets of genes to assay the strength of selection acting on the genomes of Rickettsia, Anaplasmataceae, and free-living γ-proteobacteria. As a general trend, Ka/Ks were found to decrease with increasing divergence between genomes. The high Ka/Ks for closely related genomes are probably due to a lag in the removal of slightly deleterious nonsynonymous mutations by natural selection. Interestingly, we also observed a decrease of the rate of gene loss with increasing divergence, suggesting a similar lag in the removal of slightly deleterious pseudogene alleles. For larger divergence (Ks > 0.2), Ka/Ks converge toward similar values indicating that the levels of selection are roughly equivalent between intracellular α-proteobacteria and their free-living relatives. This contrasts with the view that obligate endocellular microorganisms tend to evolve faster as a consequence of reduced effectiveness of selection, and suggests a major role of enhanced background mutation rates on the fast protein divergence in the obligate intracellular α-proteobacteria.

Genome downsizing and fast sequence divergence are frequently observed in bacteria living exclusively within the cells of higher eukaryotes. However, the driving forces and contributions of these processes to the genome diversity of the microorganisms remain poorly understood. The genus Rickettsia, a group of small obligate intracellular pathogens of humans, provides a fascinating model to study the genome downsizing process. In this article, we used seven Rickettsia genomes to reconstruct the genome of their ancestor and inferred the origin and fate of the genes found in today's species. We identify the process of gene loss as the main cause of genome diversification within the genus and show that the rate of gene loss, sequence divergence, and genome rearrangements are highly variable across the various Rickettsia lineages. This heterogeneity likely reflects the intricate effects of specialization to distinct arthropod hosts and critical alterations of the gene repertoire, such as the losses of DNA repair genes and the amplification of mobile genes. In contrast, we did not find evidence for the role of reduced population sizes on the long-term acceleration of sequence evolution. Overall, the data presented in this article shed new light on the fundamental evolutionary processes that drive the evolution of obligate intracellular bacteria.

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.

Dynamics of the IncW genetic backbone imply general trends in conjugative plasmid evolution

Plasmids cannot be understood as mere tools for genetic exchange: they are themselves subject to the forces of evolution. Their genomic and phylogenetic features have been less studied in this respect. Focusing on the IncW incompatibility group, which includes the smallest known conjugative plasmids, we attempt to unveil some common trends in plasmid evolution. The functional modules of IncW genetic backbone are described, with emphasis on their architecture and relationships to other plasmid groups. Some plasmid regions exhibit strong phylogenetic mosaicism, in striking contrast to others of unusual synteny conservation. The presence of genes of unknown function that are widely distributed in plasmid genomes is also emphasized, exposing the existence of ill-defined yet conserved plasmid functions. Conjugation is an essential hallmark of IncW plasmid biology and special attention is given to the organization and evolution of its transfer modules. Genetic exchange between plasmids and their hosts is analysed by following the evolution of the type IV secretion system. Adaptation of the trw conjugative machinery to pathogenicity functions in Bartonella is discussed as an example of how plasmids can change their host modus vivendi. Starting from the phage paradigm, our analysis articulates novel concepts that apply to plasmid evolution.

Tempo and mode of early gene loss in endosymbiotic bacteria from insects

BACKGROUND

Understanding evolutionary processes that drive genome reduction requires determining the tempo (rate) and the mode (size and types of deletions) of gene losses. In this study, we analysed five endosymbiotic genome sequences of the gamma-proteobacteria (three different Buchnera aphidicola strains, Wigglesworthia glossinidia, Blochmannia floridanus) to test if gene loss could be driven by the selective importance of genes. We used a parsimony method to reconstruct a minimal ancestral genome of insect endosymbionts and quantified gene loss along the branches of the phylogenetic tree. To evaluate the selective or functional importance of genes, we used a parameter that measures the level of adaptive codon bias in E. coli (i.e. codon adaptive index, or CAI), and also estimates of evolutionary rates (Ka) between pairs of orthologs either in free-living bacteria or in pairs of symbionts.

RESULTS

Our results demonstrate that genes lost in the early stages of symbiosis were on average less selectively constrained than genes conserved in any of the extant symbiotic strains studied. These results also extend to more recent events of gene losses (i.e. among Buchnera strains) that still tend to concentrate on genes with low adaptive bias in E. coli and high evolutionary rates both in free-living and in symbiotic lineages. In addition, we analyzed the physical organization of gene losses for early steps of symbiosis acquisition under the hypothesis of a common origin of different symbioses. In contrast with previous findings we show that gene losses mostly occurred through loss of rather small blocks and mostly in syntenic regions between at least one of the symbionts and present-day E. coli.

CONCLUSION

At both ancient and recent stages of symbiosis evolution, gene loss was at least partially influenced by selection, highly conserved genes being retained more readily than lowly conserved genes: although losses might result from drift due to the bottlenecking of endosymbiontic populations, we demonstrated that purifying selection also acted by retaining genes of greater selective importance.

The fate of laterally transferred genes: life in the fast lane to adaptation or death

Large-scale genome arrangement plays an important role in bacterial genome evolution. A substantial number of genes can be inserted into, deleted from, or rearranged within genomes during evolution. Detecting or inferring gene insertions/deletions is of interest because such information provides insights into bacterial genome evolution and speciation. However, efficient inference of genome events is difficult because genome comparisons alone do not generally supply enough information to distinguish insertions, deletions, and other rearrangements. In this study, homologous genes from the complete genomes of 13 closely related bacteria were examined. The presence or absence of genes from each genome was cataloged, and a maximum likelihood method was used to infer insertion/deletion rates according to the phylogenetic history of the taxa. It was found that whole gene insertions/deletions in genomes occur at rates comparable to or greater than the rate of nucleotide substitution and that higher insertion/deletion rates are often inferred to be present at the tips of the phylogeny with lower rates on more ancient interior branches. Recently transferred genes are under faster and relaxed evolution compared with more ancient genes. Together, this implies that many of the lineage-specific insertions are lost quickly during evolution and that perhaps a few of the genes inserted by lateral transfer are niche specific.

Comparative rates of evolution in endosymbiotic nuclear genomes

BACKGROUND

The nucleomorphs associated with secondary plastids of cryptomonads and chlorarachniophytes are the sole examples of organelles with eukaryotic nuclear genomes. Although not as widespread as their prokaryotic equivalents in mitochondria and plastids, nucleomorph genomes share similarities in terms of reduction and compaction. They also differ in several aspects, not least in that they encode proteins that target to the plastid, and so function in a different compartment from that in which they are encoded.

RESULTS

Here, we test whether the phylogenetically distinct nucleomorph genomes of the cryptomonad, Guillardia theta, and the chlorarachniophyte, Bigelowiella natans, have experienced similar evolutionary pressures during their transformation to reduced organelles. We compared the evolutionary rates of genes from nuclear, nucleomorph, and plastid genomes, all of which encode proteins that function in the same cellular compartment, the plastid, and are thus subject to similar selection pressures. Furthermore, we investigated the divergence of nucleomorphs within cryptomonads by comparing G. theta and Rhodomonas salina.

CONCLUSION

Chlorarachniophyte nucleomorph genes have accumulated errors at a faster rate than other genomes within the same cell, regardless of the compartment where the gene product functions. In contrast, most nucleomorph genes in cryptomonads have evolved faster than genes in other genomes on average, but genes for plastid-targeted proteins are not overly divergent, and it appears that cryptomonad nucleomorphs are not presently evolving rapidly and have therefore stabilized. Overall, these analyses suggest that the forces at work in the two lineages are different, despite the similarities between the structures of their genomes.

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.

Comparative genomics reveals what makes an enterobacterial plant pathogen

The bacterial family Enterobacteriaceae contains some of the most devastating human and animal pathogens, including Escherichia coli, Salmonella enterica and species of Yersinia and Shigella. These are among the best-studied of any organisms, yet there is much to be learned about the nature and evolution of interactions with their hosts and with the wider environment. Comparative and functional genomics have fundamentally improved our understanding of their modes of adaptation to different ecological niches and the genes that determine their pathogenicity. In addition to animal pathogens, Enterobacteriaceae include important plant pathogens, such as Erwinia carotovora subsp. atroseptica (Eca), the first plant-pathogenic enterobacterium to be sequenced. This review focuses on genomic comparisons between Eca and other enterobacteria, with particular emphasis on the differences that exemplify or explain the plant-associated lifestyle(s) of Eca. Horizontal gene transfer in Eca may directly have led to the acquisition of a number of determinants that mediate its interactions, pathogenic or otherwise, with plants, offering a glimpse into its evolutionary divergence from animal-pathogenic enterobacteria.

Inference and analysis of the relative stability of bacterial chromosomes

The stability of genomes is highly variable, both in terms of gene content and gene order. Here I calibrate the loss of gene order conservation (GOC) through time by fitting a simple probabilistic model on pairwise comparisons involving 126 bacterial genomes. The model computes the probability of separation of pairs of contiguous genes per unit of time and fits the data better than previous ones while allowing a mechanistic interpretation for the loss of GOC with time. Although the information on operons is not used in the model, I observe, as expected, that most highly conserved pairs of genes are indeed within operons. However, even the other pairs are much more conserved than expected given the observed experimental rearrangement rates. After 500 Myr, about 50% of the originally contiguous orthologues remain so in the average genome. Hence, the large majority of rearrangements must be deleterious and random genome rearrangements are unlikely to provide for positively selected structural changes. I then use the deviations from the model to define an intrinsic measure of genome stability that allowed the comparison of distantly related genomes and the inference of ancestral states. This shows that clades differ in genome stability, with cyanobacteria being the least stable and gamma-proteobacteria the most stable. Without correction for phylogeny, free-living bacteria are the least stable group of genomes, followed by pathogens, and then endomutualists. However, after correction for phylogenetic inertia (or the removal of cyanobacteria from the analysis), there is no significant association between genome stability and lifestyle or genome size. Hence, although this method has allowed uncovering some of mechanisms leading to rearrangements, we still ignore the forces that differentially shape selection upon genome stability in different species.

Chromosomal stasis versus plasmid plasticity in aphid endosymbiont Buchnera aphidicola

The study of three genomes of the aphid endosymbiont Buchnera aphidicola has revealed an extraordinary stasis: conservation of gene order and genetic composition of the chromosome, while the chromosome size and number of genes has reduced. The reduction in genome size appears to be ongoing since some lineages we now know to have even smaller chromosomes than the first B. aphidicola analysed. The current sequencing by our group of one of these smaller genomes with an estimated size of 450 kb, and its comparison with the other three available genomes provide insights into the nature of processes involved in shrinkage. We discuss whether B. aphidicola might be driven to extinction and be replaced by secondary aphid endosymbionts. In some lineages, genes encoding key enzymes in the pathways leading to tryptophan and leucine biosynthesis (trpEG and leuABCD, respectively) are located on plasmids, rather than the main chromosome. In contrast to the stasis of the main chromosome, plasmid genes have frequently been transferred to the main chromosome and undergone other gene rearrangements. We propose a two-step scenario to explain these contrasting modes of evolution. Essential genes may have escaped regulation by moving to plasmids in a moving B. aphidicola ancestor. B. aphidicola became polyploidy at a given stage of its evolution and plasmid genes have been transferred to the main chromosome through several independent events.