Identification of 3 phylogenetically related HIV-1 BG intersubtype circulating recombinant forms in Cuba

Identification of CRF89_BF, a new member of an HIV-1 circulating BF intersubtype recombinant form family widely spread in South America

Circulating recombinant forms (CRFs) contribute substantially to the HIV-1 pandemic. Among 105 CRFs described in the literature, 16 are BF intersubtype recombinants, most of South American origin, of which CRF12_BF is the most widely spread. A BF recombinant cluster identified in Bolivia was suggested to represent a new CRF_BF. Here we find that it belongs to a larger cluster incorporating 39 viruses collected in 7 countries from 3 continents, 22 of them in Spain, most from Bolivian or Peruvian individuals, and 12 in South America (Bolivia, Argentina, and Peru). This BF cluster comprises three major subclusters, two associated with Bolivian and one with Peruvian individuals. Near full-length genome sequence analyses of nine viruses, collected in Spain, Bolivia, and Peru, revealed coincident BF mosaic structures, with 13 breakpoints, 6 and 7 of which coincided with CRF12_BF and CRF17_BF, respectively. In a phylogenetic tree, they grouped in a clade closely related to these CRFs, and more distantly to CRF38_BF and CRF44_BF, all circulating in South America. These results allowed to identify a new HIV-1 CRF, designated CRF89_BF. Through phylodynamic analyses, CRF89_BF emergence was estimated in Bolivia around 1986. CRF89_BF is the fifth CRF member of the HIV-1 recombinant family related to CRF12_BF.

Characterization of New Recombinant Forms of HIV-1 From the Comunitat Valenciana (Spain) by Phylogenetic Incongruence

Recombination is one of the main processes shaping the evolution of HIV-1, with relevant consequences for its epidemiology. In fact, Circulating and Unique Recombinant Forms (CRFs and URFs) cause 23% of current infections. The routine analyses of antiretroviral resistance yield partial pol gene sequences that can be exploited for molecular epidemiology surveillance but also to study viral diversity and to detect potential recombinant samples. Among the pol sequences derived from a large sample dataset from the Comunitat Valenciana (Spain), we identified nine putative recombinant samples. We aimed at fully characterizing these samples and performing a detailed analysis of the corresponding recombination events. We obtained nearly full-genome sequences and used jpHMM and RDP4 to detect and characterize recombinant fragments. We assessed the confidence of these inferences by likelihood mapping and phylogenetic placement with topology congruence tests. Next, we performed a phylogenetic analysis of each putative recombinant fragment to determine its relationships to previously described recombinant forms. We found that two samples related to CRF44_BF whereas the rest corresponded to new URFs (two URF_AD, one URF_BG that can constitute a new CRF resulting from subtype B and CRF24_BG, and two URF_cpx composed of A, G, K, H, and J subtypes). These URFs have a complex recombination pattern that cannot be determined accurately. They seem to have arisen by successive recombination events among lineages, including other CRFs. Our results highlight the usefulness of routine surveillance analysis for the detection of new HIV-1 recombination forms and, at the same time, the need for full-genome sequencing and recombination detection guidelines to properly characterize this complex process.

Genetic Characterization of a Panel of Diverse HIV-1 Isolates at Seven International Sites

