Rapid Sequencing of Complete env Genes from Primary HIV-1 Samples

Optimized SMRT-UMI protocol produces highly accurate sequence datasets from diverse populations-Application to HIV-1 quasispecies

Pathogen diversity resulting in quasispecies can enable persistence and adaptation to host defenses and therapies. However, accurate quasispecies characterization can be impeded by errors introduced during sample handling and sequencing, which can require extensive optimizations to overcome. We present complete laboratory and bioinformatics workflows to overcome many of these hurdles. The Pacific Biosciences single molecule real-time platform was used to sequence polymerase-chain reaction (PCR) amplicons derived from cDNA templates tagged with unique molecular identifiers (SMRT-UMI). Optimized laboratory protocols were developed through extensive testing of different sample preparation conditions to minimize between-template recombination during PCR. The use of UMI allowed accurate template quantitation as well as removal of point mutations introduced during PCR and sequencing to produce a highly accurate consensus sequence from each template. Production of highly accurate sequences from the large datasets produced from SMRT-UMI sequencing is facilitated by a novel bioinformatic pipeline, Probabilistic Offspring Resolver for Primer IDs (PORPIDpipeline). PORPIDpipeline automatically filters and parses circular consensus reads by sample, identifies and discards reads with UMIs likely created from PCR and sequencing errors, generates consensus sequences, checks for contamination within the dataset, and removes any sequence with evidence of PCR recombination, heteroduplex formation, or early cycle PCR errors. The optimized SMRT-UMI sequencing and PORPIDpipeline methods presented here represent a highly adaptable and established starting point for accurate sequencing of diverse pathogens. These methods are illustrated through characterization of human immunodeficiency virus quasispecies in a virus transmitter-recipient pair of individuals.

HIV-PULSE: a long-read sequencing assay for high-throughput near full-length HIV-1 proviral genome characterization

A deep understanding of the composition of the HIV-1 reservoir is necessary for the development of targeted therapies and the evaluation of curative efforts. However, current near full-length (NFL) HIV-1 proviral genome sequencing assays are based on labor-intensive and costly principles of repeated PCRs at limiting dilution, restricting their scalability. To address this, we developed a high-throughput, long-read sequencing assay called HIV-PULSE (HIV Proviral UMI-mediated Long-read Sequencing). This assay uses unique molecular identifiers (UMIs) to tag individual HIV-1 genomes, allowing for the omission of the limiting dilution step and enabling long-range PCR amplification of many NFL genomes in a single PCR reaction, while simultaneously overcoming poor single-read accuracy. We optimized the assay using HIV-infected cell lines and then applied it to blood samples from 18 individuals living with HIV on antiretroviral therapy, yielding a total of 1308 distinct HIV-1 genomes. Benchmarking against the widely applied Full-Length Individual Proviral Sequencing assay revealed similar sensitivity (11 vs 18%) and overall good concordance, although at a significantly higher throughput. In conclusion, HIV-PULSE is a cost-efficient and scalable assay that allows for the characterization of the HIV-1 proviral landscape, making it an attractive method to study the HIV-1 reservoir composition and dynamics.

Evaluation of Archival HIV DNA in Brain and Lymphoid Tissues

HIV reservoirs persist in anatomic compartments despite antiretroviral therapy (ART). Characterizing archival HIV DNA in the central nervous system (CNS) and other tissues is crucial to inform cure strategies. We evaluated paired autopsy brain-frontal cortex (FC), occipital cortex (OCC), and basal ganglia (BG)-and peripheral lymphoid tissues from 63 people with HIV. Participants passed away while virally suppressed on ART at the last visit and without evidence of CNS opportunistic disease. We quantified total HIV DNA in all participants and obtained full-length HIV-envelope (FL HIV-env) sequences from a subset of 14 participants. We detected HIV DNA (gag) in most brain (65.1%) and all lymphoid tissues. Lymphoid tissues had higher HIV DNA levels than the brain (P < 0.01). Levels of HIV gag between BG and FC were similar (P > 0.2), while OCC had the lowest levels (P = 0.01). Females had higher HIV DNA levels in tissues than males (gag, P = 0.03; 2-LTR, P = 0.05), suggesting possible sex-associated mechanisms for HIV reservoir persistence. Most FL HIV-env sequences (n = 143) were intact, while 42 were defective. Clonal sequences were found in 8 out of 14 participants, and 1 participant had clonal defective sequences in the brain and spleen, suggestive of cell migration. From 10 donors with paired brain and lymphoid sequences, we observed evidence of compartmentalized sequences in 2 donors. Our data further the idea that the brain is a site for archival HIV DNA during ART where compartmentalized provirus may occur in a subset of people. Future studies assessing FL HIV-provirus and replication competence are needed to further evaluate the HIV reservoirs in tissues. IMPORTANCE HIV infection of the brain is associated with adverse neuropsychiatric outcomes, despite efficient antiretroviral treatment. HIV may persist in reservoirs in the brain and other tissues, which can seed virus replication if treatment is interrupted, representing a major challenge to cure HIV. We evaluated reservoirs and genetic features in postmortem brain and lymphoid tissues from people with HIV who passed away during suppressed HIV replication. We found a differential distribution of HIV reservoirs across brain regions which was lower than that in lymphoid tissues. We observed that most HIV reservoirs in tissues had intact envelope sequences, suggesting they could potentially generate replicative viruses. We found that women had higher HIV reservoir levels in brain and lymphoid tissues than men, suggesting possible sex-based mechanisms of maintenance of HIV reservoirs in tissues, warranting further investigation. Characterizing the archival HIV DNA in tissues is important to inform future HIV cure strategies.

