Next Article in Journal
Complex Virome in a Mesenteric Lymph Node from a Californian Sea Lion (Zalophus californianus) with Polyserositis and Steatitis
Next Article in Special Issue
Employing Molecular Phylodynamic Methods to Identify and Forecast HIV Transmission Clusters in Public Health Settings: A Qualitative Study
Previous Article in Journal
Endogenous Retroviruses Walk a Fine Line between Priming and Silencing
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 8
10.3390/v12080791
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of HIV Transmission Clusters among Natives and Foreigners Living in Italy

by
Lavinia Fabeni
1,†,
Maria Mercedes Santoro
2,*,†,
Patrizia Lorenzini
3,
Stefano Rusconi
4,
Nicola Gianotti
5,
Andrea Costantini
6,
Loredana Sarmati
7,
Andrea Antinori
3,
Francesca Ceccherini-Silberstein
2,
Antonella d’Arminio Monforte
8,
Annalisa Saracino
9,
Enrico Girardi
10 and
on behalf of the Icona Foundation Study Cohort
‡
1
Laboratory of Virology, National Institute for Infectious Diseases “Lazzaro Spallanzani” IRCCS, 00149 Rome, Italy
2
Department of Experimental Medicine, University of Rome “Tor Vergata”, 00133 Rome, Italy
3
Clinical Department, National Institute for Infectious Diseases “Lazzaro Spallanzani” IRCCS, 00149 Rome, Italy
4
III Division of Infectious Diseases, ASST Fatebenefratelli Sacco, Luigi Sacco University Hospital, 20157 Milan, Italy
5
Infectious Diseases, San Raffaele Scientific Institute, 20127 Milan, Italy
6
Department of Clinical and Molecular Science, Marche Polytechnic University, 60121 Ancona, Italy
7
Clinical Infectious Diseases, Policlinic Foundation Tor Vergata, 00133 Rome, Italy
8
Department of Health Sciences, University of Milan, 00161 Milan, Italy
9
Department of Biochemical Sciences and Human Oncology, University of Bari, 70121 Bari, Italy
10
Clinical Epidemiology Unit, National Institute for Infectious Diseases “Lazzaro Spallanzani” IRCCS, 00149 Rome, Italy
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
‡
Membership of the Icona Foundation Study Cohort is provided in the Acknowledgments.
Viruses 2020, 12(8), 791; https://doi.org/10.3390/v12080791
Submission received: 10 June 2020 / Revised: 6 July 2020 / Accepted: 8 July 2020 / Published: 23 July 2020
(This article belongs to the Special Issue HIV Molecular Epidemiology for Prevention 2020)
Browse Figure
Review Reports Versions Notes

Abstract

:
We aimed at evaluating the characteristics of HIV-1 molecular transmission clusters (MTCs) among natives and migrants living in Italy, diagnosed between 1998 and 2018. Phylogenetic analyses were performed on HIV-1 polymerase (pol) sequences to characterise subtypes and identify MTCs, divided into small (SMTCs, 2–3 sequences), medium (MMTCs, 4–9 sequences) and large (LMTCs, ≥10 sequences). Among 3499 drug-naïve individuals enrolled in the Italian Cohort Naive Antiretroviral (ICONA) cohort (2804 natives; 695 migrants), 726 (20.8%; 644 natives, 82 migrants) were involved in 228 MTCs (6 LMTCs, 36 MMTCs, 186 SMTCs). Migrants contributed 14.4% to SMTCs, 7.6% to MMTCs and 7.1% to LMTCs, respectively. HIV-1 non-B subtypes were found in 51 MTCs; noteworthy was that non-B infections involved in MTCs were more commonly found in natives (n = 47) than in migrants (n = 4). Factors such as Italian origin, being men who have sex with men (MSM), younger age, more recent diagnosis and a higher CD4 count were significantly associated with MTCs. Our findings show that HIV-1 clustering transmission among newly diagnosed individuals living in Italy is prevalently driven by natives, mainly MSM, with a more recent diagnosis and frequently infected with HIV-1 non-B subtypes. These results can contribute to monitoring of the HIV epidemic and guiding the public health response to prevent new HIV infections.

1. Introduction

Although HIV-1 is preventable through effective public health measures, substantial HIV transmission continues to occur worldwide. In 2018, 1.7 million HIV-1 newly infected individuals were reported [1]. In particular, 141,552 newly diagnosed HIV infections were reported in 50 of the 53 countries in the European Region of the World Health Organization (WHO) [2], while 2847 new cases were reported in Italy [3].
Among the HIV-1-diagnosed individuals in the EU/EEA in 2018, 42% were migrants, defined as originating from outside of the country in which they were diagnosed [2]. In Italy, migrants accounted for nearly 30% of all newly diagnosed HIV infections in recent years [2,3]. As a consequence of the high percentage of migrants in several European countries including Italy, the prevalence of numerous HIV-1 non-B subtypes and circulating recombinant forms (CRFs) has rapidly increased in several previously B-restricted areas [4,5,6,7,8,9]. However, the increasing prevalence of non-B subtypes is only partially due to migrants, since a substantial amount of such increase is attributable to the European-born population, including Italians [4,5,7,8,9,10].
A considerable proportion of the new HIV diagnoses is involved in well-defined clusters prevalently characterized by high-risk sexual behaviours [6,8,11,12,13,14]. Accurate identification of these transmission clusters and an understanding of HIV-transmission dynamics are very important to efficiently target public health interventions. In this context, phylogenetic analysis represents one of the most important strategies to better describe and monitor the local HIV epidemics. When combined with epidemiological and clinical data, the results of such analysis can be of public health relevance, for example by identifying how virus lineages are restricted to, or mix among, different demographic and behavioural subgroups [15,16].
The aim of this study was to evaluate molecular transmission clusters (MTCs) among native and foreign individuals living in Italy diagnosed between 1998 and 2018. Characteristics and drivers of the HIV-1 epidemic in Italy were also evaluated by combining MTCs analysis with demographic, epidemiological and laboratory data.

2. Materials and Methods

2.1. Cohort Description

The Italian Cohort Naive Antiretroviral (ICONA) Foundation Study is a multi-centre prospective observational cohort of HIV-1-infected patients, which was set up in 1997 [17]. Eligible patients are those starting combined antiretroviral therapy when they are naive to antiretrovirals, regardless of the reason for which they had never been previously treated and of the stage of their disease. All patients sign a consent form to participate in ICONA, in accordance with the ethical standards of the committee on human experimentation and the Helsinki Declaration (1983 revision). Demographic, clinical and laboratory data and information on therapy are collected for all participants and recorded using electronic data collection; sensitive data are collected only in an anonymous format. Subjects enrolled in the ICONA cohort are broadly representative of the Italian HIV-diagnosed population [18].

2.2. Study Population

For this specific study, 3499 newly diagnosed patients enrolled in the ICONA cohort during the period between 1998 and 2018, for whom an HIV-1 pol sequence was available at the time of the diagnosis, were considered for the analysis. Individuals born in a country other than Italy were classified as migrants. Pol sequences (containing the full-length protease (PR) and the first 240/335 reverse transcriptase (RT) codons) were produced from plasma samples through the classical Sanger sequencing performed according to the protocols used at the laboratories at each Clinical Center contributing to the cohort.

2.3. Phylogenetic Analyses

Molecular phylogeny was used to determine both HIV-1 subtype for each sequence, as previously described [19], and for the cluster analysis. MTCs were first deduced using the HIV-TRACE (Transmission Cluster Engine) web tool [20,21]. Pairwise genetic distances were obtained using HIV-TRACE under the Tamura-Nei 93 (TN93) nucleotide substitution model. This is the most general nucleotide substitution model for which distances can be estimated directly from counts of nucleotide pairs in aligned sequences [21]. Putative transmission links were inferred when two sequences had TN93 genetic distance ≤0.010 substitutions/site. This stringent genetic distance threshold has been chosen according to the suggestions by Wertheim and colleagues [22], based on the consideration that within a single person, HIV pol sequences tend not to diverge more than 0.01 substitutions/site from the baseline sequence in the first 10 years of infection [23], with the total sequence divergence tending to be less than 0.02 substitutions/site [24]. Therefore, one would expect a meaningful epidemiologically genetic distance threshold for identifying transmission partners to fall between 0.01 and 0.02 substitutions/site.
To avoid the influence of convergence evolution at antiretroviral drug resistance mutations, sequences were stripped at positions related to drug resistance.
The robustness of the MTCs was further tested using the maximum likelihood (ML) method and a Bayesian analysis. MTCs observed by HIV-TRACE analysis were not considered if not retrieved by ML and/or Bayesian phylogenetic analysis. The ML tree was inferred as previously described, using the MEGA version 6.06 software [25]. MTCs were recognised by bootstrap values >90% and a <0.015 average genetic distance. Finally, the Bayesian phylogenetic tree was reconstructed with MrBayes version 3.2.2, using a GTR + I + Г5. The Monte Carlo Markov Chain search was run for 5 × 106 generations, with the trees sampled every 100th generation (with a burn-in of 50%) [26]. Statistical support was obtained by calculating the posterior probability of each monophyletic clade, and a posterior consensus tree was generated after 50% burn-in. Clades were considered epidemiological clusters only if a posterior probability of ≥0.90 was inferred.
MTCs were then divided into small clusters (SMTCs, 2–3 sequences), medium clusters (MMTCs, 4–9 sequences) and large clusters (LMTCs, ≥10 sequences).

