Mutation rate and genome reduction in endosymbiotic and free-living bacteria

A new lysine biosynthetic enzyme from a bacterial endosymbiont shaped by genetic drift and genome reduction

The effect of population bottlenecks and genome reduction on enzyme function is poorly understood. Candidatus Liberibacter solanacearum is a bacterium with a reduced genome that is transmitted vertically to the egg of an infected psyllid-a population bottleneck that imposes genetic drift and is predicted to affect protein structure and function. Here, we define the function of Ca. L. solanacearum dihydrodipicolinate synthase (CLsoDHDPS), which catalyzes the committed branchpoint reaction in diaminopimelate and lysine biosynthesis. We demonstrate that CLsoDHDPS is expressed in Ca. L. solanacearum and expression is increased ~2-fold in the insect host compared to in planta. CLsoDHDPS has decreased thermal stability and increased aggregation propensity, implying mutations have destabilized the enzyme but are compensated for through elevated chaperone expression and a stabilized oligomeric state. CLsoDHDPS uses a ternary-complex kinetic mechanism, which is to date unique among DHDPS enzymes, has unusually low catalytic ability, but an unusually high substrate affinity. Structural studies demonstrate that the active site is more open, and the structure of CLsoDHDPS with both pyruvate and the substrate analogue succinic-semialdehyde reveals that the product is both structurally and energetically different and therefore evolution has in this case fashioned a new enzyme. Our study suggests the effects of genome reduction and genetic drift on the function of essential enzymes and provides insights on bacteria-host co-evolutionary associations. We propose that bacteria with endosymbiotic lifestyles present a rich vein of interesting enzymes useful for understanding enzyme function and/or informing protein engineering efforts.

Genome reduction occurred in early Prochlorococcus with an unusually low effective population size

In the oligotrophic sunlit ocean, the most abundant free-living planktonic bacterial lineages evolve convergently through genome reduction. The cyanobacterium Prochlorococcus responsible for 10% global oxygen production is a prominent example. The dominant theory known as "genome streamlining" posits that they have extremely large effective population sizes (Ne) such that selection for metabolic efficiency acts to drive genome reduction. Because genome reduction largely took place anciently, this theory builds on the assumption that their ancestors' Ne was similarly large. Constraining Ne for ancient ancestors is challenging because experimental measurements of extinct organisms are impossible and alternatively reconstructing ancestral Ne with phylogenetic models gives large uncertainties. Here, we develop a new strategy that leverages agent-based modeling to simulate the changes in the genome-wide ratio of radical to conservative nonsynonymous nucleotide substitution rate (dR/dC) in a possible range of Ne in ancestral populations. This proxy shows expected increases with decreases of Ne only when Ne falls to about 10 k - 100 k or lower, magnitudes characteristic of Ne of obligate endosymbiont species where drift drives genome reduction. Our simulations therefore strongly support a scenario where the primary force of Prochlorococcus genome reduction is drift rather than selection.

Pan-genome analysis of six Paracoccus type strain genomes reveal lifestyle traits

The genus Paracoccus capable of inhabiting a variety of different ecological niches both, marine and terrestrial, is globally distributed. In addition, Paracoccus is taxonomically, metabolically and regarding lifestyle highly diverse. Until now, little is known on how Paracoccus can adapt to such a range of different ecological niches and lifestyles. In the present study, the genus Paracoccus was phylogenomically analyzed (n = 160) and revisited, allowing species level classification of 16 so far unclassified Paracoccus sp. strains and detection of five misclassifications. Moreover, we performed pan-genome analysis of Paracoccus-type strains, isolated from a variety of ecological niches, including different soils, tidal flat sediment, host association such as the bluespotted cornetfish, Bugula plumosa, and the reef-building coral Stylophora pistillata to elucidate either i) the importance of lifestyle and adaptation potential, and ii) the role of the genomic equipment and niche adaptation potential. Six complete genomes were de novo hybrid assembled using a combination of short and long-read technologies. These Paracoccus genomes increase the number of completely closed high-quality genomes of type strains from 15 to 21. Pan-genome analysis revealed an open pan-genome composed of 13,819 genes with a minimal chromosomal core (8.84%) highlighting the genomic adaptation potential and the huge impact of extra-chromosomal elements. All genomes are shaped by the acquisition of various mobile genetic elements including genomic islands, prophages, transposases, and insertion sequences emphasizing their genomic plasticity. In terms of lifestyle, each mobile genetic elements should be evaluated separately with respect to the ecological context. Free-living genomes, in contrast to host-associated, tend to comprise (1) larger genomes, or the highest number of extra-chromosomal elements, (2) higher number of genomic islands and insertion sequence elements, and (3) a lower number of intact prophage regions. Regarding lifestyle adaptations, free-living genomes share genes linked to genetic exchange via T4SS, especially relevant for Paracoccus, known for their numerous extrachromosomal elements, enabling adaptation to dynamic environments. Conversely, host-associated genomes feature diverse genes involved in molecule transport, cell wall modification, attachment, stress protection, DNA repair, carbon, and nitrogen metabolism. Due to the vast number of adaptive genes, Paracoccus can quickly adapt to changing environmental conditions.

Single-amplified genomes reveal most streamlined free-living marine bacteria

Evolutionary adaptations of prokaryotes to the environment sometimes result in genome reduction. Our knowledge of this phenomenon among free-living bacteria remains scarce. We address the dynamics and limits of genome reduction by examining one of the most abundant bacteria in the ocean, the SAR86 clade. Despite its abundance, comparative genomics has been limited by the absence of pure cultures and the poor representation in metagenome-assembled genomes. We co-assembled multiple previously available single-amplified genomes to obtain the first complete genomes from members of the four families. All families showed a convergent evolutionary trajectory with characteristic features of streamlined genomes, most pronounced in the TMED112 family. This family has a genome size of ca. 1 Mb and only 1 bp as median intergenic distance, exceeding values found in other abundant microbes such as SAR11, OM43 and Prochlorococcus. This genomic simplification led to a reduction in the biosynthesis of essential molecules, DNA repair-related genes, and the ability to sense and respond to environmental factors, which could suggest an evolutionary dependence on other co-occurring microbes for survival (Black Queen hypothesis). Therefore, these reconstructed genomes within the SAR86 clade provide new insights into the limits of genome reduction in free-living marine bacteria.

Identification of shared bacterial strains in the vaginal microbiota of related and unrelated reproductive-age mothers and daughters using genome-resolved metagenomics

It has been suggested that the human microbiome might be vertically transmitted from mother to offspring and that early colonizers may play a critical role in development of the immune system. Studies have shown limited support for the vertical transmission of the intestinal microbiota but the derivation of the vaginal microbiota remains largely unknown. Although the vaginal microbiota of children and reproductive age women differ in composition, the vaginal microbiota could be vertically transmitted. To determine whether there was any support for this hypothesis, we examined the vaginal microbiota of daughter-mother pairs from the Baltimore metropolitan area (ages 14-27, 32-51; n = 39). We assessed whether the daughter's microbiota was similar in composition to their mother's using metataxonomics. Permutation tests revealed that while some pairs did have similar vaginal microbiota, the degree of similarity did not exceed that expected by chance. Genome-resolved metagenomics was used to identify shared bacterial strains in a subset of the families (n = 22). We found a small number of bacterial strains that were shared between mother-daughter pairs but identified more shared strains between individuals from different families, indicating that vaginal bacteria may display biogeographic patterns. Earlier-in-life studies are needed to demonstrate vertical transmission of the vaginal microbiota.

Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales

The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel’s lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.

The evolutionary forces driving genome size in bacteria and archaea have been subject to debate during the last decades. Typically, independent comparative analyses have suggested that unique variables, such as the strength of selection, environmental complexity, and mutation rate, are the main drivers of this trait, without considering for potential biases derived from shared ancestry. Here, we applied a phylogeny-based statistical approach to assess how tightly genome size in bacteria and archaea is linked to evolutionary history. Moreover, we also evaluated the predictability of genome size from the strength of purifying selection and ecological strategy on a broad diversity of bacteria and archaea genomes under a phylogenetic comparative framework. Our approach indicates that despite the ability of bacteria and archaea to rapidly exchange genes, a strong phylogenetic signal to genome size distributions can be recovered at broad phylogenetic scales.

Prochlorococcus have low global mutation rate and small effective population size

Prochlorococcus are the most abundant free-living photosynthetic carbon-fixing organisms in the ocean. Prochlorococcus show small genome sizes, low genomic G+C content, reduced DNA repair gene pool and fast evolutionary rates, which are typical features of endosymbiotic bacteria. Nevertheless, their evolutionary mechanisms are believed to be different. Evolution of endosymbiotic bacteria is dominated by genetic drift owing to repeated population bottlenecks, whereas Prochlorococcus are postulated to have extremely large effective population sizes (N_e) and thus drift has rarely been considered. However, accurately extrapolating N_e requires measuring an unbiased global mutation rate through mutation accumulation, which is challenging for Prochlorococcus. Here, we managed this experiment over 1,065 days using Prochlorococcus marinus AS9601, sequenced genomes of 141 mutant lines and determined its mutation rate to be 3.50 × 10^-10 per site per generation. Extrapolating N_e additionally requires identifying population boundaries, which we defined using PopCOGenT and over 400 genomes related to AS9601. Accordingly, we calculated its N_e to be 1.68 × 10⁷, which is only reasonably greater than that of endosymbiotic bacteria but surprisingly smaller than that of many free-living bacteria extrapolated using the same approach. Our results therefore suggest that genetic drift is a key driver of Prochlorococcus evolution.

Enhanced Mutation Rate, Relaxed Selection, and the "Domino Effect" are associated with Gene Loss in Blattabacterium, A Cockroach Endosymbiont

Intracellular endosymbionts have reduced genomes that progressively lose genes at a timescale of tens of million years. We previously reported that gene loss rate is linked to mutation rate in Blattabacterium, however, the mechanisms causing gene loss are not yet fully understood. Here, we carried out comparative genomic analyses on the complete genome sequences of a representative set of 67 Blattabacterium strains, with sizes ranging between 511 and 645 kb. We found that 200 of the 566 analyzed protein-coding genes were lost in at least one lineage of Blattabacterium, with the most extreme case being one gene that was lost independently in 24 lineages. We found evidence for three mechanisms influencing gene loss in Blattabacterium. First, gene loss rates were found to increase exponentially with the accumulation of substitutions. Second, genes involved in vitamin and amino acid metabolism experienced relaxed selection in Cryptocercus and Mastotermes, possibly triggered by their vertically inherited gut symbionts. Third, we found evidence of epistatic interactions among genes leading to a "domino effect" of gene loss within pathways. Our results highlight the complexity of the process of genome erosion in an endosymbiont.

Defining life and evolution: Essay on the origin, expansion, and evolution of living matter

This essay aims to define the origin, expansion, and evolution of living matter. The first formations, identified as remains, fossils, traces etc. of life are almost as old as the Earth itself. During four billion years, life on the Earth has continuously existed and been implemented in the range of conditions, ensuring the liquid state of water. During the entire period of life existence, its evolution was proceeding with the tendency of multidirectionality, after each catastrophe tending to the diversity and vastness of distribution, and all the currently living species, regardless of their complexity, have the same evolutionary age. The property of reproductive surplus (multiplication) is inherent in all the living matter. The reproduction of all the living matter is implemented via the "development" - a process of continuous occurrence of something new that did not exist in the previous moment in the reproduced individual at each specific moment of time with the tendency towards the reproduction of a "copy". In its fundamental basis, Life is based on a programme, its material support is implemented and exists not in the field of causative-consecutive events, but in the field of programmed-causative-consecutive events. This predetermines the "biology laws", the behaviour of the material constituent of Life at each time period, and the future of the material constituent of life.

Mechanisms driving genome reduction of a novel Roseobacter lineage

Members of the marine Roseobacter group are key players in the global carbon and sulfur cycles. While over 300 species have been described, only 2% possess reduced genomes (mostly 3-3.5 Mbp) compared to an average roseobacter (>4 Mbp). These taxonomic minorities are phylogenetically diverse but form a Pelagic Roseobacter Cluster (PRC) at the genome content level. Here, we cultivated eight isolates constituting a novel Roseobacter lineage which we named 'CHUG'. Metagenomic and metatranscriptomic read recruitment analyses showed that CHUG members are globally distributed and active in marine pelagic environments. CHUG members possess some of the smallest genomes (~2.6 Mb) among all known roseobacters, but they do not exhibit canonical features of typical bacterioplankton lineages theorized to have undergone genome streamlining processes, like higher coding density, fewer paralogues and rarer pseudogenes. While CHUG members form a genome content cluster with traditional PRC members, they show important differences. Unlike other PRC members, neither the relative abundances of CHUG members nor their relative gene expression levels are correlated with chlorophyll a concentration across the global samples. CHUG members cannot utilize most phytoplankton-derived metabolites or synthesize vitamin B₁₂, a key metabolite mediating the roseobacter-phytoplankton interactions. This combination of features is evidence for the hypothesis that CHUG members may have evolved a free-living lifestyle decoupled from phytoplankton. This ecological transition was accompanied by the loss of signature genes involved in roseobacter-phytoplankton symbiosis, suggesting that relaxation of purifying selection owing to lifestyle shift is likely an important driver of genome reduction in CHUG.

Genome Evolution: Mutation Is the Main Driver of Genome Size in Prokaryotes

Despite intense research on genome architecture since the 2000s, genome-size evolution in prokaryotes has remained puzzling. Using a phylogenetic approach, a new study found that increased mutation rate is associated with gene loss and reduced genome size in prokaryotes.

The microbiome: Composition and locations

The human body is home to a diverse and functionally important assemblage of symbiotic microbes that varies predictably over different spatial scales, both within and across body sites. The composition of these spatially distinct microbial consortia can be impacted by a variety of stochastic and deterministic forces, including dispersal from different source communities, and selection by regionally-specific host processes for the enrichment of physiologically significant taxa. In this chapter, we review the composition, function, and assembly of the healthy human gastrointestinal, skin, vaginal, and respiratory microbiomes, with special emphasis on the regional distribution of microbes throughout the gastrointestinal tract.

