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

Defective HIV-1 DNA pol sequences are not associated with HIV-1 DNA levels and drive most APOBEC-context drug resistance mutations

INTRODUCTION

The specificity of HIV-1 DNA genotypic resistance tests (GRTs) is hampered by the detection of the APOBEC-context drug resistance mutations (AC DRMs), usually harboured by replication-incompetent proviruses. We sought factors associated with defective sequences in the HIV-1 pol region. In addition, AC DRMs and their link with defective sequences were investigated.

METHODS

We included ART-treated patients with viral suppression or plasma viral load (VL) lower than 200 copies/mL, who underwent HIV-1 DNA genotyping, with successful sequencing of protease (PR), reverse transcriptase (RT) and integrase (IN) regions. Sequencing was performed using either the Sanger method or the Sentosa® NGS approach with a 20% cut-off. All hypermutated sequences and/or those containing at least one stop codon were considered defective.

RESULTS

A total of 613 HIV-1 DNA GRTs were analysed. Defective sequences were identified for 186 samples (30.3%) including 65 PR sequences, 92 RT sequences and 65 IN sequences. No association, including HIV-1 DNA levels, was found with the detection of defective pol sequences. A total of 226 AC DRMs were recorded in all sequences. Most of these mutations (78%) were harboured by defective sequences. AC DRMs did not emerge in the plasma viral population, and likely do not impact the virological response to ART.

CONCLUSIONS

Defectives pol sequences were not associated with HIV-1 DNA levels and harboured most of the AC DRMs. Such mutations likely have no clinical impact, and should not be reported in routine practice. Consensus guidelines for reporting HIV-1 DNA GRTs are needed, especially for the assessment and management of AC DRMs.

Question and answers on long-acting therapy with cabotegravir and rilpivirine in people with HIV

Modern ART is evolving, allowing the use of new drug formulations and alternative routes of administration to oral therapy. Long-acting (LA) cabotegravir and rilpivirine, the first fully injectable antiretroviral regimen approved for clinical use, is a test case for this new route of administration, and an innovation with implications for the quality of life of people with HIV (PWH). However, its use requires a reorganization of outpatient clinics and outpatient services, and a number of issues remain to be defined regarding the management of PWH on LA drugs, including the correct selection of people who can be treated with LA cabotegravir and rilpivirine. There is also ongoing debate about the best way to monitor both efficacy and tolerability of LA treatment and whether the management of virological failures and blips should be different from that reserved for oral regimens. The present article reviews the data on the use and management of LA cabotegravir and rilpivirine in different settings, with a review of clinical trial data and also the first available real-life experiences. The article focuses on the following: the reasons for the use of LA drugs; the implementation of their use in clinical practice; and the monitoring of treated people over time.

Proviral HIV DNA sequencing to inform suitability for long-acting injectable cabotegravir/rilpivirine

Increasingly, people attending HIV services are requesting long-acting injectable (LAI) antiretroviral treatment (ART). However, without HIV RNA resistance-associated mutation results, individuals are considered unsuitable for LAI ART. We present our experience of sequencing proviral HIV DNA in 30 individuals to inform suitability for LAI ART, of whom 23 were considered suitable. In conclusion, optimization of diagnostic tools such as proviral HIV DNA sequencing to confirm suitability for LAI ART would be a welcome addition.

Positive selection, genetic recombination, and intra-host evolution in novel equine coronavirus genomes and other members of the Embecovirus subgenus

There are several examples of coronaviruses in the Betacoronavirus subgenus Embecovirus that have jumped from an animal to the human host. Studying how evolutionary factors shape coronaviruses in non-human hosts may provide insight into the coronavirus host-switching potential. Equids, such as horses and donkeys, are susceptible to equine coronaviruses (ECoVs). With increased testing prevalence, several ECoV genome sequences have become available for molecular evolutionary analyses, especially those from the United States of America (USA). To date, no analyses have been performed to characterize evolution within coding regions of the ECoV genome. Here, we obtain and describe four new ECoV genome sequences from infected equines from across the USA presenting clinical symptoms of ECoV, and infer ECoV-specific and Embecovirus-wide patterns of molecular evolution. Within two of the four data sets analyzed, we find evidence of intra-host evolution within the nucleocapsid (N) gene, suggestive of quasispecies development. We also identify 12 putative genetic recombination events within the ECoV genome, 11 of which fall in ORF1ab. Finally, we infer and compare sites subject to positive selection on the ancestral branch of each major Embecovirus member clade. Specifically, for the two currently identified human coronavirus (HCoV) embecoviruses that have spilled from animals to humans (HCoV-OC43 and HCoV-HKU1), we find that there are 42 and 2 such sites, respectively, perhaps reflective of the more complex ancestral evolutionary history of HCoV-OC43, which involves several different animal hosts.IMPORTANCEThe Betacoronavirus subgenus Embecovirus contains coronaviruses that not only pose a health threat to animals and humans, but also have jumped from animal to human host. Equids, such as horses and donkeys are susceptible to equine coronavirus (ECoV) infections. No studies have systematically examined evolutionary patterns within ECoV genomes. Our study addresses this gap and provides insight into intra-host ECoV evolution from infected horses. Further, we identify and report natural selection pattern differences between two embecoviruses that have jumped from animals to humans [human coronavirus OC43 and HKU1 (HCoV-OC43 and HCoV-HKU1, respectively)], and hypothesize that the differences observed may be due to the different animal host(s) that each virus circulated in prior to its jump into humans. Finally, we contribute four novel, high-quality ECoV genomes to the scientific community.

