Get 20M+ Full-Text Papers For Less Than $1.50/day. Subscribe now for You or Your Team.

Learn More →

High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain: A phylogeographic approach

High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain:... OPENACCESS Phylogenetic studies are a valuable tool to understand viral transmission patterns and the Citation: PeÂrez-Parra S, Chueca N, Alvarez M, role of immigration in HIV-1 spread. We analyzed the spatio-temporal relationship of differ- Pasquau J, Omar M, Collado A, et al. (2017) High ent HIV-1 non-B subtype variants over time using phylogenetic analysis techniques. We col- prevalence and diversity of HIV-1 non-B genetic lected 693 pol (PR+RT) sequences that were sampled from 2005 to 2012 from naïve forms due to immigration in southern Spain: A patients in different hospitals in southern Spain. We used REGA v3.0 to classify them into phylogeographic approach. PLoS ONE 12(10): e0186928. https://doi.org/10.1371/journal. subtypes and recombinant forms, which were confirmed by phylogenetic analysis through pone.0186928 maximum likelihood (ML) using RAxML. For the main HIV-1 non-B variants, publicly avail- Editor: Jason Blackard, University of Cincinnati able, genetically similar sequences were sought using HIV-BLAST. The presence of HIV-1 College of Medicine, UNITED STATES lineages circulating in our study population was established using ML and Bayesian infer- Received: May 1, 2017 ence (BEAST v1.7.5) and transmission networks were identified. We detected 165 (23.4%) patients infected with HIV-1 non-B variants: 104 (63%) with recombinant viruses in pol: Accepted: October 10, 2017 CRF02_AG (71, 43%), CRF14_BG (8, 4.8%), CRF06_cpx (5, 3%) and nine other recombi- Published: October 30, 2017 nant forms (11, 6.7%) and unique recombinants (9, 5.5%). The rest (61, 37%) were infected Copyright:© 2017 Perez-Parra et al. This is an with non-recombinant subtypes: A1 (30, 18.2%), C (7, [4.2%]), D (3, [1.8%]), F1 (9, 5.5%) open access article distributed under the terms of and G (12, 7.3%). Most patients infected with HIV-1 non-B variants were men (63%, p< the Creative Commons Attribution License, which permits unrestricted use, distribution, and 0.001) aged over 35 (73.5%, p< 0.001), heterosexuals (92.2%, p< 0.001), from Africa reproduction in any medium, provided the original (59.5%, p< 0.001) and living in the El Ejido area (62.4%, p<0.001). We found lineages of author and source are credited. epidemiological relevance (mainly within Subtype A1), imported primarily through female Data Availability Statement: All relevant data are sex workers from East Europe. We detected 11 transmission clusters of HIV-1 non-B Sub- within the paper and its Supporting Information types, which included patients born in Spain in half of them. We present the phylogenetic files. profiles of the HIV-1 non-B variants detected in southern Spain, and explore their putative Funding: This work was supported by National geographical origins. Our data reveals a high HIV-1 genetic diversity likely due to the import Plan of I+D+i (RD12/0017/006), granted to FGG; of viral lineages that circulate in other countries. The highly immigrated El Ejido area acts as General Ministry of Evaluation and European, granted to FGG; Fund of Regional Development a gateway through which different subtypes are introduced into other regions, hence the (FEDER), granted to FGG; Health Institute Carlos III importance of setting up epidemiological control measures to prevent future outbreaks. (ISCIII), PI12/01053, granted to FGG; and Ministry of Health of Andalusia (AC-0082-2013), granted to PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 1 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain FGG. The funders had no role in study design, data Introduction collection and analysis, decision to publish, or The high evolutionary rate and recombination capacity of human immunodeficiency virus preparation of the manuscript. type 1 (HIV-1) determine the existence of an array of subtypes and recombinant forms circu- Competing interests: The authors have declared lating worldwide [1±4]. HIV-1 non-subtype B (ªnon-Bº) variants cause around 90% of infec- that no competing interests exist. tions worldwide, and largely predominate in African or Eastern European countries with generalized HIV-1 epidemics. Subtypes C and A, and circulating recombinant forms (CRF) CRF01_AE and CRF02_AG, are responsible alone for 70% of the world's infections [5]. Nowa- days the proportion of infections by HIV-1 non-B variants in Spain lies at 12±15%, depending on the study and technique used to characterize variants [6,7]. Nonetheless, the predominance of HIV-1 subtype B in developed countries (where antiretroviral therapy is more widespread), implies that this is the most widely studied subtype from the genetic, biological and therapeutic viewpoints. The full biological meaning of the genetic variability of HIV-1 is still not completely understood. However, several major differences between the biological properties of certain genetic subtypes in have been described; e.g., virulence, tropism and transmissibility [8,9], use of chemokine co-receptors [10], disease progression [11], susceptibility to some anti- retroviral drugs [12,13], sensitivity to viral load quantification methods [14,15] and detection [16]. These findings evidence the importance of epidemiological information about different subtypes. Eastern Andalusia is located in south-eastern Spain, and includes the provinces of Almerõ Âa, Granada and Jae Ân. Given its location and geographic closeness to the African continent, this region has received a notable foreign migratory influx in the last decade. Andalusia is the fourth Spanish Autonomous Community in number of foreign population, only surpassed by Catalonia, Madrid and the Valencian Community. The main source of immigration in Eastern Andalusia stems from its intensive farming practices, mainly in the El Ejido area (located in the province of Almerõ Âa), where one in every four citizens is an immigrant. Phylogenetic analyses, in conjunction with geographical data, can assess the existing rela- tionship between migratory events and spread of HIV-1 on a local scale [17±20], and to study HIV-1 transmission networks locally [21±24]. As in previous studies [25], our center collects the HIV-1 pol gene sequences linked to the patients' clinical data to monitor baseline drug resistance in naïve individuals in Eastern Andalusia. Our aims were to describe the molecular epidemiology and evolutionary history of non-B forms in Eastern Andalusia over the 2005± 2012 period, and to explore their putative geographical origin prior to their arrival to our region. Methods Study population During the study period (2005±2012), 693 pol gene sequences of patients newly diagnosed with HIV-1 in different Eastern Andalusian hospitals were collected from routine drug resis- tance analyses. These hospitals were distributed in 3 provinces: Granada (which included its    capital city of Granada and Motril), Jaen, and Almerõa (including its capital city of Almerõa and El Ejido). The pol sequences (protease (PR), codons 4±99; reverse transcriptase (RT), codons 38±247) obtained by the Trugene1 HIV Genotyping kit (Siemens, NAD), were linked to demographic (risk group, age, gender, country of origin, sampling year, and attending hos- pital), clinical (CD4+ T-cell count) and virological (plasma viral load) information. Demo- graphic information was voluntarily collected during clinical interviews. This study was approved by the San Cecilio Hospital's Ethics Committee, and no consent information was required as patient information remained anonymous and was de-identified prior to analyses. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 2 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain HIV-1 pol sequencing and subtype assignment All the sequences were trimmed to 883 nucleotides (nt) and aligned using ClustalW [26]. The viral subtype was studied with the REGA v3.0 subtyping tool (http://dbpartners.stanford. edu:8080/RegaSubtyping/stanford-hiv/typingtool/), and was confirmed by phylogenetic analy- sis through maximum likelihood (ML) using the randomized Accelerated Maximum Likeli- hood (RAxML) program, accessible on the CIPRES Science Gateway [27]. The general time- reversible (GTR) model with a gamma-distributed heterogeneity rate across sites was employed, applying 1000 bootstrap iterations. A representative dataset of HIV-1 group M sequences, including non- recombinant subtypes (A-K) and recombinant forms (at least four representative sequences of each non-recombinant subtype and the CRF currently available from the analysis) were downloaded from the Los Alamos HIV sequence database (http:// www.hiv.lanl.gov) was used as a reference dataset (S1 Table). The assignment to any subtype/CRF was considered definitive if the query sequence was included with the reference sequences corresponding to that viral variant in a monophyletic cluster supported by high bootstrap values (>70%) [28]. Any genetic form not associated with reference subtypes/CRFs was classified as a unique recombinant form (URF), whose recombi- nation pattern was further studied by a Bootscan analysis using the SimPlot v3.5.1 software [29]. The bootscanning method in SimPlot consists of a sliding-window phylogenetic boot- strap analysis of the query sequence aligned against a set of reference strains to reveal break- points. The Neighbor-Joining algorithm was selected, with the Kimura 2-parameter substitution model. We employed a window size of 200nt moving in 10nt increments. We used a minimum cutoff for the bootstrap value of 70% to reliably assign each of the breakpoint segments to a parental variant. We have submitted to GenBank the major groups of HIV-1 non-B variants under accession numbers MF628109 to MF628250. These were defined as those found in at least five patients. With the aim of protecting the identity of patients infected with rare genetic forms of HIV-1, and for similar scientific and ethical reasons as explained in other HIV cohorts [30±32], we decided not to submit to GenBank those sequences corresponding to the less frequent variants. Inference of the putative geographical origins of the HIV-1 non-B variants circulating in Andalusia To further characterize the relationships among the major groups of HIV-1 non-B variants, we interrogated GenBank for genetically related sequences to our major subtypes/recombinant forms using HIV-BLAST (http://www.hiv.lanl.gov/content/sequence/BASIC_BLAST/basic_ blast.html). The 10 most closely related GenBank sequences to each of our study sequences, were downloaded and included in each dataset. We also included all the pol sequences (start: 2293 and end: 3290, HXB2 coordinates), available in the HIV Los Alamos database sampled in Spain for each dataset: subtype A1 (n = 60), subtype C (n = 52), subtype F (n = 143), subtype G (n = 64), CRF14_BG (n = 25), and CRF02_AG (n = 265). Since very few sequences for CRF06_cpx were available in public databases (http://www.hiv.lanl.gov/content/sequence/ HIV/mainpage.html), we included them all (n = 110). All these individual sequence datasets were put together (n = 970) and a global phylogenetic analysis was performed using RAxML (GTR + Gamma model) and 1000 bootstrap iterations for this analysis. The phylogenetic relatedness between the sequences was studied, and a 70% bootstrap value was taken as a significantly reliable value [28]. Thresholds for low genetic dis- tance, which are commonly used as a proxy for divergence time, were not applied to the cluster definition in the ML trees since these clusters were further confirmed and analyzed using a PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 3 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain time-stamped Bayesian phylogenetic analysis with BEAST, as described below. International non-B lineages (defined as phylogenetic associations of at least one sequence from our cohort clustered with sequences from different countries), and `Andalusian clusters' (monophyletic associations of sequences in our cohort alone), were identified in the global ML tree. A Bayesian Markov Chain Monte Carlo (MCMC) approach was applied to each of the indi- vidual HIV-1 non-B subtype/CRFs datasets described above, which included the most geneti- cally similar sequences found with HIV-1 BLAST, as implemented in BEAST v1.7.5 [33]. The Shapiro-Rambaut-Drummond-2006 (SRD06) substitution model was used, together with a relaxed uncorrelated lognormal clock (UCLN)[34] and a demographic non parametric model, Bayesian Skyline Plot (BSP) [35]. This model combination was chosen because it best fits the analysis of the HIV-1 pol data run in the majority of studies [36]. The MCMC was run for 250 million states sampling every 50000. The evolutionary rate (μ, nucleotide substitutions per site per year, subst./site/year) for the different HIV-1 non-B subtypes/CRFs (S2 Table), and the most recent common ancestors (MRCA) of the different HIV-1 non-B clusters, were esti- mated. Only traces with an effective sample size (ESS) > 200 for all the parameters, after excluding an initial 10% burnin, were accepted as visualised in TRACER, v1.6 (http://tree.bio. ed.ac.uk/software/tracer/). Maximum Clades Credibility (MCC) trees were constructed in each case to summarise the posterior tree distributions. In these MCC trees, the more epidemiologically relevant clusters and lineages, previously identified in the global ML tree, were studied; and a node support cut- off (posterior probability (pp) above 0.9) was applied for their confirmation. Trees were viewed and edited in FigTree, v. 1.4.0 (http://tree.bio.ed.ac.uk/software/figtree). Analysis of the antiretroviral drug resistance mutations Drug resistance mutations were identified in the pol sequences using the HIVseq program, which is available in the HIV Drug Resistance Database of Stanford University (https://hivdb. stanford.edu/hivseq/by-sequences/), and also using the WHO surveillance drug resistance mutation list (last updated in 2009 by Bennett and colleagues) [37]. Statistical analyses A multivariate logistic regression analysis was performed to determine the predictive effect of the demographic, clinical and virological characteristics on the adscription to each subtype/ CRF. The statistical significance of these characteristics, compared to the total proportion of infected patients, was studied by a hypothesis contrast using a z-test. The statistical analysis was performed with SPSS 22.0. Results Epidemiological surveillance of the non-B HIV-1 genetics forms Of the 693 total included patients, 165 (23.8%) were infected with different genetic forms of HIV-1 non-B variants. Most of them (n = 104, 63%) were recombinant viruses in pol: 95 (57.6%) corresponded to 12 different CRFs and nine (5.5%) were URFs. The other patients (n = 61, 37%) were infected with five non-recombinant subtypes: A1, D, C, F1 and G (see Fig 1). The demographic, clinical and virological characteristics of the patients according to the genetic HIV-1 non-B forms are provided in Table 1. Most of the patients infected with non-B variants were men (63%, p < 0.001) aged over 35 (73.5%, p < 0.001), heterosexual (92.2%, p < 0.001), African (58.2%, p < 0.001), and living in the El Ejido area (62.4%, p<0.001). The PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 4 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 1. Distribution of the HIV-1 non-subtype B genetic forms detected in Eastern Andalusia over the 2005±2012 period. https://doi.org/10.1371/journal.pone.0186928.g001 full list of countries of origin for patients infected with non-B forms and born abroad (n = 127 [77%]) were: Argentina, n = 1, Brazil, n = 5, Burkina Faso, n = 1, Cameroon, n = 1, Colombia, n = 1, Congo, n = 5, Ivory Coast, n = 1, Cuba, n = 1, Gambia, n = 2, Ghana, n = 13, Guinea, n = 10, Guinea-Bissau, n = 11, Equatorial Guinea, n = 4, Lithuania, n = 1, Mali, n = 9, Morocco, n = 3, Mauritania, n = 1, Nigeria, n = 17, Dominican Republic, n = 1, Romania, n = 6, Russia, n = 14, Senegal, n = 15, Sierra Leone, n = 1 and South Africa, n = 3. The rest of subjects (n = 36, 22%) had been born in Spain. The multivariate logistic regression analyses demonstrated a higher risk of carrying HIV-1 subtype A for females (OR = 6.17, p = 0.026) and non Africans (OR = 0.08, p = 0.008; S3 Table). The other HIV-1 non-B genetic forms showed no predictive effect of the demographic, clinical and virological characteristics (data not shown). Twenty-three patients were infected with unusual HIV-1 non-B variants (i.e., those variants found in four patients or fewer). Of them, 10 (44.4%) were observed in Spanish patients (Table 2). The recombination patterns for the different URFs obtained according to the Boot- scan analysis are presented in Fig 2. Geographical distribution of the various HIV-1 non B genetic forms The geographic distribution of the different HIV-1 non B subtypes and recombinant forms are represented on the map of Eastern Andalusia (Fig 3). Most of the patients infected with HIV-1 non-B variants were sampled in El Ejido (62.4%) or in the city of Granada, (22.4%), whereas   non-B variants were less frequent in the cities of Almerõa (10.9%), Jaen (2.4%) and Motril (1.8%). PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 5 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Analysis of the putative geographical origins of the main HIV-1 non-B genetic forms found in Eastern Andalusia In order to characterize the phylogenetic relationship of the patients infected with the most frequently found HIV-1 non-B variants (those found in 5 patients), the global ML tree (Fig 4) revealed the existence of 13 international lineages in Eastern Andalusia (Table 3) and 11 Andalusian clusters (Table 4) that involved patients in our cohort. The Bayesian analyses (Figs 5 and 6) showed that most of these Andalusian clusters origi- nated in the first decade of this century, and mainly included patients sampled in El Ejido. The low CD4 count of the patients included in most of these transmission networks suggests a late HIV diagnosis in a high proportion of patients (Table 4). Table 1. Demographic, clinical and virological characteristics of the patients infected with HIV-1 non B variants sampled over the 2005±2012 period. NON-RECOMBINANT SUBTYPES RECOMBINANT FORMS CHARACTERISTICS A1 C F1 G CRF02_AG CRF14_BG CRF06_cpx Others Total Gender Male 10 (33.3) 4 (57.1) 8 (88.9) 9 (75) 48 (69) 6 (75) 3 (60) 15 (65.2) 103 (62.4) Female 20 (66.7) 3 (42.9) 1 (11.1) 3 (25) 23 (31) 2 (25) 2 (40) 8 (34.8) 62 (37.6) Age (159)* <35 12 (40) 3 (42.9) 4 (50) 3 (25) 16 (23.2) 1 (20) 3 (15) 42 (26.4) 35±45 12 (40) 3 (42.9) 2 (25) 4 (33.3) 31 (44.9) 5 (62.5) 4 (80) 13 (65) 74 (46.5) >45 6 (20) 1 (14.3) 2 (25) 5 (41.7) 22 (31.9) 3 (37.5) 4 (20) 43 (27) Risk factor (153)* HTX 30 (100) 7 (100) 8 (88.9) 12 (100) 61 (95.3) 6 (75) 3 (75) 14 (73.7) 141 (92.2) MSM 1 (1.6) 1 (12.5) 1 (25) 5 (26.3) 8 (5.2) IVDU 1 (11.1) 2 (3.1) 1 (12.5) 4 (2.6) Country of origin (163)* Spanish 11 (36.7) 4 (44.4) 1 (8.3) 6 (8.6) 4 (50) 10 (47.6) 36 (22) East Europe 15 (50) 2 (28.6) 2 (22.2) 2 (9.5) 21 (12.9) North Africa 3 (4.2) 3 (1.8) West Africa 2 (6.7) 2 (28.6) 1 (11.1) 11 (91.7) 54 (76) 4 (50) 5 (100) 3 (14.3) 82 (50.3) Central and South Africa 1 (3.3) 1(14.3) 8 (11.3) 2 (9.5) 12 (7.4) Central and South America 1 (3.3) 2 (28.6) 2 (22.2) 4 (19) 9 (5.5) Viral Load (158)* <10000 8 (26.7) 1 (14.3) 1 (12.5) 1 (8.3) 15 (22.1) 1 (12.5) 5 (23.8) 32 (20.3) 10000±100000 10 (33.3) 2 (28.6) 1 (12.5) 6 (50) 24 (35.3) 5 (62.5) 2 (50) 9 (42.9) 59 (37.3) >100000 12 (40) 4 (57.1) 6 (75) 5 (41.7) 29 (42.6) 2 (25) 2 (50) 7 (33.3) 67 (42.4) CD4 count (149)* <200 4 (16) 4 (50) 6 (50) 35 (52.2) 2 (25) 1 (20) 1 (5.3) 53 (35.6) 201±350 6 (24) 3 (60) 3 (37.5) 3 (25) 19 (28.4) 3 (37.5) 2 (40) 4 (21.1) 43 (28.9) >350 15 (60) 2 (40) 1 (12.5) 3 (25) 13 (19.4) 3 (37.5) 2 (40) 14 (73.7) 53 (35.6) Sampling date interval 2005±2007 4 (13.3) 3 (42.9) 6 (50) 11 (15.5) 1 (12.5) 3 (13) 28 (17) 2007±2009 3 (10) 1 (14.3) 1 (11.1) 2 (16.7) 30 (42.3) 5 (62.5) 4 (80) 3 (13) 49 (29.7) 2009±2011 10 (33.3) 2 (28.6) 1 (11.1) 3 (25) 12 (16.9) 1 (12.5) 1 (20) 8 (34.8) 38 (23) 2011±2013 13 (43.3) 1 (14.3) 7 (77.8) 1 (8.3) 18 (25.4) 1 (12.5) 9 (39.1) 50 (30.3) Total 30 7 9 12 71 8 5 23 165 *Date available for the number of indicated patients. https://doi.org/10.1371/journal.pone.0186928.t001 PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 6 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Table 2. Clinical, demographic and virological characteristics of the patients infected with infrequent HIV-1 non B genetic variants over the 2005± 2012 period. Patient ID Genetic form Gender (M/F) Risk Factor Age Year of diagnosis CD4 count Viral Load Country of origin 1 CRF01_AE M HTX 38 2008 764 421 Spain 2 CRF03_AB M MSM 28 2011 715 404000 Spain 3 CRF20_BG M HTX 65 2005 521 5910 Spain 4 CRF20_BG F HTX 40 2005 476 77700 Cuba 5 CRF22_01_A1 M HTX 37 2009 474 11700 Equatorial Guinea 6 CRF25_CPX F UNK NA 2012 NA 36300 Spain 7 CRF26_AU F HTX 25 2012 15 1180000 Spain 8 CRF37_CPX M HTX 39 2007 405 135000 Equatorial Guinea 9 CRF47_BF M HTX 39 2009 389 12279 Spain 10 CRF47_BF F HTX 32 2011 392 4122 Spain 11 CRF49_CPX M HTX 36 2009 288 2250000 Senegal 12 D F HTX 39 2011 443 22233 Equatorial Guinea 13 D M MSM 48 2009 526 55723 Spain 14 D M MSM 41 2009 416 417477 Colombia 15 URF BC F HTX 35 2011 723 27911 Romania 16 URF BC M HTX 40 2011 450 211933 Romania 17 URF BF M MSM 38 2011 246 143039 Spain 18 URF BF F HTX 38 2012 NA 44100 Brazil 19 URF BF M HTX 52 2006 462 3200 Argentina 20 URF BG F NA NA 2008 NA NA NA 21 URF BG M NA NA 2009 NA NA NA 22 URF BG M HTX 38 2009 NA NA Guinea-Bissau 23 URF BD M NA 49 2009 259 2800 Spain NA: Not Available; M: Male; F: Female. https://doi.org/10.1371/journal.pone.0186928.t002 In order to provide more information about the scale of the trees shown, we provide in the S4 Table the distribution of patristic (uncorrected) pairwise genetic distances between sequences included in each of the ML and Bayesian trees generated in this article. Non-recombinant subtypes Thirty (18.2%) patients were infected with HIV-1 subtype A1. The viral sequences were geneti- cally similar according to HIV-BLAST to 21 GenBank sequences from Bulgaria, the Demo- cratic Republic of Congo, Croatia and Greece with 13, 6, 1 and 1 cases, respectively. The ML analysis (Fig 4) detected a large international lineage that involved sequences from Eastern Europe (lineage L1.A1 in Table 3) and grouped 21 patients from our cohort: 16 women born abroad (Eastern Europe (n = 14), the Dominican Republic (n = 1) and Lithuania (n = 1)) and 5 Spanish men. This lineage also included 23 GenBank sequences, also originating from Eastern Europe: Bulgaria, n = 10, Russia, n = 5, Poland, n = 1 and the Ukraine, n = 1. Within this line- age, we found two clusters (A.1 and A.2), formed exclusively by Spanish men and female sex workers born in Russia, all being patients sampled in Eastern Andalusia. The A.1 local cluster involved 4 sequences from Spanish patients living in the capital of Granada, its origin was esti- mated to be 2008.5 (95%CI: 2006.6±2010.3), and the sequences presented the resistance muta- tion K103N in the RT gene. This cluster was also phylogenetically related to viruses that circulate in Eastern Europe. Unlike most of the HIV-1 non-B clusters, patients in the A.1 clus- ter showed a high CD4 count (mean = 590, range = 534±701). Moreover, the Bayesian PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 7 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 2. Bootscan analysis of the unique recombinant forms (URF) found in Eastern Andalusia. The analysis was applied to the concatenated sequences that corresponded to HXB2 coordinates 2283±2549 (PR) and 2661±3290 (RT). https://doi.org/10.1371/journal.pone.0186928.g002 phylogenetic tree revealed short internode branches, which may indicate short times between infections. Finally, five subtype A1 sequences from our cohort corresponded to patients from Africa: Mali, n = 2, Equatorial Guinea, n = 1 and Spain, n = 2, not clustered in transmission cluster. Seven (4.2%) sequences corresponded to HIV-1 subtype C, and showed high genetic simi- larity to 23 GenBank sequences sampled in South Africa (n = 12), Brazil (n = 6) and Bulgaria (n = 4). We thus found two main ways of subtype C entrance to our area: South Africa and Brazil: a Brazilian male patient from our cohort grouped with 6 GenBank sequences from Bra- zil; and a South African male patient grouped with GenBank sequences from South Africa (n = 2) and Somalia (n = 1) (Fig 4). Within this subtype, we also found a single cluster (C.1) formed by patients from Brazil (n = 1) and Romania (n = 2). Nine (5.5%) sequences corresponded to HIV-1 subtype F1 and showed a high genetic simi- larity to 21 GenBank sequences from Brazil (n = 14), Bulgaria (n = 4) and the Democratic Republic of Congo (n = 3). We found only one Andalusian F1 cluster: a sequence pair (cluster F.1), that originated in 2010.2 (95%CI: 2010±2011) and was formed by two male injection drug users sampled in Jaen and who were of Brazilian and Spanish origins. This sequence pair was included among GenBank sequences from Brazil in the ML tree. However, we found 2 international lineages: L1.F, which grouped two Romanian heterosexual patients from our cohort with GenBank sequences sampled in Eastern Europe, mainly Romania (n = 6) and Bul- garia (n = 2). The second F1 subtype lineage (L3.F) included Spanish men who have sex with men (MSM) sampled in North Spain, and also a MSM from our cohort. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 8 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 3. Geographical distribution of the patients infected with HIV-1 non-B variants over the 2005±2012 period. The percentage of each subtype/CRF in relation to all the HIV-1 non-B genetic forms is shown in each region. https://doi.org/10.1371/journal.pone.0186928.g003 Twelve (7.3%) patients of our cohort were infected with HIV-1 subtype G, who came from different western and central African countries: Mali (n = 1), Nigeria (n = 6), Ghana (n = 3) and Guinea-Bissau (n = 1). They presented high genetic similarity to 5 GenBank sequences from the Republic of Congo (n = 4) and Bulgaria (n = 1). None of these sequences was epide- miologically related according to our data. We found only one Nigerian patient whose sequence grouped with another one of the same country of origin (lineage L1.G). HIV-1 recombinant forms Eight (4.8%) patients in our cohort were infected with the recombinant CRF14_BG form, which in the pol analyses typically forms a monophyletic cluster within the subtype G crown. These eight patients came from Spain (n = 4), Guinea (n = 2), and Guinea-Bissau (n = 2). We also found a sin- gle small Andalusian cluster (cluster 14BG.1), which originated in 2004.4 (95% CI:2003.8±2005), and was formed by two Spanish patients. Finally, two patients from Guinea and Guinea Bissau grouped with sequences from Equatorial Guinea (lineages L1.14BG and L2.14BG). We found 71 (43%) patients, mainly from western African countries (77.5%), infected with CRF02_AG. Of these, 11 (14%) were grouped into five small Andalusian Cluster: 4 clusters of two patients and one with three patients. We detected 5 different lineages (L1.02AG-L5.02AG) of viruses sampled in other countries, with patients from our cohort who came mainly from Western Africa. To study the phylogenetic profile of variant CRF06_cpx, we used all the sequences available in Los Alamos HIV given their small number, n = 110 (see Fig 4). We found 5 patients in our cohort (3%) to be infected with variant CRF06_cpx, who came from different western African countries: Nigeria (n = 3), Ghana (n = 1) and Senegal (n = 1). These sequences grouped with GenBank sequences from the neighboring Western African countries of Burkina Faso, Togo and Nigeria. However, we found no significant association among the patients infected with PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 9 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 4. Global ML phylogenetic tree inferred for the main HIV-1 non-B genetic forms sampled in Eastern Andalusia. The phylogenetic tree was constructed by the general time-reversible with gamma-distributed rate heterogeneity across sites model of substitution implemented into RAxML. Branches are drawn on scale with the bar at the bottom, which represents 0.04 nucleotide substitution per site. Statistically highly supported nodes (bootstrap values>70%) are indicated by an asterisk (*). Andalusian clusters and international lineages are highlighted in yellow and blue, respectively. The Andalusian sequence names contain a three-part code: Sequence number, sampling site (AL: Almerõ Âa, EJ: El Ejido GR: Granada, JA: Jae  n, MO: Motril) and the code of the most likely country of infection. https://doi.org/10.1371/journal.pone.0186928.g004 this genetic form, and the CRF06_cpx sequences sampled in our cohort were interspersed in the tree. Discussion In Eastern Andalusia, most HIV-1 non-B subtype genetic forms were found among immigrant heterosexual population, mainly African males or Eastern European females. These patients were living preferentially in El Ejido, an area that potentially acts as a gateway for diverse HIV- PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 10 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Table 3. HIV-1 non-B international lineages involving sequences sampled in Eastern Andalusia and sequences from different countries. Lineages Sequences Country of origin of the Node Country of Most likely countries of infection for Articles with (EA;OR) patients sampled in our support collection for other other patients (Accession ID) related cohort (n) (BT) sequences (n) sequences Subtype A L1.A 21;23 Spanish (5), Russian 97 Spain (11) Ukraine (KC340416) [20,38,39] (14), Lithuanian (1) Bulgaria (10) Russia (KC340417, KC340407, Dominican (1) Poland (1) KC340486, KC340436, KC340495) Croatia (1) Subtype C L1.C 1;1 Senegalese (1) 100 Spain (1) Morocco (JQ351995) [40] Subtype F L1.F 2:21 Romanian (2) 70 Spain (19) Romania (KC340133, KC340134, [20,39,41] Bulgaria (2) KC340419, KC340378, KC340643, KC340132) L2.F 1;1 Spanish (1) 100 Spain (1) Sub-Saharan Africa (KC340440) L3.F 1;112 Spanish (1) 100 Spain (112) Spain (KJ883085-KJ883089, [42,43] KJ883091-KJ883108, KJ883110-KJ883152, JN010216-JN859590) Subtype G L1.G 1;1 Nigerian (1) 100 Spain (1) Nigeria (FJ481667) [44] CRF_14BG L1.14BG 1;1 Guinean-Bissau (1) 100 Spain (1) Equatorial-Guinea (JX428555) [45] L2.14BG 1;1 Guinean (1) 100 Spain (1) Equatorial-Guinea (EU255306) [45] Form CRF02_AG L1.02AG 1;2 Guinean (1) 100 Spain (2) Ghana (KC340379) [23] L2.02AG 2;1 Moroccan (1) 86 Spain (1) Cameroon (KC340644) Spanish (1) L3.02AG 1;1 Nigerian (1) 86 Spain (1) Spain (HF567877) [46] L4.02AG 1;6 Ghanaian (1) 96 Spain (6) Bolivia (EU255444) [44,47] Ecuador (KC340123, KC340124, FJ481711) L5.02AG 1;5 Equatorial-Guinean (1) 82 Spain (5) Equatorial-Guinea (EU255527) [47,48] EA: sequences sampled in Eastern Andalusia; OR: sequences sampled in other geographic regions; BT: bootstrap (in the global maximum-likelihood tree). https://doi.org/10.1371/journal.pone.0186928.t003 1 variants to enter the Eastern Andalusian region. These findings are explained by the fact the El Ejido's economy is mainly based on greenhouse farming, for which a large industry has emerged in recent years thanks to immigrant labor, made up of people mainly from Africa. The prevalence of HIV-1 non-B variants in eastern Andalusia is similar to that reported in a study performed in the nearby Western areas of Andalusia (23%) [49], but is still much higher than that found elsewhere in Spain [6,7]. An increased prevalence has been noted for HIV-1 non-B variants and their genetic diversity in Eastern Andalusia in recent years: 22% of autochthonous patients were infected with HIV-1 non-B forms between 2005 and 2012, as opposed to the 12.8% reported in former studies conducted between 1997 and 2001 [50]. We also detected 12 different CRFs and nine URFs, a variability that is probably related to the increased migration rate reported in southern Spain in the last decade [49,51]. The least frequent HIV-1 non-B variants were detected often among Spanish patients (43%, [10/23]), and most of the clusters formed by these variants included at least one Spanish patient (55%, [6/11]). These data suggest that although these HIV-1 non-B variants seem to be PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 11 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Table 4. Demographical, clinical, virological and phylogenetic characteristics of the Andalusian clusters found for the main HIV-1 non B variants. Cluster No. Sampling Support Risk Country of Viral Load CD4 count Location area tMRCA (95% ID Patients interval (BT:PP) factor origin (median, IQR) (median, HPD) Log IQR) Subtype A A.1 4 2011 100;0.9 HTX Spanish 5 (4.6±6) 590 (534± Granada 2008.5 (2006.6± 701) 2010.3) A.2 3 2007±2011 96;0.9 HTX Spanish- 4.7 (4.4±5.3) 301 (152± Granada-El 2003.9 (1998.4± Russian 517) Ejido 2004.7) A.3 2 2009±2010 100;0.9 HTX Spanish- 4.9 (4±5) 391 (325± El Ejido 2008 (2007.4± Russian 456) 2009) Subtype C C.1 3 2005±2012 100;0.9 HTX Romanian- 5.2 (5.1±5.3) 499 (253± Granada-Ejido 1998.3 (1992.3± Brazilian 739) 2001.9) Subtype F F. 2 2011±2012 90;0.9 HTX Brazilian 5.7(5.4±5.9) 62 (50±73) Jaen 2010.2 (2010± 2011) Form CRF14_BG 14BG.1 2 2005±2007 100;0.9 HTX Spanish 5 425 (267± Granada 2004.4 (2003.8± 583) 2005) Form CRF02AG 02AG.1 2 2007 100;0.9 HTX Guinean 5.6(4±6) 309 (180± El Ejido 2003.8 (2001± 438) 2006) 02AG.2 2 2007±2010 100;0.9 HTX Senegalese 3(4±5) 391 (381± El Ejido 2006.4 (2006± 400) 2007) 02AG.3 2 2011 99;0.9 HTX Spanish- 5.9(3±6) 326 (100± Granada- 2008.2 (2007.2± Moroccan 551) Almerõ Âa 2009.9) 02AG.4 3 2006±2011 98;0.9 HTX Spanish- Mali- 5.5(5.3±6.2) 140 (78± El Ejido 2003.6 (2000.3± 306) 2005.3) 02AG.5 2 2008 100;0.9 HTX Ghana 6 126 (45± El Ejido 2007.7 (2007.2± 200) 2009) https://doi.org/10.1371/journal.pone.0186928.t004 due to imported cases in most cases, they have also gradually penetrated the autochthonous population in recent years. The phylogenetic and epidemiological study of the HIV-1 non-B variants in our region showed that these variants account for high proportion of infections among migrant patients, and that these viruses were genetically close to those circulating in these subjects' countries of origin. This indicates that many patients were infected before they arrived in Spain. These sequences sampled in other countries, and available in public databases, act as a control to avoid overestimating the local transmission clusters that include patients who are most proba- bly unrelated in epidemiological terms. As previously shown in a national study [7], CRF02_AG was the most frequent HIV non-B variant in our population (43%). Nonetheless, the small proportion of their phylogenetic asso- ciation is surprising (14%, [11/71]). This clustering rate was much higher for other HIV-1 non-B subtypes, such as subtype A1 (27.3%, [9/33]), where we discovered an international lineage (L1.A1) that mostly included a particularly vulnerable group of Russian female sex workers and potentially their local customers. According to our analysis, it would appear that most of the non-B cases detected in Eastern Andalusia were generally imported cases as most were identified in immigrant populations. Our analyses suggest that many of these cases form part of international HIV-1 lineages that originated in Eastern Europe, South America and sub-Saharan Africa. However, we also PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 12 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 5. Bayesian phylogenetic tree inferred for the subtype A1, C, F1 and G/CRF14_BG pol sequences sampled in Eastern Andalusia and genetically similar sequences from GenBank. Red branches correspond to the sequences sampled in eastern Andalusia from 2005 to 2012. Statistically highly supported nodes (posterior probability values above 0.9) are indicated with an asterisk (*). Andalusian clusters are highlighted in yellow. Andalusian sequences names contain a three-part code: Sequence number, sampling site (AL: Almerõ Âa, EJ: El Ejido GR: Granada, JA: Jae  n, MO: Motril) and the code of the most likely country of infection. https://doi.org/10.1371/journal.pone.0186928.g005 PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 13 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 6. Bayesian phylogenetic tree inferred for the CRF02_AG pol sequences sampled in Eastern Andalusia and genetically similar sequences from GenBank. Red branches correspond to the sequences sampled in eastern Andalusia from 2005 to 2012. Statistically highly supported nodes (posterior probability values above 0.9) are indicated with an asterisk (*). Andalusian clusters are highlighted in yellow. Andalusian sequences names contain a three-part code: Sequence number, sampling site (AL: Almerõ Âa, EJ: El Ejido GR: Granada, JA: Jae  n, MO: Motril) and the code of the most likely country of infection. https://doi.org/10.1371/journal.pone.0186928.g006 PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 14 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain identified 11 intra-region clusters, which might suggest the local dissemination of some non-B variants, particularly those which involve autochthonous Spanish subjects (6/11) and recent emergence times according to the phylogenetic reconstruction. On the other hand, clusters formed by foreign subjects with old common ancestors most likely reflect imported infections. The methods used herein involve a number of sampling limitations that affect this and many other similar studies. Since we relied on a BLAST search to identify the genetically clos- est sequences (from both Spain and abroad) that could form part of the same transmission net- works as our sequences, we depended on the sequences deposited in databases. Unfortunately, this availability is sometimes very low, particularly for non-B variants. Therefore, we cannot rule out that close and more informative sequences were not captured as they have not been sampled. This was the reason why we added all the sequences available in HIV Los Alamos col- lected in Spain. We also demonstrated the presence of 13 different lineages of viruses that cir- culated in our region, which grouped with other patients from different Spanish cohorts, mainly foreign patients. Fortunately, very few sequences included in transmission clusters persented resistance to first-line antiretroviral drugs. This information agrees with the common conception that viruses with resistance mutations present a biological disadvantage against wild strains, which weakens their transmission efficacy. Likewise, the drug resistance mutations detected affect mainly reverse transcriptase inhibitor drugs. We detected transmitted resistance mutations in four of the five patients grouped in Cluster A.4, which would cause high level resistance to nevirapine and efavirenz. We decided to study only the resistance mutations present in trans- mission clusters, which would have a stronger epidemiological impact. Further detailed infor- mation will be provided in future works. The constant epidemiological surveillance in our population, for which phylogenetic analysis tools are used, is a particularly important measure to study past outbreaks of genetic HIV-1 non- B subtype variants, and to prevent future ones. Likewise, as transmission cluster size seems to predict its expansion in time [52], we could expect some transmission chains of HIV-1 non-sub- type variants to become larger in size in forthcoming years, and more Spanish individuals to be included. We herein detected the presence of one patient from our cohort related to a fast spread- ing cluster among Spanish MSM infected with subtype F in Galicia (NW Spain) [43], a transmis- sion cluster which, as Delgado et al. suggest, would probably be closely linked to viruses that circulate in Eastern Europe [42]. These authors [53] have also described a subtype A cluster that is being transmitted among individuals in different areas of Spain. Finally, Patiño et al. [54] have warned about the novel appearance of variant CRF19_cpx among Spanish MSM individuals. Adequate knowledge about the characteristics of local epidemics, the study of risk groups and the prevalence of different viral subtypes are all fundamental aspects to successfully design HIV-1 prevention campaigns. In the present study, we demonstrate that phylogenetic studies which combine demographic, clinical and geographical data from different HIV-1 non-B subtypes in Eastern Andalusia provide very useful information to epidemiologically monitor and control HIV-1 spread and its origin in imported cases. Its use will help to reinforce and implement effi- cient actions to prevent HIV-1 from spreading between autochthonous and migrant populations. Supporting information S1 Table. HIV-1 reference sequence dataset used in the phylogenetic analysis. (DOCX) S2 Table. Evolutionary rates for each of the main HIV-1 lineages found in this study obtained through Bayesian phylogenetic inference. (DOCX) PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 15 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain S3 Table. Multivariate logistic regression analysis performed for the HIV-1 subtype A1 infections. (DOCX) S4 Table. Distribution of patristic (uncorrected) genetic distances among the HIV-1 pol sequences in each dataset included in this study. (DOCX) Author Contributions Conceptualization: Santiago Pe Ârez-Parra, Natalia Chueca, Marta Alvarez, Gonzalo Yebra. Data curation: Santiago Pe Ârez-Parra, Juan Pasquau, Mohamed Omar, Antonio Collado, David Vinuesa, Ana Belen Lozano, Federico Garcõ Âa. Formal analysis: Santiago Pe Ârez-Parra, Federico Garcõ Âa. Funding acquisition: Federico Garcõ Âa. Investigation: Santiago Pe Ârez-Parra, Natalia Chueca, Marta Alvarez, Federico Garcõ Âa. Methodology: Santiago Pe Ârez-Parra, Natalia Chueca, Gonzalo Yebra, Federico Garcõ Âa. Project administration: Federico Garcõ Âa. Resources: Federico Garcõ Âa. Software: Santiago Pe Ârez-Parra. Supervision: Santiago Pe Ârez-Parra, Juan Pasquau, Mohamed Omar, Antonio Collado, David Vinuesa, Ana Belen Lozano, Gonzalo Yebra, Federico Garcõ Âa. Validation: Santiago Pe Ârez-Parra, Federico Garcõ Âa. Visualization: Santiago Pe Ârez-Parra, Federico Garcõ Âa. Writing ± original draft: Santiago Pe Ârez-Parra, Gonzalo Yebra, Federico Garcõ Âa. Writing ± review & editing: Santiago Pe Ârez-Parra, Gonzalo Yebra, Federico Garcõ Âa. References 1. Vidal N, Mulanga C, Bazepeo SE, Mwamba JK, Tshimpaka J-W, Kashi M, et al. Distribution of HIV-1 variants in the Democratic Republic of Congo suggests increase of subtype C in Kinshasa between 1997 and 2002. J Acquir Immune Defic Syndr 1999. 2005 Dec 1; 40(4):456±62. 2. Taniguchi Y, Takehisa J, Bikandou B, Mboudjeka I, N'Doundou-N'Kodia M-Y, Obengui null, et al. Genetic subtypes of HIV type 1 based on the vpu/env sequences in the Republic of Congo. AIDS Res Hum Retroviruses. 2002 Jan 1; 18(1):79±83. https://doi.org/10.1089/088922202753394745 PMID: 3. Shankarappa R, Margolick JB, Gange SJ, Rodrigo AG, Upchurch D, Farzadegan H, et al. Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol. 1999 Dec; 73(12):10489±502. PMID: 10559367 4. Robertson DL, Anderson JP, Bradac JA, Carr JK, Foley B, Funkhouser RK, et al. HIV-1 nomenclature proposal. Science. 2000 Apr 7; 288(5463):55±6. PMID: 10766634 5. Hemelaar J, Gouws E, Ghys PD, Osmanov S, WHO-UNAIDS Network for HIV Isolation and Characteri- sation. Global trends in molecular epidemiology of HIV-1 during 2000±2007. AIDS Lond Engl. 2011 Mar 13; 25(5):679±89. 6. Yebra G, de Mulder M, Martõ Ân L, Rodrõ Âguez C, Labarga P, Viciana I, et al. Most HIV type 1 non-B infec- tions in the Spanish cohort of antiretroviral treatment-naïve HIV-infected patients (CoRIS) are due to recombinant viruses. J Clin Microbiol. 2012 Feb; 50(2):407±13. https://doi.org/10.1128/JCM.05798-11 PMID: 22162552 PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 16 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain 7. Garcõ Âa F, Pe  rez-Cachafeiro S, Guillot V, Alvarez M, Pe  rez-Romero P, Pe  rez-Elõ Âas MJ, et al. Transmis- sion of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treat- ment naïve HIV-infected individuals (CoRIS). Antiviral Res. 2011 Aug; 91(2):150±3. https://doi.org/10. 1016/j.antiviral.2011.05.010 PMID: 21663768 8. Artenstein AW, VanCott TC, Mascola JR, Carr JK, Hegerich PA, Gaywee J, et al. Dual infection with human immunodeficiency virus type 1 of distinct envelope subtypes in humans. J Infect Dis. 1995 Apr; 171(4):805±10. PMID: 7706806 9. van Harmelen J, Wood R, Lambrick M, Rybicki EP, Williamson AL, Williamson C. An association between HIV-1 subtypes and mode of transmission in Cape Town, South Africa. AIDS Lond Engl. 1997 Jan; 11(1):81±7. 10. Tscherning C, Alaeus A, Fredriksson R, Bjo È rndal A, Deng H, Littman DR, et al. Differences in chemo- kine coreceptor usage between genetic subtypes of HIV-1. Virology. 1998 Feb 15; 241(2):181±8. https://doi.org/10.1006/viro.1997.8980 PMID: 9499793 11. Kanki PJ, Hamel DJ, Sankale  JL, Hsieh C c, Thior I, Barin F, et al. Human immunodeficiency virus type 1 subtypes differ in disease progression. J Infect Dis. 1999 Jan; 179(1):68±73. https://doi.org/10.1086/ 314557 PMID: 9841824 12. Apetrei C, Descamps D, Collin G, Loussert-Ajaka I, Damond F, Duca M, et al. Human immunodefi- ciency virus type 1 subtype F reverse transcriptase sequence and drug susceptibility. J Virol. 1998 May; 72(5):3534±8. PMID: 9557632 13. Taylor BS, Sobieszczyk ME, McCutchan FE, Hammer SM. The challenge of HIV-1 subtype diversity. N Engl J Med. 2008 Apr 10; 358(15):1590±602. https://doi.org/10.1056/NEJMra0706737 PMID: 14. Alaeus A, Lidman K, So È nnerborg A, Albert J. Subtype-specific problems with quantification of plasma HIV-1 RNA. AIDS Lond Engl. 1997 Jun; 11(7):859±65. 15. Parekh B, Phillips S, Granade TC, Baggs J, Hu DJ, Respess R. Impact of HIV type 1 subtype variation on viral RNA quantitation. AIDS Res Hum Retroviruses. 1999 Jan 20; 15(2):133±42. https://doi.org/10. 1089/088922299311556 PMID: 10029245 16. Brennan CA, Bodelle P, Coffey R, Harris B, Holzmayer V, Luk K-C, et al. HIV global surveillance: foun- dation for retroviral discovery and assay development. J Med Virol. 2006; 78 Suppl 1:S24±9. 17. Gilbert MTP, Rambaut A, Wlasiuk G, Spira TJ, Pitchenik AE, Worobey M. The emergence of HIV/AIDS in the Americas and beyond. Proc Natl Acad Sci U S A. 2007 Nov 20; 104(47):18566±70. https://doi. org/10.1073/pnas.0705329104 PMID: 17978186 18. Gray RR, Tatem AJ, Lamers S, Hou W, Laeyendecker O, Serwadda D, et al. Spatial phylodynamics of HIV-1 epidemic emergence in east Africa. AIDS Lond Engl. 2009 Sep 10; 23(14):F9±17. 19. Esbjo È rnsson J, Mild M, Månsson F, Norrgren H, Medstrand P. HIV-1 molecular epidemiology in Guinea- Bissau, West Africa: origin, demography and migrations. PloS One. 2011; 6(2):e17025. https://doi.org/ 10.1371/journal.pone.0017025 PMID: 21365013 20. Gonza  lez-Alba JM, Holguõ Ân A, Garcia R, Garcõ Âa-Bujalance S, Alonso R, Sua  rez A, et al. Molecular Sur- veillance of HIV-1 in Madrid, Spain: a Phylogeographic Analysis !. J Virol. 2011 Oct; 85(20):10755±63. https://doi.org/10.1128/JVI.00454-11 PMID: 21795343 21. Aldous JL, Pond SK, Poon A, Jain S, Qin H, Kahn JS, et al. Characterizing HIV transmission networks across the United States. Clin Infect Dis Off Publ Infect Dis Soc Am. 2012 Oct; 55(8):1135±43. 22. Grabowski MK, Redd AD. Molecular tools for studying HIV transmission in sexual networks. Curr Opin HIV AIDS. 2014 Mar; 9(2):126±33. https://doi.org/10.1097/COH.0000000000000040 PMID: 24384502 23. Cuevas MT, Muñoz-Nieto M, Thomson MM, Delgado E, Iribarren JA, Cilla G, et al. HIV-1 Transmission Cluster With T215D Revertant Mutation Among Newly Diagnosed Patients From the Basque Country, Spain: JAIDS J Acquir Immune Defic Syndr. 2009 May; 51(1):99±103. https://doi.org/10.1097/QAI. 0b013e318199063e PMID: 19282784 24. Vega Y, Delgado E, Ferna  ndez-Garcõ Âa A, Cuevas MT, Thomson MM, Montero V, et al. Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004±2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations. PLOS ONE. 2015 May 26; 10(5):e0125699. https://doi.org/10.1371/journal.pone.0125699 PMID: 26010948 25. Perez-Parra S, Chueca-Porcuna N, Alvarez-Estevez M, Pasquau J, Omar M, Collado A, et al. [Study of human immunodeficiency virus transmission chains in Andalusia: Analysis from baseline antiretroviral resistance sequences.]. Enferm Infecc Microbiol Clin. 2015 Jan 31; 26. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinforma Oxf Engl. 2007 Nov 1; 23(21):2947±8. 27. Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylo- genetic trees. In: Gateway Computing Environments Workshop (GCE), 2010. 2010. p. 1±8. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 17 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain 28. Hillis DM, Bull JJ. An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phyloge- netic Analysis. Syst Biol. 1993 Jun 1; 42(2):182±92. 29. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999 Jan; 73(1):152±60. PMID: 9847317 30. Leigh Brown AJ, Lycett SJ, Weinert L, Hughes GJ, Fearnhill E, Dunn DT, et al. Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J Infect Dis. 2011 Nov; 204 (9):1463±9. https://doi.org/10.1093/infdis/jir550 PMID: 21921202 È Â 31. Kouyos RD, von Wyl V, Yerly S, Boni J, Taffe P, Shah C, et al. Molecular epidemiology reveals long- term changes in HIV type 1 subtype B transmission in Switzerland. J Infect Dis. 2010 May 15; 201 (10):1488±97. https://doi.org/10.1086/651951 PMID: 20384495 32. Esbjo È rnsson J, Mild M, Audelin A, Fonager J, Skar H, Bruun Jørgensen L, et al. HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries. Virus Evol. 2016 Jan; 2(1):vew010. https://doi.org/10.1093/ve/vew010 PMID: 33. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012 Aug; 29(8):1969±73. https://doi.org/10.1093/molbev/mss075 PMID: 22367748 34. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed Phylogenetics and Dating with Confidence. PLoS Biol. 2006 Mar 14; 4(5):e88. https://doi.org/10.1371/journal.pbio.0040088 PMID: 16683862 35. Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005 May; 22(5):1185±92. https://doi.org/10.1093/ molbev/msi103 PMID: 15703244 36. Hue S, Brown AE, Ragonnet-Cronin M, Lycett SJ, Dunn DT, Fearnhill E, et al. Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions. AIDS Lond Engl. 2014 Aug 24; 28(13):1967±75. 37. Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, Kiuchi M, et al. Drug Resistance Muta- tions for Surveillance of Transmitted HIV-1 Drug-Resistance: 2009 Update. PLoS ONE [Internet]. 2009 Mar 6 [cited 2016 Sep 3]; 4(3). Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2648874/ 38. Alexiev I, Shankar A, Wensing AMJ, Beshkov D, Elenkov I, Stoycheva M, et al. Low HIV-1 transmitted drug resistance in Bulgaria against a background of high clade diversity. J Antimicrob Chemother. 2015; 70(6):1874±80. https://doi.org/10.1093/jac/dkv011 PMID: 25652746 39. De Mendoza C, Garrido C, Poveda E, Corral A, Zahonero N, Treviño A, et al. Changes in drug resis- tance patterns following the introduction of HIV type 1 non-B subtypes in Spain. AIDS Res Hum Retrovi- ruses. 2009 Oct; 25(10):967±72. https://doi.org/10.1089/aid.2008.0166 PMID: 19842792 40. Trends in Drug Resistance Prevalence in HIV-1±infected Children in Madrid (PDF Download Available) [Internet]. ResearchGate. [cited 2017 Mar 15]. Available from: https://www.researchgate.net/ publication/229074263_Trends_in_Drug_Resistance_Prevalence_in_HIV-1-infected_Children_in_ Madrid 41. Ferna  ndez-Garc õa A, Cuevas MT, Vinogradova A, Rakhmanova A, Pe  rez-Alvarez L, de Castro RO, et al. Near full-length genome characterization of a newly identified HIV type 1 subtype F variant circu- lating in St. Petersburg, Russia. AIDS Res Hum Retroviruses. 2009 Nov; 25(11):1187±91. https://doi. org/10.1089/aid.2009.0140 PMID: 19943791    42. Delgado E, Cuevas MT, Domõnguez F, Vega Y, Cabello M, Fernandez-Garcõa A, et al. 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. PloS One. 2015; 10 (11):e0143325. https://doi.org/10.1371/journal.pone.0143325 PMID: 26599410   43. Thomson MM, Fernandez-Garcõ Âa A, Delgado E, Vega Y, Dõ Âez-Fuertes F, Sanchez-Martõ Ânez M, et al. Rapid expansion of a HIV-1 subtype F cluster of recent origin among men who have sex with men in Galicia, Spain. J Acquir Immune Defic Syndr 1999. 2012 Mar 1; 59(3):e49±51. 44. Yebra G, de Mulder M, del Romero J, Rodrõ Âguez C, Holguõ Ân A. HIV-1 non-B subtypes: High transmitted NNRTI-resistance in Spain and impaired genotypic resistance interpretation due to variability. Antiviral Res. 2010 Feb; 85(2):409±17. https://doi.org/10.1016/j.antiviral.2009.11.010 PMID: 20004217 45. Yebra G, de Mulder M, Pe  rez-Elõ Âas MJ, Pe  rez-Molina JA, Gala  n JC, Llenas-Garcõ Âa J, et al. Increase of transmitted drug resistance among HIV-infected sub-Saharan Africans residing in Spain in contrast to the native population. PloS One. 2011; 6(10):e26757. https://doi.org/10.1371/journal.pone.0026757 PMID: 22046345 46. Bracho MA, Sentandreu V, Alastrue  I, Belda J, Juan A, Ferna  ndez-Garcõ Âa E, et al. Emerging trends in CRF02_AG variants transmission among men who have sex with men in Spain. J Acquir Immune Defic Syndr 1999. 2014 Mar 1; 65(3):e130±3. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 18 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain 47. Holguõ Ân A, de Mulder M, Yebra G, Lo  pez M, Soriano V. Increase of non-B subtypes and recombinants among newly diagnosed HIV-1 native Spaniards and immigrants in Spain. Curr HIV Res. 2008 Jun; 6 (4):327±34. PMID: 18691031 48. Ferna  ndez-Garc õa A, Cuevas MT, Muñoz-Nieto M, Ocampo A, Pinilla M, Garcõ Âa V, et al. Development of a panel of well-characterized human immunodeficiency virus type 1 isolates from newly diagnosed patients including acute and recent infections. AIDS Res Hum Retroviruses. 2009 Jan; 25(1):93±102. https://doi.org/10.1089/aid.2008.0174 PMID: 19113978   49. de Felipe B, Perez-Romero P, Abad-Fernandez M, Fernandez-Cuenca F, Martinez-Fernandez FJ, Trastoy M, et al. Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000±2010. Virol J. 2011; 8:416. https://doi.org/10.1186/1743-422X- 8-416 PMID: 21871090 50. Alvarez M, Garcõ Âa F, Martõ Ânez NM, Garcõ Âa F, Bernal C, Vela CM, et al. Introduction of HIV type 1 non-B subtypes into Eastern Andalusia through immigration. J Med Virol. 2003 May; 70(1):10±3. https://doi. org/10.1002/jmv.10368 PMID: 12629637 51. Instituto Nacional de Estadõ Âstica (INE) [Internet]. Available from: www.ine.es 52. Brenner BG, Roger M, Stephens D, Moisi D, Hardy I, Weinberg J, et al. Transmission Clustering Drives the Onward Spread of the HIV Epidemic Among Men Who Have Sex With Men in Quebec. J Infect Dis. 2011 Oct 1; 204(7):1115±9. https://doi.org/10.1093/infdis/jir468 PMID: 21881127 53. Delgado E, Cuevas MT, Vega Y, Montero V, Sa  nchez M, Carrera C, et al. Identificacio  n de un cluster de subtipo A que se transmite entre hombres que tienen relaciones sexuales con hombres en diversas comunidades auto  nomas de España. In Ma  laga: VI Congreso Nacional de GESIDA y 8.a Reunio  n Doc- ente de la RIS (SEIMC); 2014. p. 12. Available from: http://www.gesida-seimc.org/contenidos/ congresos/anteriores/2014/gesida2014-VIcongresocomunicaciones.pdf 54. Patiño Galindo JA, Torres-Puente M, Gimeno C, Ortega E, Navarro D, Galindo MJ, et al. Expansion of the CRF19_cpx Variant in Spain. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2015 Aug; 69:146±9. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 19 / 19 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png PLoS ONE Pubmed Central

High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain: A phylogeographic approach

Loading next page...
 
/lp/pubmed-central/high-prevalence-and-diversity-of-hiv-1-non-b-genetic-forms-due-to-ude1c0N0zN

References (132)

  • J. Harmelen, R. Wood, M. Lambrick, E. Rybicki, A. Williamson, C. Williamson (1997)

    An association between HIV‐1 subtypes and mode of transmission in Cape Town, South Africa

    AIDS, 11

  • J. González-Alba, Á. Holguín, R. Garcia, S. García-Bujalance, R. Alonso, A. Suárez, R. Delgado, L. Cardeñoso, R. González, I. García-Bermejo, F. Portero, C. Mendoza, F. González-Candelas, J. Galán (2011)

    Molecular Surveillance of HIV-1 in Madrid, Spain: a Phylogeographic Analysis

    Journal of Virology, 85

  • (YebraG, de MulderM, Pérez-ElíasMJ, Pérez-MolinaJA, GalánJC, Llenas-GarcíaJ, et al Increase of transmitted drug resistance among HIV-infected sub-Saharan Africans residing in Spain in contrast to the native population. PloS One. 2011;6(10):e26757 doi: 10.1371/journal.pone.0026757 22046345)

    YebraG, de MulderM, P�rez-El�asMJ, P�rez-MolinaJA, Gal�nJC, Llenas-Garc�aJ, et al

    YebraG, de MulderM, Pérez-ElíasMJ, Pérez-MolinaJA, GalánJC, Llenas-GarcíaJ, et al Increase of transmitted drug resistance among HIV-infected sub-Saharan Africans residing in Spain in contrast to the native population. PloS One. 2011;6(10):e26757 doi: 10.1371/journal.pone.0026757 22046345, YebraG, de MulderM, Pérez-ElíasMJ, Pérez-MolinaJA, GalánJC, Llenas-GarcíaJ, et al Increase of transmitted drug resistance among HIV-infected sub-Saharan Africans residing in Spain in contrast to the native population. PloS One. 2011;6(10):e26757 doi: 10.1371/journal.pone.0026757 22046345

  • (AlaeusA, LidmanK, SönnerborgA, AlbertJ. Subtype-specific problems with quantification of plasma HIV-1 RNA. AIDS Lond Engl. 1997 6;11(7):859–65.)

    AlaeusA, LidmanK, S�nnerborgA, AlbertJ. Subtype-specific problems with quantification of plasma HIV-1 RNA. AIDS Lond Engl. 1997

    AlaeusA, LidmanK, SönnerborgA, AlbertJ. Subtype-specific problems with quantification of plasma HIV-1 RNA. AIDS Lond Engl. 1997 6;11(7):859–65., AlaeusA, LidmanK, SönnerborgA, AlbertJ. Subtype-specific problems with quantification of plasma HIV-1 RNA. AIDS Lond Engl. 1997 6;11(7):859–65.

  • Trends in Drug Resistance Prevalence in HIV-1-infected Children in Madrid (PDF Download Available

  • MTP Gilbert (2007)

    18566

    Proc Natl Acad Sci U S A, 104

  • Á. Holguín, M. Mulder, G. Yebra, Marisa López, V. Soriano (2008)

    Increase of non-B subtypes and recombinants among newly diagnosed HIV-1 native Spaniards and immigrants in Spain.

    Current HIV research, 6 4

  • J. Esbjörnsson, Mattias Mild, F. Månsson, H. Norrgren, P. Medstrand (2011)

    HIV-1 Molecular Epidemiology in Guinea-Bissau, West Africa: Origin, Demography and Migrations

    PLoS ONE, 6

  • I. Alexiev, A. Shankar, A. Wensing, D. Beshkov, I. Elenkov, M. Stoycheva, D. Nikolova, M. Nikolova, W. Switzer (2015)

    Low HIV-1 transmitted drug resistance in Bulgaria against a background of high clade diversity.

    The Journal of antimicrobial chemotherapy, 70 6

  • N Vidal (2005)

    456

    J Acquir Immune Defic Syndr 1999, 40

  • (EsbjörnssonJ, MildM, AudelinA, FonagerJ, SkarH, Bruun JørgensenL, et al HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries. Virus Evol. 2016 1;2(1):vew010 doi: 10.1093/ve/vew010 27774303)

    Esbj�rnssonJ, MildM, AudelinA, FonagerJ, SkarH, Bruun J�rgensenL, et al

    EsbjörnssonJ, MildM, AudelinA, FonagerJ, SkarH, Bruun JørgensenL, et al HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries. Virus Evol. 2016 1;2(1):vew010 doi: 10.1093/ve/vew010 27774303, EsbjörnssonJ, MildM, AudelinA, FonagerJ, SkarH, Bruun JørgensenL, et al HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries. Virus Evol. 2016 1;2(1):vew010 doi: 10.1093/ve/vew010 27774303

  • Jeannette Aldous, Sergei Pond, A. Poon, Sonia Jain, Huifang Qin, James Kahn, M. Kitahata, B. Rodriguez, Ann Dennis, S. Boswell, R. Haubrich, Davey Smith (2012)

    Characterizing HIV transmission networks across the United States.

    Clinical infectious diseases : an official publication of the Infectious Diseases Society of America, 55 8

  • A. Drummond, M. Suchard, Dong Xie, A. Rambaut (2012)

    Bayesian Phylogenetics with BEAUti and the BEAST 1.7

    Molecular Biology and Evolution, 29

  • B. Taylor, M. Sobieszczyk, F. McCutchan, S. Hammer (2008)

    The challenge of HIV-1 subtype diversity.

    The New England journal of medicine, 358 15

  • C. Tscherning, A. Alaeus, R. Fredriksson, Å. Björndal, H. Deng, D. Littman, E. Fenyö, J. Albert (1998)

    Differences in chemokine coreceptor usage between genetic subtypes of HIV-1.

    Virology, 241 2

  • (Leigh BrownAJ, LycettSJ, WeinertL, HughesGJ, FearnhillE, DunnDT, et al Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J Infect Dis. 2011 11;204(9):1463–9. doi: 10.1093/infdis/jir550 21921202)

    Leigh BrownAJ, LycettSJ, WeinertL, HughesGJ, FearnhillE, DunnDT, et al

    Leigh BrownAJ, LycettSJ, WeinertL, HughesGJ, FearnhillE, DunnDT, et al Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J Infect Dis. 2011 11;204(9):1463–9. doi: 10.1093/infdis/jir550 21921202, Leigh BrownAJ, LycettSJ, WeinertL, HughesGJ, FearnhillE, DunnDT, et al Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J Infect Dis. 2011 11;204(9):1463–9. doi: 10.1093/infdis/jir550 21921202

  • B. Felipe, P. Pérez-Romero, M. Abad-Fernández, F. Fernández-Cuenca, Francisco Martinez-Fernandez, M. Trastoy, R. Mata, L. López-Cortés, M. Leal, P. Viciana, A. Vallejo (2011)

    Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000-2010

    Virology Journal, 8

  • Federico García, Santiago Pérez-Cachafeiro, V. Guillot, M. Álvarez, P. Pérez-Romero, M. Pérez-Elías, I. Viciana, J. Blanco, María López-Diéguez, C. Mendoza (2011)

    Transmission of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treatment naïve HIV-infected individuals (CoRIS).

    Antiviral research, 91 2

  • (van HarmelenJ, WoodR, LambrickM, RybickiEP, WilliamsonAL, WilliamsonC. An association between HIV-1 subtypes and mode of transmission in Cape Town, South Africa. AIDS Lond Engl. 1997 1;11(1):81–7.)

    van HarmelenJ, WoodR, LambrickM, RybickiEP, WilliamsonAL, WilliamsonC. An association between HIV-1 subtypes and mode of transmission in Cape Town, South Africa. AIDS Lond Engl. 1997

    van HarmelenJ, WoodR, LambrickM, RybickiEP, WilliamsonAL, WilliamsonC. An association between HIV-1 subtypes and mode of transmission in Cape Town, South Africa. AIDS Lond Engl. 1997 1;11(1):81–7., van HarmelenJ, WoodR, LambrickM, RybickiEP, WilliamsonAL, WilliamsonC. An association between HIV-1 subtypes and mode of transmission in Cape Town, South Africa. AIDS Lond Engl. 1997 1;11(1):81–7.

  • Rebecca Gray, A. Tatem, S. Lamers, W. Hou, O. Laeyendecker, D. Serwadda, N. Sewankambo, Ronald Gray, M. Wawer, T. Quinn, M. Goodenow, M. Salemi (2009)

    Spatial phylodynamics of HIV-1 epidemic emergence in east Africa

    AIDS, 23

  • Y Vega (2015)

    e0125699

    PLOS ONE, 10

  • (TscherningC, AlaeusA, FredrikssonR, BjörndalA, DengH, LittmanDR, et al Differences in chemokine coreceptor usage between genetic subtypes of HIV-1. Virology. 1998 2 15;241(2):181–8. doi: 10.1006/viro.1997.8980 9499793)

    TscherningC, AlaeusA, FredrikssonR, Bj�rndalA, DengH, LittmanDR, et al

    TscherningC, AlaeusA, FredrikssonR, BjörndalA, DengH, LittmanDR, et al Differences in chemokine coreceptor usage between genetic subtypes of HIV-1. Virology. 1998 2 15;241(2):181–8. doi: 10.1006/viro.1997.8980 9499793, TscherningC, AlaeusA, FredrikssonR, BjörndalA, DengH, LittmanDR, et al Differences in chemokine coreceptor usage between genetic subtypes of HIV-1. Virology. 1998 2 15;241(2):181–8. doi: 10.1006/viro.1997.8980 9499793

  • A. Fernández-García, M. Cuevas, M. Muñoz-Nieto, A. Ocampo, M. Pinilla, V. García, E. Serrano-Bengoechea, M. Lezaún, E. Delgado, M. Thomson, M. González-Galeano, G. Contreras, R. Nájera, L. Pérez-Alvárez (2009)

    Development of a panel of well-characterized human immunodeficiency virus type 1 isolates from newly diagnosed patients including acute and recent infections.

    AIDS research and human retroviruses, 25 1

  • (EsbjörnssonJ, MildM, MånssonF, NorrgrenH, MedstrandP. HIV-1 molecular epidemiology in Guinea-Bissau, West Africa: origin, demography and migrations. PloS One. 2011;6(2):e17025 doi: 10.1371/journal.pone.0017025 21365013)

    Esbj�rnssonJ, MildM, M�nssonF, NorrgrenH, MedstrandP. HIV-1 molecular epidemiology in Guinea-Bissau, West Africa: origin, demography and migrations. PloS One. 2011;6(2):e17025

    EsbjörnssonJ, MildM, MånssonF, NorrgrenH, MedstrandP. HIV-1 molecular epidemiology in Guinea-Bissau, West Africa: origin, demography and migrations. PloS One. 2011;6(2):e17025 doi: 10.1371/journal.pone.0017025 21365013, EsbjörnssonJ, MildM, MånssonF, NorrgrenH, MedstrandP. HIV-1 molecular epidemiology in Guinea-Bissau, West Africa: origin, demography and migrations. PloS One. 2011;6(2):e17025 doi: 10.1371/journal.pone.0017025 21365013

  • (Instituto Nacional de Estadística (INE) [Internet]. Available from: www.ine.es)

    Instituto Nacional de Estad�stica (INE) [Internet]. Available from: www.ine.es

    Instituto Nacional de Estadística (INE) [Internet]. Available from: www.ine.es, Instituto Nacional de Estadística (INE) [Internet]. Available from: www.ine.es

  • (Trends in Drug Resistance Prevalence in HIV-1–infected Children in Madrid (PDF Download Available) [Internet]. ResearchGate. [cited 2017 Mar 15]. Available from: https://www.researchgate.net/publication/229074263_Trends_in_Drug_Resistance_Prevalence_in_HIV-1-infected_Children_in_Madrid)

    Trends in Drug Resistance Prevalence in HIV-1?infected Children in Madrid (PDF Download Available) [Internet]. ResearchGate. [cited 2017 Mar 15]. Available from: https://www.researchgate.net/publication/229074263_Trends_in_Drug_Resistance_Prevalence_in_HIV-1-infected_Children_in_Madrid

    Trends in Drug Resistance Prevalence in HIV-1–infected Children in Madrid (PDF Download Available) [Internet]. ResearchGate. [cited 2017 Mar 15]. Available from: https://www.researchgate.net/publication/229074263_Trends_in_Drug_Resistance_Prevalence_in_HIV-1-infected_Children_in_Madrid, Trends in Drug Resistance Prevalence in HIV-1–infected Children in Madrid (PDF Download Available) [Internet]. ResearchGate. [cited 2017 Mar 15]. Available from: https://www.researchgate.net/publication/229074263_Trends_in_Drug_Resistance_Prevalence_in_HIV-1-infected_Children_in_Madrid

  • A Alaeus (1997)

    859

    AIDS Lond Engl, 11

  • G. Yebra, M. Mulder, J. Romero, C. Rodríguez, Á. Holguín (2010)

    HIV-1 non-B subtypes: High transmitted NNRTI-resistance in Spain and impaired genotypic resistance interpretation due to variability.

    Antiviral research, 85 2

  • R. Shankarappa, J. Margolick, S. Gange, A. Rodrigo, D. Upchurch, H. Farzadegan, Phalguni Gupta, C. Rinaldo, G. Learn, Xi He, Xiao-li Huang, J. Mullins (1999)

    Consistent Viral Evolutionary Changes Associated with the Progression of Human Immunodeficiency Virus Type 1 Infection

    Journal of Virology, 73

  • (González-AlbaJM, HolguínÁ, GarciaR, García-BujalanceS, AlonsoR, SuárezA, et al Molecular Surveillance of HIV-1 in Madrid, Spain: a Phylogeographic Analysis ▿. J Virol. 2011 10;85(20):10755–63. doi: 10.1128/JVI.00454-11 21795343)

    Gonz�lez-AlbaJM, Holgu�n�, GarciaR, Garc�a-BujalanceS, AlonsoR, Su�rezA, et al

    González-AlbaJM, HolguínÁ, GarciaR, García-BujalanceS, AlonsoR, SuárezA, et al Molecular Surveillance of HIV-1 in Madrid, Spain: a Phylogeographic Analysis ▿. J Virol. 2011 10;85(20):10755–63. doi: 10.1128/JVI.00454-11 21795343, González-AlbaJM, HolguínÁ, GarciaR, García-BujalanceS, AlonsoR, SuárezA, et al Molecular Surveillance of HIV-1 in Madrid, Spain: a Phylogeographic Analysis ▿. J Virol. 2011 10;85(20):10755–63. doi: 10.1128/JVI.00454-11 21795343

  • (TaylorBS, SobieszczykME, McCutchanFE, HammerSM. The challenge of HIV-1 subtype diversity. N Engl J Med. 2008 4 10;358(15):1590–602. doi: 10.1056/NEJMra0706737 18403767)

    TaylorBS, SobieszczykME, McCutchanFE, HammerSM. The challenge of HIV-1 subtype diversity. N Engl J Med. 2008

    TaylorBS, SobieszczykME, McCutchanFE, HammerSM. The challenge of HIV-1 subtype diversity. N Engl J Med. 2008 4 10;358(15):1590–602. doi: 10.1056/NEJMra0706737 18403767, TaylorBS, SobieszczykME, McCutchanFE, HammerSM. The challenge of HIV-1 subtype diversity. N Engl J Med. 2008 4 10;358(15):1590–602. doi: 10.1056/NEJMra0706737 18403767

  • D. Robertson, J. Anderson, J. Bradac, J. Carr, B. Foley, R. Funkhouser, F. Gao, B. Hahn, M. Kalish, C. Kuiken, G. Learn, T. Leitner, F. McCutchan, S. Osmanov, M. Peeters, D. Pieniążek, M. Salminen, P. Sharp, Steven Wolinsky, B. Korber (2000)

    HIV-1 Nomenclature Proposal

    Science, 288

  • (Fernández-GarcíaA, CuevasMT, VinogradovaA, RakhmanovaA, Pérez-AlvarezL, de CastroRO, et al Near full-length genome characterization of a newly identified HIV type 1 subtype F variant circulating in St. Petersburg, Russia. AIDS Res Hum Retroviruses. 2009 11;25(11):1187–91. doi: 10.1089/aid.2009.0140 19943791)

    Fern�ndez-Garc�aA, CuevasMT, VinogradovaA, RakhmanovaA, P�rez-AlvarezL, de CastroRO, et al

    Fernández-GarcíaA, CuevasMT, VinogradovaA, RakhmanovaA, Pérez-AlvarezL, de CastroRO, et al Near full-length genome characterization of a newly identified HIV type 1 subtype F variant circulating in St. Petersburg, Russia. AIDS Res Hum Retroviruses. 2009 11;25(11):1187–91. doi: 10.1089/aid.2009.0140 19943791, Fernández-GarcíaA, CuevasMT, VinogradovaA, RakhmanovaA, Pérez-AlvarezL, de CastroRO, et al Near full-length genome characterization of a newly identified HIV type 1 subtype F variant circulating in St. Petersburg, Russia. AIDS Res Hum Retroviruses. 2009 11;25(11):1187–91. doi: 10.1089/aid.2009.0140 19943791

  • DL Robertson (2000)

    55

    Science, 288

  • J Hemelaar (2011)

    679

    AIDS Lond Engl, 25

  • (Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Gateway Computing Environments Workshop (GCE), 2010. 2010. p. 1–8.)

    Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Gateway Computing Environments Workshop (GCE), 2010. 2010. p. 1?8.

    Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Gateway Computing Environments Workshop (GCE), 2010. 2010. p. 1–8., Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In: Gateway Computing Environments Workshop (GCE), 2010. 2010. p. 1–8.

  • (ParekhB, PhillipsS, GranadeTC, BaggsJ, HuDJ, RespessR. Impact of HIV type 1 subtype variation on viral RNA quantitation. AIDS Res Hum Retroviruses. 1999 1 20;15(2):133–42. doi: 10.1089/088922299311556 10029245)

    ParekhB, PhillipsS, GranadeTC, BaggsJ, HuDJ, RespessR. Impact of HIV type 1 subtype variation on viral RNA quantitation. AIDS Res Hum Retroviruses. 1999

    ParekhB, PhillipsS, GranadeTC, BaggsJ, HuDJ, RespessR. Impact of HIV type 1 subtype variation on viral RNA quantitation. AIDS Res Hum Retroviruses. 1999 1 20;15(2):133–42. doi: 10.1089/088922299311556 10029245, ParekhB, PhillipsS, GranadeTC, BaggsJ, HuDJ, RespessR. Impact of HIV type 1 subtype variation on viral RNA quantitation. AIDS Res Hum Retroviruses. 1999 1 20;15(2):133–42. doi: 10.1089/088922299311556 10029245

  • (de FelipeB, Pérez-RomeroP, Abad-FernándezM, Fernandez-CuencaF, Martinez-FernandezFJ, TrastoyM, et al Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000–2010. Virol J. 2011;8:416 doi: 10.1186/1743-422X-8-416 21871090)

    de FelipeB, P�rez-RomeroP, Abad-Fern�ndezM, Fernandez-CuencaF, Martinez-FernandezFJ, TrastoyM, et al

    de FelipeB, Pérez-RomeroP, Abad-FernándezM, Fernandez-CuencaF, Martinez-FernandezFJ, TrastoyM, et al Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000–2010. Virol J. 2011;8:416 doi: 10.1186/1743-422X-8-416 21871090, de FelipeB, Pérez-RomeroP, Abad-FernándezM, Fernandez-CuencaF, Martinez-FernandezFJ, TrastoyM, et al Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000–2010. Virol J. 2011;8:416 doi: 10.1186/1743-422X-8-416 21871090

  • (DrummondAJ, HoSYW, PhillipsMJ, RambautA. Relaxed Phylogenetics and Dating with Confidence. PLoS Biol. 2006 3 14;4(5):e88 doi: 10.1371/journal.pbio.0040088 16683862)

    DrummondAJ, HoSYW, PhillipsMJ, RambautA. Relaxed Phylogenetics and Dating with Confidence. PLoS Biol. 2006

    DrummondAJ, HoSYW, PhillipsMJ, RambautA. Relaxed Phylogenetics and Dating with Confidence. PLoS Biol. 2006 3 14;4(5):e88 doi: 10.1371/journal.pbio.0040088 16683862, DrummondAJ, HoSYW, PhillipsMJ, RambautA. Relaxed Phylogenetics and Dating with Confidence. PLoS Biol. 2006 3 14;4(5):e88 doi: 10.1371/journal.pbio.0040088 16683862

  • M. Grabowski, A. Redd (2014)

    Molecular tools for studying HIV transmission in sexual networks

    Current Opinion in HIV and AIDS, 9

  • (BrennanCA, BodelleP, CoffeyR, HarrisB, HolzmayerV, LukK-C, et al HIV global surveillance: foundation for retroviral discovery and assay development. J Med Virol. 2006;78 Suppl 1:S24–9.16622874)

    BrennanCA, BodelleP, CoffeyR, HarrisB, HolzmayerV, LukK-C, et al

    BrennanCA, BodelleP, CoffeyR, HarrisB, HolzmayerV, LukK-C, et al HIV global surveillance: foundation for retroviral discovery and assay development. J Med Virol. 2006;78 Suppl 1:S24–9.16622874, BrennanCA, BodelleP, CoffeyR, HarrisB, HolzmayerV, LukK-C, et al HIV global surveillance: foundation for retroviral discovery and assay development. J Med Virol. 2006;78 Suppl 1:S24–9.16622874

  • (2014)

    Identificación de un cluster de subtipo A que se transmite entre hombres que tienen relaciones sexuales con hombres en diversas comunidades autónomas de España

  • (GarcíaF, Pérez-CachafeiroS, GuillotV, AlvarezM, Pérez-RomeroP, Pérez-ElíasMJ, et al Transmission of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treatment naïve HIV-infected individuals (CoRIS). Antiviral Res. 2011 8;91(2):150–3. doi: 10.1016/j.antiviral.2011.05.010 21663768)

    Garc�aF, P�rez-CachafeiroS, GuillotV, AlvarezM, P�rez-RomeroP, P�rez-El�asMJ, et al

    GarcíaF, Pérez-CachafeiroS, GuillotV, AlvarezM, Pérez-RomeroP, Pérez-ElíasMJ, et al Transmission of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treatment naïve HIV-infected individuals (CoRIS). Antiviral Res. 2011 8;91(2):150–3. doi: 10.1016/j.antiviral.2011.05.010 21663768, GarcíaF, Pérez-CachafeiroS, GuillotV, AlvarezM, Pérez-RomeroP, Pérez-ElíasMJ, et al Transmission of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treatment naïve HIV-infected individuals (CoRIS). Antiviral Res. 2011 8;91(2):150–3. doi: 10.1016/j.antiviral.2011.05.010 21663768

  • (AlexievI, ShankarA, WensingAMJ, BeshkovD, ElenkovI, StoychevaM, et al Low HIV-1 transmitted drug resistance in Bulgaria against a background of high clade diversity. J Antimicrob Chemother. 2015;70(6):1874–80. doi: 10.1093/jac/dkv011 25652746)

    AlexievI, ShankarA, WensingAMJ, BeshkovD, ElenkovI, StoychevaM, et al

    AlexievI, ShankarA, WensingAMJ, BeshkovD, ElenkovI, StoychevaM, et al Low HIV-1 transmitted drug resistance in Bulgaria against a background of high clade diversity. J Antimicrob Chemother. 2015;70(6):1874–80. doi: 10.1093/jac/dkv011 25652746, AlexievI, ShankarA, WensingAMJ, BeshkovD, ElenkovI, StoychevaM, et al Low HIV-1 transmitted drug resistance in Bulgaria against a background of high clade diversity. J Antimicrob Chemother. 2015;70(6):1874–80. doi: 10.1093/jac/dkv011 25652746

  • J. Esbjörnsson, Mattias Mild, A. Audelin, J. Fonager, H. Skar, Louise Jørgensen, K. Liitsola, P. Björkman, G. Bratt, M. Gisslén, A. Sönnerborg, C. Nielsen, P. Medstrand, J. Albert (2016)

    HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries

    Virus Evolution, 2

  • D. Bennett, R. Camacho, D. Otelea, D. Kuritzkes, H. Fleury, M. Kiuchi, W. Heneine, R. Kantor, M. Jordan, J. Schapiro, A. Vandamme, P. Sandstrom, C. Boucher, D. Vijver, S. Rhee, Tommy Liu, D. Pillay, R. Shafer (2009)

    Drug Resistance Mutations for Surveillance of Transmitted HIV-1 Drug-Resistance: 2009 Update

    PLoS ONE, 4

  • Juan Galindo, M. Torres-Puente, C. Gimeno, E. Ortega, D. Navarro, M. Galindo, Laura Navarro, V. Navarro, A. Juan, J. Belda, M. Bracho, F. González-Candelas (2015)

    Expansion of the CRF19_cpx Variant in Spain.

    Journal of clinical virology : the official publication of the Pan American Society for Clinical Virology, 69

  • (BennettDE, CamachoRJ, OteleaD, KuritzkesDR, FleuryH, KiuchiM, et al Drug Resistance Mutations for Surveillance of Transmitted HIV-1 Drug-Resistance: 2009 Update. PLoS ONE [Internet]. 2009 3 6 [cited 2016 Sep 3];4(3). Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2648874/)

    BennettDE, CamachoRJ, OteleaD, KuritzkesDR, FleuryH, KiuchiM, et al

    BennettDE, CamachoRJ, OteleaD, KuritzkesDR, FleuryH, KiuchiM, et al Drug Resistance Mutations for Surveillance of Transmitted HIV-1 Drug-Resistance: 2009 Update. PLoS ONE [Internet]. 2009 3 6 [cited 2016 Sep 3];4(3). Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2648874/, BennettDE, CamachoRJ, OteleaD, KuritzkesDR, FleuryH, KiuchiM, et al Drug Resistance Mutations for Surveillance of Transmitted HIV-1 Drug-Resistance: 2009 Update. PLoS ONE [Internet]. 2009 3 6 [cited 2016 Sep 3];4(3). Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2648874/

  • S. Pérez-Parra, N. Chueca-Porcuna, Marta Álvarez-Estévez, J. Pasquau, Mohamed Omar, A. Collado, D. Vinuesa, Ana Lozano, F. García-García (2015)

    [Study of human immunodeficiency virus transmission chains in Andalusia: analysis from baseline antiretroviral resistance sequences].

    Enfermedades infecciosas y microbiologia clinica, 33 9

  • C. Mendoza, C. Garrido, E. Poveda, A. Corral, N. Zahonero, A. Treviño, L. Anta, V. Soriano (2009)

    Changes in drug resistance patterns following the introduction of HIV type 1 non-B subtypes in Spain.

    AIDS research and human retroviruses, 25 10

  • (BrennerBG, RogerM, StephensD, MoisiD, HardyI, WeinbergJ, et al Transmission Clustering Drives the Onward Spread of the HIV Epidemic Among Men Who Have Sex With Men in Quebec. J Infect Dis. 2011 10 1;204(7):1115–9. doi: 10.1093/infdis/jir468 21881127)

    BrennerBG, RogerM, StephensD, MoisiD, HardyI, WeinbergJ, et al

    BrennerBG, RogerM, StephensD, MoisiD, HardyI, WeinbergJ, et al Transmission Clustering Drives the Onward Spread of the HIV Epidemic Among Men Who Have Sex With Men in Quebec. J Infect Dis. 2011 10 1;204(7):1115–9. doi: 10.1093/infdis/jir468 21881127, BrennerBG, RogerM, StephensD, MoisiD, HardyI, WeinbergJ, et al Transmission Clustering Drives the Onward Spread of the HIV Epidemic Among Men Who Have Sex With Men in Quebec. J Infect Dis. 2011 10 1;204(7):1115–9. doi: 10.1093/infdis/jir468 21881127

  • Mark Larkin, G. Blackshields, N. Brown, R. Chenna, P. McGettigan, H. McWilliam, F. Valentin, I. Wallace, A. Wilm, R. Lopez, J. Thompson, T. Gibson, D. Higgins (2007)

    Clustal W and Clustal X version 2.0

    Bioinformatics, 23 21

  • M. Bracho, V. Sentandreu, I. Alastrué, J. Belda, A. Juan, E. Fernández-García, C. Santos, T. Zafra, T. Tasa, S. Colomina, F. González-Candelas (2014)

    Emerging trends in CRF02_AG variants transmission among men who have sex with men in Spain.

    Journal of acquired immune deficiency syndromes, 65 3

  • MA Larkin (2007)

    2947

    Bioinforma Oxf Engl, 23

  • F García (2011)

    150

    Antiviral Res, 91

  • PJ Kanki (1999)

    68

    J Infect Dis, 179

  • C Apetrei (1998)

    3534

    J Virol, 72

  • (De MendozaC, GarridoC, PovedaE, CorralA, ZahoneroN, TreviñoA, et al Changes in drug resistance patterns following the introduction of HIV type 1 non-B subtypes in Spain. AIDS Res Hum Retroviruses. 2009 10;25(10):967–72. doi: 10.1089/aid.2008.0166 19842792)

    De MendozaC, GarridoC, PovedaE, CorralA, ZahoneroN, Trevi�oA, et al

    De MendozaC, GarridoC, PovedaE, CorralA, ZahoneroN, TreviñoA, et al Changes in drug resistance patterns following the introduction of HIV type 1 non-B subtypes in Spain. AIDS Res Hum Retroviruses. 2009 10;25(10):967–72. doi: 10.1089/aid.2008.0166 19842792, De MendozaC, GarridoC, PovedaE, CorralA, ZahoneroN, TreviñoA, et al Changes in drug resistance patterns following the introduction of HIV type 1 non-B subtypes in Spain. AIDS Res Hum Retroviruses. 2009 10;25(10):967–72. doi: 10.1089/aid.2008.0166 19842792

  • (Fernández-GarcíaA, CuevasMT, Muñoz-NietoM, OcampoA, PinillaM, GarcíaV, et al Development of a panel of well-characterized human immunodeficiency virus type 1 isolates from newly diagnosed patients including acute and recent infections. AIDS Res Hum Retroviruses. 2009 1;25(1):93–102. doi: 10.1089/aid.2008.0174 19113978)

    Fern�ndez-Garc�aA, CuevasMT, Mu�oz-NietoM, OcampoA, PinillaM, Garc�aV, et al

    Fernández-GarcíaA, CuevasMT, Muñoz-NietoM, OcampoA, PinillaM, GarcíaV, et al Development of a panel of well-characterized human immunodeficiency virus type 1 isolates from newly diagnosed patients including acute and recent infections. AIDS Res Hum Retroviruses. 2009 1;25(1):93–102. doi: 10.1089/aid.2008.0174 19113978, Fernández-GarcíaA, CuevasMT, Muñoz-NietoM, OcampoA, PinillaM, GarcíaV, et al Development of a panel of well-characterized human immunodeficiency virus type 1 isolates from newly diagnosed patients including acute and recent infections. AIDS Res Hum Retroviruses. 2009 1;25(1):93–102. doi: 10.1089/aid.2008.0174 19113978

  • R Shankarappa (1999)

    10489

    J Virol, 73

  • (HuéS, BrownAE, Ragonnet-CroninM, LycettSJ, DunnDT, FearnhillE, et al Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions. AIDS Lond Engl. 2014 8 24;28(13):1967–75.)

    Hu�S, BrownAE, Ragonnet-CroninM, LycettSJ, DunnDT, FearnhillE, et al

    HuéS, BrownAE, Ragonnet-CroninM, LycettSJ, DunnDT, FearnhillE, et al Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions. AIDS Lond Engl. 2014 8 24;28(13):1967–75., HuéS, BrownAE, Ragonnet-CroninM, LycettSJ, DunnDT, FearnhillE, et al Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions. AIDS Lond Engl. 2014 8 24;28(13):1967–75.

  • Y Taniguchi (2002)

    79

    AIDS Res Hum Retroviruses, 18

  • AW Artenstein (1995)

    805

    J Infect Dis, 171

  • S. Hué, A. Brown, M. Ragonnet-Cronin, S. Lycett, D. Dunn, E. Fearnhill, D. Dolling, A. Pozniak, D. Pillay, V. Delpech, A. Brown (2014)

    Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions

    AIDS, 28

  • R. Kouyos, V. Wyl, S. Yerly, J. Böni, P. Taffé, C. Shah, P. Bürgisser, T. Klimkait, R. Weber, B. Hirschel, M. Cavassini, H. Furrer, M. Battegay, P. Vernazza, E. Bernasconi, M. Rickenbach, B. Ledergerber, S. Bonhoeffer, H. Günthard (2010)

    Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland.

    The Journal of infectious diseases, 201 10

  • B Parekh (1999)

    133

    AIDS Res Hum Retroviruses, 15

  • A. Alaeus, K. Lidman, A. Sönnerborg, J. Albert (1997)

    Subtype‐specific problems with quantification of plasma HIV‐1 RNA

    AIDS, 11

  • A. Brown, S. Lycett, L. Weinert, G. Hughes, E. Fearnhill, D. Dunn (2011)

    Transmission Network Parameters Estimated From HIV Sequences for a Nationwide Epidemic

    The Journal of Infectious Diseases, 204

  • (Patiño GalindoJA, Torres-PuenteM, GimenoC, OrtegaE, NavarroD, GalindoMJ, et al Expansion of the CRF19_cpx Variant in Spain. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2015 8;69:146–9.)

    Pati�o GalindoJA, Torres-PuenteM, GimenoC, OrtegaE, NavarroD, GalindoMJ, et al

    Patiño GalindoJA, Torres-PuenteM, GimenoC, OrtegaE, NavarroD, GalindoMJ, et al Expansion of the CRF19_cpx Variant in Spain. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2015 8;69:146–9., Patiño GalindoJA, Torres-PuenteM, GimenoC, OrtegaE, NavarroD, GalindoMJ, et al Expansion of the CRF19_cpx Variant in Spain. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2015 8;69:146–9.

  • (KankiPJ, HamelDJ, SankaléJL, HsiehC c, ThiorI, BarinF, et al Human immunodeficiency virus type 1 subtypes differ in disease progression. J Infect Dis. 1999 1;179(1):68–73. doi: 10.1086/314557 9841824)

    KankiPJ, HamelDJ, Sankal�JL, HsiehC c, ThiorI, BarinF, et al

    KankiPJ, HamelDJ, SankaléJL, HsiehC c, ThiorI, BarinF, et al Human immunodeficiency virus type 1 subtypes differ in disease progression. J Infect Dis. 1999 1;179(1):68–73. doi: 10.1086/314557 9841824, KankiPJ, HamelDJ, SankaléJL, HsiehC c, ThiorI, BarinF, et al Human immunodeficiency virus type 1 subtypes differ in disease progression. J Infect Dis. 1999 1;179(1):68–73. doi: 10.1086/314557 9841824

  • (HillisDM, BullJJ. An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phylogenetic Analysis. Syst Biol. 1993 6 1;42(2):182–92.)

    HillisDM, BullJJ. An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phylogenetic Analysis. Syst Biol. 1993

    HillisDM, BullJJ. An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phylogenetic Analysis. Syst Biol. 1993 6 1;42(2):182–92., HillisDM, BullJJ. An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phylogenetic Analysis. Syst Biol. 1993 6 1;42(2):182–92.

  • C. Apetrei, D. Descamps, G. Collin, Ibtisam Loussert-Ajaka, F. Damond, M. Duca, F. Simon, F. Brun-vezinet (1998)

    Human Immunodeficiency Virus Type 1 Subtype F Reverse Transcriptase Sequence and Drug Susceptibility

    Journal of Virology, 72

  • (YebraG, de MulderM, MartínL, RodríguezC, LabargaP, VicianaI, et al Most HIV type 1 non-B infections in the Spanish cohort of antiretroviral treatment-naïve HIV-infected patients (CoRIS) are due to recombinant viruses. J Clin Microbiol. 2012 2;50(2):407–13. doi: 10.1128/JCM.05798-11 22162552)

    YebraG, de MulderM, Mart�nL, Rodr�guezC, LabargaP, VicianaI, et al

    YebraG, de MulderM, MartínL, RodríguezC, LabargaP, VicianaI, et al Most HIV type 1 non-B infections in the Spanish cohort of antiretroviral treatment-naïve HIV-infected patients (CoRIS) are due to recombinant viruses. J Clin Microbiol. 2012 2;50(2):407–13. doi: 10.1128/JCM.05798-11 22162552, YebraG, de MulderM, MartínL, RodríguezC, LabargaP, VicianaI, et al Most HIV type 1 non-B infections in the Spanish cohort of antiretroviral treatment-naïve HIV-infected patients (CoRIS) are due to recombinant viruses. J Clin Microbiol. 2012 2;50(2):407–13. doi: 10.1128/JCM.05798-11 22162552

  • (CuevasMT, Muñoz-NietoM, ThomsonMM, DelgadoE, IribarrenJA, CillaG, et al HIV-1 Transmission Cluster With T215D Revertant Mutation Among Newly Diagnosed Patients From the Basque Country, Spain: JAIDS J Acquir Immune Defic Syndr. 2009 5;51(1):99–103. doi: 10.1097/QAI.0b013e318199063e 19282784)

    CuevasMT, Mu�oz-NietoM, ThomsonMM, DelgadoE, IribarrenJA, CillaG, et al

    CuevasMT, Muñoz-NietoM, ThomsonMM, DelgadoE, IribarrenJA, CillaG, et al HIV-1 Transmission Cluster With T215D Revertant Mutation Among Newly Diagnosed Patients From the Basque Country, Spain: JAIDS J Acquir Immune Defic Syndr. 2009 5;51(1):99–103. doi: 10.1097/QAI.0b013e318199063e 19282784, CuevasMT, Muñoz-NietoM, ThomsonMM, DelgadoE, IribarrenJA, CillaG, et al HIV-1 Transmission Cluster With T215D Revertant Mutation Among Newly Diagnosed Patients From the Basque Country, Spain: JAIDS J Acquir Immune Defic Syndr. 2009 5;51(1):99–103. doi: 10.1097/QAI.0b013e318199063e 19282784

  • M. Gilbert, A. Rambaut, G. Wlasiuk, T. Spira, A. Pitchenik, M. Worobey (2007)

    The emergence of HIV/AIDS in the Americas and beyond

    Proceedings of the National Academy of Sciences, 104

  • (ThomsonMM, Fernández-GarcíaA, DelgadoE, VegaY, Díez-FuertesF, Sánchez-MartínezM, et al Rapid expansion of a HIV-1 subtype F cluster of recent origin among men who have sex with men in Galicia, Spain. J Acquir Immune Defic Syndr 1999. 2012 3 1;59(3):e49–51.)

    ThomsonMM, Fern�ndez-Garc�aA, DelgadoE, VegaY, D�ez-FuertesF, S�nchez-Mart�nezM, et al

    ThomsonMM, Fernández-GarcíaA, DelgadoE, VegaY, Díez-FuertesF, Sánchez-MartínezM, et al Rapid expansion of a HIV-1 subtype F cluster of recent origin among men who have sex with men in Galicia, Spain. J Acquir Immune Defic Syndr 1999. 2012 3 1;59(3):e49–51., ThomsonMM, Fernández-GarcíaA, DelgadoE, VegaY, Díez-FuertesF, Sánchez-MartínezM, et al Rapid expansion of a HIV-1 subtype F cluster of recent origin among men who have sex with men in Galicia, Spain. J Acquir Immune Defic Syndr 1999. 2012 3 1;59(3):e49–51.

  • (ApetreiC, DescampsD, CollinG, Loussert-AjakaI, DamondF, DucaM, et al Human immunodeficiency virus type 1 subtype F reverse transcriptase sequence and drug susceptibility. J Virol. 1998 5;72(5):3534–8. 9557632)

    ApetreiC, DescampsD, CollinG, Loussert-AjakaI, DamondF, DucaM, et al

    ApetreiC, DescampsD, CollinG, Loussert-AjakaI, DamondF, DucaM, et al Human immunodeficiency virus type 1 subtype F reverse transcriptase sequence and drug susceptibility. J Virol. 1998 5;72(5):3534–8. 9557632, ApetreiC, DescampsD, CollinG, Loussert-AjakaI, DamondF, DucaM, et al Human immunodeficiency virus type 1 subtype F reverse transcriptase sequence and drug susceptibility. J Virol. 1998 5;72(5):3534–8. 9557632

  • G. Yebra, M. Mulder, M. Pérez-Elías, J. Pérez-Molina, J. Galán, J. Llenas-García, S. Moreno, Á. Holguín (2011)

    Increase of Transmitted Drug Resistance among HIV-Infected Sub-Saharan Africans Residing in Spain in Contrast to the Native Population

    PLoS ONE, 6

  • (YebraG, de MulderM, del RomeroJ, RodríguezC, HolguínA. HIV-1 non-B subtypes: High transmitted NNRTI-resistance in Spain and impaired genotypic resistance interpretation due to variability. Antiviral Res. 2010 2;85(2):409–17. doi: 10.1016/j.antiviral.2009.11.010 20004217)

    YebraG, de MulderM, del RomeroJ, Rodr�guezC, Holgu�nA. HIV-1 non-B subtypes: High transmitted NNRTI-resistance in Spain and impaired genotypic resistance interpretation due to variability. Antiviral Res. 2010

    YebraG, de MulderM, del RomeroJ, RodríguezC, HolguínA. HIV-1 non-B subtypes: High transmitted NNRTI-resistance in Spain and impaired genotypic resistance interpretation due to variability. Antiviral Res. 2010 2;85(2):409–17. doi: 10.1016/j.antiviral.2009.11.010 20004217, YebraG, de MulderM, del RomeroJ, RodríguezC, HolguínA. HIV-1 non-B subtypes: High transmitted NNRTI-resistance in Spain and impaired genotypic resistance interpretation due to variability. Antiviral Res. 2010 2;85(2):409–17. doi: 10.1016/j.antiviral.2009.11.010 20004217

  • Mark Miller, W. Pfeiffer, Terri Schwartz (2010)

    Creating the CIPRES Science Gateway for inference of large phylogenetic trees

    2010 Gateway Computing Environments Workshop (GCE)

  • (VegaY, DelgadoE, Fernández-GarcíaA, CuevasMT, ThomsonMM, MonteroV, et al Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004–2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations. PLOS ONE. 2015 5 26;10(5):e0125699 doi: 10.1371/journal.pone.0125699 26010948)

    VegaY, DelgadoE, Fern�ndez-Garc�aA, CuevasMT, ThomsonMM, MonteroV, et al

    VegaY, DelgadoE, Fernández-GarcíaA, CuevasMT, ThomsonMM, MonteroV, et al Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004–2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations. PLOS ONE. 2015 5 26;10(5):e0125699 doi: 10.1371/journal.pone.0125699 26010948, VegaY, DelgadoE, Fernández-GarcíaA, CuevasMT, ThomsonMM, MonteroV, et al Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004–2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations. PLOS ONE. 2015 5 26;10(5):e0125699 doi: 10.1371/journal.pone.0125699 26010948

  • B. Brenner, M. Roger, David Stephens, D. Moisi, I. Hardy, J. Weinberg, Reuven Turgel, H. Charest, J. Koopman, M. Wainberg (2011)

    Transmission clustering drives the onward spread of the HIV epidemic among men who have sex with men in Quebec.

    The Journal of infectious diseases, 204 7

  • G Yebra (2012)

    407

    J Clin Microbiol, 50

  • N. Vidal, C. Mulanga, S. Bazepeo, J. Mwamba, J. Tshimpaka, M. Kashi, N. Mama, C. Laurent, F. Lepira, E. Delaporte, M. Peeters (2005)

    Distribution of HIV-1 Variants in the Democratic Republic of Congo Suggests Increase of Subtype C in Kinshasa Between 1997 and 2002

    JAIDS Journal of Acquired Immune Deficiency Syndromes, 40

  • (GilbertMTP, RambautA, WlasiukG, SpiraTJ, PitchenikAE, WorobeyM. The emergence of HIV/AIDS in the Americas and beyond. Proc Natl Acad Sci U S A. 2007 11 20;104(47):18566–70. doi: 10.1073/pnas.0705329104 17978186)

    GilbertMTP, RambautA, WlasiukG, SpiraTJ, PitchenikAE, WorobeyM. The emergence of HIV/AIDS in the Americas and beyond. Proc Natl Acad Sci U S A. 2007

    GilbertMTP, RambautA, WlasiukG, SpiraTJ, PitchenikAE, WorobeyM. The emergence of HIV/AIDS in the Americas and beyond. Proc Natl Acad Sci U S A. 2007 11 20;104(47):18566–70. doi: 10.1073/pnas.0705329104 17978186, GilbertMTP, RambautA, WlasiukG, SpiraTJ, PitchenikAE, WorobeyM. The emergence of HIV/AIDS in the Americas and beyond. Proc Natl Acad Sci U S A. 2007 11 20;104(47):18566–70. doi: 10.1073/pnas.0705329104 17978186

  • (BrachoMA, SentandreuV, AlastruéI, BeldaJ, JuanA, Fernández-GarcíaE, et al Emerging trends in CRF02_AG variants transmission among men who have sex with men in Spain. J Acquir Immune Defic Syndr 1999. 2014 3 1;65(3):e130–3.)

    BrachoMA, SentandreuV, Alastru�I, BeldaJ, JuanA, Fern�ndez-Garc�aE, et al

    BrachoMA, SentandreuV, AlastruéI, BeldaJ, JuanA, Fernández-GarcíaE, et al Emerging trends in CRF02_AG variants transmission among men who have sex with men in Spain. J Acquir Immune Defic Syndr 1999. 2014 3 1;65(3):e130–3., BrachoMA, SentandreuV, AlastruéI, BeldaJ, JuanA, Fernández-GarcíaE, et al Emerging trends in CRF02_AG variants transmission among men who have sex with men in Spain. J Acquir Immune Defic Syndr 1999. 2014 3 1;65(3):e130–3.

  • B. Parekh, S. Phillips, T. Granade, J. Baggs, D. Hu, R. Respess (1999)

    Impact of HIV type 1 subtype variation on viral RNA quantitation.

    AIDS research and human retroviruses, 15 2

  • J. Hemelaar, E. Gouws, P. Ghys, S. Osmanov (2011)

    Global trends in molecular epidemiology of HIV-1 during 2000–2007

    AIDS, 25

  • Y. Taniguchi, J. Takehisa, B. Bikandou, I. Mboudjeka, Marie-Yvonne N'Doundou-N'Kodia, Obengui, Michel M'Pandi, P. M'pele, Y. Harada, E. Ido, M. Hayami, H. Ichimura, H. Parra (2002)

    Genetic subtypes of HIV type 1 based on the vpu/env sequences in the Republic of Congo.

    AIDS research and human retroviruses, 18 1

  • (TaniguchiY, TakehisaJ, BikandouB, MboudjekaI, N’Doundou-N’KodiaM-Y, Obenguinull, et al Genetic subtypes of HIV type 1 based on the vpu/env sequences in the Republic of Congo. AIDS Res Hum Retroviruses. 2002 1 1;18(1):79–83. doi: 10.1089/088922202753394745 11804559)

    TaniguchiY, TakehisaJ, BikandouB, MboudjekaI, N?Doundou-N?KodiaM-Y, Obenguinull, et al

    TaniguchiY, TakehisaJ, BikandouB, MboudjekaI, N’Doundou-N’KodiaM-Y, Obenguinull, et al Genetic subtypes of HIV type 1 based on the vpu/env sequences in the Republic of Congo. AIDS Res Hum Retroviruses. 2002 1 1;18(1):79–83. doi: 10.1089/088922202753394745 11804559, TaniguchiY, TakehisaJ, BikandouB, MboudjekaI, N’Doundou-N’KodiaM-Y, Obenguinull, et al Genetic subtypes of HIV type 1 based on the vpu/env sequences in the Republic of Congo. AIDS Res Hum Retroviruses. 2002 1 1;18(1):79–83. doi: 10.1089/088922202753394745 11804559

  • G. Yebra, M. Mulder, Leticia Martín, C. Rodríguez, P. Labarga, I. Viciana, J. Berenguer, M. Alemán, J. Pineda, Federico García, Á. Holguín (2011)

    Most HIV Type 1 Non-B Infections in the Spanish Cohort of Antiretroviral Treatment-Naïve HIV-Infected Patients (CoRIS) Are Due to Recombinant Viruses

    Journal of Clinical Microbiology, 50

  • MT Cuevas (2009)

    99

    JAIDS J Acquir Immune Defic Syndr, 51

  • E. Delgado, M. Cuevas, Francisco Domínguez, Y. Vega, M. Cabello, A. Fernández-García, M. Pérez‐Losada, M. Castro, V. Montero, Mónica Sánchez, A. Mariño, H. Álvarez, Patricia Ordóñez, A. Ocampo, C. Miralles, S. Pérez-Castro, M. López-Álvarez, R. Rodríguez, M. Trigo, J. Diz-Arén, C. Hinojosa, P. Bachiller, S. Hernáez-Crespo, R. Cisterna, E. Garduño, L. Pérez-Alvárez, M. Thomson (2015)

    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

    PLoS ONE, 10

  • J Esbjörnsson (2011)

    e17025

    PloS One, 6

  • (GrayRR, TatemAJ, LamersS, HouW, LaeyendeckerO, SerwaddaD, et al Spatial phylodynamics of HIV-1 epidemic emergence in east Africa. AIDS Lond Engl. 2009 9 10;23(14):F9–17.)

    GrayRR, TatemAJ, LamersS, HouW, LaeyendeckerO, SerwaddaD, et al

    GrayRR, TatemAJ, LamersS, HouW, LaeyendeckerO, SerwaddaD, et al Spatial phylodynamics of HIV-1 epidemic emergence in east Africa. AIDS Lond Engl. 2009 9 10;23(14):F9–17., GrayRR, TatemAJ, LamersS, HouW, LaeyendeckerO, SerwaddaD, et al Spatial phylodynamics of HIV-1 epidemic emergence in east Africa. AIDS Lond Engl. 2009 9 10;23(14):F9–17.

  • (RobertsonDL, AndersonJP, BradacJA, CarrJK, FoleyB, FunkhouserRK, et al HIV-1 nomenclature proposal. Science. 2000 4 7;288(5463):55–6. 10766634)

    RobertsonDL, AndersonJP, BradacJA, CarrJK, FoleyB, FunkhouserRK, et al

    RobertsonDL, AndersonJP, BradacJA, CarrJK, FoleyB, FunkhouserRK, et al HIV-1 nomenclature proposal. Science. 2000 4 7;288(5463):55–6. 10766634, RobertsonDL, AndersonJP, BradacJA, CarrJK, FoleyB, FunkhouserRK, et al HIV-1 nomenclature proposal. Science. 2000 4 7;288(5463):55–6. 10766634

  • (DrummondAJ, SuchardMA, XieD, RambautA. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012 8;29(8):1969–73. doi: 10.1093/molbev/mss075 22367748)

    DrummondAJ, SuchardMA, XieD, RambautA. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012

    DrummondAJ, SuchardMA, XieD, RambautA. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012 8;29(8):1969–73. doi: 10.1093/molbev/mss075 22367748, DrummondAJ, SuchardMA, XieD, RambautA. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012 8;29(8):1969–73. doi: 10.1093/molbev/mss075 22367748

  • A. Drummond, A. Rambaut, B. Shapiro, O. Pybus (2005)

    Bayesian coalescent inference of past population dynamics from molecular sequences.

    Molecular biology and evolution, 22 5

  • D. Hillis, J. Bull (1993)

    An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phylogenetic Analysis

    Systematic Biology, 42

  • C Tscherning (1998)

    181

    Virology, 241

  • (ArtensteinAW, VanCottTC, MascolaJR, CarrJK, HegerichPA, GayweeJ, et al Dual infection with human immunodeficiency virus type 1 of distinct envelope subtypes in humans. J Infect Dis. 1995 4;171(4):805–10. 7706806)

    ArtensteinAW, VanCottTC, MascolaJR, CarrJK, HegerichPA, GayweeJ, et al

    ArtensteinAW, VanCottTC, MascolaJR, CarrJK, HegerichPA, GayweeJ, et al Dual infection with human immunodeficiency virus type 1 of distinct envelope subtypes in humans. J Infect Dis. 1995 4;171(4):805–10. 7706806, ArtensteinAW, VanCottTC, MascolaJR, CarrJK, HegerichPA, GayweeJ, et al Dual infection with human immunodeficiency virus type 1 of distinct envelope subtypes in humans. J Infect Dis. 1995 4;171(4):805–10. 7706806

  • C. Brennan, P. Bodelle, Rebecca Coffey, B. Harris, V. Holzmayer, K. Luk, P. Swanson, J. Yamaguchi, A. Vallari, S. Devare, G. Schochetman, J. Hackett (2006)

    HIV global surveillance: Foundation for retroviral discovery and assay development

    Journal of Medical Virology, 78

  • M. Cuevas, M. Muñoz-Nieto, M. Thomson, E. Delgado, J. Iribarren, G. Cilla, A. Fernández-García, J. Santamaría, M. Lezaún, Laura Jiménez, L. López-Soria, M. Sota, G. Contreras, R. Nájera, L. Pérez-Alvárez (2009)

    HIV-1 Transmission Cluster With T215D Revertant Mutation Among Newly Diagnosed Patients From the Basque Country, Spain

    JAIDS Journal of Acquired Immune Deficiency Syndromes, 51

  • A. Drummond, S. Ho, M. Phillips, A. Rambaut (2006)

    Relaxed Phylogenetics and Dating with Confidence

    PLoS Biology, 4

  • M. Álvarez, Federico García, N. Martínez, Fernando García, C. Bernal, Carmen Vela, G. Angulo, J. Quero (2003)

    Introduction of HIV type 1 non‐B subtypes into Eastern Andalusia through immigration

    Journal of Medical Virology, 70

  • CA Brennan (2006)

    S24

    J Med Virol, 78

  • JL Aldous (2012)

    1135

    Clin Infect Dis Off Publ Infect Dis Soc Am, 55

  • (AldousJL, PondSK, PoonA, JainS, QinH, KahnJS, et al Characterizing HIV transmission networks across the United States. Clin Infect Dis Off Publ Infect Dis Soc Am. 2012 10;55(8):1135–43.)

    AldousJL, PondSK, PoonA, JainS, QinH, KahnJS, et al

    AldousJL, PondSK, PoonA, JainS, QinH, KahnJS, et al Characterizing HIV transmission networks across the United States. Clin Infect Dis Off Publ Infect Dis Soc Am. 2012 10;55(8):1135–43., AldousJL, PondSK, PoonA, JainS, QinH, KahnJS, et al Characterizing HIV transmission networks across the United States. Clin Infect Dis Off Publ Infect Dis Soc Am. 2012 10;55(8):1135–43.

  • P. Kanki, D. Hamel, J. Sankalé, Chung-cheng Hsieh, I. Thior, Francis Barin, Stephen Woodcock, A. Guéye-Ndiaye, Er Zhang, Monty Montano, T. Siby, R. Marlink, Ibra Ndoye, Myron Essex, S. Mboup (1999)

    Human immunodeficiency virus type 1 subtypes differ in disease progression.

    The Journal of infectious diseases, 179 1

  • (HemelaarJ, GouwsE, GhysPD, OsmanovS, WHO-UNAIDS Network for HIV Isolation and Characterisation. Global trends in molecular epidemiology of HIV-1 during 2000–2007. AIDS Lond Engl. 2011 3 13;25(5):679–89.)

    HemelaarJ, GouwsE, GhysPD, OsmanovS, WHO-UNAIDS Network for HIV Isolation and Characterisation. Global trends in molecular epidemiology of HIV-1 during 2000?2007. AIDS Lond Engl. 2011

    HemelaarJ, GouwsE, GhysPD, OsmanovS, WHO-UNAIDS Network for HIV Isolation and Characterisation. Global trends in molecular epidemiology of HIV-1 during 2000–2007. AIDS Lond Engl. 2011 3 13;25(5):679–89., HemelaarJ, GouwsE, GhysPD, OsmanovS, WHO-UNAIDS Network for HIV Isolation and Characterisation. Global trends in molecular epidemiology of HIV-1 during 2000–2007. AIDS Lond Engl. 2011 3 13;25(5):679–89.

  • K. Lole, R. Bollinger, R. Paranjape, D. Gadkari, S. Kulkarni, Nicole Novak, Roxann Ingersoll, H. Sheppard, S. Ray (1999)

    Full-Length Human Immunodeficiency Virus Type 1 Genomes from Subtype C-Infected Seroconverters in India, with Evidence of Intersubtype Recombination

    Journal of Virology, 73

  • M. Thomson, A. Fernández-García, E. Delgado, Y. Vega, F. Díez-Fuertes, M. Sánchez-Martínez, M. Pinilla, M. Castro, A. Mariño, Patricia Ordóñez, A. Ocampo, Alfredo Silva, S. Pérez-Castro, M. López-Álvarez, M. Trigo, L. Pérez-Alvárez (2012)

    Rapid expansion of a HIV-1 subtype F cluster of recent origin among men who have sex with men in Galicia, Spain.

    Journal of acquired immune deficiency syndromes, 59 3

  • BS Taylor (2008)

    1590

    N Engl J Med, 358

  • MK Grabowski (2014)

    126

    Curr Opin HIV AIDS, 9

  • JM González-Alba (2011)

    10755

    J Virol, 85

  • (DelgadoE, CuevasMT, DomínguezF, VegaY, CabelloM, Fernández-GarcíaA, et al 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. PloS One. 2015;10(11):e0143325 doi: 10.1371/journal.pone.0143325 26599410)

    DelgadoE, CuevasMT, Dom�nguezF, VegaY, CabelloM, Fern�ndez-Garc�aA, et al

    DelgadoE, CuevasMT, DomínguezF, VegaY, CabelloM, Fernández-GarcíaA, et al 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. PloS One. 2015;10(11):e0143325 doi: 10.1371/journal.pone.0143325 26599410, DelgadoE, CuevasMT, DomínguezF, VegaY, CabelloM, Fernández-GarcíaA, et al 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. PloS One. 2015;10(11):e0143325 doi: 10.1371/journal.pone.0143325 26599410

  • (AlvarezM, GarcíaF, MartínezNM, GarcíaF, BernalC, VelaCM, et al Introduction of HIV type 1 non-B subtypes into Eastern Andalusia through immigration. J Med Virol. 2003 5;70(1):10–3. doi: 10.1002/jmv.10368 12629637)

    AlvarezM, Garc�aF, Mart�nezNM, Garc�aF, BernalC, VelaCM, et al

    AlvarezM, GarcíaF, MartínezNM, GarcíaF, BernalC, VelaCM, et al Introduction of HIV type 1 non-B subtypes into Eastern Andalusia through immigration. J Med Virol. 2003 5;70(1):10–3. doi: 10.1002/jmv.10368 12629637, AlvarezM, GarcíaF, MartínezNM, GarcíaF, BernalC, VelaCM, et al Introduction of HIV type 1 non-B subtypes into Eastern Andalusia through immigration. J Med Virol. 2003 5;70(1):10–3. doi: 10.1002/jmv.10368 12629637

  • (ShankarappaR, MargolickJB, GangeSJ, RodrigoAG, UpchurchD, FarzadeganH, et al Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol. 1999 12;73(12):10489–502. 10559367)

    ShankarappaR, MargolickJB, GangeSJ, RodrigoAG, UpchurchD, FarzadeganH, et al

    ShankarappaR, MargolickJB, GangeSJ, RodrigoAG, UpchurchD, FarzadeganH, et al Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol. 1999 12;73(12):10489–502. 10559367, ShankarappaR, MargolickJB, GangeSJ, RodrigoAG, UpchurchD, FarzadeganH, et al Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol. 1999 12;73(12):10489–502. 10559367

  • (HolguínA, de MulderM, YebraG, LópezM, SorianoV. Increase of non-B subtypes and recombinants among newly diagnosed HIV-1 native Spaniards and immigrants in Spain. Curr HIV Res. 2008 6;6(4):327–34. 18691031)

    Holgu�nA, de MulderM, YebraG, L�pezM, SorianoV. Increase of non-B subtypes and recombinants among newly diagnosed HIV-1 native Spaniards and immigrants in Spain. Curr HIV Res. 2008

    HolguínA, de MulderM, YebraG, LópezM, SorianoV. Increase of non-B subtypes and recombinants among newly diagnosed HIV-1 native Spaniards and immigrants in Spain. Curr HIV Res. 2008 6;6(4):327–34. 18691031, HolguínA, de MulderM, YebraG, LópezM, SorianoV. Increase of non-B subtypes and recombinants among newly diagnosed HIV-1 native Spaniards and immigrants in Spain. Curr HIV Res. 2008 6;6(4):327–34. 18691031

  • A. Fernández-García, M. Cuevas, Anna Vinogradova, A. Rakhmanova, L. Pérez-Alvárez, R. Castro, S. Osmanov, M. Thomson (2009)

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

    AIDS research and human retroviruses, 25 11

  • (GrabowskiMK, ReddAD. Molecular tools for studying HIV transmission in sexual networks. Curr Opin HIV AIDS. 2014 3;9(2):126–33. doi: 10.1097/COH.0000000000000040 24384502)

    GrabowskiMK, ReddAD. Molecular tools for studying HIV transmission in sexual networks. Curr Opin HIV AIDS. 2014

    GrabowskiMK, ReddAD. Molecular tools for studying HIV transmission in sexual networks. Curr Opin HIV AIDS. 2014 3;9(2):126–33. doi: 10.1097/COH.0000000000000040 24384502, GrabowskiMK, ReddAD. Molecular tools for studying HIV transmission in sexual networks. Curr Opin HIV AIDS. 2014 3;9(2):126–33. doi: 10.1097/COH.0000000000000040 24384502

  • (LarkinMA, BlackshieldsG, BrownNP, ChennaR, McGettiganPA, McWilliamH, et al Clustal W and Clustal X version 2.0. Bioinforma Oxf Engl. 2007 11 1;23(21):2947–8.)

    LarkinMA, BlackshieldsG, BrownNP, ChennaR, McGettiganPA, McWilliamH, et al

    LarkinMA, BlackshieldsG, BrownNP, ChennaR, McGettiganPA, McWilliamH, et al Clustal W and Clustal X version 2.0. Bioinforma Oxf Engl. 2007 11 1;23(21):2947–8., LarkinMA, BlackshieldsG, BrownNP, ChennaR, McGettiganPA, McWilliamH, et al Clustal W and Clustal X version 2.0. Bioinforma Oxf Engl. 2007 11 1;23(21):2947–8.

  • (LoleKS, BollingerRC, ParanjapeRS, GadkariD, KulkarniSS, NovakNG, et al Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999 1;73(1):152–60. 9847317)

    LoleKS, BollingerRC, ParanjapeRS, GadkariD, KulkarniSS, NovakNG, et al

    LoleKS, BollingerRC, ParanjapeRS, GadkariD, KulkarniSS, NovakNG, et al Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999 1;73(1):152–60. 9847317, LoleKS, BollingerRC, ParanjapeRS, GadkariD, KulkarniSS, NovakNG, et al Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999 1;73(1):152–60. 9847317

  • (DrummondAJ, RambautA, ShapiroB, PybusOG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005 5;22(5):1185–92. doi: 10.1093/molbev/msi103 15703244)

    DrummondAJ, RambautA, ShapiroB, PybusOG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005

    DrummondAJ, RambautA, ShapiroB, PybusOG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005 5;22(5):1185–92. doi: 10.1093/molbev/msi103 15703244, DrummondAJ, RambautA, ShapiroB, PybusOG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005 5;22(5):1185–92. doi: 10.1093/molbev/msi103 15703244

  • (Pérez-ParraS, Chueca-PorcunaN, Álvarez-EstevezM, PasquauJ, OmarM, ColladoA, et al [Study of human immunodeficiency virus transmission chains in Andalusia: Analysis from baseline antiretroviral resistance sequences.]. Enferm Infecc Microbiol Clin. 2015 1 31;)

    P�rez-ParraS, Chueca-PorcunaN, �lvarez-EstevezM, PasquauJ, OmarM, ColladoA, et al [Study of human immunodeficiency virus transmission chains in Andalusia: Analysis from baseline antiretroviral resistance sequences.]. Enferm Infecc Microbiol Clin. 2015

    Pérez-ParraS, Chueca-PorcunaN, Álvarez-EstevezM, PasquauJ, OmarM, ColladoA, et al [Study of human immunodeficiency virus transmission chains in Andalusia: Analysis from baseline antiretroviral resistance sequences.]. Enferm Infecc Microbiol Clin. 2015 1 31;, Pérez-ParraS, Chueca-PorcunaN, Álvarez-EstevezM, PasquauJ, OmarM, ColladoA, et al [Study of human immunodeficiency virus transmission chains in Andalusia: Analysis from baseline antiretroviral resistance sequences.]. Enferm Infecc Microbiol Clin. 2015 1 31;

  • (Delgado E, Cuevas MT, Vega Y, Montero V, Sánchez M, Carrera C, et al. Identificación de un cluster de subtipo A que se transmite entre hombres que tienen relaciones sexuales con hombres en diversas comunidades autónomas de España. In Málaga: VI Congreso Nacional de GESIDA y 8.a Reunión Docente de la RIS (SEIMC); 2014. p. 12. Available from: http://www.gesida-seimc.org/contenidos/congresos/anteriores/2014/gesida2014-VIcongresocomunicaciones.pdf)

    Delgado E, Cuevas MT, Vega Y, Montero V, S�nchez M, Carrera C, et al. Identificaci�n de un cluster de subtipo A que se transmite entre hombres que tienen relaciones sexuales con hombres en diversas comunidades aut�nomas de Espa�a. In M�laga: VI Congreso Nacional de GESIDA y 8.a Reuni�n Docente de la RIS (SEIMC); 2014. p. 12. Available from: http://www.gesida-seimc.org/contenidos/congresos/anteriores/2014/gesida2014-VIcongresocomunicaciones.pdf

    Delgado E, Cuevas MT, Vega Y, Montero V, Sánchez M, Carrera C, et al. Identificación de un cluster de subtipo A que se transmite entre hombres que tienen relaciones sexuales con hombres en diversas comunidades autónomas de España. In Málaga: VI Congreso Nacional de GESIDA y 8.a Reunión Docente de la RIS (SEIMC); 2014. p. 12. Available from: http://www.gesida-seimc.org/contenidos/congresos/anteriores/2014/gesida2014-VIcongresocomunicaciones.pdf, Delgado E, Cuevas MT, Vega Y, Montero V, Sánchez M, Carrera C, et al. Identificación de un cluster de subtipo A que se transmite entre hombres que tienen relaciones sexuales con hombres en diversas comunidades autónomas de España. In Málaga: VI Congreso Nacional de GESIDA y 8.a Reunión Docente de la RIS (SEIMC); 2014. p. 12. Available from: http://www.gesida-seimc.org/contenidos/congresos/anteriores/2014/gesida2014-VIcongresocomunicaciones.pdf

  • A. Artenstein, T. Vancott, J. Mascola, J. Carr, P. Hegerich, J. Gaywee, E. Sanders-Buell, M. Robb, D. Dayhoff, S. Thitivichianlert, S. Nitayaphan, J. McNeil, D. Birx, R. Michael, D. Burke, F. McCutchan (1995)

    Dual infection with human immunodeficiency virus type 1 of distinct envelope subtypes in humans.

    The Journal of infectious diseases, 171 4

  • RR Gray (2009)

    F9

    AIDS Lond Engl, 23

  • Y. Vega, E. Delgado, A. Fernández-García, M. Cuevas, M. Thomson, V. Montero, Mónica Sánchez, A. Sánchez, L. Pérez-Alvárez (2015)

    Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004-2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations

    PLoS ONE, 10

  • (VidalN, MulangaC, BazepeoSE, MwambaJK, TshimpakaJ-W, KashiM, et al Distribution of HIV-1 variants in the Democratic Republic of Congo suggests increase of subtype C in Kinshasa between 1997 and 2002. J Acquir Immune Defic Syndr 1999. 2005 12 1;40(4):456–62.)

    VidalN, MulangaC, BazepeoSE, MwambaJK, TshimpakaJ-W, KashiM, et al

    VidalN, MulangaC, BazepeoSE, MwambaJK, TshimpakaJ-W, KashiM, et al Distribution of HIV-1 variants in the Democratic Republic of Congo suggests increase of subtype C in Kinshasa between 1997 and 2002. J Acquir Immune Defic Syndr 1999. 2005 12 1;40(4):456–62., VidalN, MulangaC, BazepeoSE, MwambaJK, TshimpakaJ-W, KashiM, et al Distribution of HIV-1 variants in the Democratic Republic of Congo suggests increase of subtype C in Kinshasa between 1997 and 2002. J Acquir Immune Defic Syndr 1999. 2005 12 1;40(4):456–62.

  • (KouyosRD, von WylV, YerlyS, BöniJ, TafféP, ShahC, et al Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland. J Infect Dis. 2010 5 15;201(10):1488–97. doi: 10.1086/651951 20384495)

    KouyosRD, von WylV, YerlyS, B�niJ, Taff�P, ShahC, et al

    KouyosRD, von WylV, YerlyS, BöniJ, TafféP, ShahC, et al Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland. J Infect Dis. 2010 5 15;201(10):1488–97. doi: 10.1086/651951 20384495, KouyosRD, von WylV, YerlyS, BöniJ, TafféP, ShahC, et al Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland. J Infect Dis. 2010 5 15;201(10):1488–97. doi: 10.1086/651951 20384495

  • J van Harmelen (1997)

    81

    AIDS Lond Engl, 11

Publisher
Pubmed Central
Copyright
© 2017 Pérez-Parra et al
ISSN
1932-6203
eISSN
1932-6203
DOI
10.1371/journal.pone.0186928
Publisher site
See Article on Publisher Site

Abstract

OPENACCESS Phylogenetic studies are a valuable tool to understand viral transmission patterns and the Citation: PeÂrez-Parra S, Chueca N, Alvarez M, role of immigration in HIV-1 spread. We analyzed the spatio-temporal relationship of differ- Pasquau J, Omar M, Collado A, et al. (2017) High ent HIV-1 non-B subtype variants over time using phylogenetic analysis techniques. We col- prevalence and diversity of HIV-1 non-B genetic lected 693 pol (PR+RT) sequences that were sampled from 2005 to 2012 from naïve forms due to immigration in southern Spain: A patients in different hospitals in southern Spain. We used REGA v3.0 to classify them into phylogeographic approach. PLoS ONE 12(10): e0186928. https://doi.org/10.1371/journal. subtypes and recombinant forms, which were confirmed by phylogenetic analysis through pone.0186928 maximum likelihood (ML) using RAxML. For the main HIV-1 non-B variants, publicly avail- Editor: Jason Blackard, University of Cincinnati able, genetically similar sequences were sought using HIV-BLAST. The presence of HIV-1 College of Medicine, UNITED STATES lineages circulating in our study population was established using ML and Bayesian infer- Received: May 1, 2017 ence (BEAST v1.7.5) and transmission networks were identified. We detected 165 (23.4%) patients infected with HIV-1 non-B variants: 104 (63%) with recombinant viruses in pol: Accepted: October 10, 2017 CRF02_AG (71, 43%), CRF14_BG (8, 4.8%), CRF06_cpx (5, 3%) and nine other recombi- Published: October 30, 2017 nant forms (11, 6.7%) and unique recombinants (9, 5.5%). The rest (61, 37%) were infected Copyright:© 2017 Perez-Parra et al. This is an with non-recombinant subtypes: A1 (30, 18.2%), C (7, [4.2%]), D (3, [1.8%]), F1 (9, 5.5%) open access article distributed under the terms of and G (12, 7.3%). Most patients infected with HIV-1 non-B variants were men (63%, p< the Creative Commons Attribution License, which permits unrestricted use, distribution, and 0.001) aged over 35 (73.5%, p< 0.001), heterosexuals (92.2%, p< 0.001), from Africa reproduction in any medium, provided the original (59.5%, p< 0.001) and living in the El Ejido area (62.4%, p<0.001). We found lineages of author and source are credited. epidemiological relevance (mainly within Subtype A1), imported primarily through female Data Availability Statement: All relevant data are sex workers from East Europe. We detected 11 transmission clusters of HIV-1 non-B Sub- within the paper and its Supporting Information types, which included patients born in Spain in half of them. We present the phylogenetic files. profiles of the HIV-1 non-B variants detected in southern Spain, and explore their putative Funding: This work was supported by National geographical origins. Our data reveals a high HIV-1 genetic diversity likely due to the import Plan of I+D+i (RD12/0017/006), granted to FGG; of viral lineages that circulate in other countries. The highly immigrated El Ejido area acts as General Ministry of Evaluation and European, granted to FGG; Fund of Regional Development a gateway through which different subtypes are introduced into other regions, hence the (FEDER), granted to FGG; Health Institute Carlos III importance of setting up epidemiological control measures to prevent future outbreaks. (ISCIII), PI12/01053, granted to FGG; and Ministry of Health of Andalusia (AC-0082-2013), granted to PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 1 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain FGG. The funders had no role in study design, data Introduction collection and analysis, decision to publish, or The high evolutionary rate and recombination capacity of human immunodeficiency virus preparation of the manuscript. type 1 (HIV-1) determine the existence of an array of subtypes and recombinant forms circu- Competing interests: The authors have declared lating worldwide [1±4]. HIV-1 non-subtype B (ªnon-Bº) variants cause around 90% of infec- that no competing interests exist. tions worldwide, and largely predominate in African or Eastern European countries with generalized HIV-1 epidemics. Subtypes C and A, and circulating recombinant forms (CRF) CRF01_AE and CRF02_AG, are responsible alone for 70% of the world's infections [5]. Nowa- days the proportion of infections by HIV-1 non-B variants in Spain lies at 12±15%, depending on the study and technique used to characterize variants [6,7]. Nonetheless, the predominance of HIV-1 subtype B in developed countries (where antiretroviral therapy is more widespread), implies that this is the most widely studied subtype from the genetic, biological and therapeutic viewpoints. The full biological meaning of the genetic variability of HIV-1 is still not completely understood. However, several major differences between the biological properties of certain genetic subtypes in have been described; e.g., virulence, tropism and transmissibility [8,9], use of chemokine co-receptors [10], disease progression [11], susceptibility to some anti- retroviral drugs [12,13], sensitivity to viral load quantification methods [14,15] and detection [16]. These findings evidence the importance of epidemiological information about different subtypes. Eastern Andalusia is located in south-eastern Spain, and includes the provinces of Almerõ Âa, Granada and Jae Ân. Given its location and geographic closeness to the African continent, this region has received a notable foreign migratory influx in the last decade. Andalusia is the fourth Spanish Autonomous Community in number of foreign population, only surpassed by Catalonia, Madrid and the Valencian Community. The main source of immigration in Eastern Andalusia stems from its intensive farming practices, mainly in the El Ejido area (located in the province of Almerõ Âa), where one in every four citizens is an immigrant. Phylogenetic analyses, in conjunction with geographical data, can assess the existing rela- tionship between migratory events and spread of HIV-1 on a local scale [17±20], and to study HIV-1 transmission networks locally [21±24]. As in previous studies [25], our center collects the HIV-1 pol gene sequences linked to the patients' clinical data to monitor baseline drug resistance in naïve individuals in Eastern Andalusia. Our aims were to describe the molecular epidemiology and evolutionary history of non-B forms in Eastern Andalusia over the 2005± 2012 period, and to explore their putative geographical origin prior to their arrival to our region. Methods Study population During the study period (2005±2012), 693 pol gene sequences of patients newly diagnosed with HIV-1 in different Eastern Andalusian hospitals were collected from routine drug resis- tance analyses. These hospitals were distributed in 3 provinces: Granada (which included its    capital city of Granada and Motril), Jaen, and Almerõa (including its capital city of Almerõa and El Ejido). The pol sequences (protease (PR), codons 4±99; reverse transcriptase (RT), codons 38±247) obtained by the Trugene1 HIV Genotyping kit (Siemens, NAD), were linked to demographic (risk group, age, gender, country of origin, sampling year, and attending hos- pital), clinical (CD4+ T-cell count) and virological (plasma viral load) information. Demo- graphic information was voluntarily collected during clinical interviews. This study was approved by the San Cecilio Hospital's Ethics Committee, and no consent information was required as patient information remained anonymous and was de-identified prior to analyses. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 2 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain HIV-1 pol sequencing and subtype assignment All the sequences were trimmed to 883 nucleotides (nt) and aligned using ClustalW [26]. The viral subtype was studied with the REGA v3.0 subtyping tool (http://dbpartners.stanford. edu:8080/RegaSubtyping/stanford-hiv/typingtool/), and was confirmed by phylogenetic analy- sis through maximum likelihood (ML) using the randomized Accelerated Maximum Likeli- hood (RAxML) program, accessible on the CIPRES Science Gateway [27]. The general time- reversible (GTR) model with a gamma-distributed heterogeneity rate across sites was employed, applying 1000 bootstrap iterations. A representative dataset of HIV-1 group M sequences, including non- recombinant subtypes (A-K) and recombinant forms (at least four representative sequences of each non-recombinant subtype and the CRF currently available from the analysis) were downloaded from the Los Alamos HIV sequence database (http:// www.hiv.lanl.gov) was used as a reference dataset (S1 Table). The assignment to any subtype/CRF was considered definitive if the query sequence was included with the reference sequences corresponding to that viral variant in a monophyletic cluster supported by high bootstrap values (>70%) [28]. Any genetic form not associated with reference subtypes/CRFs was classified as a unique recombinant form (URF), whose recombi- nation pattern was further studied by a Bootscan analysis using the SimPlot v3.5.1 software [29]. The bootscanning method in SimPlot consists of a sliding-window phylogenetic boot- strap analysis of the query sequence aligned against a set of reference strains to reveal break- points. The Neighbor-Joining algorithm was selected, with the Kimura 2-parameter substitution model. We employed a window size of 200nt moving in 10nt increments. We used a minimum cutoff for the bootstrap value of 70% to reliably assign each of the breakpoint segments to a parental variant. We have submitted to GenBank the major groups of HIV-1 non-B variants under accession numbers MF628109 to MF628250. These were defined as those found in at least five patients. With the aim of protecting the identity of patients infected with rare genetic forms of HIV-1, and for similar scientific and ethical reasons as explained in other HIV cohorts [30±32], we decided not to submit to GenBank those sequences corresponding to the less frequent variants. Inference of the putative geographical origins of the HIV-1 non-B variants circulating in Andalusia To further characterize the relationships among the major groups of HIV-1 non-B variants, we interrogated GenBank for genetically related sequences to our major subtypes/recombinant forms using HIV-BLAST (http://www.hiv.lanl.gov/content/sequence/BASIC_BLAST/basic_ blast.html). The 10 most closely related GenBank sequences to each of our study sequences, were downloaded and included in each dataset. We also included all the pol sequences (start: 2293 and end: 3290, HXB2 coordinates), available in the HIV Los Alamos database sampled in Spain for each dataset: subtype A1 (n = 60), subtype C (n = 52), subtype F (n = 143), subtype G (n = 64), CRF14_BG (n = 25), and CRF02_AG (n = 265). Since very few sequences for CRF06_cpx were available in public databases (http://www.hiv.lanl.gov/content/sequence/ HIV/mainpage.html), we included them all (n = 110). All these individual sequence datasets were put together (n = 970) and a global phylogenetic analysis was performed using RAxML (GTR + Gamma model) and 1000 bootstrap iterations for this analysis. The phylogenetic relatedness between the sequences was studied, and a 70% bootstrap value was taken as a significantly reliable value [28]. Thresholds for low genetic dis- tance, which are commonly used as a proxy for divergence time, were not applied to the cluster definition in the ML trees since these clusters were further confirmed and analyzed using a PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 3 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain time-stamped Bayesian phylogenetic analysis with BEAST, as described below. International non-B lineages (defined as phylogenetic associations of at least one sequence from our cohort clustered with sequences from different countries), and `Andalusian clusters' (monophyletic associations of sequences in our cohort alone), were identified in the global ML tree. A Bayesian Markov Chain Monte Carlo (MCMC) approach was applied to each of the indi- vidual HIV-1 non-B subtype/CRFs datasets described above, which included the most geneti- cally similar sequences found with HIV-1 BLAST, as implemented in BEAST v1.7.5 [33]. The Shapiro-Rambaut-Drummond-2006 (SRD06) substitution model was used, together with a relaxed uncorrelated lognormal clock (UCLN)[34] and a demographic non parametric model, Bayesian Skyline Plot (BSP) [35]. This model combination was chosen because it best fits the analysis of the HIV-1 pol data run in the majority of studies [36]. The MCMC was run for 250 million states sampling every 50000. The evolutionary rate (μ, nucleotide substitutions per site per year, subst./site/year) for the different HIV-1 non-B subtypes/CRFs (S2 Table), and the most recent common ancestors (MRCA) of the different HIV-1 non-B clusters, were esti- mated. Only traces with an effective sample size (ESS) > 200 for all the parameters, after excluding an initial 10% burnin, were accepted as visualised in TRACER, v1.6 (http://tree.bio. ed.ac.uk/software/tracer/). Maximum Clades Credibility (MCC) trees were constructed in each case to summarise the posterior tree distributions. In these MCC trees, the more epidemiologically relevant clusters and lineages, previously identified in the global ML tree, were studied; and a node support cut- off (posterior probability (pp) above 0.9) was applied for their confirmation. Trees were viewed and edited in FigTree, v. 1.4.0 (http://tree.bio.ed.ac.uk/software/figtree). Analysis of the antiretroviral drug resistance mutations Drug resistance mutations were identified in the pol sequences using the HIVseq program, which is available in the HIV Drug Resistance Database of Stanford University (https://hivdb. stanford.edu/hivseq/by-sequences/), and also using the WHO surveillance drug resistance mutation list (last updated in 2009 by Bennett and colleagues) [37]. Statistical analyses A multivariate logistic regression analysis was performed to determine the predictive effect of the demographic, clinical and virological characteristics on the adscription to each subtype/ CRF. The statistical significance of these characteristics, compared to the total proportion of infected patients, was studied by a hypothesis contrast using a z-test. The statistical analysis was performed with SPSS 22.0. Results Epidemiological surveillance of the non-B HIV-1 genetics forms Of the 693 total included patients, 165 (23.8%) were infected with different genetic forms of HIV-1 non-B variants. Most of them (n = 104, 63%) were recombinant viruses in pol: 95 (57.6%) corresponded to 12 different CRFs and nine (5.5%) were URFs. The other patients (n = 61, 37%) were infected with five non-recombinant subtypes: A1, D, C, F1 and G (see Fig 1). The demographic, clinical and virological characteristics of the patients according to the genetic HIV-1 non-B forms are provided in Table 1. Most of the patients infected with non-B variants were men (63%, p < 0.001) aged over 35 (73.5%, p < 0.001), heterosexual (92.2%, p < 0.001), African (58.2%, p < 0.001), and living in the El Ejido area (62.4%, p<0.001). The PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 4 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 1. Distribution of the HIV-1 non-subtype B genetic forms detected in Eastern Andalusia over the 2005±2012 period. https://doi.org/10.1371/journal.pone.0186928.g001 full list of countries of origin for patients infected with non-B forms and born abroad (n = 127 [77%]) were: Argentina, n = 1, Brazil, n = 5, Burkina Faso, n = 1, Cameroon, n = 1, Colombia, n = 1, Congo, n = 5, Ivory Coast, n = 1, Cuba, n = 1, Gambia, n = 2, Ghana, n = 13, Guinea, n = 10, Guinea-Bissau, n = 11, Equatorial Guinea, n = 4, Lithuania, n = 1, Mali, n = 9, Morocco, n = 3, Mauritania, n = 1, Nigeria, n = 17, Dominican Republic, n = 1, Romania, n = 6, Russia, n = 14, Senegal, n = 15, Sierra Leone, n = 1 and South Africa, n = 3. The rest of subjects (n = 36, 22%) had been born in Spain. The multivariate logistic regression analyses demonstrated a higher risk of carrying HIV-1 subtype A for females (OR = 6.17, p = 0.026) and non Africans (OR = 0.08, p = 0.008; S3 Table). The other HIV-1 non-B genetic forms showed no predictive effect of the demographic, clinical and virological characteristics (data not shown). Twenty-three patients were infected with unusual HIV-1 non-B variants (i.e., those variants found in four patients or fewer). Of them, 10 (44.4%) were observed in Spanish patients (Table 2). The recombination patterns for the different URFs obtained according to the Boot- scan analysis are presented in Fig 2. Geographical distribution of the various HIV-1 non B genetic forms The geographic distribution of the different HIV-1 non B subtypes and recombinant forms are represented on the map of Eastern Andalusia (Fig 3). Most of the patients infected with HIV-1 non-B variants were sampled in El Ejido (62.4%) or in the city of Granada, (22.4%), whereas   non-B variants were less frequent in the cities of Almerõa (10.9%), Jaen (2.4%) and Motril (1.8%). PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 5 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Analysis of the putative geographical origins of the main HIV-1 non-B genetic forms found in Eastern Andalusia In order to characterize the phylogenetic relationship of the patients infected with the most frequently found HIV-1 non-B variants (those found in 5 patients), the global ML tree (Fig 4) revealed the existence of 13 international lineages in Eastern Andalusia (Table 3) and 11 Andalusian clusters (Table 4) that involved patients in our cohort. The Bayesian analyses (Figs 5 and 6) showed that most of these Andalusian clusters origi- nated in the first decade of this century, and mainly included patients sampled in El Ejido. The low CD4 count of the patients included in most of these transmission networks suggests a late HIV diagnosis in a high proportion of patients (Table 4). Table 1. Demographic, clinical and virological characteristics of the patients infected with HIV-1 non B variants sampled over the 2005±2012 period. NON-RECOMBINANT SUBTYPES RECOMBINANT FORMS CHARACTERISTICS A1 C F1 G CRF02_AG CRF14_BG CRF06_cpx Others Total Gender Male 10 (33.3) 4 (57.1) 8 (88.9) 9 (75) 48 (69) 6 (75) 3 (60) 15 (65.2) 103 (62.4) Female 20 (66.7) 3 (42.9) 1 (11.1) 3 (25) 23 (31) 2 (25) 2 (40) 8 (34.8) 62 (37.6) Age (159)* <35 12 (40) 3 (42.9) 4 (50) 3 (25) 16 (23.2) 1 (20) 3 (15) 42 (26.4) 35±45 12 (40) 3 (42.9) 2 (25) 4 (33.3) 31 (44.9) 5 (62.5) 4 (80) 13 (65) 74 (46.5) >45 6 (20) 1 (14.3) 2 (25) 5 (41.7) 22 (31.9) 3 (37.5) 4 (20) 43 (27) Risk factor (153)* HTX 30 (100) 7 (100) 8 (88.9) 12 (100) 61 (95.3) 6 (75) 3 (75) 14 (73.7) 141 (92.2) MSM 1 (1.6) 1 (12.5) 1 (25) 5 (26.3) 8 (5.2) IVDU 1 (11.1) 2 (3.1) 1 (12.5) 4 (2.6) Country of origin (163)* Spanish 11 (36.7) 4 (44.4) 1 (8.3) 6 (8.6) 4 (50) 10 (47.6) 36 (22) East Europe 15 (50) 2 (28.6) 2 (22.2) 2 (9.5) 21 (12.9) North Africa 3 (4.2) 3 (1.8) West Africa 2 (6.7) 2 (28.6) 1 (11.1) 11 (91.7) 54 (76) 4 (50) 5 (100) 3 (14.3) 82 (50.3) Central and South Africa 1 (3.3) 1(14.3) 8 (11.3) 2 (9.5) 12 (7.4) Central and South America 1 (3.3) 2 (28.6) 2 (22.2) 4 (19) 9 (5.5) Viral Load (158)* <10000 8 (26.7) 1 (14.3) 1 (12.5) 1 (8.3) 15 (22.1) 1 (12.5) 5 (23.8) 32 (20.3) 10000±100000 10 (33.3) 2 (28.6) 1 (12.5) 6 (50) 24 (35.3) 5 (62.5) 2 (50) 9 (42.9) 59 (37.3) >100000 12 (40) 4 (57.1) 6 (75) 5 (41.7) 29 (42.6) 2 (25) 2 (50) 7 (33.3) 67 (42.4) CD4 count (149)* <200 4 (16) 4 (50) 6 (50) 35 (52.2) 2 (25) 1 (20) 1 (5.3) 53 (35.6) 201±350 6 (24) 3 (60) 3 (37.5) 3 (25) 19 (28.4) 3 (37.5) 2 (40) 4 (21.1) 43 (28.9) >350 15 (60) 2 (40) 1 (12.5) 3 (25) 13 (19.4) 3 (37.5) 2 (40) 14 (73.7) 53 (35.6) Sampling date interval 2005±2007 4 (13.3) 3 (42.9) 6 (50) 11 (15.5) 1 (12.5) 3 (13) 28 (17) 2007±2009 3 (10) 1 (14.3) 1 (11.1) 2 (16.7) 30 (42.3) 5 (62.5) 4 (80) 3 (13) 49 (29.7) 2009±2011 10 (33.3) 2 (28.6) 1 (11.1) 3 (25) 12 (16.9) 1 (12.5) 1 (20) 8 (34.8) 38 (23) 2011±2013 13 (43.3) 1 (14.3) 7 (77.8) 1 (8.3) 18 (25.4) 1 (12.5) 9 (39.1) 50 (30.3) Total 30 7 9 12 71 8 5 23 165 *Date available for the number of indicated patients. https://doi.org/10.1371/journal.pone.0186928.t001 PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 6 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Table 2. Clinical, demographic and virological characteristics of the patients infected with infrequent HIV-1 non B genetic variants over the 2005± 2012 period. Patient ID Genetic form Gender (M/F) Risk Factor Age Year of diagnosis CD4 count Viral Load Country of origin 1 CRF01_AE M HTX 38 2008 764 421 Spain 2 CRF03_AB M MSM 28 2011 715 404000 Spain 3 CRF20_BG M HTX 65 2005 521 5910 Spain 4 CRF20_BG F HTX 40 2005 476 77700 Cuba 5 CRF22_01_A1 M HTX 37 2009 474 11700 Equatorial Guinea 6 CRF25_CPX F UNK NA 2012 NA 36300 Spain 7 CRF26_AU F HTX 25 2012 15 1180000 Spain 8 CRF37_CPX M HTX 39 2007 405 135000 Equatorial Guinea 9 CRF47_BF M HTX 39 2009 389 12279 Spain 10 CRF47_BF F HTX 32 2011 392 4122 Spain 11 CRF49_CPX M HTX 36 2009 288 2250000 Senegal 12 D F HTX 39 2011 443 22233 Equatorial Guinea 13 D M MSM 48 2009 526 55723 Spain 14 D M MSM 41 2009 416 417477 Colombia 15 URF BC F HTX 35 2011 723 27911 Romania 16 URF BC M HTX 40 2011 450 211933 Romania 17 URF BF M MSM 38 2011 246 143039 Spain 18 URF BF F HTX 38 2012 NA 44100 Brazil 19 URF BF M HTX 52 2006 462 3200 Argentina 20 URF BG F NA NA 2008 NA NA NA 21 URF BG M NA NA 2009 NA NA NA 22 URF BG M HTX 38 2009 NA NA Guinea-Bissau 23 URF BD M NA 49 2009 259 2800 Spain NA: Not Available; M: Male; F: Female. https://doi.org/10.1371/journal.pone.0186928.t002 In order to provide more information about the scale of the trees shown, we provide in the S4 Table the distribution of patristic (uncorrected) pairwise genetic distances between sequences included in each of the ML and Bayesian trees generated in this article. Non-recombinant subtypes Thirty (18.2%) patients were infected with HIV-1 subtype A1. The viral sequences were geneti- cally similar according to HIV-BLAST to 21 GenBank sequences from Bulgaria, the Demo- cratic Republic of Congo, Croatia and Greece with 13, 6, 1 and 1 cases, respectively. The ML analysis (Fig 4) detected a large international lineage that involved sequences from Eastern Europe (lineage L1.A1 in Table 3) and grouped 21 patients from our cohort: 16 women born abroad (Eastern Europe (n = 14), the Dominican Republic (n = 1) and Lithuania (n = 1)) and 5 Spanish men. This lineage also included 23 GenBank sequences, also originating from Eastern Europe: Bulgaria, n = 10, Russia, n = 5, Poland, n = 1 and the Ukraine, n = 1. Within this line- age, we found two clusters (A.1 and A.2), formed exclusively by Spanish men and female sex workers born in Russia, all being patients sampled in Eastern Andalusia. The A.1 local cluster involved 4 sequences from Spanish patients living in the capital of Granada, its origin was esti- mated to be 2008.5 (95%CI: 2006.6±2010.3), and the sequences presented the resistance muta- tion K103N in the RT gene. This cluster was also phylogenetically related to viruses that circulate in Eastern Europe. Unlike most of the HIV-1 non-B clusters, patients in the A.1 clus- ter showed a high CD4 count (mean = 590, range = 534±701). Moreover, the Bayesian PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 7 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 2. Bootscan analysis of the unique recombinant forms (URF) found in Eastern Andalusia. The analysis was applied to the concatenated sequences that corresponded to HXB2 coordinates 2283±2549 (PR) and 2661±3290 (RT). https://doi.org/10.1371/journal.pone.0186928.g002 phylogenetic tree revealed short internode branches, which may indicate short times between infections. Finally, five subtype A1 sequences from our cohort corresponded to patients from Africa: Mali, n = 2, Equatorial Guinea, n = 1 and Spain, n = 2, not clustered in transmission cluster. Seven (4.2%) sequences corresponded to HIV-1 subtype C, and showed high genetic simi- larity to 23 GenBank sequences sampled in South Africa (n = 12), Brazil (n = 6) and Bulgaria (n = 4). We thus found two main ways of subtype C entrance to our area: South Africa and Brazil: a Brazilian male patient from our cohort grouped with 6 GenBank sequences from Bra- zil; and a South African male patient grouped with GenBank sequences from South Africa (n = 2) and Somalia (n = 1) (Fig 4). Within this subtype, we also found a single cluster (C.1) formed by patients from Brazil (n = 1) and Romania (n = 2). Nine (5.5%) sequences corresponded to HIV-1 subtype F1 and showed a high genetic simi- larity to 21 GenBank sequences from Brazil (n = 14), Bulgaria (n = 4) and the Democratic Republic of Congo (n = 3). We found only one Andalusian F1 cluster: a sequence pair (cluster F.1), that originated in 2010.2 (95%CI: 2010±2011) and was formed by two male injection drug users sampled in Jaen and who were of Brazilian and Spanish origins. This sequence pair was included among GenBank sequences from Brazil in the ML tree. However, we found 2 international lineages: L1.F, which grouped two Romanian heterosexual patients from our cohort with GenBank sequences sampled in Eastern Europe, mainly Romania (n = 6) and Bul- garia (n = 2). The second F1 subtype lineage (L3.F) included Spanish men who have sex with men (MSM) sampled in North Spain, and also a MSM from our cohort. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 8 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 3. Geographical distribution of the patients infected with HIV-1 non-B variants over the 2005±2012 period. The percentage of each subtype/CRF in relation to all the HIV-1 non-B genetic forms is shown in each region. https://doi.org/10.1371/journal.pone.0186928.g003 Twelve (7.3%) patients of our cohort were infected with HIV-1 subtype G, who came from different western and central African countries: Mali (n = 1), Nigeria (n = 6), Ghana (n = 3) and Guinea-Bissau (n = 1). They presented high genetic similarity to 5 GenBank sequences from the Republic of Congo (n = 4) and Bulgaria (n = 1). None of these sequences was epide- miologically related according to our data. We found only one Nigerian patient whose sequence grouped with another one of the same country of origin (lineage L1.G). HIV-1 recombinant forms Eight (4.8%) patients in our cohort were infected with the recombinant CRF14_BG form, which in the pol analyses typically forms a monophyletic cluster within the subtype G crown. These eight patients came from Spain (n = 4), Guinea (n = 2), and Guinea-Bissau (n = 2). We also found a sin- gle small Andalusian cluster (cluster 14BG.1), which originated in 2004.4 (95% CI:2003.8±2005), and was formed by two Spanish patients. Finally, two patients from Guinea and Guinea Bissau grouped with sequences from Equatorial Guinea (lineages L1.14BG and L2.14BG). We found 71 (43%) patients, mainly from western African countries (77.5%), infected with CRF02_AG. Of these, 11 (14%) were grouped into five small Andalusian Cluster: 4 clusters of two patients and one with three patients. We detected 5 different lineages (L1.02AG-L5.02AG) of viruses sampled in other countries, with patients from our cohort who came mainly from Western Africa. To study the phylogenetic profile of variant CRF06_cpx, we used all the sequences available in Los Alamos HIV given their small number, n = 110 (see Fig 4). We found 5 patients in our cohort (3%) to be infected with variant CRF06_cpx, who came from different western African countries: Nigeria (n = 3), Ghana (n = 1) and Senegal (n = 1). These sequences grouped with GenBank sequences from the neighboring Western African countries of Burkina Faso, Togo and Nigeria. However, we found no significant association among the patients infected with PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 9 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 4. Global ML phylogenetic tree inferred for the main HIV-1 non-B genetic forms sampled in Eastern Andalusia. The phylogenetic tree was constructed by the general time-reversible with gamma-distributed rate heterogeneity across sites model of substitution implemented into RAxML. Branches are drawn on scale with the bar at the bottom, which represents 0.04 nucleotide substitution per site. Statistically highly supported nodes (bootstrap values>70%) are indicated by an asterisk (*). Andalusian clusters and international lineages are highlighted in yellow and blue, respectively. The Andalusian sequence names contain a three-part code: Sequence number, sampling site (AL: Almerõ Âa, EJ: El Ejido GR: Granada, JA: Jae  n, MO: Motril) and the code of the most likely country of infection. https://doi.org/10.1371/journal.pone.0186928.g004 this genetic form, and the CRF06_cpx sequences sampled in our cohort were interspersed in the tree. Discussion In Eastern Andalusia, most HIV-1 non-B subtype genetic forms were found among immigrant heterosexual population, mainly African males or Eastern European females. These patients were living preferentially in El Ejido, an area that potentially acts as a gateway for diverse HIV- PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 10 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Table 3. HIV-1 non-B international lineages involving sequences sampled in Eastern Andalusia and sequences from different countries. Lineages Sequences Country of origin of the Node Country of Most likely countries of infection for Articles with (EA;OR) patients sampled in our support collection for other other patients (Accession ID) related cohort (n) (BT) sequences (n) sequences Subtype A L1.A 21;23 Spanish (5), Russian 97 Spain (11) Ukraine (KC340416) [20,38,39] (14), Lithuanian (1) Bulgaria (10) Russia (KC340417, KC340407, Dominican (1) Poland (1) KC340486, KC340436, KC340495) Croatia (1) Subtype C L1.C 1;1 Senegalese (1) 100 Spain (1) Morocco (JQ351995) [40] Subtype F L1.F 2:21 Romanian (2) 70 Spain (19) Romania (KC340133, KC340134, [20,39,41] Bulgaria (2) KC340419, KC340378, KC340643, KC340132) L2.F 1;1 Spanish (1) 100 Spain (1) Sub-Saharan Africa (KC340440) L3.F 1;112 Spanish (1) 100 Spain (112) Spain (KJ883085-KJ883089, [42,43] KJ883091-KJ883108, KJ883110-KJ883152, JN010216-JN859590) Subtype G L1.G 1;1 Nigerian (1) 100 Spain (1) Nigeria (FJ481667) [44] CRF_14BG L1.14BG 1;1 Guinean-Bissau (1) 100 Spain (1) Equatorial-Guinea (JX428555) [45] L2.14BG 1;1 Guinean (1) 100 Spain (1) Equatorial-Guinea (EU255306) [45] Form CRF02_AG L1.02AG 1;2 Guinean (1) 100 Spain (2) Ghana (KC340379) [23] L2.02AG 2;1 Moroccan (1) 86 Spain (1) Cameroon (KC340644) Spanish (1) L3.02AG 1;1 Nigerian (1) 86 Spain (1) Spain (HF567877) [46] L4.02AG 1;6 Ghanaian (1) 96 Spain (6) Bolivia (EU255444) [44,47] Ecuador (KC340123, KC340124, FJ481711) L5.02AG 1;5 Equatorial-Guinean (1) 82 Spain (5) Equatorial-Guinea (EU255527) [47,48] EA: sequences sampled in Eastern Andalusia; OR: sequences sampled in other geographic regions; BT: bootstrap (in the global maximum-likelihood tree). https://doi.org/10.1371/journal.pone.0186928.t003 1 variants to enter the Eastern Andalusian region. These findings are explained by the fact the El Ejido's economy is mainly based on greenhouse farming, for which a large industry has emerged in recent years thanks to immigrant labor, made up of people mainly from Africa. The prevalence of HIV-1 non-B variants in eastern Andalusia is similar to that reported in a study performed in the nearby Western areas of Andalusia (23%) [49], but is still much higher than that found elsewhere in Spain [6,7]. An increased prevalence has been noted for HIV-1 non-B variants and their genetic diversity in Eastern Andalusia in recent years: 22% of autochthonous patients were infected with HIV-1 non-B forms between 2005 and 2012, as opposed to the 12.8% reported in former studies conducted between 1997 and 2001 [50]. We also detected 12 different CRFs and nine URFs, a variability that is probably related to the increased migration rate reported in southern Spain in the last decade [49,51]. The least frequent HIV-1 non-B variants were detected often among Spanish patients (43%, [10/23]), and most of the clusters formed by these variants included at least one Spanish patient (55%, [6/11]). These data suggest that although these HIV-1 non-B variants seem to be PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 11 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Table 4. Demographical, clinical, virological and phylogenetic characteristics of the Andalusian clusters found for the main HIV-1 non B variants. Cluster No. Sampling Support Risk Country of Viral Load CD4 count Location area tMRCA (95% ID Patients interval (BT:PP) factor origin (median, IQR) (median, HPD) Log IQR) Subtype A A.1 4 2011 100;0.9 HTX Spanish 5 (4.6±6) 590 (534± Granada 2008.5 (2006.6± 701) 2010.3) A.2 3 2007±2011 96;0.9 HTX Spanish- 4.7 (4.4±5.3) 301 (152± Granada-El 2003.9 (1998.4± Russian 517) Ejido 2004.7) A.3 2 2009±2010 100;0.9 HTX Spanish- 4.9 (4±5) 391 (325± El Ejido 2008 (2007.4± Russian 456) 2009) Subtype C C.1 3 2005±2012 100;0.9 HTX Romanian- 5.2 (5.1±5.3) 499 (253± Granada-Ejido 1998.3 (1992.3± Brazilian 739) 2001.9) Subtype F F. 2 2011±2012 90;0.9 HTX Brazilian 5.7(5.4±5.9) 62 (50±73) Jaen 2010.2 (2010± 2011) Form CRF14_BG 14BG.1 2 2005±2007 100;0.9 HTX Spanish 5 425 (267± Granada 2004.4 (2003.8± 583) 2005) Form CRF02AG 02AG.1 2 2007 100;0.9 HTX Guinean 5.6(4±6) 309 (180± El Ejido 2003.8 (2001± 438) 2006) 02AG.2 2 2007±2010 100;0.9 HTX Senegalese 3(4±5) 391 (381± El Ejido 2006.4 (2006± 400) 2007) 02AG.3 2 2011 99;0.9 HTX Spanish- 5.9(3±6) 326 (100± Granada- 2008.2 (2007.2± Moroccan 551) Almerõ Âa 2009.9) 02AG.4 3 2006±2011 98;0.9 HTX Spanish- Mali- 5.5(5.3±6.2) 140 (78± El Ejido 2003.6 (2000.3± 306) 2005.3) 02AG.5 2 2008 100;0.9 HTX Ghana 6 126 (45± El Ejido 2007.7 (2007.2± 200) 2009) https://doi.org/10.1371/journal.pone.0186928.t004 due to imported cases in most cases, they have also gradually penetrated the autochthonous population in recent years. The phylogenetic and epidemiological study of the HIV-1 non-B variants in our region showed that these variants account for high proportion of infections among migrant patients, and that these viruses were genetically close to those circulating in these subjects' countries of origin. This indicates that many patients were infected before they arrived in Spain. These sequences sampled in other countries, and available in public databases, act as a control to avoid overestimating the local transmission clusters that include patients who are most proba- bly unrelated in epidemiological terms. As previously shown in a national study [7], CRF02_AG was the most frequent HIV non-B variant in our population (43%). Nonetheless, the small proportion of their phylogenetic asso- ciation is surprising (14%, [11/71]). This clustering rate was much higher for other HIV-1 non-B subtypes, such as subtype A1 (27.3%, [9/33]), where we discovered an international lineage (L1.A1) that mostly included a particularly vulnerable group of Russian female sex workers and potentially their local customers. According to our analysis, it would appear that most of the non-B cases detected in Eastern Andalusia were generally imported cases as most were identified in immigrant populations. Our analyses suggest that many of these cases form part of international HIV-1 lineages that originated in Eastern Europe, South America and sub-Saharan Africa. However, we also PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 12 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 5. Bayesian phylogenetic tree inferred for the subtype A1, C, F1 and G/CRF14_BG pol sequences sampled in Eastern Andalusia and genetically similar sequences from GenBank. Red branches correspond to the sequences sampled in eastern Andalusia from 2005 to 2012. Statistically highly supported nodes (posterior probability values above 0.9) are indicated with an asterisk (*). Andalusian clusters are highlighted in yellow. Andalusian sequences names contain a three-part code: Sequence number, sampling site (AL: Almerõ Âa, EJ: El Ejido GR: Granada, JA: Jae  n, MO: Motril) and the code of the most likely country of infection. https://doi.org/10.1371/journal.pone.0186928.g005 PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 13 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain Fig 6. Bayesian phylogenetic tree inferred for the CRF02_AG pol sequences sampled in Eastern Andalusia and genetically similar sequences from GenBank. Red branches correspond to the sequences sampled in eastern Andalusia from 2005 to 2012. Statistically highly supported nodes (posterior probability values above 0.9) are indicated with an asterisk (*). Andalusian clusters are highlighted in yellow. Andalusian sequences names contain a three-part code: Sequence number, sampling site (AL: Almerõ Âa, EJ: El Ejido GR: Granada, JA: Jae  n, MO: Motril) and the code of the most likely country of infection. https://doi.org/10.1371/journal.pone.0186928.g006 PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 14 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain identified 11 intra-region clusters, which might suggest the local dissemination of some non-B variants, particularly those which involve autochthonous Spanish subjects (6/11) and recent emergence times according to the phylogenetic reconstruction. On the other hand, clusters formed by foreign subjects with old common ancestors most likely reflect imported infections. The methods used herein involve a number of sampling limitations that affect this and many other similar studies. Since we relied on a BLAST search to identify the genetically clos- est sequences (from both Spain and abroad) that could form part of the same transmission net- works as our sequences, we depended on the sequences deposited in databases. Unfortunately, this availability is sometimes very low, particularly for non-B variants. Therefore, we cannot rule out that close and more informative sequences were not captured as they have not been sampled. This was the reason why we added all the sequences available in HIV Los Alamos col- lected in Spain. We also demonstrated the presence of 13 different lineages of viruses that cir- culated in our region, which grouped with other patients from different Spanish cohorts, mainly foreign patients. Fortunately, very few sequences included in transmission clusters persented resistance to first-line antiretroviral drugs. This information agrees with the common conception that viruses with resistance mutations present a biological disadvantage against wild strains, which weakens their transmission efficacy. Likewise, the drug resistance mutations detected affect mainly reverse transcriptase inhibitor drugs. We detected transmitted resistance mutations in four of the five patients grouped in Cluster A.4, which would cause high level resistance to nevirapine and efavirenz. We decided to study only the resistance mutations present in trans- mission clusters, which would have a stronger epidemiological impact. Further detailed infor- mation will be provided in future works. The constant epidemiological surveillance in our population, for which phylogenetic analysis tools are used, is a particularly important measure to study past outbreaks of genetic HIV-1 non- B subtype variants, and to prevent future ones. Likewise, as transmission cluster size seems to predict its expansion in time [52], we could expect some transmission chains of HIV-1 non-sub- type variants to become larger in size in forthcoming years, and more Spanish individuals to be included. We herein detected the presence of one patient from our cohort related to a fast spread- ing cluster among Spanish MSM infected with subtype F in Galicia (NW Spain) [43], a transmis- sion cluster which, as Delgado et al. suggest, would probably be closely linked to viruses that circulate in Eastern Europe [42]. These authors [53] have also described a subtype A cluster that is being transmitted among individuals in different areas of Spain. Finally, Patiño et al. [54] have warned about the novel appearance of variant CRF19_cpx among Spanish MSM individuals. Adequate knowledge about the characteristics of local epidemics, the study of risk groups and the prevalence of different viral subtypes are all fundamental aspects to successfully design HIV-1 prevention campaigns. In the present study, we demonstrate that phylogenetic studies which combine demographic, clinical and geographical data from different HIV-1 non-B subtypes in Eastern Andalusia provide very useful information to epidemiologically monitor and control HIV-1 spread and its origin in imported cases. Its use will help to reinforce and implement effi- cient actions to prevent HIV-1 from spreading between autochthonous and migrant populations. Supporting information S1 Table. HIV-1 reference sequence dataset used in the phylogenetic analysis. (DOCX) S2 Table. Evolutionary rates for each of the main HIV-1 lineages found in this study obtained through Bayesian phylogenetic inference. (DOCX) PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 15 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain S3 Table. Multivariate logistic regression analysis performed for the HIV-1 subtype A1 infections. (DOCX) S4 Table. Distribution of patristic (uncorrected) genetic distances among the HIV-1 pol sequences in each dataset included in this study. (DOCX) Author Contributions Conceptualization: Santiago Pe Ârez-Parra, Natalia Chueca, Marta Alvarez, Gonzalo Yebra. Data curation: Santiago Pe Ârez-Parra, Juan Pasquau, Mohamed Omar, Antonio Collado, David Vinuesa, Ana Belen Lozano, Federico Garcõ Âa. Formal analysis: Santiago Pe Ârez-Parra, Federico Garcõ Âa. Funding acquisition: Federico Garcõ Âa. Investigation: Santiago Pe Ârez-Parra, Natalia Chueca, Marta Alvarez, Federico Garcõ Âa. Methodology: Santiago Pe Ârez-Parra, Natalia Chueca, Gonzalo Yebra, Federico Garcõ Âa. Project administration: Federico Garcõ Âa. Resources: Federico Garcõ Âa. Software: Santiago Pe Ârez-Parra. Supervision: Santiago Pe Ârez-Parra, Juan Pasquau, Mohamed Omar, Antonio Collado, David Vinuesa, Ana Belen Lozano, Gonzalo Yebra, Federico Garcõ Âa. Validation: Santiago Pe Ârez-Parra, Federico Garcõ Âa. Visualization: Santiago Pe Ârez-Parra, Federico Garcõ Âa. Writing ± original draft: Santiago Pe Ârez-Parra, Gonzalo Yebra, Federico Garcõ Âa. Writing ± review & editing: Santiago Pe Ârez-Parra, Gonzalo Yebra, Federico Garcõ Âa. References 1. Vidal N, Mulanga C, Bazepeo SE, Mwamba JK, Tshimpaka J-W, Kashi M, et al. Distribution of HIV-1 variants in the Democratic Republic of Congo suggests increase of subtype C in Kinshasa between 1997 and 2002. J Acquir Immune Defic Syndr 1999. 2005 Dec 1; 40(4):456±62. 2. Taniguchi Y, Takehisa J, Bikandou B, Mboudjeka I, N'Doundou-N'Kodia M-Y, Obengui null, et al. Genetic subtypes of HIV type 1 based on the vpu/env sequences in the Republic of Congo. AIDS Res Hum Retroviruses. 2002 Jan 1; 18(1):79±83. https://doi.org/10.1089/088922202753394745 PMID: 3. Shankarappa R, Margolick JB, Gange SJ, Rodrigo AG, Upchurch D, Farzadegan H, et al. Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol. 1999 Dec; 73(12):10489±502. PMID: 10559367 4. Robertson DL, Anderson JP, Bradac JA, Carr JK, Foley B, Funkhouser RK, et al. HIV-1 nomenclature proposal. Science. 2000 Apr 7; 288(5463):55±6. PMID: 10766634 5. Hemelaar J, Gouws E, Ghys PD, Osmanov S, WHO-UNAIDS Network for HIV Isolation and Characteri- sation. Global trends in molecular epidemiology of HIV-1 during 2000±2007. AIDS Lond Engl. 2011 Mar 13; 25(5):679±89. 6. Yebra G, de Mulder M, Martõ Ân L, Rodrõ Âguez C, Labarga P, Viciana I, et al. Most HIV type 1 non-B infec- tions in the Spanish cohort of antiretroviral treatment-naïve HIV-infected patients (CoRIS) are due to recombinant viruses. J Clin Microbiol. 2012 Feb; 50(2):407±13. https://doi.org/10.1128/JCM.05798-11 PMID: 22162552 PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 16 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain 7. Garcõ Âa F, Pe  rez-Cachafeiro S, Guillot V, Alvarez M, Pe  rez-Romero P, Pe  rez-Elõ Âas MJ, et al. Transmis- sion of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treat- ment naïve HIV-infected individuals (CoRIS). Antiviral Res. 2011 Aug; 91(2):150±3. https://doi.org/10. 1016/j.antiviral.2011.05.010 PMID: 21663768 8. Artenstein AW, VanCott TC, Mascola JR, Carr JK, Hegerich PA, Gaywee J, et al. Dual infection with human immunodeficiency virus type 1 of distinct envelope subtypes in humans. J Infect Dis. 1995 Apr; 171(4):805±10. PMID: 7706806 9. van Harmelen J, Wood R, Lambrick M, Rybicki EP, Williamson AL, Williamson C. An association between HIV-1 subtypes and mode of transmission in Cape Town, South Africa. AIDS Lond Engl. 1997 Jan; 11(1):81±7. 10. Tscherning C, Alaeus A, Fredriksson R, Bjo È rndal A, Deng H, Littman DR, et al. Differences in chemo- kine coreceptor usage between genetic subtypes of HIV-1. Virology. 1998 Feb 15; 241(2):181±8. https://doi.org/10.1006/viro.1997.8980 PMID: 9499793 11. Kanki PJ, Hamel DJ, Sankale  JL, Hsieh C c, Thior I, Barin F, et al. Human immunodeficiency virus type 1 subtypes differ in disease progression. J Infect Dis. 1999 Jan; 179(1):68±73. https://doi.org/10.1086/ 314557 PMID: 9841824 12. Apetrei C, Descamps D, Collin G, Loussert-Ajaka I, Damond F, Duca M, et al. Human immunodefi- ciency virus type 1 subtype F reverse transcriptase sequence and drug susceptibility. J Virol. 1998 May; 72(5):3534±8. PMID: 9557632 13. Taylor BS, Sobieszczyk ME, McCutchan FE, Hammer SM. The challenge of HIV-1 subtype diversity. N Engl J Med. 2008 Apr 10; 358(15):1590±602. https://doi.org/10.1056/NEJMra0706737 PMID: 14. Alaeus A, Lidman K, So È nnerborg A, Albert J. Subtype-specific problems with quantification of plasma HIV-1 RNA. AIDS Lond Engl. 1997 Jun; 11(7):859±65. 15. Parekh B, Phillips S, Granade TC, Baggs J, Hu DJ, Respess R. Impact of HIV type 1 subtype variation on viral RNA quantitation. AIDS Res Hum Retroviruses. 1999 Jan 20; 15(2):133±42. https://doi.org/10. 1089/088922299311556 PMID: 10029245 16. Brennan CA, Bodelle P, Coffey R, Harris B, Holzmayer V, Luk K-C, et al. HIV global surveillance: foun- dation for retroviral discovery and assay development. J Med Virol. 2006; 78 Suppl 1:S24±9. 17. Gilbert MTP, Rambaut A, Wlasiuk G, Spira TJ, Pitchenik AE, Worobey M. The emergence of HIV/AIDS in the Americas and beyond. Proc Natl Acad Sci U S A. 2007 Nov 20; 104(47):18566±70. https://doi. org/10.1073/pnas.0705329104 PMID: 17978186 18. Gray RR, Tatem AJ, Lamers S, Hou W, Laeyendecker O, Serwadda D, et al. Spatial phylodynamics of HIV-1 epidemic emergence in east Africa. AIDS Lond Engl. 2009 Sep 10; 23(14):F9±17. 19. Esbjo È rnsson J, Mild M, Månsson F, Norrgren H, Medstrand P. HIV-1 molecular epidemiology in Guinea- Bissau, West Africa: origin, demography and migrations. PloS One. 2011; 6(2):e17025. https://doi.org/ 10.1371/journal.pone.0017025 PMID: 21365013 20. Gonza  lez-Alba JM, Holguõ Ân A, Garcia R, Garcõ Âa-Bujalance S, Alonso R, Sua  rez A, et al. Molecular Sur- veillance of HIV-1 in Madrid, Spain: a Phylogeographic Analysis !. J Virol. 2011 Oct; 85(20):10755±63. https://doi.org/10.1128/JVI.00454-11 PMID: 21795343 21. Aldous JL, Pond SK, Poon A, Jain S, Qin H, Kahn JS, et al. Characterizing HIV transmission networks across the United States. Clin Infect Dis Off Publ Infect Dis Soc Am. 2012 Oct; 55(8):1135±43. 22. Grabowski MK, Redd AD. Molecular tools for studying HIV transmission in sexual networks. Curr Opin HIV AIDS. 2014 Mar; 9(2):126±33. https://doi.org/10.1097/COH.0000000000000040 PMID: 24384502 23. Cuevas MT, Muñoz-Nieto M, Thomson MM, Delgado E, Iribarren JA, Cilla G, et al. HIV-1 Transmission Cluster With T215D Revertant Mutation Among Newly Diagnosed Patients From the Basque Country, Spain: JAIDS J Acquir Immune Defic Syndr. 2009 May; 51(1):99±103. https://doi.org/10.1097/QAI. 0b013e318199063e PMID: 19282784 24. Vega Y, Delgado E, Ferna  ndez-Garcõ Âa A, Cuevas MT, Thomson MM, Montero V, et al. Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004±2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations. PLOS ONE. 2015 May 26; 10(5):e0125699. https://doi.org/10.1371/journal.pone.0125699 PMID: 26010948 25. Perez-Parra S, Chueca-Porcuna N, Alvarez-Estevez M, Pasquau J, Omar M, Collado A, et al. [Study of human immunodeficiency virus transmission chains in Andalusia: Analysis from baseline antiretroviral resistance sequences.]. Enferm Infecc Microbiol Clin. 2015 Jan 31; 26. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinforma Oxf Engl. 2007 Nov 1; 23(21):2947±8. 27. Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylo- genetic trees. In: Gateway Computing Environments Workshop (GCE), 2010. 2010. p. 1±8. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 17 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain 28. Hillis DM, Bull JJ. An Empirical Test of Bootstrapping as a Method for Assessing Confidence in Phyloge- netic Analysis. Syst Biol. 1993 Jun 1; 42(2):182±92. 29. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999 Jan; 73(1):152±60. PMID: 9847317 30. Leigh Brown AJ, Lycett SJ, Weinert L, Hughes GJ, Fearnhill E, Dunn DT, et al. Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J Infect Dis. 2011 Nov; 204 (9):1463±9. https://doi.org/10.1093/infdis/jir550 PMID: 21921202 È Â 31. Kouyos RD, von Wyl V, Yerly S, Boni J, Taffe P, Shah C, et al. Molecular epidemiology reveals long- term changes in HIV type 1 subtype B transmission in Switzerland. J Infect Dis. 2010 May 15; 201 (10):1488±97. https://doi.org/10.1086/651951 PMID: 20384495 32. Esbjo È rnsson J, Mild M, Audelin A, Fonager J, Skar H, Bruun Jørgensen L, et al. HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries. Virus Evol. 2016 Jan; 2(1):vew010. https://doi.org/10.1093/ve/vew010 PMID: 33. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012 Aug; 29(8):1969±73. https://doi.org/10.1093/molbev/mss075 PMID: 22367748 34. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed Phylogenetics and Dating with Confidence. PLoS Biol. 2006 Mar 14; 4(5):e88. https://doi.org/10.1371/journal.pbio.0040088 PMID: 16683862 35. Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005 May; 22(5):1185±92. https://doi.org/10.1093/ molbev/msi103 PMID: 15703244 36. Hue S, Brown AE, Ragonnet-Cronin M, Lycett SJ, Dunn DT, Fearnhill E, et al. Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions. AIDS Lond Engl. 2014 Aug 24; 28(13):1967±75. 37. Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, Kiuchi M, et al. Drug Resistance Muta- tions for Surveillance of Transmitted HIV-1 Drug-Resistance: 2009 Update. PLoS ONE [Internet]. 2009 Mar 6 [cited 2016 Sep 3]; 4(3). Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2648874/ 38. Alexiev I, Shankar A, Wensing AMJ, Beshkov D, Elenkov I, Stoycheva M, et al. Low HIV-1 transmitted drug resistance in Bulgaria against a background of high clade diversity. J Antimicrob Chemother. 2015; 70(6):1874±80. https://doi.org/10.1093/jac/dkv011 PMID: 25652746 39. De Mendoza C, Garrido C, Poveda E, Corral A, Zahonero N, Treviño A, et al. Changes in drug resis- tance patterns following the introduction of HIV type 1 non-B subtypes in Spain. AIDS Res Hum Retrovi- ruses. 2009 Oct; 25(10):967±72. https://doi.org/10.1089/aid.2008.0166 PMID: 19842792 40. Trends in Drug Resistance Prevalence in HIV-1±infected Children in Madrid (PDF Download Available) [Internet]. ResearchGate. [cited 2017 Mar 15]. Available from: https://www.researchgate.net/ publication/229074263_Trends_in_Drug_Resistance_Prevalence_in_HIV-1-infected_Children_in_ Madrid 41. Ferna  ndez-Garc õa A, Cuevas MT, Vinogradova A, Rakhmanova A, Pe  rez-Alvarez L, de Castro RO, et al. Near full-length genome characterization of a newly identified HIV type 1 subtype F variant circu- lating in St. Petersburg, Russia. AIDS Res Hum Retroviruses. 2009 Nov; 25(11):1187±91. https://doi. org/10.1089/aid.2009.0140 PMID: 19943791    42. Delgado E, Cuevas MT, Domõnguez F, Vega Y, Cabello M, Fernandez-Garcõa A, et al. 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. PloS One. 2015; 10 (11):e0143325. https://doi.org/10.1371/journal.pone.0143325 PMID: 26599410   43. Thomson MM, Fernandez-Garcõ Âa A, Delgado E, Vega Y, Dõ Âez-Fuertes F, Sanchez-Martõ Ânez M, et al. Rapid expansion of a HIV-1 subtype F cluster of recent origin among men who have sex with men in Galicia, Spain. J Acquir Immune Defic Syndr 1999. 2012 Mar 1; 59(3):e49±51. 44. Yebra G, de Mulder M, del Romero J, Rodrõ Âguez C, Holguõ Ân A. HIV-1 non-B subtypes: High transmitted NNRTI-resistance in Spain and impaired genotypic resistance interpretation due to variability. Antiviral Res. 2010 Feb; 85(2):409±17. https://doi.org/10.1016/j.antiviral.2009.11.010 PMID: 20004217 45. Yebra G, de Mulder M, Pe  rez-Elõ Âas MJ, Pe  rez-Molina JA, Gala  n JC, Llenas-Garcõ Âa J, et al. Increase of transmitted drug resistance among HIV-infected sub-Saharan Africans residing in Spain in contrast to the native population. PloS One. 2011; 6(10):e26757. https://doi.org/10.1371/journal.pone.0026757 PMID: 22046345 46. Bracho MA, Sentandreu V, Alastrue  I, Belda J, Juan A, Ferna  ndez-Garcõ Âa E, et al. Emerging trends in CRF02_AG variants transmission among men who have sex with men in Spain. J Acquir Immune Defic Syndr 1999. 2014 Mar 1; 65(3):e130±3. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 18 / 19 High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain 47. Holguõ Ân A, de Mulder M, Yebra G, Lo  pez M, Soriano V. Increase of non-B subtypes and recombinants among newly diagnosed HIV-1 native Spaniards and immigrants in Spain. Curr HIV Res. 2008 Jun; 6 (4):327±34. PMID: 18691031 48. Ferna  ndez-Garc õa A, Cuevas MT, Muñoz-Nieto M, Ocampo A, Pinilla M, Garcõ Âa V, et al. Development of a panel of well-characterized human immunodeficiency virus type 1 isolates from newly diagnosed patients including acute and recent infections. AIDS Res Hum Retroviruses. 2009 Jan; 25(1):93±102. https://doi.org/10.1089/aid.2008.0174 PMID: 19113978   49. de Felipe B, Perez-Romero P, Abad-Fernandez M, Fernandez-Cuenca F, Martinez-Fernandez FJ, Trastoy M, et al. Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000±2010. Virol J. 2011; 8:416. https://doi.org/10.1186/1743-422X- 8-416 PMID: 21871090 50. Alvarez M, Garcõ Âa F, Martõ Ânez NM, Garcõ Âa F, Bernal C, Vela CM, et al. Introduction of HIV type 1 non-B subtypes into Eastern Andalusia through immigration. J Med Virol. 2003 May; 70(1):10±3. https://doi. org/10.1002/jmv.10368 PMID: 12629637 51. Instituto Nacional de Estadõ Âstica (INE) [Internet]. Available from: www.ine.es 52. Brenner BG, Roger M, Stephens D, Moisi D, Hardy I, Weinberg J, et al. Transmission Clustering Drives the Onward Spread of the HIV Epidemic Among Men Who Have Sex With Men in Quebec. J Infect Dis. 2011 Oct 1; 204(7):1115±9. https://doi.org/10.1093/infdis/jir468 PMID: 21881127 53. Delgado E, Cuevas MT, Vega Y, Montero V, Sa  nchez M, Carrera C, et al. Identificacio  n de un cluster de subtipo A que se transmite entre hombres que tienen relaciones sexuales con hombres en diversas comunidades auto  nomas de España. In Ma  laga: VI Congreso Nacional de GESIDA y 8.a Reunio  n Doc- ente de la RIS (SEIMC); 2014. p. 12. Available from: http://www.gesida-seimc.org/contenidos/ congresos/anteriores/2014/gesida2014-VIcongresocomunicaciones.pdf 54. Patiño Galindo JA, Torres-Puente M, Gimeno C, Ortega E, Navarro D, Galindo MJ, et al. Expansion of the CRF19_cpx Variant in Spain. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2015 Aug; 69:146±9. PLOS ONE | https://doi.org/10.1371/journal.pone.0186928 October 30, 2017 19 / 19

Journal

PLoS ONEPubmed Central

Published: Oct 30, 2017

There are no references for this article.