Genetic association of the antiviral restriction factor TRIM5alpha with human immunodeficiency virus type 1 infection

The innate antiviral factor TRIM5alpha restricts the replication of some retroviruses through its interaction with the viral capsid protein, leading to abortive infection. While overexpression of human TRIM5alpha results in modest restriction of human immunodeficiency virus type 1 (HIV-1), this inhibition is insufficient to block productive infection of human cells. We hypothesized that polymorphisms within TRIM5 may result in increased restriction of HIV-1 infection. We sequenced the TRIM5 gene (excluding exon 5) and the 4.8-kb 5' putative regulatory region in genomic DNA from 110 HIV-1-infected subjects and 96 exposed seronegative persons, along with targeted gene sequencing in a further 30 HIV-1-infected individuals. Forty-eight single nucleotide polymorphisms (SNPs), including 20 with allele frequencies of >1.0%, were identified. Among these were two synonymous and eight nonsynonymous coding polymorphisms. We observed no association between TRIM5 polymorphism in HIV-1-infected subjects and their set-point viral load after acute infection, although one TRIM5 haplotype was weakly associated with more rapid CD4(+) T-cell loss. Importantly, a TRIM5 haplotype containing the nonsynonymous SNP R136Q showed increased frequency among HIV-1-infected subjects relative to exposed seronegative persons, with an odds ratio of 5.49 (95% confidence interval = 1.83 to 16.45; P = 0.002). Nonetheless, we observed no effect of individual TRIM5alpha nonsynonymous mutations on the in vitro HIV-1 susceptibility of CD4(+) T cells. Therefore, any effect of TRIM5alpha polymorphism on HIV-1 infection in primary lymphocytes may depend on combinations of SNPs or on DNA sequences in linkage disequilibrium with the TRIM5alpha coding sequence.     

Host-specific driving force in human immunodeficiency virus type 1 evolution in vivo.

To investigate the process of human immunodeficiency virus type 1 (HIV-1) evolution in vivo, a total of 179 HIV-1 V3 sequences derived from cell-free plasma were determined from serial samples in three epidemiologically linked individuals (one infected blood donor and two transfusion recipients) over a maximum period of 8 years. A systematic analysis of pairwise comparisons of intrapatient sequences, both within and between each sample time point, revealed a preponderance and accumulation of nonsynonymous rather than synonymous substitutions in the V3 loop and flanking regions as they diverged over time. This strongly argues for the dominant role that positive selection for amino acid change plays in governing the pattern and process of HIV-1 env V3 evolution in vivo and nullifies hypotheses of purely neutral or mutation-driven evolution or completely chance events. In addition, different rates of evolution of HIV-1 were observed in these three different individuals infected with the same viral strain, suggesting that the degree of positive pressure for HIV-1 amino acid change is host dependent. Finally, the observed similar rate of accumulation in divergence within and between infected individuals suggests that the process of genetic divergence in the HIV epidemic proceeds regardless of host-to-host transmission events, i.e., that transmission does not reset the evolutionary clock.     

Molecular epidemiology of human immunodeficiency virus type 1 sub-subtype A3 in Senegal from 1988 to 2001

The global human immunodeficiency virus (HIV)epidemic is characterized by significant genetic diversity in circulating viruses. We have recently characterized a group of viruses that form a distinct sub-subtype within the subtype A radiation, which we have designated HIV type 1 (HIV-1) sub-subtype A, circulating in West Africa. A prospective study of a cohort of female sex workers (FSW) in Dakar, Senegal over an 18-year period indicated that an A3-specific sequence in the C2-V3 region of the env gene was found in 46 HIV-1-infected women. HIV-1 sub-subtype A3 appeared in the FSW population as early as 1988 and continued to be transmitted as of 2001. We also found that HIV-1 A3 is not confined to the FSW cohort in Senegal but is also circulating in the general population in Dakar. Furthermore, analyses of viral sequences from a few other West and Central African countries also demonstrated evidence of HIV-1 A3 sequence in isolates from HIV-1-infected people in Ivory Coast, Nigeria, Niger, Guinea Bissau, Benin, and Equatorial Guinea. Overall, because of the evidence of sub-subtype A3 in the general population in Senegal, as well as in a few neighboring West and Central African countries, along with the increasing incidence of infection with A3-containing viruses in the Dakar high-risk FSW population, we feel that HIV-1 sub-subtype A3 viruses are important to distinguish and monitor.     