HIV-1 subtypes and drug resistance are routinely tested by many international surveillance groups. However, results from different sites often vary. A systematic comparison of results from multiple sites is needed to determine whether a standardized protocol is required for consistent and accurate data analysis. A panel of well-characterized HIV-1 isolates (N = 50) from the External Quality Assurance Program Oversight Laboratory (EQAPOL) was assembled for evaluation at seven international sites. This virus panel included seven subtypes, six circulating recombinant forms (CRFs), nine unique recombinant forms (URFs) and three group O viruses. Seven viruses contained 10 major drug resistance mutations (DRMs). HIV-1 isolates were prepared at a concentration of 107 copies/ml and compiled into blinded panels. Subtypes and DRMs were determined with partial or full pol gene sequences by conventional Sanger sequencing and/or Next Generation Sequencing (NGS). Subtype and DRM results were reported and decoded for comparison with full-length genome sequences generated by EQAPOL. The partial pol gene was amplified by RT-PCR and sequenced for 89.4%-100% of group M viruses at six sites. Subtyping results of majority of the viruses (83%-97.9%) were correctly determined for the partial pol sequences. All 10 major DRMs in seven isolates were detected at these six sites. The complete pol gene sequence was also obtained by NGS at one site. However, this method missed six group M viruses and sequences contained host chromosome fragments. Three group O viruses were only characterized with additional group O-specific RT-PCR primers employed by one site. These results indicate that PCR protocols and subtyping tools should be standardized to efficiently amplify diverse viruses and more consistently assign virus genotypes, which is critical for accurate global subtype and drug resistance surveillance. Targeted NGS analysis of partial pol sequences can serve as an alternative approach, especially for detection of low-abundance DRMs.

High Degree of HIV-1 Group M (HIV-1M) Genetic Diversity within Circulating Recombinant Forms: Insight into the Early Events of HIV-1M Evolution

