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.     

Complex positive selection pressures drive the evolution of HIV-1 with different co-receptor tropisms

HIV-1 co-receptor tropism is central for understanding the transmission and pathogenesis of HIV-1 infection. We performed a genome-wide comparison between the adaptive evolution of R5 and X4 variants from HIV-1 subtypes B and C. The results showed that R5 and X4 variants experienced differential evolutionary patterns and different HIV-1 genes encountered various positive selection pressures, suggesting that complex selection pressures are driving HIV-1 evolution. Compared with other hypervariable regions of Gp120, significantly more positively selected sites were detected in the V3 region of subtype B X4 variants, V2 region of subtype B R5 variants, and V1 and V4 regions of subtype C X4 variants, indicating an association of positive selection with co-receptor recognition/binding. Intriguingly, a significantly higher proportion (33.3% and 55.6%, P<0.05) of positively selected sites were identified in the C3 region than other conserved regions of Gp120 in all the analyzed HIV-1 variants, indicating that the C3 region might be more important to HIV-1 adaptation than previously thought. Approximately half of the positively selected sites identified in the env gene were identical between R5 and X4 variants. There were three common positively selected sites (96, 113 and 281) identified in Gp41 of all X4 and R5 variants from subtypes B and C. These sites might not only suggest a functional importance in viral survival and adaptation, but also imply a potential cross-immunogenicity between HIV-1 R5 and X4 variants, which has important implications for AIDS vaccine development.     

A study of the coevolutionary patterns operating within the env gene of the HIV-1 group M subtypes

The env gene of human immunodeficiency virus (HIV) is a functionally important gene responsible for the production of protein products (gp120 and gp41) involved in host cell recognition, binding, and entry. This occurs through a complex and, as yet, not fully understood process of protein-protein interaction and within and between protein functional communication. Exposure on the surface of active HIV virions means the gp120-gp41 complexes are subjected to intense immune system pressure and have, therefore, evolved mechanisms to avoid neutralization. Using protein-coding sequences representing all the HIV type-1 (HIV-1) group M subtypes, we have identified amino acids within the env gene whose evolution is inextricably linked over the entire HIV-1 group M epidemic. We identified 848 pairs of coevolving residues (involving 263 out of 764 amino acid sites), which represent 0.29% of all possible pairs. Of the coevolving pairs, 68% were significantly correlated by hydrophobicity, molecular weight, or both hydrophobicity and molecular weight. Subsequent grouping of coevolving pairs resulted in the identification of 290 groups of amino acid residues, with the size of these groups ranging from 2 to 10 amino acid residues. Many of these dependencies are correlated by function including CD4 binding, coreceptor binding, glycosylation, and protein-protein interaction. This analysis provides important information regarding the functional dependencies observed within all the HIV-1 group M subtypes and may assist in the identification of functional protein domains and therapeutic targets within the HIV-1 env gene.     

QUASI analysis of the HIV-1 envelope sequences in the Los Alamos National Laboratory HIV sequence database: pattern and distribution of positive selection sites and their frequencies over years.

The envelope (env) protein of human immunodeficiency virus type 1 (HIV-1) plays a crucial role in virus entry and is a central target for HIV vaccine design. Using the QUASI program, we analyzed the conserved regions of all currently available env sequences in the Los Alamos National Laboratory HIV Sequence Database and identified positive selection (PS) sites that are likely to be restricted by host immune responses. We found that PS sites are dispersed across conserved regions of env sequence, and that the C3, C4, and C5 regions were the most targeted. Several regions were identified as being PS free and were mainly distributed in the C1 and C2 regions. When comparing individual QUASI PS site frequencies across clades and geographical regions with the overall frequency of the entire env database, the env sequences from North America showed significantly lower PS site frequency, while those from Asia were significantly higher using Student's t test. The QUASI PS site frequency of env proteins from viruses isolated from different years showed that the PS site frequencies of the env population increased over time. Our study provides an overview of PS sites across the conserved regions of HIV-1 env sequences.     

Evolution and biological characterization of human immunodeficiency virus type 1 subtype E gp120 V3 sequences following horizontal and vertical virus transmission in a single family

