COMET: adaptive context-based modeling for ultrafast HIV-1 subtype identification

Viral sequence classification has wide applications in clinical, epidemiological, structural and functional categorization studies. Most existing approaches rely on an initial alignment step followed by classification based on phylogenetic or statistical algorithms. Here we present an ultrafast alignment-free subtyping tool for human immunodeficiency virus type one (HIV-1) adapted from Prediction by Partial Matching compression. This tool, named COMET, was compared to the widely used phylogeny-based REGA and SCUEAL tools using synthetic and clinical HIV data sets (1 090 698 and 10 625 sequences, respectively). COMET''s sensitivity and specificity were comparable to or higher than the two other subtyping tools on both data sets for known subtypes. COMET also excelled in detecting and identifying new recombinant forms, a frequent feature of the HIV epidemic. Runtime comparisons showed that COMET was almost as fast as USEARCH. This study demonstrates the advantages of alignment-free classification of viral sequences, which feature high rates of variation, recombination and insertions/deletions. COMET is free to use via an online interface.     

Testing the hypothesis of a recombinant origin of human immunodeficiency virus type 1 subtype E

The human immunodeficiency virus type 1 (HIV-1) epidemic in Southeast Asia has been largely due to the emergence of clade E (HIV-1E). It has been suggested that HIV-1E is derived from a recombinant lineage of subtype A (HIV-1A) and subtype E, with multiple breakpoints along the E genome. We obtained complete genome sequences of clade E viruses from Thailand (93TH057 and 93TH065) and from the Central African Republic (90CF11697 and 90CF4071), increasing the total number of HIV-1E complete genome sequences available to seven. Phylogenetic analysis of complete genomes showed that subtypes A and E are themselves monophyletic, although together they also form a larger monophyletic group. The apparent phylogenetic incongruence at different regions of the genome that was previously taken as evidence of recombination is shown to be not statistically significant. Furthermore, simulations indicate that bootscanning and pairwise distance results, previously used as evidence for recombination, can be misleading, particularly when there are differences in substitution or evolutionary rates across the genomes of different subtypes. Taken jointly, our analyses suggest that there is inadequate support for the hypothesis that subtype E variants are derived from a recombinant lineage. In contrast, many other HIV strains claimed to have a recombinant origin, including viruses for which only a single parental strain was employed for analysis, do indeed satisfy the statistical criteria we propose. Thus, while intersubtype recombinant HIV strains are indeed circulating, the criteria for assigning a recombinant origin to viral structures should include statistical testing of alternative hypotheses to avoid inappropriate assignments that would obscure the true evolutionary properties of these viruses.     

Comparing mutational pathways to lopinavir resistance in HIV-1 subtypes B versus C

Although combination antiretroviral therapies seem to be effective at controlling HIV-1 infections regardless of the viral subtype, there is increasing evidence for subtype-specific drug resistance mutations. The order and rates at which resistance mutations accumulate in different subtypes also remain poorly understood. Most of this knowledge is derived from studies of subtype B genotypes, despite not being the most abundant subtype worldwide. Here, we present a methodology for the comparison of mutational networks in different HIV-1 subtypes, based on Hidden Conjunctive Bayesian Networks (H-CBN), a probabilistic model for inferring mutational networks from cross-sectional genotype data. We introduce a Monte Carlo sampling scheme for learning H-CBN models for a larger number of resistance mutations and develop a statistical test to assess differences in the inferred mutational networks between two groups. We apply this method to infer the temporal progression of mutations conferring resistance to the protease inhibitor lopinavir in a large cross-sectional cohort of HIV-1 subtype C genotypes from South Africa, as well as to a data set of subtype B genotypes obtained from the Stanford HIV Drug Resistance Database and the Swiss HIV Cohort Study. We find strong support for different initial mutational events in the protease, namely at residue 46 in subtype B and at residue 82 in subtype C. The inferred mutational networks for subtype B versus C are significantly different sharing only five constraints on the order of accumulating mutations with mutation at residue 54 as the parental event. The results also suggest that mutations can accumulate along various alternative paths within subtypes, as opposed to a unique total temporal ordering. Beyond HIV drug resistance, the statistical methodology is applicable more generally for the comparison of inferred mutational networks between any two groups.     

Frequent Intra-Subtype Recombination among HIV-1 Circulating in Tanzania

The study estimated the prevalence of HIV-1 intra-subtype recombinant variants among female bar and hotel workers in Tanzania. While intra-subtype recombination occurs in HIV-1, it is generally underestimated. HIV-1 env gp120 V1-C5 quasispecies from 45 subjects were generated by single-genome amplification and sequencing (median (IQR) of 38 (28–50) sequences per subject). Recombination analysis was performed using seven methods implemented within the recombination detection program version 3, RDP3. HIV-1 sequences were considered recombinant if recombination signals were detected by at least three methods with p-values of ≤0.05 after Bonferroni correction for multiple comparisons. HIV-1 in 38 (84%) subjects showed evidence for intra-subtype recombination including 22 with HIV-1 subtype A1, 13 with HIV-1 subtype C, and 3 with HIV-1 subtype D. The distribution of intra-patient recombination breakpoints suggested ongoing recombination and showed selective enrichment of recombinant variants in 23 (60%) subjects. The number of subjects with evidence of intra-subtype recombination increased from 29 (69%) to 36 (82%) over one year of follow-up, although the increase did not reach statistical significance. Adjustment for intra-subtype recombination is important for the analysis of multiplicity of HIV infection. This is the first report of high prevalence of intra-subtype recombination in the HIV/AIDS epidemic in Tanzania, a region where multiple HIV-1 subtypes co-circulate. HIV-1 intra-subtype recombination increases viral diversity and presents additional challenges for HIV-1 vaccine design.     

Classification of HIV-1 sequences using profile Hidden Markov Models

Accurate classification of HIV-1 subtypes is essential for studying the dynamic spatial distribution pattern of HIV-1 subtypes and also for developing effective methods of treatment that can be targeted to attack specific subtypes. We propose a classification method based on profile Hidden Markov Model that can accurately identify an unknown strain. We show that a standard method that relies on the construction of a positive training set only, to capture unique features associated with a particular subtype, can accurately classify sequences belonging to all subtypes except B and D. We point out the drawbacks of the standard method; namely, an arbitrary choice of threshold to distinguish between true positives and true negatives, and the inability to discriminate between closely related subtypes. We then propose an improved classification method based on construction of a positive as well as a negative training set to improve discriminating ability between closely related subtypes like B and D. Finally, we show how the improved method can be used to accurately determine the subtype composition of Common Recombinant Forms of the virus that are made up of two or more subtypes. Our method provides a simple and highly accurate alternative to other classification methods and will be useful in accurately annotating newly sequenced HIV-1 strains.     

Recombination confounds the early evolutionary history of human immunodeficiency virus type 1: subtype G is a circulating recombinant form

Human immunodeficiency virus type 1 (HIV-1) is classified in nine subtypes (A to D, F, G, H, J, and K), a number of subsubtypes, and several circulating recombinant forms (CRFs). Due to the high level of genetic diversity within HIV-1 and to its worldwide distribution, this classification system is widely used in fields as diverse as vaccine development, evolution, epidemiology, viral fitness, and drug resistance. Here, we demonstrate how the high recombination rates of HIV-1 may confound the study of its evolutionary history and classification. Our data show that subtype G, currently classified as a pure subtype, has in fact a recombinant history, having evolved following recombination between subtypes A and J and a putative subtype G parent. In addition, we find no evidence for recombination within one of the lineages currently classified as a CRF, CRF02_AG. Our analysis indicates that CRF02_AG was the parent of the recombinant subtype G, rather than the two having the opposite evolutionary relationship, as is currently proposed. Our results imply that the current classification of HIV-1 subtypes and CRFs is an artifact of sampling history, rather than reflecting the evolutionary history of the virus. We suggest a reanalysis of all pure subtypes and CRFs in order to better understand how high rates of recombination have influenced HIV-1 evolutionary history.     

HIV-1 genetic diversity in Western Brittany,France

We describe human immunodeficiency virus type 1 (HIV-1) diversity in Western Brittany, France, and trace the dissemination of HIV-1 non-B subtype infection. The strategy for HIV-1 subtyping used involved subtype specific enzyme immunoassays, heteroduplex mobility assays and phylogenetic analysis of the sequences of env encoding the V3 loop region. Samples were obtained from 567 patients: 465 (82%) were of subtype B and 66 (11.6%) were not (20 were subtype A, 11 subtype C, four subtype D, seven subtype F, five subtype G and 19 others with circulating recombinant forms: 4CRF01_AE, 11CRF02_AG, 1H, 3CRF11_cpx). These findings are consistent with other studies of French populations. There is an epidemiological correlation between subtype B and homosexual or heterosexual contamination in France and between non-B subtype and heterosexual contamination in Africa.     

Similar Replicative Fitness Is Shared by the Subtype B and Unique BF Recombinant HIV-1 Isolates that Dominate the Epidemic in Argentina

The HIV-1 epidemic in South America is dominated by pure subtypes (mostly B and C) and more than 7 BF and BC recombinant forms. In Argentina, circulating recombinant forms (CRFs) comprised of subtypes B and F make up more than 50% of HIV infections. For this study, 28 HIV-1 primary isolates were obtained from patients in Buenos Aires, Argentina and initially classified into subtype B (n = 9, 32.1%), C (n = 1, 3.6%), and CRFs (n = 18, 64.3%) using partial pol and vpu-env sequences, which proved to be inconsistent and inaccurate for these phylogenetic analyses. Near full length genome sequences of these primary HIV-1 isolates revealed that nearly all intersubtype BF recombination sites were unique and countered previous “CRF” B/F classifications. The majority of these Argentinean HIV-1 isolates were CCR5-using but 4 had a dual/mixed tropism as predicted by both phenotypic and genotypic assays. Comparison of the replicative fitness of these BF primary HIV-1 isolates to circulating B, F, and C HIV-1 using pairwise competitions in peripheral blood mononuclear cells (PBMCs) indicated a similarity in fitness of these BF recombinants to subtypes B and F HIV-1 (of the same co-receptor usage) whereas subtype C HIV-1 was significantly less fit than all as previously reported. These results suggest that the multitude of BF HIV-1 strains present within the Argentinean population do not appear to have gained replicative fitness following recent B and F recombination events.     

Subtypes of human immunodeficiency virus type 1 and disease stage among women in Nairobi, Kenya

In sub-Saharan Africa, where the effects of human immunodeficiency virus type 1 (HIV-1) have been most devastating, there are multiple subtypes of this virus. The distribution of different subtypes within African populations is generally not linked to particular risk behaviors. Thus, Africa is an ideal setting in which to examine the diversity and mixing of viruses from different subtypes on a population basis. In this setting, it is also possible to address whether infection with a particular subtype is associated with differences in disease stage. To address these questions, we analyzed the HIV-1 subtype, plasma viral loads, and CD4 lymphocyte levels in 320 women from Nairobi, Kenya. Subtype was determined by a combination of heteroduplex mobility assays and sequence analyses of envelope genes, using geographically diverse subtype reference sequences as well as envelope sequences of known subtype from Kenya. The distribution of subtypes in this population was as follows: subtype A, 225 (70.3%); subtype D, 65 (20.5%); subtype C, 22 (6.9%); and subtype G, 1 (0.3%). Intersubtype recombinant envelope genes were detected in 2.2% of the sequences analyzed. Given that the sequences analyzed represented only a small fraction of the proviral genome, this suggests that intersubtype recombinant viral genomes may be very common in Kenya and in other parts of Africa where there are multiple subtypes. The plasma viral RNA levels were highest in women infected with subtype C virus, and women infected with subtype C virus had significantly lower CD4 lymphocyte levels than women infected with the other subtypes. Together, these data suggest that women in Kenya who are infected with subtype C viruses are at more advanced stages of immunosuppression than women infected with subtype A or D. There are at least two models to explain the data from this cross-sectional study; one is that infection with subtype C is associated with a more rapid disease progression, and the second is that subtype C represents an older epidemic in Kenya. Discriminating between these possibilities in a longitudinal study will be important for increasing our understanding of the role of specific subtypes in the transmission and pathogenesis of HIV-1.     

HIV-1 nef protein structures associated with brain infection and dementia pathogenesis

