Molecular analysis of gp41 sequences of HIV type 1 subtype C from India

Analysis of sequence diversity and selection pressure in HIV-1 clade C gp41 from India

Evaluation of viral diversity is critical for the rational design of treatment modalities against Human immunodeficiency virus (HIV). Predominated by HIV-1 clade C (HIV-1C), the epidemic in India represents the third largest population infected with HIV-1 globally. Glycoprotein 41 (gp41) is critical for viral replication and is a target for the design of therapeutic strategies. However, documentation of viral diversity of gp41 gene in infected individuals from India remains limited. Present study employed high throughput sequencing to examine variation in gp41 amplicons generated from blood derived viruses in 24 HIV-1C infected individuals from Mumbai, India. Sequence diversity profiles were documented in different functional domains of gp41. Furthermore, through a meta-analysis approach, all reported gp41 sequences from India (N = 70) were compared with those from South Africa (N = 126), country with the largest HIV epidemic globally, also predominated by HIV-1C. A total of 44 positions displayed statistically significant differential (p < 0.05) Shannon entropy in the two regions. This comparison also identified 11 codon sites undergoing distinct selection, 8 of which remained differentially selected in an extended comparison of data from Asia (N = 137) and Africa(N = 383). Assessment of correlated mutation networks associated with differentially selected residues revealed common as well as distinct interaction networks. Furthermore, codon usage analysis revealed 17 differentially selected codons (Mann-Whitney test, p < 0.001) in Asia and Africa. Dissimilar trends in GC content across codon positions were also observed. In depth understanding of these divergent evolutionary signatures through extended analysis with larger data-sets would assist development of effective interventions being considered for HIV-1C.

Effect of diversity in gp41 membrane proximal external region of primary HIV-1 Indian subtype C sequences on interaction with broadly neutralizing antibodies 4E10 and 10E8

Human Immunodeficiency Virus-1 Clade C (HIV-1C) dominates the AIDS epidemic in India, afflicting 2.1 million individuals within the country and more than 15 million people worldwide. Membrane proximal external region (MPER) is an attractive target for broadly neutralizing antibody (bNAb) based therapies. However, information on MPER sequence diversity from India is meagre due to limited sampling of primary viral sequences. In the present study, we examined the variation in MPER of HIV-1C from 24 individuals in Mumbai, India by high throughput sequencing of uncultured viral sequences. Deep sequencing of MPER (662-683; HXB2 envelope amino acid numbering) allowed quantification of intra-individual variation up to 65% at positions 662, 665, 668, 674 and 677 within this region. These variable positions included contact sites targeted by bNAbs 2F5, Z13e1, 4E10 as well as 10E8. Both major and minor epitope variants i.e. 'haplotypes' were generated for each sample dataset. A total of 23, 34 and 25 unique epitope haplotypes could be identified for bNAbs 2F5, Z13e1 and 4E10/10E8 respectively. Further analysis of 4E10 and 10E8 epitopes from our dataset and meta-analysis of previously reported HIV-1 sequences from India revealed 26 epitopes (7 India-specific), heretofore untested for neutralization sensitivity. Peptide-Ab docking predicted 13 of these to be non-binding to 10E8. ELISA, Surface Plasmon Resonance and peptide inhibition of HIV-1 neutralization assays were then performed which validated predicted weak/non-binding interactions for peptides corresponding to six of these epitopes. These results highlight the under-representation of 10E8 non-binding HIV-1C MPER sequences from India. Our study thus underscores the need for increased surveillance of primary circulating envelope sequences for development of efficacious bNAb-based interventions in India.

CodABC: a computational framework to coestimate recombination, substitution, and molecular adaptation rates by approximate Bayesian computation

The estimation of substitution and recombination rates can provide important insights into the molecular evolution of protein-coding sequences. Here, we present a new computational framework, called "CodABC," to jointly estimate recombination, substitution and synonymous and nonsynonymous rates from coding data. CodABC uses approximate Bayesian computation with and without regression adjustment and implements a variety of codon models, intracodon recombination, and longitudinal sampling. CodABC can provide accurate joint parameter estimates from recombining coding sequences, often outperforming maximum-likelihood methods based on more approximate models. In addition, CodABC allows for the inclusion of several nuisance parameters such as those representing codon frequencies, transition matrices, heterogeneity across sites or invariable sites. CodABC is freely available from http://code.google.com/p/codabc/, includes a GUI, extensive documentation and ready-to-use examples, and can run in parallel on multicore machines.

