The genetic diversity of HIV-1 is very broad in Cameroon where all group M clades and several circulating recombinant forms CRFs (in particular CRF01_AE, CRF02_AG, CRF06_cpx, CRF09_cpx, CRF11_cpx, CRF13_cpx, CRF22_01A1, CRF36_cpx, and CRF37_cpx), in addition to groups O, N, and P viruses, have been identified.1–10 Among the broad genetic diversity of HIV-1 strains, CRF02_AG represents more than 65% of HIV infections with an additional 26% classified as unique recombinant forms (URFs).8,10 Cocirculation of different subtypes and CRFs in Cameroon results in the continuing emergence of new intersubtype recombinants or even M/O recombinants.11 Therefore, Cameroon is an ideal area to investigate the genetic diversity of HIV-1 and its impact on the global pandemic. In 2002, we initiated a molecular and serologic epidemiology survey in blood donors in Douala and Yaoundé, in Cameroon to investigate HIV-1 genetic diversity, study virologic and immunologic characteristics of the viruses, evaluate the performance of US Food and Drug Administration licensed HIV-1 assays in the presence of numerous HIV-1 variants, and identify samples to serve as candidate reference reagents for diagnostics. We found that several group M subtypes (A, B, C, D, F2, and G) and CRFs (CRF02_AG, CRF06_cpx, CRF11_cpx, CRF13_cpx, CRF19_cpx, and CRF22_01A1) were circulating in this population. A new URF strain, 02CAMLT04, had CRF22_01A1gag-CRF02_AGenv genotype,12 and a novel strain, 02CAMLT72 isolated in Douala, was designated as pure CRF22_01A1.13
CRF22_01A1 strain was initially identified in Cameroon in 20014 and circulating in this country for many years. The studies of Brennan et al14 reported that the percentage of CRF22_01A1 in concordant specimens was 6.6% (1996–1999), 4.7% (2000–2002), and 6.6% (2003–2004) over the 9-year period, however, molecular, epidemiological, and evolutionary study of CRF22_01A1 phylogenetic association with different HIV-1 subtypes was not described. Here, we report 14 full-length genomic sequences of HIV-1 specimens isolated in Cameroon from 2002 to 2003 and 2006 to 2010. Of these 14 specimens, 5 were pure CRF22_01A1, 6 were determined to be recombinants of CRF22_01A1 with CRF02_AG. These new data support involvement of CRF22_01A1 in the emerging diversity of HIV-1 in Cameroon.
For HIV characterization studies, 441 blood samples from HIV-1–seropositive individuals were collected at Douala (LT), Bamenda (LPH), Buea (BDHS), Limbe (LB), Yaounde (ARC), and a few villages in Cameroon (NYU) in 2002–2010. A total of 208 viruses were isolated and genotyped based on partial sequence analysis, 63% of them were CRF02_AG, 26% were URFs, and 6% were CRF22_01A1-containing recombinants.12,15 In this study, 13 viruses containing CRF22_01A1 genome and 3 URF strains were selected for full-length sequence analysis, however, nucleic acid was available for only 14 of them which were completely sequenced and characterized.
RNA Extraction, Reverse Transcriptase—Polymerase Chain Reaction, and Nearly Full-Length Genomic Sequencing
Viral RNA was extracted from the plasma samples using QIAamp Viral RNA Mini kit (Qiagen, Valencia, CA). Then, reverse transcription, polymerase chain reaction, and sequencing were performed as previously described.6,7 The nearly full-length genomic sequences reported in this study have been deposited in the GenBank database (accession numbers: EU743964, JN864047-JN864059).
Phylogenetic and Recombination Sequence Analysis
Nucleotide sequences were aligned with HIV-1 reference strains of different subtypes, subsubtypes, and CRFs from Los Alamos HIV database (A1_SE8538, A1_PS1044, A1_92UG037, A1_Q23 17, A1SE7535, A1_92RW008; A2_97CDKTB48, A2_94CY017; B_HXB2 LAI IIIB, B_671_00T36; C_BR025, C_ETH2220; D_ELI, D_94UG114; F1_VI850, F1_93BR020_1; F2_MP255, F2_02CM.0016BBY; G_DRCBL, G_92NG083; H_VI997, H_VI991; J_SE7887, J_SE7022; K_EQTB11C, K_MP535; 01_90CF4071, 01_90CF11697, 01_93TH051, 01_93TH253, 01_CM240; 02_pBD6_15, 02_IBNG)16 using Clustal W method17 as implemented in Vector NTI software (Invitrogen, Carlsbad, CA). The phylogenetic reconstructions were also performed with MEGA 4 software package using the neighbor-joining method.18 To analyze the recombinant structure of the new viruses, the SimPlot 3.5.1 software was used to determine the percentage of similarity between selected pairs of sequences and to calculate bootscan plots.19 The regions that did not cluster with any of the known subtypes were submitted to BLAST analysis (http://hiv-web.lanl.gov/BASIC-BLAST) to find the closely related sequences of other HIV-1 strains. The viral nucleotide positions correspond to the reference strain HXB2 (GenBank Accession Number K03455) used throughout this article.
Near full-length genome sequencing was completed for 14 viruses of interest. One virus, NYU488, was collected from Bapile located in the East Province of Cameroon, and 13 viruses were collected from the Southwest Province of Cameroon (Fig. 1A). Phylogenetic analysis of these sequences confirmed their relatedness, as they all clustered together in the CRF22_01A1, CRF02_AG, or subtype G radiation with a high bootstrap value (BV). As shown in Figure 2, 6 strains (02CAMLT04, LPH27MF, LB005, LB011, LB013, and LB054) were deeply clustered in the CRF22_01A1 radiation (BV, 100%) revealing a close relationship, 3 strains (LB045, NYU488, and ARC087) clustered between CRF22_01A1 and CRF01_AE suggesting they are CRF22_01A1 recombinants. Four strains (ARC007, LB052, BDSH129, and LT66) clustered with CRF02_AG (BV, 93%) and LT31 clustered with subtype G (BV, 100%). These sequences were further analyzed by bootscanning to determine the phylogenetic mosaic structure. The 02CAMLT04 sequence was plotted against references restricted to the subsubtype A1, F2, subtype B, C, H, K, CRF22_01A1, and CRF02_AG. The genome of 02CAMLT04 was divided into 3 segments at the 2 breakpoints of nt 6951 and nt 9002, and the neighbor-joining trees were further built for each segment (Fig. 3). The first segment extended from the gag to the beginning of the env gene (nt 789-6951), and clustered with CRF22_01A1 (BV, 100%). The second segment representing most of the env gene (nt 6952-9002) clustered with CRF02_AG radiation (BV, 100%) and the rest of the segments covering the nef gene (nt 9003-9408) clustered with CRF22_01A1. A significant portion (76%) of the 02CAMLT04 genome was in perfect alignment with CRF22_01A1 reference sequences except that most of the env gene clustered with CRF02_AG. These results demonstrated that the 02CAMLT04 virus was generated by the recombination of CRF22_01A1 and CRF02_AG. To our knowledge, this is the first report of a novel URF of CRF22_01A1 recombinant with CRF02_AG in full-length sequence.
Bootscan analysis of these viruses revealed the presence of CRF22_01A1 homology, the overall structures of these CRF22_01A1 containing viruses are diagrammed in Figure 4. Five sequences (LPH27MF, LB005, LB011, LB013, and LB054) clustered phylogenetically with CRF22_01A1 references throughout the entire genome demonstrating that these viruses share the same genomic mosaic and can be designated as pure CRF22_01A1. Two viruses (ARC087 and NYU488) shared identical recombinant breakpoints and were largely CRF22_01A1 (80% of whole genome spanning mostly pol to nef genes) and 20% for CRF02_AG. LB045 virus was divided into 4 segments with a complex CRF02/CRF22/CRF02/CRF22 recombinant, about 65% of the LB045 genome was CRF22_01A1. BDSH129 virus was almost entirely CRF02_AG (84%) with a small piece of CRF22_01A1 (16%) spanning the gag gene and the beginning of the pol gene. LB052, a triple CRF02/CRF22/CRF02 recombinant, was also predominantly CRF02_AG (76%) with a small portion of CRF22_01A1 fragment spanning the most of env gene from the beginning. ARC007 virus was a subtype B/CRF02 recombinant strain. LT66 and LT31 viruses were a complex F2/CRF02/F2/CRF02/F2 and CRF01/G/CRF02/G recombinant, respectively (data not shown).
To further understand the involvement of CRF22_01A1 in genetic recombination, the blast search of HIV sequence database against 6162 nt of 02CAMLT04 segment (nt 789-6951) yielded 4 high search score HIV-1 strains 01CM.0130NY, 01CM.1152NG, 02CM.3163MN, and 01CM.4008HAN which were originally identified in Cameroon and classified as 01AF2U, A1U, AF2, and 01DU, respectively.20 Phylogenetic analysis showed that 01CM.0130NY virus clustered with high confidence to CRF22_01A1 radiation (BV, 100%, Fig. 2), 02CM.3163MN and 01CN.1152NG clustered with other CRF22_01A1-containing recombinant (BV, 85%). 01CM.4008HAN was a unique strain clustered between CRF01_AE and subsubtype A1. Further bootscanning analysis showed that they were recombinants of CRF22_01A1 with HIV-1 subtypes F, D, subsubtype A1, F2, or CRFs, all of these strains had a unique mosaic pattern combining CRF22_01A1 (Fig. 4). The detailed phylogenetic composition of CRF22_01A1-containing recombinants described above is listed in Supplemental Digital Content 1 (see Table 1, http://links.lww.com/QAI/A315).
We identified 5 pure CRF22_01A1 and 6 CRF22_01A1-containing URFs strains in Cameroon. We also reclassified 4 CRF22_01A1-containing recombinants from the HIV Sequence GenBank (Fig. 4). Nine unique recombination patterns were identified; breakpoints were found throughout the genome, the most definitive breakpoints were at the beginning of the pol, accessory or beginning of the env region. CRF22_01A1 fragments may be identified in any region across the HIV whole genome. Table 1 summarizes the subtype/CRF assignments of recombinants containing CRF22_01A1 fragments evaluated in samples collected from 2000 to 2010. Since the identification of CRF22_01A1 in Cameroon, it may be emerging as a dominant HIV-1 variant in this population that could act as a parental subtype for recombination events. The emergence of more complex intersubtype recombinants containing fragments of CRF22_01A1, such as CRF36_cpx, have been reported in Cameroon.7 Reanalysis of the CRF36 genome including CRF22_01A1 reference sequences showed that the CRF36_cpx contains the fragment of CRF22_01A1, CRF01_AE, and CRF02_AG. Brennan et al14 found that 22.4% of HIV-1 URFs were CRF22_01A1 in the blood donor samples collected in 1996–2004 in Cameroon. Our recent study of samples collected from 2006 to 2008 in Cameroon also identified 3 CRF02gag/CRF22pol/CRF22env strains.12,13 These studies highlight the complexity of recombinants between CRF22_01A1 and CRF02_AG in Cameroon (Table 1). CRF02_AG appeared to be the dominant virus in the population and accounted for 60%–68% of HIV-1 infections with an additional 26% classified as URFs.2,21–23 This high degree of recombination could be due to preferential founder effects, and the subsequent introduction of the other genetic form, such as CRF22_01A1 or CRF02_AG, frequently resulting in recombinants containing fragments of CRF22_01A1 or CRF02_AG.
One interesting finding in this study is that all of the recombinants of HIV-1 characterized displayed a fairly similar breakpoint at the accessory gene and env region, which spans from vif to gp120 gene. For example, a breakpoint of 02CAMLT04 genome in nucleotide position nt 6911 in gp120 was also observed around similar position in 02CAMLT72 (nt 7015), 01CM.1152NG (nt 6963), LB052 (nt 6099), LB045 (nt 6261), 01CM.4008HAN (nt 6088), and 01CM.0130NY (nt 6668). The segment from one-third of the env to the 5´ end of the nef region (nt 7016-8895) in 02CAMLT72 showed a CRF22_01A1 virus derived from subsubtype A1 but was replaced by CRF02_AG in a similar region of 02CAMLT04 virus, suggesting that this region may be a hotspot for recombination between CRF22_01A1 and other subtypes/CRFs. The exact prevalence of the CRF22_01A1 intersubtype recombinants is still unclear. However, according to the current study, given that the CRF22_01A1-containing recombinants were possibly introduced to the population as early as 2000, CRF22_01A1 could be more broadly distributed in Cameroon even outside this country (Fig. 1B). Although the CRF22_01A1 strain was prevalent over the 10-year period, the numbers of reported CRF22_01A1 cases were limited since reference strains of CRF22_01A1 were not well classified until now. The complexity of CRF22_01A1 strains in Cameroon seemed to increase during the 2002–2010 period. The CRF22_01A1 strain has also been identified in Saudi Arabia,24,25 Equatorial Guinea,26,27 Central African patients,28 and United States29 indicating that the CRF22 strain was spreading to different geographic regions.
It has been well known that HIV-1 CRFs behave like a pure subtype of HIV-1 and are able to recombine with other subtypes or CRFs during HIV-1 viral reverse transcription when co-infection or super-infection occurs. We document here that CRF22_01A1 is circulating in Cameroon and is involved in complex intersubtype recombination events with subtypes (D or F), subsubtypes (A1 or F2), and CRFs (CRF01_AE or CRF02_AG). CRF22_01A1 recombinants seem to cocirculate at higher proportions than its pure prototype. Our findings further demonstrate the dynamic evolution of emerging variants in Cameroon, which could potentially impact the phlyogenetic nature of the epidemic in this region in the future.
The authors wish to acknowledge Drs Mingjie Zhang, Krishnakumar Devadas, and Robin Biswas for review of the article. The findings and conclusions in this article have not been formally disseminated by the Food and Drug Administration and should not be construed to represent any Agency determination or policy.
1. Gao F, Bailes E, Robertson DL, et al.. Origin of HIV-1 in the chimpanzee Pan troglodytes troglodytes. Nature. 1999;397:436–441.
2. Hahn BH, Shaw GM, De Cock KM, et al.. AIDS as a zoonosis: scientific and public health implications. Science. 2000;287:607–614.
3. Keele BF, Van Heuverswyn F, Li Y, et al.. Chimpanzee reservoirs of pandemic and nonpandemic HIV-1. Science. 2006;313:523–526.
4. Carr JK, Torimiro JN, Wolfe ND, et al.. The AG recombinant IbNG and novel strains of group M HIV-1 are common in Cameroon. Virology. 2001;286:168–181.
5. Plantier JC, Leoz M, Dickerson JE, et al.. A new human immunodeficiency virus derived from gorillas. Nat Med. 2009;15:871–872.
6. Powell RL, Zhao J, Konings FA, et al.. Circulating recombinant form (CRF) 37_cpx: an old strain in Cameroon composed of diverse, genetically distant lineages of subtypes A and G. AIDS Res Hum Retroviruses. 2007;23:923–933.
7. Powell RL, Zhao J, Konings FA, et al.. Identification of a novel circulating recombinant form (CRF) 36_cpx in Cameroon that combines two CRFs (01_AE and 02_AG) with ancestral lineages of subtypes A and G. AIDS Res Hum Retroviruses. 2007;23:1008–1019.
8. Vergne L, Bourgeois A, Mpoudi-Ngole E, et al.. Biological and genetic characteristics of HIV infections in Cameroon reveals dual group M and O infections and a correlation between SI-inducing phenotype of the predominant CRF02_AG variant and disease stage. Virology. 2003;310:254–266.
9. Wilbe K, Casper C, Albert J, et al.. Identification of two CRF11-cpx genomes and two preliminary representatives of a new circulating recombinant form (CRF13-cpx) of HIV type 1 in Cameroon. AIDS Res Hum Retroviruses. 2002;18:849–856.
10. Zhong P, Burda S, Urbanski M, et al.. HIV type 1 group M clades infecting subjects from rural villages in equatorial rain forests of Cameroon. J Acquir Immune Defic Syndr. 2002;31:495–505.
11. Yamaguchi J, Bodelle P, Vallari AS, et al.. HIV infections in northwestern Cameroon: identification of HIV type 1 group O and dual HIV type 1 group M and group O infections. AIDS Res Hum Retroviruses. 2004;20:944–957.
12. Machuca A, Tang S, Hu J, et al.. Increased genetic diversity and intersubtype recombinants of HIV-1 in blood donors from urban Cameroon. J Acquir Immune Defic Syndr. 2007;45:361–363.
13. Zhao J, Tang S, Ragupathy V, et al.. Identification and genetic characterization of a novel CRF22_01A1 recombinant form of HIV type 1 in Cameroon. AIDS Res Hum Retroviruses. 2010;26:1033–1045.
14. Brennan C, Bodelle P, Coffey R, et al.. The prevalence of diverse HIV-1 strains was stable in Cameroonian blood donors from 1996 to 2004. J Acquir Immune Defic Syndr. 2008;49:432–439.
15. Ragupathy V, Zhao J, Wood O, et al.. Identification of new, emerging HIV-1 unique recombinant forms and drug resistant viruses circulating in Cameroon. Virol J. 2011;8:185.
17. Thompson JD, Gibson TJ, Plewniak F, et al.. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25:4876–4882.
18. Tamura K, Dudley J, Nei M, et al.. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–1599.
19. Lole KS, Bollinger RC, Paranjape RS, et al.. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73:152–160.
20. Kijak GH, Sanders-Buell E, Wolfe ND, et al.. Development and application of a high-throughput HIV type 1 genotyping assay to identify CRF02_AG in West/West Central Africa. AIDS Res Hum Retroviruses. 2004;20:521–530.
21. McCutchan FE. Understanding the genetic diversity of HIV-1. AIDS. 2000;14(suppl 3):S31–S44.
22. Robertson DL, Anderson JP, Bradac JA, et al.. HIV-1 nomenclature proposal. Science. 2000;288:55–56.
23. Swanson P, de Mendoza C, Joshi Y, et al.. Impact of human immunodeficiency virus type 1 (HIV-1) genetic diversity on performance of four commercial viral load assays: LCx HIV RNA quantitative, AMPLICOR HIV-1 MONITOR v1.5, VERSANT HIV-1 RNA 3.0, and NucliSens HIV-1 QT. J Clin Microbiol. 2005;43:3860–3868.
24. Yamaguchi J, Badreddine S, Swanson P, et al.. Identification of new CRF43_02G and CRF25_cpx in Saudi Arabia based on full genome sequence analysis of six HIV type 1 isolates. AIDS Res Hum Retroviruses. 2008;24:1327–1335.
25. Badreddine S, Smith K, van Zyl H, et al.. Identification and characterization of HIV type 1 subtypes present in the Kingdom of Saudi Arabia: high level of genetic diversity found. AIDS Res Hum Retroviruses. 2007;23:667–674.
26. Yebra G, Rivas P, Herrero MD, et al.. Clinical differences and viral diversity between newly HIV type 1-diagnosed African and non-African patients in Spain (2005–2007). AIDS Res Hum Retroviruses. 2009;25:37–44.
27. Djoko CF, Wolfe ND, Vidal N, et al.. HIV type 1 pol gene diversity and genotypic antiretroviral drug resistance mutations in Malabo, Equatorial Guinea. AIDS Res Hum Retroviruses. 2010;26:1027–1031.
28. Moussa S, Pinson P, Pelembi P, et al.. First data on HIV-1 resistance mutations to antiretroviral drugs in Central African Republic. AIDS Res Hum Retroviruses. 2010;26:1247–1248.
29. Pyne MT, Holzmayer V, John Hackett J, et al.. Analysis of HIV subtype prevalence and geographic distribution in a large U.S. population. Presented at: 16th Conference on Retroviruses & Opportunistic Infections; February 8–11, 2009; Montreal, Canada. 2009;Poster #292.