A cluster of HIV type 1 subtype C sequences from Ethiopia, observed in full genome analysis, is not sustained in subgenomic regions

Collaborative Mining of Whole Genome Sequences for Intelligent HIV-1 Sub-Strain(s) Discovery

BACKGROUND

Effective global antiretroviral vaccines and therapeutic strategies depend on the diversity, evolution, and epidemiology of their various strains as well as their transmission and pathogenesis. Most viral disease-causing particles are clustered into a taxonomy of subtypes to suggest pointers toward nucleotide-specific vaccines or therapeutic applications of clinical significance sufficient for sequence-specific diagnosis and homologous viral studies. These are very useful to formulate predictors to induce cross-resistance to some retroviral control drugs being used across study areas.

OBJECTIVE

This research proposed a collaborative framework of hybridized (Machine Learning and Natural Language Processing) techniques to discover hidden genome patterns and feature predictors, for HIV-1 genome sequences mining.

METHOD

630 human HIV-1 genome sequences above 8500 bps were excavated from the National Center for Biotechnology Information (NCBI) database (https://www.ncbi.nlm.nih.gov) for 21 countries across different continents, Antarctica exempt. These sequences were transformed and learned using a self-organizing map (SOM). To discriminate emerging/new sub-strain(s), the HIV-1 reference genome was included as part of the input isolates/samples during the training. After training the SOM, component planes defining pattern clusters of the input datasets were generated, for cognitive knowledge mining and subsequent labelling of the datasets. Additional genome features including dinucleotide transmission recurrences, codon recurrences, and mutation recurrences, were finally extracted from the raw genomes to construct output classification targets for supervised learning.

RESULTS

SOM training explains the inherent pattern diversity of HIV-1 genomes as well as inter- and intra-country transmissions in which mobility might play an active role, as corroborated by literature. Nine sub-strains were discovered after disassembling the SOM correlation hunting matrix space attributed to disparate clusters. Cognitive knowledge mining separated similar pattern clusters bounded by a certain degree of correlation range, discovered by the SOM. A Kruskal-Wallis rank-sum test and Wilcoxon rank-sum test showed statistically significant variations in dinucleotide, codon, and mutation patterns.

CONCLUSION

Results of the discovered sub-strains and response clusters visualizations corroborate existing literature, with significant haplotype variations. The proposed framework would assist in the development of decision support systems for easy contact tracing, infectious disease surveillance, and studying the progressive evolution of the reference HIV-1 genome.

Factors influencing HIV-1 phylogenetic clustering

PURPOSE OF REVIEW

A major goal of public health in relation to HIV/AIDS is to prevent new transmissions in communities. Phylogenetic techniques have improved our understanding of the structure and dynamics of HIV transmissions. However, there is still no consensus about phylogenetic methodology, sampling coverage, gene target and/or minimum fragment size.

RECENT FINDINGS

Several studies use a combined methodology, which includes both a genetic or patristic distance cut-off and a branching support threshold to identify phylogenetic clusters. However, the choice about these thresholds remains an inherently subjective process, which affects the results of these studies. There is still a lack of consensus about the genomic region and the size of fragments that should be used, although there seems to be emerging a consensus that using longer segments, allied with the use of a realistic model of evolution and a codon alignment, increases the likelihood of inferring true transmission clusters. The pol gene is still the most used genomic region, but recent studies have suggested that whole genomes and/or sequences from nef and gp41 are also good targets for cluster reconstruction.

SUMMARY

The development and application of standard methodologies for phylogenetic clustering analysis will advance our understanding of factors associated with HIV transmission. This will lead to the design of more precise public health interventions.

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.

Importance of Viral Sequence Length and Number of Variable and Informative Sites in Analysis of HIV Clustering

To improve the methodology of HIV cluster analysis, we addressed how analysis of HIV clustering is associated with parameters that can affect the outcome of viral clustering. The extent of HIV clustering and tree certainty was compared between 401 HIV-1C near full-length genome sequences and subgenomic regions retrieved from the LANL HIV Database. Sliding window analysis was based on 99 windows of 1,000 bp and 45 windows of 2,000 bp. Potential associations between the extent of HIV clustering and sequence length and the number of variable and informative sites were evaluated. The near full-length genome HIV sequences showed the highest extent of HIV clustering and the highest tree certainty. At the bootstrap threshold of 0.80 in maximum likelihood (ML) analysis, 58.9% of near full-length HIV-1C sequences but only 15.5% of partial pol sequences (ViroSeq) were found in clusters. Among HIV-1 structural genes, pol showed the highest extent of clustering (38.9% at a bootstrap threshold of 0.80), although it was significantly lower than in the near full-length genome sequences. The extent of HIV clustering was significantly higher for sliding windows of 2,000 bp than 1,000 bp. We found a strong association between the sequence length and proportion of HIV sequences in clusters, and a moderate association between the number of variable and informative sites and the proportion of HIV sequences in clusters. In HIV cluster analysis, the extent of detectable HIV clustering is directly associated with the length of viral sequences used, as well as the number of variable and informative sites. Near full-length genome sequences could provide the most informative HIV cluster analysis. Selected subgenomic regions with a high extent of HIV clustering and high tree certainty could also be considered as a second choice.

Antiviral activity and in vitro mutation development pathways of MK-6186, a novel nonnucleoside reverse transcriptase inhibitor

MK-6186 is a novel nonnucleoside reverse transcriptase inhibitor (NNRTI) which displays subnanomolar potency against wild-type (WT) virus and the two most prevalent NNRTI-resistant RT mutants (K103N and Y181C) in biochemical assays. In addition, it showed excellent antiviral potency against K103N and Y181C mutant viruses, with fold changes (FCs) of less than 2 and 5, respectively. When a panel of 12 common NNRTI-associated mutant viruses was tested with MK-6186, only 2 relatively rare mutants (Y188L and V106I/Y188L) were highly resistant, with FCs of >100, and the remaining viruses showed FCs of <10. Furthermore, a panel of 96 clinical virus isolates with NNRTI resistance mutations was evaluated for susceptibility to NNRTIs. The majority (70%) of viruses tested displayed resistance to efavirenz (EFV), with FCs of >10, whereas only 29% of the mutant viruses displayed greater than 10-fold resistance to MK-6186. To determine whether MK-6186 selects for novel resistance mutations, in vitro resistance selections were conducted with one isolate each from subtypes A, B, and C under low-multiplicity-of-infection (MOI) conditions. The results showed a unique mutation development pattern in which L234I was the first mutation to emerge in the majority of the experiments. In resistance selection under high-MOI conditions with subtype B virus, V106A was the dominant mutation detected in the breakthrough viruses. More importantly, mutant viruses selected by MK-6186 showed FCs of <10 against EFV or etravirine (ETR), and the mutant viruses containing mutations selected by EFV or ETR were sensitive to MK-6186 (FCs of <10).

Antiviral activity of MK-4965, a novel nonnucleoside reverse transcriptase inhibitor

Nonnucleoside reverse transcriptase inhibitors (NNRTIs) are the mainstays of therapy for the treatment of human immunodeficiency virus type 1 (HIV-1) infections. However, the effectiveness of NNRTIs can be hampered by the development of resistance mutations which confer cross-resistance to drugs in the same class. Extensive efforts have been made to identify new NNRTIs that can suppress the replication of the prevalent NNRTI-resistant viruses. MK-4965 is a novel NNRTI that possesses both diaryl ether and indazole moieties. The compound displays potency at subnanomolar concentrations against wild-type (WT), K103N, and Y181C reverse transcriptase (RT) in biochemical assays. MK-4965 is also highly potent against the WT virus and two most prevalent NNRTI-resistant viruses (viruses that harbor the K103N or the Y181C mutation), against which it had 95% effective concentrations (EC(95)s) of <30 nM in the presence of 10% fetal bovine serum. The antiviral EC(95) of MK-4965 was reduced approximately four- to sixfold when it was tested in 50% human serum. Moreover, MK-4965 was evaluated with a panel of 15 viruses with NNRTI resistance-associated mutations and showed a superior mutant profile to that of efavirenz but not to that of etravirine. MK-4965 was similarly effective against various HIV-1 subtypes and viruses containing nucleoside reverse transcriptase inhibitor or protease inhibitor resistance-conferring mutations. A two-drug combination study showed that the antiviral activity of MK-4965 was nonantagonistic with each of the 18 FDA-licensed drugs tested vice versa in the present study. Taken together, these in vitro data show that MK-4965 possesses the desired properties for further development as a new NNRTI for the treatment of HIV-1 infection.

High replication fitness and transmission efficiency of HIV-1 subtype C from India: Implications for subtype C predominance

HIV-1 subtype C has been the predominant subtype throughout the course of the HIV-1 epidemic in India regardless of the geographic region of the country. In an effort to understand the mechanism of subtype C predominance in this country, we have investigated the in vitro replication fitness and transmission efficiency of HIV-1 subtypes A and C from India. Using a dual infection growth competition assay, we found that primary HIV-1 subtype C isolates had higher overall relative fitness in PBMC than subtype A primary isolates. Moreover, in an ex vivo cervical tissue derived organ culture, subtype C isolates displayed higher transmission efficiency across cervical mucosa than subtype A isolates. We found that higher fitness of subtype C was not due to a trans effect exerted by subtype C infected PBMC. A half genome A/C recombinant clone in which the 3' half of the viral genome of subtype A was replaced with the corresponding subtype C3' half, had similar replicative fitness as the parental subtype A. These results suggest that the higher replication fitness and transmission efficiency of subtype C virus compared to subtype A virus from India is most probably not due to the envelope gene alone and may be due to genes present within the 5' half of the viral genome or to a more complex interaction between the genes located within the two halves of the viral genome. These data provide a model to explain the asymmetric distribution of subtype C over other subtypes in India.

Full-length HIV type 1 proviral sequencing of 10 highly exposed women from Nairobi, Kenya reveals a high proportion of intersubtype recombinants

Phylogenetic analysis has revealed that the current HIV/AIDS pandemic consists of a multitude of different viral clades and recombinant viruses. The predominant circulating HIV-1 clade in Kenya is A1; however, Kenya borders countries where different subtypes are prominent, making Kenya a likely location for recombination. Previous studies have reported significant differences in the proportions of sequences in Kenya that are intersubtype recombinants. Studies that performed sequence-based typing on multiple HIV-1 genomic regions or full-length sequences found higher rates of recombination than those that examined a single gene or gene fragment. In this study, we describe full-length HIV-1 proviral sequence-based genotyping after limited peripheral blood mononuclear cell (PBMC) coculture. Ten subjects from a highly exposed cohort located in Nairobi, Kenya were examined. Pairwise comparison found minimal difference between sequences generated directly from patient PBMC DNA compared to sequences from cocultured PBMC DNA. Of the 10 full-length HIV-1 sequences examined, five were nonrecombinant clade A1, while the other five were unique intersubtype recombinants. Although this frequency of recombination is higher than previously described in Kenya, this finding is in agreement with previous full-length sequence data. Interestingly, although all the nonrecombinant sequences were clade A1, not all the recombinant sequences contained a clade A1 sequence.

Variable contexts and levels of hypermutation in HIV-1 proviral genomes recovered from primary peripheral blood mononuclear cells

APOBEC-mediated cytidine deamination of HIV-1 genomes during reverse transcription has been shown to be a potent mechanism of host restriction for HIV-1 infection ex vivo and in vitro. However, this defense system can be overcome by the viral protein Vif. Unlike other mechanisms of host restriction, the APOCEC-Vif interaction leaves an imprint on integrated proviruses in the form of G-->A hypermutation. In the current work we systematically studied levels, contexts, and patterns of HIV-1 hypermutation in vivo. The analysis of 24 full-genome HIV-1 sequences retrieved from primary PBMCs, representing infections with several HIV-1 clades, and the inclusion of 7 cognate pairs of hypermutated/non-hypermutated sequences derived from the same patient sample, provided a comprehensive view of the characteristics of APOBEC-mediated restriction in vivo. Levels of hypermutation varied nearly 5-fold among the studied proviruses. GpG motifs were most frequently affected (22/24 proviruses). Levels of hypermutation varied across the genome. The reported "twin peak" pattern of hypermutation was observed in 18/24 hypermutants, but the remainder exhibited singular non-conforming patterns. These data suggest considerable complexity in the interplay of host restriction and viral defense during HIV-1 infection.

Extensive intrasubtype recombination in South African human immunodeficiency virus type 1 subtype C infections

Recombinant human immunodeficiency virus type 1 (HIV-1) strains containing sequences from different viral genetic subtypes (intersubtype) and different lineages from within the same subtype (intrasubtype) have been observed. A consequence of recombination can be the distortion of the phylogenetic signal. Several intersubtype recombinants have been identified; however, less is known about the frequency of intrasubtype recombination. For this study, near-full-length HIV-1 subtype C genomes from 270 individuals were evaluated for the presence of intrasubtype recombination. A sliding window schema (window, 2 kb; step, 385 bp) was used to partition the aligned sequences. The Shimodaira-Hasegawa test detected significant topological incongruence in 99.6% of the comparisons of the maximum-likelihood trees generated from each sequence partition, a result that could be explained by recombination. Using RECOMBINE, we detected significant levels of recombination using five random subsets of the sequences. With a set of 23 topologically consistent sequences used as references, bootscanning followed by the interactive informative site test defined recombination breakpoints. Using two multiple-comparison correction methods, 47% of the sequences showed significant evidence of recombination in both analyses. Estimated evolutionary rates were revised from 0.51%/year (95% confidence interval [CI], 0.39 to 0.53%) with all sequences to 0.46%/year (95% CI, 0.38 to 0.48%) with the putative recombinants removed. The timing of the subtype C epidemic origin was revised from 1961 (95% CI, 1947 to 1962) with all sequences to 1958 (95% CI, 1949 to 1960) with the putative recombinants removed. Thus, intrasubtype recombinants are common within the subtype C epidemic and these impact analyses of HIV-1 evolution.

Biologic and genetic characterization of a panel of 60 human immunodeficiency virus type 1 isolates, representing clades A, B, C, D, CRF01_AE, and CRF02_AG, for the development and assessment of candidate vaccines

A critical priority for human immunodeficiency virus type 1 (HIV-1) vaccine development is standardization of reagents and assays for evaluation of immune responses elicited by candidate vaccines. To provide a panel of viral reagents from multiple vaccine trial sites, 60 international HIV-1 isolates were expanded in peripheral blood mononuclear cells and characterized both genetically and biologically. Ten isolates each from clades A, B, C, and D and 10 isolates each from CRF01_AE and CRF02_AG were prepared from individuals whose HIV-1 infection was evaluated by complete genome sequencing. The main criterion for selection was that the candidate isolate was pure clade or pure circulating recombinant. After expansion in culture, the complete envelope (gp160) of each isolate was verified by sequencing. The 50% tissue culture infectious dose and p24 antigen concentration for each viral stock were determined; no correlation between these two biologic parameters was found. Syncytium formation in MT-2 cells and CCR5 or CXCR4 coreceptor usage were determined for all isolates. Isolates were also screened for neutralization by soluble CD4, a cocktail of monoclonal antibodies, and a pool of HIV-1-positive patient sera. The panel consists of 49 nonsyncytium-inducing isolates that use CCR5 as a major coreceptor and 11 syncytium-inducing isolates that use only CXCR4 or both coreceptors. Neutralization profiles suggest that the panel contains both neutralization-sensitive and -resistant isolates. This collection of HIV-1 isolates represents the six major globally prevalent strains, is exceptionally large and well characterized, and provides an important resource for standardization of immunogenicity assessment in HIV-1 vaccine trials.

HLA-A and -B allele expression and ability to develop anti-Gag cross-clade responses in subtype C HIV-1–infected Ethiopians

A cohort of 35 human immunodeficiency virus type 1 (HIV-1) subtype C-infected Ethiopians was studied to define the HLA phenotype in all 35 subjects and highly conserved Gag protein regions involved in cross-clade cell-mediated immunity. Full-length Gag virus sequences were determined in 15 individuals. CD8 cell-mediated immune responses were detected by interferon-gamma ELISpot assay. HLA-A*03, -B*49, and -B*57 allelic frequencies were relatively higher than in other African populations. Anti-p17 (aa 1-60) CD8+ were detectable in the highest number of individuals. Anti-p17 (aa 1-60 and 51-110) cross-clade responses against subtype B and C were detected in 50% of the tested subjects. The p24 KF11 (aa 162-172) epitope was found to be immunodominant among the HLA-B*5703--positive individuals. These data represent the first report of correlating HLA phenotype and HIV-specific cell-mediated immune responses among infected Ethiopians and may be useful in designing cytotoxic T lymphocyte-inducing vaccines for this part of Africa.