Effects of random mutations in the human immunodeficiency virus type 1 transcriptional promoter on viral fitness in different host cell environments

Dynamic Pneumococcal Genetic Adaptations Support Bacterial Growth and Inflammation during Coinfection with Influenza

Streptococcus pneumoniae (pneumococcus) is one of the primary bacterial pathogens that complicates influenza virus infections. These bacterial coinfections increase influenza-associated morbidity and mortality through a number of immunological and viral-mediated mechanisms, but the specific bacterial genes that contribute to postinfluenza pathogenicity are not known. Here, we used genome-wide transposon mutagenesis (Tn-Seq) to reveal bacterial genes that confer improved fitness in influenza virus-infected hosts. The majority of the 32 genes identified are involved in bacterial metabolism, including nucleotide biosynthesis, amino acid biosynthesis, protein translation, and membrane transport. We generated mutants with single-gene deletions (SGD) of five of the genes identified, SPD1414, SPD2047 (cbiO1), SPD0058 (purD), SPD1098, and SPD0822 (proB), to investigate their effects on in vivo fitness, disease severity, and host immune responses. The growth of the SGD mutants was slightly attenuated in vitro and in vivo, but each still grew to high titers in the lungs of mock- and influenza virus-infected hosts. Despite high bacterial loads, mortality was significantly reduced or delayed with all SGD mutants. Time-dependent reductions in pulmonary neutrophils, inflammatory macrophages, and select proinflammatory cytokines and chemokines were also observed. Immunohistochemical staining further revealed altered neutrophil distribution with reduced degeneration in the lungs of influenza virus-SGD mutant-coinfected animals. These studies demonstrate a critical role for specific bacterial genes and for bacterial metabolism in driving virulence and modulating immune function during influenza-associated bacterial pneumonia.

Effects of SARS-CoV-2 mutations on protein structures and intraviral protein-protein interactions

Since 2019, severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) causing coronavirus disease 2019 (COVID‐19) has infected 10 millions of people across the globe, and massive mutations in virus genome have occurred during the rapid spread of this novel coronavirus. Variance in protein sequence might lead to a change in protein structure and interaction, then further affect the viral physiological characteristics, which could bring tremendous influence on the pandemic. In this study, we investigated 20 nonsynonymous mutations in the SARS‐CoV‐2 genome in which incidence rates were all ≥ 1% as of September 1st, 2020, and then modeled and analyzed the mutant protein structures. The results showed that four types of mutations caused dramatic changes in protein structures (RMSD ≥ 5.0 Å), which were Q57H and G251V in open‐reading frames 3a (ORF3a), S194L, and R203K/G204R in nucleocapsid (N). Next, we found that these mutations also affected the binding affinity of intraviral protein interactions. In addition, the hot spots within these docking mutant complexes were altered, among which the mutation Q57H was involved in both Orf3a–S and Orf3a–Orf8 protein interactions. Besides, these mutations were widely distributed all over the world, and their occurrences fluctuated as time went on. Notably, the incidences of R203K/G204R in N and Q57H in Orf3a were both over 50% in some countries. Overall, our findings suggest that SARS‐CoV‐2 mutations could change viral protein structure, binding affinity, and hot spots of the interface, thereby might have impacts on SARS‐CoV‐2 transmission, diagnosis, and treatment of COVID‐19.

Transposon Insertion Sequencing, a Global Measure of Gene Function

The goal of genomics and systems biology is to understand how complex systems of factors assemble into pathways and structures that combine to form living organisms. Great advances in understanding biological processes result from determining the function of individual genes, a process that has classically relied on characterizing single mutations. Advances in DNA sequencing has made available the complete set of genetic instructions for an astonishing and growing number of species. To understand the function of this ever-increasing number of genes, a high-throughput method was developed that in a single experiment can measure the function of genes across the genome of an organism. This occurred approximately 10 years ago, when high-throughput DNA sequencing was combined with advances in transposon-mediated mutagenesis in a method termed transposon insertion sequencing (TIS). In the subsequent years, TIS succeeded in addressing fundamental questions regarding the genes of bacteria, many of which have been shown to play central roles in bacterial infections that result in major human diseases. The field of TIS has matured and resulted in studies of hundreds of species that include significant innovations with a number of transposons. Here, we summarize a number of TIS experiments to provide an understanding of the method and explanation of approaches that are instructive when designing a study. Importantly, we emphasize critical aspects of a TIS experiment and highlight the extension and applicability of TIS into nonbacterial species such as yeast.

Droplet Tn-Seq combines microfluidics with Tn-Seq for identifying complex single-cell phenotypes

