A Founder Effect Led Early SARS-CoV-2 Transmission in Spain.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) whole-genome analysis has identified five large clades worldwide which emerged in 2019 (19A and 19B) and in 2020 (20A, 20B, and 20C). This study aimed to analyze the diffusion of SARS-CoV-2 in Spain using maximum-likelihood phylogenetic and Bayesian phylodynamic analyses. The most recent common ancestor (MRCA) of the SARS-CoV-2 pandemic was estimated to have emerged in Wuhan, China, around 24 November 2019. Phylogenetic analyses of the first 12,511 SARS-CoV-2 whole-genome sequences obtained worldwide, including 290 from 11 different regions of Spain, revealed 62 independent introductions of the virus in the country. Most sequences from Spain were distributed in clades characterized by a D614G substitution in the S gene (20A, 20B, and 20C) and an L84S substitution in ORF8 (19B) with 163 and 118 sequences, respectively, with the remaining sequences branching in 19A. A total of 110 (38%) sequences from Spain grouped in four different monophyletic clusters of clade 20A (20A-Sp1 and 20A-Sp2) and 19B clade (19B-Sp1 and 19B-Sp2) along with sequences from 29 countries worldwide. The MRCAs of clusters 19A-Sp1, 20A-Sp1, 19A-Sp2, and 20A-Sp2 were estimated to have occurred in Spain around 21 and 29 January and 6 and 17 February 2020, respectively. The prevalence of clade 19B in Spain (40%) was by far higher than in any other European country during the first weeks of the epidemic, probably as a result of a founder effect. However, this variant was replaced by G614-bearing viruses in April. In vitro assays showed an enhanced infectivity of pseudotyped virions displaying the G614 substitution compared with those having D614, suggesting a fitness advantage of D614G.IMPORTANCE Multiple SARS-CoV-2 introductions have been detected in Spain, and at least four resulted in the emergence of locally transmitted clusters that originated not later than mid-February, with further dissemination to many other countries around the world, and a few weeks before the explosion of COVID-19 cases detected in Spain during the first week of March. The majority of the earliest variants detected in Spain branched in the clade 19B (D614 viruses), which was the most prevalent clade during the first weeks of March, pointing to a founder effect. However, from mid-March to June 2020, G614-bearing viruses (clades 20A, 20B, and 20C) overcame D614 variants in Spain, probably as a consequence of an evolutionary advantage of this substitution in the spike protein. A higher infectivity of G614-bearing viruses than D614 variants was detected, suggesting that this substitution in SARS-CoV-2 spike protein could be behind the variant shift observed in Spain.

Impact of clade, geography, and age of the epidemic on HIV-1 neutralization by antibodies.

UNLABELLED: Neutralizing antibodies (nAbs) are a high priority for vaccines that aim to prevent the acquisition of HIV-1 infection. Vaccine effectiveness will depend on the extent to which induced antibodies neutralize the global diversity of circulating HIV-1 variants. Using large panels of genetically and geographically diverse HIV-1 Env-pseudotyped viruses and chronic infection plasma samples, we unambiguously show that cross-clade nAb responses are commonly induced in response to infection by any virus clade. Nonetheless, neutralization was significantly greater when the plasma clade matched the clade of the virus being tested. This within-clade advantage was diminished in older, more-diverse epidemics in southern Africa, the United States, and Europe compared to more recent epidemics in Asia. It was most pronounced for circulating recombinant form (CRF) 07_BC, which is common in China and is the least-divergent lineage studied; this was followed by the slightly more diverse Asian CRF01_AE. We found no evidence that transmitted/founder viruses are generally more susceptible to neutralization and are therefore easier targets for vaccination than chronic viruses. Features of the gp120 V1V2 loop, in particular, length, net charge, and number of N-linked glycans, were associated with Env susceptibility and plasma neutralization potency in a manner consistent with neutralization escape being a force that drives viral diversification and plasma neutralization breadth. The overall susceptibility of Envs and potencies of plasma samples were highly predictive of the neutralization outcome of any single virus-plasma combination. These findings highlight important considerations for the design and testing of candidate HIV-1 vaccines that aim to elicit effective nAbs. IMPORTANCE: An effective HIV-1 vaccine will need to overcome the extraordinary variability of the virus, which is most pronounced in the envelope glycoproteins (Env), which are the sole targets for neutralizing antibodies (nAbs). Distinct genetic lineages, or clades, of HIV-1 occur in different locales that may require special consideration when designing and testing vaccines candidates. We show that nAb responses to HIV-1 infection are generally active across clades but are most potent within clades. Because effective vaccine-induced nAbs are likely to share these properties, optimal coverage of a particular clade or combination of clades may require clade-matched immunogens. Optimal within-clade coverage might be easier to achieve in regions such as China and Thailand, where the epidemic is more recent and the virus less diverse than in southern Africa, the United States, and Europe. Finally, features of the first and second hypervariable regions of gp120 (V1V2) may be critical for optimal vaccine design.