Antigen pressure from two founder viruses induces multiple insertions at a single antibody position to generate broadly neutralizing HIV antibodies

Vaccination strategies aimed at maturing broadly neutralizing antibodies (bnAbs) from naïve precursors are hindered by unusual features that characterize these Abs, including insertions and deletions (indels). Longitudinal studies of natural HIV infection cases shed light on the complex processes underlying bnAb development and have suggested a role for superinfection as a potential enhancer of neutralization breadth. Here we describe the development of a potent bnAb lineage that was elicited by two founder viruses to inform vaccine design. The V3-glycan targeting bnAb lineage (PC39-1) was isolated from subtype C-infected IAVI Protocol C elite neutralizer, donor PC39, and is defined by the presence of multiple independent insertions in CDRH1 that range from 1-11 amino acids in length. Memory B cell members of this lineage are predominantly atypical in phenotype yet also span the class-switched and antibody-secreting cell compartments. Development of neutralization breadth occurred concomitantly with extensive recombination between founder viruses before each virus separated into two distinct population "arms" that evolved independently to escape the PC39-1 lineage. Ab crystal structures show an extended CDRH1 that can help stabilize the CDRH3. Overall, these findings suggest that early exposure of the humoral system to multiple related Env molecules could promote the induction of bnAbs by focusing Ab responses to conserved epitopes.

Optimized SMRT-UMI protocol produces highly accurate sequence datasets from diverse populations - application to HIV-1 quasispecies

Pathogen diversity resulting in quasispecies can enable persistence and adaptation to host defenses and therapies. However, accurate quasispecies characterization can be impeded by errors introduced during sample handling and sequencing which can require extensive optimizations to overcome. We present complete laboratory and bioinformatics workflows to overcome many of these hurdles. The Pacific Biosciences single molecule real-time platform was used to sequence PCR amplicons derived from cDNA templates tagged with universal molecular identifiers (SMRT-UMI). Optimized laboratory protocols were developed through extensive testing of different sample preparation conditions to minimize between-template recombination during PCR and the use of UMI allowed accurate template quantitation as well as removal of point mutations introduced during PCR and sequencing to produce a highly accurate consensus sequence from each template. Handling of the large datasets produced from SMRT-UMI sequencing was facilitated by a novel bioinformatic pipeline, Probabilistic Offspring Resolver for Primer IDs (PORPIDpipeline), that automatically filters and parses reads by sample, identifies and discards reads with UMIs likely created from PCR and sequencing errors, generates consensus sequences, checks for contamination within the dataset, and removes any sequence with evidence of PCR recombination or early cycle PCR errors, resulting in highly accurate sequence datasets. The optimized SMRT-UMI sequencing method presented here represents a highly adaptable and established starting point for accurate sequencing of diverse pathogens. These methods are illustrated through characterization of human immunodeficiency virus (HIV) quasispecies.

HIV-PULSE: A long-read sequencing assay for high-throughput near full-length HIV-1 proviral genome characterization

A deep understanding of the composition of the HIV-1 reservoir is necessary for the development of targeted therapies and the evaluation of curative efforts. However, current near full-length (NFL) HIV-1 proviral genome sequencing assays are based on labor-intensive and costly principles of repeated PCRs at limiting dilution, restricting their scalability. To address this, we developed a high-throughput, long-read sequencing assay called HIV-PULSE (HIV P roviral U MI-mediated L ong-read Se quencing). This assay uses unique molecular identifiers (UMIs) to tag individual HIV-1 genomes, allowing for the omission of the limiting dilution step and enabling long-range PCR amplification of many NFL genomes in a single PCR reaction, while simultaneously overcoming poor single-read accuracy. We optimized the assay using HIV-infected cell lines and then applied it to blood samples from 18 individuals living with HIV on antiretroviral therapy, yielding a total of 1,308 distinct HIV-1 genomes. Benchmarking against the widely applied Full-Length Individual Proviral Sequencing assay revealed similar sensitivity (11% vs 18%) and overall good concordance, though at a significantly higher throughput. In conclusion, HIV-PULSE is a cost-efficient and scalable assay that allows for the characterization of the HIV-1 proviral landscape, making it an attractive method to study the HIV-1 reservoir composition and dynamics.

