A multiple-site-specific heteroduplex tracking assay as a tool for the study of viral population dynamics

Short Communication: A Recombinant Variant with Increased Envelope Entry Efficiency Emerged During Early Infection of an HIV-1 Subtype C Dual Infected Rapid Progressor

Mutations in functionally constrained sites of the HIV envelope (Env) can affect entry efficiency and are potential targets for vaccine and drug design. We investigated Du151, a dual-infected individual with rapid disease progression. At her death 19 months postinfection (mpi), she was infected with a recombinant variant, which outgrew both parental viruses. We aimed to determine whether the recombinant virus had enhanced Env entry efficiency compared to the parental viruses and to identify the functional determinant. We generated 15 env clones at 1, 2, 8, and 19 mpi. Pseudovirus carrying a recombinant Env clone (PSV clone), C18 (19 mpi), had significantly higher entry efficiency compared to the parents, suggesting that the recombinant virus had enhanced fitness. To identify the functional determinant, we compared two recombinant PSV clones (C18 and C63)-differing in entry efficiency (2-fold) and by four and three amino acids in gp120 and gp41, respectively. The increased entry efficiency of a C18-gp41 PSV chimera indicated that the three amino acids in the C18 gp41 region were involved (K658, G671, and F717). Site-directed mutagenesis of the three amino acids of C63 showed that a single amino acid mutation, R658K, increased pseudovirion entry efficiency. The introduction of R658 into two PSV clones (C1 and C18) decreased their entry efficiency, suggesting that R658 carries a fitness cost. Thus, our data suggest that a recombinant virus emerged at 19 mpi with enhanced Env entry efficiency. Therefore, K658 in gp41 could in part be a contributing factor to the increased viral load and rapid disease progression of Du151.

Distribution, persistence and interchange of Epstein-Barr virus strains among PBMC, plasma and saliva of primary infection subjects

Our study aimed at investigating the distribution, persistence and interchange of viral strains among peripheral blood mononuclear cells (PBMC), plasma and saliva of primary Epstein-Barr virus (EBV) infection subjects. Twelve infectious mononucleosis (IM) patients and eight asymptomatic individuals (AS) with primary EBV infection were followed longitudinally at several time points for one year from the time of diagnosis, when blood and saliva samples were collected and separated into PBMC, plasma and saliva, representing circulating B cell, plasma and epithelial cell compartments, respectively. To survey the viral strains, genotyping assays for the natural polymorphisms in two latent EBV genes, EBNA2 and LMP1, were performed and consisted of real-time PCR on EBNA2 to distinguish type 1 and 2 viruses, fluorescent-based 30-bp typing assay on LMP1 to distinguish deletion and wild type LMP1, and fluorescent-based heteroduplex tracking assays on both EBNA2 and LMP1 to distinguish defined polymorphic variants. No discernible differences were observed between IM patients and AS. Multiple viral strains were acquired early at the start of infection. Stable persistence of dominant EBV strains in the same tissue compartment was observed throughout the longitudinal samples. LMP1-defined strains, China 1, China 2 and Mediterranean+, were the most common strains observed. EBNA2-defined groups 1 and 3e predominated the PBMC and saliva compartments. Concordance of EBNA2 and LMP1 strains between PBMC and saliva suggested ready interchange of viruses between circulating B cell and epithelial cell pools, whilst discordance of viral strains observed between plasma and PBMC/saliva indicated presence of viral pools in other undetermined tissue compartments. Taken together, the results indicated that the distribution, persistence and interchange of viral strains among the tissue compartments are more complex than those proposed by the current model of EBV life cycle.

Inferring viral population structures using heteroduplex mobility and DNA sequence analyses

Heteroduplex mobility (HMA) and tracking assays (HTA) are used to assess genetic relationships between DNA molecules. While distinguishing relationships between clonal or nearly clonal molecules is relatively straightforward, inferring population structures is more complex. To address this issue, HIV-1 quasispecies with varying levels of diversity were studied using both HTA and DNA sequencing. Viral diversity estimates and the temporal features of virus evolution were found to be generally concordant between HTA and DNA sequencing. In addition, the distribution of pairwise differences and the rates of virus divergence were similar between the two methods. These findings support the use of HTA to characterize variant populations of DNA and strengthen previous inferences concerning the evolution of HIV-1 over the course of infection.

Minority HIV-1 resistant variants in recent infection and in patients who failed first-line antiretroviral therapy with no detectable resistance-associated mutations in Thailand

Through the Thai National AIDS Program, approximately 200,000 patients infected with HIV are on antiretroviral (ARV) therapy. Although studies have shown low prevalence of primary HIV-1 resistance transmission in Thailand and in Southeast Asia where subtype CRF01_AE is predominant, minority HIV-1 drug resistance has not been studied. Two groups of patients, whose conventional genotyping results showed no drug resistance-associated mutations, were investigated: 104 homosexual men recently infected with HIV-1, naïve to ARV treatment and 22 first-line non-nucleoside reverse transcriptase inhibitor (NNRTI)-based failure patients. Pyrosequencing (PSQ) assay was developed to detect and quantify minority Y181C and M184V variants from the patients' plasma samples. The sensitivity of PSQ to detect minority Y181C and M184V variants was approximately 1%. 1/104 (0.5%) and 3/101 (3%) samples were found harboring Y181C and M184V in the group of homosexual men recently infected with HIV-1. In patients with first-line treatment failure, one had a minority M184V mutation (4.5%). The prevalence of Y181C and M184V minority variants in homosexual men recently infected and naïve to treatment was low in Thailand. Systematic monitoring of primary resistance transmission in Thailand and this region is essential to guide whether genotypic resistance test is required prior to commencing the first-line highly active antiretroviral therapy (HAART).

CAPRISA 004 tenofovir microbicide trial: no impact of tenofovir gel on the HIV transmission bottleneck