While Tn-Seq is a powerful tool to determine genome-wide bacterial fitness in high-throughput, culturing transposon-mutant libraries in pools can mask community or other complex single-cell phenotypes. Droplet Tn-Seq (dTn-Seq) solves this problem by microfluidics facilitated encapsulation of individual transposon mutants into growth medium-in-oil droplets, thereby enabling isolated growth, free from the influence of the population. Here we describe and validate microfluidic chip design, production, encapsulation, and dTn-Seq sample preparation. We determine that 1-3% of mutants in Streptococcus pneumoniae have a different fitness when grown in isolation and show how dTn-Seq can help identify leads for gene function, including those involved in hyper-competence, processing of alpha-1-acid glycoprotein, sensitivity against the human leukocyte elastase and microcolony formation. Additionally, we show dTn-Seq compatibility with microscopy, FACS and investigations of bacterial cell-to-cell and bacteria-host cell interactions. dTn-Seq reduces costs and retains the advantages of Tn-Seq, while expanding the method's original applicability.

On the stability of sequences inserted into viral genomes

Viruses are widely used as vectors for heterologous gene expression in cultured cells or natural hosts, and therefore a large number of viruses with exogenous sequences inserted into their genomes have been engineered. Many of these engineered viruses are viable and express heterologous proteins at high levels, but the inserted sequences often prove to be unstable over time and are rapidly lost, limiting heterologous protein expression. Although virologists are aware that inserted sequences can be unstable, processes leading to insert instability are rarely considered from an evolutionary perspective. Here, we review experimental work on the stability of inserted sequences over a broad range of viruses, and we present some theoretical considerations concerning insert stability. Different virus genome organizations strongly impact insert stability, and factors such as the position of insertion can have a strong effect. In addition, we argue that insert stability not only depends on the characteristics of a particular genome, but that it will also depend on the host environment and the demography of a virus population. The interplay between all factors affecting stability is complex, which makes it challenging to develop a general model to predict the stability of genomic insertions. We highlight key questions and future directions, finding that insert stability is a surprisingly complex problem and that there is need for mechanism-based, predictive models. Combining theoretical models with experimental tests for stability under varying conditions can lead to improved engineering of viral modified genomes, which is a valuable tool for understanding genome evolution as well as for biotechnological applications, such as gene therapy.

Cat and Mouse: HIV Transcription in Latency, Immune Evasion and Cure/Remission Strategies

There is broad scientific and societal consensus that finding a cure for HIV infection must be pursued. The major barrier to achieving a cure for HIV/AIDS is the capacity of the HIV virus to avoid both immune surveillance and current antiretroviral therapy (ART) by rapidly establishing latently infected cell populations, termed latent reservoirs. Here, we provide an overview of the rapidly evolving field of HIV cure/remission research, highlighting recent progress and ongoing challenges in the understanding of HIV reservoirs, the role of HIV transcription in latency and immune evasion. We review the major approaches towards a cure that are currently being explored and further argue that small molecules that inhibit HIV transcription, and therefore uncouple HIV gene expression from signals sent by the host immune response, might be a particularly promising approach to attain a cure or remission. We emphasize that a better understanding of the game of "cat and mouse" between the host immune system and the HIV virus is a crucial knowledge gap to be filled in both cure and vaccine research.

Antibiotics Disrupt Coordination between Transcriptional and Phenotypic Stress Responses in Pathogenic Bacteria

Bacterial genes that change in expression upon environmental disturbance have commonly been seen as those that must also phenotypically matter. However, several studies suggest that differentially expressed genes are rarely phenotypically important. We demonstrate, for Gram-positive and Gram-negative bacteria, that these seemingly uncoordinated gene sets are involved in responses that can be linked through topological network analysis. However, the level of coordination is stress dependent. While a well-coordinated response is triggered in response to nutrient stress, antibiotics trigger an uncoordinated response in which transcriptionally and phenotypically important genes are neither linked spatially nor in their magnitude. Moreover, a gene expression meta-analysis reveals that genes with large fitness changes during stress have low transcriptional variation across hundreds of other conditions, and vice versa. Our work suggests that cellular responses can be understood through network models that incorporate regulatory and genetic relationships, which could aid drug target predictions and genetic network engineering.

Strain Dependent Genetic Networks for Antibiotic-Sensitivity in a Bacterial Pathogen with a Large Pan-Genome

