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

Evaluating whole HIV-1 genome sequence for estimation of incidence and migration in a rural South African community

Background: South Africa has the largest number of people living with HIV (PLWHIV) in the world, with HIV prevalence and transmission patterns varying greatly between provinces. Transmission between regions is still poorly understood, but phylodynamics of HIV-1 evolution can reveal how many infections are attributable to contacts outside a given community. We analysed whole genome HIV-1 genetic sequences to estimate incidence and the proportion of transmissions between communities in Hlabisa, a rural South African community. Methods: We separately analysed HIV-1 for gag, pol, and env genes sampled from 2,503 PLWHIV. We estimated time-scaled phylogenies by maximum likelihood under a molecular clock model. Phylodynamic models were fitted to time-scaled trees to estimate transmission rates, effective number of infections, incidence through time, and the proportion of infections imported to Hlabisa. We also partitioned time-scaled phylogenies with significantly different distributions of coalescent times. Results: Phylodynamic analyses showed similar trends in epidemic growth rates between 1980 and 1990. Model-based estimates of incidence and effective number of infections were consistent across genes. Parameter estimates with gag were generally smaller than those estimated with pol and env. When estimating the proportions of new infections in Hlabisa from immigration or transmission from external sources, our posterior median estimates were 85% (95% credible interval (CI) = 78%–92%) for gag, 62% (CI = 40%–78%) for pol, and 77% (CI = 58%–90%) for env in 2015. Analysis of phylogenetic partitions by gene showed that most close global reference sequences clustered within a single partition. This suggests local evolving epidemics or potential unmeasured heterogeneity in the population. Conclusions: We estimated consistent epidemic dynamic trends for gag, pol and env genes using phylodynamic models. There was a high probability that new infections were not attributable to endogenous transmission within Hlabisa, suggesting high inter-connectedness between communities in rural South Africa.

Genetic Diversity and Acquired Drug Resistance Mutations Detected by Deep Sequencing in Virologic Failures among Antiretroviral Treatment Experienced Human Immunodeficiency Virus-1 Patients in a Pastoralist Region of Ethiopia

Purpose

This study was conducted to investigate the drug resistance mutations and genetic diversity of HIV-1 in ART experienced patients in South Omo, Ethiopia.

Patients and Methods

A cross-sectional study conducted on 253 adult patients attending ART clinics for ≥6 months in South Omo. Samples with VL ≥1000 copies/mL were considered as virological failures (VF) and their reverse transcriptase gene codons 90–234 were sequenced using Illumina MiSeq. MinVar was used for the identification of the subtypes and drug resistance mutations. Phylogenetic tree was constructed by neighbor-joining method using the maximum likelihood model.

Results

The median duration of ART was 51 months and 18.6% (47/253) of the patients exhibited VF. Of 47 viraemic patients, the genome of 41 were sequenced and subtype C was dominant (87.8%) followed by recombinant subtype BC (4.9%), M-09-CPX (4.9) and BF1 (2.4%). Of 41 genotyped subjects, 85.4% (35/41) had at least one ADR mutation. Eighty-one percent (33/41) of viraemic patients harbored NRTI resistance mutations, and 48.8% (20/41) were positive for NNRTI resistance mutations, with 43.9% dual resistance mutations. Among NRTI resistance mutations, M184V (73.2%), K219Q (63.4%) and T215 (56.1%) complex were the most mutated positions, while the most common NNRTI resistance mutations were K103N (24.4%), K101E, P225H and V108I 7.5% each. Active tuberculosis (aOR=13, 95% CI= 3.46–29.69), immunological failure (aOR=3.61, 95% CI=1.26–10.39), opportunistic infections (aOR=8.39, 95% CI= 1.75–40.19), and poor adherence were significantly associated with virological failure, while rural residence (aOR 2.37; 95% CI: 1.62–9.10, P= 0.05), immunological failures (aOR 2.37; 95% CI: 1.62–9.10, P= 0.05) and high viral load (aOR 16; 95% CI: 5.35 51.59, P <0.001) were predictors of ADR mutation among the ART experienced and viraemic study subjects.

