Technical and regulatory shortcomings of the TaqMan version 1 HIV viral load assay.

BACKGROUND: The lower limit of detection of the original Roche Amplicor HIV plasma viral load (pVL) assay (50 copies/mL) has defined HIV treatment success. The Amplicor assay, however, has been replaced by the Roche TaqMan assay(s). Changes to the limits of detection and calibration have not been validated for clinical utility. Sudden increases in the number of patients with detectable pVL have been reported following the introduction of the TaqMan version 1 assay. METHODS: Between October 2009 and April 2010 all routine pVL samples from British Columbia, Canada, with 40-250 copies/mL by TaqMan were re-tested by Amplicor (N = 1198). Subsequent short-term virological and resistance outcomes were followed in patients with unchanged therapy (N = 279; median 3.2 months follow-up). RESULTS: TaqMan and Amplicor values correlated poorly at low pVL values. Low-level pVL by TaqMan was not associated with impending short-term virological failure; only 17% of patients with 40-250 copies/mL by TaqMan had detectable pVL by Amplicor at follow-up. During the follow-up period only 20% of patients had an increase in pVL by TaqMan (median [IQR]: 80 [36-283] copies/mL). In addition, in ~2.4% of samples pVL was dramatically underestimated by TaqMan due to poor binding of the proprietary TaqMan primers. CONCLUSIONS: The replacement of Amplicor with the TaqMan assay has altered the previously accepted definition of HIV treatment failure without any evidence to support the clinical relevance of the new definition. Given the systematic differences in measurement in the low pVL range the British Columbia HIV treatment guidelines now use a threshold of >250 copies/mL by TaqMan to define treatment failure.

Automating HIV drug resistance genotyping with RECall, a freely accessible sequence analysis tool.

Genotypic HIV drug resistance testing is routinely used to guide clinical decisions. While genotyping methods can be standardized, a slow, labor-intensive, and subjective manual sequence interpretation step is required. We therefore performed external validation of our custom software RECall, a fully automated sequence analysis pipeline. HIV-1 drug resistance genotyping was performed on 981 clinical samples at the Stanford Diagnostic Virology Laboratory. Sequencing trace files were first interpreted manually by a laboratory technician and subsequently reanalyzed by RECall, without intervention. The relative performances of the two methods were assessed by determination of the concordance of nucleotide base calls, identification of key resistance-associated substitutions, and HIV drug resistance susceptibility scoring by the Stanford Sierra algorithm. RECall is freely available at http://pssm.cfenet.ubc.ca. In total, 875 of 981 sequences were analyzed by both human and RECall interpretation. RECall analysis required minimal hands-on time and resulted in a 25-fold improvement in processing speed (∼150 technician-hours versus ∼6 computation-hours). Excellent concordance was obtained between human and automated RECall interpretation (99.7% agreement for >1,000,000 bases compared). Nearly all discordances (99.4%) were due to nucleotide mixtures being called by one method but not the other. Similarly, 98.6% of key antiretroviral resistance-associated mutations observed were identified by both methods, resulting in 98.5% concordance of resistance susceptibility interpretations. This automated sequence analysis tool provides both standardization of analysis and a significant improvement in data workflow. The time-consuming, error-prone, and dreadfully boring manual sequence analysis step is replaced with a fully automated system without compromising the accuracy of reported HIV drug resistance data.

Increasingly successful highly active antiretroviral therapy delays the emergence of new HLA class I-associated escape mutations in HIV-1.

BACKGROUND: HLA class I-restricted cytotoxic T lymphocytes and highly active antiretroviral therapy (HAART) exert strong selective pressures on human immunodeficiency virus type 1 (HIV-1), leading to escape mutations compromising virologic control. Immune responses continue to shape HIV-1 evolution after HAART initiation, but the extent and rate at which this occurs remain incompletely quantified. Here, we characterize the incidence and clinical correlates of HLA-associated evolution in HIV-1 Pol after HAART initiation in a large, population-based observational cohort. METHODS: British Columbia HAART Observational, Medical Evaluation and Research cohort participants with available HLA class I types and longitudinal posttherapy protease/reverse transcriptase sequences were studied (n = 619; median, 5 samples per patient and 5.2 years of follow-up). HLA-associated polymorphisms were defined according to published reference lists. Rates and correlates of immune-mediated HIV-1 evolution were investigated using multivariate Cox proportional hazard models incorporating baseline and time-dependent plasma viral load and CD4 response data. RESULTS: New HLA-associated escape events were observed in 269 (43%) patients during HAART and occurred at 49 of 63 (78%) investigated immune-associated sites in Pol. In time-dependent analyses adjusting for baseline factors, poorer virologic, but not immunologic, response to HAART was associated with increased risk of immune escape of 1.9-fold per log(10) viral load increment (P < .0001). Reversion of escape mutations following HAART initiation was extremely rare. CONCLUSIONS: HLA-associated HIV-1 evolution continues during HAART to an extent that is inversely related to the virologic success of therapy. Minimizing the degree of immune escape could represent a secondary benefit of effective HAART.

