Next Article in Journal
Exploring HIV-1 Maturation: A New Frontier in Antiviral Development
Next Article in Special Issue
A Comparison of Sanger Sequencing and Amplicon-Based Next Generation Sequencing Approaches for the Detection of HIV-1 Drug Resistance Mutations
Previous Article in Journal
Topical Protease Inhibitor Increases Tumor-Free and Overall Survival in CD4-Depleted Mouse Model of Anal Cancer
Previous Article in Special Issue
Development and Optimization of Oligonucleotide Ligation Assay (OLA) Probes for Detection of HIV-1 Resistance to Dolutegravir
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Daniele Armenia
1,*,
Luca Carioti
2,
Valeria Micheli
3,
Isabella Bon
4,
Tiziano Allice
5,
Celestino Bonura
6,
Bianca Bruzzone
7,
Fiorenza Bracchitta
3,
Francesco Cerutti
5,
Giovanni Maurizio Giammanco
6,
Federica Stefanelli
7,
Maria Addolorata Bonifacio
8,
Ada Bertoli
9,
Marialinda Vatteroni
10,
Gabriele Ibba
11,
Federica Novazzi
12,
Maria Rosaria Lipsi
13,
Nunzia Cuomo
14,
Ilaria Vicenti
15,
Francesca Ceccherini-Silberstein
2,
Barbara Rossetti
16,
Antonia Bezenchek
17,18,
Francesco Saladini
15,
Maurizio Zazzi
15 and
Maria Mercedes Santoro
2,* on behalf of the Italian HIV NGS network and the ARCA Cohort
add Show full author list remove Hide full author list
1
Departmental Faculty, UniCamillus, Saint Camillus International University of Health Sciences, 00131 Rome, Italy
2
Department of Experimental Medicine, University of Rome “Tor Vergata”, 00133 Rome, Italy
3
Laboratory of Clinical Microbiology, Virology and Bioemergencies, ASST Fatebenefratelli Sacco-University of Milan, 20157 Milan, Italy
4
Microbiology Unit, IRCCS Azienda Ospedaliero-Universitaria di Bologna, 40138 Bologna, Italy
5
Laboratory of Microbiology and Virology, Amedeo di Savoia Hospital, 10149 Turin, Italy
6
Dipartimento di Promozione della Salute, Materno-Infantile, di Medicina Interna e Specialistica di Eccellenza “G. D’Alessandro” (PROSAMI), Azienda Ospedaliera Universitaria Policlinico “P. Giaccone”-University of Palermo, 90127 Palermo, Italy
7
Hygiene Unit, Ospedale Policlinico San Martino, 16132 Genoa, Italy
8
Section of Experimental and Clinical Pathology, Department of Precision and Regenerative Medicine and Jonic Area, University of Bari, 70121 Bari, Italy
9
Virology Unit, Polyclinic of “Tor Vergata”, 00133 Rome, Italy
10
Virology Unit, AOU Pisana, 56126 Pisa, Italy
11
Microbiology and Virology Unit, Diagnostic Department, AOU Sassari, 07100 Sassari, Italy
12
Department of Medicine and Technological Innovation, University of Insubria, 21100 Varese, Italy
13
Microbiology and Virology Unit, Policlinico Riuniti Foggia Hospital, 71121 Foggia, Italy
14
U.O.C. Microbiologia e Virologia, P.O. “D.Cotugno”-AO dei Colli, 80100 Napoli, Italy
15
Department of Medical Biotechnologies, University of Siena, 53100 Siena, Italy
16
Infectious Disease Department, USL SUDEST, Toscana, Misericordia Hospital, 58100 Grosseto, Italy
17
IPRO-InformaPRO S.r.l., 00152 Rome, Italy
18
EuResist Network GEIE, 00152 Rome, Italy
*
Authors to whom correspondence should be addressed.
Viruses 2024, 16(9), 1422; https://doi.org/10.3390/v16091422
Submission received: 7 August 2024 / Revised: 2 September 2024 / Accepted: 3 September 2024 / Published: 6 September 2024
(This article belongs to the Special Issue Antiviral Resistance Mutations)

Abstract

:
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.

1. Introduction

