Evolving patterns of HIV-1 transmitted drug resistance in Poland in the years 2000-2008

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

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.

Molecular epidemiology of HIV-1 infection in Europe: An overview

Human Immunodeficiency Virus type 1 (HIV-1) is characterised by vast genetic diversity. Globally circulating HIV-1 viruses are classified into distinct phylogenetic strains (subtypes, sub-subtypes) and several recombinant forms. Here we describe the characteristics and evolution of European HIV-1 epidemic over time through a review of published literature and updated queries of existing HIV-1 sequence databases. HIV-1 in Western and Central Europe was introduced in the early-1980s in the form of subtype B, which is still the predominant clade. However, in Eastern Europe (Former Soviet Union (FSU) countries and Russia) the predominant strain, introduced into Ukraine in the mid-1990s, is subtype A (A_FSU) with transmission mostly occurring in People Who Inject Drugs (PWID). In recent years, the epidemic is evolving towards a complex tapestry with an increase in the prevalence of non-B subtypes and recombinants in Western and Central Europe. Non-B epidemics are mainly associated with immigrants, heterosexuals and females but more recently, non-B clades have also spread amongst groups where non-B strains were previously absent - non-immigrant European populations and amongst men having sex with men (MSM). In some countries, non-B clades have spread amongst the native population, for example subtype G in Portugal and subtype A in Greece, Albania and Cyprus. Romania provides a unique case where sub-subtype F1 has predominated throughout the epidemic. In contrast, HIV-1 epidemic in FSU countries remains more homogeneous with A_FSU clade predominating in all countries. The differences between the evolution of the Western epidemic and the Eastern epidemic may be attributable to differences in transmission risk behaviours, lifestyle and the patterns of human mobility. The study of HIV-1 epidemic diversity provides a useful tool by which we can understand the history of the pandemic in addition to allowing us to monitor the spread and growth of the epidemic over time.

Distribution and time trends of HIV-1 variants in Poland: Characteristics of non-B clades and recombinant viruses

The spread of HIV-1 subtypes varies considerably both worldwide and within Europe, with non-B variants commonly found across various exposure groups. This study aimed to analyse the distribution and temporal trends in HIV-1 subtype variability across Poland. For analysis of the subtype distribution, 1219 partial pol sequences obtained from patients followed up in 9 of 17 Polish HIV treatment centres were used. Subtyping was inferred using the maximum likelihood method; recombination was assessed using the bootscanning and jumping profile hidden Markov model methods. Subtype B dominated in the studied group (n=1059, 86.9%); in 160 (13.1%) sequences, non-B variants were present [A1 (n=63, 5.2%), D (n=43, 3.5%), C (n=22, 1.8%), and F1 (n=2, 0.2%)]. In 25 (2.1%) cases circulating recombinant forms (CRFs) were found. Five A1 variants (0.4%) were unique AB recombinant forms (URF) not previously identified in Poland. Non-B clades were notably more common among females (n=73, 45.6%, p<0.001) and heterosexual individuals (n=103, 66.5%, p<0.001) and less frequent among men who have sex with men (MSM) (n=27, 17.42%, p<0.001). HIV-1 viral load at diagnosis was higher among non-B cases [median: 5.0 (IQR: 4.4-5.6)] vs. [median: 4.8 (IQR: 4.3-5.4) log copies/ml for subtype B (p<0.001)] with a lower CD4(+) lymphocyte count at baseline [median: 248 (IQR: 75-503) for non-B vs. median: 320 (IQR: 125-497) cells/μl for subtype B; p<0.001]. The frequency of the non-B subtypes proved stable from 2008 (11.5%) to 2014 (8.0%) [OR: 0.95 (95% CI: 0.84-1.07), p=0.4], with no temporal differences for exposure groups, gender, age and AIDS. Despite the predominance of subtype B, the variability of HIV in Poland is notable; both CRFs and URFs are present in the analysed population. Non-B variants are associated with heterosexual transmission, more advanced HIV disease and have stable temporal frequencies.

Transmitted HIV drug resistance in antiretroviral-treatment-naive patients from Poland differs by transmission category and subtype

OBJECTIVES

The surveillance of HIV-transmitted drug resistance mutations (t-DRMs), including temporal trends across subtypes and exposure groups, remains a priority in the current management of the epidemic worldwide.

METHODS

A cross-sectional analysis of 833 treatment-naive patients from 9 of 17 Polish HIV treatment centres. Partial pol sequences were used to analyse drug resistance with a general time reversible (GTR)-based maximum likelihood algorithm used for cluster/pair identification. Mutation frequencies and temporal trends were investigated.

RESULTS

t-DRMs were observed in 9% of cases (5.8% for NRTI, 1.2% NNRTI and 2.0% PI mutations) and were more common among heterosexually infected (HET) individuals (13.4%) compared with MSM (8.3%, P = 0.03) or injection drug users (IDUs; 2.9%, P = 0.001) and in MSM compared with IDUs (P = 0.046). t-DRMs were more frequent in cases infected with the non-B variant (21.6%) compared with subtype B (6.6%, P < 0.001). With subtype B a higher mutation frequency was found in MSM compared with non-MSM cases (8.3% versus 1.8% for IDU + HET, P = 0.038), while non-B variants were associated with heterosexual exposure (30.4% for HET versus 4.8% for MSM, P = 0.019; versus 0 for IDU, P = 0.016). Trends in t-DRM frequencies were stable over time except for a decrease in NNRTI t-DRMs among MSM (P = 0.0662) and an NRTI t-DRM decrease in HET individuals (P = 0.077). With subtype B a higher frequency of sequence pairs/clusters in MSM (50.4%) was found compared with HET (P < 0.001) and IDUs (P = 0.015).

CONCLUSIONS

Despite stable trends over time, patterns of t-DRMs differed notably between transmission categories and subtypes: subtype B was associated with MSM transmission and clustering while in non-B clades t-DRMs were more common and were associated with heterosexual infections.

Burden of nonnucleoside reverse transcriptase inhibitor resistance in HIV-1-infected patients: a systematic review and meta-analysis

The prevalence of HIV drug resistance varies with geographic location, year, and treatment exposure. This study generated yearly estimates of nonnucleoside reverse transcriptase inhibitor (NNRTI) resistance in treatment-naive (TN) and treatment-experienced (TE) patients in the United States (US), Europe (EU), and Canada. Studies reporting NNRTI resistance identified in electronic databases and 11 conferences were analyzed in three groups: (1) TN patients in one of four geographic regions [US, Canada, EU countries with larger surveillance networks ("EU1"), and EU countries with fewer data ("EU2")]; (2) TE patients from any region; and (3) TN patients failing NNRTI-based treatments in clinical trials. Analysis data included 158 unique studies from 22 countries representing 84 cohorts of TN patients, 21 cohorts of TE patients, and 8 trials reporting resistance at failure. From 1995 to 2000, resistance prevalence in TN patients increased in US and EU1 from 3.1% to 7.5% and 0.8% to 3.6%, respectively. Resistance in both regions stabilized in 2006 onward. Little resistance was identified in EU2 before 2000, and increased from 2006 (5.0%) to 2010 (13.7%). One TN Canadian study was identified and reported resistance of 8.1% in 2006. Half of TN clinical trial patients had resistance after treatment failure at weeks 48-144. Resistance in TE patients increased from 1998 (10.1%) to 2001 (44.0%), then decreased after 2004. Trends in NNRTI resistance among TN patients show an increased burden in the US and some EU countries compared to others. These findings signify a need for alternate first-line treatments in some regions.

Transmitted drug resistance to rilpivirine among antiretroviral-naïve patients living with HIV from northern Poland

Introduction

Rilpivirine (RPV) is a second-generation non-nucleoside reverse transcriptase inhibitor (NNRTI) that was recently approved for the treatment of antiretroviral-naïve individuals with HIV-1 viral load of <100,000 copies/ml. As transmission of the drug resistance mutations to this NNRTI may affect treatment outcomes, the frequency of primary, RPV-associated drug resistance mutations was assessed in this study.