"Dynamic range" of inferred phenotypic HIV drug resistance values in clinical practice.

BACKGROUND: 'Virtual' or inferred phenotypes (vPhenotypes) are commonly used to assess resistance to antiretroviral agents in patients failing therapy. In this study, we provide a clinical context for understanding vPhenotype values. METHODS: All HIV-infected persons enrolled in the British Columbia Drug Treatment Program with a baseline plasma viral load (pVL) and follow-up genotypic resistance and pVL results were included up to October 29, 2008 (N = 5,277). Change from baseline pVL was determined as a function of Virco vPhenotype, and the "dynamic range" (defined here by the 10th and 90th percentiles for fold-change in IC50 amongst all patients) was estimated from the distribution of vPhenotye fold-changes across the cohort. RESULTS: The distribution of vPhenotypes from a large cohort of HIV patients who have failed therapy are presented for all available antiretroviral agents. A maximum change in IC50 of at least 13-fold was observed for all drugs. The dideoxy drugs, tenofovir and most PIs exhibited small "dynamic ranges" with values of <4-fold change observed in > 99% of samples. In contrast, zidovudine, lamivudine, emtricitabine and the non-nucleoside reverse transcriptase inihibitors (excluding etravirine) had large dynamic ranges. CONCLUSION: We describe the populational distribution of vPhenotypes such that vPhenotype results can be interpreted relative to other patients in a drug-specific manner.

Gene discovery in EST sequences from the wheat leaf rust fungus Puccinia triticina sexual spores, asexual spores and haustoria, compared to other rust and corn smut fungi.

BACKGROUND: Rust fungi are biotrophic basidiomycete plant pathogens that cause major diseases on plants and trees world-wide, affecting agriculture and forestry. Their biotrophic nature precludes many established molecular genetic manipulations and lines of research. The generation of genomic resources for these microbes is leading to novel insights into biology such as interactions with the hosts and guiding directions for breakthrough research in plant pathology. RESULTS: To support gene discovery and gene model verification in the genome of the wheat leaf rust fungus, Puccinia triticina (Pt), we have generated Expressed Sequence Tags (ESTs) by sampling several life cycle stages. We focused on several spore stages and isolated haustorial structures from infected wheat, generating 17,684 ESTs. We produced sequences from both the sexual (pycniospores, aeciospores and teliospores) and asexual (germinated urediniospores) stages of the life cycle. From pycniospores and aeciospores, produced by infecting the alternate host, meadow rue (Thalictrum speciosissimum), 4,869 and 1,292 reads were generated, respectively. We generated 3,703 ESTs from teliospores produced on the senescent primary wheat host. Finally, we generated 6,817 reads from haustoria isolated from infected wheat as well as 1,003 sequences from germinated urediniospores. Along with 25,558 previously generated ESTs, we compiled a database of 13,328 non-redundant sequences (4,506 singlets and 8,822 contigs). Fungal genes were predicted using the EST version of the self-training GeneMarkS algorithm. To refine the EST database, we compared EST sequences by BLASTN to a set of 454 pyrosequencing-generated contigs and Sanger BAC-end sequences derived both from the Pt genome, and to ESTs and genome reads from wheat. A collection of 6,308 fungal genes was identified and compared to sequences of the cereal rusts, Puccinia graminis f. sp. tritici (Pgt) and stripe rust, P. striiformis f. sp. tritici (Pst), and poplar leaf rust Melampsora species, and the corn smut fungus, Ustilago maydis (Um). While extensive homologies were found, many genes appeared novel and species-specific; over 40% of genes did not match any known sequence in existing databases. Focusing on spore stages, direct comparison to Um identified potential functional homologs, possibly allowing heterologous functional analysis in that model fungus. Many potentially secreted protein genes were identified by similarity searches against genes and proteins of Pgt and Melampsora spp., revealing apparent orthologs. CONCLUSIONS: The current set of Pt unigenes contributes to gene discovery in this major cereal pathogen and will be invaluable for gene model verification in the genome sequence.

