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

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

Introduction

The HIV-1 epidemic is still one of the major public health problems in Spain. Around 3,500–4,000 new diagnoses (NDs) of HIV-1 infection are reported every year, with an estimated incidence of NDs of 7.5 per 100,000 population in 2019 (HIV STI and Hepatitis Surveillance Unit, 2020). A decreasing trend in the incidence of NDs has been observed since 2010, although it is still higher than the average rate found in the EU/EEA (European Centre for Disease Prevention and Control and WHO Regional Office for Europe, 2020). In 2019, most (56%) of the NDs were diagnosed among men who had sex with men (MSM) and 36% of all reported NDs were in people born outside of Spain (HIV STI and Hepatitis Surveillance Unit, 2020).

Molecular epidemiology is an important tool for describing the HIV-1 epidemic (Brenner et al., 2013; Paraskevis et al., 2016; Wertheim et al., 2017; Board et al., 2020). Individuals whose transmission is related group phylogenetically in clades named transmission clusters (TCs). Due to the high genetic variability of HIV-1, phylogenetic analysis allows to reconstruct transmission events through the identified TCs and infer the history of the HIV-1 epidemic (Hué et al., 2004, 2005).

National databases of protease and reverse transcriptase (Pr-RT) sequences, primarily obtained for antiretroviral drug resistance testing, contain valuable data about HIV-1 expansion and have been used in molecular epidemiology studies (Ragonnet-Cronin et al., 2016a; Parczewski et al., 2017; Oster et al., 2018; Petersen et al., 2018; Verhofstede et al., 2018; Pineda-Peña et al., 2019; Fabeni et al., 2021). The phylogenetic analyses performed in these studies combined with clinical and epidemiological data of the patients provide relevant public health information for the implementation of control measures and for monitoring the HIV-1 epidemic (Brenner et al., 2013; Vasylyeva et al., 2016, 2020; Paraskevis et al., 2019; Board et al., 2020; Campbell et al., 2020).

Spanish clinical guidelines recommend to perform a genotypic drug resistance test before starting antiretroviral therapy (ART) in all HIV-1-diagnosed patients and in ART-failing cases [AIDS study group (GESIDA) of the Spanish Society of Infectious Diseases and Clinical Microbiology and the National AIDS Plan, 2020]. Pr-RT sequences obtained for these tests have been analyzed in different molecular epidemiology studies, describing the genetic features of the HIV-1 epidemic in different regions of Spain (Thomson et al., 2001; Holguín et al., 2007; Cuevas et al., 2009a,b; González-Alba et al., 2011; Yebra et al., 2013; Pérez-Parra et al., 2015, 2016; Patiño-Galindo et al., 2016, 2017b; González-Domenech et al., 2020). The phylogenetic studies have also allowed the identification of large TCs among MSM which are actively growing in Spain (Delgado et al., 2015, 2019; Patiño-Galindo et al., 2017a; González-Domenech et al., 2018), as well as transmitted drug resistance mutations which are spreading in TCs (Cuevas et al., 2009b; Vega et al., 2015; Viciana et al., 2016; González-Domenech et al., 2018, 2020), highlighting the value of molecular epidemiology as a tool for HIV surveillance.

The Basque Country and Galicia are two regions located in northern Spain, comparable in terms of population and HIV-1 diagnosis rates. During 2013–2018, the mean resident population and rates of new HIV diagnoses were estimated in 2,2 million and 6,9 per 100,000 inhabitants for the Basque Country and 2,7 million and 5,6 per 100,000 inhabitants for Galicia (HIV STI and Hepatitis Surveillance Unit, 2020; Instituto Nacional de Estadistica, 2021). In this study, we analyzed the frequency of grouping in TCs among patients diagnosed in 2013–2018 in these northern Spanish regions, identifying the different features associated with TCs in the HIV-1 epidemics in these regions.

Materials and Methods

Patients

Samples from individuals newly diagnosed of HIV-1 infection in the northern Spanish regions of Galicia and Basque Country during 2013–2018, sent to the HIV Biology and Variability Unit, Centro Nacional de Microbiología, Instituto de Salud Carlos III, were included in the study.

The representativeness of the NDs of our cohort was estimated comparing them to the NDs reported to the Spanish information system on new HIV diagnoses (SINIVIH by its Spanish acronym; HIV STI and Hepatitis Surveillance Unit, 2020), from the same period and studied regions. As patients in both databases cannot be linked, gender, age group, transmission route and region of birth frequencies were compared between both groups of individuals using the Chi-squared test.

The use of anonymized, de-identified clinical/demographic and sequence data was reviewed and approved under an exempt protocol by the Bioethics and Animal Well-being Committee of Instituto de Salud Carlos III, with report numbers CEI PI 38_2016-v3 (dated 20 June 2016) and CEI PI 31_2019-v5 (dated 6 November 2019). This study did not require written informed consent by study participants, except for those participants who required additional samples different from the ones obtained in the routine clinical practice.

Nucleic Acid Extraction, Amplification, and Sequencing

Nucleic acid was extracted from plasma or whole blood samples. RNA was extracted from 1 ml plasma using NUCLISENS® easyMAG® (BioMérieux, Marcy l’Etoile, France) and DNA was extracted from 200 μl whole blood using QIAamp® DNA DSP blood mini kit (Qiagen, Hilden, Germany), following the manufacturer’s instructions. A Pr-RT fragment of pol (HXB2 positions 2253–3629) was amplified by RT-PCR followed by nested PCR from RNA or by nested PCR from DNA. Reagents, PCR thermal profiles, and primers are described in Supplementary Tables 1, 2.

Population sequencing was performed with ABI Prism BigDye Terminator Cycle Sequencing kit and ABI 3730 XL sequencer (Applied Biosystems, Foster City, CA, United States) in the Genomic Unit of Instituto de Salud Carlos III. Sequences were assembled with SeqMan Pro v.12.2.1 (DNA STAR Lasergene, Madison, WI, United States) and edited with BioEdit v.7.2.5 (Hall, 1999).

Phylogenetic Analyses and Transmission Cluster Identification

Partial pol sequences from the patients included in this study and sequences obtained by us from another 10,687 HIV-1-infected individuals attended in clinical centers from Spain whose samples were collected in the period 1999–2019 were included in the analyses. A single sequence per patient was used in the phylogenetic analysis. In cases where several sequences were available for one patient, that closer to the date of HIV-1 diagnosis was selected.

Reference sequences from the different subtypes and CRFs retrieved from the Los Alamos HIV Sequence Database¹ were included for the phylogenetic and bootscanning analysis. In addition, we conducted BLAST (Altschul et al., 1990) searches of each of the sequences obtained by us against the Los Alamos database, including in the phylogenetic analysis up to 10 most similar sequences from each search.

Sequences were analyzed phylogenetically by an approximately maximum-likelihood method using FastTree2 (Price et al., 2010), running it in a local desktop computer. In these analyses, the general time reversible model of nucleotide substitution with CAT approximation to account for among-site heterogeneity in substitution rates was used, and the reliability of nodes was assessed with Shimodaira–Hasegawa (SH)-like local support values (Shimodaira and Hasegawa, 1999; Guindon et al., 2010). Classification of sequences in subtypes and recombinant forms was based on clustering with clade references in approximately maximum-likelihood trees. Sequences suspected of intersubtype recombination were subsequently analyzed by bootscanning with SimPlot v3.5 (Lole et al., 1999).

