PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2017; 12(2): e0172473.
Published online 2017 February 24. doi:  10.1371/journal.pone.0172473
PMCID: PMC5325290

Expanding HIV-1 subtype B transmission networks among men who have sex with men in Poland

Zhefeng Meng, Editor

Abstract

Introduction

Reconstruction of HIV transmission links allows to trace the spread and dynamics of infection and guide epidemiological interventions. The aim of this study was to characterize transmission networks among subtype B infected patients from Poland.

Material and methods

Maximum likelihood phylogenenetic trees were inferred from 966 HIV-1 subtype B protease/reverse transcriptase sequences from patients followed up in nine Polish HIV centers. Monophyletic clusters were identified using 3% within-cluster distance and 0.9 bootstrap values. Interregional links for the clusters were investigated and time from infection to onward transmission estimated using Bayesian dated MCMC phylogeny.

Results

Three hundred twenty one (33.2%) sequences formed 109 clusters, including ten clusters of ≥5 sequences (n = 81, 8.4%). Transmission networks were more common among MSM (234 sequences, 68.6%) compared to other infection routes (injection drug use: 28 (8.2%) and heterosexual transmissions: 59 (17.3%) cases, respectively [OR:3.5 (95%CI:2.6–4.6),p<0.001]. Frequency of clustering increased from 26.92% in 2009 to 50.6% in 2014 [OR:1.18 (95%CI:1.06–1.31),p = 0.0026; slope +2.8%/year] with median time to onward transmission within clusters of 1.38 (IQR:0.59–2.52) years. In multivariate models clustering was associated with both MSM transmission route [OR:2.24 (95%CI:1.38–3.65),p<0.001] and asymptomatic stage of HIV infection [OR:1.93 (95%CI:1.4–2.64),p<0.0001]. Additionally, interregional networks were linked to MSM transmissions [OR:4.7 (95%CI:2.55–8.96),p<0.001].

Conclusions

Reconstruction of the HIV-1 subtype B transmission patterns reveals increasing degree of clustering and existence of interregional networks among Polish MSM. Dated phylogeny confirms the association between onward transmission and recent infections. High transmission dynamics among Polish MSM emphasizes the necessity for active testing and early treatment in this group.

Introduction

Scale up of antiretroviral treatment (ART), mother-to-child transmission reduction programs and earlier introduction of ART have resulted in decreased transmission of new HIV infections worldwide, especially in Africa and Asia [1,2]. However, in the recent years increase in the number of the new cases has been observed in high and middle income European countries.

This increase was associated with injection drug use (IDU) in Eastern Europe including countries such as Russia, Ukraine and Belarus and upsurge of new infections among men-who-have-sex-with-men (MSM) in West European countries [3,4]. HIV-1 epidemics in Poland was historically associated with IDU, which was the predominant route of infection until 2005 [5]. However, in recent years a shift in the dominant transmission route was observed with a 14-fold increase in the infection frequency among MSM—from 2.5 new diagnoses per million of male inhabitants in 2000 to 33.8 in 2011 and beyond [6] as well as a 90% decline in the IDU related transmissions [7]. Currently more than 70% of new diagnoses are amongst MSM aged 25–34 years [8]. Also, recent models have indicated a high proportion (68.3%, 95%CI 53.9–76.1%) of undiagnosed HIV among MSM, which is related to low testing rates in Poland and lack of uniform preventive strategy focused on this group, despite higher testing rates associated with MSM exposure [9,10].

It has been reported that new infections were associated with inconsistent condom use and concentrated in the large cities, however, factors driving this epidemics in Poland remain largely understudied [6]. Across HIV infected populations subtype B remained stably the most common HIV-1 variant (95.8–86.1%) [11,12] and associated with male gender, MSM transmissions, higher CD4 lymphocyte counts and lower HIV plasma viremia at care entry [13,14].

Phylogenetics provides an insight into the HIV network structure including its transmission dynamics [15]. Use of dated sequence data, coupled with demographic and clinical information allows to model characteristics of the epidemics including population-based patterns of HIV transmission [16]. Dated sequence analyses from cohorts with known and unknown infection time, showed that early infection was associated with onward transmissions in 25–50% of cases [17,18,19,20]. For above-mentioned purposes HIV-1 partial pol sequences obtained from routine genotyping assays have been successfully used. Characterizing populations forming transmission networks allows to target interventions to the ones at risk. Using phylogeny relationship between various viral strains may be established with identification of closely interrelated subsets—phylogenetic clusters [19]. Typically, a cluster contains all taxa—in case of HIV analyses patient derived sequences—descending from the similar lineage or ancestry [21]. Using these clusters chronology of the spread viral clades spread may be visualized. Typically clusters are defined by the genetic discane between analyzed sequences and/or statistical support (bootstrap) within the inferred tree [22,23,24]. In HIV medicine clustering dynamics may reflect changing patient mobility and allows to assess the efficacy of ART as prevention strategy [19,25,26].

Having previously identified common clustering among MSM [13], in this study we wished to investigate interregional links and time-trends associated with transmission networks for the HIV-1 subtype B. The aim of this study was to provide insight, based on the phylogenetic data, into the characteristics and trends in the dynamics of HIV transmission in Poland. We also aimed to investigate network patterns on the national level to guide HIV prevention.

Materials and methods

Study population and sequencing

For this study 966 HIV-1 subtype B samples collected from patients with defined transmission route for the purpose of genotypic drug resistance analyses were used. Patients were linked to care in 9 of 17 Polish HIV treatment centres (Białystok, Bydgoszcz, Gdańsk, Łódź, Kraków, Poznan, Szczecin, Wroclaw, Zielona Gora). Sequencing was performed at the Clinical Laboratory at the Department of Infectious, Tropical Diseases and Acquired Immune Deficiency, Pomeranian Medical University in Szczecin, Poland.

The samples included sequences collected both from treatment-experienced (n = 176, 18.21%) and antiretroviral-naive (n = 790, 81.87%) patients diagnosed in Poland from 1989 to 2014, prior to the introduction of the antiretroviral treatment. Data for subtype B treatment-naive patients were partially collected as a part of the Northern-Poland Resistance Transmission Tracking for HIV-1 (NoRTTH) project funded by the Polish AIDS Society, and the resistance data were previously published elsewhere [13]. For the molecular clock analyses sequence sampling dates were used. For individuals with multiple sequences available, the earliest sequence after the positive confirmation HIV test, was included.

The analysed data included: gender, CD4 lymphocyte count and HIV viral load at care entry, CD4 lymphocyte nadir, age at diagnosis (the age of first positive confirmation test), HCV co-infection (defined as presence of anti-HCV antibody or HCV-RNA), transmission route (self-defined by the patient) and CDC clinical stage according to the 1993 classification system for HIV infection [27] at care entry. As seroconversion time point was often unavailable in the source documentation data for the acute HIV infections were not collected—date of infection was based on the date of the first positive and confirmed HIV test.

The study protocol was approved by the institutional review board named Bioethical Committee of Pomeranian Medical University in Szczecin (Komisja Bioetyczna Pomorskiego Uniwersytetu Medycznego w Szczecinie), Poland (approval number KB-0012/08/12). Research was conducted in accordance with the Declaration of Helsinki. As the patient data were coded and anonymous, while all the procedures, including genotyping, are the part of the medical care on the people living with HIV, no separate written consent was obtained but physicians informed all subjects on the planned research and checked for verbal non-opposition from their patients.

Sequencing

HIV RNA extraction, as well as reverse transcriptase and protease genotyping, were performed using a genotyping assay (Viroseq 2.8, Abbott Molecular, Abbott Park, Illinois, USA) according to the manufacturer’s protocol. Protease/reverse transcriptase sequences obtained were 1302 b.p. long; location from the start of HXB2 genome: positions 2253 to 3525. Identification of subtype B was performed based on the partial pol sequences using two methods: on-line software (REGA genotyping 2.0 tool; http://bioafrica.mrc.ac.za/rega-genotype/html/subtypinghiv.html) and phylogenetic analyses using maximum likelihood (ML) method under the GTR+γ model with approximate likelihood ratio test (aLRT) (online PHYMLv.3.0-http://www.atgc-montpellier.fr/phyml) [28]. For reference, sequences with known subtype from the 2014 version of the HIV sequence compendium (Los Alamos National Laboratory) were used, supplemented with CRFs from the HIV sequence database (http://www.hiv.lanl.gov/components/sequence/HIV/search/search.html). All sutype B sequences included in the analysis were assessed for recombination using two methods: a simplot bootscan (2-parameter Kimura model, window size 200 b.p., 20.p. step) [29] and jumping-profile hidden Markov model (jpHMM) [30]. Recombinant sequences were removed from the final dataset. Thirty-seven codon sites associated with protease and reverse transcriptase resistance mutations were manually deleted from the sequence alignments. The final alignment was 1191 b.p. long.

Identification of clustering

First, clusters were identified using the Maximum Likelihood (ML) method with the NNI-SPR sub-tree algorithm under the GTR+γ with the PHYML v 3.0 using automatic model selection [28]. Clustering was assessed by Cluster Picker software with the maximum genetic distances calculated by the program [31]. Monophyletic clusters were identified with Cluster-picker software using 3% within-cluster distance (substitutions/site) threshold and 0.9 aLRT bootstrap values. Phylogeographic interregional links within the clustered sequences were investigated.

Additionally, for the 321 clustered sequences Bayesian Monte Carlo Marcov Chain (MCMC) was used. Pooling of clustered sequences was performed for computational reasons. Three replicates of 35 million generations were run in BEAST v2.3.1 (2015) [32] using a coalescent constant population prior and a GTR model with uncorrelated relaxed exponential molecular clock [33]. This model was cosidered the most fit using phyML modeling. The molecular clock was calibrated using parameters derived from the phyML analysis, namely gamma shape: 0.33, estimated proportion of invariant sites with the initial value of 0.5, and the following GTR substitution rates: AC: 2.621, AG: 7.422, AT: 0.922, GC: 0.943, CT: 14.36, GT: 1.0. Lognormal priors were applied to the coalescent population size, proportion invariants sites for partitions and clock mean on partitions. All prior and posterior effective sample size (ESS) values exceeded 200. A consensus tree with posterior probabilities for branch support was obtained and annotated with TreeAnnotator v 1.5.4 using a 10% burnin. All trees were visualized in Figtree v.1.2.2.

Time between transmission events was calculated from the dated MCMC tree using distances between MRCAs for each cluster (internode intervals)—for each identified cluster/pair the time to the most recent common ancestor (tMRCA) was calculated. Medians and interquartile ranges for these intervals were calculated separately both for all clusters and clusters ≥5 sequences.

Statistics

Statistical comparisons were performed using Fisher’s exact and Chi2 tests for nominal variables. Continuous variables were analysed using the Mann-Whitney U-test for nonparametric statistics (HIV viral load, age at diagnosis, baseline and nadir CD4 lymphocyte counts, internode intervals). Commercial software (Statistica 10PL, Statasoft, Warsaw, Poland) was used for above statistical calculations. Time trends for the frequency of transmission pairs and clusters were calculated using the logistic regression (R statistical platform version 3.1.0, package MASS). Two separate analyses (pairs + clusters ≥3 sequences and clusters ≥3 sequences only) were performed, with presence (coded as 1) or lack of clusters (coded as 0) investigated for every year. For analyses of time trends in the clustering the 707 (73.2% of the total cohort) patients diagnosed from 2009–2014 were included, as the sample size for the years 1989–2008 did not exceed 50 sequences. This sample (707 cases) represents 10.8% of the new HIV-1 diagnoses in Poland for the years 2009–2014 (6560 cases).

Sequence data

Sequences from this study have been submitted to GenBank and may be accessed with the following IDs: KC409158-KC409222, KM057325-KM057372, KM283892-KM284680, KT340108- KT340205.

Results

Clustering among subtype B infections

Based on the constructed a maximum likelihood phylogenetic tree a total of 109 clusters and pairs containing 321 (33.2%) sequences were found, including 36 clusters of ≥3 individuals incorporating 165 (17.1%) sequences (S1 Fig). There were ten clusters of ≥5 sequences (n = 81, 8.4%), including two larger ones comprising 13 and 20 patients. All these clusters contained only individuals infected by the MSM route (shaded in grey at S1 Fig). Frequency of transmission pairs and clusters notably increased over time: from 26.92% in 2009 to 50.6% in 2014 [OR: 1.18 (95%CI: 1.06–1.31), p = 0.0026; slope for the frequency increase 2.8%/year]. This trend was similar when sequence pairs were excluded from the clustering analyses (increase from 9.62% in 2009 to 25.3% in 2014 [OR: 1.21 (95%CI: 1.06–1.38), p = 0.004; slope for the frequency increase 2.2%/year] (Fig 1).

Fig 1
Proportions and logistic regression curves over time (2009–2014) for clusters/pairs (A, left) and clusters ≥3 sequences (B, right).

Clinical characteristics of clustered cases

Univariate associations between clustering and group characteristics are presented in Table 1. Sequence clusters/pairs were notably more frequent among MSM (229 cases, 71.3%) compared to other transmission routes (injection drug use (IDU): 32 (10.0%) and heterosexual transmissions: 60 (18.7%) cases, respectively [p<0.001, OR: 3.5 (95%CI: 2.6–4.6)]. Furthermore, clusters (groups of ≥ 3 cases) predominantly contained MSM-derived sequences (n = 131, 79.4%) with only 24 (14.5%) ones from patients infected by heterosexual contact and 10 (6.1%) by IDUs route. Moreover, patients with sequences forming pairs and clusters were younger at HIV diagnosis [median: 31 (IQR: 26–37) compared to 32 (27–39) years for unclustered sequences, p = 0.023], mostly asymptomatic at the time of analysis (67.8%, p<0.001), predominantly antiretroviral treatment naive (90.0%, p<0.001), and free of hepatitis C virus coinfection (83.3%, p<0.001). Additionally, CD4 lymphocyte counts at genotyping and CD4 lymphocyte nadirs were significantly higher for clustered/paired patients [baseline CD4 median: 407 (IQR: 209–565) cells/μl versus 290 (111–467) cells/μl, p<0.001 and nadir CD4 median 356 (IQR: 187–493) versus 223 (IQR: 88–383) cells/μl, p<0.001 for paired/clustered and unclustered sequences, respectively], while HIV viral load was similar across these categories. Those clinical, immunologic and virologic patterns were similar for clusters of ≥ 3 cases (after removal of sequence pairs from the analysis).

Table 1
Demographic, clinical and laboratory cohort characteristics associated with pairs and clusters ≥ 3 cases.

In multivariate analyses of clusters/pairs adjusted for transmission route, gender, clinical category at diagnosis, antiretroviral treatment history (treated vs. naive) and HCV status the factors associated with clustering were MSM transmission route [OR: 2.24 (95%CI: 1.38–3.65), p<0.001] and CDC category A [OR: 1.93 (95%CI: 1.4–2.64), p<0.0001]. Similar results were obtained for clusters of ≥ 3 cases with OR: 2.093 (95%CI: 1.14–3.85), p = 0.018 for MSM transmissions and OR: 2.64 (95%CI: 1.71–4.04), p<0.0001 for asymptomatic infection at HIV diagnosis.

Transmission networks

To reflect the dynamics of the HIV transmissions in Poland in recent years we analysed the transmissions networks among the sites participating in the study. For this purpose the clusters and pairs containing sequences from patients followed-up at different centers were identified. Such interregional networks were found in 31 clusters/pairs containing 126 (13%) sequences from various centers (Fig 2) and were more common among MSM compared to other, transmission routes. In clusters/pairs containing patients from ≥2 centres 112 (11.6%) sequences were derived from MSM patients compared to four (0.4%) from IDU and 10 (1%) from heterosexually infected cases [p<0.001, OR:4.7 (95%CI: 2.55–8.96))]. As expected, frequency of clusters containing interregional sequences increased notably with the cluster size (OR for trend: 2.94 (95% CI: 1.51–5.64), p<0.001) (Fig 3).

Fig 2
Interregional transmission networks identified in the Bayesian MCMC analysis.
Fig 3
Cluster frequencies by size and type of clustering.

Time intervals between clustered transmission events

To further investigate phylogenetic relationships between the clusters and pairs dated Bayesian MCMC phylogeny using a subset of 321 closely related sequences forming clusters and pairs was reconstructed (marked in blue on the S1 Fig). Median genetic distances for clusters/pairs were 0.0199 (IQR: 0.013–0.0237) and 0.023 (IQR: 0.0199–0.0276) for clusters ≥ 3 sequences.

To reflect the time between transmission events we have inferred a dated tree and calculated distances between tMRCA and previous nodes [20] for clusters of ≥ 3 sequences (Fig 4). Separate analysis was performed for clusters of ≥ 5 cases (exact topology for the time-scaled phylogenies for the clusters ≥ 5 shown in Fig 4). Time between the estimated times of infection for the clusters, reflecting time to onward transmission, ranged from 0.05 to 6.08 years [median 1.38 (IQR: 0.59–2.52) years]. For clusters of ≥5 sequences the internode intervals were significantly shorter compared to the smaller clusters [median 1.1 (IQR: 0.43–1.9) years versus median 2.13 (IQR: 1.1–3.66) years, p = 0.0026, S2 Fig]. Distribution of the internode distances for the clusters indicates that 22.4% of clustered onward transmissions occurred within the first 6 months, while 37.6% within the first year from infection; Analyzing only MSM clusters, the model suggests that the 23.1% and 35.9% of clustered transmissions occur within the first 6 months and one year from index case infection, respectively. This indicates that at least 3.9% (7.0% for MSM) of transmissions in the studied group occurred within 6 months and 6.6% (10.9% for MSM) within the first year of infection.

Fig 4
Bayesian, time-annotated MCMC phylogenetic tree of the subtype B pol sequences derived from clustered sequences coloured by cluster size.

Discussion

In the light of the rising number of new HIV diagnoses, clustered transmissions among MSM in Europe [4,34] and the shift in the character of epidemics in Poland from injection drug use related to concentrated, sexually-driven epidemics [6] we focused this study on the phylodynamic analysis of transmission networks. Data from different centers were analysed to minimize the sampling bias related to the localized epidemics and to reflect countrywide trends in the molecular epidemiology of HIV-1. The dataset included HIV sequences collected for genotyping purposes, and was limited to the subtype B as this clade remains the most prevalent in the country [11,12,35]. We have shown that the frequency (31%) of clustering is associated with MSM transmissions and early stage of HIV infection with common (13%) interregional transmission networks reflecting patient mobility Moreover, it was demonstrated that clustering is becoming more frequent over time, which indicates the uncontrolled, expanding epidemics and forward infections from individuals probably unaware of their HIV status, especially among MSM populations. This is in line with the surveillance data in Poland with the rapid increase in new HIV diagnoses in MSM, especially aged 24–34 years. In this group the number of newly diagnosed infections increased >20 fold, from 13 in the year 2000, to 313 in 2011 [6]. As of the most recent annual report (2013) 68.4% of newly diagnosed infections with disclosed transmission category were observed in MSM [7]. Clustering may only partially be explained by the earlier diagnosis during the course of infection over time, as the percentage of the late presenters reported from the epidemiological data was increasing [7].

Association between MSM transmission route and clustering within the subtype B infected cases has been previously reported in various cohorts [18,19,25,36,37,38]. Clinical characteristics of the clustered sequences followed the patterns observed in other studies including lack of prior antiretroviral treatment exposure, higher CD4 lymphocyte counts, earlier stage of HIV disease and more recent sampling [16,21]. Also observed clustering frequencies were within the 18.2–60% range observed across European and North-American datasets [16,36,37,39,40] and similar to the 31.2% observed in SPREAD cohort using the same internode genetic distance threshold of 3% [41]. Transmission networks observed in our study were not limited to the single geographical site with 13% of the total analyzed sequences (39.25% of the clustered ones) dispersed across HIV treatment centers. Of note, this may be associated with clustering networks, but also reflects the patient migration for personal and professional reasons.

It should be noted, that in the aforementioned SPREAD report, within a much smaller dataset of 193 Polish sequences from different geographical sites than the current study, a large cluster of 15 sequences was found. This is consistent with the identification of large, interregional clusters in our analysis including the ones comprising of 13 and 20 sequences. These large clusters may reflect rapidly growing sub-epidemics [19] and have been found throughout subtype B infected populations worldwide [38,42,43,44,45], however the exact transmission pattern within the transmission clusters is difficult to analyze in the cohort of chronically infected individuals. As expansion of the clusters over time has been reported [19,46] this finding supports the need for early introduction of antiretroviral therapy in populations with frequent clustering aiming to break the transmission chains [47].

Additionally, we have modeled the onward transmission time for the clustered infections indicating that 50% of clustered transmissions occur within 16.5 months after infection, with 25% within 7 months and 75% within 30 months.

Although interval modeling was performed using data inferred from samples with unknown duration of HIV infection, the results are highly consistent with the ones obtained for the Dutch observational cohort Athena [18] with reported median time to onward transmissions of 17 months (IQR: 7–17 months). Moreover, our data are in accordance to the study by Lewis et al., indicating that 22.4% (23.1% among MSM) of new infections occurred within the first 6 months from the index case [20]. Also, our observation that in the clusters ≥5 sequences 25% percent of transmissions occurred during the first 5.3 months and 50% of transmissions within the 13.2 months is consistent with the report by Brenner et al., modeling that 25%–30% of clustered transmissions occur within 6-months and 50% over a 14- to 17-month period [19].

These data, coupled with the temporal increase in the number of small and larger clusters, confirms the role of undiagnosed, recently infected individuals in fuelling of the current HIV epidemics in Poland, similarly to Canada, UK, Netherlands or France [18,19,25,48]. Our study is expanding the evidence on the enlarging MSM epidemics in Poland with the increasing importance of transmission networks and shows that current molecular dynamics of HIV spread reflects that of the West-European countries. One of the consequences of expanding transmission networks may be the increase in frequency of primary drug resistance [42,44], however, we have shown previously that trends in drug resistance were stable over time, with infrequent transmission within clusters and median prevalence of 8.3% [13,49]. Given the fact that estimated tMRCA infection times are within the timeframe of increased frequency of sexually-driven HIV infections observed over the last 10 years we suggest that clustered transmissions may significantly contribute to the current increase in the prevalence of HIV in Poland. These data are relevant for the design of prophylactic interventions and emphasize the need for prevention programs to focus on early detection and treatment among MSMs as it may reduce the frequency of onward transmissions within the networks [47,50]. Of note, prevention campaigns targeted to MSM in Poland are currently limited to anonymous counseling and testing, with no pre-exposure prophylaxis available. Scarce programs such as "Test and Keep in care"–program targeted at identification and referral to care for the newly diagnosed patients—do not seem sufficient as testing rates among populations at risk remain low [51]. The strategy should be "Seek, Test, Link, Treat, and Retain” as suggested by Hull at al. [52] and may be implemented across all populations, not only MSM.

Limitation of the study is related to the convenience of sampling—certain proportion of cases representing intermediate transmissions may have not been indentified and therefore not included in the analysis. This may result in the overestimation of the transmission intervals. Additionally, the data are limited only to the centers participating in the study—mostly Western and Southern provinces, however this is the largest sequence dataset available for the multiple sites in the country. Also, it should be emphasized that sequence datasets are not sufficient to determine the specific transmission events and ideally should be coupled with contact information [53]. It should also be noted, that definitions of genetic distance based clustering are inconsistent across the published studies ranging from 1.5% (most widely used for recent infections and local cohorts) to 4.5% (chronically infected individuals from nationwide analyses) [19,25,36,40,54]. We have used the 3% threshold for the intercluster genetic distances proposed by Kaye et al. [44] and used by the European Spread cohort [41].

Conclusions

Clustering networks reflect expanding HIV epidemics among MSM in Poland. Interregional clustering reflects high dynamics of the nationwide HIV epidemics in the MSM group. Presented molecular data allow to emphasize the need for the evidence-based interventions which should focus on the proactive identification of populations at the highest risk, testing and early referral to antiretroviral treatment. Pre-exposure antiretroviral drug use should be implemented as preventive strategy among MSM in Poland as this group contributes significantly to onward HIV transmissions.

Supporting information

S1 Fig

Maximum likelihood tree showing the indentified clusters marked in blue.

Clusters ≥ 5 sequences shaded in grey. All clustered sequences were used for the MCMC analysis using beast.

(PDF)

S2 Fig

Histograms presenting distributions in time between the dates of infection for the clustered sequences for all (a) sequences/pairs, (b) clusters ≥ 5 sequences, (c) clusters of 2–4.

(TIF)

Acknowledgments

The authors would like to express the gratitude to all of the patients participating in the study.

Funding Statement

The study was funded partially by the internal funds of Pomeranian Medical University, Szczecin, Poland, and partially by the Polish Scientific AIDS Society funds. MP, MLP, MWJ and UA received the salary from Pomeranian Medical University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

Sequences from this study have been submitted to GenBank and may be accessed with the following IDs: GU906857-GU906877, KC409134-KC409222, KM057325-KM057372, KM283892-KM284680, KT340108-KT340205.

References

1. Hill A, Pozniak A (2015) HIV treatment cascades: how can all countries reach the UNAIDS 90-90-90 target? AIDS 29: 2523–2525. doi: 10.1097/QAD.0000000000000864 [PubMed]
2. Roehr B (2015) UNAIDS celebrates success in halting spread of HIV and sets goals for 2030. BMJ 351: h3832 doi: 10.1136/bmj.h3832 [PubMed]
4. Beyrer C, Sullivan P, Sanchez J, Baral SD, Collins C, Wirtz A,.et al. (2013) The increase in global HIV epidemics in MSM. AIDS 27: 2665–2678. doi: 10.1097/01.aids.0000432449.30239.fe [PubMed]
5. Rosinska M (2006) Current trends in HIV/ AIDS epidemiology in Poland, 1999–2004. Euro surveillance: bulletin europeen sur les maladies transmissibles = European communicable disease bulletin 11.
6. Rosinska M, Janiec J, Niedzwiedzka-Stadnik M (2013) Increase of new HIV diagnoses among men who have sex with men in Poland, 2000 to 2011. Euro surveillance: bulletin europeen sur les maladies transmissibles = European communicable disease bulletin 18: 20642. [PubMed]
7. Niedzwiedzka-Stadnik M, Pielacha M, Rosinska M (2015) HIV and AIDS in Poland in 2013. Przeglad epidemiologiczny 69: 267–272, 381–264. [PubMed]
8. Niedzwiedzka-Stadnik M, Pielacha M, Rosinska M (2014) HIV and AIDS in Poland in 2012. Przeglad epidemiologiczny 68: 283–289, 383–286. [PubMed]
9. Rosinska M, Gwiazda P, DEA D, Presanis AM (2015) Bayesian evidence synthesis to estimate HIV prevalence in men who have sex with men in Poland at the end of 2009. Epidemiology and infection: 1–17. [PMC free article] [PubMed]
10. Rosinska M, Simmons R, Marzec-Boguslawska A, Janiec J, Porter K (2016) Relating HIV testing patterns in Poland to risky and protective behaviour. AIDS care 28: 423–431. doi: 10.1080/09540121.2015.1100702 [PubMed]
11. Smolen-Dzirba J, Rosinska M, Kruszynski P, Bratosiewicz-Wasik J, Janiec J, Beniowski M, et al. (2012) Molecular epidemiology of recent HIV-1 infections in southern Poland. Journal of medical virology 84: 1857–1868. doi: 10.1002/jmv.23395 [PubMed]
12. Stanczak GP, Stanczak JJ, Marczynska M, Firlag-Burkacka E, Wiercinska-Drapalo A, Leszczyszyn-Pynka M, et al. (2010) Evolving patterns of HIV-1 transmitted drug resistance in Poland in the years 2000–2008. Journal of medical virology 82: 1291–1294. doi: 10.1002/jmv.21782 [PubMed]
13. Parczewski M, Leszczyszyn-Pynka M, Witak-Jedra M, Maciejewska K, Rymer W, Szymczak A, et al. (2015) Transmitted HIV drug resistance in antiretroviral-treatment-naive patients from Poland differs by transmission category and subtype. The Journal of antimicrobial chemotherapy 70: 233–242. doi: 10.1093/jac/dku372 [PubMed]
14. Parczewski M, Leszczyszyn-Pynka M, Bander D, Urbanska A, Stanczak G, Boron-Kaczmarska A.. (2010) Characteristics of HIV-1 non-B subtype infections in Northwest Poland. Journal of medical virology 82: 1306–1313. doi: 10.1002/jmv.21797 [PubMed]
15. Grabowski MK, Redd AD (2014) Molecular tools for studying HIV transmission in sexual networks. Current opinion in HIV and AIDS 9: 126–133. doi: 10.1097/COH.0000000000000040 [PMC free article] [PubMed]
16. Aldous JL, Pond SK, Poon A, Jain S, Qin H, Kahn JS,. (2012) Characterizing HIV transmission networks across the United States. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America 55: 1135–1143. [PMC free article] [PubMed]
17. Brenner BG, Roger M, Routy JP, Moisi D, Ntemgwa M, Matte C, et al. (2007) High rates of forward transmission events after acute/early HIV-1 infection. The Journal of infectious diseases 195: 951–959. doi: 10.1086/512088 [PubMed]
18. Bezemer D, van Sighem A, Lukashov VV, van der Hoek L, Back N, Schuurman R. et al. (2010) Transmission networks of HIV-1 among men having sex with men in the Netherlands. AIDS 24: 271–282. doi: 10.1097/QAD.0b013e328333ddee [PubMed]
19. Brenner BG, Wainberg MA (2013) Future of phylogeny in HIV prevention. Journal of acquired immune deficiency syndromes 63 Suppl 2: S248–254. [PMC free article] [PubMed]
20. Lewis F, Hughes GJ, Rambaut A, Pozniak A, Leigh Brown AJ (2008) Episodic sexual transmission of HIV revealed by molecular phylodynamics. PLoS medicine 5: e50 doi: 10.1371/journal.pmed.0050050 [PubMed]
21. Volz EM, Koopman JS, Ward MJ, Brown AL, Frost SD (2012) Simple epidemiological dynamics explain phylogenetic clustering of HIV from patients with recent infection. PLoS computational biology 8: e1002552 doi: 10.1371/journal.pcbi.1002552 [PMC free article] [PubMed]
22. Buonaguro L, Petrizzo A, Tagliamonte M, Vitone F, Re MC, Pilotti E, et al. (2008) Molecular and phylogenetic analysis of HIV-1 variants circulating in Italy. Infectious agents and cancer 3: 13 doi: 10.1186/1750-9378-3-13 [PMC free article] [PubMed]
23. Hue S, Pillay D, Clewley JP, Pybus OG (2005) Genetic analysis reveals the complex structure of HIV-1 transmission within defined risk groups. Proceedings of the National Academy of Sciences of the United States of America 102: 4425–4429. doi: 10.1073/pnas.0407534102 [PubMed]
24. Leitner T, Escanilla D, Franzen C, Uhlen M, Albert J (1996) Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis. Proceedings of the National Academy of Sciences of the United States of America 93: 10864–10869. [PubMed]
25. Leigh Brown AJ, Lycett SJ, Weinert L, Hughes GJ, Fearnhill E, Dunn D. (2011) Transmission network parameters estimated from HIV sequences for a nationwide epidemic. The Journal of infectious diseases 204: 1463–1469. doi: 10.1093/infdis/jir550 [PMC free article] [PubMed]
26. Paraskevis D, Pybus O, Magiorkinis G, Hatzakis A, Wensing AM, van de Vijver D A,. (2009) Tracing the HIV-1 subtype B mobility in Europe: a phylogeographic approach. Retrovirology 6: 49 doi: 10.1186/1742-4690-6-49 [PMC free article] [PubMed]
27. (1993) From the Centers for Disease Control and prevention. 1993 revised classification system for HIV infection and expanded surveillance case definition for AIDS among adolescents and adults. JAMA: the journal of the American Medical Association 269: 460. [PubMed]
28. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic biology 59: 307–321. doi: 10.1093/sysbio/syq010 [PubMed]
29. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. (1999) Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. Journal of virology 73: 152–160. [PMC free article] [PubMed]
30. Chan PA, Reitsma MB, DeLong A, Boucek B, Nunn A, Salemi M, et al. (2014) Phylogenetic and geospatial evaluation of HIV-1 subtype diversity at the largest HIV center in Rhode Island. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases 28: 358–366. doi: 10.1016/j.meegid.2014.03.027 [PMC free article] [PubMed]
31. Ragonnet-Cronin M, Lee BE, Plitt SS, Zahariadis G, Merks H, et al. (2013) Baseline clinical HIV genotypes are a valid measure of transmitted drug resistance within the treatment-naive population. Journal of acquired immune deficiency syndromes 64: 443–447. doi: 10.1097/QAI.0b013e3182a4b991 [PubMed]
32. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC evolutionary biology 7: 214 doi: 10.1186/1471-2148-7-214 [PMC free article] [PubMed]
33. Drummond AJ, Ho SY, Phillips MJ, Rambaut A (2006) Relaxed phylogenetics and dating with confidence. PLoS biology 4: e88 doi: 10.1371/journal.pbio.0040088 [PubMed]
34. Vermund SH, Leigh-Brown AJ (2012) The HIV Epidemic: High-Income Countries. Cold Spring Harbor perspectives in medicine 2: a007195 doi: 10.1101/cshperspect.a007195 [PMC free article] [PubMed]
35. Parczewski M, Witak-Jedra M, Maciejewska K, Bociaga-Jasik M, Skwara P, Garlicki A, et al. (2014) Time trends in HIV-1 transmitted drug resistance mutation frequency in Poland. Journal of the International AIDS Society 17: 19753 doi: 10.7448/IAS.17.4.19753 [PMC free article] [PubMed]
36. Chalmet K, Staelens D, Blot S, Dinakis S, Pelgrom J, Pelgrom J, et al. (2010) Epidemiological study of phylogenetic transmission clusters in a local HIV-1 epidemic reveals distinct differences between subtype B and non-B infections. BMC infectious diseases 10: 262 doi: 10.1186/1471-2334-10-262 [PMC free article] [PubMed]
37. Chan PA, Hogan JW, Huang A, DeLong A, Salemi M, Mayer KH, et al. (2015) Phylogenetic Investigation of a Statewide HIV-1 Epidemic Reveals Ongoing and Active Transmission Networks Among Men Who Have Sex With Men. Journal of acquired immune deficiency syndromes 70: 428–435. doi: 10.1097/QAI.0000000000000786 [PMC free article] [PubMed]
38. Grossman Z, Avidor B, Mor Z, Chowers M, Levy I, Shahar E, et al. (2015) A Population-Structured HIV Epidemic in Israel: Roles of Risk and Ethnicity. PloS one 10. [PMC free article] [PubMed]
39. Cuevas M, Fernandez-Garcia A, Sanchez-Garcia A, Gonzalez-Galeano M, Pinilla M, Sanchez-Martinez M, et al. (2009) Incidence of non-B subtypes of HIV-1 in Galicia, Spain: high frequency and diversity of HIV-1 among men who have sex with men. Euro surveillance: bulletin europeen sur les maladies transmissibles = European communicable disease bulletin 14. [PubMed]
40. Yebra G, Holguin A, Pillay D, Hue S (2013) Phylogenetic and demographic characterization of HIV-1 transmission in Madrid, Spain. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases 14: 232–239. doi: 10.1016/j.meegid.2012.12.006 [PubMed]
41. Frentz D, Wensing AM, Albert J, Paraskevis D, Abecasis AB, Hamouda O, et al. (2013) Limited cross-border infections in patients newly diagnosed with HIV in Europe. Retrovirology 10: 36 doi: 10.1186/1742-4690-10-36 [PMC free article] [PubMed]
42. Brenner BG, Roger M, Moisi DD, Oliveira M, Hardy I, Turgel R, et al. (2008) Transmission networks of drug resistance acquired in primary/early stage HIV infection. AIDS 22: 2509–2515. doi: 10.1097/QAD.0b013e3283121c90 [PMC free article] [PubMed]
43. Brenner BG, Roger M, Stephens D, Moisi D, Hardy I, Weinberg J, et al. (2011) Transmission clustering drives the onward spread of the HIV epidemic among men who have sex with men in Quebec. The Journal of infectious diseases 204: 1115–1119. doi: 10.1093/infdis/jir468 [PMC free article] [PubMed]
44. Kaye M, Chibo D, Birch C (2008) Phylogenetic investigation of transmission pathways of drug-resistant HIV-1 utilizing pol sequences derived from resistance genotyping. Journal of acquired immune deficiency syndromes 49: 9–16. doi: 10.1097/QAI.0b013e318180c8af [PubMed]
45. Little SJ, Kosakovsky Pond SL, Anderson CM, Young JA, Wertheim JO, Mehta SR, et al. (2014) Using HIV networks to inform real time prevention interventions. PloS one 9. [PMC free article] [PubMed]
46. Vega Y, Delgado E, Fernandez-Garcia A, Cuevas MT, Thomson MM, Montero V, et al. (2015) Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004–2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations. PloS one 10. [PMC free article] [PubMed]
47. Wertheim JO, Kosakovsky Pond SL, Little SJ, De Gruttola V (2011) Using HIV transmission networks to investigate community effects in HIV prevention trials. PloS one 6: e27775 doi: 10.1371/journal.pone.0027775 [PMC free article] [PubMed]
48. Frange P, Meyer L, Deveau C, Tran L, Goujard C, Ghosn J, et al. (2012) Recent HIV-1 infection contributes to the viral diffusion over the French territory with a recent increasing frequency. PloS one 7: e31695 doi: 10.1371/journal.pone.0031695 [PMC free article] [PubMed]
49. Parczewski M, Urbanska A, Maciejewska K, Witak-Jedra M, Leszczyszyn-Pynka M (2014) Transmitted drug resistance to rilpivirine among antiretroviral-naive patients living with HIV from northern Poland. Journal of the International AIDS Society 17: 18929 doi: 10.7448/IAS.17.1.18929 [PMC free article] [PubMed]
50. Yeghiazarian L, Cumberland WG, Yang OO (2013) A stochastic multi-scale model of HIV-1 transmission for decision-making: application to a MSM population. PloS one 8. [PMC free article] [PubMed]
51. Rosinska M, Simmons R, Marzec-Boguslawska A, Janiec J, Porter K (2015) Relating HIV testing patterns in Poland to risky and protective behaviour. AIDS care: 1–9. [PubMed]
52. Hull MW, Wu Z, Montaner JS (2012) Optimizing the engagement of care cascade: a critical step to maximize the impact of HIV treatment as prevention. Current opinion in HIV and AIDS 7: 579–586. doi: 10.1097/COH.0b013e3283590617 [PubMed]
53. Pennings PS, Holmes SP, Shafer RW (2014) HIV-1 transmission networks in a small world. The Journal of infectious diseases 209: 180–182. doi: 10.1093/infdis/jit525 [PMC free article] [PubMed]
54. Smith DM, May SJ, Tweeten S, Drumright L, Pacold ME, Kosakovsky Pond SL, et al. (2009) A public health model for the molecular surveillance of HIV transmission in San Diego, California. AIDS 23: 225–232. doi: 10.1097/QAD.0b013e32831d2a81 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science