HIV-1 Subtype D Infections among Caucasians from Northwestern Poland—Phylogenetic and Clinical Analysis

Abstract

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.

About the Authors:

Miłosz Parczewski

* E-mail: [email protected]

Affiliation: Department of Infectious Diseases and Hepatology, Pomeranian Medical University, Szczecin, Poland

Magdalena Leszczyszyn-Pynka

Affiliation: Department of Infectious Diseases and Hepatology, Pomeranian Medical University, Szczecin, Poland

Dorota Bander

Affiliation: Department of Infectious Diseases and Hepatology, Pomeranian Medical University, Szczecin, Poland

Anna Urbanska

Affiliation: Department of Infectious Diseases and Hepatology, Pomeranian Medical University, Szczecin, Poland

Anna Boroń-Kaczmarska

Affiliation: Department of Infectious Diseases and Hepatology, Pomeranian Medical University, Szczecin, Poland

Introduction

Phylogenetic analyses have been widely implemented in molecular epidemiology of HIV infection and allowed not only to understand genetic sequence variability, but also to trace the flow of the variants and drug resistant clades as well as to infer a rate of molecular evolution to date the ancestry [1]. Clusters of sexual transmissions have been reconstructed and dated using sequence data obtained for drug resistance genotyping, providing information on the transmission patterns and allowing to identify distinct epidemiological entities [2]. Tracking changes in HIV epidemics related to subtype diversity and drug resistance is important in Europe, where genetic diversity of HIV is rising, with a trend to increase the prevalence of non-subtype B infections [3]. This variation is related either to immigration or the founder effect with subsequent rapid spread within the populations at risk. Migration from Africa and South America has been identified as a major source of diversity in Western Europe [4] with multiple, continuous inflow of the new clades into the population and common association with heterosexual transmissions [5], [6], [7]. Onward infections with the non-subtype B clades reported in Europe were mostly acquired from the partners originating from countries of the higher HIV prevalence [8], [9], however unique epidemics among drug users have also been noted [10]. Such explosive outbreaks in countries where drug injection is, or used to be common, were observed both for subtype B and non-B clades. This often results in high prevalence of a variant uncommon elsewhere, as was the case in Estonia and Kaliningrad region of Russia [11], [12]. Studies have also revealed a nosocomial epidemics with subtype F, of the Angolan ancestry, among children in Romania due to a distribution of contaminated needles and blood products [13], [14].

Acquisition of non-B subtypes may be associated with high prevalence of minor drug resistance mutations in untreated individuals and changes in virus tropism [15], [16], [17], [18]. Prevalence of major protease and reverse transcriptase mutations is similar across the variants, however both virologic and immunologic treatment efficacy is similar to the subtype B [19], [20]. In general, frequency of primary drug resistance differs locally in Europe, often being associated with transmission clusters, fluctuates over time and is usually limited to a single class of antiretrovirals [21], [22], [23], [24]. In Poland the prevalence of transmitted drug resistance in the years 2006–2009 ranged from 3.9% to 7.4% [25]. The predominant subtype B was observed in 95.8% of the newly diagnosed HIV infections [25].Of note, sexual transmissions now tend to be the most common in Poland, which is a significant change of epidemiological pattern from the earlier dominant injection drug use related transmissions [26]. We have previously identified higher prevalence of non-B subtypes, with a number of subtype D clades in the local population from Northwestern Poland (21%). These variants have been associated with more advanced HIV disease and more common heterosexual transmission route if compared to the subtype B [27].

In this study the sequences and clinical data from a cluster of locally diagnosed individuals infected with subtype D were analyzed. Firstly, we used molecular phylogeography to trace the import of subtype D into Poland, infer between-host phylogenies, indentify transmission clusters and determine ancestry with molecular clock methods. Secondly, baseline drug susceptibility (protease, reverse transcriptase and integrase) and V3 based predicted viral tropism were assessed. Finally, treatment efficacy among patients infected with this variant and secondary drug resistance were summarized.

Results

Group characteristics

Of the total of 190 patient sequences sampled so far in at the Department of Infectious Diseases and Hepatology, Szczecin, Poland HIV-1 subtype D was found in 24 (12.63%) cases. The first case of the infection with this variant was confirmed in a year 1994 in a young female, who had never travelled outside Poland. All cases were of Caucasian origin with the predominance of females (n = 18, 75%) and a median age at diagnosis of 49.5 years (IQR: 29–56 years). For only one patient presumed infection from a partner travelling overseas (sailor) was recorded, however the source sample was not available. Travel history for this contact is not specific, with frequent worldwide travels and no clear African connection. In the remaining cases no history of travel outside Europe was noted.

