The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages.

The independent emergence late in 2020 of the B.1.1.7, B.1.351, and P.1 lineages of SARS-CoV-2 prompted renewed concerns about the evolutionary capacity of this virus to overcome public health interventions and rising population immunity. Here, by examining patterns of synonymous and non-synonymous mutations that have accumulated in SARS-CoV-2 genomes since the pandemic began, we find that the emergence of these three "501Y lineages" coincided with a major global shift in the selective forces acting on various SARS-CoV-2 genes. Following their emergence, the adaptive evolution of 501Y lineage viruses has involved repeated selectively favored convergent mutations at 35 genome sites, mutations we refer to as the 501Y meta-signature. The ongoing convergence of viruses in many other lineages on this meta-signature suggests that it includes multiple mutation combinations capable of promoting the persistence of diverse SARS-CoV-2 lineages in the face of mounting host immune recognition.

Omicron extensively but incompletely escapes Pfizer BNT162b2 neutralization.

The emergence of the SARS-CoV-2 variant of concern Omicron (Pango lineage B.1.1.529), first identified in Botswana and South Africa, may compromise vaccine effectiveness and lead to re-infections1. Here we investigated Omicron escape from neutralization by antibodies from South African individuals vaccinated with Pfizer BNT162b2. We used blood samples taken soon after vaccination from individuals who were vaccinated and previously infected with SARS-CoV-2 or vaccinated with no evidence of previous infection. We isolated and sequence-confirmed live Omicron virus from an infected person and observed that Omicron requires the angiotensin-converting enzyme 2 (ACE2) receptor to infect cells. We compared plasma neutralization of Omicron relative to an ancestral SARS-CoV-2 strain and found that neutralization of ancestral virus was much higher in infected and vaccinated individuals compared with the vaccinated-only participants. However, both groups showed a 22-fold reduction in vaccine-elicited neutralization by the Omicron variant. Participants who were vaccinated and had previously been infected exhibited residual neutralization of Omicron similar to the level of neutralization of the ancestral virus observed in the vaccination-only group. These data support the notion that reasonable protection against Omicron may be maintained using vaccination approaches.

Omicron infection enhances Delta antibody immunity in vaccinated persons.

The extent to which Omicron infection1-9, with or without previous vaccination, elicits protection against the previously dominant Delta (B.1.617.2) variant is unclear. Here we measured the neutralization capacity against variants of severe acute respiratory syndrome coronavirus 2 in 39 individuals in South Africa infected with the Omicron sublineage BA.1 starting at a median of 6 (interquartile range 3-9) days post symptom onset and continuing until last follow-up sample available, a median of 23 (interquartile range 19-27) days post symptoms to allow BA.1-elicited neutralizing immunity time to develop. Fifteen participants were vaccinated with Pfizer's BNT162b2 or Johnson & Johnson's Ad26.CoV2.S and had BA.1 breakthrough infections, and 24 were unvaccinated. BA.1 neutralization increased from a geometric mean 50% focus reduction neutralization test titre of 42 at enrolment to 575 at the last follow-up time point (13.6-fold) in vaccinated participants and from 46 to 272 (6.0-fold) in unvaccinated participants. Delta virus neutralization also increased, from 192 to 1,091 (5.7-fold) in vaccinated participants and from 28 to 91 (3.0-fold) in unvaccinated participants. At the last time point, unvaccinated individuals infected with BA.1 had low absolute levels of neutralization for the non-BA.1 viruses and 2.2-fold lower BA.1 neutralization, 12.0-fold lower Delta neutralization, 9.6-fold lower Beta variant neutralization, 17.9-fold lower ancestral virus neutralization and 4.8-fold lower Omicron sublineage BA.2 neutralization relative to vaccinated individuals infected with BA.1. These results indicate that hybrid immunity formed by vaccination and Omicron BA.1 infection should be protective against Delta and other variants. By contrast, infection with Omicron BA.1 alone offers limited cross-protection despite moderate enhancement.

