Genetic diversity of HIV type 1 along the coastal strip of Kenya

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.

HIV-1 Transmission Patterns Within and Between Risk Groups in Coastal Kenya

HIV-1 transmission patterns within and between populations at different risk of HIV-1 acquisition in Kenya are not well understood. We investigated HIV-1 transmission networks in men who have sex with men (MSM), injecting drug users (IDU), female sex workers (FSW) and heterosexuals (HET) in coastal Kenya. We used maximum-likelihood and Bayesian phylogenetics to analyse new (N = 163) and previously published (N = 495) HIV-1 polymerase sequences collected during 2005-2019. Of the 658 sequences, 131 (20%) were from MSM, 58 (9%) IDU, 109 (17%) FSW, and 360 (55%) HET. Overall, 206 (31%) sequences formed 61 clusters. Most clusters (85%) consisted of sequences from the same risk group, suggesting frequent within-group transmission. The remaining clusters were mixed between HET/MSM (7%), HET/FSW (5%), and MSM/FSW (3%) sequences. One large IDU-exclusive cluster was found, indicating an independent sub-epidemic among this group. Phylodynamic analysis of this cluster revealed a steady increase in HIV-1 infections among IDU since the estimated origin of the cluster in 1987. Our results suggest mixing between high-risk groups and heterosexual populations and could be relevant for the development of targeted HIV-1 prevention programmes in coastal Kenya.

Older age at infection and nulliparity are associated with long-term non-progression in female sex workers infected with non-subtype B HIV-1

Studies have reported on HIV-infected, antiretroviral therapy (ART)-naïve individuals who show minimal disease progression despite prolonged infection. The characteristics of these long-term non-progressors (LTNPs) are not well-characterized in populations predominantly infected with non-subtype B HIV-1. Female sex workers in Mombasa, Kenya who acquired HIV-1 were studied to ascertain immunologic disease progression. Long-term non-progression was defined as an ART-naïve duration of infection ≥7 years and a majority of CD4+ cell counts ≥600 cells/µl with a non-declining CD4+ trend. Correlates of long-term non-progression were determined using multivariable logistic regression. Between February 1993 and March 2014, 332 women acquired HIV-1. Of these, 77 (23%) had ≥7 years of follow-up and 13 (17%) were categorized as LTNPs. Factors associated with long-term non-progression included age >30 years at infection (aOR = 9.41, 95% CI: 1.48–59.86, P = 0.005) and nulliparity (aOR = 20.19, 95% CI: 1.36–299.90, P = 0.03). Each log10 copies/ml increase in viral load (VL) set point was associated with a lower likelihood of being a LTNP (aOR = 0.31, 95% CI: 0.12–0.79, P = 0.01). These findings suggest that age and parity may influence the likelihood of long-term non-progression through mechanisms that are not mediated by the effects of these variables on VL. Future studies should seek to determine whether the associations presented are reproducible.

HIV-1 subtype diversity, transmission networks and transmitted drug resistance amongst acute and early infected MSM populations from Coastal Kenya

Background

HIV-1 molecular epidemiology amongst men who have sex with men (MSM) in sub-Saharan Africa remains not well characterized. We aimed to determine HIV-1 subtype distribution, transmission clusters and transmitted drug resistance (TDR) in acute and early infected MSM from Coastal Kenya.

Methods

Analysis of HIV-1 partial pol sequences from MSM recruited 2005–2017 and sampled within six months of the estimated date of infection. Volunteers were classified as men who have sex with men exclusively (MSME) or with both men and women (MSMW). HIV-1 subtype and transmission clusters were determined by maximum-likelihood phylogenetics. TDR mutations were determined using the Stanford HIV drug resistance database.

Results

Of the 97 volunteers, majority (69%) were MSMW; 74%, 16%, 9% and 1% had HIV-1 subtypes A1, D, C or G, respectively. Overall, 65% formed transmission clusters, with substantial mixing between MSME and MSMW. Majority of volunteer sequences were either not linked to any reference sequence (56%) or clustered exclusively with sequences of Kenyan origin (19%). Eight (8% [95% CI: 4–16]) had at least one TDR mutation against nucleoside (n = 2 [2%]) and/or non-nucleoside (n = 7 [7%]) reverse transcriptase inhibitors. The most prevalent TDR mutation was K103N (n = 5), with sequences forming transmission clusters of two and three taxa each. There were no significant differences in HIV-1 subtype distribution and TDR between MSME and MSMW.

Conclusions

This HIV-1 MSM epidemic was predominantly sub-subtype A1, of Kenyan origin, with many transmission clusters and having intermediate level of TDR. Targeted HIV-1 prevention, early identification and care interventions are warranted to break the transmission cycle amongst MSM from Coastal Kenya.

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.

Prevalence of Transmitted Drug Resistance Mutations in HIV-1-Infected Drug-Naive Patients from Urban and Suburban Regions of Kenya

HIV was first described in Kenya in 1984-1985. Currently, Kenya has an estimated HIV-1 prevalence of 6.2%. With the introduction of antiretroviral drugs, the survival of most HIV patients has been prolonged markedly. However, this is greatly threatened by increasing rates of antiretroviral dug resistance, which may eventually lead to suboptimal treatment outcomes. The objective of this study was to characterize currently occurring antiretroviral drug resistance mutations among drug-naive patients visiting two referral hospitals in Kenya. Using polymerase chain reaction, the HIV protease gene was amplified from blood samples of 63 study participants. The sequences were used to determine HIV-1 subtype and presence/prevalence of mutations associated with resistance to protease inhibitors. Finally, the protease gene was variably measured using Shannon entropy analysis. Analysis of frequency of HIV-1 subtypes revealed subtype A to be the predominant subtype, while the analysis of drug resistance mutations revealed the presence of four minor drug resistance mutations associated weakly with resistance to protease inhibitors. Among these mutations, L33I was the most prevalent mutation. Shannon entropy analysis revealed high genomic variability, especially in region spanning nucleotides 1-55, 113-170, and 205-240. This study warrants the need for dedicated efforts to improve compliance to antiretroviral therapy and reduce transmitted resistance rates, which will greatly ensure the therapeutic efficacy of antiretroviral drugs.

Molecular Epidemiology and Transmission Dynamics of Recent and Long-Term HIV-1 Infections in Rural Western Kenya

Objective

To identify unique characteristics of recent versus established HIV infections and describe sexual transmission networks, we characterized circulating HIV-1 strains from two randomly selected populations of ART-naïve participants in rural western Kenya.

Methods

Recent HIV infections were identified by the HIV-1 subtype B, E and D, immunoglobulin G capture immunoassay (IgG BED-CEIA) and BioRad avidity assays. Genotypic and phylogenetic analyses were performed on the pol gene to identify transmitted drug resistance (TDR) mutations, characterize HIV subtypes and potential transmission clusters. Factors associated with recent infection and clustering were assessed by logistic regression.

Results

Of the 320 specimens, 40 (12.5%) were concordantly identified by the two assays as recent infections. Factors independently associated with being recently infected were age ≤19 years (P = 0.001) and history of sexually transmitted infections (STIs) in the past six months (P = 0.004). HIV subtype distribution differed in recently versus chronically infected participants, with subtype A observed among 53% recent vs. 68% chronic infections (p = 0.04) and subtype D among 26% recent vs. 12% chronic infections (p = 0.012). Overall, the prevalence of primary drug resistance was 1.16%. Of the 258 sequences, 11.2% were in monophyletic clusters of between 2–4 individuals. In multivariate analysis factors associated with clustering included having recent HIV infection P = 0.043 and being from Gem region P = 0.002.

Conclusions

Recent HIV-1 infection was more frequent among 13–19 year olds compared with older age groups, underscoring the ongoing risk and susceptibility of younger persons for acquiring HIV infection. Our findings also provide evidence of sexual networks. The association of recent infections with clustering suggests that early infections may be contributing significant proportions of onward transmission highlighting the need for early diagnosis and treatment as prevention for ongoing prevention. Larger studies are needed to better understand the structure of these networks and subsequently implement and evaluate targeted interventions.

Genotypic Variability of HIV-1 Reverse Transcriptase Gene from Long-Term Antiretroviral-Experienced Patients in Kenya

There is continuous need to track genetic profiles of HIV strains circulating in different geographic settings to hasten vaccine discovery and inform public health and intervention policies. We partially sequenced the reverse transcriptase region of the HIV-1 pol gene from a total of 54 Kenyan patients aged 18-56 years who continued highly active antiretroviral treatment (HAART) for between 8 and 102 months. Subtyping was done using both the JPHMM tool and phylogenetic method. HIV-1 subtype A1 was the predominant strain in circulation, representing 57.4% and 70.4% of all isolates as determined by JPHMM and phylogenetic methods, respectively. Subtypes D (14.8%, 7.4%), C (5.6%, 9.3%), and A2 (0%, 5.6%) were determined at respective prevalence by both methods. JPHMM identified 22.2% of the isolates as recombinants. This surveillance focused on the RT gene and reaffirms the predominance of subtype A and an increasing proportion of recombinant strains in the Kenyan epidemic.

HIV-1 subtype diversity based on envelope C2V3 sequences from Kenyan patients on highly active antiretroviral therapy

There is a continuous need to genetically characterize the HIV strains in circulation in order to assess interventions and inform vaccine discovery. We partially sequenced the envelope C2V3 gene from a total of 59 Kenyan patients on highly active antiretroviral treatment (HAART) and determined HIV subtypes using both the JPHMM subtyping tool and the phylogenetic method. HIV-1 subtype A1 was the predominant strain in circulation, representing 65.5% and 74.5% of all isolates as determined by JPHMM and phylogenetic methods, respectively. Subtypes C and D were the next most prevalent pure strains at 9.1% each by both methods. JPHMM identified 9.1% of the isolates as recombinant. Four isolates had short sequences not covering the entire C2V3 region and were thus not subtyped. From this study, subtype A viruses are still the predominant HIV-1 strains in local circulation in Kenya. Constant surveillance is needed to update molecular trends under continuing HAART scale-up.

HIV diversity and drug resistance from plasma and non-plasma analytes in a large treatment programme in western Kenya

INTRODUCTION

Antiretroviral resistance leads to treatment failure and resistance transmission. Resistance data in western Kenya are limited. Collection of non-plasma analytes may provide additional resistance information.

METHODS

We assessed HIV diversity using the REGA tool, transmitted resistance by the WHO mutation list and acquired resistance upon first-line failure by the IAS-USA mutation list, at the Academic Model Providing Access to Healthcare (AMPATH), a major treatment programme in western Kenya. Plasma and four non-plasma analytes, dried blood-spots (DBS), dried plasma-spots (DPS), ViveST(TM)-plasma (STP) and ViveST-blood (STB), were compared to identify diversity and evaluate sequence concordance.

RESULTS

Among 122 patients, 62 were treatment-naïve and 60 treatment-experienced; 61% were female, median age 35 years, median CD4 182 cells/µL, median viral-load 4.6 log10 copies/mL. One hundred and ninety-six sequences were available for 107/122 (88%) patients, 58/62 (94%) treatment-naïve and 49/60 (82%) treated; 100/122 (82%) plasma, 37/78 (47%) attempted DBS, 16/45 (36%) attempted DPS, 14/44 (32%) attempted STP from fresh plasma and 23/34 (68%) from frozen plasma, and 5/42 (12%) attempted STB. Plasma and DBS genotyping success increased at higher VL and shorter shipment-to-genotyping time. Main subtypes were A (62%), D (15%) and C (6%). Transmitted resistance was found in 1.8% of plasma sequences, and 7% combining analytes. Plasma resistance mutations were identified in 91% of treated patients, 76% NRTI, 91% NNRTI; 76% dual-class; 60% with intermediate-high predicted resistance to future treatment options; with novel mutation co-occurrence patterns. Nearly 88% of plasma mutations were identified in DBS, 89% in DPS and 94% in STP. Of 23 discordant mutations, 92% in plasma and 60% in non-plasma analytes were mixtures. Mean whole-sequence discordance from frozen plasma reference was 1.1% for plasma-DBS, 1.2% plasma-DPS, 2.0% plasma-STP and 2.3% plasma-STB. Of 23 plasma-STP discordances, one mutation was identified in plasma and 22 in STP (p<0.05). Discordance was inversely significantly related to VL for DBS.

CONCLUSIONS

In a large treatment programme in western Kenya, we report high HIV-1 subtype diversity; low plasma transmitted resistance, increasing when multiple analytes were combined; and high-acquired resistance with unique mutation patterns. Resistance surveillance may be augmented by using non-plasma analytes for lower-cost genotyping in resource-limited settings.

Phylodynamics of HIV-1 subtype C epidemic in east Africa

The HIV-1 subtype C accounts for an important fraction of HIV infections in east Africa, but little is known about the genetic characteristics and evolutionary history of this epidemic. Here we reconstruct the origin and spatiotemporal dynamics of the major HIV-1 subtype C clades circulating in east Africa. A large number (n = 1,981) of subtype C pol sequences were retrieved from public databases to explore relationships between strains from the east, southern and central African regions. Maximum-likelihood phylogenetic analysis of those sequences revealed that most (>70%) strains from east Africa segregated in a single regional-specific monophyletic group, here called C_EA. A second major Ethiopian subtype C lineage and a large collection of minor Kenyan and Tanzanian subtype C clades of southern African origin were also detected. A Bayesian coalescent-based method was then used to reconstruct evolutionary parameters and migration pathways of the C_EA African lineage. This analysis indicates that the C_EA clade most probably originated in Burundi around the early 1960s, and later spread to Ethiopia, Kenya, Tanzania and Uganda, giving rise to major country-specific monophyletic sub-clusters between the early 1970s and early 1980s. The results presented here demonstrate that a substantial proportion of subtype C infections in east Africa resulted from dissemination of a single HIV local variant, probably originated in Burundi during the 1960s. Burundi was the most important hub of dissemination of that subtype C clade in east Africa, fueling the origin of new local epidemics in Ethiopia, Kenya, Tanzania and Uganda. Subtype C lineages of southern African origin have also been introduced in east Africa, but seem to have had a much more restricted spread.

HIV type 1 subtype surveillance in central Kenya

Human immunodeficiency virus 1 (HIV-1) infection is characterized by genetic diversity such that specific viral subtypes are predominant in specific geographic areas. To determine circulating subtypes of HIV-1 in different parts of central Kenya, a cross-sectional study was carried out on HIV-1-positive blood samples collected from consenting individuals in eight hospitals of Kenya's central province. Proviral DNA was extracted from peripheral blood mononuclear cells. Polymerase chain reaction and direct sequencing using primers generated from a highly conserved region of HIV-1 env gp41 were carried out. Ninety-six samples were successfully amplified and sequenced. Analysis of the sequences showed that a majority of them belonged to subtype A1 (67/96, 69.8%), followed by subtypes D (18, 18.7%) and C (11/96, 11.5%). Consistent with findings in other parts of Kenya, HIV-1 subtype A1 was the most dominant virus in circulation. Continued surveillance of circulating subtypes of HIV-1 in Kenya is important in determining the evolution of the HIV/AIDS epidemic in Kenya.

HIV type 1 in a rural coastal town in Kenya shows multiple introductions with many subtypes and much recombination

The extent of HIV-1 diversity was examined among patients attending a rural district hospital in a coastal area of Kenya. The pol gene was sequenced in samples from 153 patients. Subtypes were designated using the REGA, SCUEAL, and jpHMM programs. The most common subtype was A1, followed by C and D; A2 and G were also detected. However, a large proportion of the samples was found to be recombinants, which clustered within the pure subtype branches. Phylogeographic analysis of Kilifi sequences compared with those from other regions of Africa showed that while many sequences were closely related to sequences from Kenya, others were most closely related to known sequences from other parts of Africa, including West Africa. Overall, these data indicate that there have been multiple introductions of HIV-1 into this small rural town and surroundings with ongoing diversity being generated by recombination.

Molecular characterization of HIV type 1 among HIV-infected respondents in a cohort being prepared for HIV Phase III vaccine clinical trials, Western Kenya

Kenya is one of the sub-Saharan African countries affected by HIV-1 infection and AIDS. We investigated HIV-1 genetic diversity in 130 individuals from Busia, Bungoma, and Kakamega in western Kenya as part of an HIV-1 vaccine feasibility study in preparation for Phase III efficacy clinical trials. After RNA extraction the partial gag (484 bp) and env (1297 bp) regions were amplified and directly sequenced. Phylogenetic analysis was done using MEGA version 4 and recombinants were identified using the jpHMM tool and phylogenetic analysis. HIV-1 sequences were amplified from 122 of the 130 samples, 118 (90.8%) from the gag region and 78 (60 %) from the env region and 74 samples (56.9%) from both the gag and env regions. Of these sequenced on both regions, 51.4% were subtype A, 9.4% subtype D, 1.4% subtype C, 4.1% subtype G, and 33.7% were discordant and thus possible recombinants, including A1/C, A1/D, A1/A2, and A2/C. The jpHMM tool indicated a further two samples with CD and BD breakpoints within the env gene and one within the gag gene (A1C). An additional sample had an A1D breakpoint in the gag gene, but the envelope was not amplified. HIV-1 subtype diversity in western Kenya should be considered in vaccines designed for clinical trials in this region and this genetic diversity should be continuously monitored.