Ancestral niche separation and evolutionary rate differentiation between sister marine flavobacteria lineages

Marine flavobacteria are specialists for polysaccharide degradation. They dominate in habitats enriched with polysaccharides, but are also prevalent in pelagic environments where polysaccharides are less available. These niches are likely occupied by distinct lineages, but evolutionary processes underlying their niche differentiation remain elusive. Here, genomic analyses and physiological assays indicate that the sister flavobacteria lineages Leeuwenhoekiella and Nonlabens likely explore polysaccharide-rich macroalgae and polysaccharide-poor pelagic niches respectively. Phylogenomic analyses inferred that the niche separation likely occurred anciently and coincided with increased sequence evolutionary rate in Nonlabens compared with Leeuwenhoekiella. Further analyses ruled out the known mechanisms likely driving evolutionary rate acceleration, including reduced selection efficiency, decreased generation time and increased mutation rate. In particular, the mutation rates were determined using an unbiased experimental method, which measures the present-day populations and may not reflect ancestral populations. These data collectively lead to a new hypothesis that an ancestral and transient mutation rate increase resulted in evolutionary rate increase in Nonlabens. This hypothesis was supported by inferring that gains and losses of genes involved in SOS response, a mechanism known to drive transiently increased mutation rate, coincided with evolutionary rate acceleration. Our analyses highlight the evolutionary mechanisms underlying niche differentiation of flavobacteria lineages.

Strong Purifying Selection Is Associated with Genome Streamlining in Epipelagic Marinimicrobia

Marine microorganisms inhabiting nutrient-depleted waters play critical roles in global biogeochemical cycles due to their abundance and broad distribution. Many of these microbes share similar genomic features including small genome size, low % G + C content, short intergenic regions, and low nitrogen content in encoded amino acid residue side chains (N-ARSC), but the evolutionary drivers of these characteristics are unclear. Here, we compared the strength of purifying selection across the Marinimicrobia, a candidate phylum which encompasses a broad range of phylogenetic groups with disparate genomic features, by estimating the ratio of nonsynonymous and synonymous substitutions (dN/dS) in conserved marker genes. Our analysis reveals that epipelagic Marinimicrobia that exhibit features consistent with genome streamlining have significantly lower dN/dS values when compared with their mesopelagic counterparts. We also found a significant positive correlation between median dN/dS values and % G + C content, N-ARSC, and intergenic region length. We did not identify a significant correlation between dN/dS ratios and estimated genome size, suggesting the strength of selection is not a primary factor shaping genome size in this group. Our findings are generally consistent with genome streamlining theory, which postulates that many genomic features of abundant epipelagic bacteria are the result of adaptation to oligotrophic nutrient conditions. Our results are also in agreement with previous findings that genome streamlining is common in epipelagic waters, suggesting that microbes inhabiting this region of the ocean have been shaped by strong selection together with prevalent nutritional constraints characteristic of this environment.

Adapting the engine to the fuel: mutator populations can reduce the mutational load by reorganizing their genome structure

Background

Mutators are common in bacterial populations, both in natural isolates and in the lab. The fate of these lineages, which mutation rate is increased up to 100 ×, has long been studied using population genetics models, showing that they can spread in a population following an environmental change. However in stable conditions, they suffer from the increased mutational load, hence being overcome by non-mutators. However, these results don’t take into account the fact that an elevated mutation rate can impact the genetic structure, hence changing the sensitivity of the population to mutations. Here we used Aevol, an in silico experimental evolution platform in which genomic structures are free to evolve, in order to study the fate of mutator populations evolving for a long time in constant conditions.

Results

Starting from wild-types that were pre-evolved for 300,000 generations, we let 100 mutator populations (point mutation rate ×100) evolve for 100,000 further generations in constant conditions. As expected all populations initially undergo a fitness loss. However, after that the mutator populations started to recover. Most populations ultimately recovered their ancestors fitness, and a significant fraction became even fitter than the non-mutator control clones that evolved in parallel. By analyzing the genomes of the mutators, we show that the fitness recovery is due to two mechanisms: i. an increase in robustness through compaction of the coding part of the mutator genomes, ii. an increase of the selection coefficient that decreases the mean-fitness of the population. Strikingly the latter is due to the accumulation of non-coding sequences in the mutators genomes.

Conclusion

Our results show that the mutational burden that is classically thought to be associated with mutator phenotype is escapable. On the long run mutators adapted their genomes and reshaped the distribution of mutation effects. Therewith the lineage is able to recover fitness even though the population still suffers the elevated mutation rate. Overall these results change our view of mutator dynamics: by being able to reduce the deleterious effect of the elevated mutation rate, mutator populations may be able to last for a very long time; A situation commonly observed in nature.

Genome Rearrangement Shapes Prochlorococcus Ecological Adaptation

Prochlorococcus, the most abundant and smallest known free-living photosynthetic microorganism, plays a key role in marine ecosystems and biogeochemical cycles. Prochlorococcus genome evolution is a fundamental issue related to how Prochlorococcus clades adapted to different ecological niches. Recent studies revealed that the gene gain and loss is crucial to the clade differentiation. The significance of our research is that we interpreted the Prochlorococcus genome evolution from the perspective of genome structure and associated the genome rearrangement with the Prochlorococcus clade differentiation and subsequent ecological adaptation.

Prochlorococcus is the most abundant and smallest known free-living photosynthetic microorganism and is a key player in marine ecosystems and biogeochemical cycles. Prochlorococcus can be broadly divided into high-light-adapted (HL) and low-light-adapted (LL) clades. In this study, we isolated two low-light-adapted clade I (LLI) strains from the western Pacific Ocean and obtained their genomic data. We reconstructed Prochlorococcus evolution based on genome rearrangement. Our results showed that genome rearrangement might have played an important role in Prochlorococcus evolution. We also found that the Prochlorococcus clades with streamlined genomes maintained relatively high synteny throughout most of their genomes, and several regions served as rearrangement hotspots. Backbone analysis showed that different clades shared a conserved backbone but also had clade-specific regions, and the genes in these regions were associated with ecological adaptations.

IMPORTANCEProchlorococcus, the most abundant and smallest known free-living photosynthetic microorganism, plays a key role in marine ecosystems and biogeochemical cycles. Prochlorococcus genome evolution is a fundamental issue related to how Prochlorococcus clades adapted to different ecological niches. Recent studies revealed that the gene gain and loss is crucial to the clade differentiation. The significance of our research is that we interpreted the Prochlorococcus genome evolution from the perspective of genome structure and associated the genome rearrangement with the Prochlorococcus clade differentiation and subsequent ecological adaptation.

Comparative Genomic Analysis of Holospora spp., Intranuclear Symbionts of Paramecia

While most endosymbiotic bacteria are transmitted only vertically, Holospora spp., an alphaproteobacterium from the Rickettsiales order, can desert its host and invade a new one. All bacteria from the genus Holospora are intranuclear symbionts of ciliates Paramecium spp. with strict species and nuclear specificity. Comparative metabolic reconstruction based on the newly sequenced genome of Holospora curviuscula, a macronuclear symbiont of Paramecium bursaria, and known genomes of other Holospora species shows that even though all Holospora spp. can persist outside the host, they cannot synthesize most of the essential small molecules, such as amino acids, and lack some central energy metabolic pathways, including glycolysis and the citric acid cycle. As the main energy source, Holospora spp. likely rely on nucleotides pirated from the host. Holospora-specific genes absent from other Rickettsiales are possibly involved in the lifestyle switch from the infectious to the reproductive form and in cell invasion.

Why Might Bacterial Pathogens Have Small Genomes?

Bacteria that cause serious disease often have smaller genomes, and fewer genes, than their nonpathogenic, or less pathogenic relatives. Here, we review evidence for the generality of this association, and summarise the various reasons why the association might hold. We focus on the population genetic processes that might lead to reductive genome evolution, and show how several of these could be connected to pathogenicity. We find some evidence for most of the processes having acted in bacterial pathogens, including several different modes of genome reduction acting in the same lineage. We argue that predictable processes of genome evolution might not reflect any common underlying process.

When the ribosomal DNA does not tell the truth: The case of the taxonomic position of Kurtia argillacea, an ericoid mycorrhizal fungus residing among Hymenochaetales

The nuclear ribosomal DNA (nuc-rDNA) is widely used for the identification and phylogenetic reconstruction of Agaricomycetes. However, nuc-rDNA-based phylogenies may sometimes be in conflict with phylogenetic relationships derived from protein coding genes. In this study, the taxonomic position of the basidiomycetous mycobiont that forms the recently discovered sheathed ericoid mycorrhiza was investigated, because its nuc-rDNA is highly dissimilar to any other available fungal sequences in terms of nucleotide composition and length, and its nuc-rDNA-based phylogeny is inconclusive and significantly disagrees with protein coding sequences and morphological data. In the present work, this mycobiont was identified as Kurtia argillacea (= Hyphoderma argillaceum) residing in the order Hymenochaetales (Basidiomycota). Bioinformatic screening of the Kurtia ribosomal DNA sequence indicates that it represents a gene with a non-standard substitution rate or nucleotide composition heterogeneity rather than a deep paralogue or a pseudogene. Such a phenomenon probably also occurs in other lineages of the Fungi and should be taken into consideration when nuc-rDNA (especially that with unusual nucleotide composition) is used as a sole marker for phylogenetic reconstructions. Kurtia argillacea so far represents the only confirmed non-sebacinoid ericoid mycorrhizal fungus in the Basidiomycota and its intriguing placement among mostly saprobic and parasitic Hymenochaetales begs further investigation of its eco-physiology.

Mutator genomes decay, despite sustained fitness gains, in a long-term experiment with bacteria

Understanding the extreme variation among bacterial genomes remains an unsolved challenge in evolutionary biology, despite long-standing debate about the relative importance of natural selection, mutation, and random drift. A potentially important confounding factor is the variation in mutation rates between lineages and over evolutionary history, which has been documented in several species. Mutation accumulation experiments have shown that hypermutability can erode genomes over short timescales. These results, however, were obtained under conditions of extremely weak selection, casting doubt on their general relevance. Here, we circumvent this limitation by analyzing genomes from mutator populations that arose during a long-term experiment with Escherichia coli, in which populations have been adaptively evolving for >50,000 generations. We develop an analytical framework to quantify the relative contributions of mutation and selection in shaping genomic characteristics, and we validate it using genomes evolved under regimes of high mutation rates with weak selection (mutation accumulation experiments) and low mutation rates with strong selection (natural isolates). Our results show that, despite sustained adaptive evolution in the long-term experiment, the signature of selection is much weaker than that of mutational biases in mutator genomes. This finding suggests that relatively brief periods of hypermutability can play an outsized role in shaping extant bacterial genomes. Overall, these results highlight the importance of genomic draft, in which strong linkage limits the ability of selection to purge deleterious mutations. These insights are also relevant to other biological systems evolving under strong linkage and high mutation rates, including viruses and cancer cells.

Excess of non-conservative amino acid changes in marine bacterioplankton lineages with reduced genomes

Surface ocean waters are dominated by planktonic bacterial lineages with highly reduced genomes. The best examples are the cyanobacterial genus Prochlorococcus, the alphaproteobacterial clade SAR11 and the gammaproteobacterial clade SAR86, which together represent over 50% of the cells in surface oceans. Several studies have identified signatures of selection on these lineages in today's ocean and have postulated selection as the primary force throughout their evolutionary history. However, massive loss of genomic DNA in these lineages often occurred in the distant past, and the selective pressures underlying these ancient events have not been assessed. Here, we probe ancient selective pressures by computing %GC-corrected rates of conservative and radical nonsynonymous nucleotide substitutions. Surprisingly, we found an excess of radical changes in several of these lineages in comparison to their relatives with larger genomes. Furthermore, analyses of allelic genome sequences of several populations within these lineages consistently supported that radical replacements are more likely to be deleterious than conservative changes. Our results suggest coincidence of massive genomic DNA losses and increased power of genetic drift, but we also suggest that additional evidence independent of the nucleotide substitution analyses is needed to support a primary role of genetic drift driving ancient genome reduction of marine bacterioplankton lineages.

Natural selection drove metabolic specialization of the chromatophore in Paulinella chromatophora

Background

Genome degradation of host-restricted mutualistic endosymbionts has been attributed to inactivating mutations and genetic drift while genes coding for host-relevant functions are conserved by purifying selection. Unlike their free-living relatives, the metabolism of mutualistic endosymbionts and endosymbiont-originated organelles is specialized in the production of metabolites which are released to the host. This specialization suggests that natural selection crafted these metabolic adaptations. In this work, we analyzed the evolution of the metabolism of the chromatophore of Paulinella chromatophora by in silico modeling. We asked whether genome reduction is driven by metabolic engineering strategies resulted from the interaction with the host. As its widely known, the loss of enzyme coding genes leads to metabolic network restructuring sometimes improving the production rates. In this case, the production rate of reduced-carbon in the metabolism of the chromatophore.

Results

We reconstructed the metabolic networks of the chromatophore of P. chromatophora CCAC 0185 and a close free-living relative, the cyanobacterium Synechococcus sp. WH 5701. We found that the evolution of free-living to host-restricted lifestyle rendered a fragile metabolic network where >80% of genes in the chromatophore are essential for metabolic functionality. Despite the lack of experimental information, the metabolic reconstruction of the chromatophore suggests that the host provides several metabolites to the endosymbiont. By using these metabolites as intracellular conditions, in silico simulations of genome evolution by gene lose recover with 77% accuracy the actual metabolic gene content of the chromatophore. Also, the metabolic model of the chromatophore allowed us to predict by flux balance analysis a maximum rate of reduced-carbon released by the endosymbiont to the host. By inspecting the central metabolism of the chromatophore and the free-living cyanobacteria we found that by improvements in the gluconeogenic pathway the metabolism of the endosymbiont uses more efficiently the carbon source for reduced-carbon production. In addition, our in silico simulations of the evolutionary process leading to the reduced metabolic network of the chromatophore showed that the predicted rate of released reduced-carbon is obtained in less than 5% of the times under a process guided by random gene deletion and genetic drift. We interpret previous findings as evidence that natural selection at holobiont level shaped the rate at which reduced-carbon is exported to the host. Finally, our model also predicts that the ABC phosphate transporter (pstSACB) which is conserved in the genome of the chromatophore of P. chromatophora strain CCAC 0185 is a necessary component to release reduced-carbon molecules to the host.