Improved detection of CXCR4-using HIV by V3 genotyping: application of population-based and "deep" sequencing to plasma RNA and proviral DNA.

BACKGROUND: Tropism testing should rule out CXCR4-using HIV before treatment with CCR5 antagonists. Currently, the recombinant phenotypic Trofile assay (Monogram) is most widely utilized; however, genotypic tests may represent alternative methods. METHODS: Independent triplicate amplifications of the HIV gp120 V3 region were made from either plasma HIV RNA or proviral DNA. These underwent standard, population-based sequencing with an ABI3730 (RNA n = 63; DNA n = 40), or "deep" sequencing with a Roche/454 Genome Sequencer-FLX (RNA n = 12; DNA n = 12). Position-specific scoring matrices (PSSMX4/R5) (-6.96 cutoff) and geno2pheno[coreceptor] (5% false-positive rate) inferred tropism from V3 sequence. These methods were then independently validated with a separate, blinded dataset (n = 278) of screening samples from the maraviroc MOTIVATE trials. RESULTS: Standard sequencing of HIV RNA with PSSM yielded 69% sensitivity and 91% specificity, relative to Trofile. The validation dataset gave 75% sensitivity and 83% specificity. Proviral DNA plus PSSM gave 77% sensitivity and 71% specificity. "Deep" sequencing of HIV RNA detected >2% inferred-CXCR4-using virus in 8/8 samples called non-R5 by Trofile, and <2% in 4/4 samples called R5. CONCLUSIONS: Triplicate analyses of V3 standard sequence data detect greater proportions of CXCR4-using samples than previously achieved. Sequencing proviral DNA and "deep" V3 sequencing may also be useful tools for assessing tropism.

Improved virological outcomes in British Columbia concomitant with decreasing incidence of HIV type 1 drug resistance detection.

BACKGROUND: There have been limited studies evaluating temporal changes in the incidence of detection of drug resistance among human immunodeficiency virus type 1 (HIV-1) isolates and concomitant changes in plasma HIV load for treated individuals in a population-wide setting. METHODS: Longitudinal plasma viral load and genotypic resistance data were obtained from patients receiving antiretroviral therapy from the British Columbia Drug Treatment Program from July 1996 through December 2008. A total of 24,652 resistance tests were available from 5422 individuals. The incidence of successful plasma viral load suppression and of resistance to each of 3 antiretroviral categories (nucleoside/nucleotide reverse-transcriptase inhibitors, nonnucleoside reverse transcriptase inhibitors, and protease inhibitors) was calculated for the population receiving therapy. RESULTS: There has been a drastic decrease in the incidence of new cases of HIV-1 drug resistance in individuals followed during 1996-2008. In 1997, the incidence rate of any newly detected resistance was 1.73 cases per 100 person-months of therapy, and by 2008, the incidence rate had decreased >12-fold, to 0.13 cases per 100 person-months of therapy. This decrease in the incidence of resistance has occurred at an exponential rate, with half-times on the order of 2-3 years. Concomitantly, the proportion of individuals with plasma viral load suppression has increased linearly over time (from 64.7% with HIV RNA levels <50 copies/mL in 2000 to 87.0% in 2008; R2=0.97; P<.001). CONCLUSIONS: Our results suggest an increasing effectiveness of highly active antiretroviral therapy at the populational level. The vast majority of treated patients in British Columbia now have either suppressed plasma viral load or drug-susceptible HIV-1, according to their most recent test results.

Silent mutations are selected in HIV-1 reverse transcriptase and affect enzymatic efficiency.

OBJECTIVE: Missense mutations in HIV-1 reverse transcriptase are frequently selected in response to therapy; we examined whether silent mutations were also selected for by HIV therapy. DESIGN: Retrospective, observational analysis. Biochemical assays. METHODS: A comparison of the reverse transcriptase gene, from antiretroviral- naive (N = 812) and experienced individuals (N = 2212), reveals two silent mutations (K65K and K66K) that are strongly associated with treatment experience. To assess reverse transcription efficiency, steady-state kinetic assays were carried out using recombinant purified HIV-1 reverse transcriptase and a series of synthetic RNA/DNA template/primer substrates. The RNA templates spanned codons 60-77 in the reverse transcriptase and included different combinations of mutations at codons 65, 66, 67, and 70. RESULTS: Silent AAG mutations (or mixtures) at reverse transcriptase codons 65 and/or 66 were observed in 812 samples from 351 patients and 2129 samples from 829 patients, respectively. In clade B samples, there was a very strong relationship between the silent mutations and the thymidine analogue mutations, in particular D67N. Steady-state kinetic experiments demonstrated that HIV-1 reverse transcriptase exhibited a strong tendency to pause and/or dissociate at codons 65 and 66 on RNA templates that contained the D67N and K70R mutations. However, when the K66 or K66 AAA to AAG mutations were added to the background of the 67 and 70 mutational changes, these pausing and/or dissociation events were largely alleviated. CONCLUSION: Silent mutations at codons 65 and/or 66 are strongly coselected with thymidine analogue mutations. These data provide the first evidence for an RNA-level mechanism of direct relevance to drug resistance.

