High Level of Pre-Treatment HIV-1 Drug Resistance and Its Association with HLA Class I-Mediated Restriction in the Pumwani Sex Worker Cohort.

BACKGROUND: We analyzed the prevalence of pre-antiretroviral therapy (ART) drug resistance mutations (DRMs) in a Kenyan population. We also examined whether host HLA class I genes influence the development of pre-ART DRMs. METHODS: The HIV-1 proviral DNAs were amplified from blood samples of 266 ART-naïve women from the Pumwani Sex Worker cohort of Nairobi, Kenya using a nested PCR method. The amplified HIV genomes were sequenced using next-generation sequencing technology. The prevalence of pre-ART DRMs was investigated. Correlation studies were performed between HLA class I alleles and HIV-1 DRMs. RESULTS: Ninety-eight percent of participants had at least one DRM, while 38% had at least one WHO surveillance DRM. M184I was the most prevalent clinically important variant, seen in 37% of participants. The DRMs conferring resistance to one or more integrase strand transfer inhibitors were also found in up to 10% of participants. Eighteen potentially relevant (p < 0.05) positive correlations were found between HLA class 1 alleles and HIV drug-resistant variants. CONCLUSIONS: High levels of HIV drug resistance were found in all classes of antiretroviral drugs included in the current first-line ART regimens in Africa. The development of DRMs may be influenced by host HLA class I-restricted immunity.

Deep sequencing of near full-length HIV-1 genomes from plasma identifies circulating subtype C and infrequent occurrence of AC recombinant form in Southern India.

India has the third largest number of HIV-1-infected individuals accounting for approximately 2.1 million people, with a predominance of circulating subtype C strains and a low prevalence of subtype A and A1C and BC recombinant forms, identified over the past two decades. Recovery of near full-length HIV-1 genomes from a plasma source coupled with advances in next generation sequencing (NGS) technologies and development of universal methods for amplifying whole genomes of HIV-1 circulating in a target geography or population provides the opportunity for a detailed analysis of HIV-1 strain identification, evolution and dynamics. Here we describe the development and implementation of approaches for HIV-1 NGS analysis in a southern Indian cohort. Plasma samples (n = 20) were obtained from HIV-1-confirmed individuals living in and around the city of Bengaluru. Near full-length genome recovery was obtained for 9 Indian HIV-1 patients, with recovery of full-length gag and env genes for 10 and 2 additional subjects, respectively. Phylogenetic analyses indicate the majority of sequences to be represented by subtype C viruses branching within a monophyletic clade, comprising viruses from India, Nepal, Myanmar and China and closely related to a southern African cluster, with a low prevalence of the A1C recombinant form also present. Development of algorithms for bespoke recovery and analysis at a local level will further aid clinical management of HIV-1 infected Indian subjects and delineate the progress of the HIV-1 pandemic in this and other geographical regions.

HIV-1 Subtypes and 5'LTR-Leader Sequence Variants Correlate with Seroconversion Status in Pumwani Sex Worker Cohort.

Within the Pumwani sex worker cohort, a subgroup remains seronegative, despite frequent exposure to HIV-1; some of them seroconverted several years later. This study attempts to identify viral variations in 5'LTR-leader sequences (5'LTR-LS) that might contribute to the late seroconversion. The 5'LTR-LS contains sites essential for replication and genome packaging, viz, primer binding site (PBS), major splice donor (SD), and major packaging signal (PS). The 5'LTR-LS of 20 late seroconverters (LSC) and 122 early seroconverters (EC) were amplified, cloned, and sequenced. HelixTree 6.4.3 was employed to classify HIV subtypes and sequence variants based on seroconversion status. We find that HIV-1 subtypes A1.UG and D.UG were overrepresented in the viruses infecting the LSC (P < 0.0001). Specific variants of PBS (Pc < 0.0001), SD1 (Pc < 0.0001), and PS (Pc < 0.0001) were present only in the viral population from EC or LSC. Combinations of PBS [PBS-2 (Pc < 0.0001) and PBS-3 (Pc < 0.0001)] variants with specific SD sequences were only seen in LSC or EC. Combinations of A1.KE or D with specific PBS and SD variants were only present in LSC or EC (Pc < 0.0001). Furthermore, PBS variants only present in LSC co-clustered with PBS references utilizing tRNAArg; whereas, the PBS variants identified only in EC co-clustered with PBS references using tRNALys3 and its variants. This is the first report that specific PBS, SD1, and PS sequence variants within 5'LTR-LS are associated with HIV-1 seroconversion, and it could aid designing effective anti-HIV strategies.

