Phylodynamics of HIV-1 from a phase-III AIDS vaccine trial in North America

Evaluation of haplotype callers for next-generation sequencing of viruses

Currently, the standard practice for assembling next-generation sequencing (NGS) reads of viral genomes is to summarize thousands of individual short reads into a single consensus sequence, thus confounding useful intra-host diversity information for molecular phylodynamic inference. It is hypothesized that a few viral strains may dominate the intra-host genetic diversity with a variety of lower frequency strains comprising the rest of the population. Several software tools currently exist to convert NGS sequence variants into haplotypes. Previous benchmarks of viral haplotype reconstruction programs used simulation scenarios that are useful from a mathematical perspective but do not reflect viral evolution and epidemiology. Here, we tested twelve NGS haplotype reconstruction methods using viral populations simulated under realistic evolutionary dynamics. We simulated coalescent-based populations that spanned known levels of viral genetic diversity, including mutation rates, sample size and effective population size, to test the limits of the haplotype reconstruction methods and to ensure coverage of predicted intra-host viral diversity levels (especially HIV-1). All twelve investigated haplotype callers showed variable performance and produced drastically different results that were mainly driven by differences in mutation rate and, to a lesser extent, in effective population size. Most methods were able to accurately reconstruct haplotypes when genetic diversity was low. However, under higher levels of diversity (e.g., those seen intra-host HIV-1 infections), haplotype reconstruction quality was highly variable and, on average, poor. All haplotype reconstruction tools, except QuasiRecomb and ShoRAH, greatly underestimated intra-host diversity and the true number of haplotypes. PredictHaplo outperformed, in regard to highest precision, recall, and lowest UniFrac distance values, the other haplotype reconstruction tools followed by CliqueSNV, which, given more computational time, may have outperformed PredictHaplo. Here, we present an extensive comparison of available viral haplotype reconstruction tools and provide insights for future improvements in haplotype reconstruction tools using both short-read and long-read technologies.

A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing

Washington, DC continues to experience a generalized HIV-1 epidemic. We characterized the local phylodynamics of HIV-1 in DC using next-generation sequencing (NGS) data. Viral samples from 68 participants from 2016 through 2017 were sequenced and paired with epidemiological data. Phylogenetic and network inferences, drug resistant mutations (DRMs), subtypes and HIV-1 diversity estimations were completed. Haplotypes were reconstructed to infer transmission clusters. Phylodynamic inferences based on the HIV-1 polymerase (pol) and envelope genes (env) were compared. Higher HIV-1 diversity (n.s.) was seen in men who have sex with men, heterosexual, and male participants in DC. 54.0% of the participants contained at least one DRM. The 40-49 year-olds showed the highest prevalence of DRMs (22.9%). Phylogenetic analysis of pol and env sequences grouped 31.9-33.8% of the participants into clusters. HIV-TRACE grouped 2.9-12.8% of participants when using consensus sequences and 9.0-64.2% when using haplotypes. NGS allowed us to characterize the local phylodynamics of HIV-1 in DC more broadly and accurately, given a better representation of its diversity and dynamics. Reconstructed haplotypes provided novel and deeper phylodynamic insights, which led to networks linking a higher number of participants. Our understanding of the HIV-1 epidemic was expanded with the powerful coupling of HIV-1 NGS data with epidemiological data.

A 28-Year History of HIV-1 Drug Resistance and Transmission in Washington, DC

Washington, DC consistently has one of the highest annual rates of new HIV-1 diagnoses in the United States over the last 10 years. To guide intervention and prevention strategies to combat DC HIV infection, it is helpful to understand HIV transmission dynamics in a historical context. Toward this aim, we conducted a retrospective study (years 1987-2015) of 3,349 HIV pol sequences (1,026 bp) from 1,996 individuals living in the DC area belonging to three different cohorts. We coupled HIV sequence data with clinical information (sex, risk factor, race/ethnicity, viral load, subtype, anti-retroviral regimen) to identify circulating drug resistant mutations (DRM) and transmission clusters and assess their persistence over time. Of the transmission clusters identified in the DC area, 78.0 and 31.7% involved MSM and heterosexuals, respectively. The longest spread of time for a single cluster was 5 years (2007-2012) using a distance-based network inference approach and 27 years (1987-2014) using a maximum likelihood phylogenetic approach. We found eight subtypes and nine recombinants. Genetic diversity increased steadily over time with a slight peak in 2009 and remained constant thereafter until 2015. Nucleotide diversity also increased over time while relative genetic diversity (BEAST) remained relatively steady over the last 28 years with slight increases since 2000 in subtypes B and C. Sequences from individuals on drug therapy contained the highest total number of DRMs (1,104-1,600) and unique DRMs (63-97) and the highest proportion (>20%) of resistant individuals. Heterosexuals (43.94%), MSM (40.13%), and unknown (44.26%) risk factors showed similar prevalence of DRMs, while injection drug users had a lower prevalence (33.33%). Finally, there was a 60% spike in the number of codons with DRMs between 2007 and 2010. Past patterns of HIV transmission and DRM accumulation over time described here will help to predict future efficacy of ART drugs based on DRMs persisting over time and identify risk groups of interest for prevention and intervention efforts within the DC population. Our results show how longitudinal data can help to understand the temporal dynamics of HIV-1 at the local level.

Virus-Like Particle, Liposome, and Polymeric Particle-Based Vaccines against HIV-1

It is acknowledged that vaccines remain the best hope for eliminating the HIV-1 epidemic. However, the failure to produce effective vaccine immunogens and the inability of conventional delivery strategies to elicit the desired immune responses remains a central theme and has ultimately led to a significant roadblock in HIV vaccine development. Consequently, significant efforts have been applied to generate novel vaccine antigens and delivery agents, which mimic viral structures for optimal immune induction. Here, we review the latest developments that have occurred in the nanoparticle vaccine field, with special emphasis on strategies that are being utilized to attain highly immunogenic, systemic, and mucosal anti-HIV humoral and cellular immune responses. This includes the design of novel immunogens, the central role of antigen-presenting cells, delivery routes, and biodistribution of nanoparticles to lymph nodes. In particular, we will focus on virus-like-particle formulations and their preclinical uses within the HIV prophylactic vaccine setting.

Characterization of HIV diversity, phylodynamics and drug resistance in Washington, DC

BACKGROUND

Washington DC has a high burden of HIV with a 2.0% HIV prevalence. The city is a national and international hub potentially containing a broad diversity of HIV variants; yet few sequences from DC are available on GenBank to assess the evolutionary history of HIV in the US capital. Towards this general goal, here we analyze extensive sequence data and investigate HIV diversity, phylodynamics, and drug resistant mutations (DRM) in DC.

METHODS

Molecular HIV-1 sequences were collected from participants infected through 2015 as part of the DC Cohort, a longitudinal observational study of HIV+ patients receiving care at 13 DC clinics. Sequences were paired with Cohort demographic, risk, and clinical data and analyzed using maximum likelihood, Bayesian and coalescent approaches of phylogenetic, network and population genetic inference. We analyzed 601 sequences from 223 participants for int (~864 bp) and 2,810 sequences from 1,659 participants for PR/RT (~1497 bp).

RESULTS

Ninety-nine and 94% of the int and PR/RT sequences, respectively, were identified as subtype B, with 14 non-B subtypes also detected. Phylodynamic analyses of US born infected individuals showed that HIV population size varied little over time with no significant decline in diversity. Phylogenetic analyses grouped 13.5% of the int sequences into 14 clusters of 2 or 3 sequences, and 39.0% of the PR/RT sequences into 203 clusters of 2-32 sequences. Network analyses grouped 3.6% of the int sequences into 4 clusters of 2 sequences, and 10.6% of the PR/RT sequences into 76 clusters of 2-7 sequences. All network clusters were detected in our phylogenetic analyses. Higher proportions of clustered sequences were found in zip codes where HIV prevalence is highest (r = 0.607; P<0.00001). We detected a high prevalence of DRM for both int (17.1%) and PR/RT (39.1%), but only 8 int and 12 PR/RT amino acids were identified as under adaptive selection. We observed a significant (P<0.0001) association between main risk factors (men who have sex with men and heterosexuals) and genotypes in the five well-supported clusters with sufficient sample size for testing.

DISCUSSION

Pairing molecular data with clinical and demographic data provided novel insights into HIV population dynamics in Washington, DC. Identification of populations and geographic locations where clustering occurs can inform and complement active surveillance efforts to interrupt HIV transmission.

Genetic and epidemiological insights into the emergence of peste des petits ruminants virus (PPRV) across Asia and Africa

Small ruminants are important components in the livelihood of millions of households in many parts of the world. The spread of the highly contagious peste des petits ruminants (PPR) disease, which is caused by an RNA virus, PPRV, across Asia and Africa remains a major concern. The present study explored the evolutionary and epidemiological dynamics of PPRV through the analyses of partial N-gene and F-gene sequences of the virus. All the four previously described PPRV lineages (I-IV) diverged from their common ancestor during the late-19(th) to early-20(th) century. Among the four lineages, PPRV-IV showed pronounced genetic structuring across the region; however, haplotype sharing among the geographic regions, together with the presence of multiple genetic clusters within a country, indicates the possibility of frequent mobility of the diseased individuals across the region. The gradual decline in the effective number of infections suggests a limited genetic variation, which could be attributed to the effective vaccination that has been practiced since 1990s. However, the movement of infected animals across the region likely contributes to the spread of PPRV-IV. No evidence of positive selection was identified from this study.

Phylogenetic relatedness of circulating HIV-1C variants in Mochudi, Botswana

BACKGROUND

Determining patterns of HIV transmission is increasingly important for the most efficient use of modern prevention interventions. HIV phylogeny can provide a better understanding of the mechanisms underlying HIV transmission networks in communities.

METHODS

To reconstruct the structure and dynamics of a local HIV/AIDS epidemic, the phylogenetic relatedness of HIV-1 subtype C env sequences obtained from 785 HIV-infected community residents in the northeastern sector of Mochudi, Botswana, during 2010-2013 was estimated. The genotyping coverage was estimated at 44%. Clusters were defined based on relatedness of HIV-1C env sequences and bootstrap support of splits.

RESULTS

The overall proportion of clustered HIV-1C env sequences was 19.1% (95% CI 17.5% to 20.8%). The proportion of clustered sequences from Mochudi was significantly higher than the proportion of non-Mochudi sequences that clustered, 27.0% vs. 14.7% (p = 5.8E-12; Fisher exact test). The majority of clustered Mochudi sequences (90.1%; 95% CI 85.1% to 93.6%) were found in the Mochudi-unique clusters. None of the sequences from Mochudi clustered with any of the 1,244 non-Botswana HIV-1C sequences. At least 83 distinct HIV-1C variants, or chains of HIV transmission, in Mochudi were enumerated, and their sequence signatures were reconstructed. Seven of 20 genotyped seroconverters were found in 7 distinct clusters.

CONCLUSIONS

The study provides essential characteristics of the HIV transmission network in a community in Botswana, suggests the importance of high sampling coverage, and highlights the need for broad HIV genotyping to determine the spread of community-unique and community-mixed viral variants circulating in local epidemics. The proposed methodology of cluster analysis enumerates circulating HIV variants and can work well for surveillance of HIV transmission networks. HIV genotyping at the community level can help to optimize and balance HIV prevention strategies in trials and combined intervention packages.

Sequences in glycoprotein gp41, the CD4 binding site, and the V2 domain regulate sensitivity and resistance of HIV-1 to broadly neutralizing antibodies

The swarm of quasispecies that evolves in each HIV-1-infected individual represents a source of closely related Env protein variants that can be used to explore various aspects of HIV-1 biology. In this study, we made use of these variants to identify mutations that confer sensitivity and resistance to the broadly neutralizing antibodies found in the sera of selected HIV-1-infected individuals. For these studies, libraries of Env proteins were cloned from infected subjects and screened for infectivity and neutralization sensitivity. The nucleotide sequences of the Env proteins were then compared for pairs of neutralization-sensitive and -resistant viruses. In vitro mutagenesis was used to identify the specific amino acids responsible for the neutralization phenotype. All of the mutations altering neutralization sensitivity/resistance appeared to induce conformational changes that simultaneously enhanced the exposure of two or more epitopes located in different regions of gp160. These mutations appeared to occur at unique positions required to maintain the quaternary structure of the gp160 trimer, as well as conformational masking of epitopes targeted by neutralizing antibodies. Our results show that sequences in gp41, the CD4 binding site, and the V2 domain all have the ability to act as global regulators of neutralization sensitivity. Our results also suggest that neutralization assays designed to support the development of vaccines and therapeutics targeting the HIV-1 Env protein should consider virus variation within individuals as well as virus variation between individuals.

Monoclonal antibodies to the V2 domain of MN-rgp120: fine mapping of epitopes and inhibition of α4β7 binding

BACKGROUND

Recombinant gp120 (MN-rgp120) was a major component of the AIDSVAX B/E vaccine used in the RV144 trial. This was the first clinical trial to show that vaccination could prevent HIV infection in humans. A recent RV144 correlates of protection study found that protection correlated with the presence of antibodies to the V2 domain. It has been proposed that antibodies to the α4β7 binding site in the V2 domain might prevent HIV-1 infection by blocking the ability of virions to recognize α4β7 on activated T-cells. In this study we investigated the specificity of monoclonal antibodies (MAbs) to the V2 domain of MN-rgp120 and examined the possibility that these antibodies could inhibit the binding of MN-rgp120 to the α4β7 integrin.