Parallel evolution of drug resistance in HIV: failure of nonsynonymous/synonymous substitution rate ratio to detect selection.

Parallel or convergent evolution at the molecular level has been difficult to demonstrate especially when rigorous statistical criteria are applied. We present sequence data from the protease gene from eight patients infected with the human immunodeficiency virus (HIV-1). These patients have been on multiple drug therapies for at least 2 years. We present sequence data from two timepoints: time zero--the initiation of drug therapy--and a subsequent timepoint between 59 and 104 weeks after the initiation of drug therapy. In addition to the sequence data, we present viral load data from both initial and final timepoints. Our phylogenetic analyses indicate significant evolution of virus from initial to final time points, even in three of eight patients who show low viral loads. Of the five patients who escaped drug therapy, identical amino acid replacements were seen in all five patients at two different codon positions, an indication of parallel evolution. We also measured genetic diversity for these patients and found no correlation between genetic diversity and viral load. Finally, we calculated the nonsynonymous and synonymous substitution rates and showed that the ratio of nonsynonymous to synonymous substitution compared to the value of one may be a poor indicator of natural selection.     

Distinct human immunodeficiency virus type 1 subtype A virus circulating in West Africa: sub-subtype A3

Phylogenetic analyses demonstrate significant diversity in worldwide circulating strains of human immunodeficiency virus type 1 (HIV-1). Detailed studies have revealed a complex pattern of intersubtype recombinations, as well as evidence of sub-subtypes circulating in various populations. In this study, we characterized an HIV-1 strain that had previously been identified as a distinct subcluster within the subtype A radiation based on partial sequence data. These viruses were of particular interest given that we recently found that their prevalence was significantly higher in dually infected individuals compared to women who were singly infected with HIV-1. Five viruses isolated from commercial sex workers in Dakar, Senegal, were full-length PCR amplified and sequenced. Phylogenetic analyses indicated that, whereas three of these viruses were closely related and clustered overall within the HIV-1 subtype A radiation, they were distinct from previously characterized sub-subtype A1 and A2 viruses. The clustering pattern was maintained in the individual gag, pol, and env regions of the genome. Distance calculations between these viruses, which we termed A3, and other reference sub-subtype A1 and A2 viruses fell in the range of distances between previously characterized sub-subtype groups. In addition, we found evidence of two A3-containing recombinants in our cohort. These recombinants are mosaics composed of sequence from both sub-subtype A3 and CRF02_AG, the major circulating recombinant form in this West African population. Based on phylogenetic analyses, we propose that the group of viruses found in the Dakar sex worker cohort, previously referred to as HIV-1 A subcluster 2, be referred to as HIV-1 sub-subtype A3.     

Baseline CD4+ T Cell Counts Correlates with HIV-1 Synonymous Rate in HLA-B*5701 Subjects with Different Risk of Disease Progression

HLA-B*5701 is the host factor most strongly associated with slow HIV-1 disease progression, although risk of progression may vary among patients carrying this allele. The interplay between HIV-1 evolutionary rate variation and risk of progression to AIDS in HLA-B*5701 subjects was studied using longitudinal viral sequences from high-risk progressors (HRPs) and low-risk progressors (LRPs). Posterior distributions of HIV-1 genealogies assuming a Bayesian relaxed molecular clock were used to estimate the absolute rates of nonsynonymous and synonymous substitutions for different set of branches. Rates of viral evolution, as well as in vitro viral replication capacity assessed using a novel phenotypic assay, were correlated with various clinical parameters. HIV-1 synonymous substitution rates were significantly lower in LRPs than HRPs, especially for sets of internal branches. The viral population infecting LRPs was also characterized by a slower increase in synonymous divergence over time. This pattern did not correlate to differences in viral fitness, as measured by in vitro replication capacity, nor could be explained by differences among subjects in T cell activation or selection pressure. Interestingly, a significant inverse correlation was found between baseline CD4⁺ T cell counts and mean HIV-1 synonymous rate (which is proportional to the viral replication rate) along branches representing viral lineages successfully propagating through time up to the last sampled time point. The observed lower replication rate in HLA-B*5701 subjects with higher baseline CD4⁺ T cell counts provides a potential model to explain differences in risk of disease progression among individuals carrying this allele.     

