Accurately measuring recombination between closely related HIV-1 genomes

Deep mutational scanning reveals functional constraints and antibody-escape potential of Lassa virus glycoprotein complex

Lassa virus is estimated to cause thousands of human deaths per year, primarily due to spillovers from its natural host, Mastomys rodents. Efforts to create vaccines and antibody therapeutics must account for the evolutionary variability of the Lassa virus's glycoprotein complex (GPC), which mediates viral entry into cells and is the target of neutralizing antibodies. To map the evolutionary space accessible to GPC, we used pseudovirus deep mutational scanning to measure how nearly all GPC amino-acid mutations affected cell entry and antibody neutralization. Our experiments defined functional constraints throughout GPC. We quantified how GPC mutations affected neutralization with a panel of monoclonal antibodies. All antibodies tested were escaped by mutations that existed among natural Lassa virus lineages. Overall, our work describes a biosafety-level-2 method to elucidate the mutational space accessible to GPC and shows how prospective characterization of antigenic variation could aid the design of therapeutics and vaccines.

CROPseq-multi: a versatile solution for multiplexed perturbation and decoding in pooled CRISPR screens

Forward genetic screens seek to dissect complex biological systems by systematically perturbing genetic elements and observing the resulting phenotypes. While standard screening methodologies introduce individual perturbations, multiplexing perturbations improves the performance of single-target screens and enables combinatorial screens for the study of genetic interactions. Current tools for multiplexing perturbations are incompatible with pooled screening methodologies that require mRNA-embedded barcodes, including some microscopy and single cell sequencing approaches. Here, we report the development of CROPseq-multi, a CROPseq1-inspired lentiviral system to multiplex Streptococcus pyogenes (Sp) Cas9-based perturbations with mRNA-embedded barcodes. CROPseq-multi has equivalent per-guide activity to CROPseq and low lentiviral recombination frequencies. CROPseq-multi is compatible with enrichment screening methodologies and optical pooled screens, and is extensible to screens with single-cell sequencing readouts. For optical pooled screens, an optimized and multiplexed in situ detection protocol improves barcode detection efficiency 10-fold, enables detection of recombination events, and increases decoding efficiency 3-fold relative to CROPseq. CROPseq-multi is a widely applicable multiplexing solution for diverse SpCas9-based genetic screening approaches.

Deep mutational scanning reveals functional constraints and antigenic variability of Lassa virus glycoprotein complex

Lassa virus is estimated to cause thousands of human deaths per year, primarily due to spillovers from its natural host, Mastomys rodents. Efforts to create vaccines and antibody therapeutics must account for the evolutionary variability of Lassa virus's glycoprotein complex (GPC), which mediates viral entry into cells and is the target of neutralizing antibodies. To map the evolutionary space accessible to GPC, we use pseudovirus deep mutational scanning to measure how nearly all GPC amino-acid mutations affect cell entry and antibody neutralization. Our experiments define functional constraints throughout GPC. We quantify how GPC mutations affect neutralization by a panel of monoclonal antibodies and show that all antibodies are escaped by mutations that exist among natural Lassa virus lineages. Overall, our work describes a biosafety-level-2 method to elucidate the mutational space accessible to GPC and shows how prospective characterization of antigenic variation could aid design of therapeutics and vaccines.

Elevated HIV Viral Load is Associated with Higher Recombination Rate In Vivo

HIV's exceptionally high recombination rate drives its intrahost diversification, enabling immune escape and multidrug resistance within people living with HIV. While we know that HIV's recombination rate varies by genomic position, we have little understanding of how recombination varies throughout infection or between individuals as a function of the rate of cellular coinfection. We hypothesize that denser intrahost populations may have higher rates of coinfection and therefore recombination. To test this hypothesis, we develop a new approach (recombination analysis via time series linkage decay or RATS-LD) to quantify recombination using autocorrelation of linkage between mutations across time points. We validate RATS-LD on simulated data under short read sequencing conditions and then apply it to longitudinal, high-throughput intrahost viral sequencing data, stratifying populations by viral load (a proxy for density). Among sampled viral populations with the lowest viral loads (<26,800 copies/mL), we estimate a recombination rate of 1.5×10-5 events/bp/generation (95% CI: 7×10-6 to 2.9×10-5), similar to existing estimates. However, among samples with the highest viral loads (>82,000 copies/mL), our median estimate is approximately 6 times higher. In addition to co-varying across individuals, we also find that recombination rate and viral load are associated within single individuals across different time points. Our findings suggest that rather than acting as a constant, uniform force, recombination can vary dynamically and drastically across intrahost viral populations and within them over time. More broadly, we hypothesize that this phenomenon may affect other facultatively asexual populations where spatial co-localization varies.

Elevated HIV viral load is associated with higher recombination rate in vivo

HIV's exceptionally high recombination rate drives its intra-host diversification, enabling immune escape and multi-drug resistance within people living with HIV. While we know that HIV's recombination rate varies by genomic position, we have little understanding of how recombination varies throughout infection or between individuals as a function of the rate of cellular coinfection. We hypothesize that denser intra-host populations may have higher rates of coinfection and therefore recombination. To test this hypothesis, we develop a new approach (Recombination Analysis via Time Series Linkage Decay, or RATS-LD) to quantify recombination using autocorrelation of linkage between mutations across time points. We validate RATS-LD on simulated data under short read sequencing conditions and then apply it to longitudinal, high-throughput intra-host viral sequencing data, stratifying populations by viral load (a proxy for density). Among sampled viral populations with the lowest viral loads (< 26,800 copies/mL), we estimate a recombination rate of 1.5 × 10-5 events/bp/generation (95% CI: 7 × 10-6 - 2.9 × 10-5), similar to existing estimates. However, among samples with the highest viral loads (> 82,000 copies/mL), our median estimate is approximately 6 times higher. In addition to co-varying across individuals, we also find that recombination rate and viral load are associated within single individuals across different time points. Our findings suggest that rather than acting as a constant, uniform force, recombination can vary dynamically and drastically across intra-host viral populations and within them over time. More broadly, we hypothesize that this phenomenon may affect other facultatively asexual populations where spatial co-localization varies.

HIV drug resistance in the era of contemporary antiretroviral therapy: A clinical perspective

Contemporary antiretroviral therapy (ART) regimens have high barriers to the development of drug resistance. However, resistance to earlier antiretrovirals and uncommon cases of resistance to contemporary ART illustrate the continued need for good clinical management of HIV drug resistance. Here, we describe HIV drug-resistance mechanisms, the interaction of HIV drug-resistant mutations and the patterns of drug resistance to contemporary ART. We then provide guidance on the management of HIV drug resistance, including how to limit the development of resistance and manage virologic failure that is complicated by resistance. To complement this, links to resources and treatment guidelines are provided that can assist with the interpretation of HIV drug resistance test results and optimal ART selection in the clinic.

Mapping the neutralizing specificity of human anti-HIV serum by deep mutational scanning

Understanding the specificities of human serum antibodies that broadly neutralize HIV can inform prevention and treatment strategies. Here, we describe a deep mutational scanning system that can measure the effects of combinations of mutations to HIV envelope (Env) on neutralization by antibodies and polyclonal serum. We first show that this system can accurately map how all functionally tolerated mutations to Env affect neutralization by monoclonal antibodies. We then comprehensively map Env mutations that affect neutralization by a set of human polyclonal sera that neutralize diverse strains of HIV and target the site engaging the host receptor CD4. The neutralizing activities of these sera target different epitopes, with most sera having specificities reminiscent of individual characterized monoclonal antibodies, but one serum targeting two epitopes within the CD4-binding site. Mapping the specificity of the neutralizing activity in polyclonal human serum will aid in assessing anti-HIV immune responses to inform prevention strategies.

Massively parallel evaluation and computational prediction of the activities and specificities of 17 small Cas9s

Recently, various small Cas9 orthologs and variants have been reported for use in in vivo delivery applications. Although small Cas9s are particularly suited for this purpose, selecting the most optimal small Cas9 for use at a specific target sequence continues to be challenging. Here, to this end, we have systematically compared the activities of 17 small Cas9s for thousands of target sequences. For each small Cas9, we have characterized the protospacer adjacent motif and determined optimal single guide RNA expression formats and scaffold sequence. High-throughput comparative analyses revealed distinct high- and low-activity groups of small Cas9s. We also developed DeepSmallCas9, a set of computational models predicting the activities of the small Cas9s at matched and mismatched target sequences. Together, this analysis and these computational models provide a useful guide for researchers to select the most suitable small Cas9 for specific applications.

Mapping the neutralizing specificity of human anti-HIV serum by deep mutational scanning

Understanding the specificities of human serum antibodies that broadly neutralize HIV can inform prevention and treatment strategies. Here we describe a deep mutational scanning system that can measure the effects of combinations of mutations to HIV envelope (Env) on neutralization by antibodies and polyclonal serum. We first show that this system can accurately map how all functionally tolerated mutations to Env affect neutralization by monoclonal antibodies. We then comprehensively map Env mutations that affect neutralization by a set of human polyclonal sera known to target the CD4-binding site that neutralize diverse strains of HIV. The neutralizing activities of these sera target different epitopes, with most sera having specificities reminiscent of individual characterized monoclonal antibodies, but one sera targeting two epitopes within the CD4 binding site. Mapping the specificity of the neutralizing activity in polyclonal human serum will aid in assessing anti-HIV immune responses to inform prevention strategies.

A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spike

A major challenge in understanding SARS-CoV-2 evolution is interpreting the antigenic and functional effects of emerging mutations in the viral spike protein. Here, we describe a deep mutational scanning platform based on non-replicative pseudotyped lentiviruses that directly quantifies how large numbers of spike mutations impact antibody neutralization and pseudovirus infection. We apply this platform to produce libraries of the Omicron BA.1 and Delta spikes. These libraries each contain ∼7,000 distinct amino acid mutations in the context of up to ∼135,000 unique mutation combinations. We use these libraries to map escape mutations from neutralizing antibodies targeting the receptor-binding domain, N-terminal domain, and S2 subunit of spike. Overall, this work establishes a high-throughput and safe approach to measure how ∼105 combinations of mutations affect antibody neutralization and spike-mediated infection. Notably, the platform described here can be extended to the entry proteins of many other viruses.

A pseudovirus system enables deep mutational scanning of the full SARS-CoV-2 spike

A major challenge in understanding SARS-CoV-2 evolution is interpreting the antigenic and functional effects of emerging mutations in the viral spike protein. Here we describe a new deep mutational scanning platform based on non-replicative pseudotyped lentiviruses that directly quantifies how large numbers of spike mutations impact antibody neutralization and pseudovirus infection. We demonstrate this new platform by making libraries of the Omicron BA.1 and Delta spikes. These libraries each contain ~7000 distinct amino-acid mutations in the context of up to ~135,000 unique mutation combinations. We use these libraries to map escape mutations from neutralizing antibodies targeting the receptor binding domain, N-terminal domain, and S2 subunit of spike. Overall, this work establishes a high-throughput and safe approach to measure how ~10 5 combinations of mutations affect antibody neutralization and spike-mediated infection. Notably, the platform described here can be extended to the entry proteins of many other viruses.

Detection of SARS-CoV-2 intra-host recombination during superinfection with Alpha and Epsilon variants in New York City

Recombination is an evolutionary process by which many pathogens generate diversity and acquire novel functions. Although a common occurrence during coronavirus replication, detection of recombination is only feasible when genetically distinct viruses contemporaneously infect the same host. Here, we identify an instance of SARS-CoV-2 superinfection, whereby an individual was infected with two distinct viral variants: Alpha (B.1.1.7) and Epsilon (B.1.429). This superinfection was first noted when an Alpha genome sequence failed to exhibit the classic S gene target failure behavior used to track this variant. Full genome sequencing from four independent extracts reveals that Alpha variant alleles comprise around 75% of the genomes, whereas the Epsilon variant alleles comprise around 20% of the sample. Further investigation reveals the presence of numerous recombinant haplotypes spanning the genome, specifically in the spike, nucleocapsid, and ORF 8 coding regions. These findings support the potential for recombination to reshape SARS-CoV-2 genetic diversity.

RNA Structures and Their Role in Selective Genome Packaging

To generate infectious viral particles, viruses must specifically select their genomic RNA from milieu that contains a complex mixture of cellular or non-genomic viral RNAs. In this review, we focus on the role of viral encoded RNA structures in genome packaging. We first discuss how packaging signals are constructed from local and long-range base pairings within viral genomes, as well as inter-molecular interactions between viral and host RNAs. Then, how genome packaging is regulated by the biophysical properties of RNA. Finally, we examine the impact of RNA packaging signals on viral evolution.

Adeno-Associated Viral Vector Mediated Expression of Broadly- Neutralizing Antibodies Against HIV-Hitting a Fast-Moving Target

The vast genetic variability of HIV has impeded efforts towards a cure for HIV. Lifelong administration of combined antiretroviral therapy (cART) is highly effective against HIV and has markedly increased the life expectancy of HIV infected individuals. However, the long-term usage of cART is associated with co-morbidities and the emergence of multidrug-resistant escape mutants necessitating the development of alternative approaches to combat HIV/AIDS. In the past decade, the development of single-cell antibody cloning methods has facilitated the characterization of a diverse array of highly potent neutralizing antibodies against a broad range of HIV strains. Although the passive transfer of these broadly neutralizing antibodies (bnAbs) in both animal models and humans has been shown to elicit significant antiviral effects, long term virologic suppression requires repeated administration of these antibodies. Adeno-associated virus (AAV) mediated antibody gene transfer provides a long-term expression of these antibodies from a single administration of the recombinant vector. Therefore, this vectored approach holds promises in the treatment and prevention of a chronic disease like HIV infection. Here, we provide an overview of HIV genetic diversity, AAV vectorology, and anti-HIV bnAbs and summarize the promises and challenges of the application of AAV in the delivery of bnAbs for HIV prevention and therapy.

Predicting the efficiency of prime editing guide RNAs in human cells

Prime editing enables the introduction of virtually any small-sized genetic change without requiring donor DNA or double-strand breaks. However, evaluation of prime editing efficiency requires time-consuming experiments, and the factors that affect efficiency have not been extensively investigated. In this study, we performed high-throughput evaluation of prime editor 2 (PE2) activities in human cells using 54,836 pairs of prime editing guide RNAs (pegRNAs) and their target sequences. The resulting data sets allowed us to identify factors affecting PE2 efficiency and to develop three computational models to predict pegRNA efficiency. For a given target sequence, the computational models predict efficiencies of pegRNAs with different lengths of primer binding sites and reverse transcriptase templates for edits of various types and positions. Testing the accuracy of the predictions using test data sets that were not used for training, we found Spearman's correlations between 0.47 and 0.81. Our computational models and information about factors affecting PE2 efficiency will facilitate practical application of prime editing.

Prediction of the sequence-specific cleavage activity of Cas9 variants

Several Streptococcus pyogenes Cas9 (SpCas9) variants have been developed to improve an enzyme's specificity or to alter or broaden its protospacer-adjacent motif (PAM) compatibility, but selecting the optimal variant for a given target sequence and application remains difficult. To build computational models to predict the sequence-specific activity of 13 SpCas9 variants, we first assessed their cleavage efficiency at 26,891 target sequences. We found that, of the 256 possible four-nucleotide NNNN sequences, 156 can be used as a PAM by at least one of the SpCas9 variants. For the high-fidelity variants, overall activity could be ranked as SpCas9 ≥ Sniper-Cas9 > eSpCas9(1.1) > SpCas9-HF1 > HypaCas9 ≈ xCas9 >> evoCas9, whereas their overall specificities could be ranked as evoCas9 >> HypaCas9 ≥ SpCas9-HF1 ≈ eSpCas9(1.1) > xCas9 > Sniper-Cas9 > SpCas9. Using these data, we developed 16 deep-learning-based computational models that accurately predict the activity of these variants at any target sequence.

Imaging-based pooled CRISPR screening reveals regulators of lncRNA localization

Pooled-library CRISPR screening provides a powerful means to discover genetic factors involved in cellular processes in a high-throughput manner. However, the phenotypes accessible to pooled-library screening are limited. Complex phenotypes, such as cellular morphology and subcellular molecular organization, as well as their dynamics, require imaging-based readout and are currently beyond the reach of pooled-library CRISPR screening. Here we report an all imaging-based pooled-library CRISPR screening approach that combines high-content phenotype imaging with high-throughput single guide RNA (sgRNA) identification in individual cells. In this approach, sgRNAs are codelivered to cells with corresponding barcodes placed at the 3' untranslated region of a reporter gene using a lentiviral delivery system with reduced recombination-induced sgRNA-barcode mispairing. Multiplexed error-robust fluorescence in situ hybridization (MERFISH) is used to read out the barcodes and hence identify the sgRNAs with high accuracy. We used this approach to screen 162 sgRNAs targeting 54 RNA-binding proteins for their effects on RNA localization to nuclear compartments and uncovered previously unknown regulatory factors for nuclear RNA localization. Notably, our screen revealed both positive and negative regulators for the nuclear speckle localization of a long noncoding RNA, MALAT1, suggesting a dynamic regulation of lncRNA localization in subcellular compartments.

Dynamics of HIV-1 coinfection in different susceptible target cell populations during cell-free infection

HIV-1 mutations rapidly accumulate through genetic recombination events, which require the infection of a single cell by two virions (coinfection). Accumulation of mutations in the viral population may lead to immune escape and high-level drug resistance. The existence of cell subpopulations characterized by different susceptibility to HIV-1 infection has been proposed as an important parameter driving coinfection (Dang et al., 2004). While the mechanism and the quantification of HIV-1 coinfection have been recently investigated by mathematical models, the detailed dynamics of this process during cell-free infection remains elusive. In this study, we constructed ordinary differential equations considering the heterogeneity of target cell populations during cell-free infection in cell culture, and reproduced the cell culture experimental data. Our mathematical analyses showed that the presence of two differently susceptible target cell subpopulations could explain our experimental datasets, while increasing the number of subpopulations did not improve the fitting. In addition, we quantitatively demonstrated that cells infected by multiple viruses mainly accumulated from one cell subpopulation under cell-free infection conditions. In particular, the frequency of infection events in the more susceptible subpopulation was 6.11-higher than that from the other subpopulation, and 98.3% of coinfected cells emerged from the more susceptible subpopulation. Our mathematical-experimental approach is able to extract such a quantitative information, and can be easily applied to other virus infections.

Frequent sgRNA-barcode recombination in single-cell perturbation assays

Simultaneously detecting CRISPR-based perturbations and induced transcriptional changes in the same cell is a powerful approach to unraveling genome function. Several lentiviral approaches have been developed, some of which rely on the detection of distally located genetic barcodes as an indirect proxy of sgRNA identity. Since barcodes are often several kilobases from their corresponding sgRNAs, viral recombination-mediated swapping of barcodes and sgRNAs is feasible. Using a self-circularization-based sgRNA-barcode library preparation protocol, we estimate the recombination rate to be ~50% and we trace this phenomenon to the pooled viral packaging step. Recombination is random, and decreases the signal-to-noise ratio of the assay. Our results suggest that alternative approaches can increase the throughput and sensitivity of single-cell perturbation assays.

Molecular Epidemiology of HIV-1 Virus in Puerto Rico: Novel Cases of HIV-1 Subtype C, D, and CRF-24BG

HIV-1 subtype B virus is the most prevalent subtype in Puerto Rico (PR), accounting for about 90% of infection in the island. Recently, other subtypes and circulating recombinant forms (CRFs), including F(12_BF), A (01_BF), and CRF-39 BF-like, have been identified. The purpose of this study is to assess the distribution of drug resistance mutations and subtypes in PR. A total of 846 nucleotide sequences from the period comprising 2013 through 2017 were obtained from our "HIV Genotyping" test file. Phylogenetic and molecular epidemiology analyses were performed to evaluate the evolutionary dynamics and prevalence of drug resistance mutations. According to our results, we detected a decrease in the prevalence of protease inhibitor, nucleoside reverse transcriptase inhibitor (NRTI), and non-NRTI (NNRTI) resistance mutations over time. In addition, we also detected recombinant forms and, for the first time, identified subtypes C, D, and CRF-24BG in PR. Recent studies suggest that non-subtypes B are associated with a high risk of treatment failure and disease progression. The constant monitoring of viral evolution and drug resistance mutation dynamics is important to establish appropriate efforts for controlling viral expansion.

Neutral Theory and Rapidly Evolving Viral Pathogens

The evolution of viral pathogens is shaped by strong selective forces that are exerted during jumps to new hosts, confrontations with host immune responses and antiviral drugs, and numerous other processes. However, while undeniably strong and frequent, adaptive evolution is largely confined to small parts of information-packed viral genomes, and the majority of observed variation is effectively neutral. The predictions and implications of the neutral theory have proven immensely useful in this context, with applications spanning understanding within-host population structure, tracing the origins and spread of viral pathogens, predicting evolutionary dynamics, and modeling the emergence of drug resistance. We highlight the multiple ways in which the neutral theory has had an impact, which has been accelerated in the age of high-throughput, high-resolution genomics.

Tracking HIV-1 recombination to resolve its contribution to HIV-1 evolution in natural infection

Recombination in HIV-1 is well documented, but its importance in the low-diversity setting of within-host diversification is less understood. Here we develop a novel computational tool (RAPR (Recombination Analysis PRogram)) to enable a detailed view of in vivo viral recombination during early infection, and we apply it to near-full-length HIV-1 genome sequences from longitudinal samples. Recombinant genomes rapidly replace transmitted/founder (T/F) lineages, with a median half-time of 27 days, increasing the genetic complexity of the viral population. We identify recombination hot and cold spots that differ from those observed in inter-subtype recombinants. Furthermore, RAPR analysis of longitudinal samples from an individual with well-characterized neutralizing antibody responses shows that recombination helps carry forward resistance-conferring mutations in the diversifying quasispecies. These findings provide insight into molecular mechanisms by which viral recombination contributes to HIV-1 persistence and immunopathogenesis and have implications for studies of HIV transmission and evolution in vivo.

On the design of CRISPR-based single-cell molecular screens

Several groups recently coupled CRISPR perturbations and single-cell RNA-seq for pooled genetic screens. We demonstrate that vector designs of these studies are susceptible to ∼50% swapping of guide RNA-barcode associations because of lentiviral template switching. We optimized a published alternative, CROP-seq, in which the guide RNA also serves as the barcode, and here confirm that this strategy performs robustly and doubled the rate at which guides are assigned to cells to 94%.

Retroviral RNA Dimerization: From Structure to Functions

The genome of the retroviruses is a dimer composed by two homologous copies of genomic RNA (gRNA) molecules of positive polarity. The dimerization process allows two gRNA molecules to be non-covalently linked together through intermolecular base-pairing. This step is critical for the viral life cycle and is highly conserved among retroviruses with the exception of spumaretroviruses. Furthermore, packaging of two gRNA copies into viral particles presents an important evolutionary advantage for immune system evasion and drug resistance. Recent studies reported RNA switches models regulating not only gRNA dimerization, but also translation and packaging, and a spatio-temporal characterization of viral gRNA dimerization within cells are now at hand. This review summarizes our current understanding on the structural features of the dimerization signals for a variety of retroviruses (HIVs, MLV, RSV, BLV, MMTV, MPMV…), the mechanisms of RNA dimer formation and functional implications in the retroviral cycle.

A platform for functional assessment of large variant libraries in mammalian cells

Sequencing-based, massively parallel genetic assays have revolutionized our ability to quantify the relationship between many genotypes and a phenotype of interest. Unfortunately, variant library expression platforms in mammalian cells are far from ideal, hindering the study of human gene variants in their physiologically relevant cellular contexts. Here, we describe a platform for phenotyping variant libraries in transfectable mammalian cell lines in two steps. First, a landing pad cell line with a genomically integrated, Tet-inducible cassette containing a Bxb1 recombination site is created. Second, a single variant from a library of transfected, promoter-less plasmids is recombined into the landing pad in each cell. Thus, every cell in the recombined pool expresses a single variant, allowing for parallel, sequencing-based assessment of variant effect. We describe a method for incorporating a single landing pad into a defined site of a cell line of interest, and show that our approach can be used generate more than 20 000 recombinant cells in a single experiment. Finally, we use our platform in combination with a sequencing-based assay to explore the N-end rule by simultaneously measuring the effects of all possible N-terminal amino acids on protein expression.

Effect of the Latent Reservoir on the Evolution of HIV at the Within- and Between-Host Levels

The existence of long-lived reservoirs of latently infected CD4+ T cells is the major barrier to curing HIV, and has been extensively studied in this light. However, the effect of these reservoirs on the evolutionary dynamics of the virus has received little attention. Here, we present a within-host quasispecies model that incorporates a long-lived reservoir, which we then nest into an epidemiological model of HIV dynamics. For biologically plausible parameter values, we find that the presence of a latent reservoir can severely delay evolutionary dynamics within a single host, with longer delays associated with larger relative reservoir sizes and/or homeostatic proliferation of cells within the reservoir. These delays can fundamentally change the dynamics of the virus at the epidemiological scale. In particular, the delay in within-host evolutionary dynamics can be sufficient for the virus to evolve intermediate viral loads consistent with maximising transmission, as is observed, and not the very high viral loads that previous models have predicted, an effect that can be further enhanced if viruses similar to those that initiate infection are preferentially transmitted. These results depend strongly on within-host characteristics such as the relative reservoir size, with the evolution of intermediate viral loads observed only when the within-host dynamics are sufficiently delayed. In conclusion, we argue that the latent reservoir has important, and hitherto under-appreciated, roles in both within- and between-host viral evolution.

Viral reverse transcriptases

Reverse transcriptases (RTs) play a major role in the replication of Retroviridae, Metaviridae, Pseudoviridae, Hepadnaviridae and Caulimoviridae. RTs are enzymes that are able to synthesize DNA using RNA or DNA as templates (DNA polymerase activity), and degrade RNA when forming RNA/DNA hybrids (ribonuclease H activity). In retroviruses and LTR retrotransposons (Metaviridae and Pseudoviridae), the coordinated action of both enzymatic activities converts single-stranded RNA into a double-stranded DNA that is flanked by identical sequences known as long terminal repeats (LTRs). RTs of retroviruses and LTR retrotransposons are active as monomers (e.g. murine leukemia virus RT), homodimers (e.g. Ty3 RT) or heterodimers (e.g. human immunodeficiency virus type 1 (HIV-1) RT). RTs lack proofreading activity and display high intrinsic error rates. Besides, high recombination rates observed in retroviruses are promoted by poor processivity that causes template switching, a hallmark of reverse transcription. HIV-1 RT inhibitors acting on its polymerase activity constitute the backbone of current antiretroviral therapies, although novel drugs, including ribonuclease H inhibitors, are still necessary to fight HIV infections. In Hepadnaviridae and Caulimoviridae, reverse transcription leads to the formation of nicked circular DNAs that will be converted into episomal DNA in the host cell nucleus. Structural and biochemical information on their polymerases is limited, although several drugs inhibiting HIV-1 RT are known to be effective against the human hepatitis B virus polymerase. In this review, we summarize current knowledge on reverse transcription in the five virus families and discuss available biochemical and structural information on RTs, including their biosynthesis, enzymatic activities, and potential inhibition.

Antiviral drug resistance as an adaptive process

Antiviral drug resistance is a matter of great clinical importance that, historically, has been investigated mostly from a virological perspective. Although the proximate mechanisms of resistance can be readily uncovered using these methods, larger evolutionary trends often remain elusive. Recent interest by population geneticists in studies of antiviral resistance has spurred new metrics for evaluating mutation and recombination rates, demographic histories of transmission and compartmentalization, and selective forces incurred during viral adaptation to antiviral drug treatment. We present up-to-date summaries on antiviral resistance for a range of drugs and viral types, and review recent advances for studying their evolutionary histories. We conclude that information imparted by demographic and selective histories, as revealed through population genomic inference, is integral to assessing the evolution of antiviral resistance as it pertains to human health.

HIV-1 Mutation and Recombination Rates Are Different in Macrophages and T-cells

High rates of mutation and recombination help human immunodeficiency virus (HIV) to evade the immune system and develop resistance to antiretroviral therapy. Macrophages and T-cells are the natural target cells of HIV-1 infection. A consensus has not been reached as to whether HIV replication results in differential recombination between primary T-cells and macrophages. Here, we used HIV with silent mutation markers along with next generation sequencing to compare the mutation and the recombination rates of HIV directly in T lymphocytes and macrophages. We observed a more than four-fold higher recombination rate of HIV in macrophages compared to T-cells (p < 0.001) and demonstrated that this difference is not due to different reliance on C-X-C chemokine receptor type 4 (CXCR4) and C-C chemokine receptor type 5 (CCR5) co-receptors between T-cells and macrophages. We also found that the pattern of recombination across the HIV genome (hot and cold spots) remains constant between T-cells and macrophages despite a three-fold increase in the overall recombination rate. This indicates that the difference in rates is a general feature of HIV DNA synthesis during macrophage infection. In contrast to HIV recombination, we found that T-cells have a 30% higher mutation rate than macrophages (p < 0.001) and that the mutational profile is similar between these cell types. Unexpectedly, we found no association between mutation and recombination in macrophages, in contrast to T-cells. Our data highlights some of the fundamental difference of HIV recombination and mutation amongst these two major target cells of infection. Understanding these differences will provide invaluable insights toward HIV evolution and how the virus evades immune surveillance and anti-retroviral therapeutics.

Estimating the in-vivo HIV template switching and recombination rate

BACKGROUND

HIV recombination has been estimated in vitro using a variety of approaches, and shows a high rate of template switching per reverse transcription event. In-vivo studies of recombination generally measure the accumulation of recombinant strains over time, and thus do not directly estimate a comparable template switching rate.

METHOD

To examine whether the estimated in-vitro template switching rate is representative of the rate that occurs during HIV infection in vivo, we adopted a novel approach, analysing single genome sequences from early founder viruses to study the in-vivo template switching rate in the env region of HIV.

RESULTS

We estimated the in-vivo per cycle template switching rate to be between 0.5 and 1.5/1000 nt, or approximately 5-14 recombination events over the length of the HIV genome.

CONCLUSION

The in-vivo estimated template switching rate is close to the in-vitro estimated rate found in primary T lymphocytes but not macrophages, which is consistent with the majority of HIV infection occurring in T lymphocytes.

A general method to eliminate laboratory induced recombinants during massive, parallel sequencing of cDNA library

BACKGROUND

Massive, parallel sequencing is a potent tool for dissecting the regulation of biological processes by revealing the dynamics of the cellular RNA profile under different conditions. Similarly, massive, parallel sequencing can be used to reveal the complexity of viral quasispecies that are often found in the RNA virus infected host. However, the production of cDNA libraries for next-generation sequencing (NGS) necessitates the reverse transcription of RNA into cDNA and the amplification of the cDNA template using PCR, which may introduce artefact in the form of phantom nucleic acids species that can bias the composition and interpretation of original RNA profiles.

METHOD

Using HIV as a model we have characterised the major sources of error during the conversion of viral RNA to cDNA, namely excess RNA template and the RNaseH activity of the polymerase enzyme, reverse transcriptase. In addition we have analysed the effect of PCR cycle on detection of recombinants and assessed the contribution of transfection of highly similar plasmid DNA to the formation of recombinant species during the production of our control viruses.

RESULTS

We have identified RNA template concentrations, RNaseH activity of reverse transcriptase, and PCR conditions as key parameters that must be carefully optimised to minimise chimeric artefacts.

CONCLUSIONS

Using our optimised RT-PCR conditions, in combination with our modified PCR amplification procedure, we have developed a reliable technique for accurate determination of RNA species using NGS technology.

Recombination in viruses: mechanisms, methods of study, and evolutionary consequences

Recombination is a pervasive process generating diversity in most viruses. It joins variants that arise independently within the same molecule, creating new opportunities for viruses to overcome selective pressures and to adapt to new environments and hosts. Consequently, the analysis of viral recombination attracts the interest of clinicians, epidemiologists, molecular biologists and evolutionary biologists. In this review we present an overview of three major areas related to viral recombination: (i) the molecular mechanisms that underlie recombination in model viruses, including DNA-viruses (Herpesvirus) and RNA-viruses (Human Influenza Virus and Human Immunodeficiency Virus), (ii) the analytical procedures to detect recombination in viral sequences and to determine the recombination breakpoints, along with the conceptual and methodological tools currently used and a brief overview of the impact of new sequencing technologies on the detection of recombination, and (iii) the major areas in the evolutionary analysis of viral populations on which recombination has an impact. These include the evaluation of selective pressures acting on viral populations, the application of evolutionary reconstructions in the characterization of centralized genes for vaccine design, and the evaluation of linkage disequilibrium and population structure.

Retroviral vectors for analysis of viral mutagenesis and recombination

Retrovirus population diversity within infected hosts is commonly high due in part to elevated rates of replication, mutation, and recombination. This high genetic diversity often complicates the development of effective diagnostics, vaccines, and antiviral drugs. This review highlights the diverse vectors and approaches that have been used to examine mutation and recombination in retroviruses. Retroviral vectors for these purposes can broadly be divided into two categories: those that utilize reporter genes as mutation or recombination targets and those that utilize viral genes as targets of mutation or recombination. Reporter gene vectors greatly facilitate the detection, quantification, and characterization of mutants and/or recombinants, but may not fully recapitulate the patterns of mutagenesis or recombination observed in native viral gene sequences. In contrast, the detection of mutations or recombination events directly in viral genes is more biologically relevant but also typically more challenging and inefficient. We will highlight the advantages and disadvantages of the various vectors and approaches used as well as propose ways in which they could be improved.

Fifteen to twenty percent of HIV substitution mutations are associated with recombination

HIV undergoes high rates of mutation and recombination during reverse transcription, but it is not known whether these events occur independently or are linked mechanistically. Here we used a system of silent marker mutations in HIV and a single round of infection in primary T lymphocytes combined with a high-throughput sequencing and mathematical modeling approach to directly estimate the viral recombination and mutation rates. From >7 million nucleotides (nt) of sequences from HIV infection, we observed 4,801 recombination events and 859 substitution mutations (≈1.51 and 0.12 events per 1,000 nt, respectively). We used experimental controls to account for PCR-induced and transfection-induced recombination and sequencing error. We found that the single-cycle virus-induced mutation rate is 4.6 × 10(-5) mutations per nt after correction. By sorting of our data into recombined and nonrecombined sequences, we found a significantly higher mutation rate in recombined regions (P = 0.003 by Fisher's exact test). We used a permutation approach to eliminate a number of potential confounding factors and confirm that mutation occurs around the site of recombination and is not simply colocated in the genome. By comparing mutation rates in recombined and nonrecombined regions, we found that recombination-associated mutations account for 15 to 20% of all mutations occurring during reverse transcription.

Identifying recombination hot spots in the HIV-1 genome

HIV-1 infection is characterized by the rapid generation of genetic diversity that facilitates viral escape from immune selection and antiretroviral therapy. Despite recombination's crucial role in viral diversity and evolution, little is known about the genomic factors that influence recombination between highly similar genomes. In this study, we use a minimally modified full-length HIV-1 genome and high-throughput sequence analysis to study recombination in gag and pol in T cells. We find that recombination is favored at a number of recombination hot spots, where recombination occurs six times more frequently than at corresponding cold spots. Interestingly, these hot spots occur near important features of the HIV-1 genome but do not occur at sites immediately around protease inhibitor or reverse transcriptase inhibitor drug resistance mutations. We show that the recombination hot and cold spots are consistent across five blood donors and are independent of coreceptor-mediated entry. Finally, we check common experimental confounders and find that these are not driving the location of recombination hot spots. This is the first study to identify the location of recombination hot spots between two similar viral genomes with great statistical power and under conditions that closely reflect natural recombination events among HIV-1 quasispecies.

IMPORTANCE

The ability of HIV-1 to evade the immune system and antiretroviral therapy depends on genetic diversity within the viral quasispecies. Retroviral recombination is an important mechanism that helps to generate and maintain this genetic diversity, but little is known about how recombination rates vary within the HIV-1 genome. We measured recombination rates in gag and pol and identified recombination hot and cold spots, demonstrating that recombination is not random but depends on the underlying gene sequence. The strength and location of these recombination hot and cold spots can be used to improve models of viral dynamics and evolution, which will be useful for the design of robust antiretroviral therapies.

Next-generation sequencing of HIV-1 RNA genomes: determination of error rates and minimizing artificial recombination

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

HIV-1 Env C2-V4 diversification in a slow-progressor infant reveals a flat but rugged fitness landscape

Human immunodeficiency virus type-1 (HIV-1) fitness has been associated with virus entry, a process mediated by the envelope glycoprotein (Env). We previously described Env genetic diversification in a Zambian, subtype C infected, slow-progressor child (1157i) in parallel with an evolving neutralizing antibody response. Because of the role the Variable-3 loop (V3) plays in transmission, cell tropism, neutralization sensitivity, and fitness, longitudinally isolated 1157i C2-V4 alleles were cloned into HIV-1NL4-3-eGFP and -DsRed2 infectious molecular clones. The fluorescent reporters allowed for dual-infection competitions between all patient-derived C2-V4 chimeras to quantify the effect of V3 diversification and selection on fitness. 'Winners' and 'losers' were readily discriminated among the C2-V4 alleles. Exceptional sensitivity for detection of subtle fitness differences was revealed through analysis of two alleles differing in a single synonymous amino acid. However, when the outcomes of N = 33 competitions were averaged for each chimera, the aggregate analysis showed that despite increasing diversification and divergence with time, natural selection of C2-V4 sequences in this individual did not appear to be producing a 'survival of the fittest' evolutionary pattern. Rather, we detected a relatively flat fitness landscape consistent with mutational robustness. Fitness outcomes were then correlated with individual components of the entry process. Env incorporation into particles correlated best with fitness, suggesting a role for Env avidity, as opposed to receptor/coreceptor affinity, in defining fitness. Nevertheless, biochemical analyses did not identify any step in HIV-1 entry as a dominant determinant of fitness. Our results lead us to conclude that multiple aspects of entry contribute to maintaining adequate HIV-1 fitness, and there is no surrogate analysis for determining fitness. The capacity for subtle polymorphisms in Env to nevertheless significantly impact viral fitness suggests fitness is best defined by head-to-head competition.

Suboptimal provirus expression explains apparent nonrandom cell coinfection with HIV-1

Despite the ability of primate lentiviruses to prevent infected cells from being reinfected, cell coinfection has occurred in the past and has shaped virus evolution by promoting the biogenesis of heterozygous virions and recombination during reverse transcription. In vitro experiments have shown that cell coinfection with HIV is more frequent than would be expected if coinfection were a random process. A possible explanation for this bias is the heterogeneity of target cells and the preferred infection of a subpopulation. To address this question, we compared the frequency of double-positive cells measured following coincubation with green fluorescent protein (GFP) and DsRed HIV reporter viruses with that of stochastic coinfection calculated as the product of the frequencies of GFP- and DsRed-positive cells upon incubation with either reporter virus. Coinfection was more frequent than would be expected on the grounds of stochastic infection, due to the underestimation of single-infection frequencies, which mathematically decreased the calculated frequency. Indeed, when cells were incubated with either reporter virus, a fraction of the cells were scored as uninfected yet harbored a silent provirus that was reactivated upon coinfection through cross talk between viral elements. When such cross talk was avoided, experimental and calculated coinfection frequencies matched, indicating random coinfection. The proportion of infected cells harboring a silent provirus was estimated from coinfection experiments and was shown to be cell type dependent but independent of the virus entry route.

The origin of genetic diversity in HIV-1

One of the hallmarks of HIV infection is the rapid development of a genetically complex population (quasispecies) from an initially limited number of infectious particles. Genetic diversity remains one of the major obstacles to eradication of HIV. The viral quasispecies can respond rapidly to selective pressures, such as that imposed by the immune system and antiretroviral therapy, and frustrates vaccine design efforts. Two unique features of retroviral replication are responsible for the unprecedented variation generated during infection. First, mutations are frequently introduced into the viral genome by the error prone viral reverse transcriptase and through the actions of host cellular factors, such as the APOBEC family of nucleic acid editing enzymes. Second, the HIV reverse transcriptase can utilize both copies of the co-packaged viral genome in a process termed retroviral recombination. When the co-packaged viral genomes are genetically different, retroviral recombination can lead to the shuffling of mutations between viral genomes in the quasispecies. This review outlines the stages of the retroviral life cycle where genetic variation is introduced, focusing on the principal mechanisms of mutation and recombination. Understanding the mechanistic origin of genetic diversity is essential to combating HIV.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

The dynamics of HIV-1 adaptation in early infection

Human immunodeficiency virus type 1 (HIV-1) undergoes a severe population bottleneck during sexual transmission and yet adapts extremely rapidly to the earliest immune responses. The bottleneck has been inferred to typically consist of a single genome, and typically eight amino acid mutations in viral proteins spread to fixation by the end of the early chronic phase of infection in response to selection by CD8(+) T cells. Stochastic simulation was used to examine the effects of the transmission bottleneck and of potential interference among spreading immune-escape mutations on the adaptive dynamics of the virus in early infection. If major viral population genetic parameters are assigned realistic values that permit rapid adaptive evolution, then a bottleneck of a single genome is not inconsistent with the observed pattern of adaptive fixations. One requirement is strong selection by CD8(+) T cells that decreases over time. Such selection may reduce effective population sizes at linked loci through genetic hitchhiking. However, this effect is predicted to be minor in early infection because the transmission bottleneck reduces the effective population size to such an extent that the resulting strong selection and weak mutation cause beneficial mutations to fix sequentially and thus avoid interference.

The evolution of HIV: inferences using phylogenetics

Molecular phylogenetics has revolutionized the study of not only evolution but also disparate fields such as genomics, bioinformatics, epidemiology, ecology, microbiology, molecular biology and biochemistry. Particularly significant are its achievements in population genetics as a result of the development of coalescent theory, which have contributed to more accurate model-based parameter estimation and explicit hypothesis testing. The study of the evolution of many microorganisms, and HIV in particular, have benefited from these new methodologies. HIV is well suited for such sophisticated population analyses because of its large population sizes, short generation times, high substitution rates and relatively small genomes. All these factors make HIV an ideal and fascinating model to study molecular evolution in real time. Here we review the significant advances made in HIV evolution through the application of phylogenetic approaches. We first examine the relative roles of mutation and recombination on the molecular evolution of HIV and its adaptive response to drug therapy and tissue allocation. We then review some of the fundamental questions in HIV evolution in relation to its origin and diversification and describe some of the insights gained using phylogenies. Finally, we show how phylogenetic analysis has advanced our knowledge of HIV dynamics (i.e., phylodynamics).

Evolution of viral genomes: interplay between selection, recombination, and other forces

RNA viruses evolve very rapidly, often recombine, and are subject to strong host (immune response) and anthropogenic (antiretroviral drugs) selective forces. Given their compact and extensively sequenced genomes, comparative analysis of RNA viral data can provide important insights into the molecular mechanisms of adaptation, pathogenicity, immune evasion, and drug resistance. In this chapter, we present an example-based overview of recent advances in evolutionary models and statistical approaches that enable screening viral alignments for evidence of adaptive change in the presence of recombination, detecting bursts of directional adaptive evolution associated with the phenotypic changes, and detecting of coevolving sites in viral genes.

Analysing recombination in nucleotide sequences

Throughout the living world, genetic recombination and nucleotide substitution are the primary processes that create the genetic variation upon which natural selection acts. Just as analyses of substitution patterns can reveal a great deal about evolution, so too can analyses of recombination. Evidence of genetic recombination within the genomes of apparently asexual species can equate with evidence of cryptic sexuality. In sexually reproducing species, nonrandom patterns of sequence exchange can provide direct evidence of population subdivisions that prevent certain individuals from mating. Although an interesting topic in its own right, an important reason for analysing recombination is to account for its potentially disruptive influences on various phylogenetic-based molecular evolution analyses. Specifically, the evolutionary histories of recombinant sequences cannot be accurately described by standard bifurcating phylogenetic trees. Taking recombination into account can therefore be pivotal to the success of selection, molecular clock and various other analyses that require adequate modelling of shared ancestry and draw increased power from accurately inferred phylogenetic trees. Here, we review various computational approaches to studying recombination and provide guidelines both on how to gain insights into this important evolutionary process and on how it can be properly accounted for during molecular evolution studies.

Reducing chimera formation during PCR amplification to ensure accurate genotyping

Measurements of population diversity are fundamental to the reconstruction of the evolutionary and epidemiological history of organisms. Commonly used protocols to measure population diversity using the polymerase chain reaction (PCR) are prone to the introduction of artificial chimeras. These are often difficult to detect and can confound the correct interpretation of results due to the false generation of recombinants when the underlying DNA sample contains multiple distinct templates. This study presents a standardised procedure to suppress the formation of artificial chimeras during PCR amplification. The solution is based on the accurate determination of the efficiency and end point of the log-linear phase of a PCR. This procedure will facilitate the generation of data sets that more accurately reflect the underlying population diversity rather than artifacts introduced by the process itself.