Nationwide quality assurance of high-throughput diagnostic molecular testing during the SARS-CoV-2 pandemic: role of the Belgian National Reference Centre.

Since the onset of the coronavirus disease (COVID-19) pandemic in Belgium, UZ/KU Leuven has played a crucial role as the National Reference Centre (NRC) for respiratory pathogens, to be the first Belgian laboratory to develop and implement laboratory developed diagnostic assays for SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) and later to assess the quality of commercial kits. To meet the growing demand for decentralised testing, both clinical laboratories and government-supported high-throughput platforms were gradually deployed across Belgium. Consequently, the role of the NRC transitioned from a specialised testing laboratory to strengthening capacity and coordinating quality assurance. Here, we outline the measures taken by the NRC, the national public health institute Sciensano and the executing clinical laboratories to ensure effective quality management of molecular testing throughout the initial two years of the pandemic (March 2020 to March 2022).

Chronic and Early Antiretroviral Therapy Impact Human Immunodeficiency Virus (HIV) Serological Assay Sensitivity, Leading to More False-Negative Test Results in HIV Diagnosis.

This retrospective study evaluated the reactivity of 3 human immunodeficiency virus (HIV) confirmatory assays (INNO-LIA, Geenius, and MP) and 7 HIV rapid tests on samples from 2 different study populations in Belgium. For the early-treated cohort (83 HIV-1 adult patients treated within 3 months after infection), HIV-1 diagnosis was not obtained in at least 1 confirmatory assay in 12.0% (10/83) and in an HIV rapid test in 31.3% (26/83). Confirmation assay sensitivities ranged from 87.5% to 95.2%, whereas rapid test assay sensitivities ranged from 75.9% to 100%. The time to treatment initiation or the length of time on treatment did not have a statistical influence on the probability to obtain a false-negative test result. The fastest reversion was demonstrated after 4 months of treatment. Among the long-term treated cohort (390 HIV-1 patients with ≥ 9 years of undetectable viral load), false-negative test results were found in at least 1 HIV confirmatory assay for 2.1% (8/390) of the patients and in a HIV rapid test for 4.9% (19/390). Confirmation assay sensitivities ranged from 98.1% to 99.5%, whereas rapid test sensitivities ranged from 96.2% to 100%. Longer treatment increased nonreactivity of the HIV rapid tests (P = .033). Undetectable viral load decreases the sensitivities of HIV diagnostic tests, and further monitoring of the performance of serological assays is advised.

An Evolutionary Model-Based Approach To Quantify the Genetic Barrier to Drug Resistance in Fast-Evolving Viruses and Its Application to HIV-1 Subtypes and Integrase Inhibitors.

Viral pathogens causing global disease burdens are often characterized by high rates of evolutionary changes. The extensive viral diversity at baseline can shorten the time to escape from therapeutic or immune selective pressure and alter mutational pathways. The impact of genotypic background on the barrier to resistance can be difficult to capture, particularly for agents in experimental stages or that are recently approved or expanded into new patient populations. We developed an evolutionary model-based counting method to quickly quantify the population genetic potential to resistance and assess population differences. We demonstrate its applicability to HIV-1 integrase inhibitors, as their increasing use globally contrasts with limited availability of non-B subtype resistant sequence data and corresponding knowledge gap. A large sequence data set encompassing most prevailing HIV-1 subtypes and resistance-associated mutations of currently approved integrase inhibitors was investigated. A complex interplay between codon predominance, polymorphisms, and associated evolutionary costs resulted in a subtype-dependent varied genetic potential for 15 resistance mutations against integrase inhibitors. While we confirm the lower genetic barrier of subtype B for G140S, we convincingly discard a similar effect previously suggested for G140C. A supplementary analysis for HIV-1 reverse transcriptase inhibitors identified a lower genetic barrier for K65R in subtype C through differential codon usage not reported before. To aid evolutionary interpretations of genomic differences for antiviral strategies, we advanced existing counting methods with increased sensitivity to identify subtype dependencies of resistance emergence. Future applications include novel HIV-1 drug classes or vaccines, as well as other viral pathogens.