Predominance of CXCR4 tropism in HIV-1 CRF14_BG strains from newly diagnosed infections.

OBJECTIVES: R5-tropic viruses are associated with HIV-1 transmission and predominate during the early stages of infection. X4-tropic populations have been detected in ~50% of patients with late-stage disease infected with subtype B viruses. In this study, we compared the frequency of X4 tropism in individuals infected with HIV-1 CRF14_BG viruses, which have a V3 loop of subtype B, with a control group of individuals infected with subtype B viruses. METHODS: Sixty-three individuals infected with HIV-1 CRF14_BG (n = 31) or subtype B (n = 32) were studied. Similar proportions of newly diagnosed and chronically infected individuals were included in the subtype B and CRF14_BG groups. V3 sequences were obtained and coreceptor tropism was predicted using the Geno2pheno[coreceptor] algorithm. V3 net charge and 11/25 rules were also used for coreceptor prediction. RESULTS: Overall, X4 tropism was more frequent among individuals infected with CRF14_BG viruses (87.1%) than subtype B viruses (34.3%), a difference that was statistically highly significant (P = 0.00001). Importantly, the frequencies among newly diagnosed individuals were 90% and 13.3%, respectively (P = 0.0007). Characteristic amino acids in the V3 loop (T13, M14, V19 and W20) were identified at higher frequencies in CRF14_BG viruses (54%) than subtype B viruses (0%; P < 0.000001). CONCLUSIONS: CRF14_BG is the genetic form with the highest proportion of X4-tropic viruses reported to date in newly diagnosed and chronic infections. This suggests high pathogenicity for CRF14_BG viruses, potentially leading to rapid disease progression. CCR5 antagonists will be ineffective in most CRF14_BG-infected patients, even at early stages of infection.

Improvement of HIV-1 coreceptor tropism prediction by employing selected nucleotide positions of the env gene in a Bayesian network classifier.

OBJECTIVES: This study aimed to develop a genotypic method to predict HIV-1 coreceptor usage by employing the nucleotide sequence of the env gene in a tree-augmented naive Bayes (TAN) classifier, and to evaluate its accuracy in prediction compared with other available tools. METHODS: A wrapper data-mining strategy interleaved with a TAN algorithm was employed to evaluate the predictor value of every single-nucleotide position throughout the HIV-1 env gene. Based on these results, different nucleotide positions were selected to develop a TAN classifier, which was employed to predict the coreceptor tropism of all the full-length env gene sequences with information on coreceptor tropism currently available at the Los Alamos HIV Sequence Database. RESULTS: Employing 26 nucleotide positions in the TAN classifier, an accuracy of 95.6%, a specificity (identification of CCR5-tropic viruses) of 99.4% and a sensitivity (identification of CXCR4/dual-tropic viruses) of 80.5% were achieved for the in silico cross-validation. Compared with the phenotypic determination of coreceptor usage, the TAN algorithm achieved more accurate predictions than WebPSSM and Geno2pheno [coreceptor] (P<0.05). CONCLUSIONS: The use of the methodology presented in this work constitutes a robust strategy to identify genetic patterns throughout the HIV-1 env gene differently present in CCR5-tropic and CXCR4/dual-tropic viruses. Moreover, the TAN classifier can be used as a genotypic tool to predict the coreceptor usage of HIV-1 isolates reaching more accurate predictions than with other widely used genotypic tools. The use of this algorithm could improve the correct prescribing of CCR5 antagonist drugs to HIV-1-infected patients.

Identification of a HIV-1 circulating BF1 recombinant form (CRF75_BF1) of Brazilian origin that also circulates in Southwestern Europe.

Introduction: The high recombinogenic potential of HIV-1 has resulted in the generation of countless unique recombinant forms (URFs) and around 120 reported circulating recombinant forms (CRFs). Here we identify through analyses of near full-length genomes (NFLG) a new HIV-1 CRF derived from subtypes B and F1. Methods: HIV-1 protease-reverse transcriptase (Pr-RT) sequences were obtained by RT-PCR amplification from plasma RNA. Near full-length genome sequences were obtained after amplification by RT-PCR in 5 overlapping fragments. Phylogenetic sequence analyses were performed via maximum likelihood. Mosaic structures were analyzed by bootscanning and phylogenetic analyses of genome segments. Temporal and geographical estimations of clade emergence were performed with a Bayesian coalescent method. Results: Through phylogenetic analyses of HIV-1 Pr-RT sequences obtained by us from samples collected in Spain and downloaded from databases, we identified a BF1 recombinant cluster segregating from previously reported CRFs comprising 52 viruses, most from Brazil (n = 26), Spain (n = 11), and Italy (n = 9). The analyses of NFLG genomes of 4 viruses of the cluster, 2 from Spain and 2 from Italy, allowed to identify a new CRF, designated CRF75_BF1, which exhibits a complex mosaic structure with 20 breakpoints. All 4 patients harboring CRF75_BF1 viruses studied by us had CD4+ T-cell lymphocyte counts below 220/mm3 less than one year after diagnosis, a proportion significantly higher (p = 0.0074) than the 29% found in other patients studied in Spain by us during the same period. The origin of the clade comprising CRF75_BF1 and related viruses was estimated around 1984 in Brazil, with subsequent introduction of CRF75_BF1 in Italy around 1992, and migration from Italy to Spain around 1999. Conclusion: A new HIV-1 CRF, designated CRF75_BF1, has been identified. CRF75_BF1 is the 6th CRF of South American origin initially identified in Western Europe, reflecting the increasing relationship of South American and European HIV-1 epidemics. The finding of low CD4+ T-cell lymphocyte counts early after diagnosis in patients harboring CRF75_BF1 viruses warrants further investigation on the virulence of this variant.

Evaluation of genotypic tropism prediction tests compared with in vitro co-receptor usage in HIV-1 primary isolates of diverse subtypes.

OBJECTIVES: To evaluate the sensitivity and specificity of genotypic methods for predicting the co-receptor usage of subtypes B and non-B HIV-1 primary isolates, using as gold standard the infectivity of each primary isolate in GHOST cells stably expressing HIV-1 co-receptors. METHODS: Primary isolates were obtained by co-culturing either patient's peripheral blood mononuclear cells (PBMCs) or ultracentrifuged plasma with donor-activated PBMCs. In vitro co-receptor usage was determined by infecting GHOST cells. Tropism prediction, based on V3 sequences, was determined with simple rules and bioinformatic tools (Geno2pheno[coreceptor] and WebPSSM). RESULTS: This study includes 102 HIV-1 primary isolates; 23 (22.5%) subtype B and 79 (77.5%) non-B genetic forms. V3 sequences were classified into six subtypes (A-G), although 32 (31.4%) were circulating recombinant forms and 21 (20.6%) were unique recombinant forms. Sixty-nine isolates were R5, 27 R5X4 and 6 X4. The highest levels of sensitivity and specificity for the detection of X4 strains among V3 sequences, between 91% and 100%, were obtained by using PSSM(x4r5), PSSM(si/nsi) and the 11/25 rule for sequences of subtypes A, B and G, but not for subtype F. Establishing the recommended cut-off for clinical settings of a 10% false positive rate for Geno2pheno, we obtained 93% specificity and 97% sensitivity. CONCLUSIONS: Comparing genotypic assays for HIV-1 co-receptor use with a cell-culture phenotypic assay could provide more reliable results of sensitivity and specificity for the detection of X4 strains than comparing them with recombinant assays, considered as gold standard. In general, except for subtype F isolates, there is a good correlation for tropism prediction.

Factors associated with HIV-1 resistance to integrase strand transfer inhibitors in Spain: Implications for dolutegravir-containing regimens.

Integrase strand transfer inhibitor (INSTI)-containing regimens in HIV-1-infected patients have experienced a global increase. Recently, WHO has emphasized the need to fast-track the transition to dolutegravir (DTG)-based antiretroviral (ARV) treatments. However, continued surveillance of INSTI resistance is recommended. In this study, clinical, epidemiological, and virological features associated with INSTI resistance diagnosed in Spain were analyzed. Samples collected between 2008 and 2021 from HIV-1-infected patients were analyzed in integrase, protease, and reverse transcriptase using Sanger population sequencing. ARV drug resistance was evaluated with the Stanford University HIVdb program. Among 2,696 patients, 174 (6.5%) had INSTI resistance, all of them to first-generation INSTIs, and 71 (2.6%) had also resistance to second-generation INSTIs. Of these, only 5 individuals were exposed to DTG as the only INSTI, in whom resistance development was associated with poor treatment adherence and/or resistance to other ARV classes. Of newly HIV-1-diagnosed individuals, 0.92% harbored INSTI-resistant viruses, with low prevalences maintained along time, and only one had low-level resistance to DTG. Persons who inject drugs, age over 39 years, resistance to other ARV classes, and longer time from diagnosis were associated with INSTI resistance (p < 0.001). Non-subtype B INSTI-resistant viruses lacked the Q148H + G140S resistance pathway and showed lower INSTI resistance levels than subtype B viruses. In conclusion, INSTI resistance is uncommon and associated with long-term infections, older age and additional resistance to other ARV drug classes, and is rare in newly diagnosed HIV-1 infections. Our results also support the preferential use of DTG-containing regimens in first-line treatments, although surveillance of INSTI resistance is encouraged.

Transmission Clusters, Predominantly Associated With Men Who Have Sex With Men, Play a Main Role in the Propagation of HIV-1 in Northern Spain (2013-2018).

Viruses of HIV-1-infected individuals whose transmission is related group phylogenetically in transmission clusters (TCs). The study of the phylogenetic relations of these viruses and the factors associated with these individuals is essential to analyze the HIV-1 epidemic. In this study, we examine the role of TCs in the epidemiology of HIV-1 infection in Galicia and the Basque County, two regions of northern Spain. A total of 1,158 HIV-1-infected patients from both regions with new diagnoses (NDs) in 2013-2018 were included in the study. Partial HIV-1 pol sequences were analyzed phylogenetically by approximately maximum-likelihood with FastTree 2. In this analysis, 10,687 additional sequences from samples from HIV-1-infected individuals collected in Spain in 1999-2019 were also included to assign TC membership and to determine TCs' sizes. TCs were defined as those which included viruses from ≥4 individuals, at least 50% of them Spaniards, and with ≥0.95 Shimodaira-Hasegawa-like node support in the phylogenetic tree. Factors associated to TCs were evaluated using odds ratios (OR) and their 95% CI. Fifty-one percent of NDs grouped in 162 TCs. Male patients (OR: 2.6; 95% CI: 1.5-4.7) and men having sex with men (MSM; OR: 2.1; 95% CI: 1.4-3.2) had higher odds of belonging to a TC compared to female and heterosexual patients, respectively. Individuals from Latin America (OR: 0.3; 95% CI: 0.2-0.4), North Africa (OR: 0.4; 95% CI: 0.2-1.0), and especially Sub-Saharan Africa (OR: 0.02; 95% CI: 0.003-0.2) were inversely associated to belonging to TCs compared to native Spaniards. Our results show that TCs are important components of the HIV-1 epidemics in the two Spanish regions studied, where transmission between MSM is predominant. The majority of migrants were infected with viruses not belonging to TCs that expand in Spain. Molecular epidemiology is essential to identify local peculiarities of HIV-1 propagation. The early detection of TCs and prevention of their expansion, implementing effective control measures, could reduce HIV-1 infections.

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 gp120env. 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.

The Origin, Epidemiology, and Phylodynamics of Human Immunodeficiency Virus Type 1 CRF47_BF.

CRF47_BF is a circulating recombinant form (CRF) of the human immunodeficiency virus type 1 (HIV-1), the etiological agent of AIDS. CRF47_BF represents one of 19 CRFx_BFs and has a geographic focus in Spain, where it was first identified in 2010. Since its discovery, CRF47_BF has expanded considerably in Spain, predominantly through heterosexual contact (∼56% of the infections). Little is known, however, about the origin and diversity of this CRF or its epidemiological correlates, as very few samples have been available so far. This study conducts a phylogenetic analysis with representatives of all CRFx_BF sequence types along with HIV-1 M Group subtypes to validate that the CRF47_BF sequences share a unique evolutionary history. The CRFx_BF sequences cluster into a single, not well supported, clade that includes their dominant parent subtypes (B and F). This clade also includes subtype D and excludes sub-subtype F2. However, the CRF47_BF sequences all share a most recent common ancestor. Further analysis of this clade couples CRF47_BF protease-reverse transcriptase sequences and epidemiological data from an additional 87 samples collected throughout Spain, as well as additional CRF47_BF database sequences from Brazil and Spain to investigate the origin and phylodynamics of CRF47_BF. The Spanish region with the highest proportion of CRF47_BF samples in the data set was the Basque Country (43.7%) with Navarre next highest at 19.5%. We include in our analysis epidemiological data on host sex, mode of transmission, time of collection, and geographic region. The phylodynamic analysis indicates that CRF47_BF originated in Brazil around 1999-2000 and spread to Spain from Brazil in 2002-2003. The virus spread rapidly throughout Spain with an increase in population size from 2011 to 2015 and leveling off more recently. Three strongly supported clusters associated with Spanish regions (Basque Country, Navarre, and Aragon), together comprising 60.8% of the Spanish samples, were identified, one of which was also associated with transmission among men who have sex with men. The expansion in Spain of CRF47_BF, together with that of other CRFs and subtype variants of South American origin, previously reported, reflects the increasing relationship between the South American and European HIV-1 epidemics.

Genetic Diversity and Drug Resistance Mutations in Reverse Transcriptase and Protease Genes of HIV-1 Isolates from Southwestern Siberia.

The analysis of a pol gene fragment encoding protease and part of reverse transcriptase was carried out for 55 sera collected in 2016 and 2018 from HIV-1-infected patients diagnosed in 2014-2018 living in the south of Western Siberia, Russia: Altai Territory (n = 11), Republic of Altai (n = 15), Kemerovo region (n = 18), and Novosibirsk region (n = 11). CRF63_02A was the dominant genetic form (>70%) in the Altai Territory and Kemerovo and Novosibirsk regions, with subsubtype A6 comprising <30% of samples. In the Altai Republic, subsubtype A6 was predominant (53%), with 33% of viruses belonging to CRF63_02A. Four CRF63_02A/A6 unique recombinant forms were identified in the Altai Territory, Kemerovo Region, and the Altai Republic. A majority (11 of 15) of CRF63_02A viruses from Kemerovo were grouped in a cluster. Antiretroviral (ARV) drug resistance mutations were found in 6 (14%) of 43 drug-naive patients. This study provides new insights in HIV-1 molecular epidemiology and prevalence of transmitted ARV drug resistance mutations in Southwestern Siberia.

Identification of a New HIV-1 BC Intersubtype Circulating Recombinant Form (CRF108_BC) in Spain.

The extraordinary genetic variability of human immunodeficiency virus type 1 (HIV-1) group M has led to the identification of 10 subtypes, 102 circulating recombinant forms (CRFs) and numerous unique recombinant forms. Among CRFs, 11 derived from subtypes B and C have been identified in China, Brazil, and Italy. Here we identify a new HIV-1 CRF_BC in Northern Spain. Originally, a phylogenetic cluster of 15 viruses of subtype C in protease-reverse transcriptase was identified in an HIV-1 molecular surveillance study in Spain, most of them from individuals from the Basque Country and heterosexually transmitted. Analyses of near full-length genome sequences from six viruses from three cities revealed that they were BC recombinant with coincident mosaic structures different from known CRFs. This allowed the definition of a new HIV-1 CRF designated CRF108_BC, whose genome is predominantly of subtype C, with four short subtype B fragments. Phylogenetic analyses with database sequences supported a Brazilian ancestry of the parental subtype C strain. Coalescent Bayesian analyses estimated the most recent common ancestor of CRF108_BC in the city of Vitoria, Basque Country, around 2000. CRF108_BC is the first CRF_BC identified in Spain and the second in Europe, after CRF60_BC, both phylogenetically related to Brazilian subtype C strains.

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.

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.

Diverse Large HIV-1 Non-subtype B Clusters Are Spreading Among Men Who Have Sex With Men in Spain.

In Western Europe, the HIV-1 epidemic among men who have sex with men (MSM) is dominated by subtype B. However, recently, other genetic forms have been reported to circulate in this population, as evidenced by their grouping in clusters predominantly comprising European individuals. Here we describe four large HIV-1 non-subtype B clusters spreading among MSM in Spain. Samples were collected in 9 regions. A pol fragment was amplified from plasma RNA or blood-extracted DNA. Phylogenetic analyses were performed via maximum likelihood, including database sequences of the same genetic forms as the identified clusters. Times and locations of the most recent common ancestors (MRCA) of clusters were estimated with a Bayesian method. Five large non-subtype B clusters associated with MSM were identified. The largest one, of F1 subtype, was reported previously. The other four were of CRF02_AG (CRF02_1; n = 115) and subtypes A1 (A1_1; n = 66), F1 (F1_3; n = 36), and C (C_7; n = 17). Most individuals belonging to them had been diagnosed of HIV-1 infection in the last 10 years. Each cluster comprised viruses from 3 to 8 Spanish regions and also comprised or was related to viruses from other countries: CRF02_1 comprised a Japanese subcluster and viruses from 8 other countries from Western Europe, Asia, and South America; A1_1 comprised viruses from Portugal, United Kingom, and United States, and was related to the A1 strain circulating in Greece, Albania and Cyprus; F1_3 was related to viruses from Romania; and C_7 comprised viruses from Portugal and was related to a virus from Mozambique. A subcluster within CRF02_1 was associated with heterosexual transmission. Near full-length genomes of each cluster were of uniform genetic form. Times of MRCAs of CRF02_1, A1_1, F1_3, and C_7 were estimated around 1986, 1989, 2013, and 1983, respectively. MRCA locations for CRF02_1 and A1_1 were uncertain (however initial expansions in Spain in Madrid and Vigo, respectively, were estimated) and were most probable in Bilbao, Spain, for F1_3 and Portugal for C_7. These results show that the HIV-1 epidemic among MSM in Spain is becoming increasingly diverse through the expansion of diverse non-subtype B clusters, comprising or related to viruses circulating in other countries.

IL-7 is a potent and proviral strain-specific inducer of latent HIV-1 cellular reservoirs of infected individuals on virally suppressive HAART.

The persistence of HIV-1 in virally suppressed infected individuals on highly active antiretroviral therapy (HAART) remains a major therapeutic problem. The use of cytokines has been envisioned as an additional therapeutic strategy to stimulate latent proviruses in these individuals. Immune activation therapy using IL-2 has shown some promise. In the present study, we found that IL-7 was significantly more effective at enhancing HIV-1 proviral reactivation than either IL-2 alone or IL-2 combined with phytohemagglutinin (PHA) in CD8-depleted PBMCs. IL-7 also showed a positive trend for inducing proviral reactivation from resting CD4(+) T lymphocytes from HIV-1-infected patients on suppressive HAART. Moreover, the phylogenetic analyses of viral envelope gp120 genes from induced viruses indicated that distinct proviral quasispecies had been activated by IL-7, as compared with those activated by the PHA/IL-2 treatment. These studies thus demonstrate that different activators of proviral latency may perturb and potentially deplete only selected, specific portions of the proviral archive in virally suppressed individuals. The known immunomodulatory effects of IL-7 could be combined with its ability to stimulate HIV-1 replication from resting CD4(+) T lymphocytes, in addition to other moieties, to potentially deplete HIV-1 reservoirs and lead to the rational design of immune-antiretroviral approaches.

Genome-scale analysis of evolutionary rate and selection in a fast-expanding Spanish cluster of HIV-1 subtype F1.

This work is aimed at assessing the presence of positive selection and/or shifts of the evolutionary rate in a fast-expanding HIV-1 subtype F1 transmission cluster affecting men who have sex with men in Spain. We applied Bayesian coalescent phylogenetics and selection analyses to 23 full-coding region sequences from patients belonging to that cluster, along with other 19 F1 epidemiologically-unrelated sequences. A shift in the overall evolutionary rate of the virus, explained by positively selected sites in the cluster, was detected. We also found one substitution in Nef (H89F) that was specific to the cluster and experienced positive selection. These results suggest that fast transmission could have been facilitated by some inherent genetic properties of this HIV-1 variant.

HIV-1 Genetic Diversity in Recently Diagnosed Infections in Moscow: Predominance of AFSU, Frequent Branching in Clusters, and Circulation of the Iberian Subtype G Variant.

HIV-1 protease-reverse transcriptase sequences from 62 HIV-1-infected individuals recently diagnosed in Moscow were analyzed. Subtype A former Soviet Union (FSU) (AFSU) variant was the predominant clade (62.9%), followed by subtype B (22.6%), unique recombinants (6.5%), subtype G (6.5%), and CRF01_AE (1.6%). AFSU predominated among people who inject drugs (88.9%) and heterosexually acquired infections (77.8%), while subtype B was the most prevalent genetic form among men who have sex with men (44%), although AFSU was also frequent in this population (36%). Forty-eight (77.4%) viruses branched within intrasubtype clusters, three of which, of subtype B, had a majority of viruses collected outside of FSU. The four subtype G viruses identified in this study belonged to the Portuguese-Spanish (Iberian) variant and, together with three from databases, formed a Russian cluster closely related to viruses from Denmark. This is the first report of the circulation of the Iberian subtype G variant in Russia.

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.

New insights into the origin of the HIV type 1 subtype A epidemic in former Soviet Union's countries derived from sequence analyses of preepidemically transmitted viruses.

The HIV-1 subtype A epidemic affecting injecting drug users (IDU) in former Soviet Union (FSU) countries started dramatically in Odessa, southern Ukraine, in 1995, and is caused by a variant of monophyletic origin, often designated IDU-A. We phylogenetically analyzed one near full-length genome and two partial sequences of three HIV-1 subtype A viruses collected in St. Petersburg, Russia, heterosexually transmitted in 1992-1994. The sequences branched basally to the IDU-A clade, together with eight viruses from Odessa collected in 1993, all presumably acquired heterosexually, and two viruses from the Democratic Republic of Congo. Of all other FSU sequences in databases, only those from three recently collected viruses, one from Ukraine and two from northwestern Russia, at least one of them acquired heterosexually, branched basally to the IDU-A cluster. The results indicate that the FSU IDU-A variant derives from a strain that initially propagated heterosexually in Ukraine and originated in central Africa.