Analysis of the origin and evolutionary history of HIV-1 CRF28_BF and CRF29_BF reveals a decreasing prevalence in the AIDS epidemic of Brazil

Viruses Previously Identified in Brazil as Belonging to HIV-1 CRF72_BF1 Represent Two Closely Related Circulating Recombinant Forms, One of Which, Designated CRF122_BF1, Is Also Circulating in Spain

Circulating recombinant forms (CRFs) are important components of the HIV-1 pandemic. Those derived from recombination between subtype B and subsubtype F1, with 18 reported, most of them of South American origin, are among the most diverse. In this study, we identified a HIV-1 BF1 recombinant cluster that is expanding in Spain, transmitted mainly via heterosexual contact, which, analyzed in near full-length genomes in four viruses, exhibited a coincident BF1 mosaic structure, with 12 breakpoints, that fully coincided with that of two viruses (10BR_MG003 and 10BR_MG005) from Brazil, previously classified as CRF72_BF1. The three remaining Brazilian viruses (10BR_MG002, 10BR_MG004, and 10BR_MG008) previously identified as CRF72_BF1 exhibited mosaic structures highly similar, but not identical, to that of the Spanish viruses and to 10BR_MG003 and 10BR_MG005, with discrepant subtypes in two short genome segments, located in pol and gp120^env. Based on these results, we propose that the five viruses from Brazil previously identified as CRF72_BF1 actually belong to two closely related CRFs, one comprising 10BR_MG002, 10BR_MG004, and 10BR_MG008, which keep their CRF72_BF1 designation, and the other, designated CRF122_BF1, comprising 10BR_MG003, 10BR_MG005, and the viruses of the identified Spanish cluster. Three other BF1 recombinant genomes, two from Brazil and one from Italy, previously identified as unique recombinant forms, were classified as CRF72_BF1. CRF122_BF1, but not CRF72_BF1, was associated with protease L89M substitution, which was reported to contribute to antiretroviral drug resistance. Phylodynamic analyses estimate the emergence of CRF122_BF1 in Brazil around 1987. Given their close phylogenetic relationship and similar structures, the grouping of CRF72_BF1 and CRF122_BF1 in a CRF family is proposed.

Social, Epidemiological, and Virological Characteristics from Peruvian Subjects Living with HIV-1/AIDS with Different Sexual Risk Behavior

HIV-1 genetic diversity and resistance profile might change according to the risky sexual behavior of the host. To show this, we recruited 134 individuals between the years 2015 and 2017 identified as transgender women sex workers (TWSW, n = 73) and Heterosexual Military Officers (HET-MO, n = 61). After obtaining informed consent, we collected a blood sample to perform the HIV genotyping, CD4 cell count, and viral load. We used bioinformatics approaches for detecting resistance mutations and recombination events. Epidemiological data showed that both groups reported sexually transmitted diseases and they were widespread among TWSW, especially syphilis and herpes virus (35.6%). Illegal drugs consumption was higher among TWSW (71.2%), whereas condom use was inconsistent for both HET-MO (57.4%) and TWSW (74.0%). TWSW showed the shortest time exposition to antiretroviral therapy (ART) (3.5 years) and the lowest access to ART (34.2%) that conducted treatment failure (>4 logs). HIV-1 sequences from TWSW and HET-MO were analyzed to determine the genetic diversity and antiretroviral drug resistance. Phylogeny analysis revealed 125 (93%) cases of subtype B, 01 subtype A (0.76%), 07 (5.30%) BF recombinants, and 01 (0.76%) AG recombinant. Also, TWSW showed a higher recombination index (9.5%, 7/73) than HET-MO (1.5%, 1/68). HET-MO only showed acquired resistance (26.23%, 16/61), whereas TWSW showed both acquired as transmitted resistance (9.59% for each). In conclusion, TWSW and HET-MO showed significant differences considering the epidemiological characteristics, genetic diversity, recombination events, and HIV resistance profile.

Identification of CRF66_BF, a New HIV-1 Circulating Recombinant Form of South American Origin

Circulating recombinant forms (CRFs) are important components of the HIV-1 pandemic. Among 110 reported in the literature, 17 are BF1 intersubtype recombinant, most of which are of South American origin. Among these, all 5 identified in the Southern Cone and neighboring countries, except Brazil, derive from a common recombinant ancestor related to CRF12_BF, which circulates widely in Argentina, as deduced from coincident breakpoints and clustering in phylogenetic trees. In a HIV-1 molecular epidemiological study in Spain, we identified a phylogenetic cluster of 20 samples from 3 separate regions which were of F1 subsubtype, related to the Brazilian strain, in protease-reverse transcriptase (Pr-RT) and of subtype B in integrase. Remarkably, 14 individuals from this cluster (designated BF9) were Paraguayans and only 4 were native Spaniards. HIV-1 transmission was predominantly heterosexual, except for a subcluster of 6 individuals, 5 of which were men who have sex with men. Ten additional database sequences, from Argentina (n = 4), Spain (n = 3), Paraguay (n = 1), Brazil (n = 1), and Italy (n = 1), branched within the BF9 cluster. To determine whether it represents a new CRF, near full-length genome (NFLG) sequences were obtained for 6 viruses from 3 Spanish regions. Bootscan analyses showed a coincident BF1 recombinant structure, with 5 breakpoints, located in p17^gag, integrase, gp120, gp41-rev overlap, and nef, which was identical to that of two BF1 recombinant viruses from Paraguay previously sequenced in NFLGs. Interestingly, none of the breakpoints coincided with those of CRF12_BF. In a maximum likelihood phylogenetic tree, all 8 NFLG sequences grouped in a strongly supported clade segregating from previously identified CRFs and from the CRF12_BF “family” clade. These results allow us to identify a new HIV-1 CRF, designated CRF66_BF. Through a Bayesian coalescent analysis, the most recent common ancestor of CRF66_BF was estimated around 1984 in South America, either in Paraguay or Argentina. Among Pr-RT sequences obtained by us from HIV-1-infected Paraguayans living in Spain, 14 (20.9%) of 67 were of CRF66_BF, suggesting that CRF66_BF may be one of the major HIV-1 genetic forms circulating in Paraguay. CRF66_BF is the first reported non-Brazilian South American HIV-1 CRF_BF unrelated to CRF12_BF.

Geographically-stratified HIV-1 group M pol subtype and circulating recombinant form sequences

Accurate classification of HIV-1 group M lineages, henceforth referred to as subtyping, is essential for understanding global HIV-1 molecular epidemiology. Because most HIV-1 sequencing is done for genotypic resistance testing pol gene, we sought to develop a set of geographically-stratified pol sequences that represent HIV-1 group M sequence diversity. Representative pol sequences differ from representative complete genome sequences because not all CRFs have pol recombination points and because complete genome sequences may not faithfully reflect HIV-1 pol diversity. We developed a software pipeline that compiled 6,034 one-per-person complete HIV-1 pol sequences annotated by country and year belonging to 11 pure subtypes and 70 CRFs and selected a set of sequences whose average distance to the remaining sequences is minimized for each subtype/CRF and country to generate a Geographically-Stratified set of 716 Pol Subtype/CRF (GSPS) reference sequences. We provide extensive data on pol diversity within each subtype/CRF and country combination. The GSPS reference set will also be useful for HIV-1 pol subtyping.

Genomic characterization and molecular evolution analysis of subtype B and BF recombinant HIV-1 strains among Argentinean men who have sex with men reveal a complex scenario

Currently, data on HIV-1 circulating strains among men who have sex with men (MSM) in Argentina is scarce. In South America, the distribution and the prevalence of BF recombinants are dissimilar and exhibit an underappreciated heterogeneity of recombinant structures. Here, we studied for the first time the genetic diversity of HIV-1 BF recombinants and their evolution over time through in-depth phylogenetic analysis and multiple recombination detection methods involving 337 HIV-1 nucleotide sequences (25 near full-length (NFL) and 312 partial pol gene) obtained from Argentinean MSM. The recombination profiles were studied using multiple in silico tools to characterize the genetic mosaicism, and phylogenetic approaches to infer their relationships. The evolutionary history of BF recombinants and subtype B sequences was reconstructed by a Bayesian coalescent-based method. By phylogenetic inference, 81/312 pol sequences clustered within BF clade. Of them, 46 sequences showed a genetic mosaic with CRF12_BF-like patterns, including plausible second-generation recombinants. Other CRFs_BF like (CRF17, 28, 29, 39, 42, 44, 47) and probable URFs_BF were less frequently found. Phylogenetic and recombination analyses on NFL sequences allowed a meticulous definition of new BF mosaics of genomic patterns. The Bayesian analyses pointed out quite consistent onset dates for the CRFs_BF clade based on B and F gene datasets (~1986 and ~1991 respectively). These results indicate that the CRFs_BF variants have been circulating among Argentinean MSM for about 30 years. This study reveals, through growing evidence showing the importance of MSM in the dynamics of the HIV-1 epidemic in Argentina, the coexistence of CRF12_BF-like and high diversity of strains exhibiting several BF mosaic patterns, including non-reported URFs that may reflect active clusters as potential intervention targets to hinder HIV-1 transmission.

Characterization of HIV-1 CRF90_BF1 and putative novel CRFs_BF1 in Central West, North and Northeast Brazilian regions

The Brazilian AIDS epidemic has been characterized by an increasing rate of BF1 recombinants and so far eight circulating recombinant forms/CRFs_BF1 have been described countrywide. In this study, pol sequences (protease/PR, reverse transcriptase/RT) of 87 BF1 mosaic isolates identified among 828 patients living in six Brazilian States from three geographic regions (Central West, North, Northeast) were analyzed. Phylogenetic and bootscan analyses were performed to investigate the evolutionary relationship and mosaic structure of BF1 isolates. Those analyses showed that 20.7% of mosaics (18 out of 87) were CRFs-like isolates, mostly represented by CRF28/CRF29_BF-like viruses (14 out of 18). We also identified five highly supported clusters that together comprise 42 out of 87 (48.3%) BF1 sequences, each cluster containing at least five sequences sharing a similar mosaic structure, suggesting possible new unidentified CRFs_BF1. The divergence time of these five potential new CRFs_BF1 clusters was estimated using a Bayesian approach and indicate that they probably originated between the middle 1980s and the middle 1990s. DNA was extracted from whole blood and four overlapping fragments were amplified by PCR providing full/near full length genomes (FLG/NFLG) and partial genomes. Eleven HIV-1 isolates from Cluster # 5 identified in epidemiologically unlinked individuals living in Central West and North regions provided FLG/NFLG/partial genome sequences with identical mosaic structure. These viruses differ from any known CRF_BF1 reported to date and were named CRF90_BF1 by the Los Alamos National Laboratory. This is the 9^th CRF_BF1 described in Brazil and the first one identified in Central West and North regions. Our results highlight the importance of continued molecular screening and surveillance studies, especially of full genome sequences to understand the evolutionary dynamics of the HIV-1 epidemic in a country of continental dimensions as Brazil.

Assessing the HIV-1 Epidemic in Brazilian Drug Users: A Molecular Epidemiology Approach

Person who inject illicit substances have an important role in HIV-1 blood and sexual transmission and together with person who uses heavy non-injecting drugs may have less than optimal adherence to anti-retroviral treatment and eventually could transmit resistant HIV variants. Unfortunately, molecular biology data on such key population remain fragmentary in most low and middle-income countries. The aim of the present study was to assess HIV infection rates, evaluate HIV-1 genetic diversity, drug resistance, and to identify HIV transmission clusters in heavy drug users (DUs). For this purpose, DUs were recruited in the context of a Respondent-Driven Sampling (RDS) study in different Brazilian cities during 2009. Overall, 2,812 individuals were tested for HIV, and 168 (6%) of them were positive, of which 19 (11.3%) were classified as recent seroconverters, corresponding to an estimated incidence rate of 1.58%/year (95% CI 0.92–2.43%). Neighbor joining phylogenetic trees from env and pol regions and bootscan analyses were employed to subtype the virus from132 HIV-1-infected individuals. HIV-1 subtype B was prevalent in most of the cities under analysis, followed by BF recombinants (9%-35%). HIV-1 subtype C was the most prevalent in Curitiba (46%) and Itajaí (86%) and was also detected in Brasília (9%) and Campo Grande (20%). Pure HIV-1F infections were detected in Rio de Janeiro (9%), Recife (6%), Salvador (6%) and Brasília (9%). Clusters of HIV transmission were assessed by Maximum likelihood analyses and were cross-compared with the RDS network structure. Drug resistance mutations were verified in 12.2% of DUs. Our findings reinforce the importance of the permanent HIV-1 surveillance in distinct Brazilian cities due to viral resistance and increasing subtype heterogeneity all over Brazil, with relevant implications in terms of treatment monitoring, prophylaxis and vaccine development.

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.

Origin and evolution of the unique hepatitis C virus circulating recombinant form 2k/1b

Since its initial identification in St. Petersburg, Russia, the recombinant hepatitis C virus (HCV) 2k/1b has been isolated from several countries throughout Eurasia. The 2k/1b strain is the only recombinant HCV to have spread widely, raising questions about the epidemiological background in which it first appeared. In order to further understand the circumstances by which HCV recombinants might be formed and spread, we estimated the date of the recombination event that generated the 2k/1b strain using a Bayesian phylogenetic approach. Our study incorporates newly isolated 2k/1b strains from Amsterdam, The Netherlands, and has employed a hierarchical Bayesian framework to combine information from different genomic regions. We estimate that 2k/1b originated sometime between 1923 and 1956, substantially before the first detection of the strain in 1999. The timescale and the geographic spread of 2k/1b suggest that it originated in the former Soviet Union at about the time that the world's first centralized national blood transfusion and storage service was being established. We also reconstructed the epidemic history of 2k/1b using coalescent theory-based methods, matching patterns previously reported for other epidemic HCV subtypes. This study demonstrates the practicality of jointly estimating dates of recombination from flanking regions of the breakpoint and further illustrates that rare genetic-exchange events can be particularly informative about the underlying epidemiological processes.

The Use of Bioinformatics for Studying HIV Evolutionary and Epidemiological History in South America

The South American human immunodeficiency virus type 1 (HIV-1) epidemic is driven by several subtypes (B, C, and F1) and circulating and unique recombinant forms derived from those subtypes. Those variants are heterogeneously distributed around the continent in a country-specific manner. Despite some inconsistencies mainly derived from sampling biases and analytical constrains, most of studies carried out in the area agreed in pointing out specificities in the evolutionary dynamics of the circulating HIV-1 lineages. In this paper, we covered the theoretical basis, and the application of bioinformatics methods to reconstruct the HIV spatial-temporal dynamics, unveiling relevant information to understand the origin, geographical dissemination and the current molecular scenario of the HIV epidemic in the continent, particularly in the countries of Southern Cone.