Developing high-throughput HIV incidence assay with pyrosequencing platform

Precision detection of recent HIV infections using high-throughput genomic incidence assay

HIV incidence is a key measure for tracking disease spread and identifying populations and geographic regions where new infections are most concentrated. The HIV sequence population provides a robust signal for the stage of infection. Large-scale and high-precision HIV sequencing is crucial for effective genomic incidence surveillance. We produced 1,034 full-length envelope gene sequences from a seroconversion cohort by conducting HIV microdrop sequencing and measuring the genomic incidence assay's genome similarity index (GSI) dynamics. The measured dynamics of 9 of 12 individuals aligned with the GSI distribution estimated independently using 417 publicly available incident samples. We enhanced the capacity to identify individuals with recent infections, achieving predicted detection accuracies of 92% (89%-94%) for cases within 6 months and 81% (74%-87%) for cases within 9 months. These accuracy levels agreed with the observed detection accuracy intervals of an independent validation data set. Additionally, we produced 131 full-length envelope gene sequences from eight individuals with chronic HIV infection. This analysis confirmed a false recency rate (FRR) of 0%, which was consistent with 162 publicly available chronic samples. The mean duration of recent infection (MDRI) was 238 (209-267) days, indicating an 83% improvement in performance compared to current recent infection testing algorithms. The shifted Poisson mixture model was then used to estimate the time since infection, and the model estimates showed an 88% consistency with the days post infection derived from HIV RNA test dates and/or seroconversion dates. HIV microdrop sequencing provides unique prospects for large-scale incidence surveillance using high-throughput sequencing. IMPORTANCE Accurate identification of recently infected individuals is vital for prioritizing specific populations for interventions, reducing onward transmission risks, and optimizing public health services. However, current HIV-specific antibody-based methods have not been satisfactory in accurately identifying incident cases, hindering the use of HIV recency testing for prevention efforts and partner protection. Genomic incidence assays offer a promising alternative for identifying recent infections. In our study, we used microdroplet technologies to produce a large number of complete HIV envelope gene sequences, enabling the accurate detection of early infection signs. We assessed the dynamics of the incidence assay's metrics and compared them with statistical models. Our approach demonstrated high accuracy in identifying individuals with recent infections, achieving predicted detection rates exceeding 90% within 6 months and over 80% within 9 months of infection. This high-resolution method holds significant potential for enhancing the effectiveness of HIV incidence screening for case-based surveillance in public health initiatives.

Full-spectrum HIV drug resistance mutation detection by high-resolution complete pol gene sequencing

BACKGROUND

Drug resistance mutation testing is a key element for HIV clinical management, informing effective treatment regimens. However, resistance screening in current clinical practice is limited in reporting linked cross-class resistance mutations and minority variants, both of which may increase the risk of virological failure.

METHODS

To address these limitations, we obtained 358 full-length pol gene sequences from 52 specimens of 20 HIV infected individuals by combining microdroplet amplification, unique molecular identifier (UMI) labeling, and long-read high-throughput sequencing.

RESULTS

We conducted a rigorous assessment of the accuracy of our pipeline for precision drug resistance mutation detection, verifying that a sequencing depth of 35 high-throughput reads achieved complete, error-free pol gene sequencing. We detected 26 distinct drug resistance mutations to Protease Inhibitors (PIs), Nucleoside Reverse Transcriptase Inhibitors (NRTIs), Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs), and Integrase Strand Transfer Inhibitors (INSTIs). We detected linked cross-class drug resistance mutations (PI+NRTI, PI+NNRTI, and NRTI+NNRTI) that confer cross-resistance to multiple drugs in different classes. Fourteen different types of minority mutations were also detected with frequencies ranging from 3.2% to 19%, and the presence of these mutations was verified by Sanger reference sequencing. We detected a putative transmitted drug resistance mutation (TDRM) in one individual that persisted for over seven months from the first sample collected at the acute stage of infection prior to seroconversion.

CONCLUSIONS

Our comprehensive drug resistance mutation profiling can advance clinical practice by reporting mutation linkage and minority variants to better guide antiretroviral therapy options.

Microdrop Human Immunodeficiency Virus Sequencing for Incidence and Drug Resistance Surveillance

BACKGROUND