It has been suggested that immune-pressure-mediated positive selection operates to maintain the antigenic polymorphism on the third variable (V3) loop of the gp120 of human immunodeficiency virus type 1 (HIV-1). Here we present evidence, on the basis of sequencing 147 independently cloned env C2/V3 segments from a single family (father, mother, and their child), that the intensity of positive selection is related to the V3 lineage. Phylogenetic analysis and amino acid comparison of env C2/V3 and gag p17/24 regions indicated that a single HIV-1 subtype E source had infected the family. The analyses of unique env C2/V3 clones revealed that two V3 lineage groups had evolved in the parents. Group 1 was maintained with low variation in all three family members regardless of the clinical state or the length of infection, whereas group 2 was only present in symptomatic individuals and was more positively charged and diverse than group 1. Only virus isolates carrying the group 2 V3 sequences infected and induced syncytia in MT2 cells, a transformed CD4(+)-T-cell line. A statistically significant excess of nonsynonymous substitutions versus synonymous substitutions was demonstrated only for the group 2 V3 region. The data suggest that HIV-1 variants, possessing the more homogeneous group 1 V3 element and exhibiting the non-syncytium-inducing phenotype, persist in infected individuals independent of clinical status and appear to be more resistant to positive selection pressure.     

Human Immunodeficiency Virus Type 1 Intergroup (M/O) Recombination in Cameroon

Here we describe, for the first time, recombinants between two highly divergent major groups of human immunodeficiency virus type 1 (HIV-1), M and O, within a Cameroonian woman infected with three different HIV-1 strains, a group O virus, a subtype D virus, and a recently reported IBNG (A/G)-like recombinant virus. Using nested extra-long PCR amplification, we sequenced from the pol region to the env region including accessory genes of the viral genome obtained from the patient's uncultured peripheral blood mononuclear cells and examined the phylogenetic position of each gene. Compared with sequential blood samples obtained in 1995 and 1996, there were multiple segmental exchanges between three HIV-1 strains (O, D, and IBNG) and all the recombinants appeared to be derived from a common M/O ancestor. Importantly, recombination between groups M and O occurred, even though the homology between these two groups is 69, 76, 68, and 55% in the gag, pol, vif-vpr, and env regions, respectively. Recombination between strains with such distant lineages may contribute substantially to generating new HIV-1 variants.     

Genealogical evidence for positive selection in the nef gene of HIV-1.

The pattern and process of evolution in the nef gene of HIV-1 was analyzed within and among patients. Using a maximum likelihood method that allows for variable intensity of selection pressure among codons, strong positive selection was detected in a hemophiliac patient over 30 mo of infection. By reconstructing the process of allele substitution in this patient using parsimony, the synapomorphic amino acid changes separating each time point were found to have high probabilities of being under positive selection, with selective coefficients of at least 3.6%. Positive selection was also detected among 39 nef sequences from HIV-1 subtype B. In contrast, multiple pairwise comparisons of nonsynonymous and synonymous substitution rates provided no good evidence for positive selection and sliding window analyses failed to detect most positively selected sites. These findings demonstrate that positive selection is an important determinant of nef gene evolution and that genealogy-based methods outperform pairwise methods in the detection of adaptive evolution. Mapping the locations of positively selected sites may also be of use in identifying targets of the immune response and hence aid vaccine design.     

Semen-specific genetic characteristics of human immunodeficiency virus type 1 env

Human immunodeficiency virus type 1 (HIV-1) in the male genital tract may comprise virus produced locally in addition to virus transported from the circulation. Virus produced in the male genital tract may be genetically distinct, due to tissue-specific cellular characteristics and immunological pressures. HIV-1 env sequences derived from paired blood and semen samples from the Los Alamos HIV Sequence Database were analyzed to ascertain a male genital tract-specific viral signature. Machine learning algorithms could predict seminal tropism based on env sequences with accuracies exceeding 90%, suggesting that a strong genetic signature does exist for virus replicating in the male genital tract. Additionally, semen-derived viral populations exhibited constrained diversity (P < 0.05), decreased levels of positive selection (P < 0.025), decreased CXCR4 coreceptor utilization, and altered glycosylation patterns. Our analysis suggests that the male genital tract represents a distinct selective environment that contributes to the apparent genetic bottlenecks associated with the sexual transmission of HIV-1.     

Evidence of selection pressures of neuraminidase gene (NA) of influenza A virus subtype H5N1 on different hosts in Guangxi Province of China

In the present study, the possible evidence of positive selection was analyzed for the neuraminidase (NA) sequences of Guangxi H5N1 strains of China. Based on an overall site-specific positive selection analysis, it was found that NA gene of H5N1 Guangxi strains underwent purifying selection and no significant positively selected sites were identified. For the branch-specific positive selection analysis, there was no positive selection evidence for the branches leading to different poultry hosts (chicken, duck and goose). Conclusively, positive selection seems not possible (if not rare) for the NA gene in influenza H5N1 subtype, at least for the samples found in Guangxi Province of China.     

