JOURNAL OF VIROLOGY, Nov. 2004, p. 12438–12445 0022-538X/04/$08.00⫹0 DOI: 10.1128/JVI.78.22.12438–12445.2004 Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Vol. 78, No. 22
Distinct Human Immunodeficiency Virus Type 1 Subtype A Virus Circulating in West Africa: Sub-Subtype A3 Seema Thakore Meloni,1 Bohye Kim,2 Jean-Louis Sankale´,1 Donald J. Hamel,1 Sodsai Tovanabutra,2 Souleymane Mboup,3 Francine E. McCutchan,2 and Phyllis J. Kanki1* Department of Immunology and Infectious Diseases, Harvard School of Public Health, Boston, Massachusetts1; and Henry M. Jackson Foundation, Rockville, Maryland2; and Laboratoire de Bacte´riologie-Virologie, Universite´ Cheikh Anta Diop, Dakar, Senegal3 Received 16 February 2004/Accepted 22 June 2004
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. Subtype nomenclature provides us with an important tool that enables us to monitor the spread of human immunodeficiency virus type 1 (HIV-1), and the creation of an international HIV genetic sequence database has facilitated better classification of globally circulating viral strains. Phylogenetic analyses of the env and gag genes of HIV-1 indicate that there are three distinct groups (M, N, and O) of the virus circulating worldwide. Within the major group M, there are currently nine defined subtypes, or clades (A to D, F to H, J, and K) (3, 6, 18, 19, 23, 24, 34). In addition, there are a number of circulating recombinant forms (CRFs) in the world, which are associated with significant transmission and disease in various populations worldwide (23). Detailed phylogenetic studies have revealed evidence of sub-subtypes, which are characterized as regional subclusters composed of distinct lineages that are closely related to a given subtype (11, 21, 36–38). Sub-subtypes F2 (36) and A2 (11) have been identified and characterized. There have also been suggestions that certain previously identified subtypes could have been categorized as sub-subtypes (19, 26, 36–38). For example, there is evidence suggesting that the previously identified subtype D viruses might have been more appropriately named sub-subtype B2 (26, 38). Due to the complexity of the nomenclature issues, an international group of experts convened to define methods of clas* Corresponding author. Mailing address: Department of Immunology and Infectious Diseases, Harvard School of Public Health, 651 Huntington Ave., Boston, MA 02115. Phone: (617) 432-1267. Fax: (617) 432-3575. E-mail:
[email protected].
sifying HIV viral sequences (23, 26, 27). Emphasis was placed on improving precision with regard to subtype classifications. The group pointed out the importance of recognizing the presence of sub-subtypes. In their recommendations, they mapped out criteria for defining sub-subtypes by genetic distance estimates. Within each of the subtypes and sub-subtypes identified, a range of genetic diversity has been noted (1, 11, 13, 18, 30). Subtypes are genetically equidistant from one another, with the intersubtype nucleotide distances ranging from 15 to 22% in the gag gene and from 20 to 30% in the env gene (7, 11, 12, 26). Intrasubtype distances range from 3 to 10% in the gag gene and 5 to 12% in the env gene (7, 11, 26). Various reports indicate that intersub-subtype diversity ranges from 7 to 12% in the gag gene from 11 to 16% in the env gene (11, 26, 27). Subtype A viruses are interesting to study not only because they are the second most prevalent strain of HIV-1 circulating globally but also because portions of the subtype are found in a majority of the CRFs, particularly in those of epidemiologic consequence (11, 23, 25, 26). This group of viruses is estimated to account for 27% of global HIV-1 infections, with a majority of the cases represented by the CRF02㛭AG strains (23). HIV-1 infection in West Africa is principally due to CRF02㛭AG (2, 4, 6, 19, 22, 23, 25). Recent work conducted in our laboratory revealed that a distinct group of viruses within the HIV-1A subtype radiation is also present in Senegal (30, 32, 35). Phylogenetic analyses indicated that this group of viruses, which we previously referred to as HIV-1A subcluster 2, differed from sub-subtypes A1 and A2 (30, 32). Since the categorization of these viruses was primarily based on sequencing of the
12438
VOL. 78, 2004
HIV-1 SUB-SUBTYPE A3
12439
TABLE 1. Characteristics of samples from Dakar female sex worker cohort selected for full-length sequencinga Sample identification
HIV-1 case type
DDJ360 DDI579 DDJ369 DDJ362 DDJ364
Incident Incident Incident Prevalent Prevalent
Serohistory
HIV-1 HIV-1 HIV-1 HIV-1 HIV-1
Yr of HIV-1 seroconversion
Sample yr
env sub-subtype
gag sub-subtype
1991 1994 2001 1989ⴱ 1989ⴱ
1996 2001 2001 1996 1996
A3ⴱⴱ A3ⴱⴱ A3ⴱⴱ A3ⴱⴱ A3ⴱⴱ
A3ⴱⴱ A3ⴱⴱ A3ⴱⴱ A3ⴱⴱ A3ⴱⴱ
only only only and HIV-2 and HIV-2
a ⴱ, Since these women were already seroprevalent at the time they entered the study, the dates correspond to the year of first bleed; ⴱⴱ, at the time these samples were first identified for full-length sequencing, we had not yet verified that they could be referred to as sub-subtype A3 (and were referring to them as HIV-1 A subcluster 2 viruses), but for the purposes of this table we are designating the samples as A3.
C2-V3 region of the env gene and a small region of the gag gene, further analysis was warranted. In addition, we were particularly interested in further characterization of this virus after we found that it was actually overrepresented in women dually infected with HIV-1 and HIV-2 compared to women infected with HIV-1 alone (32). The goal of this investigation was to obtain additional evidence to determine the genetic relatedness and relationship of this subcluster of subtype A viruses with respect to reference HIV-1 A1 and HIV-1 A2 viruses; we therefore analyzed full-length HIV-1 sequences. Sequences from five different subcluster two viruses were obtained and analyzed. Pairwise genetic distances were evaluated and revealed that three of the five viruses could be designated a distinct sub-subtype within the HIV-1 A radiation, which we have named A3.
Sequence analysis and statistical methods. A multiple alignment of the virtually full-length genome with reference sequences of all major subtypes and sub-subtypes (A1, A2, B-D, F1, F2, G, H, J, and K), as well as the full-length sequences we had obtained, was manually generated in a Genetic Data Environment window (33). A maximum-likelihood tree of nearly full-length sequences was produced by using the DNAML program in the PHYLIP package version 3.6 (9) using a transition/transversion ratio of 1.42 (16) and empirical base frequencies. Sequences were gap stripped. Maximum-likelihood trees for gag, pol, and env subregions were also generated. The sequences were also analyzed for hypermutation by using the HYPERMUT program (http://www.hiv .lanl.gov). Pairwise genetic distances were calculated by using the DNADIST program in the PHYLIP version 3.6 package (9), with an F84 model of evolution and a transition/transversion ratio of 1.42. Calculations were done with sequences from the following subtype and sub-subtype categories: A1 (KEQ23, UG037), A2 (94cy017, 97cdktb48), B (HXB2, SF2), C (ETH2220, 93IN905), D (NDK, KE2059), F1 (VI850, F9363), F2 (CAM53657, MP255), G (NG093, DRC.BL), H (VI991, CF056), J (SE9173㛭3, SE9280), and K (EGTB117, MP535c). The alignments were gap stripped.
MATERIALS AND METHODS Study population and sample collection. Since 1985, we have conducted a prospective study of registered female sex workers (FSWs) in Dakar, Senegal; blood samples and questionnaire data were collected after obtaining informed consent. The specific details of the study recruitment procedures and methods have been previously described elsewhere (14, 15). Serostatus was determined by immunoblot on whole virus lysates and recombinant envelope peptides and by diagnostic HIV-1 and/or HIV-2 PCR. Time of infection for women who converted to HIV-positive serostatus while in the study was estimated as the midpoint between their last seronegative and first seropositive samples. DNA amplification and sequencing. Proviral DNA was extracted from cocultured peripheral blood mononuclear cells by using a kit from Qiagen, Inc. (Chatsworth, Calif.). Preliminary subtype and sub-subtype designation was done based on the C2-V3 region of the env gene (350 bp). For C2-V3 sequence, DNA samples were amplified by nested PCR with primers and reaction conditions that have been described previously (31, 39). The PCR products were purified and directly sequenced by using the second-round primers. When sequence data could not be obtained through direct sequencing, purified products were cloned into the pCR2.1 vector (T/A Cloning; Invitrogen, San Diego, Calif.), and sequence data were generated from clones. For women who appeared to have evidence of the HIV-1 subtype A subcluster 2 virus in the C2-V3 region of the env (including AF085293, AF085294, AF085299, AF085309, AF085310, AF085319, and AF085326), sequencing of the 3⬘-p24/5⬘-p7 region of the gag (495 bp) was also performed. The gag sequences were amplified by nested PCR with two sets of primers: p108-p109 for the first round and p91-p92 for the second round, as previously described (30). These products were directly sequenced by using second-round primers. All sequences for both strands of DNA were determined by dye terminator cycle sequencing, using Taq polymerase (Perkin-Elmer Applied Biosystems Division, Foster City, Calif.) and an automatic sequencer (ABI 377; Perkin-Elmer Applied Biosystems Division). Based on preliminary env and gag data, samples from five HIV-1-infected subjects were selected for full-length sequencing. Nearly the entire genome was obtained in one continuous segment by using the nested primers MSF12bOFMR1, for round 1, and F2nst-ofm19, for round 2, utilizing methods previously described (4, 29). The amplicon was directly sequenced by using Big Dye terminators on an ABI 3100 capillary sequencer (Applied Biosystems), according to manufacturer’s instructions. Sequences were assembled by using Sequencher software version 3.1 (Genecodes, Inc., Ann Arbor, Mich.).
FIG. 1. Maximum-likelihood tree of full-length sequences including three A3 sequences (DDI579, DDJ360, and DDJ369) and reference sequences from all known subtypes and sub-subtypes (A1, A2, B, C, D, F1, F2, G, H, J, and K; Los Alamos HIV Sequence database [http://hiv-web.lanl.gov]). The alignments were gap stripped, and a transition/transversion ratio of 1.42 was used. The scale bar indicates 1% nucleotide sequence divergence. The A3 sequences are indicated in boldface.
12440
MELONI ET AL.
J. VIROL.
FIG. 2. Maximum-likelihood subregion trees. (a) gag; (b) pol; (c) env. Sequences maintain unique clustering patterns in individual gene regions. The alignments were gap stripped, and these trees were generated by using the DNAML program in the PHYLIP version 3.6 package (9). A transition/transversion ratio of 1.42 was used for all trees. The scale bar indicates 1% nucleotide sequence divergence. For these trees, reference sequences from all major subtypes and sub-subtypes (A1, A2, B, C, D, F1, F2, G, H, J, and K) were used. The A3 sequences are indicated in boldface.
Sequence alignments were also used to generate bootscan graphs for recombination analyses by using the SIMPLOT version 3.4 beta program (17). Bootscan analyses were performed on neighbor-joining (NJ) trees, using a 400-bp moving window along the aligned sequences in 50-bp increments. A transition/transversion ratio of 1.42 was used with an F84 model of evolution. Gaps were stripped from the alignments. In preliminary analyses, reference sequences from all major subtypes, sub-subtypes, and CRF02㛭AG were used. Based on recombination patterns from the initial analyses, a subset of sequences was selected for the final bootscan graphs; the subset of reference sequences used included sub-subtype A3 (DDI578, DDJ360, and DDJ369), CRF02㛭AG (DJ264, IbNG, SE7812, MP1211, and MP807), subtype G (G6165, HH8793, NG083, and DRC.BL), and subtype B (HXB2, RF, SF2, and RL42) as the outgroup. The reference sequences were grouped and the final bootscan plot was generated by using the consensus for the reference sequence groups. One hundred bootstrap replicates were analyzed for the plots. Bootstrap values were plotted at the midpoint of each window along the sequence. Breakpoints were determined by using informative site analyses and by maximizing 2 values (28). Using the breakpoints, sequence alignments were clipped into appropriate subregions for further analysis. NJ trees for each of the recombined subregions were generated to confirm the results of the bootscan. All NJ trees were generated by the SEQBOOT, DNADIST, NEIGHBOR, and CONSENSE programs in the PHYLIP version 3.6 package (9). For all NJ trees, the F84 model of evolution and a transition/transversion ratio of 1.42 was used.
Sequence information. The full-length sequences have been deposited into the GenBank database under accession numbers AY521629 to AY521633.
RESULTS We amplified nearly full-length genome of viruses from five women in our cohort that we had previously identified as having evidence of the HIV-1 subtype A subcluster 2 virus in the C2-V3 region of the env gene and the 3⬘-p24/5⬘-p7 region of the gag gene (Table 1); in the table, we referred to these sequences as A3. Three of the women (DDI579, DDJ360, and DDJ369) had seroincident infections and were infected with HIV-1 only. The remaining two women (DDJ362 and DDJ364) we selected for sequencing were seroprevalent upon entry into the cohort and were dually infected with HIV-1 and HIV-2. Characterization of HIV-1 A subcluster 2 as sub-subtype A3. To understand the phylogenetic relationships of the newly characterized viruses, we constructed trees from the near fulllength genome. Initial phylogenetic tree analyses with se-
VOL. 78, 2004
quences from all five women revealed that viruses from two of the women (DDJ362 and DDJ364) had branching patterns that suggested they might be recombinants; the details of these analyses are discussed in greater detail below. Figure 1 presents a maximum-likelihood tree with the three remaining sequences (DDI579, DDJ360, and DDJ369) that contained unique sequence across the length of the entire genome. These sequences clustered overall with HIV-1 subtype A viruses but formed a subcluster that was distinct from A1 sequences as well as from the recently described A2 viruses. It should be noted that because of its unique branching pattern within the subcluster, DDJ369 was also analyzed as a possible recombinant. Bootscan analyses revealed that DDJ369 did not appear to be a recombinant and clustered with the other A3 sequences across the entire length of the genome (data not shown). Finally, for the three sequences with distinct sequence across the length of the genome, all reading frames were open and of complete length. None of the genomes had major insertions or deletions. For sequence DDJ369, the predicted start codon of the env gene had a transition from ATG to ACG, thus changing the amino acid from methionine (M) to a threonine (T); DDJ369 has a methionine three codons downstream, which might serve as the actual start codon for the env. Overall, when the sequences were trimmed into subregions based on the open reading frames for major genes, we found that the sequences maintained their clustering pattern (Fig. 2). Finally, the sequences were analyzed by using the HYPERMUT program (http://www.hiv.lanl.gov), and no evidence of hypermutation was found. Pairwise distance analysis has previously been used to distinguish between subtypes and sub-subtypes (1, 11, 26, 37). Intersubtype and intrasubtype distances were calculated for the full-length sequences based on reference strains (Fig. 3). These values are useful in determining whether a lineage is more likely a subtype or a sub-subtype (27). By convention, distance values that were less than the distance between different subtypes but greater than the distance within subtype categories were considered to be suggestive of a sub-subtype (4). The genetic distances between sub-subtypes A3 and A1, as well as A3 and A2, fall in the range of 9 to 13%, which is within the range of distances seen between other sub-subtypes (A1:A2 and F1:F2) and subtypes that should have been designated sub-subtypes (B:D; Fig. 3). Overall, all of the sub-subtype genetic distance values fell between the values for the intersubtype and intrasubtype comparisons. Identification of recombinant genomes. As mentioned above, two of the five samples that we had full-length sequenced revealed distinct clustering patterns when originally put into a phylogenetic tree. We generated two separate maximum-likelihood trees to highlight the clustering pattern of sequences DDJ362 and DDJ364 (Fig. 4a and b). For the first (Fig. 4a), the maximum-likelihood tree was generated by using major subtype and sub-subtype reference sequences, but no CRF02㛭AG reference sequences. In this tree, DDJ362 and DDJ364 clustered together overall with A3 sequences but formed a small subgroup. However, when CRF02㛭AG was added to the tree (Fig. 4b), the DDJ364 sequence clustered with CRF02㛭AG sequences. Although the DDJ362 sequence continued to cluster with the A3 sequences, the branching pattern suggested that it might be a recombinant sequence.
HIV-1 SUB-SUBTYPE A3
12441
FIG. 3. Plots of intrasubtype, intersubtype, and inter-sub-subtype genetic distances for full-length sequences. To calculate these distances, a pairwise distance matrix was constructed by using DNADIST from the PHYLIP version 3.6 package (9), with the F84 model of evolution, a transition/transversion ratio of 1.42, and empirical base frequencies. All alignments were gap stripped. CRFs were excluded from these calculations. Subtypes B and D were evaluated separately from all of the other intersubtype comparisons. Inter-sub-subtype distances were computed between sub-subtypes A1 and A2, sub-subtypes F1 and F2, and subtypes B and D. The value 0.01 on the x axis indicates a 1% nucleotide sequence divergence. The y axis shows the percentage of pairwise comparisons that have the same nucleotide sequence divergence.
Accordingly, we investigated these sequences for possible recombination events. Bootscan plots were generated, and sequence breakpoints were determined by using the bootscan results (17, 28). Initial similarity plot and bootscan analyses included representatives from all of the major subtypes and sub-subtypes (A1, A2, A3, B, C, D, F1, F2, G, J, and K), as well as CRF02. The final bootscan plots were generated by using representative CRF02㛭AG (DJ264, IbNG, SE7812, MP1211, and MP807), sub-subtype A3 (DDI579, DDJ360, and DDJ369), and subtype G (G6165, HH8793, NG083, and DRC.BL) sequences. Sequences from HIV-1 subtype B (HXB2, RF, SF2, and RL42) were also incorporated in the analyses as the outgroup. This assessment indicated that the sequences of DDJ362 and DDJ364 were a mosaic of CRF02㛭AG and HIV-1 A3. The analyses indicated that the sequence from subject DDJ362 grouped with high bootstrap values to A3 sequence in the gag gene, the central region of the pol gene and the gp120 portion of the env gene. The remaining portions of the genome clustered with high bootstrap values with CRF02㛭AG sequence
12442
MELONI ET AL.
J. VIROL.
FIG. 4. Recombination analyses for two sequences previously designated as containing A3 in the C2-V3 region of the env. An alignment of nearly full-length strains was used; long terminal repeats were not analyzed. Alignments were gap stripped prior to analysis. (a) Maximumlikelihood tree containing references from all subtypes and sub-subtypes (A1, A2, B, C, D, F1, F2, G, H, J, and K), A3 sequences (DDI579, DDJ360, and DDJ369) and the two sequences assumed to be recombinants (DDJ362 and DDJ364). (b) Maximum-likelihood tree containing references from all subtypes and sub-subtypes (A1, A2, B, C, D, F1, F2, G, H, J, and K), A3 sequences (DDI579, DDJ360, and DDJ369), CRF02 sequences and the two sequences assumed to be recombinants (DDJ362 and DDJ364); both maximum-likelihood trees were generated by using DNAML from the PHYLIP version 3.6 package (9), using a transition/transversion ratio of 1.42. (c and d) Recombination analyses for sequences from subjects DDJ362 and DDJ364, respectively. The recombinant regions and corresponding breakpoints were determined by using bootscan plots (28). Bootscans were generated on NJ trees of gap-stripped sequences, using a transition/transversion ratio of 1.42. Final plots for both recombinants contain reference sequences from sub-subtype A3 (DDI578, DDJ360, and DDJ369), CRF02㛭AG (DJ264, IbNG, SE7812, MP1211, and MP807), subtype G (G6165, HH8793, NG083, and DRC.BL), and subtype B (HXB2, RF, SF2, and RL42) as the outgroup. The reference sequences were grouped, and the final bootscan plot was generated by using the consensus for the reference sequence groups. Bootscans were performed with a moving window of 400 bp and a step size of 50 bp. Breakpoints were determined by using informative site analyses and by maximizing 2 values. For both samples, bootscan plots are shown. Below the plots, are schematics of the predicted recombinant sequence patterns, as well as the NJ phylogenetic trees for each of the recombined subregions within the plot. For each of the sequences, the roman numerals on the bootscan correspond to the roman numerals above the phylogenetic tree. All NJ trees were generated by using the SEQBOOT, DNADIST, NEIGHBOR, and CONSENSE programs in the PHYLIP version 3.6 package (9). For all NJ trees, the F84 model of evolution and a transition/transversion ratio of 1.42 was used.
(Fig. 4c). Similarly, the viral sequence from DDJ364 was composed of a mosaic of A3 and CRF02㛭AG. As shown in Fig. 4d, sequence from this subject clustered with high bootstrap values with reference A3 in the V1 to V5 regions of env. In addition, the sequence clustered with A3 in portions of the gag, pol, and nef genes. DISCUSSION Through phylogenetic analyses, we confirmed that the subcluster of viruses circulating in the Dakar sex worker cohort differs from the previously characterized HIV-1 sub-subtypes A1 and A2. The subjects chosen for more detailed analyses were previously characterized as being infected with HIV-1 subtype A subcluster 2 viruses (30, 32). Examination of fulllength sequences revealed that three subjects had sequences that form a distinct subcluster. Furthermore, phylogenetic analyses of subregions verified that these viruses maintained their distinct clustering patterns in the major genes. Genetic
distance analysis has previously been identified as a valid method of distinguishing between the subtypes and sub-subtypes of HIV-1. An analysis of nucleotide distances revealed that our A3 sequences differed from previously identified and characterized A1 and A2 sequences by approximately 9 to 13% across the length of the entire genome. This difference was comparable to those found between other sub-subtypes, such as F1:F2, and subtypes B:D. Based on phylogenetic information, we believe that the group of viruses in our cohort previously designated a subcluster within HIV-1 A can be considered a sub-subtype of HIV-1 A. We propose that the subcluster be referred to as HIV-1 sub-subtype A3. When looking at our three sequences that clustered together as sub-subtype A3, we found that the virus from DDJ369 was more distant from those from DDI579 and DDJ360, respectively, than they were from each other. We conducted bootscan analyses and ruled out that the sequence from subject DDJ369 was a recombinant. This finding suggests that either diversity
VOL. 78, 2004
HIV-1 SUB-SUBTYPE A3
FIG. 4—Continued.
12443
12444
MELONI ET AL.
has accrued in the A3 sequences in Senegal over the years or that there might have been more than one introduction of A3 into Senegal, such that the virus diversified prior to its importation into Senegal. Viral recombination is generally seen in areas that have several subtypes or recombinant forms of the virus circulating in the population (5). In areas where multiple subtypes cocirculate, it is estimated that intersubtype recombinant sequences (ISRs) account for between 8 and 24% of infections (23). Senegal has a high prevalence of CRF02㛭AG infections (4, 24, 25, 30, 32) and recombination with this form of the virus could have been predicted. Given that subtype A viruses have been identified in a majority of the circulating intersubtype recombinants, as well as both intergroup (M/O) recombinants (10, 20, 24), it is also not surprising to find that the HIV-1 A3 viruses are involved in recombination events. In the present study, we found two viruses with recombinant genomes, where the HIV-1 A3 sequence was maintained in the env gene. Interestingly, the two subjects that had recombinant viruses also happened to be dually infected, with HIV-1 and HIV-2. However, given that both subjects DDJ362 and DDJ364 were seroprevalent with both HIV-1 and HIV-2 upon entry into the cohort, it is difficult to speculate on the implications of infection with these HIV-1 recombinants. Another point to note is that these women were infected with A3-containing recombinant viruses as early as 1989. That suggests that A3 might have been in the cohort prior to 1989 for recombination events to take place. Alternatively, the virus might have recombined in another country prior to 1989 and then been imported into the Dakar FSW cohort already a recombinant. However, given we have no information regarding the actual time of infection for these two women, we cannot make assumptions as to when A3 was first introduced into the population. Ultimately, the discovery that the viruses from these two women were recombinants reemphasizes the importance of full-length sequencing in subtype and sub-subtype characterization; had we based our sub-subtype designations on env and gag results only, these two women would have originally been characterized as infected with A3. A clearer picture of sequence diversity in a given population will lead to better understanding and insights into the biologic consequences of viral diversity. Identification of the different subtypes and sub-subtypes of HIV-1 circulating in a given region will aid in molecular surveillance and development of relevant vaccine candidates (4). Since sequence changes are thought to represent, in part, escape from selective pressure, genetic analyses and molecular surveillance are critical to understanding host immune response (21). In addition, as subtype-specific vaccines are being produced and tested, it is imperative that we understand the scope of diversity. Phylogenetic analyses have revealed a significant amount of sequence divergence within each of the subtypes, and it may well be that understanding regional differences is requisite for subtype-specific vaccines to be effective (8). In summary, our data indicate that the women previously identified as being infected by HIV-1 subtype A subcluster 2 viruses can be classified as infected with HIV-1 sub-subtype A3. The phylogenetic analyses of the full-length genomes corroborate the findings from the original studies (30, 32). Continued monitoring will indicate the future role of this sub-
J. VIROL.
subtype in the epidemic and will demonstrate the contribution of HIV-1 sub-subtype A3 to recombinant viruses. Further molecular epidemiologic studies will require additional discrimination of this particular sub-subtype from other circulating subtype A viruses to identify its role in the epidemic. ACKNOWLEDGMENTS We thank Abdoulaye Dieng Sarr, Christopher Mullins, Aissatou Gue`ye-Ndiaye, Ibrahima Traore´, Mamadou Cire´ Dia, and Ousmane Diouf for their continued assistance and support on this project. We are extremely grateful for all of the assistance provided by the McCutchan lab in generating the full-length sequences. This study was supported in part by U.S. Army grant DAMD 1795-C-5005 and by NIH grants AI43879 and AI467274. REFERENCES 1. Anderson, J., A. Rodrigo, G. Learn, A. Madan, C. Delahunty, M. Coon, M. Girard, S. Osmanov, L. Hood, and J. Mullins. 2000. Testing the hypothesis of a recombinant origin of human immunodeficiency virus type 1 subtype E. J. Virol. 74:10752–10765. 2. Andersson, S., H. Norrgren, F. Dias, G. Biberfeld, and J. Albert. 1999. Molecular characterization of human immunodeficiency virus (HIV)-1 and -2 individuals from Guinea-Bissau with single or dual infections: predominance of a distinct HIV-1 subtype A/G recombinant in West Africa. Virology 267:312–330. 3. Burke, D., and F. McCutchan. 1997. Global distribution of human immunodeficiency virus 1 clades, p. 119–126. In T. Vincent, J. DeVita, S. Hellman, and S. Rosenberg (ed.), AIDS: biology, diagnosis, treatment, and preventions, 4th ed. Lippincott-Raven Publishers, Philadelphia, Pa. 4. Carr, J., T. Laukkanen, M. Salminen, J. Albert, A. Alaeus, B. Kim, E. Sanders-Buell, D. Birx, and F. McCutchan. 1999. Characterization of subtype A HIV-1 from Africa by full genome sequencing. AIDS 13:18189– 18226. 5. Cornelissen, M., G. Kampiga, F. Zodrager, and J. Goudsmit. 1996. UNAIDS network for HIV isolation and characterization: human immunodeficiency virus type 1 subtypes defined by env show high frequency of recombinant gag genes. J. Virol. 70:8209–8212. 6. Cornelissen, M., R. vandenBurg, F. Zorgdrager, and J. Goudsmit. 2000. Spread of distinct human immunodeficiency virus type 1 AG recombinant lineages in Africa. J. Gen. Virol. 81:515–523. 7. Dowling, W. E., B. Kim, C. J. Mason, K. M. Wasunna, U. Alam, L. Elson, D. L. Birx, M. L. Robb, F. E. McCutchan, and J. K. Carr. 2002. Forty-one near full-length HIV-1 sequences from Kenya reveal an epidemic of subtype A and A-containing recombinants. AIDS 16:1809–1820. 8. Ellenberger, D. L., B. Li, L. D. Lupo, S. M. Owen, J. Nkengasong, M. S. Kadio-Morokro, J. Smith, H. Robinson, M. Ackers, A. Greenberg, T. Folks, and S. Butera. 2002. Generation of a consensus sequence from prevalent and incident HIV-1 infections in West Africa to guide AIDS vaccine development. Virology 302:156–163. 9. Felsenstein, J. 1992. PHYLIP (phylogenetic inference package), version 3.5c ed. Department of Genetics, University of Washington, Seattle. 10. Gao, F., D. Robertson, S. Morrison, H. Hi, S. Graig, J. Decker, P. Fultz, M. Girard, G. Shaw, B. Hahn, and P. Sharp. 1996. The heterosexual human immunodeficiency virus type 1 epidemic in Thailand is caused by an intersubtype (A/E) recombinant of African origin. J. Virol. 70:7013–7029. 11. Gao, F., N. Vidal, Y. Li, S. A. Trask, Y. Chen, L. G. Kostrikis, D. D. Ho, J. Kim, M. D. Oh, K. Choe, M. Salminen, D. L. Robertson, G. M. Shaw, B. H. Hahn, and M. Peeters. 2001. Evidence of two distinct subsubtypes within the HIV-1 subtype A radiation. AIDS Res. Hum. Retrovir. 17:675–688. 12. Janssens, W., A. Buve, and J. N. Nkengasong. 1997. The puzzle of HIV-1 subtypes in Africa. AIDS 11:705–712. 13. Kaleebu, P., A. Ross, D. Morgan, D. Yirrell, J. Oram, A. Rutebemberwa, F. Lyagoba, L. Hamilton, B. Biryahwaho, and J. Whitworth. 2001. Relationship between HIV-1 env subtypes A and D and disease progression in a rural Ugandan cohort. AIDS 15:293–299. 14. Kanki, P., S. Mboup, R. Marlink, K. Travers, C. C. Hsieh, A. Gueye, C. Boye, J. L. Sankale, C. Donnelly, W. Leisenring, T. Siby, I. Thior, M. Dia, E. Gueye, I. Ndoye, and M. Essex. 1992. Prevalence and risk determinants of human immunodeficiency virus type 2 (HIV-2) and human immunodeficiency virus type 1 (HIV-1) in West African female prostitutes. Am. J. Epidemiol. 136:895–907. 15. Kanki, P. J., D. J. Hamel, J. L. Sankale, C. Hsieh, I. Thior, F. Barin, S. A. Woodcock, A. Gueye-Ndiaye, E. Zhang, M. Montano, T. Siby, R. Marlink, I. Noloye, M. E. Essex, and S. Mboup. 1999. Human immunodeficiency virus type 1 subtypes differ in disease progression. J. Infect. Dis. 179:68–73. 16. Leitner, T., S. Kumar, and J. Albert. 1997. Tempo and mode of nucleotide substitutions in gag and env gene fragments in human immunodeficiency
VOL. 78, 2004
17.
18.
19. 20.
21. 22.
23. 24.
25.
26.
27.
virus type 1 populations with a known transmission history. J. Virol. 71:4761– 4770. Lole, K. S., R. C. Bollinger, R. S. Paranjape, D. Gadkari, S. S. Kulkarni, N. G. Novak, R. Ingersoll, H. W. Sheppard, and S. C. Ray. 1999. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J. Virol. 73:152–160. Mboudjeka, I., B. Bikandou, L. Zekeng, J. Takehisa, Y. Harada, Y. Yamaguchi-Kabata, Y. Taniguchi, E. Ido, L. Kaptue, P. M’pelle, H. Parra, M. Ikeda, M. Hayami, and T. Miura. 1999. Genetic diversity of HIV-1 group M from Cameroon and Republic of Congo. Arch. Virol. 144:2291–2311. McCutchan, F. E. 2000. Understanding the genetic diversity of HIV-1. AIDS 14:S31–S44. McCutchan, F. E., J. K. Carr, M. Bajani, E. Sanders-Buell, T. O. Harry, T. C. Stoeckli, K. E. Robbins, W. Gashau, A. Nasidi, W. Janssens, and M. L. Kalish. 1999. Subtype G and multiple forms of A/G intersubtype recombinant human immunodeficiency virus type 1 in Nigeria. Virology 254:226–234. McGrath, K. M., N. G. Hoffman, W. Resch, J. A. E. Nelson, and R. Swanstrom. 2001. Using HIV-1 sequence variability to explore virus biology. Virus Res. 76:137–160. Montavon, C., C. Toure-kane, F. Liegeois, E. Mpoudi, A. Bourgeois, L. Vergne, J. L. Perret, A. Boumah, E. Saman, and S. Mboup. 2000. Most env and gag subtype A HIV-1 viruses circulating in West and West Central Africa are similar to the prototype AG recombinant virus IBNG. J. Acquir. Immune Defic. Syndr. 23:363–374. Osmanov, S., C. Pattou, N. Walker, B. Schwardlander, and J. Esparza. 2002. Estimated global distribution and regional spread of HIV-1 genetic subtypes in the year 2000. J. Acquir. Immune Defic. Syndr. 29:184–190. Peeters, M. 2000. Recombinant HIV sequences: their role in the global epidemic, p. I39–I54. In C. Kuiken, B. Foley, B. Hahn, B. Korber, F. McCutchan, P. Marx, J. Mellors, J. Mullins, J. Sidroski, and J. Wolinsky (ed.), Human Retroviruses and AIDS 2000: a compilation and analysis of nucleic acid and amino acid sequences. Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, N.Mex. Peeters, M., E. Esu-Williams, L. Vergne, C. Montavon, C. Mulanga-Kabeya, T. Harry, A. Ibironke, D. Lesage, D. Patrel, and E. Delaporte. 2000. Predominance of subtype A and G HIV type 1 in Nigeria, with geographical differences in their distribution. AIDS Res. Hum. Retrovir. 16:315–325. Robertson, D. L., J. P. Anderson, J. A. Bradac, J. K. Carr, B. Foley, R. K. Funkhouser, F. Gao, B. H. Hahn, M. L. Kalish, C. Kuiken, G. H. Learn, T. Leitner, F. McCutchan, S. Osmanov, M. Peeters, D. Pieniazek, M. Salminen, P. M. Sharp, S. Wolinsky, and B. Korber. 1999. HIV-1 nomenclature proposal, p. 492–505. In C. L. Kuiken, B. Foley, B. Hahn, B. Korber, F. McCutchan, P. A. Marx, J. W. Mellors, J. I. Mullins, J. Sodroski, and S. Wolinsky (ed.), Human Retroviruses and AIDS 1999. Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, N.Mex. Robertson, D. L., J. P. Anderson, J. A. Bradac, J. K. Carr, B. Foley, R. K. Funkhouser, F. Gao, B. H. Hahn, M. L. Kalish, C. Kuiken, G. H. Learn, T.
HIV-1 SUB-SUBTYPE A3
28.
29.
30.
31.
32.
33. 34. 35.
36.
37.
38.
39.
12445
Leitner, F. McCutchan, S. Osmanov, M. Peeters, D. Pieniazek, M. Salminen, P. M. Sharp, S. Wolinsky, and B. Korber. 2000. HIV-1 nomenclature proposal. Science 288:55–56. Salminen, M. O., J. K. Carr, D. S. Burke, and F. E. McCutchan. 1995. Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res. Hum. Retrovir. 11:1423–1425. Salminen, M. O., C. Koch, E. Sanders-Buell, P. K. Ehrenberg, N. L. Michael, J. K. Carr, D. S. Burke, and F. E. McCutchan. 1995. Recovery of virtually full-length HIV-1 provirus of diverse subtypes from primary virus cultures using the polymerase chain reaction. Virology 213:80–86. Sankale, J., D. Hamel, A. Woolsey, I. Traore, M. Dia, A. Gueye-Ndiaye, M. Essex, S. Mboup, and P. Kanki. 2000. Molecular evolution of human immunodeficiency virus type 1 subtype A in Senegal: 1988–1997. J. Hum. Virol. 3:157–164. Sankale, J.-L., S. Mboup, M. E. Essex, and P. J. Kanki. 1998. Genetic characterization of viral quasispecies in blood and cervical secretions of HIV-1- and HIV-2-infected women. AIDS Res. Hum. Retrovir. 14:1473– 1481. Sarr, A. D., J. L. Sankale, D. J. Hamel, K. U. Travers, A. Gueye-Ndiaye, M. Essex, S. Mboup, and P. J. Kanki. 2000. Interaction with human immunodeficiency virus (HIV) type 2 predicts HIV type 1 genotype. Virology 268: 402–410. Smith, S. 1993. Genetic Data Environment (GDE) v2.2 Manual. Millipore Imaging Systems, Ann Arbor, Mich. Subbarao, S., and G. Schochetman. 1996. Genetic variability of HIV-1. AIDS 10:S13–S23. Toure-Kane, C., C. Montavon, M. A. Faye, P. M. Gueye, P. S. Sow, I. Ndoye, A. Gaye-Diallo, E. Delaporte, M. Peeters, and S. Mboup. 2000. Identification of all HIV type 1 group M subtypes in Senegal, a country with low and stable seroprevalence. AIDS Res. Hum. Retrovir. 16:603–609. Triques, K., A. Bourgeois, S. Saragosti, N. Vidal, E. Mpoudi-Ngole, N. Nzilambi, C. Apetrei, M. Ekwalanga, E. Delaporte, and M. Peeters. 1999. High diversity of HIV-1 subtype F strains in Central Africa. Virology 259: 99–109. Triques, K., A. Bourgeois, N. Vidal, E. Mpoudi-Ngole, C. Mulanga-Kabeya, N. Nzilambi, N. Torimiro, E. Saman, E. Delaporte, and M. Peeters. 2000. Near-full-length genome sequencing of divergent African HIV type 1 subtype F viruses leads to the identification of a new HIV type 1 subtype designated K. AIDS Res. Hum. Retrovir. 16:139–151. VanDerAuwera, G., W. Janssens, L. Heyndrickx, and G. VanDerGroen. 2001. Reanalysis of full-length HIV type 1 group M subtype K and subsubtype F2 with an MS-DOS bootscanning program. AIDS Res. Hum. Retrovir. 17:185–189. Wolinsky, S. M., C. M. Wike, B. T. M. Korber, C. Hutto, W. P. Parks, L. L. Rosenblum, K. J. Kunstman, M. R. Furtado, and J. L. Munoz. 1992. Selective transmission of human immunodeficiency virus type-1 variants from mothers to infants. Science 255:1134–1137.