Precise and cost-efficient human immunodeficiency virus (HIV) incidence and drug resistance surveillances are in high demand for the advancement of the 90-90-90 "treatment for all" target.

METHODS

We developed microdrop HIV sequencing for the HIV incidence and drug resistance assay (HIDA), a single-blood-draw surveillance tool for incidence and drug resistance mutation (DRM) detection. We amplified full-length HIV envelope and pol gene sequences within microdroplets, and this compartmental amplification with long-read high-throughput sequencing enabled us to recover multiple unique sequences.

RESULTS

We achieved greater precision in determining the stage of infection than current incidence assays, with a 1.2% false recency rate (proportion of misclassified chronic infections) and a 262-day mean duration of recent infection (average time span of recent infection classification) from 83 recently infected and 81 chronically infected individuals. Microdrop HIV sequencing demonstrated an increased capacity to detect minority variants and linked DRMs. By screening all 93 World Health Organization surveillance DRMs, we detected 6 pretreatment drug resistance mutations with 2.6%-13.2% prevalence and cross-linked mutations.

CONCLUSIONS

HIDA with microdrop HIV sequencing may promote global HIV real-time surveillance by serving as a precise and high-throughput cross-sectional survey tool that can be generalized for surveillance of other pathogens.

HIITE: HIV-1 incidence and infection time estimator

Motivation

Around 2.1 million new HIV-1 infections were reported in 2015, alerting that the HIV-1 epidemic remains a significant global health challenge. Precise incidence assessment strengthens epidemic monitoring efforts and guides strategy optimization for prevention programs. Estimating the onset time of HIV-1 infection can facilitate optimal clinical management and identify key populations largely responsible for epidemic spread and thereby infer HIV-1 transmission chains. Our goal is to develop a genomic assay estimating the incidence and infection time in a single cross-sectional survey setting.

Results

We created a web-based platform, HIV-1 incidence and infection time estimator (HIITE), which processes envelope gene sequences using hierarchical clustering algorithms and informs the stage of infection, along with time since infection for incident cases. HIITE's performance was evaluated using 585 incident and 305 chronic specimens' envelope gene sequences collected from global cohorts including HIV-1 vaccine trial participants. HIITE precisely identified chronically infected individuals as being chronic with an error less than 1% and correctly classified 94% of recently infected individuals as being incident. Using a mixed-effect model, an incident specimen's time since infection was estimated from its single lineage diversity, showing 14% prediction error for time since infection. HIITE is the first algorithm to inform two key metrics from a single time point sequence sample. HIITE has the capacity for assessing not only population-level epidemic spread but also individual-level transmission events from a single survey, advancing HIV prevention and intervention programs.

Availability and implementation

Web-based HIITE and source code of HIITE are available at http://www.hayounlee.org/software.html.

Supplementary information

Supplementary data are available at Bioinformatics online.

Detection of IgG3 antibodies specific to the human immunodeficiency virus type 1 (HIV-1) p24 protein as marker for recently acquired infection

Reducing the risk of human immunodeficiency virus type 1 (HIV-1) transmission is still a public health priority. The development of effective control strategies relies on the quantification of the effects of prophylactic and therapeutic measures in disease incidence. Although several assays can be used to estimate HIV incidence, these estimates are limited by the poor performance of these assays in distinguishing recent from long-standing infections. To address such limitation, we have developed an assay to titrate p24-specific IgG3 antibodies as a marker of recent infection. The assay is based on a recombinant p24 protein capable to detect total IgG antibodies in sera using a liquid micro array and enzyme-linked immunosorbent assay. Subsequently, the assay was optimised to detect and titrate anti-p24 IgG3 responses in a panel of sequential specimens from seroconverters over 24 months. The kinetics of p24-specific IgG3 titres revealed a transient peak in the 4 to 5-month period after seroconversion. It was followed by a sharp decline, allowing infections with less than 6 months to be distinguished from older ones. The developed assay exhibited a mean duration of recent infection of 144 days and a false-recent rate of ca. 14%. Our findings show that HIV-1 p24-specific IgG3 titres can be used as a tool to evaluate HIV incidence in serosurveys and to monitor the efficacy of vaccines and other transmission control strategies.

The HIV Genomic Incidence Assay Meets False Recency Rate and Mean Duration of Recency Infection Performance Standards