Conclusion

Our evolutionary analysis suggests that in the case of Paulinella chromatophora natural selection at the holobiont level played a prominent role in shaping the metabolic specialization of the chromatophore. We propose that natural selection acted as a “metabolic engineer” by favoring metabolic restructurings that led to an increased release of reduced-carbon to the host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0947-6) contains supplementary material, which is available to authorized users.

Gene Loss and Horizontal Gene Transfer Contributed to the Genome Evolution of the Extreme Acidophile "Ferrovum"

Acid mine drainage (AMD), associated with active and abandoned mining sites, is a habitat for acidophilic microorganisms that gain energy from the oxidation of reduced sulfur compounds and ferrous iron and that thrive at pH below 4. Members of the recently proposed genus “Ferrovum” are the first acidophilic iron oxidizers to be described within the Betaproteobacteria. Although they have been detected as typical community members in AMD habitats worldwide, knowledge of their phylogenetic and metabolic diversity is scarce. Genomics approaches appear to be most promising in addressing this lacuna since isolation and cultivation of “Ferrovum” has proven to be extremely difficult and has so far only been successful for the designated type strain “Ferrovum myxofaciens” P3G. In this study, the genomes of two novel strains of “Ferrovum” (PN-J185 and Z-31) derived from water samples of a mine water treatment plant were sequenced. These genomes were compared with those of “Ferrovum” sp. JA12 that also originated from the mine water treatment plant, and of the type strain (P3G). Phylogenomic scrutiny suggests that the four strains represent three “Ferrovum” species that cluster in two groups (1 and 2). Comprehensive analysis of their predicted metabolic pathways revealed that these groups harbor characteristic metabolic profiles, notably with respect to motility, chemotaxis, nitrogen metabolism, biofilm formation and their potential strategies to cope with the acidic environment. For example, while the “F. myxofaciens” strains (group 1) appear to be motile and diazotrophic, the non-motile group 2 strains have the predicted potential to use a greater variety of fixed nitrogen sources. Furthermore, analysis of their genome synteny provides first insights into their genome evolution, suggesting that horizontal gene transfer and genome reduction in the group 2 strains by loss of genes encoding complete metabolic pathways or physiological features contributed to the observed diversification.

DNA Repair Is Associated with Information Content in Bacteria, Archaea, and DNA Viruses

The concept of a "proteomic constraint" proposes that DNA repair capacity is positively correlated with the information content of a genome, which can be approximated to the size of the proteome (P). This in turn implies that DNA repair genes are more likely to be present in genomes with larger values of P. This stands in contrast to the common assumption that informational genes have a core function and so are evenly distributed across organisms. We examined the presence/absence of 18 DNA repair genes in bacterial genomes. A positive relationship between gene presence and P was observed for 17 genes in the total dataset, and 16 genes when only nonintracellular bacteria were examined. A marked reduction of DNA repair genes was observed in intracellular bacteria, consistent with their reduced value of P. We also examined archaeal and DNA virus genomes, and show that the presence of DNA repair genes is likewise related to a larger value of P. In addition, the products of the bacterial genes mutY, vsr, and ndk, involved in the correction of GC/AT mutations, are strongly associated with reduced genome GC content. We therefore propose that a reduction in information content leads to a loss of DNA repair genes and indirectly to a reduction in genome GC content in bacteria by exposure to the underlying AT mutation bias. The reduction in P may also indirectly lead to the increase in substitution rates observed in intracellular bacteria via loss of DNA repair genes.

The evolution of sex: A new hypothesis based on mitochondrial mutational erosion: Mitochondrial mutational erosion in ancestral eukaryotes would favor the evolution of sex, harnessing nuclear recombination to optimize compensatory nuclear coadaptation

The evolution of sex in eukaryotes represents a paradox, given the "twofold" fitness cost it incurs. We hypothesize that the mutational dynamics of the mitochondrial genome would have favored the evolution of sexual reproduction. Mitochondrial DNA (mtDNA) exhibits a high-mutation rate across most eukaryote taxa, and several lines of evidence suggest that this high rate is an ancestral character. This seems inexplicable given that mtDNA-encoded genes underlie the expression of life's most salient functions, including energy conversion. We propose that negative metabolic effects linked to mitochondrial mutation accumulation would have invoked selection for sexual recombination between divergent host nuclear genomes in early eukaryote lineages. This would provide a mechanism by which recombinant host genotypes could be rapidly shuffled and screened for the presence of compensatory modifiers that offset mtDNA-induced harm. Under this hypothesis, recombination provides the genetic variation necessary for compensatory nuclear coadaptation to keep pace with mitochondrial mutation accumulation.

Strong genome-wide selection early in the evolution of Prochlorococcus resulted in a reduced genome through the loss of a large number of small effect genes

The smallest genomes of any photosynthetic organisms are found in a group of free-living marine cyanobacteria, Prochlorococcus. To determine the underlying evolutionary mechanisms, we developed a new method to reconstruct the steps leading to the Prochlorococcus genome reduction using 12 Prochlorococcus and 6 marine Synechococcus genomes. Our results reveal that small genome sizes within Prochlorococcus were largely determined shortly after the split of Prochlorococcus and Synechococcus (an early big shrink) and thus for the first time decouple the genome reduction from Prochlorococcus diversification. A maximum likelihood approach was then used to estimate changes of nucleotide substitution rate and selection strength along Prochlorococcus evolution in a phylogenetic framework. Strong genome wide purifying selection was associated with the loss of many genes in the early evolutionary stage. The deleted genes were distributed around the genome, participated in many different functional categories and in general had been under relaxed selection pressure. We propose that shortly after Prochlorococcus diverged from its common ancestor with marine Synechococcus, its population size increased quickly thus increasing efficacy of selection. Due to limited nutrients and a relatively constant environment, selection favored a streamlined genome for maximum economy. Strong genome wide selection subsequently caused the loss of genes with small functional effect including the loss of some DNA repair genes. In summary, genome reduction in Prochlorococcus resulted in genome features that are similar to symbiotic bacteria and pathogens, however, the small genome sizes resulted from an increase in genome wide selection rather than a consequence of a reduced ecological niche or relaxed selection due to genetic drift.

The transcriptome landscape of Prochlorococcus MED4 and the factors for stabilizing the core genome

BACKGROUND

Gene gain and loss frequently occurs in the cyanobacterium Prochlorococcus, a phototroph that numerically dominates tropical and subtropical open oceans. However, little is known about the stabilization of its core genome, which contains approximately 1250 genes, in the context of genome streamlining. Using Prochlorococcus MED4 as a model organism, we investigated the constraints on core genome stabilization using transcriptome profiling.

RESULTS