Methods

For the study, 244 viral genome sequences from antiretroviral-naïve individuals were obtained by bulk sequencing. RPV-associated mutations were divided into RPV resistance mutations (K101E/P, E138A/G/K/Q/R, V179L, Y181C/I/V, Y188L, H221Y, F227C and M230I/L) according to the International AIDS Society-USA (IAS-USA) mutation list and variants potentially affecting RPV susceptibility (L100I, K101H/T, E138S, V179F/D/G/T, G190A/E/S, F227L and M230V) based on the in vitro and in vivo data.

Results

IAS-USA RPV drug resistance mutations were found in 5.3% sequences, with E138A and E138G being the most common (3.7 and 0.8%, respectively), followed by K101E (0.4%) and Y181C (0.4%), with no significant differences in the frequency between subtype B and non-B clades. Mutations potentially reducing RPV susceptibility were found in 2.5% of sequences, and they included V179D (1.6%) and G190A (0.8%), with equal distribution among non-B (n=2, 2.5%) and subtype B (n=4, 2.5%) clades. Clustering of RPV mutations was infrequent.

Conclusions

Prevalence of RPV-associated drug resistance mutations was low in the analysed sample and did not vary across the subtypes. The frequency of variants with potential influence on RPV susceptibility was similar among non-B variants if compared to B clades. Transmitted drug resistance to RPV is uncommon, which makes this a good option for the treatment of ARV-naïve patients; however, genotype resistance testing should remain compulsory before starting an RPV-based regimen.

The emerging and re-emerging human immunodeficiency virus epidemics in Europe

We provide a summary of the current status of the human immunodeficiency virus (HIV) epidemic, including data on the transmission of drug-resistant virus in the European Region of the WHO. The review was conducted by searching the reports of the European Centre for Disease Prevention and Control and the United Nations General Assembly Special Session country reports to identify the number of HIV cases reported in 2002-2011, the number of HIV tests performed, and the results of the most recent HIV surveys in at-risk groups. In the West, a 5-year linear trend analysis suggests an increase in the number of reported HIV cases in Belgium, Greece, and Iceland, and a decline in The Netherlands, Ireland, Portugal, Switzerland, and the UK. In the Centre, the number of reported cases increased in Bulgaria and Turkey, and in the East, the number of reported cases increased in Armenia, Georgia, and Ukraine. Estonia and Latvia reported a significant downward trend. HIV transmission in men who have sex with men (MSM) and heterosexual transmission accounted for, respectively, 40% and 38% of newly reported HIV cases in the West in 2011, whereas the epidemic in the Centre is predominantly concentrated in MSM. Although injection drug use remains the key driving force of the epidemic in the East, there is increasing evidence of sexual transmission from the core groups of injection drug users, and increasing spread within MSM. The pattern of transmitted drug resistance (TDR) is heterogeneous accross Europe; a significant number of clusters of TDR have been recently observed in several European countries.

Molecular epidemiology of recent HIV-1 infections in southern Poland

The genetic diversity of human immunodeficiency virus type 1 (HIV-1) offers an opportunity to track the development of the epidemic across different populations. Viral pol gene fragments from 55 individuals of Polish origin with recent HIV-1 infection identified in 2008-2010 in four Polish cities were analyzed. Viral sequences were compared with sequences from 100 individuals (reference group) infected before 2004. Viral spread among groups with different HIV transmission categories was compared using a phylogenetic approach. The majority of sequences from individuals with recent infection were subtype B (93%) within which four transmission clusters (18% of samples) were detected. Samples from men infected through sex between men and from persons infected through injecting drugs were broadly separated (P < 0.0001), while samples from individuals infected by heterosexual contacts were dispersed uniformly within phylogenetic tree (P = 0.244) inferred from viral sequences derived from individuals infected recently and the reference group. The percentage of samples from persons infected by heterosexual contacts which clustered with samples from men infected through sex between men was not significantly higher for those with recent infection (47%), compared to the reference group (36%). In conclusion, men infected by sex between men and individuals infected through injecting drugs appear to form separate HIV transmission networks in Poland. The recent spread of HIV-1 among persons infected with subtype B by heterosexual contacts appears to be linked to both these groups.