Distinct HIV-1 Population Structure across Meningeal and Peripheral T Cells and Macrophage Lineage Cells

HIV-1 sequence population structure among brain and nonbrain cellular compartments is incompletely understood. Here, we compared proviral pol and env high-quality consensus single-molecule real-time (SMRT) sequences derived from CD3⁺ T cells and CD14⁺ macrophage lineage cells from meningeal or peripheral (spleen, blood) tissues obtained at autopsy from two individuals with viral suppression on antiretroviral therapy (ART). Phylogenetic analyses showed strong evidence of population structure between CD3⁺ and CD14⁺ virus populations. Distinct env variable-region characteristics were also found between CD3⁺ and CD14⁺ viruses. Furthermore, shared macrophage-tropic amino acid residues (env) and drug resistance mutations (pol) between meningeal and peripheral virus populations were consistent with the meninges playing a role in viral gene flow across the blood-brain barrier. Overall, our results point toward potential functional differences among meningeal and peripheral CD3⁺ and CD14⁺ virus populations and a complex evolutionary history driven by distinct selection pressures and/or viral gene flow. IMPORTANCE Different cell types and/or tissues may serve as a reservoir for HIV-1 during ART-induced viral suppression. We compared proviral pol and env sequences from CD3⁺ T cells and CD14⁺ macrophage lineage cells from brain and nonbrain tissues from two virally suppressed individuals. We found strong evidence of viral population structure among cells/tissues, which may result from distinct selective pressures across cell types and anatomic sites.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.

Ultradeep HIV-1 Proviral Envelope Sequencing Reveals Complex Population Structure within and between Brain and Splenic Tissues

Understanding tissue-based HIV-1 proviral population structure is important for improving treatment strategies for individuals with HIV-associated neurological disorders (HAND). Previous analyses have revealed HIV-1 envelope (env) population structure between brain and peripheral tissues as well as Env functional differences, especially in individuals with HAND. Furthermore, population structure has been detected among different anatomical locations in the brain itself, although such patterns are inconsistent across individuals and less strongly associated with the presence/absence of HAND. Here, we utilized the Pacific Biosciences single-molecule real-time (SMRT) high-throughput technology to generate thousands of sequences for each tissue, along with phylogenetic and distance-based analyses, to investigate env sequences from paired brain and spleen samples from eight individuals with/without HAND. To account for the high error rate associated with SMRT sequencing, we used a clustering approach to identify high-quality consensus sequences representative of ≥10 reads ("HQCS10"). In parallel, we characterized variable regions from nonclustered sequences to identify potential functional differences. We found evidence for significant population structure between brain and spleen tissues, as well as among brain tissues and within the same brain tissue, in individuals both with and without HAND. Variable region analysis showed differences in length and charge among brain and nonbrain tissues as well as within the brain, suggesting possible functional differences. Our results demonstrate the complexity of HIV-1 env structure/gene flow among tissues and support the concept that selective pressures in different tissue microenvironments drive viral evolution and adaptation. IMPORTANCE Understanding the evolution of HIV-1 in the brain compared to other tissues is important for improving treatment strategies for individuals with HIV-associated neurological disorders (HAND). We utilized high-throughput sequencing technology to generate thousands of full-length env sequences from paired brain and spleen samples from eight individuals with/without HAND. We found significant viral population structure for participants both with and without HAND, providing robust evidence for the brain as a compartmentalized tissue and potentially a viral reservoir. We also found striking genetic differences between virus populations, even from the same tissue, suggesting the potential for functional differences and the possibility for multiple evolutionary pathways that result in similar tropisms and/or other tissue-adapted characteristics. Our results demonstrate the complexity of viral population structure within the brain and suggest that analysis of peripheral blood samples alone may not be fully informative with respect to improving strategies to treat or eradicate HIV-1.

Fc-mediated effector function contributes to the in vivo antiviral effect of an HIV neutralizing antibody

Treatment of HIV infection with either antiretroviral (ARV) therapy or neutralizing monoclonal antibodies (NAbs) leads to a reduction in HIV plasma virus. Both ARVs and NAbs prevent new rounds of viral infection, but NAbs may have the additional capacity to accelerate the loss of virus-infected cells through Fc gamma receptor (FcγR)-mediated effector functions, which should affect the kinetics of plasma-virus decline. Here, we formally test the role of effector function in vivo by comparing the rate and timing of plasma-virus clearance in response to a single-dose treatment with either unmodified NAb or those with either reduced or augmented Fc function. When infused into viremic simian HIV (SHIV)-infected rhesus macaques, there was a 21% difference in slope of plasma-virus decline between NAb and NAb with reduced Fc function. NAb engineered to increase FcγRIII binding and improve antibody-dependent cellular cytotoxicity (ADCC) in vitro resulted in arming of effector cells in vivo, yet led to viral-decay kinetics similar to NAbs with reduced Fc function. These studies show that the predominant mechanism of antiviral activity of HIV NAbs is through inhibition of viral entry, but that Fc function can contribute to the overall antiviral activity, making them distinct from standard ARVs.

VRC34-Antibody Lineage Development Reveals How a Required Rare Mutation Shapes the Maturation of a Broad HIV-Neutralizing Lineage

Rare mutations have been proposed to restrict the development of broadly neutralizing antibodies against HIV-1, but this has not been explicitly demonstrated. We hypothesized that such rare mutations might be identified by comparing broadly neutralizing and non-broadly neutralizing branches of an antibody-developmental tree. Because sequences of antibodies isolated from the fusion peptide (FP)-targeting VRC34-antibody lineage suggested it might be suitable for such rare mutation analysis, we carried out next-generation sequencing (NGS) on B cell transcripts from donor N123, the source of the VRC34 lineage, and functionally and structurally characterized inferred intermediates along broadly neutralizing and poorly neutralizing developmental branches. The broadly neutralizing VRC34.01 branch required the rare heavy-chain mutation Y33P to bind FP, whereas the early bifurcated VRC34.05 branch did not require this rare mutation and evolved less breadth. Our results demonstrate how a required rare mutation can restrict development and shape the maturation of a broad HIV-1-neutralizing antibody lineage.

RIFRAF: a frame-resolving consensus algorithm

Motivation

Protein coding genes can be studied using long-read next generation sequencing. However, high rates of indel sequencing errors are problematic, corrupting the reading frame. Even the consensus of multiple independent sequence reads retains indel errors. To solve this problem, we introduce Reference-Informed Frame-Resolving multiple-Alignment Free template inference algorithm (RIFRAF), a sequence consensus algorithm that takes a set of error-prone reads and a reference sequence and infers an accurate in-frame consensus. RIFRAF uses a novel structure, analogous to a two-layer hidden Markov model: the consensus is optimized to maximize alignment scores with both the set of noisy reads and with a reference. The template-to-reads component of the model encodes the preponderance of indels, and is sensitive to the per-base quality scores, giving greater weight to more accurate bases. The reference-to-template component of the model penalizes frame-destroying indels. A local search algorithm proceeds in stages to find the best consensus sequence for both objectives.

Results

Using Pacific Biosciences SMRT sequences from an HIV-1 env clone, NL4-3, we compare our approach to other consensus and frame correction methods. RIFRAF consistently finds a consensus sequence that is more accurate and in-frame, especially with small numbers of reads. It was able to perfectly reconstruct over 80% of consensus sequences from as few as three reads, whereas the best alternative required twice as many. RIFRAF is able to achieve these results and keep the consensus in-frame even with a distantly related reference sequence. Moreover, unlike other frame correction methods, RIFRAF can detect and keep true indels while removing erroneous ones.

Availability and implementation

RIFRAF is implemented in Julia, and source code is publicly available at https://github.com/MurrellGroup/Rifraf.jl.

Supplementary information

Supplementary data are available at Bioinformatics online.

Long-read amplicon denoising

Long-read next-generation amplicon sequencing shows promise for studying complete genes or genomes from complex and diverse populations. Current long-read sequencing technologies have challenging error profiles, hindering data processing and incorporation into downstream analyses. Here we consider the problem of how to reconstruct, free of sequencing error, the true sequence variants and their associated frequencies from PacBio reads. Called 'amplicon denoising', this problem has been extensively studied for short-read sequencing technologies, but current solutions do not always successfully generalize to long reads with high indel error rates. We introduce two methods: one that runs nearly instantly and is very accurate for medium length reads and high template coverage, and another, slower method that is more robust when reads are very long or coverage is lower. On two Mock Virus Community datasets with ground truth, each sequenced on a different PacBio instrument, and on a number of simulated datasets, we compare our two approaches to each other and to existing algorithms. We outperform all tested methods in accuracy, with competitive run times even for our slower method, successfully discriminating templates that differ by a just single nucleotide. Julia implementations of Fast Amplicon Denoising (FAD) and Robust Amplicon Denoising (RAD), and a webserver interface, are freely available.

Rapid and Focused Maturation of a VRC01-Class HIV Broadly Neutralizing Antibody Lineage Involves Both Binding and Accommodation of the N276-Glycan

