Within-host Fitness Research Articles

Within-Host Viral Growth and Immune Response Rates Predict Foot-and-Mouth Disease Virus Transmission Dynamics for African Buffalo.

AbstractInfectious disease dynamics operate across biological scales: pathogens replicate within hosts but transmit among populations. Functional changes in the pathogen-host interaction thus generate cascading effects across organizational scales. We investigated within-host dynamics and among-host transmission of three strains (SAT-1, -2, -3) of foot-and-mouth disease viruses (FMDVs) in their wildlife host, African buffalo. We combined data on viral dynamics and host immune responses with mathematical models to ask the following questions: How do viral and immune dynamics vary among strains? Which viral and immune parameters determine viral fitness within hosts? And how do within-host dynamics relate to virus transmission? Our data reveal contrasting within-host dynamics among viral strains, with SAT-2 eliciting more rapid and effective immune responses than SAT-1 and SAT-3. Within-host viral fitness was overwhelmingly determined by variation among hosts in immune response activation rates but not by variation among individual hosts in viral growth rate. Our analyses investigating across-scale linkages indicate that viral replication rate in the host correlates with transmission rates among buffalo and that adaptive immune activation rate determines the infectious period. These parameters define the virus's relative basic reproductive number (ℛ0), suggesting that viral invasion potential may be predictable from within-host dynamics.

Holobiont Evolution: Population Theory for the Hologenome.

AbstractThis article develops mathematical theory for the population dynamics of microbiomes with their hosts and for holobiont evolution caused by holobiont selection. The objective is to account for the formation of microbiome-host integration. Microbial population dynamic parameters must mesh with the host's for coexistence. A horizontally transmitted microbiome is a genetic system with "collective inheritance." The microbial source pool in the environment corresponds to the gamete pool for nuclear genes. Poisson sampling of the microbial source pool corresponds to binomial sampling of the gamete pool. However, holobiont selection on the microbiome does not lead to a counterpart of the Hardy-Weinberg law or to directional selection that always fixes microbial genes conferring the highest holobiont fitness. A microbe might strike an optimal fitness balance between lowering its within-host fitness while increasing holobiont fitness. Such microbes are replaced by otherwise identical microbes that contribute nothing to holobiont fitness. This replacement can be reversed by hosts that initiate immune responses to nonhelpful microbes. This discrimination leads to microbial species sorting. Host-orchestrated species sorting followed by microbial competition, rather than coevolution or multilevel selection, is predicted to be the cause of microbiome-host integration.

Variation in within-host replication kinetics among virus genotypes provides evidence of specialist and generalist infection strategies across three salmonid host species.

Theory of the evolution of pathogen specialization suggests that a specialist pathogen gains high fitness in one host, but this comes with fitness loss in other hosts. By contrast, a generalist pathogen does not achieve high fitness in any host, but gains ecological fitness by exploiting different hosts, and has higher fitness than specialists in nonspecialized hosts. As a result, specialist pathogens are predicted to have greater variation in fitness across hosts, and generalists would have lower fitness variation across hosts. We test these hypotheses by measuring pathogen replicative fitness as within-host viral loads from the onset of infection to the beginning of virus clearance, using the rhabdovirus infectious hematopoietic necrosis virus (IHNV) in salmonid fish. Based on field prevalence and virulence studies, the IHNV subgroups UP, MD, and L are specialists, causing infection and mortality in sockeye salmon, steelhead, and Chinook salmon juveniles, respectively. The UC subgroup evolved naturally from a UP ancestor and is a generalist infecting all three host species but without causing severe disease. We show that the specialist subgroups had the highest peak and mean viral loads in the hosts in which they are specialized, and they had low viral loads in nonspecialized hosts, resulting in large variation in viral load across hosts. Viral kinetics show that the mechanisms of specialization involve the ability to both maximize early virus replication and avoid clearance at later times, with different mechanisms of specialization evident in different host-virus combinations. Additional nuances in the data included different fitness levels for nonspecialist interactions, reflecting different trade-offs for specialist viruses in other hosts. The generalist UC subgroup reached intermediate viral loads in all hosts and showed the smallest variation in fitness across hosts. The evolution of the UC generalist from an ancestral UP sockeye specialist was associated with fitness increases in steelhead and Chinook salmon, but only slight decreases in fitness in sockeye salmon, consistent with low- or no-cost generalism. Our results support major elements of the specialist-generalist theory, providing evidence of a specialist-generalist continuum in a vertebrate pathogen. These results also quantify within-host replicative fitness trade-offs resulting from the natural evolution of specialist and generalist virus lineages in multi-host ecosystems.

