Inference of global HIV-1 sequence patterns and preliminary feature analysis

The epidemiology of HIV-1 varies in different areas of the world, and it is possible that this complexity may leave unique footprints in the viral genome. Thus, we attempted to find significant patterns in global HIV-1 genome sequences. By applying the rule inference algorithm RIPPER (Repeated Incremental Pruning to Produce Error Reduction) to multiple sequence alignments of Env sequences from four classes of compiled datasets, we generated four sets of signature patterns. We found that these patterns were able to distinguish southeastern Asian from nonsoutheastern Asian sequences with 97.5% accuracy, Chinese from non-Chinese sequences with 98.3% accuracy, African from non-African sequences with 88.4% accuracy, and southern African from non-southern African sequences with 91.2% accuracy. These patterns showed different associations with subtypes and with amino acid positions. In addition, some signature patterns were characteristic of the geographic area from which the sample was taken. Amino acid features corresponding to the phylogenetic clustering of HIV-1 sequences were consistent with some of the deduced patterns. Using a combination of patterns inferred from subtypes B, C, and all subtypes chimeric with CRF01_AE worldwide, we found that signature patterns of subtype C were extremely common in some sampled countries (for example, Zambia in southern Africa), which may hint at the origin of this HIV-1 subtype and the need to pay special attention to this area of Africa. Signature patterns of subtype B sequences were associated with different countries. Even more, there are distinct patterns at single position 21 with glycine, leucine and isoleucine corresponding to subtype C, B and all possible recombination forms chimeric with CRF01_AE, which also indicate distinct geographic features. Our method widens the scope of inference of signature from geographic, genetic, and genomic viewpoints. These findings may provide a valuable reference for epidemiological research or vaccine design.     

Variability at human immunodeficiency virus type 1 subtype C protease cleavage sites: an indication of viral fitness?

Naturally occurring polymorphisms in the protease of human immunodeficiency virus type 1 (HIV-1) subtype C would be expected to lead to adaptive (compensatory) changes in protease cleavage sites. To test this hypothesis, we examined the prevalences and patterns of cleavage site polymorphisms in the Gag, Gag-Pol, and Nef cleavage sites of C compared to those in non-C subtypes. Codon-based maximum-likelihood methods were used to assess the natural selection and evolutionary history of individual cleavage sites. Seven cleavage sites (p17/p24, p24/p2, NC/p1, NC/TFP, PR/RT, RT/p66, and p66/IN) were well conserved over time and in all HIV-1 subtypes. One site (p1/p6(gag)) exhibited moderate variation, and four sites (p2/NC, TFP/p6(pol), p6(pol)/PR, and Nef) were highly variable, both within and between subtypes. Three of the variable sites are known to be major determinants of polyprotein processing and virion production. P2/NC controls the rate and order of cleavage, p6(gag) is an important phosphoprotein required for virion release, and TFP/p6(pol), a novel cleavage site in the transframe domain, influences the specificity of Gag-Pol processing and the activation of protease. Overall, 58.3% of the 12 HIV-1 cleavage sites were significantly more diverse in C than in B viruses. When analyzed as a single concatenated fragment of 360 bp, 96.0% of group M cleavage site sequences fell into subtype-specific phylogenetic clusters, suggesting that they coevolved with the virus. Natural variation at C cleavage sites may play an important role, not only in regulation of the viral cycle but also in disease progression and response to therapy.     

Comparative magnitude of cross-strain conservation of HIV variable loop neutralization epitopes

Although the sequence variable loops of the human immunodeficiency virus' (HIV-1) surface envelope glycoprotein (gp120) can exhibit good immunogenicity, characterizing conserved (invariant) cross-strain neutralization epitopes within these loops has proven difficult. We recently developed a method to derive sensitive and specific signature motifs for the three-dimensional (3D) shapes of the HIV-1 neutralization epitopes in the third variable (V3) loop of gp120 that are recognized by human monoclonal antibodies (mAbs). We used the signature motif method to estimate the conservation of these epitopes across circulating worldwide HIV-1 strains. The epitope targeted by the anti-V3 loop neutralizing mAb 3074 is present in 87% of circulating strains, distributed nearly evenly among all subtypes. The results for other anti-V3 Abs are: 3791, present in 63% of primarily non-B subtypes; 2219, present in 56% of strains across all subtypes; 2557, present in 52% across all subtypes; 447-52D, present in 11% of primarily subtype B strains; 537-10D, present in 9% of primarily subtype B strains; and 268-D, present in 5% of primarily subtype B strains. The estimates correlate with in vitro tests of these mAbs against diverse viral panels. The mAb 3074 thus targets an epitope that is nearly completely conserved among circulating HIV-1 strains, demonstrating the presence of an invariant structure hidden in the dynamic and sequence-variable V3 loop in gp120. Since some variable loop regions are naturally immunogenic, designing immunogens to mimic their conserved epitopes may be a promising vaccine discovery approach. Our results suggest one way to quantify and compare the magnitude of the conservation.     

An Evolutionary Model-Based Algorithm for Accurate Phylogenetic Breakpoint Mapping and Subtype Prediction in HIV-1

Genetically diverse pathogens (such as Human Immunodeficiency virus type 1, HIV-1) are frequently stratified into phylogenetically or immunologically defined subtypes for classification purposes. Computational identification of such subtypes is helpful in surveillance, epidemiological analysis and detection of novel variants, e.g., circulating recombinant forms in HIV-1. A number of conceptually and technically different techniques have been proposed for determining the subtype of a query sequence, but there is not a universally optimal approach. We present a model-based phylogenetic method for automatically subtyping an HIV-1 (or other viral or bacterial) sequence, mapping the location of breakpoints and assigning parental sequences in recombinant strains as well as computing confidence levels for the inferred quantities. Our Subtype Classification Using Evolutionary ALgorithms (SCUEAL) procedure is shown to perform very well in a variety of simulation scenarios, runs in parallel when multiple sequences are being screened, and matches or exceeds the performance of existing approaches on typical empirical cases. We applied SCUEAL to all available polymerase (pol) sequences from two large databases, the Stanford Drug Resistance database and the UK HIV Drug Resistance Database. Comparing with subtypes which had previously been assigned revealed that a minor but substantial (≈5%) fraction of pure subtype sequences may in fact be within- or inter-subtype recombinants. A free implementation of SCUEAL is provided as a module for the HyPhy package and the Datamonkey web server. Our method is especially useful when an accurate automatic classification of an unknown strain is desired, and is positioned to complement and extend faster but less accurate methods. Given the increasingly frequent use of HIV subtype information in studies focusing on the effect of subtype on treatment, clinical outcome, pathogenicity and vaccine design, the importance of accurate, robust and extensible subtyping procedures is clear.     

Mapping protease inhibitor resistance to human immunodeficiency virus type 1 sequence polymorphisms within patients

Resistance genotyping provides an important resource for the clinical management of patients infected with human immunodeficiency virus type 1 (HIV-1). However, resistance to protease (PR) inhibitors (PIs) is a complex phenotype shaped by interactions among nearly half of the residues in HIV-1 PR. Previous studies of the genetic basis of PI resistance focused on fixed substitutions among populations of HIV-1, i.e., host-specific adaptations. Consequently, they are susceptible to a high false discovery rate due to founder effects. Here, we employ sequencing “mixtures” (i.e., ambiguous base calls) as a site-specific marker of genetic variation within patients that is independent of the phylogeny. We demonstrate that the transient response to selection by PIs is manifested as an excess of nonsynonymous mixtures. Using a sample of 5,651 PR sequences isolated from both PI-naive and -treated patients, we analyze the joint distribution of mixtures and eight PIs as a Bayesian network, which distinguishes residue-residue interactions from direct associations with PIs. We find that selection for resistance is associated with the emergence of nonsynonymous mixtures in two distinct groups of codon sites clustered along the substrate cleft and distal regions of PR, respectively. Within-patient evolution at several positions is independent of PIs, including those formerly postulated to be involved in resistance. These positions are under strong positive selection in the PI-naive patient population, implying that other factors can produce spurious associations with resistance, e.g., mutational escape from the immune response.     

Genetic diversity of human immunodeficiency virus type 2: evidence for distinct sequence subtypes with differences in virus biology.

The virulence properties of human immunodeficiency virus type 2 (HIV-2) are known to vary significantly and to range from relative attenuation in certain individuals to high-level pathogenicity in others. These differences in clinical manifestations may, at least in part, be determined by genetic differences among infecting virus strains. Evaluation of the full spectrum of HIV-2 genetic diversity is thus a necessary first step towards understanding its molecular epidemiology, natural history of infection, and biological diversity. In this study, we have used nested PCR techniques to amplify viral sequences from the DNA of uncultured peripheral blood mononuclear cells from 12 patients with HIV-2 seroreactivity. Sequence analysis of four nonoverlapping genomic regions allowed a comprehensive analysis of HIV-2 phylogeny. The results revealed (i) the existence of five distinct and roughly equidistant evolutionary lineages of HIV-2 which, by analogy with HIV-1, have been termed sequence subtypes A to E; (ii) evidence for a mosaic HIV-2 genome, indicating that coinfection with genetically divergent strains and recombination can occur in HIV-2-infected individuals; and (iii) evidence supporting the conclusion that some of the HIV-2 subtypes may have arisen from independent introductions of genetically diverse sooty mangabey viruses into the human population. Importantly, only a subset of HIV-2 strains replicated in culture: all subtype A viruses grew to high titers, but attempts to isolate representatives of subtypes C, D, and E, as well as the majority of subtype B viruses, remained unsuccessful. Infection with all five viral subtypes was detectable by commercially available serological (Western immunoblot) assays, despite intersubtype sequence differences of up to 25% in the gag, pol, and env regions. These results indicate that the genetic and biological diversity of HIV-2 is far greater than previously appreciated and suggest that there may be subtype-specific differences in virus biology. Systematic natural history studies are needed to determine whether this heterogeneity has clinical relevance and whether the various HIV-2 subtypes differ in their in vivo pathogenicity.     