Implications of hepatitis C virus subtype 1a migration patterns for virus genetic sequencing policies in Italy.

BACKGROUND: In-depth phylogeographic analysis can reveal migration patterns relevant for public health planning. Here, as a model, we focused on the provenance, in the current Italian HCV subtype 1a epidemic, of the NS3 resistance-associated variant (RAV) Q80K, known to interfere with the action of NS3/4A protease inhibitor simeprevir. HCV1a migration patterns were analysed using Bayesian phylodynamic tools, capitalising on newly generated and publicly available time and geo-referenced NS3 encoding virus genetic sequence data. RESULTS: Our results showed that both immigration and local circulation fuel the current Italian HCV1a epidemic. The United States and European continental lineages dominate import into Italy, with the latter taking the lead from the 1970s onwards. Since similar migration patterns were found for Q80K and other lineages, no clear differentiation of the risk for failing simeprevir can be made between patients based on their migration and travel history. Importantly, since HCV only occasionally recombines, these results are readily transferable to the genetic sequencing policy concerning NS5A RAVs. CONCLUSIONS: The patient migration and travel history cannot be used to target only part of the HCV1a infected population for drug resistance testing before start of antiviral therapy. Consequently, it may be cost-effective to expand genotyping efforts to all HCV1a infected patients eligible for simeprevir-based therapies.

Phylogenetic evidence for underreporting of male-to-male sex among human immunodeficiency virus-infected donors in the Netherlands and Flanders.

BACKGROUND: Separate transmission networks for human immunodeficiency virus (HIV) coexist. Molecular typing of viral genomes can provide insight in HIV transmission routes in donors for whom risk behavior-based donor selection failed. STUDY DESIGN AND METHODS: This study includes all HIV-infected Dutch and Flemish donors in the period 2005 to 2014 (n = 55). Part of the HIV polymerase (pol) gene was amplified, sequenced, and compared with more than 10,000 HIV strains obtained from HIV-infected Dutch and Flemish patients. The most likely transmission route was determined based on HIV phylogeny and the donor's self-reported risk behavior during the exit interview. RESULTS: HIV-infected donors were predominantly male (69%), were repeat donors (73%), were born in the Netherlands or Belgium (95%), and harbored HIV Subtype B (68%). Seventy-five percent of HIV-infected male donors were part of robust phylogenetic clusters linked to male-to-male sex, while only 24% of HIV-infected male donors reported male-to-male sex during posttest counseling. Sex between men and women accounted for 13% of HIV infections in male donors and 93% of HIV infections in female donors based on phylogenetic analysis. Only 40% of HIV-infected female donors had HIV Subtype B; 65% of female donors reported a foreign partner and indeed HIV sequences interspersed with sequences from HIV-endemic areas abroad, in particular sub-Saharan Africa. CONCLUSION: HIV typing helps to understand HIV transmission routes in donor populations. We found substantial underreporting of male-to-male sex among HIV-infected male donors. Donor education on HIV risk factors and the danger of window-period donations and a donor environment that encourages frank disclosure of sexual behavior will contribute to a decrease of HIV-infected donors.

Impact of HCV genotype on treatment regimens and drug resistance: a snapshot in time.

The introduction of highly potent direct-acting antivirals (DAAs) has revolutionized hepatitis C virus treatment. Nevertheless, viral eradication worldwide remains a challenge also in the era of DAA treatment, because of the high associated costs, high numbers of undiagnosed patients, high re-infection rates in some risk groups and suboptimal drug efficacies associated with host and viral factors as well as advanced stages of liver disease. A correct determination of the HCV genotype allows administration of the most appropriate antiviral regimen. Additionally, HCV genetic sequencing improves our understanding of resistance-associated variants, either naturally occurring before treatment, acquired by transmission at HCV infection, or emerging after virological failure. Because treatment response rates, and the prevalence and development of drug resistance variants differ for each DAA regimen and HCV genotype, this review summarizes treatment opportunities per HCV genotype, and focuses on viral genetic sequencing to guide clinical decision making. Although approval of the first pan-genotypic DAA-only regimen is expected soon, HCV genetic sequencing will remain important because when DAA therapies fail, genotyping and resistance testing to select a new active DAA combination will be essential. Copyright © 2016 John Wiley & Sons, Ltd.

