Next Article in Journal
Deletion of CD2-Like (CD2v) and C-Type Lectin-Like (EP153R) Genes from African Swine Fever Virus Georgia-∆9GL Abrogates Its Effectiveness as an Experimental Vaccine
Next Article in Special Issue
Increasing Prevalence of HIV-1 Transmitted Drug Resistance in Portugal: Implications for First Line Treatment Recommendations
Previous Article in Journal
Synthesis, Structure, and Function of Human Adenovirus Small Non-Coding RNAs
Previous Article in Special Issue
What Should Health Departments Do with HIV Sequence Data?
 
 
viruses-logo
Article Menu

Article Menu

Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Nationwide Study about the Dispersal Patterns of the Predominant HIV-1 Subtypes A1 and B in Greece: Inference of the Molecular Transmission Clusters

by Evangelia Georgia Kostaki 1, Maria Gova 1, Georgios Adamis 2, Georgios Xylomenos 2, Maria Chini 3, Nikos Mangafas 3, Marios Lazanas 3, Simeon Metallidis 4, Olga Tsachouridou 4, Vasileios Papastamopoulos 5, Dimitrios Chatzidimitriou 6, Eleni Kakalou 5, Anastasia Antoniadou 7, Antonios Papadopoulos 7, Mina Psichogiou 8, Dimitrios Basoulis 8, Dimitrios Pilalas 6, Ifigeneia Papageorgiou 1, Dimitra Paraskeva 9, Georgios Chrysos 9, Vasileios Paparizos 10, Sofia Kourkounti 10, Helen Sambatakou 11, Vasileios Bolanos 11, Nikolaos V. Sipsas 12, Malvina Lada 13, Emmanouil Barbounakis 14, Evrikleia Kantzilaki 14, Periklis Panagopoulos 15, Vasilis Petrakis 15, Stelios Drimis 9, Charalambos Gogos 16, Angelos Hatzakis 1, Apostolos Beloukas 17,18,*, Lemonia Skoura 6,† and Dimitrios Paraskevis 1,*,†add Show full author list remove Hide full author list
1
Department of Hygiene, Epidemiology and Medical Statistics, Medical School, National and Kapodistrian University of Athens, 11527 Athens, Greece
2
1st Department of Internal Medicine, G. Gennimatas General Hospital, 11527 Athens, Greece
3
3rd Department of Internal Medicine-Infectious Diseases Unit, “Korgialeneio-Benakeio” Red Cross General Hospital, 11526 Athens, Greece
4
1st Department of Internal Medicine, AHEPA University Hospital, Medical School, Aristotle University of Thessaloniki, 54636 Thessaloniki, Greece
5
5th Department of Internal Medicine and Infectious Diseases, Evaggelismos General Hospital, 10676 Athens, Greece
6
National AIDS Reference Centre of Northern Greece, Department of Microbiology, Medical School, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
7
4th Department of Medicine, Attikon General Hospital, Medical School, National and Kapodistrian University of Athens, 12462 Athens, Greece
8
1st Department of Medicine, Laikon General Hospital, Medical School, National and Kapodistrian University of Athens, 11527 Athens, Greece
9
Department of Internal Medicine, Tzaneio General Hospital, 18536 Piraeus, Greece
10
HIV/AIDS Unit, A. Syngros Hospital of Dermatology and Venereology, 16121 Athens, Greece
11
HIV Unit, 2nd Department of Internal Medicine, Hippokration General Hospital, Medical School, National and Kapodistrian University of Athens, 11527 Athens, Greece
12
Department of Pathophysiology, Laikon General Hospital, Medical School, National and Kapodistrian University of Athens, 11527 Athens, Greece
13
2nd Department of Internal Medicine, Sismanogleion General Hospital, 15126 Marousi, Greece
14
Department of Internal Medicine, University Hospital of Heraklion “PAGNI”, Medical School, University of Crete, 71110 Heraklion, Greece
15
Department of Internal Medicine, University General Hospital, Democritus University of Thrace, 68100 Alexandroupolis, Greece
16
Department of Internal Medicine and Infectious Diseases, University Hospital of Patras, 26504 Rio, Greece
17
Department of Clinical Infection, Microbiology and Immunology, Institute of Infection and Global Health, University of Liverpool, Liverpool L697BE, UK
18
Department of Biomedical Sciences, School of Health Sciences, University of West Attica, 12243 Athens, Greece
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 20 September 2020 / Revised: 15 October 2020 / Accepted: 16 October 2020 / Published: 19 October 2020
(This article belongs to the Special Issue HIV Molecular Epidemiology for Prevention 2020)

Abstract

:
Our aim was to investigate the dispersal patterns and parameters associated with local molecular transmission clusters (MTCs) of subtypes A1 and B in Greece (predominant HIV-1 subtypes). The analysis focused on 1751 (28.4%) and 2575 (41.8%) sequences of subtype A1 and B, respectively. Identification of MTCs was based on phylogenetic analysis. The analyses identified 38 MTCs including 2–1518 subtype A1 sequences and 168 MTCs in the range of 2–218 subtype B sequences. The proportion of sequences within MTCs was 93.8% (1642/1751) and 77.0% (1982/2575) for subtype A1 and B, respectively. Transmissions within MTCs for subtype A1 were associated with risk group (Men having Sex with Men vs. heterosexuals, OR = 5.34, p < 0.001) and Greek origin (Greek vs. non-Greek origin, OR = 6.05, p < 0.001) and for subtype B, they were associated with Greek origin (Greek vs. non-Greek origin, OR = 1.57, p = 0.019), younger age (OR = 0.96, p < 0.001), and more recent sampling (time period: 2011–2015 vs. 1999–2005, OR = 3.83, p < 0.001). Our findings about the patterns of across and within country dispersal as well as the parameters associated with transmission within MTCs provide a framework for the application of the study of molecular clusters for HIV prevention.

1. Introduction