The interaction between an antibiotic and bacterium is not merely restricted to the drug and its direct target, rather antibiotic induced stress seems to resonate through the bacterium, creating selective pressures that drive the emergence of adaptive mutations not only in the direct target, but in genes involved in many different fundamental processes as well. Surprisingly, it has been shown that adaptive mutations do not necessarily have the same effect in all species, indicating that the genetic background influences how phenotypes are manifested. However, to what extent the genetic background affects the manner in which a bacterium experiences antibiotic stress, and how this stress is processed is unclear. Here we employ the genome-wide tool Tn-Seq to construct daptomycin-sensitivity profiles for two strains of the bacterial pathogen Streptococcus pneumoniae. Remarkably, over half of the genes that are important for dealing with antibiotic-induced stress in one strain are dispensable in another. By confirming over 100 genotype-phenotype relationships, probing potassium-loss, employing genetic interaction mapping as well as temporal gene-expression experiments we reveal genome-wide conditionally important/essential genes, we discover roles for genes with unknown function, and uncover parts of the antibiotic's mode-of-action. Moreover, by mapping the underlying genomic network for two query genes we encounter little conservation in network connectivity between strains as well as profound differences in regulatory relationships. Our approach uniquely enables genome-wide fitness comparisons across strains, facilitating the discovery that antibiotic responses are complex events that can vary widely between strains, which suggests that in some cases the emergence of resistance could be strain specific and at least for species with a large pan-genome less predictable.

Distribution of mutational fitness effects and of epistasis in the 5' untranslated region of a plant RNA virus

BACKGROUND

Understanding the causes and consequences of phenotypic variability is a central topic of evolutionary biology. Mutations within non-coding cis-regulatory regions are thought to be of major effect since they affect the expression of downstream genes. To address the evolutionary potential of mutations affecting such regions in RNA viruses, we explored the fitness properties of mutations affecting the 5'-untranslated region (UTR) of a prototypical member of the picorna-like superfamily, Tobacco etch virus (TEV). This 5' UTR acts as an internal ribosomal entry site (IRES) and is essential for expression of all viral genes.

RESULTS

We determined in vitro the folding of 5' UTR using the selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) technique. Then, we created a collection of single-nucleotide substitutions on this region and evaluated the statistical properties of their fitness effects in vivo. We found that, compared to random mutations affecting coding sequences, mutations at the 5' UTR were of weaker effect. We also created double mutants by combining pairs of these single mutations and found variation in the magnitude and sign of epistatic interactions, with an enrichment of cases of positive epistasis. A correlation exists between the magnitude of fitness effects and the size of the perturbation made in the RNA folding structure, suggesting that the larger the departure from the predicted fold, the more negative impact in viral fitness.

CONCLUSIONS

Evidence that mutational fitness effects on the short 5' UTR regulatory sequence of TEV are weaker than those affecting its coding sequences have been found. Epistasis among pairs of mutations on the 5' UTR ranged between the extreme cases of synthetic lethal and compensatory. A plausible hypothesis to explain all these observations is that the interaction between the 5' UTR and the host translational machinery was shaped by natural selection to be robust to mutations, thus ensuring the homeostatic expression of viral genes even at high mutation rates.

Extensive recombination-induced disruption of genetic interactions is highly deleterious but can be partially reversed by small numbers of secondary recombination events

Although homologous recombination can potentially provide viruses with vastly more evolutionary options than are available through mutation alone, there are considerable limits on the adaptive potential of this important evolutionary process. Primary among these is the disruption of favorable coevolved genetic interactions that can occur following the transfer of foreign genetic material into a genome. Although the fitness costs of such disruptions can be severe, in some cases they can be rapidly recouped by either compensatory mutations or secondary recombination events. Here, we used a maize streak virus (MSV) experimental model to explore both the extremes of recombination-induced genetic disruption and the capacity of secondary recombination to adaptively reverse almost lethal recombination events. Starting with two naturally occurring parental viruses, we synthesized two of the most extreme conceivable MSV chimeras, each effectively carrying 182 recombination breakpoints and containing thorough reciprocal mixtures of parental polymorphisms. Although both chimeras were severely defective and apparently noninfectious, neither had individual movement-, encapsidation-, or replication-associated genome regions that were on their own "lethally recombinant." Surprisingly, mixed inoculations of the chimeras yielded symptomatic infections with viruses with secondary recombination events. These recombinants had only 2 to 6 breakpoints, had predominantly inherited the least defective of the chimeric parental genome fragments, and were obviously far more fit than their synthetic parents. It is clearly evident, therefore, that even when recombinationally disrupted virus genomes have extremely low fitness and there are no easily accessible routes to full recovery, small numbers of secondary recombination events can still yield tremendous fitness gains. Importance: Recombination between viruses can generate strains with enhanced pathological properties but also runs the risk of producing hybrid genomes with decreased fitness due to the disruption of favorable genetic interactions. Using two synthetic maize streak virus genome chimeras containing alternating genome segments derived from two natural viral strains, we examined both the fitness costs of extreme degrees of recombination (both chimeras had 182 recombination breakpoints) and the capacity of secondary recombination events to recoup these costs. After the severely defective chimeras were introduced together into a suitable host, viruses with between 1 and 3 secondary recombination events arose, which had greatly increased replication and infective capacities. This indicates that even in extreme cases where recombination-induced genetic disruptions are almost lethal, and 91 consecutive secondary recombination events would be required to reconstitute either one of the parental viruses, moderate degrees of fitness recovery can be achieved through relatively small numbers of secondary recombination events.