Alterations of the genital mucosal barrier may influence the number of viruses transmitted from a human immunodeficiency virus-infected source host to the newly infected individual. We used heteroduplex tracking assay and single-genome sequencing to investigate the effect of a tenofovir-based microbicide gel on the transmission bottleneck in women who seroconverted during the CAPRISA 004 microbicide trial. Seventy-seven percent (17 of 22; 95% confidence interval [CI], 56%-90%) of women in the tenofovir gel arm were infected with a single virus compared with 92% (13 of 14; 95% CI, 67%->99%) in the placebo arm (P = .37). Tenofovir gel had no discernable impact on the transmission bottleneck.

Interplay between single resistance-associated mutations in the HIV-1 protease and viral infectivity, protease activity, and inhibitor sensitivity

Resistance-associated mutations in the HIV-1 protease modify viral fitness through changes in the catalytic activity and altered binding affinity for substrates and inhibitors. In this report, we examine the effects of 31 mutations at 26 amino acid positions in protease to determine their impact on infectivity and protease inhibitor sensitivity. We found that primary resistance mutations individually decrease fitness and generally increase sensitivity to protease inhibitors, indicating that reduced virion-associated protease activity reduces virion infectivity and the reduced level of per virion protease activity is then more easily titrated by a protease inhibitor. Conversely, mutations at more variable positions (compensatory mutations) confer low-level decreases in sensitivity to all protease inhibitors with little effect on infectivity. We found significant differences in the observed effect on infectivity with a pseudotype virus assay that requires the protease to cleave the cytoplasmic tail of the amphotropic murine leukemia virus (MuLV) Env protein. Additionally, we were able to mimic the fitness loss associated with resistance mutations by directly reducing the level of virion-associated protease activity. Virions containing 50% of a D25A mutant protease were 3- to 5-fold more sensitive to protease inhibitors. This level of reduction in protease activity also resulted in a 2-fold increase in sensitivity to nonnucleoside inhibitors of reverse transcriptase and a similar increase in sensitivity to zidovudine (AZT), indicating a pleiotropic effect associated with reduced protease activity. These results highlight the interplay between enzyme activity, viral fitness, and inhibitor mechanism and sensitivity in the closed system of the viral replication complex.

Minority variants of drug-resistant HIV

Minor drug-resistant variants exist in every patient infected with human immunodeficiency virus (HIV). Because these minority variants are usually present at very low levels, they cannot be detected and quantified using conventional genotypic and phenotypic tests. Recently, several assays have been developed to characterize these low-abundance drug-resistant variants in the large, genetically complex population that is present in every HIV-infected individual. The most important issue is what results generated by these assays can predict clinical or treatment outcomes and might guide the management of patients in clinical practice. Cutoff values for the detection of these low-abundance viral variants that predict an increased risk of treatment failure should be determined. These thresholds may be specific for each mutation and treatment regimen. In this review, we summarize the attributes and limitations of the currently available detection assays and review the existing information about both acquired and transmitted drug-resistant minority variants.

Misclassification of drug failure in Plasmodium falciparum clinical trials in southeast Asia

Most trials of antimalarials occur in areas in which reinfections are possible. For Plasmodium falciparum, reinfections are distinguished from recrudescences by polymerase chain reaction analysis of 3 polymorphic genes. However, the validity of this approach has never been rigorously tested. We tested for misclassification in 6 patients from clinical trials in Thailand and Cambodia who were classified as being reinfected by the standard polymerase chain reaction protocol. Using heteroduplex tracking assays and direct DNA sequencing, we found that 5 (83%) of 6 patients were misclassified. Misclassification in this manner overestimates the efficacy of antimalarials and delays the recognition of decreasing therapeutic efficacy, thus delaying potential changes in policy.

Compartmentalized human immunodeficiency virus type 1 originates from long-lived cells in some subjects with HIV-1-associated dementia

Human immunodeficiency virus type 1 (HIV-1) invades the central nervous system (CNS) shortly after systemic infection and can result in the subsequent development of HIV-1–associated dementia (HAD) in a subset of infected individuals. Genetically compartmentalized virus in the CNS is associated with HAD, suggesting autonomous viral replication as a factor in the disease process. We examined the source of compartmentalized HIV-1 in the CNS of subjects with HIV-1–associated neurological disease and in asymptomatic subjects who were initiating antiretroviral therapy. The heteroduplex tracking assay (HTA), targeting the variable regions of env, was used to determine which HIV-1 genetic variants in the cerebrospinal fluid (CSF) were compartmentalized and which variants were shared with the blood plasma. We then measured the viral decay kinetics of individual variants after the initiation of antiretroviral therapy. Compartmentalized HIV-1 variants in the CSF of asymptomatic subjects decayed rapidly after the initiation of antiretroviral therapy, with a mean half-life of 1.57 days. Rapid viral decay was also measured for CSF-compartmentalized variants in four HAD subjects (t_1/2 mean = 2.27 days). However, slow viral decay was measured for CSF-compartmentalized variants from an additional four subjects with neurological disease (t_1/2 range = 9.85 days to no initial decay). The slow decay detected for CSF-compartmentalized variants was not associated with poor CNS drug penetration, drug resistant virus in the CSF, or the presence of X4 virus genotypes. We found that the slow decay measured for CSF-compartmentalized variants in subjects with neurological disease was correlated with low peripheral CD4 cell count and reduced CSF pleocytosis. We propose a model in which infiltrating macrophages replace CD4⁺ T cells as the primary source of productive viral replication in the CNS to maintain high viral loads in the CSF in a substantial subset of subjects with HAD.

Infection of the central nervous system (CNS) with human immunodeficiency virus type 1 (HIV-1) can lead to the development of HIV-1–associated dementia, a severe neurological disease that results in cognitive and motor impairment. Individuals that are chronically infected with HIV-1 sometimes display unique viral variants in their cerebrospinal fluid (CSF) that are not detected in the blood virus population, termed CSF-compartmentalized variants. The cell type that produces CSF-compartmentalized virus throughout the course of infection has not been determined. We used a sensitive assay to detect compartmentalized variants in the CSF of subjects with and without neurological disease, and then measured the decay kinetics of compartmentalized virus when subjects were starting antiretroviral therapy. We found that compartmentalized virus decays rapidly in asymptomatic subjects. Additionally, we detected differential decay (i.e. rapid or slow) in subjects with neurological disease, and this was associated with the number of white blood cells in the CSF. Our data supports a model of HIV-1 infection in the CNS where compartmentalized virus is produced by a long-lived cell type (slow decay), and this virus can be amplified by short-lived cells (rapid decay) that traffic into the CNS, but is increasingly produced from long-lived cells in the immunodeficient state.