Evolution of the human immunodeficiency virus type 1 subtype-specific V3 domain is confined to a sequence space with a fixed distance to the subtype consensus.

Human immunodeficiency virus type 1 (HIV-1) strains can be separated into genetic subtypes based on phylogenetic analysis of the envelope gene. Once it had been shown that population-wide intrasubtype genetic variation of HIV-1 strains increases in the course of the AIDS epidemic, it remained uncertain whether HIV-1 subtypes are phenotypic entities spreading as distinct virus populations. To examine this, we applied Eigen's concepts of sequence geometry and fitness topography to the analysis of intrasubtype evolution of the gp120 V3 domain of HIV-1 subtypes A, B, C, and D in the course of the global AIDS epidemic. We observed that despite the high evolution rate of HIV-1, the nonsynonymous distances to the subtype consensus of sequences obtained early in the epidemic are similar to those obtained more than 10 years later, in contrast to the synonymous distances, which increased steadily over time. For HIV-1 subtype B, we observed that the evolution rate of the individual sequences is independent of their distance from the subtype B consensus, but for the individual sequences most distant from the consensus evolution away from the consensus is constrained. As a result, individual HIV-1 genomes fluctuate within a sequence space with fixed distance to the subtype consensus. Our findings suggest that the evolution of the V3 domain of HIV-1 subtypes A, B, C, and D is confined to an area in sequence space within a fixed distance to the consensus of a respective subtype. This in turn indicates that each HIV-1 subtype is a distinct viral quasispecies that is well adapted to the present environment, able to maintain its identity in the V3 region over time, and unlikely to merge during progression of the AIDS epidemic.     

Evidence for heterogeneous selective pressures in the evolution of the env gene in different human immunodeficiency virus type 1 subtypes

Recent studies have demonstrated the emergence of human immunodeficiency virus type 1 (HIV-1) subtypes with various levels of fitness. Using heterogeneous maximum-likelihood models of adaptive evolution implemented in the PAML software package, with env sequences representing each HIV-1 group M subtype, we examined the various intersubtype selective pressures operating across the env gene. We found heterogeneity of evolutionary mechanisms between the different subtypes with a category of amino acid sites observed that had undergone positive selection for subtypes C, F1, and G, while these sites had undergone purifying selection in all other subtypes. Also, amino acid sites within subtypes A and K that had undergone purifying selection were observed, while these sites had undergone positive selection in all other subtypes. The presence of such sites indicates heterogeneity of selective pressures within HIV-1 group M subtype evolution that may account for the various levels of fitness of the subtypes.     

HIV-1 Genetic Diversity and Transmitted Drug Resistance Mutations among Patients from the North, Central and South Regions of Angola

Background

Angola presents a very complex HIV-1 epidemic characterized by the co-circulation of several HIV-1 group M subtypes, intersubtype recombinants and unclassified (U) variants. The viral diversity outside the major metropolitan regions (Luanda and Cabinda) and the prevalence of transmitted drug resistance mutations (DRM) since the introduction of HAART in 2004, however, has been barely studied.

Methods

One hundred and one individuals from the Central (n = 44), North (n = 35), and South (n = 22) regions of Angola were diagnosed as HIV-1 positive and had their blood collected between 2008 and 2010, at one of the National Referral Centers for HIV diagnosis, the Kifangondo Medical Center, located in the border between the Luanda and Bengo provinces. Angolan samples were genotyped based on phylogenetic and bootscanning analyses of the pol (PR/RT) gene and their drug resistance profile was analyzed.

Results

Among the 101 samples analyzed, 51% clustered within a pure group M subtype, 42% were classified as intersubtype recombinants, and 7% were denoted as U. We observed an important variation in the prevalence of different HIV-1 genetic variants among country regions, with high frequency of subtype F1 in the North (20%), intersubtype recombinants in the Central (42%), and subtype C in the South (45%). Statistically significant difference in HIV-1 clade distribution was only observed in subtype C prevalence between North vs South (p = 0.0005) and Central vs South (p = 0.0012) regions. DRM to NRTI and/or NNRTI were detected in 16.3% of patients analyzed.

Conclusions

These results demonstrate a heterogeneous distribution of HIV-1 genetic variants across different regions in Angola and also revealed an unexpected high frequency of DRM to RT inhibitors in patients that have reported no antiretroviral usage, which may decrease the efficiency of the standard first-line antiretroviral regimens currently used in the country.     

Comparative study of adaptive molecular evolution in different human immunodeficiency virus groups and subtypes

Molecular adaptation, as characterized by the detection of positive selection, was quantified in a number of genes from different human immunodeficiency virus type 1 (HIV-1) group M subtypes, group O, and an HIV-2 subtype using the codon-based maximum-likelihood method of Yang and coworkers (Z. H. Yang, R. Nielsen, N. Goldman, and A. M. K. Pedersen, Genetics 155:431-449, 2000). The env gene was investigated further since it exhibited the strongest signal for positive selection compared to those of the other two major HIV genes (gag and pol). In order to investigate the pattern of adaptive evolution across env, the location and strength of positive selection in different HIV-1 sequence alignments was compared. The number of sites having a significant probability of being positively selected varied among these different alignment data sets, ranging from 25 in HIV-1 group M subtype A to 40 in HIV-1 group O. Strikingly, there was a significant tendency for positively selected sites to be located at the same position in different HIV-1 alignments, ranging from 10 to 16 shared sites for the group M intersubtype comparisons and from 6 to 8 for the group O to M comparisons, suggesting that all HIV-1 variants are subject to similar selective forces. As the host immune response is believed to be the dominant driving force of adaptive evolution in HIV, this result would suggest that the same sites are contributing to viral persistence in diverse HIV infections. Thus, the positions of the positively selected sites were investigated in reference to the inferred locations of different epitope types (antibody, T helper, and cytotoxic T lymphocytes) and the positions of N and O glycosylation sites. We found a significant tendency for positively selected sites to fall outside T-helper epitopes and for positively selected sites to be strongly associated with N glycosylation sites.     

Selection promotes organ compartmentalization in HIV-1: evidence from gag and pol genes

The existence of organ-specific human immunodeficiency virus type 1 (HIV-1) populations within infected hosts has been long lasting studied. Previous work established that population subdivision by organs occurs at the envelope env gene, but less is known about other genomic regions. Here, we used a population genetics approach to detect organ compartmentalization in proviral sequences of HIV-1 gag and pol genes. Significant population structure was found in pol (100% of cases) and gag (33%) pair-wise organ comparisons. The degree of compartmentalization positively correlated with the ratio of nonsynonymous to synonymous substitutions, and codons showing organ compartmentalization were more likely to be under significantly positive selection. This suggests that HIV-1 populations dynamically adapt to locally variable intra-host environments. In the case of pol gene, differential penetration of antiretroviral drugs might account for the observed pattern, whereas for gag gene, local selective pressures remain unexplored.     

Genetic analysis of Indian HIV-1 nef: subtyping, variability and implications

HIV-1 subtype C is predominant in India and globally. In the present study, we analyze HIV-1 subtype C regulatory protein Nef sequences from five recent Indian seroconverters and five long-term survivors (LTSs) for variability at crucial functional domains. Sequence analysis suggested the possibility of using regulatory gene sequences for viral subtyping and evolutionary studies apart from structural genes. In the phylogenetic tree, Indian nef sequences segregated away from other reported subtype C sequences, forming an Indian subclade within subtype C. Our studies also suggested no evidence for the association of truncated Nef with slow progression of disease, as all LTSs had intact Nef. We could identify some variations in juxtapositions to crucial functional domains, especially in seroconverter sequences, when comparing them with others. In phylogenetic analysis, specifically for the base-pair regions 411-428 and 478-525, our seroconverter sequences segregated away from those reported earlier in the literature, indicating specific evolutionary changes in these conserved regions of nef in currently circulating viruses. But the dN/dS ratio for our samples was less than one on comparing them with reported subtype C and representative sequences of different clades, strongly emphasizing the necessity of sequence conservation at different disease stages and even across clades. HLA-I binding epitope predictions for common Indian HLAs indicated that specific mutations in seroconverter Nef may alter the intensity of epitope binding, which may alter the outcome of the immune response. Hence these data would be useful in designing Nef epitopes to be included in multi-epitope HIV-1 vaccine for the Indian population and would also be of immense help in HIV-1 evolutionary studies.