T cell responses to SARS-CoV-2 spike cross-recognize Omicron.

The SARS-CoV-2 Omicron variant (B.1.1.529) has multiple spike protein mutations1,2 that contribute to viral escape from antibody neutralization3-6 and reduce vaccine protection from infection7,8. The extent to which other components of the adaptive response such as T cells may still target Omicron and contribute to protection from severe outcomes is unknown. Here we assessed the ability of T cells to react to Omicron spike protein in participants who were vaccinated with Ad26.CoV2.S or BNT162b2, or unvaccinated convalescent COVID-19 patients (n = 70). Between 70% and 80% of the CD4+ and CD8+ T cell response to spike was maintained across study groups. Moreover, the magnitude of Omicron cross-reactive T cells was similar for Beta (B.1.351) and Delta (B.1.617.2) variants, despite Omicron harbouring considerably more mutations. In patients who were hospitalized with Omicron infections (n = 19), there were comparable T cell responses to ancestral spike, nucleocapsid and membrane proteins to those in patients hospitalized in previous waves dominated by the ancestral, Beta or Delta variants (n = 49). Thus, despite extensive mutations and reduced susceptibility to neutralizing antibodies of Omicron, the majority of T cell responses induced by vaccination or infection cross-recognize the variant. It remains to be determined whether well-preserved T cell immunity to Omicron contributes to protection from severe COVID-19 and is linked to early clinical observations from South Africa and elsewhere9-12.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa.

The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

Escape of SARS-CoV-2 501Y.V2 from neutralization by convalescent plasma.

SARS-CoV-2 variants of concern (VOC) have arisen independently at multiple locations1,2 and may reduce the efficacy of current vaccines that target the spike glycoprotein of SARS-CoV-23. Here, using a live-virus neutralization assay, we compared the neutralization of a non-VOC variant with the 501Y.V2 VOC (also known as B.1.351) using plasma collected from adults who were hospitalized with COVID-19 during the two waves of infection in South Africa, the second wave of which was dominated by infections with the 501Y.V2 variant. Sequencing demonstrated that infections of plasma donors from the first wave were with viruses that did not contain the mutations associated with 501Y.V2, except for one infection that contained the E484K substitution in the receptor-binding domain. The 501Y.V2 virus variant was effectively neutralized by plasma from individuals who were infected during the second wave. The first-wave virus variant was effectively neutralized by plasma from first-wave infections. However, the 501Y.V2 variant was poorly cross-neutralized by plasma from individuals with first-wave infections; the efficacy was reduced by 15.1-fold relative to neutralization of 501Y.V2 by plasma from individuals infected in the second wave. By contrast, cross-neutralization of first-wave virus variants using plasma from individuals with second-wave infections was more effective, showing only a 2.3-fold decrease relative to neutralization of first-wave virus variants by plasma from individuals infected in the first wave. Although we tested only one plasma sample from an individual infected with a SARS-CoV-2 variant with only the E484K substitution, this plasma sample potently neutralized both variants. The observed effective neutralization of first-wave virus by plasma from individuals infected with 501Y.V2 provides preliminary evidence that vaccines based on VOC sequences could retain activity against other circulating SARS-CoV-2 lineages.

Detection of a SARS-CoV-2 variant of concern in South Africa.

Continued uncontrolled transmission of SARS-CoV-2 in many parts of the world is creating conditions for substantial evolutionary changes to the virus1,2. Here we describe a newly arisen lineage of SARS-CoV-2 (designated 501Y.V2; also known as B.1.351 or 20H) that is defined by eight mutations in the spike protein, including three substitutions (K417N, E484K and N501Y) at residues in its receptor-binding domain that may have functional importance3-5. This lineage was identified in South Africa after the first wave of the epidemic in a severely affected metropolitan area (Nelson Mandela Bay) that is located on the coast of the Eastern Cape province. This lineage spread rapidly, and became dominant in Eastern Cape, Western Cape and KwaZulu-Natal provinces within weeks. Although the full import of the mutations is yet to be determined, the genomic data-which show rapid expansion and displacement of other lineages in several regions-suggest that this lineage is associated with a selection advantage that most plausibly results from increased transmissibility or immune escape6-8.