A Sensitivity and Consistency Comparison Between Next-Generation Sequencing and Sanger Sequencing in HIV-1 Pretreatment Drug Resistance Testing

Next-generation sequencing (NGS) for HIV drug resistance (DR) testing has an increasing number of applications for the detection of low-abundance drug-resistant variants (LA-DRVs) in regard to its features as a quasi-species. However, there is less information on its detection performance in DR detection with NGS. To determine the feasibility of using NGS technology in LA-DRV detection for HIV-1 pretreatment drug resistance, 80 HIV-infected individuals who had never undergone antiretroviral therapy were subjected to both NGS and Sanger sequencing (SS) in HIV-1 drug resistance testing. The results reported in this study show that NGS exhibits higher sensitivity for drug resistance identification than SS at a 5% detection threshold. NGS showed a better consistency compared with that of SS for both protease inhibitors (PIs) and integrase inhibitors (INSTIs), with a figure amounting to more than 90%, but worse consistency in nucleotide reverse transcriptase inhibitors (NRTIs), with a consistency ranging from only 61.25% to 87.50%. The consistency of non-nucleotide reverse transcriptase inhibitors (NNRTIs) between NGS and SS was around 85%. NGS showed the highest sensitivity of 87.0% at a 5% threshold. The application of NGS technology in HIV-1 genotype resistance detection in different populations infected with HIV requires further documentation and validation.

Archived HIV-1 Drug Resistance Mutations: Role of Proviral HIV-1 DNA Genotype for the Management of Virological Responder People Living with HIV

Despite its effectiveness in controlling plasma viremia, antiretroviral therapy (ART) cannot target proviral DNA, which remains an obstacle to HIV-1 eradication. When treatment is interrupted, the reservoirs can act as a source of viral rebound, highlighting the value of proviral DNA as an additional source of information on an individual's overall resistance burden. In cases where the viral load is too low for successful HIV-1 RNA genotyping, HIV-1 DNA can help identify resistance mutations in treated individuals. The absence of treatment history, the need to adjust ART despite undetectable viremia, or the presence of LLV further support the use of genotypic resistance tests (GRTs) on HIV-1 DNA. Conventionally, GRTs have been achieved through Sanger sequencing, but the advances in NGS are leading to an increase in its use, allowing the detection of minority variants present in less than 20% of the viral population. The clinical significance of these mutations remains under debate, with interpretations varying based on context. Additionally, proviral DNA is subject to APOBEC3-induced hypermutation, which can lead to defective, nonviable viral genomes, a factor that must be considered when performing GRTs on HIV-1 DNA.

Genetic Diversity From Proviral DNA as a Proxy for Time Since HIV-1 Infection

HIV-1 RNA genetic diversity predicts time since infection, which is important for clinical care and research. It is unclear, however, whether proviral DNA genetic diversity sampled under suppressive antiretroviral therapy can be used for this purpose. We tested whether proviral genetic diversity from next-generation sequencing predicts time since infection and recency in 221 people with HIV-1 with known infection time. Proviral diversity was significantly associated with time since infection (P < 5×10-7, R2 up to 25%) and predictive of treatment initiation during recent infection (area under the curve-receiver operating characteristic up to 0.85). This shows the utility of proviral genetic diversity as a proxy for time since infection.

Comparison of Different HIV-1 Resistance Interpretation Tools for Next-Generation Sequencing in Italy

BACKGROUND

Next-generation sequencing (NGS) is gradually replacing Sanger sequencing for HIV genotypic drug resistance testing (GRT). This work evaluated the concordance among different NGS-GRT interpretation tools in a real-life setting.

METHODS

Routine NGS-GRT data were generated from viral RNA at 11 Italian laboratories with the AD4SEQ HIV-1 Solution v2 commercial kit. NGS results were interpreted by the SmartVir system provided by the kit and by two online tools (HyDRA Web and Stanford HIVdb). NGS-GRT was considered valid when the coverage was >100 reads (100×) at each PR/RT/IN resistance-associated position listed in the HIVdb 9.5.1 algorithm.

RESULTS

Among 629 NGS-GRT, 75.2%, 74.2%, and 70.9% were valid according to SmartVir, HyDRA Web, and HIVdb. Considering at least two interpretation tools, 463 (73.6%) NGS-GRT had a valid coverage for resistance analyses. The proportion of valid samples was affected by viremia <10,000-1000 copies/mL and non-B subtypes. Mutations at an NGS frequency >10% showed fair concordance among different interpretation tools.

CONCLUSION

This Italian survey on NGS resistance testing suggests that viremia levels and HIV subtype affect NGS-GRT coverage. Within the current routine method for NGS-GRT, only mutations with frequency >10% seem reliably detected across different interpretation tools.

Use of next-generation sequencing on HIV-1 DNA to assess archived resistance in highly treatment-experienced people with multidrug-resistant HIV under virological control: data from the PRESTIGIO Registry

