A validated in-house assay for HIV drug resistance mutation surveillance from dried blood spot specimens.

Despite increasing scale-up of antiretroviral therapy (ART) coverage, challenges related to adherence and HIV drug resistance (HIVDR) remain. The high cost of HIVDR surveillance is a persistent challenge with implementation in resource-constrained settings. Dried blood spot (DBS) specimens have been demonstrated to be a feasible alternative to plasma or serum for HIVDR genotyping and are more suitable for lower resource settings. There is a need for affordable HIVDR genotyping assays which can amplify HIV-1 sequences from DBS specimens, particularly those with low viral loads, at a low cost. Here, we present an in-house assay capable of reliably amplifying HIV-1 protease and partial reverse transcriptase genes from DBS specimens, which covers the complete World Health Organization 2009 list of drug resistance mutations under surveillance. DBS specimens were prepared using whole blood spiked with HIV-1 at concentrations of 10,000, 5000, 1000, and 500 copies/mL (n=30 for each concentration). Specimens were tested in triplicate. A two-step approach was used consisting of cDNA synthesis followed by nested PCR. The limit of detection of the assay was calculated to be approximately 5000 (95% CI: 3200-10,700) copies/mL for the protease gene and 3600 (95% CI: 2200-10,000) copies/mL for reverse transcriptase. The assay was observed to be most sensitive with higher viral load specimens (97.8% [95% CI: 92.2-99.7]) for both protease and reverse transcriptase at 10,000 copies/mL with performance decreasing with the use of specimens with lower viral loads (46.7% [36.1-57.5] and 60.0% [49.1-70.2] at 500 copies/mL for protease and reverse transcriptase, respectively). Ultimately, this assay presents a promising opportunity for use in resource-constrained settings. Future work should involve validation under field conditions including sub-optimal storage conditions and preparation of DBS with fingerprick blood in order to accurately reflect real-world collection scenarios.

HIV acquisition prior to entry into formal sex work: inference from next-generation viral sequencing.

