Person:

Boutwell, C

Viruses diversify over time within hosts, often undercutting the effectiveness of host defenses and therapeutic interventions. To design successful vaccines and therapeutics, it is critical to better understand viral diversification, including comprehensively characterizing the genetic variants in viral intra-host populations and modeling changes from transmission through the course of infection. Massively parallel sequencing technologies can overcome the cost constraints of older sequencing methods and obtain the high sequence coverage needed to detect rare genetic variants (<1%) within an infected host, and to assay variants without prior knowledge. Critical to interpreting deep sequence data sets is the ability to distinguish biological variants from process errors with high sensitivity and specificity. To address this challenge, we describe V-Phaser, an algorithm able to recognize rare biological variants in mixed populations. V-Phaser uses covariation (i.e. phasing) between observed variants to increase sensitivity and an expectation maximization algorithm that iteratively recalibrates base quality scores to increase specificity. Overall, V-Phaser achieved >97% sensitivity and >97% specificity on control read sets. On data derived from a patient after four years of HIV-1 infection, V-Phaser detected 2,015 variants across the ∼10 kb genome, including 603 rare variants (<1% frequency) detected only using phase information. V-Phaser identified variants at frequencies down to 0.2%, comparable to the detection threshold of allele-specific PCR, a method that requires prior knowledge of the variants. The high sensitivity and specificity of V-Phaser enables identifying and tracking changes in low frequency variants in mixed populations such as RNA viruses.

Deep sequencing technologies have the potential to transform the study of highly variable viral pathogens by providing a rapid and cost-effective approach to sensitively characterize rapidly evolving viral quasispecies. Here, we report on a high-throughput whole HIV-1 genome deep sequencing platform that combines 454 pyrosequencing with novel assembly and variant detection algorithms. In one subject we combined these genetic data with detailed immunological analyses to comprehensively evaluate viral evolution and immune escape during the acute phase of HIV-1 infection. The majority of early, low frequency mutations represented viral adaptation to host CD8+ T cell responses, evidence of strong immune selection pressure occurring during the early decline from peak viremia. CD8+ T cell responses capable of recognizing these low frequency escape variants coincided with the selection and evolution of more effective secondary HLA-anchor escape mutations. Frequent, and in some cases rapid, reversion of transmitted mutations was also observed across the viral genome. When located within restricted CD8 epitopes these low frequency reverting mutations were sufficient to prime de novo responses to these epitopes, again illustrating the capacity of the immune response to recognize and respond to low frequency variants. More importantly, rapid viral escape from the most immunodominant CD8+ T cell responses coincided with plateauing of the initial viral load decline in this subject, suggestive of a potential link between maintenance of effective, dominant CD8 responses and the degree of early viremia reduction. We conclude that the early control of HIV-1 replication by immunodominant CD8+ T cell responses may be substantially influenced by rapid, low frequency viral adaptations not detected by conventional sequencing approaches, which warrants further investigation. These data support the critical need for vaccine-induced CD8+ T cell responses to target more highly constrained regions of the virus in order to ensure the maintenance of immunodominant CD8 responses and the sustained decline of early viremia.

Genome editing via CRISPR/Cas9 has rapidly become the tool of choice by virtue of its efficacy and ease of use. However, CRISPR/Cas9-mediated genome editing in clinically relevant human somatic cells remains untested. Here, we report CRISPR/Cas9 targeting of two clinically relevant genes, B2M and CCR5, in primary human CD4+ T cells and CD34+ hematopoietic stem and progenitor cells (HSPCs). Use of single RNA guides led to highly efficient mutagenesis in HSPCs but not in T cells. A dual guide approach improved gene deletion efficacy in both cell types. HSPCs that had undergone genome editing with CRISPR/Cas9 retained multilineage potential. We examined predicted on- and off-target mutations via target capture sequencing in HSPCs and observed low levels of off-target mutagenesis at only one site. These results demonstrate that CRISPR/Cas9 can efficiently ablate genes in HSPCs with minimal off-target mutagenesis, which could have broad applicability for hematopoietic cell-based therapy.

Cytotoxic-T-lymphocyte (CTL) escape mutations in human immunodeficiency viruses encode amino acid substitutions in positions that disrupt CTL targeting, thereby increasing virus survival and conferring a relative fitness benefit. However, it is now clear that CTL escape mutations can also confer a fitness cost, and there is increasing evidence to suggest that in some cases, e.g., escape from HLA-B57/B5801-restricted responses, the costs to the escape virus may affect the clinical course of infection. To quantify the magnitude of the costs of HLA-B57/B5801 escape, a highly sensitive dual-infection assay that uses synonymous nucleotide sequence tags to quantify viral relative replication capacity (RRC) was developed. We then asked whether such CTL escape mutations had an impact equivalent to that seen for a benchmark mutation, the M184V antiretroviral drug resistance mutation of reverse transcriptase (RRC(V184) = 0.86). To answer the question, the RRCs were quantified for escape mutations in three immunodominant HLA-B57/B5801 epitopes in capsid: A146P in IW9 (RRC(P146) = 0.91), A163G in KF11 (RRC(G163) = 0.89), and T242N in TW10 (RRC(N242) = 0.86). Individually, the impact of the escape mutations on RRC was comparable to that of M184V, while coexpression of the mutations resulted in substantial further reductions, with the maximum impact observed for the triple mutant (RRC(P146-G163-N242) = 0.62). By comparison to M184V, the magnitude of the reductions in RRC caused by the escape mutations, particularly when coexpressed, suggests that the costs of escape are sufficient to affect in vivo viral dynamics and may thus play a role in the protective effect associated with HLA-B57/B5801.