Connection domain mutations N348I and A360V in HIV-1 reverse transcriptase enhance resistance to 3'-azido-3'-deoxythymidine through both RNase H-dependent and -independent mechanisms.

Thymidine analogue-associated mutations (TAMs) in reverse transcriptase (RT) of the human immunodeficiency virus type 1 (HIV-1) cause resistance to 3'-azido-3'-deoxythymidine (AZT) through excision of the incorporated monophosphate. Mutations in the connection domain of HIV-1 RT can augment AZT resistance. It has been suggested that these mutations compromise RNase H cleavage, providing more time for AZT excision to occur. However, the underlying mechanism remains elusive. Here, we focused on connection mutations N348I and A360V that are frequently observed in clinical samples of treatment-experienced patients. We show that both N348I and A360V, in combination with TAMs, decrease the efficiency of RNase H cleavage and increase excision of AZT in the presence of the pyrophosphate donor ATP. The TAMs/N348I/A360V mutant accumulates transiently formed, shorter hybrids that can rebind to RT before the template is irreversibly degraded. These hybrids dissociate selectively from the RNase H-competent complex, whereas binding in the polymerase-competent mode is either not affected with N348I or modestly improved with A360V. Both connection domain mutations can compensate for TAM-mediated deficits in processive DNA synthesis, and experiments with RNase H negative mutant enzymes confirm an RNase H-independent contribution to increased levels of resistance to AZT. Moreover, the combination of diminished RNase H cleavage and increased processivity renders the use of both PP(i) and ATP advantageous, whereas classic TAMs solely enhance the ATP-dependent reaction. Taken together, our findings demonstrate that distinct, complementary mechanisms can contribute to higher levels of excision of AZT, which in turn can amplify resistance to this drug.

N348I in the connection domain of HIV-1 reverse transcriptase confers zidovudine and nevirapine resistance.

BACKGROUND: The catalytically active 66-kDa subunit of the human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) consists of DNA polymerase, connection, and ribonuclease H (RNase H) domains. Almost all known RT inhibitor resistance mutations identified to date map to the polymerase domain of the enzyme. However, the connection and RNase H domains are not routinely analysed in clinical samples and none of the genotyping assays available for patient management sequence the entire RT coding region. The British Columbia Centre for Excellence in HIV/AIDS (the Centre) genotypes clinical isolates up to codon 400 in RT, and our retrospective statistical analyses of the Centre's database have identified an N348I mutation in the RT connection domain in treatment-experienced individuals. The objective of this multidisciplinary study was to establish the in vivo relevance of this mutation and its role in drug resistance. METHODS AND FINDINGS: The prevalence of N348I in clinical isolates, the time taken for it to emerge under selective drug pressure, and its association with changes in viral load, specific drug treatment, and known drug resistance mutations was analysed from genotypes, viral loads, and treatment histories from the Centre's database. N348I increased in prevalence from below 1% in 368 treatment-naïve individuals to 12.1% in 1,009 treatment-experienced patients (p = 7.7 x 10(-12)). N348I appeared early in therapy and was highly associated with thymidine analogue mutations (TAMs) M41L and T215Y/F (p < 0.001), the lamivudine resistance mutations M184V/I (p < 0.001), and non-nucleoside RTI (NNRTI) resistance mutations K103N and Y181C/I (p < 0.001). The association with TAMs and NNRTI resistance mutations was consistent with the selection of N348I in patients treated with regimens that included both zidovudine and nevirapine (odds ratio 2.62, 95% confidence interval 1.43-4.81). The appearance of N348I was associated with a significant increase in viral load (p < 0.001), which was as large as the viral load increases observed for any of the TAMs. However, this analysis did not account for the simultaneous selection of other RT or protease inhibitor resistance mutations on viral load. To delineate the role of this mutation in RT inhibitor resistance, N348I was introduced into HIV-1 molecular clones containing different genetic backbones. N348I decreased zidovudine susceptibility 2- to 4-fold in the context of wild-type HIV-1 or when combined with TAMs. N348I also decreased susceptibility to nevirapine (7.4-fold) and efavirenz (2.5-fold) and significantly potentiated resistance to these drugs when combined with K103N. Biochemical analyses of recombinant RT containing N348I provide supporting evidence for the role of this mutation in zidovudine and NNRTI resistance and give some insight into the molecular mechanism of resistance. CONCLUSIONS: This study provides the first in vivo evidence that treatment with RT inhibitors can select a mutation (i.e., N348I) outside the polymerase domain of the HIV-1 RT that confers dual-class resistance. Its emergence, which can happen early during therapy, may significantly impact on a patient's response to antiretroviral therapies containing zidovudine and nevirapine. This study also provides compelling evidence for investigating the role of other mutations in the connection and RNase H domains in virological failure.