Long-term adaptation following influenza A virus host shifts results in increased within-host viral fitness due to higher replication rates, broader dissemination within the respiratory epithelium and reduced tissue damage

The mechanisms and consequences of genome evolution on viral fitness following host shifts are poorly understood. In addition, viral fitness -the ability of an organism to reproduce and survive- is multifactorial and thus difficult to quantify. Influenza A viruses (IAVs) circulate broadly among wild birds and have jumped into and become endemic in multiple mammalian hosts, including humans, pigs, dogs, seals, and horses. H3N8 equine influenza virus (EIV) is an endemic virus of horses that originated in birds and has been circulating uninterruptedly in equine populations since the early 1960s. Here, we used EIV to quantify changes in infection phenotype associated to viral fitness due to genome-wide changes acquired during long-term adaptation. We performed experimental infections of two mammalian cell lines and equine tracheal explants using the earliest H3N8 EIV isolated (A/equine/Uruguay/63 [EIV/63]), and A/equine/Ohio/2003 (EIV/2003), a monophyletic descendant of EIV/63 isolated 40 years after the emergence of H3N8 EIV. We show that EIV/2003 exhibits increased resistance to interferon, enhanced viral replication, and a more efficient cell-to-cell spread in cells and tissues. Transcriptomics analyses revealed virus-specific responses to each virus, mainly affecting host immunity and inflammation. Image analyses of infected equine respiratory explants showed that despite replicating at higher levels and spreading over larger areas of the respiratory epithelium, EIV/2003 induced milder lesions compared to EIV/63, suggesting that adaptation led to reduced tissue pathogenicity. Our results reveal previously unknown links between virus genotype and the host response to infection, providing new insights on the relationship between virus evolution and fitness.

Fitness Estimation for Viral Variants in the Context of Cellular Coinfection.

Animal models are frequently used to characterize the within-host dynamics of emerging zoonotic viruses. More recent studies have also deep-sequenced longitudinal viral samples originating from experimental challenges to gain a better understanding of how these viruses may evolve in vivo and between transmission events. These studies have often identified nucleotide variants that can replicate more efficiently within hosts and also transmit more effectively between hosts. Quantifying the degree to which a mutation impacts viral fitness within a host can improve identification of variants that are of particular epidemiological concern and our ability to anticipate viral adaptation at the population level. While methods have been developed to quantify the fitness effects of mutations using observed changes in allele frequencies over the course of a host’s infection, none of the existing methods account for the possibility of cellular coinfection. Here, we develop mathematical models to project variant allele frequency changes in the context of cellular coinfection and, further, integrate these models with statistical inference approaches to demonstrate how variant fitness can be estimated alongside cellular multiplicity of infection. We apply our approaches to empirical longitudinally sampled H5N1 sequence data from ferrets. Our results indicate that previous studies may have significantly underestimated the within-host fitness advantage of viral variants. These findings underscore the importance of considering the process of cellular coinfection when studying within-host viral evolutionary dynamics.

Evaluating the fitness of PA/I38T-substituted influenza A viruses with reduced baloxavir susceptibility in a competitive mixtures ferret model