RNA-Seq technique was used to obtain the transcriptome map of Prochlorococcus MED4, including operons, untranslated regions, non-coding RNAs, and novel genes. Genome-wide expression profiles revealed that three factors contribute to core genome stabilization. First, a negative correlation between gene expression levels and protein evolutionary rates was observed. Highly expressed genes were overrepresented in the core genome but not in the flexible genome. Gene necessity was determined as a second powerful constraint on genome evolution through functional enrichment analysis. Third, quick mRNA turnover may increase corresponding proteins' fidelity among genes that were abundantly expressed. Together, these factors influence core genome stabilization during MED4 genome evolution.

CONCLUSIONS

Gene expression, gene necessity, and mRNA turnover contribute to core genome maintenance during cyanobacterium Prochlorococcus genus evolution.

Evolution of small prokaryotic genomes

As revealed by genome sequencing, the biology of prokaryotes with reduced genomes is strikingly diverse. These include free-living prokaryotes with ∼800 genes as well as endosymbiotic bacteria with as few as ∼140 genes. Comparative genomics is revealing the evolutionary mechanisms that led to these small genomes. In the case of free-living prokaryotes, natural selection directly favored genome reduction, while in the case of endosymbiotic prokaryotes neutral processes played a more prominent role. However, new experimental data suggest that selective processes may be at operation as well for endosymbiotic prokaryotes at least during the first stages of genome reduction. Endosymbiotic prokaryotes have evolved diverse strategies for living with reduced gene sets inside a host-defined medium. These include utilization of host-encoded functions (some of them coded by genes acquired by gene transfer from the endosymbiont and/or other bacteria); metabolic complementation between co-symbionts; and forming consortiums with other bacteria within the host. Recent genome sequencing projects of intracellular mutualistic bacteria showed that previously believed universal evolutionary trends like reduced G+C content and conservation of genome synteny are not always present in highly reduced genomes. Finally, the simplified molecular machinery of some of these organisms with small genomes may be used to aid in the design of artificial minimal cells. Here we review recent genomic discoveries of the biology of prokaryotes endowed with small gene sets and discuss the evolutionary mechanisms that have been proposed to explain their peculiar nature.

Ecological and evolutionary genomics of marine photosynthetic organisms

Environmental (ecological) genomics aims to understand the genetic basis of relationships between organisms and their abiotic and biotic environments. It is a rapidly progressing field of research largely due to recent advances in the speed and volume of genomic data being produced by next generation sequencing (NGS) technologies. Building on information generated by NGS-based approaches, functional genomic methodologies are being applied to identify and characterize genes and gene systems of both environmental and evolutionary relevance. Marine photosynthetic organisms (MPOs) were poorly represented amongst the early genomic models, but this situation is changing rapidly. Here we provide an overview of the recent advances in the application of ecological genomic approaches to both prokaryotic and eukaryotic MPOs. We describe how these approaches are being used to explore the biology and ecology of marine cyanobacteria and algae, particularly with regard to their functions in a broad range of marine ecosystems. Specifically, we review the ecological and evolutionary insights gained from whole genome and transcriptome sequencing projects applied to MPOs and illustrate how their genomes are yielding information on the specific features of these organisms.

Gene copy number variation and its significance in cyanobacterial phylogeny

Background

In eukaryotes, variation in gene copy numbers is often associated with deleterious effects, but may also have positive effects. For prokaryotes, studies on gene copy number variation are rare. Previous studies have suggested that high numbers of rRNA gene copies can be advantageous in environments with changing resource availability, but further association of gene copies and phenotypic traits are not documented. We used one of the morphologically most diverse prokaryotic phyla to test whether numbers of gene copies are associated with levels of cell differentiation.

Results

We implemented a search algorithm that identified 44 genes with highly conserved copies across 22 fully sequenced cyanobacterial taxa. For two very basal cyanobacterial species, Gloeobacter violaceus and a thermophilic Synechococcus species, distinct phylogenetic positions previously found were supported by identical protein coding gene copy numbers. Furthermore, we found that increased ribosomal gene copy numbers showed a strong correlation to cyanobacteria capable of terminal cell differentiation. Additionally, we detected extremely low variation of 16S rRNA sequence copies within the cyanobacteria. We compared our results for 16S rRNA to three other eubacterial phyla (Chroroflexi, Spirochaetes and Bacteroidetes). Based on Bayesian phylogenetic inference and the comparisons of genetic distances, we could confirm that cyanobacterial 16S rRNA paralogs and orthologs show significantly stronger conservation than found in other eubacterial phyla.

Conclusions

A higher number of ribosomal operons could potentially provide an advantage to terminally differentiated cyanobacteria. Furthermore, we suggest that 16S rRNA gene copies in cyanobacteria are homogenized by both concerted evolution and purifying selection. In addition, the small ribosomal subunit in cyanobacteria appears to evolve at extraordinary slow evolutionary rates, an observation that has been made previously for morphological characteristics of cyanobacteria.

The minimal genome: a metabolic and environmental comparison

The field of Synthetic Biology seeks to apply engineering principles to biology in order to produce novel biological systems. One approach to accomplish this goal is the genome-driven cell engineering approach, which searches for functioning minimal genomes in naturally occurring microorganisms, which can then be used as a template for future systems. Currently a prototypical minimal genome has not been discovered. This review analyzes the organisms Mycoplasma pneumoniae, Pelagibacter ubique, Vesicomyosocius okutanii and Prochlorococcus marinus as models of heterotrophic symbiont, heterotrophic free-living, autotrophic symbiont and autotrophic free-living organisms respectively and compares them to the current minimal cell model in order to determine which most closely resembles a true minimal genome. M. pneumoniae possesses a genome of 816 394 base pairs (bp) with 688 open reading frames (ORF) and a severely limited metabolism. Pelagibacter ubique possesses a 1 308 000 bp genome with 1354 ORF and has a fully functional metabolism but requires a reduced form of sulphur. Vesicomyosocius okutanii possesses a 1 020 000 bp genome with 975 ORF and is deficient in the production of threonine, isoleucine and ubiquinone. Prochlorococcus marinus possesses a 1 751 080 bp genome with 1884 ORF and has a complete metabolism with no deficiencies. The current minimal cell model requires a genome to be of limited size, culturalble and having minimal media requirements as such it is the conclusion of this review that P. marinus best fits this model. Further, future research should concentrate on genome reduction experiments using P. marinus and the search for additional minimal genomes should concentrate on autotrophic free-living organisms.

Microbial lifestyle and genome signatures

Microbes are known for their unique ability to adapt to varying lifestyle and environment, even to the extreme or adverse ones. The genomic architecture of a microbe may bear the signatures not only of its phylogenetic position, but also of the kind of lifestyle to which it is adapted. The present review aims to provide an account of the specific genome signatures observed in microbes acclimatized to distinct lifestyles or ecological niches. Niche-specific signatures identified at different levels of microbial genome organization like base composition, GC-skew, purine-pyrimidine ratio, dinucleotide abundance, codon bias, oligonucleotide composition etc. have been discussed. Among the specific cases highlighted in the review are the phenomena of genome shrinkage in obligatory host-restricted microbes, genome expansion in strictly intra-amoebal pathogens, strand-specific codon usage in intracellular species, acquisition of genome islands in pathogenic or symbiotic organisms, discriminatory genomic traits of marine microbes with distinct trophic strategies, and conspicuous sequence features of certain extremophiles like those adapted to high temperature or high salinity.

