ABSTRACT
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 1158 newly HIV-1-diagnosed patients (NDs) from both regions diagnosed 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 and determine the 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% confidence intervals (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. Differences in distribution and sizes of local TCs were observed in both regions reflecting a different epidemic pattern. Our results show that TCs play an important role in the spread of HIV-1 infection 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.
INTRODUCCION
The HIV-1 epidemic is still one of the major public health problems in Spain. Around 3,500 - 4,000 new diagnoses of HIV-1 infection (NDs) 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 (Board et al., 2020; Brenner et al., 2013; Paraskevis et al., 2016; Wertheim et al., 2017). 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 reconstructing transmission events through the identified TCs and infer the history of the HIV-1 epidemic (Hué et al., 2004; Hué et al., 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 (Fabeni et al., 2021; Oster et al., 2018; Parczewski et al., 2017; Petersen et al., 2018; Pineda-Peña et al., 2019; Verhofstede et al., 2018). 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 (Board et al., 2020; Brenner et al., 2013; Campbell et al., 2020; Paraskevis et al., 2019; Vasylyeva et al., 2016; Vasylyeva 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 (Cuevas et al., 2009a; Cuevas et al., 2009b; González-Alba et al., 2011; González-Domenech et al., 2020; Holguín et al., 2007; Patiño-Galindo et al., 2016; Patiño-Galindo et al., 2017b; Pérez-Parra et al., 2016; Pérez-Parra et al., 2015; Thomson et al., 2001; Yebra et al., 2013). The phylogenetic studies have also allowed the identification of large TCs among MSM which are actively growing in Spain (Delgado et al., 2015; Delgado et al., 2019; González-Domenech et al., 2018; Patiño-Galindo et al., 2017a), as well as transmitted drug resistance mutations which are spreading in TCs (Cuevas et al., 2009b; González-Domenech et al., 2018; González-Domenech et al., 2020; Vega et al., 2015; Viciana et al., 2016), 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 diagnose 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 analyze the role of TCs in the transmission of HIV-1 infection 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.
MATERIAL 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.
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 and 2.
Population sequencing was performed with ABI Prism BigDye Terminator Cycle Sequencing kit and ABI 3730 XL sequencer (Applied Biosystems, Foster City, CA, U.S.A.) 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, USA) and edited with BioEdit v.7.2.5 (Hall, 1999) (http://www.mbio.ncsu.edu).
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 1999-2019 were included in the analyses. Sequences were analyzed phylogenetically by an approximately maximum-likelihood method using FastTree2 (Price et al., 2010). 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. 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)
TCs 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.
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 p-value <0.05. A multivariable logistic regression model was performed to identify factors associated to TCs. The model was adjusted by sex, 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% confidence interval (CI). Data analyses were performed using the STATA statistical software package Version 16 (Stata Corporation, College Station, TX, US).
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).
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, being higher in Sub-Saharan Africans than in Latin Americans: female [52% (44/85) vs. 26% (37/180)] and heterosexual transmission [88% (65/74) vs. 33%. (54/165)], respectively.
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%). 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 with 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 (Table 2), 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 84%, than in Galicia, where it was 41% 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 (Table 2).
New HIV-1 diagnoses from Basque Country and Galicia (2013-18) included in the study and reported to the HIV National Surveillance System.
TCs identified in the study
Fifty one percent (594/1158) of the individuals belonged to 162 different TCs. Regarding the 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 3). Eighteen (11%) TCs comprised patients from both regions. The sizes of the TCs ranged from 4 to 204 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 3), although the difference was not statistically significant (p=0.207).
Characteristic of the identified TCs of the study.
A total of 15 large TCs (here defined as those with ≥30 individuals) were identified in the study (Table 4), 10 of them of subtype B. MSM was the associated transmission route in 14 (93%) of them. 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 4). 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 4).
Characteristics of the transmission clusters comprising ≥ 30 patients in the Basque Country and Galicia (2013-2018).
Factors associated with TCs
The percentage of patients in TCs was higher in males (57%) than in females (22%) and in MSM (67%) compared to heterosexuals (31%). This percentage was also higher in Spaniards (64%) than in Latin Americans (35%) or in Sub-Saharan Africans (1.2%) and in individuals infected with viruses of subtypes F (76%) or B (58%) compared to 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 5). 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
Characteristics of the new HIV diagnoses included in TCs (2013-2018).
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 6). 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 6)
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 6). 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 6).
DISCUSSION
The phylogenetic analysis of 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 84%.
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 4 documented infections. Therefore, our more strict TC definition, compared to some other studies, has allowed us to identify the structure and factors associated to HIV-1 strains whose transmission is established in Spain.
TCs 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 (Duran-Ramirez et al., 2021; Parczewski et al., 2017; Wertheim et al., 2017). Also, studies in other Spanish regions have found this association of TCs with MSM (González-Domenech et al., 2020; Pérez-Parra et al., 2015). 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 2000 (Nuñez et al., 2018).
Spaniards showed higher odds than Latin Americans or Africans of belonging to TCs. 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. Similarly to our results, a low percentage of Sub-Saharan Africans in TCs have been also reported in Italy (Fabeni et al., 2020) and Belgium (Verhofstede et al., 2018). 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.
In Spain, 36% of HIV-1 infections diagnosed in 2019 were in migrants (HIV, STI and hepatitis surveillance unit, 2020). Determining the origin of HIV-1 infections in the migrant population would allow proper allocation of resources for preventive measures (Álvarez del Arco et al., 2017; Fakoya et al., 2015). A previous study based on models using clinical and behavioral data found a high level of post-migration HIV acquisition, up to 71% of the migrants in Spain (Álvarez del Arco et al., 2017). However, phylogenetic analyses can provide additional information on the geographical area of HIV-1 acquisition, decreasing the bias due to data collection or failure to consider the infections acquired by migrants when visiting their countries (Álvarez del Arco et al., 2017), as migrants’ mobility was associated with increased risk for HIV acquisition (Dias et al., 2020). 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.
Subtype B was the most frequent genetic form associated to TCs. Only subtype F infections (76% of NDs in TCs) showed a higher probability 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; Delgado et al., 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 (Delgado et al., 2019). This phenomenon has been also documented in other European countries (Dauwe et al., 2015; Duran Ramirez et al., 2021; Fabeni et al., 2019; Ragonnet-Cronin et al., 2016; Verhofstede et al., 2018), where viruses introduced from abroad have expanded successfully among the MSM population.
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; Fabeni et al., 2019; González-Domenech et al., 2020; Paraskevis et al., 2019; Verhofstede et al., 2018), 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 regional level.
In this study we have shown that MSM TCs are a keystone of the HIV-1 epidemic at regional and country levels. 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 Developmenty, World Bank, 2015) are recommended. In addition, new strategies and other tools for reaching the MSM population, such as the apps for sex dating suggested in other studies (Duncan et al., 2018; Wirden et al., 2019) could be considered for preventing HIV-1 infections.
Implementation of measures targeting MSM could reduce indirectly the spread of HIV-1, which is mainly MSM-driven, in other population groups. National surveillance data in the USA 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 TCs is frequent (Esbjörnsson et al., 2016; Hué et al., 2014; Ragonnet-Cronin et al., 2018; Verhofstede et al., 2018), as we have also observed in large TCs in Spain (Delgado et al., 2019). Analyses to evaluate the contribution of MSM in 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. New approaches and strategies are needed to reach the UNAIDS 90-90-90 target proposed by 2020 (UNAIDS, 2014). To achieve this goal, the implementation of a molecular HIV-1 surveillance system to monitor the evolution of the epidemic in Spain could be of great help. Such a surveillance system 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.
Data Availability
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation. In
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, Ma 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: Ma José López de Goicoechea, José Mayo; Hospital Universitario Donostia: Gustavo Cilla, Julio Arrizabalaga, José Antonio Iribarren, Ma Jesús Bustinduy, Ma Julia Echevarría. Ma 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, Ma Ángeles Castro. Hospital Universitario Lucus Augusti, Lugo: Ramón Rabuñal-Rey, María José García-País, Ma 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, Ma Á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.
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).
ETHICS STATEMENT
Sequences were derived from antiretroviral resistance genotypic 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.
AUTHOR CONTRIBUTIONS
HG, ED, and MT conceived the study and supervised the experimental work. ED, HG, MT, and JC 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 read and approved the text.
ACKNOWLEDGEMENTS
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, Centro Nacional de Microbiología, Instituto de Salud Carlos III, for technical assistance in sequencing.
REFERENCES
