Reliable genotypic tropism tests for the major HIV-1 subtypes

Net charge and position 22 of the V3 loop are associated with HIV-1 tropism in recently infected female sex workers in Nairobi, Kenya

Human immunodeficiency virus (HIV) infection affects around 37 million people worldwide, and in Kenya, key populations especially female sex workers (FSW), are thought to play a substantial role in the wider, mostly heterosexual HIV-1 transmission structure. Notably, HIV tropism has been found to correlate with HIV-1 transmission and disease progression in HIV-infected patients. In this study, recently infected FSWs from Nairobi, Kenya, were assessed for HIV tropism and the factors related to it. We used a cross-sectional study design to analyze 76 HIV-1 positive plasma samples obtained from FSWs enrolled in sex worker outreach program clinics in Nairobi between November 2020 and April 2021. The effects of clinical, demographic, and viral genetic characteristics were determined using multivariable logistic regression. HIV-1 subtype A1 accounted for 89.5% of all cases, with a prevalence of CXCR4-tropic viruses of 26.3%. WebPSSMR5X4 and Geno2Pheno [G2P:10-15% false positive rate] showed high concordance of 88%. Subjects infected with CXCR4-tropic viruses had statistically significant lower baseline CD4+T-cell counts than those infected with CCR5-tropic viruses (P = .044). Using multivariable logistic regression and adjusting for potential confounders, we found that net charge, the amino acid at position 22 of the V3 loop, and the geographic location of the subject were associated with tropism. A unit increase in V3 loop's net-charge increased the odds of a virus being CXCR4-tropic by 2.4 times (OR = 2.40, 95%CI = 1.35-5.00, P = .007). Second, amino acid threonine at position 22 of V3 loop increased the odds of a strain being X4 by 55.7 times compared to the alanine which occurred in CCR5-tropic strains (OR = 55.7, 95%CI = 4.04-84.1, P < .003). The Kawangware sex worker outreach program clinic was associated with CXCR4-tropic strains (P = .034), but there was there was no evidence of a distinct CXCR4-tropic transmission cluster. In conclusion, this study revealed a high concordance of WebPSSMR5X4 and Geno2Pheno in predicting HIV tropism. The most striking finding was that amino acid position 22 of the V3 loop is linked to tropism in HIV-1 subtype A1. Additional studies with a large dataset are warranted to confirm our findings.

Evolution of Multiple Domains of the HIV-1 Envelope Glycoprotein during Coreceptor Switch with CCR5 Antagonist Therapy

HIV-1 uses CD4 as a receptor and chemokine receptors CCR5 and/or CXCR4 as coreceptors. CCR5 antagonists are a class of antiretrovirals used to inhibit viral entry. Phenotypic prediction algorithms such as Geno2Pheno are used to assess CCR5 antagonist eligibility, for which the V3 region is screened. However, there exist scenarios where the algorithm cannot give an accurate prediction of tropism. The current study examined coreceptor shift of HIV-1 from CCR5-tropic strains to CXCR4-tropic or dual-tropic strains among five subjects in a clinical trial of the CCR5 antagonist vicriviroc. Envelope gene amplicon libraries were constructed and subjected to next-generation sequencing, as well as single-clone sequencing and functional analyses. Approximately half of the amplified full-length single envelope-encoding clones had no significant activity for infection of cells expressing high levels of CD4 and CCR5 or CXCR4. Functional analysis of 9 to 21 individual infectious clones at baseline and at the time of VF were used to construct phylogenetic trees and sequence alignments. These studies confirmed that specific residues and the overall charge of the V3 loop were the major determinants of coreceptor use, in addition to specific residues in other domains of the envelope protein in V1/V2, V4, C3, and C4 domains that may be important for coreceptor shift. These results provide greater insight into the viral genetic determinants of coreceptor shift. IMPORTANCE This study is novel in combining single-genome sequence analysis and next-generation sequencing to characterize HIV-1 quasispecies. The work highlights the importance of mutants present at frequencies of 1% or less in development of drug resistance. This study highlights a critical role of specific amino acid substitutions outside V3 that contribute to coreceptor shift as well as important roles of the V1/V2, V4, C3, and C4 domain residues.

Genetic source completeness of HIV-1 circulating recombinant forms (CRFs) predicted by multi-label learning

MOTIVATION

Infection with strains of different subtypes and the subsequent crossover reading between the two strands of genomic RNAs by host cells' reverse transcriptase are the main causes of the vast HIV-1 sequence diversity. Such inter-subtype genomic recombinants can become circulating recombinant forms (CRFs) after widespread transmissions in a population. Complete prediction of all the subtype sources of a CRF strain is a complicated machine learning problem. It is also difficult to understand whether a strain is an emerging new subtype and if so, how to accurately identify the new components of the genetic source.

RESULTS

We introduce a multi-label learning algorithm for the complete prediction of multiple sources of a CRF sequence as well as the prediction of its chronological number. The prediction is strengthened by a voting of various multi-label learning methods to avoid biased decisions. In our steps, frequency and position features of the sequences are both extracted to capture signature patterns of pure subtypes and CRFs. The method was applied to 7185 HIV-1 sequences, comprising 5530 pure subtype sequences and 1655 CRF sequences. Results have demonstrated that the method can achieve very high accuracy (reaching 99%) in the prediction of the complete set of labels of HIV-1 recombinant forms. A few wrong predictions are actually incomplete predictions, very close to the complete set of genuine labels.

AVAILABILITY AND IMPLEMENTATION

https://github.com/Runbin-tang/The-source-of-HIV-CRFs-prediction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Epidemiological data analysis of viral quasispecies in the next-generation sequencing era

The unprecedented coverage offered by next-generation sequencing (NGS) technology has facilitated the assessment of the population complexity of intra-host RNA viral populations at an unprecedented level of detail. Consequently, analysis of NGS datasets could be used to extract and infer crucial epidemiological and biomedical information on the levels of both infected individuals and susceptible populations, thus enabling the development of more effective prevention strategies and antiviral therapeutics. Such information includes drug resistance, infection stage, transmission clusters and structures of transmission networks. However, NGS data require sophisticated analysis dealing with millions of error-prone short reads per patient. Prior to the NGS era, epidemiological and phylogenetic analyses were geared toward Sanger sequencing technology; now, they must be redesigned to handle the large-scale NGS datasets and properly model the evolution of heterogeneous rapidly mutating viral populations. Additionally, dedicated epidemiological surveillance systems require big data analytics to handle millions of reads obtained from thousands of patients for rapid outbreak investigation and management. We survey bioinformatics tools analyzing NGS data for (i) characterization of intra-host viral population complexity including single nucleotide variant and haplotype calling; (ii) downstream epidemiological analysis and inference of drug-resistant mutations, age of infection and linkage between patients; and (iii) data collection and analytics in surveillance systems for fast response and control of outbreaks.