The VH1-2 restricted VRC01-class of antibodies targeting the HIV envelope CD4 binding site are a major focus of HIV vaccine strategies. However, a detailed analysis of VRC01-class antibody development has been limited by the rare nature of these responses during natural infection and the lack of longitudinal sampling of such responses. To inform vaccine strategies, we mapped the development of a VRC01-class antibody lineage (PCIN63) in the subtype C infected IAVI Protocol C neutralizer PC063. PCIN63 monoclonal antibodies had the hallmark VRC01-class features and demonstrated neutralization breadth similar to the prototype VRC01 antibody, but were 2- to 3-fold less mutated. Maturation occurred rapidly within ∼24 months of emergence of the lineage and somatic hypermutations accumulated at key contact residues. This longitudinal study of broadly neutralizing VRC01-class antibody lineage reveals early binding to the N276-glycan during affinity maturation, which may have implications for vaccine design.

The Neutralizing Antibody Response to the HIV-1 Env Protein

BACKGROUND

A vaccine able to elicit broadly neutralizing antibodies capable of blocking infection by global viruses has not been achieved, and remains a key public health challenge.

OBJECTIVE

During infection, a robust strain-specific neutralizing response develops in most people, but only a subset of infected people develop broadly neutralizing antibodies. Understanding how and why these broadly neutralizing antibodies develop has been a focus of the HIV-1 vaccine field for many years, and has generated extraordinary insights into the neutralizing response to HIV-1 infection.

RESULTS

This review describes the features, targets and developmental pathways of early strainspecific antibodies and later broadly neutralizing antibodies, and explores the reasons such broad antibodies are not more commonly elicited during infection.

CONCLUSION

The insights from these studies have been harnessed for the development of pioneering new vaccine approaches that seek to drive B cell maturation towards breadth. Overall, this review describes how findings from infected donors have impacted on active and passive immunization approaches that seek to prevent HIV-1 infection.

Next-generation sequencing of HIV-1 single genome amplicons

The analysis of HIV-1 sequences has helped understand the viral molecular epidemiology, monitor the development of antiretroviral drug resistance, and design candidate vaccines. The introduction of single genome amplification (SGA) has been a major advancement in the field, allowing for the characterization of multiple sequences per patient while preserving linkage among polymorphisms in the same viral genome copy. Sequencing of SGA amplicons is performed by capillary Sanger sequencing, which presents low throughput, requires a high amount of template, and is highly sensitive to template/primer mismatching. In order to meet the increasing demand for HIV-1 SGA amplicon sequencing, we have developed a platform based on benchtop next-generation sequencing (NGS) (IonTorrent) accompanied by a bioinformatics pipeline capable of running on computer resources commonly available at research laboratories. During assay validation, the NGS-based sequencing of 10 HIV-1 env SGA amplicons was fully concordant with Sanger sequencing. The field test was conducted on plasma samples from 10 US Navy and Marine service members with recent HIV-1 infection (sampling interval: 2005–2010; plasma viral load: 5,884–194,984 copies/ml). The NGS analysis of 101 SGA amplicons (median: 10 amplicons/individual) showed within-individual viral sequence profiles expected in individuals at this disease stage, including individuals with highly homogeneous quasispecies, individuals with two highly homogeneous viral lineages, and individuals with heterogeneous viral populations. In a scalability assessment using the Ion Chef automated system, 41/43 tested env SGA amplicons (95%) multiplexed on a single Ion 318 chip showed consistent gene-wide coverage >50×. With lower sample requirements and higher throughput, this approach is suitable to support the increasing demand for high-quality and cost-effective HIV-1 sequences in fields such as molecular epidemiology, and development of preventive and therapeutic strategies.

Full-Length Envelope Analyzer (FLEA): A tool for longitudinal analysis of viral amplicons

Next generation sequencing of viral populations has advanced our understanding of viral population dynamics, the development of drug resistance, and escape from host immune responses. Many applications require complete gene sequences, which can be impossible to reconstruct from short reads. HIV env, the protein of interest for HIV vaccine studies, is exceptionally challenging for long-read sequencing and analysis due to its length, high substitution rate, and extensive indel variation. While long-read sequencing is attractive in this setting, the analysis of such data is not well handled by existing methods. To address this, we introduce FLEA (Full-Length Envelope Analyzer), which performs end-to-end analysis and visualization of long-read sequencing data. FLEA consists of both a pipeline (optionally run on a high-performance cluster), and a client-side web application that provides interactive results. The pipeline transforms FASTQ reads into high-quality consensus sequences (HQCSs) and uses them to build a codon-aware multiple sequence alignment. The resulting alignment is then used to infer phylogenies, selection pressure, and evolutionary dynamics. The web application provides publication-quality plots and interactive visualizations, including an annotated viral alignment browser, time series plots of evolutionary dynamics, visualizations of gene-wide selective pressures (such as dN/dS) across time and across protein structure, and a phylogenetic tree browser. We demonstrate how FLEA may be used to process Pacific Biosciences HIV env data and describe recent examples of its use. Simulations show how FLEA dramatically reduces the error rate of this sequencing platform, providing an accurate portrait of complex and variable HIV env populations. A public instance of FLEA is hosted at http://flea.datamonkey.org. The Python source code for the FLEA pipeline can be found at https://github.com/veg/flea-pipeline. The client-side application is available at https://github.com/veg/flea-web-app. A live demo of the P018 results can be found at http://flea.murrell.group/view/P018.

Viral populations constantly evolve and diversify. In this article we introduce a method, FLEA, for reconstructing and visualizing the details of evolutionary changes. FLEA specifically processes data from sequencing platforms that generate reads that are long, but error-prone. To study the evolutionary dynamics of entire genes during viral infection, data is collected via long-read sequencing at discrete time points, allowing us to understand how the virus changes over time. However, the experimental and sequencing process is imperfect, so the resulting data contain not only real evolutionary changes, but also mutations and other genetic artifacts caused by sequencing errors. Our method corrects most of these errors by combining thousands of erroneous sequences into a much smaller number of unique consensus sequences that represent biologically meaningful variation. The resulting high-quality sequences are used for further analysis, such as building an evolutionary tree that tracks and interprets the genetic changes in the viral population over time. FLEA is open source, and is freely available online.

Vaccine-Induced Protection from Homologous Tier 2 SHIV Challenge in Nonhuman Primates Depends on Serum-Neutralizing Antibody Titers

Passive administration of HIV neutralizing antibodies (nAbs) can protect macaques from hard-to-neutralize (tier 2) chimeric simian-human immunodeficiency virus (SHIV) challenge. However, conditions for nAb-mediated protection after vaccination have not been established. Here, we selected groups of 6 rhesus macaques with either high or low serum nAb titers from a total of 78 animals immunized with recombinant native-like (SOSIP) Env trimers. Repeat intrarectal challenge with homologous tier 2 SHIV_BG505 led to rapid infection in unimmunized and low-titer animals. High-titer animals, however, demonstrated protection that was gradually lost as nAb titers waned over time. An autologous serum ID₅₀ nAb titer of ∼1:500 afforded more than 90% protection from medium-dose SHIV infection. In contrast, antibody-dependent cellular cytotoxicity and T cell activity did not correlate with protection. Therefore, Env protein-based vaccination strategies can protect against hard-to-neutralize SHIV challenge in rhesus macaques by inducing tier 2 nAbs, provided appropriate neutralizing titers can be reached and maintained.

Maintenance and reappearance of extremely divergent intra-host HIV-1 variants

Understanding genetic variation in human immunodeficiency virus (HIV) is clinically and immunologically important for patient treatment and vaccine development. We investigated the longitudinal intra-host genetic variation of HIV in over 3,000 individuals in the US National HIV Surveillance System with at least four reported HIV-1 polymerase (pol) sequences. In this population, we identified 149 putative instances of superinfection (i.e. an individual sequentially infected with genetically divergent, polyphyletic viruses). Unexpectedly, we discovered a group of 240 individuals with consecutively sampled viral strains that were >0.015 substitutions/site divergent, despite remaining monophyletic in the phylogeny. Viruses in some of these individuals had a maximum genetic divergence approaching that found between two random, unrelated HIV-1 subtype-B pol sequences within the US population. Individuals with these highly divergent viruses tended to be diagnosed nearly a decade earlier in the epidemic than people with superinfection or virus with less intra-host genetic variation, and they had distinct transmission risk factor profiles. To better understand this genetic variation in cases with extremely divergent, monophyletic viruses, we performed molecular clock phylogenetic analysis. Our findings suggest that, like Hepatitis C virus, extremely divergent HIV lineages can be maintained within an individual and reemerge over a period of years.

Development of broadly neutralizing antibodies in HIV-1 infected elite neutralizers

Broadly neutralizing antibodies (bNAbs), able to prevent viral entry by diverse global viruses, are a major focus of HIV vaccine design, with data from animal studies confirming their ability to prevent HIV infection. However, traditional vaccine approaches have failed to elicit these types of antibodies. During chronic HIV infection, a subset of individuals develops bNAbs, some of which are extremely broad and potent. This review describes the immunological and virological factors leading to the development of bNAbs in such "elite neutralizers". The features, targets and developmental pathways of bNAbs from their precursors have been defined through extraordinarily detailed within-donor studies. These have enabled the identification of epitope-specific commonalities in bNAb precursors, their intermediates and Env escape patterns, providing a template for vaccine discovery. The unusual features of bNAbs, such as high levels of somatic hypermutation, and precursors with unusually short or long antigen-binding loops, present significant challenges in vaccine design. However, the use of new technologies has led to the isolation of more than 200 bNAbs, including some with genetic profiles more representative of the normal immunoglobulin repertoire, suggesting alternate and shorter pathways to breadth. The insights from these studies have been harnessed for the development of optimized immunogens, novel vaccine regimens and improved delivery schedules, which are providing encouraging data that an HIV vaccine may soon be a realistic possibility.