Decision tree for accurate infection timing in individuals newly diagnosed with HIV-1 infection.

BACKGROUND: There is today no gold standard method to accurately define the time passed since infection at HIV diagnosis. Infection timing and incidence measurement is however essential to better monitor the dynamics of local epidemics and the effect of prevention initiatives. METHODS: Three methods for infection timing were evaluated using 237 serial samples from documented seroconversions and 566 cross sectional samples from newly diagnosed patients: identification of antibodies against the HIV p31 protein in INNO-LIA, SediaTM BED CEIA and SediaTM LAg-Avidity EIA. A multi-assay decision tree for infection timing was developed. RESULTS: Clear differences in recency window between BED CEIA, LAg-Avidity EIA and p31 antibody presence were observed with a switch from recent to long term infection a median of 169.5, 108.0 and 64.5 days after collection of the pre-seroconversion sample respectively. BED showed high reliability for identification of long term infections while LAg-Avidity is highly accurate for identification of recent infections. Using BED as initial assay to identify the long term infections and LAg-Avidity as a confirmatory assay for those classified as recent infection by BED, explores the strengths of both while reduces the workload. The short recency window of p31 antibodies allows to discriminate very early from early infections based on this marker. BED recent infection results not confirmed by LAg-Avidity are considered to reflect a period more distant from the infection time. False recency predictions in this group can be minimized by elimination of patients with a CD4 count of less than 100 cells/mm3 or without no p31 antibodies. For 566 cross sectional sample the outcome of the decision tree confirmed the infection timing based on the results of all 3 markers but reduced the overall cost from 13.2 USD to 5.2 USD per sample. CONCLUSIONS: A step-wise multi assay decision tree allows accurate timing of the HIV infection at diagnosis at affordable effort and cost and can be an important new tool in studies analyzing the dynamics of local epidemics or the effects of prevention strategies.

Transmission of HIV Drug Resistance and the Predicted Effect on Current First-line Regimens in Europe.

BACKGROUND: Numerous studies have shown that baseline drug resistance patterns may influence the outcome of antiretroviral therapy. Therefore, guidelines recommend drug resistance testing to guide the choice of initial regimen. In addition to optimizing individual patient management, these baseline resistance data enable transmitted drug resistance (TDR) to be surveyed for public health purposes. The SPREAD program systematically collects data to gain insight into TDR occurring in Europe since 2001. METHODS: Demographic, clinical, and virological data from 4140 antiretroviral-naive human immunodeficiency virus (HIV)-infected individuals from 26 countries who were newly diagnosed between 2008 and 2010 were analyzed. Evidence of TDR was defined using the WHO list for surveillance of drug resistance mutations. Prevalence of TDR was assessed over time by comparing the results to SPREAD data from 2002 to 2007. Baseline susceptibility to antiretroviral drugs was predicted using the Stanford HIVdb program version 7.0. RESULTS: The overall prevalence of TDR did not change significantly over time and was 8.3% (95% confidence interval, 7.2%-9.5%) in 2008-2010. The most frequent indicators of TDR were nucleoside reverse transcriptase inhibitor (NRTI) mutations (4.5%), followed by nonnucleoside reverse transcriptase inhibitor (NNRTI) mutations (2.9%) and protease inhibitor mutations (2.0%). Baseline mutations were most predictive of reduced susceptibility to initial NNRTI-based regimens: 4.5% and 6.5% of patient isolates were predicted to have resistance to regimens containing efavirenz or rilpivirine, respectively, independent of current NRTI backbones. CONCLUSIONS: Although TDR was highest for NRTIs, the impact of baseline drug resistance patterns on susceptibility was largest for NNRTIs. The prevalence of TDR assessed by epidemiological surveys does not clearly indicate to what degree susceptibility to different drug classes is affected.

Possible footprints of APOBEC3F and/or other APOBEC3 deaminases, but not APOBEC3G, on HIV-1 from patients with acute/early and chronic infections.