Population Effectiveness of Dolutegravir Implementation in Uganda: A Prospective Observational Cohort Study (DISCO), 48-Week Results.

BACKGROUND: Tenofovir/lamivudine/dolutegravir (TLD) is the preferred first-line antiretroviral therapy (ART) regimen for people with HIV (PWH), including those who were previously virologically suppressed on nonnucleoside reverse transcriptase inhibitors (NNRTIs). We sought to estimate the real-world effectiveness of the TLD transition in Ugandan public-sector clinics. METHODS: We conducted a prospective cohort study of PWH aged ≥18 years who were transitioned from NNRTI-based ART to TLD. Study visits were conducted on the day of TLD transition and 24 and 48 weeks later. The primary end point was viral suppression (<200 copies/mL) at 48 weeks. We collected blood for retrospective viral load (VL) assessment and conducted genotypic resistance tests for specimens with VL >500 copies/mL. RESULTS: We enrolled 500 participants (median age 47 years; 41% women). At 48 weeks after TLD transition, 94% of participants were in care with a VL <200 copies/mL (n = 469/500); 2% (n = 11/500) were lost from care or died; and only 2% (n = 9/500) had a VL >500 copies/mL. No incident resistance to DTG was identified. Few participants (2%, n = 9/500) discontinued TLD due to adverse events. CONCLUSIONS: High rates of viral suppression, high tolerability, and lack of emergent drug resistance support use of TLD as the preferred first-line regimen in the region. CLINICAL TRIALS REGISTRATION: NCT04066036.

Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function.

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Clinical Evaluation of Severe Acute Respiratory Syndrome Coronavirus 2 Rapid Antigen Tests During the Omicron Wave in South Africa.

We evaluated the performance of nasal and nasopharyngeal Standard Q COVID-19 [coronavirus disease 2019] Ag tests (SD Biosensor) and the Panbio COVID-19 Ag Rapid Test Device (nasal; Abbott) against the Abbott RealTime severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) assay during the Omicron (clades 21M, 21K, and 21L) wave in South Africa. Overall, all evaluated tests performed well, with high sensitivity (range, 77.78%-81.42%) and excellent specificity values (>99%). The sensitivity of rapid antigen tests increased above 90% in samples with cycle threshold <20, and all 3 tests performed best within the first week after symptom onset.

Comparison of SARS-CoV-2 sequencing using the ONT GridION and the Illumina MiSeq.

BACKGROUND: Over 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the SARS-CoV-2 consensus genomes obtained from the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs. RESULTS: Our results show that the MiSeq had a significantly higher number of consensus genomes classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher genome coverage and mutation counts than the GridION. CONCLUSION: Due to the low genome coverage, high number of indels, and sensitivity to SARS-CoV-2 viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for SARS-CoV-2 genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole consensus genomes.

Rapid Replacement of SARS-CoV-2 Variants by Delta and Subsequent Arrival of Omicron, Uganda, 2021.

Genomic surveillance in Uganda showed rapid replacement of severe acute respiratory syndrome coronavirus 2 over time by variants, dominated by Delta. However, detection of the more transmissible Omicron variant among travelers and increasing community transmission highlight the need for near-real-time genomic surveillance and adherence to infection control measures to prevent future pandemic waves.

Identification of SARS-CoV-2 Omicron variant using spike gene target failure and genotyping assays, Gauteng, South Africa, 2021.