HIV incidence is a primary metric for epidemic surveillance and prevention efficacy assessment. HIV incidence assay performance is evaluated via false recency rate (FRR) and mean duration of recent infection (MDRI). We conducted a meta-analysis of 438 incident and 305 chronic specimens' HIV envelope genes from a diverse global cohort. The genome similarity index (GSI) accurately characterized infection stage across diverse host and viral factors. All except one chronic specimen had GSIs below 0.67, yielding a FRR of 0.33 [0-0.98] %. We modeled the incidence assay biomarker dynamics with a logistic link function assuming individual variabilities in a Beta distribution. The GSI probability density function peaked close to 1 in early infection and 0 around two years post infection, yielding MDRI of 420 [361, 467] days. We tested the assay by newly sequencing 744 envelope genes from 59 specimens of 21 subjects who followed from HIV negative status. Both standardized residuals and Anderson-Darling tests showed that the test dataset was statistically consistent with the model biomarker dynamics. This is the first reported incidence assay meeting the optimal FRR and MDRI performance standards. Signatures of HIV gene diversification can allow precise cross-sectional surveillance with a desirable temporal range of incidence detection.

Pairwise diversity and tMRCA as potential markers for HIV infection recency

Intrahost human immunodeficiency virus (HIV)-1 diversity increases linearly over time. We assessed the extent to which mean pairwise distances and the time to the most recent common ancestor (tMRCA) inferred from intrahost HIV-1C env sequences were associated with the estimated time of HIV infection. Data from a primary HIV-1C infection study in Botswana were used for this analysis (N = 42). A total of 2540 HIV-1C env gp120 variable loop region 1 to conserved region 5 (V1C5) of the HIV-1 envelope gp120 viral sequences were generated by single genome amplification and sequencing, with an average of 61 viral sequences per participant and 11 sequences per time point per participant. Raw pairwise distances were calculated for each time point and participant using the ape package in R software. The tMRCA was estimated using phylogenetic inference implemented in Bayesian Evolutionary Analysis by Sampling Trees v1.8.2. Pairwise distances and tMRCA were significantly associated with the estimated time since HIV infection (both P < 0.001). Taking into account multiplicity of HIV infection strengthened these associations. HIV-1C env-based pairwise distances and tMRCA can be used as potential markers for HIV recency. However, the tMRCA estimates demonstrated no advantage over the pairwise distances estimates.

Moving towards a reliable HIV incidence test - current status, resources available, future directions and challenges ahead

In 2011 the Incidence Assay Critical Path Working Group reviewed the current state of HIV incidence assays and helped to determine a critical path to the introduction of an HIV incidence assay. At that time the Consortium for Evaluation and Performance of HIV Incidence Assays (CEPHIA) was formed to spur progress and raise standards among assay developers, scientists and laboratories involved in HIV incidence measurement and to structure and conduct a direct independent comparative evaluation of the performance of 10 existing HIV incidence assays, to be considered singly and in combinations as recent infection test algorithms. In this paper we report on a new framework for HIV incidence assay evaluation that has emerged from this effort over the past 5 years, which includes a preliminary target product profile for an incidence assay, a consensus around key performance metrics along with analytical tools and deployment of a standardized approach for incidence assay evaluation. The specimen panels for this evaluation have been collected in large volumes, characterized using a novel approach for infection dating rules and assembled into panels designed to assess the impact of important sources of measurement error with incidence assays such as viral subtype, elite host control of viraemia and antiretroviral treatment. We present the specific rationale for several of these innovations, and discuss important resources for assay developers and researchers that have recently become available. Finally, we summarize the key remaining steps on the path to development and implementation of reliable assays for monitoring HIV incidence at a population level.

Enhancement of viral escape in HIV-1 Nef by STEP vaccination

OBJECTIVE

Properly priming cytotoxic T-lymphocyte (CTL) responses is an important task in HIV-1 vaccination. However, the STEP trial showed no efficacy even though the vaccine elicited HIV-specific CTL responses. Our study is to investigate whether or not the STEP vaccine enhanced viral escape in infected volunteers.

METHODS

The signature of viral escape, the presence of multiple escape variants, could be falsely represented by the existence of multiple founder viruses. Therefore, we use a mathematical model to designate STEP study patients with infections from a single founder virus. We then conduct permutation tests on each of 9988 Gag, Pol, and Nef overlapping peptides to identify epitopes with significant differences in diversity between the vaccine and placebo groups using previously published STEP trial sequence data.