BACKGROUND

To clarify whether next-generation sequencing (NGS) can be useful for resistance assessment in virologically suppressed highly treatment-experienced (HTE) individuals with MDR HIV.

METHODS

Ninety-one participants from the PRESTIGIO Registry were included. NGS was performed on HIV-DNA at 1%, 5% and 20% cut-offs; major drug resistance mutations (DRMs) were evaluated and compared with those detected in historical plasma genotypic resistance testing (h-GRT). APOBEC editing was also characterized.

RESULTS

Participants had a complex and long treatment history [median 23 (IQR 21-25) years of ART exposure) and had been virologically suppressed since a median of 3 (IQR 2-5) years. Among all major DRMs detected by HIV-DNA NGS and/or h-GRT, 30% were exclusively found through NGS. The highest detection rate of historical major DRMs was reached with NGS set at 1%, but unusual substitutions and extensive APOBEC hypermutations suggest technical issues and poor clinical relevance in the 1%-5% interval. At NGS set at 5%, 67.2% of historical major DRMs were detected. The number of major DRMs detected exclusively by DNA-NGS as minority variants (frequency 5%-20%) was significantly higher in individuals who later experienced virological rebound compared with those who maintained virological control [median 2 (IQR 1-3) versus 1 (0-2), P = 0.030] and positively correlated with viraemia levels at rebound (rho = 0.474, P = 0.030).

CONCLUSIONS

In non-viraemic people with an MDR virus, HIV-1 DNA NGS set at 5% is an acceptable technical cut-off that might help to reveal mutations with a potential clinical relevance. Moreover, the number of minority resistance mutations additionally detected by NGS might be associated with loss of virological control.

Performance of the Applied Biosystems HIV-1 Genotyping Kit with Integrase

HIV genotyping is used to assess HIV susceptibility to antiretroviral drugs. The Applied Biosystems HIV-1 Genotyping Kit with Integrase (AB kit, Thermo Fisher Scientific) detects resistance-associated mutations (RAMs) in HIV protease (PR), reverse transcriptase (RT), and integrase (IN). We compared results from the AB kit with results obtained previously with the ViroSeq HIV-1 Genotyping System. DNA amplicons from the AB kit were also analyzed using next-generation sequencing (NGS). HIV RNA was extracted using the MagNA Pure 24 instrument (Roche Diagnostics; 96 plasma samples, HIV subtype B, viral load range: 530-737,741 copies/mL). FASTA files were generated from AB kit data using Exatype (Hyrax Biosciences). DNA amplicons from the AB kit were also analyzed by NGS using the Nextera XT kit (Illumina). Drug resistance was predicted using the Stanford HIV Drug Resistance Database. The mean genetic distance for sequences from ViroSeq and the AB kit was 0.02% for PR/RT and 0.04% for IN; 103 major RAMs were detected by both methods. Four additional major RAMs were detected by the AB kit only. These four major RAMs were also detected by NGS (detected in 18.1%-38.2% of NGS reads). NGS detected 27 major RAMs that were not detected with either of the Sanger sequencing-based kits. All major RAMs detected with ViroSeq were detected with the AB kit; additional RAMs were detected with the AB kit only. DNA amplicons from the AB kit can be used for NGS for more sensitive detection of RAMs.

Population-based nanopore sequencing of the HIV-1 pangenome to identify drug resistance mutations

HIV-1 drug resistance genotypic tests have primarily been performed by Sanger sequencing of gene segments encoding different drug target proteins. Since the number of targets has increased with the addition of a new class of antiretroviral drugs, a simple high-throughput system for assessing nucleotide sequences throughout the HIV-1 genome is required. Here, we developed a new solution using nanopore sequencing of viral pangenomes amplified by PCR. Benchmark tests using HIV-1 molecular clones demonstrated an accuracy of up to 99.9%. In addition, validation tests of our protocol in 106 clinical samples demonstrated high concordance of drug resistance and tropism genotypes (92.5% and 98.1%, respectively) between the nanopore sequencing-based results and archived clinical determinations made based on Sanger sequencing data. These results suggest that our new approach will be a powerful solution for the comprehensive survey of HIV-1 drug resistance mutations in clinical settings.

Low-level viremia episodes appear to affect the provirus composition of the circulating cellular HIV reservoir during antiretroviral therapy

Low-level viremia (LLV) ranging from 50 to 1,000 copies/ml is common in most HIV-1-infected patients receiving antiretroviral therapy (ART). However, the source of LLV and the impact of LLV on the HIV-1 reservoir during ART remain uncertain. We hypothesized that LLV may arise from the HIV reservoir and its occurrence affect the composition of the reservoir after LLV episodes. Accordingly, we investigated the genetic linkage of sequences obtained from plasma at LLV and pre-ART time points and from peripheral blood mononuclear cells (PBMCs) at pre-ART, pre-LLV, LLV, and post-LLV time points. We found that LLV sequences were populated with a predominant viral quasispecies that accounted for 67.29%∼100% of all sequences. Two episodes of LLV in subject 1, spaced 6 months apart, appeared to have originated from the stochastic reactivation of latently HIV-1-infected cells. Moreover, 3.77% of pre-ART plasma sequences were identical to 67.29% of LLV-3 plasma sequences in subject 1, suggesting that LLV may have arisen from a subset of cells that were infected before ART was initiated. No direct evidence of sequence linkage was found between LLV viruses and circulating cellular reservoirs in all subjects. The reservoir size, diversity, and divergence of the PBMC DNA did not differ significantly between the pre- and post-LLV sampling points (P > 0.05), but the composition of viral reservoir quasispecies shifted markedly before and after LLV episodes. Indeed, subjects with LLV had a higher total PBMC DNA level, greater viral diversity, a lower proportion of variants with identical sequences detected at two or more time points, and a shorter variant duration during ART compared with subjects without LLV. Overall, our findings suggested that LLV viruses may stem from an unidentified source other than circulating cellular reservoirs. LLV episodes may introduce great complexity into the HIV reservoir, which brings challenges to the development of treatment strategies.