Conclusion

The study revealed considerable prevalence of VF and ADR mutation with the associated risk indicators. Regular virological monitoring and drug resistance genotyping methods should be implemented for better ART treatment outcomes of the nation.

Phenotypic and Genotypic Co-receptor Tropism Testing in HIV-1 Epidemic Region of Tanzania Where Multiple Non-B Subtypes Co-circulate

HIV human immunodeficiency virus type I (HIV-1) entry inhibitor potency is dependent on viral co-receptor tropisms and thereby tropism determination is clinically important. However, phenotypic tropisms of HIV-1 non-B subtypes have been poorly investigated and the genotypic prediction algorithms remain insufficiently validated. To clarify this issue, we recruited 52 treatment-naïve, HIV-1-infected patients in Tanzania, where multiple HIV-1 non-B subtypes co-circulate. Sequence analysis of 93 infectious envelope clones isolated from their plasma viral RNA revealed the co-circulation of subtypes A1, C, D, and inter-subtype recombinant forms (isRFs). Phenotypic tropism assays revealed that lentivirus reporters pseudotyped with 75 (80.6%) and 5 (5.4%) envelope clones could establish infection toward U87.CD4 cells expressing CCR5 (R5) and CXCR4 (X4), respectively; whereas the remaining 13 (14%) clones could infect both cells. Genotypic analyses by widely used algorithms including V3 net charge, Geno2pheno, WebPSSM, and PhenoSeq showed that almost all phenotypic X4-tropic clones and only 15 of 75 phenotypic R5-tropic clones were concordantly predicted. However, the remaining 60 phenotypic R5-tropic clones were discordantly predicted by at least one algorithm. In particular, 2 phenotypic R5-tropic clones were discordantly predicted by all algorithms tested. Taken together, the results demonstrate the limitation of currently available genotypic algorithms for predicting co-receptor inference among co-circulating multiple non-B subtypes and emerging isRFs. Also, the phenotypic tropism dataset presented here could be valuable for retraining of the widely used genotypic prediction algorithms to enhance their performance.

Cytotoxic Lymphocytes Target HIV-1 Gag Through Granzyme M-Mediated Cleavage

Untreated HIV-1 infection leads to a slow decrease in CD4⁺ T cell lymphocytes over time resulting in increased susceptibility to opportunistic infections (acquired immunodeficiency syndrome, AIDS) and ultimately death of the infected individual. Initially, the host's immune response controls the infection, but cannot eliminate the HIV-1 from the host. Cytotoxic lymphocytes are the key effector cells in this response and can mediate crucial antiviral responses through the release of a set of proteases called granzymes towards HIV-1-infected cells. However, little is known about the immunological molecular mechanisms by which granzymes could control HIV-1. Since we noted that HIV-1 subtype C (HIV-1C) Gag with the tetrapeptide insertion PYKE contains a putative granzyme M (GrM) cleavage site (KEPL) that overlaps with the PYKE insertion, we analyzed the proteolytic activity of GrM towards Gag. Immunoblot analysis showed that GrM could cleave Gag proteins from HIV-1B and variants from HIV-1C of which the Gag-PYKE variant was cleaved with extremely high efficiency. The main cleavage site was directly after the insertion after leucine residue 483. GrM-mediated cleavage of Gag was also observed in co-cultures using cytotoxic lymphocytes as effector cells and this cleavage could be inhibited by a GrM inhibitor peptide. Altogether, our data indicate towards a noncytotoxic immunological mechanism by which GrM-positive cytotoxic lymphocytes target the HIV-1 Gag protein within infected cells to potentially control HIV-1 infection. This mechanism could be exploited in new therapeutic strategies to treat HIV-1-infected patients to improve immunological control of the infection.