RESULTS

We identify signatures of vaccine-enhanced viral escape within HIV-1 Nef from the STEP trial. Vaccine-treated patients showed a greater level of epitope diversity in one of the immunodomiant epitopes, EVGFPVRPQVPL (Nef65-76), compared with placebo-treated patients (P = 0.0038). In the other three Nef epitopes, there is a marginally significant difference in the epitope diversity between the vaccine and placebo group (P < 0.1). This greater epitope diversity was neither due to any difference in infection duration nor overall nef gene diversity between the two groups, suggesting that the increase in viral escape was likely mediated by vaccine-induced T-cell responses.

CONCLUSION

Viral escape in Nef is elevated preferentially in STEP vaccine-treated individuals, suggesting that vaccination primarily modulated initial CTL responses. Our observations provide important insights into improving vaccine-primed first immune control.

Analysis of Viral Diversity in Relation to the Recency of HIV-1C Infection in Botswana

BACKGROUND

Cross-sectional, biomarker methods to determine HIV infection recency present a promising and cost-effective alternative to the repeated testing of uninfected individuals. We evaluate a viral-based assay that uses a measure of pairwise distances (PwD) to identify HIV infection recency, and compare its performance with two serologic incidence assays, BED and LAg. In addition, we assess whether combination BED plus PwD or LAg plus PwD screening can improve predictive accuracy by reducing the likelihood of a false-recent result.

METHODS

The data comes from 854 time-points and 42 participants enrolled in a primary HIV-1C infection study in Botswana. Time points after treatment initiation or with evidence of multiplicity of infection were excluded from the final analysis. PwD was calculated from quasispecies generated using single genome amplification and sequencing. We evaluated the ability of PwD to correctly classify HIV infection recency within <130, <180 and <360 days post-seroconversion using Receiver Operator Characteristics (ROC) methods. Following a secondary PwD screening, we quantified the reduction in the relative false-recency rate (rFRR) of the BED and LAg assays while maintaining a sensitivity of either 75, 80, 85 or 90%.

RESULTS

The final analytic sample consisted of 758 time-points from 40 participants. The PwD assay was more accurate in classifying infection recency for the 130 and 180-day cut-offs when compared with the recommended LAg and BED thresholds. A higher AUC statistic confirmed the superior predictive performance of the PwD assay for the three cut-offs. When used for combination screening, the PwD assay reduced the rFRR of the LAg assay by 52% and the BED assay by 57.8% while maintaining a 90% sensitivity for the 130 and 180-day cut-offs respectively.

CONCLUSION

PwD can accurately determine HIV infection recency. A secondary PwD screening reduces misclassification and increases the accuracy of serologic-based assays.

Molecular clock of HIV-1 envelope genes under early immune selection

Background

The molecular clock hypothesis that genes or proteins evolve at a constant rate is a key tool to reveal phylogenetic relationships among species. Using the molecular clock, we can trace an infection back to transmission using HIV-1 sequences from a single time point. Whether or not a strict molecular clock applies to HIV-1’s early evolution in the presence of immune selection has not yet been fully examined.

Results

We identified molecular clock signatures from 1587 previously published HIV-1 full envelope gene sequences obtained since acute infection in 15 subjects. Each subject’s sequence diversity linearly increased during the first 150 days post infection, with rates ranging from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.54 \times 10^{ - 5}$$\end{document}1.54×10-5 to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3.91 \times 10^{ - 5}$$\end{document}3.91×10-5 with a mean of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2.69 \times 10^{ - 5}$$\end{document}2.69×10-5 per base per day. The rate of diversification for 12 out of the 15 subjects was comparable to the neutral evolution rate. While temporal diversification was consistent with evolution patterns in the absence of selection, mutations from the founder virus were highly clustered on statistically identified selection sites, which diversified more than 65 times faster than non-selection sites. By mathematically quantifying deviations from the molecular clock under various selection scenarios, we demonstrate that the deviation from a constant clock becomes negligible as multiple escape lineages emerge. The most recent common ancestor of a virus pair from distinct escape lineages is most likely the transmitted founder virus, indicating that HIV-1 molecular dating is feasible even after the founder viruses are no longer detectable.

Conclusions

