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

Comparative genomics of a vertically transmitted thiotrophic bacterial ectosymbiont and its close free-living relative

Thiotrophic symbioses between sulphur-oxidizing bacteria and various unicellular and metazoan eukaryotes are widespread in reducing marine environments. The giant colonial ciliate Zoothamnium niveum, however, is the only host of thioautotrophic symbionts that has been cultivated along with its symbiont, the vertically transmitted ectosymbiont Candidatus Thiobius zoothamnicola (short Thiobius). Because theoretical predictions posit a smaller genome in vertically transmitted endosymbionts compared to free-living relatives, we investigated whether this is true also for an ectosymbiont. We used metagenomics to recover the high-quality draft genome of this bacterial symbiont. For comparison we have also sequenced a closely related free-living cultured but not formally described strain Milos ODIII6 (short ODIII6). We then performed comparative genomics to assess the functional capabilities at gene, metabolic pathway and trait level. 16S rRNA gene trees and average amino acid identity confirmed the close phylogenetic relationship of both bacteria. Indeed, Thiobius has about a third smaller genome than its free-living relative ODIII6, with reduced metabolic capabilities and fewer functional traits. The functional capabilities of Thiobius were a subset of those of the more versatile ODIII6, which possessed additional genes for oxygen, sulphur and hydrogen utilization and for the acquisition of phosphorus illustrating features that may be adaptive for the unstable environmental conditions at hydrothermal vents. In contrast, Thiobius possesses genes potentially enabling it to utilize lactate and acetate heterotrophically, compounds that may be provided as byproducts by the host. The present study illustrates the effect of strict host-dependence of a bacterial ectosymbiont on genome evolution and host adaptation.

Functional divergence in the proteins encoded by ARO80 from S. uvarum, S. kudriavzevii and S. cerevisiae explain differences in the aroma production during wine fermentation

Phenylethanol (PE) and phenylethyl acetate (PEA) are commonly desired compounds in wine because of their rose-like aroma. The yeast S. cerevisiae produces the PE either through de novo biosynthesis by shikimate pathway followed by the Ehrlich pathway or the direct phenylalanine catabolism via Ehrlich pathway, and then converted into PEA. Previous work demonstrated that, compared to S. cerevisiae, other Saccharomyces species, such as S. kudriavzevii and S. uvarum, produce higher concentrations of PE and PEA from the precursor phenylalanine, which indicates differential activities of the biosynthetic-involved enzymes. A previous in-silico analysis suggested that the transcriptional activator Aro80p is one of the best candidates to explain these differences. An improved functional analysis identified significant radical amino acid changes in the S. uvarum and S. kudriavzevii Aro80p that could impact the expression of the catabolic genes ARO9 and ARO10, and hence, the production of PE from phenylalanine. Indeed, wine S. cerevisiae strains carrying the S. uvarum and S. kudriavzevii ARO80 alleles increased the production of both compounds in the presence of phenylalanine by increasing the expression of ARO9 and ARO10. This study provides novel insights of the unidentified Aro80p regulatory region and the potential usage of alternatives ARO80 alleles to enhance the PE and PEA concentration in wine.

Bacterial Group II Intron Genomic Neighborhoods Reflect Survival Strategies: Hiding and Hijacking

Group II (gII) introns are mobile retroelements that can spread to new DNA sites through retrotransposition, which can be influenced by a variety of host factors. To determine if these host factors bear any relationship to the genomic location of gII introns, we developed a bioinformatic pipeline wherein we focused on the genomic neighborhoods of bacterial gII introns within their native contexts and sought to determine global relationships between introns and their surrounding genes. We found that, although gII introns inhabit diverse regions, these neighborhoods are often functionally enriched for genes that could promote gII intron retention or proliferation. On one hand, we observe that gII introns are frequently found hiding in mobile elements or after transcription terminators. On the other hand, gII introns are enriched in locations in which they could hijack host functions for their movement, potentially timing expression of the intron with genes that produce favorable conditions for retrotransposition. Thus, we propose that gII intron distributions have been shaped by relationships with their surrounding genomic neighbors.

Methylation of rRNA as a host defense against rampant group II intron retrotransposition

BACKGROUND

Group II introns are mobile retroelements, capable of invading new sites in DNA. They are self-splicing ribozymes that complex with an intron-encoded protein to form a ribonucleoprotein that targets DNA after splicing. These molecules can invade DNA site-specifically, through a process known as retrohoming, or can invade ectopic sites through retrotransposition. Retrotransposition, in particular, can be strongly influenced by both environmental and cellular factors.

RESULTS

To investigate host factors that influence retrotransposition, we performed random insertional mutagenesis using the ISS1 transposon to generate a library of over 1000 mutants in Lactococcus lactis, the native host of the Ll.LtrB group II intron. By screening this library, we identified 92 mutants with increased retrotransposition frequencies (RTP-ups). We found that mutations in amino acid transport and metabolism tended to have increased retrotransposition frequencies. We further explored a subset of these RTP-up mutants, the most striking of which is a mutant in the ribosomal RNA methyltransferase rlmH, which exhibited a reproducible 20-fold increase in retrotransposition frequency. In vitro and in vivo experiments revealed that ribosomes in the rlmH mutant were defective in the m3Ψ modification and exhibited reduced binding to the intron RNA.

CONCLUSIONS

Taken together, our results reinforce the importance of the native host organism in regulating group II intron retrotransposition. In particular, the evidence from the rlmH mutant suggests a role for ribosome modification in limiting rampant retrotransposition.

Comparative Genomics Between Saccharomyces kudriavzevii and S. cerevisiae Applied to Identify Mechanisms Involved in Adaptation

Yeasts belonging to the Saccharomyces genus play an important role in human-driven fermentations. The species S. cerevisiae has been widely studied because it is the dominant yeast in most fermentations and it has been widely used as a model eukaryotic organism. Recently, other species of the Saccharomyces genus are gaining interest to solve the new challenges that the fermentation industry are facing. One of these species is S. kudriavzevii, which exhibits interesting physiological properties compared to S. cerevisiae, such as a better adaptation to grow at low temperatures, a higher glycerol synthesis and lower ethanol production. The aim of this study is to understand the molecular basis behind these phenotypic differences of biotechnological interest by using a species-based comparative genomics approach. In this work, we sequenced, assembled and annotated two new genomes of S. kudriavzevii. We used a combination of different statistical methods to identify functional divergence, signatures of positive selection and acceleration of substitution rates at specific amino acid sites of proteins in S. kudriavzevii when compared to S. cerevisiae, and vice versa. We provide a list of candidate genes in which positive selection could be acting during the evolution of both S. cerevisiae and S. kudriavzevii clades. Some of them could be related to certain important differences in metabolism previously reported by other authors such us DAL3 and ARO4, involved in nitrogen assimilation and amino acid biosynthesis. In addition, three of those genes (FBA1, ZIP1, and RQC2) showed accelerated evolutionary rates in Sk branch. Finally, genes of the riboflavin biosynthesis were also among those genes with a significant higher rate of nucleotide substitution and those proteins have amino acid positions contributing to functional divergence.

Proteome-wide analysis of functional divergence in bacteria: exploring a host of ecological adaptations

Functional divergence is the process by which new genes and functions originate through the modification of existing ones. Both genetic and environmental factors influence the evolution of new functions, including gene duplication or changes in the ecological requirements of an organism. Novel functions emerge at the expense of ancestral ones and are generally accompanied by changes in the selective forces at constrained protein regions. We present software capable of analyzing whole proteomes, identifying putative amino acid replacements leading to functional change in each protein and performing statistical tests on all tabulated data. We apply this method to 750 complete bacterial proteomes to identify high-level patterns of functional divergence and link these patterns to ecological adaptations. Proteome-wide analyses of functional divergence in bacteria with different ecologies reveal a separation between proteins involved in information processing (Ribosome biogenesis etc.) and those which are dependent on the environment (energy metabolism, defense etc.). We show that the evolution of pathogenic and symbiotic bacteria is constrained by their association with the host, and also identify unusual events of functional divergence even in well-studied bacteria such as Escherichia coli. We present a description of the roles of phylogeny and ecology in functional divergence at the level of entire proteomes in bacteria.

Purifying selection and molecular adaptation in the genome of Verminephrobacter, the heritable symbiotic bacteria of earthworms

While genomic erosion is common among intracellular symbionts, patterns of genome evolution in heritable extracellular endosymbionts remain elusive. We study vertically transmitted extracellular endosymbionts (Verminephrobacter, Betaproteobacteria) that form a beneficial, species-specific, and evolutionarily old (60–130 Myr) association with earthworms. We assembled a draft genome of Verminephrobacter aporrectodeae and compared it with the genomes of Verminephrobacter eiseniae and two nonsymbiotic close relatives (Acidovorax). Similar to V. eiseniae, the V. aporrectodeae genome was not markedly reduced in size and showed no A–T bias. We characterized the strength of purifying selection (ω = dN/dS) and codon usage bias in 876 orthologous genes. Symbiont genomes exhibited strong purifying selection (ω = 0.09 ± 0.07), although transition to symbiosis entailed relaxation of purifying selection as evidenced by 50% higher ω values and less codon usage bias in symbiont compared with reference genomes. Relaxation was not evenly distributed among functional gene categories but was overrepresented in genes involved in signal transduction and cell envelope biogenesis. The same gene categories also harbored instances of positive selection in the Verminephrobacter clade. In total, positive selection was detected in 89 genes, including also genes involved in DNA metabolism, tRNA modification, and TonB-dependent iron uptake, potentially highlighting functions important in symbiosis. Our results suggest that the transition to symbiosis was accompanied by molecular adaptation, while purifying selection was only moderately relaxed, despite the evolutionary age and stability of the host association. We hypothesize that biparental transmission of symbionts and rare genetic mixing during transmission can prevent genome erosion in heritable symbionts.

Evolutionary microbial genomics: insights into bacterial host adaptation

Host-adapted bacteria include mutualists and pathogens of animals, plants and insects. Their study is therefore important for biotechnology, biodiversity and human health. The recent rapid expansion in bacterial genome data has provided insights into the adaptive, diversifying and reductive evolutionary processes that occur during host adaptation. The results have challenged many pre-existing concepts built from studies of laboratory bacterial strains. Furthermore, recent studies have revealed genetic changes associated with transitions from parasitism to mutualism and opened new research avenues to understand the functional reshaping of bacteria as they adapt to growth in the cytoplasm of a eukaryotic host.

New paradigms for the evolution of beneficial infections

A longstanding paradigm predicts that microbial parasites and mutualists exhibit disparate evolutionary patterns. Parasites are predicted to promote arms races with hosts, rapid evolution and sexual recombination. By contrast, mutualists have been linked with beneficial coadaptation, evolutionary stasis and asexuality. In this review we discuss the recent surge of molecular data on microbes that are being used to test and reshape these ideas. New analyses reveal that beneficial microbes often share mechanisms of infection and defense with parasites, and can also exhibit rapid evolution and extensive genetic exchange. To explain these patterns, new paradigms must take into account the varied population biology of beneficial microbes, their potential conflicts with hosts, and the mosaic nature of genome evolution that requires locus-based tests to analyze the genetics of host adaptation.

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.

Structural calibration of the rates of amino acid evolution in a search for Darwin in drifting biological systems

In the last two decades, many reports of proteins under positive selection have brought the neutral theory into question. However, the methods used to detect selection have ignored the evolvability of amino acids within proteins, which is fundamental to distinguishing positive selection from the relaxed constraints caused by genetic drift. Disentangling these two counterbalancing forces is essential to test the neutral theory. Here, we calibrate rates of amino acid divergence by using structural information from the full set of crystallized proteins in bacteria. In agreement with previous reports, we show that rates of amino acid evolution correlate negatively with the number of per-amino acid atomic interactions. Calibration of the rates of evolution allows identifying signatures of selection in biological systems that evolve under strong genetic drift, such as endosymbiotic bacteria. Application of this method identifies different rates and dynamics of evolution for highly connected amino acids in the structure compared with sparsely connected ones. We also unearth patterns of Darwinian selection in fundamental cellular proteins in endosymbiotic bacteria including the cochaperonin GroES, ribosomal proteins, proteins involved in cell cycle control, DNA-binding proteins, and proteins involved in DNA replication and repair. This is, to our knowledge, the first attempt to distinguish adaptive evolution from relaxed constraints in biological systems under genetic drift.

Two chaperonin systems in bacterial genomes with distinct ecological roles

Bacterial chaperonins are essential to cell viability and have a role in endosymbiosis, which leads to increased biological complexity. However, the extent to which chaperonins promote ecological innovation is unknown. We screened 622 bacterial genomes for genes encoding chaperonins, and found archaeal-like chaperonins in bacteria that inhabit archaeal ecological niches. We found that chaperonins encoded in pathogenic bacteria are the most functionally divergent. We identified the molecular basis of the dramatic structural changes in mitochondrial GROEL, a highly derived chaperonin gene. Our analysis suggests that chaperonins are important capacitors of evolutionary and ecological change.