Hypermutation and stress adaptation in bacteria

Hypermutability is a phenotype characterized by a moderate to high elevation of spontaneous mutation rates and could result from DNA replication errors, defects in error correction mechanisms and many other causes. The elevated mutation rates are helpful to organisms to adapt to sudden and unforeseen threats to survival. At the same time hypermutability also leads to the generation of many deleterious mutations which offset its adaptive value and therefore disadvantageous. Nevertheless, it is very common in nature, especially among clinical isolates of pathogens. Hypermutability is inherited by indirect (second order) selection along with the beneficial mutations generated. At large population sizes and high mutation rates many cells in the population could concurrently acquire beneficial mutations of varying adaptive (fitness) values. These lineages compete with the ancestral cells and also among themselves for fixation. The one with the 'fittest' mutation gets fixed ultimately while the others are lost. This has been called 'clonal interference' which puts a speed limit on adaptation. The original clonal interference hypothesis has been modified recently. Nonheritable (transient) hypermtability conferring significant adaptive benefits also occur during stress response although its molecular basis remains controversial. The adaptive benefits of heritable hypermutability are discussed with emphasis on host-pathogen interactions.

The spontaneous mutation frequencies of Prochlorococcus strains are commensurate with those of other bacteria

The marine cyanobacterium Prochlorococcus, the smallest and most abundant oxygenic phototroph, has an extremely streamlined genome and a high rate of protein evolution. High-light adapted strains of Prochlorococcus in particular have seemingly inadequate DNA repair systems, raising the possibility that inadequate repair may lead to high mutation rates. Prochlorococcus mutation rates have been difficult to determine, in part because traditional methods involving quantifying colonies on solid selective media are not straightforward for this organism. Here we used a liquid dilution method to measure the approximate number of antibiotic-resistant mutants in liquid cultures of Prochlorococcus strains previously unexposed to antibiotic selection. Several antibiotics for which resistance in other bacteria is known to result from a single base pair change were used. The resulting frequencies of antibiotic resistance in Prochlorococcus cultures allowed us to then estimate maximum spontaneous mutation rates, which were similar to those in organisms such as E. coli (∼5.4 × 10(-7) per gene per generation). Therefore, despite the lack of some DNA repair genes, it appears unlikely that the Prochlorcoccus genomes studied here are currently being shaped by unusually high mutation rates.

Ecological and temporal constraints in the evolution of bacterial genomes

Studies on the experimental evolution of microorganisms, on their in vivo evolution (mainly in the case of bacteria producing chronic infections), as well as the availability of multiple full genomic sequences, are placing bacteria in the playground of evolutionary studies. In the present article we review the differential contribution to the evolution of bacterial genomes that processes such as gene modification, gene acquisition and gene loss may have when bacteria colonize different habitats that present characteristic ecological features. In particular, we review how the different processes contribute to evolution in microbial communities, in free-living bacteria or in bacteria living in isolation. In addition, we discuss the temporal constraints in the evolution of bacterial genomes, considering bacterial evolution from the perspective of processes of short-sighted evolution and punctual acquisition of evolutionary novelties followed by long stasis periods.

Genome reduction by deletion of paralogs in the marine cyanobacterium Prochlorococcus

Several isolates of the marine cyanobacterial genus Prochlorococcus have smaller genome sizes than those of the closely related genus Synechococcus. In order to test whether loss of protein-coding genes has contributed to genome size reduction in Prochlorococcus, we reconstructed events of gene family evolution over a strongly supported phylogeny of 12 Prochlorococcus genomes and 9 Synechococcus genomes. Significantly, more events both of loss of paralogs within gene families and of loss of entire gene families occurred in Prochlorococcus than in Synechococcus. The number of nonancestral gene families in genomes of both genera was positively correlated with the extent of genomic islands (GIs), consistent with the hypothesis that horizontal gene transfer (HGT) is associated with GIs. However, even when only isolates with comparable extents of GIs were compared, significantly more events of gene family loss and of paralog loss were seen in Prochlorococcus than in Synechococcus, implying that HGT is not the primary reason for the genome size difference between the two genera.

The UlaG protein family defines novel structural and functional motifs grafted on an ancient RNase fold

Background

Bacterial populations are highly successful at colonizing new habitats and adapting to changing environmental conditions, partly due to their capacity to evolve novel virulence and metabolic pathways in response to stress conditions and to shuffle them by horizontal gene transfer (HGT). A common theme in the evolution of new functions consists of gene duplication followed by functional divergence. UlaG, a unique manganese-dependent metallo-β-lactamase (MBL) enzyme involved in L-ascorbate metabolism by commensal and symbiotic enterobacteria, provides a model for the study of the emergence of new catalytic activities from the modification of an ancient fold. Furthermore, UlaG is the founding member of the so-called UlaG-like (UlaGL) protein family, a recently established and poorly characterized family comprising divalent (and perhaps trivalent) metal-binding MBLs that catalyze transformations on phosphorylated sugars and nucleotides.

Results

Here we combined protein structure-guided and sequence-only molecular phylogenetic analyses to dissect the molecular evolution of UlaG and to study its phylogenomic distribution, its relatedness with present-day UlaGL protein sequences and functional conservation. Phylogenetic analyses indicate that UlaGL sequences are present in Bacteria and Archaea, with bona fide orthologs found mainly in mammalian and plant-associated Gram-negative and Gram-positive bacteria. The incongruence between the UlaGL tree and known species trees indicates exchange by HGT and suggests that the UlaGL-encoding genes provided a growth advantage under changing conditions. Our search for more distantly related protein sequences aided by structural homology has uncovered that UlaGL sequences have a common evolutionary origin with present-day RNA processing and metabolizing MBL enzymes widespread in Bacteria, Archaea, and Eukarya. This observation suggests an ancient origin for the UlaGL family within the broader trunk of the MBL superfamily by duplication, neofunctionalization and fixation.

Conclusions

Our results suggest that the forerunner of UlaG was present as an RNA metabolizing enzyme in the last common ancestor, and that the modern descendants of that ancestral gene have a wide phylogenetic distribution and functional roles. We propose that the UlaGL family evolved new metabolic roles among bacterial and possibly archeal phyla in the setting of a close association with metazoans, such as in the mammalian gastrointestinal tract or in animal and plant pathogens, as well as in environmental settings. Accordingly, the major evolutionary forces shaping the UlaGL family include vertical inheritance and lineage-specific duplication and acquisition of novel metabolic functions, followed by HGT and numerous lineage-specific gene loss events.

Loss of genetic redundancy in reductive genome evolution