Variability in mutational fitness effects prevents full lethal transitions in large quasispecies populations

The distribution of mutational fitness effects (DMFE) is crucial to the evolutionary fate of quasispecies. In this article we analyze the effect of the DMFE on the dynamics of a large quasispecies by means of a phenotypic version of the classic Eigen's model that incorporates beneficial, neutral, deleterious, and lethal mutations. By parameterizing the model with available experimental data on the DMFE of Vesicular stomatitis virus (VSV) and Tobacco etch virus (TEV), we found that increasing mutation does not totally push the entire viral quasispecies towards deleterious or lethal regions of the phenotypic sequence space. The probability of finding regions in the parameter space of the general model that results in a quasispecies only composed by lethal phenotypes is extremely small at equilibrium and in transient times. The implications of our findings can be extended to other scenarios, such as lethal mutagenesis or genomically unstable cancer, where increased mutagenesis has been suggested as a potential therapy.

Transposon insertion sequencing: a new tool for systems-level analysis of microorganisms

Our knowledge of gene function has increasingly lagged behind gene discovery, hindering our understanding of the genetic basis of microbial phenotypes. Recently, however, massively parallel sequencing has been combined with traditional transposon mutagenesis in techniques referred to as transposon sequencing (Tn-seq), high-throughput insertion tracking by deep sequencing (HITS), insertion sequencing (INSeq) and transposon-directed insertion site sequencing (TraDIS), making it possible to identify putative gene functions in a high-throughput manner. Here, we describe the similarities and differences of these related techniques and discuss their application to the probing of gene function and higher-order genome organization.

Extreme genetic fragility of the HIV-1 capsid

Genetic robustness, or fragility, is defined as the ability, or lack thereof, of a biological entity to maintain function in the face of mutations. Viruses that replicate via RNA intermediates exhibit high mutation rates, and robustness should be particularly advantageous to them. The capsid (CA) domain of the HIV-1 Gag protein is under strong pressure to conserve functional roles in viral assembly, maturation, uncoating, and nuclear import. However, CA is also under strong immunological pressure to diversify. Therefore, it would be particularly advantageous for CA to evolve genetic robustness. To measure the genetic robustness of HIV-1 CA, we generated a library of single amino acid substitution mutants, encompassing almost half the residues in CA. Strikingly, we found HIV-1 CA to be the most genetically fragile protein that has been analyzed using such an approach, with 70% of mutations yielding replication-defective viruses. Although CA participates in several steps in HIV-1 replication, analysis of conditionally (temperature sensitive) and constitutively non-viable mutants revealed that the biological basis for its genetic fragility was primarily the need to coordinate the accurate and efficient assembly of mature virions. All mutations that exist in naturally occurring HIV-1 subtype B populations at a frequency >3%, and were also present in the mutant library, had fitness levels that were >40% of WT. However, a substantial fraction of mutations with high fitness did not occur in natural populations, suggesting another form of selection pressure limiting variation in vivo. Additionally, known protective CTL epitopes occurred preferentially in domains of the HIV-1 CA that were even more genetically fragile than HIV-1 CA as a whole. The extreme genetic fragility of HIV-1 CA may be one reason why cell-mediated immune responses to Gag correlate with better prognosis in HIV-1 infection, and suggests that CA is a good target for therapy and vaccination strategies.

The HIV-1 capsid protein (CA) is absolutely essential for viral replication and there is, therefore, intense evolutionary pressure for HIV-1 CA to conserve its functions. However, HIV-1 CA is also a key target of the host immune response, which should provide evolutionary pressure to diversify CA sequence. Genetic robustness, or fragility, is defined as the ability, or lack thereof, of a biological entity to preserve function in the face of sequence changes. Thus, it should be advantageous to HIV-1 CA to evolve genetic robustness. Here, we present the results of extensive, random mutagenesis of single amino acids in CA that reveal an extreme genetic fragility. Although CA participates in several steps in HIV-1 replication, the biological basis for its genetic fragility was primarily the need to participate in the efficient and proper assembly of mature virion particles. The extreme genetic fragility of HIV-1 CA may be one reason why immune responses to it correlate with better prognosis in HIV-1 infection, and suggests that CA is a good target for therapy and vaccination strategies.

Magnitude and sign epistasis among deleterious mutations in a positive-sense plant RNA virus

How epistatic interactions between mutations determine the genetic architecture of fitness is of central importance in evolution. The study of epistasis is particularly interesting for RNA viruses because of their genomic compactness, lack of genetic redundancy, and apparent low complexity. Moreover, interactions between mutations in viral genomes determine traits such as resistance to antiviral drugs, virulence and host range. In this study we generated 53 Tobacco etch potyvirus genotypes carrying pairs of single-nucleotide substitutions and measured their separated and combined deleterious fitness effects. We found that up to 38% of pairs had significant epistasis for fitness, including both positive and negative deviations from the null hypothesis of multiplicative effects. Interestingly, the sign of epistasis was correlated with viral protein-protein interactions in a model network, being predominantly positive between linked pairs of proteins and negative between unlinked ones. Furthermore, 55% of significant interactions were cases of reciprocal sign epistasis (RSE), indicating that adaptive landscapes for RNA viruses maybe highly rugged. Finally, we found that the magnitude of epistasis correlated negatively with the average effect of mutations. Overall, our results are in good agreement to those previously reported for other viruses and further consolidate the view that positive epistasis is the norm for small and compact genomes that lack genetic robustness.

A fine scale phenotype-genotype virulence map of a bacterial pathogen

A large fraction of the genes from sequenced organisms are of unknown function. This limits biological insight, and for pathogenic microorganisms hampers the development of new approaches to battle infections. There is thus a great need for novel strategies that link genotypes to phenotypes for microorganisms. We describe a high-throughput strategy based on the method Tn-seq that can be applied to any genetically manipulatable microorganism. By screening 17 in vitro and two in vivo (carriage and infection) conditions for the pathogen Streptococcus pneumoniae, we create a resource consisting of >1800 interactions that is rich in new genotype–phenotype relationships. We describe genes that are involved in differential carbon source utilization in the host, as well as genes that are involved both in virulence and in resistance against specific in vitro stresses, thereby revealing selection pressures that the pathogen experiences in vivo. We reveal the secondary response to an antibiotic, including a dual role efflux pump also involved in resistance to pH stress. Through genetic-interaction mapping and gene-expression analysis we define the mechanism of attenuation and the regulatory relationship between a two-component system and a core biosynthetic pathway specific to microorganisms. Thus, we have generated a resource that provides detailed insight into the biology and virulence of S. pneumoniae and provided a road map for similar discovery in other microorganisms.

Control of virulence by small RNAs in Streptococcus pneumoniae

Small noncoding RNAs (sRNAs) play important roles in gene regulation in both prokaryotes and eukaryotes. Thus far, no sRNA has been assigned a definitive role in virulence in the major human pathogen Streptococcus pneumoniae. Based on the potential coding capacity of intergenic regions, we hypothesized that the pneumococcus produces many sRNAs and that they would play an important role in pathogenesis. We describe the application of whole-genome transcriptional sequencing to systematically identify the sRNAs of Streptococcus pneumoniae. Using this approach, we have identified 89 putative sRNAs, 56 of which are newly identified. Furthermore, using targeted genetic approaches and Tn-seq transposon screening, we demonstrate that many of the identified sRNAs have important global and niche-specific roles in virulence. These data constitute the most comprehensive analysis of pneumococcal sRNAs and provide the first evidence of the extensive roles of sRNAs in pneumococcal pathogenesis.

Effect of host species on the distribution of mutational fitness effects for an RNA virus

Knowledge about the distribution of mutational fitness effects (DMFE) is essential for many evolutionary models. In recent years, the properties of the DMFE have been carefully described for some microorganisms. In most cases, however, this information has been obtained only for a single environment, and very few studies have explored the effect that environmental variation may have on the DMFE. Environmental effects are particularly relevant for the evolution of multi-host parasites and thus for the emergence of new pathogens. Here we characterize the DMFE for a collection of twenty single-nucleotide substitution mutants of Tobacco etch potyvirus (TEV) across a set of eight host environments. Five of these host species were naturally infected by TEV, all belonging to family Solanaceae, whereas the other three were partially susceptible hosts belonging to three other plant families. First, we found a significant virus genotype-by-host species interaction, which was sustained by differences in genetic variance for fitness and the pleiotropic effect of mutations among hosts. Second, we found that the DMFEs were markedly different between Solanaceae and non-Solanaceae hosts. Exposure of TEV genotypes to non-Solanaceae hosts led to a large reduction of mean viral fitness, while the variance remained constant and skewness increased towards the right tail. Within the Solanaceae hosts, the distribution contained an excess of deleterious mutations, whereas for the non-Solanaceae the fraction of beneficial mutations was significantly larger. All together, this result suggests that TEV may easily broaden its host range and improve fitness in new hosts, and that knowledge about the DMFE in the natural host does not allow for making predictions about its properties in an alternative host.

Quasispecies spatial models for RNA viruses with different replication modes and infection strategies