The circulation of Omicron BA.1 led to the rapid increase in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cases in South Africa in November 2021, which warranted the use of more rapid detection methods. We, therefore, assessed the ability to detect Omicron BA.1 using genotyping assays to identify specific mutations in SARS-CoV-2 positive samples, Gauteng province, South Africa. The TaqPath™ COVID-19 real-time polymerase chain reaction assay was performed on all samples selected to identify spike gene target failure (SGTF). SARS-CoV-2 genotyping assays were used for the detection of del69/70 and K417N mutation. Whole-genome sequencing was performed on a subset of genotyped samples to confirm these findings. Of the positive samples received, 11.0% (175/1589) were randomly selected to assess if SGTF and genotyping assays, that detect del69/70 and K417N mutations, could identify Omicron BA.1. We identified SGTF in 98.9% (173/175) of samples, of which 88.0% (154/175) had both the del69/70 and K417N mutation. The genotyped samples (45.7%; 80/175) that were sequenced confirmed Omicron BA.1 (97.5%; 78/80). Our data show that genotyping for the detection of the del69/70 and K417N coupled with SGTF is efficient to exclude Alpha and Beta variants and rapidly detect Omicron BA.1. However, we still require assays for the detection of unique mutations that will allow for the differentiation between other Omicron sublineages. Therefore, the use of genotyping assays to detect new dominant or emerging lineages of SARS-CoV-2 will be beneficial in limited-resource settings.

Acquired HIV drug resistance and virologic monitoring in a HIV hyper-endemic setting in KwaZulu-Natal Province, South Africa.

BACKGROUND: Introduction of tenofovir (TDF) plus lamivudine (3TC) and dolutegravir (DTG) in first- and second-line HIV treatment regimens in South Africa warrants characterization of acquired HIV-1 drug resistance (ADR) mutations that could impact DTG-based antiretroviral therapy (ART). In this study, we sought to determine prevalence of ADR mutations and their potential impact on susceptibility to drugs used in combination with DTG among HIV-positive adults (≥ 18 years) accessing routine care at a selected ART facility in KwaZulu-Natal, South Africa. METHODS: We enrolled adult participants in a cross-sectional study between May and September 2019. Eligible participants had a most recent documented viral load (VL) ≥ 1000 copies/mL after at least 6 months on ART. We genotyped HIV-1 reverse transcriptase and protease genes by Sanger sequencing and assessed ADR. We characterized the effect of ADR mutations on the predicted susceptibility to drugs used in combination with DTG. RESULTS: From 143 participants enrolled, we obtained sequence data for 115 (80%), and 92.2% (95% CI 85.7-96.4) had ADR. The proportion with ADR was similar for participants on first-line ART (65/70, 92.9%, 95% CI 84.1-97.6) and those on second-line ART (40/44, 90.9%, 95% CI 78.3-97.5), and was present for the single participant on third-line ART. Approximately 89% (62/70) of those on first-line ART had dual class NRTI and NNRTI resistance and only six (13.6%) of those on second-line ART had major PI mutations. Most participants (82%) with first-line viraemia maintained susceptibility to Zidovudine (AZT), and the majority of them had lost susceptibility to TDF (71%) and 3TC (84%). Approximately two in every five TDF-treated individuals had thymidine analogue mutations (TAMs). CONCLUSIONS: Susceptibility to AZT among most participants with first-line viraemia suggests that a new second-line regimen of AZT + 3TC + DTG could be effective. However, atypical occurrence of TAMs in TDF-treated individuals suggests a less effective AZT + 3TC + DTG regimen in a subpopulation of patients. As most patients with first-line viraemia had at least low-level resistance to TDF and 3TC, identifying viraemia before switch to TDF + 3TC + DTG is important to avoid DTG functional monotherapy. These findings highlight a need for close monitoring of outcomes on new standardized treatment regimens.

Unlocking the efficiency of genomics laboratories with robotic liquid-handling.

In research and clinical genomics laboratories today, sample preparation is the bottleneck of experiments, particularly when it comes to high-throughput next generation sequencing (NGS). More genomics laboratories are now considering liquid-handling automation to make the sequencing workflow more efficient and cost effective. The question remains as to its suitability and return on investment. A number of points need to be carefully considered before introducing robots into biological laboratories. Here, we describe the state-of-the-art technology of both sophisticated and do-it-yourself (DIY) robotic liquid-handlers and provide a practical review of the motivation, implications and requirements of laboratory automation for genome sequencing experiments.

Impact of pretreatment low-abundance HIV-1 drug-resistant variants on virological failure among HIV-1/TB-co-infected individuals.

OBJECTIVES: To determine the impact of pretreatment low-abundance HIV-1 drug-resistant variants (LA-DRVs) on virological failure (VF) among HIV-1/TB-co-infected individuals treated with NNRTI first-line ART. METHODS: We conducted a case-control study of 170 adults with HIV-1/TB co-infection. Cases had at least one viral load (VL) ≥1000 RNA copies/mL after ≥6 months on NNRTI-based ART, and controls had sustained VLs <1000 copies/mL. We sequenced plasma viruses by Sanger and MiSeq next-generation sequencing (NGS). We assessed drug resistance mutations (DRMs) using the Stanford drug resistance database, and analysed NGS data for DRMs at ≥20%, 10%, 5% and 2% thresholds. We assessed the effect of pretreatment drug resistance (PDR) on VF. RESULTS: We analysed sequences from 45 cases and 125 controls. Overall prevalence of PDR detected at a ≥20% threshold was 4.7% (8/170) and was higher in cases than in controls (8.9% versus 3.2%), P = 0.210. Participants with PDR at ≥20% had almost 4-fold higher odds of VF (adjusted OR 3.7, 95% CI 0.8-18.3) compared with those without, P = 0.104. PDR prevalence increased to 18.2% (31/170) when LA-DRVs at ≥2% were included. Participants with pretreatment LA-DRVs only had 1.6-fold higher odds of VF (adjusted OR 1.6, 95% CI 0.6-4.3) compared with those without, P = 0.398. CONCLUSIONS: Pretreatment DRMs and LA-DRVs increased the odds of developing VF on NNRTI-based ART, although without statistical significance. NGS increased detection of DRMs but provided no additional benefit in identifying participants at risk of VF at lower thresholds. More studies assessing mutation thresholds predictive of VF are required to inform use of NGS in treatment decisions.

Predicted antiviral activity of tenofovir versus abacavir in combination with a cytosine analogue and the integrase inhibitor dolutegravir in HIV-1-infected South African patients initiating or failing first-line ART.

Objectives: The WHO recently recommended the use of a new first-line ART containing dolutegravir. We investigated the efficacy of NRTI backbones (tenofovir or abacavir with a cytosine analogue) in low- and middle-income countries where there is significant prior exposure to antiretrovirals and drug resistance to NRTIs. Methods: Within the treatment-as-prevention study in South Africa, we selected participants with available next-generation sequencing (NGS) data for the HIV-1 pol gene at trial entry; they were either ART initiators (n = 1193) or already established on ART (n = 94). NGS of the HIV-1 pol gene was carried out using MiSeq technology; reverse transcriptase drug resistance mutations (DRMs) were detected at 5% (DRM5%) and 20% (DRM20%) for all 1287 participants. Genotypic susceptibility was assessed using the Stanford HIVDB resistance interpretation algorithm. Results: NRTI DRM20% and DRM5% were detected among 5/1193 (0.4%) and 9/1193 (0.8%) of ART initiators, respectively. There was tenofovir exposure in 73/94 (77.7%) of those established on ART, with full susceptibility to abacavir in 57/94 (60.6%) and 56/94 (59.6%) for DRM20% and DRM5%, respectively, while 67/94 (71.3%) and 64/94 (68.1%) were fully susceptible to tenofovir, respectively. The differences between tenofovir and abacavir were not statistically significant at the 20% or 5% variant level (P = 0.16 and 0.29, respectively). NGS detection of variants at the 5% level increased detection of K65R in both naive and treated groups. One of 607 integrase sequences carried a DRM20% (Q148R). Conclusions: Dolutegravir with a cytosine analogue plus tenofovir or abacavir appears to have similar efficacy in South Africans naive to ART. NGS should be considered in HIV drug resistance surveillance.

