HIV-1 pol gene variation is sufficient for reconstruction of transmissions in the era of antiretroviral therapy

Genetic characterization of human immunodeficiency virus type 1 transmission in the Middle East and North Africa

BACKGROUND

The HIV-1 spread in the Middle East and North Africa (MENA) has not been previously characterised using the phylogenetic approach. The aim of the current study was to investigate the genetic diversity and domestic transmission of HIV-1 in the MENA.

METHODS

A total of 2036 HIV-1 sequences available in Genbank and collected in the MENA during 1988-2016 were used together with 715 HIV-1 reference sequences that were retrieved from Genbank based on genetic similarity with the MENA sequences. The REGA and COMET tools were used to determine HIV-1 subtypes and circulating recombinant forms. Maximum Likelihood and Bayesian phylogenetic analyses were used to identify and date HIV-1 transmission clusters.

RESULTS

At least 21 HIV-1 subtypes and recombinant forms were prevalent in the MENA. Subtype B was the most common variant (39%), followed by CRF35_AD (19%) and CRF02_AG (14%). The most common genetic region was pol, and 675 partial pol sequences (average of 1005 bp) were eligible for detailed phylogenetic analysis. Fifty-four percent of the MENA sequences formed HIV-1 transmission clusters. Whereas numerous clusters were country-specific, some clusters indicated transmission links between countries for subtypes B, C and CRF02_AG. This was more common in North Africa compared with the Middle East (p < 0.001). Recombinant forms had a larger proportion of clustering compared to pure subtypes (p < 0.001). The largest MENA clusters dated back to 1991 (an Algerian CRF06_cpx cluster of 43 sequences) and 2002 (a Tunisian CRF02_AG cluster of 48 sequences).

CONCLUSIONS

We found an extensive HIV-1 diversity in the MENA and a high proportion of sequences in transmission clusters. This study highlights the need for preventive measures in the MENA to limit HIV-1 spread in this region.

Defining HIV-1 transmission clusters based on sequence data

Understanding HIV-1 transmission dynamics is relevant to both screening and intervention strategies of HIV-1 infection. Commonly, HIV-1 transmission chains are determined based on sequence similarity assessed either directly from a sequence alignment or by inferring a phylogenetic tree. This review is aimed at both nonexperts interested in understanding and interpreting studies of HIV-1 transmission, and experts interested in finding the most appropriate cluster definition for a specific dataset and research question. We start by introducing the concepts and methodologies of how HIV-1 transmission clusters usually have been defined. We then present the results of a systematic review of 105 HIV-1 molecular epidemiology studies summarizing the most common methods and definitions in the literature. Finally, we offer our perspectives on how HIV-1 transmission clusters can be defined and provide some guidance based on examples from real life datasets.

Increase in human immunodeficiency virus 1 diversity and detection of various subtypes and recombinants in north-eastern Brazil

PURPOSE

Diverse human immunodeficiency virus 1 (HIV-1) subtypes and circulating recombinant forms are found in Brazil. The majority of HIV-1 molecular epidemiological studies in Brazil have been conducted in the southern and south-eastern regions of the country, although several recent studies in the north-eastern region have addressed this issue. The objective of this study was to molecularly characterize HIV-1 circulating in Pernambuco, north-eastern Brazil.

METHODOLOGY

A total of 64 samples were collected from 2002 to 2003, and another 103 were collected from 2007 to 2009. The protease and partial reverse transcriptase regions of the HIV-1 polymerase-encoding (pol) gene were sequenced, and subtyping, recombination and phylogenetic analyses were performed.Results/Key findings. Subtype B (60.9 %) was found to be predominant, followed by HIV-1 F (31.4 %). Several BF recombinants (4.2 %), and BC and AG recombinants were also identified. The intra-subtype genetic diversity was estimated to be 0.065 (sd±0.004) for HIV-1 B and 0.055 (sd±0.004) for HIV-1 F, reflecting a greater accumulation of mutations in subtype B (P<0.01). More codons were found to be under positive selective pressure in samples collected from 2007 to 2009, from individuals with a T-cell count≥200 cells mm-3 and from women. Coalescence data indicated that the subtype F population has been continuously expanding.

CONCLUSIONS

HIV-1 shows high genetic diversity in the state of Pernambuco. Thus, additional molecular evaluations of circulating strains will provide a better understanding of the epidemic and may lead to more effective preventive strategies.

The multi-faceted dynamics of HIV-1 transmission in Northern Alberta: A combined analysis of virus genetic and public health data

Molecular epidemiology has become a key tool for tracking infectious disease epidemics. Here, the spread of the most prevalent HIV-1 subtypes in Northern Alberta, Canada, was characterized with a Bayesian phylogenetic approach using 1146 HIV-1 pol sequences collected between 2007 and 2013 for routine clinical management purposes. Available patient metadata were qualitatively interpreted and correlated with onwards transmission using Fisher exact tests and logistic regression. Most infections were from subtypes A (n=36), B (n=815) and C (n=211). Africa is the dominant origin location for subtypes A and C while the subtype B epidemic was seeded from the USA and Middle America and, from the early 1990s onwards, mostly by interprovincial spread. Subtypes A (77.8%) and C (74.0%) were usually heterosexually transmitted and circulate predominantly among Blacks (61.1% and 85% respectively). Subtype B was mostly found among Caucasians (48.6%) and First Nations (36.8%), and its modes of transmission were stratified by ethnic origin. Compared to subtypes A (5.6%) and C (3.8-10.0%), a larger portion of subtype B patients were found within putative provincial transmission networks (20.3-29.5%), and this almost doubled when focusing on nationwide transmission clusters (37.9-57.5%). No clear association between cluster membership and particular patient characteristics was found. This study reveals complex and multi-faceted transmission dynamics of the HIV-1 epidemic in this otherwise low HIV prevalence population in Northern Alberta, Canada. These findings can aid public health planning.

Short Communication: Phylogenetic Evidence of HIV-1 Transmission Between Adult and Adolescent Men Who Have Sex with Men

HIV-1 incidence among youth, especially men who have sex with men (MSM), is increasing in the United States. We aimed to better understand the patterns of adolescent HIV-1 acquisition, to help guide future prevention interventions. We conducted a study combining epidemiologic and HIV-1 pol sequence data from a retrospective cohort of HIV-infected adults and adolescents in Seattle, WA between 2000 and 2013. Adolescents were defined as 13-24 years of age at the time of first HIV-1 care. Maximum-likelihood phylogenetic trees were reconstructed to identify putative viral transmission clusters of two or more individuals, followed by multivariable regression tests of associations between clustering and demographic and clinical parameters. The dataset included 3,102 sequences from 1,953 individuals; 72 putative transmission clusters were identified, representing 168 individuals (8.6%). MSM and MSM/intravenous drug use (IDU) were positively associated with clustering, with aOR 3.18 (95% CI: 1.34-7.55) and 2.59 (95% CI: 1.04-6.49), respectively. African American race was negatively associated with clustering (aOR 0.54 95% CI: 0.32-0.91). Twenty-five clusters contained one adolescent and five clusters contained two adolescents. Other individuals who clustered with adolescents were predominantly male (95%), white (85%), and either MSM (66%) or MSM/IDU (16%), with a greater mean age (34 years vs. 22 years; p < .01). In this Seattle cohort, HIV-1 transmission linkages were identified between white male adolescents and older MSM adults. Interventions aimed at age-discrepant pairs may reduce HIV-1 infections in adolescent males.

Identifying Transmission Clusters with Cluster Picker and HIV-TRACE

We compared the behavior of two approaches (Cluster Picker and HIV-TRACE) at varying genetic distances to identify transmission clusters. We used three HIV gp41 sequence datasets originating from the Rakai Community Cohort Study: (1) next-generation sequence (NGS) data from nine linked couples; (2) NGS data from longitudinal sampling of 14 individuals; and (3) Sanger consensus sequences from a cross-sectional dataset (n = 1,022) containing 91 epidemiologically linked heterosexual couples. We calculated the optimal genetic distance threshold to separate linked versus unlinked NGS datasets using a receiver operating curve analysis. We evaluated the number, size, and composition of clusters detected by Cluster Picker and HIV-TRACE at six genetic distance thresholds (1%-5.3%) on all three datasets. We further tested the effect of using all NGS, versus only a single variant for each patient/time point, for datasets (1) and (2). The optimal gp41 genetic distance threshold to distinguish linked and unlinked couples and individuals was 5.3% and 4%, respectively. HIV-TRACE tended to detect larger and fewer clusters, whereas Cluster Picker detected more clusters containing only two sequences. For NGS datasets (1) and (2), HIV-TRACE and Cluster Picker detected all linked pairs at 3% and 4% genetic distances, respectively. However, at 5.3% genetic distance, 20% of couples in dataset (3) did not cluster using either program, and for >1/3 of couples cluster assignment were discordant. We suggest caution in choosing thresholds for clustering analyses in a generalized epidemic.

