GPT-4 performance on querying scientific publications: reproducibility, accuracy, and impact of an instruction sheet.

BACKGROUND: Large language models (LLMs) that can efficiently screen and identify studies meeting specific criteria would streamline literature reviews. Additionally, those capable of extracting data from publications would enhance knowledge discovery by reducing the burden on human reviewers. METHODS: We created an automated pipeline utilizing OpenAI GPT-4 32 K API version "2023-05-15" to evaluate the accuracy of the LLM GPT-4 responses to queries about published papers on HIV drug resistance (HIVDR) with and without an instruction sheet. The instruction sheet contained specialized knowledge designed to assist a person trying to answer questions about an HIVDR paper. We designed 60 questions pertaining to HIVDR and created markdown versions of 60 published HIVDR papers in PubMed. We presented the 60 papers to GPT-4 in four configurations: (1) all 60 questions simultaneously; (2) all 60 questions simultaneously with the instruction sheet; (3) each of the 60 questions individually; and (4) each of the 60 questions individually with the instruction sheet. RESULTS: GPT-4 achieved a mean accuracy of 86.9% - 24.0% higher than when the answers to papers were permuted. The overall recall and precision were 72.5% and 87.4%, respectively. The standard deviation of three replicates for the 60 questions ranged from 0 to 5.3% with a median of 1.2%. The instruction sheet did not significantly increase GPT-4's accuracy, recall, or precision. GPT-4 was more likely to provide false positive answers when the 60 questions were submitted individually compared to when they were submitted together. CONCLUSIONS: GPT-4 reproducibly answered 3600 questions about 60 papers on HIVDR with moderately high accuracy, recall, and precision. The instruction sheet's failure to improve these metrics suggests that more sophisticated approaches are necessary. Either enhanced prompt engineering or finetuning an open-source model could further improve an LLM's ability to answer questions about highly specialized HIVDR papers.

HIV-1 5'-Leader Mutations in Plasma Viruses Before and After the Development of Reverse Transcriptase Inhibitor-Resistance Mutations.