Major coexisting human immunodeficiency virus type 1 env gene subpopulations in the peripheral blood are produced by cells with similar turnover rates and show little evidence of genetic compartmentalization

A distinctive feature of chronic human immunodeficiency virus type 1 (HIV-1) infection is the presence of multiple coexisting genetic variants, or subpopulations, that comprise the HIV-1 population detected in the peripheral blood. Analysis of HIV-1 RNA decay dynamics during the initiation of highly active antiretroviral therapy (HAART) has been a valuable tool for modeling the life span of infected cells that produce the bulk HIV-1 population. However, different HIV-1 target cells may have different turnover rates, and it is not clear whether the bulk HIV-1 RNA decay rate actually represents a composite of the decay rates of viral subpopulations compartmentalized in different cellular subsets with different life spans. Using heteroduplex tracking assays targeting the highly variable V3 or V4-V5 regions of the HIV-1 env gene in eight subjects, we found that all detectable coexisting HIV-1 variants in the peripheral blood generally decayed at similar rates during the initiation of HAART, suggesting that all of the variants were produced by cells with similar life spans. Furthermore, single genome amplification and coreceptor phenotyping revealed that in two subjects coexisting HIV-1 variants with distinct CXCR4 or CCR5 coreceptor phenotypes decayed with similar rates. Also, in nine additional subjects, recombination and a lack of genetic compartmentalization between X4 and R5 variants were observed, suggesting an overlap in host cell range. Our results suggest that the HIV-1 env subpopulations detectable in the peripheral blood are produced by cells with similar life spans and are not genetically isolated within particular cell types.

Identification and recovery of minor HIV-1 variants using the heteroduplex tracking assay and biotinylated probes

We describe a method to identify and recover minor human immunodeficiency virus type 1 (HIV-1) sequence variants from a complex population. The original heteroduplex tracking assay (HTA) was modified by incorporating a biotin tag into the probe to allow for direct sequence determination of the query strand. We used this approach to recover sequences from minor HIV-1 variants in the V3 region of the env gene, and to identify minor drug-resistant variants in pro. The biotin-HTA targeting of the V3 region of env allowed us to detect minor V3 variants, of which 45% were classified as CXCR4-using viruses. In addition, the biotin-protease HTA was able to detect mixtures of wild-type sequence and drug-resistance mutations in four subjects that were not detected by bulk sequence analysis. The biotin-HTA is a robust assay that first separates genetic variants then allows direct sequence analysis of major and minor variants.

Multiplex mutagenically separated polymerase chain reaction assay for rapid detection of human mitochondrial DNA variations in coding area

Recently, it has been recognized that information in the mitochondrial DNA (mtDNA) coding region can provide additional forensic discrimination with respect to the standard typing of the D-loop region, increasing the forensic power of mtDNA testing, which is sometimes rather limited. In the present study, we simultaneously typed ten single nucleotide polymorphisms (SNP) in the coding region by use of mutagenically separated polymerase chain reaction (MS-PCR) in the Chinese Chengdu population. This technique, in which different-size allele-specific primers were used, specifically amplified both alleles of mtDNA in the same reaction. Subsequent gel electrophoresis showed ten of the allelic products of different loci. Using multiplex MS-PCR, 30 primers were added simultaneously into one reaction tube to identify ten SNPs. The mtDNA variations of 160 individuals from the Chinese Chengdu population were examined and classified into 18 haplotypes. The multiplex MS-PCR method is suitable for large-scale screening studies of mtDNA variability because it is both rapid and economical.

Characterization of minority subpopulations in the mutant spectrum of HIV-1 quasispecies by successive specific amplifications

RNA viruses do not replicate as defined genomic nucleotide sequences but rather as complex distributions of mutant genomes termed viral quasispecies. Quasispecies dynamics has a number of relevant biological consequences in ribo- and retroviruses, among these the possible presence of memory genomes as minority components of their mutant spectra. Minority memory genomes reflect those viral subpopulations that were dominant at an earlier phase of viral evolution, and can quickly re-emerge to react to certain selective pressures, as it was documented with HIV-1 in vivo. Therefore, an adequate clinical management of HIV-1 requires the development of experimental methods for the detection and quantification of minority viral subpopulations, even at levels of less than 1% of the total quasispecies. We describe a new approach based on successive, highly specific PCR amplifications, which allows the genetic characterization of minority genomes present in increasingly smaller proportion in viral populations. We have coined the term 'quasispecies diving' to reflect the progressive draw on minority or 'deeper' genomes in the mutant spectrum of the quasispecies. In the case of the multidrug-resistant HIV-1 strain analyzed here, quasispecies diving allowed the detection of mutant minority genomes at an unprecedented level of 0.0054% of the amplified viral population. This approach represents a general strategy for the genetic characterization of smaller minority genomes in complex molecular populations.

Minority-variant pfcrt K76T mutations and chloroquine resistance, Malawi

Genotyping of the chloroquine-resistance biomarker pfcrt (Plasmodium falciparum chloroquine resistance transporter gene) suggests that, in the absence of chloroquine pressure, Plasmodium falciparum parasites in Malawi have reverted to chloroquine sensitivity. However, malaria infections in Africa are commonly polyclonal, and standard PCRs cannot detect minority genotypes if present in <20% of the parasites in an individual host. We have developed a multiple site-specific heteroduplex tracking assay (MSS-HTA) that can detect pfcrt 76T mutant parasites consisting of as little as 1% of the parasite population. In clinical samples, no pfcrt 76T was detected in 87 pregnant Malawian women by standard PCR. However, 22 (25%) contained minority-variant resistant genotypes detected by the MSS-HTA. These results were confirmed by subcloning and sequencing. This finding suggests that the chloroquine-resistant genotype remains common in Malawians and that PCR-undetectable drug-resistant genotypes may be present in disease-endemic populations. Surveillance for minority-variant drug resistant mutations may be useful in making antimalarial drug policy.