Identification of a large, fast-expanding HIV-1 subtype B transmission cluster among MSM in Valencia, Spain

We describe and characterize an exceptionally large HIV-1 subtype B transmission cluster occurring in the Comunidad Valenciana (CV, Spain). A total of 1806 HIV-1 protease-reverse transcriptase (PR/RT) sequences from different patients were obtained in the CV between 2004 and 2014. After subtyping and generating a phylogenetic tree with additional HIV-1 subtype B sequences, a very large transmission cluster which included almost exclusively sequences from the CV was detected (n = 143 patients). This cluster was then validated and characterized with further maximum-likelihood phylogenetic analyses and Bayesian coalescent reconstructions. With these analyses, the CV cluster was delimited to 113 patients, predominately men who have sex with men (MSM). Although it was significantly located in the city of Valencia (n = 105), phylogenetic analyses suggested this cluster derives from a larger HIV lineage affecting other Spanish localities (n = 194). Coalescent analyses estimated its expansion in Valencia to have started between 1998 and 2004. From 2004 to 2009, members of this cluster represented only 1.46% of the HIV-1 subtype B samples studied in Valencia (n = 5/143), whereas from 2010 onwards its prevalence raised to 12.64% (n = 100/791). In conclusion, we have detected a very large transmission cluster in the CV where it has experienced a very fast growth in the recent years in the city of Valencia, thus contributing significantly to the HIV epidemic in this locality. Its transmission efficiency evidences shortcomings in HIV control measures in Spain and particularly in Valencia.

Using nearly full-genome HIV sequence data improves phylogeny reconstruction in a simulated epidemic

HIV molecular epidemiology studies analyse viral pol gene sequences due to their availability, but whole genome sequencing allows to use other genes. We aimed to determine what gene(s) provide(s) the best approximation to the real phylogeny by analysing a simulated epidemic (created as part of the PANGEA_HIV project) with a known transmission tree. We sub-sampled a simulated dataset of 4662 sequences into different combinations of genes (gag-pol-env, gag-pol, gag, pol, env and partial pol) and sampling depths (100%, 60%, 20% and 5%), generating 100 replicates for each case. We built maximum-likelihood trees for each combination using RAxML (GTR + Γ), and compared their topologies to the corresponding true tree’s using CompareTree. The accuracy of the trees was significantly proportional to the length of the sequences used, with the gag-pol-env datasets showing the best performance and gag and partial pol sequences showing the worst. The lowest sampling depths (20% and 5%) greatly reduced the accuracy of tree reconstruction and showed high variability among replicates, especially when using the shortest gene datasets. In conclusion, using longer sequences derived from nearly whole genomes will improve the reliability of phylogenetic reconstruction. With low sample coverage, results can be highly variable, particularly when based on short sequences.

Concurrency and HIV transmission network characteristics among MSM with recent HIV infection

OBJECTIVE

Sexual partner concurrency is common among MSM and may increase the probability of HIV transmission during recent (acute or early) infection. We examined the relationship between concurrency and HIV transmission network characteristics (proxies for HIV transmission) among MSM with recent HIV infection.

DESIGN

Observational study integrating behavioral, clinical, and molecular epidemiology.

METHODS

We inferred a partial HIV transmission network using 986 HIV-1 pol sequences obtained from HIV-infected individuals in San Diego, California (1996-2015). We further analyzed data from 285 recently HIV-infected MSM in the network who provided information on up to three sexual partners in the past 3 months, including the timing of intercourse with each partner. Concurrency was defined as sexual partners overlapping in time. Logistic and negative binomial regressions were used to investigate the link between concurrency and HIV transmission network characteristics (i.e. clustering and degree or number of connections to others in the network) among these MSM.

RESULTS

Of recently HIV-infected MSM (n = 285), 54% reported concurrent partnerships and 54% were connected by at least one putative transmission link to others (i.e. clustered) in the network (median degree = 1.0; interquartile range: 0.0-3.0). Concurrency was positively associated with HIV transmission network clustering (adjusted odds ratio = 1.83, 95% confidence interval: 1.08, 3.10) and degree (adjusted incidence rate ratio = 1.48, 95% confidence interval: 1.02, 2.15).

CONCLUSION

Our findings provide empirical evidence consistent with the hypothesis that concurrency facilitates HIV transmission during recent infection. Interventions to mitigate the impact of concurrency on HIV transmission may help curb the HIV epidemic among MSM.

Impacts and shortcomings of genetic clustering methods for infectious disease outbreaks

For infectious diseases, a genetic cluster is a group of closely related infections that is usually interpreted as representing a recent outbreak of transmission. Genetic clustering methods are becoming increasingly popular for molecular epidemiology, especially in the context of HIV where there is now considerable interest in applying these methods to prioritize groups for public health resources such as pre-exposure prophylaxis. To date, genetic clustering has generally been performed with ad hoc algorithms, only some of which have since been encoded and distributed as free software. These algorithms have seldom been validated on simulated data where clusters are known, and their interpretation and similarities are not transparent to users outside of the field. Here, I provide a brief overview on the development and inter-relationships of genetic clustering methods, and an evaluation of six methods on data simulated under an epidemic model in a risk-structured population. The simulation analysis demonstrates that the majority of clustering methods are systematically biased to detect variation in sampling rates among subpopulations, not variation in transmission rates. I discuss these results in the context of previous work and the implications for public health applications of genetic clustering.

Coreceptor usage of Chinese HIV-1 and impact of X4/DM transmission clusters among recently infected men who have sex with men

To characterize the current frequency of HIV-1 coreceptor usage in China and assess the candidacy of CCR5 antagonists for treatment of HIV infections. In addition, we aimed to evaluate the potential of X4/DM virus transmission in recently infected men who have sex with men (MSM) individuals.Viral tropism testing was performed on samples from 399 MSM individuals and on 2408 available Chinese HIV-1 V3 sequences downloaded from the Los Alamos database using Geno2pheno and WebPSSM in combination. The transmission clusters were evaluated using pol sequences from 291 recently infected MSM with a maximum likelihood, maximum pairwise distance, and Bayesian inference.A higher prevalence of X4/DM viruses was observed in individuals infected with CRF01_AE strains than with subtype B (27.8% vs 12.2%, P < 0.001) and CRF07_BC/CRF08_BC/C (27.8% vs 1.0%, P < 0.001). Seven clusters containing only X4/DM viruses were detected in 40 transmission clusters. No significant difference in proportions between clustered X4/DM viruses and R5 viruses was found (P = 0.683).The high proportion of CXCR4 usage for CRF01_AE strains may result in the loss of susceptibility to maraviroc since CRF01_AE has become the most prevalent strains in China. The high prevalence of X4/DM viruses among recently CRF01_AE-infected individuals may be attributed to the stochasticity of HIV transmission, which implied that the early viral tropism screening and treatment would be the key for controlling the epidemic of CRF01_AE strains in China.

Deep sequencing increases hepatitis C virus phylogenetic cluster detection compared to Sanger sequencing

Effective surveillance and treatment strategies are required to control the hepatitis C virus (HCV) epidemic. Phylogenetic analyses are powerful tools for reconstructing the evolutionary history of viral outbreaks and identifying transmission clusters. These studies often rely on Sanger sequencing which typically generates a single consensus sequence for each infected individual. For rapidly mutating viruses such as HCV, consensus sequencing underestimates the complexity of the viral quasispecies population and could therefore generate different phylogenetic tree topologies. Although deep sequencing provides a more detailed quasispecies characterization, in-depth phylogenetic analyses are challenging due to dataset complexity and computational limitations. Here, we apply deep sequencing to a characterized population to assess its ability to identify phylogenetic clusters compared with consensus Sanger sequencing. For deep sequencing, a sample specific threshold determined by the 50th percentile of the patristic distance distribution for all variants within each individual was used to identify clusters. Among seven patristic distance thresholds tested for the Sanger sequence phylogeny ranging from 0.005-0.06, a threshold of 0.03 was found to provide the maximum balance between positive agreement (samples in a cluster) and negative agreement (samples not in a cluster) relative to the deep sequencing dataset. From 77 HCV seroconverters, 10 individuals were identified in phylogenetic clusters using both methods. Deep sequencing analysis identified an additional 4 individuals and excluded 8 other individuals relative to Sanger sequencing. The application of this deep sequencing approach could be a more effective tool to understand onward HCV transmission dynamics compared with Sanger sequencing, since the incorporation of minority sequence variants improves the discrimination of phylogenetically linked clusters.

