Estimating relative fitness in viral competition experiments

The relative fitness of viral variants has previously been defined as the slope of the logarithmic ratio of the genotype or phenotype frequencies in time plots of pairwise competition experiments. Developing mathematical models for such experiments by employing the conventional coefficient of selection s, we demonstrate that this logarithmic ratio gives the fitness difference, rather than the relative fitness. This fitness difference remains proportional to the actual replication rate realized in the particular experimental setup and hence cannot be extrapolated to other situations. Conversely, the conventional relative fitness (1 + s) should be more generic. We develop an approach to compute the generic relative fitness in conventional competition experiments. This involves an estimation of the total viral replication during the experiment and requires an estimate of the average lifetime of productively infected cells. The novel approach is illustrated by estimating the relative fitness, i.e., the relative replication rate, of a set of zidovudine-resistant human immunodeficiency virus type 1 variants. A tool for calculating the relative fitness from observed changes in viral load and genotype (or phenotype) frequencies is publically available on the website at http://www-binf.bio.uu.nl/( approximately )rdb/fitness.html.     

vFitness: a web-based computing tool for improving estimation of <Emphasis Type="Italic">in vitro</Emphasis> HIV-1 fitness experiments

Background  

The replication rate (or fitness) between viral variants has been investigated in vivo and in vitro for human immunodeficiency virus (HIV). HIV fitness plays an important role in the development and persistence of drug resistance. The accurate estimation of viral fitness relies on complicated computations based on statistical methods. This calls for tools that are easy to access and intuitive to use for various experiments of viral fitness.     

Estimating Costs and Benefits of CTL Escape Mutations in SIV/HIV Infection

Mutations that allow SIV/HIV to avoid the cytotoxic T lymphocyte (CTL) response are well documented. Recently, there have been a few attempts at estimating the costs of CTL escape mutations in terms of the reduction in viral fitness and the killing rate at which the CTL response specific to one viral epitope clears virus-infected cells. Using a mathematical model we show that estimation of both parameters depends critically on the underlying changes in the replication rate of the virus and the changes in the killing rate over time (which in previous studies were assumed to be constant). We provide a theoretical basis for estimation of these parameters using in vivo data. In particular, we show that 1) by assuming unlimited virus growth one can obtain a minimal estimate of the fitness cost of the escape mutation, and 2) by assuming no virus growth during the escape, one can obtain a minimal estimate of the average killing rate. We also discuss the conditions under which better estimates of the average killing rate can be obtained.     

Near-Neighbor Interactions in the NS3-4A Protease of HCV Impact Replicative Fitness of Drug-Resistant Viral Variants

A variety of amino acid substitutions in the NS3-4A protease of the hepatitis C virus lead to protease inhibitor (PI) resistance. Many of these significantly impair the replication fitness of the resistant variants in a genotype- and subtype-dependent manner, a critical factor in determining the probability with which resistant variants will persist. However, the underlying molecular mechanisms are unknown. Here, we present a novel residue-interaction network approach to determine how near-neighbor interactions of PI resistance mutations in NS3-4A can impact protease functional sites dependent on their genomic background. We constructed subtype-specific consensus residue networks for subtypes 1a and 1b from protease structure ensembles combined with biological properties of protein residues and evolutionary amino acid conservation. By applying local and global network topology analysis and visual exploration, we characterize PI resistance-associated sites and outline differences in near-neighbor interactions. We find local residue-interaction patterns and features at protease functional sites that are subtype specific. The noncovalent bonding patterns indicate higher fitness costs conferred by PI resistance mutations in a subtype 1b genomic background and explain the prevalence of Q80K and R155K in subtype 1a. Based on local residue interactions, we predict a subtype-specific role for the protease residue NS3–Q80 in molecular mechanisms related to the assembly of infectious virus particles that is supported by experimental data on the capacity of Q80K variants to replicate and produce infectious virus in subtype 1a and 1b cell culture.     