For all cases heterosexual transmissions were recorded, and patients were HCV antibody negative at HIV diagnosis. In one case HCV infection with antibody production was observed in the third year of observation. Patients tended to be diagnosed at symptomatic stages of the HIV infection (n = 21, 87.5%) with a low baseline and nadir lymphocyte CD4 counts (median 62 cells/µl [IQR: 23–246 cells/µl] and 35 cells/µl [IQR: 10–144 cells/µl, respectively) and a high baseline HIV-1 viraemia (median of 5.6 log HIV-1 copies/ml [IQR: 4.9–6 log HIV-1 copies/ml]).

Total observation period in the group was 116 person-years with a mean observation time of 57.82 months (95% CI: 35.37–80.27). In fourteen cases (58.3%) AIDS was diagnosed within the period of observation. Mortality rate was 3.44 per 100 person-years. Of the four deaths recorded, three were AIDS-related and occurred within the first year of observation; two cases of large B-cell lymphoma, one Pneumocystis jirovecii pneumonia and the fourth related to the cardiac event eight years after HIV diagnosis.

Phylogeography and dated phylogeny of local subtype D

For the phylogeographic and dated phylogeny the earliest possible sample from every patient was selected and sequenced. The majority of samples were obtained from treatment-naive patients (n = 22, 91.7%), the remaining two samples were from patients on ARV.

In the maximum likelihood tree all sequences from Poland formed a monophyletic cluster nested within the Ugandan sequences (figure 1, shown in red). Of note, subtype D sequences found in various European countries (figure 1, marked in yellow) were interspersed among various lineages of African origin, with no clustering with Polish sequences. Branch lengths of our sequences indicated independent evolution of the imported variant without admixture of African variants or transmission to other European countries.

[Figure omitted. See PDF.]

Figure 1. Relationships between the subtype D sequences inferred using maximum likelihood phylogeny (GenBank deposited pol sequences).

Country-specific sequences are marked with the same color: red – Poland, blue – Uganda, green – Tanzania, yellow – Europe (except Poland), brown – Cameroon, magenta - Senegal, cyan- Sudan, dark green – other African countries, violet – South America, grey – Asia, orange – North America.

https://doi.org/10.1371/journal.pone.0031674.g001

In the reconstruction of the subtype D phylogeny, the maximum likelihood tree (figure 2a) and the Bayesian clade credibility tree (figure 2b) with dated ancestry had a similar topology. Time to the most recent common ancestor inferred using the relaxed clock model with GTR+γ+Γ was in accordance to the transmission dates from the patients' documentation. Four transmission sub-clusters may be observed with both maximum likelihood bootstrap values and Bayesian posterior probability >90% (figure 2b, sub-clusters numbered 1–4, blue boxes). Using Bayesian Monte Carlo Markov Chain, dating to the time to the most recent common ancestor (tMRCA) for all regional sequences was the year 1989 (table 1). By the year 1994, which was the date of the first confirmed infection with the subtype D, the virus was probably already circulating locally. Bayesian skyline plot (figure 3) indicated that a population size has been fairly stable over time with a plateau in the last three years.

[Figure omitted. See PDF.]

Figure 2. Phylogenetic trees of the subtype D sequences from Northwestern Poland.

Figure a - maximum likelihood tree with bootstrap values for 1000 replicates drawn at the branches. Figure b – time scaled Bayesian MCMC tree. On the tree branches estimated time to the most recent common ancestor (tMRCA) and posterior probabilities expressed as percentage are shown. For both figures clustered sequences are marked in red and four identified clusters indicated as blue boxes and numbered are drawn on the right. Drug resistance mutations are marked at the tip nodes after the sequence identifier. *source patient for the transmission of the drug resistance within the cluster.

https://doi.org/10.1371/journal.pone.0031674.g002

[Figure omitted. See PDF.]

Figure 3. Bayesian skyline plot for estimation of the number of subtype D HIV-1 cases in the local population.

95% CI are marked in blue. Y-axis: predicted number of cases (log scale), X-axis: timescale (years).

https://doi.org/10.1371/journal.pone.0031674.g003

[Figure omitted. See PDF.]

Table 1. Time to the most recent common ancestors (tMRCA) for the entire cluster and sub-clustered sequences.

https://doi.org/10.1371/journal.pone.0031674.t001

Clusters of drug resistance

Four observed clustered sequences were linked epidemiologically, which was reflected in the patients' records – cases within the sub-clusters 1–3 were sexual partners. Baseline drug resistance was observed in 12 of the 22 antiretroviral naive patients (54.5%) with pretreatment genotyping available. Primary drug resistance was the most common within the fourth, the youngest (tMRCA in the 2002 [95% HPD: 2005.36–1998.52]) sub-cluster (10/11 cases, 90.9%), with nucleoside reverse transcriptase (NRTI) mutations observed. Ten patients within this sub-cluster were female and one male; epidemiological data from patients' records suggested possible sexual contacts between the man and at least three women. Median age in this subgroup was 56 years (IQR: 53–60 years); all patients were diagnosed with clinical features of immunodeficiency (category B and C according to the CDC 1993).

In two patients non-nucleoside reverse transcriptase (NNRTI) drug resistant virus was observed in pretreatment samples (as shown in figure 2b). It was possible to identify the source of infection for one of these cases (sub-cluster 3, source marked with asterisk, data from the sequencing of the sample taken after virological failure, sequence not included in the tree). Frequency of the baseline drug resistance has been shown in figure 4. No baseline integrase or protease resistance was observed in the entire group. Non-R5 tropism was assigned in three (13.6%) baseline, treatment-naive samples with geno2pheno false positive rate (FPR) threshold set to 10% (one if FPR of 5.75% was used), sequences with FPR<10% were not clustered.

[Figure omitted. See PDF.]

Figure 4. Frequency of baseline drug resistance mutations among treatment-naive patients.

https://doi.org/10.1371/journal.pone.0031674.g004

Response to antiretroviral treatment and secondary drug resistance

The group of 20 (83.3%) patients were commenced on stable antiretroviral therapy. Of the remaining one patient refused consistently, one was lost to follow-up within a month following the HIV diagnosis, one died before the treatment was introduced and one died within the first month of treatment. The median for the highest recorded lymphocyte CD4 count on cART was 369.5 (IQR: 294.25–596.5) cells/µl, with only five (25%) patients maintaining stable CD4 count >500 cells/µl for at least six months during the follow-up period. In 10 (50%) cases clinician recorded poor adherence, while pharmacy records confirmed <90% adherence in 9 (45%) patients. In 12 (60%) cases the primary therapy was virologically successful (table 2), and as expected, failure was associated with poor or <90% adherence (p<0.001) and development of drug resistance (p = 0.008, OR: 20 (95%CI: 1.7–290). Secondary drug resistance was noted in all eight failing patients with NRTI (M184V, T215 complex) and NNRTI (K103N) being the most common mutations (figure 5). No secondary protease drug resistance mutations were found in the group. For virologically failing patients non-R5 tropism was assigned in two of these cases (25%, geno2pheno algorithm with FPR 10%); if FPR 5.75% was used the non-R5 tropism would be called in one (12.5%) patient.

[Figure omitted. See PDF.]

Figure 5. Frequency of secondary drug resistance mutations in patients failing antiretroviral therapy.

https://doi.org/10.1371/journal.pone.0031674.g005

[Figure omitted. See PDF.]

Table 2. Drug resistance and treatment efficacy in the group infected with subtype D.

https://doi.org/10.1371/journal.pone.0031674.t002

Interestingly, the treated patients with the baseline drug resistance tended to present with higher rates of virologic success if compared to the ones with drug susceptible viruses (90%, versus 30%, p = 0.02997, 2-sided Fisher's exact test), however this phenomenon was associated with higher adherence.

Discussion

In this study we have evaluated the subtype D infections occurring among Caucasian individuals from Northwestern Poland. Differently to the cases of infection with this subtype observed across Europe [15], [17], [28], [29], the variant was found among native, non-immigrant and non-travelling Poles, and acquired locally. As is common in our country, HIV was diagnosed late, usually with symptoms of immunodeficiency and low baseline lymphocyte CD4 count [30]. Moreover, subtype D diagnosis in the advanced stages of the disease may be related to a faster clinical progression observed among patients infected with this subtype, lower lymphocyte CD4 count during follow-up and faster rates of lymphocyte CD4 decline, as previously described [31], [32], [33]. This may be the result of the increased degree of apoptosis in lymphocyte CD4 populations prior to the start of the antiretroviral treatment and a higher number of inhibitory programmed cell death receptors on the CD4 cells [34]. Mortality in the group was slightly higher if compared to the entire cohort (majority infected with the subtype B) from the region (3.44 versus 2.66 deaths per 100 person-years), which is in accordance to the previous reports [25], [35], [36], [37], [38]. Of note, all infections were associated with heterosexual transmissions with a percentage of infected women and mean age at diagnosis higher than among typical HIV positive population observed in Poland (75% versus 24.4%, respectively [26]). This may partially explain late diagnosis and as a consequence higher mortality and lower lymphocyte CD4 counts as heterosexual females, especially in peri- and postmenopausal period, are infrequently tested for HIV.

Observed infections with subtype D infections formed a network with probably one index case and subsequent local spread. This cluster of heterosexually acquired infections seemed independent from subtype B ones spreading among men who have sex with men (MSM),injection drug users (IDU) and heterosexual populations [27]. Dated phylogeny and phylogeographic analyses suggested that the tMRCA of the subtype D found in Poland is 1989.2 (95% HPD: 1994.46–1968.83) and is of Ugandan origin. The sequences cluster monophyletically within the sequences from Uganda, which confirmed a single transmission event in the past with the subsequent local spread. Local character of the transmissions was further confirmed by a lack of clustering with sequences from other European countries. Networks of transmission had been investigated both for B and non-B subtype infections using Bayesian MCMC inference methods, with identification of transmission clusters among MSM [2], [39] and heterosexual individuals [40]. Despite evidence of the slower dynamics among heterosexually acquired non-B subtypes, both small and large clusters with epidemiologically linked patients were found in the studies cited above. Similarly, in our study three small (two sequences) and one large (eleven) clusters were identified. Moreover, the Bayesian skyline plot indicates a stable population size. Dated phylogeny for the large cluster revealed that the majority of these transmissions occurred within the short period of time, probably in the early stage of infection.

Additionally, the high frequency of primary drug resistance mutations (DRMs) observed in this study was associated to the observed clustering within the network of transmission and presence of M41L and T215 revertant mutations. Similar results were reported in various countries for different subtypes, with most common transmission of NRTI DRMs [22], [41], [42], [43]. As expected, no integrase inhibitor (InI) transmitted resistance was found, which is consistent to the fact that all infections were acquired before introduction of InI in the clinical practice. We noted, that transmission of the drug resistance occurred both among individuals newly diagnosed, treatment-naive patients (cluster 4) and from a chronically infected, heavily treated patient (cluster 3). Both phenomena have been described previously [43], [44], [45], [46], [47], however this has been the first report on clustered drug resistance in Caucasians with subtype D. Non-R5 viruses were found infrequently in the group (in 13.6% when FPR 10% was used for the tropism assignment and 4.5% for 5.75 FPR), similarly to the result (14.8%) from a phenotypic assay performed in French patients [48]. Of note, our finding is in contrast to the studies demonstrating a high frequency of an ×4 tropism in Ugandan clades [49].

Immunologic and virologic antiretroviral therapy efficacy in the observed, stably treated group was poorer than expected for individuals on cART with no specific drug combination clearly superior in the group. Virologic failure was noted in 40% of cases, compared to the 24.7% in the European observational studies [50]. Association between poor adherence, treatment failure and development of secondary drug resistance was obvious, however the fact that the primary drug resistance did not affect the virologic treatment efficacy in the setting of good adherence must be acknowledged. It may be hypothesized that the poor adherence may be associated with the higher probability of cognitive impairment observed among individuals infected with this subtype [51]. Recently a pathogenetic mechanism of this phenomenon was proposed, as a certain features of the subtype D nef protein 3D structure was associated with long-term progressive dementia in HIV patients, suggesting that this cognitive deficit may be related to the altered folding or binding potential of this protein [52].

It has to be noted, that the limitations of this study can be related to the fact that the certain samples with subtype D might have remained undetected, as it was not possible to sequence all the samples from the patients followed-up in the centre. Neither index case, nor the certain source patient for the cluster 4 were found.

To sum up, local spread of HIV-1 subtype D infections described here represented an independent cluster, detected among heterosexually exposed individuals in parallel to the widely distributed subtype B infections. In phylogenetic inference analyses Ugandan origin of the virus was identified with probable single transmission event. Frequency of drug resistance in this group was high, especially within sub-clustered sequences, while the treatment efficacy was poor and independent of the primary drug resistance. Late diagnoses and association with heterosexual mode of transmission might fuel spread of these infections, especially in the light of evidence for the epidemiological clustering which suggests that new transmissions of HIV subtype D have been under diagnosed.

Methods

Population and ethics statement

For the study participants followed up at the Department of Infectious Diseases and Hepatology, Szczecin, Poland were recruited. The centre provides care for HIV infected patients, mostly from the Northwest part of Poland from the early 90-ties, and so far 960 patients have been observed. Clinical, epidemiological and laboratory data have been collected as a part of a routine patient care. The study has been approved by the local ethical committee of Pomeranian Medical University, Szczecin, Poland. All patients for whom sequence data (n = 190) were available have been analyzed. Treatment adherence was assessed based on the patient records (number of months of medications dispensed by the number of months of follow-up, expressed as a percentage) as well as clinical assessment based on the patient's statement regarding missed doses at home or treatment interruptions. Antiretroviral therapy was commenced according to the clinical standard at the time of introduction, stable therapy was defined as >30 consecutive days of uninterrupted treatment. Virologic success of the therapy was defined as stable serum viral load <400 copies/ml in the first available analysis following 6 months of treatment, and as a consequence no evidence of secondary drug resistance.

Genotyping and subtyping

Plasma samples collected from patients between the years 1996 and 2011 were stored at −80 degrees C and analyzed retrospectively. HIV RNA extraction as well as reverse transcriptase and protease genotyping was performed using Viroseq 2.7 and 2.8 genotyping assays (Abbott Molecular, Abbott Park, IL) according to manufacturer's protocol. Initial subtyping was performed using REGA genotyping 2.0 tool (http://bioafrica.mrc.ac.za/rega-genotype/html/subtypinghiv.html) based on the partial pol sequence obtained by Viroseq methodology. Obtained protease/reverse transcriptase sequences were 1302 b.p. long; location from the start of HXB2 genome: positions 2253–3525. Subsequently, the subtype was confirmed in phylogenetic analysis with a selection of reference sequences listed in the HIV Sequence Compendium 2010 (Los Alamos National Laboratory, Los Alamos, U.S.A. http://www.hiv.lanl.gov). For this purpose bootstrapped (1000 replicates) neighbor-joining tree with Kimura 2 parameter model was used (MEGA version 4.0 [53]). To exclude recombinant sequences, bootscanning (200 base pair (b.p.) window, step 20 b.p., Kimura 2 parameter model) using Simplot software (S. Ray, John Hopkins University, Baltimore, USA) with a set of reference sequences was performed [54]. When necessary, samples from treatment-failing subtype D patients (genotyped for the purpose of drug resistance monitoring) were resequenced from the baseline (first available) sample. V3 sequences were obtained from both baseline and on-treatment samples while integrase was sequenced from the first samples only. HIV-1 integrase was amplified and sequenced using reagents and conditions specified by Laethem et al. [55], with the following amplification and sequencing primers: AGGAGCAGAAACTTWCTATGTAGATGG (outer forward), TTCTTCCTGCCATAGGARATGCCTAAG (outer reverse), TTCRGGATYAGAAGTAAAYATAGTAACAG (inner forward), TCCTGTATGCARACCCCAATATG (inner reverse and sequencing reverse), GCACAYAAAGGRATTGGAGGAAATGAAC (sequencing, forward), GGVATTCCCTACAATCCCCAAAG (sequencing, forward), GAATACTGCCATTTGTACTGCTG (sequencing, reverse). Obtained integrase sequences were 866 b.p. long, location from the start of HXB2 genome: positions 4230–5096. V3 loop sequencing was carried out following the protocol provided by the HIV Centre of Excellence (personal communication, prof. Richard Harrigan), with the following amplification and sequencing primers: GAGCCAATTCCCATACATTATTGT (outer forward), GCCCATAGTGCTTCCTGCTGCTCCCAAGAACC (outer reverse), TGTGCCCCAGCTGGTTTTGCGAT (inner forward), TATAATTCACTTCTCCAATTGTCC (inner reverse), AATGTCAGYACAGTACAATGTACAC (sequencing, forward), GAAAAATTCCCTTCCACAATTAAA (sequencing, reverse). Obtained integrase sequences were 107 b.p. long, location from the start of HXB2 genome: positions 7110–7217. Briefly, nested PCR was carried out following reverse transcription of extracted HIV-1 RNA as stated in the methodologies cited above. Amplicons were used for sequencing by standard techniques with BigDye technology using an ABI 3500 platform (Applied Biosystems, Foster City, CA, USA). V3 loop analyses were done in triplicate. Two overlapping sequencing reactions (forward and reverse) were performed for each sample. Reverse transcriptase and protease sequences were assembled with a software provided with a Viroseq 2.8 kit, while V3 and integrase sequences using the Recall online tool (http://pssm.cfenet.ubc.ca) [56]. Integrase and V3 sequences were obtained for all patients. Obtained V3 sequence data were interpreted using the geno2pheno tool (www.geno2pheno.org) with two FPR thresholds: 10% (as suggested by the European Guidelines on HIV Tropism testing) and 5.75% (from the Merit trial) [57], [58], [59].Drug resistance interpretation was performed using the Stanford DB database (hivdb.stanford.edu) [60] with transmitted HIV-1 drug resistance interpreted according to the surveillance recommendations [61], [62].

Phylogenetic inference

To analyze the origin of our sequences, we have downloaded all HIV-1 group M subtype D pol sequences of at least 900 b.p. length, stored in the GenBank, using Los Alamos National Laboratory HIV sequence database (www.hiv.lanl.gov accessed on 16.09.2011). This dataset of 741 sequences was aligned with our partial pol sequences using Clustal ×2.0.10 (www.clustal.org). Bootstrapped (1000 replicates) maximum likelihood phylogeny using a GTR+γ model with empirical nucleotide frequencies was inferred using PAUP* v. 4.0 (http://paup.csit.fsu.edu/order.html). Separate analyses were performed for the dataset of 24 local sequences with maximum-likelihood methodology (bootstrapped with 1000 replicates) using PHYML v. 3.0 software with subtree pruning and regrafting (SPR) algorithm under GTR substitution model [63]. Moreover, to estimate the time to the most recent common ancestor (tMRCA) for local sequences we used a Bayesian Monte Carlo Marcov Chain (MCMC). Two replicates of 10 million generations were run in BEAST v. 1.5.3 [64], under the scenarios of the constant and exponential population size as well as Bayesian skyline plot using a GTR+γ+Γ and HKY models with uncorrelated lognormal relaxed molecular clock [1]. Log likelihoods for the models presented were as follows: −5790.27 for the relaxed a GTR+γ+Γ model with constant population size, −5790.27 for the relaxed a GTR+γ+Γ model with an exponential prior, −5833.38 for relaxed a GTR+2 codon positions with constant population size, −5837.73 for the HKY model. Models were compared using Chi-square distribution test, allowing to identify that GTR+2 codon position and HKY models were significantly less fit than the GTR+γ+Γ. As a GTR+γ+Γ models with constant population size and exponential priors were not statistically different, we have selected the scenario with the lowest marginal likelihood (constant population size). A consensus tree with posterior probabilities for branch support was obtained and annotated with TreeAnnotator v. 1.5.4. All trees were visualized in Figtree v.1.2.2.

Statistics

Statistical comparisons were performed using the Fisher's exact t tests for nominal variables while continuous variables were analyzed using the U-Mann Whitney test using Statistica 8.0 PL software (Statasoft, Poland).

Sequence data

HIV-1 subtype D sequences from this study sequences were submitted to the GenBank and may be accessed with the following IDs: GU906860, GU906864, GU906871, GU906872, GU906873, GU906874, JQ305750–JQ305792.

Acknowledgments

The authors would like to express their gratitude to the staff of the Clinic for Acquired Immunodeficiency, Regional Hospital, Szczecin, Poland, for the help with data collection; and to the staff of the Research Laboratory, BC Centre for Excellence in HIV/AIDS, Vancouver, Canada, especially Richard Harrigan, who provided the V3 sequencing methodology and access to the Recall software.

Author Contributions

Conceived and designed the experiments: MP. Performed the experiments: MP AU. Analyzed the data: MP. Contributed reagents/materials/analysis tools: MP AU. Wrote the paper: MP MLP DB AU ABK.

Citation: Parczewski M, Leszczyszyn-Pynka M, Bander D, Urbanska A, Boroń-Kaczmarska A (2012) HIV-1 Subtype D Infections among Caucasians from Northwestern Poland—Phylogenetic and Clinical Analysis. PLoS ONE 7(2): e31674. https://doi.org/10.1371/journal.pone.0031674

References

1. Drummond AJ, Ho SY, Phillips MJ, Rambaut A (2006) Relaxed phylogenetics and dating with confidence. PLoS biology 4: e88.AJ DrummondSY HoMJ PhillipsA. Rambaut2006Relaxed phylogenetics and dating with confidence.PLoS biology4e88

2. 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.F. LewisGJ HughesA. RambautA. PozniakAJ Leigh Brown2008Episodic sexual transmission of HIV revealed by molecular phylodynamics.PLoS medicine5e50

3. Paraskevis D, Magiorkinis E, Magiorkinis G, Sypsa V, Paparizos V, et al. (2007) Increasing prevalence of HIV-1 subtype A in Greece: estimating epidemic history and origin. The Journal of infectious diseases 196: 1167–1176.D. ParaskevisE. MagiorkinisG. MagiorkinisV. SypsaV. Paparizos2007Increasing prevalence of HIV-1 subtype A in Greece: estimating epidemic history and origin.The Journal of infectious diseases19611671176

4. Geretti AM (2006) HIV-1 subtypes: epidemiology and significance for HIV management. Current opinion in infectious diseases 19: 1–7.AM Geretti2006HIV-1 subtypes: epidemiology and significance for HIV management.Current opinion in infectious diseases1917

5. Buonaguro L, Petrizzo A, Tagliamonte M, Vitone F, Re MC, et al. (2008) Molecular and phylogenetic analysis of HIV-1 variants circulating in Italy. Infectious agents and cancer 3: 13.L. BuonaguroA. PetrizzoM. TagliamonteF. VitoneMC Re2008Molecular and phylogenetic analysis of HIV-1 variants circulating in Italy.Infectious agents and cancer313

6. Vachot L, Ataman-Onal Y, Terrat C, Durand PY, Ponceau B, et al. (2004) Short communication: retrospective study to time the introduction of HIV type 1 non-B subtypes in Lyon, France, using env genes obtained from primary infection samples. AIDS research and human retroviruses 20: 687–691.L. VachotY. Ataman-OnalC. TerratPY DurandB. Ponceau2004Short communication: retrospective study to time the introduction of HIV type 1 non-B subtypes in Lyon, France, using env genes obtained from primary infection samples.AIDS research and human retroviruses20687691

7. Barlow KL, Tatt ID, Cane PA, Pillay D, Clewley JP (2001) Recombinant strains of HIV type 1 in the United Kingdom. AIDS research and human retroviruses 17: 467–474.KL BarlowID TattPA CaneD. PillayJP Clewley2001Recombinant strains of HIV type 1 in the United Kingdom.AIDS research and human retroviruses17467474

8. Holguin A, Alvarez A, Soriano V (2002) HIV-1 subtype J recombinant viruses in Spain. AIDS research and human retroviruses 18: 523–529.A. HolguinA. AlvarezV. Soriano2002HIV-1 subtype J recombinant viruses in Spain.AIDS research and human retroviruses18523529

9. Holguin A, Alvarez A, Soriano V (2002) High prevalence of HIV-1 subtype G and natural polymorphisms at the protease gene among HIV-infected immigrants in Madrid. AIDS 16: 1163–1170.A. HolguinA. AlvarezV. Soriano2002High prevalence of HIV-1 subtype G and natural polymorphisms at the protease gene among HIV-infected immigrants in Madrid.AIDS1611631170

10. Esteves A, Parreira R, Piedade J, Venenno T, Franco M, et al. (2003) Spreading of HIV-1 subtype G and envB/gagG recombinant strains among injecting drug users in Lisbon, Portugal. AIDS research and human retroviruses 19: 511–517.A. EstevesR. ParreiraJ. PiedadeT. VenennoM. Franco2003Spreading of HIV-1 subtype G and envB/gagG recombinant strains among injecting drug users in Lisbon, Portugal.AIDS research and human retroviruses19511517

11. Zetterberg V, Ustina V, Liitsola K, Zilmer K, Kalikova N, et al. (2004) Two viral strains and a possible novel recombinant are responsible for the explosive injecting drug use-associated HIV type 1 epidemic in Estonia. AIDS research and human retroviruses 20: 1148–1156.V. ZetterbergV. UstinaK. LiitsolaK. ZilmerN. Kalikova2004Two viral strains and a possible novel recombinant are responsible for the explosive injecting drug use-associated HIV type 1 epidemic in Estonia.AIDS research and human retroviruses2011481156

12. Liitsola K, Tashkinova I, Laukkanen T, Korovina G, Smolskaja T, et al. (1998) HIV-1 genetic subtype A/B recombinant strain causing an explosive epidemic in injecting drug users in Kaliningrad. AIDS 12: 1907–1919.K. LiitsolaI. TashkinovaT. LaukkanenG. KorovinaT. Smolskaja1998HIV-1 genetic subtype A/B recombinant strain causing an explosive epidemic in injecting drug users in Kaliningrad.AIDS1219071919

13. Patrascu IV, Dumitrescu O (1993) The epidemic of human immunodeficiency virus infection in Romanian children. AIDS research and human retroviruses 9: 99–104.IV PatrascuO. Dumitrescu1993The epidemic of human immunodeficiency virus infection in Romanian children.AIDS research and human retroviruses999104

14. Mehta SR, Wertheim JO, Delport W, Ene L, Tardei G, et al. (2011) Using phylogeography to characterize the origins of the HIV-1 subtype F epidemic in Romania. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases 11: 975–979.SR MehtaJO WertheimW. DelportL. EneG. Tardei2011Using phylogeography to characterize the origins of the HIV-1 subtype F epidemic in Romania.Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases11975979

15. Vergne L, Peeters M, Mpoudi-Ngole E, Bourgeois A, Liegeois F, et al. (2000) Genetic diversity of protease and reverse transcriptase sequences in non-subtype-B human immunodeficiency virus type 1 strains: evidence of many minor drug resistance mutations in treatment-naive patients. Journal of clinical microbiology 38: 3919–3925.L. VergneM. PeetersE. Mpoudi-NgoleA. BourgeoisF. Liegeois2000Genetic diversity of protease and reverse transcriptase sequences in non-subtype-B human immunodeficiency virus type 1 strains: evidence of many minor drug resistance mutations in treatment-naive patients.Journal of clinical microbiology3839193925

16. Frange P, Chaix ML, Raymond S, Galimand J, Deveau C, et al. (2010) Low frequency of CXCR4-using viruses in patients at the time of primary non-subtype-B HIV-1 infection. Journal of clinical microbiology 48: 3487–3491.P. FrangeML ChaixS. RaymondJ. GalimandC. Deveau2010Low frequency of CXCR4-using viruses in patients at the time of primary non-subtype-B HIV-1 infection.Journal of clinical microbiology4834873491

17. Garcia F, Perez-Cachafeiro S, Guillot V, Alvarez M, Perez-Romero P, et al. (2011) Transmission of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treatment naive HIV-infected individuals (CoRIS). Antiviral research 91: 150–153.F. GarciaS. Perez-CachafeiroV. GuillotM. AlvarezP. Perez-Romero2011Transmission of HIV drug resistance and non-B subtype distribution in the Spanish cohort of antiretroviral treatment naive HIV-infected individuals (CoRIS).Antiviral research91150153

18. Huang W, Eshleman SH, Toma J, Fransen S, Stawiski E, et al. (2007) Coreceptor tropism in human immunodeficiency virus type 1 subtype D: high prevalence of CXCR4 tropism and heterogeneous composition of viral populations. Journal of virology 81: 7885–7893.W. HuangSH EshlemanJ. TomaS. FransenE. Stawiski2007Coreceptor tropism in human immunodeficiency virus type 1 subtype D: high prevalence of CXCR4 tropism and heterogeneous composition of viral populations.Journal of virology8178857893

19. De Felipe B, Perez-Romero P, Abad-Fernandez M, Fernandez-Cuenca F, Martinez-Fernandez FJ, et al. (2011) Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000–2010. Virology journal 8: 416.B. De FelipeP. Perez-RomeroM. Abad-FernandezF. Fernandez-CuencaFJ Martinez-Fernandez2011Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000–2010.Virology journal8416

20. Geretti AM, Harrison L, Green H, Sabin C, Hill T, et al. (2009) Effect of HIV-1 subtype on virologic and immunologic response to starting highly active antiretroviral therapy. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America 48: 1296–1305.AM GerettiL. HarrisonH. GreenC. SabinT. Hill2009Effect of HIV-1 subtype on virologic and immunologic response to starting highly active antiretroviral therapy.Clinical infectious diseases: an official publication of the Infectious Diseases Society of America4812961305

21. Audelin AM, Gerstoft J, Obel N, Mathiesen L, Laursen A, et al. (2011) Molecular Phylogenetics of Transmitted Drug Resistance in Newly Diagnosed HIV Type 1 Individuals in Denmark, a Nation-Wide Study. AIDS research and human retroviruses 27(12): 1283–1290.AM AudelinJ. GerstoftN. ObelL. MathiesenA. Laursen2011Molecular Phylogenetics of Transmitted Drug Resistance in Newly Diagnosed HIV Type 1 Individuals in Denmark, a Nation-Wide Study.AIDS research and human retroviruses271212831290

22. Skoura L, Metallidis S, Buckton AJ, Mbisa JL, Pilalas D, et al. (2011) Molecular and epidemiological characterization of HIV-1 infection networks involving transmitted drug resistance mutations in Northern Greece. The Journal of antimicrobial chemotherapy 66(12): 2831–2837.L. SkouraS. MetallidisAJ BucktonJL MbisaD. Pilalas2011Molecular and epidemiological characterization of HIV-1 infection networks involving transmitted drug resistance mutations in Northern Greece.The Journal of antimicrobial chemotherapy661228312837

23. Vercauteren J, Wensing AM, van de Vijver DA, Albert J, Balotta C, et al. (2009) Transmission of drug-resistant HIV-1 is stabilizing in Europe. The Journal of infectious diseases 200: 1503–1508.J. VercauterenAM WensingDA van de VijverJ. AlbertC. Balotta2009Transmission of drug-resistant HIV-1 is stabilizing in Europe.The Journal of infectious diseases20015031508

24. SPREAD programme (2008) Transmission of drug-resistant HIV-1 in Europe remains limited to single classes. AIDS 22: 625–635.SPREAD programme2008Transmission of drug-resistant HIV-1 in Europe remains limited to single classes.AIDS22625635

25. Stanczak GP, Stanczak JJ, Firlag-Burkacka E, Wiercinska-Drapalo A, Leszczyszyn-Pynkad M, et al. (2007) [Transmission of HIV-1 drug resistance among newly diagnosed patients in Poland] Article in Polish. Przeglad epidemiologiczny 61: 29–34.GP StanczakJJ StanczakE. Firlag-BurkackaA. Wiercinska-DrapaloM. Leszczyszyn-Pynkad2007[Transmission of HIV-1 drug resistance among newly diagnosed patients in Poland] Article in Polish.Przeglad epidemiologiczny612934

26. Werbinska-Sienkiewicz B, Rosinska M, Furman S (2011) [HIV and AIDS in Poland in 2009]. Article in Polish. Przeglad epidemiologiczny 65: 289–295.B. Werbinska-SienkiewiczM. RosinskaS. Furman2011[HIV and AIDS in Poland in 2009]. Article in Polish.Przeglad epidemiologiczny65289295

27. Parczewski M, Leszczyszyn-Pynka M, Bander D, Urbanska A, Stanczak G, et al. (2010) Characteristics of HIV-1 non-B subtype infections in Northwest Poland. Journal of medical virology 82: 1306–1313.M. ParczewskiM. Leszczyszyn-PynkaD. BanderA. UrbanskaG. Stanczak2010Characteristics of HIV-1 non-B subtype infections in Northwest Poland.Journal of medical virology8213061313

28. Easterbrook PJ, Smith M, Mullen J, O'Shea S, Chrystie I, et al. (2010) Impact of HIV-1 viral subtype on disease progression and response to antiretroviral therapy. Journal of the International AIDS Society 13: 4.PJ EasterbrookM. SmithJ. MullenS. O'SheaI. Chrystie2010Impact of HIV-1 viral subtype on disease progression and response to antiretroviral therapy.Journal of the International AIDS Society134

29. Thomson MM, Najera R (2001) Travel and the introduction of human immunodeficiency virus type 1 non-B subtype genetic forms into Western countries. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America 32: 1732–1737.MM ThomsonR. Najera2001Travel and the introduction of human immunodeficiency virus type 1 non-B subtype genetic forms into Western countries.Clinical infectious diseases: an official publication of the Infectious Diseases Society of America3217321737

30. Bander D, Leszczyszyn-Pynka M, Boron-Kaczmarska A (2009) [Late AIDS diagnosis in patients hospitalized in Clinic of Infectious Diseases and Hepatology PAM in years 2003–2007]. Przeglad epidemiologiczny 63: 61–66.D. BanderM. Leszczyszyn-PynkaA. Boron-Kaczmarska2009[Late AIDS diagnosis in patients hospitalized in Clinic of Infectious Diseases and Hepatology PAM in years 2003–2007].Przeglad epidemiologiczny636166

31. Baeten JM, Chohan B, Lavreys L, Chohan V, McClelland RS, et al. (2007) HIV-1 subtype D infection is associated with faster disease progression than subtype A in spite of similar plasma HIV-1 loads. The Journal of infectious diseases 195: 1177–1180.JM BaetenB. ChohanL. LavreysV. ChohanRS McClelland2007HIV-1 subtype D infection is associated with faster disease progression than subtype A in spite of similar plasma HIV-1 loads.The Journal of infectious diseases19511771180

32. Kiwanuka N, Robb M, Laeyendecker O, Kigozi G, Wabwire-Mangen F, et al. (2010) HIV-1 viral subtype differences in the rate of CD4+ T-cell decline among HIV seroincident antiretroviral naive persons in Rakai district, Uganda. Journal of acquired immune deficiency syndromes 54: 180–184.N. KiwanukaM. RobbO. LaeyendeckerG. KigoziF. Wabwire-Mangen2010HIV-1 viral subtype differences in the rate of CD4+ T-cell decline among HIV seroincident antiretroviral naive persons in Rakai district, Uganda.Journal of acquired immune deficiency syndromes54180184

33. Kaleebu P, French N, Mahe C, Yirrell D, Watera C, et al. (2002) Effect of human immunodeficiency virus (HIV) type 1 envelope subtypes A and D on disease progression in a large cohort of HIV-1-positive persons in Uganda. The Journal of infectious diseases 185: 1244–1250.P. KaleebuN. FrenchC. MaheD. YirrellC. Watera2002Effect of human immunodeficiency virus (HIV) type 1 envelope subtypes A and D on disease progression in a large cohort of HIV-1-positive persons in Uganda.The Journal of infectious diseases18512441250

34. Bousheri S, Burke C, Ssewanyana I, Harrigan R, Martin J, et al. (2009) Infection with different hiv subtypes is associated with CD4 activation-associated dysfunction and apoptosis. Journal of acquired immune deficiency syndromes 52: 548–552.S. BousheriC. BurkeI. SsewanyanaR. HarriganJ. Martin2009Infection with different hiv subtypes is associated with CD4 activation-associated dysfunction and apoptosis.Journal of acquired immune deficiency syndromes52548552

35. Lutalo T, Gray RH, Wawer M, Sewankambo N, Serwadda D, et al. (2007) Survival of HIV-infected treatment-naive individuals with documented dates of seroconversion in Rakai, Uganda. AIDS 21:: Suppl 6S15–9.T. LutaloRH GrayM. WawerN. SewankamboD. Serwadda2007Survival of HIV-infected treatment-naive individuals with documented dates of seroconversion in Rakai, Uganda.AIDS21:Suppl 6S159

36. Kiwanuka N, Laeyendecker O, Robb M, Kigozi G, Arroyo M, et al. (2008) Effect of human immunodeficiency virus Type 1 (HIV-1) subtype on disease progression in persons from Rakai, Uganda, with incident HIV-1 infection. The Journal of infectious diseases 197: 707–713.N. KiwanukaO. LaeyendeckerM. RobbG. KigoziM. Arroyo2008Effect of human immunodeficiency virus Type 1 (HIV-1) subtype on disease progression in persons from Rakai, Uganda, with incident HIV-1 infection.The Journal of infectious diseases197707713

37. Parczewski M, Bander D, Leszczyszyn-Pynka M, Urbanska A, Kaczmarczyk M, et al. (2011) Risk of all-cause mortality in HIV infected patients is associated with clinical, immunologic predictors and the CCR5 Delta32 deletion. PloS one 6: e22215.M. ParczewskiD. BanderM. Leszczyszyn-PynkaA. UrbanskaM. Kaczmarczyk2011Risk of all-cause mortality in HIV infected patients is associated with clinical, immunologic predictors and the CCR5 Delta32 deletion.PloS one6e22215

38. Stanczak GP, Stanczak JJ, Marczynska M, Firlag-Burkacka E, Wiercinska-Drapalo A, 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.GP StanczakJJ StanczakM. MarczynskaE. Firlag-BurkackaA. Wiercinska-Drapalo2010Evolving patterns of HIV-1 transmitted drug resistance in Poland in the years 2000–2008.Journal of medical virology8212911294

39. 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.S. HueD. PillayJP ClewleyOG Pybus2005Genetic 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 America10244254429

40. Hughes GJ, Fearnhill E, Dunn D, Lycett SJ, Rambaut A, et al. (2009) Molecular phylodynamics of the heterosexual HIV epidemic in the United Kingdom. PLoS pathogens 5: e1000590.GJ HughesE. FearnhillD. DunnSJ LycettA. Rambaut2009Molecular phylodynamics of the heterosexual HIV epidemic in the United Kingdom.PLoS pathogens5e1000590

41. 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.M. KayeD. ChiboC. Birch2008Phylogenetic investigation of transmission pathways of drug-resistant HIV-1 utilizing pol sequences derived from resistance genotyping.Journal of acquired immune deficiency syndromes49916

42. Cuevas MT, Munoz-Nieto M, Thomson MM, Delgado E, Iribarren JA, et al. (2009) HIV-1 transmission cluster with T215D revertant mutation among newly diagnosed patients from the Basque Country, Spain. Journal of acquired immune deficiency syndromes 51: 99–103.MT CuevasM. Munoz-NietoMM ThomsonE. DelgadoJA Iribarren2009HIV-1 transmission cluster with T215D revertant mutation among newly diagnosed patients from the Basque Country, Spain.Journal of acquired immune deficiency syndromes5199103

43. Brenner BG, Roger M, Moisi DD, Oliveira M, Hardy I, et al. (2008) Transmission networks of drug resistance acquired in primary/early stage HIV infection. AIDS 22: 2509–2515.BG BrennerM. RogerDD MoisiM. OliveiraI. Hardy2008Transmission networks of drug resistance acquired in primary/early stage HIV infection.AIDS2225092515

44. Chalmet K, Staelens D, Blot S, Dinakis S, 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.K. ChalmetD. StaelensS. BlotS. DinakisJ. Pelgrom2010Epidemiological study of phylogenetic transmission clusters in a local HIV-1 epidemic reveals distinct differences between subtype B and non-B infections.BMC infectious diseases10262

45. Hue S, Gifford RJ, Dunn D, Fernhill E, Pillay D (2009) Demonstration of sustained drug-resistant human immunodeficiency virus type 1 lineages circulating among treatment-naive individuals. Journal of virology 83: 2645–2654.S. HueRJ GiffordD. DunnE. FernhillD. Pillay2009Demonstration of sustained drug-resistant human immunodeficiency virus type 1 lineages circulating among treatment-naive individuals.Journal of virology8326452654

46. Xiridou M, Geskus R, de Wit J, Coutinho R, Kretzschmar M (2004) Primary HIV infection as source of HIV transmission within steady and casual partnerships among homosexual men. AIDS 18: 1311–1320.M. XiridouR. GeskusJ. de WitR. CoutinhoM. Kretzschmar2004Primary HIV infection as source of HIV transmission within steady and casual partnerships among homosexual men.AIDS1813111320

47. Lai A, Violin M, Ebranati E, Franzetti M, Micheli V, et al. (2011) Transmission of Resistant Hiv-1 Variants and Epidemiological Chains in Italian Newly Diagnosed Individuals. AIDS research and human retroviruses. A. LaiM. ViolinE. EbranatiM. FranzettiV. Micheli2011Transmission of Resistant Hiv-1 Variants and Epidemiological Chains in Italian Newly Diagnosed Individuals.AIDS research and human retroviruses

48. Raymond S, Delobel P, Chaix ML, Cazabat M, Encinas S, et al. (2011) Genotypic prediction of HIV-1 subtype D tropism. Retrovirology 8: 56.S. RaymondP. DelobelML ChaixM. CazabatS. Encinas2011Genotypic prediction of HIV-1 subtype D tropism.Retrovirology856

49. Kaleebu P, Nankya IL, Yirrell DL, Shafer LA, Kyosiimire-Lugemwa J, et al. (2007) Relation between chemokine receptor use, disease stage, and HIV-1 subtypes A and D: results from a rural Ugandan cohort. Journal of acquired immune deficiency syndromes 45: 28–33.P. KaleebuIL NankyaDL YirrellLA ShaferJ. Kyosiimire-Lugemwa2007Relation between chemokine receptor use, disease stage, and HIV-1 subtypes A and D: results from a rural Ugandan cohort.Journal of acquired immune deficiency syndromes452833

50. Reekie J, Mocroft A, Ledergerber B, Beniowski M, Clotet B, et al. (2010) History of viral suppression on combination antiretroviral therapy as a predictor of virological failure after a treatment change. HIV medicine 11: 469–478.J. ReekieA. MocroftB. LedergerberM. BeniowskiB. Clotet2010History of viral suppression on combination antiretroviral therapy as a predictor of virological failure after a treatment change.HIV medicine11469478

51. Sacktor N, Nakasujja N, Skolasky RL, Rezapour M, Robertson K, et al. (2009) HIV subtype D is associated with dementia, compared with subtype A, in immunosuppressed individuals at risk of cognitive impairment in Kampala, Uganda. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America 49: 780–786.N. SacktorN. NakasujjaRL SkolaskyM. RezapourK. Robertson2009HIV subtype D is associated with dementia, compared with subtype A, in immunosuppressed individuals at risk of cognitive impairment in Kampala, Uganda.Clinical infectious diseases: an official publication of the Infectious Diseases Society of America49780786

52. Lamers SL, Poon AF, McGrath MS (2011) HIV-1 nef protein structures associated with brain infection and dementia pathogenesis. PloS one 6: e16659.SL LamersAF PoonMS McGrath2011HIV-1 nef protein structures associated with brain infection and dementia pathogenesis.PloS one6e16659

53. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular biology and evolution 24: 1596–1599.K. TamuraJ. DudleyM. NeiS. Kumar2007MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0.Molecular biology and evolution2415961599

54. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, 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.KS LoleRC BollingerRS ParanjapeD. GadkariSS Kulkarni1999Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination.Journal of virology73152160

55. Van Laethem K, Schrooten Y, Covens K, Dekeersmaeker N, De Munter P, et al. (2008) A genotypic assay for the amplification and sequencing of integrase from diverse HIV-1 group M subtypes. Journal of virological methods 153: 176–181.K. Van LaethemY. SchrootenK. CovensN. DekeersmaekerP. De Munter2008A genotypic assay for the amplification and sequencing of integrase from diverse HIV-1 group M subtypes.Journal of virological methods153176181

56. Harrigan PRDW, Wynhoven B, Mo T, Hall T, Galli R (2002) The Performance of ReCall basecalling software for high-throughput HIV drug resistance basecalling using “in-house” methods. PRDW HarriganB. WynhovenT. MoT. HallR. Galli2002The Performance of ReCall basecalling software for high-throughput HIV drug resistance basecalling using “in-house” methods.XIV International AIDS Conference, Spain, Barcelona Abstract #TuPeB4598. XIV International AIDS Conference, Spain, Barcelona Abstract #TuPeB4598.

57. Vandekerckhove LP, Wensing AM, Kaiser R, Brun-Vezinet F, Clotet B, et al. (2011) European guidelines on the clinical management of HIV-1 tropism testing. The Lancet infectious diseases 11: 394–407.LP VandekerckhoveAM WensingR. KaiserF. Brun-VezinetB. Clotet2011European guidelines on the clinical management of HIV-1 tropism testing.The Lancet infectious diseases11394407

58. Swenson LC, Mo T, Dong WW, Zhong X, Woods CK, et al. (2011) Deep Third Variable Sequencing for HIV Type 1 Tropism in Treatment-Naive Patients: A Reanalysis of the MERIT Trial of Maraviroc. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America 53: 732–742.LC SwensonT. MoWW DongX. ZhongCK Woods2011Deep Third Variable Sequencing for HIV Type 1 Tropism in Treatment-Naive Patients: A Reanalysis of the MERIT Trial of Maraviroc.Clinical infectious diseases: an official publication of the Infectious Diseases Society of America53732742

59. Lengauer T, Sander O, Sierra S, Thielen A, Kaiser R (2007) Bioinformatics prediction of HIV coreceptor usage. Nature biotechnology 25: 1407–1410.T. LengauerO. SanderS. SierraA. ThielenR. Kaiser2007Bioinformatics prediction of HIV coreceptor usage.Nature biotechnology2514071410

60. Liu TF, Shafer RW (2006) Web resources for HIV type 1 genotypic-resistance test interpretation. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America 42: 1608–1618.TF LiuRW Shafer2006Web resources for HIV type 1 genotypic-resistance test interpretation.Clinical infectious diseases: an official publication of the Infectious Diseases Society of America4216081618

61. Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, et al. (2009) Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update. PloS one 4: e4724.DE BennettRJ CamachoD. OteleaDR KuritzkesH. Fleury2009Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update.PloS one4e4724

62. Shafer RW, Rhee SY, Pillay D, Miller V, Sandstrom P, et al. (2007) HIV-1 protease and reverse transcriptase mutations for drug resistance surveillance. AIDS 21: 215–223.RW ShaferSY RheeD. PillayV. MillerP. Sandstrom2007HIV-1 protease and reverse transcriptase mutations for drug resistance surveillance.AIDS21215223

63. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, et al. (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic biology 59: 307–321.S. GuindonJF DufayardV. LefortM. AnisimovaW. Hordijk2010New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0.Systematic biology59307321

64. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC evolutionary biology 7: 214.AJ DrummondA. Rambaut2007BEAST: Bayesian evolutionary analysis by sampling trees.BMC evolutionary biology7214