Genetic diversity and transmission networks of HIV-1 strains among men having sex with men (MSM) in Lomé, Togo

Understanding the HIV epidemic in key populations is important. Today only scarce information is available on HIV-1 strains that circulate in men having sex with men (MSM) in sub-Saharan Africa. Here, we studied for the first time the genetic diversity of HIV-1 strains circulating in the MSM population in Lomé, the capital city from Togo. The overall subtype/CRF distribution in pol (protease and/or partial reverse transcriptase (RT)) among the 79 HIV-1 strains from MSM was as follows: CRF02_AG (72%, n=57), subtype G (2.5%, n=2), sub-subtype A3 (1.3%, n=1), and unique recombinant forms (URF) (24%, n=19). Among the 19 URFs four different mosaic structures were observed, annotated as URF1 to URF4. Fifteen sequences (URF1) had the same mosaic structure in pol (G/CRF02_AG) and could represent a new circulating recombinant form (CRF). Phylogenetic analysis of the RT sequences showed that there were several introductions of CRF02_AG strains in the MSM population, however half of the CRF02_AG and all URF1 strains formed a separate, well-supported cluster suggesting one major introduction of CRF02_AG in the MSM population followed by efficient transmission and emergence of a possible new CRF. At least 40% of the strains fell into recent transmission chains involving two to seven MSM. Comparison with >950 HIV-1 sequences from previous studies in Togo showed intermixing of the HIV-1 epidemics between MSM and the general population. Moreover, an HIV-1 strain from a recently HIV-1 infected male patient from Germany, fell within a cluster of HIV-1 strains from MSM from Togo, illustrating recent exchange between MSM from Africa and people from other geographic regions. With growing evidence of the importance of MSM in the dynamic of the HIV epidemic in Africa there is an urgent need for appropriate interventions to limit HIV transmission in this population group.

Epidemiological, Clinical and Antiretroviral Susceptibility Characterization of Human Immunodeficiency Virus Subtypes B and Non-B in Pernambuco, Northeast Brazil

Background

HIV-1 diversity causes important differences in the virus’ biological properties and their interactions with hosts, such as cell tropism, responses to antiretroviral therapy, drug-resistance, and disease progression.

Objectives

We evaluated the interrelationship of phylogenetic inference with epidemiological and laboratory data for HIV-1 isolates circulating in Pernambuco, Northeast Region—Brazil.

Study design

A total of 168 HIV-1 pol sequences were analysed, 64 were obtained from 2002–2003, and 104, from 2007–2009. Socio-demographic, clinical, and behavioural data were obtained from medical records. Laboratory testing enabled the determination of recent HIV-1 infections and co-infections with HBV, HCV, HTLV, or syphilis. Surveillance drug-resistance mutation analysis and antiretroviral susceptibility profiling were performed using HIV Drug-Resistance Database.

Results

HIV-1 non-B was associated with female, lower education, lower viral loads, and higher T cell counts mean. Frequencies of co-infection HIV-HBV, HIV-HCV, and HIV-syphilis were 27.8% (95% CI: 19.8–37.7), 1.04% (95% CI: 0.05–5.00) and 14.7% (95% CI: 8.6–23.0), respectively. Drug-resistant mutations rate was 2.98% (95% CI: 1.10–6.47). HIV-HBV subtype B co-infection was associated with men who have sex with men (MSM), higher education, higher viral loads and males. HIV-syphilis subtype non-B co-infection was associated with MSM status, lower T cell counts and males.

Conclusions

Data showed the importance of molecular characterisations of the HIV-1 epidemic and its relation with epidemiological and clinical characteristics of the population, as well as its association with other infectious diseases, so they can effort to improve preventive measures for health services and more information about the progress and effects of the epidemic in Northeastern–Brazil.

Near real-time monitoring of HIV transmission hotspots from routine HIV genotyping: an implementation case study

BACKGROUND

HIV evolves rapidly and therefore infections with similar genetic sequences are likely linked by recent transmission events. Clusters of related infections can represent subpopulations with high rates of transmission. We describe the implementation of an automated near real-time system to monitor and characterise HIV transmission hotspots in British Columbia, Canada.

METHODS

In this implementation case study, we applied a monitoring system to the British Columbia drug treatment database, which holds more than 32 000 anonymised HIV genotypes for nearly 9000 residents of British Columbia living with HIV. On average, five to six new HIV genotypes are deposited in the database every day, which triggers an automated reanalysis of the entire database. We extracted clusters of five or more individuals with short phylogenetic distances between their respective HIV sequences. The system generated monthly reports of the growth and characteristics of clusters that were distributed to public health officers.

FINDINGS

In June, 2014, the monitoring system detected the expansion of a cluster by 11 new cases during 3 months, including eight cases with transmitted drug resistance. This cluster generally comprised young men who have sex with men. The subsequent report precipitated an enhanced public health follow-up to ensure linkage to care and treatment initiation in the affected subpopulation. Of the nine cases associated with this follow-up, all had already been linked to care and five cases had started treatment. Subsequent to the follow-up, three additional cases started treatment and most cases achieved suppressed viral loads. During the next 12 months, we detected 12 new cases in this cluster with reduction in the onward transmission of drug resistance.

INTERPRETATION

Our findings show the first application of an automated phylogenetic system monitoring a clinical database to detect a recent HIV outbreak and support the ensuing public health response. By making secondary use of routinely collected HIV genotypes, this approach is cost-effective, attains near real-time monitoring of new cases, and can be implemented in all settings in which HIV genotyping is the standard of care.

FUNDING

BC Centre for Excellence in HIV/AIDS, the Canadian Institutes for Health Research, the Genome Canada-CIHR Partnership in Genomics and Personalized Health, and the US National Institute on Drug Abuse.

Transmission dynamics of HIV-1 subtype B in the Basque Country, Spain

This work was aimed to study the HIV-1 subtype B epidemics in the Basque Country, Spain. 1727 HIV-1 subtype B sequences comprising protease and reverse transcriptase (PR/RT) coding regions, sampled between 2001 and 2008, were analyzed. 156 transmission clusters were detected by means of phylogenetic analyses. Most of them comprised less than 4 individuals and, in total, they included 441 patients. Six clusters comprised 10 or more patients and were further analyzed in order to study their origin and diversification. Four clusters included men who had unprotected homosexual sex (MSM), one group was formed by intravenous drug users (IDUs), and another included both IDUs and people infected through unprotected heterosexual sex (HTs). Most of these clusters originated from the mid-1980s to the mid-1990s. Only one cluster, formed by MSM, originated after 2000. The time between infections was significantly lower in MSM groups than in those containing IDUs (P-value <0.0001). Nucleoside RT and non-nucleoside RT inhibitor (NRTI and NNRTI)-resistance mutations to antiretroviral treatment were found in these six clusters except the most recent MSM group, but only the IDU clusters presented protease inhibitor (PI)-resistance mutations. The most prevalent mutations for each inhibitor class were PI L90M, NRTI T215D/Y/F, and NNRTI K103N, which were also among the most prevalent resistant variants in the whole dataset. In conclusion, while most infections occur as isolated introductions into the population, the number of infections found to be epidemiologically related within the Basque Country is significant. Public health control measures should be reinforced to prevent the further expansion of transmission clusters and resistant mutations occurring within them.

Prevalence of different HIV-1 subtypes in sexual transmission in China: a systematic review and meta-analysis

Sexual transmission has become the primary route of HIV transmission in China. Therefore, a comprehensive overview of HIV-1 subtype distribution is necessary for the prevention and control of the HIV epidemic. The present study aimed to provide a comprehensive prevalence estimate of different HIV-1 subtypes in sexual transmission in China. We conducted a systematic literature review for studies of HIV-1 subtypes in English and Chinese through several databases. Eligible articles were screened and selected by two authors independently. Random-effects model were applied to calculate the pooled prevalence of different HIV-1 subtypes, and subgroup analyses examined prevalence estimates across time, locations, and populations. A total of 130 eligible studies were identified, including 18 752 successfully genotyped samples. The pooled prevalence of CRF01_AE, subtype B, CRF07_BC, CRF08_BC, and subtype C were 44·54% (95% CI 40·81-48·30), 18·31% (95% CI 14·71-22·17), 16·45% (95% CI 13·82-19·25), 2·55% (95% CI 1·56-3·73), 0·37% (95% CI 0·11-0·72), respectively. The prevalence of subtype B in sexual transmission decreased, while the prevalence of CRF01_AE and CRF07_BC in sexual transmission, and CRF08_BC in heterosexual transmission increased. There is significant variation in HIV-1 subtype distribution between regions. The distribution of HIV-1 subtypes and circulating recombinant forms have changed significantly. The high genetic variability of HIV-1 poses a significant challenge for disease control and surveillance in China.