Changes in human immunodeficiency virus type 1 fitness and genetic diversity during disease progression

This study examined the relationship between ex vivo human immunodeficiency virus type 1 (HIV-1) fitness and viral genetic diversity during the course of HIV-1 disease. Primary HIV-1 isolates from 10 patients at different time points were competed against control HIV-1 strains in peripheral blood mononuclear cell (PBMC) cultures to determine relative fitness values. Patient HIV-1 isolates sequentially gained fitness during disease at a significant rate that directly correlated with viral load and HIV-1 env C2V3 diversity. A loss in both fitness and viral diversity was observed upon the initiation of antiretroviral therapy. A possible relationship between genotype and phenotype (virus replication efficiency) is supported by the parallel increases in ex vivo fitness and viral diversity during disease, of which the correlation is largely based on specific V3 sequences. Syncytium-inducing, CXCR4-tropic HIV-1 isolates did have higher relative fitness values than non-syncytium-inducing, CCR5-tropic HIV-1 isolates, as determined by dual virus competitions in PBMC, but increases in fitness during disease were not solely powered by a gradual switch in coreceptor usage. These data provide in vivo evidence that increasing HIV-1 replication efficiency may be related to a concomitant increase in HIV-1 diversity, which in turn may be a determining factor in disease progression.     

Modeling and estimation of replication fitness of human immunodeficiency virus type 1 in vitro experiments by using a growth competition assay

Growth competition assays have been developed to quantify the relative fitnesses of human immunodeficiency virus (HIV-1) mutants. In this article we develop mathematical models to describe viral/cellular dynamic interactions in the assay experiment, from which new competitive fitness indices or parameters are defined. These indices include the log fitness ratio (LFR), the log relative fitness (LRF), and the production rate ratio (PRR). From the population genetics perspective, we clarify the confusion and correct the inconsistency in the definition of relative fitness in the literature of HIV-1 viral fitness. The LFR and LRF are easier to estimate from the experimental data than the PRR, which was misleadingly defined as the relative fitness in recent HIV-1 research literature. Calculation and estimation methods based on two data points and multiple data points were proposed and were carefully studied. In particular, we suggest using both standard linear regression (method of least squares) and a measurement error model approach for more-accurate estimates of competitive fitness parameters from multiple data points. The developed methodologies are generally applicable to any growth competition assays. A user-friendly computational tool also has been developed and is publicly available on the World Wide Web at http://www.urmc.rochester.edu/bstools/vfitness/virusfitness.htm.     

Simian immunodeficiency virus variants that differ in pathogenicity differ in fitness under rapid cell turnover conditions

Simian immunodeficiency virus (SIV) has been shown to progress through a number of changes that lead to the emergence of pathogenic viral variants in macaques initially infected with a mildly cytopathic variant, SIVMneCL8. One of these late-stage isolates, SIVMne170, replicates to high levels in vivo and causes a rapid disease course when reintroduced into na?ve macaques, resulting in a viral set point up to 3,000-fold higher than the set point of the parental virus, SIVMneCL8. However, in cell culture both viruses replicate with similar kinetics. One major difference between in vivo and in vitro cultures is the life span of the infected cells. Here, we manipulated the life span of infected cells in vitro, and we show that the fitness of SIVMne170 in cultures with a limited cell life span dramatically increased compared to its fitness in cultures with a nonlimited life span of cells. The increase in fitness was at least partially due to the fact that the rapid turnover system eliminates the negative influence of the cytopathic effects associated with replication of SIVMne170. Because the relative fitness of SIVMneCL8 and SIVMne170 observed in the rapid turnover system more accurately reflects their fitness in vivo, the system represents an improved approach to comparing relative fitness of viruses.     

HIV-1 evolution under pressure of protease inhibitors: Climbing the stairs of viral fitness

The human immunodeficiency virus (HIV-1) has evolved into a viral quasispecies with a high replication capacity or fitness. Antiretroviral drugs potently inhibit replication of the wild-type virus, but HIV-1 responds by selection of drug-resistant variants. Here we review, in brief, the evolution of resistance to protease inhibitors that is characterized by severe fitness losses and an abundance of subsequent repair strategies. The possibility to restrict HIV-1 fitness is discussed in relation to the control of HIV-1 pathogenesis.     

A Multi-Variant,Viral Dynamic Model of Genotype 1 HCV to Assess the in vivo Evolution of Protease-Inhibitor Resistant Variants

Variants resistant to compounds specifically targeting HCV are observed in clinical trials. A multi-variant viral dynamic model was developed to quantify the evolution and in vivo fitness of variants in subjects dosed with monotherapy of an HCV protease inhibitor, telaprevir. Variant fitness was estimated using a model in which variants were selected by competition for shared limited replication space. Fitness was represented in the absence of telaprevir by different variant production rate constants and in the presence of telaprevir by additional antiviral blockage by telaprevir. Model parameters, including rate constants for viral production, clearance, and effective telaprevir concentration, were estimated from 1) plasma HCV RNA levels of subjects before, during, and after dosing, 2) post-dosing prevalence of plasma variants from subjects, and 3) sensitivity of variants to telaprevir in the HCV replicon. The model provided a good fit to plasma HCV RNA levels observed both during and after telaprevir dosing, as well as to variant prevalence observed after telaprevir dosing. After an initial sharp decline in HCV RNA levels during dosing with telaprevir, HCV RNA levels increased in some subjects. The model predicted this increase to be caused by pre-existing variants with sufficient fitness to expand once available replication space increased due to rapid clearance of wild-type (WT) virus. The average replicative fitness estimates in the absence of telaprevir ranged from 1% to 68% of WT fitness. Compared to the relative fitness method, the in vivo estimates from the viral dynamic model corresponded more closely to in vitro replicon data, as well as to qualitative behaviors observed in both on-dosing and long-term post-dosing clinical data. The modeling fitness estimates were robust in sensitivity analyses in which the restoration dynamics of replication space and assumptions of HCV mutation rates were varied.     

Genetic variation and quasi-species.

During the past year the relative fitness, that is, overall replication ability, and the fitness gain of animal virus variants have been quantified, providing new insight into the dynamics for the generation of RNA virus quasi-species. Measurements of mutant frequencies and rates of genetic diversification have confirmed the extreme complexity of RNA virus and retrovirus populations.     

Building of an experimental cline with Arabidopsis thaliana to estimate herbicide fitness cost

Various management strategies aim at maintaining pesticide resistance frequency under a threshold value by taking advantage of the benefit of the fitness penalty (the cost) expressed by the resistance allele outside the treated area or during the pesticide selection "off years." One method to estimate a fitness cost is to analyze the resistance allele frequency along transects across treated and untreated areas. On the basis of the shape of the cline, this method gives the relative contributions of both gene flow and the fitness difference between genotypes in the treated and untreated areas. Taking advantage of the properties of such migration-selection balance, an artificial cline was built up to optimize the conditions where the fitness cost of two herbicide-resistant mutants (acetolactate synthase and auxin-induced target genes) in the model species Arabidopsis thaliana could be more accurately measured. The analysis of the microevolutionary dynamics in these experimental populations indicated mean fitness costs of approximately 15 and 92% for the csr1-1 and axr2-1 resistances, respectively. In addition, negative frequency dependence for the fitness cost was also detected for the axr2-1 resistance. The advantages and disadvantages of the cline approach are discussed in regard to other methods of cost estimation. This comparison highlights the powerful ability of an experimental cline to measure low fitness costs and detect sensibility to frequency-dependent variations.     

In vivo emergence of vicriviroc resistance in a human immunodeficiency virus type 1 subtype C-infected subject