High rate of occult hepatitis B virus infection in hemodialysis units of KwaZulu-Natal, South Africa.

Occult hepatitis B virus (HBV) infection (OBI) is defined as the presence of HBV DNA in the liver with or without detectable HBV DNA in the serum of individuals testing HBV surface antigen (HBsAg) negative using currently available assays. The prevalence of OBI among patients receiving hemodialysis (HD) treatment remains poorly characterized in South Africa despite the high prevalence of HBV. We sought to determine the prevalence of OBI in HD units in tertiary hospitals of KwaZulu-Natal and to characterize the HBV S gene mutations potentially responsible for OBI. A cross-sectional descriptive study of residual diagnostic plasma samples from 85 HBsAg-negative patients receiving HD treatment was included. The PreS/S gene was amplified with a nested HBV polymerase chain reaction for downstream next-generation sequencing, to determine the viral genotype and identify S gene mutations associated with OBI. Nine of the 85 samples had OBI, based on detectable HBV DNA. The point prevalence of OBI was 10.6% (95% control interval: 5.5%-19.1%). Phylogenetic analysis of the samples with OBI showed that all belonged to genotype A. Three (~33%) samples had mutations in the major hydrophilic region (MHR) within the S gene, three (~33%) had mutations within the S gene but outside the MHR, whilst the remaining three had no mutations observed. The prevalence of OBI in HBsAg-negative patients undergoing HD was 10.6%, suggesting that OBI is a clinically significant problem in patients with HD in this region. The screening methods for HBV infection need to be revised to include nucleic acid testing.

Contribution of Gag and Protease to HIV-1 Phenotypic Drug Resistance in Pediatric Patients Failing Protease Inhibitor-Based Therapy.

Protease inhibitors (PIs) are used as a first-line regimen in HIV-1-infected children. Here we investigated the phenotypic consequences of amino acid changes in Gag and protease on lopinavir (LPV) and ritonavir (RTV) susceptibility among pediatric patients failing PI therapy. The Gag-protease from isolates from 20 HIV-1 subtype C-infected pediatric patients failing an LPV and/or RTV-based regimen was phenotyped using a nonreplicativein vitroassay. Changes in sensitivity to LPV and RTV relative to that of the matched baseline (pretherapy) sample were calculated. Gag and protease amino acid substitutions associated with PI failure were created in a reference clone by site-directed mutagenesis and assessed. Predicted phenotypes were determined using the Stanford drug resistance algorithm. Phenotypic resistance or reduced susceptibility to RTV and/or LPV was observed in isolates from 10 (50%) patients, all of whom had been treated with RTV. In most cases, this was associated with protease resistance mutations, but substitutions at Gag cleavage and noncleavage sites were also detected. Gag amino acid substitutions were also found in isolates from three patients with reduced drug susceptibilities who had wild-type protease. Site-directed mutagenesis confirmed that some amino acid changes in Gag contributed to PI resistance but only in the presence of major protease resistance-associated substitutions. The isolates from all patients who received LPV exclusively were phenotypically susceptible. Baseline isolates from the 20 patients showed a large (47-fold) range in the 50% effective concentration of LPV, which accounted for most of the discordance seen between the experimentally determined and the predicted phenotypes. Overall, the inclusion of thegaggene and the use of matched baseline samples provided a more comprehensive assessment of the effect of PI-induced amino acid changes on PI resistance. The lack of phenotypic resistance to LPV supports the continued use of this drug in pediatric patients.