THETA: a new genotypic approach for predicting HIV-1 CRF02-AG coreceptor usage

MOTIVATION

The circulating recombinant form of HIV-1 CRF02-AG is the most frequent non-B subtype in Europe. Anti-HIV therapy and pathophysiological studies on the impact of HIV-1 tropism require genotypic determination of HIV-1 tropism for non-B subtypes. But genotypic approaches based on analysis of the V3 envelope region perform poorly when used to determine the tropism of CRF02-AG. We, therefore, designed an algorithm based on information from the gp120 and gp41 ectodomain that better predicts the tropism of HIV-1 subtype CRF02-AG.

RESULTS

We used a bio-statistical method to identify the genotypic determinants of CRF02-AG coreceptor use. Toulouse HIV Extended Tropism Algorithm (THETA), based on a Least Absolute Shrinkage and Selection Operator method, uses HIV envelope sequence from phenotypically characterized clones. Prediction of R5X4/X4 viruses was 86% sensitive and that of R5 viruses was 89% specific with our model. The overall accuracy of THETA was 88%, making it sufficiently reliable for predicting the tropism of subtype CRF02-AG sequences.

AVAILABILITY AND IMPLEMENTATION

Binaries are freely available for download at https://github.com/viro-tls/THETA. It was implemented in Matlab and supported on MS Windows platform. The sequence data used in this work are available from GenBank under the accession numbers MK618182-MK618417.

Longitudinal analysis of subtype C envelope tropism for memory CD4+ T cell subsets over the first 3 years of untreated HIV-1 infection

Background

HIV-1 infects a wide range of CD4⁺ T cells with different phenotypic properties and differing expression levels of entry coreceptors. We sought to determine the viral tropism of subtype C (C-HIV) Envelope (Env) clones for different CD4⁺ T cell subsets and whether tropism changes during acute to chronic disease progression. HIV-1 envs were amplified from the plasma of five C-HIV infected women from three untreated time points; less than 2 months, 1-year and 3-years post-infection. Pseudoviruses were generated from Env clones, phenotyped for coreceptor usage and CD4⁺ T cell subset tropism was measured by flow cytometry.

Results

A total of 50 C-HIV envs were cloned and screened for functionality in pseudovirus infection assays. Phylogenetic and variable region characteristic analysis demonstrated evolution in envs between time points. We found 45 pseudoviruses were functional and all used CCR5 to mediate entry into NP2/CD4/CCR5 cells. In vitro infection assays showed transitional memory (TM) and effector memory (EM) CD4⁺ T cells were more frequently infected (median: 46% and 25% of total infected CD4⁺ T cells respectively) than naïve, stem cell memory, central memory and terminally differentiated cells. This was not due to these subsets contributing a higher proportion of the CD4⁺ T cell pool, rather these subsets were more susceptible to infection (median: 5.38% EM and 2.15% TM cells infected), consistent with heightened CCR5 expression on EM and TM cells. No inter- or intra-participant changes in CD4⁺ T cell subset tropism were observed across the three-time points.

Conclusions

CD4⁺ T cell subsets that express more CCR5 were more susceptible to infection with C-HIV Envs, suggesting that these may be the major cellular targets during the first 3 years of infection. Moreover, we found that viral tropism for different CD4⁺ T cell subsets in vitro did not change between Envs cloned from acute to chronic disease stages. Finally, central memory, naïve and stem cell memory CD4⁺ T cell subsets were susceptible to infection, albeit inefficiently by Envs from all time-points, suggesting that direct infection of these cells may help establish the latent reservoir early in infection.

Identification of novel molecular determinants of co-receptor usage in HIV-1 subtype F V3 envelope sequences

HIV-1 determinants of coreceptor usage within the gp120 V3 loop have been broadly studied over the past years. This information has led to the development of state-of the-art bioinformatic tools that are useful to predict co-receptor usage based on the V3 loop sequence mainly of subtypes B, C and A. However, these methods show a poor performance for subtype F V3 loops, which are found in an increasing number of HIV-1 strains worldwide. In the present work we investigated determinants of viral tropisms in the understudied subtype F by looking at genotypic and structural information of coreceptor:V3 loop interactions in a novel group of 40 subtype F V3 loops obtained from HIV-1 strains phenotypically characterized either as syncytium inducing or non-syncytium inducing by the MT-2 assay. We provide novel information about estimated interactions energies between a set of V3 loops with known tropism in subtype F, that allowed us to improve predictions of the coreceptor usage for this subtype. Understanding genetic and structural features underlying HIV coreceptor usage across different subtypes is relevant for the rational design of preventive and therapeutic strategies aimed at limiting the HIV-1 epidemic worldwide.

Comparative analyses of error handling strategies for next-generation sequencing in precision medicine

Next-generation sequencing (NGS) offers the opportunity to sequence millions and billions of DNA sequences in a short period, leading to novel applications in personalized medicine, such as cancer diagnostics or antiviral therapy. Nevertheless, sequencing technologies have different error rates, which occur during the sequencing process. If the NGS data is used for diagnostics, these sequences with errors are typically neglected or a worst-case scenario is assumed. In the current study, we focused on the impact of ambiguous bases on therapy recommendations for Human Immunodeficiency Virus 1 (HIV-1) patients. Concretely, we analyzed the treatment recommendation with entry blockers based on prediction models for co-receptor tropism. We compared three different error handling strategies that have been used in the literature, namely (i) neglection, (ii) worst-case assumption, and (iii) deconvolution with a majority vote. We could show that for two or more ambiguous positions per sequence a reliable prediction is generally no longer possible. Moreover, also the position of ambiguity plays a crucial role. Thus, we analyzed the error probability distributions of existing sequencing technologies, e.g., Illumina MiSeq or PacBio, with respect to the aforementioned error handling strategies and it turned out that neglection outperforms the other strategies in the case where no systematic errors are present. In other cases, the deconvolution strategy with the majority vote should be preferred.

Understanding the mechanisms driving the spread of subtype C HIV-1

Human immunodeficiency virus type 1 (HIV-1) subtype C (C-HIV) is the most prevalent form of HIV-1 globally, accounting for approximately 50% of infections worldwide. C-HIV is the predominant and near-exclusive subtype in the low resource regions of India and Southern Africa. Given the vast diversity of HIV-1 subtypes, it is curious as to why C-HIV constitutes such a large proportion of global infections. This enriched prevalence may be due to phenotypic differences between C-HIV isolates and other viral strains that permit enhanced transmission efficiency or, pathogenicity, or might due to the socio-demographics of the regions where C-HIV is endemic. Here, we compare the mechanisms of C-HIV pathogenesis to less prominent HIV-1 subtypes, including viral genetic and phenotypic characteristics, and host genetic variability, to understand whether evolutionary factors drove C-HIV to predominance.

SCOTCH: subtype A coreceptor tropism classification in HIV-1

Motivation

The V3 loop of the gp120 glycoprotein of the Human Immunodeficiency Virus 1 (HIV-1) is considered to be responsible for viral coreceptor tropism. gp120 interacts with the CD4 receptor of the host cell and subsequently V3 binds either CCR5 or CXCR4. Due to the fact that the CCR5 coreceptor is targeted by entry inhibitors, a reliable prediction of the coreceptor usage of HIV-1 is of great interest for antiretroviral therapy. Although several methods for the prediction of coreceptor tropism are available, almost all of them have been developed based on only subtype B sequences, and it has been shown in several studies that the prediction of non-B sequences, in particular subtype A sequences, are less reliable. Thus, the aim of the current study was to develop a reliable prediction model for subtype A viruses.

Results

Our new model SCOTCH is based on a stacking approach of classifier ensembles and shows a significantly better performance for subtype A sequences compared to other available models. In particular for low false positive rates (between 0.05 and 0.2, i.e. recommendation in the German and European Guidelines for tropism prediction), SCOTCH shows significantly better prediction performances in terms of partial area under the curves and diagnostic odds ratios compared to existing tools, and thus can be used to reliably predict coreceptor tropism for subtype A sequences.

Availability and implementation

SCOTCH can be downloaded/accessed at http://www.heiderlab.de.

Near full genome characterization of HIV-1 unique recombinant forms in Cameroon reveals dominant CRF02_AG and F2 recombination patterns

INTRODUCTION

In Cameroon, a manifold diversity of HIV strains exists with CRF02_AG and unique recombinant forms (URFs) being the predominant strains. In recent years, a steady increase in URFs and clade F2 viruses has been monitored through partial genome sequencing. There is an information gap in the characterization of emerging URFs along the full genome, which is needed to address the challenges URFs pose towards diagnosis, treatment and HIV-1 vaccine design.

METHOD

Eighteen Cameroonian URFs from samples collected between the years 2000 and 2015 were studied using a newly developed near full genome sequencing (NFGS) protocol based on variable nested RT-PCRs with a versatile primer set. Near full genomes were characterized for recombination patterns and sequence signatures with possible impact on antiretroviral treatment or Env-directed immune responses. Third-generation sequencing (3GS) of near full or half genomes (HGs) gave insight into intra-patient URF diversity.

RESULTS

The characterized URFs were composed of a broad variety of subtypes and recombinants including A, F, G, CRF01_AE, CRF02_AG and CRF22_01A1. Phylogenetic analysis unveiled dominant CRF02_AG and F2 recombination patterns. 3GS indicated a high intra-patient URF diversity with up to four distinct viral sub-populations present in plasma at the same time. URF pol genomic analysis revealed a number of accessory drug resistance mutations (DRMs) in the ART-naïve participants. Genotypic env analysis suggests CCR5 usage in 14/18 samples and identified deviations at residues, critical for gp120/gp41 interphase and CD4 binding site broadly neutralizing antibodies in more than half of the studied URFs. V1V2 sites of immune pressure in the human RV144 vaccine study varied in more than a third of URFs.

CONCLUSIONS

This study identified novel mosaic patterns in URFs in Cameroon. In line with the regional predominance of CRF_02AG and the increased prevalence of clade F2, prominent CRF_02AG and F2 background patterns were observed underlying the URFs. In the context of the novel mosaic genomes, the impact of the identified accessory DRMs and Env epitope variations on treatment and immune control remains elusive. The evolving diversity of HIV-1 URFs in Cameroon requires continuous monitoring to respond to the increasing challenges for diagnosis, antiretroviral treatment and prevention.

HIVCoR: A sequence-based tool for predicting HIV-1 CRF01_AE coreceptor usage

Determination of HIV-1 coreceptor usage is strongly recommended before starting the coreceptor-specific inhibitors for HIV treatment. Currently, the genotypic assays are the most interesting tools due to they are more feasible than phenotypic assays. However, most of prediction models were developed and validated by data set of HIV-1 subtype B and C. The present study aims to develop a powerful and reliable model to accurately predict HIV-1 coreceptor usage for CRF01_AE subtype called HIVCoR. HIVCoR utilized random forest and support vector machine as the prediction model, together with amino acid compositions, pseudo amino acid compositions and relative synonymous codon usage frequencies as the input feature. The overall success rate of 93.79% was achieved from the external validation test on the objective benchmark dataset. Comparison results indicated that HIVCoR was superior to other bioinformatics tools and genotypic predictors. For the convenience of experimental scientists, a user-friendly webserver has been established at http://codes.bio/hivcor/.

Coreceptor Tropism and Maraviroc Sensitivity of Clonally Derived Ethiopian HIV-1C Strains Using an in-house Phenotypic Assay and Commonly Used Genotypic Methods

OBJECTIVES

Genotypic Tropism Testing (GTT) tools are generally developed based on HIV-1 subtype B (HIV-1B) and used for HIV-1C as well but with a large discordance of prediction between different methods. We used an established phenotypic assay for comparison with GTT methods and for the determination of in vitro maraviroc sensitivity of pure R5-tropic and dual-tropic HIV-1C.

METHODS

Plasma was obtained from 58 HIV-1C infected Ethiopians. Envgp120 was cloned into a luciferase tagged NL4-3 plasmid. Phenotypic tropism was determined by in house method and the V3 sequences were analysed by five GTT methods. In vitro maraviroc sensitivity of R5-tropic and dual-tropic isolates were compared in the TZMbl cell-line.

RESULTS

The phenotypes were classified as R5 in 92.4% and dual tropic (R5X4) in 7.6% of 79 clones. The concordance between phenotype and genotype ranged from 64.7% to 84.3% depending on the GTT method. Only 46.9% of the R5 phenotypes were predicted as R5 by all GTT tools while R5X4 phenotypes were predicted as X4 by four methods, but not by Raymond's method. All six tested phenotypic R5 clones, as well as five of six of dual tropic clones, showed a dose response to maraviroc.

CONCLUSION

There is a high discordance between GTT methods, which underestimates the presence of R5 and overestimates X4 strains compared to a phenotypic assay. Currently available GTT algorithms should be further improved for tropism prediction in HIV-1C. Maraviroc has an in vitro activity against most HIV-1C viruses and could be considered as an alternative regimen in individuals infected with CCR5-tropic HIV-1C viruses.

Transmission patterns of HIV-1 non-R5 strains in Poland

HIV-1 env sequencing enables predictions of viral coreceptor tropism and phylogenetic investigations of transmission events. The aim of the study was to estimate the contribution of non-R5 strains to the viral spread in Poland. Partial proviral env sequences were retrieved from baseline blood samples of patients with newly diagnosed HIV-1 infection between 2008-2014, including 46 patients with recent HIV-1 infection (RHI), and 246 individuals with long-term infection (LTHI). These sequences were subjected to the genotypic coreceptor tropism predictions and phylogenetic analyses to identify transmission clusters. Overall, 27 clusters with 57 sequences (19.5%) were detected, including 15 sequences (26.3%) from patients with RHI. The proportion of non-R5 strains among all study participants was 23.3% (68/292), and was comparable between patients with RHI and LTHI (11/46, 23.9% vs 57/246, 23.2%; p = 1.000). All 11 patients with non-R5 strains and RHI were men having sex with men (MSM). Among these patients, 4 had viral sequences grouped within phylogenetic cluster with another sequence of non-R5 strain obtained from patient with LTHI, indicating potential acquisition of non-R5 HIV-1 for at least 4/46 (8.7%) patients with RHI. We were unable to confirm the contribution of patients with RHI to the forward transmission of non-R5 strains, but a relatively high proportion of non-R5 strains among them deserves attention due to the limited susceptibility to CCR5 antagonists.

The Molecular Basis of pH-Modulated HIV gp120 Binding Revealed

Decades of research has yet to provide a vaccine for HIV, the virus which causes AIDS. Recent theoretical research has turned attention to mucosa pH levels over systemic pH levels. Previous research in this field developed a computational approach for determining pH sensitivity that indicated higher potential for transmission at mucosa pH levels present during intercourse. The process was extended to incorporate a principal component analysis (PCA)-based machine learning technique for classification of gp120 proteins against a known transmitted variant called Biomolecular Electro-Static Indexing (BESI). The original process has since been extended to the residue level by a process we termed Electrostatic Variance Masking (EVM) and used in conjunction with BESI to determine structural differences present among various subspecies across Clades A1 and C. Results indicate that structures outside of the core selected by EVM may be responsible for binding affinity observed in many other studies and that pH modulation of select substructures indicated by EVM may influence specific regions of the viral envelope protein (Env) involved in protein-protein interactions.

Individuals with HIV-1 Subtype C Infection and Cryptococcal Meningitis Exhibit Viral Genetic Intermixing of HIV-1 Between Plasma and Cerebrospinal Fluid and a High Prevalence of CXCR4-Using Variants

The genotypic properties of human immunodeficiency virus type 1 (HIV-1) subtype C in individuals presenting with cryptococcal meningitis (CM) are not well established. Employing single-genome amplification as well as bulk PCR, cloning and sequencing strategies, we evaluated the genetic properties of HIV-1 subtype C env in 16 antiretroviral therapy-naive study participants with CM. Eleven of the 16 participants had matched blood plasma and cerebrospinal fluid (CSF) evaluated, with the rest having either a plasma or CSF sample evaluated. Before antiretroviral therapy initiation, matched plasma and CSF-derived env sequences of all 11 participants displayed genetic intermixing between the two compartments. Overall, 7 of the 16 (∼43.8%) participants harbored CXCR4-using variants in plasma and/or CSF, according to coreceptor usage prediction algorithms. This study suggests that HIV-1 subtype C genetic intermixing between peripheral blood and the central nervous system is common in individuals presenting with CM, and that CXCR4 usage is present in one or both compartments in approximately 44% of individuals.

Novel Strategy To Adapt Simian-Human Immunodeficiency Virus E1 Carrying env from an RV144 Volunteer to Rhesus Macaques: Coreceptor Switch and Final Recovery of a Pathogenic Virus with Exclusive R5 Tropism

The phase III RV144 human immunodeficiency virus (HIV) vaccine trial conducted in Thailand remains the only study to show efficacy in decreasing the HIV acquisition risk. In Thailand, circulating recombinant forms of HIV clade A/E (CRF01_AE) predominate; in such viruses, env originates from clade E (HIV-E). We constructed a simian-human immunodeficiency virus (SHIV) chimera carrying env isolated from an RV144 placebo recipient in the SHIV-1157ipd3N4 backbone. The latter contains long terminal repeats (LTRs) with duplicated NF-κB sites, thus resembling HIV LTRs. We devised a novel strategy to adapt the parental infectious molecular clone (IMC), R5 SHIV-E1, to rhesus macaques: the simultaneous depletion of B and CD8⁺ cells followed by the intramuscular inoculation of proviral DNA and repeated administrations of cell-free virus. High-level viremia and CD4⁺ T-cell depletion ensued. Passage 3 virus unexpectedly caused acute, irreversible CD4⁺ T-cell loss; the partially adapted SHIV had become dual tropic. Virus and IMCs with exclusive R5 tropism were reisolated from earlier passages, combined, and used to complete adaptation through additional macaques. The final isolate, SHIV-E1p5, remained solely R5 tropic. It had a tier 2 neutralization phenotype, was mucosally transmissible, and was pathogenic. Deep sequencing revealed 99% Env amino acid sequence conservation; X4-only and dual-tropic strains had evolved independently from an early branch of parental SHIV-E1. To conclude, our primate model data reveal that SHIV-E1p5 recapitulates important aspects of HIV transmission and pathobiology in humans.IMPORTANCE Understanding the protective principles that lead to a safe, effective vaccine against HIV in nonhuman primate (NHP) models requires test viruses that allow the evaluation of anti-HIV envelope responses. Reduced HIV acquisition risk in RV144 has been linked to nonneutralizing IgG antibodies with a range of effector activities. Definitive experiments to decipher the mechanisms of the partial protection observed in RV144 require passive-immunization studies in NHPs with a relevant test virus. We have generated such a virus by inserting env from an RV144 placebo recipient into a SHIV backbone with HIV-like LTRs. The final SHIV-E1p5 isolate, grown in rhesus monkey peripheral blood mononuclear cells, was mucosally transmissible and pathogenic. Earlier SHIV-E passages showed a coreceptor switch, again mimicking HIV biology in humans. Thus, our series of SHIV-E strains mirrors HIV transmission and disease progression in humans. SHIV-E1p5 represents a biologically relevant tool to assess prevention strategies.