The difference between regional rates of HIV-associated dementia (HAD) in patients infected with different subtypes of HIV suggests that genetic determinants exist within HIV that influence the ability of the virus to replicate in the central nervous system (in Uganda, Africa, subtype D HAD rate is 89%, while subtype A HAD rate is 24%). HIV-1 nef is a multifunctional protein with known toxic effects in the brain compartment. The goal of the current study was to identify if specific three-dimensional nef structures may be linked to patients who developed HAD. HIV-1 nef structures were computationally derived for consensus brain and non-brain sequences from a panel of patients infected with subtype B who died due to varied disease pathologies and consensus subtype A and subtype D sequences from Uganda. Site directed mutation analysis identified signatures in brain structures that appear to change binding potentials and could affect folding conformations of brain-associated structures. Despite the large sequence variation between HIV subtypes, structural alignments confirmed that viral structures derived from patients with HAD were more similar to subtype D structures than to structures derived from patient sequences without HAD. Furthermore, structures derived from brain sequences of patients with HAD were more similar to subtype D structures than they were to their own non-brain structures. The potential finding of a brain-specific nef structure indicates that HAD may result from genetic alterations that alter the folding or binding potential of the protein.     

Genetic analysis of HIV-1 subtypes in Nairobi, Kenya

Background

Genetic analysis of a viral infection helps in following its spread in a given population, in tracking the routes of infection and, where applicable, in vaccine design. Additionally, sequence analysis of the viral genome provides information about patterns of genetic divergence that may have occurred during viral evolution.

Objective

In this study we have analyzed the subtypes of Human Immunodeficiency Virus -1 (HIV-1) circulating in a diverse sample population of Nairobi, Kenya.

Methodology

69 blood samples were collected from a diverse subject population attending the Aga Khan University Hospital in Nairobi, Kenya. Total DNA was extracted from peripheral blood mononuclear cells (PBMCs), and used in a Polymerase Chain Reaction (PCR) to amplify the HIV gag gene. The PCR amplimers were partially sequenced, and alignment and phylogenetic analysis of these sequences was performed using the Los Alamos HIV Database.

Results

Blood samples from 69 HIV-1 infected subjects from varying ethnic backgrounds were analyzed. Sequence alignment and phylogenetic analysis showed 39 isolates to be subtype A, 13 subtype D, 7 subtype C, 3 subtype AD and CRF01_AE, 2 subtype G and 1 subtype AC and 1 AG. Deeper phylogenetic analysis revealed HIV subtype A sequences to be highly divergent as compared to subtypes D and C.

Conclusion

Our analysis indicates that HIV-1 subtypes in the Nairobi province of Kenya are dominated by a genetically diverse clade A. Additionally, the prevalence of highly divergent, complex subtypes, intersubtypes, and the recombinant forms indicates viral mixing in Kenyan population, possibly as a result of dual infections.     

The distribution of HIV-1 subtypes in Fukuoka, Japan.

To determine the incidence of human immunodeficiency virus type-1 (HIV-1) subtypes in Fukuoka, Japan, viruses from 41 HIV-1 infected individuals were subtyped. Subtyping by V3-loop enzyme-linked immunosorbent assay (ELISA) showed 31 of the 41 subjects as subtype B (MN type), one as subtype A, one as subtype C, and eight untypable. The subject infected with subtype C was identified as a foreigner; the subtype A subject was Japanese. A phylogenetic analysis of nucleic acid sequences from the env C2-V3 region was also conducted. Genetic subtyping was successful for 25 samples: 23 samples were determined as subtype B, one subtype A and one subtype E. One of the individuals infected with subtype B, as well as the subtype A and subtype E subjects, were not Japanese. This study indicated that subtype B (USA and European type) is still dominant among HIV-1 infections in Fukuoka. Further, no Japanese were subtype E positive, which is increasing in the Kanto region. It is notable, however, that subtype A and subtype C infections, which are rare in Japan, were found in Fukuoka, located far from the metropolitan area of Tokyo.     

Identification of three feline immunodeficiency virus (FIV) env gene subtypes and comparison of the FIV and human immunodeficiency virus type 1 evolutionary patterns.

Feline immunodeficiency virus (FIV) is a lentivirus associated with AIDS-like illnesses in cats. As such, FIV appears to be a feline analog of human immunodeficiency virus (HIV). A hallmark of HIV infection is the large degree of viral genetic diversity that can develop within an infected individual and the even greater and continually increasing level of diversity among virus isolates from different individuals. Our goal in this study was to determine patterns of FIV genetic diversity by focusing on a 684-nucleotide region encompassing variable regions V3, V4, and V5 of the FIV env gene in order to establish parallels and distinctions between FIV and HIV type 1 (HIV-1). Our data demonstrate that, like HIV-1, FIV can be separated into distinct envelope sequence subtypes (three are described here). Similar to that found for HIV-1, the pairwise sequence divergence within an FIV subtype ranged from 2.5 to 15.0%, whereas that between subtypes ranged from 17.8 to 26.2%. However, the high number of synonymous nucleotide changes among FIV V3 to V5 env sequences may also include a significant number of back mutations and suggests that the evolutionary distances among FIV subtypes are underestimated. Although only a few subtype B viruses were available for examination, the pattern of diversity between the FIV A and B subtypes was found to be significantly distinct; subtype B sequences had proportionally fewer mutations that changed amino acids, compared with silent changes, suggesting a more advanced state of adaptation to the host. No similar distinction was evident for HIV-1 subtypes. The diversity of FIV genomes within individual infected cats was found to be as high as 3.7% yet twofold lower than that within HIV-1-infected people over a comparable region of the env gene. Despite these differences, significant parallels between patterns of FIV evolution and HIV-1 evolution exist, indicating that a wide array of potentially divergent virus challenges need to be considered in FIV vaccine and pathogenesis studies.     

Characterization update of HIV-1 M subtypes diversity and proposal for subtypes A and D sub-subtypes reclassification

Background

The large and constantly evolving HIV-1 pandemic has led to an increasingly complex diversity. Because of some taxonomic difficulties among the most diverse HIV-1 subtypes, and taking advantage of the large amount of sequence data generated in the recent years, we investigated novel lineage patterns among the main HIV-1 subtypes.

Results

All HIV full-length genomes available in public databases were analysed (n?=?2017). Maximum likelihood phylogenies and pairwise genetic distance were obtained. Clustering patterns and mean distributions of genetic distances were compared within and across the current groups, subtypes and sub-subtypes of HIV-1 to detect and analyse any divergent lineages within previously defined HIV lineages. The level of genetic similarity observed between most HIV clades was deeply consistent with the current classification. However, both subtypes A and D showed evidence of further intra-subtype diversification not fully described by the nomenclature system at the time and could be divided into several distinct sub-subtypes.

Conclusions

With this work, we propose an updated nomenclature of sub-types A and D better reflecting their current genetic diversity and evolutionary patterns. Allowing a more accurate nomenclature and classification system is a necessary step for easier subtyping of HIV strains and a better detection or follow-up of viral epidemiology shifts.

     

Estimating the Rate of Intersubtype Recombination in Early HIV-1 Group M Strains

West Central Africa has been implicated as the epicenter of the HIV-1 epidemic, and almost all group M subtypes can be found there. Previous analysis of early HIV-1 group M sequences from Kinshasa in the Democratic Republic of Congo, formerly Zaire, revealed that isolates from a number of individuals fall in different positions in phylogenetic trees constructed from sequences from opposite ends of the genome as a result of recombination between viruses of different subtypes. Here, we use discrete ancestral trait mapping to develop a procedure for quantifying HIV-1 group M intersubtype recombination across phylogenies, using individuals'' gag (p17) and env (gp41) subtypes. The method was applied to previously described HIV-1 group M sequences from samples obtained in Kinshasa early in the global radiation of HIV. Nine different p17 and gp41 intersubtype recombinant combinations were present in the data set. The mean number of excess ancestral subtype transitions (NEST) required to map individuals'' p17 subtypes onto the gp14 phylogeny samples, compared to the number required to map them onto the p17 phylogenies, and vice versa, indicated that excess subtype transitions occurred at a rate of approximately 7 × 10⁻³ to 8 × 10⁻³ per lineage per year as a result of intersubtype recombination. Our results imply that intersubtype recombination may have occurred in approximately 20% of lineages evolving over a period of 30 years and confirm intersubtype recombination as a substantial force in generating HIV-1 group M diversity.