HIV-1 p24Gag adaptation to modern and archaic HLA-allele frequency differences in ethnic groups contributes to viral subtype diversification

Pathogen-driven selection and past interbreeding with archaic human lineages have resulted in differences in human leukocyte antigen (HLA)-allele frequencies between modern human populations. Whether or not this variation affects pathogen subtype diversification is unknown. Here we show a strong positive correlation between ethnic diversity in African countries and both human immunodeficiency virus (HIV)-1 p24gag and subtype diversity. We demonstrate that ethnic HLA-allele differences between populations have influenced HIV-1 subtype diversification as the virus adapted to escape common antiviral immune responses. The evolution of HIV Subtype B (HIV-B), which does not appear to be indigenous to Africa, is strongly affected by immune responses associated with Eurasian HLA variants acquired through adaptive introgression from Neanderthals and Denisovans. Furthermore, we show that the increasing and disproportionate number of HIV-infections among African Americans in the USA drive HIV-B evolution towards an Africa-centric HIV-1 state. Similar adaptation of other pathogens to HLA variants common in affected populations is likely.

Increased HIV-1 pretreatment drug resistance with consistent clade homogeneity among ART-naive HIV-1 infected individuals in Ethiopia

BACKGROUND

The development of pretreatment drug resistance (PDR) is becoming an obstacle to the success of antiretroviral therapy (ART). Besides, data from developing settings including Ethiopia is still limited. Therefore, this study was aimed to assess HIV-1 genetic diversity and PDR mutations among ART-naive recently diagnosed HIV-1 infected individuals in Addis Ababa, Ethiopia.

METHODS

Institutional based cross-sectional study was conducted from June to December 2018 in Addis Ababa among ART-naive recently diagnosed individuals. Partial HIV-1 pol region covering the entire protease (PR) and partial reverse transcriptase (RT) regions of 51 samples were amplified and sequenced using an in-house assay. Drug resistance mutations were examined using calibrated population resistance (CPR) tool version 6.0 from the Stanford HIV drug resistance database and the International Antiviral Society-USA (IAS-USA) 2019 mutation list.

RESULTS

According to both algorithms used, 9.8% (5/51) of analyzed samples had at least one PDR Mutation. PDR mutations to Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs) were the most frequently detected (7.8% and 9.8%, according to the CPR tool and IAS-USA algorithm, respectively). The most frequently observed NNRTIs-associated mutations common to both algorithms were K103N (2%), Y188L (2%), K101E (2%), and V106A (2%), while E138A (2%) was observed according to IAS-USA only. Y115F and M184V (mutations that confer resistance to NRTIs) dual mutations were detected according to both criteria in a single study participant (2%). PDR mutation to protease inhibitors was found to be low (only G73S; 2% according to the CPR tool). Phylogenetic analysis showed that 98% (50/51) of the study participants were infected with HIV-1C virus while one individual (2%) was infected with HIV-1A1 virus.

CONCLUSIONS

This study showed an increased level of PDR and persistence HIV-1C clade homogeneity after 15 years of the rollout of ART and 3 decades of HIV-1C circulation in Ethiopia, respectively. Therefore, we recommend routine baseline genotypic drug resistance testing for all newly diagnosed HIV infected patients before initiating treatment. This will aid the selection of appropriate therapy in achieving the 90% of patients having an undetectable viral load in consonance with the UN target.

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

MiDRMpol: A High-Throughput Multiplexed Amplicon Sequencing Workflow to Quantify HIV-1 Drug Resistance Mutations against Protease, Reverse Transcriptase, and Integrase Inhibitors

