vCOMBAT: a novel tool to create and visualize a computational model of bacterial antibiotic target-binding.

BACKGROUND: As antibiotic resistance creates a significant global health threat, we need not only to accelerate the development of novel antibiotics but also to develop better treatment strategies using existing drugs to improve their efficacy and prevent the selection of further resistance. We require new tools to rationally design dosing regimens from data collected in early phases of antibiotic and dosing development. Mathematical models such as mechanistic pharmacodynamic drug-target binding explain mechanistic details of how the given drug concentration affects its targeted bacteria. However, there are no available tools in the literature that allow non-quantitative scientists to develop computational models to simulate antibiotic-target binding and its effects on bacteria. RESULTS: In this work, we have devised an extension of a mechanistic binding-kinetic model to incorporate clinical drug concentration data. Based on the extended model, we develop a novel and interactive web-based tool that allows non-quantitative scientists to create and visualize their own computational models of bacterial antibiotic target-binding based on their considered drugs and bacteria. We also demonstrate how Rifampicin affects bacterial populations of Tuberculosis bacteria using our vCOMBAT tool. CONCLUSIONS: The vCOMBAT online tool is publicly available at https://combat-bacteria.org/ .

Shigella sonnei infection of zebrafish reveals that O-antigen mediates neutrophil tolerance and dysentery incidence.

Shigella flexneri is historically regarded as the primary agent of bacillary dysentery, yet the closely-related Shigella sonnei is replacing S. flexneri, especially in developing countries. The underlying reasons for this dramatic shift are mostly unknown. Using a zebrafish (Danio rerio) model of Shigella infection, we discover that S. sonnei is more virulent than S. flexneri in vivo. Whole animal dual-RNAseq and testing of bacterial mutants suggest that S. sonnei virulence depends on its O-antigen oligosaccharide (which is unique among Shigella species). We show in vivo using zebrafish and ex vivo using human neutrophils that S. sonnei O-antigen can mediate neutrophil tolerance. Consistent with this, we demonstrate that O-antigen enables S. sonnei to resist phagolysosome acidification and promotes neutrophil cell death. Chemical inhibition or promotion of phagolysosome maturation respectively decreases and increases neutrophil control of S. sonnei and zebrafish survival. Strikingly, larvae primed with a sublethal dose of S. sonnei are protected against a secondary lethal dose of S. sonnei in an O-antigen-dependent manner, indicating that exposure to O-antigen can train the innate immune system against S. sonnei. Collectively, these findings reveal O-antigen as an important therapeutic target against bacillary dysentery, and may explain the rapidly increasing S. sonnei burden in developing countries.

Horizontal DNA Transfer Mechanisms of Bacteria as Weapons of Intragenomic Conflict.

Horizontal DNA transfer (HDT) is a pervasive mechanism of diversification in many microbial species, but its primary evolutionary role remains controversial. Much recent research has emphasised the adaptive benefit of acquiring novel DNA, but here we argue instead that intragenomic conflict provides a coherent framework for understanding the evolutionary origins of HDT. To test this hypothesis, we developed a mathematical model of a clonally descended bacterial population undergoing HDT through transmission of mobile genetic elements (MGEs) and genetic transformation. Including the known bias of transformation toward the acquisition of shorter alleles into the model suggested it could be an effective means of counteracting the spread of MGEs. Both constitutive and transient competence for transformation were found to provide an effective defence against parasitic MGEs; transient competence could also be effective at permitting the selective spread of MGEs conferring a benefit on their host bacterium. The coordination of transient competence with cell-cell killing, observed in multiple species, was found to result in synergistic blocking of MGE transmission through releasing genomic DNA for homologous recombination while simultaneously reducing horizontal MGE spread by lowering the local cell density. To evaluate the feasibility of the functions suggested by the modelling analysis, we analysed genomic data from longitudinal sampling of individuals carrying Streptococcus pneumoniae. This revealed the frequent within-host coexistence of clonally descended cells that differed in their MGE infection status, a necessary condition for the proposed mechanism to operate. Additionally, we found multiple examples of MGEs inhibiting transformation through integrative disruption of genes encoding the competence machinery across many species, providing evidence of an ongoing "arms race." Reduced rates of transformation have also been observed in cells infected by MGEs that reduce the concentration of extracellular DNA through secretion of DNases. Simulations predicted that either mechanism of limiting transformation would benefit individual MGEs, but also that this tactic's effectiveness was limited by competition with other MGEs coinfecting the same cell. A further observed behaviour we hypothesised to reduce elimination by transformation was MGE activation when cells become competent. Our model predicted that this response was effective at counteracting transformation independently of competing MGEs. Therefore, this framework is able to explain both common properties of MGEs, and the seemingly paradoxical bacterial behaviours of transformation and cell-cell killing within clonally related populations, as the consequences of intragenomic conflict between self-replicating chromosomes and parasitic MGEs. The antagonistic nature of the different mechanisms of HDT over short timescales means their contribution to bacterial evolution is likely to be substantially greater than previously appreciated.

Efficient Inference of Recent and Ancestral Recombination within Bacterial Populations.

Prokaryotic evolution is affected by horizontal transfer of genetic material through recombination. Inference of an evolutionary tree of bacteria thus relies on accurate identification of the population genetic structure and recombination-derived mosaicism. Rapidly growing databases represent a challenge for computational methods to detect recombinations in bacterial genomes. We introduce a novel algorithm called fastGEAR which identifies lineages in diverse microbial alignments, and recombinations between them and from external origins. The algorithm detects both recent recombinations (affecting a few isolates) and ancestral recombinations between detected lineages (affecting entire lineages), thus providing insight into recombinations affecting deep branches of the phylogenetic tree. In simulations, fastGEAR had comparable power to detect recent recombinations and outstanding power to detect the ancestral ones, compared with state-of-the-art methods, often with a fraction of computational cost. We demonstrate the utility of the method by analyzing a collection of 616 whole-genomes of a recombinogenic pathogen Streptococcus pneumoniae, for which the method provided a high-resolution view of recombination across the genome. We examined in detail the penicillin-binding genes across the Streptococcus genus, demonstrating previously undetected genetic exchanges between different species at these three loci. Hence, fastGEAR can be readily applied to investigate mosaicism in bacterial genes across multiple species. Finally, fastGEAR correctly identified many known recombination hotspots and pointed to potential new ones. Matlab code and Linux/Windows executables are available at https://users.ics.aalto.fi/~pemartti/fastGEAR/ (last accessed February 6, 2017).

Pneumococcal Capsule Synthesis Locus cps as Evolutionary Hotspot with Potential to Generate Novel Serotypes by Recombination.

Diversity of the polysaccharide capsule in Streptococcus pneumoniae-main surface antigen and the target of the currently used pneumococcal vaccines-constitutes a major obstacle in eliminating pneumococcal disease. Such diversity is genetically encoded by almost 100 variants of the capsule biosynthesis locus, cps. However, the evolutionary dynamics of the capsule remains not fully understood. Here, using genetic data from 4,519 bacterial isolates, we found cps to be an evolutionary hotspot with elevated substitution and recombination rates. These rates were a consequence of relaxed purifying selection and positive, diversifying selection acting at this locus, supporting the hypothesis that the capsule has an increased potential to generate novel diversity compared with the rest of the genome. Diversifying selection was particularly evident in the region of wzd/wze genes, which are known to regulate capsule expression and hence the bacterium's ability to cause disease. Using a novel, capsule-centered approach, we analyzed the evolutionary history of 12 major serogroups. Such analysis revealed their complex diversification scenarios, which were principally driven by recombination with other serogroups and other streptococci. Patterns of recombinational exchanges between serogroups could not be explained by serotype frequency alone, thus pointing to nonrandom associations between co-colonizing serotypes. Finally, we discovered a previously unobserved mosaic serotype 39X, which was confirmed to carry a viable and structurally novel capsule. Adding to previous discoveries of other mosaic capsules in densely sampled collections, these results emphasize the strong adaptive potential of the bacterium by its ability to generate novel antigenic diversity by recombination.

Heterogeneity in the frequency and characteristics of homologous recombination in pneumococcal evolution.

The bacterium Streptococcus pneumoniae (pneumococcus) is one of the most important human bacterial pathogens, and a leading cause of morbidity and mortality worldwide. The pneumococcus is also known for undergoing extensive homologous recombination via transformation with exogenous DNA. It has been shown that recombination has a major impact on the evolution of the pathogen, including acquisition of antibiotic resistance and serotype-switching. Nevertheless, the mechanism and the rates of recombination in an epidemiological context remain poorly understood. Here, we proposed several mathematical models to describe the rate and size of recombination in the evolutionary history of two very distinct pneumococcal lineages, PMEN1 and CC180. We found that, in both lineages, the process of homologous recombination was best described by a heterogeneous model of recombination with single, short, frequent replacements, which we call micro-recombinations, and rarer, multi-fragment, saltational replacements, which we call macro-recombinations. Macro-recombination was associated with major phenotypic changes, including serotype-switching events, and thus was a major driver of the diversification of the pathogen. We critically evaluate biological and epidemiological processes that could give rise to the micro-recombination and macro-recombination processes.

Evolution of stress response in the face of unreliable environmental signals.

Most organisms live in ever-changing environments, and have to cope with a range of different conditions. Often, the set of biological traits that are needed to grow, reproduce, and survive varies between conditions. As a consequence, organisms have evolved sensory systems to detect environmental signals, and to modify the expression of biological traits in response. However, there are limits to the ability of such plastic responses to cope with changing environments. Sometimes, environmental shifts might occur suddenly, and without preceding signals, so that organisms might not have time to react. Other times, signals might be unreliable, causing organisms to prepare themselves for changes that then do not occur. Here, we focus on such unreliable signals that indicate the onset of adverse conditions. We use analytical and individual-based models to investigate the evolution of simple rules that organisms use to decide whether or not to switch to a protective state. We find evolutionary transitions towards organisms that use a combination of random switching and switching in response to the signal. We also observe that, in spatially heterogeneous environments, selection on the switching strategy depends on the composition of the population, and on population size. These results are in line with recent experiments that showed that many unicellular organisms can attain different phenotypic states in a probabilistic manner, and lead to testable predictions about how this could help organisms cope with unreliable signals.

Estimating the fitness cost of escape from HLA presentation in HIV-1 protease and reverse transcriptase.

Human immunodeficiency virus (HIV-1) is, like most pathogens, under selective pressure to escape the immune system of its host. In particular, HIV-1 can avoid recognition by cytotoxic T lymphocytes (CTLs) by altering the binding affinity of viral peptides to human leukocyte antigen (HLA) molecules, the role of which is to present those peptides to the immune system. It is generally assumed that HLA escape mutations carry a replicative fitness cost, but these costs have not been quantified. In this study, we assess the replicative cost of mutations which are likely to escape presentation by HLA molecules in the region of HIV-1 protease and reverse transcriptase. Specifically, we combine computational approaches for prediction of in vitro replicative fitness and peptide binding affinity to HLA molecules. We find that mutations which impair binding to HLA-A molecules tend to have lower in vitro replicative fitness than mutations which do not impair binding to HLA-A molecules, suggesting that HLA-A escape mutations carry higher fitness costs than non-escape mutations. We argue that the association between fitness and HLA-A binding impairment is probably due to an intrinsic cost of escape from HLA-A molecules, and these costs are particularly strong for HLA-A alleles associated with efficient virus control. Counter-intuitively, we do not observe a significant effect in the case of HLA-B, but, as discussed, this does not argue against the relevance of HLA-B in virus control. Overall, this article points to the intriguing possibility that HLA-A molecules preferentially target more conserved regions of HIV-1, emphasizing the importance of HLA-A genes in the evolution of HIV-1 and RNA viruses in general.

Ongoing shuffling of protein fragments diversifies core viral functions linked to interactions with bacterial hosts.

Biological modularity enhances evolutionary adaptability. This principle is vividly exemplified by bacterial viruses (phages), which display extensive genomic modularity. Phage genomes are composed of independent functional modules that evolve separately and recombine in various configurations. While genomic modularity in phages has been extensively studied, less attention has been paid to protein modularity-proteins consisting of distinct building blocks that can evolve and recombine, enhancing functional and genetic diversity. Here, we use a set of 133,574 representative phage proteins and highly sensitive homology detection to capture instances of domain mosaicism, defined as fragment sharing between two otherwise unrelated proteins, and to understand its relationship with functional diversity in phage genomes. We discover that unrelated proteins from diverse functional classes frequently share homologous domains. This phenomenon is particularly pronounced within receptor-binding proteins, endolysins, and DNA polymerases. We also identify multiple instances of recent diversification via domain shuffling in receptor-binding proteins, neck passage structures, endolysins and some members of the core replication machinery, often transcending distant taxonomic and ecological boundaries. Our findings suggest that ongoing diversification via domain shuffling is reflective of a co-evolutionary arms race, driven by the need to overcome various bacterial resistance mechanisms against phages.

Non-genetic inheritance and the patterns of antagonistic coevolution.

BACKGROUND: Antagonistic species interactions can lead to coevolutionary genotype or phenotype frequency oscillations, with important implications for ecological and evolutionary processes. However, direct empirical evidence of such oscillations is rare. The rarity of observations is generally attributed to inherent difficulties of ecological and evolutionary long-term studies, to weak or absent interaction between species, or to the absence of negative frequency-dependence. RESULTS: Here, we show that another factor - non-genetic inheritance, mediated for example by epigenetic mechanisms - can completely eliminate oscillations in the presence of such negative frequency dependence, even if only a small fraction of offspring are affected. We analytically derive the threshold value of this fraction at which the dynamics change from oscillatory to stable, and investigate how selection, mutation and generation times differences between the two species affect the threshold value. These results strongly suggest that the lack of phenotype frequency oscillations should not be attributed to the lack of strong interactions between antagonistic species. CONCLUSIONS: Given increasing evidence of non-genetic effects on the outcomes of antagonistic species interactions, we suggest that these effects should be incorporated into ecological and evolutionary models of interacting species.

The impact of environmental change on host-parasite coevolutionary dynamics.

Environmental factors are known to affect the strength and the specificity of interactions between hosts and parasites. However, how this shapes patterns of coevolutionary dynamics is not clear. Here, we construct a simple mathematical model to study the effect of environmental change on host-parasite coevolutionary outcome when interactions are of the matching-alleles or the gene-for-gene type. Environmental changes may effectively alter the selective pressure and the level of specialism in the population. Our results suggest that environmental change altering the specificity of selection in antagonistic interactions can produce alternating time windows of cyclical allele-frequency dynamics and cessation thereof. This type of environmental impact can also explain the maintenance of polymorphism in gene-for-gene interactions without costs. Overall, our study points to the potential consequences of environmental variation in coevolution, and thus the importance of characterizing genotype-by-genotype-by-environment interactions in natural host-parasite systems, especially those that change the direction of selection acting between the two species.

Diversity and evolution of surface polysaccharide synthesis loci in Enterobacteriales.

Bacterial capsules and lipopolysaccharides are diverse surface polysaccharides (SPs) that serve as the frontline for interactions with the outside world. While SPs can evolve rapidly, their diversity and evolutionary dynamics across different taxonomic scales has not been investigated in detail. Here, we focused on the bacterial order Enterobacteriales (including the medically relevant Enterobacteriaceae), to carry out comparative genomics of two SP locus synthesis regions, cps and kps, using 27,334 genomes from 45 genera. We identified high-quality cps loci in 22 genera and kps in 11 genera, around 4% of which were detected in multiple species. We found SP loci to be highly dynamic genetic entities: their evolution was driven by high rates of horizontal gene transfer (HGT), both of whole loci and component genes, and relaxed purifying selection, yielding large repertoires of SP diversity. In spite of that, we found the presence of (near-)identical locus structures in distant taxonomic backgrounds that could not be explained by recent exchange, pointing to long-term selective preservation of locus structures in some populations. Our results reveal differences in evolutionary dynamics driving SP diversity within different bacterial species, with lineages of Escherichia coli, Enterobacter hormaechei and Klebsiella aerogenes most likely to share SP loci via recent exchange; and lineages of Salmonella enterica, Citrobacter sakazakii and Serratia marcescens most likely to share SP loci via other mechanisms such as long-term preservation. Overall, the evolution of SP loci in Enterobacteriales is driven by a range of evolutionary forces and their dynamics and relative importance varies between different species.

Horizontal gene transfer rate is not the primary determinant of observed antibiotic resistance frequencies in Streptococcus pneumoniae.

The extent to which evolution is constrained by the rate at which horizontal gene transfer (HGT) allows DNA to move between genetic lineages is an open question, which we address in the context of antibiotic resistance in Streptococcus pneumoniae. We analyze microbiological, genomic, and epidemiological data from the largest-to-date sequenced pneumococcal carriage study in 955 infants from a refugee camp on the Thailand-Myanmar border. Using a unified framework, we simultaneously test prior hypotheses on rates of HGT and a key evolutionary covariate (duration of carriage) as determinants of resistance frequencies. We conclude that in this setting, there is little evidence of HGT playing a major role in determining resistance frequencies. Instead, observed resistance frequencies are best explained as the outcome of selection acting on a pool of variants, irrespective of the rate at which resistance determinants move between genetic lineages.

Diversity-Generating Machines: Genetics of Bacterial Sugar-Coating.

Bacterial pathogens and commensals are surrounded by diverse surface polysaccharides which include capsules and lipopolysaccharides. These carbohydrates play a vital role in bacterial ecology and interactions with the environment. Here, we review recent rapid advancements in this field, which have improved our understanding of the roles, structures, and genetics of bacterial polysaccharide antigens. Genetic loci encoding the biosynthesis of these antigens may have evolved as bacterial diversity-generating machines, driven by selection from a variety of forces, including host immunity, bacteriophages, and cell-cell interactions. We argue that the high adaptive potential of polysaccharide antigens should be taken into account in the design of polysaccharide-targeting medical interventions like conjugate vaccines and phage-based therapies.

Evolution of the Staphylococcus argenteus ST2250 Clone in Northeastern Thailand Is Linked with the Acquisition of Livestock-Associated Staphylococcal Genes.

Staphylococcus argenteus is a newly named species previously described as a divergent lineage of Staphylococcus aureus that has recently been shown to have a global distribution. Despite growing evidence of the clinical importance of this species, knowledge about its population epidemiology and genomic architecture is limited. We used whole-genome sequencing to evaluate and compare S. aureus (n = 251) and S. argenteus (n = 68) isolates from adults with staphylococcal sepsis at several hospitals in northeastern Thailand between 2006 and 2013. The majority (82%) of the S. argenteus isolates were of multilocus sequence type 2250 (ST2250). S. aureus was more diverse, although 43% of the isolates belonged to ST121. Bayesian analysis suggested an S. argenteus ST2250 substitution rate of 4.66 (95% confidence interval [CI], 3.12 to 6.38) mutations per genome per year, which was comparable to the S. aureus ST121 substitution rate of 4.07 (95% CI, 2.61 to 5.55). S. argenteus ST2250 emerged in Thailand an estimated 15 years ago, which contrasts with the S. aureus ST1, ST88, and ST121 clades that emerged around 100 to 150 years ago. Comparison of S. argenteus ST2250 genomes from Thailand and a global collection indicated a single introduction into Thailand, followed by transmission to local and more distant countries in Southeast Asia and further afield. S. argenteus and S. aureus shared around half of their core gene repertoire, indicating a high level of divergence and providing strong support for their classification as separate species. Several gene clusters were present in ST2250 isolates but absent from the other S. argenteus and S. aureus study isolates. These included multiple exotoxins and antibiotic resistance genes that have been linked previously with livestock-associated S. aureus, consistent with a livestock reservoir for S. argenteus These genes appeared to be associated with plasmids and mobile genetic elements and may have contributed to the biological success of ST2250.IMPORTANCE In this study, we used whole-genome sequencing to understand the genome evolution and population structure of a systematic collection of ST2250 S. argenteus isolates. A newly identified ancestral species of S. aureus, S. argenteus has become increasingly known as a clinically important species that has been reported recently across various countries. Our results indicate that S. argenteus has spread at a relatively rapid pace over the past 2 decades across northeastern Thailand and acquired multiple exotoxin and antibiotic resistance genes that have been linked previously with livestock-associated S. aureus Our findings highlight the clinical importance and potential pathogenicity of S. argenteus as a recently emerging pathogen.

The role of recombination for the coevolutionary dynamics of HIV and the immune response.

The evolutionary implications of recombination in HIV remain not fully understood. A plausible effect could be an enhancement of immune escape from cytotoxic T lymphocytes (CTLs). In order to test this hypothesis, we constructed a population dynamic model of immune escape in HIV and examined the viral-immune dynamics with and without recombination. Our model shows that recombination (i) increases the genetic diversity of the viral population, (ii) accelerates the emergence of escape mutations with and without compensatory mutations, and (iii) accelerates the acquisition of immune escape mutations in the early stage of viral infection. We see a particularly strong impact of recombination in systems with broad, non-immunodominant CTL responses. Overall, our study argues for the importance of recombination in HIV in allowing the virus to adapt to changing selective pressures as imposed by the immune system and shows that the effect of recombination depends on the immunodominance pattern of effector T cell responses.

On the evolution of sexual reproduction in hosts coevolving with multiple parasites.

Host-parasite coevolution has been studied extensively in the context of the evolution of sex. Although hosts typically coevolve with several parasites, most studies considered one-host/one-parasite interactions. Here, we study population-genetic models in which hosts interact with two parasites. We find that host/multiple-parasite models differ nontrivially from host/single-parasite models. Selection for sex resulting from interactions with a single parasite is often outweighed by detrimental effects due to the interaction between parasites if coinfection affects the host more severely than expected based on single infections, and/or if double infections are more common than expected based on single infections. The resulting selection against sex is caused by strong linkage-disequilibria of constant sign that arise between host loci interacting with different parasites. In contrast, if coinfection affects hosts less severely than expected and double infections are less common than expected, selection for sex due to interactions with individual parasites can now be reinforced by additional rapid linkage-disequilibrium oscillations with changing sign. Thus, our findings indicate that the presence of an additional parasite can strongly affect the evolution of sex in ways that cannot be predicted from single-parasite models, and that thus host/multiparasite models are an important extension of the Red Queen Hypothesis.