The existence of various highly divergent HIV-1 lineages and of recombination-derived sequence tracts of indeterminate origin within established circulating recombinant forms (CRFs) strongly suggests that HIV-1 group M (HIV-1M) diversity is not fully represented under the current classification system. Here we used a fully exploratory screen for recombination on a set of 480 near-full-length genomes representing the full known diversity of HIV-1M. We decomposed recombinant sequences into their constituent parts and then used maximum-likelihood phylogenetic analyses of this mostly recombination-free data set to identify rare divergent sequence lineages that fall outside the major named HIV-1M taxonomic groupings. We found that many of the sequence fragments occurring within CRFs (including CRF04_cpx, CRF06_cpx, CRF11_cpx, CRF18_cpx, CRF25_cpx, CRF27_cpx, and CRF49_cpx) are in fact likely derived from divergent unclassified parental lineages that may predate the current subtypes, even though they are presently identified as derived from currently defined HIV-1M subtypes. Our evidence suggests that some of these CRFs are descended predominantly from what were or are major previously unidentified HIV-1M lineages that were likely epidemiologically relevant during the early stages of the HIV-1M epidemic. The restriction of these divergent lineages to the Congo basin suggests that they were less infectious and/or simply not present at the time and place of the initial migratory wave that triggered the global epidemic.IMPORTANCE HIV-1 group M (HIV-1M) likely spread to the rest of the world from the Congo basin in the mid-1900s (N. R. Faria et al., Science 346:56-61, 2014, http://dx.doi.org/10.1126/science.1256739) and is today the principal cause of the AIDS pandemic. Here, we show that large sequence fragments from several HIV-1M circulating recombinant forms (CRFs) are derived from divergent parental lineages that cannot reasonably be classified within the nine established HIV-1M subtypes. These lineages are likely to have been epidemiologically relevant in the Congo basin at the onset of the epidemic. Nonetheless, they appear not to have undergone the same explosive global spread as other HIV-1M subtypes, perhaps because they were less transmissible. Concerted efforts to characterize more of these divergent lineages could allow the accurate inference and chemical synthesis of epidemiologically key ancestral HIV-1M variants so as to directly test competing hypotheses relating to the viral genetic factors that enabled the present pandemic.

CRF19_cpx is an Evolutionary fit HIV-1 Variant Strongly Associated With Rapid Progression to AIDS in Cuba

Background

Clinicians reported an increasing trend of rapid progression (RP) (AIDS within 3 years of infection) in Cuba.

Methods

Recently infected patients were prospectively sampled, 52 RP at AIDS diagnosis (AIDS-RP) and 21 without AIDS in the same time frame (non-AIDS). 22 patients were sampled at AIDS diagnosis (chronic-AIDS) retrospectively assessed as > 3 years infected. Clinical, demographic, virological, epidemiological and immunological data were collected. Pol and env sequences were used for subtyping, transmission cluster analysis, and prediction of resistance, co-receptor use and evolutionary fitness. Host, immunological and viral predictors of RP were explored through data mining.

Findings

Subtyping revealed 26 subtype B strains, 6 C, 6 CRF18_cpx, 9 CRF19_cpx, 29 BG-recombinants and other subtypes/URFs. All patients infected with CRF19 belonged to the AIDS-RP group. Data mining identified CRF19, oral candidiasis and RANTES levels as the strongest predictors of AIDS-RP. CRF19 was more frequently predicted to use the CXCR4 co-receptor, had higher fitness scores in the protease region, and patients had higher viral load at diagnosis.

Interpretation

CRF19 is a recombinant of subtype D (C-part of Gag, PR, RT and nef), subtype A (N-part of Gag, Integrase, Env) and subtype G (Vif, Vpr, Vpu and C-part of Env). Since subtypes D and A have been associated with respectively faster and slower disease progression, our findings might indicate a fit PR driving high viral load, which in combination with co-infections may boost RANTES levels and thus CXCR4 use, potentially explaining the fast progression. We propose that CRF19 is evolutionary very fit and causing rapid progression to AIDS in many newly infected patients in Cuba.

•

We propose that CRF19 is evolutionary very fit, causing rapid progression to AIDS in many newly infected patients in Cuba.

•

CRF19 is a recombinant of subtype D, subtype A and subtype G, with a subtype D protease estimated to be particularly fit.

•

A fit protease with high viral load and co-infections, may boost RANTES levels and thus CXCR4 use, hence fast progression.

Recombination pattern reanalysis of some HIV-1 circulating recombination forms suggest the necessity and difficulty of revision

BACKGROUND

Recombination is one of the major mechanisms underlying the generation of HIV-1 variability. Currently 61 circulating recombinant forms of HIV-1 have been identified. With the development of recombination detection techniques and accumulation of HIV-1 reference stains, more accurate mosaic structures of circulating recombinant forms (CRFs), like CRF04 and CRF06, have undergone repeated analysis and upgrades. Such revisions may also be necessary for other CRFs. Unlike previous studies, whose results are based primarily on a single recombination detection program, the current study was based on multiple recombination analysis, which may have produced more impartial results.

METHODS

Representative references of 3 categories of intersubtype recombinants were selected, including BC recombinants (CRF07 and CRF08), BG recombinants (CRF23 and CRF24), and BF recombinants (CRF38 and CRF44). They were reanalyzed in detail using both the jumping profile hidden Markov model and RDP3.

RESULTS

The results indicate that revisions and upgrades are very necessary and the entire re-analysis suggested 2 types of revision: (i) length of inserted fragments; and (ii) number of inserted fragments. The reanalysis also indicated that determination of small regions of about 200 bases or fewer should be performed with more caution.

CONCLUSION

Results indicated that the involvement of multiple recombination detection programs is very necessary. Additionally, results suggested two major challenges, one involving the difficulty of accurately determining the locations of breakpoints and the second involving identification of small regions of about 200 bases or fewer with greater caution. Both indicate the complexity of HIV-1 recombination. The resolution would depend critically on development of a recombination analysis algorithm, accumulation of HIV-1 stains, and a higher sequencing quality. With the changes in recombination pattern, phylogenetic relationships of some CRFs may also change. All these results may be critical to understand the role of recombination in a complex and dynamic HIV evolution.

Molecular epidemiology of HIV-1 in Panama: origin of non-B subtypes in samples collected from 2007 to 2013

Phylogenetic studies have suggested that the HIV-1 epidemic in the Americas is mainly dominated by HIV subtype B. However, countries of South America and the Caribbean have recently reported changes in their circulating HIV-1 genetic profiles. The aim of this study was to characterize the molecular profile of the HIV-1 epidemic in Panama by the analysis of 655 polymerase gene (pol) sequences that were obtained from HIV-infected Panamanians diagnosed between 1987 and 2013. Blood samples were collected from recently infected, antiretroviral drug-naïve and treatment-experienced subjects since mid-2007 to 2013. Viral RNA from plasma was extracted and sequences of HIV protease and reverse transcriptase genes were obtained. Bootscanning and phylogenetic methods were used for HIV subtyping and to trace the putative origin of non-B subtype strains. Our results showed that HIV-1 infections in Panama are dominated by subtype B (98.9%). The remaining 1.1% is represented by a diverse collection of recombinant variants including: three URFs_BC, one CRF20_BG, and one CRF28/29_BF, in addition to one subtype F1 and one subtype C, none of which were previously reported in Panama. The non-B subtype variants detected in Panama were probably introduced from Brazil (subtype F1 and CRF28/29_BF), Cuba (CRF20_BG), Dominican Republic (URFs_BC) and India (subtype C). Panama is the geographical vertex that connects the North with South America and the Caribbean through trade and cultural relations, which may explain the observed introductions of non-B subtype HIV-1 variants from both the Caribbean and South America into this Central American country.

Phylodynamics of the HIV-1 epidemic in Cuba

Previous studies have shown that the HIV-1 epidemic in Cuba displayed a complex molecular epidemiologic profile with circulation of several subtypes and circulating recombinant forms (CRF); but the evolutionary and population history of those viral variants remains unknown. HIV-1 pol sequences of the most prevalent Cuban lineages (subtypes B, C and G, CRF18_cpx, CRF19_cpx, and CRFs20/23/24_BG) isolated between 1999 and 2011 were analyzed. Maximum-likelihood analyses revealed multiple introductions of subtype B (n≥66), subtype C (n≥10), subtype G (n≥8) and CRF18_cpx (n≥2) viruses in Cuba. The bulk of HIV-1 infections in this country, however, was caused by dissemination of a few founder strains probably introduced from North America/Europe (clades B_CU-I and B_CU-II), east Africa (clade C_CU-I) and central Africa (clades G_CU, CRF18_CU and CRF19_CU), or locally generated (clades CRFs20/23/24_BG). Bayesian-coalescent analyses show that the major HIV-1 founder strains were introduced into Cuba during 1985–1995; whereas the CRFs_BG strains emerged in the second half of the 1990s. Most HIV-1 Cuban clades appear to have experienced an initial period of fast exponential spread during the 1990s and early 2000s, followed by a more recent decline in growth rate. The median initial growth rate of HIV-1 Cuban clades ranged from 0.4 year⁻¹ to 1.6 year⁻¹. Thus, the HIV-1 epidemic in Cuba has been a result of the successful introduction of a few viral strains that began to circulate at a rather late time of the AIDS pandemic, but then were rapidly disseminated through local transmission networks.

Inference of global HIV-1 sequence patterns and preliminary feature analysis

The epidemiology of HIV-1 varies in different areas of the world, and it is possible that this complexity may leave unique footprints in the viral genome. Thus, we attempted to find significant patterns in global HIV-1 genome sequences. By applying the rule inference algorithm RIPPER (Repeated Incremental Pruning to Produce Error Reduction) to multiple sequence alignments of Env sequences from four classes of compiled datasets, we generated four sets of signature patterns. We found that these patterns were able to distinguish southeastern Asian from nonsoutheastern Asian sequences with 97.5% accuracy, Chinese from non-Chinese sequences with 98.3% accuracy, African from non-African sequences with 88.4% accuracy, and southern African from non-southern African sequences with 91.2% accuracy. These patterns showed different associations with subtypes and with amino acid positions. In addition, some signature patterns were characteristic of the geographic area from which the sample was taken. Amino acid features corresponding to the phylogenetic clustering of HIV-1 sequences were consistent with some of the deduced patterns. Using a combination of patterns inferred from subtypes B, C, and all subtypes chimeric with CRF01_AE worldwide, we found that signature patterns of subtype C were extremely common in some sampled countries (for example, Zambia in southern Africa), which may hint at the origin of this HIV-1 subtype and the need to pay special attention to this area of Africa. Signature patterns of subtype B sequences were associated with different countries. Even more, there are distinct patterns at single position 21 with glycine, leucine and isoleucine corresponding to subtype C, B and all possible recombination forms chimeric with CRF01_AE, which also indicate distinct geographic features. Our method widens the scope of inference of signature from geographic, genetic, and genomic viewpoints. These findings may provide a valuable reference for epidemiological research or vaccine design.

Transmitted drug-resistance in human immunodeficiency virus-infected adult population in El Salvador, Central America

El Salvador harbours one of the largest Central American human immunodeficiency virus (HIV) epidemics, but few studies have analysed it in depth. Here, we describe the presence of transmitted drug resistance (TDR) and HIV variants in the HIV-infected adult population in El Salvador. Dried blood spots from 119 HIV-infected antiretroviral-naive adults attended in El Salvador were collected in 2011. The TDR was assessed according to the list recommended by the WHO. HIV-1 variants were described using phylogeny. Pol sequences could be amplified in 88 patients (50.6% men), with a mean age of 35 years. Almost all (96.7%) were infected with HIV through sexual practice and 58.7% were recently diagnosed. The mean CD4(+) count was 474 cells/mm(3) and 43.1% and 15.5% of patients showed moderate (<500 CD4 cells) or severe (<200) immune suppression, respectively. HIV-1 viral load was >100 000 copies/mL in 24.7% of patients and <2000 copies/mL in 9.1%. Five samples (5.7%) harboured any TDR mutation: 2.3% for nucleoside reverse transcriptase inhibitor (NRTI) and non-nucleoside reverse transcriptase inhibitor (NNRTI), and 1.4% for protease inhibitor (PI). All showed only one TDR single-class resistance mutation: M184I (two cases) for NRTI, K101E and K103N for NNRTI and L23I for PI. All viruses excepting one (URF_BG) belonged to subtype B. No phylogenetic TDR networks were found. In conclusion, we report a TDR prevalence of 5.7% in El Salvador, lower than in other Central American studies. Periodical studies are essential to monitor and prevent TDR emergence in low-income and middle-income regions. Also, more efforts are needed to promote early diagnosis and prevention of infection in El Salvador.

Construction and phenotypic characterization of HIV type 1 functional envelope clones of subtypes G and F

Subtype G has been estimated to represent the fourth most prevalent clade in the HIV-1 pandemic and subtype F is widely circulating in parts of South America (frequently within BF recombinant forms) and in Romania. However, functional envelope clones of these subtypes are lacking, which are needed for studies on antibody-mediated neutralization, coreceptor usage, and efficiency of viral entry inhibitor drugs. Here we report the construction, neutralization properties, and coreceptor usage of HIV-1 functional envelope clones of subtypes G (n = 15) and F (n = 7). These clones were obtained through RT-PCR amplification of HIV-1 gp160 from plasma RNA, and were used for pseudovirus production. All 15 subtype G-enveloped pseudoviruses were resistant to neutralization by gp120-targeted broadly neutralizing monoclonal antibodies (MAbs) b12 and 2G12, while a majority were neutralized by gp41-targeted MAbs 2F5 and 4E10. With regard to the subtype F envelopes, all seven pseudoviruses were resistant to 2F5 and b12, six were resistant to G12, and six were neutralized by 4E10. Coreceptor usage testing revealed that 21 of 22 envelopes were CCR5-tropic, including all 15 subtype G envelopes, seven of which were from patients with CD4(+) T cell counts <200/ml. These results confirm the broadly neutralizing activity of 4E10 on envelope clones across all tested group M clades, including subtypes G and F, reveal the resistance of most subtype F-enveloped pseudoviruses to broadly neutralizing MAbs b12, 2G12, and 2F5, and suggest that, similarly to subtype C, CXCR4 tropism is uncommon in subtype G, even at advanced stages of infection.

Molecular surveillance of HIV-1 in Madrid, Spain: a phylogeographic analysis

The molecular epidemiology of HIV-1 is constantly changing, mainly as a result of human migratory flows and the high adaptive ability of the virus. In recent years, Spain has become one of Europe's main destinations for immigrants and one of the western European countries with the highest rates of HIV-positive patients. Using a phylogeographic approach, we have analyzed the relationship between HIV-1 variants detected in immigrant and native populations of the urban area of Madrid. Our project was based on two coincidental facts. First, resistance tests were extended to naïve and newly diagnosed patients, and second, the Spanish government legislated the provision of legal status to many immigrants. This allowed us to obtain a large data set (n = 2,792) from 11 Madrid hospitals of viral pol sequences from the two populations, and with this unique material, we explored the impact of immigration in the epidemiological trends of HIV-1 variants circulating in the largest Spanish city. The prevalence of infections by non-B HIV-1 variants in the studied cohort was 9%, rising to 25% among native Spanish patients. Multiple transmission events involving different lineages and subsubtypes were observed in all the subtypes and recombinant forms studied. Our results also revealed strong social clustering among the most recent immigrant groups, such as Russians and Romanians, but not in those groups who have lived in Madrid for many years. Additionally, we document for the first time the presence of CRF47_BF and CRF38_BF in Europe, and a new BG recombinant form found in Spaniards and Africans is tentatively proposed. These results suggest that the HIV-1 epidemic will evolve toward a more complex epidemiological landscape.

Molecular epidemiology of HIV type 1 in Mexico: emergence of BG and BF intersubtype recombinants

The molecular epidemiology of subtypes and intersubtype recombinants (IRs) of human immunodeficiency virus type 1 (HIV-1) in Mexico has not been characterized fully. Understanding its regional distribution, prevalence, adaptability, viral fitness, pathogenicity, and immunogenicity is decisive for any design of an effective HIV vaccine. The aim of this study was to describe the presence of IRs types BG and BF in a Mexican population. Protease and reverse transcriptase regions of the pol gene were sequenced using an automated sequencing system. A phylogenic tree was constructed and genetic distances were calculated using MEGA 3.1. Recombination analysis was done by bootscan using SimPlot software. Two hundred and twenty-three HIV-1-positive individuals were enrolled in the study. At baseline, the mean plasma viral load was 285,500 HIV-1 RNA copies/ml and the mean CD4 cell count was 213 cells/ml. Subtype B was found in 220 (98.6%) samples, whereas IRs were found in three patients (1.4%): two (0.9%) with BG and one (0.45%) with BF. IRs were observed in 2/124 (1.6%) samples from treated patients and in 1/99 (1.0%) from naive patients. The presence of these HIV forms at low frequency points to the need for research on the diversity, geographic distribution, and evolution of other subtypes including circulating recombinant forms and IRs to understand the molecular epidemiology and tendencies of the HIV infection in Mexico.

Heterogeneity and penetration of HIV-1 non-subtype B viruses in an Italian province: public health implications

This study assessed changes in prevalence and distribution of HIV-1 non-subtype B viruses in Italian and immigrant patients over two decades in a province in Italy. All HIV-positive patients who underwent genotypic resistance testing were selected. Prevalence of non-subtype B viruses in 3-year periods was calculated. All sequences of non-subtype B and those provided by REGA as unassigned were analysed for phylogenetic relationships. In total, 250/1563 (16%) individuals were infected with a non-subtype B virus. Prevalence increased over time, reaching a peak (31.5%) in 2004-2006. In Italian patients, the most frequent subtypes were B (92.5%) and F1 (4%). F1 subtype was also prevalent in patients from South America (13.6%); in patients of African origin, CRF02_AG (54.9%) and G (12.3%) were the most frequent. HIV-1 non-subtype B infections in Italians were mostly found in patients who acquired HIV sexually. A phylogenetic relationship between F subtypes in Italian and representative HIV-1 sequences from Brazil was found. C subtypes in Italians were phylogenetically related to subtypes circulating in Brazil. Inter-subtype recombinants were also found in the latest years. The HIV-1 epidemic in Brescia province evolved to the point where about 1/3 patients recently diagnosed harboured non-B HIV subtypes. The distribution of HIV-1 non-B subtypes in Italian patients resembled that in South American patients and phylogenetic relatedness between some Italian and South American HIV-1 strains was found. The possible epidemiological link between these two populations would have been missed by looking only at risk factors for HIV acquisition declared by patients. The evidence of inter-subtype recombinants points to significant genetic assortment. Overall our results support phylogenetic analysis as a tool for epidemiological investigation in order to guide targeted prevention strategies.

Temporal and spatial dynamics of human immunodeficiency virus type 1 circulating recombinant forms 08_BC and 07_BC in Asia

Human immunodeficiency virus type 1 (HIV-1) CRF08_BC and CRF07_BC are two major recombinants descended from subtypes B' and C. Despite their massive epidemic impact in China, their migration patterns and divergence times remain unknown. Phylogenetic and population genetic analyses were performed on 228 HIV-1 sequences representing CRF08_BC, CRF07_BC, and subtype C strains from different locations across China, India, and Myanmar. Genome-specific rates of evolution and divergence times were estimated using a Bayesian Markov chain Monte Carlo framework under various evolutionary models. CRF08_BC originated in 1990.3 (95% credible region [CR], 1988.6 to 1991.9) in Yunnan province before spreading to Guangxi (south) and Liaoning (northeast) around 1995. Inside Guangxi region, the eastward expansion of CRF08_BC continued from Baise city (west) to Binyang (central) between 1997 and 1998 and later spread into Pingxiang around 1999 in the south, mainly through injecting drug users. Additionally, CRF07_BC diverged from its common ancestor in 1993.3 (95% CR, 1991.2 to 1995.2) before crossing the border into southern Taiwan in late 1990s. Phylogenetic analysis indicates that both CRF08_BC and CRF07_BC can trace their origins to Yunnan. The parental Indian subtype C lineage likely entered China around 1981.2 (95% CR, 1976.7 to 1985.9). Using a multiple unlinked locus model, we also showed that the dates of divergence calculated in this study may not be significantly affected by intrasubtype recombination among different lineages. This is the first phylodynamic study depicting the spatiotemporal dynamics of HIV/AIDS in East Asia.

Long-range recombination gradient between HIV-1 subtypes B and C variants caused by sequence differences in the dimerization initiation signal region

HIV-1 intersubtype recombinants have an increasingly important role in shaping the AIDS pandemic. We sought to understand the molecular mechanisms that generate intersubtype HIV-1 recombinants. We analyzed recombinants of HIV-1 subtypes B and C, and identified their crossover junctions in the viral genome from the 5' long terminal repeat (LTR) to the end of pol. We identified 56 recombination events in 56 proviruses; the distribution of these events indicated an apparent recombination gradient: there were significantly more crossover junctions in the 3' half than in the 5' half of the region analyzed. HIV-1 subtypes B and C have different dimerization initiation signal (DIS). We hypothesized that the inability of subtype B and C RNAs to form perfect base-pairing of the DIS affects the dimeric RNA structure and causes a decrease in recombination events at the 5' end of the viral genome. To test this hypothesis, we examined recombinants generated from a subtype C virus and a modified subtype B virus containing a subtype C DIS. In the 56 proviruses analyzed, we identified 96 recombination events, which are significantly more frequent than in the B/C recombinants. Furthermore, these crossover junctions were distributed evenly throughout the region analyzed, indicating that the recombination gradient was corrected by matching the DIS. Therefore, base-pairing at the DIS has an important function during HIV-1 reverse transcription, most likely in maintaining nucleic-acid structure in the complex. These findings reveal elements important to retroviral recombination and provide insights into the generation of HIV-1 intersubtype recombinants that are important to the AIDS epidemic.