Transmission clusters were defined as those comprising viruses from four or more individuals, at least 50% of them Spanish, and whose sequences grouped in the phylogenetic tree with a SH-like node support value ≥ 0.95. This limit in the percentage of migrants included in TCs was used to minimize the effect of viruses which are circulating in the countries of origin of the migrants, grouping phylogenetically in TCs in our sequence dataset due to the relatedness of the viruses of those countries, although they are not actually spreading in Spain.

Statistical Analysis

To evaluate the association between variables, Chi-squared and Fisher’s exact tests were used, with associations being considered statistically significant at a value of p < 0.05. A multivariable logistic regression model was performed to identify factors associated to TCs. The model was adjusted by gender, age group, transmission route, country of origin of the patient, Spanish region of sample collection and HIV-1 genetic form. Associations were measured using the odds ratio (OR) and its 95% CI. Data analyses were performed using the STATA statistical software package Version 16 (Stata Corporation, College Station, TX, United States).

Sequence Accession Numbers

To avoid the identification of transmission networks and the potential breech of patient confidentiality, only 19% of the sequences have been deposited in GenBank, including at least two sequences from each TC comprising five or more patients with new HIV-1 diagnoses in the period of study. The sequences are under the following accession numbers: KT276255, KT276260, KT276263-KT276264, KT276266-KT276267, KT276270-KT276271, KU685562-KU685565, KU685567, KU685569-KU685575, KU685577, KU685581-KU685582, KU685586, KU685588-KU685590, KU685592, KX534325, KX534329, KX534331, KY465968, KY496624, KY514084, KY989950, MF999250-MF999256, MF999258-MF999261, MK177721-MK177733, MK177735-MK177752, MK177754-MK177757, MK177761-MK177772, MK177783-MK17785, MK177787-MK177790, MK177792, MK177795-MK177796, MK177798-MK177799, MK177801-MK177803, MK177805-MK177807, MK298150, MT436242, MT436244, MT436246-MT436247, MT436249-MT436250, MW344920-MW344921, MW584217-MW584224, OK011532, OK011542, OK011545-OK011546, OK011549, OK148912, OK148914, OK148917-OK148919, OK148921, OK148941-OK148942, OK148953-OK148957, OK148961, OK148965, OK148968-OK148972, OK148974, OL982314-OL982315, and OM914651-OM914711.

Results

Study Population

A total of 1,158 HIV-1 NDs from Galicia and Basque Country diagnosed in 2013–2018 were included in the study (Table 1). The majority of the patients were male (82%), with MSM being the most common transmission route (46%). Heterosexuals represented 34% and 15% of the individuals were male whose transmission route was reported as sexual without specifying whether they were in the MSM or heterosexual category. Subtype B (63%) was the most frequent genetic form among the individuals included in the study, followed by unique recombinant forms (URFs; 8.2%), CRF02_AG (7.7%), and subtype F (7.6%; Table 1).

TABLE 1

Table 1. Characteristics of the HIV-1 newly diagnosed individuals included in the study.

Seventy percent of the patients were native Spanish (Table 1), followed by Latin Americans (17%) and Sub-Saharan Africans (8.0%). There were important differences in gender and transmission route proportions between these two migrant populations, with female [52% (44/85) vs. 26% (37/180)] and heterosexual transmission [88% (65/74) vs. 33% (54/165)] being higher in Sub-Saharan Africans than in Latin Americans.

The distribution of patients by gender and transmission route was similar in both studied regions (Table 1), although the frequency of people who inject drugs (PWID) was slightly higher in Galicia (5.3 vs. 3.2%), but the difference was not statistically significant (p = 0.13). The percentage of patients of non-Spanish origin was 35% in the Basque Country (13% corresponding to Africans), a proportion which is statistically higher (p < 0.001) compared to the 17% found in Galicia (2.2% Africans). Statistical differences in the distribution of genetic forms were also observed between both regions (p < 0.001), with the frequencies of subtype F (13%) and CRF02_AG (9.9%) being higher in Galicia and the Basque Country, respectively (Table 1).

Representativeness of the Patients Included in the Study

Our cohort of newly diagnosed HIV-1-infected patients represents 64% of the notified cases in 2013–2018 (Supplementary Table 3), assuming that all NDs included in our study have been reported to the SINIVIH. This percentage was higher in the Basque Country, where it reached 86%, than in Galicia, where it was 42% of the notified cases. No statistical differences were found between both patient groups regarding the distribution of individuals by gender, transmission route, region of origin, or age group (Supplementary Table 3). A map with the studied regions, with the numbers of sequences included in our dataset and the representativeness in each region of the NDs analyzed in this study is shown in Supplementary Figure 1.

TCs Identified in the Study

Fifty-one percent (594/1158) of the individuals belonged to 162 different TCs. Regarding regions, 45% (351/776) of the patients were included in 98 TCs in the Basque Country and 64% (243/382) belonged to 82 TCs in Galicia (Table 2). Eighteen (11%) TCs comprised patients from both regions. The sizes of the TCs ranged from 4 to 205 patients. MSM was the main transmission route associated with 59% of the TCs (Table 3). Transmission routes other than MSM were more frequently associated with TCs in Galicia (44%) than in the Basque County (35%; Table 2), although the difference was not statistically significant (p = 0.207).

TABLE 2

Table 2. Characteristic of the identified transmission clusters (TCs) of the study.

TABLE 3

Table 3. Characteristics of the transmission clusters comprising ≥30 patients in the Basque Country and Galicia (2013–2018).

Non-subtype B genetic forms were found in 33 (21%) of the 162 TCs identified in the study. Sixteen of them derived from genetic forms circulating in Latin America (BF1 recombinants, Brazilian subtype C and F1 variants, CRF19_cpx and CRF20_BG) and 10 from genetic forms circulating in Africa (CRF02_AG, CRF02/A3, CRF02/G, CRF06_cpx and subtype C strains of African origin). Twenty-nine percent (138/474) of the Spaniards belonging to TCs were members of non-subtype B TCs.

A total of 15 large TCs (here defined as those with ≥30 individuals) were identified in the study (Table 3), 10 of them of subtype B. MSM was the associated transmission route in 14 (93%) of them. Seventy-six percent (67/88) of subtype F infections belonged to the MSM-associated TCs F1_1 and F1_3. Most of these large TCs have spread mainly in a specific region. Thus, TCs B07, B08, B09, B10, B12, B70, and F1_3 were found mainly in the Basque Country, while TCs F1_1, B13, BG_2, and B05 were spread mainly in Galicia (Table 3). In addition, A1_1, CRF02_1, B50 and B31 have a global expansion in Spain, with a high number of patients outside of both studied regions. Four TCs had greater than 50% increase in NDs during the period 2013–2018: F1_3 (94%), A1_1 (76%), B70 (72%), and CRF02_1 (53%; Table 3).

A phylogenetic tree which includes all the TCs comprising five or more patients newly diagnosed in the studied period in Galicia and Basque Country is shown in Supplementary Figure 2.