Little is known about the in vivo development of resistance to human immunodeficiency virus type 1 (HIV-1) CCR5 antagonists. We studied 29 subjects with virologic failure from a phase IIb study of the CCR5 antagonist vicriviroc (VCV) and identified one individual with HIV-1 subtype C who developed VCV resistance. Studies with chimeric envelopes demonstrated that changes within the V3 loop were sufficient to confer VCV resistance. Resistant virus showed VCV-enhanced replication, cross-resistance to another CCR5 antagonist, TAK779, and increased sensitivity to aminooxypentane-RANTES and the CCR5 monoclonal antibody HGS004. Pretreatment V3 loop sequences reemerged following VCV discontinuation, implying that VCV resistance has associated fitness costs.     

Modeling Quasispecies and Drug Resistance in Hepatitis C Patients Treated with a Protease Inhibitor

Telaprevir, a novel hepatitis C virus (HCV) NS3-4A serine protease inhibitor, has demonstrated substantial antiviral activity in patients infected with HCV. However, drug-resistant HCV variants were detected in vivo at relatively high frequency a few days after drug administration. Here we use a two-strain mathematical model to explain the rapid emergence of drug resistance in HCV patients treated with telaprevir monotherapy. We examine the effects of backward mutation and liver cell proliferation on the preexistence of the mutant virus and the competition between wild-type and drug-resistant virus during therapy. We also extend the two-strain model to a general model with multiple viral strains. Mutations during therapy only have a minor effect on the dynamics of various viral strains, although they are capable of generating low levels of HCV variants that would otherwise be completely suppressed because of fitness disadvantages. Liver cell proliferation may not affect the pretreatment frequency of mutant variants, but is able to influence the quasispecies dynamics during therapy. It is the relative fitness of each mutant strain compared with wild-type that determines which strain(s) will dominate the virus population. This study provides a theoretical framework for exploring the prevalence of preexisting mutant variants and the evolution of drug resistance during treatment with other HCV protease inhibitors or polymerase inhibitors.     

Parasitism increases and decreases the costs of insecticide resistance in mosquitoes

Adaptations conferring resistance to xenobiotics (antibiotics, insecticides, herbicides, etc.) are often costly to the organism's fitness in the absence of the selecting agent. In such conditions, and unless other mutations compensate for the costs of resistance, sensitive individuals are expected to out-reproduce resistant individuals and drive resistance alleles to a low frequency, with the rate and magnitude of this decline being proportional to the costs of resistance. However, this evolutionary dynamic is open to modification by other sources of selection acting on the relative fitness of susceptible and resistant individuals. Here we show parasitism not only as a source of selection capable of modifying the costs of organophosphate insecticide resistance in mosquitoes, but also that qualitatively different interactions (increasing or decreasing the relative fitness of resistant individuals) occurred depending on the particular form of resistance involved. As estimates of the parasite's fitness also varied according to its host's form of resistance, our data illustrate the potential for epidemiological feedbacks to influence the strength and direction of selection acting on resistance mutations in untreated environments.     

Decreased infectivity of a neutralization-resistant equine infectious anemia virus variant can be overcome by efficient cell-to-cell spread

Two variants of equine infectious anemia virus (EIAV) that differed in sensitivity to broadly neutralizing antibody were tested in direct competition assays. No differences were observed in the growth curves and relative fitness scores of EIAVs of principal neutralizing domain variants of groups 1 (EIAV(PND-1)) and 5 (EIAV(PND-5)), respectively; however, the neutralization-resistant EIAV(PND-5) variant was less infectious in single-round replication assays. Infectious center assays indicated similar rates of cell-to-cell spread, which was approximately 1,000-fold more efficient than cell-free infectivity. These data indicate that efficient cell-to-cell spread can overcome the decreased infectivity that may accompany immune escape and should be considered in studies assessing the relative levels of fitness among lentivirus variants, including HIV-1.