HIV Type 1 molecular epidemiology in cuba: high genetic diversity, frequent mosaicism, and recent expansion of BG intersubtype recombinant forms

Cuban history of CRF19 recombinant subtype of HIV-1

CRF19 is a recombinant form of HIV-1 subtypes D, A1 and G, which was first sampled in Cuba in 1999, but was already present there in 1980s. CRF19 was reported almost uniquely in Cuba, where it accounts for ∼25% of new HIV-positive patients and causes rapid progression to AIDS (∼3 years). We analyzed a large data set comprising ∼350 pol and env sequences sampled in Cuba over the last 15 years and ∼350 from Los Alamos database. This data set contained both CRF19 (∼315), and A1, D and G sequences. We performed and combined analyses for the three A1, G and D regions, using fast maximum likelihood approaches, including: (1) phylogeny reconstruction, (2) spatio-temporal analysis of the virus spread, and ancestral character reconstruction for (3) transmission mode and (4) drug resistance mutations (DRMs). We verified these results with a Bayesian approach. This allowed us to acquire new insights on the CRF19 origin and transmission patterns. We showed that CRF19 recombined between 1966 and 1977, most likely in Cuban community stationed in Congo region. We further investigated CRF19 spread on the Cuban province level, and discovered that the epidemic started in 1970s, most probably in Villa Clara, that it was at first carried by heterosexual transmissions, and then quickly spread in the 1980s within the "men having sex with men" (MSM) community, with multiple transmissions back to heterosexuals. The analysis of the transmission patterns of common DRMs found very few resistance transmission clusters. Our results show a very early introduction of CRF19 in Cuba, which could explain its local epidemiological success. Ignited by a major founder event, the epidemic then followed a similar pattern as other subtypes and CRFs in Cuba. The reason for the short time to AIDS remains to be understood and requires specific surveillance, in Cuba and elsewhere.

Increasing prevalence and local transmission of non-B HIV-1 subtypes in the French Antilles and French Guiana between 1995 and 2018

The Caribbean and South American French Overseas Territories (CSAFOT) are the regions most heavily affected by the Human Immunodeficiency Virus type 1 (HIV-1) epidemic in France. Although dominated by HIV-1 subtype B, the detection of non-B subtypes and the great proportion of HIV-positive persons born abroad demonstrated the potential for local spread of non-B subtype strains in CSAFOT. To reconstruct the epidemiologic dynamics of major non-B subtype clusters spreading in CSAFOT, we conducted phylogenetic and evolutionary analyses of 2,523 HIV-1 pol sequences collected from patients living in Martinique, Guadeloupe, and French Guiana from 1995 to 2018. A large variety of HIV-1 non-B subtype strains (eight subtypes, twelve CRFs, and multiple URFs) have been introduced in CSAFOT and their prevalence significantly increases over time in Martinique and Guadeloupe. We identified twelve major transmission networks of non-B subtypes (CRF02_AG and subtypes A3, C, D, and F1) that probably arose in Guadeloupe, Martinique, French Guiana, and mainland France between the late 1970s and the middle 2000s. Phylogeographic analyses support frequent non-B subtype viral transmissions within CSAFOT as well as transatlantic transmission between CSAFOT and mainland France. Domestic transmission networks of non-B subtype variants in CSAFOT comprise both men having sex with men and heterosexual individuals from different age groups. Different HIV-1 non-B subtype variants were sequentially introduced in CSAFOT between the late 1970s and the middle 2000s and are currently spreading through domestic, regional, and/or transatlantic networks of individuals from different age and risk groups.

Emergence as an outbreak of the HIV-1 CRF19_cpx variant in treatment-naïve patients in southern Spain

Background

CRF19_cpx is a complex circulating recombination form (CRF) of HIV-1. We describe the characteristics of an outbreak of the CRF19_cpx variant among treatment-naïve patients in southern Spain.

Methods

The study was undertaken at the Virgen de la Victoria Hospital, a reference centre for the analysis of HIV-1 genotype in Malaga (Spain). Subtyping was performed through REGA v3.0 and the relationship of our CRF19_cpx sequences, among themselves and regarding other reference sequences from the same variant, was defined by phylogenetic analysis. We used PhyML program to perform a reconstruction of the phylogeny by Maximum Likelihood method as well as further confirmation of the transmission clusters by Bayesian inference. Additionally, we collected demographic, clinical and immunovirological data.

Results

Between 2011 and 2016, we detected 57 treatment-naïve patients with the CRF19_cpx variant. Of these, 55 conformed a very well-defined transmission cluster, phylogenetically close to CRF19_cpx sequences from the United Kingdom. The origin of this subtype in Malaga was dated between 2007 and 2010. Over 50% of the patients presented the non-nucleoside reverse transcriptase inhibitor G190A resistance mutation. This variant was mostly represented by young adult Spanish men who had sex with men. Almost half of them were recent seroconverters, though a similar percentage was diagnosed at a late state of HIV infection. Five cases of AIDS and one non-AIDS defined death occurred during follow-up. The majority of patients treated with first-line combination antiretroviral therapy (ART) responded.

Conclusions

We report the largest HIV-1 CRF19_cpx cohort of treatment-naïve patients outside Cuba, almost all emerging as an outbreak in the South of Spain. Half the cases had the G190A resistance mutation. Unlike previous studies, the variant from Malaga seems less pathogenic, with few AIDS events and an excellent response to ART.

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.