Antigenicity and immunogenicity of a synthetic human immunodeficiency virus type 1 group m consensus envelope glycoprotein

Genetic variation of human immunodeficiency virus (HIV-1) represents a major obstacle for AIDS vaccine development. To decrease the genetic distances between candidate immunogens and field virus strains, we have designed and synthesized an artificial group M consensus env gene (CON6 gene) to be equidistant from contemporary HIV-1 subtypes and recombinants. This novel envelope gene expresses a glycoprotein that binds soluble CD4, utilizes CCR5 but not CXCR4 as a coreceptor, and mediates HIV-1 entry. Key linear, conformational, and glycan-dependent monoclonal antibody epitopes are preserved in CON6, and the glycoprotein is recognized equally well by sera from individuals infected with different HIV-1 subtypes. When used as a DNA vaccine followed by a recombinant vaccinia virus boost in BALB/c mice, CON6 env gp120 and gp140CF elicited gamma interferon-producing T-cell responses that recognized epitopes within overlapping peptide pools from three HIV-1 Env proteins, CON6, MN (subtype B), and Chn19 (subtype C). Sera from guinea pigs immunized with recombinant CON6 Env gp120 and gp140CF glycoproteins weakly neutralized selected HIV-1 primary isolates. Thus, the computer-generated "consensus" env genes are capable of expressing envelope glycoproteins that retain the structural, functional, and immunogenic properties of wild-type HIV-1 envelopes.     

A new subtype of human immunodeficiency virus type 1 (MVP-5180) from Cameroon.

A new subtype (MVP-5180) of human immunodeficiency virus type 1 (HIV-1) was isolated from a Cameroonian AIDS patient. MVP-5180 was grown in several human T-cell lines and the monocytic U937 line. MVP-5180 DNA could not be amplified by nested primer PCR with conventional env primers and could be only very faintly amplified with gag and pol primers. Most German, Ivoirian, and Malawian anti-HIV-1 sera reacted faintly or moderately with Env proteins in an MVP-5180 immunoblot, whereas some Cameroonian sera reacted strongly. Of HIV-1-infected Cameroonians, 8% were identified by serological methods as infected with MVP-5180; 7% were positive when MVP-5180-specific PCR env primers were used. DNA sequence analysis of MVP-5180 showed that its genetic organization was that of HIV-1, with 65% similarity to HIV-1 and 56% similarity to HIV-2 consensus sequences. The env gene of MVP-5180 had similarities to HIV-1 and HIV-2 of 53 and of 49%, respectively. V3 loop analysis identified a crown of Gly-Pro-Met-Arg by using cloned DNA and Gly-Pro-Leu-Arg by using PCR-amplified DNA, neither of which configuration has been described for other HIV strains. In an analysis of relationships, MVP-5180 occupied a position distant to all other HIV-1 strains, including the chimpanzee simian immunodeficiency virus type 1 SIVcpz and the Uganda virus U455, and closer to the HIV-1/HIV-2 divergence node. MVP-5180, together with another Cameroonian isolate, ANT-70, constitutes a group subtype O of the most divergent HIV-1 isolates yet identified. Characterization of MVP-5180 is important for understanding the natural history of the primate immunodeficiency viruses and for the development of vaccines and diagnostics.     

The replicative fitness of primary human immunodeficiency virus type 1 (HIV-1) group M, HIV-1 group O, and HIV-2 isolates

The main (M) group of human immunodeficiency virus type 1 (HIV-1) is responsible for the global AIDS epidemic while HIV-1 group O (outlier) and HIV type 2 are endemic only in west and central Africa. The failure of HIV-2 and especially HIV-1 group O to spread following the initial zoonotic jumps is not well understood. This study was designed to examine the relative replicative capacities between these human lentiviruses. A pairwise competition experiment was performed with peripheral blood mononuclear cells with eight HIV-2 isolates, 6 group O viruses, and 15 group M viruses of subtype A (2 viruses), B (5 viruses), C (4 viruses), D (2 viruses) and CRF01_AE (2 viruses). HIV-1 group M isolates of any subtype were typically 100-fold-more fit than group O or HIV-2 strains when competed in peripheral blood mononuclear cells from various humans. This order in replicative fitness was also observed when virus pairs were added to human dendritic cells and then cocultured with primary, quiescent T cells, which is the model for HIV-1 transmission. These results suggest that reduced replicative and transmission fitness may be contributing to the low prevalence and limited geographical spread of HIV-2 and group O HIV-1 in the human population.     