HIV-1 Genetic Diversity and Its Impact on Baseline CD4+T Cells and Viral Loads among Recently Infected Men Who Have Sex with Men in Shanghai,China

The HIV-1 epidemic among men who have sex with men (MSM) has been spreading throughout China. Shanghai, a central gathering place for MSM, is facing a continuously increasing incidence of HIV-1 infection. In order to better understand the dynamics of HIV-1 diversity and its influence on patient’s immune status at baseline on diagnosis, 1265 newly HIV-1-infected MSM collected from January 2009 to December 2013 in Shanghai were retrospectively analyzed for genetic subtyping, CD4+T cell counts, and viral loads. HIV-1 phylogenetic analysis revealed a broad viral diversity including CRF01_AE (62.13%), CRF07_BC (24.51%), subtype B (8.06%), CRF55_01B (3.24%), CER67_01B (0.95%), CRF68_01B (0.4%), CRF08_BC (0.08%) and CRF59_01B (0.08%). Twenty-four unique recombination forms (URFs) (1.98%) were identified as well. Bayesian inference analysis indicated that the introduction of CRF01_AE strain (1997) was earlier than CRF07_BC strain (2001) into MSM population in Shanghai based on the time of the most recent common ancestor (tMRCA). Three epidemic clusters and five sub-clusters were found in CRF01_AE. Significantly lower CD4+T cell count was found in individuals infected with CRF01_AE than in those infected with CRF07_BC infection (P<0.01), whereas viral load was significantly higher those infected with CRF01_AE than with CRF07_BC (P<0.01). In addition, the patients with >45 years of age were found to have lower CD4+T cell counts and higher viral loads than the patients with <25 years of age (P<0.05). This study reveals the presence of HIV-1 subtype diversity in Shanghai and its remarkable influence on clinical outcome. A real-time surveillance of HIV-1 viral diversity and phylodynamics of epidemic cluster, patient’s baseline CD4+T cell count and viral load would be of great value to monitoring of disease progression, intervention for transmission, improvement of antiretroviral therapy strategy and design of vaccines.     

Quantifying Selection against Synonymous Mutations in HIV-1 env Evolution

Intrapatient evolution of human immunodeficiency virus type 1 (HIV-1) is driven by the adaptive immune system resulting in rapid change of HIV-1 proteins. When cytotoxic CD8⁺ T cells or neutralizing antibodies target a new epitope, the virus often escapes via nonsynonymous mutations that impair recognition. Synonymous mutations do not affect this interplay and are often assumed to be neutral. We test this assumption by tracking synonymous mutations in longitudinal intrapatient data from the C2-V5 part of the env gene. We find that most synonymous variants are lost even though they often reach high frequencies in the viral population, suggesting a cost to the virus. Using published data from SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) assays, we find that synonymous mutations that disrupt base pairs in RNA stems flanking the variable loops of gp120 are more likely to be lost than other synonymous changes: these RNA hairpins might be important for HIV-1. Computational modeling indicates that, to be consistent with the data, a large fraction of synonymous mutations in this genomic region need to be deleterious with a cost on the order of 0.002 per day. This weak selection against synonymous substitutions does not result in a strong pattern of conservation in cross-sectional data but slows down the rate of evolution considerably. Our findings are consistent with the notion that large-scale patterns of RNA structure are functionally relevant, whereas the precise base pairing pattern is not.     

Direct evidence of lower viral replication rates in vivo in human immunodeficiency virus type 2 (HIV-2) infection than in HIV-1 infection