METHODOLOGY/PRINCIPAL FINDINGS

Nine MAbs to the V2 domain were isolated from mice immunized with recombinant envelope proteins. The ability of these MAbs to inhibit HIV infection, block the binding of gp120 to CD4, and block the binding of MN-rgp120 to the α4β7 integrin was measured. Mutational analysis showed that eight of the MAbs recognized two immunodominant clusters of amino acids (166-168 and 178-183) located at either end of the C strand within the four-strand anti-parallel sheet structure comprising the V1/V2 domain.

CONCLUSIONS/SIGNIFICANCE

These studies showed that the antigenic structure of the V2 domain is exceedingly complex and that MAbs isolated from mice immunized with MN-rgp120 exhibited a high level of strain specificity compared to MAbs to the V2 domain isolated from HIV-infected humans. We found that immunization with MN-rgp120 readily elicits antibodies to the V2 domain and some of these were able to block the binding of MN-rgp120 to the α4β7 integrin.

Genetic diversity and molecular epidemiology of HIV transmission

The high genetic diversity of HIV is one of its most significant features, as it has consequences in global distribution, vaccine design, therapy success, disease progression, transmissibility and viral load testing. Studying HIV diversity helps to understand its origins, migration patterns, current distribution and transmission events. New advances in sequencing technologies based on the parallel acquisition of data are now used to characterize within-host and population processes in depth. Additionally, we have seen similar advances in statistical methods designed to model the past history of lineages (the phylodynamic framework) to ultimately gain better insights into the evolutionary history of HIV. We can, for example, estimate population size changes, lineage dispersion over geographic areas and epidemiological parameters solely from sequence data. In this article, we review some of the evolutionary approaches used to study transmission patterns and processes in HIV and the insights gained from such studies.

The evolution of HIV: inferences using phylogenetics

Molecular phylogenetics has revolutionized the study of not only evolution but also disparate fields such as genomics, bioinformatics, epidemiology, ecology, microbiology, molecular biology and biochemistry. Particularly significant are its achievements in population genetics as a result of the development of coalescent theory, which have contributed to more accurate model-based parameter estimation and explicit hypothesis testing. The study of the evolution of many microorganisms, and HIV in particular, have benefited from these new methodologies. HIV is well suited for such sophisticated population analyses because of its large population sizes, short generation times, high substitution rates and relatively small genomes. All these factors make HIV an ideal and fascinating model to study molecular evolution in real time. Here we review the significant advances made in HIV evolution through the application of phylogenetic approaches. We first examine the relative roles of mutation and recombination on the molecular evolution of HIV and its adaptive response to drug therapy and tissue allocation. We then review some of the fundamental questions in HIV evolution in relation to its origin and diversification and describe some of the insights gained using phylogenies. Finally, we show how phylogenetic analysis has advanced our knowledge of HIV dynamics (i.e., phylodynamics).

Phylodynamics of HIV-1 from a phase III AIDS vaccine trial in Bangkok, Thailand

BACKGROUND

In 2003, a phase III placebo-controlled trial (VAX003) was completed in Bangkok, Thailand. Of the 2,546 individuals enrolled in the trial based on high risk for infection through injection drug use (IDU), we obtained clinical samples and HIV-1 sequence data (envelope glycoprotein gene gp120) from 215 individuals who became infected during the trial. Here, we used these data in combination with other publicly available gp120 sequences to perform a molecular surveillance and phylodynamic analysis of HIV-1 in Thailand.

METHODOLOGY AND FINDINGS

Phylogenetic and population genetic estimators were used to assess HIV-1 gp120 diversity as a function of vaccination treatment, viral load (VL) and CD4(+) counts, to identify transmission clusters and to investigate the timescale and demographics of HIV-1 in Thailand. Three HIV-1 subtypes were identified: CRF01_AE (85% of the infections), subtype B (13%) and CRF15_AE (2%). The Bangkok IDU cohort showed more gp120 diversity than other Asian IDU cohorts and similar diversity to that observed in sexually infected individuals. Moreover, significant differences (P<0.02) in genetic diversity were observed in CRF01_AE IDU with different VL and CD4(+) counts. No phylogenetic structure was detected regarding any of the epidemiological and clinical factors tested, although high proportions (35% to 50%) of early infections fell into clusters, which suggests that transmission chains associated with acute infection play a key role on HIV-1 spread among IDU. CRF01_AE was estimated to have emerged in Thailand in 1984.5 (1983-1986), 3-6 years before the first recognition of symptomatic patients (1989). The relative genetic diversity of the HIV-1 population has remained high despite decreasing prevalence rates since the mid 1990s.