UNLABELLED: Members of the apolipoprotein B mRNA-editing enzyme-catalytic polypeptide-like-3 (APOBEC3) innate cellular cytidine deaminase family, particularly APOBEC3F and APOBEC3G, can cause extensive and lethal G-to-A mutations in HIV-1 plus-strand DNA (termed hypermutation). It is unclear if APOBEC3-induced mutations in vivo are always lethal or can occur at sublethal levels that increase HIV-1 diversification and viral adaptation to the host. The viral accessory protein Vif counteracts APOBEC3 activity by binding to APOBEC3 and promoting proteasome degradation; however, the efficiency of this interaction varies, since a range of hypermutation frequencies are observed in HIV-1 patient DNA. Therefore, we examined "footprints" of APOBEC3G and APOBEC3F activity in longitudinal HIV-1 RNA pol sequences from approximately 3,000 chronically infected patients by determining whether G-to-A mutations occurred in motifs that were favored or disfavored by these deaminases. G-to-A mutations were more frequent in APOBEC3G-disfavored than in APOBEC3G-favored contexts. In contrast, mutations in APOBEC3F-disfavored contexts were relatively rare, whereas mutations in contexts favoring APOBEC3F (and possibly other deaminases) occurred 16% more often than average G-to-A mutations. These results were supported by analyses of >500 HIV-1 env sequences from acute/early infection. IMPORTANCE: Collectively, our results suggest that APOBEC3G-induced mutagenesis is lethal to HIV-1, whereas mutagenesis caused by APOBEC3F and/or other deaminases may result in sublethal mutations that might facilitate viral diversification. Therefore, Vif-specific cytotoxic T lymphocyte (CTL) responses and drugs that manipulate the interplay between Vif and APOBEC3 may have beneficial or detrimental clinical effects depending on how they affect the binding of Vif to various members of the APOBEC3 family.

The genealogical population dynamics of HIV-1 in a large transmission chain: bridging within and among host evolutionary rates.

Transmission lies at the interface of human immunodeficiency virus type 1 (HIV-1) evolution within and among hosts and separates distinct selective pressures that impose differences in both the mode of diversification and the tempo of evolution. In the absence of comprehensive direct comparative analyses of the evolutionary processes at different biological scales, our understanding of how fast within-host HIV-1 evolutionary rates translate to lower rates at the between host level remains incomplete. Here, we address this by analyzing pol and env data from a large HIV-1 subtype C transmission chain for which both the timing and the direction is known for most transmission events. To this purpose, we develop a new transmission model in a Bayesian genealogical inference framework and demonstrate how to constrain the viral evolutionary history to be compatible with the transmission history while simultaneously inferring the within-host evolutionary and population dynamics. We show that accommodating a transmission bottleneck affords the best fit our data, but the sparse within-host HIV-1 sampling prevents accurate quantification of the concomitant loss in genetic diversity. We draw inference under the transmission model to estimate HIV-1 evolutionary rates among epidemiologically-related patients and demonstrate that they lie in between fast intra-host rates and lower rates among epidemiologically unrelated individuals infected with HIV subtype C. Using a new molecular clock approach, we quantify and find support for a lower evolutionary rate along branches that accommodate a transmission event or branches that represent the entire backbone of transmitted lineages in our transmission history. Finally, we recover the rate differences at the different biological scales for both synonymous and non-synonymous substitution rates, which is only compatible with the 'store and retrieve' hypothesis positing that viruses stored early in latently infected cells preferentially transmit or establish new infections upon reactivation.

HIV-1 Gag C-terminal amino acid substitutions emerging under selective pressure of protease inhibitors in patient populations infected with different HIV-1 subtypes.

HIV-1 Gag amino acid substitutions associated with protease inhibitor (PI) treatment have mainly been reported in subtype B, while information on other subtypes is scarce. Using sequences from 11613 patients infected with different HIV-1 subtypes, we evaluated the prevalence of 93 Gag amino acid substitutions and their association with genotypic PI resistance. A significant association was found for 13 Gag substitutions, including A431V in both subtype B and CRF01_AE. K415R in subtype C and S451G in subtype B were newly identified. Most PI-associated Gag substitutions are located in the flexible C-terminal domain, revealing the key role this region plays in PI resistance.

Comparative performances of HIV-1 RNA load assays at low viral load levels: results of an international collaboration.

Low-level viremia during antiretroviral therapy and its accurate measurement are increasingly relevant. Here, we present an international collaboration of 4,221 paired blood plasma viral load (pVL) results from four commercial assays, emphasizing the data with low pVL. The assays compared were the Abbott RealTime assay, the Roche Amplicor assay, and the Roche TaqMan version 1 and version 2 assays. The correlation between the assays was 0.90 to 0.97. However, at a low pVL, the correlation fell to 0.45 to 0.85. The observed interassay concordance was higher when detectability was defined as 200 copies/ml than when it was defined as 50 copies/ml. A pVL of ∼100 to 125 copies/ml by the TaqMan version 1 and version 2 assays corresponded best to a 50-copies/ml threshold with the Amplicor assay. Correlation and concordance between the viral load assays were lower at a low pVL. Clear guidelines are needed on the clinical significance of low-level viremia.

RegaDB: community-driven data management and analysis for infectious diseases.

SUMMARY: RegaDB is a free and open source data management and analysis environment for infectious diseases. RegaDB allows clinicians to store, manage and analyse patient data, including viral genetic sequences. Moreover, RegaDB provides researchers with a mechanism to collect data in a uniform format and offers them a canvas to make newly developed bioinformatics tools available to clinicians and virologists through a user friendly interface. AVAILABILITY AND IMPLEMENTATION: Source code, binaries and documentation are available on http://rega.kuleuven.be/cev/regadb. RegaDB is written in the Java programming language, using a web-service-oriented architecture.

Declining prevalence of HIV-1 drug resistance in antiretroviral treatment-exposed individuals in Western Europe.

HIV-1 drug resistance represents a major obstacle to infection and disease control. This retrospective study analyzes trends and determinants of resistance in antiretroviral treatment (ART)-exposed individuals across 7 countries in Europe. Of 20 323 cases, 80% carried at least one resistance mutation: these declined from 81% in 1997 to 71% in 2008. Predicted extensive 3-class resistance was rare (3.2% considering the cumulative genotype) and peaked at 4.5% in 2005, decreasing thereafter. The proportion of cases exhausting available drug options dropped from 32% in 2000 to 1% in 2008. Reduced risk of resistance over calendar years was confirmed by multivariable analysis.

Mutations selected in HIV-2-infected patients failing a regimen including atazanavir.

OBJECTIVES: To investigate mutations selected in viruses from HIV-2-infected patients failing a highly active antiretroviral treatment (HAART) regimen including atazanavir/ritonavir. METHODS: Twenty-eight HIV-2-infected patients previously exposed to atazanavir/ritonavir and failing therapy were studied. The protease (PR) gene was amplified and sequenced, and mutations emerging under atazanavir/ritonavir selective pressure were reported. RESULTS: The I50L mutation emerged in 4 out of 28 HIV-2-infected patients failing a HAART regimen including atazanavir/ritonavir. Besides I50L, four PR mutations previously associated with protease inhibitor resistance (I54L, I64V, V71I and I82F) and six PR mutations of unknown impact (V10I, E37D, S43T, K45R, I75V and F85L) in HIV-2 were also identified in this small group of patients. CONCLUSIONS: Several mutations were associated with virological failure of a regimen including atazanavir/ritonavir in HIV-2-infected patients, including I50L for the first time. It should be included in HIV-2 algorithms for interpretation of genotypic resistance data, and taken into account when making therapeutic decisions for HIV-2-infected patients.

Superinfection with drug-resistant HIV is rare and does not contribute substantially to therapy failure in a large European cohort.