OBJECTIVE: To infer the timing of HIV acquisition in relation to self-reported events in the sexual life course of adolescent girls and young women (AGYW) who self-identify as female sex workers (FSW) in Mombasa, Kenya. DESIGN: Next-generation viral sequencing of samples of AGYW living with HIV in the Transitions study, a cross-sectional bio-behavioural survey of AGYW aged 14-24 years in Mombasa, Kenya. METHOD: Dried blood spot specimens were collected from study participants ( n  = 37, all FSW). A portion of the HIV pol gene was sequenced using an in-house next-generation sequencing assay for HIV drug resistance mutation genotyping. Estimated time since infection (ETI) was inferred using the HIV EVO web-based tool ( https://hiv.biozentrum.unibas.ch/ETI/ ), and data on self-reported events were obtained from the survey. RESULTS: The median ETI among FSW was 3.4 (interquartile range = 1.7, 6.3) years, with a median ETI of 1.5 years prior to entry into formal sex work. We estimated that 74.1% (95% confidence interval = 53.7-88.9%) of participants living with HIV and who self-identified as FSW likely acquired HIV prior to self-identification as a sex worker. CONCLUSIONS: Findings suggest a large fraction of prevalent HIV infection among AGYW engaged in sex work stems from acquisition prior to entry into formal sex work. Current HIV prevention programs tailored for sex workers may miss key opportunities for HIV prevention as they are designed to reach women after entry into formal sex work, signaling a need for tailored programs to reach high-risk AGYW earlier on in their sexual life course.

Point-of-Care Tests for HIV Drug Resistance Monitoring: Advances and Potentials.

HIV/AIDS is a global public health crisis that is yet to be contained. Effective management of HIV drug resistance (HIVDR) supported by close resistance monitoring is essential in achieving the WHO 95-95-95 targets, aiming to end the AIDS epidemic by 2030. Point-of-care tests (POCT) enable decentralized HIVDR testing with a short turnaround time and minimal instrumental requirement, allowing timely initiation of effective antiretroviral therapy (ART) and regimen adjustment as needed. HIVDR POCT is of particular significance in an era when ART access is scaling up at a global level and enhanced HIVDR monitoring is urgently needed, especially for low-to-middle-income countries. This article provides an overview of the currently available technologies that have been applied or potentially used in HIVDR POCT. It may also benefit the continued research and development efforts toward more innovative HIVDR diagnostics.

Application of a Sanger-Based External Quality Assurance Strategy for the Transition of HIV-1 Drug Resistance Assays to Next Generation Sequencing.

The National Institute of Allergy and Infectious Diseases (NIAID) Virology Quality Assurance (VQA) established a robust proficiency testing program for Sanger sequencing (SS)-based HIV-1 drug resistance (HIVDR) testing in 2001. While many of the lessons learned during the development of such programs may also apply to next generation sequencing (NGS)-based HIVDR assays, challenges remain for the ongoing evaluation of NGS-based testing. These challenges include a proper assessment of assay accuracy and the reproducibility of low abundance variant detection, intra- and inter-assay performance comparisons among laboratories using lab-defined tests, and different data analysis pipelines designed for NGS. In collaboration with the World Health Organization (WHO) Global HIVDR Laboratory Network and the Public Health Agency of Canada, the Rush VQA program distributed archived proficiency testing panels to ten laboratories to evaluate internally developed NGS assays. Consensus FASTA files were submitted using 5%, 10%, and 20% variant detection thresholds, and scored based on the same criteria used for SS. This small study showed that the SS External Quality Assurance (EQA) approach can be used as a transitional strategy for using NGS to generate SS-like data and for ongoing performance while using NGS data from the same quality control materials to further evaluate NGS assay performance.

HIV-1 phylodynamic analysis among people who inject drugs in Pakistan correlates with trends in illicit opioid trade.

Pakistan is considered by the World Health Organization to currently have a "concentrated" HIV-1 epidemic due to a rapid rise in infections among people who inject drugs (PWID). Prevalence among the country's nearly 105,000 PWID is estimated to be 37.8% but has been shown to be higher in several large urban centers. A lack of public health resources, the common use of professional injectors and unsafe injection practices are believed to have fueled the outbreak. Here we evaluate the molecular characteristics of HIV-1 sequences (n = 290) from PWID in several Pakistani cities to examine transmission dynamics and the association between rates of HIV-1 transmission with regards to regional trends in opioid trafficking. Tip-to-tip (patristic) distance based phylogenetic cluster inferences and BEAST2 Bayesian phylodynamic analyses of time-stamped data were performed on HIV-1 pol sequences generated from dried blood spots collected from 1,453 PWID as part of a cross-sectional survey conducted in Pakistan during 2014/2015. Overall, subtype A1 strains were dominant (75.2%) followed by CRF02_AG (14.1%), recombinants/unassigned (7.2%), CRF35_AD (2.1%), G (1.0%) and C (0.3%). Nearly three quarters of the PWID HIV-1 sequences belonged to one of five distinct phylogenetic clusters. Just below half (44.4%) of individuals in the largest cluster (n = 118) did seek help injecting from professional injectors which was previously identified as a strong correlate of HIV-1 infection. Spikes in estimated HIV-1 effective population sizes coincided with increases in opium poppy cultivation in Afghanistan, Pakistan's western neighbor. Structured coalescent analysis was undertaken in order to investigate the spatial relationship of HIV-1 transmission among the various cities under study. In general terms, our analysis placed the city of Larkana at the center of the PWID HIV-1 epidemic in Pakistan which is consistent with previous epidemiological data.

Quality Control of Next-Generation Sequencing-Based HIV-1 Drug Resistance Data in Clinical Laboratory Information Systems Framework.

Next-generation sequencing (NGS) in HIV drug resistance (HIVDR) testing has the potential to improve both clinical and public health settings, however it challenges the normal operations of quality management systems to be more flexible due to its complexity, massive data generation, and rapidly evolving protocols. While guidelines for quality management in NGS data have previously been outlined, little guidance has been implemented for NGS-based HIVDR testing. This document summarizes quality control procedures for NGS-based HIVDR testing laboratories using a laboratory information systems (LIS) framework. Here, we focus in particular on the quality control measures applied on the final sequencing product aligned with the recommendations from the World Health Organization HIV Drug Resistance Laboratory Network.

Performance comparison of next generation sequencing analysis pipelines for HIV-1 drug resistance testing.

Next generation sequencing (NGS) is a trending new standard for genotypic HIV-1 drug resistance (HIVDR) testing. Many NGS HIVDR data analysis pipelines have been independently developed, each with variable outputs and data management protocols. Standardization of such analytical methods and comparison of available pipelines are lacking, yet may impact subsequent HIVDR interpretation and other downstream applications. Here we compared the performance of five NGS HIVDR pipelines using proficiency panel samples from NIAID Virology Quality Assurance (VQA) program. Ten VQA panel specimens were genotyped by each of six international laboratories using their own in-house NGS assays. Raw NGS data were then processed using each of the five different pipelines including HyDRA, MiCall, PASeq, Hivmmer and DEEPGEN. All pipelines detected amino acid variants (AAVs) at full range of frequencies (1~100%) and demonstrated good linearity as compared to the reference frequency values. While the sensitivity in detecting low abundance AAVs, with frequencies between 1~20%, is less a concern for all pipelines, their specificity dramatically decreased at AAV frequencies <2%, suggesting that 2% threshold may be a more reliable reporting threshold for ensured specificity in AAV calling and reporting. More variations were observed among the pipelines when low abundance AAVs are concerned, likely due to differences in their NGS read quality control strategies. Findings from this study highlight the need for standardized strategies for NGS HIVDR data analysis, especially for the detection of minority HIVDR variants.

A MiSeq-HyDRA platform for enhanced HIV drug resistance genotyping and surveillance.

Conventional HIV drug resistance (HIVDR) genotyping utilizes Sanger sequencing (SS) methods, which are limited by low data throughput and the inability of detecting low abundant drug resistant variants (LADRVs). Here we present a next generation sequencing (NGS)-based HIVDR typing platform that leverages the advantages of Illumina MiSeq and HyDRA Web. The platform consists of a fully validated sample processing protocol and HyDRA web, an open web portal that allows automated customizable NGS-based HIVDR data processing. This platform was characterized and validated using a panel of HIV-spiked plasma representing all major HIV-1 subtypes, pedigreed plasmids, HIVDR proficiency specimens and clinical specimens. All examined major HIV-1 subtypes were consistently amplified at viral loads of ≥1,000 copies/ml. The gross error rate of this platform was determined at 0.21%, and minor variations were reliably detected down to 0.50% in plasmid mixtures. All HIVDR mutations identifiable by SS were detected by the MiSeq-HyDRA protocol, while LADRVs at frequencies of 1~15% were detected by MiSeq-HyDRA only. As compared to SS approaches, the MiSeq-HyDRA platform has several notable advantages including reduced cost and labour, and increased sensitivity for LADRVs, making it suitable for routine HIVDR monitoring for both patient care and surveillance purposes.

Bioinformatic data processing pipelines in support of next-generation sequencing-based HIV drug resistance testing: the Winnipeg Consensus.

INTRODUCTION: Next-generation sequencing (NGS) has several advantages over conventional Sanger sequencing for HIV drug resistance (HIVDR) genotyping, including detection and quantitation of low-abundance variants bearing drug resistance mutations (DRMs). However, the high HIV genomic diversity, unprecedented large volume of data, complexity of analysis and potential for error pose significant challenges for data processing. Several NGS analysis pipelines have been developed and used in HIVDR research; however, the absence of uniformity in data processing strategies results in lack of consistency and comparability of outputs from different pipelines. To fill this gap, an international symposium on bioinformatic strategies for NGS-based HIVDR testing was held in February 2018 in Winnipeg, Canada, convening laboratory scientists, bioinformaticians and clinicians involved in four recently developed, publicly available NGS HIVDR pipelines. The goal of this symposium was to establish a consensus on effective bioinformatic strategies for NGS data management and its use for HIVDR reporting. DISCUSSION: Essential functionalities of an NGS HIVDR pipeline were divided into five analytic blocks: (1) NGS read quality control (QC)/quality assurance (QA); (2) NGS read alignment and reference mapping; (3) HIV variant calling and variant QC; (4) NGS HIVDR reporting; and (5) extended data applications and additional considerations for data management. The consensuses reached among the participants on all major aspects of these blocks are summarized here. They encompass not only recommended data management and analysis strategies, but also detailed bioinformatic approaches that help ensure accuracy of the derived HIVDR analysis outputs for both research and potential clinical use. CONCLUSIONS: While NGS is being adopted more broadly in HIVDR testing laboratories, data processing is often a bottleneck hindering its generalized application. The proposed standardization of NGS read QC/QA, read alignment and reference mapping, variant calling and QC, HIVDR reporting and relevant data management strategies in this "Winnipeg Consensus" may serve as a starting guideline for NGS HIVDR data processing that informs the refinement of existing pipelines and those yet to be developed. Moreover, the bioinformatic strategies presented here may apply more broadly to NGS data analysis of microbes harbouring significant genomic diversity.

Silencing the buzz: a new approach to population suppression of mosquitoes by feeding larvae double-stranded RNAs.

BACKGROUND: Mosquito-borne diseases threaten over half the world's human population, making the need for environmentally-safe mosquito population control tools critical. The sterile insect technique (SIT) is a biological control method that can reduce pest insect populations by releasing a large number of sterile males to compete with wild males for female mates to reduce the number of progeny produced. Typically, males are sterilized using radiation, but such methods can reduce their mating competitiveness. The method is also most effective if only males are produced, but this requires the development of effective sex-sorting methods. Recent efforts to use transgenic methods to produce sterile male mosquitoes have increased interest in using SIT to control some of our most serious disease vectors, but the release of genetically modified mosquitoes will undoubtedly encounter considerable delays as regulatory agencies deal with safety issues and public concerns. METHODS: Testis genes in the dengue vector Aedes aegypti were identified using a suppression subtractive hybridization technique. Mosquito larvae were fed double-stranded RNAs (dsRNAs) that targeted both the testis genes and a female sex determination gene (doublesex) to induce RNA interference (RNAi) -mediated sterility and inhibition of female development. Fertility and mating competiveness of the treated males were assessed in small-scale mating competition experiments. RESULTS: Feeding mosquito larvae dsRNAs targeting testis genes produced adult males with greatly reduced fertility; several dsRNAs produced males that were highly effective in competing for mates. RNAi-mediated knockdown of the female-specific isoform of doublesex was also effective in producing a highly male-biased population of mosquitoes, thereby overcoming the need to sex-sort insects before release. CONCLUSIONS: The sequence-specific gene-silencing mechanism of this RNAi technology renders it adaptable for species-specific application across numerous insect species. We envisage its use for traditional large-scale reared releases of mosquitoes and other pest insects, although the technology might also have potential for field-based control of mosquitoes where eggs deposited into a spiked larval site lead to the release of new sterile males.

Immunogenicity of sequences around HIV-1 protease cleavage sites: potential targets and population coverage analysis for a HIV vaccine targeting protease cleavage sites.

Developing an effective preventative vaccine against HIV-1 has proved to be a great challenge. The classical and proven vaccine approach has failed so far or produced a modest effect, new approaches are needed. In this study we evaluated the immunogenicity of the sequences around the protease cleavage sites (PCS) and the population coverage of a vaccine targeting HIV-1 PCS. The sequence conservation was evaluated by comparing entropy score of sequences around PCS with Gag and Pol. The immunogenicity of sequences around the 12 PCS (+10/-10 amino acids) was analyzed by identifying epitopes of HLA class I alleles in PCS region using four approaches: (1) identification of previously reported HLA class I allele epitopes around PCS region; (2) screening and validating epitopes of 8 HLA class I alleles common to most world populations using iTopia Epitope Discovery system and IFN-γ ELISpot assays; (3) screening of 151 patients of Pumwani cohort for PBMC IFN-γ ELISPOT responses to the subtype A and D consensus around PCS region; and (4) prediction of HLA alleles with epitopes around the PCS using NetMHCpan. Population coverage was calculated using the web-based analysis tool of the Immune Epitope Database based on HLA class I genotype frequencies from dbMHC database. The results showed that many HLA class I alleles have multiple epitopes in the 12 PCS regions, indicating sequence immunogenicity around PCS. Multiple epitopes of many HLA class I alleles common to >95% world populations have been identified around the 12 PCS region. Targeting these sites is a feasible vaccine approach.

HIV-1 clade D is associated with increased rates of CD4 decline in a Kenyan cohort.

HIV-1 is grouped phylogenetically into clades, which may impact rates of HIV-1 disease progression. Clade D infection in particular has been shown to be more pathogenic. Here we confirm in a Nairobi-based prospective female sex worker cohort (1985-2004) that Clade D (n = 54) is associated with a more rapid CD4 decline than clade A1 (n = 150, 20.6% vs 13.4% decline per year, 1.53-fold increase, p = 0.015). This was independent of "protective" HLA and country of origin (p = 0.053), which in turn were also independent predictors of the rate of CD4 decline (p = 0.026 and 0.005, respectively). These data confirm that clade D is more pathogenic than clade A1. The precise reason for this difference is currently unclear, and requires further study. This is first study to demonstrate difference in HIV-1 disease progression between clades while controlling for protective HLA alleles.

For protection from HIV-1 infection, more might not be better: a systematic analysis of HIV Gag epitopes of two alleles associated with different outcomes of HIV-1 infection.

A subset of women in the Pumwani Sex Worker Cohort, established in 1985 in Nairobi, Kenya, remains uninfected despite repeated high-risk exposure (HIV-exposed, seronegative [HESN]) through active sex work. This HESN phenotype is associated with several alleles of human leukocyte antigens (HLAs) and specific CD8(+) and CD4(+) T cell responses to HIV-1. The associations of HLA alleles with differential HIV-1 infection are most likely due to their different abilities to present antigen and the different immune responses they induce. The characteristics of epitopes of HLA alleles associated with different outcomes of HIV-1 infection might therefore point to a vital clue for developing an effective vaccine. In this study, we systematically analyzed HIV-1 clade A and D Gag CD8(+) T cell epitopes of two HLA class I alleles associated with different outcomes of HIV-1 infection. Binding affinity and off-rates of the identified epitopes were determined. Gamma interferon (IFN-γ) enzyme-linked immunospot (ELISpot) assays with patient peripheral blood mononuclear cells (PBMCs) validated the epitopes. Epitope-specific CD8(+) T cells were further phenotyped for memory markers with tetramer staining. Our study showed that the protective allele A*01:01 recognizes only three Gag epitopes. By contrast, B*07:02, the allele associated with susceptibility, binds 30 epitope variants. These two alleles differ most importantly in the spectrum of Gag epitopes they can present and not in affinity, off-rates, the location of the epitopes, or epitope-specific Tem/Tcm frequencies. The binding of more epitopes and strong IFN-gamma ELISpot responses are associated with susceptibility to HIV-1 infection, while more focused antigen recognition of multiple subtypes is protective. Rational vaccine design should take these observations into account.

A comparison of parallel pyrosequencing and sanger clone-based sequencing and its impact on the characterization of the genetic diversity of HIV-1.

BACKGROUND: Pyrosequencing technology has the potential to rapidly sequence HIV-1 viral quasispecies without requiring the traditional approach of cloning. In this study, we investigated the utility of ultra-deep pyrosequencing to characterize genetic diversity of the HIV-1 gag quasispecies and assessed the possible contribution of pyrosequencing technology in studying HIV-1 biology and evolution. METHODOLOGY/PRINCIPAL FINDINGS: HIV-1 gag gene was amplified from 96 patients using nested PCR. The PCR products were cloned and sequenced using capillary based Sanger fluorescent dideoxy termination sequencing. The same PCR products were also directly sequenced using the 454 pyrosequencing technology. The two sequencing methods were evaluated for their ability to characterize quasispecies variation, and to reveal sites under host immune pressure for their putative functional significance. A total of 14,034 variations were identified by 454 pyrosequencing versus 3,632 variations by Sanger clone-based (SCB) sequencing. 11,050 of these variations were detected only by pyrosequencing. These undetected variations were located in the HIV-1 Gag region which is known to contain putative cytotoxic T lymphocyte (CTL) and neutralizing antibody epitopes, and sites related to virus assembly and packaging. Analysis of the positively selected sites derived by the two sequencing methods identified several differences. All of them were located within the CTL epitope regions. CONCLUSIONS/SIGNIFICANCE: Ultra-deep pyrosequencing has proven to be a powerful tool for characterization of HIV-1 genetic diversity with enhanced sensitivity, efficiency, and accuracy. It also improved reliability of downstream evolutionary and functional analysis of HIV-1 quasispecies.

Identification and characterization of HLA-A*0301 epitopes in HIV-1 gag proteins using a novel approach.

Identification of CTL epitopes correlated to immune protection is important for the development of vaccines that enhance T-cell mediated immune responses. The correlation of positively selected amino acids (PS) of HIV-1 with host HLA alleles can identify regions containing potential T-cell epitopes. However, the specific epitopes have to be identified and characterized using overlapping peptides through T-cell functional assays. In this study we used a new approach to identify and characterize potential epitopes in the gag region containing PS mutations that significantly correlated with HLA-A*0301. The iTopia Epitope Discovery System was used to rapidly screen a panel of peptides overlapping the regions containing PS mutations and the peptides identified were assessed for relative affinity and complex stability. The potential epitopes were then validated by interferon gamma (IFN-gamma) ELISpot assays with patient PBMCs. Using this approach we identified/confirmed the predicted HLA-A*0301 epitopes in two regions of gag containing PS mutations V7I and K403R, one previously reported and the other novel. Five of the seven peptides that bound to A*0301 contained the K403R mutation and corresponded to the documented LARNCRAPRK-A3 supertype epitope. Two epitope variants, RASVLSGGK and RASILSGGK containing the V7I mutation, were identified using the iTopia Epitope Discovery System, however only the consensus variant (RAK9C) was confirmed using the ELISpot assay and it represents a novel A*0301 epitope.

Clade-specific evolution mediated by HLA-B*57/5801 in human immunodeficiency virus type 1 clade A1 p24.

HLA-B*57-mediated selection pressure leads to a typical escape pathway in human immunodeficiency virus type 1 (HIV-1) CD8 epitopes such as TW10. Whether this T242N pathway is shared by all clades remains unknown. We therefore assessed the nature of HLA-B*57 selection in a large, observational Kenyan cohort where clades A1 and D predominate. While T242N was ubiquitous in clade D HLA-B*57(+) subjects, this mutation was rare (15%) in clade A1. Instead, P243T and I247L were selected by clade A1-infected HLA-B*57 subjects but not by HLA-B*5801(+) subjects. Our data suggest that clade A1 consensus proline at Gag residue 243 might represent an inherent block to T242N escape in clade A1. We confirmed immunologically that P243T and I247L likely represent escape mutations. HLA-B*57 evolution also differed between clades in the KF11 and IW9 epitopes. A better understanding of clade-specific evolution is important for the development of HIV vaccines in regions with multiple clades.

An integrative bioinformatic approach for studying escape mutations in human immunodeficiency virus type 1 gag in the Pumwani Sex Worker Cohort.

Human immunodeficiency virus type 1 (HIV-1) is able to evade the host cytotoxic T-lymphocyte (CTL) response through a variety of escape avenues. Epitopes that are presented to CTLs are first processed in the presenting cell in several steps, including proteasomal cleavage, transport to the endoplasmic reticulum, binding by the HLA molecule, and finally presentation to the T-cell receptor. An understanding of the potential of the virus to escape CTL responses can aid in designing an effective vaccine. To investigate such a potential, we analyzed HIV-1 gag from 468 HIV-1-positive Kenyan women by using several bioinformatic approaches that allowed the identification of positively selected amino acids in the HIV-1 gag region and study of the effects that these mutations could have on the various stages of antigen processing. Correlations between positively selected residues and mean CD4 counts also allowed study of the effect of mutation on HIV disease progression. A number of mutations that could create or destroy proteasomal cleavage sites or reduce binding affinity of the transport antigen processing protein, effectively hindering epitope presentation, were identified. Many mutations correlated with the presence of specific HLA alleles and with lower or higher CD4 counts. For instance, the mutation V190I in subtype A1-infected individuals is associated with HLA-B*5802 (P = 4.73 x 10(-4)), a rapid-progression allele according to other studies, and also to a decreased mean CD4 count (P = 0.019). Thus, V190I is a possible HLA escape mutant. This method classifies many positively selected mutations across the entire gag region according to their potential for immune escape and their effect on disease progression.