Ultradeep single-molecule real-time sequencing of HIV envelope reveals complete compartmentalization of highly macrophage-tropic R5 proviral variants in brain and CXCR4-using variants in immune and peripheral tissues

Despite combined antiretroviral therapy (cART), HIV+ patients still develop neurological disorders, which may be due to persistent HIV infection and selective evolution in brain tissues. Single-molecule real-time (SMRT) sequencing technology offers an improved opportunity to study the relationship among HIV isolates in the brain and lymphoid tissues because it is capable of generating thousands of long sequence reads in a single run. Here, we used SMRT sequencing to generate ~ 50,000 high-quality full-length HIV envelope sequences (> 2200 bp) from seven autopsy tissues from an HIV+/cART+ subject, including three brain and four non-brain sites. Sanger sequencing was used for comparison with SMRT data and to clone functional pseudoviruses for in vitro tropism assays. Phylogenetic analysis demonstrated that brain-derived HIV was compartmentalized from HIV outside the brain and that the variants from each of the three brain tissues grouped independently. Variants from all peripheral tissues were intermixed on the tree but independent of the brain clades. Due to the large number of sequences, a clustering analysis at three similarity thresholds (99, 99.5, and 99.9%) was also performed. All brain sequences clustered exclusive of any non-brain sequences at all thresholds; however, frontal lobe sequences clustered independently of occipital and parietal lobes. Translated sequences revealed potentially functional differences between brain and non-brain sequences in the location of putative N-linked glycosylation sites (N-sites), V1 length, V3 charge, and the number of V4 N-sites. All brain sequences were predicted to use the CCR5 co-receptor, while most non-brain sequences were predicted to use CXCR4 co-receptor. Tropism results were confirmed by in vitro infection assays. The study is the first to use a SMRT sequencing approach to study HIV compartmentalization in tissues and supports other reports of limited trafficking between brain and non-brain sequences during cART. Due to the long sequence length, we could observe changes along the entire envelope gene, likely caused by differential selective pressure in the brain that may contribute to neurological disease.

Coupling of Single Molecule, Long Read Sequencing with IMGT/HighV-QUEST Analysis Expedites Identification of SIV gp140-Specific Antibodies from scFv Phage Display Libraries

The simian immunodeficiency virus (SIV)/macaque model of human immunodeficiency virus (HIV)/acquired immunodeficiency syndrome pathogenesis is critical for furthering our understanding of the role of antibody responses in the prevention of HIV infection, and will only increase in importance as macaque immunoglobulin (IG) gene databases are expanded. We have previously reported the construction of a phage display library from a SIV-infected rhesus macaque (Macaca mulatta) using oligonucleotide primers based on human IG gene sequences. Our previous screening relied on Sanger sequencing, which was inefficient and generated only a few dozen sequences. Here, we re-analyzed this library using single molecule, real-time (SMRT) sequencing on the Pacific Biosciences (PacBio) platform to generate thousands of highly accurate circular consensus sequencing (CCS) reads corresponding to full length single chain fragment variable. CCS data were then analyzed through the international ImMunoGeneTics information system^® (IMGT^®)/HighV-QUEST (www.imgt.org) to identify variable genes and perform statistical analyses. Overall the library was very diverse, with 2,569 different IMGT clonotypes called for the 5,238 IGHV sequences assigned to an IMGT clonotype. Within the library, SIV-specific antibodies represented a relatively limited number of clones, with only 135 different IMGT clonotypes called from 4,594 IGHV-assigned sequences. Our data did confirm that the IGHV4 and IGHV3 gene usage was the most abundant within the rhesus antibodies screened, and that these genes were even more enriched among SIV gp140-specific antibodies. Although a broad range of VH CDR3 amino acid (AA) lengths was observed in the unpanned library, the vast majority of SIV gp140-specific antibodies demonstrated a more uniform VH CDR3 length (20 AA). This uniformity was far less apparent when VH CDR3 were classified according to their clonotype (range: 9–25 AA), which we believe is more relevant for specific antibody identification. Only 174 IGKV and 588 IGLV clonotypes were identified within the VL sequences associated with SIV gp140-specific VH. Together, these data strongly suggest that the combination of SMRT sequencing with the IMGT/HighV-QUEST querying tool will facilitate and expedite our understanding of polyclonal antibody responses during SIV infection and may serve to rapidly expand the known scope of macaque V genes utilized during these responses.

HIV Envelope Glycoform Heterogeneity and Localized Diversity Govern the Initiation and Maturation of a V2 Apex Broadly Neutralizing Antibody Lineage