Transmission of drug-resistant HIV-1 from an infected individual to a caregiver.

OBJECTIVE: To describe a suspected case of intrafamilial transmission of drug-resistant HIV-1 from an infected individual to his caregiver. METHODS: HIV-1 protease, reverse transcriptase, nef and gp41 DNA sequences were determined from plasma samples after RNA extraction and nested RT-PCR. Phylogenetic trees were generated using neighbour-joining methods. RESULTS: HIV-1 infection was diagnosed in an individual (index) following the death of her HIV-infected adopted son (source). The index case cared for her son in the 9 months prior to his death. No obvious instances of contact with bodily secretions or other traditional risk factors for HIV transmission were identified, although the index case had mild eczema. At time of death, the source case's HIV plasma viral load was > 1,000,000 HIV copies/ml. Genotyping of the protease and reverse transcriptase of both the source and index samples revealed similar drug-resistance mutations despite the index case being antiretroviral-naive. Both source and index V3 genotypes were consistent with syncytium-inducing virus. The mean pairwise amino acid identity between the source and index samples was 97.0%, 99.2%, 92.6% and 94.6% in protease, reverse transcriptase, nef and gp41, respectively. The source and index sequences clustered together on phylogenetic trees when compared with 50 randomly selected sequences from British Columbia. CONCLUSION: Although exceptionally rare, our case demonstrates transmission of drug-resistant HIV-1 from an infected individual to his caregiver in the absence of traditional risk factors. It should be noted that the source had a very high viral load during the terminal phase of his illness.

Effects of human leukocyte antigen class I genetic parameters on clinical outcomes and survival after initiation of highly active antiretroviral therapy.

BACKGROUND: Human leukocyte antigen (HLA) class I variation influences the progression of untreated human immunodeficiency virus (HIV) disease; however, it is not known whether HLA class I variation may influence clinical outcomes after initiation of highly active antiretroviral therapy (HAART). METHODS: Associations between HLA class I genotypes and pretherapy clinical parameters were investigated in a cohort of 765 antiretroviral-naive adults initiating HAART. Cox proportional hazards regression was used to investigate the effects of HLA class I genotypes on time to suppression of the viral load to <500 HIV RNA copies/mL, time to an increase in the CD4 cell count to >100 cells/mm(3) above the baseline count, and time to nonaccidental death over a >5-year period after initiation of HAART. RESULTS: Homozygosity at any HLA class I locus and possession of common HLA alleles were associated with a higher pretherapy viral load (P<.05). In multivariate analyses controlling for sociodemographic and clinical parameters at baseline, HLA class I homozygosity was significantly associated with a poorer CD4 cell response (P=.008), whereas possession of uncommon HLA alleles was associated with slower virologic suppression after initiation of HAART (P=.02). We observed no significant association between HLA parameters and time to nonaccidental death after initiation of HAART (P>.05, univariate analysis). CONCLUSION: HLA class I zygosity-dependent and frequency-dependent effects may influence short-term HAART outcomes, and, thus, they deserve further investigation. No effects of these HLA parameters on survival after initiation of HAART were observed.

Genomics of hybrid poplar (Populus trichocarpax deltoides) interacting with forest tent caterpillars (Malacosoma disstria): normalized and full-length cDNA libraries, expressed sequence tags, and a cDNA microarray for the study of insect-induced defences in poplar.