The ability of HIV-1 to escape from immune surveillance in many different directions is the driving force of molecular clock persistence. This finding advances our understanding of the robustness of HIV-1’s molecular clock under immune selection, implying the potential for molecular dating.

Electronic supplementary material

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

Fecal bacterial microbiome diversity in chronic HIV-infected patients in China

The purpose of this study was to identify fecal bacterial microbiome changes in patients with chronic human immunodeficiency virus (HIV) infection in China. Bacterial 16S rRNA genes were amplified, sequenced (454 pyrosequencing), and clustered into operational taxonomic units using the QIIME software. Relative abundance at the phylum and genus levels were calculated. Alpha diversity was determined by Chao 1 and observed-species indices, and beta diversity was determined by double principal component analysis using the estimated phylogeny-based unweighted Unifrac distance matrices. Fecal samples of the patients with chronic HIV-infection tended to be enriched with bacteria of the phyla Firmicutes (47.20%±0.43 relative abundance) and Proteobacteria (37.21%±0.36) compared with those of the non-HIV infected controls (17.95%±0.06 and 3.81%±0.02, respectively). Members of the genus Bilophila were exclusively detected in samples of the non-HIV infected controls. Bacteroides and arabacteroides were more abundant in the chronic HIV-infected patients. Our study indicated that chronic HIV-infected patients in China have a fecal bacterial microbiome composition that is largely different from that found in non-HIV infected controls, and further study is needed to evaluate whether microbiome changes play a role in disease complications in the distal gut, including opportunistic infections.

SPMM: estimating infection duration of multivariant HIV-1 infections

MOTIVATION

Illustrating how HIV-1 is transmitted and how it evolves in the following weeks is an important step for developing effective vaccination and prevention strategies. It is currently possible through DNA sequencing to account for the diverse array of viral strains within an infected individual. This provides an unprecedented opportunity to pinpoint when each patient was infected and which viruses were transmitted.

RESULTS

Here we develop a mathematical tool for early HIV-1 evolution within a subject whose infection originates either from a single or multiple viral variants. The shifted Poisson mixture model (SPMM) provides a quantitative guideline for segregating viral lineages, which in turn enables us to assess when a subject was infected. The infection duration estimated by SPMM showed a statistically significant linear relationship with that by Fiebig laboratory staging (P = 0.00059) among 37 acutely infected subjects. Our tool provides a functional approach to understanding early genetic diversity, one of the most important parameters for deciphering HIV-1 transmission and predicting the rate of disease progression.

AVAILABILITY AND IMPLEMENTATION

SPMM, webserver, is available at http://www.hayounlee.org/web-tools.html.

CONTACT

hayoun@usc.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Quasispecies Analyses of the HIV-1 Near-full-length Genome With Illumina MiSeq

Human immunodeficiency virus type-1 (HIV-1) exhibits high between-host genetic diversity and within-host heterogeneity, recognized as quasispecies. Because HIV-1 quasispecies fluctuate in terms of multiple factors, such as antiretroviral exposure and host immunity, analyzing the HIV-1 genome is critical for selecting effective antiretroviral therapy and understanding within-host viral coevolution mechanisms. Here, to obtain HIV-1 genome sequence information that includes minority variants, we sought to develop a method for evaluating quasispecies throughout the HIV-1 near-full-length genome using the Illumina MiSeq benchtop deep sequencer. To ensure the reliability of minority mutation detection, we applied an analysis method of sequence read mapping onto a consensus sequence derived from de novo assembly followed by iterative mapping and subsequent unique error correction. Deep sequencing analyses of aHIV-1 clone showed that the analysis method reduced erroneous base prevalence below 1% in each sequence position and discarded only < 1% of all collected nucleotides, maximizing the usage of the collected genome sequences. Further, we designed primer sets to amplify the HIV-1 near-full-length genome from clinical plasma samples. Deep sequencing of 92 samples in combination with the primer sets and our analysis method provided sufficient coverage to identify >1%-frequency sequences throughout the genome. When we evaluated sequences of pol genes from 18 treatment-naïve patients' samples, the deep sequencing results were in agreement with Sanger sequencing and identified numerous additional minority mutations. The results suggest that our deep sequencing method would be suitable for identifying within-host viral population dynamics throughout the genome.

Identifying Recent HIV Infections: From Serological Assays to Genomics