Background: HIV-1 RT initiation depends on interaction between viral 5'-leader RNA, RT, and host tRNA3Lys. We therefore sought to identify co-evolutionary changes between the 5'-leader and RT in viruses developing RT-inhibitor resistance mutations. Methods: We sequenced 5'-leader positions 37-356 of paired plasma virus samples from 29 individuals developing the NRTI-resistance mutation M184V, 19 developing an NNRTI-resistance mutation, and 32 untreated controls. 5'-leader variants were defined as positions where ≥20% of NGS reads differed from the HXB2 sequence. Emergent mutations were defined as nucleotides undergoing ≥4-fold change in proportion between baseline and follow-up. Mixtures were defined as positions containing ≥2 nucleotides each present in ≥20% of NGS reads. Results: Among 80 baseline sequences, 87 positions (27.2%) contained a variant; 52 contained a mixture. Position 201 was the only position more likely to develop a mutation in the M184V (9/29 vs. 0/32; p=0.0006) or NNRTI-resistance (4/19 vs. 0/32; p=0.02; Fisher's Exact Test) groups than the control group. Mixtures at positions 200 and 201 occurred in 45.0% and 28.8%, respectively, of baseline samples. Because of the high proportion of mixtures at these positions, we analyzed 5'-leader mixture frequencies in two additional datasets: five publications reporting 294 dideoxyterminator clonal GenBank sequences from 42 individuals and six NCBI BioProjects reporting NGS datasets from 295 individuals. These analyses demonstrated position 200 and 201 mixtures at proportions similar to those in our samples and at frequencies several times higher than at all other 5'-leader positions. Conclusions: Although we did not convincingly document co-evolutionary changes between RT and 5'-leader sequences, we identified a novel phenomenon, wherein positions 200 and 201, immediately downstream of the HIV-1 primer binding site exhibited an extraordinarily high likelihood of containing a nucleotide mixture. Possible explanations for the high mixture rates are that these positions are particularly error-prone or provide a viral fitness advantage.

HIV-1 Group M Capsid Amino Acid Variability: Implications for Sequence Quality Control of Genotypic Resistance Testing.

BACKGROUND: With the approval of the HIV-1 capsid inhibitor, lenacapavir, capsid sequencing will be required for managing lenacapavir-experienced individuals with detectable viremia. Successful sequence interpretation will require examining new capsid sequences in the context of previously published sequence data. METHODS: We analyzed published HIV-1 group M capsid sequences from 21,012 capsid-inhibitor naïve individuals to characterize amino acid variability at each position and influence of subtype and cytotoxic T lymphocyte (CTL) selection pressure. We determined the distributions of usual mutations, defined as amino acid differences from the group M consensus, with a prevalence ≥ 0.1%. Co-evolving mutations were identified using a phylogenetically-informed Bayesian graphical model method. RESULTS: 162 (70.1%) positions had no usual mutations (45.9%) or only conservative usual mutations with a positive BLOSUM62 score (24.2%). Variability correlated independently with subtype-specific amino acid occurrence (Spearman rho = 0.83; p < 1 × 10-9) and the number of times positions were reported to contain an HLA-associated polymorphism, an indicator of CTL pressure (rho = 0.43; p = 0.0002). CONCLUSIONS: Knowing the distribution of usual capsid mutations is essential for sequence quality control. Comparing capsid sequences from lenacapavir-treated and lenacapavir-naïve individuals will enable the identification of additional mutations potentially associated with lenacapavir therapy.

Highly Ambiguous HIV-1 Pol Positions Encoding Multiple Amino Acids Usually Result from Antiviral or Immune Selection Pressure.

HIV-1 pol nucleotide ambiguities encoding amino acid mixtures occur commonly during population-based genotypic drug resistance testing. However, few studies have addressed the validity of sequences with fully ambiguous codons (FACs) containing codons translatable to more than four amino acids. We identified 839 published HIV-1 pol sequences with 846 FACs at 131 positions and determined their distribution relative to 215 HLA-associated pol positions (HAPs) and 84 drug-resistance positions. Among HIV-1 reverse transcriptase (RT) and protease sequences from antiretroviral therapy (ART)-naive and -experienced persons, there was a strong correlation between the likelihood a position was a FAC and that it was an HAP (Spearman's correlation coefficient rho >0.40; p < 1e-6). Among HIV-1 RT sequences from ART-experienced persons, there was a correlation between the likelihood that a position was a FAC and that it was a drug-resistance position (rho = 0.2; p = 8e-4). In the context of population-based genotypic resistance testing, FACs usually result from antiviral or immune selection pressure.

Sierra SARS-CoV-2 sequence and antiviral resistance analysis program.

INTRODUCTION: Although most laboratories are capable of employing established protocols to perform full-genome SARS-CoV-2 sequencing, many are unable to assess sequence quality, select appropriate mutation-detection thresholds, or report on the potential clinical significance of mutations in the targets of antiviral therapy METHODS: We describe the technical aspects and benchmark the performance of Sierra SARS-CoV-2, a program designed to perform these functions on user-submitted FASTQ and FASTA sequence files and lists of Spike mutations. Sierra SARS-CoV-2 indicates which sequences contain an unexpectedly large number of unusual mutations and which mutations are associated with reduced susceptibility to clinical stage mAbs, the RdRP inhibitor remdesivir, or the Mpro inhibitor nirmatrelvir RESULTS: To assess the performance of Sierra SARS-CoV-2 on FASTQ files, we applied it to 600 representative FASTQ sequences and compared the results to the COVID-19 EDGE program. To assess its performance on FASTA files, we applied it to nearly one million representative FASTA sequences and compared the results to the GISAID mutation annotation. To assess its performance on mutations lists, we applied it to 13,578 distinct Spike RBD mutation patterns and showed that exactly or partially matching annotations were available for 88% of patterns CONCLUSION: Sierra SARS-CoV-2 leverages previously published data to improve the quality control of submitted viral genomic data and to provide functional annotation on the impact of mutations in the targets of antiviral SARS-CoV-2 therapy. The program can be found at https://covdb.stanford.edu/sierra/sars2/ and its source code at https://github.com/hivdb/sierra-sars2.

Susceptibility of SARS-CoV-2 Omicron Variants to Therapeutic Monoclonal Antibodies: Systematic Review and Meta-analysis.

SARS-CoV-2 Omicron variants contain many mutations in its spike receptor-binding domain, the target of all authorized monoclonal antibodies (MAbs). Determining the extent to which Omicron variants reduced MAb susceptibility is critical to preventing and treating COVID-19. We systematically reviewed PubMed and three preprint servers, last updated 11 April 2022, for the in vitro activity of authorized MAbs against the Omicron variants. Fifty-one studies were eligible, including 50 containing Omicron BA.1 susceptibility data and 17 containing Omicron BA.2 susceptibility data. The first two authorized MAb combinations, bamlanivimab/etesevimab and casirivimab/imdevimab, were largely inactive against the Omicron BA.1 and BA.2 variants. In 34 studies, sotrovimab displayed a median 4.0-fold (interquartile range [IQR]: 2.6 to 6.9) reduction in activity against Omicron BA.1, and in 12 studies, it displayed a median 17-fold (IQR: 13 to 30) reduction in activity against Omicron BA.2. In 15 studies, the combination cilgavimab/tixagevimab displayed a median 86-fold (IQR: 27 to 151) reduction in activity against Omicron BA.1, and in six studies, it displayed a median 5.4-fold (IQR: 3.7 to 6.9) reduction in activity against Omicron BA.2. In eight studies against Omicron BA.1 and six studies against Omicron BA.2, bebtelovimab displayed no reduction in activity. Disparate results between assays were common. For authorized MAbs, 51/268 (19.0%) results for wild-type control variants and 78/348 (22.4%) results for Omicron BA.1 and BA.2 variants were more than 4-fold below or 4-fold above the median result for that MAb. Highly disparate results between published assays indicate a need for improved MAb susceptibility test standardization or interassay calibration. IMPORTANCE Monoclonal antibodies (MAbs) targeting the SARS-CoV-2 spike protein are among the most effective measures for preventing and treating COVID-19. However, SARS-CoV-2 Omicron variants contain many mutations in their spike receptor-binding domains, the target of all authorized MAbs. Therefore, determining the extent to which Omicron variants reduced MAb susceptibility is critical to preventing and treating COVID-19. We identified 51 studies that reported the in vitro susceptibility of the two main Omicron variants BA.1 and BA.2 to therapeutic MAbs in advanced clinical development, including eight authorized individual MAbs and three authorized MAb combinations. We estimated the degree to which different MAbs displayed reduced activity against Omicron variants. The marked loss of activity of many MAbs against Omicron variants underscores the importance of developing MAbs that target conserved regions of spike. Highly disparate results between assays indicate the need for improved MAb susceptibility test standardization.

Coronavirus Resistance Database (CoV-RDB): SARS-CoV-2 susceptibility to monoclonal antibodies, convalescent plasma, and plasma from vaccinated persons.

As novel SARS-CoV-2 variants with different patterns of spike protein mutations have emerged, the susceptibility of these variants to neutralization by antibodies has been rapidly assessed. However, neutralization data are generated using different approaches and are scattered across different publications making it difficult for these data to be located and synthesized. The Stanford Coronavirus Resistance Database (CoV-RDB; https://covdb.stanford.edu) is designed to house comprehensively curated published data on the neutralizing susceptibility of SARS-CoV-2 variants and spike mutations to monoclonal antibodies (mAbs), convalescent plasma (CP), and vaccinee plasma (VP). As of December 31, 2021, CoV-RDB encompassed 257 publications including 91 (35%) containing 9,070 neutralizing mAb susceptibility results, 131 (51%) containing 16,773 neutralizing CP susceptibility results, and 178 (69%) containing 33,540 neutralizing VP results. The database also records which spike mutations are selected during in vitro passage of SARS-CoV-2 in the presence of mAbs and which emerge in persons receiving mAbs as treatment. The CoV-RDB interface interactively displays neutralizing susceptibility data at different levels of granularity by filtering and/or aggregating query results according to one or more experimental conditions. The CoV-RDB website provides a companion sequence analysis program that outputs information about mutations present in a submitted sequence and that also assists users in determining the appropriate mutation-detection thresholds for identifying non-consensus amino acids. The most recent data underlying the CoV-RDB can be downloaded in its entirety from a GitHub repository in a documented machine-readable format.

The biological and clinical significance of emerging SARS-CoV-2 variants.

The past several months have witnessed the emergence of SARS-CoV-2 variants with novel spike protein mutations that are influencing the epidemiological and clinical aspects of the COVID-19 pandemic. These variants can increase rates of virus transmission and/or increase the risk of reinfection and reduce the protection afforded by neutralizing monoclonal antibodies and vaccination. These variants can therefore enable SARS-CoV-2 to continue its spread in the face of rising population immunity while maintaining or increasing its replication fitness. The identification of four rapidly expanding virus lineages since December 2020, designated variants of concern, has ushered in a new stage of the pandemic. The four variants of concern, the Alpha variant (originally identified in the UK), the Beta variant (originally identified in South Africa), the Gamma variant (originally identified in Brazil) and the Delta variant (originally identified in India), share several mutations with one another as well as with an increasing number of other recently identified SARS-CoV-2 variants. Collectively, these SARS-CoV-2 variants complicate the COVID-19 research agenda and necessitate additional avenues of laboratory, epidemiological and clinical research.

SARS-CoV-2 Antiviral Therapy.

The development of effective antiviral therapy for COVID-19 is critical for those awaiting vaccination, as well as for those who do not respond robustly to vaccination. This review summarizes 1 year of progress in the race to develop antiviral therapies for COVID-19, including research spanning preclinical and clinical drug development efforts, with an emphasis on antiviral compounds that are in clinical development or that are high priorities for clinical development. The review is divided into sections on compounds that inhibit SARS-CoV-2 enzymes, including its polymerase and proteases; compounds that inhibit virus entry, including monoclonal antibodies; interferons; and repurposed drugs that inhibit host processes required for SARS-CoV-2 replication. The review concludes with a summary of the lessons to be learned from SARS-CoV-2 drug development efforts and the challenges to continued progress.

Temporal Trends in HIV-1 Mutations Used for the Surveillance of Transmitted Drug Resistance.

In 2009, a list of nonpolymorphic HIV-1 drug resistance mutations (DRMs), called surveillance DRMs (SDRMs), was created to monitor transmitted drug resistance (TDR). Since 2009, TDR increased and antiretroviral therapy (ART) practices changed. We examined the changing prevalence of SDRMs and identified candidate SDRMs defined as nonpolymorphic DRMs present on ≥ 1 expert DRM list and in ≥0.1% of ART-experienced persons. Candidate DRMs were further characterized according to their association with antiretrovirals and changing prevalence. Among NRTI-SDRMs, tenofovir-associated mutations increased in prevalence while thymidine analog mutations decreased in prevalence. Among candidate NRTI-SDRMs, there were six tenofovir-associated mutations including three which increased in prevalence (K65N, T69deletion, K70G/N/Q/T). Among candidate NNRTI-SDRMs, six that increased in prevalence were associated with rilpivirine (E138K/Q, V179L, H221Y) or doravirine (F227C/L) resistance. With the notable exceptions of I47A and I50L, most PI-SDRMs decreased in prevalence. Three candidate PI-SDRMs were accessory darunavir-resistance mutations (L10F, T74P, L89V). Adding the candidate SDRMs listed above was estimated to increase NRTI, NNRTI, and PI TDR prevalence by 0.1%, 0.3%, and 0.3%, respectively. We describe trends in the prevalence of nonpolymorphic HIV-1 DRMs in ART-experienced persons. These data should be considered in decisions regarding SDRM list updates and TDR monitoring.

Coronavirus Antiviral Research Database (CoV-RDB): An Online Database Designed to Facilitate Comparisons between Candidate Anti-Coronavirus Compounds.

BACKGROUND: To prioritize the development of antiviral compounds, it is necessary to compare their relative preclinical activity and clinical efficacy. METHODS: We reviewed in vitro, animal model, and clinical studies of candidate anti-coronavirus compounds and placed extracted data in an online relational database. RESULTS: As of August 2020, the Coronavirus Antiviral Research Database (CoV-RDB; covdb.stanford.edu) contained over 2800 cell culture, entry assay, and biochemical experiments, 259 animal model studies, and 73 clinical studies from over 400 published papers. SARS-CoV-2, SARS-CoV, and MERS-CoV account for 85% of the data. Approximately 75% of experiments involved compounds with known or likely mechanisms of action, including monoclonal antibodies and receptor binding inhibitors (21%), viral protease inhibitors (17%), miscellaneous host-acting inhibitors (10%), polymerase inhibitors (9%), interferons (7%), fusion inhibitors (5%), and host protease inhibitors (5%). Of 975 compounds with known or likely mechanism, 135 (14%) are licensed in the U.S. for other indications, 197 (20%) are licensed outside the U.S. or are in human trials, and 595 (61%) are pre-clinical investigational compounds. CONCLUSION: CoV-RDB facilitates comparisons between different candidate antiviral compounds, thereby helping scientists, clinical investigators, public health officials, and funding agencies prioritize the most promising compounds and repurposed drugs for further development.

Analysis of unusual and signature APOBEC-mutations in HIV-1 pol next-generation sequences.

INTRODUCTION: At low mutation-detection thresholds, next generation sequencing (NGS) for HIV-1 genotypic resistance testing is susceptible to artifactual detection of mutations arising from PCR error and APOBEC-mediated G-to-A hypermutation. METHODS: We analyzed published HIV-1 pol Illumina NGS data to characterize the distribution of mutations at eight NGS mutation detection thresholds: 20%, 10%, 5%, 2%, 1%, 0.5%, 0.2%, and 0.1%. At each threshold, we determined proportions of amino acid mutations that were unusual (defined as having a prevalence <0.01% in HIV-1 group M sequences) or signature APOBEC mutations. RESULTS: Eight studies, containing 855 samples, in the NCBI Sequence Read Archive were analyzed. As detection thresholds were lowered, there was a progressive increase in the proportion of positions with usual and unusual mutations and in the proportion of all mutations that were unusual. The median proportion of positions with an unusual mutation increased gradually from 0% at the 20% threshold to 0.3% at the 1% threshold and then exponentially to 1.3% (0.5% threshold), 6.9% (0.2% threshold), and 23.2% (0.1% threshold). In two of three studies with available plasma HIV-1 RNA levels, the proportion of positions with unusual mutations was negatively associated with virus levels. Although the complete set of signature APOBEC mutations was much smaller than that of unusual mutations, the former outnumbered the latter in one-sixth of samples at the 0.5%, 1%, and 2% thresholds. CONCLUSIONS: The marked increase in the proportion of positions with unusual mutations at thresholds below 1% and in samples with lower virus loads suggests that, at low thresholds, many unusual mutations are artifactual, reflecting PCR error or G-to-A hypermutation. Profiling the numbers of unusual and signature APOBEC pol mutations at different NGS mutation detection thresholds may be useful to avoid selecting a threshold that is too low and poses an unacceptable risk of identifying artifactual mutations.

Expanded Spectrum of Antiretroviral-Selected Mutations in Human Immunodeficiency Virus Type 2.

BACKGROUND: HIV-1 and HIV-2 differ in their antiretroviral (ARV) susceptibilities and drug resistance mutations (DRMs). METHODS: We analyzed published HIV-2 pol sequences to identify HIV-2 treatment-selected mutations (TSMs). Mutation prevalences were determined by HIV-2 group and ARV status. Nonpolymorphic mutations were those in <1% of ARV-naive persons. TSMs were those associated with ARV therapy after multiple comparisons adjustment. RESULTS: We analyzed protease (PR) sequences from 483 PR inhibitor (PI)-naive and 232 PI-treated persons; RT sequences from 333 nucleoside RT inhibitor (NRTI)-naive and 252 NRTI-treated persons; and integrase (IN) sequences from 236 IN inhibitor (INSTI)-naive and 60 INSTI-treated persons. In PR, 12 nonpolymorphic TSMs occurred in ≥11 persons: V33I, K45R, V47A, I50V, I54M, T56V, V62A, A73G, I82F, I84V, F85L, L90M. In RT, 9 nonpolymorphic TSMs occurred in ≥10 persons: K40R, A62V, K70R, Y115F, Q151M, M184VI, S215Y. In IN, 11 nonpolymorphic TSMs occurred in ≥4 persons: Q91R, E92AQ, T97A, G140S, Y143G, Q148R, A153G, N155H, H156R, R231 5-amino acid insertions. Nine of 32 nonpolymorphic TSMs were previously unreported. CONCLUSIONS: This meta-analysis confirmed the ARV association of previously reported HIV-2 DRMs and identified novel TSMs. Genotypic and phenotypic studies of HIV-2 TSMs will improve approaches to predicting HIV-2 ARV susceptibility and treating HIV-2-infected persons.

Integrase strand transfer inhibitor (INSTI)-resistance mutations for the surveillance of transmitted HIV-1 drug resistance.

BACKGROUND: Integrase strand transfer inhibitors (INSTIs) are expected to be widely adopted globally, requiring surveillance of resistance emergence and transmission. OBJECTIVES: We therefore sought to develop a standardized list of INSTI-resistance mutations suitable for the surveillance of transmitted INSTI resistance. METHODS: To characterize the suitability of the INSTI-resistance mutations for transmitted HIV-1 drug resistance (TDR) surveillance, we classified them according to their presence on published expert lists, conservation in INSTI-naive persons, frequency in INSTI-treated persons and contribution to reduced in vitro susceptibility. Mutation prevalences were determined using integrase sequences from 17302 INSTI-naive and 2450 INSTI-treated persons; 53.3% of the INSTI-naive sequences and 20.0% of INSTI-treated sequences were from non-B subtypes. Approximately 10% of sequences were from persons who received dolutegravir alone or a first-generation INSTI followed by dolutegravir. RESULTS: Fifty-nine previously recognized (or established) INSTI-resistance mutations were present on one or more of four published expert lists. They were classified into three main non-overlapping groups: 29 relatively common non-polymorphic mutations, occurring in five or more individuals and significantly selected by INSTI treatment; 8 polymorphic mutations; and 22 rare mutations. Among the 29 relatively common INSTI-selected mutations, 24 emerged as candidates for inclusion on a list of INSTI surveillance drug-resistance mutations: T66A/I/K, E92G/Q, G118R, F121Y, E138A/K/T, G140A/C/S, Y143C/H/R/S, S147G, Q148H/R/K, N155H, S230R and R263K. CONCLUSIONS: A set of 24 non-polymorphic INSTI-selected mutations is likely to be useful for quantifying INSTI-associated TDR. This list may require updating as more sequences become available from INSTI-experienced persons infected with HIV-1 non-subtype B viruses and/or receiving dolutegravir.

Prospective Evaluation of the Vela Diagnostics Next-Generation Sequencing Platform for HIV-1 Genotypic Resistance Testing.

Genotypic antiretroviral drug resistance testing is a critical component of the global efforts to control the HIV-1 epidemic. This study investigates the semiautomated, next-generation sequencing (NGS)-based Vela Diagnostics Sentosa SQ HIV-1 Genotyping Assay in a prospective cohort of HIV-1-infected patients. Two-hundred sixty-nine samples were successfully sequenced by both NGS and Sanger sequencing. Among the 261 protease/reverse transcriptase (PR/RT) sequences, a mean of 0.37 drug resistance mutations were identified by both Sanger and NGS, 0.08 by NGS alone, and 0.03 by Sanger alone. Among the 50 integrase sequences, a mean of 0.3 drug resistance mutations were detected by both Sanger and NGS, and 0.08 by NGS alone. NGS estimated higher levels of drug resistance to one or more antiretroviral drugs for 6.5% of PR/RT sequences and 4.0% of integrase sequences, whereas Sanger estimated higher levels of drug resistance for 3.8% of PR/RT sequences. Although the samples successfully sequenced by the Sentosa SQ HIV Genotyping Assay demonstrated similar predicted resistance compared with Sanger, 44% of Sentosa runs failed quality control requiring 17 additional runs. This semi-automated NGS-based assay may aid in HIV-1 genotypic drug resistance testing, though numerous quality control issues were observed when this platform was used in a clinical laboratory setting. With additional refinement, the Sentosa SQ HIV-1 Genotyping Assay may contribute to the global efforts to control HIV-1.

A systematic review of the genetic mechanisms of dolutegravir resistance.

BACKGROUND: Characterizing the mutations selected by the integrase strand transfer inhibitor (INSTI) dolutegravir and their effects on susceptibility is essential for identifying viruses less likely to respond to dolutegravir therapy and for monitoring persons with virological failure (VF) on dolutegravir therapy. METHODS: We systematically reviewed dolutegravir resistance studies to identify mutations emerging under dolutegravir selection pressure, the effect of INSTI resistance mutations on in vitro dolutegravir susceptibility, and the virological efficacy of dolutegravir in antiretroviral-experienced persons. RESULTS AND CONCLUSIONS: We analysed 14 studies describing 84 in vitro passage experiments, 26 studies describing 63 persons developing VF plus INSTI resistance mutations on a dolutegravir-containing regimen, 41 studies describing dolutegravir susceptibility results, and 22 clinical trials and 16 cohort studies of dolutegravir-containing regimens. The most common INSTI resistance mutations in persons with VF on a dolutegravir-containing regimen were R263K, G118R, N155H and Q148H/R, with R263K and G118R predominating in previously INSTI-naive persons. R263K reduced dolutegravir susceptibility ∼2-fold. G118R generally reduced dolutegravir susceptibility >5-fold. The highest levels of reduced susceptibility occurred in viruses containing Q148 mutations in combination with G140 and/or E138 mutations. Dolutegravir two-drug regimens were highly effective for first-line therapy and for virologically suppressed persons provided dolutegravir's companion drug was fully active. Dolutegravir three-drug regimens were highly effective for salvage therapy in INSTI-naive persons provided one or more of dolutegravir's companion drugs was fully active. However, dolutegravir monotherapy in virologically suppressed persons and functional dolutegravir monotherapy in persons with active viral replication were associated with a non-trivial risk of VF plus INSTI resistance mutations.

Amino Acid Prevalence of HIV-1 pol Mutations by Direct Polymerase Chain Reaction and Single Genome Sequencing.

The presence of many highly unusual HIV-1 mutations at a minority variant threshold by next-generation sequence (NGS) may indicate that a high proportion of variants at or just above the threshold represent PCR errors. The validity of this hypothesis depends on the concept that highly unusual mutations detected by population-based sequencing are also highly unusual within a person's virus population. Highly unusual mutations were defined as mutations with a prevalence <0.01% in group M HIV-1 direct PCR population-based Sanger sequences in the Stanford HIV Drug Resistance Database. Single genome Sanger sequences [single genome sequences (SGSs)] were analyzed because they are not subject to PCR error. Permutation analyses compared the proportion of highly unusual mutations in SGSs with the empirical frequencies of these mutations in repeated random selections of population-based sequences. We created a database of 11,258 pol SGSs in 963 plasma samples from 345 persons with active virus replication and analyzed the subset of samples containing 10 or more SGSs. Highly unusual mutations occurred more commonly in samples undergoing SGS compared with population-based sequencing in protease (3.9% vs. 0.8%; p < .001), reverse transcriptase (6.5% vs. 1.5%; p < .001), and integrase (5.0% vs. 1.8%; p < .001). Highly unusual mutations occur more commonly in SGSs than in population-based sequences. However, they comprise a small proportion of all SGS mutations supporting the concept that the presence of many highly unusual mutations just above an NGS threshold suggests that the threshold is too low.

Multiplex Solid-Phase Melt Curve Analysis for the Point-of-Care Detection of HIV-1 Drug Resistance.

A point-of-care HIV-1 genotypic resistance assay that could be performed during a clinic visit would enable care providers to make informed treatment decisions for patients starting therapy or experiencing virologic failure on therapy. The main challenge for such an assay is the genetic variability at and surrounding each drug-resistance mutation (DRM). We analyzed a database of diverse global HIV sequences and used thermodynamic simulations to design an array of surface-bound oligonucleotide probe sets with each set sharing distinct 5' and 3' flanking sequences but having different centrally located nucleotides complementary to six codons at HIV-1 DRM reverse transcriptase position 103: AAA, AAC, AAG, AAT, AGA, and AGC. We then performed in vitro experiments using 80-mer oligonucleotides and PCR-amplified DNA from clinical plasma HIV-1 samples and culture supernatants that contained subtype A, B, C, D, CRF01_AE, and CRF02_AG viruses. Multiplexed solid-phase melt curve analysis discriminated perfectly among each of the six reported reverse transcriptase position 103 codons in both 80-mers and clinical samples. The sensitivity and specificity for detecting targets that contained AAC mixed with targets that contained AAA were >98% when AAC was present at a proportion of ≥10%. Multiplexed solid-phase melt curve analysis is a promising approach for developing point-of-care assays to distinguish between different codons in genetically variable regions such as those surrounding HIV-1 DRMs.

Selection analyses of paired HIV-1 gag and gp41 sequences obtained before and after antiretroviral therapy.

Most HIV-1-infected individuals with virological failure on a pharmacologically-boosted protease inhibitor (PI) regimen do not develop PI-resistance protease mutations. One proposed explanation is that HIV-1 gag or gp41 cytoplasmic domain mutations might also reduce PI susceptibility. In a recent study of paired gag and gp41 sequences from individuals with virological failure on a PI regimen, we did not identify PI-selected mutations and concluded that if such mutations existed, larger numbers of paired sequences from multiple studies would be needed for their identification. In this study, we generated site-specific amino acid profiles using gag and gp41 published sequences from 5,338 and 4,242 ART-naïve individuals, respectively, to assist researchers identify unusual mutations arising during therapy and to provide scripts for performing established and novel maximal likelihood estimates of dN/dS substitution rates in paired sequences. The pipelines used to generate the curated sequences, amino acid profiles, and dN/dS analyses will facilitate the application of consistent methods to paired gag and gp41 sequence datasets and expedite the identification of potential sites under PI-selection pressure.

Comparison of an In Vitro Diagnostic Next-Generation Sequencing Assay with Sanger Sequencing for HIV-1 Genotypic Resistance Testing.

The ability of next-generation sequencing (NGS) technologies to detect low frequency HIV-1 drug resistance mutations (DRMs) not detected by dideoxynucleotide Sanger sequencing has potential advantages for improved patient outcomes. We compared the performance of an in vitro diagnostic (IVD) NGS assay, the Sentosa SQ HIV genotyping assay for HIV-1 genotypic resistance testing, with Sanger sequencing on 138 protease/reverse transcriptase (RT) and 39 integrase sequences. The NGS assay used a 5% threshold for reporting low-frequency variants. The level of complete plus partial nucleotide sequence concordance between Sanger sequencing and NGS was 99.9%. Among the 138 protease/RT sequences, a mean of 6.4 DRMs was identified by both Sanger and NGS, a mean of 0.5 DRM was detected by NGS alone, and a mean of 0.1 DRM was detected by Sanger sequencing alone. Among the 39 integrase sequences, a mean of 1.6 DRMs was detected by both Sanger sequencing and NGS and a mean of 0.15 DRM was detected by NGS alone. Compared with Sanger sequencing, NGS estimated higher levels of resistance to one or more antiretroviral drugs for 18.2% of protease/RT sequences and 5.1% of integrase sequences. There was little evidence for technical artifacts in the NGS sequences, but the G-to-A hypermutation was detected in three samples. In conclusion, the IVD NGS assay evaluated in this study was highly concordant with Sanger sequencing. At the 5% threshold for reporting minority variants, NGS appeared to attain a modestly increased sensitivity for detecting low-frequency DRMs without compromising sequence accuracy.