Studies have shown that human immunodeficiency virus type 2 (HIV-2) is less pathogenic than HIV-1, with a lower rate of disease progression. Similarly, plasma viral loads are lower in HIV-2 infection, suggesting that HIV-2 replication is restricted in vivo in comparison to that of HIV-1. However, to date, in vivo studies characterizing replication intermediates in the viral life cycle of HIV-2 have been limited. In order to test the hypothesis that HIV-2 has a lower replication rate in vivo than HIV-1 does, we quantified total viral DNA, integrated proviral DNA, cell-associated viral mRNA, and plasma viral loads in peripheral blood samples from groups of therapy-na?ve HIV-1-infected (n = 21) and HIV-2-infected (n = 18) individuals from Dakar, Senegal, with CD4(+) T-cell counts of >200/microl. Consistent with our previous findings, total viral DNA loads were similar between HIV-1 and HIV-2 and plasma viral loads were higher among HIV-1-infected individuals. Proportions of DNA in the integrated form were also similar between these viruses. In contrast, levels of viral mRNA were lower in HIV-2 infection. Our study indicates that HIV-2 is able to establish a stable, integrated proviral infection in vivo, but that accumulation of viral mRNA is attenuated in HIV-2 infection relative to that in HIV-1 infection. The differences in viral mRNA are consistent with the differences in plasma viral loads between HIV-1 and HIV-2 and suggest that lower plasma viral loads, and possibly the attenuated pathogenesis of HIV-2, can be explained by lower rates of viral replication in vivo.     

Development of an HIV-1 Subtype Panel in China: Isolation and Characterization of 30 HIV-1 Primary Strains Circulating in China

BackgroundThe complex epidemic and significant diversity of HIV-1 strains in China pose serious challenges for surveillance and diagnostic assays, vaccine development and clinical management. There is a lack of HIV-1 isolates in current canonical HIV-1 subtype panels that can represent HIV-1 diversity in China; an HIV-1 subtype panel for China is urgently needed.MethodsBlood samples were collected from HIV-1 infected patients participating in the drug-resistance surveillance program in China. The samples were isolated, cultured and stored as neat culture supernatant. The HIV-1 isolates were fully characterized. The panel was used to compare 2 viral load assays and 2 p24 assays as the examples of how this panel could be used.ResultsAn HIV-1 subtype panel for China composed of 30 HIV-1 primary strains of four subtypes (B [including Thai-B], CRF01_AE, CRF07_BC and G) was established. The samples were isolated and cultured to a high-titer (10⁶-10⁹ copies/ml)/high-volume (40ml). The HIV-1 isolates were fully characterized by the final viral load, p24 concentration, gag-pol and envC2V3 sequencing, co-receptor prediction, determination of the four amino acids at the tip of the env V3-loop, glycosylation sites in the V3 loop and the drug-resistance mutations. The comparison of two p24 assays and two viral load assays on the isolates illustrated how this panel may be used for the evaluation of diagnostic assay performance. The Pearson value between p24 assays were 0.938. The viral load results showed excellent concordance and agreement for samples of Thai-B, but lower correlations for samples of CRF01_AE.ConclusionThe current panel of 30 HIV-1 isolates served as a basis for the development of a comprehensive panel of fully characterized viral isolates, which could reflect the current dynamic and complex HIV-1 epidemic in China. This panel will be available to support HIV-1 research, assay evaluation, vaccine and drug development.     

High Prevalence of Severe Food Insecurity and Malnutrition among HIV-Infected Adults in Senegal,West Africa

Background

Malnutrition and food insecurity are associated with increased mortality and poor clinical outcomes among people living with HIV/AIDS; however, the prevalence of malnutrition and food insecurity among people living with HIV/AIDS in Senegal, West Africa is unknown. The objective of this study was to determine the prevalence and severity of food insecurity and malnutrition among HIV-infected adults in Senegal, and to identify associations between food insecurity, malnutrition, and HIV outcomes.

Methods

We conducted a cross-sectional study at outpatient clinics in Dakar and Ziguinchor, Senegal. Data were collected using participant interviews, anthropometry, the Household Food Insecurity Access Scale, the Individual Dietary Diversity Scale, and chart review.

Results

One hundred and nine HIV-1 and/or HIV-2 participants were enrolled. The prevalence of food insecurity was 84.6% in Dakar and 89.5% in Ziguinchor. The prevalence of severe food insecurity was 59.6% in Dakar and 75.4% in Ziguinchor. The prevalence of malnutrition (BMI <18.5) was 19.2% in Dakar and 26.3% in Ziguinchor. Severe food insecurity was associated with missing clinic appointments (p = 0.01) and not taking antiretroviral therapy due to hunger (p = 0.02). Malnutrition was associated with lower CD4 cell counts (p = 0.01).

Conclusions

Severe food insecurity and malnutrition are highly prevalent among HIV-infected adults in both Dakar and Ziguinchor, and are associated with poor HIV outcomes. Our findings warrant further studies to determine the root causes of malnutrition and food insecurity in Senegal, and the short- and long-term impacts of malnutrition and food insecurity on HIV care. Urgent interventions are needed to address the unacceptably high rates of malnutrition and food insecurity in this population.     

Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection

To understand the high variability of the asymptomatic interval between primary human immunodeficiency virus type 1 (HIV-1) infection and the development of AIDS, we studied the evolution of the C2-V5 region of the HIV-1 env gene and of T-cell subsets in nine men with a moderate or slow rate of disease progression. They were monitored from the time of seroconversion for a period of 6 to 12 years until the development of advanced disease in seven men. Based on the analysis of viral divergence from the founder strain, viral population diversity within sequential time points, and the outgrowth of viruses capable of utilizing the CXCR4 receptor (X4 viruses), the existence of three distinct phases within the asymptomatic interval is suggested: an early phase of variable duration during which linear increases ( approximately 1% per year) in both divergence and diversity were observed; an intermediate phase lasting an average of 1.8 years, characterized by a continued increase in divergence but with stabilization or decline in diversity; and a late phase characterized by a slowdown or stabilization of divergence and continued stability or decline in diversity. X4 variants emerged around the time of the early- to intermediate-phase transition and then achieved peak representation and began a decline around the transition between the intermediate and late phases. The late-phase transition was also associated with failure of T-cell homeostasis (defined by a downward inflection in CD3(+) T cells) and decline of CD4(+) T cells to 相似文献   

Genetic heterogeneity of hypervariable region 1 of the hepatitis C virus (HCV) genome and sensitivity of HCV to alpha interferon therapy