Baloxavir is approved in several countries for the treatment of uncomplicated influenza in otherwise-healthy and high-risk patients. Treatment-emergent viruses with reduced susceptibility to baloxavir have been detected in clinical trials, but the likelihood of widespread occurrence depends on replication capacity and onward transmission. We evaluated the fitness of A/H3N2 and A/H1N1pdm09 viruses with the polymerase acidic (PA) I38T-variant conferring reduced susceptibility to baloxavir relative to wild-type (WT) viruses, using a competitive mixture ferret model, recombinant viruses and patient-derived virus isolates. The A/H3N2 PA/I38T virus showed a reduction in within-host fitness but comparable between-host fitness to the WT virus, while the A/H1N1pdm09 PA/I38T virus had broadly similar within-host fitness but substantially lower between-host fitness. Although PA/I38T viruses replicate and transmit between ferrets, our data suggest that viruses with this amino acid substitution have lower fitness relative to WT and this relative fitness cost was greater in A/H1N1pdm09 viruses than in A/H3N2 viruses.

The within-host fitness of HIV-1 increases with age in ART-na\xefve HIV-1 subtype C infected children

As the immune system develops with age, children combat infections better. HIV-1, however, targets an activated immune system, potentially rendering children increasingly permissive to HIV-1 infection as they grow. How HIV-1 fitness changes with age in children is unknown. Here, we estimated the within-host basic reproductive ratio, R0, a marker of viral fitness, in HIV-1 subtype C-infected children in India, aged between 84 days and 17 years. We measured serial viral load and CD4 T cell counts in 171 children who initiated first-line ART. For 25 children, regular and frequent measurements provided adequate data points for analysis using a mathematical model of viral dynamics to estimate R0. For the rest, we used CD4 counts for approximate estimation of R0. The viral load decline during therapy was biphasic. The mean lifespans of productively and long-lived infected cells were 1.4 and 27.8 days, respectively. The mean R0 was 1.5 in children aged < 5 years, increased with age, and approached 6.0 at 18 years, close to 5.8 estimated previously for adults. The tolerogenic immune environment thus compromises HIV-1 fitness in young children. Early treatment initiation, when the R0 is small, will likely improve viral control, in addition to suppressing the latent reservoir.

A de novo approach to inferring within-host fitness effects during untreated HIV-1 infection.

In the absence of effective antiviral therapy, HIV-1 evolves in response to the within-host environment, of which the immune system is an important aspect. During the earliest stages of infection, this process of evolution is very rapid, driven by a small number of CTL escape mutations. As the infection progresses, immune escape variants evolve under reduced magnitudes of selection, while competition between an increasing number of polymorphic alleles (i.e., clonal interference) makes it difficult to quantify the magnitude of selection acting upon specific variant alleles. To tackle this complex problem, we developed a novel multi-locus inference method to evaluate the role of selection during the chronic stage of within-host infection. We applied this method to targeted sequence data from the p24 and gp41 regions of HIV-1 collected from 34 patients with long-term untreated HIV-1 infection. We identify a broad distribution of beneficial fitness effects during infection, with a small number of variants evolving under strong selection and very many variants evolving under weaker selection. The uniquely large number of infections analysed granted a previously unparalleled statistical power to identify loci at which selection could be inferred to act with statistical confidence. Our model makes no prior assumptions about the nature of alleles under selection, such that any synonymous or non-synonymous variant may be inferred to evolve under selection. However, the majority of variants inferred with confidence to be under selection were non-synonymous in nature, and in most cases were have previously been associated with either CTL escape in p24 or neutralising antibody escape in gp41. We also identified a putative new CTL escape site (residue 286 in gag), and a region of gp41 (including residues 644, 648, 655 in env) likely to be associated with immune escape. Sites inferred to be under selection in multiple hosts have high within-host and between-host diversity although not all sites with high between-host diversity were inferred to be under selection at the within-host level. Our identification of selection at sites associated with resistance to broadly neutralising antibodies (bNAbs) highlights the need to fully understand the role of selection in untreated individuals when designing bNAb based therapies.

