Hostname: page-component-586b7cd67f-dlnhk Total loading time: 0 Render date: 2024-11-24T09:51:45.254Z Has data issue: false hasContentIssue false

Molecular epidemiology of HIV-1 infection in immigrant population in northern Italy

Published online by Cambridge University Press:  05 February 2020

Caterina Sagnelli
Affiliation:
Department of Mental Health and Public Medicine, Campania University Luigi Vanvitelli, Naples, Italy
Caterina Uberti-Foppa
Affiliation:
Department of Infectious Diseases, Vita-Salute University, San Raffaele Scientific Institute, Milan, Italy
Sabrina Bagaglio
Affiliation:
Department of Infectious Diseases, Vita-Salute University, San Raffaele Scientific Institute, Milan, Italy
Eleonora Cella
Affiliation:
Medical Statistics and molecular Epidemiology Unit, Campus Bio-Medico University, Rome, Italy
Vittoria Scolamacchia
Affiliation:
Medical Statistics and molecular Epidemiology Unit, Campus Bio-Medico University, Rome, Italy
Hamid Hasson
Affiliation:
Department of Infectious Diseases, Vita-Salute University, San Raffaele Scientific Institute, Milan, Italy
Stefania Salpietro
Affiliation:
Department of Infectious Diseases, Vita-Salute University, San Raffaele Scientific Institute, Milan, Italy
Emanuela Messina
Affiliation:
Department of Infectious Diseases, Vita-Salute University, San Raffaele Scientific Institute, Milan, Italy
Giulia Morsica
Affiliation:
Department of Infectious Diseases, Vita-Salute University, San Raffaele Scientific Institute, Milan, Italy
Silvia Angeletti
Affiliation:
Laboratory clinical Science, Campus Bio-Medico University, Rome, Italy
Massimo Ciccozzi
Affiliation:
Medical Statistics and molecular Epidemiology Unit, Campus Bio-Medico University, Rome, Italy
Adriano Lazzarin
Affiliation:
Department of Infectious Diseases, Vita-Salute University, San Raffaele Scientific Institute, Milan, Italy
Evangelista Sagnelli*
Affiliation:
Department of Mental Health and Public Medicine, Campania University Luigi Vanvitelli, Naples, Italy
*
Author for correspondence: Evangelista Sagnelli, E-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Human immunodeficiency virus-1 (HIV-1) is characterised by a vast genetic diversity classified into distinct phylogenetic strains and recombinant forms. We describe the HIV-1 molecular epidemiology and evolution of 129 consecutive HIV-1 positive migrants living in Milan (northern Italy). Polymerase gene sequences of 116 HIV-1 subtype-B positive patients were aligned with HIV-1 reference sequences (https://www.ncbi.nlm.nih.gov/) by using MAFFT alignment and edited by using Bioedit software. A maximum likelihood (ML) phylogenetic tree was performed by MEGA7 and was visualised by using FigTree v1.4.3. Of 129 migrants, 35 were born in Europe (28 in Eastern Europe), 70 in the Americas (67 in South America), 15 in Africa and nine in Asia; 76.4% were men who have sex with men (MSM). The serotype HIV-1-B prevailed (89.9%), followed by -C, -F1, -D and -A. Compared with 116 HIV-B patients, the 13 with HIV-non-B showed lower Nadir of CD4+ cell/mmc (P = 0.043), more frequently had sub Saharan origin (38.5 vs. 1.72%, P = 0.0001) and less frequently were MSM (40 vs. 74.5%, P = 0.02). The ML phylogenetic tree of the 116 HIV-1 subtype-B positive patients showed 13 statistically supported nodes (bootstrap > 70%). Most of the sequences included in these nodes have been isolated from male patients from the Americas and the most common risk factor was MSM. The low number of HIV-1 non-B subtype patients did not allow to perform this analysis. These results suggest a shift of HIV-1 prevention projects' focus and a continuous monitoring of HIV-1 molecular epidemiology among entry populations. Prevention efforts based on HIV molecular epidemiology may improve public health surveillance setting.

Type
Original Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s) 2020

Introduction

Two major types of human immunodeficiency virus (HIV) have been identified, HIV type 1 and HIV type 2 (HIV-2), with HIV-2, differing genetically from HIV-1 by nearly 55%. The rapid worldwide spread of HIV-1 has been favoured by its enormous genetic variability and rapid evolution, making the virus highly adaptable to new hosts. HIV-1 sequences differ by up to 10% in single individuals [Reference Korber1] due to its error-prone high replication and substitution rate consequent to the lack of proof-reading activity of the reverse transcriptase enzyme [Reference Seillier-Moiseiwitsch2, Reference Peeters and Sharp3]. HIV is classified into four distinct groups: M (major), O (outlier), N (non-M/non-O) and P [Reference Van Heuverswyn4Reference D'Arc8]. Group M viruses are responsible for the HIV-1 global pandemic, further classified in distinct subtypes (A-D, F\H, J and K), sub-subtypes (A1, A2, F1 and F2) and circulating recombinant forms (CRFs) [‘www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html'].

Subtype B is the dominant one in high-income countries with sexual transmissions among men who have sex with men (MSM), injection drug users (IDUs) and transfusions in hemophiliacs, with non-B subtypes/CRFs rising steadily among heterosexuals [Reference Paraskevis9Reference Lai15].

In Italy, the HIV-1 epidemic is mainly sustained by the subtype B genetic form, although other and novel subtypes and CRFs have been identified [Reference Balotta16Reference Buonaguro19].

The percentage of infection with HIV-1 non-B subtypes ranged from 2.4% to 19.4% in different Italian studies reaching 63.0% among migrants [Reference Balotta20Reference Ciccozzi22], with an overtime significant increase [Reference Ciccozzi17, Reference Lai23] paralleling the increase in migratory flows from geographic areas at high prevalence of HIV-1 infection [Reference Ciccozzi24Reference Abecasis26]. Accordingly, the increased prevalence in HIV-1 non-B subtypes in Italy has been associated with the increase of the immigration flow [Reference Abecasis26, Reference Hemelaar27]. Indeed, the proportion of migrants among the new diagnoses of HIV infection has increased from 11% in 1992 to 35.8% in 2016 [‘http://old.iss.it/binary/publ/cont/COA.pdf’].

Ciccozzi et al. reported a 90.7% prevalence of HIV-1 subtype B in 215 HIV-1 pol gene sequences collected between 1992 and 2010, in Italy. Three main clusters were detected which roots dated to 1987. Most of the observed viral gene flow events occurred from heterosexual to intravenous drug users. Phylogenetic and molecular clock analysis showed an early HIV-1 subtype B introduction in the mid-1980 and dissemination within local risk-specific clusters [Reference Ciccozzi28].

Lo Presti et al. analysed 883 HIV-1 subtype pol gene sequences obtained from patients living in the geographical area of Bergamo (a chief town in northern Italy) declaring different transmission risks. In this study, the 25% of the observed gene flow was from people living in the north valleys to lowland and 40.5% from a heterosexual risk group to injecting drug users. The latter were identified as the central link and the mercenary sex as the most common route of transmission and gene flow between this group and both heterosexual and homosexual individuals [Reference Presti A29].

The history of human being is history of migration which have determined the current ethnic composition of countries. At present, most of migrants come to north-western nations from their native countries in the south of the world for a variety of reasons, including family re-unification, desire for economic prosperity, natural disaster, dictatorial policy or escaping conflict. Due to the socioeconomic and political crises or tribal war in Sub-Saharan Africa counties, Western countries have become lands of immigration from this geographic area at high HIV endemicity. The migrant population living in Italy are prevalently male, young, sexually active [Reference Mabey and Mayaud30] and with broken family ties. They often live in overcrowded residences, are not socially integrated and have limited access to healthcare services [Reference Tiedje31].

A high proportion of migrants living in European countries was born in geographic areas at high HIV endemicity and, consequently, is infected with HIV more frequently than the subjects born in the host country. The European Centre for Disease Prevention and Control reported that from 2007 to 2011 about 40% of all new cases of HIV infection diagnosed in Europe were registered in foreign patients [32], who nevertheless constituted about the 10% of the total resident population. In Italy, the median incidence of HIV infection at the end of 2013 was 19.1 per 100 000 for migrants and 4.9 per 100 000 for Italian-born subjects, reflecting the frequent arrival of migrants from sub Saharan countries [33].

Compared with the native populations of Italy and of other European countries, HIV infection is diagnosed in a more advanced stage in migrant populations [33Reference Sulis37], as was also confirmed in Italy by a recent survey [34], due to various psychological and bureaucratic barriers delaying HIV testing and antiretroviral therapy administration [Reference Fakoya38, Reference Alvarez-del Arco39].

The study of HIV-1 epidemic diversity provides a useful tool through which we can understand the history of the pandemic and monitor its spread and epidemic growth over time. The objective of this study was to analyse in detail the epidemiology, characteristics and evolution of HIV in 129 migrants residing in Milan (northern Italy). Since 116 of the 129 subjects were infected with HIV-1B and only 13 with HIV-1-non B, the phylogenetic and evolutionary analysis was focused on the HIV-1 subtype B sequences only.

Material and methods

Study population

One-hundred twenty-nine consecutive anti-HIV-1 positive adult migrants were enrolled from 2008 to 2018 at the Department of Infectious Diseases, Vita-Salute University, San Raffaele Scientific Institute, Milan, Italy.

The following criteria for including a patient into the study should have be fulfilled: (1) availability of a serum sample obtained at the first observation, stored frozen at −70 °C, never towed and used only for this investigation and (2) availability of registered demographic data, self-reported risk factors for the acquisition of HIV infection, geographic area of origin and biochemical virologic and immunologic data obtained at the time of first observation when all patients were naïve for anti-retroviral therapy. All the data mentioned above have been filled in a pre-coded questionnaire.

Criteria for exclusion were an ongoing or previous antiretroviral therapy and a not even slight compliance with the criteria required for inclusion (see above at points 1 and 2).

Serological determinations

Serum HBsAg was sought by a commercial immunoenzymatic assay (Abbott Laboratories, North Chicago, IL, USA) and anti-HCV antibody by a 3rd generation commercial immunoenzymatic assay (Ortho Diagnostic Systems, Neckargemund, Germany). Antibodies to HIV 1 and 2 were sought using a commercial ELISA (Abbott Lab., North Chicago, Ill, USA) and positive results were always confirmed by Western blot analysis (Genelabs Diagnostics, Science Park Drive, Singapore), according to the Italian law.

Lymphocyte subsets (CD4+ and CD8+) were evaluated by flow cytofluorimetry using monoclonal antibodies and a fluorescence-activated cell sorter scan (Becton Dickinson, Mountain View, USA). Liver function tests, triglyceride and cholesterol levels were determined according to routine methods. The body mass index (BMI) was determined by standard procedures.

Statistical analysis

Continuous variables were not normally distributed and were summarised as median and interquartile range, and categorical variables as absolute and relative frequencies. The Mann–Whitney and χ 2 non parametric test have been used, the P value less than 5% was considered statistically significant. The statistical analysis was performed using Stata software version 14.1 (StataCorp Texas 77845 USA).

Molecular epidemiology and phylogenetic procedures

For the phylogenic analyses only HIV-1 subtype B were considered (116 sequences on 129), because mostly represented.

Two different datasets were built to investigate the phylogenetic relationships and the genetic variability of the HIV-1 virus pol gene.

The first dataset was built using 116 HIV-1 virus pol gene subtype-B sequences plus 76 reference sequences, downloaded from NCBI [‘http://www.ncbi.nlm.nih.gov/genbank/’]. The second dataset were built using 33 HIV-1 virus pol gene subtype-B sequences. This dataset included sequences with patients known information and included in statistically supported clusters. These two datasets were used to perform, respectively, maximum likelihood (ML) and Bayesian dated trees, evolutionary and phylodynamic analyses.

The reference sequences were selected based on the following criteria: (1) sequences already published in peer-reviewed journals; (2) known sampling date and location and (3) all the available sequences as always described [Reference Sagnelli40]. All sequences were aligned using MAFFT [Reference Katoh41] and manual editing was performed with Bioedit, removing gaps and cutting to identical sequence lengths.

MEGA7 was used to select the simplest evolutionary model that adequately fitted the sequence data for the datasets by using the ‘Models’ tool. The phylogenetic signal was tested with TreePuzzle [Reference Strimmer and von Haeseler42, Reference Salemi43], DAMBE [Reference Xia44] and MEGA7 [Reference Kumar45] software.

The ML tree was inferred on the first dataset; the statistical robustness and reliability of the branching order was confirmed with the bootstrap analysis (bootstrap values >70%) [Reference Sagnelli46, Reference Sagnelli47, ‘http://beast.bio.ed.ac.uk’].

The evolutionary rate of the HIV-1 virus pol gene subtype-B (second dataset) was estimated by calibrating a molecular clock using known sequence sampling times with the Bayesian Markov Chain Monte Carlo (MCMC) method implemented in BEAST v. 1.10.1 [Reference Drummond and Rambaut48, Reference Kass and Raftey49]. To investigate the demographic history, independent MCMC runs were carried out enforcing both a strict and relaxed clock with an uncorrelated log normal rate distribution and one of the following coalescent priors: constant population size, exponential growth, non-parametric smooth skyride plot Gaussian Markov Random Field and non-parametric Bayesian skyline plot (BSP) [Reference Drummond and Rambaut48Reference Baele and Lemey55]. Chains were conducted for at least 50 × 106 generations and sampled every 5000 steps for each molecular clock model. Convergence of the MCMC was assessed by calculating the ESS for each parameter.

Results

The initial characteristics of the 129 immigrants with HIV infection are reported in Tables 1 and 2. They aged 35.12 ± 8.65 years, have been followed-up for 9.59 ± 4.79 years and were prevalently males (82.94%). Reliable information on the main risk factor for acquiring HIV infection was available only for 110 patients; risky heterosexuality was declared by 21 (19.09%) of which two being sexual partner of an HIV-positive subject, male homosexuality by 84 (76.36%) of which four bisexual, drug addiction by four (3.64%) and being born of an anti-HIV positive mother by one (0.91%) (Table 1).

Table 1. Demographics and epidemiological characteristics of all the 129 HIV-1 positive immigrants enrolled, and according to HIV-1 serotype

Table 2. Biochemical and immunovirological characteristics of the 129 HIV-1 positive immigrants, and according to HIV-1 serotype.

Of the 129 migrants, 35 were born in Europe (28 in Eastern Europe), 70 in the Americas (67 in South America, of whom 21 in Brazil and 21 in Peru) and 15 in Africa (of whom seven born in sub Saharan Africa and nine in Asia) (Table 1).

The serotype HIV-1 B prevailed (89.9%), followed by serotypes HIV-1 C (4.6%), HIV-1 F1 (3.9%), HIV-1 D (0.78) and HIV-1 serotypes A (0.78%) (Table 1).

The 129 HIV-1 positive immigrants showed a quite good immunological condition at enrolment since the mean CD4+ cells per mm3 count was 472.00 ± 257.00 (M ± s.d.), the percentage of CD4% 21.46 ± 9.3 and the nadir median value through the observation 326.58 ± 197.42 CD4+ cell/mmc (Table 2). The mean HIV-RNA, cps/ml was 79 013.30 ± 171 001.37, indicating an active HIV replication (Table 2). In addition, the mean values of initial serum of liver enzymes, metabolic indexes and creatininemia and the infrequent coinfection with HBV (three cases) and with HCV (five cases) testify that the 129 patients as a whole had a satisfactory clinical condition (Table 2).

Compared with the 116 patients with serotype HIV-B, the 13 with serotype HIV-non-B showed a significantly lower Nadir of CD4+ cell/mmc (P = 0.043), more frequently were born in Sub-Saharan Africa (38.46 vs. 1.72%, P = 0.0001) and less frequently were MSM (40.0 vs. 74.5%, P = 0.02); no other significant demographic, epidemiological and laboratory difference was observed between these two subgroups of patients (Tables 1 and 2).

Phylogenetic analysis

The best evolutionary model was K2 + G + I for the first and HKY + G for the second dataset. The phylogenetic noise for the datasets was investigated by means of likelihood mapping. The percentage of dots falling in the middle area of the triangles was 8.6% for the first dataset and 6.7% for the second dataset, as none of the dataset showed more than 30% of noise, contained sufficient phylogenetic signals.

The percentage of the Parsimony-Info sites ranged from 42% (the first) to 66.7% (second dataset); the percentage of the constant sites ranged from 23.84% (second dataset) to 46% (the first dataset).

The phylogenetic signal analysis using a transition/transversion vs. divergence graph and the Xia's test (P < 0.001) did not show evidence for substitution saturation. This indicated that enough signal for phylogenetic inference existed (data not shown).

All the analyses performed showed a sufficient phylogenetic signal for further phylogenetic analysis.

ML phylogenetic tree of the first dataset was reported in Figure 1. Many statistically supported clusters have been found. The first cluster showed a correlation between a male patient born in Brazil tested for HIV-1 virus in 2008 and a heterosexual female patient born in Peru tested for HIV-1 virus in 2011. The second cluster underlined a relationship among a female patient born in Albania tested for HIV-1 virus in 2011, a female patient born in Albania, arrived in Italy in 2011, tested for HIV-1 virus in 2013 and a male patient born in Albania tested for HIV-virus in 2016. The third cluster is divided in two sub-clusters. The first sub-cluster showed a correlation between a homosexual male patient born in Montenegro tested for HIV-1 virus in 2008 and a homosexual male patient born in Germany tested for HIV-1 virus in 2011. The second sub-cluster showed a relationship among a homosexual male patient born in El Salvador tested for HIV-1 virus in 2012, a homosexual male patient born in Ecuador tested for HIV-1 virus in 2016 and a homosexual male patient born in Peru tested for HIV-1 virus in 2015. The fourth cluster observed a correlation between a male patient born in Spain tested for HIV-virus in 2016 and a male patient born in Peru, arrived in Italy in 2016, tested for HIV-1 virus in 2016. The fifth cluster showed a relationship between a homosexual male patient born in Ireland, tested for HIV-1 virus in 2008 and a homosexual male patient born in Philippines tested for HIV-1 virus in 2016. The sixth cluster evidenced a correlation between a homosexual male patient born in Brazil, arrived in Italy in 2009, tested for HIV-1 virus in 2012 and a homosexual male patient born in USA, texted for HIV-1 virus in 2016. The seventh cluster correlated a drug addicted female patient born in Nigeria tested for HIV-1 in 2013 and a HIV-1 virus subtype-B reference sequences isolated from Peru in 2007. The eighth cluster showed a relationship between a drug addicted female patient born in Nigeria tested for HIV-1 virus in 2013 and a male patient born in Tunisia tested for HIV-1 virus in 2014. The ninth cluster correlated the sequence 92726_RT and the sequence 92754_RT, two sequences with no information found. The tenth cluster showed a correlation between a heterosexual female patient born in Brazil tested for HIV-1 virus in 2014 and a homosexual male patient born in Colombia, arrived in Italy in 2008, tested for HIV-1 virus in 2009. The eleventh cluster supported a relationship between a homosexual male patient born in Turkey tested for HIV-1 in 2014, a homosexual male patient born in Ecuador tested for HIV-1 in 2009 and a homosexual male patient born in Peru for HIV-1 in 2011. The twelfth cluster showed a correlation among a homosexual male patient born in Peru tested for HIV-1 virus in 2008, the sequence 92779_RT with no information found and a homosexual male patient born in Philippines tested for HIV-1 in 2015.

Fig. 1. Phylogenetic relationship of the HIV-1 subtype-B sequences with the subtype-B reference sequences downloaded from the NCBI sequence database (https://www.ncbi.nlm.nih.gov/). The asterisk along a branch represents significant statistical support for the clade subtending that branch (bootstrap support 70%). The collapsed cluster in the figure were what not statistically significant.

The thirteenth cluster supported a relationship between a homosexual male patient born in Argentine tested for HIV-1 virus in 2010 and a homosexual male patient born in Bulgaria tested for HIV-1 virus in 2010. This cluster showed also a sub-cluster composed of a homosexual male patient born in Brazil tested for HIV-1 in 2012 and a homosexual male patient born in Spain tested for HIV-1 in 2014. The fourteenth cluster correlated a male patient born in Israel tested for HIV-1 virus in 2008 and a male patient born in Hungary and tested for HIV-1 in 2008.

Bayesian reconstruction of the time-scaled phylogeny

Figure 2 shows the maximum clade credibility (MCC) tree of second dataset.

Fig. 2. MCC tree of the second dataset. Branches are scaled in time. Significant posterior probability support (pp ≥ 0.9) as indicated by an asterisk. Clusters are highlighted.

The BSP as a demographic model with a relaxed molecular clock was selected as the most appropriate to describe the evolutionary history of HIV-1 B pol sequences (lnBF > 6). Molecular clock calibration estimated the evolutionary rate of the HIV-1 B pol sequences at 2.17 × 10−3 substitutions site per year (95% HPD 1.37 × 10−3–2.9 × 10−3).

There were ten statistically supported clusters. Cluster A includes three sequences, isolated from male homosexual patients isolated in Turkey Ecuador and Peru. The most recent common ancestor (tMRCA) corresponding to 1996, inside this cluster two sequences formed a group dated back to January 2008.

Cluster B dated back to 2001 and included four sequences isolated from male MSM patients from Brazil, Argentina, Bulgaria and Spain.

Cluster C included two sequences, it dated back to June 2007. The virus of these two sequences was isolated in the same day and seems to be the same virus circulating in these two male patients.

Cluster D dated back 1989 and included three sequences isolated from one female and two males. Interestingly sequence 11288_RT was isolated from a male patient from Peru arrived in Italy in 2016 and the tMRCA of the cluster including this strain and one another sequences dated back to May 2016.

Cluster E included two sequences from male MSM patients, dating back to 2002.

Cluster F included two sequences from male and female patients, dating back to 2006.

Cluster G included three sequences from Albanian patients, dating back to 2008. Only of one strain was known the date of arrival in Italy, 9804_RT, the group that include this strain dated back to December 2011 and this patient arrived in Italy in 2011.

Cluster H included five sequences from male MSM patients, dating back to 2004.

Cluster I included two sequences from South-American patients, dating back to 2004, one of this two sequences belongs to a Colombian patient arrived in Italy in 2008.

Cluster J included three sequences from male patients, dating back to 2007; 92779_RT sequence was isolated from an Egyptian patient arrived in 2013 in Italy and the group that included this strain dated back in July 2011.

Conclusion

It is commonly accepted that human mobility is linked to the transmission of, and to the susceptibility to infectious diseases. Migrants, defined as individuals who move from their country of origin to another, account for 40% of newly-diagnosed cases of HIV infection in the European Union/European Economic Area. Populations at high risk for HIV infection include migrants from countries where the HIV prevalence is high and those participating in high risk behaviour. Travellers contribute to the spread of HIV-1 genetic diversity worldwide, and in the developing world migration of rural populations and civil war are additional contributing factors.

Here we describe the characteristics and evolution of HIV-1 in immigrant population in northern Italy observed between 2008–2018. Phylogenetic analysis revealed different clusters with different characteristic factors probably causing the transmission of the virus within these groups. People mostly arrive in Italy already infected in their country of origin, sometime the transmission having occurred within the familiar marriage and sometimes, as in MSM and sex workers, because of frequent unsafe sexual intercourses. Unfortunately, in some cluster we did not catch any information to deduce an epidemiological link between phylogenetic tree results and the possible source of infection.

Migration is a humanitarian problem requiring political decisions including specific measures to reduce the number of susceptible people and improving their living condition. In this regard, analysing the measles model in a migrant population, some authors suggested screening and vaccination campaigns as the key points to achieve an adequate immunisation coverage and disease control [Reference Ceccarelli56]. A clinical surveillance performed in 2016 proved the good state of health of a migrant population, but the associated microbiological surveillance performed by blood, rectal, pharyngeal and nasal swabs traced unusual microorganisms such as Pseudomonas (P.) putida, P. monteilii, P. fulva, P. mosselii; Aeromonas (A.) veronii, A. caviae, A. hydrophila, A. guillouiae; Acinetobacter (Ac.) Ac. lwoffii, Ac. johnsonii, Ac. tjernberg; Pantoea (Pa.) agglomerans and Pa. calida. Among these microorganisms, some strains were resistant to carbapenems or ESBL producers or methicillin resistant [Reference Ciccozzi57]. Therefore, there is an indication to apply to migrant population measures useful to reduce the chances of transmitting infectious agents (screening and prophylaxis) from migrants to the native population, from this to migrants and among migrants.

The data presented in this study underline the importance of molecular and epidemiological screening for HIV infection for all newly arrived migrants, being they asylum seekers or economic migrants, since the spreading of HIV-1 variants has implications for diagnostic, treatment and vaccine development. These results suggest a shift of HIV-1 prevention projects' focus and a continuous monitoring of HIV-1 molecular epidemiology among entry populations. Prevention efforts based on HIV molecular epidemiology may improve public health surveillance setting.

Conflict of interest

All the authors of the manuscript declare that they have no conflict of interest regarding this paper.

Footnotes

*

Equally contributed to this work, designed the study and wrote the manuscript.

References

1.Korber, B et al. (2001) Evolutionary and immunological implications of contemporary HIV-1 variation. British Medical Bulletin 58, 1942.CrossRefGoogle ScholarPubMed
2.Seillier-Moiseiwitsch, F et al. (1994) Genetic variability of the human immunodeficiency virus: statistical and biological issues. Annual Review of Genetics 28, 559596.CrossRefGoogle ScholarPubMed
3.Peeters, M and Sharp, PM (2000) Genetic diversity of HIV-1: the moving target. AIDS (London, England) 14, S129S140.Google ScholarPubMed
4.Van Heuverswyn, F et al. (2006) Human immunodeficiency viruses: SIV infection in wild gorillas. Nature 444, 164.CrossRefGoogle ScholarPubMed
5.Gurtler, L et al. (1994) A new subtype of human immunodeficiency virus type 1 (MVP-5180) from Cameroon. Journal of Virology 68, 15811585.CrossRefGoogle ScholarPubMed
6.Simon, F et al. (1998) Identification of a new human immunodeficiency virus type 1 distinct from group M and group O. Nature Medicine 4, 10321037.CrossRefGoogle Scholar
7.Ayouba, A et al. (2000) HIV1 group N among HIV-1-seropositive individuals in Cameroon. AIDS (London, England) 14, 26232625.CrossRefGoogle Scholar
8.D'Arc, M et al. (2015) Origin of the HIV-1 group O epidemic in western lowland gorillas. Proceedings of the National Academy of Sciences of the United States of America 112, E1343E1352.CrossRefGoogle Scholar
9.Paraskevis, D et al. (2009) Tracing the HIV-1 subtype B mobility in Europe: a phylogeographic approach. Retrovirology 6, 49.CrossRefGoogle ScholarPubMed
10.Kantor, R et al. (2005) Impact of HIV-1 subtype and antiretroviral therapy on protease and reverse transcriptase genotype: results of a global collaboration. PLOS Medicine 2, e112.CrossRefGoogle ScholarPubMed
11.Vermund, SH and Leigh-Brown, AJ (2012) The HIV epidemic: high-income countries. Cold Spring Harbor Perspectives in Medicine 2, a007195.CrossRefGoogle ScholarPubMed
12.Esbjörnsson, J et al. (2016) HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries. Virus Evolution 2, vew010.CrossRefGoogle ScholarPubMed
13.Pyne, MT et al. (2013) Large-scale analysis of the prevalence and geographic distribution of HIV-1 non-B variants in the United States. Journal of Clinical Microbiology 51, 26622669.CrossRefGoogle ScholarPubMed
14.Chaix, ML et al. (2013) Increasing HIV-1 non-B subtype primary infections in patients in France and effect of HIV subtypes on virological and immunological responses to combined antiretroviral therapy. Clinical Infectious Diseases 56, 880887.CrossRefGoogle ScholarPubMed
15.Lai, A et al. (2019) Advances on human immunodeficiency virus type 2 infection. EC Microbiology 15, 964973.Google Scholar
16.Balotta, C et al. (2001) Increasing prevalence of non-clade B HIV-1 strains in heterosexual men and women, as monitored by analysis of reverse transcriptase and protease sequences. Journal of Acquired Immune Deficiency Syndromes 27, 499505.CrossRefGoogle ScholarPubMed
17.Ciccozzi, M et al. (2012) HIV-1 non-B subtypes in Italy: a growing trend. New Microbiologica 5, 377386.Google Scholar
18.Baldanti, F et al. (2008) Changes in circulation of B and non-B HIV strains: spotlight on a reference centre for infectious diseases in Northern Italy. Journal of Medical Virology 80, 947952.CrossRefGoogle ScholarPubMed
19.Buonaguro, L et al. (2008) Molecular and phylogenetic analysis of HIV-1 variants circulating in Italy. Infectious Agents and Cancer 3, 13.CrossRefGoogle ScholarPubMed
20.Balotta, C et al. (1997) Comparable biological and molecular determinants in HIV type 1-infected longterm nonprogressors and recently infected individuals. AIDS Research and Human Retroviruses 13, 337341.CrossRefGoogle Scholar
21.Buonaguro, L et al. (2004) Screening of HIV-1 isolates by reverse heteroduplex mobility assay and identification of non-B subtypes in Italy. Journal of Acquired Immune Deficiency Syndromes 37, 12951306.CrossRefGoogle ScholarPubMed
22.Ciccozzi, M et al. (2007) An outbreak of HIV-1 subtype G among Italian injecting drug users. AIDS (London, England) 21, 12131215.CrossRefGoogle Scholar
23.Lai, A et al. (2010) Changing patterns in HIV-1 non-B clade prevalence and diversity in Italy over three decades. HIV Medicine 11, 593602.CrossRefGoogle ScholarPubMed
24.Ciccozzi, M et al. (2019) The phylogenetic approach for viral infectious disease evolution and epidemiology: an updating review. Journal of Medical Virology 91, 17071724.CrossRefGoogle ScholarPubMed
25.Hemelaar, J et al. (2011) Global trends in molecular epidemiology of HIV-1 during 2000-2007. AIDS (London, England) 25, 679689.CrossRefGoogle ScholarPubMed
26.Abecasis, AB et al. (2013) HIV-1 subtype distribution and its demographic determinants in newly diagnosed patients in Europe suggest highly compartmentalized epidemics. Retrovirology 10, 7.CrossRefGoogle ScholarPubMed
27.Hemelaar, J et al. (2011) Global trends in molecular epidemiology of HIV-1 during 2000–2007. AIDS (London, England) 25, 679689.CrossRefGoogle ScholarPubMed
28.Ciccozzi, M et al. (2013) HIV type 1 origin and transmission dynamics among different risk groups in sardinia: molecular epidemiology within the close boundaries of an Italian island. AIDS Research and Human Retroviruses 9, 404410.CrossRefGoogle Scholar
29.Presti A, L et al. (2013) Migration patterns of HIV-1 subtype B virus in Northern Italy. New Microbiologica 36, 7579.Google ScholarPubMed
30.Mabey, D and Mayaud, P (1997) Sexually transmitted diseases in mobile populations. Genitourinary Medicine 73, 1822.Google ScholarPubMed
31.Tiedje, K et al. (2014) A focus group study of healthy eating knowledge, practices, and barriers among adult and adolescent immigrants and refugees in the United States. International Journal of Behavioral Nutrition and Physical Activity 16, 63.CrossRefGoogle Scholar
32.European Center for Disease Prevention and Control (ECDC) (2014) ECDC TECHNICAL REPORT Assessing the burden of key infectious diseases affecting migrant populations. Available at https://ecdc.europa.eu/sites/portal/files/media/en/publications/Publications/assessing-burden-disease-migrant-populations.pdf.Google Scholar
33.Istituto Superiore Di Sanità (ISS)–Centro Operativo AIDS (2014) CoA Aggiornamento delle nuove diagnosi di infezione da HIV e dei casi di AIDS in Italia al 31 Dicembre 2013. Not Ist Super Sanità 27(Suppl. 1), 148.Google Scholar
34.European Centre for Disease Prevention and Control (ECDC) (2013) Migrant health: sexual transmission of HIV within migrant groups in the EU/EEA and implications for effective interventions. Available at https://ecdc.europa.eu/sites/portal/files/media/en/publications/Publications/Migrant-health-sexual-transmission.pdf.Google Scholar
35.Del Amo, J et al. (2011) The epidemiology of HIV and AIDS reports in migrants in the 27 European Union countries, Norway and Iceland: 1999–2006. European Journal of Public Health 21, 620626.CrossRefGoogle ScholarPubMed
36.Boyd, AE et al. (2005) Ethnic differences in stage of presentation of adults newly diagnosed with HIV-1 infection in south London. HIV Medicine 6, 5965.CrossRefGoogle ScholarPubMed
37.Sulis, G et al. (2014) The HIV/Migrants Study Group Clinical and epidemiological features of HIV/AIDS infection among migrants at first access to healthcare services as compared to Italian patients in Italy: a retrospective multicentre study, 2000–2010. Infection 42, 859867.CrossRefGoogle Scholar
38.Fakoya, I et al. (2008) Barriers to HIV testing for migrant black Africans in Western Europe. HIV Medicine 9(Suppl. 2), 2325.CrossRefGoogle ScholarPubMed
39.Alvarez-del Arco, D et al. (2013) HIV testing and counselling for migrant populations living in high-income countries: a systematic review. European Journal of Public Health 6, 10391045.CrossRefGoogle Scholar
40.Sagnelli, C et al. (2014) Molecular epidemiology of hepatitis B virus genotypes circulating in acute hepatitis B patients in the Campania region. Journal of Medical Virology 86, 16831693.CrossRefGoogle ScholarPubMed
41.Katoh, K et al. (2005) MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Research 33, 511518.CrossRefGoogle ScholarPubMed
42.Strimmer, K and von Haeseler, A (1997) Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proceedings of the National Academy of Sciences of the United States of America 94, 68156819.CrossRefGoogle ScholarPubMed
43.Salemi, M et al. (2008) High-resolution molecular epidemiology and evolutionary history of HIV-1 subtypes in Albania. PLoS One 3, e1390.CrossRefGoogle ScholarPubMed
44.Xia, X (2013) DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution. Molecular Biology and Evolution 30, 17201728.CrossRefGoogle ScholarPubMed
45.Kumar, S et al. (2016) MEGA7: molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Molecular Biology and Evolution 33, 18701874.CrossRefGoogle ScholarPubMed
46.Sagnelli, C et al. (2017) Molecular diversity in irregular or refugee immigrant patients with HBV-genotype-E infection living in the metropolitan area of Naples. Journal of Medical Virology 89, 10151024.CrossRefGoogle ScholarPubMed
47.Sagnelli, C et al. (2018) HBV molecular epidemiology and clinical condition of immigrants living in Italy C. Infection 46, 523531.CrossRefGoogle Scholar
48.Drummond, AJ and Rambaut, A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology 7, 214.CrossRefGoogle ScholarPubMed
49.Kass, R and Raftey, A (1995) Bayes factors. Journal of the American Statistical Association 90, 773795.CrossRefGoogle Scholar
50.Drummond, AJ et al. (2005) Bayesian coalescent inference of past population dynamics from molecular sequences. Molecular Biology and Evolution 22, 11851192.CrossRefGoogle ScholarPubMed
51.Drummond, AJ et al. (2002) Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics 61, 13071320.Google Scholar
52.Minin VN, B et al. (2008) Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics. Molecular Biology and Evolution 25, 14591471.CrossRefGoogle ScholarPubMed
53.Baele, G et al. (2013) Make the most of your samples: Bayes factor estimators for high-dimensional models of sequence evolution. BMC Bioinformatics 14, 85.CrossRefGoogle ScholarPubMed
54.Baele, G et al. (2013) Accurate model selection of relaxed molecular clocks in Bayesian phylogenetics. Molecular Biology and Evolution 30, 239243.CrossRefGoogle ScholarPubMed
55.Baele, G and Lemey, P (2008) Bayesian evolutionary model testing in the phylogenomics era: matching model complexity with computational efficiency. Bioinformatics (Oxford, England) 29, 19701979.CrossRefGoogle Scholar
56.Ceccarelli, G et al. (2018) Susceptibility to measles in migrant population: implication for policy makers. Journal of Travel Medicine 25. doi: 10.1093/jtm/tax080.CrossRefGoogle ScholarPubMed
57.Ciccozzi, M et al. (2018) Sentinel surveillance data from Eritrean migrants in Italy: the theory of “Healthy Migrants”. Travel Medicine and Infectious Disease 22, 5865.CrossRefGoogle Scholar
Figure 0

Table 1. Demographics and epidemiological characteristics of all the 129 HIV-1 positive immigrants enrolled, and according to HIV-1 serotype

Figure 1

Table 2. Biochemical and immunovirological characteristics of the 129 HIV-1 positive immigrants, and according to HIV-1 serotype.

Figure 2

Fig. 1. Phylogenetic relationship of the HIV-1 subtype-B sequences with the subtype-B reference sequences downloaded from the NCBI sequence database (https://www.ncbi.nlm.nih.gov/). The asterisk along a branch represents significant statistical support for the clade subtending that branch (bootstrap support 70%). The collapsed cluster in the figure were what not statistically significant.

Figure 3

Fig. 2. MCC tree of the second dataset. Branches are scaled in time. Significant posterior probability support (pp ≥ 0.9) as indicated by an asterisk. Clusters are highlighted.