Dramatic advances have been made in HIV management, thanks to the continuous improvement of antiretroviral therapy (ART) and extensive use of resistance testing from the time when treatment starts, and beyond [1,2]. Information about drug resistance helps in selecting more effective antiretroviral regimens and contributing to high rates of virologic success. Genotypic resistance testing (GRT) through Sanger bulk sequencing of plasma HIV RNA has been long and effectively supporting ART. However, this procedure only provides information on the tip of the iceberg, allowing the detection of resistant variants with >15–20% frequency [3,4,5]. The assessment of low-abundance drug-resistance mutations is now possible through Next-Generation Sequencing (NGS) approaches that are becoming affordable in terms of costs and turnaround time [6,7]. Moreover, several studies have shown that NGS-GRT is highly concordant with Sanger sequencing at a 20% threshold [8,9,10,11,12,13,14]. For these reasons, NGS-GRT is gradually replacing Sanger sequencing, but harmonization is lagging behind due to the heterogeneity of currently available lab platforms and NGS data analysis. In this context, virology laboratories are making a lot of effort to adapt their sequencing routine to the new technologies and interpret NGS data to provide an understandable clinical report. However, several issues still need to be addressed including the harmonization of NGS-based HIV drug resistance testing protocols and subsequent data processing and reporting, both of which may benefit from improved automation to minimize artificial errors. While Sanger sequences can be easily generated and analyzed by drug resistance interpretation algorithms such as the Stanford HIV Drug Resistance Database HIVdb algorithm (https://hivdb.stanford.edu/, accessed on 31 January 2024), the harmonization of NGS-GRT is more complex and includes three main steps: (1) wet-lab steps to generate PCR amplicons that cover the pol region and prepare NGS libraries; (2) loading products on NGS platforms; and (3) bioinformatics pipelines, which convert NGS data into user-interpretable HIV drug resistance results [15,16,17]. For the third step, several freely available web tools or standalone software embedded with commercial NGS kits are currently available, providing a straightforward assessment of resistance without any bioinformatics skills [10,18]. Among them, the NGS adaptation of HIVdb and HyDRA Web are well-consolidated free web tools widely used both for research and diagnostic purposes [13,19]. Among standalone tools, SmartVir is provided with the AD4SEQ HIV-1 Solution v2 kit, the most widely used commercial kit for HIV-1 NGS-GRT in Italy (https://www.arrowdiagnostics.it/, accessed on 31 January 2024). To our knowledge, no data on the concordance of NGS-GRT interpretation are available on a large scale. This work aimed at comparing HIVdb, HyDRA Web, and SmartVir for the analysis of NGS-GRT data obtained from a large dataset from Italian real-life settings.

2. Materials and Methods

2.1. Study Design

This is an observational study conducted on NGS data retrieved from 11 Italian virology laboratories, involved in the Italian collaborative HIV NGS Network created within the ARCA cohort (https://db.dbarca.net/, accessed on 31 January 2024) to share information, protocols, and data for improving the harmonization of NGS usage in routine HIV diagnostics. At present, the network connects 36 Italian centers involved in HIV GRT for both diagnostic and research.
All NGS-GRT were performed using plasma HIV-1 RNA samples of viremic people living with HIV (PLWH) processed by using the commercial kit AD4SEQ HIV-1 Solution v2 (Arrow Diagnostics S.r.l., Genova, Italy). According to the documentation provided by the manufacturer, this CE-IVD kit can reliably detect mutations (associated or not with drug resistance) in samples with viremia values >500 copies/mL for the diagnostic routine in different HIV-1 viral subtypes. For each NGS-GRT, information about contextual viremia, treatment status (ART-naïve or -experienced), and the NGS platform used (Illumina MiSeq or iSeq100) were retrieved. NGS data were interpreted through the standalone software provided with the commercial kit (SmartVir) and the two HyDRA Web and HIVdb free online tools.

2.2. RNA Extraction

All blood samples were centrifuged at 2400× g for 15 min at 4 °C for plasma separation. After centrifugation, RNA was extracted using different extraction kits (ELITE InGenius SP200/SP1000 kit, ELITechGroup; QIAamp UltraSens Virus Kit, QIAGEN; NucliSENS® easyMag®/eMAG®, bioMèrieux), according to the manufacturer’s recommendations. The extraction and elution volumes were adjusted according to the viremia levels of samples: <1000 copies/mL (ultracentrifugation for 2 h at 4 °C—input volume of 400 μL−1 mL eluted in 50 μL), 1000–5000 (input volume of 400 μL−1 mL eluted in 50 μL), and >5000 (input volume of 200 μL−1 mL eluted in 50–100 μL).

2.3. Library Preparation for Illumina iSeq100/MiSeq Platforms

The amplicon-based library was prepared according to the manufacturer’s recommendations, allowing us to amplify and sequence protease (PR: 1–99 aa), reverse transcriptase (RT: 1–440 aa), and integrase (INT: 1–288 aa) regions, with a theoretical analytical sensitivity of 500 copies/mL. The final pool was loaded with 10% PhiX control at 9 pM on the MiSeq V2 reagent kit 500-cycle cartridge and MiSeq Reagent Nano kit v2 flow cell (Illumina Inc., San Diego, CA, USA), or at 110 pM on iSeq100 i1 reagent 300-cycle cartridge v2 (Illumina Inc., San Diego, CA, USA).

2.4. Sequence Data Analysis and Coverage Evaluation

The raw FastQ files obtained after sequencing were analyzed using the standalone software SmartVir (SmartSeq s.r.l., Alessandria, Italy) and the web tools HyDRA Web (Ver v1.7.0; available at https://hydra.canada.ca/pages/home, accessed on 31 January 2024) and HIVdb (ver 9.5.1; available at: https://hivdb.stanford.edu/hivdb/by-reads/, accessed on 31 January 2024) set at default parameters for filtering low-quality reads and at 1% cut-off for variant detection. Subtype was assessed according to the HIVdb ver 9.5.1 algorithm. Frequency cut-offs of 5% and 20% were also evaluated. Coverage >100 reads (100×) per position was considered the minimum cut-off for the valid detection of HIV drug resistance as previously reported for HIV NGS-GRT [5,20] and according to manufacturer instructions. Specifically, an NGS-GRT was considered valid for resistance analyses when 100× at each PR/RT/IN resistance-associated position listed in the HIVdb algorithm ver 9.5.1 was obtained. The reliability rate ([number of valid NGS-GRT] × 100/[total NGS-GRT]) was calculated for each tool and for at least two tools. The reliability rate was stratified according to subtype, viremia, platform used, and treatment status. In particular, viremia levels were stratified according to the following strata: ≤1000; 1001–10,000; 10,001–100,000; 100,001–1,000,000; >1,000,000 HIV-1 RNA copies/mL. Predictors of reliability for resistance interpretation were assessed through univariable and multivariable logistic regression.

2.5. Resistance Analysis and Variant Detection Concordance

The flow chart for resistance analysis is reported in Figure 1. All the substitutions detected at resistance-associated positions with a frequency >1% were evaluated. Mutations were stratified according to the following frequency detection groups: (i) frequency ranging from 1% to 5%: low-level minority variants; (ii) frequency ranging from 5% to 20%: minority variants; and (iii) frequency >20%: majority variants. The detection of unusual mutations (defined as having a prevalence <0.01% in HIV-1 group M sequences) according to the HIVdb algorithm or stop codons was considered as part of the assessment of sequencing accuracy. Mutations concordantly detected from all the interpretation tools as low-level minority variants were excluded and not considered valid for resistance interpretation according to already published studies on NGS data [19,21]. To assess the concordance in resistance detection among different interpretation systems, the frequency of each mutation was considered consistent when it fell within the same frequency detection group across the interpretation tools. Concordance was assessed overall and per each position of PR, RT, and IN regions analyzed.

2.6. Statistical Analyses

Differences between dichotomous or categorical variables were assessed using Chi-Squared or Fisher exact tests as appropriate. Differences between continuous variables were determined with the Mann–Whitney or Kruskal–Wallis test as appropriate. Univariable and multivariable logistic regression models were built to assess whether subtype, viremia levels, NGS platform, and treatment status were predictors of coverage reliability. The coefficient of variation (CV), defined as the ratio of the standard deviation to the mean of mutation’s frequency, was calculated to estimate the variability in detecting mutations across interpretation tools. The maximum frequency of mutations with discordant frequency detection group was calculated and the third quartile (3rdQ) of the distributions was considered the threshold of frequency above which misclassifications are less likely to be present.

3. Results

3.1. Samples’ Characteristics and Coverage Reliability Rate

The characteristics of the 629 samples analyzed are reported in Table 1. Most of them had a contextual viremia >1000 copies/mL (92.5%) and were processed through the iSeq100 platform (77.1%). A considerable proportion of samples were from individuals infected with HIV-1 non-B subtypes (41.2%).
The reliability rate was 75.2%, 74.2%, and 70.9 according to SmartVir, HyDRA Web, and HIVdb. Four-hundred and sixty-three (73.6%) NGS-GRT had valid coverage for resistance analyses with at least two interpretation tools. These proportions were affected by HIV-1 subtype and viremia (Figure 2). Namely, the two online interpretation tools were significantly affected by non-B subtypes (Figure 2C,E,G) while subtype G was associated with the lowest rate of coverage reliability with all the interpretation tools. A slight decrease in the coverage reliability rate (below 73%) was observed for all three tools even at plasma HIV-1 RNA levels <10,000 copies/mL (Figure 2B,D,F,H). However, the coverage reliability rate considerably decreased at viremia levels ≤1000 copies/mL, specifically for HyDRA Web and SmartVir tools (Figure 2B,F).
When evaluating valid coverage for resistance interpretation by at least two tools, multivariable logistic regression analyses showed that viremia levels below 10,000 copies /mL (and especially below 1000 copies/mL) and subtypes CRF02_AG, G, and/or heterogeneous recombinant forms were negatively associated with a valid coverage for resistance analyses (Table 2). When each tool was evaluated separately, results varied and only the subtype transversally affected coverage reliability (Table 2).

3.2. Resistance Interpretation

Among 463 GRT with valid coverage in PR/RT/IN for at least two tools, 9596 mutations were detected. Of these, 7483 (78.0%) were low-level minority variants. Concerning unusual mutations, among the 626 mutations identified, the highest proportion (561, 89.6%) was detected as low-level minority variants, followed by discordant minority variants ranging from 5–20% (55, 8.8%) and majority variants (10, 1.6%). They were never found among mutations concordantly detected by at least two tools (0%, p < 0.001).
Overall, excluding low-level minority variants and unusual mutations, 2049 (690 in PR; 768 in RT; 591 in IN) mutations with a frequency >5% for at least one tool were identified, most concordantly detected with a median (IQR) CV of 1.4% (0.5–8.2%). Specifically, 1739 (84.9%) mutations were concordantly detected by at least two tools (majority variants: 64.0%; minority variants: 20.9%). Discordant mutations (N = 310, 15.1%) were mostly minority variants (frequency range 5–20%: 12.4%; frequency >20%: 2.6%) and were more likely observed in the IN region. In fact, the highest median CV was observed among IN mutations (median [IQR] CV IN: 6.1% [1.2–20.4%]; RT: 1.2% [0.4–6.6%]; PR: 0.7% [0.4–2.9%], p < 0.001).
Mutation reports were then evaluated by comparing each tool to the other (Figure 3). Good concordance was observed above a 10% frequency. Discordant minority variants were mainly detected below 10% frequency regardless of the tool used. In general, discordant majority variants were less likely to be observed.

3.3. Concordance in Detecting Mutations at Resistance-Associated Positions of PR/RT/IN

The frequency distribution of mutations detected at each resistance-associated position in the PR/RT/IN regions is reported in Figure 4. In the majority of positions, most mutations were concordantly detected. Concerning discordant mutations, the maximum frequency among tools of minority variants with discordant detection was rarely observed above 10%. In particular, discordant minority variants with a maximum frequency above 10% were found, especially in IN, as follows: (i) 4 spanning PR at positions 10, 47, 48, and 82 (Figure 4A); (ii) 17 spanning RT at positions 100, 101, 103, 106, 108, 179, 225, and 238; and (iii) 28 spanning IN at position 66, 74, 95, 118, 121, 138, 140, 147, 151, 153, 157, and 232 (Figure 4B–D). In general, discordant majority variants (orange crosses) were more rarely detected.

4. Discussion

Considering the ongoing transition from Sanger to NGS technologies for routine HIV-1 resistance testing, the collection and evaluation of information retrieved from NGS in real life is crucial. In the present study, for the first time at a nationwide level in Italy, a considerable amount of NGS data (including several HIV-1 non-B subtypes) retrieved from real life was collected and explored. Among the 629 NGS-GRT tested with three different interpretation tools, around 70–75% had valid coverage (at least 100×) to evaluate drug resistance mutations in the HIV-1 PR/RT/IN regions. The coverage rate observed in this study was similar to or even higher than that observed in other studies evaluating in-house or commercial kits [14,20,22].
A key question is whether around 75% of NGS-GRT success might be satisfactory in the clinical setting. Notably, certain HIV-1 subtypes and moderate-low viremia levels may affect this success rate. In fact, we found that subtype G and heterogeneous recombinant forms were associated with the lowest rate (<60%) of valid coverage regardless of the interpretation tool used. This was not observed in previous studies evaluating in-house or commercial NGS methods [20,23], likely because of the lower HIV-1 subtype heterogeneity among samples tested. Indeed, in the present study, more than 40% of clinical samples had non-B subtype strains, including around 10% of heterogeneous CRF. Thus, technical adjustments should be implemented to reduce the impact of genetic variability on the coverage issue.
Concerning viremia levels, a slight decrease in coverage reliability was generally observed at levels below 10,000 copies/mL, thus raising the concern about the proper identification of resistance mutations at viral loads commonly linked to viral failure during antiviral therapy and an increased likelihood of drug resistance emergence [24,25,26].
Nevertheless, a considerable decrease in the coverage reliability rate was observed for samples with contextual viremia below 1000 copies/mL and this has been emphasized by the inclusion of several samples with viremia levels below 500 copies/mL (median [IQR]: 276 [187–510] copies/mL), which is the analytical sensitivity threshold reported for the kit used in this study.
Taken together, these results confirm that NGS for routine HIV GRT is affordable even if there are still some issues related to HIV-1 non-B subtypes and low viral loads. Unfortunately, these limitations are clinically relevant given the increasing number of non-B subtypes in many geographic areas (including Italy) [27,28,29] and the frequent GRT requests for individuals with low-level viremia [30].
Another important aspect evaluated in this study is the concordance in detecting mutations using different NGS data processing tools. This is an important point considering that virologists may use multiple interpretation tools for making resistance reports, especially to verify the presence of resistance mutations as minority variants. Several studies suggest that minority variants with a frequency below 3–5% may result from sequencing artifacts [19,21]. Indeed, we observed the highest proportion of mutations recognized as unusual and/or affected by discordant detection at a frequency in the range of 1 to 5%, supporting the technical cut-off for HIV NGS-GRT at 5%. This finding is in line with the manufacturer’s recommendations about the frequency cut-off to be considered. However, at a frequency >5%, variation among tools in detecting minority variants remained quite common, affecting around 15% of the total number of mutations detected, particularly in the IN region. In fact, the highest median CV was observed for IN mutations and discordant mutations were observed widely across the IN region. Again, this is relevant considering that INSTI are the most currently used drugs, recommended as part of first-line regimens from guidelines [1,2], and uncertainty in detecting resistance minority variants might become challenging when making GRT reports. Despite this, most of the discordant mutations were generally detected below 10%. While this is important information to guide resistance interpretation among virology labs that use commercial kits, differences may arise with other NGS and analysis pipelines. For example, another study evaluating five NGS HIV drug resistance pipelines [17] concluded that the 2% threshold may be appropriate to detect resistance mutations with high specificity. Given the wide heterogeneity in wet-lab procedures, platforms, and interpretation systems, a universal technical cut-off may remain challenging to set.
This study might have some limitations. First of all, the real performance of the whole NGS-GRT system was not evaluated because we compared only the NGS data interpretation tools. Since all the samples evaluated were processed with the same commercial kit for NGS, the technical cut-offs proposed cannot be extrapolated to other commercial or in-house approaches. Finally, no information about the clinical role of minority mutations can be concluded considering the mere technical scope of the present study.

5. Conclusions

In conclusion, the usage of NGS-GRT in Italy is increasing, mostly based on the availability of a standardized kit. However, technical improvements are required. At present, a technical frequency cut-off of 10% should be advisable to process the NGS-GRT data and guide a cautious transition from the well-established Sanger sequencing to NGS in the clinical setting. Clearly, given the extraordinary pace of technical progress in the field, recommendations are subject to frequent updates to ensure that clinical applications benefit from the latest technology.

Author Contributions

Conceptualization, D.A. and M.M.S.; methodology, D.A.; software, D.A. and L.C.; validation, D.A.; formal analysis, D.A.; investigation, D.A. and M.M.S.; resources, V.M., I.B., T.A., C.B., B.B., F.B., F.C., G.M.G., F.S. (Federica Stefanelli), M.A.B., A.B. (Ada Bertoli), M.V., G.I., F.N., M.R.L., N.C., I.V. and F.S. (Francesco Saladini); data curation, D.A., L.C. and A.B. (Antonia Bezenchek); writing—original draft preparation, D.A. and M.M.S.; writing—review and editing, M.Z., V.M. and I.B.; visualization L.C., V.M., I.B., T.A., C.B., B.B., F.B., F.C., G.M.G., F.S. (Federica Stefanelli), M.A.B., A.B. (Ada Bertoli), M.V., G.I., F.N., M.R.L., N.C., I.V., F.C.-S., B.R. and F.S. (Francesco Saladini); supervision, F.C.-S., M.Z. and M.M.S.; project administration, M.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

No specific funding was provided for this study. The current work was a joint activity of the Italian HIV NGS Network and the ARCA database cohort. The ARCA Database and initiatives are supported by unconditional and educational grants from Arrow Diagnostics, MSD, Gilead Sciences, Janssen Infectious Diseases & Vaccines, and ViiV Healthcare.

Institutional Review Board Statement

The ARCA database was approved by the central Ethics Committee (COMITATO ETICO REGIONE TOSCANA-AREA VASTA SUD EST-C.E.A.V.S.E.-21 July 2014; approval code: progetto ARCA) and by the Ethics Committees of the individual clinic centers. The study was performed in accordance with the ethical principles of the Declaration of Helsinki and the Good Clinical Practice guidelines of the International Conference on Harmonization.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author because they contain personal sensitive information subject to GDPR; the related privacy and ethical issues are managed by following the ARCA Board rules on sensitive data sharing for research purposes.

Acknowledgments

A part of this work was presented at the 21st European Meeting on HIV & Hepatitis Congress, Rome, Italy, 7–9 June 2023 (Abstract 8). We wish to thank all the clinicians, virologists, and data managers throughout Italy who contributed with their work to develop, expand, and maintain the Italian HIV NGS network and the ARCA Cohort. We also thank Debra Mandatori for having revised and edited the manuscript. The following centers and people are involved in the Italian HIV NGS network: Hygiene Unit, Ospedale Policlinico San Martino, Genoa: Bianca Bruzzone, Federica Stefanelli, Nadia Randazzo; Laboratory of Microbiology and Virology, Amedeo di Savoia Hospital, Turin: Tiziano Allice, Francesco Cerutti; Immunopathology and Cancer Biomarkers, National Cancer Institute, Centro di Riferimento Oncologico, Istituto di Ricovero e Cura a Carattere Scientifico di Aviano, Aviano: Stefania Zanussi; Department of Medicine and Technological innovation, University of Insubria, Varese: Nicasio Mancini, Federica Novazzi; Clinical Immunology, Allergy and Advanced Biotechnologies Unit, Azienda USL-IRCCS, Reggio Emilia: Chiara Carone, Romina Corsini, Ilaria Giovanelli, Alessandro Zerbini; Department of Medicine and Surgery, University of Parma, Parma: Mirko Buttrini, Flora De Conto; Molecular Microbiology and Virology Unit, Azienda Ospedaliero Universitaria di Modena, Modena: William Gennari; Microbiology Unit, Hospital of Legnano, Legnano: Pierangelo Clerici, Massimo De Paschale; Infectious Diseases Unit, Hospital of Legnano, Legnano: Marco Franzetti; Information Technology Unit, Hospital of Legnano, Legnano: Marco Andriolo; S.C. Microbiologia, Fondazione IRCCS San Gerardo dei Tintori, Monza: Sergio Malandrin, Monica Tallarita; Microbiology and Virology Unit, ASST Papa Giovanni XXIII, Bergamo: Marco Enrico Giovanni Arosio, Claudio Francesco Farina, Davide Guarneri; Department of Biomedical and Clinical Sciences Luigi Sacco, University of Milan, Milan: Alessia Lai; Laboratory of Clinical Microbiology, Virology and Bioemergencies, ASST Fatebenefratelli Sacco-University of Milan, Milan: Fiorenza Bracchitta, Federica De Poli, Valeria Micheli; Department of Oncology and Hemato-Oncology, University of Milan, Milan: Claudia Alteri; Department of Health Sciences, University of Milan, Milan: Elisa Borghi; Laboratory of Microbiology and Virology, IRCCS San Raffaele Scientific Institute, Milan: Nicola Clementi; Complex structure of Microbiology and Virology, Policlinic of Milan, Milan: Anna Paola Callegaro; Microbiology Unit, IRCCS Azienda Ospedaliero-Universitaria di Bologna: Alessia Bertoldi, Isabella Bon; Infectious Disease Department, USL SUDEST, Toscana, Misericordia Hospital, Grosseto: Barbara Rossetti; Microbiology and Virology Unit, Careggi University Hospital, Florence: Maria Grazia Colao, Paola Pace, Seble Tekle; Virology Unit, AOU Pisana, Pisa: Maria Linda Vatteroni; Department of Medical Biotechnologies, University of Siena, Siena: Francesco Saladini, Ilaria Vicenti, Maurizio Zazzi; IPRO-InformaPRO S.r.l.-Rome EuResist Network GEIE, Rome: Antonia Bezenchek; Virology Unit, National Institute for Infectious Diseases L Spallanzani, IRCCS, Rome: Giulia Berno, Lavinia Fabeni, Federica Forbici, Fabrizio Maggi, Martina Rueca; Institute of Clinical Infectious Diseases, Fondazione Policlinico Universitario A. Gemelli IRCCS: Francesca Lombardi; Department of Laboratory and Hematological Sciences, Fondazione Policlinico Universitario A. Gemelli IRCCS: Maria Rosaria Santangelo; Department of Molecular Medicine, Sapienza University of Rome: Guido Antonelli, Ombretta Turriziani; Microbiology and Diagnostics of Immunology Unit, Bambino Gesù Children Hospital IRCCS, Rome: Luna Colagrossi, Carlo Federico Perno, Rossana Scutari; UniCamillus, Saint Camillus International University of Health Sciences, Rome: Daniele Armenia, Bouba Yagai; Department of Experimental Medicine, University of Rome Tor Vergata, Rome: Maria Concetta Bellocchi, Luca Carioti, Francesca Ceccherini Silberstein, Eizadi Moghadam Hossein, Omar El Khalili, Greta Marchegiani, Maria Santoro, Daniele Spalletta; Virology Unit, Polyclinic of Tor Vergata, Rome: Sergio Bernardini, Ada Bertoli, Sandro Grelli, Vita Petrone, Lorenzo Piermatteo, Giulia Torre; U.O.C. Microbiologia e Virologia, P.O. “D.Cotugno”-AO dei Colli–Napoli: Nunzia Cuomo, Adriana Raddi; Section of Experimental and Clinical Pathology, Department of Precision and Regenerative Medicine and Jonic Area, University of Bari, Bari: Maria Addolorata Bonifacio, Maria Addolorata Mariggiò, Anna Volpe; Microbiology and Virology Unit, Policlinico Riuniti Foggia Hospital, Foggia: Maria Rosaria Lipsi; Microbiology and Virology Unit, Diagnostic Department, AOU Sassari, Sassari: Flavia Angioj, Gabriele Ibba, Sergio Uzzau; Azienda Ospedaliera Universitaria Policlinico “P. Giaccone”-University of Palermo, Palermo: Celestino Bonura, Sara Cannella, Giovanni M. Giammanco, Chiara Mascarella; Unit of Infectious Diseases, School of Medicine and Surgery, Kore University of Enna, Enna: Manuela Ceccarelli; Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina: Teresa Pollicino.

Conflicts of Interest

The authors declare no conflicts of interest with respect to the research, authorship, and/or publication of this article.

References

  1. European AIDS Clinical Society. Guidelines—Version 12.0. 2023. Available online: https://www.eacsociety.org/media/guidelines-12.0.pdf (accessed on 31 July 2024).
  2. Panel on Antiretroviral Guidelines for Adults and Adolescents, Guidelines for the Use of Antiretroviral Agents in Adults and Adolescents with HIV, Department of Health and Human Services. 2024. Available online: https://clinicalinfo.hiv.gov/sites/default/files/guidelines/documents/adult-adolescent-arv/guidelines-adult-adolescent-arv.pdf (accessed on 31 July 2024).
  3. Parkin, N.T.; Avila-Rios, S.; Bibby, D.F.; Brumme, C.J.; Eshleman, S.H.; Harrigan, P.R.; Howison, M.; Hunt, G.; Ji, H.; Kantor, R.; et al. Multi-Laboratory Comparison of Next-Generation to Sanger-Based Sequencing for HIV-1 Drug Resistance Genotyping. Viruses 2020, 12, 694. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, N.Y.; Kao, S.W.; Liu, Z.H.; Wu, T.S.; Tsai, C.L.; Lin, H.H.; Wong, W.W.; Chang, Y.Y.; Chen, S.S.; Ku, S.W.W. Shall I Trust the Report? Variable Performance of Sanger Sequencing Revealed by Deep Sequencing on HIV Drug Resistance Mutation Detection. Int. J. Infect. Dis. 2020, 93, 182–191. [Google Scholar] [CrossRef]
  5. Labate, L.; Bruzzone, B.; Spagnuolo, V.; Zazzi, M.; Santoro, M.; di Biagio, A.; Castagna, A. PRESTIGIO RING: “A 59-year-old HIV-1 positive, highly treatment-experienced woman failing darunavir/ritonavir plus raltegravir”. New Microbiol. 2023, 46, 226–230. [Google Scholar] [CrossRef] [PubMed]
  6. Chabria, S.B.; Gupta, S.; Kozal, M.J. Deep Sequencing of HIV: Clinical and Research Applications. Annu. Rev. Genom. Hum. Genet. 2014, 15, 295–325. [Google Scholar] [CrossRef]
  7. Ávila-Ríos, S.; Parkin, N.; Swanstrom, R.; Paredes, R.; Shafer, R.; Ji, H.; Kantor, R. Next-Generation Sequencing for HIV Drug Resistance Testing: Laboratory, Clinical, and Implementation Considerations. Viruses 2020, 12, 617. [Google Scholar] [CrossRef] [PubMed]
  8. Tzou, P.L.; Ariyaratne, P.; Varghese, V.; Lee, C.; Rakhmanaliev, E.; Villy, C.; Yee, M.; Tan, K.; Michel, G.; Pinsky, B.A.; et al. Comparison of an In Vitro Diagnostic Next-Generation Sequencing Assay with Sanger Sequencing for HIV-1 Genotypic Resistance Testing. J. Clin. Microbiol. 2018, 56, e00105-18. [Google Scholar] [CrossRef] [PubMed]
  9. Alidjinou, E.K.; Deldalle, J.; Hallaert, C.; Robineau, O.; Ajana, F.; Choisy, P.; Hober, D.; Bocket, L. RNA and DNA Sanger Sequencing versus Next-Generation Sequencing for HIV-1 Drug Resistance Testing in Treatment-Naive Patients. J. Antimicrob. Chemother. 2017, 72, 2823–2830. [Google Scholar] [CrossRef]
  10. Nicot, F.; Jeanne, N.; Raymond, S.; Delfour, O.; Carcenac, R.; Lefebvre, C.; Sauné, K.; Delobel, P.; Izopet, J. Performance Comparison of Deep Sequencing Platforms for Detecting HIV-1 Variants in the Pol Gene. J. Med. Virol. 2018, 90, 1486–1492. [Google Scholar] [CrossRef]
  11. Dudley, D.M.; Bailey, A.L.; Mehta, S.H.; Hughes, A.L.; Kirk, G.D.; Westergaard, R.P.; O’Connor, D.H. Cross-Clade Simultaneous HIV Drug Resistance Genotyping for Reverse Transcriptase, Protease, and Integrase Inhibitor Mutations by Illumina MiSeq. Retrovirology 2014, 11, 122. [Google Scholar] [CrossRef]
  12. Ram, D.; Leshkowitz, D.; Gonzalez, D.; Forer, R.; Levy, I.; Chowers, M.; Lorber, M.; Hindiyeh, M.; Mendelson, E.; Mor, O. Evaluation of GS Junior and MiSeq Next-Generation Sequencing Technologies as an Alternative to Trugene Population Sequencing in the Clinical HIV Laboratory. J. Virol. Methods 2015, 212, 12–16. [Google Scholar] [CrossRef]
  13. Taylor, T.; Lee, E.R.; Nykoluk, M.; Enns, E.; Liang, B.; Capina, R.; Gauthier, M.K.; Van Domselaar, G.; Sandstrom, P.; Brooks, J.; et al. A MiSeq-HyDRA Platform for Enhanced HIV Drug Resistance Genotyping and Surveillance. Sci. Rep. 2019, 9, 8970. [Google Scholar] [CrossRef] [PubMed]
  14. Bonifacio, M.A.; Genchi, C.; Lagioia, A.; Talamo, V.; Volpe, A.; Mariggiò, M.A. Analytical Assessment of the Vela Diagnostics NGS Assay for HIV Genotyping and Resistance Testing: The Apulian Experience. Int. J. Mol. Sci. 2022, 23, 2727. [Google Scholar] [CrossRef]
  15. Brumme, C.J.; Poon, A.F.Y. Promises and Pitfalls of Illumina Sequencing for HIV Resistance Genotyping. Virus Res. 2017, 239, 97–105. [Google Scholar] [CrossRef] [PubMed]
  16. Casadellà, M.; Paredes, R. Deep Sequencing for HIV-1 Clinical Management. Virus Res. 2017, 239, 69–81. [Google Scholar] [CrossRef]
  17. Lee, E.R.; Parkin, N.; Jennings, C.; Brumme, C.J.; Enns, E.; Casadellà, M.; Howison, M.; Coetzer, M.; Avila-Rios, S.; Capina, R.; et al. Performance Comparison of next Generation Sequencing Analysis Pipelines for HIV-1 Drug Resistance Testing. Sci. Rep. 2020, 10, 1634. [Google Scholar] [CrossRef]
  18. Noguera-Julian, M.; Lee, E.R.; Shafer, R.W.; Kantor, R.; Ji, H. Dry Panels Supporting External Quality Assessment Programs for Next Generation Sequencing-Based HIV Drug Resistance Testing. Viruses 2020, 12, 666. [Google Scholar] [CrossRef] [PubMed]
  19. Tzou, P.L.; Kosakovsky Pond, S.L.; Avila-Rios, S.; Holmes, S.P.; Kantor, R.; Shafer, R.W. Analysis of Unusual and Signature APOBEC-Mutations in HIV-1 Pol Next-Generation Sequences. PLoS ONE 2020, 15, e0225352. [Google Scholar] [CrossRef]
  20. Pyne, M.T.; Simmon, K.E.; Mallory, M.A.; Hymas, W.C.; Stevenson, J.; Barker, A.P.; Hillyard, D.R. HIV-1 Drug Resistance Assay Using Ion Torrent Next Generation Sequencing and On-Instrument End-to-End Analysis Software. J. Clin. Microbiol. 2022, 60, e00253-22. [Google Scholar] [CrossRef]
  21. Balakrishna, S.; Loosli, T.; Zaheri, M.; Frischknecht, P.; Huber, M.; Kusejko, K.; Yerly, S.; Leuzinger, K.; Perreau, M.; Ramette, A.; et al. Frequency Matters: Comparison of Drug Resistance Mutation Detection by Sanger and Next-Generation Sequencing in HIV-1. J. Antimicrob. Chemother. 2023, 78, 656–664. [Google Scholar] [CrossRef]
  22. Weber, J.; Volkova, I.; Sahoo, M.K.; Tzou, P.L.; Shafer, R.W.; Pinsky, B.A. Prospective Evaluation of the Vela Diagnostics Next-Generation Sequencing Platform for HIV-1 Genotypic Resistance Testing. J. Mol. Diagn. 2019, 21, 961. [Google Scholar] [CrossRef]
  23. Mohamed, S.; Boulmé, R.; Sayada, C. From Capillary Electrophoresis to Deep Sequencing: An Improved HIV-1 Drug Resistance Assessment Solution Using In Vitro Diagnostic (IVD) Assays and Software. Viruses 2023, 15, 571. [Google Scholar] [CrossRef] [PubMed]
  24. Gonzalez-Serna, A.; Min, J.E.; Woods, C.; Chan, D.; Lima, V.D.; Montaner, J.S.; Harrigan, P.R.; Swenson, L.C. Performance of HIV-1 drug resistance testing at low-level viremia and its ability to predict future virologic outcomes and viral evolution in treatment-naive individuals. Clin. Infect. Dis. 2014, 58, 1165–1173. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  25. Santoro, M.M.; Fabeni, L.; Armenia, D.; Alteri, C.; Di Pinto, D.; Forbici, F.; Bertoli, A.; Di Carlo, D.; Gori, C.; Carta, S.; et al. Reliability and clinical relevance of the HIV-1 drug resistance test in patients with low viremia levels. Clin. Infect. Dis. 2014, 58, 1156–1164. [Google Scholar] [CrossRef] [PubMed]
  26. Swenson, L.C.; Min, J.E.; Woods, C.K.; Cai, E.; Li, J.Z.; Montaner, J.S.; Harrigan, P.R.; Gonzalez-Serna, A. HIV drug resistance detected during low-level viraemia is associated with subsequent virologic failure. AIDS 2014, 28, 1125–1134. [Google Scholar] [CrossRef] [PubMed] [PubMed Central]
  27. Williams, A.; Menon, S.; Crowe, M.; Agarwal, N.; Biccler, J.; Bbosa, N.; Ssemwanga, D.; Adungo, F.; Moecklinghoff, C.; Macartney, M.; et al. Geographic and Population Distributions of Human Immunodeficiency Virus (HIV)-1 and HIV-2 Circulating Subtypes: A Systematic Literature Review and Meta-Analysis (2010–2021). J. Infect. Dis. 2023, 228, 1583–1591. [Google Scholar] [CrossRef]
  28. Fabeni, L.; Armenia, D.; Abbate, I.; Gagliardini, R.; Mazzotta, V.; Bertoli, A.; Gennari, W.; Forbici, F.; Berno, G.; Piermatteo, L.; et al. HIV-1 Transmitted Drug Resistance in Newly Diagnosed Individuals in Italy over the Period 2015–21. J. Antimicrob. Chemother. 2024, dkae189. [Google Scholar] [CrossRef]
  29. Bbosa, N.; Kaleebu, P.; Ssemwanga, D. HIV Subtype Diversity Worldwide. Curr. Opin. HIV AIDS 2019, 14, 153–160. [Google Scholar] [CrossRef]
  30. Wu, F.; Simonetti, F.R. Learning from Persistent Viremia: Mechanisms and Implications for Clinical Care and HIV-1 Cure. Curr. HIV/AIDS Rep. 2023, 20, 428–439. [Google Scholar] [CrossRef]
Figure 1. Flow chart for the evaluation of concordance among tools for HIV-1 NGS data interpretation in resistance assessment.
Figure 1. Flow chart for the evaluation of concordance among tools for HIV-1 NGS data interpretation in resistance assessment.
Viruses 16 01422 g001
Figure 2. Proportion of samples with valid coverage for resistance evaluation obtained through different NGS data interpretation tools, according to subtypes and contextual viremia levels. Bars represent the proportion of samples with 100× coverage among all PR/RT/IN positions associated with resistance according to the Stanford drug resistance algorithm (HIVdb 9.5.0). Black and white bars represent proportions according to subtypes and viremia levels, respectively. (A,B) SmartVir standalone tool; (C,D) free web tool HIVdb; (E,F) free web tool HyDRA Web; (G,H) at least 2 tools among SmartVir, HyDRA Web, and HIVdb.
Figure 2. Proportion of samples with valid coverage for resistance evaluation obtained through different NGS data interpretation tools, according to subtypes and contextual viremia levels. Bars represent the proportion of samples with 100× coverage among all PR/RT/IN positions associated with resistance according to the Stanford drug resistance algorithm (HIVdb 9.5.0). Black and white bars represent proportions according to subtypes and viremia levels, respectively. (A,B) SmartVir standalone tool; (C,D) free web tool HIVdb; (E,F) free web tool HyDRA Web; (G,H) at least 2 tools among SmartVir, HyDRA Web, and HIVdb.
Viruses 16 01422 g002
Figure 3. Concordance in detecting mutations between different NGS resistance interpretation tools. (A) Comparison between SmartVir and HIVdb; (B) comparison between SmartVir and HyDRA Web; (C) comparison between HyDRA Web and HIVdb. Each symbol represents a substitution detected with a valid coverage per each tool. Blue and green symbols represent mutations concordantly detected as majority (frequency >20%) and minority variants (frequency 1–5%), respectively. Red crosses represent mutations discordantly detected as low-level minority mutations (frequency 1–5%) or minority mutations (5–20%). Orange crosses represent mutations discordantly detected as minority mutations (frequency 5–20%) or majority mutations (5–20%). In the table on the right of scatter plots, the distribution of maximum frequency of discordant mutations detected between tools is reported; the third quartile of the distribution is highlighted in bold. Dotted lines represent 10% of frequency cut-off.
Figure 3. Concordance in detecting mutations between different NGS resistance interpretation tools. (A) Comparison between SmartVir and HIVdb; (B) comparison between SmartVir and HyDRA Web; (C) comparison between HyDRA Web and HIVdb. Each symbol represents a substitution detected with a valid coverage per each tool. Blue and green symbols represent mutations concordantly detected as majority (frequency >20%) and minority variants (frequency 1–5%), respectively. Red crosses represent mutations discordantly detected as low-level minority mutations (frequency 1–5%) or minority mutations (5–20%). Orange crosses represent mutations discordantly detected as minority mutations (frequency 5–20%) or majority mutations (5–20%). In the table on the right of scatter plots, the distribution of maximum frequency of discordant mutations detected between tools is reported; the third quartile of the distribution is highlighted in bold. Dotted lines represent 10% of frequency cut-off.
Viruses 16 01422 g003
Figure 4. Concordance in detecting mutations among at least two NGS resistance interpretation tools at each position associated with resistance spanning PR/RT/IN. (A) Protease positions associated with resistance to PI. (B) Reverse transcriptase positions associated with resistance to NRTI. (C) Reverse transcriptase positions associated with resistance to NNRTI. (D) Integrase positions associated with resistance to INSTI. X axis represents the positions associated with drug resistance according to the Stanford drug resistance algorithm (HIVdb 9.5.0); Y axis represents the maximum frequency detected in the reports retrieved from the three interpretation tools. Each symbol represents a substitution detected with a valid coverage in at least two tools. Blue and green symbols represent mutations concordantly detected as majority (frequency > 20%) and minority variants (frequency 1–5%), respectively. Red crosses represent mutations discordantly detected as low-level minority mutations (frequency 1–5%) or minority mutations (5–20%). Orange crosses represent mutations discordantly detected as minority mutations (frequency 5–20%) or majority mutations (5–20%).
Figure 4. Concordance in detecting mutations among at least two NGS resistance interpretation tools at each position associated with resistance spanning PR/RT/IN. (A) Protease positions associated with resistance to PI. (B) Reverse transcriptase positions associated with resistance to NRTI. (C) Reverse transcriptase positions associated with resistance to NNRTI. (D) Integrase positions associated with resistance to INSTI. X axis represents the positions associated with drug resistance according to the Stanford drug resistance algorithm (HIVdb 9.5.0); Y axis represents the maximum frequency detected in the reports retrieved from the three interpretation tools. Each symbol represents a substitution detected with a valid coverage in at least two tools. Blue and green symbols represent mutations concordantly detected as majority (frequency > 20%) and minority variants (frequency 1–5%), respectively. Red crosses represent mutations discordantly detected as low-level minority mutations (frequency 1–5%) or minority mutations (5–20%). Orange crosses represent mutations discordantly detected as minority mutations (frequency 5–20%) or majority mutations (5–20%).
Viruses 16 01422 g004
Table 1. Sample characteristics.
Table 1. Sample characteristics.
CharacteristicsOverall (n = 629)
Age, years, Median (IQR)44 (34–53)
Year of genotyping, n (%)
<202028 (4.4)
2020–2021104 (16.5)
2022–2023497 (79.1)
Plasma HIV-1 RNA, log (copies/mL), Median (IQR)4.9 (4.2–5.6)
≤1000 copies/mL, n (%)41 (6.5)
1001–10,000 copies/mL, n (%)90 (14.3)
10,001–100,000 copies/mL, n (%)207 (32.9)
100,001–1,000,000 copies/mL, n (%)212 (33.7)
>1,000,000 copies/mL, n (%)79 (12.6)
HIV-1Subtype, n (%) 1
B370 (58.8)
CRF02_AG88 (14.0)
F37 (5.9)
C36 (5.7)
A38 (6.0)
G19 (3.8)
Others 241 (6.5)
Treatment status, n (%)
Naive313 (49.8)
Experienced152 (24.2)
Unknown164 (26.1)
Platform used for NGS
iSeq100485 (77.1)
MiSeq144 (22.9)
1 Subtype assessed through HIVdb algorithm 9.5.1. 2 Other subtypes: B + F (N = 10), CRF01_AE (N = 9), B + G (N = 5), CRF18_cpx (N = 4), CRF06_cpx (N = 2), CRF09_cpx (N = 2), CRF25_cpx (N = 2), D (N = 2), A + G (N = 1), B + C (N = 1), B + CRF01_AE (N = 1), CRF12_BF (N = 1), CRF29_BF (N = 1).
Table 2. Factors associated with coverage reliability rate for resistance interpretation of NGS-GRT performed for clinical routine.
Table 2. Factors associated with coverage reliability rate for resistance interpretation of NGS-GRT performed for clinical routine.
VariablesOdd Ratio to Have Reliable Coverage for Resistance Interpretation (100X Coverage at All PR/RT/IN Resistance Positions)
At Least 2 ToolsSmartVirHIVdbHydra
CrudeAdjusted 1CrudeAdjusted 1CrudeAdjusted 1CrudeAdjusted 1
OR
(95% C.I.)
p ValueOR
(95% C.I.)
p ValueOR
(95% C.I.)
p ValueOR
(95% C.I.)
p ValueOR
(95% C.I.)
p ValueOR
(95% C.I.)
p ValueOR
(95% C.I.)
p ValueOR
(95% C.I.)
p Value
Subtype
B 21 1 1 1 1 1 1
CRF02_AG0.4 (0.2–0.7)<0.0010.4 (0.2–0.6)<0.0010.7 (0.4–1.1)0.1320.6 (0.3–1.0)0.0510.4 (0.2–0.6)<0.0010.4 (0.2–0.6)<0.0010.5 (0.3–0.8)0.0020.4 (0.2–0.7)0.001
F0.6 (0.3–1.3)0.1690.5 (0.3–1.2)0.111.1 (0.5–2.5)0.8521.1 (0.5–2.5)0.8920.5 (0.2–0.9)0.0320.5 (0.2–1)0.0430.7 (0.3–1.5)0.3810.7 (0.3–1.4)0.293
A0.8 (0.3–1.7)0.4790.7 (0.3–1.6)0.4141.5 (0.6–3.7)0.3891.3 (0.5–3.4)0.5410.3 (0.2–0.7)0.0030.3 (0.2–0.7)0.0030.9 (0.4–2.1)0.8430.9 (0.4–2)0.764
C0.5 (0.2–1)0.0450.5 (0.2–1)0.0521.1 (0.5–2.5)0.7881.1 (0.5–2.7)0.7520.4 (0.2–0.8)0.0060.4 (0.2–0.8)0.0130.5 (0.2–0.9)0.0270.4 (0.2–0.9)0.025
G0.3 (0.1–0.7)0.0070.3 (0.1–0.7)0.0050.4 (0.2–1.1)0.0640.4 (0.1–1.0)0.0400.2 (0.1–0.6)0.0020.2 (0.1–0.5)0.0010.3 (0.1–0.7)0.0100.3 (0.1–0.7)0.009
Others0.4 (0.2–0.8)0.0070.3 (0.2–0.7)0.0020.6 (0.3–1.1)0.1160.5 (0.2–0.9)0.0280.2 (0.1–0.5)<0.0010.2 (0.1–0.5)<0.0010.6 (0.3–1.3)0.2150.5 (0.2–1.1)0.079
Platform used
Illumina Iseq100 21 1 1 1 1 1
Illumina Miseq1.5 (0.9–2.3)0.0861.6 (1–2.5)0.0650.9 (0.6–1.4)0.615--2.1 (1.3–3.3)0.0022.1 (1.2–3.6)0.0081.2 (0.8–1.9)0.376--
Viremia (copies/mL)
≤10000.5 (0.3–1)0.0580.4 (0.2–0.7)0.0050.3 (0.1–0.6)<0.0010.3 (0.1–0.5)<0.0010.7 (0.4–1.5)0.379--0.4 (0.2–0.9)0.0160.3 (0.2–0.7)0.004
1001–10,0000.6 (0.3–1)0.0520.5 (0.3–0.9)0.0140.8 (0.4–1.3)0.3410.7 (0.4–1.3)0.2790.7 (0.4–1.1)0.108--0.5 (0.3–0.9)0.0270.5 (0.3–0.9)0.012
10,001–100,0001.1 (0.7–1.8)0.5771.1 (0.7–1.7)0.8211.2 (0.7–1.9)0.4781.2 (0.7–1.9)0.5591.1 (0.7–1.6)0.769--1.1 (0.7–1.8)0.6471 (0.6–1.7)0.901
100,001–1,000,000 21 1 1 1 1 1 1
>1,000,0001 (0.5–1.7)0.891 (0.5–1.8)0.8740.8 (0.4–1.5)0.4830.8 (0.4–1.4)0.3830.9 (0.5–1.5)0.6110.9 (0.5–1.6)0.6940.9 (0.5–1.6)0.629
Treatment status
cART naïve 21 1 1 1
cART experienced1.3 (0.8–2)0.249--1.0 (0.6–1.5)0.8731.1 (0.7–1.9)0.6221.7 (1.1–2.7)0.0151.4 (0.9–2.2)0.1921.2 (0.8–1.9)0.377--
Unknown1.3 (0.8–1.9)0.31--0.6 (0.4–1.0)0.0380.7 (0.4–1.1)0.1011.6 (1.0–2.4)0.0400.9 (0.6–1.6)0.8291 (0.7–1.6)0.885--
1 Variables with p < 0.1 at univariable analyses were included in multivariable models. 2 Reference (dummy). Reported variables significantly associated with reliable coverage are in bold.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Armenia, D.; Carioti, L.; Micheli, V.; Bon, I.; Allice, T.; Bonura, C.; Bruzzone, B.; Bracchitta, F.; Cerutti, F.; Giammanco, G.M.; et al. Comparison of Different HIV-1 Resistance Interpretation Tools for Next-Generation Sequencing in Italy. Viruses 2024, 16, 1422. https://doi.org/10.3390/v16091422

AMA Style

Armenia D, Carioti L, Micheli V, Bon I, Allice T, Bonura C, Bruzzone B, Bracchitta F, Cerutti F, Giammanco GM, et al. Comparison of Different HIV-1 Resistance Interpretation Tools for Next-Generation Sequencing in Italy. Viruses. 2024; 16(9):1422. https://doi.org/10.3390/v16091422

Chicago/Turabian Style

Armenia, Daniele, Luca Carioti, Valeria Micheli, Isabella Bon, Tiziano Allice, Celestino Bonura, Bianca Bruzzone, Fiorenza Bracchitta, Francesco Cerutti, Giovanni Maurizio Giammanco, and et al. 2024. "Comparison of Different HIV-1 Resistance Interpretation Tools for Next-Generation Sequencing in Italy" Viruses 16, no. 9: 1422. https://doi.org/10.3390/v16091422

APA Style

Armenia, D., Carioti, L., Micheli, V., Bon, I., Allice, T., Bonura, C., Bruzzone, B., Bracchitta, F., Cerutti, F., Giammanco, G. M., Stefanelli, F., Bonifacio, M. A., Bertoli, A., Vatteroni, M., Ibba, G., Novazzi, F., Lipsi, M. R., Cuomo, N., Vicenti, I., ... Santoro, M. M., on behalf of the Italian HIV NGS network and the ARCA Cohort. (2024). Comparison of Different HIV-1 Resistance Interpretation Tools for Next-Generation Sequencing in Italy. Viruses, 16(9), 1422. https://doi.org/10.3390/v16091422

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop