Volume 10, Number 7—July 2004
Recombinant Viruses and Early Global HIV-1 Epidemic
Central Africa was the epicenter of the HIV type 1 (HIV-1) pandemic. Understanding the early epidemic in the Democratic Republic of the Congo, formerly Zaire, could provide insight into how HIV evolved and assist vaccine design and intervention efforts. Using enzyme immunosorbent assay, we tested 3,988 serum samples collected in Kinshasa in the mid-1980s and confirmed seroreactivity by Western blot. Polymerase chain reaction of gag p17, env C2V3C3, and/or gp41; DNA sequencing; and genetic analyses were performed. Gene regions representing all the HIV-1 group M clades and unclassifiable sequences were found. From two or three short gene regions, 37% of the strains represented recombinant viruses, multiple infections, or both, which suggests that if whole genome sequences were available, most of these strains would have mosaic genomes. We propose that the HIV epidemic was established in Central Africa by the early 1980s and that some recombinant viruses most likely seeded the early global epidemic.
Almost 70% of all HIV type 1 infections worldwide are found in sub-Saharan Africa, and Central Africa is the only region where all HIV-1 groups and group M subtypes have been identified (1–3). This broad range of diversity suggested Central Africa as the epicenter of the HIV-1 pandemic (1). Phylogenetic reconstruction has shown that HIV-1 appears in three distinct lineages, groups M, N, and O (4–7), and each is believed to have arisen through separate zoonotic infections with chimpanzee simian immunodeficiency virus strains in Central Africa (8,9). HIV-1 group M viruses are primarily responsible for the current global epidemic, while group O infections are far fewer and generally found in west-central Africa. Rare strains of group N viruses have been identified in Cameroon, also located in western equatorial Africa. The group M viruses can be subdivided into 9 subtypes, A–D, F–H, J, and K, and at least 14 circulating recombinant forms (CRFs) (http://hiv-web.lanl.gov). These phylogenetic clusterings of HIV strains, subtypes, and CRFs are important records of the epidemic histories of HIV, and they highlight the differences in their geographic distributions and areas where they are endemic.
Some of the earliest characterized whole-genome HIV-1 sequences were collected in Zaire, currently called the Democratic Republic of the Congo (DRC), in the late 1970s and early 1980s. During this time, Kinshasa, the capital of Zaire, was the largest urban area in Central Africa, with a population estimated at 2.5 million (10). From 1984 through 1991, the Zaire Department of Public Health conducted a long-term collaborative HIV research and surveillance program, Project SIDA, with the U.S. Department of Health and Human Services and the Institute of Tropical Medicine, Antwerp, Belgium. In 1991, Project SIDA was abruptly terminated because of civil unrest.
A low and stable HIV prevalence in Zaire/DRC from 1976–1997 (11,12) has been documented, despite social and political upheavals. A genetic survey of HIV strains performed in three regions of the DRC in 1997 demonstrated a high diversity of HIV-1 group M subtypes (1). Serum samples, available from Project SIDA, offered us the opportunity to characterize HIV strains collected in Kinshasa, Zaire, in 1984 and 1986 from employees at Mama Yemo Hospital (13,14). The results provide new insights into the dynamics of HIV infections in a low-prevalence area where multiple subtypes cocirculate, early in the global HIV-1 epidemic.
Specimens and Serologic Analyses
A total of 3,988 serum samples from a collection of samples acquired through Project SIDA from a 1984 and 1986 cross-sectional study of HIV infection among hospital employees at Mama Yemo Hospital in Kinshasa, Zaire; these samples were sent to the National Institutes of Health (NIH) and stored at –20°C until shipment to the Centers for Disease Control and Prevention for serologic and genetic analysis. Specimens numbered 30,000 through 32,000 were collected in 1984. Samples numbered 33,000 through 35,585 were collected among the same population in 1986. The specimens tested represent a convenience sampling from the U.S. repository at NIH and did not include all or only the samples in the original study. The volumes of the specimens ranged from approximately 200 μL to 4 mL. All samples were tested in two separate whole viral-lysate enzyme immunosorbent assay (EIA) for antibodies to HIV-1 and HIV-2 (Genetic Systems, Redmond, WA). Samples with insufficient volumes of serum were excluded from the analysis. Serum specimens reactive by EIA were further tested by using an HIV-1/2 Western blot assay (Genelabs Diagnostics, Singapore, version 2.2).
RNA Extraction, RT, and PCR
Samples with sufficient volumes of serum for reverse transcription–polymerase chain reaction (RT-PCR) amplification with three different sets of primers were selected for further study. RNA was extracted from serum by either the NucliSens nucleic acid manual or automated protocols (Organon Teknika, Boxtel, the Netherlands). RT-PCR was performed by using the Promega RT kit (Promega, Madison, WI) following the manufacturer’s protocol. PCR reactions were placed in a GeneAMP 9600 thermocycler (Perkin-Elmer Cetus, Norwalk, CT) for 35–40 cycles as follows: 94°C for 30s, 55°C for 30s, and 72°C for 60s for all but env gp41, which was at 50°C for 30s and 72°C for 60s, with the final extension at 72°C for 5 min. Nested PCR products were electrophoresed in 1.5% agarose gels (Gibco, Grand Island, NY) and visualized with ethidium bromide staining. To avoid contamination, sample processing and pre-PCR set-up were performed in different rooms than post-PCR manipulations. All samples with discordant phylogenies in different gene regions were verified by having a different person reamplify and resequence the gene regions from a second, unopened vial of serum, if available, or from the previously used vial.
For nested PCR amplification of 380-bp env C2V3C3 fragments, two sets were used: JH44F and JH35MR (outside) and JH33F and JH48R (nested) (15), and MK369 and MK616 (outside) and MK650 and CO601 (nested) (15). Primers for PCR amplification of a 475-bp gag p17 fragment were CL1028 and AB1033 (outside) and CL1029 and AB1032 (nested) (16); those used for a 460-bp fragment of env gp41 were GP40F1 and GP41R1 (outside) and GP46F2 and GP48R2 or GP47R2 (nested) (17).
Sequencing PCR Products
The PCR-amplified products were purified with Qiagen PCR purification kits (Qiagen Inc, Chatsworth, CA) and directly sequenced by using the ABI PRISM Big Dye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems/Perkin Elmer, Foster City, CA) and both forward and reverse nested PCR primers. The sequencing reactions were purified by using the DyeEx Spin Kit (Qiagen–USA, Valencia, CA) and resolved on a 377 ABI Automated DNA Cycle Sequencer (Applied Biosystems).
Sequences were edited by using Sequencher Software v.3.1 (Gene Codes Corp., Madison, WI) and aligned with the Se-Al sequence alignment editor v.1.0 (http://evolve.zoo.ox.ac.uk/software.html?id=seal). Sequences from DRC collected in 1997 (1) were truncated to align with our C2V3C3 sequences, and gaps and ambiguous positions were removed, resulting in an alignment of 304 nt sequences; the first 60 nt represented the C2 region. For the outgroup in our tree, a consensus sequence representing the most common nucleotide at each position for each subtype represented in the Zaire/DRC dataset was selected, and then an overall DRC consensus sequence (DRCcons) was inferred. No positions were undefined in any of the consensus sequences, and no unclassified or unresolved (subtype assignment) sequences were included in the consensus sequence calculation. The Modeltest program, v.3.1 (18) was used to test for a statistically justified model of DNA substitution. The model of evolutionary change used in the tree was Transversion Model (TVM) + g (0.8575), where g is the shape parameter of the γ distribution (heterogeneity among sites), and TVM is the model of substitution whereby A↔G = C↔T and the other four rates are unique. To search for the best tree with such a large dataset (280 taxa) using the DRCcons sequence as the outgroup, several relatively fast tree-building algorithms were employed: Weighbor (http://www.t10.lanl.gov/billb/weighbor/index.html), neighbor-joining (PAUP), BioNJ (PAUP), and Fitch (PAUP). Trees derived from these programs were input to the PAUP* (19) tree scores program to compare the generated log likelihood scores. The neighbor-joining tree program was selected for our phylogenetic analysis since it had the best log likelihood score, although it was not substantially better (Shimodaira-Hasegawa test) than the BioNJ or Weighbor trees.
Genetic Distance Analysis
Subtype and CRF designations were determined from the phylogenetic analysis, and sequences grouping in subtypes were input into the MEGA program version 2.1 (http://www.megasoftware.net) to calculate the means and standard errors of the intrasubtype/intraCRF diversity, using the Kimura 2P distance model.
All 3,988 serum samples were tested by EIA, and 209 seroreactive samples were further tested by Western blot: 140 (3.5%) were HIV-1 seropositive, and 69 (1.7%) had indeterminate Western blot patterns. Samples from infected persons who have not yet seroconverted can display indeterminate Western blot patterns. If we were to assume that all 69 persons with indeterminate Western blot results were seroconverting at the time of sample collection, then the upper limit for the frequency of HIV infections within these early samples was still low at 5.2%. No HIV-2 EIA-positive specimens were confirmed by HIV-2 Western blot analysis.
The serum samples had been stored at –20°C for almost 2 decades, a temperature at which RNases are still active. Because we were concerned that the quality of the RNA would be compromised, we attempted to amplify three relatively short gene regions: the p17 region of gag (474 bp) and two env fragments, C2V3C3 (400 bp) and gp41 (460 bp). All serum specimens could not be amplified, most likely because of the poor quality of the samples. A comparison of intrasubtype diversity was determined for our 50 C2V3C3 sequences and for 181 DRC C2V3C3 sequences from 1997 (1, GenBank, http://www.ncbi.nlm.nih.gov/) with known subtype assignments (Table 1). Unclassifiable sequences were omitted from this analysis. We found the intrasubtype variability was already high in the mid-1980s (9.6%–18.7%) yet was significantly greater for each 1997 env subtype with sufficient numbers for comparisons (A, C, D, E, F1, G). No significant change was found in the frequency of each subtype (Table 1).
Figure 1 represents a neighbor-joining tree containing the following classifiable and unclassifiable C2V3C3 sequences: 56 of our early sequences from Kinshasa, 197 from the DRC collected in 1997 (subtype classifications taken from the V3–V5 sequences published in GenBank), and subtype-specific reference sequences. Some sequences that were difficult to classify in trees containing just our 1980s Kinshasa sequences and subtype-specific reference strains (data not shown) were easier to classify in the context of the larger combined Zaire/DRC phylogenetic tree. At least one of our early strains clustered with subtypes A, B/D, C, D, F (F1), G, H, J, and K, and the circulating recombinant form (CRF)01-AE clustered with 1997 CRF01 sequences. Also, some of our unclassifiable sequences clustered with 1997 unclassifiable DRC strains. One major clade contained subtype A sequences; however, at least four other distinct lineages branched independently in the tree, yet also contained 1997 sequences previously identified as subtype A viruses. One of these lineages that contained 1997 sequences previously designated as subtype A or unclassifiable also contained CRF01-AE strains in the apical portion of the cluster. Three additional unique lineages are present in the tree composed of unclassifiable sequences, perhaps representing new subtypes. Unclassifiable 1997 sequences also branched between the F1 and F2 lineages along with one of our 1980s strains, which suggests a continuum of diversity within subtype F strains instead of distinct sub-subtypes.
Even the node connecting the subtype B and D sequences was no longer distinct because of sequences falling basal to the two lineages, which made the subtype B strains appear to be part of a larger subtype D group. The 1997 sequence labeled D.3 in the tree [97DC.KS26 (accession no. AJ404096)] possibly contained a recombinant breakpoint, which caused it to branch very deep along the D lineage. However, our sequence 30884 also fell outside the node connecting subtypes B and D. In fact, this strain significantly clustered within the node connecting subtypes B and D in p17, C2V3C3, and gp41, yet did not significantly cluster with either subtype. A combined Zaire/DRC consensus sequence was constructed and used as an outgroup in this tree. The branch representing this sequence was very short and was located close to the center of the tree.
To determine an estimate of potentially “pure” subtypes, we analyzed 66 of our mid-1980s HIV strains from which we were able to amplify and sequence two or three gene regions: 53 (80.3%) viruses were amplified in all three gene regions (p17, C2V3C3, and gp41). Figure 2 shows the viruses with concordant phylogenies (possibly pure subtypes) in the sequenced gene regions as follows: A (27%), C (3%), D (18%), F1 (2%), G (8%), K (3%), and our unusual B/D virus (2%). The remaining 37% of the strains represented recombinant viruses; 32% of these appeared to be unique recombinants, except where indicated in Figure 2, and 5% were CRF01_AE, a mosaic lineage containing predominantly subtype A and a unique env lineage called E. Not only was this circulating recombinant form present by the mid-1980s, but only env subtypes A, D, and G were more prevalent. Some of the discrepant gene regions possibly indicate the presence of dual/multiple infections, and the sets of specific gene region primers selectively amplified different strains. However, since our hospital worker population does not represent a high-risk group because the prevalence of infections in Kinshasa was low and recombinant viruses are generally rapidly observed following dual infections (20), we henceforth refer to samples with discordant gene regions as recombinant viruses.
Despite a low frequency of HIV infections in Zaire in the mid-1980s compared to the level of infections currently observed in Africa, we found a remarkably high diversity of HIV-1 strains, with gene regions representing all the group M clades as well as unclassifiable regions. Intrasubtype diversity was already high among our samples, indicating that the HIV epidemic in Kinshasa was already mature by the mid-1980s. No significant change was seen in the frequency of group M subtypes between our early Kinshasa strains and a set of DRC samples collected in 1997 (1), despite political and social upheavals from 1986 to 1997.
We hypothesized that if the HIV epidemic originated in central Africa, phylogenetic analysis of our 1984/1986 Kinshasa and 1997 DRC C2V3C3 sequences might provide a more detailed history of HIV-1 evolution. We found a spectrum of diversity both within and between currently recognized HIV-1 group M subtypes and sub-subtypes (Figure 1). Some lineages, thought to represent discrete subtypes (or sub-subtypes) within the current global epidemic, no longer appeared distinct when analyzed in the context of the large diversity within Zaire/DRC sequences, i.e., multiple, distinct C2V3C3 subtype A lineages; new lineages containing unclassifiable strains; a continuum of diversity within and between the F1 and F2 sub-subtypes; and deep, unclassifiable branches. These findings demonstrate that the HIV diversity in the mid-1980s in Kinshasa was far more complex than in strains currently found in other parts of the world.
Another finding was the placement of C2V3C3 sequences of CRF01-AE in the apical branches of a lineage containing subtype A and unclassifiable sequences as basal branches. CRF01-AE is generally believed to be a recombinant virus with gag and pol genes sharing a common ancestry with subtype A, while the env and vpu genes are derived from a currently unknown parental strain, subtype E. However, others hypothesize that the unique subtype E portion of envelope is not a result of recombination (21) but may be because a higher evolutionary rate in the env and vpu genes make it unrecognizable as subtype A. Our data may be the first demonstration of subtype E env sequences grouping within a cluster of subtype A and previously unclassifiable sequences.
A combined Zaire/DRC consensus sequence used as an outgroup was very short and located near the center of the tree, close to the expected location of the hypothetical ancestor of the HIV-1 group M subtypes. This finding, along with the complexity of the combined phylogenetic tree, supports the theory that the group M subtypes evolved from a common ancestor, which may have originated in this part of Africa.
Since less than 1/10 of the genome was sequenced and analyzed in the present study, and recombination can occur anywhere along the genome, 37% recombinant viruses is very likely a significant underestimate of the actual frequency of chimeric strains cocirculating in Kinshasa by the mid-1980s. In fact, most HIV-1 strains may have contained recombinant genomes by that time. Unfortunately, our data do not allow us to discriminate whether this level of recombination through superinfection intensified the fitness of HIV strains at that time or merely contributed to viral diversity.
Although the precise ancestry of HIV-1 is still uncertain, it appears to have a zoonotic origin (8,9). Ecologic factors that would have allowed human exposure to a natural nonhuman primate host carrying a precursor virus to HIV were therefore instrumental in introducing the virus into humans. This fact suggests that zoonotic introductions of HIV may have occurred in isolated populations but went unnoticed as long as the recipients remained isolated (11) (Figure 3A). However, with increasing movements from rural areas to cities, which occurred in sub-Saharan Africa during the 1960s and 1970s, such isolation became, or was becoming, increasingly rare. The number of sub-Saharan African cities with >500,000 inhabitants rapidly increased from 3 in 1960 to 28 by 1980 (22). From 1965 to 1985, the proportion of the total population living in urban areas of central Africa rose from 21% to 35% and in western Africa from 17% to 29%.
Transmission of HIV would have amplified among high-risk populations, such as commercial sex workers, where HIV superinfections generating recombinant viruses could have resulted (Figure 3B). These recombinants could then have been spread to other, lower risk persons, such as men who visit female sex workers, their spouses, other heterosexual partners, and persons at highest risk for sexually transmitted diseases (23). Other transmission routes could also have helped create and spread recombinant viruses. For example, in 1986, studies from Kinshasa showed a strong association between receiving medical injections and HIV seropositivity among HIV-seropositive infants born to seronegative mothers (24), among hospitalized children 2–14 years of age (25), and among healthcare workers (26). Demand for blood transfusions was also high because of malaria-associated anemia, pregnancy-related complications, and sickle-cell anemia (27,28); at this time, blood was generally not screened for HIV infection (29). Therefore, many factors associated with population migration and urbanization, patterns of sexual activity, and medical practices could have played roles in producing and spreading high numbers of recombinant viruses in Kinshasa by the mid-1980s.
When HIV subtypes were initially genetically characterized in the early 1990s, the first identified viruses were assumed to represent pure subtypes, and viruses found afterward were compared to these prototypic strains. However, our data demonstrate that substantial intersubtype recombination had already occurred at the time when HIV-1 viruses were initially classified (Figure 3C). With so many recombinant viruses present at the cusp of the global HIV epidemic, at least some of the recombinant viruses in central Africa were likely classified as pure subtypes after being exported from Africa and establishing regional epidemics in other parts of the world. This realization could affect our current understanding of the range of diversity within the HIV epidemic (http://hiv-web.lanl.gov/content/hiv-db/REVIEWS/PEETERS2000/Peeters.html); of immunogen selection for vaccine design (30,31), especially in light of superinfections involving multiple subtypes (32); and of mathematical models based on molecular clocks and nonrecombinant strains (33–35) for dating the introduction of HIV into humans.
Dr. Kalish is Team Leader of the Viral Evolution and Transmission Laboratory, HIV and Retrovirology Branch, Division of AIDS, STD, and TB Laboratory Research, National Center for HIV, STD, and TB Prevention, Centers for Disease Control and Prevention. Her research interests include the evolution and molecular epidemiology of HIV, the study of its unusual variants, and investigations of atypical forms of HIV transmission of public health importance.
We thank Andrew Leigh Brown for critical review of the manuscript and Kimberly Distel for editorial review.
- Vidal N, Peeters M, Mulanga-Kabeya C, Nzilambi N, Robertson D, Ilunga W, Unprecedented degree of human immunodeficiency virus type 1 (HIV-1) group M genetic diversity in the Democratic Republic of Congo suggests that the HIV-1 pandemic originated in Central Africa. J Virol. 2000;74:10498–507.
- Murphy E, Korber B, Georges-Courbot MC, You B, Pinter A, Cook D, Diversity of V3 region sequences of human immunodeficiency viruses type 1 from the central African Republic. AIDS Res Hum Retroviruses. 1993;9:997–1006.
- Yang C, Dash B, Hanna SL, Frances HS, Nzilambi N, Colebunders RC, Predominance of HIV type 1 subtype G among commercial sex workers from Kinshasa, Democratic Republic of Congo. AIDS Res Hum Retroviruses. 2001;17:361–5.
- de Leys R, Vanderborght B, vanden Haesevelde M, Heyndrickx L, van Geel A, Wauters C, et al. Isolation and partial characterization of an unusual human immunodeficiency retrovirus from two persons of west-central Africa origin. J Virol. 1990;64:1207–16.
- Gurtler LG, Hauser PH, Eberle J, von Brunn A, Knapp S, Zekeng L, A new subtype of human immunodeficiency virus type 1 (MVP-5180) from Cameroon. J Virol. 1994;68:1581–5.
- Simon F, Mauclere P, Roques P, Loussert-Ajaka I, Muller-Trutwin MC, Saragosti S, Identification of a new human immunodeficiency virus type 1 distinct from group M and group O. Nat Med. 1998;4:1032–7.
- Charneau P, Borman AM, Quillent C, Guetard D, Chamaret S, Cohen J, Isolation and envelope sequence of a highly divergent HIV-1 isolate: definition of a new HIV-1 group. Virology. 1994;205:247–53.
- Gao F, Bailes E, Robertson DL, Chen Y, Rodenburg CM, Michael SF, Origin of HIV-1 in chimpanzee Pan troglodytes troglodytes. Nature. 1999;397:436–41.
- Hahn GH, Shaw GM, De Cock KM, Sharp PM. AIDS as a zoonosis: scientific and public health implications. Science. 2000;287:607–14.
- Mann JM, Francis H, Quinn T, Asila PK, Bosenge N, Nzilambi N, Surveillance for AIDS in a central African city, Kinshasa, Zaire. JAMA. 1986;255:3255–9.
- Nzilambi N, De Cock KM, Forthal DN, Francis H, Ryder RW, Malbe I, The prevalence of infection with human immunodeficiency virus over a 10-year period in rural Zaire. N Engl J Med. 1988;318:276–9.
- Mulanga-Kabeya C, Nzilambi N, Edidi B, Minlangu M, Tshimpaka T, Kambembo L, Evidence of stable HIV seroprevalences in selected populations in the Democratic Republic of the Congo. AIDS. 1998;12:905–10.
- Mann JM, Nzilambi N, Piot P, Bosenge N, Kalala M, Francis H, HIV infection and associated risk factors in female prostitutes in Kinshasa, Zaire. AIDS. 1988;2:249–54.
- N’Galy B, Ryder RW, Bila K, Mwandagalirwa G, Colebunders RL, Francis H, Human immunodeficiency virus infection among employees in an African hospital. N Engl J Med. 1988;319:1123–7.
- Ellenberger DL, Pieniazek D, Nkengasong J, Luo C-C, Maurice C, Janini M, Genetic analysis of human immunodeficiency virus in Abidjan, Cote d’Ivoire reveals predominance of HIV-1 subtype A and introduction of subtype G. AIDS Res Hum Retroviruses. 1999;15:3–9.
- Schochetman GS, Subbarao S, Kalish ML. Methods for studying genetic variation of the human immunodeficiency virus (HIV). In: Adolph KA, editor. Viral genome methods. Boca Raton (FL): CRC Press; 1996. p. 25–41.
- Yang C, Pieniazek D, Owen SM, Fridlund C, Nkengasong J, Mastro TD, Detection of phylogenetically diverse human immunodeficiency virus type 1 groups M and O from plasma by using highly sensitive and specific generic primers. J Clin Microbiol. 1999;37:2581–6.
- Posada D, Crandall KA. MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998;14:817–8.
- Swofford DL. PAUP* Phylogenetic Analysis Using Parsimony (*and Other Methods) Version 4. Sunderland (MA): Sinauer Associates; 2003.
- Iglesias-Sanchez MJ, Lopez-Galindez C. Analysis, quantification, and evolutionary consequences of HIV-1 in vitro recombination. Virology. 2002;304:392–402.
- Anderson JP, Rodrigo AG, Learn GH, Madan A, Delahunty C, Coon M, Testing the hypothesis of a recombinant origin of human immunodeficiency virus type 1 subtype E. J Virol. 2000;74:10752–65.
- Oucho JO, Gould WTS. Chapter 7: internal migration, urbanization, and population distribution. In: Foote KA, Hill KH, Martin LG, editors. Demographic change in sub-Saharan Africa. Washington (DC): National Academy Press; 1993. p. 256–96.
- Quinn TC. The global burden of the HIV pandemic. Lancet. 1996;348:99–106.
- Mann JM, Francis H, Davachi F, Baudoux P, Quinn TC, Nzilambi N, Risk factors for human immunodeficiency virus seropositivity among children 1 to 24 months old in Kinshasa, Zaire. Lancet. 1986;2:654–7.
- Mann JM, Francis H, Davachi F, Baudoux P, Quinn TC, Nzilambi N, Human immunodeficiency virus seroprevalence in pediatric patients 2 to 14 years of age at Mama Yemo Hospital, Kinshasa, Zaire. Pediatrics. 1986;78:673–7.
- Mann JM, Francis H, Quinn TC, Bila K, Asila PK, Bosenge N, HIV seroprevalence among hospital workers in Kinshasa, Zaire. JAMA. 1986;256:3099–102.
- Greenberg AE, Nguyen-Dinh P, Mann JM, Kobote N, Colebunders RL, Francis H, The association between malaria, blood transfusions, and HIV seropositivity in a pediatric population in Kinshasa, Zaire. JAMA. 1988;259:545–9.
- Jager H, N’Galy B, Perriens J, Nseka K, Davachi F, Kabeya CM, Prevention of transfusion-associated HIV transmission in Kinshasa, Zaire: HIV screening is not enough. AIDS. 1990;4:571–4.
- N’tita I, Mulanga K, Dulat C, Lusamba D, Rehle T, Korte R, Risk of transfusion-associated HIV transmission in Kinshasa, Zaire. AIDS. 1991;5:437–9.
- Gaschen B, Taylor J, Yusim K, Foley B, Gao F, Lang D, Diversity considerations in HIV-1 vaccine selection. Science. 2002;296:2354–60.
- Nickle DC, Jensen MA, Gottlieb GS, Shriner D, Learn GH, Rodrigo AG, Consensus and ancestral state HIV vaccines. Science. 2003;299:1515–8.
- Ramos A, Hu DJ, Nguyen L, Phan KO, Vanichseni S, Promadej N, Intersubtype human immunodeficiency virus type 1 superinfection following seroconversion to primary infection in two injection drug users. J Virol. 2002;76:7444–52.
- Schierup MH, Hein J. Consequences of recombination on traditional phylogenetic analysis. Genetics. 2000;156:879–91.
- Wain-Hobson S, Renoux-Elbe C, Vartanian J-P, Meyerhans A. Network analysis of human and simian immunodeficiency virus sequence sets reveals massive recombination resulting in shorter pathways. J Gen Virol. 2003;84:885–95.
- Korber B, Muldoon M, Theiler J, Gao F, Gupta R, Lapedes A, Timing the ancestor of the HIV-1 pandemic strains. Science. 2000;288:1789–96.
Suggested citation for this article: Kalish ML, Robbins KE, Pieniazek D, Schaefer A, Nzilambi N, Quinn TC, et al. Recombinant viruses and early global HIV-1 epidemic. Emerg Infect Dis [serial on the Internet]. 2004 Jul [date cited]. Available from: http://wwwnc.cdc.gov/eid/article/10/7/03-0904.htm