1. Introduction
HIV remains a major health challenge to date. Although HIV transmissions have declined over the last few years, incidence remains high (
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) (
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 (
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] (
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
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.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figures and Tables
Enlarge this image.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.
Enlarge this image.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.
Enlarge this image.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.
Demographic characteristics of a subset of the study population (i.e., PLHIV infected with subtypes A1 and B).
| Characteristic | Subtype A1 | Subtype B |
|---|---|---|
| Age, years, median (IQR) | 37 (30–46) 1 | 36 (30–44) 2 |
| Gender, n (%) | ||
| Male | 1.435 (81.9) | 2217 (86.1) |
| Female | 232 (13.3) | 282 (10.9) |
| Unknown | 84 (4.8) | 76 (3.0) |
| Transmission risk group, n (%) | ||
| MSM | 1035 (59.1) | 1718 (66.7) |
| PWID | 150 (8.6) | 182 (7.1) |
| Heterosexuals | 387 (22.1) | 438 (17.0) |
| Other | 20 (1.1) | 88 (3.4) |
| Unknown | 159 (9.1) | 149 (5.8) |
| Origin, n (%) | ||
| Greek | 1321 (75.4) | 2117 (82.2) |
| Non-Greek | 197 (11.3) | 186 (7.2) |
| Unknown | 233 (13.3) | 272 (10.6) |
| Sampling period | ||
| 1999–2005 | 365 (20.8) | 1020 (39.6) |
| 2006–2010 | 503 (28.7) | 618 (24.0) |
| 2011–2015 | 880 (50.3) | 934 (36.3) |
| Unknown | 3 (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 2Demographic characteristics of PLHIV who were infected with subtypes A1 and B within and outside molecular transmission clusters (MTCs).
| Subtype A1 | Subtype B | |||
|---|---|---|---|---|
| Characteristic | Within |
Outside MTCs | Within |
Outside MTCs |
| Age, years, median (IQR) | 37 (29–46) 1 | 37 (31–48) 2 | 35 (29–42) 3 | 40 (34–48) 4 |
| Gender, n (%) | ||||
| Male | 1373 (83.6) | 62 (56.9) | 1729 (87.2) | 488 (82.3) |
| Female | 193 (11.8) | 39 (35.8) | 196 (9.9) | 86 (14.5) |
| Unknown | 76 (4.6) | 8 (7.3) | 57 (2.9) | 19 (3.2) |
| Transmission risk group, n (%) | ||||
| MSM | 1022 (62.2) | 13 (11.9) | 1373 (69.3) | 345 (58.2) |
| PWID | 124 (7.6) | 26 (23.9) | 148 (7.5) | 34 (5.7) |
| Heterosexuals | 339 (20.6) | 48 (44.0) | 321 (16.2) | 117 (19.7) |
| Other | 18 (1.1) | 2 (1.8) | 23 (1.1) | 65 (11.0) |
| Unknown | 139 (8.5) | 20 (18.4) | 117 (5.9) | 32 (5.4) |
| Origin, n (%) | ||||
| Greek | 1283 (78.1) | 38 (34.9) | 1635 (82.5) | 482 (81.3) |
| Non-Greek | 150 (9.2) | 47 (43.1) | 139 (7.0) | 47 (7.9) |
| Unknown | 209 (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) |
| Unknown | 3 (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 3Multivariable logistic regression estimates using the presence in molecular transmission clusters (MTCs) as the binary outcome variable.
| Subtype A1 1 | Subtype B 1 | |||||
|---|---|---|---|---|---|---|
| Explanatory Variable | Odds Ratio | 95% CI | p-Value | Odds Ratio | 95% CI | p-Value |
| Age | 0.99 | 0.97–1.01 | 0.341 | 0.96 | 0.95–0.97 | <0.001 |
| Gender | ||||||
| Male | 1.52 | 0.87–2.67 | 0.140 | 1.33 | 0.91–1.97 | 0.141 |
| Female 2 | ||||||
| Unknown | 2.80 | 0.64–12.3 | 0.172 | 0.84 | 0.22–3.16 | 0.796 |
| Transmission risk group | ||||||
| MSM | 5.34 | 2.52–11.31 | <0.001 | 0.98 | 0.69–1.38 | 0.893 |
| PWID | 0.42 | 0.22–0.80 | 0.008 | 0.80 | 0.49–1.30 | 0.364 |
| Heterosexuals 2 | ||||||
| Other | 0.86 | 0.17–4.47 | 0.861 | 0.10 | 0.05–0.18 | <0.001 |
| Unknown | 0.81 | 0.37–1.80 | 0.610 | 0.96 | 0.48–1.95 | 0.919 |
| Origin | ||||||
| Greek | 6.05 | 3.62–10.10 | <0.001 | 1.57 | 1.08–2.29 | 0.019 |
| Non-Greek 2 | ||||||
| Unknown | 1.27 | 0.62–2.60 | 0.516 | 1.26 | 0.74–2.12 | 0.394 |
| Sampling period | ||||||
| [1999–2005] 2 | ||||||
| [2006–2010] | 0.51 | 0.25–1.01 | 0.054 | 2.33 | 1.80–3.02 | <0.001 |
| [2011–2015] | 0.75 | 0.38–1.48 | 0.411 | 3.83 | 2.95–4.96 | <0.001 |
| Unknown | 1 | - | - | 1 | - | - |
| Sampling area | ||||||
| Central/Southern Greece 2 | ||||||
| Northern Greece | 0.91 | 0.54–1.55 | 0.736 | 1.32 | 0.96–1.81 | 0.092 |
| Crete | 1 | - | - | 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 4Migration 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 from | Central/Southern Greece | Northern Greece | Crete |
| Central/Southern Greece | 90.9 | 11.6 | |
| Northern Greece | 35.4 1 | 0.2 | |
| Crete | 0.1 1 | 0 | |
| Subtype B | |||
| Central/Southern Greece | 67.6 | 15.0 | |
| Northern Greece | 17.4 1 | 1.4 1 | |
| Crete | 1.2 1 | 0.5 1 | |
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.
© 2020 by the authors.
Details
; Gova, Maria 1 ; Adamis, Georgios 2 ; Xylomenos, Georgios 2 ; Chini, Maria 3 ; Mangafas, Nikos 3 ; Lazanas, Marios 3 ; Metallidis, Simeon 4
; Tsachouridou, Olga 4 ; Papastamopoulos, Vasileios 5 ; Chatzidimitriou, Dimitrios 6
; Kakalou, Eleni 5 ; Antoniadou, Anastasia 7 ; Papadopoulos, Antonios 7 ; Psichogiou, Mina 8 ; Basoulis, Dimitrios 8 ; Pilalas, Dimitrios 6
; Papageorgiou, Ifigeneia 1 ; Paraskeva, Dimitra 9 ; Chrysos, Georgios 9
; Paparizos, Vasileios 10 ; Kourkounti, Sofia 10 ; Sambatakou, Helen 11 ; Bolanos, Vasileios 11 ; Sipsas, Nikolaos V 12
; Lada, Malvina 13
; Barbounakis, Emmanouil 14 ; Kantzilaki, Evrikleia 14 ; Panagopoulos, Periklis 15
; Petrakis, Vasilis 15 ; Drimis, Stelios 9 ; Gogos, Charalambos 16 ; Hatzakis, Angelos 1 ; Beloukas, Apostolos 17
; Skoura, Lemonia 6 ; Paraskevis, Dimitrios 1
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; Department of Biomedical Sciences, School of Health Sciences, University of West Attica, 12243 Athens, Greece