Hepatitis C virus (HCV) populations persist in vivo as a mixture of heterogeneous viruses called quasispecies. The relationship between the genetic heterogeneity of these variants and their responses to antiviral treatment remains to be elucidated. We have studied 26 virus strains to determine the influence of hypervariable region 1 (HVR-1) of the HCV genome on the effectiveness of alpha interferon (IFN-alpha) therapy. Following PCR amplification, we cloned and sequenced HVR-1. Pretreatment serum samples from 13 individuals with chronic hepatitis C whose virus was subsequently eradicated by treatment were compared with samples from 13 nonresponders matched according to the major factors known to influence the response, i.e., sex, genotype, and pretreatment serum HCV RNA concentration. The degree of virus variation was assessed by analyzing 20 clones per sample and by calculating nucleotide sequence entropy (complexity) and genetic distances (diversity). Types of mutational changes were also determined by calculating nonsynonymous substitutions per nonsynonymous site (K(a)) and synonymous substitutions per synonymous site (K(s)). The paired-comparison analysis of the nucleotide sequence entropy and genetic distance showed no statistical differences between responders and nonresponders. By contrast, nonsynonymous substitutions were more frequent than synonymous substitutions (P 相似文献   

Independent variation and positive selection in env V1 and V2 domains within maternal-infant strains of human immunodeficiency virus type 1 in vivo.

Multiple targets for immune recognition and cellular tropism are localized to the V1 and V2 hypervariable regions in the amino portion of human immunodeficiency virus type 1 (HIV-1) gp120env. We have assessed genetic diversity in env V1 and V2 hypervariable domains in vivo within epidemiologically related strains of HIV-1. Our strategy was to analyze longitudinal samples from two seropositive mothers and multiple children infected by perinatal transmission. Although the V1 and V2 domains are closely linked in the HIV-1 genome, nucleotide sequences in V1 and in V2 evolved independently in maternal-infant viruses in vivo. A high proportion of the nucleotide substitutions would introduce amino acid diversity in V1 and in V2. A significant excess of nonsynonymous over synonymous substitutions was identified in HIV-1 env V1 and V2 peptides in the mothers and in two older children but was not generally apparent in HIV-1 sequences in infants. An excess of nonsynonymous over synonymous substitutions indicated that there is positive selection for independent genetic variation in the V1 and V2 domains in vivo. It is likely that there are host responses to complex determinants in the V1 or V2 hypervariable domain of HIV-1 gp120.     

Founder virus population related to route of virus transmission: a determinant of intrahost human immunodeficiency virus type 1 evolution?

We and others have shown that in individual human immunodeficiency virus type 1 (HIV-1) infection, the adaptive evolution of HIV-1 is influenced by host immune competence. In this study, we tested the hypothesis that in addition to selective forces operating within the host, transmission bottlenecks have an impact on HIV-1 intrahost evolution. Therefore, we studied the intrahost evolution of the V3 region of the external glycoprotein gp120 of HIV-1 during the 3- and 5-year periods following seroconversion after parenteral versus sexual (male-to-male) transmission in 41 participants of the Amsterdam prospective cohorts of homosexual men (n = 31) and intravenous drug users (IVDUs; n = 10) who were AIDS free and had comparable numbers of CD4+ cells. We observed that HIV-1 strains in homosexual men accumulated over 5 years more nonsynonymous substitutions within the V3 loop than HIV-1 strains in IVDUs as a result of lower rates of nonsynonymous evolution in both the initial 3-year period from seroconversion and the following 2-year period as well as a larger proportion of nonsynonymous back substitutions in IVDUs. The mean numbers of synonymous substitutions did not differ between the two risk groups. Since HIV-1 strains in IVDUs could be distinguished from the viruses of homosexual men based on several nucleotide substitutions of which the most conserved is a synonymous substitution at the tip of the V3 loop (GGC pattern), we studied whether the founder virus population itself has an impact on the intrahost evolution of HIV-1. The mean number of nonsynonymous substitutions accumulated over 5 years within the V3 loop was lower in 10 IVDUs infected by the HIV-1 strains with the GGC signature than in 4 IVDUs infected by HIV-1 strains lacking this pattern, while the mean numbers of synonymous substitutions were similar in the two groups.     

Infection with multiple human immunodeficiency virus type 1 variants is associated with faster disease progression

Human immunodeficiency virus type 1 (HIV-1)-infected individuals develop a genetically diverse virus population over time, but often only a limited number of viral variants are transmitted from a chronic carrier to a newly infected person. Interestingly, many women but few men are infected by multiple HIV-1 variants from a single partner. To determine whether the complexity of the infecting virus population influences clinical outcome, we examined viral diversity in the HIV-1 envelope sequences present at primary infection in 156 women from Kenya for whom we had follow-up data on viral RNA levels and CD4 T-cell counts. Eighty-nine women had multiple viral genotypes, while 67 women had a single genotype at primary infection. Women who acquired multiple viral genotypes had a significantly higher viral load (median, 4.84 versus 4.64 log(10) copies/ml, P = 0.04) and a significantly lower CD4(+)-T-cell count (median, 416 versus 617 cells/mm(3), P = 0.01) 4 to 24 months after infection compared to women who were infected with a single viral genotype. These studies suggest that early HIV-1 genetic diversity is linked to faster disease progression.     

Host-specific modulation of the selective constraints driving human immunodeficiency virus type 1 env gene evolution

To address the evolution of human immunodeficiency virus type 1 (HIV-1) within a single host, we analyzed the HIV-1 C2-V5 env regions of both cell-free genomic-RNA- and proviral-DNA-derived clones. Sequential samples were collected over a period of 3 years from six untreated subjects (three typical progressors [TPs] and three slow progressors [SPs], all with a comparable length of infection except one. The evolutionary analysis of the C2-V5 env sequences performed on 506 molecular clones (253 RNA- and 253 DNA-derived sequences) highlighted a series of differences between TPs and SPs. In particular, (i) clonal sequences from SPs (DNA and RNA) showed lower nucleotide similarity than those from TPs (P = 0. 0001), (ii) DNA clones from SPs showed higher intra- and intersample nucleotide divergence than those from TPs (P < 0.05), (iii) higher host-selective pressure was generally detectable in SPs (DNA and RNA sequences), and (iv) the increase in the genetic distance of DNA and RNA sequences over time was paralleled by an increase in both synonymous (Ks) and nonsynonymous (Ka) substitutions in TPs but only in nonsynonymous substitutions in SPs. Several individual peculiarities of the HIV-1 evolutionary dynamics emerged when the V3, V4, and V5 env regions of both TPs and SPs were evaluated separately. These peculiarities, probably reflecting host-specific features of selective constraints and their continuous modulation, are documented by the dynamics of Ka/Ks ratios of hypervariable env domains.     

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.