This article has Open Peer Review reports available.
A considerable proportion of CRF01_AE strains in China originated from circulating intrasubtype recombinant forms (CIRF)
© Jia et al. 2015
Received: 22 January 2015
Accepted: 5 November 2015
Published: 16 November 2015
In this study, the prevalence of HIV-1 CRF01_AE intrasubtype recombinants in China is estimated and their contributions to the epidemic are explored.
Available HIV-1 complete genomes of CRF01_AE were retrieved from the HIV database. The two alignments were evaluated with RDP3. Recombinants were defined as cases in which the recombination signal was supported by at least 3 methods with P-values of ≤0.05 after Bonferroni correction for multiple comparisons implemented in RDP3. Phylogenetic analysis was performed to further investigate the role of intrasubtype recombinants in epidemics.
Here, 124 out of the 339 sequences from around the world (36.6 %) showed significant evidence of recombination. Here, 84 of these recombinants were from China, accounting for 54.9 % of local total sequences (84 out of 153). The results indicated non-negligible levels of intrasubtype recombination. Subsequent phylogenetic analysis indicated that a considerable proportion of CRF01_AE strains in China originated from circulating intrasubtype recombinant forms. Three large, well-supported intrasubtype recombinants clusters were identified here. Through a survey of risk factors and sampling cities and provinces, cluster I and cluster II were found to be prevalent primarily among men who have sex with men in major northern cities. Cluster III was prevalent among heterosexuals and intravenous drug users in southern and southwestern provinces.
The current work highlighted the remarkable prevalence of intrasubtype recombination within the CRF01_AE epidemic and emphasized the value of intrasubtype recombinants, which came to circulate in the same manner as intersubtype recombinants.
Recombination is one of the two most important causes of the high prevalence of variation of human immunodeficiency virus type 1 (HIV-1) . Analyses of recombination can reveal a great deal about evolution in general, much like analyses of nucleotide substitution can. The presence of multiple viral variants can produce HIV-1 recombinant strains that contain sequences either from different genetic subtypes (intersubtype recombination) or different viruses within the same subtype (intrasubtype recombination). Intersubtype recombination of HIV-1 occurs frequently and can produce many intersubtype recombinants, currently known to include 72 circulating recombinant forms (CRFs) of HIV-1 (http://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html) and numerous unique recombinant forms (URFs) of HIV-1. Intersubtype recombination is one major contributing cause of HIV-1 variability, and it facilitates the rapid generation of viral variants with high replicative capacity, drug resistance, or different expression of antigenic epitopes (summarized in [1, 2]). HIV-1 intersubtype recombination can be detected relatively easily . However, HIV-1 intrasubtype recombination has been difficult to detect . To date, most previous studies have focused on HIV-1 intersubtype recombination [5–7]. They have provided us with a very deep understanding of the evolutionary, clinical, and biological relevance, but only a few have addressed intrasubtype recombination in HIV-1. For example, two independent researches revealed frequent intrasubtype recombination among HIV-1 circulating in both South African and Tanzania [4, 8]. In this way, the HIV epidemic and other issues related to HIV-1 intrasubtype recombination have been overlooked. Historically, intersubtype recombinants of HIV-1 have been identified using subtype-specific reference sequences. However, due to the lack of such reference sequences, identification of HIV-1 intrasubtype recombinants has been limited to cases in which known sequences of transmitted multiple viral variants could serve as references for recombination analysis . The introduction of recombination detection software RDP3 eliminates the need for reference sequences and allows us to estimate the prevalence of intrasubtype recombination of CRF01_AE in China [8, 9].
The HIV-1 CRF01_AE was initially named subtype E before it was found to be a recombinant virus [6, 10]. Phylogenetic analysis have indicated that CRF01_AE originated in central Africa in the 1970s and spread to Thailand in the 1980s through heterosexual transmission [10, 11]. Later it was confirmed as the first large-scale epidemic of an intersubtype recombinant around the world [6, 12]. In China, CRF01_AE strains were first identified in the early 1990s. They were found among persons at risk of sexual transmission and intravenous drug users (IDUs) in Yunnan and Guangxi, two provinces in southwestern China that border Myanmar and Vietnam [12–14]. According to the latest nationwide molecular epidemiological survey of newly reported cases in 2006, CRF01_AE strains quickly became the most widespread strains of HIV-1 according to geographic and risk group distributions, accounting for 28 % of nationwide HIV infections during the period . Additionally, CRF01_AE have become the principal strain among men who have sex with men (MSM) in China [16, 17]. A recent study demonstrated that seven distinct phylogenetic clusters of CRF01_AE can be identified among CRF01_AE epidemics in China . Molecular clock analysis indicated that multiple genetically distinct lineages were independently introduced to China from Southeast Asia during the mid-to-late 1990s and subsequently spread into different risk groups and geographic regions within China .
Based on the Los Alamos HIV database (http://www.hiv.lanl.gov/content/index) in which many CRF01_AE sequences have been collected, this study first addresses the prevalence of intrasubtype recombinants among near-full-length CRF01_AE genomes at the host population level in China and then evaluates their contributions to viral transmission and epidemics. The results showed significantly high prevalence of intrasubtype recombinants in circulation and indicated that a considerable proportion of CRF01_AE strains in China originated from circulating intrasubtype recombinant forms.
Selection of sequences
Geographic origin of sequences and sampling year
No. of sequences
CENTRAL AFRICAN REPUBLIC (CF)
HONG KONG (HK)
11 (8 sequences have no sampling year)
133 (1 sequences have no sampling year)
UNITED STATES (US)
Identification of recombination
RDP3, a recombination analysis tool for statistical identification and characterization of recombination events in nucleotide sequences, was used to perform the analysis [8, 9]. RDP3 simultaneously uses a range of non-parametric recombination detection methods: RDP, Geneconv , Bootscan [20, 21], Maxchi [22, 23], Chimaera , SiScan , and 3Seq . RDP3 treats every sequence within the analyzed sequence alignment as a potential recombinant. It systematically screens sequence triplets or quartets to identify a recombinant sequence and two specific sequences that could serve as parents. At the same time, RDP3 performs a statistical evaluation of recombination signals [8, 9]. This approach eliminates the need for reference sequences [8, 9, 26]. The main strength of RDP3 is that it simultaneously uses several different methods to both detect and characterize the recombination events within a sequence alignment without any prior user indication of a non-recombinant set of reference sequences [8, 9]. The sequences are set to linear. The highest acceptable P-value is set to 0.05. Within bootscan analysis, step-size was set to 50 and bootstrap replicates were set to 80 due to file size restriction. The other parameters are default RDP3 settings. An HIV-1 sequence was considered to be recombinant when the recombination signal was supported by at least 3 methods with P-values of ≤0.05 after Bonferroni correction for multiple comparisons implemented in RDP3 [8, 9, 27]. The breakpoint positions predicted were manually checked using recombination signal analysis implemented in RDP3. Although “list events detected by 3 methods” was selected, only 9 out of 22 events were supported by only 3 methods in these results. All others were supported by >3 methods. The subsequently identified cluster I (corresponding to event 19) and cluster II (corresponding to event 20) were both supported by 4 methods. Cluster III (corresponding to event 22) was supported by 5 methods.
jModelTest v0.1.1 was used to find the best-fitting model of nucleotide substitution for the dataset 2 . The general time reversible (GTR) model with g-distributed (G) among-site rate heterogeneity and a proportion of invariant sites (I) (the GTR + I + G model) were chosen as the most appropriate model on the basis of the standard Akaike information criterion (AIC) and Bayesian information criterion (BIC). PhyML 3.0 was used to estimate a maximum likelihood phylogenetic tree for near full-length genome (NFLG) sequences in China using the GTR + I + G nucleotide substitution model . Tree topologies were searched heuristically using the subtree pruning and regrafting procedure (SPR). The confidence of each node in phylogenetic trees was determined using the fast-likelihood-based method of aLRT SH-like. The final maximum likelihood tree was visualized using Mega5 . The topology of the tree was further confirmed by another round of phylogenetic analysis using a Beast tool based on Bayesian inference (data not shown).
Prevalence of intrasubtype recombination of CRF01_AE among populations in China
The alignment of local NFLG sequences from 153 subjects from China infected with HIV-1 CRF01_AE (dataset 2) were first analyzed for evidence of intrasubtype recombination. The purpose of this round of intrasubtype recombination analysis under the background of mainland China was to compare its results to those of the subsequent rounds of intrasubtype recombination analysis of CRF01_AE sequences in mainland China under the background of all over the world, i.e., the background of 339 sequences.
Comparison of results of Chinese CRF01_AE strains intrasubtype analysis between different background areas (mainland China vs all over the world)
No. of sequences from China
Sequences inferred as recombinant
Further analysis was performed on the alignment of 339 CRF01_AE sequences around the world (dataset 1). Results suggested that there are 124 recombinant sequences, accounting for 36.6 % of the total. The putative recombinant sequences corresponding to specific recombination event are listed in Additional file 1. The recombinants were found to be associated with 22 separate recombination events. The most shocking thing was that 84 out of the 124 were from China. The results under the background of mainland China indicated that only 16 were recombinant, but when evaluated against the background of the whole world, many more Chinese sequences were proven to be recombinants (84 versus16). The updated prevalence was as high as 54.9 %, reflecting a significantly high prevalence of intrasubtype recombinants in circulation (Table 2). The apparent discrepancy between the 2 analyses indicated that, when placed in a much wider scope, more strains were connected with either minor or major parents, and therefore much more accurate indication of the amount of intrasubtype recombination can be found. Such results also indicate how background areas might best be selected when performing similar analyses.
Phylogenetic and recombination analysis indicated that considerable CRF01_AE strains in China originated from circulating intrasubtype recombinant forms
Among these 124 recombinants, 1 was from Afghanistan, and all the others were from China (84 recombinants), Thailand (14 recombinants), or Vietnam (25 recombinants), which are geographically close. Such clustering in geographic distribution and the fact that distinct CRF01_AE lineages were independently introduced to China from Southeast Asia suggest that these intrasubtype recombinants very likely originated within this region and that there have probably been epidemics caused by intrasubtype recombinants from the same recombination event.
CRF01_AE was confirmed as the first large-scale epidemic of an intersubtype recombinant HIV-1 strain in the world [6, 12]. The current analysis revealed that intrasubtype recombination of CRF01_AE also contributed greatly to the circulating among HIV-1 population in China. In addition to these 3 big recombinant clusters, there were also exits several much smaller recombinant clusters including only 2 or 3 recombinant sequences (indicated by triangles). Unlike large clusters, they underwent intrasubtype recombination but did not seem to have caused the transmission rate to increase.
In cluster IV, only 3 of 37 sequences were detected as recombinants (accession number: JX112841, JX112817, and GQ845126). These three sequences corresponded to different recombinant events. The first 2 recombinants were associated with recombinant event 6 and the third with recombinant event 17. This is markedly different from recombinant clusters I and II, which were traced to a specific intrasubtype recombinant event, respectively. This case indicated that there first was an extensive epidemic caused by the pure subtype. The scattered intrasubtype recombinants were generated during this period of widespread circulation. The current work revealed there to have almost the same characteristics and rules within intrasubtype recombination and within intersubtype recombination. Take these three intrasubtype recombinants for example, they were the same as the unique recombinant forms (URF), in intersubtype recombination.
It is important to point out, besides the common Recombinant Event 22 (with 5 methods supported, the breakpoint positions of which were HXB2 nt 7156–7549), 5 sequences in cluster III also contained secondary recombinations: GU564230, GU564221, AY008714, and JX112810 are also involved in Recombinant Event 11. GQ845125 is also involved in Recombinant Event 10. In addition, the non-clustered, unique putative recombinant sequence of EF036533 was found to contain three types of segment insertion, so it has three Recombination Event Numbers, including Event 1, Event 2, and Event 20. These phenomena suggest that these 6 strains have undergone multiple recombination events.
Through a survey on the risk factor and sampling cities and provinces, the current cluster I and the current cluster II of intrasubtype recombinants were both found primarily among MSM in major northern cities while the current cluster III was prevalent among heterosexuals and IDUs in southern and southwestern provinces (Additional file 4).
Due to the very short distances between strains, intrasubtype recombination is much more difficult to detect with sequence analysis techniques than intersubtype recombination is. Additionally, recombination analysis also always requires reference sequences. However, such sequences are difficult to find for intrasubtype recombination. Therefore, unlike intersubtype recombination, which has attracted a great deal of attention, research on intrasubtype recombination is relatively sparse. RDP3 treats each sequence within the analyzed sequence alignment as a potential recombinant. Then it systematically screens sequence triplets or quartets to identify a recombinant and two specific sequences that could serve as parents during statistical evaluation of recombination signals [8, 9]. This approach can eliminate the need for reference sequences. Moreover, RDP3 also uses a range of different recombination detection methods to both detect and characterize the recombination events that are evident within the sequence alignment [8, 9]. This can increase the sensitivity of the detection process considerably. This produces more accurate results.
Results in this study first demonstrated that intrasubtype recombination can be detected and is frequent among near-full-length HIV-1 subtype CRF01_AE genomes. In the CRF01_AE epidemic in China, as many as 54.9 % of local total sequences (84 out of 153) are indicated as intrasubtype recombinants. This indicates that the epidemic in China merits further study. The present work demonstrated for the first time that, like intersubtype recombinants that have caused epidemics, such as CRF02_AG, CRF07_BC, CRF08_BC, and BF [31–34], intrasubtype recombinants from the same recombination event can also lead to the prevalence of HIV-1 and AIDS. As shown in Fig. 2 and Additional file 4, intrasubtype recombinants from cluster I and cluster II were found primarily among MSM in major northern cities. Intrasubtype recombinants within cluster III were prevalent among heterosexuals and IDUs in southern and southwestern provinces. Such strains are described in the present work as circulating intrasubtype recombinant forms, echoing circulating recombinant forms. This second term is used to describe intersubtype recombinants that have caused epidemic.
Several previous studies of intrasubtype recombination have concentrated on individual patients [35–37]. Others have focused on populations but within a single geographic location . This work is the first to estimate intrasubtype recombination of viral sequences of CRF01_AE at the host population level using NFLG sequences around the world. As the data show, when placed in a much wider scope, the major and/or minor parents of more strains can be found, and this provides a much more accurate amount of intrasubtype recombination.
The prevalence of intrasubtype recombination provides important evidence of frequent intrasubtype dual infections at the level of populations because for detectable recombination, a host must be first infected by two different viral strains and then these strains must infect a common cell, become copackaged, and generate an infectious new progeny.
There still are several limitations of the work. First, the study was based on sequences that have been submitted to HIV database (the most recent sequences were submitted in 2011). In this way, the results may not reflect the most recent state of HIV infection. Second, although RDP3 can provide a much better performance for intrasubtype recombination detection via simultaneously using seven methods, there still exists limitation in sensitivity due to the very short genetic distances within sequences from the same subtype. With the development of high-throughput sequencing and the availability of more advanced recombination detection software, more optimized analysis may be possible. This would provide even more accurate information regarding intrasubtype recombination.
In summary, the current work highlighted the remarkable prevalence of intrasubtype recombination within the CRF01_AE epidemic and emphasized the value of intrasubtype recombinants, which came to circulate in the same manner as intersubtype recombinants. All these results indicated that much more attention should be paid to intrasubtype recombination.
We would like to thank all the participants for their contributions and cooperation.
This study was supported by the National Science and Technology Special Projects on Major Infectious Diseases (Grant No. 2012ZX10001-002).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Onafuwa-Nuga A, Telesnitsky A. The remarkable frequency of human immunodeficiency virus type 1 genetic recombination. Microbiol Mol Biol Rev. 2009;73:451–80. Table of Contents.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomson MM, Delgado E, Herrero I, Villahermosa ML, Vazquez-de Parga E, Cuevas MT, et al. Diversity of mosaic structures and common ancestry of human immunodeficiency virus type 1 BF intersubtype recombinant viruses from Argentina revealed by analysis of near full-length genome sequences. J Gen Virol. 2002;83:107–19.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA. The effect of recombination on the accuracy of phylogeny estimation. J Mol Evol. 2002;54:396–402.View ArticlePubMedGoogle Scholar
- Rousseau CM, Learn GH, Bhattacharya T, Nickle DC, Heckerman D, Chetty S, et al. Extensive intrasubtype recombination in South African human immunodeficiency virus type 1 subtype C infections. J Virol. 2007;81:4492–500.View ArticlePubMedPubMed CentralGoogle Scholar
- Jia L, Li L, Li H, Liu S, Wang X, Bao Z, et al. Recombination pattern reanalysis of some HIV-1 circulating recombination forms suggest the necessity and difficulty of revision. PLoS One. 2014;9, e107349.View ArticlePubMedPubMed CentralGoogle Scholar
- Robertson DL, Sharp PM, McCutchan FE, Hahn BH. Recombination in HIV-1. Nature. 1995;374:124–6.View ArticlePubMedGoogle Scholar
- Salminen MO, Carr JK, Robertson DL, Hegerich P, Gotte D, Koch C, et al. Evolution and probable transmission of intersubtype recombinant human immunodeficiency virus type 1 in a Zambian couple. J Virol. 1997;71:2647–55.PubMedPubMed CentralGoogle Scholar
- Kiwelu IE, Novitsky V, Margolin L, Baca J, Manongi R, Sam N, et al. Frequent Intra-Subtype Recombination among HIV-1 Circulating in Tanzania. PLoS One. 2013;8, e71131.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. 2010;26:2462–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Gao F, Robertson DL, Morrison SG, Hui H, Craig S, Decker J, et al. The heterosexual human immunodeficiency virus type 1 epidemic in Thailand is caused by an intersubtype (A/E) recombinant of African origin. J Virol. 1996;70:7013–29.PubMedPubMed CentralGoogle Scholar
- McCutchan FE, Artenstein AW, Sanders-Buell E, Salminen MO, Carr JK, Mascola JR, et al. Diversity of the envelope glycoprotein among human immunodeficiency virus type 1 isolates of clade E from Asia and Africa. J Virol. 1996;70:3331–8.PubMedPubMed CentralGoogle Scholar
- Liao H, Tee KK, Hase S, Uenishi R, Li XJ, Kusagawa S, et al. Phylodynamic analysis of the dissemination of HIV-1 CRF01_AE in Vietnam. Virology. 2009;391:51–6.View ArticlePubMedGoogle Scholar
- Yu XF, Chen J, Shao Y, Beyrer C, Lai S. Two subtypes of HIV-1 among injection-drug users in southern China. Lancet. 1998;351:1250.View ArticlePubMedGoogle Scholar
- Cheng H, Zhang J, Capizzi J, Young NL, Mastro TD. HIV-1 subtype E in Yunnan, China. Lancet. 1994;344:953–4.View ArticlePubMedGoogle Scholar
- He X, Xing H, Ruan Y, Hong K, Cheng C, Hu Y, et al. A comprehensive mapping of HIV-1 genotypes in various risk groups and regions across China based on a nationwide molecular epidemiologic survey. PLoS One. 2012;7, e47289.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao B, Han X, Dai D, Liu J, Ding H, Xu J, et al. New trends of primary drug resistance among HIV type 1-infected men who have sex with men in Liaoning Province, China. AIDS Res Hum Retroviruses. 2011;27:1047–53.View ArticlePubMedGoogle Scholar
- Zhang X, Li S, Li X, Xu J, Li D, Ruan Y, et al. Characterization of HIV-1 subtypes and viral antiretroviral drug resistance in men who have sex with men in Beijing, China. AIDS. 2007;21 Suppl 8:S59–65.View ArticlePubMedGoogle Scholar
- Feng Y, He X, Hsi JH, Li F, Li X, Wang Q, et al. The rapidly expanding CRF01_AE epidemic in China is driven by multiple lineages of HIV-1 viruses introduced in the 1990s. AIDS. 2013;27:1793–802.View ArticlePubMedPubMed CentralGoogle Scholar
- Padidam M, Sawyer S, Fauquet CM. Possible emergence of new geminiviruses by frequent recombination. Virology. 1999;265:218–25.View ArticlePubMedGoogle Scholar
- Martin DP, Posada D, Crandall KA, Williamson C. A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retroviruses. 2005;21:98–102.View ArticlePubMedGoogle Scholar
- Salminen MO, Carr JK, Burke DS, McCutchan FE. Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res Hum Retroviruses. 1995;11:1423–5.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA. Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci U S A. 2001;98:13757–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith JM. Analyzing the mosaic structure of genes. J Mol Evol. 1992;34:126–9.PubMedGoogle Scholar
- Gibbs MJ, Armstrong JS, Gibbs AJ. Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000;16:573–82.View ArticlePubMedGoogle Scholar
- Boni MF, Posada D, Feldman MW. An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics. 2007;176:1035–47.View ArticlePubMedPubMed CentralGoogle Scholar
- Novitsky V, Wang R, Margolin L, Baca J, Rossenkhan R, Moyo S, et al. Transmission of single and multiple viral variants in primary HIV-1 subtype C infection. PLoS One. 2011;6, e16714.View ArticlePubMedPubMed CentralGoogle Scholar
- Sentandreu V, Jimenez-Hernandez N, Torres-Puente M, Bracho MA, Valero A, Gosalbes MJ, et al. Evidence of recombination in intrapatient populations of hepatitis C virus. PLoS One. 2008;3, e3239.View ArticlePubMedPubMed CentralGoogle Scholar
- Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.View ArticlePubMedGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lihana RW, Ssemwanga D, Abimiku A, Ndembi N. Update on HIV-1 diversity in Africa: a decade in review. AIDS Rev. 2012;14:83–100.PubMedGoogle Scholar
- Lau KA, Wang B, Saksena NK. Emerging trends of HIV epidemiology in Asia. AIDS Rev. 2007;9:218–29.PubMedGoogle Scholar
- De Sa Filho DJ, Sucupira MC, Caseiro MM, Sabino EC, Diaz RS, Janini LM. Identification of two HIV type 1 circulating recombinant forms in Brazil. AIDS Res Hum Retroviruses. 2006;22:1–13.View ArticlePubMedGoogle Scholar
- Aulicino PC, Kopka J, Mangano AM, Rocco C, Iacono M, Bologna R, et al. Circulation of novel HIV type 1 A, B/C, and F subtypes in Argentina. AIDS Res Hum Retroviruses. 2005;21:158–64.View ArticlePubMedGoogle Scholar
- Philpott S, Burger H, Tsoukas C, Foley B, Anastos K, Kitchen C, et al. Human immunodeficiency virus type 1 genomic RNA sequences in the female genital tract and blood: compartmentalization and intrapatient recombination. J Virol. 2005;79:353–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Shriner D, Rodrigo AG, Nickle DC, Mullins JI. Pervasive genomic recombination of HIV-1 in vivo. Genetics. 2004;167:1573–83.View ArticlePubMedPubMed CentralGoogle Scholar
- van Rij RP, Worobey M, Visser JA, Schuitemaker H. Evolution of R5 and X4 human immunodeficiency virus type 1 gag sequences in vivo: evidence for recombination. Virology. 2003;314:451–9.View ArticlePubMedGoogle Scholar