In this paper, we review serological and molecular based methods to identify HIV infection recency. The accurate identification of recent HIV infection continues to be an important research area and has implications for HIV prevention and treatment interventions. Longitudinal cohorts that follow HIV negative individuals over time are the current gold standard approach, but they are logistically challenging, time consuming and an expensive enterprise. Methods that utilize cross-sectional testing and biomarker information have become an affordable alternative to the longitudinal approach. These methods use well-characterized biological makers to differentiate between recent and established HIV infections. However, recent results have identified a number of limitations in serological based assays that are sensitive to the variability in immune responses modulated by HIV subtypes, viral load and antiretroviral therapy. Molecular methods that explore the dynamics between the timing of infection and viral evolution are now emerging as a promising approach. The combination of serological and molecular methods may provide a good solution to identify recent HIV infection in cross-sectional data. As part of this review, we present the advantages and limitations of serological and molecular based methods and their potential complementary role for the identification of HIV infection recency.

Assessing Biases in the Evaluation of Classification Assays for HIV Infection Recency

Identifying recent HIV infection cases has important public health and clinical implications. It is essential for estimating incidence rates to monitor epidemic trends and evaluate the effectiveness of interventions. Detecting recent cases is also important for HIV prevention given the crucial role that recently infected individuals play in disease transmission, and because early treatment onset can improve the clinical outlook of patients while reducing transmission risk. Critical to this enterprise is the development and proper assessment of accurate classification assays that, based on cross-sectional samples of viral sequences, help determine infection recency status. In this work we assess some of the biases present in the evaluation of HIV recency classification algorithms that rely on measures of within-host viral diversity. Particularly, we examine how the time since infection (TSI) distribution of the infected subjects from which viral samples are drawn affect performance metrics (e.g., area under the ROC curve, sensitivity, specificity, accuracy and precision), potentially leading to misguided conclusions about the efficacy of classification assays. By comparing the performance of a given HIV recency assay using six different TSI distributions (four simulated TSI distributions representing different epidemic scenarios, and two empirical TSI distributions), we show that conclusions about the overall efficacy of the assay depend critically on properties of the TSI distribution. Moreover, we demonstrate that an assay with high overall classification accuracy, mainly due to properly sorting members of the well-represented groups in the validation dataset, can still perform notoriously poorly when sorting members of the less represented groups. This is an inherent issue of classification and diagnostics procedures that is often underappreciated. Thus, this work underscores the importance of acknowledging and properly addressing evaluation biases when proposing new HIV recency assays.

Differentiation of acute from chronic hepatitis C virus infection by nonstructural 5B deep sequencing: a population-level tool for incidence estimation

The ability to classify acute versus chronic hepatitis C virus (HCV) infections at the time of diagnosis is desirable to improve the quality of surveillance information. The aim of this study was to differentiate acute from chronic HCV infections utilizing deep sequencing. HCV nonstructural 5B (NS5B) amplicons (n = 94) were generated from 77 individuals (13 acute and 64 chronic HCV infections) in British Columbia, Canada, with documented seroconversion time frames. Amplicons were deep sequenced and HCV genomic diversity was measured by Shannon entropy (SE) and a single nucleotide variant (SNV) analysis. The relationship between each diversity measure and the estimated days since infection was assessed using linear mixed models, and the ability of each diversity measure to differentiate acute from chronic infections was assessed using generalized estimating equations. Both SE and the SNV diversity measures were significantly different for acute versus chronic infections (P < 0.009). NS5B nucleotide diversity continued to increase for at least 3 years postinfection. Among individuals with the least uncertainty with regard to duration of infection (n = 39), the area under the receiver operating characteristic curve (AUROC) was high (0.96 for SE; 0.98 for SNV). Although the AUROCs were lower (0.86 for SE; 0.80 for SNV) when data for all individuals were included, they remain sufficiently high for epidemiological purposes. Synonymous mutations were the primary discriminatory variable accounting for over 78% of the measured genetic diversity.

CONCLUSIONS

NS5B sequence diversity assessed by deep sequencing can differentiate acute from chronic HCV infections and, with further validation, could become a powerful population-level surveillance tool for incidence estimation.

A comparison of two measures of HIV diversity in multi-assay algorithms for HIV incidence estimation

BACKGROUND

