Abstract
In Western Europe, the HIV-1 epidemic among men who have sex with men (MSM) is dominated by subtype B. However, recently, other genetic forms have been reported to circulate in this population, as evidenced by their grouping in clusters predominantly comprising European individuals. Here we describe four large HIV-1 non-subtype B clusters spreading among MSM in Spain. Samples were collected in 9 regions. A pol fragment was amplified from plasma RNA or blood-extracted DNA. Phylogenetic analyses were performed via maximum likelihood, including database sequences of the same genetic forms as the identified clusters. Times and locations of the most recent common ancestors (MRCA) of clusters were estimated with a Bayesian method. Five large non-subtype B clusters associated with MSM were identified. The largest one, of F1 subtype, was reported previously. The other four were of CRF02_AG (CRF02_1; n = 115) and subtypes A1 (A1_1; n = 66), F1 (F1_3; n = 36), and C (C_7; n = 17). Most individuals belonging to them had been diagnosed of HIV-1 infection in the last 10 years. Each cluster comprised viruses from 3 to 8 Spanish regions and also comprised or was related to viruses from other countries: CRF02_1 comprised a Japanese subcluster and viruses from 8 other countries from Western Europe, Asia, and South America; A1_1 comprised viruses from Portugal, United Kingom, and United States, and was related to the A1 strain circulating in Greece, Albania and Cyprus; F1_3 was related to viruses from Romania; and C_7 comprised viruses from Portugal and was related to a virus from Mozambique. A subcluster within CRF02_1 was associated with heterosexual transmission. Near full-length genomes of each cluster were of uniform genetic form. Times of MRCAs of CRF02_1, A1_1, F1_3, and C_7 were estimated around 1986, 1989, 2013, and 1983, respectively. MRCA locations for CRF02_1 and A1_1 were uncertain (however initial expansions in Spain in Madrid and Vigo, respectively, were estimated) and were most probable in Bilbao, Spain, for F1_3 and Portugal for C_7. These results show that the HIV-1 epidemic among MSM in Spain is becoming increasingly diverse through the expansion of diverse non-subtype B clusters, comprising or related to viruses circulating in other countries.
Introduction
HIV-1 exhibits a characteristic high genetic variability, derived from elevated mutation and recombination rates. Through these mechanisms, the main (M) HIV-1 group, causative of the pandemic, has evolved into multiple genetic forms, designated subtypes, of which nine have been identified, subsubtypes, circulating recombinant forms (CRFs), of which 93 are currently recognized (HIV Sequence Database, 2019; Reis et al., 2019), and unique recombinant forms (URFs). The most globally prevalent HIV-1 genetic form is subtype C, estimated to represent around 47% worldwide infections, followed, in this order, by subtype B, subtype A, CRF02_AG, CRF01_AE, subtype G, and subtype D, with each of the remaining genetic forms estimated to represent less than 1% of global infections (Hemelaar et al., 2018).
In Western Europe, the predominant HIV-1 genetic form is subtype B, which was initially introduced among MSM and persons who inject drugs (PWID) (Lukashov et al., 1996; Casado et al., 2000; Kuiken et al., 2000; Hué et al., 2005; Beloukas et al., 2016). In early descriptions of HIV-1 genetic diversity in Western Europe, non-subtype B genetic forms were restricted to heterosexually infected immigrants coming from areas where those clades predominate, mainly sub-Saharan Africans, and European individuals epidemiologically linked to people from such areas (Fransen et al., 1996; Op de Coul et al., 1998; Thomson and Nájera, 2001). The first reports of HIV-1 non-subtype B genetic forms circulating in Western Europe among individuals without known epidemiological links to other geographic areas described the circulation of CRF01_AE among PWID in Finland (Liitsola et al., 2000) and of subtype G and CRF14_BG among a minority of HIV-1-infected PWID in the region of Galicia, Northwest Spain (Thomson et al., 2001; Delgado et al., 2002). Subsequent studies showed that the genetic forms circulating in Galicia derived from a subtype G variant widely circulating in Portugal, transmitted both through sexual contact and among PWID (Esteves et al., 2002, 2003; Palma et al., 2007; Carvalho et al., 2015). In recent years, an increasing prevalence of non-subtype B infections has been observed in Western Europe, reflecting both their importation from other geographical areas and their circulation among the local population (Abecasis et al., 2013; Beloukas et al., 2016; Hemelaar et al., 2018; Paraskevis et al., 2019).
In the current HIV-1 epidemic in Western Europe, the predominant propagation mode is through sexual contact among MSM (European Centre for Disease Prevention and Control, 2017 and WHO Regional Office for Europe 2017; Núñez et al., 2018), a population in which a resurgence of the HIV-1 epidemic has been observed since the 2000s, which is part of a global phenomenon (Bezemer et al., 2008; Beyrer et al., 2012). This has been accompanied by the emergence of phylogenetically identifiable transmission clusters, whose expansion is mostly driven by individuals with recent infection who are unaware of their HIV status (Lewis et al., 2008; Cuevas et al., 2009; Bezemer et al., 2010; Chalmet et al., 2010; Fisher et al., 2010; Zehender et al., 2010; Ambrosioni et al., 2012; Frange et al., 2012; Thomson et al., 2012; Audelin et al., 2013; Delgado et al., 2015; Esbjörnsson et al., 2016; Hoenigl et al., 2016; Chaillon et al., 2017; Patiño-Galindo et al., 2017; Parczewski et al., 2017; Verhofstede et al., 2018; Paraskevis et al., 2019). As expected, most clusters are of subtype B, but multiple instances of propagation of other HIV-1 genetic forms among European MSM have also been reported. These include subtypes A1 (Lai et al., 2016; Ragonnet-Cronin et al., 2016), C (de Oliveira et al., 2010; Lai et al., 2014; Ragonnet-Cronin et al., 2016), and F1 (Castro et al., 2010; Lai et al., 2012; Thomson et al., 2012; Delgado et al., 2015; Vinken et al., 2017; Verhofstede et al., 2018); CRF01_AE (von Wyl et al., 2011), CRF02_AG (Giuliani et al., 2013; Brand et al., 2014; Dauwe et al., 2015; Tamalet et al., 2015; Chaillon et al., 2017; Verhofstede et al., 2018), CRF17_BF (Fabeni et al., 2015), CRF19_cpx (Patiño-Galindo et al., 2015; González-Domenech et al., 2018; Pérez-Parra et al., 2018), CRF50_A1D (Foster et al., 2014), CRF56_cpx (Leoz et al., 2013), and CRF60_BC (Monno et al., 2012; Simonetti et al., 2014). However, the expansion of these clades has had a limited impact on the overall genetic composition of the HIV-1 epidemic among MSM in Western Europe, which is still largely dominated by subtype B. The only exception, though geographically restricted, is a large F1 subtype cluster of Brazilian ancestry, which represented a substantial proportion of new HIV-1 diagnoses among MSM in Northwest Spain (Thomson et al., 2012; Delgado et al., 2015). Here we describe four additional large non-subtype B clusters expanding among MSM in Spain, of CRF02_AG and of subtypes A1, F1 and C, each circulating in several Spanish regions and related to viruses from other countries.
Materials and Methods
Samples
Plasma or whole blood samples were collected from 1999 to 2018 from HIV-1-infected individuals attended at hospitals from 15 provinces from 9 regions of Spain (Basque Country, Galicia, Navarre, Castilla y León, La Rioja, Madrid, Castilla-La Mancha, Aragón, and Comunidad Valenciana). The regional sample representativeness is variable, being the greatest in the regions of Basque Country, where all public hospitals participated, and Galicia, where all but one public hospitals participated. The study was approved by the Bioethics and Animal Well-being Committee of Instituto de Salud Carlos III, Majadahonda, Madrid, Spain. Written informed consent was obtained from all participants in the study.
RNA and DNA Extraction, RT-PCR Amplification, and Sequencing
Amplification of HIV-1 fragments was done either from plasma RNA or from DNA extracted from whole blood. RNA was extracted from 1 ml plasma using Nuclisens EasyMAG kit (bioMérieux, Marcy l’Etoile, France), following the manufacturer’s instructions. DNA was extracted from 200 μl blood using QIAmp DNA Blood mini kit (Qiagen, Hilden, Germany), following the manufacturer’s instructions. An HIV-1 PR-RT fragment (approximately 1.4 kb) was amplified by RT-PCR followed by nested PCR, in the case of RNA, or by nested PCR, in the case of DNA, using previously reported primers (Delgado et al., 2015). In selected samples, NFLG amplification was done in four overlapping segments, as described (Delgado et al., 2002; Sierra et al., 2005), using RNA extracted either from plasma or from the primary isolate’s culture supernatant grown from plasma using a previously described protocol (Delgado et al., 2012). Direct sequencing of the amplified products was done using an automated capillary sequencer. Sequence electropherograms were assembled and edited with Seqman (DNASTAR, Madison, WI, United States). Newly obtained sequences are deposited in GenBank under accessions MK177651-MK177824 (PR-RT sequences) and MK177825-MK177829, KT276258, KY496622, and KY989952 (NFLG sequences).
Phylogenetic Sequence Analyses
Sequences were aligned with MAFFT v.7 (Katoh and Standley, 2013). Initial trees with all sequences were constructed with the approximate ML method implemented in FastTree v2.1.10 (Price et al., 2010) using the general time reversible (GTR)+CAT evolutionary model and assessment of node support with Shimodaira-Hasegawa (SH)-like local support values. Transmission clusters were defined as those supported by SH-like values ≥ 0.95 comprising four or more individuals, with a majority of them being native Spanish. To determine whether the identified clusters were still supported when globally sampled sequences were included and to identify viruses from other areas belonging or phylogenetically related to them, the analyses were repeated including all PR-RT sequences of 1 kb or longer of the same genetic form available at the Los Alamos HIV Sequence Database (HIV Sequence Database, 2019), downloaded with the option “one sequence per patient.” Support for clusters thus identified with FastTree was subsequently assessed with phylogenetic trees constructed via ML with PhyML v3.0 (Guindon et al., 2005), using the best-fit evolutionary model selected by Smart Model Selection (SMS) program (Lefort et al., 2017) and heuristic searches based on subtree pruning and regrafting (SPR) moves, with estimation of node support by the approximate likelihood ratio test (aLRT), SH-like procedure (Anisimova and Gascuel, 2006; Guindon et al., 2010). To keep computational times reasonable, in the analyses with PhyML only several hundred sequences (around 200–400) branching most closely to the clusters of interest in the previous FastTree analyses, together with PR-RT sequences from NFLG from databases, were included. Clusters were confirmed if in the analyses including database sequences with both FastTree and PhyML their node support values were ≥ 0.95. Trees were visualized with Dendroscope 3 (Huson and Scornavacca, 2012) or FigTree v1.4.21. Intersubtype recombination was analyzed by bootscanning with Simplot v3.5 (Lole et al., 1999), with tree construction with the neighbor-joining method, using the Kimura 2-parameter substitution model and windows of 400 or 600 nt moving in 20 nt steps.
Antiretroviral Drug Resistance Determination
Antiretroviral (ARV) drug resistance was analyzed with the Calibrated Population Resistance Tool (Gifford et al., 2009).
Temporal and Geographic Estimations of Cluster Origins
To estimate times of the most recent common ancestors (tMRCA) of clusters and their most probable geographical locations, a Bayesian MCMC coalescent method as implemented in BEAST v1.8.1 (Drummond et al., 2012) was used. Prior to these analyses, the existence of temporal signal in the datasets was assessed by an analysis of root-to-tip distances against dates of sampling using TempEst (Rambaut et al., 2016). For the BEAST analyses, we used all PR-RT sequences ≥ 1 kb from each cluster and related sequences, as determined in the ML phylogenetic analyses, excluding sequences without data on year or location of sample collection. PR-RT sequences derived from NFLG sequences of the corresponding genetic form downloaded from the HIV Sequence Database were also included in these analyses, using not more than 5 sequences per country of sampling. In the case of the CRF02_1 cluster, 40 randomly selected CRF02_AG PR-RT sequences lacking drug resistance mutations downloaded from the Los Alamos HIV Sequences database were included, since using those derived from NFLG resulted in relatively low r2 values in the TempEst analysis. We chose an HKY substitution model with gamma-distributed among-site rate heterogeneity and two partitions in codon positions (1st+2nd; 3rd) (Shapiro et al., 2006); uniform priors were used for absolute substitution rates (0-0.02 sub/site/year) and relative substitution rates in codon positions 1st+2nd and 3rd (0-10); we also used an uncorrelated lognormal relaxed clock model and a Bayesian skyline plot demographic model (Drummond et al., 2005). Each MCMC chain was run for 100 million to 200 million generations, sampling every 5,000 generations. MCMC convergence and effective sample sizes (ESS) were checked with Tracer v.1.62, ensuring that the ESS of each parameter was > 200. Results were summarized with a maximum clade credibility (MCC) tree, using TreeAnnotator v1.8.1, after removal of a 50% burn-in. The MCC trees were visualized with FigTree v1.4.2.1 Parameter uncertainty was summarized in the 95% highest posterior density (HPD) intervals.
Statistical Analyses
Differences in clustering frequency between MSM and heterosexually infected individuals and changes in proportions of non-subtype B infections along time in newly diagnosed sexually infected individuals were analyzed with Fisher’s exact test. Only native Spish individuals were included in these analyses in order to focus on locally circulating strains, minimizing the confounding effect of imported HIV-1 variants acquired in other countries. Numbers of individuals used in these analyses were 7080 in the first and 2060 in the second.
Results
Using HIV-1 PR-RT sequences from 10,506 individuals obtained by us (whose data are summarized in Supplementary Table 1), we identified 320 phylogenetic clusters comprising 4 or more individuals, 247 of subtype B and 73 of other genetic forms. Belonging to a cluster was more frequent among Spanish MSM than among heterosexually-infected Spanish individuals (68% vs. 31%; p = 7.7 × 10-6) (Supplementary Figure 1). Differences were also significant when subtype B and non-subtype B clusters were analyzed separately. An increase along time of non-subtype B infections among newly diagnosed Spanish individuals was also observed. Statistically significant increases were observed in Spanish MSM from 2005–2009 to 2010–2014 (from 10.9% to 26.7%; p = 1.4 × 10-6) and among heterosexually infected Spanish individuals from 2010–2014 to 2015–2018 (from 28.5% to 39.7%; p = 0.0097) (Supplementary Figure 2).
Of the non-subtype B clusters, 5 large ones (here defined as those comprising 10 or more individuals) were associated with MSM. The largest one, of F1 subtype (currently comprising 187 individuals), was reported previously (Thomson et al., 2012; Delgado et al., 2015). The other four were of CRF02_AG and of subtypes A1, F1, and C, henceforth designated CRF02_1, A1_1, F1_3, and C_7, respectively. All four clusters were well supported when the analyses were repeated with FastTree including all PR-RT sequences > 1 kb of the respective genetic forms available at the Los Alamos HIV Sequence database (with numbers from 1,318 of subtype F1 to 22,762 of subtype C) and with PhyML including several hundred (from 312 for subtype C to 513 for CRF02_AG) database sequences branching closer to each cluster in the FastTree analysis and derived from NFLG sequences. These analyses also allowed to identify viruses from databases belonging or closely related to the clusters. Epidemiological data of samples studied by us belonging to these clusters and to subclusters within them are shown in Table 1.
Table 1
| Cluster/ Subcluster | N | Gender | Transmission route | Region | Country of origin† | Years of HIV-1 diagnosis | |||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| M | F | ND | MSM | MSx | HT | PWID | ND | BC | GA | NA | MD | AR | LR | CV | CL | CM | Spain | LAm | Other | ND | |||
| CRF02_1 | 67 | 58 | 8 | 1 | 30 | 12 | 18 | 1 | 6 | 31 | 3 | 3 | 10 | 15 | 3 | 1 | 1 | 53 | 10 | 4 | 2008–2018 | ||
| 1 | 49 | 47 | 1 | 1 | 28 | 12 | 5 | 1 | 3 | 26 | 3 | 1 | 9 | 6 | 3 | 1 | 40 | 7 | 2 | 2008–2018 | |||
| 1_1 | 17 | 17 | 11 | 4 | 1 | 1 | 15 | 2 | 14 | 2 | 1 | 2008–2018 | |||||||||||
| 1_2 | 7 | 7 | 3 | 1 | 3 | 1 | 6 | 6 | 1 | 2016–2018 | |||||||||||||
| 3 | 10 | 5 | 5 | 10 | 9 | 1 | 8 | 1 | 1 | 2016–2018 | |||||||||||||
| A1_1 | 54 | 50 | 3 | 1 | 36 | 6 | 6 | 6 | 6 | 22 | 18 | 3 | 5 | 37 | 10 | 5 | 2 | 2006–2018 | |||||
| 1 | 30 | 27 | 3 | 18 | 4 | 3 | 5 | 6 | 18 | 1 | 5 | 19 | 5 | 5 | 1 | 2013–2017 | |||||||
| 2 | 16 | 15 | 1 | 12 | 1 | 2 | 1 | 15 | 1 | 13 | 2 | 1 | 2009–2018 | ||||||||||
| 2_1 | 10 | 10 | 8 | 1 | 1 | 9 | 1 | 2016–2018 | |||||||||||||||
| F1_3 | 36 | 36 | 17 | 10 | 3 | 6 | 30 | 3 | 3 | 29 | 6 | 1 | 2014–2018 | ||||||||||
| C_7 | 14 | 14 | 5 | 6 | 3 | 4 | 9 | 1 | 10 | 1 | 3 | 2009–2017 | |||||||||||
| 1 | 8 | 8 | 5 | 2 | 1 | 7 | 6 | 1 | 1 | 2012–2017 | |||||||||||||
| 2 | 4 | 4 | 3 | 1 | 4 | 3 | 1 | 2009–2015 | |||||||||||||||
Epidemiological data of clusters∗.
∗Only data from samples processed by the authors are included. †Country of which the individual is native. ND, no data. HT, heterosexual. MSx, male, unspecified sexual. BC, Basque Country. MD, Madrid. GA, Galicia. NA, Navarre. AR, Aragon. CL, Castilla y León. CM, Castilla-La Mancha. CV, Comunidad Valenciana. LR, La Rioja. LAm, Latin American country.
CRF02_1 Cluster
CRF02_1 comprised 115 individuals, 67 studied by us and 48 whose sequences were retrieved from databases (Figure 1). Most samples were collected in Spain, but there were also samples from Japan, Switzerland, United Kingdom, Ecuador, Netherlands, Sweden, Germany, Malaysia, and Hong Kong. Samples from Spain were from 8 regions, mainly from Basque Country, Madrid, and Aragon.
FIGURE 1
Years of HIV-1 diagnosis were 2008–2018.
CRF02_1 comprised three subclusters, designated with numerals 1-3. CRF02_1_1 comprised 60 individuals from 7 Spanish regions, two from United Kingdom and one from Germany. Within it, two subsubclusters were distinguished, associated with the Basque Country and the city of Zaragoza, Aragón, respectively; CRF02_1_2 comprised 13 individuals from Japan; and CRF02_1_3 comprised 11 individuals, 10 of them from Zaragoza.
Among samples studied by us, 88% were from men. Transmission route was sexual in the great majority, and, among those infected via sexual contact, 51% were self-declared MSM. Interestingly, all 10 individuals from subcluster CRF02_1_3 studied by us (5 women and 5 men, all from Zaragoza) were infected via heterosexual contact.
Most individuals in this and in the other clusters here described were native Spanish (Table 1).
A1_1 Cluster
A1_1 comprised 66 individuals, 54 studied by us and 12 with sequences deposited in databases (7 from United Kingdom, 4 from Portugal, and 1 from United States) (Figure 2). A1_1 was also related to viruses from the A1 subtype lineage circulating in Greece, Albania, and Cyprus (Ciccozzi et al., 2005; Paraskevis et al., 2007; Pineda-Peña et al., 2018). Samples from Spain were collected in 5 regions, mainly in Galicia and Navarre.
FIGURE 2
Forty five of 47 individuals with available data were diagnosed in 2012–2018.
Within A1_1, there were two main subclusters: A1_1_1, comprising all individuals from Navarre, Basque Country, and La Rioja, and 1 from Madrid; and, A1_1_2 comprising individuals mostly from Galicia, with a majority grouping in a subsubcluster.
Most (93%) individuals in the cluster were men, with predominance of MSM.
F1_3 Cluster
F1_3 comprised 36 individuals. There were no sequences from databases belonging to F1_3, but viruses from Romania, most of them transmitted sexually (Niculescu et al., 2015), were closely related to it (Figure 3). Five viruses from Spain (4 of them sequenced by us), 3 from United Kingdom, and 1 from Poland also branched close to F1_3, interspersed among the Romanian samples. Most F1_3 samples were from Basque Country.
FIGURE 3
All infections were diagnosed in 2014 or later. All individuals in F1_3 were men. Transmission route was sexual in all for which data were available, with a majority of MSM.
C_7 Cluster
C_7 comprised 17 individuals, including 14 studied by us and 3 from Portugal, whose sequences were retrieved from databases. One sequence from Mozambique branched basally to C_7 (Figure 4). C_7 comprised two subclusters, C_7_1, comprising all but one samples from Galicia, and C_7_2, comprising all 4 samples from Basque Country. All but one had been diagnosed in 2012–2017.
FIGURE 4
All but one Spanish samples were from Galicia or Basque Country. All individuals in C_7 were sexually infected men, with 5 being self-declared MSM.
Bayesian Analyses and Temporal and Geographic Estimations
To estimate the temporal and geographic origins of clusters and subclusters, Bayesian coalescent analyses were performed with PR-RT sequences, summarizing the posterior distribution of trees with MCC trees. Prior to these analyses, temporal signal was analyzed, revealing a clock-like structure in all four datasets used for subsequent analyses (r2= 0.4395 in CRF02_1, r2= 0.5399 in A1_1, r2= 0.7175 in F1_3, and r2= 0.3518 in C_7), indicating that the datasets contained sufficient temporal structure for reliable estimation of divergence times.
In the Bayesian analyses, all four clusters previously defined via ML were supported by node PP values > 0.98.
The tMRCA of the entire CRF02_1 cluster was estimated around 1986, but the location of the MRCA was uncertain, since the most probable one was supported by a PP < 0.5. However, the location of the MRCA of the subcluster comprising all but the 8 most basal sequences (collected in Switzerland, Sweden, Germany, Malaysia, and Hong Kong), with tMRCA around 1988, had a strong support in Madrid, with a PP of 0.975 (Figure 5). Subcluster CRF02_1_1 emerged around 2002 in Madrid, CRF02_1_2 around 2006 in Japan, and CRF02_1_3 around 2006 in Zaragoza.
FIGURE 5
The estimated tMRCA of the entire A1_1 cluster was around 1989, but the location of the MRCA was uncertain, since the most probable one was supported by a PP < 0.5. However, the location of the MRCA of the subcluster comprising all but the 6 most basal sequences, with tMRCA around 1994, was supported by a PP of 0.814 in the city of Vigo, Galicia. The emergence of subcluster A1_1_1 was around 2010, with highest location PP in Pamplona, and that of A1_1_2 was around 2004 in Vigo, Galicia (Figure 6). It should be pointed out that 2 samples from Portugal and 6 from United Kingdom could not be used in the Bayesian analyses, since no collection year was available for them.
FIGURE 6
The estimated tMRCA of F1_3 was around 2013 in the city of Bilbao, with a strongly supported ancestry in Romania (Figure 7).
FIGURE 7
Finally, the estimated tMRCA of C_7 was around 1983. Its most probable origin was Portugal, but with a PP support for the entire cluster of only 0.6. However, the location of the MRCA of the subcluster comprising all but the most basal sample had a strong support in Portugal (PP = 0.91) (Figure 8). Subcluster C_7_1 emerged around 2008 in Vigo and C_7_2 around 2004 in Vitoria, Basque Country.
FIGURE 8
ARV Drug Resistance Mutations
In CRF02_1, 5 sequences (4 from Spain and one from Germany) had ARV drug resistance mutations. One was from a drug-experienced individual in therapeutic failure with multiple drug resistance mutations, and two, with K101E and K103N, respectively, mutations of resistance to non-nucleoside reverse transcriptase inhibitors (NNRTI), were from drug-naïve individuals. The other two, with K103N and K101E mutations, respectively, were from database sequences without data on drug treatment. In A1_1, one database sequence from United Kingdom had Y188C and G190A NNRTI resistance mutations. In C_7, all but 2 sequences had L90M mutation associated with protease inhibitor drug resistance; all 13 Spanish sequences with this mutation were from drug-naïve individuals.
Near Full-Length Genome Sequences
To determine whether the viruses from the clusters were of uniform genetic form all along their genomes or were interclade recombinants, two NFLG sequences were obtained for each cluster, either from plasma RNA (P2648, P3075, P4496, P4346, and P4476) or from RNA extracted from culture supernatant (NA0048_2, X3303_2, and X3988). An additional NFLG from the A1_1 cluster (X2110, GenBank accession FJ670523) had been obtained previously by us (Cuevas et al., 2010). Bootscan analyses showed that all were of uniform genetic form along their genomes (Figure 9). We note that NFLG sequences of two CRF02_1 viruses from United Kingdom (Figure 1) are also available at sequence databases (Yebra et al., 2018; HIV Sequence Database, 2019) (GenBank accessions MF109381, MF109550).
FIGURE 9
Discussion
In Western Europe, subtype B has been largely predominant among MSM since the early HIV-1 epidemic, but in recent years other genetic forms have been reported to be circulating in this population, as evidenced by their grouping in phylogenetic clusters comprising mostly European individuals. In this study, based on a large dataset from Spain, we found an increase in proportions of non-subtype B infections among MSM in recent years (Supplementary Figure 2) and higher clustering frequency in MSM compared to heterosexually infected individuals (Supplementary Figure 1). Among clusters associated with MSM, five large ones were of non-subtype B genetic forms, one of which, of F1 subtype, was reported previously (Thomson et al., 2012; Delgado et al., 2015). The other four were of CRF02_AG and subtypes A1, F1, and C, for which here we analyze epidemiological correlations, estimated emergence times and places, NFLGs, and drug resistance mutations.
The CRF02_AG cluster (CRF02_1) comprised 115 individuals, including 67 studied by us and 48 whose sequences were retrieved from databases, making it one of the largest non-subtype B clusters circulating among MSM reported to date in Western Europe (Delgado et al., 2015; Vinken et al., 2017). CRF02_AG is the predominant HIV-1 genetic form in most West African countries (Montavon et al., 2000; Hemelaar et al., 2018) and is common in West-Central Africa. It also propagates as a minor form in several Western European countries (Giuliani et al., 2013; Brand et al., 2014; Tamalet et al., 2015; Beloukas et al., 2016; Chaillon et al., 2017; Verhofstede et al., 2018), Tunisia (El Moussi et al., 2017), and Brazil (Delatorre et al., 2012), and in 2002 caused an outbreak among PWID in Uzbekistan (Carr et al., 2005), with subsequent dissemination to Kazakhstan (Eyzaguirre et al., 2007; Lapovok et al., 2014) and Russia (Moskaleychik et al., 2015), giving rise to CRF63_02A1 through recombination with the former Soviet Union subtype A variant (Baryshev et al., 2014; Shcherbakova et al., 2014). CRF02_AG has been reported to be one of the most common non-subtype B genetic forms in Western Europe (together with subtypes A1 and C) (Abecasis et al., 2013; Beloukas et al., 2016; Hemelaar et al., 2018) and in Spain (Yebra et al., 2012). The CRF02_AG cluster here described is not completely new, since a cluster of four individuals from the region of Valencia belonging to it was reported by other authors, who also noted that 9 sequences from databases, 7 from Spain and 2 from Ecuador, were related to the Valencian cluster (Bracho et al., 2014). However, the data here presented considerably enlarge the size and the geographic range of the cluster. CRF02_1 comprises viruses from 8 Spanish regions and from 9 other countries, from Western Europe, Asia and South America, with 13 Japanese viruses grouping in a monophyletic subcluster, indicating that it is circulating in this country. Although CRF02_1 is mainly associated with MSM, a subcluster comprising 11 individuals, 10 of them from the city of Zaragoza, propagates via heterosexual contact (Table 1). The origin of CRF02_1 is not recent, with a tMRCA estimated around 1986, with uncertain location, for the entire cluster, or 1988 in Madrid for the subcluster excluding the 8 most basal sequences, but its three major subclusters emerged in recent years, with tMRCAs in the 2000s.
The A1 subtype cluster (A1_1) is the second largest cluster here described, with 66 individuals, 54 studied by us and 12 with sequences in databases. A1 subtype circulates mainly in Eastern, Central and Western Africa (Hemelaar et al., 2018), all former Soviet Union (FSU) countries (Bobkova, 2013), Greece (Paraskevis et al., 2007), Albania (Ciccozzi et al., 2005), and Cyprus (Pineda-Peña et al., 2018), and as a minor form in India (Pandey et al., 2016), although some authors designate the variants circulating in Western Africa and FSU as distinct subsubtypes (A3 and A6, respectively) (Meloni et al., 2004; Foley et al., 2016). The lineage circulating in Greece and Albania, also detected in Cyprus, is of monophyletic origin, with estimated tMRCA around 1978 (Paraskevis et al., 2007). A1 subtype clusters have been reported in United Kingdom (Gifford et al., 2007; Hughes et al., 2009; Ragonnet-Cronin et al., 2016), Italy (Lai et al., 2016), Switzerland (von Wyl et al., 2011), and Portugal (Carvalho et al., 2015), but the one here reported is the largest reported to date in Western Europe. A1_1 comprises individuals from 5 Spanish regions and 3 other countries (United Kingdom, Portugal, and United States), and is related to the Greek-Albanian A1 lineage (Figure 2). Its origin is not recent, with estimated tMRCA around 1989, with uncertain location for the entire cluster, or around 1994 in Vigo, Galicia, excluding the 6 most basal sequences, but its two major subclusters are of recent origin, with tMRCAs around 2004 and 2010, respectively.
The F1 cluster (F1_3) comprises 36 individuals, all resident in Spain, most of them in Bilbao, Basque Country. Of the clusters here described, this is the one with the most recent origin, with estimated tMRCA around 2013 in Bilbao. It is currently increasing in size, with 6 individuals newly diagnosed in 2018. Subtype F1 is circulating in Central Africa, Brazil and Romania (Dumitrescu et al., 1994; Louwagie et al., 1994; Bandea et al., 1995; Apetrei et al., 1997; Op de Coul et al., 2000; Hemelaar et al., 2018), and F1 subtype clusters have been recently identified in Spain (Thomson et al., 2012; Delgado et al., 2015), Belgium (Vinken et al., 2017; Verhofstede et al., 2018), Switzerland (Castro et al., 2010), Italy (Lai et al., 2012), and Portugal (Carvalho et al., 2015). F1_3 is most closely related to F1 viruses from Romania (Figure 4), which are related to viruses circulating in Angola (Guimarães et al., 2009) and initially propagated among adults via sexual contact, with subsequent propagation among institutionalized children through contaminated injection equipment (Op de Coul et al., 2000; Bello et al., 2012), and more recently among PWID (Temereanca et al., 2013; Niculescu et al., 2015). The expansion of an F1 subtype cluster of Romanian ancestry in Spain has its counterpart in the recent expansion of CRF14_BG, originally described in Spain (Thomson et al., 2001; Delgado et al., 2002) and Portugal (Esteves et al., 2003; Abecasis et al., 2011; Bartolo et al., 2011), among PWID in Romania (Niculescu et al., 2015). The exchange of HIV-1 genetic forms between Romania and Spain reflects the presence of a large Romanian immigrant population in Spain (Instituto Nacional de Estadística, 2018), frequently traveling between both countries.
The subtype C cluster (C_7) comprises 17 viruses, 14 from three Spanish regions, predominantly from Galicia, and three from Portugal. It comprises a Galician and a Basque subclusters (Figure 4). This is not the first report of this cluster, as it was described by us when it comprised only 7 individuals, within the context of the description of HIV-1 clusters bearing ARV drug resistance mutations, noting that viruses belonging to it carry the L90M mutation of resistance to protease inhibitors (Vega et al., 2015). This mutation was found in 15 of 17 C_7 viruses. Subtype C is the most prevalent clade in the HIV-1 pandemic, circulating mainly in Southern and Eastern Africa, Southern Brazil, and South Asia (Hemelaar et al., 2018). In Western Europe, subtype C clusters have been reported in the United Kingdom (Hughes et al., 2009; de Oliveira et al., 2010; Ragonnet-Cronin et al., 2016), Italy (Lai et al., 2014), and Portugal (Carvalho et al., 2015).
The origin of C_7 is not recent, with a tMRCA around 1983, with a most probable origin in Portugal, but the tMRCAs of its subclusters are relatively recent, in 2008 and 2004, respectively. C_7 is related to a virus from Mozambique, a former Portuguese colony where subtype C is the predominant HIV-1 genetic form (Bellocchi et al., 2005).
It should be noted that 3 individuals in each of clusters F1_3 and C_7, composed entirely of men, declared being heterosexual, and that in clusters CRF02_1 (excluding the heterosexual-associated CRF02_1_3 subcluster) and A1_1, among sexually infected individuals the number of self-declared heterosexual men exceeds the number of women (8 vs. 3 in CRF02_1 and 6 vs. 3 in A1_1). This suggests that at least some of the self-declared heterosexual men within the clusters could in fact be MSM, who, due to social stigma and discrimination, do not declare their real sexual behaviors, as suggested by other authors who found similar discrepancies between self-reported heterosexual behavior and phylogenetic clustering with sequences from MSM (Hué et al., 2014; Hoenigl et al., 2016; Ragonnet-Cronin et al., 2018).
The descriptions of the four clusters here analyzed is in line with those of other non-subtype B clusters reported to have expanded among MSM in Western Europe. However there are some salient features of the clusters here described that should be highlighted: first, the relatively large size of CRF02_1_1 and A1_1, greater than most non-subtype B clusters reported among MSM in Western Europe; second, their wide geographic distribution among different countries, which contrasts with the predominantly within-country clustering found by other authors in Europe (Frentz et al., 2013; Paraskevis et al., 2019); and, third, the rapid expansion of F1_3, with 36 diagnoses in only 4 years.
The expansion of large clusters among MSM in recent years, as here reported, reflects the existence of high risk sexual behaviors, which should prompt public health authorities to implement public health measures aimed at preventing HIV-1 transmission in this population, including behavioral interventions to reduce risky practices, preexposure prophylaxis (Grant et al., 2010; Volk et al., 2015; McCormack et al., 2016), and early diagnosis and treatment of HIV-1 infections (European Centre for Disease Prevention and Control, 2015; United Nations Population Fund, the Global Forum on MSM and HIV, United Nations Development Programme, World Health Organization, United States Agency for International Development, and World Bank, 2015). Prevention of HIV-1 transmission among MSM could also result in a reduction of heterosexually-transmitted infections, which can have their source in MSM networks, as seen in subcluster CRF02_1_3, associated with heterosexual transmission, and as reported by other authors (Oster et al., 2015; Esbjörnsson et al., 2016).
Continued HIV-1 molecular surveillance will be necessary to gain insight in real time on the dynamics of expansion of transmission networks, which will allow to focus prophylactic efforts in populations with the highest risk of HIV-1 acquisition and ongoing transmission (Little et al., 2014; Wang et al., 2015; Poon et al., 2016; Ratmann et al., 2016; Brenner et al., 2017; Chaillon et al., 2017; German et al., 2017; Oster et al., 2018; Wertheim et al., 2018) and to monitor the efficacy of public health interventions aimed at controlling the epidemic (Wertheim et al., 2011; Magiorkinis et al., 2018). HIV-1 molecular surveillance can also provide important information for the design of vaccine immunogens adapted to the major HIV-1 variants actively propagating in different areas, considering the correlation of HIV-1 clades to susceptibility to protective immune responses (Cao et al., 2000; Thomson et al., 2002; Binley et al., 2004; Geldmacher et al., 2007; Seaman et al., 2010; Hraber et al., 2014), and with potential to induce broadly neutralizing antibody responses (Kouyos et al., 2018), and will allow to obtain reagents derived from these variants for use in vaccine-related research (Cuevas et al., 2010; Revilla et al., 2011; Hora et al., 2016). These reagents will also be useful for studies on the biological basis of increased pathogenicity (Baeten et al., 2007; Kiwanuka et al., 2010; Li et al., 2014; Pérez-Álvarez et al., 2014; Kouri et al., 2015; Venner et al., 2016) and transmissibility (Kiwanuka et al., 2009) and diminished response to ARV drugs (Pernas et al., 2014; Cid-Silva et al., 2018) exhibited by some HIV-1 variants.
Statements
Data availability statement
The datasets generated for this study can be found in GenBank, MK177651–MK177824, MK177825–MK177829, KT276258, KY496622, and KY989952.
Ethics statement
This study was carried out in accordance with the recommendations of ’name of guidelines, name of committee’ with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Bioethics and Animal Well-being Committee of Instituto de Salud Carlos III, Majadahonda, Madrid, Spain.
Author contributions
MT, ED, and LP-Á conceived the study and supervised the experimental work. ED, MT, MC, AF-G, FD-F, JC, JM-L, and MS processed sequences and performed phylogenetic analyses. MT and FD-F performed phylodynamic analyses. HG performed data curation and phylogenetic analyses. SB, VM, AF-G, MS-M, EG-B, and CC performed experimental work. The members of the Spanish Group for the Study of New HIV Diagnoses recruited patients and obtained epidemiological data. MT wrote the manuscript with contributions from the other authors.
Funding
This work was funded through Acción Estratégica en Salud Intramural (AESI), Instituto de Salud Carlos III, project “Estudios sobre vigilancia epidemiológica molecular del VIH-1 en España,” PI16CIII/00033; Red de Investigación en SIDA (RIS), Instituto de Salud Carlos III, Subdirección General de Evaluación y Fondo Europeo de Desarrollo Regional (FEDER), Plan Nacional I+D+I, project RD16ISCIII/0002/0004; scientific agreements with Consellería de Sanidade, Government of Galicia (MVI 1004/16) and Osakidetza-Servicio Vasco de Salud, Government of Basque Country (MVI 1001/16); European Research Infrastructures for Poverty Related Diseases (EURIPRED). Seventh Framework Program: FP7-Capacities-infrastructures-2012-1, grant agreement 312661; and Dirección General de Farmacia, Ministerio de Sanidad, Servicios Sociales e Igualdad, Government of Spain (grant EC11-272).
Acknowledgments
We would like to thank José Antonio Taboada, from Consellería de Sanidade, Xunta de Galicia, and Daniel Zulaika, from Osakidetza-Servicio Vasco de Salud, for their support of this study, and the personnel at the Genomic Unit, Instituto de Salud Carlos III, for technical assistance in sequencing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00655/full#supplementary-material
Abbreviations
- MCC
maximum clade credibility
- MCMC
Markov chain Monte Carlo
- ML
maximum likelihood
- MRCA
most recent common ancestor
- MSM
men who have sex with men
- NFLG
near full-length genome
- PP
posterior probability
- PR-RT
protease-reverse transcriptase
- PWID
people who inject drugs
- tMRCA
time of most recent common ancestor
References
1
AbecasisA. B.MartinsA.CostaI.CarvalhoA. P.DiogoI.GomesP.et al (2011). Molecular epidemiological analysis of paired pol/env sequences from Portuguese HIV type 1 patients.AIDS Res. Hum. Retroviruses27803–805. 10.1089/AID.2010.0312
2
AbecasisA. B.WensingA. M.ParaskevisD.VercauterenJ.TheysK.van de VijverD. A.et al (2013). HIV-1 subtype distribution and its demographic determinants in newly diagnosed patients in Europe suggest highly compartmentalized epidemics.Retrovirology10:7. 10.1186/1742-4690-10-17
3
AmbrosioniJ.JunierT.DelhumeauC.CalmyA.HirschelB.ZdobnovE.et al (2012). Impact of highly active antiretroviral therapy on the molecular epidemiology of newly diagnosed HIV infections.AIDS262079–2086. 10.1097/QAD.0b013e32835805b6
4
AnisimovaM.GascuelO. (2006). Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative.Syst. Biol.55539–552. 10.1080/10635150600755453
5
ApetreiC.Loussert-AjakaI.CollinG.LetourneurF.DucaM.SaragostiS.et al (1997). HIV type 1 subtype F sequences in Romanian children and adults.AIDS Res. Hum. Retroviruses13363–365. 10.1089/aid.1997.13.363
6
AudelinA. M.CowanS. A.ObelN.NielsenC.JorgensenL. B.GerstoftJ. (2013). Phylogenetics of the Danish HIV epidemic: the role of very late presenters in sustaining the epidemic.J. Acquir. Immune Defic. Syndr.62102–108. 10.1097/QAI.0b013e318276becc
7
BaetenJ. M.ChohanB.LavreysL.ChohanV.McClellandR. S.CertainL.et al (2007). HIV-1 subtype D infection is associated with faster disease progression than subtype A in spite of similar plasma HIV-1 loads.J. Infect. Dis.1951177–1180. 10.1086/512682
8
BandeaC. I.RamosA.PieniazekD.PascuR.TanuriA.SchochetmanG.et al (1995). Epidemiologic and evolutionary relationships between Romanian and Brazilian HIV-subtype F strains.Emerg. Infect. Dis.191–93. 10.3201/eid0103.950305
9
BartoloI.AbecasisA. B.BorregoP.BarrosoH.McCutchanF.GomesP.et al (2011). Origin and epidemiological history of HIV-1 CRF14_BG.PLoS One6:e24130. 10.1371/journal.pone.0024130
10
BaryshevP. B.BogachevV. V.GashnikovaN. M. (2014). HIV-1 genetic diversity in Russia: CRF63_02A1, a new HIV type 1 genetic variant spreading in Siberia.AIDS Res. Hum. Retroviruses30592–597. 10.1089/AID.2013.0196
11
BelloG.AfonsoJ. M.MorgadoM. G. (2012). Phylodynamics of HIV-1 subtype F1 in Angola, Brazil and Romania.Infect. Genet. Evol.121079–1086. 10.1016/j.meegid.2012.03.014
12
BellocchiM. C.ForbiciF.PalombiL.GoriC.CoelhoE.SvicherV.et al (2005). Subtype analysis and mutations to antiviral drugs in HIV-1-infected patients from Mozambique before initiation of antiretroviral therapy: results from the DREAM programme.J. Med. Virol.76452–458. 10.1002/jmv.20382
13
BeloukasA.PsarrisA.GiannelouP.KostakiE.HatzakisA.ParaskevisD. (2016). Molecular epidemiology of HIV-1 infection in Europe: an overview.Infect. Genet. Evol.46180–189. 10.1016/j.meegid.2016.06.033
14
BeyrerC.BaralS. D.van GriensvenF.GoodreauS. M.ChariyalertsakS.WirtzA. L.et al (2012). Global epidemiology of HIV infection in men who have sex with men.Lancet380367–377. 10.1016/S0140-6736(12)60821-6
15
BezemerD.de WolfF.BoerlijstM. C.van SighemA.HollingsworthT. D.PrinsM.et al (2008). A resurgent HIV-1 epidemic among men who have sex with men in the era of potent antiretroviral therapy.AIDS221071–1077. 10.1097/QAD.0b013e3282fd167c
16
BezemerD.van SighemA.LukashovV. V.van der HoekL.BackN.SchuurmanR.et al (2010). Transmission networks of HIV-1 among men having sex with men in the Netherlands.AIDS24271–282. 10.1097/QAD.0b013e328333ddee
17
BinleyJ. M.WrinT.KorberB.ZwickM. B.WangM.ChappeyC.et al (2004). Comprehensive cross-clade neutralization analysis of a panel of anti-human immunodeficiency virus type 1 monoclonal antibodies.J. Virol.7813232–13252. 10.1128/JVI.78.23.13232-13252.2004
18
BobkovaM. (2013). Current status of HIV-1 diversity and drug resistance monitoring in the former USSR.AIDS Rev.15204–212.
19
BrachoM. A.SentandreuV.AlastruéI.BeldaJ.JuanA.Fernández-GarcíaE.et al (2014). Emerging trends in CRF02_AG variants transmission among men who have sex with men in Spain.J. Acquir. Immune Defic. Syndr.65e130–e133. 10.1097/01.qai.0000435602.73469.56
20
BrandD.MoreauA.CazeinF.LotF.PillonelJ.BrunetS.et al (2014). Characteristics of patients recently infected with HIV-1 non-B subtypes in France: a nested study within the mandatory notification system for new HIV diagnoses.J. Clin. Microbiol.524010–4016. 10.1128/JCM.01141-14
21
BrennerB. G.IbanescuR. I.HardyI.RogerM. (2017). Genotypic and phylogenetic insights on prevention of the spread of HIV-1 and drug resistance in “real-world” Settings.Viruses10:E10. 10.3390/v10010010
22
CaoH.ManiI.VincentR.MugerwaR.MugyenyiP.KankiP.et al (2000). Cellular immunity to human immunodeficiency virus type 1 (HIV-1) clades: relevance to HIV-1 vaccine trials in Uganda.J. Infect. Dis.1821350–1356. 10.1086/315868
23
CarrJ. K.NadaiY.EyzaguirreL.SaadM. D.KhakimovM. M.YakubovS. K.et al (2005). Outbreak of a West African recombinant of HIV-1 in Tashkent, Uzbekistan.J. Acquir. Immune Defic. Syndr.39570–575.
24
CarvalhoA.CostaP.TriunfanteV.BrancaF.RodriguesF.SantosC. L.et al (2015). Analysis of a local HIV-1 epidemic in Portugal highlights established transmission of non-B and non-G subtypes.J. Clin. Microbiol.531506–1514. 10.1128/JCM.03611-14
25
CasadoC.UrtasunI.SaragostiS.ChaixM. L.de RossiA.CattelanA. M.et al (2000). Different distribution of HIV type 1 genetic variants in European patients with distinct risk practices.AIDS Res. Hum. Retroviruses16299–304. 10.1089/088922200309403
26
CastroE.KhonkarlyM.CiuffredaD.BurgisserP.CavassiniM.YerlyS.et al (2010). HIV-1 drug resistance transmission networks in southwest Switzerland.AIDS Res. Hum. Retroviruses261233–1238. 10.1089/aid.2010.0083
27
ChaillonA.EssatA.FrangeP.SmithD. M.DelaugerreC.BarinF.et al (2017). Spatiotemporal dynamics of HIV-1 transmission in France (1999-2014) and impact of targeted prevention strategies.Retrovirology14:15. 10.1186/s12977-017-0339-4
28
ChalmetK.StaelensD.BlotS.DinakisS.PelgromJ.PlumJ.et al (2010). Epidemiological study of phylogenetic transmission clusters in a local HIV-1 epidemic reveals distinct differences between subtype B and non-B infections.BMC Infect. Dis.10:262. 10.1186/1471-2334-10-262
29
CiccozziM.GoriC.BorosS.Ruiz-AlvarezM. J.HarxhiA.DervishiM.et al (2005). Molecular diversity of HIV in Albania.J. Infect. Dis.192475–479. 10.1086/431599
30
Cid-SilvaP.Margusino-FraminanL.Balboa-BarreiroV.Martín-HerranzI.Castro-IglesiasA.Pernas-SoutoB.et al (2018). Initial treatment response among HIV subtype F infected patients who started antiretroviral therapy based on integrase inhibitors.AIDS32121–125. 10.1097/QAD.0000000000001679
31
CuevasM. T.Fernández-GarcíaA.PinillaM.García-ÁlvarezV.ThomsonM.DelgadoE.et al (2010). Short communication: biological and genetic characterization of HIV type 1 subtype B and nonsubtype B transmitted viruses: usefulness for vaccine candidate assessment.AIDS Res. Hum. Retroviruses261019–1025. 10.1089/aid.2010.0018
32
CuevasM. T.Muñoz-NietoM.ThomsonM. M.DelgadoE.IribarrenJ. A.CillaG.et al (2009). HIV-1 transmission cluster with T215D revertant mutation among newly diagnosed patients from the Basque Country, Spain.J. Acquir. Immune Defic. Syndr.5199–103. 10.1097/QAI.0b013e318199063e
33
DauweK.MortierV.SchauvliegeM.Van Den HeuvelA.FransenK.ServaisJ. Y.et al (2015). Characteristics and spread to the native population of HIV-1 non-B subtypes in two European countries with high migration rate.BMC Infect. Dis.15:524. 10.1186/s12879-015-1217-0
34
de OliveiraT.PillayD.GiffordR. J. (2010). The HIV-1 subtype C epidemic in South America is linked to the United Kingdom.PLoS One5:e9311. 10.1371/journal.pone.0009311
35
DelatorreE. O.BelloG.Eyer-SilvaW. A.Chequer-FernandezS. L.MorgadoM. G.Couto-FernandezJ. C. (2012). Evidence of multiple introductions and autochthonous transmission of the HIV type 1 CRF02_AG clade in Brazil.AIDS Res. Hum. Retroviruses281369–1372. 10.1089/aid.2011.0381
36
DelgadoE.CuevasM. T.DomínguezF.VegaY.CabelloM.Fernández-GarcíaA.et al (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 One10:e0143325. 10.1371/journal.pone.0143325
37
DelgadoE.Fernández-GarcíaA.VegaY.CuevasT.PinillaM.GarcíaV.et al (2012). Evaluation of genotypic tropism prediction tests compared with in vitro co-receptor usage in HIV-1 primary isolates of diverse subtypes.J. Antimicrob. Chemother.6725–31. 10.1093/jac/dkr438
38
DelgadoE.ThomsonM. M.VillahermosaM. L.SierraM.OcampoA.MirallesC.et al (2002). Identification of a newly characterized HIV-1 BG intersubtype circulating recombinant form in Galicia, Spain, which exhibits a pseudotype-like virion structure.J. Acquir. Immune Defic. Syndr.29536–543. 10.1097/00126334-200204150-00016
39
DrummondA. J.RambautA.ShapiroB.PybusO. G. (2005). Bayesian coalescent inference of past population dynamics from molecular sequences.Mol. Biol. Evol.221185–1192. 10.1093/molbev/msi103
40
DrummondA. J.SuchardM. A.XieD.RambautA. (2012). Bayesian phylogenetics with BEAUti and the BEAST 1.7.Mol. Biol. Evol.291969–1973. 10.1093/molbev/mss075
41
DumitrescuO.KalishM. L.KliksS. C.BandeaC. I.LevyJ. A. (1994). Characterization of human immunodeficiency virus type 1 isolates from children in Romania: identification of a new envelope subtype.J. Infect. Dis.169281–288. 10.1093/infdis/169.2.281
42
El MoussiA.ThomsonM. M.DelgadoE.CuevasM. T.NasrM.AbidS.et al (2017). Genetic diversity of HIV-1 in Tunisia.AIDS Res. Hum. Retroviruses3377–81. 10.1089/AID.2016.0164
43
EsbjörnssonJ.MildM.AudelinA.FonagerJ.SkarH.BruunJ. L.et al (2016). HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries.Virus Evol.2:vew010. 10.1093/ve/vew010
44
EstevesA.ParreiraR.PiedadeJ.VenennoT.FrancoM.Germanoet al (2003). Spreading of HIV-1 subtype G and envB/gagG recombinant strains among injecting drug users in Lisbon, Portugal.AIDS Res. Hum. Retroviruses19511–517. 10.1089/088922203766774568
45
EstevesA.ParreiraR.VenennoT.FrancoM.PiedadeJ.Germanoet al (2002). Molecular epidemiology of HIV type 1 infection in Portugal: high prevalence of non-B subtypes.AIDS Res. Hum. Retroviruses18313–325. 10.1089/088922202753519089
46
European Centre for Disease Prevention and Control (2015). ECDC Guidance: HIV and STI Prevention Among Men Who Have Sex with Men.Available at: https://ecdc.europa.eu/sites/portal/files/media/en/publications/Publications/hiv-sti-prevention-among-men-who-have-sex-with-men-guidance.pdf (accessed November 18 2018).
47
European Centre for Disease Prevention and Control (2017). HIV/AIDS Surveillance in Europe.Available at: https://ecdc.europa.eu/en/publications-data/hivaids-surveillance-europe-2017-2016-data (accessed November 18 2018).
48
EyzaguirreL. M.ErasilovaI. B.NadaiY.SaadM. D.KovtunenkoN. G.GomatosP. J.et al (2007). Genetic characterization of HIV-1 strains circulating in Kazakhstan.J. Acquir. Immune Defic. Syndr.4619–23. 10.1097/QAI.0b013e318073c620
49
FabeniL.AlteriC.OrchiN.GoriC.BertoliA.ForbiciF.et al (2015). Recent transmission clustering of HIV-1 C and CRF17_BF strains characterized by NNRTI-related mutations among newly diagnosed men in central Italy.PLoS One10:e0135325. 10.1371/journal.pone.0135325
50
FisherM.PaoD.BrownA. E.SudarshiD.GillO. N.CaneP.et al (2010). Determinants of HIV-1 transmission in men who have sex with men: a combined clinical, epidemiological and phylogenetic approach.AIDS241739–1747. 10.1097/QAD.0b013e32833ac9e6
51
FoleyB. T.LeitnerT.ParaskevisD.PeetersM. (2016). Primate immunodeficiency virus classification and nomenclature: review.Infect. Genet. Evol.46150–158. 10.1016/j.meegid.2016.10.018
52
FosterG. M.AmbroseJ. C.HuéS.DelpechV. C.FearnhillE.AbecasisA. B.et al (2014). Novel HIV-1 recombinants spreading across multiple risk groups in the United Kingdom: the identification and phylogeography of Circulating Recombinant Form (CRF) 50_A1D.PLoS One9:e83337. 10.1371/journal.pone.0083337
53
FrangeP.MeyerL.DeveauC.TranL.GoujardC.GhosnJ.et al (2012). Recent HIV-1 infection contributes to the viral diffusion over the French territory with a recent increasing frequency.PLoS One7:e31695. 10.1371/journal.pone.0031695
54
FransenK.BuveA.NkengasongJ. N.LagaM.van der GroenG. (1996). Longstanding presence in Belgians of multiple non-B HIV-1 subtypes.Lancet347:1403. 10.1016/S0140-6736(96)91042-9
55
FrentzD.WensingA. M.AlbertJ.ParaskevisD.AbecasisA. B.HamoudaO.et al (2013). Limited cross-border infections in patients newly diagnosed with HIV in Europe.Retrovirology10:36. 10.1186/1742-4690-10-36
56
GeldmacherC.CurrierJ. R.GerhardtM.HauleA.MabokoL.BirxD.et al (2007). In a mixed subtype epidemic, the HIV-1 Gag-specific T-cell response is biased towards the infecting subtype.AIDS21135–143. 10.1097/01.aids.0000247589.77061.f7
57
GermanD.GrabowskiM. K.BeyrerC. (2017). Enhanced use of phylogenetic data to inform public health approaches to HIV among men who have sex with men.Sex Health1489–96. 10.1071/SH16056
58
GiffordR. J.de OliveiraT.RambautA.PybusO. G.DunnD.VandammeA. M.et al (2007). Phylogenetic surveillance of viral genetic diversity and the evolving molecular epidemiology of human immunodeficiency virus type 1.J. Virol.8113050–13056. 10.1128/JVI.00889-07
59
GiffordR. J.LiuT. F.RheeS. Y.KiuchiM.HuéS.PillayD.et al (2009). The calibrated population resistance tool: standardized genotypic estimation of transmitted HIV-1 drug resistance.Bioinformatics251197–1198. 10.1093/bioinformatics/btp134
60
GiulianiM.SantoroM. M.LoP. A.CellaE.ScognamiglioP.LaiA.et al (2013). Circulation of HIV-1 CRF02_AG among MSM population in central Italy: a molecular epidemiology-based study.Biomed. Res. Int.2013:810617. 10.1155/2013/810617
61
González-DomenechC. M.VicianaI.DelayeL.MayorgaM. L.PalaciosR.de la TorreJ.et al (2018). Emergence as an outbreak of the HIV-1 CRF19_cpx variant in treatment-naive patients in southern Spain.PLoS One13:e0190544. 10.1371/journal.pone.0190544
62
GrantR. M.LamaJ. R.AndersonP. L.McMahanV.LiuA. Y.VargasL.et al (2010). Preexposure chemoprophylaxis for HIV prevention in men who have sex with men.N. Engl. J. Med.3632587–2599. 10.1056/NEJMoa1011205
63
GuimarãesM. L.VicenteA. C.OtsukiK.da SilvaR. F.FranciscoM.da SilvaF. G.et al (2009). Close phylogenetic relationship between Angolan and Romanian HIV-1 subtype F1 isolates.Retrovirology6:39. 10.1186/1742-4690-6-39
64
GuindonS.DufayardJ. F.LefortV.AnisimovaM.HordijkW.GascuelO. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0.Syst. Biol.59307–321. 10.1093/sysbio/syq010
65
GuindonS.LethiecF.DurouxP.GascuelO. (2005). PHYML Online–a web server for fast maximum likelihood-based phylogenetic inference.Nucleic Acids Res.33W557–W559. 10.1093/nar/gki352
66
HemelaarJ.ElangovanH.YunJ.Dickson-TettehL.FlemingerI.KirtleyS.et al (2018). Global and regional molecular epidemiology of HIV-1, 1990-2015: a systematic review, global survey, and trend analysis.Lancet Infect. Dis.19143–155. 10.1016/S1473-3099(18)30647-9
67
HIV Sequence Database (2019). HIV Sequence Database. Available at https://www.hiv.lanl.gov/content/sequence/HIV/mainpage.html (accessed February 15 2019).
68
HoeniglM.ChaillonA.KesslerH. H.HaasB.StelzlE.WeningerK.et al (2016). Characterization of HIV transmission in South-East Austria.PLoS One11:e0151478. 10.1371/journal.pone.0151478
69
HoraB.KeatingS. M.ChenY.SanchezA. M.SabinoE.HuntG.et al (2016). Genetic characterization of a panel of diverse HIV-1 isolates at seven international sites.PLoS One11:e0157340. 10.1371/journal.pone.0157340
70
HraberP.KorberB. T.LapedesA. S.BailerR. T.SeamanM. S.GaoH.et al (2014). Impact of clade, geography, and age of the epidemic on HIV-1 neutralization by antibodies.J. Virol.8812623–12643. 10.1128/JVI.01705-14
71
HuéS.BrownA. E.Ragonnet-CroninM.LycettS. J.DunnD. T.FearnhillE.et al (2014). Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions.AIDS281967–1975. 10.1097/QAD.0000000000000383
72
HuéS.PillayD.ClewleyJ. P.PybusO. G. (2005). Genetic analysis reveals the complex structure of HIV-1 transmission within defined risk groups.Proc. Natl. Acad. Sci. U.S.A.1024425–4429. 10.1073/pnas.0407534102
73
HughesG. J.FearnhillE.DunnD.LycettS. J.RambautA.Leigh BrownA. J. (2009). Molecular phylodynamics of the heterosexual HIV epidemic in the United Kingdom.PLoS Pathog.5:e1000590. 10.1371/journal.ppat.1000590
74
HusonD. H.ScornavaccaC. (2012). Dendroscope 3: an interactive tool for rooted phylogenetic trees and networks.Syst. Biol.611061–1067. 10.1093/sysbio/sys062
75
Instituto Nacional de Estadística (2018). Estadística del Padrón Continuo. Datos Provisionales a 1 de Enero de.Available at: http://www.ine.es/jaxi/Datos.htm?path=/t20/e245/p04/provi/l0/&file=00000010.px (accessed November 18 2018).
76
KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability.Mol. Biol. Evol.30772–780. 10.1093/molbev/mst010
77
KiwanukaN.LaeyendeckerO.QuinnT. C.WawerM. J.ShepherdJ.RobbM.et al (2009). HIV-1 subtypes and differences in heterosexual HIV transmission among HIV-discordant couples in Rakai, Uganda.AIDS232479–2484. 10.1097/QAD.0b013e328330cc08
78
KiwanukaN.RobbM.LaeyendeckerO.KigoziG.Wabwire-MangenF.MakumbiF. E.et al (2010). HIV-1 viral subtype differences in the rate of CD4+ T-cell decline among HIV seroincident antiretroviral naive persons in Rakai district, Uganda.J. Acquir. Immune Defic. Syndr.54180–184. 10.1097/QAI.0b013e3181c98fc0
79
KouriV.KhouriR.AlemánY.AbrahantesY.VercauterenJ.Pineda-PeñaA. C.et al (2015). CRF19_cpx is an evolutionary fit HIV-1 variant strongly associated with rapid progression to AIDS in Cuba.EBioMedicine2244–254. 10.1016/j.ebiom.2015.01.015
80
KouyosR. D.RusertP.KadelkaC.HuberM.MarzelA.EbnerH.et al (2018). Tracing HIV-1 strains that imprint broadly neutralizing antibody responses.Nature561406–410. 10.1038/s41586-018-0517-0
81
KuikenC.ThakallapalliR.EsklidA.de RondeA. (2000). Genetic analysis reveals epidemiologic patterns in the spread of human immunodeficiency virus.Am. J. Epidemiol.152814–822. 10.1093/aje/152.9.814
82
LaiA.BozziG.FranzettiM.BindaF.SimonettiF. R.de LucaA.et al (2016). HIV-1 A1 subtype epidemic in Italy originated from Africa and Eastern Europe and shows a high frequency of transmission chains involving intravenous drug users.PLoS One11:e0146097. 10.1371/journal.pone.0146097
83
LaiA.BozziG.FranzettiM.BindaF.SimonettiF. R.MicheliV.et al (2014). Phylogenetic analysis provides evidence of interactions between Italian heterosexual and South American homosexual males as the main source of national HIV-1 subtype C epidemics.J. Med. Virol.86729–736. 10.1002/jmv.23891
84
LaiA.SimonettiF. R.ZehenderG.de LucaA.MicheliV.MeravigliaP.et al (2012). HIV-1 subtype F1 epidemiological networks among Italian heterosexual males are associated with introduction events from South America.PLoS One7:e42223. 10.1371/journal.pone.0042223
85
LapovokI.KazennovaE.LagaV.VasilyevA.UtegenovaA.AbishevA.et al (2014). Short communication: molecular epidemiology of HIV type 1 infection in Kazakhstan: CRF02_AG prevalence is increasing in the southeastern provinces.AIDS Res. Hum. Retroviruses30769–774. 10.1089/AID.2013.0291
86
LefortV.LonguevilleJ. E.GascuelO. (2017). SMS: smart model selection in PhyML.Mol. Biol. Evol.342422–2424. 10.1093/molbev/msx149
87
LeozM.FeyertagF.CharpentierC.DelaugerreC.WirdenM.LemeeV.et al (2013). Characterization of CRF56_cpx, a new circulating B/CRF02/G recombinant form identified in MSM in France.AIDS272309–2312. 10.1097/QAD.0b013e3283632e0c
88
LewisF.HughesG. J.RambautA.PozniakA.Leigh BrownA. J. (2008). Episodic sexual transmission of HIV revealed by molecular phylodynamics.PLoS Med.5:e50. 10.1371/journal.pmed.0050050
89
LiY.HanY.XieJ.GuL.LiW.WangH.et al (2014). CRF01_AE subtype is associated with X4 tropism and fast HIV progression in Chinese patients infected through sexual transmission.AIDS28521–530. 10.1371/journal.pmed.0050050
90
LiitsolaK.RistolaM.HolmstromP.SalminenM.Brummer-KorvenkontioH.SimolaS.et al (2000). An outbreak of the circulating recombinant form AECM240 HIV-1 in the Finnish injection drug user population.AIDS142613–2615. 10.1097/00002030-200011100-00028
91
LittleS. J.Kosakovsky PondS. L.AndersonC. M.YoungJ. A.WertheimJ. O.MehtaS. R.et al (2014). Using HIV networks to inform real time prevention interventions.PLoS One9:e98443. 10.1371/journal.pone.0098443
92
LoleK. S.BollingerR. C.ParanjapeR. S.GadkariD.KulkarniS. S.NovakN. G.et al (1999). Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination.J. Virol.73152–160.
93
LouwagieJ.DelwartE. L.MullinsJ. I.McCutchanF. E.EddyG.BurkeD. S. (1994). Genetic analysis of HIV-1 isolates from Brazil reveals presence of two distinct genetic subtypes.AIDS Res. Hum. Retroviruses10561–567. 10.1089/aid.1994.10.561
94
LukashovV. V.KuikenC. L.VlahovD.CoutinhoR. A.GoudsmitJ. (1996). Evidence for HIV type 1 strains of U.S. intravenous drug users as founders of AIDS epidemic among intravenous drug users in northern Europe.AIDS Res. Hum. Retroviruses121179–1183. 10.1089/aid.1996.12.1179
95
MagiorkinisG.KaramitrosT.VasylyevaT. I.WilliamsL. D.MbisaJ. L.HatzakisA.et al (2018). An innovative study design to assess the community effect of interventions to mitigate HIV epidemics using transmission-chain phylodynamics.Am. J. Epidemiol.10.1093/aje/kwy160[Epub ahead of print].
96
McCormackS.DunnD. T.DesaiM.DollingD. I.GafosM.GilsonR.et al (2016). Pre-exposure prophylaxis to prevent the acquisition of HIV-1 infection (PROUD): effectiveness results from the pilot phase of a pragmatic open-label randomised trial.Lancet38753–60. 10.1016/S0140-6736(15)00056-2
97
MeloniS. T.KimB.SankaleJ. L.HamelD. J.TovanabutraS.MboupS.et al (2004). Distinct human immunodeficiency virus type 1 subtype A virus circulating in West Africa: sub-subtype A3.J. Virol.7812438–12445. 10.1128/JVI.78.22.12438-12445.2004
98
MonnoL.BrindicciG.LaiA.PunziG.AltamuraM.SimonettiF. R.et al (2012). An outbreak of HIV-1 BC recombinants in Southern Italy.J. Clin. Virol.55370–373. 10.1016/j.jcv.2012.08.014
99
MontavonC.Toure-KaneC.LiegeoisF.MpoudiE.BourgeoisA.VergneL.et al (2000). Most env and gag subtype A HIV-1 viruses circulating in West and West Central Africa are similar to the prototype AG recombinant virus IBNG.J. Acquir. Immune Defic. Syndr.23363–374. 10.1097/00126334-200004150-00001
100
MoskaleychikF. F.LagaV. Y.DelgadoE.VegaY.Fernández-GarcíaA.PérezA.et al (2015). [Rapid spread of the HIV-1 circular recombinant CRF02-AG in Russia and neighboring countries].Vopr. Virusol.6014–19.
101
NiculescuI.ParaschivS.ParaskevisD.AbagiuA.BatanI.BanicaL.et al (2015). Recent HIV-1 outbreak among intravenous drug users in Romania: evidence for cocirculation of CRF14_BG and subtype F1 Strains.AIDS Res. Hum. Retroviruses31488–495. 10.1089/aid.2014.0189
102
NúñezO.HernandoV.DíazA. (2018). Estimating the number of people living with HIV and the undiagnosed fraction in Spain in 2013.AIDS322573–2581. 10.1097/QAD.0000000000001989
103
Op de CoulE.van den BurgR.AsjöB.GoudsmitJ.CupsaA.et al (2000). Genetic evidence of multiple transmissions of HIV type 1 subtype F within Romania from adult blood donors to children.AIDS Res. Hum. Retroviruses16327–336. 10.1089/088922200309205
104
Op de CoulE. L.LukashovV. V.van DoornumG. J.GoudsmitJ.CoutinhoR. A. (1998). Multiple HIV-1 subtypes present amongst heterosexuals in Amsterdam 1988-1996: no evidence for spread of non-B subtypes.AIDS121253–1255. 10.1097/00002030-199810000-00024
105
OsterA. M.FranceA. M.PanneerN.Bañez OcfemiaM. C.CampbellE.DasguptaS.et al (2018). Identifying clusters of recent and rapid HIV transmission through analysis of molecular surveillance data.J. Acquir. Immune Defic. Syndr.79543–550. 10.1097/QAI.0000000000001856
106
OsterA. M.WertheimJ. O.HernandezA. L.Bañez OcfemiaM. C.SaduvalaN.HallH. I. (2015). Using molecular HIV surveillance data to understand transmission between subpopulations in the United States.J. Acquir. Immune Defic. Syndr.70444–451. 10.1097/QAI.0000000000000809
107
PalmaA. C.AraújoF.DuqueV.BorgesF.PaixãoM. T.CamachoR. (2007). Molecular epidemiology and prevalence of drug resistance-associated mutations in newly diagnosed HIV-1 patients in Portugal.Infect. Genet. Evol.7391–398. 10.1016/j.meegid.2007.01.009
108
PandeyS. S.CherianS.ThakarM.ParanjapeR. S. (2016). Short communication: phylogenetic and molecular characterization of six full-length HIV-1 genomes from India reveals a monophyletic lineage of Indian sub-subtype A1.AIDS Res. Hum. Retroviruses32489–502. 10.1089/AID.2015.0207
109
ParaskevisD.BeloukasA.StasinosK.PantazisN.de MendozaC.BannertN. (2019). HIV-1 molecular transmission clusters in nine European countries and Canada: association with demographic and clinical factors.BMC Med.17:4. 10.1186/s12916-018-1241-1
110
ParaskevisD.MagiorkinisE.MagiorkinisG.SypsaV.PaparizosV.LazanasM.et al (2007). Increasing prevalence of HIV-1 subtype A in Greece: estimating epidemic history and origin.J. Infect. Dis.1961167–1176. 10.1086/521677
111
ParczewskiM.Leszczyszyn-PynkaM.Witak-JedraM.SzetelaB.GasiorowskiJ.KnyszB.et al (2017). Expanding HIV-1 subtype B transmission networks among men who have sex with men in Poland.PLoS One12:e0172473. 10.1371/journal.pone.0172473
112
Patiño-GalindoJ. A.Torres-PuenteM.BrachoM. A.AlastruéI.JuanA.NavarroD.et al (2017). The molecular epidemiology of HIV-1 in the Comunidad Valenciana (Spain): analysis of transmission clusters.Sci. Rep.7:11584. 10.1038/s41598-017-10286-1
113
Patiño-GalindoJ. A.Torres-PuenteM.GimenoC.OrtegaE.NavarroD.GalindoM. J.et al (2015). Expansion of the CRF19_cpx variant in Spain.J. Clin. Virol.69146–149. 10.1016/j.jcv.2015.06.094
114
Pérez-ÁlvarezL.DelgadoE.VegaY.MonteroV.CuevasT.Fernández-GarcíaA.et al (2014). Predominance of CXCR4 tropism in HIV-1 CRF14_BG strains from newly diagnosed infections.J. Antimicrob. Chemother.69246–253. 10.1093/jac/dkt305
115
Pérez-ParraS.ÁlvarezM.Fernández-CaballeroJ. A.PérezA. B.SantosJ.BisbalO.et al (2018). Continued propagation of the CRF19_cpx variant among HIV-positive MSM patients in Spain.J. Antimicrob. Chemother.731031–1038. 10.1093/jac/dkx474
116
PernasB.GrandalM.MenaA.Castro-IglesiasA.CañizaresA.Wyleset al (2014). High prevalence of subtype F in newly diagnosed HIV-1 persons in northwest Spain and evidence for impaired treatment response.AIDS281837–1840. 10.1097/QAD.0000000000000326
117
Pineda-PeñaA. C.TheysK.StylianouD. C.DemetriadesI.AbecasisA. B.KostrikisL. G. (2018). HIV-1 Infection in Cyprus, the Eastern Mediterranean European frontier: a densely sampled transmission dynamics analysis from 1986 to 2012.Sci. Rep.8:1702. 10.1038/s41598-017-19080-5
118
PoonA. F.GustafsonR.DalyP.ZerrL.DemlowS. E.WongJ.et al (2016). Near real-time monitoring of HIV transmission hotspots from routine HIV genotyping: an implementation case study.Lancet HIV3e231–e238. 10.1016/S2352-3018(16)00046-1
119
PriceM. N.DehalP. S.ArkinA. P. (2010). FastTree 2–approximately maximum-likelihood trees for large alignments.PLoS One5:e9490. 10.1371/journal.pone.0009490
120
Ragonnet-CroninM.HuéS.HodcroftE. B.TostevinA.DunnD.FawcettT.et al (2018). Non-disclosed men who have sex with men in UK HIV transmission networks: phylogenetic analysis of surveillance data.Lancet HIV5e309–e316. 10.1016/S2352-3018(18)30062-6
121
Ragonnet-CroninM.LycettS. J.HodcroftE. B.HuéS.FearnhillE.BrownA. E.et al (2016). Transmission of non-B HIV subtypes in the United Kingdom is increasingly driven by large non-heterosexual transmission clusters.J. Infect. Dis.2131410–1418. 10.1093/infdis/jiv758
122
RambautA.LamT. T.MaxC. L.PybusO. G. (2016). Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen).Virus Evol.2:vew007. 10.1093/ve/vew007
123
RatmannO.van SighemA.BezemerD.GavryushkinaA.JurriaansS.WensingA.et al (2016). Sources of HIV infection among men having sex with men and implications for prevention.Sci. Transl. Med.8:320ra2. 10.1126/scitranslmed.aad1863
124
ReisM. N. GGuimarñesM. L.BelloG.StefaniM. M. A. (2019). Identification of new HIV-1 circulating recombinant forms CRF81_cpx and CRF99_BF1 in central western brazil and of unique BF1 recombinant forms.Front. Microbiol.10:97. 10.3389/fmicb.2019.00097
125
RevillaA.DelgadoE.ChristianE. C.DalrympleJ.VegaY.CarreraC.et al (2011). Construction and phenotypic characterization of HIV type 1 functional envelope clones of subtypes G and F.AIDS Res. Hum. Retroviruses27889–901. 10.1089/AID.2010.0177
126
SeamanM. S.JanesH.HawkinsN.GrandpreL. E.DevoyC.GiriA.et al (2010). Tiered categorization of a diverse panel of HIV-1 Env pseudoviruses for assessment of neutralizing antibodies.J. Virol.841439–1452. 10.1128/JVI.02108-09
127
ShapiroB.RambautA.DrummondA. J. (2006). Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences.Mol. Biol. Evol.237–9. 10.1093/molbev/msj021
128
ShcherbakovaN. S.ShalamovaL. A.DelgadoE.Fernández-GarcíaA.VegaY.KarpenkoL. I.et al (2014). Short communication: molecular epidemiology, phylogeny, and phylodynamics of CRF63_02A1, a recently originated HIV-1 circulating recombinant form spreading in Siberia.AIDS Res. Hum. Retroviruses30912–919. 10.1089/AID.2014.0075
129
SierraM.ThomsonM. M.RíosM.CasadoG.de CastroR. O.DelgadoE.et al (2005). The analysis of near full-length genome sequences of human immunodeficiency virus type 1 BF intersubtype recombinant viruses from Chile, Venezuela and Spain reveals their relationship to diverse lineages of recombinant viruses related to CRF12_BF.Infect. Genet. Evol.5209–217. 10.1016/j.meegid.2004.07.010
130
SimonettiF. R.LaiA.MonnoL.BindaF.BrindicciG.PunziG.et al (2014). Identification of a new HIV-1 BC circulating recombinant form (CRF60_BC) in Italian young men having sex with men.Infect. Genet. Evol.23176–181. 10.1016/j.meegid.2014.02.007
131
TamaletC.RavauxI.MoreauJ.BregigeonS.TourresC.RichetH.et al (2015). Emergence of clusters of CRF02_AG and B human immunodeficiency viral strains among men having sex with men exhibiting HIV primary infection in southeastern France.J. Med. Virol.871327–1333. 10.1002/jmv.24184
132
TemereancaA.EneL.MehtaS.ManolescuL.DuiculescuD.RutaS. (2013). Transmitted HIV drug resistance in treatment-naive Romanian patients.J. Med. Virol.851139–1147. 10.1002/jmv.23572
133
ThomsonM. M.DelgadoE.ManjónN.OcampoA.VillahermosaM. L.MariñoA.et al (2001). HIV-1 genetic diversity in Galicia Spain: BG intersubtype recombinant viruses circulating among injecting drug users.AIDS15509–516. 10.1097/00002030-200103090-00010
134
ThomsonM. M.Fernández-GarcíaA.DelgadoE.VegaY.Díez-FuertesF.Sánchez-MartínezM.et al (2012). 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.59e49–e51. 10.1097/QAI.0b013e3182400fc4
135
ThomsonM. M.NájeraR. (2001). Travel and the introduction of human immunodeficiency virus type 1 non-B subtype genetic forms into Western countries.Clin. Infect. Dis.321732–1737. 10.1086/320764
136
ThomsonM. M.Pérez-ÁlvarezL.NájeraR. (2002). Molecular epidemiology of HIV-1 genetic forms and its significance for vaccine development and therapy.Lancet Infect. Dis.2461–471. 10.1016/S1473-3099(02)00343-2
137
United Nations Population Fund the Global Forum on MSM and HIV United Nations Development Programme World Health Organization United States Agency for International Development and World Bank (2015). Implementing Comprehensive HIV and STI Programmes with Men Who Have Sex with Men.Available at: https://www.unfpa.org/es/node/13155
138
VegaY.DelgadoE.Fernández-GarcíaA.CuevasM. T.ThomsonM. M.MonteroV.et al (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 One10:e0125699. 10.1371/journal.pone.0125699
139
VennerC. M.NankyaI.KyeyuneF.DemersK.KwokC.ChenP. L.et al (2016). Infecting HIV-1 subtype predicts disease progression in women of sub-Saharan Africa.EBioMedicine13305–314. 10.1016/j.ebiom.2016.10.014
140
VerhofstedeC.DauweK.FransenK.Van LaethemK.Van den WijngaertS.RuelleJ.et al (2018). Phylogenetic analysis of the Belgian HIV-1 epidemic reveals that local transmission is almost exclusively driven by men having sex with men despite presence of large African migrant communities.Infect. Genet. Evol.6136–44. 10.1016/j.meegid.2018.03.002
141
VinkenL.FransenK.Pineda-PeñaA. C.AlexievI.BalottaC.DebaisieuxL.et al (2017). HIV-1 sub-subtype F1 outbreak among MSM in Belgium.Virus Evol.3(Suppl. 1):vew036.020. 10.1093/ve/vew036.020
142
VolkJ. E.MarcusJ. L.PhengrasamyT.BlechingerD.NguyenD. P.FollansbeeS.et al (2015). No new HIV infections with increasing use of HIV preexposure prophylaxis in a clinical practice setting.Clin. Infect. Dis.611601–1603. 10.1093/cid/civ778
143
von WylV.KouyosR. D.YerlyS.BoniJ.ShahC.BurgisserP.et al (2011). The role of migration and domestic transmission in the spread of HIV-1 non-B subtypes in Switzerland.J. Infect. Dis.2041095–1103. 10.1093/infdis/jir491
144
WangX.WuY.MaoL.XiaW.ZhangW.DaiL.et al (2015). Targeting HIV prevention based on molecular epidemiology among deeply sampled subnetworks of men who have sex with men.Clin. Infect. Dis.611462–1468. 10.1093/cid/civ526
145
WertheimJ. O.Kosakovsky PondS. L.LittleS. J.De GruttolaV. (2011). Using HIV transmission networks to investigate community effects in HIV prevention trials.PLoS One6:e27775. 10.1371/journal.pone.0027775
146
WertheimJ. O.MurrellB.MehtaS. R.ForgioneL. A.Kosakovsky PondS. L.SmithD. M.et al (2018). Growth of HIV-1 molecular transmission clusters in New York City.J. Infect. Dis.2181943–1953. 10.1093/infdis/jiy431
147
YebraG.de MulderM.MartínL.RodríguezC.LabargaP.VicianaI.et al (2012). Most HIV type 1 non-B infections in the Spanish cohort of antiretroviral treatment-naive HIV-infected patients (CoRIS) are due to recombinant viruses.J. Clin. Microbiol.50407–413. 10.1128/JCM.05798-11
148
YebraG.FramptonD.Gallo CassarinoT.RaffleJ.HubbJ.FernsR. B.et al (2018). A high HIV-1 strain variability in London, UK, revealed by full-genome analysis: results from the ICONIC project.PLoS One13:e0192081. 10.1371/journal.pone.0192081
149
ZehenderG.EbranatiE.LaiA.SantoroM. M.AlteriC.GiulianiM.et al (2010). Population dynamics of HIV-1 subtype B in a cohort of men-having-sex-with-men in Rome, Italy.J. Acquir. Immune Defic. Syndr.55156–160. 10.1097/QAI.0b013e3181eb3002
Summary
Keywords
HIV-1, molecular epidemiology, phylogeny, phylodynamics, men who have sex with men, subtypes, circulating recombinant forms, clusters
Citation
Delgado E, Benito S, Montero V, Cuevas MT, Fernández-García A, Sánchez-Martínez M, García-Bodas E, Díez-Fuertes F, Gil H, Cañada J, Carrera C, Martínez-López J, Sintes M, Pérez-Álvarez L, Thomson MM and The Spanish Group for the Study of New HIV Diagnoses (2019) Diverse Large HIV-1 Non-subtype B Clusters Are Spreading Among Men Who Have Sex With Men in Spain. Front. Microbiol. 10:655. doi: 10.3389/fmicb.2019.00655
Received
25 November 2018
Accepted
15 March 2019
Published
03 April 2019
Volume
10 - 2019
Edited by
Gkikas Magiorkinis, National and Kapodistrian University of Athens, Greece
Reviewed by
Apostolos Beloukas, University of West Attica, Greece; Juan Angel Patino Galindo, Columbia University, United States
Updates
Copyright
© 2019 Delgado, Benito, Montero, Cuevas, Fernández-García, Sánchez-Martínez, García-Bodas, Díez-Fuertes, Gil, Cañada, Carrera, Martínez-López, Sintes, Pérez-Álvarez, Thomson and the Spanish Group for the Study of New HIV Diagnoses.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael M. Thomson, mthomson@isciii.es
This article was submitted to Virology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.