HIV remains a major health challenge to date. Although HIV transmissions have declined over the last few years, incidence remains high (http://www.UNAIDS.org). According to the UNAIDS/WHO, 1.7 million people became newly infected with HIV in 2019 (http://www.UNAIDS.org). HIV-1 genetic diversity is a hallmark of the virus and since the beginning of the epidemic, group M of HIV-1 has been genetically classified into pure subtypes (A–D, F–H, J–L), sub-subtypes, recombinant forms, and unclassified viruses [1,2]. The recombinant forms of the virus are further divided into circulating recombinant forms (CRFs) and unique recombinant forms (URFs) [1]. The prevalence of HIV-1 subtypes and recombinants varies greatly with the subtype C being predominant, followed by subtypes B, A, CRF02_AG, and CRF01_AE [3]. In the absence of an effective vaccine and cure, prevention of HIV-1 transmission provides the only strategy to mitigate its impact and reduce incidence (http://www.UNAIDS.org).
Phylogenetic analyses of HIV-1 have several applications in public health, including investigation of transmission dynamics and dispersal and outbreak investigation [4,5,6,7,8,9,10,11,12,13,14,15,16]. Recently, the Centers for Disease Control and Prevention (CDC) has discussed the benefits of HIV cluster analysis for public health action [17]. Investigation on HIV molecular clusters can provide important data about the dispersal patterns of the epidemic, as well as the characteristics of people living with HIV (PLHIV) infected within clusters [4,5,6,7,8,9,11,12,13,14,15,16,18]. Moreover, by investigating the molecular clusters, we can identify those at higher risk for HIV infection (e.g., the social links of recently infected PLHIV through contact tracing programs) and intervene with targeted strategies to prevent HIV [17].
HIV cluster analysis has been performed in many instances and there is increasing evidence that this knowledge can be useful for public health purposes. The benefits of cluster analysis include the ability to identify clustered infections that occurred recently and, therefore, focus on populations at higher risk of getting or transmitting HIV [17]. This knowledge can be extracted from data generated for HIV drug resistance testing, which is being performed widely.
Our previous studies have shown that subtypes A1 and B are the most prevalent, with the former following an increasing trend after 2000 in Greece [19]. The outbreak among people who inject drugs (PWID) in the Athens metropolitan area that ignited in 2011 was a different situation, where most viruses were classified as CRF14_BG followed by CRF35_AD and subtypes A1 and B [15,20]. The local clusters were of considerable importance for the high prevalence of non-nucleoside reverse transcriptase inhibitors (NNRTI) resistance mutations in Greece, where the majority of the most prevalent resistance mutations (e.g., E138A/G, K103N, Y181C) were due to onward transmissions within molecular transmission clusters (MTCs) [21].
Our aim was to provide a framework of HIV cluster analysis for public health action in Greece, by means of molecular epidemiology. Specifically, we aimed to identify the MTCs of subtypes A1 and B, which are the most prevalent clades in Greece, and to investigate the parameters associated with infections within MTCs. We also aimed to infer the within country patterns of HIV-1 dispersal.

2. Materials and Methods

2.1. Study Population

The study dataset included 6166 HIV-1 sequences available in the protease and partial reverse transcriptase regions (pol gene) obtained from PLHIV during 1999–2015 in Greece. This number accounts for 57.2% of all people diagnosed with HIV-1 during the same period in Greece (n = 10,787). Surveillance data were available from the National Public Health Organization (formerly, Hellenic Centre for Disease Control and Prevention) [22]. Specifically, our study sample consisted of 4790 (77.7%) sequences generated in Athens (sampled from Central/Southern Greece) and 1376 (22.3%) sequences generated in Thessaloniki (1298 (94.3%) sampled from Northern Greece and 78 (5.7%) from Crete). In the case of multiple available sequences per patient (sequences collected at different time points during the patient’s follow-up), only the earliest available sequence was included in the analysis. A unique serial number was used to specify each patient, while duplicate patients were detected and excluded from the study population. The current project has been approved by the Ethics Committee of the Medical School, National and Kapodistrian University of Athens.

2.2. Subtyping Analysis

HIV-1 subtyping was carried out by using the online automated subtyping tools COMET (COntext-based Modeling for Expeditious Typing) (http://comet.retrovirology.lu) and REGA (http://dbpartners.stanford.edu/RegaSubtyping/html/indexhbv.html). In addition, the study sequences (n = 6166) were analyzed phylogenetically along with a set of 216 reference sequences representative of all pure subtypes, sub-subtypes, and the majority of CRFs, collected from the HIV-1 sequence database (http://www.hiv.lanl.gov/). Sequence alignment and editing were performed on MEGA X [23]. Editing was conducted manually and according to the coding reading frame. Phylogenetic analysis was performed by the approximate maximum likelihood (ML) method, as implemented in FastTree v2.1 [24], using Generalized Time Reversible (GTR) as a substitution model with CAT approximation. Tree visualization and annotation were performed on the programs Dendroscope v3.5.7 (http://dendroscope.org) and FigTree v1.4 (http://tree.bio.ed.ac.uk/software/figtree/). Furthermore, all the study sequences which remained unclassified from the previous procedures were tested for the presence of recombination by using the bootscanning approach as implemented in SimPlot v3.5.1 [25]. HIV-1 subtyping confirmation by phylogenetic and recombination analysis has focused on the characterization of subtypes A1 and B, which are the most prevalent clades in Greece, as described previously [19,21].

2.3. Phylogenetic Analysis

The dispersal patterns of HIV-1 subtypes A1 and B in Greece were investigated by means of phylogenetic analysis. Analysis was performed separately for each subtype and was repeated in five replicates using a different set of randomly selected references, in order to minimize the possibility of overrepresentation of countries or geographic regions with a higher number of available sequences due to more frequent HIV drug resistance testing. Specifically, sequences from our study population (A1: 1751; B: 2575) were analyzed along with a randomly selected globally sampled dataset of sequences (A1: 1500; B: 2000), used as references. Reference sequences were available in the HIV-1 sequence database and their subtype was checked and confirmed by COMET. The phylogenetic trees were inferred by the ML method under the GTR + G4 substitution model as implemented in RAxML v8.2.10 [26].
Thereafter, phylogenetic clusters consisting of at least two sequences sampled from Greece at a proportion greater than 70% compared to the total number of sequences within the cluster (geographic criterion) were defined as tentative MTCs [27,28]. Tentative MTCs that were present and fulfilled the geographic criterion in all five rounds of phylogenetic analysis were defined as MTCs. Robustness of all the large MTCs (i.e., MTCs including more than 10 sequences) was confirmed by ML transfer bootstrap observation (TBE) and Bayesian phylogenetic analysis [29]. Specifically, phylogenetic analysis using either ML or Bayesian methods were performed for each individual MTC using as references the most closely related sequences identified after a BLAST search using the HIV BLAST tool (http://www.hiv.lanl.gov). Only unique sequences identified using 10 BLAST matches to each query were included in the analysis. ML phylogenetic analysis was conducted by the RAxML v8.2.10 program, as specified in the above paragraph. Bayesian analysis was performed by using MrBayes [30] with GTR + G and run for 10 × 106 generations (burnin: 10%). The convergence of all parameters of the Markov chain Monte Carlo (MCMC) was assessed by Tracer [31], if the effective sample size (ESS) was higher than 200. MTCs receiving bootstrap support > 75% or posterior probability support > 0.85 were considered as significant.

2.4. Phylogeographic Analysis

In order to infer the within country patterns of HIV-1 dispersal, we performed phylogeographic analysis on all the large MTCs of subtypes A1 and B. Specifically, phylogeographic analysis was performed separately on the MTCs of each subtype over a high number of bootstrap trees (n = 408), reconstructed by ML phylogenetic analysis that was conducted on the RAxML v8.2.10 program. Initially, we tested if the HIV transmissions (viral migration events) occurred more frequently among PLHIV from the same geographic area (Central/Southern Greece, Northern Greece, Crete) in a significant manner. The hypothesis testing was performed by character reconstruction using the criterion of parsimony on Mesquite v3.61 [32]. Afterwards, we quantified the viral migration events between the different geographic areas by using PAUP*4.0 [33], and we assessed if they were higher than those expected by chance (panmixis). Further details about the way of estimating viral migration events has been described in detail elsewhere [34,35,36].

2.5. Statistical Analysis

Mean and standard deviation (SD)/median and interquartile range (IQR) for continuous variables and absolute and relative frequencies for categorical variables were used to summarize the demographic data of the study. The non-parametric Mann–Whitney U test was used for the comparisons of the relevant distributions, by defining the level of significance at 5% and by adjusting the level of significance according to the Bonferroni correction for multiple comparisons. The demographic parameters associated with the local transmission of subtypes A1 and B in Greece were estimated by multivariable logistic regression models that were fit to a subset of the original data consisting of complete observations (A: 1635; B: 2453). Presence in MTCs was the binary outcome variable, while age, gender, transmission risk group, origin, sampling period, and area of sampling were chosen as possible explanatory variables in both models. Statistical analyses were performed on STATA 13.1 (StataCorp LP).

2.6. Sequences Availability

The whole dataset corresponds to a dense sampling of PLHIV in Greece. Therefore, to avoid the risk of PLHIV identification, sequences will be available upon reasonable request.

3. Results

3.1. Subtyping Results

Subtyping analysis performed by COMET and REGA on our study sequences (n = 6166) revealed that the vast majority of the sequences were classified into subtypes B (2575; 41.8%) and A1 (1751; 28.4%). CRF14_BG (469; 7.6%), subtype C (175; 2.8%), CRF35_AD (160; 2.6%), CRF56_cpx (148; 2.4%), and CRF02_AG (144; 2.3%) were also detected at low prevalence. In addition, 153 (2.5%) sequences were found to belong to other subtypes and 173 (2.8%) sequences to other CRFs, while the rest of them (418; 6.8%) were classified as recombinant or unclassified forms of the virus.
Our analysis has focused on the detection of the two most prevalent clades (A1 and B) in Greece. The demographic characteristics of the PLHIV infected with subtypes A1 and B are presented in detail in Table 1. In both subtypes, most of the PLHIV were males (B: 86.1%; A1: 81.9%) of Greek origin (B: 82.2%; A1: 75.4%) with similar age (in years, median values: B: 36; A1: 37). In addition, the predominant mode of infection in this population was through sex between men (Men having Sex with Men—MSM) (B: 66.7%; A: 59.1%). A higher proportion of heterosexuals were found among those infected with subtype A1 (22.1%) versus subtype B (17.0%). Furthermore, a significant decreasing trend was observed for the prevalence of subtype B over the study period (1999–2015) (p = 0.001), using the sampling years of the sequences as a proxy of the diagnosis years of the treatment-naive PLHIV infected with subtype B. The decreasing trend of subtype B was also found to be significant over the time period between 2005 and 2015, when HIV drug resistance testing was performed routinely on treatment-naïve PLHIV and the sampling coverage was higher (p = 0.012). No significant trend was observed over time for subtype A1 (1999–2015: p = 0.410; 2005–2015: p = 0.139).

3.2. Subtypes A1 and B Dispersal Patterns and Regional Transmission

Phylogenetic analysis using large sets of randomly selected reference sequences was performed in five replicates for each one of the most prevalent subtypes (A1 and B) in order to characterize the spread of the HIV-1 epidemic in Greece. Further analysis confirmed the phylogenetic credibility of MTCs with more than 10 sequences (large MTCs). A new version of phylogenetic bootstrap (TBE) was selected since with large phylogenies, the classical Felsenstein’s bootstrap proportion (FBP) tends to give low support, especially at deep branches [29]. Estimation of MTCs revealed the levels of regional transmission per subtype by dividing the number of HIV-1 sequences from Greece found within MTCs with the total number of available sequences from Greece. Levels of regional transmission corresponds to PLHIV who have been infected locally.
Analyses revealed that the subtype A1 and B sequences from Greece clustered at different points in the ML trees (Figure 1a and Figure 2a). In addition, it was found that a considerable number of clusters fulfilled the criteria to be considered as MTCs for both subtypes. Specifically, 93.8% (1642 of 1751) of subtype A1 sequences formed 38 MTCs, with a range of 2–1518 sequences (Figure 1b and Figure S1a). For subtype A1, a large MTC was identified (Figure 1b), while for subtype B, several smaller MTCs were detected across the tree (Figure 2b and Figure S1b). Even though the median number of subtype A1 sequences per MTC was 2 (IQR: 2–4), 86.7% (1518) of the subtype A1 sequences fell within a large single MTC (Figure 1b). This finding indicates that the subtype A1 epidemic in Greece is highly monophyletic and corresponds to transmissions that occurred locally. This large single MTC received high TBE bootstrap support (99%), and according to the sensitivity analysis that was performed previously [19], it is a single outbreak. The evolutionary distances in this MTC were quite large since this strain has been introduced at the early stage of the HIV-1 epidemic in Greece. Specifically, the date of the most recent common ancestor of the subtype A1 in Greece was estimated to be 1977.9 (95% highest posterior density interval: 1973.7–1981.9) [19]. In addition, most of the PLHIV infected within the large single MTC of subtype A1 were males (n = 1305; 86.0%) of Greek origin (n = 1219; 80.3%) with a median age of 36 years (IQR: 29–46). The predominant mode of infection in this population was through sex between men (MSM) (n = 999; 65.8%), while most of the sequences from these PLHIV were sampled from Central/Southern Greece (n= 1066; 70.2%). Furthermore, phylogenetic analysis of subtype B sequences revealed the existence of 168 MTCs with a range of 2–218 sequences (Figure 2b and Figure S1b). The proportion of subtype B sequences within MTCs was 77.0% (1982 of 2575). The high levels of dispersal of subtype B (Figure 2a) and the existence of a large number of MTCs (Figure 2b) indicate that the introduction of this subtype into Greece has occurred from multiple sources.
The number of subtype A1 and B introductions into Greece was calculated separately for each subtype by adding the number of sequences found outside the MTCs (unclustered sequences—each corresponding to an independent introduction) and the number of MTCs (e.g., for two MTCs including 5 and 10 sequences, we counted only two introductions—as the number of founder strains within each MTC). The number of introductions was found equal to 147 (38 MTCs + 109 unclustered sequences) and 761 (168 MTCs + 593 unclustered sequences) for subtypes A1 and B, respectively. The proportion of founder strains transmitted locally (MTCs) was similar for subtype A1 (25.8%; 38/147) and B (22.1%; 168/759).

3.3. Characteristics of Transmissions Within and Outside MTCs

Our analysis revealed that PLHIV whose sequences found within the MTCs (infected locally) were mostly male (A1: 1373; 83.6%/B: 1729; 87.2%), MSM (A1: 1022; 62.2%/B: 1373; 69.3%), and of Greek origin (A1: 1283; 78.1%/B: 1635; 82.5%) (percentages correspond to PLHIV infected within MTCs) (Table 2). For both subtypes, PLHIV with a more recent diagnosis were more frequently found within MTCs (A1: 2011–2015: 50.5%; 2006–2010: 28.0%; 1999–2005: 21.3%/B: 2011–2015: 41.5%; 2006–2010: 25.7%; 1999–2005: 32.6%) (Table 2).
PLHIV infected with subtype A1 outside MTCs (n = 109) were reported mostly as heterosexuals (48; 44.0%), PWID (26; 23.9%), of non-Greek origin (47; 43.1%), and with an available sequence with a more recent sampling time (2011–2015: 45.9%; 2006–2010: 40.4%; 1999–2005: 13.7%) (percentages were estimated for PLHIV infected outside MTCs) (Table 2). For subtype B, infections outside MTCs were more frequent for males (488; 82.3%), MSM (345; 58.2%), PLHIV of Greek origin (482; 81.3%), and PLHIV with an available sequence with an earlier sampling time (2011–2015: 18.9%; 2006–2010: 18.2%; 1999–2005: 62.9%).

3.4. Factors Associated with Local Transmission of Subtypes A1 and B

The results of the statistical analysis (multivariable logistic regression models) for the identification of demographic factors associated with local transmission of subtypes A1 and B in Greece are presented in detail in Table 3. Analysis revealed that transmissions within MTCs of subtype A1 were associated with the MSM risk group (MSM vs. heterosexuals, OR = 5.34, p < 0.001) and Greek origin (Greek vs. non-Greek origin, OR = 6.05, p < 0.001), having been adjusted for the rest of the incorporated variables (Table 3). For subtype B, local transmissions were associated with Greek origin (Greek vs. non-Greek origin, OR = 1.57 p = 0.019) and more recent sampling (sampling period: 2011–2015 vs. 1999–2005, OR = 3.83, p < 0.001) (Table 3). In addition, it was found that younger age was also associated with an increased probability of local transmission of subtype B (OR = 0.96, p < 0.001) (Table 3).

3.5. Patterns of Geographical Dispersal

To investigate the patterns of virus spread among different geographic areas in Greece, we performed phylogeographic analysis on the sequences found within MTCs that were generated in Athens and sampled from most areas in Greece (area 1), were generated in Thessaloniki and sampled from Northern Greece (area 2), and a small sample of sequences from Crete also generated in Thessaloniki (area 3). Phylogeographic analyses suggested that the total number of viral migration events corresponding to transmissions among the three different areas was 181 (median value; IQR: 178–184) and 136 (median value; IQR: 134–138) for subtypes A1 and B, respectively. These numbers were significantly lower than the total number of viral migrations expected by chance under the simulated scenario of panmixis (i.e., random mixing of sequences assuming that the level of regional transmission is negligible). Specifically, the total number of viral migrations under panmixis were 398 (median value; IQR: 394–402) and 259 (median value; IQR: 257–262) for subtypes A1 and B, respectively, and were significantly higher than the estimated values for the bootstrap trees for both subtypes (p < 0.001). These findings suggest that transmissions occur at the regional level; a conclusion that can also be drawn by the visual inspection of the trees, where the existence of several regional clusters is obvious (Figure 3).
Additionally, we quantified the viral migration events among the different areas and we assessed if they were higher than those expected by chance under panmixis. According to the migration matrix (Table 4), the number of transmissions that were exported from Central/Southern Greece to the other localities (Northern Greece: subtype A1—90.9; subtype B—67.6/Crete: subtype A1—11.6; subtype B—15.0) was lower than those expected by chance, and only transmissions exported from Northern Greece to Central/Southern Greece (subtype A1: 35.4; subtype B: 17.4) were significantly higher than panmixis. Similar results were found between Crete and Central/Southern Greece; however, the number of transmissions was very low (subtype A1: 0.1; subtype B: 1.2) (Table 4).

4. Discussion

In the current study, we investigated the patterns of HIV-1 transmission across and among different geographic areas, and we estimated the parameters associated with regional transmission for the two most prevalent subtypes in Greece. We found that a high proportion of subtypes A1 and B, which are the predominant clades, occur locally at 93.8% and 77.0%, respectively. These provide some of the largest proportions of regional transmissions compared to approximately 50% that was reported previously [18,37,38,39,40]. Possible explanations are the dense sample used in our analysis (our study population represents the 57.2% of the total population of PLHIV diagnosed in Greece during 1999–2015) and the method implemented for the identification of MTCs that was based on phylogenetic analysis using a novel bootstrap approach, which is more appropriate for large datasets [29]. Moreover, we did not implement a genetic threshold as a criterion for the identification of MTCs, since our aim was to identify not only PLHIV with close transmission links, but all individuals infected locally. For subtype B, where the dispersal pattern was not monophyletic as was found for subtype A1, the regional transmissions were estimated equal to 77.0%, which is much higher than the estimations of previous studies (reasons mentioned above). Notably, Greece provides one of the few countries in Europe where subtype A1 is found at high prevalence and is associated with the local population and not with the PWID risk group [19,41]. Similarly, it has been reported for Albania and Cyprus [39,42], while for most countries in Western and Central Europe, subtype C and recombinant forms are the most prevailing viruses [3]. The subtype A1 epidemic in Greece differs from the A6 epidemic in Russia and Eastern European countries, where the majority of transmission have been found among the PWID and the origin of infection is from a distinct route [19,43,44,45,46,47]. Furthermore, a recent study showed that Greece provides the main source of subtype A1 transmissions in Europe and remote locations, such as the USA and Australia [48].
The proportion of subtype A1 and B introductions leading to MTCs were similar, suggesting that approximately one-quarter of cases were associated with different levels of onward transmissions in MTCs. These findings indicate that one of the reasons why the non-B subtypes are less frequently associated with transmission chains is probably that they are less abundant in European countries than subtype B [49] (http://www.hiv.lanl.gov). Although the estimated proportion of strains leading to MTCs can be extrapolated as a proportion of onward transmission for the whole population infected within the MTCs, it provides a proxy of the proportion of cases providing as sources for HIV infection. Given that our study population is a subsample of the total number of PLHIV in Greece, the previous proportion should be considered as the minimum estimate.
Our analysis suggested that regional transmission was associated for both subtypes with Greek origin, while the MSM risk group was identified as an additional associated parameter only for subtype A1 and more recent diagnosis and younger age for subtype B, respectively. These findings are in accordance with findings from Europe and the USA, where the MSM risk group, male gender, younger age, and recent diagnosis were the main factors associated with HIV-1 infections within MTCs [18,37,38,50,51,52], suggesting that prevention efforts should focus on populations with these characteristics. The exact effect of the time of sampling to MTCs cannot be estimated, since for sequences sampled earlier, identifying their transmission links is less likely and therefore, the likelihood is lower to be placed in MTCs [18].
Regarding the dispersal patterns of the virus, our analyses focused on sequences found within MTCs and we found that regional dispersal dominated within the geographic areas of study. This finding was supported by the significantly lower number of migration events found among the different areas as compared to the migration events expected under the scenario of random mixing of the sequences. The quantification of viral migration events among the different localities suggested a limited number of transmissions from the largest area to the two others, whereas transmissions from Northern Greece to Central/Southern Greece were higher than expected by chance. We should note that the number of transmissions from Central/Southern Greece to the other areas was the largest, but according to statistical phylogeography, they were less than those expected under the random mixing scenario. The high numbers of transmissions from Central/Southern Greece were probably due to the larger sample size of this area compared to Northern Greece or Crete; however, these numbers were lower than what was expected by chance. These findings indicate that Northern Greece, although it is smaller in sample size than Central/Southern Greece, transmits the virus at high rates to the latter area. Therefore, besides the major local dispersal pattern observed within the three different areas, some cross-regional transmissions occur, with Northern Greece acting as a major source for this type of viral dispersal. A similar pattern has been observed for the dispersal of the most prevalent mutations associated with resistance to antiretroviral drugs (ARV), where the major pattern included regional dispersal with some cases of cross-regional transmissions [21,53]. Similar findings were reported elsewhere, highlighting the important effect of transmission clusters in the increasing prevalence of NNRTI resistance mutations among treatment-naïve populations [37,51,54,55,56].
In conclusion, our study provides a detailed picture of the characteristics of the HIV-1 epidemic regarding the dispersal patterns of the most prevalent subtypes in Greece. Since we are one step before the implementation of molecular surveillance as a method to identify PLHIV and their links to higher risk for HIV infection, our aim was to perform a comprehensive analysis of MTCs and their characteristics. We identified the MTCs using a novel approach for bootstrapping and we showed that HIV-1 regional transmission dominates in Greece and, also, it is associated with the MSM risk group (subtype A1) and Greek origin (both subtypes). Overall, our findings provide some novel insights that can be important for national policy against HIV. Moreover, we showed that approximately one-quarter of introductions leads to regional dispersal. This estimate might provide a proxy of the percentage of PLHIV contributed to onward transmissions.

Supplementary Materials

The following is available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/10/1183/s1, Figure S1: Distribution of number of sequences per molecular transmission cluster (MTC) for: (a) subtype A1 and (b) subtype B.

Author Contributions

Conceptualization, E.G.K., A.H., A.B., L.S. and D.P. (Dimitrios Paraskevis); methodology, E.G.K., M.G. and D.P. (Dimitrios Paraskevis); formal analysis, E.G.K., M.G., I.P., L.S. and D.P. (Dimitrios Paraskevis); investigation, G.A., G.X., M.C., N.M., M.L. (Marios Lazanas), S.M., O.T., V.P. (Vasileios Papastamopoulos), D.C., E.K. (Eleni Kakalou), A.A., A.P., M.P., D.B., D.P. (Dimitrios Pilalas), D.P. (Dimitra Paraskeva), G.C., V.P. (Vasileios Paparizos), S.K., H.S., V.B., N.V.S., M.L. (Malvina Lada), E.B., E.K. (Evrikleia Kantzilaki), P.P., V.P. (Vasilis Petrakis), S.D., C.G. and L.S.; data curation, E.G.K. and M.G.; writing—original draft preparation, E.G.K. and D.P. (Dimitrios Paraskevis); writing—review and editing, E.G.K., M.G., G.A., G.X., M.C., N.M., M.L. (Marios Lazanas), S.M., O.T., V.P. (Vasileios Papastamopoulos), D.C., E.K. (Eleni Kakalou), A.A., A.P., M.P., D.B., D.P. (Dimitrios Pilalas), I.P., D.P. (Dimitra Paraskeva), G.C., V.P. (Vasileios Paparizos), S.K., H.S., V.B., N.V.S., M.L. (Malvina Lada), E.B., E.K. (Evrikleia Kantzilaki), P.P., V.P. (Vasilis Petrakis), S.D., C.G., A.H., A.B., L.S. and D.P. (Dimitrios Paraskevis); visualization, E.G.K.; supervision, D.P. (Dimitrios Paraskevis); project administration, L.S. and D.P. (Dimitrios Paraskevis); funding acquisition, L.S. and D.P. (Dimitrios Paraskevis); All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Gilead Sciences Hellas, project reference number 4205/22-01-2018. Article Processing Charge (APC) was covered by the University of Liverpool Gold Open Access scheme.

Acknowledgments

A part of this work was presented at the HIV Drug Therapy Glasgow 2018 (HIV Glasgow 2018), Glasgow, United Kingdom, 28–31 October 2018 (Poster ID P037).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Foley, B.T.; Leitner, T.; Paraskevis, D.; Peeters, M. Primate immunodeficiency virus classification and nomenclature: Review. Infect. Genet. Evol. 2016, 46, 150–158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Yamaguchi, J.; Vallari, A.; McArthur, C.; Sthreshley, L.; Cloherty, G.A.; Berg, M.G.; Rodgers, M.A. Complete Genome Sequence of CG-0018a-01 Establishes HIV-1 Subtype L. J. Acquir Immune Defic. Syndr. 2020, 83, 319–322. [Google Scholar] [CrossRef] [PubMed]
  3. Hemelaar, J.; Elangovan, R.; Yun, J.; Dickson-Tetteh, L.; Fleminger, I.; Kirtley, S.; Williams, B.; Gouws-Williams, E.; Ghys, P.D.; Abimiku, A.l.G.; et al. Global and regional molecular epidemiology of HIV-1, 1990–2015: A systematic review, global survey, and trend analysis. Lancet Infect. Dis. 2019, 19, 143–155. [Google Scholar] [CrossRef]
  4. Abeler-Dorner, L.; Grabowski, M.K.; Rambaut, A.; Pillay, D.; Fraser, C.; PANGEA Consortium. PANGEA-HIV 2: Phylogenetics And Networks for Generalised Epidemics in Africa. Curr. Opin. HIV AIDS 2019, 14, 173–180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Van de Vijver, D.; Boucher, C.A.B. Insights on transmission of HIV from phylogenetic analysis to locally optimize HIV prevention strategies. Curr. Opin. HIV AIDS 2018, 13, 95–101. [Google Scholar] [CrossRef] [PubMed]
  6. Brenner, B.G.; Ibanescu, R.I.; Hardy, I.; Roger, M. Genotypic and Phylogenetic Insights on Prevention of the Spread of HIV-1 and Drug Resistance in “Real-World” Settings. Viruses 2017, 10. [Google Scholar] [CrossRef] [Green Version]
  7. 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. Retrovirology 2017, 14, 15. [Google Scholar] [CrossRef] [Green Version]
  8. Grabowski, M.K.; Lessler, J. Phylogenetic insights into age-disparate partnerships and HIV. Lancet HIV 2017, 4, e8–e9. [Google Scholar] [CrossRef]
  9. Bezemer, D.; Cori, A.; Ratmann, O.; van Sighem, A.; Hermanides, H.S.; Dutilh, B.E.; Gras, L.; Rodrigues Faria, N.; van den Hengel, R.; Duits, A.J.; et al. Dispersion of the HIV-1 Epidemic in Men Who Have Sex with Men in the Netherlands: A Combined Mathematical Model and Phylogenetic Analysis. PLoS Med. 2015, 12, e1001898. [Google Scholar] [CrossRef] [Green Version]
  10. Vasylyeva, T.I.; Zarebski, A.; Smyrnov, P.; Williams, L.D.; Korobchuk, A.; Liulchuk, M.; Zadorozhna, V.; Nikolopoulos, G.; Paraskevis, D.; Schneider, J.; et al. Phylodynamics Helps to Evaluate the Impact of an HIV Prevention Intervention. Viruses 2020, 12, 469. [Google Scholar] [CrossRef] [Green Version]
  11. Paraskevis, D.; Nikolopoulos, G.K.; Magiorkinis, G.; Hodges-Mameletzis, I.; Hatzakis, A. The application of HIV molecular epidemiology to public health. Infect. Genet. Evol. 2016, 46, 159–168. [Google Scholar] [CrossRef] [PubMed]
  12. Brenner, B.G.; Wainberg, M.A. Future of phylogeny in HIV prevention. J. Acquir. Immune Defic. Syndr. 2013, 63, S248–S254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sivay, M.V.; Grabowski, M.K.; Zhang, Y.; Palumbo, P.J.; Guo, X.; Piwowar-Manning, E.; Hamilton, E.L.; Viet Ha, T.; Antonyak, S.; Imran, D.; et al. Phylogenetic Analysis of Human Immunodeficiency Virus from People Who Inject Drugs in Indonesia, Ukraine, and Vietnam: HPTN 074. Clin. Infect. Dis. 2019. [Google Scholar] [CrossRef]
  14. Alpren, C.; Dawson, E.L.; John, B.; Cranston, K.; Panneer, N.; Fukuda, H.D.; Roosevelt, K.; Klevens, R.M.; Bryant, J.; Peters, P.J.; et al. Opioid Use Fueling HIV Transmission in an Urban Setting: An Outbreak of HIV Infection Among People Who Inject Drugs-Massachusetts, 2015–2018. Am. J. Public Health 2020, 110, 37–44. [Google Scholar] [CrossRef]
  15. Kostaki, E.G.; Frampton, D.; Paraskevis, D.; Pantavou, K.; Ferns, B.; Raffle, J.; Grant, P.; Kozlakidis, Z.; Hadjikou, A.; Pavlitina, E.; et al. Near Full-length Genomic Sequencing and Molecular Analysis of HIV-Infected Individuals in a Network-based Intervention (TRIP) in Athens, Greece: Evidence that Transmissions Occur More Frequently from those with High HIV-RNA. Curr. HIV Res. 2018, 16, 345–353. [Google Scholar] [CrossRef] [PubMed]
  16. Peters, P.J.; Pontones, P.; Hoover, K.W.; Patel, M.R.; Galang, R.R.; Shields, J.; Blosser, S.J.; Spiller, M.W.; Combs, B.; Switzer, W.M.; et al. HIV Infection Linked to Injection Use of Oxymorphone in Indiana, 2014–2015. N. Engl. J. Med. 2016, 375, 229–239. [Google Scholar] [CrossRef] [Green Version]
  17. CDC: Detecting and Responding to HIV Transmission Clusters. A Guide for Health Departments. Available online: https://www.cdc.gov/hiv/pdf/funding/announcements/ps18-1802/CDC-HIV-PS18-1802-AttachmentE-Detecting-Investigating-and-Responding-to-HIV-Transmission-Clusters.pdf (accessed on 19 September 2020).
  18. Paraskevis, D.; Beloukas, A.; Stasinos, K.; Pantazis, N.; de Mendoza, C.; Bannert, N.; Meyer, L.; Zangerle, R.; Gill, J.; Prins, M.; et al. HIV-1 molecular transmission clusters in nine European countries and Canada: Association with demographic and clinical factors. BMC Med. 2019, 17, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Paraskevis, D.; Magiorkinis, E.; Magiorkinis, G.; Sypsa, V.; Paparizos, V.; Lazanas, M.; Gargalianos, P.; Antoniadou, A.; Panos, G.; Chrysos, G.; et al. Increasing prevalence of HIV-1 subtype A in Greece: Estimating epidemic history and origin. J. Infect. Dis. 2007, 196, 1167–1176. [Google Scholar] [CrossRef] [Green Version]
  20. Paraskevis, D.; Paraschiv, S.; Sypsa, V.; Nikolopoulos, G.; Tsiara, C.; Magiorkinis, G.; Psichogiou, M.; Flampouris, A.; Mardarescu, M.; Niculescu, I.; et al. Enhanced HIV-1 surveillance using molecular epidemiology to study and monitor HIV-1 outbreaks among intravenous drug users (IDUs) in Athens and Bucharest. Infect. Genet. Evol. 2015, 35, 109–121. [Google Scholar] [CrossRef]
  21. Paraskevis, D.; Kostaki, E.; Magiorkinis, G.; Gargalianos, P.; Xylomenos, G.; Magiorkinis, E.; Lazanas, M.; Chini, M.; Nikolopoulos, G.; Skoutelis, A.; et al. Prevalence of drug resistance among HIV-1 treatment-naive patients in Greece during 2003–2015: Transmitted drug resistance is due to onward transmissions. Infect. Genet. Evol. 2017, 54, 183–191. [Google Scholar] [CrossRef]
  22. NPHO: HIV/AIDS Surveillance in Greece Diagnoses through 31/12/2019. Available online: https://eody.gov.gr/wp-content/uploads/2018/12/epidimiologiko-deltio-hiv-2019.pdf (accessed on 19 September 2020).
  23. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  24. Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree 2—Approximately maximum-likelihood trees for large alignments. PLoS ONE 2010, 5, e9490. [Google Scholar] [CrossRef] [PubMed]
  25. Lole, K.S.; Bollinger, R.C.; Paranjape, R.S.; Gadkari, D.; Kulkarni, S.S.; Novak, N.G.; Ingersoll, R.; Sheppard, H.W.; Ray, S.C. Full-Length Human Immunodeficiency Virus Type 1 Genomes from Subtype C-Infected Seroconverters in India, with Evidence of Intersubtype Recombination. J. Virol. 1999, 73, 152–160. [Google Scholar] [CrossRef] [Green Version]
  26. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef]
  27. Kostaki, E.G.; Flampouris, A.; Karamitros, T.; Chueca, N.; Alvarez, M.; Casas, P.; Alejos, B.; Hatzakis, A.; Garcia, F.; Paraskevis, D.; et al. Spatiotemporal Characteristics of the Largest HIV-1 CRF02_AG Outbreak in Spain: Evidence for Onward Transmissions. Front. Microbiol. 2019, 10, 370. [Google Scholar] [CrossRef] [PubMed]
  28. Von Wyl, V.; Kouyos, R.D.; Yerly, S.; Boni, J.; Shah, C.; Burgisser, P.; Klimkait, T.; Weber, R.; Hirschel, B.; Cavassini, M.; et al. The role of migration and domestic transmission in the spread of HIV-1 non-B subtypes in Switzerland. J. Infect. Dis. 2011, 204, 1095–1103. [Google Scholar] [CrossRef] [PubMed]
  29. Lemoine, F.; Domelevo Entfellner, J.B.; Wilkinson, E.; Correia, D.; Davila Felipe, M.; De Oliveira, T.; Gascuel, O. Renewing Felsenstein’s phylogenetic bootstrap in the era of big data. Nature 2018, 556, 452–456. [Google Scholar] [CrossRef]
  30. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Hohna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Maddison, W.P.; Maddison, D.R. Mesquite: A Modular System for Evolutionary Analysis. Version 3.61. Available online: http://www.mesquiteproject.org (accessed on 19 September 2020).
  33. Wilgenbusch, J.C.; Swofford, D. Inferring evolutionary trees with PAUP*. Curr. Protoc. Bioinform. 2003. [Google Scholar] [CrossRef]
  34. Angelis, K.; Albert, J.; Mamais, I.; Magiorkinis, G.; Hatzakis, A.; Hamouda, O.; Struck, D.; Vercauteren, J.; Wensing, A.M.; Alexiev, I.; et al. Global Dispersal Pattern of HIV Type 1 Subtype CRF01_AE: A Genetic Trace of Human Mobility Related to Heterosexual Sexual Activities Centralized in Southeast Asia. J. Infect. Dis. 2015, 211, 1735–1744. [Google Scholar] [CrossRef] [PubMed]
  35. Magiorkinis, G.; Angelis, K.; Mamais, I.; Katzourakis, A.; Hatzakis, A.; Albert, J.; Lawyer, G.; Hamouda, O.; Struck, D.; Vercauteren, J.; et al. The global spread of HIV-1 subtype B epidemic. Infect. Genet. Evol. 2016, 46, 169–179. [Google Scholar] [CrossRef] [PubMed]
  36. Paraskevis, D.; Kostaki, E.; Nikolopoulos, G.K.; Sypsa, V.; Psichogiou, M.; Del Amo, J.; Hodges-Mameletzis, I.; Paraskeva, D.; Skoutelis, A.; Malliori, M.; et al. Molecular Tracing of the Geographical Origin of Human Immunodeficiency Virus Type 1 Infection and Patterns of Epidemic Spread Among Migrants Who Inject Drugs in Athens. Clin. Infect. Dis. 2017, 65, 2078–2084. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. 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. [Google Scholar] [CrossRef] [PubMed]
  38. Verhofstede, C.; Mortier, V.; Dauwe, K.; Callens, S.; Deblonde, J.; Dessilly, G.; Delforge, M.L.; Fransen, K.; Sasse, A.; Stoffels, K.; et al. Exploring HIV-1 Transmission Dynamics by Combining Phylogenetic Analysis and Infection Timing. Viruses 2019, 11, 1096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Kostrikis, L.G.; Hezka, J.; Stylianou, D.C.; Kostaki, E.; Andreou, M.; Kousiappa, I.; Paraskevis, D.; Demetriades, I. HIV-1 transmission networks across Cyprus (2010–2012). PLoS ONE 2018, 13, e0195660. [Google Scholar] [CrossRef] [Green Version]
  40. Lewis, F.; Hughes, G.J.; Rambaut, A.; Pozniak, A.; Brown, A.J.L. Episodic Sexual Transmission of HIV Revealed by Molecular Phylodynamics. PLoS Med. 2008, 5, e50. [Google Scholar] [CrossRef] [Green Version]
  41. Beloukas, A.; Psarris, A.; Giannelou, P.; Kostaki, E.; Hatzakis, A.; Paraskevis, D. Molecular epidemiology of HIV-1 infection in Europe: An overview. Infect. Genet. Evol. 2016, 46, 180–189. [Google Scholar] [CrossRef] [Green Version]
  42. Ciccozzi, M.; Gori, C.; Boros, S.; Ruiz-Alvarez, M.J.; Harxhi, A.; Dervishi, M.; Qyra, S.; Schinaia, N.; D’Arrigo, R.; Ceccherini-Silberstein, F.; et al. Molecular Diversity of HIV in Albania. J. Infect. Dis. 2005, 192, 475–479. [Google Scholar] [CrossRef] [Green Version]
  43. Aibekova, L.; Foley, B.; Hortelano, G.; Raees, M.; Abdraimov, S.; Toichuev, R.; Ali, S. Molecular epidemiology of HIV-1 subtype A in former Soviet Union countries. PLoS ONE 2018, 13, e0191891. [Google Scholar] [CrossRef]
  44. Schlosser, M.; Kartashev, V.V.; Mikkola, V.H.; Shemshura, A.; Saukhat, S.; Kolpakov, D.; Suladze, A.; Tverdokhlebova, T.; Hutt, K.; Heger, E.; et al. HIV-1 Sub-Subtype A6: Settings for Normalised Identification and Molecular Epidemiology in the Southern Federal District, Russia. Viruses 2020, 12, 475. [Google Scholar] [CrossRef] [Green Version]
  45. Bobkov, A.; Cheingsong-Popov, R.; Selimova, L.; Ladnaya, N.; Kazennova, E.; Kravchenko, A.; Fedotov, E.; Saukhat, S.; Zverev, S.; Pokrovsky, V.; et al. An HIV Type 1 Epidemic among Injecting Drug Users in the Former Soviet Union Caused by a Homogeneous Subtype A Strain. AIDS Res. Hum. Retrovir. 1997, 13, 1195–1201. [Google Scholar] [CrossRef] [PubMed]
  46. Thomson, M.M.; de Parga, E.V.; Vinogradova, A.; Sierra, M.; Yakovlev, A.; Rakhmanova, A.; Delgado, E.; Casado, G.; Munoz, M.; Carmona, R.; et al. New insights into the origin of the HIV type 1 subtype A epidemic in former Soviet Union’s countries derived from sequence analyses of preepidemically transmitted viruses. AIDS Res. Hum. Retrovir. 2007, 23, 1599–1604. [Google Scholar] [CrossRef]
  47. Diez-Fuertes, F.; Cabello, M.; Thomson, M.M. Bayesian phylogeographic analyses clarify the origin of the HIV-1 subtype A variant circulating in former Soviet Union’s countries. Infect. Genet. Evol. 2015, 33, 197–205. [Google Scholar] [CrossRef] [PubMed]
  48. Araujo, P.M.M.; Carvalho, A.; Pingarilho, M.; Abecasis, A.B.; Osorio, N.S. Characterization of a large cluster of HIV-1 A1 infections detected in Portugal and connected to several Western European countries. Sci. Rep. 2019, 9, 7223. [Google Scholar] [CrossRef] [Green Version]
  49. Abecasis, A.B.; Wensing, A.M.J.; Paraskevis, D.; Vercauteren, J.; Theys, K.; Van de Vijver, D.A.M.C.; Albert, J.; Asjö, B.; Balotta, C.; Beshkov, D.; et al. HIV-1 subtype distribution and its demographic determinants in newly diagnosed patients in Europe suggest highly compartmentalized epidemics. Retrovirology 2013, 10. [Google Scholar] [CrossRef] [PubMed]
  50. 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] [Green Version]
  51. Drescher, S.M.; von Wyl, V.; Yang, W.L.; Boni, J.; Yerly, S.; Shah, C.; Aubert, V.; Klimkait, T.; Taffe, P.; Furrer, H.; et al. 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] [PubMed] [Green Version]
  52. Kostaki, E.G.; Soulie, C.; Visseaux, B.; Storto, A.; Charpentier, C.; Wirden, M.; Landman, R.; Katlama, C.; Calvez, V.; Descamps, D.; et al. Molecular Analysis Suggests Post-Migration HIV-1 Acquisition among Migrants in Paris. In Proceedings of the Virtual Conference on Retroviruses and Opportunistic Infections, Boston, MA, USA, 8–11 March 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 19 September 2020).
  53. Skoura, L.; Metallidis, S.; Buckton, A.J.; Mbisa, J.L.; Pilalas, D.; Papadimitriou, E.; Papoutsi, A.; Haidich, A.B.; Chrysanthidis, T.; Tsachouridou, O.; et al. Molecular and epidemiological characterization of HIV-1 infection networks involving transmitted drug resistance mutations in Northern Greece. J. Antimicrob. Chemother. 2011, 66, 2831–2837. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. 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]
  55. Rhee, S.Y.; Blanco, J.L.; Jordan, M.R.; Taylor, J.; Lemey, P.; Varghese, V.; Hamers, R.L.; Bertagnolio, S.; Rinke de Wit, T.F.; Aghokeng, A.F.; et al. 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, e1001810. [Google Scholar] [CrossRef] [Green Version]
  56. Kleyn, T.J.; Liedtke, M.D.; Harrison, D.L.; Lockhart, S.M.; Salvaggio, M.R.; Ripley, T.L.; Rathbun, R.C. Incidence of transmitted antiretroviral drug resistance in treatment-naive HIV-1-infected persons in a large South Central United States clinic. Ann. Pharm. 2014, 48, 470–475. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Unrooted phylogenetic trees estimated by RAxML v8.2.10 of HIV-1 subtype A1 sequences from Greece and a global reference dataset. (a) Sequences from Greece are marked in light green (unclustered sequences) and red (sequences found within local molecular transmission clusters-MTCs) in contrast with sequences from other countries marked in blue. (b) The MTCs are indicated as triangles.