HIV-1 transmission between MSM and heterosexuals, and increasing proportions of circulating recombinant forms in the Nordic Countries

Increased knowledge about HIV-1 transmission dynamics in different transmission groups and geographical regions is fundamental for assessing and designing prevention efforts against HIV-1 spread. Since the first reported cases of HIV infection during the early 1980s, the HIV-1 epidemic in the Nordic countries has been dominated by HIV-1 subtype B and MSM transmission. HIV-1 pol sequences and clinical data of 51 per cent of all newly diagnosed HIV-1 infections in Sweden, Denmark, and Finland in the period 2000-2012 (N = 3,802) were analysed together with a large reference sequence dataset (N = 4,537) by trend analysis and phylogenetics. Analysis of the eight dominating subtypes and CRFs in the Nordic countries (A, B, C, D, G, CRF01_AE, CRF02_AG, and CRF06_cpx) showed that the subtype B proportion decreased while the CRF proportion increased over the study period. A majority (57 per cent) of the Nordic sequences formed transmission clusters, with evidence of mixing both geographically and between transmission groups. Detailed analyses showed multiple occasions of transmissions from MSM to heterosexuals and that active transmission clusters more often involved single than multiple Nordic countries. The strongest geographical link was between Denmark and Sweden. Finally, Denmark had a larger proportion of heterosexual domestic spread of HIV-1 subtype B (75 per cent) compared with Sweden (49 per cent) and Finland (57 per cent). We describe different HIV-1 transmission patterns between countries and transmission groups in a large geographical region. Our results may have implications for public health interventions in targeting HIV-1 transmission networks and identifying where to introduce such interventions.

Phylodynamic analysis of HIV sub-epidemics in Mochudi, Botswana

Southern Africa continues to be the epicenter of the HIV/AIDS epidemic. This HIV-1 subtype C epidemic has a predominantly heterosexual mode of virus transmission and high (>15%) HIV prevalence among adults. The epidemiological dynamics of the HIV-1C epidemic in southern Africa are still poorly understood. Here, we aim at a better understanding of HIV transmission dynamics by analyzing HIV-1 subtype C sequences from Mochudi, a peri-urban village in Botswana. HIV-1C env gene sequences (gp120 V1C5) were obtained through enhanced household-based HIV testing and counseling in Mochudi. More than 1200 sequences were generated and phylogenetically distinct sub-epidemics within Mochudi identified. The Bayesian birth-death skyline plot was used to estimate the effective reproductive number, R, and the timing of virus transmission, to classify sub-epidemics as "acute" (those with recent viral transmissions) or "historic" (those without recent viral transmissions). We identified two of the 15 sub-epidemics as "acute." The median estimates of R among the clusters ranged from 0.72 to 1.77. The majority of HIV lineages, 11 out of 15 clusters with 5+ members, appear to have been introduced to Mochudi between 1996 and 2002. The median peak duration of viral transmissions was 7.1 years (range 2.9-9.7 years). The median life span of identified HIV sub-epidemics, i.e., the time between the inferred epidemic origin and its most recent sample, was 13.1 years (range 10.2-22.1 years). Most viral transmissions within the sub-epidemics occurred between 1997 and 2007. The time period during which infected people are infectious appears to have decreased since the introduction of the national ART program in Botswana. Real-time HIV genotyping and breaking down local HIV epidemics into phylogenetically distinct sub-epidemics may help to reveal the structure and dynamics of HIV transmission networks in communities, and aid in the design of targeted interventions for members of the acute sub-epidemics that likely fuel local HIV/AIDS epidemics.

Co-infections and transmission networks of HCV, HIV-1 and HPgV among people who inject drugs

Co-infections with human immunodeficiency virus type 1 (HIV-1) and human pegivirus (HPgV) are common in hepatitis C virus (HCV)-infected individuals. However, analysis on the evolutionary dynamics and transmission network profiles of these viruses among individuals with multiple infections remains limited. A total of 228 injecting drug users (IDUs), either HCV- and/or HIV-1-infected, were recruited in Kuala Lumpur, Malaysia. HCV, HIV-1 and HPgV genes were sequenced, with epidemic growth rates assessed by the Bayesian coalescent method. Based on the sequence data, mono-, dual- and triple-infection were detected in 38.8%, 40.6% and 20.6% of the subjects, respectively. Fifteen transmission networks involving HCV (subtype 1a, 1b, 3a and 3b), HIV-1 (CRF33_01B) and HPgV (genotype 2) were identified and characterized. Genealogical estimates indicated that the predominant HCV, HIV-1 and HPgV genotypes were introduced into the IDUs population through multiple sub-epidemics that emerged as early as 1950s (HCV), 1980s (HIV-1) and 1990s (HPgV). By determining the difference in divergence times between viral lineages (ΔtMRCA), we also showed that the frequency of viral co-transmission is low among these IDUs. Despite increased access to therapy and other harm reduction interventions, the continuous emergence and coexistence of new transmission networks suggest persistent multiple viral transmissions among IDUs.

Combining the Estimated Date of HIV Infection with a Phylogenetic Cluster Study to Better Understand HIV Spread: Application in a Paris Neighbourhood

Objectives

To relate socio-demographic and virological information to phylogenetic clustering in HIV infected patients in a limited geographical area and to evaluate the role of recently infected individuals in the spread of HIV.

Methods

HIV-1 pol sequences from newly diagnosed and treatment-naive patients receiving follow-up between 2008 and 2011 by physicians belonging to a health network in Paris were used to build a phylogenetic tree using neighbour-joining analysis. Time since infection was estimated by immunoassay to define recently infected patients (very early infected presenters, VEP). Data on socio-demographic, clinical and biological features in clustered and non-clustered patients were compared. Chains of infection structure was also analysed.

Results

547 patients were included, 49 chains of infection containing 108 (20%) patients were identified by phylogenetic analysis. analysis. Eighty individuals formed pairs and 28 individuals were belonging to larger clusters. The median time between two successive HIV diagnoses in the same chain of infection was 248 days [CI = 176–320]. 34.7% of individuals were considered as VEP, and 27% of them were included in chains of infection. Multivariable analysis showed that belonging to a cluster was more frequent in VEP and those under 30 years old (OR: 3.65, 95 CI 1.49–8.95, p = 0.005 and OR: 2.42, 95% CI 1.05–5.85, p = 0.04 respectively). The prevalence of drug resistance was not associated with belonging to a pair or a cluster. Within chains, VEP were not grouped together more than chance predicted (p = 0.97).

Conclusions

Most newly diagnosed patients did not belong to a chain of infection, confirming the importance of undiagnosed or untreated HIV infected individuals in transmission. Furthermore, clusters involving both recently infected individuals and longstanding infected individuals support a substantial role in transmission of the latter before diagnosis.

Importance of Viral Sequence Length and Number of Variable and Informative Sites in Analysis of HIV Clustering

To improve the methodology of HIV cluster analysis, we addressed how analysis of HIV clustering is associated with parameters that can affect the outcome of viral clustering. The extent of HIV clustering and tree certainty was compared between 401 HIV-1C near full-length genome sequences and subgenomic regions retrieved from the LANL HIV Database. Sliding window analysis was based on 99 windows of 1,000 bp and 45 windows of 2,000 bp. Potential associations between the extent of HIV clustering and sequence length and the number of variable and informative sites were evaluated. The near full-length genome HIV sequences showed the highest extent of HIV clustering and the highest tree certainty. At the bootstrap threshold of 0.80 in maximum likelihood (ML) analysis, 58.9% of near full-length HIV-1C sequences but only 15.5% of partial pol sequences (ViroSeq) were found in clusters. Among HIV-1 structural genes, pol showed the highest extent of clustering (38.9% at a bootstrap threshold of 0.80), although it was significantly lower than in the near full-length genome sequences. The extent of HIV clustering was significantly higher for sliding windows of 2,000 bp than 1,000 bp. We found a strong association between the sequence length and proportion of HIV sequences in clusters, and a moderate association between the number of variable and informative sites and the proportion of HIV sequences in clusters. In HIV cluster analysis, the extent of detectable HIV clustering is directly associated with the length of viral sequences used, as well as the number of variable and informative sites. Near full-length genome sequences could provide the most informative HIV cluster analysis. Selected subgenomic regions with a high extent of HIV clustering and high tree certainty could also be considered as a second choice.

Emergence of clusters of CRF02_AG and B human immunodeficiency viral strains among men having sex with men exhibiting HIV primary infection in southeastern France

The number of new HIV diagnoses is increasing in the western world and transmission clusters have been recently identified among men having sex with men despite Highly Active Antiretroviral Therapy efficacy. The objective of this study was to assess temporal trends, epidemiological, clinical and virological characteristics of primary HIV infections. A retrospective analysis of 79 patients presenting primary HIV infections from 2005 to 2012 was performed in Marseille University Hospitals, southeastern France. Clinical, epidemiological and immunovirological data including phylogeny based on the polymerase gene were collected. 65 males and 14 females were enrolled. The main transmission route was homosexual contact (60.8%). Patients were mostly infected with subtype B (73.4%) and CRF02_AG (21.5%) HIV-1 strains. An increase in the annual number of HIV seroconversions among new HIV diagnoses from 5% in 2005 to 11.2% in 2012 (P = 0.06) and of the proportion of CRF02_AG HIV strains among primary HIV infections in 2011-2012 as compared to 2005-2010 (P = 0.055) was observed. Phylogenetic analysis revealed four transmission clusters including three transmission clusters among men having sex with men: two large clusters of nine CRF02_AG, six B HIV strains; and one small cluster of three B HIV strains. Clusters involved more frequently men (P = 0.01) belonging to caucasian ethicity (P = 0.05), with a higher HIV RNA load at inclusion (P = 0.03). These data highlight the importance of improving epidemiological surveillance and of implementing suitable prevention strategies to control the spread of HIV transmission among men having sex with men.

Transmission clustering among newly diagnosed HIV patients in Chicago, 2008 to 2011: using phylogenetics to expand knowledge of regional HIV transmission patterns

INTRODUCTION

HIV transmission cluster analyses can inform HIV prevention efforts. We describe the first such assessment for transmission clustering among HIV patients in Chicago.

METHODS

We performed transmission cluster analyses using HIV pol sequences from newly diagnosed patients presenting to Chicago's largest HIV clinic between 2008 and 2011. We compared sequences through progressive pairwise alignment, using neighbor joining to construct an unrooted phylogenetic tree. We defined clusters as >2 sequences among which each sequence had at least 1 partner within a genetic distance of ≤1.5%. We used multivariable regression to examine factors associated with clustering and used geospatial analysis to assess geographic proximity of phylogenetically clustered patients.

RESULTS

We compared sequences from 920 patients, median age of 35 years, 75% male, 67% black, 23% Hispanic, and 8% had a rapid plasma reagin titer ≥1:16 concurrent with their HIV diagnosis. We had HIV transmission risk data for 54%; 43% identified as men who have sex with men (MSM). Phylogenetic analysis demonstrated 123 patients (13%) grouped into 26 clusters, the largest having 20 members. In multivariable regression, age <25, black race, MSM status, male gender, higher HIV viral load, and rapid plasma reagin ≥1:16 associated with clustering. We did not observe geographic grouping of genetically clustered patients.

DISCUSSION

Our results demonstrate high rates of HIV transmission clustering, without local geographic foci, among young black MSM in Chicago. Applied prospectively, phylogenetic analyses could guide prevention efforts and help break the cycle of transmission.

Analysis of a local HIV-1 epidemic in portugal highlights established transmission of non-B and non-G subtypes

The existing data support Portugal as the western European country with the highest HIV-1 subtype diversity. However, detailed phylogenetic studies of Portuguese HIV-1 epidemics are still scarce. Thus, our main goal was to analyze the phylodynamics of a local HIV-1 infection in the Portuguese region of Minho. Molecular epidemiological analysis was applied to data from 289 HIV-1-infected individuals followed at the reference hospital of the province of Minho, Portugal, at which isolated viruses had been sequenced between 2000 and 2012. Viruses of the G (29.1%) and B (27.0%) subtypes were the most frequent, followed by recombinant forms (17.6%) and the C (14.5%), F1 (7.3%), and A1 (4.2%) subtypes. Multinomial logistic regression revealed that the odds of being infected with the A1 and F1 subtypes increased over the years compared with those with B, G, or C subtypes or recombinant viruses. As expected, polyphyletic patterns suggesting multiple and old introductions of the B and G subtypes were found. However, transmission clusters of non-B and non-G viruses among native individuals were also found, with the dates of the most recent common ancestor estimated to be in the early 2000s. Our study supports that the HIV-1 subtype diversity in the Portuguese region of Minho is high and has been increasing in a manner that is apparently driven by factors other than immigration and international travel. Infections with A1 and F1 viruses in the region of Minho are becoming established and are mainly found in sexually transmitted clusters, reinforcing the need for more efficacious control measures targeting this infection route.

Impact of sampling density on the extent of HIV clustering

Identifying and monitoring HIV clusters could be useful in tracking the leading edge of HIV transmission in epidemics. Currently, greater specificity in the definition of HIV clusters is needed to reduce confusion in the interpretation of HIV clustering results. We address sampling density as one of the key aspects of HIV cluster analysis. The proportion of viral sequences in clusters was estimated at sampling densities from 1.0% to 70%. A set of 1,248 HIV-1C env gp120 V1C5 sequences from a single community in Botswana was utilized in simulation studies. Matching numbers of HIV-1C V1C5 sequences from the LANL HIV Database were used as comparators. HIV clusters were identified by phylogenetic inference under bootstrapped maximum likelihood and pairwise distance cut-offs. Sampling density below 10% was associated with stochastic HIV clustering with broad confidence intervals. HIV clustering increased linearly at sampling density >10%, and was accompanied by narrowing confidence intervals. Patterns of HIV clustering were similar at bootstrap thresholds 0.7 to 1.0, but the extent of HIV clustering decreased with higher bootstrap thresholds. The origin of sampling (local concentrated vs. scattered global) had a substantial impact on HIV clustering at sampling densities ≥10%. Pairwise distances at 10% were estimated as a threshold for cluster analysis of HIV-1 V1C5 sequences. The node bootstrap support distribution provided additional evidence for 10% sampling density as the threshold for HIV cluster analysis. The detectability of HIV clusters is substantially affected by sampling density. A minimal genotyping density of 10% and sampling density of 50-70% are suggested for HIV-1 V1C5 cluster analysis.

Detection of transmission clusters of HIV-1 subtype C over a 21-year period in Cape Town, South Africa

Introduction

Despite recent breakthroughs in the fight against the HIV/AIDS epidemic within South Africa, the transmission of the virus continues at alarmingly high rates. It is possible, with the use of phylogenetic methods, to uncover transmission events of HIV amongst local communities in order to identify factors that may contribute to the sustained transmission of the virus. The aim of this study was to uncover transmission events of HIV amongst the infected population of Cape Town.

Methods and Results

We analysed gag p24 and RT-pol sequences which were generated from samples spanning over 21-years with advanced phylogenetic techniques. We identified two transmission clusters over a 21-year period amongst randomly sampled patients from Cape Town and the surrounding areas. We also estimated the origin of each of the identified transmission clusters with the oldest cluster dating back, on average, 30 years and the youngest dating back roughly 20 years.

Discussion and Conclusion

These transmission clusters represent the first identified transmission events among the heterosexual population in Cape Town. By increasing the number of randomly sampled specimens within a dataset over time, it is possible to start to uncover transmission events of HIV amongst local communities in generalized epidemics. This information can be used to produce targeted interventions to decrease transmission of HIV in Africa.

The impact of clinical, demographic and risk factors on rates of HIV transmission: a population-based phylogenetic analysis in British Columbia, Canada

BACKGROUND

The diversification of human immunodeficiency virus (HIV) is shaped by its transmission history. We therefore used a population based province wide HIV drug resistance database in British Columbia (BC), Canada, to evaluate the impact of clinical, demographic, and behavioral factors on rates of HIV transmission.

METHODS

We reconstructed molecular phylogenies from 27,296 anonymized bulk HIV pol sequences representing 7747 individuals in BC-about half the estimated HIV prevalence in BC. Infections were grouped into clusters based on phylogenetic distances, as a proxy for variation in transmission rates. Rates of cluster expansion were reconstructed from estimated dates of HIV seroconversion.

RESULTS

Our criteria grouped 4431 individuals into 744 clusters largely separated with respect to risk factors, including large established clusters predominated by injection drug users and more-recently emerging clusters comprising men who have sex with men. The mean log10 viral load of an individual's phylogenetic neighborhood (composed of 5 other individuals with shortest phylogenetic distances) increased their odds of appearing in a cluster by >2-fold per log10 viruses per milliliter.

CONCLUSIONS

Hotspots of ongoing HIV transmission can be characterized in near real time by the secondary analysis of HIV resistance genotypes, providing an important potential resource for targeting public health initiatives for HIV prevention.

Phylogenetic studies of transmission dynamics in generalized HIV epidemics: an essential tool where the burden is greatest?

Efficient and effective HIV prevention measures for generalized epidemics in sub-Saharan Africa have not yet been validated at the population level. Design and impact evaluation of such measures requires fine-scale understanding of local HIV transmission dynamics. The novel tools of HIV phylogenetics and molecular epidemiology may elucidate these transmission dynamics. Such methods have been incorporated into studies of concentrated HIV epidemics to identify proximate and determinant traits associated with ongoing transmission. However, applying similar phylogenetic analyses to generalized epidemics, including the design and evaluation of prevention trials, presents additional challenges. Here we review the scope of these methods and present examples of their use in concentrated epidemics in the context of prevention. Next, we describe the current uses for phylogenetics in generalized epidemics and discuss their promise for elucidating transmission patterns and informing prevention trials. Finally, we review logistic and technical challenges inherent to large-scale molecular epidemiological studies of generalized epidemics and suggest potential solutions.

Evolving molecular epidemiological profile of human immunodeficiency virus 1 in the southwest border of China

Background

We have previously reported in Xishuangbanna (Banna) Dai Autonomous Prefecture, a well-developed tourist destination in the southwest border of China, that HIV-1 transmitted dominantly through heterosexual contact with less divergent genotypes and few drug resistant mutations [1]. Due to the rapid increase of newly diagnosed HIV-1 cases per year in Banna in recent years, it’s important to evaluate the evolution of HIV-1 molecular epidemiology for the better understanding of ongoing HIV-1 outbreak in this region.

Methodology/Principal Findings

By sequencing of HIV-1 pol genes and phylogenetic analysis, we conducted a molecular epidemiologic study in 352 HIV-1-seropositive highly active antiretroviral treatment (HAART)-naïve individuals newly diagnosed at the Banna Center for Disease Control and Prevention between 2009 and 2011. Of 283 samples (84.1% taken from heterosexually acquired adults, 10.6% from needle-sharing drug users, 2.8% from men who have sex with men, 0.4% from children born from HIV-1-infected mothers, and 2.1% remained unknown) with successful sequencing for pol gene, we identified 108 (38.2%) HIV-1 subtype CRF08_BC, 101 (35.7%) CRF01_AE, 49 (17.3%) CRF07_BC, 5 (1.8%) C/CRF57_BC, 3 (1.1%) B’, 1 (0.4%) B/CRF51_01B, and 16 (5.7%) unique recombinants forms. Among these infected individuals, 104 (36.7%) cases showed drug resistant or resistance-relevant mutations, and 4 of them conferring high-level resistance to 3TC/FTC, EFV/NVP or NFV. Phylogenetic analysis revealed 21 clusters (2–7 sequences) with only 21.2% (60/283) sequences involved.

Conclusion/Significance

In contrast to our previous findings, CRF08_BC, replaced CRF01_AE, became the dominant genotype of HIV-1 in Banna prefecture. The viral strains with drug resistance mutations were detected frequently in newly diagnosed HIV-1-infected individuals in this region.

Trends and predictors of transmitted drug resistance (TDR) and clusters with TDR in a local Belgian HIV-1 epidemic

We aimed to study epidemic trends and predictors for transmitted drug resistance (TDR) in our region, its clinical impact and its association with transmission clusters. We included 778 patients from the AIDS Reference Center in Leuven (Belgium) diagnosed from 1998 to 2012. Resistance testing was performed using population-based sequencing and TDR was estimated using the WHO-2009 surveillance list. Phylogenetic analysis was performed using maximum likelihood and Bayesian techniques. The cohort was predominantly Belgian (58.4%), men who have sex with men (MSM) (42.8%), and chronically infected (86.5%). The overall TDR prevalence was 9.6% (95% confidence interval (CI): 7.7-11.9), 6.5% (CI: 5.0-8.5) for nucleoside reverse transcriptase inhibitors (NRTI), 2.2% (CI: 1.4-3.5) for non-NRTI (NNRTI), and 2.2% (CI: 1.4-3.5) for protease inhibitors. A significant parabolic trend of NNRTI-TDR was found (p = 0.019). Factors significantly associated with TDR in univariate analysis were male gender, Belgian origin, MSM, recent infection, transmission clusters and subtype B, while multivariate and Bayesian network analysis singled out subtype B as the most predictive factor of TDR. Subtype B was related with transmission clusters with TDR that included 42.6% of the TDR patients. Thanks to resistance testing, 83% of the patients with TDR who started therapy had undetectable viral load whereas half of the patients would likely have received a suboptimal therapy without this test. In conclusion, TDR remained stable and a NNRTI up-and-down trend was observed. While the presence of clusters with TDR is worrying, we could not identify an independent, non-sequence based predictor for TDR or transmission clusters with TDR that could help with guidelines or public health measures.

Discrepant amplification results during the development of an assay leads to reclassification of two AIDS reagent repository HIV-2 isolates as HIV-1

The development and verification of HIV-2 assays depends in part on the availability of well-characterized samples, including those from reagent repositories. During the development of an HIV-2 RNA quantification assay, two HIV-2 viral isolates (CDC 301340 and CDC 301342) obtained from the NIAID AIDS Reagent and Reference Repository were not detected leading to an investigation. Two HIV-2 primers/probe sets of known performance in real-time viral RNA quantification assays, targeting different regions of the virus, also failed to generate RT-PCR products for these two isolates. These isolates were tested in the HIV-1 specific COBAS AmpliPrep/COBAS TaqMan HIV-1 Test v2.0 (Roche Molecular Diagnostics) and were quantified at high copy number. Other HIV-2 isolates tested were not amplified in the COBAS HIV-1 TaqMan assay. Furthermore, the discrepant isolates were highly reactive in an HIV-1 p24 antigen test while the other HIV-2 isolates showed very weak, if any, cross-reactivity with the HIV-1 p24 assay. Phylogenetic tree analysis of sequences from the protease-reverse transcriptase regions of the discrepant HIV-2 isolates mapped with HIV-1 Group M, Subtype CRF02_AG confirming these isolates were of HIV-1 origin and had been misclassified as HIV-2. The use of misclassified isolates in the verification of molecular and immunological assays can lead to misinterpretation of test results, misdirection of efforts into assay redesign and increased development costs. The results of this study were shared with the NIAID AIDS Reagent Program, leading to the reclassification of the two discrepant isolates as HIV-1.

Phylogenetic and geospatial evaluation of HIV-1 subtype diversity at the largest HIV center in Rhode Island

Individuals infected with HIV-1 non-B subtypes are understudied in the United States. Their characterization may augment prevention and treatment interventions. We examined the regional molecular epidemiology of non-B subtypes using a combined phylogenetic and geospatial approach. HIV-1 pol sequences and clinical data obtained for routine clinical care were aggregated from 2004 to 2011 at the largest HIV center in Rhode Island. Subtyping was performed by neighbor-joining and maximum-likelihood phylogeny and compared across eight commonly used tools (HIVdb, REGA, RIP, NCBI, Geno2Pheno, EuResist, jpHMM and STAR) using proportional odds ordinal regression. Individuals with non-B subtypes were characterized according to demographics and risk factors for infection, intra-subtype clustering by maximum-likelihood phylogeny, and geospatial hotspot analysis using Getis-Ord Gi(∗) statistics. Of 1277 unique sequences, phylogenetic subtyping demonstrated 8.3% (N=106, 95% CI 6.8-10%) non-B subtypes and circulating recombinant forms (CRFs): CRF02_AG=46; A=15; C=15; CRF01_AE=6; CRF06_CPX=5; CRF14_BG=5; G=3; CRF43_02G=3; D=3; CRF24_BG=3; CRF11_CPX=1; F1=1. Compared to phylogeny, Geno2Pheno was the most concordant (86% exact match) followed by REGA (85%), EuResist (85%) and STAR (82%). Of 106 individuals with non-B subtypes, 50% were male, 71% acquired infection through heterosexual transmission; 76%, were born in Africa, 6% Southeast Asia, 5% the United States, 3% Central America, 1% Europe, and 9% unknown. Eighty percent of CRF02_AG, 93% of A and 87% of C sequences were from African-born individuals. Twenty-two percent of non-B subtypes formed transmission clusters, including a significant number of younger individuals with perinatally-acquired infection. Geospatial analyses revealed hotspots of B and non-B subtypes in the state capital with a more concentrated focus among non-B subtypes. Molecular examination of regional HIV diversity revealed a larger than expected non-subtype B infected population, mostly born in Africa, with low ongoing regional transmission. Phylogenetic and geospatial characterization of infection clusters is helpful to identify targets for treatment and prevention interventions.

Evidence that hepatitis C virus genome partly controls infection outcome

Infection by hepatitis C virus (HCV) leads to one of two outcomes; either the infection resolves within approximately 6 months or the virus can persist indefinitely. Host genetics are known to affect the likelihood of clearance or persistence. By contrast, the importance of the virus genotype in determining infection outcome is unknown, as quantifying this effect traditionally requires well-characterized transmission networks, which are rare. Extending phylogenetic approaches previously developed to estimate the virus control over set-point viral load in HIV-1 infections, we simulate inheritance of a binary trait along a phylogenetic tree, use this data to quantify how infection outcomes cluster and ascertain the effect of virus genotype on these. We apply our method to the Hepatitis C Incidence and Transmission Study in prisons (HITS-p) data set from Australia, as this cohort prospectively identified incident cases including viraemic subjects who ultimately clear the virus, thus providing us with a unique collection of sequences from clearing infections. We detect significant correlations between infection outcome and virus distance in the phylogeny for viruses of Genotype 1, with estimates lying at around 67%. No statistically significant estimates were obtained for viruses of Genotype 3a.

Toward an endgame: finding and engaging people unaware of their HIV-1 infection in treatment and prevention

Epidemic modeling suggests that a major scale-up in HIV treatment could have a dramatic impact on HIV incidence. This has led both researchers and policymakers to set a goal of an "AIDS-Free Generation." One of the greatest obstacles to achieving this objective is the number of people with undiagnosed HIV infection. Despite recent innovations, new research strategies are needed to identify, engage, and successfully treat people who are unaware of their infection.

HIV Type 1 transmission networks among men having sex with men and heterosexuals in Kenya

We performed a molecular phylogenetic study on HIV-1 polymerase sequences of men who have sex with men (MSM) and heterosexual patient samples in Kenya to characterize any observed HIV-1 transmission networks. HIV-1 polymerase sequences were obtained from samples in Nairobi and coastal Kenya from 84 MSM, 226 other men, and 364 women from 2005 to 2010. Using Bayesian phylogenetics, we tested whether sequences clustered by sexual orientation and geographic location. In addition, we used trait diffusion analyses to identify significant epidemiological links and to quantify the number of transmissions between risk groups. Finally, we compared 84 MSM sequences with all HIV-1 sequences available online at GenBank. Significant clustering of sequences from MSM at both coastal Kenya and Nairobi was found, with evidence of HIV-1 transmission between both locations. Although a transmission pair between a coastal MSM and woman was confirmed, no significant HIV-1 transmission was evident between MSM and the comparison population for the predominant subtype A (60%). However, a weak but significant link was evident when studying all subtypes together. GenBank comparison did not reveal other important transmission links. Our data suggest infrequent intermingling of MSM and heterosexual HIV-1 epidemics in Kenya.

Superinfection with drug-resistant HIV is rare and does not contribute substantially to therapy failure in a large European cohort

Background

Superinfection with drug resistant HIV strains could potentially contribute to compromised therapy in patients initially infected with drug-sensitive virus and receiving antiretroviral therapy. To investigate the importance of this potential route to drug resistance, we developed a bioinformatics pipeline to detect superinfection from routinely collected genotyping data, and assessed whether superinfection contributed to increased drug resistance in a large European cohort of viremic, drug treated patients.

Methods

We used sequence data from routine genotypic tests spanning the protease and partial reverse transcriptase regions in the Virolab and EuResist databases that collated data from five European countries. Superinfection was indicated when sequences of a patient failed to cluster together in phylogenetic trees constructed with selected sets of control sequences. A subset of the indicated cases was validated by re-sequencing pol and env regions from the original samples.

Results

4425 patients had at least two sequences in the database, with a total of 13816 distinct sequence entries (of which 86% belonged to subtype B). We identified 107 patients with phylogenetic evidence for superinfection. In 14 of these cases, we analyzed newly amplified sequences from the original samples for validation purposes: only 2 cases were verified as superinfections in the repeated analyses, the other 12 cases turned out to involve sample or sequence misidentification. Resistance to drugs used at the time of strain replacement did not change in these two patients. A third case could not be validated by re-sequencing, but was supported as superinfection by an intermediate sequence with high degenerate base pair count within the time frame of strain switching. Drug resistance increased in this single patient.

Conclusions

Routine genotyping data are informative for the detection of HIV superinfection; however, most cases of non-monophyletic clustering in patient phylogenies arise from sample or sequence mix-up rather than from superinfection, which emphasizes the importance of validation. Non-transient superinfection was rare in our mainly treatment experienced cohort, and we found a single case of possible transmitted drug resistance by this route. We therefore conclude that in our large cohort, superinfection with drug resistant HIV did not compromise the efficiency of antiretroviral treatment.

Automated analysis of phylogenetic clusters

BACKGROUND

As sequence data sets used for the investigation of pathogen transmission patterns increase in size, automated tools and standardized methods for cluster analysis have become necessary. We have developed an automated Cluster Picker which identifies monophyletic clades meeting user-input criteria for bootstrap support and maximum genetic distance within large phylogenetic trees. A second tool, the Cluster Matcher, automates the process of linking genetic data to epidemiological or clinical data, and matches clusters between runs of the Cluster Picker.

RESULTS

We explore the effect of different bootstrap and genetic distance thresholds on clusters identified in a data set of publicly available HIV sequences, and compare these results to those of a previously published tool for cluster identification. To demonstrate their utility, we then use the Cluster Picker and Cluster Matcher together to investigate how clusters in the data set changed over time. We find that clusters containing sequences from more than one UK location at the first time point (multiple origin) were significantly more likely to grow than those representing only a single location.

CONCLUSIONS

The Cluster Picker and Cluster Matcher can rapidly process phylogenetic trees containing tens of thousands of sequences. Together these tools will facilitate comparisons of pathogen transmission dynamics between studies and countries.

Treatment-naive individuals are the major source of transmitted HIV-1 drug resistance in men who have sex with men in the Swiss HIV Cohort Study

BACKGROUND

Human immunodeficiency virus type 1 (HIV-1) transmitted drug resistance (TDR) can compromise antiretroviral therapy (ART) and thus represents an important public health concern. Typically, sources of TDR remain unknown, but they can be characterized with molecular epidemiologic approaches. We used the highly representative Swiss HIV Cohort Study (SHCS) and linked drug resistance database (SHCS-DRDB) to analyze sources of TDR.

METHODS

ART-naive men who have sex with men with infection date estimates between 1996 and 2009 were chosen for surveillance of TDR in HIV-1 subtype B (N = 1674), as the SHCS-DRDB contains pre-ART genotypic resistance tests for >69% of this surveillance population. A phylogeny was inferred using pol sequences from surveillance patients and all subtype B sequences from the SHCS-DRDB (6934 additional patients). Potential sources of TDR were identified based on phylogenetic clustering, shared resistance mutations, genetic distance, and estimated infection dates.

RESULTS

One hundred forty of 1674 (8.4%) surveillance patients carried virus with TDR; 86 of 140 (61.4%) were assigned to clusters. Potential sources of TDR were found for 50 of 86 (58.1%) of these patients. ART-naive patients constitute 56 of 66 (84.8%) potential sources and were significantly overrepresented among sources (odds ratio, 6.43 [95% confidence interval, 3.22-12.82]; P < .001). Particularly large transmission clusters were observed for the L90M mutation, and the spread of L90M continued even after the near cessation of antiretroviral use selecting for that mutation. Three clusters showed evidence of reversion of K103N or T215Y/F.

CONCLUSIONS

Many individuals harboring viral TDR belonged to transmission clusters with other Swiss patients, indicating substantial domestic transmission of TDR in Switzerland. Most TDR in clusters could be linked to sources, indicating good surveillance of TDR in the SHCS-DRDB. Most TDR sources were ART naive. This, and the presence of long TDR transmission chains, suggests that resistance mutations are frequently transmitted among untreated individuals, highlighting the importance of early diagnosis and treatment.

The power and pitfalls of HIV phylogenetics in public health

Phylogenetics is the application of comparative studies of genetic sequences in order to infer evolutionary relationships among organisms. This tool can be used as a form of molecular epidemiology to enhance traditional population-level communicable disease surveillance. Phylogenetic study has resulted in new paradigms being created in the field of communicable diseases and this commentary aims to provide the reader with an explanation of how phylogenetics can be used in tracking infectious diseases. Special emphasis will be placed upon the application of phylogenetics as a tool to help elucidate HIV transmission patterns and the limitations to these methods when applied to forensic analysis. Understanding infectious disease epidemiology in order to prevent new transmissions is the sine qua non of public health. However, with increasing epidemiological resolution, there may be an associated potential loss of privacy to the individual. It is within this context that we aim to promote the discussion on how to use phylogenetics to achieve important public health goals, while at the same time protecting the rights of the individual.