Biological systems evolved to be functionally robust in uncertain environments, but also highly adaptable. Such robustness is partly achieved by genetic redundancy, where the failure of a specific component through mutation or environmental challenge can be compensated by duplicate components capable of performing, to a limited extent, the same function. Highly variable environments require very robust systems. Conversely, predictable environments should not place a high selective value on robustness. Here we test this hypothesis by investigating the evolutionary dynamics of genetic redundancy in extremely reduced genomes, found mostly in intracellular parasites and endosymbionts. By combining data analysis with simulations of genome evolution we show that in the extensive gene loss suffered by reduced genomes there is a selective drive to keep the diversity of protein families while sacrificing paralogy. We show that this is not a by-product of the known drivers of genome reduction and that there is very limited convergence to a common core of families, indicating that the repertoire of protein families in reduced genomes is the result of historical contingency and niche-specific adaptations. We propose that our observations reflect a loss of genetic redundancy due to a decreased selection for robustness in a predictable environment.

Dependence of the cyanobacterium Prochlorococcus on hydrogen peroxide scavenging microbes for growth at the ocean's surface

The phytoplankton community in the oligotrophic open ocean is numerically dominated by the cyanobacterium Prochlorococcus, accounting for approximately half of all photosynthesis. In the illuminated euphotic zone where Prochlorococcus grows, reactive oxygen species are continuously generated via photochemical reactions with dissolved organic matter. However, Prochlorococcus genomes lack catalase and additional protective mechanisms common in other aerobes, and this genus is highly susceptible to oxidative damage from hydrogen peroxide (HOOH). In this study we showed that the extant microbial community plays a vital, previously unrecognized role in cross-protecting Prochlorococcus from oxidative damage in the surface mixed layer of the oligotrophic ocean. Microbes are the primary HOOH sink in marine systems, and in the absence of the microbial community, surface waters in the Atlantic and Pacific Ocean accumulated HOOH to concentrations that were lethal for Prochlorococcus cultures. In laboratory experiments with the marine heterotroph Alteromonas sp., serving as a proxy for the natural community of HOOH-degrading microbes, bacterial depletion of HOOH from the extracellular milieu prevented oxidative damage to the cell envelope and photosystems of co-cultured Prochlorococcus, and facilitated the growth of Prochlorococcus at ecologically-relevant cell concentrations. Curiously, the more recently evolved lineages of Prochlorococcus that exploit the surface mixed layer niche were also the most sensitive to HOOH. The genomic streamlining of these evolved lineages during adaptation to the high-light exposed upper euphotic zone thus appears to be coincident with an acquired dependency on the extant HOOH-consuming community. These results underscore the importance of (indirect) biotic interactions in establishing niche boundaries, and highlight the impacts that community-level responses to stress may have in the ecological and evolutionary outcomes for co-existing species.

Differential gene retention in plastids of common recent origin

The cyanobacterium-derived plastids of algae and plants have supported the diversification of much of extant eukaryotic life. Inferences about early events in plastid evolution must rely on reconstructing events that occurred over a billion years ago. In contrast, the photosynthetic amoeba Paulinella chromatophora provides an exceptional model to study organelle evolution in a prokaryote-eukaryote (primary) endosymbiosis that occurred approximately 60 mya. Here we sequenced the plastid genome (0.977 Mb) from the recently described Paulinella FK01 and compared the sequence with the existing data from the sister taxon Paulinella M0880/a. Alignment of the two plastid genomes shows significant conservation of gene order and only a handful of minor gene rearrangements. Analysis of gene content reveals 66 differential gene losses that appear to be outright gene deletions rather than endosymbiotic gene transfers to the host nuclear genome. Phylogenomic analysis validates the plastid ancestor as a member of the Synechococcus-Prochlorococcus group, and the cyanobacterial provenance of all plastid genes suggests that these organelles were not targets of interphylum gene transfers after endosymbiosis. Inspection of 681 DNA alignments of protein-encoding genes shows that the vast majority have dN/dS ratios <<1, providing evidence for purifying selection. Our study demonstrates that plastid genomes in sister taxa are strongly constrained by selection but follow distinct trajectories during the earlier phases of organelle evolution.

Distinct, ecotype-specific genome and proteome signatures in the marine cyanobacteria Prochlorococcus

BACKGROUND

The marine cyanobacterium Prochlorococcus marinus, having multiple ecotypes of distinct genotypic/phenotypic traits and being the first documented example of genome shrinkage in free-living organisms, offers an ideal system for studying niche-driven molecular micro-diversity in closely related microbes. The present study, through an extensive comparative analysis of various genomic/proteomic features of 6 high light (HL) and 6 low light (LL) adapted strains, makes an attempt to identify molecular determinants associated with their vertical niche partitioning.

RESULTS

Pronounced strand-specific asymmetry in synonymous codon usage is observed exclusively in LL strains. Distinct dinucleotide abundance profiles are exhibited by 2 LL strains with larger genomes and G+C-content approximately 50% (group LLa), 4 LL strains having reduced genomes and G+C-content approximately 35-37% (group LLb), and 6 HL strains. Taking into account the emergence of LLa, LLb and HL strains (based on 16S rRNA phylogeny), a gradual increase in average aromaticity, pI values and beta- & coil-forming propensities and a decrease in mean hydrophobicity, instability indices and helix-forming propensities of core proteins are observed. Greater variations in orthologous gene repertoire are found between LLa and LLb strains, while higher number of positively selected genes exist between LL and HL strains.

CONCLUSION

Strains of different Prochlorococcus groups are characterized by distinct compositional, physicochemical and structural traits that are not mere remnants of a continuous genetic drift, but are potential outcomes of a grand scheme of niche-oriented stepwise diversification, that might have driven them chronologically towards greater stability/fidelity and invoked upon them a special ability to inhabit diverse oceanic environments.

Prochlorococcus: advantages and limits of minimalism

Prochlorococcus is the key phytoplanktonic organism of tropical gyres, large ocean regions that are depleted of the essential macronutrients needed for photosynthesis and cell growth. This cyanobacterium has adapted itself to oligotrophy by minimizing the resources necessary for life through a drastic reduction of cell and genome sizes. This rarely observed strategy in free-living organisms has conferred on Prochlorococcus a considerable advantage over other phototrophs, including its closest relative Synechococcus, for life in this vast yet little variable ecosystem. However, this strategy seems to reach its limits in the upper layer of the S Pacific gyre, the most oligotrophic region of the world ocean. By losing some important genes and/or functions during evolution, Prochlorococcus has seemingly become dependent on co-occurring microorganisms. In this review, we present some of the recent advances in the ecology, biology, and evolution of Prochlorococcus, which because of its ecological importance and tiny genome is rapidly imposing itself as a model organism in environmental microbiology.

The life and death of gene families

One of the unique insights provided by the growing number of fully sequenced genomes is the pervasiveness of gene duplication and gene loss. Indeed, several metrics now suggest that rates of gene birth and death per gene are only 10-40% lower than nucleotide substitutions per site, and that per nucleotide, the consequent lineage-specific expansion and contraction of gene families may play at least as large a role in adaptation as changes in orthologous sequences. While gene family evolution is pervasive, it may be especially important in our own evolution since it appears that the "revolving door" of gene duplication and loss has undergone multiple accelerations in the lineage leading to humans. In this paper, we review current understanding of gene family evolution including: methods for inferring copy number change, evidence for adaptive expansion and adaptive contraction of gene families, the origins of new families and deaths of previously established ones, and finally we conclude with a perspective on challenges and promising directions for future research.