Use of genotypic HIV DNA testing: a DELPHI-type consensus

OBJECTIVES

As many disparities in the clinical use of HIV DNA sequencing are observed, a DELPHI-type consensus was initiated in France to homogenize use, techniques and interpretation of results.

METHODS

Based on a literature review and clinical experience, a steering committee (SC) of eight virologists and one infectious disease specialist formulated statements. Statements were submitted to an independent and anonymous electronic vote of virologists and HIV clinicians in France, between October 2022 and December 2022.

RESULTS

The SC developed 20 statements grouped into six categories: clinical situations for the use of HIV DNA genotyping; techniques for performing HIV DNA genotyping; consideration of apolipoprotein B mRNA editing enzyme (APOBEC) mutations; genotyping results reporting; recycling of antiretrovirals; and availability of HIV DNA genotyping tests and delays. Twenty-one virologists and 47 clinicians participated in two voting rounds and 18/20 (90%) assertions reached a 'strong' consensus. For example, that prior genotyping on HIV DNA is useful for clinical decision-making when considering switching to some long-acting regimens or to reduce the number of antiretroviral agents in virologically suppressed patients for whom RNA data are unavailable/not exploitable/not sufficiently informative. Two statements achieved no consensus: reporting any detected viral minority population for discussion in multidisciplinary meetings (virologists), and possible risk of virological failure when using a second-generation InSTI plus lamivudine or emtricitabine regimen in patients with undetectable viral load within ≥1 year and in the presence of a documented M184V mutation within the last 5 years (clinicians).

CONCLUSIONS

This DELPHI-type consensus will facilitate the strengthening and harmonization of good practice when performing HIV DNA sequencing.

HIV-1 Drug Resistance Detected by Next-Generation Sequencing among ART-Naïve Individuals: A Systematic Review and Meta-Analysis

BACKGROUND

There are an increasing number of articles focused on the prevalence and clinical impact of pretreatment HIV drug resistance (PDR) detected by Sanger sequencing (SGS). PDR may contribute to the increased likelihood of virologic failure and the emergence of new resistance mutations. As SGS is gradually replaced by next-generation sequencing (NGS), it is necessary to assess the levels of PDR using NGS in ART-naïve patients systematically. NGS can detect the viral variants (low-abundance drug-resistant HIV-1 variants (LA-DRVs)) of virus quasi-species at levels below 20% that SGS may fail to detect. NGS has the potential to optimize current HIV drug resistance surveillance methods and inform future research directions. As the NGS technique has high sensitivity, it is highly likely that the level of pretreatment resistance would be underestimated using conventional techniques.

METHODS

For the systematic review and meta-analysis, we searched for original studies published in PubMed, Web of Science, Scopus, and Embase before 30 March 2023 that focused exclusively on the application of NGS in the detection of HIV drug resistance. Pooled prevalence estimates were calculated using a random effects model using the 'meta' package in R (version 4.2.3). We described drug resistance detected at five thresholds (>1%, 2%, 5%, 10%, and 20% of virus quasi-species). Chi-squared tests were used to analyze differences between the overall prevalence of PDR reported by SGS and NGS.

RESULTS

A total of 39 eligible studies were selected. The studies included a total of 15,242 ART-naïve individuals living with HIV. The prevalence of PDR was inversely correlated with the mutation detection threshold. The overall prevalence of PDR was 29.74% at the 1% threshold, 22.43% at the 2% threshold, 15.47% at the 5% threshold, 12.95% at the 10% threshold, and 11.08% at the 20% threshold. The prevalence of PDR to INSTIs was 1.22% (95%CI: 0.58-2.57), which is the lowest among the values for all antiretroviral drugs. The prevalence of LA-DRVs was 9.45%. At the 2% and 20% detection threshold, the prevalence of PDR was 22.43% and 11.08%, respectively. Resistance to PIs and INSTIs increased 5.52-fold and 7.08-fold, respectively, in those with a PDR threshold of 2% compared with those with PDR at 20%. However, resistance to NRTIs and NNRTIs increased 2.50-fold and 2.37-fold, respectively. There was a significant difference between the 2% and 5% threshold for detecting HIV drug resistance. There was no statistically significant difference between the results reported by SGS and NGS when using the 20% threshold for reporting resistance mutations.

CONCLUSION

In this study, we found that next-generation sequencing facilitates a more sensitive detection of HIV-1 drug resistance than SGS. The high prevalence of PDR emphasizes the importance of baseline resistance and assessing the threshold for optimal clinical detection using NGS.

Inflammatory profile of vertically HIV-1 infected adolescents receiving ART in Cameroon: a contribution toward optimal pediatric HIV control strategies

Antiretroviral therapy (ART) has improved the lifespan of people living with HIV. However, their immune system remains in a state of sustained activation/inflammation, which favors viral replication and depletion of helper T-cells with varying profiles according to ART-response. We herein sought to ascertain the inflammatory profile of adolescents living with perinatal HIV-1 infection (ALPHI) receiving ART in an African context. In this cross-sectional and comparative study among ART-experienced ALPHI in Yaoundé-Cameroon, HIV-1 RNA was measured by Abbott Real-time PCR; CD4 cells were enumerated using flow cytometry; serum cytokines were measured by ELISA; HIV-1 proviral DNA was genotyped by Sanger-sequencing; and archived drug resistance mutations (ADRMs) were interpreted using Stanford HIVdb.v9.0.1. Overall, 73 adolescents were enrolled (60 ALPHI and 13 HIV-1 negative peers) aged 15 (13-18) years; 60.00% were female. ART median duration was 92 (46-123) months; median viral load was 3.99 (3.17-4.66) RNA Log10 (copies)/mL and median CD4 count was 326 (201-654) cells/mm3. As compared to HIV-negative adolescents, TNFα was highly expressed among ALPHI (p<0.01). Following a virological response, inflammatory cytokines (IFNγ and IL-12), anti-inflammatory cytokines (IL-4 and IL-10) and inflammation-related cytokines (IL-6 and IL-1β) were highly expressed with viral suppression (VS) vs. virological failure (VF), while the chemokine CCL3 was highly expressed with VF (p<0.01). Regarding the immune response, the inflammatory cytokine TNFα was highly expressed in those that are immunocompetent (CD4≥500 cell/mm3) vs. immunocompromised (CD4<500 cell/mm3), p ≤ 0.01; while chemokine CCL2 was highly expressed in the immunocompromised (p<0.05). In the presence of ADRMs, IL-4 and CCL3 were highly expressed (p=0.027 and p=0.043 respectively). Among ART-experienced ALPHI in Cameroon, the TNFα cytokine was found to be an inflammatory marker of HIV infection; IFNγ, IL-1β, IL-6, and IL-12 are potential immunological markers of VS and targeting these cytokines in addition to antiretroviral drugs may improve management. Moreover, CCL3 and CCL2 are possible predictors of VF and/or being immunocompromised and could serve as surrogates of poor ART response.

Brief Report: Comparative Analysis of Pre-existing HIV Drug Resistance Mutations in Proviral DNA Using Next-Generation Sequencing and Routine HIV RNA Genotyping

BACKGROUND

We investigated whether deep sequencing of archived HIV DNA of antiretroviral-naive persons with acute/early HIV infection could identify transmitted drug resistance mutations (DRM), per the IAS drug resistance algorithm, which are not detected by routine bulk (consensus) sequencing.

METHODS

Deep sequencing of HIV DNA from peripheral blood mononuclear cells and consensus sequencing from concurrent blood plasma (BP) was performed from antiretroviral (ART)-naive adults with recent infection. We compared the prevalence of low-frequency (2%-20%) and high-frequency (>20%) nonnucleoside reverse transcriptase inhibitor (NNRTI), nucleoside reverse transcriptase inhibitor (NRTI), and protease inhibitor (PI) DRM.

RESULTS

Overall, 190 individuals were included, 72 (37.9%) with acute, 20 (10.5%) with very early, and 98 (51.6%) with recent HIV infection. Although all DRM detected in plasma appeared in archived proviral DNA, 9 high-frequency mutations were only detected in HIV DNA. These included 3 NRTI mutations, 4 NNRTI mutations, 1 PI mutation, and 1 H221Y (associated rilpivirine resistance) mutation. When considering DRM <20%, 11 NNRTI, 7 NRTI, 6 PI, and 3 F227L (associated doravirine resistance) mutations were found exclusively in HIV DNA. Interestingly, although 2 high-frequency M184V appeared in both DNA and RNA, low-frequency M184I were exclusive to HIV DNA (n = 6). No participants experienced virologic failure after initiating ART during the median 25.39 ± 3.13 months of follow-up on treatment.

CONCLUSION

Although most high-frequency DRMs were consistently detected in HIV RNA and HIV DNA, the presence of low-frequency DRM in proviral DNA may be relevant for clinicians because these mutations could become dominant under drug selection pressure.

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.

Evaluation of Circulating and Archived HIV-1 Integrase Drug-Resistance Variants among Patients on Third-Line ART in Cameroon: Implications for Dolutegravir-Containing Regimens in Resource-Limited Settings