Reconstructing the timing and dispersion routes of HIV-1 subtype B epidemics in the Caribbean and Central America: a phylogenetic story

The Caribbean and Central America are among the regions with highest HIV-1B prevalence worldwide. Despite of this high virus burden, little is known about the timing and the migration patterns of HIV-1B in these regions. Migration is one of the major processes shaping the genetic structure of virus populations. Thus, reconstruction of epidemiological network may contribute to understand HIV-1B evolution and reduce virus prevalence. We have investigated the spatio-temporal dynamics of the HIV-1B epidemic in The Caribbean and Central America using 1,610 HIV-1B partial pol sequences from 13 Caribbean and 5 Central American countries. Timing of HIV-1B introduction and virus evolutionary rates, as well as the spatial genetic structure of the HIV-1B populations and the virus migration patterns were inferred. Results revealed that in The Caribbean and Central America most of the HIV-1B variability was generated since the 80 s. At odds with previous data suggesting that Haiti was the origin of the epidemic in The Caribbean, our reconstruction indicated that the virus could have been disseminated from Puerto Rico and Antigua. These two countries connected two distinguishable migration areas corresponding to the (mainly Spanish-colonized) Easter and (mainly British-colonized) Western islands, which indicates that virus migration patterns are determined by geographical barriers and by the movement of human populations among culturally related countries. Similar factors shaped the migration of HIV-1B in Central America. The HIV-1B population was significantly structured according to the country of origin, and the genetic diversity in each country was associated with the virus prevalence in both regions, which suggests that virus populations evolve mainly through genetic drift. Thus, our work contributes to the understanding of HIV-1B evolution and dispersion pattern in the Americas, and its relationship with the geography of the area and the movements of human populations.

Future of phylogeny in HIV prevention

The success of the HIV Prevention Trials Network 052 trial has led to revisions in HIV-1 treatment guidelines. Antiretroviral therapy may reduce the risk of HIV-1 transmissions at the population level. The design of successful treatment as prevention interventions will be predicated on a comprehensive understanding of the spatial, temporal, and biological dynamics of heterosexual men who have sex with men and intravenous drug user epidemics. Viral phylogenetics can capture the underlying structure of transmission networks based on the genetic interrelatedness of viral sequences and cluster networks that could not be otherwise identified. This article describes the phylogenetic expansion of the Montreal men who have sex with men epidemic over the last decade. High rates of coclustering of primary infections are associated with 1 infection leading to 13 onward transmissions. Phylogeny substantiates the role of primary and recent stage infection in transmission dynamics, underlying the importance of timely diagnosis and immediate antiretroviral therapy initiation to avert transmission cascades.

HIV-1 clade B pol evolution following primary infection

OBJECTIVE

Characterize intra-individual HIV-1 subtype B pol evolution in antiretroviral naive individuals.

DESIGN

Longitudinal cohort study of individuals enrolled during primary infection.

METHODS

Eligible individuals were antiretroviral naïve participants enrolled in the cohort from December 1997-December 2005 and having at least two blood samples available with the first one collected within a year of their estimated date of infection. Population-based pol sequences were generated from collected blood samples and analyzed for genetic divergence over time in respect to dual infection status, HLA, CD4 count and viral load.

RESULTS

93 participants were observed for a median of 1.8 years (Mean = 2.2 years, SD =1.9 years). All participants classified as mono-infected had less than 0.7% divergence between any two of their pol sequences using the Tamura-Nei model (TN93), while individuals with dual infection had up to 7.0% divergence. The global substitution rates (substitutions/nucleotide/year) for mono and dually infected individuals were significantly different (p<0.001); however, substitution rates were not associated with HLA haplotype, CD4 or viral load.

CONCLUSIONS

Even after a maximum of almost 9 years of follow-up, all mono-infected participants had less than 1% divergence between baseline and longitudinal sequences, while participants with dual infection had 10 times greater divergence. These data support the use of HIV-1 pol sequence data to evaluate transmission events, networks and HIV-1 dual infection.

Social network-based recruitment successfully reveals HIV-1 transmission networks among high-risk individuals in El Salvador

OBJECTIVE

HIV in Central America is concentrated among certain groups such as men who have sex with men (MSM) and female sex workers (FSWs). We compared social recruitment chains and HIV transmission clusters from 699 MSM and 787 FSWs to better understand factors contributing to ongoing HIV transmission in El Salvador.

METHODS

Phylogenies were reconstructed using pol sequences from 119 HIV-positive individuals recruited by respondent-driven sampling (RDS) and compared with RDS chains in 3 cities in El Salvador. Transmission clusters with a mean pairwise genetic distance ≤ 0.015 and Bayesian posterior probabilities =1 were identified. Factors associated with cluster membership were evaluated among MSM.

RESULTS

Sequences from 34 (43%) MSM and 4 (10%) FSW grouped in 14 transmission clusters. Clusters were defined by risk group (12 MSM clusters) and geographic residence (only 1 spanned separate cities). In 4 MSM clusters (all n = 2), individuals were also members of the same RDS chain, but only 2 had members directly linked through recruitment. All large clusters (n ≥ 3) spanned >1 RDS chain. Among MSM, factors independently associated with cluster membership included recent infection by BED assay (P = 0.02), sex with stable male partners (P = 0.02), and sex with ≥ 3 male partners in the past year (P = 0.04).

CONCLUSIONS

We found few HIV transmissions corresponding directly with the social recruitment. However, we identified clustering in nearly one-half of MSM suggesting that RDS recruitment was indirectly but successfully uncovering transmission networks, particularly among recent infections. Interrogating RDS chains with phylogenetic analyses may help refine methods for identifying transmission clusters.

Molecular typing of the local HIV-1 epidemic in Serbia

Worldwide HIV-1 pandemic is becoming increasingly complex, with growing heterogeneity of subtypes and recombinant viruses. Previous studies have documented HIV-1 subtype B as the predominant one in Serbia, with limited presence and genetic diversity of non B subtypes. In recent years, MSM transmission has become the most frequently reported risk for HIV infection among newly diagnosed patients in Serbia, but very little is known of the network structure and dynamics of viral transmission in this and other risk groups. To gain insight about the HIV-1 subtypes distribution pattern as well as characteristics of HIV-1 transmission clusters in Serbia, we analyzed the genetic diversity of the pol gene segment in 221 HIV-1-infected patients sampled during 2002-2011. Subtype B was found to still be the most prevalent one in Serbia, accounting for over 90% of samples, while greater diversity of other subtypes was found than previously reported, including subtypes G, C, A, F, CRF01 and CRF02. In total, 41.3% of analyzed subtype B sequences were found associated in transmission clusters/network, that are highly related with MSM transmission route.

Application of a new informatics tool for contamination screening in the HIV sequencing laboratory

BACKGROUND

Current HIV-1 sequencing-based methods for detecting drug resistance-associated mutations are open and susceptible to contamination. Informatic identification of clinical sequences that are nearly identical to one another may indicate specimen-to-specimen contamination or another laboratory-associated issue.

OBJECTIVES

To design an informatic tool to rapidly identify potential contamination in the clinical laboratory using sequence analysis and to establish reference ranges for sequence variation in the HIV-1 protease and reverse transcriptase regions among a U.S. patient population.

STUDY DESIGN

We developed an open-source tool named HIV Contamination Detection (HIVCD). HIVCD was utilized to make pairwise comparisons of nearly 8000 partial HIV-1 pol gene sequences from patients across the United States and to calculate percent identities (PIDs) for each pair. ROC analysis and standard deviations of PID data were used to determine reference ranges for between-patient and within-patient comparisons and to guide selection of a threshold for identifying abnormally high PID between two unrelated sequences.

RESULTS

The PID reference range for between-patient comparisons ranged from 83.8 to 95.7% while within-patient comparisons ranged from 96 to 100%. Interestingly, 48% of between-patient sequence pairs with a PID>96.5 were geographically related. The selected threshold for abnormally high PIDs was 96 (AUC=0.993, sensitivity=0.980, specificity=0.999). During routine use, HIVCD identified a specimen mix-up and the source of contamination of a negative control.

CONCLUSIONS

In our experience, HIVCD is easily incorporated into laboratory workflow, useful for identifying potential laboratory errors, and contributes to quality testing. This type of analysis should be incorporated into routine laboratory practice.