Factors Associated With TCs

The percentage of patients in TCs was higher in males (57%) than in females (22%), and in MSM (67%) than in heterosexuals (31%). This percentage was also higher in Spaniards (64%) than in Latin Americans (35%) or Sub-Saharan Africans (1.2%), and in individuals infected with viruses of subtypes F (76%) or B (58%) than in those infected with CRF02_AG viruses (17%). The different distribution of TCs according to gender, mode of HIV-1 transmission, country of origin of the patient, and HIV-1 genetic form were statistically significant (Table 4). The age group was also associated to TCs (p = 0.036). The age range of 20–29 years showed an increased percentage of patients in TCs (59%), which was especially high in Galicia (81%, p = 0.003) but was not statistical significant in the Basque Country, and no statistical differences were found in the distribution of main variables.

TABLE 4

Table 4. Characteristics of the new HIV diagnoses included in TCs (2013–2018).

Along the studied period, there was no clear trend in the proportion of patients associated with TCs (Supplementary Figure 3). In Galicia, a decrease of notified HIV-1 NDs (Supplementary Figure 3, panel A) with a slight reduction of the number of TCs and an increase in the ratio of patients per TC (Supplementary Figure 3, panel D) was observed during the study period.

In the multivariate analysis, males (OR: 2.6; 95% CI: 1.5–4.7) and MSM (OR: 2.1; 95% CI: 1.4–3.2) had higher odds of belonging to a TC compared to females and heterosexuals, respectively (Table 5). Patients from Galicia (OR: 1.5; 95% CI: 1.1–2.2) and subtype F infections (OR: 2.5; 95% CI: 1.3–5.0) were also factors associated with TCs (Table 5).

TABLE 5

Table 5. Factors associated with transmission clusters, univariate/multivariate analysis.

Individuals from Latin America (OR: 0.3; 95% CI: 0.2–0.4), North Africa (OR: 0.4; 95% CI: 0.2–1.0) and especially from Sub-Saharan Africa (OR: 0.02; 95% CI: 0.003–0.2) were inversely associated to belonging to TCs compared to Spanish patients (Table 5). Similarly, patients over 29 years old (OR: 0.6; 95% CI: 0.4–0.9) and infections with a virus of a CRF_BF (OR: 0.1; 95% CI: 0.03–0.4) or of a genetic form classified in the “other” category (OR: 0.3; 95% CI: 0.2–0.8) had lower odds of belonging to TCs than patients in the age range of 20–29 years and with subtype B infections, respectively (Table 5). These associations, although with slightly lower OR, where also observed when the criterion on the country of origin of the patients was excluded from the TC definition (Supplementary Table 4).

Discussion

The phylogenetic analysis of sequences from almost 12,000 HIV-1-infected patients from our database has allowed us to investigate the transmission events occurring during the 2013–2018 period in two regions of northern Spain, Galicia and Basque Country, identifying TCs and their associated factors. Our ND cohort had a good representativeness, corresponding to 64% of the NDs reported to the SINIVIH from these regions during that period, especially in the Basque Country where it reached 86%.

Our analysis has assigned 51% of the NDs from Galicia and Basque Country to 162 different TCs, indicating that TCs are playing an important role in the spread of HIV-1 infections in both regions. There is no consensus on the criteria for defining a TC, as these should depend on the aim of the study (Hassan et al., 2017). We considered only TCs including at least 50% Spaniards and at least four documented infections, because they may be more epidemiologically relevant to describe HIV-1 strains spreading in an area. Smaller clusters comprising two (“pairs”) or three (“triplets”) individuals are much more numerous than larger ones (Hughes et al., 2009; Brown et al., 2011; Yebra et al., 2013) and in some studies are distinguished from larger clusters (Hoenigl et al., 2016; Fabeni et al., 2020). They could represent sexual couples with or without a second relationship by one of the members of the couple, and their epidemiologic relevance is uncertain. The inclusion of these requirements in the TC definition has allowed us to better identify the structure and factors associated to HIV-1 strains whose transmission is established in Spain.

Transmission clusters were found associated with male gender and MSM. In fact, almost all the largest TCs (≥30 patients) identified in our study were associated with MSM. HIV-1 molecular epidemiological studies in different countries, independently of the methodology used, have found a similar association (Parczewski et al., 2017; Wertheim et al., 2017; Ramirez et al., 2021). Also, studies in other Spanish regions have found this association of TCs with MSM (Pérez-Parra et al., 2015; González-Domenech et al., 2020). These findings indicate that MSM TCs are playing a main role in spreading HIV-1 infections and are consistent with the fact that MSM is the main frequent mode of HIV-1 transmission in Spain since the early 2000s (Nuñez et al., 2018).

In Spain, 36% of HIV-1 infections diagnosed in 2019 were in migrants (HIV STI and Hepatitis Surveillance Unit, 2020). Migrants play an important role in HIV-1 epidemic spread, contributing to bridging between networks from different geographic areas (Grossman et al., 2015; Su et al., 2018; Grabowski et al., 2020). In this study, we have examined the frequency of membership in Spanish TCs of migrants from different geographic areas. We have found that Latin Americans and Africans showed lower odds of belonging to TCs than Spaniards. However, there were important differences between these migrant populations, as 35% of Latin Americans were included in TCs, compared to only 1.2% of Sub-Saharan Africans, suggesting that HIV-1 transmission between Latin Americans and Spaniards is relatively common, presumably derived from cultural and linguistic affinities (Osorno-González de León et al., 2021). Similarly to our results, a low percentage of Sub-Saharan Africans in TCs have been also reported in Italy (Fabeni et al., 2020), Belgium (Verhofstede et al., 2018), and Madrid, Spain (Yebra et al., 2013), while in Israel, cross-ethnic HIV-1 spread from Israeli-born individuals to Ethiopian-born immigrants was rare (Grossman et al., 2015). In our study, TCs are mainly MSM-driven. In contrast, we have found high frequencies of females (52%) and heterosexuals (88%) in Sub-Saharan Africans, which, in addition to the likely HIV-1 acquisition in their countries of origin, can explain the low proportion of individuals in this migrant population belonging to TCs. Indeed, we have observed that Sub-Saharan Africans are infected frequently with genetic forms circulating in Africa.

Determining the origin of HIV-1 infections in the migrant population would allow proper allocation of resources for preventive measures (Fakoya et al., 2015; Álvarez del Arco et al., 2017). Previous studies based on models using clinical and behavioral data estimated high levels of post-migration HIV acquisition in Europe (Desgrées-du-Loâ et al., 2015; Álvarez del Arco et al., 2017), up to 71% among migrants in Spain (Álvarez del Arco et al., 2017). This figure could be an overestimate due to biases in data collection and failure to consider infections acquired by migrants when visiting their countries, as migrants’ mobility is associated with increased risk for HIV acquisition (Fenton et al., 2001; Dias et al., 2020). Phylogenetic analyses can provide additional information on the geographical area of HIV-1 acquisition. In our study, the lower percentage of foreigners found in TCs could be due to infections with virus strains which are not circulating in Spain, suggestive of imported infections, which would be especially frequent in Sub-Saharan Africans. Further studies using a combined phylogenetic and epidemiological approach can address the accuracy of the origin of the HIV-1 infections in the migrant populations in Spain.

With regard to the role of immigration on the origin of HIV-1 TCs in Spain, it is also interesting to note that 21% TCs identified in this study are of non-subtype B genetic forms, with their ancestry showing some correlation with the predominant geographic origins of the immigrant population in Spain. Thus, the fact that 16 TCs derive from genetic forms circulating in Latin America and Africa may correlate to the high number of immigrants from these geographic areas living in Spain, representing around 36 and 18%, respectively, of the approximately 7.3 million immigrants (Instituto Nacional de Estadistica, 2021). Similarly, the expansion of a large F1 subsubtype cluster of Romanian origin (F1_3) in Spain (Delgado et al., 2019) may relate to the fact that among European immigrants in Spain, Romanians are the most numerous, representing around 8% of all immigrants (Instituto Nacional de Estadistica, 2021).

Subtype B was the most represented genetic form associated to TCs. Only subtype F infections (76% of NDs in TCs) showed a higher frequency than subtype B infections (56% of NDs in TCs) of belonging to TCs. This was mainly due to the presence of two large TCs of F1 subsubtype, designated F1_1 (205 individuals) and F1_3 (36 individuals), which have spread among MSM mainly in Galicia and Basque Country, respectively (Delgado et al., 2015, 2019). Interestingly, two other large non-subtype B TCs associated with MSM, of CRF02_AG and A1 subsubtype, are also present in these regions (Delgado et al., 2019). In Western Europe, the HIV-1 epidemic among MSM is dominated by subtype B. However, the epidemic among MSM in Spain is becoming increasingly diverse through the expansion of multiple non-subtype B TCs, comprising or related to viruses circulating in other countries (González-Domenech et al., 2018; Delgado et al., 2019). This phenomenon has been also documented in other European countries (Dauwe et al., 2015; Ragonnet-Cronin et al., 2016b; Verhofstede et al., 2018; Fabeni et al., 2019; Ramirez et al., 2021), where viruses introduced from abroad have expanded successfully among the MSM population.

In general, there was not a clear temporal trend in the proportion of patients belonging to TC, although in Galicia our data show an increase in the proportion of patients in TCs in the last years. However, trends are difficult to evaluate due to the presence of fast-growing TCs in some periods.

A different distribution of TCs is responsible for the spread of HIV-1 in Galicia and the Basque Country. Large MSM TCs have spread locally in both regions and the majority of the identified TCs were associated with this population group. However, we found a higher frequency of TCs associated with PWID or heterosexual transmission in Galicia. Galician patients have 1.5 higher odds of belonging to TCs than individuals from the Basque Country. This can be related to the low percentage of migrants in Galicia, who have a reduced association to TCs, and the high number of patients belonging to the large F1_1 TC, currently comprising 205 individuals, which has successfully spread in Galicia and other Spanish regions (Delgado et al., 2015), even in other European countries (Delgado et al., 2015; Vinken et al., 2019). We have found an association of young patients (age range 20–29 years) to TCs in Galicia, similarly to the findings in other studies (Brenner et al., 2017; Dennis et al., 2018; Verhofstede et al., 2018; Fabeni et al., 2019; Paraskevis et al., 2019; González-Domenech et al., 2020), but not in the Basque Country. Public health interventions for reducing high risk behavior and control measures in this age group in Galicia should be implemented to decrease HIV-1 transmissions in the region. Moreover, our results highlight that even regions with comparable demographic features can show a different HIV-1 epidemic pattern, likely due to specific social and cultural differences at the regional level.

In this study, we have shown that MSM TCs are a keystone of the HIV-1 epidemic at regional and probably also at the country level, according to the propagation of the large MSM TCs identified in this study in different regions of Spain. Consequently, reinforcement of public health measures for preventing HIV-1 infections in MSM, such as behavioral interventions to reduce risky practices, pre-exposure prophylaxis (Grant et al., 2010; McCormack et al., 2016; Riddell et al., 2018), and early diagnosis and treatment of HIV-1 infections (European Centre for Disease Prevention and Control, 2015; United Nations Population Fund Global Forum on MSM & HIV United Nations Development Programme World Health Organization United States Agency for International Development World Bank, 2015) are recommended.

Implementation of preventive public health measures targeting MSM could reduce indirectly the spread of HIV-1 in Spain, which is mainly MSM-driven, in other population groups. National surveillance data in the United States suggest that infections among heterosexual women predominantly originate from MSM (Oster et al., 2015). Moreover, the presence of female and self-declared heterosexual male individuals in MSM-associated TCs is frequent (Hué et al., 2014; Esbjörnsson et al., 2016; Ragonnet-Cronin et al., 2018; Verhofstede et al., 2018), as we have also observed in the large TCs in our study (Table 3). Mixed demographics in TCs allows the spread of HIV-1 infections between different population groups, as we have seen in the expansion of MSM-associated TCs to heterosexuals in Spain (Delgado et al., 2019), indicating the importance of monitoring the evolution of TCs for implementing proper control measures in each population group. Studies to assess the contribution of MSMs to HIV-1 transmission to other population groups are needed in Spain.

The high representativeness of our cohort, especially in the Basque Country, suggests that the largest and most spread TCs in these regions have been identified, which supports the conclusions of the study. However, our data may have some limitations as some of the NDs from Basque Country and Galicia during the studied period were not analyzed and not all the TCs propagating in Spain were detected in our phylogenetic analysis.

Our study has provided valuable data on the HIV-1 epidemic in Spain, identifying specific features at geographical and population levels for tailoring more specific public health interventions. The implementation of a molecular HIV-1 surveillance system to monitor the evolution of the epidemic in Spain could allow to promptly detect rapidly expanding TCs needing urgent investigation and the implementation of additional public health measures to prevent their spread (Centers for Disease Control and Prevention, 2018), as well as to readjust the control measures in place according to the evolution of the TCs.

Members of the Spanish Group for the Study of New HIV Diagnoses

Basque Country: Hospital Universitario de Araba, Vitoria: Andrés Canut-Blasco, José Joaquín Portu, Carmen Gómez-González; Hospital Universitario de Basurto, Bilbao: Josefa Muñoz, Mª Carmen Nieto, María Zuriñe Zubero, Silvia Hernáez-Crespo, Estibaliz Ugalde; Hospital Universitario de Cruces, Barakaldo: Luis Elorduy, Elena Bereciartua, Leyre López Soria; Hospital de Galdakao: Mª José López de Goicoechea, José Mayo; Hospital Universitario Donostia: Gustavo Cilla, Julio Arrizabalaga, José Antonio Iribarren, Mª Jesús Bustinduy, Mª Julia Echevarría. Mª Yolanda Salicio, David Grandioso. Galicia: Área sanitaria de Ferrol: Ana Mariño, Patricia Ordóñez, Hortensia Álvarez, Nieves Valcarce; Complejo Hospitalario Universitario de A Coruña: Ángeles Cañizares, Mª Ángeles Castro. Hospital Universitario Lucus Augusti, Lugo: Ramón Rabuñal-Rey, María José García-País, Mª José Gude-González, Pilar Alonso-García, Antonio Moreno-Flores; Complejo Hospitalario Universitario de Ourense: Juan García Costa, Ricardo Fernández-Rodríguez, Raúl Rodríguez-Pérez, Jorge Guitián, María Dolores Díaz-López, María Genoveva Naval-Calviño; Complejo Hospitalario Universitario de Vigo: Celia Miralles, Antonio Ocampo, Sonia Pérez-Castro, Jorge Julio Cabrera; Complejo Hospitalario Universitario de Pontevedra: Julio Díz-Arén, Matilde Trigo, Mª Ángeles Pallarés. Navarra: Complejo Hospitalario de Navarra, Pamplona: Carmen Ezpeleta Baquedano, Carmen Martín Salas, Irati Arregui García, María Gracia Ruiz de Alda. Madrid: Centro Sanitario Sandoval, Madrid: Jorge del Romero, Carmen Rodríguez, Mar Vera, Óskar Ayerdi, Eva Orviz; Hospital Universitario de Fuenlabrada: María Isabel García-Arata, Santiago Prieto-Menchero; Hospital Clínico Universitario San Carlos, Madrid: Esther Culebras, Icíar Rodríguez-Avial; Hospital Universitario Fundación Jiménez Díaz, Madrid: Raquel Téllez-Pérez, Olalla Calabia-González, Alfonso Cabello-Úbeda, Miguel Górgolas Hernández-Mora; Hospital Universitario Severo Ochoa, Leganés: Sara María Quevedo, Lucía Puente, Manuel del Álamo; Hospital Fundación Alcorcón, Madrid: Carolina Campelo Gutiérrez, María José Goyanes Galán; Castilla y León: Hospital Clínico Universitario de Valladolid: Carmen Hinojosa, Carlos Dueñas, Begoña Monteagudo, Edita Sánchez; Hospital Río Hortega, Valladolid: Jessica Abadía, Belén Lorenzo Vidal; Hospital Virgen de la Concha, Zamora: Teresa Martín-Domínguez, Rosa Martínez-González. La Rioja: Hospital San Pedro, Logroño: José Ramón Blanco, Miriam Blasco. Aragón: Hospital Universitario Miguel Servet, Zaragoza: Ana María Martínez-Sapiña, Piedad Arazo. Castilla-La Mancha: Hospital Universitario de Guadalajara: Alejandro González Praetorius, Complejo Hospitalario de Toledo: César Gómez-Hernando, José Largo-Pau; Comunitat Valenciana: Hospital Universitari Sant Joan d’Alacant: Fernando Buñuel, Ana Infante.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/genbank/, KT276255, KT276260, KT276263-KT276264, KT276266-KT276267, KT276270-KT276271, KU685562-KU685565, KU685567, KU685569-KU685575, KU685577, KU685581-KU685582, KU685586, KU685588-KU685590, KU685592, KX534325, KX534329, KX534331, KY465968, KY496624, KY514084, KY989950, MF999250-MF999256, MF999258-MF999261, MK177721-MK177733, MK177735-MK177752, MK177754-MK177757, MK177761-MK177772, MK177783-MK17785, MK177787-MK177790, MK177792, MK177795-MK177796, MK177798-MK177799, MK177801-MK177803, MK177805-MK177807, MK298150, MT436242, MT436244, MT436246-MT436247, MT436249-MT436250, MW344920-MW344921, MW584217-MW584224, OK011532, OK011542, OK011545-OK011546, OK011549, OK148912, OK148914, OK148917-OK148919, OK148921, OK148941-OK148942, OK148953-OK148957, OK148961, OK148965, OK148968-OK148972, OK148974, OL982314-OL982315, and OM914651-OM914711.

Ethics Statement

The studies involving human participants were reviewed and approved by Bioethics and Animal Well-being Committee of Instituto de Salud Carlos III. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

HG, ED, and MT conceived the study and supervised the experimental work. ED, HG, MT, and JC-G processed sequences and performed phylogenetic analyses. HG and ED performed data curation. HG, AD, and LG performed statistical analyses. SB, MS, VM, and EG-B performed experimental work. The members of the Spanish Group for the Study of New HIV Diagnoses recruited patients and obtained epidemiological and clinical data. HG wrote the manuscript draft, with contributions to the text by AD, MT, and ED. All authors contributed to the article and approved the submitted version.

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 and Project “Epidemiología molecular del VIH-1 en España y su utilidad para investigaciones biológicas y en vacunas“PI19CIII/0042; 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 ICDCI, project RD16ISCIII/0002/0004; and 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).

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.

Publisher’s Note

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.

Acknowledgments

The authors 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, Centro Nacional de Microbiología, Instituto de Salud Carlos III, for technical assistance in sequencing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.782609/full#supplementary-material

Footnotes

1. ^https://www.hiv.lanl.gov/content/sequence/HIV/mainpage.html

References

AIDS Study Group (GESIDA) of the Spanish Society of Infectious Diseases and Clinical Microbiology and the National AIDS Plan. (2020). Consensus document on antirretroviral therapy in adults infected by the human inmunodeficiency virus (Updated 2020). Available at: https://gesida-seimc.org/wp-content/uploads/2020/07/TAR_GUIA_GESIDA_2020_COMPLETA_Julio.pdf (Accessed July 16, 2020).

Google Scholar

Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Álvarez del Arco, D., Fakoya, I., Thomadakis, C., Pantazis, N., Touloumi, G., Gennotte, A. F., et al. (2017). High levels of postmigration HIV acquisition within nine European countries. AIDS 31, 1979–1988. doi: 10.1097/QAD.0000000000001571

PubMed Abstract | CrossRef Full Text | Google Scholar

Board, A. R., Oster, A. M., Song, R., Gant, Z., Linley, L., Watson, M., et al. (2020). Geographic distribution of HIV transmission networks in the United States. J. Acquir. Immune Defic. Syndr. 85, e32–e40. doi: 10.1097/QAI.0000000000002448

PubMed Abstract | CrossRef Full Text | Google Scholar

Brenner, B. G., Ibanescu, R. I., Hardy, I., Stephens, D., Otis, J., Moodie, E., et al. (2017). Large cluster outbreaks sustain the HIV epidemic among MSM in Quebec. AIDS 31, 707–717. doi: 10.1097/QAD.0000000000001383

PubMed Abstract | CrossRef Full Text | Google Scholar

Brenner, B., Wainberg, M. A., and Roger, M. (2013). Phylogenetic inferences on HIV-1 transmission: implications for the design of prevention and treatment interventions. AIDS 27, 1045–1057. doi: 10.1097/QAD.0b013e32835cffd9

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, A. J. L., Lycett, S. J., Weinert, L., Hughes, G. J., Fearnhill, E., and Dunn, D. T. (2011). Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J. Infect. Dis. 204, 1463–1469. doi: 10.1093/infdis/jir550

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, E. M., Patala, A., Shankar, A., Li, J. F., Johnson, J. A., Westheimer, E., et al. (2020). Phylodynamic analysis complements partner services by identifying acute and unreported HIV transmission. Viruses 12:145. doi: 10.3390/v12020145

PubMed Abstract | CrossRef Full Text | Google Scholar