BACKGROUND: Superinfection with drug resistant HIV strains could potentially contribute to compromised therapy in patients initially infected with drug-sensitive virus and receiving antiretroviral therapy. To investigate the importance of this potential route to drug resistance, we developed a bioinformatics pipeline to detect superinfection from routinely collected genotyping data, and assessed whether superinfection contributed to increased drug resistance in a large European cohort of viremic, drug treated patients. METHODS: We used sequence data from routine genotypic tests spanning the protease and partial reverse transcriptase regions in the Virolab and EuResist databases that collated data from five European countries. Superinfection was indicated when sequences of a patient failed to cluster together in phylogenetic trees constructed with selected sets of control sequences. A subset of the indicated cases was validated by re-sequencing pol and env regions from the original samples. RESULTS: 4425 patients had at least two sequences in the database, with a total of 13816 distinct sequence entries (of which 86% belonged to subtype B). We identified 107 patients with phylogenetic evidence for superinfection. In 14 of these cases, we analyzed newly amplified sequences from the original samples for validation purposes: only 2 cases were verified as superinfections in the repeated analyses, the other 12 cases turned out to involve sample or sequence misidentification. Resistance to drugs used at the time of strain replacement did not change in these two patients. A third case could not be validated by re-sequencing, but was supported as superinfection by an intermediate sequence with high degenerate base pair count within the time frame of strain switching. Drug resistance increased in this single patient. CONCLUSIONS: Routine genotyping data are informative for the detection of HIV superinfection; however, most cases of non-monophyletic clustering in patient phylogenies arise from sample or sequence mix-up rather than from superinfection, which emphasizes the importance of validation. Non-transient superinfection was rare in our mainly treatment experienced cohort, and we found a single case of possible transmitted drug resistance by this route. We therefore conclude that in our large cohort, superinfection with drug resistant HIV did not compromise the efficiency of antiretroviral treatment.

Phylogenetic Analysis of Hepatitis C Virus Infections in a Large Belgian Cohort Using Next-Generation Sequencing of Full-Length Genomes.

The hepatitis C virus (HCV) epidemic in Western countries is primarily perpetuated by the sub-populations of men who have sex with men (MSM) and people who inject drugs (PWID). Understanding the dynamics of transmission in these communities is crucial for removing the remaining hurdles towards HCV elimination. We sequenced 269 annotated HCV plasma samples using probe enrichment and next-generation sequencing, obtaining 224 open reading frames of HCV (OR497849-OR498072). Maximum likelihood phylogenies were generated on the four most prevalent subtypes in this study (HCV1a, 1b, 3a, 4d) with a subsequent transmission cluster analysis. The highest rate of clustering was observed for HCV4d samples (13/17 (76.47%)). The second highest rate of clustering was observed in HCV1a samples (42/78 (53.85%)) with significant association with HIV-positive MSM. HCV1b and HCV3a had very low rates of clustering (2/83 (2.41%) and (0/29)). The spread of the prevalent subtype HCV1b appears to have been largely curtailed, and we demonstrate the onwards transmission of HCV1a and HCV4d in the HIV-positive MSM population across municipal borders. More systematic data collection and sequencing is needed to allow a better understanding of the HCV transmission among the community of PWID and overcome the remaining barriers for HCV elimination in Belgium.

Treatment-associated polymorphisms in protease are significantly associated with higher viral load and lower CD4 count in newly diagnosed drug-naive HIV-1 infected patients.