Figure 1. Unrooted phylogenetic trees estimated by RAxML v8.2.10 of HIV-1 subtype A1 sequences from Greece and a global reference dataset. (a) Sequences from Greece are marked in light green (unclustered sequences) and red (sequences found within local molecular transmission clusters-MTCs) in contrast with sequences from other countries marked in blue. (b) The MTCs are indicated as triangles.
Viruses 12 01183 g001
Figure 2. Unrooted phylogenetic trees estimated by RAxML v8.2.10 of HIV-1 subtype B sequences from Greece and a global reference dataset. (a) Sequences from Greece are marked in light green (unclustered sequences) and red (sequences found within local molecular transmission clusters-MTCs) in contrast with sequences from other countries marked in blue. (b) The MTCs are indicated as triangles.
Figure 2. Unrooted phylogenetic trees estimated by RAxML v8.2.10 of HIV-1 subtype B sequences from Greece and a global reference dataset. (a) Sequences from Greece are marked in light green (unclustered sequences) and red (sequences found within local molecular transmission clusters-MTCs) in contrast with sequences from other countries marked in blue. (b) The MTCs are indicated as triangles.
Viruses 12 01183 g002
Figure 3. Unrooted phylogenetic tree estimated by RAxML v8.2.10 of HIV-1 subtype A1 sequences from Greece that fell within the large single molecular transmission cluster (MTC). Sequences from Central/Southern Greece are marked in light blue, sequences from Northern Greece in orange, and sequences from Crete in dark green. Clusters consisting of sequences from Northern Greece and Crete are indicated as triangles.
Figure 3. Unrooted phylogenetic tree estimated by RAxML v8.2.10 of HIV-1 subtype A1 sequences from Greece that fell within the large single molecular transmission cluster (MTC). Sequences from Central/Southern Greece are marked in light blue, sequences from Northern Greece in orange, and sequences from Crete in dark green. Clusters consisting of sequences from Northern Greece and Crete are indicated as triangles.
Viruses 12 01183 g003
Table 1. Demographic characteristics of a subset of the study population (i.e., PLHIV infected with subtypes A1 and B).
Table 1. Demographic characteristics of a subset of the study population (i.e., PLHIV infected with subtypes A1 and B).
CharacteristicSubtype A1Subtype B
Age, years, median (IQR)37 (30–46) 136 (30–44) 2
Gender, n (%)
Male1.435 (81.9)2217 (86.1)
Female232 (13.3)282 (10.9)
Unknown84 (4.8)76 (3.0)
Transmission risk group, n (%)
MSM1035 (59.1)1718 (66.7)
PWID150 (8.6)182 (7.1)
Heterosexuals387 (22.1)438 (17.0)
Other20 (1.1)88 (3.4)
Unknown159 (9.1)149 (5.8)
Origin, n (%)
Greek1321 (75.4)2117 (82.2)
Non-Greek197 (11.3)186 (7.2)
Unknown233 (13.3)272 (10.6)
Sampling period
1999–2005365 (20.8)1020 (39.6)
2006–2010503 (28.7)618 (24.0)
2011–2015880 (50.3)934 (36.3)
Unknown3 (0.2)3 (0.1)
Total, n, (%)1751 (100.0)2575 (100.0)
1n = 1642 (of 1751 PLHIV); 2 n = 2457 (of 2575 PLHIV). PLHIV—people living with HIV; IQR—interquartile range; MSM—men having sex with men; PWID—people who inject drugs.
Table 2. Demographic characteristics of PLHIV who were infected with subtypes A1 and B within and outside molecular transmission clusters (MTCs).
Table 2. Demographic characteristics of PLHIV who were infected with subtypes A1 and B within and outside molecular transmission clusters (MTCs).
Subtype A1Subtype B
CharacteristicWithin
MTCs
Outside MTCsWithin
MTCs
Outside MTCs
Age, years, median (IQR)37 (29–46) 137 (31–48) 235 (29–42) 340 (34–48) 4
Gender, n (%)
Male1373 (83.6)62 (56.9)1729 (87.2)488 (82.3)
Female193 (11.8)39 (35.8)196 (9.9)86 (14.5)
Unknown76 (4.6)8 (7.3)57 (2.9)19 (3.2)
Transmission risk group, n (%)
MSM1022 (62.2)13 (11.9)1373 (69.3)345 (58.2)
PWID124 (7.6)26 (23.9)148 (7.5)34 (5.7)
Heterosexuals339 (20.6)48 (44.0)321 (16.2)117 (19.7)
Other18 (1.1)2 (1.8)23 (1.1)65 (11.0)
Unknown139 (8.5)20 (18.4)117 (5.9)32 (5.4)
Origin, n (%)
Greek1283 (78.1)38 (34.9)1635 (82.5)482 (81.3)
Non-Greek150 (9.2)47 (43.1)139 (7.0)47 (7.9)
Unknown209 (12.7)24 (22.0)208 (10.5)64 (10.8)
Sampling period
[1999–2005]350 (21.3)15 (13.7)647 (32.6)373 (62.9)
[2006–2010]459 (28.0)44 (40.4)510 (25.7)108 (18.2)
[2011–2015]830 (50.5)50 (45.9)822 (41.5)112 (18.9)
Unknown3 (0.2)0 (0.0)3 (0.2)0 (0.0)
Total, n, (%)1642 (100.0)109 (100.0)1982 (100.0)593 (100.0)
1n = 1540 (of 1642 PLHIV); 2 n = 102 (of 109 PLHIV); 3 n = 1891 (of 1982 PLHIV); 4 n = 566 (of 593 PLHIV). PLHIV—people living with HIV; IQR—interquartile range; MSM—men having sex with men; PWID—people who inject drugs.
Table 3. Multivariable logistic regression estimates using the presence in molecular transmission clusters (MTCs) as the binary outcome variable.
Table 3. Multivariable logistic regression estimates using the presence in molecular transmission clusters (MTCs) as the binary outcome variable.
Subtype A1 1Subtype B 1
Explanatory VariableOdds Ratio95% CIp-ValueOdds Ratio95% CIp-Value
Age0.990.97–1.010.3410.960.95–0.97<0.001
Gender
Male1.520.87–2.670.1401.330.91–1.970.141
Female 2
Unknown2.800.64–12.30.1720.840.22–3.160.796
Transmission risk group
MSM5.342.52–11.31<0.0010.980.69–1.380.893
PWID0.420.22–0.800.0080.800.49–1.300.364
Heterosexuals 2
Other0.860.17–4.470.8610.100.05–0.18<0.001
Unknown0.810.37–1.800.6100.960.48–1.950.919
Origin
Greek6.053.62–10.10<0.0011.571.08–2.290.019
Non-Greek 2
Unknown1.270.62–2.600.5161.260.74–2.120.394
Sampling period
[1999–2005] 2
[2006–2010]0.510.25–1.010.0542.331.80–3.02<0.001
[2011–2015]0.750.38–1.480.4113.832.95–4.96<0.001
Unknown1--1--
Sampling area
Central/Southern Greece 2
Northern Greece0.910.54–1.550.7361.320.96–1.810.092
Crete1--1--
1 Each model was fit to a subset of the original data (A1: 1751; B: 2575) consisting of 1635 subtype A1 and 2453 subtype B complete observations; 2 Reference category. CI—confidence interval; MSM—men having sex with men; PWID—people who inject drugs.
Table 4. Migration matrix representing the mean number of viral migrations events (transmissions) among the different geographic areas in Greece across all the phylogenetic bootstrap trees (n = 408), separately for each HIV-1 subtype under study.
Table 4. Migration matrix representing the mean number of viral migrations events (transmissions) among the different geographic areas in Greece across all the phylogenetic bootstrap trees (n = 408), separately for each HIV-1 subtype under study.
Subtype A1
Importing to
Exporting fromCentral/Southern GreeceNorthern GreeceCrete
Central/Southern Greece 90.911.6
Northern Greece35.41 0.2
Crete0.110
Subtype B
Central/Southern Greece 67.615.0
Northern Greece17.41 1.41
Crete1.210.51
1 Statistically significant pathways (compared with the null hypothesis representing the simulated scenario of panmixis). The level of significance has been adjusted according to the Bonferroni correction for multiple comparisons.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kostaki, E.G.; Gova, M.; Adamis, G.; Xylomenos, G.; Chini, M.; Mangafas, N.; Lazanas, M.; Metallidis, S.; Tsachouridou, O.; Papastamopoulos, V.; et al. A Nationwide Study about the Dispersal Patterns of the Predominant HIV-1 Subtypes A1 and B in Greece: Inference of the Molecular Transmission Clusters. Viruses 2020, 12, 1183. https://0-doi-org.brum.beds.ac.uk/10.3390/v12101183

AMA Style

Kostaki EG, Gova M, Adamis G, Xylomenos G, Chini M, Mangafas N, Lazanas M, Metallidis S, Tsachouridou O, Papastamopoulos V, et al. A Nationwide Study about the Dispersal Patterns of the Predominant HIV-1 Subtypes A1 and B in Greece: Inference of the Molecular Transmission Clusters. Viruses. 2020; 12(10):1183. https://0-doi-org.brum.beds.ac.uk/10.3390/v12101183

Chicago/Turabian Style

Kostaki, Evangelia Georgia, Maria Gova, Georgios Adamis, Georgios Xylomenos, Maria Chini, Nikos Mangafas, Marios Lazanas, Simeon Metallidis, Olga Tsachouridou, Vasileios Papastamopoulos, and et al. 2020. "A Nationwide Study about the Dispersal Patterns of the Predominant HIV-1 Subtypes A1 and B in Greece: Inference of the Molecular Transmission Clusters" Viruses 12, no. 10: 1183. https://0-doi-org.brum.beds.ac.uk/10.3390/v12101183

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