Broadly Neutralizing Antibody Responses in a Subset of HIV-1-Infected Individuals in Chennai, India

Identification of broadly neutralizing antibodies (NAbs) generated during the course of HIV-1 infection is essential for effective HIV-1 vaccine design. The magnitude and breadth of neutralizing activity in the sera from 46 antiretroviral treatment-naive HIV-1 clade C-infected individuals was measured in a single round infection assay using TZM-bl cells and multisubtype panel of env-pseudotyped viruses. Higher levels of NAb response (NAb titer 500 to >40 000) were measured in these patients against tier 1 and tier 2 viruses. The average magnitude of the NAb responses of chronically infected individuals against heterologous viruses was consistently higher than the response observed from individuals with long-term nonprogressor ( P = .086). To conclude, high titers of HIV-1 cross-neutralizing activity were observed in the sera from a subset of HIV-1-infected individuals in Chennai, India. Additional studies of the epitopes recognized by these antibodies may facilitate the discovery of an effective vaccine immunogen.

Variations in autologous neutralization and CD4 dependence of b12 resistant HIV-1 clade C env clones obtained at different time points from antiretroviral naïve Indian patients with recent infection

Background

Limited information is available on HIV-1 Indian clade C sensitivities to autologous antibodies during the course of natural infection. In the present study, a total of 37 complete envelope clones (Env) were amplified at different time points predominantly from the plasma of five Indian patients with recent HIV-1 infection and envelope-pseudotyped viruses were examined for their magnitude of sensitivity to autologous plasma antibodies during natural course of infection.

Results

Variable low levels of neutralization were consistently detected with contemporaneous autologous plasma. In contrast to clade B and African clade C HIV-1 envelopes, Env clones obtained from four patients were found to be resistant to IgG1b12. The majority of the Env clones were resistant to 2G12 and 2F5 due to the absence of the minimal motifs required for antibody recognition, but were sensitive to 4E10. Nonetheless, Env clones from one patient were found to be sensitive to 2G12, atypical for clade C, and one Env clone exhibited unusual sensitivity to 17b, suggesting spontaneous exposure of CD4i epitopes. Phylogenetic analysis revealed that Env clones were closely clustered within patients. Variation in the potential N-linked glycosylation pattern also appeared to be different in patients over the course of infection. Interestingly, we found that the sensitivity of Envs to contemporaneous autologous NAbs correlated positively with increased sensitivity to soluble CD4 and inversely with anti-CD4 antibody and Envs with increased NAb sensitivity were able to efficiently infect HeLa cells expressing low CD4.

Conclusion

Our data showed considerable variations in autologous neutralization of these early HIV-1 clade C Envs in each of these patients and indicate greater exposure to CD4 of Envs that showed increased autologous neutralization. Interestingly, Env clones obtained from a single patient at different time points were found to retain sensitivity to b12 antibody that binds to CD4 binding site in Env in contrast to Envs obtained from other patients. However, we did not find any association between increased b12 sensitivity of Envs obtained from this particular patient with their degree of exposure to CD4.

Highly complex neutralization determinants on a monophyletic lineage of newly transmitted subtype C HIV-1 Env clones from India

Little is known about the neutralization properties of HIV-1 in India to optimally design and test vaccines. For this reason, a functional Env clone was obtained from each of ten newly acquired, heterosexually transmitted HIV-1 infections in Pune, Maharashtra. These clones formed a phylogenetically distinct genetic lineage within subtype C. As Env-pseudotyped viruses the clones were mostly resistant to IgG1b12, 2G12 and 2F5 but all were sensitive to 4E10. When compared to a large multi-subtype panel of Env-pseudotyped viruses (subtypes B, C and CRF02_AG) in neutralization assays with a multi-subtype panel of HIV-1-positive plasma samples, the Indian Envs were remarkably complex. With the exception of the Indian Envs, results of a hierarchical clustering analysis showed a strong subtype association with the patterns of neutralization susceptibility. From these patterns we were able to identify 19 neutralization cluster-associated amino acid signatures in gp120 and 14 signatures in the ectodomain and cytoplasmic tail of gp41. We conclude that newly transmitted Indian Envs are antigenically complex in spite of close genetic similarity. Delineation of neutralization-associated amino acid signatures provides a deeper understanding of the antigenic structure of HIV-1 Env.

High replication fitness and transmission efficiency of HIV-1 subtype C from India: Implications for subtype C predominance