CONCLUSIONS

Our study and recent epidemiological reports indicate that HIV-1 is still a major threat in Thailand and suggest that HIV awareness and prevention needs to be strengthened to avoid AIDS resurgence.

Magnitude and breadth of a nonprotective neutralizing antibody response in an efficacy trial of a candidate HIV-1 gp120 vaccine

BACKGROUND

A candidate vaccine consisting of human immunodeficiency virus type 1 (HIV-1) subunit gp120 protein was found previously to be nonprotective in an efficacy trial (Vax004) despite strong antibody responses against the vaccine antigens. Here we assessed the magnitude and breadth of neutralizing antibody responses in Vax004.

METHODS

Neutralizing antibodies were measured against highly sensitive (tier 1) and moderately sensitive (tier 2) strains of HIV-1 subtype B in 2 independent assays. Vaccine recipients were stratified by sex, race, and high versus low behavioral risk of HIV-1 acquisition.

RESULTS

Most vaccine recipients mounted potent neutralizing antibody responses against HIV-1(MN) and other tier 1 viruses. Occasional weak neutralizing activity was detected against tier 2 viruses. The response against tier 1 and tier 2 viruses was significantly stronger in women than in men. Race and behavioral risk of HIV-1 acquisition had no significant effect on the response. Prior vaccination had little effect on the neutralizing antibody response that arose after infection.

CONCLUSIONS

Weak overall neutralizing antibody responses against tier 2 viruses is consistent with a lack of protection in this trial. The magnitude and breadth of neutralization reported here should be useful for identifying improved vaccines.

Viral phylodynamics and the search for an 'effective number of infections'

Information on the dynamics of the effective population size over time can be obtained from the analysis of phylogenies, through the application of time-varying coalescent models. This approach has been used to study the dynamics of many different viruses, and has demonstrated a wide variety of patterns, which have been interpreted in the context of changes over time in the ‘effective number of infections’, a quantity proportional to the number of infected individuals. However, for infectious diseases, the rate of coalescence is driven primarily by new transmissions i.e. the incidence, and only indirectly by the number of infected individuals through sampling effects. Using commonly used epidemiological models, we show that the coalescence rate may indeed reflect the number of infected individuals during the initial phase of exponential growth when time is scaled by infectivity, but in general, a single change in time scale cannot be used to estimate the number of infected individuals. This has important implications when integrating phylogenetic data in the context of other epidemiological data.

Phase III HIV vaccine trial in Thailand: a step toward a protective vaccine for HIV

The large human efficacy trail in Thailand, RV144, was concluded in the summer of 2009. This is the first Phase III trial to show limited, but significant, efficacy in preventing HIV acquisition. This trial represents the first sign that a preventive vaccine for HIV may be feasible. The vaccine regimen tested in Thailand consisted of priming with a Canarypox vector carrying three synthetic HIV genes. The priming was followed by booster inoculations with two recombinant envelope proteins from HIV, clade B and E. The need to understand the role in protection from HIV acquisition of the new responses, induced by this vaccine combination, has brought together many researchers with the common goal of improving the development of a safe and effective vaccine for HIV.

Mutation at a single position in the V2 domain of the HIV-1 envelope protein confers neutralization sensitivity to a highly neutralization-resistant virus

Understanding the determinants of neutralization sensitivity and resistance is important for the development of an effective human immunodeficiency virus type 1 (HIV-1) vaccine. In these studies, we have made use of the swarm of closely related envelope protein variants (quasispecies) from an extremely neutralization-resistant clinical isolate in order to identify mutations that conferred neutralization sensitivity to antibodies in sera from HIV-1-infected individuals. Here, we describe a virus with a rare mutation at position 179 in the V2 domain of gp120, where replacement of aspartic acid (D) by asparagine (N) converts a virus that is highly resistant to neutralization by multiple polyclonal and monoclonal antibodies, as well as antiviral entry inhibitors, to one that is sensitive to neutralization. Although the V2 domain sequence is highly variable, D at position 179 is highly conserved in HIV-1 and simian immunodeficiency virus (SIV) and is located within the LDI/V recognition motif of the recently described α4β7 receptor binding site. Our results suggest that the D179N mutation induces a conformational change that exposes epitopes in both the gp120 and the gp41 portions of the envelope protein, such as the CD4 binding site and the MPER, that are normally concealed by conformational masking. Our results suggest that D179 plays a central role in maintaining the conformation and infectivity of HIV-1 as well as mediating binding to α4β7.