HIV-1 subtype B: Traces of a pandemic

Human migration is a major process that shaped the origin and dissemination of HIV. Within HIV-1, subtype B (HIV-1B) is the most disseminated variant and it is assumed to be the causative agent in approximately 11% of all cases of HIV worldwide. Phylogenetic studies have revealed that HIV-1B emerged in Kinshasa (Africa) and was introduced into the Caribbean region via Haiti in or around 1966 by human migration. After localized dispersion, the virus was brought to the United States of America via homosexual/bisexual contact around 1969. Inside USA, the incidence of HIV-1B infection increased exponentially and it became established in the population, affecting not only homosexual individuals but also heterosexual individuals and injecting drug users. Soon after, the virus was disseminated and became established in other regions, including Europe, Asia, Latin America, and Australia. Recent studies suggest that, in addition to this pandemic clade, several lineages have emerged from Haiti and reached other Caribbean and Latin American countries via short-distance dissemination. Different subtype B genetic variants have also been detected in these epidemics. Four genetic variants have been described to date: subtype B', which mainly circulates in Thailand and other Asian countries; a specific variant mainly found in Trinidad and Tobago; the GPGS variant, which is primarily detected in Korea; and the GWGR variant, which is mainly detected in Brazil. This paper reviews the evolution of HIV-1B and its impact on the human population.

Co-infection with two strains of Brome mosaic bromovirus reveals common RNA recombination sites in different hosts

We have previously reported intra-segmental crossovers in Brome mosaic virus (BMV) RNAs. In this work, we studied the homologous recombination of BMV RNA in three different hosts: barley (Hordeum vulgare), Chenopodium quinoa, and Nicotiana benthamiana that were co-infected with two strains of BMV: Russian (R) and Fescue (F). Our work aimed at (1) establishing the frequency of recombination, (2) mapping the recombination hot spots, and (3) addressing host effects. The F and R nucleotide sequences differ from each other at many translationally silent nucleotide substitutions. We exploited this natural variability to track the crossover sites. Sequencing of a large number of cDNA clones revealed multiple homologous crossovers in each BMV RNA segment, in both the whole plants and protoplasts. Some recombination hot spots mapped at similar locations in different hosts, suggesting a role for viral factors, but other sites depended on the host. Our results demonstrate the chimeric ('mosaic') nature of the BMV RNA genome.

Phylogeny and Phylogeography of a Recent HIV-1 Subtype F Outbreak among Men Who Have Sex with Men in Spain Deriving from a Cluster with a Wide Geographic Circulation in Western Europe

We recently reported the rapid expansion of an HIV-1 subtype F cluster among men who have sex with men (MSM) in the region of Galicia, Northwest Spain. Here we update this outbreak, analyze near full-length genomes, determine phylogenetic relationships, and estimate its origin. For this study, we used sequences of HIV-1 protease-reverse transcriptase and env V3 region, and for 17 samples, near full-length genome sequences were obtained. Phylogenetic analyses were performed via maximum likelihood. Locations and times of most recent common ancestors were estimated using Bayesian inference. Among samples analyzed by us, 100 HIV-1 F1 subsubtype infections of monophyletic origin were diagnosed in Spain, including 88 in Galicia and 12 in four other regions. Most viruses (n = 90) grouped in a subcluster (Galician subcluster), while 7 from Valladolid (Central Spain) grouped in another subcluster. At least 94 individuals were sexually-infected males and at least 71 were MSM. Seventeen near full-length genomes were uniformly of F1 subsubtype. Through similarity searches and phylogenetic analyses, we identified 18 viruses from four other Western European countries [Switzerland (n = 8), Belgium (n = 5), France (n = 3), and United Kingdom (n = 2)] and one from Brazil, from samples collected in 2005–2011, which branched within the subtype F cluster, outside of both Spanish subclusters, most of them corresponding to recently infected individuals. The most probable geographic origin and age of the Galician subcluster was Ferrol, Northwest Galicia, around 2007, while the Western European cluster probably emerged in Switzerland around 2002. In conclusion, a recently expanded HIV-1 subtype F cluster, the largest non-subtype B cluster reported in Western Europe, continues to spread among MSM in Spain; this cluster is part of a larger cluster with a wide geographic circulation in diverse Western European countries.

Expansion of the CRF19_cpx Variant in Spain

BACKGROUND

HIV-1 CRF19_cpx, is a recombinant variant found almost exclusively in Cuba and recently associated to a faster AIDS onset. Infection with this variant leads to higher viral loads and levels of RANTES and CXCR4 co-receptor use.

OBJECTIVES

The goal of this study was to assess the presence of CRF19_cpx in the Spanish province of Valencia, given its high pathogenicity.

STUDY DESIGN

1294 HIV-1 protease-reverse transcriptase (PR/RT) sequences were obtained in Valencia (Spain), between 2005 and 2014. After subtyping, the detected CRF19_cpx sequences were aligned with 201 CRF19_cpx and 66 subtype D sequences retrieved from LANL, and subjected to maximum-likelihood phylogenetic analyses and Bayesian coalescent reconstructions. The presence of resistance mutations in the PR/RT region of these sequences was also analyzed.

RESULTS

Among the 9 CRF19_cpx sequences from different patients found (prevalence <0.1%), 7 grouped in two well-supported clades (groups A, n=4, and B, n=3), suggesting the existence of at least two independent introductions which subsequently started to expand in the studied Spanish region. Unprotected sex between men was the only known transmission route. Coalescent analyses suggested that the introductions in Valencia occurred between 2008 and 2010. Resistance mutations in the RT region were found in all sequences from group A (V139D) and in two sequences from group B (E138A).

CONCLUSIONS

This study reports for the first time the recent expansion of CRF19_cpx outside Cuba. Our results suggest that CRF19_cpx might become an emerging HIV variant in Spain, affecting Spanish native MSM and not only Cuban migrants.

Performance of an in-house human immunodeficiency virus type 1 genotyping system for assessment of drug resistance in Cuba

As commercial human immunodeficiency virus type 1 drug resistance assays are expensive, they are not commonly used in resource-limited settings. Hence, a more affordable in-house procedure was set up taking into account the specific epidemiological and economic circumstances of Cuba. The performance characteristics of the in-house assay were evaluated using clinical samples with various subtypes and resistance patterns. The lower limit of amplification was determined on dilutions series of 20 clinical isolates and ranged from 84 to 529 RNA copies/mL. For the assessment of trueness, 14 clinical samples were analyzed and the ViroSeq HIV-1 Genotyping System v2.0 was used as the reference standard. The mean nucleotide sequence identity between the two assays was 98.7% ± 1.0. Additionally, 99.0% of the amino acids at drug resistance positions were identical. The sensitivity and specificity in detecting drug resistance mutations was respectively 94.1% and 99.5%. Only few discordances in drug resistance interpretation patterns were observed. The repeatability and reproducibility were evaluated using 10 clinical samples with 3 replicates per sample. The in-house test was very precise as nucleotide sequence identity among paired nucleotide sequences ranged from 98.7% to 99.9%. The acceptance criteria were met by the in-house test for all performance characteristics, demonstrating a high degree of accuracy. Subsequently, the applicability in routine clinical practice was evaluated on 380 plasma samples. The amplification success rate was 91% and good quality consensus sequences encoding the entire protease and the first 335 codons in reverse transcriptase could be obtained for 99% of the successful amplicons. The reagent cost per sample using the in-house procedure was around € 80 per genotyping attempt. Overall, the in-house assay provided good results, was feasible with equipment and reagents available in Cuba and was half as expensive as commercial assays.

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.

Short communication: Molecular epidemiology, phylogeny, and phylodynamics of CRF63_02A1, a recently originated HIV-1 circulating recombinant form spreading in Siberia

The HIV-1 epidemic in Russia is dominated by the former Soviet Union subtype A (A(FSU)) variant, but other genetic forms are circulating in the country. One is the recently described CRF63_02A1, derived from recombination between a CRF02_AG variant circulating in Central Asia and A(FSU), which has spread in the Novosibirsk region, Siberia. Here we phylogenetically analyze pol and env segments from 24 HIV-1 samples from the Novosibirsk region collected in 2013, with characterization of three new near full-length genome CRF63_02A1 sequences, and estimate the time of the most recent common ancestor (tMRCA) and the demographic growth of CRF63_02A1 using a Bayesian method. The analyses revealed that CRF63_02A1 is highly predominant in the Novosibirsk region (81.2% in pol sequences) and is transmitted both among injecting drug users and by heterosexual contact. Similarity searches with database sequences combined with phylogenetic analyses show that CRF63_02A1 is circulating in East Kazakhstan and the Eastern area of Russia bordering China. The analyses of near full-length genome sequences show that its mosaic structure is more complex than reported, with 18 breakpoints. The tMRCA of CRF63_02A1 was estimated around 2006, with exponential growth in 2008-2009 and subsequent stabilization. These results provide new insights into the molecular epidemiology, phylogeny, and phylodynamics of CRF63_02A1.

Spatiotemporal dynamics of dissemination of non-pandemic HIV-1 subtype B clades in the Caribbean region

The Human immunodeficiency virus type-1 (HIV-1) epidemic in the Caribbean region is mostly driven by subtype B; but information about the pattern of viral spread in this geographic region is scarce and different studies point to quite divergent models of viral dissemination. In this study, we reconstructed the spatiotemporal and population dynamics of the HIV-1 subtype B epidemic in the Caribbean. A total of 1,806 HIV-1 subtype B pol sequences collected from 17 different Caribbean islands between 1996 and 2011 were analyzed together with sequences from the United States (n = 525) and France (n = 340) included as control. Maximum Likelihood phylogenetic analyses revealed that HIV-1 subtype B infections in the Caribbean are driven by dissemination of the pandemic clade (BPANDEMIC) responsible for most subtype B infections across the world, and older non-pandemic lineages (BCAR) characteristics of the Caribbean region. The non-pandemic BCAR strains account for >40% of HIV-1 infections in most Caribbean islands; with exception of Cuba and Puerto Rico. Bayesian phylogeographic analyses indicate that BCAR strains probably arose in the island of Hispaniola (Haiti/Dominican Republic) around the middle 1960s and were later disseminated to Trinidad and Tobago and to Jamaica between the late 1960s and the early 1970s. In the following years, the BCAR strains were also disseminated from Hispaniola and Trinidad and Tobago to other Lesser Antilles islands at multiple times. The BCAR clades circulating in Hispaniola, Jamaica and Trinidad and Tobago appear to have experienced an initial phase of exponential growth, with mean estimated growth rates of 0.35-0.45 year(-1), followed by a more recent stabilization since the middle 1990s. These results demonstrate that non-pandemic subtype B lineages have been widely disseminated through the Caribbean since the late 1960s and account for an important fraction of current HIV-1 infections in the region.

Prevalence of HIV-1 subtypes among men who have sex with men in China: a systematic review

The distribution of prevalent HIV-1 strains are still complex in China. Men who have sex with men (MSM) play an important bridging role in spreading HIV. The aim of our study was to quantitatively evaluate the prevalence of HIV-1 subtypes among the MSM population in China from published studies. Relevant studies were searched by selection criteria from CNKI, CBM, Pubmed, etc. We computed the estimates of the pooled proportion of HIV-1 subtypes. Heterogeneity between studies was investigated and measured using Cochran's Q statistic and the I (2) statistic. All analyses were conducted by the R statistical package version 2.13.1. A meta-analysis was performed, which included 19 articles. For comprehensive analysis of env, gag and pol genes, the pooled estimates for the prevalence of subtype B was 28.25% (95% CI: 18.10-39.66%), CRF01_AE was 53.46% (95% CI: 46.11-60.74%), CRF07_BC was 18.66% (95% CI: 13.06-25.01%) and CRF08_BC was 5.85% (95% CI: 2.73-10.07%), respectively. In subgroup analysis, the proportion of subtype B decreased, while the proportion of CRF01_AE and CRF07_BC showed an increasing tendency. Beijing, Guangdong and Henan provinces had high proportions of subtype CRF01_AE while Guangdong and Hebei provinces had the highest proportions of subtype B and CRF07_BC, respectively. A high genetic variability of HIV-1 presents a serious challenge for HIV prevention and treatment strategies among MSM in China.

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.

Current Trends of HIV Recombination Worldwide

One of the major characteristics of HIV-1 is its high genetic variability and extensive heterogeneity. This characteristic is due to its molecular traits, which in turn allows it to vary, recombine, and diversify at a high frequency. As such, it generates complex molecular forms, termed recombinants, which evade the human immune system and so survive. There is no sequence constraint to the recombination pattern as it appears to occur at inter-group (between groups M and O), as well as interand intra-subtype within group M. Rapid emergence and active global transmission of HIV-1 recombinants, known as circulating recombinant forms (CRFs) and unique recombinant forms (URFs), requires urgent attention. To date, 55 CRFs have been reported around the world. The first CRF01_AE originated from Central Africa but spread widely in Asia. The most recent CRF; CRF55_01B is a recombinant form of CRF01_AE and subtype B, although its origin is yet to be publicly disclosed. HIV-1 recombination is an ongoing event and plays an indispensable role in HIV epidemics in different regions. Africa, Asia and South America are identified as recombination hot-spots. They are affected by continual emergence and cocirculation of newly emerging CRFs and URFs, which are now responsible for almost 20% of HIV-1 infections worldwide. Better understanding of recombinants is necessary to determine their biological and molecular attributes.

Antiretroviral drug resistance in HIV-1 therapy-naive patients in Cuba

In Cuba, antiretroviral therapy rollout started in 2001 and antiretroviral therapy coverage has reached almost 40% since then. The objectives of this study were therefore to analyze subtype distribution, and level and patterns of drug resistance in therapy-naive HIV-1 patients. Four hundred and one plasma samples were collected from HIV-1 therapy-naive patients in 2003 and in 2007-2011. HIV-1 drug resistance genotyping was performed in the pol gene and drug resistance was interpreted according to the WHO surveillance drug-resistance mutations list, version 2009. Potential impact on first-line therapy response was estimated using genotypic drug resistance interpretation systems HIVdb version 6.2.0 and Rega version 8.0.2. Phylogenetic analysis was performed using Neighbor-Joining. The majority of patients were male (84.5%), men who have sex with men (78.1%) and from Havana City (73.6%). Subtype B was the most prevalent subtype (39.3%), followed by CRF20-23-24_BG (19.5%), CRF19_cpx (18.0%) and CRF18_cpx (10.3%). Overall, 29 patients (7.2%) had evidence of drug resistance, with 4.0% (CI 1.6%-4.8%) in 2003 versus 12.5% (CI 7.2%-14.5%) in 2007-2011. A significant increase in drug resistance was observed in recently HIV-1 diagnosed patients, i.e. 14.8% (CI 8.0%-17.0%) in 2007-2011 versus 3.8% (CI 0.9%-4.7%) in 2003 (OR 3.9, CI 1.5-17.0, p=0.02). The majority of drug resistance was restricted to a single drug class (75.8%), with 55.2% patients displaying nucleoside reverse transcriptase inhibitor (NRTI), 10.3% non-NRTI (NNRTI) and 10.3% protease inhibitor (PI) resistance mutations. Respectively, 20.7% and 3.4% patients carried viruses containing drug resistance mutations against NRTI+NNRTI and NRTI+NNRTI+PI. The first cases of resistance towards other drug classes than NRTI were only detected from 2008 onwards. The most frequent resistance mutations were T215Y/rev (44.8%), M41L (31.0%), M184V (17.2%) and K103N (13.8%). The median genotypic susceptibility score for the commonly prescribed first-line therapies was 2.5. This analysis emphasizes the need to perform additional surveillance studies to accurately assess the level of transmitted drug resistance in Cuba, as the extent of drug resistance might jeopardize effectiveness of first-line regimens prescribed in Cuba and might necessitate the implementation of baseline drug resistance testing.

Generation of a replication-competent chimeric simian-human immunodeficiency virus carrying env from subtype C clinical isolate through intracellular homologous recombination

A new simian-human immunodeficiency virus (SHIV), carrying env from an uncloned HIV-1 subtype C clinical isolate (97ZA012), was generated through intracellular homologous recombination, a DNA repair mechanism of the host cell. PCR fragments amplified from an existing SHIV plasmid (a 7-kb fragment from the 5' end and a 1.5-kb fragment from the 3' end) and a 4-kb fragment amplified from 97ZA012 cDNA containing env were co-transfected to human lymphoid cells. The resulting recombinant was subjected to serial passage in rhesus peripheral blood mononuclear cells (RhPBMCs). The resulting SHIV 97ZA012 was replication competent in RhPBMCs and monkey alveolar macrophages, and possessed CCR5 preference as an entry co-receptor. Experimental infection of rhesus macaques with SHIV 97ZA012 caused high titers of plasma viremia and a transient but profound depletion of CD4(+) T lymphocytes in the lung. Animal-to-animal passage was shown to be a promising measure for further adaptation of the virus in monkeys.

HIV type 1 genetic diversity in newly diagnosed Cuban patients

Knowledge of the genetic diversity of HIV-1 constitutes a fundamental premise in the epidemiological surveillance. In the present study, the HIV-1 genetic variability from 142 Cuban patients who were diagnosed with HIV-1 infection during 2009 and 2010 was determined. HIV-1 subtypes were determined by partial RT-PCR and sequencing of the HIV-1 pol gene. The phylogenetic analysis showed that 47 (33.1 %) samples were subtypes B and 95 (66.9 %) were non-B subtypes, where G, H, and C subtypes, as well as the recombinant forms CRF19_cpx, CRF18_cpx, and CRFs BG, were included. The circulation of CRF05_DF was detected for the first time in Cuba. The analyses of recombinants showed the presence of recombinant CRF18_cpx/CRF19_cpx. The study confirms the high genetic diversity of HIV-1 and the circulation of new genetic variants in the studied population, which indicates the importance of maintaining constant epidemiological surveillance in Cuba.

Increasing diversity of Human Immunodeficiency Virus type 1 subtypes circulating in Australia

Characterization of HIV subtypes can provide a more comprehensive understanding of the epidemic within a distinct region, and when combined with notification data, may also be helpful in enhancing current HIV prevention strategies. In this study, we characterized 1056 HIV-positive individuals (948 males and 108 females) living in Victoria and whose infection was detected for the first time between 2005 and 2010 inclusive. HIV-1 strains were subtyped based on pol gene sequence. Phylogenetic analysis was performed on all non-B subtype sequences identified. Of the 1056 sequences analyzed, 825 were subtype B and 231 were non-B. Overall 6 HIV-1 subtypes, 6 circulating recombinant forms (CRFs), and 12 unique recombinant forms (URFs) were identified. Regardless of gender, the majority of individuals were infected with a subtype B virus (78%). Subtype B was dominant in males (n=806, 85%). In contrast, the majority of females were infected with non-B subtypes (n=89, 82%), in particular subtype C (n=48, 45%). Phylogenetic analysis of the non-B subtypes revealed that the majority of clustering, and thereby transmission, occurred with CRF01_AE strains. Despite the relatively high numbers identified in females there was very little clustering of subtype C viruses. Subtypes C and A1 both historically associated with heterosexual transmission, and CRF01_AE often associated with IVDU, were also associated with transmission within the MSM population, demonstrating the potential for non-B subtypes to expand into the MSM population. The observation of increasing numbers of females and heterosexual males infected with non-subtype B viruses, the majority imported through migration and travel to countries where there is a high prevalence of HIV, suggests a targeted public health message may be required to prevent further increases within these two groups.

Quality of life of people with HIV/AIDS receiving antiretroviral therapy in Cuba: a cross-sectional study of the national population

OBJECTIVES

We studied the effect of antiretroviral therapy (ART) on the quality of life (QOL) of Cubans with HIV/AIDS.

METHODS

We conducted a cross-sectional study including administration of the Medical Outcomes Study-HIV Health Survey Questionnaire to a representative sample of the 1592 Cubans receiving ART in 2004. For univariate analyses, we compared mean HIV scale scores. We used logistic regression models to estimate the association between role function and year of diagnosis, between pain and sex, and between health transition and region of diagnosis, with adjustment for demographics, ART regimen, and clinical status.

RESULTS

There were 354 participants (73 women, 281 men). Scores for all functional activities showed means higher than 80 out of 100. Pain interfered more in women than in men (73.2 vs 81.9; P = .01). When HIV diagnosis occurred after 2001, the probability of experiencing difficulties performing work (odds ratio [OR] = 4.42; 95% CI = 1.83, 10.73) and pain (OR = 1.70; 95% CI = 1.01, 2.88) increased compared with earlier diagnosis. People treated with indinavir showed a greater perception of general health (58.9 vs 52.4; P = .045) and greater health improvement (78.6 vs 67.8; P = .002).

CONCLUSIONS

Although Cubans receiving ART are maintaining a high QOL, we observed significant differences by sex and time of diagnosis. QOL assessment can serve as a health outcome and may allow identification of QOL reductions potentially related to ART side effects.

The origin and diversity of the HIV-1 pandemic

This review examines the enormous progress that has been made in the past decade in understanding the origin of HIV, HIV genetic variability, and the impact of global HIV diversity on the pandemic. Multiple zoonotic transmissions of simian immunodeficiency virus (SIV) have resulted in different HIV lineages in humans. In addition, the high mutation and recombination rates during viral replication result in a great genetic variability of HIV within individuals, as well as within populations, upon which evolutionary selection pressures act. The global HIV pandemic is examined in the context of HIV evolution, and the global diversity of HIV subtypes and recombinants is discussed in detail. Finally, the impact of HIV diversity on pathogenesis, transmission, diagnosis, treatment, the immune response, and vaccine development is reviewed.

Origin and epidemiological history of HIV-1 CRF14_BG

Background

CRF14_BG isolates, originally found in Spain, are characterized by CXCR4 tropism and rapid disease progression. This study aimed to identify the origin of CRF14_BG and reconstruct its epidemiological history based on new isolates from Portugal.

Methodology/Principal Findings

C2V3C3 env gene sequences were obtained from 62 samples collected in 1993–1998 from Portuguese HIV-1 patients. Full-length genomic sequences were obtained from three patients. Viral subtypes, diversity, divergence rate and positive selection were investigated by phylogenetic analysis. The molecular structure of the genomes was determined by bootscanning. A relaxed molecular clock model was used to date the origin of CRF14_BG. Geno2pheno was used to predict viral tropism. Subtype B was the most prevalent subtype (45 sequences; 73%) followed by CRF14_BG (8; 13%), G (4; 6%), F1 (2; 3%), C (2; 3%) and CRF02_AG (1; 2%). Three CRF14_BG sequences were derived from 1993 samples. Near full-length genomic sequences were strongly related to the CRF14_BG isolates from Spain. Genetic diversity of the Portuguese isolates was significantly higher than the Spanish isolates (0.044 vs 0.014, P<0.0001). The mean date of origin of the CRF14_BG cluster was estimated to be 1992 (range, 1989 and 1996) based on the subtype G genomic region and 1989 (range, 1984–1993) based on the subtype B genomic region. Most CRF14_BG strains (78.9%) were predicted to be CXCR4. Finally, up to five amino acids were under selective pressure in subtype B V3 loop whereas only one was found in the CRF14_BG cluster.

Conclusions

CRF14_BG emerged in Portugal in the early 1990 s soon after the beginning of the HIV-1 epidemics, spread to Spain in late 1990 s as a consequence of IVDUs migration and then to the rest of Europe. CXCR4 tropism is a general characteristic of this CRF that may have been selected for by escape from neutralizing antibody response.

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.

Short communication: Molecular epidemiology of HIV type 1 in the Republic of Dagestan, Russian Federation: virtually uniform circulation of subtype A, former Soviet Union variant, with predominance of the V77I(PR) subvariant

We examine the distribution of viral genetic forms and the presence of antiretroviral drug resistance mutations in HIV-1 infections in the Republic of Dagestan, in the North Caucasus area of Russia, where a recent large increase in HIV-1 infections has been documented. Samples were collected from 41 HIV-1-infected individuals from Dagestan, most of them from the cities of Derbent (n = 21) and Mahachkala (n = 18). Thirty six were injecting drug users and five were infected by heterosexual contact. None was on antiretroviral drug treatment. HIV-1 protease and a segment of reverse transcriptase were amplified by RT-PCR from plasma RNA and sequenced, and phylogenetic trees were constructed via maximum likelihood. Forty (97.6%) of 41 samples were of subtype A, former Soviet Union variant (A(FSU)), of which 27 (67.5%) clustered with the subvariant containing the V77I substitution in protease (V77I(PR)). Within this cluster, 13 viruses formed a local subcluster, 10 of which were from Derbent. Four viruses clustered with the A(SP2) subcluster, recently identified in St. Petersburg, two with a virus from Georgia and one with a virus from Azerbaijan. No mutations associated with antiretroviral drug resistance were detected. The results, therefore, show the relationship of the HIV-1 epidemic in Dagestan with that of other areas of Russia and of neighboring countries, and reveal the spread of the A(FSU) V77I(PR) variant in the North Caucasus area.

Evaluation of HIV type 1 strains in men having sex with men and in female sex workers in Mombasa, Kenya

We compared HIV-1 strains in incident and prevalent infections in a cohort of men having sex with men (MSM) and female sex workers (FSW) near Mombasa, Kenya and conducted a cross-sectional study of viral isolates from a sample of HIV-1-infected MSM and FSW in Kilifi, Coast Province, Kenya. RNA extracted from plasma of 13 MSM, 9 FSW, and one heterosexual male was amplified by nested RT-PCR and the products were directly sequenced. HIV-1 strains from 21 individuals were characterized with one or more complete genome sequences, and two were sequenced in the Nef gene. The envelope quasispecies was also studied in one individual. Among MSM, eight strains were subtype A and five were recombinant. There were two epidemiologically linked pairs of sequences; one pair was subtype A and the other pair was a complex AA2CD recombinant of identical structure. Another MSM was dually infected with DG recombinant strains of related, but nonidentical, structure. MSM also harbored AC and AD recombinant strains. The FSW harbored seven subtype A strains, an AD recombinant, and an AA2D strain related to CRF16_A2D. The one heterosexual male studied had a subtype A infection. This MSM epidemic in Kenya appears to be of local origin, harboring many strains typical of the broader Kenyan epidemic. Characteristics of a close social network were identified, with extended chains of transmission, novel recombinant strains possibly generated within the network, and a relatively high proportion of recombinant and dual infections.