To ensure the long-term efficacy of dolutegravir (DTG), we evaluated the genotypic profile in viral reservoirs among patients on third-line (3L) antiretroviral therapy (ART) in Cameroon, according to prior exposure to raltegravir (RAL). A facility-based study was conducted from May through December 2021 among patients on 3L ART from HIV treatment centers in Yaoundé and Douala. Viral load was measured, and genotyping was performed on plasma RNA and proviral DNA. HIV-1 drug resistance mutations were interpreted using HIVdb.v9.1 and phylogeny analysis was performed using MEGA.v7, with P < 0.05 considered significant. Of the 12,093 patients on ART, 53 fully met our inclusion criteria. The median (IQR) age was 51 years (40 to 55 years), and the male/female ratio was 4/5. The median duration on integrase strand-transfer inhibitors (INSTI)-containing regimens was 18 months (12 to 32 months), and 15.09% (8/53) were exposed to RAL. The most administered 3L ART was TDF+3TC+DTG+DRV/r (33.96%, 18/53). Only 5.66% (3/53) had unsuppressed viremia (>1000 copies/mL). Resistance testing in proviral DNA was successful for 18/22 participants and revealed 1/18 patients (5.56%, in the RAL-arm) with archived mutations at major resistance positions (G140R and G163R). Five subtypes were identified, CRF02_AG (12/18), CRF22_01AE (3/18), A1 (1/18), G (1/18), and F2 (1/18). In Cameroon, 3L-experienced patients had a good virological response with a low level of archived mutations in the integrase. This finding underscored the use of DTG-containing ART for heavily treated patients in similar programmatic settings. However, patients with prior exposure to RAL should be closely monitored following a stratified or personalized approach to mitigate risks of INSTI-resistance, alongside pharmacovigilance. IMPORTANCE We described the analysis of the genotypes of the population within third-line antiviral therapy in Cameroon, with a focus on defining the effects of prior raltegravir (RAL) treatment and resistance mutations for current dolutegravir (DTG) treatment. While supporting the current transition to DTG-containing ART in resource-limited settings toward the achievement of the UNAIDS' goal of HIV elimination by 2030, our findings suggested that RAL-exposed patients may need a specific monitoring approach either in a stratified or personalized model of third-line ART to ensure the long-term success of DTG-containing regimens.

Genotypic Resistance Testing of HIV-1 DNA in Peripheral Blood Mononuclear Cells

HIV-1 DNA exists in nonintegrated linear and circular episomal forms and as integrated proviruses. In patients with plasma viremia, most peripheral blood mononuclear cell (PBMC) HIV-1 DNA consists of recently produced nonintegrated virus DNA while in patients with prolonged virological suppression (VS) on antiretroviral therapy (ART), most PBMC HIV-1 DNA consists of proviral DNA produced months to years earlier. Drug-resistance mutations (DRMs) in PBMCs are more likely to coexist with ancestral wild-type virus populations than they are in plasma, explaining why next-generation sequencing is particularly useful for the detection of PBMC-associated DRMs. In patients with ongoing high levels of active virus replication, the DRMs detected in PBMCs and in plasma are usually highly concordant. However, in patients with lower levels of virus replication, it may take several months for plasma virus DRMs to reach detectable levels in PBMCs. This time lag explains why, in patients with VS, PBMC genotypic resistance testing (GRT) is less sensitive than historical plasma virus GRT, if previous episodes of virological failure and emergent DRMs were either not prolonged or not associated with high levels of plasma viremia. Despite the increasing use of PBMC GRT in patients with VS, few studies have examined the predictive value of DRMs on the response to a simplified ART regimen. In this review, we summarize what is known about PBMC HIV-1 DNA dynamics, particularly in patients with suppressed plasma viremia, the methods used for PBMC HIV-1 GRT, and the scenarios in which PBMC GRT has been used clinically.

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.

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.

Viral resistance burden and APOBEC editing correlate with virological response in heavily treatment-experienced people living with multi-drug resistant HIV

BACKGROUND

The impact of drug resistance mutational load and APOBEC editing in heavily treatment-experienced (HTE) people living with multidrug-resistant HIV has not been investigated.

MATERIAL AND METHODS

This study explored the HIV-DNA and HIV-RNA mutational load of drug resistance and APOBEC-related mutations through next-generation sequencing (NGS, Illumina MiSeq) in 20 failing HTE participants enrolled in the PRESTIGIO registry.

RESULTS

The patients showed high levels of both HIV-DNA (4.5 [4.0-5.2] log₁₀ copies/10⁶ T-CD4+ cell) and HIV-RNA (4.5 [4.1-5.0] log₁₀ copies/mL) with complex resistance patterns in both compartments. Among the 255 drug-resistant mutations found, 66.3% were concordantly detected in both HIV-DNA and HIV-RNA; 71.3% of mutations were already present in historical Sanger genotypes. At an intra-patient frequency > 5%, a considerable proportion of mutations detected through DNA-NGS were found in historical genotypes but not through RNA-NGS, and few patients had APOBEC-related mutations. Of 14 patients who switched therapy, the five who failed treatment had DNA resistance with higher intra-patient frequency and higher DNA/RNA mutational load in a context of tendentially less pronounced APOBEC editing compared with those who responded.

CONCLUSIONS

Using NGS in HIV-DNA and HIV-RNA together with APOBEC editing evaluation might help to identify HTE individuals with MDR who are more prone to experience virological failure.

Patterns of pretreatment drug resistance mutations of very early diagnosed and treated infants in Botswana

OBJECTIVE

To describe the occurrence of HIV drug resistance mutations (DRMs) in both intact and defective HIV-1 cell-associated DNA (HIV-1 CAD) among early-treated infants.

DESIGN

The Botswana EIT Study (ClinicalTrials.gov NCT02369406) initiated antiretroviral therapy (ART) in the first week of life and evaluated HIV-1 in plasma and peripheral blood mononuclear cells (PBMCs).

METHODOLOGY

We analyzed 257 near-HIV-1 full-length sequences (nFLS) obtained by Illumina next-generation sequencing from infants near birth. Sanger sequencing of pol was performed for mothers at delivery and children with clinical failure through 96 weeks. DRMs were identified using the Stanford HIV Drug Resistance Database.

RESULTS

In 27 infants, median PBMC HIV-1 proviral load was 492 copies/ml [IQR: 78-1246 copies/ml] at a median of 2 days (range 1-32); 18 (66.7%) had no DRMs detected; six (22.2%) had DRMs detected in defective DNA only, and three (11.1%) had DRMs in both defective and intact DNA (P = 0.09). A total of 60 of 151 (37.7%) defective sequences had at least one DRM: 31.8% NNRTI, 15.2% NRTI, 5.3% protease inhibitor, and 15.5% INSTI-associated mutations. In intact sequences, 33 of 106 (31.1%) had at least 1 DRM: 29.2% NNRTI, 7.5% NRTI, 0.9% protease inhibitor, and no INSTI-associated mutations. For all three infants with intact sequence DRMs, corresponding DRMs occurred in maternal plasma at delivery. Archived DRMs were detectable at a later clinical rebound on only one occasion.

CONCLUSION

Defective HIV-1 cell-associated DNA sequences may overestimate the prevalence of drug resistance among early-treated children. The impact of DRMs from intact proviruses on long-term treatment outcomes warrants further investigation.

Diversity of HIV-1 subtypes and transmitted drug-resistance mutations among minority HIV-1 variants in a Turkish cohort

BACKGROUND

The World Health Organization (WHO) recommends the surveillance of transmitted drug resistance mutations (TDRMs) to ensure the effectiveness and sustainability of HIV treatment programs.

OBJECTIVE

Our aim was to determine the TDRMs and evaluate the distribution of HIV-1 subtypes using and compared next-generation sequencing (NGS) and Sanger-based sequencing (SBS) in a cohort of 44 antiretroviral treatment-naïve patients.

METHODS

All samples that were referred to the microbiology laboratory for HIV drug resistance analysis between December 2016 and February 2018 were included in the study. After exclusions, 44 treatment-naive adult patients with a viral load of >1000 copies/mL were analyzed. DNA sequencing for reverse transcriptase and protease regions was performed using both DeepChek ABL single round kit and Sanger-based ViroSeq HIV-1 Genotyping System. The mutations and HIV-1 subtypes were analyzed using the Stanford HIVdb version 8.6.1 Genotypic Resistance software, and TDRMs were assessed using the WHO surveillance drug-resistance mutation database. HIV-1 subtypes were confirmed by constructing a maximum-likelihood phylogenetic tree using Los Alamos IQ-Tree software.

RESULTS

NGS identified nucleos(t)ide reverse transcriptase inhibitor (NRTI)-TDRMs in 9.1% of the patients, non-nucleos(t)ide reverse transcriptase inhibitor (NNRTI)-TDRMs in 6.8% of the patients, and protease inhibitor (PI)-TDRMs in 18.2% of the patients at a detection threshold of ≥1%. Using SBS, 2.3% and 6.8% of the patients were found to have NRTI- and NNRTI-TDRMs, respectively, but no major PI mutations were detected. M41L, L74I, K65R, M184V, and M184I related to NRTI, K103N to NNRTI, and N83D, M46I, I84V, V82A, L24I, L90M, I54V to the PI sites were identified using NGS. Most mutations were found in low-abundance (frequency range: 1.0% - 4.7%) HIV-1 variants, except M41L and K103N. The subtypes of the isolates were found as follows; 61.4% subtype B, 18.2% subtype B/CRF02_AG recombinant, 13.6% subtype A, 4.5% CRF43_02G, and 2.3% CRF02_AG. All TDRMs, except K65R, were detected in HIV-1 subtype B isolates.

CONCLUSION

The high diversity of protease site TDRMs in the minority HIV-1 variants and prevalence of CRFs were remarkable in this study. All minority HIV-1 variants were missed by conventional sequencing. TDRM prevalence among minority variants appears to be decreasing over time at our center.

Human APOBEC3 Variations and Viral Infection

Human APOBEC3 (apolipoprotein B mRNA-editing catalytic polypeptide-like 3) enzymes are capable of inhibiting a wide range of endogenous and exogenous viruses using deaminase and deaminase-independent mechanisms. These enzymes are essential components of our innate immune system, as evidenced by (a) their strong positive selection and expansion in primates, (b) the evolution of viral counter-defense mechanisms, such as proteasomal degradation mediated by HIV Vif, and (c) hypermutation and inactivation of a large number of integrated HIV-1 proviruses. Numerous APOBEC3 single nucleotide polymorphisms, haplotypes, and splice variants have been identified in humans. Several of these variants have been reported to be associated with differential antiviral immunity. This review focuses on the current knowledge in the field about these natural variations and their roles in infectious diseases.

Multi-Laboratory Comparison of Next-Generation to Sanger-Based Sequencing for HIV-1 Drug Resistance Genotyping

Next-generation sequencing (NGS) is increasingly used for HIV-1 drug resistance genotyping. NGS methods have the potential for a more sensitive detection of low-abundance variants (LAV) compared to standard Sanger sequencing (SS) methods. A standardized threshold for reporting LAV that generates data comparable to those derived from SS is needed to allow for the comparability of data from laboratories using NGS and SS. Ten HIV-1 specimens were tested in ten laboratories using Illumina MiSeq-based methods. The consensus sequences for each specimen using LAV thresholds of 5%, 10%, 15%, and 20% were compared to each other and to the consensus of the SS sequences (protease 4-99; reverse transcriptase 38-247). The concordance among laboratories' sequences at different thresholds was evaluated by pairwise sequence comparisons. NGS sequences generated using the 20% threshold were the most similar to the SS consensus (average 99.6% identity, range 96.1-100%), compared to 15% (99.4%, 88.5-100%), 10% (99.2%, 87.4-100%), or 5% (98.5%, 86.4-100%). The average sequence identity between laboratories using thresholds of 20%, 15%, 10%, and 5% was 99.1%, 98.7%, 98.3%, and 97.3%, respectively. Using the 20% threshold, we observed an excellent agreement between NGS and SS, but significant differences at lower thresholds. Understanding how variation in NGS methods influences sequence quality is essential for NGS-based HIV-1 drug resistance genotyping.

Dry Panels Supporting External Quality Assessment Programs for Next Generation Sequencing-Based HIV Drug Resistance Testing

External quality assessment (EQA) is a keystone element in the validation and implementation of next generation sequencing (NGS)-based HIV drug resistance testing (DRT). Software validation and evaluation is a critical element in NGS EQA programs. While the development, sharing, and adoption of wet lab protocols is coupled with the increasing access to NGS technology worldwide, rendering it easy to produce NGS data for HIV-DRT, bioinformatic data analysis remains a bottleneck for most of the diagnostic laboratories. Several computational tools have been made available, via free or commercial sources, to automate the conversion of raw NGS data into an actionable clinical report. Although different software platforms yield equivalent results when identical raw NGS datasets are analyzed for variations at higher abundance, discrepancies arise when variations at lower frequencies are considered. This implies that validation and performance assessment of the bioinformatics tools applied in NGS HIV-DRT is critical, and the origins of the observed discrepancies should be determined. Well-characterized reference NGS datasets with ground truth on the genotype composition at all examined loci and the exact frequencies of HIV variations they may harbor, so-called dry panels, would be essential in such cases. The strategic design and construction of such panels are challenging but imperative tasks in support of EQA programs for NGS-based HIV-DRT and the validation of relevant bioinformatics tools. Here, we present criteria that can guide the design of such dry panels, which were discussed in the Second International Winnipeg Symposium themed for EQA strategies for NGS HIVDR assays.

Quality Control of Next-Generation Sequencing-Based HIV-1 Drug Resistance Data in Clinical Laboratory Information Systems Framework

Next-generation sequencing (NGS) in HIV drug resistance (HIVDR) testing has the potential to improve both clinical and public health settings, however it challenges the normal operations of quality management systems to be more flexible due to its complexity, massive data generation, and rapidly evolving protocols. While guidelines for quality management in NGS data have previously been outlined, little guidance has been implemented for NGS-based HIVDR testing. This document summarizes quality control procedures for NGS-based HIVDR testing laboratories using a laboratory information systems (LIS) framework. Here, we focus in particular on the quality control measures applied on the final sequencing product aligned with the recommendations from the World Health Organization HIV Drug Resistance Laboratory Network.

Next-Generation Sequencing for HIV Drug Resistance Testing: Laboratory, Clinical, and Implementation Considerations

Higher accessibility and decreasing costs of next generation sequencing (NGS), availability of commercial kits, and development of dedicated analysis pipelines, have allowed an increasing number of laboratories to adopt this technology for HIV drug resistance (HIVDR) genotyping. Conventional HIVDR genotyping is traditionally carried out using population-based Sanger sequencing, which has a limited capacity for reliable detection of variants present at intra-host frequencies below a threshold of approximately 20%. NGS has the potential to improve sensitivity and quantitatively identify low-abundance variants, improving efficiency and lowering costs. However, some challenges exist for the standardization and quality assurance of NGS-based HIVDR genotyping. In this paper, we highlight considerations of these challenges as related to laboratory, clinical, and implementation of NGS for HIV drug resistance testing. Several sources of variation and bias occur in each step of the general NGS workflow, i.e., starting material, sample type, PCR amplification, library preparation method, instrument and sequencing chemistry-inherent errors, and data analysis options and limitations. Additionally, adoption of NGS-based HIVDR genotyping, especially for clinical care, poses pressing challenges, especially for resource-poor settings, including infrastructure and equipment requirements and cost, logistic and supply chains, instrument service availability, personnel training, validated laboratory protocols, and standardized analysis outputs. The establishment of external quality assessment programs may help to address some of these challenges and is needed to proceed with NGS-based HIVDR genotyping adoption.