The detection of drug resistance mutations (DRMs) in minor viral populations is of potential clinical importance. However, sophisticated computational infrastructure and competence for analysis of high-throughput sequencing (HTS) data lack at most diagnostic laboratories. Thus, we have proposed a new pipeline, MiDRMpol, to quantify DRM from the HIV-1 pol region. The gag-vpu region of 87 plasma samples from HIV-infected individuals from three cohorts was amplified and sequenced by Illumina HiSeq2500. The sequence reads were adapter-trimmed, followed by analysis using in-house scripts. Samples from Swedish and Ethiopian cohorts were also sequenced by Sanger sequencing. The pipeline was validated against the online tool PASeq (Polymorphism Analysis by Sequencing). Based on an error rate of <1%, a value of >1% was set as reliable to consider a minor variant. Both pipelines detected the mutations in the dominant viral populations, while discrepancies were observed in minor viral populations. In five HIV-1 subtype C samples, minor mutations were detected at the <5% level by MiDRMpol but not by PASeq. MiDRMpol is a computationally as well as labor efficient bioinformatics pipeline for the detection of DRM from HTS data. It identifies minor viral populations (<20%) of DRMs. Our method can be incorporated into large-scale surveillance of HIV-1 DRM.

A viral genome wide association study and genotypic resistance testing in patients failing first line antiretroviral therapy in the first large countrywide Ethiopian HIV cohort

BACKGROUND

Antiretroviral therapy (ART) was rolled-out in Ethiopia in 2005, but there are no reports on outcome of ART and human immunodeficiency virus drug resistance (HIVDR) at national level. We described acquired drug resistance mutations in pol gene and performed a viral genome wide association study in virologic treatment failure patients who started first line ART during 2009-2011 in the first large countrywide HIV cohort in Ethiopia.

METHODS

The outcome of tenofovir (TDF)- and zidovudine (ZDV)-based ART was defined in 874 ART naïve patients using the on-treatment (OT) and intention-to-treat (ITT) analyses. Genotypic resistance testing was done in patients failing ART (> 1000 copies/ml) at month 6 and 12. Near full-length genome sequencing (NFLG) was used to assess amino acid changes in HIV-1 gag, pol, vif, vpr, tat, vpu, and nef genes between paired baseline and month 6 samples.

RESULTS

High failure rates were found in ITT analysis at month 6 and 12 (23.3%; 33.9% respectively). Major nucleoside and non-nucleoside reverse transcriptase (NRTI/NNRTI) drug resistance mutations were detected in most failure patients at month 6 (36/47; 77%) and month 12 (20/30; 67%). A high rate of K65R was identified only in TDF treated patients (35.7%; 50.0%, respectively). No significant difference was found in failure rate or extent of HIVDR between TDF- and ZDV- treated patients. All target regions of interest for HIVDR were described by NFLG in 16 patients tested before initiation of ART and at month 6.

CONCLUSION

In this first Ethiopian national cohort, a high degree of HIVDR was seen among ART failure patients, independent on whether TDF- or ZDV was given. However, the major reason to ART failure was lost-to-follow-up rather than virologic failure. Our NFLG assay covered all relevant target genes for antiretrovirals and is an attractive alternative for HIVDR surveillance.

HIV-1 Subtype C with PYxE Insertion Has Enhanced Binding of Gag-p6 to Host Cell Protein ALIX and Increased Replication Fitness