Resistance and viral subtypes: how important are the differences and why do they occur?

PURPOSE OF REVIEW

The global HIV-1 pandemic has evolved to include 11 subtypes and 34 circulating recombinant forms. Our knowledge of HIV-1 response to antiretroviral drugs and emergent drug resistance has, however, been limited to subtype B infections circulating in Europe and North America, with little comparative information on non-B subtypes representing approximately 90% of worldwide epidemics. This review summarizes publications in the past year that highlight intersubtype differences influencing viral susceptibility to antiretroviral drugs and emergent drug resistance.

RECENT FINDINGS

Cumulative findings from clinical studies suggest that antiretroviral therapy will be of benefit in the overall treatment of non-B subtype infections, and result in drug-resistance profiles comparable to those observed for subtype B infections. Nevertheless, the 10-15% sequence diversity in the Pol region contributes to intersubtype differences in response to particular nucleoside and non-nucleoside analogues, as well as protease inhibitors. Distinct signature mutations and mutational pathways are identified for specific non-B subtypes. The implications of subtype on clinical outcome and interpretative algorithms are described.

SUMMARY

Understanding intersubtype differences in drug resistance is important in optimizing treatment strategies in resource-poor settings. Hopefully, this may assist in the design of prophylactic approaches to prevent HIV-1 horizontal and vertical HIV-1 transmission.

Near full-length genome characterization of a newly identified HIV type 1 subtype F variant circulating in St. Petersburg, Russia

We report the near full-length genome characterization of an HIV-1 subtype F virus (D88_845) collected in St. Petersburg, Russia, from a 25-year-old Russian woman perinatally infected in 1982. In a Bayesian phylogenetic analysis, the genome sequence branched basally to the subsubtype F1 clade. In partial sequences, D88_845 clustered with 13 other subtype F sequences from Russia, corresponding to gag (n = 2), pol (n = 3), and env (n = 8) segments. At least 11 of these sequences are from samples collected in St. Petersburg from heterosexually infected Russian individuals. In each of these segments, the Russian viruses formed a monophyletic cluster that branched as a sister clade of the F1 subsubtype. One sequence from Belgium branched with D88_845 with a posterior probability of 0.99. This is the first report on the identification and near full-length genome characterization of the subtype F variant circulating in St. Petersburg, which is closely related to, but distinct from, the F1 subsubtype.

Molecular epidemiology of HIV-1 in St Petersburg, Russia: predominance of subtype A, former Soviet Union variant, and identification of intrasubtype subclusters

OBJECTIVES

To examine HIV-1 genetic diversity in St. Petersburg.

METHODS

Partial HIV-1 pol sequences from 102 plasma samples collected in 2006 were analyzed with a Bayesian phylogeny inference method.

RESULTS

Subtype A, former Soviet Union (FSU) variant (AFSU), was the predominant clade (89.3%); other clades were subtypes B (9.7%) and F1 (1%). AFSU was predominant both among injecting drug users (98.2%) and heterosexually infected individuals (91.4%), whereas subtype B was more prevalent among homosexual men (75%). Within the AFSU variant, most sequences (93.5%) branched within 1 of 4 strongly supported subclusters. The largest comprised 63% AFSU viruses and was uncommon outside St Petersburg. A second subcluster (17.4% AFSU viruses) corresponds to the variant with the V77I substitution in protease, which is widely circulating in different FSU countries. Two minor subclusters comprised 8.7% and 6.5% AFSU viruses, respectively. There was no correlation between risk exposure and AFSU subclusters. Six of 8 subtype B sequences, 4 of them from homosexual men, grouped in a monophyletic subcluster.

CONCLUSIONS

The results of this study show a great predominance of AFSU viruses in St Petersburg and point to a few phylogenetically identifiable introductions as the origin of most current HIV-1 AFSU infections in the city.

HIV diversity in Venezuela: predominance of HIV type 1 subtype B and genomic characterization of non-B variants

The aim of this study was the analysis of human immunodeficiency virus (HIV) diversity in Venezuela, and the characterization of variants other than subtype B. A total of 425 HIV isolates, collected between 2003 and 2008, were analyzed. The sequence of at least one genomic region (Pol, Env, Vif, or Nef ) was available for all of them and at least two genomic regions were analyzed in 46% of them. From the 425 HIV isolates analyzed, 421 (99.1%) were classified as HIV-1 subtype B. The four non-subtype B isolates correspond to one subtype C, one recombinant AG, and two HIV-2 isolates. This study shows that HIV-1 subtype B is still highly predominant in Venezuela. Whereas some sporadic cases of other HIV types can be found, they do not seem to have disseminated to the present.

Gene expression analysis from human immunodeficiency virus type 1 subtype C promoter and construction of bicistronic reporter vectors