Next-generation sequencing provides an added value in determining drug resistance and viral tropism in Cameroonian HIV-1 vertically infected children

With limited and low-genetic barrier drugs used for the prevention of mother-to-child transmission (PMTCT) of HIV in sub-Saharan Africa, vertically transmitted HIV-1 drug-resistance (HIVDR) is concerning and might prompt optimal pediatric strategies.The aim of this study was to ascertain HIVDR and viral-tropism in majority and minority populations among Cameroonian vertically infected children.A comparative analysis among 18 HIV-infected children (7 from PMTCT-exposed mothers and 11 from mothers without PMTCT-exposure) was performed. HIVDR and HIV-1 co-receptor usage was evaluated by analyzing sequences obtained by both Sanger sequencing and ultra-deep 454-pyrosequencing (UDPS), set at 1% threshold.Overall, median (interquartile range) age, viremia, and CD4 count were 6 (4-10) years, 5.5 (4.9-6.0) log10 copies/mL, and 526 (282-645) cells/mm, respectively. All children had wild-type viruses through both Sanger sequencing and UDPS, except for 1 PMTCT-exposed infant harboring minority K103N (8.31%), born to a mother exposed to AZT+3TC+NVP. X4-tropic viruses were found in 5 of 15 (33.3%) children (including 2 cases detected only by UDPS). Rate of X4-tropic viruses was 0% (0/6) below 5 years (also as minority species), and became relatively high above 5 years (55.6% [5/9], P = .040. X4-tropic viruses were higher with CD4 ≤15% (4/9 [44.4%]) versus CD4 >15% (1/6 [16.7%], P = .580); similarly for CD4 ≤200 (3/4 [75%]) versus CD4 >200 (2/11 [18.2%] cells/mm, P = .077.NGS has the ability of excluding NRTI- and NNRTI-mutations as minority species in all but 1 children, thus supporting the safe use of these drug-classes in those without such mutations, henceforth sparing ritonavir-boosted protease inhibitors or integrase inhibitors for the few remaining cases. In children under five years, X4-tropic variants would be rare, suggesting vertical-transmission with CCR5-tropic viruses and possible maraviroc usage at younger ages.

Next generation sequencing reveals a high frequency of CXCR4 utilizing viruses in HIV-1 chronically infected drug experienced individuals in South Africa

BACKGROUND

Entry inhibitors, such as Maraviroc, bind to CCR5 inhibiting entry of CCR5 utilizing viruses (R5 viruses). In the course of HIV infection, CXCR4 utilizing viruses (X4 viruses) may emerge and outgrow R5 viruses, and potentially limit the effectiveness of Maraviroc. The use of Maraviroc is reserved for salvage therapy in South Africa.

OBJECTIVE

In this study, we examined the frequency of R5 and X4 viruses, using next generation sequencing, in patients under treatment to draw inferences on the utility of Maraviroc in a South African population.

STUDY DESIGN

Proviral DNA was isolated from peripheral blood mononuclear cells (PBMC) of 72 chronically HIV infected patients on antiretroviral treatment. HIV V3 loop gene was amplified and sequenced on an Illumina MiniSeq platform. Viral subtypes were determined by the jumping profile Hidden Markov Model (jpHMM) and REGA genotyping tools. De Novo consensus sequences were derived for the majority and minority populations for each patient using Geneious^® software version 8.1.5. HIV-1 tropism was inferred using PSSMsinsi, Geno2pheno and Phenoseq-C web-based tools.

RESULTS

Quality V3 loop sequences were obtained from 72 patients, with 5 years (range: 0-16) median duration on treatment. Subtypes A1, B and C viruses were identified at frequencies of 4% (3/72), 4% (3/72) and 92% (66/72) respectively. Fifty four percent (39/72) of patients exclusively harboured R5 viral quasispecies; and 21% (15/72) exclusively harbored X4 viral quasispecies. Twenty five percent of patients (18/72) harbored dual/mixture of R5X4 quasispecies. Of these 18 patients, about 28% (5/18) harbored the R5+X4, a mixture with a majority R5 and minority X4 viruses, while about 72% (13/18) harbored the R5X4+ mixture with a majority X4 and minority R5 viruses. The proportion of all patients who harbored X4 viruses either exclusively or dual/mixture was 46% (33/72). Thirty-five percent (23/66) of the patients who were of HIV-1 subtype C harboured X4 viruses (χ² = 3.58; p = .058), and 57% of these (13/23) harbored X4 viruses exclusively. CD4⁺ cell count less than 350 cell/μl was associated with the presence of X4 viruses (χ² = 4.99; p = .008).

CONCLUSION

The effectiveness of Maraviroc as a component in salvage therapy may be compromised for a significant number of chronically infected patients harboring CXCR4 utilizing viruses.

Phenotypic co-receptor tropism and Maraviroc sensitivity in HIV-1 subtype C from East Africa

Genotypic tropism testing (GTT) for co-receptor usage is a recommended tool for clinical practice before administration of the CCR5-antagonist maraviroc. For some isolates, phenotypic tropism testing (PTT) revealed discordant results with GTT. In this study, we performed a comparative study between GTT and PTT in HIV-1C from East Africa (HIV-1C_EA) and compared the data with HIV-1B and 01_AE and described the maraviroc susceptibility in the CCR5-tropic strains. Patient-derived HIV-1 envgp120 region was cloned into a modified pNL4-3 plasmid expressing the luciferase gene. rPhenotyping dissected single clones from 31 HIV-1C_EA infected patients and four strains with known phenotype. Additionally, 68 clones from 18 patients (HIV-1B: 5, 01_AE: 7, HIV-1C_EA: 6) were used to determine the PTT in GHOST cell line. The respective V3-sequences were used for GTT. R5-tropic strains from HIV-1C_EA (n = 20) and non-C (n = 12) were tested for maraviroc sensitivity in TZMbl cell line. The GTT falsely called a higher proportion of X4-tropic strains in HIV-1C_ET compared to PTT by both rPhenotyping and the GHOST-cell assay. When multiple clones were tested in a subset of patients’ samples, both dual-tropic and R5-tropic strains were identified for HIV-1C. Relatively higher EC₅₀ values were observed in HIV-1C strains than the non-C strains (p = 0.002).

Prediction of coreceptor usage by five bioinformatics tools in a large Ethiopian HIV-1 subtype C cohort

Background

Genotypic tropism testing (GTT) has been developed largely on HIV-1 subtype B. Although a few reports have analysed the utility of GTT in other subtypes, more studies using HIV-1 subtype C (HIV-1C) are needed, considering the huge contribution of HIV-1C to the global epidemic.

Methods

Plasma was obtained from 420 treatment-naïve HIV-1C infected Ethiopians recruited 2009–2011. The V3 region was sequenced and the coreceptor usage was predicted by five tools: Geno2Pheno clinical–and clonal–models, PhenoSeq-C, C-PSSM and Raymond’s algorithm. The impact of baseline tropism on antiretroviral treatment (ART) outcome was evaluated.

Results

Of 352 patients with successful baseline V3 sequences, the proportion of predicted R5 virus varied between the methods by 12.5% (78.1%-90.6%). However, only 58.2% of the predictions were concordant and only 1.7% were predicted to be X4-tropic across the five methods. Compared pairwise, the highest concordance was between C-PSSM and Geno2Pheno clonal (86.4%). In bivariate intention to treat (ITT) analysis, R5 infected patients achieved treatment success more frequently than X4 infected at month six as predicted by Geno2Pheno clinical (77.8% vs 58.7%, P = 0.004) and at month 12 by C-PSSM (61.9% vs 46.6%, P = 0.038). However, in the multivariable analysis adjusted for age, gender, baseline CD4 and viral load, only tropism as predicted by C-PSSM showed an impact on month 12 (P = 0.04, OR 2.47, 95% CI 1.06–5.79).

Conclusion

Each of the bioinformatics models predicted R5 tropism with comparable frequency but there was a large discordance between the methods. Baseline tropism had an impact on outcome of first line ART at month 12 in multivariable ITT analysis but only based on prediction by C-PSSM which thus possibly could be used for predicting outcome of ART in HIV-1C infected Ethiopians.

HIV-1 and SIV Predominantly Use CCR5 Expressed on a Precursor Population to Establish Infection in T Follicular Helper Cells

Background

T follicular helper (Tfh) cells are increasingly recognized as a major reservoir of HIV infection that will likely need to be addressed in approaches to curing HIV. However, Tfh express minimal CCR5, the major coreceptor for HIV-1, and the mechanism by which they are infected is unclear. We have previously shown that macaque Tfh lack CCR5, but are infected in vivo with CCR5-using SIV at levels comparable to other memory CD4⁺ T cells. Similarly, human splenic Tfh cells are highly infected with HIV-1 DNA. Therefore, we set out to examine the mechanism of infection of Tfh cells.

Methodology

Tfh and other CD4⁺ T cell subsets from macaque lymph nodes and spleens, splenic Tfh from HIV⁺ subjects, and tonsillar Tfh from HIV-uninfected subjects were isolated by cell sorting prior to cell surface and molecular characterization. HIV proviral gp120 sequences were submitted to genotypic and phenotypic tropism assays. Entry of CCR5- and CXCR4-using viruses into Tfh from uninfected tonsillar tissue was measured using a fusion assay.

Results

Phylogenetic analysis, genotypic, and phenotypic analysis showed that splenic Tfh cells from chronic HIV⁺ subjects were predominantly infected with CCR5-using viruses. In macaques, purified CCR5⁺PD-1^{intermediate(int)+} memory CD4⁺ T cells were shown to include pre-Tfh cells capable of differentiating in vitro to Tfh by upregulation of PD-1 and Bcl6, confirmed by qRT-PCR and single-cell multiplex PCR. Infected PD-1^int cells survive, carry SIV provirus, and differentiate into PD-1^hi Tfh after T cell receptor stimulation, suggesting a pathway for SIV infection of Tfh. In addition, a small subset of macaque and human PD-1^hi Tfh can express low levels of CCR5, which makes them susceptible to infection. Fusion assays demonstrated CCR5-using HIV-1 entry into CCR5⁺ Tfh and pre-Tfh cells from human tonsils.

Conclusion

The major route of infection of Tfh in macaques and humans appears to be via a CCR5-expressing pre-Tfh population. As the generation of Tfh are important for establishing effective immune responses during primary infections, Tfh are likely to be an early target of HIV-1 following transmission, creating an important component of the reservoir that has the potential to expand over time.

Partial HIV C2V3 envelope sequence analysis reveals association of coreceptor tropism, envelope glycosylation and viral genotypic variability among Kenyan patients on HAART

Background

HIV-1 is highly variable genetically and at protein level, a property it uses to subvert antiviral immunity and treatment. The aim of this study was to assess if HIV subtype differences were associated with variations in glycosylation patterns and co-receptor tropism among HAART patients experiencing different virologic treatment outcomes.

Methods

A total of 118 HIV env C2V3 sequence isolates generated previously from 59 Kenyan patients receiving highly active antiretroviral therapy (HAART) were examined for tropism and glycosylation patterns. For analysis of Potential N-linked glycosylation sites (PNGs), amino acid sequences generated by the NCBI’s Translate tool were applied to the HIVAlign and the N-glycosite tool within the Los Alamos Database. Viral tropism was assessed using Geno2Pheno (G2P), WebPSSM and Phenoseq platforms as well as using Raymond’s and Esbjörnsson’s rules. Chi square test was used to determine independent variables association and ANOVA applied on scale variables.

Results

At respective False Positive Rate (FPR) cut-offs of 5% (p = 0.045), 10% (p = 0.016) and 20% (p = 0.005) for CXCR4 usage within the Geno2Pheno platform, HIV-1 subtype and viral tropism were significantly associated in a chi square test. Raymond’s rule (p = 0.024) and WebPSSM (p = 0.05), but not Phenoseq or Esbjörnsson showed significant associations between subtype and tropism. Relative to other platforms used, Raymond’s and Esbjörnsson’s rules showed higher proportions of X4 variants, while WebPSSM resulted in lower proportions of X4 variants across subtypes. The mean glycosylation density differed significantly between subtypes at positions, N277 (p = 0.034), N296 (p = 0.036), N302 (p = 0.034) and N366 (p = 0.004), with HIV-1D most heavily glycosylated of the subtypes. R5 isolates had fewer PNGs than X4 isolates, but these differences were not significant except at position N262 (p = 0.040). Cell-associated isolates from virologic treatment success subjects were more glycosylated than cell-free isolates from virologic treatment failures both for the NXT (p = 0.016), and for all the patterns (p = 0.011).

Conclusion

These data reveal significant associations of HIV-1 subtype diversity, viral co-receptor tropism, viral suppression and envelope glycosylation. These associations have important implications for designing therapy and vaccines against HIV. Heavy glycosylation and preference for CXCR4 usage of HIV-1D may explain rapid disease progression in patients infected with these strains.

Monophylogenetic HIV-1C epidemic in Ethiopia is dominated by CCR5-tropic viruses-an analysis of a prospective country-wide cohort

Background

CCR5 coreceptor using HIV-1 subtype C (HIV-1C) has been reported to dominate the Ethiopian epidemic. However, almost all data have been obtained from two large cities in the central and north-west regions and recent data is lacking.

Methods

Plasma were obtained from 420 treatment-naïve patients recruited 2009–2011 to a large country-wide Ethiopian cohort. The V3 region was sequenced and the co-receptor tropism was predicted by the clinical and clonal models of the geno2pheno tool at different false positive rates (fpr) and for subtype. In an intention to treat analysis the impact of baseline tropism on outcome of antiretroviral therapy was evaluated.

Results

V3 loop sequencing was successful in 352 (84%) patients. HIV-1C was found in 350 (99.4%) and HIV-1A in two (0.6%) patients. When comparing the geno2pheno fpr10% clonal and clinical models, 24.4% predictions were discordant. X4-virus was predicted in 17.0 and 19.0%, respectively, but the predictions were concordant in only 6%. At fpr5%, concordant X4-virus predictions were obtained in 3.1%. The proportion of X4-tropic virus (clonal fpr10%) increased from 5.6 to 17.3% (p < 0.001) when 387 Ethiopian V3 loop sequences dated from 1984 to 2003 were compared with ours. In an intention to treat analysis, 67.9% reached treatment success at month 6 and only 50% at month 12. Only age and not tropism predicted therapy outcome and no difference was found in CD4+ cell gain between R5-tropic and X4-tropic infected patients. At viral failure, R5 to X4 switch was rare while X4 to R5 switch occurred more frequently (month 6: p = 0.006; month 12: p = 0.078).

Conclusion

The HIV-1C epidemic is monophylogenetic in all regions of Ethiopia and R5-tropic virus dominates, even in patients with advanced immunodeficiency, although the proportion of X4-tropic virus seems to have increased over the last two decades. Geno2pheno clinical and clonal prediction models show a large discrepancy at fpr10%, but not at fpr5%. Hence further studies are needed to assess the utility of genotypic tropism testing in HIV-1C. In ITT analysis only age and not tropism influenced the outcome.

Non-R5-tropic HIV-1 in subtype A1 and D infections were associated with lower pretherapy CD4+ cell count but not with PI/(N)NRTI therapy outcomes in Mbarara, Uganda

BACKGROUND

Previous studies suggest that infection with non-R5-tropic subtype B HIV-1, compared with R5, is associated with a more rapid decline in CD4 cell count, but does not affect PI/(N)NRTI therapy outcome. Here, we explored clinical correlates associated with viral tropism in subtype A1 and D infections.

METHODS

HIV-1 subtype A1 (n = 196) and D (n = 143) pretherapy plasma samples and up to 7.5 years of posttherapy virologic and CD4 data were collected from a cross-sectional cohort in Mbarara, Uganda. Tropism and subtype were inferred using env V3 (geno2pheno) and gp41 (RIP) Sanger sequences. For each subtype, R5 infection was compared with non-R5 in terms of: pretherapy viral load and CD4 cell count (Mann-Whitney tests), and therapy outcomes, including time to virologic suppression, postsuppression virologic rebound, CD4 decline and CD4 recovery (log-rank tests).

RESULTS

A 94% of all patients in this study achieved virologic suppression within median 3 months posttherapy. In both subtypes, non-R5 infection was associated with lower pretherapy CD4 cell count (non-R5 vs. R5; A: median 57 vs. 147 cells/μl P = 0.005; D: 80 vs. 128 cells/μl P = 0.006). Multivariable linear regression confirmed that tropism, not subtype nor the interaction between subtype and tropism, was a significant predictor of pretherapy CD4 cell count (P < 0.0001). None of pretherapy viral load, time to virologic suppression, virologic rebound, CD4 decline nor CD4 recovery was significantly different (all P > 0.09).

CONCLUSION

Regardless of HIV-1 subtype or tropism, the majority of patients in this Ugandan cohort responded to therapy, even though non-R5 infection was associated with lower pretherapy CD4 cell count.

Sequential CCR5-Tropic HIV-1 Reactivation from Distinct Cellular Reservoirs following Perturbation of Elite Control

Background

HIV Elite Controllers may reveal insights into virus persistence given they harbour small reservoir sizes, akin to HIV non-controllers treated early with combination antiretroviral therapy. Both groups of patients represent the most promising candidates for interventions aimed at sustained remission or ‘cure’. Analytic treatment interruption (ATI) in the latter group leads to stochastic rebound of virus, though it is unclear whether loss of elite control is also associated with similar rebound characteristics.

Methods

We studied three discrete periods of virus rebound during myeloma related immune disruption over 2.5 years in an elite controller who previously underwent autologous stem cell transplantation (ASCT) in the absence of any antiretroviral therapy. Single genome sequencing of the V1-V4 region of env in PBMC and plasma was performed and phylogenies reconstructed. Average pairwise distance (APD) was calculated and non-parametric methods used to assess compartmentalisation. Coreceptor usage was predicted based on genotypic algorithms.

Results

122 single genome sequences were obtained (median 26 sequences per rebound). The initial rebounding plasma env sequences following ASCT represented two distinct lineages, and clustered with proviral DNA sequences isolated prior to ASCT. One of the lineages was monophyletic, possibly indicating reactivation from clonally expanded cells. The second rebound occurred 470 days after spontaneous control of the first rebound and was phylogenetically distinct from the first, confirmed by compartmentalisation analysis, with a different cellular origin rather than ongoing replication. By contrast, third rebound viruses clustered with second rebound viruses, with evidence for ongoing evolution that was associated with lymphopenia and myeloma progression. Following ASCT a shift in tropism from CXCR4-tropic viruses to a CCR5-tropic population was observed to persist through to the third rebound.

Conclusions

Our data highlight similarities in the viral reservoir between elite and non-controllers undergoing ATI following allogeneic transplantation. The lack of propagation of CXCR4 using viruses following transplantation warrants further study.

Phylogenetic Analysis of Ethiopian HIV-1 Subtype C Near Full-Length Genomes Reveals High Intrasubtype Diversity and a Strong Geographical Cluster

In this study, we characterize HIV-1 subtype C (HIV-1C) strains at the near full-length genome (NFLG) level and perform genotypic drug resistance testing (GRT) and genotypic tropism testing (GTT) from Ethiopia (HIV-1CET). Plasma samples (n = 150) were obtained from therapy-naive individuals residing in Addis Ababa, Ethiopia in 2008. HIV-NFLG was performed in a subset of patients (n = 30). GRT (pol) and GTT (V3 env) were performed using in-house methods. GTT was analyzed by PhenoSeq-C. The phylogenetic analysis of the NLFG identified two separate clusters of HIV-1CET, although all strains formed one large overarching cluster together. At NFLG, greater diversity was found among HIV-1CET strains compared to HIV-1C strains from other geographical locations. The geographic clustering was weak in the small subgenomic (pol and env) regions. The primary drug-resistant mutations were identified at a low level (<5%). GTT identified that 12% (12/102) of the patients were predicted to be harboring X4-tropic or both R5/X4-tropic viruses.

Genotypic Prediction of Co-receptor Tropism of HIV-1 Subtypes A and C

Antiretroviral treatment of Human Immunodeficiency Virus type-1 (HIV-1) infections with CCR5-antagonists requires the co-receptor usage prediction of viral strains. Currently available tools are mostly designed based on subtype B strains and thus are in general not applicable to non-B subtypes. However, HIV-1 infections caused by subtype B only account for approximately 11% of infections worldwide. We evaluated the performance of several sequence-based algorithms for co-receptor usage prediction employed on subtype A V3 sequences including circulating recombinant forms (CRFs) and subtype C strains. We further analysed sequence profiles of gp120 regions of subtype A, B and C to explore functional relationships to entry phenotypes. Our analyses clearly demonstrate that state-of-the-art algorithms are not useful for predicting co-receptor tropism of subtype A and its CRFs. Sequence profile analysis of gp120 revealed molecular variability in subtype A viruses. Especially, the V2 loop region could be associated with co-receptor tropism, which might indicate a unique pattern that determines co-receptor tropism in subtype A strains compared to subtype B and C strains. Thus, our study demonstrates that there is a need for the development of novel algorithms facilitating tropism prediction of HIV-1 subtype A to improve effective antiretroviral treatment in patients.

Highly Accurate Structure-Based Prediction of HIV-1 Coreceptor Usage Suggests Intermolecular Interactions Driving Tropism

HIV-1 entry into host cells is mediated by interactions between the V3-loop of viral glycoprotein gp120 and chemokine receptor CCR5 or CXCR4, collectively known as HIV-1 coreceptors. Accurate genotypic prediction of coreceptor usage is of significant clinical interest and determination of the factors driving tropism has been the focus of extensive study. We have developed a method based on nonlinear support vector machines to elucidate the interacting residue pairs driving coreceptor usage and provide highly accurate coreceptor usage predictions. Our models utilize centroid-centroid interaction energies from computationally derived structures of the V3-loop:coreceptor complexes as primary features, while additional features based on established rules regarding V3-loop sequences are also investigated. We tested our method on 2455 V3-loop sequences of various lengths and subtypes, and produce a median area under the receiver operator curve of 0.977 based on 500 runs of 10-fold cross validation. Our study is the first to elucidate a small set of specific interacting residue pairs between the V3-loop and coreceptors capable of predicting coreceptor usage with high accuracy across major HIV-1 subtypes. The developed method has been implemented as a web tool named CRUSH, CoReceptor USage prediction for HIV-1, which is available at http://ares.tamu.edu/CRUSH/.

Characteristics and spread to the native population of HIV-1 non-B subtypes in two European countries with high migration rate

Background

Non-B subtypes account for at least 50 % of HIV-1 infections diagnosed in Belgium and Luxembourg. They are considered to be acquired through heterosexual contacts and infect primarily individuals of foreign origin. Information on the extent to which non-B subtypes spread to the local population is incomplete.

Methods

Pol and env gene sequences were collected from 410 non-subtype B infections. Profound subtyping was performed using 5 subtyping tools and sequences of both pol and env. Demographic information, disease markers (viral load, CD4 count) and viral characteristics (co-receptor tropism) were compared between subtypes. Maximum likelihood phylogenetic trees were constructed and examined for clustering.

Results

The majority of non-B infections were diagnosed in patients originating from Africa (55.8 %), individuals born in Western Europe represented 30.5 %. Heterosexual transmission was the most frequently reported transmission route (79.9 %), MSM transmission accounted for 12.2 % and was significantly more frequently reported for Western Europeans (25.7 % versus 4.3 % for individuals originating from other regions; p < 0.001). Subtypes A and C and the circulating recombinant forms CRF01_AE and CRF02_AG were the most represented and were included in the comparative analysis. Native Western Europeans were underrepresented for subtype A (14.5 %) and overrepresented for CRF01_AE (38.6 %). The frequency of MSM transmission was the highest for CRF01_AE (18.2 %) and the lowest for subtype A (0 %). No differences in age, gender, viral load or CD4 count were observed. Prevalence of CXCR4-use differed between subtypes but largely depended on the tropism prediction algorithm applied. Indications for novel intersubtype recombinants were found in 20 patients (6.3 %). Phylogenetic analysis revealed only few and small clusters of local transmission but could document one cluster of CRF02_AG transmission among Belgian MSM.

Conclusions

The extent to which non-B subtypes spread in the native Belgian-Luxembourg population is higher than expected, with 30.5 % of the non-B infections diagnosed in native Western Europeans. These infections resulted from hetero- as well as homosexual transmission. Introduction of non-B variants in the local high at risk population of MSM may lead to new sub-epidemics and/or increased genetic variability and is an evolution that needs to be closely monitored.