Understanding how broadly neutralizing antibodies (bnAbs) to HIV envelope (Env) develop during natural infection can help guide the rational design of an HIV vaccine. Here, we described a bnAb lineage targeting the Env V2 apex and the Ab-Env co-evolution that led to development of neutralization breadth. The lineage Abs bore an anionic heavy chain complementarity-determining region 3 (CDRH3) of 25 amino acids, among the shortest known for this class of Abs, and achieved breadth with only 10% nucleotide somatic hypermutation and no insertions or deletions. The data suggested a role for Env glycoform heterogeneity in the activation of the lineage germline B cell. Finally, we showed that localized diversity at key V2 epitope residues drove bnAb maturation toward breadth, mirroring the Env evolution pattern described for another donor who developed V2-apex targeting bnAbs. Overall, these findings suggest potential strategies for vaccine approaches based on germline-targeting and serial immunogen design.

Characterization of Hepatitis C Virus (HCV) Envelope Diversification from Acute to Chronic Infection within a Sexually Transmitted HCV Cluster by Using Single-Molecule, Real-Time Sequencing

In contrast to other available next-generation sequencing platforms, PacBio single-molecule, real-time (SMRT) sequencing has the advantage of generating long reads albeit with a relatively higher error rate in unprocessed data. Using this platform, we longitudinally sampled and sequenced the hepatitis C virus (HCV) envelope genome region (1,680 nucleotides [nt]) from individuals belonging to a cluster of sexually transmitted cases. All five subjects were coinfected with HIV-1 and a closely related strain of HCV genotype 4d. In total, 50 samples were analyzed by using SMRT sequencing. By using 7 passes of circular consensus sequencing, the error rate was reduced to 0.37%, and the median number of sequences was 612 per sample. A further reduction of insertions was achieved by alignment against a sample-specific reference sequence. However, in vitro recombination during PCR amplification could not be excluded. Phylogenetic analysis supported close relationships among HCV sequences from the four male subjects and subsequent transmission from one subject to his female partner. Transmission was characterized by a strong genetic bottleneck. Viral genetic diversity was low during acute infection and increased upon progression to chronicity but subsequently fluctuated during chronic infection, caused by the alternate detection of distinct coexisting lineages. SMRT sequencing combines long reads with sufficient depth for many phylogenetic analyses and can therefore provide insights into within-host HCV evolutionary dynamics without the need for haplotype reconstruction using statistical algorithms.IMPORTANCE Next-generation sequencing has revolutionized the study of genetically variable RNA virus populations, but for phylogenetic and evolutionary analyses, longer sequences than those generated by most available platforms, while minimizing the intrinsic error rate, are desired. Here, we demonstrate for the first time that PacBio SMRT sequencing technology can be used to generate full-length HCV envelope sequences at the single-molecule level, providing a data set with large sequencing depth for the characterization of intrahost viral dynamics. The selection of consensus reads derived from at least 7 full circular consensus sequencing rounds significantly reduced the intrinsic high error rate of this method. We used this method to genetically characterize a unique transmission cluster of sexually transmitted HCV infections, providing insight into the distinct evolutionary pathways in each patient over time and identifying the transmission-associated genetic bottleneck as well as fluctuations in viral genetic diversity over time, accompanied by dynamic shifts in viral subpopulations.

Antibody 10-1074 suppresses viremia in HIV-1-infected individuals

Monoclonal antibody 10-1074 targets the V3 glycan supersite on the HIV-1 envelope (Env) protein. It is among the most potent anti-HIV-1 neutralizing antibodies isolated so far. Here we report on its safety and activity in 33 individuals who received a single intravenous infusion of the antibody. 10-1074 was well tolerated and had a half-life of 24.0 d in participants without HIV-1 infection and 12.8 d in individuals with HIV-1 infection. Thirteen individuals with viremia received the highest dose of 30 mg/kg 10-1074. Eleven of these participants were 10-1074-sensitive and showed a rapid decline in viremia by a mean of 1.52 log10 copies/ml. Virologic analysis revealed the emergence of multiple independent 10-1074-resistant viruses in the first weeks after infusion. Emerging escape variants were generally resistant to the related V3-specific antibody PGT121, but remained sensitive to antibodies targeting nonoverlapping epitopes, such as the anti-CD4-binding-site antibodies 3BNC117 and VRC01. The results demonstrate the safety and activity of 10-1074 in humans and support the idea that antibodies targeting the V3 glycan supersite might be useful for the treatment and prevention of HIV-1 infection.

Error rates, PCR recombination, and sampling depth in HIV-1 whole genome deep sequencing

Deep sequencing is a powerful and cost-effective tool to characterize the genetic diversity and evolution of virus populations. While modern sequencing instruments readily cover viral genomes many thousand fold and very rare variants can in principle be detected, sequencing errors, amplification biases, and other artifacts can limit sensitivity and complicate data interpretation. For this reason, the number of studies using whole genome deep sequencing to characterize viral quasi-species in clinical samples is still limited. We have previously undertaken a large scale whole genome deep sequencing study of HIV-1 populations. Here we discuss the challenges, error profiles, control experiments, and computational test we developed to quantify the accuracy of variant frequency estimation.