Human immunodeficiency virus type 1 subtype C (HIV-1C) has a natural deletion of a YPxL motif in its Gag-p6 late domain. This domain mediates the binding of Gag to host cell protein ALIX and subsequently facilitates viral budding. In a subset of HIV-1C-infected individuals, the tetrapeptide insertion PYxE has been identified at the deleted YPxL motif site. Here, we report the consequences of PYxE insertion on the interaction with ALIX and the relevance regarding replication fitness and drug sensitivity. In our three HIV-1C cohorts, PYKE and PYQE were most prevalent among PYxE variants. Through in silico predictions and in vitro experiments, we showed that HIV-1C Gag has an increased binding to ALIX when the PYxE motif is present. To go more into the clinical relevance of the PYxE insertion, we obtained patient-derived gag-pol sequences from HIV-1C_PYxEi viruses and inserted them in a reference HIV-1 sequence. Viral growth was increased, and the sensitivity to the protease inhibitor (PI) lopinavir (LPV) and nucleoside reverse transcriptase inhibitor tenofovir alafenamide (TAF) was decreased for some of the HIV-1C PYxE variants compared to that of wild-type variants. Our data suggest that PYxE insertion in Gag restores the ability of Gag to bind ALIX and correlates with enhanced viral fitness in the absence or presence of LPV and TAF. The high prevalence and increased replication fitness of the HIV-1C virus with PYxE insertion indicates the clinical importance of these viral variants.IMPORTANCE Genomic differences within HIV-1 subtypes is associated with various degrees of viral spread, disease progression, and clinical outcome. Viral budding is essential in the HIV-1 life cycle and mainly mediated through the interaction of Gag with host proteins. Two motifs within Gag-p6 mediate binding of host cell proteins and facilitate budding. HIV-1C has a natural deletion of one of these two motifs, resulting in an inability to bind to host cell protein ALIX. Previously, we have identified a tetrapeptide (PYxE) insertion at this deleted motif site in a subset of HIV-1C patients. Here, we report the incidence of PYxE insertions in three different HIV-1C cohorts, and the insertion restores the binding of Gag to ALIX. It also increases viral growth even in the presence of the antiretroviral drugs lopinavir and tenofovir alafenamide. Hence, PYxE insertion in HIV-1C might be biologically relevant for viruses and clinically significant among patients.

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

Recent increased identification and transmission of HIV-1 unique recombinant forms in Sweden

A temporal increase in non-B subtypes has earlier been described in Sweden by us and we hypothesized that this increased viral heterogeneity may become a hotspot for the development of more complex and unique recombinant forms (URFs) if the epidemics converge. In the present study, we performed subtyping using four automated tools and phylogenetic analysis by RAxML of pol gene sequences (n = 5246) and HIV-1 near full-length genome (HIV-NFLG) sequences (n = 104). A CD4⁺ T-cell decline trajectory algorithm was used to estimate time of HIV infection. Transmission clusters were identified using the family-joining method. The analysis of HIV-NFLG and pol gene described 10.6% (11/104) and 2.6% (137/5246) of the strains as URFs, respectively. An increasing trend of URFs was observed in recent years by both approaches (p = 0·0082; p < 0·0001). Transmission cluster analysis using the pol gene of all URFs identified 14 clusters with two to eight sequences. Larger transmission clusters of URFs (BF1 and 01B) were observed among MSM who mostly were sero-diagnosed in recent time. Understanding the increased appearance and transmission of URFs in recent years could have importance for public health interventions and the use of HIV-NFLG would provide better statistical support for such assessments.

Using nearly full-genome HIV sequence data improves phylogeny reconstruction in a simulated epidemic

HIV molecular epidemiology studies analyse viral pol gene sequences due to their availability, but whole genome sequencing allows to use other genes. We aimed to determine what gene(s) provide(s) the best approximation to the real phylogeny by analysing a simulated epidemic (created as part of the PANGEA_HIV project) with a known transmission tree. We sub-sampled a simulated dataset of 4662 sequences into different combinations of genes (gag-pol-env, gag-pol, gag, pol, env and partial pol) and sampling depths (100%, 60%, 20% and 5%), generating 100 replicates for each case. We built maximum-likelihood trees for each combination using RAxML (GTR + Γ), and compared their topologies to the corresponding true tree’s using CompareTree. The accuracy of the trees was significantly proportional to the length of the sequences used, with the gag-pol-env datasets showing the best performance and gag and partial pol sequences showing the worst. The lowest sampling depths (20% and 5%) greatly reduced the accuracy of tree reconstruction and showed high variability among replicates, especially when using the shortest gene datasets. In conclusion, using longer sequences derived from nearly whole genomes will improve the reliability of phylogenetic reconstruction. With low sample coverage, results can be highly variable, particularly when based on short sequences.