Multi-assay algorithms (MAAs) can be used to estimate HIV incidence in cross-sectional surveys. We compared the performance of two MAAs that use HIV diversity as one of four biomarkers for analysis of HIV incidence.

METHODS

Both MAAs included two serologic assays (LAg-Avidity assay and BioRad-Avidity assay), HIV viral load, and an HIV diversity assay. HIV diversity was quantified using either a high resolution melting (HRM) diversity assay that does not require HIV sequencing (HRM score for a 239 base pair env region) or sequence ambiguity (the percentage of ambiguous bases in a 1,302 base pair pol region). Samples were classified as MAA positive (likely from individuals with recent HIV infection) if they met the criteria for all of the assays in the MAA. The following performance characteristics were assessed: (1) the proportion of samples classified as MAA positive as a function of duration of infection, (2) the mean window period, (3) the shadow (the time period before sample collection that is being assessed by the MAA), and (4) the accuracy of cross-sectional incidence estimates for three cohort studies.

RESULTS

The proportion of samples classified as MAA positive as a function of duration of infection was nearly identical for the two MAAs. The mean window period was 141 days for the HRM-based MAA and 131 days for the sequence ambiguity-based MAA. The shadows for both MAAs were <1 year. Both MAAs provided cross-sectional HIV incidence estimates that were very similar to longitudinal incidence estimates based on HIV seroconversion.

CONCLUSIONS

MAAs that include the LAg-Avidity assay, the BioRad-Avidity assay, HIV viral load, and HIV diversity can provide accurate HIV incidence estimates. Sequence ambiguity measures obtained using a commercially-available HIV genotyping system can be used as an alternative to HRM scores in MAAs for cross-sectional HIV incidence estimation.

Deep sequencing: becoming a critical tool in clinical virology

Population (Sanger) sequencing has been the standard method in basic and clinical DNA sequencing for almost 40 years; however, next-generation (deep) sequencing methodologies are now revolutionizing the field of genomics, and clinical virology is no exception. Deep sequencing is highly efficient, producing an enormous amount of information at low cost in a relatively short period of time. High-throughput sequencing techniques have enabled significant contributions to multiples areas in virology, including virus discovery and metagenomics (viromes), molecular epidemiology, pathogenesis, and studies of how viruses to escape the host immune system and antiviral pressures. In addition, new and more affordable deep sequencing-based assays are now being implemented in clinical laboratories. Here, we review the use of the current deep sequencing platforms in virology, focusing on three of the most studied viruses: human immunodeficiency virus (HIV), hepatitis C virus (HCV), and influenza virus.

Assessment of ambiguous base calls in HIV-1 pol population sequences as a biomarker for identification of recent infections in HIV-1 incidence studies

An increase in the proportion of ambiguous base calls in HIV-1 pol population sequences during the course of infection has been demonstrated in different study populations, and sequence ambiguity thresholds to classify infections as recent or nonrecent have been suggested. The aim of our study was to evaluate sequence ambiguities as a candidate biomarker for use in an HIV-1 incidence assay using samples from antiretroviral treatment-naive seroconverters with known durations of infection (German HIV-1 Seroconverter Study). We used 2,203 HIV-1 pol population sequences derived from 1,334 seroconverters to assess the sequence ambiguity method (SAM). We then compared the serological incidence BED capture enzyme immunoassay (BED-CEIA) with the SAM for a subset of 723 samples from 495 seroconverters and evaluated a multianalyte algorithm that includes BED-CEIA results, SAM results, viral loads, and CD4 cell counts for 453 samples from 325 seroconverters. We observed a significant increase in the proportion of sequence ambiguities with the duration of infection. A sequence ambiguity threshold of 0.5% best identified recent infections with 76.7% accuracy. The mean duration of recency was determined to be 208 (95% confidence interval, 196 to 221) days. In the subset analysis, BED-CEIA achieved a significantly higher accuracy than the SAM (84.6 versus 75.5%, P < 0.001) and results were concordant for 64.2% (464/723) of the samples. Also, the multianalyte algorithm did not show better accuracy than the BED-CEIA (83.4 versus 84.3%, P = 0.786). In conclusion, the SAM and the multianalyte algorithm including SAM were inferior to the BED-CEIA, and the proportion of sequence ambiguities is therefore not a preferable biomarker for HIV-1 incidence testing.