We report the cloning and sequence analysis of the long terminal repeat (LTR) of several primary HIV-1 subtype C strains of India. Phylogenetically, all the LTRs and the paired env sequences clustered with subtype C reference strains. The LTRs demonstrated extensive polymorphism in the transcription factor binding sites (TFBS) within the enhancer and the modulator regions. We generated reporter vectors under the control of a select subset of the subtype C LTRs. The reporter vectors are distinguished by the simultaneous expression of two independent reporter genes, secreted alkaline phosphatase (SEAP) and enhanced green fluorescence protein (EGFP), in response to Tat. Expression of EGFP was facilitated by engineering an internal ribosome entry site (IRES) into the expression cassette. Although subtype C strains cause a large majority of the global infections, and important differences in the transcription factor binding sites have been identified in the subtype C promoter, few reporter vectors containing subtype C-LTR have been described. We analyzed gene expression from the C-LTR reporter vectors in different cell lines under diverse experimental conditions and compared it to the B-LTR reporter vector. The reporter vectors were responsive to Tat derived from diverse viral subtypes. Furthermore, a positive correlation was observed between the expression of the reporter genes and the viral structural protein p24 when the cells were infected with viral molecular clones. The LTR reporters we developed could be of significant use in the study of viral transactivation, in the evaluation of biological properties of viral subtypes, and in the screening for antiviral inhibitors.

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

BG intersubtype recombinants represented 11.6% of HIV-1 isolates in a recent survey in Cuba based on pol sequences, most of them forming a single clade further subdivided into 3 subclades. Here, we analyze 8 near full-length genomes and 1 gag-pol sequence from epidemiologically unlinked Cuban BG recombinants from these 3 subclades (3 from each). Near full-length sequences were also obtained from 3 subtype G and 2 subtype B Cuban viruses. Phylogenetic relationships were estimated via maximum likelihood, and mosaic structures of the recombinants were inferred with the bootscanning, MaxChi, Genconv, and GARD methods. For the near full-length genomes, all recombinants formed a strongly supported clade further subdivided into the same subclades previously defined in pol. Mosaic structures were identical within each subclade and different among subclades, although 5 breakpoints were coincident among all recombinants. Individual phylogenetic trees for nonrecombinant fragments (concatenated B and G subtype segments) indicated a common ancestry for the parental viruses and their relationships to local subtype B and G strains. These results allow us to identify 3 new BG intersubtype circulating recombinant forms in Cuba derived from a common recombinant ancestor, which originated from B and G subtype parental strains circulating in Cuba.

Genotypic resistance to antiretroviral drugs in patients infected with several HIV type 1 genetic forms in Cuba

The main objective of this study is to evaluate the prevalence of resistance-associated mutations to reverse transcriptase (RT) and protease (PR) inhibitors (I) 2 years after the introduction of antiretroviral treatment in Cuba, analyzing the mutations corresponding to different HIV-1 genetic forms circulating in Cuba. A total of 425 plasma samples were collected in 2003, corresponding to 175 (41.2%) subtype B and 250 (58.8%) non-B genetic forms, including 56 (22.4 %) non-B subtypes, 112 (44.8%) circulating recombinant forms (CRFs), and 82 (32.8%) unique RFs (URFs). Of these, 175 (41.2%) patients were under highly active antiretroviral therapy (HAART) and 250 (58.8%) were treatment-naive. The presence of RT and PR resistance-associated mutations was established by sequencing. Levels of resistance were evaluated according to the Stanford Database program (http://hivdb.stanford.edu). The prevalence of resistance to RTI was 52.2% among RTI-treated patients, 51.5% for subtype B, and 53.2% for non-B genetic forms, including CRF18_cpx, CRF19_cpx, subtype C, and BG URF. In treatment-naive patients it was 6.4% in subtype B and 4.2% in non-B subtypes and RFs. The prevalence of resistance to PRI was 30% among PRI-treated patients, 28% in subtype B and 31% in non-B genetic forms, and 3.2% among treatment-naive subjects, mostly BG recombinants. In conclusion, significant differences in the prevalence of resistance to RTI and PRI were not detected among the most frequent genetic forms from treated patients, suggesting that the genetic diversity of HIV-1 in Cuba does not play a main role in the development of resistance to antiretroviral drugs. The presence of transmitted resistance mutations supports the study of resistance at baseline of treatment.

Antiretroviral drug resistance and phylogenetic diversity of HIV-1 in Chile

This study reports the analysis of human immunodeficiency virus type 1 (HIV-1) protease (PR) and reverse transcriptase (RT) coding sequences from 136 HIV-1-infected subjects from Chile, 66 (49%) of them under antiretroviral (ARV) treatment. The prevalence of mutations conferring high or intermediate resistance levels to ARVs was 77% among treated patients and 2.5% among drug-naïve subjects. The distribution of resistance prevalence in treated patients by drug class was 61% to nucleoside RT inhibitors, 84% to nonnucleoside RT inhibitors, and 46% to PR inhibitors. Phylogenetic analysis revealed that 115 (85%) subjects were infected with subtype B viruses, 1 with a subtype F1 virus, and 20 (15%) carried BF intersubtype recombinants. Most BF recombinants grouped into two clusters, one related to CRF12_BF, while the other could represent a new circulating recombinant form (CRF). In conclusion, this is the first report analysing the prevalence of ARV resistance which includes patients under HAART from Chile. Additionally, phylogenetic analysis of the PR-RT coding sequences reveals the presence of BF intersubtype recombinants.