Empirical observations and theoretical studies suggest that viruses may use different replication strategies to amplify their genomes, which impact the dynamics of mutation accumulation in viral populations and therefore, their fitness and virulence. Similarly, during natural infections, viruses replicate and infect cells that are rarely in suspension but spatially organized. Surprisingly, most quasispecies models of virus replication have ignored these two phenomena. In order to study these two viral characteristics, we have developed stochastic cellular automata models that simulate two different modes of replication (geometric vs stamping machine) for quasispecies replicating and spreading on a two-dimensional space. Furthermore, we explored these two replication models considering epistatic fitness landscapes (antagonistic vs synergistic) and different scenarios for cell-to-cell spread, one with free superinfection and another with superinfection inhibition. We found that the master sequences for populations replicating geometrically and with antagonistic fitness effects vanished at low critical mutation rates. By contrast, the highest critical mutation rate was observed for populations replicating geometrically but with a synergistic fitness landscape. Our simulations also showed that for stamping machine replication and antagonistic epistasis, a combination that appears to be common among plant viruses, populations further increased their robustness by inhibiting superinfection. We have also shown that the mode of replication strongly influenced the linkage between viral loci, which rapidly reached linkage equilibrium at increasing mutations for geometric replication. We also found that the strategy that minimized the time required to spread over the whole space was the stamping machine with antagonistic epistasis among mutations. Finally, our simulations revealed that the multiplicity of infection fluctuated but generically increased along time.

Pervasive sign epistasis between conjugative plasmids and drug-resistance chromosomal mutations

Multidrug-resistant bacteria arise mostly by the accumulation of plasmids and chromosomal mutations. Typically, these resistant determinants are costly to the bacterial cell. Yet, recently, it has been found that, in Escherichia coli bacterial cells, a mutation conferring resistance to an antibiotic can be advantageous to the bacterial cell if another antibiotic-resistance mutation is already present, a phenomenon called sign epistasis. Here we study the interaction between antibiotic-resistance chromosomal mutations and conjugative (i.e., self-transmissible) plasmids and find many cases of sign epistasis (40%)—including one of reciprocal sign epistasis where the strain carrying both resistance determinants is fitter than the two strains carrying only one of the determinants. This implies that the acquisition of an additional resistance plasmid or of a resistance mutation often increases the fitness of a bacterial strain already resistant to antibiotics. We further show that there is an overall antagonistic interaction between mutations and plasmids (52%). These results further complicate expectations of resistance reversal by interdiction of antibiotic use.

Bacteria can become resistant to antibiotics by spontaneous mutation of chromosomal genes or through the acquisition of horizontally mobile genetic elements, mainly conjugative plasmids. Plasmid-borne resistance is widespread among bacterial pathogens. Plasmids generally entail a cost to the host, associated with the replication and maintenance of the genetic element and with the expression of its genes. Therefore, in the absence of antibiotic, both plasmids and resistance mutations are often deleterious and confer a fitness cost to the cell. Here we studied epistatic interactions between five natural conjugative plasmids and ten chromosomal mutations conferring resistance to three types of antibiotics, making a total of 50 different combinations of chromosomal mutations and conjugative plasmids. We show that sometimes plasmids confer an advantage to bacterial strains carrying resistance mutations in their chromosome. This occurs in 32% (16 out of 50) of tested combinations. Furthermore, in 5 out of 50 plasmid-mutations combinations studied (10%), we observed an increased fitness when a plasmid-bearing bacterial cell acquires a drug-resistant mutation. These examples of sign epistasis are highly unexpected. This work explains, at least in part, how multidrug resistance evolved so rapidly.

Fitness ranking of individual mutants drives patterns of epistatic interactions in HIV-1

Fitness interactions between mutations, referred to as epistasis, can strongly impact evolution. For RNA viruses and retroviruses with their high mutation rates, epistasis may be particularly important to overcome fitness losses due to the accumulation of deleterious mutations and thus could influence the frequency of mutants in a viral population. As human immunodeficiency virus type 1 (HIV-1) resistance to azidothymidine (AZT) requires selection of sequential mutations, it is a good system to study the impact of epistasis. Here we present a thorough analysis of a classical AZT-resistance pathway (the 41–215 cluster) of HIV-1 variants by fitness measurements in single round infection assays covering physiological drug concentrations ex vivo. The sign and value of epistasis varied and did not predict the epistatic effect on the mutant frequency. This complex behavior is explained by the fitness ranking of the variants that strongly depends on environmental factors, i.e., the presence and absence of drugs and the host cells used. Although some interactions compensate fitness losses, the observed small effect on the relative mutant frequencies suggests that epistasis might be inefficient as a buffering mechanism for fitness losses in vivo. While the use of epistasis-based hypotheses to make general assumptions on the evolutionary dynamics of viral populations is appealing, our data caution their interpretation without further knowledge on the characteristics of the viral mutant spectrum under different environmental conditions.

Varied sensitivity to therapy of HIV-1 strains in CD4+ lymphocyte sub-populations upon ART initiation

Background

Although antiretroviral therapy (ART) has proven its success against HIV-1, the long lifespan of infected cells and viral latency prevent eradication. In this study we analyzed the sensitivity to ART of HIV-1 strains in naïve, central memory and effector memory CD4⁺ lymphocyte subsets.

Methods

From five patients cellular HIV-1 infection levels were quantified before and after initiation of therapy (2-5 weeks). Through sequencing the C2V3 region of the HIV-1 gp120 envelope, we studied the effect of short-term therapy on virus variants derived from naïve, central memory and effector memory CD4⁺ lymphocyte subsets.

Results

During short-term ART, HIV-1 infection levels declined in all lymphocyte subsets but not as much as RNA levels in serum. Virus diversity in the naïve and central memory lymphocyte populations remained unchanged, whilst diversity decreased in serum and the effector memory lymphocytes. ART differentially affected the virus populations co-circulating in one individual harboring a dual HIV-1 infection. Changes in V3 charge were found in all individuals after ART initiation with increases within the effector memory subset and decreases found in the naïve cell population.

Conclusions

During early ART virus diversity is affected mainly in the serum and effector memory cell compartments. Differential alterations in V3 charge were observed between effector memory and naïve populations. While certain cell populations can be targeted preferentially during early ART, some virus strains demonstrate varied sensitivity to therapy, as shown from studying two strains within a dual HIV-1 infected individual.

Mutational fitness effects in RNA and single-stranded DNA viruses: common patterns revealed by site-directed mutagenesis studies

The fitness effects of mutations are central to evolution, yet have begun to be characterized in detail only recently. Site-directed mutagenesis is a powerful tool for achieving this goal, which is particularly suited for viruses because of their small genomes. Here, I discuss the evolutionary relevance of mutational fitness effects and critically review previous site-directed mutagenesis studies. The effects of single-nucleotide substitutions are standardized and compared for five RNA or single-stranded DNA viruses infecting bacteria, plants or animals. All viruses examined show very low tolerance to mutation when compared with cellular organisms. Moreover, for non-lethal mutations, the mean fitness reduction caused by single mutations is remarkably constant (0.10-0.13), whereas the fraction of lethals varies only modestly (0.20-0.41). Other summary statistics are provided. These generalizations about the distribution of mutational fitness effects can help us to better understand the evolution of RNA and single-stranded DNA viruses.

Simple genomes, complex interactions: epistasis in RNA virus

Owed to their reduced size and low number of proteins encoded, RNA viruses and other subviral pathogens are often considered as being genetically too simple. However, this structural simplicity also creates the necessity for viral RNA sequences to encode for more than one protein and for proteins to carry out multiple functions, all together resulting in complex patterns of genetic interactions. In this work we will first review the experimental studies revealing that the architecture of viral genomes is dominated by antagonistic interactions among loci. Second, we will also review mathematical models and provide a description of computational tools for the study of RNA virus dynamics and evolution. As an application of these tools, we will finish this review article by analyzing a stochastic bit-string model of in silico virus replication. This model analyzes the interplay between epistasis and the mode of replication on determining the population load of deleterious mutations. The model suggests that, for a given mutation rate, the deleterious mutational load is always larger when epistasis is predominantly antagonistic than when synergism is the rule. However, the magnitude of this effect is larger if replication occurs geometrically than if it proceeds linearly.

Replication mode and landscape topology differentially affect RNA virus mutational load and robustness

Regardless of genome polarity, intermediaries of complementary sense must be synthesized and used as templates for the production of new genomic strands. Depending on whether these new genomic strands become themselves templates for producing extra antigenomic ones, thus giving rise to geometric growth, or only the firstly synthesized antigenomic strands can be used to this end, thus following Luria's stamping machine model, the abundances and distributions of mutant genomes will be different. Here we propose mathematical and bit string models that allow distinguishing between stamping machine and geometric replication. We have observed that, regardless the topology of the fitness landscape, the critical mutation rate at which the master sequence disappears increases as the mechanism of replication switches from purely geometric to stamping machine. We also found that, for a wide range of mutation rates, large-effect mutations do not accumulate regardless the scheme of replication. However, mild mutations accumulate more in the geometric model. Furthermore, at high mutation rates, geometric growth leads to a population collapse for intermediate values of mutational effects at which the stamping machine still produces master genomes. We observed that the critical mutation rate was weakly dependent on the strength of antagonistic epistasis but strongly dependent on synergistic epistasis. In conclusion, we have shown that RNA viruses may increase their robustness against the accumulation of deleterious mutations by replicating as stamping machines and that the magnitude of this benefit depends on the topology of the fitness landscape assumed.

Tn-seq: high-throughput parallel sequencing for fitness and genetic interaction studies in microorganisms

Biological pathways are structured in complex networks of interacting genes. Solving the architecture of such networks may provide valuable information, such as how microorganisms cause disease. Here we present a method (Tn-seq) for accurately determining quantitative genetic interactions on a genome-wide scale in microorganisms. Tn-seq is based on the assembly of a saturated Mariner transposon insertion library. After library selection, changes in frequency of each insertion mutant are determined by sequencing of the flanking regions en masse. These changes are used to calculate each mutant’s fitness. Fitness was determined for each gene of the gram-positive bacterium Streptococcus pneumoniae, a causative agent of pneumonia and meningitis. A genome-wide screen for genetic interactions identified both alleviating and aggravating interactions that could be further divided into seven distinct categories. Due to the wide activity of the Mariner transposon, Tn-seq has the potential to contribute to the exploration of complex pathways across many different species.

Clinical significance of human immunodeficiency virus type 1 replication fitness

The relative fitness of a variant, according to population genetics theory, is that variant's relative contribution to successive generations. Most drug-resistant human immunodeficiency virus type 1 (HIV-1) variants have reduced replication fitness, but at least some of these deficits can be compensated for by the accumulation of second-site mutations. HIV-1 replication fitness also appears to influence the likelihood of a drug-resistant mutant emerging during treatment failure and is postulated to influence clinical outcomes. A variety of assays are available to measure HIV-1 replication fitness in cell culture; however, there is no agreement regarding which assays best correlate with clinical outcomes. A major limitation is that there is no high-throughput assay that incorporates an internal reference strain as a control and utilizes intact virus isolates. Some retrospective studies have demonstrated statistically significant correlations between HIV-1 replication fitness and clinical outcomes in some patient populations. However, different studies disagree as to which clinical outcomes are most closely associated with fitness. This may be in part due to assay design, sample size limitations, and differences in patient populations. In addition, the strength of the correlations between fitness and clinical outcomes is modest, suggesting that, at present, it would be difficult to utilize these assays for clinical management.

Modeling Within-Host Evolution of HIV: Mutation, Competition and Strain Replacement

Virus evolution during infection of a single individual is a well-known feature of disease progression in chronic viral diseases. However, the simplest models of virus competition for host resources show the existence of a single dominant strain that grows most rapidly during the initial period of infection and competitively excludes all other virus strains. Here, we examine the dynamics of strain replacement in a simple model that includes a convex trade-off between rapid virus reproduction and long-term host cell survival. Strains are structured according to their within-cell replication rate. Over the course of infection, we find a progression in the dominant strain from fast- to moderately-replicating virus strains featuring distinct jumps in the replication rate of the dominant strain over time. We completely analyze the model and provide estimates for the replication rate of the initial dominant strain and its successors. Our model lays the groundwork for more detailed models of HIV selection and mutation. We outline future directions and application of related models to other biological situations.

The evolution of sex: empirical insights into the roles of epistasis and drift

Despite many years of theoretical and experimental work, the explanation for why sex is so common as a reproductive strategy continues to resist understanding. Recent empirical work has addressed key questions in this field, especially regarding rates of mutation accumulation in sexual and asexual organisms, and the roles of negative epistasis and drift as sources of adaptive constraint in asexually reproducing organisms. At the same time, new ideas about the evolution of sexual recombination are being tested, including intriguing suggestions of an important interplay between sex and genetic architecture, which indicate that sex and recombination could have affected their own evolution.

Adaptation of HIV-1 depends on the host-cell environment

Many viruses have the ability to rapidly develop resistance against antiviral drugs and escape from the host immune system. To which extent the host environment affects this adaptive potential of viruses is largely unknown. Here we show that for HIV-1, the host-cell environment is key to the adaptive potential of the virus. We performed a large-scale selection experiment with two HIV-1 strains in two different T-cell lines (MT4 and C8166). Over 110 days of culture, both virus strains adapted rapidly to the MT4 T-cell line. In contrast, when cultured on the C8166 T-cell line, the same strains did not show any increase in fitness. By sequence analyses and infections with viruses expressing either yellow or cyan fluorescent protein, we were able to show that the absence of adaptation was linked to a lower recombination rate in the C8166 T-cell line. Our findings suggest that if we can manipulate the host-cellular factors that mediate viral evolution, we may be able to significantly retard viral adaptability.

Epistasis between deleterious mutations and the evolution of recombination

Epistasis and the evolution of recombination are closely intertwined: epistasis generates linkage disequilibria (i.e. statistical associations between alleles), whereas recombination breaks them up. The mutational deterministic hypothesis (MDH) states that high recombination rates are maintained because the breaking up of linkage disequilibria generated by negative epistasis enables more efficient purging of deleterious mutations. However, recent theoretical and experimental work challenges the MDH. Experimental evidence suggests that negative epistasis, required by the MDH, is relatively uncommon. On the theoretical side, population genetic models suggest that, compared with the combined effects of drift and selection, epistasis generates a negligible amount of linkage disequilibria. Here, we assess these criticisms and discuss to what extent they invalidate the MDH as an explanation for the evolution of recombination.