Centers for Disease Control and Prevention (2018). Detecting, investigating, and responding to HIV transmission clusters. Available at: https://www.cdc.gov/hiv/pdf/funding/announcements/ps18-1802/CDC-HIV-PS18-1802-AttachmentE-Detecting-Investigating-and-Responding-to-HIV-Transmission-Clusters.pdf (Accessed September 13, 2018).

Google Scholar

Cuevas, M. T., Fernández-García, A., Sánchez-García, A., González-Galeano, M., Pinilla, M., Sánchez-Martínez, M., et al. (2009a). Incidence of non-B subtypes of HIV-1 in Galicia, Spain: high frequency and diversity of HIV-1 among men who have sex with men. Euro Surveill. 14:19413. doi: 10.2807/ese.14.47.19413-en

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuevas, M. T., Muñoz-Nieto, M., Thomson, M. M., Delgado, E., Iribarren, J. A., Cilla, G., et al. (2009b). HIV-1 transmission cluster with T215D revertant mutation among newly diagnosed patients from the Basque Country, Spain. J. Acquir. Immune Defic. Syndr. 51, 99–103. doi: 10.1097/QAI.0b013e318199063e

PubMed Abstract | CrossRef Full Text | Google Scholar

Dauwe, K., Mortier, V., Schauvliege, M., Van Den Heuvel, A., Fransen, K., Servais, J. 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. doi: 10.1186/s12879-015-1217-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Delgado, E., Benito, S., Montero, V., Cuevas, M. T., Fernández-García, A., Sánchez-Martínez, M., et al. (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

PubMed Abstract | CrossRef Full Text | Google Scholar

Delgado, E., Cuevas, M. T., Domínguez, F., Vega, Y., Cabello, M., Fernández-García, A., 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 One 10:e0143325. doi: 10.1371/journal.pone.0143325

PubMed Abstract | CrossRef Full Text | Google Scholar

Dennis, A. M., Volz, E., Frost, A. S. M. S., Hossain, M., Poon, A. F. Y., Rebeiro, P. F., et al. (2018). HIV-1 transmission clustering and phylodynamics highlight the important role of young men who have sex with men. AIDS Res. Hum. Retrovir. 34, 879–888. doi: 10.1089/AID.2018.0039

PubMed Abstract | CrossRef Full Text | Google Scholar

Desgrées-du-Loâ, A., Pannetier, J., Ravalihasy, A., Gosselin, A., Supervie, V., Panjo, H., et al. (2015). Sub-Saharan African migrants living with HIV acquired after migration, France, ANRS PARCOURS study, 2012 to 2013. Euro Surveill. 20:30065. doi: 10.2807/1560-7917.ES.2015.20.46.30065

PubMed Abstract | CrossRef Full Text | Google Scholar

Dias, S., Gama, A., Loos, J., Roxo, L., Simões, D., and Nöstlinger, C. (2020). The role of mobility in sexual risk behaviour and HIV acquisition among sub-Saharan African migrants residing in two European cities. PLoS One 15:e0228584. doi: 10.1371/journal.pone.0228584

PubMed Abstract | CrossRef Full Text | Google Scholar

Esbjörnsson, J., Mild, M., Audelin, A., Fonager, J., Skar, H., Bruun-Jørgensen, 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. doi: 10.1093/ve/vew010

PubMed Abstract | CrossRef Full Text | Google Scholar

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 June 17, 2015).

Google Scholar

European Centre for Disease Prevention and Control and WHO Regional Office for Europe (2020). HIV/AIDS surveillance in Europe 2020–2019 data. Available at: https://www.ecdc.europa.eu/sites/default/files/documents/hiv-surveillance-report-2020.pdf (Accessed November 26, 2020).

Google Scholar

Fabeni, L., Alteri, C., Berno, G., Scutari, R., Orchi, N., De, C. G., et al. (2019). Characterisation of HIV-1 molecular transmission clusters among newly diagnosed individuals infected with non-B subtypes in Italy. Sex. Transm. Infect. 95, 619–625. doi: 10.1136/sextrans-2019-054017

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabeni, L., Rozera, G., Berno, G., Giombini, E., Gori, C., Orchi, N., et al. (2021). Molecular transmission dynamics of primary HIV infections in Lazio region, years 2013-2020. Viruses 13:176. doi: 10.3390/v13020176

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabeni, L., Santoro, M. M., Lorenzini, P., Rusconi, S., Gianotti, N., Costantini, A., et al. (2020). Evaluation of HIV transmission clusters among natives and foreigners living in Italy. Viruses 12:791. doi: 10.3390/v12080791

PubMed Abstract | CrossRef Full Text | Google Scholar

Fakoya, I., Ávarez del Arco, D., Woode-Owusu, M., Monge, S., Rivero-Montesdeoca, Y., Delpech, V., et al. (2015). A systematic review of post-migration acquisition of HIV among migrants from countries with generalised HIV epidemics living in Europe: implications for effectively managing HIV prevention programmes and policy. BMC Public Health 15:561. doi: 10.1186/s12889-015-1852-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Fenton, K. A., Chinouya, M., Davidson, O., and Copas, A. (2001). HIV transmission risk among sub-Saharan Africans in London travelling to their countries of origin. AIDS 15, 1442–1445. doi: 10.1097/00002030-200107270-00017

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Alba, J. M., Holguín, A., García, R., García-Bujalance, S., Alonso, R., Suáez, A., et al. (2011). Molecular surveillance of HIV-1 in Madrid, Spain: a phylogeographic analysis. J. Virol. 85, 10755–10763. doi: 10.1128/JVI.00454-11

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Domenech, C. M., Sena-Corrales, G., Viciana-Ramos, I., Palacios-Muñoz, R., Mora-Navas, L., Clavijo-Frutos, E., et al. (2020). High prevalence of sequences included in transmission clusters within newly diagnosed HIV-1 patients in southern Spain (2004-2015). Microb. Drug Resist. 26, 1090–1097. doi: 10.1089/mdr.2019.0344

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Domenech, C. M., Viciana, I., Delaye, L., Mayorga, M. L., Palacios, R., de la Torre, J., et al. (2018). Emergence as an outbreak of the HIV-1 CRF19_cpx variant in treatment-naïve patients in southern Spain. PLoS One 13:e0190544. doi: 10.1371/journal.pone.0190544

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabowski, M. K., Lessler, J., Bazaale, J., Nabukalu, D., Nankinga, J., Nantume, B., et al. (2020). Migration, hotspots, and dispersal of HIV infection in Rakai, Ugand. Nat. Commun. 11:976. doi: 10.1038/s41467-020-14636-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Grant, R. M., Lama, J. R., Anderson, P. L., McMahan, V., Liu, A. Y., Vargas, L., et al. (2010). Preexposure chemoprophylaxis for HIV prevention in men who have sex with men. N. Engl. J. Med. 363, 2587–2599. doi: 10.1056/NEJMoa1011205

PubMed Abstract | CrossRef Full Text | Google Scholar

Grossman, Z., Avidor, B., Mor, Z., Chowers, M., Levy, I., Shahar, E., et al. (2015). A population-structured HIV epidemic in Israel: roles of risk and ethnicity. PLoS One 10:e0135061. doi: 10.1371/journal.pone.0135061

PubMed Abstract | CrossRef Full Text | Google Scholar

Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, T. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp. Ser. 41, 95–98.

Google Scholar

Hassan, A. S., Pybus, O. G., Sanders, E. J., Albert, J., and Esbjornsson, J. (2017). Defining HIV-1 transmission clusters based on sequence data. AIDS 31, 1211–1222. doi: 10.1097/QAD.0000000000001470

PubMed Abstract | CrossRef Full Text | Google Scholar

HIV STI and Hepatitis Surveillance Unit (2020). Vigilancia del VIH y sida en España: Sistema de información sobre nuevos diagnósticos de VIH y registro nacional de casos de sida. Plan Nacional sobre el Sida. Available at: https://www.mscbs.gob.es/ca//ciudadanos/enfLesiones/enfTransmisibles/sida/vigilancia/Informe_VIH_SIDA_20201130.pdf (Accessed November 27, 2020).

Google Scholar

Hoenigl, M., Chaillon, A., Kessler, H. H., Haas, B., Stelzl, E., Weninger, K., et al. (2016). Characterization of HIV transmission in south-East Austria. PLoS One 11:e0151478. doi: 10.1371/journal.pone.0151478

PubMed Abstract | CrossRef Full Text | Google Scholar

Holguín, A., Pena, M. J., Troncoso, F., and Soriano, V. (2007). Introduction of non-B subtypes among Spaniards newly diagnosed with HIV type 1 in the Canary Islands. AIDS Res. Hum. Retrovir. 23, 498–502. doi: 10.1089/aid.2006.0191

PubMed Abstract | CrossRef Full Text | Google Scholar

Hué, S., Brown, A. E., Ragonnet-Cronin, M., Lycett, S. J., Dunn, D. T., Fearnhill, E., et al. (2014). Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions. AIDS 28, 1967–1975. doi: 10.1097/QAD.0000000000000383

PubMed Abstract | CrossRef Full Text | Google Scholar

Hué, S., Clewley, J. P., Cane, P. A., and Pillay, D. (2004). HIV-1 pol gene variation is sufficient for reconstruction of transmissions in the era of antiretroviral therapy. AIDS 18, 719–728. doi: 10.1097/00002030-200403260-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Hué, S., Pillay, D., Clewley, J. P., and Pybus, O. G. (2005). Genetic analysis reveals the complex structure of HIV-1 transmission within defined risk groups. Proc. Natl. Acad. Sci. U. S. A. 102, 4425–4429. doi: 10.1073/pnas.0407534102

PubMed Abstract | CrossRef Full Text | Google Scholar

Hughes, G. J., Fearnhill, E., Dunn, D., Lycett, S. J., Rambaut, A., and Leigh Brown, A. J. (2009). Molecular phylodynamics of the heterosexual HIV epidemic in the United Kingdom. PLoS Pathog. 5:e1000590. doi: 10.1371/journal.ppat.1000590

PubMed Abstract | CrossRef Full Text | Google Scholar

Instituto Nacional de Estadistica (2021). Población residente. Available at: https://www.ine.es/jaxiT3/Datos.htm?t=9681 (Accessed September 8, 2021).

Google Scholar

Lole, K. S., Bollinger, R. C., Paranjape, R. S., Gadkari, D., Kulkarni, S. S., Novak, N. 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. 73, 152–160. doi: 10.1128/JVI.73.1.152-160.1999

PubMed Abstract | CrossRef Full Text | Google Scholar

McCormack, S., Dunn, D. T., Desai, M., Dolling, D. I., Gafos, M., Gilson, R., 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. Lancet 387, 53–60. doi: 10.1016/S0140-6736(15)00056-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Nuñez, O., Hernando, V., and Díaz, A. (2018). Estimating the number of people living with HIV and the undiagnosed fraction in Spain in 2013. AIDS 32, 2573–2581. doi: 10.1097/QAD.0000000000001989

PubMed Abstract | CrossRef Full Text | Google Scholar

Osorno-González de León, M. F., Delgado, E., Benito, S., Sánchez, M., Gil, H., Montero, V., et al. (2021). 11th Conference on HIV Science, abstract PEC252. Available at: https://ias2021.org/wp-content/uploads/2021/07/IAS2021_Abstracts_web.pdf (Accessed July 18–21, 2021).

Google Scholar

Oster, A. M., France, A. M., Panneer, N., Bañez-Ocfemia, M. C., Campbell, E., Dasgupta, S., et al. (2018). Identifying clusters of recent and rapid HIV transmission through analysis of molecular surveillance data. J. Acquir. Immune Defic. Syndr. 79, 543–550. doi: 10.1097/QAI.0000000000001856

PubMed Abstract | CrossRef Full Text | Google Scholar

Oster, A. M., Wertheim, J. O., Hernandez, A. L., Ocfemia, M. C., Saduvala, N., and Hall, H. I. (2015). Using molecular HIV surveillance data to understand transmission between subpopulations in the United States. J. Acquir. Immune Defic. Syndr. 70, 444–451. doi: 10.1097/QAI.0000000000000809

PubMed Abstract | CrossRef Full Text | Google Scholar

Paraskevis, D., Beloukas, A., Stasinos, K., Pantazis, N., de Mendoza, C., Bannert, N., et al. (2019). HIV-1 molecular transmission clusters in nine European countries and Canada: association with demographic and clinical factors. BMC Med. 17:4. doi: 10.1186/s12916-018-1241-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Paraskevis, D., Nikolopoulos, G. K., Magiorkinis, G., Hodges-Mameletzis, I., and Hatzakis, A. (2016). The application of HIV molecular epidemiology to public health. Infect. Genet. Evol. 46, 159–168. doi: 10.1016/j.meegid.2016.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Parczewski, M., Leszczyszyn-Pynka, M., Witak-Jedra, M., Szetela, B., Gasiorowski, J., Knysz, B., et al. (2017). Expanding HIV-1 subtype B transmission networks among men who have sex with men in Poland. PLoS One 12:e0172473. doi: 10.1371/journal.pone.0172473

PubMed Abstract | CrossRef Full Text | Google Scholar

Patiño-Galindo, J. A., Thomson, M. M., Pérez-Álvarez, L., Delgado, E., Cuevas, M. T., Fernández-García, A., et al. (2016). Transmission dynamics of HIV-1 subtype B in the Basque Country, Spain. Infect. Genet. Evol. 40, 91–97. doi: 10.1016/j.meegid.2016.02.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Patiño-Galindo, J. A., Torres-Puente, M., Bracho, M. A., Alastrué, I., Juan, A., Navarro, D., et al. (2017a). Identification of a large, fast-expanding HIV-1 subtype B transmission cluster among MSM in Valencia, Spain. PLoS One 12:e0171062. doi: 10.1371/journal.pone.0171062

PubMed Abstract | CrossRef Full Text | Google Scholar

Patiño-Galindo, J. A., Torres-Puente, M., Bracho, M. A., Alastrué, I., Juan, A., Navarro, D., et al. (2017b). The molecular epidemiology of HIV-1 in the Comunidad Valenciana (Spain): analysis of transmission clusters. Sci. Rep. 7:11584. doi: 10.1038/s41598-017-10286-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Pérez-Parra, S., Chueca, N., Ávarez, M., Pasquau, J., Omar, M., Collado, A., et al. (2016). Phylodynamic and phylogeographic profiles of subtype B HIV-1 epidemics in South Spain. PLoS One 11:e0168099. doi: 10.1371/journal.pone.0168099

