Development of a panel of well-characterized human immunodeficiency virus type 1 isolates from newly diagnosed patients including acute and recent infections

High prevalence and diversity of HIV-1 non-B genetic forms due to immigration in southern Spain: A phylogeographic approach

Phylogenetic studies are a valuable tool to understand viral transmission patterns and the role of immigration in HIV-1 spread. We analyzed the spatio-temporal relationship of different HIV-1 non-B subtype variants over time using phylogenetic analysis techniques. We collected 693 pol (PR+RT) sequences that were sampled from 2005 to 2012 from naïve patients in different hospitals in southern Spain. We used REGA v3.0 to classify them into subtypes and recombinant forms, which were confirmed by phylogenetic analysis through maximum likelihood (ML) using RAxML. For the main HIV-1 non-B variants, publicly available, genetically similar sequences were sought using HIV-BLAST. The presence of HIV-1 lineages circulating in our study population was established using ML and Bayesian inference (BEAST v1.7.5) and transmission networks were identified. We detected 165 (23.4%) patients infected with HIV-1 non-B variants: 104 (63%) with recombinant viruses in pol: CRF02_AG (71, 43%), CRF14_BG (8, 4.8%), CRF06_cpx (5, 3%) and nine other recombinant forms (11, 6.7%) and unique recombinants (9, 5.5%). The rest (61, 37%) were infected with non-recombinant subtypes: A1 (30, 18.2%), C (7, [4.2%]), D (3, [1.8%]), F1 (9, 5.5%) and G (12, 7.3%). Most patients infected with HIV-1 non-B variants were men (63%, p < 0.001) aged over 35 (73.5%, p < 0.001), heterosexuals (92.2%, p < 0.001), from Africa (59.5%, p < 0.001) and living in the El Ejido area (62.4%, p<0.001). We found lineages of epidemiological relevance (mainly within Subtype A1), imported primarily through female sex workers from East Europe. We detected 11 transmission clusters of HIV-1 non-B Subtypes, which included patients born in Spain in half of them. We present the phylogenetic profiles of the HIV-1 non-B variants detected in southern Spain, and explore their putative geographical origins. Our data reveals a high HIV-1 genetic diversity likely due to the import of viral lineages that circulate in other countries. The highly immigrated El Ejido area acts as a gateway through which different subtypes are introduced into other regions, hence the importance of setting up epidemiological control measures to prevent future outbreaks.

Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004-2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations

Our objectives were to carry out an epidemiological surveillance study on transmitted drug resistance (TDR) among individuals newly diagnosed of HIV-1 infection during a nine year period in Spain and to assess the role of transmission clusters (TC) in the propagation of resistant strains. An overall of 1614 newly diagnosed individuals were included in the study from January 2004 through December 2012. Individuals come from two different Spanish regions: Galicia and the Basque Country. Resistance mutations to reverse transcriptase inhibitors (RTI) and protease inhibitors (PI) were analyzed according to mutations included in the surveillance drug-resistance mutations list updated in 2009. TC were defined as those comprising viruses from five or more individuals whose sequences clustered in maximum likelihood phylogenetic trees with a bootstrap value ≥90%. The overall prevalence of TDR to any drug was 9.9%: 4.9% to nucleoside RTIs (NRTIs), 3.6% to non-nucleoside RTIs (NNRTIs), and 2.7% to PIs. A significant decrease of TDR to NRTIs over time was observed [from 10% in 2004 to 2% in 2012 (p=0.01)]. Sixty eight (42.2%) of 161 sequences with TDR were included in 25 TC composed of 5 or more individuals. Of them, 9 clusters harbored TDR associated with high level resistance to antiretroviral drugs. T215D revertant mutation was transmitted in a large cluster comprising 25 individuals. The impact of epidemiological networks on TDR frequency may explain its persistence in newly diagnosed individuals. The knowledge of the populations involved in TC would facilitate the design of prevention programs and public health interventions.

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.

An integrated map of HIV genome-wide variation from a population perspective

Background

The HIV pandemic is characterized by extensive genetic variability, which has challenged the development of HIV drugs and vaccines. Although HIV genomes have been classified into different types, groups, subtypes and recombinants, a comprehensive study that maps HIV genome-wide diversity at the population level is still lacking to date. This study aims to characterize HIV genomic diversity in large-scale sequence populations, and to identify driving factors that shape HIV genome diversity.

Results

A total of 2996 full-length genomic sequences from 1705 patients infected with 16 major HIV groups, subtypes and circulating recombinant forms (CRFs) were analyzed along with structural, immunological and peptide inhibitor information. Average nucleotide diversity of HIV genomes was almost 50% between HIV-1 and HIV-2 types, 37.5% between HIV-1 groups, 14.7% between HIV-1 subtypes, 8.2% within individual HIV-1 subtypes and less than 1% within single patients. Along the HIV genome, diversity patterns and compositions of nucleotides and amino acids were highly similar across different groups, subtypes and CRFs. Current HIV-derived peptide inhibitors were predominantly derived from conserved, solvent accessible and intrinsically ordered structures in the HIV-1 subtype B genome. We identified these conserved regions in Capsid, Nucleocapsid, Protease, Integrase, Reverse transcriptase, Vpr and the GP41 N terminus as potential drug targets. In the analysis of factors that impact HIV-1 genomic diversity, we focused on protein multimerization, immunological constraints and HIV-human protein interactions. We found that amino acid diversity in monomeric proteins was higher than in multimeric proteins, and diversified positions were preferably located within human CD4 T cell and antibody epitopes. Moreover, intrinsic disorder regions in HIV-1 proteins coincided with high levels of amino acid diversity, facilitating a large number of interactions between HIV-1 and human proteins.

Conclusions

This first large-scale analysis provided a detailed mapping of HIV genomic diversity and highlighted drug-target regions conserved across different groups, subtypes and CRFs. Our findings suggest that, in addition to the impact of protein multimerization and immune selective pressure on HIV-1 diversity, HIV-human protein interactions are facilitated by high variability within intrinsically disordered structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0148-6) contains supplementary material, which is available to authorized users.

Identification of new and unusual rev and nef transcripts expressed by an HIV type 1 primary isolate

We analyzed RNA splice site usage in three HIV-1 subtype B primary isolates through reverse transcriptase polymerase chain reaction (RT-PCR) amplification of spliced RNAs using a fluorescently labeled primer, with computerized size determination and quantification of PCR products, which were also identified by clone sequencing. In one isolate, P2149-3, unusual and unreported spliced transcripts were detected. This isolate preferentially used for rev RNA generation a 3' splice site (3'ss) located five nucleotides upstream of A4a, previously identified only in a T cell line-adapted virus and in a group O isolate, and designated A4d. P2149-3 also used an unreported 3'ss for rev RNA generation, designated A4h, located 20 nucleotides upstream of 3'ss A4c. Additionally, unusual nef RNAs using 3'ss A5a and A7a and with exon composition 1.3.7 were identified. The identification of several unusual and unreported spliced transcripts in an HIV-1 primary isolate suggests a greater diversity of splice site usage in HIV-1 than previously appreciated.

Improving HIV coreceptor usage prediction in the clinic using hints from next-generation sequencing data

MOTIVATION

Due to the high mutation rate of human immunodeficiency virus (HIV), drug-resistant-variants emerge frequently. Therefore, researchers are constantly searching for new ways to attack the virus. One new class of anti-HIV drugs is the class of coreceptor antagonists that block cell entry by occupying a coreceptor on CD4 cells. This type of drug just has an effect on the subset of HIVs that use the inhibited coreceptor. A good prediction of whether the viral population inside a patient is susceptible to the treatment is hence very important for therapy decisions and pre-requisite to administering the respective drug. The first prediction models were based on data from Sanger sequencing of the V3 loop of HIV. Recently, a method based on next-generation sequencing (NGS) data was introduced that predicts labels for each read separately and decides on the patient label through a percentage threshold for the resistant viral minority.

RESULTS

We model the prediction problem on the patient level taking the information of all reads from NGS data jointly into account. This enables us to improve prediction performance for NGS data, but we can also use the trained model to improve predictions based on Sanger sequencing data. Therefore, also laboratories without NGS capabilities can benefit from the improvements. Furthermore, we show which amino acids at which position are important for prediction success, giving clues on how the interaction mechanism between the V3 loop and the particular coreceptors might be influenced.

AVAILABILITY

A webserver is available at http://coreceptor.bioinf.mpi-inf.mpg.de.

CONTACT

nico.pfeifer@mpi-inf.mpg.de.

Evaluation of heterologous vaginal SHIV SF162p4 infection following vaccination with a polyvalent Clade B virus-like particle vaccine

The vast diversity of HIV-1 infections has greatly impeded the development of a successful HIV-1/AIDS vaccine. Previous vaccine work has demonstrated limited levels of protection against SHIV/SIV infection, but protection was observed only when the challenge virus was directly matched to the vaccine strain. As it is likely impossible to directly match the vaccine strain to all infecting strains in nature, it is necessary to develop an HIV-1 vaccine that can protect against a heterologous viral challenge. In this study we investigated the ability of polyvalent and consensus vaccines to protect against a heterologous clade B challenge. Rhesus macaques were vaccinated with ConB or PolyB virus-like particle vaccines. All vaccines were highly immunogenic with high titers of antibody found in all vaccinated groups against SIV Gag. Antibody responses were also observed against a diverse panel of clade B envelopes. Following vaccination nonhuman primates (NHPs) were challenged via the vaginal route with SHIV(SF162p4). The PolyB vaccine induced a 66.7% reduction in the rate of infection as well as causing a two log reduction in viral burden if infection was not blocked. ConB vaccination had no effect on either the infection rate or viral burden. These results indicate that a polyvalent clade-matched vaccine is better able to protect against a heterologous challenge as compared to a consensus vaccine.

Identification of new splice sites used for generation of rev transcripts in human immunodeficiency virus type 1 subtype C primary isolates

The HIV-1 primary transcript undergoes a complex splicing process by which more than 40 different spliced RNAs are generated. One of the factors contributing to HIV-1 splicing complexity is the multiplicity of 3′ splice sites (3'ss) used for generation of rev RNAs, with two 3'ss, A4a and A4b, being most commonly used, a third site, A4c, used less frequently, and two additional sites, A4d and A4e, reported in only two and one isolates, respectively. HIV-1 splicing has been analyzed mostly in subtype B isolates, and data on other group M clades are lacking. Here we examine splice site usage in three primary isolates of subtype C, the most prevalent clade in the HIV-1 pandemic, by using an in vitro infection assay of peripheral blood mononuclear cells. Viral spliced RNAs were identified by RT-PCR amplification using a fluorescently-labeled primer and software analyses and by cloning and sequencing the amplified products. The results revealed that splice site usage for generation of rev transcripts in subtype C differs from that reported for subtype B, with most rev RNAs using two previously unreported 3'ss, one located 7 nucleotides upstream of 3'ss A4a, designated A4f, preferentially used by two isolates, and another located 14 nucleotides upstream of 3'ss A4c, designated A4g, preferentially used by the third isolate. A new 5′ splice site, designated D2a, was also identified in one virus. Usage of the newly identified splice sites is consistent with sequence features commonly found in subtype C viruses. These results show that splice site usage may differ between HIV-1 subtypes.

Evaluation of genotypic tropism prediction tests compared with in vitro co-receptor usage in HIV-1 primary isolates of diverse subtypes

OBJECTIVES

To evaluate the sensitivity and specificity of genotypic methods for predicting the co-receptor usage of subtypes B and non-B HIV-1 primary isolates, using as gold standard the infectivity of each primary isolate in GHOST cells stably expressing HIV-1 co-receptors.

METHODS

Primary isolates were obtained by co-culturing either patient's peripheral blood mononuclear cells (PBMCs) or ultracentrifuged plasma with donor-activated PBMCs. In vitro co-receptor usage was determined by infecting GHOST cells. Tropism prediction, based on V3 sequences, was determined with simple rules and bioinformatic tools (Geno2pheno[coreceptor] and WebPSSM).

RESULTS

This study includes 102 HIV-1 primary isolates; 23 (22.5%) subtype B and 79 (77.5%) non-B genetic forms. V3 sequences were classified into six subtypes (A-G), although 32 (31.4%) were circulating recombinant forms and 21 (20.6%) were unique recombinant forms. Sixty-nine isolates were R5, 27 R5X4 and 6 X4. The highest levels of sensitivity and specificity for the detection of X4 strains among V3 sequences, between 91% and 100%, were obtained by using PSSM(x4r5), PSSM(si/nsi) and the 11/25 rule for sequences of subtypes A, B and G, but not for subtype F. Establishing the recommended cut-off for clinical settings of a 10% false positive rate for Geno2pheno, we obtained 93% specificity and 97% sensitivity.

CONCLUSIONS

Comparing genotypic assays for HIV-1 co-receptor use with a cell-culture phenotypic assay could provide more reliable results of sensitivity and specificity for the detection of X4 strains than comparing them with recombinant assays, considered as gold standard. In general, except for subtype F isolates, there is a good correlation for tropism prediction.

Characterization of partial and near full-length genomes of HIV-1 strains sampled from recently infected individuals in São Paulo, Brazil

Background

Genetic variability is a major feature of human immunodeficiency virus type 1 (HIV-1) and is considered the key factor frustrating efforts to halt the HIV epidemic. A proper understanding of HIV-1 genomic diversity is a fundamental prerequisite for proper epidemiology, genetic diagnosis, and successful drugs and vaccines design. Here, we report on the partial and near full-length genomic (NFLG) variability of HIV-1 isolates from a well-characterized cohort of recently infected patients in São Paul, Brazil.

Methodology

HIV-1 proviral DNA was extracted from the peripheral blood mononuclear cells of 113 participants. The NFLG and partial fragments were determined by overlapping nested PCR and direct sequencing. The data were phylogenetically analyzed.

Results

Of the 113 samples (90.3% male; median age 31 years; 79.6% homosexual men) studied, 77 (68.1%) NFLGs and 32 (29.3%) partial fragments were successfully subtyped. Of the successfully subtyped sequences, 88 (80.7%) were subtype B sequences, 12 (11%) BF1 recombinants, 3 (2.8%) subtype C sequences, 2 (1.8%) BC recombinants and subclade F1 each, 1 (0.9%) CRF02 AG, and 1 (0.9%) CRF31 BC. Primary drug resistance mutations were observed in 14/101 (13.9%) of samples, with 5.9% being resistant to protease inhibitors and nucleoside reverse transcriptase inhibitors (NRTI) and 4.9% resistant to non-NRTIs. Predictions of viral tropism were determined for 86 individuals. X4 or X4 dual or mixed-tropic viruses (X4/DM) were seen in 26 (30.2%) of subjects. The proportion of X4 viruses in homosexuals was detected in 19/69 (27.5%).

Conclusions

Our results confirm the existence of various HIV-1 subtypes circulating in São Paulo, and indicate that subtype B account for the majority of infections. Antiretroviral (ARV) drug resistance is relatively common among recently infected patients. The proportion of X4 viruses in homosexuals was significantly higher than the proportion seen in other study populations.

Construction and phenotypic characterization of HIV type 1 functional envelope clones of subtypes G and F

Subtype G has been estimated to represent the fourth most prevalent clade in the HIV-1 pandemic and subtype F is widely circulating in parts of South America (frequently within BF recombinant forms) and in Romania. However, functional envelope clones of these subtypes are lacking, which are needed for studies on antibody-mediated neutralization, coreceptor usage, and efficiency of viral entry inhibitor drugs. Here we report the construction, neutralization properties, and coreceptor usage of HIV-1 functional envelope clones of subtypes G (n = 15) and F (n = 7). These clones were obtained through RT-PCR amplification of HIV-1 gp160 from plasma RNA, and were used for pseudovirus production. All 15 subtype G-enveloped pseudoviruses were resistant to neutralization by gp120-targeted broadly neutralizing monoclonal antibodies (MAbs) b12 and 2G12, while a majority were neutralized by gp41-targeted MAbs 2F5 and 4E10. With regard to the subtype F envelopes, all seven pseudoviruses were resistant to 2F5 and b12, six were resistant to G12, and six were neutralized by 4E10. Coreceptor usage testing revealed that 21 of 22 envelopes were CCR5-tropic, including all 15 subtype G envelopes, seven of which were from patients with CD4(+) T cell counts <200/ml. These results confirm the broadly neutralizing activity of 4E10 on envelope clones across all tested group M clades, including subtypes G and F, reveal the resistance of most subtype F-enveloped pseudoviruses to broadly neutralizing MAbs b12, 2G12, and 2F5, and suggest that, similarly to subtype C, CXCR4 tropism is uncommon in subtype G, even at advanced stages of infection.

Short communication: Biological and genetic characterization of HIV type 1 subtype B and nonsubtype B transmitted viruses: usefulness for vaccine candidate assessment

Due to the extraordinary degree of genetic diversity of HIV-1 and the structural complexity of its envelope glycoproteins, designing an effective vaccine is difficult, requiring the development of viral reagents to assess vaccine-elicited neutralizing antibodies. The aim of this study was to improve on our previously developed panel of HIV-1 strains of different genetic forms, focusing on strains from acute and recently acquired infections as the most representative of the transmitted viruses. HIV-1 primary isolates were expanded in peripheral blood mononuclear cells. Viral stocks of 40 ml each were produced. Syncytium-inducing (SI) phenotype, coreceptor use, and TCID(50)/ml were determined. Near full-length HIV-1 genomes were amplified by RT-nested PCR in four overlapping segments. Phylogenetic analyses were performed with neighbor-joining trees and bootscanning. Forty-four HIV-1 strains were included in the panel. Twenty-four (54.1%) strains were from early infections (16 acute and 8 recent); of them, 21 (87%) were sexually transmitted. NSI/R5 phenotype was detected in 37 (84.1%) viruses and SI/R5,X4 in another 7 (15.9%). TCID(50)/ml ranged between 10(4) and 10(6.6). Twelve different genetic forms constituted this panel: subtypes A1, B, C, F1, and G; circulating recombinant forms CRF02_AG, CRF14_BG, and CRF24_BG; and unique recombinant forms CRF02_AG/A3, BF1, CRF12_BF/B, and DF1G. In conclusion, in this study, we report the development of a comprehensive and well-characterized panel of HIV-1 isolates for assessing neutralization in HIV vaccine research. This panel is available for distribution through the Programme EVA Centre for AIDS Reagents, National Institute for Biological Standard and Control (NIBSC).

Identification of a new HIV type 1 circulating BF intersubtype recombinant form (CRF47_BF) in Spain

We report the identification of a new HIV-1 circulating recombinant form (CRF47_BF) derived from subtypes B and F. It was initially identified in protease-reverse transcriptase sequences from nine individuals from three separate regions of Spain who acquired HIV-1 infection via sexual contact. All nine sequences formed a strongly supported phylogenetic cluster, branching apart from all known CRFs, and in bootscan analyses were BF mosaics with two coincident breakpoints. Two epidemiologically unlinked viruses were sequenced in near full-length genomes, which exhibited identical mosaic structures, with 16 intersubtype breakpoints in a genome predominantly of subtype B. Subtype F segments of the new CRF failed to cluster with any of the near full-length genome subtype F sequences available in public databases. Recent dates of HIV-1 diagnoses and short genetic distances suggest a recent origin of this CRF. This is the tenth reported CRF_BF, the first apparently having originated outside of South America.