Insect herbivory reshapes a native leaf microbiome.

Insect herbivory is pervasive in plant communities, but its impact on microbial plant colonizers is not well-studied in natural systems. By calibrating sequencing-based bacterial detection to absolute bacterial load, we find that the within-host abundance of most leaf microbiome (phyllosphere) taxa colonizing a native forb is amplified within leaves affected by insect herbivory. Herbivore-associated bacterial amplification reflects community-wide compositional shifts towards lower ecological diversity, but the extent and direction of such compositional shifts can be interpreted only by quantifying absolute abundance. Experimentally eliciting anti-herbivore defences reshaped within-host fitness ranks among Pseudomonas spp. field isolates and amplified a subset of putatively phytopathogenic P. syringae in a manner causally consistent with observed field-scale patterns. Herbivore damage was inversely correlated with plant reproductive success and was highly clustered across plants, which predicts tight co-clustering with putative phytopathogens across hosts. Insect herbivory may thus drive the epidemiology of plant-infecting bacteria as well as the structure of a native plant microbiome by generating variation in within-host bacterial fitness at multiple phylogenetic and spatial scales. This study emphasizes that 'non-focal' biotic interactions between hosts and other organisms in their ecological settings can be crucial drivers of the population and community dynamics of host-associated microbiomes.

Mutagenesis Scanning Uncovers Evolutionary Constraints on Tobacco Etch Potyvirus Membrane-Associated 6K2 Protein.

RNA virus high mutation rate is a double-edged sword. At the one side, most mutations jeopardize proteins functions; at the other side, mutations are needed to fuel adaptation. The relevant question then is the ratio between beneficial and deleterious mutations. To evaluate this ratio, we created a mutant library of the 6K2 gene of tobacco etch potyvirus that contains every possible single-nucleotide substitution. 6K2 protein anchors the virus replication complex to the network of endoplasmic reticulum membranes. The library was inoculated into the natural host Nicotiana tabacum, allowing competition among all these mutants and selection of those that are potentially viable. We identified 11 nonsynonymous mutations that remain in the viral population at measurable frequencies and evaluated their fitness. Some had fitness values higher than the wild-type and some were deleterious. The effect of these mutations in the structure, transmembrane properties, and function of 6K2 was evaluated in silico. In parallel, the effect of these mutations in infectivity, virus accumulation, symptoms development, and subcellular localization was evaluated in the natural host. The α-helix H1 in the N-terminal part of 6K2 turned out to be under purifying selection, while most observed mutations affect the link between transmembrane α-helices H2 and H3, fusing them into a longer helix and increasing its rigidity. In general, these changes are associated with higher within-host fitness and development of milder or no symptoms. This finding suggests that in nature selection upon 6K2 may result from a tradeoff between within-host accumulation and severity of symptoms.

Characterization of Influenza B Virus Variants with Reduced Neuraminidase Inhibitor Susceptibility.

Treatment options for influenza B virus infections are limited to neuraminidase inhibitors (NAIs), which block the neuraminidase (NA) glycoprotein on the virion surface. The development of NAI resistance would therefore result in a loss of antiviral treatment options for influenza B virus infections. This study characterized two contemporary influenza B viruses with known resistance-conferring NA amino acid substitutions, D197N and H273Y, detected during routine surveillance. The D197N and H273Y variants were characterized in vitro by assessing NA enzyme activity and affinity, as well as replication in cell culture compared to those of NAI-sensitive wild-type viruses. In vivo studies were also performed in ferrets to assess the replication and transmissibility of each variant. Mathematical models were used to analyze within-host and between-host fitness of variants relative to wild-type viruses. The data revealed that the H273Y variant had NA enzyme function similar to that of its wild type but had slightly reduced replication and transmission efficiency in vivo The D197N variant had impaired NA enzyme function, but there was no evidence of reduction in replication or transmission efficiency in ferrets. Our data suggest that the influenza B virus variant with the H273Y NA substitution had a more notable reduction in fitness compared to wild-type viruses than the influenza B variant with the D197N NA substitution. Although a D197N variant is yet to become widespread, it is the most commonly detected NAI-resistant influenza B virus in surveillance studies. Our results highlight the need to carefully monitor circulating viruses for the spread of influenza B viruses with the D197N NA substitution.