HIV-1 subtype C has been the predominant subtype throughout the course of the HIV-1 epidemic in India regardless of the geographic region of the country. In an effort to understand the mechanism of subtype C predominance in this country, we have investigated the in vitro replication fitness and transmission efficiency of HIV-1 subtypes A and C from India. Using a dual infection growth competition assay, we found that primary HIV-1 subtype C isolates had higher overall relative fitness in PBMC than subtype A primary isolates. Moreover, in an ex vivo cervical tissue derived organ culture, subtype C isolates displayed higher transmission efficiency across cervical mucosa than subtype A isolates. We found that higher fitness of subtype C was not due to a trans effect exerted by subtype C infected PBMC. A half genome A/C recombinant clone in which the 3' half of the viral genome of subtype A was replaced with the corresponding subtype C3' half, had similar replicative fitness as the parental subtype A. These results suggest that the higher replication fitness and transmission efficiency of subtype C virus compared to subtype A virus from India is most probably not due to the envelope gene alone and may be due to genes present within the 5' half of the viral genome or to a more complex interaction between the genes located within the two halves of the viral genome. These data provide a model to explain the asymmetric distribution of subtype C over other subtypes in India.

Genetic diversity of HIV type 1 subtype C env gene sequences from India

Considering the severity of the HIV-1 subtype C epidemic, data on the epidemiology and distribution of HIV subtypes in India are relatively sparse. Keeping this in view, 28 env gene sequences from patients were sequenced and analyzed. The samples were collected over a period of 10 years from 1995 to 2004. Assessment of the interisolate genetic distances of the study isolates, which were all subtype C, showed interisolate distances varying from 2 to 19% (mean: 14%) with the maximum diversity observed in the samples collected in 2003-2004. Analysis of the phylogenetic relationships among subtype C env sequences from six different countries and our study isolates revealed an overall star-like phylogeny with almost all sequences from India forming a monophyletic lineage. A lower diversity within the immunodominant epitopes was found. The data generated from this study should prove valuable for the production of vaccine against subtype C.

Human immunodeficiency virus type 1 drug resistance mutations in peripheral blood mononuclear cell proviral DNA among antiretroviral treatment-naive and treatment-experienced patients from Pune, India

The prevalence of HIV drug resistance (HIVDR) mutations in the HIV protease (PR) and reverse transcriptase (RT) genes was estimated from peripheral blood mononuclear cells (PBMCs) in a study population of 25 antiretroviral (ARV) therapy-naive and 50 ARV-experienced chronically infected patients from Pune city, Maharashtra State, western India. Of the 75 study HIV-1 sequences, 73 belonged to subtype C and 2 to subtype A1. On phylogenetic analysis, the study subtype C sequences sub clustered randomly with different Indian and non-Indian subtype C sequences, emphasizing the HIV-1 subtype C pol gene diversity. The heterosexual route was the most common route of transmission (74.67%). There were no observable HIVDR mutations in ARV-naive patients. The ARV-experienced patients had a history of exposure to nucleoside reverse transcriptase inhibitor and nonnucleoside reverse transcriptase inhibitor combinations. At least one HIVDR mutation in RT was observed in 29 (80.55%) of ARV-experienced patients with evidence of failing therapy. M184V was the most common observed HIVDR mutation. No PR major mutations were observed among ARV-experienced patients. A higher prevalence of proviral HIVDR mutations in PBMCs was associated with irregular adherence to therapy (p < 0.05) and HIV-1 RNA levels > 1000 copies/ml (p < 0.001).

Extensive cross-reactive neutralizing antibody response in Indian patients with limited genetic diversity of HIV-1

Genome sequence analysis of HIV-1 subtype C viruses from India shows monophyletic lineage and relatively limited genetic diversity. To understand its immunological implication, cross-reactivity of neutralizing antibody response was investigated. In primary screening, neutralizing antibody response to single heterologous primary HIV-1 subtype C isolate was assessed in plasma samples from 235 HIV-1 infected, anti-retroviral treatment naive individuals from Pune, India. Plasma samples that showed > or =90% neutralization and two randomly selected plasma samples that showed 50-60% neutralization were tested against a panel of primary HIV-1 subtype C isolates obtained from epidemiologically unlinked individuals from India. The neutralizing antibody response showed extensive cross-neutralization, suggesting presence of shared neutralization determinants among circulating HIV-1 subtype C viruses in India.