2.4. Evaluation of HIV-1 Transmitted Drug Resistance

Transmitted drug resistance (TDR) was evaluated by considering the WHO 2009 list [27] with the additional RT mutations K65E/N, E138A/G/K/Q/R, V179L, G190Q, T215N, H221Y, F227C and M230I, reported in the International AIDS Society list [28] and/or the HIVdb V.8.9-1 (https://hivdb.stanford.edu/, last updated 25/10/2019). HIV-1 strains were defined as resistant if carrying at least one TDR mutation.

2.5. Statistical Analysis

Comparisons between patients belonging or not to clusters were evaluated using the Mann–Whitney test for quantitative variables and Chi-square or Fisher’s exact test for categorical variables. Factors potentially associated with MTCs were evaluated using multivariable logistic regression analysis. In particular, the following factors were included in the model: gender, mode of HIV transmission, nation of birth, level of education, type of employment, plasma HIV RNA and CD4 cell count at sequencing and calendar year of HIV diagnosis. All analyses were performed using STATA version 15.1 (StataCorp College Station, Texas, USA). p-values <0.05 were considered as statistically significant.

2.6. Sequences Deposited to GenBank

All 726 sequences involved in the 228 molecular transmission clusters have been deposited into GenBank (accession numbers: MT787750-MT788351; MT788352-MT788382; MT788383-MT788475). All other sequences are available on request.

3. Results

3.1. Patients’ Characteristics

A total of 3499 HIV-1 newly diagnosed participants in the ICONA cohort from 1997 through 2018 were included (Table 1). The majority of these individuals (2872; 82.1%) were male, Italian (2804; 80.1%) and MSM (1789; 51.1%). Migrants accounted for 19.9%, had lived in Italy for a median (interquartile range, IQR) of 5 (1–9) years and were represented mainly by subjects from Central and South America (241, 6.9%) and Africa (219, 6.3%), followed by Europe (187, 5.3%), Asia (38, 1.1%) and other areas (10, 0.3%) (Table 1 and Table S1). Phylogenetic analysis revealed that most individuals (n = 2556; 73.1%) were infected with B subtype, followed by CRF02_AG (187; 5.3%), F1 (179; 5.1%), C (148; 4.2%), A1 (104; 3.0%), CRF60_BC (64; 51.8%) and other subtypes/CRFs (261; 7.5%).
The overall prevalence of TDR was 14.3% (n = 501), mainly imputable to resistance associated with non-nucleoside RT inhibitors (NNRTI) (9.6%).

3.2. Cluster Identification

Overall, 228 MTCs including 2–52 individuals were identified by both HIV-TRACE and phylogenetic analyses (Figure 1, Figures S1 and S2). No differences between the two methods were identified. MTCs involved 726 individuals (20.7%; 644 natives and 82 migrants). Of the 228 MTCs identified, 186 were SMTCs (2–3 sequences) involving 402 individuals (11.5%), 36 MMTCs (4–9 sequences) involving 184 individuals (5.3%) and 6 LMTCs (≥10 sequences) involving 140 individuals (4.0%) (Table 2). The presence of both natives and migrants was found in 66.7% of LMTCs, 33.3% of MMTCs and 23.1% of SMTCs. Migrants were involved mainly in SMTCs (n = 58; 70.7%). Individuals infected with HIV-1 B subtype in MTCs were present with a higher prevalence than those infected with non-B subtype (71.4% vs. 28.6%). Individuals infected with non-B subtypes were part of 51 MTCs (22.4%) (SMTCs: 37; MMTCs: 11; LMTCs: 3; Table 2); only four of these MTCs were entirely composed of migrants (7.8%).
By analyzing the date of diagnosis of individuals involved in MTCs, it was found that before 2010 the proportion of individuals involved in the clusters was 34.6% (n = 251: SMTCs, n = 156, 38.8%; MMTCs/LMTCs = 95, 29.3%). Only a few individuals were involved in MTCs before 2007 (n = 38, 5.2%: SMTCS, n = 30, 7.5%; MMTCs/LMTCs = 8, 2.5%). Starting from 2011, the proportion of individuals involved in MTCs was 65.4% (n = 475: SMTCs, n = 246, 61.2%; MMTCs/LMTCs = 229, 70.7%). Migrants, on the other hand, were only involved in MTCs from 2007. In particular, the proportion of migrants in MTCs between 2007 and 2010 was overall 26.8% (n = 22: SMTCs, n = 18, 31.0%; MMTCs/LMTCs = 4, 16.7%), whereas after 2011 the proportion of migrants in MTCs was 73.2% (n = 60: SMTCs, n = 40, 69.0%; MMTCs/LMTCs = 20, 83.3%).
Regarding TDR, 81 (11.2%) individuals involved in MTCs carried a resistant virus. This proportion was lower compared to that found in individuals out of MTCs (15.1%, p = 0.006). Individuals with a resistant virus were involved in 36 SMTCs (29 B, 2 C, 2 CRF01_AE, 1 CRF02_AG, 1 CRF24_BG, 1 F1), five MMTCs (three B, one CRF02_AG, and one CRF12_BF) and one LMTC (CRF60_BC).

3.3. Epidemiological Characteristics and Factors Associated with MTCs

With regard to the characteristics of individuals involved in MTCs in comparison with those out of MTCs, male gender (94.3% vs. 78.9%, p < 0.001), being MSM (75.8% vs. 45.0%, p < 0.001) and natives (88.7% vs. 77.8%, p < 0.001) were significantly associated with being in MTCs (Table 1). Individuals in MTCs were also younger (median (IQR) age: 32 (27–40) vs. 38 (31–46) years, p < 0.001), more recently diagnosed (median (IQR) year of diagnosis: 2012 (2009–2014) vs. 2011 (2007–2014), p < 0.001), and had higher CD4 cell count on diagnosis (median (IQR) cells/mm3: 459 (322–624) vs. 353 (177–523), p < 0.001) in comparison with individuals out of MTCs. Moreover, a lower prevalence of TDR was found among individuals in MTCs compared to those out of MTCs (11.2% vs. 15.1%, p < 0.001). Multivariable logistic regression confirmed that Italian origin, being MSM, of younger age, a more recent diagnosis, and higher CD4 cell count were factors significantly associated with MTCs (Table 1).

3.4. Characteristics of the Medium/Large Clusters

Information of MMTCs (n = 36) and LMTCs (n = 6) is reported in Table 2 and Table 3. The majority of MTCs were composed only of males (MMTCs: 32/36, 88.9%; LMTCs: 5/6, 83.3%). The most common risk factor was same-gender sex (49.8% of the MSM in the overall MTCs), followed by heterosexual contacts (24.8% of the heterosexuals in the overall MTCs). As regards to migrants, they contributed 7.6% to MMTCs and 7.1% to LMTCs, respectively. The 24 migrants involved in these medium/large MTCs were mainly from Central/South America or other European countries (Table 3).
Some of the individuals involved in these MTCs carried the following pathways of mutations associated with drug resistance: I85V in PR (n = 1), and A62V (n = 3), E138A (n = 6) or K101E (n = 5) in RT (Table 3). The median genetic distance was ≤0.016 substitution/site in each MTC, highlighting the tight correlation among individuals involved in the MTC.

3.5. Characteristics of the Small Clusters

Focusing attention on the 186 SMTCs, 156 (83.9%) were composed of two persons, and 30 (16.1%) of three persons (Table 2). The contribution of females in these SMTCs was minimal (34/402 (8.5%): 31 in pairs and 3 in MTCs composed of three individuals). Overall, these SMTCs were composed mainly of MSM (54.3%). Interestingly, 32.8% included individuals with mixed transmission categories, mainly among MSM and heterosexual contacts (19.9%) predominantly found in TCs (86.5%) pairs. The great majority of these TCs were exclusively composed of Italians (73.7%), while about 23.1% of them were composed of mixed nationality. Only six (3.2%) of the SMTCs were totally composed of migrant couples (Table 2) originating from Africa (n = 3), Asia (n = 1), America (n = 1) and Eastern Europe (n = 1) (data not shown); four of them (67%) were composed of heterosexual individuals that came from the same country in Africa (n = 3) and Asia (n = 1) (Figure 1), while the other two were composed of MSM.
Regarding TDR, 66 individuals (involved in 36 SMTCs, 19.3%) carried drug resistance mutations.
In particular, seven SMTCs (12 individuals) carried only PI resistance mutations (K23I, M46I/L, F53Y, I85V, L90M); twenty-four SMTCs (43 individuals) carried NNRTI resistance mutations (L100I, K101E, K103N/S, V106I, E138A/G/K, G190A), alone or in combination with other RTI resistance mutations. Specifically, one SMTC, composed of two individuals, carried the combination K103N+E138A, and one SMTC had the combination L100I + K103N + L210W + T215D.
Five SMTCs (11 individuals) carried NRTI resistance mutations (M41L, D67N, L210W, T215S/L/D, K219Q), alone or in combination with other NRTI mutations. In particular, the following combinations were found in all the individuals involved in the MTCs: M41L+T215L (three individuals), D67N+K219Q (two individuals), while one pair was composed of one individual with the M41L and one with the M41L+T215D combination.

4. Discussion

This study provides molecular evidence of HIV-1 epidemiological clusters circulating in Italy among 3499 native and foreign individuals diagnosed between 1998 and 2018 enrolled in the ICONA cohort. Around 21% of the overall individuals included in the analysis were involved in clusters. Clustering analysis showed the presence of a high number of MTCs; in particular, 228 MTCs including 6 LMTCs, 36 MMTCs and 186 SMTCs were identified.
In general, clustered transmission was prevalently driven by native individuals (644, 88.7%); in fact, 222 MTCs (97.4%) were composed mainly of Italian individuals, and more in particular, 157 MTCs (70.7%) were composed exclusively of Italian persons and 65 MTCs (29.3%) showed an intermixing between migrants and Italians. Only six SMTCs of the overall MTCs were composed exclusively of a couple of migrant individuals—67% were composed by heterosexual individuals and came from the same country in Africa and Asia. This highlights how some migrants probably acquired HIV in their home country and how in certain geographical areas the transmission route is mainly represented by heterosexual contacts [29]. A minimal contribution of migrants in MTCs was also found in Belgium, despite the high prevalence of foreign individuals living in this country [10]. In particular, in the Belgian cohort, the percentage of foreign people was much higher than in our cohort (53.0% vs. 19.9%). Furthermore, foreign individuals infected with HIV (for example those from Europe) might look slightly less actively involved in cluster formation in Italy than in Belgium, although both cohorts showed similar tendencies. Despite these differences between the two cohorts, overall findings support the fact that clustering is mainly driven by natives.
Individuals involved in MTCs were younger, mainly male and MSM, with a recent diagnosis and higher CD4 count compared to those out of MTCs. These findings are in agreement with several studies on transmission clusters from the USA, Asia and Europe [8,10,14,19,30,31,32,33,34,35,36] and led us to hypothesise that MSM are more aware than others of having both high-risk behaviour and a high self-perception of risk for acquiring HIV, and for this reason, they test more frequently than others and are thus diagnosed earlier.
Even though MSM was the most prevalent risk factor in MTCs, different risk behaviours were observed among individuals belonging to the same MTCs. In particular, mixtures of MSM and heterosexual contacts were observed in most MTCs. In fact, at least one heterosexual contact was observed in 19.9%, 27.8%, and 83.0% of each SMTCs, MMTCs, and LMTCs, respectively. Noteworthy was that the majority of MTCs found the number of MSM clearly outweighed the number of heterosexuals, and in many, only men were represented. This observation is in line with that reported in previous studies and suggests that misreporting of transmission risk may be an important reason for the presence of mixed clusters [10,37,38]. The most important reason for misreporting could be fear or embarrassment to disclose sex with same-gender sex. Another explanation for the presence of mixed transmission risk clusters may be due to the fact that homosexual orientation is frequently masked by bisexual behaviours in older individuals, and men who have sex with men and women are a potential bridge population for transmitting HIV to heterosexual women [6,10,39,40]. Along the same line, intravenous drug users may serve as a go-between risk category.
Of note, about a quarter of MTCs were characterized by HIV-1 non-B subtypes; in particular, 50% of LMTCs and 30% of MMTCs were driven by these viral strains. By considering the fact that only a small portion of migrants was present in these non-B MTCs, in the most part of cases interspersed with natives, this finding confirms the fact that HIV-1 non-B variants have entered and are to date circulating in newly diagnosed individuals originating from Italy. This is in line with the important increase of non-B pure subtypes and CRFs reported in several European countries during the last decade. This increase in non-B subtypes occurred in conjunction with associated epidemiological changes, such as (1) increase in immigration and tourism from areas with non-B subtypes, (2) travel-associated infections and (3) increased internal circulation in Italy of imported subtypes [5,7,8,9,10,12,14,41,42,43,44,45,46,47].
By specifically considering the date of diagnosis, most individuals in MTCs were diagnosed after 2010, thus suggesting that the expansion of epidemics in Italy seems mainly associated with recent MTCs, rather than the growth of older, established ones. This has also been confirmed by multivariable logistic regression analysis, showing that a more recent diagnosis is a factor positively associated with MTCs. The involvement of migrants in MTCs has a similar temporal trend of the overall population. From the public health point of view, all this information is important for designing prevention and public health intervention strategies.
Regarding TDR, 14.3% of the overall individuals harboured a resistant virus, mainly imputable to NNRTI resistance mutations. This TDR prevalence is slightly higher than in other studies [7,39,46,48] because in the analysis, we considered the E138A mutation, associated with resistance to etravirine and rilpivirine [49,50,51,52]. In fact, by removing this mutation, the overall prevalence of TDR in our study population was found to be 9.8%, similar to that already reported in literature [7,10,42,46,48,53], remaining mainly attributed to NNRTI resistance (4.8%). As E138A is a polymorphic mutation that ranges in prevalence from about 1% to 5% depending on subtype [54], its high prevalence in our population might not necessarily be related to TDR.
By evaluating the prevalence of TDR in and out of MTCs, similarly to that observed in previous studies [8,46], a resistant virus was less frequent among individuals involved in MTCs compared with those out (11.2% vs. 15.1%), and this was also confirmed by excluding the E138A mutation (data not shown). On the other hand, different results have been found by Verhofstede et al. [10] describing a slightly higher TDR prevalence within the clusters. These discrepancies might be explained by different characteristics of the study populations analyzed. Resistance to more than one class was found only in one SMTCs, where mutations for both NRTIs and NNRTIs were present. By assessing TDR prevalence in MTCs according to drug class, we mainly found the mutations associated with resistance to NNRTIs K101E and E138A, which are at low fitness cost [55,56].
As with any other observational study, our data may have some limitations. Firstly, our study population might not be representative of the overall population living in Italy. In particular, it has been estimated that to date, migrant individuals with new HIV-1 diagnoses accounted for nearly 30% of all the newly diagnosed HIV infections in Italy in recent years [3]. Despite the fact that in the overall ICONA cohort this percentage is confirmed [17], in our specific study population, the overall percentage of migrants was lower (around 20%). This is due to the fact that the sequences at diagnosis were not available for all the individuals enrolled in the cohort. Thus the contribution of migrants in MTCs could be underestimated; in fact, unsampled individuals might limit the correct idea of the migrant circulation in our country. However, our results are in line with other previous Italian studies in different study populations that showed that the contribution of migrants in the cluster is minimal [6,8,46].
Moreover, we should take into account that phylogenetic analyses are limited by their inability to infer direction of transmission: the absence of the contribution of sequences from drug-treated individuals and unsampled individuals might limit correct transmission dynamics. Finally, another significant limitation to this study is the lack of information on the country of infection, because these are important data to establish the geographical origin of the MTCs. This could explain, for example, the presence of couples in our study population composed only of heterosexual individuals of the same gender, even though we cannot exclude a misclassification of the risk factor for these cases. Nonetheless, the detection of MTCs and the associated demographic and viroimmunological parameters can help us to focus proven HIV prevention tools where they are needed most.
In conclusion, this study shows that HIV-1 newly diagnosed subjects are involved in several MTCs in the past two decades in Italy. Clustered transmission, especially for large clusters, is prevalently driven by natives, mainly MSM, with a recent diagnosis and frequently infected with HIV-1 non-B subtype. Participation of migrants in clustered transmission is rare.
These results reinforce the fact that phylodynamics represent one of the most important tools to better describe and monitor local HIV epidemics, by correlating the genetic relationship of the viruses with information on demographics, transmission mode and new infections.
Overall, our findings can contribute to monitoring of the HIV epidemic and guiding the public health response in Italy to prevent new infections.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/8/791/s1, Table S1: Characteristics of foreign individuals with HIV living in Italy, Figure S1: Phylogenetic tree of the 3499 pol sequences analysed and coloured according to the most representative subtypes, Figure S2. Phylogenetic tree of the 726 pol sequences involved in 228 molecular transmission clusters (MTCs).

Author Contributions

Conceptualization, M.M.S., F.C.-S., A.d.M., A.S. and E.G.; methodology, L.F., M.M.S. and P.L.; software, L.F. and P.L.; validation, L.F., M.M.S. and P.L.; formal analysis, L.F. and P.L.; investigation, L.F. and M.M.S.; data curation, P.L.; writing—original draft preparation, L.F. and M.M.S.; writing—review and editing, L.F. and M.M.S.; visualization, S.R., N.G., A.C., L.S., A.A.; supervision, L.F., M.M.S., F.C.-S., A.d.M., A.S. and E.G. All authors have read and agreed to the published version of the manuscript.

Funding

Icona Foundation has received unrestricted grants and grants for specific projects by Gilead, ViiV Healthcare, MSD. No grants for this specific study.

Acknowledgments

A part of this work was presented at the Conference on Retroviruses and Opportunistic Infections (CROI). Boston, Massachusetts, 8-11 March, 2020 (Abstract ID 920). We gratefully thank Debra Mandatori for having revised and edited the manuscript. Icona Fundation Study Group: BOARD OF DIRECTORS: A d’Arminio Monforte (President), A Antinori (Vice-President), M Andreoni, A Castagna, F Castelli, R Cauda, G Di Perri, M Galli, R Iardino, G Ippolito, A Lazzarin, GC Marchetti, G Rezza, F von Schloesser, P Viale. SCIENTIFIC SECRETARY: A d’Arminio Monforte, A Antinori, A Castagna, F Ceccherini-Silberstein, A Cozzi-Lepri, E Girardi, A Gori, S Lo Caputo, F Maggiolo, C Mussini, M Puoti, CF Perno. STEERING COMMITTEE: A Antinori, F Bai, A Bandera, S Bonora, M Borderi, A Calcagno, MR Capobianchi, A Castagna, F Ceccherini-Silberstein, S Cicalini, A Cingolani, P Cinque, A Cozzi-Lepri, A d’Arminio Monforte, A Di Biagio, R Gagliardini, E Girardi, N Gianotti, A Gori, G Guaraldi, G Lapadula, M Lichtner, A Lai, S Lo Caputo, G Madeddu, F Maggiolo, G Marchetti, E Merlini, C Mussini, S Nozza, CF Perno, S Piconi, C Pinnetti, M Puoti, E Quiros Roldan, R Rossotti, S Rusconi, MM Santoro, A Saracino, L Sarmati, V Spagnuolo, V Svicher, L Taramasso. STATISTICAL AND MONITORING TEAM: A Cozzi-Lepri, I Fanti, L Galli, P Lorenzini, A Rodano’, M Macchia, A Tavelli. COMMUNITY ADVISORY BOARD: A Bove, A Camposeragna, M Errico, M Manfredini, A Perziano, V Calvino. BIOLOGICAL BANK INMI: F Carletti, S Carrara, A Di Caro, S Graziano, F Petroni, G Prota, S Truffa. PARTICIPATING PHYSICIANS AND CENTERS: Italy A Giacometti, A Costantini, V Barocci (Ancona); G Angarano, L Monno, E Milano (Bari); F Maggiolo, C Suardi (Bergamo); P Viale, V Donati, G Verucchi (Bologna); F Castelnuovo, C Minardi, E Quiros Roldan (Brescia); B Menzaghi, C Abeli (Busto Arsizio); L Chessa, F Pes (Cagliarti); B Cacopardo, B Celesia (Catania); J Vecchiet, K Falasca (Chieti); A Pan, S Lorenzotti (Cremona); L Sighinolfi, D Segala (Ferrara); P Blanc, F Vichi (Firenze); G Cassola, M Bassetti, A Alessandrini, N Bobbio, G Mazzarello (Genova); M Lichtner, L Fondaco, (Latina); P Bonfanti, C Molteni (Lecco); A Chiodera, P Milini (Macerata); G Nunnari, G Pellicanò (Messina); A d’Arminio Monforte, M Galli, A Lazzarin, G Rizzardini, M Puoti, A Castagna, ES Cannizzo, MC Moioli, R Piolini, D Bernacchia, A Poli, C Tincati, (Milano); C Mussini, C Puzzolante (Modena); C Migliorino, G Lapadula (Monza); V Sangiovanni, G Borgia, V Esposito, G Di Flumeri, I Gentile, V Rizzo (Napoli); AM Cattelan, S Marinello (Padova); A Cascio, M Trizzino (Palermo); D Francisci, E Schiaroli (Perugia); G Parruti, F Sozio (Pescara); C Lazzaretti, R Corsini (Reggio Emilia); M Andreoni, A Antinori, R Cauda, A Cristaudo, V Vullo, R Acinapura, S Lamonica, M Capozzi, A Mondi, A Cingolani, M Rivano Capparuccia, G Iaiani, A Latini, G Onnelli, MM Plazzi, G De Girolamo, A Vergori (Roma); M Cecchetto, F Viviani (Rovigo); G Madeddu, A De Vito(Sassari); B Rossetti, F Montagnani (Siena); A Franco, R Fontana Del Vecchio (Siracusa); C Di Giuli (Terni); P Caramello, G Di Perri, S Bonora, GC Orofino, M Sciandra (Torino); A Londero (Udine); V Manfrin, G Battagin (Vicenza); G Starnini, A Ialungo (Viterbo).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. The Joint United Nations Program on HIV/AIDS (UNAIDS). Available online: http://www.unaids.org/en (accessed on 28 April 2020).
  2. ECDC: Surveillance Report. HIV/AIDS Surveillance in Europe 2019 – 2018 Data. Available online: https://www.ecdc.europa.eu/sites/default/files/documents/hiv-surveillance-report-2019.pdf (accessed on 28 April 2020).
  3. Regine, V.; Pugliese, L.; Boros, S.; Santaquilani, M.; Ferri, M.; Suligoi, B. Aggiornamento delle nuove diagnosi di infezione da HIV e dei casi di AIDS in Italia al 31 dicembre 2018. Notiziario ISS October 2019, 32:10. Available online: http://www.salute.gov.it/imgs/C_17_notizie_3963_0_file.pdf (accessed on 28 April 2020).
  4. Monno, L.; Brindicci, G.; Lo Caputo, S.; Punzi, G.; Scarabaggio, T.; Riva, C.; Di Bari, C.; Pierotti, P.; Saracino, A.; Lagioia, A.; et al. HIV-1 subtypes and circulating recombinant forms (CRFs) from HIV-infected patients residing in two regions of central and southern Italy. J. Med. Virol. 2005, 75, 483–490. [Google Scholar] [CrossRef] [PubMed]
  5. Dauwe, K.; Mortier, V.; Schauvliege, M.; Van Den Heuvel, A.; Fransen, K.; Servais, J.Y.; Bercoff, D.P.; Seguin-Devaux, C.; Verhofstede, C. Characteristics and spread to the native population of HIV-1 non-B subtypes in two European countries with high migration rate. BMC Infect. Dis. 2015, 15, 524. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Fabeni, L.; Alteri, C.; Orchi, N.; Gori, C.; Bertoli, A.; Forbici, F.; Montella, F.; Pennica, A.; De Carli, G.; Giuliani, M.; et al. Recent transmission clustering of HIV-1 C and CRF17_BF strains characterized by NNRTI -Related mutations among newly diagnosed men in central Italy. PLoS ONE 2015, 10, e0135325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Rossetti, B.; Di Giambenedetto, S.; Torti, C.; Postorino, M.C.; Punzi, G.; Saladini, F.; Gennari, W.; Borghi, V.; Monno, L.; Pignataro, A.R.; et al. Evolution of transmitted HIV-1 drug resistance and viral subtypes circulation in Italy from 2006 to 2016. HIV Med. 2018, 19, 619–628. [Google Scholar] [CrossRef] [PubMed]
  8. Fabeni, L.; Alteri, C.; Berno, G.; Scutari, R.; Orchi, N.; De Carli, G.; Bertoli, A.; Carioti, L.; Gori, C.; Forbici, F.; et al. Characterisation of HIV-1 molecular transmission clusters among newly diagnosed individuals infected with non-B subtypes in Italy. Sex. Transm. Infect. 2019, 95, 619–625. [Google Scholar] [CrossRef] [PubMed]
  9. Sagnelli, C.; Uberti-Foppa, C.; Bagaglio, S.; Cella, E.; Scolamacchia, V.; Hasson, H.; Salpietro, S.; Messina, E.; Morsica, G.; Angeletti, S.; et al. Molecular epidemiology of HIV-1 infection in immigrant population in northern Italy. Epidemiol. Infect. 2020, 148, e19. [Google Scholar] [CrossRef] [Green Version]
  10. Verhofstede, C.; Dauwe, K.; Fransen, K.; Van Laethem, K.; Van den Wijngaert, S.; Ruelle, J.; Delforge, M.L.; Vancutsem, E.; Vaira, D.; Stoffels, K.; et al. 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. 2018, 61, 36–44. [Google Scholar] [CrossRef]
  11. Brenner, B.G.; Roger, M.; Routy, J.P.; Moisi, D.; Ntemgwa, M.; Matte, C.; Baril, J.G.; Thomas, R.; Rouleau, D.; Bruneau, J.; et al. High rates of forward transmission events after acute/early HIV-1 infection. J. Infect. Dis. 2007, 195, 951–959. [Google Scholar] [CrossRef]
  12. Lai, A.; Bozzi, G.; Franzetti, M.; Binda, F.; Simonetti, F.R.; Micheli, V.; Meraviglia, P.; Corsi, P.; Bagnarelli, P.; De Luca, A.; et al. Phylogenetic analysis provides evidence of interactions between Italian heterosexual and South American homosexual males as the main source of national HIV-1 subtype C epidemics. J. Med. Virol. 2014, 86, 729–736. [Google Scholar] [CrossRef]
  13. Ragonnet-Cronin, M.; Lycett, S.J.; Hodcroft, E.B.; Hué, S.; Fearnhill, E.; Brown, A.E.; Delpech, V.; Dunn, D.; Leigh Brown, A.J.; United Kingdom HIV Drug Resistance Database. Transmission of Non-B HIV Subtypes in the United Kingdom Is Increasingly Driven by Large Non-Heterosexual Transmission Clusters. J. Infect. Dis. 2016, 213, 1410–1418. [Google Scholar] [CrossRef] [Green Version]
  14. Kostaki, E.G.; Soulie, C.; Visseaux, B.; Storto, A.; Charpentier, C.; Wirden, W.; Landman, R.; Katlama, C.; Calvez, V.; Descamps, D.; et al. Molecular Analysis Suggests Post-Migration HIV-1 Acquisition among Migrants in Paris. Virtual Conference on Retroviruses and Opportunistic Infections. Boston, Massachusetts, 2020. Abstract Number 921. Available online: https://www.croiconference.org/abstract/molecular-analysis-suggests-post-migration-hiv-1-acquisition-among-migrants-in-paris/ (accessed on 28 April 2020).
  15. Brenner, B.; Wainberg, M.A.; Roger, M. Phylogenetic inferences on HIV-1 transmission: Implications for the design of prevention and treatment interventions. AIDS 2013, 27, 1045–1057. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Beltrán, C.G.; Kostaki, E.G.; Carioti, L.; Alvarez, M.; Olalla, J.; Vidal Ampurdanes, M.C.; Montero, M.; García-Bujalance, S.; Blanco, J.R.; Rivero, M.; et al. HIV trace vs phylogenetic Analysis: Unraveling Transmission Clusters in Spain. Virtual Conference on Retroviruses and Opportunistic Infections. Boston, Massachusetts, 2020, Abstract Number 927. Available online: https://www.croiconference.org/abstract/hiv-trace-vs-phylogenetic-analysis-unraveling-transmission-clusters-in-spain/ (accessed on 28 April 2020).
  17. The Italian Cohort Naive Antiretroviral (ICONA) Foundation Study. Available online: http://www.fondazioneicona.org/ (accessed on 12 July 2020).
  18. Vourli, G.; Pharris, A.; Cazein, F.; Costagliola, D.; Dabis, F.; Del Amo, J.; Delpech, V.; Díaz, A.; Girardi, E.; Gourlay, A.; et al. Are European HIV cohort data within EuroCoord representative of the diagnosed HIV population? AIDS 2019, 33, 133–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Fabeni, L.; Berno, G.; Fokam, J.; Bertoli, A.; Alteri, C.; Gori, C.; Forbici, F.; Takou, D.; Vergori, A.; Zaccarelli, M.; et al. Comparative Evaluation of Subtyping Tools for Surveillance of Newly Emerging HIV-1 Strains. J. Clin. Microbiol. 2017, 55, 2827–2837. [Google Scholar] [CrossRef] [Green Version]
  20. HIV-TRACE (Transmission Cluster Engine) Web Tool. Available online: http://hivtrace.datamonkey.org/hivtrace (accessed on 12 July 2020).
  21. Kosakovsky Pond, S.L.; Weaver, S.; Leigh Brown, A.J.; Wertheim, J.O. HIV-TRACE (TRAnsmission Cluster Engine): A Tool for Large Scale Molecular Epidemiology of HIV-1 and Other Rapidly Evolving Pathogens. Mol. Biol. Evol. 2018, 35, 1812–1819. [Google Scholar] [CrossRef] [Green Version]
  22. Wertheim, J.O.; Kosakovsky Pond, S.L.; Forgione, L.A.; Mehta, S.R.; Murrell, B.; Shah, S.; Smith, D.M.; Scheffler, K.; Torian, L.V. Social and Genetic Networks of HIV-1 Transmission in New York City. PLoS Pathog. 2017, 13, e1006000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hightower, G.K.; May, S.J.; Pérez-Santiago, J.; Pacold, M.E.; Wagner, G.A.; Little, S.J.; Richman, D.D.; Mehta, S.R.; Smith, D.M.; Pond, S.L. HIV-1 clade B pol evolution following primary infection. PloS ONE 2013, 8, e68188. [Google Scholar] [CrossRef] [Green Version]
  24. Poon, A.F.; Joy, J.B.; Woods, C.K.; Shurgold, S.; Colley, G.; Brumme, C.J.; Hogg, R.S.; Montaner, J.S.; Harrigan, P.R. The impact of clinical, demographic and risk factors on rates of HIV transmission: A population-based phylogenetic analysis in British Columbia, Canada. J. Infect. Dis. 2015, 211, 926–935. [Google Scholar] [CrossRef] [Green Version]
  25. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  26. Ronquist, F.; Huelsenbeck, J.P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef] [Green Version]
  27. Bennett, D.E.; Camacho, R.J.; Otelea, D.; Kuritzkes, D.R.; Fleury, H.; Kiuchi, M.; Heneine, W.; Kantor, R.; Jordan, M.R.; Schapiro, J.M.; et al. Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update. PLoS ONE 2009, 4, e4724. [Google Scholar] [CrossRef] [Green Version]
  28. Wensing, A.M.; Calvez, V.; Ceccherini-Silberstein, F.; Charpentier, C.; Günthard, H.F.; Paredes, R.; Shafer, R.W.; Richman, D.D. 2019 Update of the Drug Resistance Mutations in HIV-1. Top. Antivir. Med. 2019, 27, 111–121. [Google Scholar] [PubMed]
  29. Fakoya, I.; Álvarez-del Arco, D.; Woode-Owusu, M.; Monge, S.; Rivero-Montesdeoca, Y.; Delpech, V.; Rice, B.; Noori, T.; Pharris, A.; Amato-Gauci, A.J.; et al. 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. 2015, 15, 561. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Aldous, J.L.; Pond, S.K.; Poon, A.; Jain, S.; Qin, H.; Kahn, J.S.; Kitahata, M.; Rodriguez, B.; Dennis, A.M.; Boswell, S.L.; et al. Characterizing HIV transmission networks across the United States. Clin. Infect. Dis. 2012, 55, 1135–1143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Drescher, S.M.; von Wyl, V.; Yang, W.L.; Böni, J.; Yerly, S.; Shah, C.; Aubert, V.; Klimkait, T.; Taffé, P.; Furrer, H.; et al. Swiss HIV Cohort Study. Treatment-naive individuals are the major source of transmitted HIV-1 drug resistance in men who have sex with men in the Swiss HIV Cohort Study. Clin. Infect. Dis. 2014, 58, 285–294. [Google Scholar] [CrossRef] [Green Version]
  32. Brenner, B.G.; Ibanescu, R.I.; Hardy, I.; Stephens, D.; Otis, J.; Moodie, E.; Grossman, Z.; Vandamme, A.M.; Roger, M.; Wainberg, M.A. The Montreal PHI, SPOT cohorts. Large cluster outbreaks sustain the HIV epidemic among MSM in Quebec. AIDS 2017, 31, 707–717. [Google Scholar] [CrossRef]
  33. Lunar, M.M.; Mlakar, J.; Poljak, M. Increase in the numbers of HIV-1 non-B subtypes and potential recombinant forms circulating among Slovenian MSM in recent years. Virus Evol. 2017, 3 (Suppl. 1). [Google Scholar] [CrossRef] [Green Version]
  34. Patiño-Galindo, J.Á.; Torres-Puente, M.; Bracho, M.A.; Alastrué, I.; Juan, A.; Navarro, D.; Galindo, M.J.; Gimeno, C.; Ortega, E.; González-Candelas, F. Identification of a large, fast-expanding HIV-1 subtype B transmission cluster among MSM in Valencia, Spain. PLoS ONE 2017, 12, e0171062. [Google Scholar] [CrossRef] [Green Version]
  35. Parczewski, M.; Leszczyszyn-Pynka, M.; Witak-Jędra, M.; Szetela, B.; Gąsiorowski, J.; Knysz, B.; Bociąga-Jasik, M.; Skwara, P.; Grzeszczuk, A.; Jankowska, M.; et al. Expanding HIV-1 subtype B transmission networks among men who have sex with men in Poland. PLoS ONE 2017, 12, e0172473. [Google Scholar] [CrossRef] [Green Version]
  36. Stecher, M.; Chaillon, A.; Eberle, J.; Behrens, G.M.N.; Eis-Hübinger, A.M.; Lehmann, C.; Jablonka, A.; Bogner, J.; Fätkenheuer, G.; Spinner, C.D.; et al. Molecular Epidemiology of the HIV Epidemic in Three German Metropolitan Regions - Cologne/Bonn, Munich and Hannover, 1999-2016. Sci. Rep. 2018, 8, 6799. [Google Scholar] [CrossRef]
  37. Hué, S.; Brown, A.E.; Ragonnet-Cronin, M.; Lycett, S.J.; Dunn, D.T.; Fearnhill, E.; Dolling, D.I.; Pozniak, A.; Pillay, D.; Delpech, V.C.; et al. UK Collaboration on HIV Drug Resistance and the Collaborative HIV, Anti-HIV Drug Resistance Network (CHAIN). Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions. AIDS 2014, 28, 1967–1975. [Google Scholar] [CrossRef] [PubMed]
  38. Van de Laar, T.J.; Bezemer, D.; van Laethem, K.; Vandewalle, G.; de Smet, A.; van Wijngaerden, E.; Claas, E.C.; van Sighem, A.I.; Vandamme, A.M.; Compernolle, V.; et al. Phylogenetic evidence for underreporting of male-to-male sex among human immunodeficiency virus-infected donors in the Netherlands and Flanders. Transfusion 2017, 57, 1235–1247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Tao, J.; Ruan, Y.; Yin, L.; Vermund, S.H.; Shepherd, B.E.; Shao, Y.; Qian, H.Z. Sex with women among men who have sex with men in China: Prevalence and sexual practices. AIDS Patient Care STDS 2013, 27, 524–528. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Bruhn, C.A.W.; Audelin, A.M.; Helleberg, M.; Bjorn-Mortensen, K.; Obel, N.; Gerstoft, J.; Nielsen, C.; Melbye, M.; Medstrand, P.; Gilbert, M.T.P.; et al. The origin and emergence of an HIV-1 epidemic: From introduction to endemicity. AIDS 2014, 28, 1031–1040. [Google Scholar] [CrossRef] [PubMed]
  41. Lai, A.; Riva, C.; Marconi, A.; Balestrieri, M.; Razzolini, F.; Meini, G.; Vicenti, I.; Rosi, A.; Saladini, F.; Caramma, I.; et al. Changing patterns in HIV-1 non-B clade prevalence and diversity in Italy over three decades. HIV Med. 2010, 11, 593–602. [Google Scholar] [CrossRef] [PubMed]
  42. García, F.; Pérez-Cachafeiro, S.; Guillot, V.; Alvarez, M.; Pérez-Romero, P.; Pérez-Elías, M.J.; Viciana, I.; Blanco, J.R.; López-Dieguez, M.; de Mendoza, C.; et al. Transmission of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treatment naïve HIV-infected individuals (CoRIS). Antiviral. Res. 2011, 91, 150–153. [Google Scholar] [CrossRef]
  43. Patiño Galindo, J.A.; Torres-Puente, M.; Gimeno, C.; Ortega, E.; Navarro, D.; Galindo, M.J.; Navarro, L.; Navarro, V.; Juan, A.; Belda, J.; et al. CRIVIH. Expansion of the CRF19_cpx variant in Spain. J. Clin. Virol. 2015, 69, 146–149. [Google Scholar] [CrossRef]
  44. Esbjörnsson, J.; Mild, M.; Audelin, A.; Fonager, J.; Skar, H.; Bruun Jørgensen, L.; Liitsola, K.; Björkman, P.; Bratt, G.; Gisslén, M.; et al. SPREAD/ESAR Programme. HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries. Virus Evol. 2016, 2, vew010. [Google Scholar] [CrossRef] [Green Version]
  45. Chaillon, A.; Essat, A.; Frange, P.; Smith, D.M.; Delaugerre, C.; Barin, F.; Ghosn, J.; Pialoux, G.; Robineau, O.; Rouzioux, C.; et al. Spatiotemporal dynamics of HIV-1 transmission in France (1999-2014) and impact of targeted prevention strategies. Retrovirol. 2017, 14, 15. [Google Scholar] [CrossRef] [Green Version]
  46. Fabeni, L.; Alteri, C.; Di Carlo, D.; Orchi, N.; Carioti, L.; Bertoli, A.; Gori, C.; Forbici, F.; Continenza, F.; Maffongelli, G.; et al. Dynamics and phylogenetic relationships of HIV-1 transmitted drug resistance according to subtype in Italy over the years 2000–14. J. Antimicrob. Chemother. 2017, 72, 2837–2845. [Google Scholar] [CrossRef]
  47. Stecher, M.; Chaillon, A.; Stephan, C.; Wasmuth, J.C.; Eberle, J.; Behrens, G.; Roider, J.; Spinner, C.D.; Müller, M.C.; Knops, E.; et al. The Translational Platform HIV (TPHIV). Non-B Subtypes HIV 1 Infections in Germany before and after the European Migrant Crisis. Virtual Conference on Retroviruses and Opportunistic Infections. Boston, Massachusetts, 2020, Abstract Number 919. Available online: https://www.croiconference.org/abstract/non-b-subtype-hiv-infections-in-germany-before-and-after-the-european-migrant-crisis/ (accessed on 28 April 2020).
  48. Pineda-Peña, A.C.; Schrooten, Y.; Vinken, L.; Ferreira, F.; Li, G.; Trovão, N.S.; Khouri, R.; Derdelinckx, I.; De Munter, P.; Kücherer, C.; et al. Trends and predictors of transmitted drug resistance (TDR) and clusters with TDR in a local Belgian HIV-1 epidemic. PLoS ONE 2014, 9, e101738. [Google Scholar] [CrossRef]
  49. Haddad, M.; Stawiski, E.; Benhamida, J.; Coakley, E. Improved Genotypic Algorithm for Predicting Etravirine Susceptibility: Comprehensive List of Mutations Identified through Correlation with Matched Phenotype. 17th Conference on Retroviruses and Opportunistic Infections. San Francisco, California, 2010, Abstract Number 574. Available online: http://hivandhepatitis.com/2010_conference/croi/posters/Haddad.pdf (accessed on 28 April 2020).
  50. Haddad, M.; Napolitano, L.A.; Paquet, A.C.; Evans, M.D.; Petropoulos, C.J.; Whitcomb, J.; Rimsky, L.; Vingerhoets, J.; Picchio, G.; Coakley, E. Impact of HIV-1 reverse transcriptase E138 mutations on rilpivirine drug susceptibility. International Workshop on HIV and Hepatitis Virus Drug Resistance and Curative Strategies. Los Cabos, Mexico. Antivir. Ther. 2011, 16 (Suppl. 1), A18. [Google Scholar]
  51. Tambuyzer, L.; Nijs, S.; Daems, B.; Picchio, G.; Vingerhoets, J. Effect of mutations at position E138 in HIV-1 reverse transcriptase on phenotypic susceptibility and virologic response to etravirine. J. Acquir. Immune. Defic. Syndr. 2011, 58, 18–22. [Google Scholar] [CrossRef] [PubMed]
  52. Giannini, A.; Vicenti, I.; Materazzi, A.; Boccuto, A.; Dragoni, F.; Zazzi, M.; Saladini, F. The HIV-1 reverse transcriptase E138A natural polymorphism decreases the genetic barrier to resistance to etravirine in vitro. J. Antimicrob. Chemother. 2019, 74, 607–613. [Google Scholar] [CrossRef] [PubMed]
  53. Bavaro, D.F.; Di Carlo, D.; Rossetti, B.; Bruzzone, B.; Vicenti, I.; Pontali, E.; Zoncada, A.; Lombardi, F.; Di Giambenedetto, S.; Borghi, V.; et al. Pretreatment HIV drug resistance and treatment failure in non-Italian HIV-1-infected patients enrolled in ARCA. Antivir. Ther. 2020. [Google Scholar] [CrossRef]
  54. Sluis-Cremer, N.; Jordan, M.R.; Huber, K.; Wallis, C.L.; Bertagnolio, S.; Mellors, J.W.; Parkin, N.T.; Harrigan, P.R. E138A in HIV-1 reverse transcriptase is more common in subtype C than B: Implications for rilpivirine use in resource-limited settings. Antiviral Res. 2014, 107, 31–34. [Google Scholar] [CrossRef] [Green Version]
  55. Rhee, S.Y.; Blanco, J.L.; Jordan, M.R.; Taylor, J.; Lemey, P.; Varghese, V.; Hamers, R.L.; Bertagnolio, S.; de Wit, T.F.; Aghokeng, A.F.; et al. Correction: Geographic and Temporal Trends in the Molecular Epidemiology and Genetic Mechanisms of Transmitted HIV-1 Drug Resistance: An Individual-Patient- and Sequence-Level Meta-Analysis. PLoS Med. 2015, 12, e1001845. [Google Scholar] [CrossRef] [Green Version]
  56. Machnowska, P.; Meixenberger, K.; Schmidt, D.; Jessen, H.; Hillenbrand, H.; Gunsenheimer-Bartmeyer, B.; Hamouda, O.; Kücherer, C.; Bannert, N.; German HIV-1 Seroconverter Study Group. Prevalence and persistence of transmitted drug resistance mutations in the German HIV-1 Seroconverter Study Cohort. PLoS ONE 2019, 14, e0209605. [Google Scholar] [CrossRef]
Viruses 12 00791 g001 550
Figure 1. Clusters’ population by HIV-TRACE. The coloured geometric shapes identify the nationality and risk factors of the 3499 individuals involved in the study. Shapes without connections represent individuals not involved in molecular transmission clusters. Shapes connected with lines represent individuals involved in molecular transmission clusters.
Figure 1. Clusters’ population by HIV-TRACE. The coloured geometric shapes identify the nationality and risk factors of the 3499 individuals involved in the study. Shapes without connections represent individuals not involved in molecular transmission clusters. Shapes connected with lines represent individuals involved in molecular transmission clusters.
Viruses 12 00791 g001
Table
Table 1. Patient’s characteristics and factors associated with HIV-1 molecular transmission clusters.
Table 1. Patient’s characteristics and factors associated with HIV-1 molecular transmission clusters.
VariablesOverallOut of ClusterIn Clusterp-Value aAdjusted Model b
n = 34992773 (79.3%)726 (20.7%)OR(95% CI)p-Value
Male gender, n (%)2872 (82.1%)2187 (78.9%)685 (94.3%)<0.001---
Age, years, median (IQR)37 (30–45)38 (31–46)32 (27–40)<0.0010.650.59–0.72<0.001
Mode of HIV transmission, n (%)
F heterosexual553 (15.8%)513 (18.5%)40 (5.5%)<0.0011.00--
F IVDU32 (0.9%)31 (1.1%)1 (0.1%) 0.490.06–3.800.497
M heterosexual713 (20.4%)628 (22.7%)85 (11.7%) 1.821.19–2.780.006
M IVDU161 (4.6%)145 (5.2%)16 (2.2%) 1.520.80–2.900.204
MSM1789 (51.1%)1247 (45.0%)550 (75.8%) 3.462.39–5.03<0.001
Other/unknown251 (7.2%)209 (7.5%)34 (4.7%) 2.731.66–4.48<0.001
Nation of birth, n (%)
Italy2804 (80.1%)2160 (77.8%)644 (88.7%)<0.0011.00--
Africa219 (6.3%)212 (7.7%)7 (0.9%) 0.180.08–0.39<0.001
Central and South America241 (6.9%)201 (7.3%)40 (5.5%) 0.490.33–0.710.001
Europe187 (5.3%)159 (5.7%)28 (3.9%) 0.290.08–0.970.045
Asia38 (1.1%)35 (1.3%)3 (0.4%) 0.620.40–0.970.035
Other10 (0.3%)6 (0.2%)4 (0.6%) 2.610.62–10.970.189
Education, n (%)
Primary school169 (4.8%)158 (5.7%)11 (1.5%)<0.0010.870.45–1.710.691
Secondary school585 (16.7%)505 (18.2%)80 (11.0%) 0.910.68–1.210.518
College/University1762 (50.4%)1329 (47.9%)433 (59.6%) 1.00-
Unknown983 (28.1%)781 (28.2%)202 (27.8%) 1.030.83–1.280.773
Employment, n (%)
Employed1476 (42.2%)1148 (41.4%)328 (45.2%)<0.0011.00--
Unemployed461 (13.2%)389 (14.0%)72 (9.9%) 0.910.67–1.250.565
Self-employed526 (15.0%)413 (14.9%)113 (15.6) 0.970.75–1.260.840
Student146 (4.2%)93 (3.4%)53 (7.3%) 0.830.56–1.240.360
Housewife94 (2.7%)88 (3.2%)6 (0.8%) 1.180.47–2.940.723
Other278 (7.9%)244 (8.8%)34 (4.7%) 0.690.46–1.050.083
Unknown518 (14.8%)398 (14.3%)120 (16.5%) 0.990.75–1.300.929
HIV RNA, copies/mL, n (%)
<1000122 (3.5%)99 (3.6%)23 (3.2%)0.0050.660.39–1.120.127
1000–10,000559 (16.0%)445 (16.1%)114 (15.7%) 0.810.61–1.080.147
10,000–100,0001470 (42.0%)1126 (40.6%)344 (47.4%) 0.940.75–1.170.562
>100,0001118 (32.0%)905 (32.6%)213 (29.3%) 1.00-
Unknown230 (6.6%)198 (7.1%)32 (4.4%) 0.760.38–1.550.451
CD4, cells/mm3, n (%)
≤200791 (22.6%)721 (26.0%)70 (9.6%)<0.0011.00--
201–5001485 (42.4%)1165 (42.0%)320 (44.1%) 2.221.64–2.99<0.001
>5001003 (28.7%)703 (25.4%)300 (41.3%) 3.012.20–4.13<0.001
Unknown220 (6.3%)184 (6.6%)36 (5.0%) 1.900.93–3.900.078
Year of diagnosis,
median (IQR)		2011
(2008–2014)		2011
(2007–2014)		2012
(2009–2014)		<0.0011.091.06–1.11<0.001
Subtype, n (%)
A1104 (3.0%)85 (3.1%)19 (2.6%)<0.001---
B2556 (73.1%)2038 (73.5%)518 (71.4%) ---
C148 (4.2%)119 (4.3%)29 (4.0%) ---
CRF02_AG187 (5.3%)141 (5.1%)46 (6.3%) ---
CRF60_BC64 (1.8%)12 (0.4%)52 (7.2%) ---
F1179 (5.1%)157 (5.7%)22 (3.0%) ---
Other261 (7.5%)221 (8.0%)40(5.5%) ---
Transmitted drug resistance, n (%)
Any drug class501 (14.3)420 (15.1)81 (11.2)0.006---
NNRTI335 (9.6)281 (10.1)54 (7.4)0.028---
NRTI150 (4.3)134 (4.8)16 (2.2)0.002---
PI66 (1.9)53 (1.9)13 (1.8)0.832---
a By Mann–Whitney test (for quantitative variables) and χ2 test or Fisher’s exact test (for categorical variables), as appropriate. p-values <0.05 were considered statistically significant and were reported in bold. b Adjusted for: Sex, age, mode of HIV transmission, nation of birth, education, employment, plasma HIV-RNA, CD4 cell count, year of diagnosis. Variables that were significant in univariable analysis (p < 0.05) were considered for the multivariable model. F: female; IVDU: intravenous drug user; M: male; MSM: men who have sex with men; NNRTI: non-nucleoside reverse transcriptase inhibitor; NRTI: nucleos(t)ide reverse transcriptase inhibitor; PI: protease inhibitor.
Table
Table 2. Characteristics of molecular transmission clusters found among 3499 drug-naïve individuals living in Italy.
Table 2. Characteristics of molecular transmission clusters found among 3499 drug-naïve individuals living in Italy.
Cluster’s CharacteristicsSmall ClustersMedium Clusters
 