As part of a genomics strategy to characterize inducible defences against insect herbivory in poplar, we developed a comprehensive suite of functional genomics resources including cDNA libraries, expressed sequence tags (ESTs) and a cDNA microarray platform. These resources are designed to complement the existing poplar genome sequence and poplar (Populus spp.) ESTs by focusing on herbivore- and elicitor-treated tissues and incorporating normalization methods to capture rare transcripts. From a set of 15 standard, normalized or full-length cDNA libraries, we generated 139,007 3'- or 5'-end sequenced ESTs, representing more than one-third of the c. 385,000 publicly available Populus ESTs. Clustering and assembly of 107,519 3'-end ESTs resulted in 14,451 contigs and 20,560 singletons, altogether representing 35,011 putative unique transcripts, or potentially more than three-quarters of the predicted c. 45,000 genes in the poplar genome. Using this EST resource, we developed a cDNA microarray containing 15,496 unique genes, which was utilized to monitor gene expression in poplar leaves in response to herbivory by forest tent caterpillars (Malacosoma disstria). After 24 h of feeding, 1191 genes were classified as up-regulated, compared to only 537 down-regulated. Functional classification of this induced gene set revealed genes with roles in plant defence (e.g. endochitinases, Kunitz protease inhibitors), octadecanoid and ethylene signalling (e.g. lipoxygenase, allene oxide synthase, 1-aminocyclopropane-1-carboxylate oxidase), transport (e.g. ABC proteins, calreticulin), secondary metabolism [e.g. polyphenol oxidase, isoflavone reductase, (-)-germacrene D synthase] and transcriptional regulation [e.g. leucine-rich repeat transmembrane kinase, several transcription factor classes (zinc finger C3H type, AP2/EREBP, WRKY, bHLH)]. This study provides the first genome-scale approach to characterize insect-induced defences in a woody perennial providing a solid platform for functional investigation of plant-insect interactions in poplar.

Discordances between interpretation algorithms for genotypic resistance to protease and reverse transcriptase inhibitors of human immunodeficiency virus are subtype dependent.

The major limitation of drug resistance genotyping for human immunodeficiency virus remains the interpretation of the results. We evaluated the concordance in predicting therapy response between four different interpretation algorithms (Rega 6.3, HIVDB-08/04, ANRS [07/04], and VGI 8.0). Sequences were gathered through a worldwide effort to establish a database of non-B subtype sequences, and demographic and clinical information about the patients was gathered. The most concordant results were found for nonnucleoside reverse transcriptase (RT) inhibitors (93%), followed by protease inhibitors (84%) and nucleoside RT inhibitor (NRTIs) (76%). For therapy-naive patients, for nelfinavir, especially for subtypes C and G, the discordances were driven mainly by the protease (PRO) mutational pattern 82I/V + 63P + 36I/V for subtype C and 82I + 63P + 36I + 20I for subtype G. Subtype F displayed more discordances for ritonavir in untreated patients due to the combined presence of PRO 20R and 10I/V. In therapy-experienced patients, subtype G displayed a lot of discordances for saquinavir and indinavir due to mutational patterns involving PRO 90 M and 82I. Subtype F had more discordance for nelfinavir attributable to the presence of PRO 88S and 82A + 54V. For the NRTIs lamivudine and emtricitabine, CRF01_AE had more discordances than subtype B due to the presence of RT mutational patterns 65R + 115 M and 118I + 215Y, respectively. Overall, the different algorithms agreed well on the level of resistance scored, but some of the discordances could be attributed to specific (subtype-dependent) combinations of mutations. It is not yet known whether therapy response is subtype dependent, but the advice given to clinicians based on a genotypic interpretation algorithm differs according to the subtype.

Molecular and clinical epidemiology of CXCR4-using HIV-1 in a large population of antiretroviral-naive individuals.

OBJECTIVE: We wished to characterize the epidemiological and clinical correlates of CXCR4-using human immunodeficiency virus type 1 (HIV-1) ("X4 variants") in a cross-sectional analysis of a large population of antiretroviral-naive individuals. METHODS: HIV-1 coreceptor use was determined in the last pretherapy plasma sample for 1191 individuals initiating triple-combination therapy in British Columbia, Canada. Baseline variables investigated included sociodemographic characteristics, plasma viral load (pVL), CD4 cell count, AIDS diagnosis, HIV-1 V3 loop sequence, and human CCR5 Delta 32 genotype. RESULTS: Individuals harboring X4 variants (n = 178 of 979 phenotyped samples; 18.2%) displayed a poorer baseline clinical profile than individuals harboring exclusively CCR5-using HIV-1 ("R5 variants") (median pVL, 175,000 vs. 120,000 copies of HIV-1 RNA/mL [P = .0006]; median CD4 cell count, 110 vs. 290 cells/mm(3) [P < .0001]). Individuals heterozygous for the CCR5 Delta 32 deletion (n = 128 of 967; 13.2%) were at 2.5 times higher risk of harboring X4 variants, compared with those without the deletion (multivariate P = .0005). The presence of basic amino acids at codon 11 and/or codon 25 of HIV-1 V3 (n = 109 of 955; 11.4%) was associated with a 9.1 times higher risk of harboring X4 variants (multivariate P < .0001), regardless of CCR5 Delta 32 genotype. In multivariate analyses adjusting for baseline parameters, HIV-1 coreceptor use was not found to be a significant predictor of survival or treatment response. CONCLUSION: Baseline CD4 cell count, pVL, HIV-1 V3 sequence, and CCR5 Delta 32 genotype were the strongest determinants of CXCR4-using HIV-1 in this population. After adjustment for baseline parameters, the presence of X4 variants before initiation of highly active antiretroviral therapy was not independently associated with a poorer outcome of therapy.