The prevalence of transmitted drug resistance in newly diagnosed HIV-infected individuals in Croatia: the role of transmission clusters of men who have sex with men carrying the T215S surveillance drug resistance mutation

The aim of this study was to determine the prevalence of transmitted drug resistance (TDR) in newly diagnosed and treatment-naive HIV-infected patients from Croatia and evaluate a possible contribution of transmission clusters to the spread of resistant virus. The study enrolled treatment-naive HIV-infected patients that entered clinical care at the Croatian Reference Center for HIV/AIDS between 2006 and 2008. The protease gene and a part of the reverse transcriptase gene of the HIV-1 genome were sequenced by using the Trugene HIV-1 Genotyping System. The prevalence of transmitted drug resistance was analyzed by using the surveillance drug resistance mutations (SDRM) list recommended by the WHO in 2009. We report findings for 118 of 180 eligible patients (65.6% coverage). SDRM were detected in 26 of 118 patients (22.0%) who were infected with subtype B and belonged mostly to the men having sex with men (MSM). The majority of patients with primary resistance carried SDRM associated with resistance to nucleoside analogues reverse transcriptase inhibitors (NRTIs, 23 of 118 patients, 19.5%). The most frequently found NRTI SDRM was T215S (17 of 118 patients, 14.4%). SDRM associated with resistance to nonnucleoside reverse transcriptase inhibitors were detected in three (2.5%) patients and primary resistance to protease inhibitors was not detected. Non-B subtypes were detected in 13/118 patients (11%). A total of 12 transmission pairs and eight distinct transmission clusters were identified with the largest cluster harboring sequences from 19 patients; among them all but two were carrying the T215S mutation. This study showed a high prevalence of TDR in newly diagnosed MSM from Croatia and is an important contribution concerning the relationship between local transmission clusters and the spread of resistant virus.

HIV-1 subtype D infections among Caucasians from Northwestern Poland--phylogenetic and clinical analysis

BACKGROUND

HIV-1 subtype D infections, which are associated with a faster rate of progression and lymphocyte CD4 decline, cognitive deficit and higher mortality, have rarely been found in native Europeans. In Northwestern Poland, however, infections with this subtype had been identified. This study aimed to analyze the sequence and clinical data for patients with subtype D using molecular phylogeography and identify transmission clusters and ancestry, as well as drug resistance, baseline HIV tropism and antiretroviral treatment efficacy.

METHODS

Phylogenetic analyses of local HIV-1 subtype D sequences were performed, with time to the most recent common ancestor inferred using bayesian modeling. Sequence and drug resistance data were linked with the clinical and epidemiological information.

RESULTS

Subtype D was found in 24 non-immigrant Caucasian, heterosexually infected patients (75% of females, median age at diagnosis of 49.5 years; IQR: 29-56 years). Partial pol sequences clustered monophyletically with the clades of Ugandan origin and no evidence of transmission from other European countries was found. Time to the most common recent ancestor was 1989.24 (95% HPD: 1968.83-1994.46). Baseline drug resistance to nucleoside reverse transcriptase inhibitors was observed in 54.5% of cases (mutations: M41L, K103N, T215S/D) with evidence of clustering, no baseline integrase or protease resistance and infrequent non-R5 tropism (13.6%). Virologic failure was observed in 60% of cases and was associated with poor adherence (p<0.001) and subsequent development of drug resistance (p = 0.008, OR: 20 (95%CI: 1.7-290).

CONCLUSIONS

Local subtype D represented an independently transmitted network with probably single index case, high frequency of primary drug resistance and evidence of transmission clusters.