Tobamovirus 3'-Terminal Gene Overlap May be a Mechanism for within-Host Fitness Improvement.

Overlapping genes (OGs) are a universal phenomenon in all kingdoms, and viruses display a high content of OGs combined with a high rate of evolution. It is believed that the mechanism of gene overlap is based on overprinting of an existing gene. OGs help virus genes compress a maximum amount of information into short sequences, conferring viral proteins with novel features and thereby increasing their within-host fitness. Analysis of tobamovirus 3′-terminal genes reveals at least two modes of OG organization and mechanisms of interaction with the host. Originally isolated from Solanaceae species, viruses (referred to as Solanaceae-infecting) such as tobacco mosaic virus do not show 3′-terminal overlap between movement protein (MP) and coat protein (CP) genes but do contain open reading frame 6 (ORF6), which overlaps with both genes. Conversely, tobamoviruses, originally isolated from Brassicaceae species (referred to as Brassicaceae-infecting) and also able to infect Solanaceae plants, have no ORF6 but are characterized by overlapping MP and CP genes. Our analysis showed that the MP/CP overlap of Brassicaceae-infecting tobamoviruses results in the following: (i) genome compression and strengthening of subgenomic promoters; (ii) CP gene early expression directly from genomic and dicistronic MP subgenomic mRNA using an internal ribosome entry site (IRES) and a stable hairpin structure in the overlapping region; (iii) loss of ORF6, which influences the symptomatology of Solanaceae-infecting tobamoviruses; and (iv) acquisition of an IRES polypurine-rich region encoding an MP nuclear localization signal. We believe that MP/CP gene overlap may constitute a mechanism for host range expansion and virus adjustment to Brassicaceae plants.

Fast-killing parasites can be favoured in spatially structured populations.

It is becoming increasingly clear that the evolution of infectious disease is influenced by host population structure. Theory predicts that parasites should be more ‘prudent’—less transmissible—in spatially structured host populations. However, here we (i) highlight how low transmission, the phenotype being selected for in this in context, may also be achieved by rapacious host exploitation, if fast host exploitation confers a local, within-host competitive advantage and (ii) test this novel concept in a bacteria–virus system. We found that limited host availability and, to a lesser extent, low relatedness favour faster-killing parasites with reduced transmission. By contrast, high host availability and high relatedness favour slower-killing, more transmissible parasites. Our results suggest high, rather than low, virulence may be selected in spatially structured host–parasite communities where local competition and hence selection for a within-host fitness advantage is high.This article is part of the themed issue ‘Opening the black box: re-examining the ecology and evolution of parasite transmission’.

Feline Immunodeficiency Virus Cross-Species Transmission: Implications for Emergence of New Lentiviral Infections.