Rates of antiretroviral resistance among HIV-infected patients with and without a history of injection drug use.

BACKGROUND: There exist concerns regarding the potential for elevated rates of antiretroviral resistance among HIV-infected injection drug users (IDUs) prescribed highly active antiretroviral therapy (HAART), however, no population-based study has examined if IDUs have elevated rates of antiretroviral resistance in comparison to non-IDUs. OBJECTIVE: To evaluate the time to the development of antiretroviral resistance among antiretroviral-naive patients with and without a history of injection drug use. METHODS: In British Columbia there is a province-wide HIV/AIDS treatment program that provides antiretrovirals free of charge. We examined all antiretroviral-naive patients initiating HAART between 1 August 1996 and 30 September 2000 and who were followed to 31 March 2002. The main outcome measure was the time to class-specific antiretroviral resistance. Cumulative antiretroviral resistance rates among IDUs and non-IDUs were evaluated using Kaplan-Meier methods and relative hazards were estimated using Cox regression. RESULTS: Overall, 1191 antiretroviral-naive patients initiated HAART during the study period. Resistance mutations were observed in 298 (25%) subjects during the first 30 months of HAART. In comparison with non-IDUs, the risk of protease inhibitor resistance [relative hazard (RH), 0.9; 95% confidence interval (CI), 0.5-1.6] and non-nucleoside reverse transcriptase inhibitor resistance (RH, 1.5; 95% CI, 1.0-2.2) were similar among IDUs, and there were no differences in the rates of resistance to the sub-classes of nucleoside reverse transcriptase inhibitors. CONCLUSIONS: Resistance to all major classes of antiretrovirals were similar among IDUs and non-IDUs after 30 months of follow-up. These findings should help to allay fears that prescribing HAART to IDUs may result in elevated rates of resistance.

Rare mutations at codon 103 of HIV-1 reverse transcriptase can confer resistance to non-nucleoside reverse transcriptase inhibitors.

BACKGROUND: The K103N mutation in HIV-1 reverse transcriptase (RT) confers high-level resistance to current non-nucleoside reverse transcriptase inhibitors (NNRTI). The prevalence and resistance profile of HIV-1 with other substitutions at RT codon 103 is less well documented. METHODS: K103 substitutions among over 70,000 clinical samples submitted for routine antiretroviral resistance testing at two independent centres were examined. Phenotypic resistance profiles of isolates harboring rare K103 variants in the absence of known NNRTI-associated resistance mutations were retrieved from Virco's correlative genotype/phenotype database. Genotyped samples with known treatment histories were retrieved from the British Columbia Centre for Excellence in HIV/AIDS database. Site-directed mutants containing K103 variants were constructed and phenotyped. RESULTS: K103N, R and S were observed in 29, 1.8, and 0.9% of Virco isolates and in 16, 1.5 and 0.4% of British Columbia isolates. K103T/Q/H substitutions were observed only rarely (<0.2%). The prevalence of unusual codon 103 substitutions remained stable over 5 years, except K103S, which increased over fourfold in both datasets. K103R/Q-containing clinical isolates remained phenotypically susceptible to NNRTI, whereas K103S/T/H-containing isolates showed over 10-fold decreased NNRTI susceptibility. Among patients with a known treatment history, K103S/T/H were observed primarily in individuals failing NNRTI-containing regimens. Site-directed mutants confirmed decreased susceptibility to NNRTI in K103S/T/H-containing recombinants. CONCLUSION: Variants at HIV RT codon 103 other than K103N are observed relatively rarely in clinical isolates, but K103 S, T and H confer decreased susceptibility to NNRTI. These data are relevant for interpretive genotype algorithms and in the design of assays specific to RT codon 103 mutations.

Impact of HIV-1 subtype and antiretroviral therapy on protease and reverse transcriptase genotype: results of a global collaboration.

BACKGROUND: The genetic differences among HIV-1 subtypes may be critical to clinical management and drug resistance surveillance as antiretroviral treatment is expanded to regions of the world where diverse non-subtype-B viruses predominate. METHODS AND FINDINGS: To assess the impact of HIV-1 subtype and antiretroviral treatment on the distribution of mutations in protease and reverse transcriptase, a binomial response model using subtype and treatment as explanatory variables was used to analyze a large compiled dataset of non-subtype-B HIV-1 sequences. Non-subtype-B sequences from 3,686 persons with well characterized antiretroviral treatment histories were analyzed in comparison to subtype B sequences from 4,769 persons. The non-subtype-B sequences included 461 with subtype A, 1,185 with C, 331 with D, 245 with F, 293 with G, 513 with CRF01_AE, and 618 with CRF02_AG. Each of the 55 known subtype B drug-resistance mutations occurred in at least one non-B isolate, and 44 (80%) of these mutations were significantly associated with antiretroviral treatment in at least one non-B subtype. Conversely, of 67 mutations found to be associated with antiretroviral therapy in at least one non-B subtype, 61 were also associated with antiretroviral therapy in subtype B isolates. CONCLUSION: Global surveillance and genotypic assessment of drug resistance should focus primarily on the known subtype B drug-resistance mutations.

HIV-1 drug resistance: degree of underestimation by a cross-sectional versus a longitudinal testing approach.

Genotyping of human immunodeficiency virus type 1 (HIV-1) for antiretroviral drug resistance is routinely used both in clinical practice, to guide the selection of options for an individual's antiretroviral therapy, and in epidemiological studies, to estimate levels of antiretroviral drug resistance in a patient population. However, reliance on results of a single test can result in an underestimation of antiretroviral drug resistance. In the present study, we quantified the prevalence of resistance-associated mutations found in recent genotypic tests of 1734 HIV-1-infected, treatment-experienced subjects who had at least 3 genotypic tests (n = 11,404 genotypic tests total; median, 5 tests/subject) and compared it with that of resistance-associated mutations ever detected in these subjects between 1996 and 2004. Single-point analyses underestimated antiretroviral drug resistance, particularly for nucleoside analogues, in both individuals and patient populations. For example, the prevalence of resistance-associated mutation M184V/I was 25.5% in the most recent genotypes and 58.8% in available historical genotypes. Our results suggest that analysis of a combined historical genotype rather than of a cross-sectional genotype may lead to more accurate estimates of antiretroviral drug resistance in individual patients and in patient populations.

Predictors of HIV drug-resistance mutations in a large antiretroviral-naive cohort initiating triple antiretroviral therapy.

OBJECTIVE: The objective of this study was to systematically characterize the incidence and determinants of antiretroviral resistance in the HOMER (Highly Active Antiretroviral Therapy [HAART] Observational Medical Evaluation and Research) cohort of 1191 human immunodeficiency virus-infected, antiretroviral-naive adults initiating HAART in British Columbia, Canada. METHODS: All plasma samples with plasma virus loads (pVLs) >1000 copies/mL collected during the first 30 months of follow-up were genotyped for drug resistance. The primary outcome measure was time to the first detection of major drug-resistance mutation(s). Cox proportional hazard regression was used to identify factors significantly associated with the detection of drug-resistance mutations. RESULTS: Drug-resistance mutations were detected in 298 subjects (25%). Factors significantly associated with detection of drug-resistance mutations included high baseline pVL (multivariate hazard ratio [HR], 1.59; P<.001) and adherence (estimated using prescription-refill data and/or untimed plasma drug-concentration measurements). When compared with subjects with low (0%-<20%) prescription-refill percentages, subjects at an elevated risk of harboring drug-resistance mutations were those with relatively high but imperfect prescription-refill percentages (80%-<90%; multivariate HR, 4.15; P<.001) and those with essentially perfect (>/=95%) refill percentages but with 2 plasma drug concentrations below the steady-state trough concentration minus 1 standard deviation (multivariate HR, 4.57; P<.001). Initial use of nonnucleoside reverse-transcriptase inhibitor-based HAART was significantly associated with multiclass drug resistance (multivariate HR, 1.84; P=.001). CONCLUSION: High baseline pVLs and substantial but imperfect levels of adherence were major predictors of antiretroviral resistance.