Influence of HLA class I haplotypes on HIV-1 seroconversion and disease progression in Pumwani sex worker cohort.

We examined the effect of HLA class I haplotypes on HIV-1 seroconversion and disease progression in the Pumwani sex worker cohort. This study included 595 HIV-1 positive patients and 176 HIV negative individuals. HLA-A, -B, and -C were typed to 4-digit resolution using sequence-based typing method. HLA class I haplotype frequencies were estimated using PyPop 32-0.6.0. The influence of haplotypes on time to seroconversion and CD4+ T cell decline to <200 cells/mm3 were analyzed by Kaplan-Meier analysis using SPSS 13.0. Before corrections for multiple comparisons, three 2-loci haplotypes were significantly associated with faster seroconversion, including A*23∶01-C*02∶02 (p = 0.014, log rank(LR) = 6.06, false-discovery rate (FDR) = 0.056), B*42∶01-C*17∶01 (p = 0.01, LR = 6.60, FDR = 0.08) and B*07∶02-C*07∶02 (p = 0.013, LR = 6.14, FDR = 0.069). Two A*74∶01 containing haplotypes, A*74∶01-B*15∶03 (p = 0.047, LR = 3.942, FDR = 0.068) and A*74∶01-B*15∶03-C*02∶02 (p = 0.045, LR = 4.01, FDR = 0.072) and B*14∶02-C*08∶02 (p = 0.021, LR = 5.36, FDR = 0.056) were associated with slower disease progression. Five haplotypes, including A*30∶02-B*45∶01 (p = 0.0008, LR = 11.183, FDR = 0.013), A*30∶02-C*16∶01 (p = 0.015, LR = 5.97, FDR = 0.048), B*53∶01-C*04∶01 (p = 0.010, LR = 6.61, FDR = 0.08), B*15∶10-C*03∶04 (p = 0.031, LR = 4.65, FDR = 0.062), and B*58∶01-C*03∶02 (p = 0.037, LR = 4.35, FDR = 0.066) were associated with faster progression to AIDS. After FDR corrections, only the associations of A*30∶02-B*45∶01 and A*30∶02-C*16∶01 with faster disease progression remained significant. Cox regression and deconstructed Kaplan-Meier survival analysis showed that the associations of haplotypes of A*23∶01-C*02∶02, B*07∶02-C*07∶02, A*74∶01-B*15∶03, A*74∶01-B*15∶03-C*02∶02, B*14∶02-C*08∶02 and B*58∶01-C*03∶02 with differential seroconversion or disease progression are due to the dominant effect of a single allele within the haplotypes. The true haplotype effect was observed with A*30∶02-B*45∶01, A*30∶02-C*16∶02, B*53∶01-C*04∶01 B*15∶10-C*03∶04, and B*42∶01-C*17∶01. In these cases, the presence of both alleles accelerated the disease progression or seroconversion than any of the single allele within the haplotypes. Our study showed that the true effects of HLA class I haplotypes on HIV seroconversion and disease progression exist and the associations of HLA class I haplotype can also be due to the dominant effect of a single allele within the haplotype.

Interplay between HIV-1 and Host Genetic Variation: A Snapshot into Its Impact on AIDS and Therapy Response.

As of February 2012, 50 circulating recombinant forms (CRFs) have been reported for HIV-1 while one CRF for HIV-2. Also according to HIV sequence compendium 2011, the HIV sequence database is replete with 414,398 sequences. The fact that there are CRFs, which are an amalgamation of sequences derived from six or more subtypes (CRF27_cpx (cpx refers to complex) is a mosaic with sequences from 6 different subtypes besides an unclassified fragment), serves as a testimony to the continual divergent evolution of the virus with its approximate 1% per year rate of evolution, and this phenomena per se poses tremendous challenge for vaccine development against HIV/AIDS, a devastating disease that has killed 1.8 million patients in 2010. Here, we explore the interaction between HIV-1 and host genetic variation in the context of HIV/AIDS and antiretroviral therapy response.