Owing to a complex history of host-parasite coevolution, lentiviruses exhibit a high degree of species specificity. Given the well-documented viral archeology of human immunodeficiency virus (HIV) emergence following human exposures to simian immunodeficiency virus (SIV), an understanding of processes that promote successful cross-species lentiviral transmissions is highly relevant. We previously reported natural cross-species transmission of a subtype of feline immunodeficiency virus, puma lentivirus A (PLVA), between bobcats (Lynx rufus) and mountain lions (Puma concolor) for a small number of animals in California and Florida. In this study, we investigate host-specific selection pressures, within-host viral fitness, and inter- versus intraspecies transmission patterns among a larger collection of PLV isolates from free-ranging bobcats and mountain lions. Analyses of proviral and viral RNA levels demonstrate that PLVA fitness is severely restricted in mountain lions compared to that in bobcats. We document evidence of diversifying selection in three of six PLVA genomes from mountain lions, but we did not detect selection among 20 PLVA isolates from bobcats. These findings support the hypothesis that PLVA is a bobcat-adapted virus which is less fit in mountain lions and under intense selection pressure in the novel host. Ancestral reconstruction of transmission events reveals that intraspecific PLVA transmission has occurred among panthers (Puma concolor coryi) in Florida following the initial cross-species infection from bobcats. In contrast, interspecific transmission from bobcats to mountain lions predominates in California. These findings document outcomes of cross-species lentiviral transmission events among felids that compare to the emergence of HIV from nonhuman primates.IMPORTANCE Cross-species transmission episodes can be singular, dead-end events or can result in viral replication and spread in the new species. The factors that determine which outcome will occur are complex, and the risk of new virus emergence is therefore difficult to predict. We used molecular techniques to evaluate the transmission, fitness, and adaptation of puma lentivirus A (PLVA) between bobcats and mountain lions in two geographic regions. Our findings illustrate that mountain lion exposure to PLVA is relatively common but does not routinely result in communicable infections in the new host. This is attributed to efficient species barriers that largely prevent lentiviral adaptation. However, the evolutionary capacity for lentiviruses to adapt to novel environments may ultimately overcome host restriction mechanisms over time and under certain ecological circumstances. This phenomenon provides a unique opportunity to examine cross-species transmission events leading to new lentiviral emergence.

Bacterial evolution: Making a host-adapted bacterium.

Genomic analysis of Salmonella Enteritidis isolates taken from a single individual with a chronic and relapsing infection reveals how these bacteria have adapted to their host surroundings. Increased within-host fitness comes at the expense of ancillary traits such as virulence.

Multiple Barriers to the Evolution of Alternative Gene Orders in a Positive-Strand RNA Virus.

The order in which genes are organized within a genome is generally not conserved between distantly related species. However, within virus orders and families, strong conservation of gene order is observed. The factors that constrain or promote gene-order diversity are largely unknown, although the regulation of gene expression is one important constraint for viruses. Here we investigate why gene order is conserved for a positive-strand RNA virus encoding a single polyprotein in the context of its authentic multicellular host. Initially, we identified the most plausible trajectory by which alternative gene orders could evolve. Subsequently, we studied the accessibility of key steps along this evolutionary trajectory by constructing two virus intermediates: (1) duplication of a gene followed by (2) loss of the ancestral gene. We identified five barriers to the evolution of alternative gene orders. First, the number of viable positions for reordering is limited. Second, the within-host fitness of viruses with gene duplications is low compared to the wild-type virus. Third, after duplication, the ancestral gene copy is always maintained and never the duplicated one. Fourth, viruses with an alternative gene order have even lower fitness than viruses with gene duplications. Fifth, after more than half a year of evolution in isolation, viruses with an alternative gene order are still vastly inferior to the wild-type virus. Our results show that all steps along plausible evolutionary trajectories to alternative gene orders are highly unlikely. Hence, the inaccessibility of these trajectories probably contributes to the conservation of gene order in present-day viruses.

Within-Host Niche Differences and Fitness Trade-offs Promote Coexistence of Plant Viruses.

Pathogens live in diverse, competitive communities, yet the processes that maintain pathogen diversity remain elusive. Here, we use a species-rich, well-studied plant virus system, the barley yellow dwarf viruses, to examine the mechanisms that regulate pathogen diversity. We empirically parameterized models of three viruses, their two aphid vectors, and one perennial grass host. We found that high densities of both aphids maximized virus diversity and that competition limited the coexistence of two closely related viruses. Even limited ability to simultaneously infect (coinfect) host individuals strongly promoted virus coexistence; preventing coinfection led to priority effects. Coinfection generated stabilizing niche differences by allowing viruses to share hosts. However, coexistence also required trade-offs between vector generalist and specialist life-history strategies. Our predicted outcomes broadly concur with previous field observations. These results show how competition within individual hosts and vectors may lead to unexpected population-level outcomes between pathogens, including coexistence, competitive exclusion, and priority effects, and how contemporary coexistence theory can help to predict these outcomes.

Maintenance of phenotypic diversity within a set of virulence encoding genes of the malaria parasite Plasmodium falciparum.

Infection by the human malaria parasite Plasmodium falciparum results in a broad spectrum of clinical outcomes, ranging from severe and potentially life-threatening malaria to asymptomatic carriage. In a process of naturally acquired immunity, individuals living in malaria-endemic regions build up a level of clinical protection, which attenuates infection severity in an exposure-dependent manner. Underlying this shift in the immunoepidemiology as well as the observed range in malaria pathogenesis is the var multigene family and the phenotypic diversity embedded within. The var gene-encoded surface proteins Plasmodium falciparum erythrocyte membrane protein 1 mediate variant-specific binding of infected red blood cells to a diverse set of host receptors that has been linked to specific disease manifestations, including cerebral and pregnancy-associated malaria. Here, we show that cross-reactive immune responses, which minimize the within-host benefit of each additionally expressed gene during infection, can cause selection for maximum phenotypic diversity at the genome level. We further show that differential functional constraints on protein diversification stably maintain uneven ratios between phenotypic groups, in line with empirical observation. Our results thus suggest that the maintenance of phenotypic diversity within P. falciparum is driven by an evolutionary trade-off that optimizes between within-host parasite fitness and between-host selection pressure.

The impact of high-order epistasis in the within-host fitness of a positive-sense plant RNA virus.

RNA viruses are the main source of emerging infectious diseases because of the evolutionary potential bestowed by their fast replication, large population sizes and high mutation and recombination rates. However, an equally important property, which is usually neglected, is the topography of the fitness landscape. How many fitness maxima exist and how well they are connected is especially interesting, as this determines the number of accessible evolutionary pathways. To address this question, we have reconstructed a region of the fitness landscape of tobacco etch potyvirus constituted by mutations observed during the experimental adaptation of the virus to the novel host Arabidopsis thaliana. Fitness was measured for many genotypes and showed the existence of multiple peaks and holes in the landscape. We found prevailing epistatic effects between mutations, with cases of reciprocal sign epistasis being common among pairs of mutations. We also found that high-order epistasis was as important as pairwise epistasis in their contribution to fitness. Therefore, results suggest that the landscape was rugged due to the existence of holes caused by lethal genotypes, that a very limited number of potential neutral paths exist and that it contained a single adaptive peak.

Evolutionary Trade-Offs Underlie the Multi-faceted Virulence of Staphylococcus aureus.

Bacterial virulence is a multifaceted trait where the interactions between pathogen and host factors affect the severity and outcome of the infection. Toxin secretion is central to the biology of many bacterial pathogens and is widely accepted as playing a crucial role in disease pathology. To understand the relationship between toxicity and bacterial virulence in greater depth, we studied two sequenced collections of the major human pathogen Staphylococcus aureus and found an unexpected inverse correlation between bacterial toxicity and disease severity. By applying a functional genomics approach, we identified several novel toxicity-affecting loci responsible for the wide range in toxic phenotypes observed within these collections. To understand the apparent higher propensity of low toxicity isolates to cause bacteraemia, we performed several functional assays, and our findings suggest that within-host fitness differences between high- and low-toxicity isolates in human serum is a contributing factor. As invasive infections, such as bacteraemia, limit the opportunities for onward transmission, highly toxic strains could gain an additional between-host fitness advantage, potentially contributing to the maintenance of toxicity at the population level. Our results clearly demonstrate how evolutionary trade-offs between toxicity, relative fitness, and transmissibility are critical for understanding the multifaceted nature of bacterial virulence.