Uneven genetic robustness of HIV-1 integrase

Structural Impact of Ex Vivo Resistance Mutations on HIV-1 Integrase Polymers Induced by Allosteric Inhibitors

HIV-1 integrase (IN) is targeted by two classes of antivirals: integrase strand transfer inhibitors (INSTIs), which bind to the active site within the catalytic core domain (CCD), and allosteric integrase inhibitors (ALLINIs), which bind at the CCD dimer interface. ALLINIs were initially designed to disrupt interactions with the cellular cofactor LEDGF/p75, but it has become clear that ALLINIs primarily act by promoting formation of aberrant integrase polymers. The ALLINIs achieve this by stabilizing ectopic intermolecular interactions between the CCD dimer and the integrase carboxy-terminal domain (CTD), which disrupts viral maturation. Previously, we determined the structure of full-length HIV-1 IN bound to the ALLINI GSK1264 at 4.4 Å resolution, revealing its polymerization mechanism. More recently, we reported the X-ray crystal structure of a minimal ternary complex between CCD, CTD, and the ALLINI BI-224436 at a higher resolution. In this study, we improve the original 4.4 Å structure using this higher-resolution information and report two new structures of full-length HIV-1 IN harboring escape mutations in the CCD (Trp131Cys) or CTD (Asn222Lys) bound with the prototype ALLINI BI-D at 4.5 Å. These structures reveal perturbations to the tertiary organization associated with escape substitutions, which correlate with their reduced ability to form ectopic ALLINI-induced polymers in vitro. These findings suggest a general structural mechanism of ALLINI resistance and provide insights for the design of improved ALLINIs.

Altered Host microRNAomics in HIV Infections: Therapeutic Potentials and Limitations

microRNAs have emerged as essential regulators of health and disease, attracting significant attention from researchers across diverse disciplines. Following their identification as noncoding oligonucleotides intricately involved in post-transcriptional regulation of protein expression, extensive efforts were devoted to elucidating and validating their roles in fundamental metabolic pathways and multiple pathologies. Viral infections are significant modifiers of the host microRNAome. Specifically, the Human Immunodeficiency Virus (HIV), which affects approximately 39 million people worldwide and has no definitive cure, was reported to induce significant changes in host cell miRNA profiles. Identifying and understanding the effects of the aberrant microRNAome holds potential for early detection and therapeutic designs. This review presents a comprehensive overview of the impact of HIV on host microRNAome. We aim to review the cause-and-effect relationship between the HIV-induced aberrant microRNAome that underscores miRNA's therapeutic potential and acknowledge its limitations.

Multimodal Functionalities of HIV-1 Integrase

Integrase is the retroviral protein responsible for integrating reverse transcripts into cellular genomes. Co-packaged with viral RNA and reverse transcriptase into capsid-encased viral cores, human immunodeficiency virus 1 (HIV-1) integrase has long been implicated in reverse transcription and virion maturation. However, the underlying mechanisms of integrase in these non-catalytic-related viral replication steps have remained elusive. Recent results have shown that integrase binds genomic RNA in virions, and that mutational or pharmacological disruption of integrase-RNA binding yields eccentric virion particles with ribonucleoprotein complexes situated outside of the capsid shell. Such viruses are defective for reverse transcription due to preferential loss of integrase and viral RNA from infected target cells. Parallel research has revealed defective integrase-RNA binding and eccentric particle formation as common features of class II integrase mutant viruses, a phenotypic grouping of viruses that display defects at steps beyond integration. In light of these new findings, we propose three new subclasses of class II mutant viruses (a, b, and c), all of which are defective for integrase-RNA binding and particle morphogenesis, but differ based on distinct underlying mechanisms exhibited by the associated integrase mutant proteins. We also assess how these findings inform the role of integrase in HIV-1 particle maturation.

Vulnerable targets in HIV-1 Pol for attenuation-based vaccine design

Identification of viral immune escape mutations that compromise HIV's ability to replicate may aid rational attenuation-based vaccine design. Previously we reported amino acids associated with altered viral replication capacity (RC) from a sequence-function analysis of 487 patient-derived RT-integrase sequences. In this study, site-directed mutagenesis experiments were performed to validate the effect of these mutations on RC. Viral reverse transcripts were measured by quantitative PCR and structural modelling was performed to gain further insight into the effect of reverse transcriptase (RT) mutations on reverse transcription. RT-integrase variants in or flanking cytotoxic T cell epitopes in the RT palm (158S), RT thumb (241I and 257V) and integrase catalytic core domain (124N) were confirmed to significantly reduce RC. RT mutants showed a delayed initiation of viral DNA synthesis. Structural models provide insight into how these attenuating RT mutations may affect amino acid interactions in the helix clamp, primer grip and catalytic site regions.

CRISPR Genome Editing Applied to the Pathogenic Retrovirus HTLV-1

CRISPR editing of retroviral proviruses has been limited to HIV-1. We propose human T-cell leukemia virus type 1 (HTLV-1) as an excellent model to advance CRISPR/Cas9 genome editing technologies against actively expressing and latent retroviral proviruses. HTLV-1 is a tumorigenic human retrovirus responsible for the development of both leukemia/lymphoma (ATL) and a neurological disease (HAM/TSP). The virus immortalizes and persists in CD4⁺ T lymphocytes that survive for the lifetime of the host. The most important drivers of HTLV-1-mediated transformation and proliferation are the tax and hbz viral genes. Tax, transcribed from the plus-sense or genome strand, is essential for de novo infection and cellular immortalization. Hbz, transcribed from the minus-strand, supports proliferation and survival of infected cells in both its protein and mRNA forms. Abrogating the function or expression of tax and/or hbz by genome editing and mutagenic double-strand break repair may disable HTLV-1-infected cell growth/survival and prevent immune modulatory effects and ultimately HTLV-1-associated disease. In addition, the HTLV-1 viral genome is highly conserved with remarkable sequence homogeneity, both within the same host and even among different HTLV isolates. This offers more focused guide RNA targeting. In addition, there are several well-established animal models for studying HTLV-1 infection in vivo as well as cell immortalization in vitro. Therefore, studies with HTLV-1 may provide a better basis to assess and advance in vivo genome editing against retroviral infections.

A de novo approach to inferring within-host fitness effects during untreated HIV-1 infection

In the absence of effective antiviral therapy, HIV-1 evolves in response to the within-host environment, of which the immune system is an important aspect. During the earliest stages of infection, this process of evolution is very rapid, driven by a small number of CTL escape mutations. As the infection progresses, immune escape variants evolve under reduced magnitudes of selection, while competition between an increasing number of polymorphic alleles (i.e., clonal interference) makes it difficult to quantify the magnitude of selection acting upon specific variant alleles. To tackle this complex problem, we developed a novel multi-locus inference method to evaluate the role of selection during the chronic stage of within-host infection. We applied this method to targeted sequence data from the p24 and gp41 regions of HIV-1 collected from 34 patients with long-term untreated HIV-1 infection. We identify a broad distribution of beneficial fitness effects during infection, with a small number of variants evolving under strong selection and very many variants evolving under weaker selection. The uniquely large number of infections analysed granted a previously unparalleled statistical power to identify loci at which selection could be inferred to act with statistical confidence. Our model makes no prior assumptions about the nature of alleles under selection, such that any synonymous or non-synonymous variant may be inferred to evolve under selection. However, the majority of variants inferred with confidence to be under selection were non-synonymous in nature, and in most cases were have previously been associated with either CTL escape in p24 or neutralising antibody escape in gp41. We also identified a putative new CTL escape site (residue 286 in gag), and a region of gp41 (including residues 644, 648, 655 in env) likely to be associated with immune escape. Sites inferred to be under selection in multiple hosts have high within-host and between-host diversity although not all sites with high between-host diversity were inferred to be under selection at the within-host level. Our identification of selection at sites associated with resistance to broadly neutralising antibodies (bNAbs) highlights the need to fully understand the role of selection in untreated individuals when designing bNAb based therapies.

A CRISPR/Cas9 library to map the HIV-1 provirus genetic fitness

The integrated proviral genome is the major barrier to a cure for HIV-1 infection. Genome editing technologies, such as CRISPR/Cas9, may disable or remove the HIV-1 provirus by introducing DNA double strand breaks at sequence specific sites in the viral genome. Host DNA repair by the error-prone non-homologous end joining pathway generates mutagenic insertions or deletions at the break. CRISPR/Cas9 editing has been shown to reduce replication competent viral genomes in cell culture, but only a minority of possible genome editing targets have been assayed. Currently there is no map of double strand break genetic fitness for HIV-1 to inform the choice of editing targets. However, CRISPR/Cas9 genome editing makes it possible to target double strand breaks along the length of the provirus to generate a double strand break genetic fitness map. We identified all possible HIV-1 targets with different bacterial species of CRISPR/Cas9. This library of guide RNAs was evaluated for GC content and potential off-target sites in the human genome. Complexity of the library was reduced by eliminating duplicate guide RNA targets in the HIV-1 long terminal repeats and targets in the env gene. Although the HIV-1 genome is AT-rich, the S. pyogenes CRISPR/Cas9 with the proto-spacer adjacent motif NGG offers the most HIV-1 guide RNAs. This library of HIV-1 guide RNAs may be used to generate a double strand break genetic fragility map to be further applied to any genome editing technology designed for the HIV-1 provirus. Keywords: HIV-1; genome editing; CRISPR; genetic fitness; guide RNAs.

Effects of Arbovirus Multi-Host Life Cycles on Dinucleotide and Codon Usage Patterns

Arthropod-borne viruses (arboviruses) of vertebrates including dengue, zika, chikungunya, Rift Valley fever, and blue tongue viruses cause extensive morbidity and mortality in humans, agricultural animals, and wildlife across the globe. As obligate intercellular pathogens, arboviruses must be well adapted to the cellular and molecular environment of both their arthropod (invertebrate) and vertebrate hosts, which are vastly different due to hundreds of millions of years of separate evolution. Here we discuss the comparative pressures on arbovirus RNA genomes as a result of a dual host life cycle, focusing on pressures that do not alter amino acids. We summarize what is currently known about arboviral genetic composition, such as dinucleotide and codon usage, and how cyclical infection of vertebrate and invertebrate hosts results in different genetic profiles compared with single-host viruses. To serve as a comparison, we compile what is known about arthropod tRNA, dinucleotide, and codon usages and compare this with vertebrates. Additionally, we discuss the potential roles of genetic robustness in arboviral evolution and how it may vary from other viruses. Overall, both arthropod and vertebrate hosts influence the resulting genetic composition of arboviruses, but a great deal remains to be investigated.

Evidence for the Selective Basis of Transition-to-Transversion Substitution Bias in Two RNA Viruses

The substitution rates of transitions are higher than expected by chance relative to those of transversions. Many have argued that selection disfavors transversions, as nonsynonymous transversions are less likely to conserve biochemical properties of the original amino acid. Only recently has it become feasible to directly test this selective hypothesis by comparing the fitness effects of a large number of transition and transversion mutations. For example, a recent study of six viruses and one beta-lactamase gene did not find evidence supporting the selective hypothesis. Here, we analyze the relative fitness effects of transition and transversion mutations from our recently published genome-wide study of mutational fitness effects in influenza virus. In contrast to prior work, we find that transversions are significantly more detrimental than transitions. Using what we believe to be an improved statistical framework, we also identify a similar trend in two HIV data sets. We further demonstrate a fitness difference in transition and transversion mutations using four deep mutational scanning data sets of influenza virus and HIV, which provided adequate statistical power. We find that three of the most commonly cited radical/conservative amino acid categories are predictive of fitness, supporting their utility in studies of positive selection and codon usage bias. We conclude that selection is a major contributor to the transition:transversion substitution bias in viruses and that this effect is only partially explained by the greater likelihood of transversion mutations to cause radical as opposed to conservative amino acid changes.

Cryo-EM structures and atomic model of the HIV-1 strand transfer complex intasome

Like all retroviruses, HIV-1 irreversibly inserts a viral DNA (vDNA) copy of its RNA genome into host target DNA (tDNA). The intasome, a higher-order nucleoprotein complex composed of viral integrase (IN) and the ends of linear vDNA, mediates integration. Productive integration into host chromatin results in the formation of the strand transfer complex (STC) containing catalytically joined vDNA and tDNA. HIV-1 intasomes have been refractory to high-resolution structural studies. We used a soluble IN fusion protein to facilitate structural studies, through which we present a high-resolution cryo-electron microscopy (cryo-EM) structure of the core tetrameric HIV-1 STC and a higher-order form that adopts carboxyl-terminal domain rearrangements. The distinct STC structures highlight how HIV-1 can use the common retroviral intasome core architecture to accommodate different IN domain modules for assembly.

The Envelope Gene of Transmitted HIV-1 Resists a Late Interferon Gamma-Induced Block

Type I interferon (IFN) signaling engenders an antiviral state that likely plays an important role in constraining HIV-1 transmission and contributes to defining subsequent AIDS pathogenesis. Type II IFN (IFN-γ) also induces an antiviral state but is often primarily considered to be an immunomodulatory cytokine. We report that IFN-γ stimulation can induce an antiviral state that can be both distinct from that of type I interferon and can potently inhibit HIV-1 in primary CD4⁺ T cells and a number of human cell lines. Strikingly, we find that transmitted/founder (TF) HIV-1 viruses can resist a late block that is induced by type II IFN, and the use of chimeric IFN-γ-sensitive/resistant viruses indicates that interferon resistance maps to the env gene. Simultaneously, in vitro evolution also revealed that just a single amino acid substitution in the envelope can confer substantial resistance to IFN-mediated inhibition. Thus, the env gene of transmitted HIV-1 confers resistance to a late block that is phenotypically distinct from blocks previously described to be resisted by env and is therefore mediated by unknown IFN-γ-stimulated factor(s) in human CD4⁺ T cells and cell lines. This important unidentified block could play a key role in constraining HIV-1 transmission.

IMPORTANCE The human immune system can hinder invading pathogens through interferon (IFN) signaling. One consequence of this signaling is that cells enter an antiviral state, increasing the levels of hundreds of defenses that can inhibit the replication and spread of viruses. The majority of HIV-1 infections result from a single virus particle (the transmitted/founder) that makes it past these defenses and colonizes the host. Thus, the founder virus is hypothesized to be a relatively interferon-resistant entity. Here, we show that certain HIV-1 envelope genes have the unanticipated ability to resist specific human defenses mediated by different types of interferons. Strikingly, the envelope gene from a founder HIV-1 virus is far better at evading these defenses than the corresponding gene from a common HIV-1 lab strain. Thus, these defenses could play a role in constraining the transmission of HIV-1 and may select for transmitted viruses that are resistant to this IFN-mediated inhibition.

In vivo mutation rates and the landscape of fitness costs of HIV-1

Mutation rates and fitness costs of deleterious mutations are difficult to measure in vivo but essential for a quantitative understanding of evolution. Using whole genome deep sequencing data from longitudinal samples during untreated HIV-1 infection, we estimated mutation rates and fitness costs in HIV-1 from the dynamics of genetic variation. At approximately neutral sites, mutations accumulate with a rate of 1.2 × 10-5 per site per day, in agreement with the rate measured in cell cultures. We estimated the rate from G to A to be the largest, followed by the other transitions C to T, T to C, and A to G, while transversions are less frequent. At other sites, mutations tend to reduce virus replication. We estimated the fitness cost of mutations at every site in the HIV-1 genome using a model of mutation selection balance. About half of all non-synonymous mutations have large fitness costs (>10 percent), while most synonymous mutations have costs <1 percent. The cost of synonymous mutations is especially low in most of pol where we could not detect measurable costs for the majority of synonymous mutations. In contrast, we find high costs for synonymous mutations in important RNA structures and regulatory regions. The intra-patient fitness cost estimates are consistent across multiple patients, indicating that the deleterious part of the fitness landscape is universal and explains a large fraction of global HIV-1 group M diversity.

HIV-1 Integrase Binds the Viral RNA Genome and Is Essential during Virion Morphogenesis

While an essential role of HIV-1 integrase (IN) for integration of viral cDNA into human chromosome is established, studies with IN mutants and allosteric IN inhibitors (ALLINIs) have suggested that IN can also influence viral particle maturation. However, it has remained enigmatic as to how IN contributes to virion morphogenesis. Here, we demonstrate that IN directly binds the viral RNA genome in virions. These interactions have specificity, as IN exhibits distinct preference for select viral RNA structural elements. We show that IN substitutions that selectively impair its binding to viral RNA result in eccentric, non-infectious virions without affecting nucleocapsid-RNA interactions. Likewise, ALLINIs impair IN binding to viral RNA in virions of wild-type, but not escape mutant, virus. These results reveal an unexpected biological role of IN binding to the viral RNA genome during virion morphogenesis and elucidate the mode of action of ALLINIs.

A Highly Conserved Residue in HIV-1 Nef Alpha Helix 2 Modulates Protein Expression

The HIV-1 Nef protein has been established as a key pathogenic determinant of HIV/AIDS, but there is little knowledge of how the extensive genetic diversity of HIV-1 affects Nef function. Upon compiling a set of subtype-specific reference strains, we identified a subtype C reference strain, C.BR92025, that contained natural polymorphisms at otherwise highly conserved residues 13, 84, and 92. Interestingly, strain C.BR92025 Nef displayed impaired Nef function and had decreased protein expression. We have demonstrated that strain C.BR92025 Nef has a higher rate of protein turnover than highly expressed Nef proteins and that this higher rate of protein turnover is due to an alanine-to-valine substitution at Nef residue 84. These findings highlight residue A84 as a major determinant of HIV-1 Nef expression.

Extensive genetic diversity is a defining characteristic of human immunodeficiency virus type 1 (HIV-1) and poses a significant barrier to the development of an effective vaccine. To better understand the impact of this genetic diversity on the HIV-1 pathogenic factor Nef, we compiled a panel of reference strains from the NIH Los Alamos HIV Database. Initial sequence analysis identified point mutations at Nef residues 13, 84, and 92 in subtype C reference strain C.BR92025 from Brazil. Functional analysis revealed impaired major histocompatibility complex class I and CD4 downregulation of strain C.BR92025 Nef, which corresponded to decreased protein expression. Metabolic labeling demonstrated that strain C.BR92025 Nef has a greater rate of protein turnover than subtype B reference strain B.JRFL that, on the basis of mutational analysis, is related to Nef residue A84. An alanine-to-valine substitution at position 84, located in alpha helix 2 of Nef, was sufficient to alter the rate of turnover of an otherwise highly expressed Nef protein. In conclusion, these findings highlight HIV-1 Nef residue A84 as a major determinant of protein expression that may offer an additional avenue to disrupt or mediate the effects of this key HIV-1 pathogenic factor.

IMPORTANCE The HIV-1 Nef protein has been established as a key pathogenic determinant of HIV/AIDS, but there is little knowledge of how the extensive genetic diversity of HIV-1 affects Nef function. Upon compiling a set of subtype-specific reference strains, we identified a subtype C reference strain, C.BR92025, that contained natural polymorphisms at otherwise highly conserved residues 13, 84, and 92. Interestingly, strain C.BR92025 Nef displayed impaired Nef function and had decreased protein expression. We have demonstrated that strain C.BR92025 Nef has a higher rate of protein turnover than highly expressed Nef proteins and that this higher rate of protein turnover is due to an alanine-to-valine substitution at Nef residue 84. These findings highlight residue A84 as a major determinant of HIV-1 Nef expression.

HIV-1 Protease, Reverse Transcriptase, and Integrase Variation

HIV-1 protease (PR), reverse transcriptase (RT), and integrase (IN) variability presents a challenge to laboratories performing genotypic resistance testing. This challenge will grow with increased sequencing of samples enriched for proviral DNA such as dried blood spots and increased use of next-generation sequencing (NGS) to detect low-abundance HIV-1 variants. We analyzed PR and RT sequences from >100,000 individuals and IN sequences from >10,000 individuals to characterize variation at each amino acid position, identify mutations indicating APOBEC-mediated G-to-A editing, and identify mutations resulting from selective drug pressure. Forty-seven percent of PR, 37% of RT, and 34% of IN positions had one or more amino acid variants with a prevalence of ≥1%. Seventy percent of PR, 60% of RT, and 60% of IN positions had one or more variants with a prevalence of ≥0.1%. Overall 201 PR, 636 RT, and 346 IN variants had a prevalence of ≥0.1%. The median intersubtype prevalence ratios were 2.9-, 2.1-, and 1.9-fold for these PR, RT, and IN variants, respectively. Only 5.0% of PR, 3.7% of RT, and 2.0% of IN variants had a median intersubtype prevalence ratio of ≥10-fold. Variants at lower prevalences were more likely to differ biochemically and to be part of an electrophoretic mixture compared to high-prevalence variants. There were 209 mutations indicative of APOBEC-mediated G-to-A editing and 326 mutations nonpolymorphic treatment selected. Identification of viruses with a high number of APOBEC-associated mutations will facilitate the quality control of dried blood spot sequencing. Identifying sequences with a high proportion of rare mutations will facilitate the quality control of NGS.

IMPORTANCE Most antiretroviral drugs target three HIV-1 proteins: PR, RT, and IN. These proteins are highly variable: many different amino acids can be present at the same position in viruses from different individuals. Some of the amino acid variants cause drug resistance and occur mainly in individuals receiving antiretroviral drugs. Some variants result from a human cellular defense mechanism called APOBEC-mediated hypermutation. Many variants result from naturally occurring mutation. Some variants may represent technical artifacts. We studied PR and RT sequences from >100,000 individuals and IN sequences from >10,000 individuals to quantify variation at each amino acid position in these three HIV-1 proteins. We performed analyses to determine which amino acid variants resulted from antiretroviral drug selection pressure, APOBEC-mediated editing, and naturally occurring variation. Our results provide information essential to clinical, research, and public health laboratories performing genotypic resistance testing by sequencing HIV-1 PR, RT, and IN.

On the Causes of Evolutionary Transition:Transversion Bias

A pattern in which nucleotide transitions are favored several fold over transversions is common in molecular evolution. When this pattern occurs among amino acid replacements, explanations often invoke an effect of selection, on the grounds that transitions are more conservative in their effects on proteins. However, the underlying hypothesis of conservative transitions has never been tested directly. Here we assess support for this hypothesis using direct evidence: the fitness effects of mutations in actual proteins measured via individual or paired growth experiments. We assembled data from 8 published studies, ranging in size from 24 to 757 single-nucleotide mutations that change an amino acid. Every study has the statistical power to reveal significant effects of amino acid exchangeability, and most studies have the power to discern a binary conservative-vs-radical distinction. However, only one study suggests that transitions are significantly more conservative than transversions. In the combined set of 1,239 replacements (544 transitions, 695 transversions), the chance that a transition is more conservative than a transversion is 53 % (95 % confidence interval 50 to 56) compared with the null expectation of 50 %. We show that this effect is not large compared with that of most biochemical factors, and is not large enough to explain the several-fold bias observed in evolution. In short, the available data have the power to verify the “conservative transitions” hypothesis if true, but suggest instead that selection on proteins plays at best a minor role in the observed bias.

Extremely High Mutation Rate of HIV-1 In Vivo

Rates of spontaneous mutation critically determine the genetic diversity and evolution of RNA viruses. Although these rates have been characterized in vitro and in cell culture models, they have seldom been determined in vivo for human viruses. Here, we use the intrapatient frequency of premature stop codons to quantify the HIV-1 genome-wide rate of spontaneous mutation in DNA sequences from peripheral blood mononuclear cells. This reveals an extremely high mutation rate of (4.1 ± 1.7) × 10⁻³ per base per cell, the highest reported for any biological entity. Sequencing of plasma-derived sequences yielded a mutation frequency 44 times lower, indicating that a large fraction of viral genomes are lethally mutated and fail to reach plasma. We show that the HIV-1 reverse transcriptase contributes only 2% of mutations, whereas 98% result from editing by host cytidine deaminases of the A3 family. Hypermutated viral sequences are less abundant in patients showing rapid disease progression compared to normal progressors, highlighting the antiviral role of A3 proteins. However, the amount of A3-mediated editing varies broadly, and we find that low-edited sequences are more abundant among rapid progressors, suggesting that suboptimal A3 activity might enhance HIV-1 genetic diversity and pathogenesis.

The rate of spontaneous mutation of the HIV-1 genome within its human host is exceptionally high, is mostly driven by host cytidine deaminases, and probably plays a role in disease progression.

The high levels of genetic diversity of the HIV-1 virus grant it the ability to escape the immune system, to rapidly evolve drug resistance, and to circumvent vaccination strategies. However, our knowledge of HIV-1 mutation rates has been largely restricted to in vitro and cell culture studies because of the inherent complexity of measuring these rates in vivo. Here, by analyzing the frequency of premature stop codons, we show that the HIV-1 mutation rate in vivo is two orders of magnitude higher than that predicted by in vitro studies, making it the highest reported mutation rate for any biological system. A large component of this rate is from host cellular cytidine deaminases, which induce mutations in the viral DNA as a defense mechanism. While the HIV-1 genome is hypermutated in blood cells, only a very small fraction of these mutations reach the plasma, indicating that many viruses are defective as a result of the extremely high mutation load. In addition, we find that the HIV-1 mutation rate tends to be higher in patients showing normal disease progression than in those undergoing rapid progression, emphasizing the negative impact on viral fitness of hypermutation by host cytidine deaminases. However, we also observe subpopulations of weakly-mutated viral genomes whose sequence diversity may influence viral pathogenesis. Our work highlights the fine balance for HIV-1 between enough mutation to evade host responses and too much mutation that can inactivate the virus.

Epistasis as a determinant of the HIV-1 protease's robustness to mutation

The robustness of phenotypes to mutation is critical to protein evolution; robustness may be an adaptive trait if it promotes evolution. We hypothesised that native proteins subjected to natural selection in vivo should be more robust than proteins generated in vitro in the absence of natural selection. We compared the mutational robustness of two human immunodeficiency virus type 1 (HIV-1) proteases with comparable catalytic efficiencies, one isolated from an infected individual and the second generated in vitro via random mutagenesis. Single mutations in the protease (82 and 60 in the wild-type and mutant backgrounds, respectively) were randomly generated in vitro and the catalytic efficiency of each mutant was determined. No differences were observed between these two protease variants when lethal, neutral, and deleterious mutations were compared (P = 0.8025, chi-squared test). Similarly, average catalytic efficiency (−72.6% and −64.5%, respectively) did not significantly differ between protease mutant libraries (P = 0.3414, Mann Whitney test). Overall, the two parental proteins displayed similar mutational robustness. Importantly, strong and widespread epistatic interactions were observed when the effect of the same mutation was compared in both proteases, suggesting that epistasis can be a key determinant of the robustness displayed by the in vitro generated protease.