Independent evolution of human immunodeficiency virus type 1 env V1/V2 and V4/V5 hypervariable regions during chronic infection

Using DNA heteroduplex tracking assays, we characterized human immunodeficiency virus type 1 env V4/V5 genetic populations in multiple blood plasma samples collected over an average of 7 months from 24 chronically infected human subjects. We observed complex and dynamic V4/V5 genetic populations in most subjects. Comparisons of V4/V5 and V1/V2 population changes over the course of the study showed that major shifts in genetic populations frequently occurred in one region but not the other, and these observations were independently confirmed in one subject by single-genome sequencing. These results suggest that the V1/V2 and V4/V5 regions of env often evolve independently during chronic infection.

Sequential turnover of human immunodeficiency virus type 1 env throughout the course of infection

We examined the rates of variant population turnover of the V1-V2 and V4-V5 hypervariable domains of the human immunodeficiency virus type 1 (HIV-1) gp120 molecule in longitudinal plasma samples from 14 men with chronic HIV-1 infection using heteroduplex tracking assays (HTA). Six men had high rates of CD4+ T-cell loss, and eight men had low rates of CD4+ T-cell loss over 2.5 to 8 years of infection. We found that V1-V2 and V4-V5 env populations changed dramatically over time in all 14 subjects; the changes in these regions were significantly correlated with each another over time. The subjects with rapid CD4 loss had significantly less change in their env populations than the subjects with slow CD4 loss. The two subjects with rapid CD4 loss and sustained low CD4 counts (<150/microl for at least 2 years) showed stabilization of their V1-V2 and V4-V5 populations as reflected by low levels of total change in HTA pattern and low HTA indices (a novel measure of the emergence of new bands and band distribution); this stabilization was not observed in other subjects. The stabilization of env variant populations at low CD4 counts following periods of rapid viral evolution suggests that selective pressure on env, likely from new immune responses, is minimal when CD4 counts drop dramatically and remain low for extended periods of time.

Use of a novel assay based on intact recombinant viruses expressing green (EGFP) or red (DsRed2) fluorescent proteins to examine the contribution of pol and env genes to overall HIV-1 replicative fitness

Multiple studies have described a reduction in the replicative fitness of HIV-1 isolates harboring mutations that confer resistance to antiretroviral drugs. Contradictory results, however, have been obtained depending on the methodology used in each study (Quinones-Mateu, M.E., Arts, E.J., 2002. Fitness of drug resistant HIV-I: methodology and clinical implications. Drug Resist. Update 5, 224-233), affecting our understanding of the potential relationship of viral replicative fitness with HIV-1 disease. It has been demonstrated previously that both pol and env genes play a major role in HIV-1 replicative fitness of clinical isolates. Therefore, measuring clinically relevant replicative fitness using recombinant viruses where a single mutation and/or viral gene have been introduced does not seem like a reasonable approach in this era of multi-target antiretroviral therapy. A novel method was developed to measure HIV-1 replicative fitness based on recombinant viruses expressing the enhanced green fluorescent (EGFP) or the Discosoma sp. red fluorescent (DsRed2) proteins in a HIV-1NL4-3 backbone. Contrary to previous designs to analyze HIV-1 fitness, these replication competent viruses were created in an intact viral genetic background (without deleting or affecting the expression of any viral gene). This new system was used to evaluate the contribution of drug-resistance mutations in the pol and env genes to overall viral replicative fitness (in the presence and absence of drug pressure) using direct growth competition experiments. Mutations in pol showed a stronger effect on HIV-1 replicative fitness than mutations in the env gene associated with resistance to enfuvirtide, corroborating the plasticity of the later gene to accept mutations and the sensibility of the protease and reverse transcriptase enzymes to drug-associated primary mutations. In conclusion, a new protocol was used to measure HIV-1 replicative fitness in either the presence or absence of antiretroviral drugs, which may be used as a high-throughput assay to help us understand the clinical significance of viral fitness.

Antiretroviral therapy 2006: Pharmacology, applications, and special situations

As we approach the completion of the first 25 years of the human immunodeficiency virus (HIV) epidemic, there have been dramatic improvements in the care of patients with HIV infection. These have prolonged life and decreased morbidity. There are twenty currently available antiretrovirals approved in the United States for the treatment of this infection. The medications, including their pharmacokinetic properties, side effects, and dosing are reviewed. In addition, the current approach to the use of these medicines is discussed. We have included a section addressing common comorbid conditions including hepatitis B and C along with tuberculosis.

Virus fitness: concept, quantification, and application to HIV population dynamics

Viral fitness has been broadly studied during the past three decades, mainly to test evolutionary models and population theories difficult to analyze and interpret with more complex organisms. More recent studies, however, are focused in the role of fitness on viral transmission, pathogenesis, and drug resistance. Here, we used human immunodeficiency virus (HIV) as one of the most relevant models to evaluate the importance of viral quasispecies and fitness in HIV evolution, population dynamics, disease progression, and potential clinical implications.

Quantifying mixed populations of drug-resistant human immunodeficiency virus type 1

In order to survive prolonged treatment with antiretroviral nucleoside analogs, the human immunodeficiency virus type 1 (HIV-1) is selectively forced to acquire mutations in the reverse transcriptase (RT) gene. Some of these mutations are more common than others and have become markers for antiretroviral resistance. For the early detection of these markers, a novel MultiCode-RTx one-step testing system to rapidly and simultaneously characterize mixtures of HIV-1 targets was designed. For cDNA, nucleotide polymorphisms for codon M184V (ATG to GTG) and K65R (AAA to AGA) could be differentiated and quantified even when the population mixture varied as much as 1 to 10,000. Standard mixed-population curves using 1 to 100% of the mutant or wild type generated over 4 logs of total viral particle input did not affect the overall curves, making the method robust. The system was also applied to a small set of samples extracted from infected individuals on nucleoside reverse transcriptase inhibitor therapy. Of 13 samples tested, all were positive for HIV and 10 of the 13 genotypes determined were concordant with the line probe assay. MultiCode-RTx could be applied to other drug-selected mutations in the viral genome or for applications where single-base changes in DNA or RNA occur at frequencies reaching 0.01% to 1%, respectively.