Positive selection for gains of N-linked glycosylation sites in hemagglutinin during evolution of H3N2 human influenza A virus

The number of N-linked glycosylation sites in the globular head of hemagglutinin (HA) has increased during evolution of H3N2 human influenza A virus. Here natural selection operating on the gains of N-linked glycosylation sites was examined by using the single-site analysis and the single-substitution analysis. In the single-site analysis, positive selection was not inferred at the amino acid sites where the substitutions generating N-linked glycosylation sites were observed, but was detected at antigenic sites. In contrast, in the single-substitution analysis, positive selection was detected for the amino acid substitutions generating N-linked glycosylation sites. The single-site analysis and the single-substitution analysis appeared to be suitable for detecting recurrent and episodic natural selection, respectively. The gains of N-linked glycosylation sites were likely to be positively selected for the function of shielding antigenic sites from immune responses. At the antigenic sites, positive selection appeared to have operated not only on the radical substitution but also on the conservative substitution in terms of the charge of amino acids, suggesting that the antigenic drift is not a by-product of the evolution of receptor binding avidity in HA of human H3N2 virus.     

SIVcpz from a naturally infected Cameroonian chimpanzee: Biological and genetic comparison with HIV-1 N

Thus far, simian immunodeficiency virus from chimpanzees (SIVcpz) genomes have been characterized as Pan troglodytes troglodytes and show a strong relation with human immunodeficiency virus (HIV)-1 N in their env genes. We fully characterized another SIVcpz from P. t. troglodytes . This chimpanzee (Cam5) was, as was also the host of SIVcpz-cam3, wild born in Cameroon, a region where all three groups of HIV-1 (M, N and O) co-occur. In contrast to other SIVcpz, SIVcpz-cam5 was isolated immediately after the rescue of the animal. Our data demonstrate that SIVcpz-cam5, like SIVcpz-cam3, grows easily on human peripheral blood mononuclear cells (PBMCs) and uses CCR5 as a co-receptor similar to HIV-1 N YBF30. Phylogenetic analysis based on the entire env gene shows that SIVcpz-cam5 falls into the same unique subcluster as HIV-1 N YBF30, SIVcpz-cam3 and SIVcpz-US. A phylogenetic relationship was also found with the vif gene of HIV-1 N. This study provides proof that HIV-1 N related viruses circulate in wild P. t. troglodytes .     

SIVcpz from a naturally infected Cameroonian chimpanzee: biological and genetic comparison with HIV-1 N

Thus far, simian immunodeficiency virus from chimpanzees (SIVcpz) genomes have been characterized as Pan troglodytes troglodytes and show a strong relation with human immunodeficiency virus (HIV)-1 N in their env genes. We fully characterized another SIVcpz from P. t. troglodytes. This chimpanzee (Cam5) was, as was also the host of SIVcpz-cam3, wild born in Cameroon, a region where all three groups of HIV-1 (M, N and O) co-occur. In contrast to other SIVcpz, SIVcpz-cam5 was isolated immediately after the rescue of the animal. Our data demonstrate that SIVcpz-cam5, like SIVcpz-cam3, grows easily on human peripheral blood mononuclear cells (PBMCs) and uses CCR5 as a co-receptor similar to HIV-1 N YBF30. Phylogenetic analysis based on the entire env gene shows that SIVcpz-cam5 falls into the same unique subcluster as HIV-1 N YBF30, SIVcpz-cam3 and SIVcpz-US. A phylogenetic relationship was also found with the vif gene of HIV-1 N. This study provides proof that HIV-1 N related viruses circulate in wild P. t. troglodytes.     

Molecular epidemiology of human immunodeficiency virus-1 in the state of Ceará, Northeast,Brazil

We analyzed, by env and gag heteroduplex mobility assay, 149 human immunodeficiency virus (HIV-1) positive samples collected in Ceará during the year 2000. The prevalence of subtype B was 81.2% and the prevalence of subtype F and B/F recombinants were both 2.7%. Eight (5.4%) and 12 (8%) out of 149 samples showed indeterminate results in the env and gag analysis respectively. By FokI restriction fragment length polymorphism, 34% of the subtype B samples were identified as the typical Brazilian subtype B.In the present study, we identified HIV-1 subtype F and B/F in Ceará for the first time. Our results contribute to the understanding of HIV in Brazil, and may prove useful for the development of vaccine candidates.     

Selection for human immunodeficiency virus type 1 envelope glycosylation variants with shorter V1-V2 loop sequences occurs during transmission of certain genetic subtypes and may impact viral RNA levels

Designing an effective human immunodeficiency virus type 1 (HIV-1) vaccine will rely on understanding which variants, from among the myriad of circulating HIV-1 strains, are most commonly transmitted and determining whether such variants have an Achilles heel. Here we show that heterosexually acquired subtype A HIV-1 envelopes have signature sequences that include shorter V1-V2 loop sequences and fewer predicted N-linked glycosylation sites relative to the overall population of circulating variants. In contrast, recently transmitted subtype B variants did not, and this was true for cases where the major risk factor was homosexual contact, as well as for cases where it was heterosexual contact. This suggests that selection during HIV-1 transmission may vary depending on the infecting subtype. There was evidence from 23 subtype A-infected women for whom there was longitudinal data that those who were infected with viruses with fewer potential N-linked glycosylation sites in V1-V2 had lower viral set point levels. Thus, our study also suggests that the extent of glycosylation in the infecting virus could impact disease progression.     

Epidemiology of foot-and-mouth disease in Landhi Dairy Colony, Pakistan, the world largest Buffalo colony

Background

The high diversity of HIV variants driving the global AIDS epidemic has caused many to doubt whether an effective vaccine against the virus is possible. However, by identifying the selective forces that are driving the ongoing diversification of HIV and characterising their genetic consequences, it may be possible to design vaccines that pre-empt some of the virus' more common evasion tactics. One component of such vaccines might be the envelope protein, gp41. Besides being targeted by both the humoral and cellular arms of the immune system this protein mediates fusion between viral and target cell membranes and is likely to be a primary determinant of HIV transmissibility.

Results

Using recombination aware analysis tools we compared site specific signals of selection in gp41 sequences from different HIV-1 M subtypes and circulating recombinant forms and identified twelve sites evolving under positive selection across multiple major HIV-1 lineages. To identify evidence of selection operating during transmission our analysis included two matched datasets sampled from patients with acute or chronic subtype C infections. We identified six gp41 sites apparently evolving under different selection pressures during acute and chronic HIV-1 infections. These sites mostly fell within functional gp41 domains, with one site located within the epitope recognised by the broadly neutralizing antibody, 4E10.

Conclusion

Whereas these six sites are potentially determinants of fitness and are therefore good candidate targets for subtype-C specific vaccines, the twelve sites evolving under diversifying selection across multiple subtypes might make good candidate targets for broadly protective vaccines.     

Codon-substitution models for heterogeneous selection pressure at amino acid sites

Comparison of relative fixation rates of synonymous (silent) and nonsynonymous (amino acid-altering) mutations provides a means for understanding the mechanisms of molecular sequence evolution. The nonsynonymous/synonymous rate ratio (omega = d(N)d(S)) is an important indicator of selective pressure at the protein level, with omega = 1 meaning neutral mutations, omega < 1 purifying selection, and omega > 1 diversifying positive selection. Amino acid sites in a protein are expected to be under different selective pressures and have different underlying omega ratios. We develop models that account for heterogeneous omega ratios among amino acid sites and apply them to phylogenetic analyses of protein-coding DNA sequences. These models are useful for testing for adaptive molecular evolution and identifying amino acid sites under diversifying selection. Ten data sets of genes from nuclear, mitochondrial, and viral genomes are analyzed to estimate the distributions of omega among sites. In all data sets analyzed, the selective pressure indicated by the omega ratio is found to be highly heterogeneous among sites. Previously unsuspected Darwinian selection is detected in several genes in which the average omega ratio across sites is <1, but in which some sites are clearly under diversifying selection with omega > 1. Genes undergoing positive selection include the beta-globin gene from vertebrates, mitochondrial protein-coding genes from hominoids, the hemagglutinin (HA) gene from human influenza virus A, and HIV-1 env, vif, and pol genes. Tests for the presence of positively selected sites and their subsequent identification appear quite robust to the specific distributional form assumed for omega and can be achieved using any of several models we implement. However, we encountered difficulties in estimating the precise distribution of omega among sites from real data sets.