n = 36		Large Clusters
 
n = 6
OverallClusters of 2 PersonsClusters of 3 Persons
n = 186n = 156n = 30
Sampling interval (min–max)2–3--4–910–52
Age interval (min–max)19–6519–6520–5618–5318–58
Mode of HIV transmission, n (%)
Only MSM101 (54.3%)82 (52.5%)19 (63.3%)18 (50.0%)-
Only heterosexual (F+M)13 (7.0%)13 (8.3%)---
Only heterosexual (only F)4 (2.2%)4 (2.6%)---
Only heterosexual (only M)5 (2.7%)5 (3.2%)---
Only IVDU (F+M)1 (0.5%)1 (0.6%)---
Only IVDU (only M)2 (1.1%)2 (1.3%)---
Mixed (F+M)18 (9.7%)9 (5.8%)4 (13.3%)4 (11.1%)1 (16.7%)
Mixed (only M)42 (22.5%)40 (25.5%)7 (23.3%)14 (38.9%)5 (83.3%)
Nationality, n (%)
Only Italians137 73.7%)122 (78.2%)15 (50.0%)24 (66.7%)2 (33.3%)
Italians and migrants43 (23.1%)28 (18.0%)15 (50.0%)12 (33.3%)4 (66.7%)
Only migrants6 (3.2%)6 (3.8%)---
Subtype, n (%)
A13 (1.6%)2 (1.2%)1 (3.3%)2 (5.6%)-
B149 (80.1%)125 (80.1%)24 (80.0%)25 (69.3%)3 (50.0%)
C5 (2.7%)4 (2.6%)1 (3.3%)2 (5.6%)1 (16.7%)
CRF02_AG5 (2.7%)3 (1.9%)2 (6.7%)5 (13.9%)1 (16.7%)
CRF60_BC----1 (16.7%)
F111 (5.9%)11 (7.1%)---
Other13 (7.0%)11 (7.1%)2 (6.7%)2 (5.6%)-
Resistance, n (%)
None150 (80.7)123 (78.8)27 (90.0)31 (86.1)5 (83.3)
Any drug class36 (19.3)33 (21.2)3 (10.0)5 (13.9)1 (16.7)
NNRTI24 (12.9)22 (14.1)2 (6.7)3 (8.3)1 (16.7)
NRTI6 (3.2)5 (3.2)1 (3.3)1 (2.8)0 (0.0)
PI7 (3.8)7 (4.5)0 (0.0)1 (2.8)0 (0.0)
IVDU: intravenous drug user; MSM: men who have sex with men; NNRTI: non-nucleoside reverse transcriptase inhibitor; NRTI: nucleos(t)ide reverse transcriptase inhibitor; PI: protease inhibitor.
Table
Table 3. Clusters characteristics.
Table 3. Clusters characteristics.
Medium MTCs (4–9 Sequences)
IDSubtypeCluster Size a
n		Sampling Interval YearHIV Diagnosis Interval YearMigrants
n (%)		Migrants’ Nationality
n		Age, Years
Median (IQR)		Risk Factor
n		Genetic Distance
Mean (SE)		Resistance
1B42006–20072006–20070 (0.0)-46 (44–46)4 MSM0.008 (0.002)No
2B42007–20122004–20120 (0.0)-46 (42–50)1 F Het, 2 M Het, 1 M Unk0.010 (0.002)No
3B42008–20142008–20140 (0.0)-52 (50–53)4 MSM0.008 (0.002)No
4B42008–20162008–20160 (0.0)-39 (37–41)3 MSM, 1 M Unk0.006 (0.002)No
5B42009–20122009–20120 (0.0)-54 (51–56)3 MSM, 1 F Het0.009 (0.002)No
6B42009–20162009–20160 (0.0)-36 (33–41)3 MSM, 1 M Unk0.009 (0.002)Yes (n = 4; RT: E138A)
7B42010–201120101 (25.0)1 NA50 (46–52)4 MSM0.004 (0.002)No
8B42011–20132011–20130 (0.0)-37 (35–38)4 MSM0.008 (0.002)No
9B42011–20172011–20170 (0.0)-53 (48–56)4 MSM0.007 (0.002)No
10B42013–20142010–20140 (0.0)-41 (34–48)4 MSM0.006 (0.002)No
11B42013–20182013–20181 (25.0)1 EE41 (37–43)3 MSM, 1 M Het0.008 (0.002)No
12B420162015–20160 (0.0)-38 (28–49)4 MSM0.004 (0.001)Yes (n = 1; PR: I85V)
13B52007–20122007–20111 (20.0)1 NA38 (37–45)5 MSM0.006 (0.001)No
14B52008–20142008–20140 (0.0)-36 (36–37)5 MSM0.011 (0.002)No
15B52010–20122009–20122 (40.0)1 CA, 1 EE41 (37–43)5 MSM0.008 (0.002)No
16B52011–20172008–20170 (0.0)-38 (37–45)4 MSM, 1 M Het0.009 (0.002)No
17B52011–20172011–20170 (0.0)-33 (32–35)5 MSM0.009 (0.002)Yes (n = 5; RT: K101E)
18B52013–20182013–20181 (20.0)1 SA38 (36–43)3 MSM, 2 M Het0.011 (0.002)No
19B62007–20112006–20110 (0.0)-41 (36–48)6 MSM0.010 (0.002)No
20B62007–20112007–20110 (0.0)-42 (37–45)5 MSM, 1 M Unk0.011 (0.002)No
21B62010–20142010–20140 (0.0)-37 (35–43)5 MSM, 1 M Het0.011 (0.002)No
22B72006–20162006–20091 (14.3)1 SA43 (39–48)7 MSM0.010 (0.002)No
23B72009–20142008–20142 (28.6)2 SA51 (40–53)7 MSM0.011 (0.002)No
24B92010–20162010–20140 (0.0)-40 (34–52)7 MSM, 1 M Het, 1 M Unk0.009 (0.002)No
25B92017–20182017–20180 (0.0)-34 (31–41)8 MSM, 1 M Unk0.006 (0.001)No
26A142013–20152013–20150 (0.0)-47 (45–50)4 MSM0.010 (0.002)No
27A182014–20162013–20161 (12.5)1 WE37 (33–39)8 MSM0.008 (0.002)No
28C42007–20112006–20111 (25.0)1 SA46 (44–47)3 MSM, 1 M Het0.007 (0.002)No
29C42012–20162012–20160 (0.0)-50 (45–52)3 M Het, 1 M MSM0.009 (0.002)No
30CRF02_AG42009–20132009–20120 (0.0)-35 (35–37)4 MSM0.007 (0.002)No
31CRF02_AG42009–20132009–20130 (0.0)-38 (38–40)2 MSM, 1 F Het, 1 M Unk0.004 (0.001)No
32CRF02_AG4201420141 (25.0)1 EE34 (30–39)3 MSM, 1 M Het0.009 (0.002)Yes (n = 1; RT: E138A)
33CRF02_AG52013–20142012–20140 (0.0)-34 (33–35)4 MSM, 1 M Unk0.006 (0.002)No
34CRF02_AG72010–20172010–20160 (0.0)-41 (36–51)6 MSM, 1 M Het0.005 (0.001)No
35CRF12_BF52014–20152013–20151 (20.0)1 SA34 (27–34)2 F Het; 2 M Het, 1 M Unk0.010 (0.002)Yes (n = 3; RT: A62V)
36CRF20_BG72013–20172013–20171 (14.3)1 Aus33 (32–37)7 MSM0.005 (0.001)No
Large MTCs (≥10 Sequences)
1B142008–20162008–20160 (0.0)-40 (34–49)11 MSM, 2 M Het, 1 M IVDU0.010 (0.001)No
2B192007–20152007–20152 (10.5)1 EE, 1 SA37 (32–43)18 MSM, 1 M Unk0.016 (0.002)No
3B352009–20172009–20172 (5.7)1 WE, 1 SA34 (33–43)33 MSM, 1 M Het, 1 M Unk0.013 (0.002)No
4C102011–20162011–20160 (0.0)-39 (35–42)8 MSM, 1 M Het, 1 M Unk0.009 (0.001)No
5CRF02_AG102006–20162006–20161 (10.0)1 CA42 (35–46)9 MSM, 1 M Het0.011 (0.002)No
6CRF60_BC522008–20182008–20185 (9.6)4 EE, 1 Unk34 (31–37)41 MSM, 2 F Het, 3 M Het, 6 M Unk0.012 (0.001)Yes (n = 1; RT: E138A)
a Number of individuals involved in a specific MTC. Het: heterosexual; MSM: men who have sex with men; Unk: unknown; PR: protease; RT: reverse transcriptase. Aus: Australian, CA: Central American; EE: East European; F: female; M: male; NA: North American; SA: South American; WE: West European.

Share and Cite

MDPI and ACS Style

Fabeni, L.; Santoro, M.M.; Lorenzini, P.; Rusconi, S.; Gianotti, N.; Costantini, A.; Sarmati, L.; Antinori, A.; Ceccherini-Silberstein, F.; d’Arminio Monforte, A.; et al. Evaluation of HIV Transmission Clusters among Natives and Foreigners Living in Italy. Viruses 2020, 12, 791. https://doi.org/10.3390/v12080791

AMA Style

Fabeni L, Santoro MM, Lorenzini P, Rusconi S, Gianotti N, Costantini A, Sarmati L, Antinori A, Ceccherini-Silberstein F, d’Arminio Monforte A, et al. Evaluation of HIV Transmission Clusters among Natives and Foreigners Living in Italy. Viruses. 2020; 12(8):791. https://doi.org/10.3390/v12080791

Chicago/Turabian Style

Fabeni, Lavinia, Maria Mercedes Santoro, Patrizia Lorenzini, Stefano Rusconi, Nicola Gianotti, Andrea Costantini, Loredana Sarmati, Andrea Antinori, Francesca Ceccherini-Silberstein, Antonella d’Arminio Monforte, and et al. 2020. "Evaluation of HIV Transmission Clusters among Natives and Foreigners Living in Italy" Viruses 12, no. 8: 791. https://doi.org/10.3390/v12080791

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop