Rapid evolution of the sequences and gene repertoires of secreted proteins in bacteria

Freshwater genome-reduced bacteria exhibit pervasive episodes of adaptive stasis

The emergence of bacterial species is rooted in their inherent potential for continuous evolution and adaptation to an ever-changing ecological landscape. The adaptive capacity of most species frequently resides within the repertoire of genes encoding the secreted proteome (SP), as it serves as a primary interface used to regulate survival/reproduction strategies. Here, by applying evolutionary genomics approaches to metagenomics data, we show that abundant freshwater bacteria exhibit biphasic adaptation states linked to the eco-evolutionary processes governing their genome sizes. While species with average to large genomes adhere to the dominant paradigm of evolution through niche adaptation by reducing the evolutionary pressure on their SPs (via the augmentation of functionally redundant genes that buffer mutational fitness loss) and increasing the phylogenetic distance of recombination events, most of the genome-reduced species exhibit a nonconforming state. In contrast, their SPs reflect a combination of low functional redundancy and high selection pressure, resulting in significantly higher levels of conservation and invariance. Our findings indicate that although niche adaptation is the principal mechanism driving speciation, freshwater genome-reduced bacteria often experience extended periods of adaptive stasis. Understanding the adaptive state of microbial species will lead to a better comprehension of their spatiotemporal dynamics, biogeography, and resilience to global change.

SOCfinder: a genomic tool for identifying social genes in bacteria

Bacteria cooperate by working collaboratively to defend their colonies, share nutrients, and resist antibiotics. Nevertheless, our understanding of these remarkable behaviours primarily comes from studying a few well-characterized species. Consequently, there is a significant gap in our understanding of microbial social traits, particularly in natural environments. To address this gap, we can use bioinformatic tools to identify genes that control cooperative or otherwise social traits. Existing tools address this challenge through two approaches. One approach is to identify genes that encode extracellular proteins, which can provide benefits to neighbouring cells. An alternative approach is to predict gene function using annotation tools. However, these tools have several limitations. Not all extracellular proteins are cooperative, and not all cooperative behaviours are controlled by extracellular proteins. Furthermore, existing functional annotation methods frequently miss known cooperative genes. We introduce SOCfinder as a new tool to find bacterial genes that control cooperative or otherwise social traits. SOCfinder combines information from several methods, considering if a gene is likely to [1] code for an extracellular protein [2], have a cooperative functional annotation, or [3] be part of the biosynthesis of a cooperative secondary metabolite. We use data on two extensively-studied species (P. aeruginosa and B. subtilis) to show that SOCfinder is better at finding known cooperative genes than existing tools. We also use theory from population genetics to identify a signature of kin selection in SOCfinder cooperative genes, which is lacking in genes identified by existing tools. SOCfinder opens up a number of exciting directions for future research, and is available to download from https://github.com/lauriebelch/SOCfinder.

Gut immune responses and evolution of the gut microbiome-a hypothesis

The gut microbiome is an assemblage of microbes that have profound effects on their hosts. The composition of the microbiome is affected by bottom-up, among-taxa interactions and by top-down, host effects, which includes the host immune response. While the high-level composition of the microbiome is generally stable over time, component strains and genotypes will constantly be evolving, with both bottom-up and top-down effects acting as selection pressures, driving microbial evolution. Secretory IgA is a major feature of the gut's adaptive immune response, and a substantial proportion of gut bacteria are coated with IgA, though the effect of this on bacteria is unclear. Here we hypothesize that IgA binding to gut bacteria is a selection pressure that will drive the evolution of IgA-bound bacteria, so that they will have a different evolutionary trajectory than those bacteria not bound by IgA. We know very little about the microbiome of wild animals and even less about their gut immune responses, but it must be a priority to investigate this hypothesis to understand if and how host immune responses contribute to microbiome evolution.

Signatures of kin selection in a natural population of the bacteria Bacillus subtilis

Laboratory experiments have suggested that bacteria perform a range of cooperative behaviors, which are favored because they are directed toward relatives (kin selection). However, there is a lack of evidence for cooperation and kin selection in natural bacterial populations. Molecular population genetics offers a promising method to study natural populations because the theory predicts that kin selection will lead to relaxed selection, which will result in increased polymorphism and divergence at cooperative genes. Examining a natural population of Bacillus subtilis, we found consistent evidence that putatively cooperative traits have higher polymorphism and greater divergence than putatively private traits expressed at the same rate. In addition, we were able to eliminate alternative explanations for these patterns and found more deleterious mutations in genes controlling putatively cooperative traits. Overall, our results suggest that cooperation is favored by kin selection, with an average relatedness of r = .79 between interacting individuals.

Is cooperation favored by horizontal gene transfer?

It has been hypothesized that horizontal gene transfer on plasmids can facilitate the evolution of cooperation, by allowing genes to jump between bacteria, and hence increase genetic relatedness at the cooperative loci. However, we show theoretically that horizontal gene transfer only appreciably increases relatedness when plasmids are rare, where there are many plasmid-free cells available to infect (many opportunities for horizontal gene transfer). In contrast, when plasmids are common, there are few opportunities for horizontal gene transfer, meaning relatedness is not appreciably increased, and so cooperation is not favored. Plasmids, therefore, evolve to be rare and cooperative, or common and noncooperative, meaning plasmid frequency and cooperativeness are never simultaneously high. The overall level of plasmid-mediated cooperation, given by the product of plasmid frequency and cooperativeness, is therefore consistently negligible or low.

Characterization and Spatial Mapping of the Human Gut Metasecretome

Bacterially secreted proteins play an important role in microbial physiology and ecology in many environments, including the mammalian gut. While gut microbes have been extensively studied over the past decades, little is known about the proteins that they secrete into the gastrointestinal tract. In this study, we developed and applied a computational pipeline to a comprehensive catalog of human-associated metagenome-assembled genomes in order to predict and analyze the bacterial metasecretome of the human gut, i.e., the collection of proteins secreted out of the cytoplasm by human gut bacteria. We identified the presence of large and diverse families of secreted carbohydrate-active enzymes and assessed their phylogenetic distributions across different taxonomic groups, which revealed an enrichment in Bacteroidetes and Verrucomicrobia. By mapping secreted proteins to available metagenomic data from endoscopic sampling of the human gastrointestinal tract, we specifically pinpointed regions in the upper and lower intestinal tract along the lumen and mucosa where specific glycosidases are secreted by resident microbes. The metasecretome analyzed in this study constitutes the most comprehensive list of secreted proteins produced by human gut bacteria reported to date and serves as a useful resource for the microbiome research community. IMPORTANCE Bacterially secreted proteins are necessary for the proper functioning of bacterial cells and communities. Secreted proteins provide bacterial cells with the ability to harvest resources from the exterior, import these resources into the cell, and signal to other bacteria. In the human gut microbiome, these actions impact host health and allow the maintenance of a healthy gut bacterial community. We utilized computational tools to identify the major components of human gut bacterially secreted proteins and determined their spatial distribution in the gastrointestinal tract. Our analysis of human gut bacterial secreted proteins will allow a better understanding of the impact of gut bacteria on human health and represents a step toward identifying new protein functions with interesting applications in biomedicine and industry.

Gene transferability and sociality do not correlate with gene connectivity

The connectivity of a gene, defined as the number of interactions a gene's product has with other genes' products, is a key characteristic of a gene. In prokaryotes, the complexity hypothesis predicts that genes which undergo more frequent horizontal transfer will be less connected than genes which are only very rarely transferred. We tested the role of horizontal gene transfer, and other potentially important factors, by examining the connectivity of chromosomal and plasmid genes, across 134 diverse prokaryotic species. We found that (i) genes on plasmids were less connected than genes on chromosomes; (ii) connectivity of plasmid genes was not correlated with plasmid mobility; and (iii) the sociality of genes (cooperative or private) was not correlated with gene connectivity.

Evolution-proof inhibitors of public good cooperation: a screening strategy inspired by social evolution theory

Interference with public good cooperation provides a promising novel antimicrobial strategy since social evolution theory predicts that resistant mutants will be counter-selected if they share the public benefits of their resistance with sensitive cells in the population. Although this hypothesis is supported by a limited number of pioneering studies, an extensive body of more fundamental work on social evolution describes a multitude of mechanisms and conditions that can stabilize public behaviour, thus potentially allowing resistant mutants to thrive. In this paper we theorize on how these different mechanisms can influence the evolution of resistance against public good inhibitors. Based hereon, we propose an innovative 5-step screening strategy to identify novel evolution-proof public good inhibitors, which involves a systematic evaluation of the exploitability of public goods under the most relevant experimental conditions, as well as a careful assessment of the most optimal way to interfere with their action. Overall, this opinion paper is aimed to contribute to long-term solutions to fight bacterial infections.

Impact of horizontal gene transfer on emergence and stability of cooperative virulence in Salmonella Typhimurium

Intestinal inflammation fuels the transmission of Salmonella Typhimurium (S.Tm). However, a substantial fitness cost is associated with virulence expression. Mutations inactivating transcriptional virulence regulators generate attenuated variants profiting from inflammation without enduring virulence cost. Such variants interfere with the transmission of fully virulent clones. Horizontal transfer of functional regulatory genes (HGT) into attenuated variants could nevertheless favor virulence evolution. To address this hypothesis, we cloned hilD, coding for the master regulator of virulence, into a conjugative plasmid that is highly transferrable during intestinal colonization. The resulting mobile hilD allele allows virulence to emerge from avirulent populations, and to be restored in attenuated mutants competing against virulent clones within-host. However, mutations inactivating the mobile hilD allele quickly arise. The stability of virulence mediated by HGT is strongly limited by its cost, which depends on the hilD expression level, and by the timing of transmission. We conclude that robust evolution of costly virulence expression requires additional selective forces such as narrow population bottlenecks during transmission.

Salmonella Typhimurium virulence is costly and can be lost by mutation during infection. Bakkeren et al. show that virulence restoration via horizontal gene transfer is only transient while transmission bottlenecks promote long-term virulence stability.

Kin selection for cooperation in natural bacterial populations

Bacteria secrete many molecules outside the cell, where they provide benefits to other cells. One potential reason for producing these “public goods” is that they benefit closely related cells that share the gene for cooperation (kin selection). While many laboratory studies have supported this hypothesis, there is a lack of evidence that kin selection favors cooperation in natural populations. We examined bacterial genomes from the environment and used population genetics theory to analyze the DNA sequences. Our analyses suggest that public goods cooperation has indeed been favored by kin selection in natural populations.

Bacteria produce a range of molecules that are secreted from the cell and can provide a benefit to the local population of cells. Laboratory experiments have suggested that these “public goods” molecules represent a form of cooperation, favored because they benefit closely related cells (kin selection). However, there is a relative lack of data demonstrating kin selection for cooperation in natural populations of bacteria. We used molecular population genetics to test for signatures of kin selection at the genomic level in natural populations of the opportunistic pathogen Pseudomonas aeruginosa. We found consistent evidence from multiple traits that genes controlling putatively cooperative traits have higher polymorphism and greater divergence and are more likely to harbor deleterious mutations relative to genes controlling putatively private traits, which are expressed at similar rates. These patterns suggest that cooperative traits are controlled by kin selection, and we estimate that the relatedness for social interactions in P. aeruginosa is r = 0.84. More generally, our results demonstrate how molecular population genetics can be used to study the evolution of cooperation in natural populations.

The evolution and role of eukaryotic-like domains in environmental intracellular bacteria: the battle with a eukaryotic cell

Intracellular pathogens that are able to thrive in different environments, such as Legionella spp. which preferentially live in protozoa in aquatic environments or environmental Chlamydiae which replicate either within protozoa or a range of animals, possess a plethora of cellular biology tools to influence their eukaryotic host. The host manipulation tools that evolved in the interaction with protozoa, confer these bacteria the capacity to also infect phylogenetically distinct eukaryotic cells, such as macrophages and thus they can also be human pathogens. To manipulate the host cell, bacteria use protein secretion systems and molecular effectors. Although these molecular effectors are encoded in bacteria, they are expressed and function in a eukaryotic context often mimicking or inhibiting eukaryotic proteins. Indeed, many of these effectors have eukaryotic-like domains. In this review we propose that the main pathways environmental intracellular bacteria need to subvert in order to establish the host eukaryotic cell as a replication niche are chromatin remodelling, ubiquitination signalling, and modulation of protein-protein interactions via tandem repeat domains. We then provide mechanistic insight into how these proteins might have evolved as molecular weapons. Finally, we highlight that in environmental intracellular bacteria the number of eukaryotic-like domains and proteins is considerably higher than in intracellular bacteria specialised to an isolated niche, such as obligate intracellular human pathogens. As mimics of eukaryotic proteins are critical components of host pathogen interactions, this distribution of eukaryotic-like domains suggests that the environment has selected them.

Plasmids do not consistently stabilize cooperation across bacteria but may promote broad pathogen host-range

Horizontal gene transfer via plasmids could favour cooperation in bacteria, because transfer of a cooperative gene turns non-cooperative cheats into cooperators. This hypothesis has received support from theoretical, genomic and experimental analyses. In contrast, we show here, with a comparative analysis across 51 diverse species, that genes for extracellular proteins, which are likely to act as cooperative ‘public goods’, were not more likely to be carried on either: (i) plasmids compared to chromosomes; or (ii) plasmids that transfer at higher rates. Our results were supported by theoretical modelling which showed that while horizontal gene transfer can help cooperative genes initially invade a population, it has less influence on the longer-term maintenance of cooperation. Instead, we found that genes for extracellular proteins were more likely to be on plasmids when they coded for pathogenic virulence traits, in pathogenic bacteria with a broad host-range.

Metagenomics and Other Omics Approaches to Bacterial Communities and Antimicrobial Resistance Assessment in Aquacultures

The shortage of wild fishery resources and the rising demand for human nutrition has driven a great expansion in aquaculture during the last decades in terms of production and economic value. As such, sustainable aquaculture production is one of the main priorities of the European Union's 2030 agenda. However, the intensification of seafood farming has resulted in higher risks of disease outbreaks and in the increased use of antimicrobials to control them. The selective pressure exerted by these drugs provides the ideal conditions for the emergence of antimicrobial resistance hotspots in aquaculture facilities. Omics technology is an umbrella term for modern technologies such as genomics, metagenomics, transcriptomics, proteomics, culturomics, and metabolomics. These techniques have received increasing recognition because of their potential to unravel novel mechanisms in biological science. Metagenomics allows the study of genomes in microbial communities contained within a certain environment. The potential uses of metagenomics in aquaculture environments include the study of microbial diversity, microbial functions, and antibiotic resistance genes. A snapshot of these high throughput technologies applied to microbial diversity and antimicrobial resistance studies in aquacultures will be presented in this review.

Strain Structure and Dynamics Revealed by Targeted Deep Sequencing of the Honey Bee Gut Microbiome

The factors driving fine-scale composition and dynamics of gut microbial communities are poorly understood. In this study, we used metagenomic amplicon deep sequencing to decipher the strain dynamics of two key members of the honey bee gut microbiome. Using this high-throughput and cost-effective approach, we were able to confirm results from previous large-scale whole-genome shotgun (WGS) metagenomic sequencing studies while also gaining additional insights into the community dynamics of two core members of the honey bee gut microbiome. Moreover, we were able to show that cryptic strains are not responsible for the observed variations in microbiome composition across bees.

Host-associated microbiomes can be critical for the health and proper development of animals and plants. The answers to many fundamental questions regarding the modes of acquisition and microevolution of microbiome communities remain to be established. Deciphering strain-level dynamics is essential to fully understand how microbial communities evolve, but the forces shaping the strain-level dynamics of microbial communities remain largely unexplored, mostly because of methodological issues and cost. Here, we used targeted strain-level deep sequencing to uncover the strain dynamics within a host-associated microbial community using the honey bee gut microbiome as a model system. Our results revealed that amplicon sequencing of conserved protein-coding gene regions using species-specific primers is a cost-effective and accurate method for exploring strain-level diversity. In fact, using this method we were able to confirm strain-level results that have been obtained from whole-genome shotgun sequencing of the honey bee gut microbiome but with a much higher resolution. Importantly, our deep sequencing approach allowed us to explore the impact of low-frequency strains (i.e., cryptic strains) on microbiome dynamics. Results show that cryptic strain diversity is not responsible for the observed variations in microbiome composition across bees. Altogether, the findings revealed new fundamental insights regarding strain dynamics of host-associated microbiomes.

IMPORTANCE The factors driving fine-scale composition and dynamics of gut microbial communities are poorly understood. In this study, we used metagenomic amplicon deep sequencing to decipher the strain dynamics of two key members of the honey bee gut microbiome. Using this high-throughput and cost-effective approach, we were able to confirm results from previous large-scale whole-genome shotgun (WGS) metagenomic sequencing studies while also gaining additional insights into the community dynamics of two core members of the honey bee gut microbiome. Moreover, we were able to show that cryptic strains are not responsible for the observed variations in microbiome composition across bees.

Phylogenetic profiling, an untapped resource for the prediction of secreted proteins and its complementation with sequence-based classifiers in bacterial type III, IV and VI secretion systems

In the establishment and maintenance of the interaction between pathogenic or symbiotic bacteria with a eukaryotic organism, protein substrates of specialized bacterial secretion systems called effectors play a critical role once translocated into the host cell. Proteins are also secreted to the extracellular medium by free-living bacteria or directly injected into other competing organisms to hinder or kill. In this work, we explore an approach based on the evolutionary dependence that most of the effectors maintain with their specific secretion system that analyzes the co-occurrence of any orthologous protein group and their corresponding secretion system across multiple genomes. We compared and complemented our methodology with sequence-based machine learning prediction tools for the type III, IV and VI secretion systems. Finally, we provide the predictive results for the three secretion systems in 1606 complete genomes at http://www.iib.unsam.edu.ar/orgsissec/.

Community diversity and habitat structure shape the repertoire of extracellular proteins in bacteria

We test the hypothesis that the frequency and cost of extracellular proteins produced by bacteria, which often depend on cooperative processes, vary with habitat structure and community diversity. The integration of the environmental distribution of bacteria (using 16S datasets) and their genomes shows that bacteria living in more structured habitats encode more extracellular proteins. In contrast, the effect of community diversity depends on protein function: it’s positive for proteins implicated in antagonistic interactions and negative for those involved in nutrient acquisition. Extracellular proteins are costly and endure stronger selective pressure for low cost and for low diffusivity in less structured habitats and in more diverse communities. Finally, Bacteria found in multiple types of habitats, including host-associated generalists, encode more extracellular proteins than niche-restricted bacteria. These results show that ecological variables, notably habitat structure and community diversity, shape the evolution of the repertoires of genes encoding extracellular proteins and thus affect the ability of bacteria to manipulate their environment.

Microbes secrete a repertoire of extracellular proteins to serve various functions depending on the ecological context. Here the authors examine how bacterial community composition and habitat structure affect the extracellular proteins, showing that generalist species and those living in more structured environments produce more extracellular proteins, and that costs of production are lower in more diverse communities.

ChlamDB: a comparative genomics database of the phylum Chlamydiae and other members of the Planctomycetes-Verrucomicrobiae-Chlamydiae superphylum

ChlamDB is a comparative genomics database containing 277 genomes covering the entire Chlamydiae phylum as well as their closest relatives belonging to the Planctomycetes-Verrucomicrobiae-Chlamydiae (PVC) superphylum. Genomes can be compared, analyzed and retrieved using accessions numbers of the most widely used databases including COG, KEGG ortholog, KEGG pathway, KEGG module, Pfam and InterPro. Gene annotations from multiple databases including UniProt (curated and automated protein annotations), KEGG (annotation of pathways), COG (orthology), TCDB (transporters), STRING (protein-protein interactions) and InterPro (domains and signatures) can be accessed in a comprehensive overview page. Candidate effectors of the Type III secretion system (T3SS) were identified using four in silico methods. The identification of orthologs among all PVC genomes allows users to perform large-scale comparative analyses and to identify orthologs of any protein in all genomes integrated in the database. Phylogenetic relationships of PVC proteins and their closest homologs in RefSeq, comparison of transmembrane domains and Pfam domains, conservation of gene neighborhood and taxonomic profiles can be visualized using dynamically generated graphs, available for download. As a central resource for researchers working on chlamydia, chlamydia-related bacteria, verrucomicrobia and planctomyces, ChlamDB facilitates the access to comprehensive annotations, integrates multiple tools for comparative genomic analyses and is freely available at https://chlamdb.ch/. Database URL: https://chlamdb.ch/.

Antibiotic Resistance Gene Diversity and Virulence Gene Diversity Are Correlated in Human Gut and Environmental Microbiomes

Human beings have used large amounts of antibiotics, not only in medical contexts but also, for example, as growth factors in agriculture and livestock, resulting in the contamination of the environment. Even when pathogenic bacteria are the targets of antibiotics, hundreds of nonpathogenic bacterial species are affected as well. Therefore, both pathogenic and nonpathogenic bacteria have gradually become resistant to antibiotics. We tested whether there is still cooccurrence of resistance and virulence determinants. We performed a comparative study of environmental and human gut metagenomes from different individuals and from distinct human populations across the world. We found a great diversity of antibiotic resistance determinants (AR diversity [ARd]) and virulence factors (VF diversity [VFd]) in metagenomes. Importantly there is a correlation between ARd and VFd, even after correcting for protein family richness. In the human gut, there are less ARd and VFd than in more diversified environments, and yet correlations between the ARd and VFd are stronger. They can vary from very high in Malawi, where antibiotic consumption is unattended, to nonexistent in the uncontacted Amerindian population. We conclude that there is cooccurrence of resistance and virulence determinants in human gut microbiomes, suggesting a possible coselective mechanism.IMPORTANCE Every year, thousands of tons of antibiotics are used, not only in human and animal health but also as growth promoters in livestock. Consequently, during the last 75 years, antibiotic-resistant bacterial strains have been selected in human and environmental microbial communities. This implies that, even when pathogenic bacteria are the targets of antibiotics, hundreds of nonpathogenic bacterial species are also affected. Here, we performed a comparative study of environmental and human gut microbial communities issuing from different individuals and from distinct human populations across the world. We found that antibiotic resistance and pathogenicity are correlated and speculate that, by selecting for resistant bacteria, we may be selecting for more virulent strains as a side effect of antimicrobial therapy.

Antibiotics as both friends and foes of the human gut microbiome: The microbial community approach

The exposure of the human gut to antibiotics can have a great impact on human health. Antibiotics pertain to the preservation of human health and are useful tools for fighting bacterial infections. They can be used for curing infections and can play a critical role in immunocompromised or chronic patients, or in fighting childhood severe malnutrition. Yet, the genomic and phylogenetic diversity of the human gut changes under antibiotic exposure. Antibiotics can also have severe side effects on human gut health, due to the spreading of potential antibiotic resistance genetic traits and to their correlation with virulence of some bacterial pathogens. They can shape, and even disrupt, the composition and functioning diversity of the human gut microbiome. Traditionally bacterial antibiotic resistances have been evaluated at clone or population level. However, the understanding of these two apparently disparate perspectives as both friends and foes may come from the study of microbiomes as a whole and from the evaluation of both positive and negative effects of antibiotics on microbial community dynamics and diversity. In this review we present some metagenomic tools and databases that enable the studying of antibiotic resistance in human gut metagenomes, promoting the development of personalized medicine strategies, new antimicrobial therapy protocols and patient follow-up.

Position Matters: Network Centrality Considerably Impacts Rates of Protein Evolution in the Human Protein-Protein Interaction Network

The proteins of any organism evolve at disparate rates. A long list of factors affecting rates of protein evolution have been identified. However, the relative importance of each factor in determining rates of protein evolution remains unresolved. The prevailing view is that evolutionary rates are dominantly determined by gene expression, and that other factors such as network centrality have only a marginal effect, if any. However, this view is largely based on analyses in yeasts, and accurately measuring the importance of the determinants of rates of protein evolution is complicated by the fact that the different factors are often correlated with each other, and by the relatively poor quality of available functional genomics data sets. Here, we use correlation, partial correlation and principal component regression analyses to measure the contributions of several factors to the variability of the rates of evolution of human proteins. For this purpose, we analyzed the entire human protein–protein interaction data set and the human signal transduction network—a network data set of exceptionally high quality, obtained by manual curation, which is expected to be virtually free from false positives. In contrast with the prevailing view, we observe that network centrality (measured as the number of physical and nonphysical interactions, betweenness, and closeness) has a considerable impact on rates of protein evolution. Surprisingly, the impact of centrality on rates of protein evolution seems to be comparable, or even superior according to some analyses, to that of gene expression. Our observations seem to be independent of potentially confounding factors and from the limitations (biases and errors) of interactomic data sets.

Secreted Proteins Defy the Expression Level-Evolutionary Rate Anticorrelation

The rates of evolution of the proteins of any organism vary across orders of magnitude. A primary factor influencing rates of protein evolution is expression. A strong negative correlation between expression levels and evolutionary rates (the so-called E-R anticorrelation) has been observed in virtually all studied organisms. This effect is currently attributed to the abundance-dependent fitness costs of misfolding and unspecific protein-protein interactions, among other factors. Secreted proteins are folded in the endoplasmic reticulum, a compartment where chaperones, folding catalysts, and stringent quality control mechanisms promote their correct folding and may reduce the fitness costs of misfolding. In addition, confinement of secreted proteins to the extracellular space may reduce misinteractions and their deleterious effects. We hypothesize that each of these factors (the secretory pathway quality control and extracellular location) may reduce the strength of the E-R anticorrelation. Indeed, here we show that among human proteins that are secreted to the extracellular space, rates of evolution do not correlate with protein abundances. This trend is robust to controlling for several potentially confounding factors and is also observed when analyzing protein abundance data for 6 human tissues. In addition, analysis of mRNA abundance data for 32 human tissues shows that the E-R correlation is always less negative, and sometimes nonsignificant, in secreted proteins. Similar observations were made in Caenorhabditis elegans and in Escherichia coli, and to a lesser extent in Drosophila melanogaster, Saccharomyces cerevisiae and Arabidopsis thaliana. Our observations contribute to understand the causes of the E-R anticorrelation.

Abundance and co-occurrence of extracellular capsules increase environmental breadth: Implications for the emergence of pathogens

Extracellular capsules constitute the outermost layer of many bacteria, are major virulence factors, and affect antimicrobial therapies. They have been used as epidemiological markers and recently became vaccination targets. Despite the efforts to biochemically serotype capsules in a few model pathogens, little is known of their taxonomic and environmental distribution. We developed, validated, and made available a computational tool, CapsuleFinder, to identify capsules in genomes. The analysis of over 2500 prokaryotic genomes, accessible in a database, revealed that ca. 50% of them—including Archaea—encode a capsule. The Wzx/Wzy-dependent capsular group was by far the most abundant. Surprisingly, a fifth of the genomes encode more than one capsule system—often from different groups—and their non-random co-occurrence suggests the existence of negative and positive epistatic interactions. To understand the role of multiple capsules, we queried more than 6700 metagenomes for the presence of species encoding capsules and showed that their distribution varied between environmental categories and, within the human microbiome, between body locations. Species encoding capsules, and especially those encoding multiple capsules, had larger environmental breadths than the other species. Accordingly, capsules were more frequent in environmental bacteria than in pathogens and, within the latter, they were more frequent among facultative pathogens. Nevertheless, capsules were frequent in clinical samples, and were usually associated with fast-growing bacteria with high infectious doses. Our results suggest that capsules increase the environmental range of bacteria and make them more resilient to environmental perturbations. Capsules might allow opportunistic pathogens to profit from empty ecological niches or environmental perturbations, such as those resulting from antibiotic therapy, to colonize the host. Capsule-associated virulence might thus be a by-product of environmental adaptation. Understanding the role of capsules in natural environments might enlighten their function in pathogenesis.

Extracellular capsules protect bacterial cells from external aggressions such as antibiotics or desiccation, but can also be targeted by vaccines. Since little was known about their frequency across Prokaryotes, we created and made freely available a computational tool, CapsuleFinder, to identify them from genomic data. Surprisingly, its use showed that many bacterial strains, especially those with the largest genomes, encode several capsules. The frequencies of the different combinations of capsule groups depended strongly on the phyla and the groups themselves, suggesting the existence of epistatic interactions between capsules. Bacteria encoding capsule systems were found in many natural environments, and were frequent in the human microbiome. In contrast to their frequent association with virulence, we found many more capsules in non-pathogens or facultative pathogens than among obligatory pathogens. We suggest that capsules increase the environmental breadth of bacteria thereby facilitating host colonization by opportunistic pathogens.

Unexpected genomic features in widespread intracellular bacteria: evidence for motility of marine chlamydiae

Chlamydiae are obligate intracellular bacteria comprising important human pathogens and symbionts of protists. Molecular evidence indicates a tremendous diversity of chlamydiae particularly in marine environments, yet our current knowledge is based mainly on terrestrial representatives. Here we provide first insights into the biology of marine chlamydiae representing three divergent clades. Our analysis of single-cell amplified genomes revealed hallmarks of the chlamydial lifestyle, supporting the ancient origin of their characteristic developmental cycle and major virulence mechanisms. Surprisingly, these chlamydial genomes encode a complete flagellar apparatus, a previously unreported feature. We show that flagella are an ancient trait that was subject to differential gene loss among extant chlamydiae. Together with a chemotaxis system, these marine chlamydiae are likely motile, with flagella potentially playing a role during host cell infection. This study broadens our view on chlamydial biology and indicates a largely underestimated potential to adapt to different hosts and environments.

Identification of protein secretion systems in bacterial genomes

Bacteria with two cell membranes (diderms) have evolved complex systems for protein secretion. These systems were extensively studied in some model bacteria, but the characterisation of their diversity has lagged behind due to lack of standard annotation tools. We built online and standalone computational tools to accurately predict protein secretion systems and related appendages in bacteria with LPS-containing outer membranes. They consist of models describing the systems’ components and genetic organization to be used with MacSyFinder to search for T1SS-T6SS, T9SS, flagella, Type IV pili and Tad pili. We identified ~10,000 candidate systems in bacterial genomes, where T1SS and T5SS were by far the most abundant and widespread. All these data are made available in a public database. The recently described T6SSⁱⁱⁱ and T9SS were restricted to Bacteroidetes, and T6SSⁱⁱ to Francisella. The T2SS, T3SS, and T4SS were frequently encoded in single-copy in one locus, whereas most T1SS were encoded in two loci. The secretion systems of diderm Firmicutes were similar to those found in other diderms. Novel systems may remain to be discovered, since some clades of environmental bacteria lacked all known protein secretion systems. Our models can be fully customized, which should facilitate the identification of novel systems.

Comparative genomic analysis of evolutionarily conserved but functionally uncharacterized membrane proteins in archaea: Prediction of novel components of secretion, membrane remodeling and glycosylation systems

A systematic comparative genomic analysis of all archaeal membrane proteins that have been projected to the last archaeal common ancestor gene set led to the identification of several novel components of predicted secretion, membrane remodeling, and protein glycosylation systems. Among other findings, most crenarchaea have been shown to encode highly diverged orthologs of the membrane insertase YidC, which is nearly universal in bacteria, eukaryotes, and euryarchaea. We also identified a vast family of archaeal proteins, including the C-terminal domain of N-glycosylation protein AglD, as membrane flippases homologous to the flippase domain of bacterial multipeptide resistance factor MprF, a bifunctional lysylphosphatidylglycerol synthase and flippase. Additionally, several proteins were predicted to function as membrane transporters. The results of this work, combined with our previous analyses, reveal an unexpected diversity of putative archaeal membrane-associated functional systems that remain to be functionally characterized. A more general conclusion from this work is that the currently available collection of archaeal (and bacterial) genomes could be sufficient to identify (almost) all widespread functional modules and develop experimentally testable predictions of their functions.

Genes associated with ant social behavior show distinct transcriptional and evolutionary patterns

Studies of the genetic basis and evolution of complex social behavior emphasize either conserved or novel genes. To begin to reconcile these perspectives, we studied how the evolutionary conservation of genes associated with social behavior depends on regulatory context, and whether genes associated with social behavior exist in distinct regulatory and evolutionary contexts. We identified modules of co-expressed genes associated with age-based division of labor between nurses and foragers in the ant Monomorium pharaonis, and we studied the relationship between molecular evolution, connectivity, and expression. Highly connected and expressed genes were more evolutionarily conserved, as expected. However, compared to the rest of the genome, forager-upregulated genes were much more highly connected and conserved, while nurse-upregulated genes were less connected and more evolutionarily labile. Our results indicate that the genetic architecture of social behavior includes both highly connected and conserved components as well as loosely connected and evolutionarily labile components.

DOI:http://dx.doi.org/10.7554/eLife.04775.001

Animal species vary widely in their degree of social behavior. Some species live solitarily, and others, such as ants and humans, form large societies. Many researchers have tried to understand the genetic changes underlying the evolution of social behavior. Some researchers suggest that it involves recycling existing genes that also have other conserved functions. Others propose that the evolution of social behavior involves completely new genes that are not found in related but solitary species.

Ants are one of the best-studied social animals. An established colony can contain many 1000s of individuals that live and work together and perform different roles. The queen's job is to lay eggs, while the worker ants do everything else, including collecting food, caring for the young, and protecting the colony. In some species of ant—including the pharaoh ant—a worker's role changes as it ages. Younger workers tend to stay in the nest and nurse the brood, while older workers tend to leave the nest and forage for food.

Mikheyev and Linksvayer asked: which genes are responsible for this age-based division of labor? And how did this aspect of social behavior evolve? First, after observing pharaoh ants from two colonies set up in the laboratory, they confirmed that workers nursing the brood were on average almost a week younger than those seen collecting food. Next Mikheyev and Linksvayer identified which genes were expressed in ants of different ages, or ants engaged in different tasks. Specific sets of genes were expressed more (or ‘up-regulated’) in nurse workers, while others were up-regulated in foraging workers.

Mikheyev and Linksvayer then investigated how rapidly these genes had evolved by comparing them to related genes found in other social insects (fire ants and honey bees). They also determined the ‘connectivity’ of these genes by asking how many other genes showed similar expression patterns. In many organisms, how rapidly a gene evolves depends on how tightly connected its expression is to the expression of other genes; highly connected genes evolve more slowly.

The genes that were expressed more in the older foraging workers were both more highly connected and more evolutionarily conserved in the other social insects. Genes that were up-regulated in the younger nurse workers were more loosely connected and rapidly evolving.

Mikheyev and Linksvayer's findings show that the evolution of social behavior in animals involves both new genes, which tend to be loosely connected, and conserved genes, which tend to be more highly connected.

DOI:http://dx.doi.org/10.7554/eLife.04775.002

Genome comparison of three serovar 5 pathogenic strains of Haemophilus parasuis: insights into an evolving swine pathogen

Haemophilus parasuis is the causative agent of Glässer’s disease, a systemic disorder characterized by polyarthritis, polyserositis and meningitis in pigs. Although it is well known that H. parasuis serovar 5 is the most prevalent serovar associated with the disease, the genetic differences among strains are only now being discovered. Genomes from two serovar 5 strains, SH0165 and 29755, are already available. Here, we present the draft genome of a third H. parasuis serovar 5 strain, the formal serovar 5 reference strain Nagasaki. An in silico genome subtractive analysis with full-length predicted genes of the three H. parasuis serovar 5 strains detected 95, 127 and 95 strain-specific genes (SSGs) for Nagasaki, SH0165 and 29755, respectively. We found that the genomic diversity within these three strains was high, in part because of a high number of mobile elements. Furthermore, a detailed analysis of large sequence polymorphisms (LSPs), encompassing regions ranging from 2 to 16 kb, revealed LSPs in virulence-related elements, such as a Toll-IL receptor, the AcrA multidrug efflux protein, an ATP-binding cassette (ABC) transporter, lipopolysaccharide-synthetizing enzymes and a tripartite ATP-independent periplasmic (TRAP) transporter. The whole-genome codon adaptation index (CAI) was also calculated and revealed values similar to other well-known bacterial pathogens. In addition, whole-genome SNP analysis indicated that nucleotide changes tended to be increased in membrane-related genes. This analysis provides further evidence that the genome of H. parasuis has been subjected to multiple lateral gene transfers (LGTs) and to fine-tuning of virulence factors, and has the potential for accelerated genome evolution.

Massive expansion of Ubiquitination-related gene families within the Chlamydiae

Gene loss, gain, and transfer play an important role in shaping the genomes of all organisms; however, the interplay of these processes in isolated populations, such as in obligate intracellular bacteria, is less understood. Despite a general trend towards genome reduction in these microbes, our phylogenomic analysis of the phylum Chlamydiae revealed that within the family Parachlamydiaceae, gene family expansions have had pronounced effects on gene content. We discovered that the largest gene families within the phylum are the result of rapid gene birth-and-death evolution. These large gene families are comprised of members harboring eukaryotic-like ubiquitination-related domains, such as F-box and BTB-box domains, marking the largest reservoir of these proteins found among bacteria. A heterologous type III secretion system assay suggests that these proteins function as effectors manipulating the host cell. The large disparity in copy number of members in these families between closely related organisms suggests that nonadaptive processes might contribute to the evolution of these gene families. Gene birth-and-death evolution in concert with genomic drift might represent a previously undescribed mechanism by which isolated bacterial populations diversify.

Targeting virulence: can we make evolution-proof drugs?

Antivirulence drugs are a new type of therapeutic drug that target virulence factors, potentially revitalising the drug-development pipeline with new targets. As antivirulence drugs disarm the pathogen, rather than kill or halt pathogen growth, it has been hypothesized that they will generate much weaker selection for resistance than traditional antibiotics. However, recent studies have shown that mechanisms of resistance to antivirulence drugs exist, seemingly damaging the 'evolution-proof' claim. In this Opinion article, we highlight a crucial distinction between whether resistance can emerge and whether it will spread to a high frequency under drug selection. We argue that selection for resistance can be reduced, or even reversed, using appropriate combinations of target and treatment environment, opening a path towards the development of evolutionarily robust novel therapeutics.

Transcriptional abundance is not the single force driving the evolution of bacterial proteins

Background

Despite rapid progress in understanding the mechanisms that shape the evolution of proteins, the relative importance of various factors remain to be elucidated. In this study, we have assessed the effects of 16 different biological features on the evolutionary rates (ERs) of protein-coding sequences in bacterial genomes.

Results

Our analysis of 18 bacterial species revealed new correlations between ERs and constraining factors. Previous studies have suggested that transcriptional abundance overwhelmingly constrains the evolution of yeast protein sequences. This transcriptional abundance leads to selection against misfolding or misinteractions. In this study we found that there was no single factor in determining the evolution of bacterial proteins. Not only transcriptional abundance (codon adaptation index and expression level), but also protein-protein associations (PPAs), essentiality (ESS), subcellular localization of cytoplasmic membrane (SLM), transmembrane helices (TMH) and hydropathicity score (HS) independently and significantly affected the ERs of bacterial proteins. In some species, PPA and ESS demonstrate higher correlations with ER than transcriptional abundance.

Conclusions

Different forces drive the evolution of protein sequences in yeast and bacteria. In bacteria, the constraints are involved in avoiding a build-up of toxic molecules caused by misfolding/misinteraction (transcriptional abundance), while retaining important functions (ESS, PPA) and maintaining the cell membrane (SLM, TMH and HS). Each of these independently contributes to the variation in protein evolution.

Insight into the evolution of the histidine triad protein (HTP) family in Streptococcus

The Histidine Triad Proteins (HTPs), also known as Pht proteins in Streptococcus pneumoniae, constitute a family of surface-exposed proteins that exist in many pathogenic streptococcal species. Although many studies have revealed the importance of HTPs in streptococcal physiology and pathogenicity, little is known about their origin and evolution. In this study, after identifying all htp homologs from 105 streptococcal genomes representing 38 different species/subspecies, we analyzed their domain structures, positions in genome, and most importantly, their evolutionary histories. By further projecting this information onto the streptococcal phylogeny, we made several major findings. First, htp genes originated earlier than the Streptococcus genus and gene-loss events have occurred among three streptococcal groups, resulting in the absence of the htp gene in the Bovis, Mutans and Salivarius groups. Second, the copy number of htp genes in other groups of Streptococcus is variable, ranging from one to four functional copies. Third, both phylogenetic evidence and domain structure analyses support the division of two htp subfamilies, designated as htp I and htp II. Although present mainly in the pyogenic group and in Streptococcus suis, htp II members are distinct from htp I due to the presence of an additional leucine-rich-repeat domain at the C-terminus. Finally, htp genes exhibit a faster nucleotide substitution rate than do housekeeping genes. Specifically, the regions outside the HTP domains are under strong positive selection. This distinct evolutionary pattern likely helped Streptococcus to easily escape from recognition by host immunity.