BACKGROUND: The effect of drug resistance transmission on disease progression in the newly infected patient is not well understood. Major drug resistance mutations severely impair viral fitness in a drug free environment, and therefore are expected to revert quickly. Compensatory mutations, often already polymorphic in wild-type viruses, do not tend to revert after transmission. While compensatory mutations increase fitness during treatment, their presence may also modulate viral fitness and virulence in absence of therapy and major resistance mutations. We previously designed a modeling technique that quantifies genotypic footprints of in vivo treatment selective pressure, including both drug resistance mutations and polymorphic compensatory mutations, through the quantitative description of a fitness landscape from virus genetic sequences. RESULTS: Genotypic correlates of viral load and CD4 cell count were evaluated in subtype B sequences from recently diagnosed treatment-naive patients enrolled in the SPREAD programme. The association of surveillance drug resistance mutations, reported compensatory mutations and fitness estimated from drug selective pressure fitness landscapes with baseline viral load and CD4 cell count was evaluated using regression techniques. Protease genotypic variability estimated to increase fitness during treatment was associated with higher viral load and lower CD4 cell counts also in treatment-naive patients, which could primarily be attributed to well-known compensatory mutations at highly polymorphic positions. By contrast, treatment-related mutations in reverse transcriptase could not explain viral load or CD4 cell count variability. CONCLUSIONS: These results suggest that polymorphic compensatory mutations in protease, reported to be selected during treatment, may improve the replicative capacity of HIV-1 even in absence of drug selective pressure or major resistance mutations. The presence of this polymorphic variation may either reflect a history of drug selective pressure, i.e. transmission from a treated patient, or merely be a result of diversity in wild-type virus. Our findings suggest that transmitted drug resistance has the potential to contribute to faster disease progression in the newly infected host and to shape the HIV-1 epidemic at a population level.

HIV-1 protease mutation 82M contributes to phenotypic resistance to protease inhibitors in subtype G.

OBJECTIVES: The purpose of this study was the qualitative and quantitative assessment of the in vitro effect of HIV-1 protease (PR) mutation 82M on replication capacity and susceptibility to the eight clinically available PR inhibitors (PIs). METHODS: The 82M substitution was introduced by site-directed mutagenesis in wild-type subtype B and G strains, as well as reverted back to wild-type in a therapy-failing strain. The recombinant viruses were evaluated for their replication capacity and susceptibility to PIs. RESULTS: The single 82M mutation within a wild-type subtype B or G background did not result in drug resistance. However, the in vitro effect of single PR mutations on PI susceptibility is not always distinguishable from wild-type virus, and particular background mutations and polymorphisms are required to detect significant differences in the drug susceptibility profile. Consequently, reverting the 82M mutation back to wild-type (82I) in a subtype G isolate from a patient that failed therapy with multiple other PR mutations did result in significant increases in susceptibility towards indinavir and lopinavir and minor increases in susceptibility towards amprenavir and atazanavir. The presence of the 82M mutation also slightly decreased viral replication, whether it was in the genetic background of subtype B or subtype G. CONCLUSIONS: Our results suggest that 82M has an impact on PI susceptibility and that this effect is not due to a compensatory effect on the replication capacity. Because 82M is not observed as a polymorphism in any subtype, these observations support the inclusion of 82M in drug resistance interpretation systems and PI mutation lists.

Prevalence and Evolution of Transmitted Human Immunodeficiency Virus Drug Resistance in Belgium Between 2013 and 2019.

Background: To assess the prevalence and evolution of transmitted drug resistance (TDR) in Belgium, a total of 3708 baseline human immunodeficiency virus (HIV)-1 polymerase sequences from patients diagnosed between 2013 and 2019 were analyzed. Methods: Protease and reverse-transcriptase HIV-1 sequences were collected from the 7 national Aids Reference Laboratories. Subtype determination and drug resistance scoring were performed using the Stanford HIV Drug Resistance Database. Trends over time were assessed using linear regression, and the maximum likelihood approach was used for phylogenetic analysis. Results: A total of 17.9% of the patients showed evidence of TDR resulting in at least low-level resistance to 1 drug (Stanford score ≥15). If only the high-level mutations (Stanford score ≥60) were considered, TDR prevalence dropped to 6.3%. The majority of observed resistance mutations impacted the sensitivity for nonnucleoside reverse-transcriptase inhibitors (NNRTIs) (11.4%), followed by nucleoside reverse-transcriptase inhibitors (6.2%) and protease inhibitors (2.4%). Multiclass resistance was observed in 2.4%. Clustered onward transmission was evidenced for 257 of 635 patients (40.5%), spread over 25 phylogenetic clusters. Conclusions: The TDR prevalence remained stable between 2013 and 2019 and is comparable to the prevalence in other Western European countries. The high frequency of NNRTI mutations requires special attention and follow-up. Phylogenetic analysis provided evidence for local clustered onward transmission of some frequently detected mutations.