How infection-triggered pathobionts influence virulence evolution.

Host-pathogen interactions can be influenced by the host microbiota, as the microbiota can facilitate or prevent pathogen infections. In addition, members of the microbiota can become virulent. Such pathobionts can cause co-infections when a pathogen infection alters the host immune system and triggers dysbiosis. Here we performed a theoretical investigation of how pathobiont co-infections affect the evolution of pathogen virulence. We explored the possibility that the likelihood of pathobiont co-infection depends on the evolving virulence of the pathogen. We found that, in contrast to the expectation from classical theory, increased virulence is not always selected for. For an increasing likelihood of co-infection with increasing pathogen virulence, we found scenario-specific selection for either increased or decreased virulence. Evolutionary changes, however, in pathogen virulence do not always translate into similar changes in combined virulence of the pathogen and the pathobiont. Only in one of the scenarios where pathobiont co-infection is triggered above a specific virulence level we found a reduction in combined virulence. This was not the case when the probability of pathobiont co-infection linearly increased with pathogen virulence. Taken together, our study draws attention to the possibility that host-microbiota interactions can be both the driver and the target of pathogen evolution. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

Correcting for Antibody Waning in Cumulative Incidence Estimation from Sequential Serosurveys.

Serosurveys are a widely used tool to estimate the cumulative incidence, i.e. the fraction of a population that have been infected by a given pathogen. These surveys rely on serological assays that measure the level of pathogen-specific antibodies. Because antibody levels are waning, the fraction of previously infected individuals that have sero-reverted increases with time past infection. To avoid underestimating the true cumulative incidence, it is therefore essential to correct for waning antibody levels. We present an empirically-supported approach for sero-reversion correction in cumulative incidence estimation when sequential serosurveys are conducted in the context of a newly emerging infectious disease. The correction is based on the observed dynamics of antibody titers in sero-positive cases and validated using several in silico test scenarios. Furthermore, through this approach we revise a previous cumulative incidence estimate, which relies on the assumption of an exponentially-declining probability of sero-reversion over time, of SARS-CoV-2 of 76% in Manaus, Brazil, by October 2020 to 47.6% (43.5% - 53.5%). This estimate has implications e.g. for the proximity to herd immunity in Manaus in late 2020.

The pharmacokinetic-pharmacodynamic modelling framework as a tool to predict drug resistance evolution.

Pharmacokinetic-pharmacodynamic (PKPD) models, which describe how drug concentrations change over time and how that affects pathogen growth, have proven highly valuable in designing optimal drug treatments aimed at bacterial eradication. However, the fast rise of antimicrobial resistance calls for increased focus on an additional treatment optimization criterion: avoidance of resistance evolution. We demonstrate here how coupling PKPD and population genetics models can be used to determine treatment regimens that minimize the potential for antimicrobial resistance evolution. Importantly, the resulting modelling framework enables the assessment of resistance evolution in response to dynamic selection pressures, including changes in antimicrobial concentration and the emergence of adaptive phenotypes. Using antibiotics and antimicrobial peptides as an example, we discuss the empirical evidence and intuition behind individual model parameters. We further suggest several extensions of this framework that allow a more comprehensive and realistic prediction of bacterial escape from antimicrobials through various phenotypic and genetic mechanisms.

Virulence decomposition for bifurcating infections.

A fundamental goal in infection biology is to understand the emergence of variation in pathogen virulence-here defined as the decrease in host fitness caused by a pathogen. To uncover the sources of such variation, virulence can be decomposed into both host- and pathogen-associated components. However, decomposing virulence can be challenging owing to complex within-host pathogen dynamics such as bifurcating infections, which recently received increased empirical and theoretical attention. Bifurcating infections are characterized by the emergence of two distinct infection types: (i) terminal infections with high pathogen loads resulting in rapid host death, and (ii) persistent infections with lower loads and delayed host death. Here, we propose to use discrete mixture models to perform separate virulence decompositions for each infection type. Using this approach, we reanalysed a recently published experimental dataset on bacterial load and survival in Drosophila melanogaster. This analysis revealed several advantages of the new approach, most importantly the generation of a more comprehensive picture of the varying sources of virulence in different bacterial species. Beyond this application, our approach could provide valuable information for ground-truthing and improving theoretical models of within-host infection dynamics, which are developed to predict variation in infection outcome and pathogen virulence.

Per-pathogen virulence of HIV-1 subtypes A, C and D.

HIV-1 subtypes differ in their clinical manifestations and the speed in which they spread. In particular, the frequency of subtype C is increasing relative to subtypes A and D. We investigate whether HIV-1 subtypes A, C and D differ in their per-pathogen virulence and to what extend this explains the difference in spread between these subtypes. We use data from the hormonal contraception and HIV-1 genital shedding and disease progression among women with primary HIV infection study. For each study participant, we determine the set-point viral load value, CD4+ T cell level after primary infection and CD4+ T cell decline. Based on both the CD4+ T cell count after primary infection and CD4+ T cell decline, we estimate the time until AIDS. We then obtain our newly introduced measure of virulence as the inverse of the estimated time until AIDS. After fitting a model to the measured virulence and set-point viral load values, we tested if this relation varies per subtype. We found that subtype C has a significantly higher per-pathogen virulence than subtype A. Based on an evolutionary model, we then hypothesize that differences in the primary length of infection period cause the observed variation in the speed of spread of the subtypes.

Assessing the relative importance of bacterial resistance, persistence and hyper-mutation for antibiotic treatment failure.

To curb the rising threat of antimicrobial resistance, we need to understand the routes to antimicrobial treatment failure. Bacteria can survive treatment by using both genetic and phenotypic mechanisms to diminish the effect of antimicrobials. We assemble empirical data showing that, for example, Pseudomonas aeruginosa infections frequently contain persisters, transiently non-growing cells unaffected by antibiotics (AB) and hyper-mutators, mutants with elevated mutation rates, and thus higher probability of genetic resistance emergence. Resistance, persistence and hyper-mutation dynamics are difficult to disentangle experimentally. Hence, we use stochastic population modelling and deterministic fitness calculations to investigate the relative importance of genetic and phenotypic mechanisms for immediate treatment failure and establishment of prolonged, chronic infections. We find that persistence causes 'hidden' treatment failure with very low cell numbers if antimicrobial concentrations prevent growth of genetically resistant cells. Persister cells can regrow after treatment is discontinued and allow for resistance evolution in the absence of AB. This leads to different mutational routes during treatment and relapse of an infection. By contrast, hyper-mutation facilitates resistance evolution during treatment, but rarely contributes to treatment failure. Our findings highlight the time and concentration dependence of different bacterial mechanisms to escape AB killing, which should be considered when designing 'failure-proof' treatments.

Decomposing virulence to understand bacterial clearance in persistent infections.

Following an infection, hosts cannot always clear the pathogen, instead either dying or surviving with a persistent infection. Such variation is ecologically and evolutionarily important because it can affect infection prevalence and transmission, and virulence evolution. However, the factors causing variation in infection outcomes, and the relationship between clearance and virulence are not well understood. Here we show that sustained persistent infection and clearance are both possible outcomes across bacterial species showing a range of virulence in Drosophila melanogaster. Variation in virulence arises because of differences in the two components of virulence: bacterial infection intensity inside the host (exploitation), and the amount of damage caused per bacterium (per parasite pathogenicity). As early-phase exploitation increased, clearance rates later in the infection decreased, whereas there was no apparent effect of per parasite pathogenicity on clearance rates. Variation in infection outcomes is thereby determined by how virulence - and its components - relate to the rate of pathogen clearance. Taken together we demonstrate that the virulence decomposition framework is broadly applicable and can provide valuable insights into host-pathogen interactions.

Epidemiological and evolutionary consequences of different types of CRISPR-Cas systems.

Bacteria have adaptive immunity against viruses (phages) in the form of CRISPR-Cas immune systems. Currently, 6 types of CRISPR-Cas systems are known and the molecular study of three of these has revealed important molecular differences. It is unknown if and how these molecular differences change the outcome of phage infection and the evolutionary pressure the CRISPR-Cas systems faces. To determine the importance of these molecular differences, we model a phage outbreak entering a population defending exclusively with a type I/II or a type III CRISPR-Cas system. We show that for type III CRISPR-Cas systems, rapid phage extinction is driven by the probability to acquire at least one resistance spacer. However, for type I/II CRISPR-Cas systems, rapid phage extinction is characterized by an a threshold-like behaviour: any acquisition probability below this threshold leads to phage survival whereas any acquisition probability above it, results in phage extinction. We also show that in the absence of autoimmunity, high acquisition rates evolve. However, when CRISPR-Cas systems are prone to autoimmunity, intermediate levels of acquisition are optimal during a phage outbreak. As we predict an optimal probability of spacer acquisition 2 factors of magnitude above the one that has been measured, we discuss the origin of such a discrepancy. Finally, we show that in a biologically relevant parameter range, a type III CRISPR-Cas system can outcompete a type I/II CRISPR-Cas system with a slightly higher probability of acquisition.

Modeling Polygenic Antibiotic Resistance Evolution in Biofilms.

The recalcitrance of biofilms to antimicrobials is a multi-factorial phenomenon, including genetic, physical, and physiological changes. Individually, they often cannot account for biofilm recalcitrance. However, their combination can increase the minimal inhibitory concentration of antibiotics needed to kill bacterial cells by three orders of magnitude, explaining bacterial survival under otherwise lethal drug treatment. The relative contributions of these factors depend on the specific antibiotics, bacterial strain, as well as environmental and growth conditions. An emerging population genetic property-increased biofilm genetic diversity-further enhances biofilm recalcitrance. Here, we develop a polygenic model of biofilm recalcitrance accounting for multiple phenotypic mechanisms proposed to explain biofilm recalcitrance. The model can be used to generate predictions about the emergence of resistance-its timing and population genetic consequences. We use the model to simulate various treatments and experimental setups. Our simulations predict that the evolution of resistance is impaired in biofilms at low antimicrobial concentrations while it is facilitated at higher concentrations. In scenarios that allow bacteria exchange between planktonic and biofilm compartments, the evolution of resistance is further facilitated compared to scenarios without exchange. We compare these predictions to published experimental observations.

Applying mixture model methods to SARS-CoV-2 serosurvey data from Geneva.

Serosurveys are an important tool to estimate the true extent of the current SARS-CoV-2 pandemic. So far, most serosurvey data have been analyzed with cutoff-based methods, which dichotomize individual measurements into sero-positives or negatives based on a predefined cutoff. However, mixture model methods can gain additional information from the same serosurvey data. Such methods refrain from dichotomizing individual values and instead use the full distribution of the serological measurements from pre-pandemic and COVID-19 controls to estimate the cumulative incidence. This study presents an application of mixture model methods to SARS-CoV-2 serosurvey data from the SEROCoV-POP study from April and May 2020 in Geneva (2766 individuals). Besides estimating the total cumulative incidence in these data (8.1% (95% CI: 6.8%-9.9%)), we applied extended mixture model methods to estimate an indirect indicator of disease severity, which is the fraction of cases with a distribution of antibody levels similar to hospitalized COVID-19 patients. This fraction is 51.2% (95% CI: 15.2%-79.5%) across the full serosurvey, but differs between three age classes: 21.4% (95% CI: 0%-59.6%) for individuals between 5 and 40 years old, 60.2% (95% CI: 21.5%-100%) for individuals between 41 and 65 years old and 100% (95% CI: 20.1%-100%) for individuals between 66 and 90 years old. Additionally, we find a mismatch between the inferred negative distribution of the serosurvey and the validation data of pre-pandemic controls. Overall, this study illustrates that mixture model methods can provide additional insights from serosurvey data.

Population genetics, biofilm recalcitrance, and antibiotic resistance evolution.

Biofilms are communities of bacteria forming high-density sessile colonies. Such a lifestyle comes associated with costs and benefits: while the growth rate of biofilms is often lower than that of their free-living counterparts, this cost is readily repaid once the colony is subjected to antibiotics. Biofilms can grow in antibiotic concentrations a thousand times higher than planktonic bacteria. While numerous mechanisms have been proposed to explain biofilm recalcitrance towards antibiotics, little is yet known about their effect on the evolution of resistance. We synthesize the current understanding of biofilm recalcitrance from a pharmacodynamic and a population genetics perspective. Using the pharmacodynamic framework, we discuss the effects of various mechanisms and show that biofilms can either promote or impede resistance evolution.

Effector function does not contribute to protection from virus challenge by a highly potent HIV broadly neutralizing antibody in nonhuman primates.

Protection from immunodeficiency virus challenge in nonhuman primates (NHPs) by a first-generation HIV broadly neutralizing antibody (bnAb) b12 has previously been shown to benefit from interaction between the bnAb and Fcγ receptors (FcγRs) on immune cells. To investigate the mechanism of protection for a more potent second-generation bnAb currently in clinical trials, PGT121, we carried out a series of NHP studies. These studies included treating with PGT121 at a concentration at which only half of the animals were protected to avoid potential masking of FcγR effector function benefits by dominant neutralization and using a new variant that more completely eliminated all rhesus FcγR binding than earlier variants. In contrast to b12, which required FcγR binding for optimal protection, we concluded that PGT121-mediated protection is not augmented by FcγR interaction. Thus, for HIV-passive antibody prophylaxis, these results, together with existing literature, emphasize the importance of neutralization potency for clinical antibodies, with effector function requiring evaluation for individual antibodies.

Bacteria primed by antimicrobial peptides develop tolerance and persist.

Antimicrobial peptides (AMPs) are key components of innate immune defenses. Because of the antibiotic crisis, AMPs have also come into focus as new drugs. Here, we explore whether prior exposure to sub-lethal doses of AMPs increases bacterial survival and abets the evolution of resistance. We show that Escherichia coli primed by sub-lethal doses of AMPs develop tolerance and increase persistence by producing curli or colanic acid, responses linked to biofilm formation. We develop a population dynamic model that predicts that priming delays the clearance of infections and fuels the evolution of resistance. The effects we describe should apply to many AMPs and other drugs that target the cell surface. The optimal strategy to tackle tolerant or persistent cells requires high concentrations of AMPs and fast and long-lasting expression. Our findings also offer a new understanding of non-inherited drug resistance as an adaptive response and could lead to measures that slow the evolution of resistance.

Estimating the cumulative incidence of SARS-CoV-2 with imperfect serological tests: Exploiting cutoff-free approaches.

Large-scale serological testing in the population is essential to determine the true extent of the current SARS-CoV-2 pandemic. Serological tests measure antibody responses against pathogens and use predefined cutoff levels that dichotomize the quantitative test measures into sero-positives and negatives and use this as a proxy for past infection. With the imperfect assays that are currently available to test for past SARS-CoV-2 infection, the fraction of seropositive individuals in serosurveys is a biased estimator of the cumulative incidence and is usually corrected to account for the sensitivity and specificity. Here we use an inference method-referred to as mixture-model approach-for the estimation of the cumulative incidence that does not require to define cutoffs by integrating the quantitative test measures directly into the statistical inference procedure. We confirm that the mixture model outperforms the methods based on cutoffs, leading to less bias and error in estimates of the cumulative incidence. We illustrate how the mixture model can be used to optimize the design of serosurveys with imperfect serological tests. We also provide guidance on the number of control and case sera that are required to quantify the test's ambiguity sufficiently to enable the reliable estimation of the cumulative incidence. Lastly, we show how this approach can be used to estimate the cumulative incidence of classes of infections with an unknown distribution of quantitative test measures. This is a very promising application of the mixture-model approach that could identify the elusive fraction of asymptomatic SARS-CoV-2 infections. An R-package implementing the inference methods used in this paper is provided. Our study advocates using serological tests without cutoffs, especially if they are used to determine parameters characterizing populations rather than individuals. This approach circumvents some of the shortcomings of cutoff-based methods at exactly the low cumulative incidence levels and test accuracies that we are currently facing in SARS-CoV-2 serosurveys.

Long-term experimental evolution of HIV-1 reveals effects of environment and mutational history.

An often-returning question for not only HIV-1, but also other organisms, is how predictable evolutionary paths are. The environment, mutational history, and random processes can all impact the exact evolutionary paths, but to which extent these factors contribute to the evolutionary dynamics of a particular system is an open question. Especially in a virus like HIV-1, with a large mutation rate and large population sizes, evolution is expected to be highly predictable if the impact of environment and history is low, and evolution is not neutral. We investigated the effect of environment and mutational history by analyzing sequences from a long-term evolution experiment, in which HIV-1 was passaged on 2 different cell types in 8 independent evolutionary lines and 8 derived lines, 4 of which involved a switch of the environment. The experiments lasted for 240-300 passages, corresponding to approximately 400-600 generations or almost 3 years. The sequences show signs of extensive parallel evolution-the majority of mutations that are shared between independent lines appear in both cell types, but we also find that both environment and mutational history significantly impact the evolutionary paths. We conclude that HIV-1 evolution is robust to small changes in the environment, similar to a transmission event in the absence of an immune response or drug pressure. We also find that the fitness landscape of HIV-1 is largely smooth, although we find some evidence for both positive and negative epistatic interactions between mutations.

Number of HIV-1 founder variants is determined by the recency of the source partner infection.

During sexual transmission, the high genetic diversity of HIV-1 within an individual is frequently reduced to one founder variant that initiates infection. Understanding the drivers of this bottleneck is crucial to developing effective infection control strategies. Little is known about the importance of the source partner during this bottleneck. To test the hypothesis that the source partner affects the number of HIV founder variants, we developed a phylodynamic model calibrated using genetic and epidemiological data on all existing transmission pairs for whom the direction of transmission and the infection stage of the source partner are known. Our results suggest that acquiring infection from someone in the acute (early) stage of infection increases the risk of multiple-founder variant transmission compared with acquiring infection from someone in the chronic (later) stage of infection. This study provides the first direct test of source partner characteristics to explain the low frequency of multiple-founder strain infections.

Disentangling non-specific and specific transgenerational immune priming components in host-parasite interactions.

Exposure to a pathogen primes many organisms to respond faster or more efficiently to subsequent exposures. Such priming can be non-specific or specific, and has been found to extend across generations. Disentangling and quantifying specific and non-specific effects is essential for understanding the genetic epidemiology of a system. By combining a large infection experiment and mathematical modelling, we disentangle different transgenerational effects in the crustacean model Daphnia magna exposed to different strains of the bacterial parasite Pasteuria ramosa. In the experiment, we exposed hosts to a high dose of one of three parasite strains, and subsequently challenged their offspring with multiple doses of the same (homologous) or a different (heterologous) strain. We find that exposure of Daphnia to Pasteuria decreases the susceptibility of their offspring by approximately 50%. This transgenerational protection is not larger for homologous than for heterologous parasite challenges. Methodologically, our work represents an important contribution not only to the analysis of immune priming in ecological systems but also to the experimental assessment of vaccines. We present, for the first time, an inference framework to investigate specific and non-specific effects of immune priming on the susceptibility distribution of hosts-effects that are central to understanding immunity and the effect of vaccines.

Intestinal epithelial NAIP/NLRC4 restricts systemic dissemination of the adapted pathogen Salmonella Typhimurium due to site-specific bacterial PAMP expression.

Inflammasomes can prevent systemic dissemination of enteropathogenic bacteria. As adapted pathogens including Salmonella Typhimurium (S. Tm) have evolved evasion strategies, it has remained unclear when and where inflammasomes restrict their dissemination. Bacterial population dynamics establish that the NAIP/NLRC4 inflammasome specifically restricts S. Tm migration from the gut to draining lymph nodes. This is solely attributable to NAIP/NLRC4 within intestinal epithelial cells (IECs), while S. Tm evades restriction by phagocyte NAIP/NLRC4. NLRP3 and Caspase-11 also fail to restrict S. Tm mucosa traversal, migration to lymph nodes, and systemic pathogen growth. The ability of IECs (not phagocytes) to mount a NAIP/NLRC4 defense in vivo is explained by particularly high NAIP/NLRC4 expression in IECs and the necessity for epithelium-invading S. Tm to express the NAIP1-6 ligands-flagella and type-III-secretion-system-1. Imaging reveals both ligands to be promptly downregulated following IEC-traversal. These results highlight the importance of intestinal epithelial NAIP/NLRC4 in blocking bacterial dissemination in vivo, and explain why this constitutes a uniquely evasion-proof defense against the adapted enteropathogen S. Tm.