Increased human immunodeficiency virus type 1 (HIV-1) env compartmentalization in the presence of HIV-1-associated dementia

The human immunodeficiency virus type 1 (HIV-1) surface Env protein has been implicated in the development of HIV-1-associated dementia (HAD). HIV-1 env diversity was analyzed by heteroduplex tracking assay in 27 infected subjects with various neurological statuses. env compartmentalization between the blood and cerebral spinal fluid (CSF) was apparent with all neurological categories. However, in subjects with HAD, significantly more CSF virus was represented by CNS-unique env variants. Variants specialized for replication in the CNS may play a larger role in the development of HAD. Alternatively, HAD may be associated with a more pronounced state of immunosuppression that permits more extensive replication and independent evolution within the CNS compartment.

Evolution of human immunodeficiency virus type 1 protease genotypes and phenotypes in vivo under selective pressure of the protease inhibitor ritonavir

We examined the population dynamics of human immunodeficiency virus type 1 pro variants during the evolution of resistance to the protease inhibitor ritonavir (RTV) in vivo. pro variants were followed in subjects who had added RTV to their previously failed reverse transcriptase inhibitor therapy using a heteroduplex tracking assay designed to detect common resistance-associated mutations. In most cases the initial variant appeared rapidly within 2 to 3 months followed by one or more subsequent population turnovers. Some of the subsequent transitions between variants were rapid, and some were prolonged with the coexistence of multiple variants. In several cases variants without resistance mutations persisted despite the emergence of new variants with an increasing number of resistance-associated mutations. Based on the rate of turnover of pro variants in the RTV-treated subjects we estimated that the mean fitness of newly emerging variants was increased 1.2-fold (range, 1.02 to 1.8) relative to their predecessors. A subset of pro genes was introduced into infectious molecular clones. The corresponding viruses displayed impaired replication capacity and reduced susceptibility to RTV. A subset of these clones also showed increased susceptibility to two nonnucleoside reverse transcriptase inhibitors and the protease inhibitor saquinavir. Finally, a significant correlation between the reduced replication capacity and reduced processing at the gag NC-p1 processing site was noted. Our results reveal a complexity of patterns in the evolution of resistance to a protease inhibitor. In addition, these results suggest that selection for resistance to one protease inhibitor can have pleiotropic effects that can affect fitness and susceptibility to other drugs.

Effect of a protease inhibitor-induced genetic bottleneck on human immunodeficiency virus type 1 env gene populations

The initiation of drug therapy or the addition of a new drug to preexisting therapy can have a significant impact on human immunodeficiency virus type 1 (HIV-1) populations within a person. Drug therapy directed at reverse transcriptase and protease can result in dramatic decreases in virus load, causing a contraction in the virus population that represents a potential genetic bottleneck as a subset of virus with genomes carrying resistance mutations repopulate the host. While this bottleneck exerts an effect directly on the region that is being targeted by the drugs, it also affects other regions of the viral genome. We have applied heteroduplex tracking assays (HTA) specific to variable regions 1 and 2 (V1/V2) and variable region 3 (V3) of the HIV-1 env gene to analyze the effect of a genetic bottleneck created by the selection of resistance to ritonavir, a protease inhibitor. Subjects were classified into groups on the basis of the extent of the initial drop in virus load and the duration of virus load reduction. Subjects with a strong initial drop in virus load exhibited a loss of heterogeneity in the env region at virus load rebound; in contrast, subjects with a weak initial drop in virus load exhibited little to no loss of heterogeneity at virus load rebound in either region of env examined. The duration of virus load reduction also affected env populations. Subjects that had prolonged reductions exhibited slower population diversification and the appearance of new V1/V2 species after rebound. The longer reduction of virus load in these subjects may have allowed for improved immune system function, which in turn could have selected for new escape mutants. Subjects with rapid rebound quickly reequilibrated the entry env variants back into the resistant population. When the pro gene developed further resistance mutations subsequent to virus load rebound, no changes were observed in V1/V2 or V3 populations, suggesting that the high virus loads allowed the env populations to reequilibrate rapidly. The rapid equilibration of env variants during pro gene sequence transitions at high virus load suggests that recombination is active in defining the HIV-1 sequence population. Conversely, part of the success of suppressive antiviral therapy may be to limit the potential for evolution through recombination, which requires dually infected cells.

Compartmentalized human immunodeficiency virus type 1 present in cerebrospinal fluid is produced by short-lived cells

Human immunodeficiency virus type 1 (HIV-1) invades the central nervous system (CNS) during primary infection and persists in this compartment by unknown mechanisms over the course of infection. In this study, we examined viral population dynamics in four asymptomatic subjects commencing antiretroviral therapy to characterize cellular sources of HIV-1 in the CNS. The inability to monitor viruses directly in the brain poses a major challenge in studying HIV-1 dynamics in the CNS. Studies of HIV-1 in cerebrospinal fluid (CSF) provide a useful surrogate for the sampling of virus in the CNS, but they are complicated by the fact that infected cells in local CNS tissues and in the periphery contribute to the population pool of HIV-1 in CSF. We utilized heteroduplex tracking assays to differentiate CSF HIV-1 variants that were shared with peripheral blood plasma from those that were compartmentalized in CSF and therefore presumably derived from local CNS tissues. We then tracked the relative decline of individual viral variants during the initial days of antiretroviral therapy. We found that HIV-1 variants compartmentalized in CSF declined rapidly during therapy, with maximum half-lives of approximately 1 to 3 days. These kinetics emulate the decline in HIV-1 produced from short-lived CD4+ T cells in the periphery, suggesting that a similarly short-lived, HIV-infected cell population exists within the CNS. We propose that short-lived CD4+ T cells trafficking between the CNS and the periphery play an important role in amplifying and maintaining HIV-1 populations in the CNS during the asymptomatic phase of infection.

Multiple V1/V2 env variants are frequently present during primary infection with human immunodeficiency virus type 1

Human immunodeficiency virus type 1 (HIV-1) exists as a complex population of multiple genotypic variants in persons with chronic infection. However, acute HIV-1 infection via sexual transmission is a low-probability event in which there is thought to be low genetic complexity in the initial inoculum. In order to assess the viral complexity present during primary HIV-1 infection, the V1/V2 and V3 variable regions of the env gene were examined by using a heteroduplex tracking assay (HTA) capable of resolving these genotypic variants. Blood plasma samples from 26 primary HIV-1-infected subjects were analyzed for their level of diversity. Half of the subjects had more than one V1/V2 viral variant during primary infection, indicating the frequent transmission of multiple variants. This observation is inconsistent with the idea of infrequent transmission based on a small transmitting inoculum of cell-free virus. In chronically infected subjects, the complexity of the viral populations was even greater in both the V1/V2 and the V3 regions than in acutely infected subjects, indicating that in spite of the presence of multiple variants in acute infection, the virus does pass through a genetic bottleneck during transmission. We also examined how well the infecting virus penetrated different anatomical compartments by using the HTA. Viral variants detected in blood plasma were compared to those detected in seminal plasma and/or cerebral spinal fluid of six individuals. The virus in each of these compartments was to a large extent identical to virus in blood plasma, a finding consistent with rapid penetration of the infecting variant(s). The low-probability transmission of multiple variants could be the result of transient periods of hyperinfectiousness or hypersusceptibility. Alternatively, the inefficient transfer of a multiply infected cell could account for both the low probability of transmission and the transfer of multiple variants.

Heterosexual transmission of multiple highly conserved viral variants in HIV-1 subtype C-infected seronegative women

Viral population diversity was assessed in samples collected from five HIV-infected women who were RNA positive and antibody negative. Similar to studies in men, highly conserved viral variants were detected (mean nucleotide diversity of 0.11% for p17p24, 0.32% for C2C3). In two individuals diversity was uncharacteristically lower in env C2C3 than in gag pl7p24, suggesting selection in env at this very early stage of infection.

Sequencing-based detection of low-frequency human immunodeficiency virus type 1 drug-resistant mutants by an RNA/DNA heteroduplex generator-tracking assay

Drug-resistant viruses may be present as minority variants during early treatment failures or following discontinuation of failed antiretroviral regimens. A limitation of the traditional direct PCR population sequencing method is its inability to detect human immunodeficiency virus type 1 (HIV-1) variants present at frequencies lower than 20%. A drug resistance genotyping assay based on the isolation and DNA sequencing of minority HIV protease variants is presented here. A multiple-codon-specific heteroduplex generator probe was constructed to improve the separation of HIV protease genes varying in sequence at 12 codons associated with resistance to protease inhibitors. Using an RNA molecule as probe allowed the simple sequencing of protease variants isolated as RNA/DNA heteroduplexes with different electrophoretic mobilities. The protease gene RNA heteroduplex generator-tracking assay (RNA-HTA) was tested on plasma quasispecies from 21 HIV-1-infected persons in whom one or more protease resistance mutations emerged during therapy or following initiation of salvage regimens. In 11 of 21 cases, RNA-HTA testing of virus from the first episode of virologic failure identified protease resistance mutations not seen by population-based PCR sequencing. In 8 of these 11 cases, all of the low-frequency drug resistance mutations detected exclusively by RNA-HTA during the first episode became detectable by population-based PCR sequencing at the later time point. Distinct sets of protease mutations could be linked on different genomes in patients with high-frequency protease gene lineages. The enhanced detection of minority drug resistance variants using a sequencing-based assay may improve the efficacy of genotype-assisted salvage therapies.

Role of minority populations of human immunodeficiency virus type 1 in the evolution of viral resistance to protease inhibitors

Human immunodeficiency virus type 1 (HIV-1) drug resistance results from the accumulation of mutations in the viral genes targeted by the drugs. These genetic changes, however, are commonly detected and monitored by techniques that only take into account the dominant population of plasma virus. Because HIV-1-infected patients harbor a complex and diverse mixture of virus populations, the mechanisms underlying the emergence and the evolution of resistance are not fully elucidated. Using techniques that allow the quantification of resistance mutations in minority virus species, we have monitored the evolution of resistance in plasma virus populations from patients failing protease inhibitor treatment. Minority populations with distinct resistance genotypes were detected in all patients throughout the evolution of resistance. The emergence of new dominant genotypes followed two possible mechanisms: (i) emergence of a new mutation in a currently dominant genotype and (ii) emergence of a new genotype derived from a minority virus species. In most cases, these population changes were associated with an increase in resistance at the expense of a reduction in replication capacity. Our findings provide a preliminary indication that minority viral species, which evolve independently of the majority virus population, can eventually become dominant populations, thereby serving as a reservoir of diversity and possibly accelerating the development of drug resistance.

Genetic and phenotypic analyses of human immunodeficiency virus type 1 escape from a small-molecule CCR5 inhibitor

We have described previously the generation of an escape variant of human immunodeficiency virus type 1 (HIV-1), under the selection pressure of AD101, a small molecule inhibitor that binds the CCR5 coreceptor (A. Trkola, S. E. Kuhmann, J. M. Strizki, E. Maxwell, T. Ketas, T. Morgan, P. Pugach, S. X. L. Wojcik, J. Tagat, A. Palani, S. Shapiro, J. W. Clader, S. McCombie, G. R. Reyes, B. M. Baroudy, and J. P. Moore, Proc. Natl. Acad. Sci. USA 99:395-400, 2002). The escape mutant, CC101.19, continued to use CCR5 for entry, but it was at least 20,000-fold more resistant to AD101 than the parental virus, CC1/85. We have now cloned the env genes from the the parental and escape mutant isolates and made chimeric infectious molecular clones that fully recapitulate the phenotypes of the corresponding isolates. Sequence analysis of the evolution of the escape mutants suggested that the most relevant changes were likely to be in the V3 loop of the gp120 glycoprotein. We therefore made a series of mutant viruses and found that full AD101 resistance was conferred by four amino acid changes in V3. Each change individually caused partial resistance when they were introduced into the V3 loop of a CC1/85 clone, but their impact was dependent on the gp120 context in which they were made. We assume that these amino acid changes alter how the HIV-1 Env complex interacts with CCR5. Perhaps unexpectedly, given the complete dependence of the escape mutant on CCR5 for entry, monomeric gp120 proteins expressed from clones of the fully resistant isolate failed to bind to CCR5 on the surface of L1.2-CCR5 cells under conditions where gp120 proteins from the parental virus and a partially AD101-resistant virus bound strongly. Hence, the full impact of the V3 substitutions may only be apparent at the level of the native Env complex.

Highly uneven distribution of tenofovir-selected simian immunodeficiency virus in different anatomical sites of rhesus macaques

Antiviral tenofovir monotherapy was used to determine whether drug-selected simian immunodeficiency virus (SIV) variants replaced their wild-type progenitors at the same rate in different tissues of six rhesus macaques. The relative frequencies of drug-resistant and wild-type genotypes were measured longitudinally in blood and in 23 lymphoid and nonlymphoid tissues collected at necropsy. The mutant/wild-type genotype ratio was measured using a heteroduplex tracking assay targeting tenofovir-selected SIV reverse transcriptase codons. After the initiation of tenofovir treatment in animals with high steady-state viremia levels, resistant genotypes emerged in the plasma within 1 to 8 weeks and in five of six cases reached frequencies of nearly 100% within 4 to 25 weeks. The appearance of tenofovir-resistant genotypes in peripheral blood mononuclear cell (PBMC) DNA was generally delayed by 1 to 2 weeks and in one case was completely absent. Necropsies performed 8 to 55 weeks after the initiation of tenofovir treatment showed the frequency of resistant SIV genotypes to be generally higher in tissue RNA than DNA fractions. The frequency of drug-resistant genotypes varied widely between anatomical sites, including different lymph nodes of the same animal. Except for the epidydimis, the tissues with the lowest rates of proviral replacement by tenofovir-resistant genotypes differed between animals. The highly uneven distribution of tenofovir-resistant genotypes in different tissues seen shortly after the initiation of tenofovir monotherapy may reflect differences in local antiviral drug selection pressures and/or the stochastic effect of small effective populations of drug-resistant variants randomly seeding different anatomical sites early in therapy.

Turnover of env variable region 1 and 2 genotypes in subjects with late-stage human immunodeficiency virus type 1 infection

The env gene of human immunodeficiency virus type 1 (HIV-1) includes some of the most genetically diverse regions of the viral genome, which are called variable regions 1 through 5 (V1 through V5). We have developed a heteroduplex tracking assay to detect changes in variable regions 1 and 2 of env (V1/V2-HTA). Using sequences from two molecular clones as probes, we have studied the nature of longitudinal virus population changes in a cohort of HIV-1-infected subjects. Viral sequences present in 21 subjects with late-stage HIV-1 infection were initially screened for stability of the virus population by V1/V2-HTA. The virus populations at entry comprised an average of five coexisting V1/V2 genotypic variants (as identified by HTA). Eight of the 21 subjects were examined in detail because of the dynamic behavior of their env variants over an approximately 9-month period. In each of these cases we detected a single discrete transition of V1/V2 genotypes based on monthly sampling. The major V1/V2 genotypes (those present at >10% abundance) from the eight subjects were cloned and sequenced to define the nature of V1/V2 variability associated with a discrete transition. Based on a comparison of V1/V2 genotypic variants present at entry with the newly emerged variants we categorized the newly emerged variants into two groups: variants without length differences and variants with length differences. Variants without length differences had fewer nucleotide substitutions, with the changes biased to either V1 or V2, suggestive of recent evolutionary events. Variants with length differences included ones with larger numbers of changes that were distributed, suggestive of recall of older genotypes. Most length differences were located in domains where the codon motif AVT (V = A, G, C) had become enriched and fixed. Finally, recombination events were detected in two subjects, one of which resulted in the reassortment of V1 and V2 regions. We suggest that turnover in V1/V2 populations was largely driven by selection on either V1 or V2 and that escape was accomplished either through changes focused in the region under selection or by the appearance of a highly divergent variant.

Selection of high-level resistance to human immunodeficiency virus type 1 protease inhibitors

Protease inhibitors represent some of the most potent agents available for therapeutic strategies designed to inhibit human immunodeficiency virus type 1 (HIV-1) replication. Under certain circumstances the virus develops resistance to the inhibitor, thereby negating the benefits of this therapy. We have carried out selections for high-level resistance to each of three protease inhibitors (indinavir, ritonavir, and saquinavir) in cell culture. Mutations accumulated over most of the course of the increasing selective pressure. There was significant overlap in the identity of the mutations selected with the different inhibitors, and this gave rise to high levels of cross-resistance. Virus particles from the resistant variants all showed defects in processing at the NC/p1 protease cleavage site in Gag. Selections with pairs of inhibitors yielded similar patterns of resistance mutations. A virus that could replicate at near-toxic levels of the three protease inhibitors combined was selected. The pro sequence of this virus was similar to that of the viruses that had been selected for high-level resistance to each of the drugs singly. Finally, a molecular clone carrying the eight most common resistance mutations seen in these selections was characterized. The sequence of this virus was relatively stable during selection for revertants in spite of displaying poor processing at the NC/p1 site and having significantly reduced fitness. These results reveal patterns of drug resistance that extend to near the limits of attainable selective pressure with these inhibitors and confirm the patterns of cross-resistance for these three inhibitors and the attenuation of virion protein processing and fitness that accompanies high-level resistance.

Fitness of drug resistant HIV-1: methodology and clinical implications

Recent studies of human immunodeficiency virus type 1 (HIV-1) fitness have examined the potential relationship with plasma viral load, drug resistance, and disease progression. For example, treatment of HIV-1 infected individuals with antiretroviral drugs may result in the selection and emergence of inhibitor-resistant variants with reduced replicative capacity. However, it is still unclear whether in vitro HIV-1 fitness has any direct relationship to in vivo disease progression or treatment success. A related question is which in vitro assay of viral fitness is the most appropriate for comparison with in vivo HIV-1 fitness. Characterization of the relative viral fitness of drug-resistant HIV-1 strains may lead to a better understanding of whether or not less fit viruses pose a clinical benefit to the patient.

Mutagenically separated PCR assay for rapid detection of M41L and K70R zidovudine resistance mutations in CRF01_AE (subtype E) human immunodeficiency virus type 1

A rapid zidovudine (ZDV) resistance genotypic assay was developed based on the mutagenically separated PCR (MS-PCR) technique to detect two ZDV-resistant mutations, M41L and K70R in CRF01_AE (subtype E). Endpoint dilution analysis revealed that the newly constructed MS-PCR assay could successfully detect three to nine copies of human immunodeficiency virus type 1 template RNA. The test against wild-type and mutant template mixtures in different ratios demonstrated that the assay could detect 10% minor population, at least. Fifty-one subtype E clinical samples were analyzed by the newly constructed MS-PCR assay and direct nucleotide sequencing. The concordance of the two assays was 92 and 100% in codons 41 and 70, respectively. The MS-PCR assay is a rapid, simple, and inexpensive assay that is highly sensitive in detecting mutant targets, including minor populations. Thus, it could be used as a powerful tool for epidemiological surveillance of drug-resistant mutations in developing countries.

Identification of Epstein-Barr virus strain variants in hairy leukoplakia and peripheral blood by use of a heteroduplex tracking assay

Epstein-Barr virus (EBV) strains can be distinguished by specific sequence variations in the LMP1 gene. In this study, a heteroduplex tracking assay (HTA) specific for LMP1 was developed to precisely identify the prototypic undeleted strain B958, other undeleted strains (Ch2, AL, NC, and Med-), and strains with the 30-bp deletion (Med+ and Ch1). This technique also provides an estimate of the relative abundance of strains in patient samples. In this study, EBV strains were identified in 25 hairy leukoplakia (HLP) biopsies and six matched peripheral blood samples and throat washes with the LMP1-HTA. To investigate the relationship of the virus found in the peripheral blood to that in the HLP lesion, the strain variants in the peripheral blood B lymphocytes and those present within the epithelial cells in the HLP lesion and in throat washes were identified. In many of the subjects, compartmental differences in the EBV strain profiles in the oral cavity and peripheral blood were readily apparent. The throat wash specimens usually had a strain profile similar to that within the corresponding HLP sample, which was distinct from the strain profile detected in the peripheral blood. These analyses reveal that the nature of EBV infection can be very dynamic, with changes in relative strain abundance over time as well as the appearance of new strains. The patterns of abundance in the blood and oral cavity provide evidence for compartmentalization and for the transmission of strains between the blood and oropharynx.

Nelfinavir-resistant, amprenavir-hypersusceptible strains of human immunodeficiency virus type 1 carrying an N88S mutation in protease have reduced infectivity, reduced replication capacity, and reduced fitness and process the Gag polyprotein precursor aberrantly

The evolution of human immunodeficiency virus type 1 (HIV-1) strains with reduced susceptibility to protease inhibitors (PIs) is a major cause of PI treatment failure. A subset of subjects failing a therapy regimen containing the PI nelfinavir developed mutations at position 88 in the protease region. The N88S mutation occurring in some of these subjects induces amprenavir hypersusceptibility and a reduction of fitness and replication capacity. Here we demonstrate that substitutions L63P and V77I in protease, in combination, partially compensate for the loss of fitness, loss of replication capacity, loss of specific infectivity, and aberrant Gag processing induced by the N88S mutation. In addition, these mutations partially ablate amprenavir hypersusceptibility. Addition of mutation M46L to a strain harboring mutations L63P, V77I, and N88S resulted in a reduction of fitness and infectivity without changing Gag-processing efficiency, while amprenavir hypersusceptibility was further diminished. The ratio of reverse transcriptase activity to p24 protein was reduced in this strain compared to that in the other variants, suggesting that the M46L effect on fitness occurred through a mechanism different from a Gag-processing defect. We utilized these mutant strains to undertake a systematic comparison of indirect, single, cycle-based measures of fitness with direct, replication-based fitness assays and demonstrated that both yield consistent results. However, we observed that the magnitude of the fitness loss for one of the mutants varied depending on the assay used.

Population genetic analysis of the protease locus of human immunodeficiency virus type 1 quasispecies undergoing drug selection, using a denaturing gradient-heteroduplex tracking assay

Monitoring the evolution of human immunodeficiency virus type 1 (HIV-1) drug resistance requires measuring the frequency of closely related genetic variants making up the complex viral quasispecies found in vivo. In order to resolve both major and minor (>/=2%) protease gene variants differing by one or more nucleotide substitutions, we analyzed PCR products derived from plasma viral quasispecies by using a combination of denaturing gradient gel electrophoresis and DNA heteroduplex tracking assays. Correct population sampling of the high level of genetic diversity present within viral quasispecies could be documented by parallel analysis of duplicate, independently generated PCR products. The composition of genetically complex protease gene quasispecies remained constant over short periods of time in the absence of treatment and while plasma viremia fell >100-fold following the initiation of protease inhibitor ritonavir monotherapy. Within a month of initiating therapy, a strong reduction in the genetic diversity of plasma viral populations at the selected protease locus was associated with rising plasma viremia and the emergence of drug resistance. The high levels of protease genetic diversity seen before treatment reemerged only months later. In one patient, reduction in genetic diversity at the protease gene was observed concomitantly with an increase in diversity at the envelope gene (E. L. Delwart, P. Heng, A. Neumann, and M. Markowitz, J. Virol. 72:2416-2421, 1998), indicating that opposite population genetic changes can take place in different HIV-1 loci. The rapid emergence of drug-resistant HIV-1 was therefore associated with a strong, although only transient, reduction in genetic diversity at the selected locus. The denaturing gradient-heteroduplex tracking assay is a simple method for the separation and quantitation of very closely related, low-frequency, genetic variants within complex viral populations.