PubMed Abstract | CrossRef Full Text | Google Scholar

Pérez-Parra, S., Chueca-Porcuna, N., Ávarez-Estevez, M., Pasquau, J., Omar, M., Collado, A., et al. (2015). Los estudios de resistencia a antirretrovirales como herramienta para el análsis de los clusters de transmisión del virus de la inmunodefiencia humana. Enferm. Infecc. Microbiol. Clin. 33, 603–608. doi: 10.1016/j.eimc.2014.11.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Petersen, A., Cowan, S. A., Nielsen, J., Fischer, T. K., and Fonager, J. (2018). Characterisation of HIV-1 transmission clusters and drug-resistant mutations in Denmark, 2004 to 2016. Euro Surveill. 23:1700633. doi: 10.2807/1560-7917.ES.2018.23.44.1700633

PubMed Abstract | CrossRef Full Text | Google Scholar

Pineda-Peña, A. C., Pingarilho, M., Li, G., Vrancken, B., Libin, P., Gomes, P., et al. (2019). Drivers of HIV-1 transmission: the Portuguese case. PLoS One 14:e0218226. doi: 10.1371/journal.pone.0218226

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2--approximately maximum-likelihood trees for large alignments. PLoS One 5:e9490. doi: 10.1371/journal.pone.0009490

PubMed Abstract | CrossRef Full Text | Google Scholar

Ragonnet-Cronin, M., Hué, S., Hodcroft, E. B., Tostevin, A., Dunn, D., Fawcett, T., et al. (2018). Non-disclosed men who have sex with men in UK HIV transmission networks: phylogenetic analysis of surveillance data. Lancet HIV 5, e309–e316. doi: 10.1016/S2352-3018(18)30062-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ragonnet-Cronin, M., Lycett, S. J., Hodcroft, E. B., Hué, S., Fearnhill, E., Brown, A. E., et al. (2016a). Transmission of non-B HIV subtypes in the United Kingdom is increasingly driven by large non-heterosexual transmission clusters. J. Infect. Dis. 213, 1410–1418. doi: 10.1093/infdis/jiv758

PubMed Abstract | CrossRef Full Text | Google Scholar

Ragonnet-Cronin, M., Shilaih, M., Günthard, H. F., Hodcroft, E. B., Böni, J., FearnHill, E., et al. (2016b). A direct comparision of two densenly sampled HIV epidemics: the UK and Switzerland. Sci. Rep. 6:32251. doi: 10.1038/SREP32251

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramirez, J. J. D., Ballouz, T., Nguyen, H., Kusejko, K., Chaudron, S. E., Huber, M., et al. (2021). Increasing frequency and transmission of HIV-1 non-B subtypes among men who have sex with men in the Swiss HIV cohort study. J. Infect. Dis. 225, 306–316. doi: 10.1093/infdis/jiab360

PubMed Abstract | CrossRef Full Text | Google Scholar

Riddell, J., Amico, K. R., and Mayer, K. H. (2018). HIV preexposure prophylaxis: a review. JAMA 319, 1261–1268. doi: 10.1001/jama.2018.1917

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimodaira, H., and Hasegawa, M. (1999). Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol. Biol. Evol. 16, 1114–1116. doi: 10.1093/oxfordjournals.molbev.a026201

CrossRef Full Text | Google Scholar

Su, L., Liang, S., Hou, X., Zhong, P., Wei, D., Fu, Y., et al. (2018). Impact of worker emigration on HIV epidemics in labour export areas: a molecular epidemiology investigation in Guangyuan, China. Sci. Rep. 8:16046. doi: 10.1038/s41598-018-33996-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, M. M., Delgado, E., Manjón, N., Ocampo, A., Villahermosa, M. L., Mariño, A., et al. (2001). HIV-1 genetic diversity in Galicia Spain: BG intersubtype recombinant viruses circulating among injecting drug users. AIDS 15, 509–516. doi: 10.1097/00002030-200103090-00010

PubMed Abstract | CrossRef Full Text | Google Scholar

United Nations Population Fund Global Forum on MSM & HIV United Nations Development Programme World Health Organization United States Agency for International Development World Bank (2015). Implementing comprehensive HIV and STI programmes with men who have sex with men: practical guidance for collaborative interventions. Available at: https://www.unfpa.org/es/node/13155 (Accessed September 8, 2015).

Google Scholar

Vasylyeva, T. I., Friedman, S. R., Paraskevis, D., and Magiorkinis, G. (2016). Integrating molecular epidemiology and social network analysis to study infectious diseases: towards a socio-molecular era for public health. Infect. Genet. Evol. 46, 248–255. doi: 10.1016/j.meegid.2016.05.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Vasylyeva, T. I., Zarebski, A., Smyrnov, P., Williams, L. D., Korobchuk, A., Liulchuk, M., et al. (2020). Phylodynamics helps to evaluate the impact of an HIV prevention intervention. Viruses 12:469. doi: 10.3390/v12040469

PubMed Abstract | CrossRef Full Text | Google Scholar

Vega, Y., Delgado, E., Fernández-García, A., Cuevas, M. T., Thomson, M. M., Montero, V., 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 One 10:e0125699. doi: 10.1371/journal.pone.0125699

PubMed Abstract | CrossRef Full Text | Google Scholar

Verhofstede, C., Dauwe, K., Fransen, K., Van, L. K., Van den Wijngaert, S., Ruelle, J., 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. 61, 36–44. doi: 10.1016/j.meegid.2018.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Viciana, I., González-Domenech, C. M., Palacios, R., Delgado, M., Del, A. A., Tellez, F., et al. (2016). Clinical, virological and phylogenetic characterization of a multiresistant HIV-1 strain outbreak in naive patients in southern Spain. J. Antimicrob. Chemother. 71, 357–361. doi: 10.1093/jac/dkv332

PubMed Abstract | CrossRef Full Text | Google Scholar

Vinken, L., Fransen, K., Cuypers, L., Alexiev, I., Balotta, C., Debaisieux, L., et al. (2019). Earlier initiation of antiretroviral treatment coincides with an initial control of the HIV-1 sub-subtype F1 outbreak among men-having-sex-with-men in Flanders, Belgium. Front. Microbiol. 10:613. doi: 10.3389/fmicb.2019.00613

PubMed Abstract | CrossRef Full Text | Google Scholar

Wertheim, J. O., Kosakovsky Pond, S. L., Forgione, L. A., Mehta, S. R., Murrell, B., Shah, S., et al. (2017). Social and genetic networks of HIV-1 transmission in New York City. PLoS Pathog. 13:e1006000. doi: 10.1371/journal.ppat.1006000

PubMed Abstract | CrossRef Full Text | Google Scholar

Yebra, G., Holguín, A., Pillay, D., and Hué, S. (2013). Phylogenetic and demographic characterization of HIV-1 transmission in Madrid, Spain. Infect. Genet. Evol. 14, 232–239. doi: 10.1016/j.meegid.2012.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar